nr0b1 (DAX1) mutation in zebrafish causes female-to-male sex reversal through abnormal gonadal proliferation and differentiation

Tilapia, a good model for studying reproductive endocrinology

The Nile tilapia (Oreochromis niloticus), with a system of XX/XY sex determination, is a worldwide farmed fish with a shorter sexual maturation time than that of most cultured fish. Tilapia show a spawning cycle of approximately 14 days and can be artificially propagated in the laboratory all year round to obtain genetically all female (XX) and all male (XY) fry. Its genome sequence has been opened, and a perfect gene editing platform has been established. With a moderate body size, it is convenient for taking enough blood to measure hormone level. In recent years, using tilapia as animal model, we have confirmed that estrogen is crucial for female development because 1) mutation of star2, cyp17a1 or cyp19a1a (encoding aromatase, the key enzyme for estrogen synthesis) results in sex reversal (SR) due to estrogen deficiency in XX tilapia, while mutation of star1, cyp11a1, cyp17a2, cyp19a1b or cyp11c1 affects fertility due to abnormal androgen, cortisol and DHP levels in XY tilapia; 2) when the estrogen receptors (esr2a/esr2b) are mutated, the sex is reversed from female to male, while when the androgen receptors are mutated, the sex cannot be reversed; 3) the differentiated ovary can be transdifferentiated into functional testis by inhibition of estrogen synthesis, and the differentiated testis can be transdifferentiated into ovary by simultaneous addition of exogenous estrogen and androgen synthase inhibitor; 4) loss of male pathway genes amhy, dmrt1, gsdf causes SR with upregulation of cyp19a1a in XY tilapia. Disruption of estrogen synthesis rescues the male to female SR of amhy and gsdf but not dmrt1 mutants; 5) mutation of female pathway genes foxl2 and sf-1 causes SR with downregulation of cyp19a1a in XX tilapia; 6) the germ cell SR of foxl3 mutants fails to be rescued by estrogen treatment, indicating that estrogen determines female germ cell fate through foxl3. This review also summarized the effects of deficiency of other steroid hormones, such as androgen, DHP and cortisol, on fish reproduction. Overall, these studies demonstrate that tilapia is an excellent animal model for studying reproductive endocrinology of fish.

Nuclear receptor superfamily structural diversity in pacific oyster: In silico identification of estradiol binding candidates

The increasing presence of anthropogenic contaminants in aquatic environments poses challenges for species inhabiting contaminated sites. Due to their structural binding characteristics to ligands that inhibit or activate gene transcription, these xenobiotic compounds frequently target the nuclear receptor superfamily. The present work aims to understand the potential interaction between the hormone 17-β-estradiol, an environmental contaminant, and the nuclear receptors of Crassostrea gigas, the Pacific oyster. This filter-feeding, sessile oyster species is subject to environmental changes and exposure to contaminants. In the Pacific oyster, the estrogen-binding nuclear receptor is not able to bind this hormone as it does in vertebrates. However, another receptor may exhibit responsiveness to estrogen-like molecules and derivatives. We employed high-performance in silico methodologies, including three-dimensional modeling, molecular docking and atomistic molecular dynamics to identify likely binding candidates with the target moecule. Our approach revealed that among the C. gigas nuclear receptor superfamily, candidates with the most favorable interaction with the molecule of interest belonged to the NR1D, NR1H, NR1P, NR2E, NHR42, and NR0B groups. Interestingly, NR1H and NR0B were associated with planktonic/larval life cycle stages, while NR1P, NR2E, and NR0B were associated with sessile/adult life stages. The application of this computational methodological strategy demonstrated high performance in the virtual screening of candidates for binding with the target xenobiotic molecule and can be employed in other studies in the field of ecotoxicology in non-model organisms.

De novo assembly, characterization and comparative transcriptome analysis of gonads reveals sex-biased genes in Coreoperca whiteheadi

The wild Coreoperca whiteheadi is considered as the primordial species in sinipercine fish, which has valuable genetic information. Unfortunately, C. whiteheadi was listed as a near-threatened species because of the environmental pollution, over-exploitation and species invasion. Therefore, more genetic information is needed to have a better understanding of gonadal development in C. whiteheadi. Here, the first gonadal transcriptomes analysis of C. whiteheadi was conducted and 277.14 million clean reads were generated. A total of 96,753 unigenes were successfully annotated. By comparing ovary and testis transcriptomes, a total of 21,741 differentially expressed genes (DEGs) were identified, of which 12,057 were upregulated and 9684 were downregulated in testes. Among them, we also identified about 53 differentially expressed sex-biased genes. Subsequently, the expression of twenty-four DEGs were confirmed by real-time fluorescence quantitative PCR. Furthermore, the histological analysis was conducted on ovaries and testes of one-year-old C. whiteheadi. Our results provided basic support for further studies on the function of sex-biased genes and the molecular mechanism of sex determination and reproduction in C. whiteheadi.

De Novo Assembly, Characterization and Comparative Transcriptome Analysis of the Gonads of Jade Perch (Scortum barcoo)

Due to the high meat yield and rich nutritional content, jade perch (Scortum barcoo) has become an important commercial aquaculture species in China. Jade perch has a slow growth rate, taking 3-4 years to reach sexual maturity, and has almost no difference in body size between males and females. However, the study of its gonad development and reproduction regulation is still blank, which limited the yield increase. Herein, the gonad transcriptomes of juvenile males and females of S. barcoo were identified for the first time. A total of 107,060 unigenes were successfully annotated. By comparing male and female gonad transcriptomes, a total of 23,849 differentially expressed genes (DEGs) were identified, of which 9517 were downregulated, and 14,332 were upregulated in the testis. In addition, a large number of DEGs involved in sex differentiation, gonadal development and differentiation and gametogenesis were identified, and the differential expression patterns of some genes were further verified using real-time fluorescence quantitative PCR. The results of this study will provide a valuable resource for further studies on sex determination and gonadal development of S. barcoo.

Sex blind: bridging the gap between drug exposure and sex-related gene expression in Danio rerio using next-generation sequencing (NGS) data and a literature review to find the missing links in pharmaceutical and environmental toxicology studies

The sex of both humans and Danio rerio has previously been shown to affect the way individuals respond to drug exposure. Genes which allow identification of sex in juvenile zebrafish show potential to reveal these confounding variables between sex in toxicological and preclinical trials but the link between these is so far missing. These sex-specific, early expressed genes where expression is not altered by drug exposure must be carefully selected for this purpose. We aimed to discover genes which can be used in pharmaceutical trials and environmental toxicology studies to uncover sex-related variations in gene expression with drug application using the model organism Danio rerio. Previously published early sex determining genes from King et al. were evaluated as well as additional genes selected from our zebrafish Next-generation sequencing (NGS) data which are known from previously published works not to be susceptible to changes in expression with drug exposure. NGS revealed a further ten female-specific genes (vtg1, cyp17a1, cyp19a1a, igf3, ftz-f1, gdf9, foxl2a, Nr0b1, ipo4, lhcgr) and five male related candidate genes (FKBP5, apobb1, hbaa1, dmrt1, spata6) which are also expressed in juvenile zebrafish, 28 days post fertilisation (dpf). Following this, a literature review was performed to classify which of these early-expressed sex specific genes are already known to be affected by drug exposure in order to determine candidate genes to be used in pharmaceutical trials or environmental toxicology testing studies. Discovery of these early sex-determining genes in Danio rerio will allow identification of sex-related responses to drug testing to improve sex-specific healthcare and the medical treatment of human patients.

Becoming female: Ovarian differentiation from an evolutionary perspective

Differentiation of the bipotential gonadal primordium into ovaries and testes is a common process among vertebrate species. While vertebrate ovaries eventually share the same functions of producing oocytes and estrogens, ovarian differentiation relies on different morphogenetic, cellular, and molecular cues depending on species. The aim of this review is to highlight the conserved and divergent features of ovarian differentiation through an evolutionary perspective. From teleosts to mammals, each clade or species has a different story to tell. For this purpose, this review focuses on three specific aspects of ovarian differentiation: ovarian morphogenesis, the evolution of the role of estrogens on ovarian differentiation and the molecular pathways involved in granulosa cell determination and maintenance.

Expression and Transcript Localization of star, sf-1, and dax-1 in the Early Brain of the Orange-Spotted Grouper Epinephelus coioides

We investigated the developmental expression and localization of sf-1 and dax-1 transcripts in the brain of the juvenile orange-spotted grouper in response to steroidogenic enzyme gene at various developmental ages in relation to gonadal sex differentiation. The sf-1 transcripts were significantly higher from 110-dah (day after hatching) and gradually increased up to 150-dah. The dax-1 mRNA, on the other hand, showed a decreased expression during this period, in contrast to sf-1 expression. At the same time, the early brain had increased levels of steroidogenic gene (star). sf-1 and star hybridization signals were found to be increased in the ventromedial hypothalamus at 110-dah; however, dax-1 mRNA signals decreased in the early brain toward 150-dah. Furthermore, the exogenous estradiol upregulated star and sf-1 transcripts in the early brain of the grouper. These findings suggest that sf-1 and dax-1 may have an antagonistic expression pattern in the early brain during gonadal sex differentiation. Increased expression of steroidogenic gene together with sf-1 during gonadal differentiation strongly suggests that sf-1 may play an important role in the juvenile grouper brain steroidogenesis and brain development.

Sexual determination in zebrafish

Zebrafish have emerged as a major model organism to study vertebrate reproduction due to their high fecundity and external development of eggs and embryos. The mechanisms through which zebrafish determine their sex have come under extensive investigation, as they lack a definite sex-determining chromosome and appear to have a highly complex method of sex determination. Single-gene mutagenesis has been employed to isolate the function of genes that determine zebrafish sex and regulate sex-specific differentiation, and to explore the interactions of genes that promote female or male sexual fate. In this review, we focus on recent advances in understanding of the mechanisms, including genetic and environmental factors, governing zebrafish sex development with comparisons to gene functions in other species to highlight conserved and potentially species-specific mechanisms for specifying and maintaining sexual fate.

nr0b1 (DAX1) loss of function in zebrafish causes hypothalamic defects via abnormal progenitor proliferation and differentiation

The nuclear receptor DAX-1 (encoded by the NR0B1 gene) is presented in the hypothalamic tissues in humans and other vertebrates. Human patients with NR0B1 mutations often have hypothalamic-pituitary defects, but the involvement of NR0B1 in hypothalamic development and function is not well understood. Here, we report the disruption of the nr0b1 gene in zebrafish causes abnormal expression of gonadotropins, a reduction in fertilization rate, and an increase in post-fasting food intake, which is indicative of abnormal hypothalamic functions. We find that loss of nr0b1 increases the number of prodynorphin (pdyn)-expressing neurons but decreases the number of pro-opiomelanocortin (pomcb)-expressing neurons in the zebrafish hypothalamic arcuate region (ARC). Further examination reveals that the proliferation of progenitor cells is reduced in the hypothalamus of nr0b1 mutant embryos accompanying with the decreased expression of genes in the Notch signaling pathway. Additionally, the inhibition of Notch signaling in wild-type embryos increases the number of pdyn neurons, mimicking the nr0b1 mutant phenotype. In contrast, ectopic activation of Notch signaling in nr0b1 mutant embryos decreases the number of pdyn neurons. Taken together, our results suggest that nr0b1 regulates neural progenitor proliferation and maintenance to ensure normal hypothalamic neuron development.

Nonsense variant of NR0B1 causes hormone disorders associated with congenital adrenal hyperplasia

Congenital adrenal hyperplasia (CAH) is a rare X-linked recessive inherited disease that is considered a major cause of steroidogenesis disorder and is associated with variants or complete deletion of the NR0B1 gene. The DAX-1 protein (encoded by NR0B1) is a vertebrate-specific orphan nuclear receptor and is also a transcriptional factor for adrenal and reproductive development. CAH usually causes adrenal insufficiency in infancy and early childhood, leading to hypogonadotropic hypogonadism in adulthood; however, few adult cases have been reported to date. In this study, we examined a Chinese family with one adult patient with CAH, and identified a putative variant of NR0B1 gene via next-generation sequencing (NGS), which was confirmed with Sanger sequencing. A novel nonsense variant (c.265C>T) was identified in the NR0B1 gene, which caused the premature termination of DAX-1 at residue 89 (p.G89*). Furthermore, mutant NR0B1 gene displayed a partial DAX-1 function, which may explain the late pathogenesis in our case. Additionally, qPCR revealed the abnormal expression of four important genes identified from ChIP-seq, which were associated with energy homeostasis and steroidogenesis, and were influenced by the DAX-1 mutant. In addition, hormone disorders can be caused by DAX-1 mutant and partially recovered by siRNA of PPARGC1A. Herein, we identified a novel nonsense variant (c.265C>T) of NR0B1 in a 24-year-old Chinese male who was suffering from CAH. This mutant DAX-1 protein was found to have disordered energy homeostasis and steroidogenesis based on in vitro studies, which was clinically consistent with the patient’s phenotypic features.

Characterization, expression, and regulatory effects of nr0b1a and nr0b1b in spotted scat (Scatophagus argus)

Nuclear receptor subfamily 0 group B member 1 (Nr0b1) belongs to the nuclear receptor (NR) superfamily. It plays critical roles in sex determination, sex differentiation, and gonadal development in mammals. In this study, the duplicated genes nr0b1a and nr0b1b were identified in spotted scat (Scatophagus argus). Phylogenetic and synteny analyses revealed that, unlike nr0b1a, nr0b1b was retained in several species of teleosts after an nr0b1 gene duplication event but was secondarily lost in other fish species, amphibians, reptiles, birds, and mammals. In a sequence analysis, only 1.5 LXXLL-related repeat motifs were identified in spotted scat Nr0b1a, Nr0b1b, and non-mammalian Nr0b1a/Nr0b1, different from the 3.5 repeat motifs in mammalian Nr0b1. By qPCR, nr0b1a and nr0b1b were highly expressed in testes from stages IV to V and in ovaries from stages II to IV, respectively. Male-to-female sex reversal was induced in XY spotted scat by the administration of exogenous E2. A qPCR analysis showed that nr0b1b mRNA expression was higher in sex-reversed XY fish than in control XY fish, with no difference in nr0b1a. A luciferase assay showed that spotted scat Nr0b1a and Nr0b1b did not individually activate cyp19a1a gene transcription. As in mammals, spotted scat Nr0b1a suppressed Nr5a1-mediated cyp19a1a expression, despite containing only 1.5 LXXLL-related repeat motifs in its N-terminal region, while Nr0b1b stimulated Nr5a1-mediated cyp19a1a transcription. These results demonstrated that nr0b1a and nr0b1b in spotted scat have distinct expression patterns and regulatory effects and further indicate that nr0b1b might be involved in ovarian development by regulating Nr5a1-mediated cyp19a1a expression.

A fish with no sex: gonadal and adrenal functions partition between zebrafish NR5A1 co-orthologs

People with NR5A1 mutations experience testicular dysgenesis, ovotestes, or adrenal insufficiency, but we do not completely understand the origin of this phenotypic diversity. NR5A1 is expressed in gonadal soma precursor cells before expression of the sex-determining gene SRY. Many fish have two co-orthologs of NR5A1 that likely partitioned ancestral gene subfunctions between them. To explore ancestral roles of NR5A1, we knocked out nr5a1a and nr5a1b in zebrafish. Single-cell RNA-seq identified nr5a1a-expressing cells that co-expressed genes for steroid biosynthesis and the chemokine receptor Cxcl12a in 1-day postfertilization (dpf) embryos, as does the mammalian adrenal-gonadal (interrenal-gonadal) primordium. In 2dpf embryos, nr5a1a was expressed stronger in the interrenal-gonadal primordium than in the early hypothalamus but nr5a1b showed the reverse. Adult Leydig cells expressed both ohnologs and granulosa cells expressed nr5a1a stronger than nr5a1b. Mutants for nr5a1a lacked the interrenal, formed incompletely differentiated testes, had no Leydig cells, and grew far larger than normal fish. Mutants for nr5a1b formed a disorganized interrenal and their gonads completely disappeared. All homozygous mutant genotypes lacked secondary sex characteristics, including male breeding tubercles and female sex papillae, and had exceedingly low levels of estradiol, 11-ketotestosterone, and cortisol. RNA-seq showed that at 21dpf, some animals were developing as females and others were not, independent of nr5a1 genotype. By 35dpf, all mutant genotypes greatly under-expressed ovary-biased genes. Because adult nr5a1a mutants form gonads but lack an interrenal and conversely, adult nr5a1b mutants lack a gonad but have an interrenal, the adrenal, and gonadal functions of the ancestral nr5a1 gene partitioned between ohnologs after the teleost genome duplication, likely owing to reciprocal loss of ancestral tissue-specific regulatory elements. Identifying such elements could provide hints to otherwise unexplained cases of Differences in Sex Development.

Igf3 is essential for ovary differentiation in zebrafish†

Zebrafish gonadal sexual differentiation is an important but poorly understood subject. Previously, we have identified a novel insulin-like growth factor (Igf) named insulin-like growth factor 3 (Igf3) in teleosts. The importance of Igf3 in oocyte maturation and ovulation has been recently demonstrated by us in zebrafish. In this study, we have further found the essential role of Igf3 in gonadal sexual differentiation of zebrafish. A differential expression pattern of igf3 between ovary and testis during sex differentiation (higher level in ovary than in testis) was found in zebrafish. An igf3 knockout zebrafish line was established using TALENs-mediated gene knockout technique. Intriguingly, all igf3 homozygous mutants were males due to the female-to-male sex reversal occurred during sex differentiation. Further analysis showed that Igf3 did not seem to affect the formation of so-called juvenile ovary and oocyte-like germ cells. Oocyte development was arrested at primary growth stage, and the ovary was gradually sex-reversed to testis before 60 day post fertilization (dpf). Such sex reversal was likely due to decreased germ cell proliferation by suppressing PI3K/Akt pathway in early ovaries of igf3 mutants. Estrogen is considered as a master regulator in fish sex differentiation. Here, we found that igf3 expression could be upregulated by estrogen in early stages of ovarian follicles as evidenced in in vitro treatment assays and cyp19a1a mutant zebrafish, and E2 failed to rescue the defects of igf3 mutants in ovarian development, suggesting that Igf3 may serve as a downstream factor of estrogen signaling in sex differentiation. Taken together, we demonstrated that Igf3 is essential for ovary differentiation in zebrafish.

Sex determination, gonadal sex differentiation, and plasticity in vertebrate species

A diverse array of sex determination (SD) mechanisms, encompassing environmental to genetic, have been found to exist among vertebrates, covering a spectrum from fixed SD mechanisms (mammals) to functional sex change in fishes (sequential hermaphroditic fishes). A major landmark in vertebrate SD was the discovery of the SRY gene in 1990. Since that time, many attempts to clone an SRY ortholog from nonmammalian vertebrates remained unsuccessful, until 2002, when DMY/dmrt1by was discovered as the SD gene of a small fish, medaka. Surprisingly, however, DMY/dmrt1by was found in only 2 species among more than 20 species of medaka, suggesting a large diversity of SD genes among vertebrates. Considerable progress has been made over the last 3 decades, such that it is now possible to formulate reasonable paradigms of how SD and gonadal sex differentiation may work in some model vertebrate species. This review outlines our current understanding of vertebrate SD and gonadal sex differentiation, with a focus on the molecular and cellular mechanisms involved. An impressive number of genes and factors have been discovered that play important roles in testicular and ovarian differentiation. An antagonism between the male and female pathway genes exists in gonads during both sex differentiation and, surprisingly, even as adults, suggesting that, in addition to sex-changing fishes, gonochoristic vertebrates including mice maintain some degree of gonadal sexual plasticity into adulthood. Importantly, a review of various SD mechanisms among vertebrates suggests that this is the ideal biological event that can make us understand the evolutionary conundrums underlying speciation and species diversity.

TCDD-induced multi- and transgenerational changes in the methylome of male zebrafish gonads

The legacy endocrine disrupting chemical and aryl hydrocarbon receptor agonist, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), is produced as a byproduct of industrial processes and causes adverse health effects ranging from skin irritation to cancer. TCDD endpoints are also observed in subsequent, unexposed generations; however, the mechanisms of these multi- and transgenerational effects are unknown. We hypothesized an epigenetic mechanism, specifically DNA methylation for the transgenerational, male-mediated reproductive effects of developmental TCDD exposure. Using whole genome bisulfite sequencing, we evaluated DNA methylation changes in three generations of zebrafish, the first of which was exposed to TCDD during sexual development at 50 ppt for 1 h at both 3- and 7-week post-fertilization. We discovered that TCDD induces multi- and transgenerational methylomic changes in testicular tissue from zebrafish with decreased reproductive capacity, but most significantly in the indirectly exposed F1 generation. In comparing differentially methylated genes to concurrent transcriptomic changes, we identified several genes and pathways through which transgenerational effects of low level TCDD exposure are likely inherited. These include significant differential methylation of genes involved in reproduction, endocrine function, xenobiotic metabolism, and epigenetic processing. Notably, a number of histone modification genes were both differentially methylated and expressed in all generations, and many differentially methylated genes overlapped between multiple generations. Collectively, our results suggest that DNA methylation is a promising mechanism to explain male-mediated transgenerational reproductive effects of TCDD exposure in zebrafish, and these effects are likely inherited through integration of multiple epigenetic pathways.

Molecular and morphological sex differentiation in sablefish (Anoplopoma fimbria), a marine teleost with XX/XY sex determination

Phenotypic sex of an organism is determined by molecular changes in the gonads, so-called molecular sex differentiation, which should precede the rise of cellular or anatomical sex-distinguishing features. This study characterized molecular and morphological sex differentiation in sablefish (Anoplopoma fimbria), a marine teleost with established XX/XY genotypic sex determination. Next generation sequencing was conducted on sablefish ovarian and testicular mRNAs to obtain sequences for transcripts associated with vertebrate sex determination and differentiation and early reproductive development. Gene-specific PCRs were developed to determine the distribution and ontogenetic gonadal expression of transcription, growth, steroidogenic and germline factors, as well as gonadotropin and steroid receptors. Molecular changes associated with sex differentiation were first apparent in both XY- and XX-genotype sablefish at ~ 60 mm in body length and prior to histological signs of sex differentiation. The earliest and most robust markers of testicular differentiation were gsdf, amh, dmrt1, cyp11b, star, sox9a, and fshr. Markedly elevated mRNA levels of several steroidogenesis-related genes and ar2 in differentiating testes suggested that androgens play a role in sablefish testicular differentiation. The earliest markers of ovarian differentiation were cyp19a1a, lhcgr, foxl2, nr0b1, and igf3. Other transcripts such as figla, zp3, and pou5f3 were expressed predominantly in XX-genotype fish and significantly increased with the first appearance and subsequent development of primary oocytes. This study provides valuable insight to the developmental sequence of events associated with gonadal sex differentiation in marine teleosts with XX/XY sex determination. It also implicates particular genes in processes of male and female development and establishes robust molecular markers for phenotypic sex in sablefish, useful for ongoing work related to sex control and reproductive sterilization.

Sex-Dependent RNA Editing and N6-adenosine RNA Methylation Profiling in the Gonads of a Fish, the Olive Flounder (Paralichthys olivaceus)

Adenosine-to-inosine (A-to-I) editing and N6-methyladenosine (m6A) are two of the most abundant RNA modifications. Here, we examined the characteristics of the RNA editing and transcriptome-wide m6A modification profile in the gonads of the olive flounder, Paralichthys olivaceus, an important maricultured fish in Asia. The gonadal differentiation and development of the flounder are controlled by genetic as well as environmental factors, and the epigenetic mechanism may play an important role. In total, 742 RNA editing events were identified, 459 of which caused A to I conversion. Most A-to-I sites were located in 3′UTRs, while 61 were detected in coding regions (CDs). The number of editing sites in the testis was higher than that in the ovary. Transcriptome-wide analyses showed that more than one-half of the transcribed genes presented an m6A modification in the flounder gonads, and approximately 60% of the differentially expressed genes (DEGs) between the testis and ovary appeared to be negatively correlated with m6A methylation enrichment. Further analyses revealed that the mRNA expression of some sex-related genes (e.g., dmrt1 and amh) in the gonads may be regulated by changes in mRNA m6A enrichment. Functional enrichment analysis indicated that the RNA editing and m6A modifications were enriched in several canonical pathways (e.g., Wnt and MAPK signaling pathways) in fish gonads and in some pathways whose roles have not been investigated in relation to fish sex differentiation and gonadal development (e.g., PPAR and RNA degradation pathways). There were 125 genes that were modified by both A-to-I editing and m6A, but the two types of modifications mostly occurred at different sites. Our results suggested that the presence of sex-specific RNA modifications may be involved in the regulation of gonadal development and gametogenesis.

Zebrafish as a model for studying ovarian development: Recent advances from targeted gene knockout studies

Ovarian development is a complex process controlled by precise coordination of multiple factors. The targeted gene knockout technique is a powerful tool to study the functions of these factors. The successful application of this technique in mice in the past three decades has significantly enhanced our understanding on the molecular mechanism of ovarian development. Recently, with the advent of genome editing techniques, targeted gene knockout research can be carried out in many species. Zebrafish has emerged as an excellent model system to study the control of ovarian development. Dozens of genes related to ovarian development have been knocked out in zebrafish in recent years. Much new information and perspectives on the molecular mechanism of ovarian development have been obtained from these mutant zebrafish. Some findings have challenged conventional views. Several genes have been identified for the first time in vertebrates to control ovarian development. Focusing on ovarian development, the purpose of this review is to briefly summarize recent findings using these gene knockout zebrafish models, and compare these findings with mammalian models. These established mutants and rapid development of gene knockout techniques have prompted zebrafish as an ideal animal model for studying ovarian development.

A Hormone That Lost Its Receptor: Anti-Müllerian Hormone (AMH) in Zebrafish Gonad Development and Sex Determination

Fetal mammalian testes secrete Anti-Müllerian hormone (Amh), which inhibits female reproductive tract (Müllerian duct) development. Amh also derives from mature mammalian ovarian follicles, which marks oocyte reserve and characterizes polycystic ovarian syndrome. Zebrafish (Danio rerio) lacks Müllerian ducts and the Amh receptor gene amhr2 but, curiously, retains amh To discover the roles of Amh in the absence of Müllerian ducts and the ancestral receptor gene, we made amh null alleles in zebrafish. Results showed that normal amh prevents female-biased sex ratios. Adult male amh mutants had enormous testes, half of which contained immature oocytes, demonstrating that Amh regulates male germ cell accumulation and inhibits oocyte development or survival. Mutant males formed sperm ducts and some produced a few offspring. Young female mutants laid a few fertile eggs, so they also had functional sex ducts. Older amh mutants accumulated nonvitellogenic follicles in exceedingly large but sterile ovaries, showing that Amh helps control ovarian follicle maturation and proliferation. RNA-sequencing data partitioned juveniles at 21 days postfertilization (dpf) into two groups that each contained mutant and wild-type fish. Group21-1 upregulated ovary genes compared to Group21-2, which were likely developing as males. By 35 dpf, transcriptomes distinguished males from females and, within each sex, mutants from wild types. In adult mutants, ovaries greatly underexpressed granulosa and theca genes, and testes underexpressed Leydig cell genes. These results show that ancestral Amh functions included development of the gonadal soma in ovaries and testes and regulation of gamete proliferation and maturation. A major gap in our understanding is the identity of the gene encoding a zebrafish Amh receptor; we show here that the loss of amhr2 is associated with the breakpoint of a chromosome rearrangement shared among cyprinid fishes.

Genetic regulation of sex determination and maintenance in zebrafish (Danio rerio)

Over the last several decades zebrafish (Danio rerio) has become a major model organism for the study of vertebrate development and physiology. Given this, it may be surprising how little is known about the mechanism that zebrafish use to determine sex. While zebrafish are a gonochoristic species (having two sexes) that do not switch sex as adults, it was appreciated early on that sex ratios obtained from breeding lab domesticated lines were not typically a 1:1 ratio of male and female, suggesting that sex was not determined by a strict chromosomal mechanism. Here we will review the recent progress toward defining the genetic mechanism for sex determination in both wild and domesticated zebrafish.

Genetic regulation of sex determination and maintenance in zebrafish (Danio rerio)

Over the last several decades zebrafish (Danio rerio) has become a major model organism for the study of vertebrate development and physiology. Given this, it may be surprising how little is known about the mechanism that zebrafish use to determine sex. While zebrafish are a gonochoristic species (having two sexes) that do not switch sex as adults, it was appreciated early on that sex ratios obtained from breeding lab domesticated lines were not typically a 1:1 ratio of male and female, suggesting that sex was not determined by a strict chromosomal mechanism. Here we will review the recent progress toward defining the genetic mechanism for sex determination in both wild and domesticated zebrafish.

Loss-of-function of sox3 causes follicle development retardation and reduces fecundity in zebrafish

Folliculogenesis is essential for production of female gametes in vertebrates. However, the molecular mechanisms underlying follicle development, particularly apoptosis regulation in ovary, remain elusive. Here, we generated sox3 knockout zebrafish lines using CRISPR/Cas9. sox3 knockout led to follicle development retardation and a reduced fecundity in females. Comparative analysis of transcriptome between sox3^−/− and wild-type ovaries revealed that Sox3 was involved in pathways of ovarian steroidogenesis and apoptosis. Knockout of sox3 promoted follicle apoptosis and obvious apoptosis signals were detected in somatic cells of stages III and IV follicles of sox3^−/− ovaries. Moreover, Sox3 can bind to and activate the promoter of cyp19a1a. Up-regulation of Cyp19a1a expression promoted 17β-estradiol synthesis, which inhibited apoptosis in follicle development. Thus, Sox3 functions as a regulator of Cyp19a1a expression, via 17β-E2 linking apoptosis suppression, which is implicated in improving female fecundity.

Cloning, characterization and function analysis of DAX1 in Chinese loach (Paramisgurnus dabryanus)

The mechanisms of sex determination and differentiation have not been elucidated in most fish species. In this study, the full-length cDNAs of DAX1 was cloned and characterized in aquaculture fish Chinese loach (Paramisgurnus dabryanus), designated as Pd-DAX1. The cDNA sequence of Pd-DAX1 was 1261 bp, including 795 bp open reading frame (ORF) encoding 264 amino acids. Pd-DAX1 shares highly identical sequence with DAX1 homologues from different species. The expression profiles of Pd-DAX1 in different developmental stages and diverse adult tissues were analyzed by quantitative real-time RT-PCR and in situ hybridization (ISH). Pd-DAX1 was continuously expressed during embryogenesis, with the extensive distribution in the development of the central nervous system. Tissue distribution analysis revealed that Pd-DAX1 expressed widely in adult tissues, with the highest expression level found in testis, moderate level in ovary, showing a sex-dimorphic expression pattern. Pd-DAX1 mainly located in spermatogonia cells, spermatocytes, primary oocytes and previtellogenic oocyte cells, implying that Pd-DAX1 may involve in gametogenesis. These preliminary findings suggest that Pd-DAX1 gene is highly conserved during vertebrate evolution and involved in a wide range of developmental processes including embryogenesis, central nervous system development and gonad development.

Characteristics of Cyp11a during Gonad Differentiation of the Olive Flounder Paralichthys olivaceus

The P450 side-chain cleavage enzyme, P450scc (Cyp11a) catalyzes the first enzymatic step for the synthesis of all steroid hormones in fish. To study its roles in gonads of the olive flounder Paralichthys olivaceus, an important maricultured fish species, we isolated the cyp11a genomic DNA sequence of 1396 bp, which consists of 5 exons and 4 introns. Semi-quantitative reverse transcription polymerase chain reaction (RT-PCR) results indicated that the flounder cyp11a was exclusively expressed in gonad and head kidney tissues. Its expression level in the testis was higher than that in the ovary. According to the in situ hybridization patterns, cyp11a was mainly expressed in the Leydig cells of the testis, and the thecal cells of the ovary. Immunofluorescence analysis showed that Cyp11a was located in the cytoplasm of the cultured flounder testis cells. Further quantitative real-time PCR results presented the cyp11a differential expression patterns during gonad differentiation. Among different sampling points of the 17β-estradiol (E2, 5 ppm) treatment group, cyp11a expression levels were relatively high in the differentiating ovary (30 and 40 mm total length, TL), and then significantly decreased in the differentiated ovary (80, 100 and 120 mm TL, p < 0.05). The pregnenolone level also dropped in the differentiated ovary. In the high temperature treatment group (HT group, 28 ± 0.5 °C), the cyp11a expression level fluctuated remarkably in the differentiating testis (60 mm TL), and then decreased in the differentiated testis (80, 100 mm TL, p < 0.05). In the testosterone (T, 5 ppm) treatment group, the cyp11a was expressed highly in undifferentiated gonads and the differentiating testis, and then dropped in the differentiated testis. Moreover, the levels of cholesterol and pregnenolone of the differentiating testis in the HT and T groups increased. The expression level of cyp11a was significantly down-regulated after the cultured flounder testis cells were treated with 75 and 150 μM cyclic adenosine monophosphate (cAMP), respectively (p < 0.05), and significantly up-regulated after treatment with 300 μM cAMP (p < 0.05). Both nuclear receptors NR5a2 and NR0b1 could significantly up-regulate the cyp11a gene expression in a dosage dependent way in the testis cells detected by cell transfection analysis (p < 0.05). The above data provides evidence that cyp11a would be involved in the flounder gonad differentiation and development.

Zebrafish androgen receptor is required for spermatogenesis and maintenance of ovarian function

The androgen receptor (AR) is a nuclear receptor protein family member and inducible transcription factor that modulates androgen target gene expression. Studies using a mouse model confirmed the need for ar in reproductive development, particularly spermatogenesis. Here, we investigated the role of ar in zebrafish using CRISPR/Cas9 gene targeting technology. Targeted disruption of ar in zebrafish increases the number of female offspring and increases offspring weight. In addition, ar-null male zebrafish have female secondary sex characteristics. More importantly, targeted disruption of ar in zebrafish causes male infertility via defective spermatogenesis and female premature ovarian failure during growth. Mechanistic assays suggest that these effects are caused by fewer proliferated cells and more apoptotic cells in ar-null testes. Moreover, genes involved in reproductive development, estradiol induction and hormone synthesis were dys-regulated in testes and ovaries and the reproductive-endocrine axis was disordered. Our data thus suggest that the zebrafish ar is required for spermatogenesis and maintenance of ovarian function, which confirms evolutionarily conserved functions of ar in vertebrates, as well as indicates that ar-null zebrafish are a suitable model for studying pathologic mechanisms related to androgen disorders.

Distinct and Cooperative Roles of amh and dmrt1 in Self-Renewal and Differentiation of Male Germ Cells in Zebrafish

Spermatogenesis is a fundamental process in male reproductive biology and depends on precise balance between self-renewal and differentiation of male germ cells. However, the regulative factors for controlling the balance are poorly understood. In this study, we examined the roles of amh and dmrt1 in male germ cell development by generating their mutants with Crispr/Cas9 technology in zebrafish. Amh mutant zebrafish displayed a female-biased sex ratio, and both male and female amh mutants developed hypertrophic gonads due to uncontrolled proliferation and impaired differentiation of germ cells. A large number of proliferating spermatogonium-like cells were observed within testicular lobules of the amh-mutated testes, and they were demonstrated to be both Vasa- and PH3-positive. Moreover, the average number of Sycp3- and Vasa-positive cells in the amh mutants was significantly lower than in wild-type testes, suggesting a severely impaired differentiation of male germ cells. Conversely, all the dmrt1-mutated testes displayed severe testicular developmental defects and gradual loss of all Vasa-positive germ cells by inhibiting their self-renewal and inducing apoptosis. In addition, several germ cell and Sertoli cell marker genes were significantly downregulated, whereas a prominent increase of Insl3-positive Leydig cells was revealed by immunohistochemical analysis in the disorganized dmrt1-mutated testes. Our data suggest that amh might act as a guardian to control the balance between proliferation and differentiation of male germ cells, whereas dmrt1 might be required for the maintenance, self-renewal, and differentiation of male germ cells. Significantly, this study unravels novel functions of amh gene in fish.

Characterization and expression of StAR2a and StAR2b in the olive flounder Paralichthys olivaceus

Steroidogenic acute regulatory protein 2 (StAR2) is a key protein in transporting cholesterol from the outer mitochondria membrane to the inner mitochondria membrane for sex steroid synthesis. In this study, two StAR2 gene isoforms, StAR2a and StAR2b, were isolated from the olive flounder Paralichthys olivaceus gonads. Semi-quantitative RT-PCR results indicated that their expression levels were higher in testis than those in ovary. StAR2a was mainly expressed in the thecal cells and ooplasm of ovary, and Leydig cells and spermatid of testis according to the results of in situ hybridization. The quantitative real-time PCR results showed that the expressions of StAR2a and StAR2b were high in undifferentiation gonads and differentiating testis, and then decreased in differentiated testis in the high temperature (28°C) and exogenous testosterone treatment groups. While, in the exogenous 17β-estradiol treatment group, both genes' expression levels were high in differentiating ovary, and then significantly decreased in differentiated ovary (P<0.05). StAR2a and StAR2b expression levels were significantly down-regulated in the cultured testis cells treated with the 75 and 150μM cAMP, but significantly up-regulated in the cultured testis cells treated with the 300μM cAMP (P<0.05). Moreover, their expression levels were significantly up-regulated by transfecting the cultured testis cells with pcDNA3.1-NR5a2 and pcDNA3.1-NR0b1 (P<0.05). Above study showed that expression of StAR2 was regulated by cAMP and the transcription factors, NR5a2 and NR0b1, indicating that StAR2 may have functions in flounder gonadal differentiation and maintenance.

Dissection of Larval Zebrafish Gonadal Tissue

Although wild zebrafish possess a ZZ/ZW sex-determination system, domesticated zebrafish have lost the sex chromosome. They utilize a polygenic sex determination system, where several genes distributed throughout the genome collectively determine the sex identities of individual fish. Currently, the genes involved in regulating gonad development and how they work remain elusive. Normally, isolating gonadal tissue is the first step to examine the sex developmental processes. Here, we present a procedure to isolate gonadal tissue from 17 dpf (days post fertilization) and 25 dpf zebrafish larvae. The isolated gonadal tissue may be subsequently examined by morphology and gene expression profiling.

Gene editing nuclease and its application in tilapia

Gene editing nucleases including zinc-finger nucleases (ZFNs), transcription activator like effector nucleases (TALENs) and clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated (Cas) system (CRISPR/Cas9) provide powerful tools that improve our ability to understand the physiological processes and their underlying mechanisms. To date, these approaches have already been widely used to generate knockout and knockin models in a large number of species. Fishes comprise nearly half of extant vertebrate species and provide excellent models for studying many aspects of biology. In this review, we present an overview of recent advances in the use of gene editing nucleases for studies of fish species. We focus particularly on the use of TALENs and CRISPR/Cas9 genome editing for studying sex determination in tilapia.