Skip to main content

Auxin controls circadian flower opening and closure in the waterlily

Abstract

Background

Flowers open at sunrise and close at sunset, establishing a circadian floral movement rhythm to facilitate pollination as part of reproduction. By the coordination of endogenous factors and environmental stimuli, such as circadian clock, photoperiod, light and temperature, an appropriate floral movement rhythm has been established; however, the underlying mechanisms remain unclear.

Results

In our study, we use waterlily as a model which represents an early-diverging grade of flowering plants, and we aim to reveal the general mechanism of flower actions. We found that the intermediate segment of petal cells of waterlily are highly flexible, followed by a circadian cell expansion upon photoperiod stimuli. Auxin causes constitutively flower opening while auxin inhibitor suppresses opening event. Subsequent transcriptome profiles generated from waterlily’s intermediate segment of petals at different day-time points showed that auxin is a crucial phytohormone required for floral movement rhythm via the coordination of YUCCA-controlled auxin synthesis, GH3-mediated auxin homeostasis, PIN and ABCB-dependent auxin efflux as well as TIR/AFB-AUX/IAA- and SAUR-triggered auxin signaling. Genes involved in cell wall organization were downstream of auxin events, resulting in the output phenotypes of rapid cell expansion during flower opening and cell shrinkage at flower closure stage.

Conclusions

Collectively, our data demonstrate a central regulatory role of auxin in floral movement rhythm and provide a global understanding of flower action in waterlily, which could be a conserved feature of angiosperms.

Background

In flowering plants, many flowers open during the day and close at night, exhibiting a 24-h circadian rhythm, which is known as Linné’s floral clock [1]. This circadian movement pattern serves as an environmental adaptation characteristic for creating an opportunity for pollination [1]; in Ipomoea nil, the flower-opening time is determined by the length of dark period [2]. However, without a light-dark switch period, flowers continue to display a rhythm in opening and closing when in a complete light or dark environment using their “internal clock”, i.e. circadian clock [3]. Yon et al. found that silencing the core circadian clock genes Late Elongated Hypocotyl (LHY) and ZEITLUPE (ZTL) in Nicotiana attenuate strongly altered the floral rhythms [4]. Therefore, the internal circadian clock can recognize a 24-h rhythmic oscillation, establishing an endogenously sensing system to adapt to seasonal environmental cycles. Due to the daily changes in solar direction and irradiation, light becomes the key connection between the circadian clock and other environmental signaling pathways. So far, sophisticated light responsive networks have been intensively established in model plant developmental processes (summarized in [5]), and the output models are conceptually useful for explaining the connections between light and circadian rhythm in flowers. For example, sunflowers track the sun from east to west, reorienting during the night so that their leaves and apices face east before sunrise [6]. We term this phenomenon as phototropism, which is a typical plant response to light identified during floral rhythms [7]. Phototropism has been well characterized as being controlled by an asymmetric auxin distribution [8]. Correspondingly, application of the auxins indole acetic acid (IAA) and naphthaleneacetic acid (NAA) strongly promotes the opening of iris flowers [9, 10]. Thus, auxin is also involved in the regulation of circadian floral movement.

Auxin, one of the major plant hormones, is essential to regulate diverse processes of plant development by coordination of cell expansion and cell division, such as embryogenesis, tropic response, shoot branching, vascular development, apical dominance, flowering and fruit ripening [11, 12]. The effects of auxin on plant development have been largely attributed to the transcriptional and translational regulations governed by the tryptophan aminotransferase (TAA) and flavin monooxygenase (YUC)-dependent synthesis pathway, Transport Inhibitor Response1/Auxin-Related F-BOX (TIR1/AFB)-Aux/IAA-regulated signaling pathway and PIN-FORMED (PIN)-mediated transport pathway [13,14,15,16,17,18]. Cell wall biogenesis has been implicated as a core event downstream of auxin regulations that determines anisotropic cell expansion during plant development including floral circadian movement. In Ipomoea nil, transcripts of xyloglucan endotransglucosylase/hydrolases (XTHs) encoding the genes InXTH1-InXTH4 in petals, which are involved in cell wall modification, were closely correlated with the rate of flower opening, further controlling plastic petal growth [2]. In the ‘Mitchell’ petunia, knockdown of cell-wall-associated β-galactosidases, which determine the galactose level among cell wall polysaccharides, severely reduces the flower-opening angle due to the disruption of petal integrity [19]. Hence, a combination of transcriptional, translational and post-translational regulations leads to coordination among various processes including the circadian clock, hormone biosynthesis, light signaling, and cell wall modification to establish a proper opening-closing period in flowers. Although the phenomenon of floral movements have been broadly described during the past 20 years’ research [9, 20], the underlying regulatory mechanisms are largely elusive.

Flower is a characteristic feature of the angiosperms, and floral rhythm has long been discussed regarding the evolutionary origin of flower and subsequent diversification. Nymphaeales (also called waterlily), one of the most ancient angiosperm lineages, represents an early-diverging grade of flowering plants [21]. Nymphaeales consists of five genera with approximately 60 species [22, 23], and most waterlily species have flowers that open during the day and close at night. Thus, owing to the special phylogenetic position of Nymphaeales, evaluating flower actions in Nymphaeales is highly useful to reconstruct the molecular characteristics of angiosperms.

In this study, Nymphaeales flowers were used as a model to study the rhythm of floral movement. We observed that the intermediate segment of petal cells of the Nymphaeales flower are severely flexible and respond to the photoperiod. Exogenous application of auxin pronouncedly blocks flower closure; whereas, treatment of auxin inhibitor suppresses flower opening. Transcriptome analysis generated from the intermediate segment of petals at different day-time points provide a global understanding that auxin metabolism, signaling and transport coordinate the movement rhythm formation of flowers. In addition, cell wall modification which acts as a downstream cascade of auxin is synchronously oscillating with auxin events, resulting in circadian cell expansion to accomplish floral opening and closing. This study offers an integrated view for understanding floral circadian rhythms and provides useful information for studying the traits of flowering plants across angiosperms.

Results

The petal is responsible for sensing the circadian signal for movement

Plants keep track of the day-night photoperiod, ensuring reproductive success. To study this typical floral rhythm process and understand the mechanism by which the flower opening period is initiated, we adopted flowers from one of the most ancient angiosperm lineages, Nymphaeales (henceforth called waterlily), as the study model system. In the late summer of 2016 (photoperiod: 6:00–18:00 gradual light/18:00–6:00 gradual darkness; temperature: 25–32 °C), we tracked the flower movement of Nymphaea colorata via time-lapse photography for every 0.5 h. As shown in Additional file 1: Figure S1A, the flower remained closed at 7:30 and rapidly opened 0.5 h later (8:00); at 9:00, the opened angle of the flower reached the maximum value; and after 16:30, the flowers started closing continuously. Thus, the flowers of Nymphaea colorata showed a natural opening-closure pattern. Due to the shortage of Nymphaea colorata samples, we purchased the fresh-cut flowers of the waterlily cultivar Nymphaea nouchali as experimental materials. To test if the cut flowers exhibit a similar circadian movement pattern as natural flowers, the cut flowers were cultivated and tracked in a growth chamber at a controlled temperature at 25 °C and a 16-h light/8-h dark photoperiod (7:00–23:00 light/23:00–7:00 dark). Accompanied with the light exposure at 7:00, the flowers gradually opened. At 10:00, the opening angle rapidly increased until 14:00, when the opened angle reached the maximum of 75°. After 14:00, the flowers gradually closed, and at 23:00, the flowers were almost completely closed (Fig. 1a). Thus, with the manual manipulation of growth conditions, the cut flowers also displayed a stable opening-closure rhythm movement pattern, as shown by the key time points: 7:00 (initiation), 10:00 (rapid opening), 14:00 (maximal opening) and 18:00 (closure).

Fig. 1
figure 1

Intermediate petal determines flower opening and closure movement rhythm. a Floral opening and closure rhythm in Nymphaea nouchali cut flowers was tracked from 7:00 to 23:00. b Floral opening-closure rhythm without sepal. c Without sepal, the floral opening and closure rhythm was tracked when only 1 cm of the petal was left. White arrows highlighted the movement rhythm of the petals. d Without the sepal and stamen, the floral opening and closure rhythm was tracked when only 1 cm of the petal was left. White arrows highlighted the movement rhythm of the petals. e-f Floral opening-closure angles were measured according to the above pictures. Blue rectangular boxes below the curves highlight the maximum opening period. Scale bar, 25 mm (a-d). Error bar = SEM (e, f)

To explore the most sensitive part of flower that is responsive for period regulation, we first cut off the sepal of the flowers and tested the responsiveness between flowers with and without sepals. Without the sepals, the flowers continued to display an opening-closure pattern following period of (7:00–10:00–14:00–18:00) identical to that of the flowers with sepals. Moreover, without the sepals, the flowers exhibited a greater opening angle (maximum angle = 130°), suggesting that the sepals might act as a shield to restrict petal movement while protecting the internal flower organs (Fig. 1b, e). The sepals did not influence the floral rhythm; thus, we removed the sepals for the following experiments, owing to the convenience of study on internal flower organs.

To further evaluate the contribution of petals, the petals were divided into three sections: the top segment (> 1 cm distance to the bottom end), the intermediate segment (0.5–1 cm distance to the bottom end) and the bottom segment (< 0.5 cm to the bottom end). Petals without the top sections continued to response to the photoperiod and displayed a floral movement pattern similar to that of a complete flower (Fig. 1c, e), whereas the petals could not expand properly when only the bottom section remained (Fig. 1e, Additional file 1: Figure S1B, according to measurement standard of Additional file 1: Figure S1C). Meanwhile, the stamens exhibited an irreversible opening action (Fig. 1c, Additional file 1: Figure S1B). These observations suggest that the intermediate segment of petals is very sensitive to receive environmental signals and reacts. Furthermore, we even removed the stamens and left only 1 cm of the bottom segment of petals to verify if the stamens provide a thrust to push petal opening. However, without the stamens, the basal 1 cm segment of the petal was still able to respond to the photoperiod and display the 7:00–10:00–14:00–18:00 movement pattern (Fig. 1d, f). Thus, the physiological data indicated that the petal is the most important organ for receiving circadian signals and displays plastic curvature, and the 1 cm segment of the basal petal that connects to the sepal and the receptacle is primarily responsibility for floral movement.

Cell morphogenesis of the intermediate petal is circadian-dependent

To determine whether the intermediate segment of petal is responsive for circadian rhythm, we tracked longitudinal sections of waterlily flowers at different time points (7:00, 10:00 and 14:00 were the flower-opening points, 18:00 and 21:00 were the flower-closure points). As shown in Fig. 2a, the intermediate segments of all four petal layers responded to the photoperiod, exhibiting gradually expanded curvatures. Particularly, the intermediate segment showed a strikingly plastic change accompanied with the circadian petal movement (Fig. 2a). To better visualize the cellular change, we used Scanning Electronic Microscopy (SEM) to track the adaxial and abaxial epidermis of the intermediate petal. In both the adaxial and abaxial epidermis, long and narrow cells quickly expanded or shrank along the horizontal direction at different time points. As we observed, from 7:00 to 14:00, the cells were continuously swelling, reaching a maximum cell width of ~ 22 μm; later, from 14:00 till 21:00, the expanded cells shrank back to 19 μm (Fig. 2b-c). Obviously, the amplitude of cell length alteration in adaxial side was higher than abaxial side. Moreover, flower opening (at 10:00) elevated adaxial cell width whereas did not influence abaxial cells, implying that adaxial cells were more flexible than abaxial ones. Hence, the cellular change of intermediate segment of petal was coordinated by the alteration of cell width and cell length (Fig. 2c). We thereby calculated a ratio to compare the cell width with its individual cell length; the abaxial cells apparently exhibited a circadian oscillation of width/length ratio, as shown by the highest at 10:00 and 14:00, when the flowers rapidly opened, and the lowest at 18:00 and 21:00, when the flowers were closing (Fig. 2c). This oscillation pattern of cellular morphogenesis is tightly correlated with petal curvature. Therefore, the rapid remodeling of cells confirmed the above speculation that the intermediate petal cells receive and respond to circadian movement signals, subsequently bending via anisotropic cell expansion.

Fig. 2
figure 2

Petal intermediate cells are responsible for floral movement rhythm. a Longitudinal section of flowers showed that several petal layers were all responsive to the photoperiod. The adaxial and abaxial sides of the petals layers were marked, and white dot lines marked the intermediate and bottom segments of petals. b-c Cell morphology of the adaxial and abaxial epidermis at different time points was tracked by scanning electronic microscopy. Individual cell length and width, as well as the ratio of the cell width to the cell length, were quantified in the C chart. More than 85 cells were randomly selected and quantified in the individual samples. Scale bar, 2.5 mm (a), 25 μm (b). Error bar = SEM (c)

Flower opening-closure pattern is influenced by exogenous auxin

The plant hormone auxin is well-known to stimulate cell expansion and division via the promotion of cell wall flexibility and extensibility, which depends on the cell wall modification [24, 25]. Cell expansion at flower opening stage while shrinkage at flower closure stage are tightly correlated with the dynamic property regulated by auxin. It intrigues us to further understand if auxin is required for floral rhythm movement. To verify our speculation, we sprayed the synthetic auxin analog NAA and the auxin transport inhibitor 2,3,5-triiodobenzoic acid (TIBA) [26] on the intermediate petals before the phototracking. NAA significantly enhanced the rate and angle of flower opening (Fig. 3). Moreover, compared with the DMSO-control samples, whose closed angle was approximately 30°, the flowers could not be completely closed after NAA application, as shown by 70° closure angles at 20:00 and 80° at 23:00 (Fig. 3). It appears that NAA triggered the second rhythm of flower movement after 20:00. In contrast, TIBA suppressed flower opening. The maximum opening angle was significantly smaller than that in control, and the opening schedule was delayed for 1 h following TIBA treatment (Fig. 3). Thus, consistent with the role of auxin in the opening of the iris flower [10], our finding further unveils that auxin participates in the floral rhythm of opening and closure.

Fig. 3
figure 3

Flower opening-closure pattern is regulated by auxin. a-c Floral movement angles were measured after treatment with NAA and TIBA (10 μM NAA/TIBA was pretreated for 8 h before time-lapse photography, and then 10 μM NAA/TIBA was sprayed on the intermediate segment and bottom segments of petal for every 2 h treatment). Scale bar, 25 mm (a). Error bar = SEM (b, c). Significant analysis of the maximum opening angles among the DMSO-, NAA- or TIBA-treated samples was performed by Student’s t-test. ****, p < 0.0001

Flower opening-closure pattern is coordinated by endogenous auxin metabolism

In order to further understand the mechanism of auxin regulation and global transcription-based regulatory network occurred during the flower opening-closure processes, we performed a de-novo RNA-seq analysis based on the above described key time points of intermediate segment of petals. Samples at 6:00 (T1) were used as the control when the flowers were still under the dark period; and then, petal samples were individually collected at 7:00 (T2), 10:00 (T3), 14:00 (T4) and 18:00 (T5) (Fig. 4a). The differential transcriptome profiles were analyzed based on the standard of q-value< 0.05, fold change> 1.5, and FPKM (expected number of Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced) ≥ 1 by comparison of different couples (such as T2 vs. T1). For annotation, in total 33,333 waterlily unigenes were searched and identified against the NCBI non-redundant (Nr) database using BlastX. According to GO enrichment and KEGG pathway analysis, in the T2 vs. T1 comparison group when the flowers started to receive the light signal, most of the downregulated genes (in total 422) were involved in oxidative stress and cell wall modification, and most of the upregulated genes participated in photosynthesis or lipid biosynthesis pathways (in total 166) (Additional file 2: Figure S2). At T3, the genes related to cell wall biogenesis started to be upregulated, corresponding to rapid flower opening action. At T4 when the flower opened maximally, the expression level of photosynthesis- and lipid-related genes decreased. When the flower was closing at T5, peroxidase-related genes were activated (Additional file 2: Figure S2). Thus, cells from the intermediate area of petals undergo a series of cellular changes in conjunction with the transcriptional regulation of photosynthesis, cell wall reorganization, oxidation, and lipid biosynthesis, among others.

Fig. 4
figure 4

Global transcriptome description of flower opening-closure. a Description of RNA-seq transcriptome samples collected from the intermediate and bottom section of petals. b-i An overview of different clusters. The differentially expressed genes during the opening and closure of the flower at different time points (6:00, 7:00, 10:00, 14:00 and 18:00) were grouped into 8 clusters. The gene expression patterns of these 8 clusters is shown. Red arrows mark the clusters with upregulated genes, and blue arrows mark the clusters with downregulated genes. The global expression pattern curves were simulated in the upper right frames

To further understand the biological correlation during floral movement rhythm, we clustered all differential expressed genes into 8 clusters based on their expression pattern (Fig. 4b-i). At the onset of light stimulation at 7:00, the responsive genes of cluster 1 were immediately upregulated compared with their expression at 6:00 (Additional file 3: Figure S3A). According to a false discovery rate (FDR) < 0.05, 70 GO items of cluster 1 were found, and these items mainly belonged to photosynthesis and transmembrane transporters. Photosynthesis is a fundamental process used by plants to convert light energy to sugar and organic compounds [27]; therefore, photosynthesis I/II reaction center subunits and chlorophyll a-b binding proteins were greatly stimulated upon light stimulation (Additional file 3: Figure S3B, Additional file 4 Table S1).

Based on the physiological data, opening of flower is regulated by exogenous auxin or auxin inhibitor treatment (Fig. 3), we try to understand if endogenous auxin level fine-tunes the regulation of floral movement rhythm. Indole-3-acetic acid (IAA) is the major natural auxin in plants and primarily synthesized from the substrate, Tryptophan (Trp). Trp is firstly converted to indole-3-pyruvate (IPA) by the Tryptophan aminotransferase (TAA) enzymes, subsequently IAA is produced from IPA by flavin monooxygenase (YUC) proteins [16, 28, 29]. Cytochrome P450 enzymes such as CYP72B2, CYP79B3 and CYP83B1 are also involved in auxin synthesis [15, 16]. Using annotated genes as query input, we tracked the DEGs of auxin metabolism events which were mainly present in cluster 3 (Fig. 5a-b). In the transcriptome profile, YUC homologues of YUC8, YUC9 and YUC10 were pronouncedly upregulated at the rapid opening stage of flowers and YUC5 was activated at flower closure stage (Table 1, Fig. 5b). Similar as the oscillation of YUCs, CYP83B1 was also induced at rapid opening stage. However, TAA family showed a different expression pattern, as shown by upregulation of TAA3 but downregulation of TAA2 during flower opening (Table 1, Fig. 5b). Nevertheless, most of auxin synthesis-related genes are activated at flower opening stage, indicating that floral movement rhythm requires the participation of new synthesized auxin (Fig. 5b).

Fig. 5
figure 5

Flower opening-closure pattern is coordinated by auxin synthesis, signaling and transport. a GO items provided the transcriptome profiles of cluster 3, including vesicle trafficking, auxin signaling and cell wall modification-related genes. b-d Heatmap of individual regulators involved in auxin biosynthesis (b), auxin transport (c) and auxin signaling (d) provided the global transcriptome regulation of auxin during flower opening and closure processes (the corresponding genes were listed in Tables 1, 2 and 3 according to the labelled number “1” to “21” in each cluster) (modified from [89])

Table 1 Auxin metabolism (corresponding to Fig. 5b)

Free IAA comprises no more than 25% of the total amount of IAA, depending on the plant tissue and species. Thus, in addition to free auxin, IAA is present in a variety of modified forms, including ester-linked indole-3-acetic acid (IAA)-sugar conjugates and amide-linked IAA-amino acid conjugates [30]. Free IAA is biologically active in diverse plant developmental processes; however, most IAA conjugates are inactive, and they serve as a reservoir to maintain appropriate balance of free IAA and conjugated IAA level [30]. The conversion of IAA to IAA-amino acids is determined by the Gretchen Hagen3 (GH3) family of auxin-inducible acyl amido synthetases [31]. In contract, IAA-Leu-resistant1/ILR1-like (ILR1/ILL) encodes an amidohydrolase that cleaves IAA-amino acid conjugates to release active IAA [32]. The transcriptome profile showed that ILL5 and ILL6 were downregulated along flower opening, whereas 11 out of 16 GH3 homologues were particularly upregulated at flower opening stage (Table 1, Fig. 5b), suggesting that flower opening action triggers the frequent conversion of free IAA to conjugated IAA, which serves a storage of auxin for the further reuse. Hence, the influence of auxin synthesis and auxin homeostasis together demonstrated that flower opening triggers the transcriptional regulation of endogenous auxin metabolism, which coordinates the opening-closure rhythm formation of flower.

Auxin flow toward the cell wall determines the plastic cell expansion

Auxin is synthesized in the tissues with active cell division and expansion, and is directionally transported via a cell-to-cell manner to the entire plant body [18]. The transportation of auxin is controlled by auxin influx and efflux carriers. PINFORMEDs (PINs) are the well-known auxin efflux transporters, which are constitutively cycling and recycling between the plasma membrane (PM) and endosomal compartments to control the polar auxin flow out of the cell [17]. Expression of PIN1 homologues was relatively high at the flower opening stage while decreased when the flowers were closing (Fig. 5c, Table 2), suggesting that PIN-mediated auxin efflux is required for the rapid opening action. PP2A phosphatase is an important regulator for polar-targeting of PINs via dephosphorylation [33]. One PP2A homologue was downregulated while two PP2A homologues were upregulated during flower opening (Fig. 5c, Table 2), indicating that a fine-tune regulation of PIN phosphorylation occurs along the opening action. In addition, subcellular dynamics of PINs is coordinated by the endosomal regulators of guanine-nucleotide exchange factors for ADP-ribosylation factor GTPases (ARF-GEF) and Brefeldin A-inhibited guanine nucleotide-exchange protein (BIG) which mediate vesicle budding process [34,35,36], as well as Sorting Nexin1 (SNX1) and Vacuolar protein sorting proteins (VPS) which control vacuolar sorting process [37,38,39]. Expression of BIGs was stimulated associated with opening action, whereas SNX1 and VPSs were downregulated at maximal opening stage (Fig. 5c, Table 2), implying that the opening of flower triggers frequent endocytosis events but suppresses degradation of PINs. Altogether, compared with flower closure stage, opening action induces PIN cycling between PM and endosomal compartments, thereby resulting in more auxin flow toward the cell wall.

Table 2 Auxin transport (corresponding to Fig. 5c)

P-glycoprotein (PGP/ABCB) proteins are defined as another type of auxin efflux transporters, which interact with PINs and maintain the stability of PINs at the PM [40]. Ten ABCB homologues were found in DEGs, and 6 of them were upregulated at rapid opening stage (Fig. 5c, Table 2), suggesting that ABCBs might assist PINs to maintain high auxin flowing out of the cell during rapid flower opening. Due to the sequence similarity with PIN proteins, PIN-LIKES (PILS) proteins are defined as the putative auxin carriers which regulate intracellular auxin accumulation and homeostasis at the endoplasmic reticulum (ER), subsequently affecting nuclear auxin signaling [41]. PILS2 and PILS6 were suppressed at maximal opening of flower, while PILS3 and PILS7 were promoted along flower opening (Fig. 5c, Table 2), indicating that floral movement stimulates a feedback auxin flow among the intercellular compartments.

Among all auxin transporters, the PIN and ABCB proteins are the major auxin efflux carriers, whereas AUXIN1/LIKE-AUX1 (AUX/LAX) proteins represent the major auxin influx carriers, which are responsive for auxin uptake within the cells [42]. Interestingly, 5 out of 7 AUX1/LAX homologue was upregulated at 14:00 and 18:00 of flower closing stage (Fig. 5c, Table 2), implying that auxin uptake is required for closure action. It is apparent that flower opening promotes auxin efflux while flower closing enhanced auxin influx. The above transcriptome data supported a scenario that flower opening action pumps more auxin to the cell wall, presumably leading to the further cell wall remodeling; while flower closing action keeps auxin stay within the cell, ready for the second round of floral movement.

Auxin signaling is highly active throughout the entire floral movement rhythm

Auxin is extraordinarily multifunctional, with different cells responding very differently to changes in auxin levels, which is required for the coordination of auxin signaling at transcriptional level [43]. The core components of the nuclear auxin signaling machinery consist of three protein families: the F-box Transport Inhibitor Response 1/Auxin Signaling F-box Protein (TIR1/AFB) auxin co-receptors, the Auxin/Indole-3-Acetic Acid (Aux/IAA) transcriptional repressors, and the Auxin Response Factor (ARF) transcription factors [44]. Briefly, the small molecule auxin targets to its nuclear receptor Transport Inhibitor Response 1/Auxin Signaling F-BOX (TIR1/AFB) proteins, bringing TIR1/AFBs and transcriptional repressors Aux/IAAs together and subsequently leading to the ubiquitination and degradation of this complex. In the transcriptome profile, expression of TIR1 and AFBs was gradually upregulated associated with floral movement rhythm, and reached to the maximal level at 18:00 (Fig. 5d, Table 3), suggesting that auxin signaling is required throughout the entire processes of flower opening and closure. Aux/IAAs act as transcriptional repressors to dimerize with Auxin Response Factor (ARF), allowing the rapid changes of transcription in response to auxin [28, 43]. In DEGs of waterlily Aux/IAA homologues, 13 out of 17 IAAs maintained high transcript during the opening process, but was suppressed at closing stage (Fig. 5d, Table 3). TOPLESS (TPL) which physically interacts with IAA12 and acts as a co-repressor to regulate auxin transcription [45], was consistently oscillating with TIR1/AFBs (Fig. 5d, Table 3). Moreover, the downstream components of auxin signaling, ARFs were also generally high expressed at flower opening stage (Fig. 5d, Table 3). Thus, expression of auxin signaling regulators was generally induced associated with flower opening, suggesting that opening action requires very active auxin perception network.

Table 3 Auxin Signaling (corresponding to Fig. 5d)

Particularly, DEGs of small auxin up RNA (SAUR) family are greatly found in our transcriptome profile, which represent the largest family of early auxin response genes with elusive auxin-related functions. Recently studies have implicated the key roles of SAURs which integrate auxin and environmental signals in a wide range of plant development, including auxin-dependent tropic growth, apical hook development, leaf senescence, root growth, seed germination etc. [46]. There are 81 SAURs in Arabidopsis, 58 SAURs in rice, 18 SAURs in moss, 134 SAURs in potato and 79 SAURs in maize [46]. Due to the gene evolution and differentiation, tandem and segmental duplication events are usually happened which contributed to the expansion of the SAUR family genes. Interestingly, 25 SAUR homologues were differentially expressed during floral opening-closure process: SAUR32, SAUR36 and SAUR72 were highly expressed at the onset of flower opening (7:00), while being suppressed at the later flower movement processes; 13 out of 25 SAURs were particularly induced associated with rapid opening process (10:00); and the rest of SAURs were stimulated at the closing stage (18:00) (Fig. 5d, Table 3). Thus, different SAUR members displayed distinct responsiveness to flower opening or closing actions, supporting the underlying roles of SAURs involved in auxin-mediated anisotropic cell expansion during floral movement.

In summary, through the regulation of auxin metabolism, auxin transportation and auxin signaling, auxin serves as a key phytohormone to determine the circadian movement rhythm of flowers.

Cell wall remodeling is downstream of auxin events and controlled by auxin oscillation during floral movement

In addition to auxin regulations, plenty of vesicle trafficking-related genes were found in cluster 3, which maintained high expression levels from 7:00 till 10:00 (Figs. 4d, 5a, Additional file 5: Figure S4A), including multiple COPI vesicle coat genes such as Coatomer Subunits, Conserved oligomeric Golgi complex subunits (COG) and Protein transport Sec genes (Additional file 4: Table S2), suggesting that vesicle trafficking is the earliest responding process during the rapid opening stage. Vesicle trafficking is a dynamics and rapid movement process, including endocytosis and exocytosis, packaging the signal molecules from the synthesis locations, such as endoplasmic reticulum (ER) to Golgi body; the molecules are finally released to their locations either on the plasma membrane, secreted to the extracellular space or targeted to the vacuoles for degradation [47]. Particularly, vesicle flow is crucial for cellular morphogenesis as it delivers membrane cargoes and materials for cell wall formation. Thus, promotion of vesicle trafficking might activate reorganization of the cell wall components [48].

Plant cell walls are rigid, flexible and extensible layers composed of polysaccharides, proteins, and lignin. Among the wall polysaccharides, cellulose, made of β-1,4-linked glucan microfibril, is the major wall component synthesized by cellulose synthase (CesA) complexes [49]. In the transcriptome profiles, multiple cell wall synthesis- or modification-related genes were differently responding to floral movement at distinct stages. For example, CESAs were gradually upregulated at 7:00 till 18:00, suggesting that cell wall is continuously synthesized along the entire floral movement process (Fig. 6c, Table 4). The synthesized cell wall components can be loosened either by cell wall-loosening proteins Expansin (EXP) or xyloglucan endotransglucosylase (XTH/XET). EXPs promote cell wall acidification to increase wall extensibility, and XETs break the xyloglucan cross-links between cellulose microfibrils to loosen cell wall [50, 51]. Several EXPs were highly expressed from 6:00 to 10:00 of the early stage of flower opening, whereas XETs showed upregulation from 10:00 till 18:00 (Fig. 6c, Table 4). Regarding the previous study showed that auxin triggers rapid cell wall acidification [52], the complementary responsiveness of EXPs and XETs to floral movement rhythm in our study also implies that EXPs-mediated wall loosening is important for early flower opening, probably in conjunction of auxin-induced rapid wall acidification; while the later stage of floral movement requires XETs-dependent cellulose microfibrils rearrangement.

Fig. 6
figure 6

Cell wall remodeling is controlled by auxin during floral movement. a-b GO items provided the transcriptome profiles of cluster 5 and 8, including cell wall organization and modification, transmembrane receptor and kinase-mediated signaling cascades. c Heatmap of components involved in cell wall organization and modification provided the global transcriptome regulation of cell wall remodeling event, in conjunction of auxin regulation and floral movement rhythm (the corresponding genes were listed in Table 4 according to the labelled number “22” to “27” in each cluster)

Table 4 Cell wall synthesis (corresponding to Fig. 6C)

In addition, pectin is another major component of the plant cell wall, and it contains a highly heterogeneous group of polymers including homogalacturonans and rhamnogalacturonans I and II [53]. Pectinesterase (PME) catalyzes the de-esterification of pectin into pectate and methanol to alter the integrity of the cell wall [54]. Consistent with XETs expression, PMEs were additionally activated from 10:00 to 18:00 (Fig. 6c, Table 4), implying that a dramatic degradation and re-establishment of the cell wall occurred during the floral movement. Apparently, activation of these genes involved in cell wall remodeling was sustained throughout the entire process of floral movement, suggesting the importance of cell wall organization during flower opening and closure (Fig. 6c, Table 4).

Strikingly, a large number of plasma membrane-associated kinases, such as Serine/threonine-protein kinase, Brassinosteroid-Signaling Kinase 3 (BSK3) and NIMA (never in mitosis gene a)-related kinase 5 (NEK5) [55, 56]; Leucine-rich repeat receptor (LRR)-like serine/threonine-protein kinase, Receptor-like protein kinase 2 (RPK2) and Barely Any Meristem 1 (BAM1) [57, 58]; Receptor-like protein kinase FERONIA and transmembrane kinase (TMK) [59, 60], were found in cluster 8 (Additional file 5: Figure S4B, Additional file 4: Table S3). Most of them belong to plasma membrane-localized receptor-like kinases (RLKs). Ligands-activated RLKs are involved in a wide range of plant growth processes, such as plant immunity, male-female interaction and cell wall modification, which perceive and transmit extracellular signals into the cytosolic, leading to signaling transduction by phosphorylation [61,62,63]. For example, THESEUS1 (THE1) is a cell wall integrity sensor to mediate the response of plant cells to the perturbation of cellulose synthesis [64]. FERONIA (FER) elicits calcium transients to maintain cell wall integrity during salt stress [65]. In the transcriptome profile, THE1 and FER were upregulated from 10:00 till 18:00 (Fig. 6c, Additional file 4: Table S3), which raised a possibility that RLKs might be the signal sensors of cell wall to coordinate the floral movement rhythm.

It has been well described that auxin induces cell elongation and expansion through its effect on modification of the cell wall [24]. The above transcriptome data further confirm that auxin signaling and cell wall remodeling were synchronously oscillating, leading to the output phenotype of cell expansion at flower opening stage and cell shrinkage at flower closure stage. Thus, we would speculate a scenario of floral movement rhythm: I. at flower opening stage when cells anisotropically expand, auxin is highly synthesized and transported to the cell wall, subsequently leading to the dramatic loosening and reorganization of cell wall; II. At flower closure stage when cell expansion are restricted, auxin synthesis is relatively delayed with low auxin flow to the cell wall (Fig. 6c). Thus, our finding further unveils that auxin-mediated cell wall remodeling is essential for the movement rhythm formation of flower opening-closure.

Discussion

For the purpose of pollination, flowers open for a considerable period, and then abscise immediately after pollination [66]. Several factors are involved in this process, including transcriptional factors, hormonal regulation, cell wall modification and circadian rhythm control [9, 20]. In our study, one of the most ancient flowering plants, Nymphaeales, is used to explore the mechanism of floral rhythm. This work globally described auxin-dependent transcriptional signatures of flower opening-closure movement in Nymphaeales, which could be conserved in most of flowering plants.

Van Doorn et al. have summarized the major mechanisms during flower opening and closure [9, 20]. They concluded that turgor pressure and water transport are the major restrictive factors to limit fast movement. During the flower opening of the Asiatic Lily, the most influential factor is the change in angle between the midribs and the pedicel, where an increasing curvature of the midribs pulls the tepals apart [67]. Additionally, Liang and Mahadevana indicated that the edges of lily tepals are wrinkled as the flower opens [68]. Plant structures overcome these limits of turgor pressure, such as sepals of the waterlily flower, resulting in a very rapid cell expansion. Our research revealed that the intermediate petal is the major tissue responding to endogenous change and environmental stimuli. The realistic petal curvature determines the efficiency of insect-propagated pollination. Thus, plastic intermediate cells of petal offer a good platform to study the global signaling events during the circadian floral movement.

Consistent with the previous studies on light signaling which plays fundamental roles in diverse aspects of plant development, including floral movement rhythm, the crucial light signaling components were rapidly responding at the onset of flower opening (Additional file 3: Figure S3C). Using Arabidopsis homologues as input, the transcript oscillations of waterlily homologues of red/far-red light receptors Phytochromes (PHYs) [69], blue light receptor Cryptochromes (CRYs) [70] and Phototropins (PHOTs) [71], Ultraviolet-B receptor UV-B resistance 8 (UVR8) [72], and their downstream effectors including Phytochrome-Interacting Factor 3 (PIF3), Constitutive Photomorphogenic 1 (COP1), Suppressor Of Phytochrome A1 (SPA1), Elongated Hypocotyl 5 (HY5) [73] were depicted in a light signaling network (Additional file 3: Figure S3C, Additional file 4 Table S4). Moreover, the transcript level of the core circadian clock regulators, downstream of light signaling, such as Pseudo-Response Regulators (APRR) [74], Late Elongated Hypocotyl (LHY) [75] and Casein kinase (CKL) [76], were additionally triggered associated with flower opening (Additional file 3: Figure S3C, Additional file 4: Table S4). Therefore, the consistent upregulation of photosynthesis, light signaling and circadian clock events support the fundamental role of light at the beginning of flower opening.

Beside of the regulations of light signaling, auxin and cell wall, meanwhile other signaling pathways are significantly influenced associated with floral movement, as seen by downregulations of transcription factors and ribonucleoproteins in cluster 2, 4 and 6 (Additional file 6: Figure S5, Additional file 4: Tables S5-S7). In model plants, increasing numbers of TFs are shown to participate in flower development. For example, the expression of ERF TFs is highly correlated with flower senescence, because increased ethylene is produced during floral organ senescence [77, 78]. WRKY12 and WRKY13 oppositely regulate flowering under a short-day condition by controlling the plant hormone gibberellin [61]. Interestingly, RhNAC100 transcriptional factor was identified in rose petal, which is rapidly induced by ethylene in the petals. Silencing of RhNAC100 significantly increased the petal size and promoted cell expansion in the petal abaxial sub-epidermis [79]; this result is consistent with our finding that NAC TFs were downregulated during flower opening and petal cell expansion. Therefore, floral opening was associated with lower transcriptional cascades, suggesting that these transcriptional factors might act as negative regulators of cell expansion, which deserve further investigation.

Conclusion

Similar to all touch responses, which can turn plants away from aggressors or enable flowers to move toward sunrise, morphogenesis during the flower opening-closure process is also influenced by a combination of endogenous and exogenous perturbations [80]. In our study, transcriptome profiles of the intermediate cells of petal were unearthed in the waterlily, which confirm that the signaling cascades of auxin metabolism, signaling and transport play center roles in flower opening and closure, and cell wall remodeling event regulated by auxin realizes circadian cell expansion and shrinkage associated with flower movement rhythm. Additionally, photosynthesis, light signaling and transcriptional regulations are individually incorporated into different stages of floral movement. Our study provided a preliminary understanding of the opening-closure processes of the waterlily flower by dividing them into three stages (Fig. 7).

Fig. 7
figure 7

Preliminary model during flower opening and closure processes. During 7:00–23:00 light/23:00–6:00 dark photoperiod, the flowers gradually open from 7:00 to 14:00, and then followed by the closure process from 14:00 to 23:00. We divide the process into three stages: I. Rapid-response stage (7:00–10:00), involved in the upregulation of auxin signaling, cell wall modification and photoreaction events; II. Continuously blooming stage (10:00–14:00), involved in activation of auxin synthesis, auxin signaling, auxin efflux and vesicle trafficking; and suppression of transcriptional and translational events; III. Late-response stage (14:00–18:00), involved in the stimulation of processes including auxin signaling, auxin influx, cell wall synthesis and remodeling; and suppression of transcription. The upregulated pathways are marked red, and the downregulated pathways are marked blue

Stage I (7:00–10:00): rapid-response stage, when the flowers are initially triggered for opening. At this stage, photosynthesis and light perception are activated, while gene transcription is blocked. Photosynthesis is the first reactive event once the flowers are exposed to light, which utilizes light energy to synthesize carbohydrate molecules (sucrose) for consumption or storage. Photoreception-mediated extension of the cell volume of the petal intermediate segment may be initiated by an unloading of sucrose from the phloem, possibly resulting in the lower apoplastic water potential [80], subsequently causing a change in cell turgor pressure. Active auxin perception system creates a complex but flexible environment, being able to re-arrange cell wall components (Fig. 7).

Stage II (10:00–14:00): flower-opening stage, when the flowers are quickly blooming until the maximal angle. At this stage, vesicle trafficking, auxin synthesis, auxin signaling and efflux, and cell wall remodeling processes occupy dominant roles during flowers blooming, accompanied with downregulation of transcription. Frequent vesicles trafficking events enhance the endocytosis and exocytosis of PM-located cargoes, such as auxin transporters PIN proteins, facilitating active auxin flow during organogenesis [81]. Apparently, auxin-related genes, such as GH3 family, AUX/IAAs, ARFs, SAURs are enriched in cluster 3 which particularly responded to rapid flower opening. YUC-dependent auxin synthesis, GH3-regulated auxin homeostasis, TIR1/AFBs and SAUR-mediated auxin signaling, as well as PIN- and PGP-regulated auxin efflux are primarily induced during stage II. A high level of auxin is synthesized and then being pumped to the cell wall; consequently, cell walls are acidified by high auxin and loosened by expansins [82].

Stage III (14:00–18:00): late response stage, when flowers start to be closed. The waterlily represents a typical flowering plant with a cross-pollination trait. The opened flower waits for insects that bring pollens from other flowers, crawling down towards its stigma for further pollination [80]. Apparently, the completely opened flowers with expanded stamens facilitate insects in touching the bottom of the stigma. Therefore, the fertilization event occurs usually after the maximum opening (after 14:00). In our transcriptome profiles, several RLKs are activated during stage III, possibly leading to fertilization. On the other hand, RLKs might function as cell-wall signaling receptors, in conjunction with cell-wall modification enzymes, to coordinate cell wall reorganization during flower closure. Auxin influx and cell wall modification events are still active in stage III, which confirms the crucial contributions of auxin and cell wall throughout the entire processes of flower opening and closure.

Methods

Plant materials and growth condition

Nymphaea colorata (Peter, Gustav Albert, 1928, Abhandlungen der Königlichen Gesellschaft der Wissenschaften, Göttingen. Mathematisch-Physikalische Klasse, n. s. 13(2): 58, 68, f. 10. 1928, website: http://www.tropicos.org/Name/22600193 or https://en.wikipedia.org/wiki/Nymphaea_colorata, preserved in Botanischen Garten, Berlin.) plants were cultivated under a natural environment, and the flower opening-closure pictures of Nymphaea colorata were captured in September of 2016 (Additional file 1: Figure S1A). For all the remaining experimental tests, cut flowers (with long stalks) of the waterlily cultivar Nymphaea nouchali (Burman, Nicolaas Laurens (Nicolaus Laurent), 1768, Flora Indica. .. nec non Prodromus Florae Capensis 120. 1768, website: http://www.ipni.org/ipni/idPlantNameSearch.do?id=135096-3&output_format=lsid-metadata&show_history=true, or https://en.wikipedia.org/wiki/Nymphaea_nouchali, preserved in Gandhi Park, East Fort, Thiruvananthapuram, India), which were commercially ordered from the market and hydroponically grown in a chamber controlled at 25 °C with a 16-h light/8-h dark photoperiod (7:00 am-23:00 pm light/23:00 pm-7:00 am dark), were used.

Growth condition and chemical treatment

NAA (Biosharp) and TIBA (Macklin) were dissolved in dimethyl sulfoxide (DMSO). For the chemical treatment, 10 μM NAA and 10 μM TIBA were pre-sprayed on the intermediate and bottom section of the petals for 8 h before phototracking, and then the floral movement was tracked for the following 24 h by capturing images using an automatic camera.

Time-lapse photography and measurement of flowering angles

The flower opening-closure patterns were tracked using time-lapse photography every 1 h for a total of 24 h. Deviated angles of petal opening and closure were measured using ImageJ (Additional file 1: Figure S1C). In each experiment, at least 10 flowers were analyzed, or at least 85 cells were measured. The experiments were repeated three times. Statistical analyses were performed using unpaired two-tailed t-test, where ns, *, **, *** and **** correspond to no significant difference, p-value < 0.05, < 0.01, < 0.001, < 0.0001, respectively.

RNA-seq, expression annotation, GO and KEGG pathway enrichment analyses

The sepals of each flowers were removed, and the extracted samples were collected from the intermediate and bottom segments of petal (≤1 cm distance to the bottom end) at different times of 6:00, 7:00, 10:00, 14:00 and 18:00 (Fig. 3a). Each sample was collected from at least three flowers, and the samples at each time point were collected in triplicated and ground into fine powder with liquid nitrogen for RNA extraction. The total RNA of the waterlily flowers was extracted using the RNAprep Pure Plant Kit (Tiangen). RNA-Seq libraries were sequenced on an Illumina HiSeq X platform. After demultiplexing, the adapter sequences were trimmed, and poly-N and low-quality reads were removed to obtain the high-quality, clean reads. The Next Generation Sequencing (NGS) reads for the samples were listed in Additional file 4: Table S8. Trinity V2.06 was used to assemble the Paired-end clean reads and the assembled results indicated the contig Nx statistics (eg. the contig N50 value) (Additional file 4: Table S9). TransDecoder V2.0.1 was applied to predict the protein-coding regions of the transcripts. The NCBI/nr, SWISS Prot, and Phytozome V12 databases were searched to annotate the peptide functions [83]. To quantify the expression levels of genes, Trinity reference was annotated, and paired-end reads were mapped to the assembled transcripts and assigned to genes using RSEM V1.2.20. Differential expression analysis was performed using R package DEseq2 V1.18.1. The FPKM (expected number of Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced) of each gene was calculated based on the length of the gene and the reads counts mapped to this gene. Genes with a fold change > 1.5 and a q value< 0.05 were assigned as differentially expressed genes (DEGs). The DEGs were mapped to GO terms in the GO database (http://www.geneontology.org/) in order to calculate the gene numbers for every term. KEGG (http://www.genome.jp/kegg/) is used to perform the pathway analysis of the DEGs. Statistical enrichment of the DEGs for the GO terms implemented by the agriGO v2.0. Significantly enriched GO terms and KEGG pathways (corrected p < 0.05) were identified based on a hypergeometric test [84, 85].

Cluster analysis was performed using the MeV4.8.1 software. Firstly, data were normalized by genes/rows. Secondly, significant DEGs were selected according to significance analysis for microarrays Delta = 0.26. The k-means method was used with the following settings: Pearson correlation, number of clusters = 8, and maximum iterations = 50 [86].

References

  1. Frund J, Dormann CF, Tscharntke T. Linne’s floral clock is slow without pollinators--flower closure and plant-pollinator interaction webs. Ecol Lett. 2011;14(9):896–904.

    Article  PubMed  Google Scholar 

  2. Shinozaki Y, Tanaka R, Ono H, Ogiwara I, Kanekatsu M, van Doorn WG, et al. Length of the dark period affects flower opening and the expression of circadian-clock associated genes as well as xyloglucan endotransglucosylase/hydrolase genes in petals of morning glory (Ipomoea nil). Plant Cell Rep. 2014;33(7):1121–31.

    PubMed  CAS  Google Scholar 

  3. Mas P. Circadian clock signaling in Arabidopsis thaliana: from gene expression to physiology and development. Int J Dev Biol. 2005;49(5–6):491–500.

    Article  PubMed  CAS  Google Scholar 

  4. Yon F, Joo Y, Cortes Llorca L, Rothe E, Baldwin IT, Kim SG. Silencing Nicotiana attenuata LHY and ZTL alters circadian rhythms in flowers. New Phytol. 2016;209(3):1058–66.

    Article  PubMed  CAS  Google Scholar 

  5. Jiao Y, Lau OS, Deng XW. Light-regulated transcriptional networks in higher plants. Nat Rev Genet. 2007;8(3):217–30.

    Article  PubMed  CAS  Google Scholar 

  6. Atamian HS, Creux NM, Brown EA, Garner AG, Blackman BK, Harmer SL. Circadian regulation of sunflower heliotropism, floral orientation, and pollinator visits. Science. 2016;353(6299):587–90.

    Article  PubMed  CAS  Google Scholar 

  7. Harmer SL, Brooks CJ. Growth-mediated plant movements: hidden in plain sight. Curr Opin Plant Biol. 2018;41:89–94.

    Article  PubMed  CAS  Google Scholar 

  8. Goyal A, Szarzynska B, Fankhauser C. Phototropism: at the crossroads of light-signaling pathways. Trends Plant Sci. 2013;18(7):393–401.

    Article  PubMed  CAS  Google Scholar 

  9. van Doorn WG, Van Meeteren U. Flower opening and closure: a review. J Exp Bot. 2003;54(389):1801–12.

    Article  PubMed  CAS  Google Scholar 

  10. van Doorn WG, Dole I, Celikel FG, Harkema H. Opening of Iris flowers is regulated by endogenous auxins. J Plant Physiol. 2013;170(2):161–4.

    Article  PubMed  CAS  Google Scholar 

  11. Mockaitis K, Estelle M. Auxin receptors and plant development: a new signaling paradigm. Annu Rev Cell Dev Biol. 2008;24:55–80.

    Article  PubMed  CAS  Google Scholar 

  12. Chapman EJ, Estelle M. Mechanism of auxin-regulated gene expression in plants. Annu Rev Genet. 2009;43:265–85.

    Article  PubMed  CAS  Google Scholar 

  13. Dharmasiri N, Dharmasiri S, Estelle M. The F-box protein TIR1 is an auxin receptor. Nature. 2005;435(7041):441–5.

    Article  PubMed  CAS  Google Scholar 

  14. Kepinski S, Leyser O. The Arabidopsis F-box protein TIR1 is an auxin receptor. Nature. 2005;435(7041):446–51.

    Article  PubMed  CAS  Google Scholar 

  15. Zhao Y. Auxin biosynthesis and its role in plant development. Annu Rev Plant Biol. 2010;61:49–64.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Zhao Y. Auxin biosynthesis: a simple two-step pathway converts tryptophan to indole-3-acetic acid in plants. Mol Plant. 2012;5(2):334–8.

    Article  PubMed  CAS  Google Scholar 

  17. Adamowski M, Friml J. PIN-dependent auxin transport: action, regulation, and evolution. Plant Cell. 2015;27(1):20–32.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Habets ME, Offringa R. PIN-driven polar auxin transport in plant developmental plasticity: a key target for environmental and endogenous signals. New Phytol. 2014;203(2):362–77.

    Article  PubMed  CAS  Google Scholar 

  19. O'Donoghue EM, Somerfield SD, Deroles SC, Sutherland PW, Hallett IC, Erridge ZA, et al. Simultaneous knock-down of six beta-galactosidase genes in petunia petals prevents loss of pectic galactan but decreases petal strength. Plant Physiol Biochem. 2017;113:208–21.

    Article  PubMed  CAS  Google Scholar 

  20. van Doorn WG, Kamdee C. Flower opening and closure: an update. J Exp Bot. 2014;65(20):5749–57.

    Article  PubMed  CAS  Google Scholar 

  21. Donoghue MJ. Phylogenies and the analysis of evolutionary sequences, with examples from seed plants. Evolution. 1989;43(6):1137–56.

    Article  PubMed  Google Scholar 

  22. Friis EM, Pedersen KR, Crane PR. Fossil evidence of water lilies (Nymphaeales) in the early cretaceous. Nature. 2001;410(6826):357–60.

    Article  PubMed  CAS  Google Scholar 

  23. Iles WJ, Lee C, Sokoloff DD, Remizowa MV, Yadav SR, Barrett MD, et al. Reconstructing the age and historical biogeography of the ancient flowering-plant family Hydatellaceae (Nymphaeales). BMC Evol Biol. 2014;14:102.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Lehman TA, Smertenko A, Sanguinet KA. Auxin, microtubules, and vesicle trafficking: conspirators behind the cell wall. J Exp Bot. 2017;68(13):3321–9.

    Article  PubMed  CAS  Google Scholar 

  25. Majda M, Robert S. The Role of Auxin in Cell Wall Expansion. Int J Mol Sci. 2018;19(4).

  26. Amijima M, Iwata Y, Koizumi N, Mishiba K. The polar auxin transport inhibitor TIBA inhibits endoreduplication in dark grown spinach hypocotyls. Plant Sci. 2014;225:45–51.

    Article  PubMed  CAS  Google Scholar 

  27. Olson JM. Photosynthesis in the Archean era. Photosynth Res. 2006;88(2):109–17.

    Article  PubMed  CAS  Google Scholar 

  28. Woodward AW, Bartel B. Auxin: regulation, action, and interaction. Ann Bot. 2005;95(5):707–35.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Yue J, Hu X, Huang J. Origin of plant auxin biosynthesis. Trends Plant Sci. 2014;19(12):764–70.

    Article  PubMed  CAS  Google Scholar 

  30. Ludwig-Muller J. Auxin conjugates: their role for plant development and in the evolution of land plants. J Exp Bot. 2011;62(6):1757–73.

    Article  PubMed  CAS  Google Scholar 

  31. Staswick PE, Serban B, Rowe M, Tiryaki I, Maldonado MT, Maldonado MC, et al. Characterization of an Arabidopsis enzyme family that conjugates amino acids to indole-3-acetic acid. Plant Cell. 2005;17(2):616–27.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. Bartel B, Fink GR. ILR1, an amidohydrolase that releases active indole-3-acetic acid from conjugates. Science. 1995;268(5218):1745–8.

    Article  PubMed  CAS  Google Scholar 

  33. Michniewicz M, Zago MK, Abas L, Weijers D, Schweighofer A, Meskiene I, et al. Antagonistic regulation of PIN phosphorylation by PP2A and PINOID directs auxin flux. Cell. 2007;130(6):1044–56.

    Article  PubMed  CAS  Google Scholar 

  34. Teh OK, Moore I. An ARF-GEF acting at the Golgi and in selective endocytosis in polarized plant cells. Nature. 2007;448(7152):493–6.

    Article  PubMed  CAS  Google Scholar 

  35. Naramoto S, Otegui MS, Kutsuna N, de Rycke R, Dainobu T, Karampelias M, et al. Insights into the localization and function of the membrane trafficking regulator GNOM ARF-GEF at the Golgi apparatus in Arabidopsis. Plant Cell. 2014;26(7):3062–76.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Kitakura S, Adamowski M, Matsuura Y, Santuari L, Kouno H, Arima K, et al. BEN3/BIG2 ARF GEF is involved in Brefeldin A-sensitive trafficking at the trans-Golgi network/early endosome in Arabidopsis thaliana. Plant Cell Physiol. 2017;58(10):1801–11.

    Article  PubMed  CAS  Google Scholar 

  37. Hanzawa T, Shibasaki K, Numata T, Kawamura Y, Gaude T, Rahman A. Cellular auxin homeostasis under high temperature is regulated through a sorting NEXIN1-dependent endosomal trafficking pathway. Plant Cell. 2013;25(9):3424–33.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Lee HN, Zarza X, Kim JH, Yoon MJ, Kim SH, Lee JH, et al. Vacuolar trafficking protein VPS38 is dispensable for autophagy. Plant Physiol. 2018;176(2):1559–72.

    Article  PubMed  CAS  Google Scholar 

  39. Wang HJ, Hsu YW, Guo CL, Jane WN, Wang H, Jiang L, et al. VPS36-dependent multivesicular bodies are critical for Plasmamembrane protein turnover and vacuolar biogenesis. Plant Physiol. 2017;173(1):566–81.

    Article  PubMed  CAS  Google Scholar 

  40. Blakeslee JJ, Bandyopadhyay A, Lee OR, Mravec J, Titapiwatanakun B, Sauer M, et al. Interactions among PIN-FORMED and P-glycoprotein auxin transporters in Arabidopsis. Plant Cell. 2007;19(1):131–47.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Barbez E, Kubes M, Rolcik J, Beziat C, Pencik A, Wang B, et al. A novel putative auxin carrier family regulates intracellular auxin homeostasis in plants. Nature. 2012;485(7396):119–22.

    Article  PubMed  CAS  Google Scholar 

  42. Swarup R, Peret B. AUX/LAX family of auxin influx carriers-an overview. Front Plant Sci. 2012;3:225.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Leyser O. Auxin Signaling. Plant Physiol. 2018;176(1):465–79.

    Article  PubMed  CAS  Google Scholar 

  44. Lavy M, Estelle M. Mechanisms of auxin signaling. Development. 2016;143(18):3226–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  45. Szemenyei H, Hannon M, Long JA. TOPLESS mediates auxin-dependent transcriptional repression during Arabidopsis embryogenesis. Science. 2008;319(5868):1384–6.

    Article  PubMed  CAS  Google Scholar 

  46. Ren H, Gray WM. SAUR proteins as effectors of hormonal and environmental signals in plant growth. Mol Plant. 2015;8(8):1153–64.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Chen X, Irani NG, Friml J. Clathrin-mediated endocytosis: the gateway into plant cells. Curr Opin Plant Biol. 2011;14(6):674–82.

    Article  PubMed  CAS  Google Scholar 

  48. Lunn D, Phan TD, Tucker GA, Lycett GW. Cell wall composition of tomato fruit changes during development and inhibition of vesicle trafficking is associated with reduced pectin levels and reduced softening. Plant Physiol Biochem. 2013;66:91–7.

    Article  PubMed  CAS  Google Scholar 

  49. Endler A, Persson S. Cellulose synthases and synthesis in Arabidopsis. Mol Plant. 2011;4(2):199–211.

    Article  PubMed  CAS  Google Scholar 

  50. Park YB, Cosgrove DJ. Xyloglucan and its interactions with other components of the growing cell wall. Plant Cell Physiol. 2015;56(2):180–94.

    Article  PubMed  CAS  Google Scholar 

  51. Cosgrove DJ. Plant expansins: diversity and interactions with plant cell walls. Curr Opin Plant Biol. 2015;25:162–72.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  52. Fendrych M, Leung J, Friml J. TIR1/AFB-aux/IAA auxin perception mediates rapid cell wall acidification and growth of Arabidopsis hypocotyls. elife. 2016;5:e19048.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Micheli F. Pectin methylesterases: cell wall enzymes with important roles in plant physiology. Trends Plant Sci. 2001;6(9):414–9.

    Article  PubMed  CAS  Google Scholar 

  54. Jolie RP, Duvetter T, Van Loey AM, Hendrickx ME. Pectin methylesterase and its proteinaceous inhibitor: a review. Carbohydr Res. 2010;345(18):2583–95.

    Article  PubMed  CAS  Google Scholar 

  55. Tang W, Kim TW, Oses-Prieto JA, Sun Y, Deng Z, Zhu S, et al. BSKs mediate signal transduction from the receptor kinase BRI1 in Arabidopsis. Science. 2008;321(5888):557–60.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  56. Zhang B, Chen HW, Mu RL, Zhang WK, Zhao MY, Wei W, et al. NIMA-related kinase NEK6 affects plant growth and stress response in Arabidopsis. Plant J. 2011;68(5):830–43.

    Article  PubMed  CAS  Google Scholar 

  57. DeYoung BJ, Bickle KL, Schrage KJ, Muskett P, Patel K, Clark SE. The CLAVATA1-related BAM1, BAM2 and BAM3 receptor kinase-like proteins are required for meristem function in Arabidopsis. Plant J. 2006;45(1):1–16.

    Article  PubMed  CAS  Google Scholar 

  58. Kinoshita A, Betsuyaku S, Osakabe Y, Mizuno S, Nagawa S, Stahl Y, et al. RPK2 is an essential receptor-like kinase that transmits the CLV3 signal in Arabidopsis. Development. 2010;137(22):3911–20.

    Article  PubMed  CAS  Google Scholar 

  59. Escobar-Restrepo JM, Huck N, Kessler S, Gagliardini V, Gheyselinck J, Yang WC, et al. The FERONIA receptor-like kinase mediates male-female interactions during pollen tube reception. Science. 2007;317(5838):656–60.

    Article  PubMed  CAS  Google Scholar 

  60. Xu T, Dai N, Chen J, Nagawa S, Cao M, Li H, et al. Cell surface ABP1-TMK auxin-sensing complex activates ROP GTPase signaling. Science. 2014;343(6174):1025–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  61. Li H, Yang WC. RLKs orchestrate the signaling in plant male-female interaction. Sci China Life Sci. 2016;59(9):867–77.

    Article  PubMed  CAS  Google Scholar 

  62. Wang J, Qu B, Dou S, Li L, Yin D, Pang Z, et al. The E3 ligase OsPUB15 interacts with the receptor-like kinase PID2 and regulates plant cell death and innate immunity. BMC Plant Biol. 2015;15:49.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  63. Oh E, Seo PJ, Kim J. Signaling peptides and receptors coordinating plant root development. Trends Plant Sci. 2018;23(4):337-51.

  64. Hematy K, Sado PE, Van Tuinen A, Rochange S, Desnos T, Balzergue S, et al. A receptor-like kinase mediates the response of Arabidopsis cells to the inhibition of cellulose synthesis. Curr Biol. 2007;17(11):922–31.

    Article  PubMed  CAS  Google Scholar 

  65. Feng W, Kita D, Peaucelle A, Cartwright HN, Doan V, Duan Q, et al. The FERONIA receptor kinase maintains Cell-Wall integrity during salt stress through ca(2+) signaling. Curr Biol. 2018;28(5):666–675 e5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  66. Rogers HJ. From models to ornamentals: how is flower senescence regulated? Plant Mol Biol. 2013;82(6):563–74.

    Article  PubMed  CAS  Google Scholar 

  67. Bieleski R, Elgar J, Heyes J. Mechanical aspects of rapid flower opening in Asiatic lily. Ann Bot. 2000;86(6):1175–83.

    Article  Google Scholar 

  68. Liang H, Mahadevan L. Growth, geometry, and mechanics of a blooming lily. Proc Natl Acad Sci U S A. 2011;108(14):5516–21.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Wang H, Wang H. Phytochrome signaling: time to tighten up the loose ends. Mol Plant. 2015;8(4):540–51.

    Article  PubMed  CAS  Google Scholar 

  70. Li QH, Yang HQ. Cryptochrome signaling in plants. Photochem Photobiol. 2007;83(1):94–101.

    Article  PubMed  CAS  Google Scholar 

  71. Hohm T, Preuten T, Fankhauser C. Phototropism: translating light into directional growth. Am J Bot. 2013;100(1):47–59.

    Article  PubMed  CAS  Google Scholar 

  72. Yin R, Ulm R. How plants cope with UV-B: from perception to response. Curr Opin Plant Biol. 2017;37:42–8.

    Article  PubMed  CAS  Google Scholar 

  73. Liu H, Liu B, Zhao C, Pepper M, Lin C. The action mechanisms of plant cryptochromes. Trends Plant Sci. 2011;16(12):684–91.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  74. Matsushika A, Makino S, Kojima M, Yamashino T, Mizuno T. The APRR1/TOC1 quintet implicated in circadian rhythms of Arabidopsis thaliana: II. Characterization with CCA1-overexpressing plants. Plant Cell Physiol. 2002;43(1):118–22.

    Article  PubMed  CAS  Google Scholar 

  75. Green RM, Tobin EM. The role of CCA1 and LHY in the plant circadian clock. Dev Cell. 2002;2(5):516–8.

    Article  PubMed  CAS  Google Scholar 

  76. Lu SX, Liu H, Knowles SM, Li J, Ma L, Tobin EM, et al. A role for protein kinase casein kinase2 alpha-subunits in the Arabidopsis circadian clock. Plant Physiol. 2011;157(3):1537–45.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  77. Liu J, Li J, Wang H, Fu Z, Liu J, Yu Y. Identification and expression analysis of ERF transcription factor genes in petunia during flower senescence and in response to hormone treatments. J Exp Bot. 2011;62(2):825–40.

    Article  PubMed  CAS  Google Scholar 

  78. Chen WH, Li PF, Chen MK, Lee YI, Yang CH. FOREVER YOUNG FLOWER negatively regulates ethylene response DNA-binding factors by activating an ethylene-responsive factor to control Arabidopsis floral organ senescence and abscission. Plant Physiol. 2015;168(4):1666–83.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  79. Pei H, Ma N, Tian J, Luo J, Chen J, Li J, et al. An NAC transcription factor controls ethylene-regulated cell expansion in flower petals. Plant Physiol. 2013;163(2):775–91.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  80. Braam J. In touch: plant responses to mechanical stimuli. New Phytol. 2005;165(2):373–89.

    Article  PubMed  Google Scholar 

  81. Kitakura S, Vanneste S, Robert S, Lofke C, Teichmann T, Tanaka H, et al. Clathrin mediates endocytosis and polar distribution of PIN auxin transporters in Arabidopsis. Plant Cell. 2011;23(5):1920–31.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  82. Sanchez-Rodriguez C, Rubio-Somoza I, Sibout R, Persson S. Phytohormones and the cell wall in Arabidopsis during seedling growth. Trends Plant Sci. 2010;15(5):291–301.

    Article  PubMed  CAS  Google Scholar 

  83. Amin S, Prentis PJ, Gilding EK, Pavasovic A. Assembly and annotation of a non-model gastropod (Nerita melanotragus) transcriptome: a comparison of de novo assemblers. BMC Res Notes. 2014;7:488.

    Article  PubMed  PubMed Central  Google Scholar 

  84. Du Z, Zhou X, Ling Y, Zhang Z, Su Z. agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010;38(Web Server issue):W64–70.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  85. Tian T, Liu Y, Yan H, You Q, Yi X, Du Z, et al. agriGO v2.0: a GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Res. 2017;45(W1):W122–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  86. Ulitsky I, Maron-Katz A, Shavit S, Sagir D, Linhart C, Elkon R, et al. Expander: from expression microarrays to networks and functions. Nat Protoc. 2010;5(2):303–22.

    Article  PubMed  CAS  Google Scholar 

  87. Pham VN, Kathare PK, Huq E. Phytochromes and Phytochrome interacting factors. Plant Physiol. 2018;176(2):1025–38.

    Article  PubMed  CAS  Google Scholar 

  88. Inoue K, Araki T, Endo M. Circadian clock during plant development. J Plant Res. 2018;131(1):59–66.

    Article  PubMed  Google Scholar 

  89. Weijers D, Friml J. SnapShot: auxin signaling and transport. Cell. 2009;136(6):1172.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

We thank Deshu Lin for microscopy assistance.

Funding

This work was supported by the National Key Research and Development Program of China (2017YFA0506100), Program of FAFU-Horticulture College Young Scholar, National Young 1000 Talent Grant and Fujian 100 Talent Program.

Availability of data and materials

The datasets generated or analyzed during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations

Authors

Contributions

MK: performed the experiments of time-lapse photography, scanning electronic microscopy and RNA-seq transcriptome analysis; ZG: assistance of RNA-seq transcriptome analysis, drawing the cartoon models of Figs. 5 and 6; JC: assistance of RNA-seq transcriptome analysis; YQ: assistance of scanning electronic microscopy observation and quantification; LZ and XC: contributed waterlily materials and initiate the project; ZG and XC: designed the project and wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xu Chen.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

The intermediate petal segment determines the flower opening and closure rhythm.Corresponding to Fig. 1. (A) The opening and closure rhythm of Nymphaea colorata Peter’s flower was tracked from 7:30 to 18:00. (B) Without the sepal, the floral opening-closure movement was tracked when only 0.5 cm of the petal was left. White arrows highlighted the petal movement. (C) Measurement methods of floral opening angles. Scale bar, 25 mm (A-C). (TIF 6190 kb)

Additional file 2:

Global transcriptome description of the flower opening and closure. Corresponding to Fig. 4. Global analysis of GO items and KEGG pathways during the flower opening and closure processes are shown. The transcriptome at 6:00 (T1) was used as the control, and the transcriptomes at 7:00 (T2), 10:00(T3), 14:00 (T4) and 18:00 (T5) were individually compared with each other (T/T comparisons). The total changed numbers are listed below. Red represents the upregulated groups, and blue marks the downregulated groups (TIF 1165 kb)

Additional file 3:

The involvement of light signaling during flower opening. (A-B) GO items and heatmaps provided the transcriptome profiles of cluster 1, including genes involved in photosynthesis (list in Additional file 4: Table S1). LHCB, Chlorophyll a-b binding protein; PSB, Oxygen-evolving enhancer protein; PSA, Photosystem I reaction center subunit V; CAB, Chlorophyll a-b binding protein; CAP, Chlorophyll a-b binding protein CP24; PPD, PsbP domain-containing protein; (C) Transcript profiling of light signaling network. The basic red/far-red light signaling cascade is based on the PHYs-PIF signaling module, as referred in [87]. Blue light signaling primarily consists of the CRY-COP1/SPA1-HY5 signaling module for photomorphogenesis, CRY-COP1/SPA1-GI-CO-FT module for photoperiodic flowering, PHOTs-mediated signaling for phototropism, as referred to [8, 70, 73]. UVB-triggered signaling is established by UVR8-COP1/SPA1-HY5 module, as referred in [72]. The input light signals additionally stimulate the endogenous circadian oscillator, which is feedback regulated by the transcriptional loop of the APRRs-TOCCCA1-LHY-ELF4-FKF1-ZTL clock signaling module, as referred to [88]. PHY, Phytochrome; CRY, Cryptochrome; PHOT, Phototropin; UVR8, UV-B resistance 8; PIF3, Phytochrome-interacting factor 3; ELF, Early flowering; COP1, Constitutive photomorphogenic 1; SPA1, Suppressor of phytochrome A1; HY5, Elongated hypocotyl 5; FKF1, Flavin Binding, Kelch Repeat, F-BOX1; ZTL, ZEITLUPE; APRR, Pseudo-Response Regulators; TOC1, Timing of CAB expression1; LHY, Late Elongated Hypocotyl; CCA1, Circadian Clock Associated 1; GI, GIGANTEA; CO, CONSTANS; FT, Protein Flowering Locus T; CKL, Casein kinase 1-like; EXP, Expansin. The expression patterns of these genes are marked as heatmaps (the corresponding genes were listed in Additional file 4: Table S4 according to the labelled number “28” to “42” in each cluster). (TIF 2690 kb)

Additional file 4:

Table S1. Transcriptome profile of photosynthesis-regulators, corresponding to Additional file 3: Figure S3B. Table S2. Transcriptome profile of vesicle trafficking-regulators, corresponding to Additional file 5: Figure S4A. Table S3. Transcriptome profile of receptor kinases, corresponding to Additional file 5: Figure S4B. Table S4. Transcriptome profile of light signaling-regulators, corresponding to Additional file 3: Figure S3C. Table S5. Transcriptome profile of transcriptional factors, corresponding to Additional file 6: Figure S5B. Table S6. Transcriptome profile of protein synthesis-related genes, corresponding to Additional file 6: Figure S5D. Table S7. Transcriptome profile of transcriptional factors, corresponding to Additional file 6: Figure S5E. Table S8. NGS reads for waterlily RNA. Table S9. Statistical result of waterlily RNA transcriptome assembly. (DOCX 65 kb)

Additional file 5:

Transcriptome description of vesicle trafficking and receptor kinase-mediated cascades. (A) Heatmap provided the transcriptome profiles of vesicle trafficking-related genes in cluster 3. The corresponding genes were listed in Additional file 4: Table S2. (B) Heatmap provided the transcriptome profiles of RLKs in cluster 8. The corresponding genes were listed in Additional file 4: Table S3. (TIF 640 kb)

Additional file 6:

Description of flower opening-associated signaling events with downregulations. (A-F) GO items and heatmaps provided the transcriptome profiles of cluster 2, 4 and 6, mainly including transcriptional factors and protein synthesis-related genes (list in Additional file 4: Tables S5-S7). In B, ERF, Ethylene-responsive transcription factor; WRKY, WRKY transcription factor; NAC, NAC domain-containing protein; MADS, MADS-box transcription factor; HSFA, Heat stress transcription factor; BZIP, Basic leucine zipper; CDF, Cyclic dof factor. In D, RPL, ribosomal protein L; RPS, ribosomal protein S; TIF, Eukaryotic translation initiation factor. In E, Homeo-box, Homeobox-leucine zipper protein; Zinc finger, Dof zinc finger protein; ERF, Ethylene-responsive transcription factor; RAP, Ethylene-responsive transcription factor RAP; CRF, Ethylene-responsive transcription factor CRF; BBM, AP2-like ethylene-responsive transcription factor BBM; IDD, indeterminate-domain; NAC, NAC domain-containing protein; WRKY, WRKY transcription factor; GATA, GATA transcription factor; TGA, Transcription factor TGA. (TIF 3371 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ke, M., Gao, Z., Chen, J. et al. Auxin controls circadian flower opening and closure in the waterlily. BMC Plant Biol 18, 143 (2018). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-018-1357-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-018-1357-7

Keywords