Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Using Co-Expression Analysis and Stress-Based Screens to Uncover Arabidopsis Peroxisomal Proteins Involved in Drought Response

  • Jiying Li,

    Affiliations Department of Energy Plant Research Laboratory, Michigan State University, East Lansing, Michigan, United States of America, Genetics Graduate Program, Michigan State University, East Lansing, Michigan, United States of America

  • Jianping Hu

    huji@msu.edu

    Affiliations Department of Energy Plant Research Laboratory, Michigan State University, East Lansing, Michigan, United States of America, Plant Biology Department, Michigan State University, East Lansing, Michigan, United States of America

Abstract

Peroxisomes are essential organelles that house a wide array of metabolic reactions important for plant growth and development. However, our knowledge regarding the role of peroxisomal proteins in various biological processes, including plant stress response, is still incomplete. Recent proteomic studies of plant peroxisomes significantly increased the number of known peroxisomal proteins and greatly facilitated the study of peroxisomes at the systems level. The objectives of this study were to determine whether genes that encode peroxisomal proteins with related functions are co-expressed in Arabidopsis and identify peroxisomal proteins involved in stress response using in silico analysis and mutant screens. Using microarray data from online databases, we performed hierarchical clustering analysis to generate a comprehensive view of transcript level changes for Arabidopsis peroxisomal genes during development and under abiotic and biotic stress conditions. Many genes involved in the same metabolic pathways exhibited co-expression, some genes known to be involved in stress response are regulated by the corresponding stress conditions, and function of some peroxisomal proteins could be predicted based on their co-expression pattern. Since drought caused expression changes to the highest number of genes that encode peroxisomal proteins, we subjected a subset of Arabidopsis peroxisomal mutants to a drought stress assay. Mutants of the LON2 protease and the photorespiratory enzyme hydroxypyruvate reductase 1 (HPR1) showed enhanced susceptibility to drought, suggesting the involvement of peroxisomal quality control and photorespiration in drought resistance. Our study provided a global view of how genes that encode peroxisomal proteins respond to developmental and environmental cues and began to reveal additional peroxisomal proteins involved in stress response, thus opening up new avenues to investigate the role of peroxisomes in plant adaptation to environmental stresses.

Introduction

Peroxisomes are small and single membrane-delimited organelles that house numerous oxidative reactions connected to metabolism and development. These organelles are dynamic in nature, as their abundance, morphology and protein composition can be remodeled in response to developmental and environmental cues to adapt to the need of the organism [1,2,3]. Plant peroxisomes perform conserved functions such as β-oxidation of fatty acids and related metabolites and detoxification of reactive oxygen species (ROS), as well as plant-specific functions including photorespiration and metabolism of hormones such as jasmonate (JA) and auxin. Peroxisomes are crucial to virtually every developmental stage in plants, from embryogenesis, seedling development, vegetative and reproductive development, to senescence, and were recently shown to be involved in plant response to biotic and abiotic stresses [2,4]. The number of known peroxisomal proteins has risen to ~170 in Arabidopsis, largely due to recent peroxisomal proteome analyses followed by in vivo protein targeting verifications [5].

Peroxisomes possess many oxidative reactions that produce H2O2, as well as ROS-scavenging enzymes such as catalase and ascorbate-glutathione cycle enzymes [4,6]. ROS is a key component in stress responses [7]. Suppression of catalase 1 in tobacco resulted in necrotic lesions in high light and increased susceptibility to paraquat, salt and ozone [8]. Mutants of Arabidopsis catalase 2 develop photoperiod-dependent leaf lesions [9]. Evidence from melon, Arabidopsis and tobacco suggested the involvement of several peroxisomal photorespiratory enzymes, e.g., hydroxypyruvate reductase (HPR), serine:glyoxylate aminotransferase (SGT), alanine:glyoxylate aminotransferase (AGT), and glycolate oxidase (GOX) in immune response, possibly through ROS production [10,11,12].

Peroxisomes are also involved in stress response through mechanisms other than ROS homeostasis. Arabidopsis Ca2+-dependent protein kinase CPK1 is physically associated with peroxisomes and functions in a SA-dependent signaling pathway that leads to plant resistance to both fungal and bacterial pathogens [13,14]. Arabidopsis PEN2 is a peroxisome-associated myrosinase involved in callose deposition and glucosinolate hydrolysis necessary to generate antimicrobial products, thus is required for plant resistance against a broad spectrum of nonhost fungal pathogens, [15,16,17,18,19]. Furthermore, JA biosynthetic enzymes, some of which reside in peroxisomes, have been shown to affect systemic acquired resistance (SAR) to varying degrees [20]. It was suggested that the final step of SA biosynthesis, i.e., cinnamate to SA via the reduction of two carbons, may occur through β–oxidation in the peroxisome [21], thus making the peroxisome a potential player in SAR signaling. Interestingly, some virus species can hijack peroxisomes for viral RNA replication, causing the proliferation of peroxisome-like vesicle structures and leading to plant necrosis [22,23], which adds another layer of peroxisomal involvement in plant-pathogen interaction.

Despite these findings, there are still substantial knowledge gaps in the role of peroxisomes in stress response and how the functions of these peroxisomal proteins may be connected. To further identify peroxisomal proteins involved in plant response to various stress conditions, a peroxisome-centered systematic approach is needed. Recent advances in genome-wide transcriptomic and gene ontology enrichment analyses have provided valuable information on gene functions and mechanisms of biological processes. An important finding from these analyses is that genes functioning in the same pathway are often co-regulated by shared transcriptional regulatory systems and thus co-express across development and/or under many stress conditions [24]. To this end, we performed a genome-wide transcriptomic analysis of genes that encode peroxisomal proteins in Arabidopsis, trying to determine whether peroxisomal genes involved in the same biochemical pathways are co-expressed and whether we could identify new peroxisomal proteins involved in stress response using this type of in silico analysis. We followed up the in silico analysis with a pilot drought-based mutant screen, which identified the role of the peroxisomal LON2 protease and the photorespiratory enzyme hydroxypyruvate reductase 1 (HPR1) in drought resistance. Our study marks the beginning of systematic identifications of peroxisomal proteins involved in plant adaptation to stresses.

Materials and Methods

Plant materials and growth conditions

Arabidopsis thaliana ecotype Col-0 was used as wild type (WT). T-DNA insertion mutant lines were obtained from the Arabidopsis Biological Resource Center (ABRC; http://www.arabidopsis.org/) and confirmed by PCR genotyping. Seeds were sown in the soil, stratified in the dark at 4°C for 3 days, and plants were grown in a controlled growth chamber at 22°C under long-day conditions (16 hrs white light at 100 μmol photons m-2 s-1 and 8 hrs dark) for 3.5 weeks before drought treatment.

Microarray data analysis and heatmap visualization

Microarray datasets containing expression data of Arabidopsis peroxisomal genes from various tissues at different developmental stages were obtained from the AtGenExpress database, and expression data under biotic and abiotic stresses were downloaded from NCBI Gene Expression Omnibus (GEO) database (Table A in S1 File). Peroxisomal gene expression data obtained from various developmental stages were directly extracted from the whole-genome data and used for generating the heatmap. For data on biotic and abiotic stresses, log2-normalized data were extracted for peroxisomal genes from the whole-genome expression profile, using methods previously described [25]. Analysis was performed using the Bioconductor software [26] with the statistical computing language R (version 2.15.2). Normalization of gene expression values was carried out with the robust multi-array average (RMA) algorithm [27] implemented in the Affy package of Bioconductor. Statistical significance of the differential expression values were assessed with Linear models for microarray (limma) package [28]. Hierarchical clustering of the differentially expressed genes was visualized by creating heatmaps using the color palette package RColorBrewer and the gplots package [29].

Chlorophyll fluorescence measurements

Chlorophyll fluorescence images of intact plants were obtained from a custom-designed plant imager chamber, using a previously described method [30]. Plants in the pots were placed in the imaging chamber in the dark for 20 min for dark adaptation before minimal chlorophyll fluorescence Fo was measured. Later, maximal fluorescence Fm was measured when a saturating pulse of light was applied. Fv/Fm = (Fm-Fo)/Fm. Fluorescence images were analyzed by ImageJ [31].

Drought stress assays

For the drought tolerance screen, each selected mutant (two plants) and two WT plants were grown in the same pot under long day conditions (specified above) for 3.5 weeks, after which point plants stopped receiving water for 18 days before Fv/Fm measurement was conducted. For the follow-up analysis of the lon2 and hpr1 mutants, Fv/Fm measurement was repeated in the same way as in the primary screen, and watered plants were added as the control. Leaf samples from the drought-treated and control plants were harvested for chlorophyll content measurement, relative water content (RWC) and anthocyanin quantification as described previously [32].

For chlorophyll measurement, rosette leaves were weighed and placed into 2 ml 80% acetone in the dark for 3 days. Absorbance at 645 nm and 663 nm was measured using a spectrophotometer. Total chlorophyll content = (22.22 x A645 + 9.05 x A663) μg/ml x 2 ml /leaf fresh weight in mg.

To measure relative water content, rosette leaves were cut and immediately weighed as fresh weight (FW), and then placed in distilled deionized water at 4°C in the dark for 24 hrs, and the weight was recorded as turgid weight (TW). Then the rosette leaf sample was placed at 60°C for 2 days and the weight was recorded as dry weight (DW). Relative water content = (FW-DW) / (TW-DW) X 100%.

For anthocyanin measurement, rosette leaves were weighed, frozen by liquid nitrogen, and ground to powder. After adding 2 ml extraction buffer (1% HCl in methanol), the samples were placed at 4°C overnight. Later, an equal amount of chloroform was added, and the mixture was centrifuged for 5 min. After the top supernatant was transferred to a new tube, equal volume of 60% extraction buffer was added. Absorbance of each tube at 530 nm and 657 nm were measured with a spectrophotometer. Anthocyanin content = (A530-A657) /weight.

Results

Co-expression analysis of genes that encode peroxisomal proteins during development and in response to stresses

Microarray datasets containing expression data of Arabidopsis peroxisomal genes from various tissues at different developmental stages and under biotic and abiotic stresses were downloaded from the AtGenExpress database and NCBI Gene Expression Omnibus (GEO) database, respectively (Table A in S1 File). Developmental data were obtained from different tissues from seedlings, adult and senescing leaves, flowers, and siliques and seeds at various maturation stages. Abiotic stress conditions included high light, cold, hypoxia, drought, salt and the major stress hormone abscisic acid (ABA). Biotic stresses included the bacterial pathogen Pseudomonas syringae pv. Tomato (pst) DC3000, fungal pathogen Botrytis cinerea, and Pathogen-Associated Molecular Patterns (PAMPs) such as the bacterial flagellin 22 (flg22), bacterial elongation factor (elf18) and the fungal elicitor chitin (Table A in S1 File). Expression profiles of 160 peroxisomal genes (Table B in S1 File) were extracted from the whole-genome expression profile, clustered by hierarchical clustering analysis based on the extent of co-expression, and visualized by heatmaps.

Not surprisingly, many peroxisomal genes that function in the same metabolic pathways are co-regulated during development (Fig 1). For example, genes that encode glyoxylate cycle enzymes isocitrate lyase (ICL), malate synthase (MLS), and citrate synthase 1 (CSY1) are clustered together and co-up-regulated during seed maturation but co-repressed in other developmental stages (Fig 1). This is consistent with the fact that the glyoxylate cycle, a pathway that converts the β–oxidation product acetyl-CoA to succinate and malate to be used for gluconeogenesis, is primarily if not exclusively active in seeds and early seedling development [2,33]. In agreement with their roles in photorespiration, which recycles 2-phosphoglycolate produced by the oxygenase activity of RuBisCO back to the Calvin-Benson cycle [34], the expression of genes for the peroxisomal photorespiratory enzymes hydroxypyruvate reductase 1 (HPR1), glycolate oxidase 1 and 2 (GOX1 and GOX2, which are indistinguishable in microarrays due to high sequence identity), and peroxisomal malate dehydrogenase MDH2 is high in vegetative tissues but diminished during seed development. In contrast to these three clustered genes, genes that encode two other photorespiratory enzymes, glutamate: glyoxylate aminotransferase 1 (GGT1) and alanine: glyoxylate aminotransferase 1 (AGT1), are expressed during seed development as well, indicating that photorespiration may not be the only process that these two enzymes participate. Interestingly, genes that encode the JA biosynthetic enzymes 12-oxophytodienoate reductase 3 (OPR3) and OPC-8:0 CoA ligase 1 (OPCL1), disease related protein PEN2, and the small heat shock protein ACD31.2 are clustered together and also co-expressed with the photorespiration genes HPR1, GOX1/2 and MDH2 throughout development. This pattern indicates that these stress-related and development-related (in the case of JA biosynthesis) genes might be under similar regulatory circuitry as those photorespiration genes. As the major H2O2 detoxification enzymes, the three catalases are mostly constitutively expressed throughout development (Fig 1).

thumbnail
Fig 1. Heatmap of transcript levels of peroxisomal genes in various developmental stages.

Absolute gene expression values downloaded from the AtGenExpress database were used for heatmap generation. Genes discussed in the text are in red. Developmental stages include: 1. seedling_cotyledons; 2. seedling_hypocotyl; 3. seedling_leaves1+2; 4. adult_leaves; 5. senescing leaves; 6. flower; 7. silique_stage3; 8. silique_stage4; 9. silique_stage5; 10. seed_stage6; 11. seed_stage7; 12. seed_stage8; 13. seed_stage9; 14. seed_stage10.

https://doi.org/10.1371/journal.pone.0137762.g001

Transcriptional reprogramming during stresses is an important mechanism to confer stress tolerance [35,36]. Many genes that encode peroxisomal proteins are significantly regulated by abiotic stresses (Fig 2A). Among them, the small heat shock protein-encoding gene AtHsp15.7 showed >150-fold increase in expression under high light (Figure A in S1 File) and had to be removed from the heatmap in Fig 2A to prevent it from masking the changes in other genes in the heatmap. Genes encoding the peroxisomal proliferation factors PEX11b, PEX11c and PEX11d are all up-regulated by hypoxia (Fig 2A), which is consistent with a previous finding that hypoxia stress can rapidly stimulate peroxisomal extension over endoplasmic reticulum [37]. CAT2 and CAT3 expressions are also significantly up-regulated during drought stress (Fig 2A), consistent with their role as major ROS detoxification enzymes in stress response [38]. The glyoxylate cycle genes ICL, MLS and CSY1 are again clustered (Fig 2A), suggesting the tight regulation of this pathway by abiotic stress factors. Genes encoding the photorespiration enzymes GOX1/2, HPR1, MDH2, AGT1 and the chloroplast/peroxisome dual localized organelle division protein dynamic related protein 5B (DRP5B) are clustered (Fig 2A), raising the interesting possibility that photorespiration and the proliferation of peroxisomes are co-regulated during plant adaptation to abiotic stresses.

thumbnail
Fig 2. Heatmaps showing peroxisomal gene expression under stress conditions.

Expression values are log2 normalized fold changes from untreated plants. Genes discussed in the text are in red. (A) Gene expression under abiotic stresses. Genes subjected to mutant analysis are underscored. (B) Gene expression under biotic stresses.

https://doi.org/10.1371/journal.pone.0137762.g002

In response to biotic stresses, genes for some peroxisomal proteins previously shown to be involved in defense exhibited strong transcriptional reprogramming. For example, PEN2 (Penetration 2) is induced by two PAMPs, flg22 and chitin, but not by elf18 or the two pathogens (Fig 2B), supporting its major role as a myrosinase in PAMP-triggered immunity [15,16]. The two JA biosynthetic genes OPR3 and OPCL1 are co-up-regulated by flg22, chitin, P. syringae and B. cinerea (Fig 2B), consistent with JA’s role as an important defense hormone [39]. Interestingly, photorespiratory genes such as HPR1, CAT2, GOX1/2, MDH2, AGT1 and GGT2 are co-down-regulated by elf18, P. syringae and B. cinerea (Fig 2B), which is in agreement with the idea that photorespiration may play a defense role against pathogens through H2O2-dependent and -independent metabolisms [40]. The peroxisomal elongation factor gene PEX11b is again co-expressed with several photorespiratory genes during biotic stresses (Fig 2B), which is in accordance with its co-expression with photorespiratory genes in response to light [41] and the role of PEX11b in inducing peroxisomal proliferation during dark-to-light transition [42]. This data also indicated a potential need to increase peroxisomal abundance during pathogen defense.

A drought tolerance mutant screen revealed the role of the LON2 protease and the photorespiratory enzyme hydroxypyruvate reductase 1 (HPR1) in drought response

Based on the rule of “guilt-by-association” [24], those peroxisomal genes that showed significant up-regulation of transcript levels by some stresses have the potential to play a role in these specific conditions. To test this hypothesis, we decided to choose a stress condition, under which significant regulation of expression is seen for the highest number of peroxisomal genes, to screen for mutants with altered response.

Among the abiotic and biotic stress conditions examined, drought and the bacterial PAMP elf18 trigger expression changes to the highest number of peroxisomal genes (Figure B in S1 File). Drought is one of the most common environmental stresses that limit plant growth and development. Plants have evolved sophisticated adaptive drought tolerance mechanisms, including increased level of water transporting capacity, decrease of evaporative water, up-regulation of osmolytes and chaperone proteins, activation of Ca2+-dependent, ABA-dependent and other signaling pathways, and regulation of the transcript levels of the genes involved [43]. Mutant plants defective in these processes may exhibit increased drought sensitivity, such as increased water loss and ion leakage, decrease of photosynthesis rate, degradation of chlorophyll, and eventually cell death and plant withering [43]. As such, we used a drought tolerance assay as an initial screen to test the prediction from in silico analysis.

We have a collection of Arabidopsis mutants, which has facilitated us in discovering functions of newly identified peroxisomal proteins in previous studies [25,44,45]. To identify peroxisomal proteins involved in drought stress response, we first selected 26 mutants for 18 genes, most of which showed transcript level changes under drought or the drought stress hormone ABA. These included the up-regulated genes CAT2, GOX3, Hsp15.7, CSY3, Macrophage Migration Inhibitory Factor 1 (MIF1) and LON2 protease, and the down-regulated genes polyamine oxidase PAO2, thiolase KAT5 and acyl-CoA activating enzyme AAE14. We also included mutants of proteins involved in major peroxisomal pathways (photorespiration and fatty acids β-oxidation) but do not show obvious changes in transcript levels under drought, i.e. hpr1, gox1, acx3, acx6, and acx1 acx5. Mutants of mildly regulated genes, such as peroxisomal NAD+ transporter PXN and beta-hydroxyisobutyryl-CoA hydrolase CHY1 were assayed as well (Table C in S1 File).

For an efficient and quantitative drought tolerance assessment, we used the photosynthetic efficiency Fv/Fm as a drought susceptibility indicator in our screen. Each mutant was grown in the same pot with the wild type plant for 3.5 weeks with periodic irrigation, followed by an 18-day drought period, at the end of which Fv/Fm was measured. The positive control, ABA biosynthetic mutant aba1 [46], and mutants for the LON2 protease and photorespiratory enzyme HPR1, showed statistically significant decrease in Fv/Fm after drought treatment (Fig 3).

thumbnail
Fig 3. Fv/Fm comparison between the selected peroxisomal mutants and wild type plants grown in the same pot.

The aba1 mutant (salk_059469) was used as a positive control. Two biological replicates were used for each genotype. Asterisk indicates p-value < 0.01 in Student’s t test.

https://doi.org/10.1371/journal.pone.0137762.g003

The lon2 and hpr1 mutants were further analyzed to assess defects in drought resistance. Prior to drought treatment, lon2 and hpr1 mutants exhibited similar Fv/Fm values to that of the wild type (Fig 4A and 4B), suggesting that the drought sensitive phenotypes we observed were specific to drought stress and not a result from general growth defect. Compared with watered plants and drought-treated wild type plants, drought-treated lon2-2, hpr1-1 and hpr1-2 displayed an age-dependent gradient of photosynthetic defect, which was stronger in older leaves and milder in young leaves (Fig 4C and 4D). These mutants also showed defects in other drought stress indicators, including reduced anthocyanin induction (Fig 4F), accelerated chlorophyll degradation (Fig 4G), and lower relative water content (Fig 4H).

thumbnail
Fig 4. Drought resistance phenotypes of lon2 and hpr1 mutants.

(A) Images of plants (left) and color-coded chlorophyll fluorescence that indicates Fv/Fm values (right). (B) Fv/Fm comparison between mutants and wild type. Four biological replicates of each genotype were used. No significant difference in Fv/Fm was observed. (C) Plant images and color-coded chlorophyll fluorescence images that indicate Fv/Fm values. (D-E) Fv/Fm comparison between mutants and wild type. Two biological replicates were used for each genotype under each condition. Asterisk, p < 0.01 in Students’ t test. (F-H) Quantification of anthocyanin (F), chlorophyll (G), and relative water content (H) in mutants and wild type plants. Three biological replicates were used for each genotype under each condition. Asterisk indicates p < 0.01 in Students’ t test; N.D, not detectable; D, drought; W, watered.

https://doi.org/10.1371/journal.pone.0137762.g004

Discussion

We have constructed peroxisome-centered transcriptomic heatmaps using Arabidopsis microarray data from various developmental stages and under biotic and abiotic stresses. Results obtained from the in silico analysis not only showed correlation between protein function and expression regulation for many proteins with known function, but also provided information from which previously unknown roles may be inferred for some peroxisomal proteins. For example, out of the two peroxisomal MDH isoforms, MDH2 is tightly co-expressed with HPR1 and GOX1/2, suggesting that MDH2 is the major MDH isoform that functions in photorespiration. Among all the genes analyzed, AtHsp15.7, shows the strongest up-regulation by high light, suggesting that the small heat shock protein Hsp15.7, which was shown to be a stress-inducible constituent of the peroxisome [47], may facilitate the re-folding of proteins that have been partially unfolded or damaged under high light stress. Up-regulation of peroxisome elongation factors such as PEX11b, PEX11c and PEX11d under hypoxial and biotic stresses suggests that increased volume of peroxisome might be a mechanism for the plant to deal with enhanced oxidative stress. Finally, the fact that unknown protein UP6, which was previously shown to play a minor role in β–oxidation [45,48,49], is strongly up-regulated during late seed developmental stages, indicates a possible role for this protein in seed maturation.

Following the transcriptome analysis, we used a drought stress assay to test promising gene candidates, and identified LON2 and HPR1 as contributors to drought tolerance. HPR1 converts hydroxypyruvate to glycerate during photorespiration [2]. It is possible that drought induces stomatal closure, which limits the atmospheric uptake of CO2, thus activating the oxygenase activity of RuBisCO and subsequently, photorespiration. In hpr1, the accumulated photorespiratory metabolites may inhibit RuBisCO activity and slow down the Calvin-Benson cycle due to decreased supply of glycerate, thus leading to accumulated NADPH and ROS that cause a series of oxidative damages. Mutants for the other two genes directly (GOX1) or indirectly (CAT2) involved in photorespiration did not show a drought phenotype, possibly due to their functional redundancy with GOX2 and CAT1/CAT3, respectively.

LON2 is a protease with unknown substrates and a role in peroxisomal matrix protein import and degradation [50,51]. Our transcriptomic analysis found LON2 to be up-regulated by ABA by 4-fold, which is consistent with its 8-fold induction by ABA in guard cells [52], suggesting that LON2 may play a role in drought response through ABA signaling and peroxisomal protein quality control pathways. Since peroxisomal degradation via autophagy was shown to be enhanced in the lon2 mutant, especially in older leaves [50,53], it is possible that there are insufficient peroxisomes in the lon2 mutant to carry out photorespiration, which is critical for plant survival under drought conditions. This may also explain why lon2 and hpr1 display stronger phenotypes than the ABA biosynthetic mutant aba1 in the initial Fv/Fm screen, because our screen measured photosynthetic efficiency, which is directly impacted by photorespiration deficiencies in these two peroxisomal mutants.

Although in silico analysis is powerful for function prediction, mutants for many genes whose transcript levels are regulated by drought did not exhibit obvious drought tolerance defects. Given the manageable size of the peroxisomal proteome and available mutants, stress-based mutant screens should be a more direct way to identify peroxisomal proteins involved in stress response. In this initial screen, we identified strong drought sensitive phenotypes in the knockout mutants of the LON2 protease and the photorespiratory enzyme HPR1, suggesting that future larger-scale screens would be promising to investigate the role of peroxisomes in plant adaptation to environmental stresses comprehensively. The next step would be to link these peroxisomal proteins with the global stress response networks.

Supporting Information

S1 File. Supporting information.

This file includes the following: Figure A. Expression of the small heat shock protein gene AtHsp15.7 in response to abiotic stresses; Figure B. Total number of Arabidopsis peroxisomal genes with significantly changed expression levels in response to stresses; Table A. Microarray datasets used in this study; Table B. Arabidopsis peroxisomal gene list; Table C. Mutants used in the primary screen for drought tolerance.

https://doi.org/10.1371/journal.pone.0137762.s001

(PDF)

Acknowledgments

We would like to thank David Hall for assistance with Fv/Fm measurement, Dr. Jin Chen and Sahra Uygun for help with microarray data analysis, the Arabidopsis Biological Resource Center for providing the T-DNA insertion mutant lines, Dr. David Kramer for sharing the gox1 mutant and Dr. Gregg Howe for sharing the acx1 acx5 double mutant. No conflict of interest declared.

Author Contributions

Conceived and designed the experiments: JL JH. Performed the experiments: JL. Analyzed the data: JL JH. Contributed reagents/materials/analysis tools: JL JH. Wrote the paper: JL JH.

References

  1. 1. Pieuchot L, Jedd G (2012) Peroxisome assembly and functional diversity in eukaryotic microorganisms. Annu Rev Microbiol 66: 237–263. pmid:22994494
  2. 2. Hu J, Baker A, Bartel B, Linka N, Mullen RT, Reumann S, et al. (2012) Plant peroxisomes: biogenesis and function. Plant Cell 24: 2279–2303. pmid:22669882
  3. 3. Smith JJ, Aitchison JD (2013) Peroxisomes take shape. Nat Rev Mol Cell Biol 14: 803–817. pmid:24263361
  4. 4. Kaur N, Reumann S, Hu J (2009) Peroxisome biogenesis and function. The Arabidopsis Book 7: e0123. pmid:22303249
  5. 5. Kaur N, Hu J (2011) Defining the plant peroxisomal proteome: from Arabidopsis to rice. Front Plant Sci 2: 103. pmid:22645559
  6. 6. Del Rio LA (2011) Peroxisomes as a cellular source of reactive nitrogen species signal molecules. Arch Biochem Biophys 506: 1–11. pmid:21055384
  7. 7. Kotchoni SO, Gachomo EW (2006) The reactive oxygen species network pathways:an essential prerequisite for perception of pathogen attack and the acquired disease resistance in plants. J Biosci 31: 389–404. pmid:17006022
  8. 8. Willekens H, Chamnongpol S, Davey M, Schraudner M, Langebartels C, Van Montagu M, et al. (1997) Catalase is a sink for H2O2 and is indispensable for stress defence in C3 plants. Embo J 16: 4806–4816. pmid:9305623
  9. 9. Queval G, Issakidis-Bourguet E, Hoeberichts FA, Vandorpe M, Gakiere B, Vanacker H, et al. (2007) Conditional oxidative stress responses in the Arabidopsis photorespiratory mutant cat2 demonstrate that redox state is a key modulator of daylength-dependent gene expression, and define photoperiod as a crucial factor in the regulation of H2O2-induced cell death. Plant J 52: 640–657. pmid:17877712
  10. 10. Okinaka Y, Yang CH, Herman E, Kinney A, Keen NT (2002) The P34 syringolide elicitor receptor interacts with a soybean photorespiration enzyme, NADH-dependent hydroxypyruvate reductase. Mol Plant Microbe Interact 15: 1213–1218. pmid:12481993
  11. 11. Taler D, Galperin M, Benjamin I, Cohen Y, Kenigsbuch D (2004) Plant eR genes that encode photorespiratory enzymes confer resistance against disease. Plant Cell 16: 172–184. pmid:14688292
  12. 12. Rojas CM, Senthil-Kumar M, Wang K, Ryu CM, Kaundal A, Mysore KS. (2012) Glycolate oxidase modulates reactive oxygen species-mediated signal transduction during nonhost resistance in Nicotiana benthamiana and Arabidopsis. Plant Cell 24: 336–352. pmid:22286136
  13. 13. Dammann C, Ichida A, Hong B, Romanowsky SM, Hrabak EM, Harmon AC, et al. (2003) Subcellular targeting of nine calcium-dependent protein kinase isoforms from Arabidopsis. Plant Physiol 132: 1840–1848. pmid:12913141
  14. 14. Coca M, San Segundo B (2010) AtCPK1 calcium-dependent protein kinase mediates pathogen resistance in Arabidopsis. Plant J 63: 526–540. pmid:20497373
  15. 15. Lipka V, Dittgen J, Bednarek P, Bhat R, Wiermer M, Stein M, et al. (2005) Pre- and postinvasion defenses both contribute to nonhost resistance in Arabidopsis. Science 310: 1180–1183. pmid:16293760
  16. 16. Clay NK, Adio AM, Denoux C, Jander G, Ausubel FM (2009) Glucosinolate metabolites required for an Arabidopsis innate immune response. Science 323: 95–101. pmid:19095898
  17. 17. Maeda K, Houjyou Y, Komatsu T, Hori H, Kodaira T, Ishikawa A. (2009) AGB1 and PMR5 contribute to PEN2-mediated preinvasion resistance to Magnaporthe oryzae in Arabidopsis thaliana. Mol Plant Microbe Interact 22: 1331–1340. pmid:19810803
  18. 18. Bednarek P, Pislewska-Bednarek M, Svatos A, Schneider B, Doubsky J, Mansurova M, et al. (2009) A glucosinolate metabolism pathway in living plant cells mediates broad-spectrum antifungal defense. Science 323: 101–106. pmid:19095900
  19. 19. Westphal L, Scheel D, Rosahl S (2008) The coi1-16 mutant harbors a second site mutation rendering PEN2 nonfunctional. Plant Cell 20: 824–826. pmid:18430804
  20. 20. Spoel SH, Dong X (2012) How do plants achieve immunity? Defense without specialized immune cells. Nat Rev Immunol 12: 89–100. pmid:22273771
  21. 21. Reumann S (2004) Specification of the peroxisome targeting signals type 1 and type 2 of plant peroxisomes by bioinformatics analyses. Plant Physiol 135: 783–800. pmid:15208424
  22. 22. Rochon D, Singh B, Reade R, Theilmann J, Ghoshal K, Alam SB, et al. (2014) The p33 auxiliary replicase protein of Cucumber necrosis virus targets peroxisomes and infection induces de novo peroxisome formation from the endoplasmic reticulum. Virology 452: 133–142. pmid:24606690
  23. 23. McCartney AW, Greenwood JS, Fabian MR, White KA, Mullen RT (2005) Localization of the tomato bushy stunt virus replication protein p33 reveals a peroxisome-to-endoplasmic reticulum sorting pathway. Plant Cell 17: 3513–3531. pmid:16284309
  24. 24. Schmid M, Davison TS, Henz SR, Pape UJ, Demar M, Vingron M, et al. (2005) A gene expression map of Arabidopsis thaliana development. Nat Genet 37: 501–506. pmid:15806101
  25. 25. Quan S, Yang P, Cassin-Ross G, Kaur N, Switzenberg R, Aung K, et al. (2013) Proteome analysis of peroxisomes from etiolated Arabidopsis seedlings identifies a peroxisomal protease involved in beta-oxidation and development. Plant Physiol 163: 1518–1538. pmid:24130194
  26. 26. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. (2004) Bioconductor: open software development for computational biology and bioinformatics. Genome Biol 5: R80. pmid:15461798
  27. 27. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, et al. (2003) Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 4: 249–264. pmid:12925520
  28. 28. Smyth GK (2004) Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol 3: Article3. pmid:16646809
  29. 29. Warnes GR, Bolker B, Bonebakker L, Gentleman R, Liaw WHA, Lumley T, et al. (2015) gplots: Various R Programming Tools for Plotting Data. http://cran.r-project.org/web/packages/gplots/index.html.
  30. 30. Attaran E, Major IT, Cruz JA, Rosa BA, Koo AJ, Chen J, et al. (2014) Temporal dynamics of growth and photosynthesis suppression in response to jasmonate signaling. Plant Physiol 165: 1302–1314. pmid:24820026
  31. 31. Schneider CA, Rasband WS, Eliceiri KW (2012) NIH Image to ImageJ: 25 years of image analysis. Nat Methods 9: 671–675. pmid:22930834
  32. 32. Pandey N, Ranjan A, Pant P, Tripathi RK, Ateek F, Pandey HP, et al. (2013) CAMTA 1 regulates drought responses in Arabidopsis thaliana. BMC Genomics 14: 216. pmid:23547968
  33. 33. Pracharoenwattana I, Smith SM (2008) When is a peroxisome not a peroxisome? Trends Plant Sci 13: 522–525. pmid:18768343
  34. 34. Bauwe H, Hagemann M, Fernie AR (2010) Photorespiration: players, partners and origin. Trends Plant Sci 15: 330–336. pmid:20403720
  35. 35. Gehan MA, Greenham K, Mockler TC, McClung CR (2015) Transcriptional networks—crops, clocks, and abiotic stress. Curr Opin Plant Biol 101016/jpbi200908001 24: 39–46. pmid:25646668
  36. 36. Yamaguchi-Shinozaki K, Shinozaki K (2006) Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses. Annu Rev Plant Biol 57: 781–803. pmid:16669782
  37. 37. Sinclair AM, Trobacher CP, Mathur N, Greenwood JS, Mathur J (2009) Peroxule extension over ER-defined paths constitutes a rapid subcellular response to hydroxyl stress. Plant J 59: 231–242. pmid:19292761
  38. 38. Kaur N, Reumann S, Hu J (2009) Peroxisome biogenesis and function. Arabidopsis Book 7: e0123. pmid:22303249
  39. 39. Browse J (2009) Jasmonate passes muster: a receptor and targets for the defense hormone. Annu Rev Plant Biol 60: 183–205. pmid:19025383
  40. 40. Sorhagen K, Laxa M, Peterhansel C, Reumann S (2013) The emerging role of photorespiration and non-photorespiratory peroxisomal metabolism in pathogen defence. Plant Biology 15: 723–736. pmid:23506300
  41. 41. Kaur N, Li J, Hu J (2013) Peroxisomes and photomorphogenesis. Subcell Biochem 69: 195–211. pmid:23821150
  42. 42. Desai M, Hu J (2008) Light induces peroxisome proliferation in Arabidopsis seedlings through the photoreceptor phytochrome A, the transcription factor HY5 HOMOLOG, and the peroxisomal protein PEROXIN11b. Plant Physiol 146: 1117–1127. pmid:18203870
  43. 43. Zhu JK (2002) Salt and drought stress signal transduction in plants. Annu Revi Plant Biol 53: 247–273.
  44. 44. Cassin-Ross G, Hu J (2014) A simple assay to identify peroxisomal proteins involved in 12-oxo-phytodienoic acid metabolism. Plant Signal Behav 9: e29464.
  45. 45. Cassin-Ross G, Hu J (2014) Systematic phenotypic screen of Arabidopsis peroxisomal mutants identifies proteins involved in beta-oxidation. Plant Physiol 166: 1546–1559. pmid:25253886
  46. 46. Xiong L, Lee H, Ishitani M, Zhu JK (2002) Regulation of osmotic stress-responsive gene expression by the LOS6/ABA1 locus in Arabidopsis. J Biol Chem 277: 8588–8596. pmid:11779861
  47. 47. Ma C, Haslbeck M, Babujee L, Jahn O, Reumann S (2006) Identification and characterization of a stress-inducible and a constitutive small heat-shock protein targeted to the matrix of plant peroxisomes. Plant Physiol 141: 47–60. pmid:16531488
  48. 48. Quan S, Switzenberg R, Reumann S, Hu J (2010) In vivo subcellular targeting analysis validates a novel peroxisome targeting signal type 2 and the peroxisomal localization of two proteins with putative functions in defense in Arabidopsis. Plant Signal Behav 5: 151–153. pmid:20009535
  49. 49. Reumann S, Quan S, Aung K, Yang P, Manandhar-Shrestha K, Holbrook D, et al. (2009) In-depth proteome analysis of Arabidopsis leaf peroxisomes combined with in vivo subcellular targeting verification indicates novel metabolic and regulatory functions of peroxisomes. Plant Physiol 150: 125–143. pmid:19329564
  50. 50. Farmer LM, Rinaldi MA, Young PG, Danan CH, Burkhart SE, Bartel B. (2013) Disrupting autophagy restores peroxisome function to an Arabidopsis lon2 mutant and reveals a role for the LON2 protease in peroxisomal matrix protein degradation. Plant Cell 25: 4085–4100. pmid:24179123
  51. 51. Lingard MJ, Bartel B (2009) Arabidopsis LON2 is necessary for peroxisomal function and sustained matrix protein import. Plant Physiol 151: 1354–1365. pmid:19748917
  52. 52. Leonhardt N, Kwak JM, Robert N, Waner D, Leonhardt G, Schroeder JI. (2004) Microarray expression analyses of Arabidopsis guard cells and isolation of a recessive abscisic acid hypersensitive protein phosphatase 2C mutant. Plant Cell 16: 596–615. pmid:14973164
  53. 53. Goto-Yamada S, Mano S, Nakamori C, Kondo M, Yamawaki R, Kato A, et al. (2014) Chaperone and protease functions of LON protease 2 modulate the peroxisomal transition and degradation with autophagy. Plant Cell Physiol 55: 482–496. pmid:24492254