Improper excess light energy dissipation in Arabidopsis results in a metabolic reprogramming
- Martin Frenkel†2,
- Carsten Külheim†1,
- Hanna Johansson Jänkänpää1,
- Oskar Skogström1,
- Luca Dall'Osto3,
- Jon Ågren4,
- Roberto Bassi3, 5,
- Thomas Moritz6,
- Jon Moen2 and
- Stefan Jansson1Email author
© Frenkel et al; licensee BioMed Central Ltd. 2009
Received: 21 May 2008
Accepted: 26 January 2009
Published: 26 January 2009
Plant performance is affected by the level of expression of PsbS, a key photoprotective protein involved in the process of feedback de-excitation (FDE), or the qE component of non-photochemical quenching, NPQ.
In studies presented here, under constant laboratory conditions the metabolite profiles of leaves of wild-type Arabidopsis thaliana and plants lacking or overexpressing PsbS were very similar, but under natural conditions their differences in levels of PsbS expression were associated with major changes in metabolite profiles. Some carbohydrates and amino acids differed ten-fold in abundance between PsbS-lacking mutants and over-expressers, with wild-type plants having intermediate amounts, showing that a metabolic shift had occurred. The transcriptomes of the genotypes also varied under field conditions, and the genes induced in plants lacking PsbS were similar to those reportedly induced in plants exposed to ozone stress or treated with methyl jasmonate (MeJA). Genes involved in the biosynthesis of JA were up-regulated, and enzymes involved in this pathway accumulated. JA levels in the undamaged leaves of field-grown plants did not differ between wild-type and PsbS-lacking mutants, but they were higher in the mutants when they were exposed to herbivory.
These findings suggest that lack of FDE results in increased photooxidative stress in the chloroplasts of Arabidopsis plants grown in the field, which elicits a response at the transcriptome level, causing a redirection of metabolism from growth towards defence that resembles a MeJA/JA response.
The photosynthetic light reaction provides plants with the energy required for growth and metabolism. However, due to the inherent properties of the excited pigment molecules and high-energy electron carriers, harmful side reactions can occur that can lead to the destruction of pigments, proteins and lipids and, unless these side reactions are controlled, death of the plant. Several mechanisms for regulating the photosynthetic light reaction have evolved, and the fact that all photosynthetic organisms have photoprotective processes, operating in different ways and over varying timescales, suggests that they are essential for survival. Feedback de-excitation (FDE), also called the qE – or energy-dependent – component of non-photochemical quenching (NPQ), is a photoprotective process of crucial importance for the plant. PsbS, a 22 kDa protein of the LHC superfamily, is necessary for the functioning of FDE . Shortly after an increase in light intensity, plants reduce the excitation pressure on photosystem II (PSII) by inducing FDE, changing the conformation of the photosynthetic antenna from an "energy transfer state" to a "quenched state", allowing excess energy to be dissipated as heat. Feedback de-excitation is reversible within a few minutes, whereas state transition, a second mechanism that can reduce excitation pressure, is reversible over slightly longer time-scales (ca. 15 minutes). Plants lacking FDE through mutations in the PsbS gene show no marked phenotypic deviations from their respective wild-type in controlled growth chamber conditions [1, 2], but exhibit reductions in fitness when grown under fluctuating light or field conditions . Other short-term responses to increases in light intensities include increases in cyclic electron transfer rates, activation of the Calvin cycle, and photorespiration. Long-term responses include reductions in effective light intensities by thickening and tilting of leaves, accumulation of anthocyanins, and movements of the chloroplasts [reviewed by: ]. A reduction in light harvesting can be achieved by reducing the size of the effective light-harvesting antenna, while energy utilization can be increased by up-regulation of the "dark" reactions. In addition, plants have scavenging systems that neutralize light-induced reactive oxygen species (ROS). ROS are generated, to some degree, at all light intensities, but especially when the intensity exceeds the capacity of the plant's photon utilization mechanisms.
Previously we have attempted to elucidate the mechanistic background of the reduced Darwinian fitness seen in field-grown plants lacking FDE by measuring selected physiological and metabolic parameters . In the study reported here we have used DNA microarrays and metabolic profiling to obtain more information on ways in which plants respond to variations in FDE capacity, using mutants with no or overexpressed levels of PsbS. Furthermore, we have examined the responses of genotypes varying in PsbS expression, both in the field and laboratory, to the effects of herbivory. Since plants grown under natural conditions show larger plant-to-plant variation than those grown under controlled conditions a tightly controlled experimental design was used in order to obtain reliable data. If the amount of PsbS in wild-type plants is sufficient to saturate relevant responses, PsbS over-expressers would be expected to have an identical phenotype to that of wild-type plants (unlike PsbS-deficient plants). However, it is also possible that increasing FDE capacity to levels seen in the PsbS over-expressers may be beneficial in some respects. Thus, comparing the three types of plants in the field was expected to give information that could help assess this possibility, and provide further insights into plant responses affected by differences in FDE capacity. Here we present data that provide new indications regarding the role of FDE at the whole plant level, including evidence of an overlap between the stress responses resulting from increased excitation of photosystem II and those induced by the jasmonate pathway.
PsbS overexpression did not have a significant impact on fitness
Arabidopsis thaliana genotypes Col-0 and oePsbSwere grown in an experimental garden during the summer of 2003, according to methods described by . The weather was drier and warmer, the overall seed yield of wild-type plants was lower and their mortality rates were higher – suggesting that the plants were under greater stress – during 2003 than during the cited study (2000 and 2001. Measurements on field-grown plants confirmed that the two genotypes had FDE levels corresponding to expectations based on their PsbS levels (data not shown). OePsbS plants produced significantly more siliques, but fewer seeds per silique, than wild-type plants and consequently there was no significant difference in total seed production per plant between the two genotypes; oePsbS plants producing on average 0.3 (± 8.3)% more seeds than wild type seed plants. In contrast, npq4 plants produced on average 36% less seeds than wild type plants in 2000/2001  and 2005 .
PsbS levels influence plant metabolism under field conditions but not under standard lab conditions
In our earlier experiments  the fitness of the npq4 mutant was found to be similar to that of wild-type plants under constant environmental conditions. Differences only occurred when plants were grown under variable conditions, with the most notable effects arising in the field where they were exposed to irregular and unpredictable changes in, for instance, light intensity. Therefore, we performed a detailed metabolomic investigation of wild-type, npq4 and oePsbS plants grown under natural conditions in a similar manner to that described above. We also used DNA microarrays to document the total pattern of gene expression for the three genotypes.
Metabolic profiling of leaves from field-grown plants was performed in the same way as for plants grown under controlled conditions. As before, PLS-DA was used in an attempt to classify the three genotypes, and the first two components from the resulting model (R2X = 0.37, R2Y = 0.89, Q2Y = 0.19) successfully separated the genotypes (data not shown). To validate these results another PLS-DA model was calculated, using only the npq4 and oePsbS sample-sets, and predictions were made for the wild-type data based on this model. This model predicted that the wild-type plants would be located between the npq4 and oePsbS genotypes in most respects (Fig. 1B). The new model showed that the main differences between the genotypes appeared to be trends relating to PsbS expression, with npq4 and oePsbS at opposite ends and the wild-type in the middle, strongly implying that PsbS levels influence the composition of the plant metabolome.
Thus, there were distinct between-genotype differences in the metabolic profiles of the plants grown in the field, but not in those of plants grown under controlled conditions in the lab, even though variations between replicates were greater among plants grown in the field (Fig. 1).
Major changes occurred in amino acid and carbohydrate metabolism
Metabolites with significantly differing levels in wild-type, npq4 and oePsbS
Similarly, carbohydrate metabolism was strongly influenced. Levels of galactose, raffinose and isomaltose were 30-, 9- and 8-fold higher, respectively in npq4 than in oePsbS plants. The levels of fructose and glucose also differed in a similar way, but not as markedly. Sucrose was the only carbohydrate detected that was more abundant in oePsbS than in npq4 plants. The trends in the amounts of glucose, fructose and sucrose between npq4 and wild-type were similar to those observed in plants grown in the field since germination, as detected by traditional metabolite analyses . These findings indicate that plants grown in climate chambers and transferred to the field after five days have strong metabolic similarities, at least in their carbohydrate metabolism, to those grown in the field since germination.
Major changes in composition were also detected amongst lipid metabolites. For example, the fatty acids diacylgalactosylglycerol and monogalactosyl-glycerol were more abundant in npq4, while concentrations of certain other lipids were higher in oePsbS plants. However, it should be noted that our metabolite profiling approach is not optimized for analysing non-polar compounds, and specific changes in lipids should be analyzed by dedicated lipid analysis techniques. Notable differences detected in the abundance of metabolites of other classes included findings that concentrations of two compounds involved in detoxification of ROS, alpha-tocopherol and dehydroascorbic acid (tentatively identified), were higher in the npq4 plants, while other metabolites with central roles in primary metabolism, such as malate, fumarate and citrate, were more abundant in oePsbS, as was spermidine.
Overall, these data show that the "core metabolism" of Arabidopsis, as well as some of its secondary metabolism, responded quite dramatically to variations in levels of PsbS in plants grown in the field.
PsbS levels influenced the transcriptome
When comparing three samples, three different pair-wise comparisons can be made or one can attempt to integrate data from all samples into a single analysis. There were strong ordinal relationships between PsbS expression and the concentrations of various metabolites, with the wild-type in an intermediate position between npq4 and oePsbS along the x-axis, and similar patterns were observed in the transcriptome data (Fig. 2). Therefore, in the following section we will mainly discuss the transcriptomal differences observed between npq4 and wild-type plants, since the same general trends were found in the other pair-wise comparisons.
Out of the 25392 spots on the full CATMA cDNA array, expression data were obtained from 15167 (the remaining gave too weak hybridization signals), with 526 showing statistically significant differential expression (False Discovery Rate, FDR, corrected p < 0.05) between npq4 and oePsbS. The 23 genes that had the most significantly different expression levels are listed in Additional file 3. The increased level of mRNA encoding galactinol synthase, which catalyzes the first committed step in the synthesis of raffinose-family oligosaccharides (RFOs) , in npq4 plants is consistent with the high accumulation of raffinose in this genotype. Other genes that were significantly up-regulated in npq4 include, inter alia, those encoding a dormancy-associated protein, a lipoxygenase, adenylosuccinate lyase (which is involved in purine biosynthesis), arginase (which catalyzes the conversion of arginine to urea and ornithine), and 1,4-alpha-D-glucan maltohydrolase (beta-amylase).
Expression levels of photosynthesis-related genes were higher in oePsbS
The development of Gene Ontology (GO) classifications and the introduction of MAPMAN software [14, 15] have made it easier to study coordinated gene responses involved in specific biochemical pathways or physiological processes. Here, MAPMAN was used to visualize changes in transcription, since we believe it provides the most useful interface and allows more readily comprehensible classifications than other more complex categorization methods.
MAPMAN software includes a list of "BINs", which contain genes that have similar functions. For example, BIN 1 contains all nuclear genes involved in photosynthesis, and BIN 1.1 contains those encoding products involved in the photosynthetic light reactions. MAPMAN can graphically display the results from a microarray experiment, illustrating the overall expression trends of genes from specific BINs, and apply statistical tools to determine the BINs and subBINs that have been significantly up- or downregulated. When classes of genes represented by BINs are analyzed quality filtering can be applied less stringently than when individual differentially expressed genes are analyzed, because a single unregulated gene may show differential expression purely through chance when large numbers of genes are analysed, as they are in DNA microarray experiments. However, these events are stochastic, and analysis of a class containing several genes is less likely to show a false positive result. Hence, if the filtering is too stringent it may result in the failure to detect classes of genes that, for example, are all only slightly down-regulated, but appropriate statistical methods can be used to detect and quantify the differential expression of groups of genes.
List of MAPMAN "BINs" showing evidence for differential expression
npq4 – wt
npq4 – oe
oe – wt
High exp in
PS lightreaction photosystem II
PS lightreaction photosystem II PSII polypeptide subunits
PS lightreaction photosystem I
PS lightreaction photosystem I PSIpolypeptide subunits
PS calvin cyle
major CHO metabolism
minor CHO metabolism raffinose family raffinose synthases
hormone metabolism jasmonate
hormone metabolism jasmonate synthesis-degradation
protein synthesis chloroplast/mito – plastid ribosomal protein
protein synthesis chloroplast/mito – plastid ribosomal protein plastid
transport peptides and oligopeptides
not assigned no ontology pentatricopeptide (PPR) repeat
BINs involved in the metabolism of diverse secondary metabolites contained many genes with differential transcript abundance, but the patterns were not consistent within whole BINs. For example, many genes in the BIN "Anthocyanin metabolism" were up-regulated in npq4, in accordance with our previous finding that under many conditions npq4 mutants accumulate more anthocyanins than wild-type plants , but the BIN as a whole was not significantly up-regulated. Similarly, the expression levels of many genes involved in both the anabolism and catabolism of amino acids also differed significantly between genotypes, although there were no overall genotype-related trends in these BINs.
The differential expression of many genes with functions relating to plastid metabolism was also corroborated by a different analysis, in which we extracted expression data for all the genes corresponding to the 1590 nuclear-encoded genes with plastid-related functions on the cDNA macroarray, as performed by . Analysis of the genes in this set (excluding those with A ≤ 7) showed that expression levels of genes representing seven of the 23 regulons, as identified by , differed significantly between genotypes (data not shown). For example, we obtained data for 60 of the 80 genes assigned to regulon 1 (which mainly consists of genes encoding photosystem subunits) and the expression levels of these genes were on average 35% higher in oePsbS than in npq4 plants.
Several genes in BINs related to redox regulation (e.g. the "Ascorbate/glutathione", "Glutaredoxin", "Peroredoxin" and "Dismutase/catalase" BINs) are involved in protection against photo-oxidative stress, and thus we expected their expression to be affected by changes in PsbS levels in field-grown plants. However, even though there were tendencies for some genes in the "Glutaredoxin" BIN to be up-regulated, and "Heme-oxidation proteins" to be down-regulated (data not shown), patterns of mRNA levels relating to genes in these BINs showed no overall differences between genotypes. The changes in PsbS levels may not have induced adjustments in the overall expression of genes encoding ROS protection systems, since the systems may have been sufficiently strong to cope with the oxidative stress experienced by plants of all three genotypes during the study period. However, these findings indicate that the main plant response to changes in PsbS expression, at least in their acclimation state at the time we sampled them, was not a localized increase in the capacity to cope with photo-oxidative stress, but rather a change in primary metabolism. Interestingly, many BINs relating to jasmonate (JA) synthesis were most strongly expressed in npq4 plants.
The metabolic reprogramming may be mediated, at least in part, by jasmonate
Overlap between microarray results
It also seemed plausible that such an effect on plant core metabolism could be mediated by plant hormone signalling. Therefore, we compared our gene list with sets of genes induced by various hormone treatments in the same way as above. Treatments with auxin, cytokinin, giberellin, ethylene or brassinolide do not apparently result in responses similar to that seen in npq4. The effect of abscissic acid showed some overlap (10% up, 2% down), whereas the response to salicylic acid was dissimilar (4% up, 10% down). However, plant responses to treatment with methyl jasmonate (MeJA) appear to be closest to those observed in field-grown plants lacking PsbS (17% up, 4% down), suggesting that at least part of the mechanism involved was a JA response.
An SDS-PAGE gel displaying peptides from wild-type and npq4 thylakoids showed the presence of LOX-C and AOS bands in npq4 samples even in NL conditions, while in the wild-type samples they were present in lower levels even when grown under HL. In each genotype the intensity of the LOX-C and AOS bands, and that of a lower molecular weight band (34 kDa), was proportional to the level of stress experienced by the plant as determined by increases in zeaxanthin, tocopherol and lipid peroxides. We conclude that LOX-C and AOS are found as thylakoid-associated proteins in the npq4 genotype and/or in HL conditions, and that their amount seems to be proportional to the level of light stress experienced by the plants.
Increased JA signalling in npq4plants may influence herbivore preferences
This study shows that under natural conditions Arabidopsis plants respond to changes in their FDE capacity quite dramatically. Our initial assumption was that a lack of PsbS, and hence FDE, would reduce the plants' Darwinian fitness as a consequence of changes in their photosynthetic performance. This may be true in part , but this study demonstrates that the plant's "core metabolism" (carbohydrate and amino acid metabolism) is also strongly affected, with large changes to the transcriptome occurring in plants lacking PsbS. Several amino acids (e.g. glutamine, glycine, aspartate and threonine) were highly depleted in plants lacking PsbS, whereas levels of others (e.g. leucine) were higher. Various carbohydrates also accumulated in plants lacking PsbS including (inter alia) galactose, raffinose and isomaltose, while sucrose levels were reduced. Knowledge of metabolic changes induced in Arabidopsis by different treatments is rather limited, but several of these compounds have been reported to accumulate under various kinds of stresses. For instance, there have been suggestions that the raffinose family oligosaccharides (RFOs) provide tolerance against drought, high salinity and cold stress . Hence, the metabolic reprogramming that occurred in plants lacking PsbS could be described in general terms as a stress response. Since transcriptional changes induced in Arabidopsis by various treatments have been intensively examined in several large-scale DNA microarray projects, they are far better characterized than overall metabolome changes. This allowed us to compare our transcriptome data with responses induced by a variety of biotic and abiotic stresses, as well as hormone treatments.
Comparing microarray data obtained using different platforms is far from straightforward but there are examples showing that cross-platform comparisons can give reliable results. For example, a benchmarking of GST-based CATMA arrays (which we used), and Affymetrics arrays indicated that the results were largely consistent . However, no study, including our own, has provided data on the degree of overlap between sets of transcripts in identical RNA samples identified as being upregulated in both CATMA and Affymetrics analyses. Hence, it is not easy to evaluate the absolute level of overlap between the set of genes induced in npq4 plants and those induced by other treatments, especially since our data were obtained under field conditions, while those presented in all other studies published to data were obtained from plants grown under controlled conditions in the lab. Nevertheless, the response seen in the npq4 genotype was most similar to those reportedly induced by ozone exposure or wounding, and also similar to a JA/MeJ response. The latter finding was further corroborated in four ways: comparison with other array data; detailed analysis of genes and transcripts associated with the JA/MeJa pathway in our plants; proteomic analysis; and measuring the JA levels within the plants. We hypothesize that a lack of PsbS in Arabidopsis plants leads to increased photooxidative stress and consequently induces a stress response. The signals involved are likely to be complex and to involve several hormones/signalling substances. However, since the response we found most closely resembled a JA/MeJA response in terms of the upregulated genes in microarray experiments, JA is probably one of the substances involved, and perhaps the most important. Additional studies using mutants with disturbances in different signalling pathways will hopefully provide further information regarding the nature of the signalling pathways involved.
It should be pointed out that plants probably have several mechanisms for sensing high light, and the signalling pathways involved in the responses observed in our study are unlikely to be identical to those that mediate general high light responses. For example, the effect of cryptochrome on high light-induced changes in gene expression  is apparently mediated by a different pathway from the response reported here, and the connection between ABA and high-light responses reported by  may involve different pathways too. However, it should be noted that there was some overlap between the genes induced in npq4 and those reportedly induced by ABA, although to a lower degree than the overlap seen with the JA response. Plants exposed to excessive light conditions need to modify many aspects of their growth and metabolism, since high light may be accompanied with any one or more of various other stressors, so multiple pathways could allow plants to respond to ecological challenges flexibly, with strong selective advantages.
In contrast to the abovementioned studies, here we compared genotypes that differed only in their capacity to dissipate excess light in the photosynthetic light-harvesting antenna, grown under identical environmental conditions. Thus, we focused on a single aspect of high-light signalling, mediated by signals directly generated at the photosynthetic reaction site. The results show that increased photooxidative stress at photosystem II is somehow sensed and a nuclear signal is generated that modifies gene expression. This in turn leads to a reprogramming of plant metabolism, with the implication that JA could be involved in the process. JA might even be the direct signal, and our finding that two enzymes involved in the JA biosynthesis pathway were enriched in the thylakoid membranes of PsbS-deficient plants may provide clues regarding the mechanism behind the response. Activation of the JA pathway is believed to occur at the site of JA biosynthesis; the chloroplastic membrane . Therefore, JA is locally synthesized in response to light stress, and is also amplified by a positive feedback loop . However, the level of JA was unaffected in npq4 plants unless they were provoked by herbivory. This does not exclude the direct involvement of the octadecanoid pathway in the response; increases could be induced in the flux through the pathway or in levels of either MeJA or JA-Ile, which may be the active species . However, it is also possible that plants lacking PsbS could be more "primed" and respond more actively to herbivory by increasing JA levels. At this point it is not possible to distinguish between these possibilities, but we hope to address these possibilities in future studies.
Regardless of the exact role that JA plays in the process, the metabolic reprogramming may be sufficient to explain the lower fitness observed in plants lacking PsbS, since JA-induced responses are known to reduce seed production . Furthermore, Arabidopsis plants treated with JA have an appearance that, in our experience, closely resembles that of plants exposed to moderate light stress, with anthocyanin accumulation, changes in leaf morphology and growth retardation (unpublished results). Overall, these findings suggest that plants lacking PsbS allocate resources away from reproduction towards defence. It has become increasingly apparent that stress responses are very complex. For instance, different herbivores may elicit specific responses in attacked plants and jasmonate signalling may affect different herbivores in different ways . Furthermore, although responses to various abiotic stresses share several common features, they also exhibit many highly specific features. Since most stress signalling has been shown to involve ROS in some way, one speculative possibility is that the changes we have observed could correspond to a more basal metabolic shift from reproduction and/or growth towards defence that increases the plant's ability to respond competently to stimuli elicited by specific stresses.
PsbS levels in wild-type plants do not appear to be "saturated", suggesting that selection acts against the production of higher levels of PsbS, even though it would increase the FDE capacity and hence potentially provide more photoprotection. Although this must theoretically hold true for the expression levels of all genes, it is not easy to demonstrate experimentally. An interesting implication of our data is that too much photoprotection could compromise "stress signalling", either in a more general sense, as discussed above, or more specifically through the involvement of JA/MeJA. Hence, plants with elevated PsbS levels could have less competence to cope with other stresses, and recent data indicate that this is indeed the case (Frenkel et al, in preparation). Overall, these findings suggest that there might be a tradeoff between photoprotection and tolerance to other stresses, as indicated by our data showing differences in herbivore preferences between wild-type plants and plants lacking PsbS.
Many mutants show marked phenotypic deviations from corresponding wild types, but mutations in most genes lead to more subtle changes in knock-out or antisense plants, with the lack of observable phenotypic deviations being the rule rather than the exception (see e.g. [28, 29]). Functional redundancy, of course, could be one reason for this, but it is clear that in many cases no observable phenotypic change may occur because the plants are not grown under suitable conditions for such changes to appear. Plant metabolic networks have highly complex architectures and the ability to compensate for enzyme deficiencies, or even an insufficiency in an entire pathway, using alternative pathways or substances that can perform similar functions, is likely to be of great adaptive value. Growing plants under naturally stressful conditions, rather than under controlled conditions in the lab, inevitably generates data that are complex and difficult to interpret. To our knowledge, the present study provides the first demonstration that metabolomic and transcriptomic analyses of mutants and over-expressers grown in the field can provide new and unexpected conclusions about individual gene products. We have been able to separate the genotypic differences from uncontrolled variation, and metabolomic analyses performed in the field and the laboratory clearly demonstrated that genotypic differences in expression patterns observed in field-grown material were not necessarily present under controlled conditions. It should be stressed that the responses described here were not detected under standard growth conditions in the lab. Therefore, we believe that studies like this, using plants grown under natural conditions, are essential for appreciating the full complexity of plants' genetic and metabolic composition and their interactions with the environment.
We have performed metabolomic, transcriptomic, proteomic, hormone and herbivore preference analyses of plants containing differing amounts of PsbS, both in the lab and the field. The results show that the level of PsbS is a determinant of plant fitness, metabolism and gene expression, perhaps through influencing the JA response.
Plant material and growth conditions
The plants used in the experiments represented three genotypes: Arabidopsis thaliana cv Columbia-0 (referred to as wild-type, wt), oePsbS  and npq4-1 . The npq4-1 mutant has a fast-neutron induced deletion of the gene encoding PsbS. Plants with the npq4-1 genotype had been backcrossed for three generations after the mutagen treatment. Plants overexpressing the PsbS gene were generated by Agrobacterium tumefaciens-mediated transformation of Arabidopsis thaliana (Col-0) plants, using the plasmid pXPL7 containing a 3198-bp fragment of the BAC clone F9J23. Plants of the line used here (designated line 5) were backcrossed for one generation after homozygous lines were identified (F3 generation).
For the fitness experiment, 40 plants representing each of the three genotypes were grown as described in . The plants were grown in an experimental garden under natural conditions; they received no fertilization and were watered ad lib when dry. Since oePsbS is a transgenic line we abided by the terms of our growing permit, which stipulates that all transgenic Arabidopsis plants grown in the field must be covered by a net from the time the first plant bolts. This ensures insects are excluded from the site and thus prevents accidental insect-mediated cross-pollination with wild plants. After flowering had ceased all the siliques were counted, and seeds were also counted in three randomly selected siliques from each plant. The mean number of seeds per silique was then multiplied by the number of siliques produced per plant to estimate total seed production.
For the metabolomic and transcriptomic measurements, the three genotypes (wild-type, npq4 and oePsbS) were pre-grown in a growth chamber under short day conditions (8-hour light periods with 400 μmol quanta m-2 s-1) with a 23°C/18°C temperature regime and relative humidity of 70%. All plants were moved to an experimental garden after five weeks of growth and then left there for a further five days. At the field site the photon flux density (PPFD) ranged from very low levels up to 600 W/m2 (ca. 2300 μmol quanta m-2 s-1) during the (ca. 15-hour) photoperiods and was highly variable. The temperature varied between 16° and 28°C, and the relative humidity (RH) between 30 and 100%. At 11 am on the fifth day, one leaf of each of 10 plants representing each genotype was harvested for metabolite measurements, and two pools of 3–4 leaf rosettes from each genotype were harvested for transcript analysis. Samples were frozen immediately after harvesting and stored at -80°C until further use.
For the metabolomic measurements from the constant conditions, the three genotypes were grown in a growth chamber under short day conditions (8-hour light periods with 200 μmol quanta m-2 s-1) at ca 23°C temperature and relative humidity of about 30%.
Approximately 10 mg of each sample was extracted and analyzed for metabolite profiling, according to  with some minor modifications. To the chloroform:MeOH:H2O (6:2:2) extraction mixture (1 ml) used for each sample, the following stable isotope reference compounds were added: [2H7]-cholesterol, [13C3]-myristic acid, [13C4]-hexadecanoic acid, [2H4]-succinic acid, [13C5, 15N]-glutamic acid, [13C5]-proline, [13C4]-disodium α-ketoglutarate, [13C12]-sucrose, [2H4]-putrescine, [2H6]-salicylic acid and [13C6]-glucose. The extraction was performed using an MM 301 Vibration Mill (Retsch GmbH & Co. KG, Haan, Germany) at a frequency of 30 Hz s-1 for 3 min after adding 3 mm tungsten carbide beads (Retsch GmbH & Co. KG, Haan, Germany) to each tube to increase the extraction efficiency. The extraction was followed by centrifugation in an Eppendorf centrifuge (Model 5417C) for 10 min at 14 000 rpm, after which 200 μl of the supernatant was transferred to a GC-vial and evaporated to dryness. The samples were then derivatized using 30 μL of methoxyamine hydrochloride (15 mg mL-1) in pyridine for 1 h at 70°C, followed by 16 h at room temperature, and subsequently trimethylsilylated by adding 30 μL of MSTFA with 1% TMCS and incubating for 1 h at room temperature. After silylation, 30 μL of heptane containing 45 ng μl-1 methyl octadecanoate was added. The samples were analyzed, together with blank control samples and a series of n-alkanes, according to , by GC/TOFMS (C12-C40), to facilitate calculation of retention indices .
One μL of each derivatized sample was injected splitless by an Agilent 7683 autosampler (Agilent, Atlanta, GA, USA) into an Agilent 6890 gas chromatograph equipped with a 10 m, 0.18 mm i.d. fused silica capillary column with a chemically bonded 0.18 μm DB 5-MS stationary phase (J&W Scientific, Folsom, CA, USA). The injector temperature was 270°C, the septum purge flow rate was 20 ml min-1 and the purge was turned on after 60 s. The gas flow rate through the column was 1 ml min-1, the column temperature was held at 70°C for 2 min, then increased by 40°C min-1 to 320°C and held there for 2 min. The column effluent was introduced into the ion source of a Pegasus III time-of-flight mass spectrometer (Leco Corp, St Joseph, MI, USA). The transfer line and the ion source temperatures were 250°C and 200°C, respectively. Ions were generated by a 70 eV electron beam at an ionization current of 2.0 mA, and 30 spectra s-1 were recorded in the mass range 50 to 800 m/z. The acceleration voltage was turned on after a solvent delay of 170 s. The detector voltage was 1660 V.
All non-processed MS-files from the metabolic analysis were exported from the ChromaTOF software in NetCDF format into MATLAB™ 6.5 software (Mathworks, Natick, MA, USA), where all data pre-treatment procedures, such as base-line correction, chromatogram alignment, data compression and Hierarchical Multivariate Curve Resolution (H-MCR) were performed using custom scripts according to . All manual integration was performed using ChromaTOF 2.00 software (Leco Corp, St Joseph, MI, USA) and all multivariate statistical investigations (PCA, PLS-DA) were performed using Simca software 10.5.0.0 (Umetrics, Umeå, Sweden). The following statistics for the PLS-DA models are discussed at various points throughout this paper: R2X is the cumulative modelled variation in X, R2Y is the cumulative modelled variation in Y, and Q2Y is the cumulative predicted variation in Y, according to cross-validation. The range of these parameters is 0–1, where 1 indicates a perfect model.
Significantly differing metabolites were detected by interpreting the first weight vector (w1), as described by , from the npq4/oePsbS PLS-DA model, together with the 95% confidence intervals calculated using jack-knifing .
Raw data were extracted from the Tiff images generated by ScanArray Express 2.1 using GenePix Pro 5.1 software (Axon Instruments, CA, USA). The threshold settings for the spot diameter resize feature were < 75% and > 150%, and the composite pixel intensity (CPI) was set to 300. The data were uploaded and normalized, filtered and analyzed within UPSC-BASE . The complete set of data can be found at http://www.upscbase.db.umu.se, experiment UMA-0060. In UPSC-BASE, UmeaSAMED restricted linear scaling (RLS) was performed to combine the different intensity scan data, thereby increasing the linear range of the signals http://www.umu.se/climi/bact/Microarray/R-libraries.htm. B-statistics were used to describe the quality of the technical and biological replicates. Two biological and two technical replicates were used and a total of 12 slides were hybridized with probes from the three genotypes (Fig 8). Before averaging the data for the biological replicates multivariate statistical investigations (PCA) were performed using Simca 10.5 software (Umetrics, Umeå, Sweden). For analyses involving MAPMAN software, datasets were exported from UPSC-BASE with a MAPMAN/AraCyc export tool. Multiple testing corrections were done in MAPMAN using the Wilcoxon Rank Sum Test and Benjamini Hochberg correction.
Thylakoid membranes were purified from wild-type and npq4 leaves grown at different light and temperature conditions as previously described . 1D and 2D SDS-PAGE was performed using a Tris-sulphate buffer system  in a 12–18% acrylamide gradient. Gels were stained with Coomassie R. Spots were excised, vacuum-dried and subjected to mass spectroscopy after trypsin digestion according to standard protocols. cDNAs for LOX-C and AOS were obtained from the Nottingham Arabidopsis Centre and expressed in E. coli using a pQE-50 vector including a 6-HIS tag at the N-terminus. The recombinant proteins were isolated using a Ni++ column and further purified by SDS-PAGE and electro-elution of Coomassie-stained bands. Purified proteins were injected into rabbits and antisera were obtained. Immunoblotting was performed following 1D or 2D SDS-PAGE using alkaline phosphatase-coupled antibodies for detection. Zeaxanthin and tocopherol contents were determined by HPLC analysis of acetone extracts, and levels of peroxidised lipids were determined by measuring malonyl-dialdehyde (MDA) derivatives, separated and quantified by HPLC according to methods described by .
We performed a cafeteria experiment to test for differences in herbivore preferences for the genotypes, as follows. We collected netted slugs (Deroceras reticulatum, Mullusca, Gastropoda) from the experimental site. Plants for the cafeteria experiments (pair-wise comparisons between the genotypes) were sown in 5 cm wide pots and then vernalised for 4–5 days at 4°C. After the vernalisation they were transferred to a growth chamber with short day conditions (8-hour light periods with 150 μmol quanta m-2 s-1), a 23°C/18°C temperature regime and constant relative humidity of 70% for ca. six weeks. After reaching a rosette size large enough to punch out several leaf discs of 1 cm diameter, they were brought to the experimental field site and exposed to natural, varying light conditions. Plants were taken for the experiments after 4–7 days in the field.
Eight leaf discs (1 cm diameter) from each of the two genotypes were taken and placed on wet filter paper in a 14.5 cm Petri dish, arranged in chess-board pattern, prior to introducing an individual of Deroceras. Pictures of the Petri dishes were taken after 15 hours, and the leaf area eaten during the experiment was quantified by image analysis using Scion Image http://www.scioncorp.com. Every cafeteria experiment consisted of 10 replicates in paired combinations of the genotypes, and repeated on 7 different days giving a total of 70 runs.
During the summer of 2007, we scored naturally occurring herbivory in our experimental garden. To examine whether age (and size) of the plant affected herbivory levels, we used plants that were both 6 and 21 days old. The plants that were 21 days old were transferred to the field seven days after the younger plants. All plants were scored for leaf damage four times during the summer (after 18, 25, 29 and 34 days in the field according to the ages of the younger plants). For the analysis, we defined plants that had lost at least one leaf or more as damaged by herbivores.
We thank Krishna Niyogi for providing seeds for the oePsbS line. This work was supported by grants from the Swedish Research Council (VR), the Swedish Research Council for Environment, Agricultural Sciences and Spatial Planning (FORMAS), the Wallenberg Consortium North and the Kempe Foundation. R.B. thanks GENEFUN for financial support.
- Li XP, Bjorkman O, Shih C, Grossman AR, Rosenquist M, Jansson S, Niyogi KK: A pigment-binding protein essential for regulation of photosynthetic light harvesting. Nature. 2000, 403 (6768): 391-395. 10.1038/35000131.PubMedView ArticleGoogle Scholar
- Havaux M, Niyogi KK: The violaxanthin cycle protects plants from photooxidative damage by more than one mechanism. Proc Natl Acad Sci USA. 1999, 96 (15): 8762-8767. 10.1073/pnas.96.15.8762.PubMedPubMed CentralView ArticleGoogle Scholar
- Külheim C, Ågren J, Jansson S: Rapid regulation of light harvesting and plant fitness in the field. Science. 2002, 297 (5578): 91-93. 10.1126/science.1072359.PubMedView ArticleGoogle Scholar
- Horton P, Ruban AV, Walters RG: Regulation of light harvesting in green plants. Annual Review of Plant Physiology and Plant Molecular Biology. 1996, 47: 655-684. 10.1146/annurev.arplant.47.1.655.PubMedView ArticleGoogle Scholar
- Kulheim C, Jansson S: What leads to reduced fitness in non-photochemical quenching mutants?. Physiologia Plantarum. 2005, 125 (2): 202-211. 10.1111/j.1399-3054.2005.00547.x.View ArticleGoogle Scholar
- Frenkel M, Bellafiore S, Rochaix JD, Jansson S: Hierarchy amongst photosynthetic acclimation responses for plant fitness. Physiologia Plantarum. 2007, 129 (2): 455-459. 10.1111/j.1399-3054.2006.00831.x.View ArticleGoogle Scholar
- Gullberg J, Jonsson P, Nordstrom A, Sjostrom M, Moritz T: Design of experiments: an efficient strategy to identify factors influencing extraction and derivatization of Arabidopsis thaliana samples in metabolomic studies with gas chromatography/mass spectrometry. Analytical Biochemistry. 2004, 331 (2): 283-295. 10.1016/j.ab.2004.04.037.PubMedView ArticleGoogle Scholar
- Jonsson P, Bruce SJ, Moritz T, Trygg J, Sjöström M, Plumb R, Granger J, Maibaum E, Nicholson JK, Holmes E, Antti H: Extraction interpretation and validation of information for comparing samples in metabolic LC/MS data sets. Analyst. 2005, 130 (5): 701-707. 10.1039/b501890k.PubMedView ArticleGoogle Scholar
- Trygg J, Wold S: O2-PLS, a two-block (X-Y) latent variable regression (LVR) method with an integral OSC filter. Journal of Chmeometrics. 2003, 17 (1): 53-64. 10.1002/cem.775.View ArticleGoogle Scholar
- Martens H, Martens M: Modified Jack-knife estimation of parameter uncertainty in bilinear modelling by partial least squares regression (PLSR). Food Quality and Preference. 2000, 11 (1–2): 5-16. 10.1016/S0950-3293(99)00039-7.View ArticleGoogle Scholar
- Schauer N, Steinhauser D, Strelkov S, Schomburg D, Allison G, Moritz T, Lundgren K, Roessner-Tunali U, Forbes MG, Willmitzer L, Fernie AR, Kopka J: GC-MS libraries for the rapid identification of metabolites in complex biological samples. FEBS Lett. 2005, 579 (6): 1332-1337. 10.1016/j.febslet.2005.01.029.PubMedView ArticleGoogle Scholar
- Fouad WM, Rathinasabapathi B: Expression of bacterial L-aspartate-alpha-decarboxylase in tobacco increases beta-alanine and pantothenate levels and improves thermotolerance. Plant Molecular Biology. 2006, 60 (4): 495-505. 10.1007/s11103-005-4844-9.PubMedView ArticleGoogle Scholar
- Peterbauer T, Richter A: Biochemistry and physiology of raffinose family oligosacharides and galactosyl cyclitols in seeds. Seed Science research. 2001, 11 (3): 185-197.Google Scholar
- Thimm O, Blasing O, Gibon Y, Nagel A, Meyer S, Kruger P, Selbig J, Muller LA, Rhee SY, Stitt M: MAPMAN: a user-driven tool to display genomic data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004, 37 (6): 914-939. 10.1111/j.1365-313X.2004.02016.x.PubMedView ArticleGoogle Scholar
- Usadel B, Nagel A, Thimm O, Redestig H, Blaesing OE, Palacios-Rojas N, Selbig J, Hannemann J, Piques MC, Steinhauser D, Scheible WR, Gibon Y, Morcuende R, Weicht D, Meyer S, Stitt M: Extension of the visualization tool mapman to allow statistical analysis of arrays, display of corresponding genes, and comparison with known responses. Plant Physiol. 2005, 138 (3): 1195-1204. 10.1104/pp.105.060459.PubMedPubMed CentralView ArticleGoogle Scholar
- Noren H, Svensson P, Stegmark R, Funk C, Adamska I, Andersson B: Expression of the early light-induced protein but not the PsbS protein is influenced by low temperature and depends on the developmental stage of the plant in field-grown pea cultivars. Plant Cell Environ. 2003, 26 (2): 245-253. 10.1046/j.1365-3040.2003.00954.x.View ArticleGoogle Scholar
- Biehl A, Richly E, Noutsos C, Salamini F, Leister D: Analysis of 101 nuclear transcriptomes reveals 23 distinct regulons and their relationship to metabolism, chromosomal gene distribution and co-ordination of nuclear and plastid gene expression. Gene. 2005, 344: 33-41. 10.1016/j.gene.2004.09.009.PubMedView ArticleGoogle Scholar
- Farmer EE, Ryan CA: Octadecanoid precursors of jasmonic acid activate the synthesis of wound inducible proteinase inhibitors. Plant Cell. 1992, 4 (2): 129-134. 10.1105/tpc.4.2.129.PubMedPubMed CentralView ArticleGoogle Scholar
- Taji T, Ohsumi C, Seki M, Kasuga M, Kobayashi M, Yamaguchi-Shinozaki K, Shinozaki K: Important roles of drought- and cold-inducible genes for galactinol synhtase in stress tolerance in Arabidopsis thaliana. Plant J. 2002, 29 (4): 417-426. 10.1046/j.0960-7412.2001.01227.x.PubMedView ArticleGoogle Scholar
- Allemeersch J, Durinck S, Vanderhaeghen R, Alard P, Maes R, Seeuws K, Bogaert T, Coddens K, Deschouwer K, Van Hummelen P, Vuylsteke M, Moreau Y, Kwekkeboom J, Wijfjes AH, May S, Beynon J, Hilson P, Kuiper MT: Benchmarking the CATMA microarray. A novel tool for Arabidopsis transcriptome analysis. Plant Physiology. 2005, 137 (2): 588-601. 10.1104/pp.104.051300.PubMedPubMed CentralView ArticleGoogle Scholar
- Kleine T, Kindgren P, Benedict C, Hendrickson L, Strand A: Genome-wide gene expression analysis reveals a critical role for CRYPTOCHROME1 in the response of arabidopsis to high irradiance. Plant Physiology. 2007, 144 (3): 1391-1406. 10.1104/pp.107.098293.PubMedPubMed CentralView ArticleGoogle Scholar
- Rossel JB, Walter PB, Hendrickson L, Chow WS, Poole A, Mullineaux PM, Pogson BJ: A mutation affecting ASCORBATE PEROXIDASE 2 gene expression reveals a link between responses to high light and drought tolerance. Plant Cell and Environment. 2006, 29 (2): 269-281. 10.1111/j.1365-3040.2005.01419.x.View ArticleGoogle Scholar
- Turner JG, Ellis C, Devoto A: The jasmonate signal pathway. Plant Cell. 2002, 14: S153-S164.PubMedPubMed CentralGoogle Scholar
- Sasaki Y, Asamizu E, Shibata D, Nakamura Y, Kaneko T, Awai K, Amagai M, Kuwata C, Tsugane T, Masuda T, Shimada H, Takamiya K, Ohta H, Tabata S: Monitoring of methyl jasmonate-responsive genes in Arabidopsis by cDNA macroarray: Self-activation of jasmonic acid biosynthesis and crosstalk with other phytohormone signaling pathways. DNA Res. 2001, 8 (4): 153-161. 10.1093/dnares/8.4.153.PubMedView ArticleGoogle Scholar
- Staswick PE, Tiryaki I: The oxylipin signal jasmonic acid is activated by an enzyme that conjugates it to isoleucine in Arabidopsis. Plant Cell. 2004, 16 (8): 2117-2127. 10.1105/tpc.104.023549.PubMedPubMed CentralView ArticleGoogle Scholar
- Baldwin IT: Jasmonate-induced responses are costly but benefit plants under attack in native populations. Proceedings of the National Academy of Sciences of the United States of America. 1998, 95 (14): 8113-8118. 10.1073/pnas.95.14.8113.PubMedPubMed CentralView ArticleGoogle Scholar
- Kessler A, Halitschke R, Baldwin IT: Silencing the jasmonate cascade: Induced plant defenses and insect populations. Science. 2004, 305 (5684): 665-668. 10.1126/science.1096931.PubMedView ArticleGoogle Scholar
- Ganeteg U, Klimmek F, Jansson S: Lhca5 – an LHC-type protein associated with photosystem I. Plant Molecular Biology. 2004, 54 (5): 641-651. 10.1023/B:PLAN.0000040813.05224.94.PubMedView ArticleGoogle Scholar
- Ganeteg U, Strand A, Gustafsson P, Jansson S: The properties of the chlorophyll a/b-binding proteins Lhca2 and Lhca3 studied in vivo using antisense inhibition. Plant Physiology. 2001, 127 (1): 150-158. 10.1104/pp.127.1.150.PubMedPubMed CentralView ArticleGoogle Scholar
- Li X-P, Müller-Moulé P, Gilmore AM, Niyogi KK: PsbS-dependant enhancement of feedback de-excitation protect photosystem II from photoinhibition. Proc Natl Acad Sci USA. 2002, 99 (23): 15222-15227. 10.1073/pnas.232447699.PubMedPubMed CentralView ArticleGoogle Scholar
- Hilson P, Allemeersch J, Altmann T, Aubourg S, Avon A, Beyon J, Bhalerao R, Bitton F, Caboche M, Cannoot B, et al: Versatile Gene-specific sequence tags for arabidopsis functional genomics: transcription profiling and reverse genomics application. Genome Res. 2004, 14: 2176-2189. 10.1101/gr.2544504.PubMedPubMed CentralView ArticleGoogle Scholar
- Sjodin A, Bylesjo M, Skogstrom O, Eriksson D, Nilsson P, Ryden P, Jansson S, Karlsson J: UPSC-BASE – Populus transcriptomics online. Plant Journal. 2006, 48 (5): 806-817. 10.1111/j.1365-313X.2006.02920.x.PubMedView ArticleGoogle Scholar
- Dall'Osto L, Caffarri S, Bassi R: A mechanism of nonphotochemical energy dissipation, independent from PsbS, revealed by a conformational change in the antenna protein CP26. Plant Cell. 2005, 17 (4): 1217-1232. 10.1105/tpc.104.030601.PubMedPubMed CentralView ArticleGoogle Scholar
- Bassi R, Dainese P: A Supramolecular Light-Harvesting Complex from Chloroplast Photosystem-Ii Membranes. European Journal of Biochemistry. 1992, 204 (1): 317-326. 10.1111/j.1432-1033.1992.tb16640.x.PubMedView ArticleGoogle Scholar
- Havaux M, Eymery F, Porfirova S, Rey P, Dormann P: Vitamin E protects against photoinhibition and photooxidative stress in Arabidopsis thaliana. Plant Cell. 2005, 17 (12): 3451-3469. 10.1105/tpc.105.037036.PubMedPubMed CentralView ArticleGoogle Scholar