- © 2011 American Society of Plant Biologists. All rights reserved.
Abstract
Abscisic acid-, stress-, and ripening-induced (ASR) proteins were first described about 15 years ago as accumulating to high levels during plant developmental processes and in response to diverse stresses. Currently, the effects of ASRs on water deficit tolerance and the ways in which their physiological and biochemical functions lead to this stress tolerance remain poorly understood. Here, we characterized the ASR gene family from maize (Zea mays), which contains nine paralogous genes, and showed that maize ASR1 (ZmASR1) was encoded by one of the most highly expressed paralogs. Ectopic expression of ZmASR1 had a large overall impact on maize yield that was maintained under water-limited stress conditions in the field. Comparative transcriptomic and proteomic analyses of wild-type and ZmASR1-overexpressing leaves led to the identification of three transcripts and 16 proteins up- or down-regulated by ZmASR1. The majority of them were involved in primary and/or cellular metabolic processes, including branched-chain amino acid (BCAA) biosynthesis. Metabolomic and transcript analyses further indicated that ZmASR1-overexpressing plants showed a decrease in BCAA compounds and changes in BCAA-related gene expression in comparison with wild-type plants. Interestingly, within-group correlation matrix analysis revealed a close link between 13 decreased metabolites in ZmASR1-overexpressing leaves, including two BCAAs. Among these 13 metabolites, six were previously shown to be negatively correlated to biomass, suggesting that ZmASR1-dependent regulation of these 13 metabolites might contribute to regulate leaf growth, resulting in improvement in kernel yield.
Since one-third of the world’s food is produced on irrigated land (Food and Agriculture Organization of the United Nations, 2008), the likely impacts of recurrent heat waves and water stress episodes on global food production are numerous. In the European heat wave of 2003, crop production was reduced by around 30% (Ciais et al., 2005). In maize (Zea mays), the most significant reductions in end-of-season kernel yields are observed when water deficit occurs during the flowering stage, either just before floral initiation or immediately after pollination, because of the impact of water deficit on ovary and young kernel abortions (Claassen, 1970; Boyer and Westgate, 2004). Water deficit during the vegetative growth phase also leads to reductions in overall productivity through reductions in kernel numbers (Boyer and Westgate, 2004). Finally, late-stage water deficit during the kernel-filling period can lead to reductions in yield, depending upon the dry matter reserves in the plant (McPherson and Boyer, 1977).
Trait-based approaches considering maize drought avoidance and dehydration tolerance mechanisms have been relatively slow to progress, as judged by the adoption of improved varieties (Salekdeh et al., 2009). Transgenic maize plants expressing the betA gene from Escherichia coli encoding choline dehydrogenase, a key enzyme in the biosynthesis of Gly betaine from choline, showed enhanced Gly betaine accumulation, resulting in greater kernel yield after drought stress in field tests (Quan et al., 2004). Expression of a maize CAAT box transcription factor, Z. mays nuclear factor Y B subunit 2 (ZmNF-YB2), has also been shown to confer drought tolerance and enhanced photosynthetic capacity under drought stress, with improvements in kernel yield in maize (Nelson et al., 2007). Additionally, two members of a family of bacterial RNA chaperones, cold shock protein A (CspA) from E. coli and CspB from Bacillus subtilis, were shown to confer vegetative tolerance and improved end-of-season kernel yield under water-limiting conditions in maize (Castiglioni et al., 2008). These studies clearly demonstrate that maize plants are amenable to improved water deficit tolerance through multiple mechanisms of action.
Having specific target traits and genes can markedly accelerate progress through the marker-assisted selection of parents and progeny in early generations (Tuberosa and Salvi, 2006). We used a proteomics strategy, with the use of a recombinant inbred line population derived from the cross between the maize inbred lines MBS847 and F2, differing in yield under water deprivation conditions, to identify candidate proteins associated with water stress response in maize (de Vienne et al., 1999). The protein quantity profiling of all individuals in the segregating population made it possible to treat the quantity of each protein as a quantitative trait and to identify a quantitative trait locus (QTL) for protein quantity, hereafter named protein quantity locus (PQL) according to Damerval et al. (1994). When the PQL colocated with a QTL related to water deficit response, a causal relationship could be inferred between the protein level and the trait variation. Furthermore, when the gene encoding the protein also colocated, it could be seen as a candidate gene. Among the identified candidate genes was a member of the abscisic acid-, stress-, and ripening-induced (ASR) gene family, ZmASR1, which colocated with its PQL and QTLs for leaf senescence and anthesis-silking interval (de Vienne et al., 1999; Jeanneau et al., 2002). These associations were validated by the observation that transgenic maize lines misexpressing ZmASR1 showed significant changes in leaf senescence under mild water deficit conditions in the field (Jeanneau et al., 2002). Interestingly, transgenic Arabidopsis (Arabidopsis thaliana) plants overexpressing an ASR gene from Lilium longiflorum (LLA23) or banana (Musa paradisiaca; MpASR) exhibited enhanced capacity to survive in water-limited conditions (Yang et al., 2005; Liu et al., 2010). Overexpression of the tomato (Solanum lycopersicum) ASR1 gene (SlASR1) in tobacco (Nicotiana tabacum) plants also resulted in decreased rates of water loss (Kalifa et al., 2004b). Together, these results suggest that water deficit tolerance is mediated at least in part through ASR proteins.
The ASR proteins are highly charged, hydrophilic, and low-Mr proteins that are widely present in the plant kingdom, except in the Brassicaceae family (for review, see Carrari et al., 2004; Battaglia et al., 2008). Since their discovery more than a decade ago, they have been described as accumulating to high levels during plant developmental processes, such as fruit ripening, pollen maturation, and Glc metabolism (Iusem et al., 1993; Wang et al., 1998; Cakir et al., 2003; Frankel et al., 2007; Yang et al., 2008), and in response to diverse stresses, including water deficit, salt, cold, limited light, and, lately, pathogen attack (Schneider et al., 1997; Riccardi et al., 1998; Vaidyanathan et al., 1999; Huang et al., 2000; Maskin et al., 2001; Kalifa et al., 2004b; Liu et al., 2010; Philippe et al., 2010; Henry et al., 2011). Nevertheless, it is noteworthy that their exact function remains enigmatic. Indeed, the molecular mechanisms underlying the biological roles of the ASR proteins cannot be deduced simply by sequence homology with other known proteins. SlASR1 and the grape (Vitis vinifera) ASR protein, VvMSA (for maturation-, stress-, and abscisic acid-induced protein), have been shown to possess DNA-binding activity (Cakir et al., 2003; Kalifa et al., 2004a; Rom et al., 2006; Maskin et al., 2007; Shkolnik and Bar-Zvi, 2008). On the basis of such findings, these ASR proteins have been proposed to regulate the transcription of sugar- and abiotic stress-regulated genes in fruit and vegetative tissues, respectively (Cakir et al., 2003; Kalifa et al., 2004b; Frankel et al., 2007; Saumonneau et al., 2008). However, a dual function of LLA23 and SlASR1 proteins as a regulator and a protective molecule upon water deficit has also been proposed (Yang et al., 2005, 2008; Konrad and Bar-Zvi, 2008).
To validate further the candidate gene ZmASR1 and gain insight into its function, we characterized the ASR gene family of maize. We showed that ZmASR1 was one of the most highly expressed ASR genes in maize and found evidence that its expression had a large overall impact on maize vegetative productivity and yield that was maintained under water-limited stress conditions in the field. We identified 25 genes that were either transcriptionally or posttranscriptionally regulated by ZmASR1, of which seven were involved in branched-chain amino acid (BCAA) biosynthesis. Our studies further revealed a close link between 13 decreased metabolites in ZmASR1-overexpressing leaves, including two BCAAs, and support the idea that ZmASR1-dependent regulation of BCAA biosynthetic genes might contribute to yield improvement.
RESULTS
ZmASR1 Belongs to the Maize ASR Gene Family, Which Consists of at Least Nine Members
To identify ZmASR genes, we first analyzed the public databases for the presence of putative ASR coding sequences in maize using the ZmASR1 protein sequence as the query (Riccardi et al., 1998). The search was then extended to maize genomic contigs (The Institute for Genomic Research), high-throughput genomic sequences, and the recently released maize genomic sequence (Schnable et al., 2009), which were matched with the EST sequences and used to design specific primers for DNA and cDNA amplification in maize lines used for mapping analysis and transgenesis. This led to the identification of nine ZmASR genes in the maize genome (Fig. 1; Supplemental Tables S1 and S2), including the published gene Bss1 (for bundle sheath strands-specific gene 1; Furumoto et al., 2000), renamed hereafter ZmASR4. ZmASR1 to ZmASR7-3 mapped to five of the 10 maize chromosomes, with chromosomes 2 and 10 harboring two loci each and chromosome 3 harboring three loci (Supplemental Table S1). ZmASR7-1 and ZmASR7-2 were separated by 3 kb, whereas ZmASR7-3 was separated from this cluster by approximately 100 kb. When the nucleotide sequences of the nine ZmASR cDNAs were aligned, sequence identity values of 15% to 97% were observed between the individual maize genes (Supplemental Table S3). All ZmASR genes possessed a conserved structure composed of two exons and one intron, except ZmASR7-1, ZmASR7-2, and ZmASR7-3, whose first and second exons were fused (Fig. 1A). This fusion between exon 1 and exon 2 was also found in sorghum (Sorghum bicolor) SbASR6 and SbASR7, suggesting that it took place before the divergence of those two species (Supplemental Table S4). In general, the length of exons 1 and 2 was conserved (approximately 200 bp), with the exception of exon 1 of ZmASR3, which was three times longer, like its barley (Hordeum vulgare), sorghum, and rice (Oryza sativa) counterparts (Supplemental Table S4).
Gene structure of ZmASR genes. A, Schematic drawing of the exon/intron structure of ZmASR genes. Black arrows and thin gray bars indicate exons (boldface numbers) and introns (lightface numbers), respectively, the sizes of which are in bp. B, Alignment of the deduced amino acid sequences of ZmASR genes. Arrows denote the two highly conserved regions of ASR proteins: a small N-terminal consensus of approximately 18 to 20 amino acids containing a stretch of six His residues (solid arrow) and a large C-terminal region containing two ABA/WDS signatures (dashed arrows). Amino acid positions determined to be the Zn2+-dependent DNA-binding activity domain and a sequence possibly hindering DNA binding of the SlASR1 protein (Rom et al., 2006) are marked with triangles and stars, respectively. Amino acid positions determined to be the A and B regions of the bipartite nuclear localization signal of the LLA23 protein (Wang et al., 2005) are overlined. Amino acids identical to LLA23 or SlASR1 proteins are shown in boldface.
The maize ASR loci encoded small proteins of 102 to 269 amino acid long that met the hydrophilin criteria (Garay-Arroyo et al., 2000), in that their GRAVY index was between −1.13 and −1.37 and 8% to 28% of their amino acid residues are Gly, and were predicted to be “natively unfolded” in solution, given their enrichment in disorder-promoting charged amino acids (Supplemental Table S4; Goldgur et al., 2007). They all contained the two ASR-specific conserved motifs (e.g. a small N-terminal consensus of 18–20 amino acids containing a typical stretch of six His residues and a large C-terminal region of at least 80 amino acids containing two abscisic acid [ABA]/water deficit stress [WDS] signatures; Canel et al., 1995; Padmanabhan et al., 1997), except ZmASR3 and ZmASR4, which did not possess the N-terminal His-rich domain (Fig. 1B; Supplemental Table S4). The first exon of ZmASR3 encoded an additional domain that was not conserved in other ZmASR proteins. Additionally, ZmASR7-1 presented a single nucleotide insertion in a position corresponding to amino acid 57 within the second ABA/WDS signature, leading to a frameshift and consequently to a unique and maize-specific deduced amino acid sequence at the C terminus (Fig. 1B). With the exception of ZmASR4 and ZmASR5, ZmASR proteins contained a sequence motif at the C terminus (Fig. 1B; Supplemental Table S4), presenting highest identity (71%–86%) with the Zn2+-dependent DNA-binding activity domain of the SlASR1 protein (Rom et al., 2006). Close examination of the C-terminal extremity of ZmASR1, ZmASR2, and ZmASR3 revealed the presence of a region structurally similar to the bipartite nuclear localization signal of LLA23 (Fig. 1B; Supplemental Table S4; Wang et al., 2005), while five other ZmASR proteins presented only the first motif (A in Fig. 1B) but lacked the second motif (B in Fig. 1B).
To assess the diversity of the ASR gene family beyond maize, a phylogenetic tree was constructed that included the nine maize genes and 91 related genes from representatives of different groups of Spermatophyta found in the data banks (Supplemental Fig. S1). In the case of rice, we adopted the nomenclature published during the time course of this study (Philippe et al., 2010). Maximum likelihood and Bayesian reconstruction methods were congruent and confirmed the previously established division of plant ASR proteins into three distinct clades corresponding to conifer, eudicot, and monocot sequences (Supplemental Fig. S1; Frankel et al., 2007). By adding all available maize sequences, we showed that within monocots, the Commelinid sequences form a monophyletic group, with Poaceae sequences falling into two distinct clades, the I-1/I-2/I-3 clade and the II-1/II-2 clade (Supplemental Fig. S1). Within the I-1/I-2/I-3 clade, ZmASR1, ZmASR2, ZmASR3, and ZmASR4 were related to SbASR1, SbASR2, SbASR3, and OsASR6, respectively (Supplemental Fig. S1). Within the II-1/II-2 clade, ZmASR5 and ZmASR6 were related to SbASR5 and SbASR4, respectively, while ZmASR7-1, ZmASR7-2, and ZmASR7-3 appeared as three recent paralogs that were related to SbASR7 (Supplemental Fig. S1).
Response of ZmASR Transcript Levels to Water Deficit and Water Stress-Related Treatments
To investigate the role of the different ZmASR genes in the maize plant and under water stress conditions, quantitative reverse transcription (qRT)-PCR was used to monitor their transcript levels in leaf and kernel tissues at different developmental stages, the growth of both tissues being very sensitive to water deficit, and root tissue, which partly maintains growth under water deficit (Westgate and Steudle, 1985). Transcript levels of ZmASR6, ZmASR7-2, and ZmASR7-3 were not quantified because it was not possible to design gene-specific primers of sufficient quality for qRT-PCR. The relative abundance of the other six ZmASR amplicons could be compared, as we showed that the PCR efficiency for each gene was very similar (see “Materials and Methods”).
Under well-watered conditions, all the ZmASR genes tested were transcribed in leaf, kernel, and root tissues (Fig. 2). In leaves, the expression of ZmASR1 displayed a bell-shaped expression pattern (Fig. 2A), peaking at 10 cm from the leaf insertion point (e.g. beyond the growing zone, which has a constant length of 7–8 cm in maize; Tardieu et al., 2000). The transcript levels of ZmASR3, ZmASR4, ZmASR5, and ZmASR7-1 increased along the leaf, while the transcript levels of ZmASR2 were rather stable along the segments from 2 to 10 cm and decreased in the mature part of the leaf (Fig. 2A). In unfertilized and fertilized caryopses, the transcript levels of ZmASR1, ZmASR2, and ZmASR3 displayed a bell-shaped expression pattern, which peaked at 0, 7, and 9 d after pollination (DAP), respectively, and remained rather stable after 12 DAP (Fig. 2, F and G). The transcript levels of ZmASR4 increased moderately during the filling stage, whereas expression of ZmASR5 and ZmASR7-1 peaked at 7 DAP and, after a drop at 12 DAP, remained rather stable throughout development (Fig. 2, F and G). Under water deficit conditions, statistically significant transcript level increases were observed for ZmASR1, ZmASR2, ZmASR4, and ZmASR7-1 in leaves (Fig. 2A) but only for ZmASR4 in kernels (Fig. 2G).
Transcript levels of ZmASR genes in leaf, root, and kernel tissues. A, ZmASR transcript levels in different leaf sections of B73 plants grown under well-watered (control) or water deficit (deficit) conditions. L5, Mature leaf 5; L6, growing leaf 6; M, middle of the fifth leaf blade. B to E, ZmASR transcript levels in leaves (B and C) and roots (D and E) from 10-d-old F2 plantlets supplied in continuous light during 24 h with culture medium (control) complemented with 3.7 μm ABA (B and D), 1% (w/v) Glc (C and E), or 7% (w/v) PEG 8000 (C and E). F and G, ZmASR transcript levels in developing kernels from A188 (F) or MBS847 (G) plants grown under well-watered (control) or water deficit (deficit; water interruption 7 d before pollination) conditions. The light gray, gray, and dark gray boxes indicate the early, filling, and desiccation phases of the maize kernel, respectively. qRT-PCR was performed on total RNA of the indicated tissues using gene-specific primers (Supplemental Table S2). Transcript levels were normalized against the stable endogenous ZmGRP2 gene and shown relative to ZmASR5 transcript levels in the middle of the fifth leaf blade in control conditions. Values represent means of three biological replicates ± se.
For plants, the responses to water deficit include both ABA-dependent and ABA-independent pathways (Nakashima et al., 2009). Moreover, ABA induces some ASR genes in a range of plant species, including tomato, tobacco, grape, banana, sugarcane (Saccharum officinarum), and rice (Rossi et al., 1998; Sugiharto et al., 2002; Cakir et al., 2003; Takasaki et al., 2008; Henry et al., 2011). Sugars such as Glc also exhibit interactions with ABA in controlling seedling development (Rolland et al., 2006), and grape and tomato ASR proteins may be involved in Glc metabolism (Cakir et al., 2003; Frankel et al., 2007). Additionally, polyethylene glycol (PEG)-induced water stress increases tomato, potato (Solanum tuberosum), and sugarcane ASR transcript levels (Amitai-Zeigerson et al., 1995; Huang et al., 2000; Doczi et al., 2002; Sugiharto et al., 2002). Therefore, we also studied ZmASR transcript levels in response to exogenous ABA (3.7 μm), Glc (1%), and PEG (7%) treatments. Transcript levels of ZmASR1, ZmASR2, ZmASR3, and ZmASR4 were up-regulated between 4 and 18 h in control leaves and roots, while those of ZmASR5 were slightly down-regulated after 4 h and up-regulated at 12 h in control leaves and roots, respectively (Fig. 2, B and D). Comparatively, transcript levels of ZmASR7-1 were up-regulated at 12 h and stable in control leaves and roots (Fig. 2, B and D). These expression data might be due to the memory of a circadian cycle. In response to ABA treatment, we found that in leaves, transcript levels of ZmASR5 were down-regulated, while those of ZmASR1, ZmASR3, and ZmASR4 were up-regulated (Fig. 2B). In roots, up-regulation by ABA appeared a common feature to ZmASR1, ZmASR2, ZmASR4, and ZmASR7-1 (Fig. 2D). Supplying Glc barely altered ZmASR1 transcript levels in leaves (Fig. 2C). In contrast, in roots, up-regulation by Glc appeared a common feature to ZmASR1 and ZmASR5 (Fig. 2E). This Glc response was unlikely to be due to an osmotic effect, since the PEG treatment, at an equal osmotic potential, had no effect on both ZmASR1 and ZmASR5 expression (Fig. 2E). Nevertheless, supplying PEG decreased ZmASR5 transcript levels in leaves as well as ZmASR2 and ZmASR7-1 transcript levels in roots (Fig. 2, C and E).
ZmASR1 Is the Primary Detectable ZmASR Protein Responding to Water Deficit in Leaves
Expression at the transcript level is not always reflected at the protein level, so we combined separation of soluble proteins by two-dimensional gel electrophoresis (2-DE) with quantitative analysis of silver-stained gels and identification by liquid chromatography-tandem mass spectrometry (LC-MS/MS) to investigate the protein expression pattern of the ZmASR proteins in ear leaf and kernel tissues harvested at different DAPs from plants subjected to water deprivation 7 d before pollination. Our 2-DE condition allowed the analysis of five ZmASRs (ZmASR1, ZmASR2, ZmASR3, ZmASR4, and ZmASR5), the pI of the remaining ZmASRs (ZmASR6, ZmASR7-1, ZmASR7-2, and ZmASR7-3) being higher than that of the immobilized pH gradient strips used. 2-DE analysis revealed ZmASR1, ZmASR2, and ZmASR3 in ear leaves (Fig. 3, A–G) and only ZmASR1 in kernels (Fig. 3, I–K). A gel-free analysis, which has the advantage of being independent of the protein pI, gave the same results (data not shown). Interestingly, ZmASR1, encoded by one of the most highly expressed ZmASR genes, was the most abundant detected ZmASR isoform (Fig. 3H). Under well-watered conditions, its expression was rather stable during the early phase in kernels (Fig. 3H), paralleling the slight decrease observed at the transcript levels only after 7 DAP (Fig. 2G). Under water deficit conditions, its expression gradually increased in ear leaves (Fig. 3H), while it did not show any significant change in kernels (Fig. 3L). In contrast, ZmASR2 showed a relatively moderate increase in ear leaves in response to water deficit, while ZmASR3 did not show any change (Fig. 3H). These data suggest that ZmASR1 encodes the major ASR isoform for water stress responses in maize.
Expression profile of ZmASR1, ZmASR2, and ZmASR3 proteins in ear leaf and kernel tissues. 2-DE gels of ZmASR1, ZmASR2, and ZmASR3 proteins extracted from 0-, 4-, 7-, and 12-DAP ear leaf (150 μg of protein per strip) or kernel (30 μg of protein per strip) of MBS847 plants daily irrigated (control) or subjected to water deprivation (deficit) 7 d before pollination are shown. A to G, 2-DE gel images of 4-DAP ear leaf in well-watered (B, D, and F) and water deficit (A, C, E, and G) conditions with the identified ZmASR1, ZmASR2, and ZmASR3 isoforms. H, Quantification of ZmASR1, ZmASR2, and ZmASR3 isoforms in the ear leaf. I to K, 2-DE gel images of 4-DAP kernel (4K) in well-watered (J) and water deficit (I and K) conditions with the identified ZmASR1 isoform. L, Quantification of the ZmASR1 isoform in the kernel. Values represent means of three biological replicates ± se.
Ectopic Expression of ZmASR1 Maintains Kernel Yield under Water-Limited Conditions
To address ZmASR1 function, transgenic maize plants overexpressing the ZmASR1 coding sequence (ZmASR1-OE) under the control of the constitutive cassava vein mosaic virus promoter (Verdaguer et al., 1998) were generated. After selection and regeneration, followed by two backcrosses with the maize line F2, nine independent hemizygous lines that carried a unique T-DNA insertion and expressed ZmASR1 protein (data not shown) and their wild-type sister lines were selected and evaluated in a field site under well-watered and water-limited conditions (see “Materials and Methods”; Supplemental Fig. S2, A and B). qRT-PCR experiments on leaf 11 (below the ear leaf) 5 d before silk emergence showed that transgenic lines expressed ZmASR1 about 13- and 4-fold more than wild-type lines under well-watered and water-limited conditions, respectively (Supplemental Fig. S2C). An across-event analysis demonstrated that ZmASR1-OE plants exhibited a significant 3% and 4% increase in ear leaf area relative to wild-type plants under both well-watered and water-limited stress conditions, respectively, and a significant 12% increase in shoot biomass yield in well-watered conditions (Fig. 4, A and B; Supplemental Table S5). ZmASR1-OE plants also demonstrated an increase (P < 0.10) in dry leaf weight and total chlorophyll content under both well-watered and water-limited stress conditions (Supplemental Table S5). Mean yield of wild-type plants at the water deficit block was 2.5 tons ha−1, representing an approximately 43% reduction in yield compared with wild-type plants under well-watered conditions (Fig. 4C; Supplemental Table S5). Yield averages of ZmASR1-OE plants as a group were significantly greater than in wild-type plants, by 7% and 17% under well-watered and water-limited conditions, respectively (Fig. 4, C and D; Supplemental Table S5). Together, these results provided evidence that ectopic expression of ZmASR1 improved and maintained maize kernel yield under well-watered and water-limited conditions in the field.
ZmASR1-OE plants demonstrate improved vegetative and reproductive performance under water-limiting conditions. Ear leaf area (A), shoot biomass yield (B), kernel dry weight (C), and kernel number (D) values of wild-type (WT) and ZmASR1-OE plants are shown for groupings of nine individual ZmASR1 events. Values are means ± se (n = 52–60). When two samples show different letters above the bar, the difference between them is significant (normal letters, P < 0.05). When both genotype and treatment effects are significant, a, b, a+α, and b+α are indicated (see “Materials and Methods”; and Supplemental Table S5).
Transcriptomic, Proteomic, and Metabolic Adjustments in ZmASR1-OE Plants
To obtain additional cues with regard to ZmASR1 function, we used the 46K array constructed by the Maize Oligonucleotide Array Project (http://www.maizearray.org) followed by qRT-PCR experiments, 2-DE, and gas chromatography coupled to time-of-flight MS analysis to contrast the transcriptome, the proteome, and the metabolome, respectively, of ZmASR1-OE leaves with that of wild-type leaves 5 d before silk emergence under well-watered and water deficit conditions.
A first gene list of 16 differentially expressed genes was established from the transcriptome comparison between ZmASR1-OE and wild-type leaves under both treatment conditions based on P < 0.01 by two-way ANOVA and the criterion of a 1.5-fold change in transcript abundance. To confirm the differential expression, qRT-PCR experiments were performed based on the same samples that had been used for the initial microarray analysis. Because it was not possible to design gene-specific primers of sufficient quality for qRT-PCR for five target genes, the list was shortened to 11 candidates. Only three of the 11 potential target genes were significantly affected by ZmASR1-OE with opposite effects: (1) up-regulation of MZ00000739 and MZ00024643 by ZmASR1-OE under well-watered conditions; (2) down-regulation of MZ00022082 by ZmASR1-OE under both well-watered and water deficit conditions (Fig. 5; Supplemental Table S2). MZ00000739 (ZmSPL14) likely encoded a nuclear SQUAMOSA promoter-binding protein (SBP) box transcription factor of the plant-specific subfamily IIa, most closely related to AtSPL14 from Arabidopsis (Stone et al., 2005; Guo et al., 2008). MZ00022082 encoded a protein with unknown function carrying a WD40 repeat-like-containing domain and was accordingly hereafter named ZmWD40.1. Finally, MZ00024643 shared highest identity with At5g14680 (75%) and At3g01520 (71%) sequences encoding Arabidopsis vacuolar USPA (for universal stress protein A of E. coli) domain proteins (Kerk et al., 2003; Endler et al., 2006) and was accordingly hereafter named ZmUSPA.1.
Transcript levels of confirmed ZmASR1 target genes identified by microarray. qRT-PCR was performed on total RNA prepared from the 11th leaves of wild-type (WT) and ZmASR1-OE plants that had been used for the initial microarray analysis using gene-specific primers (Supplemental Table S2). Transcript levels were normalized against the stable endogenous ZmGRP2 gene and shown relative to ZmBCAT2 transcript levels in the 11th leaf in well-watered (control) conditions. Values represent means of biological duplicates ± se. When two samples show different letters above the bar, the difference between them is significant (normal letters, P < 0.05).
The proteome comparison revealed 22 protein spots that were consistently and significantly more abundant (seven spots) or less abundant (15 spots) in ZmASR1-OE leaves than in wild-type leaves in one or both treatment conditions (Table I). Subsequent analyses by LC-MS/MS led to the identification of 16 of them and showed that they were encoded by 15 distinct genes (Table I; Supplemental Table S6). These comprised five up-regulated proteins, annotated as class II aspartyl-tRNA synthetase, β-d-glucosidase, spermidine synthase 1, glucan endo-1,3-β-glucosidase 5, and NAD-dependent epimerase/dehydratase, and 11 down-regulated proteins, annotated as trigger factor-like protein, adenylosuccinate synthetase (two spots encoded by the same gene), plastid transcriptionally active 16, 3-isopropylmalate dehydrogenase (IPMDH; two isoforms, ZmIPMDH1 and ZmIPMDH2), pyruvate dehydrogenase subunit E1β, sedoheptulose-1,7-bisphosphatase, PSII stability/assembly factor, C2-domain containing protein, and thiazole biosynthetic enzyme 1-1. Interestingly, none of the corresponding 15 genes generated expression data with significant P value thresholds in the transcriptomic analysis (data not shown).
Proteins that showed significant changes in ZmASR1-OE leaves compared with wild-type leaves
Nontargeted metabolite profiling identified 84 unique compounds (Supplemental Table S7). Among them, 20 were consistently and significantly affected by ZmASR1-OE with opposite effects: increase in lactate and urea in water deficit conditions and decrease in 18 other metabolites, including 10 amino acids (Ile, Leu, Val, Phe, Trp, Asn, Gln, Pro, Ala, and Gly), two of their derivatives or potential precursors (benzoate and citramalate), and three sugars (galactinol, Glc, and Suc), in both well-watered and water deficit conditions (Table II).
Metabolites that showed significant changes in ZmASR1-OE leaves compared with wild-type leaves
ZmASR1-OE Influences the Expression of Additional BCAA- and Pro-Related Target Genes and a Specific Combination of Metabolites
Since the decrease in the three BCAAs Ile, Leu, and Val was in reasonable agreement with that of the ZmIPMDH1 and ZmIPMDH2 isoforms involved in Leu biosynthesis (Fig. 6), we hypothesized that ZmASR1 could regulate the transcription of key genes involved in drought-responsive pathways that were modified by ZmASR1-OE. Therefore, we identified those genes in the maize genome based on sequence homology to Arabidopsis genes and examined the expression level of 32 of them in wild-type and ZmASR1-OE leaves by qRT-PCR experiments (Fig. 6; Supplemental Fig. S3; Supplemental Table S2; The Arabidopsis Information Resource and AraCyc; Fang et al., 1992; Taji et al., 2002; Martin et al., 2006; Less and Galili, 2008; Urano et al., 2009). We did not find evidence for transcriptional regulation of ZmIPMDH1 and ZmIPMDH2/3 genes, in agreement with the transcriptomic analysis (Supplemental Fig. S3). However, three other BCAA-related transcripts, one for acetohydroxyacid synthase (ZmAHAS1) and two for ketolacid reductoisomerase (ZmKAR1 and ZmKARI2), showed significant increases in ZmASR1-OE leaves in comparison with wild-type leaves under well-watered conditions (Fig. 6). Furthermore, one transcript for branched-chain aminotransferase (ZmBCAT4; BCAA biosynthesis) and two transcripts for Δ1-pyrroline-5-carboxylate synthetase (ZmP5CS2 and ZmP5CS3; Pro biosynthesis) were up-regulated under both treatment conditions, while another branched-chain aminotransferase paralog (ZmBCAT2) was down-regulated in ZmASR1-OE leaves under both treatment conditions (P < 0.10; Fig. 6; Supplemental Fig. S3). These seven biosynthetic genes were induced in response to water deficit, but only one, ZmBCAT2, had an expression pattern consistent with the metabolic changes in ZmASR1-OE leaves compared with wild-type leaves (Fig. 6).
ZmASR1-OE influences BCAA-related gene expression and BCAA composition. A, The BCAA biosynthetic pathway. Substrates and products are in lightface and boldface text, respectively. Enzymes are in gray, and the number of isoforms identified in maize is indicated in parentheses. The enzymes in boldface are the rate-limiting enzymes during the water deficit response. Thick arrows indicate common steps in BCAA biosynthesis. Enzymes in boxes are ZmASR1 targets. B to D, Relative transcript (B), protein (C), and metabolite (D) levels in the 11th leaves of wild-type (WT) and ZmASR1-OE plants according to Figure 5, Table I, and Table II. Values represent means of biological duplicates ± se. When two samples show different letters above the bar, the difference between them is significant (roman letters, P < 0.05; italic letters, P < 0.10). When both genotype and treatment effects are significant, a, b, a+α, and b+α are indicated (see “Materials and Methods”). AHAS, Acetohydroxyacid synthase; DHAD, dihydroxyacid dehydratase; DW, dry weight; IPMI, isopropylmalate isomerase; IPMS, isopropylmalate synthase; TD, Thr deaminase.
To provide further insight into the impact of ZmASR1-OE on the transcriptome, proteome, and metabolome of the maize leaf, a within-group correlation matrix analysis was performed. Among 10,153 pairwise correlations computed between 36 transcripts (ZmASR1, ZmSPL14, ZmWD40.1, ZmUSPA.1, and the 32 genes related to metabolic pathways modified by ZmASR1-OE), 22 protein spots, the total chlorophyll content, and 84 metabolites, we obtained 765, 170, and 23 individual significant correlations at P < 0.05, P < 0.01, and P < 0.001, respectively (Supplemental Tables S8 and S9). These numbers were nearly twice that expected under the null hypothesis (Supplemental Table S9). Furthermore, the number of positive and negative individual significant correlations did not follow a symmetric binomial distribution (Supplemental Table S9). Similar differences were observed for the transcript-transcript, metabolite-metabolite, and transcript-metabolite submatrices, suggesting that the individual significant correlations obtained for these submatrices were most likely valuable (Supplemental Table S9). The analysis of the metabolite-metabolite submatrix revealed that 13 of the 18 metabolites that were less abundant in ZmASR1-OE leaves than in wild-type leaves showed high positive correlation to each other (Table III). Interestingly, six of them were previously reported to be negatively correlated to biomass in Arabidopsis (Table III; Meyer et al., 2007; Sulpice et al., 2009). Furthermore, six others showed high positive correlation with metabolites negatively correlated to biomass in Arabidopsis (Table III; Meyer et al., 2007; Sulpice et al., 2009). The analysis of the transcript-metabolite submatrix showed that ZmASR1 target transcripts involved in BCAA biosynthesis were positively correlated with each other and/or BCAAs and revealed a close link between six ZmASR1 target transcripts (ZmBCAT4, ZmKARI2, ZmP5CS2, ZmP5CS3, ZmUSPA.1, and ZmWD40.1) and the specific combination of 13 decreased metabolites (Supplemental Tables S10 and S11).
Pairwise correlations of ZmASR1 target metabolites against decreased and/or biomass-related metabolites
DISCUSSION
Expression of ASR Genes in Response to Water Deficit in Maize and Other Poaceae Species
Here, we identified nine ASR genes in the maize genome, making the ZmASR gene family the biggest ASR gene subfamily identified to date. Phylogenetic analyses indicated several duplication events after the divergence of Liliales from Commelinids, giving rise to five paralogous clades of ASR genes in the Poaceae family. Anchoring of the ASR genes on the maize physical map shows that all Poaceae sequences can be assigned to one of the Poaceae ancestor protochromosomes previously defined by Salse et al. (2008). Thus, sequences from the subclades I-1, I-2, I-3, II-1, and II-2 trace back to chromosomes A11/A12, A4, A1, A2/A4, and A1, respectively. Consequently, the data support the idea that the Poaceae intermediate ancestor with 12 protochromosomes (Murat et al., 2010) already had seven ASR genes. The sequence of events raising this small gene family in grasses might have taken place during the diploidization process following the Poaceae-specific whole-genome duplication (Murat et al., 2010). A larger sample of sequence data from monocots other than Poaceae might help to further elucidate this point.
Our gene-specific qRT-PCR expression data combined with 2-DE analysis and the recent work of Philippe et al. (2010) in rice strongly suggest that members of the Poaceae subclade I-1, including ZmASR1 and OsASR5, are the most prevalent ASR proteins in major plant tissues. It is noteworthy that recently released large-scale quantitative proteomics studies only revealed ZmASR1 and ZmASR2 in maize leaves (Friso et al., 2010; Majeran et al., 2010). In C4 Poaceae species, such as maize, a key component of mature leaves is the partitioning of photosynthetic processes between the bundle sheath and the mesophyll. Among the 25,800 transcripts (about 80% of the predicted maize transcriptome) that were revealed along the maize leaf gradient, Li et al. (2010) identified only five (ZmASR1, ZmASR2, ZmASR3, ZmASR4, and ZmASR5) of the nine ZmASR genes characterized in this study. It is striking that they all belong to the suite of genes (21%) that were determined to be differentially expressed between bundle sheath and mesophyll cells (Li et al., 2010). Thus, ZmASR4 showed enriched expression in bundle sheath cells, in agreement with a previous work that used differential screening analysis (Furumoto et al., 2000), while ZmASR1, ZmASR2, ZmASR3, and ZmASR5 showed enriched expression in mesophyll cells. Such a pattern of expression suggests that the ZmASR gene family shows subfunctionalization. Interestingly, almost all Poaceae sequences in subclades I-1, I-2, and II-2, including ZmASR1, ZmASR2, ZmASR3, ZmASR7-2, and ZmASR7-3, share a putative nuclear localization signal and substantial sequence identity with the Zn2+-dependent DNA-binding activity domain of LLA23 and SlASR1 proteins, raising the interesting possibility that the ASR proteins belonging to these subclades act as regulators modulating gene expression. Unlike the other members of the maize ASR gene family, ZmASR7-1, ZmASR7-2, and ZmASR7-3 do not contain any introns. Because ZmASR7-1 was expressed at a low level, whereas the expression of ZmASR7-2 and ZmASR7-3 was not determined by qRT-PCR, conventional RT-PCR was used to monitor the transcript levels of ZmASR7-2 and ZmASR7-3 and to rule out the possibility that these genes might be pseudogenes. Both transcripts were detected in leaf and root tissues, although at very low levels, regardless of the water stress-related conditions (Supplemental Fig. S4).
Additionally, our transcript and protein expression data confirmed and extended findings of up-regulation by water deficit of ZmASR1 in maize leaves (Riccardi et al., 1998, 2004). Indeed, we showed up-regulation of ZmASR1, ZmASR2, ZmASR4, and ZmASR7-1 by water deficit in leaves, although direct treatment with ABA in leaves induced elevations of ZmASR1, ZmASR3, and ZmASR4. Up-regulation by water deficit and ABA treatment was previously shown for the sugarcane and rice ZmASR1 orthologs SoDIP22 and OsASR5 (Sugiharto et al., 2002; Rabbani et al., 2003; Takasaki et al., 2008). These data indicate that members of the Poaceae subclade I-1, except ZmASR2, may play important roles in ABA-dependent pathways in the water deficit response in Poaceae leaves. In contrast, ZmASR2 and ZmASR7-1 are regulated by ABA-independent pathways under water deficit in mature leaves, as reported previously for potato DS2 genes (Silhavy et al., 1995; Doczi et al., 2002). It is worth noting that maize and rice (Philippe et al., 2010) ASR genes belonging to Poaceae subclades I-3 (ZmASR4 and OsASR6) and II-2 (ZmASR7-1 and OsASR1) were all up-regulated by water deficit in mature leaves. In contrast, maize and rice ASR genes belonging to Poaceae subclades I-2 (ZmASR3 and OsASR4) and II-1 (ZmASR5 and OsASR3) were not significantly altered and down-regulated by water deficit in mature leaves, respectively. Therefore, these data suggest that the ASR expression in response to water deficit represents an evolutionarily conserved regulatory mechanism in the Poaceae and several unrelated eudicot plants. A better understanding of the physiological roles of the maize ASR genes will require more exhaustive definition of their expression profiles in response to water deficit at the cellular level. It would also be important to determine whether these genes are affected by water deficit at the posttranslational level.
ZmASR1 Maintained Kernel Yield under Water Deficit Conditions through a Set of Low-Mr Metabolites Related to Growth Rate
ASR overexpression studies in plants are exemplified by works carried out on Solanoideae (Kalifa et al., 2004b; Frankel et al., 2007) and Arabidopsis, although the latter does not include any ASR homologous genes (Yang et al., 2005; Liu et al., 2010). Our study demonstrated that ZmASR1-OE improved the kernel yield of transgenic maize in the field under well-watered and water-limited conditions. The positive yield impact of ZmASR1 was both on kernel number and weight. These results were unexpected based on our previous PQL study, which did not reveal any QTL for maize yield and its component in the vicinity of the ZmASR1 chromosomal region (de Vienne et al., 1999; Jeanneau et al., 2002). Nevertheless, it may be pointed out that recently identified QTLs for ear length under water-stressed conditions (Lu et al., 2006) and female flowering time (Marino et al., 2009) are centered on the ZmASR1 locus. The finding that ZmASR1-OE plants demonstrated significant improvements in ear leaf area and shoot biomass yield, as well as increases in leaf dry weight and total chlorophyll content, implies that these improvements in vegetative productivity translated into improvements in reproductive performance and kernel yield. Enhanced vegetative tolerance under well-watered conditions is surprising, as we would have expected no significant differences, as described previously under mild water deficit conditions (Jeanneau et al., 2002). Nevertheless, transgenic maize plants overexpressing ZmNF-YB2 were darker green and flowered 1 to 3 d earlier than the control plants under conditions of ample water supply (Nelson et al., 2007). Additionally, enhanced vegetative tolerance and improved kernel yield under water deficit conditions were produced by the overexpression of ZmNF-YB2, betA, or Csp genes in maize (Quan et al., 2004; Nelson et al., 2007; Castiglioni et al., 2008). Continued evaluation of ZmASR1-OE across different hybrid genetic backgrounds adapted for the marketplace will be important ongoing work. Furthermore, testing the ability of other water deficit-inducible ZmASR genes for their ability to confer water deficit tolerance will further define the structure-function relationships and improve our understanding of the mode of action of the ZmASR proteins.
A noteworthy observation from our data is that 18 low-Mr metabolites decrease in mature ZmASR1-OE leaf. Of high relevance are 12 standard or related amino acids and three sugars: (1) Ala, which is synthesized from pyruvate; (2 and 3) Asn and Gln, which are the central regulators of carbon/nitrogen metabolism; (4) Gly, whose metabolism occurs by several different pathways, including Thr catabolism; (5–7) Ile (Thr-derived pathway) and the two other BCAAs, Leu and Val (pyruvate-derived pathway), that show the highest hydrophobicity among the proteinogenic amino acids and are accordingly the major constituents of transmembrane regions of membrane proteins; (8 and 9) Phe and Trp, whose synthesis initiates from chorismate and that serve in plants as precursors for a wide range of secondary metabolites having multiple functions; (10) Pro, which is synthesized from either Glu or Orn and has multiple functions in stress adaptation, recovery, and signaling; (11) benzoate, which derives from Phe and is a precursor to several important benzenoid compounds, including the defense signaling compound salicylic acid; (12) citramalate, which may serve as a precursor for Ile biosynthesis in plants as in microorganisms (de Kraker et al., 2007; Joshi et al., 2010; Less et al., 2010; Szabados and Savouré, 2010); (13) Suc, which is the major end product of photosynthesis and the major transport form of carbon from source and sink; (14) its cleavage product Glc, which enhances the expression of genes implied in carbohydrate metabolism (Koch, 1996; Trouverie et al., 2004); (15) galactinol (raffinose family oligosaccharide [RFO]), which plays a novel role in the protection of cellular metabolism, in particular, the photosynthesis of chloroplast, from oxidative damage caused by several types of abiotic stresses (Nishizawa et al., 2008). Our finding that ZmASR1-OE triggers an adjustment in the Glc level is in good agreement with the observation that a QTL for hexose content colocates with the ZmASR1 locus (Pelleschi et al., 2006) and a previous report in potato in which SlASR1 was overexpressed (Frankel et al., 2007). Additionally, the lower levels of Pro in ZmASR1-OE plants is in close agreement with previous reports in tobacco and potato plants in which SlASR1 was overexpressed and down-regulated, respectively (Kalifa et al., 2004b; Frankel et al., 2007). Lower levels of Phe, Val, and Suc had also been described for potato plants overexpressing SlASR1 (Frankel et al., 2007), although the remaining other decreased metabolites had never been identified before. Interestingly, among the 18 decreased metabolites, 13 were positively correlated to each other. Previous reports documented that biomass can be predicted by a specific metabolite status (Meyer et al., 2007; Sulpice et al., 2009). It is intriguing that among the 13 positively correlated metabolites, six (Ala, Phe, Val, Suc, benzoate, and citramalate) were previously shown to be negatively correlated to biomass (Meyer et al., 2007; Sulpice et al., 2009). It is also noteworthy that six others (transaconitate, transcaffeoylquinate, galactinol, Glc, monomethylphosphate, and Pro) showed positive correlation to metabolites negatively correlated to biomass (Meyer et al., 2007; Sulpice et al., 2009). Additionally, it was found that Pro levels accumulated under water deficit conditions were negatively correlated with maize yield (Udomprasert et al., 1999). Since ZmASR1-OE plants demonstrated improvement and maintenance in shoot biomass yield under well-watered and water-limited conditions, respectively, as well as decreases in this specific set of 13 metabolites under both treatment conditions, we hypothesize that ZmASR1 might negatively impact the level of these 13 metabolites, which would contribute to biomass yield and subsequently result in improvements and maintenance of kernel yield under well-watered and water deficit conditions, respectively.
Influence of ZmASR1 on the Expression of Genes Involved in BCAA Biosynthesis
What is the relationship between ZmASR1 and the level of these 13 metabolites that may impact on biomass and kernel yield? Our global comparative study led to the identification of 25 genes that were either transcriptionally (10 genes) or posttranscriptionally (16 proteins encoded by 15 distinct genes) up-regulated (eight transcripts and five proteins) or down-regulated (two transcripts and 11 proteins) by ZmASR1-OE in mature leaves. Three-quarters of them encoded proteins normally localized in plastids or accumulating in these organelles under osmotic stress and involved in primary metabolic pathways. Importantly, the microarray analysis and the additional qRT-PCR experiments revealed different target genes from the proteomic analysis. Therefore, the two approaches were complementary. It is quite evident that some genes with transcript changes in ZmASR1-OE leaves compared with wild-type leaves might not have been revealed by the proteomic analysis because of their low abundance. Nevertheless, none of the 16 identified proteins that showed substantial changes in abundance in ZmASR1-OE leaves compared with wild-type leaves were transcriptionally regulated in ZmASR1-OE leaves. This is not necessarily surprising, since it has been shown recently that posttranscriptional regulation can affect the protein abundance without affecting transcript levels (Böhmer and Schroeder, 2011). These data, therefore, strengthen the idea that ZmASR1, similar to other ASR proteins, may act both as a transcriptional regulator and a chaperone-like protein protector (Saumonneau et al., 2008; Urtasun et al., 2010). A synergistic relationship was found between the osmolyte Gly betaine and the SlASR1 protein, suggesting a combined mechanism of action of these two proteins in response to abiotic stresses (Konrad and Bar-Zvi, 2008). However, the enhanced Gly betaine accumulation in transgenic betA-overexpressing maize lines led to greater increases in total soluble sugars and free amino acids (Quan et al., 2004), in contrast to the metabolite effects we report in ZmASR1-OE plants. Thus, the mechanism of action of ZmASR1 is likely to be different from that of Gly betaine. Just recently, a wheat (Triticum aestivum) group 2 LEA protein named DHN-5, which is closely related to the maize LEA protein RAB17 and involved in salt and osmotic tolerance (Brini et al., 2007), has been shown to enhance the thermostability and activity of β-d-glucosidase (Brini et al., 2010), a target protein of ZmASR1. It is striking that β-d-glucosidase and 12 other ZmASR1 target genes belong to the list of established and potential thioredoxin (Trx) targets (He et al., 2009; Montrichard et al., 2009). Closer analysis shows that seven additional ZmASR1 target genes may contain Cys residues involved in disulfide bridges, making a total of 20 ZmASR1 target genes potentially linked to Trx. Thus, our results open the way to the very speculative hypothesis that ZmASR1 exhibits a chaperone-like activity on selectively redox-regulated proteins.
Among the established and potential Trx targets with expression changes in ZmASR1-OE plants, we found seven genes involved in BCAA biosynthesis. A unique feature of BCAA biosynthesis, which seems to occur exclusively in plastids, is that Val and Ile are formed in two parallel pathways using the same enzymes, namely AHAS, KARI, dihydroxyacid dehydratase, and BCAT, while Leu formation branches off from 2-oxoisovalerate, the last intermediates of the Val biosynthetic pathway, to follow a three-step chain elongation catalyzed by isopropylmalate synthase, isopropylmalate isomerase, and IPMDH, which ends with a transamination step catalyzed by a BCAT (Binder et al., 2007; Joshi et al., 2010). Decreases in ZmBCAT2 transcripts and ZmIPMDH1 and ZmIPMDH2 proteins in ZmASR1-OE leaves in comparison with wild-type leaves fit well with the finding that ZmASR1-OE triggers significant decreases in BCAAs. Furthermore, transcript levels of ZmAHAS1, ZmKARI1, ZmKARI2, and ZmBCAT4 were correlated with each other and/or the levels of the specific set of 13 decreased metabolites. Importantly, ZmKARI1 and ZmKARI2 showed enriched expression in mesophyll cells in comparison with bundle sheath cells (Majeran et al., 2005; Friso et al., 2010), similar to that was observed for ZmASR1 (Li et al., 2010). There are contradictory reports in the literature on the distribution of ZmIPMDH1 and ZmIPMDH2 across mesophyll and bundle sheath chloroplasts, with a 2-fold enrichment in either mesophyll (Majeran et al., 2005) or bundle sheath (Friso et al., 2010) chloroplasts. Distribution data for ZmAHAS1, ZmBCAT2, and ZmBCAT4 are as yet unavailable, although it is noteworthy that ZmBCAT1 and ZmBCAT3 showed enriched expression in mesophyll cells (Friso et al., 2010). Additionally, it may be noted that Arabidopsis Trx-m1 is the most effective Trx on AtIPMDH1 in comparison with Arabidopsis Trx-f1 or bacterial Trx (He et al., 2009). Its maize ortholog (Trx-m4) also showed enriched expression in mesophyll cells (Majeran et al., 2005; Friso et al., 2010). Taken together, our data demonstrate the causal relationship between the overexpression of ZmASR1 and genes involved in BCAA biosynthesis and support the idea that the improvement in biomass associated with ZmASR1-OE involves the coordinated transcriptional regulation of BCAA biosynthetic genes. The resolution of the other ZmASR1 target genes is severely hampered by the fact that only limited functional information is available on them. In this respect, a noteworthy observation is that Cañas et al. (2010) have just recently shown that increasing the expression of the enzyme P5CS may have a beneficial effect on yield. Additional experiments will be required to reveal the role of the other ZmASR1 gene targets in water deficit tolerance.
MATERIALS AND METHODS
Plant Materials and Cultures
Stress Treatments at the Vegetative Stage
To carry out the water deficit experiment, maize (Zea mays) B73 plants grown in a greenhouse (photoperiod of 16 h supplemented with artificial lighting, 25°C/18°C [day/night], and 60%/55% [day/night] relative humidity) with appropriate daily watering (50% soil water content) until the fifth leaf reached a 35-cm length (about 21 d on average) were deprived of water (WDS test) or daily irrigated (50% soil water content; control condition). The fifth and sixth leaves were harvested when growth of the sixth leaf stopped (0.2 cm or less; WDS test) or after 3 d of constant leaf elongation rate (control condition), 5 to 7 h after the onset of the light period. Four-centimeter-long fragments of the sixth leaves and a 10-cm-long fragment from the middle of the fifth leaf blade were immediately frozen in liquid nitrogen. For ABA, Glc, and PEG treatments, we used leaves and roots described in detail by Trouverie et al. (2004). All samples were ground in a 2-mL microcentrifuge tube containing both 5- and 7-mm steel beads for two to three series of 20 s using a TissueLyser II at 20 Hz (Qiagen) and stored at −80°C until analysis.
Stress Treatment at the Reproductive Stage
Tissues used to establish expression profiles during maize ovule and kernel development (−5 to 70 DAP) were described in detail by Massonneau et al. (2005). To carry out the water deficit experiment, maize MBS847 plants were grown in the greenhouse as described by Qin et al. (2004). Seven days before pollination, plants were deprived of water or irrigated daily. Sampling took place at 0, 4, 7, and 12 DAP, 5 to 7 h after the onset of the light period. The ears and 10-cm-long fragments from the bottom thirds of the corresponding axillary leaves were harvested from plants subjected to each treatment. Mature ovules and kernels were collected from the central part of the ears, separated from the adhering perianth and the pedicel, and stored in liquid nitrogen. All plant materials were ground in a mortar in liquid nitrogen and stored at −80°C until analysis.
Stress Treatment and Yield Trial in the Field
To evaluate the ZmASR1-OE plants in the field under limited water treatment, nine independently integrated ZmASR1-OE events were grown at Magneraud, France, in May 2005 using two randomized block designs (one to four rows per event; 45 plants per row in each block, transgene-positive and transgene-negative plants in each row identified by PCR during the late vegetative stage of development). One block was maintained in a well-watered condition and the other one placed under limited water supply using mobile devices preventing watering by rainfall (Supplemental Fig. S1A) and three sets of two Watermark probes dispatched at the top third, middle third, and bottom third of each block. For each location, there were two probes placed vertically in the soil, 30 and 60 cm below the surface, and connected to a data logger, allowing daily collection of soil water potential data. These data and growth stage of plants were integrated into the “grapher Watermark” software designed for a specific quality of soil, allowing the monitoring of the trial for signs of water deficit throughout the season. Consequently, the limiting-water condition was achieved at the late vegetative stage of development (e.g. from 19 d before silk emergence; Supplemental Fig. S1, B and C). Among the nine events, the three with strongest expression of the transgene (two to three rows per event in each block) were selected for further analysis. Two-centimeter-long fragments in the bottom third of the 11th leaves below the ear leaf were collected 57 d after sowing (5 d before silk emergence), 5 to 7 h after the onset of the light period, and immediately frozen in liquid nitrogen. They were then ground in a 2-mL microcentrifuge tube containing both 5- and 7-mm steel beads for two to four series of 20 s, using a TissueLyser II at 20 Hz (Qiagen). Biological duplicates were formed by dividing each row in two groups and pooling the three selected events, without changing the proportion of each row: samples included 46 plants in total, with 14 to 18 plants per event and two to three rows per event. All samples were stored at −80°C until analysis.
Plant Transformation
The plasmid used for the production of ZmASR1-OE plants contains the backbone of vector pSB12 (Komari et al., 1996), a Basta resistance cassette (rice Actin promoter and intron, Bar and Nos terminators) next to the right border, and the ZmASR1 cDNA of maize MBS847 (GenBank accession no. AX297905), under the control of the cassava vein mosaic virus promoter, next to the left border. Agrobacterium tumefaciens-mediated transformation of maize inbred line A188 was based on a published protocol (Ishida et al., 1996). For each transformation event, the number of T-DNA insertions was evaluated by Southern blot, the integrity of the transgene was verified by PCR, and the expression level of the transgene was evaluated by RT-PCR and 2-DE as described (Jeanneau et al., 2002; Depège-Fargeix et al., 2011). Seventeen independent transformation events were crossed with the maize inbred line F2, in which the ZmASR1 protein is not detected (Riccardi et al., 1998; Jeanneau et al., 2002). T1 events were then selected and backcrossed one more time with the maize inbred line F2. Among the 17 independently integrated T2 events, the nine with strongest expression of the transgene compared with the endogenous ZmASR1 gene based on 2-DE were evaluated in the field. In the T2 generation containing 50% transgenics, the wild-type plants were used as control plants.
Database Search, Gene Structure Determination, and Chromosomal Location of ZmASR Genes
To search for ZmASR1 homologs in maize, we used the ZmASR1 cDNA in BLAST analysis using EST assemblies in The Gene Index database (http://compbio.dfci.harvard.edu/tgi/plant.html), assemblies of genomic DNA fragments at The Plant Genomics Resources (http://blast.jcvi.org/tgi_maize/index.cgi), high-throughput genomic sequences from bacterial artificial chromosome clones at the National Center for Biotechnology Information (NCBI; http://www.ncbi.nlm.nih.gov/HTGS), and the recent release of the maize genome sequence (http://maizesequence.org). The members of the maize ASR gene family suggested from these combined in silico analyses were confirmed by sequencing of amplified cDNA and genomic DNA of maize B73 plus one to five additional maize inbred lines (Mo17, F2, F252, MBS847, and/or A188). They were then mapped on the intermated B73 × Mo17 (IBM) or F2 × F252 (LHRF) mapping panels using the REFMAP050110 (IBM_Gnp2004) or LHRF_Gnp2004 framework map, respectively (Falque et al., 2005), using gene-specific primers (Supplemental Table S2). Alternatively, BLAST analysis on MaizeSequence.org was used in the absence of polymorphism between the maize inbred lines that led to the IBP or LHRF mapping panels. Deduced amino acid sequences were annotated by BLAST against known ASR proteins at NCBI and query against InterPro (http://www.ebi.ac.uk/Tools/pfa/iprscan/). The ProtParam tool (Gasteigler et al., 2005; http://web.expasy.org/protparam) was also used to determine the physical and chemical properties of the deduced amino acid sequences.
Phylogenetic Analysis
Sequence alignments were managed using BioEdit version 7.0.0 (Hall, 1999) and visually refined on the basis of the amino acid sequences. Bayesian phylogenetic analysis was carried out using MrBayes version 3.1.1 (Huelsenbeck and Ronquist, 2001) with a general time-reversible model, with a proportion of invariable sites and a γ-distribution for site-specific rates and a partition according to the three codon positions on the four most conserved regions only (272-nucleotide alignment matrix). Three chains (two-heated) were run twice for 107 generations, with a sampling of the cold chain parameters every 100 generations and a burn-in of 25,000 samples. Convergence was followed with potential scale reduction factor and average sd of split frequencies. A majority rule consensus tree was built with posterior probabilities of nodes above 0.85 × 100 indicated.
qRT-PCR and RT-PCR Analyses
The total RNA isolation from maize tissues collected in the stress experiments (leaves, roots, mature ovules, and kernels) with the TRIzol Reagent (Invitrogen) and ethanol precipitation was carried out as described (Trouverie et al., 2004). Subsequent DNase treatment and DNase inactivation were carried out according to the instructions of the supplier (Ambion). Two micrograms of total RNA was reverse transcribed using random hexamers (Invitrogen), 100 units of SuperScript II (Invitrogen), and 40 units of recombinant Rnasin RNase inhibitor (Promega) in a final volume of 20 μL. For expression profiles during maize ovule and kernel development, the total RNA isolation with phenol/chloroform extraction and acetic acid/ethanol precipitation, as well as subsequent DNase treatment and RT, were carried out as described (Massonneau et al., 2005).
Primer sequences were designed using Primer Express 2.0 (Applied Biosystem). Subsequent qRT-PCR analysis was carried out as described (Capelle et al., 2010). The calibration step of the experiment checked for equivalent PCR efficiency of the different genes (to allow comparison and normalization). Standard curves (log of cDNA dilution versus cycle threshold) using serial 10-fold dilution of cDNA were built for each pair of selected primers, a 100% efficiency corresponding to a slope of −3.3 (Marino et al., 2003). Practically, only pairs of primers yielding a slope of −3.3 ± 0.1 were selected. The specificity of the amplification (checked by dissociation curve analysis, gel electrophoresis, and sequencing of the PCR product) and the use of appropriate control genes were also assessed. Normalization of the results was achieved using the Gly-rich RNA-binding protein2 (ZmGRP2) gene. After an initial test of five genes (actin, polyubiquitin, α-tubulin, H+-ATPase, and ZmGRP2) on the range of tissues and treatments we wished to compare, only ZmGRP2 was found to be stably expressed, in agreement with recent genome-wide observations in 55 tissues of maize B73 line (Sekhon et al., 2011). Conventional RT-PCR analysis was used to provide evidence that ZmASR7-2 and ZmASR7-3 genes are expressed. It was performed using the real-time PCR system and the SYBR Green PCR master mix as described (Capelle et al., 2010) to allow comparison with the qRT-PCR experiments. Products of the PCR were loaded on agarose gels and stained with ethidium bromide.
Microarray Analysis
Total RNA was extracted using the Nucleospin 8 RNA Extraction Kit (Macherey-Nagel) and in vitro amplified based on suggestions from the Maize Oligonucleotide Array Project (http://www.maizearray.org/files/cRNA_Target_Production_Using_RNA_Amplification.pdf). Hybridization was carried out as described (http://www.maizearray.org/files/Hybridization_Protocol_For_cRNA_Targets.pdf) using the 46K array constructed by the Maize Oligonucleotide Array Project. Quantile normalization of the raw data was carried out using the R version 2.4.1 software with a log10 transformation of the data initially performed. The criteria for the inclusion of a gene in the list of differentially expressed genes were a 1.5-fold increase or decrease and P < 0.01 by two-way ANOVA. The cDNA and genomic sequences corresponding to the 70-mers present on the microarray were established by BLAST of the 70-mers against the high-throughput genomic sequences database and the recent release of the maize genome sequence and regularly updated. Deduced amino acid sequences were annotated by BLAST against the Arabidopsis (Arabidopsis thaliana) genome at NCBI and screened for known conserved domains using the Center for Biological Sequence Analysis database (http://www.cbs.dtu.dk/services/).
Protein Extraction, Gel Staining, and 2-DE Analysis
Protein extraction with TCA/acetone precipitation was carried out as described (Méchin et al., 2007). Proteins were solubilized in UKS (for urea-K2CO3-sodium dodecyl sulfate) buffer (Méchin et al., 2007) and quantified using the 2-D Quant Kit (GE Healthcare Life Sciences). Isoelectrofocusing was carried out using 24-cm-long, pH 4 to 7 Immobiline Dry Strips (GE Healthcare Life Sciences) to which 50 μg (for protein quantification; 150 and 30 μg for ear leaf and kernel tissues, respectively) or 300 μg (for protein identification) of dissolved protein was added. Active rehydration and focusing were achieved at 20°C in a Protean IEF Cell (Bio-Rad) by increasing the voltage step by step from 50 to 10,000 V (13 h at 50 V, 0.5 h at 200 V, 0.5 h at 500 V, 1 h at 1,000 V, then increase to 10,000 V) and stopped when 84,000 Vh was reached. Strips were equilibrated as described (Görg et al., 1987) and sealed at the top of a 1-mm-thick two-dimensional gel (24 × 22 cm) with 1% (w/v) low-melting agarose in SDS electrophoresis buffer (25 mm Trizma base, 0.2 m Gly, and 0.1% [v/v] SDS). Separation of continuous gels (11% [w/v] acrylamide, 2.9% [w/v] piperazine diacrylamide) was performed at 14°C (20 V for 1 h, 140 V for 15 h, and 20 V for 2 h) in a Protean Plus Dodeca Cell Electrophoresis Chamber (Bio-Rad) until the bromphenol blue front reached the end of the gel. For protein quantification, a silver-staining procedure was performed as described (Méchin et al., 2003). Scanning was carried out at 300 dpi with a 16-bit grayscale pixel depth using an image scanner (Amersham Biosciences) and analyzed with Progenesis software (Nonlinear Dynamics) according to Zivy (2007). For protein identification, gels were stained in colloidal Coomassie blue according to Yan et al. (2000) and spots were excised for MS. Excised gels were restained with silver nitrate as described (Méchin et al., 2003) and compared with the nonexcised silver-stained gels.
LC-MS/MS and Protein Identification
Two-dimensional gel spot digestion and LC-MS/MS analysis were performed as described (Page et al., 2010). A database search was performed with X!Tandem (version 2010.01.01.4; http://www.thegpm.org/tandem/). Enzymatic cleavage was declared as a trypsin digestion with one possible miscleavage. Cys carboxyamidomethylation and Met oxidation were set to static and possible modifications, respectively. Precursor mass and fragment mass tolerance were 2.0 and 0.8, respectively. A refinement search was added with similar parameters, except that semitryptic peptide and possible N-terminal protein acetylation were searched. The maize genome sequence (release 4a.53; 53,764 entries), the UniProt database restricted to maize (http://www.uniprot.org; release 15.11; 43,694 entries), and a contaminant database including in particular trypsin and keratins were used. Only peptides with a E value smaller than 0.1 were reported. Identified proteins were filtered and grouped using the X!Tandem pipeline (http://pappso.inra.fr/bioinfo/xtandempipeline/) with the following criteria: (1) a minimum of two different peptides was required with an E value smaller than 0.05; (2) a protein E value (calculated as the product of unique peptide E values) smaller than 10−4 was required. In the case of identification with only two or three MS/MS spectra, similarity between the experimental and the theoretical MS/MS spectra was visually checked. To take redundancy into account, proteins with at least one peptide in common were grouped. This allowed grouping of proteins of similar function. Within each group, proteins with at least one specific peptide relative to other members of the group were reported as subgroups.
Chlorophyll Content Measurement and Metabolite Analysis
Total chlorophyll content was measured after extraction of leaf material into 80% acetone and assay at 663 and 645 nm as described (Arnon, 1949). Gas chromatography coupled to time-of-flight MS analysis was carried out as described (Noctor et al., 2007). Because automated peak integration was occasionally erroneous, integration was verified manually for each compound in all analyses.
Two-Way ANOVA and Within-Group Correlation Matrix Analysis
A log10 transformation of the data was initially performed when the se values showed that the se increased in proportion to the treatment. Either nontransformed or log10-transformed values were then used for two-way ANOVA using R software version 2.8.1 (http://cran.at.r-project.org). The two possible means models for two-way ANOVA are the additive model and the interaction model (Scheffé, 1999; Seltman, 2010). An F test for interactions was performed to determine whether the additive model could be retained. When the interaction model was needed, the Bonferroni method was applied for pairwise comparisons. The significance was placed at the 0.05 level.
In a multivariate ANOVA model, the within-group correlation matrix is calculated from the correlation matrix of residuals after regressing out the influence of the fixed effects by linear regression (Hair et al., 2010). Functions of the R software were used. We were unable to test the association between two variables independently of the others because we did not have enough data to determine partial correlation coefficients. Therefore, the null hypothesis of a pairwise correlation coefficient was evaluated by a test derived from an F test and extended to all correlation coefficients using the Bonferroni method. The significance was placed at the 0.05 level.
Supplemental Data
The following materials are available in the online version of this article.
Supplemental Figure S1. Majority rule consensus tree obtained using Bayesian interference analysis of 100 ASR sequences.
Supplemental Figure S2. Stress treatment and yield trial of ZmASR1-OE plants.
Supplemental Figure S3. Relative expression profile of selected genes encoding ZmASR1 targets and/or rate-limiting enzymes in leaves 11 of ZmASR1-OE and wild-type plants.
Supplemental Figure S4. Transcript levels of ZmASR7-1, ZmASR7-2, and ZmASR7-3 genes in leaf and root tissues.
Supplemental Table S1. ZmASR gene mapping.
Supplemental Table S2. Sequences of the primers used for mapping and qRT-PCR analyses.
Supplemental Table S3. Sequence identity between individual ZmASR genes.
Supplemental Table S4. Summary of exon lengths of various ASR genes from Poaceae and properties of their deduced amino acid sequences.
Supplemental Table S5. ZmASR1-OE plants maintain kernel yield under water-limited conditions.
Supplemental Table S6. Proteins that showed significant changes in ZmASR1-OE leaves compared with wild-type leaves, their annotated function, location, and quantification.
Supplemental Table S7. Unique metabolites identified in ZmASR1-OE and wild-type leaves.
Supplemental Table S8. Within-group correlation matrix.
Supplemental Table S9. Number of individual significant correlations for the whole matrix and its submatrices.
Supplemental Table S10. Pairwise correlations of ZmASR1 targets against BCAA-related ZmASR1 targets.
Supplemental Table S11. Pairwise correlations of ZmASR1 target transcripts against decreased and/or biomass-related metabolites.
Acknowledgments
We thank Ingrid Doridant and Carine Hécart for maize transformation and culture, Carine Remoué for help with physical mapping, the Biogemma transcriptomic team for the hybridization and analysis of the microarray, and Jean Coursol for guiding us through the statistical analyses. Caroline Lelarge-Trouverie is acknowledged for preliminary metabolomic analysis on ZmASR1-OE.
Footnotes
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: Sylvie Coursol (coursol{at}moulon.inra.fr).
↵1 This work was supported by the Génoplante projects B06 and WaterLess (grant no. ANR05–GPLA–034–05) and by a Ph.D. research fellowship of the Université Paris-Sud 11 to L.V.
↵2 Present address: Université Picardie Jules Verne, EA3900 Biologie des Plantes et Contrôle des Insectes Ravageurs, F–80039 Amiens, France.
↵3 Present address: INRA, UMR 950 Ecophysiologie Végétale, Agronomie et Nutritions N, C, S, F–14032 Caen, France.
↵4 Present address: Field Crop Division, Limagrain Europe, Site de la Garenne, F–63720 Chappes, France.
↵[OA] Open Access articles can be viewed online without a subscription.
↵[W] The online version of this article contains Web-only data.
- Received March 22, 2011.
- Accepted August 4, 2011.
- Published August 18, 2011.