Transcriptomic responses to heat stress in invasive and native blue mussels (genus Mytilus): molecular correlates of invasive success

Stability and expression patterns of housekeeping genes in Mediterranean mussels (Mytilus galloprovincialis) under field investigations

The use of marine mussels as biological models encompasses a broad range of research fields, in which the application of RNA analyses disclosed novel biomarkers of environmental stress and investigated biochemical mechanisms of action. Quantitative real-time PCR (qPCR) is the gold standard for these studies, and despite its wide use and available protocols, it may be affected by technical flaws requiring reference gene data normalization. In this study, stability of housekeeping genes commonly employed as reference genes in qPCR analyses with Mytilus galloprovincialis was explored under field conditions. Mussels were collected from farms in the Northwestern Adriatic Sea. The sampling strategy considered latitudinal gradients of environmental parameters (proxied by location), gender, and their interactions with seasonality. Analyses of gene stability were performed using different algorithms. BestKeeper and geNorm agreed that combination of the ribosomal genes 18S ribosomal RNA (18S) and 28S ribosomal RNA (28S) was the best normalization strategy in the conditions tested, which agrees with available evidence. NormFinder provided different normalization strategies, involving combinations of tubulin (TUB)/28S (Gender/Season effect) or TUB/helicase (HEL) (Location/Season effect). Since NormFinder considers data grouping and computes both intra- and inter-group stability variations, it should work better with complex experimental designs and dataset structuring. Under the selected normalization strategies, expressions of the variable housekeeping genes actin (ACT) and elongation factor-1α (EF1) correlated with seasonal and latitudinal changes of abiotic environmental factors and mussel physiological status. Results point to consider ACT and EF1 expressions as molecular biomarkers of mussel general physiological status in field studies.

Heat Hardening Ameliorates Apoptotic and Inflammatory Effects Through Increased Autophagy in Mussels

The severity, frequency, and duration of extreme events, in the context of global warming, have placed many marine ecosystems at high risk. Therefore, the application of methods that can mediate the impacts of global warming on marine organisms seems to be an emerging necessity in the near term. In this context, enhancing the thermal resilience of marine organisms may be crucial for their sustainability. It has been shown that the repeated time-limited exposure of an organism to an environmental stimulus modifies its response mode, thus enhancing resilience and allowing adaptation of the physiological and developmental phenotype to environmental stress. In the present study, we investigated the "stress memory" effect caused by heat hardening on Mytilus galloprovincialis cellular pathways to identify the underlying biochemical mechanisms that enhance mussel thermal tolerance. Heat hardening resulted in increased ETS activity and ATP production and increased autophagic performance at all elevated temperatures (24 °C, 26 °C, and 28 °C). Furthermore, at these increased temperatures, apoptosis and inflammation remain at significantly lower levels in pregnant individuals than in nonhardened individuals. Autophagy, as a negative regulator of apoptosis, may lead to decreased damage to surrounding cells, which in turn alleviates inflammatory effects. In conclusion, the exposure of mussels to heat hardening seems to provide a physiological response that enhances heat tolerance and increases cell survival through increased energy production and reduced cell death and inflammatory responses. The latter can be utilized for the management and conservation of aquatic species of economic value or endangered status.

Influence of salinity on the thermal tolerance of aquatic organisms

Aquatic organisms are challenged by changes in their external environment, such as temperature and salinity fluctuations. If these variables interacted with each other, the response of organisms to temperature changes would be modified by salinity and vice versa. We tested for potential interaction between temperature and salinity effects on freshwater, brackish, and marine organisms, including algae, macrophytes, heterotrophic protists, parasites, invertebrates, and fish. We performed a meta-analysis that compared the thermal tolerance (characterised by the temperature optimum, lower and upper temperature limits, and thermal breadth) at various salinities. The meta-analysis was based on 90 articles (algae: 15; heterotrophic protists: 1; invertebrates: 43; and fish: 31). Studies on macrophytes and parasites were lacking. We found that decreasing salinity significantly increased and decreased the lower and upper temperature limits, respectively, in all groups. Thus, a lowered salinity increased the thermal sensitivity of organisms. These findings mainly reflect the response of brackish and marine organisms to salinity changes, which dominated our database. The few studies on freshwater species showed that their lower thermal limits increased and the upper thermal limits decreased with increasing salinity, albeit statistically nonsignificant. Although non-significant, the response of thermal tolerance to salinity changes differed between various organism groups. It generally decreased in the order of: algae > invertebrates > fish. Overall, our findings indicate adverse effects of salinity changes on the temperature tolerance of aquatic organisms. For freshwater species, studies are comparatively scarce and further studies on their thermal performance at various salinity gradients are required to obtain more robust evidence for interactions between salinity and temperature tolerance. Considering test conditions such as acclimation temperature and potential infection with parasites in future studies may decrease the variability in the relationship between salinity and thermal tolerance.

Can Transcriptomics Elucidate the Role of Regulation in Invasion Success?

When a species invades a novel environment, it must bridge the environment-phenotype mismatch in its new range to persist. Contemporary invasion biology research has focused on the role that trait variation and adaptation, and their underlying genomic factors, play in a species' adaptive potential, and thus facilitating invasion. Empirical studies have provided valuable insights into phenotypes that persist and arise in novel environments, coupled with 'omics tools that further the understanding of the contributions of genomic architecture in species establishment. Particularly, the use of transcriptomics to explore the role of plasticity in the initial stages of an invasion is growing. Here, we assess the role of various mechanisms relating to regulation and functional adaptation (often measured via the transcriptome) that support trait-specific plasticity in invasive species, allowing phenotypic variability without directly altering genomic diversity. First, we present a comprehensive review of the studies utilising transcriptomics in invasion biology. Second, we collate the evidence for and against the role of a range of regulatory processes in contributing to invasive species plasticity. Finally, we pose open questions in invasion biology where the use of transcriptome data may be valuable, as well as discuss the methodological limitations.

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.

Elucidating responses of the intertidal clam Ruditapes philippinarum to compound extreme oceanic events

Ocean acidification and heatwaves caused by rising CO2 affect bivalves and other coastal organisms. Intertidal bivalves are vital to benthic ecosystems, but their physiological and metabolic responses to compound catastrophic climate events are unknown. Here, we examined Manila clam (Ruditapes philippinarum) responses to low pH and heatwaves. Biochemical and gene expression demonstrated that pH and heatwaves greatly affect physiological energy enzymes and genes expression. In the presence of heatwaves, Manila clams expressed more enzymes and genes involved in physiological energetics regardless of acidity, even more so than in the presence of both. In this study, calcifying organisms' biochemical and molecular reactions are more susceptible to temperature rises than acidity. Acclimation under harsh weather conditions was consistent with thermal stress increase at lower biological organization levels. These substantial temporal biochemical and molecular patterns illuminate clam tipping points. This study helps us understand how compound extreme weather and climate events affect coastal bivalves for future conservation efforts.

Molecular and behavioural responses of the mussel Mytilus edulis exposed to a marine heatwave

Marine heatwaves (MHW) threaten marine organisms and tend to increase in frequency and intensity. We exposed the blue mussel Mytilus edulis to a MHW lasting 23 days, including two 10-d periods of thermal intensity increase of +5 °C (20 °C-25 °C) interspersed by 1 day back to 20 °C, followed by a 4-d recovery period. We investigated behaviour responses of mussels and gene expression changes relative to the circadian rhythm (Per), oxidative stress (SOD), cellular apoptosis (CASP3), energy production (ATPs), and general stress response (hsp70). Results showed that the MHW disturbed the valve activity of mussels. Particularly, mussels increased the number of valve micro-closures, showing a stressful state of organisms. Mussels also decreased Per, CASP3, ATPs, and Hsp70 gene expression. Some behavioural and molecular effects persisted after the MHW, suggesting a limited recovery capacity of individuals. This work highlighted the vulnerability of M. edulis to a realistic MHW.

The mechanisms involved in byssogenesis in Pteria penguin under different temperatures

Byssus is important for marine bivalves to adhere robustly to diverse substrates and resist environmental impacts. The winged pearl oyster, Pteria penguin, can reattach or not reattach to the same environment, which leaves the development and survival of the oyster population at risk. In this study, diverse methods were employed to evaluate the byssus quality and explore the mechanism of byssus secretion at different temperatures. The results demonstrated that oysters maintained their byssus properties at different temperatures through polyphenol oxidase (PPO) and reactive oxygen species (ROS) variation. They were both higher at 27 °C than at 21 °C. Furthermore, PPO activities of WB27 (31.78 U/g ± 1.50 U/g) were significantly higher than NB27, WB21, and NB21. Sectional observation revealed three types of vesicles, from which a novel vesicle might participate in byssogenesis as a putative metal storage particle. Moreover, cytoskeletal proteins may cooperate with cilia to transport byssal proteins, which then facilitate byssus formation under the regulation of upstream signals. Transcriptome analysis demonstrated that protein quality control, ubiquitin-mediated proteolysis, and cytoskeletal reorganization-related genes contributed to adaptation to temperature changes and byssus fabrication, and protection-related genes play a critical role in byssogenesis, byssus toughness, and durability. These results were utilized to create a byssogenesis mechanism model, to reveal the foot gland and vesicle types of P. penguin and provide new insights into adaptation to temperature changes and byssus fabrication in sessile bivalves.

Population Genetic Divergence among Worldwide Gene Pools of the Mediterranean Mussel Mytilus galloprovincialis

The Mediterranean mussel Mytilus galloprovincialis is distributed in both hemispheres either natively or introduced. The updated population genetic distribution of this species provides a useful knowledge against which future distribution shifts could be assessed. This study, performed with seven microsatellite markers and three reference species (M. edulis, M. chilensis and M. trossulus), aimed to determine the scenario of genetic divergence between 15 samples of M. galloprovincialis from 10 localities in Europe, Africa, Asia, Australia, North America and South America. In agreement with previous data, M. trossulus was the most divergent taxon of the genus, but M. chilensis appeared as an intermediate taxon between M. edulis and M. galloprovincialis, though closer to this latter. M. galloprovincialis from the Atlantic Northeast appears as the most likely source of worldwide exotic settlements instead of the previously thought Mediterranean population. The successful worldwide establishment of M. galloprovincialis suggests it is a flexible evolutionary species (FES), i.e., a species or population whose genetic background allows it to rapidly adapt to changing environments. This natural endowed plastic adaptation makes it a candidate resilient species amidst the ongoing climatic change.

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.

The environmental cellular stress response: the intertidal as a multistressor model

The wild poses a multifaceted challenge to the maintenance of cellular function. Therefore, a multistressor approach is essential to predict the cellular mechanisms which promote homeostasis and underpin whole-organism tolerance. The intertidal zone is particularly dynamic, and thus, its inhabitants provide excellent models to assess mechanisms underpinning multistressor tolerance. Here, we critically review our current understanding of the regulation of the cellular stress response (CSR) under multiple abiotic stressors in intertidal organisms and consider to what extent a multistressor approach brings us closer to understanding responses in the wild. The function of the CSR has been well documented in laboratory and field exposures with a view to understanding single-stressor thermal effects. Multistressor studies still remain relatively limited in comparison but have applied three main approaches: (i) laboratory application of multiple stressors in isolation, (ii) multiple stressors applied in combination, and (iii) field-based correlation of multiple stressors against the CSR. The application of multiple stressors in isolation has allowed the identification of putative, shared stress pathways but overlooks non-additive stressor interactions on the CSR. Combined stressor studies are relatively limited in number but already highlight variable effects on the CSR dependent upon stressor type, timing, and magnitude. Field studies have allowed the identification of responsive components of the CSR to various stressors in situ but are correlative, not causative. A combined approach involving laboratory multistressor studies linking the CSR to whole-organism tolerance as well as field studies is required if we are to understand the role of the CSR in the natural environment.

Cardiac performance and heart gene network provide dynamic responses of bay scallop Argopecten irradians irradians exposure to marine heatwaves

The increased frequency of marine heat waves (MHWs) caused by global climate change is predicted to threaten the survival of economic bivalves, therefore having severely adverse effects on local ecological communities and aquaculture production. However, the study of scallops facing MHWs is still scarce, particularly in the scallop Argopecten irradians irradians, which has a significant share of "blue foods" in northern China. In the present study, bay scallop heart was selected to detect its cardiac performance, oxidative impairment and dynamic molecular responses, accompanied by assessing survival variations of individuals in the simulated scenario of MWHs (32 °C) with different time points (0 h, 6 h, 12 h, 24 h, 3 d, 6 d and 10 d). Notably, cardiac indices heart rate (HR), heart amplitude (HA), rate-amplitude product (RAP) and antioxidant enzyme activities superoxide dismutase (SOD) and catalase (CAT) all peaked at 24 h but sharply dropped on 3 d, coinciding with mortality. Transcriptome analysis revealed that the heart actively defended against heat stress at the acute stage (<24 h) via energy supply, misfolded proteins correction and enhanced signal transduction, whereas regulation of the defense response and apoptotic process combined with twice transcription initiation were the dominant responses at the chronic stage (3-10 d). In particular, HSP70 (heat shock protein 70), HSP90 and CALR (calreticulin) in the endoplasmic reticulum were identified as the hub genes (top 5 %) in the HR-associated module via WGCNA (weighted gene co-expression network analysis) trait-module analysis, followed by characterization of their family members and diverse expression patterns under heat exposure. Furthermore, RNAi-mediated knockdown of CALR expression (after 24 h) significantly weakened the thermotolerance of scallops, as evidenced by a drop of 1.31 °C in ABT (Arrhenius break temperature) between the siRNA-injected group and the control group. Our findings elucidated the dynamic molecular responses at the transcriptome level and verified the cardiac functions of CALR in bay scallops confronted with stimulated MHWs.

Chromosome-Level Genome Assembly of the Blue Mussel Mytilus chilensis Reveals Molecular Signatures Facing the Marine Environment

The blue mussel Mytilus chilensis is an endemic and key socioeconomic species inhabiting the southern coast of Chile. This bivalve species supports a booming aquaculture industry, which entirely relies on artificially collected seeds from natural beds that are translocated to diverse physical-chemical ocean farming conditions. Furthermore, mussel production is threatened by a broad range of microorganisms, pollution, and environmental stressors that eventually impact its survival and growth. Herein, understanding the genomic basis of the local adaption is pivotal to developing sustainable shellfish aquaculture. We present a high-quality reference genome of M. chilensis, which is the first chromosome-level genome for a Mytilidae member in South America. The assembled genome size was 1.93 Gb, with a contig N50 of 134 Mb. Through Hi-C proximity ligation, 11,868 contigs were clustered, ordered, and assembled into 14 chromosomes in congruence with the karyological evidence. The M. chilensis genome comprises 34,530 genes and 4795 non-coding RNAs. A total of 57% of the genome contains repetitive sequences with predominancy of LTR-retrotransposons and unknown elements. Comparative genome analysis of M. chilensis and M. coruscus was conducted, revealing genic rearrangements distributed into the whole genome. Notably, transposable Steamer-like elements associated with horizontal transmissible cancer were explored in reference genomes, suggesting putative relationships at the chromosome level in Bivalvia. Genome expression analysis was also conducted, showing putative genomic differences between two ecologically different mussel populations. The evidence suggests that local genome adaptation and physiological plasticity can be analyzed to develop sustainable mussel production. The genome of M. chilensis provides pivotal molecular knowledge for the Mytilus complex.

Elevated aerial temperature modulates digestive enzyme activities in Mytilus californianus

The marine intertidal mussel Mytilus californianus aggregates to form beds along the Pacific shores of North America. As a sessile organism it must cope with fluctuations in temperature during low-tide aerial exposure, which elevates maintenance costs and negatively affects its overall energy budget. The function of its digestive gland is to release enzymes that break apart ingested polymers for subsequent nutrient absorption. The effects of elevated aerial warming acclimation on the functioning of digestive gland enzymes are not well studied. In this study we asked whether digestive gland carbohydases and proteases could be overstimulated in warm condition to possibly mitigate the costs related to the heat-shock response. We compared mussels acclimated to a + 9 °C heat-shock during daily low-tide aerial exposure to mussels acclimated to isothermal tidal conditions in a simulated intertidal system. The results showed fairly consistent activities of cellulase, trypsin, and amino-peptidase across tidal variation and between thermal treatments; however, amylase activity was lower in warmed versus cool mussels across low and high-tide. We also observed the expression of heat-shock genes in gill tissue during warm tidal conditions, suggestive that moderate temperatures during aerial exposure can induce a stress response.

Cold-adapted amphipod species upon heat stress: Proteomic responses and their correlation with transcriptomic responses

The cellular heat shock response (HSR) comprises transcriptomic and proteomic reactions to thermal stress. It was here addressed, how the proteomic, together with the transcriptomic HSR, relate to the thermal sensitivities of three cold-adapted but differently thermo-sensitive freshwater amphipod species. The proteomes of thermosensitive Eulimnogammarus verrucosus and thermotolerant Eulimnogammarus cyaneus, both endemic to Lake Baikal, and of thermotolerant Holarctic Gammarus lacustris were investigated upon 24 h exposure to the species-specific 10 % lethal temperatures (LT10). Furthermore, correlations of heat stress induced changes in proteomes (this study) and transcriptomes (previous study with identical experimental design) were examined. Proteomes indicated that the HSR activated processes encompassed (i) proteostasis maintenance, (ii) maintenance of cell adhesion, (iii) oxygen transport, (iv) antioxidant response, and (v) regulation of protein synthesis. Thermo-sensitive E. verrucosus showed the most pronounced proteomic HSR and the lowest correlation of transcriptomic and proteomic HSRs. For proteins related to translation (ribosomal proteins, elongation factors), transcriptomic and proteomic changes were inconsistent: transcripts were downregulated in many cases, with levels of corresponding proteins remaining unchanged. In the Eulimnogammarus species, levels of hemocyanin protein but not transcript were increased upon heat stress, suggesting a HSR also directed to enhance oxygen transport. Thermosensitive E. verrucosus showed the most pronounced relocation of transcription/translation activity to proteostasis maintenance, which may indicate that the general species-specific stability of protein structure could be a fundamental determinant of thermotolerance. By combining transcriptomic and proteomic response data, this study provides a comprehensive picture of the cellular HSR components in the studied amphipods.

Integrated transcriptomic and metabolomic analysis sheds new light on adaptation of Pinctada fucata martensii to short-term hypoxic stress

Analyses of the transcriptome and metabolome were conducted to clarify alterations of key genes and metabolites in pearl oysters following exposure to short-term hypoxic treatment. We totally detected 209 DEGs between the control and hypoxia groups. Enrichment analysis indicated the enrichment of GO terms including "oxidation-reduction process", "ECM organization", "chaperone cofactor-dependent protein refolding", and "ECM-receptor interaction" KEGG pathway by the DEGs. In addition, between the two groups, a total of 28 SDMs were identified, which were implicated in 13 metabolic pathways, such as "phenylalanine metabolism", "D-amino acid metabolism", and "aminoacyl-tRNA biosynthesis". Results suggest that pearl oysters are exposed to oxidative stress and apoptosis under short-term hypoxia. Also, pearl oysters might adapt to short-term hypoxic treatment by increasing antioxidant activity, modulating immune and biomineralization activities, maintaining protein homeostasis, and reorganizing the cytoskeleton. The results of our study help unveil the mechanisms by which pearl oysters respond adaptively to short-term hypoxia.

Antibiotic-Induced Treatments Reveal Stress-Responsive Gene Expression in the Endangered Lichen Lobaria pulmonaria

Antibiotics are primarily found in the environment due to human activity, which has been reported to influence the structure of biotic communities and the ecological functions of soil and water ecosystems. Nonetheless, their effects in other terrestrial ecosystems have not been well studied. As a result of oxidative stress in organisms exposed to high levels of antibiotics, genotoxicity can lead to DNA damage and, potentially, cell death. In addition, in symbiotic organisms, removal of the associated microbiome by antibiotic treatment has been observed to have a big impact on the host, e.g., corals. The lung lichen Lobaria pulmonaria has more than 800 associated bacterial species, a microbiome which has been hypothesized to increase the lichen's fitness. We artificially exposed samples of L. pulmonaria to antibiotics and a stepwise temperature increase to determine the relative effects of antibiotic treatments vs. temperature on the mycobiont and photobiont gene expression and the viability and on the community structure of the lichen-associated bacteria. We found that the mycobiont and photobiont highly reacted to different antibiotics, independently of temperature exposure. We did not find major differences in bacterial community composition or alpha diversity between antibiotic treatments and controls. For these reasons, the upregulation of stress-related genes in antibiotic-treated samples could be caused by genotoxicity in L. pulmonaria and its photobiont caused by exposure to antibiotics, and the observed stress responses are reactions of the symbiotic partners to reduce damage to their cells. Our study is of great interest for the community of researchers studying symbiotic organisms as it represents one of the first steps to understanding gene expression in an endangered lichen in response to exposure to toxic environments, along with dynamics in its associated bacterial communities.

Species- and site-specific circulating bacterial DNA in Subantarctic sentinel mussels Aulacomya atra and Mytilus platensis

Impacts of climate changes are particularly severe in polar regions where warmer temperatures and reductions in sea-ice covers threaten the ecological integrity of marine coastal ecosystems. Because of their wide distribution and their ecological importance, mussels are currently used as sentinel organisms in monitoring programs of coastal ecosystems around the world. In the present study, we exploited the concept of liquid biopsy combined to a logistically friendly sampling method to study the hemolymphatic bacterial microbiome in two mussel species (Aulacomya atra and Mytilus platensis) in Kerguelen Islands, a remote Subantarctic volcanic archipelago. We found that the circulating microbiome signatures of both species differ significantly even though their share the same mussel beds. We also found that the microbiome differs significantly between sampling sites, often correlating with the particularity of the ecosystem. Predictive models also revealed that both species have distinct functional microbiota, and that the circulating microbiome of Aulacomya atra was more sensitive to changes induced by acute thermal stress when compared to Mytilus platensis. Taken together, our study suggests that defining circulating microbiome is a useful tool to assess the health status of marine ecosystems and to better understand the interactions between the sentinel species and their habitat.

Interactive effects of multiple stressors on the physiological performance of the invasive mussel Mytilus galloprovincialis

The invasive mussel Mytilus galloprovincialis is a heat-tolerant species relative to its competative congener M. trossulus, that dominates warm seawater environments but it is unknown how multiple stressors (MS) may affect its physiology. Our study determined the effects of MS on the metabolic rate (MR), superoxide dismutase (SOD) antioxidant enzyme activity, and clearance rate (CR) of M. galloprovincialis. Mussels were exposed for 7 d to hyposalinity (20, 28, 34 ppt) then to heat shock (17, 20, 25 °C) after which MR and SOD activity were determined. CR was quantified following a 30 min MS exposure. We found a significant influence of MS on MR, SOD, and CR. We identified synergistic effects on MR under the most extreme treatment. SOD activity was the greatest under 20 °C exposure while CR declined under heat shock. Thus, our study suggests that mussels experiencing MS may become energy limited as MR increases and feeding rates decrease.

Autophagic event and metabolomic disorders unveil cellular toxicity of environmental microplastics on marine polychaete Hediste diversicolor

Although the hazards of microplastics (MPs) have been quite well explored, the aberrant metabolism and the involvement of the autophagy pathway as an adverse response to environmental MPs in benthic organisms are still unclear. The present work aims to assess the impact of different environmental MPs collected from the south coast of the Mediterranean Sea, composed by polyethylene (PE), polyethylene vinyl acetate (PEVA), low-density polyethylene (LDPE), high-density polyethylene (HDPE), polypropylene (PP) and polyamide (PA) on the metabolome and proteome of the marine polychaete Hediste diversicolor. As a result, all the microplastic types were detected with Raman microspectroscopy in polychaetes tissues, causing cytoskeleton damage and induced autophagy pathway manifested by immunohistochemical labeling of specific targeted proteins, through Tubulin (Tub), Microtubule-associated protein light chain 3 (LC3), and p62 (also named Sequestosome 1). Metabolomics was conducted to further investigate the metabolic alterations induced by the environmental MPs-mixture in polychaetes. A total of 28 metabolites were differentially expressed between control and MPs-treated polychaetes, which showed elevated levels of amino acids, glucose, ATP/ADP, osmolytes, glutathione, choline and phosphocholine, and reduced concentration of aspartate. These novel findings extend our understanding given the toxicity of environmental microplastics and unravel their underlying mechanisms.

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.

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.

Does the environmental history of mussels have an effect on the physiological response to additional stress under experimental conditions?

Expected effects on marine biota of the ongoing elevation of water temperature and high latitudes is of major concern when considering the reliability of coastal ecosystem production. To compare the capacity of coastal organisms to cope with a temperature increase depending on their environmental history, responses of adult blue mussels (Mytilus spp.) taken from two sites differentially exposed to chemical pollution were investigated during an experimental exposure to a thermal stress. Immune parameters were notably altered by extreme warming and transcriptional changes for a broad selection of genes were associated to the temperature increase following a two-step response pattern. Site-specific responses suggested an influence of environmental history and support the possibility of a genetic basis in the physiological response. However no meaningful difference was detected between the response of hybrids and M galloprovincialis. This study brings new information about the capacity of mussels to cope with the ongoing elevation of water temperature in these coastal ecosystems.

Molecular Responses to Thermal and Osmotic Stress in Arctic Intertidal Mussels (Mytilus edulis): The Limits of Resilience

Increases in Arctic temperatures have accelerated melting of the Greenland icesheet, exposing intertidal organisms, such as the blue mussel Mytilus edulis, to high air temperatures and low salinities in summer. However, the interaction of these combined stressors is poorly described at the transcriptional level. Comparing expression profiles of M. edulis from experimentally warmed (30 °C and 33 °C) animals kept at control (23‰) and low salinities (15‰) revealed a significant lack of enrichment for Gene Ontology terms (GO), indicating that similar processes were active under all conditions. However, there was a progressive increase in the abundance of upregulated genes as each stressor was applied, with synergistic increases at 33 °C and 15‰, suggesting combined stressors push the animal towards their tolerance thresholds. Further analyses comparing the effects of salinity alone (23‰, 15‰ and 5‰) showed high expression of stress and osmoregulatory marker genes at the lowest salinity, implying that the cell is carrying out intracellular osmoregulation to maintain the cytosol as hyperosmotic. Identification of aquaporins and vacuolar-type ATPase transcripts suggested the cell may use fluid-filled cavities to excrete excess intracellular water, as previously identified in embryonic freshwater mussels. These results indicate that M. edulis has considerable resilience to heat stress and highly efficient mechanisms to acclimatise to lowered salinity in a changing world.

Differences in chemical contaminants bioaccumulation and ecotoxicology biomarkers in Mytilus edulis and Mytilus galloprovincialis and their hybrids

The Mytilus mussels are spread all over the world and many related species coexist in several areas and can produce hybrid offspring. Mussels have been used for decades in national and international programs to monitor chemical contamination in the environment. Differences in bioaccumulation and biotransformation abilities between species and their hybrids should be evaluated to assess the comparability of the results obtained within the international biomonitoring programs. The objective of this study was to characterize bioaccumulation abilities and biomarker responses in Mytilus edulis, Mytilus galloprovincialis and their hybrids via an in situ transplantation experimentation on their progenies. Four mussel groups (M. edulis, M. galloprovincialis and two hybrids batches) issued from genetically characterized parents were transplanted for one year in Charente Maritime (France) to ensure their exposure to identical sources of contamination. The bioaccumulation of several families of contaminants (trace metals, polycyclic aromatic hydrocarbons, polybrominated diphenyl ethers, polychlorinated biphenyls), the response of several biomarkers (DNA strand breaks level, lysosomal membrane stability, metallothionein content, acetylcholine esterase activity) and some physiological parameters (growth, mortality, gonadal development), were analyzed. Differences were observed between species, however they were contaminant-specific. Variations in contaminants levels were observed between progenies, with higher levels of Cu, PBDE, PCB in M. edulis, and higher levels of Cd, Hg, Zn in M galloprovincialis. This study demonstrated that variations in contaminant bioaccumulation and different biomarker responses exist between Mytilus species in the field. Data on species or the presence of hybrid individuals (or introgression) is an important additional parameter to add to biomonitoring programs databases.

Gene expression responses to thermal shifts in the endangered lichen Lobaria pulmonaria

Anthropogenic climate change has led to unprecedented shifts in temperature across many ecosystems. In a context of rapid environmental changes, acclimation is an important process as it may influence the capacity of organisms to survive under novel thermal conditions. Mechanisms of acclimation could involve upregulation of stress response genes involved in protein folding, DNA damage repair and the regulation of signal transduction genes, along with a simultaneous downregulation of genes involved in growth or the cell cycle, in order to maintain cellular functions and equilibria. We transplanted Lobaria pulmonaria lichens originating from different forests to determine the relative effects of long-term acclimation and genetic factors on the variability in expression of mycobiont and photobiont genes. We found a strong response of the mycobiont and photobiont to high temperatures, regardless of sample origin. The green-algal photobiont had an overall lower response than the mycobiont. Gene expression of both symbionts was also influenced by acclimation to transplantation sites and by genetic factors. L. pulmonaria seems to have evolved powerful molecular pathways to deal with environmental fluctuations and stress and can acclimate to new habitats by transcriptomic convergence. Although L. pulmonaria has the molecular machinery to counteract short-term thermal stress, survival of lichens such as L. pulmonaria depends mostly on their long-term positive carbon balance, which can be compromised by higher temperatures and reduced precipitation, and both these outcomes have been predicted for Central Europe in connection with global climate change.

Transcriptomic Profile of the Cockle Cerastoderma edule Exposed to Seasonal Diarrhetic Shellfish Toxin Contamination

Bivalves constitute an important source of proteins for human consumption, but some accumulate biotoxins such as diarrhetic shellfish toxins (DSTs), constituting a risk to human health. The cockle Cerastoderma edule is one of the most important species harvested in the Portuguese coast but also one of the most affected species due to recurrent DSTs exposure. However, little is known regarding the effects of the toxins produced by blooming dinoflagellates on C. edule. Herein, we explore the Differentially Expressed Genes (DEGs) of two tissues (gills and digestive gland) from wild cockles sampled in Portugal, through their whole transcriptomic response in two different seasons (exposed and not exposed to DSTs). The de novo transcriptome assembly returned 684,723 contigs, N50 of 1049, and 98.53% completeness. Altogether, 1098 DEGs were identified, of which 353 DEGs were exclusive for the digestive gland, 536 unique for the gills and 209 DEGs were common. Among DEGs were identified known DSTs-biomarkers including glutathione peroxidase, glutathione S-transferase, superoxide dismutase, cytochrome P450, ABC transporters, actin and tubulin-related proteins, Heat shock proteins and complement C1Q-like proteins. This study provides the first transcriptomic profile of C. edule, giving new insights about its molecular responses under different environmental conditions of DSTs exposure.

An integrated, multi-level analysis of thermal effects on intertidal molluscs for understanding species distribution patterns

Elucidating the physiological mechanisms that underlie thermal stress and discovering how species differ in capacities for phenotypic acclimatization and evolutionary adaptation to this stress is critical for understanding current latitudinal and vertical distribution patterns of species and for predicting their future state in a warming world. Such mechanistic analyses require careful choice of study systems (species and temperature-sensitive traits) and design of laboratory experiments that reflect the complexities of in situ conditions. Here, we critically review a wide range of studies of intertidal molluscs that provide mechanistic accounts of thermal effects across all levels of biological organization - behavioural, organismal, organ level, cellular, molecular, and genomic - and show how temperature-sensitive traits govern distribution patterns and capacities for coping with thermal stress. Comparisons of congeners from different thermal habitats are especially effective means for identifying adaptive variation. We employ these mechanistic analyses to illustrate how species differ in the severity of threats posed by rising temperature. Counterintuitively, we show that some of the most heat-tolerant species may be most threatened by increases in temperatures because of their small thermal safety margins and minimal abilities to acclimatize to higher temperatures. We discuss recent molecular biological and genomic studies that provide critical foundations for understanding the types of evolutionary changes in protein structure, RNA secondary structure, genome content, and gene expression capacities that underlie adaptation to temperature. Duplication of stress-related genes, as found in heat-tolerant molluscs, may provide enhanced capacity for coping with higher temperatures. We propose that the anatomical, behavioural, physiological, and genomic diversity found among intertidal molluscs, which commonly are of critical importance and high abundance in these ecosystems, makes this group of animals a highly appropriate study system for addressing questions about the mechanistic determinants of current and future distribution patterns of intertidal organisms.

Molecular response of a sub-antarctic population of the blue mussel (Mytilus edulis platensis) to a moderate thermal stress

The Kerguelen Islands (49°26'S, 69°50'E) represent a unique environment due to their geographical isolation, which protects them from anthropogenic pollution. The ability of the endemic mussel, part of the Mytilus complex, to cope with moderate heat stress was explored using omic tools. Transcripts involved in six major metabolic functions were selected and the qRT-PCR data indicated mainly changes in aerobic and anaerobic energy metabolism and stress response. Proteomic comparisons revealed a typical stress response pattern with cytoskeleton modifications and elements suggesting increased energy metabolism. Results also suggest conservation of protein homeostasis by the long-lasting presence of HSP while a general decrease in transcription is observed. The overall findings are consistent with an adaptive response to moderate stresses in mussels in good physiological condition, i.e. living in a low-impact site, and with the literature concerning this model species. Therefore, local blue mussels could be advantageously integrated into biomonitoring strategies, especially in the context of Global Change.

Treatment with ascorbic acid normalizes the aerobic capacity, antioxidant defence, and cell death pathways in thermally stressed Mytilus galloprovincialis

Considering temperature's upcoming increase due to climate change, combined with the fact that Mediterranean mussels Mytilus galloprovincialis (Lamarck, 1819) live at their upper limits [critical temperatures (Tc) beyond 25 °C], we cannot be sure of this species' sustainable future in the Mediterranean Sea. Deviation from optimum temperatures leads to cellular damage due to oxidative stress. Although ascorbic acid (AA) is a major scavenger of reactive oxygen species (ROS), its capacity to minimize oxidative stress effects is scarcely studied in aquatic organisms. Thus, treatment with 5 mM and 10 mM AA of thermally stressed molluscs had been employed in order to examine its antioxidant capacity. While 5 mM had no effect, 10 mM normalized COX1 and ND2 relative mRNA levels, and superoxide dismutase (SOD), catalase, and glutathione reductase (GR) enzymatic activity levels in both examined tissues: posterior adductor muscle (PAM) and mantle. ATP levels, probably providing the adequate energy for antioxidant defence in thermally stressed mussels, is also normalized under 10 mM AA treatment. Moreover, autophagic indicators such as LC3 II/I and SQSTM1/p62 levels are normalized, indicating autophagy amelioration. Apoptosis also seems to be inhibited since both Bax/Bcl-2 and cleaved caspase substrate levels decrease with 10 mM AA treatment. Therefore, treatment of mussels with AA seems to produce threshold effects, although the precise underlying mechanisms must be elucidated in future studies. These findings show that treatment of mussels with effective antioxidants can be useful as a strategic approach for the reduction of the deleterious effects on mussels' summer mortality in aquaculture zones.

Resilience in Greenland intertidal Mytilus: The hidden stress defense

The Arctic is experiencing particularly rapid rates of warming, consequently invasive boreal species are now able to survive the less extreme Arctic winter temperatures. Whilst persistence of intertidal and terrestrial species in the Arctic is primarily determined by their ability to tolerate the freezing winters, air temperatures in the Arctic summer can reach 36 °C in the intertidal, which is beyond the upper thermal limits of many marine species. This is normally lethal for the conspicuous ecosystem engineer Mytilus edulis. Transcriptomic analyses were undertaken on both in situ collected and experimentally warmed animals to understand whether M. edulis is able to tolerate these very high summer temperatures. Surprisingly there was no significant enrichment for Gene Ontology terms (GO) when comparing the inner and outer fjord intertidal animals with outer fjord subtidal (control) animals, representing animals collected at 27 °C, 19 °C and 3 °C respectively. This lack of differentiation indicated a wide acclimation ability in this species. Conversely, significant enrichment for processes such as signal transduction, cytoskeleton and cellular protein modification was identified in the expression profiles of the 22 °C and 32 °C experimentally heated animals. This difference in gene expression between in situ collected and experimentally warmed animals was almost certainly due to the former being acclimated to a fluctuating, but predictable, temperature regime, which has increased their thermal tolerances. Interestingly, there was no evidence for enrichment of the classical cellular stress response in any of the animals sampled. Identification of a massive expansion of the HSPA12 heat shock protein 70 kDa gene family presented the possibility of these genes acting as intertidal regulators underpinning thermal resilience. This expansion has resulted in a modified cellular stress response, as an evolutionary adaptation to the rigour of the invasive intertidal life style. Thus, M. edulis appear to have considerable capacity to withstand the current rates of Arctic warming, and the very large attendant thermal variation.

New challenges of marine ecotoxicology in a global change context

Currently, research agenda in marine ecotoxicology is facing new challenges with the emergence of newly and complex synthetized chemicals. The study of the fate and adverse effects of toxicants remains increasingly complicated with global change events. Ecotoxicology had provided for a decades, precious scientific data and knowledge but also technical and management tools for the environmental community. Regarding those, it is necessary to update methodologies dealing with these issues such as combined effect of conventional and emergent stressors and global changes. In this point of view article, we discuss one hand the new challenges of ecotoxicology in this context, and in the other hand, the need of updating agenda and methodologies currently used in monitoring programs and finally recommendations and future research needs. Among recommendations, it could be cited the necessity to perform long-term experiments, the standardization of sentinel species and taking benefit from baseline studies and omics technologies.

Freshening increases the susceptibility to heat stress in intertidal mussels (Mytilus edulis) from the Arctic

Temperatures in the Arctic are increasing at a faster pace than at lower latitudes resulting in range expansion of boreal species. In Greenland, the warming also drives accelerating melt of the Greenland Ice Sheet resulting in more meltwater entering Greenland fjords in summer. Our aim was to determine if increasing summer temperatures combined with lower salinity can induce the expression of stress-related proteins, for example, heat shock protein, in boreal intertidal mussels in Greenland, and whether low salinity reduces the upper thermal limit at which mortality occurs. We conducted a mortality experiment, using 12 different combinations of salinity and air temperature treatments during a simulated tidal regime, and quantified the change in mRNA levels of five stress-related genes (hsp24, hsp70, hsp90, sod and p38) in surviving mussels to discern the level of sublethal stress. Heat-induced mortality occurred in mussels exposed to an air temperature of 30°C and mortality was higher in treatments with lowered salinity (5 and 15‰), which confirms that low habitat salinity decreases the upper thermal limit of Mytilus edulis. The gene expression analysis supported the mortality results, with the highest gene expression found at combinations of high temperature and low salinity. Combined with seasonal measurements of intertidal temperatures in Greenland, we suggest heat stress occurs in low salinity intertidal area, and that further lowered salinity in coastal water due to increased run-off can make intertidal bivalves more susceptible to summer heat stress. This study thus provides an example of how different impacts of climate warming can work synergistically to stress marine organisms.

Transcriptomic analyses provide insights into the adaptive responses to heat stress in the ark shells, Scapharca subcrenata

The ark shell, Scapharca subcrenata, is susceptible to high temperature which may lead to mass mortality in hot summers. Herein, we conducted the transcriptomic analyses of haemocytes in ark shells under thermal stress, to reveal the underlying molecular mechanisms of heat resistance in these animals. The results showed that a total of 7773, 11,500 and 13,046 unigenes were expressed differentially at 12, 24 and 48 h post thermal stress, respectively. The expression levels of key DEGs as revealed by RNA-seq were confirmed by quantitative real-time PCR. GO and KEGG enrichment analyses showed that the DEGs were mainly associated with apoptosis, NF-kappa B signaling pathway, TNF signaling pathway and RIG-I-like receptor signaling pathway. Among the DEGs, 40 were candidate heat stress response-related genes and 169 were identified to be involved in antioxidant defense, cell detoxification, protein metabolism and endoplasmic reticulum stress responses. It seemed that ark shells may adapt to short term thermal stress through regulation of protein metabolism, DNA replication and anti-apoptotic system. However, if the stress sustains, it may cause irreparable injury gradually in the animals due to oxygen limitation and metabolic dysregulation. Noteworthily, the expression of DEGs involved in protein biosynthesis and proteolysis was significantly elevated in ark shells under heat stress. These findings may provide preliminary insights into the molecular response of ark shells to acute thermal stress and lay the groundwork for marker-assisted selection of heat-resistant strains in S. subcrenata.

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.

Omics approaches for conservation biology research on the bivalve Chamelea gallina

The striped venus (Chamelea gallina) is an important economic resource in the Mediterranean Basin; this species has exhibited a strong quantitative decline in the Adriatic Sea. The aim of this work was to provide a comprehensive view of the biological status of C. gallina to elucidate the bioecological characteristics and genetic diversity of wild populations. To the best of our knowledge, this investigation is the first to perform a multidisciplinary study on C. gallina based on two omics approaches integrated with histological, ecotoxicological, and chemical analyses and with the assessment of environmental parameters. The results obtained through RNA sequencing indicated that the striped venus has a notable ability to adapt to different environmental conditions. Moreover, the stock reduction exhibited by this species in the last 2 decades seems not to have negatively affected its genetic diversity. Indeed, the high level of genetic diversity that emerged from our ddRAD dataset analyses is ascribable to the high larval dispersal rate, which might have played a “compensatory role” on local fluctuations, conferring to this species a good adaptive potential to face the environmental perturbations. These findings may facilitate the efforts of conservation biologists to adopt ad hoc management plans for this fishery resource.

The effects of cold stress on Mytilus species in the natural environment

Environmental stressors induce changes in marine mussels from molecular (e.g., neurotransmitter and chaperone concentration, and expression of immune- and heat-shock protein-related genes) to physiological (e.g., filtration and heart rates, the number of circulating hemocytes) levels. Temperature directly affects the biogeographic distribution of mussels. Chaperones might form an essential part of endogenous protective mechanisms for the adaptation of these animals to low temperatures in nature. Here, we review the available studies dealing with cold stress responses of Mytilidae family members in their natural environment.

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.

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

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

Expected ocean warming conditions significantly alter the transcriptome of developing postlarval American lobsters (Homarus americanus): Implications for energetic trade-offs

The American lobster (Homarus americanus) is one of the most iconic and economically valuable fishery species in the Northwestern Atlantic. Surface ocean temperatures are rapidly increasing across much of the species' range, raising concern about resiliency in the face of environmental change. Warmer temperatures accelerate rates of larval development and enhance survival to the postlarval stage, but the potential costs at the molecular level have rarely been addressed. We explored how exposure to current summer temperatures (16 °C) or temperature regimes mimicking projected moderate or extreme warming scenarios (18 °C and 22 °C, respectively) for the Gulf of Maine during development influences the postlarval lobster transcriptome. After de novo assembling the transcriptome, we identified 2542 differentially expressed (DE; adjusted p < 0.05) transcripts in postlarvae exposed to 16 °C vs. 22 °C, and 422 DE transcripts in postlarvae reared at 16 °C vs. 18 °C. Lobsters reared at 16 °C significantly over-expressed transcripts related to cuticle formation and the immune response up to 14.4- and 8.5-fold respectively, relative to those reared at both 18 °C and 22 °C. In contrast, the expression of transcripts affiliated with metabolism increased up to 7.1-fold as treatment temperature increased. These results suggest that lobsters exposed to projected warming scenarios during development experience a shift in the transcriptome that reflects a potential trade-off between maintaining immune defenses and sustaining increased physiological rates under a warming environment. This could have major implications for post-settlement survival through increased risk of mortality due to disease and/or starvation if energetic demands cannot be met.

Integrating GWAS and Transcriptomics to Identify the Molecular Underpinnings of Thermal Stress Responses in Drosophila melanogaster

Thermal tolerance of an organism depends on both the ability to dynamically adjust to a thermal stress and preparatory developmental processes that enhance thermal resistance. However, the extent to which standing genetic variation in thermal tolerance alleles influence dynamic stress responses vs. preparatory processes is unknown. Here, using the model species Drosophila melanogaster, we used a combination of Genome Wide Association mapping (GWAS) and transcriptomic profiling to characterize whether genes associated with thermal tolerance are primarily involved in dynamic stress responses or preparatory processes that influence physiological condition at the time of thermal stress. To test our hypotheses, we measured the critical thermal minimum (CTmin) and critical thermal maximum (CTmax) of 100 lines of the Drosophila Genetic Reference Panel (DGRP) and used GWAS to identify loci that explain variation in thermal limits. We observed greater variation in lower thermal limits, with CTmin ranging from 1.81 to 8.60°C, while CTmax ranged from 38.74 to 40.64°C. We identified 151 and 99 distinct genes associated with CTmin and CTmax, respectively, and there was strong support that these genes are involved in both dynamic responses to thermal stress and preparatory processes that increase thermal resistance. Many of the genes identified by GWAS were involved in the direct transcriptional response to thermal stress (72/151 for cold; 59/99 for heat), and overall GWAS candidates were more likely to be differentially expressed than other genes. Further, several GWAS candidates were regulatory genes that may participate in the regulation of stress responses, and gene ontologies related to development and morphogenesis were enriched, suggesting many of these genes influence thermal tolerance through effects on development and physiological status. Overall, our results suggest that thermal tolerance alleles can influence both dynamic plastic responses to thermal stress and preparatory processes that improve thermal resistance. These results also have utility for directly comparing GWAS and transcriptomic approaches for identifying candidate genes associated with thermal tolerance.

Effects of balancing selection and microhabitat temperature variations on heat tolerance of the intertidal black mussel Septifer virgatus

Realistic assessments of the impacts of global warming on population extinction risk are likely to require an integrated analysis of the roles of standing genetic variation, microhabitat thermal complexity, and the inter-individual variation of heat tolerance due to both genetic differences and seasonal acclimatization effects. Here, we examine whether balancing selection and microhabitat temperature heterogeneity can interact to enhance the population persistence to thermal stress for the black mussel Septifer virgatus. We deployed biomimetic data loggers on the shore to measure the microhabitat-specific thermal variation from June 2014 to April 2016. Thermal tolerance of specimens was indexed by measuring effects of temperature on heart rate. Genotyping of specimens was performed using double digestion restriction association RADSeq (ddRADseq). Our results show that inter-individual variations in thermal tolerance correlate significantly with genetic differences at some specific gene loci, and that heterozygotes have higher thermal tolerances than homozygotes. The observed seasonal changes in genotype frequency suggest that these loci are under balancing selection. The ability of thermally resistant heterozygotes to survive in sun-exposed microhabitats acts to balance the loss of homozygotes during summer and enable the persistence of genetic polymorphisms. Population persistence of the mussel is also facilitated by the micro-scale variation in temperature, which provides refugia from thermal stress. Our results emphasize that inter-individual variation in thermal tolerance and in microhabitat heterogeneity in temperature are important for the persistence of populations in rocky shore habitats.

Gene regulation underpinning increased thermal tolerance in a laboratory-evolved coral photosymbiont

Small increases in ocean temperature can disrupt the obligate symbiosis between corals and dinoflagellate microalgae, resulting in coral bleaching. Little is known about the genes that drive the physiological and bleaching response of algal symbionts to elevated temperature. Moreover, many studies to-date have compared highly divergent strains, making it challenging to accredit specific genes to contrasting traits. Here, we compare transcriptional responses at ambient (27°C) and bleaching-relevant (31°C) temperatures in a monoclonal, wild-type (WT) strain of Symbiodiniaceae to those of a selected-strain (SS), derived from the same monoclonal culture and experimentally evolved to elevated temperature over 80 generations (2.5 years). Thousands of genes were differentially expressed at a log fold-change of >8 between the WT and SS over a 35 days temperature treatment period. At 31°C, WT cells exhibited a temporally unstable transcriptomic response upregulating genes involved in the universal stress response such as molecular chaperoning, protein repair, protein degradation and DNA repair. Comparatively, SS cells exhibited a temporally stable transcriptomic response and downregulated many stress response genes that were upregulated by the WT. Among the most highly upregulated genes in the SS at 31°C were algal transcription factors and a gene probably of bacterial origin that encodes a type II secretion system protein, suggesting interactions with bacteria may contribute to the increased thermal tolerance of the SS. Genes and functional pathways conferring thermal tolerance in the SS could be targeted in future genetic engineering experiments designed to develop thermally resilient algal symbionts for use in coral restoration and conservation.

Molecular mechanisms underlying the effects of temperature increase on Mytilus sp. and their hybrids at early larval stages

The present work aims to investigate the effects of water temperature increase on Mytilus galloprovincilis and Mytilus edulis pure larvae (PG, PE) and their hybrids (HFG, HFE). D-larvae were maintained at 18 °C or exposed to a higher temperature of 22 °C for 48 h. Initially, Embryotoxicity test was evaluated. Second, a transcriptomic analysis using a recently developed microarray platform was applied to determine the main biological processes involved in early life stages responses to temperature increase. Finally, an immunofluorescence investigation was performed to bridge the gap between transcriptomic regulation and the real changes at cellular/tissue levels. Embryotoxicity test revealed a higher sensitivity of M. edulis (PE) D-larvae as well as hybrids from females M. edulis (HFE) to temperature increase, with the highest rate of larval malformations. Transcriptomic results indicated a lack of an adequate heat shock protein (Hsp) response in PE and HFE larvae (the high expression was observed in PG larvae); the differential expression of gene involved in translation, energy metabolism and oxidative stress response may contribute to explain the observed complex alterations in the studied conditions. As revealed by immunohistochemistry, cytoskeleton proteins changes associated with a drastic decrease of Histidine-Rich Glycoprotein (HRG) may elucidate the larval abnormalities in shell development observed for PE and HFE larvae. Overall, the results indicate that each type of pure larva (PG and PE) and their respective female hybrid (HFG and HFE) react similarly to the temperature increase. Our data should be carefully considered in view of the water temperature increase in marine ecosystems and especially for the mussel's species in confluence zones.

High-throughput quantification of protein structural change reveals potential mechanisms of temperature adaptation in Mytilus mussels

Background

Temperature exerts a strong influence on protein evolution: species living in thermally distinct environments often exhibit adaptive differences in protein structure and function. However, previous research on protein temperature adaptation has focused on small numbers of proteins and on proteins adapted to extreme temperatures. Consequently, less is known about the types and quantity of evolutionary change that occurs to proteins when organisms adapt to small shifts in environmental temperature. In this study, these uncertainties were addressed by developing software that enabled comparison of structural changes associated with temperature adaptation (hydrogen bonding, salt bridge formation, and amino acid use) among large numbers of proteins from warm- and cold-adapted species of marine mussels, Mytilus galloprovincialis and Mytilus trossulus, respectively.

Results

Small differences in habitat temperature that characterize the evolutionary history of Mytilus mussels were sufficient to cause protein structural changes consistent with temperature adaptation. Hydrogen bonds and salt bridges that increase stability and protect against heat-induced denaturation were more abundant in proteins from warm-adapted M. galloprovincialis compared with proteins from cold-adapted M. trossulus. These structural changes were related to deviations in the use of polar and charged amino acids that facilitate formation of hydrogen bonds and salt bridges within proteins, respectively. Enzymes, in particular those within antioxidant and cell death pathways, were over-represented among proteins with the most hydrogen bonds and salt bridges in warm-adapted M. galloprovincialis. Unlike extremophile proteins, temperature adaptation in Mytilus proteins did not involve substantial changes in the number of hydrophobic or large volume amino acids, nor in the content of glycine or proline.

Conclusions

Small shifts in organism temperature tolerance, such as that needed to cope with climate warming, may result from structural and functional changes to a small percentage of the proteome. Proteins in which function is dependent on large conformational change, notably enzymes, may be particularly sensitive to temperature perturbation and represent foci for natural selection. Protein temperature adaptation can occur through different types and frequencies of structural change, and adaptive mechanisms used to cope with small shifts in habitat temperature appear different from mechanisms used to retain protein function at temperature extremes.

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.

Mitochondrial thermo-sensitivity in invasive and native freshwater mussels

Climate change is impacting many, if not all, forms of life. Increases in extreme temperature fluctuations and average temperatures can cause stress, particularly in aquatic sessile ectotherms such as freshwater mussels. However, some species seem to thrive more than others in face of temperature-related stressors. Thermal tolerance may, for example, explain the success of invasive species. It is also known that mitochondria can play a key role in setting an ectothermic species' thermal tolerance. In this study, we aimed to characterize the mitochondrial thermo-tolerance in invasive and endemic freshwater mussels. With the use of high-resolution respirometry, we analyzed the mitochondrial respiration of two freshwater bivalve species exposed to a broad range of temperatures. We noticed that the invasive dreissenid Dreissena bugensis possessed a less thermo-tolerant mitochondrial metabolism than the endemic unionid Elliptio complanata This lack of tolerance was linked with a more noticeable aerobic metabolic depression at elevated temperatures. This decrease in mitochondrial metabolic activity was also linked with an increase in leak oxygen consumption as well as a stable maintenance of the activity of cytochrome c oxidase in both species. These findings may be associated both with the species' life history characteristics, as D. bugensis is more adapted to unstable habitats, in which selection pressures for resistance adaptations are reduced. Our findings add to the growing body of literature characterizing the mitochondrial metabolism of many aquatic ectotherms in our changing world.

The cellular stress response and temperature: Function, regulation, and evolution

The cellular stress response (CSR) is critical for enabling organisms to cope with thermal damage to proteins, nucleic acids, and membranes. It is a graded response whose properties vary with the degree of cellular damage. Molecular damage has positive, as well as negative, function-perturbing effects. Positive effects include crucial regulatory interactions that orchestrate involvement of the different components of the CSR. Thermally unfolded proteins signal for rapid initiation of transcription of genes encoding heat shock proteins (HSPs), central elements of the heat shock response (HSR). Thermal disruption of messenger RNA (mRNA) secondary structures in untranslated regions leads to the culling of the mRNA pool: thermally labile mRNAs for housekeeping proteins are degraded by exonucleases; heat-resistant mRNAs for stress proteins like HSPs then can monopolize the translational apparatus. Thus, proteins and RNA function as "cellular thermometers," and evolved differences in their thermal stabilities enable rapid initiation of the CSR whenever cell temperature rises significantly above the normal thermal range of a species. Covalent DNA damage, which may result from increased production of reactive oxygen species, is temperature-dependent; its extent may determine cellular survival. High levels of stress that exceed capacities for molecular repair can lead to proteolysis, inhibition of cell division, and programmed cell death (apoptosis). Onset of these processes may occur later in the stress period, after initiation of the HSR, to allow HSPs opportunity to restore protein homeostasis. Delay of these energy costly processes may also result from shortfalls in availability of adenosine triphosphate and reducing power during times of peak stress.