Molecular characterization and functional divergence of two Gadd45g homologs in sex determination in half-smooth tongue sole (Cynoglossus semilaevis)

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.

Global expression response of genes in sex-undifferentiated Nile tilapia gonads after exposure to trace letrozole

The aromatase inhibitor letrozole can be found in rivers, effluents, and even drinking water. Studies have demonstrated that letrozole affects various metabolic pathways and may cause reproductive toxicity, especially in fish exposed during development. However, studies on the effect of a low concentration of letrozole at the whole-gonad transcriptomic level in the early stage of fish sexual development have not been investigated. The aim of our study was to explore the potential effects of a low concentration of letrozole on the gonad transcriptome of Nile tilapia at an early stage of sexual development. In this study, 9 dpf (days postfertilization) Nile tilapia were exposed to trace letrozole for 12 days. Letrozole exposure from 9 dpf to 21 dpf persistently altered phenotypic sex development and induced the male-biased sex ratio. The transcriptome results showed that 1173 differentially expressed genes (DEGs) were present in the female control vs 1.5 μg/L letrozole-treated female comparison group and that 1576 DEGs were present in the 1.5 μg/L letrozole-treated female vs male control comparison group. Differentially expressed gene enrichment analysis revealed several crucial pathways, including the drug metabolism-cytochrome P450 pathway, the ErbB-PI3K/Akt/mTOR pathway, and the calcium signalling pathway. Further analysis of these identified DEGs indicated that some key genes correlated with metabolism and epigenetic regulation were significantly affected by letrozole, such as UDP-glucuronosyltransferase (Ugt), glutathione S-transferase omega-1 (Gsto1), lysine-specific demethylase 6bb (Kdm6bb, original name is Kdm6a), jumonji and AT-rich interaction domain containing 2 (Jarid2b, original name is Jarid2), growth arrest and DNA damage inducible gamma (Gadd45g), and chromobox protein 7 (Cbx7). The qRT-PCR validation results for twelve DEGs showed that the Pearson's correlation of the log₁₀^{fold change} values between the qPCR and RNA-Seq results was 0.90, indicating the accuracy and reliability of the RNA-Seq results. Our study is the first to report the effect of letrozole on the transcriptome of gonads from fish during early-stage sexual development. These findings will be useful for understanding the toxic effects and molecular mechanisms of letrozole exposure at the early stage of gonad development on the sexual development of aquatic organisms.

A combination of genome-wide association study screening and SNaPshot for detecting sex-related SNPs and genes in Cynoglossus semilaevis

Chinese tongue sole (Cynoglossus semilaevis) males and females exhibit great differences in growth rate and appearance. The species is heterogametic (ZW/ZZ) and has sex-reversed "pseudomales" that are genetically female and physiologically male. In this study, we identified eight sex-specific single nucleotide polymorphism (SNP) markers for the sex identification of C. semilaevis by using a combination of genome-wide association study (GWAS) screening and SnaPshot validation. Candidate SNPs were screened using genotyping by sequencing to perform GWAS of the differential SNPs between the sexes of C. semilaevis. The SNP loci were amplified using a multiplex PCR system and detected via SNaPshot, which enables multiplexing of up to 30-40 SNPs in a single assay and ensures high accuracy of the results. The molecular markers detected in our study were used to successfully identify normal males and pseudomales from 45 caught and 40 cultured C. semilaevis specimens. Linkage disequilibrium analysis showed that the eight SNP loci were related to each other, with a strong linkage. Moreover, we investigated the expression of prdm6 mRNA containing a missense SNP and confirmed that the gene is differentially expressed in the gonads of the different sexes of C. semilaevis; the expression of prdm6 mRNA was significantly higher in the males than in the females and pseudomales. This means prdm6 may be related to sex differentiation in C. semilaevis.

Molecular characterization and expression profiles provide new insights into GATA5 functions in tongue sole (Cynoglossus semilaevis)

GATA5 is a member of the GATA transcription factor family, which serves essential roles in varieties of cellular functions and biological processes. In this study, we have accomplished the molecular cloning, bioinformatic analysis and preliminary function study of C. semilaevis GATA5. The full-length cDNA nucleotide sequence is 1955 bp, with a coding sequence of 1167 bp, which encodes a polypeptide of 388 amino acids. Homology, phylogenetic, gene structure and synteny analysis showed that C. semilaevis GATA5 was highly conserved among vertebrates. Tissue distribution pattern exhibited that C. semilaevis GATA5 was significantly expressed in heart, intestine, liver, kidney and gonad, with a sexual dimorphic feature observed in testis and ovary. Embryonic development expression profiles showed that C. semilaevis GATA5 transcripts increased at the blastula stage, and peaked at the heat-beating period. Strong signals were detected at spermatids of male testis and stage III oocytes of female ovary by ISH. The expression of C. semilaevis GATA5 was regulated by 17α-MT and E2 after hormone stimulation to the ovary. Together, all the results pointed out that GATA5 might play a vital role during gonadal maturation and the reproductive cycle of C. semilaevis. This study lays the foundation for further researches on the sex control breeding in tongue sole.

Molecular characterization of Pod1 during sex development in Chinese tongue sole (Cynoglossus semilaevis)

Pod1 encodes a Class II bHLH transcription factor involved in the development of a number of tissues such as gonad, spleen, lungs and heart. However, to date, little is known about its function in teleosts. In this study, we cloned and characterized Pod1 gene from Cynoglossus semilaevis. This gene contains three exons and two introns, with the full-length cDNA of 918 nucleotides that encodes a 183 amino acid protein with a conserved bHLH domain. Realtime quantitative PCR revealed that Pod1 was predominantly expressed in the testes of C. semilaevis. In different stages of testes development, Pod1 expression was undetectable up to 120 days after hatching (dah), and then increased at 210 dah and 1 year after hatching (yah). Furthermore, in situ hybridization (ISH) analysis revealed that Pod1 was mainly localized in the germ cells of testes, but was not detected in ovarian cells; which suggested its possible functions in spermatogenesis of C. semilaevis. The methylation profile analysis of Pod1 genomic sequence in the gonads showed that the differences in their putative promoter regions of Pod1 among ovary, male and pseudo-male testes were not obvious. Thus, further research might be needed to evaluate whether Pod1 expression is regulated by epigenetic level.

Characterization and expression pattern of r-spondin1 in Cynoglossus semilaevis

r-spondin1 (rspo1) encodes a secreted protein that is involved in the determination and differentiation of the mammalian ovary. However, little information is yet available for teleosts. Here, we identified a homologue of rspo1 in Cynoglossus semilaevis. The full-length cDNA of rspo1 had a length of 2,703 bp with an open reading frame of 834 bp, encoding a protein with a length of 277 amino-acids. rspo1 expression was detected via qRT-PCR in various tissues, and significant sexually dimorphic expression was observed in the gonads. Furthermore, ISH located rspo1 in germ cells such as spermatogonia, spermatocytes, spermatids, spermatozoa, and oocytes, as well as in somatic cells of the gonads. Following knockdown of rspo1 in an ovarian cell line, the expressions of wnt4a, β-catenin, foxl2, and StAR were highly affected; wnt4a and β-catenin were significantly downregulated, whereas foxl2 and StAR were significantly upregulated. In summary, these data suggest that rspo1 may be involved in the regulation of ovarian development and differentiation through a conserved pathway, while the function of the gene in the testis remains elusive.

Locus Mapping, Molecular Cloning, and Expression Analysis of rps6kb2, a Novel Metamorphosis-Related Gene in Chinese Tongue Sole (Cynoglossus semilaevis)

Flatfish metamorphosis denotes the extraordinary transformation of a symmetric pelagic larva into an asymmetric benthic juvenile. This unique process involves eye migration, a 90° rotation in posture, and asymmetrical pigmentation for adaptation to a benthic lifestyle. In the present study, we used genetics to map a metamorphosis-related locus (q-10M) in the male linkage group (LG10M), a small interval of 0.9 cM corresponding to a 1.8 M-bp physical area in chromosome 9 in the Chinese tongue sole (Cynoglossus semilaevis). Combined with single-marker analysis, ribosomal protein S6 kinase 2 (rps6kb2) a member of the family of AGC kinases was identified as a novel metamorphosis-related candidate gene. Its expression pattern during metamorphosis was determined by quantitative RT-PCR and whole-mount in situ hybridization analysis. rps6kb2 gene was significantly expressed in metamorphic climax stage larvae and distributed in all the tissues transforming during metamorphosis, including tail, jaw, eye and skin of larvae. The results suggest that rps6kb2 has a general role in tissue transformations during flatfish metamorphosis including tail changes, skull remodeling, eye migration, and asymmetrical pigmentation.

Two Figla homologues have disparate functions during sex differentiation in half-smooth tongue sole (Cynoglossus semilaevis)

Figla is a germ-cell-specific transcription factor associated with ovary development and differentiation. In vertebrates, one transcriptional form of Figla is commonly found. However, besides the common form of this gene (named Figla_tv1), a new transcriptional form (named Figla_tv2) was identified in half-smooth tongue sole (Cynoglossus semilaevis). The full-length cDNA of Figla_tv1 was 1057 bp long with a 591-bp open reading frame encoding a predicted 196 amino acid protein, while Figla_tv2 encoded a 125 amino acid protein. Figla_tv1 and Figla_tv2 expression in various tissues was detected by qRT-PCR. Figla_tv1 was expressed mainly in ovary, skin and liver, while Figla_tv2 was expressed in all examined tissues. In the gonads, Figla_tv1 was expressed in ovary, while Figla_tv2 was predominately expressed in testis of pseudomales. Further, in situ hybridization located Figla_tv1 only in oocytes and Figla_tv2 mainly in germ cells of pseudomale testis. After knocking down Figla_tv2 in a pseudomale testis cell line, the expression of two steroid hormone-encoding genes, StAR and P450scc, was significantly up-regulated (P < 0.05). Our findings suggest that Figla_tv1 has a conserved function in folliculogenesis, as in other vertebrates, and that Figla_tv2 may have a role in the spermatogenesis of pseudomales by regulating the synthesis and metabolism of steroid hormones.