Evaluating the Hypoxia Response of Ruffe and Flounder Gills by a Combined Proteome and Transcriptome Approach

Developing molecular classifiers to detect environmental stressors, smolt stages and morbidity in coho salmon, Oncorhynchus kisutch

Aquatic species are increasingly confronted with environmental stressors because of climate change. Although molecular technologies have advanced our understanding of how organisms respond to stressors in laboratory settings, the ability to detect physiological responses to specific stressors under complex field conditions remains underdeveloped. This research applied multi-stressor challenge trials on coho salmon, employing the "Salmon Fit-Chips" genomic tool and a random forest-based classification model to develop classifiers predictive for chronic thermal and hypoxic stress, as well as salinity acclimation, smolt stage and morbidity status. The study also examined how smolts and de-smolts (smolts not having entered SW during the smolt window) responded transcriptionally to exposure to saltwater. Using RF classifiers optimized with 4 to 12 biomarkers, we identified transcriptional signatures that accurately predicted the presence of each stressor and physiological state, achieving prediction accuracy rates between 86.8 % and 100 %, regardless of other background stressors present. Stressor recovery time was established by placing fish back into non-stressor conditions after stress exposure, providing important context to stressor detections in field applications. Recovery from thermal and hypoxic stress requires about 3 and 2 days, respectively, with >3 days needed for re-acclimation to freshwater for seawater acclimated fish. The study also found non-additive (synergistic) effects of multiple stressors on mortality risk. Importantly, osmotic stress associated with de-smolts was the most important predictor of mortality. In saltwater, de-smolts exposed to salinity, high temperature, and hypoxia experienced a 9-fold increase in mortality compared to those only exposed to saltwater, suggesting a synergistic response to multiple stressors. These findings suggest that delays in hatchery releases to support release of larger fish need to be carefully scrutinized to ensure fish are not being released as de-smolts, which are highly susceptible to additional climate-induced stressors like rising temperatures and reduced dissolved oxygen levels in the marine environment.

Network-based integration of omics, physiological and environmental data in real-world Elbe estuarine Zander

Coastal and estuarine environments are under endogenic and exogenic pressures jeopardizing survival and diversity of inhabiting biota. Information of possible synergistic effects of multiple (a)biotic stressors and holobiont interaction are largely missing in estuaries like the Elbe but are of importance to estimate unforeseen effects on animals' physiology. Here, we seek to leverage host-transcriptional RNA-seq and gill mucus microbial 16S rRNA metabarcoding data coupled with physiological and abiotic measurements in a network analysis approach to decipher the impact of multiple stressors on the health of juvenile Sander lucioperca along one of the largest European estuaries. We find mesohaline areas characterized by gill tissue specific transcriptional responses matching osmosensing and tissue remodeling. Liver transcriptomes instead emphasized that zander from highly turbid areas were undergoing starvation which was supported by compromised body condition. Potential pathogenic bacteria, including Shewanella, Acinetobacter, Aeromonas and Chryseobacterium, dominated the gill microbiome along the freshwater transition and oxygen minimum zone. Their occurrence coincided with a strong adaptive and innate transcriptional immune response in host gill and enhanced energy demand in liver tissue supporting their potential pathogenicity. Taken together, we show physiological responses of a fish species and its microbiome to abiotic factors whose impact is expected to increase with consequences of climate change. We further present a method for the close-meshed detection of the main stressors and bacterial species with disease potential in a highly productive ecosystem.

Are the effects of catch-and-release angling evident in changes to mRNA abundances related to metabolism, acid-base regulation and stress in lake trout (Salvelinus namaycush) gills?

Catch-and-release (C&R) angling is a conservation-oriented practice intended to reduce the impact recreational angling has on fish populations. Even though most recreationally angled fish are released, little is known about how C&R angling impacts fish at the cellular or tissue level. As the first to explore the impacts of C&R angling on mRNA abundances, our study aimed to identify how the stress of angling influenced metabolism, acid-base regulation and cellular stress in the gills of lake trout (Salvelinus namaycush). Because gills are responsible for metabolic gas exchange, are crucial sites of acid-base homeostasis and respond to stressors quickly, we hypothesized that the relative mRNA abundance of genes related to these three physiological processes would be altered after angling. We took gill samples of live lake trout at 0, 2 or 48 h after fish were angled by rod and reel, and then used quantitative PCR (qPCR) to measure the relative abundance of nine candidate mRNA transcripts. Heat shock protein 70 (hsp70) mRNA levels significantly increased over 5-fold 2 h after angling, indicating a potential activation of a cytoprotective response. However, contrary to our hypothesis, we observed no change in the relative mRNA abundance of genes related to metabolism or acid-base regulation in response to C&R angling within a 48-h period. As C&R angling can negatively impact fish populations, further use of transcript-level studies will allow us to understand the impact C&R has on specific tissues and improve our knowledge of how C&R influences overall fish health.

Transcriptomic and proteomic analyses reveal the common and unique pathway(s) underlying different skin colors of leopard coral grouper (Plectropomus leopardus)

To gain a comprehensive and unbiased molecular understanding of the different skin colors of P. leopardus, we used Illumina HiSeq 2500 and TMT (Tandem Mass Tag) to compare transcription and protein levels between red and black skin of P. leopardus. We identified 797 upregulated and 314 downregulated genes (differentially expressed genes; DEGs) in red (RG) compared with black (BG) skin of P. leopardus. We also identified 377 differentially abundant proteins (DAPs), including 314 upregulated and 63 downregulated proteins. These DEGs and DAPs were significantly enriched in melanin synthesis (e.g., pyrimidine metabolism, Phenylalanine, tyrosine, and tryptophan biosynthesis, melanogenesis, phenylalanine metabolism, and tyrosine metabolism), oxidative phosphorylation (e.g., phosphonate and phosphinate metabolism, and oxidative phosphorylation), energy metabolism (e.g., HIF-1, glycolysis/gluconeogenesis, fatty acid biosynthesis, and fatty acid degradation), and signal transduction (e.g., Wnt, calcium, MAPK, and cGMP-PKG signaling pathways), etc. Further analysis of MAPKs showed that the activation levels of its main members JNK1 and ERK1/2 differed significantly between red and black skin colors. After RNAi was used to interfere with ERK1/2, it was found that the local skin of the tail of P. leopardus would turn black. Combined transcriptome and proteome analysis showed that most DEGs-DAPs in red skin were higher than in black skin (58 were upregulated, 1 was downregulated, and 4 were opposite). These DEGs-DAPs showed that the differences between red and black skin tissues of P. leopardus were related primarily to energy metabolism, signal transduction and cytoskeleton. These findings are not only conducive to understand the skin color regulation mechanism of P. leopardus and other coral reef fish, but also provide an important descriptive to the breeding of color strains. SIGNIFICANCE OF THE STUDY: The skin color of P. leopardus gradually darkens or blackens due to environmental factors such as changes in light intensity and human activities, and this directly affects its ornamental and economic value. In this study, RNAseq and TMT were used to conduct comparative quantitative transcriptomics and proteomics and analyze differences between red and black P. leopardus skin. The results showed that energy metabolism, signal transduction and cytoskeleton were the main metabolic pathways causing their skin color differences. These findings contribute to existing data describing fish skin color, and provide information about protein levels, which are of great significance to a deeper understanding of the skin color regulation mechanism in P. leopardus and other coral reef fishes.

Identification of Hypoxia-Specific Biomarkers in Salmonids Using RNA-Sequencing and Validation Using High-Throughput qPCR

Identifying early gene expression responses to hypoxia (i.e., low dissolved oxygen) as a tool to assess the degree of exposure to this stressor is crucial for salmonids, because they are increasingly exposed to hypoxic stress due to anthropogenic habitat change, e.g., global warming, excessive nutrient loading, and persistent algal blooms. Our goal was to discover and validate gill gene expression biomarkers specific to the hypoxia response in salmonids across multi-stressor conditions. Gill tissue was collected from 24 freshwater juvenile Chinook salmon (Oncorhynchus tshawytscha), held in normoxia [dissolved oxygen (DO) > 8 mg L^-1] and hypoxia (DO = 4‒5 mg L^-1) in 10 and 18° temperatures for up to six days. RNA-sequencing (RNA-seq) was then used to discover 240 differentially expressed genes between hypoxic and normoxic conditions, but not affected by temperature. The most significantly differentially expressed genes had functional roles in the cell cycle and suppression of cell proliferation associated with hypoxic conditions. The most significant genes (n = 30) were selected for real-time qPCR assay development. These assays demonstrated a strong correlation (r = 0.88; P < 0.001) between the expression values from RNA-seq and the fold changes from qPCR. Further, qPCR of the 30 candidate hypoxia biomarkers was applied to an additional 322 Chinook salmon exposed to hypoxic and normoxic conditions to reveal the top biomarkers to define hypoxic stress. Multivariate analyses revealed that smolt stage, water salinity, and morbidity status were relevant factors to consider with the expression of these genes in relation to hypoxic stress. These hypoxia candidate genes will be put into application screening Chinook salmon to determine the identity of stressors impacting the fish.

Transcriptome analysis reveals molecular strategies in gills and heart of large yellow croaker (Larimichthys crocea) under hypoxia stress

The gills and heart are two major targets of hypoxia in fish. However, the molecular responses in fish gills and heart to hypoxia challenge remain unclear. Here, RNA-Seq technology was used to study the gene expression profiles in gills and heart of large yellow croaker (Larimichthys crocea) at 6, 24, and 48 h after hypoxia stress. A total of 1,546 and 2,746 differentially expressed genes (DEGs) were identified in gills and heart, respectively. Expression changes of nine genes in each tissue were further validated by the qPCR. Based on KEGG and Gene ontology enrichments, we found that various innate immunity-related genes, such as complement components (C1qs, C2, C3, C6, and C7), chemokines (CCL3, CCL17, CCL19, CCL25, and CXCL8_L3), chemokine receptors (CCR9, CXCR1, and CXCR3), and nitric oxide synthase (NOS), were significantly down-regulated in gills and/or heart, suggesting that innate immune processes mediated by these genes may be inhibited by hypoxia. The genes involved in both glycolysis pathway (LDHA) and tricarboxylic acid cycle (IDH2 and OGDH) were up-regulated in gills and heart of hypoxic large yellow croakers, possibly because gill and heart tissues need enough energy to accelerate gas exchange and blood circulation. Hypoxia also affected the ion transport in gills of large yellow croaker, through down-regulating the expression levels of numerous classical ion transporters, including HVCN1, SLC20A2, SLC4A4, RHBG, RHCG, and SCN4A, suggesting an energy conservation strategy to hypoxia stress. All these results indicate that the immune processes, glycolytic pathways, and ion transport were significantly altered in gills and/or heart of large yellow croaker under hypoxia, possibly contributing to maintain cellular energy balance during hypoxia. Our data, therefore, afford new information to understand the tissue-specific molecular responses of bony fish to hypoxia stress.

De novo transcriptome assembly, annotation and comparison of four ecological and evolutionary model salmonid fish species

Background

Salmonid fishes exhibit high levels of phenotypic and ecological variation and are thus ideal model systems for studying evolutionary processes of adaptive divergence and speciation. Furthermore, salmonids are of major interest in fisheries, aquaculture, and conservation research. Improving understanding of the genetic mechanisms underlying traits in these species would significantly progress research in these fields. Here we generate high quality de novo transcriptomes for four salmonid species: Atlantic salmon (Salmo salar), brown trout (Salmo trutta), Arctic charr (Salvelinus alpinus), and European whitefish (Coregonus lavaretus). All species except Atlantic salmon have no reference genome publicly available and few if any genomic studies to date.

Results

We used paired-end RNA-seq on Illumina to generate high coverage sequencing of multiple individuals, yielding between 180 and 210 M reads per species. After initial assembly, strict filtering was used to remove duplicated, redundant, and low confidence transcripts. The final assemblies consisted of 36,505 protein-coding transcripts for Atlantic salmon, 35,736 for brown trout, 33,126 for Arctic charr, and 33,697 for European whitefish and are made publicly available. Assembly completeness was assessed using three approaches, all of which supported high quality of the assemblies: 1) ~78% of Actinopterygian single-copy orthologs were successfully captured in our assemblies, 2) orthogroup inference identified high overlap in the protein sequences present across all four species (40% shared across all four and 84% shared by at least two), and 3) comparison with the published Atlantic salmon genome suggests that our assemblies represent well covered (~98%) protein-coding transcriptomes. Thorough comparison of the generated assemblies found that 84-90% of transcripts in each assembly were orthologous with at least one of the other three species. We also identified 34-37% of transcripts in each assembly as paralogs. We further compare completeness and annotation statistics of our new assemblies to available related species.

Conclusion

New, high-confidence protein-coding transcriptomes were generated for four ecologically and economically important species of salmonids. This offers a high quality pipeline for such complex genomes, represents a valuable contribution to the existing genomic resources for these species and provides robust tools for future investigation of gene expression and sequence evolution in these and other salmonid species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4379-x) contains supplementary material, which is available to authorized users.

The expression of hypoxia-inducible factor-1α gene is not affected by low-oxygen conditions in yellow perch (Perca flavescens) juveniles

Hypoxia can affect various fish populations, including yellow perch Perca flavescens, which is an economically and ecologically important species in Lake Erie, a freshwater system that often experiences hypoxia in the hypolimnetic part of the lake. Fish, similarly to mammals, possess molecular oxygen sensor-hypoxia-inducible factor-1 (HIF-1), a transcription factor that can affect expression of many downstream genes related to animal growth and locomotion, protein synthesis, as well as ATP and amino acid metabolism. HIF-1 is a heterodimer, which consists of two subunits: oxygen-sensitive and oxygen-insensitive subunits, α and β, respectively. In this study, we report first on the molecular cloning and sequencing of P. flavescens HIF-1α. The full-length complementary DNA (cDNA) was isolated and submitted to the GenBank with accession number KT783483. It consists of 3529 base pairs (bp) carrying a single open-reading frame that encompasses 2250 bp of the coding region, 247 bp of the 5' untranslated region (UTR), and 1032 bp of the 3' UTR. The "de novo" prediction of the 3D structure of HIF-1α protein, which consists of 749 amino acids, is presented, too. We then utilized One-Step Taqman® real-time RT-PCR technology to monitor changes in HIF-1α messenger RNA (mRNA) copies in response to chronic hypoxic stress. An experiment was conducted using 14-day post-swim-up stage yellow perch larvae with uninflated swim bladders. This experiment included three treatment groups: hypoxia, mid-hypoxia, and normoxia, in four replicates (four tanks per treatment) with the following dissolved oxygen levels: 3, 4, and >7 mg O₂/L, respectively. At the end (2 weeks) and in the middle (1 week) of the experiment, fish from each tank were sampled for body measurements and molecular biology analysis. The results showed no differences in survival (∼90%) between treatment groups. Oxygen concentration was lowered to 3.02 ± 0.15 (mean ± SE) mg O₂/L with no adverse effect on fish survival. The highest growth rate was observed in the normoxic group. A similar trend was observed with fish body length. The growth rate of fish declined with decreasing water-dissolved oxygen. The number of HIF-1α mRNA copies was not significantly different between hypoxic, mid-hypoxic, and normoxic conditions, and this was true for fish obtained in the middle and at the end of the experiment. Graphical abstract.

Proteomic responses of European flounder to temperature and hypoxia as interacting stressors: Differential sensitivities of populations

In the context of global change, ectotherms are increasingly impacted by abiotic perturbations. Along the distribution area of a species, the populations at low latitudes are particularly exposed to temperature increase and hypoxic events. In this study, we have compared the proteomic responses in the liver of European flounder populations, by using 2-D electrophoresis. One southern peripheral population from Portugal vs two northern core populations from France, were reared in a common garden experiment. Most of the proteomic differences were observed between the two experimental conditions, a cold vs a warm and hypoxic conditions. Consistent differentiations between populations were observed in accumulation of proteins involved in the bioenergetics- and methionine-metabolisms, fatty acids transport, and amino-acid catabolism. The specific regulation of crucial enzymes like ATP-synthase and G6PDH, in the liver of the southern population, could be related to a possible local adaptation. This southern peripheral population is spatially distant from northern core populations and has experienced dissimilar ecological conditions; thus it may contain genotypes that confer resilience to climate changes.