Transcriptome analysis of male and female Sebastiscus marmoratus

Comparative transcriptome analysis of gonads in male and female Pseudobagrus ussuriensis (Bagridae, Siluriformes)

As an important aquaculture fish in the Heilongjiang River Basin, Pseudobagrus ussuriensis has high economic value, and all-male culture is beneficial to the economic development of this fish. In this study, the transcriptomes of gonads in males and females were analyzed, and some genes related to gonad development were found. A total of 82,931 unigenes were found (average length 1504 bp, N50 1829 bp). In addition, 4689 differentially expressed genes (DEGs; including 1424 genes upregulated and 3265 genes downregulated in males) were identified. Some genes associated with testis development (such as Dmrt1 and Ropn1l) were significantly upregulated in males, while genes related to ovary development (such as Wnt2, PLC, Cyp19a, ZP3) were significantly downregulated in males, demonstrating that these genes have a crucial influence on gonad development in P. ussuriensis. Some signaling pathways related to gonad development were found, such as the Wnt pathway and oocyte meiosis. The results of RNA-seq obtained in this study provide theoretical data for elucidating the potential mechanism of gonad development of P. ussuriensis and reliable genomic data for the establishment of mono-sex breeding of P. ussuriensis.

The impact of benzo[a]pyrene on murine allergic airway inflammation via epigenetic remodeling

Simultaneous exposure to both BaP and house dust mites (HDM) has been shown to exacerbate pulmonary inflammation and hyperresponsiveness in a murine asthma model. The mechanistic insight into epigenetic inheritance for this effect, however, remains to be clarified. As such, in this study, we explore the molecular basis for the enhancement of asthma. Female BAL/C mice were intranasally administered HDM (25 µg in 25 μL saline) and/or BaP (10 μg/kg) every other day for 9 weeks. RNA sequencing and DNA methylation assessment were used to explore the underlying mechanism. Following simultaneous exposure to HDM and BaP, mice exhibited pulmonary inflammation and the transcript level of IL4i1b, muc4 and IL22ra2 that were associated with altered DNA methylation, suggesting that there may be an epigenetic basis for BaP-induced asthma exacerbation. Our data suggest that DNA methylation is a major epigenetic modification that accompanies airway remodeling associated with changes in the allergic mice.

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.

Comparative Transcriptome Analysis Revealed Genes Involved in Sexual and Polyploid Growth Dimorphisms in Loach (Misgurnus anguillicaudatus)

Simple Summary

Misgurnus anguillicaudatus not only exhibits sexual size dimorphism, but also shows polyploid size dimorphism. Here, we performed comparative transcriptome integration analysis of multiple tissues of diploid and tetraploid M. anguillicaudatus of both sexes. We found that differences in energy metabolism and steroid hormone synthesis levels may be the main causes of sexual and polyploidy growth dimorphisms of M. anguillicaudatus. Fast-growing M. anguillicaudatus (tetraploids, females) have higher levels of energy metabolism and lower steroid hormone synthesis and fatty acid degradation abilities than slow-growing M. anguillicaudatus (diploids, males).

Abstract

Sexual and polyploidy size dimorphisms are widespread phenomena in fish, but the molecular mechanisms remain unclear. Loach (Misgurnus anguillicaudatus) displays both sexual and polyploid growth dimorphism phenomena, and are therefore ideal models to study these two phenomena. In this study, RNA-seq was used for the first time to explore the differentially expressed genes (DEGs) between both sexes of diploid and tetraploid loaches in four tissues (brain, gonad, liver, and muscle). Results showed that 21,003, 17, and 1 DEGs were identified in gonad, liver, and muscle tissues, respectively, between females and males in both diploids and tetraploids. Regarding the ploidy levels, 4956, 1496, 2187, and 1726 DEGs were identified in the brain, gonad, liver, and muscle tissues, respectively, between tetraploids and diploids of the same sex. When both sexual and polyploid size dimorphisms were considered simultaneously in the four tissues, only 424 DEGs were found in the gonads, indicating that these gonadal DEGs may play an important regulatory role in regulating sexual and polyploid size dimorphisms. Regardless of the sex or ploidy comparison, the significant DEGs involved in glycolysis/gluconeogenesis and oxidative phosphorylation pathways were upregulated in faster-growing individuals, while steroid hormone biosynthesis-related genes and fatty acid degradation and elongation-related genes were downregulated. This suggests that fast-growing loaches (tetraploids, females) have higher energy metabolism levels and lower steroid hormone synthesis and fatty acid degradation abilities than slow-growing loaches (diploids, males). Our findings provide an archive for future systematic research on fish sexual and polyploid dimorphisms.

Dimorphism of sex and gonad-development-related genes in male and female lined seahorse, Hippocampus erectus, based on transcriptome analyses

Seahorse is characterized by its male pregnancy and sex-role reversal. To better understand the sexual dimorphism of male and female seahorses based on essential genes, we performed systematic transcriptome studies for both genders. A total of 157,834,590 cleaned reads were obtained and assembled into 129,268 transcripts and 31,764 could be annotated. Results showed that 176 up-regulated and 391 down-regulated transcripts were identified in the male seahorses compared with those in females. Genes involved in sex differentiation, such as dmrt1, sox9, fem1 and vasa, were identified and characterized. Moreover, the essential genes involved in reproductive molecular pathway were identified and analyzed in seahorses. In conclusion, the present study provides an archive for the future systematic research on seahorse sex differentiation.

Identification of genes involved in gonadal sex differentiation and the dimorphic expression pattern in Takifugu rubripes gonad at the early stage of sex differentiation

Quantifying the expression of mRNAs in the gonads at the critical stage of molecular sex differentiation stage might help to clarify the regulatory network during early sex differentiation and provide new information on the role of sex-related genes in gonadal function. In this study, transcriptomic analysis of sex-related genes expression profiles in fugu gonads at 60 and 90 days after hatching (dah) was conducted firstly, and a total of 112,504,991 clean reads, encompassing 28.35 Gb of sequences were retrieved. Twenty-three thousand eight hundred ten genes were found to be expressed in juvenile fugu gonads, and we mainly focused on the differentially expressed genes that have the potential to be involved in the gonadal sex differentiation. For 60-dah juveniles, we identified 1014 genes that were upregulated in the ovary and 1570 that were upregulated in the testis. For 90-dah juveniles, we identified 1287 genes that were upregulated in the ovary and 1500 that were upregulated in the testis. The dimorphic expression patterns of 15 genes in gonads at 30 and 40 dah were further investigate using qPCR. Cyp11b and star were expressed at higher levels in XY than in XX, while cyp11a1 and cyp19a1a were expressed at higher levels in XX than in XY at 30 dah. At 40 dah, the levels of gsdf, dmrt1, dmrt3, cyp11c1, star, and hsd3b expression were higher in XY, while the levels of foxl2, cyp19a1a, wnt9b, and foxD4 expression were higher in XX. Sox9, cyp11a1, cyp17a1, cyp17a2, and nr5a2 were expressed at similar levels in XX and XY at 40 dah. This is the first report of gonadal transcriptome of fugu at early sex differentiation stage, and our results provide an archive for further study on molecular mechanism underlying sex differentiation in this species.

Transcriptome characterization of HPG axis from Chinese sea perch Lateolabrax maculatus

Here the transcriptome and differential gene expression in the adult brain and gonads of the Chinese sea perch Lateolabrax maculatus were reported. A total of 78 256 909 clean reads were generated from the adult brain, ovary and testis by using the Illumina HiSeq2000 platform and assembled into 274 909 contigs. A total of 31 683 unigenes were annotated based on sequence similarity and 20 702 unigenes were found to exhibit 8237 gene ontology terms and 3888 signal pathways. Transcripts of 26 623 unigenes were present in all of the tissues, whereas pairwise comparisons revealed that 671/367, 496/315 and 1668/580 unigenes were up-down regulated by at least two-fold between the brain and ovary, ovary and testis and brain and testis, respectively. Homology search led to the identification of reproduction-associated genes of the brain-gonad axis, including those involved in sex differentiation and maintenance. The data provided an integrated and comprehensive transcriptome resource for L. maculatus, which could be used for further research on hypothalamus-pituitary-gonad axis gene function, reproduction regulation and sex-biased gene expression.

Widespread patterns of sexually dimorphic gene expression in an avian hypothalamic-pituitary-gonadal (HPG) axis

The hypothalamic-pituitary-gonadal (HPG) axis is a key biological system required for reproduction and associated sexual behaviors to occur. In the avian reproductive model of the rock dove (Columba livia), we characterized the transcript community of each tissue of the HPG axis in both sexes, thereby significantly expanding our mechanistic insight into HPG activity. We report greater sex-biased differential expression in the pituitary as compared to the hypothalamus, with multiple genes more highly expressed in the male pituitary being related to secretory function, and multiple genes more highly expressed in the female pituitary being related to reproduction, growth, and development. We report tissue-specific and sex-biased expression in genes commonly investigated when studying reproduction, highlighting the need for sex parity in future studies. In addition, we uncover new targets of investigation in both sexes, which could potentially change our understanding of HPG function.

De novo transcriptome sequencing and analysis of male, pseudo-male and female yellow perch, Perca flavescens

Transcriptome sequencing could facilitate discovery of sex-biased genes, biological pathways and molecular markers, which could help clarify the molecular mechanism of sex determination and sexual dimorphism, and assist with selective breeding in aquaculture. Yellow perch has unique gonad system and sexual dimorphism and is an alternative model to study mechanism of sex determination, sexual dimorphism and sexual selection. In this study, we performed the de novo assembly of yellow perch gonads and muscle transcriptomes by high throughput Illumina sequencing. A total of 212,180 contigs were obtained, ranging from 127 to 64,876 bp, and N50 of 1,066 bp. The assembly RNA-Seq contigs (≥200bp) were then used for subsequent analyses, including annotation, pathway analysis, and microsatellites discovery. No female- and pseudo-male-biased genes were involved in any pathways while male-biased genes were involved in 29 pathways, and neuroactive ligand receptor interaction and enzyme of trypsin (enzyme code, EC: 3.4.21.4) was highly involved. Pyruvate kinase (enzyme code, EC: 2.7.1.40), which plays important roles in cell proliferation, was highly expressed in muscles. In addition, a total of 183,939 SNPs, 11,286 InDels and 41,479 microsatellites were identified. This study is the first report on transcriptome information in Percids, and provides rich resources for conducting further studies on understanding the molecular basis of sex determinations, sexual dimorphism, and sexual selection in fish, and for population studies and marker-assisted selection in Percids.

The testis and ovary transcriptomes of the rock bream (Oplegnathus fasciatus): A bony fish with a unique neo Y chromosome

The rock bream (Oplegnathus fasciatus) is considerably one of the most economically important marine fish in East Asia and has a unique neo-Y chromosome system that is a good model to study the sex determination and differentiation in fish. In the present study, we used Illumina sequencing technology (HiSeq2000) to sequence, assemble and annotate the transcriptome of the testis and ovary tissues of rock bream. A total of 40,004,378 (NCBI SRA database SRX1406649) and 53,108,992 (NCBI SRA database SRX1406648) high quality reads were obtained from testis and ovary RNA sequencing, respectively, and 60,421 contigs (with average length of 1301 bp) were obtained after de novo assembling with Trinity software. Digital gene expression analysis reveals 14,036 contigs that show gender-enriched expressional profile with either testis-enriched (237 contigs) or ovary-enriched (581 contigs) with RPKM > 100. There are 237 male- and 582 female-abundant expressed genes that show sex dimorphic expression. We hope that the gonad transcriptome and those gender-enriched transcripts of rock bream can provide some insight into the understanding of genome-wide transcriptome profile of teleost gonad tissue and give useful information in fish gonad development.

Transcriptome Profile Analysis of Ovarian Tissues from Diploid and Tetraploid Loaches Misgurnus anguillicaudatus

RNA sequencing and short-read assembly was utilized to produce a transcriptome of ovarian tissues from three-year-old diploid and tetraploid loaches (Misgurnus anguillicaudatus). A total of 28,369 unigenes were obtained, comprising 10,546 unigenes with length longer than 1000 bp. More than 73% of the unigenes were annotated through sequence comparison with databases. The RNA-seq data revealed that 2253 genes were differentially expressed between diploid and tetraploid loaches, including 1263 up-regulated and 990 down-regulated genes in tetraploid loach. Some differentially expressed genes, such as vitellogenin (Vtg), gonadotropin releasing hormone receptor type A (GnRHRA), steroidogenic acute regulatory protein (StAR), mitogen-activated protein kinase 14a (MAPK14a), ATP synthase subunit alpha (atp5a), and synaptonemal complex protein 1 (Scp1), were involved in regulation of cell proliferation, division, gene transcription, ovarian development and energy metabolism, suggesting that these genes were related to egg diameter of the loach. Results of transcriptome profiling here were validated using real time quantitative PCR in ten selected genes. The present study provided insights into the transcriptome profile of ovarian tissues from diploid and tetraploid loaches Misgurnus anguillicaudatus, which was made available to the research community for functional genomics, comparative genomics, polyploidy evolution and molecular breeding of this loach and other related species.

Sequencing and De Novo Assembly of the Gonadal Transcriptome of the Endangered Chinese Sturgeon (Acipenser sinensis)

Background

The Chinese sturgeon (Acipenser sinensis) is endangered through anthropogenic activities including over-fishing, damming, shipping, and pollution. Controlled reproduction has been adopted and successfully conducted for conservation. However, little information is available on the reproductive regulation of the species. In this study, we conducted de novo transcriptome assembly of the gonad tissue to create a comprehensive dataset for A. sinensis.

Results

The Illumina sequencing platform was adopted to obtain 47,333,701 and 47,229,705 high quality reads from testis and ovary cDNA libraries generated from three-year-old A. sinensis. We identified 86,027 unigenes of which 30,268 were annotated in the NCBI non-redundant protein database and 28,281 were annotated in the Swiss-prot database. Among the annotated unigenes, 26,152 and 7,734 unigenes, respectively, were assigned to gene ontology categories and clusters of orthologous groups. In addition, 12,557 unigenes were mapped to 231 pathways in the Kyoto Encyclopedia of Genes and Genomes Pathway database. A total of 1,896 unigenes, potentially differentially expressed between the two gonad types, were found, with 1,894 predicted to be up-regulated in ovary and only two in testis. Fifty-five potential gametogenesis-related genes were screened in the transcriptome and 34 genes with significant matches were found. Besides, more paralogs of 11 genes in three gene families (sox, apolipoprotein and cyclin) were found in A. sinensis compared to their orthologs in the diploid Danio rerio. In addition, 12,151 putative simple sequence repeats (SSRs) were detected.

Conclusions

This study provides the first de novo transcriptome analysis currently available for A. sinensis. The transcriptomic data represents the fundamental resource for future research on the mechanism of early gametogenesis in sturgeons. The SSRs identified in this work will be valuable for assessment of genetic diversity of wild fish and genealogy management of cultured fish.

De novo characterization of fall dormant and nondormant alfalfa (Medicago sativa L.) leaf transcriptome and identification of candidate genes related to fall dormancy

Alfalfa (Medicago sativa L.) is one of the most widely cultivated perennial forage legumes worldwide. Fall dormancy is an adaptive character related to the biomass production and winter survival in alfalfa. The physiological, biochemical and molecular mechanisms causing fall dormancy and the related genes have not been well studied. In this study, we sequenced two standard varieties of alfalfa (dormant and non-dormant) at two time points and generated approximately 160 million high quality paired-end sequence reads using sequencing by synthesis (SBS) technology. The de novo transcriptome assembly generated a set of 192,875 transcripts with an average length of 856 bp representing about 165.1 Mb of the alfalfa leaf transcriptome. After assembly, 111,062 (57.6%) transcripts were annotated against the NCBI non-redundant database. A total of 30,165 (15.6%) transcripts were mapped to 323 Kyoto Encyclopedia of Genes and Genomes pathways. We also identified 41,973 simple sequence repeats, which can be used to generate markers for alfalfa, and 1,541 transcription factors were identified across 1,350 transcripts. Gene expression between dormant and non-dormant alfalfa at different time points were performed, and we identified several differentially expressed genes potentially related to fall dormancy. The Gene Ontology and pathways information were also identified. We sequenced and assembled the leaf transcriptome of alfalfa related to fall dormancy, and also identified some genes of interest involved in the fall dormancy mechanism. Thus, our research focused on studying fall dormancy in alfalfa through transcriptome sequencing. The sequencing and gene expression data generated in this study may be used further to elucidate the complete mechanisms governing fall dormancy in alfalfa.

Gonadal transcriptome analysis of male and female olive flounder (Paralichthys olivaceus)

Olive flounder (Paralichthys olivaceus) is an important commercially cultured marine flatfish in China, Korea, and Japan, of which female grows faster than male. In order to explore the molecular mechanism of flounder sex determination and development, we used RNA-seq technology to investigate transcriptomes of flounder gonads. This produced 22,253,217 and 19,777,841 qualified reads from ovary and testes, which were jointly assembled into 97,233 contigs. Among them, 23,223 contigs were mapped to known genes, of which 2,193 were predicted to be differentially expressed in ovary and 887 in testes. According to annotation information, several sex-related biological pathways including ovarian steroidogenesis and estrogen signaling pathways were firstly found in flounder. The dimorphic expression of overall sex-related genes provides further insights into sex determination and gonadal development. Our study also provides an archive for further studies of molecular mechanism of fish sex determination.