Redirecting tropane alkaloid metabolism reveals pyrrolidine alkaloid diversity in Atropa belladonna.

Plant-specialized metabolism is complex, with frequent examples of highly branched biosynthetic pathways, and shared chemical intermediates. As such, many plant-specialized metabolic networks are poorly characterized. The N-methyl Δ1 -pyrrolinium cation is a simple pyrrolidine alkaloid and precursor of pharmacologically important tropane alkaloids. Silencing of pyrrolidine ketide synthase (AbPyKS) in the roots of Atropa belladonna (Deadly Nightshade) reduces tropane alkaloid abundance and causes high N-methyl Δ1 -pyrrolinium cation accumulation. The consequences of this metabolic shift on alkaloid metabolism are unknown. In this study, we utilized discovery metabolomics coupled with AbPyKS silencing to reveal major changes in the root alkaloid metabolome of A. belladonna. We discovered and annotated almost 40 pyrrolidine alkaloids that increase when AbPyKS activity is reduced. Suppression of phenyllactate biosynthesis, combined with metabolic engineering in planta, and chemical synthesis indicates several of these pyrrolidines share a core structure formed through the nonenzymatic Mannich-like decarboxylative condensation of the N-methyl Δ1 -pyrrolinium cation with 2-O-malonylphenyllactate. Decoration of this core scaffold through hydroxylation and glycosylation leads to mono- and dipyrrolidine alkaloid diversity. This study reveals the previously unknown complexity of the A. belladonna root metabolome and creates a foundation for future investigation into the biosynthesis, function, and potential utility of these novel alkaloids.

Fruity, sticky, stinky, spicy, bitter, addictive, and deadly: evolutionary signatures of metabolic complexity in the Solanaceae.

Covering: 2000-2022Plants collectively synthesize a huge repertoire of metabolites. General metabolites, also referred to as primary metabolites, are conserved across the plant kingdom and are required for processes essential to growth and development. These include amino acids, sugars, lipids, and organic acids. In contrast, specialized metabolites, historically termed secondary metabolites, are structurally diverse, exhibit lineage-specific distribution and provide selective advantage to host species to facilitate reproduction and environmental adaptation. Due to their potent bioactivities, plant specialized metabolites attract considerable attention for use as flavorings, fragrances, pharmaceuticals, and bio-pesticides. The Solanaceae (Nightshade family) consists of approximately 2700 species and includes crops of significant economic, cultural, and scientific importance: these include potato, tomato, pepper, eggplant, tobacco, and petunia. The Solanaceae has emerged as a model family for studying the biochemical evolution of plant specialized metabolism and multiple examples exist of lineage-specific metabolites that influence the senses and physiology of commensal and harmful organisms, including humans. These include, alcohols, phenylpropanoids, and carotenoids that contribute to fruit aroma and color in tomato (fruity), glandular trichome-derived terpenoids and acylsugars that contribute to plant defense (stinky & sticky, respectively), capsaicinoids in chilli-peppers that influence seed dispersal (spicy), and steroidal glycoalkaloids (bitter) from Solanum, nicotine (addictive) from tobacco, as well as tropane alkaloids (deadly) from Deadly Nightshade that deter herbivory. Advances in genomics and metabolomics, coupled with the adoption of comparative phylogenetic approaches, resulted in deeper knowledge of the biosynthesis and evolution of these metabolites. This review highlights recent progress in this area and outlines opportunities for - and challenges of-developing a more comprehensive understanding of Solanaceae metabolism.

Correction: Fruity, sticky, stinky, spicy, bitter, addictive, and deadly: evolutionary signatures of metabolic complexity in the Solanaceae.

Correction for 'Fruity, sticky, stinky, spicy, bitter, addictive, and deadly: evolutionary signatures of metabolic complexity in the Solanaceae' by Paul D. Fiesel et al., Nat. Prod. Rep., 2022, 39, 1438-1464, https://doi.org/10.1039/D2NP00003B.

Optimising the use of gene expression data to predict plant metabolic pathway memberships.

Plant metabolites from diverse pathways are important for plant survival, human nutrition and medicine. The pathway memberships of most plant enzyme genes are unknown. While co-expression is useful for assigning genes to pathways, expression correlation may exist only under specific spatiotemporal and conditional contexts. Utilising > 600 tomato (Solanum lycopersicum) expression data combinations, three strategies for predicting memberships in 85 pathways were explored. Optimal predictions for different pathways require distinct data combinations indicative of pathway functions. Naive prediction (i.e. identifying pathways with the most similarly expressed genes) is error prone. In 52 pathways, unsupervised learning performed better than supervised approaches, possibly due to limited training data availability. Using gene-to-pathway expression similarities led to prediction models that outperformed those based simply on expression levels. Using 36 experimental validated genes, the pathway-best model prediction accuracy is 58.3%, significantly better compared with that for predicting annotated genes without experimental evidence (37.0%) or random guess (1.2%), demonstrating the importance of data quality. Our study highlights the need to extensively explore expression-based features and prediction strategies to maximise the accuracy of metabolic pathway membership assignment. The prediction framework outlined here can be applied to other species and serves as a baseline model for future comparisons.

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.

Characterization of Trichome-Expressed BAHD Acyltransferases in Petunia axillaris Reveals Distinct Acylsugar Assembly Mechanisms within the Solanaceae.

Acylsugars are synthesized in the glandular trichomes of the Solanaceae family and are implicated in protection against abiotic and biotic stress. Acylsugars are composed of either sucrose or glucose esterified with varying numbers of acyl chains of differing length. In tomato (Solanum lycopersicum), acylsugar assembly requires four acylsugar acyltransferases (ASATs) of the BAHD superfamily. Tomato ASATs catalyze the sequential esterification of acyl-coenzyme A thioesters to the R4, R3, R3', and R2 positions of sucrose, yielding a tetra-acylsucrose. Petunia spp. synthesize acylsugars that are structurally distinct from those of tomato. To explore the mechanisms underlying this chemical diversity, a Petuniaaxillaris transcriptome was mined for trichome preferentially expressed BAHDs. A combination of phylogenetic analyses, gene silencing, and biochemical analyses coupled with structural elucidation of metabolites revealed that acylsugar assembly is not conserved between tomato and petunia. In P. axillaris, tetra-acylsucrose assembly occurs through the action of four ASATs, which catalyze sequential addition of acyl groups to the R2, R4, R3, and R6 positions. Notably, in P. axillaris, PaxASAT1 and PaxASAT4 catalyze the acylation of the R2 and R6 positions of sucrose, respectively, and no clear orthologs exist in tomato. Similarly, petunia acylsugars lack an acyl group at the R3' position, and congruently, an ortholog of SlASAT3, which catalyzes acylation at the R3' position in tomato, is absent in P. axillaris Furthermore, where putative orthologous relationships of ASATs are predicted between tomato and petunia, these are not supported by biochemical assays. Overall, these data demonstrate the considerable evolutionary plasticity of acylsugar biosynthesis.

A root-expressed L-phenylalanine:4-hydroxyphenylpyruvate aminotransferase is required for tropane alkaloid biosynthesis in Atropa belladonna.

The tropane alkaloids, hyoscyamine and scopolamine, are medicinal compounds that are the active components of several therapeutics. Hyoscyamine and scopolamine are synthesized in the roots of specific genera of the Solanaceae in a multistep pathway that is only partially elucidated. To facilitate greater understanding of tropane alkaloid biosynthesis, a de novo transcriptome assembly was developed for Deadly Nightshade (Atropa belladonna). Littorine is a key intermediate in hyoscyamine and scopolamine biosynthesis that is produced by the condensation of tropine and phenyllactic acid. Phenyllactic acid is derived from phenylalanine via its transamination to phenylpyruvate, and mining of the transcriptome identified a phylogenetically distinct aromatic amino acid aminotransferase (ArAT), designated Ab-ArAT4, that is coexpressed with known tropane alkaloid biosynthesis genes in the roots of A. belladonna. Silencing of Ab-ArAT4 disrupted synthesis of hyoscyamine and scopolamine through reduction of phenyllactic acid levels. Recombinant Ab-ArAT4 preferentially catalyzes the first step in phenyllactic acid synthesis, the transamination of phenylalanine to phenylpyruvate. However, rather than utilizing the typical keto-acid cosubstrates, 2-oxoglutarate, pyruvate, and oxaloacetate, Ab-ArAT4 possesses strong substrate preference and highest activity with the aromatic keto-acid, 4-hydroxyphenylpyruvate. Thus, Ab-ArAT4 operates at the interface between primary and specialized metabolism, contributing to both tropane alkaloid biosynthesis and the direct conversion of phenylalanine to tyrosine.

Evolution of a complex locus for terpene biosynthesis in solanum.

Functional gene clusters, containing two or more genes encoding different enzymes for the same pathway, are sometimes observed in plant genomes, most often when the genes specify the synthesis of specialized defensive metabolites. Here, we show that a cluster of genes in tomato (Solanum lycopersicum; Solanaceae) contains genes for terpene synthases (TPSs) that specify the synthesis of monoterpenes and diterpenes from cis-prenyl diphosphates, substrates that are synthesized by enzymes encoded by cis-prenyl transferase (CPT) genes also located within the same cluster. The monoterpene synthase genes in the cluster likely evolved from a diterpene synthase gene in the cluster by duplication and divergence. In the orthologous cluster in Solanum habrochaites, a new sesquiterpene synthase gene was created by a duplication event of a monoterpene synthase followed by a localized gene conversion event directed by a diterpene synthase gene. The TPS genes in the Solanum cluster encoding cis-prenyl diphosphate-utilizing enzymes are closely related to a tobacco (Nicotiana tabacum; Solanaceae) diterpene synthase encoding Z-abienol synthase (Nt-ABS). Nt-ABS uses the substrate copal-8-ol diphosphate, which is made from the all-trans geranylgeranyl diphosphate by copal-8-ol diphosphate synthase (Nt-CPS2). The Solanum gene cluster also contains an ortholog of Nt-CPS2, but it appears to encode a nonfunctional protein. Thus, the Solanum functional gene cluster evolved by duplication and divergence of TPS genes, together with alterations in substrate specificity to utilize cis-prenyl diphosphates and through the acquisition of CPT genes.

KNOX genes influence a gradient of fruit chloroplast development through regulation of GOLDEN2-LIKE expression in tomato.

The chlorophyll content of unripe fleshy fruits is positively correlated with the nutrient content and flavor of ripe fruit. In tomato (Solanum lycopersicum) fruit, the uniform ripening (u) locus, which encodes a GOLDEN 2-LIKE transcription factor (SlGLK2), influences a gradient of chloroplast development that extends from the stem end of the fruit surrounding the calyx to the base of the fruit. With the exception of the u locus, the factors that influence the formation of this developmental gradient are unknown. In this study, characterization and positional cloning of the uniform gray-green (ug) locus of tomato reveals a thus far unknown role for the Class I KNOTTED1-LIKE HOMEOBOX (KNOX) gene, TKN4, in specifying the formation of this chloroplast gradient. The involvement of KNOX in fruit chloroplast development was confirmed through characterization of the Curl (Cu) mutant, a dominant gain-of-function mutation of TKN2, which displays ectopic fruit chloroplast development that resembles SlGLK2 over-expression. TKN2 and TKN4 act upstream of SlGLK2 and the related gene ARABIDOPSIS PSEUDO RESPONSE REGULATOR 2-LIKE (SlAPRR2-LIKE) to establish their latitudinal gradient of expression across developing fruit that leads to a gradient of chloroplast development. Class I KNOX genes typically influence plant morphology through maintenance of meristem activity, but this study identifies a role for TKN2 and TKN4 in specifically influencing chloroplast development in fruit but not leaves, suggesting that this fundamental process is differentially regulated in these two organs.

Transcriptome-enabled marker discovery and mapping of plastochron-related genes in Petunia spp.

BACKGROUND: Petunia (Petunia × hybrida), derived from a hybrid between P. axillaris and P. integrifolia, is one of the most economically important bedding plant crops and Petunia spp. serve as model systems for investigating the mechanisms underlying diverse mating systems and pollination syndromes. In addition, we have previously described genetic variation and quantitative trait loci (QTL) related to petunia development rate and morphology, which represent important breeding targets for the floriculture industry to improve crop production and performance. Despite the importance of petunia as a crop, the floriculture industry has been slow to adopt marker assisted selection to facilitate breeding strategies and there remains a limited availability of sequences and molecular markers from the genus compared to other economically important members of the Solanaceae family such as tomato, potato and pepper. RESULTS: Here we report the de novo assembly, annotation and characterization of transcriptomes from P. axillaris, P. exserta and P. integrifolia. Each transcriptome assembly was derived from five tissue libraries (callus, 3-week old seedlings, shoot apices, flowers of mixed developmental stages, and trichomes). A total of 74,573, 54,913, and 104,739 assembled transcripts were recovered from P. axillaris, P. exserta and P. integrifolia, respectively and following removal of multiple isoforms, 32,994 P. axillaris, 30,225 P. exserta, and 33,540 P. integrifolia high quality representative transcripts were extracted for annotation and expression analysis. The transcriptome data was mined for single nucleotide polymorphisms (SNP) and simple sequence repeat (SSR) markers, yielding 89,007 high quality SNPs and 2949 SSRs, respectively. 15,701 SNPs were computationally converted into user-friendly cleaved amplified polymorphic sequence (CAPS) markers and a subset of SNP and CAPS markers were experimentally verified. CAPS markers developed from plastochron-related homologous transcripts from P. axillaris were mapped in an interspecific Petunia population and evaluated for co-localization with QTL for development rate. CONCLUSIONS: The high quality of the three Petunia spp. transcriptomes coupled with the utility of the SNP data will serve as a resource for further exploration of genetic diversity within the genus and will facilitate efforts to develop genetic and physical maps to aid the identification of QTL associated with traits of interest.

Determination of residues responsible for substrate and product specificity of Solanum habrochaites short-chain cis-prenyltransferases.

Isoprenoids are diverse compounds that have their biosynthetic origin in the initial condensation of isopentenyl diphosphate and dimethylallyl diphosphate to form C10 prenyl diphosphates that can be elongated by the addition of subsequent isopentenyl diphosphate units. These reactions are catalyzed by either cis-prenyltransferases (CPTs) or trans-prenyltransferases. The synthesis of volatile terpenes in plants typically proceeds through either geranyl diphosphate (C10) or trans-farnesyl diphosphate (C15), to yield monoterpenes and sesquiterpenes, respectively. However, terpene biosynthesis in glandular trichomes of tomato (Solanum lycopersicum) and related wild relatives also occurs via the cis-substrates neryl diphosphate (NPP) and 2Z,6Z-farnesyl diphosphate (Z,Z-FPP). NPP and Z,Z-FPP are synthesized by neryl diphosphate synthase1 (NDPS1) and Z,Z-farnesyl diphosphate synthase (zFPS), which are encoded by the orthologous CPT1 locus in tomato and Solanum habrochaites, respectively. In this study, comparative sequence analysis of NDPS1 and zFPS enzymes from S. habrochaites accessions that synthesize either monoterpenes or sesquiterpenes was performed to identify amino acid residues that correlate with the ability to synthesize NPP or Z,Z-FPP. Subsequent structural modeling, coupled with site-directed mutagenesis, highlighted the importance of four amino acids located within conserved domain II of CPT enzymes that form part of the second α-helix, for determining substrate and product specificity of these enzymes. In particular, the relative positioning of aromatic amino acid residues at positions 100 and 107 determines the ability of these enzymes to synthesize NPP or Z,Z-FPP. This study provides insight into the biochemical evolution of terpene biosynthesis in the glandular trichomes of Solanum species.

Evolution of TPS20-related terpene synthases influences chemical diversity in the glandular trichomes of the wild tomato relative Solanum habrochaites.

A systematic screen of volatile terpene production in the glandular trichomes of 79 accessions of Solanum habrochaites was conducted and revealed the presence of 21 mono- and sesquiterpenes that exhibit a range of qualitative and quantitative variation. Hierarchical clustering identified distinct terpene phenotypic modules with shared patterns of terpene accumulation across accessions. Several terpene modules could be assigned to previously identified terpene synthase (TPS) activities that included members of the TPS-e/f subfamily that utilize the unusual cis-prenyl diphosphate substrates neryl diphosphate and 2z,6z-farnesyl diphosphate. DNA sequencing and in vitro enzyme activity analysis of TPS-e/f members from S. habrochaites identified three previously unassigned enzyme activities that utilize these cisoid substrates. These produce either the monoterpenes α-pinene and limonene, or the sesquiterpene 7-epizingiberene, with the in vitro analyses that recapitulated the trichome chemistry found in planta. Comparison of the distribution of S. habrochaites accessions with terpene content revealed a strong preference for the presence of particular TPS20 alleles at distinct geographic locations. This study reveals that the unusually high intra-specific variation of volatile terpene synthesis in glandular trichomes of S. habrochaites is due at least in part to evolution at the TPS20 locus.

Differential control of ethylene responses by GREEN-RIPE and GREEN-RIPE LIKE1 provides evidence for distinct ethylene signaling modules in tomato.

The factors that mediate specific responses to the plant hormone ethylene are not fully defined. In particular, it is not known how signaling at the receptor complex can control distinct subsets of ethylene responses. Mutations at the Green-ripe (Gr) and reversion to ethylene sensitivity1 (rte1) loci, which encode homologous proteins of unknown function, influence ethylene responses in tomato (Solanum lycopersicum) and Arabidopsis (Arabidopsis thaliana), respectively. In Arabidopsis, AtRTE1 is required for function of the ETR1 ethylene receptor and acts predominantly through this receptor via direct protein-protein interaction. While most eudicot families including the Brassicaceae possess a single gene that is closely related to AtRTE1, we report that members of the Solanaceae family contain two phylogenetically distinct genes defined by GR and GREEN-RIPE LIKE1 (GRL1), creating the possibility of subfunctionalization. We also show that SlGR and SlGRL1 are differentially expressed in tomato tissues and encode proteins predominantly localized to the Golgi. A combination of overexpression in tomato and complementation of the rte1-3 mutant allele indicates that SlGR and SlGRL1 influence distinct but overlapping ethylene responses. Overexpression of SlGRL1 in the Gr mutant background provides evidence for the existence of different ethylene signaling modules in tomato that are influenced by GR, GRL1, or both. In addition, overexpression of AtRTE1 in tomato leads to reduced ethylene responsiveness in a subset of tissues but does not mimic the Gr mutant phenotype. Together, these data reveal species-specific heterogeneity in the control of ethylene responses mediated by members of the GR/RTE1 family.

Pleiotropic phenotypes of the sticky peel mutant provide new insight into the role of CUTIN DEFICIENT2 in epidermal cell function in tomato.

Plant epidermal cells have evolved specialist functions associated with adaptation to stress. These include the synthesis and deposition of specialized metabolites such as waxes and cutin together with flavonoids and anthocyanins, which have important roles in providing a barrier to water loss and protection against UV radiation, respectively. Characterization of the sticky peel (pe) mutant of tomato (Solanum lycopersicum) revealed several phenotypes indicative of a defect in epidermal cell function, including reduced anthocyanin accumulation, a lower density of glandular trichomes, and an associated reduction in trichome-derived terpenes. In addition, pe mutant fruit are glossy and peels have increased elasticity due to a severe reduction in cutin biosynthesis and altered wax deposition. Leaves of the pe mutant are also cutin deficient and the epicuticular waxes contain a lower proportion of long-chain alkanes. Direct measurements of transpiration, together with chlorophyll-leaching assays, indicate increased cuticular permeability of pe leaves. Genetic mapping revealed that the pe locus represents a new allele of CUTIN DEFICIENT2 (CD2), a member of the class IV homeodomain-leucine zipper gene family, previously only associated with cutin deficiency in tomato fruit. CD2 is preferentially expressed in epidermal cells of tomato stems and is a homolog of Arabidopsis (Arabidopsis thaliana) ANTHOCYANINLESS2 (ANL2). Analysis of cuticle composition in leaves of anl2 revealed that cutin accumulates to approximately 60% of the levels observed in wild-type Arabidopsis. Together, these data provide new insight into the role of CD2 and ANL2 in regulating diverse metabolic pathways and in particular, those associated with epidermal cells.

Altered chloroplast development and delayed fruit ripening caused by mutations in a zinc metalloprotease at the lutescent2 locus of tomato.

The chloroplast is the site of photosynthesis in higher plants but also functions as the center of synthesis for primary and specialized metabolites including amino acids, fatty acids, starch, and diverse isoprenoids. Mutants that disrupt aspects of chloroplast function represent valuable tools for defining structural and biochemical regulation of the chloroplast and its interplay with whole-plant structure and function. The lutescent1 (l1) and l2 mutants of tomato (Solanum lycopersicum) possess a range of chlorophyll-deficient phenotypes including reduced rates of chlorophyll synthesis during deetiolation and enhanced rates of chlorophyll loss in leaves and fruits as they age, particularly in response to high-light stress and darkness. In addition, the onset of fruit ripening is delayed in lutescent mutants by approximately 1 week although once ripening is initiated they ripen at a normal rate and accumulation of carotenoids is not impaired. The l2 locus was mapped to the long arm of chromosome 10 and positional cloning revealed the existence of a premature stop codon in a chloroplast-targeted zinc metalloprotease of the M50 family that is homologous to the Arabidopsis (Arabidopsis thaliana) gene ETHYLENE-DEPENDENT GRAVITROPISM DEFICIENT AND YELLOW-GREEN1. Screening of tomato germplasm identified two additional l2 mutant alleles. This study suggests a role for the chloroplast in mediating the onset of fruit ripening in tomato and indicates that chromoplast development in fruit does not depend on functional chloroplasts.

Striking natural diversity in glandular trichome acylsugar composition is shaped by variation at the Acyltransferase2 locus in the wild tomato Solanum habrochaites.

Acylsugars are polyesters of short- to medium-length acyl chains on sucrose or glucose backbones that are produced in secretory glandular trichomes of many solanaceous plants, including cultivated tomato (Solanum lycopersicum). Despite their roles in biotic stress adaptation and their wide taxonomic distribution, there is relatively little information about the diversity of these compounds and the genes responsible for their biosynthesis. In this study, acylsugar diversity was assessed for 80 accessions of the wild tomato species Solanum habrochaites from throughout the Andes Mountains. Trichome metabolites were analyzed by liquid chromatography-time of flight-mass spectrometry, revealing the presence of at least 34 structurally diverse acylsucroses and two acylglucoses. Distinct phenotypic classes were discovered that varied based on the presence of glucose or sucrose, the numbers and lengths of acyl chains, and the relative total amounts of acylsugars. The presence or absence of an acetyl chain on the acylsucrose hexose ring caused clustering of the accessions into two main groups. Analysis of the Acyltransferase2 gene (the apparent ortholog of Solyc01g105580) revealed differences in enzyme activity and gene expression correlated with polymorphism in S. habrochaites accessions that varied in acylsucrose acetylation. These results are consistent with the hypothesis that glandular trichome acylsugar acetylation is under selective pressure in some populations of S. habrochaites and that the gene mutates to inactivity in the absence of selection.

Metabolomics-guided discovery of cytochrome P450s involved in pseudotropine-dependent biosynthesis of modified tropane alkaloids.

Plant alkaloids constitute an important class of bioactive chemicals with applications in medicine and agriculture. However, the knowledge gap of the diversity and biosynthesis of phytoalkaloids prevents systematic advances in biotechnology for engineered production of these high-value compounds. In particular, the identification of cytochrome P450s driving the structural diversity of phytoalkaloids has remained challenging. Here, we use a combination of reverse genetics with discovery metabolomics and multivariate statistical analysis followed by in planta transient assays to investigate alkaloid diversity and functionally characterize two candidate cytochrome P450s genes from Atropa belladonna without a priori knowledge of their functions or information regarding the identities of key pathway intermediates. This approach uncovered a largely unexplored root localized alkaloid sub-network that relies on pseudotropine as precursor. The two cytochrome P450s catalyze N-demethylation and ring-hydroxylation reactions within the early steps in the biosynthesis of diverse N-demethylated modified tropane alkaloids.

Within- and cross-species predictions of plant specialized metabolism genes using transfer learning.

Plant specialized metabolites mediate interactions between plants and the environment and have significant agronomical/pharmaceutical value. Most genes involved in specialized metabolism (SM) are unknown because of the large number of metabolites and the challenge in differentiating SM genes from general metabolism (GM) genes. Plant models like Arabidopsis thaliana have extensive, experimentally derived annotations, whereas many non-model species do not. Here we employed a machine learning strategy, transfer learning, where knowledge from A. thaliana is transferred to predict gene functions in cultivated tomato with fewer experimentally annotated genes. The first tomato SM/GM prediction model using only tomato data performs well (F-measure = 0.74, compared with 0.5 for random and 1.0 for perfect predictions), but from manually curating 88 SM/GM genes, we found many mis-predicted entries were likely mis-annotated. When the SM/GM prediction models built with A. thaliana data were used to filter out genes where the A. thaliana-based model predictions disagreed with tomato annotations, the new tomato model trained with filtered data improved significantly (F-measure = 0.92). Our study demonstrates that SM/GM genes can be better predicted by leveraging cross-species information. Additionally, our findings provide an example for transfer learning in genomics where knowledge can be transferred from an information-rich species to an information-poor one.

QTL and Transcriptomic Analyses Implicate Cuticle Transcription Factor SHINE as a Source of Natural Variation for Epidermal Traits in Cucumber Fruit.

The fruit surface is a unique tissue with multiple roles influencing fruit development, post-harvest storage and quality, and consumer acceptability. Serving as the first line of protection against herbivores, pathogens, and abiotic stress, the surface can vary markedly among species, cultivars within species, and developmental stage. In this study we explore developmental changes and natural variation of cucumber (Cucumis sativus L.) fruit surface properties using two cucumber lines which vary greatly for these traits and for which draft genomes and a single nucleotide polymorphism (SNP) array are available: Chinese fresh market type, Chinese Long '9930' (CL9930), and pickling type, 'Gy14'. Thin-section samples were prepared from the mid-region of fruit harvested at 0, 4, 8, 12, 16, 20, 24 and 30 days post pollination (dpp), stained with Sudan IV and evaluated for cuticle thickness, depth of wax intercalation between epidermal cells, epidermal cell size and shape, and number and size of lipid droplets. 'Gy14' is characterized by columnar shaped epidermal cells, a 2-3 fold thicker cuticular layer than CL9930, increased cuticular intercalations between cells and a larger number and larger sized lipid droplets. In both lines maximal deposition of cuticle and increase in epidermal size coincided with exponential fruit growth and was largely completed by approximately 16 dpp. Phenotyping and quantitative trait locus mapping (QTL) of fruit sampled from an F7:F8 Gy14 × CL9930 recombinant inbred line (RIL) population identified QTL regions on chromosomes 1, 4 and 5. Strong QTL for epidermal cell height, cuticle thickness, intercalation depth, and diameter of lipid droplets co-localized on chromosome 1. SSR markers on chromosome 1 were used to screen for recombinants in an extended RIL population to refine the QTL region. Further fine mapping by KASP assay combined with gene expression profiling suggested a small number of candidate genes. Tissue specificity, developmental analysis of expression, allelic diversity and gene function implicate the regulatory factor CsSHINE1/WIN1 as a source of natural variation for cucumber fruit epidermal traits.

An Automated Data-Driven Pipeline for Improving Heterologous Enzyme Expression.

Enzymes are the ultimate entities responsible for chemical transformations in natural and engineered biosynthetic pathways. However, many natural enzymes suffer from suboptimal functional expression due to poor intrinsic protein stability. Further, stability enhancing mutations often come at the cost of impaired function. Here we demonstrate an automated protein engineering strategy for stabilizing enzymes while retaining catalytic function using deep mutational scanning coupled to multiple-filter based screening and combinatorial mutagenesis. We validated this strategy by improving the functional expression of a Type III polyketide synthase from the Atropa belladonna biosynthetic pathway for tropane alkaloids. The best variant had a total of 8 mutations with over 25-fold improved activity over wild-type in E. coli cell lysates, an improved melting temperature of 11.5 ± 0.6 °C, and only minimal reduction in catalytic efficiency. We show that the multiple-filter approach maintains acceptable sensitivity with homology modeling structures up to 4 Å RMS. Our results highlight an automated protein engineering tool for improving the stability and solubility of difficult to express enzymes, which has impact for biotechnological applications.