- 1Stazione Zoologica Anton Dohrn, Naples, Italy
- 2Centre of Marine Sciences, University of Algarve, Faro, Portugal
- 3GreenCoLab - Associação Oceano Verde, University of Algarve, Faro, Portugal
The circadian clock is an endogenous time-keeping mechanism that enables organisms to adapt to external environmental cycles. It produces rhythms of plant metabolism and physiology, and interacts with signaling pathways controlling daily and seasonal environmental responses through gene expression regulation. Downstream metabolic outputs, such as photosynthesis and sugar metabolism, besides being affected by the clock, can also contribute to the circadian timing itself. In marine plants, studies of circadian rhythms are still way behind in respect to terrestrial species, which strongly limits the understanding of how they coordinate their physiology and energetic metabolism with environmental signals at sea. Here, we provided a first description of daily timing of key core clock components and clock output pathways in two seagrass species, Cymodocea nodosa and Zostera marina (order Alismatales), co-occurring at the same geographic location, thus exposed to identical natural variations in photoperiod. Large differences were observed between species in the daily timing of accumulation of transcripts related to key metabolic pathways, such as photosynthesis and sucrose synthesis/transport, highlighting the importance of intrinsic biological, and likely ecological attributes of the species in determining the periodicity of functions. The two species exhibited a differential sensitivity to light-to-dark and dark-to-light transition times and could adopt different growth timing based on a differential strategy of resource allocation and mobilization throughout the day, possibly coordinated by the circadian clock. This behavior could potentially derive from divergent evolutionary adaptations of the species to their bio-geographical range of distributions.
Introduction
Earth’s rotation causes repetitive changes of day and night that are reflected in diurnal cycles of light and temperature. Most living organisms experience these changes and thus face the challenge of coordinating their lives with such rhythms. This was the primary driver for the emergence and evolution of endogenous clocks that can be set by the rising or the setting of the sun (McWatters and Devlin, 2011). In plants, circadian clocks regulate many aspects of physiology, growth, development and reproduction (Harmer et al., 2000; Harmer, 2009). These produce daily rhythms of metabolism, and interact with signaling pathways controlling daily environmental responses (Dodd et al., 2015; Kim et al., 2017). Ultimately, circadian clocks allow plants to anticipate daily and seasonal changes in the environment, conferring them an adaptive advantage (Ouyang et al., 1998).
Over the last decade, a combination of empirical research and mathematical modeling has led to the identification of several clock components in terrestrial model species and an understanding of their regulative roles within plants’ biochemical pathways (Webb and Satake, 2015). Three interconnected modules mainly compose the plants’ circadian-clock system: the receptors that perceive environmental signals and provide the inputs, the core oscillator that generates the rhythm itself, and the outputs that determine physiological and metabolic rhythms (McWatters and Devlin, 2011). In Arabidopsis thaliana, the clockwork is regulated by at least three tightly interlocked feedback loops, where clock proteins regulate their own transcription, either directly or indirectly (Barak et al., 2000; Harmer, 2009). These loops of gene expression patterns are synchronized by post-translational modifications of their own proteins.
Different photic cues (e.g., day-length, intensity, spectral composition of sun light and solar position) may contribute to the entrainment of circadian clock of plants (Millar, 2004). These signals are perceived via three regulatory families of photoreceptors including phytochromes (PHY), which are sensitive to red/far-red, and cryptochromes (CRY) and phototropins (PHOT), which adsorb in the UV-blue light (Millar, 2003), although phototropins seems not to have a direct role in light input into the circadian clock (Jones, 2009; Litthauer et al., 2016). The ZEITLUPE (ZTL) family that belongs to the LOV domain photoreceptors is the third family of blue light photoreceptors in plants, with a central role in circadian clock and photoperiodic flowering (Song et al., 2014).
Genome-wide studies in terrestrial plants have revealed that a large portion (up to 89%) of the expressed transcripts show a rhythmic expression under diverse circadian, photocycle and thermocycle growth conditions (e.g., Harmer et al., 2000; Michael et al., 2008; Zdepski et al., 2008; Filichkin et al., 2011; Jończyk et al., 2011). This is especially true for transcripts involved in main physiological pathways, such as sugar synthesis/utilization, which imply a major role of the clock system in coordinating key aspects of the plant metabolism (Harmer et al., 2000; Covington et al., 2008; Haydon et al., 2013a). Several sugar-responsive transcripts show rhythmic changes that match the diurnal variation in sugar levels (Bläsing et al., 2005). On the other hand, sugar levels themselves can entrain circadian rhythms by regulating the expression of clock-components (Haydon et al., 2013b). In general, the circadian oscillator is coupled with carbohydrate biochemistry and there is a tight connection between nuclear-encoded circadian-clock components and processes occurring in the chloroplasts (Dodd et al., 2015). For instance at dawn the phase of the circadian oscillator responds to the low light intensity detected by photoreceptors. Then, a second so-called “metabolic dawn” occurs in response to the first accumulation of sugars due to the sensitivity of clock components such as LHY/CCA1 (Haydon et al., 2013b; Dodd et al., 2015). The daily dynamic of starch production and consumption is also known to be regulated by the circadian oscillator with feedbacks from carbon availability (Haydon et al., 2013b; Dodd et al., 2015; Seki et al., 2017; Flis et al., 2019). The molecular mechanisms that coordinate the turnover of starch with the external photoperiod are currently unknown, although two different models have been proposed, mainly differing in the way the temporal and metabolic cues are integrated. In the first, it is assumed that the abundance of starch is measured (Scialdone et al., 2013), while in the other, it is assumed that sucrose is measured (Feugier and Satake, 2013; Webb and Satake, 2015; Seki et al., 2017). In the starch sensing models, the circadian clock is a passive timer used to measure the time of day (Scialdone et al., 2013). In the second, sucrose feedbacks to the circadian clock to dynamically regulate the phase of the circadian oscillator (Feugier and Satake, 2013; Seki et al., 2017).
The structure and function of the circadian system might be widely conserved across higher plant species (Serikawa et al., 2008; Filichkin et al., 2011); however, the regulation of clock-output traits can be species-specific such that predictions from model species as Arabidopsis are not always possible. For example, recent studies have highlighted fundamental differences in the circadian behavior of traits like growth rate in dicotyledonous vs. monocotyledons species (Müller et al., 2014). Such differences question if the clock oscillator has the same importance to regulate physiology, metabolism and growth across species (Müller et al., 2014). In addition, several studies have examined the potential ecological implications of endogenous differences in plant photo-perception and timing systems (Resco et al., 2009) on questions such as biogeographical distribution of species and range-shifts under climate changes (Huffeldt, 2020), and on the relevance that circadian regulation of physiology has on the interactions at ecosystem level (Resco de Dios et al., 2016; Lu et al., 2017; Hubbard et al., 2018) and biosphere–atmosphere regulations (Resco de Dios, 2013).
If the transfer of knowledge regarding circadian-regulated patterns is already complex across species inhabiting similar habitats, this is even more difficult in species experiencing completely different habitats, such as seagrasses that live submerged in marine environments. Seagrasses are a polyphyletic assemblage of ∼70 species, which evolved from land ancestors and returned to the sea in the Cretaceous period (Les et al., 1997). Seagrasses are distributed along temperate and tropical coastlines worldwide (Short et al., 2007) where they form marine forests that are among the most productive and valuable costal ecosystems (Ruiz-Frau et al., 2017) and contribute to climate change mitigation (Marbà et al., 2015). The current distribution of seagrass species across latitudes reflects their specific requirements in terms of light and temperature for growth (Short et al., 2007) and phenology (Blok et al., 2018; Yue et al., 2020), which could be influenced by circadian and/or photoperiodic responses.
Some of the biggest challenges imposed by the marine environment in respect to the terrestrial one, are the differences in light intensity and quality (Kirk, 2010), as well as the lower CO2 availability (Larkum et al., 2006, 2018). The exposure to such a different habitat could have influenced seagrasses’ circadian clock, in terms of composition of gene networks and their subsequent regulation (Olsen et al., 2016). A first evidence of rhythmic diel oscillations of transcripts related to photosynthesis and respiration in seagrasses has come from preliminary reports in Posidonia oceanica and Zostera marina (Procaccini et al., 2017; Rasmusson et al., 2017). However, studies of circadian rhythms in marine plants are still way behind in comparison to terrestrial species, notably for the lack of genomes available for most species and the inapplicability of fundamental molecular genetics approaches. This strongly limits the possibility to investigate the interconnection between energy metabolism and clock systems in marine plants, as well as how circadian and photoperiodic responses may influence their geographical distribution, productivity and survival under climate change.
Here, we want to provide a comprehensive overview of daily transcriptional timing of putative core circadian clock components and metabolic pathways such as sugar production/utilization, photoreception, photosynthesis, and cellular respiration, as well as patterns of carbohydrate accumulation, in two seagrass species, Cymodocea nodosa and Zostera marina (order Alismatales). These species are characterized by distinct biogeographical distributional ranges: C. nodosa typically spreads in warmer and temperate waters (from Mauritania to the Canary archipelago and Mediterranean Sea), while Z. marina is the dominant seagrass species in colder water of North hemisphere (Short et al., 2007). We explored daily responses of the species when co-occurring at the same geographic location (Ria Formosa lagoon, Faro, Portugal), thus being submitted to identical light and temperature cycles. This area represents the northern and southern distribution edge in the Atlantic region of C. nodosa and Z. marina, respectively. Our hypothesis is that possible inter-specific variations in daily metabolic rhythms, could reflect e.g., evolutionary adaptations of the species to divergent latitudinal distribution related to their optimal photic and thermal requirements.
Our study was conducted under natural light conditions, thus it aims at uncovering the rhythms of key pathways in the two species under naturally fluctuating conditions.
Materials and Methods
Seagrass Sampling and Experimental Design
Seagrass specimens used for this study are the same as in Ruocco et al. (2020a). In May 2019, entire shoots of Cymodocea nodosa and Zostera marina were collected by snorkeling from a shallow-water (1-2 m depth) mixed meadow near Ilha da Culatra (Faro, Portugal – 36°59′41″ N, 7°50′26″ W) with a significant portion of the sediment. Plants and sediment were then transported to the nearby Ramalhete field station (CCMAR, Centre of Marine Sciences – 37°00′22″ N, 7°58′02″ W) in darkened containers filled with seawater, limiting rhizome and roots breakage. Following removal from the field, shoots were transplanted to six 150-L cylindrical tanks located in an outdoor mesocosm facility exposed to a natural daily light regime. The bottom of each tank was covered with ca. 7 cm washed beach sand. Tanks were supplied with running seawater from the nearby area, previously filtered with sand and UV filters, in an independent open circuit configuration. Seagrass shoots (ca. 20) were randomly assigned to each tank. Plants were left to recover for 24 h following transplantation. At the end of the acclimation period, leaf sections of C. nodosa and Z. marina were collected 6 times during a 24-h cycle, i.e., at midnight (00:00), dawn (04:30), sunrise (06:30), solar noon (13:00), sunset (20:30), and dusk (22:00). All information about sunrise and sunset times were taken from https://www.timeanddate.com/. Irradiance (as Photosynthetically Active Radiation, PAR) and water temperature were monitored continuously in the tanks with a Li-192SA underwater quantum sensor connected to a Li-1400 data logger (Li-Cor), and a HOBO temperature logger, respectively. Daily patterns of PAR and temperature levels are displayed in Supplementary Figure 1. For all analyses, three independent biological replicates were collected for each species and time point (n = 3) and used for both gene expression and carbohydrate analyses. Leaf samples were rapidly cleaned of epiphytes, rinsed with distilled water and immediately submerged in RNAlater© tissue collection (Ambion, Life Technologies, Waltham, MA, United States). After one night at 4°C, leaf samples were definitely stored at –20°C.
RNA Extraction and Target Genes Selection
Total RNA from the youngest fully developed leaves of C. nodosa and Z. marina (Ruocco et al., 2019) was extracted with the Aurum™ Total RNA Mini Kit (BIO-RAD, Hercules, CA, United States), following the manufacturer’s protocol. About 5–7 cm-long leaf sections corresponding to about 70–100 mg were ground to a fine powder with a mortar and pestle containing liquid nitrogen. Samples were then homogenized through a Mixer Mill MM300 (QIAGEN, Hilden, Germany) and tungsten carbide beads (3 mm) for 3 min at 20.1 Hz. The quality and purity of the total RNA was checked using NanoDrop (ND-1000 UV-Vis spectrophotometer; NanoDrop Technologies, Wilmington, DE, United States) and 1% agarose gel electrophoresis. RNA was used when Abs260 nm/Abs280 nm and Abs260 nm/Abs230 nm ratios were > 1.8 and 1.8 < x < 2, respectively. The RNA concentration was accurately determined by the Qubit™ RNA BR Assay kit (Thermo Fisher Scientific, Waltham, MA, United States) using the Qubit 2.0 Fluorometer (Thermo Fisher Scientific). The total RNA (500 ng) from each sample was reverse-transcribed into cDNA with the iScript™ cDNA synthesis kit (BIO-RAD), according to the manufacturer’s protocol.
Primers for target genes (19) involved in key plant metabolic pathways potentially exhibiting circadian changes i.e., sugar metabolism, photosynthesis and light harvesting, photoreception and circadian clock, growth (Harmer et al., 2000), were newly designed with the primer analysis software Primer3 v. 0.4.0 (Untergasser et al., 2012) or selected from previous studies (see Table 1). The design conditions included the primer length (18–23 bp), Tm (∼60 C), GC content (50%), and product size (100–250 bp). To identify target sequences for Genes of Interests (GOIs), the corresponding protein sequences of A. thaliana were downloaded from TAIR1 and used as queries for BLASTP (Altschul et al., 1990). For Z. marina, sequences were blasted against the genome available in ORCAE2 (Olsen et al., 2016). For C. nodosa, potential homologous sequences were identified through sequence similarity search against an available transcriptome (Ruocco et al., 2017) or public databases. The BLASTP analysis was carried out setting a threshold of homology to 1e-10. The top hit for each query was selected and the homology level of each pair of target sequences of Z. marina vs. C. nodosa, was reported (Supplementary Tables 1a,b). To verify the evolutionary relationships among target genes identified in seagrasses in comparison to land plants, homologous sequences of GOIs of other three plant species (Oryza sativa, Zea mays, and Populus trichocarpa), representatives of monocots and eudicots lineages, were downloaded from public databases (i.e., KEGG, Phytozome v13, and UniProt) and included in our analysis (Supplementary Table 2). Multiple sequence alignments were performed using COBALT (Papadopoulos and Agarwala, 2007) with default parameters. A further phylogenetic analysis was carried out on GOIs for which potential orthologues sequences were identified between Z. marina and C. nodosa using the Phylemon2 toolkit (Sánchez et al., 2011; Supplementary Table 3). Poorly aligned sequence regions were removed with trimAl v1.3 (Capella-Gutiérrez et al., 2009) using the strictplus option. For each dataset, a maximum likelihood phylogenetic tree was inferred using PhyML (Guindon et al., 2010) implemented in the online tool Phyml Best AIC Tree v.1.02b. Topology optimization of best model was calculated for each independent dataset and the support of nodes was calculated using the SH-like procedure (Anisimova and Gascuel, 2006). Generated trees were visualized in UGENE v.33.0 (Okonechnikov et al., 2012) and graphically edited in Inkscape (Supplementary Figure 2).
Daily Gene-Expression Analysis
Daily patterns of transcript abundance were obtained via RT-qPCR. All reactions were carried out as outlined in Ruocco et al. (2019). Briefly, each reaction consisted in 5 μl Fast SYBR® Green Master Mix (Applied Biosystems), 1 μl cDNA (1:5 diluted) template, and 4 μl of 0.7 pmol μl–1 primers. The thermal profile of the reactions was as follows: 95°C for 20 s, 40 times 95°C for 1 s, and 60°C for 20 s. All RT-qPCR reactions were conducted in triplicate, and each assay included three no-template negative controls. To normalize target gene expression, we used the eukaryotic initiation factor 4A (eIF4A) as a reference gene (RG), which was previously demonstrated to exhibit a stable expression in Z. marina and C. nodosa along a daily cycle (Rasmusson et al., 2017; Ruocco et al., 2020a) and under a range of different conditions (Ransbotyn and Reusch, 2006; Winters et al., 2011; Olivé et al., 2017). RT-qPCR efficiencies for all primer pairs were calculated from the slopes of the standard curves of the threshold cycle (CT) vs. cDNA concentration with the equation E = 10–1/slope. The relative quantification of the transcript levels was obtained following (Schmittgen and Livak, 2008). In details, the negative differences in the cycles to cross the threshold value between the RG and the respective GOI (-ΔCT) were calculated according to the equation:
Mean –ΔCT values were then calculated for biological replicates at each time point (n = 3) from individual –ΔCT values. Normalized values of diurnal transcript levels of GOIs have been graphically visualized using SigmaPlot v.12.5.
Total Non-structural Carbohydrate Determination
To assess the pattern of sugar accumulation/consumption throughout the day, the total content of non-structural carbohydrates (TNC; soluble sugars and starch) was analyzed in seagrass leaf tissues (n = 3) at five sampling times [i.e., at sunrise (06:30), solar noon (13:00), sunset (20:30), midnight (0:00), dawn (04:30)]. For TNC analysis, leaves (ca. 50 mg) were finely ground in liquid nitrogen, and subsequent extraction was carried out following Ruocco et al. (2020b). Briefly, soluble sugars were solubilized by three sequential extractions with 80% (v/v) ethanol at 80°C for 15 min. After centrifugation (3,000 rpm for 10 min), the ethanol extract was used for the determination of soluble sugars’ content, while the pellet was hydrolysed for starch determination (24 h at RT) with 3 ml NaOH 0.1 M. For both soluble sugars and starch, 3% aqueous phenol (0.25 ml) and 95%–97% H2SO4 (2.5 ml) were mixed with 1 ml sample in glass tubes and the solution was allowed to rest for 30 min. Absorbance was then read at 490 and 750 nm with an Agilent 8453 UV–Vis spectrophotometer. TNC content was calculated using sucrose calibration curves (standard sucrose 99%, Biorad).
Statistical Analysis
A multivariate analysis was first used to assess the overall signal of all 19 GOIs (based on –ΔCT values). In details, a two-way permutational multivariate analysis of variance (PERMANOVA) based on the Bray-Curtis matrix, was conducted with the Primer 6 v.6.1.12 and PERMANOVA + v.1.0.2 software package. The analysis consisted of two fixed factors: Species (Sp) with two levels (i.e., C. nodosa and Z. marina) and Time (Ti) with six levels (i.e., midnight, dawn, sunrise, solar noon, sunset, and dusk). A principal component analysis (PCA) based on the Bray-Curtis matrix and clustering analysis were also performed on the multivariate gene-expression datasets with the software PAST v.3.03 (Hammer et al., 2001). Subsequently, univariate analyses (one-way ANOVAs and two-way ANOVAs), were used to assess the effects of time and species on daily transcript levels of individual GOIs and TNC levels (soluble sugars and starch). One and two-way ANOVAs were performed using the statistical package STATISTICA (StatSoft, Inc. v. 10, Brookline, MA, United States). Data normality was tested using the Shapiro–Wilk test, and the variance homogeneity was verified using Levene’s test or Cochran test. When assumptions of normality or homoscedasticity were not met, data were Box-Cox transformed. The Student–Newman–Keuls (SNK) post hoc test was used whenever significant differences were detected.
Results
Daily Cycle of Gene Expression in Cymodocea nodosa and Zostera marina
Principal component analysis and clustering analysis were used to separate samples collected at different times of the day and were based on multivariate gene-expression data. In C. nodosa, both approaches identified two main groups: a first including dawn, sunrise and solar noon samples and a second comprising sunset, dusk, and midnight samples (Figure 1A). Transcripts responsible for such separation are highlighted in the PCA biplot and their weightings are outlined in Supplementary Table 4. The most positively correlated transcripts with the PC1 (x-axis) were the circadian clock component LHY, and those involved in sucrose synthesis and transport (SPS1 and SUT2), whereas LHCA1 and FBA1 drove the separation of samples along the PC2 (y-axis). In Z. marina, PCA and clustering analysis identified three main groups: a first including sunrise, solar noon and sunset samples; a second including midnight and dawn samples; and a third comprising only dusk samples (Figure 1B). As evident from the biplot, the most positively correlated transcript with the PC1 was the transcript for the light harvesting protein LHCA1, while circadian clock and photoreception-related transcripts (i.e., LHY, CRY, and NPH1) were negatively correlated with the PC1. On the contrary, transcripts involved in photosynthesis (e.g., psbA, FD, and RBCS) and those involved in sugar metabolism (e.g., GWD and SPS1) were among the most positively correlated with the PC2. Overall, there was a significant difference in the diurnal pattern of accumulation of the targeted transcripts between the two seagrass species, as highlighted by the two-way PERMANOVA (Table 2 and Supplementary Table 5). A significant effect of time (Ti) (P < 0.05), species (Sp) (P < 0.001), as well as a significant Ti×Sp interaction (P < 0.001), was identified (Table 2). In further support of this, differences between species for each GOI based on two-way ANOVAs, are outlined in Supplementary Table 6.
Figure 1. Principal component analysis (PCA) with biplot and clustering based on multivariate gene-expression data in Cymodocea nodosa (A) and Zostera marina (B) (n = 3). Transcript loadings are outlined in Supplementary Table 4.
Table 2. Results of two-way PERMANOVA conducted on multivariate gene-expression data (–ΔCT values for all 19 Genes of Interests) in Cymodocea nodosa and Zostera marina.
Photosynthesis and Light Harvesting
Transcripts encoding for main subunits of Photosystem I (PSI) and II (PSII) (psbA and psbD), RuBisCO small subunit (RBCS) and ferredoxin (FD) in C. nodosa exhibited a very coordinate expression pattern, with max. abundance between solar noon and sunset, and a subsequent slow-down during the rest of the day (Figure 2). A similar trend was observed for transcripts encoding for light harvesting proteins, which exhibited max. abundance at noon and min. from midnight to sunrise (Figure 2). Among photosynthesis-related genes, psbD and FD were the ones significantly affected by the daily time according to the one-way ANOVA (psbD: sunset ≠ dusk; FD: solar noon ≠ dusk) (Table 3). Expression levels of LHCA1 was also significantly different between many time points during the day (Table 3). In Z. marina, photosynthesis-related genes, after a first peak at sunrise, decreased their abundance until solar noon, and then their level started to increase again until dusk (Figure 2). However, this pattern was only statistically significant for FD (Table 4). LHCA1 peaked between sunrise and sunset (Table 4).
Figure 2. Daily transcript levels (as normalized –ΔCT data) of Genes of Interests (GOIs) with a function in photosynthesis and light harvesting in C. nodosa (upper panels) and Z. marina (lower panels). Significant daily changes (P < 0.05) for each GOI based on SNK tests following one-way ANOVAs (Tables 3, 4), are indicated with different letters. Letters are coloured as the respective GOI. Values are mean ± SE for n = 3.
Table 3. Results of one-way ANOVAs based on –ΔCT data of individual Genes of Interests in Cymodocea nodosa.
Table 4. Results of one-way ANOVAs based on –ΔCT data of individual Genes of Interests in Zostera marina.
Sucrose and Starch Metabolism
In C. nodosa, transcripts encoding for proteins involved in sucrose synthesis and transport (SUT2 and SPS1) peaked at the beginning and during the light phase (i.e., at sunrise and solar noon), then their abundance decreased throughout the day, with a min. of expression at midnight (Figure 3). On the contrary, transcripts with a role in sucrose and starch metabolism globally tend to increase their expression levels at the end of the day, peaking during the night phase (i.e., from dusk to dawn) (Figure 3). In details, transcript abundance of SUT2 significantly varied between sunrise and midnight (P < 0.05), while SPS1 exhibited significant variations between sunrise and midnight (P < 0.05), and solar noon and midnight (P < 0.05) (Table 3). Sucrose synthase (SUS) displayed significant differences between solar noon and midnight (P < 0.05) (Table 3), while the other GOIs within this category did not exhibited significant changes. In Z. marina, SPS1 peaked later in the day (i.e., at sunset), while SUT2 showed no differences during the day (Figure 3 and Table 4). Genes related to sucrose and starch metabolism did not exhibited a very coordinate expression pattern, although their levels of abundance were also generally higher during dark hours. GWD exhibited significant changes at dawn, in respect to other time points (Figure 3 and Table 4).
Figure 3. Daily transcript levels (as normalized –ΔCT data) of GOIs with a function in sucrose and starch metabolism, and sucrose synthesis and transport in C. nodosa (upper panels) and Z. marina (lower panels). Significant daily changes (P < 0.05) for each GOI based on SNK tests following one-way ANOVAs (Tables 3, 4), are indicated with different letters. Letters are coloured as the respective GOI. Values are mean ± SE for n = 3.
Circadian Clock and Photoreception
In C. nodosa, the circadian clock component LHY displayed a peak of abundance in the early morning (i.e., from dawn to sunrise) and a min. of expression between sunset and midnight (Figure 4 and Table 3). ZTL, as well as blue-light photoreceptors (NPH1 and CRY1) showed no significant expression differences among the selected time points (Figure 4). In Z. marina, the circadian-clock component LHY reached its max. level slightly earlier than C. nodosa, between midnight and dawn (Figure 4 and Table 4). Blue-light photoreceptors exhibited a similar pattern, as they showed max. expression levels between midnight and dawn (Figure 4). According to one-way ANOVA, levels of LHY differed significantly at midnight and dawn in respect to the rest of the day (Table 4), while NPH1 and CRY1 varied significantly their expression between solar noon and dawn (P < 0.05 for NPH1 and P = 0.06 for CRY1) (Table 4).
Figure 4. Daily transcript levels (as normalized –ΔCT data) of GOIs with a function in circadian clock and photoreception in C. nodosa (upper panels) and Z. marina (lower panels). Significant daily changes (P < 0.05) for each GOI based on SNK tests following one-way ANOVAs (Tables 3, 4), are indicated with different letters. Letters are coloured as the respective GOI. Values are mean ± SE for n = 3.
Mitochondrial Respiration and Growth
Regarding respiration-related genes, the transcript for Cytochrome C (COX5B) showed no significant expression changes among time points in both species. The transcript encoding for AOX1A instead exhibited significant oscillation during the day in C. nodosa, with max. abundance around dusk (Figure 5 and Table 3). The transcript for Expansin B (EB) exhibited no significant daily changes in both species.
Figure 5. Daily transcript levels (as normalized –ΔCT data) of GOIs with a function in cellular respiration and growth in C. nodosa (upper panels) and Z. marina (lower panels). Significant daily changes (P < 0.05) for each GOI based on SNK tests following one-way ANOVAs (Tables 3, 4), are indicated with different letters. Letters are coloured as the respective GOI. Values are mean ± SE for n = 3.
Daily Patterns of Carbohydrate Accumulation in Cymodocea nodosa and Zostera marina
Two-way ANOVAs revealed a significant difference of soluble sugar accumulation between the two species (P < 0.05; Supplementary Table 7). Daily levels of starch accumulation were affected by both species (Sp) and time (Ti) (Supplementary Table 7). In C. nodosa, starch level peaked at sunrise (sunrise ≠ dawn, solar noon, sunset, midnight), in respect to the rest of the day (Figure 6 and Table 5). In Z. marina, starch abundance was almost constant during the day (Figure 6).
Figure 6. Daily patterns of TNC (soluble sugars and starch) accumulation in C. nodosa and Z. marina. Significant daily changes (P < 0.05) based on SNK tests following one-way ANOVAs (Table 5), are indicated with different letters. Values are mean ± SE for n = 3.
Table 5. Results of one-way ANOVAs to assess the effects of time of the day on TNC accumulation (soluble sugars and starch) in Zostera marina and Cymodocea nodosa.
Discussion
General Overview of Daily Timing in Cymodocea nodosa and Zostera marina
Here, we provided a first description of daily regulation of key pathways in two seagrass species, C. nodosa and Z. marina, co-occurring at the same geographic location, thus exposed to identical natural variations in light and temperature. In both species, genes involved in sucrose metabolism (e.g., SPS1) and core circadian clock components (i.e., LHY), in addition to those encoding for light harvesting proteins (i.e., LHCA1) or involved in redox regulation (e.g., FD) exhibited clear diurnal rhythms. Such genes were also those making a greater contribution to sample distribution in the PCA (Figure 1 and Supplementary Table 4).
Despite this common behavior, the variations we observed in their daily transcriptional timing, were species-specific (Table 2 and Supplementary Tables 5, 6), highlighting the importance of intrinsic biological, and likely ecological, attributes in determining the periodicity of functions that go beyond the (common) effects of external stimuli (Müller et al., 2014; Ruocco et al., 2020a). Five GOIs (i.e., SPS1, GWD, FD, LHY and LHCA1) showed a significant Ti×Sp interaction, indicating their daily pattern of expression was substantially different between the two seagrass species (Supplementary Table 6).
In C. nodosa, a clear separation occurred between “light” (i.e., dawn, sunrise and solar noon) and “dark” (i.e., dusk, sunset and midnight) hour responses, where solar noon and midnight represented the most different time points. Yet, similar gene-expression responses were observed during light-to-dark and dark-to-light transition times. On the other hand, in Z. marina, “light-hour responses” were shifted in time, as sunrise, solar noon and dusk clustered together, while dawn-related responses were closer to those occurring at midnight. In general, in Z. marina, larger differences were detected during transition times (dusk vs. dawn), in comparison to C. nodosa (Supplementary Table 5).
Globally, the most frequent time for the max. expression (average transcript level) of GOIs was the middle of the day, followed by midnight/dawn in C. nodosa, while the min. of expression was recorded at sunset. In Z. marina instead, max. average expression occurred later, at the end of the day (i.e., during sunset/dusk), and minimum at dawn. This might suggest a differential diurnal control of growth rates between the two species, which could provide them with physiological and/or ecological benefits (Müller et al., 2014). The regulation of plant growth generally occurs through the control of carbon metabolism and light signaling, both processes being coordinated by distinct components of the circadian clock (Nozue and Maloof, 2006; Graf and Smith, 2011; Dodd et al., 2014). In Arabidopsis, growth rates peak toward the end of the night when plants are not photosynthetically active and have to properly manage metabolic resources until the next morning (Wiese et al., 2007). This regulation is also achieved via LHY/CCA1, which could have a role in ensuring the correct timing of starch breakdown at night, and ELF3 acting to repress growth in the light, thus favoring the accumulation of carbohydrate reserves (Graf et al., 2010; Nusinow et al., 2011). In our experiment, we do not have direct measurements of diurnal growth rates in C. nodosa and Z. marina, and this has never been assessed in seagrasses. However, the shifted peaks of LHY abundance and the slightly different daily patterns of carbohydrate accumulation between species (see below), could point for a different timing of growth based on a different strategy of mobilization and allocation of resources at night, coordinated by the circadian clock (Wiese et al., 2007; Graf and Smith, 2011).
Specific Daily Regulation of Individual Genes of Interests in Cymodocea nodosa and Zostera marina
In C. nodosa, low sugar levels occurring by the end of the night, could be responsible for the induction of transcripts involved in sucrose synthesis and transport in the early morning (starting from dawn), such as SPS1 and SUT2, indicating that carbohydrates have declined to critical levels, and plants have become net consumers of carbon (Bläsing et al., 2005). Subsequent to the induction of such genes, leaf soluble sugars peaked in the light, between sunrise and solar noon, demonstrating an activation of the related biochemical pathways. On the contrary, transcripts encoding proteins responsible for sucrose breakdown e.g., SUS, were induced during the night, similar to transcripts involved in starch mobilization. Glucan water dikinase, a key enzyme regulating nocturnal starch breakdown in leaves, and other starch-mobilizing enzymes (e.g., Beta-amylase and Fructose-bisphosphate aldolase) revealed cycling transcription in A. thaliana under a LD cycle and free-running conditions (Harmer et al., 2000; Lu et al., 2005) concomitant to cycling of starch breakdown. In Z. marina, diurnal rhythms of transcripts related to sugar metabolism were shifted in time in respect to what observed in C. nodosa. For instance, SPS1, involved in sucrose synthesis, peaked later in the day (between solar noon and sunset). On the contrary, transcripts encoding enzymes responsible for sugar and starch degradation/mobilization (e.g., GWD) peaked slightly earlier than C. nodosa. Endogenous sugar clock entrainment optimizes plant growth and fitness regulating carbon homeostasis between source and sink tissues (Ohara and Satake, 2017). The observed differences in gene expression could be related to different patterns of photosynthate accumulation/translocation between species: during the day, C. nodosa could quickly use triose-phosphates produced by the Calvin cycle for sucrose synthesis to support leaf maintenance and export everything else to rhizomes, where sugars can be used for rhizome expansion or stored as starch until required; on the other hand, Z. marina seems to invest more in the leaves, accumulating (exporting less) sugars and starch during the day and using starch during the night to support sucrose synthesis and export (Geiger et al., 2000).
All the above highlights the importance that sugar levels could have for diurnal gene regulation in seagrasses (Contento et al., 2004; Bläsing et al., 2005). Yet, the rate of starch degradation at night, which is known to be under control of the circadian clock (Müller et al., 2014), could also be of fundamental importance for the carbon economy of seagrass species with large effects on their short and long-term productivity and deserve further investigation (Sulpice et al., 2009; Graf and Smith, 2011).
In C. nodosa, transcripts involved in photosynthetic carbon fixation, as well as light harvesting proteins, were triggered by light, and co-ordinately peaked in the middle of the day, thus supporting the synthesis and export of sucrose needed for plant metabolism and storage. This pattern of accumulation of photosynthesis-related transcripts is very similar to the one previously detected in shallow-growing P. oceanica plants (Procaccini et al., 2017) and matches what is observed in Arabidopsis under constant light conditions, where photosynthesis genes peaked near the middle of the subjective day (Harmer et al., 2000). However, the diurnal pattern of accumulation of such transcripts was one of the most striking differences observed between the two analyzed species. Indeed, in Z. marina, transcripts for photosynthesis-related genes were slightly induced at sunrise, and then there was a down-regulation around noon followed by a further increase between sunset and dusk. We do not have a conclusive explanation for this pattern, however photo-physiological data collected in the same experiment indicated that one of the most noticeable difference between C. nodosa and Z. marina was their susceptibility to high irradiance levels at noon (Hofman, 2019). Specifically, Z. marina was much less affected by high light than C. nodosa, whose effective quantum yield (ΔF/Fm′) dropped significantly more at noon and exhibited a much higher non-photochemical quenching (NPQ) (Hofman, 2019). Besides, the max. abundance level of AOX1A, which is involved in energy dissipation mechanisms via alternative mitochondrial respiration, also differed across species, peaking at dusk in C. nodosa and at noon in Z. marina (ns). This highlights that the two species could adopt different strategies for dissipating excess of energy throughout the day, both during photosynthesis and respiration processes, and this could be somehow linked to the differential transcriptional timing we observed (e.g., for psbA) (Phee et al., 2004). However, studies at the protein level should be conducted, as transcriptome and proteome could have very different scales of regulation during the day (Graf et al., 2017).
As discussed above, diurnal rhythms observed in C. nodosa and Z. marina could be regulated by sugar availability and the circadian clock. At this regard, LHY could play a key role, as its diurnal changes are highly significant in both species. Our data confirmed that the peak of abundance of this transcript occurs at the beginning of the day (dawn/sunrise) in C. nodosa, as observed in many terrestrial species (McWatters and Devlin, 2011; Haydon et al., 2013b; Dodd et al., 2015), while in Z. marina its peak was slightly earlier (midnight/dawn). Around dawn, the phase of the circadian oscillator is indeed adjusted in response to low light intensity detected by photoreceptors as the sun rises, entraining primary metabolism, before the so-called “metabolic dawn” that occurs later due to the increasing concentration of sugars (Haydon et al., 2013a,b; Dodd et al., 2015).
The anticipation in timing of the peak in expression of LHY in Z. marina before dawn was quite unexpected. Differences in diurnal rhythms of circadian clocks components (i.e., LHY/CCA1 and TOC1) are known to occur across ecotypes (Slotte et al., 2007). Shifts in phase of LHY/CCA1 have been tracked in A. thaliana mutants exposed to artificial light conditions (Mizoguchi et al., 2002), while changes in the transcriptional peaks of several core clock components were observed in A. thaliana plants growing under different photoperiods (Flis et al., 2016). Besides the photoperiod cues, dynamic adjustment of phase and period of the plant core clock are induced in response to sugars and several biotic and abiotic stimuli (Webb et al., 2019). Further studies under controlled conditions should be performed to investigate which are the cues that are more important for controlling the circadian responses of the different seagrass species.
Transcripts encoding blue photoreceptors CRY1 and the NPH1 significantly oscillated during the day in Z. marina, where they had a very similar expression pattern, with max. expression levels at the beginning of the light phase (i.e., dawn). This is similar to terrestrial model land plants, where diurnal rhythms of mRNA levels of photoreceptors (both cryptochromes and phytochromes) exhibit their max. in the light phase, anticipating lights-on and lights-off signals (Tóth et al., 2001), and inducing change of gene expression in both nuclear and plastid transcripts (Wang et al., 2014).
Conclusion
Here, we provided a first overview of daily transcriptional timing of key metabolic pathways in two marine plants exposed to the same photoperiod. A strength of our study is that plants were exposed to natural irradiance conditions, with irregular fluctuations and gradual LD transitions at dusk and dawn. Experiments with terrestrial species are typically conducted exposing plants in controlled environments with constant (artificial) daily irradiance and abrupt LD transitions. Such studies are often not representative of natural settings, and diurnal variations of carbon metabolism can differ significantly between sunlight and artificial light conditions (Vialet-Chabrand et al., 2016, 2017; Annunziata et al., 2017; Panter et al., 2019).
Overall, our results suggested a potential interplay between circadian clock and sugar levels in determining the periodicity of functions in both species, as demonstrated in terrestrial plants, although other factors could be at play. Large differences were observed between species in the daily regulation of transcripts related to photosynthesis and sugar metabolism, which together with the slightly different daily patterns of carbohydrate accumulation, highlight the importance of intrinsic biological and likely ecological attributes of the species. The two species could adopt e.g., different growth timing based on a differential strategy of resource allocation and mobilization, coordinated by the circadian clock. Such behavior could also derive from divergent evolutionary adaptations of the species to different photoperiodic conditions across their bio-geographical range of distributions. This could also partially explains the differential sensitivity of the two species during light-to-dark and dark-to-light transitions. Species-specific effects of twilight duration on circadian phases of animal behavior (Boulos et al., 1996) and physiology (Sosniyenko et al., 2009) have been well documented, as well as the occurrence of circadian phenotypes correlated with latitude and geographical distribution clines either in wild plant populations or within selectively bred crops (Greenham et al., 2017; Rees et al., 2021).
We are aware that our study has limitations due to a single day analysis. In addition, we can only infer, e.g., the circadian regulation of certain patterns, as this should be confirmed with ad hoc experiments under constant light or dark conditions. However, we believe our analysis addressing the diurnal response of two species along a single day is informative and broadens current knowledge of biological rhythms in seagrasses. Studies of circadian rhythms in marine plants are currently missing, although these could have wider implications, as the modulation of such rhythms could be used for improving marine plant performances and productivity (Kim et al., 2017; Steed et al., 2021) and may ultimately help in buffering the effects of climate changes on seagrass distribution and persistence (Resco et al., 2009).
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
ED and MR conceived the ideas and designed the study. ED, MR, JS, IB, JH, and KP participated in sample collection and mesocosm maintenance. MR and ED performed all gene-expression and biochemical analyses. MR led the writing of the manuscript. All authors contributed in drafting the work, revised it critically, and gave final approval for publication.
Funding
The research leading to these results received funding from the European Union’s Horizon 2020 Research and Innovation Programme under grant agreement no 730984, ASSEMBLE Plus project. This study also received Portuguese national funds from FCT - Foundation for Science and Technology through the projects GRASSMET (PTDC/MAR-EST/4257/2014) and UIDB/04326/2020. The work has been partially supported by the project Marine Hazard, PON03PE_00203_1, Italian Ministry of Education, University and Research (MIUR).
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.
Acknowledgments
We thank the staff of CCMAR (Faro, Portugal) for the hospitality and in particular all the staff at the Ramalhete Experimental Station.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2021.757187/full#supplementary-material
Footnotes
References
Altschul, S. F., Gish, W., Miller, W., Myers, E. W., and Lipman, D. J. (1990). Basic local alignment search tool. J. Mol. Biol. 215, 403–410.
Anisimova, M., and Gascuel, O. (2006). Approximate likelihood-ratio test for branches: a fast, accurate, and powerful alternative. Syst. Biol. 55, 539–552. doi: 10.1080/10635150600755453
Annunziata, M. G., Apelt, F., Carillo, P., Krause, U., Feil, R., Mengin, V., et al. (2017). Getting back to nature: a reality check for experiments in controlled environments. J. Exp. Bot. 68, 4463–4477. doi: 10.1093/jxb/erx220
Barak, S., Tobin, E. M., Green, R. M., Andronis, C., and Sugano, S. (2000). All in good time: the Arabidopsis circadian clock. Trends Plant Sci. 5, 517–522.
Bläsing, O. E., Gibon, Y., Günther, M., Höhne, M., Morcuende, R., Osuna, D., et al. (2005). Sugars and circadian regulation make major contributions to the global regulation of diurnal gene expression in Arabidopsis. Plant Cell 17, 3257–3281. doi: 10.1105/tpc.105.035261
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.
Boulos, Z., Macchi, M., and Terman, M. (1996). Effects of twilights on circadian entrainment patterns and reentrainment rates in squirrel monkeys. J. Comp. Physiol. A 179, 687–694. doi: 10.1007/BF00216132
Capella-Gutiérrez, S., Silla-Martínez, J. M., and Gabaldón, T. (2009). trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25, 1972–1973. doi: 10.1093/bioinformatics/btp348
Contento, A. L., Kim, S.-J., and Bassham, D. C. (2004). Transcriptome profiling of the response of Arabidopsis suspension culture cells to Suc starvation. Plant Physiol. 135, 2330–2347. doi: 10.1104/pp.104.044362
Covington, M. F., Maloof, J. N., Straume, M., Kay, S. A., and Harmer, S. L. (2008). Global transcriptome analysis reveals circadian regulation of key pathways in plant growth and development. Genome Biol. 9, 1–18.
Dodd, A. N., Belbin, F. E., Frank, A., and Webb, A. A. R. (2015). Interactions between circadian clocks and photosynthesis for the temporal and spatial coordination of metabolism. Front. Plant Sci. 6:245. doi: 10.3389/fpls.2015.00245
Dodd, A. N., Dalchau, N., Gardner, M. J., Baek, S. J., and Webb, A. A. (2014). The circadian clock has transient plasticity of period and is required for timing of nocturnal processes in Arabidopsis. New Phytol. 201, 168–179. doi: 10.1111/nph.12489
Feugier, F. G., and Satake, A. (2013). Dynamical feedback between circadian clock and sucrose availability explains adaptive response of starch metabolism to various photoperiods. Front. Plant Sci. 3:305. doi: 10.3389/fpls.2012.00305
Filichkin, S. A., Breton, G., Priest, H. D., Dharmawardhana, P., Jaiswal, P., Fox, S. E., et al. (2011). Global profiling of rice and poplar transcriptomes highlights key conserved circadian-controlled pathways and cis-regulatory modules. PLoS One 6:e16907. doi: 10.1371/journal.pone.0016907.g010
Flis, A., Mengin, V., Ivakov, A. A., Mugford, S. T., Hubberten, H. M., Encke, B., et al. (2019). Multiple circadian clock outputs regulate diel turnover of carbon and nitrogen reserves. Plant Cell Environ. 42, 549–573. doi: 10.1111/pce.13440
Flis, A., Sulpice, R., Seaton, D. D., Ivakov, A. A., Liput, M., Abel, C., et al. (2016). Photoperiod-dependent changes in the phase of core clock transcripts and global transcriptional outputs at dawn and dusk in Arabidopsis. Plant Cell Environ. 39, 1955–1981. doi: 10.1111/pce.12754
Geiger, D. R., Servaites, J. C., and Fuchs, M. A. (2000). Role of starch in carbon translocation and partitioning at the plant level. Funct. Plant Biol. 27, 571–582. doi: 10.1046/j.1469-8137.1997.00843.x
Graf, A., Coman, D., Uhrig, R. G., Walsh, S., Flis, A., Stitt, M., et al. (2017). Parallel analysis of Arabidopsis circadian clock mutants reveals different scales of transcriptome and proteome regulation. Open Biol. 7:160333. doi: 10.1098/rsob.160333
Graf, A., Schlereth, A., Stitt, M., and Smith, A. M. (2010). Circadian control of carbohydrate availability for growth in Arabidopsis plants at night. Proc. Natl. Acad. Sci. U.S.A. 107, 9458–9463. doi: 10.1073/pnas.0914299107
Graf, A., and Smith, A. M. (2011). Starch and the clock: the dark side of plant productivity. Trends Plant Sci. 16, 169–175. doi: 10.1016/j.tplants.2010.12.003
Greenham, K., Lou, P., Puzey, J. R., Kumar, G., Arnevik, C., Farid, H., et al. (2017). Geographic variation of plant circadian clock function in natural and agricultural settings. J. Biol. Rhythms 32, 26–34. doi: 10.1177/0748730416679307
Guindon, S., Dufayard, J. F., Lefort, V., Anisimova, M., Hordijk, W., and Gascuel, O. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307–321. doi: 10.1093/sysbio/syq010
Hammer, Ø, Harper, D., and Ryan, P. (2001). PAST: paleontological statistics software package for education and data analysis. Palaeontol. Electron. 4, 1–9.
Harmer, S., Hogenesch, J., Straume, M., Chang, H., Han, B., Zhu, T., et al. (2000). Orchestrated transcription of key pathways in Arabidopsis by the circadian clock. Science 290, 2110–2113. doi: 10.1126/science.290.5499.2110
Harmer, S. L. (2009). The circadian system in higher plants. Annu. Rev. Plant Biol. 60, 357–377. doi: 10.1146/annurev.arplant.043008.092054
Haydon, M. J., Hearn, T. J., Bell, L. J., Hannah, M. A., and Webb, A. A. (2013a). Metabolic regulation of circadian clocks. Semin. Cell Dev. Biol. 24, 414–421.
Haydon, M. J., Mielczarek, O., Robertson, F. C., Hubbard, K. E., and Webb, A. A. R. (2013b). Photosynthetic entrainment of the Arabidopsis thaliana circadian clock. Nature 502, 689–692. doi: 10.1038/nature12603
Hofman, J. D. (2019). Compared Physiological Performances of Caulerpa prolifera and Native Seagrasses of Ria Formosa. Faro: Universidade Do Algarve.
Hubbard, C. J., Brock, M. T., Van Diepen, L. T., Maignien, L., Ewers, B. E., and Weinig, C. (2018). The plant circadian clock influences rhizosphere community structure and function. ISME J. 12, 400–410.
Huffeldt, N. P. (2020). Photic barriers to poleward range-shifts. Trends Ecol. Evol. 35, 652–655. doi: 10.1016/j.tree.2020.04.011
Jończyk, M., Sobkowiak, A., Siedlecki, P., Biecek, P., Trzcinska-Danielewicz, J., Tiuryn, J., et al. (2011). Rhythmic diel pattern of gene expression in juvenile maize leaf. PLoS One 6:e23628. doi: 10.1371/journal.pone.0023628
Jones, M. A. (2009). Entrainment of the Arabidopsis circadian clock. J. Plant Biol. 52, 202–209. doi: 10.1007/s12374-009-9030-1
Kim, J. A., Kim, H.-S., Choi, S.-H., Jang, J.-Y., Jeong, M.-J., and Lee, S. I. (2017). The importance of the circadian clock in regulating plant metabolism. Int. J. Mol. Sci. 18:2680. doi: 10.3390/ijms18122680
Kirk, J. T. O. (2010). Light and Photosynthesis in Aquatic Ecosystems. Cambridge: Cambridge University Press.
Larkum, A. W., Pernice, M., Schliep, M., Davey, P., Szabo, M., Raven, J. A., et al. (2018). “Photosynthesis and metabolism of seagrasses,” in. Seagrasses of Australia, eds A. Larkum, G. Kendrick, and P. Ralph (Cham: Springer), 315–342. doi: 10.1007/978-3-319-71354-0_11
Larkum, A. W. D., Orth, R. J., and Duarte, C. M. (2006). Seagrasses: Biology, Ecology and Conservation. Amsterdam: Springer.
Les, D. H., Cleland, M. A., and Waycott, M. (1997). Phylogenetic studies in Alismatidae, ii: evolution of marine angiosperms (seagrasses) and hydrophily. Syst. Bot. 22, 443–463. doi: 10.2307/2419820
Litthauer, S., Battle, M. W., and Jones, M. A. (2016). Phototropins do not alter accumulation of evening-phased circadian transcripts under blue light. Plant Signal. Behav. 11:e1126029. doi: 10.1080/15592324.2015.1126029
Lu, H., Mcclung, C. R., and Zhang, C. (2017). Tick tock: circadian regulation of plant innate immunity. Annu. Rev. Phytopathol. 55, 287–311. doi: 10.1146/annurev-phyto-080516-035451
Lu, Y., Gehan, J. P., and Sharkey, T. D. (2005). Daylength and circadian effects on starch degradation and maltose metabolism. Plant Physiol. 138, 2280–2291. doi: 10.1104/pp.105.061903
Marbà, N., Arias-Ortiz, A., Masqué, P., Kendrick, G. A., Mazarrasa, I., Bastyan, G. R., et al. (2015). Impact of seagrass loss and subsequent revegetation on carbon sequestration and stocks. J. Ecol. 103, 296–302. doi: 10.1111/1365-2745.12370
Marín-Guirao, L., Ruiz, J. M., Dattolo, E., Garcia-Munoz, R., and Procaccini, G. (2016). Physiological and molecular evidence of differential short-term heat tolerance in Mediterranean seagrasses. Sci. Rep. 6, 1–13. doi: 10.1038/srep28615
McWatters, H. G., and Devlin, P. F. (2011). Timing in plants–a rhythmic arrangement. FEBS Lett. 585, 1474–1484. doi: 10.1016/j.febslet.2011.03.051
Michael, T., Mockler, T., Breton, G., Mcentee, C., Byer, A., Trout, J., et al. (2008). Network discovery pipeline elucidates conserved time-of-day-specific cis-regulatory modules. PLoS Genet. 4:e14. doi: 10.1371/journal.pgen.0040014
Millar, A. J. (2003). A suite of photoreceptors entrains the plant circadian clock. J. Biol. Rhythms 18, 217–226. doi: 10.1177/0748730403018003004
Millar, A. J. (2004). Input signals to the plant circadian clock. J. Exp. Bot. 55, 277–283. doi: 10.1093/jxb/erh034
Mizoguchi, T., Wheatley, K., Hanzawa, Y., Wright, L., Mizoguchi, M., Song, H. R., et al. (2002). LHY and CCA1 are partially redundant genes required to maintain circadian rhythms in Arabidopsis. Dev. Cell 2, 629–641. doi: 10.1016/s1534-5807(02)00170-3
Müller, L. M., Von Korff, M., and Davis, S. J. (2014). Connections between circadian clocks and carbon metabolism reveal species-specific effects on growth control. J. Exp. Bot. 65, 2915–2923. doi: 10.1093/jxb/eru117
Nozue, K., and Maloof, J. N. (2006). Diurnal regulation of plant growth. Plant Cell Environ. 29, 396–408.
Nusinow, D. A., Helfer, A., Hamilton, E. E., King, J. J., Imaizumi, T., Schultz, T. F., et al. (2011). The ELF4–ELF3–LUX complex links the circadian clock to diurnal control of hypocotyl growth. Nature 475, 398–402. doi: 10.1038/nature10182
Ohara, T., and Satake, A. (2017). Photosynthetic entrainment of the circadian clock facilitates plant growth under environmental fluctuations: perspectives from an integrated model of phase oscillator and phloem transportation. Front. Plant Sci. 8:1859. doi: 10.3389/fpls.2017.01859
Okonechnikov, K., Golosova, O., Fursov, M., and Ugene Team. (2012). Unipro UGENE: a unified bioinformatics toolkit. Bioinformatics 28, 1166–1167. doi: 10.1093/bioinformatics/bts091
Olivé, I., Silva, J., Lauritano, C., Costa, M. M., Ruocco, M., Procaccini, G., et al. (2017). Linking seagrass gene expression to productivity: long and short-term responses of Cymodocea nodosa to CO2 in volcanic vents. Sci. Rep. 7:42278.
Olsen, J. L., Rouzé, P., Verhelst, B., Lin, Y.-C., 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
Ouyang, Y., Andersson, C. R., Kondo, T., Golden, S. S., and Johnson, C. H. (1998). Resonating circadian clocks enhance fitness in cyanobacteria. Proc. Natl. Acad. Sci. U.S.A. 95, 8660–8664.
Panter, P. E., Muranaka, T., Cuitun-Coronado, D., Graham, C. A., Yochikawa, A., Kudoh, H., et al. (2019). Circadian regulation of the plant transcriptome under natural conditions. Front. Genet. 10:1239. doi: 10.3389/fgene.2019.01239
Papadopoulos, J. S., and Agarwala, R. (2007). COBALT: constraint-based alignment tool for multiple protein sequences. Bioinformatics 23, 1073–1079. doi: 10.1093/bioinformatics/btm076
Phee, B. K., Cho, J. H., Park, S., Jung, J. H., Lee, Y. H., Jeon, J. S., et al. (2004). Proteomic analysis of the response of Arabidopsis chloroplast proteins to high light stress. Proteomics 4, 3560–3568. doi: 10.1002/pmic.200400982
Procaccini, G., Ruocco, M., Marín-Guirao, L., Dattolo, E., Brunet, C., D’esposito, D., et al. (2017). Depth-specific fluctuations of gene expression and protein abundance modulate the photophysiology in the seagrass Posidonia oceanica. Sci. Rep. 7:42890. doi: 10.1038/srep42890
Ransbotyn, V., and Reusch, T. B. H. (2006). Housekeeping gene selection for quantitative real-time PCR assays in the seagrass Zostera marina subjected to heat stress. Limnol. Oceanogr. Methods 4, 367–373. doi: 10.4319/lom.2006.4.367
Rasmusson, L. M., Lauritano, C., Procaccini, G., Gullström, M., Buapet, P., and Björk, M. (2017). Respiratory oxygen consumption in the seagrass Zostera marina varies on a diel basis and is partly affected by light. Mar. Biol. 164:140. doi: 10.1007/s00227-017-3168-z
Rees, H., Joynson, R., Brown, J. K., and Hall, A. (2021). Naturally occurring circadian rhythm variation associated with clock gene loci in Swedish Arabidopsis accessions. Plant Cell Environ. 44, 807–820. doi: 10.1111/pce.13941
Resco, V., Hartwell, J., and Hall, A. (2009). Ecological implications of plants’ ability to tell the time. Ecol. Lett. 12, 583–592. doi: 10.1111/j.1461-0248.2009.01295.x
Resco de Dios, V. (2013). On the persistence of memory: soft clocks and terrestrial biosphere-atmosphere interactions. Plant Signal. Behav. 8, 743–752. doi: 10.4161/psb.26964
Resco de Dios, V., Gessler, A., Ferrio, J. P., Alday, J. G., Bahn, M., Del Castillo, J., et al. (2016). Circadian rhythms have significant effects on leaf-to-canopy scale gas exchange under field conditions. Gigascience 5:43. doi: 10.1186/s13742-016-0149-y
Ruiz-Frau, A., Gelcich, S., Hendriks, I., Duarte, C. M., and Marbà, N. (2017). Current state of seagrass ecosystem services: research and policy integration. Ocean Coast. Manag. 149, 107–115. doi: 10.1016/j.scitotenv.2020.142565
Ruocco, M., Ambrosino, L., Jahnke, M., Chiusano, M. L., Barrote, I., Procaccini, G., et al. (2020a). m6A RNA methylation in marine plants: first insights and relevance for biological rhythms. Int. J. Mol. Sci. 21:7508. doi: 10.3390/ijms21207508
Ruocco, M., Entrambasaguas, L., Dattolo, E., Milito, A., Marín-Guirao, L., and Procaccini, G. (2020b). A king and vassals’ tale: molecular signatures of clonal integration in Posidonia oceanica under chronic light shortage. J. Ecol. 109, 294–312.
Ruocco, M., Musacchia, F., Olivé, I., Costa, M. M., Barrote, I., Santos, R., et al. (2017). Genomewide transcriptional reprogramming in the seagrass Cymodocea nodosa under experimental ocean acidification. Mol. Ecol. 26, 4241–4259. doi: 10.1111/mec.14204
Ruocco, M. Marín-Guirao, L., and Procaccini, G. (2019). 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:24. doi: 10.1007/s00227-019-3482-8
Sánchez, R., Serra, F., Tárraga, J., Medina, I., Carbonell, J., Pulido, L., et al. (2011). Phylemon 2.0: a suite of web-tools for molecular evolution, phylogenetics, phylogenomics and hypotheses testing. Nucleic Acids Res. 39(suppl_2), W470–W474. doi: 10.1093/nar/gkr408
Schmittgen, T. D., and Livak, K. J. (2008). Analyzing real-time PCR data by the comparative CT method. Nat. Protoc. 3, 1101–1108.
Scialdone, A., Mugford, S. T., Feike, D., Skeffington, A., Borrill, P., Graf, A., et al. (2013). Arabidopsis plants perform arithmetic division to prevent starvation at night. Elife 2:e00669. doi: 10.7554/eLife.00669
Seki, M., Ohara, T., Hearn, T. J., Frank, A., Da Silva, V. C., Caldana, C., et al. (2017). Adjustment of the Arabidopsis circadian oscillator by sugar signalling dictates the regulation of starch metabolism. Sci. Rep. 7:8305. doi: 10.1038/s41598-017-08325-y
Serikawa, M., Miwa, K., Kondo, T., and Oyama, T. (2008). Functional conservation of clock-related genes in flowering plants: overexpression and RNA interference analyses of the circadian rhythm in the monocotyledon Lemna gibba. Plant Physiol. 146, 1952–1963. doi: 10.1104/pp.107.114611
Short, F., Carruthers, T., Dennison, W., and Waycott, M. (2007). Global seagrass distribution and diversity: a bioregional model. J. Exp. Mar. Biol. Ecol. 350, 3–20.
Slotte, T., Holm, K., McIntyre, L. M., Lagercrantz, U., and Lascoux, M. (2007). Differential expression of genes important for adaptation in Capsella bursa-pastoris (Brassicaceae). Plant Physiol. 145, 160–173.
Song, Y. H., Estrada, D. A., Johnson, R. S., Kim, S. K., Lee, S. Y., Maccoss, M. J., et al. (2014). Distinct roles of FKF1, Gigantea, and Zeitlupe proteins in the regulation of Constans stability in Arabidopsis photoperiodic flowering. Proc. Natl. Acad. Sci. U.S.A. 111, 17672–17677. doi: 10.1073/pnas.1415375111
Sosniyenko, S., Hut, R. A., Daan, S., and Sumová, A. (2009). Influence of photoperiod duration and light–dark transitions on entrainment of Per1 and Per2 gene and protein expression in subdivisions of the mouse suprachiasmatic nucleus. Eur. J. Neurosci. 30, 1802–1814. doi: 10.1111/j.1460-9568.2009.06945.x
Steed, G., Ramirez, D. C., Hannah, M. A., and Webb, A. A. (2021). Chronoculture, harnessing the circadian clock to improve crop yield and sustainability. Science 372:eabc9141. doi: 10.1126/science.abc9141
Sulpice, R., Pyl, E.-T., Ishihara, H., Trenkamp, S., Steinfath, M., Witucka-Wall, H., et al. (2009). Starch as a major integrator in the regulation of plant growth. Proc. Natl. Acad. Sci. U.S.A. 106, 10348–10353. doi: 10.1073/pnas.0903478106
Tóth, R., Kevei, E., Hall, A., Millar, A. J., Nagy, F., and Kozma-Bognár, L. (2001). Circadian clock-regulated expression of phytochrome and cryptochrome genes in Arabidopsis. Plant Physiol. 127, 1607–1616. doi: 10.1104/pp.010467
Tutar, O. Marín-Guirao, L., Ruiz, J. M., and Procaccini, G. (2017). Antioxidant response to heat stress in seagrasses. A gene expression study. Mar. Environ. Res. 132, 94–102. doi: 10.1016/j.marenvres.2017.10.011
Untergasser, A., Cutcutache, I., Koressaar, T., Ye, J., Faircloth, B. C., Remm, M., et al. (2012). Primer3—new capabilities and interfaces. Nucleic Acids Res. 40, e115–e115. doi: 10.1093/nar/gks596
Vialet-Chabrand, S., Matthews, J., Brendel, O., Blatt, M., Wang, Y., Hills, A., et al. (2016). Modelling water use efficiency in a dynamic environment: an example using Arabidopsis thaliana. Plant Sci. 251, 65–74. doi: 10.1016/j.plantsci.2016.06.016
Vialet-Chabrand, S., Matthews, J. S. A., Simkin, A. J., Raines, C. A., and Lawson, T. (2017). Importance of fluctuations in light on plant photosynthetic acclimation. Plant Physiol. 173, 2163. doi: 10.1104/pp.16.01767
Wang, X., Wang, Q., Nguyen, P., and Lin, C. (2014). Cryptochrome-mediated light responses in plants. Enzymes 35, 167–189. doi: 10.1016/b978-0-12-801922-1.00007-5
Webb, A. A., and Satake, A. (2015). Understanding circadian regulation of carbohydrate metabolism in Arabidopsis using mathematical models. Plant Cell Physiol. 56, 586–593. doi: 10.1093/pcp/pcv033
Webb, A. A., Seki, M., Satake, A., and Caldana, C. (2019). Continuous dynamic adjustment of the plant circadian oscillator. Nat. Commun. 10:550.
Wiese, A., Christ, M., Virnich, O., Schurr, U., and Walter, A. (2007). Spatio-temporal leaf growth patterns of Arabidopsis thaliana and evidence for sugar control of the diel leaf growth cycle. New Phytol. 174, 752–761. doi: 10.1111/j.1469-8137.2007.02053.x
Winters, G., Nelle, P., Fricke, B., Rauch, G., and Reusch, T. B. H. (2011). Effects of a simulated heat wave on photophysiology and gene expression of high- and low-latitude populations of Zostera marina. Mar. Ecol. Prog. Ser. 435, 83–95. doi: 10.3354/meps09213
Yue, S., Zhang, X., Xu, S., Zhang, Y., Zhao, P., Wang, X., et al. (2020). Reproductive strategies of the seagrass Zostera japonica under different geographic conditions in northern China. Front. Mar. Sci. 7:990. doi: 10.3389/fmars.2020.574790
Keywords: gene expression, sugars, circadian clock, marine plants, primary metabolism, photoperiod
Citation: Ruocco M, Barrote I, Hofman JD, Pes K, Costa MM, Procaccini G, Silva J and Dattolo E (2021) Daily Regulation of Key Metabolic Pathways in Two Seagrasses Under Natural Light Conditions. Front. Ecol. Evol. 9:757187. doi: 10.3389/fevo.2021.757187
Received: 11 August 2021; Accepted: 15 November 2021;
Published: 06 December 2021.
Edited by:
Isabel Marques, University of Lisbon, PortugalReviewed by:
Susan Crosthwaite, East Malling Research, United KingdomCarlos T. Hotta, University of São Paulo, Brazil
Tomoaki Muranaka, Kagoshima University, Japan
Copyright © 2021 Ruocco, Barrote, Hofman, Pes, Costa, Procaccini, Silva and Dattolo. 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.
*Correspondence: Miriam Ruocco, bWlyaWFtLnJ1b2Njb0Bzem4uaXQ=; Emanuela Dattolo, ZW1hbnVlbGEuZGF0dG9sb0Bzem4uaXQ=