Transcriptome wide analyses reveal a sustained cellular stress response in the gill tissue of Trematomus bernacchii after acclimation to multiple stressors

The hypoxia response pathway in the Antarctic fish Notothenia coriiceps is functional despite a poly Q/E insertion mutation in HIF-1α

Antarctic notothenioid fishes, inhabiting the oxygen-rich Southern Ocean, possess a polyglutamine and glutamic acid (poly Q/E) insertion mutation in the master transcriptional regulator of oxygen homeostasis, hypoxia- inducible factor-1α (HIF-1α). To determine if this mutation impairs the ability of HIF-1 to regulate gene expression in response to hypoxia, we exposed Notothenia coriiceps, with a poly Q/E insertion mutation in HIF-1α that is 9 amino acids long, to hypoxia (2.3 mg L-1 O2) or normoxia (10 mg L -1 O2) for 12 h. Heart ventricles, brain, liver, and gill tissue were harvested and changes in gene expression quantified using RNA sequencing. Levels of glycogen and lactate were also quantified to determine if anaerobic metabolism increases in response to hypoxia. Exposure to hypoxia resulted in 818 unique differentially expressed genes (DEGs) in liver tissue of N. coriiceps. Many hypoxic genes were induced, including ones involved in the MAP kinase and FoxO pathways, glycolytic metabolism, and vascular remodeling. In contrast, there were fewer than 104 unique DEGs in each of the other tissues sampled. Lactate levels significantly increased in liver in response to hypoxia, indicating that anaerobic metabolism increases in response to hypoxia in this tissue. Overall, our results indicate that the hypoxia response pathway is functional in N. coriiceps despite a poly Q/E mutation in HIF-1α, and confirm that Antarctic fishes are capable of altering gene expression in response to hypoxia.

Transcriptome profiles revealed high- and low-salinity water altered gill homeostasis in half-smooth tongue sole (Cynoglossus semilaevis)

Salinity is an important environmental factor that affects fish growth, development, and reproduction. As euryhaline fish, half-smooth tongue sole (Cynoglossus semilaevis) are a suitable species for deciphering the salinity adaptation mechanism of fish; however, the molecular mechanisms underlying low- and high-salinity responses remain unclear. In this study, RNA-seq was applied to characterize the genes and regulatory pathways involved in C. semilaevis gill responses to high- (32 ppt), low- (8 ppt), and control-salinity (24 ppt) water. Gills were rich in mitochondria-rich cells (MRCs) in high salinity. Compared with control, 2137 and 218 differentially expressed genes (DEGs) were identified in low and high salinity, respectively. The enriched functions of most DEGs were metabolism, ion transport, regulation of cell cycle, and immune response. The DEGs involved in oxidative phosphorylation, citrate cycle, and fatty acid metabolism were down-regulated in low salinity. For ion transport, high and low salinity significantly altered the expressions of prlr, ca12, and cftr. In cell cycle arrest and cellular repair, gadd45b, igfbp5, and igfbp2 were significantly upregulated in high and low salinity. For immune response, il10, il34, il12b, and crp increased in high and low salinity. Our findings suggested that alterations in material and energy metabolism, ions transport, cell cycle arrest, cellular repair, and immune response, are required to maintain C. semilaevis gill homeostasis under high and low salinity. This study provides insight into the divergence of C. semilaevis osmoregulation mechanisms acclimating to high and low salinity, which will serve as reference for the healthy culture of C. semilaevis.

Can heat shock protein 70 (HSP70) serve as biomarkers in Antarctica for future ocean acidification, warming and salinity stress?

The Antarctic Peninsula is one of the fastest-warming places on Earth. Elevated sea water temperatures cause glacier and sea ice melting. When icebergs melt into the ocean, it “freshens” the saltwater around them, reducing its salinity. The oceans absorb excess anthropogenic carbon dioxide (CO2) causing decline in ocean pH, a process known as ocean acidification. Many marine organisms are specifically affected by ocean warming, freshening and acidification. Due to the sensitivity of Antarctica to global warming, using biomarkers is the best way for scientists to predict more accurately future climate change and provide useful information or ecological risk assessments. The 70-kilodalton (kDa) heat shock protein (HSP70) chaperones have been used as biomarkers of stress in temperate and tropical environments. The induction of the HSP70 genes (Hsp70) that alter intracellular proteins in living organisms is a signal triggered by environmental temperature changes. Induction of Hsp70 has been observed both in eukaryotes and in prokaryotes as response to environmental stressors including increased and decreased temperature, salinity, pH and the combined effects of changes in temperature, acidification and salinity stress. Generally, HSP70s play critical roles in numerous complex processes of metabolism; their synthesis can usually be increased or decreased during stressful conditions. However, there is a question as to whether HSP70s may serve as excellent biomarkers in the Antarctic considering the long residence time of Antarctic organisms in a cold polar environment which appears to have greatly modified the response of heat responding transcriptional systems. This review provides insight into the vital roles of HSP70 that make them ideal candidates as biomarkers for identifying resistance and resilience in response to abiotic stressors associated with climate change, which are the effects of ocean warming, freshening and acidification in Antarctic organisms.

Morphological analysis of erythrocytes of an Antarctic teleost under heat stress: Bias of the stabling effect

The stenothermal Antarctic fish that live in the coastal waters of the Terra Nova Bay (Ross Sea) are rarely exposed to temperatures above zero during the year. We tested whether a slight temperature rise of 1.5 °C affects a sensitive biomarker such as erythrocytes morphology in sections of blood pellets of a small demersal notothen. The erythrocytes' shape descriptors showed significant or highly significant differences temporally from the capture of fish to the conclusion of the experiment. Surprisingly, the erythrocyte's morphology did not show significant differences between the two experimental conditions, returning similar results in control fish stabled at -0.9 °C and in the fish treated at +0.6 °C, although the values of the shape descriptors were often lower in the latter. This study demonstrates the critical issues of comparative physiology in the study of extremely sensitive organisms, such as the fish of the High Antarctic Zone. Moreover, the stabling effect inside the aquarium facilities appears to significantly obscure the effects of the experimental heat treatment.

Combined effects of ocean warming and acidification on marine fish and shellfish: A molecule to ecosystem perspective

It is expected that by 2050 human population will exceed nine billion leading to increased pressure on marine ecosystems. Therefore, it is conjectured various levels of ecosystem functioning starting from individual to population-level, species distribution, food webs and trophic interaction dynamics will be severely jeopardized in coming decades. Ocean warming and acidification are two prime threats to marine biota, yet studies about their cumulative effect on marine fish and shellfishes are still in its infancy. This review assesses existing information regarding the interactive effects of global environmental factors like warming and acidification in the perspective of marine capture fisheries and aquaculture industry. As climate change continues, distribution pattern of species is likely to be altered which will impact fisheries and fishing patterns. Our work is an attempt to compile the existing literatures in the biological perspective of the above-mentioned stressors and accentuate a clear outline of knowledge in this subject. We reviewed studies deciphering the biological consequences of warming and acidification on fish and shellfishes in the light of a molecule to ecosystem perspective. Here, for the first time impacts of these two global environmental drivers are discussed in a holistic manner taking into account growth, survival, behavioural response, prey predator dynamics, calcification, biomineralization, reproduction, physiology, thermal tolerance, molecular level responses as well as immune system and disease susceptibility. We suggest urgent focus on more robust, long term, comprehensive and ecologically realistic studies that will significantly contribute to the understanding of organism's response to climate change for sustainable capture fisheries and aquaculture.

Adaptations and Diversity of Antarctic Fishes: A Genomic Perspective

Antarctic notothenioid fishes are the classic example of vertebrate adaptive radiation in a marine environment. Notothenioids diversified from a single common ancestor ∼25 Mya to more than 140 species today, and they represent ∼90% of fish biomass on the continental shelf of Antarctica. As they diversified in the cold Southern Ocean, notothenioids evolved numerous traits, including osteopenia, anemia, cardiomegaly, dyslipidemia, and aglomerular kidneys, that are beneficial or tolerated in their environment but are pathological in humans. Thus, notothenioids are models for understanding adaptive radiations, physiological and biochemical adaptations to extreme environments, and genetic mechanisms of human disease. Since 2014, 16 notothenioid genomes have been published, which enable a first-pass holistic analysis of the notothenioid radiation and the genetic underpinnings of novel notothenioid traits. Here, we review the notothenioid radiation from a genomic perspective and integrate our insights with recent observations from other fish radiations. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Gill Transcriptome Analysis Revealed the Difference in Gene Expression Between Freshwater and Seawater Acclimated Guppy (Poecilia reticulata)

Guppy (Poecilia reticulata) can adapt to a wide range of salinity changes. To investigate the gene expression changes in the guppy exposed to seawater, we characterized its gill transcriptome using RNA sequencing. Experimental fish were exposed to salinity increase from 0 to 30‰ within 4 days, while control fish were cultured in freshwater (0‰ salinity). Seven days after salinity exposure, the gills were sampled and the mortality within 2 weeks was recorded. No significant difference in the cumulative mortality at the second week was found between the two groups. Transcriptomic analysis identified 3477 differentially expressed genes (DEGs), including 1067 upregulated and 2410 downregulated genes. These DEGs were enriched in several biological processes, including ion transport, ion homeostasis, ATP biosynthetic process, metabolic process, and immune system process. Oxidative phosphorylation was the most activated pathway. DEGs involved in the pathway "endoplasmic reticulum (ER)-mediated phagocytosis," "starch and sucrose metabolism," and "steroid biosynthesis" were mainly downregulated; chemokines and interleukins involved in "cytokine-cytokine receptor interaction" were differentially expressed. The present results suggested that oxidative phosphorylation had essential roles in osmoregulation in the gills of seawater acclimated guppy, during which the decline in the expression of genes encoding V-ATPases and calreticulin had a negative effect on the phagocytosis and immune response. Besides, several metabolic processes including "starch and sucrose metabolism" and "steroid biosynthesis" were affected. This study elucidates transcriptomic changes in osmotic regulation, metabolism, and immunity in seawater acclimated guppy.

Effects of Ocean Acidification on Transcriptomes in Asian Seabass Juveniles

Ocean acidification is changing the fate of marine organisms. It is essential to predict the biological responses and evolutionary processes driven by ocean acidification, to maintain the equilibrium of the marine ecosystem and to facilitate aquaculture. However, how marine organisms, particularly the marine fish species, respond to ocean acidification, is still poorly understood. Consequences of ocean acidification on finfish aquaculture are largely not well known. We studied the effects of ocean acidification for 7 days on growth, behaviour and gene expression profiles in the brain, gill and kidney of Asian seabass juveniles. Results showed that growth and behaviour were not affected by short-term ocean acidification. We found tissue-specific differentially expressed genes (DEGs) involving many molecular processes, such as organ development, growth, muscle development, ion homeostasis and neurogenesis and development, as well as behaviours. Most of the DEGs, which were functionally enriched in ion homeostasis, were related to calcium transport, followed by sodium/potassium channels. We found that genes associated with neurogenesis and development were significantly enriched, implying that ocean acidification has also adversely affected the neural regulatory mechanism. Our results indicate that although the short-term ocean acidification does not cause obvious phenotypic and behavioural changes, it causes substantial changes of gene expressions in all three analysed tissues. All these changes of gene expressions may eventually affect physiological fitness. The DEGs identified here should be further investigated to discover DNA markers associated with adaptability to ocean acidification to improve fish's capability to adapt to ocean acidification.

Tissue comparison of transcriptional response to acute acidification stress of barramundi Lates calcarifer in coastal and estuarine areas

In order to explore the common and unique physiological changes in tissues of juvenile barramundi Lates calcarifer in acidified water environment, RNA sequence analysis was used to analyze the molecular responses of liver, head kidney, and gill of juvenile barramundi in pH 7.4 and pH 8.1 seawater environment. The number of differential expression genes identified in liver, head kidney and gill were 860, 388 and 1792, respectively. Through functional enrichment analysis, the differential expression genes common to the three tissues were all related to immunity. Among the unique differential genes in the liver, pathways related to digestion, endocrine, and metabolism were enriched. Among the unique differential expression genes in gill, pathways related to genetic information processing, immunity and metabolism were enriched. The findings of the present study uncover the transcriptional changes in fish correspond to environmental pH change, and provide a better understanding on the biological process at molecular level to environmental pH adapting. This work highlights that assessments for the potential of estuarine fishes to cope with environmental pH change to develop the future conservation strategies.

Transcriptional and Catalytic Responsiveness of the Antarctic Fish Trematomus bernacchii Antioxidant System toward Multiple Stressors

Ocean-warming and acidification jeopardize Antarctic marine species, adapted to cold and constant conditions and naturally exposed to high pro-oxidant pressures and cadmium (Cd) bioavailability. The aim of this study was to investigate if projected temperature increase and pH reduction may affect the accumulation and the effects of Cd in the rockcod Trematomus bernacchii. Organisms were exposed for 14 days to six scenarios, combining environmental or increased temperature (−1 °C, +1 °C) and control or reduced pH (8.05, 7.60), either with or without Cd (40 µg/L). Responses in liver and gills were analyzed at different levels, including mRNA and functional measurements of metallothioneins and of a wide battery of antioxidants, integrated with the evaluation of the total antioxidant capacity and onset of oxidative damages. In the gills, metallothioneins and mRNA of antioxidant genes (nrf2, keap1, cat, gpx1) increased after Cd exposure, but such effects were softened by warming and acidification. Antioxidants showed slighter variations at the enzymatic level, while Cd caused glutathione increase under warming and acidified scenarios. In the liver, due to higher basal antioxidant protection, limited effects were observed. Genotoxic damage increased under the combined stressors scenario. Overall results highlighted the modulation of the oxidative stress response to Cd by multiple stressors, suggesting the vulnerability of T. bernacchii under predicted ocean change scenarios.

Antarctic ecosystems in transition - life between stresses and opportunities

Important findings from the second decade of the 21st century on the impact of environmental change on biological processes in the Antarctic were synthesised by 26 international experts. Ten key messages emerged that have stakeholder-relevance and/or a high impact for the scientific community. They address (i) altered biogeochemical cycles, (ii) ocean acidification, (iii) climate change hotspots, (iv) unexpected dynamism in seabed-dwelling populations, (v) spatial range shifts, (vi) adaptation and thermal resilience, (vii) sea ice related biological fluctuations, (viii) pollution, (ix) endangered terrestrial endemism and (x) the discovery of unknown habitats. Most Antarctic biotas are exposed to multiple stresses and considered vulnerable to environmental change due to narrow tolerance ranges, rapid change, projected circumpolar impacts, low potential for timely genetic adaptation, and migration barriers. Important ecosystem functions, such as primary production and energy transfer between trophic levels, have already changed, and biodiversity patterns have shifted. A confidence assessment of the degree of 'scientific understanding' revealed an intermediate level for most of the more detailed sub-messages, indicating that process-oriented research has been successful in the past decade. Additional efforts are necessary, however, to achieve the level of robustness in scientific knowledge that is required to inform protection measures of the unique Antarctic terrestrial and marine ecosystems, and their contributions to global biodiversity and ecosystem services.

Regulation of globin expression in Antarctic fish under thermal and hypoxic stress

In the freezing waters of the Southern Ocean, Antarctic teleost fish, the Notothenioidei, have developed unique adaptations to cope with cold, including, at the extreme, the loss of hemoglobin in icefish. As a consequence, icefish are thought to be the most vulnerable of the Antarctic fish species to ongoing ocean warming. Some icefish also fail to express myoglobin but all appear to retain neuroglobin, cytoglobin-1, cytoglobin-2, and globin-X. Despite the lack of the inducible heat shock response, Antarctic notothenioid fish are endowed with physiological plasticity to partially compensate for environmental changes, as shown by numerous physiological and genomic/transcriptomic studies over the last decade. However, the regulatory mechanisms that determine temperature/oxygen-induced changes in gene expression remain largely unexplored in these species. Proteins such as globins are susceptible to environmental changes in oxygen levels and temperature, thus playing important roles in mediating Antarctic fish adaptations. In this study, we sequenced the full-length transcripts of myoglobin, neuroglobin, cytoglobin-1, cytoglobin-2, and globin-X from the Antarctic red-blooded notothenioid Trematomus bernacchii and the white-blooded icefish Chionodraco hamatus and evaluated transcripts levels after exposure to high temperature and low oxygen levels. Basal levels of globins are similar in the two species and both stressors affect the expression of Antarctic fish globins in brain, retina and gills. Temperature up-regulates globin expression more effectively in white-blooded than in red-blooded fish while hypoxia strongly up-regulates globins in red-blooded fish, particularly in the gills. These results suggest globins function as regulators of temperature and hypoxia tolerance. This study provides the first insights into globin transcriptional changes in Antarctic fish.

Evolution of chaperome gene expression and regulatory elements in the antarctic notothenioid fishes

Confined within the cold-stable Southern Ocean, Antarctic notothenioid fishes have undergone an evolutionary loss of the inducible heat shock response (HSR), while facing perpetual low-temperature challenges to cellular proteostasis. This study examines how evolution in chronic cold has affected the shared cellular apparatus that mediates proteostasis under normal and heat stressed states. To deduce Antarctic-specific changes, we compared native expression levels across the full suite of chaperome genes and assessed the structural integrity of two crucial HSR regulators - Heat Shock Factor 1 (HSF1) that activates HSR, and heat shock elements (HSEs), the binding sites for HSF1 - between Antarctic fishes and the basal temperate notothenioid Eleginops maclovinus. Native expression levels of Antarctic fish chaperomes showed very modest changes overall, contrary to the common view of constitutive upregulation in the cold. Only a few cytosolic HSP70 genes showed greater transcription, with only the ancestrally-inducible HSPA6 strongly upregulated across all Antarctic species. Additionally, the constant cold has apparently not relaxed the selective pressures on maintaining HSF1 and HSEs in Antarctic fish. Instead, we found HSF1 experienced intensified selective pressure, with conserved sequence changes in Antarctic species suggesting optimization for non-heat-stress functional roles. HSEs of the HSP70 gene family have largely remained conserved in canonical sequence motifs and copy numbers as in E. maclovinus, showing limited impact of relaxed selective pressure. This study shows that evolution in chronic cold has led to both subtle and distinctive changes in the cellular apparatus for proteostasis and HSR, with functional consequences amenable to experimental evaluation.

Reduced stress defence responses contribute to the higher toxicity of a pesticide under warming

There is a pressing need to identify the molecular mechanisms underlying the, often magnifying, interactive effects between contaminants and natural stressors. Here we test our hypothesis that lower general stress defence responses contribute to synergistic interactions between stressors. We focus on the widespread pattern that many contaminants are more toxic at higher temperatures. Specifically, we tested the effects of an environmentally realistic low-effect and high-effect concentration of the pesticide chlorpyrifos under warming at the gene expression level in the northern house mosquito Culex pipiens molestus (Forskal, 1775). By applying the independent action model for combined stressors on RNA-sequencing data, we identified interactive gene expression patterns under combined exposure to chlorpyrifos and warming for general stress defence responses: protection of macromolecules, antioxidant processes, detoxification and energy metabolism/allocation. Most of these general stress defence response genes showed upregulated antagonistic interactions (i.e., were less upregulated than expected under the independent action model). This indicates that when pesticide exposure was combined with warming, the general stress defence responses were no longer buffering increased stress levels, which may contribute to a higher sensitivity to toxicants under warming. These upregulated antagonistic interactions were stronger for the high-effect chlorpyrifos concentration, indicating that exposure to this concentration under warming was most stressful. Our results highlight that quantitative analysis of the frequency and strength of the interaction types of general stress defence response genes, specifically focusing on antagonistic upregulations and synergistic downregulations, may advance our understanding of how natural stressors modify the toxicity of contaminants.

Understanding the Metabolic Capacity of Antarctic Fishes to Acclimate to Future Ocean Conditions

Antarctic fishes have evolved under stable, extreme cold temperatures for millions of years. Adapted to thrive in the cold environment, their specialized phenotypes will likely render them particularly susceptible to future ocean warming and acidification as a result of climate change. Moving from a period of stability to one of environmental change, species persistence will depend on maintaining energetic equilibrium, or sustaining the increased energy demand without compromising important biological functions such as growth and reproduction. Metabolic capacity to acclimate, marked by a return to metabolic equilibrium through physiological compensation of routine metabolic rate (RMR), will likely determine which species will be better poised to cope with shifts in environmental conditions. Focusing on the suborder Notothenioidei, a dominant group of Antarctic fishes, and in particular four well-studied species, Trematomus bernacchii, Pagothenia borchgrevinki, Notothenia rossii, and N. coriiceps, we discuss metabolic acclimation potential to warming and CO2-acidification using an integrative and comparative framework. There are species-specific differences in the physiological compensation of RMR during warming and the duration of acclimation time required to achieve compensation; for some species, RMR fully recovered within 3.5 weeks of exposure, such as P. borchgrevinki, while for other species, such as N. coriiceps, RMR remained significantly elevated past 9 weeks of exposure. In all instances, added exposure to increased PCO2, further compromised the ability of species to return RMR to pre-exposure levels. The period of metabolic imbalance, marked by elevated RMR, was underlined by energetic disturbance and elevated energetic costs, which shifted energy away from fitness-related functions, such as growth. In T. bernacchii and N. coriiceps, long duration of elevated RMR impacted condition factor and/or growth rate. Low growth rate can affect development and ultimately the timing of reproduction, severely compromising the species' survival potential and the biodiversity of the notothenioid lineage. Therefore, the ability to achieve full compensation of RMR, and in a short-time frame, in order to avoid long term consequences of metabolic imbalance, will likely be an important determinant in a species' capacity to persist in a changing environment. Much work is still required to develop our understanding of the bioenergetics of Antarctic fishes in the face of environmental change, and a targeted approach of nesting a mechanistic focus in an ecological and comparative framework will better aid our predictions on the effect of global climate change on species persistence in the polar regions.

Characterizing Gene Copy Number of Heat Shock Protein Gene Families in the Emerald Rockcod, Trematomus bernacchii

The suborder Notothenioidae is comprised of Antarctic fishes, several of which have lost their ability to rapidly upregulate heat shock proteins in response to thermal stress, instead adopting a pattern of expression resembling constitutive genes. Given the cold-denaturing effect that sub-zero waters have on proteins, evolution in the Southern Ocean has likely selected for increased expression of molecular chaperones. These selective pressures may have also enabled retention of gene duplicates, bolstering quantitative output of cytosolic heat shock proteins (HSPs). Given that newly duplicated genes are under more relaxed selection, it is plausible that gene duplication enabled altered regulation of such highly conserved genes. To test for evidence of gene duplication, copy number of various isoforms within major heat shock gene families were characterized via qPCR and compared between the Antarctic notothen, Trematomus bernacchii, which lost the inducible heat shock response, and the non-Antarctic notothen, Notothenia angustata, which maintains an inducible heat shock response. The results indicate duplication of isoforms within the hsp70 and hsp40 super families have occurred in the genome of T. bernacchii. The findings suggest gene duplications may have been critical in maintaining protein folding efficiency in the sub-zero waters and provided an evolutionary mechanism of alternative regulation of these conserved gene families.

Combined stress of ocean acidification and warming influence survival and drives differential gene expression patterns in the Antarctic pteropod, Limacina helicina antarctica

The ecologically important thecosome pteropods in the Limacina spp. complex have recently been the focus of studies examining the impacts global change factors - e.g., ocean acidification (OA) and ocean warming (OW) - on their performance and physiology. This focus is driven by conservation concerns where the health of pteropod populations is threatened by the high susceptibility of their shells to dissolution in low aragonite saturation states associated with OA and how coupling of these stressors may push pteropods past the limits of physiological plasticity. In this manipulation experiment, we describe changes in the transcriptome of the Antarctic pteropod, Limacina helicina antarctica, to these combined stressors. The conditions used in the laboratory treatments met or exceeded those projected for the Southern Ocean by the year 2100. We made two general observations regarding the outcome of the data: (1) Temperature was more influential than pH in terms of changing patterns of gene expression, and (2) these Antarctic pteropods appeared to have a significant degree of transcriptomic plasticity to respond to acute abiotic stress in the laboratory. In general, differential gene expression was observed amongst the treatments; here, for example, transcripts associated with maintaining protein structure and cell proliferation were up-regulated. To disentangle the effects of OA and OW, we used a weighted gene co-expression network analysis to explore patterns of change in the transcriptome. This approach identified gene networks associated with OW that were enriched for transcripts proposed to be involved in increasing membrane fluidity at warmer temperatures. Together these data provide evidence that L.h.antarctica has a limited capacity to acclimate to the combined conditions of OA and OW used in this study. This reduced scope of acclimation argues for continued study of how adaptation to polar aquatic environments may limit the plasticity of present-day populations in responding to future environmental change.

Ocean acidification promotes broad transcriptomic responses in marine metazoans: a literature survey

For nearly a decade, the metazoan-focused research community has explored the impacts of ocean acidification (OA) on marine animals, noting that changes in ocean chemistry can impact calcification, metabolism, acid-base regulation, stress response and behavior in organisms that hold high ecological and economic value. Because OA interacts with several key physiological processes in marine organisms, transcriptomics has become a widely-used method to characterize whole organism responses on a molecular level as well as inform mechanisms that explain changes in phenotypes observed in response to OA. In the past decade, there has been a notable rise in studies that examine transcriptomic responses to OA in marine metazoans, and here we attempt to summarize key findings across these studies. We find that organisms vary dramatically in their transcriptomic responses to pH although common patterns are often observed, including shifts in acid-base ion regulation, metabolic processes, calcification and stress response mechanisms. We also see a rise in transcriptomic studies examining organismal response to OA in a multi-stressor context, often reporting synergistic effects of OA and temperature. In addition, there is an increase in studies that use transcriptomics to examine the evolutionary potential of organisms to adapt to OA conditions in the future through population and transgenerational experiments. Overall, the literature reveals complex organismal responses to OA, in which some organisms will face more dramatic consequences than others. This will have wide-reaching impacts on ocean communities and ecosystems as a whole.

Accelerated evolution at chaperone promoters among Antarctic notothenioid fishes

BACKGROUND

Antarctic fishes of the Notothenioidei suborder constitutively upregulate multiple inducible chaperones, a highly derived adaptation that preserves proteostasis in extreme cold, and represent a system for studying the evolution of gene frontloading. We screened for Hsf1-binding sites, as Hsf1 is a master transcription factor of the heat shock response, and highly-conserved non-coding elements within proximal promoters of chaperone genes across 10 Antarctic notothens, 2 subpolar notothens, and 17 perciform fishes. We employed phylogenetic models of molecular evolution to determine whether (i) changes in motifs associated with Hsf1-binding and/or (ii) relaxed purifying selection or exaptation at ancestral cis-regulatory elements coincided with the evolution of chaperone frontloading in Antarctic notothens.

RESULTS

Antarctic notothens exhibited significantly fewer Hsf1-binding sites per bp at chaperone promoters than subpolar notothens and Serranoidei, the most closely-related suborder to Notothenioidei included in this study. 90% of chaperone promoters exhibited accelerated substitution rates among Antarctic notothens relative to other perciformes. The proportion of bases undergoing accelerated evolution (i) was significantly greater in Antarctic notothens than in subpolar notothens and Perciformes in 70% of chaperone genes and (ii) increased among bases that were more conserved among perciformes. Lastly, we detected evidence of relaxed purifying selection and exaptation acting on ancestrally conserved cis-regulatory elements in the Antarctic notothen lineage and its major branches.

CONCLUSION

A large degree of turnover has occurred in Notothenioidei at chaperone promoter regions that are conserved among perciform fishes following adaptation to the cooling of the Southern Ocean. Additionally, derived reductions in Hsf1-binding site frequency suggest cis-regulatory modifications to the classical heat shock response. Of note, turnover events within chaperone promoters were less frequent in the ancestral node of Antarctic notothens relative to younger Antarctic lineages. This suggests that cis-regulatory divergence at chaperone promoters may be greater between Antarctic notothen lineages than between subpolar and Antarctic clades. These findings demonstrate that strong selective forces have acted upon cis-regulatory elements of chaperone genes among Antarctic notothens.

Modeling the influence of time and temperature on the levels of fatty acids in the liver of Antarctic fish Trematomus bernacchii

Antarctic fish (Trematomus bernacchii) are an ideal group for studying the effect of ocean warming on vital physiological and biochemical mechanisms of adaptation, including changes in the fatty acid composition to higher heat tolerance in the sub-zero waters of the Southern Ocean. Despite the awareness of the impact of ocean warming on marine life, bioclimatic models describing the effect of temperature and time on fatty acid levels in marine species have not been considered yet. The objective of the present study was to investigate changes in the concentrations of fatty acids in liver from T. bernacchii in response to an increase in temperature in the Antarctic region. Changes in the concentrations of fatty acids in liver from T. bernacchii were observed after varying simultaneously and systematically the temperature and time. The fatty acid profiles were determined by gas chromatography prior to acclimation (− 1.8 °C) and after acclimation (0.0, 1.0, and 2.0 °C) at different times (1, 5, and 10 days). The observed changes were graphically visualized by expressing the fatty acid concentration in absolute units (mg g−1) as a function of the temperature and time using polynomial models. Major changes in fatty acid composition were observed at day 1 of exposition at all temperatures. At day 5, the fish seem to tolerate the new temperature condition. The concentrations of saturated fatty acids were almost constant throughout the various conditions. The concentrations of monounsaturated fatty acids (in particular 18:1n − 9) decrease at day 1 for all temperatures. In contrast, there was an increase in the concentrations of polyunsaturated fatty acids (in particular 20:5n − 3 and 22:6n − 3) with increasing temperatures after 1, 5, and 10 days of exposure. The proposed models were in agreement with reported studies on polar and temperate fish, indicating possibly similar adaptation mechanisms for teleost to cope with global warming.

Characterization of thermally sensitive miRNAs reveals a central role of the FoxO signaling pathway in regulating the cellular stress response of an extreme stenotherm, Trematomus bernacchii

Despite the lack of an inducible heat shock response (HSR), the Antarctic notothenioid fish, Trematomus bernacchii, has retained a level of physiological plasticity that can at least partially compensate for the effects of acute heat stress. Over the last decade, both physiological and transcriptomic studies have signaled these fish can mitigate the effects of acute heat stress by employing other aspects of the cellular stress response (CSR) that help confer thermotolerance as well as drive homeostatic mechanisms during long-term thermal acclimations. However, the regulatory mechanisms that determine temperature-induced changes in gene expression remain largely unexplored in this species. Therefore, this study utilized next generation sequencing coupled with an in silico approach to explore the regulatory role of microRNAs in governing the transcriptomic level response observed in this Antarctic notothenioid with respect to the CSR. Using RNAseq, we characterized the expression of 125 distinct miRNA orthologues in T. bernacchii gill tissue. Additionally, we identified 12 miRNAs that appear to be thermally responsive based on differential expression (DE) analyses performed between fish acclimated to control (^-1.5 °C) and an acute heat stress (⁺4 °C). We further characterized the functional role of these DE miRNAs using bioinformatics pipelines to identify putative gene targets of the DE miRNAs and subsequent gene set enrichment analyses, which together suggest these miRNAs are involved in regulating diverse aspects of the CSR in T. bernacchii.

Divergence, evolution and adaptation in ray-finned fish genomes

With the rapid development of next-generation sequencing technologies and bioinformatics, over 50 ray-finned fish genomes by far have been sequenced with high quality. The genomic work provides abundant genetic resources for deep understanding of divergence, evolution and adaptation in the fish genomes. They are also instructive for identification of candidate genes for functional verification, molecular breeding, and development of novel marine drugs. As an example of other omics data, the Fish-T1K project generated a big database of fish transcriptomes to integrate with these published fish genomes for potential applications. In this review, we highlight the above-mentioned recent investigations and core topics on the ray-finned fish genome research, with a main goal to obtain a deeper understanding of fish biology for theoretical and practical applications.

Transcriptome Analysis of Yamame (Oncorhynchus masou) in Normal Conditions after Heat Stress

Understanding the mechanism of high-temperature tolerance in cold-freshwater fish is crucial for predicting how certain species will cope with global warming. In this study, we investigated temperature tolerance in masu salmon (Oncorhynchus masou, known in Japan as ‘yamame’), an important aquaculture species. By selective breeding, we developed a group of yamame (F2) with high-temperature tolerance. This group was subjected to a high-temperature tolerance test and divided into two groups: High-temperature tolerant (HT) and non-high-temperature tolerant (NT). RNA was extracted from the gill and adipose fin tissues of each group, and the mRNA expression profiles were analyzed using RNA sequencing. A total of 2893 differentially expressed genes (DEGs) from the gill and 836 from the adipose fin were identified by comparing the HT and NT groups. Functional analyses were then performed to identify associated gene ontology (GO) terms and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The HT group showed a high expression of heat shock protein 70 (HSP70) gene and enriched gene expression in the extracellular matrix (ECM), cell junction, and adhesion pathways in gill tissues compared to the NT group. The HT group also exhibited highly expressed genes in glycolysis and showed lower expression of the genes in the p53 signaling pathway in adipose fin tissues. Taken together, the difference of expression of some genes in the normal condition may be responsible for the difference in heat tolerance between the HT and NT yamame in the heat stress condition.

The transcriptome of the marine calanoid copepod Temora longicornis under heat stress and recovery

Understanding the impacts of global change in zooplankton communities is crucial, as alterations in the zooplankton communities can affect entire marine ecosystems. Despite the economic and ecological importance of the calanoid copepod Temora longicornis in the Belgian part of the North Sea, molecular data is still very limited for this species. Using HiSeq Illumina sequencing, we sequenced the whole transcriptome of T. longicornis, after being exposed to realistic temperatures of 14 and 17 °C. After both an acute (1 day) and a more sustained (5 days) thermal exposure to 17 °C, we investigated gene expression differences with animals exposed to 14 °C, which may be critical for the thermal acclimation and resilience of this copepod species. We also studied the possibility of a short term stress recovery of a heat shock. A total of 179,569 transcripts were yielded, of which 44,985 putative ORF transcripts were identified. These transcripts were subsequently annotated into roughly 22,000 genes based on known sequences using Gene Ontology (GO) and KEGG databases. Temora only showed a mild response to both the temperature and the duration of the exposure. We found that the expression of 27 transcripts varied significantly with an increase in temperature of 3 °C, of which eight transcripts were differentially expressed after acute exposure only. Gene set enrichment analysis revealed that, overall, T. longicornis was more impacted by a sustained thermal exposure, rather than an immediate (acute) exposure, with two times as many enriched GO terms in the sustained treatment. We also identified several general stress responses independent of exposure time, such as modified protein synthesis, energy mobilisation, cuticle and chaperone proteins. Finally, we highlighted candidate genes of a possible recovery from heat exposure, identifying similar terms as those enriched in the heat treatments, i.e. related to for example energy metabolism, cuticle genes and extracellular matrix. The data presented in this study provides the first transcriptome available for T. longicornis which can be used for future genomic studies.

Evolution in chronic cold: varied loss of cellular response to heat in Antarctic notothenioid fish

Background

Confined within the freezing Southern Ocean, the Antarctic notothenioids have evolved to become both cold adapted and cold specialized. A marked signature of cold specialization is an apparent loss of the cellular heat shock response (HSR). As the HSR has been examined in very few notothenioid species to-date, it remains unknown whether HSR loss pervades the Antarctic radiation, or whether the broader cellular responses to heat stress has sustained similar loss. Understanding the evolutionary status of these responses in this stenothermal taxon is crucial for evaluating its adaptive potential to ocean warming under climate change.

Results

In this study, we used an acute heat stress protocol followed by RNA-Seq analyses to study the evolution of cellular-wide transcriptional responses to heat stress across three select notothenioid lineages - the basal temperate and nearest non-Antarctic sister species Eleginops maclovinus serving as ancestral proxy, the cryopelagic Pagothenia borchgrevinki and the icefish Chionodraco rastrospinosus representing cold-adapted red-blooded and hemoglobinless Antarctic notothenioids respectively. E. maclovinus displayed robust cellular stress responses including the ER Unfolded Protein Response and the cytosolic HSR, cementing the HSR as a plesiomorphy that preceded Antarctic notothenioid radiation. While the transcriptional response to heat stress was minimal in P. borchgrevinki, C. rastrospinosus exhibited robust responses in the broader cellular networks especially in inflammatory responses despite lacking the classic HSR and UPR.

Conclusion

The disparate patterns observed in these two archetypal Antarctic species indicate the evolutionary status in cellular ability to mitigate acute heat stress varies even among Antarctic lineages, which may affect their adaptive potential in coping with a warming world.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1254-6) contains supplementary material, which is available to authorized users.

Heat shock influences the fatty acid composition of the muscle of the Antarctic fish Trematomus bernacchii

In the Ross Sea region (average temperature of -1.87 °C), shelf water warming up to +0.8-+1.4 °C is predicted by 2200, so there is an urgent need to understand how organisms can respond to rising temperatures. In this study, we analyzed the effect of a heat shock on the fatty acid (FAs) composition of muscle of the Antarctic teleost Trematomus bernacchii, caught in Terra Nova Bay (Ross Sea), and held in fish tanks at 0, +1 or +2 °C, for 1, 5 and 10 days. In general, heat shock produced, beyond a reduction in total lipid content correlated to the temperature, an increase in the percentage of saturated FAs, and a decrease in mono-unsaturated FAs; however, the level of poly-unsaturated FAs did not seem to directly correlate with temperature. Principal component analysis indicated that both temperature and exposure time affect the composition of FAs in the muscle probably through an alteration of the metabolic pathways of FAs. In this study, we demonstrated that T. bernacchii was capable to rapidly acclimatize to a heat shock. This study contributes to increasing knowledge on the effect of temperature on the lipid composition of T. bernacchii and is complementary to previous studies on the gene expression and biochemistry of this species face multiple stressors.

Antarctic emerald rockcod have the capacity to compensate for warming when uncoupled from CO2 -acidification

Increases in atmospheric CO₂ levels and associated ocean changes are expected to have dramatic impacts on marine ecosystems. Although the Southern Ocean is experiencing some of the fastest rates of change, few studies have explored how Antarctic fishes may be affected by co-occurring ocean changes, and even fewer have examined early life stages. To date, no studies have characterized potential trade-offs in physiology and behavior in response to projected multiple climate change stressors (ocean acidification and warming) on Antarctic fishes. We exposed juvenile emerald rockcod Trematomus bernacchii to three PCO₂ treatments (~450, ~850, and ~1,200 μatm PCO₂ ) at two temperatures (-1 or 2°C). After 2, 7, 14, and 28 days, metrics of physiological performance including cardiorespiratory function (heart rate [f_H ] and ventilation rate [f_V ]), metabolic rate (M˙O2), and cellular enzyme activity were measured. Behavioral responses, including scototaxis, activity, exploration, and escape response were assessed after 7 and 14 days. Elevated PCO₂ independently had little impact on either physiology or behavior in juvenile rockcod, whereas warming resulted in significant changes across acclimation time. After 14 days, f_H , f_V and M˙O2 significantly increased with warming, but not with elevated PCO₂ . Increased physiological costs were accompanied by behavioral alterations including increased dark zone preference up to 14%, reduced activity by 12%, as well as reduced escape time suggesting potential trade-offs in energetics. After 28 days, juvenile rockcod demonstrated a degree of temperature compensation as f_V , M˙O2, and cellular metabolism significantly decreased following the peak at 14 days; however, temperature compensation was only evident in the absence of elevated PCO₂ . Sustained increases in f_V and M˙O2 after 28 days exposure to elevated PCO₂ indicate additive (f_V ) and synergistic (M˙O2) interactions occurred in combination with warming. Stressor-induced energetic trade-offs in physiology and behavior may be an important mechanism leading to vulnerability of Antarctic fishes to future ocean change.

Transcriptome profiling based on protein-protein interaction networks provides a core set of genes for understanding blood immune response mechanisms against Edwardsiella tarda infection in Japanese flounder (Paralichthys olivaceus)

Marine organisms are commonly under threat from various pathogens. Edwardsiella tarda is one of the fish pathogens that can infect both cultured and wild fish species. E. tarda can also infect other vertebrates, including amphibians, reptiles, and mammals. Bacteremia caused by E. tarda can be a fatal disease in humans. Blood acts as a pipeline for the fish immune system. Generating blood transcriptomic resources from fish challenged by E. tarda is crucial for understanding molecular mechanisms underlying blood immune response process. In this study, we performed transcriptome-wide gene expression profiling of Japanese flounder (Paralichthys olivaceus) challenged by 8 and 48 h E. tarda stress. An average of 37 million clean reads per library was obtained, and approximately 85.6% of these reads were successfully mapped to the reference genome. In addition, 808 and 1265 differential expression genes (DEGs) were found at 8 and 48 h post-injection, respectively. Gene Ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to search immune-related DEGs. A protein-protein interaction network was constructed to obtain the interaction relationship of immune genes during pathogens stress. Based on KEGG and protein association networks analysis, 30 hub genes were discovered and validated by quantitative RT-PCR. This study represents the first transcriptome analysis based on protein-protein interaction networks in fish and provides us with valuable gene resources for the research of fish blood immunity, which can significantly assist us to further understand the molecular mechanisms of humans and other vertebrates against E. tarda.

Biochemical, Histopathological and Molecular Responses in Gills of Leuciscus cephalus Exposed to Metals

Gills are major targets for acute metal toxicity in fish, due to their permanent contact with aquatic pollutants. To assess the effects of metals on gills of the Leuciscus cephalus (chub), fish individuals were collected from two sites in the Tur River, Romania, in upstream (site 1) and downstream (site 2) of a metal pollution source. Quantitative and hyperspectral analyses showed that Zn, Sr, and Fe concentrations were significantly higher in gills from site 2 compared with site 1. Malondialdehyde and advanced oxidation protein products levels increased 17 and 28%, respectively, whereas reduced glutathione level diminished significantly in the gills of fish collected from site 2 compared to site 1. The activities of superoxide dismutase, catalase, and glutathione-S-transferase increased significantly at 41, 21, and 28%, respectively. Proliferating cell nuclear antigen (PCNA) protein levels, as well as the amount of DNA damage, were significantly increased for site 2 compared with site 1. The induced oxidative stress generated hyperplasia, hypertrophy, and inflammation in the epithelial cells and apoptosis. Hence, this could suggest that gill cells have tried to counteract the oxidative stress-induced DNA fragmentation by PCNA up-regulation, but the PCNA expression decreased on longer time due to the low level of GSH, resulting in apoptosis.

Transcriptomic response of the Antarctic pteropod Limacina helicina antarctica to ocean acidification

Background

Ocean acidification (OA), a change in ocean chemistry due to the absorption of atmospheric CO₂ into surface oceans, challenges biogenic calcification in many marine organisms. Ocean acidification is expected to rapidly progress in polar seas, with regions of the Southern Ocean expected to experience severe OA within decades. Biologically, the consequences of OA challenge calcification processes and impose an energetic cost.

Results

In order to better characterize the response of a polar calcifier to conditions of OA, we assessed differential gene expression in the Antarctic pteropod, Limacina helicina antarctica. Experimental levels of pCO₂ were chosen to create both contemporary pH conditions, and to mimic future pH expected in OA scenarios. Significant changes in the transcriptome were observed when juvenile L. h. antarctica were acclimated for 21 days to low-pH (7.71), mid-pH (7.9) or high-pH (8.13) conditions. Differential gene expression analysis of individuals maintained in the low-pH treatment identified down-regulation of genes involved in cytoskeletal structure, lipid transport, and metabolism. High pH exposure led to increased expression and enrichment for genes involved in shell formation, calcium ion binding, and DNA binding. Significant differential gene expression was observed in four major cellular and physiological processes: shell formation, the cellular stress response, metabolism, and neural function. Across these functional groups, exposure to conditions that mimic ocean acidification led to rapid suppression of gene expression.

Conclusions

Results of this study demonstrated that the transcriptome of the juvenile pteropod, L. h. antarctica, was dynamic and changed in response to different levels of pCO₂. In a global change context, exposure of L. h. antarctica to the low pH, high pCO₂ OA conditions resulted in a suppression of transcripts for genes involved in key physiological processes: calcification, metabolism, and the cellular stress response. The transcriptomic response at both acute and longer-term acclimation time frames indicated that contemporary L. h. antarctica may not have the physiological plasticity necessary for adaptation to OA conditions expected in future decades. Lastly, the differential gene expression results further support the role of shelled pteropods such as L. h. antarctica as sentinel organisms for the impacts of ocean acidification.

Electronic supplementary material

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

Transcriptomic responses to environmental change in fishes: Insights from RNA sequencing

The need to better understand how plasticity and evolution affect organismal responses to environmental variability is paramount in the face of global climate change. The potential for using RNA sequencing (RNA-seq) to study complex responses by non-model organisms to the environment is evident in a rapidly growing body of literature. This is particularly true of fishes for which research has been motivated by their ecological importance, socioeconomic value, and increased use as model species for medical and genetic research. Here, we review studies that have used RNA-seq to study transcriptomic responses to continuous abiotic variables to which fishes have likely evolved a response and that are predicted to be affected by climate change (e.g., salinity, temperature, dissolved oxygen concentration, and pH). Field and laboratory experiments demonstrate the potential for individuals to respond plastically to short- and long-term environmental stress and reveal molecular mechanisms underlying developmental and transgenerational plasticity, as well as adaptation to different environmental regimes. We discuss experimental, analytical, and conceptual issues that have arisen from this work and suggest avenues for future study.

Regulation of gene expression is associated with tolerance of the Arctic copepod Calanus glacialis to CO2-acidified sea water

Ocean acidification is the increase in seawater pCO ₂ due to the uptake of atmospheric anthropogenic CO ₂, with the largest changes predicted to occur in the Arctic seas. For some marine organisms, this change in pCO ₂, and associated decrease in pH, represents a climate change-related stressor. In this study, we investigated the gene expression patterns of nauplii of the Arctic copepod Calanus glacialis cultured at low pH levels. We have previously shown that organismal-level performance (development, growth, respiration) of C. glacialis nauplii is unaffected by low pH. Here, we investigated the molecular-level response to lowered pH in order to elucidate the physiological processes involved in this tolerance. Nauplii from wild-caught C. glacialis were cultured at four pH levels (8.05, 7.9, 7.7, 7.5). At stage N6, mRNA was extracted and sequenced using RNA-seq. The physiological functionality of the proteins identified was categorized using Gene Ontology and KEGG pathways. We found that the expression of 151 contigs varied significantly with pH on a continuous scale (93% downregulated with decreasing pH). Gene set enrichment analysis revealed that, of the processes downregulated, many were components of the universal cellular stress response, including DNA repair, redox regulation, protein folding, and proteolysis. Sodium:proton antiporters were among the processes significantly upregulated, indicating that these ion pumps were involved in maintaining cellular pH homeostasis. C. glacialis significantly alters its gene expression at low pH, although they maintain normal larval development. Understanding what confers tolerance to some species will support our ability to predict the effects of future ocean acidification on marine organisms.

The effects of elevated temperature and ocean acidification on the metabolic pathways of notothenioid fish

The adaptations used by notothenioid fish to combat extreme cold may have left these fish poorly poised to deal with a changing environment. As such, the expected environmental perturbations brought on by global climate change have the potential to significantly affect the energetic demands and subsequent cellular processes necessary for survival. Despite recent lines of evidence demonstrating that notothenioid fish retain the ability to acclimate to elevated temperatures, the underlying mechanisms responsible for temperature acclimation in these fish remain largely unknown. Furthermore, little information exists on the capacity of Antarctic fish to respond to changes in multiple environmental variables. We have examined the effects of increased temperature and pCO₂ on the rate of oxygen consumption in three notothenioid species, Trematomus bernacchii, Pagothenia borchgrevinki, and Trematomus newnesi. We combined these measurements with analysis of changes in aerobic and anaerobic capacity, lipid reserves, fish condition, and growth rates to gain insight into the metabolic cost associated with acclimation to this dual stress. Our findings indicated that temperature is the major driver of the metabolic responses observed in these fish and that increased pCO₂ plays a small, contributing role to the energetic costs of the acclimation response. All three species displayed varying levels of energetic compensation in response to the combination of elevated temperature and pCO₂. While P. borchgrevinki showed nearly complete compensation of whole animal oxygen consumption rates and aerobic capacity, T. newnesi and T. bernacchii displayed only partial compensation in these metrics, suggesting that at least some notothenioids may require physiological trade-offs to fully offset the energetic costs of long-term acclimation to climate change related stressors.

Species-Specific Responses of Juvenile Rockfish to Elevated pCO2: From Behavior to Genomics

In the California Current ecosystem, global climate change is predicted to trigger large-scale changes in ocean chemistry within this century. Ocean acidification—which occurs when increased levels of atmospheric CO₂ dissolve into the ocean—is one of the biggest potential threats to marine life. In a coastal upwelling system, we compared the effects of chronic exposure to low pH (elevated pCO₂) at four treatment levels (i.e., pCO₂ = ambient [500], moderate [750], high [1900], and extreme [2800 μatm]) on behavior, physiology, and patterns of gene expression in white muscle tissue of juvenile rockfish (genus Sebastes), integrating responses from the transcriptome to the whole organism level. Experiments were conducted simultaneously on two closely related species that both inhabit kelp forests, yet differ in early life history traits, to compare high-CO₂ tolerance among species. Our findings indicate that these congeners express different sensitivities to elevated CO₂ levels. Copper rockfish (S. caurinus) exhibited changes in behavioral lateralization, reduced critical swimming speed, depressed aerobic scope, changes in metabolic enzyme activity, and increases in the expression of transcription factors and regulatory genes at high pCO₂ exposure. Blue rockfish (S. mystinus), in contrast, showed no significant changes in behavior, swimming physiology, or aerobic capacity, but did exhibit significant changes in the expression of muscle structural genes as a function of pCO₂, indicating acclimatization potential. The capacity of long-lived, late to mature, commercially important fish to acclimatize and adapt to changing ocean chemistry over the next 50–100 years is likely dependent on species-specific physiological tolerances.

Aquatic insects in a multistress environment: cross-tolerance to salinity and desiccation

Exposing organims to a particular stressor may enhance tolerance to a subsequent stress, when protective mechanisms against both stressors are shared. Such cross-tolerance is a common adaptive response in dynamic multivariate environments and often indicates potential co-evolution of stress traits. Many aquatic insects in inland saline waters from Mediterranean-climate regions are sequentially challenged with salinity and desiccation stress. Thus, cross-tolerance to these physiologically similar stressors could have been positively selected in insects of these regions. We used adults of the saline water beetles Enochrus jesusarribasi (Hydrophilidae) and Nebrioporus baeticus (Dytiscidae) to test cross-tolerance responses to desiccation and salinity. In independent laboratory experiments, we evaluated the effects of i) salinity stress on the subsequent resistance to desiccation and ii) desiccation stress (rapid and slow dehydration) on the subsequent tolerance to salinity. Survival, water loss and haemolymph osmolality were measured. Exposure to stressful salinity improved water control under subsequent desiccation stress in both species, with a clear cross-tolerance (enhanced performance) in N. baeticus. In contrast, general negative effects on performance were found under the inverse stress sequence. The rapid and slow dehydration produced different water loss and haemolymph osmolality dynamics that were reflected in different survival patterns. Our finding of cross-tolerance to salinity and desiccation in ecologically similar species from distant lineages, together with parallel responses between salinity and thermal stress previously found in several aquatic taxa, highlights the central role of adaption to salinity and co-occurring stressors in arid inland waters, having important implications for the species' persistence under climate change.

Seasonal immunoregulation in a naturally-occurring vertebrate

Background

Fishes show seasonal patterns of immunity, but such phenomena are imperfectly understood in vertebrates generally, even in humans and mice. As these seasonal patterns may link to infectious disease risk and individual condition, the nature of their control has real practical implications. Here we characterize seasonal dynamics in the expression of conserved vertebrate immunity genes in a naturally-occurring piscine model, the three-spined stickleback.

Results

We made genome-wide measurements (RNAseq) of whole-fish mRNA pools (n = 36) at the end of summer and winter in contrasting habitats (riverine and lacustrine) and focussed on common trends to filter habitat-specific from overarching temporal responses. We corroborated this analysis with targeted year-round whole-fish gene expression (Q-PCR) studies in a different year (n = 478). We also considered seasonal tissue-specific expression (6 tissues) (n = 15) at a third contrasting (euryhaline) locality by Q-PCR, further validating the generality of the patterns seen in whole fish analyses. Extremes of season were the dominant predictor of immune expression (compared to sex, ontogeny or habitat). Signatures of adaptive immunity were elevated in late summer. In contrast, late winter was accompanied by signatures of innate immunity (including IL-1 signalling and non-classical complement activity) and modulated toll-like receptor signalling. Negative regulators of T-cell activity were prominent amongst winter-biased genes, suggesting that adaptive immunity is actively down-regulated during winter rather than passively tracking ambient temperature. Network analyses identified a small set of immune genes that might lie close to a regulatory axis. These genes acted as hubs linking summer-biased adaptive pathways, winter-biased innate pathways and other organismal processes, including growth, metabolic dynamics and responses to stress and temperature. Seasonal change was most pronounced in the gill, which contains a considerable concentration of T-cell activity in the stickleback.

Conclusions

Our results suggest major and predictable seasonal re-adjustments of immunity. Further consideration should be given to the effects of such responses in seasonally-occurring disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2701-7) contains supplementary material, which is available to authorized users.

RNA-seq reveals a diminished acclimation response to the combined effects of ocean acidification and elevated seawater temperature in Pagothenia borchgrevinki

PURPOSE

The IPCC has reasserted the strong influence of anthropogenic CO2 contributions on global climate change and highlighted the polar-regions as highly vulnerable. With these predictions the cold adapted fauna endemic to the Southern Ocean, which is dominated by fishes of the sub-order Notothenioidei, will face considerable challenges in the near future. Recent physiological studies have demonstrated that the synergistic stressors of elevated temperature and ocean acidification have a considerable, although variable, impact on notothenioid fishes. The present study explored the transcriptomic response of Pagothenia borchgrevinki to increased temperatures and pCO2 after 7, 28 and 56days of acclimation. We compared this response to short term studies assessing heat stress alone and foretell the potential impacts of these stressors on P. borchgrevinki's ability to survive a changing Southern Ocean.

RESULTS

P. borchgrevinki did demonstrate a coordinated stress response to the dual-stressor condition, and even indicated that some level of inducible heat shock response may be conserved in this notothenioid species. However, the stress response of P. borchgrevinki was considerably less robust than that observed previously in the closely related notothenioid, Trematomus bernacchii, and varied considerably when compared across different acclimation time-points. Furthermore, the molecular response of these fish under multiple stressors displayed distinct differences compared to their response to short term heat stress alone.

CONCLUSIONS

When exposed to increased sea surface temperatures, combined with ocean acidification, P. borchgrevinki demonstrated a coordinated stress response that has already peaked by 7days of acclimation and quickly diminished over time. However, this response is less dramatic than other closely related notothenioids under identical conditions, supporting previous research suggesting that this notothenioid species is less sensitive to environmental variation.