Abstract
Leaves develop as flat lateral organs from the indeterminate shoot apical meristem. The establishment of polarity along three-dimensional axes, proximal–distal, medial–lateral, and adaxial–abaxial axes, is crucial for the growth of normal leaves. The mutations of ASYMMETRIC LEAVES1 (AS1) and AS2 of Arabidopsis thaliana cause defects in repression of the indeterminate state and the establishment of axis formation in leaves. Although many mutations have been identified that enhance the adaxial–abaxial polarity defects of as1 and as2 mutants, the roles of the causative genes in leaf development are still unknown. In this study, we found that wild-type plants treated with berberine produced pointed leaves, which are often observed in the single mutants that enhance phenotypes of as1 and as2 mutants. The berberine-treated as1 and as2 mutants formed abaxialized filamentous leaves. Berberine, an isoquinoline alkaloid compound naturally produced in various plant sources, has a growth inhibitory effect on plants that do not produce berberine. We further showed that transcript levels of meristem-specific class 1 KNOX homeobox genes and abaxial determinant genes were increased in berberine-treated as1 and as2. Berberine treated plants carrying double mutations of AS2 and the large subunit ribosomal protein gene RPL5B showed more severe defects in polarity than did the as2 single mutant plants. We suggest that berberine inhibits (a) factor(s) that might be required for leaf adaxial cell differentiation through a pathway independent of AS1 and AS2. Multiple pathways might play important roles in the formation of flat symmetric leaves.
Similar content being viewed by others
Introduction
Leaves develop as flat lateral organs from the shoot apical meristem, and the establishment of polarity along three-dimensional axes, proximal–distal, medial–lateral and adaxial–abaxial, is crucial for the growth of normal leaves (Steeves and Sussex 1989; Hudson 2000; Tsukaya 2006; Bowman and Floyd 2008; Szakonyi et al. 2010). Thus plants appear to exploit common mechanisms that are responsible for the establishment of these axes during leaf development. Genes that are involved in the formation of each of the three axes have been identified and their functions have been analyzed (Waites et al. 1998; Bowman and Smyth 1999; Sawa et al. 1999; Byrne et al. 2000; Nemhauser et al. 2000; Eshed et al. 2001; Kerstetter et al. 2001; McConnell et al. 2001; Iwakawa et al. 2002; Pekker et al. 2005; Garcia et al. 2006). The asymmetric leaves1 (as1) and as2 mutants of Arabidopsis thaliana (A. thaliana) produce leaves abnormal in axis formation, and these mutants have pleiotropic phenotypes (Rédei and Hirono 1964; Tsukaya and Uchimiya 1997; Serrano-Cartagena et al. 1999; Byrne et al. 2000, 2001, 2002; Ori et al. 2000; Semiarti et al. 2001; Iwakawa et al. 2002, 2007; Matsumura et al. 2009; Ikezaki et al. 2010). Analyses by DNA microarray and real-time quantitative RT-PCR (qRT-PCR) revealed that a mutation of either AS1 or AS2 results in an increase in the transcription levels of genes for the abaxial determinants, ETTIN/AUXIN RESPONSE FACTOR3 (ETT/ARF3), KANADI2 (KAN2) and YABBY5 (YAB5) as well as the class 1 KNOX genes, BREVIPEDICELLUS (BP), KNAT2, and KNAT6 (Iwakawa et al. 2007; Takahashi et al. 2008). The class 1 KNOX genes (BP, KNAT2, KNAT6, and SHOOT MERISTEMLESS [STM]) are expressed in the shoot apical meristem and thought to function in maintaining the indeterminate cell state (Lincoln et al. 1994; Long et al. 1996; Byrne et al. 2000; Ori et al. 2000; Semiarti et al. 2001; Ikezaki et al. 2010). AS1 encodes a myb domain protein and AS2 encodes a plant-specific protein with an AS2/LOB domain, both of which are nuclear proteins (Byrne et al. 2000; Iwakawa et al. 2002). These proteins form a complex structure (Yang et al. 2008) (hereinafter referred to as AS1–AS2) that is involved in the establishment of proper axis formation. In addition, some mutations of genes alter leaf patterning and axis formation in the absence of either AS1 or AS2: namely, ribosomal protein genes, genes involved in chromatin remodeling, histone modification, and small RNA biogenesis, and others (Ori et al. 2000; Xu et al. 2003, 2006; Kidner and Martienssen 2005; Li et al. 2005; Garcia et al. 2006; Huang et al. 2006; Yang et al. 2006; Ueno et al. 2007; Pinon et al. 2008; Yao et al. 2008; Pérez-Pérez et al. 2010; Horiguchi et al. 2011a, b; Kojima et al. 2011; Szakonyi and Byrne 2011; Ishibashi et al. 2012). Although the responsible wild-type genes also have important roles in the maintenance of the plants themselves, the roles of these wild-type genes in the establishment of leaf polarity are still unclear.
Alkaloids are nitrogen-containing low molecular weight substances that form a major class of plant-produced secondary metabolites. We screened for natural compounds that affect leaf development in A. thaliana and found that berberine, a benzylisoquinoline alkaloid, caused a pointed-leaf phenotype in wild-type plants. Berberine is naturally produced in various plant sources, such as the wood of Berberis, inner bark of Phellodendron, and rhizomes of Coptis (Grycová et al. 2007; Shitan and Yazaki 2007). The inhibitory effect of berberine on growth was reported in plant species that cannot produce berberine (Sakai et al. 2002). It has been reported that Arabidopsis plants might not produce berberine (Facchini et al. 2004); nevertheless, effects of berberine on A. thaliana have not been reported. For a half century, it has been known that berberine interacts with DNA and RNA (Klímek and Hnilica 1959; Yamagishi 1962; Kley and Hahn 1969; Wolfe et al. 1972). Berberine is used in biochemical research as a fluorescent probe of DNA and RNA (Hirakawa et al. 2005; Islam and Kumar 2008; Islam et al. 2009). Recently, it has been reported that berberine preferentially binds to triplex or quadruplex DNA (Lee et al. 1993; Ren and Chaires 1999; Franceschin et al. 2006; Arora et al. 2008; Ma et al. 2008; Bhadra and Kumar 2011). Guanine-rich sequences that presumably form quadruplex DNA are found throughout the genomes studied, including that of A. thaliana (Huppert and Balasubramanian 2005; Todd et al. 2005; Zhao et al. 2007; Mullen et al. 2010; Takahashi et al. 2012). It has been reported that sequences containing quadruplex DNA affect the transcript levels of genes located nearby (Siddiqui-Jain et al. 2002). In recent decades, berberine has been reported to possess antibacterial (Chae et al. 1999; Hu et al. 2000; Yan et al. 2008), anti-cancer (Anis et al. 2001), and antiphotooxidative activities (Kim et al. 2000). Berberine also induces apoptosis in several types of cells (Kuo et al. 1995; Hwang et al. 2006; Mantena et al. 2006) and has been used in medicines for both traditional and clinical applications. In spite of these reports and the broad uses of berberine, there have been no reports on the cause of its effects on growth inhibition and leaf morphology in plants that do not produce berberine.
In the present study, we report that wild-type Arabidopsis plants treated with berberine showed growth inhibition and changes of leaf morphology. Furthermore, we detected severe defects of leaf phenotypes related to the adaxial–abaxial polarity in as1 and as2 mutants treated with berberine. We discuss these effects of berberine that interfered with normal leaf development in A. thaliana.
Results
Phenotypes of wild-type plants treated with berberine
When we treated plants with berberine (Fig. 1Aa), they displayed obviously severe effects of growth inhibition with decreased leaf size (Fig. 1Ab), similar to those observed in plants treated with trichostatin A (Ueno et al. 2007). As shown in Fig. 1B, the sizes of the first and second rosette leaves in plants at 18 days after sowing with 20 μM berberine were comparable to those at 12 days after sowing without berberine. Furthermore, the venation pattern of berberine-treated leaves was simpler than that of untreated leaves (Fig. 1B). We also observed pointed leaves at 20 and 25 μM concentrations of berberine similar to those observed in oligocellula5-1 (oli5-1) and oligocellula7-1 (oli7-1) mutant plants, each of which has a mutation in the large subunit of ribosomal protein gene RPL5A and RPL5B, respectively (Horiguchi et al. 2006; Fujikura et al. 2009; Fig. 1Cb, D; Table 1). These oli5-1 and oli7-1 mutants were identified as showing decreased palisade leaf cell numbers and pointed-leaf phenotypes. The leaf polarity defects of oli5-1 were rather more severe than those of oli7-1 (Horiguchi et al. 2011a). We quantified the pointed-leaf phenotype by measuring the tip angles of the first and second leaves. For wild-type plants grown without berberine, the mean tip angle was 123.7°, whereas those for berberine-treated wild-type and the oli5-1 and oli7-1 mutant plants were 118.7°, 98.6°, and 104.1°, respectively (p < 0.01). In this paper, we categorized leaves with tip angles less than 110° as ‘pointed leaves’. As shown in Fig. 1Db, there were no leaves with tip angles less than 110° seen in wild-type plants without berberine; whereas, 14 % of wild-type leaves treated with berberine (6 out of 42 leaves), 92 % of oli5-1 leaves (34 out of 37 leaves), and 69 % oli7-1 leaves (34 out of 42 leaves) had leaf-tip angles smaller than 110°.
Phenotypes of as1-1 and as2-1 plants treated with berberine
Mutations that cause pointed leaves, such as oli5-1 and oli7-1, have been previously reported to enhance the phenotypes of as1 and as2 mutants (Fig. 1Da; Pinon et al. 2008; Yao et al. 2008; Horiguchi et al. 2011a; Szakonyi and Byrne 2011; Ishibashi et al. 2012). Therefore, we examined the phenotypes of wild-type and as1-1 and as2-1 mutant plants treated without and with 12.5, 20, and 25 μM berberine. The results are given in Table 1. The pointed-leaf phenotype was found in 7, 8, and 7 % of wild-type plants as well as in as1-1 and as2-1 plants treated with 12.5, 20, and 25 μM berberine, respectively (Table 1). Wild-type plants seldom showed radialized leaves (Table 1), whereas 3, 13, and 18 % of as1-1 and 8, 35, and 28 % of as2-1 plants had trumpet-like leaves and filamentous leaves indicative of a defect in adaxial–abaxial polarity at 12.5, 20, and 25 μM berberine, respectively (Fig. 2Ac, Ad, Ae, Af; Table 1). We also observed plants without true leaves. When grown with 25 μM berberine, true leaves were not seen on 17 % of wild type, 27 % of as1-1, and 34 % of as2-1 plants (Table 1). These results indicated that as1-1 and as2-1 mutant plants were more sensitive than the wild-type plants to berberine treatment.
Since trumpet-like and filamentous leaves are often generated in plants having defective adaxial–abaxial polarity determination, we examined whether these leaf phenotypes are associated with transcript levels of certain genes in the presence of berberine treatment. To examine the effect of berberine treatment on leaf adaxial–abaxial polarity, we used green fluorescent protein (GFP) under the control of the FILAMENTOUS FLOWER (FIL) promoter (FILp:GFP), which is expressed in abaxial cells of leaf primordia (Watanabe and Okada 2003). Strong signals were detected in cells located at peripheral positions in filamentous leaves in the presence of 20 μM berberine (Fig. 2Bc, Bd), whereas no signals were detected in cells located on the adaxial side of leaves of as1-1 and as2-1 plants in the absence of berberine (Fig. 2Ba, Bb).
We next examined the accumulated transcript levels of genes involved in leaf adaxial–abaxial polarity and class 1 KNOX family genes in 14-day-old shoots treated with 12.5 μM berberine. Figure 2C shows the changes in the accumulated transcript levels of representative genes. Supplementary Table S2 shows relative changes in transcript levels of all the genes analyzed in each plant. When the wild-type plants were treated with 12.5 μM berberine, the accumulated transcript levels of BP and KAN4 showed significant increases over twofold higher than those levels of the non-treated wild type (Fig. 2C). Accumulated transcript levels of BP, KAN4, ETT/ARF3, and ARF4 in either the as1-1 or as2-1 or both plants treated with berberine at 12.5 μM showed significant increases over twofold higher than those of non-treated as1-1 or as2-1 plants (Fig. 2C; Supplementary Table S2). In contrast, the accumulated transcript levels of the adaxial determinant HD-ZIPIII genes PHABULOSA (PHB), PHAVOLUTA (PHV), and REVOLUTA (REV) were not significantly changed by berberine treatment in wild-type, as1-1, and as2-1 plants (Supplementary Table S2).
Phenotypes of plants treated with other isoquinoline alkaloid compounds
To understand which structure of the isoquinoline alkaloid compounds is responsible for their function, we next examined whether two other isoquinoline alkaloid compounds, coptisine and palmatine, having chemical structures similar to that of berberine (Fig. 3Aa, Ba), affect the phenotypes of as1 and as2 leaves. The sensitivity of plants to coptisine was similar to that of berberine (Fig. 3Ab). Trumpet-like and filamentous leaves were observed at concentrations from 20 to 50 μM in as1-1 and as2-1 plants (Fig. 3Ce, 3Cf), whereas these abnormal leaves were not observed in wild-type plants (Fig. 3Cd; Table 2). Trumpet-like and filamentous leaves were observed in 28 % (27 out of 96 plants) and 43 % (42 out of 96 plants) of as1-1 and as2-1 plants at 40 μM coptisine, respectively (Table 2). Plants without true leaves were not observed with coptisine treatment, although growth was inhibited. In contrast, treatments with 100 μM palmatine had no effect on either growth or leaf morphology (Fig. 3Bb, Cg, Ch, Ci). These results suggested that some specific conformation or residues might be important for affecting plant growth and morphology.
Phenotypes of plants treated with a DNA intercalator
Since it has been reported that berberine forms complexes with duplex, triplex, and quadruplex DNA (Klímek and Hnilica 1959; Yamagishi 1962; Kley and Hahn 1969; Lee et al. 1993; Ren and Chaires 1999; Franceschin et al. 2006), we attempted to treat as1 and as2 mutant plants with one of the DNA intercalators, N-methylmesoporphyrin IX (NMM), which has been known as a major quadruplex DNA ligand (Arthanari et al. 1998). NMM treatment severely inhibited growth of plants at 10 μM and plants rarely germinated at 25 μM (Fig. 4Ab). NMM treatment strongly inhibited true leaf formation. As shown in Fig. 4C, the sensitivity to NMM was stronger in as1-1 and as2-1 plants than in wild-type plants: 73 % of as1-1 plants and 75 % of as2-1 plants were deficient in true leaf formation, whereas 61 % of the wild-type plants were deficient at 3 μM NMM. Plants treated with NMM showed filamentous leaves in all wild-type, as1-1, and as2-1 plants at low concentrations of NMM; although we did not detect a significant difference in the frequencies of filamentous leaf formation between wild-type and the as1-1 or as2-1 mutant plants (Fig. 4Bd, Be, Bf).
Increased sensitivity of berberine by a mutation in cytosolic ribosomal protein RPL5B
Most DNA intercalators are also known to bind to ribosomal RNA, a component of the ribosome complex, and result in ribosome degradation (Wolfe et al. 1972). We treated wild-type, oli7-1, and oli7-1 as2-1 double mutant plants with berberine. Berberine treatment severely inhibited the growth of oli7-1 plants at 15 and 20 μM (Fig. 5A). When we treated these plants with berberine, 4 and 18 % of oli7-1 plants did not have true leaves at 15 and 20 μM berberine, respectively, whereas very lower frequencies (1 and 2 %) of wild-type plants appeared without true leaves at 15 and 20 μM berberine (Fig. 5B). The oli7-1 plants treated with 20 μM berberine also displayed the narrower leaf phenotype (orange arrowheads in Fig. 5Ag). When treated with 15 and 20 μM berberine, the narrow-leaf phenotype was found in 4 and 1 % of oli7-1 plants, 34 and 29 % of oli7-1 as2-1 double mutant plants, and 1 and 3 % of wild-type plants, respectively (Fig. 5A, B). Trumpet-like and filamentous leaves were also observed at 15 and 20 μM berberine treatment in 5 and 1 % of oli7-1 plants, respectively, which are comparable frequencies to those of wild-type plants. As reported previously (Horiguchi et al. 2011a), 17 % of the oli7-1 as2-1 double mutants displayed a filamentous-leaves phenotype in the absence of berberine, and the frequencies of the filamentous-leaves phenotype were approximately twofold increased (37 %) in the presence of both 15 and 20 μM berberine.
The effects of berberine on transcript levels in wild-type plants
To know how berberine affects leaf cells in wild-type plants, we used Affymetrix ATH1 microarrays to analyze transcript levels of genes in Arabidopsis plants grown on MS medium with and without 12.5 μM berberine for 14 days. Microarray analysis showed that transcript levels of 3,854 out of 22,591 genes on the ATH1 gene chip (excluding the mitochondrial genome) were changed by berberine treatment (Supplementary Table S3). To date, mutations of many ribosomal protein genes and other genes have been identified that enhance as1 and or as2 leaf phenotypes (Ori et al. 2000; Xu et al. 2003, 2006; Kidner and Martienssen 2005; Li et al. 2005; Garcia et al. 2006; Huang et al. 2006; Yang et al. 2006; Ueno et al. 2007; Pinon et al. 2008; Yao et al. 2008; Pérez-Pérez et al. 2010; Horiguchi et al. 2011a, b; Kojima et al. 2011; Szakonyi and Byrne 2011; Ishibashi et al. 2012). We analyzed whether the expression of these known enhancer genes was repressed by berberine treatment; however, transcript levels of such genes were not changed significantly by berberine.
Discussion
In the present study, we found that berberine affected the establishment of leaf polarity in the as1 or as2 mutant plants of A. thaliana. Our results by treatments with low molecular weight compounds implicated the existence of a novel factor(s), which is involved together with AS1 and AS2 in the establishment of adaxial cells of leaves.
Berberine affects leaf development of Arabidopsis and upregulates the transcript levels of class 1 KNOX and leaf abaxial determinant genes
When we treated the as1-1 and as2-1 mutant plants with berberine, they produced abaxialized filamentous leaves (Fig. 2A). Berberine might inhibit a function(s) of a novel factor(s) that is involved in adaxialization independently from AS1–AS2 during leaf development. Namely, the factor(s) and AS1–AS2 might share the same target to control the adaxialization of leaves through independent pathways (Fig. 6). It has been reported that the AS1–AS2 complex directly represses the transcription of BP and KNAT2 genes (Guo et al. 2008). Transcription of the abaxial determinant genes ETT/ARF3, KAN2, and YAB5 was also repressed by AS1 and AS2 (Iwakawa et al. 2007). Our real-time qRT-PCR results also showed that berberine treatment increased transcript levels of BP and KAN4 genes in wild-type plants, and the class 1 KNOX genes BP and KNAT6 and the abaxial determinant genes KAN4, ETT/ARF3, and ARF4 in the as1 and as2 mutants. These results suggest that the novel factor(s) inhibited by berberine somehow represses BP and KAN4 genes through an unidentified pathway that would act in parallel with AS1–AS2 for the establishment of leaf polarity, resulting in the formation of flat symmetric leaves in A. thaliana (Fig. 6).
Factors inhibited by berberine and coptisine that are required for the establishment of leaf polarity in parallel with AS1 and AS2
Berberine, coptisine, and palmatine are close structural analogs of each other. They bear the same tetracyclic structure (7,8,13,13a-tetrahydroberbinium), but differ in the nature of the substituents at positions 2 and 3 on the benzo ring, being methylene dioxy for berberine and coptisine, and dimethoxy for palmatine. Similarly to berberine, coptisine and palmatine also have a buckled structure. Recently, berberine and coptisine have been reported to inhibit the receptor activator of nuclear factor-kappaB ligand (RANKL)-induced osteoclast differentiation from osteoclast precursors, the survival of mature osteoclasts, and ovariectomy-stimulated bone resorption in rats (Li et al. 1999; Hu et al. 2008; Lee et al. 2008; Lee et al. 2012). In contrast, the inhibitory effect of palmatine on osteoclast differentiation and function is markedly weaker than that of berberine (Hu et al. 2008; Lee et al. 2010, 2012). Since both berberine and coptisine affected leaf development, the methylene dioxy substituent at positions 2 and 3 on the benzo ring might be responsible for determining their targets, although their targets might be different from that in RANKL-induced osteoclast differentiation in rats.
Multiple pathways are involved in the adaxial development of as1 and as2 leaves
Our microarray analysis showed that berberine negatively affected the expression of many genes; however, we did not detect decreased transcript levels of those genes that have been known to interact genetically with AS1 and/or AS2. To seek for the means by which berberine repressed such unknown factor(s), we considered several possible mechanisms of repression (inhibition), such as via berberine intercalation with secondary structures of DNA (quadruplex and triplex) and RNA, or by binding with proteins (Meyerson et al. 1978). Berberine treated plants having double mutations of AS2 and RPL5B showed more severe defects in polarity than did those as2 single mutant plants (Fig. 5A, B). These results suggest that berberine might repress other ribosomal components in the same pathway with RPL5B/OLI7/ANG3. Our genetic analysis also suggests that a process mediated by an unknown factor(s), whose function is inhibited by berberine, is the distinct process involving RPL5B for the establishment of leaf polarity. As we previously suggested in Kojima et al. (2011), the adaxial development of leaves might be controlled independently by multiple pathways that include the products of ELO3, RDR6-AGO7, and HDT1/HDT2. Considering that plants treated with another G-quadruplex ligand NMM also showed filamentous leaves (Fig. 4B), it might also be possible that berberine represses unknown factors via intercalating secondary DNA/RNA structures such as G-quadruplex DNAs. Especially, a DNA sequence that putatively forms a G-quadruplex is found in the transcription initiation site of ribosomal DNA in both the human (Hanakahi et al. 1999) and Arabidopsis genomes (Takahashi et al. 2012). It remains to be analyzed whether berberine binds to the rDNA region or to rRNA with the resulting impaired formation of leaf polarity. Taken together, the adaxial development of leaves is likely to be controlled independently by multiple pathways on the as1 or as2 background.
Since adaxial defects of leaves caused by berberine are not obvious in wild-type plants (Fig. 2C), the wild-type AS1–AS2 function can replace the polarity establishment function of an unknown factor(s) that is inhibited by berberine. That is to say, the wild-type AS1 and AS2 genes appear to attenuate deficiencies resulting from treatments with berberine and coptisine. Further investigations of the molecular functions of AS1 and AS2 proteins will allow us to understand how they could regulate gene expression by the presently unknown mechanism(s).
Our experiments are representative of a new method to identify novel factors that take part in leaf development by using the as1 or as2 mutant. The combination of comprehensive chemical screening and genetic analysis may make it possible to elucidate new events involved in the establishment of leaf polarity.
Materials and methods
Plant strains and growth conditions
Arabidopsis thaliana ecotype Col-0 (CS1092) and the mutants as1-1 (CS3374) and as2-1 (CS3117) were obtained from the Arabidopsis Biological Resource Center (ABRC) (Columbus, OH, USA). We outcrossed as2-1 with Col-0 three times and as1-1 with Col-0 once and used the progeny for our experiments (Semiarti et al. 2001). The oli mutants, back-crossed to Col-0 more than three times, were used for our experiments (Horiguchi et al. 2006; Fujikura et al. 2009). For the analysis of plants, seeds were surface-sterilized with 0.25 % sodium hypochlorite/0.1 % Triton X-100 solution and sown on Murashige and Skoog (MS) medium (Wako Pure Chemicals, Osaka, Japan) supplemented with 2 % (w/v) sucrose and 0.8 % agar. After 2 days at 4 °C in darkness, plants were transferred to a daily regimen of 8 h of darkness and 16 h of white light at 50 μmol m−2 S−1 at 22 °C, as described previously (Semiarti et al. 2001). Berberine chloride and palmatine chloride were purchased from Sigma-Aldrich (St. Louis, Missouri, USA). Coptisine chloride was purchased from Wako Pure Chemicals. NMM was purchased from Frontier Scientific, Inc. (Logan, Utah, USA). For chemical treatments, berberine, coptisine, or palmatine dissolved in dimethyl sulfoxide, or NMM dissolved in H2O was mixed with MS agar medium at 1–100 μM concentrations and immediately dispensed into plastic multi-dishes (multidish 267061, Thermo Fisher Scientific K. K., Yokohama, Japan). Ages of plants are given in terms of the number of days after sowing.
Morphological observation of leaves
Whole leaves were observed under a stereoscopic microscope (Stemi SV11, Carl Zeiss Inc., Oberkochen, Germany). Leaves were fixed in a formalin-acetic acid-alcohol solution and cleared in a chloral hydrate solution as previously described (Tsuge et al. 1996).
Quantification of angles of leaf tips
First or second true leaves were dissected from seedlings treated with and without berberine, mounted on a glass slide, and photographed under a stereoscopic microscope (Stemi SV11, Carl Zeiss). We drew a line along the mid vein from the one-eighth distal area of a leaf blade and a vertical line on leaves of the digital image using Adobe Illustrator software (Adobe Systems Incorporated, San Jose, CA, USA). We then drew lines toward the leaf tips from the crossing points with the vertical line and the outline of the leaf blade. The resulting angles were calculated by using the software ImageJ (http://rsbweb.nih.gov/ij/). Statistical significance was calculated by Student’s T test and F test by using Excel (Microsoft Corporation, Tokyo, Japan).
Histological analyses
Thin sections and cleared specimens were prepared as described previously (Semiarti et al. 2001). Leaf sections were cut by hand, and GFP fluorescence images were obtained with a confocal laser microscope (LSM510, Carl Zeiss Inc., Oberkochen, Germany) according to Ueno et al. (2007).
Real time quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR)
Shoots were harvested 14 days after seeds were sown and immediately frozen in liquid nitrogen and stored at −80 °C. Total RNA was isolated from 14-day-old seedlings with the RNeasy Plant Mini Kit (Qiagen, Valencia, CA) according to the manufacturer’s instructions. Sample volumes were normalized as described by Iwakawa et al. (2007). The primer pairs are shown in Supplementary Table S1. For the analysis of RNA levels in Arabidopsis by real-time qRT-PCR, we prepared 10 μg of total RNA. Reverse transcription was carried out with ReverTra Ace (TOYOBO, Osaka, Japan). PCR was performed in the presence of the double-stranded DNA-specific dye SYBR Green (Applied Biosystems, Lincoln, CA). Amplification was monitored in real time with the Applied Biosystems StepOnePlus Real-Time PCR system (Applied Biosystems) according to the supplier’s recommendations. The mean value of three technical replicates was normalized by using the TUBULIN α-3 (TUA3) transcript as a control. Statistical significance was calculated by Student’s T test and F test by using Excel.
Microarray
Seedlings of wild-type plants grown with and without 12.5 μM berberine for 14 days were harvested and total RNA was extracted as described in the ‘real time qRT-PCR’ section, above. For microarray analysis, the quality and purity of the RNA were confirmed with an Ultrospec 2100 pro (GE Healthcare UK Ltd, Amersham Place, England). Total RNA samples (8 μg) were reverse-transcribed, yielding double-stranded cDNA, which was transcribed in vitro in the presence of biotin-labeled nucleotides with an IVT Labeling Kit (Affymetrix Inc., Santa Clara, CA), and purified. Labeled cRNA was fragmented and hybridized to Affymetrix ATH1 GeneChip arrays for 16 h at 45 °C according to Affymetrix protocols. Arrays were washed on an Affymetrix Fluidics Station 450 and measured for fluorescence intensity with an Affymetrix GeneChip Scanner 3000. The raw data were processed by using Affymetrix Gene Chip Operating Software (GCOS; Version 1.4.0.036).
Microarray analysis
We calculated the expression signals for two samples and the log2 ratio for berberine-treated samples against the control sample by GCOS. The detection calls and change calls were calculated based on the signal-intensity ratio between perfect-match and miss-match oligos on the array. In this experiment, we excluded those genes (6,782 genes) for which two samples showed an absent or marginal call (i.e., detection call determined by GCOS, based on the p value of the one-sided Wilcoxon signed-rank test: absent call means p ≥ 0.065 and marginal call means 0.065 > p ≥ 0.05), because absent indicates that the expression signal was undetectable. We also excluded those genes (18,737 genes) for which two samples showed no change calls (i.e., change call determined by GCOS: no change call means (1–0.006) ≥ p ≥ 0.006 for the one-sided Wilcoxon signed-rank test), because no change indicates that the expression signal is almost equal to that of the control sample.
References
Anis KV, Rajeshkumar NV, Kuttan R (2001) Inhibition of chemical carcinogenesis by berberine in rats and mice. J Pharm Pharmacol 53:763–768
Arora A, Balasubramanian C, Kumar N, Agrawal S, Ojha RP, Maiti S (2008) Binding of berberine to human telomeric quadruplex—spectroscopic, calorimetric and molecular modeling studies. FEBS J 275:3971–3983
Arthanari H, Basu S, Kawano TL, Bolton PH (1998) Fluorescent dyes specific for quadruplex DNA. Nucleic Acids Res 26:3724–3728
Bhadra K, Kumar GS (2011) Interaction of berberine, palmatine, coralyne, and sanguinarine to quadruplex DNA: a comparative spectroscopic and calorimetric study. Biochim Biophys Acta 1810:485–496
Bowman JL, Floyd SK (2008) Patterning and polarity in seed plant shoots. Annu Rev Plant Biol 59:67–88
Bowman JL, Smyth DR (1999) CRABS CLAW, a gene that regulates carpel and nectary development in Arabidopsis, encodes a novel protein with zinc finger and helix-loop-helix domains. Development 126:2387–2396
Byrne ME, Barley R, Curtis M, Arroyo JM, Dunham M, Hudson A, Martienssen RA (2000) Asymmetric leaves1 mediates leaf patterning and stem cell function in Arabidopsis. Nature 408:967–971
Byrne M, Timmermans M, Kidner C, Martienssen R (2001) Development of leaf shape. Curr Opin Plant Biol 4:38–43
Byrne ME, Simorowski J, Martienssen RA (2002) ASYMMETRIC LEAVES1 reveals knox gene redundancy in Arabidopsis. Development 129:1957–1965
Chae SH, Jeong IH, Choi DH, Oh JW, Ahn YJ (1999) Growth-inhibiting effects of Coptis japonica root-derived isoquinoline alkaloids on human intestinal bacteria. J Agric Food Chem 47:934–938
Eshed Y, Baum SF, Perea JV, Bowman JL (2001) Establishment of polarity in lateral organs of plants. Curr Biol 11:1251–1260
Facchini PJ, Bird DA, St-Pierre B (2004) Can Arabidopsis make complex alkaloids? Trends Plant Sci 9:116–122
Franceschin M, Rossetti L, D’Ambrosio A, Schirripa S, Bianco A, Ortaggi G, Savino M, Schultes C, Neidle S (2006) Natural and synthetic G-quadruplex interactive berberine derivatives. Bioorg Med Chem Lett 16:1707–1711
Fujikura U, Horiguchi G, Ponce MR, Micol JL, Tsukaya H (2009) Coordination of cell proliferation and cell expansion mediated by ribosome-related processes in the leaves of Arabidopsis thaliana. Plant J 59:499–508
Garcia D, Collier SA, Byrne ME, Martienssen RA (2006) Specification of leaf polarity in Arabidopsis via the trans-acting siRNA pathway. Curr Biol 16:933–938
Grycová L, Dostál J, Marek R (2007) Quaternary protoberberine alkaloids. Phytochemistry 68:150–175
Guo M, Thomas J, Collins G, Timmermans MC (2008) Direct repression of KNOX loci by the ASYMMETRIC LEAVES1 complex of Arabidopsis. Plant Cell 20:48–58
Hanakahi LA, Sun H, Maizels N (1999) High affinity interactions of nucleolin with G–G-paired rDNA. J Biol Chem 274:15908–15912
Hirakawa K, Kawanishi S, Hirano T (2005) The mechanism of guanine specific photooxidation in the presence of berberine and palmatine: activation of photosensitized singlet oxygen generation through DNA-binding interaction. Chem Res Toxicol 18:1545–1552
Horiguchi G, Ferjani A, Fujikura U, Tsukaya H (2006) Coordination of cell proliferation and cell expansion in the control of leaf size in Arabidopsis thaliana. J Plant Res 119:37–42
Horiguchi G, Mollá-Morales A, Pérez-Pérez JM, Kojima K, Robles P, Ponce MR, Micol JL, Tsukaya H (2011a) Differential contributions of ribosomal protein genes to Arabidopsis thaliana leaf development. Plant J 65:724–736
Horiguchi G, Nakayama H, Ishikawa N, Kubo M, Demura T, Fukuda H, Tsukaya H (2011b) ANGUSTIFOLIA3 plays roles in adaxial/abaxial patterning and growth in leaf morphogenesis. Plant Cell Physiol 52:112–124
Hu JP, Takahashi N, Yamada T (2000) Coptidis rhizoma inhibits growth and proteases of oral bacteria. Oral Dis 6:297–302
Hu JP, Nishishita K, Sakai E, Yoshida H, Kato Y, Tsukuba T, Okamoto K (2008) Berberine inhibits RANKL-induced osteoclast formation and survival through suppressing the NF-kappaB and Akt pathways. Eur J Pharmacol 580:70–79
Huang W, Pi L, Liang W, Xu B, Wang H, Cai R, Huang H (2006) The proteolytic function of the Arabidopsis 26S proteasome is required for specifying leaf adaxial identity. Plant Cell 18:2479–2492
Hudson A (2000) Development of symmetry in plants. Annu Rev Plant Physiol Plant Mol Biol 51:349–370
Huppert JL, Balasubramanian S (2005) Prevalence of quadruplexes in the human genome. Nucleic Acids Res 33:2908–2916
Hwang JM, Kuo HC, Tseng TH, Liu JY, Chu CY (2006) Berberine induces apoptosis through a mitochondria/caspases pathway in human hepatoma cells. Arch Toxicol 80:62–73
Ikezaki M, Kojima M, Sakakibara H, Kojima S, Ueno Y, Machida C, Machida Y (2010) Genetic networks regulated by ASYMMETRIC LEAVES1 (AS1) and AS2 in leaf development in Arabidopsis thaliana: KNOX genes control five morphological events. Plant J 61:70–82
Ishibashi N, Kanamaru K, Ueno Y, Kojima S, Kobayashi T, Machida C, Machida Y (2012) ASYMMETRIC-LEAVES2 and an ortholog of eukaryotic NudC domain proteins repress expression of AUXIN-RESPONSE-FACTOR and class 1 KNOX homeobox genes for development of flat symmetric leaves in Arabidopsis. Biol Open 1:197–207
Islam MM, Kumar GS (2008) RNA targeting by small molecule alkaloids: studies on the binding of palmatine and berberine to polyribonucleotides and comparison to ethidium. J Mol Struct 875:382–391
Islam MM, Chowdhury SR, Kumar GS (2009) Spectroscopic and calorimetric studies on the binding of alkaloids berberine, palmatine and coralyne to double stranded RNA polynucleotides. J Phys Chem B 113:1210–1224
Iwakawa H, Ueno Y, Semiarti E, Onouchi H, Kojima S, Tsukaya H, Hasebe M, Soma T, Ikezaki M, Machida C, Machida Y (2002) The ASYMMETRIC LEAVES2 gene of Arabidopsis thaliana, required for formation of a symmetric flat leaf lamina, encodes a member of a novel family of proteins characterized by cysteine repeats and a leucine zipper. Plant Cell Physiol 43:467–478
Iwakawa H, Iwasaki M, Kojima S, Ueno Y, Soma T, Tanaka H, Semiarti E, Machida Y, Machida C (2007) Expression of the ASYMMETRIC LEAVES2 gene in the adaxial domain of Arabidopsis leaves represses cell proliferation in this domain and is critical for the development of properly expanded leaves. Plant J 51:173–184
Kerstetter RA, Bollman K, Taylor RA, Bomblies K, Poethig RS (2001) KANADI regulates organ polarity in Arabidopsis. Nature 411:706–709
Kidner CA, Martienssen RA (2005) The role of ARGONAUTE1 (AGO1) in meristem formation and identity. Dev Biol 280:504–517
Kim JP, Jung MY, Kim JP, Kim SY (2000) Antiphotooxidative activity of protoberberines derived from Coptis japonica making in the chlorophyll-sensitized photooxidation of oil. J Agric Food Chem 48:1058–1063
Kley AK, Hahn FE (1969) Berberine: complex with DNA. Science 166:755–757
Klímek M, Hnilica L (1959) The influence of deoxyribonucleic acid on ultraviolet and visible light absorption of berberine. Arch Biochem Biophys 81:105–110
Kojima S, Iwasaki M, Takahashi H, Imai T, Matsumura Y, Fleury D, Lijsebettens MV, Machida Y, Machida C (2011) ASYMMETRIC LEAVES2 and elongator, a histone acetyltransferase complex, mediate the establishment of polarity in leaves of Arabidopsis thaliana. Plant Cell Physiol 52:1259–1273
Kuo CL, Chou CC, Yung BY (1995) Berberine complexes with DNA in the berberine-induced apoptosis in human leukemic HL-60 cells. Cancer Lett 93:193–200
Lee HW, Suh JH, Kim HN, Kim AY, Park SY, Shin CS, Choi JY, Kim JB (2008) Berberine promotes osteoblast differentiation by Runx2 activation with p38 MAPK. J Bone Miner Res 23:1227–1237
Lee JS, Latimer LJP, Hampel KJ (1993) Coralyne binds tightly to both T·A·T- and C·G·C+-containing DNA triplexes. Biochemistry 32:5591–5597
Lee JW, Mase N, Yonezawa T, Seo HJ, Jeon WB, Cha BY, Nagai K, Woo JT (2010) Palmatine attenuates osteoclast differentiation and function through inhibition of receptor activator of nuclear factor-κb ligand expression in osteoblast cells. Biol Pharm Bull 33:1733–1739
Lee JW, Iwahashi A, Hasegawa S, Yonezawa T, Jeon WB, Cha BY, Nagai K, Woo JT (2012) Coptisine inhibits RANKL-induced NF-κB phosphorylation in osteoclast precursors and suppresses function through the regulation of RANKL and OPG gene expression in osteoblastic cells. J Nat Med 66:8–16
Li H, Miyahara T, Tezuka Y, Namba T, Suzuki T, Dowaki R, Watanabe M, Nemoto N, Tonami S, Seto H, Kadota S (1999) The effect of kampo formulae on bone resorption in vitro and in vivo. II. Detailed study of berberine. Biol Pharm Bull 22:391–396
Li H, Xu L, Wang H, Yuan Z, Cao X, Yang Z, Zhang D, Xu Y, Huang H (2005) The putative RNA-dependent RNA polymerase RDR6 acts synergistically with ASYMMETRIC LEAVES1 and 2 to repress BREVIPEDICELLUS and microRNA165/166 in Arabidopsis leaf development. Plant Cell 17:2157–2171
Lincoln C, Long J, Yamaguchi J, Serikawa K, Hake S (1994) A knotted1-like homeobox gene in Arabidopsis is expressed in the vegetative meristem and dramatically alters leaf morphology when overexpressed in transgenic plants. Plant Cell 6:1859–1876
Long JA, Moan EI, Medford JI, Barton MK (1996) A member of the KNOTTED class of homeodomain proteins encoded by the STM gene of Arabidopsis. Nature 379:66–69
Ma Y, Ou TM, Hou JQ, Lu YJ, Tan JH, Gu LQ, Huang ZS (2008) 9-N-Substituted berberine derivatives: stabilization of G-quadruplex DNA and down-regulation of oncogene c-myc. Bioorg Med Chem 16:7582–7591
Mantena SK, Sharma SD, Katiyar SK (2006) Berberine, a natural product, induces G1-phase cell cycle arrest and caspase-3-dependent apoptosis in human prostate carcinoma cells. Mol Cancer Ther 5:296–308
Matsumura Y, Iwakawa H, Machida Y, Machida C (2009) Characterization of genes in the ASYMMETRIC LEAVES2/LATERAL ORGAN BOUNDARIES (AS2/LOB) family in Arabidopsis thaliana, and functional and molecular comparisons between AS2 and other family members. Plant J 58:525–537
McConnell JR, Emery J, Eshed Y, Bao N, Bowman J, Barton MK (2001) Role of PHABULOSA and PHAVOLUTA in determining radial patterning in shoots. Nature 411:709–713
Meyerson LR, McMurtrey KD, Davis VE (1978) ISOQUINOLINE ALKALOIDS inhibitory actions on cation-dependent ATP-phosphohydrolases. Neurochem Res 3:239–257
Mullen MA, Olson KJ, Dallaire P, Major F, Assmann SM, Bevilacqua PC (2010) RNA G-Quadruplexes in the model plant species Arabidopsis thaliana: prevalence and possible functional roles. Nucleic Acids Res 38:8149–8163
Nemhauser JL, Feldman LJ, Zambryski PC (2000) Auxin and ETTIN in Arabidopsis gynoecium morphogenesis. Development 127:3877–3888
Ori N, Eshed Y, Chuck G, Bowman JL, Hake S (2000) Mechanisms that control knox gene expression in the Arabidopsis shoot. Development 127:5523–5532
Pekker I, Alvarez JP, Eshed Y (2005) Auxin response factors mediate Arabidopsis organ asymmetry via modulation of KANADI activity. Plant Cell 17:2899–2910
Pérez-Pérez JM, Candela H, Robles P, López-Torrejón G, del Pozo JC, Micol JL (2010) A role for AUXIN RESISTANT3 in the coordination of leaf growth. Plant Cell Physiol 10:1661–1673
Pinon V, Etchells JP, Rossignol P, Collier SA, Arroyo JM, Martienssen RA, Byrne ME (2008) Three PIGGYBACK genes that specifically influence leaf patterning encode ribosomal proteins. Development 135:1315–1324
Rédei GP, Hirono Y (1964) Linkage studies. Arabidopsis Inf Serv 1:9–10
Ren J, Chaires JB (1999) Sequence and structural selectivity of nucleic acid binding ligands. Biochemistry 38:16067–16075
Sakai K, Shitan N, Sato F, Ueda K, Yazaki K (2002) Characterization of berberine transport into Coptis japonica cells and the involvement of ABC protein. J Exp Bot 53:1879–1886
Sawa S, Watanabe K, Goto K, Kanaya E, Morita EH, Okada K (1999) FILAMENTOUS FLOWER, a meristem and organ identity gene of Arabidopsis, encodes a protein with a zinc finger and HMG-related domains. Genes Dev 13:1079–1088
Semiarti E, Ueno Y, Tsukaya H, Iwakawa H, Machida C, Machida Y (2001) The ASYMMETRIC LEAVES2 gene of Arabidopsis thaliana regulates formation of a symmetric lamina, establishment of venation and repression of meristem-related homeobox genes in leaves. Development 128:1771–1783
Serrano-Cartagena J, Robles P, Ponce MR, Micol JL (1999) Genetic analysis of leaf form mutants from the Arabidopsis Information Service collection. Mol Gen Genet 261:725–739
Shitan N, Yazaki K (2007) Accumulation and membrane transport of plant alkaloids. Curr Pharm Biotechnol 8:244–252
Siddiqui-Jain A, Grand CL, Bearss DJ, Hurley LH (2002) Direct evidence for a G-quadruplex in a promoter region and its targeting with a small molecule to repress c-MYC transcription. Proc Natl Acad Sci U S A 99:11593–11598
Steeves TA, Sussex IM (1989) Patterns in plant development, 2nd edn. Cambridge University Press, Cambridge
Szakonyi D, Byrne ME (2011) Ribosomal protein L27a is required for growth and patterning in Arabidopsis thaliana. Plant J 65:269–281
Szakonyi D, Moschopoulos A, Byrne ME (2010) Perspectives on leaf dorsoventral polarity. J Plant Res 123:281–290
Takahashi H, Iwakawa H, Nakao S, Ojio T, Morishita R, Morikawa S, Machida Y, Machida C, Kobayashi T (2008) Knowledge-based fuzzy adaptive resonance theory and its application to the analysis of gene expression in plants. J Biosci Bioeng 106:587–593
Takahashi H, Nakagawa A, Kojima S, Takahashi A, Cha BY, Woo J-T, Nagai K, Machida Y, Machida C (2012) Discovery of novel rules for G-quadruplex-forming sequences in plants by using bioinformatics methods. J Biosci Bioeng
Todd AK, Johnston M, Neidle S (2005) Highly prevalent putative quadruplex sequence motifs in human DNA. Nucleic Acids Res 33:2901–2907
Tsuge T, Tsukaya H, Uchimiya H (1996) Two independent and polarized processes of cell elongation regulate leaf blade expansion in Arabidopsis thaliana (L.) Heynh. Development 122:1589–1600
Tsukaya H (2006) Mechanism of leaf-shape determination. Annu Rev Plant Biol 57:477–496
Tsukaya H, Uchimiya H (1997) Genetic analyses of the formation of the serrated margin of leaf blades in Arabidopsis: combination of a mutational analysis of leaf morphogenesis with the characterization of a specific marker gene expressed in hydathodes and stipules. Mol Gen Genet 256:231–238
Ueno Y, Ishikawa T, Watanabe K, Terakura S, Iwakawa H, Okada K, Machida C, Machida Y (2007) Histone deacetylases and ASYMMETRIC LEAVES2 are involved in the establishment of polarity in leaves of Arabidopsis. Plant Cell 19:445–457
Waites R, Selvadurai HR, Oliver IR, Hudson A (1998) The PHANTASTICA gene encodes a MYB transcription factor involved in growth and dorsoventrality of lateral organs in Antirrhinum. Cell 93:779–789
Watanabe K, Okada K (2003) Two discrete cis elements control the abaxial side-specific expression of the FILAMENTOUS FLOWER gene in Arabidopsis. Plant Cell 15:2592–2602
Wolfe AD, Allison RG, Hahn FE (1972) Labilizing action of intercalating drugs and dyes on bacterial ribosomes. Biochemistry 11:1569–1572
Xu L, Xu Y, Dong A, Sun Y, Pi L, Xu Y, Huang H (2003) Novel as1 and as2 defects in leaf adaxial-abaxial polarity reveal the requirement for ASYMMETRIC LEAVES1 and 2 and ERECTA functions in specifying leaf adaxial identity. Development 130:4097–4107
Xu L, Yang L, Pi L, Liu Q, Ling Q, Wang H, Poethig RS, Huang H (2006) Genetic interaction between the AS1-AS2 and RDR6-SGS3-AGO7 pathways for leaf morphogenesis. Plant Cell Physiol 47:853–863
Yamagishi H (1962) Interaction between nucleic acids and berberine sulfate. J Cell Biol 15:589–592
Yan D, Jin C, Xiao XH, Dong XP (2008) Antimicrobial properties of berberines alkaloids in Coptis chinensis Franch by microcalorimetry. J Biochem Biophys Methods 70:845–849
Yang L, Huang W, Wang H, Cai R, Xu Y, Huang H (2006) Characterizations of a hypomorphic argonaute1 mutant reveal novel AGO1 functions in Arabidopsis lateral organ development. Plant Mol Biol 61:63–78
Yang JY, Iwasaki M, Machida C, Machida Y, Zhou X, Chua NH (2008) βC1, the pathogenicity factor of TYLCCNV, interacts with AS1 to alter leaf development and suppress selective jasmonic acid responses. Genes Dev 22:2564–2577
Yao Y, Ling Q, Wang H, Huang H (2008) Ribosomal proteins promote leaf adaxial identity. Development 135:1325–1334
Zhao Y, Du Z, Li N (2007) Extensive selection for the enrichment of G4 DNA motifs in transcriptional regulatory regions of warm blooded animals. FEBS Lett 581:1951–1956
Acknowledgments
The authors thank Dr. Kazufumi Yazaki (Kyoto University), Dr. Yoko Matsumura (Nagoya University) and Dr. Sumie Keta (Chubu University) for helpful discussions, Dr. Kaname Tsutsumiuchi (Chubu University), and Dr. Yasuyuki Ishida (Chubu University) for helpful suggestions about chemical properties of isoquinoline alkaloids, and Harutaka Ichihara, Takao Yamamoto, Koki Nishimura, and Mari Takahashi (Chubu University) for technical assistance. This work was supported by Grants-in-Aid for Scientific Research on Priority Areas [no. 19060003 to Y.M.], Scientific Research (C) [no. 21570052 and 24570061 to C.M. and no. 22570053 to S.K.] and Young Scientists (B) [no. 21710211 and 24710222 to H.T.] from the Ministry of Education, Culture, Sports, Science and Technology of Japan (MEXT); and the ‘Academic Frontier’ Project for Private Universities (matching fund subsidy from MEXT, 2005–2009). This work is also supported, in part, by a grant from the Japan Foundation for Applied Enzymology.
Open Access
This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.
Author information
Authors and Affiliations
Corresponding author
Additional information
Ayami Nakagawa and Hiro Takahashi contributed equally to this work.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Nakagawa, A., Takahashi, H., Kojima, S. et al. Berberine enhances defects in the establishment of leaf polarity in asymmetric leaves1 and asymmetric leaves2 of Arabidopsis thaliana . Plant Mol Biol 79, 569–581 (2012). https://doi.org/10.1007/s11103-012-9929-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11103-012-9929-7