Sexually Dimorphic Gene Expression Associated with Growth and Reproduction of Tongue Sole (Cynoglossus semilaevis) Revealed by Brain Transcriptome Analysis

Unveiling Gene Expression Dynamics during Early Embryogenesis in Cynoglossus semilaevis: A Transcriptomic Perspective

Commencing with sperm-egg fusion, the early stages of metazoan development include the cleavage and formation of blastula and gastrula. These early embryonic events play a crucial role in ontogeny and are accompanied by a dramatic remodeling of the gene network, particularly encompassing the maternal-to-zygotic transition. Nonetheless, the gene expression dynamics governing early embryogenesis remain unclear in most metazoan lineages. We conducted transcriptomic profiling on two types of gametes (oocytes and sperms) and early embryos (ranging from the four-cell to the gastrula stage) of an economically valuable flatfish-the Chinese tongue sole Cynoglossus semilaevis (Pleuronectiformes: Cynoglossidae). Comparative transcriptome analysis revealed that large-scale zygotic genome activation (ZGA) occurs in the blastula stage, aligning with previous findings in zebrafish. Through the comparison of the most abundant transcripts identified in each sample and the functional analysis of co-expression modules, we unveiled distinct functional enrichments across different gametes/developmental stages: actin- and immune-related functions in sperms; mitosis, transcription inhibition, and mitochondrial function in oocytes and in pre-ZGA embryos (four- to 1000-cell stage); and organ development in post-ZGA embryos (blastula and gastrula). These results provide insights into the intricate transcriptional regulation of early embryonic development in Cynoglossidae fish and expand our knowledge of developmental constraints in vertebrates.

Comprehensive Transcriptome Analysis Reveals Sex-Specific Alternative Splicing Events in Zebrafish Gonads

Alternative splicing is an important way of regulating gene functions in eukaryotes. Several key genes involved in sex determination and gonadal differentiation, such as nr5a1 and ddx4, have sex-biased transcripts between males and females, suggesting a potential regulatory role of alternative splicing in gonads. Currently, the sex-specific alternative splicing events and genes have not been comprehensively studied at the genome-wide level in zebrafish. In this study, through global splicing analysis on three independent sets of RNA-seq data from matched zebrafish testes and ovaries, we identified 120 differentially spliced genes shared by the three datasets, most of which haven’t been reported before. Functional enrichment analysis showed that the GO terms of mRNA processing, mRNA metabolism and microtubule-based process were strongly enriched. The testis- and ovary-biased alternative splicing genes were identified, and part of them (tp53bp1, tpx2, mapre1a, kif2c, and ncoa5) were further validated by RT-PCR. Sequence characteristics analysis suggested that the lengths, GC contents, and splice site strengths of the alternative exons or introns may have different influences in different types of alternative splicing events. Interestingly, we identified an unexpected high proportion (over 70%) of non-frameshift exon-skipping events, suggesting that in these cases the two protein isoforms derived from alternative splicing may both have functions. Furthermore, as a representative example, we found that the alternative splicing of ncoa5 causes the loss of a conserved RRM domain in the short transcript predominantly produced in testes. Our study discovers novel sex-specific alternative splicing events and genes with high reliabilities in zebrafish testes and ovaries, which would provide attractive targets for follow-up studies to reveal the biological significances of alternative splicing events and genes in sex determination and gonadal differentiation.

Identification of crucial factors involved in Cynoglossus semilaevis sexual size dimorphism by GWAS and demonstration of zbed1 regulatory network by DAP-seq

Sexual size dimorphism (SSD), whereby females and males exhibit different body sizes, are widely documented in animals. To explore crucial regulators implicated in female-biased SSD of Chinese tongue sole (Cynoglossus semilaevis), GWAS was conducted on 350 females and 59 males. Twenty SNPs and 25 genes including zbed1, nsd3, cdc45, klhl29, and smad4 with -log(p) > 7 were screened, mainly mapping to sex chromosome. The chromosome W-linked gene zbed1 attracted particular attention because it is a master key for cell proliferation. Thus, the regulatory network of zbed1 in C. semilaevis was explored by DAP-seq and 1352 peaks were discovered in the female brain. Moreover, zbed1 potentially regulated hippo signaling pathway, cell cycle, translation, and PI3k-Akt signaling pathway in C. semilaevis. These findings identify crucial SNPs and genes associated with female-biased SSD in C. semilaevis, also provide the first genome-wide survey for the zbed1 regulatory network in fish species.

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.

Pomc Plays an Important Role in Sexual Size Dimorphism in Tilapia

Sexual dimorphism is common across the animal kingdom. Knowledge of the mechanisms of sexual size dimorphism is limited although it is important in biology and aquaculture. Tilapia is the common name for ~ 100 species of cichlid fish. Some are important aquaculture species and males outgrow females. To gain novel insights into the mechanisms underlying sexual size dimorphism, we analyzed the differences of brain transcriptomes between males and females in Mozambique tilapia and studied the function of the pro-opiomelanocortin (Pomc) gene in tilapia and zebrafish. The transcriptome analysis identified 123, 55, and 2706 sex-biased genes at 5, 30, and 90 dph (days post-hatch), respectively, indicating sexual dimorphism of gene expressions in the brain. The expression of Pomc in the tilapia brain was a female-biased at 30, 90, and 120 dph. An analysis of the DNA sequence located upstream of the tilapia Pomc transcriptional start site identified two estrogenic response elements. In vitro luciferase assay of the two elements revealed that β-estradiol significantly enhanced the expression of luciferase activity, suggesting that the expression of Pomc is mediated by estrogen. We knocked out Pomc in zebrafish using Crispr/Cas-9. The Pomc-knockout zebrafish showed faster growth and higher sensitivity to feeding as compared to the wild-type fish. Taken together, our results indicate that Pomc contributes to sexual size dimorphism and suggest that the high estrogen level in females promotes the expression of Pomc and suppresses feeding in female tilapias, which leads to the slower growth of female tilapias.

Transcriptome Profiling Reveals the Sexual Dimorphism of Gene Expression Patterns during Gonad Differentiation in the Half-Smooth Tongue Sole (Cynoglossus semilaevis)

The half-smooth tongue sole (Cynoglossus semilaevis), one of the most economically-important fish species in China, exhibits sexually dimorphic growth. An understanding of sex-related gene expression patterns in the tongue sole may inform sex regulation and breeding processes that increase fish production. However, the gene expression patterns during gonad development in the tongue sole remain unknown. In this study, transcriptome sequencing analyses were performed during gonad differentiation in the tongue sole, namely, at 62 days post-hatching (dph), 100 dph, 120 dph, and 150 dph. Differentially expressed genes associated with sex differentiation and gonad development were identified at each time point. Trend analysis showed that gene expression patterns varied over time. These expression patterns either explained common, non-sexually-dimorphic features or indicated significant sexual dimorphism. Transcript structure analyses identified both sex and time differences among samples. This study investigated the time-dependent expression patterns of several sex-related genes, including Dmrt1, Amh, Foxl2, aromatase encoding gene, Esr, and the Sox gene family, during gonad differentiation in the tongue sole. These results might clarify the significant sexual differences during early development in the tongue sole and might provide insight into the mechanisms controlling sex differentiation and development.

Transcriptomic analysis of female and male gonads in juvenile snakeskin gourami (Trichopodus pectoralis)

The snakeskin gourami (Trichopodus pectoralis) exhibits sexual dimorphism, particularly in body size. Since the snakeskin gourami is usually marketed during sexual maturation, the sexual size dimorphism has become an economically important trait. Sex-biased gene expression plays a key role in phenotypic sexual dimorphism. Therefore, using high-throughput RNA sequencing (RNA-seq) technology, we aimed to explore the differentially expressed genes (DEGs) in ovary and testis during sex differentiation in juvenile snakeskin gourami. Our results revealed a number of DEGs were demonstrated to be overexpressed in ovary (11,625 unigenes) and testis (16,120 unigenes), and the top 10 female-biased (rdh7, dnajc25, ap1s3, zp4, polb, parp12, trim39, gucy2g, rtbs, and fdxr) and male-biased (vamp3, nbl1, dnah2, ccdc11, nr2e3, spats1, pih1d2, tekt3, fbxo36, and mybl2) DEGs were suggested to be mainly associated with ovary and testis differentiation, respectively. Additionally, using real-time reverse transcription polymerase chain reaction (qRT-PCR), validation of the differential expression of 21 genes that were previously shown to be related to gonad development was performed (ar, bHLH, cyp19a1, daz, dead-end, esrb, esrrg, gnrhr, gpa, gsg1l, hsd17B, mospd1, nanos-1, nanos-2, p53, piwi-1, piwi-2, rerg, rps6ka, tgf-beta, and VgR). The results showed a significantly positive correlation (0.84; P < 0.001) between the results of RNA-seq and qRT-PCR. Therefore, RNA-seq analysis in our study identified global genes that were associated with ovary and testis differentiation in the juvenile phase of the snakeskin gourami. Our findings provide valuable transcriptomic bioinformation for further investigation of reproductive biology and applications of sex manipulation.

Integrative analysis reveals pathways associated with sex reversal in Cynoglossus semilaevis

Sex reversal is a complex biological phenomenon exhibited by Cynoglossus semilaevis. Some genetic females may irreversibly convert to pseudomales, thus increasing aquaculture costs because males grow much more slowly than females. In this study, an integrative analysis of transcriptome and proteome was performed to compare differences in gene and protein expression in females and pseudomales after gonad differentiation in C. semilaevis. Based on RNA-Seq results, 1893 genes showed differences in expression at the transcript level between females and pseudomales. Of these differentially expressed genes (DEGs), zona pellucida sperm-binding protein 4-like (LOC103393374 , ZP4), zona pellucida sperm-binding protein 4-like (LOC103396071, ZP4) and forkhead box L2 (foxl2) were highly expressed in females and doublesex and mab-3 related transcription factor 1(dmrt1) and doublesex and mab-3 related transcription factor 3 (dmrt3) were highly expressed in pseudomales. GO enrichment analysis results indicate that wnt signaling pathways and oocyte maturation are two terms enriched in female. At the protein level, Tandem Mass Tags analysis revealed that 324 proteins differed in their relative abundance between pseudomales and females. KEGG analysis found that pseudo-highly expressed proteins were enriched in the ubiquitin mediated proteolysis pathway. For integrative analysis, the Spearman correlation coefficient between the transcriptome and proteome was 0.59. Among 52 related genes, 46 DEGs (88%) were well matched in their levels of change in protein abundance. These findings reveal major active pathways in female and pseudomale gonads after sex reversal and provide new insights into molecular mechanisms associated with sex reversal regulatory network.