Single-Molecule Real-Time Sequencing for Identifying Sexual-Dimorphism-Related Transcriptomes and Genes in the Chinese Soft-Shelled Turtle (Pelodiscus sinensis)

Abstract

The Chinese soft-shelled turtle (Pelodiscus sinensis), an economically important aquatic species in China, displays considerable sexual dimorphism: the male P. sinensis is larger and, thus, more popular in the market. In this study, we obtained the full-length (FL) transcriptome data of P. sinensis by using Pacific Biosciences (PacBio)’s isoform sequencing and analyzed the transcriptome structure. In total, 1,536,849 high-quality FL transcripts were obtained through single-molecule real-time (SMRT) sequencing, which were then corrected using Illumina sequencing data. Next, 89,666 nonredundant FL transcripts were generated after mapping to the reference genome of P. sinensis; 291 fusion genes and 17,366 novel isoforms were successfully annotated using data from the nonredundant protein sequence database (NR), eukaryotic orthology groups (KOG), the Gene Ontology (GO) project, and the KEGG Orthology (KO) database. Additionally, 19,324 alternative polyadenylation sites, 101,625 alternative splicing events, 12,392 long noncoding RNAs, and 5916 transcription factors were identified. Smad4, Wif1, and 17-β-hsd were identified as female-biased genes, while Nkd2 and Prp18 held a higher expression level in males than females. In summary, we found differences between male and female P. sinensis individuals in AS, lncRNA, genes, and transcripts, which relate to the Wnt pathway, oocyte meiosis, and the TGF-β pathway. Female-biased genes such as Smad4, Wif1, and 17-β-hsd and male-biased genes such as Nkd2 and Prp18 played important roles in the sex determination of P. sinensis. FL transcripts are a precious resource for characterizing the transcriptome of P. sinensis, laying the foundation for further research on the sex-determination mechanisms of P. sinensis.