Thermo-priming triggers species-specific physiological and transcriptome responses in Mediterranean seagrasses.

Heat-priming improves plants' tolerance to a recurring heat stress event. The underlying molecular mechanisms of heat-priming are largely unknown in seagrasses. Here, ad hoc mesocosm experiments were conducted with two Mediterranean seagrass species, Posidonia oceanica and Cymodocea nodosa. Plants were first exposed to heat-priming, followed by a heat-triggering event. A comprehensive assessment of plant stress response across different levels of biological organization was performed at the end of the triggering event. Morphological and physiological results showed an improved response of heat-primed P. oceanica plants while in C. nodosa both heat- and non-primed plants enhanced their growth rates at the end of the triggering event. As resulting from whole transcriptome sequencing, molecular functions related to several cellular compartments and processes were involved in the response to warming of non-primed plants, while the response of heat-primed plants involved a limited group of processes. Our results suggest that seagrasses acquire a primed state during the priming event, that eventually gives plants the ability to induce a more energy-effective response when the thermal stress event recurs. Different species may differ in their ability to perform an improved heat stress response after priming. This study provides pioneer molecular insights into the emerging topic of seagrass stress priming and may benefit future studies in the field.

Seagrass genomes reveal ancient polyploidy and adaptations to the marine environment.

We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinum and Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marina and P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the 'savannahs of the sea' are of major concern in times of climate change and loss of biodiversity.

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.

A trait-based framework for seagrass ecology: Trends and prospects.

In the last three decades, quantitative approaches that rely on organism traits instead of taxonomy have advanced different fields of ecological research through establishing the mechanistic links between environmental drivers, functional traits, and ecosystem functions. A research subfield where trait-based approaches have been frequently used but poorly synthesized is the ecology of seagrasses; marine angiosperms that colonized the ocean 100M YA and today make up productive yet threatened coastal ecosystems globally. Here, we compiled a comprehensive trait-based response-effect framework (TBF) which builds on previous concepts and ideas, including the use of traits for the study of community assembly processes, from dispersal and response to abiotic and biotic factors, to ecosystem function and service provision. We then apply this framework to the global seagrass literature, using a systematic review to identify the strengths, gaps, and opportunities of the field. Seagrass trait research has mostly focused on the effect of environmental drivers on traits, i.e., "environmental filtering" (72%), whereas links between traits and functions are less common (26.9%). Despite the richness of trait-based data available, concepts related to TBFs are rare in the seagrass literature (15% of studies), including the relative importance of neutral and niche assembly processes, or the influence of trait dominance or complementarity in ecosystem function provision. These knowledge gaps indicate ample potential for further research, highlighting the need to understand the links between the unique traits of seagrasses and the ecosystem services they provide.

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.

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.

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.

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.

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.

m6A RNA Methylation in Marine Plants: First Insights and Relevance for Biological Rhythms.

Circadian regulations are essential for enabling organisms to synchronize physiology with environmental light-dark cycles. Post-transcriptional RNA modifications still represent an understudied level of gene expression regulation in plants, although they could play crucial roles in environmental adaptation. N6-methyl-adenosine (m6A) is the most prevalent mRNA modification, established by "writer" and "eraser" proteins. It influences the clockwork in several taxa, but only few studies have been conducted in plants and none in marine plants. Here, we provided a first inventory of m6A-related genes in seagrasses and investigated daily changes in the global RNA methylation and transcript levels of writers and erasers in Cymodocea nodosa and Zostera marina. Both species showed methylation peaks during the dark period under the same photoperiod, despite exhibiting asynchronous changes in the m6A profile and related gene expression during a 24-h cycle. At contrasting latitudes, Z. marina populations displayed overlapping daily patterns of the m6A level and related gene expression. The observed rhythms are characteristic for each species and similar in populations of the same species with different photoperiods, suggesting the existence of an endogenous circadian control. Globally, our results indicate that m6A RNA methylation could widely contribute to circadian regulation in seagrasses, potentially affecting the photo-biological behaviour of these plants.

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.

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.

Molecular level responses to chronic versus pulse nutrient loading in the seagrass Posidonia oceanica undergoing herbivore pressure.

Seagrasses are key marine foundation species, currently declining due to the compounded action of global and regional anthropogenic stressors. Eutrophication has been associated with seagrass decline, while grazing has been traditionally considered to be a natural disturbance with a relatively low impact on seagrasses. In the recent years, this assumption has been revisited. Here, by means of a 16-month field-experiment, we investigated the molecular mechanisms driving the long-term response of Posidonia oceanica to the combination of nutrient enrichment, either as a chronic (press) or pulse disturbance, and herbivory. Changes in expression levels of 19 target genes involved in key steps of photosynthesis, nutrient assimilation, chlorophyll metabolism, oxidative-stress response and plant defense were evaluated through reverse transcription-quantitative polymerase chain reaction (RT-qPCR). High herbivore pressure affected the molecular response of P. oceanica more dramatically than did enhanced nutrient levels, altering the expression of genes involved in plant tolerance and resistance traits, such as photosynthesis and defense mechanisms. Genes involved in carbon fixation and N assimilation modulated the response of plants to high nutrient levels. Availability of resources seems to modify P. oceanica response to herbivory, where the upregulation of a nitrate transporter gene was accompanied by the decline in the expression of nitrate reductase in the leaves, suggesting a change in plant-nutrient allocation strategy. Finally, press and pulse fertilizations altered nitrate uptake and reduction-related genes in opposite ways, suggesting that taking into account the temporal regime of nutrient loading is important to assess the physiological response of seagrasses to eutrophication.

Genomewide transcriptional reprogramming in the seagrass Cymodocea nodosa under experimental ocean acidification.

Here, we report the first use of massive-scale RNA-sequencing to explore seagrass response to CO2 -driven ocean acidification (OA). Large-scale gene expression changes in the seagrass Cymodocea nodosa occurred at CO2 levels projected by the end of the century. C. nodosa transcriptome was obtained using Illumina RNA-Seq technology and de novo assembly, and differential gene expression was explored in plants exposed to short-term high CO2 /low pH conditions. At high pCO2 , there was a significant increased expression of transcripts associated with photosynthesis, including light reaction functions and CO2 fixation, and also to respiratory pathways, specifically for enzymes involved in glycolysis, in the tricarboxylic acid cycle and in the energy metabolism of the mitochondrial electron transport. The upregulation of respiratory metabolism is probably supported by the increased availability of photosynthates and increased energy demand for biosynthesis and stress-related processes under elevated CO2 and low pH. The upregulation of several chaperones resembling heat stress-induced changes in gene expression highlighted the positive role these proteins play in tolerance to intracellular acid stress in seagrasses. OA further modifies C. nodosa secondary metabolism inducing the transcription of enzymes related to biosynthesis of carbon-based secondary compounds, in particular the synthesis of polyphenols and isoprenoid compounds that have a variety of biological functions including plant defence. By demonstrating which physiological processes are most sensitive to OA, this research provides a major advance in the understanding of seagrass metabolism in the context of altered seawater chemistry from global climate change.

Linking gene expression to productivity to unravel long- and short-term responses of seagrasses exposed to CO2 in volcanic vents.

Ocean acidification is a major threat for marine life but seagrasses are expected to benefit from high CO2. In situ (long-term) and transplanted (short-term) plant incubations of the seagrass Cymodocea nodosa were performed near and away the influence of volcanic CO2 vents at Vulcano Island to test the hypothesis of beneficial effects of CO2 on plant productivity. We relate, for the first time, the expression of photosynthetic, antioxidant and metal detoxification-related genes to net plant productivity (NPP). Results revealed a consistent pattern between gene expression and productivity indicating water origin as the main source of variability. However, the hypothesised beneficial effect of high CO2 around vents was not supported. We observed a consistent long- and short-term pattern of gene down-regulation and 2.5-fold NPP decrease in plants incubated in water from the vents and a generalized up-regulation and NPP increase in plants from the vent site incubated with water from the Reference site. Contrastingly, NPP of specimens experimentally exposed to a CO2 range significantly correlated with CO2 availability. The down-regulation of metal-related genes in C. nodosa leaves exposed to water from the venting site suggests that other factors than heavy metals, may be at play at Vulcano confounding the CO2 effects.

Depth-specific fluctuations of gene expression and protein abundance modulate the photophysiology in the seagrass Posidonia oceanica.

Here we present the results of a multiple organizational level analysis conceived to identify acclimative/adaptive strategies exhibited by the seagrass Posidonia oceanica to the daily fluctuations in the light environment, at contrasting depths. We assessed changes in photophysiological parameters, leaf respiration, pigments, and protein and mRNA expression levels. The results show that the diel oscillations of P. oceanica photophysiological and respiratory responses were related to transcripts and proteins expression of the genes involved in those processes and that there was a response asynchrony between shallow and deep plants probably caused by the strong differences in the light environment. The photochemical pathway of energy use was more effective in shallow plants due to higher light availability, but these plants needed more investment in photoprotection and photorepair, requiring higher translation and protein synthesis than deep plants. The genetic differentiation between deep and shallow stands suggests the existence of locally adapted genotypes to contrasting light environments. The depth-specific diel rhythms of photosynthetic and respiratory processes, from molecular to physiological levels, must be considered in the management and conservation of these key coastal ecosystems.