Proteomic and metabolic disturbances in lignin-modified Brachypodium distachyon.

Lignin biosynthesis begins with the deamination of phenylalanine and tyrosine (Tyr) as a key branch point between primary and secondary metabolism in land plants. Here, we used a systems biology approach to investigate the global metabolic responses to lignin pathway perturbations in the model grass Brachypodium distachyon. We identified the lignin biosynthetic protein families and found that ammonia-lyases (ALs) are among the most abundant proteins in lignifying tissues in grasses. Integrated metabolomic and proteomic data support a link between lignin biosynthesis and primary metabolism mediated by the ammonia released from ALs that is recycled for the synthesis of amino acids via glutamine. RNA interference knockdown of lignin genes confirmed that the route of the canonical pathway using shikimate ester intermediates is not essential for lignin formation in Brachypodium, and there is an alternative pathway from Tyr via sinapic acid for the synthesis of syringyl lignin involving yet uncharacterized enzymatic steps. Our findings support a model in which plant ALs play a central role in coordinating the allocation of carbon for lignin synthesis and the nitrogen available for plant growth. Collectively, these data also emphasize the value of integrative multiomic analyses to advance our understanding of plant metabolism.

Microevolution in the pansecondary metabolome of Aspergillus flavus and its potential macroevolutionary implications for filamentous fungi.

Fungi produce a wealth of pharmacologically bioactive secondary metabolites (SMs) from biosynthetic gene clusters (BGCs). It is common practice for drug discovery efforts to treat species' secondary metabolomes as being well represented by a single or a small number of representative genomes. However, this approach misses the possibility that intraspecific population dynamics, such as adaptation to environmental conditions or local microbiomes, may harbor novel BGCs that contribute to the overall niche breadth of species. Using 94 isolates of Aspergillus flavus, a cosmopolitan model fungus, sampled from seven states in the United States, we dereplicate 7,821 BGCs into 92 unique BGCs. We find that more than 25% of pangenomic BGCs show population-specific patterns of presence/absence or protein divergence. Population-specific BGCs make up most of the accessory-genome BGCs, suggesting that different ecological forces that maintain accessory genomes may be partially mediated by population-specific differences in secondary metabolism. We use ultra-high-performance high-resolution mass spectrometry to confirm that these genetic differences in BGCs also result in chemotypic differences in SM production in different populations, which could mediate ecological interactions and be acted on by selection. Thus, our results suggest a paradigm shift that previously unrealized population-level reservoirs of SM diversity may be of significant evolutionary, ecological, and pharmacological importance. Last, we find that several population-specific BGCs from A. flavus are present in Aspergillus parasiticus and Aspergillus minisclerotigenes and discuss how the microevolutionary patterns we uncover inform macroevolutionary inferences and help to align fungal secondary metabolism with existing evolutionary theory.

Expression quantitative trait loci mapping identified PtrXB38 as a key hub gene in adventitious root development in Populus.

Plant establishment requires the formation and development of an extensive root system with architecture modulated by complex genetic networks. Here, we report the identification of the PtrXB38 gene as an expression quantitative trait loci (eQTL) hotspot, mapped using 390 leaf and 444 xylem Populus trichocarpa transcriptomes. Among predicted targets of this trans-eQTL were genes involved in plant hormone responses and root development. Overexpression of PtrXB38 in Populus led to significant increases in callusing and formation of both stem-born roots and base-born adventitious roots. Omics studies revealed that genes and proteins controlling auxin transport and signaling were involved in PtrXB38-mediated adventitious root formation. Protein-protein interaction assays indicated that PtrXB38 interacts with components of endosomal sorting complexes required for transport machinery, implying that PtrXB38-regulated root development may be mediated by regulating endocytosis pathway. Taken together, this work identified a crucial root development regulator and sheds light on the discovery of other plant developmental regulators through combining eQTL mapping and omics approaches.

Outer membrane vesicles catabolize lignin-derived aromatic compounds in Pseudomonas putida KT2440.

Lignin is an abundant and recalcitrant component of plant cell walls. While lignin degradation in nature is typically attributed to fungi, growing evidence suggests that bacteria also catabolize this complex biopolymer. However, the spatiotemporal mechanisms for lignin catabolism remain unclear. Improved understanding of this biological process would aid in our collective knowledge of both carbon cycling and microbial strategies to valorize lignin to value-added compounds. Here, we examine lignin modifications and the exoproteome of three aromatic-catabolic bacteria: Pseudomonas putida KT2440, Rhodoccocus jostii RHA1, and Amycolatopsis sp. ATCC 39116. P. putida cultivation in lignin-rich media is characterized by an abundant exoproteome that is dynamically and selectively packaged into outer membrane vesicles (OMVs). Interestingly, many enzymes known to exhibit activity toward lignin-derived aromatic compounds are enriched in OMVs from early to late stationary phase, corresponding to the shift from bioavailable carbon to oligomeric lignin as a carbon source. In vivo and in vitro experiments demonstrate that enzymes contained in the OMVs are active and catabolize aromatic compounds. Taken together, this work supports OMV-mediated catabolism of lignin-derived aromatic compounds as an extracellular strategy for nutrient acquisition by soil bacteria and suggests that OMVs could potentially be useful tools for synthetic biology and biotechnological applications.

Bacterial Homologs of Progestin and AdipoQ Receptors (PAQRs) Affect Membrane Energetics Homeostasis but Not Fluidity.

Membrane potential homeostasis is essential for cell survival. Defects in membrane potential lead to pleiotropic phenotypes, consistent with the central role of membrane energetics in cell physiology. Homologs of the progestin and AdipoQ receptors (PAQRs) are conserved in multiple phyla of Bacteria and Eukarya. In eukaryotes, PAQRs are proposed to modulate membrane fluidity and fatty acid (FA) metabolism. The role of bacterial homologs has not been elucidated. Here, we use Escherichia coli and Bacillus subtilis to show that bacterial PAQR homologs, which we name "TrhA," have a role in membrane energetics homeostasis. Using transcriptional fusions, we show that E. coli TrhA (encoded by yqfA) is part of the unsaturated fatty acid biosynthesis regulon. Fatty acid analyses and physiological assays show that a lack of TrhA in both E. coli and B. subtilis (encoded by yplQ) provokes subtle but consistent changes in membrane fatty acid profiles that do not translate to control of membrane fluidity. Instead, membrane proteomics in E. coli suggested a disrupted energy metabolism and dysregulated membrane energetics in the mutant, though it grew similarly to its parent. These changes translated into a disturbed membrane potential in the mutant relative to its parent under various growth conditions. Similar dysregulation of membrane energetics was observed in a different E. coli strain and in the distantly related B. subtilis. Together, our findings are consistent with a role for TrhA in membrane energetics homeostasis, through a mechanism that remains to be elucidated. IMPORTANCE Eukaryotic homologs of the progestin and AdipoQ receptor family (PAQR) have been shown to regulate membrane fluidity by affecting, through unknown mechanisms, unsaturated fatty acid (FA) metabolism. The bacterial homologs studied here mediate small and consistent changes in unsaturated FA metabolism that do not seem to impact membrane fluidity but, rather, alter membrane energetics homeostasis. Together, the findings here suggest that bacterial and eukaryotic PAQRs share functions in maintaining membrane homeostasis (fluidity in eukaryotes and energetics for bacteria with TrhA homologs).

Development of an Experimental Approach to Achieve Spatially Resolved Plant Root-Associated Metaproteomics Using an Agar-Plate System.

Plant-microbe interactions in the rhizosphere play a vital role in plant health and productivity. The composition and function of root-associated microbiomes is strongly influenced by their surrounding environment, which is often customized by their host. How microbiomes change with respect to space and time across plant roots remains poorly understood, and methodologies that facilitate spatiotemporal metaproteomic studies of root-associated microbiomes are yet to be realized. Here, we developed a method that provides spatially resolved metaproteome measurements along plant roots embedded in agar-plate culture systems, which have long been used to study plants. Spatially defined agar "plugs" of interest were excised and subsequently processed using a novel peptide extraction method prior to metaproteomics, which was used to infer both microbial community composition and function. As a proof-of-principle, a previously studied 10-member community constructed from a Populus root system was grown in an agar plate with a 3-week-old Populus trichocarpa plant. Metaproteomics was performed across two time points (24 and 48 h) for three distinct locations (root base, root tip, and a region distant from the root). The spatial resolution of these measurements provides evidence that microbiome composition and expression changes across the plant root interface. Interrogation of the individual microbial proteomes revealed functional profiles related to their behavioral associations with the plant root, in which chemotaxis and augmented metabolism likely supported predominance of the most abundant member. This study demonstrated a novel peptide extraction method for studying plant agar-plate culture systems, which was previously unsuitable for (meta)proteomic measurements.

Ecosystem consequences of introducing plant growth promoting rhizobacteria to managed systems and potential legacy effects.

The rapidly growing industry of crop biostimulants leverages the application of plant growth promoting rhizobacteria (PGPR) to promote plant growth and health. However, introducing nonnative rhizobacteria may impact other aspects of ecosystem functioning and have legacy effects; these potential consequences are largely unexplored. Nontarget consequences of PGPR may include changes in resident microbiomes, nutrient cycling, pollinator services, functioning of other herbivores, disease suppression, and organic matter persistence. Importantly, we lack knowledge of whether these ecosystem effects may manifest in adjacent ecosystems. The introduced PGPR can leave a functional legacy whether they persist in the community or not. Legacy effects include shifts in resident microbiomes and their temporal dynamics, horizontal transfer of genes from the PGPR to resident taxa, and changes in resident functional groups and interaction networks. Ecosystem functions may be affected by legacies PGPR leave following niche construction, such as when PGPR alter soil pH that in turn alters biogeochemical cycling rates. Here, we highlight new research directions to elucidate how introduced PGPR impact resident microbiomes and ecosystem functions and their capacity for legacy effects.

Temporal dynamics of protein and post-translational modification abundances in Populus leaf across a diurnal period.

Populus spp. are dedicated woody biomass feedstocks for advanced biofuels and bioproducts. Proper growth and fitness of poplar as a sustainable feedstock depends on timely perception and response to environmental signals (e.g., light, temperature, water). Poplar leaves, like other C3 photosynthesis plants, have evolved oscillating or circadian rhythms that play important roles in synchronizing biological processes with external cues. To characterize this phenomenon at a molecular level, we employed bottom-up proteomics using high-resolution mass spectrometry and de novo-assisted database searching to identify abundance changes in proteins and post-translational modifications in poplar leaf tissue sampled across a 12/12-hour light/dark diurnal period.

The Moderately (D)efficient Enzyme: Catalysis-Related Damage In Vivo and Its Repair.

Enzymes have in vivo life spans. Analysis of life spans, i.e., lifetime totals of catalytic turnovers, suggests that nonsurvivable collateral chemical damage from the very reactions that enzymes catalyze is a common but underdiagnosed cause of enzyme death. Analysis also implies that many enzymes are moderately deficient in that their active-site regions are not naturally as hardened against such collateral damage as they could be, leaving room for improvement by rational design or directed evolution. Enzyme life span might also be improved by engineering systems that repair otherwise fatal active-site damage, of which a handful are known and more are inferred to exist. Unfortunately, the data needed to design and execute such improvements are lacking: there are too few measurements of in vivo life span, and existing information about the extent, nature, and mechanisms of active-site damage and repair during normal enzyme operation is too scarce, anecdotal, and speculative to act on. Fortunately, advances in proteomics, metabolomics, cheminformatics, comparative genomics, and structural biochemistry now empower a systematic, data-driven approach for identifying, predicting, and validating instances of active-site damage and its repair. These capabilities would be practically useful in enzyme redesign and improvement of in-use stability and could change our thinking about which enzymes die young in vivo, and why.

Peeling back the layers of crassulacean acid metabolism: functional differentiation between Kalanchoë fedtschenkoi epidermis and mesophyll proteomes.

Crassulacean acid metabolism (CAM) is a specialized mode of photosynthesis that offers the potential to engineer improved water-use efficiency (WUE) and drought resilience in C3 plants while sustaining productivity in the hotter and drier climates that are predicted for much of the world. CAM species show an inverted pattern of stomatal opening and closing across the diel cycle, which conserves water and provides a means of maintaining growth in hot, water-limited environments. Recent genome sequencing of the constitutive model CAM species Kalanchoë fedtschenkoi provides a platform for elucidating the ensemble of proteins that link photosynthetic metabolism with stomatal movement, and that protect CAM plants from harsh environmental conditions. We describe a large-scale proteomics analysis to characterize and compare proteins, as well as diel changes in their abundance in guard cell-enriched epidermis and mesophyll cells from leaves of K. fedtschenkoi. Proteins implicated in processes that encompass respiration, the transport of water and CO2 , stomatal regulation, and CAM biochemistry are highlighted and discussed. Diel rescheduling of guard cell starch turnover in K. fedtschenkoi compared with that observed in Arabidopsis is reported and tissue-specific localization in the epidermis and mesophyll of isozymes implicated in starch and malate turnover are discussed in line with the contrasting roles for these metabolites within the CAM mesophyll and stomatal complex. These data reveal the proteins and the biological processes enriched in each layer and provide key information for studies aiming to adapt plants to hot and dry environments by modifying leaf physiology for improved plant sustainability.

Towards engineering ectomycorrhization into switchgrass bioenergy crops via a lectin receptor-like kinase.

Soil-borne microbes can establish compatible relationships with host plants, providing a large variety of nutritive and protective compounds in exchange for photosynthesized sugars. However, the molecular mechanisms mediating the establishment of these beneficial relationships remain unclear. Our previous genetic mapping and whole-genome resequencing studies identified a gene deletion event of a Populus trichocarpa lectin receptor-like kinase gene PtLecRLK1 in Populus deltoides that was associated with poor-root colonization by the ectomycorrhizal fungus Laccaria bicolor. By introducing PtLecRLK1 into a perennial grass known to be a non-host of L. bicolor, switchgrass (Panicum virgatum L.), we found that L. bicolor colonizes ZmUbipro-PtLecRLK1 transgenic switchgrass roots, which illustrates that the introduction of PtLecRLK1 has the potential to convert a non-host to a host of L. bicolor. Furthermore, transcriptomic and proteomic analyses on inoculated-transgenic switchgrass roots revealed genes/proteins overrepresented in the compatible interaction and underrepresented in the pathogenic defence pathway, consistent with the view that pathogenic defence response is down-regulated during compatible interaction. Metabolomic profiling revealed that root colonization in the transgenic switchgrass was associated with an increase in N-containing metabolites and a decrease in organic acids, sugars, and aromatic hydroxycinnamate conjugates, which are often seen in the early steps of establishing compatible interactions. These studies illustrate that PtLecRLK1 is able to render a plant susceptible to colonization by the ectomycorrhizal fungus L. bicolor and shed light on engineering mycorrhizal symbiosis into a non-host to enhance plant productivity and fitness on marginal lands.

Metabolism of syringyl lignin-derived compounds in Pseudomonas putida enables convergent production of 2-pyrone-4,6-dicarboxylic acid.

Valorization of lignin, an abundant component of plant cell walls, is critical to enabling the lignocellulosic bioeconomy. Biological funneling using microbial biocatalysts has emerged as an attractive approach to convert complex mixtures of lignin depolymerization products to value-added compounds. Ideally, biocatalysts would convert aromatic compounds derived from the three canonical types of lignin: syringyl (S), guaiacyl (G), and p-hydroxyphenyl (H). Pseudomonas putida KT2440 (hereafter KT2440) has been developed as a biocatalyst owing in part to its native catabolic capabilities but is not known to catabolize S-type lignin-derived compounds. Here, we demonstrate that syringate, a common S-type lignin-derived compound, is utilized by KT2440 only in the presence of another energy source or when vanAB was overexpressed, as syringate was found to be O-demethylated to gallate by VanAB, a two-component monooxygenase, and further catabolized via extradiol cleavage. Unexpectedly, the specificity (kcat/KM) of VanAB for syringate was within 25% that for vanillate and O-demethylation of both substrates was well-coupled to O2 consumption. However, the native KT2440 gallate-cleaving dioxygenase, GalA, was potently inactivated by 3-O-methylgallate. To engineer a biocatalyst to simultaneously convert S-, G-, and H-type monomers, we therefore employed VanAB from Pseudomonas sp. HR199, which has lower activity for 3MGA, and LigAB, an extradiol dioxygenase able to cleave protocatechuate and 3-O-methylgallate. This strain converted 93% of a mixture of lignin monomers to 2-pyrone-4,6-dicarboxylate, a promising bio-based chemical. Overall, this study elucidates a native pathway in KT2440 for catabolizing S-type lignin-derived compounds and demonstrates the potential of this robust chassis for lignin valorization.

Metaproteomics reveals insights into microbial structure, interactions, and dynamic regulation in defined communities as they respond to environmental disturbance.

BACKGROUND: Microbe-microbe interactions between members of the plant rhizosphere are important but remain poorly understood. A more comprehensive understanding of the molecular mechanisms used by microbes to cooperate, compete, and persist has been challenging because of the complexity of natural ecosystems and the limited control over environmental factors. One strategy to address this challenge relies on studying complexity in a progressive manner, by first building a detailed understanding of relatively simple subsets of the community and then achieving high predictive power through combining different building blocks (e.g., hosts, community members) for different environments. Herein, we coupled this reductionist approach with high-resolution mass spectrometry-based metaproteomics to study molecular mechanisms driving community assembly, adaptation, and functionality for a defined community of ten taxonomically diverse bacterial members of Populus deltoides rhizosphere co-cultured either in a complex or defined medium. RESULTS: Metaproteomics showed this defined community assembled into distinct microbiomes based on growth media that eventually exhibit composition and functional stability over time. The community grown in two different media showed variation in composition, yet both were dominated by only a few microbial strains. Proteome-wide interrogation provided detailed insights into the functional behavior of each dominant member as they adjust to changing community compositions and environments. The emergence and persistence of select microbes in these communities were driven by specialization in strategies including motility, antibiotic production, altered metabolism, and dormancy. Protein-level interrogation identified post-translational modifications that provided additional insights into regulatory mechanisms influencing microbial adaptation in the changing environments. CONCLUSIONS: This study provides high-resolution proteome-level insights into our understanding of microbe-microbe interactions and highlights specialized biological processes carried out by specific members of assembled microbiomes to compete and persist in changing environmental conditions. Emergent properties observed in these lower complexity communities can then be re-evaluated as more complex systems are studied and, when a particular property becomes less relevant, higher-order interactions can be identified.

Molecular Remodeling in Populus PdKOR RNAi Roots Profiled Using LC-MS/MS Proteomics.

Plant endo-ß-1,4-glucanases belonging to the Glycoside Hydrolase Family 9 have functional roles in cell wall biosynthesis and remodeling via endohydrolysis of (1→4)-ß-d-glucosidic linkages. Modification of cell wall chemistry via RNA interference (RNAi)-mediated downregulation of Populus deltoides KORRIGAN (PdKOR), an endo-ß-1,4-glucanase familygene was shown to have functional consequences on the composition of secondary metabolome and the ability of modified roots to interact with beneficial microbes. The molecular remodeling that underlies the observed differences at metabolic, physiological, and morphological levels in roots is not well understood. Here a liquid chromatography (LC)-tandem mass spectrometry (MS/MS)-based proteome profiling approach is used to survey the molecular remodeling in root tissues of PdKOR and control plants. A total of 14316 peptides are identified and these mapped to 7139 P. deltoides proteins. Based on 90% sequence identity, the measured protein accessions represent 1187 functional protein groups. Analysis of Gene Ontology (GO) categories and specific individual proteins show differential expression of proteins relevant to plant-microbe interactions, cell wall chemistry, and metabolism. The new proteome dataset serves as a useful resource for deriving new hypotheses and empirical testing pertaining to functional roles of proteins and pathways in differential priming of plant roots to interactions with microbes.

A Viable New Strategy for the Discovery of Peptide Proteolytic Cleavage Products in Plant-Microbe Interactions.

Small peptides that are proteolytic cleavage products (PCPs) of less than 100 amino acids are emerging as key signaling molecules that mediate cell-to-cell communication and biological processes that occur between and within plants, fungi, and bacteria. Yet, the discovery and characterization of these molecules is largely overlooked. Today, selective enrichment and subsequent characterization by mass spectrometry-based sequencing offers the greatest potential for their comprehensive characterization, however qualitative and quantitative performance metrics are rarely captured. Herein, we addressed this need by benchmarking the performance of an enrichment strategy, optimized specifically for small PCPs, using state-of-the-art de novo-assisted peptide sequencing. As a case study, we implemented this approach to identify PCPs from different root and foliar tissues of the hybrid poplar Populus × canescens 717-1B4 in interaction with the ectomycorrhizal basidiomycete Laccaria bicolor. In total, we identified 1,660 and 2,870 Populus and L. bicolor unique PCPs, respectively. Qualitative results supported the identification of well-known PCPs, like the mature form of the photosystem II complex 5-kDa protein (approximately 3 kDa). A total of 157 PCPs were determined to be significantly more abundant in root tips with established ectomycorrhiza when compared with root tips without established ectomycorrhiza and extramatrical mycelium of L. bicolor. These PCPs mapped to 64 Populus proteins and 69 L. bicolor proteins in our database, with several of them previously implicated in biologically relevant associations between plant and fungus.

Correction to: Diel rewiring and positive selection of ancient plant proteins enabled evolution of CAM photosynthesis in Agave.

Following publication of the original article.

Exploiting the Dynamic Relationship between Peptide Separation Quality and Peptide Coisolation in a Multiple-Peptide Matches-per-Spectrum Approach Offers a Strategy To Optimize Bottom-Up Proteomics Throughput and Depth.

Peptide cofragmentation leads to chimeric MS/MS spectra that negatively impact traditional single-peptide match-per-spectrum (sPSM) search strategies in proteomics. The collection of chimeric spectra is influenced by peptide coelution and the width of precursor isolation windows. Although peptide cofragmentation can be reduced by advanced chromatography, such as UHPLC and 2D-HPLC separation schemes, and narrower isolation windows, chimeric spectra can still be as high as 30-50% of the total MS/MS spectra collected. Alternatively, cofragmented peptides in chimeric spectra and the use of wider isolation windows benefit multiple-peptide matches-per-spectrum (mPSM) algorithms, such as CharmeRT, which facilitate the identification of several cofragmented peptides. Considering recent advancements in LC and mPSM methodologies, we present a comprehensive examination of the levels of chimeric spectra collected in the analysis of a HeLa digest measured using different LC modes of separation and isolation windows and compare the depth of identifications obtained when these data are annotated using a sPSM or a mPSM approach. Our results demonstrate that MS/MS data derived from 1D-HPLC strategies under different gradient schemes and searched with CharmeRT yielded higher average numbers of PSMs (11%-49%), peptide analytes (10%-16%), and peptide sequences (3%-10%) compared to data derived from 1D-UHPLC runs but searched with a sPSM strategy. Interestingly, data from a 2D-HPLC separation strategy benefits more from the application of CharmeRT results when compared to a 50 cm 1D-UHPLC column employing a 500 min gradient. Overall, these results provide new insights into how to better configure LC-MS/MS measurements for improved throughput and peptide identification in complex proteomes.

The nature of the progression of drought stress drives differential metabolomic responses in Populus deltoides.

BACKGROUND AND AIMS: The use of woody crops for Quad-level (approx. 1 × 1018 J) energy production will require marginal agricultural lands that experience recurrent periods of water stress. Populus species have the capacity to increase dehydration tolerance by lowering osmotic potential via osmotic adjustment. The aim of this study was to investigate how the inherent genetic potential of a Populus clone to respond to drought interacts with the nature of the drought to determine the degree of biochemical response. METHODS: A greenhouse drought stress study was conducted on Populus deltoides 'WV94' and the resulting metabolite profiles of leaves were determined by gas chromatography-mass spectrometry following trimethylsilylation for plants subjected to cyclic mild (-0.5 MPa pre-dawn leaf water potential) drought vs. cyclic severe (-1.26 MPa) drought in contrast to well-watered controls (-0.1 MPa) after two or four drought cycles, and in contrast to plants subjected to acute drought, where plants were desiccated for up to 8 d. KEY RESULTS: The nature of drought (cyclic vs. acute), frequency of drought (number of cycles) and the severity of drought (mild vs. severe) all dictated the degree of osmotic adjustment and the nature of the organic solutes that accumulated. Whereas cyclic drought induced the largest responses in primary metabolism (soluble sugars, organic acids and amino acids), acute onset of prolonged drought induced the greatest osmotic adjustment and largest responses in secondary metabolism, especially populosides (hydroxycinnamic acid conjugates of salicin). CONCLUSIONS: The differential adaptive metabolite responses in cyclic vs. acute drought suggest that stress acclimation occurs via primary metabolism in response to cyclic drought, whereas expanded metabolic plasticity occurs via secondary metabolism following severe, acute drought. The shift in carbon partitioning to aromatic metabolism with the production of a diverse suite of higher order salicylates lowers osmotic potential and increases the probability of post-stress recovery.

Diel rewiring and positive selection of ancient plant proteins enabled evolution of CAM photosynthesis in Agave.

BACKGROUND: Crassulacean acid metabolism (CAM) enhances plant water-use efficiency through an inverse day/night pattern of stomatal closure/opening that facilitates nocturnal CO2 uptake. CAM has evolved independently in over 35 plant lineages, accounting for ~ 6% of all higher plants. Agave species are highly heat- and drought-tolerant, and have been domesticated as model CAM crops for beverage, fiber, and biofuel production in semi-arid and arid regions. However, the genomic basis of evolutionary innovation of CAM in genus Agave is largely unknown. RESULTS: Using an approach that integrated genomics, gene co-expression networks, comparative genomics and protein structure analyses, we investigated the molecular evolution of CAM as exemplified in Agave. Comparative genomics analyses among C3, C4 and CAM species revealed that core metabolic components required for CAM have ancient genomic origins traceable to non-vascular plants while regulatory proteins required for diel re-programming of metabolism have a more recent origin shared among C3, C4 and CAM species. We showed that accelerated evolution of key functional domains in proteins responsible for primary metabolism and signaling, together with a diel re-programming of the transcription of genes involved in carbon fixation, carbohydrate processing, redox homeostasis, and circadian control is required for the evolution of CAM in Agave. Furthermore, we highlighted the potential candidates contributing to the adaptation of CAM functional modules. CONCLUSIONS: This work provides evidence of adaptive evolution of CAM related pathways. We showed that the core metabolic components required for CAM are shared by non-vascular plants, but regulatory proteins involved in re-reprogramming of carbon fixation and metabolite transportation appeared more recently. We propose that the accelerated evolution of key proteins together with a diel re-programming of gene expression were required for CAM evolution from C3 ancestors in Agave.