Sex-specific differences in transcriptome profiles of brain and muscle tissue of the tropical gar

Hepatic transcriptome, transcriptional effects and antioxidant responses in Poecilia vivipara exposed to sanitary sewage

Aquatic environments are subject to threats from multiple human activities, particularly through the release of untreated sanitary sewage into the coastal environments. These effluents contain a large group of natural or synthetic compounds referred to as emerging contaminants. Monitoring the types and quantities of toxic substances in the environment, especially complex mixtures, is an exhausting and challenging task. Integrative effect-based tools, such as biomarkers, are recommended for environmental quality monitoring programs. In this study, fish Poecilia vivipara were exposed for 24 and 96 h to raw untreated sewage diluted 33 % (v/v) in order to identify hepatic genes to be used as molecular biomarkers. Through a de novo hepatic transcriptome assembly, using Illumina MiSeq, 54,285 sequences were assembled creating a reference transcriptome for this guppy species. Transcripts involved in biotransformation systems, antioxidant defenses, ABC transporters, nuclear and xenobiotic receptors were identified and evaluated by qPCR. Sanitary sewage induced transcriptional changes in AhR, PXR, CYP2K1, CYP3A30, NQO1, UGT1A1, GSTa3, GSTmu, ST1C1, SOD, ABCC1 and SOX9 genes from liver of fish, particularly after 96 h of exposure. Changes in hepatic enzyme activities were also observed. The enzymes showed differences in fish exposed to both periods, while in the gills there was a prevalence of significant results after 96 h. The observed differences were associated to gender and/or to sewage exposure. The obtained results support the use of P. vivipara as sentinel and model organism for ecotoxicological studies and evidence the importance of understanding the differential responses associated to gender.

Navigating sex and sex roles: deciphering sex-biased gene expression in a species with sex-role reversal (Syngnathus typhle)

Sexual dimorphism, the divergence in morphological traits between males and females of the same species, is often accompanied by sex-biased gene expression. However, the majority of research has focused on species with conventional sex roles, where females have the highest energy burden with both egg production and parental care, neglecting the diversity of reproductive roles found in nature. We investigated sex-biased gene expression in Syngnathus typhle, a sex-role reversed species with male pregnancy, allowing us to separate two female traits: egg production and parental care. Using RNA sequencing, we examined gene expression across organs (brain, head kidney and gonads) at various life stages, encompassing differences in age, sex and reproductive status. While some gene groups were more strongly associated with sex roles, such as stress resistance and immune defence, others were driven by biological sex, such as energy and lipid storage regulation in an organ- and age-specific manner. By investigating how genes regulate and are regulated by changing reproductive roles and resource allocation in a model system with an unconventional life-history strategy, we aim to better understand the importance of sex and sex role in regulating gene expression patterns, broadening the scope of this discussion to encompass a wide range of organisms.

First neurotranscriptome of adults Tambaquis (Colossoma macropomum) with characterization and differential expression between males and females

The Tambaqui is one of the most representative Amazon fish species, being highly exploited in fisheries, aquaculture and as a research model. Nonetheless, data about functional genome are still required to evaluate reproductive and nutrition parameters as well as resistance to pathogens. The of next-generation sequencing has allows assessing the transcriptional processes in non-model species by providing comprehensive gene collections to be used as a database in further genomic applications and increased performance of captive populations. In this study, we relied on RNAseq approach to generate the first transcriptome of the telencephalon from adult males and females of Colossoma macropomum, resulting in a reference dataset for future functional studies. We retrieved 896,238 transcripts, including the identification of 267,785 contigs and 203,790 genes. From this total, 91 transcripts were differentially expressed, being 63 and 28 of them positively regulated for females and males, respectively. The functional annotation resulted in a library of 40 candidate genes for females and 20 for males. The functional enrichment classes comprised reproductive processes (GO:0,048,609; GO:0,003,006; GO:0,044,703; GO:0,032,504; GO:0,019,953) being related to sex differentiation (e.g., SAFB) and immune response (e.g., SLC2A6, AHNAK, NLRC3, NLRP3 and IgC MHC I alpha3), thus indicating that the genes in the neurotranscriptome of Tambaqui participate in sex differentiation and homeostasis of captive specimens. These data are useful to design the selection of genes related to sex determination and animal welfare in raising systems of Tambaqui.

Transcriptomic signatures associated with underlying rapid changes in the early phase brain of bi-directional sex change in Trimma okinawae

Teleost fish exhibit remarkable sexual plasticity and divergent developmental systems, including sequential hermaphroditism. One of the more fascinating models of sexual plasticity is socially controlled sex change, which is often observed in coral reef fish. The Okinawa rubble goby, Trimma okinawae, is a bi-directional sex-changing fish. It can rapidly change sex in either direction based on social circumstances. Although behavioural and neuroendocrine sex change occurs immediately and is believed to trigger gonadal changes, the underlying mechanisms remain poorly understood. In this study, we conducted a de novo transcriptome analysis of the T. okinawae brain and identified genes that are differentially expressed between the sexes and genes that were immediately controlled by social stimulation implicating sex change. Several genes showed concordant expression shifts regardless of the sex change direction and were associated with histone modification in nerve cells. These genes are known to function in the neuroendocrine control of reproduction in nerve cells. Overall, we identified genes associated with the initiation of sex change, which provides insight into the regulation of sex change and sexual plasticity.

Identification and expression analysis of transcripts involved in taurine biosynthesis during early ontogeny of tropical gar Atractosteus tropicus

In fishes, the availability of taurine is regulated during ontogenetic development, where its endogenous synthesis capacity is species dependent. Thus, different pathways and involved enzymes have been described: pathway I (cysteine sulfinate-dependent pathway), cysteine dioxygenase type 1 (cdo1) and cysteine sulfinic acid decarboxylase (csad); pathway II (cysteic acid pathway), cdo1 and glutamic acid decarboxylase (gad); and pathway III (cysteamine pathway), 2-aminoethanethiol dioxygenase (ado); whereas taurine transporter (taut) is responsible for taurine entry into cells on the cell membrane and the mitochondria. This study determined if the tropical gar (Atractosteus tropicus), an ancient holostean fish model, has the molecular mechanism to synthesize taurine through the identification and analysis expression of transcripts coding for proteins involved in its biosynthesis and transportation, at different embryo-larvae stages and in different organs of juveniles (31 dah). We observed a fluctuating expression of all transcripts involved in the three pathways at all analyzed stages. All transcripts are expressed during the beginning of larval development; however, ado and taut show a peak expression at 9 dah, and all transcripts but csad decreased at 23 dah, when the organism ended the larval period. Furthermore, at 31 dah, we observed taut expression in all examined organs. The transcripts involved in pathways I and III are expressed differently across all organs, whereas pathway II was only observed in the brain, eye, and skin. The results suggested that taurine biosynthesis in tropical gar is regulated during its early development before first feeding, and the pathway might also be organ-type dependent.

DNA Methylation Difference between Female and Male Ussuri Catfish (Pseudobagrus ussuriensis) in Brain and Gonad Tissues

DNA methylation has been found to be involved in sex determination and differentiation in many aquaculture species. The Ussuri catfish (Pseudobagrus ussuriensis) is a popular aquaculture fish in China with high economic value in which male-biased sex dimorphism was observed in terms of body size and body weight. In this study, DNA methylation-sensitive RAD sequencing (Methyl-RAD) was used to explore the epigenetic difference between adult male and female samples in brain and gonad tissues. In brain tissues, 5,442,496 methylated cytosine sites were found and 9.94% of these sites were from symmetric CCGG or CCWGG sites. Among these sites, 321 differential DNA methylation sites (DMSs) in 171 genes were identified, while in gonad tissues, 4,043,053 methylated cytosines sites were found in total and 11.70% of them were from CCGG or CCWGG. Among these sites, 78 differential DNA methylation sites were found which were located in 64 genes. We also found several sex-determination genes among these differential methylated genes, such as amh, gsdf and hsd11b2 in brain tissues and slco3a1, socs2 and trim47 in gonad tissues. These results provided evidence for understanding the function of DNA methylation in the sex differentiation in Pseudobagrus ussuriensis, which further deepens the relationship between gene regulation and epigenetics.

Sex-biased dynamics of three-spined stickleback (Gasterosteus aculeatus) gene expression patterns

The study of the differences between sexes presents an excellent model to unravel how phenotypic variation is achieved from a similar genetic background. Sticklebacks are of particular interest since evidence of a heteromorphic chromosome pair has not always been detected. The present study investigated sex-biased mRNA and small non-coding RNA (sncRNA) expression patterns in the brain, adipose tissues, and gonads of the three-spined stickleback. The sncRNA analysis indicated that regulatory functions occurred mainly in the gonads. Alleged miRNA-mRNA interactions were established and a mapping bias of differential expressed transcripts towards chromosome 19 was observed. Key players previously shown to control sex determination and differentiation in other fish species but also genes like gapdh were among the transcripts identified. This is the first report in the three-spined stickleback demonstrating tissue-specific expression comprising both mRNA and sncRNA between sexes, emphasizing the importance of mRNA-miRNA interactions as well as new presumed genes not yet identified to have gender-specific roles.

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

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

Transcriptomics of Tasmanian Devil (Sarcophilus Harrisii) Ear Tissue Reveals Homogeneous Gene Expression Patterns across a Heterogeneous Landscape

In an era of unprecedented global change, exploring patterns of gene expression among wild populations across their geographic range is crucial for characterizing adaptive potential. RNA-sequencing studies have successfully characterized gene expression differences among populations experiencing divergent environmental conditions in a wide variety of taxa. However, few of these studies have identified transcriptomic signatures to multivariate, environmental stimuli among populations in their natural environments. Herein, we aim to identify environmental and sex-driven patterns of gene expression in the Tasmanian devil (Sarcophilus harrisii), a critically endangered species that occupies a heterogeneous environment. We performed RNA-sequencing on ear tissue biopsies from adult male and female devils from three populations at the extremes of their geographic range. There were no transcriptome-wide patterns of differential gene expression that would be suggestive of significant, environmentally-driven transcriptomic responses. The general lack of transcriptome-wide variation in gene expression levels across the devil's geographic range is consistent with previous studies that documented low levels of genetic variation in the species. However, genes previously implicated in local adaptation to abiotic environment in devils were enriched for differentially expressed genes. Additionally, three modules of co-expressed genes were significantly associated with either population of origin or sex.

Integrated analysis of mRNA and miRNA expression profiles reveals muscle growth differences between adult female and male Chinese concave-eared frogs (Odorrana tormota)

The Chinese concave-eared torrent frog (Odorrana tormota) is the first known non-mammalian vertebrate that can communicate using ultrasound. In this species, females are approximately four times as large as males, in which the female growth rate is obviously higher than that of male. Until now, the molecular mechanisms underlying muscle growth development differences between male and female frogs have not been reported. Here, we integrated mRNA and miRNA expression profiles to reveal growth differences in the hindlimb muscles of 2-year-old frogs. Among 569 differentially expressed genes (DEGs), 69 were associated with muscle growth and regeneration. Fifty-one up-regulated genes in females were potentially involved in promoting muscle growth and regeneration, whereas 18 up-regulated genes in males may lead to muscle growth inhibition and fast-twitch muscle fiber contraction. 244 DEGs were enriched in mTOR and other protein synthesis signaling pathways, and protein degradation pathways, including lysosomal protease, calpain, caspase, and ubiquitin-proteasome system pathways. It may interpret why female muscles grow faster than males. Based on expression differences of genes involved in glycolysis and oxidative metabolism, we speculated that the proportion of slow muscle fiber was higher and that of fast muscle fiber was lower in female compared with male muscle. Additionally, 767 miRNAs were identified, including 217 new miRNAs, and 6248 miRNA-negatively regulated mRNAs were predicted. The miRNA target genes were enriched in pathways related to muscle growth, protein synthesis, and degradation. Thus, in addition to the identified mRNA differential expressions, miRNAs may play other important roles in the differential regulation of hindlimb muscle growth between female and male O. tormota.

Transcriptome analysis for identification of candidate genes related to sex determination and growth in Charybdis japonica

Charybdis japonica is an important cultured crab in China and it exhibits sex differences in their growth. Growth is an important economic trait that is controlled by many genes. In order to discover the growth-related regulatory mechanisms, whole-body transcriptomic sequencing and comparative analyses in different genders of C. japonica were conducted based on Illumina RNA-seq technology. In the present study, we obtained 62,879,204 and 60,226,334 clean reads in female and male libraries, respectively. 25,000,000 clean reads of every library were randomly selected and compared with Nt database to examine the possible contamination. Results showed that all clean reads were distributed among C. japonica or other species that were closely relative to this species, indicating no-pollution. De novo assembly was performed and a total of 32,543 and 44,174 unigenes were produced in female and male of C. japonica, respectively. Among all the unigenes, 12,591 and 14,455 unigenes of female and male crabs were annotated based on protein databases. Moreover, a total of 33,926 unigenes were found to contain ORFs and 52,839 SSRs were detected. The contrast between male and female C. japonica identifying 1939 unigenes were significantly differentially expressed. In addition, we specifically discussed some gene functions and pathways potentially associated with sex determination and growth. This is the first systematic report of whole transcriptome in C. japonica. The transcriptome information provides a basic resource for further studies on understanding the molecular basis of biological processes in C. japonica and other crustaceans.