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

De Novo Assembly, Characterization, and Comparative Transcriptome Analysis of Mature Male and Female Gonads of Rabbitfish (Siganus oramin) (Bloch & Schneider, 1801)

The rabbitfish, Siganus oramin, is a commercially important table fish in southeastern China. However, there have been few studies on its gonad development and reproduction regulation. Comparative transcriptome analysis was first performed on adult male and female gonads of S. oramin. In total, 47,070 unigenes were successfully assembled and 22,737 unigenes were successfully annotated. Through comparative transcriptome analysis of male and female gonads, a total of 6722 differentially expressed genes were successfully identified, with 3528 upregulated genes and 3154 downregulated genes in the testes. In addition, 39 differentially expressed reproduction-related genes were identified. Finally, quantitative real-time PCR was used to validate the expression levels of several differentially expressed genes. These results provide important data for further studying the function of reproduction-related genes and the molecular mechanism regulating gonad development and reproduction in S. oramin.

Integration of Whole-Genome Resequencing and Transcriptome Sequencing Reveals Candidate Genes in High Glossiness of Eggshell

Eggshell gloss is an important characteristic for the manifestation of eggshell appearance. However, no study has yet identified potential candidate genes for eggshell gloss between high-gloss (HG) and low-gloss (LG) chickens. The aim of this study was to perform a preliminary investigation into the formation mechanism of eggshell gloss and to identify potential genes. The eggshell gloss of 300-day-old Rhode Island Red hens was measured from three aspects. Uterine tissues of the selected HG and LG (n = 5) hens were collected for RNA-seq. Blood samples were also collected for whole-genome resequencing (WGRS). RNA-seq analysis showed that 150 differentially expressed genes (DEGs) were identified in the uterine tissues of HG and LG hens. These DEGs were mainly enriched in the calcium signaling pathway and the neuroactive ligand-receptor interaction pathway. Importantly, these two pathways were also significantly enriched in the WGRS analysis results. Further joint analysis of WGRS and RNA-seq data revealed that 5-hydroxytryptamine receptor 1F (HTR1F), zinc finger protein 536 (ZNF536), NEDD8 ubiquitin-like modifier (NEDD8), nerve growth factor (NGF) and calmodulin 1 (CALM1) are potential candidate genes for eggshell gloss. In summary, our research provides a reference for the study of eggshell gloss and lays a foundation for improving egg glossiness in layer breeding.

Transcriptome Analysis Revealed the Early Heat Stress Response in the Brain of Chinese Tongue Sole (Cynoglossus semilaevis)

As a common influencing factor in the environment, temperature greatly influences the fish that live in the water all their life. The essential economic fish Chinese tongue sole (Cynoglossus semilaevis), a benthic fish, will experience both physiological and behavioral changes due to increases in temperature. The brain, as the central hub of fish and a crucial regulatory organ, is particularly sensitive to temperature changes and will be affected. However, previous research has mainly concentrated on the impact of temperature on the gonads of C. semilaevis. Instead, our study examines the brain using transcriptomics to investigate specific genes and pathways that can quickly respond to temperature changes. The fish were subjected to various periods of heat stress (1 h, 2 h, 3 h, and 5 h) before extracting the brain for transcriptome analysis. After conducting transcriptomic analyses, we identified distinct genes and pathways in males and females. The pathways were mainly related to cortisol synthesis and secretion, neuroactive ligand-receptor interactions, TGF beta signaling pathway, and JAK/STAT signaling pathway, while the genes included the HSP family, tshr, c-fos, c-jun, cxcr4, camk2b, and igf2. Our study offers valuable insights into the regulation mechanisms of the brain's response to temperature stress.

Transcriptome analysis of hypothalamus and pituitary tissues reveals genetic mechanisms associated with high egg production rates in Changshun green-shell laying hens

BACKGROUND

Changshun green-shell laying hens are unique to the Guizhou Province, China, and have high egg quality but relatively low yield. Egg production traits are regulated by the hypothalamus-pituitary-ovary axis. However, the underlying mechanism remains unclear. Thus, we conducted RNA sequencing of hypothalamic and pituitary tissues from low- and high-yielding Changshun green-shell laying hens to identify critical pathways and candidate genes involved in controlling the egg production rate.

RESULTS

More than 39 million clean reads per sample were obtained, and more than 82% were mapped to the Gallus gallus genome. Further analysis identified 1,817 and 1,171 differentially expressed genes (DEGs) in the hypothalamus and pituitary, respectively. Nineteen DEGs were upregulated in both the hypothalamus and pituitary of high-yielding chickens. The functions of these DEGs were mainly associated with ion transport or signal transduction. Gene set enrichment analysis revealed that the pathways enriched in the hypothalamus were mainly associated with gonadotropin-releasing hormone (GnRH) secretion, neurotransmitter release, and circadian rhythms. The pathways enriched in the pituitary were mainly associated with GnRH secretion, energy metabolism, and signal transduction. Five and four DEGs in the hypothalamus and pituitary, respectively, were selected randomly for qRT-PCR analysis. The expression trends determined via qRT-PCR were consistent with the RNA-seq results.

CONCLUSIONS

The current study identified 19 DEGs upregulated in both the hypothalamus and pituitary gland, which could provide an important reference for further studies on the molecular mechanisms underlying egg production in Changshun green-shell laying hens. In addition, enrichment analysis showed that GnRH secretion and signal transduction, especially neurotransmitter release, play crucial roles in the regulation of egg production.

Molecular biological, physiological, cytological, and epigenetic mechanisms of environmental sex differentiation in teleosts: A systematic review

Human activities have been exerting widespread stress and environmental risks in aquatic ecosystems. Environmental stress, including temperature rise, acidification, hypoxia, light pollution, and crowding, had a considerable negative impact on the life histology of aquatic animals, especially on sex differentiation (SDi) and the resulting sex ratios. Understanding how the sex of fish responds to stressful environments is of great importance for understanding the origin and maintenance of sex, the dynamics of the natural population in the changing world, and the precise application of sex control in aquaculture. This review conducted an exhaustive search of the available literature on the influence of environmental stress (ES) on SDi. Evidence has shown that all types of ES can affect SDi and universally result in an increase in males or masculinization, which has been reported in 100 fish species and 121 cases. Then, this comprehensive review aimed to summarize the molecular biology, physiology, cytology, and epigenetic mechanisms through which ES contributes to male development or masculinization. The relationship between ES and fish SDi from multiple aspects was analyzed, and it was found that environmental sex differentiation (ESDi) is the result of the combined effects of genetic and epigenetic factors, self-physiological regulation, and response to environmental signals, which involves a sophisticated network of various hormones and numerous genes at multiple levels and multiple gradations in bipotential gonads. In both normal male differentiation and ES-induced masculinization, the stress pathway and epigenetic regulation play important roles; however, how they co-regulate SDi is unclear. Evidence suggests that the universal emergence or increase in males in aquatic animals is an adaptation to moderate ES. ES-induced sex reversal should be fully investigated in more fish species and extensively in the wild. The potential aquaculture applications and difficulties associated with ESDi have also been addressed. Finally, the knowledge gaps in the ESDi are presented, which will guide the priorities of future research.

De novo assembly, characterization and comparative transcriptome analysis of gonads reveals sex-biased genes in Coreoperca whiteheadi

The wild Coreoperca whiteheadi is considered as the primordial species in sinipercine fish, which has valuable genetic information. Unfortunately, C. whiteheadi was listed as a near-threatened species because of the environmental pollution, over-exploitation and species invasion. Therefore, more genetic information is needed to have a better understanding of gonadal development in C. whiteheadi. Here, the first gonadal transcriptomes analysis of C. whiteheadi was conducted and 277.14 million clean reads were generated. A total of 96,753 unigenes were successfully annotated. By comparing ovary and testis transcriptomes, a total of 21,741 differentially expressed genes (DEGs) were identified, of which 12,057 were upregulated and 9684 were downregulated in testes. Among them, we also identified about 53 differentially expressed sex-biased genes. Subsequently, the expression of twenty-four DEGs were confirmed by real-time fluorescence quantitative PCR. Furthermore, the histological analysis was conducted on ovaries and testes of one-year-old C. whiteheadi. Our results provided basic support for further studies on the function of sex-biased genes and the molecular mechanism of sex determination and reproduction in C. whiteheadi.

De Novo Assembly, Characterization and Comparative Transcriptome Analysis of the Gonads of Jade Perch (Scortum barcoo)

Due to the high meat yield and rich nutritional content, jade perch (Scortum barcoo) has become an important commercial aquaculture species in China. Jade perch has a slow growth rate, taking 3-4 years to reach sexual maturity, and has almost no difference in body size between males and females. However, the study of its gonad development and reproduction regulation is still blank, which limited the yield increase. Herein, the gonad transcriptomes of juvenile males and females of S. barcoo were identified for the first time. A total of 107,060 unigenes were successfully annotated. By comparing male and female gonad transcriptomes, a total of 23,849 differentially expressed genes (DEGs) were identified, of which 9517 were downregulated, and 14,332 were upregulated in the testis. In addition, a large number of DEGs involved in sex differentiation, gonadal development and differentiation and gametogenesis were identified, and the differential expression patterns of some genes were further verified using real-time fluorescence quantitative PCR. The results of this study will provide a valuable resource for further studies on sex determination and gonadal development of S. barcoo.

Genome-Wide Association Study of Growth and Sex Traits Provides Insight into Heritable Mechanisms Underlying Growth Development of Macrobrachium nipponense (Oriental River Prawn)

Male hybrid oriental river prawns grow significantly faster than hybrid females. In this study, the growth and sex traits of 181 individuals of Macrobrachium nipponense were recorded, and each individual genotype was evaluated using the 2b-RAD sequencing method. The genetic parameters for growth and sex traits were estimated. A genome-wide association analysis (GWAS) of these traits was performed. In total, 18 growth-related SNPs were detected from 12 chromosomes using a mixed linear model. The most significant loci of weight are located on the position of the SNP (102638935, chromosome 13), which can explain 11.87% of the phenotypic variation. A total of 11 significant SNPs were detected on four chromosomes associated with sex trait (three on chromosome 4, one on chromosome 7 and seven on chromosome 17). The heritability of this trait is 0.8998 and belongs to the range of ultra-high heritability. Genetic correlations were prevalent among the 11 traits examined, the genetic coefficient between sex and body weight reached a significant level of −0.23. This study is the first GWAS for sex of binary and growth traits in oriental river prawn. Our results provide a set of markers for the genetic selection of growth traits and help us to further understand the genetic mechanisms of growth in Macrobrachium nipponense.

The Role of the Heat Shock Cognate Protein 70 Genes in Sex Determination and Differentiation of Chinese Tongue Sole (Cynoglossus semilaevis)

Fish sex determination can be affected by environmental temperature. This process relies on temperature-sensitive proteins such as heat shock proteins (HSPs). Our previous work found that heat shock cognate proteins (HSCs) may participate in high-temperature associated sex reversal of Chinese tongue sole (Cynoglossus semilaevis). However, the role of hsc genes in responding to high temperature and affecting sex determination/differentiation remains unclear. Here, by using C. semilaevis as model, we identified hsc70 and hsc70-like. hsc70 was abundant in the gonads with a testicular-higher expression at all gonadal development stages except for 6 months post fertilization (mpf). Intriguingly, hsc70-like showed higher expression in testes from 6 mpf on. Both long-term heat treatment during the temperature-sensitive sex-determining period and short-term heat stress at the end of this period caused different expression of hsc70/hsc70-like between sexes. The dual-luciferase assay results also suggested that these genes can respond to high temperature rapidly in vitro. Heat treatment of C. semilaevis testis cells overexpressed with hsc70/hsc70-like could affect the expression of sex-related genes sox9a and cyp19a1a. Our results indicated that hsc70 and hsc70-like were key regulators linking external high-temperature signals with sex differentiation in vivo and provide a new idea for understanding the mechanism by which high temperature affects sex determination/differentiation in teleosts.

Transcriptomic Analysis Revealed Candidate Genes Involved in Pseudomale Sperm Abnormalities in Chinese Tongue Sole (Cynoglossus semilaevis)

Chinese tongue sole (Cynoglossus semilaevis) has a ZZ/ZW sex determination system, but the genotypic female (ZW) can be sex-reversed into phenotypic males, namely, pseudomales. Pseudomale fish can produce only Z-type sperm but not W sperm. However, the molecular mechanism is unclear. To screen the key genes involved in pseudomale sperm abnormalities, we analysed the transcriptomic profiles of pseudomale and male sperm. In comparison to male sperm, 592 differentially expressed genes (DEGs) were identified in pseudomale sperm, including 499 upregulated and 93 downregulated genes. KEGG analysis indicated that the FoxO signalling pathway, especially the foxo3a and foxo6-like genes, may play an important role in spermatogenesis. The DEGs were mainly distributed on sex chromosomes, with 158 downregulated genes on the Z chromosome and 41 upregulated genes on the W chromosome. A specific area (14-15 M) on the Z chromosome was identified, which enriched eight DEGs inside the ~1 M region. In addition, there were five gene alleles on the sex chromosomes, which showed the opposite transcription pattern (upregulated for the W allele, downregulated for the Z allele). This study has provided valuable data for screening candidate genes involved in the pseudomale sperm abnormality.

Improved biomarker discovery through a plot twist in transcriptomic data analysis

Background

Transcriptomic analysis is crucial for understanding the functional elements of the genome, with the classic method consisting of screening transcriptomics datasets for differentially expressed genes (DEGs). Additionally, since 2005, weighted gene co-expression network analysis (WGCNA) has emerged as a powerful method to explore relationships between genes. However, an approach combining both methods, i.e., filtering the transcriptome dataset by DEGs or other criteria, followed by WGCNA (DEGs + WGCNA), has become common. This is of concern because such approach can affect the resulting underlying architecture of the network under analysis and lead to wrong conclusions. Here, we explore a plot twist to transcriptome data analysis: applying WGCNA to exploit entire datasets without affecting the topology of the network, followed with the strength and relative simplicity of DEG analysis (WGCNA + DEGs). We tested WGCNA + DEGs against DEGs + WGCNA to publicly available transcriptomics data in one of the most transcriptomically complex tissues and delicate processes: vertebrate gonads undergoing sex differentiation. We further validate the general applicability of our approach through analysis of datasets from three distinct model systems: European sea bass, mouse, and human.

Results

In all cases, WGCNA + DEGs clearly outperformed DEGs + WGCNA. First, the network model fit and node connectivity measures and other network statistics improved. The gene lists filtered by each method were different, the number of modules associated with the trait of interest and key genes retained increased, and GO terms of biological processes provided a more nuanced representation of the biological question under consideration. Lastly, WGCNA + DEGs facilitated biomarker discovery.

Conclusions

We propose that building a co-expression network from an entire dataset, and only thereafter filtering by DEGs, should be the method to use in transcriptomic studies, regardless of biological system, species, or question being considered.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01398-w.

Transcriptome Profiling and Expression Localization of Key Sex-Related Genes in a Socially-Controlled Hermaphroditic Clownfish, Amphiprion clarkii

Clownfish can be an excellent research model for investigating the socially-controlled sexual development of sequential hermaphrodite teleosts. However, the molecular cascades underlying the social cues that orchestrate the sexual development process remain poorly understood. Here, we performed a comparative transcriptomic analysis of gonads from females, males, and nonbreeders of Amphiprion clarkii, which constitute a complete social group, allowing us to investigate the molecular regulatory network under social control. Our analysis highlighted that the gonads of nonbreeders and males exhibited high similarities but were far from females, both in global transcriptomic profiles and histological characteristics, and identified numerous candidate genes involved in sexual development, some well-known and some novel. Significant upregulation of cyp19a1a, foxl2, nr5a1a, wnt4a, hsd3b7, and pgr in females provides strong evidence for the importance of steroidogenesis in ovarian development and maintenance, with cyp19a1a playing a central role. Amh and sox8 are two potential key factors that may regulate testicular tissue development in early and late stages, respectively, as they are expressed at higher levels in males than in females, but with slightly different expression timings. Unlike previous descriptions in other fishes, the unique expression pattern of dmrt1 in A. clarkii implied its potential function in both male and female gonads, and we speculated that it might play promoting roles in the early development of both testicular and ovarian tissues.

A Potential Role of the Translation Elongation Factor eef1a1 in Gonadal High-Temperature Perception in Chinese Tongue Sole (Cynoglossus semilaevis)

Simple Summary

The eukaryotic elongation factor 1 alpha (eef1a) gene is vital for protein translation by delivering aminoacylated tRNAs to the A/P site of the ribosome via the GTP-dependent reaction. Here, the Chinese tongue sole (Cynoglossus semilaevis) eef1a1 gene was identified, and its potential role in gonadal high-temperature perception was assessed. The full-length sequence of eef1a1 cDNA was 1809 base pair (bp) encoding a putative protein of 461 amino acids. The expression levels of eef1a1 in the ovary were significantly higher than that in the testis from 6 mpf to 3 ypf. Under high-temperature induction during sex differentiation, eef1a1 was significantly down-regulated in males, while the difference was not detected in females. Furthermore, the rapid response of eef1a1 to environmental high temperature was assessed in vitro. Our findings suggest that C. semilaevis eef1a1 might be essential for the molecular response regulatory network of external temperature affecting internal sex differentiation.

Abstract

The eukaryotic translation elongation factor 1 alpha (eef1a) gene has a well-defined role in protein synthesis. However, its role in external temperature perception and internal sex differentiation and development is still unclear. In this study, eef1a1 was identified and functionally analyzed in Chinese tongue sole (Cynoglossus semilaevis). The eef1a1 cDNA, 1809 bp in length, had a 1386 bp open reading frame (ORF) that encoded a 461 amino acid polypeptide containing one EF-1_alpha domain. eef1a1 expression levels were investigated across different tissues and during gonadal development. In the gonad, eef1a1 showed a sexually dimorphic expression pattern with a statistically higher expression level in the ovary than in the testis from 6 months postfertilization to 3 years postfertilization. Under high temperature (28 °C) treatment during C. semilaevis sex differentiation (from 30 days postfertilization to 3 months postfertilization), eef1a1 was statistically down-regulated in males, while the difference was not detected in females. In addition, the dual-luciferase assay exhibited that eef1a1 can respond to high temperature rapidly. Based on these results, C. semilaevis eef1a1 might have a dual role in the perception of external temperature changes and sex differentiation regulation.

Sex reversal explains some, but not all, climate-mediated sex ratio variation within a viviparous reptile

Evolutionary transitions in sex-determining systems have occurred frequently yet understanding how they occur remains a major challenge. In reptiles, transitions from genetic to temperature-dependent sex determination can occur if the gene products that determine sex evolve thermal sensitivity, resulting in sex-reversed individuals. However, evidence of sex reversal is limited to oviparous reptiles. Here we used thermal experiments to test whether sex reversal is responsible for differences in sex determination in a viviparous reptile, Carinascincus ocellatus, a species with XY sex chromosomes and population-specific sex ratio response to temperature. We show that sex reversal is occurring and that its frequency is related to temperature. Sex reversal was unidirectional (phenotypic males with XX genotype) and observed in both high- and low-elevation populations. We propose that XX-biased genotypic sex ratios could produce either male- or female-biased phenotypic sex ratios as observed in low-elevation C. ocellatus under variable rates of XX sex reversal. We discuss reasons why sex reversal may not influence sex ratios at high elevation. Our results suggest that the mechanism responsible for evolutionary transitions from genotypic to temperature-dependent sex determination is more complex than can be explained by a single process such as sex reversal.

Comprehensive Transcriptome Analysis of Gonadal and Somatic Tissues for Identification of Sex-Related Genes in the Largemouth Bass Micropterus salmoides

Largemouth bass (Micropterus salmoides) is an economically important fish. It can spawn many times during a breeding season, and there are no obvious morphological characteristics to distinguish male and female juvenile fish. So far, little is known about the genes regulating their sexual development in this species. Here, we performed RNA sequencing (RNA-Seq) analysis of the testis, ovary, and somatic tissue to identify sex-related genes in the largemouth bass. A total of 51,672 unigenes were obtained via the transcriptome analysis, and 5900 differential expression genes (DEGs), including 3028 up-regulated and 2872 down-regulated DEGs, were obtained in the somatic tissue, testis, and ovary. DEGs were retrieved by making comparisons: somatic tissue vs testis (1733-up and 1382-down), testis vs ovary (841-up and 807-down), and ovary vs somatic tissue (454-up and 683-down). Finally, functional annotation identified 22 key sex-related DEGs, including 13 testis-biased DEGs (dmrt1, cyp11b1, sox9, spata4, spata22, spata17, fshr, fem-1a, wt1, daz1, amh, vasa, and piwi1) and 9 ovary-biased DEGs (foxl2, gdf9, zp3, sox3, cyp19a, bmp15, fem-1b, fig. la, and piwi2). This result was further confirmed by the tissue expression detection via RT-PCR and RT-qPCR. Protein-protein interacting (PPI) network analysis revealed that the testis-specific dmrt1 interacts directly with the testis-biased DEGs (cyp11b1 and spata4) and the ovary-biased DEGs (foxl2, gdf9, zp3, sox3, cyp19a, and bmp15), suggesting that the dmrt1 as a sex-determining gene can play a dual role through inducing the testis-biased DEGs and inhibiting the ovary-biased DEGs during the testicular development. Our present results provide useful molecular data for a better understanding of sexual development in the largemouth bass.

Developmental changes in gene expression and gonad morphology during sex differentiation in Atlantic salmon (Salmo salar)

The Atlantic salmon (Salmo salar) is a globally important species for its value in fisheries and aquaculture, and as a research model. In order to characterise aspects of sex differentiation at the morphological and mRNA level in this species, the present study examined developmental changes in gonad morphology and gene expression in males and females between 0 and 79 days post hatch (dph). Morphological differentiation of the ovary (indicated by the formation of germ cell cysts) became apparent from 52 dph. By 79 dph, ovarian phenotype was evident in 100% of genotypic females. Testes remained in an undifferentiated-like state throughout the experiment, containing germ cells dispersed singularly within the gonadal region distal to the mesentery. There were no significant sex-related differences in gonad cross-section size, germ cell number or germ cell diameter during the experiment. The expression of genes involved in teleost sex differentiation (anti-müllerian hormone (amh), cytochrome P450, family 19, subfamily A, polypeptide 1a (cyp19a1a), forkhead box L2a (foxl2a), gonadal soma-derived factor (gsdf), r-spondin 1 (rspo1), sexually dimorphic on the Y chromosome (sdY)), retinoic acid-signalling (aldehyde dehydrogenase 1a2 (aldh1a2), cytochrome P450 family 26 a1 (cyp26a1), cytochrome P450 family 26 b1 (cyp26b1), t-box transcription factor 1 (tbx1a)) and neuroestrogen production (cytochrome P450, family 19, subfamily A, polypeptide 1b (cyp19a1b)) was investigated. Significant sex-related differences were observed only for the expression of amh, cyp19a1a, gsdf and sdY. In males, amh, gsdf and sdY were upregulated from 34, 59 and 44 dph respectively. In females, cyp19a1a was upregulated from 66 dph. Independent of sex, foxl2a expression was highest at 0 dph and had reduced ∼ 47-fold by the time of morphological sex differentiation at 52 dph. This study provides new insights into the timing and sequence of some physiological changes associated with sex differentiation in Atlantic salmon. These findings also reveal that some aspects of the mRNA sex differentiation pathways in Atlantic salmon are unique compared to other teleost fishes, including other salmonids.

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.

A Transcriptomic Analysis of Gonads from the Low-Temperature-Induced Masculinization of Takifugu rubripes

Simple Summary

Our study analyzed the differentiation of transcriptomes in normal and sex-reverse Takifugu rubripes, and screened out 13 differentially expressed genes related to sex differentiation. This is the first report on the gonadal transcriptome of pseudo-males in Takifugu rubripes. Our results provide an important contribution to the molecular mechanism of masculinization in a cultured fish subject to low-temperature treatment.

Abstract

The phenotypic sex of fish is usually plastic. Low-temperature treatment induces the masculinization of Takifugu rubripes, resulting in pseudo-males (PM) with the physiological sex of a male (M) and genetic sex of a female (F). For a comparison of gonadal transcriptomes, we collected gonads from three groups of T. rubripes (F, M, and PM) for high-throughput transcriptome sequencing. The results provided 467,640,218 raw reads (70.15 Gb) and a total of 436,151,088 clean reads (65.43 Gb), with an average length of 150 bp. Only 79 differentially expressed genes (DEGs) were identified between F and PM, whereas 12,041 and 11,528 DEGs were identified between F and M, and PM and M, respectively. According to the functional annotation of DEGs, 13 DEGs related to gonadal development were screened (LOC101066759, dgat1, limk1, fbxl3, col6a3, fgfr3, dusp22b, svil, abhd17b, srgap3, tmem88b, bud4, and mustn10) which might participate in formating PM. A quantitative PCR of the DEGs confirmed the reliability of the RNA-seq. Our results provide an important contribution to the genome sequence resources for T. rubripes and insight into the molecular mechanism of masculinization in a cultured fish subject to low-temperature treatment.

Transcriptomic analysis revealed gene expression profiles during the sex differentiation of Chinese tongue sole (Cynoglossus semilaevis)

Sex differentiation in aquatic fish is important both for theoretical study and practical production, as growth dimorphism frequently appears in different sexes, especially in marine fish. The deciphered genome, identification of the male-determining gene dmrt1 and established genotypic sex screening method make Chinese tongue sole (Cynoglossus semilaevis) an ideal model to study sex differentiation in fish. In this study, comparative gonadal transcriptomic analyses were conducted for genetic females and males at 48, 68, and 108 days post hatching (dph), representing pre-, during- and post-gonadal differentiation stages, although the gonad is not completely differentiated and isolable in 48 and 68 dph individuals, while it is in 108 dph individuals. Altogether, 28 libraries were constructed, and a mean of 46.64 M clean reads was obtained. Differentially expressed gene (DEG) analysis revealed that 179 genes had similar expression patterns in males and females in all three stages. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis suggested that the enriched pathways included ubiquitin-mediated proteolysis, lysosomes, and RNA degradation. Moreover, weighted gene coexpression network analyses (WGCNA) identified 14 modules, one of which was closely correlated with female differentiation, exhibiting female-biased expression in all three stages (48, 68, 108 dph). An illustrated core gene interaction network of this module identified 50 genes, most of which are on W chromosomes. Six genes, including two ubiquitin conjugating enzymes, were selected for further investigation, and their female-biased expression was confirmed in even earlier stages, at 10 and 30 dph. These data facilitate our understanding of sex differentiation in fish and provide a genomic rationale for screening candidate genes (preferentially W-linked genes) that could be involved in the female differentiation process.

Developmental aspects of the hypothalamic-pituitary network related to reproduction in teleost fish

The hypothalamic-pituitary-gonadal axis is the main system that regulates reproduction in vertebrates through a complex network that involves different neuropeptides, neurotransmitters, and pituitary hormones. Considering that this axis is established early on life, the main goal of the present work is to gather information on its development and the actions of its components during early life stages. This review focuses on fish because their neuroanatomical characteristics make them excellent models to study neuroendocrine systems. The following points are discussed: i) developmental functions of the neuroendocrine components of this network, and ii) developmental disruptions that may impact adult reproduction. The importance of the components of this network and their susceptibility to external/internal signals that can alter their specific early functions and/or even the establishment of the reproductive axis, indicate that more studies are necessary to understand this complex and dynamic network.

Gene Expression of Takifugu rubripes Gonads During AI- or MT-induced Masculinization and E2-induced Feminization

Elucidating the global molecular changes that occur during aromatase inhibitor (AI)- or 17α-methyltestosterone (MT)-induced masculinization and estradiol-17β (E2)-induced feminization is critical to understanding the roles that endocrine and genetic factors play in regulating the process of sex differentiation in fish. Here, fugu larvae were treated with AI (letrozole), MT, or E2 from 25 to 80 days after hatching (dah), and gonadal transcriptomic analysis at 80 dah was performed. The expression of dmrt1, gsdf, foxl2, and other key genes (star, hsd3b1, cyp11c1, cyp19a1a, etc.) involved in the steroid hormone biosynthesis pathway were found be altered. The expression of dmrt1, gsdf, cyp19a1a, and foxl2 was further verified by quantitative polymerase chain reaction. In the control group, the expression of dmrt1 and gsdf was significantly higher in XY larvae than in XX larvae, while the expression of foxl2 and cyp19a1a was significantly higher in XX larvae than in XY larvae (P < .05). AI treatment suppressed the expression of foxl2 and cyp19a1a, and induced the expression of dmrt1 and gsdf in XX larvae. MT treatment suppressed the expression of foxl2, cyp19a1a, dmrt1, and gsdf in XX larvae. E2 treatment suppressed the expression of dmrt1 and gsdf, but did not restore the expression of foxl2 and cyp19a1a in XY larvae. The shared response following AI, MT, and E2 treatment suggested that these genes are essential for sex differentiation. This finding offers some insight into AI or MT-induced masculinization, and E2-induced femininization in fugu.

High temperature induced masculinization of zebrafish by down-regulation of sox9b and esr1 via DNA methylation

Elevated temperature could influence the sex differentiation by altering the expression of sex-related genes in fish. However, the underlying mechanisms by which the gene expression is altered remain poorly understood. Here, we aimed to explore the role of DNA methylation in sex differentiation of zebrafish (Danio rerio) in response to elevated temperature. The results showed that high temperature (33°C) exposure of fish from 20 to 30 days post fertilization (dpf), compared to normal temperature (28°C), resulted in male-biased sex ratio and decreased expression of female-related genes including cyp19a1a, sox9b and esr1. Meanwhile, the expressions of DNA methyltransferases dnmt3a1 and dnmt3a2, and the DNA methylation levels in sox9b and esr1 promoter were significantly increased by high temperature, strongly implying that DNA methylation is involved in high temperature-induced masculinization of zebrafish. Co-treatment with 5-aza-2'-deoxycytidine (a DNA methylation inhibitor) attenuated the high temperature-induced masculinizing effect, recovered the expression of esr1 and sox9b, suppressed the transcription of dnmt3a1 and dnmt3a2, and decreased the methylation of esr1 and sox9b promoter, further confirming that DNA methylation plays an important role in high temperature-induced masculinization of zebrafish. Furthermore, the methylation of sox9b promoter decreased the enrichment of transcription factor CREB (cAMP-responsive element binding proteins). Overall, these findings suggest that high temperature induce masculinization of zebrafish by down-regulation of female-related genes via DNA methylation, providing a new insight in understanding the epigenetic mechanism of thermal-mediated sex differentiation in fish.

Temperature-Induced Sex Reversal in Reptiles: Prevalence, Discovery, and Evolutionary Implications

Sex reversal is the process by which an individual develops a phenotypic sex that is discordant with its chromosomal or genotypic sex. It occurs in many lineages of ectothermic vertebrates, such as fish, amphibians, and at least one agamid and one scincid reptile species. Sex reversal is usually triggered by an environmental cue that alters the genetically determined process of sexual differentiation, but it can also be caused by exposure to exogenous chemicals, hormones, or pollutants. Despite the occurrence of both temperature-dependent sex determination (TSD) and genetic sex determination (GSD) broadly among reptiles, only 2 species of squamates have thus far been demonstrated to possess sex reversal in nature (GSD with overriding thermal influence). The lack of species with unambiguously identified sex reversal is not necessarily a reflection of a low incidence of this trait among reptiles. Indeed, sex reversal may be relatively common in reptiles, but little is known of its prevalence, the mechanisms by which it occurs, or the consequences of sex reversal for species in the wild under a changing climate. In this review, we present a roadmap to the discovery of sex reversal in reptiles, outlining the various techniques that allow new occurrences of sex reversal to be identified, the molecular mechanisms that may be involved in sex reversal and how to identify them, and approaches for assessing the impacts of sex reversal in wild populations. We discuss the evolutionary implications of sex reversal and use the central bearded dragon (Pogona vitticeps) and the eastern three-lined skink (Bassiana duperreyi) as examples of how species with opposing patterns of sex reversal may be impacted differently by our rapidly changing climate. Ultimately, this review serves to highlight the importance of understanding sex reversal both in the laboratory and in wild populations and proposes practical solutions to foster future research.

Sex bias miRNAs in Cynoglossus semilaevis could play a role in transgenerational inheritance

Alterations of non-coding RNA profiling in spermatozoa are candidate mechanisms related to changes in paternal environment and progeny. Transgenerational inheritance of sex in pseudomales of Cynoglossus semilaevis, a fish with significant sex dimorphism, is a typical example of non-Mendelian inheritance. In the present study, miRNA profiles of spermatozoa were compared between male and pseudomale of C. semilaevis. Differential miRNAs in sperm from F0 and F1 generation also provides clues for revealing the possible role of non-coding RNA mediated transgenerational inheritance. Four sexual bias miRNAs, dre-miR-26a-5p, dre-miR-27b-3p, dre-miR-125b-5p,pol-199a-5p, were identified and verified in F0 and F1 generation of C. semilaevis. All of them were highly expressed in male sperm compared with pseudomale sperm. Function of target genes indicates that target genes of these differential RNAs are highly correlated with sex differentiation, gametogenesis and maintenance of secondary sexual characteristics. In a word, identification of epigenetic markers in gametes has great prospects in predicting susceptibility and properties in offsprings, and providing an indicator of parentalgenetic property.

A simple and rapid method for fish sex identification based on recombinase-aided amplification and its use in Cynoglossus semilaevis

Fish sex identification is a basic technique of great importance for both fish genetic studies and fisheries. Due to the sexual reversal phenomenon in many fish species, a simple and rapid molecular identification method for fish genetic sex is urgently needed to suit versatile detection scenarios, such as point-of-need applications. In this study, we took Cynoglossus semilaevis as an example, established a recombinase-aided amplification (RAA)-based method for sex identification, and combined the RAA-detection with two result visualization approaches with distinct features, capillary electrophoresis (CE) and lateral flow dipstick (LFD). Specific primers and probe were designed to specifically detect the sex chromosome W of C. semilaevis in order to distinguish the genetic sex between males, pseudo-males and females. To evaluate the performance of our methods, the genetic sex for twenty-eight males, sixty-eight pseudo-males and fifty-four females were examined with the RAA-based method and classical PCR-based genotyping method, demonstrating the consistent results of sex identification between both methods. The RAA-LFD method is operationally simple, rapid (~ 30 min) and holds great potential for point-of-need applications of fish sex identification, including fishery fields. The method presented here could be effective for identifying fish gender with the ZW karyotype.

Transcriptome analysis of ovary tissues from low- and high-yielding Changshun green-shell laying hens

Background

Changshun green-shell laying hens are unique to Guizhou Province, China, and have high egg quality. Improving egg production performance has become an important breeding task, and in recent years, the development of high-throughput sequencing technology provides a fast and exact method for genetic selection. Therefore, we aimed to use this technology to analyze the differences between the ovarian mRNA transcriptome of low and high-yield Changshun green-shell layer hens, identify critical pathways and candidate genes involved in controlling the egg production rate, and provide basic data for layer breeding.

Results

The egg production rates of the low egg production group (LP) and the high egg production group (HP) were 68.00 ± 5.56 % and 93.67 ± 7.09 %, with significant differences between the groups (p < 0.01). Moreover, the egg weight, shell thickness, strength and layer weight of the LP were significantly greater than those of the HP (p < 0.05). More than 41 million clean reads per sample were obtained, and more than 90 % of the clean reads were mapped to the Gallus gallus genome. Further analysis identified 142 differentially expressed genes (DEGs), and among them, 55 were upregulated and 87 were downregulated in the ovaries. KEGG pathway enrichment analysis identified 9 significantly enriched pathways, with the neuroactive ligand-receptor interaction pathway being the most enriched. GO enrichment analysis indicated that the GO term transmembrane receptor protein tyrosine kinase activity, and the DEGs identified in this GO term, including PRLR, NRP1, IL15, BANK1, NTRK1, CCK, and HGF may be associated with crucial roles in the regulation of egg production.

Conclusions

The above-mentioned DEGs may be relevant for the molecular breeding of Changshun green-shell laying hens. Moreover, enrichment analysis indicated that the neuroactive ligand-receptor interaction pathway and receptor protein tyrosine kinases may play crucial roles in the regulation of ovarian function and egg production.

Supplementary Information

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

Where the Ends Meet: An Overview of Sex Determination in Atheriniform Fishes

Atheriniform fishes have recently emerged as attractive models for evolutionary, ecological, and molecular/physiological studies on sex determination. Many species in this group have marked temperature-dependent sex determination (TSD) and yet many species also have a sex determinant gene that provides a strong drive for male differentiation. Thus, in these species the 2 forms of sex determination that were once considered to be mutually exclusive, environmental (ESD) and genotypic (GSD) sex determination, can coexist at environmentally relevant conditions. Here, we review the current knowledge on sex determination in atheriniform fishes with emphasis on the molecular and physiological mechanisms of ESD and GSD, the coexistence and cross-talk between these 2 mechanisms, the possibility of extragonadal transduction of environmental information and/or extragonadal onset of sex determination, and the results of field studies applying novel tools such as otolith increment analysis and molecular markers of genetic sex developed for selected New World and Old World atheriniform species. We also discuss the existence of molecular and histological mechanisms to prevent the discrepant differentiation in parts of the gonads because of ambiguous or conflicting environmental and genetic signals and particularly the possibility that the female is the default state in these species.

High elevation increases the risk of Y chromosome loss in Alpine skink populations with sex reversal

The view that has genotypic sex determination and environmental sex determination as mutually exclusive states in fishes and reptiles has been contradicted by the discovery that chromosomal sex and environmental influences can co-exist within the same species, hinting at a continuum of intermediate states. Systems where genes and the environment interact to determine sex present the opportunity for sex reversal to occur, where the phenotypic sex is the opposite of that predicted by their sex chromosome complement. The skink Bassiana duperreyi has XX/XY sex chromosomes with sex reversal of the XX genotype to a male phenotype, in laboratory experiments, and in field nests, in response to exposure to cold incubation temperatures. Here we studied the frequency of sex reversal in adult populations of B. duperreyi in response to climatic variation, using elevation as a surrogate for environmental temperatures. We demonstrate sex reversal in the wild for the first time in adults of a reptile species with XX/XY sex determination. The highest frequency of sex reversal occurred at the highest coolest elevation location, Mt Ginini (18.46%) and decreased in frequency to zero with decreasing elevation. We model the impact of this under Fisher's frequency-dependent selection to show that, at the highest elevations, populations risk the loss of the Y chromosome and a transition to temperature-dependent sex determination. This study contributes to our understanding of the risks of extinction from climate change in species subject to sex reversal by temperature, and will provide focus for future research to test on-the-ground management strategies to mitigate the effects of climate in local populations.

Two transcriptionally distinct pathways drive female development in a reptile with both genetic and temperature dependent sex determination

How temperature determines sex remains unknown. A recent hypothesis proposes that conserved cellular mechanisms (calcium and redox; 'CaRe' status) sense temperature and identify genes and regulatory pathways likely to be involved in driving sexual development. We take advantage of the unique sex determining system of the model organism, Pogona vitticeps, to assess predictions of this hypothesis. P. vitticeps has ZZ male: ZW female sex chromosomes whose influence can be overridden in genetic males by high temperatures, causing male-to-female sex reversal. We compare a developmental transcriptome series of ZWf females and temperature sex reversed ZZf females. We demonstrate that early developmental cascades differ dramatically between genetically driven and thermally driven females, later converging to produce a common outcome (ovaries). We show that genes proposed as regulators of thermosensitive sex determination play a role in temperature sex reversal. Our study greatly advances the search for the mechanisms by which temperature determines sex.

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.

High temperature stress response is not sexually dimorphic at the whole-body level and is dependent on androgens to induce sex reversal

The understanding of the molecular and endocrine mechanisms behind environmentally-induced sex reversal in fish is of great importance in the context of predicting the potential effects of climate change, especially increasing temperature. Here, we demonstrate the global effects of high temperature on genome-wide transcription in medaka (Oryzias latipes) during early development. Interestingly, data analysis did not show sexual dimorphic changes, demonstrating that thermal stress is not dependent on genotypic sex. Additionally, our results revealed significant changes in several pathways under high temperature, such as stress response from brain, steroid biosynthesis, epigenetic mechanisms, and thyroid hormone biosynthesis, among others. These microarray data raised the question of what the exact molecular and hormonal mechanisms of action are for female-to-male sex reversal under high temperatures in fish. Complementary gene expression analysis revealed that androgen-related genes increase in females (XX) experiencing high water temperature. To test the involvement of androgens in thermal-induced sex reversal, an androgen antagonist was used to treat XX medaka under a high-temperature setup. Data clearly demonstrated failure of female-to-male sex reversal when androgen action is inhibited, corroborating the importance of androgens in environmentally-induced sex reversal.

Early response to heat stress in Chinese tongue sole (Cynoglossus semilaevis): performance of different sexes, candidate genes and networks

BACKGROUND

Temperature is known to affect living organisms and alter the expression of responsive genes, which affects a series of life processes, such as development, reproduction and metabolism. Several genes and gene families have been involved in high temperature responses, such as heat shock protein (hsp) family, Jumonji family and genes related to cortisol synthesis. Gonad is a vital organ related to the existence of a species. However, the comprehensive understanding of gonadal responses to environmental temperature is limited.

RESULTS

To explore the effects of environmental temperature on genes and gene networks in gonads, we performed acute heat treatment (48 h) on Chinese tongue sole (Cynoglossus semilaevis). Gonadal transcriptome analysis was conducted on females, pseudomales and males exposed to high (28 °C) and normal (22 °C) temperatures. A total of 1226.24 million clean reads were obtained from 18 libraries. Principal component analysis (PCA) and differentially expressed gene (DEG) analysis revealed different performance of sex responses to heat stress. There were 4565, 790 and 1117 specific genes altered their expression level in females, pseudomales and males, respectively. Of these, genes related to hsp gene family, cortisol synthesis and metabolism and epigenetic regulation were involved in early heat response. Furthermore, a total of 1048 DEGs were shared among females, pesudomales and males, which may represent the inherent difference between high and normal temperatures. Genes, such as eef1akmt3, eef1akmt4, pnmt and hsp family members, were found.

CONCLUSIONS

Our results depicted for the first time the gonadal gene expression under acute high temperature treatment in Chinese tongue sole. The findings may provide a clue for understanding the responses of genes and networks to environmental temperature.

Genome-wide investigation of Dmrt gene family in large yellow croaker (Larimichthys crocea)

The Dmrt (Doublesex and Mab-3 related transcription factor) gene family is a class of crucial transcription factors characterized by a conserved DM (Doublesex/Mab-3) domain. Previous researches indicate this gene family is involved in various physiological processes, especially in sex determination/differentiation and gonad development. Despite the vital roles of the Dmrt gene family in physiological processes, the comprehensive characterization and analysis of the dmrt genes in large yellow croaker (Larimichthys crocea), one of the most commercially important marine fish in China, have not been described. In this study, we performed the first genome-wide systematic analysis of L. crocea dmrt genes through the bioinformatics method. A total of seven members of the Dmrt gene family including Lcdmrt1, Lcdmrt2a, Lcdmrt2b, Lcdmrt3, Lcdmrt4, Lcdmrt5, and Lcdmrt6 were excavated based on the genome data of L. crocea. Further analysis revealed that the dmrt genes of L. crocea were distributed unevenly across four chromosomes. There were three dmrt genes (Lcdmrt1, Lcdmrt2a, and Lcdmrt3) on 3rd chromosome, one (Lcdmrt6) on 13th chromosome, one (Lcdmrt4) on 14th chromosome, two on (Lcdmrt5 and Lcdmrt2b) 17th chromosome. The gene structure analysis indicated that the number of introns of different dmrt genes of L. crocea had some differences: Lcdmrt1 had four introns, Lcdmrt2a, Lcdmrt2b, and Lcdmrt6 had two introns, Lcdmrt3, Lcdmrt4, and Lcdmrt5 had only one intron. The expression pattern analysis with published gonad transcriptome datasets and further confirmed by qRT-PCR revealed that these members of the Dmrt gene family except for Lcdmrt4 were all sexually dimorphic and preferred expressing in testis. Furthermore, the expression pattern analysis also revealed that the expression level of Lcdmrt1 and Lcdmrt6 was significantly higher than that of other members, suggesting that these two genes may play a more important role in testis. Overall, our studies provide a comprehensive insight into the Dmrt gene family members and a basis for the further study of their biological functions in L. crocea.

Transcriptome Profiling and Analysis of Genes Associated with High Temperature-Induced Masculinization in Sex-Undifferentiated Nile Tilapia Gonad

Artificially high temperatures during critical thermosensitive periods (TSPs) can induce the sex reversal of Nile tilapia (Oreochromis niloticus) females into pseudomales; Nile tilapia is a GSD + TE (genotypic plus temperature effects) fish species. Previous studies have shown that water temperature affects the expression levels of many genes in the gonad or brain in various teleost species. However, few studies on the effect of temperature at the whole-gonad transcriptomic level in the early stage of sex differentiation have been reported in fish species exhibiting GSD + TE. In this study, RNA-Seq was performed to characterize the transcriptomic profile and identify genes exhibiting temperature- and sex-biased expressions in the Nile tilapia gonad at 21 dpf. A total of 42 genes were found to be associated with both high-temperature treatment and sex development, as the expression levels of these genes differed in both FC (female control) vs MC (male control) and FC vs FT (high temperature-treated females in the TSP). Among these genes, the transcriptional alterations of many male sex determination and differentiation genes, such as Dmrt1, Gsdf, and the DNA damage-inducible protein GADD45 alpha, suggested that the male pathway is initiated after high-temperature treatment and that its initiation may play a role in high temperature-induced masculinization in Nile tilapia. The qRT-PCR validation results for thirteen differentially expressed genes showed that the Pearson's correlation of the log₁₀ fold change values between the qPCR and RNA-Seq results was 0.70 (p < 0.001), indicating the accuracy and reliability of the RNA-Seq results. Our study provides insights into how high-temperature treatment induces the sex reversal of Nile tilapia females.