Life in the intertidal: Cellular responses, methylation and epigenetics

Genome-wide identification, evolution of DNA methyltransferases and their expression under salinity stress in Larimichthys crocea

DNA methyltransferases (Dnmts) are responsible for DNA methylation which influences patterns of gene expression and plays a crucial role in response to environmental changes. In this study, 7 LcDnmt genes were identified in the genome of large yellow croaker (Larimichthys crocea). The comprehensive analysis was conducted on gene structure, protein and location site of LcDnmts. LcDnmt proteins belonged to three groups (Dnmt1, Dnmt2, and Dnmt3) according to their conserved domains and phylogenetic analysis. Although Dnmt3 can be further divided into three sub groups (Dnmt3a, Dnmt3b, and Dnmt3l), there is no Dnmnt3l member in the large yellow croaker. Phylogenetic analysis revealed that the Dnmt family was highly conserved in teleosts. Expression patterns derived from the RNA-seq, qRT-PCR and Western blot analysis revealed that 2 LcDnmt genes (LcDnmt1 and LcDnmt3a2) significantly regulated under salinity stress in the liver, which was found to be dominantly expressed in the intestine and brain, respectively. These two genes may play an important role in the salinity stress of large yellow croaker and represent candidates for future functional analysis. Our results revealed the conservation of Dnmts during evolution and indicated a potential role of Dnmts in epigenetic regulation of response to salinity stress.

Multi-omics for studying and understanding polar life

Polar ecosystems are experiencing amongst the most rapid rates of regional warming on Earth. Here, we discuss 'omics' approaches to investigate polar biodiversity, including the current state of the art, future perspectives and recommendations. We propose a community road map to generate and more fully exploit multi-omics data from polar organisms. These data are needed for the comprehensive evaluation of polar biodiversity and to reveal how life evolved and adapted to permanently cold environments with extreme seasonality. We argue that concerted action is required to mitigate the impact of warming on polar ecosystems via conservation efforts, to sustainably manage these unique habitats and their ecosystem services, and for the sustainable bioprospecting of novel genes and compounds for societal gain.

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.

Roles of multi-level temperature-adaptive responses and microhabitat variation in establishing distributions of intertidal species

How intertidal species survive their harsh environment and how best to evaluate and forecast range shifts in species distribution are two important and closely related questions for intertidal ecologists and global change biologists. Adaptive variation in responses of organisms to environmental change across all levels of biological organization - from behavior to molecular systems - is of key importance in setting distribution patterns, yet studies often neglect the interactions of diverse types of biological variation (e.g. differences in thermal optima owing to genetic and acclimation-induced effects) with environmental variation, notably at the scale of microhabitats. Intertidal species have to cope with extreme and frequently changing thermal stress, and have shown high variation in thermal sensitivities and adaptive responses at different levels of biological organization. Here, I review the physiological and biochemical adaptations of intertidal species to environmental temperature on multiple spatial and temporal scales. With fine-scale datasets for the thermal limits of individuals and for environmental temperature variation at the microhabitat scale, we can map the thermal sensitivity for each individual in different microhabitats, and then scale up the thermal sensitivity analysis to the population level and, finally, to the species level by incorporating physiological traits into species distribution models. These more refined mechanistic models that include consideration of physiological variations have higher predictive power than models that neglect these variations, and they will be crucial to answering the questions posed above concerning adaptive mechanisms and the roles they play in governing distribution patterns in a rapidly changing world.

Intertidal limits shape covariation between metabolic plasticity, oxidative stress and telomere dynamics in Pacific oyster (Crassostrea gigas)

In intertidal zones, species such as sessile shellfish exhibit extended phenotypic plasticity to face rapid environmental changes, but whether frequent exposure to intertidal limits of the distribution range impose physiological costs for the animal remains elusive. Here, we explored how phenotypic plasticity varied along foreshore range at multiple organization levels, from molecular to cellular and whole organism acclimatization, in the Pacific oyster (Crassostrea gigas). We exposed 7-month-old individuals for up to 16 months to three foreshore levels covering the vertical range for this species, representing 20, 50 and 80% of the time spent submerged monthly. Individuals at the upper range limit produced energy more efficiently, as seen by steeper metabolic reactive norms and unaltered ATP levels despite reduced mitochondrial density. By spending most of their time emerged, oysters mounted an antioxidant shielding concomitant with lower levels of pro-oxidant proteins and postponed age-related telomere attrition. Instead, individuals exposed at the lower limit range near subtidal conditions showed lower energy efficiencies, greater oxidative stress and shorter telomere length. These results unraveled the extended acclimatization strategies and the physiological costs of living too fast in subtidal conditions for an intertidal species.

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.

Associations between DNA methylation and gene regulation depend on chromatin accessibility during transgenerational plasticity

BACKGROUND

Epigenetic processes are proposed to be a mechanism regulating gene expression during phenotypic plasticity. However, environmentally induced changes in DNA methylation exhibit little-to-no association with differential gene expression in metazoans at a transcriptome-wide level. It remains unexplored whether associations between environmentally induced differential methylation and expression are contingent upon other epigenomic processes such as chromatin accessibility. We quantified methylation and gene expression in larvae of the purple sea urchin Strongylocentrotus purpuratus exposed to different ecologically relevant conditions during gametogenesis (maternal conditioning) and modeled changes in gene expression and splicing resulting from maternal conditioning as functions of differential methylation, incorporating covariates for genomic features and chromatin accessibility. We detected significant interactions between differential methylation, chromatin accessibility, and genic feature type associated with differential expression and splicing.

RESULTS

Differential gene body methylation had significantly stronger effects on expression among genes with poorly accessible transcriptional start sites while baseline transcript abundance influenced the direction of this effect. Transcriptional responses to maternal conditioning were 4-13 × more likely when accounting for interactions between methylation and chromatin accessibility, demonstrating that the relationship between differential methylation and gene regulation is partially explained by chromatin state.

CONCLUSIONS

DNA methylation likely possesses multiple associations with gene regulation during transgenerational plasticity in S. purpuratus and potentially other metazoans, but its effects are dependent on chromatin accessibility and underlying genic features.

Epigenetic analytical approaches in ecotoxicological aquatic research

Environmental epigenetics has become a key research focus in global climate change studies and environmental pollutant investigations impacting aquatic ecosystems. Specifically, triggered by environmental stress conditions, intergenerational DNA methylation changes contribute to biological adaptive responses and survival of organisms to increase their tolerance towards these conditions. To critically review epigenetic analytical approaches in ecotoxicological aquatic research, we evaluated 78 publications reported over the past five years (2016-2021) that applied these methods to investigate the responses of aquatic organisms to environmental changes and pollution. The results show that DNA methylation appears to be the most robust epigenetic regulatory mark studied in aquatic animals. As such, multiple DNA methylation analysis methods have been developed in aquatic organisms, including enzyme restriction digestion-based and methyl-specific immunoprecipitation methods, and bisulfite (in)dependent sequencing strategies. In contrast, only a handful of aquatic studies, i.e. about 15%, have been focusing on histone variants and post-translational modifications due to the lack of species-specific affinity based immunological reagents, such as specific antibodies for chromatin immunoprecipitation applications. Similarly, ncRNA regulation remains as the least popular method used in the field of environmental epigenetics. Insights into the opportunities and challenges of the DNA methylation and histone variant analysis methods as well as decreasing costs of next generation sequencing approaches suggest that large-scale epigenetic environmental studies in model and non-model organisms will soon become available in the near future. Moreover, antibody-dependent and independent methods, such as mass spectrometry-based methods, can be used as an alternative epigenetic approach to characterize global changes of chromatin histone modifications in future aquatic research. Finally, a systematic guide for DNA methylation and histone variant methods is offered for ecotoxicological aquatic researchers to select the most relevant epigenetic analytical approach in their research.

Functional Niche Partitioning Occurs over Body Size but Not Nutrient Reserves nor Melanism in a Polar Carabid Beetle along an Altitudinal Gradient

Phenotypic plasticity can favor the emergence of different morphotypes specialized in specific ranges of environmental conditions. The existence of intraspecific partitioning confers resilience at the species scale and can ultimately determine species survival in a context of global changes. Amblystogenium pacificum is a carabid beetle endemic to the sub-Antarctic Crozet Islands, and it has two distinctive morphotypes based on body coloration. For this study, A. pacificum specimens of functional niches were sampled along an altitudinal gradient (as a proxy for temperature), and some morphological and biochemical traits were measured. We used an FAMD multivariate analysis and linear mixed-effects models to test whether these traits were related to morphotype, altitude, and sexual dimorphism. We then calculated and compared the functional niches at different altitudes and tested for niche partitioning through a hypervolume approach. We found a positive hump-shaped correlation between altitude and body size as well as higher protein and sugar reserves in females than in males. Our functional hypervolume results suggest that the main driver of niche partitioning along the altitudinal gradient is body size rather than morphotype or sex, even though darker morphotypes tended to be more functionally constrained at higher altitudes and females showed limited trait variations at the highest altitude.

Genome-wide sequencing identifies a thermal-tolerance related synonymous mutation in the mussel, Mytilisepta virgata

The roles of synonymous mutations for adapting to stressful thermal environments are of fundamental biological and ecological interests but poorly understood. To study whether synonymous mutations influence thermal adaptation at specific microhabitats, a genome-wide genotype-phenotype association analysis is carried out in the black mussels Mytilisepta virgata. A synonymous mutation of Ubiquitin-specific Peptidase 15 (MvUSP15) is significantly associated with the physiological upper thermal limit. The individuals carrying GG genotype (the G-type) at the mutant locus possess significantly lower heat tolerance compared to the individuals carrying GA and AA genotypes (the A-type). When heated to sublethal temperature, the G-type exhibit higher inter-individual variations in MvUSP15 expression, especially for the mussels on the sun-exposed microhabitats. Taken together, a synonymous mutation in MvUSP15 can affect the gene expression profile and interact with microhabitat heterogeneity to influence thermal resistance. This integrative study sheds light on the ecological importance of adaptive synonymous mutations as an underappreciated genetic buffer against heat stress and emphasizes the importance of integrative studies at a microhabitat scale for evaluating and predicting the impacts of climate change.

Pollution gradient leads to local adaptation and small-scale spatial variability of communities and functions in an urban marine environment

Urbanization of coastal habitats, of which harbors and marinas are the paragon, has led to various ecological paradigms about their functioning. Harbor infrastructures offer new hard substrata that are colonized by a wide variety of organisms (biofouling) including many introduced species. These structures also modify hydrodynamism and contaminant dispersal, leading to strong disturbance gradients within them. Differences in sessile community structure have previously been correlated to these gradients at small spatial scale (<100 m). Local adaptation might be involved to explain such results, but as correlation is not causation, the present study aims to understand the causal link between the environmental gradients and community structure through a reciprocal transplant experiment among three sites of a marina (inner, middle, entrance). Our results highlighted strong small-scale spatial variations of contaminants (trace metals, PCB, pesticides, and PAH) in sediments and animal samples which have been causally linked to changes in community composition after transplant. But historical contingency and colonization succession also play an important role. Our results provided strong evidence for local adaptation since community structure, respiration, and pollutant uptake in Bugula neritina, as well as the metabolomes of B. neritina and Ciona intestinalis were impacted by the transplant with a disadvantage for individuals transplanted from the entrance to the inner location. The here observed results may thus indicate that the disturbance gradient in marinas might constitute a staple for selecting pollutant-resistant species and populations, causing local adaptation. This highlights the importance of conducting further studies into small scale local adaptation.

Interactions Between Temperature Variability and Reproductive Physiology Across Traits in an Intertidal Crab

Thermal extremes alter population processes, which can result in part from temperature-induced movement at different spatial and temporal scales. Thermal thresholds for animal movement likely change based on underlying thermal physiology and life-history stage, a topic that requires greater study. The intertidal porcelain crab Petrolisthes cinctipes currently experiences temperatures that can reach near-lethal levels in the high-intertidal zone at low tide. However, the thermal thresholds that trigger migration to cooler microhabitats, and the extent to which crabs move in response to temperature, remain unknown. Moreover, the influence of reproductive status on these thresholds is rarely investigated. We integrated demographic, molecular, behavioral, and physiological measurements to determine if behavioral thermal limits varied due to reproductive state. Demographic data showed a trend for gravid, egg bearing, crabs to appear more often under rocks in the cooler intertidal zone where crab density is highest. In situ expression of 31 genes related to stress, metabolism, and growth in the field differed significantly based on intertidal elevation, with mid-intertidal crabs expressing the gene for the reproductive yolk protein vitellogenin (vg) earlier in the season. Furthermore, VG protein levels were shown to increase with density for female hemolymph. Testing for temperatures that elicit movement revealed that gravid females engage in heat avoidance behavior at lower temperatures (i.e., have a lower voluntary thermal maximum, VTmax) than non-gravid females. VTmax was positively correlated with the temperature of peak firing rate for distal afferent nerve fibers in the walking leg, a physiological relationship that could correspond to the mechanistic underpinning for temperature dependent movement. The vulnerability of marine organisms to global change is predicated by their ability to utilize and integrate physiological and behavioral strategies in response to temperature to maximize survival and reproduction. Interactions between fine-scale temperature variation and reproductive biology can have important consequences for the ecology of species, and is likely to influence how populations respond to ongoing climate change.

Can heat shock protein 70 (HSP70) serve as biomarkers in Antarctica for future ocean acidification, warming and salinity stress?

The Antarctic Peninsula is one of the fastest-warming places on Earth. Elevated sea water temperatures cause glacier and sea ice melting. When icebergs melt into the ocean, it “freshens” the saltwater around them, reducing its salinity. The oceans absorb excess anthropogenic carbon dioxide (CO2) causing decline in ocean pH, a process known as ocean acidification. Many marine organisms are specifically affected by ocean warming, freshening and acidification. Due to the sensitivity of Antarctica to global warming, using biomarkers is the best way for scientists to predict more accurately future climate change and provide useful information or ecological risk assessments. The 70-kilodalton (kDa) heat shock protein (HSP70) chaperones have been used as biomarkers of stress in temperate and tropical environments. The induction of the HSP70 genes (Hsp70) that alter intracellular proteins in living organisms is a signal triggered by environmental temperature changes. Induction of Hsp70 has been observed both in eukaryotes and in prokaryotes as response to environmental stressors including increased and decreased temperature, salinity, pH and the combined effects of changes in temperature, acidification and salinity stress. Generally, HSP70s play critical roles in numerous complex processes of metabolism; their synthesis can usually be increased or decreased during stressful conditions. However, there is a question as to whether HSP70s may serve as excellent biomarkers in the Antarctic considering the long residence time of Antarctic organisms in a cold polar environment which appears to have greatly modified the response of heat responding transcriptional systems. This review provides insight into the vital roles of HSP70 that make them ideal candidates as biomarkers for identifying resistance and resilience in response to abiotic stressors associated with climate change, which are the effects of ocean warming, freshening and acidification in Antarctic organisms.

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.

Parallel shape divergence between ecotypes of the limpet Nacella concinna along the Antarctic Peninsula: a new model species for parallel evolution?

Parallel phenotypic divergence is the independent differentiation between phenotypes of the same lineage or species occupying ecologically similar environments in different populations. We tested in the Antarctic limpet Nacella concinna the extent of parallel morphological divergence in littoral and sublittoral ecotypes throughout its distribution range. These ecotypes differ in morphological, behavioural and physiological characteristics. We studied the lateral and dorsal outlines of shells and the genetic variation of the mitochondrial gene Cytochrome Oxidase subunit I from both ecotypes in 17 sample sites along more than 2,000 km. The genetic data indicate that both ecotypes belong to a single evolutionary lineage. The magnitude and direction of phenotypic variation differ between ecotypes across sample sites; completely parallel ecotype-pairs (i.e., they diverge in the same magnitude and in the same direction) were detected in 84.85% of lateral and 65.15% in dorsal view comparisons. Besides, specific traits (relative shell height, position of shell apex, and elliptical/pear-shape outline variation) showed high parallelism. We observed weak morphological covariation between the two shape shell views, indicating that distinct evolutionary forces and environmental pressures could be acting on this limpet shell shape. Our results demonstrate there is a strong parallel morphological divergence pattern in N. concinna along its distribution, making this Antarctic species a suitable model for the study of different evolutionary forces shaping the shell evolution of this limpet.

Large within, and between, species differences in marine cellular responses: Unpredictability in a changing environment

Predicting the impacts of altered environments on future biodiversity requires a detailed understanding of organism responses to change. To date, studies evaluating mechanisms underlying marine organism stress responses have largely concentrated on oxygen limitation and the use of heat shock proteins as biomarkers. However, whether these biomarkers represent responses that are consistent across species and different environmental stressors remains open to question. Here we show that responses to four different thermal stresses (three rates of thermal ramping (1 °C h^-1, 1 °C day^-1 or 1 °C 3 day^-1) and a three-month acclimation to warming of 2 °C) applied to three species of Antarctic marine invertebrate produced highly individual responses in gene expression profiles, both within and between species. Mapping the gene expression profiles from each treatment for each of the three species, identified considerable difference in numbers of differentially regulated transcripts ranging from 10 to 3011. When these data were correlated across the different temperature treatments, there was no evidence for a common response with only 0-2 transcripts shared between all four treatments within any one species. There were also no shared differentially expressed genes across species, even at the same thermal ramping rates. The classical cellular stress response (CSR) i.e. up-regulation of heat shock proteins, was only strongly present in two species at the fastest ramping rate of 1 °C h^-1, albeit with different sets of stress genes expressed in each species. These data demonstrate the wide variability in response to warming at the molecular level in marine species. Therefore, identification of biodiversity stress responses engendered by changing conditions will require evaluation at the species level using targeted key members of the ecosystem, strongly correlated to the local biotic and abiotic factors.

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.

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.

Ocean acidification effects on the stress response in a calcifying antarctic coastal organism: The case of Nacella concinna ecotypes

Ocean acidification (OA) could become a serious threat for the Antarctic marine ecosystem over coming years, as the solubility of atmospheric CO₂ and CaCO₃ minerals increases at lower temperatures. We evaluated the effect of OA on the stress response of the limpet Nacella concinna by measuring gene expression levels. The experiment was performed with the two ecotypes (Littoral and Sublittoral) of the species during 54 days (IPCC, 2019 scenario RCP8.5; control, ~375 ppm; low-pH treatment, ~923 ppm). Exposure to low-pH treatment during 15 days triggered the down-regulation of two heat-shock protein genes (HSP70A, HSP70B) only in sublittoral individuals. Little variation in the relative expression values of all genes in both ecotypes was observed probably, due to a historical exposure to the substantial daily natural pH fluctuations recorded in the study area during the experiment. This study provides relevant baseline data for future OA experiments on coastal species in Antarctica.

How do life-history traits influence the fate of intertidal and subtidal Mytilus galloprovincialis in a changing climate?

Coastal organisms (i.e. intertidal or upper subtidal species) live in between the terrestrial and aquatic realms, making them particularly vulnerable to climate change. In this context, intertidal organisms may suffer from the predicted sea level rise (increasing their submerged time) while subtidal organisms may suffer from anthropically-induced hypoxia and its consequences. Although there is some knowledge on how coastal organisms adapt to environmental changes, the biochemical and physiological consequences of prolonged submergence periods have not yet been well characterized. Thus, the present study aimed to assess the biochemical alterations experienced by intertidal organisms maintained always under tidal exposure (IT); intertidal organisms maintained submersed (IS); subtidal organisms maintained always submersed (SS); subtidal organisms under tidal exposure (ST). For this, Mytilus galloprovincialis specimens from contiguous intertidal and subtidal populations were exposed to the above mentioned conditions for twenty-eight days. Results indicated that both intertidal and subtidal mussels are adapted to the oxidative stress pressure caused by tidal and submerged conditions tested. Intertidal mussels did not seem to be negatively affected by submergence while ST specimens were energetically challenged by tidal exposure. Both IT and ST mussels consumed glycogen to fuel up mechanisms aiming to maintain redox homeostasis. Overall, both intertidal and subtidal populations were capable of coping with tidal exposure, although the strategies employed differed between them. These findings indicate that although IT mussels may not significantly suffer from the longer-term submergence, hypoxic events occurring in the context of global warming and other anthropogenic impacts may have consequences on both IT and ST populations. Altogether, it is important to highlight that tides may act as a confounding factor in experiments concerning coastal organisms, as it causes additional physiological and biochemical perturbations.

Calcium homeostasis and stable fatty acid composition underpin heatwave tolerance of the keystone polychaete Hediste diversicolor

Extreme weather events, such as heatwaves, are becoming increasingly frequent, long-lasting and severe as global climate change continues, shaping marine biodiversity patterns worldwide. Increased risk of overheating and mortality across major taxa have been recurrently observed, jeopardizing the sustainability of ecosystem services. Molecular responses of species, which scale up to physiological and population responses, are determinant processes that modulate species sensitivity or tolerance to extreme weather events. Here, by integrating proteomic, fatty acid profiling and physiological approaches, we show that the tolerance of the intertidal ragworm Hediste diversicolor, a keystone species in estuarine ecosystems and an emergent blue bio-resource, to long-lasting heatwaves (24  vs 30 °C for 30 days) is shaped by calcium homeostasis, immune function and stability of fatty acid profiles. These features potentially enabled H. diversicolor to increase its thermal tolerance limit by 0.81 °C under the heatwave scenario and maintain survival. No growth trade-offs were detected, as wet weight remained stable across conditions. Biological variation of physiological parameters was lower when compared to molecular measures. Proteins showed an overall elevated coefficient of variation, although decreasing molecular variance under the heatwave scenario was observed for both proteins and fatty acids. This finding is consistent with the phenomenon of physiological canalization in extreme environments and contradicts the theory that novel conditions increase trait variation. Our results show that keystone highly valued marine polychaetes are tolerant to heatwaves, confirming the potential of H. diversicolor as a blue bio-resource and opening new avenues for sustainable marine aquaculture development.

Transcriptional frontloading contributes to cross-tolerance between stressors

The adaptive value of phenotypic plasticity for performance under single stressors is well documented. However, plasticity may only truly be adaptive in the natural multifactorial environment if it confers resilience to stressors of a different nature, a phenomenon known as cross-tolerance. An understanding of the mechanistic basis of cross-tolerance is essential to aid prediction of species resilience to future environmental change. Here, we identified mechanisms underpinning cross-tolerance between two stressors predicted to increasingly challenge aquatic ecosystems under climate change, chronic warming and hypoxia, in an ecologically-important aquatic invertebrate. Warm acclimation improved hypoxic performance through an adaptive hypometabolic strategy and changes in the expression of hundreds of genes that are important in the response to hypoxia. These 'frontloaded' genes showed a reduced reaction to hypoxia in the warm acclimated compared to the cold acclimated group. Frontloaded genes included stress indicators, immune response and protein synthesis genes that are protective at the cellular level. We conclude that increased constitutive gene expression as a result of warm acclimation reduced the requirement for inducible stress responses to hypoxia. We propose that transcriptional frontloading contributes to cross-tolerance between stressors and may promote fitness of organisms in environments increasingly challenged by multiple anthropogenic threats.

Identifying conserved polychaete molecular markers of metal exposure: Comparative analyses using the Alitta virens (Annelida, Lophotrochozoa) transcriptome

Polychaetes are vital for evaluating the effects of toxic metals in marine systems, and sensitive molecular biomarkers should be integral to monitoring efforts. However, the few polychaete markers that exist are inconsistent, even within the same species, failing to identify gene expression changes in metal-exposed animals incurring clear metabolic costs. Comparing previously characterised polychaete metal-responsive genes with those of another carefully selected species could identify biomarkers applicable across polychaetes. The ragworm Alitta virens (Sars, 1835) is particularly suited for such comparisons due to its dominance of fully saline coastal areas, widespread distribution, large biomass, and its phylogenetic position relative to other polychaete 'omic' resources. A transcriptome atlas for A. virens was generated and an RNASeq-qPCR screening approach was used to characterise the response to chronic exposures of environmentally relevant concentrations of copper and zinc in controlled mesocosms. Genes presenting dramatic expression changes in A. virens were compared with known metal-responsive genes in other polychaetes to identify new possible biomarkers and assess those currently used. This revealed some current markers should probably be abandoned (e.g. Atox1), while others, such as GST-Omega, should be used with caution, as different polychaete species appear to upregulate distinct GST-Omega orthologues. In addition, the comparisons give some indication of genes that are induced by metal exposure across phylogenetically divergent polychaetes, including a suite of haemoglobin subunits and linker chains that could play conserved roles in metal-stress response. Although such newly identified markers need further characterisation, they offer alternatives to current markers that are plainly insufficient.

Evaluating the effects of ocean warming and freshening on the physiological energetics and transcriptomic response of the Antarctic limpet Nacella concinna

In the Southern Ocean, warming and freshening are expected to be prominent signals of climate change and the reduced ability of Antarctic marine organisms to cope with changing environmental conditions could challenge their future survival. The Antarctic limpet Nacella concinna is a macroinvertebrate of rocky ecosystems, which occurs in high densities in the shallow subtidal zone. Subtidal individuals were exposed to a combination of temperatures (1, 4, 8, 11, 14 °C) and salinities (20 and 30 psu) for a 60-day period. A drastic increment in mortality was observed with seawater warming, showing that N. concinna is highly stenothermal, with limited ability to survive at temperatures warmer than 4 °C, although there was some degree of acclimation at 4 °C and ambient salinity (30 psu). This study confirmed the stenohaline characteristic of this species, with mortality reaching 50% and lower scope for growth at low salinity (20 psu) even at the control temperature (1 °C). At the sub-cellular level, limpets' low tolerance to out-of range salinity is illustrated by the activation of cell remodelling processes whereas the down-regulation of chaperones proteins and plasma membrane ATPase suggest that under the combination of warming and freshening N. concinna experiences a severe level of stress and devote much of its energy to somatic maintenance and survival. The drastic effect observed can be explained by its subtidal origin, an environment with more stable conditions. The surviving individuals at 1 °C and lowered salinity (20 psu) were either more tolerant or showing signs of acclimation after 60 days, but the combination of warming and freshening have a greater combined stress. Projections of climate change for end of the century for this part of the Antarctic can, therefore, result in a significant diminution of the subtidal population of N. concinna, affecting ecological interactions and diversity of the food web.

Deciphering mollusc shell production: the roles of genetic mechanisms through to ecology, aquaculture and biomimetics

Most molluscs possess shells, constructed from a vast array of microstructures and architectures. The fully formed shell is composed of calcite or aragonite. These CaCO₃ crystals form complex biocomposites with proteins, which although typically less than 5% of total shell mass, play significant roles in determining shell microstructure. Despite much research effort, large knowledge gaps remain in how molluscs construct and maintain their shells, and how they produce such a great diversity of forms. Here we synthesize results on how shell shape, microstructure, composition and organic content vary among, and within, species in response to numerous biotic and abiotic factors. At the local level, temperature, food supply and predation cues significantly affect shell morphology, whilst salinity has a much stronger influence across latitudes. Moreover, we emphasize how advances in genomic technologies [e.g. restriction site-associated DNA sequencing (RAD-Seq) and epigenetics] allow detailed examinations of whether morphological changes result from phenotypic plasticity or genetic adaptation, or a combination of these. RAD-Seq has already identified single nucleotide polymorphisms associated with temperature and aquaculture practices, whilst epigenetic processes have been shown significantly to modify shell construction to local conditions in, for example, Antarctica and New Zealand. We also synthesize results on the costs of shell construction and explore how these affect energetic trade-offs in animal metabolism. The cellular costs are still debated, with CaCO₃ precipitation estimates ranging from 1-2 J/mg to 17-55 J/mg depending on experimental and environmental conditions. However, organic components are more expensive (~29 J/mg) and recent data indicate transmembrane calcium ion transporters can involve considerable costs. This review emphasizes the role that molecular analyses have played in demonstrating multiple evolutionary origins of biomineralization genes. Although these are characterized by lineage-specific proteins and unique combinations of co-opted genes, a small set of protein domains have been identified as a conserved biomineralization tool box. We further highlight the use of sequence data sets in providing candidate genes for in situ localization and protein function studies. The former has elucidated gene expression modularity in mantle tissue, improving understanding of the diversity of shell morphology synthesis. RNA interference (RNAi) and clustered regularly interspersed short palindromic repeats - CRISPR-associated protein 9 (CRISPR-Cas9) experiments have provided proof of concept for use in the functional investigation of mollusc gene sequences, showing for example that Pif (aragonite-binding) protein plays a significant role in structured nacre crystal growth and that the Lsdia1 gene sets shell chirality in Lymnaea stagnalis. Much research has focused on the impacts of ocean acidification on molluscs. Initial studies were predominantly pessimistic for future molluscan biodiversity. However, more sophisticated experiments incorporating selective breeding and multiple generations are identifying subtle effects and that variability within mollusc genomes has potential for adaption to future conditions. Furthermore, we highlight recent historical studies based on museum collections that demonstrate a greater resilience of molluscs to climate change compared with experimental data. The future of mollusc research lies not solely with ecological investigations into biodiversity, and this review synthesizes knowledge across disciplines to understand biomineralization. It spans research ranging from evolution and development, through predictions of biodiversity prospects and future-proofing of aquaculture to identifying new biomimetic opportunities and societal benefits from recycling shell products.

Molecular mechanisms of biomineralization in marine invertebrates

Much recent marine research has been directed towards understanding the effects of anthropogenic-induced environmental change on marine biodiversity, particularly for those animals with heavily calcified exoskeletons, such as corals, molluscs and urchins. This is because life in our oceans is becoming more challenging for these animals with changes in temperature, pH and salinity. In the future, it will be more energetically expensive to make marine skeletons and the increasingly corrosive conditions in seawater are expected to result in the dissolution of these external skeletons. However, initial predictions of wide-scale sensitivity are changing as we understand more about the mechanisms underpinning skeletal production (biomineralization). These studies demonstrate the complexity of calcification pathways and the cellular responses of animals to these altered conditions. Factors including parental conditioning, phenotypic plasticity and epigenetics can significantly impact the production of skeletons and thus future population success. This understanding is paralleled by an increase in our knowledge of the genes and proteins involved in biomineralization, particularly in some phyla, such as urchins, molluscs and corals. This Review will provide a broad overview of our current understanding of the factors affecting skeletal production in marine invertebrates. It will focus on the molecular mechanisms underpinning biomineralization and how knowledge of these processes affects experimental design and our ability to predict responses to climate change. Understanding marine biomineralization has many tangible benefits in our changing world, including improvements in conservation and aquaculture and exploitation of natural calcified structure design using biomimicry approaches that are aimed at producing novel biocomposites.

Combined stress of ocean acidification and warming influence survival and drives differential gene expression patterns in the Antarctic pteropod, Limacina helicina antarctica

The ecologically important thecosome pteropods in the Limacina spp. complex have recently been the focus of studies examining the impacts global change factors - e.g., ocean acidification (OA) and ocean warming (OW) - on their performance and physiology. This focus is driven by conservation concerns where the health of pteropod populations is threatened by the high susceptibility of their shells to dissolution in low aragonite saturation states associated with OA and how coupling of these stressors may push pteropods past the limits of physiological plasticity. In this manipulation experiment, we describe changes in the transcriptome of the Antarctic pteropod, Limacina helicina antarctica, to these combined stressors. The conditions used in the laboratory treatments met or exceeded those projected for the Southern Ocean by the year 2100. We made two general observations regarding the outcome of the data: (1) Temperature was more influential than pH in terms of changing patterns of gene expression, and (2) these Antarctic pteropods appeared to have a significant degree of transcriptomic plasticity to respond to acute abiotic stress in the laboratory. In general, differential gene expression was observed amongst the treatments; here, for example, transcripts associated with maintaining protein structure and cell proliferation were up-regulated. To disentangle the effects of OA and OW, we used a weighted gene co-expression network analysis to explore patterns of change in the transcriptome. This approach identified gene networks associated with OW that were enriched for transcripts proposed to be involved in increasing membrane fluidity at warmer temperatures. Together these data provide evidence that L.h.antarctica has a limited capacity to acclimate to the combined conditions of OA and OW used in this study. This reduced scope of acclimation argues for continued study of how adaptation to polar aquatic environments may limit the plasticity of present-day populations in responding to future environmental change.

Antarctic environmental change and biological responses

Antarctica and the surrounding Southern Ocean are facing complex environmental change. Their native biota has adapted to the region's extreme conditions over many millions of years. This unique biota is now challenged by environmental change and the direct impacts of human activity. The terrestrial biota is characterized by considerable physiological and ecological flexibility and is expected to show increases in productivity, population sizes and ranges of individual species, and community complexity. However, the establishment of non-native organisms in both terrestrial and marine ecosystems may present an even greater threat than climate change itself. In the marine environment, much more limited response flexibility means that even small levels of warming are threatening. Changing sea ice has large impacts on ecosystem processes, while ocean acidification and coastal freshening are expected to have major impacts.

Molecular assessment of wild populations in the marine realm: Importance of taxonomic, seasonal and habitat patterns in environmental monitoring

Scientists are currently faced with the challenge of assessing the effects of anthropogenic stressors on aquatic ecosystems. Cellular stress response (CSR) biomarkers are ubiquitous and phylogenetically conserved among metazoans and have been successfully applied in environmental monitoring but they can also vary according to natural biotic and abiotic factors. The reported variability may thus limit the wide application of biomarkers in monitoring, imposing the need to identify variability levels in the field. Our aim was to carry out a comprehensive in situ assessment of the CSR (heat shock protein 70 kDa, ubiquitin, antioxidant enzymes) and oxidative damage (lipid peroxidation) in wild populations across marine taxa by collecting fish, crustaceans, mollusks and cnidarians during two different seasons (spring and summer) and two habitat types (coast and estuary). CSR end-point patterns were different between taxa with mollusks having higher biomarker levels, followed by the cnidarians, while fish and crustaceans showed lower biomarker levels. The PCA showed clear clusters related to mobility/sessile traits with sessile organisms showing greater levels (>2-fold) of CSR proteins and oxidative damage. Mean intraspecific variability in the CSR measured by the coefficient of variation (% CV) (including data from all seasons and sites) was elevated (35-94%). Overall, there was a seasonal differentiation in biomarker patterns across taxonomic groups, especially evident in fish and cnidarians. A differentiation in biomarker patterns between habitat types was also observed and associated with phenotypic plasticity or local adaptation. Overall, specimens collected in the estuary had lower biomarker levels when compared to specimens collected in the coast. This work highlights the importance of assessing baseline biomarker levels across taxa, seasons and habitats prior to applying biomarker analyses in environmental monitoring. Selecting bioindicator species, defining sampling strategies, and identifying confounding factors are crucial preliminary steps that ensure the success of biomarkers as powerful tools in biomonitoring.

Chronic exposure to copper and zinc induces DNA damage in the polychaete Alitta virens and the implications for future toxicity of coastal sites

Copper and zinc are metals that have been traditionally thought of as past contamination legacies. However, their industrial use is still extensive and current applications (e.g. nanoparticles and antifouling paints) have become additional marine environment delivery routes. Determining a pollutant's genotoxicity is an ecotoxicological priority, but in marine benthic systems putative substances responsible for sediment genotoxicity have rarely been identified. Studies that use sediment as the delivery matrix combined with exposures over life-history relevant timescales are also missing for metals. Here we assess copper and zinc's genotoxicity by exposing the ecologically important polychaete Alitta virens to sediment spiked with environmentally relevant concentrations for 9 months. Target bioavailable sediment and subsequent porewater concentrations reflect the global contamination range for coasts, whilst tissue concentrations, although elevated, were comparable with other polychaetes. Survival generally reduced as concentrations increased, but monthly analyses show that growth was not significantly different between treatments. The differential treatment mortality may have enabled the surviving worms in the high concentration treatments to capture more food thus removing any concentration treatment effects for biomass. Using the alkaline comet assay we confirm that both metals via the sediment are genotoxic at concentrations routinely found in coastal regions and this is supported by elevated DNA damage in worms from field sites. However, combined with the growth data it also highlights the tolerance of A. virens to DNA damage. Finally, using long term (decadal) monitoring data we show stable or increasing sediment concentrations of these metals for many areas. This will potentially mean coastal sediment is a significant mutagenic hazard to the benthic community for decades to come. An urgent reappraisal of the current input sources for these 'old pollutants' is, therefore, required.