Zinc-Dependent Global Transcriptional Control, Transcriptional Deregulation, and Higher Gene Copy Number for Genes in Metal Homeostasis of the Hyperaccumulator Arabidopsis halleri 1[W] Ina N. Talke 2 , Marc Hanikenne 2 , and Ute Kra ¨mer* Max Planck Institute of Molecular Plant Physiology, D–14476 Potsdam-Golm, Germany The metal hyperaccumulator Arabidopsis halleri exhibits naturally selected zinc (Zn) and cadmium (Cd) hypertolerance and accumulates extraordinarily high Zn concentrations in its leaves. With these extreme physiological traits, A. halleri phylogenet- ically belongs to the sister clade of Arabidopsis thaliana. Using a combination of genome-wide cross species microarray analysis and real-time reverse transcription-PCR, a set of candidate genes is identified for Zn hyperaccumulation, Zn and Cd hypertolerance, and the adjustment of micronutrient homeostasis in A. halleri. Eighteen putative metal homeostasis genes are newly identified to be more highly expressed in A. halleri than in A. thaliana, and 11 previously identified candidate genes are confirmed. The encoded proteins include HMA4, known to contribute to root-shoot transport of Zn in A. thaliana. Expression of either AtHMA4 or AhHMA4 confers cellular Zn and Cd tolerance to yeast (Saccharomyces cerevisiae). Among further newly implicated proteins are IRT3 and ZIP10, which have been proposed to contribute to cytoplasmic Zn influx, and FRD3 required for iron partitioning in A. thaliana. In A. halleri, the presence of more than a single genomic copy is a hallmark of several highly expressed candidate genes with possible roles in metal hyperaccumulation and metal hypertolerance. Both A. halleri and A. thaliana exert tight regulatory control over Zn homeostasis at the transcript level. Zn hyperaccumulation in A. halleri involves enhanced partitioning of Zn from roots into shoots. The transcriptional regulation of marker genes suggests that in the steady state, A. halleri roots, but not the shoots, act as physiologically Zn deficient under conditions of moderate Zn supply. Metal hyperaccumulation is a rare trait found in approximately 440 plant taxa, 11 of which have been reported to hyperaccumulate zinc (Zn; Reeves and Baker, 2000). Individuals of these taxa accumulate very high Zn concentrations of more than 10,000 mgg 21 Zn in leaf dry biomass (Baker and Brooks, 1989). The accumu- lation of higher metal concentrations in shoots than in roots has been described as a characteristic of hyper- accumulation in both hydroponic and soil culture (Baker et al., 1994; Kra ¨mer et al., 1996; Weber et al., 2004). All known metal hyperaccumulator taxa occur primarily on geologically metal-rich or metal-contaminated soils and possess extraordinarily high metal tolerance. The num- ber of Zn and nickel hyperaccumulator taxa is higher in the Brassicaceae family than in any other plant family, and includes a large number of nickel hyperaccumulator species in the section Odontarrhena of the genus Alys- sum and the two most commonly studied Zn hyper- accumulator model species Thlaspi caerulescens and Arabidopsis halleri (Reeves and Baker, 2000). Both subspecies of A. halleri, subsp. gemmifera and subsp. halleri, are hyperaccumulators of Zn, and subsp. halleri is able to grow healthily and to accumu- late very high leaf Zn concentrations of between 3,000 and 22,000 mgg 21 dry biomass within a wide dynamic range of external Zn concentrations in the field (Bert et al., 2000; Dahmani-Muller et al., 2000; Kubota and Takenaka, 2003). The leaves of individuals of some populations have additionally been reported to con- tain hyperaccumulator levels of cadmium (Cd) of more than 100 mgg 21 dry biomass (Dahmani-Muller et al., 2000). Together, A. halleri and its nonaccumulat- ing close relative Arabidopsis lyrata form the sister clade of Arabidopsis thaliana (Koch et al., 2000; Yogeeswaran et al., 2005). A. halleri and A. thaliana share on aver- age approximately 94% nucleotide sequence identity within coding regions (Becher et al., 2004). Leaves of both A. lyrata and A. thaliana commonly accumulate between 40 and 80 mgg 21 Zn, and gradually develop chlorosis when they contain significantly higher metal concentrations (Macnair et al., 1999; Desbrosses- Fonrouge et al., 2005). In hydroponic culture, the growth of A. halleri roots has been reported to tolerate at least 30-fold higher Zn and 10-fold higher Cd concentrations than root growth of A. lyrata (Macnair et al., 1999; Bert et al., 2003). These observations 1 This work was supported by the Deutsche Forschungsgemein- schaft (Kr 1967/3), the European Union (Research Training Network METALHOME, HPRN–CT–2002–00243), the German Federal Minis- try of Education and Research (Biofuture 0311877), and the Max Planck Institute of Molecular Plant Physiology. 2 These authors contributed equally to the paper. * Corresponding author; e-mail [email protected]; fax 49–331–5678–408. The author 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.plantphysiol.org) is: Ute Kra ¨mer ([email protected]). [W] The online version of this article contains Web-only data. www.plantphysiol.org/cgi/doi/10.1104/pp.105.076232 148 Plant Physiology, September 2006, Vol. 142, pp. 148–167, www.plantphysiol.org Ó 2006 American Society of Plant Biologists https://plantphysiol.org Downloaded on May 29, 2021. - Published by Copyright (c) 2020 American Society of Plant Biologists. All rights reserved.
20
Embed
Zinc-Dependent Global Transcriptional Control, Transcriptional Deregulation, and ... · Zinc-Dependent Global Transcriptional Control, Transcriptional Deregulation, and Higher Gene
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Zinc-Dependent Global Transcriptional Control,Transcriptional Deregulation, and Higher Gene CopyNumber for Genes in Metal Homeostasis of theHyperaccumulator Arabidopsis halleri1[W]
Ina N. Talke2, Marc Hanikenne2, and Ute Kramer*
Max Planck Institute of Molecular Plant Physiology, D–14476 Potsdam-Golm, Germany
The metal hyperaccumulator Arabidopsis halleri exhibits naturally selected zinc (Zn) and cadmium (Cd) hypertolerance andaccumulates extraordinarily high Zn concentrations in its leaves. With these extreme physiological traits, A. halleri phylogenet-ically belongs to the sister clade ofArabidopsis thaliana. Using a combination of genome-wide cross species microarray analysis andreal-time reverse transcription-PCR, a set of candidate genes is identified for Zn hyperaccumulation, Zn and Cd hypertolerance,and the adjustment of micronutrient homeostasis in A. halleri. Eighteen putative metal homeostasis genes are newly identified tobe more highly expressed inA. halleri than inA. thaliana, and 11 previously identified candidate genes are confirmed. The encodedproteins include HMA4, known to contribute to root-shoot transport of Zn in A. thaliana. Expression of either AtHMA4 orAhHMA4 confers cellular Zn and Cd tolerance to yeast (Saccharomyces cerevisiae). Among further newly implicated proteins areIRT3 and ZIP10, which have been proposed to contribute to cytoplasmic Zn influx, and FRD3 required for iron partitioning inA. thaliana. InA. halleri, the presence of more than a single genomic copy is a hallmark of several highly expressed candidate geneswith possible roles in metal hyperaccumulation and metal hypertolerance. Both A. halleri and A. thaliana exert tight regulatorycontrol over Zn homeostasis at the transcript level. Zn hyperaccumulation in A. halleri involves enhanced partitioning of Zn fromroots into shoots. The transcriptional regulation of marker genes suggests that in the steady state,A. halleri roots, but not the shoots,act as physiologically Zn deficient under conditions of moderate Zn supply.
Metal hyperaccumulation is a rare trait found inapproximately 440 plant taxa, 11 of which have beenreported to hyperaccumulate zinc (Zn; Reeves andBaker, 2000). Individuals of these taxa accumulate veryhigh Zn concentrations of more than 10,000mg g21 Zn inleaf dry biomass (Baker and Brooks, 1989). The accumu-lation of higher metal concentrations in shoots than inroots has been described as a characteristic of hyper-accumulation in both hydroponic and soil culture (Bakeret al., 1994; Kramer et al., 1996; Weber et al., 2004). Allknown metal hyperaccumulator taxa occur primarily ongeologically metal-rich or metal-contaminated soils andpossess extraordinarily high metal tolerance. The num-ber of Zn and nickel hyperaccumulator taxa is higher inthe Brassicaceae family than in any other plant family,and includes a large number of nickel hyperaccumulator
species in the section Odontarrhena of the genus Alys-sum and the two most commonly studied Zn hyper-accumulator model species Thlaspi caerulescens andArabidopsis halleri (Reeves and Baker, 2000).
Both subspecies of A. halleri, subsp. gemmifera andsubsp. halleri, are hyperaccumulators of Zn, andsubsp. halleri is able to grow healthily and to accumu-late very high leaf Zn concentrations of between 3,000and 22,000 mg g21 dry biomass within a wide dynamicrange of external Zn concentrations in the field (Bertet al., 2000; Dahmani-Muller et al., 2000; Kubota andTakenaka, 2003). The leaves of individuals of somepopulations have additionally been reported to con-tain hyperaccumulator levels of cadmium (Cd) ofmore than 100 mg g21 dry biomass (Dahmani-Mulleret al., 2000). Together, A. halleri and its nonaccumulat-ing close relative Arabidopsis lyrata form the sister cladeof Arabidopsis thaliana (Koch et al., 2000; Yogeeswaranet al., 2005). A. halleri and A. thaliana share on aver-age approximately 94% nucleotide sequence identitywithin coding regions (Becher et al., 2004). Leaves ofboth A. lyrata and A. thaliana commonly accumulatebetween 40 and 80 mg g21 Zn, and gradually developchlorosis when they contain significantly higher metalconcentrations (Macnair et al., 1999; Desbrosses-Fonrouge et al., 2005). In hydroponic culture, thegrowth of A. halleri roots has been reported to tolerateat least 30-fold higher Zn and 10-fold higher Cdconcentrations than root growth of A. lyrata (Macnairet al., 1999; Bert et al., 2003). These observations
1 This work was supported by the Deutsche Forschungsgemein-schaft (Kr 1967/3), the European Union (Research Training NetworkMETALHOME, HPRN–CT–2002–00243), the German Federal Minis-try of Education and Research (Biofuture 0311877), and the MaxPlanck Institute of Molecular Plant Physiology.
2 These authors contributed equally to the paper.* Corresponding author; e-mail [email protected];
fax 49–331–5678–408.The author responsible for distribution of materials integral to the
findings presented in this article in accordance with the policydescribed in the Instructions for Authors (www.plantphysiol.org) is:Ute Kramer ([email protected]).
[W] The online version of this article contains Web-only data.www.plantphysiol.org/cgi/doi/10.1104/pp.105.076232
148 Plant Physiology, September 2006, Vol. 142, pp. 148–167, www.plantphysiol.org � 2006 American Society of Plant Biologists
https://plantphysiol.orgDownloaded on May 29, 2021. - Published by Copyright (c) 2020 American Society of Plant Biologists. All rights reserved.
illustrate the dramatic alterations in metal homeostasisof A. halleri compared to both A. thaliana and A. lyrata,despite a very close phylogenetic relationship and highoverall sequence conservation among these three spe-cies (Clemens et al., 2002; Kramer and Clemens, 2005).
Researchers have begun to address the major dif-ferences between the metal homeostasis networks ofclosely related metal hyperaccumulator and nonaccu-mulator plants at the molecular level (Salt and Kramer,2000; Clemens et al., 2002; Kramer and Clemens, 2005).Using single gene and transcriptomic approaches, themost notable differences observed in hyperaccumula-tor compared to closely related nonaccumulator plantswere high relative transcript levels (RTLs) of severalcandidate genes (Pence et al., 2000; Assuncxao et al.,2001; Becher et al., 2004; Drager et al., 2004; Papoyanand Kochian, 2004; Weber et al., 2004; Hammond et al.,2006; Mirouze et al., 2006). Differences in the functionsof hyperaccumulator candidate gene products, whencompared to their A. thaliana orthologs, have not ex-tensively been documented so far (Bernard et al., 2004;Roosens et al., 2004).
The observed high transcript levels of candidategenes in A. halleri and in other hyperaccumulatormodel plants are based on a modified regulation oftranscript levels (Pence et al., 2000; Becher et al., 2004;Weber et al., 2004). This could involve at least fourpossible mechanisms. First, it is conceivable that forsome A. halleri genes, the presence or activity of prox-imal cis-elements or trans-factors controlling transcriptlevels could be altered in comparison to A. thaliana.Second, it is possible that the affinity of a Zn sensingmachinery that down-regulates transcript levels ofgenes encoding Zn acquisitory proteins may be re-duced in A. halleri. Third, Zn requirement of A. hallerimay merely be higher than in closely related non-accumulator plants. Finally, the high expression of aspecific gene of A. halleri could be an indirect conse-quence of another alteration in a different element ofthe metal homeostasis network, which is complexand tightly regulated in plants. Before dissecting themolecular basis underlying the modified regulationof expression of individual candidate genes inA. halleri,it is therefore important to understand globally howthe regulation of the metal homeostasis network isaltered.
For A. halleri MTP1, which is involved in Zn hyper-tolerance, the amplification of gene copy number con-tributes to enhanced MTP1 transcript levels, whencompared to A. thaliana or A. lyrata (Drager et al.,2004). It has not been investigated whether genomiccopy number is also increased for other genes involvedin naturally selected metal hyperaccumulation orhypertolerance of A. halleri.
Previously, root and shoot transcript profiles wereestablished in separate studies using A. halleri andA. thaliana grown under different conditions withAffymetrix 8 K chips representing only a limitednumber of genes (Becher et al., 2004; Weber et al.,2004). To obtain a comprehensive account of directly
comparable expression profiles for both roots andshoots of the two species, we employed the next gen-eration of Affymetrix Arabidopsis GeneChips, theATH1 microarray, covering more than 22,500 genes.
Here we present genome-wide transcriptional pro-files of hydroponically cultivated A. halleri and A.thaliana plants maintained under control conditions orupon a short-term exposure to high Zn concentrations.We report a set of candidate genes for metal hyper-accumulation and tolerance, for the adjustment ofmicronutrient homeostasis, and for the protectionfrom secondary abiotic stress in A. halleri. For severalof the most highly expressed candidate genes, genomiccopy number is estimated to be higher in A. halleri thanin A. thaliana by DNA gel blot. The detailed analysis oftranscriptional regulation of candidate genes by real-time reverse transcription (RT)-PCR is combined withexisting and newly generated knowledge on geneproduct function to develop a model of the alterationsin the metal homeostasis network underlying metalhyperaccumulation in A. halleri.
RESULTS
Establishment of Experimental Conditions
To compare the effect of external Zn supply onZn accumulation and partitioning in A. thaliana andA. halleri, plants were cultivated hydroponically andsupplied with different subtoxic Zn concentrations for 3weeks (Fig. 1, A and B). InA. halleri, the ratio of shoot Znconcentrations to root Zn concentrations was aboveunity in plants cultivated in a Zn-deficient hydroponicsolution containing no added Zn for 3 weeks (Fig. 1C).Continuous supply of Zn concentrations of between1 and 5 mM led to an increase in the shoot-to-root ratio ofZn concentrations, with maximum ratios approximatelytripling the ratio determined under Zn-deficient condi-tions. By contrast, A. thaliana partitioned into the shootsonly a small proportion of the Zn concentration accu-mulated in the roots (Fig. 1, A and B). This was reflectedin shoot-to-root Zn concentration ratios substantiallybelow unity, which decreased further with increasing Znsupply (Fig. 1C). These results indicated that A. halleriaccumulates Zn predominantly in the shoot even at verylow Zn supply, whereas A. thaliana primarily immobi-lizes Zn in the roots irrespective of the Zn supply. WhenA. halleri was cultivated in media containing high Znconcentrations of 30 or 300mM Zn, shoot-to-root ratios ofZn concentrations were below unity and only about one-quarter of the ratios measured at 0.3 mM Zn. In additionto the dramatic interspecies differences in Zn partition-ing between roots and shoots, the data suggest a species-specific modulation of Zn partitioning dependent on theZn supply.
To select suitable conditions and time points fortranscript profiling, transcript levels of three metalhomeostasis genes, ZIP4, ZIP9, and NAS2, were de-termined in preliminary experiments (Fig. 2). Thesegenes are known to be Zn responsive at the transcript
Zinc Homeostasis in Arabidopsis halleri and A. thaliana
Plant Physiol. Vol. 142, 2006 149
https://plantphysiol.orgDownloaded on May 29, 2021. - Published by Copyright (c) 2020 American Society of Plant Biologists. All rights reserved.
level in A. thaliana (Grotz et al., 1998; Wintz et al., 2003)and they have previously been identified as candidategenes for Zn accumulation and tolerance in A. halleri(Becher et al., 2004; Weber et al., 2004).
Plants of two different Zn steady states were sub-jected to short-term increases in external Zn supply. Inone experiment, plants grown under control (1 mM) Znconditions were exposed to excess Zn concentrationsof 30 mM Zn for A. thaliana and 300 mM Zn for A. halleri,for 2 and 8 h. The choice of concentrations reflectsnaturally selected Zn tolerance in A. halleri, whereasA. thaliana possesses only basic Zn tolerance (Clemens,2001; Becher et al., 2004). (In preliminary experiments,the chosen concentrations were the maximum Znconcentrations that did not cause a reduction in bio-mass gain in A. thaliana and A. halleri, respectively,during 4 d of exposure [M. Becher, I.N. Talke, A.N.Chardonnens, and U. Kramer, unpublished data].) Inan additional experiment, Zn-deficient plants culti-vated in a medium lacking added Zn for 3 weeks wereresupplied with a Zn concentration of 5 mM (Fig. 2).
In roots, Zn-induced changes in transcript levelswere observed as early as 2 h after increasing the Znconcentrations in the hydroponic medium (Fig. 2, toprow: ZIP4, ZIP9). In shoots, transcriptional responseswere observed 2 or 8 h after the increase in Zn concen-trations in the hydroponic medium (Fig. 2, bottom row:ZIP4, ZIP9, and NAS2). Transcripts of ZIP9 encoding atransporter of the ZIP family and NAS2 encoding anicotianamine synthase were not detectable in shoots ofeither Arabidopsis species grown in the presence ofsufficient (control) Zn, but were strongly induced uponZn deficiency (Fig. 2, ZIP9, NAS2). This suggested thatneither of the two species experienced Zn deficiencyunder the chosen control conditions.
Comparison of A. halleri and A. thaliana TranscriptProfiles with ATH1 GeneChips
According to the results from the preliminary ex-periments, microarrays were hybridized with labeledcRNAs prepared from roots and shoots of A. halleriand A. thaliana grown at 1 mM Zn and after 2 and 8 h ofexposure to excess Zn for roots and shoots, respec-tively (see ‘‘Material and Methods’’). Under controlconditions, higher signal intensities were detected fora total of 628 probe sets in roots, and for 1,739 probesets in shoots of A. halleri compared to A. thaliana (asoutlined in ‘‘Materials and Methods;’’ for visualizationsee all colored, nongray datapoints in Fig. 3, A and B;complete gene lists with annotations in SupplementalTables I and II). We also identified genes that re-sponded at the transcript level to excess Zn in A. halleri(see ‘‘Materials and Methods;’’ Supplemental Tables IIIand IV). To group the identified genes according tofunctional classes, a genome-wide annotation list wasassembled; this list includes the class of metal homeo-stasis-related genes, which contained all functionallycharacterized A. thaliana genes encoding proteins in-volved in metal transport, chelation, binding, andtrafficking and all other members of the respectiveprotein families, and genes encoding members of pro-tein families known to participate in metal homeosta-sis in other organisms (Supplemental Table XI: class‘‘Metal ion homeostasis;’’ Fig. 3, see ‘‘Material andMethods’’). The intersection of the gene list of metalhomeostasis-related genes and of the genes that areeither more highly expressed (Supplemental Tables Iand II) or Zn regulated in A. halleri (SupplementalTables III and IV) contained a total of 17 genes in theroots and 19 genes in the shoots, with seven genes in
Figure 1. Long-term Zn accumulation in A. halleri and A. thaliana roots and shoots. Hydroponically grown 4.5-week-old plantswere supplied with different Zn concentrations in the culture medium for 3 weeks before harvest. The Zn concentrations in root(A) and shoot (B) tissues were determined by inductively coupled plasma optical emission spectroscopy. White circles,A. thaliana; black circles: A. halleri. C illustrates the ratio of Zn concentrations in shoot versus root tissues. Note that at 10 mM
ZnSO4, the Zn concentration in A. thaliana roots exceeds that of A. halleri roots by more than 30-fold. In preliminaryexperiments, exposure of hydroponically grown A. thaliana to 30 mM ZnSO4 for a period of 3 weeks led to the onset of mildtoxicity symptoms; therefore, this Zn condition was excluded. Data are from one experiment representative of two independentexperiments. Values represent mean 6 SE, calculated from n 5 3 culture vessels per Zn condition. For C, data were log2transformed. Data were plotted with Origin-Pro 7 SR4 (v 7.0552; OriginLab Corporation) and logarithmic fits to the data wereperformed for A and B, a cubic B-spline connection was used for the smooth curves in C.
Talke et al.
150 Plant Physiol. Vol. 142, 2006
https://plantphysiol.orgDownloaded on May 29, 2021. - Published by Copyright (c) 2020 American Society of Plant Biologists. All rights reserved.
common between both roots and shoots (Fig. 3; Table I).We considered these genes as an initial set of majorcandidate genes for an involvement in naturally selectedmetal tolerance and hyperaccumulation in A. halleri.
The metal homeostasis genes newly identified hereas candidate genes encode proteins encompassingboth cellular functions known to be of primary im-portance in metal homeostasis: membrane transport ofmetals and metal binding/chelation. The 12 new can-didate genes encoding proteins that are likely to beinvolved in the transport of Zn or other metals acrossmembranes were Zn-regulated transporter/iron-regulatedtransporter-like protein 10 (ZIP10; roots, 111-fold), metaltolerance protein 8 (MTP8; shoots, 65-fold), heavy metal-associated domain-containing protein 4 or heavy metalATPase 4 (HMA4; shoots, 53-fold), ferric reductase defec-tive 3 (FRD3; roots, 45-fold), phosphate transporter 1;4
(PHT1;4; roots, 19-fold), cation/H1 exchanger 18 (CHX18;roots, 6.5-fold), MTP11 (shoots, 4.7-fold), P-type ATPasefrom Arabidopsis 1 (PAA1/HMA6; shoots, 4.6-fold),iron-regulated transporter 2 (IREG2; roots, 4.5-fold),iron-regulated transporter 3 (IRT3; roots, 4.5-fold), yellow-stripe 1-like protein 6 (YSL6; shoots, 4.4-fold), and ZIP3(down-regulated 2.7-fold in roots following exposure toexcess Zn; Fristedt et al., 1999; Curie et al., 2001; Maseret al., 2001; Rogers and Guerinot, 2002; Blaudez et al.,2003; Jensen et al., 2003; Lahner et al., 2003; Shikanai et al.,2003; Green and Rogers, 2004; Hussain et al., 2004; Shinet al., 2004). A total of five new candidate genes encodingmetal-binding proteins or proteins involved in the bio-synthesis of metal chelators included ferritin 2 (FER2;shoots, 11-fold), S-adenosyl-Met synthetase 1 (SAMS1;shoots, 10-fold), SAMS2 (shoots, 8.5-fold), FER1 (shoots,5.1-fold), and SAMS3 (roots and shoots, around 5-fold;
Figure 2. Expression analysis of known Zn-regulated marker genes. Real-time RT-PCR analysis of expression of ZIP4, ZIP9, andNAS2 genes in roots (top section) and shoots (bottom section) of hydroponically grownA. halleri (black bars) andA. thaliana (lightgray bars). Twodifferent basal Zn conditions uponwhich additional Znwas suppliedwere assessed: Short Zn oversupply denotes ashort-term high Zn supply to plants precultivated at control Zn concentrations; forA. thaliana, 30mMZn, and forA. halleri, 300mM
Zn were added to the culture medium for 2 and 8 h before harvest. Zn deficiency/resupply denotes a supply of 5 mM Zn to Zn-deficient plants for 2, 8, and 24 h before harvest. Note that substantial changes in root transcript levels occur as soon as 2 h after thechange of Zn concentration in the culturemediumwhereas shoot transcript levels reflect the change in external Zn supply only at alater time point. Data shown are transcript levels relative to EF1a from one experiment representative of two independentbiological experiments. Tissues from at least three culture vessels containing three plants each were pooled for each condition.Geneswere concluded to be not expressed (n.e.) whenCT valueswere above a threshold of 35 and reaction efficiencies (REs) withthe respective cDNAwere more than 5% below the mean RE for the respective primer pair (see also ‘‘Material and Methods’’).
Zinc Homeostasis in Arabidopsis halleri and A. thaliana
Plant Physiol. Vol. 142, 2006 151
https://plantphysiol.orgDownloaded on May 29, 2021. - Published by Copyright (c) 2020 American Society of Plant Biologists. All rights reserved.
Peleman et al., 1989; Petit et al., 2001). The compoundS-adenosyl-Met is the substrate for nicotianamine syn-thase, and the A. halleri nicotianamine synthase genesNAS2 andNAS3have previously been implicated in Zntolerance of A. halleri (Becher et al., 2004; Weber et al.,2004; see also Table I).
Two additional, closely related genes encoding theprotein disulfide isomerases 1 and 2 were includedamong the set of candidate genes analyzed here be-cause of their interesting regulation and potentialfunction in metal detoxification (PDI1: shoots, 39-foldhigher in A. halleri; PDI2: shoots, 17-fold; Table I;Rensing et al., 1997; Narindrasorasak et al., 2003).These proteins are the closest sequence relatives ofa previously characterized endoplasmic reticulum-localized protein from castor bean (Ricinus communis),which catalyzes the isomerization of disulfide bridgeswithin proteins. Enhanced PDI activity is thought toprotect plants from stress that results in the modifica-tion of thiol groups of endoplasmic reticulum proteins(Houston et al., 2005).
Microarray Analysis of Zn-Dependent
Transcriptional Regulation
Among candidate genes, the microarray data sug-gested a trend toward decreasing transcript levels inresponse to excess Zn in roots of both A. thaliana and A.halleri (Table I). Conversely, in A. halleri roots transcriptlevels of ZIP and NAS genes were increased under Zndeficiency, when compared to control conditions (TableI). In shoots of A. thaliana, transcript levels of severalcandidate genes were increased in response to excessZn, namely for PDI and SAMS genes, and for YSL6(Table I). In A. halleri, transcript levels of these geneswere constitutively high and unchanged in response toexcess Zn (Table I). Our data suggested that, especiallyin shoots, A. halleri responds differently to excess Znthan A. thaliana at the transcript level.
We investigated whether the species-specific pat-terns observed in metal-dependent transcriptionalregulation of candidate genes are also reflected glob-ally across the entire panel of genes represented on theATH1 microarray. In roots of A. halleri and A. thaliana,exposure to excess Zn resulted in increased transcriptlevels for 0.5% and 0.9%, respectively, of all expressedgenes, and in decreased transcript levels for only0.09% and 0.07%, respectively, of all expressed genes(Supplemental Tables III and V, P , 0.1). Thus, thepredominant directions of transcriptional responsesto excess Zn were identical in roots of A. halleri andA. thaliana, with a global trend toward up-regulationopposing a trend toward down-regulation among ourset of candidate genes in both species. In shoots ofA. thaliana, among all expressed genes, transcriptional
Figure 3. The relation between ATH1 microarray signals of GeneChipshybridized with cRNA from A. halleri and from A. thaliana grownunder control Zn conditions. For each of the 22,810 probe sets on theATH1 array, the mean signal derived from the two replicate arrays ofroots (A) and shoots (B) of A. halleri was plotted on the y axis and therespective mean signal of A. thaliana was plotted on the x axis. Signalvalues represent mean signal intensities calculated from MAS 5.0-scaled raw signals of replicate arrays, without further normalization.Symbols and colors reflect the data filtering which was applied toidentify candidate genes more highly expressed in A. halleri comparedto A. thaliana and involved in metal ion homeostasis or relatedprocesses. Gray diamonds represent all probe sets that did not passany filter. Yellow squares represent probe sets that exhibited a 4-foldhigher signal intensity in A. halleri than in A. thaliana and passed thevalue filter (for details see ‘‘Materials and Methods’’). Blue trianglesrepresent probe sets that, in addition, passed a first annotation filter (I)of genes annotated as metal ion (Zn, Fe, Cu, and Mn) binding,membrane and transport associated, oxidative stress protection asso-ciated, pathogen response related, and RNA metabolism associated(Supplemental Table XI). Red circles represent probe sets that passed afurther annotation filter (II) that included only genes with a predictedfunction in metal ion homeostasis. Candidate genes that passed the
value filter and either annotation filter and that were chosen for furtheranalysis by real-time RT-PCR are represented by green triangles (passedannotation filter I) or circles (passed annotation filter II) with a blackoutline (for details see Table I).
Talke et al.
152 Plant Physiol. Vol. 142, 2006
https://plantphysiol.orgDownloaded on May 29, 2021. - Published by Copyright (c) 2020 American Society of Plant Biologists. All rights reserved.
Table I. Summary of microarray data for A. halleri root and shoot candidate genes
ATH1 probe set identifiers and the corresponding A. thaliana locus identifiers are given for candidate genes more highly expressed in A. halleri thanin A. thaliana under control conditions (genes that passed the value filter as described in ‘‘Materials and Methods’’). In addition, genes were includedthat showed an at least 2-fold up-regulation or down-regulation upon short-term Zn supply in A. halleri. Raw signal denotes mean signal valuecalculated from signal values exported from MAS 5.0 using GeneSpring. Microarray data for short-term Zn supply and Zn deficiency are given wherethe regulation compared to controls was at least 2-fold. One sample Student’s t tests were performed in GeneSpring on ratios from cross species and onfold changes from within-species comparisons, P values are indicated as follows: * 5 P , 0.05; ** 5 P , 0.01; *** 5 P , 0.001.
ATH1 Probe
Set Identifier
AGI Locus
Identifier
Gene
Name
A. halleri
Raw Signal
A. halleri/
A. thaliana
Signal Ratio
A. halleri
Fold Change
2Zn(II)
A. halleri
Fold Change
1Zn(II),
Root: 2 h,
Shoot: 8 h
A. thaliana
Fold Change
1Zn(II),
Root: 2 h,
Shoot: 8 h
Annotation
Root CandidatesCytoplasmic metal influx
262546_at At1g31260 ZIP10a 141 111** 3.74 ZIP family of metaltransporters
253413_at At4g33020 ZIP9 617 36.2*** ZIP family of metaltransporters
260462_at At1g10970 ZIP4 722 24.6*** 2.29 0.31 ZIP family of metaltransporters
259723_at At1g60960 IRT3a 738 4.47* 2.79 0.28 0.37* ZIP family of metaltransporters
257715_at At3g12750 ZIP1 194 4.08* 2.05 ZIP family of metaltransporters
266336_at At2g32270 ZIP3a 849 1.71* 0.37 0.30* ZIP family of metaltransporters
Root-to-shoot metal transport258646_at At3g08040 FRD3a 2,890 44.6*** Multidrug and toxin
efflux family transporterOther membrane transport
266184_s_at At2g38940 PHT1;4a 822 19.1** Phosphate:H1 symporterfamily, ScPho84 like
253658_at At4g30120 HMA3 457 9.14* P1B CPX metal ATPasefamily
267266_at At2g23150 NRAMP3 685 6.61*** Natural resistance-associated macrophageprotein family
249255_at At5g41610 CHX18a 34.1 6.52* Cation/H1 exchanger ofunknown specificity,CPA2 family
up-regulation was observed in response to excessZn (0% down, 3.2% up, P , 0.1; Supplemental TableVI). By comparison, in shoots of A. halleri a smallerproportion of expressed genes responded, and tran-scriptional down-regulation was predominant (0.8%down, 0.2% up, P , 0.1; Supplemental Table IV). Thisindicated that primarily in the shoots, the two speciesrespond differently to excess Zn also at the wholetranscriptome level.
Confirmation of Candidate Genes by Real-Time RT-PCR
To assess the reliability of the cross species expres-sion differences between A. halleri and A. thalianadetected using the microarrays, we determined tran-script levels by an alternative method, quantitativereal-time RT-PCR, for a subset of 18 genes (Table II).Overall, there was good qualitative agreement be-tween the ratios of microarray signals in A. hallerirelative to A. thaliana and the ratios of transcript levels
as determined by real-time RT-PCR (Table II). The falsenegatives according to microarray hybridization,ZIP3, and HMA4 (in roots), could be attributed to ahigh divergence of the A. halleri sequence from the A.thaliana sequence in the region covered by the oligo-nucleotide probes on the microarray, in combinationwith a modest magnitude of the A. halleri to A. thalianaratio of transcript levels. For several other genes,quantitative discrepancies between microarray andreal-time RT-PCR data are likely to be attributable tocross hybridization between high transcript levels of agene in A. halleri and the probes for a closely relatedmember of the same gene family (NAS2 on NAS4 inroots, NAS3 on NAS2 in shoots, NAS2 on NAS3 inshoots under Zn deficiency, and PDI1 on PDI2). It wasreported earlier that very low signals (for example, forZIP10, MTP8, and ZIP9 in A. thaliana) are often quan-tified imprecisely on the ATH1 chip (Czechowski et al.,2004), and for these genes in particular, the real-timeRT-PCR results are likely to be more accurate.
Table I. (Continued from previous page.)
ATH1 Probe
Set Identifier
AGI Locus
Identifier
Gene
Name
A. halleri
Raw Signal
A. halleri/
A. thaliana
Signal Ratio
A. halleri
Fold Change
2Zn(II)
A. halleri
Fold Change
1Zn(II),
Root: 2 h,
Shoot: 8 h
A. thaliana
Fold Change
1Zn(II),
Root: 2 h,
Shoot: 8 h
Annotation
Other membrane transport253658_at At4g30120 HMA3 796 358** P1B CPX metal ATPase
aNewly identified candidate gene. bData for 259632_at, NAS4, and 250832_at, NAS1, in roots are given for comparison. cIdentified asfalse positive by real-time RT-PCR (Table II). dMirouze et al. (2006).
Talke et al.
154 Plant Physiol. Vol. 142, 2006
https://plantphysiol.orgDownloaded on May 29, 2021. - Published by Copyright (c) 2020 American Society of Plant Biologists. All rights reserved.
Table II. Summary of real-time RT-PCR data for selected candidate genes at control Zn concentrations
Transcript levels relative to EF1a of selected candidate genes as determined by real-time RT-PCR are given for shoot and root tissues from plants grown under control conditions. Values are meanand SE of between three and five independent experiments for A. halleri and A. thaliana. The statistical significances of the differences between A. halleri and A. thaliana expression levels wasassessed with a Student’s t test (two sample assuming unequal variances), P values are indicated as follows: *5 P, 0.05; **5 P, 0.01; ***5 P, 0.001. In real-time RT-PCR experiments, geneswere determined as not expressed (n.e.) when CT values were above a threshold of 35 and reaction efficiencies with the respective cDNAwere more than 5% below the mean RE for the respectiveprimer pair; n.d. 5 not determined.The normalized -fold difference between A. halleri and A. thaliana from ATH1 microarray data is given for comparison (see Table I for further details). Flagsdenote calls as given by MAS5: P 5 present; M 5 marginal; A 5 absent. NAS1 and NAS4 were included for comparison with the other members of the NAS family.
Gene
Name
Root Shoot Tissue with Higher
Expression in A. halleri
Compared to A. thalianaReal-Time RT-PCR Microarray Real-Time RT-PCR Microarray
A. halleri
Transcript
Level
A. thaliana
Transcript
Level
A. halleri/
A. thaliana
Transcript
Ratio
A. halleri/
A. thaliana
Signal Ratio
A. halleri
Flag
A. thaliana
Flag
A. halleri
Transcript
Level
A. thaliana
Transcript
Level
A. halleri/
A. thaliana
Transcript
Ratio
A. halleri/
A. thaliana
Signal Ratio
A. halleri
Flag
A. thaliana
Flag
Real-Time
RT-PCRMicroarray
Cytoplasmic metal influxZIP3 498 6 46.0 ** 108 6 46.1 4.62 1.71* P P 10.2 6 2.62 * 1.29 6 0.43 7.89 2.37 A A Root,
shootZIP4 155 6 19.6 *** 13.1 6 1.74 11.9 24.6*** P A 52.6 6 7.96 92.1 6 44.3 0.57 4.04* P P Root Root,
shootZIP9 654 6 189 * 50.0 6 4.61 13.1 36.2*** P P, A n.e. n.e. n.d. 6.32 A A Root RootZIP10 17.9 6 4.98 * 2.89 6 0.71 6.18 111** P, M A n.e. n.e. n.d. 3.24 A A Root RootIRT3 191 6 20.5 *** 40.9 6 14.6 4.68 4.47* P P 85.0 6 4.70 63.8 6 38.9 1.33 1.70* P P Root RootZIP6 122 6 27.9 * 26.3 6 4.78 4.64 2.16* P P 183 6 33.5 ** 18.2 6 5.35 10.1 7.34*** P M, A Root,
shootShoot
Root-to-shoot metal transportHMA4 2,092 6 401 ** 327 6 29.2 6.39 0.71 P P 1,141 6 332 * 37.9 6 31.4 30.1 52.5** P A Root,
shootShoot
FRD3 292 6 87.8 * 18.8 6 4.95 15.5 44.6*** P P 46.5 6 13.3 * 7.45 6 3.05 6.23 37.1* P A Root,shoot
Root
Other membrane transportMTP8 10.6 6 3.53 * 1.39 6 0.53 7.65 15.9* A A 32.4 6 7.84 * 2.86 6 1.45 11.3 65.0** P, A A Root,
shootShoot
PHT1;4 77.6 6 15.9 * 18.7 6 1.43 4.15 19.1** P P 108 6 10.2 *** 20.2 6 10.3 5.33 446* P A Root,shoot
Root
Nicotianamine/thiol biosynthesis pathwaysNAS1 52.0 6 10.4 49.5 6 24.9 1.05 0.81 P P 37.4 6 9.27 79.3 6 26.7 0.47 0.54 P, M PNAS2 437 6 113 * 27.5 6 14.9 15.9 48.4*** P P n.e. n.e. n.d. 3.14 P, A A Root RootNAS3 13.9 6 2.55 ** 0.20 6 0.05 71.4 5.22** P A 44.7 6 5.06 * 15.8 6 6.43 2.82 7.89** P P Root,
shootRoot,shoot
NAS4 12.3 6 1.54 6.93 6 2.12 1.78 13.7 P, M A 34.5 6 10.4 * 125 6 34.6 0.28 0.36* P POASA2 22.7 6 5.28 83.5 6 27.1 0.27 3.25* P P 13.3 6 2.08 64.7 6 30.8 0.21 10.6*** P P ShootSAMS2 881 6 140 ** 158 6 36.4 5.59 2.60* P P 422 6 65.4 ** 70.6 6 24.6 5.97 8.54*** P P Root,
shootShoot
Stress protectionPDI1 362 6 56.7 * 173 6 42.9 2.10 3.50** P P 355 6 63.3 * 91.3 6 47.1 3.89 39.5** P P Root,
shootShoot
PDI2 118 6 20.4 125 6 48.8 0.94 3.08* P P 129 6 16.9 127 6 49.3 1.01 17.5** P P Shoot
Out of all candidate genes identified here, absolutetranscript levels were highest for A. halleri HMA4 (seeTable II), which encodes a membrane transport proteinof the P1B transition metal pump family of the P-typeATPases (Axelsen and Palmgren, 2001). To determinethe cellular function of this candidate gene, which hadnot previously been implicated in metal hyperaccu-mulation or metal tolerance in A. halleri, we expressedAhHMA4 and AtHMA4 in metal hypersensitive mu-tants of budding yeast (Saccharomyces cerevisiae; Fig. 4).Both the A. thaliana and the A. halleri HMA4 proteinswere similarly able to complement Zn hypersensitivityof the zrc1 cot1 double mutant and Cd hypersensitivityof the ycf1 mutant of yeast (Fig. 4, A and B). When thetransport functions of AtHMA4 and AhHMA4 weredisrupted by converting the conserved Asp residues(D401) of the phosphorylation motifs (Axelsen and
Palmgren, 2001) into Alas, respectively, complemen-tation was no longer observed. Consequently, Zn andCd are substrates of both Arabidopsis HMA4 proteinsin yeast, and the observed complementation is attrib-utable to the transport functions of the proteins andnot to the known ability of the cytoplasmic C terminiof AtHMA4 and related proteins to bind metal ions(Bernard et al., 2004; Papoyan and Kochian, 2004).
Gene Copy Number of Highly Expressed CandidateGenes in A. halleri
Previously, high transcript levels of the Zn tolerancegene MTP1 were partly attributed to the presence ofseveral MTP1 gene copies in the genome of A. halleri,whereas the genome of the closely related nontolerantspecies A. thaliana contains only a single MTP1 genecopy (Drager et al., 2004). To estimate the genomiccopy number for other candidate genes in A. halleri, wegenerated genomic DNA blots for 12 of these genes(Kim et al., 1998; Marquardt et al., 2000; Small andWendel, 2000; Ishizaki et al., 2002). Complex bandpatterns suggested that more than one gene copy ispresent in the A. halleri genome for the genes ZIP3,ZIP6, ZIP9, and HMA4 (Fig. 5, A to D; Table III;Supplemental Fig. 1; for data evaluation see also‘‘Materials and Methods’’). For comparison, the max-imum expression of these genes was 8-, 10-, 13-, and30-fold higher in A. halleri than in A. thaliana, respec-tively, as determined by real-time RT-PCR (Table II).By contrast, the simple band patterns of the genomicDNA blots indicated that a number of other candidategenes are most likely to be present as single copies inthe A. halleri genome. These include FRD3, ZIP4,MTP8, ZIP10, IRT3, PHT1;4, and PDI1, the maximumexpression of which was 16-, 12-, 11-, 6-, 5-, 5-, and4-fold in A. halleri compared to A. thaliana, respectively(Tables II and III; Fig. 5). Among the group of candi-date genes for which DNA gel-blot analysis was per-formed, more than single gene copies were detectedfor each of the three candidate genes exhibiting thehighest transcript levels in A. halleri (normalized toEF1a transcript levels), i.e. HMA4 and ZIP9 and ZIP3(compare Table II, left section). Based on sequence dataand restriction analysis of genomic PCR products, it isunlikely that additional bands on the DNA gel blotsoriginated from uncharacterized restriction sitesrather than additional gene copies in the A. hallerigenome (see Table III; M. Hanikenne and U. Kramer,unpublished data).
Regulation of Candidate Gene Expression
Transcriptional Zn responses were further investi-gated for a subset of 18 candidate genes using real-time RT-PCR (Fig. 6). In the roots of A. halleri, short-termexposure to 300 mM Zn for 2 and 8 h resulted in thedown-regulation of transcript levels of a group ofgenes, which encode several ZIP family Zn trans-porters with likely functions in Zn influx into the
Figure 4. Functional analysis of AtHMA4 and AhHMA4 cDNAs inyeast mutant strains. Expression of Arabidopsis HMA4 cDNAs in theZn-hypersensitive zrc1 cot1 double mutant (A) and in the Cd-hyper-sensitive ycf1 mutant (B) is shown. Wild-type cells were transformedwith pFL38H-GW (ev), and mutant cells were transformed withpFL38H-GW (ev) and with wild-type and mutated (D401A) cDNAs ofboth AtHMA4 and AhHMA4 in pFL38H-GW. The phosphorylation ofthe Asp residue (here D401) of the conserved phosphorylation motif isan indispensable step of the reaction cycle in all P-type ATPases. Serialdilutions of transformants were spotted on LSP plates containing theindicated concentrations of ZnSO4 (A) or CdSO4 (B). The 100 dilutioncorresponds to an OD600 of 0.5. Plates were incubated at 30�C for 4 d.
Talke et al.
156 Plant Physiol. Vol. 142, 2006
https://plantphysiol.orgDownloaded on May 29, 2021. - Published by Copyright (c) 2020 American Society of Plant Biologists. All rights reserved.
cytoplasm (Fig. 6A). Down-regulation of transcriptlevels of this group of genes was clearly less pro-nounced in the roots of A. thaliana after short-termexposure to 30 mM Zn. It has to be kept in mind that inroots of A. halleri the Zn-induced decrease in transcriptlevels occurred from very high transcript levels undercontrol conditions. Consequently, despite a more pro-nounced down-regulation in A. halleri, final root tran-script levels 8 h after supply of excess Zn weregenerally higher in A. halleri than in A. thaliana (com-pare Fig. 2).
Different from roots, shoot transcript levels of mostZIP genes were generally as low in A. halleri as in A.thaliana under control conditions (see Fig. 2; Table II).In shoots of A. halleri, there was virtually no transcrip-tional response to a short-term oversupply of Zn,whereas A. thaliana shoots responded by a markeddown-regulation of several genes including primarilyZIP genes and PHT1;4 (Fig. 6B). To test whether theshoots of A. halleri were able to respond transcription-
ally to changes in Zn supply, we investigated geneexpression under Zn deficiency (Fig. 6). Zn deficiencyresulted in a very pronounced increase in shoot tran-script levels of genes encoding ZIP family transportersand nicotianamine synthase proteins (Fig. 6B; seeTable I). In roots of A. halleri, Zn deficiency resultedin a moderate increase in transcript levels of thesegenes, when compared to the high transcript levelsobserved under control conditions (Fig. 6A; for com-parison with A. thaliana, see Figs. 2 and 7).
To analyze the specificity of transcriptional Zn re-sponses, candidate gene regulation in A. halleri wasalso investigated following exposure to moderateexcesses of Cd (30 mM) or copper (Cu; 5 mM) andNaCl (100 mM) for 2 and 8 h (Fig. 6). The resultssuggested that overall, Zn responses were rather spe-cific. Some of the Zn-responsive genes, primarily NASand ZIP genes, were also down-regulated in responseto excess Cu, as for excess Zn. The transcriptionalresponses to excess Zn and excess Cd, however, were
Figure 5. Genomic DNA gel blot for ZIP9 (A), HMA4 (B), FRD3 (C), and MTP8 (D) in A. thaliana and in A. halleri. Blots wereperformed using genomic DNA extracted from A. thaliana (Col accession) and from two A. halleri individuals (Lan 3-1 and Lan 5of the accession Langelsheim). A to D, Five micrograms of genomic DNAwere digested with EcoRI,HindIII, orNcoI, resolved ona 0.9% (w/v) agarose gel, blotted, and hybridized with radiolabeled A. halleri-derived cDNA probes. A, An EcoRI restriction siteis present in the region spanned by the ZIP9 probe in A. thaliana only, aHindIII restriction site is present in the intron spanned bythe probe in A. halleri only, and an NcoI restriction site is present in an exon spanned by the probe in both species. B, A HindIIIrestriction site is present in one of the introns spanned by the HMA4 probe in A. halleri. C, A HindIII restriction site is present inthe region spanned by the FRD3 probe in both species. D, After restriction withNcoI, a predicted 20.9 kb fragment in A. thalianawas not detected by theMTP8 probe, probably because of the large size of this fragment (see also Table III). Note that for clarity,lanes were rearranged in the ZIP9 blot.
Zinc Homeostasis in Arabidopsis halleri and A. thaliana
Plant Physiol. Vol. 142, 2006 157
https://plantphysiol.orgDownloaded on May 29, 2021. - Published by Copyright (c) 2020 American Society of Plant Biologists. All rights reserved.
distinct. Furthermore, there was a second group ofgenes exhibiting a pattern of regulation that wasclearly different from the Zn- and Cu-repressed genes.Expression of these genes was induced in response toexposure to 100 mM NaCl and included MTP8,PHT1;4, FRD3, ZIP6, and OASA2. Transcript levels ofsubsets of these genes were also increased following
exposure to excess Cu or Cd, supporting their respon-siveness to abiotic stress. PDI transcript levels showedfew changes over the entire range of conditions, with anincrease in root PDI2 transcript levels in response toexcess Cu in A. halleri and an increase in PDI1transcript levels in response to Zn in shoots of A.thaliana.
Table III. Gene copy number estimated by genomic DNA gel blot for selected candidate genes in A. halleri
Five micrograms of genomic DNA isolated from A. thaliana (Col-0 accession) and from two A. halleri individuals (Lan 3-1 and Lan 5 of theLangelsheim accession) were digested with EcoRI, HindIII, or NcoI, resolved on a 0.9% (w/v) agarose gel, blotted, and hybridized with radiolabelledA. halleri-derived cDNA probes. The size of the detected DNA fragments is given in kilobase pair (detectable fragments were in the range of 0.2 and14 kb). Note that in A. halleri, there are a maximum of two alleles per gene copy. Blots for the four underlined genes are presented in Figure 5, andother blots are shown in Supplemental Figure 1.
Gene Copy No. in A. halleriA. halleri Lan 3.1 A. halleri Lan 5 A. thaliana Col-0
aCross hybridization of AhZIP3 probe with AtZIP5 genomic sequence. bCross hybridization of AhIRT3 and AhZIP4 probes with AtZIP4 andAtIRT3 genomic sequences, respectively. cNot detected, predicted fragment given for A. thaliana. dEnzyme cuts once in one of the intronsspanned by the cDNA probe. eEnzyme cuts once in the probe. fEnzyme cuts twice in the probe. gCross hybridization of AhHMA4 probewith AtHMA2 genomic sequence. hCross hybridization of AhPht1;4 probe with At3g54700 genomic sequence. iCross hybridization ofAhPDI1 and AhPDI2 probes with AtPDI2 and AtPDI1 genomic sequences, respectively. Note that it is likely that cross hybridization also occurs ofAhPDI1 and AhPDI2 probes with AhPDI2 and AhPDI1 genomic sequences. jAdditional band, identity not known.
Talke et al.
158 Plant Physiol. Vol. 142, 2006
https://plantphysiol.orgDownloaded on May 29, 2021. - Published by Copyright (c) 2020 American Society of Plant Biologists. All rights reserved.
Figure 6. Transcriptional responses of selected candidate genes to Zn and other ions. Transcript levels relative to EF1a wereassessed in roots (A) and shoots (B) of hydroponically grown plants subjected to short-term treatment with high Zn, long-term Zndeficiency, or short-term treatment with Cd, Cu, or Na (see ‘‘Materials and Methods’’ for details). RTLs within a treatment werenormalized to the RTL of the control treatment in the respective experiment (‘‘RTL normalized to control’’). To determinetranscriptional responses for each treatment and each gene, the mean RTL normalized to control was calculated from therespective values from two independent biological experiments. Bold letters indicate genes that exhibit higher transcript levels inA. halleri compared to A. thaliana after cultivation in the presence of control Zn concentrations. Note that some genes were notexpressed in some conditions, see Figures 2 and 7 for details.
Zinc Homeostasis in Arabidopsis halleri and A. thaliana
Plant Physiol. Vol. 142, 2006 159
https://plantphysiol.orgDownloaded on May 29, 2021. - Published by Copyright (c) 2020 American Society of Plant Biologists. All rights reserved.
Transcript levels of HMA4 appeared to be constitu-tively very high in the roots of A. halleri across alltreatments, and only slightly induced under Zn defi-ciency in A. halleri shoots (Fig. 6, A and B). A detailedanalysis confirmed the extremely high expression ofHMA4 in roots and shoots ofA. halleri (Fig. 7). Under Zndeficiency, shoot HMA4 transcript levels were slightlyincreased, and they were reduced to control levelsshortly after readdition of Zn to the medium (Fig. 7).
We investigated transcript levels of the additionalnewly identified candidate genes ZIP3 and IRT3, un-der different conditions of Zn deficiency, resupply, andoversupply (Fig. 7). There were overall similarities inshort-term Zn responses of plants maintained undercontrol (1 mM added Zn) and Zn-deficient (0 mM addedZn) conditions (Fig. 7; see also Fig. 2). In response toresupply of 5 mM Zn to Zn-deficient plants, there was adecrease in transcript levels for these Zn-deficiency-
induced genes in both A. halleri and A. thaliana. Thisdemonstrated clearly that there is no general defect inthe magnitude or time scale of transcriptional Znresponses in A. halleri (Fig. 7; see also Fig. 2: ZIP4and ZIP9 in controls and at time points 8 and 24 h). Inroots of A. halleri, after the initial down-regulation,which was strongest 2 h after resupply of Zn, a generaltrend of a gradual increase in transcript levels wasobserved at later time points (Fig. 7).
DISCUSSION
In a microarray-based cross species comparativetranscript analysis in A. thaliana and A. halleri, weidentified a set of candidate genes expressed at higherlevels in A. halleri than in A. thaliana or regulated by Znin A. halleri (Fig. 3; Table I; Supplemental Tables I–IV).
Figure 7. Zn-dependent transcriptional regulation of selected candidate genes. Real-time RT-PCR analysis of expression ofHMA4, ZIP3, and IRT3 genes in roots (top section) and shoots (bottom section) of hydroponically grown A. halleri (black bars)and A. thaliana (light gray bars). Two different basal Zn conditions upon which additional Zn was supplied were assessed: ShortZn oversupply and Zn deficiency/resupply. For a detailed description, see Figure 2 and ‘‘Materials andMethods.’’ Data shown aretranscript levels relative to EF1a from one experiment representative of two independent biological experiments. Tissues from atleast three culture vessels with three plants each were pooled for each condition. Genes were concluded to be not expressed(n.e.) when CT values were above a threshold of 35 and reaction efficiencies with the respective cDNAweremore than 5% belowthe mean RE for the respective primer pair. Note that ZIP3 was not expressed in shoots only for the condition A. thaliana, 30 mM
Zn supply, 8 h. At all other conditions, ZIP3 transcripts were detectable but RTLs were very low (A. halleri: 4.3 [control], 3.5 [2 h],and 2.6 [8 h]; A. thaliana: 1.7 [control], 1.7 [2 h], and 3.9 [control Zn deficiency]), and therefore do not appear as visible bars inthe diagram.
Talke et al.
160 Plant Physiol. Vol. 142, 2006
https://plantphysiol.orgDownloaded on May 29, 2021. - Published by Copyright (c) 2020 American Society of Plant Biologists. All rights reserved.
Microarray data were confirmed by an alternativetechnique of transcript quantification, real-time RT-PCR, for a representative subset of genes. Among thegenes identified using microarrays, combined in silicoand real-time RT-PCR analyses revealed a total of 29genes encoding putative metal homeostasis proteins,which we consider here as an initial set of majorcandidates for a role in naturally selected metal hyper-accumulation and associated metal hypertolerance inA. halleri. Of these, 18 genes have not previously beenimplicated in these traits, and 11 confirm earlier find-ings (Tables I and II; Becher et al., 2004; Weber et al.,2004; Mirouze et al., 2006).
We performed DNA gel blots to estimate genomiccopy numbers for a subset of 12 candidate genes (Fig.5; Table III; Supplemental Fig. 1). The obtained data areimportant for future studies of segregation of candi-date genes and for the cloning of genomic sequencesfrom A. halleri. Moreover, the combination of DNA gelblot with steady-state transcript level data suggestedthe presence of more than a single genomic copy inA. halleri of the genes, for which the highest transcriptlevels were detected, HMA4, ZIP9, and ZIP3 (Tables IIand III). As reported previously for MTP1 (Drageret al., 2004), increased gene copy numbers of thesegenes may contribute to the ability of A. halleri toexpress these genes at very high levels, and possiblyallow for regulatory diversification among the copiesof these genes. However, high expression of thesecandidate genes cannot be explained without addi-tionally invoking an altered regulation of transcriptlevels. Moreover, substantially higher transcript levelsin A. halleri than in A. thaliana, as also found for FRD3,ZIP4, and MTP8, were not in all cases associated withindications for more than a single genomic copy of therespective candidate gene.
Ratios of shoot-to-root Zn concentrations are gener-ally below unity in A. thaliana and other nonaccumu-lator plants, and typically above unity in A. halleri andother Zn hyperaccumulating plants (Baker et al., 1994;Dahmani-Muller et al., 2000; Becher et al., 2004; Weberet al., 2004; Fig. 1). Zn concentrations partitioned intoshoots relative to the roots were dependent on theexternal Zn supply (Fig. 1). In addition, the adjustmentof Zn partitioning as a function of Zn supply differedbetween A. thaliana and A. halleri, with maximumrelative Zn concentrations partitioned into the shoot atapproximately 5 mM in A. halleri, and under Zn defi-ciency in A. thaliana. In A. halleri the partitioning of Zninto the shoot was restricted at high external Znsupply of 30 mM Zn and above. This may contributeto Zn hypertolerance in the Langelsheim accession ofA. halleri. It has been reported that individuals fromthe most Zn-hypertolerant populations of A. halleri,which include the Langelsheim accession, accumulatelower leaf Zn concentrations than individuals fromless hypertolerant populations (Bert et al., 2000, 2002).As expected based on Zn partitioning, we observedbetween-species differences in gene expression undercontrol conditions, as well as species-specific short-
term transcriptional Zn responses of candidate genesin A. halleri and A. thaliana (Tables I and II; Figs. 1, 2, 6,and 7).
Transcript levels of genes known to be transcrip-tionally regulated by Zn, such as ZIP1, ZIP3, ZIP4,ZIP9, NAS2, and NAS3 can be used as indicators of theZn status of A. thaliana plants (Grotz et al., 1998; Wintzet al., 2003). The expression of a group of Zn deficiencymarker genes suggested that both roots and shoots ofA. thaliana are largely Zn sufficient at 1 or 5 mM Zn, andconcertedly experience Zn deficiency after prolongedgrowth in a medium lacking added Zn (Figs. 2 and 7).Similarly, the expression of the same set of genes in A.halleri plants grown at 1 or 5 mM Zn suggests that theshoots are Zn sufficient under these conditions (seeNAS2, ZIP9, and ZIP4 transcript levels in shoots inFigs. 2, 6, and 7). By contrast, the roots of A. halleri acttranscriptionally as Zn deficient after prolongedgrowth at 1 or 5 mM Zn. Thus, transcript levels of acommon set of Zn status marker genes can be inter-preted to indicate that roots and shoots of A. halleriexperience different Zn stati in the steady state.
There was a rapid and strong down-regulation oftranscript levels of ZIP3, ZIP4, and IRT3 in roots ofA. halleri following exposure to excess Zn or resupplyof Zn to roots of Zn-deficient plants grown in amedium lacking added Zn (Figs. 2 and 7). Overall,this suggests that short-term Zn-dependent down-regulation of Zn deficiency-induced transcripts inroots of A. halleri is generally as functional and efficientas in A. thaliana. Only NAS2 and ZIP9 transcript levelswere not or only partly Zn responsive in roots of A.halleri, respectively. The Zn-dependent control of tran-script abundance of these two genes in A. thalianaappears to be partially or fully inactive in roots, but notin shoots, of A. halleri. The simplest explanation for thehigh expression of Zn deficiency responsive genes inthe roots of A. halleri at normal Zn supplies could be aloss of Zn from the roots via high root-to-shoot Zntransport rates (see Fig. 1). The proteins encoded bytwo of the candidate genes have been implicated inroot-to-shoot metal translocation in A. thaliana: HMA4and FRD3.
The A. thaliana HMA4 protein is a P1B-type metalpump that localizes to the plasma membrane andmediates cellular metal efflux. Together with the re-lated AtHMA2, AtHMA4 has a role in the transport ofZn from the root to the shoot and in metal detoxifica-tion (Hussain et al., 2004; Verret et al., 2004; Mills et al.,2005). According to real-time RT-PCR analysis (Fig. 7;Table II), HMA4 transcript levels are between 4- and10-fold higher in roots and at least 30-fold higher inshoots of A. halleri than in A. thaliana. In A. halleri,HMA2 transcript levels are much lower than HMA4transcript levels (less than 0.25% in shoots and lessthan 0.2% in roots), and HMA2 transcript levels arelower in A. halleri than in A. thaliana (approximately25% in shoots and 8% in roots; I.N. Talke, L. Krall, andU. Kramer, unpublished data). Heterologous expres-sion in yeast indicated that AhHMA4 and AtHMA4
Zinc Homeostasis in Arabidopsis halleri and A. thaliana
Plant Physiol. Vol. 142, 2006 161
https://plantphysiol.orgDownloaded on May 29, 2021. - Published by Copyright (c) 2020 American Society of Plant Biologists. All rights reserved.
are similarly able to complement Zn and Cd hyper-sensitivity of yeast mutants (Fig. 4), and that comple-mentation is dependent on the transport functions ofthe proteins. Together, our data suggest that AhHMA4is functionally highly similar to AtHMA4, and thatHMA4 expression is very high in A. halleri and pre-dominant over the expression of HMA2. Thus, in thecontext of the data published on AtHMA4, the datapresented here implicate A. halleri HMA4 in Zn and Cdhypertolerance and in Zn hyperaccumulation, i.e. thepartitioning of Zn predominantly to the shoot.
The A. thaliana and A. halleri FRD3 proteins aremembers of the multidrug and toxin efflux family ofmembrane transport proteins. Although the transportfunction of AtFRD3 has not yet been directly estab-lished, the available data suggest that AtFRD3 has arole in maintaining iron (Fe) mobility between the rootand the shoot. In the frd3 mutant, root Fe deficiencyresponses are constitutively active, leading to second-ary accumulation of a number of metals includingmanganese (Mn), cobalt, Cu, Zn, and Fe in both rootsand shoots (Delhaize, 1996; Rogers and Guerinot, 2002;Lahner et al., 2003; Green and Rogers, 2004). Basedon published data, one may thus expect FRD3 to beexpressed at low levels in hyperaccumulator plants.Unexpectedly however, FRD3 is expressed at very highlevels in A. halleri (Tables I and II; Fig. 6). In agreementwith this, Zn hyperaccumulation in A. halleri is specificboth in terms of the accumulated metal and in terms ofthe accumulating tissue, which is the shoot (Becheret al., 2004). In addition, Zn accumulation is highest inthe vacuoles of leaf mesophyll cells in A. halleri (Kupperet al., 2000; Sarret et al., 2002), whereas in the A. thalianafrd3 mutant Fe and Mn appear to accumulate in the leafapoplast (Green and Rogers, 2004). In conclusion, basedon the available evidence, it can be hypothesized thatthe high expression of FRD3 in A. halleri contributesto metal homeostasis, but not specifically to the highaccumulation of Zn in shoots of A. halleri.
In addition to FRD3, several other genes found to behighly expressed in A. halleri (Table I) were previouslyassociated not with Zn or Cd homeostasis, but with thehomeostasis of other transition metals, among themFER1 (Fe), FER2 (Fe), NRAMP3 (Fe and Mn), IREG2(Fe), HMA6/PAA1 (Cu), MTP8, and MTP11 (Mn;Thomine et al., 2000; Petit et al., 2001; Delhaize et al.,2003; Shikanai et al., 2003; Thomine et al., 2003; McKieand Barlow, 2004; Lanquar et al., 2005). The encodedproteins may contribute to the tolerance and hyper-accumulation phenotypes by adjusting the homeosta-sis of Fe, Mn, and Cu to the alterations in Zn and Cdhomeostasis of A. halleri. Alternatively, these proteinsmay have direct roles in Zn or Cd homeostasis. It willbe important to establish the functions of these pro-teins in the future.
To address the specificity of Zn-dependent regula-tion of candidate gene transcript levels, we exposedplants to an excess of Cu, Cd, or NaCl in the hydro-ponic medium. Root ZIP3 and ZIP4 transcript levelsand transcript levels of all four NAS genes were
suppressed following exposure of A. halleri to excessCu (Fig. 6). Wintz et al. (2003) demonstrated that ex-pression of A. thaliana ZIP4 can complement the Cu up-take defect of the ctr1 yeast mutant. AtZIP3-dependentuptake of 65Zn21 into yeast cells was decreased by morethan 55% in the presence of a 10-fold excess of CuCl2(Grotz et al., 1998). Root-to-shoot transport of Cu in thexylem is likely to occur as a Cu-nicotianamine complex(Pich et al., 1994; Pich and Scholz, 1996; von Wirenet al., 1999). The Cu-induced down-regulation of ZIPand NAS transcript levels observed here in A. halleri(Fig. 6) may constitute a regulatory response counter-acting excessive Cu uptake and transport to the shoot.Based on the chemical similarity between Cd21 andZn21 ions, Cd toxicity may involve the disruption ofplant Zn homeostasis. However, none of the ZIP orNAS genes known to have a role in Zn acquisition anddistribution were regulated in response to an excess ofCd (Fig. 6). This may reflect the Cd hypertolerance ofA. halleri (Bert et al., 2003).
There appear to be similarities between major can-didate genes expressed at high levels in A. halleri andin the metal-tolerant Zn and Cd hyperaccumulatormodel species T. caerulescens. The A. halleri/A. lyrataclade is phylogenetically closest to A. thaliana, withan estimated divergence time of between 3.5 and5.8 million years ago (Koch et al., 2001; Yogeeswaranet al., 2005). For comparison, T. caerulescens belongs tothe clade containing Sinapis alba, which is thought tohave diverged from the A. thaliana lineage more than20 million years ago (Gao et al., 2005). In contrast toA. halleri, multiple genomic copies have not beenreported for any T. caerulescens candidate genes sofar. High transcript levels have been observed for T.caerulescens ZNT1 (Lasat et al., 2000; Pence et al., 2000;Assuncxao et al., 2001) as well as for its sequenceortholog AhZIP4 (Becher et al., 2004), and for TcZNT2(Assuncxao et al., 2001) and its sequence orthologAhIRT3 (Fig. 7), which are likely to encode cellularZn uptake systems. The sequence-orthologous TcZTP1(Assuncxao et al., 2001), Thlaspi goesingense MTP1(Persans et al., 2001), and AhMTP1 transcripts (Becheret al., 2004; Drager et al., 2004) have all been reportedto occur at very high levels in the respective speciesand encode proteins mediating metal efflux from thecytoplasm. Finally, high transcript levels have beenreported for TcHMA4 (Bernard et al., 2004; Papoyanand Kochian, 2004) and AhHMA4 (Fig. 7), and the en-coded proteins are capable of mediating cellular Zn andCd detoxification (see Fig. 4; Papoyan and Kochian,2004). These observations raise the question ofwhether metal hyperaccumulation and metal hyper-tolerance are ancestral traits that have subsequentlybeen lost in the majority of Brassicaceae species. Wefavor the alternative hypothesis that Zn hyperaccu-mulation evolved several times independently in theBrassicaceae through recent mutations leading to theconvergent overexpression of an overlapping set ofgenes responsible for the Zn hyperaccumulation andhypertolerance traits.
Talke et al.
162 Plant Physiol. Vol. 142, 2006
https://plantphysiol.orgDownloaded on May 29, 2021. - Published by Copyright (c) 2020 American Society of Plant Biologists. All rights reserved.
The comprehensive identification of candidategenes in several hyperaccumulator model plants is aprerequisite for their functional, regulatory, andsequence analysis. This will be an important steptoward a better understanding of the evolution ofnaturally selected metal hyperaccumulation and asso-ciated metal hypertolerance, and of evolutionary ge-nome dynamics in plants. Furthermore, this will allowa comparison with extreme traits in other Arabidopsisrelative model systems, such as sodium (Na) tolerancein Thellungiella halophila (Inan et al., 2004; Taji et al.,2004; Gong et al., 2005).
In summary, we have used cross species transcriptprofiling to successfully identify on a genome-widescale a number of candidate genes for naturally se-lected metal tolerance and hyperaccumulation in A.halleri. By exploring candidate genes at the genome,functional, and regulatory level we have obtained in-sights into the complex molecular basis of an extremephysiological trait.
MATERIALS AND METHODS
Plant Material, Growth Conditions,
and Experimental Treatments
Plants of Arabidopsis thaliana (L. Heynhold, accession Columbia [Col]) and
Arabidopsis halleri (L. O’Kane and Al-Shehbaz subsp. halleri, accession Langel-
sheim) were grown from seeds, which were the progeny of a single mother
plant in a pool of six individuals grown in an open greenhouse from seeds
collected at Langelsheim. Plants were cultivated in hydroponic culture as
described in Becher et al. (2004), with a photoperiod of 11 h light, 13 h dark, at
a photon flux density 145 mmol m22 s21 during the day, and day and night
temperatures of 20�C and 18�C, respectively. In all experiments, each culture
vessel of 400 mL hydroponic solution contained three individual plants, and
all treatments were harvested at the same time. Experiments for subsequent
RNA extraction were harvested by carefully separating roots and shoots of
plants, pooling tissues according to treatment, and freezing in liquid nitrogen.
Harvested tissues were stored at 280�C until further processing.
Short-term high Zn, Cd, Cu, and Na treatments were initiated when plants
were 6.5 weeks old. For short-term high Zn supply, the medium was replaced
with fresh control medium (1 mM ZnSO4) 8 h before harvest. Then nothing
(controls) or 30 mM ZnSO4 and 300 mM ZnSO4 for A. thaliana and A. halleri,
respectively, was added to the hydroponic solution 8 or 2 h before harvest.
Three biologically independent experiments, each comprising five culture
vessels per treatment (15 plants per treatment) were conducted. In experi-
ments with Cd, Cu, and Na treatments of A. halleri, media were exchanged
(1 mM ZnSO4; Cd: none; Cu: 0.5 mM CuSO4; Na: 0.1 mM Na2MO4) 3 d before
harvest; 30 mM CdCl2, 5 mM CuSO4, or 100 mM NaCl were added as treatments
8 or 2 h before harvest. Two biologically independent experiments were
conducted, with at least two culture vessels (at least six plants in total) for each
treatment.
For Zn deficiency experiments with Zn resupply, the treatment was
initiated when plants were 4.5 weeks old and lasted for 3 weeks. For
deficiency, ZnSO4 was omitted from the hydroponic solution; the sufficient
control solution contained 5 mM ZnSO4. The last change of medium was done
3 d before harvest. For Zn resupply, 5 mM ZnSO4 was added to the hydroponic
solution of Zn-deficient plants 24, 8, and 2 h before harvest. Two biologically
independent experiments were conducted, each comprising at least three
culture vessels (at least nine plants) per treatment.
In the experiments conducted to determine root and shoot Zn concentra-
tions after long-term treatment with a range of different Zn supplies, plants
were 4.5 weeks old when treatments were initiated. For A. halleri, the
hydroponic solution contained 0, 0.3, 1, 5, 10, 30, or 300 mM ZnSO4, for A.
thaliana 0, 1, 5, or 10 mM ZnSO4. Upon harvest after 3 weeks of treatment,
shoots were briefly rinsed in ultrapure water and blotted dry, and roots from
three individuals of the same species and hydroponic culture vessel were
desorbed together in 150 mL of an ice-cold solution of 5 mM CaCl2 and 1 mM
MES-KOH (pH 5.7) for 20 min, with a replacement of the solution after 5 min,
followed by two washes in 150 mL ice-cold water, each for 3 min. Tissues were
pooled by culture vessel and dried at 60�C for 3 d. Two biologically indepen-
dent experiments were conducted, each comprising three culture vessels (nine
plants) per treatment.
Determination of Zn Concentrations
Dried root and shoot material was homogenized, and for each sample, 35
to 70 mg were processed and elemental contents analyzed by inductively
coupled plasma optical emission spectroscopy using an IRIS Advantage Duo
ER/S (Thermo Jarrell Ash) as described previously (Becher et al., 2004).
RNA Extraction
For microarray expression profiling, total RNA was extracted with Trizol
(Invitrogen Life Technologies) according to the manufacturer’s instructions.
RNeasy spin columns and the RNase-free DNase set for on-column DNase
digest (Qiagen) were used according to the manufacturer’s instructions for
subsequent RNA purification and digestion of genomic DNA, and as an
alternative method of total RNA isolation when transcript levels were deter-
mined by real-time RT-PCR only (RNeasy plant mini kit and RNase-free
DNase set; Qiagen). Quality and quantity of RNA was checked visually by
denaturing gel electrophoresis and by photometric analysis (A260 and A280).
Microarray Expression Profiling and Data Analysis
Synthesis of cDNA, cRNA labeling, and the hybridization on the GeneChip
Arabidopsis ATH1 genome array were done as recommended by the man-