Dmrt1 directly regulates the transcription of the testis-biased Sox9b gene in Nile tilapia (Oreochromis niloticus)

Decoding the fish genome opens a new era in important trait research and molecular breeding in China

Aquaculture represents the fastest-growing global food production sector, as it has become an essential component of the global food supply. China has the world's largest aquaculture industry in terms of production volume. However, the sustainable development of fish culture is hindered by several concerns, including germplasm degradation and disease outbreaks. The practice of genomic breeding, which relies heavily on genome information and genotypephenotype relationships, has significant potential for increasing the efficiency of aquaculture production. In 2014, the completion of the genome sequencing and annotation of the Chinese tongue sole signified the beginning of the fish genomics era in China. Since then, domestic researchers have made dramatic progress in functional genomic studies. To date, the genomes of more than 60 species of fish in China have been assembled and annotated. Based on these reference genomes, evolutionary, comparative, and functional genomic studies have revolutionized our understanding of a wide range of biologically and economically important traits of fishes, including growth and development, sex determination, disease resistance, metamorphosis, and pigmentation. Furthermore, genomic tools and breeding techniques such as SNP arrays, genomic selection, and genome editing have greatly accelerated genetic improvement through the incorporation of functional genomic information into breeding activities. This review aims to summarize the current status, advances, and perspectives of the genome resources, genomic study of important traits, and genomic breeding techniques of fish in China. The review will provide aquaculture researchers, fish breeders, and farmers with updated information concerning fish genomic research and breeding technology. The summary will help to promote the genetic improvement of production traits and thus will support the sustainable development of fish aquaculture.

A New Investigation to Discriminate Sexes in Alive Nile Tilapia (Oreochromis niloticus) Using Cyp19a1a and Dmrt1 Gene Expression in Tail Fin Tissues

The Nile Tilapia (Oreochromis niloticus), a gonochoristic teleost fish with a XX/XY sex-determination system, is an ideal model for investigating gonadal sex differentiation. During gonadal differentiation, the expression of cyp19a1a in XX gonads and dmrt1 in XY gonads are required for undifferentiated tissues to develop into ovary or testis. In this study, quantitative real-time RT-PCR assessed the expression of cyp19a1a and dmrt1 genes in gonads and tail fin tissues. Differences in gene expression mean among sexually differentiated fish were analyzed using two-way analysis of variance (ANOVA) and validation of mixed model using discriminant analysis (DA) for morphometric traits and the gene expression in gonads and tail fin tissues used to validate and utilize them in discriminating sexes in sex-differentiated Nile Tilapia fish. The results revealed that, cyp19a1a gene expression in female ovaries was more significant than dmrt1 in male testis. In the other hand, the dmrt1 gene expression in the tail fin was higher in males than females. Both, cyp19a1a and dmrt1 genes, can discriminate fish sexes by 100% by using their expression in tail fin tissues. In conclusion, the cyp19a1a and dmrt1 genes could be used as a genetic marker to discriminate between the Nile Tilapia sexes, whereas used as an indicator for ovarian or testis differentiation in sexually differentiated Nile Tilapia using tail fin tissues. It is worth mentioning that this is the first investigation for using cyp19a1a and dmrt1 genes from Nile Tilapia tail fin tissues in sex determination.

Gonadal transcriptome analysis of sex-biased gene and genome-wide investigation of dmrt gene family in Acanthogobius ommaturus

Acanthogobius ommaturus is one of the largest goby fish, and widely distributed in the Northwestern Pacific as an annual benthic fish. This study aims to report the gonadal transcriptome of A. ommaturus and identify differentially expressed genes (DEGs) between sexes. A total of 5460 (27.94 %) DEGs were detected from genome, with 3301 (16.89 %) biased towards males and 2159 (11.05 %) towards females. Analysis of 76 known vertebrate sex-related genes revealed multiple key genes, including the male-biased genes dmrt1 (Doublesex and Mab-3 related transcription factor 1) and amh (Anti-Mullerian Hormone), and the female-biased genes foxl2 (Forkhead Box L2) and cyp19a1a (Cytochrome P450 Aromatase 19 Subfamily A1). Furthermore, a genome-wide gene family analysis focused on the most significantly differentially expressed male-biased gene, dmrt1, was conducted using the chromosomal-level genome. Six Aodmrt genes were identified and subjected to phylogenetic and protein interaction network analyses. To validate the expression pattern, quantitative real-time PCR (qRT-PCR) was performed and compared with gonadal transcriptome data. The results showed that only dmrt1 exhibited significant male-bias, while the expression levels and sex differences of other dmrt genes in the gonads were inconclusive. Interestingly, the other dmrt genes displayed higher expression levels in other tissues, suggesting currently unknown functions. In conclusion, this study provides valuable genetic information contributing to the understanding of the sex determination mechanism of A. ommaturus and bony fish.

Genome-wide investigation of the DMRT gene family sheds new insight into the regulation of sex differentiation in spotted knifejaw (Oplegnathus punctatus) with fusion chromosomes (Y)

The role of the DMRT family in male sex determination and differentiation is significant, but its regulatory role in spotted knifejaw with Y fusion chromosomes remains unclear. Through genome-wide scanning, transcriptome analysis, qPCR, FISH, and RNA interference (RNAi), we investigated the DMRT family and the dmrt1-based sex regulation network. Seven DMRTs were identified (DMRT1/2 (2a,2b)/6, DMRT4/5, DMRT3), and dmrt gene dispersion among chromosomes is possibly driven by three whole-genome duplications. Transcriptome analysis enriched genes were associated with sex regulation and constructed a network associated with dmrt1. qPCR and FISH results showed the expression dimorphism of sex-related genes in dmrt-related regulatory networks. RNAi experiments indicated a distinct sex regulation mode in spotted knifejaw. Dmrt1 knockdown upregulated male-related genes (sox9a, sox9b, dmrt1, amh, amhr2) and hsd11b2 expression, which is critical for androgen synthesis. Amhr2 is located on the heterozygous chromosome (Y) and is specifically localized in primary spermatocytes, and is extremely upregulated after dmrt1 knockdown which suggested besides the important role of dmrt1 in male differentiation, the amhr2 along with amhr2/amh system, also play important regulatory roles in maintaining high expression of the hsd11b2 and male differentiation. This study aims to further investigate sex regulatory mechanisms in species with fusion chromosomes.

Transcriptome characterization of gonadal sex differentiation in Pacific bluefin tuna, Thunnus orientalis (Temminck et Schlegel)

Tunas (genus Thunnus) are one of the most ecologically and commercially important fish worldwide. To establish a biological basis for reproduction in this globally essential species, we have recently studied crucial reproductive aspects of the Pacific bluefin tuna (T. orientalis; PBT), as a model of tuna species, based on our closed-cycle aquaculture technology. In this study, we clarified the global expression profile of the genes regulating gonadal sex differentiation in PBT, as this developmental process is vital to sexual reproduction. Based on the results of our comparative (RNA-sequencing) and temporal (qRT-PCR) transcriptome analyses using the updated genome dataset, we propose the molecular mechanisms of gonadal sex differentiation in PBT. In female gonads, foxl2 and cyp19a1a (coding aromatase) are expressed at the onset of sex differentiation. Active aromatase-mediated estrogen biosynthesis, which includes positive regulation of cyp19a1a expression by Foxl2, induces ovarian differentiation. By contrast, dmrt1 and gsdf are upregulated in differentiating male gonads lacking active estrogen synthesis. Dmrt1 and Gsdf would mainly promote testicular differentiation. Furthermore, androgen biosynthesis is upregulated in differentiating male gonad. Endogenous androgens may also be vital to testicular differentiation. This study provides the first comprehensive data clarifying the molecular basis for gonadal sex differentiation in tunas.

Polyvinyl chloride microplastics induce changes in gene expression and organ histology along the HPG axis in Cyprinus carpio var. larvae

The negative consequences of microplastics pollution on the health of aquatic species have garnered extensive attention. However, the mechanisms through which microplastics may cause harm in the reproductive processes of fish remain unknown. For this study, Cyprinus carpio var. was subjected to four treatments with various concentrations of PVC microplastics for 60 days, through food rationed diets (no plastic control, 10%, 20% and 30%). The gonadosomatic indices, gonad and brain histologies, sex hormone levels, and transcriptional and translational genes in the hypothalamic-pituitary-gonadal (HPG) axes of both sexes were observed. According to the results, the gonadosomatic indices were significantly decreased, gonadal development was delayed, and the level of estradiol (E2) in the females was significantly elevated. In addition, the expression levels of genes associated with the HPG axis in the brains and gonads (gnrh, gtha1, fshβ, cyp19b, erα, vtg1, dmrt1, sox9b, and cyp19a) and the transcription levels of apoptosis-related genes in the brains and gonads (caspase3, bax, and bcl-2) exhibited significant changes. Further investigation revealed that the translation levels of genes linked to sex differentiation and sex steroid hormone (cyp19b and dmrt1) were significantly altered. These findings indicated that PVC likely microplastics may have a negative impact on the reproductive system of Cyprinus carpio var. by inhibiting gonadal development, affecting the gonad and brain structures, and altering the levels of steroid hormones and the expression of HPG axis-related genes. This work provides new insights into the toxicity of microplastics in aquatic organisms by revealing that PVC microplastics are a potential threat against the reproduction of fish populations.

Denovo RNA-Seq analysis of ovary and testis reveals potential differentially expressed transcripts associated with gonadal unsynchronization development in Onychostoma macrolepis

The Onychostoma macrolepis (O. macrolepis) is a rare and endangered wild species. Their endangered extinction might be due to their low fertility. To further illustrate the molecular mechanism of gonad development of the male and female O. macrolepis, the present study carried out de novo testicular and ovarian transcriptome sequencing. By comparing ovary and testis, 30,869 differentially expressed unigenes (9870 in female, 20999 in male) were identified. In addition, KEGG and GO analysis suggested that the Hedgehog signaling pathway have important roles in testis maintenance and spermatogenesis, whereas the Hippo signaling pathway and Wnt signaling pathway are likely to participate in ovary maintenance. RT-qPCR analysis results were consistent with transcriptome sequencing that all of gender differentiation-related genes (FOXL2, GDF9, WNT4, CYP19A1, SOX9 and GATA4), temperature-enriched genes (NOVA1, CTGF and NR4A1), clock-related genes (PER2, PER3, CRY1, CRY2, BMAL1 and CIPC) were significantly differential expression in testis compared with ovaries. Taken together, these results revealed a potential molecular mechanism that low fertility of the O. macrolepis might strong correlate with the gonadal dyssynchrony development of the male and female, which might provide theoretical basis and technical support for artificial reproduction and germplasm resource protection of the O. macrolepis.

Master-Key Regulators of Sex Determination in Fish and Other Vertebrates-A Review

In vertebrates, mainly single genes with an allele ratio of 1:1 trigger sex-determination (SD), leading to initial equal sex-ratios. Such genes are designated master-key regulators (MKRs) and are frequently associated with DNA structural variations, such as copy-number variation and null-alleles. Most MKR knowledge comes from fish, especially cichlids, which serve as a genetic model for SD. We list 14 MKRs, of which dmrt1 has been identified in taxonomically distant species such as birds and fish. The identification of MKRs with known involvement in SD, such as amh and fshr, indicates that a common network drives SD. We illustrate a network that affects estrogen/androgen equilibrium, suggesting that structural variation may exert over-expression of the gene and thus form an MKR. However, the reason why certain factors constitute MKRs, whereas others do not is unclear. The limited number of conserved MKRs suggests that their heterologous sequences could be used as targets in future searches for MKRs of additional species. Sex-specific mortality, sex reversal, the role of temperature in SD, and multigenic SD are examined, claiming that these phenomena are often consequences of artificial hybridization. We discuss the essentiality of taxonomic authentication of species to validate purebred origin before MKR searches.

Identification and functional analysis of Dmrt1 gene and the SoxE gene in the sexual development of sea cucumber, Apostichopus japonicus

Members of the Doublesex and Mab-3-related transcription factor (Dmrt) gene family handle various vital functions in several biological processes, including sex determination/differentiation and gonad development. Dmrt1 and Sox9 (SoxE in invertebrates) exhibit a very conserved interaction function during testis formation in vertebrates. However, the dynamic expression pattern and functional roles of the Dmrt gene family and SoxE have not yet been identified in any echinoderm species. Herein, five members of the Dmrt gene family (Dmrt1, 2, 3a, 3b and 5) and the ancestor SoxE gene were identified from the genome of Apostichopus japonicus. Expression studies of Dmrt family genes and SoxE in different tissues of adult males and females revealed different expression patterns of each gene. Transcription of Dmrt2, Dmrt3a and Dmrt3b was higher expressed in the tube feet and coelomocytes instead of in gonadal tissues. The expression of Dmrt1 was found to be sustained throughout spermatogenesis. Knocking-down of Dmrt1 by means of RNA interference (RNAi) led to the downregulation of SoxE and upregulation of the ovarian regulator foxl2 in the testes. This indicates that Dmrt1 may be a positive regulator of SoxE and may play a role in the development of the testes in the sea cucumber. The expression level of SoxE was higher in the ovaries than in the testes, and knocking down of SoxE by RNAi reduced SoxE and Dmrt1 expression but conversely increased the expression of foxl2 in the testes. In summary, this study indicates that Dmrt1 and SoxE are indispensable for testicular differentiation, and SoxE might play a functional role during ovary differentiation in the sea cucumber.

Identification and Expression Analysis of Dsx and Its Positive Transcriptional Regulation of IAG in Black Tiger Shrimp (Penaeus monodon)

Doublesex (Dsx) is a polymorphic transcription factor of the DMRTs family, which is involved in male sex trait development and controls sexual dimorphism at different developmental stages in arthropods. However, the transcriptional regulation of the Dsx gene is largely unknown in decapods. In this study, we reported the cDNA sequence of PmDsx in Penaeus monodon, which encodes a 257 amino acid polypeptide. It shared many similarities with Dsx homologs and has a close relationship in the phylogeny of different species. We demonstrated that the expression of the male sex differentiation gene Dsx was predominantly expressed in the P. monodon testis, and that PmDsx dsRNA injection significantly decreased the expression of the insulin-like androgenic gland hormone (IAG) and male sex-determining gene while increasing the expression of the female sex-determining gene. We also identified a 5′-flanking region of PmIAG that had two potential cis-regulatory elements (CREs) for the PmDsx transcription. Further, the dual-luciferase reporter analysis and truncated mutagenesis revealed that PmDsx overexpression significantly promoted the transcriptional activity of the PmIAG promoter via a specific CRE. These results suggest that PmDsx is engaged in male reproductive development and positively regulates the transcription of the PmIAG by specifically binding upstream of the promoter of the PmIAG. It provides a theoretical basis for exploring the sexual regulation pathway and evolutionary dynamics of Dmrt family genes in P. monodon.

Sex determination mechanisms and sex control approaches in aquaculture animals

Aquaculture is one of the most efficient modes of animal protein production and plays an important role in global food security. Aquaculture animals exhibit extraordinarily diverse sexual phenotypes and underlying mechanisms, providing an ideal system to perform sex determination research, one of the important areas in life science. Moreover, sex is also one of the most valuable traits because sexual dimorphism in growth, size, and other economic characteristics commonly exist in aquaculture animals. Here, we synthesize current knowledge of sex determination mechanisms, sex chromosome evolution, reproduction strategies, and sexual dimorphism, and also review several approaches for sex control in aquaculture animals, including artificial gynogenesis, application of sex-specific or sex chromosome-linked markers, artificial sex reversal, as well as gene editing. We anticipate that better understanding of sex determination mechanisms and innovation of sex control approaches will facilitate sustainable development of aquaculture.

The transcription factor Sox30 is involved in Nile tilapia spermatogenesis

Spermatogenesis is a complex process in which spermatogonial stem cells differentiate and develop into mature spermatozoa. The transcriptional regulatory network involved in fish spermatogenesis remains poorly understood. Here, we demonstrate in Nile tilapia that the Sox transcription factor family member Sox30 is specifically expressed in the testes and mainly localizes to spermatocytes and spermatids. CRISPR/Cas9-mediated sox30 mutation results in abnormal spermiogenesis, reduction of sperm motility, and male subfertility. Comparative transcriptome analysis shows that sox30 mutation alters the expression of genes involved in spermatogenesis. Further chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq), ChIP-PCR, and luciferase reporter assays reveal that Sox30 positively regulates the transcription of ift140 and ptprb, two genes involved in spermiogenesis, by directly binding to their promoters. Taken together, our data indicate that Sox30 plays essential roles in Nile tilapia spermatogenesis by directly regulating the transcription of the spermiogenesis-related genes ift140 and ptprb.

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 (O_II), ovary at the yolk vesicle stage (O_IV),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 O_IV vs O_II, while there were 3653 genes that were up-regulated and 2874 genes were down-regulated in comparison of O_IV vs T, there were 3345 genes that were up-regulated and 2995 genes were down-regulated in comparison of the O_II 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 O_II and O_IV, while significantly increased during the ovotestis in the hermaphroditic period (O_T), 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 (O_T). 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.

Sex Determination and Differentiation in Teleost: Roles of Genetics, Environment, and Brain

The fish reproductive system is a complex biological system. Nonetheless, reproductive organ development is conserved, which starts with sex determination and then sex differentiation. The sex of a teleost is determined and differentiated from bipotential primordium by genetics, environmental factors, or both. These two processes are species-specific. There are several prominent genes and environmental factors involved during sex determination and differentiation. At the cellular level, most of the sex-determining genes suppress the female pathway. For environmental factors, there are temperature, density, hypoxia, pH, and social interaction. Once the sexual fate is determined, sex differentiation takes over the gonadal developmental process. Environmental factors involve activation and suppression of various male and female pathways depending on the sexual fate. Alongside these factors, the role of the brain during sex determination and differentiation remains elusive. Nonetheless, GnRH III knockout has promoted a male sex-biased population, which shows brain involvement during sex determination. During sex differentiation, LH and FSH might not affect the gonadal differentiation, but are required for regulating sex differentiation. This review discusses the role of prominent genes, environmental factors, and the brain in sex determination and differentiation across a few teleost species.

Establishment of a stem Leydig cell line capable of 11-ketotestosterone production

The deficiency or insufficiency of androgen can trigger a range of reproductive diseases as well as other symptoms. Stem Leydig cells (SLCs) are critical for the formation and maintenance of a functional androgen-producing cell (Leydig cell, LC) population throughout adult male life. However, to date, our knowledge about SLCs is poor. Here we report the derivation and characterisation of a clonal stem LC line (designated as TSL) capable of 11- ketotestosterone (11-KT) production from a 3-month-old Nile tilapia (Oreochromis niloticus) testis. The cells retained stable proliferation after 77 generations with normal karyotype and growth factor dependency. They expressed platelet-derived growth factor receptor-α (pdgfrα), nestin and chicken ovalbumin upstream promoter transcription factor II (coup-tfIIa), which are characteristic of SLCs. Upon induction in defined medium, TSLs could undergo differentiation into steroidogenically active LCs and produce 11-KT. When implanted into recipient Nile tilapia testes from which endogenous LCs had been eliminated by ethane dimethanesulphonate (EDS) treatment, the PKH26-labelled TSLs could colonise the interstitium, subsequently express steroidogenic genes and restore 11-KT production. Taken together, our data suggest that TSLs possess the ability of continuous proliferation and potential of differentiation into functional LCs invitro and invivo. To the best of our knowledge TSL might represent the first stem LC line capable of 11-KT production to date. Our study may offer new opportunities for investigating the self-renewal of SLCs and steroidogenesis invitro, and provide an invaluable invitro model for investigating endocrine disruptors.

Germline sexual fate is determined by the antagonistic action of dmrt1 and foxl3/foxl2 in tilapia

Germline sexual fate has long been believed to be determined by the somatic environment, but this idea is challenged by recent studies of foxl3 mutants in medaka. Here, we demonstrate that the sexual fate of tilapia germline is determined by the antagonistic interaction of dmrt1 and foxl3, which are transcriptionally repressed in male and female germ cells, respectively. Loss of dmrt1 rescued the germ cell sex reversal in foxl3Δ7/Δ7 XX fish, and loss of foxl3 partially rescued germ cell sex reversal but not somatic cell fate in dmrt1Δ5/Δ5 XY fish. Interestingly, germ cells lost sexual plasticity in dmrt1Δ5/Δ5 XY and foxl3Δ7/Δ7 XX single mutants, as aromatase inhibitor (AI) and estrogen treatment failed to rescue the respective phenotypes. However, recovery of germ cell sexual plasticity was observed in dmrt1/foxl3 double mutants. Importantly, mutation of somatic cell-specific foxl2 resulted in testicular development in foxl3Δ7/Δ7 or dmrt1Δ5/Δ5 mutants. Our findings demonstrate that sexual plasticity of germ cells relies on the presence of both dmrt1 and foxl3. The existence of dmrt1 and foxl3 allows environmental factors to influence the sex fate decision in vertebrates.

Comparison of Gonadal Transcriptomes Uncovers Reproduction-Related Genes with Sexually Dimorphic Expression Patterns in Diodon hystrix

Diodon hystrix is a new and emerging aquaculture species in south China. However, due to the lack of understanding of reproductive regulation, the management of breeding and reproduction under captivity remains a barrier for the commercial aquaculture of D. hystrix. More genetic information is needed to identify genes critical for gonadal development. Here, the first gonadal transcriptomes of D. hystrix were analyzed and 151.89 million clean reads were generated. All reads were assembled into 57,077 unigenes, and 24,574 could be annotated. By comparing the gonad transcriptomes, 11,487 differentially expressed genes were obtained, of which 4599 were upregulated and 6888 were downregulated in the ovaries. Using enrichment analyses, many functional pathways were found to be associated with reproduction regulation. A set of sex-biased genes putatively involved in gonad development and gametogenesis were identified and their sexually dimorphic expression patterns were characterized. The detailed transcriptomic data provide a useful resource for further research on D. hystrix reproductive manipulation.

Differential expression analysis and identification of sex-related genes by gonad transcriptome sequencing in estradiol-treated and non-treated Ussuri catfish Pseudobagrus ussuriensis

The Ussuri catfish (Pseudobagrus ussuriensis) has an XX/XY sex determination system but its sex determination gene(s) remain unknown. To better understand the molecular sex determination mechanism, transcriptome analysis was conducted to obtain sex-related gene expression profiles. Transcriptome analyses were made of male and female developing/differentiating gonads by high-throughput RNA sequencing, including gonads from fish given an estradiol-induced sex reversal treatment. A total of 81,569 unigenes were assembled and 39,904 were significantly matched to known unique proteins by comparison with public databases. Twenty specifically expressed and 142 differentially expressed sex-related genes were extracted from annotated data by comparing the treatment groups. These genes are involved in spermatogenesis (e.g., Dnali1, nectin3, klhl10, mybl1, Katnal1, Eno4, Mns1, Spag6, Tsga10, Septin7), oogenesis (e.g., Lagr5, Fmn2, Npm2, zar1, Fbxo5, Fbxo43, Prdx4, Nrip1, Lfng, Atrip), gonadal development/differentiation (e.g., Cxcr4b, Hmgb2, Cftr, Ch25h, brip1, Prdm9, Tdrd1, Star, dmrt1, Tut4, Hsd17b12a, gdf9, dnd, arf1, Spata22), and estradiol response (e.g., Mmp14, Lhcgr, vtg1, vtg2, esr2b, Piwil1, Aifm1, Hsf1, gdf9). Dmrt1 and gdf9 may play an essential role in sex determination in P. ussuriensis. The expression patterns of six random genes were validated by quantitative real-time PCR, which confirmed the reliability and accuracy of the RNA-seq results. These data provide a valuable resource for future studies of gene expression and for understanding the molecular mechanism of sex determination/differentiation and gonadal development/differentiation (including hormone-induced sexual reversal) in Ussuri catfish. This has the potential to assist in producing monosex Ussuri catfish to increase aquacultural productivity.

Molecular identification and functional analysis of MyD88 in giant freshwater prawn (Macrobrachium rosenbergii) and expression changes in response to bacterial challenge

Myeloid differentiation factor 88 (MyD88) is a crucial adaptor protein for Toll-like receptor (TLR)-mediated signaling pathways and plays an important role in immune response. In this study, the full-length cDNA of MyD88 from Macrobrachium rosenbergii (MRMyD88) was cloned. The MRMyD88 cDNA is 1758 bp long and contains a 1398-bp open reading frame. Multiple sequence alignment and phylogenetic analysis revealed that the amino acid sequence of MRMyD88 shared high identity with the known MyD88 proteins. The MRMyD88 mRNA was widely expressed in all examined tissues, with highest level in intestine, followed by gonad and pleopod. Furthermore, the MRMyD88 promoter region, spanning 1622 bp, contains several transcription factor-binding sites, including nine GATA-1 box motifs. Electrophoretic mobility shift assay showed that Gfi-1, SRF, and Oct-1 bind to the upstream region of MRMyD88. Additionally, the results showed that the expression levels of TLR1, TLR2 and TLR3 were different in response to Vibrio anguillarum, Lactobacillus plantarum and Aeromonas hydrophila infections. However, these bacteria significantly increased the expression levels of MyD88 and prophenoloxidase. These data suggest that the TLR-mediated signaling pathway is MyD88-dependent in response to pathogenic and probiotic bacteria in M. rosenbergii.

Molecular cloning, characterization, and expression analysis of a sex-biased transcriptional factor sox9 gene of mud crab Scylla paramamosain

Sox9 gene, a crucial member of the Sox gene family, is present in various organisms and involved in many physiological processes, especially in sex determination and gonad development. In this study, we reported a sox9 gene (designated as Spsox9) from Scylla paramamosain through analyzing published gonad transcriptome data. Meanwhile, the accuracy was validated by PCR technology, and the 3' sequences were cloned with 3' RACE technology. The full-length cDNA of Spsox9 is 2843 bp, consisting of a 243 bp 5' UTR, an 1124 bp 3' UTR, and a 1476 bp ORF encoding 491 amino acids. Furthermore, to better understand its conservation among crustacean species, the sox9 gene ortholog was identified in several other crustaceans species with their published transcriptome data, respectively. All of the Sox9 proteins identified in the current study had the common feature of Sox proteins (HMG domain) and were highly conserved among analyzed crustacean species. In all examined tissues, the Spsox9 was mainly expressed in the gonad (testis and ovary), eyestalk, and cerebral ganglion. During embryo development, Spsox9 was highly expressed in 5 pairs of appendages, 7 pairs of appendages, and eye-pigment formation stage. During ovary development, the expression level of Spsox9 remained stable in the first 4 stages (O1-O4) and decreased in the tertiary vitellogenesis (O5) stage. During testis development, the expression level of Spsox9 was highest in the spermatid stage (T2) and was significantly different from that in the spermatocyte stage (T1) and mature sperm stage (T3) (p < 0.05). In addition, Spsox9 exhibited a sex-biased expression pattern in T1 and O1. These present results indicated that the Spsox9 gene might play crucial roles in the gonad and embryo development of mud crab.

Comparative Genomics Studies on the dmrt Gene Family in Fish

Doublesex and mab-3-related transcription factor (dmrt) genes are widely distributed across various biological groups and play critical roles in sex determination and neural development. Here, we applied bioinformatics methods to exam cross-species changes in the dmrt family members and evolutionary relationships of the dmrt genes based on genomes of 17 fish species. All the examined fish species have dmrt1-5 while only five species contained dmrt6. Most fish harbored two dmrt2 paralogs (dmrt2a and dmrt2b), with dmrt2b being unique to fish. In the phylogenetic tree, 147 DMRT are categorized into eight groups (DMRT1-DMRT8) and then clustered in three main groups. Selective evolutionary pressure analysis indicated purifying selections on dmrt1-3 genes and the dmrt1-3-2(2a) gene cluster. Similar genomic conservation patterns of the dmrt1-dmrt3-dmrt2(2a) gene cluster with 20-kb upstream/downstream regions in fish with various sex-determination systems were observed except for three regions with remarkable diversity. Synteny analysis revealed that dmrt1, dmrt2a, dmrt2b, and dmrt3-5 were relatively conserved in fish during the evolutionary process. While dmrt6 was lost in most species during evolution. The high conservation of the dmrt1-dmrt3-dmrt2(2a) gene cluster in various fish genomes suggests their crucial biological functions while various dmrt family members and sequences across fish species suggest different biological roles during evolution. This study provides a molecular basis for fish dmrt functional analysis and may serve as a reference for in-depth phylogenomics.

High-quality chromosome-level genomes of two tilapia species reveal their evolution of repeat sequences and sex chromosomes

Tilapias are one of the most farmed fishes that are coined as "aquatic chicken" by the food industry. Nile tilapia and blue tilapia exhibit very recent transition of sex chromosome systems since their divergence approximately five million years ago, making them a great model for elucidating the molecular and evolutionary mechanisms of sex chromosome turnovers. Studies of their sex-determining pathways are also critical for developing genetic sex control in aquaculture. We report here the newly produced genomes of Nile tilapia and blue tilapia that integrate long-read sequencing and chromatin conformation data. The two nearly complete genomes have anchored over 97% of the sequences into linkage groups (LGs), and assembled majorities of complex repetitive regions including telomeres, centromeres and rDNA clusters. In particular, we inferred two episodes of repeat expansion at LG3 respectively in the ancestor of cichlids and that of tilapias. The consequential large heterochromatic region concentrated at one end of LG3 comprises tandem arrays of mRNA and small RNA genes, among which we have identified a candidate female determining gene Paics in blue tilapia. Paics shows female-specific patterns of single-nucleotide variants, copy numbers and expression patterns in gonads during early gonadogenesis. Our work provides a very important genomic resource for functional studies of cichlids, and suggested that unequal distribution of repeat content that impacts the local recombination rate might make some chromosomes more likely to become sex chromosomes.

A novel transcriptional cascade is involved in Fzr-mediated endoreplication

Endoreplication, known as endocycle, is a variant of the cell cycle that differs from mitosis and occurs in specific tissues of different organisms. Endoreplicating cells generally undergo multiple rounds of genome replication without chromosome segregation. Previous studies demonstrated that Drosophila fizzy-related protein (Fzr) and its mammalian homolog Cdh1 function as key regulators of endoreplication entrance by activating the anaphase-promoting complex/cyclosome to initiate the ubiquitination and subsequent degradation of cell cycle factors such as Cyclin B (CycB). However, the molecular mechanism underlying Fzr-mediated endoreplication is not completely understood. In this study, we demonstrated that the transcription factor Myc acts downstream of Fzr during endoreplication in Drosophila salivary gland. Mechanistically, Fzr interacts with chromatin-associated histone H2B to enhance H2B ubiquitination in the Myc promoter and promotes Myc transcription. In addition to negatively regulating CycB transcription, the Fzr-ubiquitinated H2B (H2Bub)-Myc signaling cascade also positively regulates the transcription of the MCM6 gene that is involved in DNA replication by directly binding to specific motifs within their promoters. We further found that the Fzr-H2Bub-Myc signaling cascade regulating endoreplication progression is conserved between insects and mammalian cells. Altogether, our work uncovers a novel transcriptional cascade that is involved in Fzr-mediated endoreplication.

Comparative transcriptome analysis on four types of gonadal tissues of blotched snakehead (Channa maculata)

Blotched snakehead (Channa maculata) is an economically important freshwater fish in China, of which males grow much faster than females. To illuminate the molecular mechanism of sex differentiation and gonad development, RNA-Sequencing was performed to identify sex-related genes and pathway in gonads of 6-month-old normal XX females (XX-F), normal XY males (XY-M), XY sex reversal females (XY-F) and YY super-males (YY-M). The analysis showed that many differentially expressed genes (DEGs) had similar expression patterns in XY-F and XX-F, which were different from XY-M and YY-M. qRT-PCR indicated that Amh, Dmrt1, and Sox9 had relatively high expression in testes of XY-M and YY-M. Taking Amh as an example, there was a relative fold change of 1.0 in XX-F, 2.1 fold change in XY-F, 36.1 fold change in XY-M, and 26.0 fold change in YY-M. Cyp19a1a, Figla, and Foxl2 were highly expressive in ovaries of XX-F and XY-F. Taking Figla as an example, there was a relative fold change of 557 in XX-F, 304.5 fold change in XY-F, 5.6 fold change in XY-M, and 4.4 fold change in YY-M. KEGG analysis revealed many DEGs distributed in pathways related to sex differentiation, steroid hormone synthesis and growth, etc. Significant variation and trends in relative expression levels tested by qRT-PCR were consistent with those recorded by RNA-Sequencing. This is the first time that transcriptome of snakehead has been investigated systematically and in an integrated way. Large quantities of candidate genes involved in sex differentiation, gonad development and growth dimorphism were identified. The study provides useful resources for understanding sex differentiation and growth dimorphism, potentially assisting mono-sex production of snakehead in aquaculture.

Genome-wide investigation of Dmrt gene family in large yellow croaker (Larimichthys crocea)

The Dmrt (Doublesex and Mab-3 related transcription factor) gene family is a class of crucial transcription factors characterized by a conserved DM (Doublesex/Mab-3) domain. Previous researches indicate this gene family is involved in various physiological processes, especially in sex determination/differentiation and gonad development. Despite the vital roles of the Dmrt gene family in physiological processes, the comprehensive characterization and analysis of the dmrt genes in large yellow croaker (Larimichthys crocea), one of the most commercially important marine fish in China, have not been described. In this study, we performed the first genome-wide systematic analysis of L. crocea dmrt genes through the bioinformatics method. A total of seven members of the Dmrt gene family including Lcdmrt1, Lcdmrt2a, Lcdmrt2b, Lcdmrt3, Lcdmrt4, Lcdmrt5, and Lcdmrt6 were excavated based on the genome data of L. crocea. Further analysis revealed that the dmrt genes of L. crocea were distributed unevenly across four chromosomes. There were three dmrt genes (Lcdmrt1, Lcdmrt2a, and Lcdmrt3) on 3rd chromosome, one (Lcdmrt6) on 13th chromosome, one (Lcdmrt4) on 14th chromosome, two on (Lcdmrt5 and Lcdmrt2b) 17th chromosome. The gene structure analysis indicated that the number of introns of different dmrt genes of L. crocea had some differences: Lcdmrt1 had four introns, Lcdmrt2a, Lcdmrt2b, and Lcdmrt6 had two introns, Lcdmrt3, Lcdmrt4, and Lcdmrt5 had only one intron. The expression pattern analysis with published gonad transcriptome datasets and further confirmed by qRT-PCR revealed that these members of the Dmrt gene family except for Lcdmrt4 were all sexually dimorphic and preferred expressing in testis. Furthermore, the expression pattern analysis also revealed that the expression level of Lcdmrt1 and Lcdmrt6 was significantly higher than that of other members, suggesting that these two genes may play a more important role in testis. Overall, our studies provide a comprehensive insight into the Dmrt gene family members and a basis for the further study of their biological functions in L. crocea.

Effects of 17β-estradiol on early gonadal development and expression of genes implicated in sexual differentiation of a South American teleost, Astyanax altiparanae

Gonadal sex differentiation in teleost fish shows greater plasticity as compared to other vertebrates, as it can be influenced by a variety of factors such as exogenous sex steroids. Exogenous estrogens, such as 17β-estradiol (E2), can induce feminization when administered during early embryonic development. However, the mechanisms underlying the E2-induced feminization are not fully understood, especially in Neotropical species. Therefore, the aim of this study was to evaluate the effects of E2 administration on the phenotypic sex characteristics, histological assessment of the gonads, and the expression of selected genes in Astyanax altiparanae exposed to dietary E2 prior to gonadal differentiation. At 4 days post-hatch (dph), groups of 30-40 undifferentiated larvae were fed with a diet containing varying amounts of E2 for 28 days, and fish were sampled at 90 dph. Previous studies revealed that ovary formation in A. altiparanae occurred at 58 dph, whereas the first sign of testis formation was found at 73 dph. In relation to the control, E2 exposure increased the proportion of phenotypic females in 120% and 148.4% for 4 and 6 mg E2/Kg, respectively. However, histological analysis revealed that treatments did not affect gonadal sex ratio between males and females, but induced intersex (testis-ova) in the group treated with 6 mg E2/Kg food. Treatment with E2 also altered gonadal transcript levels of a selected number of genes implicated in sexual differentiation. Males overexpressed dmrt1, sox9 and amh following E2 treatment as compared to control. Females showed increased mRNA levels of dmrt1 and sox9, which might be related to the down-regulation of cyp19a1a after E2 exposure. In summary, E2 exposure during early gonadal development affected male secondary characteristics without changing the gonadal sex ratio, and altered expression of genes implicated in sexual differentiation.

Spatial and temporal expression pattern of sex-related genes in ovo-testis of the self-fertilizing mangrove killifish (Kryptolebias marmoratus)

In vertebrates, sex determination and differentiation comprehend a fine balance between female and male factors, leading the bipotential anlage to develop towards ovary or testis, respectively. Nevertheless, the mangrove killifish, (Kryptolebias marmoratus) a simultaneous hermaphroditic species, could overcome those antagonistic pathways and evolved to develop and maintain reproductively active ovarian and testicular tissues in the same organ. Morphological and mRNA localization analyzes of developing and adult gonads demonstrate that genes related to testis (dmrt1 and amh) and ovary differentiation (foxl2 and sox9a) follow the same expression pattern observed in gonochoristic species, thus functioning as two independent organs. In addition, Amh expression patterns make it a strong candidate for initiation of the formation and maintenance of the testicular tissue in the hermaphroditic gonad. Differently from described so far, foxl3 seems to have an important role in oogenesis as well as spermatogenesis and gonadal structure.

An ex vivo Approach to Study Hormonal Control of Spermatogenesis in the Teleost Oreochromis niloticus

As the male reproductive organ, the main task of the testis is the production of fertile, haploid spermatozoa. This process, named spermatogenesis, starts with spermatogonial stem cells, which undergo a species-specific number of mitotic divisions until starting meiosis and further morphological maturation. The pituitary gonadotropins, luteinizing hormone, and follicle stimulating hormone, are indispensable for vertebrate spermatogenesis, but we are still far from fully understanding the complex regulatory networks involved in this process. Therefore, we developed an ex vivo testis cultivation system which allows evaluating the occurring changes in histology and gene expression. The experimental circulatory flow-through setup described in this work provides the possibility to study the function of the male tilapia gonads on a cellular and transcriptional level for at least 7 days. After 1 week of culture, tilapia testis slices kept their structure and all stages of spermatogenesis could be detected histologically. Without pituitary extract (tilPE) however, fibrotic structures appeared, whereas addition of tilPE preserved spermatogenic cysts and somatic interstitium completely. We could show that tilPE has a stimulatory effect on spermatogonia proliferation in our culture system. In the presence of tilPE or hCG, the gene expression of steroidogenesis related genes (cyp11b2 and stAR2) were notably increased. Other testicular genes like piwil1, amh, or dmrt1 were not expressed differentially in the presence or absence of gonadotropins or gonadotropin containing tilPE. We established a suitable system for studying tilapia spermatogenesis ex vivo with promise for future applications.

Genome-Wide identification of doublesex and Mab-3-Related transcription factor (DMRT) genes in nile tilapia (oreochromis niloticus)

Doublesex and Mab-3-related transcription factor (DMRT) gene family is extensively known for its contribution in sex determination and differentiation across phyla. Here we report the identification of five DM (doublesex and mab-3) domain genes in the Nile tilapia which includes DMRT1, DMRTa2, DMRT2a, DMRT2b and DMRT3a. The full-length sequence of DMRT genes ranges from 3526 (DMRTA2) to 1471bp (DMRT1) which encode putative proteins series from 469 to 372 amino acids. All the DMRT proteins contained at least one conserved DNA-binding DM domain. Sub-cellular localization and gene ontology revealed DMRT1 protein is maximum localized in nuclear region and gene ontology analysis showed the molecular function of 48.2%, biological process 43.6% and cellular component 25%. Chromosomal location and synteny analysis displayed that DMRT genes mostly cluster linkage group 12. Altogether, our findings provide vital genomic information for future studies of biochemical, physiological, and phylogenetic studies on DMRT genes in teleost.

Transcription of the Sox30 Gene Is Positively Regulated by Dmrt1 in Nile Tilapia

The Sox family member Sox30 is highly expressed in the testis of several vertebrate species and has been shown to play key roles in spermiogenesis. However, its transcription regulation remains unclear. Here, we analyzed the Sox30 promoter from the teleost fish Nile tilapia (Oreochromis niloticus) and predicted a putative cis-regulatory element (CRE) for doublesex and mab-3 related transcription factor 1 (Dmrt1), a male-specific transcription factor involved in male sex differentiation. Transcriptional profiling revealed that Sox30 and Dmrt1 similarly exhibited a high expression in tilapia testes from 90 days after hatching (dah) to 300 dah, and the transcription of the Sox30 gene was reduced about one-fold in the testes of male tilapia with Dmrt1 knockdown. Further dual-luciferase reporter assay confirmed that Dmrt1 overexpression significantly promoted transcriptional activity of the Sox30 promoter and this promotion was decreased following the mutation of putative CRE for Dmrt1 within the Sox30 promoter. Chromatin immunoprecipitation-based PCR (ChIP-PCR) and electrophoretic mobility shift assay (EMSA) demonstrated that Dmrt1 directly binds to putative CRE within the Sox30 promoter. These results together indicate that Dmrt1 positively regulates the transcription of the tilapia Sox30 gene by directly binding to specific CRE within the Sox30 promoter.

Transcriptomics analysis revealing candidate genes and networks for sex differentiation of yesso scallop (Patinopecten yessoensis)

Background

The Yesso scallop, Patinopecten (Mizuhopecten) yessoensis, is a commercially important bivalve in the coastal countries of Northeast Asia. It has complex modes of sex differentiation, but knowledge of the mechanisms underlying this sex determination and differentiation is limited.

Results

In this study, the gonad tissues from females and males at three developmental stages were used to investigate candidate genes and networks for sex differentiation via RNA-Req. A total of 901,980,606 high quality clean reads were obtained from 18 libraries, of which 417 expressed male-specific genes and 754 expressed female-specific genes. Totally, 10,074 genes differentially expressed in females and males were identified. Weighted gene co-expression network analysis (WGCNA) revealed that turquoise and green gene modules were significantly positively correlated with male gonads, while coral1 and black modules were significantly associated with female gonads. The most important gene for sex determination and differentiation was Pydmrt 1, which was the only gene discovered that determined the male sex phenotype during early gonadal differentiation. Enrichment analyses of GO terms and KEGG pathways revealed that genes involved in metabolism, genetic and environmental information processes or pathways are sex-biased. Forty-nine genes in the five modules involved in sex differentiation or determination were identified and selected to construct a gene co-expression network and a hypothesized sex differentiation pathway.

Conclusions

The current study focused on screening genes of sex differentiation in Yesso scallop, highlighting the potential regulatory mechanisms of gonadal development in P. yessoensis. Our data suggested that WCGNA can facilitate identification of key genes for sex differentiation and determination. Using this method, a hypothesized P. yessoensis sex determination and differentiation pathway was constructed. In this pathway, Pydmrt 1 may have a leading function.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6021-6) contains supplementary material, which is available to authorized users.