Information theory and machine learning illuminate large-scale metabolomic responses of Brachypodium distachyon to environmental change.

Plant responses to environmental change are mediated via changes in cellular metabolomes. However, <5% of signals obtained from liquid chromatography tandem mass spectrometry (LC-MS/MS) can be identified, limiting our understanding of how metabolomes change under biotic/abiotic stress. To address this challenge, we performed untargeted LC-MS/MS of leaves, roots, and other organs of Brachypodium distachyon (Poaceae) under 17 organ-condition combinations, including copper deficiency, heat stress, low phosphate, and arbuscular mycorrhizal symbiosis. We found that both leaf and root metabolomes were significantly affected by the growth medium. Leaf metabolomes were more diverse than root metabolomes, but the latter were more specialized and more responsive to environmental change. We found that 1 week of copper deficiency shielded the root, but not the leaf metabolome, from perturbation due to heat stress. Machine learning (ML)-based analysis annotated approximately 81% of the fragmented peaks versus approximately 6% using spectral matches alone. We performed one of the most extensive validations of ML-based peak annotations in plants using thousands of authentic standards, and analyzed approximately 37% of the annotated peaks based on these assessments. Analyzing responsiveness of each predicted metabolite class to environmental change revealed significant perturbations of glycerophospholipids, sphingolipids, and flavonoids. Co-accumulation analysis further identified condition-specific biomarkers. To make these results accessible, we developed a visualization platform on the Bio-Analytic Resource for Plant Biology website (https://bar.utoronto.ca/efp_brachypodium_metabolites/cgi-bin/efpWeb.cgi), where perturbed metabolite classes can be readily visualized. Overall, our study illustrates how emerging chemoinformatic methods can be applied to reveal novel insights into the dynamic plant metabolome and stress adaptation.

Orthology-based analysis helps map evolutionary diversification and predict substrate class use of BAHD acyltransferases.

Large enzyme families catalyze metabolic diversification by virtue of their ability to use diverse chemical scaffolds. How enzyme families attain such functional diversity is not clear. Furthermore, duplication and promiscuity in such enzyme families limits their functional prediction, which has produced a burgeoning set of incompletely annotated genes in plant genomes. Here, we address these challenges using BAHD acyltransferases as a model. This fast-evolving family expanded drastically in land plants, increasing from one to five copies in algae to approximately 100 copies in diploid angiosperm genomes. Compilation of >160 published activities helped visualize the chemical space occupied by this family and define eight different classes based on structural similarities between acceptor substrates. Using orthologous groups (OGs) across 52 sequenced plant genomes, we developed a method to predict BAHD acceptor substrate class utilization as well as origins of individual BAHD OGs in plant evolution. This method was validated using six novel and 28 previously characterized enzymes and helped improve putative substrate class predictions for BAHDs in the tomato genome. Our results also revealed that while cuticular wax and lignin biosynthetic activities were more ancient, anthocyanin acylation activity was fixed in BAHDs later near the origin of angiosperms. The OG-based analysis enabled identification of signature motifs in anthocyanin-acylating BAHDs, whose importance was validated via molecular dynamic simulations, site-directed mutagenesis and kinetic assays. Our results not only describe how BAHDs contributed to evolution of multiple chemical phenotypes in the plant world but also propose a biocuration-enabled approach for improved functional annotation of plant enzyme families.

A massively parallel barcoded sequencing pipeline enables generation of the first ORFeome and interactome map for rice.

Systematic mappings of protein interactome networks have provided invaluable functional information for numerous model organisms. Here we develop PCR-mediated Linkage of barcoded Adapters To nucleic acid Elements for sequencing (PLATE-seq) that serves as a general tool to rapidly sequence thousands of DNA elements. We validate its utility by generating the ORFeome for Oryza sativa covering 2,300 genes and constructing a high-quality protein-protein interactome map consisting of 322 interactions between 289 proteins, expanding the known interactions in rice by roughly 50%. Our work paves the way for high-throughput profiling of protein-protein interactions in a wide range of organisms.

Acylsugars protect Nicotiana benthamiana against insect herbivory and desiccation.

KEY MESSAGE: Nicotiana benthamiana acylsugar acyltransferase (ASAT) is required for protection against desiccation and insect herbivory. Knockout mutations provide a new resource for investigation of plant-aphid and plant-whitefly interactions. Nicotiana benthamiana is used extensively as a transient expression platform for functional analysis of genes from other species. Acylsugars, which are produced in the trichomes, are a hypothesized cause of the relatively high insect resistance that is observed in N. benthamiana. We characterized the N. benthamiana acylsugar profile, bioinformatically identified two acylsugar acyltransferase genes, ASAT1 and ASAT2, and used CRISPR/Cas9 mutagenesis to produce acylsugar-deficient plants for investigation of insect resistance and foliar water loss. Whereas asat1 mutations reduced accumulation, asat2 mutations caused almost complete depletion of foliar acylsucroses. Three hemipteran and three lepidopteran herbivores survived, gained weight, and/or reproduced significantly better on asat2 mutants than on wildtype N. benthamiana. Both asat1 and asat2 mutations reduced the water content and increased leaf temperature. Our results demonstrate the specific function of two ASAT proteins in N. benthamiana acylsugar biosynthesis, insect resistance, and desiccation tolerance. The improved growth of aphids and whiteflies on asat2 mutants will facilitate the use of N. benthamiana as a transient expression platform for the functional analysis of insect effectors and resistance genes from other plant species. Similarly, the absence of acylsugars in asat2 mutants will enable analysis of acylsugar biosynthesis genes from other Solanaceae by transient expression.

Migration through a Major Andean Ecogeographic Disruption as a Driver of Genetic and Phenotypic Diversity in a Wild Tomato Species.

Evolutionary dynamics at the population level play a central role in creating the diversity of life on our planet. In this study, we sought to understand the origins of such population-level variation in mating systems and defensive acylsugar chemistry in Solanum habrochaites-a wild tomato species found in diverse Andean habitats in Ecuador and Peru. Using Restriction-site-Associated-DNA-Sequencing (RAD-seq) of 50 S. habrochaites accessions, we identified eight population clusters generated via isolation and hybridization dynamics of 4-6 ancestral populations. Detailed characterization of mating systems of these clusters revealed emergence of multiple self-compatible (SC) groups from progenitor self-incompatible populations in the northern part of the species range. Emergence of these SC groups was also associated with fixation of deleterious alleles inactivating acylsugar acetylation. The Amotape-Huancabamba Zone-a geographical landmark in the Andes with high endemism and isolated microhabitats-was identified as a major driver of differentiation in the northern species range, whereas large geographical distances contributed to population structure and evolution of a novel SC group in the central and southern parts of the range, where the species was also inferred to have originated. Findings presented here highlight the role of the diverse ecogeography of Peru and Ecuador in generating population differentiation, and enhance our understanding of the microevolutionary processes that create biological diversity.

Specialized Metabolism in a Nonmodel Nightshade: Trichome Acylinositol Biosynthesis.

Plants make many biologically active, specialized metabolites, which vary in structure, biosynthesis, and the processes they influence. An increasing number of these compounds are documented to protect plants from insects, pathogens, or herbivores or to mediate interactions with beneficial organisms, including pollinators and nitrogen-fixing microbes. Acylsugars, one class of protective compounds, are made in glandular trichomes of plants across the Solanaceae family. While most described acylsugars are acylsucroses, published examples also include acylsugars with hexose cores. The South American fruit crop naranjilla (lulo; Solanum quitoense) produces acylsugars containing a myoinositol core. We identified an enzyme that acetylates triacylinositols, a function homologous to the last step in the acylsucrose biosynthetic pathway of tomato (Solanum lycopersicum). Our analysis reveals parallels between S. lycopersicum acylsucrose and S. quitoense acylinositol biosynthesis, suggesting a common evolutionary origin.

Coexpression Analysis Identifies Two Oxidoreductases Involved in the Biosynthesis of the Monoterpene Acid Moiety of Natural Pyrethrin Insecticides in Tanacetum cinerariifolium.

Flowers of Tanacetum cinerariifolium produce a set of compounds known collectively as pyrethrins, which are commercially important pesticides that are strongly toxic to flying insects but not to most vertebrates. A pyrethrin molecule is an ester consisting of either trans-chrysanthemic acid or its modified form, pyrethric acid, and one of three alcohols, jasmolone, pyrethrolone, and cinerolone, that appear to be derived from jasmonic acid. Chrysanthemyl diphosphate synthase (CDS), the first enzyme involved in the synthesis of trans-chrysanthemic acid, was characterized previously and its gene isolated. TcCDS produces free trans-chrysanthemol in addition to trans-chrysanthemyl diphosphate, but the enzymes responsible for the conversion of trans-chrysanthemol to the corresponding aldehyde and then to the acid have not been reported. We used an RNA sequencing-based approach and coexpression correlation analysis to identify several candidate genes encoding putative trans-chrysanthemol and trans-chrysanthemal dehydrogenases. We functionally characterized the proteins encoded by these genes using a combination of in vitro biochemical assays and heterologous expression in planta to demonstrate that TcADH2 encodes an enzyme that oxidizes trans-chrysanthemol to trans-chrysanthemal, while TcALDH1 encodes an enzyme that oxidizes trans-chrysanthemal into trans-chrysanthemic acid. Transient coexpression of TcADH2 and TcALDH1 together with TcCDS in Nicotiana benthamiana leaves results in the production of trans-chrysanthemic acid as well as several other side products. The majority (58%) of trans-chrysanthemic acid was glycosylated or otherwise modified. Overall, these data identify key steps in the biosynthesis of pyrethrins and demonstrate the feasibility of metabolic engineering to produce components of these defense compounds in a heterologous host.

Characteristics of Plant Essential Genes Allow for within- and between-Species Prediction of Lethal Mutant Phenotypes.

Essential genes represent critical cellular components whose disruption results in lethality. Characteristics shared among essential genes have been uncovered in fungal and metazoan model systems. However, features associated with plant essential genes are largely unknown and the full set of essential genes remains to be discovered in any plant species. Here, we show that essential genes in Arabidopsis thaliana have distinct features useful for constructing within- and cross-species prediction models. Essential genes in A. thaliana are often single copy or derived from older duplications, highly and broadly expressed, slow evolving, and highly connected within molecular networks compared with genes with nonlethal mutant phenotypes. These gene features allowed the application of machine learning methods that predicted known lethal genes as well as an additional 1970 likely essential genes without documented phenotypes. Prediction models from A. thaliana could also be applied to predict Oryza sativa and Saccharomyces cerevisiae essential genes. Importantly, successful predictions drew upon many features, while any single feature was not sufficient. Our findings show that essential genes can be distinguished from genes with nonlethal phenotypes using features that are similar across kingdoms and indicate the possibility for translational application of our approach to species without extensive functional genomic and phenomic resources.

Functionally divergent alleles and duplicated Loci encoding an acyltransferase contribute to acylsugar metabolite diversity in Solanum trichomes.

Glandular trichomes from tomato (Solanum lycopersicum) and other species in the Solanaceae produce and secrete a mixture of O-acylsugars (aliphatic esters of sucrose and glucose) that contribute to insect defense. Despite their phylogenetic distribution and diversity, relatively little is known about how these specialized metabolites are synthesized. Mass spectrometric profiling of acylsugars in the S. lycopersicum x Solanum pennellii introgression lines identified a chromosome 11 locus containing a cluster of BAHD acyltransferases with one gene (named Sl-ASAT3) expressed in tip cells of type I trichomes where acylsugars are made. Sl-ASAT3 was shown to encode an acyl-CoA-dependent acyltransferase that catalyzes the transfer of short (four to five carbons) branched acyl chains to the furanose ring of di-acylsucrose acceptors to produce tri-acylsucroses, which can be further acetylated by Sl-ASAT4 (previously Sl-AT2). Among the wild tomatoes, diversity in furanose ring acyl chains on acylsucroses was most striking in Solanum habrochaites. S. habrochaites accessions from Ecuador and northern Peru produced acylsucroses with short (≤C5) or no acyl chains on the furanose ring. Accessions from central and southern Peru had the ability to add short or long (up to C12) acyl chains to the furanose ring. Multiple ASAT3-like sequences were found in most accessions, and their in vitro activities correlated with observed geographical diversity in acylsugar profiles.

Consequences of Whole-Genome Triplication as Revealed by Comparative Genomic Analyses of the Wild Radish Raphanus raphanistrum and Three Other Brassicaceae Species.

Polyploidization events are frequent among flowering plants, and the duplicate genes produced via such events contribute significantly to plant evolution. We sequenced the genome of wild radish (Raphanus raphanistrum), a Brassicaceae species that experienced a whole-genome triplication event prior to diverging from Brassica rapa. Despite substantial gene gains in these two species compared with Arabidopsis thaliana and Arabidopsis lyrata, ∼70% of the orthologous groups experienced gene losses in R. raphanistrum and B. rapa, with most of the losses occurring prior to their divergence. The retained duplicates show substantial divergence in sequence and expression. Based on comparison of A. thaliana and R. raphanistrum ortholog floral expression levels, retained radish duplicates diverged primarily via maintenance of ancestral expression level in one copy and reduction of expression level in others. In addition, retained duplicates differed significantly from genes that reverted to singleton state in function, sequence composition, expression patterns, network connectivity, and rates of evolution. Using these properties, we established a statistical learning model for predicting whether a duplicate would be retained postpolyploidization. Overall, our study provides new insights into the processes of plant duplicate loss, retention, and functional divergence and highlights the need for further understanding factors controlling duplicate gene fate.

Something Old, Something New: Conserved Enzymes and the Evolution of Novelty in Plant Specialized Metabolism.

Plants produce hundreds of thousands of small molecules known as specialized metabolites, many of which are of economic and ecological importance. This remarkable variety is a consequence of the diversity and rapid evolution of specialized metabolic pathways. These novel biosynthetic pathways originate via gene duplication or by functional divergence of existing genes, and they subsequently evolve through selection and/or drift. Studies over the past two decades revealed that diverse specialized metabolic pathways have resulted from the incorporation of primary metabolic enzymes. We discuss examples of enzyme recruitment from primary metabolism and the variety of paths taken by duplicated primary metabolic enzymes toward integration into specialized metabolism. These examples provide insight into processes by which plant specialized metabolic pathways evolve and suggest approaches to discover enzymes of previously uncharacterized metabolic networks.

Molecular Evidence for Functional Divergence and Decay of a Transcription Factor Derived from Whole-Genome Duplication in Arabidopsis thaliana.

Functional divergence between duplicate transcription factors (TFs) has been linked to critical events in the evolution of land plants and can result from changes in patterns of expression, binding site divergence, and/or interactions with other proteins. Although plant TFs tend to be retained post polyploidization, many are lost within tens to hundreds of million years. Thus, it can be hypothesized that some TFs in plant genomes are in the process of becoming pseudogenes. Here, we use a pair of salt tolerance-conferring transcription factors, DWARF AND DELAYED FLOWERING1 (DDF1) and DDF2, that duplicated through paleopolyploidy 50 to 65 million years ago, as examples to illustrate potential mechanisms leading to duplicate retention and loss. We found that the expression patterns of Arabidopsis thaliana (At)DDF1 and AtDDF2 have diverged in a highly asymmetric manner, and AtDDF2 has lost most inferred ancestral stress responses. Consistent with promoter disablement, the AtDDF2 promoter has fewer predicted cis-elements and a methylated repetitive element. Through comparisons of AtDDF1, AtDDF2, and their Arabidopsis lyrata orthologs, we identified significant differences in binding affinities and binding site preference. In particular, an AtDDF2-specific substitution within the DNA-binding domain significantly reduces binding affinity. Cross-species analyses indicate that both AtDDF1 and AtDDF2 are under selective constraint, but among A. thaliana accessions, AtDDF2 has a higher level of nonsynonymous nucleotide diversity compared with AtDDF1. This may be the result of selection in different environments or may point toward the possibility of ongoing functional decay despite retention for millions of years after gene duplication.

A Feedback-Insensitive Isopropylmalate Synthase Affects Acylsugar Composition in Cultivated and Wild Tomato.

Acylsugars are insecticidal specialized metabolites produced in the glandular trichomes of plants in the Solanaceae family. In the tomato clade of the Solanum genus, acylsugars consist of aliphatic acids of different chain lengths esterified to sucrose, or less frequently to glucose. Through liquid chromatography-mass spectrometry screening of introgression lines, we previously identified a region of chromosome 8 in the Solanum pennellii LA0716 genome (IL8-1/8-1-1) that causes the cultivated tomato Solanum lycopersicum to shift from producing acylsucroses with abundant 3-methylbutanoic acid acyl chains derived from leucine metabolism to 2-methylpropanoic acid acyl chains derived from valine metabolism. We describe multiple lines of evidence implicating a trichome-expressed gene from this region as playing a role in this shift. S. lycopersicum M82 SlIPMS3 (Solyc08g014230) encodes a functional end product inhibition-insensitive version of the committing enzyme of leucine biosynthesis, isopropylmalate synthase, missing the carboxyl-terminal 160 amino acids. In contrast, the S. pennellii LA0716 IPMS3 allele found in IL8-1/8-1-1 encodes a nonfunctional truncated IPMS protein. M82 transformed with an SlIPMS3 RNA interference construct exhibited an acylsugar profile similar to that of IL8-1-1, whereas the expression of SlIPMS3 in IL8-1-1 partially restored the M82 acylsugar phenotype. These IPMS3 alleles are polymorphic in 14 S. pennellii accessions spread throughout the geographical range of occurrence for this species and are associated with acylsugars containing varying amounts of 2-methylpropanoic acid and 3-methylbutanoic acid acyl chains.

MAKER-P: a tool kit for the rapid creation, management, and quality control of plant genome annotations.

We have optimized and extended the widely used annotation engine MAKER in order to better support plant genome annotation efforts. New features include better parallelization for large repeat-rich plant genomes, noncoding RNA annotation capabilities, and support for pseudogene identification. We have benchmarked the resulting software tool kit, MAKER-P, using the Arabidopsis (Arabidopsis thaliana) and maize (Zea mays) genomes. Here, we demonstrate the ability of the MAKER-P tool kit to automatically update, extend, and revise the Arabidopsis annotations in light of newly available data and to annotate pseudogenes and noncoding RNAs absent from The Arabidopsis Informatics Resource 10 build. Our results demonstrate that MAKER-P can be used to manage and improve the annotations of even Arabidopsis, perhaps the best-annotated plant genome. We have also installed and benchmarked MAKER-P on the Texas Advanced Computing Center. We show that this public resource can de novo annotate the entire Arabidopsis and maize genomes in less than 3 h and produce annotations of comparable quality to those of the current The Arabidopsis Information Resource 10 and maize V2 annotation builds.

Characteristics and significance of intergenic polyadenylated RNA transcription in Arabidopsis.

The Arabidopsis (Arabidopsis thaliana) genome is the most well-annotated plant genome. However, transcriptome sequencing in Arabidopsis continues to suggest the presence of polyadenylated (polyA) transcripts originating from presumed intergenic regions. It is not clear whether these transcripts represent novel noncoding or protein-coding genes. To understand the nature of intergenic polyA transcription, we first assessed its abundance using multiple messenger RNA sequencing data sets. We found 6,545 intergenic transcribed fragments (ITFs) occupying 3.6% of Arabidopsis intergenic space. In contrast to transcribed fragments that map to protein-coding and RNA genes, most ITFs are significantly shorter, are expressed at significantly lower levels, and tend to be more data set specific. A surprisingly large number of ITFs (32.1%) may be protein coding based on evidence of translation. However, our results indicate that these "translated" ITFs tend to be close to and are likely associated with known genes. To investigate if ITFs are under selection and are functional, we assessed ITF conservation through cross-species as well as within-species comparisons. Our analysis reveals that 237 ITFs, including 49 with translation evidence, are under strong selective constraint and relatively distant from annotated features. These ITFs are likely parts of novel genes. However, the selective pressure imposed on most ITFs is similar to that of randomly selected, untranscribed intergenic sequences. Our findings indicate that despite the prevalence of ITFs, apart from the possibility of genomic contamination, many may be background or noisy transcripts derived from "junk" DNA, whose production may be inherent to the process of transcription and which, on rare occasions, may act as catalysts for the creation of novel genes.

NYUS.2: an automated machine learning prediction model for the large-scale real-time simulation of grapevine freezing tolerance in North America.

Accurate and real-time monitoring of grapevine freezing tolerance is crucial for the sustainability of the grape industry in cool climate viticultural regions. However, on-site data are limited due to the complexity of measurement. Current prediction models underperform under diverse climate conditions, which limits the large-scale deployment of these methods. We combined grapevine freezing tolerance data from multiple regions in North America and generated a predictive model based on hourly temperature-derived features and cultivar features using AutoGluon, an automated machine learning engine. Feature importance was quantified by AutoGluon and SHAP (SHapley Additive exPlanations) value. The final model was evaluated and compared with previous models for its performance under different climate conditions. The final model achieved an overall 1.36°C root-mean-square error during model testing and outperformed two previous models using three test cultivars at all testing regions. Two feature importance quantification methods identified five shared essential features. Detailed analysis of the features indicates that the model has adequately extracted some biological mechanisms during training. The final model, named NYUS.2, was deployed along with two previous models as an R shiny-based application in the 2022-23 dormancy season, enabling large-scale and real-time simulation of grapevine freezing tolerance in North America for the first time.

Phylogenomic analyses across land plants reveals motifs and coexpression patterns useful for functional prediction in the BAHD acyltransferase family.

The BAHD acyltransferase family is one of the largest enzyme families in flowering plants, containing dozens to hundreds of genes in individual genomes. Highly prevalent in angiosperm genomes, members of this family contribute to several pathways in primary and specialized metabolism. In this study, we performed a phylogenomic analysis of the family using 52 genomes across the plant kingdom to gain deeper insights into its functional evolution and enable function prediction. We found that BAHD expansion in land plants was associated with significant changes in various gene features. Using pre-defined BAHD clades, we identified clade expansions in different plant groups. In some groups, these expansions coincided with the prominence of metabolite classes such as anthocyanins (flowering plants) and hydroxycinnamic acid amides (monocots). Clade-wise motif-enrichment analysis revealed that some clades have novel motifs fixed on either the acceptor or the donor side, potentially reflecting historical routes of functional evolution. Co-expression analysis in rice and Arabidopsis further identified BAHDs with similar expression patterns, however, most co-expressed BAHDs belonged to different clades. Comparing BAHD paralogs, we found that gene expression diverges rapidly after duplication, suggesting that sub/neo-functionalization of duplicate genes occurs quickly via expression diversification. Analyzing co-expression patterns in Arabidopsis in conjunction with orthology-based substrate class predictions and metabolic pathway models led to the recovery of metabolic processes of most of the already-characterized BAHDs as well as definition of novel functional predictions for some uncharacterized BAHDs. Overall, this study provides new insights into the evolution of BAHD acyltransferases and sets up a foundation for their functional characterization.