- Research
- Open access
- Published:
Acclimation of circadian rhythms in woodland strawberries (Fragaria vesca L.) to Arctic and mid-latitude photoperiods
BMC Plant Biology volume 23, Article number: 483 (2023)
Abstract
Background
Though many abiotic factors are constantly changing, the photoperiod is a predictable factor that enables plants to time many physiological responses. This timing is regulated by the circadian clock, yet little is known about how the clock adapts to the differences in photoperiod between mid-latitudes and high latitudes. The primary objective of this study was to compare how clock gene expression is modified in four woodland strawberry (Fragaria vesca L.) accessions originating from two different populations in Italy (IT1: Tenno, Italy, 45°N, IT4: Salorno, Italy, 46°N) and two in Northern Norway (NOR2: Alta, Norway, 69°N, NOR13: Indre Nordnes, Norway 69°N) when grown under simulated daylength conditions of an Arctic or mid-latitude photoperiod. The second objective was to investigate whether population origin or the difference in photoperiod influenced phytohormone accumulation.
Results
The Arctic photoperiod induced lower expression in IT4 and NOR13 for six clock genes (FvLHY, FvRVE8, FvPRR9, FvPRR7, FvPRR5, and FvLUX), in IT1 for three genes (FvLHY, FvPRR9, and FvPRR5) and in NOR2 for one gene (FvPRR9). Free-running rhythms for FvLHY in IT1 and IT4 were higher after the Arctic photoperiod, while the free-running rhythm for FvLUX in IT4 was higher after the mid-latitude photoperiod. IT1 showed significantly higher expression of FvLHY and FvPRR9 than all other accessions, as well as significantly higher expression of the circadian regulated phytohormone, abscisic acid (ABA), but low levels of salicylic acid (SA). NOR13 had significantly higher expression of FvRVE8, FvTOC1, and FvLUX than all other accessions. NOR2 had extremely low levels of auxin (IAA) and high levels of the jasmonate catabolite, hydroxyjasmonic acid (OH-JA).
Conclusions
Our study shows that circadian rhythms in Fragaria vesca are driven by both the experienced photoperiod and genetic factors, while phytohormone levels are primarily determined by specific accessions’ genetic factors rather than the experienced photoperiod.
Background
The photoperiod (daylength) is a highly predictable factor that has prompted one of the most evolutionarily conserved time-keeping mechanisms in all domains of life on Earth, a circadian clock. Mid-latitudes experience moderate seasonal variations in photoperiod compared to the equator, varying from approximately 9–15 h between winter and summer. Whereas the annual changes in photoperiod in the polar regions are extreme, ranging from nearly 24 h of photosynthetic light during summer (midnight sun) to no photosynthetic light during winter. This restricts the growth season to spring and summer due to the rapid shortening in photoperiod in autumn, which affects the timing of plant growth, flowering, biomass production, and quality [1]. The circadian clock is an important factor in physiological adaptations across regional scales [2], but how the clock balances between stability and plasticity across middle and high latitudes is unknown. The woodland strawberry, Fragaria vesca, is an herbaceous perennial with a wide geographic distribution from 37°N-70°N across the northern hemisphere. With a small, fully sequenced genome (2n = 14, ~ 200 Mb), and shared sequence identity with the cultivated octoploid strawberry (Fragaria x ananassa, 2n = 56, ~ 813 Mb), this plant is an ideal model to study clock plasticity to photoperiodic differences [3,4,5].
The mechanism of the circadian clock in plants has the same transcription-translation feedback loop (TTFL) found in almost all other organisms, comprised of a core set of genes which function as transcriptional activators [6] and/or repressors [7,8,9,10] of other clock components (Fig. 1). Although Arabidopsis has received considerable attention in the study of functional characteristics of clock components in plants, it is important to recognize that the function of individual components can vary among different species [11]. Therefore, it is imperative to evaluate the plasticity of the clock in other species to understand how the clock system can adapt to new regions under future climate scenarios.
In most plants, two morning-expressed genes, CIRCADIAN CLOCK ASSOCIATED 1 (CCA1), and/or LATE ELONGATED HYPOCOTYL (LHY), are essential components to the core clock. Not all plants contain both CCA1 and LHY genes, and F. vesca contains only LHY [12]. CCA1 and LHY proteins function as transcriptional repressors of their own transcription as well as many other clock genes [9, 13]. As CCA1 and LHY levels decline throughout the day, daytime and evening phased genes are then transcribed which further represses transcription of CCA1 and LHY [7, 9]. CCA1 and LHY are two partially redundant MYB-domain containing transcription factors that belong to the small REVEILLE subfamily including the close homolog REVEILLE 8 (RVE8). RVE8 is the main transcriptional activator of evening genes by targeting an evening element (EE) motif in their promoter regions [6].
Daytime expressed genes play roles within the regulation of both morning and evening-phased complexes of the clock. Members of the PSEUDO-RESPONSE-REGULATOR gene family (PRR9, PRR7, and PRR5) are repressed by LHY and are said to express sequentially throughout the day in Arabidopsis [13, 14], however, in F. vesca PRR9, PRR7, and PRR5 have been shown to express at nearly the same time around noon [12]. These gene function as transcriptional repressors of the morning genes CCA1, LHY, and RVE8 [15]. An evening-expressed PRR, TIMING OF CHLOROPHYLL A/B BINDING PROTEIN 1 (TOC1 or PRR1), is an indispensable gene to the core clock network, which functions to repress expression of morning-expressed genes at night [16]. An Evening Complex (EC) composed of LUX ARRYTHMO (LUX) also known as PHYTOCLOCK 1 (PCL1), and EARLY FLOWERING 3 & 4 (ELF3 & ELF4) downregulates transcription of TOC1 after dusk, allowing transcription of CCA1 and LHY to restart the next dawn [13, 17]. In F. vesca, TOC1 and LUX were the only evening genes that displayed circadian rhythms [12]. Thus, the interplay of activators and repressors keeps the clock components expressed during the right phase of the day.
The key function of the core genes involved in the plant clock also tightly regulate many output pathways. More than 1/3 of the transcriptome in plants is regulated by the circadian clock [18]. Mechanisms that enhance tolerance to UV radiation or drought can provide increased resistance to abiotic stressors, while mechanisms that enhance the timing of gene-mediated defenses or elevate levels of protective hormones can boost immune responses during periods when pathogens are expected to be active [19, 20]. For example, LHY interacts with the biosynthetic pathways of abscisic acid (ABA), which function under drought and osmotic stress [21], while regulation of jasmonic acid (JA) [22], and indole-3-acetic acid or auxin (IAA) [23] are regulated by multiple levels of the circadian network to increase tolerance to biotic stress as well as regulate plant growth and development. It has been hypothesized that jasmonate biosynthesis is promoted by high-latitude summer photoperiods which increases constitutive defense against herbivory [24].
It has been shown that plants at higher latitudes have longer circadian periods of leaf movement [25, 26], however, how these very long photoperiods affect the molecular plant circadian clock, and its downstream regulation, is much less known. Here we study the effect of differences in photoperiod between Arctic latitudes and mid-latitudes on the woodland strawberry circadian clock. To do this, two F. vesca accessions from Italy (IT1: Tenno, Italy, 45°N, IT4: Salorno, Italy, 46°N) and two from Northern Norway (NOR2: Alta, Norway, 69°N, NOR13: Indre Nordnes, Norway, 69°N), were grown under simulated daylength conditions associated with the growth season at an Arctic and mid-latitude photoperiod. First, the expression patterns of seven core clock genes were analyzed, followed by analysis of phytohormone accumulation. This study furthers the understanding of how plants adapt their circadian rhythms to new photoperiods.
Results
Plant biomass grown under Arctic and mid-latitude photoperiods
One accession (NOR13) had significantly lower biomass under the Arctic photoperiod (0.62 ± 0.15 g, 0.39 ± 0.13 g, p = 0.04). All other accessions were not significantly affected by photoperiod (IT1: 1.02 ± 0.14 g, 0.96 ± 0.23 g, IT4: 0.62 ± 0.20 g, 0.62 ± 0.23 g, NOR2: 0.64 ± 0.22 g, 0.84 ± 0.37 g).
Differences in circadian clock gene expression between Arctic and mid-latitude photoperiods and the effect on free-running rhythms
Six of the seven clock genes (all but FvTOC1) had distinct phases of expression. Here, ZT 24 is the same point of the light cycle as ZT 0, however it is not grouped with ZT 0 for statistical comparison for free-running rhythms. All post-hoc results are listed in Suppl. Table 2 for significant differences between the treatments, and Suppl. Table 3 for significant differences between the accessions.
IT1, IT4, and NOR13 had significantly higher expression of FvLHY under the mid-latitude photoperiod compared to the Arctic photoperiod at ZT 0/24 (Fig. 2). Under the mid-latitude photoperiod, IT1 had significantly higher expression than all other accessions at ZT 0, and ZT 4, while NOR2 was significantly higher than IT4 at ZT 0 and ZT 4. Under the Arctic photoperiod, IT1 had significantly higher expression than all other accessions at ZT 0 and ZT 4, while NOR2 was significantly higher than IT4 and NOR13 at ZT 4 and ZT 8. Free-running rhythms after the mid-latitude photoperiod were significantly higher for IT4 at ZT 4, and NOR2 at ZT 8, but significantly lower for IT1, IT4, and NOR13 at ZT 24. Free-running rhythms after the Arctic photoperiod were significantly higher for IT1 and IT4 at ZT 4 and 8. The free-running rhythms were also significantly higher for IT1 and IT4 at ZT 8 after the Arctic photoperiod compared to after the mid-latitude photoperiod.
IT4 and NOR13 had significantly higher expression of FvRVE8 under the mid-latitude photoperiod compared to the Arctic photoperiod at ZT 0 (Fig. 3). Under the mid-latitude photoperiod, NOR13 had significantly higher expression than IT4 and NOR2 at ZT 0. Under the Arctic photoperiod, IT1 and NOR13 had significantly higher expression than NOR2 at ZT 4. Only NOR2 had a significantly higher free-running rhythm after the mid-latitude photoperiod at ZT 8.
All accessions had significantly higher expression of FvPRR9 under the mid-latitude photoperiod compared to the Arctic photoperiod at ZT 4 and ZT 8 (Fig. 4). Under the mid-latitude photoperiod, IT1 had significantly higher expression than IT4 and NOR2 at ZT 8. Under the Arctic photoperiod, IT1 had significantly higher expression than IT4 at ZT 8. Free-running rhythms after the mid-latitude photoperiod were significantly lower for IT4, NOR2, and NOR13 at ZT 4, and for all accessions at ZT 8. Free-running rhythms after the Arctic photoperiod were significantly lower for all accessions at ZT 8.
IT4 had significantly higher expression of FvPRR7 under the mid-latitude photoperiod compared to the Arctic photoperiod at ZT 12, while NOR13 had significantly higher expression under the Arctic photoperiod at ZT 20 compared to the mid-latitude photoperiod (Fig. 5). There were no significant differences in free-running rhythms after either photoperiod.
IT1, IT4, and NOR13 had significantly higher expression of FvPRR5 under the mid-latitude photoperiod compared to the Arctic photoperiod at ZT 8, while IT4 was significantly higher under the Arctic photoperiod at ZT 16 compared to the mid-latitude photoperiod (Fig. 6). Under the mid-latitude photoperiod, NOR13 and IT4 had significantly higher expression than NOR2 and IT1 at ZT 8 and ZT 12. Under the Arctic photoperiod, NOR13 had significantly higher expression than NOR2 and IT1 at ZT 12 and ZT 16, and IT4 had significantly higher expression than NOR2 and IT1 at ZT 12. Free-running rhythms after the mid-latitude photoperiod were significantly lower for IT1, IT4, and NOR13 at ZT 8.
There were no significant differences in expression of FvTOC1 within the accessions between the two treatments, or for any free-running rhythms (Fig. 7). Under the Arctic photoperiod, NOR13 had significantly higher expression than IT1 and NOR2 at ZT 0, and all other accessions at ZT 20.
IT4 and NOR13 had significantly higher expression of FvLUX under the mid-latitude photoperiod compared to the Arctic photoperiod at ZT 12, while NOR13 was significantly higher under the Arctic photoperiod at ZT 20 compared to the mid-latitude photoperiod (Fig. 8). Under the mid-latitude photoperiod, NOR13 had significantly higher expression than IT1 and NOR2 at ZT 12, and IT4 had significantly higher expression than IT1 at ZT 12. Under the Arctic photoperiod, NOR13 had significantly higher expression than all other accessions at ZT 20. Free-running rhythms after the mid-latitude photoperiod were significantly lower for IT1 at ZT 8 and ZT 12, and NOR13 at ZT 12, but significantly higher for NOR2 at ZT 20 and IT4 at ZT 24. Free-running rhythms after the Arctic photoperiod were significantly lower for IT1, IT4, and NOR13 at ZT 12, as well as for IT1 at ZT 16. The free-running rhythms were also significantly higher for IT4 at ZT 12 after the mid-latitude photoperiod compared to after the Arctic photoperiod.
Differences in phytohormone accumulation under the Arctic photoperiod, mid-latitude photoperiod, and free-running rhythms
A partial least squares discriminant analysis (PLS-DA) of all measured phytohormones together demonstrated marked separation of IT1 from the other accessions under both photoperiods (Fig. 9). Although there was no clear diurnal trend in any of the analyzed phytohormone levels (Suppl. Figure 1), significant differences in individual phytohormone level between the accessions were found.
Abscissic acid (ABA) levels were the highest in IT1 for all treatments (p = 0.001, Fig. 10A), while auxin (IAA) levels were the lowest in NOR2 (p = 0.001, Fig. 10B). Salicylic acid (SA) levels were the lowest in IT1 for all treatments (p = 0.001) while IT4 had the highest levels (p = 0.001, Fig. 10C). The jasmonate catabolite, hydroxyjasmonic acid (OH-JA), was highest in NOR2 (p = 0.001, Fig. 10D) and accumulated in much higher concentrations than JA (Fig. 10E). JA had no significant differences in accumulation at all. JA-Ile was only significantly higher when moved to constant darkness from the Arctic photoperiod (D-Arctic) in IT4 compared to IT1 (p = 0.001) and NOR13 (p = 0.04, Fig. 10F). All significant difference values for phytohormone levels are listed in Suppl. Table 4.
Discussion
This study investigated whether seven core circadian clock genes and various phytohormones and derivatives of jasmonate in four different woodland strawberry accessions were affected by the differences in photoperiod corresponding with Arctic and mid-latitude daylength conditions that differ by 3h of photosynthetic light and a lack of darkness during the night in the Arctic (midnight sun).
Circadian rhythms have lower amplitudes under Arctic photoperiods
Clock gene expression varied depending on both experienced photoperiod and endogenous factors within specific accessions. The Arctic photoperiod induced lower expression in IT4 and NOR13 for six of the seven genes (all but FvTOC1), in IT1 for three genes (FvLHY, FvPRR9, and FvPRR5), and in NOR2 for one gene (FvPRR9). IT1 had higher expression than all other accessions of FvLHY and FvPRR9 under both photoperiods, while NOR13 had higher expression of FvRVE8, FvTOC1, and FvLUX under both photoperiods. The amplitudes of free-running rhythms for some genes were also influenced by both the previously exposed photoperiod and specific accession. IT1 and IT4 initially had higher free-running rhythms of FvLHY after the Arctic photoperiod compared to the mid-latitude photoperiod, but after the mid-latitude photoperiod, free-running expression for three of the accessions began to rise again after 24 h, whereas they didn’t after the Arctic photoperiod. IT4 had a higher free-running expression of FvLUX at ZT 12 after the mid-latitude photoperiod compared to the Arctic photoperiod, though free-running expression also appeared to lose the 24-hour rhythmicity for this accession. The attenuation of overt circadian rhythms has been observed in animals adapted to Arctic environments [27], but never in plants. Given that this phenomenon was observed in at least one gene in all accessions, including those native to middle latitudes, this suggests that reduced amplitudes of circadian rhythms under Arctic photoperiods could potentially offer some functional benefits to plants. However, it is important to note that this observation was limited to a two-week period, which could simply reflect the plasticity of the clock in response to short-term changes. Further research would be necessary to understand the long-term adjustment of the clock and its impact on physiological functions. The circadian genes here were chosen for their specific roles in the plant clock system. Both FvRVE8 and FvLHY had a strong morning phase of expression, even in the Arctic daylength which had low levels of non-photosynthetic light at night, emphasizing the dependence on the predictable onset of photosynthetic light at dawn for timing of expression. While LHY is a key repressor of all other genes in the clock [13] RVE8 is a direct transcriptional activator of evening-phased genes as well as involved in temperature compensation and light signaling [28]. The transcription profiles of these genes signifies their essential role in regulating the morning phase of the circadian clock.
FvPRR9, FvPRR7, and FvPRR5 had roughly similar phases of expression at mid-day in F. vesca, which was also found in Chen et al. [12], though, here expression was lower under the Arctic photoperiod and the duration of expression was slightly longer. In Arabidopsis, in contrast, the PRR’s have been identified to express in sequential waves after dawn starting with PRR9, followed by PRR7, PRR5, PRR3, and ending with TOC1 in the evening [29]. These mid-day-phased genes act as transcriptional repressors of CCA1 and LHY in Arabidopsis [30], and coincided with the timing of decreased LHY expression in F. vesca here, confirming similarities in clock gene functions across higher plants [31].
FvTOC1’s phase of expression occurred throughout the evening and into the morning, but had the lowest expression at ZT 4 when FvLHY had the highest expression, providing supportive evidence of reciprocal repression between TOC1 and morning genes [10]. FvTOC1’s extended phase of expression from afternoon into early morning has been reported in previous studies and aligns with TOC1’s close association with ABA to maintain cellular homeostasis during drought stress via processes such as stomatal closure, which largely occur at night [29, 32,33,34]. LUX along with the other components of the Evening Complex (EC), directly regulate PRR9 in Arabidopsis, as well as play an important role in growth in the early evening, corresponding with the evening expression for FvLUX observed here [17, 35]. The EC regulates the circadian gating of hypocotyl growth in the early evening by repressing the expression of phytochrome interacting factors PIF4 and PIF5 that help regulate plant responses to differential red (R), far-red (FR), and blue light throughout the day, and are not activated again until the following morning by interacting with LHY [17, 36]. While here the evening light was only represented by 3 µmol m–2 s−1 of white light, not FR light, all clock components still behaved in ways indicating they perceived the difference between the photosynthetic light and non-photosynthetic light representative of twilight. Thus, the components in the clock system are directly dependent on the proper transcriptional timing of all other transcription factors in the TTFL to stay in sync with the external photoperiod, even when the light only fluctuates in its photosynthetic quality. This ensures proper timing of biological functions, hormone signaling, and metabolic processes in response to multiple abiotic and biotic factors [37]. However, it is yet to be determined whether the reduction in amplitude of circadian clock gene activity observed under Arctic photoperiods can be replicated in natural field conditions and warrants further studies.
Phytohormone accumulation is dependent upon endogenous factors
Considering that the connection between the circadian clock and hormone signaling is crucial for plant growth and development, and to enhancing fitness in plants [38], the differences in experienced photoperiod or population origin were hypothesized to influence phytohormone accumulation. Phytohormone accumulation was primarily determined by genetic factors of the specific accessions rather than the experienced photoperiod or latitude of origin. Here, the accession from Tenno, Italy (IT1, 45°93’N) was identified as an accession that best discriminates between all accessions, very likely due to much higher levels of abscisic acid (ABA), but low levels of salicylic acid (SA) under all treatments, including complete darkness, while the other accession from Salorno, Italy (IT4, 46°23’N), had higher levels of SA, and the two Norwegian accession had moderate levels of SA. Since LHY interacts with the biosynthesis pathways of ABA [21], this may be the reason for IT1’s consistently higher expression of both FvLHY and ABA. Increased ABA/SA ratios were observed in drought-stressed Brassica napus (rapeseed) [39], however our plants were not exposed to any known stress factors (biotic or abiotic), thus this increase may be related to adaptations from the original environment. Additionally, there was no rhythmic accumulation of any of the phytohormones, which was unexpected as ABA [21], jasmonic acid (JA) [22], jasmonates, and salicylates [40] have been reported to accumulate rhythmically in Arabidopsis.
The accession from Alta, Norway (NOR2) had extremely low levels of indole-3-acetic acid (IAA) but high levels of the jasmonate catabolite, hydroxyjasmonic acid (OH-JA). IAA signal transduction is regulated by multiple steps of the circadian clock, where high auxin inducibility was correlated with high-amplitude circadian rhythms [23]. NOR2 had relatively low amplitudes of expression for several clock genes, which may explain NOR2’s low levels of IAA.
There was no significant difference in JA accumulation, though in a recent study, JA accumulation was significantly higher in bilberry (Vaccinium myrtillus) leaves under constant light (24 h) versus a 12 h (12 L/12 D) photoperiod [24]. It can be postulated that the lack of difference here may be due to the similarities in amount of photosynthetic light. Plants grown under the Arctic daylength received only 3 h longer photosynthetic light, and still experienced a 6 h period of non-photosynthetic light (twilight), which would be more typical of the early and late growing seasons, when herbivore-resistance induced by jasmonate levels would be less necessary than in peak summer. Additionally, the low levels of jasmonic acid-isoleucine conjugate (JA-Ile) observed here are typical of plants in the absence of stress [24]. Thus, the hypothesis that biosynthesis of defense-related phytohormones is promoted by high-latitude summer photoperiods may still hold true for F. vesca and would need to be investigated under 24 h photosynthetic light. There was only one observable difference in plant growth for NOR13, which was lower under the Arctic photoperiod, but this disparity did not cause any discernable effect on clock genes or phytohormone levels. However, we cannot exclude intrinsic determinations of phytohormones level in the different accessions, which is not unlikely as actually none of the light treatments caused significant changes within an accession. Therefore, the observed differences in endogenous phytohormone levels have a much more complex role in various physiological and biochemical processes than can be elucidated in this study.
Relation between changes in photoperiod and physiological responses for future crop production
Climate change is causing a geographic shift for many species as land at northern latitudes becomes more habitable, however the experienced photoperiod may constrain physiological responses (e.g., early flowering) that are important for crop production [41]. It has been suggested that manipulation of the circadian clock may be a strategy to overcome light-related environmental changes due to its direct control of many developmental and metabolic processes [42]. However, most studies investigating the influence of photoperiod on the diurnal rhythm of gene expression in plants typically classify differences in daylength by long days (14–16 h of light) and short days (8–10 h of light) with corresponding dark periods for a 24 h cycle [43,44,45]. This is because the differences between short-days and long-days are key determining factors for circadian regulation of floral induction [31, 46]. Yet, this does not address how plants will respond to the differences in light conditions between middle and high latitudes. Since F. vesca has a natural geographic range that spans from 37°N to 70°N, and shares a significant amount of sequence identity with the economically important cultivated strawberry, this makes it useful when considering transcriptional studies that can be related to agricultural improvement [3,4,5]. Accession-specific differences have been observed for flowering responses in F. vesca, where a high-latitude population from Alta had an obligatory vernalization requirement but no vernalization was needed for flowering of other high-latitude populations [47, 48]. However, photoperiodic responses like flowering, are complex traits with many loci that have multiple regulatory factors in response to changes in daylength [49, 50]. Therefore, while it has been shown here that F. vesca accessions from Italy and Norway can acclimate their circadian clocks to both Arctic and mid-latitudinal photoperiods, the effects of the lower amplitudes of expression on physiological functioning requires further investigation. Understanding how species adapt circadian regulation to different photoperiods can reveal important selection pressures on cellular and physiological responses that may be of crucial importance for future crop production [25, 51, 52].
Conclusions
Our study shows that the circadian clock in Fragaria vesca accessions from Italy and Northern Norway can acclimate to both Arctic and mid-latitude photoperiods. Under the Arctic photoperiod, certain clock genes exhibited reduced amplitudes compared to those under the mid-latitude photoperiod. However, the extent of these variations depended on specific accessions’ endogenous factors, indicating that the amplitude of circadian rhythms is influenced by both environmental and genetic factors. In contrast, the genetic factors of specific accessions primarily determined the accumulation of phytohormones, rather than the experienced photoperiod. These differences suggest the existence of conserved adaptations from their original environments.
Materials and methods
Plant materials and treatments
Four Fragaria vesca accessions were obtained from the University of Helsinki (Dr. Timo Hytönen’s lab). Two originated from Italy; IT1 (Tenno, Ville del Monte, Lago di Tenno): 45°93’N, 10°81’E; and IT4 (Da Salorno, Pochi, Alto, Adige, Italy): 46°23’N, 11°23’ E, and two from Northern Norway; NOR2 (Alta, Leirbukta): 69°10’ N, 23°67’ E, and NOR13 (Indre Nordnes): 69°53’ N, 20°38’ E. Clones were propagated from each of the four accessions in 1:1 [v/v] peat soil and perlite under a 12 h photoperiod (18 °C) prior to light experiments. Clones were grown for 18 days at a constant temperature of 18 °C in separate growth chambers under two daylength treatments: 18 h photosynthetic light/6 h twilight for the Arctic photoperiod, or 15 h photosynthetic light/1 h twilight/8 h dark for the mid-latitude photoperiod (Fig. 7).
Light conditions were applied using fluorescent tubes (PHILIPS Cool White MASTER TL-D Super 80 58 W/840, Eindhoven, The Netherlands) to simulate photosynthetic light and light emitting diodes (LED, PHILIPS Softone 18 W) to simulate non-photosynthetic light/twilight. Total photosynthetic photon flux density (PPFD) was adjusted in each treatment to give roughly the same amount of total photosynthetic light over the two different daylengths and measured using a quantum sensor (LI-1000, LI-COR Inc., USA). The Arctic photoperiod had a PPFD of 200 µmol m–2 s– for the photosynthetic light, the mid-latitude photoperiod had a PPFD of 243 µmol m–2 s–, and the non-photosynthetic light for both treatments had a PPFD of 3 µmol m–2 s–. Leaves from three biological replicates of each accession were collected at 4 h intervals from Zeitgeber time 0 (ZT 0) to ZT 48 (N = 114 per accession), put directly into liquid N2 and stored at -80 °C until RNA extraction. At ZT 24, half of the remaining plants from each treatment (N = 36 per accession) were transferred to total darkness and concurrently sampled to observe free-running rhythms.
RNA extraction and cDNA synthesis
Total RNA was extracted from finely ground leaf tissue following a modified protocol from Jaakola, Pirttilä [53] using the E.Z.N.A Total RNA Kit I (Omega Bio-tek Inc., USA, 2020). DNA was removed by RNase-free DNase (Qiagen, US) according to manufacturer’s instructions. The purity and concentration of RNA samples were assessed using Nanodrop 2000 spectrophotometer (ThermoFisher, USA). RNA integrity was assessed using 1.2% (w/v) agarose gel electrophoresis. cDNA was synthesized from 0.75 µg total RNA in a final reaction volume of 10 µL using Superscript IV First-Strand cDNA Synthesis Reaction protocol (Invitrogen, Carlsbad, CA, USA) according to manual instructions and stored at -20 °C until use.
Primer design
Seven circadian genes, one universally conserved protein marker, and two reference genes were selected based on previous studies in F. vesca. Five gene primer sequences (FvLHY, FvPRR9, FvPRR7, FvPRR5, and FvRVE8) and one reference gene primer sequence (FvGAPDH) were selected from Chen et al. [12]. Gene ID’s for FvTOC1, FvLUX, and FvMSI1 were searched in the F. vesca genome database (F. vesca Whole Genome v1.0 (build 8) Assembly & Annotation) at the Rosaceae genome database (GDR, https:/www.rosaceae.org) to get the genomic sequences. FvMSI1 weas subject to primer design testing in a separate unpublished scientific report. NCBI Primer-BLAST (https://www.ncbi.nlm.nih.gov/tools/primer-blast/) was used to design specific primers for remaining genes according to the standard primer parameters. Primer sequences are listed in Suppl. Table 1.
Quantitative reverse transcription PCR (qRT-PCR)
qRT-PCR assays were performed using a CFX96 Touch Real-Time PCR System (Bio-Rad, Hercules, CA, USA) to analyze RNA expression. Each 15 µL reaction mixture contained 1 µL of 10−1 diluted cDNA template, 7.5 µL of SsoFast™ EvaGreen Supermix (Bio-Rad), 3.5 µL of sterile H2O, and 1.5 µL (5 µM) of each forward and reverse primer. Amplification cycling conditions were as follows: 1 denaturation cycle of 95 °C for 2 m, followed by 40 cycles of 95 °C for 5 s, 60 °C for 30 s, ending with 1 cycle for the melt curve analysis from 65 to 95 °C in 0.5 °C increments for 5 s to verify analytical specificity. The sample-maximization method [54] was used to minimize technical (run-to-run) variation. A no-template control (NTC) was included for each gene assay to account for contamination or non-specific gene products.
Expression validation
Primer pair specificity was confirmed via single peak detection with no signals on the negative controls in the melt curve analysis. All tested genes had efficiency (E%) values ranging from 89 to 117%, with regression coefficient (R2) varying from 0.96 to 1.00.
Phytohormone measurements
Remaining ground leaf tissue was lyophilized for 4 days. Between 10 and 30 mg of dry leaf material was extracted in 1.5 ml methanol containing 60 ng D4-SA (Santa Cruz Biotechnology, USA), 60 ng D6-JA (HPC Standards GmbH, Germany), 60 ng D6-ABA (Santa Cruz Biotechnology, USA), 12 ng D6-JA-Ile (HPC Standards GmbH), and D5-indolacetic acid (D5-IAA, OlChemIm s.r.o., Olomouc, Czech Republic) as internal standards. Samples were agitated on a horizontal shaker at room temperature for 10 min. The homogenate was mixed for 30 min and centrifuged at 13,000 rpm for 20 min at 4 °C and the supernatant was collected. The homogenate was re-extracted with 500 µl methanol, mixed and centrifuged and the supernatants were pooled. The combined extracts were evaporated under reduced pressure at 30 °C and dissolved in 500 µl methanol.
Phytohormone analysis was performed by LC-MS/MS as in [55, 56] on an Agilent 1260 series HPLC system (Agilent Technologies) coupled to a tandem mass spectrometer QTRAP 6500 (SCIEX, Darmstadt, Germany). Details of the instrument parameters and response factors for quantification can be found in Suppl. Table 5.
Indolacetic acid was quantified using the same LC-MS/MS system with the same chromatographic conditions but using positive mode ionization with an ion spray voltage at 5500 eV. Multiple reaction monitoring (MRM) was used to monitor analyte parent ion → product ion fragmentations as follows: m/z 176 →130 (collision energy (CE ) 19 V; declustering potential (DP) 31 V) for indolacetic acid (IAA); m/z 181 →133 + m/z 181 →134 + m/z 181 →135 (CE 19 V; DP 31 V) for D5-indolacetic acid.
Statistical analysis
The real-time qRT-PCR expression data and phytohormone data were analyzed using Jupyter Notebook. The 2− ΔCq values were normalized to the geometric mean of the two reference genes following the ΔCq method [54] and grouped by ZT 0 - ZT 24. For each gene, the mean ± std of three biological replicates over two consecutive 24 h periods were plotted for each light treatment, and the mean ± std of three biological replicates over one 24 h period were plotted for each dark treatment. A three-way analysis of variance (ANOVA) and Tukey post-hoc test was performed for each individual gene using the statsmodels package to determine at which ZT there was a significant difference between the treatments within each accession (Suppl. Table 2), and at which ZT there was a significant difference between the accessions within each treatment (Suppl. Table 3). A two-way ANOVA and Tukey post-hoc test was performed for the phytohormone results to determine which clones had significant differences in accumulation between the treatments. Both the real-time qRT-PCR results and phytohormone results were confirmed by three biological repetitions. The scripts for the qPCR data analysis and phytohormone analysis are available from the corresponding author upon request.
Availability of data and materials
The datasets generated and analyzed during the current study are available in the Python-qPCR-analysis repository, https://zenodo.org/badge/latestdoi/359723030.
References
Mølmann JAB, Dalmannsdottir S, Hykkerud AL, Hytönen T, Samkumar A, Jaakola L. Influence of Arctic light conditions on crop production and quality. Physiol Plant. 2021;172(4):1931–40.
Oravec MW, Greenham K. The adaptive nature of the plant circadian clock in natural environments. Plant Physiol. 2022;190(2):968–80.
Folta KM, Davis TM. Strawberry genes and Genomics. CRC Crit Rev Plant Sci. 2006;25(5):399–415.
Shulaev V, Sargent DJ, Crowhurst RN, Mockler TC, Folkerts O, Delcher AL, Jaiswal P, Mockaitis K, Liston A, Mane SP, Burns P, Davis TM, Slovin JP, Bassil N, Hellens RP, Evans C, Harkins T, Kodira C, Desany B, Crasta OR, Jensen RV, Allan AC, Michael TP, Setubal JC, Celton JM, Rees DJ, Williams KP, Holt SH, Ruiz Rojas JJ, Chatterjee M, Liu B, Silva H, Meisel L, Adato A, Filichkin SA, Troggio M, Viola R, Ashman TL, Wang H, Dharmawardhana P, Elser J, Raja R, Priest HD, Bryant DW Jr, Fox SE, Givan SA, Wilhelm LJ, Naithani S, Christoffels A, Salama DY, Carter J, Lopez Girona E, Zdepski A, Wang W, Kerstetter RA, Schwab W, Korban SS, Davik J, Monfort A, Denoyes-Rothan B, Arus P, Mittler R, Flinn B, Aharoni A, Bennetzen JL, Salzberg SL, Dickerman AW, Velasco R, Borodovsky M, Veilleux RE, Folta KM. The genome of woodland strawberry (Fragaria vesca). Nat Genet. 2011;43(2):109–16.
Edger PP, Poorten TJ, VanBuren R, Hardigan MA, Colle M, McKain MR, Smith RD, Teresi SJ, Nelson ADL, Wai CM, Alger EI, Bird KA, Yocca AE, Pumplin N, Ou S, Ben-Zvi G, Brodt A, Baruch K, Swale T, Shiue L, Acharya CB, Cole GS, Mower JP, Childs KL, Jiang N, Lyons E, Freeling M, Puzey JR, Knapp SJ. Origin and evolution of the octoploid strawberry genome. Nat Genet. 2019;51(3):541–7.
Rawat R, Takahashi N, Hsu PY, Jones MA, Schwartz J, Salemi MR, Phinney BS, Harmer SL. REVEILLE8 and PSEUDO-REPONSE REGULATOR5 form a negative feedback loop within the Arabidopsis circadian clock. PLoS Genet. 2011;7(3):e1001350.
Wang Z-Y, Tobin EM. Constitutive expression of the CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) gene disrupts circadian rhythms and suppresses its own expression. Cell. 1998;93(7):1207–17.
Schaffer R, Ramsay N, Samach A, Corden S, Putterill J, Carré IA, Coupland G. The late elongated hypocotyl mutation of Arabidopsis disrupts circadian rhythms and the Photoperiodic Control of Flowering. Cell. 1998;93(7):1219–29.
Alabadí D, Oyama T, Yanovsky MJ, Harmon FG, Más P, Kay SA. Reciprocal regulation between TOC1 and LHY/CCA1 within the Arabidopsis circadian clock. Science. 2001;293(5531):880–3.
Huang W, Pérez-García P, Pokhilko A, Millar AJ, Antoshechkin I, Riechmann JL, Mas P. Mapping the Core of the Arabidopsis Circadian Clock Defines the Network Structure of the Oscillator. Science. 2012;336(6077):75–9.
Laosuntisuk K, Elorriaga E, Doherty CJ. The game of timing: circadian rhythms intersect with changing environments. Annu Rev Plant Biol. 2023;74(1):511–38.
Chen XD, Wang J, Zhao MZ, Zhao F. Characterization and expression analysis of circadian clock genes in the diploid woodland strawberry Fragaria vesca. Biol Plant. 2018;62(3):451–61.
Adams S, Manfield I, Stockley P, Carre IA. Revised morning loops of the Arabidopsis Circadian Clock based on analyses of Direct Regulatory interactions. PLoS ONE. 2015;10(12):e0143943.
Matsushika A, Makino S, Kojima M, Mizuno T. Circadian waves of expression of the APRR1/TOC1 family of Pseudo-Response regulators in Arabidopsis thaliana: insight into the Plant Circadian Clock. Plant Cell Physiol. 2000;41(9):1002–12.
Nakamichi N. The Transcriptional Network in the Arabidopsis Circadian Clock System. Genes (Basel) 2020;11(11):1284.
Huang H, Alvarez S, Bindbeutel R, Shen Z, Naldrett MJ, Evans BS, Briggs SP, Hicks LM, Kay SA, Nusinow DA. Identification of Evening Complex Associated Proteins in Arabidopsis by Affinity purification and Mass Spectrometry. Mol Cell Proteomics. 2016;15(1):201–17.
Nusinow DA, Helfer A, Hamilton EE, King JJ, Imaizumi T, Schultz TF, Farre EM, Kay SA. The ELF4-ELF3-LUX complex links the circadian clock to diurnal control of hypocotyl growth. Nature. 2011;475(7356):398–402.
Covington MF, Maloof JN, Straume M, Kay SA, Harmer SL. Global transcriptome analysis reveals circadian regulation of key pathways in plant growth and development. Genome Biol. 2008;9(8):R130.
Zhang Y, Bo C, Wang L. Novel crosstalks between Circadian Clock and Jasmonic Acid Pathway Finely coordinate the tradeoff among Plant Growth, Senescence and Defense. Int J Mol Sci. 2019;20(21):5254.
Wang W, Barnaby JY, Tada Y, Li H, Tör M, Caldelari D, Lee D-u, Fu X-D, Dong X. Timing of plant immune responses by a central circadian regulator. Nature. 2011;470(7332):110–4.
Adams S, Grundy J, Veflingstad SR, Dyer NP, Hannah MA, Ott S, Carre IA. Circadian control of abscisic acid biosynthesis and signalling pathways revealed by genome-wide analysis of LHY binding targets. New Phytol. 2018;220(3):893–907.
Thines B, Parlan EV, Fulton EC. Circadian network interactions with Jasmonate Signaling and Defense. Plants (Basel Switzerland). 2019;8(8):252.
Covington MF, Harmer SL. The circadian clock regulates auxin signaling and responses in Arabidopsis. PLoS Biol. 2007;5(8):e222.
Mithöfer A, Riemann M, Faehn CA, Mrazova A, Jaakola L. Plant defense under Arctic light conditions: can plants withstand invading pests?. Front Plant Sci. 2022;13:1051107.
Greenham K, Lou P, Puzey JR, Kumar G, Arnevik C, Farid H, Willis JH, McClung CR. Geographic Variation of Plant Circadian clock function in natural and agricultural settings. J Biol Rhythms. 2017;32(1):26–34.
Michael TP, Salomé PA, Yu HJ, Spencer TR, Sharp EL, McPeek MA, Alonso JM, Ecker JR, McClung CR. Enhanced fitness conferred by naturally occurring variation in the circadian clock. Science. 2003;302(5647):1049–53.
Bloch G, Barnes BM, Gerkema MP, Helm B. Animal activity around the clock with no overt circadian rhythms: patterns, mechanisms and adaptive value. Proc Biol Sci. 2013;280(1765):20130019.
Hsu PY, Devisetty UK, Harmer SL. Accurate timekeeping is controlled by a cycling activator in Arabidopsis. Elife. 2013;2:e00473.
Makino S, Matsushika A, Kojima M, Oda Y, Mizuno T. Light response of the circadian waves of the APRR1/TOC1 quintet: when does the quintet start singing rhythmically in Arabidopsis? Plant Cell Physiol. 2001;42(3):334–9.
Nakamichi N, Kiba T, Henriques R, Mizuno T, Chua NH, Sakakibara H. PSEUDO-RESPONSE REGULATORS 9, 7, and 5 are transcriptional repressors in the Arabidopsis circadian clock. Plant Cell. 2010;22(3):594–605.
Song YH, Ito S, Imaizumi T. Similarities in the circadian clock and photoperiodism in plants. Curr Opin Plant Biol. 2010;13(5):594–603.
Strayer C, Oyama T, Schultz TF, Raman R, Somers DE, Más P, Panda S, Kreps JA, Kay SA. Cloning of the Arabidopsis clock gene TOC1, an autoregulatory response regulator homolog. Science. 2000;289(5480):768–71.
Pokhilko A, Mas P, Millar AJ. Modelling the widespread effects of TOC1 signalling on the plant circadian clock and its outputs. BMC Syst Biol. 2013;7:23.
Legnaioli T, Cuevas J, Mas P. TOC1 functions as a molecular switch connecting the circadian clock with plant responses to drought. EMBO J. 2009;28(23):3745–57.
Helfer A, Nusinow DA, Chow BY, Gehrke AR, Bulyk ML, Kay SA. LUX ARRHYTHMO encodes a nighttime repressor of circadian gene expression in the Arabidopsis core clock. Curr Biol. 2011;21(2):126–33.
Sun Q, Wang S, Xu G, Kang X, Zhang M, Ni M. SHB1 and CCA1 interaction desensitizes light responses and enhances thermomorphogenesis. Nat Commun. 2019;10(1):3110.
Webb AAR, Seki M, Satake A, Caldana C. Continuous dynamic adjustment of the plant circadian oscillator. Nat Commun. 2019;10(1):550.
Troncoso-Ponce MA, Mas P. Newly described components and regulatory mechanisms of circadian clock function in Arabidopsis thaliana. Mol Plant. 2012;5(3):545–53.
Park S-H, Lee B-r, La H, Mamun M, Bae D-W, Kim T-H. Drought Intensity-Responsive Salicylic Acid and Abscisic Acid Crosstalk with the Sugar Signaling and metabolic pathway in Brassica napus. Plants. 2021;10(3):610.
Goodspeed D, Chehab EW, Min-Venditti A, Braam J, Covington Michael F. Covington Michael. Arabidopsis synchronizes jasmonate-mediated defense with insect circadian behavior. Proc National Acad Sci. 2012;109(12):4674–7.
Ettinger AK, Buonaiuto DM, Chamberlain CJ, Morales-Castilla I, Wolkovich EM. Spatial and temporal shifts in photoperiod with climate change. New Phytol. 2021;230(2):462–74.
Kim JA, Kim HS, Choi SH, Jang JY, Jeong MJ, Lee SI. The importance of the circadian clock in regulating Plant Metabolism. Int J Mol Sci. 2017;18(12):2680.
Weng X, Lovell JT, Schwartz SL, Cheng C, Haque T, Zhang L, Razzaque S, Juenger TE. Complex interactions between day length and diurnal patterns of gene expression drive photoperiodic responses in a perennial C4 grass. Plant Cell Environ. 2019;42(7):2165–82.
Michael TP, Mockler TC, Breton G, McEntee C, Byer A, Trout JD, Hazen SP, Shen R, Priest HD, Sullivan CM, Givan SA, Yanovsky M, Hong F, Kay SA, Chory J. Network Discovery Pipeline elucidates conserved Time-of-day–specific cis-Regulatory Modules. PLoS Genet. 2008;4(2):e14.
Dalchau N, Hubbard KE, Robertson FC, Hotta CT, Briggs HM, Stan G-B, Gonçalves JM, Webb AAR. Correct biological timing in < i > Arabidopsis requires multiple light-signaling pathways. Proc National Acad Sci. 2010;107(29):13171–6.
Suárez-López P, Wheatley K, Robson F, Onouchi H, Valverde F, Coupland G. CONSTANS Mediates Between the Circadian Clock and the Control of Flowering in Arabidopsis Nature. 2001;410(6832):1116–20.
Heide OM, Sønsteby A. Interactions of temperature and photoperiod in the control of flowering of latitudinal and altitudinal populations of wild strawberry (Fragaria vesca). Physiol Plant. 2007;130(2):280–9.
Koskela EA, Kurokura T, Toivainen T, Sønsteby A, Heide OM, Sargent DJ, Isobe S, Jaakola L, Hilmarsson H, Elomaa P, Hytönen T. Altered regulation of TERMINAL FLOWER 1 causes the unique vernalisation response in an arctic woodland strawberry accession. New Phytol. 2017;216(3):841–53.
Yanovsky MJ, Kay SA. Living by the calendar: how plants know when to flower. Nat Rev Mol Cell Biol. 2003;4(4):265–75.
Giakountis A, Cremer F, Sim S, Reymond M, Schmitt J, Coupland G. Distinct patterns of genetic variation alter flowering responses of Arabidopsis accessions to different daylengths. Plant Physiol. 2010;152(1):177–91.
Hsu PY, Harmer SL. Wheels within wheels: the plant circadian system. Trends Plant Sci. 2014;19(4):240–9.
Hut RA, Paolucci S, Dor R, Kyriacou CP, Daan S. Latitudinal clines: an evolutionary view on biological rhythms. Proc Biol Sci. 2013;280(1765):20130433.
Jaakola L, Pirttilä AM, Halonen M, Hohtola A. Isolation of high quality RNA from bilberry (Vaccinium myrtillus L.) fruit. Mol Biotechnol. 2001;19(2):201–3.
Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J. qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007;8(2):R19.
Heyer M, Reichelt M, Mithöfer A. A holistic Approach to analyze systemic Jasmonate Accumulation in Individual Leaves of Arabidopsis Rosettes upon Wounding. Front Plant Sci. 2018;9:1569.
Gallinger J, Zikeli K, Zimmermann MR, Görg LM, Mithöfer A, Reichelt M, Seemüller E, Gross J, Furch ACU. Specialized 16SrX phytoplasmas induce diverse morphological and physiological changes in their respective fruit crops. PLoS Pathog. 2021;17(3):e1009459.
Acknowledgements
We thank the Climate laboratory of the Arctic University of Norway. Amos Samkumar is acknowledged for help with sampling, and Leidulf Lund is acknowledged for growing the plants.
Funding
Open access funding provided by UiT The Arctic University of Norway (incl University Hospital of North Norway) The work was financially supported by NordPlant (NordForsk grant no. 84597) as well as the Strategic Publication Scholarship from the Arctic University of Norway.
Author information
Authors and Affiliations
Contributions
Concept by LJ, JM, and CF. TH provided the plant material. CF, LJ, and MR performed the experiments. AM contributed to the implementation of the phytohormone analysis and to the analysis of the results. CF wrote the manuscript with contribution of all authors.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Fragaria vesca were used in this study. NOR2, NOR13, IT1, and IT4 accessions were kindly provided by Dr. Timo Hytönen (Department of Agricultural Sciences, Viikki Plant Science Centre, University of Helsinki, Finland). All methods were carried out in accordance with institutional, national, and international guidelines and regulations.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Supplementary Table 1.
List of genes, primers, source, and efficiency values (E%) used in this study. Supplementary Table 2. Significant differences between treatments within each accession. Supplementary Table 3. Significant differences within accessions within each treatment. Supplementary Table 4. Significant differences in phytohormone data between accessions within treatments. Supplementary Table 5. Details of analysis of phytohormones by LC-MS/MS [HPLC 1260 (Agilent Technologies)-QTRAP6500 (SCIEX)] in negative ionisation mode. Supplementary Figure 1. Average phytohormone levels within each treatment over a 24 h period.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Faehn, C., Reichelt, M., Mithöfer, A. et al. Acclimation of circadian rhythms in woodland strawberries (Fragaria vesca L.) to Arctic and mid-latitude photoperiods. BMC Plant Biol 23, 483 (2023). https://doi.org/10.1186/s12870-023-04491-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12870-023-04491-6