# Adaptation of Temperate Seagrass to Arctic Light Relies on Seasonal Acclimatization of Carbon Capture and Metabolism

Alexander Jueterbock, Bernardo Duarte, James Coyer, Jeanine L. Olsen, Martina Elisabeth Luise Kopp, Irina Smolina, Sophie Arnaud-Haond, Zi-Min Hu, Galice Hoarau

<figure><img src="image_1.png" /></figure>

## To cite this version:

Alexander Jueterbock, Bernardo Duarte, James Coyer, Jeanine L. Olsen, Martina Elisabeth Luise Kopp, et al.. Adaptation of Temperate Seagrass to Arctic Light Relies on Seasonal Acclimatization of Carbon Capture and Metabolism. Frontiers in Plant Science, 2021, 12, pp.745855. 10.3389/fpls.2021.745855. [hal-03589731](https://hal.umontpellier.fr/hal-03589731)

HAL Id: hal-03589731

https://hal.umontpellier.fr/hal-03589731

Submitted on 22 Dec 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for
updates

# Adaptation of Temperate Seagrass to Arctic Light Relies on Seasonal Acclimatization of Carbon Capture and Metabolism

Alexander Jueterbock<fnref n="1" />*, Bernardo Duarte<fnref n="2" />,<fnref n="3" />, James Coyer⁴, Jeanine L. Olsen⁵,
Martina Elisabeth Luise Kopp⁶, Irina Smolina⁶, Sophie Arnaud-Haond⁷, Zi-Min Hu⁸ and
Galice Hoarau⁶

OPEN ACCESS

**Edited by:**

Kasper Elgetti Brodersen,
University of Copenhagen,
Denmark

**Reviewed by:**

Manoj Kumar,
Indian Institute of Toxicology
Research (CSIR), India
Matthew William Fraser,
University of Western Australia,
Australia

*Correspondence:

Alexander Jueterbock
Alexander-Jueterbock@web.de

**Specialty section:**

This article was submitted to
Functional Plant Ecology,
a section of the journal
Frontiers in Plant Science

**Received:** 22 July 2021

**Accepted:** 29 October 2021

**Published:** 02 December 2021

**Citation:**

Jueterbock A, Duarte B, Coyer J, Olsen JL, Kopp MEL, Smolina I, Arnaud-Haond S, Hu Z-M and Hoarau G (2021) Adaptation of Temperate Seagrass to Arctic Light Relies on Seasonal Acclimatization of Carbon Capture and Metabolism. *Front. Plant Sci.* 12:745855. doi: 10.3389/fpls.2021.745855

Due to rising global surface temperatures, Arctic habitats are becoming thermally suitable for temperate species. Whether a temperate species can immigrate into an ice-free Arctic depends on its ability to tolerate extreme seasonal fluctuations in daylength. Thus, understanding adaptations to polar light conditions can improve the realism of models predicting poleward range expansions in response to climate change. Plant adaptations to polar light have rarely been studied and remain unknown in seagrasses. If these ecosystem engineers can migrate polewards, seagrasses will enrich biodiversity, and carbon capture potential in shallow coastal regions of the Arctic. Eelgrass (Zostera marina) is the most widely distributed seagrass in the northern hemisphere. As the only seagrass species growing as far north as 70°N, it is the most likely candidate to first immigrate into an ice-free Arctic. Here, we describe seasonal (and diurnal) changes in photosynthetic characteristics, and in genome-wide gene expression patterns under strong annual fluctuations of daylength. We compared PAM measurements and RNA-seq data between two populations at the longest and shortest day of the year: (1) a Mediterranean population exposed to moderate annual fluctuations of 10–14 h daylength and (2) an Arctic population exposed to high annual fluctuations of 0–24 h daylength. Most of the gene expression specificities of the Arctic population were found in functions of the organelles (chloroplast and mitochondrion). In winter, Arctic eelgrass conserves energy by repressing respiration and reducing photosynthetic energy fluxes. Although light-reactions, and genes involved in carbon capture and carbon storage were upregulated in summer, enzymes involved in CO₂ fixation and chlorophyll-synthesis were upregulated in winter, suggesting that winter metabolism relies not only on stored energy resources but also on active use of dim light conditions. Eelgrass is unable to use excessive amounts of light during summer and

<footnote n="1">¹Algal and Microbial Biotechnology Division, Faculty of Biosciences and Aquaculture, Nord University, Bodo, Norway,</footnote>

<footnote n="2">²Marine and Environmental Sciences Centre, Faculty of Sciences of the University of Lisbon, Lisbon, Portugal,</footnote>

<footnote n="3">³Departamento de Biologia Vegetal da Faculdade de Ciências da Universidade de Lisboa, Lisbon, Portugal, ⁴Shoals Marine Laboratory, University of New Hampshire, Durham, NH, United States, ⁵Ecological Genetics-Genomics Group, Groningen Institute of Evolutionary Life Sciences, University of Groningen, Groningen, Netherlands, ⁶Marine Molecular Ecology Group, Faculty of Biosciences and Aquaculture, Nord University, Bodo, Norway, ⁷UMR MARBEC Marine Biodiversity Exploitation and Conservation, Univ Montpellier, CNRS, IFREMER, IRD, Montpellier, France, ⁸Ocean School, Yantai University, Yantai, China</footnote>
demonstrates a significant reduction in photosynthetic performance under long daylengths, possibly to prevent photoinhibition constrains. Our study identified key mechanisms that allow eelgrass to survive under Arctic light conditions and paves the way for experimental research to predict whether and up to which latitude eelgrass can potentially migrate polewards in response to climate change.

**Keywords:** eelgrass (*Zostera marina*), climate change, Arctic light, respiration, photosynthesis, carbon capture, daylength, energy storage

INTRODUCTION

Climate-change induced increases in surface temperatures are
causing a worldwide poleward shift of temperature isotherms
(McMahon and Hays, 2006) and concomitantly, distribution
limits (Wernberg et al., 2013; Melbourne-Thomas et al., 2020),
and spring phenology of species (Clausen and Clausen, 2013;
Poloczanska et al., 2013). Furthermore, polar regions are
warming at a particularly high rate (Høye et al., 2007; Burrows
et al., 2011). For example, Arctic sea ice has melted at a
rate of 6% per decade since the late 1970 (Cai et al., 2021),
opening thermally suitable habitats with unoccupied niches
available for colonizing species from the south (Cheung
et al., 2009).

Whether a temperate species can immigrate into ice-free
Arctic regions depends not only on finding thermally suitable
conditions but also on its ability to tolerate extreme Arctic
light conditions (Huffeldt, 2020). Constant daylight or darkness
can be particularly challenging for sessile autotrophs. To what
extent constant summer daylight inhibits cellular respiration
(Kok effect; Griffin and Heskel, 2013; Heskel et al., 2013),
growth, productivity, and photosynthesis in plants (Kallio and
Valanne, 1975) remains controversial and may depend on the
species (Sysoeva et al., 2010). A poleward range shift further
implies that growth and reproduction are triggered in suboptimal
seasons if temperate species follow the same photoperiodic
changes as specific life cycle triggers when immigrating into
the Arctic (Taulavuori et al., 2010). Yet, research and, thus,
knowledge on physiological adjustments of temperate plants
to Arctic light conditions is scarce in both terrestrial, and
marine systems (Huffeldt, 2020).

Seagrasses comprise a polyphyletic group of marine flowering
plants inhabiting shallow subtidal shores from subarctic to
tropical regions. Seagrass meadows provide habitat and food
that support a diverse food web, as well as important ecosystem
services such as nutrient fixation, protection of the shoreline
from erosion, nursery grounds for commercially important
fishes, and carbon sequestration of ca. 104 Tg C year⁻¹ (Duarte
and Krause-Jensen, 2017), or 1% of the global annual carbon
emissions (ca. 9.8 Gt C year⁻¹; Friedlingstein et al., 2020).
Thus, if seagrass can grow under Arctic light conditions, its
predicted northward shift (Wilson and Lotze, 2019) will likely
enrich biodiversity and introduce drastic functional changes
in future Arctic coastal ecosystems, food webs, fish/invertebrate
stocks, and biogeochemical cycles (Krause-Jensen and
Duarte, 2014).

The eelgrass *Zostera marina* is the most widely distributed seagrass species, occurring in suitable habitats along shores of all five continents in the northern hemisphere (Green and Short, 2004). Eelgrass meadows grow in coastal waters from subtropical regions at latitudes as low as 36°N (e.g., southern Iberia; Green and Short, 2004) to polar regions at latitudes as high as ca. 70°N (e.g., Greenland at 64.5° N; Olesen et al., 2015) and northern Norway at 69.8° N (Olsen et al., 2013).

Temperature, optimally from 15-25°C (Boström et al., 2004; Staehr and Borum, 2011; Zhang et al., 2016), is considered as the most important factor controlling phases of eelgrass sexual reproduction (Blok et al., 2018). Photoperiod appears to not control eelgrass reproduction because the eelgrass genome lacks the phytochrome C gene (PHYC; Olsen et al., 2016) that is involved in the photoperiodic control of flowering in most terrestrial plants (Chen et al., 2014; Woods et al., 2014). Accordingly, the species' northern distribution limit probably is limited by cool summer water temperatures <13.5°C that prevent seed maturation (Blok et al., 2018). As its seed germination is not regulated by photoperiodic rhythms (Moore et al., 1993), eelgrass seeds are likely to develop and disperse in a warming Arctic and contribute to a seed-driven northward shift of eelgrass meadows under projected climate change (Clausen et al., 2014; Blok et al., 2018). Temperature also is one of the most important factors for eelgrass growth (Clausen et al., 2014; Zhang et al., 2016). For example, increased investment in below-ground rhizomes contributes to larger annual biomass in populations from high latitudes compared with populations from warmer waters (Clausen et al., 2014).

Below-ground rhizomes provide stored energy to fuel eelgrass
metabolism under low-light conditions during Arctic winter
(Olesen and Sand-Jensen, 1993; Burke et al., 1996) and the
concentration of stored carbohydrates in eelgrass beds increases
with latitude (Colarusso, 2007; Novak et al., 2020). However,
rising temperatures may lead to a mismatch between peak
growth - advancing to winter darkness - and available solar
radiation in spring (Clausen et al., 2014). The mismatch, in
combination with an increased annual biomass turnover under
warmer temperatures (Olesen et al., 2015), requires eelgrass
to store sufficient resources to fuel winter metabolism and the
higher and earlier production of shoots at the onset of the
growing season (Clausen et al., 2014). Thus, poleward migrating
eelgrass will have strong implications for the dynamics of
below-ground carbon sinks.
The minimum light requirement for long-term survival and growth – meeting respiratory and metabolic energy demands – ranges between 11 and 34% surface irradiance (Ochieng et al., 2010), corresponding to 1.2–12.6 mol photons m⁻² day⁻¹ (Gattuso et al., 2006). Six hours daily of saturating irradiance (Hk) was suggested to be the minimum allowing for eelgrass growth and survival (Dennison and Alberte, 1985). Furthermore, shaded eelgrass showed significantly higher mortality rates under increased temperatures (16% at 20°C vs. 0.57% at 12°C; Eriander, 2017).

Most previous research has focused on the impact of modest light reductions or depth limits of eelgrass (e.g., Ralph et al., 2007), although a period of 30 days at 10% surface irradiance (<3 mol photons m⁻² d⁻¹) limited the survival of eelgrass transplants in estuarine temperate systems (Moore et al., 1997). The potential for poleward migration into the Arctic remains unknown because eelgrass survival under extreme light periods remains understudied. In order to investigate to what extent extreme Arctic photoperiods may present an energetic constraint to projected northward range shifts of eelgrass, we compared photosynthetic characteristics and gene expression patterns between two populations at the longest and shortest day of the year: (1) a Mediterranean population experiencing moderate annual fluctuations of 10–14 h daylength and 2) an Arctic population experiencing high annual fluctuations of 0–24 h daylength. The specific objectives were to:

1. Characterize energy-saving mechanisms (in photosynthetic and metabolic pathways) that may allow eelgrass to survive during long Arctic winter nights.

2. Determine possible adaptive mechanisms of Arctic eelgrass to compensate for energy loss during winter with increased investment in energy storing mechanisms during summer.

Differences in gene expression between summer and winter or day and night samples from a single population are interpreted as resulting from seasonal or diurnal acclimatization (Box 1). Our study cannot identify whether gene-expression or photosynthetic specificities of the Norwegian population (not replicated in the French population) are acclimatizations resulting from phenotypic plasticity or whether they are genetic adaptations (Box 1) that evolved under polar light regimes. Thus, in this manuscript, we use the term "adaptation" when referring to either acclimatization or genetic adaptation (Box 1).

**BOX 1** | Glossary of adaptation-related terms

**Acclimatization** – Phenotypic change over the life span of an organism induced by epigenetic regulation in response to environmental cues.

**Adaptation** – All plastic and genetically adaptive changes that reduce the negative effects of unfamiliar and undesirable external environment. This includes both acclimatization and genetic adaptation.

**Genetic adaptation** – Evolution of genetically determined heritable traits in a population across generations in response to natural selection.

# MATERIALS AND METHODS

## Site Characteristics

Eelgrass shoots were sampled in Norway (Røvik, 67°16′06.2″ N, 15°15′38.4″ E) and France (Sète, Thau Lagoon, 43°25′08.0″ N, 3°40′03.9″ E). The youngest two shoots of each of six plants were collected at each site at noon and the following midnight around summer solstice (June 21) and winter solstice (December 21) 2017. Intact shoots were placed in RNAlater (ThermoFisher) and stored after 1–2 days at –20°C.

Defining Zeitgeber Times ZTs (time after an environmental time cue that entrains circadian rhythms, such as sunlight) that correspond between the French and Norwegian samples is complicated by the fact that the two latitudes do not only differ in diurnal light rhythms but also in light quantity and quality (Taulavuori et al., 2010) (i.e., light conditions at dawn at the French site do not correspond to dawn light conditions at the Norwegian site and, thus may mark different ZTs in the circadian cycle). We chose midday and midnight as time points to compare samples from the two locations as these represent times at which the plants from both locations had the potential to receive most/least light of the day in both seasons, and the time that is at equal temporal distance from sunset and sunrise. This implies, however, that differences between the Norwegian and French samples at corresponding time points may be partly explained by them being at different ZTs within their circadian rhythm. However, differences in ZTs are less affecting seasonal differences between the Norwegian and French, because we compared daily averages that involved measurements from both midday and midnight samplings.

Snapshots of temperature (°C) and light (Lux; Table 1) were taken at each sampling event in ca. 20 cm water depth with Pendant data loggers (Onset HOBO UA-002). Annual variation of temperature (°C), salinity, and Lux conditions at the Norwegian sampling site were recorded from spring 2016 to spring 2017 with data loggers (HOBO U24-002-C for salinity/temperature with 30 min intervals, and HOBO UA-002 for light/temperature with 10 min intervals) anchored ca. 30 cm above sandy substrate during low tide. At low tide, the light/temperature loggers were falling dry. Temperature (°C) and salinity conditions in the Thau Lagoon (43° 22′ 44.4″ N, 3° 34′ 17.4″ E) were downloaded for year 2016 from the marine environmental station database of Thau Lagoon (Messiaen et al., 2018) on July 31, 2020, which were recorded every 30 min with data loggers (WTW) at 1–2 m water depth. Light (PAR µmol m⁻² s⁻¹) conditions were recorded with dataloggers (SKUW 215 PAR Sensor, Skye Instruments) every 30 min for December 2015–July 2016 at the Thau Lagoon Observatory (43°24′53″ N, 3°41′16″ E; Trombetta et al., 2019) and were downloaded from Sea scientific open data publication (www.seanoe.org, accessed on July 31, 2020). To make the light conditions comparable between the Norwegian and French sampling sites, Lux values were converted to approximate PAR (photosynthetically active radiation) values by applying a multiplication factor of 0.019, an approximation for the spectral power distribution of sunlight (6,000 K; Creswell, 2010).
<table><caption>**TABLE 1** | Snapshots of temperature and light conditions measured at the sampling events.</caption><thead><tr><th rowspan="2">Measurement</th><th rowspan="2">Population</th><th colspan="2">June</th><th colspan="2">December</th></tr><tr><th>Night</th><th>Day</th><th>Night</th><th>Day</th></tr></thead><tbody><tr><td rowspan="2">Temperature (°C)</td><td>Norway (Røvik)</td><td>15.1</td><td>16.0</td><td>-1.9</td><td>-1.6</td></tr><tr><td>France (Sète)</td><td>22.8</td><td>21.7</td><td>7.3</td><td>8.2</td></tr><tr><td rowspan="2">Light (µmol m<sup>-2</sup> s<sup>-1</sup>)</td><td>Norway (Røvik)</td><td>6.3</td><td>235.6</td><td>0.0</td><td>0.0</td></tr><tr><td>France (Sète)</td><td>0.0</td><td>327.2</td><td>0.2</td><td>535.4</td></tr></tbody></table>

Seasonal variation in light levels was greater in Norway, where PAR levels remained close to zero from mid-November to mid-January, but reached ca. 1,000 µmol m⁻² s⁻¹ in summer, similar to France (Figures 1C,D,E; Supplementary Table 1). In summer, light levels showed a diurnal rhythm in both Norway and France, but – as expected – reached slightly higher nocturnal levels in Norway.

Summer temperatures were ca. 10°C higher in France (20–25°C) compared to Norway (10–15°C; Supplementary Figure 1; Supplementary Table 1). Stronger Norwegian diurnal temperature variations likely reflected the differences between air and water temperatures around spring tides when the dataloggers may have been shortly exposed to air during low-tides but submerged during high-tides. In contrast, the eelgrass meadow in the Thau lagoon was constantly submerged. Accordingly, the data loggers revealed relatively constant salinity throughout the year at ca. 37 PSU (Supplementary Figure 2; Supplementary Table 6), while the high fluctuation (10–30 PSU) in Norway was due to seasonal rainfall, freshwater runoff, ice formation, and melting.

## Photosynthesis

At each sampling event, we assessed leaf photochemistry using a Chlorophyll a Pulse Amplitude Modulated (PAM) Fluorometer (FluorPen FP100, Photon Systems Instruments), by measuring the increase in chlorophyll a fluorescence upon illumination after a 10-min dark period induced by leaf-clips (OJIP curve; Bussotti et al., 2010; Duarte et al., 2017), using a saturating pulse at 455 nm from 10 seagrass shoots.

When a dark-adapted leaf is illuminated with a saturating light intensity, it exhibits a so-called OJIP curve from a state where all reaction centers are opened (all quinone acceptors are oxidized and can accept electrons) to a state where all quinone acceptors are reduced and chlorophyll fluorescence reaches its maximum. The rise in fluorescence transient along this polyphasic OJIP curve depicts the rate of reduction kinetics of various components of Photosystem II (PSII). The OJIP curve provides information on the processes involved in photon absorption at the PSII antennae (ABS), transformation into an exciton and its free energy flux (RE), and transduction of its energy in the form of electronic potential via electron transport (ET) from the donor side of the PSII to the electron acceptor of the PSI. The analysis of this OJIP transient (JIP test) allows to estimate the energy fluxes, involving absorption, entrapment in the PSII, transport to the quinone pool, and dissipation in the form of heat or fluorescence that helps to

avoid photoinhibition damages. The formulae and definitions of terms derived from the JIP-test (chlorophyll a induction curves, or Kautsky curves measurements; Supplementary Table 2) are listed in Table 2.

Differences in these photosynthetic characteristics between populations (day/night/season) were tested for significance with Wilcoxon rank sum tests followed by correcting the values of *p* for multiple pairwise comparisons using the Benjamini-Hochberg method with the R package "rstatix" v0.6.0 (Kassambara, 2020). Multivariate differences in photochemistry were estimated from photochemical raw data using a canonical analysis of principal coordinates (CAP) with Euclidean distances with the Primer 6 software (Clarke and Gorley, 2005). This multivariate approach is insensitive to heterogeneous data and frequently used to compare different sample groups (Duarte et al., 2017; Cruz de Carvalho et al., 2020). Differences in the overall photochemical metabolism were estimated from normalized photochemical raw data with a canonical analysis of principal coordinates (Duarte et al., 2017). Multivariate differences were visualized as Euclidean distances along the first two principal coordinates.

## RNA Preparation

RNA was extracted from small pieces of stored samples (RNAlater; ThermoFisher) that were homogenized (PowerLyzer 24; Qiagen) and processed according to the standard protocol of the RNeasy PowerPlant Kit (Qiagen). We consistently excised the mid portion of the youngest fully mature leaf in the shoot (usually the second-rank leaf), so that tissue-age dependent within- or among-leaf variations were unlikely to explain differences in gene expression (Ruocco et al., 2019a,b). RNA concentrations were measured with Qubit™ RNA HS Assays (ThermoFisher) using a Qubit 4 fluorometer. RNA integrity was verified on TapeStation 2200 with RNA ScreenTape. RIN values were >3.6.

RNA libraries were prepared from 1 µg of extracted RNA with NEBNext® Ultra™ II Directional RNA Library Prep Kit for Illumina (New England Biolabs) following the polyA mRNA workflow of the kit protocol. Libraries were verified on TapeStation 2200 using High Sensitivity D1000 ScreenTape. Paired-end (2×150 bp) RNA sequencing (RNASeq) data were generated using a NextSeq 500/550 High Output Kit v2 (Illumina) on an Illumina NextSeq 500 platform at Nord University.

## Sequence Data Analysis

Raw reads were trimmed for low quality (Phred score Q<20, 99% base call accuracy), and TruSeq adapter sequences
A

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

C

<figure><img src="image_6.png" /></figure>

E

<figure><img src="image_7.png" /></figure>

Location • Rovika • Thau Lagoon Observatory

B

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

D

<figure><img src="image_10.png" /></figure>

**FIGURE 1** | Seasonal and diurnal gene expression changes of the Norwegian and French population of the seagrass *Zostera marina* in relation to differences in annual light conditions at the sampling sites. Number of genes upregulated in summer, and winter (A) in France and in Norway, or upregulated during day, and night (B), separately for summer and winter. The subset of genes in A and B are independent from each other but partly overlapping. Light levels in June (C), December (D), and over the entire year (E) at the Norwegian sampling site (blue) and at an observatory station in the Thau Lagoon close to the French sampling site (red). In the Thau Lagoon, values are missing from July to December. PAR: Photosynthetically active radiation.

(overlap > 3 bp) before reads <25 bp were removed with TrimGalore! v 0.6.0.<fnref n="1" /> Read quality was checked with FastQC v 0.11.8 to control for aberrant read base content, length distribution, duplication and over-representation.<fnref n="2" />

Trimmed reads were mapped to the genome of *Z. marina* (Olsen et al., 2016) v 2.2 from ORCAE (Sterck et al., 2012) and the organellar genomes (chloroplast and mitochondria) of *Z. marina* v2.1 from ORCAE (Sterck et al., 2012) with the splice-aware RNA-seq aligner STAR v 2.7.2a (sjdbOverhang 149; Dobin et al., 2013, 2016), guided by gff3 gene annotations v2.2 from ORCAE (Sterck et al., 2012), converted to gtf with

<footnote n="1">¹https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/</footnote>

<footnote n="2">²http://www.bioinformatics.babraham.ac.uk/projects/fastqc/</footnote>
<table><caption>**TABLE 2** | Formulae and definitions of terms derived from the JIP test.</caption><thead><tr><th>Term</th><th>Formula</th><th>Definition</th></tr></thead><tbody><tr><td>PI<sub>ABS</sub></td><td>[γ<sub>RC</sub>'(1 − γRC)] × [φ<sub>Po</sub><sup>2</sup>/(1 − φ<sub>Po</sub>)] × [ψ<sub>o</sub><sup>3</sup>/(1 − ψ<sub>o</sub>)]</td><td>Performance Index reflecting the overall photosynthetic performance by combining three parameters: (1) the density of reaction centers, (2) the electron transport at the onset of illumination, and (3) the maximum energy flux reaching the reaction center in PS II.</td></tr><tr><td>ABS/CS</td><td>F<sub>o</sub><sup>4</sup></td><td>Absorbed energy flux measured as the minimal fluorescence intensity in a dark adapted frond when all reaction centers are opened (all quinone acceptors are oxidized and can accept electrons)</td></tr><tr><td>ET/CS</td><td>φ<sub>Po</sub> × ψ<sub>o</sub> × (ABS/CS)</td><td>Electron transport energy flux</td></tr><tr><td>TR<sub>o</sub>/CS</td><td>φ<sub>Po</sub> × (ABS/CS)</td><td>Trapped energy flux</td></tr><tr><td>DI<sub>o</sub>/CS</td><td>(ABS/CS) − (TR<sub>o</sub>/CS)</td><td>Dissipated energy flux</td></tr><tr><td>Area</td><td>Extracted from the recorded OJIP</td><td>Size of the oxidized quinone pool; total complementary area between the fluorescence induction curve and F = F<sub>M</sub></td></tr><tr><td>RC<sub>o</sub>/CS</td><td>φ<sub>Po</sub> × (V<sub>o</sub><sup>5</sup>/M<sub>o</sub><sup>6</sup>) × ABS/CS</td><td>Density of photosystem II oxidized reaction centers per excited cross-section</td></tr><tr><td>TR<sub>o</sub>/DI<sub>o</sub></td><td></td><td>Contribution of the light reactions for primary photochemistry</td></tr><tr><td>Ψ<sub>o</sub>/(1 − Ψ<sub>o</sub>)</td><td></td><td>Contribution of the dark reactions for primary photochemistry</td></tr></tbody></table>

the program gffread.<fnref n="3" /> Alignments that contained noncanonical splice junctions were filtered out. Duplication rates were checked with the R package “dupRadar” v1.14.0 (Sayols et al., 2016).

Read counts for 21,069 genes (Supplementary Table 3) were obtained with featureCounts in the R package “Rsubread” v1.34.6 (Liao et al., 2019) with the non-default settings minOverlap = 3, MinMQS = 20, strandSpecific = 2, isPairedEnd = TRUE. Genes of low expression (library average < 5, n = 7,023) were excluded from subsequent analyses to avoid potential artefacts from sequencing errors.

Gene IDs of the Z. marina genome v2.2 (Olsen et al., 2016) were annotated with gene descriptions, Gene Ontology (GO) terms, KEGG Orthology groups (KO), and Enzyme Commission (EC) numbers (based on enzyme nomenclature) from the ORCAE database (Sterck et al., 2012). Gene ontology terms assign genes to one or more categories related to the molecular function, cellular component or biological process of the gene products. KEGG Orthology groups assign gene products to their function in a biological pathway of molecular interactions. The EC number classifies genes numerically based on the chemical reactions they catalyze.

## Multivariate Clustering of Gene Expression

To characterize overall transcriptomic differentiation, we clustered the samples hierarchically by the first five principal components of gene coverage (14,046 genes), which together explained 73.4% of the variation, using the R package “FactoMineR” v. 1.42 (Lê et al., 2008). For Principal Component Analysis (PCA), the scale.unit was set to FALSE not to scale the expression values to unit variance. Prior to PCA, we regularized log-transformed gene count values (Supplementary Table 4) with DESeq2 v 1.24.0 (Love et al., 2014) in order to assure homoscedasticity and to account for differences in sequencing depth.

## Differential Gene Expression

We tested for differential gene expression between summer and winter samples (day and night samples combined) and between day and night samples in France and Norway, respectively, with DESeq2 v1.24.0 (Love et al., 2014). The values of *p* were adjusted with Independent Hypothesis Filtering based on the R package “IHW” v1.14.0 (Ignatiadis et al., 2016). We plotted overlaps of differentially expressed genes between seasons and locations with the R package “eulerr” v6.1.0 (Larsson, 2019) and “nVennR” v0.2.1 (Pérez-Silva et al., 2018).

For each population, we tested (Fisher's exact test) for enriched gene ontology terms (biological process, cellular component, and molecular function) in genes that were differentially expressed between day and night samples or between summer and winter samples with the R package “topgo” (Alexa and Rahnenführer, 2010). Significance levels were adjusted with the Benjamini and Hochberg correction (Benjamini and Hochberg, 1995), using the p adjust function in R (R Core Team, 2019), to control for the false discovery rate when conducting multiple comparisons.

The adjusted-log10 of the values of *p* of significantly enriched GO-terms were visualized for each population in heatmaps using the heatmap.2 function of the R package “gplots” v.3.0.3 (Warnes et al., 2020). GO terms were hierarchically clustered

<footnote n="1">¹γ<sub>RC</sub> = Chl<sub>o</sub>/Chl<sub>total</sub> = RC/(ABS + RC), the probability that a PSII Chl molecule functions as reaction center (RC).</footnote>

<footnote n="2">²φ<sub>Po</sub> = TR<sub>o</sub>/ABS = [1 − (F<sub>M</sub>/F<sub>o</sub>)], the maximum quantum yield for primary photochemistry, with F<sub>M</sub> representing the maximum fluorescence intensity when all reaction centers are closed (all quinone acceptors are reduced).</footnote>

<footnote n="3">³ψ<sub>o</sub> = ET<sub>o</sub>/TR<sub>o</sub> = (1 − V<sub>o</sub>), efficiency/probability for electron transport (ET), or the probability that an electron moves further than Q<sub>h</sub> (the primary quinone acceptor).</footnote>

<footnote n="4">⁴F<sub>o</sub>, minimal fluorescence when all PSII RCs are open.</footnote>

<footnote n="5">⁵V<sub>o</sub> = (F<sub>M</sub> − F<sub>o</sub>)/(F<sub>M</sub> − F<sub>o</sub>), with F<sub>M</sub> representing the fluorescence intensity at the J-step (at 2 ms).</footnote>

<footnote n="6">⁶M<sub>o</sub>, Slope of the origin of the fluorescence rise. Maximal rate of accumulation of the fraction of closed reaction centers.</footnote>

<footnote>Based on Zhu et al. (2005) and Duarte et al., (2017).</footnote>

<footnote n="3">³https://github.com/gpertea/gffread</footnote>
with the R package “goSTAG” v1.12.0 (Bennett and Bushel, 2017) based on euclidean distances and Ward's minimum variance clustering method with Ward's (Ward, 1963) clustering criterion implemented (ward.D2 option; Murtagh and Legendre, 2014). The sample dendrogram was based on Euclidean distances and a complete linkage clustering method.

Enriched/perturbed KEGG metabolic pathways in the fold-change matrices of differentially expressed genes were identified using the R package “gage” v2.37.0 (Luo et al., 2009). Nuclear genes were mapped to KO (Kegg Orthology)-terms by translating EC numbers from the genome annotation to corresponding KO-terms based on the KEGG API list http://rest.KEGG.jp/link/ec/ko. For organellar genes, EC numbers and, thus KO-terms were not annotated. For pathways perturbed (*p* < 0.05) in at least one of the populations and seasons, we plotted the significant fold-changes of molecular functions (KO-terms) for all populations and seasons with the R package “pathview” v3.11 (Luo and Brouwer, 2013). When several genes referred to the same KO-term, we plotted their average fold-change in the respective pathway.

# RESULTS

## Photosynthesis

Overall photosynthetic performance, represented by the photochemical performance index ($P_I_{ABS}$; **Figure 2A**) showed inverse seasonal shifts in Norway and France. While the French population performed significantly better during summer, the Norwegian population performed significantly better during winter. The reduced summer performance of the Norwegian population was associated with peculiar photosynthetic characteristics. In contrast to the French population and the Norwegian winter samples, light harvesting ability increased in the Norwegian summer samples under non-saturating light levels during night-time, as represented by the increase in absorbed (ABS/CS; **Figure 2C**) and trapped ($TR_0/CS$; **Figure 2D**) energy flux, and in the number of oxidized reaction centers per leaf cross section ($RC_0/CS$; **Figure 2G**). The contribution of dark reactions to the overall photochemistry [$\Psi_0$/(1-$\Psi_0$); **Figure 2I**] decreased significantly during summer nights. However, despite readiness of the electron transport chain (ETC) for electron transport, as indicated by the significant increase in the size of the oxidized quinone pool during Norwegian summer nights (Area; **Figure 2B**), the electron transport energy flux (ET₀/CS; **Figure 2E**) did not increase significantly. During winter, despite a significant increase in the contribution of light and dark reactions to the primary photochemistry [$TR_0/DI_0$, and $\Psi_0$/(1-$\Psi_0$); **Figures 2H,I**], the Norwegian population showed a significant reduction in all energy fluxes, particularly during night-time (ABS/CS, $TR_0$/CS, and ET₀/CS; **Figures 2C–E**), including dissipated energy loss ($DI_0/CS$; **Figure 2F**).

The samples clustered well by sampling location and daylength (78 of 80 samples correctly identified) along the first two principal components showing that the photochemical traits from the chlorophyll *a* induction curves could well classify

the samples (**Figure 3**). The Norwegian summer samples grouped into one cluster, irrespective of whether they were collected at day or night. The Norwegian winter samples collected during nighttime formed another distinct cluster. Winter daytime samples from France and Norway grouped close together, indicating that winter had a prevalent effect over location.

## Gene Expression

On average, 76% of 1 million raw reads (average per library) mapped to exons (**Supplementary Table 5**). Based on the expression of 13,932 genes covered >5 times, the samples grouped into four clusters, separated by location (France/Norway) and season (Summer/Winter; **Figure 4**). The most distinct gene expression patterns were recorded in the winter samples from Norway (NW). Gene expression was most similar between winter samples from France (FS) and summer samples from Norway (NS; **Figure 4**). Strictly distinct clusters of gene expression between day and night were only observed in France during summer (FS).

Seasonal differences in gene expression were found in 8652 genes in the French population and in 8569 genes in the Norwegian population (**Supplementary Table 6**). In summer, 4,274 genes were upregulated in France, and 4,084 genes in Norway, of which 1,475 overlapped between the two populations (**Figure 1A**). In winter, 2,691 genes were upregulated in France, and 2,798 genes in Norway, of which 1,687 overlapped between the two populations (**Figure 1A**).

Diurnal differences in expression were strongest (>3,000 genes) in, and most similar (>500 genes shared) between the French summer samples and the Norwegian winter samples (**Figure 1B**). In contrast, the number of genes with diurnal expression differences was an order of magnitude lower (ca. 600 genes) in the French winter samples or the Norwegian summer samples (**Figure 1B**). All genes with diurnal differences in gene expression are characterized in **Supplementary Table 7**.

## Functional Changes

Genes involved in photosynthesis, carbon fixation (dark reactions of photosynthesis; **Supplementary Figure 3**), and sugar metabolism (fructose and mannose metabolism: **Supplementary Figure 4** and glycolysis: **Supplementary Figure 5**) showed contrasting seasonal patterns between the Norwegian and French seagrass populations. Those are likely to explain the functional basis of the strong transcriptome similarity between Norwegian winter samples and French summer samples (**Figure 4**). Summer upregulation/winter downregulation of most of the organellar functions, both related to the chloroplast and the mitochondrion, characterizes the functional transcriptome specificities of the Norwegian samples (**Figure 5**). All pathways enriched for seasonal or diurnal differentially expressed genes are listed in **Supplementary Table 8**.

## Photosynthesis and Carbon Sequestration

The Norwegian samples upregulated in summer genes related to photosynthesis and ATP metabolism (**Figure 6**;
<figure><img src="image_11.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>C</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>D</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>E</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>F</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>G</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>H</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>I</figcaption></figure>

<figure><img src="image_19.png" /></figure>

**FIGURE 2** | Photosynthetic characteristics. Snapshots of photosynthetic characteristics with standard error bars and indications for significant differences between populations, day and night, or summer and winter. Significance codes after the value of *p* correction (Benjamini-Hochberg) for multiple comparisons: "<0.05"; *, "<0.01"; **, "<0.001"; ***, "<0.0001"; ****. (A) Index of photosynthetic performance (P<sub>ABS</sub>). (B) Size of the oxidized quinone pool (Area). (C) Absorption per excited cross-section (ABS/CS). (D) Trapping per excited cross-section (TR<sub>o</sub>/CS). (E) Electron transport per excited cross-section (ET<sub>o</sub>/CS). (F) Dissipated energy flux per excited cross section (DI<sub>o</sub>/CS). (G) Oxidized reaction centers per excited cross-section (RC<sub>o</sub>/CS). (H) Contribution of the light reactions to the primary photochemistry (TR<sub>o</sub>/DI<sub>o</sub>). (I) Contribution of the dark reactions to the primary photochemistry [Ψ<sub>o</sub>/(1-Ψ<sub>o</sub>)].

<figure><img src="image_20.png" /></figure>

**FIGURE 3** | Multivariate differences among samples in photochemical data, visualized as Euclidean distances along the first two principal coordinates.

**Supplementary Figure 6), photosystem, and photosynthetic membrane (Supplementary Figure 6), especially during the day (Figure 7; Supplementary Figures 7,8). Moreover, the Norwegian summer samples showed enrichment (*p* < 0.05) of genes with significant upregulation in the fructose and mannose metabolism (ko00051; **Supplementary Figure 4**). While significantly lower expressed in winter, genes involved in photosynthesis (Figure 7), chlorophyll binding (Supplementary Figure 7), chloroplast, photosystem, and photosynthetic membrane (Supplementary Figure 8) were significantly upregulated in winter nights as compared with winter days. During summer nights, the Norwegian samples showed the upregulation of trehalose biosynthesis (Figure 7).**

The carbon fixation pathway in photosynthetic organisms (ko00710; **Supplementary Figure 3**) was enriched (*p* < 0.05) for genes that were upregulated in the Norwegian samples during summer but in the French samples during winter. The French samples upregulated during winter days genes involved in oxidoreduction (**Figure 7; Supplementary Figure 7**), as well as chloroplast genes
<figure><img src="image_21.png" /><figcaption>**FIGURE 4** | Hierarchical cluster of all 48 samples based on the first five principal components of expression in 13,932 genes. FS: France Summer, FW: France Winter, NS: Norway Summer, NW: Norway Winter.</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>**FIGURE 5** | Highly simplified overview of various functions that characterize the specificities in seasonal variation of gene expression in the Norwegian population of *Zostera marina*. Functions, enzymes and components upregulated in winter are highlighted with a blue background; such upregulated in summer are highlighted with an orange background.</figcaption></figure>

(Supplementary Figure 8), and were enriched for upregulated genes in the pentose phosphate pathway (ko00030; Supplementary Figure 9).

## Respiration

The French samples upregulated in winter genes related to respiration, including GO-terms of cellular respiration (Figure 6), mitochondria (Supplementary Figure 6), and oxidoreductase activity (Supplementary Figure 10). Moreover, the French winter samples were enriched for genes involved in glycolysis/gluconeogenesis (ko00010; Supplementary Figure 5).

In contrast, the Norwegian upregulated genes related to mitochondria (Supplementary Figure 6) and showed enrichment (*p* < 0.05) of genes with significant upregulation in glycolysis/gluconeogenesis (ko00010; Supplementary Figure 5) in summer. Only genes involved in the mitochondrial inner membrane were upregulated in the Norwegian samples during winter (Supplementary Figure 8).

## Peptide Synthesis

In winter, the Norwegian samples upregulated functions and processes involved in gene expression, including peptide synthesis, translation, RNA processing (Figure 6), binding of nucleic acids and catalytic activities on RNA and tRNA (Supplementary Figure 10), as well as genes involved in intracellular organelles such as the ribosome (Supplementary Figure 6). During winter days, the Norwegian samples had upregulated genes involved in protein maturation, peptidyl-lysine modification (Figure 7), NAD⁺ kinase, and phosphatidyl-inositol phospholipase C activity (Supplementary Figure 7).

The French samples showed diurnal changes in the expression of genes related to protein synthesis during summer, when genes involved in translation, amid biosynthesis, (Figure 7), as well as the ribosome (Supplementary Figure 8) were upregulated during the day.

## DISCUSSION

Our results provide the first assessment of seasonal (and diurnal) changes in photosynthetic and genome-wide gene expression
<figure><img src="image_23.png" /><figcaption>**FIGURE 6** | Seasonal changes in biological processes. Enriched processes are shown for each population with positive-log10 (*p*) for genes upregulated in summer (orange) and with negative-log10 (*p*) for genes upregulated in winter (cyan). The hierarchical cluster of GO terms is based on Euclidean distance and a ward.D2 clustering algorithm (Ward's minimum variance method).</figcaption></figure>

patterns specific to a population growing above the Arctic circle that contrast patterns observed in a Mediterranean (control) population. Our results indicate that adaptation to polar light regimes involves continued photosynthesis, carbon fixation and energy conservation during long winter nights, and carbon storage during the polar summer.

## Strategies to Survive Long Winter Nights
Photosynthesis and Carbon Fixation in Winter

Despite reduced rhythmicity in daylight levels during the Arctic winter, our Norwegian winter samples appeared to still follow a day-night rhythm in gene expression (Figure 1B). Indeed, the Norwegian population showed, in contrast to the French
<figure><img src="image_24.png" /></figure>

population, a far higher number of genes with diurnal expression differences during winter (3,099; **Figure 1B**) than during summer (657; **Figure 1B**). Increased night-time expression of genes related to the synthesis of chloroplasts and chlorophyll (**Figure 7**) is in line with earlier field observations of enhanced chlorophyll in winter acclimatized eelgrass (Thom et al., 2008), and suggests

that the Norwegian population actively synthesizes photosynthesis-relevant proteins to capture the smallest amounts of photosynthetically active radiation in winter.

Efficient use of very dim Arctic winter light in the Norwegian samples is suggested by a significant increase in the photochemical performance index ($PI_{ABS}$; **Figure 2A**), and in the contribution
of both light and dark reactions (Figures 2H,I) to the primary photochemistry during winter. Adjustments that may contribute to photosynthesis under low light levels include an increase in chlorophyll *b* levels in comparison to chlorophyll *a*, thus increasing photosystem II antennae size, and maximizing the photon harvesting probability (Ort et al., 2011) with a simultaneous increase (by chlorophyll *b*) in the affinity for harvesting light of long wavelengths that protrudes to greater water depths (Lee et al., 2007; Eriander, 2017). Moreover, under low light levels, eelgrass plants invest more energy into vertical leaf growth, compensating the low light levels with large photosynthetic surfaces and a reduction in canopy depth (Ochieng et al., 2010; Wong et al., 2021). While energy fluxes (ABS/CS and TR₀/CS; Figures 2C,D) during the Norwegian winter cannot match fluxes occurring in summer, they are likely to prevent the depletion of assembled carbohydrate storages. Accordingly, winter survival of eelgrass at its lower depth distribution in the Great Bay Estuary (New Hampshire, United States), depended primarily on adaptations that allowed for continued photosynthesis under extreme low light conditions (e.g., increase in leaf size, chlorophyll *a*, and maximum quantum yield of PSII; Chapter IV in Ochieng, 2008). Thus, adaptation to the Arctic appears to involve seasonal acclimatization that allows to actively use the smallest amounts of winter light *via* photosynthesis.

Accordingly, the Norwegian population appears to present the physiological potential to capture carbon in winter. For example, ribulose-bisphosphate carboxylase (RuBisCO, K01601), the initial carboxylase of the Calvin cycle that fixes atmospheric CO₂, was expressed at higher levels during winter than during summer (Figure 5; Supplementary Figure 3). This contrasts with 30% under-expression of RuBisCO in *Posidonia oceanica* leaves adapted to low light-levels (Mazzuca et al., 2009) and shows that winter gene expression in an Arctic population is not comparable with gene expression characteristics of seagrass experiencing low light levels in turbid waters year-round. In addition, the enzyme RuBisCO activase, although showing significantly higher expression in the Norwegian population during summer (Zosma13g01130, Supplementary Table 6), was still expressed during winter (Zosma294g00110, Zosma13g01130, and Zosma88g00030 in Supplementary Table 7). The presence of this activator, which removes the RuBisCO inhibitor CAP (2-carboxylarabinitol-1,5-bisphosphate), suggests that the Norwegian population could indeed maintain carbon fixation during the dark winter period. If increased RuBisCO expression in Norwegian winter samples is related to cold-acclimation, as shown in maize (Salesse-Smith et al., 2020), then warmer winters could hinder carbon assimilation and, thus, survival in the future Arctic.

Furthermore, the Norwegian winter samples had upregulated the enzyme phosphoenolpyruvate carboxylase (PEPC, K01595; Figure 5), which can catalyze CO₂ fixation with phosphoenolpyruvate to produce oxaloacetate and inorganic phosphate. This suggests that *Z. marina*, although being characterized as a C3 plant that fixes CO₂ primarily *via* RuBisCo (Koch et al., 2013), can fix some carbon during winter independent from photosynthesis, although at much lower levels

than RuBisCo would allow (Raven and Farquhar, 1990). PEPC CO₂ fixation is not uncommon in C3 plants (Chi et al., 2014), given that enzymes involved in the C4 photosynthetic biochemistry have functions in C3 plants other than in C4 photosynthesis, including the supply of carbon skeletons (*via* the intermediates oxaloacetate and malate) to the citric acid cycle (Lepiniec et al., 1994; Aubry et al., 2011; Figure 5) for the formation of new tissue. Thus, we suggest that PEP-carboxylase may present a key enzyme that allows for winter growth by decoupling photosynthesis and biomass formation not only in polar algae (Wiencke et al., 2007, 2009) but also in Arctic eelgrass.

Active growth in winter would explain the presence of freshly grown seagrass leaves in the dark season in the Norwegian population (Supplementary Figure 11). Eelgrass reaches light saturation between ca. 70 and 900 µmol m⁻² s⁻¹, and compensation above ca. 1–36 µmol m⁻² s⁻¹ (Ochieng, 2008). Given that photosynthetically active radiation exceeded 10 µmol m⁻² s⁻¹ in the Norwegian population only for short mid-day periods on two December days in 2016 (Supplementary Table 1), photosynthesis-driven carbon dioxide assimilation is expected to be very low in the Norwegian winter.

Winter growth, therefore, is expected as a result of light-independent carbon fixation and the usage of stored carbohydrates mobilized from below-ground tissues (Burke et al., 1996). For example, strong activation of sucrose synthase, of which the products can drive many metabolic pathways, has been shown in sink-tissues of eelgrass under negative carbon balances (Alcoverro et al., 1999). Sustained growth during periods of negative carbon balance (respiratory C loss exceeding photosynthetic C gain) – fueled by mobilized carbohydrate storages – has been shown in both eelgrass (Wium-Andersen and Borum, 1984; Olesen and Sand-Jensen, 1993), and the Mediterranean seagrass *P. oceanica* (Alcoverro et al., 2001). In conclusion, our results suggest that the use of low light levels and light-independent carbon fixation count among the adaptations that allow for eelgrass to grow during winter without fully depleting their stored carbohydrates.

## Energy Conservation Strategies in Winter

The Norwegian population appears to conserve energy during winter, as indicated by a significant reduction in all photosynthetic energy fluxes (ABS/CS, TR₀/CS and ET₀/CS; Figures 2C–E), resulting in a significantly reduced loss of energy *via* dissipation (DI₀/CS; Figure 2F). Such state of metabolic dormancy, limiting the enzymatic reactions of photosynthesis, may be adaptive to the combination of low light availability and low winter temperatures.

Moreover, the respiratory breakdown of carbohydrates must have been reduced in the Norwegian population during winter, given that genes related to glycolysis, ATP metabolism (Figure 6; Supplementary Figure 5), and the fructose and mannose metabolism (Supplementary Figure 4; e.g., fructose-bisphosphate-aldolase K01623; Figures 6, 7) were downregulated. Accordingly, while upregulated during winter in the French population (Figure 6), genes involved in cellular respiration were significantly downregulated in the Norwegian population during winter (e.g.,
malate dehydrogenase, pyruvate dehydrogenase, and Cytochrome c oxidase subunit 5b; **Supplementary Table 6; Figure 6**), suggesting a shutdown of catabolic pathways during the Norwegian winter. Thus, our results provide the first molecular support for winter-acclimatization in eelgrass that explains an increase in photosynthesis-to-respiration ratios, in line with reports of a greater net primary productivity in eelgrass acclimatized to winter than to summer conditions (Thom et al., 2008). Acclimatization of eelgrass to low light conditions, further involves reduced rhizome and root growth, lower production of new shoots through lateral branching, shoot density, biomass, and sheath length (Ochieng et al., 2010; Wong et al., 2021), all of which reduce the respiratory burden. Furthermore, active reduction in respiration rates likely accompanies other factors, such as increased chlorophyll a content and lower epiphytic growth, in photo-adaptive processes that are involved in lowering light-saturation and compensation points in winter-acclimatized eelgrass (Ochieng, 2008; Thom et al., 2008). Such seasonal acclimatization may to some degree counteract the respiratory burden of rising temperatures by decreasing respiration rates and critical light requirements (Zimmerman et al., 1989; Staehr and Borum, 2011).

Glycolysis (Figure 5; Supplementary Figure 5) is a catabolic pathway providing energy in the form of NADH, ATP, and pyruvate, the latter being further catabolized in the citric acid cycle as part of oxidative respiration. Its downregulation during winter in this study contrasts with upregulation of genes involved in carbohydrate catabolism (e.g., 1-fructose-bisphosphate aldolase) in low-light adapted seagrass leaves of the species *P. oceanica* (Mazzuca et al., 2009). Downregulation in our study may be explained by the ability of our low-light (winter) eelgrass to build up carbon storages during summer and, thus, to obtain carbon-storages during winter from storage tissues such as roots and rhizomes, instead of from the above-ground leaf tissue. In other words, while carbon catabolism was low in the seagrass leaves, we cannot know what processes were active in the below-ground organs, in which sucrose and starch can be depleted by 50 and 23%, respectively, after a 3-week shading period (Eriander, 2017).

That carbon storages were built up over summer was suggested by the significantly increased expression of enzymes involved in carbon fixation (Supplementary Figure 3) and by significantly increased summer expression of two genes encoding for sucrose-phosphate synthase (Zosma255g00080, Zosma87g00120, Supplementary Table 6). This enzyme is involved in the synthesis of sucrose, the dominant storage carbohydrate in seagrasses (Touchette and Burkholder, 2000). Reduced carbohydrate cleavage in the sampled winter leaves is likely to contribute to the ability of the Norwegian population to survive long dark periods as the reduced demand for metabolic energy prevents depletion of carbon storages.

## Carbon Storage and Stress Under Polar Summer Light

Upregulation of photosynthesis-related genes in the Norwegian population during summer (Figure 6) was also shown in the seagrass *P. oceanica* when exposed to high light levels (Dattolo

et al., 2014) and is likely to enable the plants to process the increased amount of captured light energy (Walters, 2005). Upregulation of genes with functions in carbon capture during summer in the Norwegian (but not the French) samples (Supplementary Figure 3) is necessary for the rapid build-up of stored carbohydrates during the short summer period, and provides a molecular explanation for increasing autumn starch content with latitude (in *Zostera noltii*; Soissons et al., 2018).

Surprisingly, the light intensity at the Norwegian sampling site showed a clear day-night rhythm during summer. Even though the summer sun does not set below the horizon, shading by mountain ridges and the low sun angle during night-time appears to create a short, but clearly noticeable drop in available photosynthetically active radiation. This diurnal light rhythm was paralleled by diurnal changes in gene expression. For example, the daytime upregulation of arginine catabolism (Figure 7) may mobilize nitrogen storages and thus, allow for active growth under nitrogen limitation (Winter et al., 2015). In addition, the nighttime upregulation of trehalose biosynthesis (Figure 7), in which the intermediate compound trehalose-6-phosphate acts as a sensor for available sucrose, is likely involved in the regulation of sucrose-levels that have built up over the day (Lunn et al., 2014; Figueroa and Lunn, 2016) and plays a role in plant growth, development, and protection against various environmental stressors (Lunn et al., 2014; Joshi et al., 2020). Thus, the diurnal changes in gene expression during summer are likely related to active growth and processing of assembled sugars. Local geography, e.g., the presence of mountain ridges, may play an important role in allowing eelgrass to process a dark phase in high latitudes.

Despite diurnal fluctuations in daylight levels, the lack of complete darkness during summer significantly lowers the photosynthetic performance (PI<sub>ABS</sub>) of the Norwegian samples (Figure 2A), in contrast to the seasonal change in photosynthetic performance in the French population (with significantly higher values during summer: Figure 2A). Particularly during summer nights, the Norwegian population showed a significant increase in absorbed (ABS/CS; Figure 2C), trapped (TR<sub>0</sub>/CS; Figure 2D), and electron transport (ET<sub>0</sub>/CS) energy fluxes, as well as the size of the oxidized quinone pool (ABS) and the number of oxidized reaction centers per leaf cross section (RC<sub>0</sub>/CS; Figure 2G). Although these photosynthetic adjustments increase the light harvesting ability and efficiency to use low-light, the Norwegian samples also showed a significant increase in energy dissipation (DI/CS), suggesting an inability to use excessive amounts of light energy in summer.

Low light intensities during the polar summer nights appear to have allowed the Norwegian plants to process a dark phase to recycle intermediates of the Krebs cycle (such as malic acid) that are consumed during the light phase of photosynthesis. The dark phase is triggered by an inhibitor (2-carboxylarabinitol-1,5-bisphosphate, CAP) that is only produced under low light when the chloroplast pH shifts from acid to basic. It was not found in any of the genes showing diurnal differential expression (Supplementary Table 7). A crucial limit to the northward shift of seagrass may be set by the latitude at which the sunlight is too strong to allow the procession of a dark phase. The
existence and identification of this latitude, however, relies on future studies involving transplant experiments to a range of increasing latitudes.

## The Metabolic Dilemma of Carbon Storage Tissue

Our findings suggest that eelgrass can store sufficient energy reserves to survive dark periods lasting longer than those at 67°N (our Norwegian sampling site). In fact, if sugar levels were depleted in mid-winter, we would have expected to see an upregulation of genes involved in protein degradation, as autophagic protein degradation has been suggested as an alternative energy source under sugar depletion in eelgrass (Mochida et al., 2019). As the expected upregulation was not observed, a poleward shift of eelgrass should not be hindered by the potential to store sufficient photosynthesates. In order to test this hypothesis, future studies must evaluate carbohydrate reserves in autumn, a good indicator of winter survival in the seagrass *Zostera noltii* (Govers et al., 2015), and reserve depletion over long dark periods directly in below-ground rhizomes.

While energy storing rhizomes allow eelgrass to survive dark periods, they also represent respiratory burdens that can increase light requirements, potentially preventing poleward migration (Olesen and Sand-Jensen, 1993; Lee et al., 2007). This is particularly true as respiration increases under rising temperatures (Lee et al., 2007; McMinn and Martin, 2013), faster than photosynthesis (higher Q10; Thom et al., 2008; Staehr and Borum, 2011). The respiratory demands of seagrass roots and rhizomes explain the relatively high amounts of light (10–30% surface irradiance; Ochieng et al., 2010) required for long-term survival, as compared with 1–5% for many species of phytoplankton and algae (Kenworthy and Fonseca, 1996). Below-ground rhizomes in eelgrass consume up to 15% of the carbon fixed, and up to 25% at the deepest range-margin (Kraemer and Alberte, 1993). Moreover, even if carbohydrate concentrations are not depleted, seagrass originating from 36°N cannot survive extreme experimental winter light limitation (2h irradiance saturated photosynthesis), presumably through the effects of anoxia on increased translocation of energy reserves (Alcoverro et al., 1999). Thus, increasing carbohydrate storage does not guarantee long-term survival in the dark.

The need of eelgrass to elevate carbon storage in warming Arctic oceans may be lowered by thermal acclimatization that increases respiratory energy efficiency, as has been shown in terrestrial Arctic tundra species (Kornfeld et al., 2013), or by increasing CO₂ levels (Zimmerman et al., 1997). Eelgrass is generally carbon limited, because it cannot efficiently use the inorganic carbon from the water column (HCO₃⁻; Lee et al., 2007). Thus, increasing CO₂ levels are expected to increase photosynthetic rates without affecting respiration (Zimmerman et al., 1997). If turbidity does not hinder long-term survival (Moore et al., 1997) in the future Arctic, increasing CO₂ levels will reduce the need for producing high below-ground biomass by shortening the seasonal period during which eelgrass experiences a negative carbon balance.

## Conclusion and Perspectives

In this study, we show that distinct seasonal expression patterns of an eelgrass population from the Norwegian subarctic are driven by genes encoding functions of the organelles (chloroplast and mitochondrion). These are likely to provide key physiological adaptations that allow eelgrass to endure the polar dark season, and suggest metabolic strategies to enhance restoration success in turbid waters. The two most prominent features were:

1. Active carbon capture in winter. Winter expression of genes encoding for functions related to carbon capture, including Rubisco and PEP-carboxylase (Figure 5), may allow for eelgrass growth and carbon fixation during the polar dark season. This resembles the season-anticipator strategy of polar algae to fix carbon and grow even under the ice during Arctic darkness when nutrient conditions are highest (Berge et al., 2015).

2. Metabolic dormancy in winter. Downregulation of genes related to cellular respiration/carbohydrate cleavage provides the first molecular support for winter-acclimatization in eelgrass that explains an increase in photosynthesis-to-respiration ratios as compared with summer conditions (Thom et al., 2008) and is likely to prevent depletion of carbohydrate storages during the dark period.

This study alone, based on two meadows and an assessment of metabolic and transcriptomic responses, does not allow us to predict whether eelgrass will be able to migrate poleward into a warming, ice-free Arctic. Although this study reveals specificities that are functionally relevant at the species' current northernmost range of distribution, its potential to migrate further poleward depends on the plastic and/or adaptive potential of the identified key functions. While genetic adaptations would rely on standing genetic diversity that allows for evolutionary change over several generations, acclimatization – which can build up within the lifetime of an organism (Kelly et al., 2012; Donelson et al., 2019) – would allow eelgrass to immediately occupy niches that open up in the Arctic. For example, in the seagrass *Thalassia testudinum*, phenotypic variation among subpopulations across an environmental gradient in the Florida Bay were mainly ascribed to phenotypic plasticity because gene flow was high within 50 km distance (Bricker et al., 2011). Further studies are needed to identify to what degree the north-south differences that we identified result from acclimatization or genetic adaptation (Box 1):

1. Reciprocal transplantation of individuals between northern to southern latitudes would allow to identify genetic adaptations as features in photosynthesis and gene expression as such remaining distinct in northern/southern populations and adaptive to northern/southern latitudes even one or more generations after the populations had been transplanted (to minimize non-genetic confounding effects; Kawecki and Ebert, 2004; Pazzaglia et al., 2021).

2. Genetic tests for outlier loci will identify genomic regions that have responded to natural selection in polar regions with exceptionally low allelic variation in 2 high-latitude populations and high genetic divergence from 2 low-latitude
populations (Luikart et al., 2003; Narum and Hess, 2011). Functions associated with genes located in these regions are likely to have adapted genetically to polar conditions. For example, outlier tests identified genes related to growth and immunity as putatively adaptive to varying temperature and light conditions along latitudinal and bathymetric gradients in the seagrass *P. oceanica* (Jahnke et al., 2019).

In order to test whether the reaction norm of eelgrass is wide enough to acclimatize to a warming Arctic, we need to expose eelgrass to a simulation of near-future polar conditions. For example, seagrass meadows in Australia (*Halodule pinifolia* and *Halophila ovalis*) were shown to die at >38 days at 0.1 mol photons m⁻² d⁻¹ (0.35% SI; Longstaff and Dennison, 1999). To determine the latitude at which eelgrass would lack sufficient daylight throughout the year, we would need an equivalent study that imitates future temperature and CO₂ levels at different daylight lengths in a common-garden setup and, ideally, includes a variety of genotypes – since these can differ in their dark tolerance (Salo et al., 2015).

## DATA AVAILABILITY STATEMENT

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at:

https://figshare.com/projects/Adaptation_of_temperate_seagrass_to_Arctic-light_relies_on seasonal acclimatization_of_carbon_capture_and_metabolism/117825

https://www.ncbi.nlm.nih.gov/, PRJNA745558

## REFERENCES

Alcoverro, T., Manzanera, M., and Romero, J. (2001). Annual metabolic carbon balance of the seagrass *Posidonia oceanica*: The importance of carbohydrate reserves. *Mar. Ecol. Prog. Ser.* 211, 105–116. doi: 10.3354/meps211105

Alcoverro, T., Zimmerman, R. C., Kohrs, D. G., and Alberte, R. S. (1999). Resource allocation and sucrose mobilization in light-limited eelgrass *Zostera marina*. *Mar. Ecol. Prog. Ser.* 187, 121–131. doi: 10.3354/meps187121

Alexa, A., and Rahnenführer, J. (2010). Enrichment analysis for Gene Ontology. R Packag. "topGO" v2.16.0. Available at: http://www.bioconductor.org/packages/2.11/bioc/html/topGO.html (Accessed November 02, 2020).

Aubry, S., Brown, N. J., and Hibberd, J. M. (2011). The role of proteins in C 3 plants prior to their recruitment into the C 4 pathway. *J. Exp. Bot.* 62, 3049–3059. doi: 10.1093/jxb/err012

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate - a practical and powerful approach to multiple testing. *J. R. Stat. Soc. B* 57, 289–300. doi: 10.2307/2346101

Bennett, B. D., and Bushel, P. R. (2017). goSTAG: gene ontology subtrees to tag and annotate genes within a set. *Source Code Biol. Med.* 12, 6–13. doi: 10.1186/s13029-017-0066-1

Berge, J., Renaud, P. E., Darnis, G., Cottier, F., Last, K., Gabrielsen, T. M., et al. (2015). In the dark: A review of ecosystem processes during the Arctic polar night. *Prog. Oceanogr.* 139, 258–271. doi: 10.1016/j.pocean.2015.08.005

Blok, S. E., Olesen, B., and Krause-jensen, D. (2018). Life history events of eelgrass *Zostera marina* L. populations across gradients of latitude and temperature. *Mar. Ecol. Prog. Ser.* 590, 79–93. doi: 10.3354/meps12479

https://dataview.ncbi.nlm.nih.gov/object/PRJNA745558?reviewer=d8qlenvsni8uvdohe3mik9n00s.

## AUTHOR CONTRIBUTIONS

GH (project leader) and AJ planned the project and designed the sampling design. AJ, GH, and SA-H collected the shoots. AJ and BD analyzed the data and wrote the manuscript. JC performed the RNA extraction. MK and IS performed the library preparation and sequencing. SA-H, JO, and Z-MH were involved in data interpretation. All authors contributed to the article and approved the submitted version.

## FUNDING

This work was funded by the Norwegian Research Council (Havkyst project 243916). JC was supported by a visiting researcher grant from the Norwegian Research Council (project 281153). Research in the Marine and Environmental Sciences Centre (MARE) was funded by “Fundação para a Ciência e Tecnologia (FCT)” UIDB/04292/2020 and Mar2020 program (VALPRAD project, MAR-01.04.02-FEAMP-0007). BD was supported by FCT throughout an Investigator Contract (CEECIND/03869/2018).

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at https://www.frontiersin.org/articles/10.3389/fpls.2021.745855/full#supplementary-material

Boström, C., Roos, C., and Rönnberg, O. (2004). Shoot morphometry and production dynamics of eelgrass in the northern Baltic Sea. *Aquat. Bot.* 79, 145–161. doi: 10.1016/j.aquabot.2004.02.002

Bricker, E., Waycott, M., Calladine, A., and Zieman, J. C. (2011). High connectivity across environmental gradients and implications for phenotypic plasticity in a marine plant. *Mar. Ecol. Prog. Ser.* 423, 57–67. doi: 10.3354/meps08962

Burke, M. K., Dennison, W. C., and Moore, K. A. (1996). Non-structural carbohydrate reserves of eelgrass *Zostera marina*. *Mar. Ecol. Prog. Ser.* 137, 195–201. doi: 10.3354/meps137195

Burrows, M. T., Schoeman, D. S., Buckley, L. B., Moore, P., Poloczanska, E. S., Brander, K. M., et al. (2011). The pace of shifting climate in marine and terrestrial ecosystems. *Science* 334, 652–655. doi: 10.1126/science.1210288

Bussotti, F., Desotgiu, R., Pollastrini, M., and Cascio, C. (2010). The JIP test: a tool to screen the capacity of plant adaptation to climate change. *Scand. J. For. Res.* 25, 43–50. doi: 10.1080/02827581.2010.485777

Cai, Q., Wang, J., Beletsky, D., Overland, J. E., Ikeda, M., and Wan, L. (2021). Summer Arctic Sea ice decline during 1850-2017 and the amplified Arctic warming during the recent decades. *Environ. Res. Lett.* 16. doi: 10.1088/1748-9326/abdb5f

Chen, A., Li, C., Hu, W., Lau, M. Y., Lin, H., Rockwell, N. C., et al. (2014). PHYTOCHROME C plays a major role in the acceleration of wheat flowering under long-day photoperiod. *Proc. Natl. Acad. Sci. U. S. A.* 111, 10037–10044. doi: 10.1073/pnas.1409795111

Cheung, W. W. L., Lam, V. W. Y., Sarmiento, J. L., Kearney, K., Watson, R., and Pauly, D. (2009). Projecting global marine biodiversity impacts under climate change scenarios. *Fish Fish.* 10, 235–251. doi: 10.1111/j.1467-2979.2008.00315.x
Chi, S., Wu, S., Yu, J., Wang, X., Tang, X., and Liu, T. (2014). Phylogeny of C4-photosynthesis enzymes based on algal transcriptomic and genomic data supports an archaeal/proteobacterial origin and multiple duplication for most C4-related genes. *PLoS One* 9, e110154. doi: 10.1371/journal.pone.0110154

Clarke, K. R., and Gorley, R. N. (2005). PRIMER 6: Getting Started with v6. researchgate.net, 1–12. Available at: https://www.researchgate.net/file. PostFileLoader.html?id=5656d8ce5dbbbdcba28b459f&assetKey=AS:299979836018691@1448532174760 (Accessed April 12, 2021).

Clausen, K. K., and Clausen, P. (2013). Earlier Arctic springs cause phenological mismatch in long-distance migrants. *Oecologia* 173, 1101–1112. doi: 10.1007/s00442-013-2681-0

Clausen, K. K., Krause-Jensen, D., Olesen, B., and Marbà, N. (2014). Seasonality of eelgrass biomass across gradients in temperature and latitude. *Mar. Ecol. Prog. Ser.* 506, 71–85. doi: 10.3354/meps10800

Colarusso, (2007). Natural and stress-induced changes in non-structural carbohydrate concentrations in eelgrass [Zostera marina L.]. Ph. D. thesis, Marine Science Center, Northwestern University, Nahant Massachusetts, USA, 88.

Creswell, L. (2010). Phytoplankton culture for aquaculture feed. *South. Reg. Aquac. Cent.*, 1–55.

Cruz de Carvalho, R., Feijão, E., Kletschkus, E., Marques, J. C., Reis-Santos, P., Fonseca, V. F., et al. (2020). Halophyte bio-optical phenotyping: A multivariate photochemical pressure index (multi-PPi) to classify salt marsh anthropogenic pressures levels. *Ecol. Indic.* 119, 106816. doi: 10.1016/j.ecolind.2020.106816

Dattolo, E., Ruocco, M., Brunet, C., Lorenti, M., Lauritano, C., D'Esposito, D., et al. (2014). Response of the seagrass *Posidonia oceanica* to different light environments: insights from a combined molecular and photo-physiological study. *Mar. Environ. Res.* 101, 225–236. doi: 10.1016/j.marenvres.2014.07.010

Dennison, W., and Alberte, R. (1985). Role of daily light period in the depth distribution of *Zostera marina* (eelgrass). *Mar. Ecol. Prog. Ser.* 25, 51–61. doi: 10.3354/meps025051

Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., et al. (2013). STAR: Ultrafast universal RNA-seq aligner. *Bioinformatics* 29, 15–21. doi: 10.1093/bioinformatics/bts635

Dobin, A., Gingeras, T. R., Spring, C., Flores, R., Sampson, J., Knight, R., et al. (2016). Mapping RNA-seq with STAR. *Curr. Protoc. Bioinforma.* 51, 586–597. doi: 10.1002/0471250953.b1114s51

Donelson, J. M., Sunday, J. M., Figueira, W. F., Gaitán-Espitia, J. D., Hobday, A. J., Johnson, C. R., et al. (2019). Understanding interactions between plasticity, adaptation and range shifts in response to marine environmental change. *Philos. Trans. R. Soc. B Biol. Sci.* 374:20180186. doi: 10.1098/rstb.2018.0186

Duarte, C. M., and Krause-Jensen, D. (2017b). Export from seagrass meadows contributes to marine carbon sequestration. *Front. Mar. Sci.* 4:13. doi: 10.3389/fmars.2017.00013

Duarte, B., Pedro, S., Marques, J. C., Adão, H., and Caçador, I. (2017a). Zostera noltii development probing using chlorophyll a transient analysis (JIP-test) under field conditions: integrating physiological insights into a photochemical stress index. *Ecol. Indic.* 76, 219–229. doi: 10.1016/j.ecolind.2017.01.023

Eriander, L. (2017). Light requirements for successful restoration of eelgrass (*Zostera marina* L.) in a high latitude environment – acclimatization, growth and carbohydrate storage. *J. Exp. Mar. Bio. Ecol.* 496, 37–48. doi: 10.1016/j.jembe.2017.07.010

Figueroa, C. M., and Lunn, J. E. (2016). A tale of two sugars: Trehalose 6-phosphate and sucrose. *Plant Physiol.* 172, 7–27. doi: 10.1104/pp.16.00417

Friedlingstein, P., O'Sullivan, M., Jones, M. W., Andrew, R. M., Hauck, J., Olsen, A., et al. (2020). Global carbon budget 2020. *Earth Syst. Sci. Data* 12, 3269–3340. doi: 10.5194/essd-12-3269-2020

Gattuso, J.-P., Gentili, B., Duarte, C. M., Kleypas, J. A., Middelburg, J. J., and Antoine, D. (2006). Light availability in the coastal ocean: impact on the distribution of benthic photosynthetic organisms and contribution to primary production. *Biogeosci. Discuss.* 3, 895–959. doi: 10.5194/bgd-3-895-2006

Govers, L. L., Suykerbuyk, W., Hoppenreijo, J. H. T., Giesen, K., Bouma, T. J., and Van Katwijk, M. M. (2015). Rhizome starch as indicator for temperate seagrass winter survival. *Ecol. Indic.* 49, 53–60. doi: 10.1016/j.ecolind.2014.10.002

Green, E. P., and Short, F. (2004). World Atlas of Seagrasses. Berkeley, USA: University of California Press.

Griffin, K. L., and Heskel, M. (2013). Breaking the cycle: how light, CO₂ and O₂ affect plant respiration. *Plant Cell Environ.* 36, 498–500. doi: 10.1111/pce.12039

Heskel, M. A., Atkin, O. K., Turnbull, M. H., and Griffin, K. L. (2013). Bringing the Kok effect to light: a review on the integration of daytime respiration

and net ecosystem exchange. *Ecosphere* 4, art98-art14. doi: 10.1890/ ES13-00120.1

Høye, T. T., Post, E., Meltofte, H., Schmidt, N. M., and Forchhammer, M. C. (2007). Rapid advancement of spring in the high Arctic. *Curr. Biol.* 17, R449–R451. doi: 10.1016/j.cub.2007.04.047

Huffeldt, N. P. (2020). Photic barriers to Poleward range-shifts. *Trends Ecol. Evol.* 35, 652–655. doi: 10.1016/j.tree.2020.04.011

Ignatiadis, N., Klaus, B., Zaugg, J. B., and Huber, W. (2016). Data-driven hypothesis weighting increases detection power in genome-scale multiple testing. *Nat. Methods* 13, 577–580. doi: 10.1038/nmeth.3885

Jahnke, M., D'Esposito, D., Orrù, L., Lamontanara, A., Dattolo, E., Badalamenti, F., et al. (2019). Adaptive responses along a depth and a latitudinal gradient in the endemic seagrass *Posidonia oceanica*. *Heredity* 122, 233–243. doi: 10.1038/s41437-018-0103-0

Joshi, R., Sahoo, K. K., Singh, A. K., Anwar, K., Pundir, P., Gautam, R. K., et al. (2020). Enhancing trehalose biosynthesis improves yield potential in marker-free transgenic rice under drought, saline, and sodic conditions. *J. Exp. Bot.* 71, 653–668. doi: 10.1093/jxb/erz462

Kallio, P., and Valanne, N. (1975). "On the Effect of Continuous Light on Photosynthesis in Mosses," in (Springer, Berlin, Heidelberg), 149–162.

Kassambara, A. (2020). rstatix: Pipe-Friendly Framework for Basic Statistical Tests. Available at: https://cran.r-project.org/package=rstatix (Accessed September 18, 2020).

Kawlecki, T. J., and Ebert, D. (2004). Conceptual issues in local adaptation. *Ecol. Lett.* 7, 1225–1241. doi: 10.1111/j.1461-0248.2004.00684.x

Kelly, S. A., Panhuis, T. M., and Stoehr, A. M. (2012). Phenotypic plasticity: molecular mechanisms and adaptive significance. *Compr. Physiol.* 2, 1417–1439. doi: 10.1002/cphy.c11008

Kenworthy, W. J., and Fonseca, M. S. (1996). Light requirements of seagrasses Halodule wrightii and Syringodium filiforme derived from the relationship between diffuse light attenuation and maximum depth distribution. *Estuaries* 19, 740–750. doi: 10.2307/1352533

Koch, M., Bowes, G., Ross, C., and Zhang, X. H. (2013). Climate change and ocean acidification effects on seagrasses and marine macroalgae. *Glob. Chang. Biol.* 19, 103–132. doi: 10.1111/j.1365-2486.2012.02791.x

Kornfeld, A., Heskel, M., Atkin, O. K., Gough, L., Griffin, K. L., Horton, T. W., et al. (2013). Respiratory flexibility and efficiency are affected by simulated global change in Arctic plants. *New Phytol.* 197, 1161–1172. doi: 10.1111/ nph.12083

Kraemer, G. P., and Alberte, R. S. (1993). Age-related patterns of metabolism and biomass in subterranean tissues of *Zostera marina* (eelgrass). *Mar. Ecol. Prog. Ser.* 95, 193–203. doi: 10.3354/meps095193

Krause-Jensen, D., and Duarte, C. M. (2014). Expansion of vegetated coastal ecosystems in the future Arctic. *Front. Mar. Sci.* 1:77. doi: 10.3389/ fmars.2014.00077

Larsson, J. (2019). eulerr: Area-Proportional Euler and Venn Diagrams with Ellipses. R package version 6.0.0. Available at: https://cran.r-project.org/package=eulerr (Accessed November 10, 2020).

Lè, S., Josse, J., and Husson, F. (2008). FactoMineR: an R package for multivariate analysis. *J. Stat. Softw.* 35, 253–258. doi: 10.1016/j.envint.2008.06.007

Lee, K. S., Park, S. R., and Kim, Y. K. (2007). Effects of irradiance, temperature, and nutrients on growth dynamics of seagrasses: a review. *J. Exp. Mar. Bio. Ecol.* 350, 144–175. doi: 10.1016/j.jembe.2007.06.016

Lepiniec, L., Vidal, J., Chollet, R., Gadal, P., and Crétin, C. (1994). Phosphoenolpyruvate carboxylase: structure, regulation and evolution. *Plant Sci.* 99, 111–124. doi: 10.1016/0168-9452(94)90168-6

Liao, Y., Smyth, G. K., and Shi, W. (2019). The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads. *Nucleic Acids Res.* 47:e47. doi: 10.1093/nar/gkz114

Longstaff, B. J., and Dennison, W. C. (1999). Seagrass survival during pulsed turbidity events: The effects of light deprivation on the seagrasses *Halodule pinifolia* and *Halophila ovalis*. *Aquat. Bot.* 65, 105–121. doi: 10.1016/ S0304-3770(99)00035-2

Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biol.* 15, 550. doi: 10.1186/s13059-014-0550-8

Luikart, G., England, P. R., Tallmon, D., Jordan, S., and Taberlet, P. (2003). The power and promise of population genomics: from genotyping to genome typing. *Nat. Rev. Genet.* 4, 981–994. doi: 10.1038/nrg1226
Lunn, J. E., Delorge, I., Figueroa, C. M., Van Dijck, P., and Stitt, M. (2014). Trehalose metabolism in plants. *Plant J.* 79, 544–567. doi: 10.1111/tpj.12509

Luo, W., and Brouwer, C. (2013). Pathview: An R/bioconductor package for pathway-based data integration and visualization. *Bioinformatics* 29, 1830–1831. doi: 10.1093/bioinformatics/btt285

Luo, W., Friedman, M. S., Shedden, K., Hankenson, K. D., and Woolf, P. J. (2009). GAGE: generally applicable gene set enrichment for pathway analysis. *BMC Bioinformatics* 10:161. doi: 10.1186/1471-2105-10-161

Mazzuca, S., Spadafora, A., Filadoro, D., Vannini, C., Marsoni, M., Cozza, R., et al. (2009). Seagrass light acclimation: 2-DE protein analysis in Posidonia leaves grown in chronic low light conditions. *J. Exp. Mar. Bio. Ecol.* 374, 113–122. doi: 10.1016/j.jembe.2009.04.010

McMahon, C. R., and Hays, G. C. (2006). Thermal niche, large-scale movements and implications of climate change for a critically endangered marine vertebrate. *Glob. Chang. Biol.* 12, 1330–1338. doi: 10.1111/j.1365-2486.2006.01174.x

McMinn, A., and Martin, A. (2013). Dark survival in a warming world. *Proc. R. Soc. Biol. Sci.* 280:20122909. doi: 10.1098/rspb.2012.2909

Melbourne-Thomas, J., Audzijonyte, A., Brasier, M. J., Cresswell, K., Fogarty, H. E., Haward, M., et al. (2020). Poleward bound: adapting to climate-driven species redistribution. *Rev. Fish Biol. Fish.* doi: 10.1007/s11160-021-09641-3 [Epub ahead of print]

Messiaen, G., Mortreux, S., Le Gall, P., Crottier, A., and Lagarde, F. (2018). Marine environmental station database of Thau lagoon. *Seanoe.* doi: 10.17882/52404

Mochida, K., Hano, T., Onduka, T., Ito, K., and Yoshida, G. (2019). Physiological responses of eelgrass (*Zostera marina*) to ambient stresses such as herbicide, insufficient light, and high water temperature. *Aquat. Toxicol.* 208, 20–28. doi: 10.1016/j.aquatox.2018.12.018

Moore, K. A., Orth, R. J., and Nowak, J. F. (1993). Environmental regulation of seed germination in *Zostera marina* L. (eelgrass) in Chesapeake Bay: effects of light, oxygen and sediment burial. *Aquat. Bot.* 45, 79–91. doi: 10.1016/0304-3770(93)90054-Z

Moore, K. A., Wetzel, R. L., and Orth, R. J. (1997). Seasonal pulses of turbidity and their relations to eelgrass (*Zostera marina* L.) survival in an estuary. *J. Exp. Mar. Bio. Ecol.* 215, 115–134. doi: 10.1016/S0022-0981(96)02774-8

Murtagh, F., and Legendre, P. (2014). Ward's hierarchical agglomerative clustering method: which algorithms implement Ward's criterion? *J. Classif.* 31, 274–295. doi: 10.1007/s00357-014-9161-z

Narum, S. R., and Hess, J. E. (2011). Comparison of FST outlier tests for SNP loci under selection. *Mol. Ecol. Resour.* 11, 184–194. doi: 10.1111/j.1755-0998.2011.02987.x

Novak, A. B., Pelletier, M. C., Colarusso, P., Simpson, J., Gutierrez, M. N., Arias-Ortiz, A., et al. (2020). Correction to: factors influencing carbon stocks and accumulation rates in eelgrass meadows Across New England, USA. *Estuar. Coasts* 43, 2183–2184. doi: 10.1007/s12237-020-00815-z (Estuaries and Coasts, (2020), 43, 8, (2076–2091), 10.1007/s12237-020-00754-9)

Ochieng, C. A. (2008). Survival Strategies of Eelgrass in Reduced Light. Vasa. Available at: http://medcontent.metapress.com/index/A65RM03P4874243N.pdf (Accessed March 08, 2021).

Ochieng, C. A., Short, F. T., and Walker, D. I. (2010). Photosynthetic and morphological responses of eelgrass (*Zostera marina* L.) to a gradient of light conditions. *J. Exp. Mar. Bio. Ecol.* 382, 117–124. doi: 10.1016/j.jembe.2009.11.007

Olesen, B., Krause-Jensen, D., Marbà, N., and Christensen, P. B. P. (2015). Eelgrass *Zostera marina* in subarctic Greenland: dense meadows with slow biomass turnover in cold waters. *Mar. Ecol. Prog. Ser.* 518, 107–121. doi: 10.3354/meps11087

Olesen, B., and Sand-Jensen, K. (1993). Seasonal acclimatization of eelgrass *Zostera marina* growth to light. *Mar. Ecol. Prog. Ser.* 94, 91–99. doi: 10.3354/meps094091

Olsen, J. L., Coyer, J. A., Stam, W. T., Moy, F. E., Christie, H., and Jorgensen, N. M. (2013). Eelgrass *Zostera marina* populations in northern Norwegian fjords are genetically isolated and diverse. *Mar. Ecol. Prog. Ser.* 486, 121–132. doi: 10.3354/meps10373

Olsen, J. L., Rouzé, P., Verhelst, B., Lin, Y., Bayer, T., Collen, J., et al. (2016). The genome of the seagrass *Zostera marina* reveals angiosperm adaptation to the sea. *Nature* 530, 331–335. doi: 10.1038/nature16548

Ort, D. R., Zhu, X., and Melis, A. (2011). Optimizing antenna size to maximize photosynthetic efficiency. *Plant Physiol.* 155, 79–85. doi: 10.1104/pp.110.165886

Pazzaglia, J., Reusch, T. B. H., Terlizzi, A., Marín-Guirao, L., and Procaccini, G. (2021). Phenotypic plasticity under rapid global changes: The intrinsic force for future seagrasses survival. *Evol. Appl.* 14, 1181–1201. doi: 10.1111/eva.13212

Pérez-Silva, J. G., Araujo-Voces, M., and Quesada, V. (2018). NVenn: Generalized, quasi-proportional Venn and Euler diagrams. *Bioinformatics* 34, 2322–2324. doi: 10.1093/bioinformatics/bty109

Poloczanska, E. S., Brown, C. J., Sydeman, W. J., Kiessling, W., Schoeman, D. S., Moore, P. J., et al. (2013). Global imprint of climate change on marine life. *Nat. Clim. Chang.* 3, 919–925. doi: 10.1038/nclimate1958

R Core Team (2019). R: A Language and Environment for Statistical Computing. Available at: https://www.r-project.org/ (Accessed Febraury 05, 2021).

Ralph, P. J., Durako, M. J., Enriquez, S., Collier, C. J., and Doblin, M. A. (2007). Impact of light limitation on seagrasses. *J. Exp. Mar. Bio. Ecol.* 350, 176–193. doi: 10.1016/j.jembe.2007.06.017

Raven, J. A., and Farquhar, G. D. (1990). The influence of N metabolism and organic acid synthesis on the natural abundance of isotopes in carbon plants. *New Phytol.* 116, 505–529. doi: 10.1111/j.1469-8137.1990.tb00536.x

Ruocco, M., De Luca, P., Marín-Guirao, L., and Procaccini, G. (2019a). Differential leaf age-dependent thermal plasticity in the keystone Seagrass *Posidonia oceanica*. *Front. Plant Sci.* 10. doi: 10.3389/fpls.2019.01556

Ruocco, M., Marín-Guirao, L., and Procaccini, G. (2019b). Within- and among-leaf variations in photo-physiological functions, gene expression and DNA methylation patterns in the large-sized seagrass *Posidonia oceanica*. *Mar. Biol.* 166, 1–18. doi: 10.1007/s00227-019-3482-8

Salesse-Smith, C. E., Sharwood, R. E., Busch, F. A., and Stern, D. B. (2020). Increased Rubisco content in maize mitigates chilling stress and speeds recovery. *Plant Biotechnol. J.* 18, 1409–1420. doi: 10.1111/pbi.13306

Salo, T., Reusch, T. B. H., and Boström, C. (2015). Genotype-specific responses to light stress in eelgrass *Zostera marina*, a marine foundation plant. *Mar. Ecol. Prog. Ser.* 519, 129–140. doi: 10.3354/meps11083

Sayols, S., Scherzinger, D., and Klein, H. (2016). dupRadar: A bioconductor package for the assessment of PCR artifacts in RNA-Seq data. *BMC Bioinformatics* 17:428. doi: 10.1186/s12859-016-1276-2

Soissons, L. M., Haanstra, E. P., Van Katwijk, M. M., Asmus, R., Auby, I., Barillé, L., et al. (2018). Latitudinal patterns in european seagrass carbon reserves: influence of seasonal fluctuations versus short-term stress and disturbance events. *Front. Plant Sci.* 9. doi: 10.3389/fpls.2018.00088

Staehr, P. A., and Borum, J. (2011). Seasonal acclimation in metabolism reduces light requirements of eelgrass (*Zostera marina*). *J. Exp. Mar. Bio. Ecol.* 407, 139–146. doi: 10.1016/j.jembe.2011.05.031

Sterck, L., Billiau, K., Abeel, T., Rouzé, P., and Van de Peer, Y. (2012). ORCAE: online resource for community annotation of eukaryotes. *Nat. Methods* 9:1041. doi: 10.1038/nmeth.2242

Sysoeva, M., Markovskaya, E., and Shibaeva, T. (2010). Plants under continuous light: A review. *Plant Stress* 4, 5–17.

Taulavuori, K., Sarala, M., and Taulavuori, E. (2010). "Growth responses of trees to arctic light environment" in *Progress in Botany*. eds. U. Lüttge, W. Beyschlag, B. Büdel and D. Francis (Berlin, Heidelberg: Springer), 157–168.

Thom, R. M., Southard, S. L., Borde, A. B., and Stoltz, P. (2008). Light requirements for growth and survival of eelgrass (*Zostera marina* L.) in Pacific northwest (USA) estuaries. *Estuar. Coasts* 31, 969–980. doi: 10.1007/s12237-008-9082-3

Touchette, B. W., and Burkholder, J. A. M. (2000). Overview of the physiological ecology of carbon metabolism in seagrasses. *J. Exp. Mar. Bio. Ecol.* 250, 169–205. doi: 10.1016/S0022-0981(00)00196-9

Trombetta, T., Vidussi, F., Mas, S., Parin, D., Simier, M., and Mostajir, B. (2019). Water temperature drives phytoplankton blooms in coastal waters. *PLoS One* 14, e0214933. doi: 10.1371/journal.pone.0214933

Walters, R. G. (2005). Towards an understanding of photosynthetic acclimation. *J. Exp. Bot.* 56, 435–447. doi: 10.1093/jxb/eri060

Ward, J. H. (1963). Hierarchical grouping to optimize an objective function. *J. Am. Stat. Assoc.* 58, 236–244. doi: 10.1080/01621459.1963.10500845

Warnes, G. R., Bolker, B., Bonebakker, L., Gentleman, R., Huber, W., Liaw, A., et al. (2020). gplots: Various R Programming Tools for Plotting Data. Available at: https://cran.r-project.org/package=gplots (Accessed June 16, 2021)..
Wernberg, T., Smale, D., and Tuya, F. (2013). An extreme climatic event alters marine ecosystem structure in a global biodiversity hotspot. *Nat. Clim.* 3, 78–82. doi: 10.1038/nclimate1627

Wiencke, C., Clayton, M. N., Gómez, I., Iken, K., Lüder, U. H., Amsler, C. D., et al. (2007). Life strategy, ecophysiology and ecology of seaweeds in polar waters. *Life Extrem. Environ.* 6, 213–244. doi: 10.1007/978-1-4020-6285-8-13

Wiencke, C., Gómez, I., and Dunton, K. (2009). Phenology and seasonal physiological performance of polar seaweeds. *Bot. Mar.* 52, 585–592. doi: 10.1515/BOT.2009.078

Wilson, K. L., and Lotze, H. K. (2019). Climate change projections reveal range shifts of eelgrass *Zostera marina* in the Northwest Atlantic. *Mar. Ecol. Prog. Ser.* 620, 47–62. doi: 10.3354/meps12973

Winter, G., Todd, C. D., Trovato, M., Forlani, G., and Funck, D. (2015). Physiological implications of arginine metabolism in plants. *Front. Plant Sci.* 6:534. doi: 10.3389/fpls.2015.00534

Wium-Andersen, G., and Borum, J. (1984). Biomass variation and autotrophic production of an epiphyte-macrophyte community in a coastal danish area: I. eelgrass (*Zostera marina* L.) biomass and net production. *Ophelia* 23, 33–46. doi: 10.1080/00785236.1984.10426603

Wong, M. C., Vercaemer, B. M., and Griffiths, G. (2021). Response and recovery of eelgrass (*Zostera marina*) to chronic and episodic light disturbance. *Estuar. Coasts* 44, 312–324. doi: 10.1007/s12237-020-00803-3

Woods, D. P., Ream, T. S., Minevich, G., Hobert, O., and Amasino, R. M. (2014). PHYTOCHROME C is an essential light receptor for photoperiodic flowering in the temperate grass, Brachypodium distachyon. *Genetics* 198, 397–408. doi: 10.1534/genetics.114.166785

Zhang, P. D., Liu, Y. S., Guo, D., Li, W. T., and Zhang, Q. (2016). Seasonal variation in growth, morphology, and reproduction of eelgrass *Zostera marina* on the eastern coast of the Shandong peninsula, China. *J. Coast. Res.* 318, 315–322. doi: 10.2112/JCOASTRES-D-14-00117.1

Zhu, X. G., Govindjee, J., Baker, N. R., DeSturler, E., Ort, D. R., and Long, S. P. (2005). Chlorophyll a fluorescence induction kinetics in leaves predicted from a model describing each discrete step of excitation energy and electron transfer associated with photosystem II. *Planta* 223, 114–133. doi: 10.1007/ s00425-005-0064-4

Zimmerman, R. C., Kohrs, D. G., Steller, D. L., and Alberte, R. S. (1997). Impacts of CO2 enrichment on productivity and light requirements of eelgrass. *Plant Physiol.* 115, 599–607. doi: 10.1104/pp.115.2.599

Zimmerman, R. C., Smith, R. D., and Alberte, R. S. (1989). Thermal acclimation and whole-plant carbon balance in *Zostera marina* L. (eelgrass). *J. Exp. Mar. Bio. Ecol.* 130, 93–109. doi: 10.1016/0022-0981(89)90197-4

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2021 Jueterbock, Duarte, Coyer, Olsen, Kopp, Smolina, Arnaud-Haond, Hu and Hoarau. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.