Induction of kidney-related gene programs through co-option of SALL1 in mole ovotestes.

Changes in gene expression represent an important source of phenotypic innovation. Yet how such changes emerge and impact the evolution of traits remains elusive. Here, we explore the molecular mechanisms associated with the development of masculinizing ovotestes in female moles. By performing integrative analyses of epigenetic and transcriptional data in mole and mouse, we identified the co-option of SALL1 expression in mole ovotestes formation. Chromosome conformation capture analyses highlight a striking conservation of the 3D organization at the SALL1 locus, but an evolutionary divergence of enhancer activity. Interspecies reporter assays support the capability of mole-specific enhancers to activate transcription in urogenital tissues. Through overexpression experiments in transgenic mice, we further demonstrate the capability of SALL1 to induce kidney-related gene programs, which are a signature of mole ovotestes. Our results highlight the co-option of gene expression, through changes in enhancer activity, as a plausible mechanism for the evolution of traits.

The single-cell chromatin landscape in gonadal cell lineage specification.

Gonad development includes sex determination and divergent maturation of the testes and ovaries. Recent advances in measuring gene expression in single cells are providing new insights into this complex process. However, the underlying epigenetic regulatory mechanisms remain unclear. Here, we profiled chromatin accessibility in mouse gonadal cells of both sexes from embryonic day 11.5 to 14.5 using single-cell assay for transposase accessible chromatin by sequencing (scATAC-seq). Our results showed that individual cell types can be inferred by the chromatin landscape, and that cells can be temporally ordered along developmental trajectories. Integrative analysis of transcriptomic and chromatin-accessibility maps identified multiple putative regulatory elements proximal to key gonadal genes Nr5a1, Sox9 and Wt1. We also uncover cell type-specific regulatory factors underlying cell type specification. Overall, our results provide a better understanding of the epigenetic landscape associated with the progressive restriction of cell fates in the gonad.

Foxp2 deficiency impairs reproduction by modulating the hypothalamic-pituitary-gonadal axis in zebrafish†.

FOXP2 was initially characterized as a transcription factor linked to speech and language disorders. Single-cell RNA sequencing reveals that Foxp2 is enriched in the gonadotrope cluster of the pituitary gland and colocalized with the hormones LHB and FSHB in chickens and mice, implying that FOXP2 might be associated with reproduction in vertebrates. Herein, we investigated the roles of foxp2 in reproduction in a Foxp2-deficient zebrafish model. The results indicated that the loss of Foxp2 inhibits courtship behavior in adult male zebrafish. Notably, Foxp2 deficiency disrupts gonad development, leading to retardation of follicle development and a decrease in oocytes in females at the full-growth stage, among other phenotypes. The transcriptome analysis (RNA-seq) also revealed that differentially expressed genes clustered into the estrogen signaling and ovarian steroidogenesis-related signaling pathways. In addition, we found that Foxp2 deficiency could modulate the hypothalamic-pituitary-gonadal axis, especially the regulation of lhb and fshb expression, in zebrafish. In contrast, the injection of human chorionic gonadotropin, a specific LH agonist, partially rescues Foxp2-impaired reproduction in zebrafish, suggesting that Foxp2 plays an important role in the regulation of reproduction via the hypothalamic-pituitary-gonadal axis in zebrafish. Thus, our findings reveal a new role for Foxp2 in the regulation of reproduction in vertebrates.

Regulation of early spermatogenesis in the giant prawn Macrobrachium rosenbergii by a GCL homolog†.

The germ cell-less gene is crucial for gonad development in various organisms. Early interventions in its expression suggested a regulatory role at the mitotic stages of spermatogenesis, and its early knockout resulted in complete sterility in Drosophila. Genomic and transcriptomic data available for the catadromous giant prawn Macrobrachium rosenbergii enabled the identification of a germ cell-less homolog for this species, which we termed MroGCL (mRNA accession number OQ533056). An open reading frame containing 494 amino acids and a typical evolutionarily conserved BTB/POZ domain suggests possible protein-protein interaction functions in keeping with the Drosophila germ cell-less protein. Genomic mapping of MroGCL showed a full length of 120 896 bases. Analysis of the temporal expression of MroGCL showed constant expression in early prawn embryonic and larval stages, but a significant increase 10 days after metamorphosis when crucial sexual differentiation processes occur in prawns. In adult animals, high expression was detected in the gonads compared to the somatic tissues. RNAi-based knock-down experiments showed that both the silenced and control groups reached advanced spermatogenic stages, but that there was a significant decrease in the yield of spermatozoa in about half of the silenced animals. This finding supports our hypothesis that MroGCL is crucial for mitosis during early stage spermatogenesis. In conclusion, this study contributes to the understanding of crustacean gonad development and provides a stepping stone in the development of environmentally valuable sterile crustacean populations.

Expression characteristics of the cyp19a1b aromatase gene and its response to 17ß-estradiol treatment in largemouth bass (Micropterus salmoides).

To investigate the regulatory role of the cyp19a1b aromatase gene in the sexual differentiation of largemouth bass (Micropterus salmoides, LMB), we obtained the full-length cDNA sequence of cyp19a1b using rapid amplification of cDNA ends technique. Tissue expression characteristics and feedback with 17-ß-estradiol (E2) were determined using quantitative real-time PCR (qRT-PCR), while gonad development was assessed through histological section observations. The cDNA sequence of LMB cyp19a1b was found to be1950 base pairs (bp) in length, including a 5' untranslated region of 145 bp, a 3' untranslated region of 278 bp, and an open reading frame encoding a protein consisting of 1527 bp that encoded 508 amino acids. The qRT-PCR results indicated that cyp19a1b abundantly expressed in the brain, followed by the gonads, and its expression in the ovaries was significantly higher than that observed in the testes (P < 0.05). After feeding fish with E2 for 30 days, the expression of cyp19a1b in the pseudo-female gonads (XY-F) was significantly higher than that in males (XY-M) (P < 0.05), whereas expression did not differ significantly between XX-F and XY-F fish (P > 0.05). Although the expression of cyp19a1b in XY-F and XX-F fish was not significantly different after 60 days (P>0.05), both exhibited significantly higher levels than that of XY-M fish (P<0.05). Histological sections analysis showed the presence of oogonia in both XY-F and XX-F fish at 30 days, while spermatogonia were observed in XY-M fish. At 60 days, primary oocytes were abundantly observed in both XY-F and XX-F fish, while a few spermatogonia were visible in XY-M fish. At 90 days, the histological sections' results showed that a large number of oocytes were visible in XY-F and XX-F fish. Additionally, the gonads of XY-M fish contained numerous spermatocytes. These results suggest that cyp19a1b plays a pivotal role in the development of ovaries and nervous system development in LMB.

Lhx2 in germ cells suppresses endothelial cell migration in the developing ovary.

LIM-homeobox genes play multiple roles in developmental processes, but their roles in gonad development are not completely understood. Herein, we report that Lhx2, Ils2, Lmx1a, and Lmx1b are expressed in a sexually dimorphic manner in mouse, rat, and human gonads during sex determination. Amongst these, Lhx2 has female biased expression in the developing gonads of species with environmental and genetic modes of sex determination. Single-cell RNAseq analysis revealed that Lhx2 is exclusively expressed in the germ cells of the developing mouse ovaries. To elucidate the roles of Lhx2 in the germ cells, we analyzed the phenotypes of Lhx2 knockout XX gonads. While the gonads developed appropriately in Lhx2 knockout mice and the somatic cells were correctly specified in the developing ovaries, transcriptome analysis revealed enrichment of genes in the angiogenesis pathway. There was an elevated expression of several pro-angiogenic factors in the Lhx2 knockout ovaries. The elevated expression of pro-angiogenic factors was associated with an increase in numbers of endothelial cells in the Lhx2-/- ovaries at E13.5. Gonad recombination assays revealed that the increased numbers of endothelial cells in the XX gonads in absence of Lhx2 was due to ectopic migration of endothelial cells in a cell non-autonomous manner. We also found that, there was increased expression of several endothelial cell-enriched male-biased genes in Lhx2 knockout ovaries. Also, in absence of Lhx2, the migrated endothelial cells formed an angiogenic network similar to that of the wild type testis, although the coelomic blood vessel did not form. Together, our results suggest that Lhx2 in the germ cells is required to suppress vascularization in the developing ovary. These results suggest a need to explore the roles of germ cells in the control of vascularization in developing gonads. Preprint version of the article is available on BioRxiv at https://doi.org/10.1101/2022.03.07.483280.

Identification, Characterization and Functional Analysis of Fibroblast Growth Factors in Black Rockfish (Sebastes schlegelii).

Fibroblast growth factors (FGFs) are short polypeptides that play essential roles in various cellular biological processes, including cell migration, proliferation, and differentiation, as well as tissue regeneration, immune response, and organogenesis. However, studies focusing on the characterization and function of FGF genes in teleost fishes are still limited. In this study, we identified and characterized expression patterns of 24 FGF genes in various tissues of embryonic and adult specimens of the black rockfish (Sebates schlegelii). Nine FGF genes were found to play essential roles in myoblast differentiation, as well as muscle development and recovery in juvelines of S. schlegelii. Moreover, sex-biased expression pattern of multiple FGF genes was recorded in the species' gonads during its development. Among them, expression of the FGF1 gene was recorded in interstitial and sertoli cells of testes, promoting germ-cell proliferation and differentiation. In sum, the obtained results enabled systematic and functional characterization of FGF genes in S. schlegelii, laying a foundation for further studies on FGF genes in other large teleost fishes.

Quantitative proteomics analysis of chicken embryos reveals key proteins that affect right gonadal degeneration in females.

In birds, embryonic gonads of females develop in a way different from mammals, with the left one develops into a functional ovary, while the right one degenerates during embryogenesis. Here, we examined the proteomics profiles of the female and male left and right gonads at embryonic day 6.5 (E6.5) with the label free tandem mass spectrometry proteomics technique. The relative protein abundance of the left and right gonads of female and male embryos was determined to identify their differential proteins. Overall, a total of 7726 proteins were identified, of which 79 and 54 proteins were significantly different in female and male right gonads compared with female left gonads and male left gonads respectively. Bioinformatics analysis showed that the proteins DMRT1, ZFPM2, TSHZ3 were potentially associated with the degeneration of the right gonads in female embryos. The proteomics in this study provide clues for further elucidation of the pathways of sex determination, sex differentiation, and right gonadal degeneration in birds.

A time-course transcriptome analysis of gonads from yellow catfish (Pelteobagrus fulvidraco) reveals genes associated with gonad development.

BACKGROUND: The yellow catfish, Pelteobagrus fulvidraco, is a commercially important fish species. It is widely distributed in the fresh water areas of China, including rivers, lakes, and reservoirs. Like many other aquaculture fish species, people have observed significant size dimorphism between male and female yellow catfish and it shows a growth advantage in males. RESULTS: Here, at the first time, the time-course transcriptome was used to explore the various expression profiles of genes in different gonad developmental stages and genders. A total of 2696 different expression genes (DEGs) were identified from different stages. Based on these DEGs, 13 gonad development related genes were identified which showed time-specific or sex biased expression patterns. CONCLUSION: This study will provide the crucial information on the molecular mechanism of gonad development of female and male yellow catfish. Especially, during the different gonad development stages, these 13 gonad development related genes exhibit various expression patterns in female and male individual respectively. These results could inspire and facilitate us to understanding the various roles of these genes play in different gonad development stages and genders.

Differential DNA methylation in Pacific oyster reproductive tissue in response to ocean acidification.

BACKGROUND: There is a need to investigate mechanisms of phenotypic plasticity in marine invertebrates as negative effects of climate change, like ocean acidification, are experienced by coastal ecosystems. Environmentally-induced changes to the methylome may regulate gene expression, but methylome responses can be species- and tissue-specific. Tissue-specificity has implications for gonad tissue, as gonad-specific methylation patterns may be inherited by offspring. We used the Pacific oyster (Crassostrea gigas) - a model for understanding pH impacts on bivalve molecular physiology due to its genomic resources and importance in global aquaculture- to assess how low pH could impact the gonad methylome. Oysters were exposed to either low pH (7.31 ± 0.02) or ambient pH (7.82 ± 0.02) conditions for 7 weeks. Whole genome bisulfite sequencing was used to identify methylated regions in female oyster gonad samples. C- > T single nucleotide polymorphisms were identified and removed to ensure accurate methylation characterization. RESULTS: Analysis of gonad methylomes revealed a total of 1284 differentially methylated loci (DML) found primarily in genes, with several genes containing multiple DML. Gene ontologies for genes containing DML were involved in development and stress response, suggesting methylation may promote gonad growth homeostasis in low pH conditions. Additionally, several of these genes were associated with cytoskeletal structure regulation, metabolism, and protein ubiquitination - commonly-observed responses to ocean acidification. Comparison of these DML with other Crassostrea spp. exposed to ocean acidification demonstrates that similar pathways, but not identical genes, are impacted by methylation. CONCLUSIONS: Our work suggests DNA methylation may have a regulatory role in gonad and larval development, which would shape adult and offspring responses to low pH stress. Combined with existing molluscan methylome research, our work further supports the need for tissue- and species-specific studies to understand the potential regulatory role of DNA methylation.

Effect of water temperature and food availability on growth performance, sex ratio and gonadal development in juvenile convict cichlid (Amatitlanianigrofasciata).

Temperature and food availability are key drivers of growth and reproductive development in fishes, but information on how they interact is poorly understood. This study investigates the effects of water temperature and food availability on growth, sex ratio and gonadal development of the convict cichlid (Amatitlania nigrofasciata) which is an ornamental fish that may be a useful lab model. For this experiment, 180 juvenile convict cichlid (0.3 ± 0.02 g) were held at three different temperatures (26, 29 and 32 °C as T1, T2 and T3) and fed to satiation (S) or a restricted diet (R: half satiation) during a 56-day experimental period. Specific growth rate was significantly higher in T2S treatment than the other groups. The highest and lowest mean oocyte sizes were recorded in T1S and T3R groups, respectively. The sex ratio of fish held at 29 °C was male biased (female, 21.0%; male, 78.9%), but this was not seen at 26 °C (female, 47.6%; male, 52.4%) or 32 °C (female, 57.1%; male, 42.9%). In T1S and T1R treatments, oocytes developed more than the other treatments and in T2S group testicular development was more advanced than other groups. These results demonstrate the complex interplay of temperature and food availability on growth and reproductive development in the convict cichlid. Appropriate food availability significantly improves growth and reproductive processes, while restricted feeding decreases growth, survival rate and reproductive performance.

Ovotestes suggest cryptic genetic influence in a reptile model for temperature-dependent sex determination.

Sex determination and differentiation in reptiles is complex. Temperature-dependent sex determination (TSD), genetic sex determination (GSD) and the interaction of both environmental and genetic cues (sex reversal) can drive the development of sexual phenotypes. The jacky dragon (Amphibolurus muricatus) is an attractive model species for the study of gene-environment interactions because it displays a form of Type II TSD, where female-biased sex ratios are observed at extreme incubation temperatures and approximately 50 : 50 sex ratios occur at intermediate temperatures. This response to temperature has been proposed to occur due to underlying sex determining loci, the influence of which is overridden at extreme temperatures. Thus, sex reversal at extreme temperatures is predicted to produce the female-biased sex ratios observed in A. muricatus. The occurrence of ovotestes during development is a cellular marker of temperature sex reversal in a closely related species Pogona vitticeps. Here, we present the first developmental data for A. muricatus, and show that ovotestes occur at frequencies consistent with a mode of sex determination that is intermediate between GSD and TSD. This is the first evidence suggestive of underlying unidentified sex determining loci in a species that has long been used as a model for TSD.

Dnmt3aa but Not Dnmt3ab Is Required for Maintenance of Gametogenesis in Nile Tilapia (Oreochromis niloticus).

Dnmt3a, a de novo methyltransferase, is essential for mammalian germ line DNA methylation. Only one Dnmt3a is identified in mammals, and homozygous mutants of Dnmt3a are lethal, while two Dnmt3a paralogs, dnmt3aa and dnmt3ab, are identified in teleosts due to the third round of genome duplication, and homozygous mutants of dnmt3aa and dnmt3ab are viable in zebrafish. The expression patterns and roles of dnmt3aa and dnmt3ab in gonadal development remain poorly understood in teleosts. In this study, we elucidated the precise expression patterns of dnmt3aa and dnmt3ab in tilapia gonads. Dnmt3aa was highly expressed in oogonia, phase I and II oocytes and granulosa cells in ovaries and spermatogonia and spermatocytes in testes, while dnmt3ab was mainly expressed in ovarian granulosa cells and testicular spermatocytes. The mutation of dnmt3aa and dnmt3ab was achieved by CRISPR/Cas9 in tilapia. Lower gonadosomatic index (GSI), increased apoptosis of oocytes and spermatocytes and significantly reduced sperm quality were observed in dnmt3aa-/- mutants, while normal gonadal development was observed in dnmt3ab-/- mutants. Consistently, the expression of apoptotic genes was significantly increased in dnmt3aa-/- mutants. In addition, the 5-methylcytosine (5-mC) level in dnmt3aa-/- gonads was decreased significantly, compared with that of dnmt3ab-/- and wild type (WT) gonads. Taken together, our results suggest that dnmt3aa, not dnmt3ab, plays important roles in maintaining gametogenesis in teleosts.

Formation of the rete testis during mouse embryonic development.

BACKGROUND: The rete testis connects seminiferous tubules of the testis with efferent ducts having a mesonephric origin. The development of the rete testis is insufficiently studied, but there is evidence suggesting that it originates from gonadal cells. Here, the formation of the rete testis was investigated from E11.5 to E16.5 using immunofluorescent staining and 3D-modeling. RESULTS: The rete testis became visible by SOX9 and PAX8 staining starting from E12.5. It was located in the mesonephros but connected with testis cords formed by Sertoli cells expressing SOX9, AMH, DMRT1. Between E13.5 and E14.5, AMH+ network of testis cords at the mesonephric side began to disintegrate in a gradient-dependent manner along the anterior-posterior axis of the gonad and connections between testis cords gradually lost AMH becoming a part of the rete. Cells combining features of Sertoli and rete cells (PAX8+ AMH+ and DMRT1+ AMH- cells) were detected starting from E14.5, suggesting that some rete cells originated from Sertoli cells. The rete ovarii, a female counterpart of the rete testis, developed in a similar way as the rete testis until E13.5. CONCLUSIONS: A part of the rete testis originates from connections between testis cords. Evidence that Sertoli cells contribute to rete cells is provided.

Comparison of differential expression genes in ovaries and testes of Pearlscale angelfish Centropyge vrolikii based on RNA-Seq analysis.

Pearlscale angelfish Centropyge vrolikii is a kind of protogynous hermaphrodite fish with a natural sexual reversion. Under appropriate social conditions, a female fish can transform into a male fish spontaneously. It is an important prerequisite for artificial breeding to understand the process of its gonadal development and sexual reversion. Gonadal development is regulated by many sex-related genes. In this study, we used unreferenced RNA-Seq technology to sequence the ovary at the perinucleolus stage (OII), ovary at the yolk vesicle stage (OIV),IV and testis (T), respectively; screened the gonadal differential expression genes (DEGs); and analyzed the expression of these genes in different developmental stages of ovary and different sex gonads. The results showed that a total of 142,589 all-unigene samples were assembled, and gene annotation was performed by COG, GO, KEGG, KOG, Pfam, Swissprot, eggNOG, and NR functional database. Comparative analysis revealed that there were 1919 genes that were up-regulated and 1289 genes were down-regulated in comparison to OIV vs OII, while there were 3653 genes that were up-regulated and 2874 genes were down-regulated in comparison of OIV vs T, there were 3345 genes that were up-regulated and 2995 genes were down-regulated in comparison of the OII vs the T. At the same time, the results verified by RT-qPCR were consistent with the variation trend of transcriptome data. Among the results, amh, sox9b, dmrt1, dmrt2, cyp11a, cyp17a, and cyp19a were significantly expressed in the testes, while sox3, sox4, sox11, sox17, and hsd3b7 were significantly expressed in the ovaries. And, the expression of the amh, sox9b, dmrt2, and dmrt1 were low in the OII and OIV, while significantly increased during the ovotestis in the hermaphroditic period (OT), and finally reached the highest level in pure testis after sex reversal. The expression of sox3, sox4, hsd3b7, sox11, and sox17 was significantly reduced during the hermaphroditic period (OT). These results suggested that these genes may play an important role in the process of sex reversal. This study is helpful to further understand the molecular regulation mechanism of gonadal development and sexual reversion in Pearlscale angelfish and also provide important clues for future studies.

Developmental and Histopathological Changes in the Gonads of Female Japanese Medaka Following Exposure to Water from the Yellow River, China.

Gonad development and histopathological changes typically associated with endocrine disruption were evaluated in female Japanese medaka (Oryzias latipes) exposed to river water from four representative cross-sections in the Yellow River (YR), China. Fish were held in the river water treatments from fertilization. Advanced ovarian development was observed in fish exposed to river water from Qinhe cross-section at 20 days post-hatch (dph) and in fish exposed to river water from all four cross-sections at 60 dph. Histopathological changes including increased oocyte atresia, perifollicular cell hyperplasia/hypertrophy, changes in ovarian staging, interstitial fibrosis and interstitial proteinaceous fluid were observed in the gonads of fish at 60 dph after exposure to river water from some cross-sections. Cytoplasmic retraction and karyoplasmic clumping were observed in fish exposed to river water from all four cross-sections at 60 dph. The results indicate that development and reproductive function in Yellow River fish is impaired, placing fish populations at risk.

Gonadal supporting cells acquire sex-specific chromatin landscapes during mammalian sex determination.

Cis-regulatory elements are critical for the precise spatiotemporal regulation of genes during development. However, identifying functional regulatory sites that drive cell differentiation in vivo has been complicated by the high numbers of cells required for whole-genome epigenetic assays. Here, we identified putative regulatory elements during sex determination by performing ATAC-seq and ChIP-seq for H3K27ac in purified XX and XY gonadal supporting cells before and after sex determination in mice. We show that XX and XY supporting cells initiate sex determination with similar chromatin landscapes and acquire sex-specific regulatory elements as they commit to the male or female fate. To validate our approach, we identified a functional gonad-specific enhancer downstream of Bmp2, an ovary-promoting gene. This work increases our understanding of the complex regulatory network underlying mammalian sex determination and provides a powerful resource for identifying non-coding regulatory elements that could harbor mutations that lead to Disorders of Sexual Development.

Dynamics of the transcriptional landscape during human fetal testis and ovary development.

STUDY QUESTION: Which transcriptional program triggers sex differentiation in bipotential gonads and downstream cellular events governing fetal testis and ovary development in humans? SUMMARY ANSWER: The characterization of a dynamically regulated protein-coding and non-coding transcriptional landscape in developing human gonads of both sexes highlights a large number of potential key regulators that show an early sexually dimorphic expression pattern. WHAT IS KNOWN ALREADY: Gonadal sex differentiation is orchestrated by a sexually dimorphic gene expression program in XX and XY developing fetal gonads. A comprehensive characterization of its non-coding counterpart offers promising perspectives for deciphering the molecular events underpinning gonad development and for a complete understanding of the etiology of disorders of sex development in humans. STUDY DESIGN, SIZE, DURATION: To further investigate the protein-coding and non-coding transcriptional landscape during gonad differentiation, we used RNA-sequencing (RNA-seq) and characterized the RNA content of human fetal testis (N = 24) and ovaries (N = 24) from 6 to 17 postconceptional week (PCW), a key period in sex determination and gonad development. PARTICIPANTS/MATERIALS, SETTING, METHODS: First trimester fetuses (6-12 PCW) and second trimester fetuses (13-14 and 17 PCW) were obtained from legally induced normally progressing terminations of pregnancy. Total RNA was extracted from whole human fetal gonads and sequenced as paired-end 2 × 50 base reads. Resulting sequences were mapped to the human genome, allowing for the assembly and quantification of corresponding transcripts. MAIN RESULTS AND THE ROLE OF CHANCE: This RNA-seq analysis of human fetal testes and ovaries at seven key developmental stages led to the reconstruction of 22 080 transcripts differentially expressed during testicular and/or ovarian development. In addition to 8935 transcripts displaying sex-independent differential expression during gonad development, the comparison of testes and ovaries enabled the discrimination of 13 145 transcripts that show a sexually dimorphic expression profile. The latter include 1479 transcripts differentially expressed as early as 6 PCW, including 39 transcription factors, 40 long non-coding RNAs and 20 novel genes. Despite the use of stringent filtration criteria (expression cut-off of at least 1 fragment per kilobase of exon model per million reads mapped, fold change of at least 2 and false discovery rate adjusted P values of less than <1%), the possibility of assembly artifacts and of false-positive differentially expressed transcripts cannot be fully ruled out. LARGE-SCALE DATA: Raw data files (fastq) and a searchable table (.xlss) containing information on genomic features and expression data for all refined transcripts have been submitted to the NCBI GEO under accession number GSE116278. LIMITATIONS, REASONS FOR CAUTION: The intrinsic nature of this bulk analysis, i.e. the sequencing of transcripts from whole gonads, does not allow direct identification of the cellular origin(s) of the transcripts characterized. Potential cellular dilution effects (e.g. as a result of distinct proliferation rates in XX and XY gonads) may account for a few of the expression profiles identified as being sexually dimorphic. Finally, transcriptome alterations that would result from exposure to pre-abortive drugs cannot be completely excluded. Although we demonstrated the high quality of the sorted cell populations used for experimental validations using quantitative RT-PCR, it cannot be totally excluded that some germline expression may correspond to cell contamination by, for example, macrophages. WIDER IMPLICATIONS OF THE FINDINGS: For the first time, this study has led to the identification of 1000 protein-coding and non-coding candidate genes showing an early, sexually dimorphic, expression pattern that have not previously been associated with sex differentiation. Collectively, these results increase our understanding of gonad development in humans, and contribute significantly to the identification of new candidate genes involved in fetal gonad differentiation. The results also provide a unique resource that may improve our understanding of the fetal origin of testicular and ovarian dysgenesis syndromes, including cryptorchidism and testicular cancers. STUDY FUNDING/COMPETING INTEREST(S): This work was supported by the French National Institute of Health and Medical Research (Inserm), the University of Rennes 1, the French School of Public Health (EHESP), the Swiss National Science Foundation [SNF n° CRS115_171007 to B.J.], the French National Research Agency [ANR n° 16-CE14-0017-02 and n° 18-CE14-0038-02 to F.C.], the Medical Research Council [MR/L010011/1 to P.A.F.] and the European Community's Seventh Framework Programme (FP7/2007-2013) [under grant agreement no 212885 to P.A.F.] and from the European Union's Horizon 2020 Research and Innovation Programme [under grant agreement no 825100 to P.A.F. and S.M.G.]. There are no competing interests related to this study.

Kiss2 but not kiss1 is involved in the regulation of social stress on the gonad development in yellowtail clownfish, Amphiprion clarkii.

The yellowtail clownfish (Amphiprion clarkii) is a hermaphrodite fish, whose sex differentiation and gonad development are closely related to its social status. The kisspeptin/KissR system is regarded as a key factor mediating social stress on reproductive regulation. In order to understand the effects of social rank stress on the yellowtail clownfish gonadal differentiation, full-length cDNAs of two paralogous genes encoding kisspeptin (kiss1 and kiss2) and KissR (kissr2 and kissr3) were cloned and characterized. The results of real-time PCR showed that kiss1 was primarily expressed in the hypothalamus, and kiss2/kissr2 were abundantly expressed in the liver, while kissr3 was almost exclusively concentrated in the cerebellum and pituitary. Moreover, both Kiss1-10 and Kiss2-10 peptides could initiate downstream signaling pathways by interacting with cognate receptors expressed in eukaryotic cells. Among the three social status groups, the mRNA levels of kiss2 in the hypothalamus and pituitary as well as kissr2 in the pituitary were significantly higher in subordinate individuals (nonbreeders) than dominate individuals (females and males); while the mRNA levels of kissr3 in the hypothalamus and gonad were low in subordinate individuals. Furthermore, the plasma estradiol (E2) and testosterone (T) levels were higher in subordinate than dominate individuals. This study shows that kiss2 is involved in the regulation of social stress on the gonad development in the yellowtail clownfish, but not kiss1.

Effects of short-time exposure to atrazine on miRNA expression profiles in the gonad of common carp (Cyprinus carpio).

BACKGROUND: Atrazine is widely used in agriculture and is a known endocrine disrupting chemical. Atrazine can seep into the water body through surface, posing a potential threat to the aquatic ecological environment and human drinking water source. In vertebrate, studies have shown that it can affect reproduction and development seriously, but its molecular mechanism for aquatic animals is unknown. Aquaculture is very common in China, especially common carp, whose females grow faster than males. However, the effects of atrazine on the reproduction of carp, especially miRNA, have not been investigated. RESULTS: In this study, common carp (Cyprinus carpio) at two key developmental stages were exposed to atrazine in vitro. Sex ratio was observed to analyze the effect of atrazine on the sex. MiRNA expression profiles were analysed to identify miRNAs related to gonad development and to reveal the atrazine mechanisms interfering with gonad differentiation. The results showed that the sex ratio was biased towards females. Atrazine exposure caused significant alteration of multiple miRNAs. Predicted targets of differently-expressed miRNAs were involved in many reproductive biology signalling pathways. CONCLUSIONS: Our results indicate that atrazine promoted the expression of female-biased genes by decreasing miRNAs in primordial gonad. In addition, our results indicate that atrazine can up-regulate aromatase expression through miRNAs, which supports the hypothesis that atrazine has endocrine-disrupting activity by altering the gene expression profile of the Hypothalamus-Pituitary-Gonad axis through its corresponding miRNAs.