Abstract
Throughout development the Arabidopsis shoot apical meristem successively undergoes several major phase transitions such as the juvenile-to-adult and floral transitions until, finally, it will produce flowers instead of leaves and shoots. Members of the Arabidopsis SBP-box gene family of transcription factors have been implicated in promoting the floral transition in dependence of miR156 and, accordingly, transgenics constitutively over-expressing this microRNA are delayed in flowering. To elaborate their roles in Arabidopsis shoot development, we analysed two of the 11 miR156 regulated Arabidopsis SBP-box genes, i.e. the likely paralogous genes SPL9 and SPL15. Single and double mutant phenotype analysis showed these genes to act redundantly in controlling the juvenile-to-adult phase transition. In addition, their loss-of-function results in a shortened plastochron during vegetative growth, altered inflorescence architecture and enhanced branching. In these aspects, the double mutant partly phenocopies constitutive MIR156b over-expressing transgenic plants and thus a major contribution to the phenotype of these transgenics as a result of the repression of SPL9 and SPL15 is strongly suggested.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
Introduction
During maturation, plants pass through several developmentally distinct growth phases in which the shoot gradually gains reproductive competence (Poethig 1990). After the transition from embryonic to postembryonic growth, plants undergo at least two further phase transitions, the vegetative as well as the reproductive phase change. During vegetative growth, rosette leaves are initiated at the flanks of the shoot apical meristem (SAM) with a certain frequency, referred to as plastochron (Erickson and Michelini 1957). After going through the reproductive phase transition, also known as the floral transition, the SAM starts to initiate floral buds instead of leaves. In Arabidopsis, as in many other plants showing day length dependent flowering, the floral transition is preceded by a transition from juvenile to adult growth. This switch, known as the vegetative phase change, is physiologically defined as achieving competence to respond to photoperiodic induction of flowering (Poethig 1990).
The transition from juvenile to adult growth is gradual and rather subtle but generally can be followed by several morphological markers. In Arabidopsis, for example, leaves produced in the juvenile phase have long petioles, are small, round and lack abaxial trichomes. In contrast, short petioles, elliptical anatomy and the development of trichomes on the abaxial side represent adult traits (Telfer et al. 1997).
Regulation of these developmental transitions is largely dependent on (changes in) environmental cues such as day length, light intensity and temperature, as well as on endogenous factors such as the plant hormone gibberellin (Telfer et al. 1997). Whereas the molecular genetic mechanisms underlying the floral transition are already worked out in increasing detail (Komeda 2004), it has only been recently that we begin to understand the molecular genetic basis of the vegetative phase change. Most genes suggested to play a role in promoting the latter phase change have been identified by the analysis of mutants showing a precocious onset of adult traits and intriguingly, link vegetative phase change to RNA silencing pathways. These genes include the Arabidopsis ortholog of exportin 5/MSN5, HASTY (HST; Telfer and Poethig 1998; Bollman et al. 2003), the zinc-finger-domain protein encoding locus SERRATE (SE; Clarke et al. 1999) and ZIPPY (ZIP), an AGO-family member (Hunter et al. 2003). More recently, screens for mutations with zip-like phenotypes resulted in alleles of SUPRESSOR OF GENE SILENCING3 (SGS3) and RNA-DEPENDENT POLYMERASE6 (RDR6), both genes required for posttranscriptional gene silencing (PTGS) and acting in the same pathways as ZIP and HST (Peragine et al. 2004). Furthermore, a precocious vegetative phase change has also been found in dicer-like 4 (dcl4) mutants (Gasciolli et al. 2005; Xie et al. 2005; Yoshikawa et al. 2005).
One explanation for the observed effects could be that target genes of this silencing pathway play a positive role on the vegetative phase change and their down-regulation consequently promotes juvenility. Hence, mutations in genes involved in this silencing pathway, as the ones described above, cause an accelerated vegetative development.
In line with this idea, members of the plant specific SBP-box gene transcription factor family have been implicated in promoting vegetative and floral phase transitions. In particular, overexpression of the Arabidopsis SBP-box gene SPL3 leads to early flowering and a significant earlier appearance of abaxial trichomes on the rosette leaves (Cardon et al. 1997; Wu and Poethig 2006). Interestingly, together with 10 of 16 other family members, SPL3 expression is post-transcriptionally controlled by miR156 and probably also by the very closely related miR157 (Rhoades et al. 2002; Schwab et al. 2005; Wu and Poethig 2006; Gandikota et al. 2007). Consistent with its role of down-regulating SPL3 and related SPL target-genes, constitutive overexpression of miR156 encoding loci has been shown to cause the production of a significantly larger number of leaves with juvenile characteristics and a delay in flowering (Schwab et al. 2005; Wu and Poethig 2006).
Although the available data clearly point to a regulatory role for the miRNA regulated SPL genes in the temporal development of the Arabidopsis shoot, the contribution of the single genes to the described phenotypes remains to be determined. Therefore, we identified and isolated mutant alleles for single SPL genes. In comparison to other miR156 targeted SPL genes, available expression data (AtGenExpress; Schmid et al. 2005) show SPL9 and SPL15 to be already quite active in the vegetative shoot apex. Accordingly, their mutant phenotypes were found to affect vegetative development. Here we report the mutant analysis of SPL9 and SPL15, two likely paralogous members of the SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factor family (Cardon et al. 1999), and discuss their redundant regulatory role on the vegetative phase change and the temporal initiation of rosette leaves.
Materials and methods
Plant material and plant growth conditions
All of the genetic stocks described in the paper were in Columbia background. The T-DNA insertion lines SALK_006573 (spl9-2), SALK_074426 (spl15-1) and SALK_138712 (spl15-2) were obtained at the Nottingham Arabidopsis Stock Centre (NASC). The T-DNA insertion lines GABI-Kat 544F04 (spl9-3) and WiscDsLox 457 (spl15-3) were obtained from GABI-Kat and the Arabidopsis Biological Research Centre (ABRC), respectively. Insertion mutant information for NASC- or ARBC-lines was obtained from the SIGNAL website at http://signal.salk.edu. Plants homozygous for the T-DNA insertions were identified by PCR using T-DNA left border- and gene-specific primers. T-DNA specific left border primers for SALK, GABI-Kat and WiscDsLox T-DNAs were 5′-GCGTGGACCGCTTGCTGCAACT-3′, 5′-ATATTGACCATCATACTCATTGC-3′ and 5′-TGGCAGGATATATTGTGGTGTAAACA-3′, respectively. In combination with the respective left border primer we used the following gene-specific primers: 5′-GCTATGGCTTAAGCCTTAAGTTAAAAGG-3′ for SALK_006573, 5′-CGTAGCTGTCGTGGACTAGTGTCAATC-3′ for SALK_074426 and SALK_138712, 5′-AACCTCTGTTCGATACCAGCCACAG-3′ for GABI-Kat 544F04 and 5′-AGCCATTGTAACCTTATCGGAGAATGAG-3′ for WiscDsLox 457. The stable En-1 insertion mutant 5ABA33-H1 (spl9-1) was obtained from the ZIGIA-population (Unte 2001). Plants homozygous for a four base pair insertion in the first exon of SPL9 caused by the excision of the En-1 transposon were backcrossed with wild type twice to obtain plants exclusively containing the four base-pair insertion without any further transposon contamination. In order to identify plants containing the mutation we used the following primer combination: 5′-AGTAAGAGGAAACCACCATGGAGATGG-3′ (forward) and 5′-AACCTTCCACTTGGCACCTTGGTATA-3′ (reverse, recognises the insertion).
All plants were grown in plastic trays or pots filled with ready-to-use commercial, pre-fertilized soil mixture (Type ED73, Werkverband eV, Sinntal-Jossa, Germany). For stratification, seeds were kept on moist paper at 4°C in the dark for 4–5 days before transferring to soil (i.e. “sowing”) in growth chambers at 22°C, 50% relative humidity. Germination and cultivation of the plants in long-day conditions (16 h light, 8 h dark) were either under approx. 70 μE/cm2/s (LD1) or 175 μE/cm2/s (LD2) light provided by fluorescent tubes (L58W/840 and L58W/25 Osram, Munich, Germany). Plants in short-day (SD) conditions (8 h light, 16 h dark) were cultivated under approx. 450 μE/cm2/s light. To determine sensitivity to photoperiodic induction of flowering, stratified seeds were germinated in a modified SD with a 9 h light period. The developing plants were kept in these conditions for 21 days before they were transferred to similar growth chambers with continuous light provided by Osram HQIT 400 W lamps. Batches of plants were returned to the modified SD conditions after 1, 3 and 5 days.
Phenotypic analysis
Flowering time was measured as the time between sowing and anthesis (opening of first flower). Bolting time was recorded when the main inflorescence had reached a height of 0.5 cm. Inflorescence height was measured between the rosette and the first flower of the main inflorescence of plants with their first siliques fully ripened. In order to count the number of site shoots, all site shoots longer than 0.5 cm were scored. Abaxial trichomes were scored using a Leica MZFLIII stereomicroscope (Wetzlar, Germany). An estimation of the rosette leaf initiation rate (L/D day−1) was obtained by dividing the number (L) of rosette leaves having reached at least 0.5 cm in length through the number of days (D) between sowing and determination. Note that this value reflects but does not equalize (average) plastochron as it should be corrected for the true start of initiation of the first leaf as well as the time needed for the first adult leaf to reach a size of 0.5 cm.
Histological analysis
Apical regions were isolated from plants grown in SD for 41 days by trimming with a razor blade. The tissue was fixed in 4% formaldehyde/0.1 M PO4, pH 7.0 for 48 h and embedded in paraffin using a Leica ASP300 tissue processor (Wetzlar, Germany). Embedded apices were cut into 8 μm thin cross sections using a Jung Autocut 2055 and photographed using a Zeiss Axiophot microscope (Göttingen, Germany) equipped with a KY-F5U 28CCD camera (JVC, Yokohama, Japan). The first cross section in which the apex was visible plus two successive sections were used to determine the diameter of the apical region. From these three measured values the average was taken for comparison. In addition, cross sectional area and circularity factor (=4π A P −2, A is area, P is perimeter) of the leaf primordia, outlined by hand on the photographs, were determined as well. The measurements were performed with help of the program ImageJ 1.35s (Wayne Rasband, National Institutes of Health, USA).
GA3 treatments
Col-0, the spl9 spl15 double mutant and the 35S::MIR156b overexpressor were grown in LD1 conditions. Immediately after germination, half of the plants were treated by spraying 100 μM GA3, 0.02% Tween 20 and this was repeated twice per week until they started flowering. The other half of the plants was similarly treated with 0.02% Tween 20.
Phylogenetic comparison
Multiple alignments of amino acid sequences were generated by the program ClustalW of the MacVector 7.2.2 software package (Accelrys Ltd., Cambridge, UK) using the BLOSUM 30 matrix with an open gap penalty of 10 and an extend gap penalty of 0.05. Only the SBP-domain was used for the phylogenetic reconstruction. The tree was constructed using the neighbour-joining algorithm of the MacVector 7.2.2 software package.
Quantitative real-time PCR analysis
To perform quantitative RT-PCR (using the iQ5 real-time PCR detection system, Bio-Rad, Munich, Germany) apical regions were collected (roots and as much of the leaves as possible were removed using tweezers) of plants cultivated 5, 9, 13, 27 and 32 days after sowing in LD1 conditions. Total RNA was extracted using the RNeasy plant mini kit (Qiagen, Hilden, Germany), including an on-column DNase digestion. First-strand cDNA was synthesized using SuperScript III RNase H reverse transcriptase (Invitrogen) starting with 2 μg of total RNA primed with an oligo(T)12–18 primer (Gibco BRL, Karlsruhe, Germany). SPL9-specific primers, 5′-AGAACATTGGATACAACAGTGATGAGG-3′ (forward) and 5′-GTTTGAGTCGCCAATTCCCTTGTAGC-3′ (reverse) as well as SPL15-specific primers, 5′-TTGGGAGATCCTACTGCGTGGTCAACC-3′ (forward) and 5′-AGCCATTGTAACCTTATCGGAGAATGAG-3′ (reverse), were designed to generate a PCR product of 171 and 300 bp, respectively. Based on the analysis of Czechowski et al. (2005) PP2A expression was used as reference for transcript normalization with the primer pair 5′-TAACGTGGCCAAAATGATGC-3′ (forward) and 5′-GTTCTCCACAACCGCTTGGT-3′ (reverse). The PCR efficiencies for SPL9, SPL15 and PP2A primers were determined to be 96, 95.5 and 96%, respectively. Quantifications, in triplicate, were performed using the Brilliant SYBRGreen QPCR kit (Stratagene, La Jolla, CA, USA), according to the manufacturer’s protocol, in a final volume of 25 μl. PCR was carried out in 250 μl optical reaction vials (Stratagene) heated for 10 min at 95°C to hot-start the Taq polymerase, followed by 40 cycles of denaturation (30 s at 95°C), annealing (30 s at 58°C) and extension (30 s at 72°C).
Semi-quantitative RT-PCR analysis
Total RNA was extracted from seedlings of the Col-0 wild-type, mutant- and transgenic lines using the RNeasy plant mini kits (Qiagen). RT-PCR with equal amounts of RNA was performed using the one-step RT-PCR kit (Qiagen). SPL9 knockout lines were identified using the following primer pair: 5′-GGTCGGGTCAGTCGGGTCAGATACC-3′ (forward) and 5′-ACTGGCCGCCTCATCACTCTTGTATCC-3′ (reverse). SPL9 mRNA is expected to yield a 415 bp fragment whereas genomic SPL9 DNA is expected to yield a 1,138 bp fragment. SPL15 knockout lines were identified using the following primer pair: 5′-AGAAGCAAGAACCGGGTCAATACC-3′ (forward) and 5′-AGCCATTGTAACCTTATCGGAGAATGAG-3′ (reverse). SPL15 mRNA is expected to yield a 666 bp fragment whereas genomic SPL15 DNA is expected to yield a 1,004 bp fragment. RT-PCR of the loading control (RAN3; At5g55190) was performed with the primer pair 5′-ACCAGCAAACCGTGGATTACCCTAGC-3′ (forward) and 5′-ATTCCACAAAGTGAAGATTAGCGTCC-3′ (reverse) to yield a fragment of 531 bp when derived from RAN3 mRNA (genomic RAN3 is expected to yield a 1,314 bp fragment).
Statistics
Graphical representations of numerical data were generated with the Microsoft Excel program (Microsoft Germany, Munich) and statistical tests were performed using the Student’s t-Test within this program. P-values lower than 0.05 were considered to be statistically relevant and the data involved to represent significant differences.
Remaining techniques and methods
Standard molecular biology techniques were performed as described by Sambrook et al. (1989). Graphical plots and digital photographic images were cropped and assembled using Adobe Photoshop (Adobe Systems, San Jose, CA, USA).
Results
Molecular characterization of the Arabidopsis SBP-box genes SPL9 and SPL15
Mutant alleles for the Arabidopsis SBP-box genes SPL9 (At2g42200) and SPL15 (At3g57920) were obtained from screening publicly available electronic databases and seed stock centres for transposon or T-DNA tagged SPL genes. For SPL9, we identified three insertion alleles designated as spl9-1 to -3 and confirmed the nature and position of their mutations (see “Materials and methods”; Fig. 1a). The first allele, spl9-1, was identified in the En-transposon mutagenised ZIGIA population (Baumann et al. 1998; Unte 2001) and most likely resulted from the excision of an inserted En-1 transposon leaving behind a 4-bp insertion footprint in the first exon. The result of this is a frame shift in the coding sequence and the generation of a stop-codon 86 base pairs after the insertion site. Both spl9-2 and spl9-3 represent T-DNA insertion mutant alleles identified within, respectively, the SALK collection (Alonso et al. 2003) and the GABI-Kat collection (Li et al. 2007).
Also three independent T-DNA insertion lines for SPL15 could be obtained and confirmed (see “Materials and methods”; Fig. 1a). Two alleles designated as spl15-1 and spl15-2 were identified within the SALK collection and one, spl15-3, within the WiscDsLox T-DNA collection.
According to data available from the AtGenExpress micro-array database (Schmid et al. 2005), both SPL9 and SPL15 transcript levels increase during development and are preferentially found in the shoot apical region and in young flowers. We confirmed this temporal expression pattern with the help of qRT-PCR (Fig. 1b). In LD1 growing conditions (see “Materials and methods”), SPL9 and SPL15 transcript levels remain comparable during the first 2–3 weeks. Thereafter, the expression level of SPL9 starts to increase followed by that of SPL15. Around 32 days after sowing (DAS), at about the time Col-0 plants have undergone their reproductive phase transition, SPL9 transcript levels have become approximately two and a half times higher in comparison to SPL15 and six times in comparison to day 5. Arabidopsis lines carrying as a transgene a genomic fragment encompassing the locus for SPL15 and with a GUS reporter gene inserted downstream of the ATG start codon, confirmed the predominantly apical expression of SPL15 (Supplementary Fig. 1).
RT-PCR performed on mRNA isolated from whole seedling plants homozygous for any of the three SPL9 or SPL15 mutant alleles (see “Materials and methods”) did not result in the detection of RNA derived of the respective genes (Fig. 1c). This strongly suggests that all mutant alleles isolated represent functional null-alleles. Accordingly, plants homozygous for any of the three spl9 mutant alleles showed highly identical phenotypes, as did all three homozygous spl15 mutants (see phenotypic analysis below; Supplementary Fig. 2). Allelic tests confirmed that the observed phenotypes are indeed due to mutation in either SPL9 or SPL15, respectively (data not shown).
With over 75% of their amino acid residues identical, SPL9 and SPL15 show high similarity on the level of their proteins. Also a phylogenetic comparison based on the SBP-box of all 17 SPL genes in Arabidopsis revealed SPL9 and SPL15 as most closely related and most likely forming a pair of paralogous genes (Fig. 2). Based on this close relationship some degree of functional redundancy could be expected and, therefore, we created double mutant lines to uncover such redundancy. To ascertain that phenotypic changes in the mutant plants are solely due to the loss-of-function of SPL9 and SPL15 we generated two different homozygous double mutant lines with the allelic combinations spl9-1 spl15-1 and spl9-2 spl15-2, respectively. Both lines exhibit the same phenotype as described in the next section. For further detailed analysis the spl9-1 spl15-1 line was chosen and in the following referred to as spl9 spl15 for simplicity.
Phenotypic analysis of spl9 and spl15 mutants
For a phenotypic analysis, we compared spl9 and spl15 single mutants, spl9 spl15 double mutants to Col-0 wild type as well as to an 35S::MIR156b transgenic line (kindly provided by D. Weigel and R. Schwab).
An interesting aspect of the MIR156b over-expressing plants, as already noticed by Schwab and co-workers (2005) is an increased rate of rosette leaf initiation which, in combination with a modest delayed flowering, results in the obvious denser rosettes of fully developed plants (Fig. 3a). In addition, advanced 35S::MIR156b plants became very bushy (Fig. 3c). We found these phenotypic aspects also displayed by the spl9 spl15 double mutant, albeit less pronounced (Fig. 3a, b). To quantify the contribution of SPL9 and SPL15 to these phenomena, we compared the number of rosette and cauline leaves of the respective single and double mutants and of the MIR156b overexpressor to wild type (Table 1). Whereas in LD2 growing conditions, the 35S::MIR156b line produced ca. eleven more rosette leaves in comparison to wild type, the single mutant lines produced, on average, only 1–2 rosette leaves more. Again, with ca. six extra rosette leaves, the spl9 spl15 double mutant differed more from wild type than the single mutants and showed a stronger tendency towards the phenotype of the MIR156b overexpressor. The number of cauline leaves remained very comparable among all mutants and wild type, although some reduction may be observed particularly in the spl9 mutants.
Also with respect to the development of side shoots, the spl9 spl15 double mutant differed more from wild type than the single mutants. In fact, spl9-1 and spl15-1 single mutants were found not to differ significantly from Col-0 plants that had formed, on average, 0.9 ± 0.6 side shoots of at least 0.5 cm in length by the time that the first siliques ripened. With an average of 2.1 ± 1.1 side shoots, the spl9 spl15 double mutant did significantly differ from wild type as did the 35S::MIR156b transgenic line with, on average, 4.1 ± 0.8 side shoots. Taken together, the phenotypic data of the spl9 spl15 double mutant clearly suggests a redundant function of SPL9 and SPL15 in shoot development and in the maintenance of apical dominance.
MiR156 is assumed to target, besides SPL9 and SPL15, exclusively other SPL genes (Rhoades et al. 2002). These too were shown to be down regulated in MIR156b over-expressing plants (Schwab et al. 2005). As in comparison to the spl9 spl15 double mutant the MIR156b over-expressor displays an even more severe aberrant phenotype, it can also be deduced that in addition to SPL9 and SPL15, other miR156-controlled SPL genes act redundantly to control shoot development and apical dominance.
In addition to the number of leaves formed before the appearance of the first flowers, we also determined for the same plants the time they needed to bolt as well as to anthesis (Table 1). On average, the spl9 and spl15 single mutants behaved similar to wild type but, as expected based on the data of Schwab et al. (2005), the 35S::MIR156b line bolted and flowered somewhat later. The spl9 spl15 double mutants showed an intermediate behaviour. Whereas for the single mutants the few more leaves formed may be accounted for by the slight delay in the transition to flowering, this delay is unlikely to explain the increased rosette leaf number of the spl9 spl15 double mutant. In line with the observation of Schwab and co-workers (2005) who reported a leaf-initiation rate per day in SD of 2.2 vs. 1.4 for the MIR156b overexpressor and the wild type, respectively, this is probably best explained by assuming a shortened plastochron during vegetative growth.
To uncover a possible cause or consequence for this increased rate of leaf initiation, we microscopically examined cross sections of the vegetative shoot apex to determine size and phyllotaxy of the spl9 spl5 double mutant and the MIR156b overexpressor and compared these to wild type. To this purpose, plants were grown for 41 days in SD conditions, whereafter the number of rosette leaves having reached at least 0.5 cm in length were recorded and their apices dissected, fixed and embedded in paraffin (see “Materials and methods”). At this age, Col-0 plants were found to have formed on average 24.5 leaves of 0.5 cm or more, the double mutant 33.8 and the MIR156b overexpressor plants already 43.4 (Fig. 4a). As the plants were of the same age, these differences most likely reflect differences in plastochron. Alternatively, one may assume large temporal differences per genotype concerning initiation of the first leaf and/or development of the last leaf recorded to have reached 0.5 cm in length. However, we obtained no indications for such discrepancies and noted increasing differences in rosette density during the entire vegetative growth phase of wild type and mutants. From these data, a relative 1.8-fold (43.4/24.5) increase in leaf initiation rate of the 35S::MIR156b transgenics over wild type can be deduced, a value that quite well matches the observation of Schwab et al. (2005). Leaf initiation rate of the double mutant seem to be increased by a factor of 1.4 (33.8/24.5) in comparison to wild type.
After sectioning the paraffin embedded material, a small but not significant difference in average SAM-diameter of the spl9 spl15 double mutant and Col-0 wild type could be observed (Fig. 4b, e–f). However, with an average diameter of 104 μm, the MIR156b overexpressor showed also a slight but yet significant (P < 0.05) decrease in its SAM size compared to Col-0 (Fig. 4b, e, g). Furthermore, both the spl9 spl15 double mutant and the MIR156b overexpressor exhibited the same phyllotaxy as wild type with rosette leaves initiated either clock- or anticlockwise with an angle of divergence of about 137.5° between successive leaves and forming a spiral lattice with a parastichy pair (3,5) (Fig. 4e–g). From these observations, it is concluded that the observed shorter plastochron is neither the result nor the cause of an altered phyllotaxy in the spl9 spl15 double mutant or the MIR156b overexpressor. The shortened plastochron, however, seems to correlate with a reduced SAM size.
As obvious from cross sections shown in Fig. 4e, g, the young leaves of the MIR156b overexpressor appear more roundish in shape in comparison to wild type leaves at similar positions. In particular, the vacuolated cells surrounding their midveins seem larger and the developing laminas reduced, i.e. represented by less small cytoplasm rich cells along their lateral margins. In addition, the stipules of the MIR156b overexpressor seem to be more prominent. In these aspects of leaf development, the spl9 spl15 double mutant seems to behave intermediate (Fig. 4e).
To quantify the difference in shape and size of the leaf primordia, we determined their circularity and cross sectional area (see “Materials and methods”) starting from the first leaf cross section found to be separated from the apical meristem. To reduce effects due to imperfect cross sectioning, i.e. not absolute perpendicular to the longitudinal axes, as well as in correlating sequentially numbered primordia of different sections, we averaged the values obtained over 10 successive primordia. As shown in Fig. 4c, circularity of the youngest 10 leaf primordia was highly similar between the different genotypes. Circularity of subsequent older primordia decreased in all, however, more rapidly in wild type such that, on average, leaf 10–20 did differ significantly between the genotypes. Interestingly, it cannot be excluded that this difference is a direct consequence of a shortened plastochron in the mutant lines. In particular as the leaf initiation rate of the MIR156b overexpressor line lies roughly one and a half times above that of wild type. Accordingly, and with respect to absolute age, leaf 11–20 of the MIR156b overexpressor may be more comparable to leaves 6–15 of wild type and to which indeed no significant difference in circularity was found. However, in cross sectional area these young leaves differed significantly. On average, 6.5 ± 2.4 × 103 μm2 for wild-type leaves 6–15 and 26.7 ± 16.2 × 103 μm2 for leaf 11–20 of the MIR156b overexpressor.
It is known that the shape and other characteristics of newly formed leaves progressively change in correlation with the vegetative phase transition (Telfer et al. 1997). Furthermore, likely due to a changed plastochron, the correlation between leaf number and flowering seems to differ for the mutants and wild type. Therefore, we further investigated the possibility that the observed differences correlate with relative altered timing of the vegetative phase transition.
Functional analysis of SPL9 and SPL15 during the vegetative phase transition
In order to determine the timing of the vegetative phase change we used the absence or presence of abaxial trichomes on rosette leaves as a morphological marker for leaves formed during the juvenile or adult growth phase, respectively (Telfer et al. 1997). On average, in our LD2 growing conditions, the first abaxial trichomes developed on rosette leaf number six of Col-0 wild-type plants (Table 1). The spl9 spl15 double mutants displayed its first abaxial trichomes on leaf number twelve and the 35S::MIR156b overexpressor on leaf number 16. Although less than the spl9 spl15 double mutant, the respective single mutants also developed significantly more juvenile leaves than wild type (Table 1).
We distinguished juvenile and adult growth phase based on a phase dimorphism, i.e. in abaxial trichomes. However, on a plant-physiological level the juvenile phase in Arabidopsis is characterized as being incompetent to respond to photoperiodic induction of flowering (Poethig 1990). To determine if this competence was indeed affected, small populations of 20–22 plants of wild type, the spl9, spl15 single and double mutants and the MIR156b overexpressor were germinated and cultivated for 3 weeks in non-inductive SD conditions (see “Materials and methods”). The plants were then brought into continuous light and batchwise shifted back to SD after either 1, 3 or 5 days. Their flowering response was recorded within a 3-week period after this inductive treatment. Plants not flowering within this period also did not flower after 2 months like plants of a control group representing all genotypes that were kept continuously in SD. As shown in Table 2, the 5-day inductive treatment caused a flowering response in 100% of the plants of all genotypes. Three days also sufficed to induce all or almost all of the wild-type and single mutant plants, whereas the response of the double mutant and particular of the MIR156b overexpressor already declined. One day of continuous light, still enough to induce half or more of the wild-type and single mutant plants, did not induce flowering in any of the MIR156b overexpressor plants and only in one-tenth of the double mutants. These results thus demonstrate that also according to physiological criteria, SPL9 and SPL15 redundantly promote the juvenile-to-adult phase transition. In addition, other miR156-regulated SPL genes are expected to contribute as well based on the behaviour of the MIR156b overexpressor.
The role of gibberellin in the function of SPL9 and SPL15
The plant hormone gibberellin (GA) is known to promote flowering in many plants and in Arabidopsis it is particularly required for flowering in SD (Wilson et al. 1992). Exogenous application of GA will induce abaxial trichomes on leaves where these are normally not present although they will not appear earlier than leaf three (Telfer et al. 1997). In order to test whether the spl9 spl15 double mutant and the MIR156b overexpressor are defective in gibberellin sensitivity or biosynthesis, we exogenously applied GA3 and compared the onset of abaxial trichome production to mock treated and wild-type plants grown in LD1 (Fig. 5). Like in wild type, the GA3 treatment strongly reduced the number of rosette leaves without abaxial trichomes, i.e. juvenile leaves, with about a factor of three. This result shows that the spl9 spl15 double mutant and the MIR156b overexpressor remained sensitive to GA3. However, in both the spl9 spl15 double mutant and the MIR156b overexpressor the amount of GA3 applied could not reduce the number of juvenile leaves to that of obtained in GA3 treated wild type.
Discussion
Transgenic plants constitutively over-expressing the plant specific miRNA156 have been described for Arabidopsis (Schwab et al. 2005; Wu and Poethig 2006) and rice (Xie et al. 2006). Recently, overexpression of a miR156 encoding locus has also been shown to be the cause of the natural maize mutant Corngrass1 phenotype (Chuck et al. 2007). Interestingly, in all three different species, overexpression of this well conserved miRNA (Axtell and Bartel 2005; Arazi et al. 2005) causes a similar phenotype suggesting an evolutionary conserved role for the function of miR156 and its SPL target genes. Generally, in comparison to the respective wild type, miR156 over-expressing plants are smaller, flower later, tend to lose apical dominance and initiate more leaves with a shorter plastochron. MiR156 targets eleven SBP-box genes in Arabidopsis but the results presented here clearly show that already simultaneous silencing of the two likely paralogous target genes, SPL9 and SPL15, well approximate the miR156 over-expressing phenotype regarding the traits mentioned above. SPL9 and SPL15 thus act as important and functionally redundant transcription factors regulating diverse processes in shoot maturation and most likely in combination with other miR156 regulated SPL genes. In agreement with this latter statement is our observation that in addition to spl9-1 and spl15-1, mutation of a third miR156 controlled gene, SPL2 (At5g43270; T-DNA insertion line SALK_022235), results in triple mutant plants showing an even better approximation to the MIR156b overexpressor phenotype (Supplementary Fig. 2; despite the absence of detectable SPL2 transcript, single homozygous spl2-1 mutant plants lack an obvious mutant phenotype, data not shown). In greenhouse LD conditions we found the triple mutant to have produced on average 17.3 ± 2.1 (n = 16) rosette leaves in comparison to 15.6 ± 2.6 for the spl9 spl15 double mutant. Col-0 wild-type and the MIR156b overexpressor plants grown in parallel produced 12.9 ± 1.7 and 22.5 ± 3.5 rosette leaves before flowering, respectively.
SPL9 and SPL15 positively regulate the juvenile-to-adult growth phase transition
It became clear with the detailed analysis of Arabidopsis MIR156b overexpressors by Wu and Poethig (2006) as well as with the description of the Corngrass1 mutant in maize by Chuck et al. (2007), that one of the major phenotypic alterations in miR156 over-expressing plants is an extended juvenile growth phase. This suggests that one of the important functions of miR156 targeted SBP-box genes is to promote the vegetative phase change. In agreement with this observation, Wu and Poethig (2006) showed that overexpression of the miRNA156 regulated gene SPL3 and its likely paralogs leads to a greatly shortened juvenile phase in Arabidopsis. Based on morphological markers (abaxial trichomes) and physiological parameters (response to photo-inductive stimulus) we found that the spl9 spl15 double mutant exhibit a delayed vegetative phase transition and, therefore, conclude that both genes are very likely involved in the positive regulation of this developmental process in a redundant fashion. Most likely, because of this redundancy, photoperiodic induction of the single spl9 and spl15 mutants is not much affected. However, the effect on the appearance of abaxial trichomes as a marker for the juvenile-to-adult phase transition appears to be stronger in the spl9 mutant in comparison to spl15. This may be due to the fact that in shoot apical development expression of SPL9 starts to increase before that of SPL15 (Fig. 1b).
As SPL9 and SPL15 promote the juvenile-to-adult growth phase transition and thus competence to respond to photoperiodic induction of flowering, it is interesting to note that both SPL9 and SPL15 themselves are strongly upregulated in the shoot apex upon such induction (Schmid et al. 2003). An additional role for these genes in establishing inflorescence or floral meristem identity may thus be suggested.
SPL9 and SPL15 negatively regulate leaf initiation rate
Our data on the leaf initiation rate suggest that SPL9 and SPL15 act negatively on leaf initiation. Silencing of both genes leads to a shorter plastochron. Other miR156 regulated SPL genes may act similarly as the plastochron is even further shortened in the MIR156b overexpressor plants. A few mutants are known to cause a shortened plastochron and most of them simultaneously affect phyllotaxy. However, we found that shortening of the plastochron due to loss of SPL gene function is neither the cause nor the result of a changed spatial distribution of leaf primordia at the shoot apex. A shorter plastochron without an altered phyllotaxy has also been reported for two rice mutants, plastochron1 and -2 (pla1, -2; Itoh et al. 1998; Miyoshi et al. 2004; Kawakatsu et al. 2006). PLA1 encodes a cytochrome P450 protein, whereas PLA2 encodes a MEI2-like RNA binding protein. In both mutants the reduction in plastochron is accompanied by an increase in the size of the SAM and a higher rate of cell division. However, although the SAM of pla2 is actually smaller than that of pla1, it has a shorter plastochron. Furthermore, higher cell division activity associated with constitutive overexpression of CyclinD shortened the plastochron in tobacco without altering SAM size (Cockcroft et al. 2000). These observations suggest, as already noticed by Kawakatsu et al. (2006), that not SAM size but rather cell division rate is decisive in plastochron duration. Also our results may lend support to this hypothesis as both the spl9 spl15 double mutant and the MIR156b overexpressor exhibit a clearly shorter plastochron than wild-type plants but their SAM sizes differ only marginally. Therefore, it will be of interest to determine if SPL9, SPL15 and other miR156 regulated SPL genes control cell division rate in the SAM and, if so, in particular if their role is mediated through the phytohormone cytokinin. Not only is cytokinin a major positive regulator of cell proliferation and division in plants (Werner et al. 2001) it is also, in mutual dependence of auxin, a major determinant in the outgrowth of lateral shoots (Sachs and Thimann 1967; Chatfield et al. 2000). This latter aspect may explain the reduced apical dominance observed for the spl9 spl15 double mutant and the MIR156b overexpressor. Finally, mutants disrupting cytokinin signalling are known to result in reduced leaf initiation rates in addition to a smaller SAM and other effects (Nishimura et al. 2004; Higuchi et al. 2004).
Do SPL9 and SPL15 negatively regulate leaf maturation rate?
Based on their observations of the pla mutants in correlation to the expression of the respective genes in leaf primordia but not in the SAM, Kawakatsu et al. (2006) proposed that the rate of leaf maturation plays a significant role in regulating the rate of leaf initiation. In addition, these authors postulated a model in which the inhibitory effect of pre-existing leaf primordia on the initiation of the next leaf is lost as they mature. Similarly, a shortened plastochron in the spl9, spl15 mutants and the MIR156b overexpressor may also be due to precocious maturation of their leaves as suggested by our comparison of cross sections through successive leaf primordia of wild type and the MIR156b overexpressor. Even if the shape, i.e. circularity, of the primordia may not significantly differ after correction for an altered plastochron by comparing primordia based on age and not on serial sequence number, their cross sectional area seems to increase more rapidly in the MIR156b overexpressor. In particular, the cells surrounding the midvein in MIR156b overexpressor leaves appear to enlarge more rapidly.
SPL9 and SPL15 do not modulate the role of GA in the vegetative phase change
Exogenous application of GA3 has been found to accelerate abaxial trichome production in Arabidopsis suggesting that gibberellins function to regulate vegetative phase change (Telfer et al. 1997). These findings are also supported by mutant analysis. For example in spindly (spy) mutants, which undergo constitutive GA response, abaxial trichomes occur on leaves initiated significantly earlier than in wild type (Jakobsen and Olszewski 1993; Telfer et al. 1997). On the other hand, Telfer et al. (1997) found that mutants blocked in GA biosynthesis as well as GA insensitive mutants are significantly delayed in the appearance of abaxial trichomes. Loss-of-function mutants for the here examined SPL9 and SPL15 genes clearly delay the appearance of abaxial trichomes.
Our treatment of spl9 spl15 double mutant and 35S::MIR156b transgenic plants with high doses of GA3 showed that, like in wild type, their number of juvenile leaves can be reduced but not to numbers equal to those found for similarly treated wild type. In fact the ratios of juvenile to adult leaves of wild type and mutant phenotypes remain highly comparable to untreated plants. From this we conclude that the role SPL9, SPL15 and other miR156 controlled SPL genes play in the vegetative phase change, is unlikely to be GA mediated although a minor contribution to GA sensitivity can not be excluded.
Outlook
MiR156 targeted members of the SBP-box family of transcription factors in both mono- and dicots appear to play an important role as positive regulators of shoot maturation and of the vegetative to reproductive phase transition in particular. Both genetic factors, i.e. miR156 and SBP-box genes, have also been suggested to be major determinants in the transition from undifferentiated to differentiated embryogenic calli of rice (Luo et al. 2006). As the interaction between SBP-box genes and miR156 is of ancient origin in land plants (Arazi et al. 2005; Riese et al. 2007) it will be interesting to learn to what extent their molecular interplay is of importance to developmental phase transitions in plants in general.
Abbreviations
- DAS:
-
Days after sowing
- GA:
-
Gibberellic acid
- LD:
-
Long day
- miRNA:
-
microRNA
- qRT-PCR:
-
Quantitative real-time PCR
- SAM:
-
Shoot apical meristem
- SBP:
-
SQUAMOSA PROMOTER BINDING PROTEIN
- SD:
-
Short day
- SPL :
-
SQUAMOSA PROMOTER BINDING PROTEIN-LIKE
References
Alonso JM, Stepanova AN, Leisse TJ, Kim CJ, Chen H, Shinn P, Stevenson DK, Zimmerman J, Barajas P, Cheuk R, Gadrinab C, Heller C, Jeske A, Koesema E, Meyers CC, Parker H, Prednis L, Ansari Y, Choy N, Deen H, Geralt M, Hazari N, Hom E, Karnes M, Mulholland C, Ndubaku R, Schmidt I, Guzman P, Aguilar-Henonin L, Schmid M, Weigel D, Carter DE, Marchand T, Risseeuw E, Brogden D, Zeko A, Crosby WL, Berry CC, Ecker JR (2003) Genome-wide insertional mutagenesis of Arabidopsis thaliana. Science 301:653–657
Arazi T, Talmor-Neiman M, Stav R, Riese M, Huijser P, Baulcombe DC (2005) Cloning and charaterization of micro-RNAs from moss. Plant J 43:837–848
Axtell MJ, Bartel DP (2005) Antiquity of microRNAs and their targets in land plants. Plant Cell 17:1658–1673
Baumann E, Lewald J, Saedler H, Schulz B, Wisman E (1998) Successful PCR-based reverse genetic screens using an En-1 mutagenized Arabidopsis thaliana population generated via single-seed descendent. Theor Appl Genet 97:729–734
Bollmann KM, Aukerman MJ, Park MY, Hunter C, Berardini TZ, Poethig RS (2003) HASTY, the Arabidopsis ortholog of exportin 5/MSN5, regulates phase change and morphogenesis. Development 130:1493–1504
Cardon G, Höhmann S, Nettesheim K, Saedler H, Huijser P (1997) Functional analysis of the Arabidopsis thaliana SBP-bx gene SPL3: a novel gene involved in the floral transition. Plant J 12:367–377
Cardon G, Höhmann S, Klein J, Nettesheim K, Saedler H, Huijser P (1999) Molecular characterization of the Arabidopsis SBP-box genes. Gene 237:91–104
Chatfield SP, Stirnberg P, Forde B, Leyser O (2000) The hormonal regulation of axillary bud growth in Arabidopsis. Plant J 24:159–169
Chuck G, Cigan AM, Saeteurn K, Hake S (2007) The heterochronic maize mutant Corngrass1 results from overexpression of a tandem microRNA. Nat Genet 39:544–549
Clarke JH, Tack D, Findlay K, Van Montagu M, Van Lijsebettens M (1999) The SERRATE locus controls the formation of the early juvenile leaves and phase length in Arabidopsis. Plant J 20:493–501
Cockcroft CE, den Boer BGW, Healy JMS, Murray JAH (2000) CyclinD control of growth rate in plants. Nature 405:575–579
Czechowski T, Stitt M, Altmann T, Udvardi MK, Scheible W-R (2005) Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiol 139:5–17
Erickson RO, Michelini FJ (1957) The plastochron index. Am J Bot 44:297–305
Gandikota M, Birkenbihl RP, Höhmann S, Cardon G, Saedler H, Huijser P (2007) The miRNA156/157 recognition element in the 3′ UTR of the Arabidopsis SBP box gene SPL3 prevents early flowering by translational inhibition in seedlings. Plant J 49:683–693
Gasciolli V, Mallory AC, Bartel DP, Vaucheret H (2005) Partially redundant functions of Arabidopis DICER-like enzymes and a role for DCL4 in producing trans-acting siRNAs. Curr Biol 15:1494–1500
Higuchi M, Pischke MS, Mähönen AP, Miyawaki K, Hashimoto Y, Seki M, Kobayashi M, Shinozaki K, Kato T, Tabata S, Helariutta Y, Sussman MR, Kakimoto T (2004) In planta functions of the Arabidopsis cytokinin receptor family. Proc Natl Acad Sci USA 101:8821–8826
Hunter C, Sun H, Poethig RS (2003) The Arabidopsis heterochronic gene ZIPPY is an ARGONAUTE family member. Curr Biol 13:1734–1739
Itoh J-I, Hasegawa A, Kitano H, Nagato Y (1998) A recessive heterochronic mutation, plastochron1, shortens the plastochron and elongates the vegetative phase in rice. Plant Cell 10:1511–1521
Jakobsen SE, Olszewski NE (1993) Mutations at the SPINDLY locus of Arabidopsis alter gibberellin signal transduction. Plant Cell 5:887–896
Kawakatsu T, Itho J-I, Miyoshi K, Kurata N, Alvarez N, Veit B, Nagato Y (2006) PLASTOCHRON2 regulates leaf initiation and maturation in rice. Plant Cell 18:612–625
Komeda Y (2004) Genetic regulation of time to flower in Arabidopsis thaliana. Annu Rev Plant Biol 55:521–535
Li Y, Rosso MG, Viehoever P, Weisshaar B (2007) GABI-Kat SimpleSearch: an Arabidopsis thaliana T-DNA mutant database with detailed information for confirmed insertions. Nucleic Acids Res 35(Database issue):D874–D878
Luo Y-C, Zhou H, Li Y, Chen J-Y, Yang J-H, Chen Y-Q, Qu L-H (2006) Rice embryogenic calli express a unique set of microRNAs, suggesting regulatory roles of microRNAs in plant post-embryogenic development. FEBS Lett 580:5111–5116
Miyoshi K, Ahn BO, Kawakatsu T, Ito Y, Itoh Y, Nagato Y, Kurata N (2004) PLASTOCHRON1 a timekeeper of leaf initiation in rice, encodes cytochrome P450. Proc Natl Acad Sci USA 101:875–880
Nishimura C, Ohashi Y, Sato S, Kato T, Tabata S, Ueguchi C (2004) Genetic analysis of Arabidopsis histidine kinase genes encoding cytokinin receptors reveals their overlapping biological functions in the regulation of shoot and root growth in Arabidopsis thaliana. Plant Cell 16:1365–1377
Peragine A, Yoshikawa M, Wu G, Albrecht HL, Poethig RS (2004) SGS3 and SGS2/SDE1/RDR6 are required for juvenile development and the production of trans-acting siRNAs in Arabidopsis. Genes Develop 18:2368–2379
Poethig RS (1990) Phase change and the regulation of shoot morphogenesis in plants. Science 250:923–930
Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B, Bartel DP (2002) Prediction of plant MicroRNA targets. Cell 110:513–520
Riese M, Höhmann S, Saedler H, Münster T, Huijser P (2007) Comparative analysis of the SBP-box gene families in P. patens and seed plants. Gene 401:28–37
Sachs T, Thimann KV (1967) The role of auxins and cytokinins in the release of buds from dominance. Am J Bot 54:136–144
Sambrook J, Fritsch EF, Maniatis T (1989) Molecular cloning: a laboratory manual. Cold Spring Harbor Laboratory Press, Cold Spring Harbor
Schmid M, Uhlenhaut NH, Godard F, Demar M, Bressan R, Weigel D, Lohmann JU (2003) Dissection of floral induction pathways using global expression analysis. Development 130:6001–6012
Schmid M, Davison TS, Henz SR, Pape UJ, Demar M, Vingron M, Schölkopf B, Weigel D, Lohmann J (2005) A gene expression map of Arabidopsis development. Nat Genet 37:501–506
Schwab R, Palatnik JF, Riester M, Schommer C, Schmid M, Weigel D (2005) Specific effects of MicroRNAs on the plant transcriptome. Dev Cell 8:517–527
Telfer A, Poethig RS (1998) HASTY: a gene that regulates the timing of shoot maturation in Arabidopsis thaliana. Development 125:1889–1898
Telfer A, Bollman KM, Poethig RS (1997) Phase change and the regulation of trichome distribution in Arabidopsis thaliana. Development 124:645–654
Unte US (2001) Funktionelle Charakterisierung von SBP-Box-Genen mit Hilfe der Modellpflanze Arabidopsis thaliana. Dissertation, University of Cologne
Werner T, Motyka V, Strnad M, Schmülling T (2001) Regulation of plant growth by cytokinin. Proc Natl Acad Sci USA 98:10487–10492
Wilson RN, Heckman JW, Somerville CR (1992) Gibberellin is required for flowering in Arabidopsis thaliana under short days. Plant Physiol 100:403–408
Wu G, Poethig RS (2006) Temporal regulation of shoot development in Arabidopsis thaliana by miR156 and its target SPL3. Development 133:3539–3547
Xie Z, Allen E, Wilken A, Carrington JC (2005) DICER-LIKE 4 functions in trans-acting small interfering RNA biogenesis and vegetative phase change in Arabidopsis thaliana. Proc Natl Acad Sci USA 102:12984–12989
Xie K, Wu C, Xiong L (2006) Genomic organization, differential expression, and interaction of SQUAMOSA promoter-binding-like transcription factors and microRNA156 in rice. Plant Physiol 142:280–293
Yoshikawa M, Peragine A, Park MY, Poethig RS (2005) A pathway for the biogenesis of trans-acting siRNAs in Arabidopsis. Genes Develop 19:2164–2175
Acknowledgements
We are grateful to Dr Detlef Weigel and Dr Rebecca Schwab (MPI, Tübingen, Germany) for the kind gift of the 35S::MIR156b line and Dr Ulrike Unte for providing the spl9-1 allele. The Salk Institute Genomic Analysis Laboratory is acknowledged for providing the sequence-indexed Arabidopsis T-DNA insertion mutants. We thank Susanne Höhmann for excellent technical support.
Open Access
This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
Author information
Authors and Affiliations
Corresponding author
Electronic supplementary material
Below is the link to the electronic supplementary material.
Supplementary Fig. 1
Expression of SPL15 as detected with a GUS reporter gene construct. A pSPL15::GUS:SPL15 reporter was constructed by subcloning a 2,793 bp genomic fragment of SPL15 beginning 1,260 bp upstream of the SPL15-ATG start codon and extending to 132 bp downstream of the SPL15-stop codon into the binary vector pGJ2148 (kindly provided by Guido Jach, MPIZ Cologne). A ß-glucuronidase (GUS) reporter was subsequently cloned in frame shortly downstream of the SPL15-ATG start codon and the whole construct stably transformed into Arabidopsis Col-0 background. (A–B) GUS staining of mature embryo gently squeezed out of the imbibed seed (A) and of a ca. 10 day old seedling grown in LD (B). Bars in A and B represent respectively 200 μm and 2 mm. (TIF 26,178 KB)
Supplementary Fig. 2
Phenotypes of spl9 and spl15 related mutants. (A) Flowering plants homozygous for spl9-1, -2 and -3 alleles in Col-0 background next to wild type. (B) Flowering plants homozygous for spl15-1, -2 and -3 alleles in Col-0 background next to wild type. (C) Flowering plants, from left to right, of Col-0 wild type, the double mutants spl9-1 spl15-1 and spl9-2 spl15-2, the triple mutant spl2-1 spl9-1 spl15-1 and the MIR156b overexpressor. (TIF 32,914 KB)
Rights and permissions
Open Access This is an open access article distributed under the terms of the Creative Commons Attribution Noncommercial License (https://creativecommons.org/licenses/by-nc/2.0), which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
About this article
Cite this article
Schwarz, S., Grande, A.V., Bujdoso, N. et al. The microRNA regulated SBP-box genes SPL9 and SPL15 control shoot maturation in Arabidopsis. Plant Mol Biol 67, 183–195 (2008). https://doi.org/10.1007/s11103-008-9310-z
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11103-008-9310-z