A trans-species missense SNP in Amhr2 is associated with sex determination in the tiger pufferfish, Takifugu rubripes (fugu)

Miniscule differences between sex chromosomes in the giant genome of a salamander

In the Mexican axolotl (Ambystoma mexicanum), sex is determined by a single Mendelian factor, yet its sex chromosomes do not exhibit morphological differentiation typical of many vertebrate taxa that possess a single sex-determining locus. As sex chromosomes are theorized to differentiate rapidly, species with undifferentiated sex chromosomes provide the opportunity to reconstruct early events in sex chromosome evolution. Whole genome sequencing of 48 salamanders, targeted chromosome sequencing and in situ hybridization were used to identify the homomorphic sex chromosome that carries an A. mexicanum sex-determining factor and sequences that are present only on the W chromosome. Altogether, these sequences cover ~300 kb of validated female-specific (W chromosome) sequence, representing ~1/100,000^th of the 32 Gb genome. Notably, a recent duplication of ATRX, a gene associated with mammalian sex-determining pathways, is one of few functional (non-repetitive) genes identified among these W-specific sequences. This duplicated gene (ATRW) was used to develop highly predictive markers for diagnosing sex and represents a strong candidate for a recently-acquired sex determining locus (or sexually antagonistic gene) in A. mexicanum.

The autosomal Gsdf gene plays a role in male gonad development in Chinese tongue sole (Cynoglossus semilaevis)

Gsdf is a key gene for testicular differentiation in teleost. However, little is known about the function of Gsdf in Chinese tongue sole (Cynoglossus semilaevis). In this study, we obtained the full-length Gsdf gene (CS-Gsdf), and functional characterization revealed its potential participation during germ cell differentiation in testes. CS-Gsdf transcription was predominantly detected in gonads, while the levels in testes were significantly higher than those in ovaries. During the different developmental stages in male gonads, the mRNA level was significantly upregulated at 86 dph, and a peak appeared at 120 dph; then, the level decreased at 1 and 2 yph. In situ hybridization revealed that CS-Gsdf mRNA was mainly localized in the Sertoli cells, spermatogonia, and spermatids in mature testes. After CS-Gsdf knockdown in the male testes cell line by RNA interference, a series of sex-related genes was influenced, including several sex differentiation genes, CS-Wnt4a, CS-Cyp19a1a and CS-Star. Based on these data, we speculated that CS-Gsdf may play a positive role in germ differentiation and proliferation via influencing genes related to sex differentiation.

The unusual rainbow trout sex determination gene hijacked the canonical vertebrate gonadal differentiation pathway

Evolutionary novelties require rewiring of transcriptional networks and/or the evolution of new gene functions. Sex determination (SD), one of the most plastic evolutionary processes, requires such novelties. Studies on the evolution of vertebrate SD revealed that new master SD genes are generally recruited from genes involved in the downstream SD regulatory genetic network. Only a single exception to this rule is currently known in vertebrates: the intriguing case of the salmonid master SD gene (sdY), which arose from duplication of an immune-related gene. This exception immediately posed the question of how a gene outside from the classical sex differentiation cascade could acquire its function as a male SD gene. Here we show that SdY became integrated in the classical vertebrate sex differentiation cascade by interacting with the Forkhead box domain of the female-determining transcription factor, Foxl2. In the presence of Foxl2, SdY is translocated to the nucleus where the SdY:Foxl2 complex prevents activation of the aromatase (cyp19a1a) promoter in cooperation with Nr5a1 (Sf1). Hence, by blocking a positive loop of regulation needed for the synthesis of estrogens in the early differentiating gonad, SdY disrupts a preset female differentiation pathway, consequently allowing testicular differentiation to proceed. These results also suggest that the evolution of unusual vertebrate master sex determination genes recruited from outside the classical pathway like sdY is strongly constrained by their ability to interact with the canonical gonadal differentiation pathway.

A novel PDZ domain-containing gene is essential for male sex differentiation and maintenance in yellow catfish (Pelteobagrus fulvidraco)

The sex-determining genes are found to be variable among different fish species. Yellow catfish (Pelteobagrus fulvidraco) is an important aquaculture fish species in China with XX/XY sex-determining type. Recently, YY super-male yellow catfish has been successfully produced by combining hormonal-induced sex reversal method with sex chromosome-linked markers. Here, we identified a novel PDZ domain-containing gene in yellow catfish designated as pfpdz1, in whose intron the sex-linked marker was located. The coding sequence of pfpdz1 in Y chromosome was identical to that in X chromosome except a missense SNP (A/T) that changes an amino acid (E8V) in the N-terminal region. Pfpdz1 displayed male-specific expression during sex differentiation. Overexpression of pfpdz1 using additive transgenesis induces XX ovary to differentiate into testis-like tissue, while the targeted inactivation of pfpdz1 in Y chromosome using CRISPR/Cas9-mediated mutagenesis triggers ovarian differentiation. Furthermore, we demonstrated that pfpdz1 initiates testicular differentiation through upregulating expression of amh, dmrt1 and sox9a1, as well as downregulating expression of cyp19a1, foxl2 and wnt4. Our data provide functional evidence that pfpdz1 is significant for male differentiation and maintenance in yellow catfish.

A comparative transcriptomic study on developmental gonads provides novel insights into sex change in the protandrous black porgy (Acanthopagrus schlegelii)

Protandrous black porgy (Acanthopagrus schlegelii) is a popular and valuable commercial marine fish in China and East Asian countries. Controlling and managing its breeding has been an imperative step towards obtaining a sustainable supply of this fish in aquaculture production systems. Therefore, study on the molecular mechanisms of sex change in black porgy has both scientific and commercial importance. Previously, we identified some candidate genes related to sex determination and differentiation from a high-quality genome assembly of the black porgy. In the present study, transcriptome sequencing of developmental gonads (including testis, ovotestis and ovary) of black porgy was performed to further investigate the sex-change mechanisms. Our results showed that the highly expressed male-related genes (dmrt1, piwi1, piwi2, sox9, sox30 and amh) at the male phase were significantly down-regulated to a substantial degree at the intersexual stage, and the female-related genes (jnk1, vasa, wnt4, figla and foxl2) were distinctly up-regulated when the fish grows into a female adult, suggesting the potential roles of these genes in sex change of the black porgy. These data also support a previous hypothesis that the femaleness will be switched on when the testis is entering the degenerated stage through the diminished dmrt1 expression. Our transcriptome data provide a very useful genomic resource for future studies on sex change and practical aquaculture in the black porgy.

Full genome survey and dynamics of gene expression in the greater amberjack Seriola dumerili

Background

Teleosts of the genus Seriola, commonly known as amberjacks, are of high commercial value in international markets due to their flesh quality and worldwide distribution. The Seriola species of interest to Mediterranean aquaculture is the greater amberjack (Seriola dumerili). This species holds great potential for the aquaculture industry, but in captivity, reproduction has proved to be challenging, and observed growth dysfunction hinders their domestication. Insights into molecular mechanisms may contribute to a better understanding of traits like growth and sex, but investigations to unravel the molecular background of amberjacks have begun only recently.

Findings

Illumina HiSeq sequencing generated a high-coverage greater amberjack genome sequence comprising 45 909 scaffolds. Comparative mapping to the Japanese yellowtail (Seriola quinqueriadiata) and to the model species medaka (Oryzias latipes) allowed the generation of in silico groups. Additional gonad transcriptome sequencing identified sex-biased transcripts, including known sex-determining and differentiation genes. Investigation of the muscle transcriptome of slow-growing individuals showed that transcripts involved in oxygen and gas transport were differentially expressed compared with fast/normal-growing individuals. On the other hand, transcripts involved in muscle functions were found to be enriched in fast/normal-growing individuals.

Conclusion

The present study provides the first insights into the molecular background of male and female amberjacks and of fast- and slow-growing fish. Therefore, valuable molecular resources have been generated in the form of a first draft genome and a reference transcriptome. Sex-biased genes, which may also have roles in sex determination or differentiation, and genes that may be responsible for slow growth are suggested.

Diverse and variable sex determination mechanisms in vertebrates

Sex is prevalent in nature and sex determination is one of the most fundamental biological processes, while the way of initiating female and male development exhibits remarkable diversity and variability across vertebrates. The knowledge on why and how sex determination mechanisms evolve unusual plasticity remains limited. Here, we summarize sex determination systems, master sex-determining genes and gene-regulatory networks among vertebrates. Recent research advancements on sex determination system transition are also introduced and discussed in some non-model animals with multiple sex determination mechanisms. This review will provide insights into the origin, transition and evolutionary adaption of different sex determination strategies in vertebrates, as well as clues for future perspectives in this field.

A draft genome of the striped catfish, Pangasianodon hypophthalmus, for comparative analysis of genes relevant to development and a resource for aquaculture improvement

Background

The striped catfish, Pangasianodon hypophthalmus, is a freshwater and benthopelagic fish common in the Mekong River delta. Catfish constitute a valuable source of dietary protein. Therefore, they are cultured worldwide, and P. hypophthalmus is a food staple in the Mekong area. However, genetic information about the culture stock, is unavailable for breeding improvement, although genetics of the channel catfish, Ictalurus punctatus, has been reported. To acquire genome sequence data as a useful resource for marker-assisted breeding, we decoded a draft genome of P. hypophthalmus and performed comparative analyses.

Results

Using the Illumina platform, we obtained both nuclear and mitochondrial DNA sequences. Molecular phylogeny using the mitochondrial genome confirmed that P. hypophthalmus is a member of the family Pangasiidae and is nested within a clade including the families Cranoglanididae and Ictaluridae. The nuclear genome was estimated at approximately 700 Mb, assembled into 568 scaffolds with an N50 of 14.29 Mbp, and was estimated to contain ~ 28,600 protein-coding genes, comparable to those of channel catfish and zebrafish. Interestingly, zebrafish produce gadusol, but genes for biosynthesis of this sunscreen compound have been lost from catfish genomes. The differences in gene contents between these two catfishes were found in genes for vitamin D-binding protein and cytosolic phospholipase A₂, which have lost only in channel catfish. The Hox cluster in catfish genomes comprised seven paralogous groups, similar to that of zebrafish, and comparative analysis clarified catfish lineage-specific losses of A5a, B10a, and A11a. Genes for insulin-like growth factor (IGF) signaling were conserved between the two catfish genomes. In addition to identification of MHC class I and sex determination-related gene loci, the hypothetical chromosomes by comparison with the channel catfish demonstrated the usefulness of the striped catfish genome as a marker resource.

Conclusions

We developed genomic resources for the striped catfish. Possible conservation of genes for development and marker candidates were confirmed by comparing the assembled genome to that of a model fish, Danio rerio, and to channel catfish. Since the catfish genomic constituent resembles that of zebrafish, it is likely that zebrafish data for gene functions is applicable to striped catfish as well.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5079-x) contains supplementary material, which is available to authorized users.

Unusual Diversity of Sex Chromosomes in African Cichlid Fishes

African cichlids display a remarkable assortment of jaw morphologies, pigmentation patterns, and mating behaviors. In addition to this previously documented diversity, recent studies have documented a rich diversity of sex chromosomes within these fishes. Here we review the known sex-determination network within vertebrates, and the extraordinary number of sex chromosomes systems segregating in African cichlids. We also propose a model for understanding the unusual number of sex chromosome systems within this clade.

Conservation of Sex-Linked Markers among Conspecific Populations of a Viviparous Skink, Niveoscincus ocellatus, Exhibiting Genetic and Temperature-Dependent Sex Determination

Sex determination systems are exceptionally diverse and have undergone multiple and independent evolutionary transitions among species, particularly reptiles. However, the mechanisms underlying these transitions have not been established. Here, we tested for differences in sex-linked markers in the only known reptile that is polymorphic for sex determination system, the spotted snow skink, Niveoscincus ocellatus, to quantify the genomic differences that have accompanied this transition. In a highland population, sex is determined genetically, whereas in a lowland population, offspring sex ratio is influenced by temperature. We found a similar number of sex-linked loci in each population, including shared loci, with genotypes consistent with male heterogamety (XY). However, population-specific linkage disequilibrium suggests greater differentiation of sex chromosomes in the highland population. Our results suggest that transitions between sex determination systems can be facilitated by subtle genetic differences.

Divergent Evolutionary Trajectories of Two Young, Homomorphic, and Closely Related Sex Chromosome Systems

There exists extraordinary variation among species in the degree and nature of sex chromosome divergence. However, much of our knowledge about sex chromosomes is based on comparisons between deeply diverged species with different ancestral sex chromosomes, making it difficult to establish how fast and why sex chromosomes acquire variable levels of divergence. To address this problem, we studied sex chromosome evolution in two species of African clawed frog (Xenopus), both of whom acquired novel systems for sex determination from a recent common ancestor, and both of whom have female (ZW/ZZ) heterogamy. Derived sex chromosomes of one species, X. laevis, have a small region of suppressed recombination that surrounds the sex determining locus, and have remained this way for millions of years. In the other species, X. borealis, a younger sex chromosome system exists on a different pair of chromosomes, but the region of suppressed recombination surrounding an unidentified sex determining gene is vast, spanning almost half of the sex chromosomes. Differences between these sex chromosome systems are also apparent in the extent of nucleotide divergence between the sex chromosomes carried by females. Our analyses also indicate that in autosomes of both of these species, recombination during oogenesis occurs more frequently and in different genomic locations than during spermatogenesis. These results demonstrate that new sex chromosomes can assume radically different evolutionary trajectories, with far-reaching genomic consequences. They also suggest that in some instances the origin of new triggers for sex determination may be coupled with rapid evolution sex chromosomes, including recombination suppression of large genomic regions.

Identification of genes involved in gonadal sex differentiation and the dimorphic expression pattern in Takifugu rubripes gonad at the early stage of sex differentiation

Quantifying the expression of mRNAs in the gonads at the critical stage of molecular sex differentiation stage might help to clarify the regulatory network during early sex differentiation and provide new information on the role of sex-related genes in gonadal function. In this study, transcriptomic analysis of sex-related genes expression profiles in fugu gonads at 60 and 90 days after hatching (dah) was conducted firstly, and a total of 112,504,991 clean reads, encompassing 28.35 Gb of sequences were retrieved. Twenty-three thousand eight hundred ten genes were found to be expressed in juvenile fugu gonads, and we mainly focused on the differentially expressed genes that have the potential to be involved in the gonadal sex differentiation. For 60-dah juveniles, we identified 1014 genes that were upregulated in the ovary and 1570 that were upregulated in the testis. For 90-dah juveniles, we identified 1287 genes that were upregulated in the ovary and 1500 that were upregulated in the testis. The dimorphic expression patterns of 15 genes in gonads at 30 and 40 dah were further investigate using qPCR. Cyp11b and star were expressed at higher levels in XY than in XX, while cyp11a1 and cyp19a1a were expressed at higher levels in XX than in XY at 30 dah. At 40 dah, the levels of gsdf, dmrt1, dmrt3, cyp11c1, star, and hsd3b expression were higher in XY, while the levels of foxl2, cyp19a1a, wnt9b, and foxD4 expression were higher in XX. Sox9, cyp11a1, cyp17a1, cyp17a2, and nr5a2 were expressed at similar levels in XX and XY at 40 dah. This is the first report of gonadal transcriptome of fugu at early sex differentiation stage, and our results provide an archive for further study on molecular mechanism underlying sex differentiation in this species.

Transitions in sex determination and sex chromosomes across vertebrate species

Despite the prevalence of sexual reproduction across eukaryotes, there is a remarkable diversity of sex-determination mechanisms. The underlying causes of this diversity remain unclear, and it is unknown whether there are convergent trends in the directionality of turnover in sex-determination mechanisms. We used the recently assembled Tree of Sex database to assess patterns in the evolution of sex-determination systems in the remarkably diverse vertebrate clades of teleost fish, squamate reptiles and amphibians. Contrary to theoretical predictions, we find no evidence that the evolution of separate sexes is irreversible, as transitions from separate sexes to hermaphroditism occur at higher rates than the reverse in fish. We also find that transitions from environmental sex determination to genetic sex determination occur at higher rates than the reverse in both squamates and fish, suggesting that genetic sex determination is more stable. However, our data are not consistent with the hypothesis that heteromorphic sex chromosomes are an "evolutionary trap." Rather, we find similar transition rates between homomorphic and heteromorphic sex chromosomes in both fish and amphibians, and to environmental sex determination from heteromorphic vs. homomorphic sex chromosome systems in fish. Finally, we find that transitions between male and female heterogamety occur at similar rates in amphibians and squamates, while transitions to male heterogamety occur at higher rates in fish. Together, these results provide the most comprehensive view to date of the evolution of vertebrate sex determination in a phylogenetic context, providing new insight into long-standing questions about the evolution of sexual reproduction.

Sex-specific dmrt1 and cyp19a1 methylation and alternative splicing in gonads of the protandrous hermaphrodite barramundi

Epigenetics is involved in sex differentiation of gonochoristic and hermaphroditic fish species, whereby two genes dmrt1 (pro-male) and cyp19a1 (pro-female) are known to play major roles. Barramundi, Lates calcarifer, is an important tropical aquaculture species that undergo natural and permanent male to female sex change, a process for which the exact underlying molecular mechanisms are still unknown. To elucidate whether DNA methylation is involved in sex control of barramundi, a next-generation bisulfite amplicon sequencing approach was used to target 146 CpG sites within proximal promoters and first exons of seven sex-related genes (dmrt1, cyp19a1, amh, foxl2, nr5a2, sox8 and sox9) of 24 testis and 18 ovaries of captive and wild adult barramundi. Moreover, comparative expression profiles of the key dmrt1 and cyp19a1 genes were further investigated using RT-qPCR and Sanger sequencing approaches, whereas expression levels of remaining targeted genes were based on available literature for the species. Results showed that cyp19a1 and amh were more methylated in males, whereas dmrt1 and nr5a2 were more methylated in females (P < 0.001), with no gender differences found for foxl2, sox8 or sox9 genes (P > 0.05). Sex-biased promoter DNA methylation was inversely related to gene expression only for dmrt1 and nr5a2, and directly related to amh expression, whereas no differences in cyp19a1 expression were found between testes and ovaries. Notably, unique sex-specific alternative splicing of dmrt1 and cyp19a1 were discovered, whereby males lacked the full-length aromatase coding cyp19a1 mRNA due to partial or total exon splicing, and females lacked the dmrt1 exon containing the DM-domain sequence. This study advances the current knowledge aiming to elucidate the genetic mechanisms within male and female gonads of this large protandrous hermaphrodite by providing the first evidence of epigenetics and alternative splicing simultaneously affecting key genes (cyp19a1 and dmrt1) central to sex differentiation pathways.

A Dense Linkage Map of Lake Victoria Cichlids Improved the Pundamilia Genome Assembly and Revealed a Major QTL for Sex-Determination

Genetic linkage maps are essential for comparative genomics, high quality genome sequence assembly and fine scale quantitative trait locus (QTL) mapping. In the present study we identified and genotyped markers via restriction-site associated DNA (RAD) sequencing and constructed a genetic linkage map based on 1,597 SNP markers of an interspecific F2 cross of two closely related Lake Victoria cichlids (Pundamilia pundamilia and P sp. 'red head'). The SNP markers were distributed on 22 linkage groups and the total map size was 1,594 cM with an average marker distance of 1.01 cM. This high-resolution genetic linkage map was used to anchor the scaffolds of the Pundamilia genome and estimate recombination rates along the genome. Via QTL mapping we identified a major QTL for sex in a ∼1.9 Mb region on Pun-LG10, which is homologous to Oreochromis niloticus LG 23 (Ore-LG23) and includes a well-known vertebrate sex-determination gene (amh).

Screening and characterisation of sex differentiation-related long non-coding RNAs in Chinese soft-shell turtle (Pelodiscus sinensis)

Long non-coding RNAs (lncRNAs) perform distinct functions in various biological processes in mammals, including sex differentiation. However, the roles of lncRNAs in other vertebrates, especially in the Chinese soft-shell turtle (Pelodiscus sinensis), remain to be clarified. In this study, we performed genome-wide analysis of the lncRNA expression profiles in gonad tissues and screened numerous sex-specific lncRNAs in the Chinese soft-shell turtle. Of the 363,310,650 clean reads obtained, 5,994 sequences were typed as lncRNAs, of which 4,463 were novel. A selection of sex-specific lncRNAs (♀ 932, ♂ 449) from female ovaries and male testis were shown to act on target genes in cis and in trans, and most were involved in gonad differentiation based on Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Furthermore, interactions among the differentially expressed lncRNA-mRNAs and protein coding genes were identified by construction of correlation networks. Overall, our systematic analysis of lncRNA expression profiles in gonad tissues revealed numerous sex-specific lncRNAs in P. sinensis. Thereby, these findings provide new insights into the function of lncRNAs in sex differentiation and highlight a group of candidate lncRNAs for future studies.

Dobzhansky-Muller incompatibilities, dominance drive, and sex-chromosome introgression at secondary contact zones: A simulation study

Dobzhansky-Muller (DM) incompatibilities involving sex chromosomes have been proposed to account for Haldane's rule (lowered fitness among hybrid offspring of the heterogametic sex) as well as Darwin's corollary (asymmetric fitness costs with respect to the direction of the cross). We performed simulation studies of a hybrid zone to investigate the effects of different types of DM incompatibilities on cline widths and positions of sex-linked markers. From our simulations, X-Y incompatibilities generate steep clines for both X-linked and Y-linked markers; random effects may produce strong noise in cline center positions when migration is high relative to fitness costs, but X- and Y-centers always coincide strictly. X-autosome and Y-autosome incompatibilities also generate steep clines, but systematic shifts in cline centers occur when migration is high relative to selection, as a result of a dominance drive linked to Darwin's corollary. Interestingly, sex-linked genes always show farther introgression than the associated autosomal genes. We discuss ways of disentangling the potentially confounding effects of sex biases in migration, we compare our results to those of a few documented contact zones, and we stress the need to study independent replicates of the same contact zone.

Discovery and functional characterization of microRNAs and their potential roles for gonadal development in spotted knifejaw, Oplegnathus punctatus

The spotted knifejaw (Oplegnathus punctatus) is a newly emerging economical fishery species in China. Studies focused on the regulation of gonadal development and gametogenesis of spotted knifejaw are still insufficient. As a key post-transcriptional regulator, miRNAs have been shown to play important roles in development and reproduction systems. In this study, small RNA deep sequencing in ovary and testis of spotted knifejaw were performed to screen miRNA expression patterns. After sequencing and bioinformatics analysis, a total of 247 conserved known miRNAs and 41 novel miRNAs were identified in spotted knifejaw gonads for the first time. In addition, 36 miRNAs were differentially expressed between testis and ovary. The putative target genes of differentially expressed (DE) miRNAs were significantly enriched in several pathways related to sexual differentiation and gonadal development, such as steroid hormone biosynthesis. Sequencing data was validated through qRT-PCR analysis of selected DE miRNAs. Dual-luciferase reporter analyses of filtered miRNA-target gene pairs confirmed that opu-miR-27b-3p targeted in piwi2 and mov10l1 3' UTRs and down-regulated their expressions in spotted knifejaw. The notion that mov10l1 and piwi2 enhance germ cells proliferation and regulate gonadal development and gametogenesis suggests that opu-miR-27b-3p may attenuated this process in the gonads of spotted knifejaw. These findings provided insights into regulatory roles of gonadal miRNAs and supplied fundamental resources for further studies on miRNA-mediated post-transcriptional regulation in reproductive system of spotted knifejaw.

Hypothalamus and pituitary transcriptome profiling of male and female Hong Kong grouper (Epinephelus akaara)

Hong Kong grouper (Epinephelus akaara) is an important commercially cultured marine fish in Asia, and a protogynous hermaphrodite with the "diandry" pattern. In order to explore the gene expression patterns of hypothalamus and pituitary between male and female Hong Kong grouper, we used RNA-seq technology to investigate transcriptomes of both tissues in immature and mature male and female adults. This produced 227,227,148 and 215,858,948 high quality reads from hypothalamus and pituitary, which were jointly assembled into 199,203 unigenes. Among them, 30,786 unigenes were mapped to known genes. Differential expression analysis revealed 275 unigenes that were differentially expressed between immature male and female adults and 561 between mature male and female adults. According to annotation and KEGG information, these differentially expressed genes (DEGs) were involved in development, metabolism, and regulation of transcription. One DEG, amino-terminal enhancer of split (AES), a member of the Groucho/transducin-like enhancer of split family of transcriptional regulators that played important roles in neurogenesis, segmentation, and sex determination, was significantly upregulated in male individuals in both immature and mature adult comparisons, indicating it may be involved in male reproductive function during development. Our report, for the first time, uses RNA-seq technology to investigate transcriptomes of both hypothalamus and pituitary in teleost fish, and provides a basis for further studies of molecular mechanism of sex determination and development in Hong Kong grouper.

Identification of Sex-associated SNPs of Greater Amberjack (Seriola dumerili)

The sex determination systems of fish are highly diverse compared with those of mammals. Thus, performing investigations using nonmodel fish species helps to understand the highly diverse sex determination systems of fish. Because greater amberjack (Seriola dumerili) is one of the most important edible fish globally and knowledge of its sex determination system is economically important in the field of aquaculture, we are interested in the mechanisms of sex determination of Seriola species. In this study, we identified sex-associated SNPs of greater amberjack using SNP information of 10 males and 10 females by an association test. We determined that the sex-associated SNPs were on chromosome 12 and mainly covered with two scaffolds (about 7.1 Mbp). Genotypes of sex-associated SNPs indicated that females are the heterogametic sex (ZZ/ZW). Furthermore, we compared the genomic structure of greater amberjack with those of Japanese amberjack (Seriola quinqueradiata), California yellowtail (Seriola dorsalis), and medaka (Oryzias latipes). Whole-genome alignments and synteny analysis indicated that the sex determination system of greater amberjack is markedly different from that of medaka and implied that the sex determination system is conserved in the Seriola species.

Transcriptome display during tilapia sex determination and differentiation as revealed by RNA-Seq analysis

Background

The factors determining sex in teleosts are diverse. Great efforts have been made to characterize the underlying genetic network in various species. However, only seven master sex-determining genes have been identified in teleosts. While the function of a few genes involved in sex determination and differentiation has been studied, we are far from fully understanding how genes interact to coordinate in this process.

Results

To enable systematic insights into fish sexual differentiation, we generated a dynamic co-expression network from tilapia gonadal transcriptomes at 5, 20, 30, 40, 90, and 180 dah (days after hatching), plus 45 and 90 dat (days after treatment) and linked gene expression profiles to both development and sexual differentiation. Transcriptomic profiles of female and male gonads at 5 and 20 dah exhibited high similarities except for a small number of genes that were involved in sex determination, while drastic changes were observed from 90 to 180 dah, with a group of differently expressed genes which were involved in gonadal differentiation and gametogenesis. Weighted gene correlation network analysis identified changes in the expression of Borealin, Gtsf1, tesk1, Zar1, Cdn15, and Rpl that were correlated with the expression of genes previously known to be involved in sex differentiation, such as Foxl2, Cyp19a1a, Gsdf, Dmrt1, and Amh.

Conclusions

Global gonadal gene expression kinetics during sex determination and differentiation have been extensively profiled in tilapia. These findings provide insights into the genetic framework underlying sex determination and sexual differentiation, and expand our current understanding of developmental pathways during teleost sex determination.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4756-0) contains supplementary material, which is available to authorized users.

Identification and Evolution of TGF-β Signaling Pathway Members in Twenty-Four Animal Species and Expression in Tilapia

Transforming growth factor β (TGF-β) signaling controls diverse cellular processes during embryogenesis as well as in mature tissues of multicellular animals. Here we carried out a comprehensive analysis of TGF-β pathway members in 24 representative animal species. The appearance of the TGF-β pathway was intrinsically linked to the emergence of metazoan. The total number of TGF-β ligands, receptors, and smads changed slightly in all invertebrates and jawless vertebrates analyzed. In contrast, expansion of the pathway members, especially ligands, was observed in jawed vertebrates most likely due to the second round of whole genome duplication (2R) and additional rounds in teleosts. Duplications of TGFB2, TGFBR2, ACVR1, SMAD4 and SMAD6, which were resulted from 2R, were first isolated. Type II receptors may be originated from the ACVR2-like ancestor. Interestingly, AMHR2 was not identified in Chimaeriformes and Cypriniformes even though they had the ligand AMH. Based on transcriptome data, TGF-β ligands exhibited a tissue-specific expression especially in the heart and gonads. However, most receptors and smads were expressed in multiple tissues indicating they were shared by different ligands. Spatial and temporal expression profiles of 8 genes in gonads of different developmental stages provided a fundamental clue for understanding their important roles in sex determination and reproduction. Taken together, our findings provided a global insight into the phylogeny and expression patterns of the TGF-β pathway genes, and hence contribute to the greater understanding of their biological roles in the organism especially in teleosts.

Sex Chromosome Evolution and Genomic Divergence in the Fish Hoplias malabaricus (Characiformes, Erythrinidae)

The Erythrinidae family (Teleostei: Characiformes) is a small Neotropical fish group with a wide distribution throughout South America, where Hoplias malabaricus corresponds to the most widespread and cytogenetically studied taxon. This species possesses significant genetic variation, as well as huge karyotype diversity among populations, as reflected by its seven major karyotype forms (i.e., karyomorphs A-G) identified up to now. Although morphological differences in their bodies are not outstanding, H. malabaricus karyomorphs are easily identified by differences in 2n, morphology and size of chromosomes, as well as by distinct evolutionary steps of sex chromosomes development. Here, we performed comparative genomic hybridization (CGH) to analyse both the intra- and inter-genomic status in terms of repetitive DNA divergence among all but one (E) H. malabaricus karyomorphs. Our results indicated that they have close relationships, but with evolutionary divergences among their genomes, yielding a range of non-overlapping karyomorph-specific signals. Besides, male-specific regions were uncovered on the sex chromosomes, confirming their differential evolutionary trajectories. In conclusion, the hypothesis that H. malabaricus karyomorphs are result of speciation events was strengthened.

The transcriptomic signature of different sexes in two protogynous hermaphrodites: Insights into the molecular network underlying sex phenotype in fish

Sex differentiation is a puzzling problem in fish due to the variety of reproductive systems and the flexibility of their sex determination mechanisms. The Sparidae, a teleost family, reflects this remarkable diversity of sexual mechanisms found in fish. Our aim was to capture the transcriptomic signature of different sexes in two protogynous hermaphrodite sparids, the common pandora Pagellus erythrinus and the red porgy Pagrus pagrus in order to shed light on the molecular network contributing to either the female or the male phenotype in these organisms. Through RNA sequencing, we investigated sex-specific differences in gene expression in both species' brains and gonads. The analysis revealed common male and female specific genes/pathways between these protogynous fish. Whereas limited sex differences found in the brain indicate a sexually plastic tissue, in contrast, the great amount of sex-biased genes observed in gonads reflects the functional divergence of the transformed tissue to either its male or female character. Α common "crew" of well-known molecular players is acting to preserve either sex identity of the gonad in these fish. Lastly, this study lays the ground for a deeper understanding of the complex process of sex differentiation in two species with an evolutionary significant reproductive system.

Gonadal development and expression of sex-specific genes during sex differentiation in the Japanese eel

The process of gonadal development and mechanism involved in sex differentiation in eels are still unclear. The objectives were to investigate the gonadal development and expression pattern of sex-related genes during sex differentiation in the Japanese eel, Anguilla japonica. For control group, the elvers of 8-10cm were reared for 8months; and for feminization, estradiol-17β (E2) was orally administered to the elvers of 8-10cm for 6months. Only males were found in the control group, suggesting a possible role of environmental factors in eel sex determination. In contrast, all differentiated eels in E2-treated group were female. Gonad histology revealed that control male eels seem to differentiate through an intersexual stage, while female eels (E2-treated) would differentiate directly from an undifferentiated gonad. Tissue distribution and sex-related genes expression during gonadal development were analyzed by qPCR. The vasa, figla and sox3 transcripts in gonads were significantly increased during sex differentiation. High vasa expression occurred in males; figla and sox3 were related to ovarian differentiation. The transcripts of dmrt1 and sox9a were significantly increased in males during testicular differentiation and development. The cyp19a1 transcripts were significantly increased in differentiating and differentiated gonads, but did not show a differential expression between the control and E2-treated eels. This suggests that cyp19a1 is involved both in testicular differentiation and development in control males, and in the early stage of ovarian differentiation in E2-treated eels. Importantly, these results also reveal that cyp19a1 is not a direct target for E2 during gonad differentiation in the eel.

Sex-specific markers developed by next-generation sequencing confirmed an XX/XY sex determination system in bighead carp (Hypophthalmichehys nobilis) and silver carp (Hypophthalmichthys molitrix)

Sex-specific markers are powerful tools for identifying sex-determination system in various animals. Bighead carp (Hypophthalmichehys nobilis) and silver carp (Hypophthalmichthys molitrix) are two of the most important edible fish in Asia, which have a long juvenility period that can lasts for 4-5 years. In this study, we found one sex-specific marker by next-generation sequencing together with bioinformatics analysis in bighead carp. The male-specific markers were used to perform molecular sexing in the progenies of artificial gynogenetic diploids and found all progenies (n = 160) were females. Meanwhile, around 1 : 1 sex ratio was observed in a total of 579 juvenile offspring from three other families. To further extend the male-specific region, we performed genome walking and got a male-specific sequence of 8,661 bp. Five pairs of primers were designed and could be used to efficiently distinguish males from females in bighead carp and silver carp. The development of these male-specific markers and results of their molecular sexing in different populations provide strong evidence for a sex determination system of female homogametry or male heterogametry (XX/XY) in bighead carp and silver carp. To the best of our knowledge, this is the first report of effective sex-specific markers in these two large carp species.

Identification of the sex-determining locus in grass puffer (Takifugu niphobles) provides evidence for sex-chromosome turnover in a subset of Takifugu species

There is increasing evidence for frequent turnover in sex chromosomes in vertebrates. Yet experimental systems suitable for tracing the detailed process of turnover are rare. In theory, homologous turnover is possible if the new sex-determining locus is established on the existing sex-chromosome. However, there is no empirical evidence for such an event. The genus Takifugu includes fugu (Takifugu rubripes) and its two closely-related species whose sex is most likely determined by a SNP at the Amhr2 locus. In these species, males are heterozygous, with G and C alleles at the SNP site, while females are homozygous for the C allele. To determine if a shift in the sex-determining locus occurred in another member of this genus, we used genetic mapping to characterize the sex-chromosome systems of Takifugu niphobles. We found that the G allele of Amhr2 is absent in T. niphobles. Nevertheless, our initial mapping suggests a linkage between the phenotypic sex and the chromosome 19, which harbors the Amhr2 locus. Subsequent high-resolution analysis using a sex-reversed fish demonstrated that the sex-determining locus maps to the proximal end of chromosome 19, far from the Amhr2 locus. Thus, it is likely that homologous turnover involving these species has occurred. The data also showed that there is a male-specific reduction of recombination around the sex-determining locus. Nevertheless, no evidence for sex-chromosome differentiation was detected: the reduced recombination depended on phenotypic sex rather than genotypic sex; no X- or Y-specific maker was obtained; the YY individual was viable. Furthermore, fine-scale mapping narrowed down the new sex-determining locus to the interval corresponding to approximately 300-kb of sequence in the fugu genome. Thus, T. niphobles is determined to have a young and small sex-determining region that is suitable for studying an early phase of sex-chromosome evolution and the mechanisms underlying turnover of sex chromosome.

Production of Tiger Puffer Takifugu rubripes Offspring from Triploid Grass Puffer Takifugu niphobles Parents

The tiger puffer Takifugu rubripes is one of the most popular aquacultural fish; however, there are two major obstacles to selective breeding. First, they have a long generation time of 2 or 3 years until maturation. Second, the parental tiger puffer has a body size (2-5 kg) much larger than average market size (0.6-1.0 kg). The grass puffer Takifugu niphobles is closely related to the tiger puffer and matures in half the time. Furthermore, grass puffer can be reared in small areas since their maturation weight is about 1/150 that of mature tiger puffer. Therefore, to overcome the obstacles of maturation size and generation time of tiger puffer, we generated surrogate grass puffer that can produce tiger puffer gametes through germ cell transplantation. Approximately 5000 tiger puffer testicular cells were transplanted into the peritoneal cavity of triploid grass puffer larvae at 1 day post hatching. When the recipient fish matured, both males and females produced donor-derived gametes. Through their insemination, we successfully produced donor-derived tiger puffer offspring presenting the same body surface dot pattern, number of dorsal fin rays, and DNA fingerprint as those of the donor tiger puffer, suggesting that the recipient grass puffer produced functional eggs and sperm derived from the donor tiger puffer. Although fine tunings are still needed to improve efficiencies, surrogate grass puffer are expected to accelerate the breeding process of tiger puffer because of their short generation time and small body size.

Brain Sexual Differentiation and Requirement of SRY: Why or Why Not?

Brain sexual differentiation is orchestrated by precise coordination of sex steroid hormones. In some species, programming of select male brain regions is dependent upon aromatization of testosterone to estrogen. In mammals, these hormones surge during the organizational and activational periods that occur during perinatal development and adulthood, respectively. In various fish and reptiles, incubation temperature during a critical embryonic period results in male or female sexual differentiation, but this can be overridden in males by early exposure to estrogenic chemicals. Testes development in mammals requires a Y chromosome and testis determining gene SRY (in humans)/Sry (all other therian mammals), although there are notable exceptions. Two species of spiny rats: Amami spiny rat (Tokudaia osimensis) and Tokunoshima spiny rat (Tokudaia tokunoshimensis) and two species of mole voles (Ellobius lutescens and Ellobius tancrei), lack a Y chromosome/Sry and possess an XO chromosome system in both sexes. Such rodent species, prototherians (monotremes, who also lack Sry), and fish and reptile species that demonstrate temperature sex determination (TSD) seemingly call into question the requirement of Sry for brain sexual differentiation. This review will consider brain regions expressing SRY/Sry in humans and rodents, respectively, and potential roles of SRY/Sry in the brain will be discussed. The evidence from various taxa disputing the requirement of Sry for brain sexual differentiation in mammals (therians and prototherians) and certain fish and reptilian species will be examined. A comparative approach to address this question may elucidate other genes, pathways, and epigenetic modifications stimulating brain sexual differentiation in vertebrate species, including humans.

Gonadal soma controls ovarian follicle proliferation through Gsdf in zebrafish

BACKGROUND

Aberrant signaling between germ cells and somatic cells can lead to reproductive disease and depends on diffusible signals, including transforming growth factor-beta (TGFB) -family proteins. The TGFB-family protein Gsdf (gonadal soma derived factor) controls sex determination in some fish and is a candidate for mediating germ cell/soma signaling.

RESULTS

Zebrafish expressed gsdf in somatic cells of bipotential gonads and expression continued in ovarian granulosa cells and testicular Sertoli cells. Homozygous gsdf knockout mutants delayed leaving the bipotential gonad state, but then became a male or a female. Mutant females ovulated a few oocytes, then became sterile, accumulating immature follicles. Female mutants stored excess lipid and down-regulated aromatase, gata4, insulin receptor, estrogen receptor, and genes for lipid metabolism, vitellogenin, and steroid biosynthesis. Mutant females contained less estrogen and more androgen than wild-types. Mutant males were fertile. Genomic analysis suggests that Gsdf, Bmp15, and Gdf9, originated as paralogs in vertebrate genome duplication events.

CONCLUSIONS

In zebrafish, gsdf regulates ovarian follicle maturation and expression of genes for steroid biosynthesis, obesity, diabetes, and female fertility, leading to ovarian and extra-ovarian phenotypes that mimic human polycystic ovarian syndrome (PCOS), suggesting a role for a related TGFB signaling molecule in the etiology of PCOS. Developmental Dynamics 246:925-945, 2017. © 2017 Wiley Periodicals, Inc.

Expression profile of amh/Amh during bi-directional sex change in the protogynous orange-spotted grouper Epinephelus coioides

Gonadal differentiation is tightly regulated by the initial sex determining gene and the downstream sex-related genes in vertebrates. However, sex change in fish can alter the sexual fate from one sex to the other. Chemical-induced maleness in the protogynous orange-spotted grouper is transient, and a reversible sex change occurs after the chemical treatment is withdrawn. We used these characteristics to study Amh signaling during bi-directional sex change in the grouper. We successfully induced the female-to-male sex change by chemical (aromatase inhibitor, AI, or methyltestosterone, MT) treatment. A dormant gonad (a low proliferation rate of early germ cells and no characteristics of both sexes) was found during the transient phase of reversible male-to-female sex change after the withdrawal of chemical administration. Our results showed that amh (anti-mullerian hormone) and its receptor amhr2 (anti-mullerian hormone receptor type 2) were significantly increased in the gonads during the process of female-to-male sex change. Amh is expressed in the Sertoli cells surrounding the type A spermatogonia in the female-to-male grouper. Male-related gene (dmrt1 and sox9) expression was immediately decreased in MT-terminated males during the reversible male-to-female sex change. However, Amh expression was found in the surrounding cells of type A spermatogonia-like cells during the transient phase of reversible male-to-female sex change. This phenomenon is correlated with the dormancy of type A spermatogonia-like cells. Thus, Amh signaling is suggested to play roles in regulating male differentiation during the female-to-male sex change and in inhibiting type-A spermatogonia-like cell proliferation/differentiation during the reversible male-to-female sex change. We suggest that Amh signaling might play dual roles during bi-directional sex change in grouper.

Transcriptome analysis of three critical periods of ovarian development in Yellow River carp (Cyprinus carpio)

Ovary development is a complex process involving numerous genes; the molecular mechanism underlying the ovary development of carp is still unknown. Here we used Illumina HiSeq™ 2500 to explore the transcriptome of undifferentiated gland (PG), juvenile ovary (OJ) and adult ovary (OA) of Yellow River carp (Cyprinus carpio). A total of 58,749 unigenes were obtained, comprising 45,707 known genes and 13,042 new genes. We identified differentially-expressed genes (DEGs) during development and characterized the functional properties of DEGs by comparison with the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes databases. qRT-PCR was used to analyze the expression of 22 DEGs and the results corresponded with those of RNA-Seq. Among DEGs between PG and OJ, some upstream regulators of gonad development were upregulated in PG, such as cyp19a and sox9, while some oocyte-specific genes were upregulated in OJ, such as nobox, bmp15 and zp2. Among DEGs between OJ and OA, many oocyte physiological function-related genes were upregulated in OA, such as fem-1 and foxl2. GO analysis showed a higher number of DEGs from PG-OJ analysis were assigned to reproduction terms. Furthermore, our investigation has also revealed DEGs identified from PG-OJ analysis were enriched in several important functional pathways, such as Fanconi anemia and the notch signal pathway. These data suggested a dynamic shift in gene expression during ovary development, and DEGs between PG and OJ provided crucial candidate gene data for the study of ovarian differentiation. Additionally, a total of 1,776,769 single nucleotide polymorphisms and 157,279 INDEs were revealed from transcriptome data. This result will contribute to knowledge of ovary differentiation of Yellow River carp.

Integrated analysis of mRNA-seq and miRNA-seq reveals the potential roles of sex-biased miRNA-mRNA pairs in gonad tissue of dark sleeper (Odontobutis potamophila)

Background

The dark sleeper (Odontobutis potamophila) is an important commercial fish species which shows a sexually dimorphic growth pattern. However, the lack of sex transcriptomic data is hindering further research and genetically selective breeding of the dark sleeper. In this study, integrated analysis of mRNA and miRNA was performed on gonad tissue to elucidate the molecular mechanisms of sex determination and differentiation in the dark sleeper.

Results

A total of 143 differentially expressed miRNAs and 16,540 differentially expressed genes were identified. Of these, 8103 mRNAs and 75 miRNAs were upregulated in testes, and 8437 mRNAs and 68 miRNAs were upregulated in ovaries. Integrated analysis of miRNA and mRNA expression profiles predicted more than 50,000 miRNA-mRNA interaction sites, and among them 27,583 negative miRNA-mRNA interactions. A number of sex related genes were targeted by sex-biased miRNAs. The relationship between 15 sex-biased genes and 15 sex-biased miRNAs verified by using qRT-PCR were described. Additionally, a number of SNPs were revealed through the transcriptome data.

Conclusions

The overall results of this study facilitate our understanding of the molecular mechanism underlying sex determination and differentiation and provide valuable genomic information for selective breeding of the dark sleeper.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3995-9) contains supplementary material, which is available to authorized users.

Transcriptomic analysis of the differentiating ovary of the protogynous ricefield eel Monopterus albus

BACKGROUND

The ricefield eel is a protogynous hermaphroditic Synbranchiform species that changes sex naturally from female to male, which offers an interesting model for studying gonadal (particularly ovarian) differentiation in vertebrates. In the present study, transcriptome sequencing of the gonad of ricefield eel larvae was performed to explore the molecular mechanisms underlying the ovarian differentiation and development.

RESULTS

A total of 301,267,988 clean reads were generated from cDNA libraries of gonadal tissues of ricefield eel larvae at 6, 9, 12, and 20 days post hatching (dph), which contained undifferentiated gonads, differentiating ovaries, ovaries with oogonia, and ovaries with meiotic oocytes, respectively. De-novo assembly of all the clean reads generated a total of 265,896 unigenes with a mean size of 720 bp and a N50 of 1107 bp. RT-qPCR analysis of the developmental expression of 13 gonadal development-related functional genes indicated that RNA-seq data are reliable. Transcriptome data suggest that high expression of female development-related genes and low expression of male development-related genes in the early gonads of ricefield eel larvae participate in the cascade of sex differentiation leading to the final female phenotype. The contrasting expression patterns of genes involved in retinoid acid (RA) synthesis and degradation might result in peak production of RA at 12 dph in the gonad of ricefield eel larvae, and induce molecular events responsible for the initiation of meiosis before the meiotic signs could be observed at 20 dph. In addition, only stra6 but not stra8 could be identified in gonadal transcriptome data of ricefield eel larvae, and the expression pattern of stra6 paralleled those of genes involved in RA synthesis, suggesting that stra6 may be a downstream target of RA and play a role in RA metabolism and/or meiotic initiation in the gonad of ricefield eel larvae.

CONCLUSIONS

The present study depicted the first large-scale RNA sequencing of the gonad of ricefield eel larvae, and identified many important functional genes, GO terms and KEGG pathways involved in gonadal development and germ cell meiosis. Results of the present study will facilitate future study on the ovarian differentiation of ricefield eels and other teleosts as well.

Dmrt1 is required for primary male sexual differentiation in Chinese soft-shelled turtle Pelodiscus sinensis

In vertebrates, the primary sex-determining signals that initiate sexual development are remarkably diverse, ranging from complete genetic to environmental cues. However, no sex determination-related genes have been functionally identified in reptiles. Here, we characterized a conserved DM domain gene, Dmrt1, in Chinese soft-shelled turtle Pelodiscus sinensis (P. sinensis), which exhibits ZZ/ZW sex chromosomes. Dmrt1 exhibited early male-specific embryonic expression, preceding the onset of gonadal sex differentiation. The expression of Dmrt1 was induced in ZW embryonic gonads that were masculinized by aromatase inhibitor treatment. Dmrt1 knockdown in ZZ embryos by RNA interference resulted in male to female sex reversal, characterized by obvious feminization of gonads, significant down-regulation of testicular markers Amh and Sox9, and remarkable up-regulation of ovarian regulators, Cyp19a1 and Foxl2. Conversely, ectopic expression of Dmrt1 led to largely masculinized genetic females, production of Amh and Sox9, and a decline in Cyp19a1 and Foxl2. These findings demonstrate that Dmrt1 is both necessary and sufficient to initiate testicular development, thereby acting as an upstream regulator of the male pathway in P. sinensis.

Identification of Sex-determining Loci in Pacific White Shrimp Litopeneaus vannamei Using Linkage and Association Analysis

The Pacific white shrimp Litopenaeus vannamei is a predominant aquaculture shrimp species in the world. Like other animals, the L. vannamei exhibited sexual dimorphism in growth trait. Mapping of the sex-determining locus will be very helpful to clarify the sex determination system and further benefit the shrimp aquaculture industry towards the production of mono-sex stocks. Based on the data used for high-density linkage map construction, linkage-mapping analysis was conducted. The sex determination region was mapped in linkage group (LG) 18. A large region from 0 to 21.205 cM in LG18 showed significant association with sex. However, none of the markers in this region showed complete association with sex in the other populations. So an association analysis was designed using the female parent, pool of female progenies, male parent, and pool of male progenies. Markers were de novo developed and those showing significant differences between female and male pools were identified. Among them, three sex-associated markers including one fully associated marker were identified. Integration of linkage and association analysis showed that the sex determination region was fine-mapped in a small region along LG18. The identified sex-associated marker can be used for the sex detection of this species at genetic level. The fine-mapped sex-determining region will contribute to the mapping of sex-determining gene and help to clarify sex determination system for L. vannamei.

Sexually dimorphic expression in developing and adult gonads shows an important role of gonadal soma-derived factor during sex differentiation in olive flounder (Paralichthys olivaceus)

Gonadal soma-derived factor (gsdf) is a new member of transforming growth factor beta (TGF-β) superfamily. As a teleost- and gonad-specific growth factor, gsdf has been indicated to play an important role in early germ cell development. However, little is known about its significance in germ cell development of olive flounder (Paralichthys olivaceus). In the present study, a 1338 bp gsdf gene was isolated from P. olivaceus for the first time. Bioinformatic analysis revealed that the genomic structure and synteny relationship of gsdf in teleosts were conserved. Quantitative real-time PCR (qRT-PCR) showed that gsdf expressed before sex gonadal differentiation, and the expression level increased rapidly after initiation of sex differentiation in males. In adult individuals, the expression of gsdf was higher in testis than that in ovary (P<0.01). In situ hybridization (ISH) indicated that gsdf mRNA was detected in the somatic cells of both males and females, and also in the cytoplasm of oocytes. These results suggested that gsdf might play an important role as initial switches to promote testis differentiation and participate in early germ cell development, such as proliferation and differentiation of spermatogonia and oogonia in P. olivaceus.

Potential contributions of heat shock proteins and related genes in sexual differentiation in yellow catfish (Pelteobagrus fulvidraco)

Sex determination and differentiation in ectotherms are very complicated affairs and usually affected by both genetic and environmental factors. Because of their temperature-sensitive expression, heat shock proteins (HSPs) are good candidates for temperature-dependent sex determination (TSD). Similar to most thermosensitive fish species, the male to female ratio increases with temperature in yellow catfish (Pelteobagrus fulvidraco). Yellow catfish is also a type of sexual size dimorphic fish, and the male individuals grow much faster than females of the same age. Therefore, research of sex differentiation in yellow catfish is important in aquiculture. In this attempt, a total of seven HSPs and related genes were identified from transcriptomes of yellow catfish by 454 pyrosequencing and Solexa sequencing that we did previously, including five genes with complete open reading frame (ORF). Phylogenetically, all these genes were compared with their counterparts from other vertebrates. All these genes were sex-biased expressed in gonads. Hspa5, Hip, and Cdc37 were expressed more highly in ovary than in testis, whereas Hsp90α, Hspb2, Hspb8, and Hspbp1 were expressed more highly in testis than in ovary. Additionally, the expression of these genes was assessed after 17α-methyltestosterone (MT) and 17α-ethinylestradiol (EE2) treatment, respectively. Our result showed that working as co-chaperones, these HSPs and related genes may regulate sex steroid receptor activities to influence gonad development in yellow catfish. Our work would help in the understanding of the mechanism of sexual differentiation in teleosts.

Sequential, Divergent, and Cooperative Requirements of Foxl2a and Foxl2b in Ovary Development and Maintenance of Zebrafish

Foxl2 is essential for mammalian ovary maintenance. Although sexually dimorphic expression of foxl2 was observed in many teleosts, its role and regulative mechanism in fish remained largely unclear. In this study, we first identified two transcript variants of foxl2a and its homologous gene foxl2b in zebrafish, and revealed their specific expression in follicular layer cells in a sequential and divergent fashion during ovary differentiation, maturation, and maintenance. Then, homozygous foxl2a mutants (foxl2a-/-) and foxl2b mutants (foxl2b-/-) were constructed and detailed comparisons, such as sex ratio, gonadal histological structure, transcriptome profiling, and dynamic expression of gonadal development-related genes, were carried out. Initial ovarian differentiation and oocyte development occur normally both in foxl2a-/- and foxl2b-/- mutants, but foxl2a and foxl2b disruptions result in premature ovarian failure and partial sex reversal, respectively, in adult females. In foxl2a-/- female mutants, sox9a-amh/cyp19a1a signaling was upregulated at 150 days postfertilization (dpf) and subsequently oocyte apoptosis was triggered after 180 dpf. In contrast, dmrt1 expression was greater at 105 dpf and increased several 100-fold in foxl2b-/- mutated ovaries at 270 dpf, along with other testis-related genes. Finally, homozygous foxl2a-/-/foxl2b-/- double mutants were constructed in which complete sex reversal occurs early and testis-differentiation genes robustly increase at 60 dpf. Given mutual compensation between foxl2a and foxl2b in foxl2b-/- and foxl2a-/- mutants, we proposed a model in which foxl2a and foxl2b cooperate to regulate zebrafish ovary development and maintenance, with foxl2b potentially having a dominant role in preventing the ovary from differentiating as testis, as compared to foxl2a.

Transcriptomic responses to environmental temperature by turtles with temperature-dependent and genotypic sex determination assessed by RNAseq inform the genetic architecture of embryonic gonadal development

Vertebrate sexual fate is decided primarily by the individual's genotype (GSD), by the environmental temperature during development (TSD), or both. Turtles exhibit TSD and GSD, making them ideal to study the evolution of sex determination. Here we analyze temperature-specific gonadal transcriptomes (RNA-sequencing validated by qPCR) of painted turtles (Chrysemys picta TSD) before and during the thermosensitive period, and at equivalent stages in soft-shell turtles (Apalone spinifera-GSD), to test whether TSD's and GSD's transcriptional circuitry is identical but deployed differently between mechanisms. Our data show that most elements of the mammalian urogenital network are active during turtle gonadogenesis, but their transcription is generally more thermoresponsive in TSD than GSD, and concordant with their sex-specific function in mammals [e.g., upregulation of Amh, Ar, Esr1, Fog2, Gata4, Igf1r, Insr, and Lhx9 at male-producing temperature, and of β-catenin, Foxl2, Aromatase (Cyp19a1), Fst, Nf-kb, Crabp2 at female-producing temperature in Chrysemys]. Notably, antagonistic elements in gonadogenesis (e.g., β-catenin and Insr) were thermosensitive only in TSD early-embryos. Cirbp showed warm-temperature upregulation in both turtles disputing its purported key TSD role. Genes that may convert thermal inputs into sex-specific development (e.g., signaling and hormonal pathways, RNA-binding and heat-shock) were differentially regulated. Jak-Stat, Nf-κB, retinoic-acid, Wnt, and Mapk-signaling (not Akt and Ras-signaling) potentially mediate TSD thermosensitivity. Numerous species-specific ncRNAs (including Xist) were differentially-expressed, mostly upregulated at colder temperatures, as were unannotated loci that constitute novel TSD candidates. Cirbp showed warm-temperature upregulation in both turtles. Consistent transcription between turtles and alligator revealed putatively-critical reptilian TSD elements for male (Sf1, Amh, Amhr2) and female (Crabp2 and Hspb1) gonadogenesis. In conclusion, while preliminary, our data helps illuminate the regulation and evolution of vertebrate sex determination, and contribute genomic resources to guide further research into this fundamental biological process.

Genome editing reveals dmrt1 as an essential male sex-determining gene in Chinese tongue sole (Cynoglossus semilaevis)

Chinese tongue sole is a marine fish with ZW sex determination. Genome sequencing suggested that the Z-linked dmrt1 is a putative male determination gene, but direct genetic evidence is still lacking. Here we show that TALEN of dmrt1 efficiently induced mutations of this gene. The ZZ dmrt1 mutant fish developed ovary-like testis, and the spermatogenesis was disrupted. The female-related genes foxl2 and cyp19a1a were significantly increased in the gonad of the ZZ dmrt1 mutant. Conversely, the male-related genes Sox9a and Amh were significantly decreased. The dmrt1 deficient ZZ fish grew much faster than ZZ male control. Notably, we obtained an intersex ZW fish with a testis on one side and an ovary on the other side. This fish was chimeric for a dmrt1 mutation in the ovary, and wild-type dmrt1 in the testis. Our data provide the first functional evidence that dmrt1 is a male determining gene in tongue sole.

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.

Comparative RNA-Seq analysis of differentially expressed genes in the testis and ovary of Takifugu rubripes

Takifugu rubripes is a classical model organism for studying the role of gonad organogenesis in such physiological processes as fertilization, sex determination, and sexual development. To explicitly investigate the mechanism associated with gonad organogenesis in T. rubripes, we obtained 44.3 million and 55.2 million raw reads from the testis and ovary, respectively, by RNA-seq and from these, 18,523 genes were identified. A total of 680 differentially expressed genes were obtained from comparison transcriptome analysis between the testis and ovary, and of these, 556 genes were up-regulated in the testis, whereas only 124 genes were upregulated in the ovary. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated that many of these genes encode proteins involved in signal transduction and gonad development. We mainly focused on the differentially expressed genes that have the potential to develop into the gonad. The generation of large scale transcriptomic data presented in this work would enrich the genetic resources of T. rubripes, which should be valuable to the comparative and evolutionary studies of teleosts.

Dmrt1 is necessary for male sexual development in zebrafish

The dmrt1 (doublesex and mab-3 related transcription factor 1) gene is a key regulator of sex determination and/or gonadal sex differentiation across metazoan animals. This is unusual given that sex determination genes are typically not well conserved. The mechanisms by which zebrafish sex is determined have remained elusive due to the lack of sex chromosomes and the complex polygenic nature of sex determination in domesticated strains. To investigate the role of dmrt1 in zebrafish sex determination and gonad development, we isolated mutations disrupting this gene. We found that the majority of dmrt1 mutant fish develop as fertile females suggesting a complete male-to-female sex reversal in mutant animals that would have otherwise developed as males. A small percentage of mutant animals became males, but were sterile and displayed testicular dysgenesis. Therefore zebrafish dmrt1 functions in male sex determination and testis development. Mutant males had aberrant gonadal development at the onset of gonadal sex-differentiation, displaying reduced oocyte apoptosis followed by development of intersex gonads and failed testis morphogenesis and spermatogenesis. By contrast, female ovaries developed normally. We found that Dmrt1 is necessary for normal transcriptional regulation of the amh (anti-Müllerian hormone) and foxl2 (forkhead box L2) genes, which are thought to be important for male or female sexual development respectively. Interestingly, we identified one dmrt1 mutant allele that co-operates with a linked segregation distorter locus to generate an apparent XY sex determination mechanism. We conclude that dmrt1 is dispensable for ovary development but necessary for testis development in zebrafish, and that dmrt1 promotes male development by transcriptionally regulating male and female genes as has been described in other animals. Furthermore, the strong sex-ratio bias caused by dmrt1 reduction-of-function points to potential mechanisms through which sex chromosomes may evolve.