Single Locus Maintains Large Variation of Sex Reversal in Half-Smooth Tongue Sole (Cynoglossus semilaevis)

Identification and Functional Analysis of Cynoglossus semilaevis Z-Linked E3 Ubiquitin Ligase rnf34

The high proportion of males in C. semilaevis hinders their industrial development. The genetic ZW individual can become a pseudomale by sex reversal. And the pseudomale can produce Z-sperm (with epigenetic information to cause sex reversal) while W-sperm is absent, which leads to an even higher male proportion in offspring. Recently, with the development of transcriptomic technologies, research on spermatogenesis in C. semilaevis has been focused on the ubiquitination pathway. In this study, we analyzed the function of the ubiquitin ligase rnf34 gene on the Z chromosome. A qPCR experiment showed that its expression level in the gonad was the highest among different tissues. In the ovary, the expression gradually increased with development from 40 days post-hatching (dph) to 1.5 years post-hatching (yph). In the testis, rnf34 showed increased expression from 40 dph to 6 months post-hatching (mpf) and stabilized up until 1.5 ypf. In situ hybridization showed that the mRNA of rnf34 was mainly distributed in the germ cells of the testis and the ovary. In vivo siRNA-mediated knockdown of the rnf34 gene in male fish affected the expression of a series of genes related to sex differentiation and spermatogenesis. These results provide genetic data on the molecular mechanisms of gonadal development and spermatogenesis in C. semilaevis.

Genome-wide association study reveals genomic loci of sex differentiation and gonadal development in Plectropomus leopardus

Introduction: Plectropomus leopardus, a commercially significant marine fish, is primarily found in the Western Pacific regions and along the coast of Southeast Asia. A thorough analysis of the molecular mechanisms involved in sex differentiation is crucial for gaining a comprehensive understanding of gonadal development and improving sex control breeding. However, the relevant fundamental studies of P. leopardus are relatively lacking. Methods: In this study, a genome-wide association study (GWAS) was conducted to investigate the genetic basis mechanism of sex differentiation and gonadal developmental traits in P. leopardus utilizing about 6,850,000 high-quality single-nucleotide polymorphisms (SNPs) derived from 168 individuals (including 126 females and 42 males) by the genome-wide efficient mixed-model association (GEMMA) algorithm. Results: The results of these single-trait GWASs showed that 46 SNP loci (-log10 p > 7) significantly associated with sex differentiation, and gonadal development traits were distributed in multiple different chromosomes, which suggested the analyzed traits were all complex traits under multi-locus control. A total of 1,838 potential candidate genes were obtained by considering a less-stringent threshold (-log10 p > 6) and ±100 kb regions surrounding the significant genomic loci. Moreover, 31 candidate genes were identified through a comprehensive analysis of significant GWAS peaks, gene ontology (GO) annotations, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, including taf7, ddx6, apoeb, sgk1, a2m, usf1, hsd3b7, dll4, xbp1, tet3, esr1, and gli3. These trait-associated genes have been shown to be involved in germline development, male sex differentiation, gonad morphogenesis, hormone receptor binding, oocyte development, male gonad development, steroidogenesis, estrogen-synthetic pathway, etc. Discussion: In the present study, multiple genomic loci of P. leopardus associated with sex differentiation and gonadal development traits were identified for the first time by using GWAS, providing a valuable resource for further research on the molecular genetic mechanism and sex control in P. leopardus. Our results also can contribute to understanding the genetic basis of the sex differentiation mechanism and gonadal development process in grouper fish.

CircRNA Identification and CircRNA-miRNA-mRNA Network in Cynoglossus semilaevis Sexual Size Dimorphism

Sexual size dimorphism (SSD), which is the sexual differences in body size, has been widely reported in various species including fishes. For Chinese tongue sole (Cynoglossus semilaevis), a flatfish exhibiting typically female-biased SSD, little is known for its epigenetic regulation mechanism, especially the role of circRNAs. Here, we identified the differently expressed abundances of circRNAs in females, males, and pseudo-males to explore the potential functions of circRNAs in Chinese tongue sole SSD. In total, 14,745 novel circRNAs were screened, among which 1461 DE circRNAs were identified from the brain, gonad, liver, and muscle in female, male, and pseudo-male individuals. The ceRNA network was subsequently constructed, including 10 circRNAs, 26 mRNAs, and 11 miRNAs. These DE mRNAs were mainly related to the mRNA surveillance pathway, metabolic pathways, and cellular senescence. Importantly, the ceRNA network has revealed that several circRNAs such as novel_circ_004374 and novel_circ_014597 may regulate homeodomain interacting protein kinase 2 (hipk2) expression by sponging miR-130-x. It is also worth exploring whether or how novel_circ_008696 regulates SET Domain Containing 2, histone lysine methyltransferase (setd2), which in turn affects the epigenetic patterns of different sexual individuals. The present study not only enriches the knowledge on the potential roles of circRNA in the physiological process, but also provides new clues for the explanation of fish SSD. In future studies, the precise function and involvement of circRNAs in female-biased SSD will require more efforts.

Environmentally-induced sex reversal in fish with chromosomal vs. polygenic sex determination

Sex ratio depends on sex determination mechanisms and is a key demographic parameter determining population viability and resilience to natural and anthropogenic stressors. There is increasing evidence that the environment can alter sex ratio even in genetically sex-determined species (GSD), as elevated temperature can cause female-to-male sex reversal (neomales). Alarmingly, neomales are being discovered in natural populations of several fish, amphibian and reptile species worldwide. Understanding the basis of neomale development is important for conservation biology. Among GSD species, it is unknown whether those with chromosomal sex determination (CSD), the most common system, will better resist the influence of high temperature than those with polygenic sex determination (PSD). Here, we compared the effects of elevated temperature in two wild zebrafish strains, Nadia (NA) and Ekkwill (EKW), which have CSD with a ZZ/ZW system, against the AB laboratory strain, which has PSD. First, we uncovered novel sex genotypes and the results showed that, at control temperature, the masculinization rate roughly doubled with the addition of each Z chromosome, while some ZW and WW fish of the wild strains became neomales. Surprisingly, we found that at elevated temperatures WW fish were just as likely as ZW fish to become neomales and that all strains were equally susceptible to masculinization. These results demonstrate that the Z chromosome is not essential for male development and that the dose of W buffers masculinization at the control temperature but not at elevated temperature. Furthermore, at the elevated temperature the testes of neomales, but not of normal males, contained more spermatozoa than at the control temperature. Our results show in an unprecedented way that, in a global warming scenario, CSD species may not necessarily be better protected against the masculinizing effect of elevated temperature than PSD species, and reveal genotype-by-temperature interactions in male sex determination and spermatogenesis.

Intensive masculinization caused by chronic heat stress in juvenile Cynoglossus semilaevis: Growth performance, gonadal histology and gene responses

The sea temperature has been observed to chronically increase during the past decades, leaving unpredictable influences to the marine biological resources. Thus, it is of vital significance to study the biological responses of ocean inhabited organisms with the artificially stimulated heat stress environment. Cynoglossus semilaevis provides us with an ideal model to study the influence of chronic heat stress on the sexual differentiation in marine teleosts for its genetic sex determination (GSD) + environmental effected (EE) sex determination system. In this study, the comparative experiment was conducted employing heated seawater (HT group) and ambient seawater (CT group) to cultivate juvenile C. semilaevis respectively. Significant differences were exhibited in growth performance and a delayed germ cell development effect was found in pseudomales formed under chronic heat stress. Using transcriptome analysis, the transcription profile of 55 days post fertilization (dpf) and 100 dpf juveniles' gonads were studied. A total of 47 libraries were constructed with an average mapping rate of 94.63% after assembling. GO and KEGG enrichment were proceeded using DEGs screened out between (1) pseudomale gonads at 55 dpf and 100 dpf in HT and CT group (2) pseudomale and female gonads at 55 dpf and 100 dpf in HT and CT group. Terms and pathways involved in steroid stimulation, reproduction ability, germ cell proliferation et al. were shed light on. The expression pattern of 29 DEGs including amh, hsp90b1, pgr et al. were also provided to supplement the results of functional enrichment. Weighted gene co-expression networks analysis (WGCNA) was constructed and hspb8-like, histone H2A.V were exhibited to play vital roles in the heat-induced masculinization. Our findings facilitate the understanding for transcriptional variations in intensive masculinization cause by chronic heat stress of C. semilaevis and provide referable study of the influences on the teleosts in elevated sea temperature.

Construction of a high-density linkage map and detection of sex-specific markers in Penaeus japonicus

Penaeus japonicus is one of the most important farmed shrimp species in many countries. Sexual dimorphism is observed in P. japonicus, in which females grow faster and larger than males; therefore, a unisexual female culture of P. japonicus could improve the efficiency of productivity. However, the genetic mechanisms underlying sex determination in P. japonicus are unclear. In this study, we constructed a high-density genetic linkage map of P. japonicus using genotyping-by-sequencing (GBS) technology in a full-sib family. The final map was 3,481.98 cM in length and contained 29,757 single nucleotide polymorphisms (SNPs). These SNPs were distributed on 41 sex-averaged linkage groups, with an average inter-marker distance of 0.123 cM. One haplotype, harboring five sex-specific SNPs, was detected in linkage group 1 (LG1), and its corresponding confidence interval ranged from 211.840 to 212.592 cM. Therefore, this high-density genetic linkage map will be informative for genome assembly and marker-assisted breeding, and the sex-linked SNPs will be helpful for further studies on molecular mechanisms of sex determination and unisexual culture of P. japonicus in the future.

Selection of functional EPHB2 genotypes from ENU mutated grass carp treated with GCRV

Background

N-ethyl-N-nitrosourea (ENU) mutagenesis is a useful method for the genetic engineering of plants, and the production of functional mutants in animal models including mice and zebrafish. Grass carp reovirus (GCRV) is a haemorrhagic disease of grass carp which has caused noteworthy losses in fingerlings over the last few years. To overcome this problem, we used ENU mutant grass carp in an attempt to identify functional resistance genes for future hereditary rearing projects in grass carp.

Results

This study used ENU-mutated grass carp to identify genetic markers associated with resistance to the haemorrhagic disease caused by GCRV. Bulked segregant analysis (BSA) was performed on two homozygous gynogenetic ENU grass carp groups who were susceptible or resistant to GCRV. This analysis identified 466,162 SNPs and 197,644 InDels within the genomes of these mixed pools with a total of 170 genes annotated in the associated region, including 49 genes with non-synonymous mutations at SNP sites and 25 genes with frame shift mutations at InDel sites. Of these 170 mutated genes, 5 randomly selected immune-related genes were shown to be more strongly expressed in the resistant group as compared to the susceptible animals. In addition, we found that one immune-related gene, EPHB2, presented with two heterozygous SNP mutations which altered the animal’s responded to GCRV disease. These SNPs were found in the intron region of EPHB2 at positions 5859 (5859^G > A) and 5968 (5968^G > A) and were significantly (p = 0.002, 0.003) associated with resistance to GCRV. These SNP sites were also shown to correlate with the GCRV-resistant phenotype in these ENU grass carp. We also evaluated the mortality of the different ENU fish genotypes in response to GCRV and the SNPs in EPHB2. The outcomes of these evaluations will be useful in future selections of GCRV-resistant genes for genetic breeding in grass carp.

Conclusion

Our results provide a proof of concept for the application of BSA-sequence analysis in detecting genes responsible for specific functional genotypes and may help to develop better methods for marker-assisted selection, especially for disease resistance in response to GCRV.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07858-x.

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.

Transcriptome of Gonads From High Temperature Induced Sex Reversal During Sex Determination and Differentiation in Chinese Tongue Sole, Cynoglossus semilaevis

The sex of Chinese tongue sole (Cynoglossus semilaevis) is determined by both genetic sex determination (GSD) and environmental sex determination (ESD), making it an ideal model to study the relationship between sex-determination and temperature. In the present study, transcriptomes of undifferentiated gonads from genetic females and males, as well as differentiated gonads from males, females, and pseudomales under high and normal temperature treatments were generated for comparative transcriptomic analysis. A mean of 68.24 M high-quality clean reads was obtained for each library. Differentially expressed genes (DEGs) between different sexes and environmental treatments were identified, revealing that the heat shock protein gene family was involved in the high temperature induced sex reversal. The Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways that were enriched in pseudomale and genetic female comparison included neuroactive ligand-receptor interaction, cortisol synthesis and secretion, and steroid hormone biosynthesis. Furthermore, weighted gene co-expression network analyses were conducted on all samples, and two modules were positive correlated with pseudomale under high temperature. An illustrated protein-protein interaction map of the module identified a hub gene, hsc70. These findings provide insights into the genetic network that is involved in sex determination and sexual differentiation, and improve our understanding of genes involved in sex reversal under high temperature.

Transcriptome Profiling Insights the Feature of Sex Reversal Induced by High Temperature in Tongue Sole Cynoglossus semilaevis

Sex reversal induced by temperature change is a common feature in fish. Usually, the sex ratio shift occurs when temperature deviates too much from normal during embryogenesis or sex differentiation stages. Despite decades of work, the mechanism of how temperature functions during early development and sex reversal remains mysterious. In this study, we used Chinese tongue sole as a model to identify features from gonad transcriptomic and epigenetic mechanisms involved in temperature induced masculinization. Some of genetic females reversed to pseudomales after high temperature treatment which caused the sex ratio imbalance. RNA-seq data showed that the expression profiles of females and males were significantly different, and set of genes showed sexually dimorphic expression. The general transcriptomic feature of pesudomales was similar with males, but the genes involved in spermatogenesis and energy metabolism were differentially expressed. In gonads, the methylation level of cyp19a1a promoter was higher in females than in males and pseudomales. Furthermore, high-temperature treatment increased the cyp19a1a promoter methylation levels of females. We observed a significant negative correlation between methylation levels and expression of cyp19ala. In vitro study showed that CpG within the cAMP response element (CRE) of the cyp19a1a promoter was hypermethylated, and DNA methylation decreased the basal and forskolin-induced activities of cyp19a1a promoter. These results suggested that epigenetic change, i.e., DNA methylation, which regulate the expression of cyp19a1a might be the mechanism for the temperature induced masculinization in tongue sole. It may be a common mechanism in teleost that can be induced sex reversal by temperature.

New locus reveals the genetic architecture of sex reversal in the Chinese tongue sole (Cynoglossus semilaevis)

Sex reversal in insects, amphibians, reptiles, and fishes is a complicated and interesting biological phenomenon. Sex reversal changes the sex ratio of populations and may complicate breeding schemes. In the Chinese tongue sole (Cynoglossus semilaevis), genetic females may change into pseudomales, thereby increasing aquaculture costs because of the lower growth rate of the males than that of the females. Here we identify a new locus associated with sex reversal; this single nucleotide polymorphism (SNP) is located in the third intron of the doublesex and mab-3 related transcription factor 1 (Dmrt1) gene on the Z chromosome (named Cyn_Z_8564889) and has two alleles, A and G. Cyn_Z_8564889 regulates sex reversal interactively with our previously detected SNP (Cyn_Z_6676874), with the genetic females simultaneously carrying the T allele of Cyn_Z_6676874 and the A allele of Cyn_Z_8564889 changing into pseudomales. Other Dmrt1 polymorphisms were detected, which formed two haplotypes. Two SNPs in the second exon of Dmrt1 result in amino acid changes, suggesting that Dmrt1 is essential in sex reversal. We also verified that pseudomales produce no or little W sperm. The interaction and linkage between Cyn_Z_6676874 and Cyn_Z_8564889 and the absence of W sperm from pseudomales unravel the genetic architecture of sex reversal in C. semilaevis.

Effects of different dietary DHA:EPA ratios on gonadal steroidogenesis in the marine teleost, tongue sole (Cynoglossus semilaevis)

The present study was conducted to investigate the effects of dietary DHA and EPA on gonadal steroidogenesis in mature females and males, with a feeding trial on tongue sole, a typical marine teleost with sexual dimorphism. Three experimental diets differing basically in DHA:EPA ratio, that is, 0·68 (diet D:E-0·68), 1·09 (D:E-1·09) and 2·05 (D:E-2·05), were randomly assigned to nine tanks of 3-year-old tongue sole (ten females and fifteen males in each tank). The feeding trail lasted for 90 d before and during the spawning season. Fish were reared in a flowing seawater system and fed to apparent satiation twice daily. Compared with diet D:E-0·68, diet D:E-1·09 significantly enhanced the oestradiol production in females, whereas diet D:E-2·05 significantly enhanced the testosterone production in males. In ovaries, diet D:E-1·09 induced highest mRNA expression of follicle-stimulating hormone receptor (FSHR), steroidogenic acute regulatory protein, 17α-hydroxylase (P450c17) and 3β-hydroxysteroid dehydrogenase (3β-HSD). In testes, diet 2·05 resulted in highest mRNA expression of FSHR, cholesterol side-chain cleavage enzyme, P450c17 and 3β-HSD. Fatty acid profiles in fish tissues reflected closely those of diets. Female fish had more gonadal EPA content but less DHA content than male fish, whereas there was a reverse observation in liver. In conclusion, the dietary DHA:EPA ratio, possibly combined with the dietary EPA:arachidonic acid ratio, differentially regulated sex steroid hormone synthesis in mature female and male tongue soles. Females seemed to require more EPA but less DHA for the gonadal steroidogenesis than males. The results are beneficial to sex-specific nutritive strategies in domestic teleost.