Expression-environment associations in transcriptomic heat stress responses for a global plant lineage

The increasing frequency and severity of heatwaves will intensify stress on plants. Given regional variation in heatwave exposure and expected differences in thermal tolerance between species it is unlikely that all plant species will be affected equally by climate change. However, little is currently known about variation in the responses of plants to heat stress, or how those responses differ among closely related species adapted to different environments. Here we quantify the response of 17 Acacia species (175 RNA-seq libraries), from across Australia's diverse biomes, to a multi-day experimental heatwave treatment to identify variation in transcriptomic and physiological responses to heat stress. Genes with known heat response functions showed consistent responses across Acacia species. Up to 10% of all genes and over 100 gene families showed significant clinal variation in the magnitude of their expression plasticity across species. Specifically, gene families linked to the temperature stress response were overrepresented among significant relationships with home range temperature conditions. Gene expression responses seen on the first day of the heatwave were more frequently associated with home range climates, while expression responses by day four were more commonly related to photosystem II acclimation. Comparative transcriptomics on non-model species has the potential to provide key information on stress response plasticity, especially when linked with our understanding of model species. Our study indicates that the pressing challenge to identifying potentially vulnerable species to climate change could be benefited by the further exploration of clinal variation in transcriptome plasticity.

Patterns of methylation and transcriptional plasticity during thermal acclimation in a reef-building coral

Phenotypic plasticity can buffer organisms against short-term environmental fluctuations. For example, previous exposure to increased temperatures can increase thermal tolerance in many species. Prior studies have found that acclimation to higher temperature can influence the magnitude of transcriptional response to subsequent acute thermal stress (hereafter, "transcriptional response modulation"). However, mechanisms mediating this gene expression response and, ultimately, phenotypic plasticity remain largely unknown. Epigenetic modifications are good candidates for modulating transcriptional response, as they broadly correlate with gene expression. Here, we investigate changes in DNA methylation as a possible mechanism controlling shifts in gene expression plasticity and thermal acclimation in the reef-building coral Acropora nana. We find that gene expression response to acute stress is altered in corals acclimated to different temperatures, with many genes exhibiting a dampened response to heat stress in corals pre-conditioned to higher temperatures. At the same time, we observe shifts in methylation during both acclimation (11 days) and acute heat stress (24 h). We observed that the acute heat stress results in shifts in gene-level methylation and elicits an acute transcriptional response in distinct gene sets. Further, acclimation-induced shifts in gene expression plasticity and differential methylation also largely occur in separate sets of genes. Counter to our initial hypothesis no overall correlation between the magnitude of differential methylation and the change in gene expression plasticity. We do find a small but statistically significant overlap in genes exhibiting both dampened expression response and shifts in methylation (14 genes), which could be candidates for further inquiry. Overall, our results suggest transcriptional response modulation occurs independently from methylation changes induced by thermal acclimation.

Environmental factors influence cross-talk between a heat shock protein and an oxidative stress protein modification in the lizard Gallotia galloti

Better understanding how organisms respond to their abiotic environment, especially at the biochemical level, is critical in predicting population trajectories under climate change. In this study, we measured constitutive stress biomarkers and protein post-translational modifications associated with oxidative stress in Gallotia galloti, an insular lizard species inhabiting highly heterogeneous environments on Tenerife. Tenerife is a small volcanic island in a relatively isolated archipelago off the West coast of Africa. We found that expression of GRP94, a molecular chaperone protein, and levels of protein carbonylation, a marker of cellular stress, change across different environments, depending on solar radiation-related variables and topology. Here, we report in a wild animal population, cross-talk between the baseline levels of the heat shock protein-like GRP94 and oxidative damage (protein carbonylation), which are influenced by a range of available temperatures, quantified through modelled operative temperature. This suggests a dynamic trade-off between cellular homeostasis and oxidative damage in lizards adapted to this thermally and topologically heterogeneous environment.

Genomics for monitoring and understanding species responses to global climate change

All life forms across the globe are experiencing drastic changes in environmental conditions as a result of global climate change. These environmental changes are happening rapidly, incur substantial socioeconomic costs, pose threats to biodiversity and diminish a species' potential to adapt to future environments. Understanding and monitoring how organisms respond to human-driven climate change is therefore a major priority for the conservation of biodiversity in a rapidly changing environment. Recent developments in genomic, transcriptomic and epigenomic technologies are enabling unprecedented insights into the evolutionary processes and molecular bases of adaptation. This Review summarizes methods that apply and integrate omics tools to experimentally investigate, monitor and predict how species and communities in the wild cope with global climate change, which is by genetically adapting to new environmental conditions, through range shifts or through phenotypic plasticity. We identify advantages and limitations of each method and discuss future research avenues that would improve our understanding of species' evolutionary responses to global climate change, highlighting the need for holistic, multi-omics approaches to ecosystem monitoring during global climate change.

Chronic exposure to environmental temperature attenuates the thermal sensitivity of salmonids

Metabolism, the biological processing of energy and materials, scales predictably with temperature and body size. Temperature effects on metabolism are normally studied via acute exposures, which overlooks the capacity for organisms to moderate their metabolism following chronic exposure to warming. Here, we conduct respirometry assays in situ and after transplanting salmonid fish among different streams to disentangle the effects of chronic and acute thermal exposure. We find a clear temperature dependence of metabolism for the transplants, but not the in-situ assays, indicating that chronic exposure to warming can attenuate salmonid thermal sensitivity. A bioenergetic model accurately captures the presence of fish in warmer streams when accounting for chronic exposure, whereas it incorrectly predicts their local extinction with warming when incorporating the acute temperature dependence of metabolism. This highlights the need to incorporate the potential for thermal acclimation or adaptation when forecasting the consequences of global warming on ecosystems.

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.

Genetic diversity and gene expression diversity shape the adaptive pattern of the aquatic plant Batrachium bungei along an altitudinal gradient on the Qinghai-Tibet plateau

It is an intriguing issue of evolutionary biology how genetic diversity and gene expression diversity shape the adaptive patterns. Comparative transcriptomic studies of wild populations in extreme environments provide critical insights into the relative contribution of genetic and expressive components. In this study, we analyzed the genetic diversity and gene expression diversity of 20 populations of the aquatic plant Batrachium bungei along elevations ranging from 2690 to 4896 m on the Qinghai-Tibet plateau (QTP). Based on single nucleotide polymorphisms (SNPs) and gene expression data from 100 individuals of B. bungei, we found that variation in genetic sequence was more sensitive to detect weak differentiation than gene expression. Using 292,613 high-quality SNPs, we documented a significant phylogeographical structure, a low within-population genetic diversity, and a high inter-population genetic differentiation in B. bungei populations. Analysis of relationship between geographic distance, genetic distance, and gene expression similarity showed that geographic isolation shaped gene flow patterns but not gene expression patterns. We observed a negative relationship between genetic diversity and gene expression diversity within and among B. bungei populations, and we demonstrated that as environmental conditions worsen with increasing altitude, genetic diversity played an increased role in maintaining the stability of populations, while the corresponding role of gene expression diversity decreased. These results suggested that genetic diversity and gene expression diversity might act as a complementary mechanism contributing to the long-term survival of B. bungei in extreme environments.

Transcriptome analysis reveals the diverse response of pearl oyster Pinctada fucata martensii after different PAMP stimulation

Bivalves have evolved effective strategies to combat different pathogens in the environment. They rely on innate immunity to deal with the invasion of various bacteria, viruses, and other microorganisms. However, the molecular mechanisms underlying the responses remain largely unknown. Herein, we constructed 21 transcriptomes of the hemocytes after lipopolysaccharide (LPS), peptidoglycan (PGN) and polyinosinic-polycytidylic acid (poly(I:C)) stimulation to investigate the molecular mechanisms underlying adaptations and plastic responses to different pathogen-related molecular patterns (PAMPs) in pearl oyster Pinctada fucata martensii. Transcriptome analysis revealed 1986-3427 responsive genes enriched in the major immune and cell cycle-related pathways at different times after PAMP stimulation, and the expression patterns of genes under these pathways are complex and diverse. Moreover, "lysosomes" were enriched 6 h after LPS and PGN stimulation, while "peroxisomes" were only enriched in poly(I:C) group. These results suggest different response strategies of pearl oyster to different PAMPs. Furthermore, we identified 261 pattern-recognition receptors (PRRs) including 4 retinoic acid-inducible gene I-like receptors, 38 NOD-like receptors, 83 Toll-like receptors, and 136 C-type lectins in the genome of P. f. martensii. The diverse expression patterns of these PRRs after different PAMP stimulation indicated that pearl oyster evolved complex and specific recognition systems due to tandem repeat and diverse domain combination, which may help pearl oyster cope with the different pathogens in the environment. The present study improved our understanding of the molecular response of pearl oyster to different PAMP stimulation.

Effects of ocean acidification and tralopyril on bivalve biomineralization and carbon cycling: A study of the Pacific Oyster (Crassostrea gigas)

The combined effects of emerging pollutants and ocean acidification (OA) on marine organisms and marine ecosystems have attracted increasing attention. However, the combined effects of tralopyril and OA on marine organisms and marine ecosystems remain unclear. In this study, Crassostrea gigas (C. gigas) were exposed to tralopyril (1 μg/L) and/or OA (PH = 7.7) for 21 days and a 14-day recovery acclimation. To investigate the stress response and potential molecular mechanisms of C. gigas to OA and tralopyril exposure alone or in combination, as well as the effects of OA and/or tralopyril on bivalve biomineralization and marine carbon cycling. The results showed that the combined toxicity was between that of acidification and tralopyril alone. Single or combined exposure activated the general stress defense responses of C. gigas mantle, affected energy metabolism and biomineralization of the organism and the carbon cycle of the marine ecosystem. Moreover, acidification-induced and tralopyril-induced toxicity showed potential recoverability at molecular and biochemical levels. This study provides a new perspective on the molecular mechanisms of tralopyril toxicity to bivalve shellfish and reveals the potential role of tralopyril and OA on marine carbon cycling.

Transcriptomic stability or lability explains sensitivity to climate stressors in coralline algae

Background

Crustose coralline algae (CCA) are calcifying red macroalgae that play important ecological roles including stabilisation of reef frameworks and provision of settlement cues for a range of marine invertebrates. Previous research into the responses of CCA to ocean warming (OW) and ocean acidification (OA) have found magnitude of effect to be species-specific. Response to OW and OA could be linked to divergent underlying molecular processes across species.

Results

Here we show Sporolithon durum, a species that exhibits low sensitivity to climate stressors, had little change in metabolic performance and did not significantly alter the expression of any genes when exposed to temperature and pH perturbations. In contrast, Porolithon onkodes, a major coral reef builder, reduced photosynthetic rates and had a labile transcriptomic response with over 400 significantly differentially expressed genes, with differential regulation of genes relating to physiological processes such as carbon acquisition and metabolism. The differential gene expression detected in P. onkodes implicates possible key metabolic pathways, including the pentose phosphate pathway, in the stress response of this species.

Conclusions

We suggest S. durum is more resistant to OW and OA than P. onkodes, which demonstrated a high sensitivity to climate stressors and may have limited ability for acclimatisation. Understanding changes in gene expression in relation to physiological processes of CCA could help us understand and predict how different species will respond to, and persist in, future ocean conditions predicted for 2100.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08931-9.

Differential transcriptomic responses to heat stress in surface and subterranean diving beetles

Subterranean habitats are generally very stable environments, and as such evolutionary transitions of organisms from surface to subterranean lifestyles may cause considerable shifts in physiology, particularly with respect to thermal tolerance. In this study we compared responses to heat shock at the molecular level in a geographically widespread, surface-dwelling water beetle to a congeneric subterranean species restricted to a single aquifer (Dytiscidae: Hydroporinae). The obligate subterranean beetle Paroster macrosturtensis is known to have a lower thermal tolerance compared to surface lineages (CT_max 38 °C cf. 42–46 °C), but the genetic basis of this physiological difference has not been characterized. We experimentally manipulated the thermal environment of 24 individuals to demonstrate that both species can mount a heat shock response at high temperatures (35 °C), as determined by comparative transcriptomics. However, genes involved in these responses differ between species and a far greater number were differentially expressed in the surface taxon, suggesting it can mount a more robust heat shock response; these data may underpin its higher thermal tolerance compared to subterranean relatives. In contrast, the subterranean species examined not only differentially expressed fewer genes in response to increasing temperatures, but also in the presence of the experimental setup employed here alone. Our results suggest P. macrosturtensis may be comparatively poorly equipped to respond to both thermally induced stress and environmental disturbances more broadly. The molecular findings presented here have conservation implications for P. macrosturtensis and contribute to a growing narrative concerning weakened thermal tolerances in obligate subterranean organisms at the molecular level.

Molecular responses of a key Antarctic species to sedimentation due to rapid climate change

Rapid regional warming causing glacial retreat and melting of ice caps in Antarctica leads benthic filter-feeders to be exposed to periods of food shortage and high respiratory impairment as a consequence of seasonal sediment discharge in the West Antarctic Peninsula coastal areas. The molecular physiological response and its fine-tuning allow species to survive acute environmental stress and are thus a prerequisite to longer-term adaptation to changing environments. Under experimental conditions, we analyzed here the metabolic response to changes in suspended sediment concentrations, through transcriptome sequencing and enzymatic measurements in a highly abundant Antarctic ascidian. We found that the mechanisms underlying short-term response to sedimentation in Cnemidocarpa verrucosa sp. A involved apoptosis, immune defense, and general metabolic depression. These mechanisms may be understood as an adaptive protection against sedimentation caused by glacial retreat. This process can strongly contribute to the structuring of future benthic filter-feeder communities in the face of climate change.

Demographic, physiological, and genetic factors linked to the poleward range expansion of the snail Nerita yoldii along the shoreline of China

Species range shift is one of the most significant consequences of climate change in the Anthropocene. A comprehensive study, including demographic, physiological, and genetic factors linked to poleward range expansion, is crucial for understanding how the expanding population occupies the new habitat. In the present study, we investigated the demographic, physiological, and genetic features of the intertidal gastropod Nerita yoldii, which has extended its northern limit by ~200 km over the former biogeographic break of the Yangtze River Estuary during recent decades. The neutral SNPs data showed that the new marginal populations formed a distinct cluster established by a few founders. Demographic modelling analysis revealed that the new marginal populations experienced a strong genetic bottleneck followed by recent demographic expansion. Successful expansion that overcame the founder effect might be attributed to its high capacity of rapid population growth and multiple introductions. According to the non-neutral SNPs under diversifying selection, there were high levels of heterozygosity in the new marginal populations, which might be beneficial for adapting to the novel thermal conditions. The common garden experiment showed that the new marginal populations have evolved divergent transcriptomic and physiological responses to heat stress, allowing them to occupy and survive in the novel environment. Lower transcriptional plasticity was observed in the new marginal populations. These results suggest a new biogeographic pattern of N. yoldii has formed with the occurrence of demographic, physiologic, and genetic changes, and emphasize the roles of adaptation of marginal populations during range expansion.

Gene expression studies of plastic and evolutionary responses to global warming

Phenotypic plasticity can be a rapid response for coping with global warming, yet may be insufficient to protect species from extinction. Evolutionary adaptation may reinforce adaptive or oppose maladaptive plastic responses. With advances in technology whole transcriptomes can provide us with an unprecedented overview of genes and functional processes underlying the interplay between plasticity and evolution. We advocate that insects provide ideal opportunities to study plasticity in non-adapted and thermally adapted populations to infer reaction norms across the whole transcriptome ('reactionomes'). This can advance our understanding of how the interplay between plasticity and evolution shapes responses to warming. So far, a limited number of studies suggest predominantly maladaptive plastic responses to novel environments that are reduced with time, but much more research is needed to infer general patterns.

Comparative Transcriptomics of the Northern Quahog Mercenaria mercenaria and Southern Quahog Mercenaria campechiensis in Response to Chronic Heat Stress

The northern quahog (Mercenaria mercenaria) supports lucrative aquaculture industries in the USA. In the southeastern USA, aquacultured M. mercenaria faces increasing risks of summer die-offs from prolonged heat waves. We used a comparative transcriptomic approach to investigate the molecular responses of M. mercenaria and its southern congener, Mercenaria campechiensis, to controlled incremental heat stress over a 4-week period. Mercenaria were exposed to temperatures from 24 to 34 °C with 2.5 °C/week, after which, gill transcriptomes were de novo assembled and annotated. During the 4 weeks of chronic heat exposure, both species had the same survival rate (96%); M. mercenaria experienced body weight gain/loss depending on the originated hatcheries while M. campechiensis experienced an average net weight loss. The upregulated genes in both species included those in chaperone-mediated protein folding and regulation of cell death pathways, while the downregulated genes in both species involved in mRNA processing and splicing pathways. Compared to M. mercenaria, M. campechiensis appears to be more sensitive to prolonged heat stress as indicated by upregulating significantly more genes in coping with oxidative stress and in the protein degradation pathways, while downregulating some inhibitors of apoptosis. We discussed this finding within their ecological and evolutionary context. Our findings highlighted the potential vulnerability of the two quahogs, especially the southern quahog, to continued ocean warming.

Insights into molecular pathways and fatty acid membrane composition during the temperature stress response in the murine C2C12 cell model

Daily and seasonal temperature fluctuations are inevitable due to climate change, which highlights the importance of studying the detrimental effects of temperature fluctuations on the health, productivity, and product quality of farm animals. Muscle membrane composition and the molecular signals are vital for muscle cell differentiation and muscle growth, but their response to temperature stress is not well characterized. Temperature changes can lead to modification of membrane components of the cell, which may affect its surroundings and intracellular signaling pathways. Using C2C12 myoblast cells as a model of skeletal muscle development, this study was designed to investigate the effects of high temperature (39 °C and 41 °C) and low temperature (35 °C) on molecular pathways in the cells as well as the cell membrane fatty acid composition. Our results show that several genes were differentially expressed in C2C12 cells cultured under heat or cold stress, and these genes were enriched important KEGG pathways including PI3K-Akt signaling pathway, lysosome and HIF- signaling pathway, Wnt signaling pathway and AMPK signaling pathway. Our analysis further reveals that several membrane transporters and genes involved in lipid metabolism and fatty acid elongation were also differentially expressed in C2C12 cells cultured under high or low temperature. Additionally, temperature stress shifts the fatty acid composition in the cell membranes, including the proportion of saturated, monounsaturated and polyunsaturated fatty acids. This study revealed an interference between fatty acid composition in the membranes and changing molecular pathways including lipid metabolism and fatty acids elongation mediated under thermal stress. These findings will reinforce a better understanding of the adaptive mechanisms in skeletal muscle under temperature stress.

Inferring the genetic responses to acute drought stress across an ecological gradient

BACKGROUND

How do xerophytic species thrive in environments that experience extreme annual drought? Although critical to the survival of many species, the genetic responses to drought stress in many non-model organisms has yet to be explored. We investigated this question in Mentzelia section Bartonia (Loasaceae), which occurs throughout western North America, including arid lands. To better understand the genetic responses to drought stress among species that occur in different habitats, the gene expression levels of three species from Mentzelia were compared across a precipitation gradient. Two de novo reference transcriptomes were generated and annotated. Leaf and root tissues were collected from control and drought shocked plants and compared to one another for differential expression. A target-gene approach was also implemented to better understand how drought-related genes from model and crop species function in non-model systems.

RESULTS

When comparing the drought-shock treatment plants to their respective control plants, we identified 165 differentially expressed clusters across all three species. Differentially expressed genes including those associated with water movement, photosynthesis, and delayed senescence. The transcriptome profiling approach was coupled with a target genes approach that measured expression of 90 genes associated with drought tolerance in model organisms. Comparing differentially expressed genes with a ≥ 2 log-fold value between species and tissue types showed significant differences in drought response. In pairwise comparisons, species that occurred in drier environments differentially expressed greater genes in leaves when drought shocked than those from wetter environments, but expression in the roots mostly produced opposite results.

CONCLUSIONS

Arid-adapted species mount greater genetic responses compared to the mesophytic species, which has likely evolved in response to consistent annual drought exposure across generations. Drought responses also depended on organ type. Xerophytes, for example, mounted a larger response in leaves to downregulate photosynthesis and senescence, while mobilizing carbon and regulating water in the roots. The complexity of drought responses in Mentzelia suggest that whole organism responses need to be considered when studying drought and, in particular, the physiological mechanisms in which plants regulate water, carbon, cell death, metabolism, and secondary metabolites.

Differential DNA methylation across environments has no effect on gene expression in the eastern oyster

It has been hypothesized that environmentally induced changes to gene body methylation could facilitate adaptive transgenerational responses to changing environments. We compared patterns of global gene expression (Tag-seq) and gene body methylation (reduced representation bisulfite sequencing) in 80 eastern oysters Crassostrea virginica from six full-sib families, common gardened for 14 months at two sites in the northern Gulf of Mexico that differed in mean salinity. At the time of sampling, oysters from the two sites differed in mass by 60% and in parasite loads by nearly two orders of magnitude. They also differentially expressed 35% of measured transcripts. However, we observed differential methylation at only 1.4% of potentially methylated loci in comparisons between individuals from these different environments, and little correspondence between differential methylation and differential gene expression. Instead, methylation patterns were largely driven by genetic differences among families, with a PERMANOVA analysis indicating nearly a two orders of magnitude greater number of genes differentially methylated between families than between environments. An analysis of CpG observed/expected values (CpG O/E) across the C. virginica genome showed a distinct bimodal distribution, with genes from the first cluster showing the lower CpG O/E values, greater methylation and higher and more stable gene expression, while genes from the second cluster showed lower methylation, and lower and more variable gene expression. Taken together, the differential methylation results suggest that only a small portion of the C. virginica genome is affected by environmentally induced changes in methylation. At this point, there is little evidence to suggest that environmentally induced methylation states would play a leading role in regulating gene expression responses to new environments.

Lack of genotype-by-environment interaction suggests limited potential for evolutionary changes in plasticity in the eastern oyster, Crassostrea virginica

Eastern oysters in the northern Gulf of Mexico are facing rapid environmental changes and can respond to this change via plasticity or evolution. Plasticity can act as an immediate buffer against environmental change, but this buffering could impact the organism's ability to evolve in subsequent generations. While plasticity and evolution are not mutually exclusive, the relative contribution and interaction between them remains unclear. In this study, we investigate the roles of plastic and evolved responses to environmental variation and Perkinsus marinus infection in Crassostrea virginica by using a common garden experiment with 80 oysters from six families outplanted at two field sites naturally differing in salinity. We use growth data, P. marinus infection intensities, 3' RNA sequencing (TagSeq) and low-coverage whole-genome sequencing to identify the effect of genotype, environment and genotype-by-environment interaction on the oyster's response to site. As one of first studies to characterize the joint effects of genotype and environment on transcriptomic and morphological profiles in a natural setting, we demonstrate that C. virginica has a highly plastic response to environment and that this response is parallel among genotypes. We also find that genes responding to genotype have distinct and opposing profiles compared to genes responding to environment with regard to expression levels, Ka/Ks ratios and nucleotide diversity. Our findings suggest that C. virginica may be able to buffer the immediate impacts of future environmental changes by altering gene expression and physiology, but the lack of genetic variation in plasticity suggests limited capacity for evolved responses.

Mechanism underlying the toxicity of the microplastic fibre transfer in the sea cucumber Apostichopus japonicus

Microscopic plastic particles (0.1 µm-5 mm) are widespread hazardous pollutants, and microfibres (MFs) are their dominant shape in habitats. Previous field and laboratory studies have demonstrated that MFs enter the coelomic fluid of sea cucumbers from the water through the respiratory tree. However, the possible mechanism underlying the toxicity of this process is not well understood. Herein, RNA-Seq was used to examine the responses of the respiratory tree during the MF transfer process in the sea cucumber Apostichopus japonicus. Polyester synthetic MFs were used, and the number of transferred MFs was controlled to the amount reported from the field. The results showed that the MFs altered gene expression as the transfer process increased. The top genes regulated by MF transfer were mainly related to metabolic processes and signal transduction pathways, with upregulated genes following low MF transfer and downregulated genes following high MF transfer. Functional enrichment analysis revealed the pathways in which differentially expressed genes were enriched under different MF transfer scenarios. The transcriptomic findings were further supported by histological observations, which revealed injury and loss of cell components. This study contributes to understanding the effects of MFs in a valuable echinoderm species through transcriptomic and histological examinations.

The use of non-lethal sampling for transcriptomics to assess the physiological status of wild fishes

Fishes respond to different abiotic and biotic stressors through changes in gene expression as a part of an integrated physiological response. Transcriptomics approaches have been used to quantify gene expression patterns as a reductionist approach to understand responses to environmental stressors in animal physiology and have become more commonly used to study wild fishes. We argue that non-lethal sampling for transcriptomics should become the norm for assessing the physiological status of wild fishes, especially when there are conservation implications. Processes at the level of the transcriptome provide a "snapshot" of the cellular conditions at a given time; however, by using a non-lethal sampling protocol, researchers can connect the transcriptome profile with fitness-relevant ecological endpoints such as reproduction, movement patterns and survival. Furthermore, telemetry is a widely used approach in fisheries to understand movement patterns in the wild, and when combined with transcriptional profiling, provides arguably the most powerful use of non-lethal sampling for transcriptomics in wild fishes. In this review, we discuss the different tissues that can be successfully incorporated into non-lethal sampling strategies, which is particularly useful in the context of the emerging field of conservation transcriptomics. We briefly describe different methods for transcriptional profiling in fishes from high-throughput qPCR to whole transcriptome approaches. Further, we discuss strategies and the limitations of using transcriptomics for non-lethally studying fishes. Lastly, as 'omics' technology continues to advance, transcriptomics paired with different omics approaches to study wild fishes will provide insight into the factors that regulate phenotypic variation and the physiological responses to changing environmental conditions in the future.

Adaptive Differences in Gene Expression in Farm-Impacted Seedbeds of the Native Blue Mussel Mytilus chilensis

The study of adaptive population differences is relevant for evolutionary biology, as it evidences the power of selective local forces relative to gene flow in maintaining adaptive phenotypes and their underlying genetic determinants. However, human-mediated hybridization through habitat translocations, a common and recurrent aquaculture practice where hybrids could eventually replace local genotypes, risk populations' ability to cope with perturbations. The endemic marine mussel Mytilus chilensis supports a booming farming industry in the inner sea of Chiloé Island, southern Chile, which entirely relies on artificially collected seeds from natural beds that are translocated to ecologically different fattening centers. A matter of concern is how farm-impacted seedbeds will potentially cope with environmental shifts and anthropogenic perturbations. This study provides the first de novo transcriptome of M. chilensis; assembled from tissue samples of mantles and gills of individuals collected in ecologically different farm-impacted seedbeds, Cochamó (41°S) and Yaldad (43°S). Both locations and tissue samples differentially expressed transcripts (DETs) in candidate adaptive genes controlling multiple fitness traits, involved with metabolism, genetic and environmental information processing, and cellular processes. From 189,743 consensus contigs assembled: 1,716 (Bonferroni p value ≤ 0.05) were DETs detected in different tissues of samples from different locations, 210 of them (fold change ≥ | 100|) in the same tissue of samples from a different location, and 665 (fold change ≥ | 4|) regardless of the tissue in samples from a different location. Site-specific DETs in Cochamó (169) and Yaldad (150) in candidate genes controlling tolerance to temperature and salinity shifts, and biomineralization exhibit a high number of nucleotide genetic variants with regular occurrence (frequency > 99%). This novel M. chilensis transcriptome should help assessing and monitoring the impact of translocations in wild and farm-impacted mussel beds in Chiloé Island. At the same time, it would help designing effective managing practices for conservation, and translocation traceability.

Adaptive Shifts in Gene Regulation Underlie a Developmental Delay in Thermogenesis in High-Altitude Deer Mice

Aerobic performance is tied to fitness as it influences an animal's ability to find food, escape predators, or survive extreme conditions. At high altitude, where low O2 availability and persistent cold prevail, maximum metabolic heat production (thermogenesis) is an aerobic performance trait that is closely linked to survival. Understanding how thermogenesis evolves to enhance survival at high altitude will yield insight into the links between physiology, performance, and fitness. Recent work in deer mice (Peromyscus maniculatus) has shown that adult mice native to high altitude have higher thermogenic capacities under hypoxia compared with lowland conspecifics, but that developing high-altitude pups delay the onset of thermogenesis. This finding suggests that natural selection on thermogenic capacity varies across life stages. To determine the mechanistic cause of this ontogenetic delay, we analyzed the transcriptomes of thermoeffector organs-brown adipose tissue and skeletal muscle-in developing deer mice native to low and high altitude. We demonstrate that the developmental delay in thermogenesis is associated with adaptive shifts in the expression of genes involved in nervous system development, fuel/O2 supply, and oxidative metabolism pathways. Our results demonstrate that selection has modified the developmental trajectory of the thermoregulatory system at high altitude and has done so by acting on the regulatory systems that control the maturation of thermoeffector tissues. We suggest that the cold and hypoxic conditions of high altitude force a resource allocation tradeoff, whereby limited energy is allocated to developmental processes such as growth, versus active thermogenesis, during early development.

Magnitude and Predictability of pH Fluctuations Shape Plastic Responses to Ocean Acidification

AbstractPhenotypic plasticity is expected to facilitate the persistence of natural populations as global change progresses. The attributes of fluctuating environments that favor the evolution of plasticity have received extensive theoretical investigation, yet empirical validation of these findings is still in its infancy. Here, we combine high-resolution environmental data with a laboratory-based experiment to explore the influence of habitat pH fluctuation dynamics on the plasticity of gene expression in two populations of the Mediterranean mussel, Mytilus galloprovincialis. We linked differences in the magnitude and predictability of pH fluctuations in two habitats to population-specific gene expression profiles in ambient and stressful pH treatments. Our results demonstrate population-based differentiation in gene expression plasticity, whereby mussels native to a habitat exhibiting a large magnitude of pH fluctuations with low predictability display reduced phenotypic plasticity between experimentally imposed pH treatments. This work validates recent theoretical findings on evolution in fluctuating environments, suggesting that the predictability of fluctuating selection pressures may play a predominant role in shaping the phenotypic variation observed across natural populations.

Evolutionary Responses to Warming

Climate change is predicted to dramatically alter biological diversity and distributions, driving extirpations, extinctions, and extensive range shifts across the globe. Warming can also, however, lead to phenotypic or behavioural plasticity, as species adapt to new conditions. Recent genomic research indicates that some species are capable of rapid evolution as selection favours adaptive responses to environmental change and altered or novel niche spaces. New advances are providing mechanistic insights into how temperature might accelerate evolution in the Anthropocene. These discoveries highlight intriguing new research directions - such as using geothermal and polar systems combined with powerful genomic tools - that will help us to understand the processes underpinning adaptive evolution and better project how ecosystems will change in a warming world.

Comparative transcriptomics of spotted seatrout (Cynoscion nebulosus) populations to cold and heat stress

Resilience to climate change depends on a species' adaptive potential and phenotypic plasticity. The latter can enhance survival of individual organisms during short periods of extreme environmental perturbations, allowing genetic adaptation to take place over generations. Along the U.S. East Coast, estuarine-dependent spotted seatrout (Cynoscion nebulosus) populations span a steep temperature gradient that provides an ideal opportunity to explore the molecular basis of phenotypic plasticity. Genetically distinct spotted seatrout sampled from a northern and a southern population were exposed to acute cold and heat stress (5 biological replicates in each treatment and control group), and their transcriptomic responses were compared using RNA-sequencing (RNA-seq). The southern population showed a larger transcriptomic response to acute cold stress, whereas the northern population showed a larger transcriptomic response to acute heat stress compared with their respective population controls. Shared transcripts showing significant differences in expression levels were predominantly enriched in pathways that included metabolism, transcriptional regulation, and immune response. In response to heat stress, only the northern population significantly upregulated genes in the apoptosis pathway, which could suggest greater vulnerability to future heat waves in this population as compared to the southern population. Genes showing population-specific patterns of expression, including hpt, acot, hspa5, and hsc71, are candidates for future studies aiming to monitor intraspecific differences in temperature stress responses in spotted seatrout. Our findings contribute to the current understanding of phenotypic plasticity and provide a basis for predicting the response of a eurythermal fish species to future extreme temperatures.

Freshwater mussels (Unionidae) brought into captivity exhibit up-regulation of genes involved in stress and energy metabolism

Approximately two thirds of freshwater mussel species in the United States and Canada are imperiled, and populations are declining rapidly. Translocation and captive management are commonly used to mitigate losses of freshwater mussel biodiversity, but these conservation tools may result in decreased growth and increased mortality. This study uses RNA-Seq to determine how translocation into captivity affects gene expression in Amblema plicata. Mussels were collected from the Muskingum River in Ohio, USA and brought into a captive holding facility. RNA was extracted from gill tissue 11 months post translocation from mussels in captivity and the Muskingum River on the same day. RNA was sequenced on an Illumina HiSeq 2500, and differential expression analysis was performed on de novo assembled transcripts. More than 1200 transcripts were up-regulated in captive mussels, and 246 were assigned functional annotations. Many up-regulated transcripts were involved in energy metabolism and the stress response, such as heat shock proteins and antioxidants. More than 500 transcripts were down-regulated in captive mussels, and 41 were assigned functional annotations. We observed an over-representation of down-regulated transcripts associated with immune response. Our work suggests that A. plicata experienced moderate levels of stress and altered energy metabolism and immune response for at least 11 months post translocation into captivity.

Differential immune gene expression associated with contemporary range expansion in two invasive rodents in Senegal

Biological invasions are major anthropogenic changes associated with threats to biodiversity and health. However, what determines the successful establishment and spread of introduced populations remains unclear. Here, we explore several hypotheses linking invasion success and immune phenotype traits, including those based on the evolution of increased competitive ability concept. We compared gene expression profiles between anciently and recently established populations of two major invading species, the house mouse Mus musculus domesticus and the black rat Rattus rattus, in Senegal (West Africa). Transcriptome analyses identified differential expression between anciently and recently established populations for 364 mouse genes and 83 rat genes. All immune-related genes displaying differential expression along the mouse invasion route were overexpressed at three of the four recently invaded sites studied. Complement activation pathway genes were overrepresented among these genes. By contrast, no particular immunological process was found to be overrepresented among the differentially expressed genes of black rat. Changes in transcriptome profiles were thus observed along invasion routes, but with different specific patterns between the two invasive species. These changes may be driven by increases in infection risks at sites recently invaded by the house mouse, and by stochastic events associated with colonization history for the black rat. These results constitute a first step toward the identification of immune eco-evolutionary processes potentially involved in the invasion success of these two rodent species.

Genetic compensation rather than genetic assimilation drives the evolution of plasticity in response to mild warming across latitudes in a damselfly

Global warming is causing plastic and evolutionary changes in the phenotypes of ectotherms. Yet, we have limited knowledge on how the interplay between plasticity and evolution shapes thermal responses and underlying gene expression patterns. We assessed thermal reaction norm patterns across the transcriptome and identified associated molecular pathways in northern and southern populations of the damselfly Ischnura elegans. Larvae were reared in a common garden experiment at the mean summer water temperatures experienced at the northern (20°C) and southern (24°C) latitudes. This allowed a space-for-time substitution where the current gene expression levels at 24°C in southern larvae are a proxy for the expected responses of northern larvae under gradual thermal evolution to the predicted 4°C warming. Most differentially expressed genes showed fixed differences across temperatures between latitudes, suggesting that thermal genetic adaptation will mainly evolve through changes in constitutive gene expression. Northern populations also frequently showed plastic responses in gene expression to mild warming, while southern populations were much less responsive to temperature. Thermal responsive genes in northern populations showed to a large extent a pattern of genetic compensation, namely gene expression that was induced at 24°C in northern populations remained at a lower constant level in southern populations, and were associated with metabolic and translation pathways. There was instead little evidence for genetic assimilation of an initial plastic response to mild warming. Our data therefore suggest that genetic compensation rather than genetic assimilation may drive the evolution of plasticity in response to mild warming in this damselfly species.

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.

Spatio-temporal patterns in the gene expression of the calanoid copepod Temora longicornis in the Belgian part of the North Sea

Marine zooplankton are increasingly being affected by recent environmental changes, such as climate change, and respond with profound spatial relocations and shifts in phenology and physiology. In order to predict whether populations are able to persist or adapt to such new conditions, it is essential to understand the molecular basis of such adaptations, which ultimately get translated into these physiological responses. To explore variation in population gene expression across time and space, we investigated transcriptome-level profiles of the calanoid copepod Temora longicornis, that were collected at four different locations in the Belgian Part of the North Sea (BPNS) on three different time points (April, June, October) in 2018. RNA-seq analysis of field collected adults identified large seasonal differences in gene expression, mainly between spring-summer and autumn samples. The largest log-fold changes occurred in a set of genes encoding for ribosomal and myosin (heavy chain) transcripts. Enrichment analysis revealed a strong seasonal pattern in vitellogenin, cuticle and glycolytic gene expression as well. We also found a positive correlation between vitellogenin expression and densities of T. longicornis. No clear spatial variation in expression patterns was found in the BPNS. This study underlines the potential of field gene expression studies for biomonitoring purposes and the significance of considering seasonal variation in future studies.

Genomic signals found using RNA sequencing show signatures of selection and subtle population differentiation in walleye (Sander vitreus) in a large freshwater ecosystem

RNA sequencing is an effective approach for studying aquatic species yielding both physiological and genomic data. However, its population genetic applications are not well-characterized. We investigate this possible role for RNA sequencing for population genomics in Lake Winnipeg, Manitoba, Canada, walleye (Sander vitreus). Lake Winnipeg walleye represent the largest component of the second-largest freshwater fishery in Canada. In the present study, large female walleye were sampled via nonlethal gill biopsy over two years at three spawning sites representing a latitudinal gradient in the lake. Genetic variation from sequenced mRNA was analyzed for neutral and adaptive markers to investigate population structure and possible adaptive variation. We find low population divergence (F ST = 0.0095), possible northward gene flow, and outlier loci that vary latitudinally in transcripts associated with cell membrane proteins and cytoskeletal function. These results indicate that Lake Winnipeg walleye may be effectively managed as a single demographically connected metapopulation with contributing subpopulations and suggest genomic differences possibly underlying observed phenotypic differences. Despite its high cost relative to other genotyping methods, RNA sequencing data can yield physiological in addition to genetic information discussed here. We therefore argue that it is useful for addressing diverse molecular questions in the conservation of freshwater species.

Population epigenetic divergence exceeds genetic divergence in the Eastern oyster Crassostrea virginica in the Northern Gulf of Mexico

Populations may respond to environmental heterogeneity via evolutionary divergence or phenotypic plasticity. While evolutionary divergence occurs through DNA sequence differences among populations, plastic divergence among populations may be generated by changes in the epigenome. Here, we present the results of a genome-wide comparison of DNA methylation patterns and genetic structure among four populations of Eastern oyster (Crassostrea virginica) in the northern Gulf of Mexico. We used a combination of restriction site-associated DNA sequencing (RADseq) and reduced representation bisulfite sequencing (RRBS) to explore population structure, gene-wide averages of F ST, and DNA methylation differences between oysters inhabiting four estuaries with unique salinity profiles. This approach identified significant population structure despite a moderately low F ST (0.02) across the freshwater boundary of the Mississippi river, a finding that may reflect recent efforts to restore oyster stock populations. Divergence between populations in CpG methylation was greater than for divergence in F ST, likely reflecting environmental effects on DNA methylation patterns. Assessment of CpG methylation patterns across all populations identified that only 26% of methylated DNA was intergenic; and, only 17% of all differentially methylated regions (DMRs) were within these same regions. DMRs within gene bodies between sites were associated with genes known to be involved in DNA damage repair, ion transport, and reproductive timing. Finally, when assessing the correlation between genomic variation and DNA methylation between these populations, we observed population-specific DNA methylation profiles that were not directly associated with single nucleotide polymorphisms or broader gene-body mean F ST trends. Our results suggest that C. virginica may use DNA methylation to generate environmentally responsive plastic phenotypes and that there is more divergence in methylation than divergence in allele frequencies.

Additive and mostly adaptive plastic responses of gene expression to multiple stress in Tribolium castaneum

Gene expression is known to be highly responsive to the environment and important for adjustment of metabolism but there is also growing evidence that differences in gene regulation contribute to species divergence and differences among locally adapted populations. However, most studies so far investigated populations when divergence had already occurred. Selection acting on expression levels at the onset of adaptation to an environmental change has not been characterized. Understanding the mechanisms is further complicated by the fact that environmental change is often multivariate, meaning that organisms are exposed to multiple stressors simultaneously with potentially interactive effects. Here we use a novel approach by combining fitness and whole-transcriptome data in a large-scale experiment to investigate responses to drought, heat and their combination in Tribolium castaneum. We found that fitness was reduced by both stressors and their combined effect was almost additive. Expression data showed that stressor responses were acting independently and did not interfere physiologically. Since we measured expression and fitness within the same individuals, we were able to estimate selection on gene expression levels. We found that variation in fitness can be attributed to gene expression variation and that selection pressures were environment dependent and opposite between control and stress conditions. We could further show that plastic responses of expression were largely adaptive, i.e. in the direction that should increase fitness.

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.

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'.

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

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

Seedling response to water stress in valley oak (Quercus lobata) is shaped by different gene networks across populations

Drought is a major stress for plants, creating a strong selection pressure for traits that enable plant growth and survival in dry environments. Many drought responses are conserved species-wide responses, while others vary among populations distributed across heterogeneous environments. We tested how six populations of the widely distributed California valley oak (Quercus lobata) sampled from contrasting climates would differ in their response to soil drying relative to well-watered controls in a common environment by measuring ecophysiological traits in 93 individuals and gene expression (RNA-seq) in 42 individuals. Populations did not differ in their adjustment of turgor loss point during soil drying, suggesting a generalized species-wide response. Differential expression analysis identified 689 genes with a common response to treatment across populations and 470 genes with population-specific responses. Weighted gene co-expression network analysis (WGCNA) identified groups of genes with similar expression patterns that may be regulated together (gene modules). Several gene modules responded differently to water stress among populations, suggesting regional differences in gene network regulation. Populations from sites with a high mean annual temperature responded to the imposed water stress with significantly greater changes in gene module expression, indicating that these populations may be locally adapted to respond to drought. We propose that this variation among valley oak populations provides a mechanism for differential tolerance to the increasingly frequent and severe droughts in California.

Transcriptome profiling reveals exposure to predicted end-of-century ocean acidification as a stealth stressor for Atlantic cod larvae

Ocean acidification (OA), a direct consequence of increasing atmospheric CO₂ concentration dissolving in ocean waters, is impacting many fish species. Little is known about the molecular mechanisms underlying the observed physiological impacts in fish. We used RNAseq to characterize the transcriptome of 3 different larval stages of Atlantic cod (Gadus morhua) exposed to simulated OA at levels (1179 µatm CO₂) representing end-of-century predictions compared to controls (503 µatm CO₂), which were shown to induce tissue damage and elevated mortality in G. morhua. Only few genes were differentially expressed in 6 and 13 days-post-hatching (dph) (3 and 16 genes, respectively), during a period when maximal mortality as a response to elevated pCO₂ occurred. At 36 dph, 1413 genes were differentially expressed, most likely caused by developmental asynchrony between the treatment groups, with individuals under OA growing faster. A target gene analysis revealed only few genes of the universal and well-defined cellular stress response to be differentially expressed. We thus suggest that predicted ocean acidification levels constitute a "stealth stress" for early Atlantic cod larvae, with a rapid breakdown of cellular homeostasis leading to organismal death that was missed even with an 8-fold replication implemented in this study.

Dissecting the Transcriptomic Basis of Phenotypic Evolution in an Aquatic Keystone Grazer

Knowledge of the molecular basis of phenotypic responses to environmental cues is key to understanding the process of adaptation. Insights to adaptation at an evolutionary time scale can be gained by observing organismal responses before and after a shift in environmental conditions, but such observations can rarely be made. Using the ecological and genomic model Daphnia, we linked transcriptomic responses and phosphorus (P)-related phenotypic traits under high and low P availability. We mapped weighted gene coexpression networks to traits previously assessed in resurrected ancient (600 years old) and modern Daphnia pulicaria from a lake with a historic shift in P-enrichment. Subsequently, we assessed evolutionary conservation or divergence in transcriptional networks of the same isolates. We discovered highly preserved gene networks shared between ancient genotypes and their modern descendants, but also detected clear evidence of transcriptional divergence between these evolutionarily separated genotypes. Our study highlights that phenotypic evolution is a result of molecular fine-tuning on different layers ranging from basic cellular responses to higher order phenotypes. In a broader context, these findings advance our understanding how populations are able to persist throughout major environmental shifts.

Differences in induced thermotolerance among populations of Olympia oysters

An organism's ability to cope with thermal stress is an important predictor of survival in a changing climate. One way in which organisms may acclimatize to thermal stress in the short-term is through induced thermotolerance, whereby exposure to a sublethal heat shock enables the organism to subsequently survive what might otherwise be a lethal event. Whether induced thermotolerance is related to basal thermotolerance is not well understood for marine organisms. Furthermore, whether populations often differ in their capacity for induced thermotolerance is also unclear. Here, we tested for differences in basal thermotolerance and induced thermotolerance among six populations of Olympia oysters (Ostrea lurida) from three California estuaries. Oysters were raised under common-garden laboratory conditions for a generation and then exposed to two treatments (control or sublethal heat shock) followed by a spectrum of temperatures that bound the upper critical temperature in order to determine LT₅₀ (temperature at which 50% of the population dies). All populations exhibited induced thermotolerance by increasing their LT₅₀ to a similar maximum temperature when extreme thermal stress was preceded by a sublethal heat shock. However, populations differed in their basal thermotolerance and their plasticity in thermotolerance. Populations with the highest basal thermotolerance were least able to modify upper critical temperature, while the population with the lowest basal thermotolerance exhibited the greatest plasticity in the upper critical temperature. Our results highlight that populations with high basal thermotolerance may be most vulnerable to climate warming because they lack the plasticity required to adjust their upper thermal limits.

Differential responses to ocean acidification between populations of Balanophyllia elegans corals from high and low upwelling environments

Ocean acidification (OA), the global decrease in surface water pH from absorption of anthropogenic CO₂ , may put many marine taxa at risk. However, populations that experience extreme localized conditions, and are adapted to these conditions predicted in the global ocean in 2,100, may be more tolerant to future OA. By identifying locally adapted populations, researchers can examine the mechanisms used to cope with decreasing pH. One oceanographic process that influences pH is wind-driven upwelling. Here we compare two Californian populations of the coral Balanophyllia elegans from distinct upwelling regimes, and test their physiological and transcriptomic responses to experimental seawater acidification. We measured respiration rates, protein and lipid content, and gene expression in corals from both populations exposed to pH levels of 7.8 and 7.4 for 29 days. Corals from the population that experiences lower pH due to high upwelling maintained the same respiration rate throughout the exposure. In contrast, corals from the low upwelling site had reduced respiration rates, protein content and lipid-class content at low pH exposure, suggesting they have depleted their energy reserves. Using RNA-Seq, we found that corals from the high upwelling site upregulated genes involved in calcium ion binding and ion transport, most likely related to pH homeostasis and calcification. In contrast, corals from the low upwelling site downregulated stress response genes at low pH exposure. Divergent population responses to low pH observed in B. elegans highlight the importance of multi-population studies for predicting a species' response to future OA.

Synergistic Effects of Temperature and Salinity on the Gene Expression and Physiology of Crassostrea virginica

The eastern oyster, Crassostrea virginica, forms reefs that provide critical services to the surrounding ecosystem. These reefs are at risk from climate change, in part because altered rainfall patterns may amplify local fluctuations in salinity, impacting oyster recruitment, survival, and growth. As in other marine organisms, warming water temperatures might interact with these changes in salinity to synergistically influence oyster physiology. In this study, we used comparative transcriptomics, measurements of physiology, and a field assessment to investigate what phenotypic changes C. virginica uses to cope with combined temperature and salinity stress in the Gulf of Mexico. Oysters from a historically low salinity site (Sister Lake, LA) were exposed to fully crossed temperature (20°C and 30°C) and salinity (25, 15, and 7 PSU) treatments. Using comparative transcriptomics on oyster gill tissue, we identified a greater number of genes that were differentially expressed (DE) in response to low salinity at warmer temperatures. Functional enrichment analysis showed low overlap between genes DE in response to thermal stress compared with hypoosmotic stress and identified enrichment for gene ontologies associated with cell adhesion, transmembrane transport, and microtubule-based process. Experiments also showed that oysters changed their physiology at elevated temperatures and lowered salinity, with significantly increased respiration rates between 20°C and 30°C. However, despite the higher energetic demands, oysters did not increase their feeding rate. To investigate transcriptional differences between populations in situ, we collected gill tissue from three locations and two time points across the Louisiana Gulf coast and used quantitative PCR to measure the expression levels of seven target genes. We found an upregulation of genes that function in osmolyte transport, oxidative stress mediation, apoptosis, and protein synthesis at our low salinity site and sampling time point. In summary, oysters altered their phenotype more in response to low salinity at higher temperatures as evidenced by a higher number of DE genes during laboratory exposure, increased respiration (higher energetic demands), and in situ differential expression by season and location. These synergistic effects of hypoosmotic stress and increased temperature suggest that climate change will exacerbate the negative effects of low salinity exposure on eastern oysters.

Single and repetitive microplastics exposures induce immune system modulation and homeostasis alteration in the edible mussel Mytilus galloprovincialis

Seashore invertebrates such as mussels are exposed to multiple bouts of pollution related to human activities. Plastic debris originating from land-based activities are a concerning issue as they may be fragmented in smaller pieces (microplastics, < 5 mm diameter) which have an excellent potential for uptake by a large variety of animals. Here, we set out to explore the whole transcriptome profiling of Mytilus galloprovincialis associated with temporal variability of microplastics concentrations. Mussels were submitted to (i) a single 18 days-exposure to a concentration of microplastics found during pollution events (4.6 E+5 polyethylene microbeads L^-1), (ii) a recovery period to investigate the reversibility of microplastics effects and (iii) a repeated exposure to microplastics to evidence acclimation to microplastics pollution events. Overall, 18 days-exposure to microplastics was mostly associated with disruption of mussel global homeostasis resulting in the production of stress and immune-related proteins and as a consequence, a diminution of energy allocated to growth. During the recovery period, a contrasting response was observed with the activation of apoptotic processes and the up-regulation of immune-receptors and stress-related proteins (glutathione peroxidase, hsp70) in mussels previously exposed to microplastics. These divergent responses, suggest that the establishment of compensatory mechanism as an attempt to recover, is not sufficient to counteract physiological stress induced by the first exposure. Finally, the differences observed in gene expression between single and repeated exposures to microplastics suggest, under the experimental conditions tested, that mussels may be able to establish a stress-memory upon microplastics exposure.

Transcriptomic responses of Serratia liquefaciens cells grown under simulated Martian conditions of low temperature, low pressure, and CO2-enriched anoxic atmosphere

Results from previous experiments indicated that the Gram-negative α-proteobacterium Serratia liquefaciens strain ATCC 27592 was capable of growth under low temperature (0 °C), low pressure (0.7 kPa), and anoxic, CO₂-dominated atmosphere–conditions intended to simulate the near-subsurface environment of Mars. To probe the response of its transcriptome to this extreme environment, S. liquefaciens ATCC 27592 was cultivated under 4 different environmental simulations: 0 °C, 0.7 kPa, CO₂ atmosphere (Condition A); 0 °C, ~101.3 kPa, CO₂ atmosphere (Condition B); 0 °C, ~101.3 kPa, ambient N₂/O₂ atmosphere (Condition C); and 30 °C, ~101.3 kPa, N₂/O₂ atmosphere (Condition D; ambient laboratory conditions). RNA-seq was performed on ribosomal RNA-depleted total RNA isolated from triplicate cultures grown under Conditions A-D and the datasets generated were subjected to transcriptome analyses. The data from Conditions A, B, or C were compared to laboratory Condition D. Significantly differentially expressed transcripts were identified belonging to a number of KEGG pathway categories. Up-regulated genes under all Conditions A, B, and C included those encoding transporters (ABC and PTS transporters); genes involved in translation (ribosomes and their biogenesis, biosynthesis of both tRNAs and aminoacyl-tRNAs); DNA repair and recombination; and non-coding RNAs. Genes down-regulated under all Conditions A, B, and C included: transporters (mostly ABC transporters); flagellar and motility proteins; genes involved in phenylalanine metabolism; transcription factors; and two-component systems. The results are discussed in the context of Mars astrobiology and planetary protection.

Transcriptomic responses to extreme low salinity among locally adapted populations of Olympia oyster (Ostrea lurida)

The Olympia oyster (Ostrea lurida) is a foundation species inhabiting estuaries along the North American west coast. In California estuaries, O. lurida is adapted to local salinity regimes and populations differ in low salinity tolerance. In this study, oysters from three California populations were reared for two generations in a laboratory common garden and subsequently exposed to low salinity seawater. Comparative transcriptomics was then used to understand species-level responses to hyposmotic stress and population-level mechanisms underlying divergent salinity tolerances. Gene expression patterns indicate Olympia oysters are sensitive to hyposmotic stress: All populations respond to low salinity by up-regulating transcripts indicative of protein unfolding, DNA damage and cell cycle arrest after sub-lethal exposure. Among O. lurida populations, transcriptomic profiles differed constitutively and in response to low salinity. Despite two generations in common-garden conditions, transcripts encoding apoptosis modulators were constitutively expressed at significantly different levels in the most tolerant population. Expression of cell death regulators may facilitate cell fate decisions when salinity declines. Following low salinity exposure, oysters from the more tolerant population expressed a small number of mRNAs at significantly higher levels than less tolerant populations. Proteins encoded by these transcripts regulate ciliary activity within the mantle cavity and may function to prolong valve closure and reduce mortality in low salinity seawater. Collectively, gene expression patterns suggest sub-lethal impacts of hyposmotic stress in Olympia oysters are considerable and that even oysters with greater low salinity tolerance may be vulnerable to future freshwater flooding events.

Coexpression network revealing the plasticity and robustness of population transcriptome during the initial stage of domesticating energy crop Miscanthus lutarioriparius

Coexpression network revealing genes with Co-variation Expression pattern (CE) and those with Top rank of Expression fold change (TE) played different roles in responding to new environment of Miscanthus lutarioriparius. Variation in gene expression level, the product of genetic and/or environmental perturbation, determines the robustness-to-plasticity spectrum of a phenotype in plants. Understanding how expression variation of plant population response to a new field is crucial to domesticate energy crops. Weighted Gene Coexpression Network Analysis (WGCNA) was used to explore the patterns of expression variation based on 72 Miscanthus lutarioriparius transcriptomes from two contrasting environments, one near the native habitat and the other in one harsh domesticating region. The 932 genes with Co-variation Expression pattern (CE) and other 932 genes with Top rank of Expression fold change (TE) were identified and the former were strongly associated with the water use efficiency (r ≥ 0.55, P ≤ 10^-7). Functional enrichment of CE genes were related to three organelles, which well matched the annotation of twelve motifs identified from their conserved noncoding sequence; while TE genes were mostly related to biotic and/or abiotic stress. The expression robustness of CE genes with high genetic diversity kept relatively stable between environments while the harsh environment reduced the expression robustness of TE genes with low genetic diversity. The expression plasticity of CE genes was increased less than that of TE genes. These results suggested that expression variation of CE genes and TE genes could account for the robustness and plasticity of acclimation ability of Miscanthus, respectively. The patterns of expression variation revealed by transcriptomic network would shed new light on breeding and domestication of energy crops.