Brain and Pituitary Transcriptome Analyses Reveal the Differential Regulation of Reproduction-Related LncRNAs and mRNAs in Cynoglossus semilaevis

Long noncoding RNAs (lncRNAs) have been identified to be involved in half-smooth tongue sole (Cynoglossus semilaevis) reproduction. However, studies of their roles in reproduction have focused mainly on the ovary, and their expression patterns and potential roles in the brain and pituitary are unclear. Thus, to explore the mRNAs and lncRNAs that are closely associated with reproduction in the brain and pituitary, we collected tongue sole brain and pituitary tissues at three stages for RNA sequencing (RNA-seq), the 5,135 and 5,630 differentially expressed (DE) mRNAs and 378 and 532 DE lncRNAs were identified in the brain and pituitary, respectively. The RNA-seq results were verified by RT-qPCR. Moreover, enrichment analyses were performed to analyze the functions of DE mRNAs and lncRNAs. Interestingly, their involvement in pathways related to metabolism, signal transduction and endocrine signaling was revealed. LncRNA-target gene interaction networks were constructed based on antisense, cis and trans regulatory mechanisms. Moreover, we constructed competing endogenous RNA (ceRNA) networks. In summary, this study provides mRNA and lncRNA expression profiles in the brain and pituitary to understand the molecular mechanisms regulating tongue sole reproduction.

Integration of Transcriptome and Methylome Highlights the Roles of Cell Cycle and Hippo Signaling Pathway in Flatfish Sexual Size Dimorphism

Sexual size dimorphism (SSD) is the difference in segments or body size between sexes prevalent in various species. Understanding the genetic architecture of SSD has remained a significant challenge owing to the complexity of growth mechanisms and the sexual influences among species. The Chinese tongue sole (Cynoglossus semilaevis), which exhibits a female-biased SSD and sex reversal from female to pseudomale, is an ideal model for exploring SSD mechanism at the molecular level. The present study aimed to integrate transcriptome and methylome analysis to unravel the genetic and epigenetic changes in female, male, and pseudomale C. semilaevis. The somatotropic and reproductive tissues (brain, liver, gonad, and muscle) transcriptomes were characterized by RNA-seq technology. Transcriptomic analysis unravelled numerous differentially expressed genes (DEGs) involved in cell growth and death-related pathways. The gonad and muscle methylomes were further employed for screening differentially methylated genes (DMGs). Relatively higher DNA methylation levels were observed in the male and pseudomale individuals. In detail, hypermethylation of the chromosome W was pronounced in the pseudomale group than in the female group. Furthermore, weighted gene co-expression network analysis showed that turquoise and brown modules positively and negatively correlated with the female-biased SSD, respectively. A combined analysis of the module genes and DMGs revealed the female-biased mRNA transcripts and hypomethylated levels in the upstream and downstream regions across the cell cycle-related genes. Moreover, the male and pseudomale-biased gene expression in the hippo signaling pathway were positively correlated with their hypermethylation levels in the gene body. These findings implied that the activation of the cell cycle and the inhibition of the hippo signaling pathway were implicated in C. semilaevis female-biased SSD. In addition, the dynamic expression pattern of the epigenetic regulatory factors, including dnmt1, dnmt3a, dnmt3b, and uhrf1, among the different sexes correspond with their distinct DNA methylation levels. Herein, we provide valuable clues for understanding female-biased SSD in C. semilaevis.

A Z-Linked E3 Ubiquitin Ligase Cs-rchy1 Is Involved in Gametogenesis in Chinese Tongue Sole, Cynoglossus semilaevis

Ubiquitin ligase (E3) plays a versatile role in gonadal development and spermatogenesis in mammals, while its function in fish is little reported. In this study, a Z-chromosome linked ubiquitin ligase rchy1 in C. semilaevis (Cs-rchy1) was cloned and characterized. The full-length cDNA was composed of 1962 bp, including 551 bp 5′UTR, 736 bp 3′UTR, and 675 bp ORF encoding a 224-amino-acid (aa) protein. Cs-rchy1 was examined among seven different tissues and found to be predominantly expressed in gonads. In testis, Cs-rchy1 could be detected from 40 days post hatching (dph) until 3 years post hatching (yph), but there was a significant increase at 6 months post hatching (mph). In comparison, the expression levels in ovary were rather stable among different developmental stages. In situ hybridization showed that Cs-rchy1 was mainly localized in germ cells, that is, spermatid and spermatozoa in testis and stage I, II and III oocytes in ovary. In vitro RNA interference found that Cs-rchy1 knockdown resulted in the decline of sox9 and igf1 in ovarian cell line and down-regulation of cyp19a in the testicular cell line. These data suggested that Cs-rchy1 might participate in gonadal differentiation and gametogenesis, via regulating steroid hormone synthesis.