Comparative transcriptome analysis reveals the sexual dimorphic expression profiles of mRNAs and non-coding RNAs in the Asian yellow pond turtle (Meauremys mutica)

Screening of temperature-responsive signalling molecules during sex differentiation in Asian yellow pond turtle (Mauremys mutica)

BACKGROUND

The Asian yellow pond turtle (Mauremys mutica) is an important commercial freshwater aquaculture species in China. This species is a highly sexually dimorphic species, with males growing at a faster rate than females and exhibits temperature-dependent sex determination (TSD), in which the incubation temperature during embryonic development determines the sexual fate. However, the mechanisms of the sex determination or sex differentiation in the Asian yellow pond turtle are remain a mystery.

RESULTS

Temperature-specific gonadal transcriptomics of the Asian yellow pond turtle were performed during the thermosensitive period (stage 15) using RNA-seq technology to identify candidate genes that initiate gonadal differentiation. We uncovered candidates that were the first to respond to temperature. These candidates were sexually dimorphic in expression, reflecting differences in gonadal (Cirbp, Runx1) and germline differentiation (Vasa, Nanos1, Piwil2), gametogenesis (Hmgb3, Zar1, Ovoinhibitor-like, Kif4), steroid hormone biosynthesis (Hsd17b5, Hsd17b6), heat shock (Dnajb6, Hsp90b1, Hsp90aa1) and transient receptor potential channel genes (Trpm1, Trpm4, Trpm6, Trpv1).

CONCLUSIONS

Our work will provide important genetic information to elucidate the mechanisms of sex control in the Asian yellow pond turtles, and will contribute important genetic resources for further studies of temperature-dependent sex determination in turtles.

Transcriptome analysis of differentially expressed circular RNAs in the testis and ovary of golden pompano (Trachinotus blochii)

The artificial breeding of golden pompano (Trachinotus blochii) has expanded greatly in recent years, and after long-term breeding efforts, clear sexual dimorphisms have been observed in T. blochii growth traits, with females growing faster. As sponges of microRNA (miRNAs), circular RNAs (CircRNAs) can alleviate miRNA inhibition of target mRNA. However, few studies have examined sex-related CircRNAs and none of those have looked at T. blochii. To further understand the role of CircRNAs in sex differentiation and sexual size dimorphism in T. blochii, six CircRNA libraries were constructed from the testes and ovaries of T. blochii. A total of 1522 CircRNAs were found distributed over all 24 chromosomes of T. blochii. 135 differentially expressed CircRNAs (DECs) were identified by screening, These DECs were then subjected to GO enrichment, which found 47 enriched pathways. A number of CircRNAs were enriched in cellular processes and metabolic processes. According to the KEGG pathway analysis, a series of sex differentiation pathways were enriched, including the GnRH, calcium, and MAPK signaling pathways. Furthermore, we selected two CircRNAs from the DECs named circ-cacna1b and circ-octc. We found that the cacna1b gene is regulated by 7 miRNAs, 3 of which were regulated by circ-cacna1b, i.e., mmu-miR-138-5p, fru-miR-138, and pma-miR-138b. In addition, the miRNA named pma-miR-138b can regulate sex-related genes, such as sox9 and dmrt1, among others. The co-expression network of CircRNA-miRNA-mRNA showed circ-cacna1b may play a crucial role in T. blochii sex differentiation by regulating pma-miR-138b to affect the expression of sex differentiation genes. The circ-octc may be one of the largest contributors to sexual size dimorphism during growth through its effect on lipid metabolism. These findings could broaden our understanding of CircRNAs and provide new insight into their function in sex differentiation and growth.

Temporal variation in DNA methylation during gonadal development in a reptile with temperature-dependent sex determination

DNA methylation plays a significant role in transducing external environmental signals to a cellular response in reptiles; however, whether the methylation patterns are conserved across species remains unclear. Here, we examined the genome-wide DNA methylation differentiation between male and female hatchling gonads of the temperature-dependent sex determination (TSD) Mauremys mutica (M. mutica) using methylation-dependent restriction-site associated DNA sequencing (MethylRAD-seq) to test differentially methylated genes underlying sexual development. Several categories, including heat shock genes (HSP90A, HSP30C), histone- (KDM8) and ubiquitin-related genes (TRIM39), kinases (WNK3) and gonad differentiation or gonadal development related genes (HSD17B8, HSD17B12), were identified as candidates for future study. Additionally, we identified several regulatory pathways potentially mediating TSD thermosensitivity such as the GnRH signaling pathway and calcium signaling pathway. These findings provide evidence that sexually dimorphic DNA methylation may be associated with sex determination or sex differentiation in TSD M. mutica.

DNA Methylation Difference between Female and Male Ussuri Catfish (Pseudobagrus ussuriensis) in Brain and Gonad Tissues

DNA methylation has been found to be involved in sex determination and differentiation in many aquaculture species. The Ussuri catfish (Pseudobagrus ussuriensis) is a popular aquaculture fish in China with high economic value in which male-biased sex dimorphism was observed in terms of body size and body weight. In this study, DNA methylation-sensitive RAD sequencing (Methyl-RAD) was used to explore the epigenetic difference between adult male and female samples in brain and gonad tissues. In brain tissues, 5,442,496 methylated cytosine sites were found and 9.94% of these sites were from symmetric CCGG or CCWGG sites. Among these sites, 321 differential DNA methylation sites (DMSs) in 171 genes were identified, while in gonad tissues, 4,043,053 methylated cytosines sites were found in total and 11.70% of them were from CCGG or CCWGG. Among these sites, 78 differential DNA methylation sites were found which were located in 64 genes. We also found several sex-determination genes among these differential methylated genes, such as amh, gsdf and hsd11b2 in brain tissues and slco3a1, socs2 and trim47 in gonad tissues. These results provided evidence for understanding the function of DNA methylation in the sex differentiation in Pseudobagrus ussuriensis, which further deepens the relationship between gene regulation and epigenetics.

Transcriptome analysis reveals key genes and pathways related to sex differentiation in the Chinese soft-shelled turtle (Pelodiscus sinensis)

Most vertebrates exhibit sexual dimorphisms in size, colour, behaviour, physiology and many others. The Chinese soft-shelled turtle (Pelodiscus sinensis) male individuals reach a larger size than females which produce significant economic implications in aquaculture. However, the mechanisms of sex determination and plastic patterns of sex differentiation in P. sinensis remain unclear. Here, comparative transcriptome analysis on male and female embryonic gonads prior to gonad formation and stages mediated gonadal differentiation of P. sinensis were performed to characterize the potential sex-related genes and their molecular pathways in P. sinensis. A total of 6369 differentially expressed genes (DEGs) were identified from day 9 and day 16 and assigned to 626 GO pathways and 161 KEGG signalling pathways, including ovarian steroidogenesis pathway, steroid hormone biosynthesis pathways, and the GnRH signalling pathway (P < 0.05). Moreover, protein interaction network analyses revealed that Akr1c3, Sult2b1, Sts, Cyp3a, Cyp1b1, Sox30 and Lhx9 might be key candidate genes for sex differentiation in P. sinensis. These data provide a genomic rationale for the sex differentiation of P. sinensis and enrich the candidate gene pool for sex differentiation.

Whole-Transcriptome Analysis Identifies Gender Dimorphic Expressions of Mrnas and Non-Coding Rnas in Chinese Soft-Shell Turtle (Pelodiscus sinensis)

In aquaculture, the Chinese soft-shelled turtle (Pelodiscus sinensis) is an economically important species with remarkable gender dimorphism in its growth patterns. However, the underlying molecular mechanisms of this phenomenon have not been elucidated well. Here, we conducted a whole-transcriptome analysis of the female and male gonads of P. sinensis. Overall, 7833 DE mRNAs, 619 DE lncRNAs, 231 DE circRNAs, and 520 DE miRNAs were identified. Some "star genes" associated with sex differentiation containing dmrt1, sox9, and foxl2 were identified. Additionally, some potential genes linked to sex differentiation, such as bmp2, ran, and sox3, were also isolated in P. sinensis. Functional analysis showed that the DE miRNAs and DE ncRNAs were enriched in the pathways related to sex differentiation, including ovarian steroidogenesis, the hippo signaling pathway, and the calcium signaling pathway. Remarkably, a lncRNA/circRNA-miRNA-mRNA interaction network was constructed, containing the key genes associated with sex differentiation, including fgf9, foxl3, and dmrta2. Collectively, we constructed a gender dimorphism profile of the female and male gonads of P. sinensis, profoundly contributing to the exploration of the major genes and potential ncRNAs involved in the sex differentiation of P. sinensis. More importantly, we highlighted the potential functions of ncRNAs for gene regulation during sex differentiation in P. sinensis as well as in other turtles.

Chromosome-level genome assembly of Asian yellow pond turtle (Mauremys mutica) with temperature-dependent sex determination system

Knowledge of sex determination has important implications in physiology, ecology and genetics, but the evolutionary mechanisms of sex determination systems in turtles have not been fully elucidated, due to a lack of reference genomes. Here, we generate a high-quality genome assembly of Asian yellow pond turtle (Mauremys mutica) using continuous long-read (PacBio platform), Illumina, and high-throughput chromatin conformation capture (Hi-C) technologies. The M. mutica haplotype has a genome size of 2.23 Gb with a contig N50 of 8.53 Mb and scaffold N50 of 141.98 Mb. 99.98% sequences of the total assembly are anchored to 26 pseudochromosomes. Comparative genomics analysis indicated that the lizard-snake-tuatara clade diverged from the bird-crocodilian-turtle clade at approximately 267.0-312.3 Mya. Intriguingly, positive selected genes are mostly enriched in the calcium signaling pathway and neuroactive ligand-receptor interaction, which are involved in the process of temperature-dependent sex determination. These findings provide important evolutionary insights into temperature-dependent sex determination system.

Comparative transcriptomic analysis reveals the gonadal development-related gene response to environmental temperature in Mauremys mutica

The Asian yellow pond turtle (Mauremys mutica) displays temperature-dependent sex determination (TSD), in which incubation temperature during embryonic development determines the sexual fate of the individual. However, the mechanism of the sex determination/differentiation of Mauremys mutica remains a mystery. Here, we first analyzed the temperature-specific gonadal transcriptomes of Mauremys mutica prior to gonad formation and gonads during the thermosensitive period. We uncovered a list of candidates that respond to temperature stimuli enriched in several categories, such as heat shock protein family members dnajb6a, dnaja4, hspa8 and hsp90aa1, temperature sensor genes mmp17 and mmp28, and putative novel temperature-responsive genes tmco6, gria3 and eif3f. Notably, striking differences were identified in the expression profiles of genes underlying sexual development, such as tex15, insr, igf1r, cirbp, esr1, dmrt2 and Serpinh1. Moreover, we analyzed the similarity and divergence of the timecourse of gene expression among Mauremys mutica and two other reported TSD turtles (Trachemys scripta and Chrysemys picta). The shared genes revealed the common gonad-specific regulatory mechanisms existing in these three TSD turtles that initiate their sexual development. Therefore, our findings could provide basic data to elucidate the mechanisms of sex determination/differentiation of M. mutica, even contributing to further understanding of these mechanisms in other TSD turtles.

MicroRNAs May Play an Important Role in Sexual Reversal Process of Chinese Soft-Shelled Turtle, Pelodiscus sinensis

The Chinese soft-shelled (Pelodiscus sinensis) turtle exhibits obvious sex dimorphism, which leads to the higher economic and nutritional value of male individuals. Exogenous hormones can cause the transformation from male to female phenotype during gonadal differentiation. However, the molecular mechanism related to the sexual reversal process is unclear. In this study, we compared the difference between the small RNAs of male, female, and pseudo-female turtles by small RNA-seq to understand the sexual reversal process of Chinese soft-shelled turtles. A certain dose of estrogen can cause the transformation of Chinese soft-shelled turtles from male to female, which are called pseudo-female individuals. The result of small RNA-seq has revealed that the characteristics of pseudo-females are very similar to females, but are strikingly different from males. The number of the microRNAs (miRNAs) of male individuals was significantly less than the number of female individuals or pseudo-female individuals, while the expression level of miRNAs of male individuals were significantly higher than the other two types. Furthermore, we found 533 differentially expressed miRNAs, including 173 up-regulated miRNAs and 360 down-regulated miRNAs, in the process of transformation from male to female phenotype. Cluster analysis of the total 602 differential miRNAs among females, males, and pseudo-females showed that miRNAs played a crucial role during the sexual differentiation. Among these differential miRNAs, we found 12 miRNAs related to gonadal development and verified their expression by qPCR. The TR-qPCR results confirmed the differential expression of 6 of the 12 miRNAs: miR-26a-5p, miR-212-5p, miR-202-5p, miR-301a, miR-181b-3p and miR-96-5p were involved in sexual reversal process, which was consistent with the results of omics. Using these six miRNAs and some of their target genes, we constructed a network diagram related to gonadal development. We suggest that these miRNAs may play an important role in the process of effective sex reversal, which would contribute to the breeding of all male strains of Chinese soft-shelled turtles.

Dynamic expression and functional analysis of circular RNA in the gonads of Chinese soft-shelled turtles (Pelodiscus sinensis)

Circular RNA (circRNA) is a noncoding RNA that can regulate a variety of biological processes. CircRNAs can regulate gene expression posttranscriptionally by acting as microRNA sponges. Many turtle species are remarkable organisms due to their reproductive processes. However, information on circRNA in the gonads of turtles is limited. In this study, 6, 121 circRNAs were identified in the testes and ovaries of Chinese soft-shelled turtles (Pelodiscus sinensis) using the Illumina platform, and 710 circRNAs were significantly differentially expressed (DE). The DE circRNAs included 541 upregulated and 169 downregulated circRNAs in the testes. GO and KEGG pathway analysis indicated that the DE circRNAs were enriched in several signaling pathways, including GnRH, Wnt, FoxO, Progesterone mediated oocyte maturation, and mTOR signaling pathways. Five DE circRNAs were randomly selected, and their relative expression levels in ovaries and testes were detected by quantitative real-time PCR. All of these circRNAs were differentially expressed. In addition, 9, 883 interactions between circRNAs and miRNAs were predicted in the turtles. Target genes of the miRNAs include a range of genes regulating gonadal development. Seven ceRNA networks (DE circRNAs-DE miRNAs-DE mRNAs), including 7 DE circRNAs, 11 DE miRNAs and 20 DE mRNAs, were constructed. The networks included Cdc6, the miR-1 family, the miR-203 family, and the miR-302 family. The expression profile of gonadal circRNAs might help to elucidate the roles of nonprotein coding RNAs in turtle gonadal development.

Vasa expression is associated with sex differentiation in the Asian yellow pond turtle, Mauremys mutica

Vasa, one of the best-studied germ cell markers plays a critical role in germ cell development and differentiation in animals. Vasa deficiency would lead to male-specific sterility in most vertebrates, but female sterility in the fly. However, the role of the vasa gene involved in germ cell differentiation is largely elusive. Here, we first characterized the expression profile of vasa products in the Asian yellow pond turtle by quantitative reverse-transcription polymerase chain reaction and fluorescence immunostaining. The results showed that vasa messenger RNA (mRNA) is initially detected in embryos at stage 16, and then dramatically increased in embryos at stage 19. In particular, like the sex-related genes, vasa mRNA exhibited differential expression in embryos between the male-producing temperature (MPT, 25°C) and the female-producing temperature (FPT, 33°C), whereas there was no difference in methylation levels of vasa promoter detected between FPT and MPT. In contrast, in the adult Asian yellow pond, the level of vasa mRNA was much higher in the testis than ovary. Moreover, the immunostaining on testicular sections and cells showed that Vasa protein was exclusively expressed in germ cells: Weak but detectable in spermatogonia, highest in spermatocytes, moderate and concentrated in chromatid bodies in spermatids and spermatozoa, and bare in somatic cells. The expression profile of Vasa protein is similar in turtle species studied so far but distinct from those in fish species in this study. The findings of this study would provide new insights into our understanding of the conservation and divergence of the vasa gene, even other germ cell genes across phyla.