Transcriptomic resilience to global warming in the seagrass Zostera marina, a marine foundation species

Exploring priming strategies to improve stress resilience of Posidonia oceanica seedlings

Seagrasses' ability to store information after exposure to stress (i.e. stress memory) and to better respond to further stress (i.e. priming) have recently been observed, although the temporal persistence of the memory and the mechanisms for priming induction remain to be defined. Here, we explored three priming strategies in Posidonia oceanica seedlings, each inducing a different level of stress, for temperature and salinity. We investigated changes in morphometry, growth rate and biomass between primed and non-primed seedlings. The results showed similar behaviour of seedlings when exposed to an acute stress event, regardless of whether they had been primed or not and of the priming strategy received. This opens the debate on the level of stress necessary for inducing a priming status and the persistence of the stress memory in P. oceanica seedlings. Although no priming-induced stress resistance was observed, seedlings showed unexpectedly high resilience to extreme levels of both abiotic stressors.

Marine heatwaves recurrence aggravates thermal stress in the surfgrass Phyllospadix scouleri

The surfgrass Phyllospadix scouleri grows in highly productive meadows along the Pacific coast of North America. This region has experienced increasingly severe marine heatwaves (MHWs) in recent years. Our study evaluated the impact of consecutive MHWs, simulated in mesocosms, on essential ecophysiological features of P. scouleri. Overall, our findings show that the plants' overall physiological status has been progressively declining. Interestingly, the indicators of physiological stress in photosynthesis only showed up once the initial heat exposure stopped (i.e., during the recovery period). The warming caused increased oxidative damage and a decrease in nitrate uptake rates. However, the levels of non-structural carbohydrates and relative growth rates were not affected. Our findings emphasize the significance of incorporating recovery periods in this type of study as they expose delayed stress responses. Furthermore, experiencing consecutive intense MHWs can harm surfgrasses over time, compromising the health of their meadows and the services they offer to the ecosystem.

Adaptive responses of eelgrass (Zostera marina L.) to ocean warming and acidification

Ocean warming (OW) and ocean acidification (OA), driven by rapid global warming accelerating at unprecedented rates, are profoundly impacting the stability of seagrass ecosystems. Yet, our current understanding of the effects of OW and OA on seagrass remains constrained. Herein, we investigated the response of eelgrass (Zostera marina L.), a representative seagrass species, to OW and OA through comprehensive transcriptomic and metabolomic analyses. The results showed notable variations in plant performance under varying conditions: OW, OA, and OWA (a combination of both conditions). Specifically, under average oceanic temperature conditions for eelgrass growth over the past 20 years -from May to November-OA promoted the production of differentially expressed genes and metabolites associated with alanine, aspartate, and glutamate metabolism, as well as starch and sucrose metabolism. Under warming condition, eelgrass was resistant to OA by accelerating galactose metabolism, along with glycine, serine, and threonine metabolism, as well as the tricarboxylic acid (TCA) cycle. Under the combined OW and OA condition, eelgrass stimulated fructose and mannose metabolism, glycolysis, and carbon fixation, in addition to galactose metabolism and the TCA cycle to face the interplay. Our findings suggest that eelgrass exhibits adaptive capacity by inducing different metabolites and associated genes, primarily connected with carbon and nitrogen metabolism, in response to varying degrees of OW and OA. The data generated here support the exploration of mechanisms underlying seagrass responses to environmental fluctuations, which hold critical significance for the future conservation and management of these ecosystems.

Signs of local adaptation by genetic selection and isolation promoted by extreme temperature and salinity in the Mediterranean seagrass Posidonia oceanica

Adaptation to local conditions is known to occur in seagrasses; however, knowledge of the genetic basis underlying this phenomenon remains scarce. Here, we analysed Posidonia oceanica from six sites within and around the Stagnone di Marsala, a semi-enclosed coastal lagoon where salinity and temperature exceed the generally described tolerance thresholds of the species. Sea surface temperatures (SSTs) were measured and plant samples were collected for the assessment of morphology, flowering rate and for screening genome-wide polymorphisms using double digest restriction-site-associated DNA sequencing. Results demonstrated more extreme SSTs and salinity levels inside the lagoon than the outer lagoon regions. Morphological results showed significantly fewer and shorter leaves and reduced rhizome growth of P. oceanica from the inner lagoon and past flowering events were recorded only for a meadow farthest away from the lagoon. Using an array of 51,329 single nucleotide polymorphisms, we revealed a clear genetic structure among the study sites and confirmed the genetic isolation and high clonality of the innermost site. In all, 14 outlier loci were identified and annotated with several proteins including those relate to plant stress response, protein transport and regulators of plant-specific developmental events. Especially, five outlier loci showed maximum allele frequency at the innermost site, likely reflecting adaptation to the extreme temperature and salinity regimes, possibly due to the selection of more resistant genotypes and the progressive restriction of gene flow. Overall, this study helps us to disentangle the genetic basis of seagrass adaptation to local environmental conditions and may support future works on assisted evolution in seagrasses.

Ocean current patterns drive the worldwide colonization of eelgrass (Zostera marina)

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

Sources of variability in seagrass fatty acid profiles and the need of identifying reliable warming descriptors

Global warming is expected to have inexorable and profound effects on marine ecosystems, particularly in foundation species such as seagrasses. Identifying responses to warming and comparing populations across natural temperature gradients can inform how future warming will impact the structure and function of ecosystems. Here, we investigated how thermal environment, intra-shoot and spatial variability modulate biochemical responses of the Mediterranean seagrass Posidonia oceanica. Through a space-for-time substitution study, Fatty acid (FA) profiles on the second and fifth leaf of the shoots were quantified at eight sites in Sardinia along a natural sea surface temperature (SST) summer gradient (about 4 °C). Higher mean SST were related to a decrease in the leaf total fatty acid content (LTFA), a reduction in polyunsaturated fatty acids (PUFA), omega-3/omega-6 PUFA and PUFA/saturated fatty acids (SFA) ratios and an increase in SFA, monounsaturated fatty acids and carbon elongation index (CEI, C18:2 n-6/C16:2 n-6) ratio. Results also revealed that FA profiles were strongly influenced by leaf age, independently of SST and spatial variability within sites. Overall, this study evidenced that the sensitive response of P. oceanica FA profiles to intra-shoot and spatial variability must not be overlooked when considering their response to temperature changes.

Gene co-expression network analysis for the selection of candidate early warning indicators of heat and nutrient stress in Posidonia oceanica

The continuous worldwide seagrasses decline calls for immediate actions in order to preserve this precious marine ecosystem. The main stressors that have been linked with decline in seagrasses are 1) the increasing ocean temperature due to climate change and 2) the continuous inputs of nutrients (eutrophication) associated with coastal human activities. To avoid the loss of seagrass populations, an "early warning" system is needed. We used Weighed Gene Co-expression Network Analysis (WGCNA), a systems biology approach, to identify potential candidate genes that can provide an early warning signal of stress in the Mediterranean iconic seagrass Posidonia oceanica, anticipating plant mortality. Plants were collected from both eutrophic (EU) and oligotrophic (OL) environments and were exposed to thermal and nutrient stress in a dedicated mesocosm. By correlating the whole-genome gene expression after 2-weeks exposure with the shoot survival percentage after 5-weeks exposure to stressors, we were able to identify several transcripts that indicated an early activation of several biological processes (BP) including: protein metabolic process, RNA metabolic process, organonitrogen compound biosynthetic process, catabolic process and response to stimulus, which were shared among OL and EU plants and among leaf and shoot apical meristem (SAM), in response to excessive heat and nutrients. Our results suggest a more dynamic and specific response of the SAM compared to the leaf, especially the SAM from plants coming from a stressful environment appeared more dynamic than the SAM from a pristine environment. A vast list of potential molecular markers is also provided that can be used as targets to assess field samples.

Ten years of marine evolutionary biology-Challenges and achievements of a multidisciplinary research initiative

The Centre for Marine Evolutionary Biology (CeMEB) at the University of Gothenburg, Sweden, was established in 2008 through a 10-year research grant of 8.7 m€ to a team of senior researchers. Today, CeMEB members have contributed >500 scientific publications, 30 PhD theses and have organised 75 meetings and courses, including 18 three-day meetings and four conferences. What are the footprints of CeMEB, and how will the centre continue to play a national and international role as an important node of marine evolutionary research? In this perspective article, we first look back over the 10 years of CeMEB activities and briefly survey some of the many achievements of CeMEB. We furthermore compare the initial goals, as formulated in the grant application, with what has been achieved, and discuss challenges and milestones along the way. Finally, we bring forward some general lessons that can be learnt from a research funding of this type, and we also look ahead, discussing how CeMEB's achievements and lessons can be used as a springboard to the future of marine evolutionary biology.

DNA methylation dynamics in a coastal foundation seagrass species under abiotic stressors

DNA methylation (DNAm) has been intensively studied in terrestrial plants in response to environmental changes, but its dynamic changes in a temporal scale remain unexplored in marine plants. The seagrass Posidonia oceanica ranks among the slowest-growing and longest-living plants on Earth, and is particularly vulnerable to sea warming and local anthropogenic pressures. Here, we analysed the dynamics of DNAm changes in plants collected from coastal areas differentially impacted by eutrophication (i.e. oligotrophic, Ol; eutrophic, Eu) and exposed to abiotic stressors (nutrients, temperature increase and their combination). Levels of global DNAm (% 5-mC) and the expression of key genes involved in DNAm were assessed after one, two and five weeks of exposure. Results revealed a clear differentiation between plants, depending on environmental stimuli, time of exposure and plants' origin. % 5-mC levels were higher during the initial stress exposure especially in Ol plants, which upregulated almost all genes involved in DNAm. Contrarily, Eu plants showed lower expression levels, which increased under chronic exposure to stressors, particularly to temperature. These findings show that DNAm is dynamic in P. oceanica during stress exposure and underlined that environmental epigenetic variations could be implicated in the regulation of acclimation and phenotypic differences depending on local conditions.

Advances in understanding multilevel responses of seagrasses to hypersalinity

Human- and nature-induced hypersaline conditions in coastal systems can lead to profound alterations of the structure and vitality of seagrass meadows and their socio-ecological benefits. In the last two decades, recent research efforts (>50 publications) have contributed significantly to unravel the physiological basis underlying the seagrass-hypersalinity interactions, although most (∼70%) are limited to few species (e.g. Posidonia oceanica, Zostera marina, Thalassia testudinum, Cymodocea nodosa). Variables related to photosynthesis and carbon metabolism are among the most prevalent in the literature, although other key metabolic processes such as plant water relations and responses at molecular (i.e. gene expression) and ultrastructure level are attracting attention. This review emphasises all these latest insights, offering an integrative perspective on the interplay among biological responses across different functional levels (from molecular to clonal structure), and their interaction with biotic/abiotic factors including those related to climate change. Other issues such as the role of salinity in driving the evolutionary trajectory of seagrasses, their acclimation mechanisms to withstand salinity increases or even the adaptive properties of populations that have historically lived under hypersaline conditions are also included. The pivotal role of the costs and limits of phenotypic plasticity in the successful acclimation of marine plants to hypersalinity is also discussed. Finally, some lines of research are proposed to fill the remaining knowledge gaps.

Single-molecule real-time sequencing of the full-length transcriptome of Halophila beccarii

Ecologically, Halophila beccarii Asch. is considered as a colonizing or a pioneer seagrass species and a “tiny but mighty” seagrass species, since it may recover quickly from disturbance generally. The use of transcriptome technology can provide a better understanding of the physiological processes of seagrasses. To date, little is known about the genome and transcriptome information of H. beccarii. In this study, we used single molecule real-time (SMRT) sequencing to obtain full-length transcriptome data and characterize the transcriptome structure. A total of 11,773 of the 15,348 transcripts were successfully annotated in seven databases. In addition, 1573 long non-coding RNAs, 8402 simple sequence repeats and 2567 transcription factors were predicted in all the transcripts. A GO analysis showed that 5843 transcripts were divided into three categories, including biological process (BP), cellular component (CC) and molecular function (MF). In these three categories, metabolic process (1603 transcripts), protein-containing complex (515 transcripts) and binding (3233 transcripts) were the primary terms in BP, CC, and MF, respectively. The major types of transcription factors were involved in MYB-related and NF-YB families. To the best of our knowledge, this is the first report of the transcriptome of H. beccarii using SMRT sequencing technology.

De novo assembly and annotation of the transcriptome of the endangered seagrass Zostera capensis: Insights from differential gene expression under thermal stress

Seagrasses are important marine ecosystem engineers but anthropogenic impacts and climate change have led to numerous population declines globally. In South Africa, Zostera capensis is endangered due to fragmented populations and heavy anthropogenic pressures on estuarine ecosystems that house the core of the populations. Addressing questions of how pressures such as climate change affect foundational species, including Z. capensis are crucial to supporting their conservation and underpin restoration efforts. Here we use ecological transcriptomics to study key functional responses of Z. capensis through quantification of gene expression after thermal stress and present the first reference transcriptome of Z. capensis. Four de novo reference assemblies (Trinity, IDBA-tran, RNAspades, SOAPdenovo) filtered through the EvidentialGene pipeline resulted in 153,755 transcripts with a BUSCO score of 66.1% for completeness. Differential expression analysis between heat stressed (32 °C for three days) and pre-warming plants identified genes involved in photosynthesis, oxidative stress, translation, metabolic and biosynthetic processes in the Z. capensis thermal stress response. This reference transcriptome is a significant contribution to the limited available genomic resources for Z. capensis and represents a vital tool for addressing questions around the species restoration and potential functional responses to warming marine environments.

A Transcriptomic Analysis of Phenotypic Plasticity in Crassostrea virginica Larvae under Experimental Acidification

Ocean acidification (OA) is a major threat to marine calcifiers, and little is known regarding acclimation to OA in bivalves. This study combined physiological assays with next-generation sequencing to assess the potential for recovery from and acclimation to OA in the eastern oyster (Crassostrea virginica) and identify molecular mechanisms associated with resilience. In a reciprocal transplant experiment, larvae transplanted from elevated pCO₂ (~1400 ppm) to ambient pCO₂ (~350 ppm) demonstrated significantly lower mortality and larger size post-transplant than oysters remaining under elevated pCO₂ and had similar mortality compared to those remaining in ambient conditions. The recovery after transplantation to ambient conditions demonstrates the ability for larvae to rebound and suggests phenotypic plasticity and acclimation. Transcriptomic analysis supported this hypothesis as genes were differentially regulated under OA stress. Transcriptomic profiles of transplanted and non-transplanted larvae terminating in the same final pCO₂ converged, further supporting the idea that acclimation underlies resilience. The functions of differentially expressed genes included cell differentiation, development, biomineralization, ion exchange, and immunity. Results suggest acclimation as a mode of resilience to OA. In addition, the identification of genes associated with resilience can serve as a valuable resource for the aquaculture industry, as these could enable marker-assisted selection of OA-resilient stocks.

Warming northward shifting southern limits of the iconic temperate seagrass (Zostera marina)

Global warming can shift the range edges of numerous species poleward. Here, eelgrass distribution was reinvestigated at its southern limits on the eastern coast of China, which indicated that there has been a northward shift in the southern limit of Z. marina. To determine if regional warming resulted in a northward shift in suitable eelgrass habitats, sixteen transplantations of adult eelgrass shoots and seeds at the historical southern distribution limit of eelgrass were conducted between 2016 and 2021. The results showed that high water temperatures in summer had negative effects on eelgrass growth, and directly triggered shoot mortality during 2016–2018. Under heat stress, antioxidant enzyme activity was initially increased, but then decreased under more stressful heat conditions; and the HSP70 protein and its molecular chaperone protein were highly expressed under heat stress. These results demonstrated that suitable eelgrass habitat was now located further north along the eastern coast of China.

•

High temperatures trigger seagrass (Zostera marina L.) restoration failure

•

None seedlings and adult shoots survived the first or second summer

•

Over-summering shoots with lower density, height, and rhizome diameter

•

Warming northward shifting eelgrass habitat range along the eastern coast of China

Local environment modulates whole-transcriptome expression in the seagrass Posidonia oceanica under warming and nutrients excess

The intensification of anomalous events of seawater warming and the co-occurrence with local anthropogenic stressors are threatening coastal marine habitats, including seagrasses, which form extensive underwater meadows. Eutrophication highly affects coastal environments, potentially summing up to the widespread effects of global climate changes. In the present study, we investigated for the first time in seagrasses, the transcriptional response of different plant organs (i.e., leaf and shoot apical meristem, SAM) of the Mediterranean seagrass Posidonia oceanica growing in environments with a different history of nutrient enrichment. To this end, a mesocosm experiment exposing plants to single (nutrient enrichment or temperature increase) and multiple stressors (nutrient enrichment plus temperature increase), was performed. Results revealed a differential transcriptome regulation of plants under single and multiple stressors, showing an organ-specific sensitivity depending on plants' origin. While leaf tissues were more responsive to nutrient stress, SAM revealed a higher sensitivity to temperature treatments, especially in plants already impacted in their native environment. The exposure to stress conditions induced the modulation of different biological processes. Plants living in an oligotrophic environment were more responsive to nutrients compared to plants from a eutrophic environment. Evidences that epigenetic mechanisms were involved in the regulation of transcriptional reprogramming were also observed in both plants' organs. These results represent a further step in the comprehension of seagrass response to abiotic stressors pointing out the importance of local pressures in a global warming scenario.

2b-RAD Genotyping of the Seagrass Cymodocea nodosa Along a Latitudinal Cline Identifies Candidate Genes for Environmental Adaptation

Plant populations distributed along broad latitudinal gradients often show patterns of clinal variation in genotype and phenotype. Differences in photoperiod and temperature cues across latitudes influence major phenological events, such as timing of flowering or seed dormancy. Here, we used an array of 4,941 SNPs derived from 2b-RAD genotyping to characterize population differentiation and levels of genetic and genotypic diversity of three populations of the seagrass Cymodocea nodosa along a latitudinal gradient extending across the Atlantic-Mediterranean boundary (i.e., Gran Canaria—Canary Islands, Faro—Portugal, and Ebro Delta—Spain). Our main goal was to search for potential outlier loci that could underlie adaptive differentiation of populations across the latitudinal distribution of the species. We hypothesized that such polymorphisms could be related to variation in photoperiod-temperature regime occurring across latitudes. The three populations were clearly differentiated and exhibited diverse levels of clonality and genetic diversity. Cymodocea nodosa from the Mediterranean displayed the highest genotypic richness, while the Portuguese population had the highest clonality values. Gran Canaria exhibited the lowest genetic diversity (as observed heterozygosity). Nine SNPs were reliably identified as outliers across the three sites by two different methods (i.e., BayeScan and pcadapt), and three SNPs could be associated to specific protein-coding genes by screening available C. nodosa transcriptomes. Two SNPs-carrying contigs encoded for transcription factors, while the other one encoded for an enzyme specifically involved in the regulation of flowering time, namely Lysine-specific histone demethylase 1 homolog 2. When analyzing biological processes enriched within the whole dataset of outlier SNPs identified by at least one method, “regulation of transcription” and “signalling” were among the most represented. Our results highlight the fundamental importance signal integration and gene-regulatory networks, as well as epigenetic regulation via DNA (de)methylation, could have for enabling adaptation of seagrass populations along environmental gradients.

Physiological and morphological effects of a marine heatwave on the seagrass Cymodocea nodosa

Marine heatwaves (MHWs) are increasing in frequency and intensity as part of climate change, yet their impact on seagrass is poorly known. The present work evaluated the physiological and morphological responses of Cymodocea nodosa to a MHW. C. nodosa shoots were transplanted into a mesocosm facility. To simulate a MHW, water temperature was raised from 20 to 28 °C, kept 7 days at 28 °C, cooled down back to 20 °C and then maintained at 20 °C during an 8-day recovery period. The potentially stressful effects of the simulated heatwave on the photosynthetic performance, antioxidative-stress level and area vs dry weight ratio of leaves were investigated. The maximum quantum yield of photosystem II (ΦPSII) increased during the heatwave, allowing the plants to maintain their photosynthetic activity at control level. Negative effects on the photosynthetic performance and leaf biomass of C. nodosa were observed during the recovery period. No significant oxidative stress was observed throughout the experiment. Overall, although C. nodosa showed a relative tolerance to MHWs compared to other species, its population in Ria Formosa is likely to be negatively affected by the forecasted climate change scenarios.

Improved chromosome-level genome assembly and annotation of the seagrass, Zostera marina (eelgrass)

Background: Seagrasses (Alismatales) are the only fully marine angiosperms.  Zostera marina (eelgrass) plays a crucial role in the functioning of coastal marine ecosystems and global carbon sequestration. It is the most widely studied seagrass and has become a marine model system for exploring adaptation under rapid climate change. The original draft genome (v.1.0) of the seagrass  Z. marina (L.) was based on a combination of Illumina mate-pair libraries and fosmid-ends. A total of 25.55 Gb of Illumina and 0.14 Gb of Sanger sequence was obtained representing 47.7× genomic coverage. The assembly resulted in ~2000 unordered scaffolds (L50 of 486 Kb), a final genome assembly size of 203MB, 20,450 protein coding genes and 63% TE content. Here, we present an upgraded chromosome-scale genome assembly and compare v.1.0 and the new v.3.1, reconfirming previous results from Olsen et al. (2016), as well as pointing out new findings.  

Methods: The same high molecular weight DNA used in the original sequencing of the Finnish clone was used. A high-quality reference genome was assembled with the MECAT assembly pipeline combining PacBio long-read sequencing and Hi-C scaffolding. 

Results: In total, 75.97 Gb PacBio data was produced. The final assembly comprises six pseudo-chromosomes and 304 unanchored scaffolds with a total length of 260.5Mb and an N50 of 34.6 MB, showing high contiguity and few gaps (~0.5%). 21,483 protein-encoding genes are annotated in this assembly, of which 20,665 (96.2%) obtained at least one functional assignment based on similarity to known proteins. 

Conclusions: As an important marine angiosperm, the improved  Z. marina genome assembly will further assist evolutionary, ecological, and comparative genomics at the chromosome level. The new genome assembly will further our understanding into the structural and physiological adaptations from land to marine life.

Limited plasticity in thermally tolerant ectotherm populations: evidence for a trade-off

Many species face extinction risks owing to climate change, and there is an urgent need to identify which species' populations will be most vulnerable. Plasticity in heat tolerance, which includes acclimation or hardening, occurs when prior exposure to a warmer temperature changes an organism's upper thermal limit. The capacity for thermal acclimation could provide protection against warming, but prior work has found few generalizable patterns to explain variation in this trait. Here, we report the results of, to our knowledge, the first meta-analysis to examine within-species variation in thermal plasticity, using results from 20 studies (19 species) that quantified thermal acclimation capacities across 78 populations. We used meta-regression to evaluate two leading hypotheses. The climate variability hypothesis predicts that populations from more thermally variable habitats will have greater plasticity, while the trade-off hypothesis predicts that populations with the lowest heat tolerance will have the greatest plasticity. Our analysis indicates strong support for the trade-off hypothesis because populations with greater thermal tolerance had reduced plasticity. These results advance our understanding of variation in populations' susceptibility to climate change and imply that populations with the highest thermal tolerance may have limited phenotypic plasticity to adjust to ongoing climate warming.

Overexpression of seagrass DnaJ gene ZjDjB1 enhances the thermotolerance of transgenic arabidopsis thaliana

Seagrass meadows are one of the most important marine resources that grow along the coast. They provide habitat and a food source for animals. They also protect the coast, fix sediment and purify seawater. In the current period of global climate change, anomalies in coastal water temperatures are increasing. A sudden increase in water temperature owing to a heat wave can have a profound effect on seagrass. Zostera japonica is a type of intertidal seagrasses, which is exposed to the air at low tide. High temperatures in the summer often lead to a decline in seagrass meadows. DnaJ proteins, also known as J proteins, are a family of conserved chaperone proteins. They are designated as J proteins because they contain a highly conserved J domain. They function as chaperones of heat shock proteins in organisms. In this study, the role of DnaJ protein (ZjDjB1) of Z. japonica under heat stress was studied. ZjDjB1 was localized to the cytoplasm and nucleus. The overexpression of ZjDjB1 in Arabidopsis thaliana results in an increase in thermotolerance and a decrease in the accumulation of reactive oxygen species and also a reduction in membrane damage. ZjDjB1 may achieve this goal by maintaining a low activity of proteolytic enzymes.

Sensitivity of Photosynthesis to Warming in Two Similar Species of the Aquatic Angiosperm Ruppia from Tropical and Temperate Habitats

Climate change-related events, such as marine heatwaves, are increasing seawater temperatures, thereby putting pressure on marine biota. The cosmopolitan distribution and significant contribution to marine primary production by the genus Ruppia makes them interesting organisms to study thermal tolerance and local adaptation. In this study, we investigated the photosynthetic responses in Ruppia to the predicted future warming in two contrasting bioregions, temperate Sweden and tropical Thailand. Through DNA barcoding, specimens were determined to Ruppia cirrhosa for Sweden and Ruppia maritima for Thailand. Photosynthetic responses were assessed using pulse amplitude-modulated fluorometry, firstly in short time incubations at 18, 23, 28, and 33 °C in the Swedish set-up and 28, 33, 38, and 43 °C in the Thai set-up. Subsequent experiments were conducted to compare the short time effects to longer, five-day incubations in 28 °C for Swedish plants and 40 °C for Thai plants. Swedish R. cirrhosa displayed minor response, while Thai R. maritima was more sensitive to both direct and prolonged temperature stress with a drastic decrease in the photosynthetic parameters leading to mortality. The results indicate that in predicted warming scenarios, Swedish R. cirrhosa may sustain an efficient photosynthesis and potentially outcompete more heat-sensitive species. However, populations of the similar R. maritima in tropical environments may suffer a decline as their productivity will be highly reduced.

Gene body DNA methylation in seagrasses: inter- and intraspecific differences and interaction with transcriptome plasticity under heat stress

The role of DNA methylation and its interaction with gene expression and transcriptome plasticity is poorly understood, and current insight comes mainly from studies in very few model plant species. Here, we study gene body DNA methylation (gbM) and gene expression patterns in ecotypes from contrasting thermal environments of two marine plants with contrasting life history strategies in order to explore the potential role epigenetic mechanisms could play in gene plasticity and responsiveness to heat stress. In silico transcriptome analysis of CpGO/E ratios suggested that the bulk of Posidonia oceanica and Cymodocea nodosa genes possess high levels of intragenic methylation. We also observed a correlation between gbM and gene expression flexibility: genes with low DNA methylation tend to show flexible gene expression and plasticity under changing conditions. Furthermore, the empirical determination of global DNA methylation (5-mC) showed patterns of intra and inter-specific divergence that suggests a link between methylation level and the plants' latitude of origin and life history. Although we cannot discern whether gbM regulates gene expression or vice versa, or if other molecular mechanisms play a role in facilitating transcriptome responsiveness, our findings point to the existence of a relationship between gene responsiveness and gbM patterns in marine plants.

Intraspecific genetic variation matters when predicting seagrass distribution under climate change

Seagrasses play a vital role in structuring coastal marine ecosystems, but their distributional range and genetic diversity have declined rapidly in recent decades. To improve conservation of seagrass species, it is important to predict how climate change may impact their ranges. Such predictions are typically made with correlative species distribution models (SDMs), which can estimate a species' potential distribution under present and future climatic scenarios given species' presence data and climatic predictor variables. However, these models are typically constructed with species-level data, and thus ignore intraspecific genetic variability, which can give rise to populations with adaptations to heterogeneous climatic conditions. Here, we explore the link between intraspecific adaptation and niche differentiation in Thalassia hemprichii, a seagrass broadly distributed in the tropical Indo-Pacific Ocean and a crucial provider of habitat for numerous marine species. By retrieving and re-analysing microsatellite data from previous studies, we delimited two distinct phylogeographical lineages within the nominal species and found an intermediate level of differentiation in their multidimensional environmental niches, suggesting the possibility for local adaptation. We then compared projections of the species' habitat suitability under climate change scenarios using species-level and lineage-level SDMs. In the Central Tropical Indo-Pacific region, models for both levels predicted considerable range contraction in the future, but the lineage-level models predicted more severe habitat loss. Importantly, the two modelling approaches predicted opposite patterns of habitat change in the Western Tropical Indo-Pacific region. Our results highlight the necessity of conserving distinct populations and genetic pools to avoid regional extinction due to climate change and have important implications for guiding future management of seagrasses.

Seagrasses in an era of ocean warming: a review

Seagrasses are valuable sources of food and habitat for marine life and are one of Earth's most efficient carbon sinks. However, they are facing a global decline due to ocean warming and eutrophication. In the last decade, with the advent of new technology and molecular advances, there has been a dramatic increase in the number of studies focusing on the effects of ocean warming on seagrasses. Here, we provide a comprehensive review of the future of seagrasses in an era of ocean warming. We have gathered information from published studies to identify potential commonalities in the effects of warming and the responses of seagrasses across four distinct levels: molecular, biochemical/physiological, morphological/population, and ecosystem/planetary. To date, we know that although warming strongly affects seagrasses at all four levels, seagrass responses diverge amongst species, populations, and over depths. Furthermore, warming alters seagrass distribution causing massive die-offs in some seagrass populations, whilst also causing tropicalization and migration of temperate species. In this review, we evaluate the combined effects of ocean warming with other environmental stressors and emphasize the need for multiple-stressor studies to provide a deeper understanding of seagrass resilience. We conclude by discussing the most significant knowledge gaps and future directions for seagrass research.

Fast and pervasive transcriptomic resilience and acclimation of extremely heat-tolerant coral holobionts from the northern Red Sea

Corals from the northern Red Sea and Gulf of Aqaba exhibit extreme thermal tolerance. To examine the underlying gene expression dynamics, we exposed Stylophora pistillata from the Gulf of Aqaba to short-term (hours) and long-term (weeks) heat stress with peak seawater temperatures ranging from their maximum monthly mean of 27 °C (baseline) to 29.5 °C, 32 °C, and 34.5 °C. Corals were sampled at the end of the heat stress as well as after a recovery period at baseline temperature. Changes in coral host and symbiotic algal gene expression were determined via RNA-sequencing (RNA-Seq). Shifts in coral microbiome composition were detected by complementary DNA (cDNA)-based 16S ribosomal RNA (rRNA) gene sequencing. In all experiments up to 32 °C, RNA-Seq revealed fast and pervasive changes in gene expression, primarily in the coral host, followed by a return to baseline gene expression for the majority of coral (>94%) and algal (>71%) genes during recovery. At 34.5 °C, large differences in gene expression were observed with minimal recovery, high coral mortality, and a microbiome dominated by opportunistic bacteria (including Vibrio species), indicating that a lethal temperature threshold had been crossed. Our results show that the S. pistillata holobiont can mount a rapid and pervasive gene expression response contingent on the amplitude and duration of the thermal stress. We propose that the transcriptomic resilience and transcriptomic acclimation observed are key to the extraordinary thermal tolerance of this holobiont and, by inference, of other northern Red Sea coral holobionts, up to seawater temperatures of at least 32 °C, that is, 5 °C above their current maximum monthly mean.

Phenotypic plasticity under rapid global changes: The intrinsic force for future seagrasses survival

Coastal oceans are particularly affected by rapid and extreme environmental changes with dramatic consequences for the entire ecosystem. Seagrasses are key ecosystem engineering or foundation species supporting diverse and productive ecosystems along the coastline that are particularly susceptible to fast environmental changes. In this context, the analysis of phenotypic plasticity could reveal important insights into seagrasses persistence, as it represents an individual property that allows species' phenotypes to accommodate and react to fast environmental changes and stress. Many studies have provided different definitions of plasticity and related processes (acclimation and adaptation) resulting in a variety of associated terminology. Here, we review different ways to define phenotypic plasticity with particular reference to seagrass responses to single and multiple stressors. We relate plasticity to the shape of reaction norms, resulting from genotype by environment interactions, and examine its role in the presence of environmental shifts. The potential role of genetic and epigenetic changes in underlying seagrasses plasticity in face of environmental changes is also discussed. Different approaches aimed to assess local acclimation and adaptation in seagrasses are explored, explaining strengths and weaknesses based on the main results obtained from the most recent literature. We conclude that the implemented experimental approaches, whether performed with controlled or field experiments, provide new insights to explore the basis of plasticity in seagrasses. However, an improvement of molecular analysis and the application of multi-factorial experiments are required to better explore genetic and epigenetic adjustments to rapid environmental shifts. These considerations revealed the potential for selecting the best phenotypes to promote assisted evolution with fundamental implications on restoration and preservation efforts.

The Genetic Component of Seagrass Restoration: What We Know and the Way Forwards

Seagrasses are marine flowering plants providing key ecological services and functions in coasts and estuaries across the globe. Increased environmental changes fueled by human activities are affecting their existence, compromising natural habitats and ecosystems’ biodiversity and functioning. In this context, restoration of disturbed seagrass environments has become a worldwide priority to reverse ecosystem degradation and to recover ecosystem functionality and associated services. Despite the proven importance of genetic research to perform successful restoration projects, this aspect has often been overlooked in seagrass restoration. Here, we aimed to provide a comprehensive perspective of genetic aspects related to seagrass restoration. To this end, we first reviewed the importance of studying the genetic diversity and population structure of target seagrass populations; then, we discussed the pros and cons of different approaches used to restore and/or reinforce degraded populations. In general, the collection of genetic information and the development of connectivity maps are critical steps for any seagrass restoration activity. Traditionally, the selection of donor population preferred the use of local gene pools, thought to be the best adapted to current conditions. However, in the face of rapid ocean changes, alternative approaches such as the use of climate-adjusted or admixture genotypes might provide more sustainable options to secure the survival of restored meadows. Also, we discussed different transplantation strategies applied in seagrasses and emphasized the importance of long-term seagrass monitoring in restoration. The newly developed information on epigenetics as well as the application of assisted evolution strategies were also explored. Finally, a view of legal and ethical issues related to national and international restoration management is included, highlighting improvements and potential new directions to integrate with the genetic assessment. We concluded that a good restoration effort should incorporate: (1) a good understanding of the genetic structure of both donors and populations being restored; (2) the analysis of local environmental conditions and disturbances that affect the site to be restored; (3) the analysis of local adaptation constraints influencing the performances of donor populations and native plants; (4) the integration of distribution/connectivity maps with genetic information and environmental factors relative to the target seagrass populations; (5) the planning of long-term monitoring programs to assess the performance of the restored populations. The inclusion of epigenetic knowledge and the development of assisted evolution programs are strongly hoped for the future.

Plasticity in parental effects confers rapid larval thermal tolerance in the estuarine anemone Nematostella vectensis

Parental effects can prepare offspring for different environments and facilitate survival across generations. We exposed parental populations of the estuarine anemone, Nematostella vectensis, from Massachusetts to elevated temperatures and quantified larval mortality across a temperature gradient. We found that parental exposure to elevated temperatures resulted in a consistent increase in larval thermal tolerance, as measured by the temperature at which 50% of larvae die (LT50), with a mean increase in LT50 of 0.3°C. Larvae from subsequent spawns returned to baseline thermal thresholds when parents were returned to normal temperatures, indicating plasticity in these parental effects. Histological analyses of gametogenesis in females suggested that these dynamic shifts in larval thermal tolerance may be facilitated by maternal effects in non-overlapping gametic cohorts. We also compared larvae from North Carolina (a genetically distinct population with higher baseline thermal tolerance) and Massachusetts parents, and observed that larvae from heat-exposed Massachusetts parents had thermal thresholds comparable to those of larvae from unexposed North Carolina parents. North Carolina parents also increased larval thermal tolerance under the same high-temperature regime, suggesting that plasticity in parental effects is an inherent trait for N. vectensis Overall, we find that larval thermal tolerance in N. vectensis shows a strong genetic basis and can be modulated by parental effects. Further understanding of the mechanisms behind these shifts can elucidate the fate of thermally sensitive ectotherms in a rapidly changing thermal environment.

A framework for understanding gene expression plasticity and its influence on stress tolerance

Phenotypic plasticity can serve as a stepping stone towards adaptation. Recently, studies have shown that gene expression contributes to emergent stress responses such as thermal tolerance, with tolerant and susceptible populations showing distinct transcriptional profiles. However, given the dynamic nature of gene expression, interpreting transcriptomic results in a way that elucidates the functional connection between gene expression and the observed stress response is challenging. Here, we present a conceptual framework to guide interpretation of gene expression reaction norms in the context of stress tolerance. We consider the evolutionary and adaptive potential of gene expression reaction norms and discuss the influence of sampling timing, transcriptomic resilience, as well as complexities related to life history when interpreting gene expression dynamics and how these patterns relate to host tolerance. We highlight corals as a case study to demonstrate the value of this framework for non-model systems. As species face rapidly changing environmental conditions, modulating gene expression can serve as a mechanistic link from genetic and cellular processes to the physiological responses that allow organisms to thrive under novel conditions. Interpreting how or whether a species can employ gene expression plasticity to ensure short-term survival will be critical for understanding the global impacts of climate change across diverse taxa.

The Seagrass Holobiont: What We Know and What We Still Need to Disclose for Its Possible Use as an Ecological Indicator

Microbes and seagrass establish symbiotic relationships constituting a functional unit called the holobiont that reacts as a whole to environmental changes. Recent studies have shown that the seagrass microbial associated community varies according to host species, environmental conditions and the host’s health status, suggesting that the microbial communities respond rapidly to environmental disturbances and changes. These changes, dynamics of which are still far from being clear, could represent a sensitive monitoring tool and ecological indicator to detect early stages of seagrass stress. In this review, the state of art on seagrass holobiont is discussed in this perspective, with the aim of disentangling the influence of different factors in shaping it. As an example, we expand on the widely studied Halophila stipulacea’s associated microbial community, highlighting the changing and the constant components of the associated microbes, in different environmental conditions. These studies represent a pivotal contribution to understanding the holobiont’s dynamics and variability pattern, and to the potential development of ecological/ecotoxicological indices. The influences of the host’s physiological and environmental status in changing the seagrass holobiont, alongside the bioinformatic tools for data analysis, are key topics that need to be deepened, in order to use the seagrass-microbial interactions as a source of ecological information.

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

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

The Seagrass Methylome Is Associated With Variation in Photosynthetic Performance Among Clonal Shoots

Evolutionary theory predicts that clonal organisms are more susceptible to extinction than sexually reproducing organisms, due to low genetic variation and slow rates of evolution. In agreement, conservation management considers genetic variation as the ultimate measure of a population's ability to survive over time. However, clonal plants are among the oldest living organisms on our planet. Here, we test the hypothesis that clonal seagrass meadows display epigenetic variation that complements genetic variation as a source of phenotypic variation. In a clonal meadow of the seagrass Zostera marina, we characterized DNA methylation among 42 shoots. We also sequenced the whole genome of 10 shoots to correlate methylation patterns with photosynthetic performance under exposure to and recovery from 27°C, while controlling for somatic mutations. Here, we show for the first time that clonal seagrass shoots display DNA methylation variation that is independent from underlying genetic variation, and associated with variation in photosynthetic performance under experimental conditions. It remains unknown to what degree this association could be influenced by epigenetic responses to transplantation-related stress, given that the methylomes showed a strong shift under acclimation to laboratory conditions. The lack of untreated control samples in the heat stress experiment did not allow us to distinguish methylome shifts induced by acclimation from such induced by heat stress. Notwithstanding, the co-variation in DNA methylation and photosynthetic performance may be linked via gene expression because methylation patterns varied in functionally relevant genes involved in photosynthesis, and in the repair and prevention of heat-induced protein damage. While genotypic diversity has been shown to enhance stress resilience in seagrass meadows, we suggest that epigenetic variation plays a similar role in meadows dominated by a single genotype. Consequently, conservation management of clonal plants should consider epigenetic variation as indicator of resilience and stability.

Root microbiomes as indicators of seagrass health

The development of early warning indicators that identify ecosystem stress is a priority for improving ecosystem management. As microbial communities respond rapidly to environmental disturbance, monitoring their composition could prove one such early indicator of environmental stress. We combined 16S rRNA gene sequencing of the seagrass root microbiome of Halophila ovalis with seagrass health metrics (biomass, productivity and Fsulphide) to develop microbial indicators for seagrass condition across the Swan-Canning Estuary and the Leschenault Estuary (south-west Western Australia); the former had experienced an unseasonal rainfall event leading to declines in seagrass health. Microbial indicators detected sites of potential stress that other seagrass health metrics failed to detect. Genera that were more abundant in 'healthy' seagrasses included putative methylotrophic bacteria (e.g. Methylotenera and Methylophaga), iron cycling bacteria (e.g. Deferrisoma and Geothermobacter) and N2 fixing bacteria (e.g. Rhizobium). Conversely, genera that were more abundant in 'stressed' seagrasses were dominated by putative sulphur-cycling bacteria, both sulphide-oxidising (e.g. Candidatus Thiodiazotropha and Candidatus Electrothrix) and sulphate-reducing (e.g. SEEP-SRB1, Desulfomonile and Desulfonema). The sensitivity of the microbial indicators developed here highlights their potential to be further developed for use in adaptive seagrass management, and emphasises their capacity to be effective early warning indicators of stress.

Elevated Temperatures Impose Transcriptional Constraints and Elicit Intraspecific Differences Between Coffee Genotypes

The projected impact of global warming on coffee production may require the heat-adapted genotypes in the next decades. To identify cellular strategies in response to warmer temperatures, we compared the effect of elevated temperature on two commercial Coffea arabica L. genotypes exploring leaf physiology, transcriptome, and carbohydrate/protein composition. Growth temperatures were 23/19°C (day/night), as optimal condition (OpT), and 30/26°C (day/night) as a possible warmer scenario (WaT). The cv. Acauã showed lower levels of leaf temperature (Tleaf) under both conditions compared to cv. Catuaí, whereas slightly or no differences for other leaf physiological parameters. Therefore, to explore temperature responsive pathways the leaf transcriptome was examined using RNAseq. Genotypes showed a marked number of differentially-expressed genes (DEGs) under OpT, however DEGs strongly decrease in both at WaT condition indicating a transcriptional constraint. DEGs responsive to WaT revealed shared and genotype-specific genes mostly related to carbohydrate metabolism. Under OpT, leaf starch content was greater in cv. Acauã and, as WaT temperature was imposed, the leaf soluble sugar did not change in contrast to cv. Catuaí, although the levels of leaf starch, sucrose, and leaf protein decreased in both genotypes. These findings revealed intraspecific differences in the underlying transcriptional and metabolic interconnected pathways responsive to warmer temperatures, which is potentially linked to thermotolerance, and thus may be useful as biomarkers in breeding for a changing climate.

Multigenerational Mitigating Effects of Ocean Acidification on In Vivo Endpoints, Antioxidant Defense, DNA Damage Response, and Epigenetic Modification in an Asexual Monogonont Rotifer

Ocean acidification (OA) is caused by changes in ocean carbon chemistry due to increased atmospheric pCO₂ and is predicted to have deleterious effects on marine ecosystems. While the potential impacts of OA on many marine species have been studied, the multigenerational effects on asexual organisms remain unknown. We found that low seawater pH induced oxidative stress and DNA damage, decreasing growth rates, fecundity, and lifespans in the parental generation, whereas deleterious effects on in vivo endpoints in F1 and F2 offspring were less evident. The findings suggest that multigenerational adaptive effects play a role in antioxidant abilities and other defense mechanisms. OA-induced DNA damage, including double-strand breaks (DSBs), was fully repaired in F1 offspring of parents exposed to OA for 7 days, indicating that an adaptation mechanism may be the major driving force behind multigenerational adaptive effects. Analysis of epigenetic modification in response to OA involved examination of histone modification of DNA repair genes and a chromatin immunoprecipitation assay, as Bombus koreanus has no methylation pattern for CpG in its genome. We conclude that DSBs, DNA repair, and histone modification play important roles in multigenerational plasticity in response to OA in an asexual monogonont rotifer.

Nitrogen sufficiency enhances thermal tolerance in habitat-forming kelp: implications for acclimation under thermal stress

Local and global changes associated with anthropogenic activities are impacting marine and terrestrial ecosystems. Macroalgae, especially habitat-forming species like kelp, play critical roles in temperate coastal ecosystems. However, their abundance and distribution patterns have been negatively affected by warming in many regions around the globe. Along with global change, coastal ecosystems are also impacted by local drivers such as eutrophication. The interaction between global and local drivers might modulate kelp responses to environmental change. This study examines the regulatory effect of NO₃⁻ on the thermal plasticity of the giant kelp Macrocystis pyrifera. To do this, thermal performance curves (TPCs) of key temperature-dependant traits–growth, photosynthesis, NO₃⁻ assimilation and chlorophyll a fluorescence–were examined under nitrate replete and deplete conditions in a short-term incubation. We found that thermal plasticity was modulated by NO₃⁻ but different thermal responses were observed among traits. Our study reveals that nitrogen, a local driver, modulates kelp responses to high seawater temperatures, ameliorating the negative impacts on physiological performance (i.e. growth and photosynthesis). However, this effect might be species-specific and vary among biogeographic regions – thus, further work is needed to determine the generality of our findings to other key temperate macroalgae that are experiencing temperatures close to their thermal tolerance due to climate change.

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

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

Stress Memory of Recurrent Environmental Challenges in Marine Invasive Species: Ciona robusta as a Case Study

Fluctuating environmental changes impose tremendous stresses on sessile organisms in marine ecosystems, in turn, organisms develop complex response mechanisms to keep adaptive homeostasis for survival. Physiological plasticity is one of the primary lines of defense against environmental challenges, and such defense often relies on the antioxidant defense system (ADS). Hence, it is imperative to understand response mechanisms of ADS to fluctuating environments. Invasive species provide excellent models to study how species cope with environmental stresses, as invasive species encounter sudden, and often recurrent, extensive environmental challenges during the whole invasion process. Here, we studied the roles of ADS on rapid response to recurrent cold challenges in a highly invasive tunicate (Ciona robusta) by simulating cold stresses during its invasion process. We assessed antioxidative indicators, including malondialdehyde (MDA), total antioxidant capacity (T-AOC), superoxide dismutase (SOD), catalase (CAT), and glutathione (GSH), as well as transcriptional changes of ADS-related genes to reveal the physiological plasticity under recurring cold stresses. Our results demonstrated that physiological homeostasis relied on the resilience of ADS, which further accordingly tuned antioxidant activity and gene expression to changing environments. The initial cold stress remodeled baselines of ADS to promote the development of stress memory, and subsequent stress memory largely decreased the physiological response to recurrent environmental challenges. All results here suggest that C. robusta could develop stress memory to maintain physiological homeostasis in changing or harsh environments. The results obtained in this study provide new insights into the mechanism of rapid physiological adaption during biological invasions.

Response of foundation macrophytes to near-natural simulated marine heatwaves

Marine heatwaves have been observed worldwide and are expected to increase in both frequency and intensity due to climate change. Such events may cause ecosystem reconfigurations arising from species range contraction or redistribution, with ecological, economic and social implications. Macrophytes such as the brown seaweed Fucus vesiculosus and the seagrass Zostera marina are foundation species in many coastal ecosystems of the temperate northern hemisphere. Hence, their response to extreme events can potentially determine the fate of associated ecosystems. Macrophyte functioning is intimately linked to the maintenance of photosynthesis, growth and reproduction, and resistance against pathogens, epibionts and grazers. We investigated morphological, physiological, pathological and chemical defence responses of western Baltic Sea F. vesiculosus and Z. marina populations to simulated near-natural marine heatwaves. Along with (a) the control, which constituted no heatwave but natural stochastic temperature variability (0HW), two treatments were applied: (b) two late-spring heatwaves (June, July) followed by a summer heatwave (August; 3HW) and (c) a summer heatwave only (1HW). The 3HW treatment was applied to test whether preconditioning events can modulate the potential sensitivity to the summer heatwave. Despite the variety of responses measured in both species, only Z. marina growth was impaired by the accumulative heat stress imposed by the 3HW treatment. Photosynthetic rate, however, remained high after the last heatwave indicating potential for recovery. Only epibacterial abundance was significantly affected in F. vesiculosus. Hence both macrophytes, and in particular F. vesiculosus, seem to be fairly tolerant to short-term marine heatwaves at least at the intensities applied in this experiment (up to 5°C above mean temperature over a period of 9 days). This may partly be due to the fact that F. vesiculosus grows in a highly variable environment, and may have a high phenotypic plasticity.

Transcriptome profiling analysis of the seagrass, Zostera muelleri under copper stress

Copper (Cu) in an essential trace metal but it can also contaminate coastal waters at high concentrations mainly from agricultural run-off and mining activities which are detrimental to marine organisms including seagrasses. The molecular mechanisms driving Cu toxicity in seagrasses are not clearly understood yet. Here, we investigated the molecular responses of the Australian seagrass, Z. muelleri at the whole transcriptomic level after 7 days of exposure to 250 μg Cu L^-1 and 500 μg Cu L^-1. The leaf-specific whole transcriptome results showed a concentration-dependent disturbance in chloroplast function, regulatory stress responses and defense mechanisms. This study provided new insights into the responses of seagrasses to trace metal stress and reports possible candidate genes which can be considered as biomarkers to improve conservation and management of seagrass meadows.

Differential Leaf Age-Dependent Thermal Plasticity in the Keystone Seagrass Posidonia oceanica

Introduction: Gene-expression patterns and their upstream regulatory mechanisms (e.g. epigenetic) are known to modulate plant acclimatability and thus tolerance to heat stress. Within species, thermal plasticity (i.e. temperature-sensitive phenotypic plasticity) and differential thermo-tolerance are recognized among different genotypes, development stages, organs or tissues. Leaf age and lifespan have been demonstrated to strongly affect photosynthetic thermo-tolerance in terrestrial species, whereas there is no information available for marine plants. Materials and Methods: Here, we investigated how an intense warming event affects molecular and photo-physiological functions in the large-sized seagrass Posidonia oceanica, at fine spatial resolution. Plants were exposed for one week at 34°C in a controlled-mesocosm system. Subsequent variations in the expression of 12 target genes and global DNA methylation level were evaluated in three leaf-age sections (i.e. basal, medium and high) established along the longitudinal axis of youngest, young and fully mature leaves of the shoot. Targeted genes were involved in photosynthesis, chlorophyll biosynthesis, energy dissipation mechanisms, stress response and programmed cell death. Molecular analyses paralleled the assessment of pigment content and photosynthetic performance of the same leaf segments, as well as of plant growth inhibition under acute warming. Results: Our data revealed, for the first time, the presence of variable leaf age-dependent stress-induced epigenetic and gene-expression changes in seagrasses, underlying photo-physiological and growth responses to heat stress. An investment in protective responses and growth arrest was observed in immature tissues; while mature leaf sections displayed a higher ability to offset gene down-regulation, possibly through the involvement of DNA methylation changes, although heat-induced damages were visible at photo-physiological level. Discussion: Overall, mature and young leaf tissues exhibited different strategies to withstand heat stress and thus a variable thermal plasticity. This should be taken in consideration when addressing seagrass response to warming and other stressors, especially in large-sized species, where sharp age differences are present within and among leaves, and other gradients of environmental factors (e.g. light) could be at play. Molecular and physiological evaluations conducted only on adult leaf tissues, as common practice in seagrass research, could give inadequate estimates of the overall plant state, and should not be considered as a proxy for the whole shoot.

SMRT sequencing of full-length transcriptome of seagrasses Zostera japonica

Seagrass meadows are among the four most productive marine ecosystems in the world. Zostera japonica (Z. japonica) is the most widely distributed species of seagrass in China. However, there is no reference genome or transcriptome available for Z. japonica, impeding progress in functional genomic and molecular ecology studies in this species. Temperature is the main factor that controls the distribution and growth of seagrass around the world, yet how seagrass responds to heat stress remains poorly understood due to the lack of genomic and transcriptomic data. In this study, we applied a combination of second- and third-generation sequencing technologies to sequence full-length transcriptomes of Z. japonica. In total, we obtained 58,134 uniform transcripts, which included 46,070 high-quality full-length transcript sequences. We identified 15,411 simple sequence repeats, 258 long non-coding RNAs and 28,038 open reading frames. Exposure to heat elicited a complex transcriptional response in genes involved in posttranslational modification, protein turnover and chaperones. Overall, our study provides the first large-scale full-length trascriptome in Zostera japonica, allowing for structural, functional and comparative genomics studies in this important seagrass species. Although previous studies have focused specifically on heat shock proteins, we found that examination of other heat stress related genes is important for studying response to heat stress in seagrass. This study provides a genetic resource for the discovery of genes related to heat stress tolerance in this species. Our transcriptome can be further utilized in future studies to understand the molecular adaptation to heat stress in Zostera japonica.

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

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

Effects of adult temperature on gene expression in a butterfly: identifying pathways associated with thermal acclimation

BACKGROUND

Phenotypic plasticity is a pervasive property of all organisms and considered to be of key importance for dealing with environmental variation. Plastic responses to temperature, which is one of the most important ecological factors, have received much attention over recent decades. A recurrent pattern of temperature-induced adaptive plasticity includes increased heat tolerance after exposure to warmer temperatures and increased cold tolerance after exposure to cooler temperatures. However, the mechanisms underlying these plastic responses are hitherto not well understood. Therefore, we here investigate effects of adult acclimation on gene expression in the tropical butterfly Bicyclus anynana, using an RNAseq approach.

RESULTS

We show that several antioxidant markers (e.g. peroxidase, cytochrome P450) were up-regulated at a higher temperature compared with a lower adult temperature, which might play an important role in the acclamatory responses subsequently providing increased heat tolerance. Furthermore, several metabolic pathways were up-regulated at the higher temperature, likely reflecting increased metabolic rates. In contrast, we found no evidence for a decisive role of the heat shock response.

CONCLUSIONS

Although the important role of antioxidant defence mechanisms in alleviating detrimental effects of oxidative stress is firmly established, we speculate that its potentially important role in mediating heat tolerance and survival under stress has been underestimated thus far and thus deserves more attention.

The transcriptome of the marine calanoid copepod Temora longicornis under heat stress and recovery

Understanding the impacts of global change in zooplankton communities is crucial, as alterations in the zooplankton communities can affect entire marine ecosystems. Despite the economic and ecological importance of the calanoid copepod Temora longicornis in the Belgian part of the North Sea, molecular data is still very limited for this species. Using HiSeq Illumina sequencing, we sequenced the whole transcriptome of T. longicornis, after being exposed to realistic temperatures of 14 and 17 °C. After both an acute (1 day) and a more sustained (5 days) thermal exposure to 17 °C, we investigated gene expression differences with animals exposed to 14 °C, which may be critical for the thermal acclimation and resilience of this copepod species. We also studied the possibility of a short term stress recovery of a heat shock. A total of 179,569 transcripts were yielded, of which 44,985 putative ORF transcripts were identified. These transcripts were subsequently annotated into roughly 22,000 genes based on known sequences using Gene Ontology (GO) and KEGG databases. Temora only showed a mild response to both the temperature and the duration of the exposure. We found that the expression of 27 transcripts varied significantly with an increase in temperature of 3 °C, of which eight transcripts were differentially expressed after acute exposure only. Gene set enrichment analysis revealed that, overall, T. longicornis was more impacted by a sustained thermal exposure, rather than an immediate (acute) exposure, with two times as many enriched GO terms in the sustained treatment. We also identified several general stress responses independent of exposure time, such as modified protein synthesis, energy mobilisation, cuticle and chaperone proteins. Finally, we highlighted candidate genes of a possible recovery from heat exposure, identifying similar terms as those enriched in the heat treatments, i.e. related to for example energy metabolism, cuticle genes and extracellular matrix. The data presented in this study provides the first transcriptome available for T. longicornis which can be used for future genomic studies.

Which Genes in a Typical Intertidal Seagrass (Zostera japonica) Indicate Copper-, Lead-, and Cadmium Pollution?

Healthy seagrasses are considered a prime indicator of estuarine and coastal ecosystem function; however, as the only group of flowering plants recolonizing the sea, seagrasses are frequently exposed to anthropogenic heavy metal pollutants, which are associated with high levels of molecular damage. To determine whether biologically relevant concentrations of heavy metals cause systematic alterations in RNA expression patterns, we performed a gene expression study using transcriptome analyses (RNA-seq). We exposed the typical intertidal seagrass Zostera japonica to 0 and 50 μM of copper (Cu), lead (Pb), and cadmium (Cd) under laboratory conditions. A total of 18,266 differentially expressed genes (DEGs) were identified, of which 2001 co-expressed genes directly related by Cu, Pb, and Cd stress. We also examined the effects of short-term heavy metal Cu, Pb, and Cd pulses on the accumulation of metals in Z. japonica and showed metal concentrations were higher in the shoots than in roots. Twelve differentially expressed genes were further analyzed for expression differences using real-time quantitative polymerase chain reaction (RT-qPCR). Our data suggest that as coastal seawater pollution worsens, the sensitive genes identified in this study may be useful biomarkers of sublethal effects and provide fundamental information for Z. japonica resistant gene engineering.

Investigating cellular stress response to heat stress in the seagrass Posidonia oceanica in a global change scenario

Posidonia oceanica meadows are facing global threats mainly due to episodic heat waves. In a mesocosm experiment, we aimed at disentangling the molecular response of P. oceanica under increasing temperature (20 °C-32 °C). The experiment was carried out in spring, when heat waves can potentially occur and plants are putatively more sensitive to heat stress, since they are deprived in carbohydrates reserves after the cold winter months. We aimed to identify the activation of different phases of the cellular stress response (CSR) reaction and the responsive genes activated or repressed in heated plants. A molecular traffic light was proposed as a response model including green (protein folding and membrane protection), yellow (ubiquitination and proteolysis) and red (DNA repair and apoptosis) categories. Additionally, we estimated phenological trait variations to complement the information obtained from the molecular proxies of stress. Despite reduced leaf growth rate, heated plants did not exhibit signs of irreversible damage, probably underlying species pre-adaptation to warm and fluctuating regimes. Gene expression analyses revealed that molecular chaperoning, DNA repair and apoptosis inhibition processes related genes were the ones that mostly responded to high thermal stress and will be target of further investigation and in situ proofing for assessing their use as indicators of P. oceanica performance under sub-lethal heat stress.

Differentiation in fitness-related traits in response to elevated temperatures between leading and trailing edge populations of marine macrophytes

The nature of species distribution boundaries is a key subject in ecology and evolution. Edge populations are potentially more exposed to climate-related environmental pressures. Despite research efforts, little is known about variability in fitness-related traits in leading (i.e., colder, high latitude) versus trailing (i.e., warmer, low latitude) edge populations. We tested whether the resilience, i.e. the resistance and recovery, of key traits differs between a distributional cold (Greenland) and warm (Portugal) range edge population of two foundation marine macrophytes, the intertidal macroalga Fucus vesiculosus and the subtidal seagrass Zostera marina. The resistance and recovery of edge populations to elevated seawater temperatures was compared under common experimental conditions using photosynthetic efficiency and expression of heat shock proteins (HSP). Cold and warm edge populations differed in their response, but this was species specific. The warm edge population of F. vesiculosus showed higher thermal resistance and recovery whereas the cold leading edge was less tolerant. The opposite was observed in Z. marina, with reduced recovery at the warm edge, while the cold edge was not markedly affected by warming. Our results confirm that differentiation of thermal stress responses can occur between leading and trailing edges, but such responses depend on local population traits and are thus not predictable just based on thermal pressures.