-
PDF
- Split View
-
Views
-
Cite
Cite
Bin-Bin Ye, Guan-Dong Shang, Yu Pan, Zhou-Geng Xu, Chuan-Miao Zhou, Ying-Bo Mao, Ning Bao, Lijun Sun, Tongda Xu, Jia-Wei Wang, AP2/ERF Transcription Factors Integrate Age and Wound Signals for Root Regeneration, The Plant Cell, Volume 32, Issue 1, January 2020, Pages 226–241, https://doi.org/10.1105/tpc.19.00378
Close - Share Icon Share
Abstract
Age and wounding are two major determinants for regeneration. In plants, the root regeneration is triggered by wound-induced auxin biosynthesis. As plants age, the root regenerative capacity gradually decreases. How wounding leads to the auxin burst and how age and wound signals collaboratively regulate root regenerative capacity are poorly understood. Here, we show that the increased levels of three closely-related miR156-targeted Arabidopsis (Arabidopsis thaliana) SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors, SPL2, SPL10, and SPL11, suppress root regeneration with age by inhibiting wound-induced auxin biosynthesis. Mechanistically, we find that a subset of APETALA2/ETHYLENE RESPONSE FACTOR (AP2/ERF) transcription factors including ABSCISIC ACID REPRESSOR1 and ERF109 is rapidly induced by wounding and serves as a proxy for wound signal to induce auxin biosynthesis. In older plants, SPL2/10/11 directly bind to the promoters of AP2/ERFs and attenuates their induction, thereby dampening auxin accumulation at the wound. Our results thus identify AP2/ERFs as a hub for integration of age and wound signal for root regeneration.
INTRODUCTION
Among the universal changes that occur with age in multicellular organisms is the decline in regenerative capacity. In mammals, for example, the heart quickly loses its capacity to regenerate within a brief period after birth (Porrello et al., 2011). Remyelination, the phenomenon by which new myelin sheaths are generated around axons in the adult central nervous system, declines with increasing age (Ruckh et al., 2012). The factors that are likely involved include accumulated DNA mutations, reduced number of stem cells, decreased cell proliferative capacity, and disturbed metabolic function (Wells and Watt, 2018). In plants, the regenerative capacity also progressively decreases with age. For example, the number of newly formed sprouts decreased significantly in Calluna vulgari plants that were more than 6 years old (Berdowski and Siepel, 1998). Similarly, the shoot regenerative capacity declines as Quercus euboica ages (Kartsonas and Papafotiou, 2007).
Plant hormones play pivotal roles in regeneration (Birnbaum and Sánchez Alvarado, 2008; Sena and Birnbaum, 2010; Perianez-Rodriguez et al., 2014; Pulianmackal et al., 2014; Su and Zhang, 2014; Xu and Huang, 2014; Ikeuchi et al., 2016; Kareem et al., 2016; Sang et al., 2018). It is well known that cytokinin induces shoot regeneration whereas auxin promotes root regeneration (Skoog and Miller, 1957). There are different ways to regenerate roots. In the tissue culture experiments, high auxin/cytokinin in root induction medium triggers the reprogramming of callus, a group of dedifferentiated cells initiated from xylem-pole pericycle cells of root explants and pericycle-like cells of aerial organs, into roots (Atta et al., 2009; Sugimoto et al., 2010; Ikeuchi et al., 2016). By contrast, root tip regeneration follows the developmental stages of embryonic patterning and is guided by spatial information provided by complementary hormone domains (Efroni et al., 2016). The detached Arabidopsis (Arabidopsis thaliana) leaf or hypocotyl can also regenerate roots on hormone-free B5 medium (Liu et al., 2014). It is proposed that the regeneration-competent cells in the pericycle around the wound undergo four sequential steps to initiate adventitious roots (Xu, 2018).
The molecular mechanism by which shoot regenerative capacity is affected by age during tissue culture has been well documented. Our previous work reveals an important role of microRNA156 (miR156), the master regulator of juvenility, in shoot regeneration (Poethig, 2009; Yu et al., 2015c). miR156 targets a group of transcription factors named SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL; Rhoades et al., 2002; Schwab et al., 2005; Wu and Poethig, 2006; Gandikota et al., 2007). Juvenile plants exhibit a high cytokinin response and regenerative capacity. As a plant ages, miR156 levels decline, alleviating the repression of its SPL targets (Wang et al., 2009; Yamaguchi et al., 2009; Xu et al., 2018). SPL directly inhibits the transcriptional activity of B-type ARABIDOPSIS RESPONSE REGULATORs and thereby impairs shoot regenerative capacity (Zhang et al., 2015). miR156 is also involved in de novo root regeneration in Arabidopsis (Xu et al., 2016; Ye et al., 2019). The hypocotyls from the plants with low miR156 activity produce fewer adventitious roots than those from wild type. Notably, this phenomenon seems to be evolutionarily conserved. The same positive relationship between miR156 level and rooting capacity has been found in Malus xiaojinensis and Medicago sativa (Aung et al., 2017; Xu et al., 2017). However, the molecular mechanism by which miR156/SPL regulates root regenerative capacity is poorly understood.
Here, we find that miR156-targeted SPLs have divergent roles in plant regeneration. Whereas the SPL9 subfamily suppressed shoot regeneration, the SPL10 subfamily specifically suppressed root regeneration by inhibiting auxin burst in response to wounding. By the combination of RNA-sequencing (RNA-seq) and chromatin immunoprecipitation followed by sequencing (ChIP-seq), we further identified a subset of APETALA2/ETHYLENE RESPONSE FACTOR (AP2/ERF) transcription factors as the proxy for wound signal to induce auxin biosynthesis, and these transcription factors act as a hub for integration of age and wound signal.
RESULTS
Divergent Roles of miR156-Targeted SPLs in Plant Regeneration
In Arabidopsis, miR156 targets 11 SPL genes (Xing et al., 2013). Whether all the SPL members contribute to root regeneration is still unknown. To address this question, we examined root regeneration capacity by culturing leaf explants on B5 medium without added exogenous hormones under long-day conditions (Liu et al., 2014). To avoid the impact of leaf age on regenerative capacity, the first/second (early/juvenile), third/fourth (mid), and sixth (late/adult) Arabidopsis leaves of the same developmental stage (3 mm in length) were used. Consistent with published results (Xu et al., 2016), the rooting capacity was progressively decreased with plant age in wild type (Supplemental Figures 1A and 1B). The first to second leaf explants exhibited higher miR156 level than the sixth leaf explants (Supplemental Figure 1C).
Blocking the function of miR156 by overexpression of MIM156 (mimicry of miR156) greatly reduced adventitious root formation (Supplemental Figures 1D and 1E). An opposite trend was observed when miR156 was overexpressed by constitutive 35S promoter (Supplemental Figures 1D and 1E). Notably, the same relationship between miR156 level and rooting capacity was found in poplar tree (Populus; Supplemental Figures 1F and 1G), indicating that the role of miR156 in rooting is evolutionarily conserved (Wang et al., 2011).
We next performed the regeneration experiments using the first to second, third to fourth, and sixth leaf explant of spl single and high-order mutants. SPL9 and SPL15 have been reported to play an important role in shoot regeneration (Zhang et al., 2015). The sixth leaf explants of spl9 spl15 exhibited the same rooting rate as wild type (Figures 1A and 1B). SPL2, SPL10, and SPL11 constitute another miR156-targeted SPL subfamily. Because SPL10 and SPL11 are located in the same genetic locus, we generated spl10 spl11 double mutants by the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 technology (Supplemental Figures 2A to 2C; Ran et al., 2017; Yin et al., 2017). Compared to wild type, spl10 spl11 showed increased rooting capacity (Figures 1A to 1C; Supplemental Figures 2D to 2I). This phenotype was further enhanced in spl2 spl10 spl11 triple mutants (Figures 1A to 1C; Supplemental Figures 2D to 2I). Intriguingly, both spl10 spl11 and spl2 spl10 spl11 triple mutant did not affect shoot regeneration capacity in tissue culture (Supplemental Figure 2J and 2K).
SPL10 Subfamily Regulates Rooting.
(A) Root regeneration assays using the sixth leaf explants from wild type and spl mutants. Pictures were taken at 20 DAC on B5 media. Scale bar = 1 cm.
(B) Quantitative analysis of root regenerative rate of the sixth leaf explants in spl single and high-order mutants at 20 DAC. The rooting rate was represented by the percentage of leaf explants with regenerated roots. Data are from five independent experiments. Bars show mean. One-way ANOVA was performed followed by a Tukey’s multiple comparisons test. **P < 0.01; ns, not significant.
(C) The root regenerative capacity of the sixth leaf explants from wild type, spl10 spl11, and spl2 spl10 spl11 at 12 DAC. The regenerative capacity was represented by the percentage of leaf explants with different number of regenerated roots. Data are from three independent experiments. Bars show mean ± sd.
(D) The rooting rate of leaf explants cultured on B5 medium without (mock) or with 10 μΜ of DEX treatment at 20 DAC. Data are from three independent experiments. Bars show mean. Two-sided Student’s t test was performed. **P < 0.01; ns, not significant.
(E) The rooting rate of leaf explants cultured on B5 medium without (mock) or with 10 μΜ of DEX treatment at 20 DAC. Data are from three independent experiments. Bars show mean. Two-sided Student’s t test was performed.**P < 0.01; ns, not significant.
(F) The rooting rate of leaf explants overexpressing rSPL9 or rSPL10 at 20 DAC. Data are from three independent experiments. Bars show mean. One-way ANOVA was performed followed by a Tukey’s multiple comparisons test. **P < 0.01; ns, not significant.
(G) and (I) Expression pattern of SPL9 and SPL10 promoter GUS reporters. The first to second leaves (G) and the sixth leaves (I) were stained for 5 h for pSPL9:GUS and 1 h 40 min for pSPL10:GUS. Scale bar = 1 mm. Three independent experiments were performed (n = 15 for each experiment). One representative picture is shown.
(H) and (J) Transverse section of pSPL9:GUS and pSPL10:GUS reporters. The first to second leaves (H) and the sixth leaves (J) were stained for the same time (24 h). xy, xylem; xp, xylem parenchyma cell; pc, procambium; ph, phloem. Scale bar = 10 μm.
We then generated SPL9- and SPL10-inducible lines (pSPL9:rSPL9-GR and pSPL10:rSPL10-GR), in which miR156-resistant forms of SPLs (rSPL) were fused to the hormone binding domain of rat glucocorticoid receptor (GR) and expressed from their own promoters. Treatment with the steroid hormone ligand dexamethasone (DEX), which leads to a translocation of the rSPL9/10-GR fusion protein from the cytoplasm to the nucleus, resulted in the same phenotype as the transgenic plants expressing rSPL9/10 (pSPL9:rSPL9/pSPL10:rSPL10; Supplemental Figure 3A). In the absence of DEX, pSPL9:rSPL9-GR and pSPL10:rSPL10-GR rooted as wild type (Figure 1D; Supplemental Figure 4A). By contrast, after application of DEX, the rooting capacity of pSPL10:rSPL10-GR was nearly abolished, whereas that of pSPL9:rSPL9-GR remained unchanged (Figure 1D; Supplemental Figure 4A).
To probe the reason why SPL10/11, but not SPL9/15, represses root regeneration, we performed a promoter swapping experiment by generating pSPL10:rSPL9, pSPL9:rSPL10, pSPL10:rSPL9-GR, and pSPL9:rSPL10-GR. Both pSPL10:rSPL9 and pSPL9:rSPL10 accelerated the juvenile-to-adult phase transition (Supplemental Figure 3B). Intriguingly, pSPL10:rSPL9 exhibited the same low rooting rate as pSPL10:rSPL10 plants (Figure 1F; Supplemental Figure 4C). Consistently, pSPL9:rSPL10 did not show altered rooting rate in comparison to wild type (Figure 1F; Supplemental Figure 4C). Similar results were obtained using pSPL10:rSPL9-GR and pSPL9:rSPL10-GR as explants (Figure 1E; Supplemental Figure 4B).
Analyses of SPL9 and SPL10 promoter β-Glucuronidase (GUS) reporter lines (two independent lines) revealed that pSPL10:GUS showed stronger staining in the procambium cells and some xylem parenchyma cells than pSPL9:GUS (Figures 1G to 1J; Supplemental Figures 3C to 3F), suggesting that the difference in promoter strength may underlie the functional divergence in root regeneration between SPL9/15 and SPL10/11. Consistent with the notion that SPL10 and SPL11 play redundant roles in rooting, SPL11 exhibited the same expression pattern as SPL10 (Supplemental Figure 3G). Taking these observations together, we conclude that the SPL9 subfamily is important for shoot regeneration, whereas SPL2, SPL10, and SPL11 play a dominant role in root regeneration.
SPL2/10/11 Regulate Root Regenerative Capacity Partially through Auxin
The formation of adventitious roots is triggered by the accumulation of auxin at the wound after leaf detachment (Xu, 2018). The impairment of root regeneration in adult wild type and p35S:MIM156 and pSPL10:rSPL10 plants, therefore, might be caused by defects in auxin biosynthesis or signal transduction. To distinguish between these two possibilities, we first examined the auxin responsiveness using DR5 as a reporter (Sabatini et al., 1999). Consistent with published results, the DR5 expression was low before wounding (Liu et al., 2014; Chen et al., 2016b). After 1 d of leaf detachment, strong DR5 signals were observed in mesophyll cells and the vasculature near the wound in the first to second leaf (Figure 2A). By contrast, DR5 expression was not greatly increased in both wound sites of the third to fourth and fifth to sixth leaves (Figure 2A; Supplemental Figure 5A). Comparison of the first to second leaves among wild type, p35S:MIR156, and p35S:MIM156 further confirmed a positive relationship between miR156 level and the enhancement of DR5 expression during root regeneration (Figures 2B and 2C; Supplemental Figures 5B and 5C).
SPL10 Subfamily Regulates Rooting Capacity through Auxin.
(A) Auxin responsiveness in wild type. Auxin responsiveness is inferred by pDR5:GFP reporter. The leaves from the plants of different ages (first to second, third to fourth, fifth to sixth) were used as explants. Four independent experiments were performed (n > 15 for each experiment). One representative result is shown. The leaf shape is marked with dashed lines. Scale bar = 500 μm.
(B) and (C) Auxin responsiveness in wild-type, p35S:MIR156, and p35S:MIM156 plants. The first to second (B) or sixth leaves (C) were used as explants. Three independent experiments were performed (n > 15 for each experiment). One representative result is shown. The leaf shape is marked with dash lines. Scale bars = 500 μm.
(D) Wound-induced IAA accumulation in wild-type leaf explants was attenuated in old plants. The ratio was calculated by IAA concentration at 2 DAC/IAA concentration at 0 DAC. The accumulation of auxin at the wound after leaf detachment at 2 DAC has been reported (Xu, 2018). The IAA ratio in the first to second wild-type leaf explants was set to 1. Four independent experiments were performed. Lines show mean. Two-sided Student’s t test was performed. **P < 0.01.
(E) and (F) Wound-induced IAA accumulation was modulated by SPL10 subfamily. The first to second (E) or sixth leaves (F) were used as explants. The IAA ratio in the first to second (E) or sixth (F) wild-type leaf explants was set to 1.0. Four (E) and three (F) independent experiments were performed, respectively. Lines show mean. Two-sided Student’s t test was performed. *P < 0.05; **P < 0.01.
(G) and (H) The expression level of ASA1 in the first to second (G) or six (H) leaf explants from wild type, spl2 spl10 spl1, and pSPL10:rSPL10 at 0, 1, and 2 DAC. The expression levels in the leaves of pSPL10:rSPL10 at 0 DAC (G) and the leaves of wild type at 0 DAC (H) were set to 1. Four (G) and two (H) independent experiments were performed, respectively. Lines show mean. One-way ANOVA followed by a Tukey’s multiple comparisons test or two-sided Student’s t test was performed. *P < 0.05; **P < 0.01; ns, not significant.
(I) and (J) Quantitative analyses of rooting rate of the first to second (I) or sixth (J) leaf explants cultured on B5 media supplemented with IAA of different concentrations. Data of spl2 spl10 spl11 and pSPL10:rSPL10 were from three independent experiments. Data of wild type were from six (I) or four (J) independent experiments.
To validate our results, we measured free indole-3-acetic acid (IAA) levels in wild type, spl2 spl10 spl1, and pSPL10:rSPL10 explants. An auxin burst was evident in the first to second wild-type leaves at 2 d after culturing (DAC; Figure 2D). This effect was greatly suppressed in the fifth to sixth leaves (Figure 2D). In agreement with this, the accumulation of IAA was attenuated in the first to second leaves of pSPL10:rSPL10 (Figure 2E) but boosted in the sixth leaf of spl2 spl10 spl11 mutant (Figure 2F; Supplemental Figure 5D).
We next examined the expression of genes involved in auxin biosynthesis. WEAK ETHYLENE INSENSITIVE2/ANTHRANILATE SYNTHASE alpha1 and WEI7/ANTHRANILATE SYNTHASE beta1 (ASB1) genes that encode alpha- and beta-subunits of a rate-limiting enzyme in Trp biosynthesis pathway have a major role in Trp-dependent auxin biosynthesis (Stepanova et al., 2005, 2008; Tao et al., 2008). TRYPTOPHAN AMINOTRANSFERASE OF ARABIDOPSIS1 (TAA1) catalyzes the conversion of Trp to indole-3-pyruvate (Stepanova et al., 2008). Compared to wild type, the expression of the ASA1, ASB1, and TAA1 was low in pSPL10:rSPL10 (Figure 2G; Supplemental Figure 5E). Upon leaf detachment, both ASA1 and ASB1 were induced in the first to second wild-type leaves within 2 d, whereas TAA1 level remained constant. The upregulation of ASA1 and ASB1 was significantly suppressed in pSPL10:rSPL10 (Figure 2G; Supplemental Figure 5E). When the sixth leaves were used as explants, the induction of ASA1 was compromised in wild type but still evident in spl2 spl10 spl11 mutant (Figure 2H). Thus, these results support the notion that the reduced rooting capacity of pSPL10:rSPL10 is likely caused by a reduced amount of auxin in leaf explants.
We then asked whether exogenous auxin is sufficient to rescue the rooting defects of pSPL10:rSPL10 plants. Compared to mock, the addition of IAA led to increased rooting rate of the first to second leaves of pSPL10:rSPL10 in a dose-dependent manner (Figure 2I; Supplemental Figure 5F). This effect was also observed in the sixth wild-type leaves but to a lesser extent in spl2 spl10 spl11 mutant background (Figure 2J; Supplemental Figure 5G). It should be noted that IAA could not fully rescue the rooting defects of pSPL10:rSPL10, implying that SPL10 does not suppress root regeneration exclusively through modulating auxin biosynthesis.
Identification of Direct Targets of SPL10 during Root Regeneration
The mechanism by which auxin sequentially induces root founder cell specification and cell fate transition has been well studied (Xu, 2018). However, the primary wounding response within hours during root regeneration is largely unknown. We harvested wild-type and pSPL10:rSPL10 leaf explants at 0 and 4 h after leaf detachment. Total RNAs were extracted and subjected to RNA-seq analyses. In wild type, a total of 3,221 genes were differentially expressed (twofold, P < 0.01). Among them, 63% of the genes (2,057/3,221) were upregulated and 37% of the genes (1,164/3,221) were downregulated (Figure 3A; Supplemental Data Sets 1 and 2). Gene ontology (GO) enrichment analysis revealed that the upregulated genes were enriched with the genes involved in jasmonic acid (JA) signaling, chitin response, and substance transportation (Supplemental Figure 6A; Supplemental Data Set 3). Further comparison with the RNA-seq data of pSPL10:rSPL10 leaf explants identified that 590 wound-induced genes (28% [590/2,057]) were repressed by SPL10 (Figure 3A; Supplemental Data Set 4, twofold, P < 0.01).
Identification of Downstream Targets of SPL10.
(A) Volcano plot. The explants were cultured on B5 media for 4 h and subjected to RNA-seq analyses. A total of 2,057 genes were induced, whereas 1,164 genes were repressed. FDR < 0.01; Log2Fold Change > 1 or < −1. Among the wound-induced genes, 590 genes were repressed by SPL10 (1 < Log2FoldChange [pSPL10:rSPL10] < Log2FoldChange [wild type]).
(B) The enrichment of GTAC motif in SPL10 ChIP-seq. P value is given at upright corner. The e represents times 10 raised to the power (×10n).
(C) Venn plot. ChIP-Seq analysis revealed 8,569 genes targeted by SPL10. Among them, 233 genes were wound-inducible.
(D) GO analysis of 233 genes that are wound-induced and repressed by SPL10. The top 20 categories are shown.
(E) Fold induction of JA-responsive genes in wild type (Col-0) and pSPL10:rSPL10 during root regeneration. The induction fold was calculated by comparison of the expression level between 0 and 4 h.
(F) Fold induction of SA-responsive genes in wild type and pSPL10:rSPL10 during root regeneration. The induction fold was calculated by comparison of the expression level between 0 and 4 h.
The RNA-seq analyses revealed that SPL10 has a global impact on wounding response by affecting expression of >3,000 genes involved in multiple plant hormone and physiological processes. However, some of these effects may be indirect. To resolve this question, we mapped the SPL10 targets in vivo by ChIP-seq using the 3xFLAG epitope-tagged version of rSPL10 expressed under its own promoter (pSPL10:3xFLAG-rSPL10). pSPL10:3xFLAG-rSPL10 showed the similar phenotype as pSPL10:rSPL10 (Supplemental Figures 6B and 6C). The explants after 4 h of leaf detachment were harvested and subjected to ChIP-seq analysis. In total, we identified ~10,506 SPL10 binding peaks (corresponding to 8,569 genes) in the genome (Supplemental Data Set 5), in which 233 genes were wound-induced and SPL10-repressed (Figure 3C; Supplemental Data Set 6). GO analysis further revealed that the genes in response to JA, salicylic acid (SA), and wounding were overrepresented (Figures 3D to 3F; Supplemental Data Set 7). Among them, 28 genes encoded transcription factors, belonging to 12 different transcription factor families (Table 1).
List of 28 Genes Encoding Transcription Factors, Induced By Wounding and Suppressed by SPL10.
| Gene Family . | Number . | Arabidopsis Genome Initiative Number and Gene Name . | |||
|---|---|---|---|---|---|
| B3 | 2 | AT3G11580 (NGAL2) | AT5G06250 (NGAL3) | — | — |
| bHLH | 3 | AT1G10585 | AT2G43010 (PIF4) | AT5G61270 (PIF7) | — |
| C2H2 | 1 | AT2G37430 (ZAT11) | — | — | — |
| DBB | 1 | AT1G78600 (BBX22) | — | — | — |
| AP2/ERF | 8 | AT1G43160 (RAP2.6) | AT3G50260 (CEJ1) | AT4G06746 (RAP2.9) | AT4G34410 (ERF109) |
| AT5G07310 (ERF115) | AT5G13330 (RAP2.6L) | AT5G53290 (CRF3) | AT5G64750 (ABR1) | ||
| LBD | 1 | AT2G28500 (LBD11) | — | — | — |
| MYB | 2 | AT2G47190 (MYB2) | AT4G21440 (MYB102) | — | — |
| MYB_related | 1 | AT3G23250 (MYB15) | — | — | — |
| NAC | 2 | AT1G52890 (NAC019) | AT3G04070 (NAC047) | — | — |
| SRS | 1 | AT5G12330 (LRP1) | — | — | — |
| Trihelix | 2 | AT2G38250 | AT5G01380 | — | — |
| WRKY | 4 | AT1G80840 (WRKY40) | AT5G07100 (WRKY26) | AT5G46350 (WRKY8) | AT5G49520 (WRKY48) |
| Gene Family . | Number . | Arabidopsis Genome Initiative Number and Gene Name . | |||
|---|---|---|---|---|---|
| B3 | 2 | AT3G11580 (NGAL2) | AT5G06250 (NGAL3) | — | — |
| bHLH | 3 | AT1G10585 | AT2G43010 (PIF4) | AT5G61270 (PIF7) | — |
| C2H2 | 1 | AT2G37430 (ZAT11) | — | — | — |
| DBB | 1 | AT1G78600 (BBX22) | — | — | — |
| AP2/ERF | 8 | AT1G43160 (RAP2.6) | AT3G50260 (CEJ1) | AT4G06746 (RAP2.9) | AT4G34410 (ERF109) |
| AT5G07310 (ERF115) | AT5G13330 (RAP2.6L) | AT5G53290 (CRF3) | AT5G64750 (ABR1) | ||
| LBD | 1 | AT2G28500 (LBD11) | — | — | — |
| MYB | 2 | AT2G47190 (MYB2) | AT4G21440 (MYB102) | — | — |
| MYB_related | 1 | AT3G23250 (MYB15) | — | — | — |
| NAC | 2 | AT1G52890 (NAC019) | AT3G04070 (NAC047) | — | — |
| SRS | 1 | AT5G12330 (LRP1) | — | — | — |
| Trihelix | 2 | AT2G38250 | AT5G01380 | — | — |
| WRKY | 4 | AT1G80840 (WRKY40) | AT5G07100 (WRKY26) | AT5G46350 (WRKY8) | AT5G49520 (WRKY48) |
The gene family, the number of the genes within the same family, and accession numbers are shown. Dashes indicate no data.
| Gene Family . | Number . | Arabidopsis Genome Initiative Number and Gene Name . | |||
|---|---|---|---|---|---|
| B3 | 2 | AT3G11580 (NGAL2) | AT5G06250 (NGAL3) | — | — |
| bHLH | 3 | AT1G10585 | AT2G43010 (PIF4) | AT5G61270 (PIF7) | — |
| C2H2 | 1 | AT2G37430 (ZAT11) | — | — | — |
| DBB | 1 | AT1G78600 (BBX22) | — | — | — |
| AP2/ERF | 8 | AT1G43160 (RAP2.6) | AT3G50260 (CEJ1) | AT4G06746 (RAP2.9) | AT4G34410 (ERF109) |
| AT5G07310 (ERF115) | AT5G13330 (RAP2.6L) | AT5G53290 (CRF3) | AT5G64750 (ABR1) | ||
| LBD | 1 | AT2G28500 (LBD11) | — | — | — |
| MYB | 2 | AT2G47190 (MYB2) | AT4G21440 (MYB102) | — | — |
| MYB_related | 1 | AT3G23250 (MYB15) | — | — | — |
| NAC | 2 | AT1G52890 (NAC019) | AT3G04070 (NAC047) | — | — |
| SRS | 1 | AT5G12330 (LRP1) | — | — | — |
| Trihelix | 2 | AT2G38250 | AT5G01380 | — | — |
| WRKY | 4 | AT1G80840 (WRKY40) | AT5G07100 (WRKY26) | AT5G46350 (WRKY8) | AT5G49520 (WRKY48) |
| Gene Family . | Number . | Arabidopsis Genome Initiative Number and Gene Name . | |||
|---|---|---|---|---|---|
| B3 | 2 | AT3G11580 (NGAL2) | AT5G06250 (NGAL3) | — | — |
| bHLH | 3 | AT1G10585 | AT2G43010 (PIF4) | AT5G61270 (PIF7) | — |
| C2H2 | 1 | AT2G37430 (ZAT11) | — | — | — |
| DBB | 1 | AT1G78600 (BBX22) | — | — | — |
| AP2/ERF | 8 | AT1G43160 (RAP2.6) | AT3G50260 (CEJ1) | AT4G06746 (RAP2.9) | AT4G34410 (ERF109) |
| AT5G07310 (ERF115) | AT5G13330 (RAP2.6L) | AT5G53290 (CRF3) | AT5G64750 (ABR1) | ||
| LBD | 1 | AT2G28500 (LBD11) | — | — | — |
| MYB | 2 | AT2G47190 (MYB2) | AT4G21440 (MYB102) | — | — |
| MYB_related | 1 | AT3G23250 (MYB15) | — | — | — |
| NAC | 2 | AT1G52890 (NAC019) | AT3G04070 (NAC047) | — | — |
| SRS | 1 | AT5G12330 (LRP1) | — | — | — |
| Trihelix | 2 | AT2G38250 | AT5G01380 | — | — |
| WRKY | 4 | AT1G80840 (WRKY40) | AT5G07100 (WRKY26) | AT5G46350 (WRKY8) | AT5G49520 (WRKY48) |
The gene family, the number of the genes within the same family, and accession numbers are shown. Dashes indicate no data.
De novo motif analyses revealed eight enriched putative transcription factor binding elements under SPL10 binding peaks. As expected, many peaks contain perfect GTAC motifs, the core binding motif for SPL transcription factors (Figure 3B; Klein et al., 1996). Interestingly, the binding motifs for basic helix-loop-helix, MYB, and NAC transcription factors were also enriched (Supplemental Figure 6E), implying that the transcription factors belonging to these three gene families may play a collaborative role with SPL10 in root regeneration.
Convergence of Age and Wound Signals on a Subset of AP2/ERFs
One of the primary events of the response to wounding that occur at the wounded site is the burst of JA (Savatin et al., 2014; Chen et al., 2016a; Hilleary and Gilroy, 2018). Our GO analyses revealed that SPL10 represses JA signaling and biosynthesis during root regeneration. Indeed, a recent article has demonstrated that miR156-targeted SPLs can stabilize JA ZIM-domain proteins, the repressor of JA signaling pathway, thereby attenuating JA responses (Mao et al., 2017).
ChIP-seq also identified eight AP2/ERFs as direct targets of SPL10 (Figure 4A; Supplemental Data Sets 1 and 5; Licausi et al., 2013). The binding of SPL10 on ABSCISIC ACID REPRESSOR1 (ABR1) promoter was verified by ChIP-PCR analysis and competitive electrophoretic mobility shift assay (EMSA; Figure 4B; Supplemental Figures 7A to 7C). Interestingly, ERF109 expression has been shown to be induced by JA (Cai et al., 2014). RT-qPCR analyses found that another seven AP2/ERFs genes also responded to JA treatment (Figures 4C to 4F; Supplemental Figures 7D to 7G). Time-course analyses further revealed that these AP2/ERFs were rapidly induced after leaf detachment and reduced to the basal levels after 1 to 2 h (Figure 4G; Supplemental Figures 7H and 8A). Consistent with RNA-seq data, the induction of AP2/ERFs was largely compromised in pSPL10:rSPL10 leaf explants (Figure 4G; Supplemental Figures 7H and 8A). Taken together, these results suggest that age and wound signals converge on AP2/ERFs during root regeneration.
AP2/ERFs Are Direct Target of SPL10.
(A) ChIP-seq data showing that SPL10 binds to the promoter region of ABR1. Schematic of ABR1 gene shows the regions for ChIP-PCR analysis and the locations of SPL10 binding motif (GTAC; inverted triangles). The data was shown by Gviz (Hahne and Ivanek, 2016).
(B) ChIP-PCR analysis. Data are from three independent experiments. Bars show mean ± sd.
(C) to (F) Expression of AP2/ERF transcription factors in response to JA treatment. For ABR1 and ERF109, explants were treated with JA for 30 min. For ERF115 and RAP2.6L, the explants were treated for 60 min. Data are from three independent experiments. Lines show mean.
(G) Time course analyses of AP2/ERFs during root regeneration. Two independent experiments were performed. One representative data are shown. Bars show mean ± sd. X-axis stands for the culturing time after leaf detachment. Result for another biological replicate is shown in Supplemental Figure 8.
(H) GUS staining of pABR1:GUS. The leaves were detached and stained at 0, 30, and 60 min. Scale bar = 200 μm. Three independent experiments were performed (n = 15 for each experiment). One representative picture is shown.
The Role of ABR1 in Root Regeneration
Previous reports have shown that ERF109 promotes lateral root initiation by activating ASA1 gene expression (Cai et al., 2014; Kong et al., 2018), whereas ERF115 contributes to root stem cell replenishment upon damage (Heyman et al., 2013). Moreover, RELATED TO AP2 6L (RAP2.6L) is involved in tissue reunion and required for the division of pith cells in the reunion process (Asahina et al., 2011). We therefore postulate that these AP2/ERF genes may play important roles in wound-induced auxin biosynthesis and root regeneration. Because the enrichment of SPL10 on the ABR1 promoter was highest among these AP2/ERFs, we characterized its expression and function in detail.
To reveal the expression pattern of ABR1 during root regeneration, we generated a reporter line, in which GUS was expressed from the regulatory sequence of ABR1 (pABR1:GUS). Under normal growth condition, GUS signals were weakly detected in leaf vascular tissues (Figure 4H). By contrast, upon wounding, the reporter expression was significantly increased at the wound after 30 min, and throughout the leaf after 1 h (Figure 4H).
To probe the role of AP2/ERFs in root regeneration, we generated transgenic plants that overexpressed ABR1 using the constitutive 35S promoter. The p35S:ABR1 plants did not exhibit growth defects during the vegetative phase (Supplemental Figure 9A). Root regeneration assays revealed that p35S:ABR1 conferred a higher root regenerative capacity than wild type (Figures 5A and 5B; Supplemental Figure 9C). Notably, p35S:ABR1 was sufficient to rescue the low root regenerative capacity in pSPL10:rSPL10 plants (Figures 5C and 5D; Supplemental Figure 9D), supporting the notion that ABR1 acts downstream of SPL10/11.
The Role of ABR1 in Root Regeneration.
(A) The rooting rate of the sixth leaf explants from wild type (Col-0) and p35S:ABR1 at 20 DAC. The rooting rate was calculated by the number of rooted explants/total number of explants. Three independent p35S:ABR1 were used. Data are from three independent experiments. One-way ANOVA was performed followed by a Tukey’s multiple comparisons test. **P < 0.001.
(B) The regenerative capacity of the sixth leaf explants from wild type (Col-0) and p35S:ABR1 at 12 DAC. The regenerative capacity was represented by the percentage of leaf explants with different number of regenerated roots. Three independent p35S:ABR1 were used. Data are from three independent experiments. Bars show mean ± sd.
(C) The rooting rate of the first to second leaf explants from pSPL10:rSPL10 and p35S:ABR1 pSPL10:rSPL10 at 20 DAC. Data are from four independent experiments. Bars show mean. One-way ANOVA was performed followed by a Tukey’s multiple comparisons test. **P < 0.01.
(D) The regenerative capacity of the first to second leaf explants from pSPL10:rSPL10 and p35S:ABR1 pSPL10:rSPL10 at 12 DAC. Data are from four independent experiments. Bars show mean ± sd.
(E) The rooting rate of the first to second leaf explants from erf mutants at 20 DAC. Data are from five independent experiments. Bars show mean. One-way ANOVA was performed followed by a Tukey’s multiple comparisons test. *P < 0.05; **P < 0.01; ns, not significant.
(F) The regenerative capacity of the first to second leaf explants from erf mutants at 12 DAC. Data are from five independent experiments. Bars show mean ± sd.
(G) The rooting rate of the first to second leaf explants from wild type (Col-0) and p35S:ABR1-SRDX at 20 DAC. Data are from three independent experiments. Bars show mean. One-way ANOVA was performed followed by a Tukey’s multiple comparisons test. **P < 0.01.
(H) The regenerative capacity of the first to second leaf explants from wild type (Col-0) and p35S:ABR1-SRDX at 12 DAC. Data are from three independent experiments. Bars show mean ± sd.
(I) Root regeneration assays using wild type and erf mutants. Pictures were taken at 20 DAC. Scale bar = 1 cm.
(J) ChIP-PCR analysis showing enrichment of ABR1-3xFLAG on the promoter of ASA1. Schematic of ASA1 genomic region shows the region for ChIP-PCR analysis (black overbar) and the location of GCC-box (inverted green triangles). The enrichment level in wild type (Col-0) was set to 1. Data are from three independent experiments. Bars show mean. Two-sided Student’s t test was performed. *P < 0.05.
(K) Wound-induced IAA accumulation in the first to second leaf explants from wild type and ERF mutant. The IAA ratio in the first to second leaf explants of p35S:ABR1-SRDX #1 was set to 1. Three independent experiments were performed. Lines show mean. One-way ANOVA was performed followed by a Tukey’s multiple comparisons test. **P < 0.01.
AP2/ERFs constitute a large family that shares highly conserved DNA binding domains, resulting in much redundancy (Asahina et al., 2011; Licausi et al., 2013; Gu et al., 2017; Phukan et al., 2017). To overcome this, we overexpressed a dominant-negative version of ABR1 by fusion ABR1 with the repressor motif SRDX (Hiratsu et al., 2003). Compared to wild type, the rooting capacity of p35S:ABR1-SRDX was completely abolished (Figures 5G and 5H; Supplemental Figure 9E). To validate our results, we obtained a T-DNA insertion mutant of abr1 and erf109, which did not exhibit growth defects during vegetative phase compared to wild type (Supplemental Figure 9B). Rooting assays revealed that erf single mutants did not exhibit strong defects in rooting, whereas abr1 erf109 double mutant showed a lower rooting rate than wild type (Figures 5E, 5F, and 5I). The defect of abr1 erf109 mutant was less severe than p35S:ABR1-SRDX, suggesting that the other six SPL10/11-regulated AP2/ERF genes may play redundant roles in root regeneration.
It has been shown that JA activates root stem cells (Chen et al., 2011; Marhava et al., 2019; Zhou et al., 2019). In addition, JA and ERF109 evoke auxin biosynthesis through ASA1 (Sun et al., 2009; Cai et al., 2014). To investigate whether ABR1 also promotes root regeneration by directly binding to the ASA1 promoter, we generated transgenic plants in which 3xFLAG-tagged ABR1 was expressed under a constitutive RIBO promoter (pRIBO:ABR1-3xFLAG, Supplemental Figure 9F). ChIP-PCR experiments found that ABR1-3xFLAG could bind to the GCC-box in the promoter of ASA1 (Figure 5J). Consistent with this, the induction of ASA1 was largely compromised in erf109 abr1 mutants (Supplemental Figure 9G), accompanied by a lower level of IAA (Figure 5K; Supplemental Figure 9H). Taken together, these results indicate that both ERF109 and ABR1 promote auxin production and rooting by activating ASA1.
DISCUSSION
Nature of the Wound Signal during Plant Regeneration
The very early responses to wounding include the changes in trans-membrane potential, increase of apoplastic glutamate, influx of Ca2+, and generation of H2O2 (Lup et al., 2016; Choi et al., 2017; Toyota et al., 2018). However, how these signals regulate gene expression in the nucleus is largely unknown. The identification of wound-induced AP2/ERFs suggests the early wound signal rapidly reprograms the transcriptome at the transcriptional level (Marhava et al., 2019; Zhang et al., 2019; Zhou et al., 2019). Intriguingly, studies from the Sugimoto lab have shown that another group of AP2/ERF-type transcriptional factor genes, named WOUND-INDUCED DEDIFFERENTIATION1 (WIND1) and its homologs WIND2, WIND3 and WIND4, is induced upon wounding (Iwase et al., 2011a; 2011b). Importantly, transient overexpression of WIND genes is able to promote callus formation at the wound. These results collectively imply that AP2/ERFs act as a proxy for wounding signal to evoke downstream gene expression and different groups of AP2/ERFs may trigger divergent regeneration processes (Heyman et al., 2018).
How the wound signal induces AP2/ERF expression remains to be elucidated. Expression analyses reveal that AP2/ERFs such as ABR1, ERF115, and ERF109 are induced not only by wounding but also by plant hormones such as abscisic acid, ethylene, JA, and SA, suggesting a common property shared by the ERF promoters. Interestingly, a recent report reveals the mechanism by which ERF1 is rapidly induced by JA (An et al., 2017). The transcriptional coactivator complex Mediator directly links JA receptor CORONATINE INSENSITIVE1 to the ERF1 promoter, thereby facilitating CORONATINE INSENSITIVE1-dependent degradation of JA ZIM-domain transcriptional repressors. In this scenario, we speculate that the acute induction of the wound-induced AP2/ERFs may be triggered by the rapid degradation of unknown transcriptional repressors that bind to the AP2/ERF promoters before wounding.
SPL10 Suppresses Root Regeneration through Multiple Pathways
The identification of AP2/ERFs as the direct of SPL10 suggests that SPL10 suppresses root regeneration by modulating auxin biosynthesis. However, we could not exclude the possibility that SPL10 regulates root regeneration via other pathways. Indeed, exogenous auxin does not fully rescue the reduced rooting defects of pSPL10:rSPL10. Moreover, our GO analysis reveals that SPL10 has a broad impact on hormone responses such as ethylene, JA, and SA. Growing evidence suggests that ethylene and SA act as rooting stimuli. For example, ethylene enhances adventitious root production in petunia (Petunia sp) cuttings and Arabidopsis hypocotyls (Druege et al., 2014; Veloccia et al., 2016). SA triggers adventitious rooting by hydrogen peroxide in mung bean (Vigna radiata) seedlings (Yang et al., 2013). Therefore, it will not be surprising that SPL10 suppresses root regeneration by attenuating ethylene and SA response in the leaf explants.
The Role of SPL10 in Adventitious Root Regeneration and Lateral Root Initiation
It has been shown that SPL10 regulates lateral root number in Arabidopsis (Yu et al., 2015b; Gao et al., 2018). Plants overexpressing miR156 produce more lateral roots whereas reducing miR156 levels leads to fewer lateral roots. The identified SPL10-AP2/ERF-ASA1 cascade is unlikely to be adopted for lateral root development because the wound-inducible AP2/ERFs are not expressed under normal growth conditions. This notion is supported by findings that, although adventitious root and lateral root converge on a similar mechanism for root primordium development, their upstream events are completely different (Verstraeten et al., 2014; Birnbaum, 2016; Druege et al., 2016; Sheng et al., 2017).
miR390, TAS3-derived trans-acting short-interfering RNAs, and AUXIN RESPONSE FACTORs form an auxin-responsive regulatory network controlling lateral root growth (Marin et al., 2010; Yoon et al., 2010; He et al., 2018). Interestingly, miR390 has a similar role in vegetative phase transition as miR156 (Adenot et al., 2006; Fahlgren et al., 2006; Hunter et al., 2006). Hence, it is possible that miR156-SPL regulates lateral root initiation through miR390-ARF pathway and its precise molecular mechanism, if true, needs further investigations.
METHODS
Plant Materials and Culture Conditions
Arabidopsis (Arabidopsis thaliana; ecotype Columbia-0) were grown at 21°C under a 16-h light/8-h dark condition with a light intensity of 80 mmol/m2/s using TLD 36W/865 and 36W/830 bulbs (Philips). The p35S:MIR156, p35S:MIM156, pSPL9:rSPL9-GR, spl9 spl15, and pDR5:GFP lines have been described in Franco-Zorrilla et al. (2007), Müller and Sheen (2008), and Wang et al. (2008, 2009). p35S:MIR156 Populus x canadensis has been described in Wang et al. (2011). The spl2 (SALK_022235), spl10 (Wiscseq_DsLox339A09), erf109 (SALK_206786), and abr1 (SALK_012151) were ordered from the Arabidopsis Biological Resource Centre. The erf109 abr1 double mutant was identified in F2 population by PCR genotyping.
Generation of Transgenic Plants
For transgenic Arabidopsis plants, the binary constructs were delivered into Agrobacterium tumefaciens GV3101 (pMP90) by the freeze-thaw method. Transgenic plants were generated by the floral dipping method (Clough and Bent, 1998) and screened with 0.05% (v/v) glufosinate (Basta; BASF) on soil, 40 mg/mL of hygromycin, or 50 mg/mL of kanamycin on 1/2 Murashige and Skoog plates.
Constructs
The oligonucleotide primers for all constructs are given in Supplemental Data Set 8.
For pSPL10:rSPL9, pSPL9:rSPL10, pSPL10:rSPL9-GR, and pSPL9:rSPL10-GR, pSPL9 and pSPL10 fragments were PCR-amplified using pSPL9:rSPL9-GR and pSPL10:rSPL10-GR as templates. ABR1 and ABR1-SRDX fragments were amplified using genomic DNA as template and then inserted into JW807. For pSPL10:3xFLAG-rSPL10, the rSPL10 fragment was inserted into pBSK:3xFLAG-X. The 3xFLAG-rSPL10 was then digested and replaced the rSPL10-GR fragment in pSPL10:rSPL10-GR. For pRIBO:ABR1-3xFLAG, ABR1 was cloned into pBSK:X-3xFLAG. The ABR1-3xFLAG was then digested and inserted into TQ120 behind RIBO promoter. pSPL9:GUS, pSPL10:GUS, and pABR1:GUS were constructed by inserting the promoters of SPL9, SPL10, and ABR1, respectively, into pBB00 in the front of the GUS coding gene.
Root Regeneration Assay
For Arabidopsis root regeneration, seeds were germinated on half strength Murashige and Skoog medium at 21°C under a 16-h light/8-h dark condition. Leaves at same developmental stage (3 mm in length) were detached and cultured on B5 medium without Suc at 21°C under a 16-h light/8-h dark condition. The rooting rate was represented by the percentage of explants with regenerated roots in a given number of explants at 20 DAC. The regenerative capacity was represented by the percentage of leaf explants with different number of regenerated roots at 12 DAC (Chen et al., 2014). For each experiment, 30 to 40 explants were used. For DEX induction experiment, 10 μΜ of DEX was added into B5 medium. Water solution with an equal volume of ethanol was used as mock. For root regeneration assay in poplar tree (Populus), the side shoots from the wild-type and p35S:MIR156 plants were detached and cultured in water. The regenerative capacity was scored after 20 d. The statistical tests are shown in Supplemental Data Set 9.
Shoot Regeneration Assay
The Arabidopsis shoot regeneration assays were performed as described by Zhang et al. (2015, 2017). The shoot regenerative capacity was represented by the number of regenerated shoots per explant. The statistical tests are shown in Supplemental Data Set 9.
Expression Analyses
Total RNA was extracted with TRIzol reagent (Invitrogen). Total RNA (1 μg) was treated with DNase I (1 unit/mL; Fermentas) and used for cDNA synthesis with oligo (dT) primer (Fermentas). The average expression levels and sd values were calculated from 2−ΔΔCt values. The RT-qPCR primers for TUB have been described in Wang et al. (2009). The oligonucleotide primers for all genes are given in Supplemental Data Set 8. The statistical tests are shown in Supplemental Data Set 9.
For hormone treatments, methyl JA (Sigma-Aldrich) was dissolved in ethanol to 50 mM and diluted d to the final concentration of 50 μΜ with water. Water solution with an equal volume of ethanol was used as mock (Mao et al., 2017). Twelve-d–old wild-type (Col-0) seedlings grown under long-day conditions were treated with mock or 50 μΜ of methyl JA. IAA was used to rescue the rooting defects in pSPL10:rSPL10 plants.
Genome Editing by CRISPR/Cas9
We used egg-cell–specific promoter-controlled CRISPR/Cas9 to generate spl10 spl11 mutant (Wang et al., 2015b). Two sgRNAs were designed and cloned into pHEE401E as described by Wang et al. (2015b). The construct was transformed into spl10 mutant by the floral dipping method. The plants carrying the mutation in SPL11 were identified by PCR.
RNA-Seq Analyses
For RNA-Seq analyses, RNA was extracted using TRIzol. Library construct and deep sequencing were performed using the HiSeq 4000 Platform (Illumina) following the manufacturer’s instructions (BGI). Three biological replicates were performed. For RNA-seq data analysis, adapter sequences and low-quality bases were filtered from the raw data using FASTP (Chen et al., 2018). All RNA-seq data were quantified by the tool Salmon (https://combine-lab.github.io/salmon/) using The Arabidopsis Information Resource 10 transcriptome as reference (Patro et al., 2017). The raw expression data were then imported into R using the tool tximport (Soneson et al., 2015). Before differential gene expression analysis, the genes whose total counts in 12 samples were <10 were prefiltered. The differential gene expression analysis was performed with the program DESeq2 (Love et al., 2014). The GO enrichment analysis was performed by the tool clusterProfiler in biological process, molecular function, and cell component, respectively (Yu et al., 2012). The RNA-seq data (BioProject PRJCA001184) were deposited with the Beijing Institute of Genomics Data Center (http://bigd.big.ac.cn). The statistical tests are shown in Supplemental Data Set 9.
Free IAA Analyses
Free IAA content measurement in Figures 2D, 2E, and 5K and Supplemental Figure 9H was performed as described in Wang et al. (2015a), with minor modifications in ultraperformance liquid chromatography-tandem mass spectrometer conditions.
IAA content measurement in Figure 2F and Supplemental Figure 5D was performed as described by Zhang et al. (2019). Briefly, the leaf explants from each sample were ground by liquid nitrogen, dissolved by 200 μL of phosphate buffered saline for 10 min on ice, and centrifuged at 12,000 rpm at 4°C for 10 min. The auxin concentration was measured with 10 μL of supernatant for each technical repeat using electrochemical detection of auxin as described by Sun et al. (2017, 2018). The statistical tests are shown in Supplemental Data Set 9.
ChIP-PCR and ChIP-seq Analyses
Briefly, the first to second leaf explants from wild type (Col-0), pSPL10:3xFLAG-rSPL10, and pRIBO:ABR1-3xFLAG were cultured on B5 medium, harvested after 4 h, and fixed according to a published protocol (Yu et al., 2013). The chromatin extract was immunoprecipitated with anti-FLAG beads (cat. nos. SA042001 and SA042005; Sigma-Aldrich). ChIP DNAs were reverse cross linked and purified with a PCR purification kit (cat. no. 28,206; Qiagen). One microliter of DNA was used for RT-qPCR analyses. The relative enrichment of 3×FLAG-rSPL10 and ABR1-3×FLAG was calculated by normalizing the amount of each immunoprecipitated fragment to input DNA, and then by normalizing the value for transgenic plants against the value for wild-type as a negative control.
For ChIP-seq data analysis, all ChIP-seq reads were first aligned to the Arabidopsis genome (The Arabidopsis Information Resource 10) with the software BowTie2 (Langmead and Salzberg, 2012). After marking duplicated reads with the program sambamba (Tarasov et al., 2015), we used the SAMtools program (Li et al., 2009) to filter duplicates (only DNA sequencing), unmapped reads, ambiguously mapped reads, and read pairs with conflicting alignments. To assess the correlation between replicates, we computed the reads coverage for consecutive bins of 10 kilobases with deepTools (Ramírez et al., 2014). Pearson correlation was then used for the analysis. All the ChIP-seq peaks were called by MACS2 with default parameters (Zhang et al., 2008). BEDTools (Quinlan and Hall, 2010) was used to identify the specific peak in pSPL10:3xFLAG-rSPL10. To find the sequence motif enriched in ChIP-seq peaks, findMotifsGenome.pl from the HOMER program was used (Heinz et al., 2010). The tool ChIPseeker was then used to retrieve the nearest genes around the peak (Yu et al., 2015a). Two biological replicates were performed. The statistical tests are shown in Supplemental Data Set 9.
GUS Staining
For pSPL9:GUS, pSPL10:GUS, and pABR1:GUS staining, explants were incubated in a GUS assay solution (100 mM of sodium phosphate buffer, 10 mM of EDTA, 0.1% [v/v] Triton X-100, and 0.1% [w/v] X-Gluc) at 37°C. The chlorophyll of the stained tissues were removed by incubating with 75% alcohol at room temperature for 12 h, and then in chloral hydrate solution (200 g of chloral hydrate, 20 g of glycerol, and 50 mL of water at 65°C for ~12 h, until tissues became transparent (Tsuge et al., 1996). The Differential Interference Contrast observations were conducted using a BX63 microscope (Olympus).
EMSA
To construct plasmid for expression of the DNA binding domain of SPL10 in Escherichia coli, the coding sequence of SPL10 was PCR-amplified and cloned into pRSF-Duet (https://www.helmholtz-muenchen.de/fileadmin/PEPF/pRSF_vectors/pRSFDuet-1_map.pdf). Oligonucleotide probes were synthesized and labeled with biotin at the 5′ end. EMSA was performed using a LightShift Chemiluminescent EMSA Kit (catalog no. 20148X; Thermo Fisher Scientific). Briefly, biotin-labeled probes were incubated in 1× binding buffer, 5 mM of MgCl2, 0.05% (v/v) NP-40, and 50 ng/μL of Poly (dI·dC) with or without proteins at room temperature for 20 min. For unlabeled probe competition, unlabeled probes were added to the reactions.
Microscopy
For confocal imaging, we used a model no. LSM880 confocal microscopic system (Zeiss) to observe GFP fluorescence.
Accession Numbers
Sequence data from this article can be found in the Arabidopsis Genome Initiative or GenBank/European Molecular Biology Laboratory databases under the following accession numbers: SPL2 (At5g43270), SPL9 (At2g42200), SPL10 (At1g27370), SPL11 (At1g27360), SPL15 (At3g57920), TUB (At5g62690), ABR1 (At5g64750), ERF109 (At4g34410), ERF115 (AT5g07310), RAP2.6L (At5g13330), RAP2.6 (At1g43160), RAP2.9 (At4g06746), CEJ1 (At3g50260), CRF3 (At5g53290), ASA1 (At5g05730), and TAA1 (At1g70560). The ChIP-seq data (BioProject PRJCA001184) were deposited with the Beijing Institute of Genomics Data Center (http://bigd.big.ac.cn).
Supplemental Data
Supplemental Figure 1. Developmental decline of rooting capacity.
Supplemental Figure 2. Generation of spl10 spl11 double mutant and phenotypic analyses.
Supplemental Figure 3. Plant phenotype and expression pattern of SPL GUS reporters.
Supplemental Figure 4. Rooting phenotype.
Supplemental Figure 5. Age and SPL10 subfamily regulate rooting capacity through auxin.
Supplemental Figure 6. RNA-seq and ChIP-seq analyses.
Supplemental Figure 7. AP2/ERF genes are the direct targets of SPL10.
Supplemental Figure 8. Time course analyses of AP2/ERFs during root regeneration.
Supplemental Figure 9. ABR1 regulates rooting capacity.
Supplemental Data Set 1. The genes upregulated after leaf detachment in wild-type leaf explants.
Supplemental Data Set 2. The genes downregulated after leaf detachment in wild-type leaf explants.
Supplemental Data Set 3. GO enrichment analysis of upregulated genes in wild-type leaf explants.
Supplemental Data Set 4. The wound-induced genes repressed by SPL10.
Supplemental Data Set 5. SPL10 binding peaks identified by ChIP-seq.
Supplemental Data Set 6. Two-hundred and thirty-three wound-induced and SPL10-repressed genes.
Supplemental Data Set 7. GO enrichment analysis of 233 candidate genes.
Supplemental Data Set 8. Oligonucleotide primer sequences.
Supplemental Data Set 9. Summary of statistical tests.
DIVE Curated Terms
The following phenotypic, genotypic, and functional terms are of significance to the work described in this paper:
Acknowledgments
We thank the Arabidopsis Biological Resource Center for seeds and Qi-Jun Chen (China Agricultural University, China) for the CRISPR/Cas9 vector; Xiao-Su Gao, Jiqin Li, Yun-Xiao He, and Shui-Ning Yang (Chinese Academy of Sciences Center for Excellence in Molecular Plant Sciences, China) for skillful technical assistance; Hongtao Liu, Lin Xu (Chinese Academy of Sciences Center for Excellence in Molecular Plant Sciences, China), and members of the J.-W. Wang lab for discussion and comments on the article; and Jijun Yan and Jinfang Chu (National Centre for Plant Gene Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, China) and Guochen Qin and Wenting Liu (Shanghai Center for Plant Stress Biology, Chinese Academy of Sciences Center for Excellence in Molecular Plant Sciences, China) for measuring free IAA contents. This work was supported National Key Research and Development Program (2016YFA0500800), the National Natural Science Foundation of China (grants 31430013, 31222029, 912173023, 31525004, and 31770399), and the Strategic Priority Research Program of the Chinese Academy of Sciences (grant XDB27030101).
AUTHOR CONTRIBUTIONS
B.-B.Y. and J.-W.W. designed the research; B.-B.Y. performed most of the experiments; G.-D.S. and Z.-G.X. analyzed the ChIP-seq and RNA-seq data; Y.P., N.B., L.S., and T.X. contributed to the auxin measurement; C.-M.Z. prepared the protein for EMSA experiment; Y.-B.M. analyzed JA-related data; B.-B.Y. and J.-W.W. analyzed the rest of the data; J.-W.W. wrote the article.
References
Author notes
The author(s) responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantcell.org) is: Jia-Wei Wang (jwwang@sippe.ac.cn).
Articles can be viewed without a subscription.




