Transcriptomic and proteomic analyses of a pale-green durum wheat mutant shows variations in photosystem components and metabolic deficiencies under drought stress
BMC Genomics volume 15, Article number: 125 (2014)
Leaf pigment content is an important trait involved in environmental interactions. In order to determine its impact on drought tolerance in wheat, we characterized a pale-green durum wheat mutant (Triticum turgidum L. var. durum) under contrasting water availability conditions.
The pale-green mutant was investigated by comparing pigment content and gene/protein expression profiles to wild-type plants at anthesis. Under well-watered (control) conditions the mutant had lower levels of chlorophylls and carotenoids, but higher levels of xanthophyll de-epoxidation compared to wild-type. Transcriptomic analysis under control conditions showed that defense genes (encoding e.g. pathogenesis-related proteins, peroxidases and chitinases) were upregulated in the mutant, suggesting the presence of mild oxidative stress that was compensated without altering the net rate of photosynthesis. Transcriptomic analysis under terminal water stress conditions, revealed the modulation of antioxidant enzymes, photosystem components, and enzymes representing carbohydrate metabolism and the tricarboxylic acid cycle, indicating that the mutant was exposed to greater oxidative stress than the wild-type plants, but had a limited capacity to respond. We also compared the two genotypes under irrigated and rain-fed field conditions over three years, finding that the greater oxidative stress and corresponding molecular changes in the pale-green mutant were associated to a yield reduction.
This study provides insight on the effect of pigment content in the molecular response to drought. Identified genes differentially expressed under terminal water stress may be valuable for further studies addressing drought resistance in wheat.
The leaf chlorophyll content of plants determines their capacity to absorb energy from sunlight and is therefore considered a measure of photosynthetic potential. Plants typically respond to water stress by reducing the leaf chlorophyll content, and hence their photosynthetic activity. Genotypes that tolerate water stress tend to have more leaf chlorophyll and lower canopy temperatures when stressed. In contrast, genotypes that are susceptible to water stress often have higher canopy temperatures, and when such plants are well-watered, with their stomata fully open, water evaporation causes leaf cooling. The ability of a genotype to remain green under water stress conditions may therefore enhance photoassimilation, thus increasing productivity and grain filling. In contrast to these general statements, a low level of leaf chlorophyll can limit the interception of light to avoid oxidative damage caused by excess radiation and could be considered a desirable trait under water stress conditions. The optimization of photosynthesis in crops under water stress therefore requires a productive balance between assimilation under favorable conditions and the avoidance of excess radiation under stress.
Chlorophylls and carotenoids play an important role in photoassimilation and the synthesis of both molecules is coupled to chloroplast development. Thylakoid organization is coordinated with the expression of nuclear genes, encoding pigment-binding proteins that are imported into the chloroplast to assemble the photosynthetic complexes and the expression of chloroplast encoded genes[8, 9]. This process is tightly regulated, so mutants deficient in chlorophyll content can arise not only through the disruption of chlorophyll biosynthesis but also via processes related to photosystem assembly and light-harvesting activity[10, 11]. Such chlorophyll mutants have been used to study the function of photosynthetic components as described in rice e.g.[12, 13], barley e.g.[14–17] and wheat e.g.[18–20].
Carotenoids are involved in photoprotection and in light harvesting, and different carotenoids play different roles in the photosynthetic complexes. For example, only β-carotene is present in the reaction centers whereas xanthophylls are present in the antennal light harvesting complex (LHC) where each protein has a specific carotenoid complement that ensures maximum photosystem efficiency. Lutein is involved in light harvesting and the quenching of 3Chl states and violaxanthin can be replaced by zeaxanthin if the de-epoxidation cycle is activated by photosystem saturation following the absorption of excess light. Zeaxanthin helps to prevent photo-inhibition by dissipating the excess of energy as heat, either by direct quenching or by inducing conformational changes in the LHC proteins. Photo-inhibition is an early consequence of excess light and other environmental stresses, reducing photosynthetic activity and leading to pigment photo-oxidation, the production of reactive oxygen species (ROS) and ultimately cell death. Chloroplasts and mitochondria are the primary sources of ROS, which damage cells by oxidizing nucleic acids, proteins and lipids. However, low concentrations of ROS act as potent signaling molecules that coordinate responses to abiotic stress, pests and pathogens. The balance between ROS production and scavenging by antioxidants can be disrupted by different forms of stress.
Wheat is a staple crop that plays a major role in global food security. Conventional breeding strategies have been used to increase grain yields and stress tolerance, but in order to feed a growing population it will be essential to include molecular breeding as well as conventional approaches to achieve further substantial gains[27, 28]. Gene expression profiling with microarrays has been used to characterize the response of wheat to heat and drought stress, cold, Fusarium infection and phase-transition, and to investigate the molecular basis of grain quality traits. In the context of drought stress, Krugman et al. compared resistant and susceptible wild emmer genotypes to investigate the overlapping regulatory and signaling processes overrepresented in the drought-resistant genotype and to identify candidate genes for drought tolerance. Proteomics has also been used to investigate the response of wheat to nitrogen, salinity stress[36, 37], low temperatures and drought[39, 40]. Poor correlation between transcriptomic and proteomic data has been reported in some stress experiments although such complementary information highlights the multilevel regulation of responses to complex environments[41, 42].
Most of the experiments discussed above involved genotypes with contrasting phenotypes and different genetic backgrounds, an approach that allows many genes and stress tolerance strategies to be identified, but makes it difficult to characterize specific genotype-phenotype relationships. This can be achieved in experiments that involve genotypes differing at a single locus e.g.. In order to control for genetic background effects, we compared a pale-green mutant (M) and its wild-type mother variety (WT) to determine the impact of leaf pigment content on gene and protein expression. We studied the chlorophyll and carotenoid content as well as the transcriptomic and proteomic profiles of the mutant and wild-type plants under well-irrigated (control) conditions and terminal water stress. We also tested the performance of the mutant and wild-type plants under irrigated (control) and rain-fed field conditions in order to assess the effect of the mutation on biomass accumulation and grain yield.
Wild-type and mutant leaf sections at the tillering stage were compared by transmission electron microscopy (TEM). As shown in Figure 1, lens-shaped chloroplasts were present in the wild-type mesophyll cells whereas those in the mutant plants were typically cylindrical. The thylakoid grana of the wild-type chloroplasts appeared to be connected by continuous intergranal lamellae parallel to the long axis of the chloroplast, whereas the mutant grana appeared to be disconnected and disorganized within the chloroplast stroma. Osmiophilic lipid globules (plastoglobules) were more abundant in the mutant chloroplasts (shown by arrows in Figure 1C and Figure 1D).
Physiological responses and field measures
Water stress was applied to plants in the glasshouse by reducing the soil water content at anthesis. This increased the negative value of the midday flag leaf water potential (ΨW), reduced the leaf conductance to water vapor values (gs) and reduced the mean net photosynthetic rate (An) but differences between genotypes were not statistically significant (Table 1). In the field, differences in aboveground biomass were significant only under rain-fed conditions, but the grain yield of the mutant plants was reduced by 10.3% under irrigated conditions and by 14.3% under rain-fed conditions, compared to wild-type plants in the same environment (Table 2).
The total carotenoid and chlorophyll levels in the flag leaves at anthesis under well-watered (control) conditions were significantly lower in the mutant than the wild-type, but differences in antheraxanthin and zeaxanthin levels were not statistically significant (Table 3). Under water stress conditions, only the antheraxanthin content differed significantly between wild-type and mutant leaves, whereas the total carotenoid and chlorophyll contents were similar in both genotypes. The xanthophyll de-epoxidation cycle appeared to be induced in the mutant plants under control conditions. Whereas under terminal water stress, wild-type plants increased the de-epoxidation rate, mutant plants did not. Therefore, there was no statistically significant difference between the genotypes under water stress conditions (Table 3).
Microarray quality analysis
Microarray data were collected using the Affymetrix Wheat Genome Array and normalized using the RMA algorithm. The average background was 5.33, well within the recommended levels. The percentage of "Present" calls ranged from 39.05% to 42.10% among the 61 k probe sets. The quality of the biological replicates was determined by calculating R2 values among replicates of the same sample and the values ranged between 0.988 and 0.997 with an average of 0.995. Principal component analysis highlighted the main sources of total variance. The two components explaining 59.3% and 15.8% of the variance represented the differences between the treatments (control vs. water stress) and between the genotypes (pale-green mutant vs. wild-type) as shown in Figure 2.
The mutation effect on gene expression was studied comparing the expression of the mutant to wild-type for each water treatment and indicated along the manuscript with the terms induced or repressed genes. In addition, the stress effect was evaluated comparing the gene expression of the water stress to control conditions for each genotype and indicated along the manuscript as up/downregulated genes.
Gene expression under control conditions
Transcriptional profiling at anthesis revealed differences in gene expression between the two genotypes under control conditions. We identified 23 induced and 15 repressed probe sets in the mutant (>2-fold difference; P < 0.05). The array data for selected probe sets were validated by quantitative real-time RT-PCR, showing strong agreement between the fold-change values observed with each method (r = 0.96, P < 0.0001; Additional file1: Figure S1).
MIPS FunCatDB was used to distribute the genes into functional categories, revealing that the induced probe sets in the mutant represented chitin catabolism, defense (particularly systemic acquired resistance) and cell wall metabolism (Table 4). The most strongly induced probe sets represented peroxidase 10 (POX10), pathogenesis-related (PR) proteins including beta-glucanases and chitinases, and a TaWRKY45-like transcription factor (Table 5). The repressed probe sets represented functions including systemic interaction with the environment, abscisic acid (ABA) responses, hormonal regulation, stress responses and DNA degradation (Table 4), e.g. sequences representing RAB (responsive to ABA) proteins, late embryogenesis abundant protein LEA1 and endonucleases (Table 5). Probe sets related to cell rescue, defense and virulence (functional category 32) were induced or repressed in the mutant when compared with the wild-type. However, none the subcategories was statistically significant at P ≤ 0.005.
Gene expression under water stress
We initially compared each genotype under stress and control conditions to identify both general and genotype-specific stress-modulated genes, which were defined as upregulated (induced by water stress) or downregulated (repressed by water stress), as mentioned above. We identified a total of 1817 probe sets that were modulated by water stress in at least one of the genotypes (Figure 3), including 1282 genes in the wild-type plants, 1422 genes in the mutant plants and 887 in both.
We also compared the wild-type and mutant plants under water stress conditions, to identify genes that were differently-expressed under stress. We identified 166 probe sets that were defined either as induced or repressed in the mutant plants, according to the level of expression under water stress conditions.
Genotype-independent stress-response genes
From the 887 probe sets modulated in both genotypes, 286 were upregulated and 601 were downregulated under water stress (Figure 3). The upregulated genes predominantly represented the following functional categories: stress response, carbon compound and carbohydrate metabolism, metabolite transport, amino acid metabolism, ABA response and proline biosynthesis (Additional file2: Table S1). Genes already known to be associated with water stress are shown in the Additional file3: Table S2. These include genes encoding sucrose synthase (SUS3, SUS4), glucose-6-phosphate dehydrogenase (G6PD4), branched-chain-amino-acid aminotransferase (BCAT3) and branched-chain alpha-keto acid dehydrogenase E1 α subunit (BCKDHE1α), delta 1-pyrroline-5-carboxylate synthetase 2 (P5CS2), pyrroline-5-carboxylate reductase (P5CR) and RAB18. The downregulated genes predominantly represented the following functional categories: phosphate and secondary metabolism, e.g. 1-aminocyclopropane-1-carboxylate synthase (ACS) and 12-oxophytodienoate reductase (OPR); interaction with the environment, e.g. chitinases, β-glucanases and pathogenesis-related proteins; immune response, e.g. transcription factor WRKY; and phenylpropanoid metabolism, e.g. phenylalanine ammonia-lyase. These genotype-independent stress-modulated genes correspond to the well-characterized drought stress response of wheat, confirming that the mutant plants preserve much of the stress-response signaling network that allows wheat plants to tolerate water stress[34, 43].
Genotype-dependent stress-response genes in wild-type plants
We identified 395 probe sets that were modulated solely in the wild-type plants or to a greater extent (above the 2-fold threshold) than in the mutant plants. This included 265 upregulated and 130 downregulated genes (Figure 3). The upregulated genes predominantly represented carbohydrate and pyruvate metabolism and the glyoxylate cycle, whereas the downregulated represented metabolism, interaction with the environment, protein modification and signaling pathways (Table 6). Selected genes specifically modulated by the wild-type plants under stress are listed in Table 7, including carbon fixation enzymes such as malate dehydrogenase (NADP-ME3) and phosphoenolpyruvate caboxykinase (PCK1), and glyoxylate cycle enzymes such as isocitrate lyase (ICL) and malate synthase (MLS).
Genotype-dependent stress-response genes in mutant plants
We identified 535 probe sets that were modulated solely in the mutant plants or to a greater extent (above the 2-fold threshold) than in the wild-type plants. This included 157 upregulated and 378 downregulated genes representing diverse functional categories (Figure 3; Table 8). The most significant categories in the upregulated group were associated with chloroplast transcription, interactions with the environment (cellular sensing, response to external stimuli, chemoperception and response), stress responses and metabolism, whereas the downregulated group was dominated by functional categories such as interaction with the environment (response to external stimuli), stress responses and plant defense (Table 8). The downregulated genes also included those induced in the mutant under control conditions.
The genes that were upregulated specifically in the mutant included those involved in the synthesis of carotenoid precursors, e.g. hydroxymethylglutaryl-CoA synthase (MVA1), 1-deoxy-D-xylulose-5-phosphate synthase (DXS) and phytoene synthase (PSY), and those involved in chloroplast gene expression, photosynthesis, chlorophyll binding and detoxification, e.g. glutamyl-tRNA reductase 1 (HEMA1), light harvesting complex II components such as LHB1B2 and ATPase subunit delta, and the Rubisco small subunit. We observed the downregulation of several genes representing glutamate metabolism such as glutamate decarboxylase (GAD), glutamate dehydrogenase (GDH), γ-glutamyltranspeptidase (GGT) and L-aspartate oxidase (AO); and several signaling pathway components, including calcineurin B-like protein 1 (CBL1), inositol phosphatase (SAL1) and MPK3 (Table 7).
Differentially-expressed stress-response genes under stress conditions
We compared the transcriptional profiles of wild-type and mutant plants under water stress and identified 166 differently-expressed probe sets (29 induced and 137 repressed in the mutant). The only functional category significantly represented among the 29 induced probe sets was the channel/pore transport class (Additional file4: Table S3), corresponding to the transcripts for two plasma membrane intrinsic proteins (PIP2;2, PIP2;5) that were downregulated in both genotypes but to a lesser extent in the mutant (Additional file3: Table S2). Other genes that were induced in the mutant under water stress included those encoding ascorbate peroxidase (APX4), sucrose: fructan 6-fructosyltransferase (SFT6), myo-inositol oxygenase (MIOX1) and dihydrodipicolinate reductase 3 (crr1) (Table 7 and Additional file3: Table S2).
The genes that were repressed in the mutant plants under stress were predominantly related to carbohydrate metabolism, peroxidase/oxidative stress responses, the glyoxylate cycle and chromosome subcellular localization (Additional file4: Table S3). Representative repressed genes included those encoding RAP2.8, GDH, GAD, CBL-interacting protein kinase 23 (CIPK23) and glyoxalase (Table 7). We also found genes encoding ICL, MLS and trehalose phosphate-synthase (TPS) which were not affected by stress in mutant plants but were upregulated in wild-type plants, thus appearing in the repressed transcript category (Table 7).
The effect of the mutation was also investigated by carrying out a comparative proteomic analysis between wild-type and mutant plants at anthesis under the same conditions as described above.
When wild-type and mutant plants were compared under well-watered (control) conditions, we matched 182 protein spots (data not shown) among which two showed a significant (>2-fold) difference in signal intensity between genotypes, although neither spot could be identified because of the low abundance of the corresponding protein. Even so, this experiment supported the transcriptomic and suggested that there were only minor differences between the genotypes under normal conditions.
When we compared each genotype under stress and control conditions, we matched 207 protein spots among both genotypes. From the 207, we found that 32 spots were modulated by water stress in the wild-type plants and 62 were modulated in the mutant plants, including 25 spots that were common to both genotypes. This left seven modulated spots specific to the wild-type plants and 37 specific to the mutant.
When the two genotypes were compared under water stress, we identified 28 spots that were more abundant in the mutant and 17 that were less abundant in the mutant. We selected 19 of the induced spots and all 17 of the repressed spots based on the likelihood of extracting useful amounts of protein, and we identified 13 proteins from the induced set and 8 from the repressed set (Additional file5: Figure S2). The proteins were identified by peptide mass fingerprinting (Additional file6: Table S4) and MS/MS analysis (Additional file7: Table S5). The induced proteins predominantly represented carbohydrate and amino acid metabolism, photosynthesis, detoxification and phytoalexin biosynthesis, whereas the repressed proteins predominantly represented photosynthesis (Table 9). Two proteins that were specifically upregulated in the mutant were the oxygen-evolving enhancer protein 2 corresponding to the PSBP subunit of photosystem II (A.n. Q00434), and catalase-1 (CATA1; A.n. Q00434). Two proteins that were specifically downregulated in the mutant were the Rubisco small subunit chain PW9 (A.n. P26667) and the photosystem I reaction center subunit IV (A.n. P13194).
Chloroplast ultrastructure and pigment content
Chlorotic mutants with affected chloroplast ultrastructure have been described in several reports[9, 14, 44–46]. The chloroplasts in these plants often contain fewer grana structures than wild-type plants, and thinner thylakoid membranes that are not stacked in the normal manner. In contrast, the chloroplasts in our pale-green mutant contained abundant grana with normally-stacked thylakoids, but their distribution within the stroma was disorganized, potentially reflecting a lack of membrane continuity or the degradation of inter-grana lamella membranes. The stacking adhesion of the thylakoids in grana structures is regulated by the phosphorylation of photosystem II proteins, which also ensures photosynthetic efficiency. We found no significant difference in the photosynthetic rate of the wild-type and mutant and no quantitative differences in transcripts or proteins representing the photosynthetic apparatus at anthesis under control conditions.
The mutant chloroplasts also contained more stromal plastoglobules than wild-type chloroplasts. These lipoprotein bodies are structurally continuous with the thylakoid membranes and facilitate content exchange with the thylakoid lumen. They may also be involved in thylakoid development and disassembly, increasing in abundance in response to oxidative stress and delivering the antioxidant molecules stored on them to the thylakoid lumen, where they scavenge ROS[49, 50]. The structure of the inter-grana lamella in the mutant chloroplasts could therefore reflect thylakoid reorganization or the impact of oxidative stress already at tillering stage under field capacity conditions.
Although the pigment content of the mutant plants was ~41% lower than in wild-type plants under control conditions, the photosynthetic rate was similar (Table 1). This seems illogical given that carotenoids and chlorophylls are required for the proper folding, assembly and stability of LHC apoproteins. However, we detected a greater xanthophyll de-epoxidation rate in the mutant as indicated by the relationship (0.5A + Z)/(VAZ). The zeaxanthin in the mutant may compensate for the lack of other xanthophylls in the antenna proteins, and may also preserve the structure of the photosystem and its membranes by dissipating photochemical energy as heat to prevent photo-oxidation. This may explain the abnormal chloroplast ultrastructure in the mutant but the near-normal photosynthetic rate. Our data suggest that, under control conditions, the mutant compensates for the lower pigment content and that minor oxidative stress has no detrimental effect on photosynthetic capacity or biomass accumulation.
However, under water stress conditions the level of de-epoxidation increases in wild-type plants, and both genotypes show similar values for the relationship (0.5A + Z)/(VAZ) (Table 3). This suggests that the mutant has a lower capacity for increasing the de-epoxidation rate. An alternative explanation is that the pool of convertible violaxanthin is rapidly depleted and only non-convertible violaxanthin remains to maintain the structural integrity of the antenna. Because xanthophylls help to prevent photo-oxidation, the mutant plants without additional de-epoxidated carotenoids to cope for the water stress damage may sense higher oxidative stress in the chloroplast under water stress. Thus, the pale-green mutant may be more susceptible since the ability of plants to withstand oxidative stress determines the level of overall stress tolerance.
Differences between the wild-type and mutant transcriptomes at anthesis (no stress)
There were only minor differences in gene expression between the wild-type and mutant plants under control conditions at anthesis, mirroring the similar net photosynthetic rates (Table 1). The few repressed genes in the mutant were predominantly related to ABA responses and the induced genes mostly represented pathogenesis-related (PR) proteins and peroxidases involved in plant defense (Table 5).
Apart from their role in defense, peroxidases may be modulated by cellular redox disturbances. In such cases, nuclear gene expression can be regulated by ROS-dependent signal transduction or low-molecular-weight antioxidants produced in the chloroplast. For example, the expression of PR proteins and peroxidases increases in ascorbate-deficient mutants. Since no lesion or pathogen attack was detected in our investigation, it is most likely that low-level of oxidative stress triggers defense gene expression in the pale-green mutant. Proteins involved in plant defense are also sensitive to hormones such as salicylic acid (SA), jasmonate (JA) and ethylene, but we did not find any differently-expressed genes related to these pathways in the mutant. Plant defense signaling pathways are also activated through WRKY transcription factors that bind to pathogen-response and wound-inducible promoter elements (W-boxes and GCC-like elements). The TaWRKY45 gene was induced in the mutant, and may help to activate some of the PR genes we also identified. This transcription factor confers resistance to Fusarium head blight (FHB) in wheat and increases the expression of PR-1 and PR-2 genes in rice to confer resistance to Pseudomonas syringae, salinity and drought tolerance, as well as reducing ABA sensitivity.
Genes encoding defense enzymes that synthesize secondary metabolites such as phenylpropanoids and flavonoids were also induced in the mutant, including naringerin, 2-oxoglutarate 3-dioxygenase and agmatine coumaroyltransferase (ACT), the latter representing the final step in the hydroxycinnamic acid amide (HCAA) synthesis pathway. Peroxidases can also couple agmatine conjugates to form antifungal compounds known as hordatines in barley and (to a lesser extent) in wheat. Finally, the blufensin-1 gene was also induced in the mutant. This is a negative regulator conferring sensitivity to powdery mildew in barley by modulating penetration resistance.
The repressed probe sets were related to systemic interactions with the environment, especially the ABA response, e.g. genes encoding dehydrins such as RAB proteins and other LEA proteins that accumulate during seed maturation and in vegetative tissues in response to abiotic stress. RAB is a LEA family dehydrin that is induced by ABA, and the Arabidopsis ortholog RAB18 is a useful marker for the loss of ABA sensitivity e.g. in plants that are deficient for the SWI/SNF chromatin-remodeling complex. The gene is upregulated when a relevant repressor is mutated, e.g. SUMO E3 ligase SIZ1. We hypothesize that the mutant induces a repressor that acts against the ABA pathway or downregulates an inducer of the same pathway. Indeed, two transcription factors were repressed in the mutant under control conditions, the Q-type C2H2 zinc finger transcription factor TaZP23 and the homeobox-leucine zipper protein ATHB-6. TaZP23 is upregulated in leaves and roots during drought stress and also following ABA treatment, suggesting a role in ABA-dependent gene regulation. ATHB-6 may act as a growth regulator in response to water deficit, and is induced by ABA, salinity and dehydration stress. It also regulates the ABA signaling pathway, acting downstream of ABI1.
The repressed probe sets also included two endonuclease genes related to the barley genes Bnuc and BEN1, which are expressed during apoptosis[67, 68]. Bnuc is also induced by salinity stress. The papain-like cysteine proteinase (HvPAP14) and TdS40 were also represented in the repressed probe sets, and these also facilitate senescence and apoptosis. HvPAP14 also generates mature forms of storage proteins in seeds, provides free amino acids during germination, and confers stress tolerance and defense against pathogens[72, 73]. TdS40 is induced during post-flowering senescence in wheat, and also by JA, SA and in chlorotic/necrotic leaf tissue following infection with Pyrenophora teres.
There were few differences between the mutant and wild-type plants in the absence of water stress at anthesis and no statistically significant differentially-expressed genes could be identified also at the tillering stage (data not shown). Therefore despite a clear difference in chlorophyll content and chloroplast structure, the mutation appears to have a minimal impact on metabolism and gene regulation, a finding that does not provide sufficient explanation for the yield loss detected under field conditions.
Water stress enhances the differences between mutant and wild-type plants
We investigated the molecular basis of the difference in yield between the two genotypes by exposing them to water stress, which is the most common environmental constraint for durum wheat. Most of the genes modulated by water stress were similarly regulated in both genotypes, suggesting a common response (Additional file3: Table S2). Genes involved in the synthesis of osmoprotectants were upregulated in both genotypes, whereas genes involved in ABA-dependent signaling, ethylene and JA synthesis were downregulated in both genotypes. Genes related to primary metabolism (amino acid and carbohydrate catabolism) were also similarly modulated, indicating that both genotypes use metabolic strategies to counter the effects of drought. In contrast, genes encoding photosystem components and enzymes representing carbohydrate metabolism and the tricarboxylic acid (TCA) cycle showed significant differences that may explain the difference in performance between wild-type and mutant plants in the field.
Photosynthesis and pigment synthesis
Both transcriptomic and proteomic analysis indicated that the components of photosystems I and II (PSI and PSII) and ATPase were differentially expressed in the wild-type and mutant plants under water stress. Several PSI proteins were underrepresented in the mutant, including reaction center (RC) protein subunit IV (PSI-E), which is required for optimal electron transport to ferredoxin and flavodoxin, and light harvesting complex LHCA1 (CAB 1B-21) (Table 9). PSI takes longer than PSII to recover from photo-inhibition because RC proteins are degraded instead of being repaired. Therefore, the loss of PSI RC proteins in the mutant is likely to exacerbate the damage caused by stress.
PSBO and PSBP maintain the structure of PSII and were specifically overrepresented in the mutant. CP29 was upregulated in both genotypes, but this was underrepresented in the mutant under stress conditions (Table 9). This protein binds LHCB proteins to the reaction center and confers photoprotection. The LHCB proteins were downregulated in both genotypes and underrepresented in the mutant. In contrast, the genes encoding LHCB proteins were upregulated in the mutant (Table 7). The beta and delta subunits of ATPase synthase were also overrepresented in the mutant, suggesting that the H+ flux across the thylakoid membrane is greater in the mutant, with a knock-on effect on ATP production.
Three genes responsible for the synthesis of carotenoid precursors were upregulated by stress specifically in the mutant: 1-deoxy-D-xylulose-5-phosphate synthase (DXS), hydroxymethylglutaryl-CoA synthase (MVA1) and phytoene synthase (PSY). One gene encoding glutamyl-tRNA reductase 1 (HEMA1), which is responsible for the synthesis of chlorophyll precursors, was also upregulated by stress specifically in the mutant. The pigment content of the mutant plants did not increase relative to the control but the activation of these genes could help to prevent a steeper decline in pigment content under water stress conditions (Table 3).
Redox changes in the chloroplast result in the stoichiometric adjustment of photosystem proteins and regulate the expression of both nuclear and plastid photosynthesis-related genes[80–82]. Because the mutant showed evidence of oxidative stress even under control conditions, the re-organization of the photosynthetic apparatus under water stress conditions could reflect the additional oxidative stress that occurs during the early stages of drought. The loss of photosynthetic proteins (especially LHCB) by the mutant in response to drought stress reduces the amount of captured energy and thus avoids the over-excitation of PSII which would generate additional ROS and cause oxidative damage.
Oxidative stress and ROS signaling
The balance between ROS synthesis and scavenging determines the level of oxidative stress in plants[54, 84]. ROS scavenging and detoxifying proteins such as peroxidases, glutathione transferases and glyoxalase were generally modulated in a similar manner in both genotypes, but more probe sets of this group were downregulated specifically in the mutant under stress. In the same direction, a glutaredoxin gene was upregulated in the wild-type plants under stress but was unaffected in the mutant (Table 7). Thylakoid ascorbate peroxidase and catalase were upregulated specifically in the mutant, suggesting a greater need for protection against ROS derived from chloroplast and peroxisome activity (Tables 7 and9). H2O2 can spread from the chloroplast to other organelles by diffusion through aquaporin-like PIP2 proteins which were downregulated in the mutant to a lesser extent than in wild-type plants (Additional file3: Table S2). These results indicate that the mutant experiences greater oxidative stress than wild-type plants under water stress conditions.
Two probe sets related to ROS and stomatal closure were specifically modulated in the mutant under water stress. The specific downregulation of glyoxalase gene suggests that the mutant is less able than wild-type plants to detoxify methylglyoxal using reduced glutathione, thus allowing the production of ROS and stomatal closure. Similarly, the specific upregulation of respiratory burst oxidase homologue D in the mutant (Table 7) suggests that the ABA-induced production of ROS in guard cells is promoted to encourage stomatal closure. However, the differential expression of genes that regulate the stomatal aperture did not result in any differences between the genotypes in terms of stomatal conductance under water stress.
ROS can induce Ca2+-dependent signaling that triggers interactions between different calcineurin B-like (CBL) sensors and their CBL-interacting protein kinase (CIPK) targets. Genotype-dependent differences were found in expression of these genes: whereas CIPK3 and CIPK10 were upregulated by stress in both genotypes, CIPK2 and CIPK6 were upregulated specifically in wild-type plants and CBL1 and CIPK23 were downregulated specifically in the mutant. These two proteins cooperate with CBL9 to activate K+-channels and regulate stomatal responses. CIPK10 and CIPK6 belong to the SnRK3-type family and mediate nutrient sensing, stress responses and ABA signaling[90, 91]. As a result of altered expression of CBLs and CIPKs different Ca2+-dependent downstream responses under water stress are expected between mutant and wild-type plants.
MPK3, which was specifically downregulated in the mutant, is also involved in ROS signaling. Once activated, it phosphorylates other stress-related proteins such as the transcription factors ZAT10 and AZF2, and enzymes in the ethylene synthesis pathway such as ACS, all of which were downregulated by stress in both genotypes. ZAT10 (STZ) is a transcriptional repressor responsive to chitin, which mediates the response to photo-oxidative stress and ROS.
Despite not showing damaged leaves or altered stomatal response, the data suggest that the pale-green mutant differs significantly from wild-type plants in terms of ROS sensing and signaling under water stress. ROS generated in the chloroplast could move more easily inside the cells and elicit a signal triggering specific responses.
Transcriptomic profiling revealed differences between the genotypes in the regulation of carbohydrate metabolism under water stress conditions (Table 7). The synthesis of osmoprotectants such as raffinose family oligosaccharides (RFOs) and trehalose was specifically upregulated in the wild-type plants, including the genes encoding galactokinase and raffinose synthases. RFOs protect membranes and scavenge ROS. Trehalose-phosphate synthase (TPS) was also specifically induced in wild-type plants and thus repressed in the mutant under stress. Trehalose and trehalose-6-phosphate act not only as osmolytes to stabilize proteins and membranes, but also as signaling molecules to regulate metabolism.
Two of the probe sets upregulated in the mutant plants encoded sucrose: sucrose 1-fructosyltransferase and sucrose:fructan 6-fructosyltransferase (SST1 and SFT6), which are responsible for synthesizing fructans from sucrose in the vacuole[95, 96]. The mutant plants therefore appeared to promote fructan synthesis for remobilization as free carbohydrates, which are directly correlated with the water-soluble carbohydrate content in wheat. Conversely, probe sets encoding cellulose synthase and cellulose synthase-like enzymes were upregulated specifically in the wild-type plants and xylosyltransferases were repressed in the mutant plants under stress (Table 7). These enzymes divert carbon to the synthesis of cell wall polysaccharides (cellulose, hemicellulose and xyloglucans) and their loss reduces growth and also inhibits photosystem activity and carbohydrate metabolism genes[97, 98]. This shows the importance of cell wall structure not only for growth and turgor changes under stress but also for carbohydrate metabolism. Furthermore, sucrose, RFOs and fructans can directly or indirectly scavenge ROS, as reviewed by Keunen et al..
The general picture that emerges from the modulation of carbohydrate metabolism suggests that the pale-green mutant is less able to deal with water stress than wild-type plants. Cell wall synthesis genes are activated in the wild-type plants whereas fructans are produced in the mutant. This correlates with the lower biomass of the mutant plants at anthesis and the lower yield under rain-fed conditions (Table 2).
TCA cycle and related pathways
The tricarboxylic acid (TCA) cycle generates energy and reducing power (ATP and NADPH) but consumes H2O and releases CO2. The substrates feeding the cycle are acetyl-CoA and other intermediates from lipid or protein catabolism. Acetyl-CoA is also produced from pyruvate (derived from glycolysis) by the action of the pyruvate dehydrogenase complex, and the subunit E1 gene representing this complex was upregulated in both genotypes under water stress. However, under conditions of carbon starvation or demand, additional strategies are used to supply TCA substrates[100, 101]. The wild-type plants appear to be more effective at these strategies given the specific upregulation of the following genes (Table 7):
Genes encoding products that convert amino acids into substrates suitable for the TCA cycle [102–104] were upregulated in both genotypes (Additional file 3: Table S2), but three were upregulated specifically in wild-type plants: methylcrotonoyl-CoA carboxylase, α and β subunits (MCCA and MCCB), and 2-oxoisovalerate dehydrogenase E2 (BCE2).
Genes encoding the NADP-malic enzyme (NADP-ME 3 and NADP-ME 4), which converts malate into pyruvate, and phosphoenolpyruvate carboxykinase (PCK1), which converts oxaloacetate into phosphoenolpyruvate, were upregulated specifically in the wild-type plants.
Genes encoding isocitrate lyase (ICL) and malate synthase (MLS) from the glyoxylate cycle were upregulated specifically in wild-type plants and thus, repressed in the mutant. This cycle produces succinate and malate acting as a bypass to avoid the loss of CO2 from the normal TCA cycle. It also metabolizes glyoxylate generated by photorespiration and prevents photo-inhibition.
Genes encoding glutamate dehydrogenase and glutamate decarboxylase (GDH and GDC) were specifically downregulated and consequently repressed in the mutant under stress. They participate in glutamate metabolism, which balances nitrogen and carbon use in the cell either by storing nitrogen as GABA, glutamate or glutamine, or by supplying carbon back to the TCA cycle. A functional GABA shunt is required for stress tolerance. In the pale-green mutant, carbon stored as glutamate and GABA may not be introduced efficiently into the TCA cycle, disrupting the nitrogen/carbon balance and reducing growth under carbon demand or starvation imposed by stress. GDH can also be inhibited by excess light or sucrose.
The feed-in to the TCA cycle and the ability to derive energy from it appear to be more limited in the mutant under stress, potentially explaining the differences between genotypes in terms of biomass and yield in the field experiments. The mutant mobilizes carbon less efficiently, with a significant impact on grain filling and yield as observed in our field experiments.
Regulators of carbon metabolism and signaling
Genes that regulate carbohydrate and nitrogen metabolism were also differentially expressed between the wild-type and mutant plants under water stress, helping to explain the basis of the metabolic differences between the genotypes (Table 7).
The BZIP63 gene was specifically downregulated and repressed in the mutant plants under stress. The transcription factor encoded by this gene is an important node in the glucose-ABA interaction network and is repressed by both ABA and glucose. It may promote starvation tolerance by regulating transcription of genes such as the raffinose synthase (RS6, upregulated in the wild-type plants).
AP2/EREBP proteins mediate hormone, sugar and redox signaling in the context of cold and drought stress. Each member has a unique expression profile and the individual proteins are differentially regulated by glucose. We found that RAP2.12 and RAP2.7 were specifically upregulated in wild-type plants and as a result appeared repressed in the mutant, and that RAP2.8 was specifically downregulated and also appeared repressed in the mutant, suggesting that the mutant and wild-type plants respond differently to the presence of sugars.
PI3P5K was upregulated specifically in the mutant plants. This enzyme increases the levels of phosphatidylinositol-5-phosphate, which functions as a signaling molecule during osmotic stress and can repress the expression of WRKY70. SAL1, which was downregulated specifically in the mutant, acts as negative regulator of drought tolerance by modulating the levels of different sugars and controlling the corresponding signaling pathways.
Since the results highlighted some modifications of signaling molecules, the enzymes of the metabolism could be fine-tuned specifically in each genotype to adjust the concentration of metabolites.
Under well-watered glasshouse conditions the reduced chlorophyll and carotenoids content of our pale-green mutant resulted on mild oxidative stress symptoms that could be compensated without altering photosynthesis and physiological performance by the induction of stress-response genes. When the availability of water was limited, we observed more extreme modulation of antioxidant enzymes and photosynthetic proteins in the mutant than in wild-type plants, suggesting that the latter were better protected against oxidative stress. Concomitant molecular changes in the mutant revealed that the regulation of carbohydrate metabolism differed from wild-type plants, particularly the energy-producing pathways and pathways feeding the TCA cycle. The modification of carbohydrate metabolism in the mutant resulted in lower biomass accumulation and hence reduced yield under field conditions. This yield reduction was significant (10%) under irrigation, but was exacerbated in rain-fed conditions (14%) in agreement with the gene expression study, thus indicating that the mutant plants were more sensitive to water scarcity. Our results support the conclusion that the pale-green mutant was less able to adapt to terminal water stress, sustaining that reduced pigment content may be disadvantageous in durum wheat under water limited conditions, although a case by case study would be required when working with other pale-green mutants.
We compared two durum wheat (Triticum turgidum L. var. durum) genotypes: the wild-type variety Borgia (from the cross IRTA-1004 x Bidi 17) and a derived pale-green mutant (MD-597). The mutant was selected in 2002 in the M2 generation after treating the wild-type variety with the mutagen sodium azide (N3Na) and backcrossing twice with the parent to remove other mutations. F1 offspring of backcrosses with the wild-type were similar to it, while the ratio of the character segregation was 3:1. Therefore, the mutation is most likely controlled by a recessive nuclear gene.The mutant plants used in our experiments were the stable eighth generation individuals (M8) derived from the isolated original mutant plant (Figure 4).
We used a split plot design with three replicates (eight pots per replicate) under controlled conditions in a glasshouse. The main plots were the water treatments (well-watered and water-stressed) and the subplots were the genotypes (wild-type and mutant). Two vernalized seedlings of the same genotype were transplanted into a 4.0-L pots (19 × 15 cm) filled with a blend of white and frozen through black sphagnum peat mixture (Klasmann-Deilmann GmbH), pH 6.5, conductivity 40 ms/m, NPK 14:16:18. Each genotype was planted in 48 pots, giving a total of 96 plants. All plants were grown with the following maximum/minimum temperatures and photoperiods: 13°C/7°C, 12 h daylight from emergence to the first node detectable stage, 18°C/10°C, 14 h daylight from then to heading, 20°C/15°C, 16 h daylight from heading to anthesis, and 25°C/20°C, 18 h daylight after anthesis. The pots were rotated twice a week to avoid temperature and light position effects.
Both genotypes were grown at two levels of soil water content: i) control, in which the soil was maintained at field capacity during the whole experiment, and ii) terminal water stress, in which irrigation was field capacity from sowing to booting (Zadoks stage 45;) and then reduced until maturity by watering the pots every second day with the amount of water necessary to maintain the pots at 1/3 of field capacity. The stress treatment was initiated by withholding water at the anthesis stage. The soil water content was monitored every second day using soil moisture probes (ECH2O dielectric aquameter, Decagon Devices, Inc. Pullman, WA, USA) in the center of the pots and was maintained by weighing pots daily and watering them accordingly. Data of the probes indicated that field capacity conditions corresponded to a soil water content around 50 ml water/100 ml of freshweight soil, while the water stress treatment reduced soil water content until ~13 ml water/100 ml soil (Table 1). The flag leaves of four plants per experimental unit were measured between 11 a.m. and 1 p.m. solar time at the booting and anthesis stages to determine (Table 1): i) the midday leaf water potential (ΨW) using a pressure chamber (Model 3005; Soil Moisture Equipment, Santa Barbara, CA); ii) leaf conductance (gs) using a steady-state porometer (LI-1600, Li-Cor, Lincoln, NE); and iii) net CO2 assimilation rate (An) using a portable infrared gas analyzer (Model LCA-4: Analytical Development, Hoddesdon, UK). For the expression profiling experiments, the flag leaf from five plants was collected for each of three biological replicates and for each water regime, between 11 a.m. and 1 p.m. solar time. The tissues were immediately frozen in liquid nitrogen and stored at -80°C. Each flag leaf was ground in liquid nitrogen and 100-mg and 40-mg aliquots were prepared for transcriptomic and proteomic analysis, respectively. Before RNA and protein extraction, five samples were pooled for each biological replicate.
Wild-type and mutant plants were grown in the field for three crop seasons (2010–2012) under irrigated (control) and rain-fed conditions in north-eastern Spain (latitude 41°38’N, longitude 0°23’ E, elevation 250 m a.s.l.) Soil type was Mesic Calcixerolic Xerochrept with a silty-clay texture and pH of 8.1. Plots were sown to a density of 500 seeds/m2 in randomized complete block designs with three replicates, each comprising eight rows (12.5 m in length and 15 cm apart). Field conditions are in Table 10. Irrigation was provided four times by flooding the first crop season (on Dec. 8th and at montly intervals from March to May) and by sprinkler irrigations the successive years in which frequencies were determined according to the environmental conditions (Additional file8: Figure S3). Fertilizer was applied as appropriate, and weeds, diseases and pests were controlled according to standard agricultural practices. Plant biomass was determined at anthesis using a sample of 0.5 m row length and crop dry weight (g/m2) following the methodology described by Villegas et al.. Plots were mechanically harvested at ripening and the yield was presented at 12% moisture level.
Three plantlets at the three fully-developed leaves stage (Zadoks stage 13;) from the glasshouse experiment, irrigated at field capacity, were collected from each genotype. Small segments (1 mm) from the central part of the third-leaf blades were fixed in 2.5% glutaraldehyde in 0.1 M phosphate buffer (pH 7.2) at 4°C. After three washes in 0.1 M phosphate buffer (pH 7.4) the samples were incubated in 1% osmium tetroxide and 2.5% potassium ferrocyanide for 2 h at 4°C. After ten washes for 10 min each with 0.1 M sodium acetate, and incubation for 30 min in 0.5% uranyl acetate followed by two washes for 10 min in 0.1 M sodium acetate at 4°C, the samples were dehydrated in increasing concentrations of ethanol followed by two washes with sodium acetate and finally embedded in epoxy EMBED-812 resin (Electron Microscopy Sciences) for three days at 60°C. Ultrathin sections (70 nm) contrasted with Reynolds’ lead citrate stain were examined using a Zeiss TEM 910 transmission electron microscope.
Pigments were quantified in four wild-type and mutant plants from each treatment regime harvested at anthesis. The blade of the flag leaf was dissected, weighed and immediately frozen in liquid nitrogen. The pigments from 100 mg lyophilized tissue were extracted with acetone for 20 min at 40°C, filtered and mixed with 9:1 hexane:ethyl ether and 2% NaCl. The organic fraction was evaporated in a nitrogen stream and argon was used for storage at -80°C. The pigments were separated on an ACQUITY Ultra Performance LC™ system (Waters, Milford, MA, USA) linked to a PDA 2996 detector (Waters, Milford, MA, USA) following the protocol of Taylor et al. but using a modified gradient.
Carotenoids were isolated using a YMC C30 carotenoid 3-μm, 2.0 × 100 mm HPLC column (Waters, Milford, MA) with a mobile phase comprising solvent A (methanol:water, 97:3, v/v) and solvent B (100% methyl tert-butylether), both solvents containing 0.05% triethylamine. The sample was warmed to 25°C and a volume of 10 μL was injected into the column, which also was thermostatically maintained at 25°C. The flow rate was set to 0.5 mL/min. The gradient program was set as follows: initial conditions 98% solvent A and 2% solvent B for 13 min; changed with a linear gradient to 62% solvent A and 38% solvent B in 1 min; hold for 2 min; changed with a linear gradient to 32% solvent A and 68% solvent B in 1 min; hold for 5 min; returned to initial conditions in 4 min, followed by equilibration for 5 min. The data analysis was done with MassLynx™ software version 4.1 (Waters, Milford, MA, USA). Carotenoids were identified by monitoring the order of elution from the column, the ultraviolet and visible spectra and the spectral fine structure and comparing them with standards and previous references. Standards trans-β-apo-8’carotenal, chlorophyll a and chlorophyll b were purchased from Sigma-Aldrich (USA), whereas lutein (L), zeaxanthin (Z), β-carotene (β-car), neoxanthin (N), violaxanthin (V) and antheraxanthin (A) were acquired from CaroteNature (Switzerland). The data were analyzed by ANOVA with the SAS-STAT package and the means were compared with a Tukey test at 5% probability level.
RNA isolation and array hybridization
Total RNA was extracted using Trizol reagent (Invitrogen, Carlsbad, CA, USA) and cleaned with RNeasy MinElute (QIAGEN) columns following the manufacturer’s instructions. Quality was determined using an Agilent Bioanalyzer 2100. The RNA samples were processed following the Affymetrix Genechip Expression Analysis Technical Manual. Single-stranded and double-stranded cDNAs were synthesized using the Affymetrix GeneChip® 3’ IVT Express Kit. Biotin-conjugated nucleotides were incorporated into aRNA by in vitro transcription to generate cRNA which was purified to remove unincorporated NTPs, salts, enzymes and inorganic phosphate before fragmentation and hybridization onto 3’ expression arrays. Specifically, we used the Affymetrix GeneChip® Wheat Genome Array. The arrays were washed and stained on an Affymetrix Fluidics Station followed by scanning with a GeneChip Scanner 3000.
GeneChip® quality analysis
The quality of microarray hybridizations was determined by checking the following parameters: i) the standard Affymetrix controls to evaluate labeling and hybridization (B2 oligonucleotides, PolyA controls such as lys, phe, thr and dap, and hybridization controls such as BioB, BioC, BioD and Cre); ii) RNA degradation plot relative to probe signal intensities of the Actin and GAPDH control genes, using the R library "simpleaffy"; iii) percentage of "Present calls" on the basis of the MAS 5.0 algorithm; and iv) R2 linear correlation coefficients among biological.
Data processing and analysis
Hybridization data were processed as recommended by Aprile et al.: i) normalization of the raw data by Robust Multi-Array Average (RMA) using the R package Affymetrix library; ii) MAS 5.0 algorithm on raw data in order to produce a detection call for each probe set "Present", "Marginal" or "Absent" (all probe sets that did not show a "Present" call in all repetitions of at least one sample were removed); iii) the filtered data were imported into the software Genespring GX7.3 (Agilent) for the analysis of differentially-expressed genes using a Welch t-test with Benjamini and Hochberg false discovery rate correction for multiple tests.
The differences in gene expression were considered to be significant at P < 0.05 and the induction or repression ratio was ≥2-fold. Principal component analysis (PCA) was then used to validate the role of the genotype and to verify the effect of water stress between the wild-type and mutant plants (Figure 2).
The terms induced or repressed were used when describing the comparison on gene expression between both genotypes independently of the water treatment (mutation effect). And the up/downregulated terms were used when comparing the two water treatments in the same genotype (stress effect).
In order to identify overrepresented gene classes within selected groups of genes compared to the entire array, the Arabidopsis thaliana best BLASTX annotations of the wheat probe sets were used as input for the MIPS Functional Catalogue Database (FunCatDB), which provides a statistical survey of the functional distribution of a given set of genes. Only the functional categories with a cut-off of P < 0.005 were considered. BLAST searches were carried out using Plant Expression Database (PLEXdb) where a complete list of additional annotation for all probe sets on the wheat genome array is available. The original microarray data have been deposited and are accessible in PLEXdb (http://www.plexdb.org) as experiment TA48: "Pale-green durum wheat mutant under terminal drought stress condition" and in NCBI's Gene Expression Omnibus (NCBI GEO;http://www.ncbi.nlm.nih.gov/geo) as experiment GSE45563.
Validation of array data with qRT-PCR
Quantitative real-time RT-PCR (qRT-PCR) was used to confirm the differential expression of the probe sets identified by microarray analysis. The total RNA isolated for the array experiment was treated with DNA-free™ DNase Treatment & Removal Reagents (Ambion, Life technologies) and cDNA was synthesized using the iScript cDNA synthesis kit (Bio-Rad, Hercules, CA, USA). PCR amplification was carried out in a 25 μL final volume containing 2× IQ SYBR Green Supermix (BioRad, Hercules, CA, USA), 500 nM of each specific forward and reverse primers and 10 ng of cDNA as the template, using the Bio-Rad CFX96 sequence detector system. The amplification program comprised a denaturation step at 95˚C for 3 min followed by 40 cycles of 95˚C for 10 s, 60˚C for 30 s and 72˚C for 20 s. The amplification products were validated by melt curve analysis using standard curves with an efficiency of 90–100%. Several repetitions were used to establish the efficiency of the standard curve for the internal reference gene. Efficiency was auto-calculated using the Bio-Rad CFX Manager software. Ct value and normalized expression were calculated using the ΔΔC(t) method. The no reverse transcriptase and no template negative controls were used to ensure absence of contamination and of nonspecific amplifications. Primer specificity was further confirmed by gel electrophoresis.
Specific qRT-PCR were used to verify nine probe sets corresponding to the following genes: Chitinase 3 (TaAffx.128418.43.S1_at), Blufensin 1 (TaAffx.26815.1.S1_at), Germin-like protein (Ta.169.1.S1_x_at), Naringenin,2-oxoglutarate 3-dioxygenase (TaAffx.6092.1.S1_at), WRKY45-like transcription factor (Ta.8614.2.S1_x_at), Agmatine coumaroyltransferase (Ta.8228.1.S1_at), RAB protein (Ta.2638.1.S1_at), Q-type C2H2 zinc finger protein (ZFP23) (TaAffx.98394.1.S1_at) and Dehydrin (WZY1-1; Ta.13255.1.S1_at). Primers were designed near the 3’, according to GenBank (NCBI) database sequences (Additional file9: Table S6). Data were calculated from the calibration curve and normalized using the expression curve of the Actin transcript (GQ339780) corresponding to probe set Ta.28253.1.S1_at, selected among several probe sets with low variation between samples as the internal reference gene (CV = 0.036). The qRT-PCR data were compared to the corresponding microarray expression values by mean of Pearson product–moment correlation coefficients (Additional file1: Figure S1).
Protein was extracted as described by Jiang et al. with modifications. The powder was resuspended in 2 mL ice-cold 10% trichloroacetic acid in acetone containing 0.07% dithiothreitol (DTT) and protease inhibitor cocktail (P9599 Sigma). After overnight incubation at -20°C, the mixture was centrifuged at 4 °C for 30 min at 14000 rpm. The pellet was washed four times with acetone containing 0.07% DTT and centrifuged for 5 min at 14000 rpm. The pellet was left to stand on ice until the acetone evaporated completely. Protein Extraction Reagent type 4 buffer (C0356 Sigma) was used to resuspend the pellet overnight at room temperature while mixing. Proteins were separated from insoluble components by centrifugation for 20 min at 14000 rpm and were stored at -80°C. The Bradford protein quantification method was used to quantify the proteins using bovine serum albumin as a standard.
2-D Electrophoresis and detection
For the first dimension separation, samples were mixed with rehydration buffer (7 M urea, 2 M thiourea, 1% C7BzO detergent, 40 mM Trizmabase, 50 mM DTT, 1% IPG buffer pH 3–10, and 0.002% bromophenol blue) to a total volume of 315 μL. Actively rehydrated (50 V/gel, 20˚C) 18-cm Bio-Rad IPG strips (pH 5–8) were loaded with protein and separated by isoelectric focusing using a Protean IEF Cell system (Bio-Rad) following the manufacturer’s instructions. The sequential gradient procedure was 500 V/linear for 30 min, 1000 V/linear for 1 h, 10000 V/linear for 1 h, and 10000 V/rapid until a total of 55,000 Vh. The current limit was 50 μA per IPG strip. After IEF separation, the gel strips were incubated for 15 min in the equilibration buffer (375 mM Tris base, 6 M urea, 20% glycerol, 2% SDS) containing 2% DTT, followed by 15 min in the same buffer containing 2.5% iodoacetamide instead of DTT. For the second-dimension separation, SDS-PAGE with 11% 22 × 20 cm polyacrylamide gels were loaded and separated in an Ettan DALTsix Electrophoresis Unit (GE Healthcare) with 0.25 M Tris–HCl (pH 8.8), 1.92 M glycine, 1% w/v SDS electrophoresis buffer, at 8 mA/gel overnight. Three replicate analytical gels, with 50 μg of proteins each, were run for each protein extract. The preparative gels were loaded with 350 μg of protein. After electrophoresis, the analytical and preparative SDS-PAGE gels were fluorescently stained with Flamingo or Oriole (Bio-Rad) respectively, and scanned with a Versadoc MP 4000 system (Bio-Rad).
Gel imaging and spot quantification were carried out using the PD-Quest Advanced 2D Gel Analysis software v8.0.1 (Bio-Rad) applying automatic spot detection and matching followed by manual/visual validation. The image was normalized by local regression and a scaling factor of 106 (ppm) was applied to avoid non-expression related variations in spot intensity. Spots showing more than 2-fold changes in intensity with a coefficient of variation <50% and with statistically significant differences by Student’s t-test (P < 0.05) were defined as differentially expressed.
Protein digestion and identification by mass spectrometry (MS)
The spots of interest were excised manually from preparative gels and digested with trypsin according to Perez-Hedó et al.. We then applied 1 μL of tryptic peptide solution to a MALDI plate, dried it at room temperature and covered it with 1 μL of saturated α-cyano-4-hydroxycinnamic acid prepared in 50% v/v ACN containing 0.1% TFA.
Proteins were identified using an AutoflexSpeed MALDI-TOF/TOF mass spectrometer (Bruker Daltonics) by peptide mass fingerprinting (PMF). For those spots lacking identification, 3–4 peptides were analyzed by MS/MS mode. Mass spectra (mode reflectron, MH+) were acquired by FlexControl v3.0 software (Bruker Daltonics), recording in the range 800–4500 Da, and the MS/MS information was obtained in LIFT (Laser-Induced Forward Transfer) mode. MS spectra were externally calibrated using Peptide Calibration Standard II (Bruker Daltonics). The peak lists obtained were compared against the Swiss-Prot, Trembl and non-redundant NCBI protein databases, using the MASCOT software package v2.3 (Matrix Sciences, UK). The search parameters were set as monoisotopic peptide masses, carbamidomethylation of cysteine and oxidation of methionine as fixed and variable modifications, respectively, one trypsin missed cleavage and a maximum of ±100 ppm for PMF peptide tolerance and ±0.4 Da for MS/MS tolerance. Protein identification was accepted for those results showing a probability-based MOWSE score with P < 0.05. We also analyzed MS and MS/MS combined spectra prepared using BioToolsv3.1 (Bruker Daltonics).
Availability of supporting data
The data set supporting the results of this article is available in the PLEXdb repository as experiment TA48.
Chlorophyll binding protein
CBL-interacting protein kinase
Late embryogenesis abundant protein
Light harvesting complex
Photosystem subunit O
Photosystem subunit P
Responsive to ABA
Related to AP2
Raffinose family oligosaccharides
Reactive oxygen species
Soil water content
Kumar S, Sehgal SK, Kumar U, Prasad PVV, Joshi AK, Gill BS: Genomic characterization of drought tolerance-related traits in spring wheat. Euphytica. 2012, 186 (1): 265-276. 10.1007/s10681-012-0675-3.
Pradhan GP, Prasad PVV, Fritz AK, Kirkham MB, Gill BS: High temperature tolerance in Aegilops species and its potential transfer to wheat. Crop Sci. 2012, 52 (1): 292-304. 10.2135/cropsci2011.04.0186.
Blum A, Shpiler L, Golan G, Mayer J: Yield stability and canopy temperature of wheat genotypes under drought-stress. Field Crops Res. 1989, 22 (4): 289-296. 10.1016/0378-4290(89)90028-2.
Araus JL, Slafer GA, Reynolds MP, Royo C: Plant breeding and drought in C-3 cereals: what should we breed for?. Annals of Botany. 2002, 89 (7): 925-940. 10.1093/aob/mcf049.
Gregersen PL, Culetic A, Boschian L, Krupinska K: Plant senescence and crop productivity. Plant Mol Biol. 2013, 82 (6): 603-622. 10.1007/s11103-013-0013-8.
Tardy F, Creach A, Havaux M: Photosynthetic pigment concentration, organization and interconversions in a pale green Syrian landrace of barley (Hordeum vulgare L., Tadmor) adapted to harsh climatic conditions. Plant Cell Environ. 1998, 21 (5): 479-489. 10.1046/j.1365-3040.1998.00293.x.
Carvalho RF, Takaki M, Azevedo RA: Plant pigments: the many faces of light perception. Acta Physiologiae Plantarum. 2011, 33 (2): 241-248. 10.1007/s11738-010-0533-7.
Pogson BJ, Albrecht V: Genetic Dissection of chloroplast biogenesis and development: an overview. Plant Physiol. 2011, 155 (4): 1545-1551. 10.1104/pp.110.170365.
Greiner S: Plastome mutants of higher plants. Genomics of Chloroplasts and Mitochondria. Edited by: Bock R, Knoop V. 2012, Netherlands: Springer, 237-266. [Govindjee and Sharkey TD (Series Editors): Advances in Photosynthesis and Respiration, vol 35.]
Leister D: Chloroplast research in the genomic age. Trends Genet. 2003, 19 (1): 47-56. 10.1016/S0168-9525(02)00003-3.
Bollivar DW: Recent advances in chlorophyll biosynthesis. Photosynthesis Res. 2006, 90 (2): 173-194.
Zhang H, Li J, Yoo J, Yoo S, Cho S, Koh H, Seo HS, Paek N: Rice chlorina-1 and chlorina-9 encode ChlD and ChlI subunits of Mg-chelatase, a key enzyme for chlorophyll synthesis and chloroplast development. Plant Mol Biol. 2006, 62 (3): 325-337. 10.1007/s11103-006-9024-z.
Chen Q, Wang LF, Su N, Qin HD, Niu HB, Wang JL, Zhai HQ, Wan JM: Photosystem 2 photochemistry and pigment composition of a yellow mutant of rice (Oryza sativa L.) under different irradiances. Photosynthetica. 2008, 46 (1): 35-39. 10.1007/s11099-008-0007-y.
Prina AR, Arias MC, Lainez V, Landau A, Maldonado S: A cytoplasmically inherited mutant controlling early chloroplast development in barley seedlings. Theor Appl Genet. 2003, 107 (8): 1410-1418. 10.1007/s00122-003-1391-0.
Morosinotto T, Bassi R, Frigerio S, Finazzi G, Morris E, Barber J: Biochemical and structural analyses of a higher plant Photosystem II supercomplex of a Photosystem I-less mutant of barley - Consequences of a chronic over-reduction of the plastoquinone pool. Febs Journal. 2006, 273 (20): 4616-4630.
Frigerio S, Campoli C, Zorzan S, Fantoni L, Crosatti C, Drepper F, Haehnel W, Cattivelli L, Morosinotto T, Bassi R: Photosynthetic antenna size in higher plants is controlled by plastoquinone redox state at post-transcriptional rather than transcriptional level. J Biol Chem. 2007, 282 (40): 29457-29469. 10.1074/jbc.M705132200.
Campoli C, Caffarri S, Svensson JT, Bassi R, Stanca AM, Cattivelli L, Crosatti C: Parallel pigment and transcriptomic analysis of four barley Albina and Xantha mutants reveals the complex network of the chloroplast-dependent metabolism. Plant Mol Biol. 2009, 71 (1–2): 173-191.
Falbel TG, Staehelin LA: Characterization of a family of chlorophyll-deficient wheat (Triticum) and barley (Hordeum vulgare) mutants with defects in the magnesium-insertion step of chlorophyll biosynthesis. Plant Physiol. 1994, 104 (2): 639-648. 10.1104/pp.104.2.639.
Falbel TG, Meehl JB, Staehelin LA: Severity of mutant phenotype in a series of chlorophyll-deficient wheat mutants depends on light intensity and the severity of the block in chlorophyll synthesis. Plant Physiol. 1996, 112 (2): 821-832. 10.1104/pp.112.2.821.
Liu X, Xu H, Zhang J, Liang G, Liu Y, Guo A: Effect of low temperature on chlorophyll biosynthesis in albinism line of wheat (Triticum aestivum) FA85. Physiol Plantarum. 2012, 145 (3): 384-394. 10.1111/j.1399-3054.2012.01604.x.
Bassi R, Caffarri S: Lhc proteins and the regulation of photosynthetic light harvesting function by xanthophylls. Photosynthesis Res. 2000, 64 (2–3): 243-256.
Jahns P, Holzwarth AR: The role of the xanthophyll cycle and of lutein in photoprotection of photosystem II. Biochim Biophys Acta-Bioenergetics. 2012, 1817 (1): 182-193. 10.1016/j.bbabio.2011.04.012.
Choudhury NK, Behera RK: Photoinhibition of photosynthesis: role of carotenoids in photoprotection of chloroplast constituents. Photosynthetica. 2001, 39 (4): 481-488. 10.1023/A:1015647708360.
Takahashi S, Badger MR: Photoprotection in plants: a new light on Photosystem II damage. Trends Plant Sci. 2011, 16 (1): 53-60. 10.1016/j.tplants.2010.10.001.
Petrov VD, Van Breusegem F: Hydrogen peroxide-a central hub for information flow in plant cells. AoB plants. 2012, pls014
Gill SS, Tuteja N: Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol Biochem. 2010, 48 (12): 909-930. 10.1016/j.plaphy.2010.08.016.
Fleury D, Jefferies S, Kuchel H, Langridge P: Genetic and genomic tools to improve drought tolerance in wheat. J Exp Bot. 2010, 61 (12): 3211-3222. 10.1093/jxb/erq152.
Leader DJ: Transcriptional analysis and functional genomics in wheat. J Cereal Sci. 2005, 41 (2): 149-163. 10.1016/j.jcs.2004.10.006.
Qin D, Wu H, Peng H, Yao Y, Ni Z, Li Z, Zhou C, Sun Q: Heat stress-responsive transcriptome analysis in heat susceptible and tolerant wheat (Triticum aestivum L.) by using Wheat Genome Array. BMC Genomics. 2008, 9: 432-10.1186/1471-2164-9-432.
Sutton F, Chen D, Ge X, Kenefick D: Cbf genes of the Fr-A2 allele are differentially regulated between long-term cold acclimated crown tissue of freeze-resistant and -susceptible, winter wheat mutant lines. BMC Plant Biol. 2009, 9: 34-10.1186/1471-2229-9-34.
Desmond OJ, Manners JM, Schenk PM, Maclean DJ, Kazan K: Gene expression analysis of the wheat response to infection by Fusarium pseudograminearum. Physiol Mol Plant Pathol. 2008, 73 (1–3): 40-47.
Winfield MO, Lu C, Wilson ID, Coghill JA, Edwards KJ: Cold- and light-induced changes in the transcriptome of wheat leading to phase transition from vegetative to reproductive growth. BMC Plant Biology. 2009, 9: 55-10.1186/1471-2229-9-55.
Wan Y, Underwood C, Toole G, Skeggs P, Zhu T, Leverington M, Griffiths S, Wheeler T, Gooding M, Poole R, Edwards KJ, Gezan S, Welham S, Snape J, Mills ENC, Mitchell RAC, Shewry PR: A novel transcriptomic approach to identify candidate genes for grain quality traits in wheat. Plant Biotech J. 2009, 7 (5): 401-410. 10.1111/j.1467-7652.2009.00405.x.
Krugman T, Chague V, Peleg Z, Balzergue S, Just J, Korol AB, Nevo E, Saranga Y, Chalhoub B, Fahima T: Multilevel regulation and signalling processes associated with adaptation to terminal drought in wild emmer wheat. Funct Integr Genomics. 2010, 10 (2): 167-186. 10.1007/s10142-010-0166-3.
Bahrman N, Le Gouis J, Negroni L, Amilhat L, Leroy P, Laine AL, Jaminon O: Differential protein expression assessed by two-dimensional gel electrophoresis for two wheat varieties grown at four nitrogen levels. Proteomics. 2004, 4 (3): 709-719. 10.1002/pmic.200300571.
Ouerghi Z, Remy R, Ouelhazi L, Ayadi A, Brulfert J: Two-dimensional electrophoresis of soluble leaf proteins, isolated from two wheat species (Triticum durum and Triticum aestivum) differing in sensitivity towards NaCl. Electrophoresis. 2000, 21 (12): 2487-2491. 10.1002/1522-2683(20000701)21:12<2487::AID-ELPS2487>3.0.CO;2-B.
Gao L, Yan X, Li X, Guo G, Hu Y, Ma W, Yan Y: Proteome analysis of wheat leaf under salt stress by two-dimensional difference gel electrophoresis (2D-DIGE). Phytochemistry. 2011, 72 (10): 1180-1191. 10.1016/j.phytochem.2010.12.008.
Xu J, Li Y, Sun J, Du L, Zhang Y, Yu Q, Liu X: Comparative physiological and proteomic response to abrupt low temperature stress between two winter wheat cultivars differing in low temperature tolerance. Plant Biology. 2013, 15 (2): 292-303. 10.1111/j.1438-8677.2012.00639.x.
Caruso G, Cavaliere C, Foglia P, Gubbiotti R, Samperi R, Lagana A: Analysis of drought responsive proteins in wheat (Triticum durum) by 2D-PAGE and MALDI-TOF mass spectrometry. Plant Sci. 2009, 177 (6): 570-576. 10.1016/j.plantsci.2009.08.007.
Ford KL, Cassin A, Bacic A: Quantitative proteomic analysis of wheat cultivars with differing drought stress tolerance. Front Plant Sci. 2011, 2: 44-
Jiang Y, Yang B, Harris NS, Deyholos MK: Comparative proteomic analysis of NaCl stress-responsive proteins in Arabidopsis roots. J Exp Bot. 2007, 58 (13): 3591-3607. 10.1093/jxb/erm207.
Peng Z, Wang M, Li F, Lv H, Li C, Xia G: A proteomic study of the response to salinity and drought stress in an introgression strain of bread wheat. Mol Cell Proteomics. 2009, 8 (12): 2676-2686. 10.1074/mcp.M900052-MCP200.
Aprile A, Mastrangelo AM, De Leonardis AM, Galiba G, Roncaglia E, Ferrari F, De Bellis L, Turchi L, Giuliano G, Cattivelli L: Transcriptional profiling in response to terminal drought stress reveals differential responses along the wheat genome. BMC Genomics. 2009, 10: 279-10.1186/1471-2164-10-279.
Landau AM, Lokstein H, Scheller HV, Lainez V, Maldonado S, Raul Prina A: A cytoplasmically inherited barley mutant is defective in Photosystem I assembly due to a temperature-sensitive defect in ycf3 splicing. Plant Physiol. 2009, 151 (4): 1802-1811. 10.1104/pp.109.147843.
Fristedt R, Willig A, Granath P, Crevecoeur M, Rochaix J, Vener AV: Phosphorylation of Photosystem II controls functional macroscopic folding of photosynthetic membranes in Arabidopsis. Plant Cell. 2009, 21 (12): 3950-3964. 10.1105/tpc.109.069435.
Simpson DJ, Vonwettstein D: Macromolecular Physiology of Plastids. 14. Viridis mutants in barley - Genetic, fluoroscopic and ultrastructural characterization. Carlsberg Res Commun. 1980, 45 (4): 283-314. 10.1007/BF02906180.
Cui Y, Jia Q, Yin Q, Lin G, Kong M, Yang Z: The GDC1 gene encodes a novel ankyrin domain-containing protein that is essential for grana formation in Arabidopsis. Plant Physiol. 2011, 155 (1): 130-141. 10.1104/pp.110.165589.
Austin JR, Frost E, Vidi P, Kessler F, Staehelin LA: Plastoglobules are lipoprotein subcompartments of the chloroplast that are permanently coupled to thylakoid membranes and contain biosynthetic enzymes. Plant Cell. 2006, 18 (7): 1693-1703. 10.1105/tpc.105.039859.
Bréhélin C, Kessler F, van Wijk KJ: Plastoglobules: versatile lipoprotein particles in plastids. Trends Plant Sci. 2007, 12 (6): 260-266. 10.1016/j.tplants.2007.04.003.
Bréhélin C, Kessler F: The plastoglobule: a bag full of lipid biochemistry tricks. Photochem Photobiol. 2008, 84 (6): 1388-1394. 10.1111/j.1751-1097.2008.00459.x.
Pogson BJ, Niyogi KK, Bjorkman O, DellaPenna D: Altered xanthophyll compositions adversely affect chlorophyll accumulation and nonphotochemical quenching in Arabidopsis mutants. Proc Natl Acad Sci USA. 1998, 95 (22): 13324-13329. 10.1073/pnas.95.22.13324.
Dall'Osto L, Lico C, Alric J, Giuliano G, Havaux M, Bassi R: Lutein is needed for efficient chlorophyll triplet quenching in the major LHCII antenna complex of higher plants and effective photoprotection in vivo under strong light. BMC Plant Biol. 2006, 6: 32-10.1186/1471-2229-6-32.
Wehner A, Grasses T, Jahns P: De-epoxidation of violaxanthin in the minor antenna proteins of photosystem II, LHCB4, LHCB5, and LHCB6. J Biol Chem. 2006, 281 (31): 21924-21933. 10.1074/jbc.M602915200.
Kar RK: Plant responses to water stress: role of reactive oxygen species. Plant Signal Behav. 2011, 6 (11): 1741-1745. 10.4161/psb.6.11.17729.
Baier M, Dietz KJ: Chloroplasts as source and target of cellular redox regulation: a discussion on chloroplast redox signals in the context of plant physiology. J Exp Bot. 2005, 56 (416): 1449-1462. 10.1093/jxb/eri161.
Colville L, Smirnoff N: Antioxidant status, peroxidase activity, and PR protein transcript levels in ascorbate-deficient Arabidopsis thaliana vtc mutants. J Exp Bot. 2008, 59 (14): 3857-3868. 10.1093/jxb/ern229.
Bahrini I, Sugisawa M, Kikuchi R, Ogawa T, Kawahigashi H, Ban T, Handa H: Characterization of a wheat transcription factor, TaWRKY45, and its effect on Fusarium head blight resistance in transgenic wheat plants. Breed Sci. 2011, 61 (2): 121-129. 10.1270/jsbbs.61.121.
Qiu Y, Yu D: Over-expression of the stress-induced OsWRKY45 enhances disease resistance and drought tolerance in Arabidopsis. Environ Exp Bot. 2009, 65 (1): 35-47. 10.1016/j.envexpbot.2008.07.002.
Nomura T, Ishizuka A, Kishida K, Islam AKMR, Endo TR, Iwamura H, Ishihara A: Chromosome arm location of the genes for the biosynthesis of hordatines in barley. Genes Genet Syst. 2007, 82 (6): 455-464. 10.1266/ggs.82.455.
Meng Y, Moscou MJ, Wise RP: Blufensin1 negatively impacts basal defense in response to barley powdery mildew. Plant Physiol. 2009, 149 (1): 271-285. 10.1104/pp.108.129031.
Battaglia M, Olvera-Carrillo Y, Garciarrubio A, Campos F, Covarrubias AA: The enigmatic LEA proteins and other hydrophilins. Plant Physiol. 2008, 148 (1): 6-24. 10.1104/pp.108.120725.
Saez A, Rodrigues A, Santiago J, Rubio S, Rodriguez PL: HAB1-SWI3B interaction reveals a link between abscisic acid signaling and putative SWI/SNF chromatin-remodeling complexes in Arabidopsis. Plant Cell. 2008, 20 (11): 2972-2988. 10.1105/tpc.107.056705.
Miura K, Lee J, Jin JB, Yoo CY, Miura T, Hasegawa PM: Sumoylation of ABI5 by the Arabidopsis SUMO E3 ligase SIZ1 negatively regulates abscisic acid signaling. Proc Natl Acad Sci USA. 2009, 106 (13): 5418-5423. 10.1073/pnas.0811088106.
Kam J, Gresshoff PM, Shorter R, Xue G: The Q-type C(2)H(2) zinc finger subfamily of transcription factors in Triticum aestivum is predominantly expressed in roots and enriched with members containing an EAR repressor motif and responsive to drought stress. Plant Mol Biol. 2008, 67 (3): 305-322. 10.1007/s11103-008-9319-3.
Soderman E, Hjellstrom M, Fahleson J, Engstrom P: The HD-Zip gene ATHB6 in Arabidopsis is expressed in developing leaves, roots and carpels and up-regulated by water deficit conditions. Plant Mol Biol. 1999, 40 (6): 1073-1083. 10.1023/A:1006267013170.
Himmelbach A, Hoffmann T, Leube M, Hohener B, Grill E: Homeodomain protein ATHB6 is a target of the protein phosphatase ABI1 and regulates hormone responses in Arabidopsis. EMBO J. 2002, 21 (12): 3029-3038. 10.1093/emboj/cdf316.
Zaina G, Morassutti C, De Amicis F, Fogher C, Marchetti S: Endonuclease genes up-regulated in tissues undergoing programmed cell death are expressed during male gametogenesis in barley. Gene. 2003, 315: 43-50.
Aoyagi S, Sugiyama M, Fukuda H: BEN1 and ZEN1 cDNAs encoding S1-type DNases that are associated with programmed cell death in plants. FEBS Lett. 1998, 429 (2): 134-138. 10.1016/S0014-5793(98)00563-8.
Muramoto Y, Watanabe A, Nakamura T, Takabe T: Enhanced expression of a nuclease gene in leaves of barley plants under salt stress. Gene. 1999, 234 (2): 315-321. 10.1016/S0378-1119(99)00193-6.
Helm M, Schmid M, Hierl G, Terneus K, Tan L, Lottspeich F, Kieliszewski MJ, Gietl C: KDEL-tailed cysteine endopeptidases involved in programmed cell death, intercalation of new cells, and dismantling of extensin scaffolds. Am J Bot. 2008, 95 (9): 1049-1062. 10.3732/ajb.2007404.
Martinez M, Cambra I, Carrillo L, Diaz-Mendoza M, Diaz I: Characterization of the entire cystatin gene family in barley and their target Cathepsin L-Like Cysteine-Proteases, partners in the hordein mobilization during seed germination. Plant Physiol. 2009, 151 (3): 1531-1545. 10.1104/pp.109.146019.
Martinez DE, Bartoli CG, Grbic V, Guiamet JJ: Vacuolar cysteine proteases of wheat (Triticum aestivum L.) are common to leaf senescence induced by different factors. J Exp Bot. 2007, 58 (5): 1099-1107. 10.1093/jxb/erl270.
Gilroy EM, Hein I, van der Hoorn R, Boevink PC, Venter E, McLellan H, Kaffarnik F, Hrubikova K, Shaw J, Holeva M, Lopez EC, Borras-Hidalgo O, Pritchard L, Loake GJ, Lacomme C, Birch PRJ: Involvement of cathepsin B in the plant disease resistance hypersensitive response. Plant J. 2007, 52 (1): 1-13. 10.1111/j.1365-313X.2007.03226.x.
Rampino P, Pataleo S, Falco V, Mita G, Perrotta C: Identification of candidate genes associated with senescence in durum wheat (Triticum turgidum subsp durum) using cDNA-AFLP. Mol Biol Rep. 2011, 38 (8): 5219-5229. 10.1007/s11033-010-0673-2.
Krupinska K, Haussuhl K, Schafer A, van der Kooij TAW, Leckband G, Lorz H, Falk J: A novel nucleus-targeted protein is expressed in barley leaves during senescence and pathogen infection. Plant Physiol. 2002, 130 (3): 1172-1180. 10.1104/pp.008565.
Rousseau F, Setif P, Lagoutte B: Evidence for the involvement of PSI-E subunit in the reduction of ferredoxin by Photosystem-I. EMBO J. 1993, 12 (5): 1755-1765.
Sonoike K: Photoinhibition of Photosystem I. physiol Plantarum. 2011, 142 (1): 56-64. 10.1111/j.1399-3054.2010.01437.x.
Yi X, Hargett SR, Frankel LK, Bricker TM: The PsbP protein, but not the PsbQ protein, is required for normal thylakoid architecture in Arabidopsis thaliana. FEBS Lett. 2009, 583 (12): 2142-2147. 10.1016/j.febslet.2009.05.048.
de Bianchi S, Betterle N, Kouril R, Cazzaniga S, Boekema E, Bassi R, Dall'Osto L: Arabidopsis mutants deleted in the light-harvesting protein Lhcb4 have a disrupted Photosystem II macrostructure and are defective in photoprotection. Plant Cell. 2011, 23 (7): 2659-2679. 10.1105/tpc.111.087320.
Pfannschmidt T, Schütze K, Brost M, Oelmüller R: A novel mechanism of nuclear photosynthesis gene regulation by redox signals from the chloroplast during photosystem stoichiometry adjustment. J Biol Chem. 2001, 276 (39): 36125-36130. 10.1074/jbc.M105701200.
Fey V, Wagner R, Bräutigam K, Pfannschmidt T: Photosynthetic redox control of nuclear gene expression. J Exp Bot. 2005, 56 (416): 1491-1498. 10.1093/jxb/eri180.
Pfannschmidt T, Bräutigam K, Wagner R, Dietzel L, Schröter Y, Steiner S, Nykytenko A: Potential regulation of gene expression in photosynthetic cells by redox and energy state: approaches towards better understanding. Ann Bot. 2009, 103 (4): 599-607.
Foyer CH, Neukermans J, Queval G, Noctor G, Harbinson J: Photosynthetic control of electron transport and the regulation of gene expression. J Exp Bot. 2012, 63 (4): 1637-1661. 10.1093/jxb/ers013.
Foyer CH, Noctor G: Redox homeostasis and antioxidant signaling: a metabolic interface between stress perception and physiological responses. Plant Cell. 2005, 17 (7): 1866-1875. 10.1105/tpc.105.033589.
Bienert GP, Moller ALB, Kristiansen KA, Schulz A, Moller IM, Schjoerring JK, Jahn TP: Specific aquaporins facilitate the diffusion of hydrogen peroxide across membranes. J Biol Chem. 2007, 282 (2): 1183-1192. 10.1074/jbc.M603761200.
Hoque TS, Uraji M, Ye W, Hossain MA, Nakamura Y, Murata Y: Methylglyoxal-induced stomatal closure accompanied by peroxidase-mediated ROS production in Arabidopsis. J Plant Physiol. 2012, 169 (10): 979-986. 10.1016/j.jplph.2012.02.007.
Kwak JM, Mori IC, Pei ZM, Leonhardt N, Torres MA, Dangl JL, Bloom RE, Bodde S, Jones JDG, Schroeder JI: NADPH oxidase AtrbohD and AtrbohF genes function in ROS-dependent ABA signaling in Arabidopsis. EMBO J. 2003, 22 (11): 2623-2633. 10.1093/emboj/cdg277.
Luan S: The CBL-CIPK network in plant calcium signaling. Trends Plant Sci. 2009, 14 (1): 37-42. 10.1016/j.tplants.2008.10.005.
Cheong YH, Pandey GK, Grant JJ, Batistic O, Li L, Kim B, Lee S, Kudla J, Luan S: Two calcineurin B-like calcium sensors, interacting with protein kinase CIPK23, regulate leaf transpiration and root potassium uptake in Arabidopsis. Plant J. 2007, 52 (2): 223-239. 10.1111/j.1365-313X.2007.03236.x.
Coello P, Hey SJ, Halford NG: The sucrose non-fermenting-1-related (SnRK) family of protein kinases: potential for manipulation to improve stress tolerance and increase yield. J Exp Bot. 2011, 62 (3): 883-893. 10.1093/jxb/erq331.
Chen L, Ren F, Zhou L, Wang Q, Zhong H, Li X: The Brassica napus Calcineurin B-Like 1/CBL-interacting protein kinase 6 (CBL1/CIPK6) component is involved in the plant response to abiotic stress and ABA signalling. J Exp Bot. 2012, 63 (17): 6211-6222. 10.1093/jxb/ers273.
Li G, Meng X, Wang R, Mao G, Han L, Liu Y, Zhang S: Dual-level regulation of ACC synthase activity by MPK3/MPK6 cascade and its downstream WRKY transcription factor during ethylene induction in Arabidopsis. Plos Genetics. 2012, 8 (6): e1002767-10.1371/journal.pgen.1002767.
Miller G, Shulaev V, Mittler R: Reactive oxygen signaling and abiotic stress. Physiol Plantarum. 2008, 133 (3): 481-489. 10.1111/j.1399-3054.2008.01090.x.
Krasensky J, Jonak C: Drought, salt, and temperature stress-induced metabolic rearrangements and regulatory networks. J Exp Bot. 2012, 63 (4): 1593-1608. 10.1093/jxb/err460.
Vijn I, Smeekens S: Fructan: more than a reserve carbohydrate?. Plant Physiol. 1999, 120 (2): 351-359. 10.1104/pp.120.2.351.
Xue G, McIntyre CL, Jenkins CLD, Glassop D, van Herwaarden AF, Shorter R: Molecular dissection of variation in carbohydrate metabolism related to water-soluble carbohydrate accumulation in stems of wheat. Plant Physiol. 2008, 146 (2): 441-454.
Park YB, Cosgrove DJ: Changes in cell wall biomechanical properties in the xyloglucan-deficient xxt1/xxt2 mutant of Arabidopsis. Plant Physiol. 2012, 158 (1): 465-475. 10.1104/pp.111.189779.
Wormit A, Butt SM, Chairam I, McKenna JF, Nunes-Nesi A, Kjaer L, O'Donnelly K, Fernie AR, Woscholski R, Barter MCL, Hamann T: Osmosensitive changes of carbohydrate metabolism in response to cellulose biosynthesis inhibition. Plant Physiol. 2012, 159 (1): 105-117. 10.1104/pp.112.195198.
Keunen E, Peshev D, Vangronsveld J, Van Den Ende W, Cuypers A: Plant sugars are crucial players in the oxidative challenge during abiotic stress: extending the traditional concept. Plant Cell Environ. 2013, doi:10.1111/pce.12061
Araújo WL, Tohge T, Ishizaki K, Leaver CJ, Fernie AR: Protein degradation - an alternative respiratory substrate for stressed plants. Trends Plant Sci. 2011, 16 (9): 489-498.
Mohammadi M, Kav NNV, Deyholos MK: Transcriptional profiling of hexaploid wheat (Triticum aestivum L.) roots identifies novel, dehydration-responsive genes. Plant Cell Environ. 2007, 30 (5): 630-645. 10.1111/j.1365-3040.2007.01645.x.
Fujiki Y, Ito M, Itoh T, Nishida I, Watanabe A: Activation of the promoters of Arabidopsis genes for the branched-chain alpha-keto acid dehydrogenase complex in transgenic tobacco BY-2 cells under sugar starvation. Plant Cell Physiol. 2002, 43 (3): 275-280. 10.1093/pcp/pcf032.
Taylor NL, Heazlewood JL, Day DA, Millar AH: Lipoic acid-dependent oxidative catabolism of alpha-keto acids in mitochondria provides evidence for branched-chain amino acid catabolism in Arabidopsis. Plant Physiol. 2004, 134 (2): 838-848. 10.1104/pp.103.035675.
Galili G: The aspartate-family pathway of plants: linking production of essential amino acids with energy and stress regulation. Plant Signal Behav. 2011, 6 (2): 192-5. 10.4161/psb.6.2.14425.
Xu H, Zhang J, Zeng J, Jiang L, Liu E, Peng C, He Z, Peng X: Inducible antisense suppression of glycolate oxidase reveals its strong regulation over photosynthesis in rice. J Exp Bot. 2009, 60 (6): 1799-1809. 10.1093/jxb/erp056.
Miyashita Y, Good AG: NAD(H)-dependent glutamate dehydrogenase is essential for the survival of Arabidopsis thaliana during dark-induced carbon starvation. J Exp Bot. 2008, 59 (3): 667-680. 10.1093/jxb/erm340.
Renault H, El Amrani A, Berger A, Mouille G, Soubigou-Taconnat L, Bouchereau A, Deleu C: γ-Aminobutyric acid transaminase deficiency impairs central carbon metabolism and leads to cell wall defects during salt stress in Arabidopsis roots. Plant Cell Environ. 2013, 36 (5): 1009-1018. 10.1111/pce.12033.
Melo-Oliveira R, Oliveira IC, Coruzzi GM: Arabidopsis mutant analysis and gene regulation define a nonredundant role for glutamate dehydrogenase in nitrogen assimilation. Proc Natl Acad Sci USA. 1996, 93 (10): 4718-4723. 10.1073/pnas.93.10.4718.
Matiolli CC, Tomaz JP, Duarte GT, Prado FM, Vieira Del Bem LE, Silveira AB, Gauer L, Guedes Correa LG, Drumond RD, Carvalho Viana AJ, Di Mascio P, Meyer C, Vincentz M: The Arabidopsis bZIP gene AtbZIP63 is a sensitive integrator of transient abscisic acid and glucose signals. Plant Physiol. 2011, 157 (2): 692-705. 10.1104/pp.111.181743.
Dietz K, Vogel MO, Viehhauser A: AP2/EREBP transcription factors are part of gene regulatory networks and integrate metabolic, hormonal and environmental signals in stress acclimation and retrograde signalling. Protoplasma. 2010, 245 (1–4): 3-14.
Price J, Laxmi A, St Martin SK, Jang JC: Global transcription profiling reveals multiple sugar signal transduction mechanisms in Arabidopsis. Plant Cell. 2004, 16 (8): 2128-2150. 10.1105/tpc.104.022616.
Ndamukong I, Jones DR, Lapko H, Divecha N, Avramova Z: Phosphatidylinositol 5-Phosphate links dehydration stress to the activity of ARABIDOPSIS TRITHORAX-LIKE Factor ATX1. Plos One. 2010, 5 (10): e13396-10.1371/journal.pone.0013396.
Wilson PB, Estavillo GM, Field KJ, Pornsiriwong W, Carroll AJ, Howell KA, Woo NS, Lake JA, Smith SM, Harvey Millar A, Von Caemmerer S, Pogson BJ: The nucleotidase/phosphatase SAL1 is a negative regulator of drought tolerance in Arabidopsis. Plant J. 2009, 58 (2): 299-317. 10.1111/j.1365-313X.2008.03780.x.
Zadoks JC, Chang TT, Konzak CF: Decimal code for growth stages of cereals. Weed Res. 1974, 14 (6): 415-421. 10.1111/j.1365-3180.1974.tb01084.x.
Villegas D, Aparicio N, Blanco R, Royo C: Biomass accumulation and main stem elongation of durum wheat grown under Mediterranean conditions. Ann Bot. 2001, 88 (4): 617-627. 10.1006/anbo.2001.1512.
Taylor KL, Brackenridge AE, Vivier MA, Oberholster A: High-performance liquid chromatography profiling of the major carotenoids in Arabidopsis thaliana leaf tissue. J Chromatogr A. 2006, 1121 (1): 83-91. 10.1016/j.chroma.2006.04.033.
Britton G, Liaaen-Jensen S, Pfander H: Carotenoids: Handbook. 2004, Birkhäuser: Basel
Aprile A, Federici C, Close TJ, De Bellis L, Cattivelli L, Roose ML: Expression of the H plus -ATPase AHA10 proton pump is associated with citric acid accumulation in lemon juice sac cells. Funct Integr Genomics. 2011, 11 (4): 551-563. 10.1007/s10142-011-0226-3.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4 (2): 249-264. 10.1093/biostatistics/4.2.249.
Gautier L, Cope L, Bolstad BM, Irizarry RA: affy - analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004, 20 (3): 307-315. 10.1093/bioinformatics/btg405.
Benjamini Y, Hochberg Y: Controlling the false discovery rate - a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995, 57 (1): 289-300.
Ruepp A, Zollner A, Maier D, Albermann K, Hani J, Mokrejs M, Tetko I, Guldener U, Mannhaupt G, Munsterkotter M, Mewes HW: The FunCat, a functional annotation scheme for systematic classification of proteins from whole genomes. Nucleic Acids Res. 2004, 32 (18): 5539-5545. 10.1093/nar/gkh894.http://mips.helmholtz-muenchen.de/proj/funcatDB,
Dash S, Van Hemert J, Hong L, Wise RP, Dickerson JA: PLEXdb: gene expression resources for plants and plant pathogens. Nucleic Acids Res. 2012, 40 (D1): D1194-D1201. 10.1093/nar/gkr938.http://www.plexdb.org,
Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29 (9): e45-10.1093/nar/29.9.e45.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3 (7): 0034.1-
Jiang Q, Chen H, Pan X, Pan Q, Shi Y, Li X, Zhang G, Wang Y, Xie S, Shen S: Proteomic analysis of wheat (Triticum aestivum L.) hybrid necrosis. Plant Sci. 2008, 175 (3): 394-401. 10.1016/j.plantsci.2008.05.017.
Bradford MM: A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein dye binding. Anal Biochem. 1976, 72 (1–2): 248-254.
Perez-Hedó M, Sanchez-Lopez I, Eizaguirre M: Comparative analysis of hemolymph proteome maps in diapausing and non-diapausing larvae of Sesamia nonagrioides. Proteome Sci. 2012, 10 (1): 58-10.1186/1477-5956-10-58.
A. Peremarti is supported by the AGROTECNIO Foundation. The support of the Efficient Use of Water Program of IRTA is acknowledged.
This study was partially funded by MICINN (Spain) under projects AGL2009-11187, AGL2012-37217 and RTA2009-00085-00, and by MIUR (Italy) under the project ISCOCEM. AP, DV, CR authors are part of the Centre CONSOLIDER INGENIO 2010 on Agrigenomics funded by the Spanish Ministry of Education and Science.
The authors declare that they have no competing interests.
AP carried out the molecular analysis and drafted the manuscript. CM participated in the transcriptomic analysis and interpretation of data. AA carried out the statistical analysis of transcriptomic data. ER carried out the microarray hybridization. CR, LC and DV helped to draft the manuscript and interpretation of data. CR released the mother variety and conceived the study. CR and DV obtained the mutant. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: Table S1: Significant functional categories represented in the probe sets modulated similarly in the wild-type and mutant plants at anthesis under water stress conditions. (DOC 90 KB)
Additional file 3: Table S2: Selected probe sets modulated similarly in the wild-type and mutant plants representing important processes occurring under water stress conditions. (DOC 128 KB)
Additional file 4: Table S3: Significant functional categories represented in the induced and repressed probe sets in the mutant plants under water stress conditions at anthesis. (DOC 34 KB)
Additional file 5: Figure S2: 2D-gel electrophoresis of leaf proteins from wild-type and mutant plants under water stress. (DOC 1 MB)
Additional file 6: Table S4: Protein identification in differently expressed spots by peptide mass fingerprinting (PMF). (DOC 195 KB)
About this article
Cite this article
Peremarti, A., Marè, C., Aprile, A. et al. Transcriptomic and proteomic analyses of a pale-green durum wheat mutant shows variations in photosystem components and metabolic deficiencies under drought stress. BMC Genomics 15, 125 (2014). https://doi.org/10.1186/1471-2164-15-125
- Chlorophyll content
- Oxidative stress
- Pale-green mutant
- Triticum turgidum L var. durum
- Water stress