Investigating the molecular basis of local adaptation to thermal stress: population differences in gene expression across the transcriptome of the copepod Tigriopus californicus

Heat tolerance, oxidative stress response tuning and robust gene activation in early-stage Drosophila melanogaster embryos

In organisms with complex life cycles, life stages that are most susceptible to environmental stress may determine species persistence in the face of climate change. Early embryos of Drosophila melanogaster are particularly sensitive to acute heat stress, yet tropical embryos have higher heat tolerance than temperate embryos, suggesting adaptive variation in embryonic heat tolerance. We compared transcriptomic responses to heat stress among tropical and temperate embryos to elucidate the gene regulatory basis of divergence in embryonic heat tolerance. The transcriptomes of tropical and temperate embryos differed in both constitutive and heat-stress-induced responses of the expression of relatively few genes, including genes involved in oxidative stress. Most of the transcriptomic response to heat stress was shared among all embryos. Embryos shifted the expression of thousands of genes, including increases in the expression of heat shock genes, suggesting robust zygotic gene activation and demonstrating that, contrary to previous reports, early embryos are not transcriptionally silent. The involvement of oxidative stress genes corroborates recent reports on the critical role of redox homeostasis in coordinating developmental transitions. By characterizing adaptive variation in the transcriptomic basis of embryonic heat tolerance, this study is a novel contribution to the literature on developmental physiology and developmental genetics.

A rapid return to normal: temporal gene expression patterns following cold exposure in the bumble bee Bombus impatiens

Bumble bees are common in cooler climates and many species likely experience periodic exposure to very cold temperatures, but little is known about the temporal dynamics of cold response mechanisms following chill exposure, especially how persistent effects of cold exposure may facilitate tolerance of future events. To investigate molecular processes involved in the temporal response by bumble bees to acute cold exposure, we compared mRNA transcript abundance in Bombus impatiens workers exposed to 0°C for 75 min (inducing chill coma) and control bees maintained at a constant ambient temperature (28°C). We sequenced the 3' end of mRNA transcripts (TagSeq) to quantify gene expression in thoracic tissue of bees at several time points (0, 10, 30, 120 and 720 min) following cold exposure. Significant differences from control bees were only detectable within 30 min after the treatment, with most occurring at the 10 min recovery time point. Genes associated with gluconeogenesis and glycolysis were most notably upregulated, while genes related to lipid and purine metabolism were downregulated. The observed patterns of expression indicate a rapid recovery after chill coma, suggesting an acute differential transcriptional response during recovery from chill coma and return to baseline expression levels within an hour, with no long-term gene expression markers of this cold exposure. Our work highlights the functions and pathways important for acute cold recovery, provides an estimated time frame for recovery from cold exposure in bumble bees, and suggests that cold hardening may be less important for these heterothermic insects.

Contrasting responses of Thermocyclops crassus and T. oithonoides (Crustacea, Copepoda) to thermal stress

Thermal tolerance is a critical factor influencing the survival of living organisms. This study focuses on the thermal resistance of copepod species, Thermocyclops crassus (Fischer, 1853) and T. oithonoides (Sars G.O., 1863), with overlapping distribution ranges in Europe. Short-term heat shock experiments were conducted to assess the thermal resistance of these copepods, considering various temperature increments and exposure durations. Additionally, the study explored the influence of heat shock on egg sac shedding, a vital indicator of population dynamics. Results indicate that widely distributed T. crassus exhibits higher thermal tolerance compared to narrowly distributed T. oithonoides, with survival rates varying under different heat shock conditions. Furthermore, T. crassus demonstrated a quicker response in dropping egg sacs in response to thermal stress, suggesting a potential adaptive mechanism for the survival of adults. However, rapid egg sac droppings pose high risks for eggs facing unfavorable conditions. T. crassus, inhabiting environments with greater temperature fluctuations such as the littoral and pelagial zones, exhibited better survival mechanisms compared to T. oithonoides, which predominantly resides in the pelagic zone. The findings have implications for understanding copepod responses to global warming and thermal pollution. This research contributes insights into the adaptive strategies of thermophilic copepod species and their ecological consequences.

Comparative Transcriptome Analyses Provide New Insights into the Evolution of Divergent Thermal Resistance in Two Eel Gobies

Adaptation to thermal conditions in tidal mudflats always involves tolerating frequent fluctuations and often extreme environmental temperatures. Regulation of gene expression plays a fundamental role in the evolution of these thermal adaptations. To identify the key gene regulatory networks associated with the thermal adaptation, we investigated the capability of cold tolerance, as well as the transcriptomic changes under cold stress in two mudflat inhabitants (Odontamblyopus lacepedii and O. rebecca) with contrasting latitude affinity. Our results revealed a remarkable divergent capacity of cold tolerance (CTmin: 0.61 °C vs. 9.57 °C) between the two gobies. Analysis of transcriptomic changes under cold stress unveiled 193 differentially expressed genes exhibiting similar expression profiles across all tissues and species, including several classic metabolic and circadian rhythm molecules such as ACOD and CIART that may represent the core cold response machinery in eel gobies. Meanwhile, some genes show a unique expression spectrum in the more cold-tolerant O. lacepedii suggesting their roles in the enhanced cold tolerance and hence the extreme thermal adaptations. In addition, a weighted gene co-expression network analysis (WGCNA) revealed a subset of metabolic hub genes including MYH11 and LIPT2 showing distinct down-regulation in O. lacepedii when exposed to cold stress which highlights the role of reduced energy consumption in the enhanced cold tolerance of eel gobies. These findings not only provide new insights into how mudflat teleosts could cope with cold stress and their potential evolutionary strategies for adapting to their thermal environment, but also have important implications for sound management and conservation of their fishery resources in a scenario of global climate warming in the marine realm.

Multi-omics reveals largely distinct transcript- and protein-level responses to the environment in an intertidal mussel

Organismal responses to stressful environments are influenced by numerous transcript- and protein-level mechanisms, and the relationships between expression changes at these levels are not always straightforward. Here, we used paired transcriptomic and proteomic datasets from two previous studies from gill of the California mussel, Mytilus californianus, to explore how simultaneous transcript and protein abundance patterns may diverge under different environmental scenarios. Field-acclimatized mussels were sampled from two disparate intertidal sites; individuals from one site were subjected to three further treatments (common garden, low-intertidal or high-intertidal outplant) that vary in temperature and feeding time. Assessing 1519 genes shared between the two datasets revealed that both transcript and protein expression patterns differentiated the treatments at a global level, despite numerous underlying discrepancies. There were far more instances of differential expression between treatments in transcript only (1451) or protein only (226) than of the two levels shifting expression concordantly (68 instances). Upregulated expression of cilium-associated transcripts (likely related to feeding) was associated with relatively benign field treatments. In the most stressful treatment, transcripts, but not proteins, for several molecular chaperones (including heat shock proteins and endoplasmic reticulum chaperones) were more abundant, consistent with a threshold model for induction of translation of constitutively available mRNAs. Overall, these results suggest that the relative importance of transcript- and protein-level regulation (translation and/or turnover) differs among cellular functions and across specific microhabitats or environmental contexts. Furthermore, the degree of concordance between transcript and protein expression can vary across benign versus acutely stressful environmental conditions.

Loss of mitochondrial performance at high temperatures is correlated with upper thermal tolerance among populations of an intertidal copepod

Environmental temperatures have pervasive effects on the performance and tolerance of ectothermic organisms, and thermal tolerance limits likely play key roles underlying biogeographic ranges and responses to environmental change. Mitochondria are central to metabolic processes in eukaryotic cells, and these metabolic functions are thermally sensitive; however, potential relationships between mitochondrial function, thermal tolerance limits and local thermal adaptation in general remain unresolved. Loss of ATP synthesis capacity at high temperatures has recently been suggested as a mechanistic link between mitochondrial function and upper thermal tolerance limits. Here we use a common-garden experiment with seven locally adapted populations of intertidal copepods (Tigriopus californicus), spanning approximately 21.5° latitude, to assess genetically based variation in the thermal performance curves of maximal ATP synthesis rates in isolated mitochondria. These thermal performance curves displayed substantial variation among populations with higher ATP synthesis rates at lower temperatures (20-25 °C) in northern populations than in southern populations. In contrast, mitochondria from southern populations maintained ATP synthesis rates at higher temperatures than the temperatures that caused loss of ATP synthesis capacity in mitochondria from northern populations. Additionally, there was a tight correlation between the thermal limits of ATP synthesis and previously determined variation in upper thermal tolerance limits among populations. This suggests that mitochondria may play an important role in latitudinal thermal adaptation in T. californicus, and supports the hypothesis that loss of mitochondrial performance at high temperatures is linked to whole-organism thermal tolerance limits in this ectotherm.

Transcriptomic and functional genetic evidence for distinct ecophysiological responses across complex life cycle stages

Organisms with complex life cycles demonstrate a remarkable ability to change their phenotypes across development, presumably as an evolutionary adaptation to developmentally variable environments. Developmental variation in environmentally sensitive performance, and thermal sensitivity in particular, has been well documented in holometabolous insects. For example, thermal performance in adults and juvenile stages exhibit little genetic correlation (genetic decoupling) and can evolve independently, resulting in divergent thermal responses. Yet, we understand very little about how this genetic decoupling occurs. We tested the hypothesis that genetic decoupling of thermal physiology is driven by fundamental differences in physiology between life stages, despite a potentially conserved Cellular Stress Response. We used RNAseq to compare transcript expression in response to a cold stressor in Drosophila melanogaster larvae and adults and used RNAi (RNA interference) to test whether knocking down nine target genes differentially affected larval and adult cold tolerance. Transcriptomic responses of whole larvae and adults during and following exposure to -5°C were largely unique both in identity of responding transcripts and in temporal dynamics. Further, we analyzed the tissue-specificity of differentially-expressed transcripts from FlyAtlas 2 data, and concluded that stage-specific differences in transcription were not simply driven by differences in tissue composition. In addition, RNAi of target genes resulted in largely stage-specific and sometimes sex-specific effects on cold tolerance. The combined evidence suggests that thermal physiology is largely stage-specific at the level of gene expression, and thus natural selection may be acting on different loci during the independent thermal adaptation of different life stages.

First transcriptome of the copepod Gladioferens pectinatus subjected to chronic contaminant exposures

Contaminants are often at low concentrations in ecosystems and their effects on exposed organisms can occur over long periods of time and across multiple generations. Alterations to subcellular mechanistic pathways in response to exposure to contaminants can provide insights into mechanisms of toxicity that methods measuring higher levels of biological may miss. Analysis of the whole transcriptome can identify novel mechanisms of action leading to impacts in exposed biota. The aim of this study was to characterise how exposures to copper, benzophenone and diclofenac across multiple generations altered molecular expression pathways in the marine copepod Gladioferens pectinatus. Results of the study demonstrated differential gene expression was observed in cultures exposure to diclofenac (569), copper (449) and benzophenone (59). Pathways linked to stress, growth, cellular and metabolic processes were altered by exposure to all three contaminants with genes associated with oxidative stress and xenobiotic regulation also impacted. Protein kinase functioning, cytochrome P450, transcription, skeletal muscle contraction/relaxation, mitochondrial phosphate translocator, protein synthesis and mitochondrial methylation were all differentially expressed with all three chemicals. The results of the study also suggested that using dimethyl sulfoxide as a dispersant influenced the transcriptome and future research may want to investigate it's use in molecular studies. Data generated in this study provides a first look at transcriptomic response of G. pectinatus exposed to contaminants across multiple generations, future research is needed to validate the identified biomarkers and link these results to apical responses such as population growth to demonstrate the predictive capacity of molecular tools.

Nature versus nurture in heat stress induced learning between inbred and outbred populations of Lymnaea stagnalis

Changing environmental conditions often lead to microevolution of traits that are adaptive under the current selection pressure. Currently, one of the major selection pressures is the rise in temperatures globally that has a severe impact on the behavioral ecology of animals. However, the role of thermal stress on neuronal plasticity and memory formation is not well understood. Thermal tolerance and sensitivity to heat stress show variation across populations of the same species experiencing different thermal regimes. We used two populations of the pond snail Lymnaea stagnalis: one lab-bred W-snails and the other wild Delta snails to test heat shock induced learning and memory formation for the Garcia effect learning paradigm. In Garcia effect, a single pairing of a heat stressor (30 °C for 1h) with a novel taste results in a taste-specific negative hedonic shift lasting 24h as long-term memory (LTM) in lab bred W-snails. In this study we used a repeated heat stress procedure to test for increased or decreased sensitivity to the heat before testing for the Garcia effect. We found that lab-bred W-snails show increased sensitivity to heat stress after repeated heat exposure for 7days, leading to enhanced LTM for Garcia effect with only 15min of heat exposure instead of standard 1h. Surprisingly, the freshly collected wild snails do not show Garcia effect. Additionally, F1 generation of wild snails raised and maintained under laboratory conditions still retain their heat stress tolerance similar to their parents and do not show a Garcia effect under standard learning paradigm or even after repeated heat stressor. Thus, we found a differential effect of heat stress on memory formation in wild and lab bred snails. Most interestingly we also show that local environmental (temperature) conditions for one generation is not enough to alter thermal sensitivity in a wild population of L. stagnalis.

Elevated Salinity Rapidly Confers Cross-Tolerance to High Temperature in a Splash-Pool Copepod

Accurate forecasting of organismal responses to climate change requires a deep mechanistic understanding of how physiology responds to present-day variation in the physical environment. However, the road to physiological enlightenment is fraught with complications: predictable environmental fluctuations of any single factor are often accompanied by substantial stochastic variation and rare extreme events, and several factors may interact to affect physiology. Lacking sufficient knowledge of temporal patterns of co-variation in multiple environmental stressors, biologists struggle to design and implement realistic and relevant laboratory experiments. In this study, we directly address these issues, using measurements of the thermal tolerance of freshly collected animals and long-term field records of environmental conditions to explore how the splash-pool copepod Tigriopus californicus adjusts its physiology as its environment changes. Salinity and daily maximum temperature-two dominant environmental stressors experienced by T. californicus-are extraordinarily variable and unpredictable more than 2-3 days in advance. However, they substantially co-vary such that when temperature is high salinity is also likely to be high. Copepods appear to take advantage of this correlation: median lethal temperature of field-collected copepods increases by 7.5°C over a roughly 120 parts-per-thousand range of ambient salinity. Complementary laboratory experiments show that exposure to a single sublethal thermal event or to an abrupt shift in salinity also elicits rapid augmentation of heat tolerance via physiological plasticity, although the effect of salinity dwarfs that of temperature. These results suggest that T. californicus's physiology keeps pace with the rapid, unpredictable fluctuations of its hypervariable physical environment by responding to the cues provided by recent sublethal stress and, more importantly, by leveraging the mechanistic cross-talk between responses to salinity and heat stress.

Modeling the impact of temperature on the population abundance of the ambrosia beetle Xyleborus affinis (Curculionidae: Scolytinae) under laboratory-reared conditions

Modeling the impact of temperature on each life stage of a beetle population represents a continuing challenge. This study evaluates the effects of five temperature treatments (20, 23, 26, 29 and 32 °C) on population abundance and timing of a colony of ambrosia beetles Xyleborus affinis reared under laboratory conditions and use this data to develop demographic and phenological models. Abundances at each life stage (eggs, larvae, pupae and adult) were examined through periodic destructive sampling; given that it was not possible to track individuals. To assess the effects of temperature on oviposition, development and survival rates we developed a novel estimation strategy based on cohorts, which does not require individual developmental data. Since oviposition was entirely unwitnessed, we assessed competing empirical ovipositional models. Rates of development were computed using a modal rate curve for each life stage, and rates were projected to cohorts in life stages assuming log-normal developmental variance. Temperature-driven survival rates were assumed to be logistic with a quadratic exponent to capture modal temperature dependence. Parameters were estimated simultaneously using minimum negative log posterior likelihood, assuming Poisson distribution of observations and using priors to inform unobserved developmental rates and enforce mechanistic constraints on oviposition models. A parabolic function best described oviposition rate. Optimal developmental temperatures were 30.5 °C, 29 °C and 27.5 °C, with maximum developmental rates of 0.26/day, 0.12/day and 0.23/day for eggs, larvae and pupae, respectively. The survival rates in the range 20-29 °C were equal to 1 in the eggs-to-larvae transition, from 0.72 to 0.35 in larvae-to-pupae transition, and from 0.2 to 0.89 in pupae-to-adults transition. This procedure effectively characterized the direct thermal effects on development and survival of each life stage in the X. affinis under laboratory conditions and would be suitable for estimating temperature dependence for other species in which individual observations are not possible.

Molecular adaptations to heat stress in the thermophilic ant genus Cataglyphis

Over the last decade, increasing attention has been paid to the molecular adaptations used by organisms to cope with thermal stress. However, to date, few studies have focused on thermophilic species living in hot, arid climates. In this study, we explored molecular adaptations to heat stress in the thermophilic ant genus Cataglyphis, one of the world's most thermotolerant animal taxa. We compared heat tolerance and gene expression patterns across six Cataglyphis species from distinct phylogenetic groups that live in different habitats and experience different thermal regimes. We found that all six species had high heat tolerance levels with critical thermal maxima (CT_max ) ranging from 43℃ to 45℃ and a median lethal temperature (LT50) ranging from 44.5℃ to 46.8℃. Transcriptome analyses revealed that, although the number of differentially expressed genes varied widely for the six species (from 54 to 1118), many were also shared. Functional annotation of the differentially expressed and co-expressed genes showed that the biological pathways involved in heat-shock responses were similar among species and were associated with four major processes: the regulation of transcriptional machinery and DNA metabolism; the preservation of proteome stability; the elimination of toxic residues; and the maintenance of cellular integrity. Overall, our results suggest that molecular responses to heat stress have been evolutionarily conserved in the ant genus Cataglyphis and that their diversity may help workers withstand temperatures close to their physiological limits.

Transcriptomic signatures of temperature adaptation in the eastern oyster Crassostrea virginica

The large geographic distribution of the eastern oyster, Crassostrea virginica, makes it an ideal species to test how populations have adapted to latitudinal gradients in temperature. Despite inhabiting distinct thermal regimes, populations of C. virginica near the species' southern and northern geographic range show no population differences in their physiological response to temperature. In this study, we used comparative transcriptomics to understand how oysters from either end of the species' range maintain enantiostasis across three acclimation temperatures (10, 20, and 30°C). With this approach, we identified genes that were differentially expressed in response to temperature between individuals of C. virginica collected from New Brunswick, Canada and Louisiana, USA. We observed a core set of genes whose expression responded to temperature in both populations, but also an even larger set of genes with expression patterns that were unique to each population. Intriguingly, the genes with population-specific responses to temperature had elevated F_ST and Ka/Ks ratios compared to the genome-wide average. In contrast, genes showing only a response to temperature were found to only have elevated F_ST values suggesting that divergent F_ST may be due to selection on linked regulatory regions rather than positive selection on protein coding regions. Taken together, our results suggest that, despite coarse-scale physiological similarities, natural selection has shaped divergent gene expression responses to temperature in geographically separated populations of this broadly eurythermal marine invertebrate.

Population differences in Chinook salmon (Oncorhynchus tshawytscha) DNA methylation: Genetic drift and environmental factors

Local adaptation and phenotypic differences among populations have been reported in many species, though most studies focus on either neutral or adaptive genetic differentiation. With the discovery of DNA methylation, questions have arisen about its contribution to individual variation in and among natural populations. Previous studies have identified differences in methylation among populations of organisms, although most to date have been in plants and model animal species. Here we obtained eyed eggs from eight populations of Chinook salmon (Oncorhynchus tshawytscha) and assayed DNA methylation at 23 genes involved in development, immune function, stress response, and metabolism using a gene-targeted PCR-based assay for next-generation sequencing. Evidence for population differences in methylation was found at eight out of 23 gene loci after controlling for developmental timing in each individual. However, we found no correlation between freshwater environmental parameters and methylation variation among populations at those eight genes. A weak correlation was identified between pairwise DNA methylation dissimilarity among populations and pairwise F ST based on 15 microsatellite loci, indicating weak effects of genetic drift or geographic distance on methylation. The weak correlation was primarily driven by two genes, GTIIBS and Nkef. However, single-gene Mantel tests comparing methylation and pairwise F ST were not significant after Bonferroni correction. Thus, population differences in DNA methylation are more likely related to unmeasured oceanic environmental conditions, local adaptation, and/or genetic drift. DNA methylation is an additional mechanism that contributes to among population variation, with potential influences on organism phenotype, adaptive potential, and population resilience.

Long-term thermal acclimation drives adaptive physiological adjustments of a marine gastropod to reduce sensitivity to climate change

Ocean warming is predicted to challenge the persistence of a variety of marine organisms, especially when combined with ocean acidification. While temperature affects virtually all physiological processes, the extent to which thermal history mediates the adaptive capacity of marine organisms to climate change has been largely overlooked. Using populations of a marine gastropod (Turbo undulatus) with different thermal histories (cool vs. warm), we compared their physiological adjustments following exposure (8-week) to ocean acidification and warming. Compared to cool-acclimated counterparts, we found that warm-acclimated individuals had a higher thermal threshold (i.e. increased CT_max by 2 °C), which was unaffected by the exposure to ocean acidification and warming. Thermal history also strongly mediated physiological effects, where warm-acclimated individuals adjusted to warming by conserving energy, suggested by lower respiration and ingestion rates, energy budget (i.e. scope for growth) and O:N ratio. After exposure to warming, warm-acclimated individuals had higher metabolic rates and greater energy budget due to boosted ingestion rates, but such compensatory feeding disappeared when combined with ocean acidification. Overall, we suggest that thermal history can be a critical mediator of physiological performance under future climatic conditions. Given the relatively gradual rate of global warming, marine organisms may be better able to adaptively adjust their physiology to future climate than what short-term experiments currently convey.

Species and population specific gene expression in blood transcriptomes of marine turtles

BACKGROUND

Transcriptomic data has demonstrated utility to advance the study of physiological diversity and organisms' responses to environmental stressors. However, a lack of genomic resources and challenges associated with collecting high-quality RNA can limit its application for many wild populations. Minimally invasive blood sampling combined with de novo transcriptomic approaches has great potential to alleviate these barriers. Here, we advance these goals for marine turtles by generating high quality de novo blood transcriptome assemblies to characterize functional diversity and compare global transcriptional profiles between tissues, species, and foraging aggregations.

RESULTS

We generated high quality blood transcriptome assemblies for hawksbill (Eretmochelys imbricata), loggerhead (Caretta caretta), green (Chelonia mydas), and leatherback (Dermochelys coriacea) turtles. The functional diversity in assembled blood transcriptomes was comparable to those from more traditionally sampled tissues. A total of 31.3% of orthogroups identified were present in all four species, representing a core set of conserved genes expressed in blood and shared across marine turtle species. We observed strong species-specific expression of these genes, as well as distinct transcriptomic profiles between green turtle foraging aggregations that inhabit areas of greater or lesser anthropogenic disturbance.

CONCLUSIONS

Obtaining global gene expression data through non-lethal, minimally invasive sampling can greatly expand the applications of RNA-sequencing in protected long-lived species such as marine turtles. The distinct differences in gene expression signatures between species and foraging aggregations provide insight into the functional genomics underlying the diversity in this ancient vertebrate lineage. The transcriptomic resources generated here can be used in further studies examining the evolutionary ecology and anthropogenic impacts on marine turtles.

Transcriptomics Analysis and Re-sequencing Reveal the Mechanism Underlying the Thermotolerance of an Artificial Selection Population of the Pacific Oyster

The Pacific oyster is a globally important aquaculture species inhabiting the intertidal environment, which experiences great temperature variation. Mass deaths in the summer pose a major challenge for the oyster industry. We initiated an artificial selection breeding program in 2017 using acute heat shock treatments of the parents to select for thermotolerance in oysters. In this study, we compared the respiration rate, summer survival rate, gene expression, and gene structure of F₂ selected oysters and non-selected wild oysters. A transcriptional analysis revealed global divergence between the selected and control groups at the larval stage, including 4764 differentially expressed genes, among which 79 genes were heat-responsive genes. Five heat shock proteins were enriched, and four of the six genes (five heat stock genes in the enriched GO terms and KEGG pathways and BAG4) were differentially expressed in 1-year-old oysters. Integration of the transcriptomic and re-sequencing data of the selected and the control groups revealed 1090 genes that differentiated in both gene structure and expression. Two SNPs (single nucleotide polymorphism) that may mediate the expression of CGI_10022585 and CGI_10024709 were validated. In addition, the respiration rate of 1-year-old oysters varied significantly between the selected group and the control group at room temperature (20°C). And the summer survival rate of the selected population was significantly improved. This study not only shows that artificial selection has a significant effect on the gene structure and expression of oysters, but it also helps reveal the mechanism underlying their tolerance of high temperature as well as the ability of oysters to adapt to climate change.

A complete enzymatic capacity for biosynthesis of docosahexaenoic acid (DHA, 22 : 6n-3) exists in the marine Harpacticoida copepod Tigriopus californicus

The long-standing paradigm establishing that global production of Omega-3 (n–3) long-chain polyunsaturated fatty acids (LC-PUFA) derived almost exclusively from marine single-cell organisms, was recently challenged by the discovery that multiple invertebrates possess methyl-end (or ωx) desaturases, critical enzymes enabling the biosynthesis of n–3 LC-PUFA. However, the question of whether animals with ωx desaturases have complete n–3 LC-PUFA biosynthetic pathways and hence can contribute to the production of these compounds in marine ecosystems remained unanswered. In the present study, we investigated the complete enzymatic complement involved in the n–3 LC-PUFA biosynthesis in Tigriopus californicus, an intertidal harpacticoid copepod. A total of two ωx desaturases, five front-end desaturases and six fatty acyl elongases were successfully isolated and functionally characterized. The T. californicus ωx desaturases enable the de novo biosynthesis of C₁₈ PUFA such as linoleic and α-linolenic acids, as well as several n–3 LC-PUFA from n–6 substrates. Functions demonstrated in front-end desaturases and fatty acyl elongases unveiled various routes through which T. californicus can biosynthesize the physiologically important arachidonic and eicosapentaenoic acids. Moreover, T. californicus possess a Δ4 desaturase, enabling the biosynthesis of docosahexaenoic acid via the ‘Δ4 pathway’. In conclusion, harpacticoid copepods such as T. californicus have complete n–3 LC-PUFA biosynthetic pathways and such capacity illustrates major roles of these invertebrates in the provision of essential fatty acids to upper trophic levels.

Tolerance Patterns and Transcriptomic Response to Extreme and Fluctuating Salinities across Populations of the Intertidal Copepod Tigriopus californicus

AbstractPopulations that tolerate extreme environmental conditions with frequent fluctuations can give valuable insights into physiological limits and adaptation. In some estuarine and marine ecosystems, organisms must adapt to extreme and fluctuating salinities, but not much is known about how varying salinities impact local adaptation across a wide geographic range. We used eight geographically and genetically divergent populations of the intertidal copepod Tigriopus californicus to test whether northern populations have greater tolerance to low salinity stresses, as they experience greater precipitation and less evaporation. We used a common-garden experiment approach and exposed all populations to acute low (1 and 3 ppt) and high (110 and 130 ppt) salinities for 24 h and to a fluctuation between baseline salinity and moderate low (7 ppt) and high (80 ppt) salinities for over 49 h. We also performed RNA sequencing at several time points during the fluctuation between baseline and salinity of 7 ppt to understand the molecular basis of divergence between two populations with differing physiological responses. We present these novel findings: (1) acute low salinity conditions caused more deaths than high salinity; (2) molecular processes that elevate proline levels increased in salinity of 7 ppt, which contrasts with other physiological studies in T. californicus that mainly associated accumulation of proline with hyperosmotic stress; and (3) tolerance to a salinity fluctuation did not follow a latitudinal trend but was instead governed by a complex interplay of factors, including population and duration of salinity stress. This highlights the importance of including a wider variety of environmental conditions in empirical studies to understand local adaptation.

Different amplitudes of temperature fluctuation induce distinct transcriptomic and metabolomic responses in the dung beetle Phanaeus vindex

Most studies exploring molecular and physiological responses to temperature have focused on constant temperature treatments. To gain a better understanding of the impact of fluctuating temperatures, we investigated the effects of increased temperature variation on Phanaeus vindex dung beetles across levels of biological organization. Specifically, we hypothesized that increased temperature variation is energetically demanding. We predicted that thermal sensitivity of metabolic rate and energetic reserves would be reduced with increasing fluctuation. To test this, we examined the responses of dung beetles to constant (20°C), low fluctuation (20±5°C), or high fluctuation (20±12°C) temperature treatments using respirometry, assessment of energetic reserves and HPLC-MS-based metabolomics. We found no significant differences in metabolic rate or energetic reserves, suggesting increased fluctuations were not energetically demanding. To understand why there was no effect of increased amplitude of temperature fluctuation on energetics, we assembled and annotated a de novo transcriptome, finding non-overlapping transcriptomic and metabolomic responses of beetles exposed to different fluctuations. We found that 58 metabolites increased in abundance in both fluctuation treatments, but 15 only did so in response to high-amplitude fluctuations. We found that 120 transcripts were significantly upregulated following acclimation to any fluctuation, but 174 were upregulated only in beetles from the high-amplitude fluctuation treatment. Several differentially expressed transcripts were associated with post-translational modifications to histones that support a more open chromatin structure. Our results demonstrate that acclimation to different temperature fluctuations is distinct and may be supported by increasing transcriptional plasticity. Our results indicate for the first time that histone modifications may underlie rapid acclimation to temperature variation.

Allele-Specific Expression and Evolution of Gene Regulation Underlying Acute Heat Stress Response and Local Adaptation in the Copepod Tigriopus californicus

Geographic variation in environmental temperature can select for local adaptation among conspecific populations. Divergence in gene expression across the transcriptome is a key mechanism for evolution of local thermal adaptation in many systems, yet the genetic mechanisms underlying this regulatory evolution remain poorly understood. Here we examine gene expression in 2 locally adapted Tigriopus californicus populations (heat tolerant San Diego, SD, and less tolerant Santa Cruz, SC) and their F1 hybrids during acute heat stress response. Allele-specific expression (ASE) in F1 hybrids was used to determine cis-regulatory divergence. We found that the number of genes showing significant allelic imbalance increased under heat stress compared to unstressed controls. This suggests that there is significant population divergence in cis-regulatory elements underlying heat stress response. Specifically, the number of genes showing an excess of transcripts from the more thermal tolerant (SD) population increased with heat stress while that number of genes with an SC excess was similar in both treatments. Inheritance patterns of gene expression also revealed that genes displaying SD-dominant expression phenotypes increase in number in response to heat stress; that is, across loci, gene expression in F1's following heat stress showed more similarity to SD than SC, a pattern that was absent in the control treatment. The observed patterns of ASE and inheritance of gene expression provide insight into the complex processes underlying local adaptation and thermal stress response.

Heat tolerance and thermal preference of the copepod Tigriopus californicus are insensitive to ecologically relevant dissolved oxygen levels

Shifting climate patterns may impose novel combinations of abiotic conditions on animals, yet understanding of the present-day interactive effects of multiple stressors remains under-developed. We tested the oxygen and capacity limited thermal tolerance (OCLTT) hypothesis and quantified environmental preference of the copepod Tigriopus californicus, which inhabits rocky-shore splashpools where diel fluctuations of temperature and dissolved oxygen (DO) are substantial. Egg-mass bearing females were exposed to a 5 h heat ramp to peak temperatures of 34.1–38.0 °C crossed with each of four oxygen levels: 22, 30, 100 and 250% saturation (4.7–5.3, 5.3–6.4, 21.2–21.3, and 50.7–53.3 kPa). Survival decreased at higher temperatures but was independent of DO. The behavioral preference of females was quantified in seven combinations of gradients of both temperature (11–37 °C) and oxygen saturation (17–206% or 3.6–43.6 kPa). Females avoided high temperatures regardless of DO levels. This pattern was more pronounced when low DO coincided with high temperature. In uniform temperature treatments, the distribution shifted toward high DO levels, especially in uniform high temperature, confirming that Tigriopus can sense environmental pO₂. These results question the ecological relevance of OCLTT for Tigriopus and raise the possibility of microhabitat selection being used within splashpool environments to avoid physiologically stressful combinations of conditions.

Sex differences in early transcriptomic responses to oxidative stress in the copepod Tigriopus californicus

Background

Patterns of gene expression can be dramatically different between males and females of the same species, in part due to genes on sex chromosomes. Here we test for sex differences in early transcriptomic response to oxidative stress in a species which lacks heteromorphic sex chromosomes, the copepod Tigriopus californicus.

Results

Male and female individuals were separately exposed to control conditions and pro-oxidant conditions (hydrogen peroxide and paraquat) for periods of 3 hours and 6 hours. Variance partitioning showed the greatest expression variance among individuals, highlighting the important information that can be obscured by the common practice of pooling individuals. Gene expression variance between sexes was greater than that among treatments, showing the profound effect of sex even when males and females share the same genome. Males exhibited a larger response to both pro-oxidants, differentially expressing more than four times as many genes, including up-regulation of more antioxidant genes, heat shock proteins and protease genes. While females differentially expressed fewer genes, the magnitudes of fold change were generally greater, indicating a more targeted response. Although females shared a smaller fraction of differentially expressed genes between stressors and time points, expression patterns of antioxidant and protease genes were more similar between stressors and more GO terms were shared between time points.

Conclusions

Early transcriptomic responses to the pro-oxidants H₂O₂ and paraquat in copepods revealed substantial variation among individuals and between sexes. The finding of such profound sex differences in oxidative stress response, even in the absence of sex chromosomes, highlights the importance of studying both sexes and the potential for developing sex-specific strategies to promote optimal health and aging in humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07179-5.

Regulation Between HSF1 Isoforms and HSPs Contributes to the Variation in Thermal Tolerance Between Two Oyster Congeners

Heat shock transcription factor 1 (HSF1) plays an important role in regulating heat shock, which can activate heat shock proteins (HSPs). HSPs can protect organisms from thermal stress. Oysters in the intertidal zone can tolerate thermal stress. The Pacific oyster (Crassostrea gigas gigas) and Fujian oyster (C. gigas angulata)—allopatric subspecies with distinct thermal tolerances—make good study specimens for analyzing and comparing thermal stress regulation. We cloned and compared HSF1 isoforms, which is highly expressed under heat shock conditions in the two subspecies. The results revealed that two isoforms (HSF1a and HSF1d) respond to heat shock in both Pacific and Fujian oysters, and different heat shock conditions led to various combinations of isoforms. Subcellular localization showed that isoforms gathered in the nucleus when exposed to heat shock. The co-immunoprecipitation revealed that HSF1d can be a dimer. In addition, we selected HSPs that are expressed under the heat shock response, according to the RNA-seq and proteomic analyses. For the HSPs, we analyzed the coding part and the promoter sequences. The result showed that the domains of HSPs are conserved in two subspecies, but the promoters are significantly different. The Dual-Luciferase assay showed that the induced expression isoform HSF1d had the highest activity in C. gigas gigas, while the constitutively-expressed HSF1a was most active in C. gigas angulata. In addition, variation in the level of HSP promoters appeared to be correlated with gene expression. We argue that this gene is regulated based on the different expression levels between the two subspecies’ responses to heat shock. In summary, various stress conditions can yield different HSF1 isoforms and respond to heat shock in both oyster subspecies. Differences in how the isoforms and promoter are activated may contribute to their differential expressions. Overall, the results comparing C. gigas gigas and C. gigas angulata suggest that these isoforms have a regulatory relationship under heat shock, providing valuable information on the thermal tolerance mechanism in these commercially important oyster species.

Biogeographic parallels in thermal tolerance and gene expression variation under temperature stress in a widespread bumble bee

Global temperature changes have emphasized the need to understand how species adapt to thermal stress across their ranges. Genetic mechanisms may contribute to variation in thermal tolerance, providing evidence for how organisms adapt to local environments. We determine physiological thermal limits and characterize genome-wide transcriptional changes at these limits in bumble bees using laboratory-reared Bombus vosnesenskii workers. We analyze bees reared from latitudinal (35.7-45.7°N) and altitudinal (7-2154 m) extremes of the species' range to correlate thermal tolerance and gene expression among populations from different climates. We find that critical thermal minima (CTMIN) exhibit strong associations with local minimums at the location of queen origin, while critical thermal maximum (CTMAX) was invariant among populations. Concordant patterns are apparent in gene expression data, with regional differentiation following cold exposure, and expression shifts invariant among populations under high temperatures. Furthermore, we identify several modules of co-expressed genes that tightly correlate with critical thermal limits and temperature at the region of origin. Our results reveal that local adaptation in thermal limits and gene expression may facilitate cold tolerance across a species range, whereas high temperature responses are likely constrained, both of which may have implications for climate change responses of bumble bees.

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

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

Transcriptomic signatures of cold adaptation and heat stress in the winter ant (Prenolepis imparis)

Climate change is a serious threat to biodiversity; it is therefore important to understand how animals will react to this stress. Ectotherms, such as ants, are especially sensitive to the climate as the environmental temperature influences myriad aspects of their biology, from optimal foraging time to developmental rate. In this study, we conducted an RNA-seq analysis to identify stress-induced genes in the winter ant (Prenolepis imparis). We quantified gene expression during heat and cold stress relative to a control temperature. From each of our conditions, we sequenced the transcriptome of three individuals. Our de novo assembly included 13,324 contigs that were annotated against the nr and SwissProt databases. We performed gene ontology and enrichment analyses to gain insight into the physiological processes involved in the stress response. We identified a total of 643 differentially expressed genes across both treatments. Of these, only seven genes were differentially expressed in the cold-stressed ants, which could indicate that the temperature we chose for trials did not induce a strong stress response, perhaps due to the cold adaptations of this species. Conversely, we found a strong response to heat: 426 upregulated genes and 210 downregulated genes. Of these, ten were expressed at a greater than ten-fold change relative to the control. The transcripts we could identify included those encoding for protein folding genes, heat shock proteins, histones, and Ca2+ ion transport. One of these transcripts, hsc70-4L was found to be under positive selection. We also characterized the functional categories of differentially expressed genes. These candidate genes may be functionally conserved and relevant for related species that will deal with rapid climate change.

Temperature-dependent life history and transcriptomic responses in heat-tolerant versus heat-sensitive Brachionus rotifers

Thermal stress response is an essential physiological trait that determines occurrence and temporal succession in nature, including response to climate change. We compared temperature-related demography in closely related heat-tolerant and heat-sensitive Brachionus rotifer species. We found significant differences in heat response, with the heat-sensitive species adopting a strategy of long survival and low population growth, while the heat-tolerant followed the opposite strategy. In both species, we examined the genetic basis of physiological variation by comparing gene expression across increasing temperatures. Comparative transcriptomic analyses identified shared and opposing responses to heat. Interestingly, expression of heat shock proteins (hsps) was strikingly different in the two species and mirrored differences in population growth rates, showing that hsp genes are likely a key component of a species' adaptation to different temperatures. Temperature induction caused opposing patterns of expression in further functional categories including energy, carbohydrate and lipid metabolism, and in genes related to ribosomal proteins. In the heat-sensitive species, elevated temperatures caused up-regulation of genes related to meiosis induction and post-translational histone modifications. This work demonstrates the sweeping reorganizations of biological functions that accompany temperature adaptation in these two species and reveals potential molecular mechanisms that might be activated for adaptation to global warming.

Latitudinal variation in norms of reaction of phenology in the greater duckweed Spirodela polyrhiza

Variable environments may result in the evolution of adaptive phenotypic plasticity when cues reliably indicate an appropriate phenotype-environment match. Although adaptive plasticity is well established for phenological traits expressed across environments, local differentiation in norms of reaction is less well studied. The switch from the production of regular fronds to overwintering 'turions' in the greater duckweed Spirodela polyrhiza is vital to fitness and is expressed as a norm of reaction induced by falling temperatures associated with the onset of winter. However, the optimal norm of reaction to temperature is expected to differ across latitudes. Here, we test the hypothesis that a gradient in the length and predictability of growing seasons across latitudes results in the evolution of reaction norms characterized by earlier turion production at higher latitudes. We test this by collecting S. polyrhiza from replicate populations across seven latitudes from Ontario to Florida and then assessing differentiation in thermal reaction norms of turion production along a common temperature gradient. As predicted, northern populations produce turions at a lower birth order and earlier; a significant latitude-by-temperature interaction suggests that reaction norm differentiation has occurred. Our results provide evidence of differentiation in reaction norms across latitudes in a phenological trait, and we discuss how the adaptive significance of this plasticity might be further tested.

Fatty acid bioconversion in harpacticoid copepods in a changing environment: a transcriptomic approach

By 2100, global warming is predicted to significantly reduce the capacity of marine primary producers for long-chain polyunsaturated fatty acid (LC-PUFA) synthesis. Primary consumers such as harpacticoid copepods (Crustacea) might mitigate the resulting adverse effects on the food web by increased LC-PUFA bioconversion. Here, we present a high-quality de novo transcriptome assembly of the copepod Platychelipus littoralis, exposed to changes in both temperature (+3°C) and dietary LC-PUFA availability. Using this transcriptome, we detected multiple transcripts putatively coding for LC-PUFA-bioconverting front-end fatty acid (FA) desaturases and elongases, and performed phylogenetic analyses to identify their relationship with sequences of other (crustacean) taxa. While temperature affected the absolute FA concentrations in copepods, LC-PUFA levels remained unaltered even when copepods were fed an LC-PUFA-deficient diet. While this suggests plasticity of LC-PUFA bioconversion within P. littoralis, none of the putative front-end desaturase or elongase transcripts was differentially expressed under the applied treatments. Nevertheless, the transcriptome presented here provides a sound basis for future ecophysiological research on harpacticoid copepods. This article is part of the theme issue 'The next horizons for lipids as 'trophic biomarkers': evidence and significance of consumer modification of dietary fatty acids'.

Effects of temperature changes on life parameters, oxidative stress, and antioxidant defense system in the monogonont marine rotifer Brachionus plicatilis

Global warming is a big concern for all organisms and many efforts have been made to reveal the potential effects of temperature elevation on aquatic organisms. However, limited studies on molecular mechanistic approaches on physiological effects due to temperature changes are available. Here, we investigated the effects of temperature changes on life parameters (e.g., population growth [total number of rotifers], and lifespan), oxidative stress levels and antioxidant activities (e.g., glutathione S-transferase [GST], catalase [CAT], superoxide dismutase [SOD]) with expression levels in the monogonont marine rotifer Brachionus plicatilis. The changes in temperatures led to significant reduction (P < 0.05) in lifespan, possibly due to significant decrease (P < 0.05) in antioxidant activities, reducing the potential to cope with significant elevation in the temperature-induced oxidative stress in B. plicatilis. To further assess the actual induction and clearance of reactive oxygen species (ROS), N-acetyl-L-cysteine was used to examine whether the temperature-induced oxidative stress could be successfully scavenged. Furthermore, expression patterns of the antioxidant-related genes (GSTs, SODs, and CATs) were down- or upregulated (P < 0.05) in response to different temperatures in B. plicatilis. Overall, these findings indicate that ROS-mediated oxidative stress led to cellular damage and antioxidant defense system, resulting in deleterious effects on life parameters in rotifer.

Paternally Inherited P-Element Copy Number Affects the Magnitude of Hybrid Dysgenesis in Drosophila simulans and D. melanogaster

Transposable elements (TEs) are repetitive regions of DNA that are able to self-replicate and reinsert themselves throughout host genomes. Since the discovery of TEs, a prevalent question has been whether increasing TE copy number has an effect on the fitness of their hosts. P-elements (PEs) in Drosophila are a well-studied TE that has strong phenotypic effects. When a female without PEs (M) is crossed to a male with them (P), the resulting females are often sterile, a phenomenon called hybrid dysgenesis (HD). Here, we used short- and long-read sequencing to infer the number of PEs in the genomes of dozens of isofemale lines from two Drosophila species and measured whether the magnitude of HD was correlated with the number of PEs in the paternal genome. Consistent with previous reports, we find evidence for a positive correlation between the paternal PE copy number and the magnitude of HD in progeny from ♀M × ♂ P crosses for both species. Other crosses are not affected by the number of PE copies. We also find that the correlation between the strength of HD and PE copy number differs between species, which suggests that there are genetic differences that might make some genomes more resilient to the potentially deleterious effects of TEs. Our results suggest that PE copy number interacts with other factors in the genome and the environment to cause HD and that the importance of these interactions is species specific.

Mitonuclear mismatch alters performance and reproductive success in naturally introgressed populations of a montane leaf beetle

Coordination between nuclear and mitochondrial genomes is critical to metabolic processes underlying animals' ability to adapt to local environments, yet consequences of mitonuclear interactions have rarely been investigated in populations where individuals with divergent mitochondrial and nuclear genomes naturally interbreed. Genetic variation in the leaf beetle Chrysomela aeneicollis was assessed along a latitudinal thermal gradient in California's Sierra Nevada. Variation at mitochondrial cytochrome oxidase II (COII) and the nuclear gene phosphoglucose isomerase (PGI) shows concordance and was significantly greater along a 65 km transect than 10 other loci. STRUCTURE analyses using neutral loci identified a southern and northern subpopulation, which interbreed in the central drainage Bishop Creek. COII and PGI were used as indicators of mitochondrial and nuclear genetic variation in field and laboratory experiments conducted on beetles from this admixed population. Fecundity, larval development rate, running speed and male mating frequency were higher for beetles with geographically "matched" than "mismatched" mitonuclear genotypes. Effects of mitonuclear mismatch were largest for individuals with northern nuclear genotypes possessing southern mitochondria and were most pronounced after heat treatment or at high elevation. These findings suggest that mitonuclear incompatibility diminishes performance and reproductive success in nature, effects that could intensify at environmental extremes.

The ghost of temperature past: interactive effects of previous and current thermal conditions on gene expression in Manduca sexta

High temperatures can negatively impact the performance and survival of organisms, particularly ectotherms. While an organism's response to high temperature stress clearly depends on current thermal conditions, its response may also be affected by the temporal pattern and duration of past temperature exposures. We used RNA sequencing of Manduca sexta larvae fat body tissue to evaluate how diurnal temperature fluctuations during development affected gene expression both independently and in conjunction with subsequent heat stress. Additionally, we compared gene expression between two M. sexta populations, a lab colony and a genetically related field population that have been separated for >300 generations and differ in their thermal sensitivities. Lab-adapted larvae were predicted to show increased expression responses to both single and repeated thermal stress, whereas recurrent exposure could decrease later stress responses for field individuals. We found large differences in overall gene expression patterns between the two populations across all treatments, as well as population-specific transcriptomic responses to temperature; more differentially expressed genes were upregulated in the field compared with lab larvae. Developmental temperature fluctuations alone had minimal effects on long-term gene expression patterns, with the exception of a somewhat elevated stress response in the lab population. Fluctuating rearing conditions did alter gene expression during exposure to later heat stress, but this effect depended on both the population and the particular temperature conditions. This study contributes to increased knowledge of molecular mechanisms underlying physiological responses of organisms to temperature fluctuations, which is needed for the development of more accurate thermal performance models.

Phenotypic Variation in Growth and Gene Expression Under Different Photoperiods in Allopatric Populations of the Copepod Tigriopus californicus

Daylength is a major environmental condition that varies seasonally and predictably along a latitudinal cline, where higher latitudes exhibit greater ranges in total daylengths. Generally, the circadian clock acts as a network of genes whose expression dynamics are known to control daily rhythms in response to daylength, and it enables the control of many physiological processes such as growth and development. While well studied in many model animals, the influence of daylength variation on phenotypic evolution is poorly examined in marine species. In this study we demonstrate that two allopatric populations of the intertidal crustacean Tigriopus californicus exhibit plastic and divergent phenotypic responses to changes in daylength. Using common-garden experiments, we discovered that shorter daylengths promoted decreased adult body size and faster growth rates in the two divergent populations, suggesting a plastic response to shortened days. In addition, the higher-latitude population exhibited a faster growth rate at any daylength condition, indicating a fixed response, possibly as a result of adaptation to respective natural light regimes. Gene expression profiles of several circadian clock genes, monitored throughout the day by quantitative polymerase chain reaction, revealed that the key core clock genes reach higher daily transcription maxima in the southern population compared to the northern population, pointing to divergent strategies used to respond to changes in daylength. Many modifier genes to the circadian clock showed similar plastic responses to the different daylengths, supporting the existence of at least some conserved gene expression across both populations. Ultimately, our results suggest that photoperiod and daylength exert a potent selective pressure underexplored in marine systems and warranting further future research.

Adaptation to climate change through genetic accommodation and assimilation of plastic phenotypes

Theory suggests that evolutionary changes in phenotypic plasticity could either hinder or facilitate evolutionary rescue in a changing climate. Nevertheless, the actual role of evolving plasticity in the responses of natural populations to climate change remains unresolved. Direct observations of evolutionary change in nature are rare, making it difficult to assess the relative contributions of changes in trait means versus changes in plasticity to climate change responses. To address this gap, this review explores several proxies that can be used to understand evolving plasticity in the context of climate change, including space for time substitutions, experimental evolution and tests for genomic divergence at environmentally responsive loci. Comparisons among populations indicate a prominent role for divergence in environmentally responsive traits in local adaptation to climatic gradients. Moreover, genomic comparisons among such populations have identified pervasive divergence in the regulatory regions of environmentally responsive loci. Taken together, these lines of evidence suggest that divergence in plasticity plays a prominent role in adaptation to climatic gradients over space, indicating that evolving plasticity is also likely to play a key role in adaptive responses to climate change through time. This suggests that genetic variation in plastic responses to the environment (G × E) might be an important predictor of species' vulnerabilities to climate-driven decline or extinction. This article is part of the theme issue 'The role of plasticity in phenotypic adaptation to rapid environmental change'.

Consequences of HSF knockdown on gene expression during the heat shock response in Tigriopus californicus

Although the existence of a cellular heat shock response is nearly universal, its relationship to organismal thermal tolerance is not completely understood. Many of the genes involved are known to be regulated by the highly conserved heat shock transcription factor-1 (HSF-1), yet the regulatory network is not fully characterized. Here, we investigated the role of HSF-1 in gene expression following thermal stress using knockdown of HSF-1 by RNA interference in the intertidal copepod Tigriopus californicus We observed some evidence for decreased transcription of heat shock protein genes following knockdown, supporting the widely acknowledged role of HSF-1 in the heat shock response. However, the majority of differentially expressed genes between the control and HSF-1 knockdown groups were upregulated, suggesting that HSF-1 normally functions to repress their expression. Differential expression observed in genes related to chitin and cuticle formation lends support to previous findings that these processes are highly regulated following heat stress. We performed a genome scan and identified a set of 396 genes associated with canonical heat shock elements. RNA-seq data did not find those genes to be more highly represented in our HSF-1 knockdown treatment, indicating that requirements for binding and interaction of HSF-1 with a given gene are not simply predicted by the presence of HSF-1 binding sites. Further study of the pathways implicated by these results and future comparisons among populations of T. californicus may help us understand the role and importance of HSF-1 in the heat shock response and, more broadly, in organismal thermal tolerance.

Constitutive and Plastic Gene Expression Variation Associated with Desiccation Resistance Differences in the Drosophila americana Species Group

Stress response mechanisms are ubiquitous and important for adaptation to heterogenous environments and could be based on constitutive or plastic responses to environmental stressors. Here we quantify constitutive and plastic gene expression differences under ambient and desiccation stress treatments, in males and females of three species of Drosophila known to differ in desiccation resistance. Drosophila novamexicana survives desiccation trials significantly longer than the two subspecies of Drosophila americana, consistent with its natural species range in the desert southwest USA. We found that desiccation stress reduces global expression differences between species—likely because many general stress response mechanisms are shared among species—but that all species showed plastic expression changes at hundreds of loci during desiccation. Nonetheless, D. novamexicana had the fewest genes with significant plastic expression changes, despite having the highest desiccation resistance. Of the genes that were significantly differentially expressed between species—either within each treatment (>200 loci), constitutively regardless of treatment (36 loci), or with different species-specific plasticity (26 loci)—GO analysis did not find significant enrichment of any major gene pathways or broader functions associated with desiccation stress. Taken together, these data indicate that if gene expression changes contribute to differential desiccation resistance between species, these differences are likely shaped by a relatively small set of influential genes rather than broad genome-wide differentiation in stress response mechanisms. Finally, among the set of genes with the greatest between-species plasticity, we identified an interesting set of immune-response genes with consistent but opposing reaction norms between sexes, whose potential functional role in sex-specific mechanisms of desiccation resistance remains to be determined.

Among-family variation in survival and gene expression uncovers adaptive genetic variation in a threatened fish

Variation in among-family transcriptional responses to different environmental conditions can help to identify adaptive genetic variation, even prior to a selective event. Coupling differential gene expression with formal survival analyses allows for the disentanglement of treatment effects, required for understanding how individuals plastically respond to environmental stressors, from the adaptive genetic variation responsible for differential survival. We combined these two approaches to investigate responses to an emerging conservation issue, thiamine (vitamin B₁ ) deficiency, in a threatened population of Atlantic salmon (Salmo salar). Thiamine is an essential vitamin that is increasingly limited in many ecosystems. In Lake Champlain, Atlantic salmon cannot acquire thiamine in sufficient quantities to support natural reproduction; fertilized eggs must be reared in hatcheries and treated with supplemental thiamine. We evaluated transcriptional responses (via RNA sequencing) to thiamine treatment across families and found 3,616 genes differentially expressed between control (no supplemental thiamine) and treatment individuals. Fewer genes changed expression equally across families (i.e., additively) than exhibited genotype × environment interactions in response to thiamine. Differentially expressed genes were related to known physiological effects of thiamine deficiency, including oxidative stress, cardiovascular irregularities and neurological abnormalities. We also identified 1,446 putatively adaptive genes that were strongly associated with among-family survival in the absence of thiamine treatment, many of which related to neurogenesis and visual perception. Our results highlight the utility of coupling RNA sequencing with formal survival analyses to identify candidate genes that underlie the among-family variation in survival required for an adaptive response to natural selection.

Is geographical variation driving the transcriptomic responses to multiple stressors in the kelp Saccharina latissima?

BACKGROUND

Kelps (Laminariales, Phaeophyceae) are brown macroalgae of utmost ecological, and increasingly economic, importance on temperate to polar rocky shores. Omics approaches in brown algae are still scarce and knowledge of their acclimation mechanisms to the changing conditions experienced in coastal environments can benefit from the application of RNA-sequencing. Despite evidence of ecotypic differentiation, transcriptomic responses from distinct geographical locations have, to our knowledge, never been studied in the sugar kelp Saccharina latissima so far.

RESULTS

In this study we investigated gene expression responses using RNA-sequencing of S. latissima from environments with contrasting temperature and salinity conditions - Roscoff, in temperate eastern Atlantic, and Spitsbergen in the Arctic. Juvenile sporophytes derived from uniparental stock cultures from both locations were pre-cultivated at 8 °C and SA 30. Sporophytes acclimated to 0 °C, 8 °C and 15 °C were exposed to a low salinity treatment (SA 20) for 24 h. Hyposalinity had a greater impact at the transcriptomic level than the temperature alone, and its effects were modulated by temperature. Namely, photosynthesis and pigment synthesis were extensively repressed by low salinity at low temperatures. Although some responses were shared among sporophytes from the different sites, marked differences were revealed by principal component analysis, differential expression and GO enrichment. The interaction between low temperature and low salinity drove the largest changes in gene expression in sporophytes from Roscoff while specimens from Spitsbergen required more metabolic adjustment at higher temperatures. Moreover, genes related to cell wall adjustment were differentially expressed between Spitsbergen and Roscoff control samples.

CONCLUSIONS

Our study reveals interactive effects of temperature and salinity on transcriptomic profiles in S. latissima. Moreover, our data suggest that under identical culture conditions sporophytes from different locations diverge in their transcriptomic responses. This is probably connected to variations in temperature and salinity in their respective environment of origin. The current transcriptomic results support the plastic response pattern in sugar kelp which is a species with several reported ecotypes. Our data provide the baseline for a better understanding of the underlying processes of physiological plasticity and may help in the future to identify strains adapted to specific environments and its genetic control.

Variation in developmental temperature alters adulthood plasticity of thermal tolerance in Tigriopus californicus

In response to environmental change, organisms rely on both genetic adaptation and phenotypic plasticity to adjust key traits that are necessary for survival and reproduction. Given the accelerating rate of climate change, plasticity may be particularly important. For organisms in warming aquatic habitats, upper thermal tolerance is likely to be a key trait, and many organisms express plasticity in this trait in response to developmental or adulthood temperatures. Although plasticity at one life stage may influence plasticity at another life stage, relatively little is known about this possibility for thermal tolerance. Here, we used locally adapted populations of the copepod Tigriopus californicus to investigate these potential effects in an intertidal ectotherm. We found that low latitude populations had greater critical thermal maxima (CT_max) than high latitude populations, and variation in developmental temperature altered CT_max plasticity in adults. After development at 25°C, CT_max was plastic in adults, whereas no adulthood plasticity in this trait was observed after development at 20°C. This pattern was identical across four populations, suggesting that local thermal adaptation has not shaped this effect among these populations. Differences in the capacities to maintain ATP synthesis rates and to induce heat shock proteins at high temperatures, two likely mechanisms of local adaptation in this species, were consistent with changes in CT_max owing to phenotypic plasticity, which suggests that there is likely mechanistic overlap between the effects of plasticity and adaptation. Together, these results indicate that developmental effects may have substantial impacts on upper thermal tolerance plasticity in adult ectotherms.

Transcriptomics of Tasmanian Devil (Sarcophilus Harrisii) Ear Tissue Reveals Homogeneous Gene Expression Patterns across a Heterogeneous Landscape

In an era of unprecedented global change, exploring patterns of gene expression among wild populations across their geographic range is crucial for characterizing adaptive potential. RNA-sequencing studies have successfully characterized gene expression differences among populations experiencing divergent environmental conditions in a wide variety of taxa. However, few of these studies have identified transcriptomic signatures to multivariate, environmental stimuli among populations in their natural environments. Herein, we aim to identify environmental and sex-driven patterns of gene expression in the Tasmanian devil (Sarcophilus harrisii), a critically endangered species that occupies a heterogeneous environment. We performed RNA-sequencing on ear tissue biopsies from adult male and female devils from three populations at the extremes of their geographic range. There were no transcriptome-wide patterns of differential gene expression that would be suggestive of significant, environmentally-driven transcriptomic responses. The general lack of transcriptome-wide variation in gene expression levels across the devil's geographic range is consistent with previous studies that documented low levels of genetic variation in the species. However, genes previously implicated in local adaptation to abiotic environment in devils were enriched for differentially expressed genes. Additionally, three modules of co-expressed genes were significantly associated with either population of origin or sex.

Within species expressed genetic variability and gene expression response to different temperatures in the rotifer Brachionus calyciflorus sensu stricto

Genetic divergence is impacted by many factors, including phylogenetic history, gene flow, genetic drift, and divergent selection. Rotifers are an important component of aquatic ecosystems, and genetic variation is essential to their ongoing adaptive diversification and local adaptation. In addition to coding sequence divergence, variation in gene expression may relate to variable heat tolerance, and can impose ecological barriers within species. Temperature plays a significant role in aquatic ecosystems by affecting species abundance, spatio-temporal distribution, and habitat colonization. Recently described (formerly cryptic) species of the Brachionus calyciflorus complex exhibit different temperature tolerance both in natural and in laboratory studies, and show that B. calyciflorus sensu stricto (s.s.) is a thermotolerant species. Even within B. calyciflorus s.s., there is a tendency for further temperature specializations. Comparison of expressed genes allows us to assess the impact of stressors on both expression and sequence divergence among disparate populations within a single species. Here, we have used RNA-seq to explore expressed genetic diversity in B. calyciflorus s.s. in two mitochondrial DNA lineages with different phylogenetic histories and differences in thermotolerance. We identify a suite of candidate genes that may underlie local adaptation, with a particular focus on the response to sustained high or low temperatures. We do not find adaptive divergence in established candidate genes for thermal adaptation. Rather, we detect divergent selection among our two lineages in genes related to metabolism (lipid metabolism, metabolism of xenobiotics).

Regional heterogeneity impacts gene expression in the subarctic zooplankter Neocalanus flemingeri in the northern Gulf of Alaska

Marine pelagic species are being increasingly challenged by environmental change. Their ability to persist will depend on their capacity for physiological acclimatization. Little is known about limits of physiological plasticity in key species at the base of the food web. Here we investigate the capacity for acclimatization in the copepod Neocalanus flemingeri, which inhabits the Gulf of Alaska, a heterogeneous and highly seasonal environment. RNA-Seq analysis of field-collected pre-adults identified large regional differences in expression of genes involved in metabolic and developmental processes and response to stressors. We found that lipid synthesis genes were up-regulated in individuals from Prince William Sound and down-regulated in the Gulf of Alaska. Up-regulation of lipid catabolic genes in offshore individuals suggests they are experiencing nutritional deficits. The expression differences demonstrate physiological plasticity in response to a steep gradient in food availability. Our transcriptional analysis reveals mechanisms of acclimatization that likely contribute to the observed resilience of this population.

Little evidence of adaptation potential to ocean acidification in sea urchins living in "Future Ocean" conditions at a CO2 vent

Ocean acidification (OA) can be detrimental to calcifying marine organisms, with stunting of invertebrate larval development one of the most consistent responses. Effects are usually measured by short-term, within-generation exposure, an approach that does not consider the potential for adaptation. We examined the genetic response to OA of larvae of the tropical sea urchin Echinometra sp. C. raised on coral reefs that were either influenced by CO2 vents (pH ~ 7.9, future OA condition) or nonvent control reefs (pH 8.2). We assembled a high quality de novo transcriptome of Echinometra embryos (8 hr) and pluteus larvae (48 hr) and identified 68,056 SNPs. We tested for outlier SNPs and functional enrichment in embryos and larvae raised from adults from the control or vent sites. Generally, highest F ST values in embryos were observed between sites (intrinsic adaptation, most representative of the gene pool in the spawned populations). This comparison also had the highest number of outlier loci (40). In the other comparisons, classical adaptation (comparing larvae with adults from the control transplanted to either the control or vent conditions) and reverse adaptation (larvae from the vent site returned to the vent or explanted at the control), we only observed modest numbers of outlier SNPs (6-19) and only enrichment in two functional pathways. Most of the outliers detected were silent substitutions without adaptive potential. We conclude that there is little evidence of realized adaptation potential during early development, while some potential (albeit relatively low) exists in the intrinsic gene pool after more than one generation of exposure.

Transgenerational Proteome Plasticity in Resilience of a Marine Copepod in Response to Environmentally Relevant Concentrations of Microplastics

Here, we examined the multigenerational effect of microplastics (6-μm polystyrene beads; with different environmentally relevant concentrations of 0.023 and 0.23 mg/L in seawater) on the marine copepod Tigriopus japonicus under two-generation exposure (F0-F1) followed by one-generation recovery (F2) in clean seawater. Also, the seven life-history traits (survival, sex ratio, developmental time of nauplius phase, developmental time to maturation, number of clutches, number of nauplii/clutch, and fecundity) were measured for each generation. Furthermore, to investigate within-generation proteomic response and transgenerational proteome plasticity, proteome profiling was conducted for the F1 and F2 copepods under the control and 0.23 mg/L microplastics treatment. The results showed successful ingestion of microplastics in F0-F1 under both exposure concentrations, while higher concentration (0.23 mg/L) of microplastics resulted in the significant reduction in survival rate, number of nauplii/clutch, and fecundity. However, the affected traits were totally restored in the recovery generation (F2). Proteomic analysis demonstrated that microplastics exposure increased several cellular biosynthesis processes and, in turn, reduced energy storage due to the trade-off, hence compromising survival and reproduction of the treated copepods in F1. Interestingly, the two-generational effect of microplastics in copepods had significant transgenerational proteome plasticity as demonstrated by increased energy metabolism and stress-related defense pathway, which accounts for regaining of the compromised phenotypic traits during recovery (i.e., F2). Overall, this study provides a molecular understanding on the effect of microplastics at a translational level under long-term multigenerational exposure in marine copepods, and also the transgenerational proteome plasticity is likely rendering the robustness of copepods in response to microplastics pollution.

Effects of oxidative stress on sex-specific gene expression in the copepod Tigriopus californicus revealed by single individual RNA-seq

Oxidative stress reflects the imbalance of pro-oxidants and antioxidants. Prolonged oxidative stress can induce cellular damage, diseases and aging, and the effects may be sex-specific. Tigriopus californicus has recently been proposed as an alternative model system for sex-specific studies due to the absence of sex chromosomes. In this study, we used comparative transcriptomic analyses to assess sex-specific transcriptional responses to oxidative stress. Male and female individuals were maintained separately in one of three treatments: 1) control conditions with an algae diet, 2) pro-oxidant (H₂O₂) conditions with an algae diet or 3) decreased antioxidant conditions (reduced carotenoids due to a yeast diet). Single individual RNA-seq was then conducted for twenty-four libraries using Ligation Mediated RNA sequencing (LM-Seq). Variance in gene expression was partitioned into 62.3% between sexes, 26.85% among individuals and 10.85% among treatments. Within each of the three treatments, expression was biased toward females. However, compared to the control treatment, males in both pro-oxidant and decreased antioxidant treatments differentially expressed more genes while females differentially expressed fewer genes but with a greater magnitude of fold change. As the first study of copepods to apply single individual RNA-seq, the findings will contribute to a better understanding of transcriptomic variation among individuals as well as sex-specific response mechanisms to oxidative stress in the absence of sex chromosomes.

Multiple stressors produce differential transcriptomic patterns in a stream-dwelling salamander

BACKGROUND

Global biodiversity is decreasing at an alarming rate and amphibians are at the forefront of this crisis. Understanding the factors that negatively impact amphibian populations and effectively monitoring their health are fundamental to addressing this epidemic. Plasma glucocorticoids are often used to assess stress in amphibians and other vertebrates, but these hormones can be extremely dynamic and impractical to quantify in small organisms. Transcriptomic responses to stress hormones in amphibians have been largely limited to laboratory models, and there have been few studies on vertebrates that have evaluated the impact of multiple stressors on patterns of gene expression. Here we examined the gene expression patterns in tail tissues of stream-dwelling salamanders (Eurycea tynerensis) chronically exposed to the stress hormone corticosterone under different temperature regimes.

RESULTS

We found unique transcriptional signatures for chronic corticosterone exposure that were independent of temperature variation. Several of the corticosterone responsive genes are known to be involved in immune system response (LY-6E), oxidative stress (GSTM2 and TRX), and tissue repair (A2M and FX). We also found many genes to be influenced by temperature (CIRBP, HSC71, HSP40, HSP90, HSP70, ZNF593). Furthermore, the expression patterns of some genes (GSTM2, LY-6E, UMOD, ZNF593, CIRBP, HSP90) show interactive effects of temperature and corticosterone exposure, compared to each treatment alone. Through a series of experiments we also showed that stressor induced patterns of expression were largely consistent across ages, life cycle modes, and tissue regeneration.

CONCLUSIONS

Outside of thermal stressors, the application of transcriptomes to monitor the health of non-human vertebrate systems has been vastly underinvestigated. Our study suggests that transcriptomic patterns harbor stressor specific signatures that can be highly informative for monitoring the diverse stressors of amphibian populations.