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.

Introduction: Bone morphogenetic proteins (BMPs) play a crucial role in bone formation and differentiation. Recent RNA-Seq results suggest that BMPs may be involved in the sex differentiation of P. sinensis, yet more relevant studies about BMPs in P. sinensis are lacking.Methods: Herein, we identified BMP gene family members, analyzed the phylogeny, collinear relationship, scaffold localization, gene structures, protein structures, transcription factors and dimorphic expression by using bioinformatic methods based on genomic and transcriptomic data of P. sinensis. Meanwhile, qRT-PCR was used to verify the RNA-Seq results and initially explore the function of the BMPs in the sex differentiation of P. sinensis.Results: A total of 11 BMP genes were identified, 10 of which were localized to their respective genomic scaffolds. Phylogenetic analysis revealed that BMP genes were divided into eight subfamilies and shared similar motifs (“WII”, “FPL”, “TNHA”, “CCVP”, and “CGC”) and domain (TGF-β superfamily). The results of the sexually dimorphic expression profile and qRT-PCR showed that Bmp2, Bmp3, Bmp15l, Bmp5, Bmp6 and Bmp8a were significantly upregulated in ovaries, while Bmp2lb, Bmp7, Bmp2bl and Bmp10 were remarkable upregulated in testes, suggesting that these genes may play a role in sex differentiation of P. sinensis.Discussion: Collectively, our comprehensive results enrich the basic date for studying the evolution and functions of BMP genes in P. sinensis.

Male and female Chinese soft-shelled turtles (Pelodiscus sinensis) have sex-dimorphic growth patterns, and males have higher commercial value because of their larger size and thicker calipash. Thus, developing sex-specific markers is beneficial to studies on all-male breeding in P. sinensis. Here, we developed an accurate and efficient workflow for the screening of sex-specific sequences with ZW or XY sex determination systems. Based on this workflow, female and male P. sinensis reference genomes of 2.23 Gb and 2.26 Gb were obtained using de novo assembly. After aligning and filtering, 4.01 Mb female-specific sequences were finally identified. Subsequently, the seven developed sex-specific primer pairs were 100% accurate in preliminary, population, and embryonic validation. The presence and absence of bands for the primers of P44, P45, P66, P67, P68, and P69, as well as two and one bands for the PB1 primer, indicate that the embryos are genetically female and male, respectively. NR and functional annotations identified several sex-determining candidate genes and related pathways, including Ran, Eif4et, and Crkl genes, and the insulin signaling pathway and the cAMP signaling pathway, respectively. Collectively, our results reveal that a ZW-type sex-determination system is present in P. sinensis and provide novel insights for the screening of sex-specific markers, sex-control breeding, and the studies of the sex determination mechanism of P. sinensis.

The Chinese soft-shelled turtle (Pelodiscus sinensis), as a type of warm-water reptile, could be induced to massive death by sharp temperature decline. Hence, the mechanism of spleen tissue responding to cold stress in the P. sinensis was investigated. The present results showed that the superoxide dismutase (SOD) activity declined from 4 to 16 days post-cold-stress (dps), while the catalase (CAT) and glutathione peroxidase (GSH-Px) activities increased, from 4 to 8 dps in the 14 °C (T14) and 7 °C (T7) stress groups. The spleen transcriptome in the T7 group and the control group (CG) at 4 dps obtained 2625 differentially expressed genes (DEGs), including 1462 upregulated and 1663 downregulated genes. The DEGs were enriched mainly in the pathways “intestinal immune network for IgA production” (Pigr, Il15ra, Tnfrsf17, Aicda, and Cd28), “toll-like receptor signaling pathway” (Mapk10, Tlr2, Tlr5, Tlr7, and Tlr8), and “cytokine–cytokine receptor interaction” (Cx3cl1, Cx3cr1, Cxcl14, Cxcr3, and Cxcr4). The metabolomic data showed that esculentic acid, tyrosol, diosgenin, heptadecanoic acid, and 7-ketodeoxycholic acid were obviously increased, while baccatin III, taurohyocholate, parthenolide, enterolactone, and tricin were decreased, in the CG vs. T7 comparison. Integrated analysis of the two omics revealed that “glycine, serine and threonine metabolism”, “FoxO signaling pathway”, and “neuroactive ligand–receptor interaction” were the main pathways responding to the cold stress. Overall, this work found that low temperature remarkably influenced the antioxidant enzyme activities, gene expression pattern, and metabolite profile in the spleen, indicating that immunity might be weakened by cold stress in P. sinensis.