Changes in zebrafish gonadal index (GSI) after feeding Aurantiochytrium sp. for 56 days. Each bar represents the mean ± SEM (n = 3). Statistically significant differences are expressed as * (* p < 0.05, ** p < 0.01).

Histological analysis of zebrafish fed 1% Aurantiochytrium sp. Extract for 56 days. (AD) is the 1% TE group, (EH) is the control group, where (D) and (H) are oil red stains, and the rest are H&E stains. Blue letters indicate oocyte development stages, including perinucleolar oocytes (PO), corticoalveolar oocytes (CAO), early vitellogenic oocytes (EVO), late and mature oocytes (LMO), and post-ovulatory follicles (POF). Scale: 100 μm.

Germinal vesicle breakdown (GVBD) in zebrafish oocytes exposed to different concentrations of Aurantiochytrium sp. active product (TE) and 17α, 20β-dihydroxy-4-progesterone-3-one (DHP). Each bar represents the mean ± SEM (n = 3). The statistically significant difference is expressed as * (* p < 0.05, ** p < 0.01).

The effect of 0.5% TE on the expression of genes related to oocyte maturation in zebrafish in vivo (AF). The data are expressed as mean ± standard error (SEM) (n = 3). The statistically significant difference is expressed as * (* p < 0.05, ** p < 0.01).

Differential expression of gonadal genes in the control group and 1% TE group. The green dots represent down-regulated differentially expressed genes, the blue dots represent up-regulated differentially expressed genes, and the black dots represent non-differentially expressed genes. The significance criteria are q < 0.05 and Fold Change ≥ 2.

Transcriptome differential gene validation and its cluster analysis (A) RT-qPCR differential expression verification. The error line indicates SEM (n = 3). (B) The consistency of log2 folds changes between RNA-Seq data (x-axis) and RT-qPCR analysis (y-axis). (C) Differentially expressed gene expression pattern cluster diagram, used to illustrate the overall pattern of gene expression in different gonadal samples.

Functional classification of gonadal differentially expressed genes GO Enrichment bubble map (A); Cell component category (B); Molecular function category; and (C) Biological process category.

Functional classification of gonadal differentially expressed genes GO Enrichment bubble map (A); Cell component category (B); Molecular function category; and (C) Biological process category.

(A) Scatter plot of KEGG pathway enrichment of differentially expressed genes, the vertical coordinate indicates the pathway name, the horizontal coordinate is the Enrichment Factor, the color represents the q-value; q-value is the p-value after correction for multiple hypothesis testing, and the size of the circle indicates the number of enriched genes in the pathway. (B) Differentially expressed genes KEGG classification diagram, the vertical axis is the name of the KEGG metabolic pathway, the left is the pathway name, and the right is the corresponding classification category of each pathway. The same column color indicates the same category.

Acknowledgments
This image is the copyrighted work of the attributed author or publisher, and ZFIN has permission only to display this image to its users. Additional permissions should be obtained from the applicable author or publisher of the image. Full text @ Animals (Basel)