Brassinosteroid signaling directs formative cell divisions and protophloem differentiation in Arabidopsis root meristems

ABSTRACT Brassinosteroids (BRs) trigger an intracellular signaling cascade through its receptors BR INSENSITIVE 1 (BRI1), BRI1-LIKE 1 (BRL1) and BRL3. Recent studies suggest that BR-independent inputs related to vascular differentiation, for instance root protophloem development, modulate downstream BR signaling components. Here, we report that protophloem sieve element differentiation is indeed impaired in bri1 brl1 brl3 mutants, although this effect might not be mediated by canonical downstream BR signaling components. We also found that their small meristem size is entirely explained by reduced cell elongation, which is, however, accompanied by supernumerary formative cell divisions in the radial dimension. Thus, reduced cell expansion in conjunction with growth retardation, because of the need to accommodate supernumerary formative divisions, can account for the overall short root phenotype of BR signaling mutants. Tissue-specific re-addition of BRI1 activity partially rescued subsets of these defects through partly cell-autonomous, partly non-cell-autonomous effects. However, protophloem-specific BRI1 expression essentially rescued all major bri1 brl1 brl3 root meristem phenotypes. Our data suggest that BR perception in the protophloem is sufficient to systemically convey BR action in the root meristem context. Highlighted article: In addition to controlling cell expansion, brassinosteroid signaling acts cell-autonomously and non-cell-autonomously to promote protophloem differentiation and restrict formative divisions in root meristems.


INTRODUCTION
Brassinosteroids (BRs) were discovered as plant cell elongation and division stimulants. Although the principally active BR, brassinolide, was isolated in 1979 (Grove et al., 1979), the BR signaling pathway was characterized much later, through genetic approaches in Arabidopsis thaliana. The key discovery was the isolation of BRASSINOSTEROID INSENSITIVE 1 (BRI1), a receptor kinase that triggers an intracellular signaling cascade upon extracellular BR perception (Li and Chory, 1997). Downstream BR signaling involves the GSK3/SHAGGY-LIKE kinase BRASSINOSTEROID-INSENSITIVE 2 (BIN2) and its substrates, the homologous transcription factors BRASSINAZOLE-RESISTANT 1 (BZR1) and BRI1-EMS-SUPPRESSOR 1 (BES1) (Li and Nam, 2002;Yin et al., 2002;Zhu et al., 2013). BR perception inactivates BIN2, thereby enabling BZR1 and BES1 to promote a transcriptional response. Among the three Arabidopsis BRI1 homologs, only BRI1-LIKE 1 (BRL1) and BRL3 are functional BR receptors (Cano-Delgado et al., 2004). brl1 or brl3 single or double loss-of-function mutants lack discernible phenotypes; however, the severe dwarfism of bri1 mutants is enhanced in bri1 brl1 brl3 triple mutants. Interestingly, tissue-specific expression of BRI1 in the epidermal cell layer largely recues bri1 dwarfism (Savaldi-Goldstein et al., 2007), but cannot complement defects in internal tissues, for example altered vascular patterning.
The role of BRs is best understood in hypocotyl elongation. Here, BR signaling is required for optimal cell expansion, which involves synergistic interaction with other hormones, for instance auxin (Nemhauser et al., 2004). BRs are also essential for primary root growth, because loss-of-function biosynthetic as well as signaling mutants have short roots. Whereas the bri1 shoot phenotype appears to primarily result from reduced cell expansion (Savaldi-Goldstein et al., 2007), the situation is more complex in the root, where cell proliferation and elongation are deeply intertwined (Scacchi et al., 2010). In root development, BR impact on cell proliferation is seemingly more prominent (Gonzalez-Garcia et al., 2011;Hacham et al., 2011). For example, it was reported that root meristem size is reduced in bri1 mutants as indicated by the number of dividing cells along cortex cell files. Although such reduction was not observed in a BR biosynthesis mutant (Chaiwanon and Wang, 2015), all mutants displayed decreased cell elongation, as indicated by, for example, mature cortex cell length. Phenotypic discrepancies could be explained by redundancies or crossregulations, for instance between BR receptors, or by alternative biosynthetic bypasses (Ohnishi et al., 2006). Moreover, analysis of tissue-specific BRI1 re-addition into receptor mutants led to the conclusion that epidermal BRI1 activity promotes root meristem growth and affects inner tissues through unknown non-cellautonomous signals (Hacham et al., 2011;Vragovićet al., 2015).
Recently, it was found that GSK3 activity, including BIN2, could be modulated by a non-BR receptor system that controls xylem vessel differentiation (Cho et al., 2014;Kondo et al., 2014). Moreover, it was reported that BIN2 is a crucial interactor of OCTOPUS (OPS) (Anne et al., 2015), a quantitative positive master regulator of protophloem differentiation in the root meristem (Rodriguez-Villalon et al., 2014;Truernit et al., 2012). Loss-offunction ops mutants display impaired protophloem sieve element differentiation, which is associated with systemic effects, such as reduced auxin activity throughout the meristem and strongly reduced root growth. Loss-of-function mutants in BREVIS RADIX (BRX), another positive regulator of protophloem differentiation, have a similar phenotype (Rodriguez-Villalon et al., 2015). Interestingly, it was reported that brx root growth could be partially rescued by BR application (Gujas et al., 2012;Mouchel et al., 2006), and ops protophloem defects could be partially suppressed by BR pathway activation (Anne et al., 2015;De Rybel et al., 2009). These observations motivated us to assess the role of BR signaling in root protophloem differentiation.

BR pathway activation cannot substitute for known protophloem differentiation factors
First, we examined whether, or to what degree, the BRX or OPS genes act through brassinosteroid signaling components. Similar to the hypocotyl, it has been suggested that threshold BR activity limits auxin activity in the root, because BR application rescues systemically reduced auxin response in brx root meristems (Gujas et al., 2012;Mouchel et al., 2006). Consistent with this, BR treatment also alleviated reduced auxin response in ops meristems (Rodriguez-Villalon et al., 2014) as indicated by the inverse auxin signaling reporter DII-VENUS (Santuari et al., 2011), whereas wild-type root meristems showed little detectable change (Fig. 1A). However, compared with brx, the effect of BR treatment in ops was less pronounced and an impact on root growth hardly detectable (Fig. 1B). To test whether BR treatment also affected brx or ops protophloem differentiation, we quantified cells that apparently failed to enter the sieve element differentiation program. These so-called 'gap' cells occur in the differentiation zone of the protophloem sieve element strands and are easily distinguished by their reduced propidium iodide cell wall staining (Scacchi et al., 2010;Truernit et al., 2012). Interestingly, BR-treated ops, but not brx mutants, showed statistically significant reduction in gap cell frequency (Fig. 1C). By contrast, no rescue was observed with respect to root growth or gap cell frequency in either genotype upon treatment with bikinin ( Fig. 1D,E), a pharmacological GSK3 (and BIN2) inhibitor (De Rybel et al., 2009). Thus, BR pathway activation could partially rescue distinct aspects of the brx or ops phenotypes. Consistent with the notion that the observed effects reflect parallel action of BR signaling and BRX, combination of a loss-of-function allele in the main BR receptor, BRI1, with a brx null allele led to an additive root phenotype in homozygous brx bri1 double mutants (Fig. S1A). Matching this observation, inspection by confocal microscopy also revealed an even more severe root meristem phenotype than in either single mutant (Fig. S1B). This included development of the protophloem, which displayed hardly any properly differentiating sieve elements. This phenotype was so severe that gap quantification became practically unfeasible. Moreover, in many brx bri1 roots (∼40%) only a single protophloem strand could be detected. bri1 brl1 brl3 triple mutants display protophloem differentiation defects To investigate further a possible involvement of the BR signaling pathway in protophloem differentiation, we next monitored the phenotype of bri1 mutants in more detail. However, although bri1 mutants displayed strongly reduced root growth, no gap cells were observed (n=44), consistent with previous reports (Anne et al., 2015). This could mean that phloem phenotypes were masked by redundancy or feedback systems. To avoid such ambiguous scenarios, we investigated bri1 brl1 brl3 triple mutants (referred to hereafter as 'triple mutant'), in which BR perception is completely shut down and therefore phenotypes are uncoupled from regulatory feedbacks on BR signal transduction or biosynthesis (Cano-Delgado et al., 2004). These experiments revealed that gap cells do indeed occur in bri1 brl1 brl3 triple mutants ( Fig. 2A), although at considerably lower frequency than in brx or ops. By contrast, brl1 brl3 double mutant protophloem were wild type in appearance (n=26), corroborating the unequally redundant action of BRL1 and BRL3 (Briggs et al., 2006), as well as the nevertheless dominant role of BRI1 in protophloem development. As expected, neither triple mutant protophloem defects nor overall root growth could be rescued by BR application. However, root growth, but not protophloem differentiation, was partially rescued by bikinin treatment (Fig. 2B,C).
Systemic propagation of local BR signaling via unknown pathways has been suggested as an explanation for non-cellautonomous rescue of BR receptor mutants by epidermis-specific re-addition of BRI1 activity (Hacham et al., 2011;Savaldi-Goldstein et al., 2007). With respect to the root, opposing, action sitedependent effects of BR signaling on root meristem size have been reported (Hacham et al., 2011;Vragovićet al., 2015). Here, we analyzed whether protophloem differentiation responds to proposed non-cell-autonomous cues by investigating triple mutants that expressed functional BRI-GFP fusions under control of either the GLABRA 2 (GL2) or the SHORT ROOT (SHR) promoter (hereafter GL2::BRI1-GFP triple and SHR::BRI1-GFP triple , respectively) (Vragovićet al., 2015) ( Fig. 2D). Whereas GL2 confers expression in epidermal non-root-hair cell files, SHR drives expression throughout the stele. In both cases, the triple mutant short root phenotype was partially rescued (Fig. 2E), but protophloem differentiation defects were not (Fig. 2F). Although this was expected for GL2::BRI1-GFP triple plants, it was surprising for the SHR::BRI1-GFP triple plants. However, we found that SHR:: BRI1-GFP was not expressed in the phloem poles ( Fig. S1C), unlike BRI1-GFP expressed under control of its native promoter (Fig.  S1D). Thus, the triple mutant protophloem defects could not be rescued by re-addition of BRI1-GFP to immediately neighboring cell files, arguing for a possibly cell-autonomous contribution of BR signaling to protophloem development.
The small root meristem of bri1 brl1 brl3 triple mutants can be entirely explained by reduced cell elongation The absence of protophloem defects in bri1 single mutants not only suggested BR receptor redundancy in protophloem development, but also that impaired protophloem differentiation is not the primary cause for the short root of bri1 or the triple mutant. In part, this phenotype has been attributed to reduced meristem size (Fig. S1E), as indicated by the number of meristematic cells along root cell files (Gonzalez-Garcia et al., 2011;Hacham et al., 2011;Vragovićet al., 2015). However, given the overall morphological irregularities in the triple mutant or its GL2::BRI1-GFP triple and SHR::BRI1-GFP triple derivatives ( Fig. 3A), assigning the border between cell proliferation and elongation could become somewhat subjective. Therefore, we measured the progression of cell expansion along meristematic cortex files instead. Plots of cumulative cell length along the cortex files indicated that cell elongation is reduced in the triple mutant at all positions. However, the corresponding wild-type and triple mutant curves were highly correlated (Pearson's r=0.96) (Fig. 3B), indicating that the small bri1 brl1 brl3 meristem can be entirely explained by reduced cell elongation, similar to the dwarf4 BR biosynthetic mutant (Chaiwanon and Wang, 2015). In wild type, cell divisions were no longer observed once cells elongated above 20 µm. In the triple mutant or its derivatives, 20 µm length was reached about ten cells later than in wild type (Fig. 3C), and cell files contained about eight more cells at 400 µm ( Fig. 3D). However, because cell elongation was generally reduced, cell divisions in the mutant lines apparently ceased earlier, before the 20 µm mark. Moreover, correlation between GL2::BRI1-GFP triple or SHR::BRI1-GFP triple and wild type was still high (r=0.80 and 0.88, respectively), and even higher when compared with bri1 brl1 brl3 (r=0.84 and 0.89, respectively), consistent with no statistically significant difference between GL2::BRI1-GFP triple and SHR::BRI1-GFP triple and the triple mutant in other parameters (Fig. 3C,D).
Interestingly, statistically significant, but small partial rescue of reduced mature cell size, as exemplified by cortex (Fig. 3E) and epidermis ( Fig. 3F) cells, was observed in SHR::BRI1-GFP triple plants (Fig. 3E,F). By contrast, GL2::BRI1-GFP triple plants displayed even further reduced mature cell length, as previously reported (Fridman et al., 2014). However, reduced root growth was still weakly complemented by GL2::BRI1-GFP, even after a longer growth period (Fig. S1F). Notably, however, the phenotype of these plants was difficult to assess because their roots twisted, which could result from the fact that GL2::BRI1-GFP is only expressed in non-root-hair cell files. The increased ratio of non-hair cell length to hair cell length might reflect this (Fig. 3G). Transverse sections showed that BRI1-GFP-expressing epidermal cells also expanded considerably more radially than their immediate hair cell neighbors (Fig. 4A), confirming earlier reports (Fridman et al., 2014), which is likely to create mechanical strain between them.
Cell division rate in the radial dimension reacts to cellautonomous as well as non-cell-autonomous cues Consistent with previous reports that BR signaling promotes columella stem cell divisions (Chaiwanon and Wang, 2015;Gonzalez-Garcia et al., 2011), bri1 brl1 brl3 mutants displayed slightly reduced columella cell layers, which was not rescued in the derivative lines (Fig. 3H). By contrast, the number of formative cell divisions was markedly increased in the radial dimension across all cell layers except the endodermis, giving rise to additional cell files ( Fig. 4A-C). Moreover, extra periclinal cortex and endodermis divisions were observed occasionally in the triple mutant, but not in wild type (Fig. 4D). This phenotype was already observed in a bri1 single mutant, although it was not observed as often as in the triple mutant (Fig. S2A). Interestingly, our observations fit with previously reported increased cyclin B1 activity in the tip of bri1 mutant root meristems (Gonzalez-Garcia et al., 2011) (with the caveat that this is not necessarily indicative of more cell divisions). Together with the reduced cell expansion described above, this observation would also explain the overall strongly reduced root growth of bri1 brl1 brl3 mutants, because accommodation of supernumerary formative divisions could presumably slow down overall root elongation.
Surprisingly, the GL2::BRI1-GFP and SHR::BRI1-GFP transgenes modulated the radial increase in cell files in a unique, opposite fashion. Whereas GL2::BRI1-GFP triple plants basically complemented, and sometimes even over-compensated the excessive radial cell proliferation, the triple mutant phenotype was further enhanced in SHR::BRI1-GFP triple plants ( Fig. 4A-C). Thus, the cell division pattern appeared to accommodate local variation in BR signaling in a partly cell-autonomous, partly non-cellautonomous fashion. The observations also suggest that overall root growth was partially rescued in both GL2::BRI1-GFP triple and SHR::BRI1-GFP triple plants for different reasons. Whereas the excess formative cell divisions, but not the reduced cell elongation, were partially compensated in GL2::BRI1-GFP triple plants, the inverse applied to SHR::BRI1-GFP triple plants.
BRI1 activity in the developing protophloem is sufficient to rescue all observed major root phenotypes of bri1 brl1 brl3 triple mutants Interestingly, in the F1 derived from crosses between GL2::BRI1-GFP triple and SHR::BRI1-GFP triple plants, the formative cell division phenotypes of the two parental lines largely compensated each other (Fig. S2B). At the same time, both mature cortex cell length and overall root growth were restored to nearly wild-type levels (Fig. S2,D). However, gap cells still occurred in these plants and, although their frequency was reduced, this was only borderline statistically significant compared with the parental lines (Fig. S2E). To evaluate directly the contribution of BR signaling to protophloem development, we therefore created transgenic lines in which we expressed a BRI1-CITRINE fusion protein under control of different, increasingly protophloem-specific promoters. The MEMBRANE-ASSOCIATED KINASE REGULATOR 5 (MAKR5) promoter confers expression in the phloem poles in the early meristem, and throughout the root vasculature later on (Kang and Hardtke, 2016) (Fig. 5A). By contrast, the BARELY ANY MERISTEM 3 (BAM3) promoter drives expression at the phloem pole from early on, with its strongest activity in the developing sieve elements and possibly some weak activity in companion cells (Depuydt et al., 2013;Rodriguez-Villalon et al., 2014) (Fig. 5B). Finally, the COTYLEDON VASCULAR PATTERN 2 (CVP2) promoter is very specific for the developing sieve element cell file (Rodriguez-Villalon et al., 2015) (Fig. 5C). Analysis of lines that carried the BAM3::BRI1-CITRINE transgene were complicated by the observation that it triggered an apparently dosage-dependent shoot phenotype (Fig. 5D). However, when MAKR5::BRI1-CITRINE or CVP2::BRI1-CITRINE transgenes were introduced into brl1 brl3 double mutant plants that were heterozygous for the bri1 mutation, surprisingly only two distinct phenotypic classes were observed in the next generation, in multiple, independent transgenic lines (Fig. S2F). The first class comprised short-rooted individuals corresponding to bri1 brl1 brl3 triple mutants (Fig. 5E), and the second class comprised individuals with fully developed roots with a wild-type appearance that were also similar to brl1 brl3 double mutant roots (Fig. 5F). The two phenotypic classes appeared in ratios of approximately 1:15, suggesting that both constructs fully (E,F) Root meristems of bri1 brl1 brl3 triple mutants (E) and brl1 brl3 double mutants (F) segregating from a bri1 +/− brl1 −/− brl3 −/− mother plant that carried a hemizygous CVP2::BRI1-CITRINE transgene. Note the autofluorescence background in A-C as exemplified by (non-transgenic) F (these images were taken with the same confocal settings). (G) Phenotypic analysis of the segregating progeny of a bri1 +/− brl1 −/− brl3 −/− mother plant that carried a hemizygous single insertion of a CVP2::BRI1-CITRINE transgene. Overall root growth, meristem appearance and protophloem gaps were scored for an entire sample of 347 seedlings, indicating rescue of the bri1 brl1 brl3 short root phenotype (G,G′) as well as the protophloem phenotype (G″,G‴) by the transgene. Asterisks indicate protophloem sieve element strands. Arrow in E indicates gap cells. Colors in G″ and G‴ refer to the subgroups indicated in G′. rescued the short root phenotype of the triple mutant. More detailed analyses confirmed this notion. For instance, in the progeny of a bri1 −/+ brl1 −/− brl3 −/− plant that carried a single CVP2::BRI1-CITRINE T-DNA insertion, 11 short root individuals and 165 long root individuals were observed. This segregation ratio was significantly different from 1/4 (P<0.0001; Fisher's exact test), but not from 1/16. Genotyping confirmed that all the 11 short root individuals were bri1 brl1 brl3 triple mutants and did not carry the transgene, consistent with the absence of CITRINE signal in their roots. Moreover, the meristems of seedlings with long roots that expressed the CVP2::BRI1-CITRINE transgene were essentially indistinguishable from the meristems of seedlings with long roots that were not transgenic (i.e. brl1 brl3 double mutants, or such double mutants heterozygous for bri1). This included 25% of complemented triple mutants, as we did not observe aborted seeds or arrested, non-germinating seedlings. In summary, in all individuals that expressed the transgenes, we observed long roots with a growth rate equal to wild type, restored cell elongation, and absence of protophloem gaps. To confirm this finding with higher statistical confidence, we analyzed an entire progeny of 347 seedlings segregating from another, independent bri1 −/+ brl1 −/− brl3 −/− plant that carried a single hemizygous CVP2::BRI1-CITRINE T-DNA insertion in confocal microscopy. Again, not a single protophloem gap was observed in meristems with wild-type appearance (Fig. 5G), demonstrating that the transgene complemented not only the reduced root growth and altered meristem morphology of bri1 brl1 brl3 triple mutants, but also their protophloem differentiation defect. However, we could not conclusively determine whether the transgenes also rescued the increase in formative cell divisions, because the presence of the transgenic BRI1 copy prevented conclusive genotyping of the (intron-less) endogenous locus except to confirm presence or absence of the mutant allele. We also could not recover triple mutants that were homozygous for the transgenes, because the shoot phenotype was not fully rescued. Therefore, with the possible exception of restricting formative cell divisions, we have shown that protophloem-specific BR signaling is sufficient to mediate the major effects of systemic BR action throughout the root meristem.

DISCUSSION
In this study, we evaluated whether BR signaling has a role in root protophloem development as recently proposed by the finding that the OPS protein can interact with, and thereby inhibit, BIN2 (Anne et al., 2015). We did indeed observe a substantial contribution of BR signaling to the progression of sieve element differentiation in BR receptor triple mutants, which appeared to be moderate, however, when compared with brx or ops mutants. Our findings suggest that OPS action cannot be entirely mediated through BIN2 and its downstream targets. In part, our results are incongruent with previous studies. For example, whereas it has been reported that bikinin treatment, but not BR treatment, could rescue ops gap cells (Anne et al., 2015), we observed the inverse. Such divergent outcomes could be explained by different assay or sampling conditions, or by small sample sizes that can lead to statistically significant yet not reproducible outcomes in individual replicate experiments. Unfortunately, this limitation is sometimes difficult to overcome because of the technical challenge and the major time investment associated with microscopic investigation of the protophloem. However, the notion emerging from our results that any OPS effects on BIN2 activity might not be mediated by the canonical downstream factors of BR signaling is also supported by the relatively small, statistically borderline significant effect of dominant bzr1-D or bes1-D mutation on ops root phenotypes (Anne et al., 2015), especially when compared with the substantial bzr1-D or bes1-D rescue of bri1 root growth (Chaiwanon and Wang, 2015). The finding that bzr1-D or bes1-D mutations have no effect on ops cotyledon vasculature phenotypes, and that OPS gain-of-function suppresses bri1 and bin2 dwarfism (Anne et al., 2015) supports this idea. Moreover, given that the link between BIN2 and OPS is so far exclusively based on observations in OPS gain-of-function scenarios, it still appears possible that BIN2 inhibition by OPS might represent a neomorphic or context-specific scenario. Indeed, our finding that BR signaling contributes to protophloem differentiation, but that GSK3 or BIN2 inhibition is not sufficient to substitute for loss of BR signaling in this context, supports this notion and is also consistent with normal protophloem development in bin2 gain-of-function mutants (Anne et al., 2015). Therefore, it still appears possible that BR perception contributes to protophloem differentiation through a GSK3/BIN2-independent route, or through one of the few GSK3 proteins that do not respond to bikinin. Alternatively, the impact of BR signaling on sieve element differentiation might be indirect. For instance, it is conceivable that the loss of BR receptors might increase the accessible pool of their co-receptors, which could then result in stoichiometric shifts towards other receptor signaling systems, such as the CLE45-BAM3 pathway, which suppresses protophloem differentiation (Depuydt et al., 2013). Our observation that the root phenotypes of brx and bri1 mutants are additive would be consistent with a parallel, possibly indirect, action of BR signaling in protophloem development.
Our observations suggest that in the root BR signaling is primarily required for proper cell elongation and to restrict the extent of formative divisions. Our finding that bri1 brl1 brl3 triple mutants produce more root cell files also matches with a previously reported increase of cyclin B1 activity in the tip of bri1 mutant root meristems (Gonzalez-Garcia et al., 2011), which already display this phenotype. Our data also indicate that meristem size in bri1 brl1 brl3 triple mutants can be entirely explained by reduced cell elongation. This might seem surprising given previous reports of altered meristematic cell number in brassinosteroid mutants (Gonzalez-Garcia et al., 2011;Hacham et al., 2011;Chaiwanon and Wang, 2015;Vragovićet al., 2015). However, the reported changes were typically rather small and not always congruent between studies. This could reflect different assay conditions and/or genetic materials, but also different quantification strategies. For instance, some authors quantified meristematic cortex cells (Hacham et al., 2011;Chaiwanon and Wang, 2015), whereas others scored epidermal cell files (Gonzalez-Garcia et al., 2011). The latter approach could be misleading if quantification is not consistently done in trichoblast or atrichoblast cell files only, because cell division and elongation rate are different in the two. In our study, we found that reduced cell elongation can by itself explain the small meristem size of bri1 brl1 brl3 triple mutants. However, it apparently cannot explain the relatively higher overall root growth reduction (i.e. ∼50% reduction in root length compared with ∼25% reduction in mature cell length). One explanation for this growth reduction would be an influence of BR signaling on the balance between cell proliferation and differentiation across cell layers, as previously proposed (Vragovićet al., 2015). Alternatively, bri1 brl1 brl3 roots might simply be growing slowly, because the supernumerary formative cell divisions have to be accommodated in space and time. In support of the latter idea, a similar short root phenotype can be created artificially by transgenic induction of formative cell divisions (De Rybel et al., 2013).
Perhaps our most surprising result was the observation that BRI1 activity that was restricted to the developing protophloem could rescue all major aspects of the bri1 brl1 brl3 triple mutant root phenotype. Such non-cell-autonomous effects of BR signaling have been described extensively before (Fridman et al., 2014;Hacham et al., 2011;Savaldi-Goldstein et al., 2007;Vragovićet al., 2015); however, their penetrance was typically found to be partial. Another example of non-cell-autonomous BR effects is our finding that increased stele cell divisions were suppressed in GL2::BRI1-GFP triple plants, whereas epidermal cell divisions were increased in SHR::BRI1-GFP triple plants, and combination of the two transgenes led to their mutual compensation and thus an intermediate phenotype. This points to flexible non-cell-autonomous accommodations of local growth effects, which could result from physiological consequences of sitespecific BR signaling action as proposed by transcriptomic analyses (Chaiwanon and Wang, 2015;Vragovićet al., 2015). Alternatively, they could reflect reactions to mechanical stress. Irrespective of the mechanism, protophloem-specific BR perception appears to be sufficient to systemically convey the major effects of BR action in the root meristem context, once more underlining the essential function of the protophloem for overall root system growth and organization (Depuydt et al., 2013;Rodriguez-Villalon et al., 2015).

Plant materials and phenotyping
The mutants and transgenic lines used in this study were all in Col-0 background and most have been described before (Geldner et al., 2007;Rodriguez-Villalon et al., 2014;Scacchi et al., 2010;Vragovićet al., 2015). Plant tissue culture, (confocal) microscopy, histology and physiological assays were performed according to standard procedures (Rodriguez-Villalon et al., 2014. Brassinolide was purchased from Sigma, bikinin from Wako Chemicals.