Differentiation In Zebrafish Research Articles

Two phases of gonadal sex differentiation in zebrafish with ZZ/ZW sex determination system

Zebrafish sex chromosomes have been identified in the wild Nadia (NA) strain, and its sex determination belongs to the female-heterogametic ZZ/ZW system. Here, we investigate the correlation between ZZ/ZW sex chromosomes in the NA strain with sex-related factors, and sort out the complicated process of sex determination in zebrafish. Two phases exist during zebrafish sex differentiation. In the first phase, ZW gonads differentiate into juvenile ovary while ZZ gonads remain indifferent. In the second phase, ZW gonads either continue ovary development or undergo female-to-male transition, while ZZ gonads undergo direct male development. The W chromosome may contribute to the first phase while the abundance of germ cells and other factors may be involved in the second phase of sex differentiation in zebrafish.

Adolescent exposure to tris(1,3-dichloro-2-propyl) phosphate (TCPP) induces reproductive toxicity in zebrafish through hypothalamic-pituitary-gonadal axis disruption

Tris(1,3-dichloro-2-propyl) phosphate (TCPP), a prevalent organophosphorus flame retardant in aquatic environments, has raised significant concerns regarding its ecological risks. This study aims to explore the impacts of TCPP on the reproductive functions of zebrafish and delineate its gender-related toxic mechanisms. By assessing the effects on zebrafish of 10 mg/L TCPP exposure from 30 to 120 days post-fertilization (dpf), we thoroughly evaluated the reproductive capability and endocrine system alterations. Our findings indicated that TCPP exposure disrupted gender differentiation in zebrafish and markedly impaired their reproductive capacity, resulting in decreased egg laying and offspring development quality. Histological analyses of gonadal tissues showed an abnormal increase in immature oocytes in females and a reduction in mature sperm count and spermatogonial structure integrity in males, collectively leading to compromised embryo quality. Additionally, molecular docking results indicated that TCPP showed a strong affinity for estrogen receptors, and TCPP-treated zebrafish exhibited imbalanced sex hormones and increased estrogen receptor expression. Alterations in genes associated with the hypothalamic-pituitary-gonadal (HPG) axis and activation of the steroidogenesis pathway suggested that TCPP targets the HPG axis to regulate sex hormone homeostasis. Tamoxifen (TAM), as a competitive inhibitor of estrogen, exhibited a biphasic effect, as evidenced by the counteraction of TCPP-induced effects in both male and female zebrafish after TAM addition. Overall, our study underscored the gender-dependent reproductive toxicity of TCPP exposure in zebrafish, characterized by diminished reproductive capacity and hormonal disturbances, likely due to interference in the HPG axis and steroidogenesis pathways. These findings emphasize the critical need to consider gender differences in chemical risk assessments for ecosystems and highlight the importance of understanding the mechanisms underlying the effects of chemical pollutants on the reproductive health of aquatic species.

New insights into the all-testis differentiation in zebrafish with compromised endogenous androgen and estrogen synthesis.

The regulatory mechanism of gonadal sex differentiation, which is complex and regulated by multiple factors, remains poorly understood in teleosts. Recently, we have shown that compromised androgen and estrogen synthesis with increased progestin leads to all-male differentiation with proper testis development and spermatogenesis in cytochrome P450 17a1 (cyp17a1)-/- zebrafish. In the present study, the phenotypes of female-biased sex ratio were positively correlated with higher Fanconi anemia complementation group L (fancl) expression in the gonads of doublesex and mab-3 related transcription factor 1 (dmrt1)-/- and cyp17a1-/-;dmrt1-/- fish. The additional depletion of fancl in cyp17a1-/-;dmrt1-/- zebrafish reversed the gonadal sex differentiation from all-ovary to all-testis (in cyp17a1-/-;dmrt1-/-;fancl-/- fish). Luciferase assay revealed a synergistic inhibitory effect of Dmrt1 and androgen signaling on fancl transcription. Furthermore, an interaction between Fancl and the apoptotic factor Tumour protein p53 (Tp53) was found in vitro. The interaction between Fancl and Tp53 was observed via the WD repeat domain (WDR) and C-terminal domain (CTD) of Fancl and the DNA binding domain (DBD) of Tp53, leading to the K48-linked polyubiquitination degradation of Tp53 activated by the ubiquitin ligase, Fancl. Our results show that testis fate in cyp17a1-/- fish is determined by Dmrt1, which is thought to stabilize Tp53 by inhibiting fancl transcription during the critical stage of sexual fate determination in zebrafish.

Anp32e protects against accumulation of H2A.Z at Sox motif containing promoters during zebrafish gastrulation

Epigenetic regulation of chromatin states is crucial for proper gene expression programs and progression during development, but precise mechanisms by which epigenetic factors influence differentiation remain poorly understood. Here we find that the histone variant H2A.Z accumulates at Sox motif-containing promoters during zebrafish gastrulation while neighboring genes become transcriptionally active. These changes coincide with reduced expression of anp32e, the H2A.Z histone removal chaperone, suggesting that loss of Anp32e may lead to increases in H2A.Z binding during differentiation. Remarkably, genetic removal of Anp32e in embryos leads to H2A.Z accumulation prior to gastrulation and developmental genes become precociously active. Accordingly, H2A.Z accumulation occurs most extensively at Sox motif-associated genes, including many which are normally activated following gastrulation. Altogether, our results provide compelling evidence for a mechanism in which Anp32e preferentially restricts H2A.Z accumulation at Sox motifs to regulate the initial phases of developmental differentiation in zebrafish.

Genetic mechanisms of multiciliated cell development: from fate choice to differentiation in zebrafish and other models.

Multiciliated cells (MCCS) form bundles of cilia and their activities are essential for the proper development and physiology of many organ systems. Not surprisingly, defects in MCCs have profound consequences and are associated with numerous disease states. Here, we discuss the current understanding of MCC formation, with a special focus on the genetic and molecular mechanisms of MCC fate choice and differentiation. Furthermore, we cast a spotlight on the use of zebrafish to study MCC ontogeny and several recent advances made in understanding MCCs using this vertebrate model to delineate mechanisms of MCC emergence in the developing kidney.

Hhex and Prox1a synergistically dictate the hepatoblast to hepatocyte differentiation in zebrafish

The specification of endoderm cells to prospective hepatoblasts is the starting point for hepatogenesis. However, how a prospective hepatoblast gains the hepatic fate remains elusive. Previous studies have shown that loss-of-function of either hhex or prox1a alone causes a small liver phenotype but without abolishing the hepatocyte differentiation, suggesting that absence of either Hhex or Prox1a alone is not sufficient to block the hepatoblast differentiation. Here, via genetic studies of the zebrafish two single (hhex−/− and prox1a−/−) and one double (hhex−/−prox1a−/−) mutants, we show that simultaneous loss-of-function of the hhex and prox1a two genes does not block the endoderm cells to gain the hepatoblast potency but abolishes the hepatic differentiation from the prospective hepatoblast. Consequently, the hhex−/−prox1a−/− double mutant displays a liverless phenotype that cannot be rescued by the injection of bmp2a mRNA. Taken together, we provide strong evidences showing that Hhex teams with Prox1a to act as a master control of the differentiation of the prospective hepatoblasts towards hepatocytes.

The osteogenic and mineralogenic potential of the microalgae Skeletonema costatum and Tetraselmis striata CTP4 in fish models

Skeletal disorders are problematic aspects for the aquaculture industry as skeletal deformities, which affect most species of farmed fish, increase production costs and affect fish welfare. Following recent findings that show the presence of osteoactive compounds in marine organisms, we evaluated the osteogenic and mineralogenic potential of commercially available microalgae strains Skeletonema costatum and Tetraselmis striata CTP4 in several fish systems. Ethanolic extracts increased extracellular matrix mineralization in gilthead seabream (Sparus aurata) bone-derived cell cultures and promoted osteoblastic differentiation in zebrafish (Danio rerio) larvae. Long-term dietary exposure to both extracts increased bone mineralization in zebrafish and upregulated the expression of genes involved in bone formation (sp7, col1a1a, oc1, and oc2), bone remodeling (acp5a), and antioxidant defenses (cat, sod1). Extracts also improved the skeletal status of zebrafish juveniles by reducing the incidence of skeletal anomalies. Our results indicate that both strains of microalgae contain osteogenic and mineralogenic compounds, and that ethanolic extracts have the potential for an application in the aquaculture sector as dietary supplements to support fish bone health. Future studies should also identify osteoactive compounds and establish whether they can be used in human health to broaden the therapeutic options for bone erosive disorders such as osteoporosis.

Norethindrone suppress the germ cell development via androgen receptor resulting in male bias

Progestins are widely used and detected in surface waters, and can affect gonad development and sexual differentiation in fish. However, the toxicological mechanisms of sexual differentiation induced by progestins are not well understood. Here, we investigated the effects of norethindrone (NET) and androgen receptor (AR) antagonist flutamide (FLU) on gonadal differentiation in zebrafish from 21 dpf (days post-fertilization) to 49 dpf. The results showed that NET caused male bias, while FLU resulted in female bias at 49 dpf. The NET and FLU mixtures significantly decreased the percentage of males compared to the NET single exposure. Molecular docking analysis showed that FLU and NET had similar docking pocket and docking posture with AR resulting in competitively forming the hydrogen bond with Thr334 of AR. These results suggested that binding to AR was the molecular initiating event of sex differentiation induced by NET. Moreover, NET strongly decreased transcription of biomarker genes (dnd1, ddx4, dazl, piwil1 and nanos1) involved in germ cell development, while FLU significantly increased transcription of these target genes. There was an increase in the number of juvenile oocytes, which was consistent with the female bias in the combined groups. The bliss independence model analysis further showed that NET and FLU had antagonistic effect on transcription and histology during gonadal differentiation. Thus, NET suppressed the germ cell development via AR, resulting in male bias. Understanding the molecular initiation of sex differentiation in progestins is essential to provide a comprehensive biological basis for ecological risk assessment.

Leukocyte Tyrosine Kinase (Ltk) Is the Mendelian Determinant of the Axolotl Melanoid Color Variant.

The great diversity of color patterns observed among amphibians is largely explained by the differentiation of relatively few pigment cell types during development. Mexican axolotls present a variety of color phenotypes that span the continuum from leucistic to highly melanistic. The melanoid axolotl is a Mendelian variant characterized by large numbers of melanophores, proportionally fewer xanthophores, and no iridophores. Early studies of melanoid were influential in developing the single-origin hypothesis of pigment cell development, wherein it has been proposed that all three pigment cell types derive from a common progenitor cell, with pigment metabolites playing potential roles in directing the development of organelles that define different pigment cell types. Specifically, these studies identified xanthine dehydrogenase (XDH) activity as a mechanism for the permissive differentiation of melanophores at the expense of xanthophores and iridophores. We used bulked segregant RNA-Seq to screen the axolotl genome for melanoid candidate genes and identify the associated locus. Dissimilar frequencies of single-nucleotide polymorphisms were identified between pooled RNA samples of wild-type and melanoid siblings for a region on chromosome 14q. This region contains gephyrin (Gphn), an enzyme that catalyzes the synthesis of the molybdenum cofactor that is required for XDH activity, and leukocyte tyrosine kinase (Ltk), a cell surface signaling receptor that is required for iridophore differentiation in zebrafish. Wild-type Ltk crispants present similar pigment phenotypes to melanoid, strongly implicating Ltk as the melanoid locus. In concert with recent findings in zebrafish, our results support the idea of direct fate specification of pigment cells and, more generally, the single-origin hypothesis of pigment cell development.

Vitamin D regulates ion regulation by affecting the ionocyte differentiation in zebrafish (Danio rerio) larvae

Freshwater teleosts frequently face the stress of varied ion and pH levels; therefore, they have developed related defense mechanisms to maintain the homeostasis of body-fluid ion and acid-base balance. The different subtypes of ionocytes expressed in the branchial epithelium of adult fish or the skin of larvae are the major sites for fish ion regulation. 1α,25-dihydroxyvitamin D3 (1α,25(OH)2D3), the bioactive form of vitamin D, is a steroid hormone that is involved in the regulation of Ca2+ uptake and acid secretion in teleosts. Our results revealed that 1α,25(OH)2D3 levels were not changed in zebrafish larvae upon exposure to low Na+ freshwater compared to normal freshwater. In contrast, 1α,25(OH)2D3 levels were substantially higher in fish exposed to acidic and low Ca2+ freshwater than in those exposed to normal freshwater. Some hormones regulate ion regulation and acid secretion by modulating ionocyte differentiation and/or proliferation in teleosts; however, the role of vitamin D in this process is unclear. Zebrafish larvae were used as a model in the present study to explore the effect of vitamin D on ionocyte proliferation and/or differentiation. The present study indicated that 1α,25(OH)2D3 treatment increased the number of foxi3a-positive cells, ionocyte progenitors, and mature ionocytes. However, the number of P63-positive epidermal stem cells did not change in the zebrafish larvae treated with 1α,25(OH)2D3. These results revealed that vitamin D exerts a positive effect on the number of ionocytes by increasing the differentiation of ionocytes. Increased ionocyte differentiation by vitamin D is suggested to elevate the capacity of ion regulation and acid secretion in zebrafish to cope with external stress. The present findings indicate the role of vitamin D in the regulation of ionocyte differentiation and provide new insights into the mechanisms of stress adaptation of fish.

The isl2a transcription factor regulates pituitary development in zebrafish.

ISL LIM homeobox 2, also known as insulin gene enhancer protein ISL-2 (ISL2), is a transcription factor gene that participates in a wide range of developmental events. However, the role of ISL2 in the hypothalamus-pituitary-thyroid axis is largely unknown. In the present study, we characterized the expression patterns of ISL2 and revealed its regulative role during embryogenesis using zebrafish. We used the CRISPR/Cas9 system to successfully establish homozygous ISL2-orthologue (isl2a and isl2b) knockout zebrafish. Moreover, we utilized these knockout zebrafish to analyze the pituitary and thyroid phenotypes in vivo. For further molecular characterization, in situ hybridization and immunofluorescence were performed. The isl2a mutant zebrafish presented with thyroid hypoplasia, reduced whole-body levels of thyroid hormones, increased early mortality, gender imbalance, and morphological retardation during maturity. Additionally, thyrotropes, a pituitary cell type, was notably decreased during development. Importantly, the transcriptional levels of pituitary-thyroid axis hormones-encoding genes, such as tshba, cga, and tg, were significantly decreased in isl2a mutants. Finally, the thyroid dysplasia in isl2a mutant larvae may be attributed to a reduction in proliferation rather than changes in apoptosis. In summary, isl2a regulates the transcriptional levels of marker genes in hypothalamus-pituitary-thyroid axis, and isl2a knockout causing low thyroid hormone levels in zebrafish. Thus, isl2a identified by the present study, is a novel regulator for pituitary cell differentiation in zebrafish, resulting in thyroid gland hypoplasia and phenotypes of hypothyroidism.

Stage-specific and cell type-specific requirements of ikzf1 during haematopoietic differentiation in zebrafish

The zinc finger transcription factor Ikaros1 (Ikzf1) is required for lymphoid development in mammals. Four zinc fingers constitute its DNA binding domain and two zinc fingers are present in the C-terminal protein interaction module. We describe the phenotypes of zebrafish homozygous for two distinct mutant ikzf1 alleles. The IT325 variant lacks the C-terminal two zinc fingers, whereas the fr105 variant retains only the first zinc finger of the DNA binding domain. An intact ikzf1 gene is required for larval T cell development, whereas low levels of adult lymphoid development recover in the mutants. By contrast, the mutants exhibit a signature of increased myelopoiesis at larval and adult stages. Both mutations stimulate erythroid differentiation in larvae, indicating that the C-terminal zinc fingers negatively regulate the extent of red blood cell production. An unexpected differential effect of the two mutants on adult erythropoiesis suggests a direct requirement of an intact DNA binding domain for entry of progenitors into the red blood cell lineage. Collectively, our results reinforce the biological differences between larval and adult haematopoiesis, indicate a stage-specific function of ikzf1 in regulating the hierarchical bifurcations of differentiation, and assign distinct functions to the DNA binding domain and the C-terminal zinc fingers.

A landscape of differentiated biological processes involved in the initiation of sex differentiation in zebrafish

Zebrafish (Danio rerio) has been used as a promising animal model to study gonadal development and gametogenesis. Although previous studies have identified critical molecules participating in zebrafish gonad differentiation, a landscape view of the biological processes involved in this process is still lacking. Here we isolated intact zebrafish differentiating gonads, at 25 days post-fertilization (dpf) and 30 dpf and conducted RNA-seq analyses on the juvenile gonads that tended to develop into ovaries or testes. Our study demonstrates that the juvenile ovary and testis at 25 dpf and 30 dpf are different at the biological process level. During ovary differentiation, the biological processes related to metabolic activities in the production of energy and maternal substances, RNA degradation, and DNA repair were enriched. During testis differentiation, the biological processes related to cell proliferation, differentiation, and morphogenesis were enriched, with a total of 15 signaling pathways. Notably, we reveal that the immune-related processes are extensively involved in the regulation of testis development. Overall, this study provides a landscape of differentiated biological processes and novel insights into the initiation of sex differentiation in zebrafish.

Using Single Cell Omics to Understand Posterior Neural Crest Differentiation in Zebrafish

Neural crest cells (NCCs) are a migratory, multipotent, and transient stem cell population essential to vertebrate development. Born along the neural tube, NCCs are regionally specified along the anterior‐posterior axis to contribute to a vast number of cell lineages; such as craniofacial cartilage, pigment cells, and large aspects of the peripheral nervous system. In our recently published single‐cell RNA‐sequencing (scRNA‐seq) analysis of posterior NCC lineages in zebrafish embryos, we discovered a combinatorial expression code of transcripts that encode for posterior Hox transcription factors, which was especially specific to NCC‐derived neural lineages. While Hox expression combinations had been characterized within pre‐migratory and recently delaminated cranial NCCs in prior studies, such a high‐resolution signature had not been described within the context of posterior NCC lineages during their differentiation. Here, we report a lifetime atlas of NCC and NCC‐derived cells along the posterior zebrafish body, greatly expanding our prior analyses. To achieve this, several published and publicly‐available zebrafish scRNA‐seq studies characterizing NCC and NCC‐derived cells were computationally merged with our prior data sets and integrated, encompassing: 24 hours post fertilization (hpf), 48‐50 hpf, 68‐70 hpf, 5 days post fertilization (dpf), juvenile, and adult tissues—producing a cohesive distribution of cell identities. Within the lifetime atlas, various Hox gene transcript signatures were observed where numerous discrete cell lineages formed over time and unique Hox expression codes resolved within these incipient lineages. The generation of a multi‐stage comprehensive transcriptional atlas from multiple scRNA‐seq studies will serve as a valuable tool for the NCC and zebrafish communities. Additionally, this atlas of cell lineages demonstrates the complexity and robustness of combinatorial Hox expression codes in posterior NCC throughout development and suggests possible functional roles therein.

Zebrafish Establish Female Germ Cell Identity by Advancing Cell Proliferation and Meiosis.

Zebrafish is a popular research model; but its mechanism of sex determination is unclear and the sex of juvenile fish cannot be distinguished. To obtain fish with defined sex, we crossed domesticated zebrafish with the Nadia strain that has a female-dominant W segment. These fish were placed on a ziwi:GFP background to facilitate sorting of fluorescent germ cells for transcriptomic analysis. We analyzed the transcriptomes of germ cells at 10–14 days postfertilization (dpf), when sex dimorphic changes started to appear. Gene ontology showed that genes upregulated in the 10-dpf presumptive females are involved in cell cycles. This correlates with our detection of increased germ cell numbers and proliferation. We also detected upregulation of meiotic genes in the presumptive females at 14 dpf. Disruption of a meiotic gene, sycp3, resulted in sex reversal to infertile males. The germ cells of sycp3 mutants could not reach diplotene and underwent apoptosis. Preventing apoptosis by disrupting tp53 restored female characteristics in sycp3 mutants, demonstrating that adequate germ cells are required for female development. Thus, our transcriptome and gene mutation demonstrate that initial germ cell proliferation followed by meiosis is the hallmark of female differentiation in zebrafish.

Perspective on automated in vivo drug screening using the chodl mutant zebrafish line.

Perspective on automated in vivo drug screening using the chodl mutant zebrafish line.

Cyp11a2 Is Essential for Oocyte Development and Spermatogonial Stem Cell Differentiation in Zebrafish.

Cytochrome P45011A1, encoded by Cyp11a1, converts cholesterol to pregnenolone (P5), the first and rate-limiting step in steroidogenesis. In zebrafish, cyp11a1 is maternally expressed and cyp11a2 is considered the ortholog of Cyp11a1 in mammals. A recent study has shown that depletion of cyp11a2 resulted in steroidogenic deficiencies and the mutants developed into males with feminized secondary sexual characteristics. Here, we independently generated cyp11a2 mutants in zebrafish and showed that the mutants can develop into males and females in the juvenile stage, but finally into infertile males with defective mating behavior in the adult stage. In the developing ovaries, the cyp11a2 mutation led to stage I oocyte apoptosis and final sex reversal, which could be partially rescued by treatment with P5 but not estradiol. In the developing testes, depletion of cyp11a2 resulted in dysfunction of Sertoli cells and lack of functional Leydig cells. Spermatogonial stem cells (SSCs) in the mutant testes underwent active self-renewal but no differentiation, resulting in a high abundance of SSCs in the testis, as revealed by immunofluorescence staining with Nanos2 antibody. The high abundance and differentiation competence of SSCs in the mutant testes were verified by a novel testicular cell transplantation method developed in this study, by transplanting mutant testicular cells into germline-depleted wild-type (WT) fish. The transplanted mutant SSCs efficiently differentiated into functional spermatids in WT hosts. Overall, our study demonstrates the functional importance of cyp11a2 in early oogenesis and differentiation of SSCs.

Loss of Polycomb Repressive Complex 2 Function Alters Digestive Organ Homeostasis and Neuronal Differentiation in Zebrafish.

Polycomb repressive complex 2 (PRC2) mediates histone H3K27me3 methylation and the stable transcriptional repression of a number of gene expression programs involved in the control of cellular identity during development and differentiation. Here, we report on the generation and on the characterization of a zebrafish line harboring a null allele of eed, a gene coding for an essential component of the PRC2. Homozygous eed-deficient mutants present a normal body plan development but display strong defects at the level of the digestive organs, such as reduced size of the pancreas, hepatic steatosis, and a loss of the intestinal structures, to die finally at around 10–12 days post fertilization. In addition, we found that PRC2 loss of function impairs neuronal differentiation in very specific and discrete areas of the brain and increases larval activity in locomotor assays. Our work highlights that zebrafish is a suited model to study human pathologies associated with PRC2 loss of function and H3K27me3 decrease.

Noncanonical protease-activated receptor 1 regulates lymphatic differentiation in zebrafish

SummaryThe differentiation of lymphatic progenitors is a crucial step in lymphangiogenesis. However, its underlying mechanism remains unclear. Here, we found that noncanonical protease-activated receptor 1 (par1) regulates the differentiation of lymphatic progenitors in zebrafish embryos. Loss of par1 function impaired lymphatic differentiation by downregulating prox1a expression in parachordal lymphangioblasts and caused compromised thoracic duct formation in zebrafish. Meanwhile, the G protein gnai2a, a par1 downstream effector, was selectively required for lymphatic development in zebrafish, and its mutation mimicked the lymphatic phenotype observed in par1 mutants. Interestingly, mmp13, but not thrombin, was required for lymphatic development in zebrafish. Furthermore, analyses of genetic interactions confirmed that mmp13b serves as a par1 upstream protease to regulate lymphatic development in zebrafish embryos. Mechanistically, par1 promotes flt4 expression and phospho-Erk1/2 activity in the posterior cardinal vein. Taken together, our findings highlight a function of par1 in the regulation of lymphatic differentiation in zebrafish embryos.

17β-Trenbolone binds to androgen receptor, decreases number of primordial germ cells, modulates expression of genes related to sexual differentiation, and affects sexual differentiation in zebrafish (Danio rerio)

Exposure to 17β-trenbolone caused a skewed sex ratio in fish. However, the molecular initiating event and key molecular event(s) remain unknown. In this study, zebrafish were exposed to 17β-trenbolone at nominal concentrations of 2 ng/L, 20 ng/L, 200 ng/L, and 2000 ng/L from fertilization to 60 days post fertilization (dpf). First, the sex ratio at 60 dpf was calculated to evaluate adverse outcomes on sexual differentiation. 17β-Trenbolone caused a skewed sex ratio toward males, with intersex individuals observed in the 20 ng/L group and all-male populations found in the 200 ng/L and 2000 ng/L groups. Then, the distribution and number of primordial germ cells, the expression of sex differentiation-related genes, and plasma vitellogenin concentrations were detected in wild-type zebrafish and the EGFP-nanos-3′UTR transgenic line using whole-mount in situ hybridization, real-time PCR, EGFP fluorescence quantification, and enzyme-linked immunosorbent assay. The results indicated that 17β-trenbolone exposure decreased the number of primordial germ cells at 1 dpf and 3 dpf, decreased expression of ovarian differentiation-related genes foxl2 and cyp19a1a at 60 dpf, increased expression of testis differentiation-related genes dmrt1, sox9a, and amh at 60 dpf, and decreased plasma vitellogenin levels at 60 dpf, revealing the key molecular events at different time points involved in affected sexual differentiation by 17β-trenbolone. Finally, molecular docking showed that 17β-trenbolone docked into ligand-binding domain of zebrafish androgen receptor with high binding energy (−3.72 kcal/mol), suggesting that binding to androgen receptor is the molecular initiating event affecting sexual differentiation by 17β-trenbolone. We found that 17β-trenbolone can bind to the zebrafish androgen receptor, decrease the number of primordial germ cells during the early embryonic stage, modulate the expression of genes related to sexual differentiation during gonadal differentiation, and eventually cause a skewed sex ratio toward males in adults.