Dmrt1 (doublesex and mab-3-related transcription factor 1) expression during gonadal development and spermatogenesis in the Japanese eel

Full-length transcriptome analysis of male and female gonads in Japanese Eel (Anguilla japonica)

BACKGROUND

The Japanese eel (Anguilla japonica) holds significant economic value in East Asia, but limitations in understanding its reproductive biology have hindered advancements in artificial breeding techniques. Previous research has primarily focused on conserved sex differentiation genes, offering limited insights into the broader molecular mechanisms driving gonadal development and sexual dimorphism. To address these limitations, this study aims to investigate key genes and pathways involved in gonadal development through a comprehensive transcriptomic analysis of male and female eel gonads.

RESULTS

PacBio Iso-Seq and Illumina RNA-Seq technologies were combined to conduct a full-length transcriptome analysis of male and female Japanese eel gonads at a post-differentiation, pre-maturation stage. A total of 24,661 unigenes were identified in ovaries and 15,023 in testes, along with genomic regulatory elements such as transcription factors, simple sequence repeats, and long non-coding RNAs. Additionally, 1,210 differentially expressed genes were detected. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses revealed significant pathways involved in cell cycle regulation, metabolic processes, apoptosis, and hormone activity. Notably, several reproductive-related genes, including bambi, ccnb1, cdc20, gdf9, prlh, ccdc39, chrebp, tspo, syce3, and ngb, demonstrated significant dimorphic expression in eel gonads.

CONCLUSIONS

This study provides valuable insights into the molecular mechanisms of gonadal differentiation and sexual dimorphism in Japanese eels. The findings expand the genetic resources available for the eel breeding industry and could facilitate the development of improved artificial breeding techniques focused on reproductive development.

Understanding Genetic Regulation of Sex Differentiation in Hermaphroditic Fish

As a fundamental taxonomic group within vertebrates, fish represent an invaluable resource for investigating the mechanisms underlying sex determination and differentiation owing to their extensive geographical distribution and rich biodiversity. Within this biological cohort, the processes of sex determination and differentiation are intricately governed by both genetic factors and the complex interplay of environmental cues. While variations in external environmental factors, particularly temperature, can exert a modulatory influence on sex differentiation in fish to a limited degree, genetic factors remain the primary determinants of sexual traits. Hermaphroditic fish display three distinct types of sexual transitions: protandry (male to female), protogyny (female-to-male), bidirectional sex change (both directions serially). These fish, characterized by their unique reproductive strategies and sexual plasticity, serve as exemplary natural models for elucidating the mechanisms of sex differentiation and sexual transitions in fish. The present review delves into the histological dynamics during gonadal development across three types of sequential hermaphroditic fish, meticulously delineating the pivotal characteristics at each stage, from the inception of primordial gonads to sexual specialization. Furthermore, it examines the regulatory genes and associated signaling pathways that orchestrate sex determination and differentiation. By systematically synthesizing these research advancements, this paper endeavors to offer a comprehensive and profound insight into the intricate mechanisms governing sex differentiation in sequential hermaphroditic fish.

The expression profiles of cyp19a1, sf-1, esrs and gths in the brain-pituitary during gonadal sex differentiation in juvenile Japanese eels

Eels are gonochoristic species whose gonadal differentiation initiates at the yellow eel stage and is influenced by environmental factors. We revealed some sex-related genes were sex dimorphically expressed in gonads during gonadal sex differentiation of Japanese eel (Anguilla japonica); however, the expression of sex-related genes in the brain-pituitary during gonadal sex differentiation in eels is still unclear. This study aimed to investigate the sex-related gene expressions in the brain-pituitary and tried to clarify their roles in the brain and gonads during gonadal sex differentiation. Based on our previous histological study, the control eels developed as males, and estradiol-17β (E2) was used for feminization. Our results showed that during testicular differentiation, the brain cyp19a1 transcripts and aromatase proteins were increased significantly; moreover, the cyp19a1, sf-1, foxl2s, and esrs (except gperb) transcripts in the midbrain/pituitary also were increased significantly. Forebrain gnrh1 transcripts increased slightly during gonadal differentiation of both sexes, but the gnrhr1b and gnrhr2 transcripts in the midbrain/pituitary were stable during gonadal differentiation. The expression levels of gths and gh in the midbrain/pituitary were significantly increased during testicular differentiation and were much higher in males than in E2-feminized females. These results implied that endogenous estrogens might play essential roles in the brain/pituitary during testicular differentiation, sf-1, foxl2s, and esrs may have roles in cyp19a1 regulation in the midbrain/pituitary of Japanese eels. For the GnRH-GTH axis, gths, especially fshb, may be regulated by esrs and involved in regulating testicular differentiation and development in Japanese eels.

A Testis-Specific DMRT1 (Double Sex and Mab-3-Related Transcription Factor 1) Plays a Role in Spermatogenesis and Gonadal Development in the Hermaphrodite Boring Giant Clam Tridacna crocea

The testis-specific double sex and mab-3-related transcription factor 1 (DMRT1) has long been recognized as a crucial player in sex determination across vertebrates, and its essential role in gonadal development and the regulation of spermatogenesis is well established. Here, we report the cloning of the key spermatogenesis-related DMRT1 cDNA, named Tc-DMRT1, from the gonads of Tridacna crocea (T. crocea), with a molecular weight of 41.93 kDa and an isoelectric point of 7.83 (pI). Our hypothesis is that DMRT1 machinery governs spermatogenesis and regulates gonadogenesis. RNAi-mediated Tc-DMRT1 knockdown revealed its critical role in hindering spermatogenesis and reducing expression levels in boring giant clams. A histological analysis showed structural changes, with normal sperm cell counts in the control group (ds-EGFP) but significantly lower concentrations of sperm cells in the experimental group (ds-DMRT1). DMRT1 transcripts during embryogenesis exhibited a significantly high expression pattern (p < 0.05) during the early zygote stage, and whole-embryo in-situ hybridization confirmed its expression pattern throughout embryogenesis. A qRT-PCR analysis of various reproductive stages revealed an abundant expression of Tc-DMRT1 in the gonads during the male reproductive stage. In-situ hybridization showed tissue-specific expression of DMRT1, with a positive signal detected in male-stage gonadal tissues comprising sperm cells, while no signal was detected in other stages. Our study findings provide an initial understanding of the DMRT1 molecular machinery controlling spermatogenesis and its specificity in male-stage gonads of the key bivalve species, Tridacna crocea, and suggest that DMRT1 predominantly functions as a key regulator of spermatogenesis in giant clams.

Dmrt1 is the only male pathway gene tested indispensable for sex determination and functional testis development in tilapia

Sex is determined by multiple factors derived from somatic and germ cells in vertebrates. We have identified amhy, dmrt1, gsdf as male and foxl2, foxl3, cyp19a1a as female sex determination pathway genes in Nile tilapia. However, the relationship among these genes is largely unclear. Here, we found that the gonads of dmrt1;cyp19a1a double mutants developed as ovaries or underdeveloped testes with no germ cells irrespective of their genetic sex. In addition, the gonads of dmrt1;cyp19a1a;cyp19a1b triple mutants still developed as ovaries. The gonads of foxl3;cyp19a1a double mutants developed as testes, while the gonads of dmrt1;cyp19a1a;foxl3 triple mutants eventually developed as ovaries. In contrast, the gonads of amhy;cyp19a1a, gsdf;cyp19a1a, amhy;foxl2, gsdf;foxl2 double and amhy;cyp19a1a;cyp19a1b, gsdf;cyp19a1a;cyp19a1b triple mutants developed as testes with spermatogenesis via up-regulation of dmrt1 in both somatic and germ cells. The gonads of amhy;foxl3 and gsdf;foxl3 double mutants developed as ovaries but with germ cells in spermatogenesis due to up-regulation of dmrt1. Taking the respective ovary and underdeveloped testis of dmrt1;foxl3 and dmrt1;foxl2 double mutants reported previously into consideration, we demonstrated that once dmrt1 mutated, the gonad could not be rescued to functional testis by mutating any female pathway gene. The sex reversal caused by mutation of male pathway genes other than dmrt1, including its upstream amhy and downstream gsdf, could be rescued by mutating female pathway gene. Overall, our data suggested that dmrt1 is the only male pathway gene tested indispensable for sex determination and functional testis development in tilapia.

Tilapia, a good model for studying reproductive endocrinology

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

Transcription factors Dmrt1a, Foxl2, and Nr5a1a potentially interact to regulate cyp19a1a transcription in ovarian follicles of ricefield eel (Monopterus albus)

Aromatase (encoded by Cyp19a1) in the ovarian follicular cells catalyzes the production of estradiol from testosterone, which plays important roles in the ovarian development of vertebrates. In the present study, the interaction of Dmrt1, Foxl2, and Nr5a1a on the regulation of cyp19a1a transcription in ovarian follicles was examined in a teleost, the ricefield eel Monopterus albus. The expression of dmrt1a, foxl2, and nr5a1a was detected in ovarian follicular cells together with cyp19a1a at the mRNA and/or protein levels. Sequence analysis identified one conserved Foxo binding site in the proximal promoter region of ricefield eel cyp19a1a. Transient transfection assay showed that Foxl2 may bind to the conserved Foxo site to activate cyp19a1a transcription and act synergistically with Nr5a1a. Mutation of either the conserved Nr5a1 site or Foxo site abolished or significantly decreased the synergistic effects of Nr5a1a and Foxl2 on cyp19a1a transcription. The sequence between Region III and I-box of Nr5a1a was critical to this synergistic effect. Dmrt1a modulated the Foxl2- and Nr5a1a-induced activation of cyp19a1a transcription and their synergistic effects in a biphasic manner, with inhibitory roles observed at lower doses (10 to 50ng) but release of the inhibition or even potentiating effects observed at higher doses (100 to 200ng). Collectively, data of the present study suggest that the interaction of Dmrt1a, Foxl2, and Nr5a1a in the ovarian follicular cells may facilitate the adequate expression of cyp19a1a and the production of estradiol, and contribute to the development and maturation of ovarian follicles in ricefield eels and other vertebrates as well.

Liposome-Encapsulated Rec8 and Dmrt1 Plasmids Induce Red-Spotted Grouper (Epinephelus akaara) Testis Maturation

In fish, the maturity of gonads plays an important role in the development and reproduction of the population, and it also dictates the success of captive breeding. Therefore, finding ways to promote gonadal maturation is an important goal in aquaculture. In this study, we injected recombinant dmrt1 and rec8 overexpression plasmids packaged in liposomes into the immature testis of red-spotted grouper (Epinephelus akaara) and measured the expression of Dmrt1 and Rec8 protein in vivo. Gonadosomatic index (GSI) and gonadal histology analyses showed that the testis developed from the immature to the mature state within 7 days after plasmid injection. Additionally, the spermatozoa concentration and motility in plasmid-injected fish was the same as that of naturally mature fish. These results provided evidence that delivery of dmrt1 and rec8 expression plasmids into the testis via injection induced testis maturation in vivo.

Noncoding Genes on Sex Chromosomes and Their Function in Sex Determination, Dosage Compensation, Male Traits, and Diseases

The mammalian Y chromosome has evolved in many species into a specialized chromosome that contributes to sex development among other male phenotypes. This function is well studied in terms of protein-coding genes. Less is known about the noncoding genome on the Y chromosome and its contribution to both sex development and other traits. Once considered junk genetic material, noncoding RNAs are now known to contribute to the regulation of gene expression and to play an important role in refining cellular functions. The prime examples are noncoding genes on the X chromosome, which mitigate the differential dosage of genes on sex chromosomes. Here, we discuss the evolution of noncoding RNAs on the Y chromosome and the emerging evidence of how micro, long, and circular noncoding RNAs transcribed from the Y chromosome contribute to sex differentiation. We briefly touch on emerging evidence that these noncoding RNAs also contribute to some other important clinical phenotypes in humans.

Involvement of Transforming Growth Factor Beta Family Genes in Gonadal Differentiation in Japanese Eel, Anguilla japonica, According to Sex-Related Gene Expressions

The gonochoristic feature with environmental sex determination that occurs during the yellow stage in the eel provides an interesting model to investigate the mechanisms of gonadal development. We previously studied various sex-related genes during gonadal sex differentiation in Japanese eels. In the present study, the members of transforming growth factor beta (TGF-β) superfamily were investigated. Transcript levels of anti-Müllerian hormone, its receptor, gonadal soma-derived factor (amh, amhr2, and gsdf, respectively) measured by real-time polymerase chain reaction (qPCR) showed a strong sexual dimorphism. Transcripts were dominantly expressed in the testis, and their levels significantly increased with testicular differentiation. In contrast, the expressions of amh, amhr2, and gsdf transcripts were low in the ovary of E2-feminized female eels. In situ hybridization detected gsdf (but not amh) transcript signals in undifferentiated gonads. amh and gsdf signals were localized to Sertoli cells and had increased significantly with testicular differentiation. Weak gsdf and no amh signals were detected in early ovaries of E2-feminized female eels. Transcript levels of amh and gsdf (not amhr2) decreased during human chorionic gonadotropin (HCG)-induced spermatogenesis in males. This study suggests that amh, amhr2, and especially gsdf might be involved in the gene pathway regulating testicular differentiation of Japanese eels.

Testis-specific expression pattern of dmrt1 and its putative regulatory region in the sea urchin (Mesocentrotus nudus)

Sea urchin (Mseocentrotus nudus) is an economically important mariculture species in several Asian countries. The growth rate and immunocompetence differ by sex in this species. However, the mechanisms of sex determination in M. nudus have remain unclear. In the present study, we focus on the dmrt1 gene of M. nudus (Mndmrt1) to investigate its dynamic expression pattern during different developmental stages. Real-time quantitative PCR (RT-qPCR) revealed that Mndmrt1 exhibits testis-specific expression and undetectable during the whole embryogenesis. With the development of ontogenetic, Mndmrt1 transcripts are first detected at 9 months post-fertilization (mpf). In addition, both the transcripts and protein of Mndmrt1 gene were specifically expressed in spermatogonia and spermatocytes, indicating that it might be a male germ cells marker in sea urchin. Significantly, the 1441 bp promoter sequence of Mndmrt1 gene was obtained by DNA walking, and one positive regulatory region at -1197/ -968 in the promoter, as well as one negative regulatory region at -1441/ -1198 have been identified by promoter activity analysis. Moreover, two regulatory regions contain multiple putative binding sites for transcription factors, including Sp1, Egr1, Sox5, CEBP, GATA and SRY. These findings suggest that Mndmrt1 may be related to testis differentiation and spermatogenesis in sea urchin and will provide an insight into understanding the regulatory mechanism of the dmrt1 gene.

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.

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.

Discovery of the Dmrt gene family members based on transcriptome analysis in mud crab Scylla paramamosain

Doublesex and mab-3 related transcription factors (Dmrts) play crucial roles in sex determination/differentiation and gonad development. The information on Dmrts and their functions are still scarce in mud crab Scylla paramamosain. In this study, 12 published transcriptome data of S. paramamosain were retrieved, pooled, and assembled. From the assembly, 7 Dmrt gene family members were identified and consisted of Spdmrt-like, Spdmrt-1a, Spdmrt-3, Spdmrt-11E, Spidmrt-1, Spdoublesex (Spdsx), and Spidmrt-2. These dmrt genes were predicted to encode 224 aa, 465 aa, 435 aa, 276 aa, 520 aa, 552 aa, and 266 aa protein precursors, respectively. The expression patterns of the dmrt genes were characterized by semi-quantitative PCR. The Spdmrt-like and Spdmrt-1a were exclusively detected in gonads, of which both expression levels in the testis were higher than that in the ovary. The Spdmrt-3, Spdmrt-11E, Spidmrt-1, Spdsx, and Spidmrt-2 were observed in various tissues; all these genes were sexually dimorphic except for dmrt-11E. Specifically, the expression level of Spdmrt-3 and Spidmrt-2 were higher in the testis than that in the ovary. On the contrary, the Spdsx and Spidmrt-1 expression level were higher in ovary than that in testis. The present study's findings provided a fundamental understanding of Dmrt gene family members involving sex determination/differentiation and gonad development in the S. paramamosain.

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.

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.

Full-length transcriptome sequencing and comparative transcriptomic analysis to uncover genes involved in early gametogenesis in the gonads of Amur sturgeon (Acipenser schrenckii)

BACKGROUND

Sturgeons (Acipenseriformes) are polyploid chondrostean fish that constitute an important model species for studying development and evolution in vertebrates. To better understand the mechanisms of reproduction regulation in sturgeon, this study combined PacBio isoform sequencing (Iso-Seq) with Illumina short-read RNA-seq methods to discover full-length genes involved in early gametogenesis of the Amur sturgeon, Acipenser schrenckii.

RESULTS

A total of 50.04 G subread bases were generated from two SMRT cells, and herein 164,618 nonredundant full-length transcripts (unigenes) were produced with an average length of 2782 bp from gonad tissues (three testes and four ovaries) from seven 3-year-old A. schrenckii individuals. The number of ovary-specific expressed unigenes was greater than those of testis (19,716 vs. 3028), and completely different KEGG pathways were significantly enriched between the ovary-biased and testis-biased DEUs. Importantly, 60 early gametogenesis-related genes (involving 755 unigenes) were successfully identified, and exactly 50% (30/60) genes of those showed significantly differential expression in testes and ovaries. Among these, the Amh and Gsdf with testis-biased expression, and the Foxl2 and Cyp19a with ovary-biased expression strongly suggested the important regulatory roles in spermatogenesis and oogenesis of A. schrenckii, respectively. We also found the four novel Sox9 transcript variants, which increase the numbers of regulatory genes and imply function complexity in early gametogenesis. Finally, a total of 236,672 AS events (involving 36,522 unigenes) were detected, and 10,556 putative long noncoding RNAs (lncRNAs) and 4339 predicted transcript factors (TFs) were also respectively identified, which were all significantly associated with the early gametogenesis of A. schrenckii.

CONCLUSIONS

Overall, our results provide new genetic resources of full-length transcription data and information as a genomic-level reference for sturgeon. Crucially, we explored the comprehensive genetic characteristics that differ between the testes and ovaries of A. schrenckii in the early gametogenesis stage, which could provide candidate genes and theoretical basis for further the mechanisms of reproduction regulation of sturgeon.