Comparative transcriptomics of spotted seatrout (Cynoscion nebulosus) populations to cold and heat stress

Transcriptomic changes behind Sparus aurata hepatic response to different aquaculture challenges: An RNA-seq study and multiomics integration

Gilthead seabream (Sparus aurata) is an important species in Mediterranean aquaculture. Rapid intensification of its production and sub-optimal husbandry practices can cause stress, impairing overall fish performance and raising issues related to sustainability, animal welfare, and food safety. The advent of next-generation sequencing technologies has greatly revolutionized the study of fish stress biology, allowing a deeper understanding of the molecular stress responses. Here, we characterized for the first time, using RNA-seq, the different hepatic transcriptome responses of gilthead seabream to common aquaculture challenges, namely overcrowding, net handling, and hypoxia, further integrating them with the liver proteome and metabolome responses. After reference-guided transcriptome assembly, annotation, and differential gene expression analysis, 7, 343, and 654 genes were differentially expressed (adjusted p-value < 0.01, log2|fold-change| >1) in the fish from the overcrowding, net handling, and hypoxia challenged groups, respectively. Gene set enrichment analysis (FDR < 0.05) suggested a scenario of challenge-specific responses, that is, net handling induced ribosomal assembly stress, whereas hypoxia induced DNA replication stress in gilthead seabream hepatocytes, consistent with proteomics and metabolomics' results. However, both responses converged upon the downregulation of insulin growth factor signalling and induction of endoplasmic reticulum stress. These results demonstrate the high phenotypic plasticity of this species and its differential responses to distinct challenging environments at the transcriptomic level. Furthermore, it provides significant resources for characterizing and identifying potentially novel genes that are important for gilthead seabream resilience and aquaculture production efficiency with regard to fish welfare.

Potential Effects of High Temperature and Heat Wave on Nanorana pleskei Based on Transcriptomic Analysis

In the context of climate change, understanding how indigenous amphibians of the Qinghai-Tibet plateau react to stresses and their coping mechanisms could be crucial for predicting their fate and successful conservation. A liver transcriptome for Nanorana pleskei was constructed using high-throughput RNA sequencing, and its gene expression was compared with frogs acclimated under either room temperature or high temperature and also heat wave exposed ones. A total of 126,465 unigenes were produced, with 66,924 (52.92%) of them being annotated. Up to 694 genes were found to be differently regulated as a result of abnormal temperature acclimatization. Notably, genes belonging to the heat shock protein (HSP) family were down-regulated in both treated groups. Long-term exposure to high-temperature stress may impair the metabolic rate of the frog and trigger the body to maintain a hypometabolic state in an effort to survive challenging times. During heat waves, unlike the high-temperature group, mitochondrial function was not impaired, and the energy supply was largely normal to support the highly energy-consuming metabolic processes. Genes were more transcriptionally suppressed when treated with high temperatures than heat waves, and the body stayed in low-energy states for combating these long-term adverse environments to survive. It might be strategic to preserve initiation to executive protein activity under heat wave stress. Under both stress conditions, compromising the protection of HSP and sluggish steroid activity occurred in frogs. Frogs were more affected by high temperatures than by heat waves.

Captivity induces large and population-dependent brain transcriptomic changes in wild-caught cane toads (Rhinella marina)

Gene expression levels are key molecular phenotypes at the interplay between genotype and environment. Mounting evidence suggests that short-term changes in environmental conditions, such as those encountered in captivity, can substantially affect gene expression levels. Yet, the exact magnitude of this effect, how general it is, and whether it results in parallel changes across populations are not well understood. Here, we take advantage of the well-studied cane toad, Rhinella marina, to examine the effect of short-term captivity on brain gene expression levels, and determine whether effects of captivity differ between long-colonized and vanguard populations of the cane toad's Australian invasion range. We compared the transcriptomes of wild-caught toads immediately assayed with those from toads captured from the same populations but maintained in captivity for seven months. We found large differences in gene expression levels between captive and wild-caught toads from the same population, with an over-representation of processes related to behaviour and the response to stress. Captivity had a much larger effect on both gene expression levels and gene expression variability in toads from vanguard populations compared to toads from long-colonized areas, potentially indicating an increased plasticity in toads at the leading edge of the invasion. Overall, our findings indicate that short-term captivity can induce large and population-specific transcriptomic changes, which has significant implications for studies comparing phenotypic traits of wild-caught organisms from different populations that have been held in captivity.

Min pig skeletal muscle response to cold stress

The increased sensitivity of pigs to ambient temperature is due to today's intensive farming. Frequent climate disasters increase the pressure on healthy pig farming. Min pigs are an indigenous pig breed in China with desirable cold resistance characteristics, and hence are ideal for obtaining cold-resistant pig breeds. Therefore, it is important to discover the molecular mechanisms that are activated in response to cold stress in the Min pig. Here, we conducted a transcriptomic analysis of the skeletal muscle of Min pigs under chronic low-temperature acclimation (group A) and acute short cold stress (group B). Cold exposure caused more genes to be upregulated. Totals of 125 and 96 differentially expressed genes (DEGs) were generated from groups A and B. Sixteen common upregulated DEGs were screened; these were concentrated in oxidative stress (SRXN1, MAFF), immune and inflammatory responses (ITPKC, AREG, MMP25, FOSL1), the nervous system (RETREG1, GADD45A, RCAN1), lipid metabolism (LRP11, LIPG, ITGA5, AMPD2), solute transport (SLC19A2, SLC28A1, SLCO4A1), and fertility (HBEGF). There were 102 and 73 genes that were specifically differentially expressed in groups A and B, respectively. The altered mRNAs were enriched in immune, endocrine, and cancer pathways. There were 186 and 91 differentially expressed lncRNAs generated from groups A and B. Analysis of the target genes suggested that they may be involved in regulating the MAPK signaling pathway for resistance to cold. The results of this study provide a comprehensive overview of cold exposure-induced transcriptional patterns in skeletal muscle of the Min pig. These results can guide future molecular studies of cold stress response in pigs for improving cold tolerance as a goal in breeding programs.

Comparative Transcriptomics of the Northern Quahog Mercenaria mercenaria and Southern Quahog Mercenaria campechiensis in Response to Chronic Heat Stress

The northern quahog (Mercenaria mercenaria) supports lucrative aquaculture industries in the USA. In the southeastern USA, aquacultured M. mercenaria faces increasing risks of summer die-offs from prolonged heat waves. We used a comparative transcriptomic approach to investigate the molecular responses of M. mercenaria and its southern congener, Mercenaria campechiensis, to controlled incremental heat stress over a 4-week period. Mercenaria were exposed to temperatures from 24 to 34 °C with 2.5 °C/week, after which, gill transcriptomes were de novo assembled and annotated. During the 4 weeks of chronic heat exposure, both species had the same survival rate (96%); M. mercenaria experienced body weight gain/loss depending on the originated hatcheries while M. campechiensis experienced an average net weight loss. The upregulated genes in both species included those in chaperone-mediated protein folding and regulation of cell death pathways, while the downregulated genes in both species involved in mRNA processing and splicing pathways. Compared to M. mercenaria, M. campechiensis appears to be more sensitive to prolonged heat stress as indicated by upregulating significantly more genes in coping with oxidative stress and in the protein degradation pathways, while downregulating some inhibitors of apoptosis. We discussed this finding within their ecological and evolutionary context. Our findings highlighted the potential vulnerability of the two quahogs, especially the southern quahog, to continued ocean warming.

Assembly-free rapid differential gene expression analysis in non-model organisms using DNA-protein alignment

BACKGROUND

RNA-seq is being increasingly adopted for gene expression studies in a panoply of non-model organisms, with applications spanning the fields of agriculture, aquaculture, ecology, and environment. For organisms that lack a well-annotated reference genome or transcriptome, a conventional RNA-seq data analysis workflow requires constructing a de-novo transcriptome assembly and annotating it against a high-confidence protein database. The assembly serves as a reference for read mapping, and the annotation is necessary for functional analysis of genes found to be differentially expressed. However, assembly is computationally expensive. It is also prone to errors that impact expression analysis, especially since sequencing depth is typically much lower for expression studies than for transcript discovery.

RESULTS

We propose a shortcut, in which we obtain counts for differential expression analysis by directly aligning RNA-seq reads to the high-confidence proteome that would have been otherwise used for annotation. By avoiding assembly, we drastically cut down computational costs - the running time on a typical dataset improves from the order of tens of hours to under half an hour, and the memory requirement is reduced from the order of tens of Gbytes to tens of Mbytes. We show through experiments on simulated and real data that our pipeline not only reduces computational costs, but has higher sensitivity and precision than a typical assembly-based pipeline. A Snakemake implementation of our workflow is available at: https://bitbucket.org/project_samar/samar .

CONCLUSIONS

The flip side of RNA-seq becoming accessible to even modestly resourced labs has been that the time, labor, and infrastructure cost of bioinformatics analysis has become a bottleneck. Assembly is one such resource-hungry process, and we show here that it can be avoided for quick and easy, yet more sensitive and precise, differential gene expression analysis in non-model organisms.