Figures
Abstract
Background
Herbivore grazing is a multiple-component process that includes wounding, defoliation, and saliva deposition. Despite the extensive published research on mechanical wounding and defoliation, no analysis to identify the genes that specify defoliation and mechanical wounding has been performed. Moreover, the influence of the expression of these genes on plant regrowth after defoliation remains poorly understood.
Results
Seven cDNA libraries for RNA samples collected from stubble tissues that had been mechanically wounded or defoliated at 2, 6 and 24 h along with the control were sequenced using the Illumina/Solexa platform. A comparative transcriptomic analysis of the sequencing data was conducted. In total, 1,836 and 3,238 genes were detected with significant differential expression levels after wounding and defoliation, respectively, during one day. GO, KOG and pathway-based enrichment analyses were performed to determine and further understand the biological functions of those differentially expressed genes (DEGs). The results demonstrated that both wounding and defoliation activated the systemic synthesis of jasmonate (JA). However, defoliation specifically reduced the expression levels of ribosomal protein genes, cell division or cell expansion-related genes, and lignin biosynthesis genes and may have negatively affected plant growth. Further analysis revealed that the regrowth of elongating leaves was significantly retarded after defoliation at 6 h through the following 7 days of measurement, suggesting that the gene expression pattern and phenotype are consistent. Fifteen genes were selected, and their expression levels were confirmed by quantitative RT-PCR (qRT-PCR). Thirteen of them exhibited expression patterns consistent with the digital gene expression (DGE) data.
Conclusions
These sequencing datasets allowed us to elucidate the common and distinct mechanisms of plant responses to defoliation and wounding. Additionally, the distinct DEGs represent a valuable resource for novel gene discovery that may improve plant resistance to defoliation from various processes.
Citation: Chen S, Cai Y, Zhang L, Yan X, Cheng L, Qi D, et al. (2014) Transcriptome Analysis Reveals Common and Distinct Mechanisms for Sheepgrass (Leymus chinensis) Responses to Defoliation Compared to Mechanical Wounding. PLoS ONE 9(2): e89495. https://doi.org/10.1371/journal.pone.0089495
Editor: Christian Schönbach, Nazarbayev University, Kazakhstan
Received: November 20, 2013; Accepted: January 21, 2014; Published: February 21, 2014
Copyright: © 2014 Chen et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by National Basic Research Program of China (“973”, 2014CB138704), the National High Technology Research and Development Program of China (“863”, 2011AA100209), the Project of Ningxia Agricultural Comprehensive Development Office (NTKJ-2013-03 (1), NTKJ-2013-04), and the Knowledge Innovation Program of the Chinese Academy of Sciences (KSCX2-EW-B-05). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Grasslands are a major part of the global ecosystem, covering 37% of the Earth's terrestrial area [1]. Herbivory is a fundamental driver of grassland plant community composition. Nearly all grasslands worldwide have been grazed by populations of wild and domestic herbivores, such as cattle [2], horses [3], sheep [4], deer [5], rabbits [6], and kangaroos [7]. Herbivore grazing is a multiple-component process that includes wounding, defoliation, and saliva deposition [8].
Wounding responses have been demonstrated in species throughout the plant kingdom, and previous research has focused on identifying systemic wound signals and the mechanisms by which they are generated, transported, and perceived [9]. Leaf wounding is a well-characterized stimulus for the rapid activation of jasmonate (JA) synthesis through distinct autonomous and non-autonomous cellular pathways [9]–[12].
Defoliation not only induces serious wounding but also removes a large fraction of photosynthetically active leaves from grass plants. Therefore, plant photosynthesis is reduced; carbohydrate reserves stored in the plant are mobilized; and the transcript levels of photosynthesis-related and carbohydrate metabolism-related genes, such as genes encoding the small subunits of Rubisco, Rubisco activase, light harvesting chlorophyll a/b binding proteins, fructan exohydrolase (FEH), sucrose transporter (SUT), invertase, β-amylase, and starch synthase, are changed [13]–[18].
We previously identified 466 genes that responded to simulated grazing using the Affymetrix Rice GeneChips. These genes encode proteins that participate in signal transduction, miRNA regulation, cell wall modification, carbohydrate metabolism, hormone synthesis, and molecular transport [8]. However, it is unclear whether the responses were due to wounding, defoliation, or other factors.
Currently, the genes that specify defoliation and mechanical wounding remain unknown, and little is understood about how expression of these genes influences plant regrowth after defoliation. In this study, the transcriptomes of sheepgrass stubble tissues collected from defoliated and mechanically wounded plants were compared using the Solexa sequencing system. By comparing the differentially expressed gene (DEG) responses to defoliation and wounding, the common and distinct pathways and genes between defoliation and wounding were identified. These sequencing datasets allowed us to elucidate the common and distinct mechanisms of the plant response to defoliation and wounding, and the distinct DEGs represent a valuable resource for novel gene discovery that may improve plant tolerance toward defoliation caused by herbivores, insects, pathogens, or mechanical stress.
Materials and Methods
Plant Materials, Growing Conditions, and Treatments
The sheepgrass variety Zhongke 3 was grown in a soil mix of peat moss and vermiculite (2∶1, v/v) in a greenhouse at 23°C with a photoperiod of 16 h light/8 h dark. Wounding and defoliation treatments were performed on eight-week-old sheepgrass seedlings. To induce defoliation, approximately two-thirds of the aboveground biomass was removed, and to induce wounding, tweezers were used to mechanically wound a site at the same position as the defoliation site. The remaining one-third of the aboveground biomass was collected 2, 6, and 24 h after defoliation, and the corresponding parts were collected for control seedlings and wounding at 2, 6, and 24 h. Each sample included 3 replicate pots and 21 seedlings. The samples were collected and immediately frozen in liquid nitrogen and stored at −80°C.
RNA Extraction and Solexa/Illumina Sequencing
Total RNA was extracted from individual control, wounded and defoliated plants at 2, 6, and 24 h using TRIzol reagent (Invitrogen, Carlsbad, CA, USA). After precipitation, the RNA was purified using Qiagen’s RNeasy Kit with on-column DNase I digestion according to the manufacturer’s instructions. RNA concentrations were measured using the Bioanalyzer 2100 (Agilent Technologies, Inc., Waldbronn, Germany), and the RNA integrity was analyzed on a 1.0% (w/v) agarose gel. Approximately 5 to 8 µg total RNA was used to construct each RNA-seq library. Poly (A) mRNAs were isolated from the total RNA using oligo (dT) magnetic beads (Illumina, San Diego, CA). RNA fragmentation, cDNA synthesis, and PCR amplification were performed according to the Illumina RNA-seq protocol (Cat # RS-100-0801). Seven cDNA libraries, which were isolated from mechanically wounded and defoliated material at 2, 6, and 24 h and the control, were sequenced using the Illumina HiSeq 2000 System at the Chinese National Human Genome Center (Shanghai, China). The sequence reads generated in this study have been deposited in the NCBI sequence read archive (SRA065691).
Analysis and Screening of Differentially Expressed Genes (DEGs)
The reads for each sample were counted and mapped back to previously generated sheepgrass reference genes [19] using the software tool Bowtie [20]. The read number of each gene was transformed into RPKM (Reads Per Kilobases per Million reads) [21], and the differently expressed genes were identified by a DEGseq package using the MARS method (MA-plot-based method with Random Sampling model) [22]. The p-value was obtained from the differential gene expression test. FDR (False Discovery Rate) manipulation was used to determine the p-value threshold in multiple tests and analyses. Both FDR ≤0.001 and the absolute value of | the log2Ratio ≥1 were used as thresholds to identify the significant DEGs [23].
Gene Ontology, EuKaryotic Orthologous Groups and Pathway Enrichment Analysis
Gene Ontology (GO) is an international standardized functional gene classification system that describes the properties of genes and their products in any organism using three ontologies: cellular component, molecular function, and biological process. GO enrichment analysis was performed using agriGO (http://bioinfo.cau.edu.cn/agriGO/, [24]). EuKaryotic Orthologous Groups (KOG) enrichment analysis was conducted through hypergeometric distribution testing using the Phyper function in the R software package (http://www.rproject.org/). The Bonferroni correction was used to adjust the p-values. The significantly enriched functional clusters were selected based on a corrected q-value less than 0.05. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was conducted using KOBAS 2.0 (http://kobas.cbi.pku.edu.cn/, [25]).
Real-time Quantitative RT-PCR (qRT-PCR) Analysis
Real-time quantitative RT-PCR (qRT-PCR) analysis was used to verify the DGE results. The RNA samples used for the qRT-PCR assays were the same as in the DGE experiments. Gene-specific primers were designed according to the reference unigene sequences using Primer Premier 6.0. Fifteen genes were selected from the DEGs under defoliation at 2 h for quantitative RT-PCR assays, and these genes were involved in JA biosynthesis, the gibberellin signaling pathway, cell wall biogenesis/degradation, and starch and sucrose metabolism (Table S1). qRT-PCR was performed according to the manufacturer’s specifications (TaKaRa SYBR Premix Ex Taq II Kit, Dalian, China). The following SYBR Green PCR cycling conditions were used: denaturation at 95°C for 10 s, followed by 40 cycles of 94°C for 5 s and 60°C for 20 s. The PCR experiments were performed using an iQ 5 Multicolor Real-time PCR Detection System (Bio-Rad, USA). The results were normalized to the sheepgrass actin gene, and the relative gene expression levels were calculated using 2−ΔΔCT.
Results
Wounding and Defoliation Exhibit Distinct DEGs
To identify the DEGs involved in wounding and defoliation, we sequenced seven cDNA libraries and generated 192 million clean reads ranged from 4,811,093 to 6,411,192 per library after filtering the low-quality and adapter sequences. These clean reads were mapped to the sheepgrass reference genes as described in the Materials and methods. In total, 36,413, 35,925, 35,113, 35,929, 34,534, 36,828, and 36,956 unigenes were detected for sequenced libraries 1–7, respectively (Table S2). The reference unigenes mapped by all clean reads were enhanced as the sequencing amount (total read number) increased. However, when the sequencing counts reached 1.5 million reads or higher, the number of detected genes became saturated (Figure S1), indicating that the sequencing depth was sufficient to cover the transcriptome.
By comparing the treated and control Solexa libraries, a great number of differentially expressed transcripts were identified. To study the subset of genes associated with wounding and defoliation, we analyzed the most differentially regulated genes with a |log2Ratio|≥1 and a FDR ≤0.001. There were 706 (356 up-regulated and 350 down-regulated), 891 (442 up-regulated and 449 down-regulated), and 803 (362 up-regulated and 441 down-regulated) genes that responded to wounding at 2, 6, and 24 h, respectively. For defoliation, we identified 2,211 (685 up-regulated and 1526 down-regulated), 997 (451 up-regulated and 546 down-regulated), and 1,111 (480 up-regulated and 631 down-regulated) genes that were significantly expressed at 2, 6, and 24 h, respectively (Figure 1A). These significantly differentially wounding- and defoliation-responsive genes were further divided into genes expressed at one time point and those expressed at overlapping time points. For example, 1,836 and 3,238 genes were differentially expressed between the wounding and defoliation treatments, respectively, during one day (Figure 1B, C), and 954 genes were uniquely expressed under only defoliation. These DEGs are listed in Table S3 and Table S4, respectively.
The categories of DEGs expressed at one time point and at overlapping time points for wounding and defoliation at 2, 6, and 24(A) A summary of the up- and down-regulated DEG numbers. (B) and (C) DEGs expressed at one time point and overlapping time points for wounding and defoliation, respectively.
GO Functional Analysis of DEGs
A total of 835 DEGs that responded to wounding and 1,560 DEGs that responded to defoliation were divided into 44 GO categories (Figure 2). In the biological process category, cellular processes (wounding 66.6%; defoliation 68.5%) and metabolic processes (wounding 61.5%; defoliation 62.5%) were the most dominant groups, followed by stimulus responses (wounding 23.9%; defoliation 19.4%) and cellular component organization (both approximately 17%). Regarding the molecular function category, 65.9% of the wounding unigenes and 63.5% of the defoliation unigenes were assigned to binding, followed by catalytic activity (both approximately 49%), transcription regulatory activity (4.6% for both), and structural molecule activity (wounding 4.0%; defoliation 7.5%). In the cellular component category, cell and cell part genes (wounding 72.1% and defoliation 75.5% for both) dominated, followed by organelle (wounding 49.9%; defoliation 52.4%) and organelle part genes (wounding 25.5%; defoliation 26.4%).
The frequency of GO terms was analyzed using the GO Slim Assignment. The x- and y-axes represent the names of clusters and the ratio of each cluster, respectively. W, mechanical wounding; D, defoliation.
To reveal significantly enriched GO terms in DEGs compared to the reference transcriptome, GO enrichment analysis of the functional significance was performed. The GO terms with p≤0.05 were considered significantly enriched among the DEGs. This analysis allowed us to determine the major molecular functions, biological processes, and cellular components that were represented among the DEGs (Figure S2). For the molecular functions, the significantly enriched GO terms of the DEGs for wounding and defoliation included DNA binding, serine-type endopeptidase inhibitor activity, oxidoreductase activity, and structural ribosome constituents specific to defoliation. For the enriched biological processes, there were 53 GO terms of DEGs for wounding and 40 GO terms of DEGs for defoliation, which represented cellular component organization and cellular component biogenesis (e.g., chromatin assembly, protein-DNA complex assembly, and nucleosome assembly), metabolic processes (fatty acid and oxylipin biosynthesis), lipid localization, and responses to stimuli (e.g. response to water and defense response to insects). For the cellular components, significantly enriched GO terms of DEGs for wounding and defoliation included chromatin, nucleosome, nucleus, extracellular region, and ribosome, and cell wall was specific to defoliation.
KOG and KEGG Pathway Enrichment Analysis of DEGs
A total of 832 wounding-responsive DEGs and 1,605 defoliation-responsive DEGs were classified functionally into 24 clusters of KOG. KOG enrichment analyses were performed by hypergeometric distribution testing, and the Bonferroni correction was applied. Significantly enriched functional clusters were selected based on a q-value ≤0.05. The results demonstrated that genes associated with chromatin structure and dynamics, energy production and conversion, and amino acid transport and metabolism were over-represented in the KOG classification of wounding (Table 1), and chromatin structure and dynamics, carbohydrate transport and metabolism, translation, ribosomal structure and biogenesis, and the cytoskeleton were over-represented in the KOG classification of defoliation (Table 2).
A total of 460 wound-responsive DEGs and 881 defoliation-responsive DEGs were classified into the KEGG pathway. For wounding, the significantly enriched KEGG pathways included valine, leucine, and isoleucine degradation; linoleic acid metabolism; ascorbate and aldarate metabolism; starch and sucrose metabolism; fatty acid elongation; arginine and proline metabolism; and alpha-linolenic acid metabolism (Table 3). For defoliation, starch and sucrose metabolism; ribosome; alpha-linolenic acid metabolism; phenylpropanoid biosynthesis; flavonoid biosynthesis; stilbenoid, diarylheptanoid and gingerol biosynthesis; linoleic acid metabolism; and phenylalanine metabolism were over-represented in the KEGG pathway classification (Table 4).
Enrichment of KEGG Pathways Common to Wounding and Defoliation
From the above KEGG pathway enrichment analysis of the DEGs, we observed that linoleic acid metabolism, alpha-linolenic acid metabolism, and starch and sucrose metabolism were enriched KEGG pathways common to both wounding and defoliation. In the alpha-linolenic acid metabolism and linoleic acid metabolism pathways, we identified genes encoding lipoxygenase (LOX), allene oxide synthase (AOS), allene oxide cyclase (AOC), 12-oxophytodienoic acid reductase, acyl-CoA oxidase (ACOX1, ACOX3), acetyl-CoA acyltransferase 1, and enoyl-CoA hydratase/3-hydroxyacyl-CoA dehydrogenase, all enzymes involved in JA biosynthesis. The gene expression levels were significantly increased for wounding and defoliation at 2 h (Figure 3). In starch and sucrose metabolism, genes involved in sucrose metabolism, including genes encoding sucrose synthase, fructokinase, and beta-fructofuranosidase (cell wall invertase), were down-regulated, while genes related to starch biosynthesis and metabolism and trehalose biosynthesis, such as genes encoding starch synthase, starch phosphorylase, beta-amylase, and trehalose 6-phosphate synthase, were up-regulated (Figure 3).
Green to bright red: down-regulation to up-regulation. W, mechanical wounding; D, defoliation.
Defoliation-enriched KEGG Pathways
For defoliation, the specifically enriched KEGG pathways were ribosome; phenylpropanoid biosynthesis; flavonoid biosynthesis; stilbenoid, diarylheptanoid and gingerol biosynthesis; and phenylalanine metabolism. Most ribosomal genes were down-regulated after defoliation at 2 h, but some were up-regulated at 24 h (Figure 4A). The expression levels of the genes involved in phenylpropanoid biosynthesis; flavonoid biosynthesis; stilbenoid, diarylheptanoid and gingerol biosynthesis; and phenylalanine metabolism were predominantly decreased (Figure 4B); these down-regulated genes included those encoding cinnamyl-alcohol dehydrogenase, peroxidase, caffeoyl-CoA O-methyltransferase, beta-glucosidase, phenylalanine/tyrosine ammonia-lyase, shikimate O-hydroxycinnamoyltransferase, chalcone synthase, anthocyanidin reductase, and flavonoid 3′,5′-hydroxylase. However, some genes involved in these pathways, including the genes encoding trans-cinnamate 4-monooxygenase, coumaroylquinate (coumaroylshikimate) 3′-monooxygenase, phenylalanine ammonia-lyase, and primary-amine oxidase, were up-regulated (Figure 4B).
(A) Ribosomal DEGs. (B) DEGs involved in phenylpropanoid biosynthesis; flavonoid biosynthesis; stilbenoid, diarylheptanoid and gingerol biosynthesis; and phenylalanine metabolism. Green to bright red: down-regulation to up-regulation. D, defoliation.
Confirmation of Read-mapped Genes by qRT-PCR
To confirm the reliability of the Solexa/Illumina sequencing technology, fifteen genes were selected as described in the Materials and methods for quantitative RT-PCR assays. The expression patterns for thirteen genes were consistent between the qRT-PCR and the DGE analyses (Figure 5). The corresponding primers are listed in Table S1. The inconsistent expression between the qPCR and DGE analyses for the two remaining genes was most likely because DGE was more sensitive in detecting low-abundance transcripts and small changes in gene expression than qRT-PCR.
Effect of Defoliation on Regrowth of Elongating Leaves
Eight-week-old plants were defoliated to a stubble height of 5 cm. The regrowth of the elongating leaves emerging above the cutting level was measured (Figure 6). The regrowth of elongating leaves was significantly restrained at 6 h after defoliation compared with the growth of the control plants, and the inhibition continued through the following seven days of measurement.
After 8 weeks of growth, sheepgrass plants were defoliated to a stubble height of 5 cm, and regrowth of elongating leaves was measured above the cutting level after defoliation at different time points.
Discussion
Distinct Global Gene Transcription Changes for Sheepgrass under Wounding and Defoliation Treatments
In this study, a transcriptome profiling analysis was performed to identify genes that are differentially expressed in wounding- and defoliation-treated plants. A sequencing depth of 4.8 to 6.4 million reads per library was reached (Table S2). The results of the saturation sequencing analyzed in the seven libraries indicated that the sequencing depth was sufficient to cover the transcriptome (Figure S1). We identified 1,836 and 3,238 wounding-responsive and defoliation-responsive DEGs, respectively, which suggested that defoliation led to nearly twice as many DEGs as wounding. Among these DEGs, 954 were uniquely expressed under defoliation and thus represent candidate genes that merit further study. Further analysis indicated that the defoliation-responsive gene number was highest at 2 h, during which the number of down-regulated defoliation-responsive DEGs (1,526) were over four times the number of down-regulated wounding-responsive DEGs (350), suggesting that plants could rapidly respond to the enormous damage of defoliation by regulating global gene expression.
Wounding and Defoliation Commonly Activate Systemic JA Synthesis
JA is derived from the unsaturated fatty acid α-linolenic acid (18∶3), an abundant octadecanoid in the cellular membranes of higher plants. Lipoxygenase (LOX) is the initial step in JA biosynthesis and catalyzes the oxygenation of α-linolenic acid, leading to the formation of hydroperoxy octadecadienoic acids, which are then converted into 12-oxo phytodienoic acid (12-OPDA) in the plastids via enzymatic reactions catalyzed first by allene oxide synthase (AOS) and subsequently by allene oxide cyclase (AOC). OPDA is then converted to JA by 12-oxo-phytodienoic acid reductase (OPR) followed by three cycles of β-oxidation in the peroxisome [26], [27]. In this study, we observed that the DEGs involved in JA biosynthesis after wounding and defoliation were significantly enriched according to the GO terms and KEGG pathway enrichment analysis (Figure S2, Table 3, 4). The transcription levels of key genes in JA biosynthesis, including genes encoding lipoxygenase (LOX), allene oxide synthase (AOS), allene oxide cyclase (AOC), 12-oxophytodienoic acid reductase (OPR), acyl-CoA oxidase (ACOX1, ACOX3), acetyl-CoA acyltransferase 1, and enoyl-CoA hydratase/3-hydroxyacyl-CoA dehydrogenase (Figure 3), were significantly increased after wounding and defoliation. Our results are consistent with the previous reports, demonstrating that wound-induced responses are triggered by the de novo synthesis of the plant hormone JA when plant tissues are injured by various means, including herbivores, pathogens, and mechanical stress [9], [28]–[29].
Defoliation Specifically Reduces the Expression Levels of Ribosomal Protein Genes, Cell Division or Cell Expansion-related Genes, and Lignin Biosynthesis Genes and may Negatively Affect Plant Growth
Enrichment analyses of GO terms, KOG, and the KEGG pathway suggested that the ribosome was specifically over-represented after defoliation. The ribosome is an essential ribonucleoprotein complex that is engaged in translation and is indispensable for growth [30]. In this study, many ribosomal protein genes, including the large 60S and small 40S subunits of the cytosolic ribosome, were identified, and their expression levels mainly decreased at 2 h after defoliation (Figure 4A). In both plants and animals, the loss of ribosomal proteins leads to reduced growth, which most likely correlated with reduced ribosome production and lower rates of protein synthesis [31]. For example, in Arabidopsis thaliana, several RP loss-of-function mutations that affect cell division or cell expansion and consequently result in a deformed leaf size and shape have been identified, indicating cell- or development-specific roles of RPs during leaf growth [32]. Notably, in this study, the expression levels of numerous genes involved in cell division or cell expansion were down-regulated after defoliation, and these genes include those encoding cyclins, cyclin-dependent kinases, and expansion (Table S4).
Lignin is a phenolic heteropolymer of the secondary cell walls that plays a major role in plant development and defense against pathogens. The monolignols represent the main component of lignin, and their synthesis involves many intermediates and enzymes. In this study, we identified several key genes involved in lignin biosynthesis that were down-regulated, including genes encoding phenylalanine ammonia-lyase (PAL), caffeoyl-CoA O-methyltransferase (CCOMT), shikimate O-hydroxycinnamoyltransferase (HCT), and cinnamyl alcohol dehydrogenase (CAD) (Figure 4B). The down-regulation of these genes has been demonstrated to reduce lignin biosynthesis and cell-wall thickness, alter xylem organization, and retard growth [33]–[36].
Notably, the regrowth of elongating leaves was significantly retarded after defoliation at 6 h through the following 7 days of measurement (Figure 6), suggesting that the above gene expression patterns correlated consistently with the phenotype. Castrillón-Arbeláez et al. demonstrated that although defoliation led to a rapid and transient reduction of non-structural carbohydrates (NSC) in grain amaranth, only a few changes in gene expression and enzyme activity could be associated with the NSC changes in an analysis of 25 genes. They also revealed that the rapid mobilization of foliar starch reserves followed by an efficient recovery of all NSC reserves after defoliation did not alter plant growth or reproductive fitness [17]. Vargas-Ortiz et al. demonstrated that defoliation reduced all NSC levels in the stems and roots of grain amaranth, and this reduction was associated with reduced sucrose synthase and cell wall invertase activity [37]. However, in this study, we observed that the expression levels of genes encoding sucrose synthase and cell wall invertase were reduced under both mechanical wounding and defoliation. Ida et al. also suggested that although defoliation in Oxytropis sericea reduced photosynthesis and nectar production, it did not alter photosynthate allocation or fruit or seed production [38]. However, Akiyama and Ågren demonstrated that defoliation in Arabidopsis thaliana reduced seed production, and this reduction correlated negatively with the removed leaf area [39]. The above contradictory results suggest that defoliation is a complex process, which includes defoliation timing, intensity, and species tolerance, and may lead to different responses. The highly complex nature of the NSC metabolic flux in plants also indicates that defoliation effects cannot be fully explained at the transcriptional level, and its effects on plant growth and reproductive fitness might be elucidated by systems biology in the future.
Supporting Information
Figure S1.
Sequencing saturation analysis of the seven libraries. The number of detected genes was enhanced as the sequencing amount (total read number) increased. W, mechanical wounding; D, defoliation.
https://doi.org/10.1371/journal.pone.0089495.s001
(TIF)
Figure S2.
Gene Ontology term enrichment analysis of wounding- and defoliation-responsive DEGs. (A–C) and (D–F) present the GO terms ‘enrichment status’ and ‘hierarchy’ for wounding and defoliation, respectively; (A and D) the biological process, (B and E) molecular function, and (C and F) cellular component branches. The classification terms and their serial numbers are represented as boxes, and the box includes the GO term, the adjusted p-value (in parentheses), the item number mapping the GO term in the query list and background, and the total number of items in the query list and background. The boxes with significant levels are indicated by color. The color scale from light to dark illustrates the p-value cutoff levels from low to high.
https://doi.org/10.1371/journal.pone.0089495.s002
(DOC)
Table S2.
Read number and mapping results for the seven independent libraries.
https://doi.org/10.1371/journal.pone.0089495.s004
(DOC)
Table S3.
DEGs expressed after mechanical wounding.
https://doi.org/10.1371/journal.pone.0089495.s005
(XLS)
Author Contributions
Conceived and designed the experiments: SC GL. Performed the experiments: SC YC LZ DQ XL. Analyzed the data: SC LZ XY LC QZ. Contributed reagents/materials/analysis tools: GL SC XL. Wrote the paper: SC.
References
- 1. O'Mara FP (2012) The role of grasslands in food security and climate change. Ann Bot 110(6): 1263–1270.
- 2. Roche LM, Latimer AM, Eastburn DJ, Tate KW (2012) Cattle grazing and conservation of a meadow-dependent amphibian species in the Sierra Nevada. PLoS One 7(4): e35734.
- 3. de Villalobos AE, Zalba SM, Peláez DV (2011) Pinus halepensis invasion in mountain pampean grassland: effects of feral horses grazing on seedling establishment. Environ Res 111(7): 953–959.
- 4. Scohier A, Dumont B (2012) How do sheep affect plant communities and arthropod populations in temperate grasslands? Animal 6(7): 1129–1138.
- 5. Lapointe L, Bussières J, Crête M, Ouellet JP (2010) Impact of growth form and carbohydrate reserves on tolerance to simulated deer herbivory and subsequent recovery in Liliaceae. Am J Bot 97(6): 913–924.
- 6. Rueda M, Rebollo S, García-Salgado G (2013) Contrasting impacts of different-sized herbivores on species richness of Mediterranean annual pastures differing in primary productivity. Oecologia 172(2): 449–459.
- 7. Meers T, Adams R (2003) The impact of grazing by Eastern Grey Kangaroos (Macropus giganteus) on vegetation recovery after fire at Reef Hills Regional Park, Victoria. Ecol Manage Rest 4: 126–132.
- 8. Chen S, Li XQ, Zhao A, Wang L, Li X, et al. (2009) Genes and pathways induced in early response to defoliation in rice seedlings. Curr Issues Mol Biol 11(2): 81–100.
- 9. Koo AJ, Howe GA (2009) The wound hormone jasmonate. Phytochemistry 70(13–14): 1571–1580.
- 10. Chung HS, Koo AJK, Gao X, Jayanty S, Thines B, et al. (2008) Regulation and function of Arabidopsis JASMONATE ZIM-domain genes in response to wounding and herbivory. Plant Physiol 146: 952–964.
- 11. Glauser G, Grata E, Dubugnon L, Rudaz S, Farmer EE, et al. (2008) Spatial and temporal dynamics of jasmonate synthesis and accumulation in Arabidopsis in response to wounding. J Biol Chem 283: 16400–16407.
- 12. Strassner J, Schaller F, Frick U, Howe GA, Weiler EW, et al. (2002) Characterization and cDNA microarray expression analysis of 12-oxophytodienoate reductases reveals differential roles for octadecanoid biosynthesis in local versus systemic wound responses. Plant J 32: 585–601.
- 13. Lee JM, Sathish P, Donaghy DJ, Roche JR (2011) Impact of defoliation severity on photosynthesis, carbon metabolism and transport gene expression in perennial ryegrass. Funct Plant Biol 38: 808–817.
- 14. Tamura KI, Sanada Y, Tase K, Komatsu T, Yoshida M (2011) Pp6-FEH1 encodes an enzyme for degradation of highly polymerized levan and is transcriptionally induced by defoliation in timothy (Phleum pratense L.). J Exp Bot 62(10): 3421–3431.
- 15. Berthier A, Desclos M, Amiard V, Morvan-Bertrand A, Demmig-Adams B, et al. (2009) Activation of sucrose transport in defoliated Lolium perenne L.: an example of apoplastic phloem loading plasticity. Plant Cell Physiol 50(7): 1329–1344.
- 16. Furet PM, Berthier A, Decau ML, Morvan-Bertrand A, Prud'homme MP, et al. (2012) Differential regulation of two sucrose transporters by defoliation and light conditions in perennial ryegrass. Plant Physiol Biochem 61: 88–96.
- 17. Castrillón-Arbeláez PA, Martínez-Gallardo N, Arnaut HA, Tiessen A, Délano-Frier JP (2012) Metabolic and enzymatic changes associated with carbon mobilization, utilization and replenishment triggered in grain amaranth (Amaranthus cruentus) in response to partial defoliation by mechanical injury or insect herbivory. BMC Plant Biol 12: 163.
- 18. Lothier J, Lasseur B, Prud’homme MP, Morvan-Bertrand A (2010) Hexokinase-dependent sugar signaling represses fructan exohydrolase activity in Lolium perenne. Funct Plant Biol 37: 1151–1160.
- 19. Chen S, Huang X, Yan X, Liang Y, Wang Y, et al. (2013) Transcriptome analysis in sheepgrass (Leymus chinensis): a dominant perennial grass of the Eurasian Steppe. PLoS One 8(7): e67974.
- 20. Langmead B, Trapnell C, Pop M, Salzberg SL (2009) Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 10(3): R25.
- 21. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B (2008) Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 5(7): 621–628.
- 22. Wang L, Feng Z, Wang X, Wang X, Zhang X (2010) DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 26(1): 136–138.
- 23. Benjamini Y, Yekutieli D (2001) The control of the false discovery rate in multiple testing under dependency. Ann Stat 29: 1165–1188.
- 24. Du Z, Zhou X, Ling Y, Zhang Z, Su Z (2010) agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res 38: W64–70.
- 25. Xie C, Mao X, Huang J, Ding Y, Wu J, et al. (2011) KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res 39: W316–322.
- 26. Lyons R, Manners JM, Kazan K (2013) Jasmonate biosynthesis and signaling in monocots: a comparative overview. Plant Cell Rep 32(6): 815–827.
- 27. Wasternack C, Hause B (2013) Jasmonates: biosynthesis, perception, signal transduction and action in plant stress response, growth and development. An update to the 2007 review in Annals of Botany. Ann Bot 111(6): 1021–1058.
- 28. Zhang Y, Turner JG (2008) Wound-induced endogenous jasmonates stunt plant growth by inhibiting mitosis. PLoS One 3(11): e3699.
- 29. Koo AJ, Gao X, Jones AD, Howe GA (2009) A rapid wound signal activates the systemic synthesis of bioactive jasmonates in Arabidopsis. Plant J 59(6): 974–986.
- 30. Horiguchi G, Van Lijsebettens M, Candela H, Micol JL, Tsukaya H (2012) Ribosomes and translation in plant developmental control. Plant Sci 191–192: 24–34.
- 31. Byrne ME (2009) A role for the ribosome in development. Trends Plant Sci 14(9): 512–529.
- 32. Schippers JHM, Mueller-Roeber B (2010) Ribosomal composition and control of leaf development. Plant Sci 179: 307–315.
- 33. Barakat A, Bagniewska-Zadworna A, Choi A, Plakkat U, DiLoreto DS, et al. (2009) The cinnamyl alcohol dehydrogenase gene family in Populus: phylogeny, organization, and expression. BMC Plant Biol 9: 26.
- 34. Guo D, Chen F, Inoue K, Blount JW, Dixon RA (2001) Downregulation of caffeic acid 3-O-methyltransferase and caffeoyl CoA 3-O-methyltransferase in transgenic alfalfa. impacts on lignin structure and implications for the biosynthesis of G and S lignin. Plant Cell 13(1): 73–88.
- 35. Day A, Neutelings G, Nolin F, Grec S, Habrant A, et al. (2009) Caffeoyl coenzyme A O-methyltransferase down-regulation is associated with modifications in lignin and cell-wall architecture in flax secondary xylem. Plant Physiol Biochem 47(1): 9–19.
- 36. Gallego-Giraldo L, Jikumaru Y, Kamiya Y, Tang Y, Dixon RA (2011) Selective lignin downregulation leads to constitutive defense response expression in alfalfa (Medicago sativa L.). New Phytol 190(3): 627–639.
- 37. Vargas-Ortiz E, Espitia-Rangel E, Tiessen A, Délano-Frier JP (2013) Grain amaranths are defoliation tolerant crop species capable of utilizing stem and root carbohydrate reserves to sustain vegetative and reproductive growth after leaf loss. PLoS One 8(7): e67879.
- 38. Ida TY, Harder LD, Kudo G (2012) Effects of defoliation and shading on the physiological cost of reproduction in silky locoweed Oxytropis sericea. Ann Bot 109(1): 237–246.
- 39. Akiyama R, Ågren J (2012) Magnitude and timing of leaf damage affect seed production in a natural population of Arabidopsis thaliana (Brassicaceae). PLoS One 7(1): e30015.