PRH1 mediates ARF7-LBD dependent auxin signaling to regulate lateral root development in Arabidopsis thaliana

Summary

The development of lateral roots in Arabidopsis thaliana is strongly dependent on signaling directed by the AUXIN RESPONSE FACTOR7 (ARF7), which in turn activates LATERAL ORGAN BOUNDARIES DOMAIN (LBD) transcription factors (LBD16, LBD18 and LBD29). Here, the product of PRH1, a PR-1 homolog annotated previously as encoding a pathogen-responsive protein, was identified as a target of ARF7-mediated auxin signaling and also as participating in the development of lateral roots. PRH1 was shown to be strongly induced by auxin treatment, and plants lacking a functional copy of PRH1 formed fewer lateral roots. The transcription of PRH1 was controlled by the binding of both ARF7 and LBDs to its promoter region.

Keywords:

Transcription factors – Seedlings – Arabidopsis thaliana – Transcriptional control – Auxins – Lateral roots – Acetic acid – Naphthalenes

Introduction

The architecture of the root system depends on the density of lateral roots (LRs) formed along with the extent of root branching. LRs are initiated from mature pericycle cells lying adjacent to the xylem pole, referred to in Arabidopsis thaliana as the xylem pole pericycle [1, 2]. A subset of these cells, namely the founder cells, undergo a series of highly organized divisions to form an LR primordium, which eventually develops into an LR. The entire process of LR development, from its initiation to its emergence, is regulated by auxin [1, 3, 4]. The accumulation of auxin in protoxylem cells results in the priming of neighboring pericycle cells to gain founder cell identity, forming an LR pre-branch site [5]. This auxin signaling is mediated by the proteins ARF7 and ARF19, which act in the IAA14/28-dependent pathway [6, 7]. LR initiation depends on the transcriptional activation, through the involvement of the transcription factor LBD16 and LBD18, of either E2Fa or CDKA1;1 and CYCB1;1 [8, 9]. Other downstream targets such as EXP14 and EXP17 and the products of certain cell wall loosening-related genes are also either directly or indirectly regulated by LBD18 to promote the emergence of an LR [10–12]. LBD18 has been also found to control LR formation through its interaction with GIP1, which results in the transcriptional activation of EXP14 [13].

LBD transcription factors are among the most well studied downstream targets of ARF7 and ARF19 [14]. ARF7 has also been found to control LR emergence through its regulation of IDA and genes which encode leucine-rich repeat receptor-like kinases such as HAE and HSL2 [15]. These interactions take place in the overlaying tissues, thereby influencing cell wall remodeling and cell wall degradation during LR emergence [15].

Here, the product of PRH1 (At2g19990), a homolog of PR-1 which has been annotated as a pathogen responsive protein [16], was identified as a novel target of ARF7-LBDs mediated auxin signaling. We provided evidence that PRH1 acted downstream of ARF7 and LBDs to control auxin-regulated LR formation.

Results

PRH1 is a target of ARF7-mediated auxin signaling

To identify additional targets involved in ARF7-mediated auxin signaling and LR development (Fig 1A and 1B), a comparison was made between the 8-day-old wild type (WT) or arf7 mutant seedlings exposed to 10 μM naphthalene acetic acid (NAA) for 4 hours respectively (Figs 1C and S1A). Based on differential comparisons (WT and arf7 with or without NAA treatment) all of the detected genes were presented in S1 Table as the raw data and the genes changed at least 2-fold (P < 0.0001, log₂fold change >1) were selected as differentially expressed genes (DEGs) shown in S2 Table. The DEGs were selected based on the comparisons (WT-NAA VS WT-Mock, arf7-NAA VS arf7-Mock, arf7-Mock VS WT-Mock, arf7-NAA VS WT-NAA), which were induced by auxin and regulated by ARF7. To narrow down the selection of candidate genes, we further adjusted the screening criteria (P < 0.0001, log₂fold change >1.5 in WT-NAA VS WT-Mock), which reduces the DEGs to 23 (Fig 1D and S1B Fig). To further explore whether these genes are involved in LR development, we then examined the expression patterns of these 23 candidate genes in LR initiation using the Arabidopsis eFP Browser website (http://bar.utoronto.ca/efp/cgi-bin/efpWeb.cgi). After this analysis, eight of the candidate DEGs were expressed in the process of LR initiation induction, and then we focused on the PRH1 (At2g19990), which displayed the highest induction by auxin while the induction was absent in arf7, for further investigations (Fig 1E, S1B and S1C Fig).

PRH1, a homolog of PR-1 has been annotated as encoding a pathogen-responsive protein and was also named as PR-1-LIKE [16]. There are 11 members in the family in Arabidopsis, and most of them are defense-related genes (S2A Fig) [17–19]. Expression of PRH1 is not salicylic acid responsive and the biochemical function is unknown [16]. Considering that PRH1 was induced by auxin (Fig 1E), the auxin induced GH3.5 was used as a positive control and the PRH1 expression was examined through time-coursed treatment with auxin (treated with 10 μM naphthalene acetic acid, NAA). The result showed that PRH1 is transcriptionally regulated by auxin and there was a significant induction after 2 hours of auxin treatment (S2B Fig). The GUS expression analysis in 8-day-old seedlings carrying a PRH1pro:GUS transgene implied that PRH1pro:GUS is expressed in the cells overlying/surrounding the lateral root primordium (LRP), especially at endodermis or cortex cells overlying the primordium in WT background (Fig 2A; S2C Fig), and the level of GUS activity was greatly enhanced following exposure of the plants to auxin (Fig 2B). When the PRH1pro:GUS transgene was expressed in the arf7 mutant background, the PRH1pro:GUS expression was strongly suppressed, and auxin treatment only slightly induced the expression of PRH1pro:GUS (Fig 2C and 2D). Consistently, PRH1 expression was completely non-responsive to auxin treatment in arf7 according to our qPCR assays (Fig 1E). The conclusion was that the auxin-responsive PRH1 is mainly expressed in the cells overlying/surrounding the LRP and its expression is dependent on ARF7-mediated auxin signaling.

PRH1 acts downstream of ARF7 to regulate LR development

To address if auxin-induced PRH1 expression is involved into LR development, we examined the LR phenotypes in both prh1 mutants and the PRH1 over-expression lines (Fig 3A and 3B; S3B Fig). The third allele of prh1 (phr1-3) generated via Crispr-cas9 system [20] was also used for the LR phenotype analysis (S3A Fig). The primordia numbers of 8-day-old seedlings at various stages were counted of all these three prh1 mutant alleles. The results showed that primordium development was affected. The numbers at stage I and Ⅳ were increased and at the later periods (stage Ⅵ, Ⅶ and emerged) were reduced (Fig 3A). Examination of LR development in prh1 mutants showed that the absence of a functional copy of PRH1 inhibited the formation of emerged LRs (Fig 3E). The expression level of PRH1 in its T-DNA insertion mutants and the over-expression lines were all verified (S3C Fig). The weak phenotype of prh1-2 in the emerged stage might be result from the down-regulation of PRH1 in this allele (Fig 3A; S3C Fig). Moreover the PRH1pro:PRH1 construct was transformed in the prh1-1 mutant, the expression of PRH1 could complement the LR defects in prh1-1 (S4 Fig). Though LR development was not greatly affected by the over-expression of PRH1 in a WT background (S3B Fig), over-expression of PRH1 could significantly increase the LR numbers of arf7 mutant plants in the presence of exogenous NAA treatment (Fig 3C–3E).

To further confirm the role of PRH1 in LR formation, we did LR analysis of prh1 mutants through a gravistimulation-based bioassay [21–23]. After synchronizing LR formation with a gravistimulus for 18 h, we observed a delayed transition from stage I to stage II of LR development which was shown by the increased LR I (lateral root stage I) ratio and decreased LR II (lateral root stage II) ratio in prh1 mutants (Fig 3B). Then we examined the phenotype of mutant seedlings after a gravistimulus for 42 h, and the lateral root initiation rates reduced both in stage VII and emerged ones, but increased in the earlier stages (stage IV or V) (Fig 3B). In summary, PRH1 acts downstream of ARF7 to regulate LR development.

ARF7 directly regulates PRH1 transcription

The observation that the abundance of PRH1 transcript was reduced in the arf7 mutant (Fig 1D) suggested that the gene was transcriptionally regulated by the auxin response factor ARF7. The gene’s promoter sequence harbors the known auxin response elements (AuxREs) TGTC (Fig 4A). A transient expression assay carried out in A. thaliana leaf protoplasts showed that ARF7 was able to activate a construct comprising the PRH1 promoter fused with LUC. The co-expression of PRH1pro:LUC and 35S::ARF7 had a large positive effect on the intensity of the luminescence signal (Fig 4B), consistent with the idea that ARF7 is able to activate the PRH1pro:LUC construct. A ChIP-qPCR assay confirmed the association of ARF7 with the AuxRE elements present on the PRH1 promoter (Fig 4A and 4C), and a yeast one hybrid assay showed that ARF7 was able to bind to the PRH1 promoter (Fig 4D). The conclusion was that the effect of ARF7 binding to the PRH1 promoter was regulating the gene’s transcription, thereby influencing LR development.

Auxin-induced PRH1 expression is regulated by LBD29

Members of the LBD family are known to function downstream of ARF7 to control LR initiation and emergence [14, 24–27]. When the transcription level of PRH1 was investigated in the three single lbd mutants lbd16 (Salk_040739), lbd18 (Salk_038125) and lbd29 (Salk_071133) combined with or without exogenous auxin treatment, the observation was that the gene was strongly repressed in the single mutants. Though the auxin induced PRH1 expression was normal in lbd16 or lbd18, the induction was strongly reduced in lbd29, indicating LBD29 plays an important role in auxin induced PRH1 expression (S5A Fig). Consistently, the over-expression of LBD29 strongly up-regulates PRH1 expression (S5B Fig). In addition, we also observed the up-regulation of PRH1 expression in LBD16, LBD18 or LBD29 over-expression lines (S5B Fig). These results suggest that LBDs regulate the transcription of PRH1, auxin induced expression of PRH1 is dependent on LBD29.

Since the auxin influx carrier LIKE-AUXIN3 (LAX3) has been reported to act downstream of ARF7-LBD29 auxin signaling module to facilitate LR emergence [27, 28], we also investigated the relationship between LAX3 and PRH1 through qPCR analysis. The results showed that the absence of PRH1 did not affect the expression of LAX3 with or without auxin treatment (S5C Fig). However, in lax3 which has defects in auxin influx, the auxin-induced PRH1 expression was reduced (S5D Fig).

LBDs regulate PRH1 by binding to its promoter

The transcriptional response of PRH1 to the loss-of-function and over-expression of certain LBDs implied that these transcription factors can also regulate PRH1. The PRH1 promoter has the typical LBD protein binding elements including CGGCG and CACGTG (S6A Fig). Although the yeast one hybrid assay produced no direct evidence for the binding of LBD proteins to the PRH1 promoter (S6B Fig), but the ChIP-qPCR assay did support the association of LBD16, LBD18 and LBD29 with the PRH1 promoter (Fig 5A). Thus a transient expression assay in A. thaliana leaf protoplasts was used to investigate whether LBD16, LBD18 and/or LBD29 was able to influence the expression of the PRH1pro:LUC transgene. The outcome of co-expressing this construct with each of 35S::LBD16, 35S::LBD18 or 35S::LBD29 was to enhance the expression of PRH1pro:LUC (Fig 5B). When a pair of effector transgenes (either 35S::LBD16 plus 35S::LBD18, 35S::LBD18 plus 35S::LBD29 or 35S::LBD16 plus 35S::LBD29) was introduced, the PRH1pro:LUC expression was even more strongly activated (Fig 5B). The overall conclusion was that the LBDs regulate PRH1 expression by directly binding to its promoter.

LBD regulated LR development is partially dependent on PRH1

To further study if LBDs regulated PRH1 transcription is involved into LBD-mediated LR development, we crossed 35S::PRH1 with lbd16, lbd18 or lbd29 mutants and analyzed LR phenotypes of the generated double mutants. Compared to lbd16, lbd18, lbd29 single mutant plants, the generated 35S::PRH1/lbd16, 35S::PRH1 lbd18 or 35S::PRH1/lbd29 double mutants have increased LR numbers, especially in the presence of auxin treatment (Fig 6), indicating that PRH1 is involved into LBD-mediated LR development.

To determine the subcellular localization of PRH1, the cytosolic marker ADH-RFP [29] or plasma membrane maker ZmCRN-mCherry [30] were coexpressed with PRH1-GFP in N. benthamiana leaf cells (S7 Fig). The green signals from PRH1-GFP were merged with the red signals from ADH-RFP, indicating that PRH1 is localized in the cytosol. Given the LR phenotype of prh1 mutants and its expression pattern, the expression levels of some cell wall formation related genes were detected in prh1-1. We found some EXPANSIN genes, including EXP3 EXP8 EXP11 and EXP14 which were all involved in cell wall loosening [12, 31], were down-regulated in prh1 (S8 Fig). Therefore, PRH1 might regulate LR emergence by influencing the EXPANSIN gene expression.

According to this study together with the previous reports, a proposed model was given in Fig 7. PRH1 is involved in ARF7-LBD dependent auxin signaling pathway regulating lateral root development. PRH1 acts as a downstream gene of both ARF7 and LBDs by the direct transcriptional regulation and then affects the expression of some EXPANSIN genes. Meanwhile LAX3 acts downstream of LBD29 to control IAA accumulation in the epidermits and cortex [32], and thus influences auxin-induced PRH1 expression indirectly.

Discussion

The product of the auxin-induced gene PRH1 controls LR formation

The role of auxin in LR development has been extensively investigated. During the initiation of LR, auxin signals are transduced through the two distinct AUX/IAA-ARFs modules IAA14-ARF7/ARF19 and IAA12-ARF5 [6, 33–36]. This process establishes LR founder cell identity, which is the initial step towards the formation of a LR. Auxin is also involved in LR emergence, since it controls the activation of IDA, HAE and HSL2 which encode cell wall remodeling enzymes and of genes encoding various expansions in the endodermal cells which lie above the LR primordia [11, 15, 37]. Members of the LBD gene family encode transcription factors, characterized by their LOB domain [38], which act downstream of ARF7-mediated auxin signaling to control LR formation [13, 24–26].

Here, the PR-1 homolog gene PRH1 was identified as a novel target of ARF7-mediated auxin signaling. Otherwise we have compared our RNA-Seq data to the previous datasets of lateral root induction, recommend in Lewis et al [39] and Himanen et al [40]. Consistently, all these studies showed that the expression of PRH1 was induced by auxin and the induction depended on ARF7. Although the over-expression of PRH1 in a WT background was not associated with any LR phenotype, its over-expression in arf7, lbd16, 1bd18 or lbd29 mutant backgrounds led to a partial rescue of the compromised LR phenotype, especially under auxin treatment (Figs 3 and 6). The more striking rescued LR formation by overexpression of PRH1 in lbd mutants than that in arf7 suggests additional factors other than LBDs act downstream of ARF7. The PRH1 promoter directs its product primarily to the cells overlying/surrounding the lateral root primordium (as shown by the site of GUS expression in plants expressing a PRH1pro:GUS transgene), which is consistent with the proposed role of PRH1 in LR formation. Auxin treatment strongly induced the expression of the PRH1pro:GUS transgene in a WT (Fig 2A and 2B), but not in an arf7 background (Fig 2C and 2D), which implied that the auxin-induced expression of PRH1 must depend on the presence of a functional copy of ARF7. LBD proteins, especially LBD29, are also involved in the induction by auxin of PRH1. Experimentally, it was shown that both ARF7 and LBDs were able to bind to the PRH1 promoter (Figs 4 and 5). As some of the LBD-binding and AuxRE elements close to each other in PRH1 promoter (Fig 5A), it will be of future interest to investigate how ARF7 coordinates with LBDs to mediate auxin signaling and therefore to control the expression of PRH1. In addition, the auxin influx carrier LIKE-AUXIN3 (LAX3) which acts downstream of ARF7-LBD29 auxin signaling module to control LR emergence [27, 28], might be involved into auxin-induced PRH1 expression through the regulation of local auxin transport and local auxin accumulations.

The homo- and heterodimerization of LBD proteins are thought to be an important determinant of LR formation [26, 41]. The AtbZIP59 transcription factor has been reported to form complexes with LBDs to direct auxin-induced callus formation and LR development [42]. All these results suggest that homodimerization and heterodimerization of signaling components are important in LR formation. Here, it was shown that PRH1 was mainly expressed in the cells overlying the primordium (Fig 2) and there was a LR emergence defect in the prh1 mutants. Several of the EXPANSIN genes were down-regulated in prh1-1(S8 Fig), indicating that PRH1 might regulate LR emergence through influencing cell wall loosening. Whether PRH1 regulates LR formation through a direct interaction with LBDs, AtZIP59 and other factors which have been shown to control LR development, is still an open question for the future.

Materials and methods

Plant materials and growing conditions

A. thaliana ecotype Col-0 was used as the wild type (WT), and all mutants and transgenes were constructed in a Col-0 background. The mutants comprised a set of T-DNA insert lines obtained from the Arabidopsis Biological Resource Center (ABRC, Columbus, OH, USA), namely prh1-1 (GK_626H08), prh1-2 (Salk_014249), lbd16 (Salk_040739), lbd18 (Salk_038125) [43] and lbd29 (Salk_071133) [27], along with arf7 (Salk_040394) [33]; lax3 [28]; the transgenic lines were Super::LBD29-GFP [44], 35S::MYC-LBD16, 35S::MYC-LBD18, 35S::PRH1, 35S::PRH1-GFP (PRH1-GFP) and PRH1pro:GUS. Seeds were surface-sterilized by fumigation with chlorine gas for 4 hours, and then plated on solidified half strength Murashige and Skoog (MS) medium. All of the seeds were held for two days at 4°C and finally grown in a greenhouse delivering a 16 h photoperiod with a constant temperature of 22°C.

Transgene constructs and the generation of transgenic lines

The 1,692 nt region upstream of the PRH1 start cordon was used as the PRH1 promoter, and was inserted, along with the PRH1 coding sequence into a pENTR^TM/SD/D-TOPO^TM plasmid (ThermoFisher). The construct was subsequently recombined with several Gateway destination vectors, namely PK7WG2 [45] for over-expression without any tags, PB7WGF2 for over-expression with an C-GFP tag, PGWB18 for over-expression with an N-Myc tag, PBI221 for over-expression, PKGWFS7.1 [45] for promoter analysis with C-GFP/GUS. Transgenesis was effected using the floral dip method [46], employing transgene constructs harbored by Agrobacterium tumefaciens strain GV3101.

Microscopy

The PRH1 subcellular localization was detected in tobacco leaf cells. The Agrobacterium tumefaciens strain GV3101 containing 35S::PRH1-GFP or the cytosolic marker ADH-RFP [29] or plasma membrane maker ZmCRN-mCherry [30] were incubated shaking until the OD₆₀₀ reaching 0.6. Centrifuge the cells at 4000g for 5 minutes, then discard the supernatant and resuspend the pellet with induction buffer (0.01mol/L MES, 0.01Mol/L MgCl₂, 150 μmol/L acetosyringone) to find OD = 0.4. Mix two kinds of agrobacterium up with a same volume and inject the mixture into 3 week old tobacco leaves. The transfected plants grown in the greenhouse for at least 48 h at 28°C before GFP and RFP caning. The fluorescent signals were visualized by using the LSM880 laser scanning confocal microscope (Zeiss). ZmCRN-mCherry and ADH-RFP signals were detected by using 561 nm laser excitation and 580–675 nm emission, PRH1-GFP signals were detected by using 561 nm laser excitation and 500–530 nm emission.

The 8-d-old seedlings of PRH1pro:GUS were used for GUS staining. The GUS staining assay followed Liu et al. [47]. GUS signal was photographed by using an Olympus BX53 microscope supported with an Olympus DP72 digital camera.

Phenotype analysis

Eight-day-old seedlings were photographed using a Epson Perfection V800 scanner for the whole seedling images and the root length was measured by Image J software. To observe lateral root primordia clearly the seedlings were treated with a serials of solutions (1. 4% HCL, 20% methanol at 57°C for 15 minutes. 2. 7% NaOH, 60% ethanol at room temperature for 15 minutes. 3. 60% ethanol at room temperature for 10 minutes. 4. 40% ethanol at room temperature for 10 minutes. 5. 20% ethanol at room temperature for 10 minutes. 6. 10% ethanol at room temperature for 10 minutes. 7. 25% glycerol, 5% ethanol at room temperature for 10 minutes) before mounted in 50% glycerol on slides. Then the lateral root primordia and the emerged lateral roots were calculated using the Olympus BX53 microscope. The differential stages of lateral roots development have been described by Malamy and Benfey [1].

RNA isolation and quantitative real time PCR (qPCR) analysis

RNA was extracted from the roots of eight-day-old seedlings (some of which had been treated with 10 μM naphthalene acetic acid (NAA), obtained from Sigma (St. Louis, MO, USA)) using a TIANGEN Biotech (Beijing) Co. RNA simple Total RNA kit (www.tiangen.com), following the manufacturer’s protocol. A 2 μg aliquot of the RNA was used as the template for synthesizing the first cDNA strand, using a Fast Quant RT Kit (TIANGEN Biotech (Beijing) Co.). Subsequent qPCRs were processed on an MyiQ Real-time PCR Detection System (Bio-Rad, Hercules, CA, USA), using the Super Real PreMix Plus SYBR Green reagent (TIANGEN Biotech (Beijing) Co.). Each sample was represented by three biological replicates, and each biological replicate by three technical replicates. The AtACTIN2 (At3g18780) sequence was used as the reference. Details of the relevant primer sequences are given in S3 Table.

RNA-Seq analysis

RNA submitted for RNA-seq analysis was isolated from the roots of eight-day-old WT and arf7 seedlings (treated with the 10 μM naphthalene acetic acid (NAA) for 4 hours and the non-treated samples), using a RNeasy Mini Kit (Qiagen, Hilden, Germany). The RNA-seq analysis were processed following the methods given by Lv et al. [48]

Chromatin immunoprecipitation coupled to quantitative PCR (ChIP-qPCR) analysis

DNA harvested from eight-day-old seedlings of WT and transgenic lines harboring either 35S::MYC-ARF7, 35S::MYC-LBD16, 35S::MYC-LBD18 (ARF7, LBD16 and LBD18 were cloned into PGWB18) or Super::LBD29-GFP [44, 49] then they were cross-linked in a solution containing 1% v/v formaldehyde. The ChIP procedure followed that given by Gendrel A-V et al. [50], and employed antibodies recognizing MYC and GFP (Abcam, UK). The immune-precipitated DNA and total input DNA were analyzed by qPCR. The enrichment in WT and transgenic plants was measured as a ratio of the bound promoter fragments over total input respectively. Details of the relevant primers are given in S3 Table.

Yeast one-hybrid (Y1H) assays

Y1H assays were carried out using the Matchmaker Gold yeast one hybrid Library Screening System (Biosciences Clontech, Palo Alto, CA, USA). The PRH1 promoter was inserted into the pAbAi vector, and the bait vector (PRH1_pro-pAbAi) was linearized, restricted by BbsI, and introduced into the Y1H Goldeast strain along with the prey vector AD-ARF7. Transformants were grown on SD/-Leu-Ura dropout plates containing different concentrations of aureobasidin A (CAT#630499, Clontech, USA). The primers used for generating the various constructs are listed in S3 Table.

Dual-luciferase transient expression assay in A. thaliana protoplasts

Protoplasts were isolated from three week old WT leaves following Yoo et al. [51]. The ligation of the PRH1 promoter to pGreen0800-LUC was realized by enzyme digestion and ligation [52]. The primers used for this procedure are listed in S3 Table. The PRH1pro:LUC reporter plasmid was transferred into Arabidopsis protoplast cells with one or more effector plasmids (ARF7, LBD16, LBD18, LBD29 –PBI221); the empty vector PBI221 was used as the negative control. The dual-luciferase assay kit (Promega, USA) allows for the quantification of both firefly and renilla luciferase activity. Signals were detected with a Synergy 2 multimode micro-plate (Centro LB960, Berthold, Germany).

Generation of prh1-3 using CRISPR technology

The experimental method was described by Yan et al [20]. The Arabidopsis plants Columbia-0 ecotype were used as wild type (WT) and were transformed with the CRISPR construct by floral dipping. The prh1-3 plants were identified at the T2 stage. The target sequence is presented in S3 Table.

Statistical analysis

Tests were applied to the data to check for normality and homogeneity before attempting an analysis of the variation. The data were analyzed using routines implemented in Prism 5 software (GraphPad Software). The Students’ t-test was applied to compare pairs of means. Values have been presented in the form means±SE; the significance thresholds were 0.05 (*), 0.01 (**) and 0.001 (***). Comparisons between multiple means were made using a one-way analysis of variance and Tukey’s test.

Accession numbers

Sequence data referred to in this paper have been deposited in both the Arabidopsis Genome Initiative and the GenBank/EMBL databases under the accession numbers: ARF7 (At5g20730, NP_851046), PRH1 (At2g19990, NP_179589), LBD16 (At2g42430, NP_565973), LBD18 (At2g45420, NP_850436), LBD29 (At3g58190, NP_191378), LAX3 (At1g77690, NP_177892). The RNA-seq data are available in the Gene Expression Omnibus database under accession number GSE122355.

Supporting information

S1 Fig [a]
Transcriptomic analysis of WT and roots before and after auxin treatment.

S2 Fig [a]
responses to auxin treatment.

S3 Fig [a]
Lateral roots (LRs) development phenotype analysis and the transcription level of in wild type (WT), mutants and its over-expression lines.

S4 Fig [a]
PRH1 could complement the defects of Lateral roots (LRs) formation in mutants.

S5 Fig [a]
The transcription of is regulated by LBD16, LBD18, LBD29 and LAX3.

S6 Fig [a]
The interaction between LBDs and promoter in yeast one-hybrid.

S7 Fig [a]
The subcellular localization of PRH1.

S8 Fig [tif]
Expression analysis of cell wall formation related genes in .

S1 Table [xlsx]
The raw data of gene expression based on the comparison between the four datasets.

S2 Table [xlsx]
Gene differential expression filter data according to the comparison between the four datasets.

S3 Table [xlsx]
Primers used in this study.