Identification of genes associated with heat tolerance in Arctic charr exposed to acute thermal stress

Tissue-specific transcriptional response of post-larval clownfish to ocean warming

Anthropogenically driven climate change is predicted to increase average sea surface temperatures, as well as the frequency and intensity of marine heatwaves in the future. This increasing temperature is predicted to have a range of negative physiological impacts on multiple life-stages of coral reef fish. Nevertheless, studies of early-life stages remain limited, and tissue-specific transcriptomic studies of post-larval coral reef fish are yet to be conducted. Here, in an aquaria-based study we investigate the tissue-specific (brain, liver, muscle, and digestive tract) transcriptomic response of post-larval (20 dph) Amphiprion ocellaris to temperatures associated with future climate change (+3 °C). Additionally, we utilized metatranscriptomic sequencing to investigate how the microbiome of the digestive tract changes at +3 °C. Our results show that the transcriptional response to elevated temperatures is highly tissue-specific, as the number of differentially expressed genes (DEGs) and gene functions varied amongst the brain (102), liver (1785), digestive tract (380), and muscle (447). All tissues displayed DEGs associated with thermal stress, as 23 heat-shock protein genes were upregulated in all tissues. Our results indicate that post-larval clownfish may experience liver fibrosis-like symptoms at +3 °C as genes associated with extracellular matrix structure, oxidative stress, inflammation, glucose transport, and metabolism were all upregulated. We also observe a shift in the digestive tract microbiome community structure, as Vibrio sp. replace Escherichia coli as the dominant bacteria. This shift is coupled with the dysregulation of various genes involved in immune response in the digestive tract. Overall, this study highlights post-larval clownfish will display tissue-specific transcriptomic responses to future increases in temperature, with many potentially harmful pathways activated at +3 °C.

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

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

Transcriptome analysis reveals genes connected to temperature adaptation in juvenile antarctic krill Euphausia superba

BACKGROUND

The Antarctic krill, Euphausia superba (E. superba), is a key organism in the Antarctic marine ecosystem and has been widely studied. However, there is a lack of transcriptome data focusing on temperature responses.

METHODS

In this study, we performed transcriptome sequencing of E. superba samples exposed to three different temperatures: -1.19 °C (low temperature, LT), - 0.37 °C (medium temperature, MT), and 3 °C (high temperature, HT).

RESULTS

Illumina sequencing generated 772,109,224 clean reads from the three temperature groups. In total, 1,623, 142, and 842 genes were differentially expressed in MT versus LT, HT versus LT, and HT versus MT, respectively. Moreover, Kyoto Encyclopedia of Genes and Genomes analysis revealed that these differentially expressed genes were mainly involved in the Hippo signaling pathway, MAPK signaling pathway, and Toll-like receptor signaling pathway. Quantitative reverse-transcription PCR revealed that ESG037073 expression was significantly upregulated in the MT group compared with the LT group, and ESG037998 expression was significantly higher in the HT group than in the LT group.

CONCLUSIONS

This is the first transcriptome analysis of E. superba exposed to three different temperatures. Our results provide valuable resources for further studies on the molecular mechanisms underlying temperature adaptation in E. superba.

Effects of juvenile thermal preconditioning on the heat-shock, immune, and stress responses of rainbow trout upon a secondary thermal challenge

Higher water temperatures and pathogens are both significant factors that negatively affect the welfare of teleost fish. In aquaculture, compared to natural populations, these problems are especially exacerbated, as the animals have relatively limited mobility, and the higher density promotes faster spread of infectious diseases. Because of the potential harm these stressors can inflict, methods that can limit the damage of these stressors are particularly valuable. As a method of interest, early-life thermal preconditioning of animals demonstrated some potential for effective improvements in thermotolerance. However, the potential effects of the method on the immune system via the heat-stress model have not been explored. In this experiment, juvenile-stage thermal preconditioned rainbow trout (Oncorhynchus mykiss) were subjected to a secondary thermal challenge, animals were collected and sampled at the time of lost equilibrium. The effects of preconditioning on the general stress response was assessed by measuring the plasma cortisol levels. In addition, we also examined hsp70 and hsc70 mRNA levels in the spleen and gill tissues, as well as IL-1β, IL-6, TNF-α, IFN-1, β2m, and MH class I transcripts via qRT-PCR. No changes in CTmax were observed between the preconditioned and control cohorts upon the second challenge. IL-1β and IL-6 transcripts were generally upregulated with increased temperature of the secondary thermal challenge, whereas IFN-1 transcripts were upregulated in the spleen, but downregulated in the gills, along with MH class I. The juvenile thermal preconditioning produced a series of changes in transcript levels for IL-1β, TNF-α, IFN-1, and hsp70 but the dynamics of these differences were inconsistent. Finally, analysis of plasma cortisol levels presented significantly lower cortisol levels in the pre-conditioned animals compared to the non-pre-conditioned control cohort.

Koala Genome Survey: An Open Data Resource to Improve Conservation Planning

Genome sequencing is a powerful tool that can inform the management of threatened species. Koalas (Phascolarctos cinereus) are a globally recognized species that captured the hearts and minds of the world during the 2019/2020 Australian megafires. In 2022, koalas were listed as ‘Endangered’ in Queensland, New South Wales, and the Australian Capital Territory. Populations have declined because of various threats such as land clearing, habitat fragmentation, and disease, all of which are exacerbated by climate change. Here, we present the Koala Genome Survey, an open data resource that was developed after the Australian megafires. A systematic review conducted in 2020 demonstrated that our understanding of genomic diversity within koala populations was scant, with only a handful of SNP studies conducted. Interrogating data showed that only 6 of 49 New South Wales areas of regional koala significance had meaningful genome-wide data, with only 7 locations in Queensland with SNP data and 4 locations in Victoria. In 2021, we launched the Koala Genome Survey to generate resequenced genomes across the Australian east coast. We have publicly released 430 koala genomes (average coverage: 32.25X, range: 11.3–66.8X) on the Amazon Web Services Open Data platform to accelerate research that can inform current and future conservation planning.

Genetic and epigenetic regulation of growth, reproduction, disease resistance and stress responses in aquaculture

Major progress has been made with genomic and genetic studies in aquaculture in the last decade. However, research on epigenetic regulation of aquaculture traits is still at an early stage. It is apparent that most, if not all, aquaculture traits are regulated at both genetic and epigenetic levels. This paper reviews recent progress in understanding of genetic and epigenetic regulation of important aquaculture traits such as growth, reproduction, disease resistance, and stress responses. Although it is challenging to make generalized statements, DNA methylation is mostly correlated with down-regulation of gene expression, especially when at promoters and enhancers. As such, methylation of growth factors and their receptors is negatively correlated with growth; hypomethylation of genes important for stress tolerance is correlated with increased stress tolerance; hypomethylation of genes important for male or female sex differentiation leads to sex differentiation into males or females, respectively. It is apparent that environmental regulation of aquaculture traits is mediated at the level of epigenetic regulation, and such environment-induced epigenetic changes appeared to be intergenerationally inherited, but evidences for transgenerational inheritance are still limited.

Genetic approaches for increasing fitness in endangered species

The global rate of wildlife extinctions is accelerating, and the persistence of many species requires conservation breeding programs. A central paradigm of these programs is to preserve the genetic diversity of the founder populations. However, this may preserve original characteristics that make them vulnerable to extinction. We introduce targeted genetic intervention (TGI) as an alternative approach that promotes traits that enable species to persist in the face of threats by changing the incidence of alleles that impact on fitness. The TGI toolkit includes methods with established efficacy in model organisms and agriculture but are largely untried for conservation, such as synthetic biology and artificial selection. We explore TGI approaches as a species-restoration tool for intractable threats including infectious disease and climate change.

The Atlantic salmon's stress- and immune-related transcriptional responses to moderate hypoxia, an incremental temperature increase, and these challenges combined

The marine environment is predicted to become warmer, and more hypoxic, and these conditions may negatively impact the health and survival of coastal fish species, including wild and farmed Atlantic salmon (Salmo salar). Thus, we examined how: (1) moderate hypoxia (∼70% air saturation) at 12°C for 3 weeks; (2) an incremental temperature increase from 12°C to 20°C (at 1°C week-1) followed by 4 weeks at 20°C; and (3) treatment "2" combined with moderate hypoxia affected transcript expression in the liver of post-smolts as compared to control conditions (normoxia, 12°C). Specifically, we assessed the expression of 45 genes related to the heat shock response, oxidative stress, apoptosis, metabolism and immunity using a high-throughput qPCR approach (Fluidigm Biomark™ HD). The expression profiles of 27 "stress"-related genes indicated that: (i) moderate hypoxia affected the expression of several stress genes at 12°C; (ii) their expression was impacted by 16°C under normoxic conditions, and this effect increased until 20°C; (iii) the effects of moderate hypoxia were not additive to those at temperatures above 16°C; and (iv) long-term (4 weeks) exposure to 20°C, with or without hypoxia, resulted in a limited acclimatory response. In contrast, the expression of 15 immune-related genes was not greatly affected until temperatures reached 20°C, and this effect was particularly evident in fish exposed to the added challenge of hypoxia. These results provide valuable information on how these two important environmental factors affect the "stress" physiology and immunology of Atlantic salmon, and we identify genes that may be useful as hypoxia and/or temperature biomarkers in salmonids and other fishes.

Intraspecific variation in tolerance of warming in fishes

Intraspecific variation in key traits such as tolerance of warming can have profound effects on ecological and evolutionary processes, notably responses to climate change. The empirical evidence for three primary elements of intraspecific variation in tolerance of warming in fishes is reviewed. The first is purely mechanistic that tolerance varies across life stages and as fishes become mature. The limited evidence indicates strongly that this is the case, possibly because of universal physiological principles. The second is intraspecific variation that is because of phenotypic plasticity, also a mechanistic phenomenon that buffers individuals' sensitivity to negative impacts of global warming in their lifetime, or to some extent through epigenetic effects over successive generations. Although the evidence for plasticity in tolerance to warming is extensive, more work is required to understand underlying mechanisms and to reveal whether there are general patterns. The third element is intraspecific variation based on heritable genetic differences in tolerance, which underlies local adaptation and may define long-term adaptability of a species in the face of ongoing global change. There is clear evidence of local adaptation and some evidence of heritability of tolerance to warming, but the knowledge base is limited with detailed information for only a few model or emblematic species. There is also strong evidence of structured variation in tolerance of warming within species, which may have ecological and evolutionary significance irrespective of whether it reflects plasticity or adaptation. Although the overwhelming consensus is that having broader intraspecific variation in tolerance should reduce species vulnerability to impacts of global warming, there are no sufficient data on fishes to provide insights into particular mechanisms by which this may occur.

Transcriptomic responses to heat stress in gill and liver of endangered Brachymystax lenok tsinlingensis

Global warming significantly affects fish, particularly cold-water fish, because increased temperature adversely impacts their abilities to grow or reproduce, and eventually influences their fitness or even causes death. To survive, fish may alter their distribution or behavior to avoid the stress, and perhaps acclimate or evolve resistance to the elevated temperature. Brachymystax lenok tsinlingensis is an endangered cold-water species in China, and it has been found to alter the altitudinal distribution, decrease swimming efficiency and develop resistance under heat exposure, which badly impact the continuing conservation work. To better protect them, it is essential to understand how they respond to thermal stress behaviorally and physiologically. Therefore, the fish were exposed to 24.5 °C and based on the time taken for them to lose equilibrium, they were separately sampled as sensitive and tolerant groups. Both gill and liver tissues were collected from both groups for transcriptome sequencing. Sequencing results demonstrated that control and tolerant groups were similar in transcriptomic patterns and sensitive groups differentially expressed more genes than tolerant ones, suggesting the gene expression of tolerant groups may return to base levels as exposure time increased. Tissue differences were the major factor affecting gene expression, and they also displayed different physiological responses to heat stress. Consistent with other studies, heat shock response, immune response, metabolic adjustment and ion transport were found to be triggered after exposed to elevated temperature. The findings would contribute to a better understanding of responding mechanisms of fish to thermal stress and provide guidance for future conservation programs.

Temperature training improves transcriptional homeostasis after heat shock in juvenile Atlantic sturgeon (Acipenser oxyrinchus)

Exposure to high temperatures can lead to thermotolerance in fish, which is hypothesized to potentially improve post-release survival in species under restocking programs, like Atlantic sturgeon. The aim of this study was to determine whether Atlantic sturgeon juveniles exposed to a 4-week temperature treatment respond differently to a subsequent heat shock than juveniles exposed to heat shock for the first time (naive fish). Response to heat shock was assessed by mapping the liver transcriptome. In total, 838 unique contigs were differentially expressed between the trained and the control group (592 downregulated, 261 upregulated, and 15 down- or upregulated, depending on the condition), corresponding to genes involved in the response to heat, tissue damage, proteolysis, and metabolism. Temperature-trained fish showed 2-4-fold fewer dysregulated contigs than naive fish, indicating their ability to maintain and recover homeostasis faster. During heat shock, hspc1 was upregulated in both experimental groups, while hspa1 and dnaja4 were exclusively upregulated in the control. Overall, compensatory mechanisms were observed in addition to the heat shock response. Only two genes, fgg and apnl, were upregulated at nearly all timepoints in both groups. Peptidases were more strongly downregulated in control fish, which also showed a reduction in lipid metabolism during recovery. Keratins, pck1, gadd45ga, and gadd45gb were differentially expressed between trained and control fish, and due to their roles in tissue protection and ER stress reduction, they might be responsible for the maintenance of the transcriptional homeostasis observed in trained fish.

Genetic Constraints, Transcriptome Plasticity, and the Evolutionary Response to Climate Change

In situ adaptation to climate change will be critical for the persistence of many ectotherm species due to their relative lack of dispersal capacity. Climate change is causing increases in both the mean and the variance of environmental temperature, each of which may act as agents of selection on different traits. Importantly, these traits may not be heritable or have the capacity to evolve independently from one another. When genetic constraints prevent the "baseline" values of thermal performance traits from evolving rapidly, phenotypic plasticity driven by gene expression might become critical. We review the literature for evidence that thermal performance traits in ectotherms are heritable and have genetic architectures that permit their unconstrained evolution. Next, we examine the relationship between gene expression and both the magnitude and duration of thermal stress. Finally, we identify genes that are likely to be important for adaptation to a changing climate and determine whether they show patterns consistent with thermal adaptation. Although few studies have measured narrow-sense heritabilities of thermal performance traits, current evidence suggests that the end points of thermal reaction norms (tolerance limits) are moderately heritable and have the potential to evolve rapidly. By contrast, performance at intermediate temperatures has substantially lower evolutionary potential. Moreover, evolution in many species appears to be constrained by genetic correlations such that populations can adapt to either increases in mean temperature or temperature variability, but not both. Finally, many species have the capacity for plastic expression of the transcriptome in response to temperature shifts, with the number of differentially expressed genes increasing with the magnitude, but not the duration, of thermal stress. We use these observations to develop a conceptual model that describes the likely trajectory of genome evolution in response to changes in environmental temperature. Our results indicate that extreme weather events, rather than gradual increases in mean temperature, are more likely to drive genetic and phenotypic change in wild ectotherms.

The synergistic interaction of thermal stress coupled with overstocking strongly modulates the transcriptomic activity and immune capacity of rainbow trout (Oncorhynchus mykiss)

The objective of the present study is to identify and evaluate informative indicators for the welfare of rainbow trout exposed to (A) a water temperature of 27 °C and (B) a stocking density of 100 kg/m³ combined with a temperature of 27 °C. The spleen-somatic and condition index, haematocrit and the concentrations of haemoglobin, plasma cortisol and glucose revealed non-significant differences between the two stress groups and the reference group 8 days after the onset of the experiments. The transcript abundance of almost 1,500 genes was modulated at least twofold in in the spleen of rainbow trout exposed to a critical temperature alone or a critical temperature combined with crowding as compared to the reference fish. The number of differentially expressed genes was four times higher in trout that were simultaneously challenged with high temperature and crowding, compared to trout challenged with high temperature alone. Based on these sets of differentially expressed genes, we identified unique and common tissue- and stress type-specific pathways. Furthermore, our subsequent immunologic analyses revealed reduced bactericidal and inflammatory activity and a significantly altered blood-cell composition in challenged versus non-challenged rainbow trout. Altogether, our data demonstrate that heat and overstocking exert synergistic effects on the rainbow trout’s physiology, especially on the immune system.

Multiplicity and Polymorphism of Fish Hemoglobins

The diversity of fish hemoglobins and the association with oxygen availability and physiological requirements during the life cycle has attracted scientists since the first report on multiple hemoglobin in fishes (Buhler and Shanks 1959). The functional heterogeneity of the fish hemoglobins enables many species to tolerate hypoxic conditions and exhausting swimming, but also to maintain the gas pressure in the swim bladder at large depths. The hemoglobin repertoire has further increased in various species displaying polymorphic hemoglobin variants differing in oxygen binding properties. The multiplicity of fish hemoglobins as particularly found in the tetraploid salmonids strongly contrasts with the complete loss of hemoglobins in Antarctic icefishes and illustrates the adaptive radiation in the oxygen transport of this successful vertebrate group.

Salmonid gene expression biomarkers indicative of physiological responses to changes in salinity and temperature, but not dissolved oxygen

An organism's ability to respond effectively to environmental change is critical to its survival. Yet, life stage and overall condition can dictate tolerance thresholds to heightened environmental stressors, such that stress may not be equally felt across individuals and at all times. Also, the transcriptional responses induced by environmental changes can reflect both generalized responses as well as others that are highly specific to the type of change being experienced. Thus, if transcriptional biomarkers specific to a stressor, even under multi-stressor conditions, can be identified, the biomarkers could then be applied in natural environments to determine when and where an individual experiences such a stressor. Here, we experimentally challenged juvenile Chinook salmon (Oncorhynchus tshawytscha) to validate candidate gill gene expression biomarkers. A sophisticated experimental design manipulated salinity (freshwater, brackish water and seawater), temperature (10, 14 and 18°C) and dissolved oxygen (normoxia and hypoxia) in all 18 possible combinations for 6 days using separate trials for three smolt statuses (pre-smolt, smolt and de-smolt). In addition, changes in juvenile behaviour, plasma variables, gill Na+/K+-ATPase activity, body size, body morphology and skin pigmentation supplemented the gene expression responses. We identified biomarkers specific to salinity and temperature that transcended the multiple stressors, smolt status and mortality (live, dead and moribund). Similar biomarkers for dissolved oxygen were not identified. This work demonstrates the unique power of gene expression biomarkers to identify a specific stressor even under multi-stressor conditions, and we discuss our next steps for hypoxia biomarkers using an RNA-seq study.

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.

Design and characterization of an 87k SNP genotyping array for Arctic charr (Salvelinus alpinus)

We have generated a high-density, high-throughput genotyping array for characterizing genome-wide variation in Arctic charr (Salvelinus alpinus). Novel single nucleotide polymorphisms (SNPs) were identified in charr from the Fraser, Nauyuk and Tree River aquaculture strains, which originated from northern Canada and fish from Iceland using high coverage sequencing, reduced representation sequencing and RNA-seq datasets. The array was designed to capture genome-wide variation from a diverse suite of Arctic charr populations. Cross validation of SNPs from various sources and comparison with previously published Arctic charr SNP data provided a set of candidate SNPs that generalize across populations. Further candidate SNPs were identified based on minor allele frequency, association with RNA transcripts, even spacing across intergenic regions and association with the sex determining (sdY) gene. The performance of the 86,503 SNP array was assessed by genotyping Fraser, Nauyuk and Tree River strain individuals, as well as wild Icelandic Arctic charr. Overall, 63,060 of the SNPs were polymorphic within at least one group and 36.8% were unique to one of the four groups, suggesting that the array design allows for characterization of both within and across population genetic diversity. The concordance between sdY markers and known phenotypic sex indicated that the array can accurately determine the sex of individuals based on genotype alone. The Salp87k genotyping array provides researchers and breeders the opportunity to analyze genetic variation in Arctic charr at a more detailed level than previously possible.

RNA-seq reveals temporal differences in the transcriptome response to acute heat stress in the Atlantic salmon (Salmo salar)

Acute heat stress is common in aquaculture and can affect diverse physiological processes in fish; however, different species of fish have various mechanisms for heat stress adaptation. In this study, we profiled the transcriptome responses of the Atlantic salmon (Salmo salar) to heat stress at 23 °C for 6 or 24 h, compared with that of fish at a normal temperature of 13 °C. The liver was selected as the target tissue for this analysis. A total of 243 and 88 genes were differentially expressed after 6 and 24 h of heat stress, respectively. Of these, only 22 were common to both time points, and most of these common genes were molecular chaperones such as heat shock cognate 71 kDa protein and heat shock protein 90-alpha. Genes such as activating transcription factor 6, calreticulin, protein disulfide isomerase A3, and protein kinase R-like endoplasmic reticulum kinase-eukaryotic initiation factor 2-alpha were only up-regulated after 6 h of heat stress; most of these genes are involved in the endoplasmic reticulum stress pathway. Indeed, endoplasmic reticulum stress was identified at 6 h but not at 24 h, suggesting that stress response plays an important role in the adaptation of Atlantic salmon to acute heat stress. Other up-regulated genes at 6 h were related to the insulin and nucleotide oligomerization domain-like receptor signaling pathways, which directly eliminate misfolded proteins and sustain sugar and lipid homeostasis. At 24 h, heat stress influenced the expression of steroid and terpenoid backbone biosynthesis, which may influence the sexual development and differentiation of Atlantic salmon. Overall, our results elucidate the transcriptome mechanisms that contribute to short-term heat tolerance in the liver of Atlantic salmon.

Transcriptomic responses to heat stress in rainbow trout Oncorhynchus mykiss head kidney

Rainbow trout (Oncorhynchus mykiss) are widely cultured throughout the word for commercial aquaculture. However, as a cold-water species, rainbow trout are highly susceptible to heat stress, which may cause pathological signs or diseases by alleviating the immune roles and then lead to mass mortality. Understanding the molecular mechanisms that occur in the rainbow trout in response to heat stress will be useful to decrease heat stress-related morbidity and mortality in trout aquaculture. In the present study, we conducted transcriptome analysis of head kidney tissue in rainbow trout under heat-stress (24 °C) and control (18 °C) conditions, to identify heat stress-induced genes and pathways. More than 281 million clean reads were generated from six head kidney libraries. Using an adjusted P-value of P < 0.05 as the threshold, a total of 443 differentially expressed genes (DEGs) were identified, including members of the HSP90, HSP70, HSP60, and HSP40 family and several cofactors or cochaperones. The RNA-seq results were confirmed by RT-qPCR. Gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis of DEGs were performed. Many genes involved in maintaining homeostasis or adapting to stress and stimuli were highly induced in response to high temperature. The most significantly enriched pathway was "Protein processing in endoplasmic reticulum (ER)", a quality control system that ensures correct protein folding or degradation of misfolded polypeptides by ER-associated degradation. Other signaling pathways involved in regulation of immune system and post-transcriptional regulation of spliceosome were also critical for thermal adaptation. These findings improve our understanding of the molecular mechanisms of heat stress responses and are useful to develop strategies for the improvement of rainbow trout survival rate during summer high-temperature period.

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.

Physiological and genomic signatures of evolutionary thermal adaptation in redband trout from extreme climates

Temperature is a master environmental factor that limits the geographical distribution of species, especially in ectotherms. To address challenges in biodiversity conservation under ongoing climate change, it is essential to characterize relevant functional limitations and adaptive genomic content at population and species levels. Here, we present evidence for adaptive divergence in cardiac function and genomic regions in redband trout (Oncorhynchus mykiss gairdneri) populations from desert and montane streams. Cardiac phenotypes of individual fish were measured in the field with a custom-built electrocardiogram apparatus. Maximum heart rate and its rate limiting temperature during acute warming were significantly higher in fish that have evolved in the extreme of a desert climate compared to a montane climate. Association mapping with 526,301 single nucleotide polymorphisms (SNPs) across the genome revealed signatures of thermal selection both within and among ecotypes. Among desert and montane populations, 435 SNPs were identified as putative outliers under natural selection and 20 of these loci showed significant association with average summer water temperatures among populations. Phenotypes for cardiac performance were variable within each ecotype, and 207 genomic regions were strongly associated with either maximum heart rate or rate limiting temperatures among individuals. Annotation of significant loci provided candidate genes that underlie thermal adaptation, including pathways associated with cardiac function (IRX5, CASQ1, CAC1D, and TITIN), neuroendocrine system (GPR17 and NOS), and stress response (SERPH). By integrating comparative physiology and population genomics, results here advance our knowledge on evolutionary processes of thermal adaptation in aquatic ectotherms.

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.

Gradual and Acute Temperature Rise Induces Crossing Endocrine, Metabolic, and Immunological Pathways in Maraena Whitefish (Coregonus maraena)

The complex and still poorly understood nature of thermoregulation in various fish species complicates the determination of the physiological status on the basis of diagnostic marker genes and indicative molecular pathways. The present study aimed to compare the physiological impacts of both gradual and acute temperature rise from 18 to 24°C on maraena whitefish in aquaculture. Microarray-based transcriptome profiles in the liver, spleen and kidney of heat-stressed maraena whitefish revealed the modulation of a significantly higher number of genes in those groups exposed to gradually rising temperatures compared with the acutely stressed groups, which might reflect early adaptation mechanisms. Moreover, we suggest a common set of 11 differentially expressed genes that indicate thermal stress induced by gradual or acute temperature rise in the three selected tissues. Besides the two pathways regulated in both data sets unfolded protein response and aldosterone signaling in epithelial cells, we identified unique tissue- and stress type-specific pathways reflecting the crossroads between signal transduction, metabolic and immunologic pathways to cope with thermal stress. In addition, comparing lists of differentially regulated genes with meta-analyzed published data sets revealed that “acute temperature rise”-responding genes that encode members of the HSP70, HSP90, and HSP40 families; their functional homologs; co-chaperones and stress-signal transducers are well-conserved across different species, tissues and/or cell types and experimental approaches.

Mechanisms of thermal adaptation and evolutionary potential of conspecific populations to changing environments

Heterogeneous and ever-changing thermal environments drive the evolution of populations and species, especially when extreme conditions increase selection pressure for traits influencing fitness. However, projections of biological diversity under scenarios of climate change rarely consider evolutionary adaptive potential of natural species. In this study, we tested for mechanistic evidence of evolutionary thermal adaptation among ecologically divergent redband trout populations (Oncorhynchus mykiss gairdneri) in cardiorespiratory function, cellular response and genomic variation. In a common garden environment, fish from an extreme desert climate had significantly higher critical thermal maximum (p < .05) and broader optimum thermal window for aerobic scope (>3°C) than fish from cooler montane climate. In addition, the desert population had the highest maximum heart rate during warming (20% greater than montane populations), indicating improved capacity to deliver oxygen to internal tissues. In response to acute heat stress, distinct sets of cardiac genes were induced among ecotypes, which helps to explain the differences in cardiorespiratory function. Candidate genomic markers and genes underlying these physiological adaptations were also pinpointed, such as genes involved in stress response and metabolic activity (hsp40, ldh-b and camkk2). These markers were developed into a multivariate model that not only accurately predicted critical thermal maxima, but also evolutionary limit of thermal adaptation in these specific redband trout populations relative to the expected limit for the species. This study demonstrates mechanisms and limitations of an aquatic species to evolve under changing environments that can be incorporated into advanced models to predict ecological consequences of climate change for natural organisms.

Differential Gene Expression Profiles and Alternative Isoform Regulations in Gill of Nile Tilapia in Response to Acute Hypoxia

Fish often encounters exposures to acute environmental hypoxia either spatially or temporally. Gill organ plays important roles in response to hypoxic stress in fish. Few studies focus on the molecular regulation mechanisms of gills under hypoxic stress. In this study, we investigated the transcriptomic response to 12-h acute hypoxia in gill of a hypoxia tolerant fish, Nile tilapia Oreochromis niloticus through RNA sequencing (RNA-Seq). We sequenced messenger RNA from three control samples and three hypoxia-treated samples. Bioinformatics analysis identified 239 differentially expressed genes (DEG) and 34 genes (DUES) that had significant differential alternative isoform regulation events in at least one exonic region in gill in response to acute hypoxia. The spatiotemporal expression analysis in five tissues (heart, liver, brain, gill, and spleen) sampled at three time points (6, 12, and 24 h) under hypoxia treatment confirmed the significant association of differential exon usages in two DUES genes (TLDC2 and SSX2IPA) with hypoxia conditions. Further functional analysis suggested several energy and immune response-related pathways, e.g., metabolic pathway and antigen processing and presentation, contained the most abundant DEG genes. We found that some GO biological processes for DEG genes were significantly enriched under hypoxic stress, such as glycolysis, metabolic process, generation of precursor metabolites and energy, and cholesterol metabolic process. Our findings suggest abundant differential gene expression changes and alternative isoform regulation events in genes involved in the hypoxia response in gill. Our results provide a basis for exploring the gene regulation mechanism under hypoxic stress in fish.

BREEDING AND GENETICS SYMPOSIUM: Breeding for resilience to heat stress effects in dairy ruminants. A comprehensive review

Selection for heat tolerant (HT) animals in dairy production has been so far linked to estimation of declines in production using milk recording and meteorological information on the day of control using reaction norms. Results from these models show that there is a reasonable amount of genetic variability in the individual response to high heat loads, which makes feasible selection of HT animals at low costs. However, the antagonistic relationship between level of production and response to heat stress (HS) implies that selection for HT animals under this approach must be done with caution so that productivity is not damaged. Decomposition of the genetic variability in principal components (PC) can provide selection criteria independent of milk production level although biological interpretation of PC is difficult. Moreover, given that response to heat stress for each animal is estimated with very sparse information collected under different physiological and management circumstances, biased (normally underestimation) and lack of accuracy may be expected. Alternative phenotypic characterization of HT can come from the use of physiological traits, which have also shown moderate heritability. However, costs of a large scale implementation based on physiological characteristics has precluded its use. Another alternative is the use of biomarkers that define heat tolerance. A review of biomarkers of HS from more recent studies is provided. Of particular interest are milk biomarkers, which together with infrared spectra prediction equations can provide useful tools for genetic selection. In the 'omics' era, genomics, transcriptomics, proteomics and metabolomics have been already used to detect genes affecting HT. A review of findings in these areas is also provided. Except for the slick hair gene, there are no other genes for which variants have been clearly associated with HT. However, integration of omics information could help in pointing at knots of the HS control network and, in the end, to a panel of markers to be used in the selection of HT animals. Overall, HT is a complex phenomenon that requires integration of fine phenotypes and omics information to provide accurate tools for selection without damaging productivity. Technological developments to make on-farm implementation feasible and with greater insight into the key biomarkers and genes involved in HT are needed.

Genomic Changes Associated with Reproductive and Migratory Ecotypes in Sockeye Salmon (Oncorhynchus nerka)

Mechanisms underlying adaptive evolution can best be explored using paired populations displaying similar phenotypic divergence, illuminating the genomic changes associated with specific life history traits. Here, we used paired migratory [anadromous vs. resident (kokanee)] and reproductive [shore- vs. stream-spawning] ecotypes of sockeye salmon (Oncorhynchus nerka) sampled from seven lakes and two rivers spanning three catchments (Columbia, Fraser, and Skeena) in British Columbia, Canada to investigate the patterns and processes underlying their divergence. Restriction-site associated DNA sequencing was used to genotype this sampling at 7,347 single nucleotide polymorphisms, 334 of which were identified as outlier loci and candidates for divergent selection within at least one ecotype comparison. Sixty-eight of these outliers were present in two or more comparisons, with 33 detected across multiple catchments. Of particular note, one locus was detected as the most significant outlier between shore and stream-spawning ecotypes in multiple comparisons and across catchments (Columbia, Fraser, and Snake). We also detected several genomic islands of divergence, some shared among comparisons, potentially showing linked signals of differential selection. The single nucleotide polymorphisms and genomic regions identified in our study offer a range of mechanistic hypotheses associated with the genetic basis of O. nerka life history variation and provide novel tools for informing fisheries management.

Genomic signatures of local adaptation reveal source-sink dynamics in a high gene flow fish species

Understanding source-sink dynamics is important for conservation management, particularly when climatic events alter species’ distributions. Following a 2011 ‘marine heatwave’ in Western Australia, we observed high recruitment of the endemic fisheries target species Choerodon rubescens, towards the cooler (southern) end of its distribution. Here, we use a genome wide set of 14 559 single-nucleotide polymorphisms (SNPs) to identify the likely source population for this recruitment event. Most loci (76%) showed low genetic divergence across the species’ range, indicating high levels of gene flow and confirming previous findings using neutral microsatellite markers. However, a small proportion of loci showed strong patterns of differentiation and exhibited patterns of population structure consistent with local adaptation. Clustering analyses based on these outlier loci indicated that recruits at the southern end of C. rubescens’ range originated 400 km to the north, at the centre of the species’ range, where average temperatures are up to 3 °C warmer. Survival of these recruits may be low because they carry alleles adapted to an environment different to the one they now reside in, but their survival is key to establishing locally adapted populations at and beyond the range edge as water temperatures increase with climate change.

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.

The expression and function of hsp30-like small heat shock protein genes in amphibians, birds, fish, and reptiles

Small heat shock proteins (sHSPs) are a superfamily of molecular chaperones with important roles in protein homeostasis and other cellular functions. Amphibians, reptiles, fish and birds have a shsp gene called hsp30, which was also referred to as hspb11 or hsp25 in some fish and bird species. Hsp30 genes, which are not found in mammals, are transcribed in response to heat shock or other stresses by means of the heat shock factor that is activated in response to an accumulation of unfolded protein. Amino acid sequence analysis revealed that representative HSP30s from different classes of non-mammalian vertebrates were distinct from other sHSPs including HSPB1/HSP27. Studies with amphibian and fish recombinant HSP30 determined that they were molecular chaperones since they inhibited heat- or chemically-induced aggregation of unfolded protein. During non-mammalian vertebrate development, hsp30 genes were differentially expressed in selected tissues. Also, heat shock-induced stage-specific expression of hsp30 genes in frog embryos was regulated at the level of chromatin structure. In adults and/or tissue culture cells, hsp30 gene expression was induced by heat shock, arsenite, cadmium or proteasomal inhibitors, all of which enhanced the production of unfolded/damaged protein. Finally, immunocytochemical analysis of frog and chicken tissue culture cells revealed that proteotoxic stress-induced HSP30 accumulation co-localized with aggresome-like inclusion bodies. The congregation of damaged protein in aggresomes minimizes the toxic effect of aggregated protein dispersed throughout the cell. The current availability of probes to detect the presence of hsp30 mRNA or encoded protein has resulted in the increased use of hsp30 gene expression as a marker of proteotoxic stress in non-mammalian vertebrates.

Characterisation of the global transcriptional response to heat shock and the impact of individual genetic variation

Background

The heat shock transcriptional response is essential to effective cellular function under stress. This is a highly heritable trait but the nature and extent of inter-individual variation in heat shock response remains unresolved.

Methods

We determined global transcription profiles of the heat shock response for a panel of lymphoblastoid cell lines established from 60 founder individuals in the Yoruba HapMap population. We explore the observed differentially expressed gene sets following heat shock, establishing functional annotations, underlying networks and nodal genes involving heat shock factor 1 recruitment. We define a multivariate phenotype for the global transcriptional response to heat shock using partial least squares regression and map this quantitative trait to associated genetic variation in search of the major genomic modulators.

Results

A comprehensive dataset of differentially expressed genes following heat shock in humans is presented. We identify nodal genes downstream of heat shock factor 1 in this gene set, notably involving ubiquitin C and small ubiquitin-like modifiers together with transcription factors. We dissect a multivariate phenotype for the global heat shock response which reveals distinct clustering of individuals in terms of variance of the heat shock response and involves differential expression of genes involved in DNA replication and cell division in some individuals. We find evidence of genetic associations for this multivariate response phenotype that involves trans effects modulating expression of genes following heat shock, including HSF1 and UBQLN1.

Conclusion

This study defines gene expression following heat shock for a cohort of individuals, establishing insights into the biology of the heat shock response and hypotheses for how variation in this may be modulated by underlying genetic diversity.

Electronic supplementary material

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

A genome-wide association study of heat stress-associated SNPs in catfish

Heat tolerance is a complex and economically important trait for catfish genetic breeding programs. With global climate change, it is becoming an increasingly important trait. To better understand the molecular basis of heat stress, a genome-wide association study (GWAS) was carried out using the 250 K catfish SNP array with interspecific backcross progenies, which derived from crossing female channel catfish with male F1 hybrid catfish (female channel catfish × male blue catfish). Three significant associated SNPs were detected by performing an EMMAX approach for GWAS. The SNP located on linkage group 14 explained 12.1% of phenotypical variation. The other two SNPs, located on linkage group 16, explained 11.3 and 11.5% of phenotypical variation respectively. A total of 14 genes with heat stress related functions were detected within the significant associated regions. Among them, five genes-TRAF2, FBXW5, ANAPC2, UBR1 and KLHL29- have known functions in the protein degradation process through the ubiquitination pathway. Other genes related to heat stress include genes involved in protein biosynthesis (PRPF4 and SYNCRIP), protein folding (DNAJC25), molecule and iron transport (SLC25A46 and CLIC5), cytoskeletal reorganization (COL12A1) and energy metabolism (COX7A2, PLCB1 and PLCB4) processes. The results provide fundamental information about genes and pathways that is useful for further investigation into the molecular mechanisms of heat stress. The associated SNPs could be promising candidates for selecting heat-tolerant catfish lines after validating their effects on larger and various catfish populations.

The pathophysiological basis and consequences of fever

There are numerous causes of a raised core temperature. A fever occurring in sepsis may be associated with a survival benefit. However, this is not the case for non-infective triggers. Where heat generation exceeds heat loss and the core temperature rises above that set by the hypothalamus, a combination of cellular, local, organ-specific, and systemic effects occurs and puts the individual at risk of both short-term and long-term dysfunction which, if severe or sustained, may lead to death. This narrative review is part of a series that will outline the pathophysiology of pyrogenic and non-pyrogenic fever, concentrating primarily on the pathophysiology of non-septic causes.

Transcriptomic responses to environmental temperature in eurythermal and stenothermal fishes

Ectothermic species like fishes differ greatly in the thermal ranges they tolerate; some eurythermal species may encounter temperature ranges in excess of 25°C, whereas stenothermal species in polar and tropical waters live at essentially constant temperatures. Thermal specialization comes with fitness trade-offs and as temperature increases due to global warming, the physiological basis of specialization and thermal plasticity has become of great interest. Over the past 50 years, comparative physiologists have studied the physiological and molecular differences between stenothermal and eurythermal fishes. It is now well known that many stenothermal fishes have lost an inducible heat shock response (HSR). Recent advances in transcriptomics have now made it possible to examine genome-wide changes in gene expression (GE) in non-model ecologically important fish, broadening our view beyond the HSR to regulation of genes involved in hundreds of other cellular processes. Here, we review the major findings from transcriptomic studies of extreme eurythermal and stenothermal fishes in response to acute and long-term exposure to temperature, both time scales being critically important for predicting climate change responses. We consider possible molecular adaptations that underlie eurythermy and stenothermy in teleosts. Furthermore, we highlight the challenges that still face the field of comparative environmental genomics and suggest fruitful paths of future investigation.

Differentially expressed genes associated with adaptation to different thermal environments in three sympatric Cuban Anolis lizards

How animals achieve evolutionary adaptation to different thermal environments is an important issue for evolutionary biology as well as for biodiversity conservation in the context of recent global warming. In Cuba, three sympatric species of Anolis lizards (Anolis allogus, A. homolechis and A. sagrei) inhabit different thermal microhabitats, thereby providing an excellent opportunity to examine how they have adapted to different environmental temperatures. Here, we performed RNA-seq on the brain, liver and skin tissues from these three species to analyse their transcriptional responses at two different temperatures. In total, we identified 400, 816 and 781 differentially expressed genes (DEGs) between the two temperatures in A. allogus, A. homolechis and A. sagrei, respectively. Only 62 of these DEGs were shared across the three species, indicating that global transcriptional responses have diverged among these species. Gene ontology (GO) analysis showed that large numbers of ribosomal protein genes were DEGs in the warm-adapted A. homolechis, suggesting that the upregulation of protein synthesis is an important physiological mechanism in the adaptation of this species to hotter environments. GO analysis also showed that GO terms associated with circadian regulation were enriched in all three species. A gene associated with circadian regulation, Nr1d1, was detected as a DEG with opposite expression patterns between the cool-adapted A. allogus and the hot-adapted A. sagrei. Because the environmental temperature fluctuates more widely in open habitats than in forests throughout the day, the circadian thermoregulation could also be important for adaptation to distinct thermal habitats.

Transcriptomic Characterization of Tambaqui (Colossoma macropomum, Cuvier, 1818) Exposed to Three Climate Change Scenarios

Climate change substantially affects biodiversity around the world, especially in the Amazon region, which is home to a significant portion of the world’s biodiversity. Freshwater fishes are susceptible to increases in water temperature and variations in the concentrations of dissolved gases, especially oxygen and carbon dioxide. It is important to understand the mechanisms underlying the physiological and biochemical abilities of fishes to survive such environmental changes. In the present study, we applied RNA-Seq and de novo transcriptome sequencing to evaluate transcriptome alterations in tambaqui when exposed to five or fifteen days of the B1, A1B and A2 climate scenarios foreseen by the IPCC. The generated ESTs were assembled into 54,206 contigs. Gene ontology analysis and the STRING tool were then used to identify candidate protein domains, genes and gene families potentially responsible for the adaptation of tambaqui to climate changes. After sequencing eight RNA-Seq libraries, 32,512 genes were identified and mapped using the Danio rerio genome as a reference. In total, 236 and 209 genes were differentially expressed at five and fifteen days, respectively, including chaperones, energetic metabolism-related genes, translation initiation factors and ribosomal genes. Gene ontology enrichment analysis revealed that mitochondrion, protein binding, protein metabolic process, metabolic processes, gene expression, structural constituent of ribosome and translation were the most represented terms. In addition, 1,202 simple sequence repeats were detected, 88 of which qualified for primer design. These results show that cellular response to climate change in tambaqui is complex, involving many genes, and it may be controlled by different cues and transcription/translation regulation mechanisms. The data generated from this study provide a valuable resource for further studies on the molecular mechanisms involved in the adaptation of tambaqui and other closely related teleost species to climate change.

Transcriptional response to heat shock in liver of snow trout (Schizothorax richardsonii)--a vulnerable Himalayan Cyprinid fish

The snow trout (Schizothorax richardsonii) belonging to family Cyprinidae, is an endemic fish of the Himalayan region. The species is tagged as vulnerable species in the IUCN red list of threatened species. The fish thrives well in snowmelt water of several streams and rivers in the region but are occasionally exposed to more than 20 °C during the summer season. Therefore, we have used deep RNA sequencing to decipher the transcriptome of snow trout and characterize the genes and molecular pathways involved in heat shock response. In this study 72,601,298 and 65,428,283 raw reads for heat-shocked and control, respectively, were obtained by Illumina paired-end sequencing technology. The de novo assembled transcriptome was tested for differential gene expression across the treatment groups. The quality of assembly was high with N75 and N50 lengths of 461 and 1274 bases, respectively. A total of 65 unique transcripts were differentially expressed in liver under heat shock and control. Annotated blast matches reveal that differentially expressed transcripts correspond to critical chaperones and molecular pathways, previously shown to be important for thermal stress in other fish species. Eight randomly selected heat-stressed responsive transcripts were also observed to be upregulated during qRT-PCR analysis. This study is the preliminary step to understanding the responses during sudden environmental changes like heat shock. The reference transcriptome database would also aid further studies on biological and physiological aspects of the snow trout under abiotic stresses.

Developmental abnormalities and differential expression of genes induced in oil and dispersant exposed Menidia beryllina embryos

Exposure of fish embryos to relatively low concentrations of oil has been implicated in sub-lethal toxicity. The objective of this study was to determine the effects of the exposure of Menidia beryllina embryos at 30-48h post-fertilization to the water accommodated fractions of oil (WAF, 200ppm, v/v), dispersants (20ppm, v/v, Corexit 9500 or 9527), and mixtures of oil and each of the dispersants to produce chemically enhanced water accommodated fractions (CEWAFs) over a 72-hour period. The polyaromatic hydrocarbon (PAH) and benzene, toluene, ethylene and xylene (BTEX) constituents of the 5X concentrated exposure solutions (control, WAF, dispersants and CEWAFs) were determined and those of the 1× exposures were derived using a dilution factor. PAH, BTEX and low molecular weight PAH constituents greater than 1ppb were observed in WAF and the dispersants, but at much higher levels in CEWAFs. The WAF and CEWAFs post-weathering were diluted at 1:5 (200ml WAF/CEWAF: 800ml 25ppt saltwater) for embryo exposures. Mortality, heartbeat, embryo normalcy, abnormality types and severities were recorded. The qPCR assay was used to quantify abundances of transcripts of target genes for sexual differentiation and sex determination (StAR, dmrt-1, amh, cyp19b, vtg and chg-L,), growth regulation (ghr) and stress response (cyp1a and Hsp90); and gapdh served as the housekeeping gene. Temperature was 21±1.5°C throughout the experimental period, while mortality was low and not significantly different (p=0.68) among treatments. Heartbeat was significantly different (0.0034) with the lowest heartbeats recorded in Corexit 9500 (67.5beats/min) and 9527 (67.1beats/min) exposed embryos compared with controls (82.7beats/min). Significantly more treated embryos were in a state of deterioration, with significantly more embryos presenting arrested tissue differentiation compared with controls (p=0.021). Exposure to WAF, dispersants and CEWAF induced aberrant expression of all the genes, with star, dmrt-1, ghr and hsp90 being significantly down-regulated in CEWAF and cyp19b in Corexit 9527. The cyp1a and cyp19b were significantly up-regulated in CEWAFs and WAF, respectively. The molecular endpoints were most sensitive, especially the expression of star, cyp19b, cyp1a, hsp90 and could therefore be used as early indicators of long term effects of Corexit 9500 and 9527 usage in oil spill management on M. beryllina, a valid sentinel for oil pollution events.

Impact of Thermal Stress on Kidney-Specific Gene Expression in Farmed Regional and Imported Rainbow Trout

Seasonal water temperatures can be stressful for fish in aquaculture and can therefore negatively influence their welfare. Although the kidney is the crucial organ associated with the primary stress response, knowledge about the stress-modulated kidney transcriptome in salmonids is limited. In the present study, we used a comparative microarray approach to characterize the general gene expression profiles of rainbow trout trunk kidney after a 2-week acclimation to mild heat (23 °C) and cold stress (8 °C). Hypothesizing that local adaptation influences stress performance, we aimed to identify differences in the temperature-induced gene expression in the regional trout strain BORN, in addition to a common imported strain. Moderate temperature challenge provoked typical stress response clusters, including heat-shock proteins or cold-inducible factors, in addition to altered energy metabolism in trout kidney. Mild cold, in particular, enhanced renal protein degradation processes, as well as mRNA and protein synthesis, while it also triggered fatty acid biosynthesis. Mild heat led to cytoskeleton-stabilizing processes and might have facilitated cell damage and infection. Furthermore, both breeding lines used different strategies for energy provision, cellular defense, and cell death/survival pathways. As a main finding, the genes involved in energy provision showed generally higher transcript levels at both temperatures in BORN trout compared to imported trout, indicating adjusted metabolic rates under local environmental conditions. Altogether, this study provides a general overview of stress-induced transcriptional patterns in rainbow trout trunk kidney, in addition to identifying genes and networks that contribute to the robustness of the BORN strain. Our analyses suggest SERPINH1 and CIRBP as general marker genes for heat stress and cold stress in trout, respectively.

Transcriptional Response to Acute Thermal Exposure in Juvenile Chinook Salmon Determined by RNAseq

Thermal exposure is a serious and growing challenge facing fish species worldwide. Chinook salmon (Oncorhynchus tshawytscha) living in the southern portion of their native range are particularly likely to encounter warmer water due to a confluence of factors. River alterations have increased the likelihood that juveniles will be exposed to warm water temperatures during their freshwater life stage, which can negatively impact survival, growth, and development and pose a threat to dwindling salmon populations. To better understand how acute thermal exposure affects the biology of salmon, we performed a transcriptional analysis of gill tissue from Chinook salmon juveniles reared at 12° and exposed acutely to water temperatures ranging from ideal to potentially lethal (12° to 25°). Reverse-transcribed RNA libraries were sequenced on the Illumina HiSeq2000 platform and a de novo reference transcriptome was created. Differentially expressed transcripts were annotated using Blast2GO and relevant gene clusters were identified. In addition to a high degree of downregulation of a wide range of genes, we found upregulation of genes involved in protein folding/rescue, protein degradation, cell death, oxidative stress, metabolism, inflammation/immunity, transcription/translation, ion transport, cell cycle/growth, cell signaling, cellular trafficking, and structure/cytoskeleton. These results demonstrate the complex multi-modal cellular response to thermal stress in juvenile salmon.

Characterization of differentially expressed genes in liver in response to the rearing temperature of rainbow trout Oncorhynchus mykiss and their heritable differences

To characterize thermal-responsive genes in fish, firstly, juvenile rainbow trout were reared in four different temperature conditions (average temperatures were 10, 14, 18, and 22 °C, respectively) and differentially expressed genes were identified. Gene expression in the liver was analyzed by the differential display method, followed by validation using real-time PCR. Subsequently, to examine whether the identified genes show heritable differences, the gene expression levels were compared among juveniles of three genetically distinct lines of rainbow trout (a strain and two closed colonies) by rearing at two different temperature conditions (average 14 and 22 °C). By rearing at 22 °C, growth retardation was observed compared with fish reared at 14 and 18 °C, and six genes were identified as differentially expressed genes in response to the rearing temperature in the gene expression analyses. With the increase in rearing temperature, gene expressions of a complement C1q and two ribosomal proteins were significantly up-regulated. On the other hand, three metabolic genes (betaine homocysteine methyltransferase, triosephosphate isomerase, and glucose-6-phosphatase) were down-regulated, indicating a metabolic depression due to high temperature. In the subsequent analyses, in response to the rearing temperature (14 and 22 °C), there was a trend that the complement C1q and glucose-6-phosphatase genes showed different expression patterns among the three rainbow trout lines, suggesting heritable differences in these genes. Our study provides information on thermal-responsive genes in fish, and we anticipate it will facilitate further investigation in the thermal biology of fish.

Facing warm temperatures during migration: cardiac mRNA responses of two adult Oncorhynchus nerka populations to warming and swimming challenges

The main findings of the current study were that exposing adult sockeye salmon Onchorhynchus nerka to a warm temperature that they regularly encounter during their river migration induced a heat shock response at an mRNA level, and this response was exacerbated with forced swimming. Similar to the heat shock response, increased immune defence-related responses were also observed after warm temperature treatment and with a swimming challenge in two different populations (Chilko and Nechako), but with some important differences. Microarray analyses revealed that 347 genes were differentially expressed between the cold (12-13° C) and warm (18-19° C) treated fish, with stress response (GO:0006950) and response to fungus (GO:0009620) elevated with warm treatment, while expression for genes involved in oxidative phosphorylation (GO:0006119) and electron transport chain (GO:0022900) elevated for cold-treated fish. Analysis of single genes with real-time quantitative PCR revealed that temperature had the most significant effect on mRNA expression levels, with swimming and population having secondary influences. Warm temperature treatment for the Chilko population induced expression of heat shock protein (hsp) 90α, hsp90β and hsp30 as well as interferon-inducible protein. The Nechako population, which is known to have a narrower thermal tolerance window than the Chilko population, showed even more pronounced stress responses to the warm treatment and there was significant interaction between population and temperature treatment for hsp90β expression. Moreover, significant interactions were noted between temperature treatment and swimming challenge for hsp90α and hsp30, and while swimming challenge alone increased expression of these hsps, the expression levels were significantly elevated in warm-treated fish swum to exhaustion. In conclusion, it seems that adult O. nerka currently encounter conditions that induce several cellular defence mechanisms during their once-in-the-lifetime migration. As river temperatures continue to increase, it remains to be seen whether or not these cellular defences provide sufficient protection for all O. nerka populations.

Detection and mapping of QTL for temperature tolerance and body size in Chinook salmon (Oncorhynchus tshawytscha) using genotyping by sequencing

Understanding how organisms interact with their environments is increasingly important for conservation efforts in many species, especially in light of highly anticipated climate changes. One method for understanding this relationship is to use genetic maps and QTL mapping to detect genomic regions linked to phenotypic traits of importance for adaptation. We used high-throughput genotyping by sequencing (GBS) to both detect and map thousands of SNPs in haploid Chinook salmon (Oncorhynchus tshawytscha). We next applied this map to detect QTL related to temperature tolerance and body size in families of diploid Chinook salmon. Using these techniques, we mapped 3534 SNPs in 34 linkage groups which is consistent with the haploid chromosome number for Chinook salmon. We successfully detected three QTL for temperature tolerance and one QTL for body size at the experiment-wide level, as well as additional QTL significant at the chromosome-wide level. The use of haploids coupled with GBS provides a robust pathway to rapidly develop genomic resources in nonmodel organisms; these QTL represent preliminary progress toward linking traits of conservation interest to regions in the Chinook salmon genome.