Linking transcriptional responses to organismal tolerance reveals mechanisms of thermal sensitivity in a mesothermal endangered fish

Non-concordant epigenetic and transcriptional responses to acute thermal stress in western mosquitofish (Gambusia affinis)

Climate change is intensifying the frequency and severity of extreme temperatures. Understanding the molecular mechanisms underlying the ability to cope with acute thermal stress is key for predicting species' responses to extreme temperature events. While many studies have focused on the individual roles of gene expression, post-transcriptional processes and epigenetic modifications in response to acute thermal stress, the relative contribution of these molecular mechanisms remains unclear. The wide range of thermal limits of western mosquitofish (Gambusia affinis) provides an opportunity to explore this interplay. Here, we quantified changes in gene expression, alternative splicing, DNA methylation and microRNA (miRNA) expression in muscle tissue dissected from mosquitofish immediately after reaching high (CTmax) or low thermal limit (CTmin). Although the numbers of genes showing expression and splicing changes in response to acute temperature stress were small, we found a possibly larger and non-redundant role of splicing compared to gene expression, with more genes being differentially spliced (DSGs) than differentially expressed (DEGs), and little overlap between DSGs and DEGs. We also identified a small proportion of CpGs showing significant methylation change (i.e. differentially methylated cytosines, DMCs) in fish at thermal limits; however, there was no overlap between DEGs and genes annotated with DMCs in both CTmax and CTmin experiments. The weak interplay between epigenetic modifications and gene expression was further supported by our discoveries of no differentially expressed miRNAs. These findings provide novel insights into the relative role of different molecular mechanisms underlying immediate responses to extreme temperatures and demonstrate non-concordant responses of epigenetic and transcriptional mechanisms to acute temperature stress.

Heatwave resilience of juvenile white sturgeon is associated with epigenetic and transcriptional alterations

Heatwaves are increasing in frequency and severity, posing a significant threat to organisms globally. In aquatic environments heatwaves are often associated with low environmental oxygen, which is a deadly combination for fish. However, surprisingly little is known about the capacity of fishes to withstand these interacting stressors. This issue is particularly critical for species of extreme conservation concern such as sturgeon. We assessed the tolerance of juvenile white sturgeon from an endangered population to heatwave exposure and investigated how this exposure affects tolerance to additional acute stressors. We measured whole-animal thermal and hypoxic performance and underlying epigenetic and transcriptional mechanisms. Sturgeon exposed to a simulated heatwave had increased thermal tolerance and exhibited complete compensation for the effects of acute hypoxia. These changes were associated with an increase in mRNA levels involved in thermal and hypoxic stress (hsp90a, hsp90b, hsp70 and hif1a) following these stressors. Global DNA methylation was sensitive to heatwave exposure and rapidly responded to acute thermal and hypoxia stress over the course of an hour. These data demonstrate that juvenile white sturgeon exhibit substantial resilience to heatwaves, associated with improved cross-tolerance to additional acute stressors and involving rapid responses in both epigenetic and transcriptional mechanisms.

Time-Course and Tissue-Specific Molecular Responses to Acute Thermal Stress in Japanese Mantis Shrimp Oratosquilla oratoria

Current understanding of adaptability to high temperatures is increasingly important as extreme weather events that can trigger immediate physiological stress in organisms have occurred more frequently. Here, we tracked transcriptomic responses of gills, hepatopancreas, and muscle to acute thermal exposure at 30 °C for 0.5, 6, and 12 h in an economically important crustacean, Oratosquilla oratoria, to gain a preliminary understanding of the tissue-specific and dynamic physiological regulation process under acute heat stress. The unique physiological responses of muscle, hepatopancreas, and gills to acute thermal stress were associated with protein degradation, lipid transport, and energy metabolism in O. oratoria, respectively. Functional enrichment analysis of differentially expressed transcripts and heat-responsive gene clusters revealed a biphasic protective responsiveness of O. oratoria developed from the early responses of signal transduction, immunity, and cytoskeleton reorganization to the responses dominated by protein turnover and energy metabolism at the mid-late stages under acute heat stress. Noteworthy, trend analysis revealed a consistently upregulated expression pattern of high molecular weight heat shock protein (HSP) family members (HSP60, HSP70, and HSP90) during the entire thermal exposure process, highlighting their importance for maintaining heat resistance in O. oratoria. Documenting the whole process of transcriptional responses at fine temporal resolution will contribute to a far-reaching comprehension of plastic responses to acute heat stress in crustaceans, which is critical in the context of a changing climate.

Effects of turbidity, temperature and predation cue on the stress response of juvenile delta smelt

The San Francisco Estuary (SFE) is one of the most degraded ecosystems in the United States, and organisms that inhabit it are exposed to a suite of environmental stressors. The delta smelt (Hypomesus transpacificus), a small semi-anadromous fish endemic to the SFE and considered an indicator species, is close to extinction in the wild. The goal of this study was to investigate how environmental alterations to the SFE, such as reductions in turbidities, higher temperatures and increased prevalence of invasive predators affect the physiology and stress response of juvenile delta smelt. Juvenile delta smelt were exposed to two temperatures (17 and 21°C) and two turbidities (1-2 and 10-11 NTU) for 2 weeks. After the first week of exposure, delta smelt were exposed to a largemouth bass (Micropterus salmoides) predator cue at the same time every day for 7 days. Fish were measured and sampled on the first (acute) and final (chronic) day of exposures to predator cues and later analyzed for whole-body cortisol, glucose, lactate, and protein. Length and mass measurements were used to calculate condition factor of fish in each treatment. Turbidity had the greatest effect on juvenile delta smelt and resulted in reduced cortisol, increased glucose and lactate, and greater condition factor. Elevated temperatures reduced available energy in delta smelt, indicated by lower glucose and total protein, whereas predator cue exposure had negligible effects on their stress response. This is the first study to show reduced cortisol in juvenile delta smelt held in turbid conditions and adds to the growing data that suggest this species performs best in moderate temperatures and turbidities. Multistressor experiments are necessary to understand the capacity of delta smelt to respond to the multivariate and dynamic changes in their natural environment, and results from this study should be considered for management-based conservation efforts.

Polygenic discrimination of migratory phenotypes in an estuarine forage fish

Migration is a complex phenotypic trait with some species containing migratory and nonmigratory individuals. Such life history variation may be attributed in part to plasticity, epigenetics, or genetics. Although considered semianadromous, recent studies using otolith geochemistry have revealed life history variation within the critically endangered Delta Smelt. Broadly categorizable as migratory or freshwater residents, we examined Restriction site Associated DNA sequencing data to test for a relationship between genetic variation and migratory behaviors. As previously shown, we found no evidence for neutral population genetic structure within Delta Smelt; however, we found significant evidence for associations between genetic variants and life history phenotypes. Furthermore, discriminant analysis of principal components, hierarchical clustering, and machine learning resulted in accurate assignment of fish into the freshwater resident or migratory classes based on their genotypes. These results suggest the presence of adaptive genetic variants relating to life history variation within a panmictic population. Mechanisms that may lead to this observation are genotype dependent habitat choice and spatially variable selection, both of which could operate each generation and are not exclusive. Given that the population of cultured Delta Smelt are being used as a refugial population for conservation, as a supply for wild population supplementation, and currently represent the majority of all living individuals of this species, we recommend that the hatchery management strategy consider the frequencies of life history-associated alleles and how to maintain this important aspect of Delta Smelt biological variation while under captive propagation.

Thermal regime during parental sexual maturation, but not during offspring rearing, modulates DNA methylation in brook charr (Salvelinus fontinalis)

Epigenetic inheritance can result in plastic responses to changing environments being faithfully transmitted to offspring. However, it remains unclear how epigenetic mechanisms such as DNA methylation can contribute to multigenerational acclimation and adaptation to environmental stressors. Brook charr (Salvelinus fontinalis), an economically important salmonid, is highly sensitive to thermal stress and is of conservation concern in the context of climate change. We studied the effects of temperature during parental sexual maturation and offspring rearing on whole-genome DNA methylation in brook charr juveniles (fry). Parents were split between warm and cold temperatures during sexual maturation, mated in controlled breeding designs, then offspring from each family were split between warm (8°C) and cold (5°C) rearing environments. Using whole-genome bisulfite sequencing, we found 188 differentially methylated regions (DMRs) due to parental maturation temperature after controlling for family structure. By contrast, offspring rearing temperature had a negligible effect on offspring methylation. Stable intergenerational inheritance of DNA methylation and minimal plasticity in progeny could result in the transmission of acclimatory epigenetic states to offspring, priming them for a warming environment. Our findings have implications pertaining to the role of intergenerational epigenetic inheritance in response to ongoing climate change.

Revelation of candidate genes and molecular mechanism of reproductive seasonality in female rohu (Labeo rohita Ham.) by RNA sequencing

BACKGROUND

Carp fish, rohu (Labeo rohita Ham.) is important freshwater aquaculture species of South-East Asia having seasonal reproductive rhythm. There is no holistic study at transcriptome level revealing key candidate genes involved in such circannual rhythm regulated by biological clock genes (BCGs). Seasonality manifestation has two contrasting phases of reproduction, i.e., post-spawning resting and initiation of gonadal activity appropriate for revealing the associated candidate genes. It can be deciphered by RNA sequencing of tissues involved in BPGL (Brain-Pituitary-Gonad-Liver) axis controlling seasonality. How far such BCGs of this fish are evolutionarily conserved across different phyla is unknown. Such study can be of further use to enhance fish productivity as seasonality restricts seed production beyond monsoon season.

RESULT

A total of ~ 150 Gb of transcriptomic data of four tissues viz., BPGL were generated using Illumina TruSeq. De-novo assembled BPGL tissues revealed 75,554 differentially expressed transcripts, 115,534 SSRs, 65,584 SNPs, 514 pathways, 5379 transcription factors, 187 mature miRNA which regulates candidate genes represented by 1576 differentially expressed transcripts are available in the form of web-genomic resources. Findings were validated by qPCR. This is the first report in carp fish having 32 BCGs, found widely conserved in fish, amphibian, reptile, birds, prototheria, marsupials and placental mammals. This is due to universal mechanism of rhythmicity in response to environment and earth rotation having adaptive and reproductive significance.

CONCLUSION

This study elucidates evolutionary conserved mechanism of photo-periodism sensing, neuroendocrine secretion, metabolism and yolk synthesis in liver, gonadal maturation, muscular growth with sensory and auditory perception in this fish. Study reveals fish as a good model for research on biological clock besides its relevance in reproductive efficiency enhancement.

The heat shock response shows plasticity in embryonic lake whitefish (Coregonus clupeaformis) exposed to repeated thermal stress

We examined the impact of repeated thermal stress on the heat shock response (HSR) of thermally sensitive lake whitefish (Coregonus clupeaformis) embryos. Our treatments were designed to mimic temperature fluctuations in the vicinity of industrial thermal effluents. Embryos were either maintained at control temperatures (3 ^oC) or exposed to a repeated thermal stress (TS) of 3 or 6 ^oC above control temperature every 3 or 6 days throughout embryonic development. At 82 days post-fertilisation, repeated TS treatments were stopped and embryos received either a high level TS of 12, 15, or 18 ^oC above ambient temperature for 1 or 4 h, or no additional TS. These treatments were carried out after a 6 h recovery from the last repeated TS. Embryos in the no repeated TS group responded, as expected, with increases in hsp70 mRNA in response to 12, 15 and 18 ^oC high-level TS. However, exposure to repeated TS of 3 or 6 ⁰C every 6 days also resulted in a significant upregulation of hsp70 mRNA relative to the controls. Importantly, these repeated TS events and the associated elevations in hsp70 attenuated the upregulation of hsp70 in response to a 1 h, high-level TS of 12 ^oC above ambient, but not to either longer (4 h) or higher (15 or 18 ^oC) TS events. Conversely, hsp90α mRNA levels were not consistently elevated in the no repeated TS groups exposed to high-level TS. In some instances, hsp90α levels appeared to decrease in embryos exposed to repeated TS followed by a high-level TS. The observed attenuation of the HSR in lake whitefish embryos demonstrates that embryos of this species have plasticity in their HSR and repeated TS may protect against high-level TS, but the response differs based on repeated TS treatment, high-level TS temperature and duration, and the gene of interest.

Mathematical modeling and analysis of the heat shock protein response during thermal stress in fish and HeLa cells

The climate change has the potential to dramatically affect species' thermal physiology and may change the ecology and evolution of species' lineages. In this work, we investigated the transition of dynamics in the heat shock response when the thermal stress approaches the upper thermal limits of species to understand how the climate change may affect the heat shock responses in ectotherms and endotherms. The heat shock protein 70, HSP70, prevents protein denaturation or misfolding under thermal stresses. When thermal stress increases, the number of misfolded proteins increases, which leads to high levels of HSP70 protein. However, when temperatures approach limits of thermal tolerance (i.e., the critical thermal maximum, CTmax, for ectotherms and the superior critical temperature, SCT, for endotherms), levels of HSP70 protein synthesis start to decrease. Thus, we hypothesized that the temperature at the first reduction of HSP abundance indicates the thermal limits of the species. In this work, we provide a mathematical model to investigate the behavior of the heat shock responses related to HSP70 protein. This model captures the dynamics of HSP70 protein and Hsp70 mRNA, in HeLa cells (i.e., representative for endotherms) and multiple species of fishes (i.e., representative for ectotherms) with different acclimation histories. Based on our hypothesis of the relationship between the HSP70 protein level and CTmax/SCT, our model provides three methods to predict the CTmax of fishes with varying acclimation temperature and the SCT of HeLa cells. The CTmax increases as the acclimation temperature increases in fishes, however the CTmax plateaus when the acclimation temperature is higher than 20°C in brook trout, a representative cool water salmonid. Our model also captures the situation that the heat shock reaction in fish is more sensitive to the heat shock temperature than HeLa cells, when the heat shock temperature is lower than the upper thermal tolerance. However, both fish and HeLa cells are sensitive to the heat shock temperature when the temperature reaches the thermal tolerance limits. Additionally, our sensitive analysis result indicates that the status of some components in the heat shock reaction changes when the temperature reaches the thermal tolerance resulting in failure in protein refolding in fish and speeding up the refolding process in HeLa cells. Mathematical analysis is also applied on a simplified mathematical model to investigate the detailed dynamics of the model, such as the steady states of the substrate, Hsp70 mRNA, and HSP70 protein, at different thermal stresses. The comparison between the original model and simplified model shows that the inhibition on HSP70 protein transcription by thermal stresses leads to the reduction in HSP70 protein at extreme temperatures.

Molecular and physiological responses predict acclimation limits in juvenile brook trout (Salvelinus fontinalis)

Understanding the resilience of ectotherms to high temperatures is essential because of the influence of climate change on aquatic ecosystems. The ability of species to acclimate to high temperatures may determine whether populations can persist in their native ranges. We examined physiological and molecular responses of juvenile brook trout (Salvelinus fontinalis) to six acclimation temperatures (5, 10, 15, 20, 23 and 25°C) that span the thermal distribution of the species to predict acclimation limits. Brook trout exhibited an upregulation of stress-related mRNA transcripts (heat shock protein 90-beta, heat shock cognate 71 kDa protein, glutathione peroxidase 1) and downregulation of transcription factors and osmoregulation-related transcripts (nuclear protein 1, Na+/K+/2Cl- co-transporter-1-a) at temperatures ≥20°C. We then examined the effects of acclimation temperature on metabolic rate (MR) and physiological parameters in fish exposed to an acute exhaustive exercise and air exposure stress. Fish acclimated to temperatures ≥20°C exhibited elevated plasma cortisol and glucose, and muscle lactate after exposure to the acute stress. Fish exhibited longer MR recovery times at 15 and 20°C compared with the 5 and 10°C groups; however, cortisol levels remained elevated at temperatures ≥20°C after 24 h. Oxygen consumption in fish acclimated to 23°C recovered quickest after exposure to acute stress. Standard MR was highest and factorial aerobic scope was lowest for fish held at temperatures ≥20°C. Our findings demonstrate how molecular and physiological responses predict acclimation limits in a freshwater fish as the brook trout in the present study had a limited ability to acclimate to temperatures beyond 20°C.

Full-length transcriptome sequences of ridgetail white prawn Exopalaemon carinicauda provide insight into gene expression dynamics during thermal stress

Marine heat waves and extreme high temperature become more frequent and intense in these years, which affected the survival of aquaculture animals. The ridgetail white prawn Exopalaemon carinicauda is an important economic species in eastern China, which has remarkable thermal tolerance. However, there has been little study of its thermal-adaptation mechanisms due to the complex genetic structure and unknown genome. To better understand the molecular mechanisms of E. carinicauda to adapt to the changing temperature, a combination of Illumina-based short reads RNA-seq and single molecule real-time-based full-length transcriptome sequencing was used in this study. In total, 17,212 unigenes from high-quality transcripts of E. carinicauda were generated and 14,663 complete ORFs were detected with an average length of 1980 bp. In addition, the transcriptome profiles of E. carinicauda treated with 34 °C heat stress for 6 and 24 h were analyzed. These differentially expressed genes were primarily enriched in oxidation-reduction process (Gene Ontology enrichment, GO) and the pathways of starch and sucrose metabolism (Kyoto Encyclopedia of Genes and Genomes enrichment, KEGG) after 6 h thermal stress, which indicated that E. carinicauda was suffering the attack by reactive oxygen species. After 24 h thermal stress, these differentially expressed genes were enriched in the pathway of lysosome, glycine, serine and threonine metabolism, fatty acid metabolism (KEGG), which indicated the oxidative stress was decreased. Interestingly, 40 genes for hemocyanin were found to be downregulated after 6 h heat stress, which indicated that the immunocompetence of E. carinicauda decreased after short term thermal stress (6 h). After 24 h thermal stress, E. carinicauda showed transcriptional adaptation to high temperature by upregulating of 11 genes encoding molecular chaperones, including HSP40 and HSP90 which were firstly reported to be related to thermal stress in E. carinicauda. These results promote a better understanding of the thermal-adaptation mechanism of E. carinicauda.

Effects of acclimation temperature on the thermal physiology in two geographically distinct populations of lake sturgeon (Acipenser fulvescens)

Temperature is one of the most important abiotic factors regulating development and biological processes in ectotherms. By 2050, climate change may result in temperature increases of 2.1-3.4°C in Manitoba, Canada. Lake sturgeon, Acipenser fulvescens, from both northern and southern populations in Manitoba were acclimated to 16, 20 and 24°C for 30 days, after which critical thermal maximum (CTmax) trials were conducted to investigate their thermal plasticity. We also examined the effects of temperature on morphological and physiological indices. Acclimation temperature significantly influenced the CTmax, body mass, hepatosomatic index, metabolic rate and the mRNA expression of transcripts involved in the cellular response to heat shock and hypoxia (HSP70, HSP90a, HSP90b, HIF-1α) in the gill of lake sturgeon. Population significantly affected the above phenotypes, as well as the mRNA expression of Na+/K+ ATPase-α1 and the hepatic glutathione peroxidase enzyme activity. The southern population had an average CTmax that was 0.71 and 0.45°C higher than the northern population at 20 and 24°C, respectively. Immediately following CTmax trials, mRNA expression of HSP90a and HIF-1α was positively correlated with individual CTmax of lake sturgeon across acclimation treatments and populations (r = 0.7, r = 0.62, respectively; P < 0.0001). Lake sturgeon acclimated to 20 and 24°C had decreased hepatosomatic indices (93 and 244% reduction, respectively; P < 0.0001) and metabolic suppression (27.7 and 42.1% reduction, respectively; P < 0.05) when compared to sturgeon acclimated to 16°C, regardless of population. Glutathione peroxidase activity and mRNA expression Na+/K+ ATPase-α1 were elevated in the northern relative to the southern population. Acclimation to 24°C also induced mortality in both populations when compared to sturgeon acclimated to 16 and 20°C. Thus, increased temperatures have wide-ranging population-specific physiological consequences for lake sturgeon across biological levels of organization.

Sub-lethal temperature thresholds indicate acclimation and physiological limits in brook trout Salvelinus fontinalis

The upper thermal tolerance of brook trout Salvelinus fontinalis was estimated using critical thermal maxima (CT_max ) experiments on fish acclimated to temperatures that span the species' thermal range (5-25°C). The CT_max increased with acclimation temperature but plateaued in fish acclimated to 20, 23 and 25°C. Plasma lactate was highest, and the hepato-somatic index (I_H ) was lowest at 23 and 25°C, which suggests additional metabolic costs at those acclimation temperatures. The results suggest that there is a sub-lethal threshold between 20 and 23°C, beyond which the fish experience reduced physiological performance.

Adaptation of plasticity to projected maximum temperatures and across climatically defined bioregions

Resilience to environmental stressors due to climate warming is influenced by local adaptations, including plastic responses. The recent literature has focused on genomic signatures of climatic adaptation, but little is known about how plastic capacity may be influenced by biogeographic and evolutionary processes. We investigate phenotypic plasticity as a target of climatic selection, hypothesizing that lineages that evolved in warmer climates will exhibit greater plastic adaptive resilience to upper thermal stress. This was experimentally tested by comparing transcriptomic responses within and among temperate, subtropical, and desert ecotypes of Australian rainbowfish subjected to contemporary and projected summer temperatures. Critical thermal maxima were estimated, and ecological niches delineated using bioclimatic modeling. A comparative phylogenetic expression variance and evolution model was used to assess plastic and evolved changes in gene expression. Although 82% of all expressed genes were found in the three ecotypes, they shared expression patterns in only 5 out of 236 genes that responded to the climate change experiment. A total of 532 genes showed signals of adaptive (i.e., genetic-based) plasticity due to ecotype-specific directional selection, and 23 of those responded to projected summer temperatures. Network analyses demonstrated centrality of these genes in thermal response pathways. The greatest adaptive resilience to upper thermal stress was shown by the subtropical ecotype, followed by the desert and temperate ecotypes. Our findings indicate that vulnerability to climate change will be highly influenced by biogeographic factors, emphasizing the value of integrative assessments of climatic adaptive traits for accurate estimation of population and ecosystem responses.

Thermal tolerance depends on season, age and body condition in imperilled redside dace Clinostomus elongatus

Urbanization tends to increase water temperatures in streams and rivers and is hypothesized to be contributing to declines of many freshwater fishes. However, factors that influence individual variation in thermal tolerance, and how these may change seasonally, are not well understood. To address this knowledge gap, we studied redside dace Clinostomus elongatus, an imperilled stream fish native to rapidly urbanizing areas of eastern North America. In wild redside dace from rural Ohio, USA, acute upper thermal tolerance (i.e. critical thermal maximum, CTmax) ranged between ~34°C in summer (stream temperature ~22°C) and 27°C in winter (stream temperature ~2°C). Juveniles had higher CTmax than adults in spring and summer, but in winter, CTmax was higher in adults. Thermal safety margins (CTmax - ambient water temperature; ~11°C) were less than the increases in peak water temperature predicted for many redside dace streams due to the combined effects of climate change and urbanization. Furthermore, behavioural agitation occurred 5-6°C below CTmax. Safety margins were larger (>20°C) in autumn and winter. In addition, redside dace were more sensitive (2.5°C lower CTmax) than southern redbelly dace Chrosomus erythrogaster, a non-imperilled sympatric cyprinid. Body condition (Fulton's K) of adult redside dace was positively correlated with CTmax, but in juveniles, this relationship was significant only in one of two summers of experiments. Next, we measured CTmax of captive redside dace fed experimentally manipulated diets. In adults, but not juveniles, CTmax was higher in fish fed a high- vs. low-ration diet, indicating a causal link between nutrition and thermal tolerance. We conclude that redside dace will be challenged by predicted future summer temperatures, especially in urbanized habitats. Thus, habitat restoration that mitigates temperature increases is likely to benefit redside dace. We also suggest habitat restoration that improves food availability may increase thermal tolerance, and thus population resilience.

Deep transcriptome analysis of the heat shock response in an Atlantic sturgeon (Acipenser oxyrinchus) cell line

Despite efforts to restore Atlantic sturgeon in European rivers, aquaculture techniques result in animals with high post-release mortality due to, among other reasons, their low tolerance to increasing water temperature. Marker genes to monitor heat stress are needed in order to identify heat-resistant fish. Therefore, an Atlantic sturgeon cell line was exposed to different heat shock protocols (30 °C and 35 °C) and differences in gene expression were investigated. In total 3020 contigs (∼1.5%) were differentially expressed. As the core of the upregulated contigs corresponded to heat shock proteins (HSP), the heat shock factor (HSF) and the HSP gene families were annotated in Atlantic sturgeon and mapped via Illumina RNA sequencing to identify heat-inducible family members. Up to 6 hsf and 76 hsp genes were identified in the Atlantic sturgeon transcriptome resources, 16 of which were significantly responsive to the applied heat shock. The previously studied hspa1 (hsp70) gene was only significantly upregulated at the highest heat shock (35 °C), while a set of 5 genes (hspc1, hsph3a, hspb1b, hspb11a, and hspb11b) was upregulated at all conditions. Although the hspc1 (hsp90a) gene was previously used as heat shock-marker in sturgeons, we found that hspb11a is the most heat-inducible gene, with up to 3296-fold higher expression in the treated cells, constituting the candidate gene markers for in vivo trials.

Sensitivities of an endemic, endangered California smelt and two non-native fishes to serial increases in temperature and salinity: implications for shifting community structure with climate change

In many aquatic systems, native fishes are in decline and the factors responsible are often elusive. In the San Francisco Estuary (SFE) in California, interactions among native and non-native species are key factors contributing to the decline in abundance of endemic, endangered Delta Smelt (Hypomesus transpacificus). Climate change and drought-related stressors are further exacerbating declines. To assess how multiple environmental changes affect the physiology of native Delta Smelt and non-native Mississippi Silverside (Menidia beryllina) and Largemouth Bass (Micropterus salmoides), fishes were exposed to serial exposures of a single stressor (elevated temperature or salinity) followed by two stressors (elevated temperature and salinity) to determine how a single stressor affects the capacity to cope with the addition of a second stressor. Critical thermal maximum (CTMax; a measure of upper temperature tolerance) was determined after 0, 2, 4 and 7 days following single and multiple stressors of elevated temperature (16°C vs. 20°C) and salinity (2.4 vs. 8-12 ppt, depending on species). Under control conditions, non-native fishes had significantly higher CTMax than the native Delta Smelt. An initial temperature or salinity stressor did not negatively affect the ability of any species to tolerate a subsequent multiple stressor. While elevated salinity had little effect on CTMax, a 4°C increase in temperature increased CTMax. Bass experienced an additive effect of increased temperature and salinity on CTMax, such that CTMax further increased under multiple stressors. In addition, Bass demonstrated physiological sensitivity to multiple stressors demonstrated by changes in hematocrit and plasma osmolality, whereas the physiology of Silversides remained unaffected. Non-native Bass and Mississippi Silversides showed consistently higher thermal tolerance limits than the native Delta Smelt, supporting their abundance in warmer SFE habitats. Continued increases in SFE water temperatures predicted with climate change may further impact endangered Delta Smelt populations directly if habitat temperatures exceed thermal limits.

Developing specific molecular biomarkers for thermal stress in salmonids

Background

Pacific salmon (Oncorhynchus spp.) serve as good biological indicators of the breadth of climate warming effects on fish because their anadromous life cycle exposes them to environmental challenges in both marine and freshwater environments. Our study sought to mine the extensive functional genomic studies in fishes to identify robust thermally-responsive biomarkers that could monitor molecular physiological signatures of chronic thermal stress in fish using non-lethal sampling of gill tissue.

Results

Candidate thermal stress biomarkers for gill tissue were identified using comparisons among microarray datasets produced in the Molecular Genetics Laboratory, Pacific Biological Station, Nanaimo, BC, six external, published microarray studies on chronic and acute temperature stress in salmon, and a comparison of significant genes across published studies in multiple fishes using deep literature mining. Eighty-two microarray features related to 39 unique gene IDs were selected as candidate chronic thermal stress biomarkers. Most of these genes were identified both in the meta-analysis of salmon microarray data and in the literature mining for thermal stress markers in salmonids and other fishes. Quantitative reverse transcription PCR (qRT-PCR) assays for 32 unique genes with good efficiencies across salmon species were developed, and their activity in response to thermally challenged sockeye salmon (O. nerka) and Chinook salmon (O. tshawytscha) (cool, 13–14 °C and warm temperatures 18–19 °C) over 5–7 days was assessed. Eight genes, including two transcripts of each SERPINH1 and HSP90AA1, FKBP10, MAP3K14, SFRS2, and EEF2 showed strong and robust chronic temperature stress response consistently in the discovery analysis and both sockeye and Chinook salmon validation studies.

Conclusions

The results of both discovery analysis and gene expression showed that a panel of genes involved in chaperoning and protein rescue, oxidative stress, and protein biosynthesis were differentially activated in gill tissue of Pacific salmon in response to elevated temperatures. While individually, some of these biomarkers may also respond to other stressors or biological processes, when expressed in concert, we argue that a biomarker panel comprised of some or all of these genes could provide a reliable means to specifically detect thermal stress in field-caught salmon.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5108-9) contains supplementary material, which is available to authorized users.

Regional differences in thermal adaptation of a cold-water fish Rhynchocypris oxycephalus revealed by thermal tolerance and transcriptomic responses

Understanding how populations adapt to different thermal environments is an important issue for biodiversity conservation in the context of recent global warming. To test the hypothesis that populations from southern region are more sensitive to climate change than northern region in cold-water species, we determined the thermal tolerance of two geographical populations of a cold-water fish, Rhynchocypris oxycephalus: the Hangzhou population from southern region and the Gaizhou population from northern region, then compared their transcriptomic responses between a control and a high temperature treatment. The results showed that the thermal tolerance range and thermal tolerance polygon area of Hangzhou population were narrower than the Gaizhou population, indicating populations from southern region were possibly more vulnerable. Further transcriptomic analysis revealed that the Gaizhou population expressed more temperature responding genes than the Hangzhou population (583 VS. 484), corresponding with their higher thermal tolerance, while some of these genes (e.g. heat shock protein) showed higher expression in the Hangzhou population under control condition, suggesting individuals from southern region possibly have already responded to the present higher environmental temperature pressure. Therefore, these results confirm the prediction that populations from southern region are more sensitive to global warming, and will be important for their future conservation.

The utility of transcriptomics in fish conservation

There is growing recognition of the need to understand the mechanisms underlying organismal resilience (i.e. tolerance, acclimatization) to environmental change to support the conservation management of sensitive and economically important species. Here, we discuss how functional genomics can be used in conservation biology to provide a cellular-level understanding of organismal responses to environmental conditions. In particular, the integration of transcriptomics with physiological and ecological research is increasingly playing an important role in identifying functional physiological thresholds predictive of compensatory responses and detrimental outcomes, transforming the way we can study issues in conservation biology. Notably, with technological advances in RNA sequencing, transcriptome-wide approaches can now be applied to species where no prior genomic sequence information is available to develop species-specific tools and investigate sublethal impacts that can contribute to population declines over generations and undermine prospects for long-term conservation success. Here, we examine the use of transcriptomics as a means of determining organismal responses to environmental stressors and use key study examples of conservation concern in fishes to highlight the added value of transcriptome-wide data to the identification of functional response pathways. Finally, we discuss the gaps between the core science and policy frameworks and how thresholds identified through transcriptomic evaluations provide evidence that can be more readily used by resource managers.

Environmental and genetic determinants of transcriptional plasticity in Chinook salmon

Variation in gene transcription is widely believed to be the mechanistic basis of phenotypically plastic traits; however, comparatively little is known about the inheritance patterns of transcriptional variation that would allow us to predict its response to selection. In addition, acclimation to different environmental conditions influences acute transcriptional responses to stress and it is unclear if these effects are heritable. To address these gaps in knowledge, we assayed levels of messenger RNA for 14 candidate genes at rest and in response to a 24-h confinement stress for 72 half-sib families of Chinook salmon reared in two different environments (hatchery and semi-natural stream channel). We observed extensive plasticity for mRNA levels of metabolic and stress response genes and demonstrated that mRNA level plasticity due to rearing environment affects mRNA level plasticity in response to stress. These effects have important implications for natural populations experiencing multiple stressors. We identified genotype-by-environment interactions for mRNA levels that were dominated by maternal effects; however, mRNA level response to challenge also exhibited a non-additive genetic basis. Our results indicate that while plasticity for mRNA levels can evolve, predicting the outcome of selection will be difficult. The inconsistency in genetic architecture among treatment groups suggests there is considerable cryptic genetic variation for gene expression.

Thermal conditions during early life influence seasonal maternal strategies in the three-spined stickleback

BACKGROUND

Conditions experienced by a female during early life may affect her reproductive strategies and maternal investment later in life. This effect of early environmental conditions is a potentially important mechanism by which animals can compensate for the negative impacts of climate change. In this study, we experimentally tested whether three-spined sticklebacks (Gasterosteus aculeatus) change their maternal strategy according to environmental temperatures experienced earlier in life. We studied maternal investment from a life-history perspective because females are expected to adjust their reproductive strategy in relation to their current and future reproductive returns as well as offspring fitness.

RESULTS

F1 families were reared in control and elevated winter temperatures and their reproductive trajectories were studied when returned to common conditions. Females that had experienced the warm winter treatment (n = 141) had a lower fecundity and reduced breeding and total lifespan compared to the control individuals (n = 159). Whereas the control females tended to produce their heaviest and largest clutches in their first reproductive attempt, the warm-acclimated females invested less in their first clutch, but then produced increasingly heavy clutches over the course of the breeding season. Egg mass increased with clutch number at a similar rate in the two groups. The warm-acclimated females increased the investment of carotenoids in the first and last clutches of the season. Thus, any transgenerational effects of the maternal thermal environment on offspring phenotype may be mediated by the allocation of antioxidants into eggs but not by egg size.

CONCLUSIONS

Our results indicate that conditions experienced by females during juvenile life have a profound effect on life-time maternal reproductive strategies. The temperature-induced changes in maternal strategy may be due to constraints imposed by the higher energetic costs of a warm environment, but it is possible that they allow the offspring to compensate for higher energetic costs and damage when they face the same thermal stress as did their mothers.

Temperature effects on gene expression and morphological development of European eel, Anguilla anguilla larvae

Temperature is important for optimization of rearing conditions in aquaculture, especially during the critical early life history stages of fish. Here, we experimentally investigated the impact of temperature (16, 18, 20, 22 and 24°C) on thermally induced phenotypic variability, from larval hatch to first-feeding, and the linked expression of targeted genes [heat shock proteins (hsp), growth hormone (gh) and insulin-like growth factors (igf)] associated to larval performance of European eel, Anguilla anguilla. Temperature effects on larval morphology and gene expression were investigated throughout early larval development (in real time from 0 to 18 days post hatch) and at specific developmental stages (hatch, jaw/teeth formation, and first-feeding). Results showed that hatch success, yolk utilization efficiency, survival, deformities, yolk utilization, and growth rates were all significantly affected by temperature. In real time, increasing temperature from 16 to 22°C accelerated larval development, while larval gene expression patterns (hsp70, hsp90, gh and igf-1) were delayed at cold temperatures (16°C) or accelerated at warm temperatures (20–22°C). All targeted genes (hsp70, hsp90, gh, igf-1, igf-2a, igf-2b) were differentially expressed during larval development. Moreover, expression of gh was highest at 16°C during the jaw/teeth formation, and the first-feeding developmental stages, while expression of hsp90 was highest at 22°C, suggesting thermal stress. Furthermore, 24°C was shown to be deleterious (resulting in 100% mortality), while 16°C and 22°C (~50 and 90% deformities respectively) represent the lower and upper thermal tolerance limits. In conclusion, the high survival, lowest incidence of deformities at hatch, high yolk utilization efficiency, high gh and low hsp expression, suggest 18°C as the optimal temperature for offspring of European eel. Furthermore, our results suggest that the still enigmatic early life history stages of European eel may inhabit the deeper layer of the Sargasso Sea and indicate vulnerability of this critically endangered species to increasing ocean temperature.

Transcriptional mechanisms underlying life-history responses to climate change in the three-spined stickleback

Phenotypic plasticity, the ability of an organism to express different phenotypes depending on the environment, provides an important mechanism by which an animal population can persist under rapid climate change. We experimentally tested both life‐history and transcriptional responses of an ecological model species, the three‐spined stickleback, to warm acclimation at the southern edge of its European range. We explored cross‐environment genetic correlations of key life‐history traits in male sticklebacks exposed to long‐term temperature changes to examine whether the plasticity pattern was variable among genotypes by using a character‐state approach. We also studied gene expression plasticity by analysing both whole‐transcriptome and candidate gene expression in brain and liver. Male sticklebacks that developed under warmer conditions during winter were smaller in size and invested less in nuptial coloration at the beginning of the breeding season, showing similar responses across different genotypes. The lack of genetic variation in life‐history responses may limit any future evolution of the thermal reaction norm in the study population. After long‐term exposure to increased winter temperatures, genes responsible for several metabolic and oxidation–reduction processes were upregulated, and some hormone genes involved in growth and reproduction were downregulated in the brain. In the liver, there was no significantly represented gene ontology by the differentially expressed genes. Since a higher temperature leads to a higher resting metabolic rate, living in warmer environments may incur higher energetic costs for ectotherms to maintain cellular homoeostasis, resulting in negative consequences for life‐history traits. The expression of genes related to metabolism, cellular homoeostasis and regulatory signalling may underlie temperature‐induced changes in life history.

Plasticity in gene transcription explains the differential performance of two invasive fish species

Phenotypic plasticity buffers organisms from environmental change and is hypothesized to aid the initial establishment of nonindigenous species in novel environments and postestablishment range expansion. The genetic mechanisms that underpin phenotypically plastic traits are generally poorly characterized; however, there is strong evidence that modulation of gene transcription is an important component of these responses. Here, we use RNA sequencing to examine the transcriptional basis of temperature tolerance for round and tubenose goby, two nonindigenous fish species that differ dramatically in the extent of their Great Lakes invasions despite similar invasion dates. We used generalized linear models of read count data to compare gene transcription responses of organisms exposed to increased and decreased water temperature from those at ambient conditions. We identify greater response in the magnitude of transcriptional changes for the more successful round goby compared with the less successful tubenose goby. Round goby transcriptional responses reflect alteration of biological function consistent with adaptive responses to maintain or regain homeostatic function in other species. In contrast, tubenose goby transcription patterns indicate a response to stressful conditions, but the pattern of change in biological functions does not match those expected for a return to homeostatic status. Transcriptional plasticity plays an important role in the acute thermal tolerance for these species; however, the impaired response to stress we demonstrate in the tubenose goby may contribute to their limited invasion success relative to the round goby. Transcriptional profiling allows the simultaneous assessment of the magnitude of transcriptional response as well as the biological functions involved in the response to environmental stress and is thus a valuable approach for evaluating invasion potential.

Sublethal salinity stress contributes to habitat limitation in an endangered estuarine fish

As global change alters multiple environmental conditions, predicting species' responses can be challenging without understanding how each environmental factor influences organismal performance. Approaches quantifying mechanistic relationships can greatly complement correlative field data, strengthening our abilities to forecast global change impacts. Substantial salinity increases are projected in the San Francisco Estuary, California, due to anthropogenic water diversion and climatic changes, where the critically endangered delta smelt (Hypomesus transpacificus) largely occurs in a low-salinity zone (LSZ), despite their ability to tolerate a much broader salinity range. In this study, we combined molecular and organismal measures to quantify the physiological mechanisms and sublethal responses involved in coping with salinity changes. Delta smelt utilize a suite of conserved molecular mechanisms to rapidly adjust their osmoregulatory physiology in response to salinity changes in estuarine environments. However, these responses can be energetically expensive, and delta smelt body condition was reduced at high salinities. Thus, acclimating to salinities outside the LSZ could impose energetic costs that constrain delta smelt's ability to exploit these habitats. By integrating data across biological levels, we provide key insight into the mechanistic relationships contributing to phenotypic plasticity and distribution limitations and advance the understanding of the molecular osmoregulatory responses in nonmodel estuarine fishes.

Effects of high temperatures on threatened estuarine fishes during periods of extreme drought

Climate change and associated increases in water temperatures may impact physiological performance in ectotherms and exacerbate endangered species declines. We used an integrative approach to assess the impact of elevated water temperature on two fishes of immediate conservation concern in a large estuary system, the threatened longfin smelt (Spirinchus thaleichthys) and endangered delta smelt (Hypomesus transpacificus). Abundances have reached record lows in California, USA, and these populations are at imminent risk of extirpation. California is currently impacted by a severe drought, resulting in high water temperatures, conditions that will become more common as a result of climate change. We exposed fish to environmentally relevant temperatures (14°C and 20°C) and used RNA sequencing to examine the transcriptome-wide responses to elevated water temperature in both species. Consistent with having a lower temperature tolerance, longfin smelt exhibited a pronounced cellular stress response, with an upregulation of heat shock proteins, after exposure to 20°C that was not observed in delta smelt. We detected an increase in metabolic rate in delta smelt at 20°C and increased expression of genes involved in metabolic processes and protein synthesis, patterns not observed in longfin smelt. Through examination of responses across multiple levels of biological organization, and by linking these responses to habitat distributions in the wild, we demonstrate that longfin smelt may be more susceptible than delta smelt to increases in temperatures, and they have little room to tolerate future warming in California. Understanding the species-specific physiological responses of sensitive species to environmental stressors is crucial for conservation efforts and managing aquatic systems globally.

Assessments at multiple levels of biological organization allow for an integrative determination of physiological tolerances to turbidity in an endangered fish species

Turbidity can influence trophic levels by altering species composition and can potentially affect fish feeding strategies and predator-prey interactions. The estuarine turbidity maximum, described as an area of increased suspended particles, phytoplankton and zooplankton, generally represents a zone with higher turbidity and enhanced food sources important for successful feeding and growth in many fish species. The delta smelt (Hypomesus transpacificus) is an endangered, pelagic fish species endemic to the San Francisco Estuary and Sacramento-San Joaquin River Delta, USA, where it is associated with turbid waters. Turbidity is known to play an important role for the completion of the species' life cycle; however, turbidity ranges in the Delta are broad, and specific requirements for this fish species are still unknown. To evaluate turbidity requirements for early life stages, late-larval delta smelt were maintained at environmentally relevant turbidity levels ranging from 5 to 250 nephelometric turbidity units (NTU) for 24 h, after which a combination of physiological endpoints (molecular biomarkers and cortisol), behavioural indices (feeding) and whole-organism measures (survival) were determined. All endpoints delivered consistent results and identified turbidities between 25 and 80 NTU as preferential. Delta smelt survival rates were highest between 12 and 80 NTU and feeding rates were highest between 25 and 80 NTU. Cortisol levels indicated minimal stress between 35 and 80 NTU and were elevated at low turbidities (5, 12 and 25 NTU). Expression of stress-related genes indicated significant responses for gst, hsp70 and glut2 in high turbidities (250 NTU), and principal component analysis on all measured genes revealed a clustering of 25, 35, 50 and 80 NTU separating the medium-turbidity treatments from low- and high-turbidity treatments. Taken together, these data demonstrate that turbidity levels that are either too low or too high affect delta smelt physiological performance, causing significant effects on overall stress, food intake and mortality. They also highlight the need for turbidity to be considered in habitat and water management decisions.

Coupled Downscaled Climate Models and Ecophysiological Metrics Forecast Habitat Compression for an Endangered Estuarine Fish

Climate change is driving rapid changes in environmental conditions and affecting population and species’ persistence across spatial and temporal scales. Integrating climate change assessments into biological resource management, such as conserving endangered species, is a substantial challenge, partly due to a mismatch between global climate forecasts and local or regional conservation planning. Here, we demonstrate how outputs of global climate change models can be downscaled to the watershed scale, and then coupled with ecophysiological metrics to assess climate change effects on organisms of conservation concern. We employed models to estimate future water temperatures (2010–2099) under several climate change scenarios within the large heterogeneous San Francisco Estuary. We then assessed the warming effects on the endangered, endemic Delta Smelt, Hypomesus transpacificus, by integrating localized projected water temperatures with thermal sensitivity metrics (tolerance, spawning and maturation windows, and sublethal stress thresholds) across life stages. Lethal temperatures occurred under several scenarios, but sublethal effects resulting from chronic stressful temperatures were more common across the estuary (median >60 days above threshold for >50% locations by the end of the century). Behavioral avoidance of such stressful temperatures would make a large portion of the potential range of Delta Smelt unavailable during the summer and fall. Since Delta Smelt are not likely to migrate to other estuaries, these changes are likely to result in substantial habitat compression. Additionally, the Delta Smelt maturation window was shortened by 18–85 days, revealing cumulative effects of stressful summer and fall temperatures with early initiation of spring spawning that may negatively impact fitness. Our findings highlight the value of integrating sublethal thresholds, life history, and in situ thermal heterogeneity into global change impact assessments. As downscaled climate models are becoming widely available, we conclude that similar assessments at management-relevant scales will improve the scientific basis for resource management decisions.