AtHB7/12 Regulate Root Growth in Response to Aluminum Stress

Aluminum (Al) stress is a major limiting factor for plant growth and crop production in acid soils. At present, only a few transcription factors involved in the regulation of Al resistance have been characterized. Here, we used reversed genetic approach through phenotype analysis of overexpressors and mutants to demonstrate that AtHB7 and AtHB12, two HD-Zip I transcription factors, participate in Al resistance. In response to Al stress, AtHB7 and AtHB12 displayed different dynamic expression patterns. Although both AtHB7 and AtHB12 positively regulate root growth in the absence of Al stress, our results showed that AtHB7 antagonizes with AtHB12 to control root growth in response to Al stress. The athb7/12 double mutant displayed a wild-type phenotype under Al stress. Consistently, our physiological analysis showed that AtHB7 and AtHB12 oppositely regulate the capacity of cell wall to bind Al. Yeast two hybrid assays showed that AtHB7 and AtHB12 could form homo-dimers and hetero-dimers in vitro, suggesting the interaction between AtHB7 and AtHB12 in the regulation of root growth. The conclusion was that AtHB7 and AtHB12 oppositely regulate Al resistance by affecting Al accumulation in root cell wall.


Introduction
A major factor constraint crop production on acidic soils worldwide is aluminum (Al) toxicity [1]. Excess Al 3+ in acid soils is toxic to plants. It has been reported that the rapid inhibition of root elongation is the earliest and most dramatic symptom of Al toxicity, resulting in a reduced and damaged root system that limits mineral nutrient and water uptake, and eventually leads to reduced production [2][3][4][5].
The root transition zone (TZ), located between the apical meristem and basal elongation zone (EZ) in root apex, is the critical site of perception of Al toxicity, which has been reported in the model plant Arabidopsis (Arabidopsis thaliana), wheat (Triticum aestivum), maize (Zea mays), sorghum (Sorghum bicolor), and common bean (Phaseolus vulgaris) [6][7][8][9][10][11]. The root TZ has also been thought to be the active OsFRDL4, OsMGT1 and OsCDT3 [34]. ART2, a homology of ART1, plays a supplementary role in Al tolerance in rice. [35]. OsWRKY22 was involved in regulation of Al-activated citrate acid secretion by regulating expression of OsFRDL4 which encodes a citrate transporter in rice [41]. ABSCISIC ACID, STRESS AND RIPENING (ASR) transcription factor ASR5 collaboratively acts with ASR1 in Al-stress response [42,43].
HD-ZIP transcription factors which are unique in plants are characterized as containing a homeodomain (HD) for binding DNA and a leucine zipper (Zip) motif for protein dimerization [44]. HD-ZIP I proteins from different species were reported to participate in various abiotic adversity responses. In rice, OsHOX22 acted as a negative regulator of drought and salt tolerance through an ABA-dependent pathway [45]. MtHB1 from Medicago truncatula controlled lateral root emergence under abiotic stresses [46]. ZmHDZ10 positively regulated drought and salt stress dependent on ABA signal pathway [47]. In Arabidopsis, AtHB6 was reported to be involved in ABA response by interacting with the PP2C phosphatase ABI1 and mediated drought stress responses acting downstream of ABI1 and ABI2 [48,49]. AtHB13 regulated response to freezing temperature via inducing genes expression which involved in cell membrane stabilization and inhibition of ice growth [50]. In addition, AtHB13 positively regulates a NAC transcription factor gene JUB1 expression to confer Arabidopsis drought tolerance [51]. AtHB7 and AtHB12 negatively module ABA signaling pathway by activating PP2C genes and repressing ABA receptors PYL5 and PYL8 [52]. Olsson et al. showed that AtHB7 and AtHB12 are negative developmental regulators in response to drought [53]. Our previous data showed that AtHB7 and AtHB12 were highly up-regulated among the differential expression genes in Col-0 wild type treated with or without Al [7].
Since AtHB7 and AtHB12 as two members of HD-ZIP I subfamily play important roles in regulation of tolerance to abiotic stress and were significantly induced by Al, we speculated that these two transcription factors might relate to regulation of Al resistance. Here, we showed that AtHB7 and AtHB12 acted in a cooperative manner under normal condition but in an antagonistic manner to regulate primary root growth under Al treatment.

AtHB7 and AtHB12 Induced by Al Stress, Which Was Independent of STOP1
AtHB7 and AtHB12, the most closed members in HD-ZIP I subfamily, were up-regulated 13 folds and 206 folds respectively by Al treatment [7]. To characterize the potential roles of AtHB7 and AtHB12 in root growth response to Al stress, we further analyzed the expression patterns of AtHB7 and AtHB12 under Al stress through a time-course analysis. Real-time qPCR results showed that AtHB7 mRNA level was dominant in early response stage and increased with Al treatment, whereas AtHB12 expression level appeared to be down-regulated initially, and then up-regulated at 6 h after Al exposure and surpassed AtHB7 accumulation level at 12 h ( Figure 1A). The different dynamic expression of AtHB7 and AtHB12 in response to Al stimuli suggests that AtHB7 and AtHB12 might play different roles in response to Al stress. To clarify the molecular mechanism of transcriptional regulation of AtHB7 and AtHB12 in response to Al stress, we firstly examined if STOP1, which is a key transcription factor during Al stress response and regulates multiple Al resistance genes expression, was involved in this process. The results showed that the expression patterns of AtHB7 and AtHB12 were not changed in stop1 mutant with or without Al treatment, indicating that Al-regulated AtHB7 and AtHB12 expression is independent of STOP1 ( Figure 1B,C).
To address the expression patterns of AtHB7 and AtHB12 in the presence of Al stress, the effect of exposure to Al stress on the spatial expression of AtHB7 and AtHB12 was analyzed by monitoring the expression of the AtHB7pro:eGFP-GUS and AtHB12pro:eGFP-GUS transgenes. After a 6 h exposure to 25 µM Al, GFP signals were highly detected in the root apex TZ, while there is no GFP signal in the root tips under non-stressed conditions ( Figure 2). To address the expression patterns of AtHB7 and AtHB12 in the presence of Al stress, the effect of exposure to Al stress on the spatial expression of AtHB7 and AtHB12 was analyzed by monitoring the expression of the AtHB7pro:eGFP-GUS and AtHB12pro:eGFP-GUS transgenes. After a 6 h exposure to 25 μM Al, GFP signals were highly detected in the root apex TZ, while there is no GFP signal in the root tips under non-stressed conditions ( Figure 2).

AtHB7 and AtHB12 Promote Root Growth by Regulating Cell Number and Cell Length
To investigate whether Al responsive AtHB7 and AtHB12 are involved in root growth, we generated 35S:AtHB7/12 overexpression transgenic plants, the dominant-negative constructs 35S:AtHB7/12-SRDX (SUPERMAN repression domain, SRDX) lines and athb7 and athb12 mutants edited by CRISPR CAS9 system (Supplementary Figure S1). Under normal condition, the roots of 35S:AtHB7 and 35S:AtHB12 were longer than those of the wild type, while 35S:AtHB7-SRDX and 35S:AtHB12-SRDX showed short-root phenotypes, which were also confirmed by root length analysis in the athb7 and athb12 mutants (Figure 3). A longitudinal zonation pattern analysis showed that root apical meristem (RAM) of 35S:AtHB7 and 35S:AtHB12 seedlings were significantly longer than that of the WT control ( Figure 4A,B). Cell numbers in the meristem zone (MZ) and cell length in the elongation zone (EZ) of 35S:AtHB7 and 35S:AtHB12 seedlings were both increased compared to

AtHB7 and AtHB12 Promote Root Growth by Regulating Cell Number and Cell Length
To investigate whether Al responsive AtHB7 and AtHB12 are involved in root growth, we generated 35S:AtHB7/12 overexpression transgenic plants, the dominant-negative constructs 35S:AtHB7/12-SRDX (SUPERMAN repression domain, SRDX) lines and athb7 and athb12 mutants edited by CRISPR CAS9 system (Supplementary Figure S1). Under normal condition, the roots of 35S:AtHB7 and 35S:AtHB12 were longer than those of the wild type, while 35S:AtHB7-SRDX and 35S:AtHB12-SRDX showed short-root phenotypes, which were also confirmed by root length analysis in the athb7 and athb12 mutants (Figure 3). A longitudinal zonation pattern analysis showed that root apical meristem (RAM) of 35S:AtHB7 and 35S:AtHB12 seedlings were significantly longer than that of the WT control ( Figure 4A,B). Cell numbers in the meristem zone (MZ) and cell length in the elongation zone (EZ) of 35S:AtHB7 and 35S:AtHB12 seedlings were both increased compared to Figure 2. Both AtHB7 and AtHB12 were induced by Al toxicity in transition zone. Expression of AtHB7pro:GFP-GUS and AtHB12pro:GFP-GUS in root tips of seedlings exposed without or with 20 µM Al for 6 h were analyzed by observation using confocal laser scanning microscopy.

AtHB7 and AtHB12 Promote Root Growth by Regulating Cell Number and Cell Length
To investigate whether Al responsive AtHB7 and AtHB12 are involved in root growth, we generated 35S:AtHB7/12 overexpression transgenic plants, the dominant-negative constructs 35S:AtHB7/12-SRDX (SUPERMAN repression domain, SRDX) lines and athb7 and athb12 mutants edited by CRISPR CAS9 system (Supplementary Figure S1). Under normal condition, the roots of 35S:AtHB7 and 35S:AtHB12 were longer than those of the wild type, while 35S:AtHB7-SRDX and 35S:AtHB12-SRDX showed short-root phenotypes, which were also confirmed by root length analysis in the athb7 and athb12 mutants ( Figure 3). A longitudinal zonation pattern analysis showed that root apical meristem (RAM) of 35S:AtHB7 and 35S:AtHB12 seedlings were significantly longer than that of the WT control ( Figure 4A,B). Cell numbers in the meristem zone (MZ) and cell length in the elongation zone (EZ) of 35S:AtHB7 and 35S:AtHB12 seedlings were both increased compared to the WT control ( Figure 4C-F). However, the athb7 mutant and the 35S:AtHB12-SRDX transgene showed a reduced size in the root MZ, and cortical cells in the EZ were shorter than those in WT, and cell numbers in the MZ were also strongly reduced ( Figure 4C-F). These results suggest that AtHB7 and AtHB12 affect both cell proliferation and cell elongation.

5
the WT control ( Figure 4C-F). However, the athb7 mutant and the 35S:AtHB12-SRDX transgene showed a reduced size in the root MZ, and cortical cells in the EZ were shorter than those in WT, and cell numbers in the MZ were also strongly reduced ( Figure 4C-F). These results suggest that AtHB7 and AtHB12 affect both cell proliferation and cell elongation.  . Asterisks indicate significant differences between wild type and transgenic lines at *: p < 0.05 and **: p < 0.01 by Tukey's test.

AtHB7 and AtHB12 Antagonistically Regulate Root Elongation under Al Stress
We further examined the root growth under Al stress. The results showed that 35S:AtHB7 transgenic lines displayed a reduced root-growth inhibition in response to Al stress, while 35S:AtHB7-SRDX transgenic lines and athb7 mutants all showed an enhanced root-growth inhibition after Al treatment ( Figure 5A). However, there were no difference between wild type and the 35S:AtHB7, 35S:AtHB7-SRDX transgenic lines and athb7 mutants when plants were exposed to solutions containing different toxic ions including sodium (Na + ), copper (Cu 2+ ), lanthanum (La 3+ ) and cadmium (Cd 2+ ), or to various pH solutions (4.5-6.5), indicating that AtHB7 responds to Al stress specifically (Supplementary Figure S2). These results indicate that AtHB7 is a positive regulator of Al resistance. In contrast to AtHB7, 35S:AtHB12 transgenic lines showed an elevated root growth inhibition in response to Al stress, while 35S:AtHB12-SRDX and athb12 mutants showed significantly less inhibition than the WT control ( Figure 5B). Similar to AtHB7, no differences in root growth were observed when WT and 35S:AtHB12, 35S:AtHB12-SRDX transgenic lines and athb12 mutant plants were exposed to other metal ions stress and different pH solutions (Supplementary Figure S3). These results indicate that AtHB12 negatively regulates Al resistance.

AtHB7 and AtHB12 Oppositely Regulate Aluminum Deposition in Cell Wall
To investigate the underlying mechanisms of how AtHB7/12 regulate Al resistance, Al content in roots of athb7 mutant, athb12 mutant and wild type was detected. Compared with wild type, athb7 mutant accumulated much more Al in roots, but athb12 mutant held less Al in roots ( Figure 6A). In depth, hematoxylin staining technique was used to visualize Al accumulation in root cell walls of wild type, AtHB7 or AtHB12 transgenic plants. The staining intensity was stronger in athb7 root tip region and decreased in AtHB7 overexpression plants compared to wild type ( Figure 6B). However, stronger intensity in AtHB12 overexpression plants and weaker stain in the athb12 mutant were observed ( Figure 6C). Taken together, these results indicate that AtHB7 and AtHB12 differently regulated Al tolerance by affecting Al accumulation in root cell wall. . Asterisks indicate significant differences between wild type and transgenic seedlings at *: p < 0.05 and **: p < 0.01 by Tukey's test.

AtHB7 and AtHB12 Antagonistically Regulate Root Elongation under Al Stress
We further examined the root growth under Al stress. The results showed that 35S:AtHB7 transgenic lines displayed a reduced root-growth inhibition in response to Al stress, while 35S:AtHB7-SRDX transgenic lines and athb7 mutants all showed an enhanced root-growth inhibition after Al treatment ( Figure 5A). However, there were no difference between wild type and the 35S:AtHB7, 35S:AtHB7-SRDX transgenic lines and athb7 mutants when plants were exposed to solutions containing different toxic ions including sodium (Na + ), copper (Cu 2+ ), lanthanum (La 3+ ) and cadmium (Cd 2+ ), or to various pH solutions (4.5-6.5), indicating that AtHB7 responds to Al stress . Asterisks indicate significant differences between wild type and transgenic seedlings at *: p < 0.05 and **: p < 0.01 by Tukey's test.

7
Al resistance. In contrast to AtHB7, 35S:AtHB12 transgenic lines showed an elevated root growth inhibition in response to Al stress, while 35S:AtHB12-SRDX and athb12 mutants showed significantly less inhibition than the WT control ( Figure 5B). Similar to AtHB7, no differences in root growth were observed when WT and 35S:AtHB12, 35S:AtHB12-SRDX transgenic lines and athb12 mutant plants were exposed to other metal ions stress and different pH solutions (Supplementary Figure S3). These results indicate that AtHB12 negatively regulates Al resistance. . Asterisks indicate significant differences between wild type and transgenic seedlings at *: p < 0.05 by Tukey's test.

AtHB7 and AtHB12 Oppositely Regulate Aluminum Deposition in Cell Wall
To investigate the underlying mechanisms of how AtHB7/12 regulate Al resistance, Al content in roots of athb7 mutant, athb12 mutant and wild type was detected. Compared with wild type, athb7 mutant accumulated much more Al in roots, but athb12 mutant held less Al in roots ( Figure 6A). In depth, hematoxylin staining technique was used to visualize Al accumulation in root cell walls of wild type, AtHB7 or AtHB12 transgenic plants. The staining intensity was stronger in athb7 root tip region and decreased in AtHB7 overexpression plants compared to wild type ( Figure 6B). However, stronger intensity in AtHB12 overexpression plants and weaker stain in the athb12 mutant were observed ( Figure 6C). Taken together, these results indicate that AtHB7 and AtHB12 differently regulated Al tolerance by affecting Al accumulation in root cell wall.

AtHB7 and AtHB12 might Act in Homodimer or Heterodimer to Control Root Growth under Al Stress
HD ZIP transcription factors are known for containing a homeodomain / leucine zipper domain, the former domain specifically binding to DNA and the latter domain mediating the formation of protein dimers. Based on the high level homology between AtHB7 and AtHB12 and the distinctively opposite function in Al-stress response, we speculated that AtHB7 and AtHB12 may interact to form a heterodimer that affects protein activity or stability to antagonistically regulate Al resistance in Arabidopsis. Yeast two-hybrid assays showed that AtHB7 could not only form homodimer with itself but also form heterodimer with AtHB12 ( Figure 7A). Further, the athb7/hb12 double mutant exhibited similar Al resistance with wild type ( Figure 7B). These results suggest that AtHB7 and AtHB12 may form a heterodimer to inhibit mutual activities in regulating downstream Al-response genes under Al stress.

AtHB7 and AtHB12 Might Act in Homodimer or Heterodimer to Control Root Growth under Al Stress
HD ZIP transcription factors are known for containing a homeodomain/leucine zipper domain, the former domain specifically binding to DNA and the latter domain mediating the formation of protein dimers. Based on the high level homology between AtHB7 and AtHB12 and the distinctively opposite function in Al-stress response, we speculated that AtHB7 and AtHB12 may interact to form a heterodimer that affects protein activity or stability to antagonistically regulate Al resistance in Arabidopsis. Yeast two-hybrid assays showed that AtHB7 could not only form homodimer with itself but also form heterodimer with AtHB12 ( Figure 7A). Further, the athb7/hb12 double mutant exhibited similar Al resistance with wild type ( Figure 7B). These results suggest that AtHB7 and AtHB12 may form a heterodimer to inhibit mutual activities in regulating downstream Al-response genes under Al stress. a heterodimer that affects protein activity or stability to antagonistically regulate Al resistance in Arabidopsis. Yeast two-hybrid assays showed that AtHB7 could not only form homodimer with itself but also form heterodimer with AtHB12 ( Figure 7A). Further, the athb7/hb12 double mutant exhibited similar Al resistance with wild type ( Figure 7B). These results suggest that AtHB7 and AtHB12 may form a heterodimer to inhibit mutual activities in regulating downstream Al-response genes under Al stress. . Asterisks indicate significant differences between wild type and mutants at *: p < 0.05 by Tukey's test.

Discussion
Al stress in one of the most limiting factors for plant growth and crop productivity in acidic soils. Identification of the molecules involved in Al tolerance mechanisms is a current challenge for Al toxicity. Previous studies identified STOP1 as a key transcription factor that regulates primary root . Asterisks indicate significant differences between wild type and mutants at *: p < 0.05 by Tukey's test.

Discussion
Al stress in one of the most limiting factors for plant growth and crop productivity in acidic soils. Identification of the molecules involved in Al tolerance mechanisms is a current challenge for Al toxicity. Previous studies identified STOP1 as a key transcription factor that regulates primary root growth in response to Al stress. In the present study, we identified two novel regulators of Al resistance. AtHB7 and AtHB12, two HD-Zip I transcription factors, oppositely regulate Al resistance in Arabidopsis. This is directly proved by the fact that AtHB7 overexpression lines and athb12 mutants are tolerant to Al and athb7 mutant, and AtHB12 overexpression lines are sensitive to Al ( Figure 5). However, AtHB7 and AtHB12 are not associated with tolerance to other metal such as La 3+ , Cd 2+ , Cu 2+ , Na + and low pH ( Supplementary Figures S2 and S3). Therefore, we conclude that AtHB7 and AtHB12 are specifically involved in Al resistance. The HD-ZIP transcription factor family consists of 47 members in Arabidopsis and several HD-ZIP I subfamily members were found to take part in response to unfavorable conditions such as salt stress, drought stress, ABA and cold stress. In the study, we firstly demonstrated that these two proteins belonging to HD-ZIP I subfamily function in Al resistance mechanisms.
Although AtHB7 and AtHB12 are the most closely related members in HD-Zip I family (Supplementary Figure S4), AtHB7 and AtHB12 showed opposite function in Al resistance. Previous study also showed that AtHB7 was repressed by water loss, while AtHB12 was induced by water loss under water stress conditions; in addition, AtHB7 and AtHB12 play different roles in certain biological process or developmental stage like stomata aperture and senescence [54].
It is interesting to note that AtHB7 and AtHB12 play opposite roles in Al resistance, while both promote root growth during normal development (Figures 3 and 5). The difference of DNA-binding homeodomains between AtHB7 and AtHB12 might lead to distinct target genes in response to Al stress; AtHB7 and AtHB12 also perhaps interact with different partners in Al-response due to the distinct leucine zipper domains between AtHB7 and AtHB12. For instance, some transcription factors interact with transcriptional repressor to negatively regulate target genes. Human transcription factor LEF1 was showed to interact with HDAC1 (Histone Deacetylase1) and its function shift from activation to repression regulated by β-Catenin-Histone Deacetylase [55]. BZR1 acts as either a transcriptional activator or a repressor partly depend on the context of promoters of target genes and trans-factors. BZR1-PIF4 (phytochrome interacting factor 4) heterodimer binds to the E-box motif to activate target expression [56]; while BZR1 interacts with TPL (Groucho/TUP1-like corepressor TOPLESS) through its EAR motif to transcriptionally repress target genes [57]. AtHB7 and AtHB12 were revealed to function as transcription activators in Arabidopsis under normal conditions by transient expression assays [44]. However, recent report showed that AtHB12 interacted with TFIIB through its AHA (aromatic large hydrophobic acidic residues) transactivation motif which locates in carboxy termini, while AtHB7 interacted with both TBP and TFIIB to activate target genes [58].
Our results showed that AtHB7 and AtHB12 were able to form heterodimers ( Figure 7A). The leucine zipper domain of HD Zip provides proper structure to allow the formation of dimers [59,60]. For example, AtHB5 was reported to be able to form homodimer and heterodimers with other HD-Zip I members [61]. OsHox1 as a member of HD-Zip II family could interact with HAT which is an Arabidopsis HD-Zip I member, whereas the corresponding overexpression transgenic plants show opposite phenotype such as developmental rate [62]. Hetero-dimerization of AtHB7 and AtHB12 seems to inhibit activities of each other under Al stress, and the athb7/hb12 double mutant show similar Al resistance with wild type ( Figure 7B).
In this study, we further demonstrated that AtHB7 and AtHB12 differently mediated Al-induced root elongation, which was attributed to distinct capacity of cell wall to bind Al ( Figure 6). It is in line with previous findings that AtHB12 regulates expression of cell wall associated genes. AtHB12 overexpression promotes leaf cell expansion by increasing expression of cell wall-related genes such as EXPANSINA10 and DWARF4 [63]. Cell wall is composed of pectin, hemicellulose and cellulose, and pectin and hemicellulose were accepted as the major targets to bind Al [64][65][66]. It is well known that compositional and structural modification of the cell wall, especially pectin methylation and modification of xyloglucan oligosaccharides, exert crucial function in Al resistance. Functional genes associated with cell wall modification were widely characterized in Arabidopsis such as WAK1 (cell wall-associated receptor kinase 1), PME46, LUH (LEUNIG_HOMOLOG), STAR1 (Sensitive to Al rhizotoxicity1), ALS3 (ALUMINUM SENSITIVE 3), and XTH17 and XTH31 [19,20,[22][23][24][25]67,68]. Next, we intend to explore whether the above genes mediated AtHB7/12 regulated Al resistance, and discover more novel Al-resistance genes controlled by AtHB7 or AtHB12 based on whole transcriptome sequencing to enrich our knowledge in mechanisms of plants Al resistance.

Plant Culture and Treatment
Arabidopsis cultured in Petri dishes with 1/2 MS solid medium or in nutrient solutions or in soils were placed in a plant growth chamber at 22 • C, 150 µmolm −2 s −1 light intensity with a long-day photoperiod (16-h light/8-h dark). For plate culture, seeds were surface disinfected with 10% sodium hypochlorite solution for 5 min, and then washed three times with sterile water [69]. The disinfected seeds were sowed onto 1/2 MS medium and incubated at 4 • C in darkness for 2 days. Then, the plates with seeds were vertically placed in the chamber. For Al and different metal treatment, Arabidopsis seeds soaked in sterile water at 4 • C overnight germinated in 2% MGRL solution (pH 5.0) without or with 6 µM AlCl 3 , 1 mM LaCl 3 , 3.5 mM CdCl 2 , 2 mM CuCl 2 , 10 mM NaCl for 6 days. For different pH treatment, the soaked Arabidopsis germinated in 2% MGRL solution with pH 6.5, 6.0, 5.5, 5.0, 4.7, 4.5 for 6 days. Treatment solution was renewed every day. Then, the primary root length was analysis by Image J (developed by National Institutes of Health) after treatment.

Materials and Transgenic Arabidopsis Construction
athb7 and athb12 mutants were obtained by YAO-promoter driven CRISPR/cas9 gene editing method according to previous reports [70]. Briefly, two sgRNA targets of each gene were designed, then cloned into middle vector in turn and eventually two pYAO:hSpCas9-target1-sgRNA-target2-sgRNA recombinant vectors for editing AtHB7 and AtHB12 respectively were constructed. The amplified AtHB7 and AtHB12 ORF were respectively integrated into pENTR ™ TOPO ® vector using pENTR™/D-TOPO™ Cloning Kit (Thermo Fisher, catalog number K2400-20), and subsequently transferred to pB7WG2 Gateway destination vector to obtain 35S:AtHB7 and 35S:AtHB12 overexpression vectors. AtHB7 and AtHB12 coding sequences were respectively fused with the EAR-repression domain (SRDX) under the control of the 35S promoter to generate 35S:AtHB7-SRDX and 35S:AtHB12-SRDX recombinant vectors. In total, 2.1 Kb DNA fragment upstream from initiation codon of AtHB7 was cloned to drive repoter genes expression using Gateway cloning technology to construct AtHB7pro:eGFP-GUS. Likewise, approximate 2.2 Kb promoter sequence of AtHB12 was amplified to obtain AtHB12pro:eGFP-GUS plasmid.

Phylogenetic Analysis
The amino acid sequences of HD-ZIP I members were retrieved from TAIR database and then were aligned with Clustal 2.0 software [71]. The phylogenetic tree was constructed using MEGA 7.0 by neighbor-joining (NJ) method with 1000 bootstrap replications [72].

RNA Isolation and qRT PCR
Total RNA was extracted from intact roots of 6-day-old plants treated without or with 20 µM Al using TaKaRa MiniBEST Plant RNA Extraction Kit (TAKARA, catalog number 9769). At first, genome DNA was removed by gDNase (TIANGEN BIOTECH Co., Ltd., catalog number KR116) (TIANGEN Biotech, Beijing, China). Then, 1 µg RNA was synthesized into the first cDNA strand with FastKing RT Kit according to manufacturers' protocols (TIANGEN BIOTECH Co., Ltd., catalog number KR116) (TIANGEN Biotech, Beijing, China). The real-time quantitative RT-PCR was performed with SuperReal PreMix Plus (SYBR Green) provided by TIANGEN company using CFX Connect Real-Time System (Bio-Rad, Hercules, CA, USA) [73]. qPCR reaction conditions were 45 cycles at 95 • C for 10 s, 56 • C for 20 s, and 72 • C for 20 s. Expression levels were normalized to the expression level of UBQ10. All qRT PCR experiments were done three repeats from different biological samples.

Confocal Microscopy
6-day-old AtHB7pro:eGFP-GUS or AtHB12pro:eGFP-GUS seedlings were treated in 2% MGRL solution without or with 20 µM Al at pH 5.0. Then, the roots were immersed in 1 µg/ml propidium iodide (PI) solution and imaged using a confocal laser scanning microscope (Zeiss LSM700). For PI and GFP channels, the 543 nm and the 488 nm wavelengths were used for excitation, respectively. Four-day-old seedlings of wild type, AtHB7 or AtHB12 transgenic plants were stain with PI for growth analysis. Cortical cell numbers in the root apical meristem (RAM) and length of elongation zone cell were determined as previous described [74]. In detail, RAM comprised of the organized cells from QC to the first elongated cell. Elongation zone boundaries were from the first dramatically elongating cell to the cell developing root hairs.

Al Content Determination
Seeds of wild type, athb7 and athb12 mutants germinated in 2% MGRL solution at pH 5.5 for 6 days and then the seedlings were exposed to 2% MGRL solution with 10 µM Al at pH 5.0. After 24 h, the roots were excised, washed three times with deionized water, dried and digested with 65% HNO 3 . Al content in roots was measured by GF-AAS (SHIMADZU, Japan).

Hematoxylin Staining
Arabidopsis seeds were germinated in 2% MGRL solution at pH 5.5 for 6 days and then the seedlings were transferred to 2% MGRL solution supplemented with 20 µM Al at pH 5.0 for 6 h. After treatment, the seedlings were stained with 0.2% hematoxylin containing 0.02% KIO 4 for 20 min and washed three times with deionized water before observations. The samples were observed and photographed using an Olympus BX53 microscope equipped with an Olympus DP72 camera system.

Yeast Two-Hybrid Assay
Analysis of protein-protein interactions in yeast were carried out using the Matchmaker Gold Yeast Two-Hybrid System provided by TAKARA. The coding sequence of AtHB7 and AtHB12 were transferred from pENTR TM to the pGBKT7 (BD) vector by Gateway LR reaction to generate the constructs AtHB7-BD and AtHB12-BD as the preys. Similarly, the amplified coding sequence of AtHB7 and AtHB12 were fused with the PGADT7 (AD), resulting in the baits. The prey and bait were co-transformed into yeast cells. The transformed yeasts were selected on dropout medium without Leu and Trp. Then analysis of interaction was indicated by the growth of yeast colonies on synthetic medium lacking Leu, Trp, His or lacking Leu, Trp, His, and Ade with 3-AT incubated at 30 • C for 3 days.

Statistical Analysis
Statistical analysis was performed by Tukey's test among treatments. Three biological replicates for each treatment were conducted for the statistical analysis in this article. Asterisks or different letters in the figures indicated significant differences as follows: * p < 0.05, ** p < 0.01, and *** p < 0.001.