Research articleHeritable variation in heat shock gene expression: a potential mechanism for adaptation to thermal stress in embryos of sea turtles
Abstract
The capacity of species to respond adaptively to warming temperatures will be key to their survival in the Anthropocene. The embryos of egg-laying species such as sea turtles have limited behavioural means for avoiding high nest temperatures, and responses at the physiological level may be critical to coping with predicted global temperature increases. Using the loggerhead sea turtle (Caretta caretta) as a model, we used quantitative PCR to characterise variation in the expression response of heat-shock genes (hsp60, hsp70 and hsp90; molecular chaperones involved in cellular stress response) to an acute non-lethal heat shock. We show significant variation in gene expression at the clutch and population levels for some, but not all hsp genes. Using pedigree information, we estimated heritabilities of the expression response of hsp genes to heat shock and demonstrated both maternal and additive genetic effects. This is the first evidence that the heat-shock response is heritable in sea turtles and operates at the embryonic stage in any reptile. The presence of heritable variation in the expression of key thermotolerance genes is necessary for sea turtles to adapt at a molecular level to warming incubation environments.
1. Introduction
Survival in a warming world depends on the ability of a species to move or remain within their current distributions and respond to changing environmental conditions by evolutionary and/or plastic responses [1–3]. Responses to changing environmental conditions include modification of life-history traits such as the timing of reproduction [3], changes in thermoregulatory behaviour and microhabitat use [4], and/or changes at a physiological level by acclimation or evolving tolerance to warmer climates [2,5]. Comparatively little attention has been paid to the capacity of embryos to survive thermal challenges, despite the fact that most embryos develop in fixed locations where stressful conditions are difficult to avoid [6]. Embryonic tolerance of thermal stress should be under strong selection, and plasticity of their physiological responses should further maximize their survival in altered environments (e.g. [7]).
Physiological responses to thermal stress begin at the molecular level, frequently with the activation of genes such as those that code for heat-shock proteins (Hsps) that mitigate damage to membranes, proteins and DNA [8–10]. Hsps in the Hsp60, Hsp70 and Hsp90 families increase expression in response to heat stress in a range of taxa (e.g. insects, fish, frogs and reptiles), and maintain protein folding and degradation and myelination of neurons to avoid apoptosis [11–17]. Because the upregulation of heat-shock protein genes (hsp) in response to thermal stress is highly conserved across taxa, it is a candidate mechanism for adaptation of embryos to higher temperatures. Consistent with this view, many studies have shown improved heat resistance after Hsp expression [8,10], and several studies have reported clinal variation in hsp genes across latitudinal [18,19] and altitudinal transects [20–22], indicating natural selection is acting on these genes. However, many factors influence Hsp expression levels (e.g. other environmental stressors, inbreeding, age) and the Hsp expression level in each species will be a balance between the costs and benefits of upregulation of Hsps [10].
The mechanisms by which sea turtle embryos respond to heat stress are of particular interest, as their lineage arose in the Mesozoic and has persisted through many instances of global heating and cooling [23–25]. The strong fidelity of nesting females to natal beaches [26] means that major shifts in rookery locations occurs slowly [27]. Consequently, rapid rises in beach temperatures due to anthropogenically forced climate change could be a novel challenge for embryos relative to more gradual thermal changes experienced throughout their evolutionary history.
Microevolution of the thermally sensitive genome of sea turtles is only possible if two key expectations are met. First, and taking heat shock genes as an example, the optimal level of hsp expression should vary among populations experiencing thermal differences in the nesting environment [28]. Second, the variation in expression of hsps in response to thermal stress must be heritable [28,29]. In general, there have been few studies on the adaptive potential of wild populations to heat stress, and laboratory evolution experiments have focused on the heritability of expression of hsp70 in most Drosophila species [30,31]. Given that Drosophila have short generation intervals (a few weeks), low levels of heritability may be sufficient for adaptation to rapid warming through evolutionary means. By contrast, sea turtles have long generation intervals of approximately 20–30 years [32], which constrains the rate of microevolution. Furthermore, if changes in female nesting behaviour cannot compensate for warmer incubation temperatures, sea turtle embryos are left with few options to avoid heat stress. The extent to which expression of hsps is heritable, or can be modified through plastic responses, will affect their ability to change their critical thermal limits.
Here, we investigated phenotypic and genetic variation of hsp gene expression in the embryos of loggerhead sea turtles (Caretta caretta), which are vulnerable to reduced fitness and higher mortality when exposed to high temperatures in terrestrial nests [33,34]. We first determined whether hsp expression differed between embryos from a temperate and a sub-tropical rookery. Although hsp sequences are generally highly conserved within and between species [11–17], we expected to find regional variation in hsp expression profiles. Secondly, we tested for and estimated the heritability of hsp expression, and the plasticity of hsp expression, in response to thermal stress. Our results provide some key parameters needed for understanding whether long-lived reptiles could adapt to the unprecedented pace of contemporary climate change through plastic and/or evolutionary adaptation.
2. Material and methods
Full methods and associated references are available in the electronic supplementary material, Materials and methods.
(a) Study species and sample collection
The Western Australia population of Caretta caretta nest at Turtle Bay on Dirk Hartog Island (DHI; 25.49827° S, 112.98719° E) and at smaller mainland rookeries including Bungelup Beach (BB; 22.282331° S, 113.831570° E) (figure 1). We collected a total of 1280 eggs from both nesting sites during peak nesting activity in January of 2013 and 2014. Eggs were collected, transported and monitored according to University of Western Australia ethics protocols. Details of these protocols can be found in the electronic supplementary material, Material and methods.
Figure 1. Dirk Hartog Island (DHI) and Bungelup Beach (BB) are two range-edge C. caretta rookeries in Western Australia, isolated by approximately 520 km of coastline [35] and 3° latitude. DHI is located within a temperate zone while BB is a located within a sub-tropical zone [36]. Map adapted from [37].
(b) Incubation and heat-shock experiments
Eggs from each female were randomly distributed among 1.5 l plastic containers to reduce clutch effects (DHI: 80 containers each with 15 eggs, BB: 20 containers each with four eggs). All eggs were held at a constant 29°C (±0.3°C) until time of heat-shock treatment. A heat shock of 36°C for 3 h was applied to the treatment group 45 days into the incubation period, following previously established protocols for inducing thermal stress in C. caretta [17]. After a 1-h cool down period, embryos were removed from the egg, weighed and given a lethal injection of MS-222 (50 mg kg−1, Sigma) [17,38] followed by decapitation to ensure death.
(c) RNA extraction and RT-qPCR
Two independent RNA extractions of heart tissue from up to five embryos per treatment group were done using the FavorPrep™ Total Tissue Mini RNA Kit (FATRK-300, Fisher Biotec) without DNase treatments, following the manufacturer's instructions. Total RNA was eluted in nuclease-free water and quantified using a Qubit v. 2.0 Fluorometer (Invitrogen). Two hundred nanograms of total RNA was reverse-transcribed into cDNA using the oligo (dT)16 primer with MultiScribe™ MuLV Reverse Transcriptase in a High Capacity RNA-to-cDNA Kit (4387406, Invitrogen), following the manufacturer's protocol. Quantitative real-time PCR was performed in triplicate using the iTaq Universal SYBR® Green Supermix (172–5121, BioRad) on a Step-One Plus PCR System (Applied Biosystems) with the following program: 95°C for 10 min; 40 cycles of 95°C for 10 s, 60°C for 60 s. RT-qPCR was performed in a 10 µl reaction with 10 ng of cDNA and final primer concentration of 200 nM. The gene-specific primers used for RT-qPCR are listed in the electronic supplementary material, table S1. Target gene mRNA expression levels were normalized to reference gene 18s mRNA expression levels [17].
(d) Statistical analyses
Cycle threshold (Cq) values were converted to relative gene expression values ΔCq and ΔΔCq using the methods described in similar studies [17,39]. Details of our calculations are provided in the electronic supplementary material, Material and methods.
A linear mixed-effects model was used to estimate variance components between and within rookeries. The model included rookery and clutch nested within rookery as random factors. Comparing total variances explained by the full model with a model having one factor removed tested the significance of each level in the analysis. Variance components were calculated using REML, and the VarCorr function in the R package nlme [40]. All ΔCq and ΔΔCq values were log-transformed with 2–X, where X is the mean ΔCq or mean ΔΔCq value, prior to analysis.
Broad- and narrow-sense heritabilities of ΔCq (within-treatment) and ΔΔCq values (across-treatment) for hsp60, hsp70 and hsp90 were estimated using an ‘animal model’ in ASReml 3.0 [41,42], fitting separate animal models to each treatment (control and heat shock), with rookery, offspring identity and dam as random factors. REML likelihood-ratio tests (REML LRT [42]) were used to find the best fitting model, by starting with a fully saturated model and then systematically reducing the number of parameters. Spearman's rank correlation was used to estimate genetic correlation between the expression levels of the genes under each treatment (see also electronic supplementary material, Materials and methods). The statistical power (B) of the correlation tests were calculated using the R package pwr [43].
3. Results and discussion
(a) Clutch is the most important component of variation in basal and increased levels of expression for hsp60, hsp70 and hsp90
We detected significant variation in the hsp expression in response to heat stress. For both the procedural control (29°C, 3 h) and the heat-shock (36°C, 3 h) treatments, the variance components analysis revealed no significant differences in expression between rookeries (table 1). However, there was a significant proportion of variation in hsp expression among clutches within rookeries for all target genes (table 1). Approximately, 30–50% of all variation in hsp expression was due to clutch effects. This suggests that the phenotypic variance explained by clutch may be due to the genetic constitution of the offspring, as well as maternal effects such as yolk quality or the thermal environment during early embryogenesis (e.g. the environment in utero immediately prior to oviposition) [44]. Female C. caretta have been shown to alternate ‘active and stay warm’ and ‘passive and stay cool’ thermoregulation strategies to optimize reproductive output [45,46]. While any flow-on effects of thermoregulation during embryogenesis are yet to be documented in sea turtles, thermoregulatory behaviours in oviparous lizards (Bassiana duperreyi) can have significant effects on hatchling phenotypes, notably on offspring sex [47,48].
| hsp60 |
hsp70 |
hsp90 |
||||
|---|---|---|---|---|---|---|
| % of total variance | p-value | % of total variance | p-value | % of total variance | p-value | |
| ΔCq 29°C | ||||||
| rookery | 22.7 | 0.215 | 0.0 | 0.999 | 0.0 | 0.999 |
| clutch | 33.3 | <0.001 | 28.0 | 0.002 | 23.0 | 0.009 |
| ΔCq 36°C | ||||||
| rookery | 14.6 | 0.446 | 3.2 | 0.792 | 28.6 | 0.143 |
| clutch | 48.6 | <0.001 | 30.4 | <0.001 | 33.9 | <0.001 |
| ΔΔCq | ||||||
| rookery | 21.3 | 0.264 | 2.9 | 0.683 | 30.4 | 0.037 |
| clutch | 39.3 | <0.001 | 4.5 | 0.442 | 11.8 | 0.041 |
(b) There is geographical variation in the plasticity of expression for hsp90 in response to heat shock
Gene expression assays were conducted on 14 clutches from a temperate rookery (DHI: offspring N = 78) and four clutches from a sub-tropical rookery (BB: offspring N = 18). All offspring increased expression for all target genes in response to an acute heat stress, consistent with a previous study [17]. Expression of hsp90 increased 50.8-fold in embryos from the temperate rookery in response to the heat shock, and 19.3-fold in embryos from the sub-tropical rookery (figure 2). These fold change differences were reflected by the significant proportion of total variance in ΔΔCq between rookeries for hsp90 (table 1). Large fold-changes were also evident for hsp60 and hsp70 (figure 2), but the differences between rookeries for these genes were non-significant (table 1). As found in the within-treatment analyses for all genes, there was a significant proportion of variance in relative gene expression between clutches within rookeries for hsp60 and hsp90.
Figure 2. Expression levels of hsp60, hsp70 and hsp90 by C. caretta embryos from a temperate (DHI, hatched bars; n = 129–130 individuals per gene) and sub-tropical (BB, grey bars; n = 30 individuals per gene) rookery in response to an acute heat shock. All mRNA data are normalized to the reference gene 18s (values are clutch mean ± s.e.m.; N = 18). Significant differences (p < 0.05) between rookeries denoted with an asterisk.
Higher relative gene expression in hsp90 in the sub-tropical rookery is consistent with the nest temperature differences between the two study sites [37,49], and suggests that embryos developing on these beaches may be locally adapted to different nest temperatures. On Ascension Island in the South Atlantic Ocean, the black sand of Northeast Bay (NEB) averages 2.6°C warmer than the white sand of Long Beach (LB), and a common garden experiment revealed that offspring of green sea turtles (Chelonia mydas) nesting on these beaches were locally adapted at a fine spatial scale [50]. Although offspring from the NEB rookery exhibited greater tolerance to warmer temperatures, hatching success was significantly lower than at the LB rookery, which led the authors to suggest that embryonic thermotolerance had not evolved in concert with the 0.5°C warming documented at these rookeries over the past 150 years [50].
In our study, stronger evidence for local adaptation would require replication at the rookery level and knowledge of the genetic relationships between nesting females [51,52]. Nevertheless, that hsp90 has significant heritable plasticity in expression may provide C. caretta embryos with a mechanism for increasing their thermotolerance as a response to climate change. Further, the degree to which plasticity allows mean phenotypes to track changes in the thermal environment will determine whether nesting populations evolve to tolerate changes in nest temperature [53]. It should also be noted that studies on Drosophila suggest that populations that evolved in high temperature environments do not upregulate Hsps to the same degree when exposed to heat stress compared to populations from cooler habitats [54]. This suggests mechanisms other than Hsp expression are more important for coping with high temperatures in permanently high temperature environments in some taxa.
(c) Variation in hsp expression is heritable in Caretta caretta
The animal models used to estimate heritability (see Material and methods) considered the effects of rookery, offspring identity and dam on hsp expression. Rookery had no significant effect on phenotype and was thereafter excluded from all models. The REML LRT showed that offspring was the best fitting pedigree factor in the animal model for hsp70, and so narrow-sense heritability was estimated for this gene. For hsp60, narrow-sense heritability was estimable in the heat shock treatment at 36°C, but not in the procedural control treatment at 29°C (table 2). The estimates of narrow-sense heritability for expression of hsp60 and hsp70 were high relative to those expected of physiological traits for ectotherms [55]. In general, physiological traits (e.g. O2 consumption, resistance to heat stress) have mid-level heritabilities (approx. 0.33), similar to those for behavioural traits (e.g. approx. 0.30), greater than those for life-history traits (e.g. approx. 0.26), but much lower than those for morphological traits (e.g. approx. 0.46) [55].
| hsp60 | hsp70 | hsp90 | |
|---|---|---|---|
| ΔCq 29°C | 0.51 ± 0.11a | 0.88 ± 0.23 | 0.32 ± 0.12a |
| ΔCq 36°C | 0.79 ± 0.28 | 0.75 ± 0.27 | 0.21 ± 0.13a |
| ΔΔCq | 0.38 ± 0.21 | 0.20 ± 0.10 | 0.37 ± 0.10a |
Although the heritability estimates of the plastic response of all three hsp genes were lower than the heritability estimates of the relative expression of each gene within treatments, they were nevertheless significant for hsp70 and hsp90 (table 2). This is not an unexpected result, as heritability estimates for plastic responses of a trait are typically lower than the heritability estimate of the trait itself (reviewed by [12,56]). Taken together, our results show that populations have the potential to adapt through selection to increase the baseline expression levels of hsps (ΔCq) and also through selection on plasticity on hsp70 and hsp90 (ΔΔCq). Such a capacity to adapt to warmer and more variable climates is not universal, even in species with large population sizes. For example, low additive genetic variance and a non-significant response to selection for upper thermal limits was found in a marine copepod [57]; and in a freshwater live-bearing fish [58]. Rainforest restricted Drosophila also appear to lack genetic variance for desiccation [59], though significant evolutionary responses in desiccation resistance has been reported under mild desiccation stress [60].
The differences in heritability estimates between relative expression and plastic expression could be due to a cost–benefit trade-off [56,61] in which expression of hsp60 and hsp70 are regulated to minimize any detrimental effects of overexpression [13]. Overexpression of hsp70 compromises the fitness and survival of Drosophila larvae [61], but the rate at which developmental temperatures rise may affect how relative expression of hsp70 evolves in natural populations [62]. Hence, it is possible that moderate increases in relative expression of hsp60 and hsp70 (figure 2) produce a short-term thermotolerance response in embryonic C. caretta, but maximal levels of expression are tightly controlled to ensure embryonic survival.
For hsp90, only broad-sense heritability was detected. Relative to hsp60 and hsp70, the low, but significant heritability estimates for hsp90 were similar to estimates for other physiological traits in other species [55,56,62]. The broad-sense heritabilities we estimated do not exclude the possibility of additive genetic variance in the expression of these genes, but do suggest significant maternal, environmental or dominance effects. Clearly, model species such as Drosophila would be more amenable to separating these effects. However, as the heritability of relative expression of hsp90 has not been previously estimated in reptiles, future studies might examine whether increased expression levels of this hsp are consistent across clutches of individual females. For example, sampling from the multiple clutches laid by females within a single season could reveal if the heritability of plasticity in expression of hsp90 is characterized more by maternal or environmental effects.
Relative gene expression in the heat-shock treatment was positively correlated with relative gene expression in the procedural control treatment for hsp60 (Spearman's rank correlation coefficient rs = 0.911, n = 18, p ≤ 0.001), hsp70 (rs = 0.624, n = 18, p ≤ 0.05) and hsp90 (rs = 0.515, n = 18, p ≤ 0.05). Thus, clutches that express highly in the control treatment also express highly in the heat-shock treatment. These correlations suggest genes that increase expression in one thermal environment will also increase expression in another.
The basal expression levels of different genes within environments were not significantly correlated in the control treatment for hsp60 and hsp70 (rs = 0.385, n = 18, p = 0.114), hsp60 and hsp90 (rs = 0.259, n = 18, p = 0.300) or hsp70 and hsp90 (rs = 0.424, n = 18, p = 0.080). In the heat-shock treatment, expression levels for hsp60 and hsp70 were significantly correlated (rs = 0.488, n = 18, p = 0.040). Increased expression levels hsp60 and hsp90 were marginally non-significant (rs = 0.451, n = 18, p = 0.060). Finally, increased expression levels of hsp70 and hsp90 were not significantly correlated (rs = 0.152, n = 18, p = 0.550). These correlations are all positive and suggest that selection for increased basal expression levels in one gene does, in general, give rise to increased expression in other genes (although here our sample sizes were too low to detect significant effects in most cases; statistical power ranged from 0.092 to 0.487 for non-significant correlations). Stress-induced transcription of hsp genes is regulated by heat-shock transcription factors (HSFs), which bind to the heat shock element sequences present in the promoters of hsp genes [63,64]. There are multiple HSFs in vertebrates (HSF1–4), which may provide redundancy and specialization of stress signals, though HSF1 appears to be the primary regulator of heat-inducible hsp gene expression in most vertebrate cells [64,65].
4. Conclusion
This is the first study to examine geographical and genetic variation in a physiological response to heat stress in sea turtle embryos, and the heritability of that response at the molecular level for any reptile species. Evolutionary adaptation to changing environments requires traits to be heritable, and we have shown heritability of expression of heat shock genes under both control and heat-stress conditions, and that the plasticity of expression of these heat shock genes is also heritable. Expression levels of hsp60 and hsp70 were attributed to both maternal and additive effects; expression levels of hsp90 were characterized by maternal and environmental effects. Additionally, we found strong correlations between expression levels of all target genes and incubation environments, suggesting that an elevated level of baseline expression results in relatively higher expression under thermal stress. While the relevance of our results to survival of embryos under field conditions are equivocal, they provide strong evidence of molecular mechanisms for tolerating and/or adapting to higher incubation temperatures. Of pressing concern is that no matter how many routes to adaption can be demonstrated, the pace of contemporary climate change may outstrip the evolutionary potential of slow-breeding species.
Ethics
All research described here adhered to local guidelines, and all appropriate ethical approval and licences were obtained.
Data accessibility
The datasets supporting this article have been uploaded onto DRYAD and as part of the electronic supplementary material.
Author's contributions
J.N.T., W.J.K. and N.J.M. conceived and designed the research, and S.W. provided substantial funding. J.N.T. carried out the molecular work and data analysis, the latter with assistance from W.J.K. J.L.T. performed heritability estimates. J.N.T., W.J.K., J.L.T., O.B., S.W., M.G.M. and N.J.M. wrote the manuscript. All authors approved the manuscript for publication.
Competing interests
The authors have no competing interests.
Funding
This project was funded by the Western Australian Department of Parks and Wildlife (to J.N.T.), the ANZ Holsworth Wildlife Trust (to N.J.M and J.N.T.) and the School of Animal Biology at the University of Western Australia. Funding from the Australian Institute of Marine Science and the Commonwealth Scientific and Industrial Research Organisation also provided scholarship support to J.N.T.).
Acknowledgements
We thank David Holley and Chris McMonagle (Western Australian Department of Parks and Wildlife, DPaW) in Shark Bay, and Peter Barnes and Keely Markovina (DPaW, Ningaloo Turtle Program in Cape Range National Park) for facilitating site access and fieldwork support. Blair Bentley, Adelaide Bevilaqua, Christie Craig, Kasey Darts, Frankie Hoult and Frances Leng provided invaluable support in the field. Additionally, we are grateful to Sherralee Lukehurst and Yvette Hitchen for assistance with the molecular work.
