Figures
Abstract
Background
Switchgrass (Panicum virgatum L.) is being developed as a bioenergy crop for many temperate regions of the world. One way to increase biomass yields is to move southern adapted lowland cultivars to more northern latitudes. However, many southerly adapted switchgrass germplasm can suffer significant winter kill in northerly climes.
Materials and Methods
Here, we have applied next-generation sequencing in combination with biochemical analyses to query the metabolism of crowns and rhizomes obtained from two contrasting switchgrass cultivars. Crowns and rhizomes from field-grown lowland (cv Kanlow) and upland (cv Summer) switchgrass cultivars were collected from three randomly selected post-flowering plants. Summer plants were senescing, whereas Kanlow plants were not at this harvest date.
Results
Principal component analysis (PCA) differentiated between both the Summer and Kanlow transcriptomes and metabolomes. Significant differences in transcript abundances were detected for 8,050 genes, including transcription factors such as WRKYs and those associated with phenylpropanoid biosynthesis. Gene-set enrichment analyses showed that a number of pathways were differentially up-regulated in the two populations. For both populations, protein levels and enzyme activities agreed well with transcript abundances for genes involved in the phenylpropanoid pathway that were up-regulated in Kanlow crowns and rhizomes. The combination of these datasets suggests that dormancy-related mechanisms had been triggered in the crowns and rhizomes of the Summer plants, whereas the crowns and rhizomes of Kanlow plants had yet to enter dormancy.
Conclusions
Delayed establishment of dormancy at more northerly latitudes could be one factor that reduces winter-survival in the high-yielding Kanlow plants. Understanding the cellular signatures that accompany the transition to dormancy can be used in the future to select plants with improved winter hardiness.
Citation: Palmer NA, Saathoff AJ, Tobias CM, Twigg P, Xia Y, Vogel KP, et al. (2014) Contrasting Metabolism in Perenniating Structures of Upland and Lowland Switchgrass Plants Late in the Growing Season. PLoS ONE 9(8): e105138. https://doi.org/10.1371/journal.pone.0105138
Editor: Marie-Anne Van Sluys, Universidade de Sao Paulo, Brazil
Received: March 25, 2014; Accepted: July 17, 2014; Published: August 18, 2014
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the paper and its Supporting Information files.
Funding: This work was supported by the Office of Science (BER), U. S. Department of Energy Grant Number DE-AI02-09ER64829, and by the USDA-ARS CRIS project 5440-21000-030-00D. 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
Switchgrass (Panicum virgatum L.), a perennial C4-grass native to the continental USA, is being developed as a major biomass feedstock for use in temperate regions [1]. Latitudinal adaptation has resulted in ecotypes with distinguishable genetic make-up [2] and differential responses to photoperiod [3]. Moving switchgrass germplasm more than one hardiness zone north or south from its native adaptation zone or latitude generally results in significantly reduced biomass yields. This is caused by early flowering in the northern adapted plants when moved to southern latitudes and reduced winter survival of southern plants when grown in northern latitudes because of delayed flowering and senescence [3].
For switchgrass, perenniality and sustainability appear to be strongly interlinked and are significantly influenced by photoperiod [3], [4]. Photoperiodic cues likely influence the orderly developmental transition from spring emergence of shoots through the senescence of the aerial parts of the plant and the imposition of dormancy in the below-ground parts of the plant post seed ripening [5]. As with other perennial warm-season grasses, the above ground parts of the switchgrass plant senesce at the end of the growing season, whereas the below-ground tissues, (comprised of the crowns, rhizomes and roots), transition to a dormant state. Tiller buds originating from the crowns and rhizomes are the source of the eventual regrowth of the above ground portions of the plant in the following growing season. Little is currently known at a molecular or cellular level on these aspects of switchgrass biology; however, if these processes can be understood it will be possible to improve both breeding and management of this crop. This knowledge can be used not only to improve the sustainability of switchgrass and other temperate perennial biomass grasses, but also potential perennial grain crops [6] which will face similar trade-offs of maximizing grain yield while maintaining crown and rhizome vitality. Seasonal senescence and dormancy has been studied in a number of other species [7]–[10], and these studies offer a scaffold to interpret data obtained from understudied species such as switchgrass.
This study focused on two tetraploid cultivars, namely, “Summer” that is adapted to the northern Great Plains of the USA and a contrasting cultivar “Kanlow” that is adapted to the southern latitudes of the USA. The cultivar Summer was developed at Brookings, South Dakota (∼44°N), USA from a native collection that originated from southeast Nebraska, USA and was first released in 1963 [11]. The cultivar Kanlow was also first released in 1963 from Manhattan, Kansas (∼39°N), USA and was developed from 200 plants originally sourced from a lowland site near Wetumka, Oklahoma, USA. Kanlow plants display adaptation to wet sites [11]. There is about a 26 cm differential in average rainfall between SE Nebraska (origin of cv Summer, drier) and Wetumka, OK (origin of Kanlow, wetter). Although we cannot rule out environmental-specific transcriptional differences between the two populations, our intention was to focus on understanding any stage specific differences. Little, if any, selection pressure was applied during the generation of the cultivars, and as with most switchgrass cultivars, both Summer and Kanlow are synthetic populations that display diversity for many plant characteristics including height, leaf morphology, phenology and genotype as has been described in other publications [1]–[3], [5], [12]. Tetraploid switchgrasses exhibit higher yields than octaploid switchgrasses, and are providing germplasm for the development of improved plant materials for the biofuel sector [13], [14]. Within the tetraploid switchgrasses, the lowland cultivars such as Kanlow out-yield upland tetraploids such as Summer, but suffer from significant winter-kill in more northerly latitudes limiting their potential to be grown in these regions [1], [3]. The aims of this study were to utilize RNA-seq with some biochemical analyses to develop knowledge into the metabolism of crown and rhizome tissues near the end of the growing season. The combination of these unbiased datasets obtained from the crowns and rhizomes from field-grown plants suggests that dormancy-related mechanisms had been triggered in the crowns and rhizomes of the Summer plants, whereas the crowns and rhizomes levels had yet to enter dormancy. Recently, RNA-seq has been used to map transcripts and discover single nucleotide polymrophisms (SNPs) in several different populations of switchgrass [15]. In this study the authors were able to efficiently map a large number of sequencing reads and identify over one million SNPs. A similar approach has been used in our study to understand the transcriptomes in the two contrasting switchgrass cultivars.
Materials and Methods
Plant Materials, Growth Conditions and Selection of Harvest Date
Stands of switchgrass cv Summer and cv Kanlow were established in the fields of the University of Nebraska-Agricultural Research Division near Mead, NE (∼41°N) in 2009, using seedlings raised in the greenhouse and transplanted to the fields. Each plot (1 m×1.2 m) contained 12 closely-spaced plants to mimic sward densities. Plots were in a field that measured approximately 35 m×30 m [16]. Agronomic management was as described earlier [17]. Plants were not irrigated and depended on ambient rainfall and stored soil moisture for growth. Climatic conditions for a few weeks prior to, and after collection are shown in Figure S1 in File S1. Daytime temperatures were generally greater than 20°C for the ∼4 weeks prior to harvest and low temperatures were generally 10° to 15°C. Three continuous days of rainfall were recorded three days prior to harvest. Soils were moist at harvest (Figure S1 in File S1).
At this location, Summer plants reach 50% anthesis generally by early July and reach physiological maturity by early-mid September. In contrast, for Kanlow plants, 50% anthesis is reached by mid-August and physiological maturity in late September/October. We did not control for plant size (although Kanlow plants were always larger than Summer plants), or measure above-ground yields. Individual plants planted in 2009 were still distinct in 2010, and the crowns and rhizomes were collected from three individual plants of each cultivar as described earlier [16] on September 29, 2010. Data were not collected on the above ground material except to visually verify seed developmental stage and if there were obvious signs of plant senescence (leaf yellowing, tiller death, etc.). Such signs of senescence were present in Summer plants, but were not observed for the Kanlow plants. This developmental state was selected as a starting point to understand the changes in the below-ground metabolism of switchgrass that impacted transition into winter dormancy. At harvest, tissues were flash frozen in liquid nitrogen and later pulverized cryogenically in order to obtain high yields of total RNA [16].
RNA-Seq
Total RNA was extracted from tissues with TRIzol Reagent (Invitrogen Corp, Carlsbad, CA), followed by further purification by Qiagen RNeasy Mini Kit (Qiagen, Valencia, CA). Sequencing libraries were prepared from poly(A)+RNA using the Illumina mRNA-Seq Sample Prep Kit according to the manufacturer's instructions. Libraries were then sequenced on the Illumina Genome Analyzer IIx using two 36-cycle sequencing kits to read 75 nucleotides of sequence from a single end of each insert, by the v8 protocol (Illumina Inc, San Diego, CA). Each sample was analyzed on a single lane. The purity, integrity, and profile of extracted total RNA was verified on an Agilent 2100 BioAnalyzer (Agilent Technologies, Inc., Santa Clara, USA), and for concentration using a NanoDrop instrument (Thermo Scientific, Waltham, USA). Once generated, quality of the libraries was verified by gel electrophoresis and on a NanoDrop spectrophotometer. The PhiX virus DNA library was used as a control in a sequencing lane on each flowcell run to monitor sequencing operation. Obtained sequences were trimmed for adapter sequences and verified for quality using the software available with the instrument.
Mapping and Differential gene expression analysis
Illumina reads were mapped to the draft switchgrass genome available at www.phytozome.net [18], using splice-junction aware Tophat2 (version 2.0.11) [19], with Bowtie2 (version 2.2.3) [20] alignment. Default parameters were used and reads flagged as having multiple alignments were not included in gene expression. Gene expression counts were calculated using featureCounts [21], a script included as part of the Subread (version 1.4.4) [22] analysis package, and the gene annotation file provided with the draft genome release. Differential expression analysis was conducted using the program edgeR (version 3.4.2) [23] in R [24]. Very low expression genes were filtered out of the data set prior to analysis by requiring 1 count per million (cpm) in at least three of the six samples. An FDR threshold of 0.05 was used to determine differentially expressed genes. Variance stabilizing transformation values, calculated using the DESeq2 (version 1.2.8) [25] package in R, were used for gene expression quantities in further downstream analyses.
Genome Functional Annotation
Draft genome nucleotide transcript sequences were functionally annotated by sequence similarity using the program Blast2GO [26] similarly to our earlier work [16]. Blast2GO annotations with an e-value threshold of 1×10−15 or lower resulting from a Blastx search against the NCBI non-redundant protein database were used in all the subsequent analyses, including Enzyme Commission (EC) (http://www.chem.qmul.ac.uk/iubmb/enzyme/) and KEGG pathways [27], [28].
Gene Set Enrichment Analysis
Gene Set Enrichment Analysis was done using GSEA from the Broad Institute [29]. Gene sets were created from Blast2GO generated KEGG pathway assignments.
Metabolite Analyses
Aliquots (100±2 mg) of ground crown and rhizome tissue were extracted for metabolites, generally following a previously described protocol [30]. Metabolites were separated and identified (where possible) by GCMS [31] using the Agilent Fiehn GC/MS Metabolomics RTL library as the basis for instrument operation and metabolite identification. Metabolite identification was accomplished using the Fiehn libraries in Chemstation and in the supplied version of AMDIS. Metabolites were extracted from 100 mg of tissues with 350 µL 100% methanol and 20 µL of 0.2 mg mL−1 ribitol in water, as an internal standard. Samples were then heated in a 70°C heat block for 15 minutes. After heating, 1 volume (370 µL) DI water (18 MΩ) was added and each tube was vigorously mixed and then centrifuged at 14,000 RPM for 10 min. The supernatants were placed into new 1.5 mL tubes, and polar and nonpolar metabolites separated by two chloroform washes (300 µL each time). The final upper phase (methanol and water) was carefully pipetted into a new microfuge tube and 25 µL aliquots from each tube were then placed into 2 mL glass vials. Next, 5 µL of 0.4 mg mL−1 myristic-d27 acid (#366889 Sigma-Aldrich Co.) in 100% methanol and 20 µL of 0.2 mg mL−1 docosane (#134457 Sigma-Aldrich Co.) in dichloromethane were added to each vial; these were used as a retention locking compound and internal standard, respectively. The vials were evaporated to dryness under vacuum (Labconco Inc.). The dried samples were derivatized using 30 µL of 40 mg mL−1 methoxyamine hydrochloride at 37°C for three hours followed by 30 min. of trimethylsilylation at 37°C using 90 µL N-methyl-N-(trimethylsilyl)trifluoroacetamide (Thermo Scientific Inc.). The derivatized sample was placed into a new 2 mL glass vial with an insert due to the small volumes. The GC column was a 30-meter, 0.250 mm I.D. HP-5MS column (Agilent Technologies, Inc.). The mass selective detector was auto-tuned regularly with perfluorotributylamine as recommended by the manufacturer and operated over a m/Q scan range of 50–600. Temperature settings were 250°C for the injector and MS source and 290°C for the transfer line. The initial oven temperature was 60°C, which was increased at a rate of 10°C per minute to 325°C. Helium was the carrier gas and the flow rate was 1.1 mL min−1. The injector was operated in splitless mode. Since only major metabolites were of interest, features that had a mean major ion count of less than 1.0×105 across both genotypes were discarded. Also, features that were missing in more than 10 of the 27 GC-MS runs per cultivar were flagged; features that were flagged this way across both cultivars were eliminated from the data set since they were not reliably detected across the biological and technical replicates. After this filtering, a total of 418 features remained in the data set, although not all of these features could be reliably identified by Chemstation or AMDIS.
Enzyme, Protein and Immunoassays
Total soluble protein was extracted by sonicating approximately 100±2 mg of pulverized crown and rhizome tissue with 600 µL of Tris-HCl buffer, pH 8.0 containing a plant protease cocktail (Sigma P9599) with 5 mM DTT for cinnamyl alcohol dehydrogenase (CAD) or with 300 µL of Tris-HCl buffer, pH 8.0, without DTT for ascorbate peroxidase (AscPx) and caffeoyl-O-methyl transferase (COMT). For phenylalanine ammonia lyase (PAL), tissues were extracted using 0.1 M borate buffer, pH 8.8. Crude homogenate was desalted using a spin column (Thermo Scientific 89889, 2 mL columns) and used directly for enzyme assays. Protein content was measured using a dye-binding assay (Pierce 660 nm Protein assay kit) using BSA as a standard. Gel electrophoresis, protein transfer and immunodetection were as described previously [32]. Chemiluminescence was detected using a BioRAD Chemidoc system. Enzyme activities were analyzed as described previously [32]–[34].
Statistical Analyses
Statistical analysis of metabolite data was conducted using SAS 9.3 (SAS Institute, Inc., Cary, NC) and JMP 9.0 (SAS Institute Inc., Cary, NC). Pairwise comparisons across the two cultivars used the Wilcoxon rank-sum test in PROC NPAR1WAY, and the raw p-values were adjusted in PROC MULTTEST to correct for multiple comparisons using the false discovery rate [35]. JMP 9.0 was used to conduct principal component analysis and generate the color map. ANOVA testing was performed on other datasets using statistical routines available in EXCEL.
Results
80% of all reads map to the switchgrass draft genome
The mapping summary for the three biological replicates obtained from field grown Summer and Kanlow plants are shown in Table S1 in File S1. Each sample yielded approximately 24–27 million reads, of which greater than 77–83% were mapped to the 0.0 early-release version of the switchgrass genome. Functional annotation of the PviDraft0 genome by Blast2GO resulted in the annotation of 45,487 gene models at an e-value threshold of 1×10−15 or lower (Table S2 in File S1).
Transcriptomes of Summer and Kanlow crowns and rhizomes are significantly different
EdgeR analysis identified a total of 8,050 genes that were differentially expressed (File S2) between Kanlow and Summer crown and rhizomes at an FDR<0.05. A principal component analysis (PCA) of this expression dataset effectively separated the three Kanlow samples from the three Summer samples (Figure 1A). Component 1 accounted for 43% of the variance and separated the two cultivars. Component 2 accounted for approximately 18% of the variance and separated the three Summer samples. A Venn diagram of the total number of genes found in the Summer and Kanlow datasets is shown in Figure 1B. In total transcripts were detected for 36,572 genes. Of this total ∼84% (30,771 genes) were common to both populations and 2,303 and 3,498 genes were unique to the Summer and Kanlow datasets respectively (Figure 1B). About 41% of the differentially expressed genes (DEGs) had two-fold or greater average expression in Kanlow, and ∼26% had two-fold or greater expression in Summer (Figure 1B).
(A) PCA of transcriptomes of the three biological replicates obtained from each ecotype. Kanlow (orange triangle) were differentiated from Summer (cyan squares) by the first component. (B) Venn diagram of numbers of common and highly expressed genes in Summer (log2FC<-1) and Kanlow (log2FC>1), including DEGs in the Kanlow (orange) and Summer (cyan) datasets. (C) Mean counts for twelve transcripts identified as WRKY transcription factors. Stars above orange bars (Kanlow) identify putative switchgrass orthologs of WRKY factors involved in defense response in other systems. The inverted triangle over a cyan bar identifies a WRKY ortholog to an Arabidopsis gene involved in senescence that is upregulated in Summer crowns and rhizomes.
A more detailed evaluation of the DEGs was performed to identify changes in transcript abundance for three plant-specific classes of transcription factors, namely MYBs, NACs and WRKYs which are known to have multiple functions in cellular processes [36]–[40]. A scan of the switchgrass genome based on Pfam PF03106 (WRKY), PF00249 (MYB) and PF02365 (NAC) domains [41], identified 178 WRKYs (∼89 WRKYs per diploid genome); 324 MYBs (∼162 MYBs per diploid genome); and 239 NACs (∼120 NACs per diploid genome). Similar searches for WRKY members in other species gave the following results: Arabidopsis – 72; rice – 94 and 82 in Brachypodium. Using a similar approach for MYBs and NACs in other species, we found: Arabidopsis – 256/112 (MYBs/NACs); rice – 233/140 (MYBs/NACs), and Brachypodium – 200/100 (MYBs/NACs). Based on the relative numbers of these genes found in model plants, we expect that approximately 80–90% of the orthologous switchgrass genes were identified by this search. A total of 199 MYBs (30 DEG) and 108 NACs (20 DEG) were found in the expression datasets. In all 119 WRKYs (∼67% of the all WRKY genes identified in the genome) were found in the Summer and Kanlow datasets. Of these 119 genes, a total of 12 were differentially expressed at an FDR<0.05 (Figure 1D). A database search indicated that orthologs to four of the WRKY genes upregulated in Kanlow (identified by a star symbol in Figure 1C) were involved in responses to pathogens. In contrast, the ortholog to one WRKY gene upregulated in Summer and has been demonstrated to be an important mediator of cellular responses to senescence (inverted triangle).
Metabolite analysis differentiates Kanlow from Summer crowns and rhizomes
Principal component analysis of a broad-based polar metabolite profiling of Kanlow and Summer crown and rhizome extracts using GC-MS revealed that Summer and Kanlow extracts were clearly differentiated by the first principal component, which explained 27% of the variance (Figure 2A). The second component differentiated among the individual genotypes within each cultivar and explained an additional 13% of the variance. Concentration ratios, approximated by major ion count, were coupled with differences that had an FDR<0.05 to generate a volcano plot of the metabolite data (Figure 2B). Concentration ratios between the cultivars were considered meaningful if the averaged ion counts for a metabolite had a ratio greater than 1.75 or less than 1.75−1. The stricter criteria resulted in 66 metabolites with sufficient FDR values and ion area ratios in Summer and 47 metabolites in Kanlow. Hierarchical two-way clustering was used to generate an initial color map of the overall metabolite profile (Figure S2 in File S1). Following this analysis, we generated the heat map shown in Figure 2C by using one-way hierarchical clustering based on cultivar, and manually reordered metabolites into more biologically related groupings which could be indicative of differences in their metabolism (Figure 2C A–G). In Summer plants, levels of several amino acids (Figure 2C group A), amine derivatives (group B) and some sugar alcohols (group C) were significantly higher. In contrast, erythrose-4-phosphate, maltotriose, sucrose and trehalose were detected in greater levels in Kanlow crowns and rhizomes. Further significant differences (based on an FDR<0.05) in the metabolism between Kanlow and Summer rhizomes were evident in the observed levels of several organic acids such as glucuronic acid, a precursor of xylans [42], and a range of secondary metabolites, including phenolic acids (Figure 2C, group D,E). Summer tissues also appeared to be in a more oxidized state as compared to Kanlow crowns and rhizomes as they contained significantly higher levels of dehydroascorbic acid, gluconic acid lactone, and lower concentrations of α-tocopherol, (Figure 2C, group D, F).
(A) PCA of overall metabolite profiles observed by GCMS for for each of three individual genotypes within each cultivar. Kanlow tissue extracts (orange triangles) are separated from Summer tissue extracts (cyan squares) by the first component. Two-way error bars are shown for each plant that was based on nine separate GCMS runs for each sample. (B) Volcano plot showing the log10FDR versus the log1.75 fold change in peak area for major ions for all metabolites detected by GCMS. Significant differences were observed for several metabolites between the two cultivars (Kanlow orange triangles) and Summer (cyan squares). Grey circles are metabolites that failed to show sufficient differences in ion area of had an FDR value of >0.05. (C) Heat map shows marked differences in tissue abundances of selected metabolites in Kanlow and Summer crowns and rhizomes. Data are the average of triplicate injections from each of three separate extractions from the three biological replicates of each cultivar. Scale is from high abundance (yellow) to low abundance (black) for each metabolite. Data were subjected to one-way hierarchical clustering based on cultivar, and metabolites were manually reordered into more biologically related groupings (A–G).
Nexus between gene-set enrichment and metabolite analyses
As a prelude to gene-set enrichment analysis (GSEA), 124 gene sets (Table S3 in File S1) were generated for the predicted transcriptome available in the PviDraft0 using Blast2GO and the Kyoto Encyclopedia of Genes and Genomes (KEGG; http://www.genome.jp/kegg/). These gene sets were then used through GSEA to identify specific pathways differentially regulated in the crown and rhizome tissues of Kanlow and Summer plants. A total of 24 metabolic pathways (FDR<0.2) were different between the two ecotypes (Table S4 in File S1). In Kanlow tissues, there was an apparent upregulation of phenylpropanoid biosynthesis (cell wall accretion), starch, sucrose and fatty acid metabolism. Within Summer plants, there was enrichment (FDR<0.001) in transcripts for enzymes involved in diterpenoid biosynthesis, and a somewhat weaker upregulation for transcripts for enzymes associated with the biosynthesis of unsaturated fatty acids (FDR<0.193). A number of metabolic pathways associated with the degradation and scavenging of carbon compounds were also enriched in Summer crowns and rhizomes (Table S4 in File S1), providing some consistency with metabolite data.
In Kanlow tissues, there was significant enrichment in the transcripts coding for enzymes in the phenylpropanoid (lignin biosynthesis) pathway as compared to Summer tissues (Figure 3). Key metabolites caffeic acid and ferulic acid were enriched in Kanlow tissues. In Summer crowns and rhizomes, levels of Phe and Tyr, precursors for the phenylpropanoid pathway were significantly greater, consistent with the transcriptomic data.
Orange boxes identify individual enzymes upregulated in Kanlow. Boxes with EC numbers and relative enrichment of transcripts (Kanlow/Summer, K/S) are shown. Boxes with numbers indicate metabolites associated with this pathway identified by GCMS. Cyan (higher levels in Summer tissues), orange (higher in Kanlow tissues) and grey boxes (not significantly different).
Enzyme activity, protein levels and enrichment of key transcription factors confirm GSEA of phenylpropanoid pathway
The phenylpropanoid pathway was targeted to understand if the changes documented for transcript abundance (Figure 3) were also evident at the protein level. Transcripts for several genes associated with the phenylpropanoid pathway were significantly enhanced in Kanlow crowns and rhizomes as compared to the Summer plants. These included PAL, 4-coumarate-CoA ligase (4-CL), CAD, COMT and caffeoyl CoA 3-O-methyltransferase (CCoAOMT) (Figure 4A). It should be noted that transcript abundance for a called gene in Figure 4A is representative of all sequences assigned to this class of gene by Blast2GO, and is likely to contain transcripts for one or more related genes within a given gene family. Abundances varied over an order of magnitude, suggesting that differential regulatory mechanisms might exist for these individual genes. Transcript levels for actin were not significantly different (Figure 4A). AscPx transcripts were slightly elevated in Summer tissues as compared to Kanlow. AscPx is involved in detoxification of ROS in cells. Next, antibodies generated to select enzymes in the lignin biosynthesis pathway were utilized and showed that apparent protein levels for PAL, 4-CL, CAD, COMT and CCoAOMT were lower in all three biological replicates of Summer crown and rhizome extracts as compared to Kanlow extracts (Figure 4B). Levels for s-adenosyl methionine synthetase (SAMS), an enzyme that produces S-adenosyl methionine, which is a substrate for COMT and CCoAOMT, was also depressed in Summer crowns and rhizomes. Signal intensity for actin (used as a loading control) was similar for the Kanlow extracts and two of the Summer extracts, and diminished in one Summer extract. However, all the other proteins probed were essentially identical to the levels detected in the two other Summer plants (Figure 4B). The consistent differences observed in the transcript abundances and protein levels were also mirrored in the activities of PAL, CAD and COMT which were significantly lower or not detected in extracts from Summer tissues as compared to extracts prepared from the Kanlow crowns and rhizomes. AscPx activities were greater in Summer plants as compared to the Kanlow plants (Figure 4C) in agreement with the transcript and protein data. These data suggested that transcription factors known to affect the phenylpropanoid pathway [43] could be similarly regulated in these plants at this harvest date. The expected relationships of the key transcription factors involved in secondary cell wall deposition and lignification, adapted from [43], is shown in Figure 4D. Individual regulatory genes are color coded in boxes and the cellular pathways they impact are shown in ovals. Arrows connect the putative relationships between these gene products in model systems and to their respective cellular processes (ovals, Figure 4D). Although MYB 26 has been depicted as a master regulator in Figure 4D, its overall role is not as well established as for the other regulatory genes [43]. The XND1 gene (pink box) is a negative regulator of xylogenesis, and inhibits programmed cell death (PCD) and related processes in Arabidopsis [44]. The best switchgrass orthologs to these individual transcription factors (color coded identically in Figure 4D and 4E) which positively regulate secondary cell growth and/or lignification exhibited greater ratios of transcript abundance (Kanlow/Summer; K/S). Most of these genes were significantly upregulated (FDR<0.05) in Kanlow crowns and rhizomes (highlighted in orange). In contrast, transcripts for the XND1 ortholog were significantly greater in Summer tissues (Figure 4E).
(A) Transcript abundances for all transcripts identified as belonging to a specific gene by Blast2GO. PAL, phenylalanine ammonia lyase; 4CL, 4-coumarate-CoA ligase; CAD, cinnamyl alcohol dehydrogenase; COMT, caffeoyl-O-methyltransferase; CCoAOMT, caffeoyl-CoA 3O-methyltransferase; SAMS, s-adenosylmethionine synthetase; AscPX, ascorbate peroxidase. Orange bars = Kanlow; Cyan bars = Summer. (B) Immunblots of crown and rhizome extracts separated by SDS-PAGE. Extracts from three biological replicates for each cultivar are shown. These were the same tissues used for transcriptome and metabolite analyses. Other designations as described for panel A. (C) Enzyme activities of PAL, COMT, AscPx and CAD. Orange bars = Kanlow; Cyan bars = Summer. Other designations as described for panel A. (D) A map of potential relationships of key transcription factors that are known to play a major role in plant cell wall developmental processes. Each specific factor is shown in a different colored box. Figure is based on work described by Zhao and Dixon (2011). (E) Best switchgrass orthologs to transcription factors shown in (D). Boxes are colored identically to Panel D above. The switchgrass loci are identified and the FDR<0.05 values are shown in the last column. For the last column, boxes in orange color (transcripts significantly upregulated in Kanlow), cyan color (significantly upregulated in Summer), uncolored, not significant.
Discussion
In this study the dynamic aspects of the crown and rhizome transcriptome and metabolome in switchgrass has been explored for the first time. Analysis at this single harvest date revealed interesting insights into the metabolism of the crowns and rhizomes obtained from lowland Kanlow and upland Summer plants which differ in their photoperiod maturity response at the field site used for this experiment [3]. The work reported here complements several recent articles have expanded the available molecular resources for switchgrass and have begun the process of understanding switchgrass cellular processed within the context of its genomic complexity [16], [45]–[50].
The use of NGS coupled to bioinformatic analyses allowed approximately 80% of all reads to be effectively mapped to the switchgrass genome. Subsequent bioinformatic analyses indicated that these transcripts arose from 36,572 genes after filtering, which is ∼56% of the total loci identified so far in the switchgrass genome, suggesting that the mapping and identification processes were relatively robust. It is likely that potentially interesting features could have been missed due to the mapping of reads to the unfinished contig-based switchgrass genome assembly. Thus, mapping of transcripts to incomplete genes, gene sequences with no discernible domains or near the end of contigs would have been excluded from our dataset. However, these data provided adequate depth to begin an understanding of the metabolism of the crowns and rhizomes from the two contrasting switchgrass populations used in this study.
RNA-Seq of the transcriptomes provided a first broad-scale means to directly assess the metabolic status of these plants. PCA showed that the transcriptomes were clearly differentiated by the first component, indicating that between population (cultivar) differences were greater than within population differences. Both switchgrass cultivars are a heterogenous collection of genotypes [11] and the data for Summer plants would support this fact. The greater variation observed between the Summer plants could be indicative of differences in the timing of aerial senescence in these genotypes and its effect on crown and rhizome metabolism. In contrast, the transcriptomes of the Kanlow plants appeared to be more tightly clustered. Indeed the metabolomes also reflected this analysis. The two types of relatively non-biased large-scale approaches (RNA-Seq and GCMS) yielded similar conclusions with regard to the divergent metabolism in the Summer plants as compared to the Kanlow plants at this specific harvest date. Environmental conditions prior to the harvest date indicated that the field site had received adequate rainfall and had not experienced freezing temperatures minimizing abiotic stresses on these plants. Although these data are supportive of the notion that the differences in metabolism observed in the crowns and rhizomes were more reflective of developmental events, it would be difficult to rule out biotic and abiotic stresses these plants may have experienced during growth and during sample handling.
Differential enrichment of transcripts, metabolites and gene-sets were observed between the two populations. Analysis of transcription factors such as WRKYs [51]–[53] which impact an extensive array of plant cellular processes was undertaken. Approximately ∼67% of the all WRKYs tentatively identified in the switchgrass genome were represented in the total transcriptomic datasets, indicative of their potential role in switchgrass crowns and rhizomes. Within this larger group, of the nine WRKY genes significantly upregulated in Kanlow plants, four genes appeared to be orthologs to Arabidopsis WRKYs with a role in the defense response to pathogens [54]–[56]. In contrast, the WRKY factor enriched in Summer tissues have been implicated as a negative regulator of leaf senescence [55]. Although the actual roles of individual switchgrass WRKYs remain to be explored, they could be reflective of other changes observed at the metabolite and GSEA levels. Summer crowns and rhizomes were enriched in some oxidized compounds and appeared to contained lower levels of α-tocopherol [57], and phenolic acids that could potentially function as antioxidants [58], [59]. Although these data show a link between the transcriptomic and metabolomic datasets, much work will be needed to unequivocally solve these interrelationships.
Transcriptomic and metabolomic data were consistent with the hypothesis that Kanlow tissues were actively utilizing carbon in contrast to Summer crowns and rhizomes. Sucrose and trehalose levels as well the GSEA results showing the enrichment of sucrose-starch metabolism pathways in Kanlow crowns were consistent with this hypothesis. Also other downstream biosynthetic pathways that rely on carbon skeletons derived from sugars were also enriched in Kanlow tissues. Most notably, there was a strong enrichment in pathways linked to cell wall accretion and growth occurring in Kanlow plants relative to the Summer plants. Data shown in Figure 4 effectively linked transcript and metabolite datasets to the underlying protein activities, at least for a portion of the lignin biosynthesis pathway which showed higher activity in Kanlow tissues. Certain phenolic metabolites, namely cinnamic acid, ferulic acid and caffeic acid that are associated with monolignol synthesis [31], [60], were also more enriched in Kanlow consistent with enhanced flow of carbon through this pathway. In contrast, the amino acid precursors for this pathway, Phe and Tyr, were significantly enriched in Summer crowns and rhizomes which suggested that down-regulation of the phenylpropanoid pathway enzymes could have led to the increased accumulation of these amino acids. The significant up-regulation in transcript abundance for the switchgrass genes identified as potential orthologs of master transcription factors that regulate secondary cell wall biogenesis and lignification, based on [43], in Kanlow plants was striking and validated the other protein and metabolite findings. In contrast, transcripts for the putative switchgrass NAC gene XND1 were significantly greater in Summer tissues, suggestive of reduction in xylogenesis, and possibly in growth of tiller initials and rhizomes. In Arabidopsis, XND1 acts as a repressor of xylogenesis and programmed cell death (PCD) [44].
A number of KEGG pathways were found in Summer crowns and rhizomes that were enriched for transcripts, although not all enzymes were populated, unlike the more clear-cut observations on the phenylpropanoid pathway in Kanlow tissues. In trying to establish the possible connections between these diverse metabolic pathways, it became apparent that most possessed a potential common metabolite in acetyl-CoA (Figure 5). Potentially, scavenging of carbon compounds to form acetyl-CoA could represent a means to meet energy demands needed for maintaining cellular health in the possible absence of photosynthate delivery from the senescent shoots in Summer plants, which are consistent with results obtained for leafy spurge [61]. Acetyl-CoA serves as a key link between cellular metabolism, energetics and the acetylation of proteins in animal cells [62]. Similar mechanisms appear to be operating in plants [63]–[65], suggesting that a deeper understanding of the flux of carbon through a acetyl-CoA hub could provide new insights into the metabolism of crowns and rhizomes accompanying the seasonal growth habit of switchgrass and related perennial plants. Unfortunately, the metabolite analysis protocol used in this study cannot detect acetyl-CoA. Future use of other platforms such as LC-ms/ms might yield information on these more unstable metabolites.
Cyan squares are KEGG pathways found to be upregulated in Summer plants by GSEA, and cyan circles are, and metabolites that were elevated in Summer plants relative to Kanlow plants. Edges connecting pathways (squares) to metabolites (circle) indicate that the given metabolite is found in the connected pathway. Pathways with direct connections in KEGG are also connected by edges. Acetyl-CoA (diamond) could be a potential linker molecule among these diverse pathways and is suggested to be a possible metabolic hub in Summer crowns and rhizomes entering dormancy.
Conclusions
Data presented in this manuscript showed that differences existed in the metabolite and transcript profiles in the crowns and rhizomes obtained from field-grown switchgrass plants belonging to two populations with divergent photoperiod responses. Crowns and rhizomes obtained from the more southern adapted cultivar Kanlow appeared to be in a “growth-mode” with enrichment in gene-sets associated with biosynthesis and accretion of tissues. In contrast, the northern adapted cultivar Summer appeared to have entered a more quiescent state with greater enrichment of transcripts and metabolites favoring channeling of carbon skeletons through an acetyl-CoA hub. Our findings provide an initial framework to understand differences in crown and rhizome metabolism which could be exploited to enhance the latitudinal adaptation of diverse switchgrass populations.
Supporting Information
File S1.
Table S1 - Mapping summary for switchgrass crown and rhizome samples analyzed by HTS; Table S2 - Annotation of all transcripts using Blast2GO; Table S3 - List of 124 KEGG Pathways populated by Blast2GO analysis of the switchgrass draft genome; Table S4 - Gene Set Enrichment Analysis Results. KEGG pathways with an FDR<0.2 are shown. Transcript observed for individual enzymes in a given pathway are tabulated (enzymes); Figure S1 - Temperature and rainfall at the field site prior to and post-harvest; Figure S2 - Heat map of all major ions detected by GCMS using two-way hierarchical clustering. Yellow = high, Black = low. Data are for three independent extractions from three plants.
https://doi.org/10.1371/journal.pone.0105138.s001
(PDF)
File S2.
Excel file with list of differentially expressed genes (DEGs).
https://doi.org/10.1371/journal.pone.0105138.s002
(XLSX)
Acknowledgments
We thank Dr. Jean-Jack Reithoven and Mei Chen at the University of Nebraska-Lincoln Genomics Core Facility for preparing, sequencing and validating sequence data. We also thank Steve Masterson and Patrick Callahan for excellent technical support. The U.S. Department of Agriculture, Agricultural Research Service, is an equal opportunity/affirmative action employer and all agency services are available without discrimination. Mention of commercial products and organizations in this manuscript is solely to provide specific information. It does not constitute endorsement by USDA-ARS over other products and organizations not mentioned. All sequences have been deposited with the NCBI under the following BioSample accessions: SAMN01125827, SAMN01125828. Temporary SubmissionID: SUB120690. The PviDraft0 sequence data used in this study were produced by the US Department of Energy Joint Genome Institute.
Author Contributions
Conceived and designed the experiments: NP AS CT PT YX KV SM SS GS. Performed the experiments: NP AS CT PT YX GS. Analyzed the data: NP AS CT PT SM SS GS. Contributed reagents/materials/analysis tools: KV CT PT YX SM SS GS. Contributed to the writing of the manuscript: NP AS CT PT YX KV SM SS GS.
References
- 1. Vogel KP, Sarath G, Saathoff AJ, Mitchell RB (2011) Switchgrass. Rsc Energy Environ S 3: 341–380.
- 2. Okada M, Lanzatella C, Saha MC, Bouton J, Wu R, et al. (2010) Complete switchgrass genetic maps reveal subgenome collinearity, preferential pairing and multilocus interactions. Genetics 185: 745–760.
- 3. Casler MD, Vogel KP, Taliaferro CM, Ehlke NJ, Berdahl JD, et al. (2007) Latitudinal and longitudinal adaptation of switchgrass populations. Crop Sci 47: 2249–2260.
- 4. Sarath G, Baird LM, Mitchell RB (2014) Senescence, dormancy and tillering in perennial C4 grasses. Plant Sci 217–218: 140–151.
- 5. Moore KJ, Moser LE, Vogel KP, Waller SS, Johnson BE, et al. (1991) Describing and quantifying growth stages of perennial forage grasses. Agron J 83: 1073–1077.
- 6. Glover JD, Reganold JP, Bell LW, Borevitz J, Brummer EC, et al. (2010) Agriculture. Increased food and ecosystem security via perennial grains. Science 328: 1638–1639.
- 7. Hemming MN, Trevaskis B (2011) Make hay when the sun shines: the role of MADS-box genes in temperature-dependant seasonal flowering responses. Plant Sci 180: 447–453.
- 8. Campoy JA, Ruiz D, Egea J (2011) Dormancy in temperate fruit trees in a global warming context: A review. Sci Hortic-Amsterdam 130: 357–372.
- 9. Anderson JV, Horvath DP, Chao WS, Foley ME (2010) Bud Dormancy in Perennial Plants: A Mechanism for Survival. Top Curr Genet 21: 69–90.
- 10. Sandve SR, Kosmala A, Rudi H, Fjellheim S, Rapacz M, et al. (2011) Molecular mechanisms underlying frost tolerance in perennial grasses adapted to cold climates. Plant Sci 180: 69–77.
- 11.
Alderson J, Sharp WC (1995) Grass Varieties in the United States. USDA-SCS, vol 170.
- 12.
Vogel KP (2004) Switchgrass. In: Moser LE, Sollenberger L, Burson B, editors. Warm-season (C4) grasses. Madison, WI: ASA-CSSA-SSSA. pp. 561–588.
- 13.
Casler MD, Mitchell RB, Vogel KP (2012) Switchgrass. In: Kole CPJC, Shonnard DR, editors. Handbook of Bioenergy Crop Plants. Boca Raton, FL: CRC Press. pp. 563–592.
- 14.
Vogel K, Mitchell R, Sarath G, Casler MD (2014) Registration of 'Liberty' switchgrass. J Plant Regist.accepted
- 15.
Childs KL, Nandety A, Hirsch CN, Góngora-Castillo E, Schmutz J, et al. Generation of Transcript Assemblies and Identification of Single Nucleotide Polymorphisms from Seven Lowland and Upland Cultivars of Switchgrass. Plant Genome. accepted
- 16. Palmer NA, Saathoff AJ, Kim J, Benson A, Tobias CM, et al. (2012) Next-Generation Sequencing of Crown and Rhizome Transcriptome from an Upland, Tetraploid Switchgrass. Bioener Res 5: 649–661.
- 17. Vogel KP, Mitchell KB (2008) Heterosis in Switchgrass: Biomass Yield in Swards. Crop Sci 48: 2159–2164.
- 18. Goodstein DM, Shu SQ, Howson R, Neupane R, Hayes RD, et al. (2012) Phytozome: a comparative platform for green plant genomics. Nucleic Acids Res 40: D1178–D1186.
- 19. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, et al. (2013) TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol 14: R36.
- 20. Langmead B, Salzberg SL (2012) Fast gapped-read alignment with Bowtie 2. Nature methods 9: 357–359.
- 21. Liao Y, Smyth GK, Shi W (2014) featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30: 923–930.
- 22. Liao Y, Smyth GK, Shi W (2013) The Subread aligner: fast, accurate and scalable read mapping by seed-and-vote. Nucleic Acids Res 41: e108.
- 23. Robinson MD, McCarthy DJ, Smyth GK (2010) edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26: 139–140.
- 24.
Team RC (2011) R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2012. Available: http://cran r-project org.
- 25.
Anders S, Huber W (2010) Differential expression analysis for sequence count data. Genome Biol 11..
- 26. Conesa A, Gotz S (2008) Blast2GO: A comprehensive suite for functional analysis in plant genomics. Int J Plant Genomics 2008: 619832.
- 27. Kanehisa M, Goto S, Kawashima S, Nakaya A (2002) The KEGG databases at GenomeNet. Nucleic Acids Res 30: 42–46.
- 28. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M (2012) KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res 40: D109–114.
- 29. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, et al. (2005) Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA 102: 15545–15550.
- 30. Roessner U, Patterson JH, Forbes MG, Fincher GB, Langridge P, et al. (2006) An investigation of boron toxicity in barley using metabolomics. Plant Physiol 142: 1087–1101.
- 31. Palmer NA, Sattler SE, Saathoff AJ, Funnell D, Pedersen JF, et al. (2008) Genetic background impacts soluble and cell wall-bound aromatics in brown midrib mutants of sorghum. Planta 229: 115–127.
- 32. Saathoff AJ, Tobias CM, Sattler SE, Haas EJ, Twigg P, et al. (2010) Switchgrass Contains Two Cinnamyl Alcohol Dehydrogenases Involved in Lignin Formation. Bioener Res 4: 1–14.
- 33. Palmer NA, Sattler SE, Saathoff AJ, Sarath G (2010) A continuous, quantitative fluorescent assay for plant caffeic acid o-methyltransferases. J Agric Food Chem 58: 5220–5226.
- 34. Gomez-Vasquez R, Day R, Buschmann H, Randles S, Beeching JR, et al. (2004) Phenylpropanoids, phenylalanine ammonia lyase and peroxidases in elicitor-challenged cassava (Manihot esculenta) suspension cells and leaves. Ann Bot (Lond) 94: 87–97.
- 35. Benjamini Y, Hochberg Y (1995) Controlling the False Discovery Rate - a Practical and Powerful Approach to Multiple Testing. J Roy Stat Soc B Met 57: 289–300.
- 36.
Chi Y, Yang Y, Zhou Y, Zhou J, Fan B, et al.. (2013) Protein-Protein Interactions in the Regulation of WRKY Transcription Factors. Mol Plant.
- 37. Agarwal P, Reddy MP, Chikara J (2011) WRKY: its structure, evolutionary relationship, DNA-binding selectivity, role in stress tolerance and development of plants. Mol Biol Rep 38: 3883–3896.
- 38. Puranik S, Sahu PP, Srivastava PS, Prasad M (2012) NAC proteins: regulation and role in stress tolerance. Trends Plant Sci 17: 369–381.
- 39. Prouse MB, Campbell MM (2012) The interaction between MYB proteins and their target DNA binding sites. Biochim Biophys Acta 1819: 67–77.
- 40. Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, et al. (2010) MYB transcription factors in Arabidopsis. Trends Plant Sci 15: 573–581.
- 41. Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, et al. (2012) The Pfam protein families database. Nucleic Acids Res 40: D290–301.
- 42. Endres S, Tenhaken R (2011) Down-regulation of the myo-inositol oxygenase gene family has no effect on cell wall composition in Arabidopsis. Planta 234: 157–169.
- 43. Zhao Q, Dixon RA (2011) Transcriptional networks for lignin biosynthesis: more complex than we thought? Trends Plant Sci 16: 227–233.
- 44. Zhao C, Avci U, Grant EH, Haigler CH, Beers EP (2008) XND1, a member of the NAC domain family in Arabidopsis thaliana, negatively regulates lignocellulose synthesis and programmed cell death in xylem. Plant J 53: 425–436.
- 45. Ersoz ES, Wright MH, Pangilinan JL, Sheehan MJ, Tobias C, et al. (2012) SNP discovery with EST and NextGen sequencing in switchgrass (Panicum virgatum L.). PLoS One 7: e44112.
- 46. Palmer NA, Saathoff AJ, Waters BM, Donze T, Heng-Moss TM, et al. (2014) Global changes in mineral transporters in tetraploid switchgrasses (Panicum virgatum L.). Front Plant Sci 4: 549.
- 47. Saathoff AJ, Donze T, Palmer NA, Bradshaw J, Heng-Moss T, et al. (2013) Towards uncovering the roles of switchgrass peroxidases in plant processes. Front Plant Sci 4: 202.
- 48. Lu F, Lipka AE, Glaubitz J, Elshire R, Cherney JH, et al. (2013) Switchgrass genomic diversity, ploidy, and evolution: novel insights from a network-based SNP discovery protocol. PLoS Genet 9: e1003215.
- 49. Zhang JY, Lee YC, Torres-Jerez I, Wang M, Yin Y, et al. (2013) Development of an integrated transcript sequence database and a gene expression atlas for gene discovery and analysis in switchgrass (Panicum virgatum L.). Plant J 74: 160–173.
- 50. Zhao K, Bartley LE (2014) Comparative genomic analysis of the R2R3 MYB secondary cell wall regulators of Arabidopsis, poplar, rice, maize, and switchgrass. Bmc Plant Biol 14: 135.
- 51. Rushton PJ, Somssich IE, Ringler P, Shen QJ (2010) WRKY transcription factors. Trends Plant Sci 15: 247–258.
- 52. Chen L, Song Y, Li S, Zhang L, Zou C, et al. (2012) The role of WRKY transcription factors in plant abiotic stresses. Biochim Biophys Acta 1819: 120–128.
- 53. Tripathi P, Rabara RC, Langum TJ, Boken AK, Rushton DL, et al. (2012) The WRKY transcription factor family in Brachypodium distachyon. Bmc Genomics 13: 270.
- 54. Ulker B, Shahid Mukhtar M, Somssich IE (2007) The WRKY70 transcription factor of Arabidopsis influences both the plant senescence and defense signaling pathways. Planta 226: 125–137.
- 55. Besseau S, Li J, Palva ET (2012) WRKY54 and WRKY70 co-operate as negative regulators of leaf senescence in Arabidopsis thaliana. J Exp Bot 63: 2667–2679.
- 56. Lai Z, Vinod K, Zheng Z, Fan B, Chen Z (2008) Roles of Arabidopsis WRKY3 and WRKY4 transcription factors in plant responses to pathogens. Bmc Plant Biol 8: 68.
- 57. Sattler SE, Gilliland LU, Magallanes-Lundback M, Pollard M, DellaPenna D (2004) Vitamin E is essential for seed longevity and for preventing lipid peroxidation during germination. Plant Cell 16: 1419–1432.
- 58. Gill SS, Tuteja N (2010) Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol Biochem 48: 909–930.
- 59. Sakihama Y, Cohen MF, Grace SC, Yamasaki H (2002) Plant phenolic antioxidant and prooxidant activities: phenolics-induced oxidative damage mediated by metals in plants. Toxicology 177: 67–80.
- 60. Boerjan W, Ralph J, Baucher M (2003) Lignin biosynthesis. Annu Rev Plant Biol 54: 519–546.
- 61. Dogramaci M, Horvath DP, Christoffers MJ, Anderson JV (2011) Dehydration and vernalization treatments identify overlapping molecular networks impacting endodormancy maintenance in leafy spurge crown buds. Funct Integr Genomic 11: 611–626.
- 62. Albaugh BN, Arnold KM, Denu JM (2011) KAT(ching) metabolism by the tail: insight into the links between lysine acetyltransferases and metabolism. Chembiochem 12: 290–298.
- 63. Berr A, Shafiq S, Shen WH (2011) Histone modifications in transcriptional activation during plant development. Biochim Biophys Acta 1809: 567–576.
- 64. Xing S, Poirier Y (2012) The protein acetylome and the regulation of metabolism. Trends Plant Sci 17: 423–430.
- 65. Wu X, Oh MH, Schwarz EM, Larue CT, Sivaguru M, et al. (2011) Lysine acetylation is a widespread protein modification for diverse proteins in Arabidopsis. Plant Physiol 155: 1769–1778.