Drivers of metabolic diversification: how dynamic genomic neighbourhoods generate new biosynthetic pathways in the Brassicaceae

Beyond co-expression: pathway discovery for plant pharmaceuticals

Plant natural products have been an important source of medicinal molecules since ancient times. To gain access to the whole diversity of these molecules for pharmaceutical applications, it is important to understand their biosynthetic origins. Whilst co-expression is a reliable tool for identifying gene candidates, a variety of complementary methods can aid in screening or refining candidate selection. Here, we review recently employed plant biosynthetic pathway discovery approaches, and highlight future directions in the field.

Classification, biosynthesis, and biological functions of triterpene esters in plants

Triterpene esters comprise a class of secondary metabolites that are synthesized by decorating triterpene skeletons with a series of oxidation, glycosylation, and acylation modifications. Many triterpene esters with important bioactivities have been isolated and identified, including those with applications in the pesticide, pharmaceutical, and cosmetic industries. They also play essential roles in plant defense against pests, diseases, physical damage (as part of the cuticle), and regulation of root microorganisms. However, there has been no recent summary of the biosynthetic pathways and biological functions of plant triterpene esters. Here, we classify triterpene esters into five categories based on their skeletons and find that C-3 oxidation may have a significant effect on triterpenoid acylation. Fatty acid and aromatic moieties are common ligands present in triterpene esters. We further analyze triterpene ester synthesis-related acyltransferases (TEsACTs) in the triterpene biosynthetic pathway. Using an evolutionary classification of BAHD acyltransferases (BAHD-ATs) and serine carboxypeptidase-like acyltransferases (SCPL-ATs) in Arabidopsis thaliana and Oryza sativa, we classify 18 TEsACTs with identified functions from 11 species. All the triterpene-skeleton-related TEsACTs belong to BAHD-AT clades IIIa and I, and the only identified TEsACT from the SCPL-AT family belongs to the CP-I subfamily. This comprehensive review of the biosynthetic pathways and bioactivities of triterpene esters provides a foundation for further study of their bioactivities and applications in industry, agricultural production, and human health.

Alcohol dehydrogenases regulated by a MYB44 transcription factor underlie Lauraceae citral biosynthesis

Lineage-specific terpenoids have arisen throughout the evolution of land plants and are believed to play a role in interactions between plants and the environment. Species-specific gene clusters in plants have provided insight on the evolution of secondary metabolism. Lauraceae is an ecologically important plant family whose members are also of considerable economic value given their monoterpene contents. However, the gene cluster responsible for the biosynthesis of monoterpenes remains yet to be elucidated. Here, a Lauraceae-specific citral biosynthetic gene cluster (CGC) was identified and investigated using a multifaceted approach that combined phylogenetic, collinearity, and biochemical analyses. The CGC comprises MYB44 as a regulator and 2 alcohol dehydrogenases (ADHs) as modifying enzymes, which derived from species-specific tandem and proximal duplication events. Activity and substrate divergence of the ADHs has resulted in the fruit of mountain pepper (Litsea cubeba), a core Lauraceae species, consisting of more than 80% citral. In addition, MYB44 negatively regulates citral biosynthesis by directly binding to the promoters of the ADH-encoding genes. The aggregation of citral biosynthetic pathways suggests that they may form the basis of important characteristics that enhance adaptability. The findings of this study provide insights into the evolution of and the regulatory mechanisms involved in plant terpene biosynthesis.

The genome of the toxic invasive species Heracleum sosnowskyi carries an increased number of genes despite absence of recent whole-genome duplications

Heracleum sosnowskyi, belonging to a group of giant hogweeds, is a plant with large effects on ecosystems and human health. It is an invasive species that contributes to the deterioration of grassland ecosystems. The ability of H. sosnowskyi to produce linear furanocoumarins (FCs), photosensitizing compounds, makes it very dangerous. At the same time, linear FCs are compounds with high pharmaceutical value used in skin disease therapies. Despite this high importance, it has not been the focus of genetic and genomic studies. Here, we report a chromosome-scale assembly of Sosnowsky's hogweed genome. Genomic analysis revealed an unusually high number of genes (55106) in the hogweed genome, in contrast to the 25-35 thousand found in most plants. However, we did not find any traces of recent whole-genome duplications not shared with its confamiliar, Daucus carota (carrot), which has approximately thirty thousand genes. The analysis of the genomic proximity of duplicated genes indicates on tandem duplications as a main reason for this increase. We performed a genome-wide search of the genes of the FC biosynthesis pathway and surveyed their expression in aboveground plant parts. Using a combination of expression data and phylogenetic analysis, we found candidate genes for psoralen synthase and experimentally showed the activity of one of them using a heterologous yeast expression system. These findings expand our knowledge on the evolution of gene space in plants and lay a foundation for further analysis of hogweed as an invasive plant and as a source of FCs.

Characterization of the horse chestnut genome reveals the evolution of aescin and aesculin biosynthesis

Horse chestnut (Aesculus chinensis) is an important medicinal tree that contains various bioactive compounds, such as aescin, barrigenol-type triterpenoid saponins (BAT), and aesculin, a glycosylated coumarin. Herein, we report a 470.02 Mb genome assembly and characterize an Aesculus-specific whole-genome duplication event, which leads to the formation and duplication of two triterpenoid biosynthesis-related gene clusters (BGCs). We also show that AcOCS6, AcCYP716A278, AcCYP716A275, and AcCSL1 genes within these two BGCs along with a seed-specific expressed AcBAHD6 are responsible for the formation of aescin. Furthermore, we identify seven Aesculus-originated coumarin glycoside biosynthetic genes and achieve the de novo synthesis of aesculin in E. coli. Collinearity analysis shows that the collinear BGC segments can be traced back to early-diverging angiosperms, and the essential gene-encoding enzymes necessary for BAT biosynthesis are recruited before the splitting of Aesculus, Acer, and Xanthoceras. These findings provide insight on the evolution of gene clusters associated with medicinal tree metabolites.

A chromosome-level genome assembly reveals that a bipartite gene cluster formed via an inverted duplication controls monoterpenoid biosynthesis in Schizonepeta tenuifolia

Biosynthetic gene clusters (BGCs) are regions of a genome where genes involved in a biosynthetic pathway are in proximity. The origin and evolution of plant BGCs as well as their role in specialized metabolism remain largely unclear. In this study, we have assembled a chromosome-scale genome of Japanese catnip (Schizonepeta tenuifolia) and discovered a BGC that contains multiple copies of genes involved in four adjacent steps in the biosynthesis of p-menthane monoterpenoids. This BGC has an unprecedented bipartite structure, with mirrored biosynthetic regions separated by 260 kilobases. This bipartite BGC includes identical copies of a gene encoding an old yellow enzyme, a type of flavin-dependent reductase. In vitro assays and virus-induced gene silencing revealed that this gene encodes the missing isopiperitenone reductase. This enzyme evolved from a completely different enzyme family to isopiperitenone reductase from closely related Mentha spp., indicating convergent evolution of this pathway step. Phylogenomic analysis revealed that this bipartite BGC has emerged uniquely in the S. tenuifolia lineage and through insertion of pathway genes into a region rich in monoterpene synthases. The cluster gained its bipartite structure via an inverted duplication. The discovered bipartite BGC for p-menthane biosynthesis in S. tenuifolia has similarities to the recently described duplicated p-menthane biosynthesis gene pairs in the Mentha longifolia genome, providing an example of the convergent evolution of gene order. This work expands our understanding of plant BGCs with respect to both form and evolution, and highlights the power of BGCs for gene discovery in plant biosynthetic pathways.

Analysis of Arabidopsis non-reference accessions reveals high diversity of metabolic gene clusters and discovers new candidate cluster members

Metabolic gene clusters (MGCs) are groups of genes involved in a common biosynthetic pathway. They are frequently formed in dynamic chromosomal regions, which may lead to intraspecies variation and cause phenotypic diversity. We examined copy number variations (CNVs) in four Arabidopsis thaliana MGCs in over one thousand accessions with experimental and bioinformatic approaches. Tirucalladienol and marneral gene clusters showed little variation, and the latter was fixed in the population. Thalianol and especially arabidiol/baruol gene clusters displayed substantial diversity. The compact version of the thalianol gene cluster was predominant and more conserved than the noncontiguous version. In the arabidiol/baruol cluster, we found a large genomic insertion containing divergent duplicates of the CYP705A2 and BARS1 genes. The BARS1 paralog, which we named BARS2, encoded a novel oxidosqualene synthase. The expression of the entire arabidiol/baruol gene cluster was altered in the accessions with the duplication. Moreover, they presented different root growth dynamics and were associated with warmer climates compared to the reference-like accessions. In the entire genome, paired genes encoding terpene synthases and cytochrome P450 oxidases were more variable than their nonpaired counterparts. Our study highlights the role of dynamically evolving MGCs in plant adaptation and phenotypic diversity.

Uncovering a miltiradiene biosynthetic gene cluster in the Lamiaceae reveals a dynamic evolutionary trajectory

The spatial organization of genes within plant genomes can drive evolution of specialized metabolic pathways. Terpenoids are important specialized metabolites in plants with diverse adaptive functions that enable environmental interactions. Here, we report the genome assemblies of Prunella vulgaris, Plectranthus barbatus, and Leonotis leonurus. We investigate the origin and subsequent evolution of a diterpenoid biosynthetic gene cluster (BGC) together with other seven species within the Lamiaceae (mint) family. Based on core genes found in the BGCs of all species examined across the Lamiaceae, we predict a simplified version of this cluster evolved in an early Lamiaceae ancestor. The current composition of the extant BGCs highlights the dynamic nature of its evolution. We elucidate the terpene backbones generated by the Callicarpa americana BGC enzymes, including miltiradiene and the terpene (+)-kaurene, and show oxidization activities of BGC cytochrome P450s. Our work reveals the fluid nature of BGC assembly and the importance of genome structure in contributing to the origin of metabolites.

Metabolomics-centered mining of plant metabolic diversity and function: Past decade and future perspectives

Plants are natural experts in organic synthesis, being able to generate large numbers of specific metabolites with widely varying structures that help them adapt to variable survival challenges. Metabolomics is a research discipline that integrates the capabilities of several types of research including analytical chemistry, statistics, and biochemistry. Its ongoing development provides strategies for gaining a systematic understanding of quantitative changes in the levels of metabolites. Metabolomics is usually performed by targeting either a specific cell, a specific tissue, or the entire organism. Considerable advances in science and technology over the last three decades have propelled us into the era of multi-omics, in which metabolomics, despite at an earlier developmental stage than genomics, transcriptomics, and proteomics, offers the distinct advantage of studying the cellular entities that have the greatest influence on end phenotype. Here, we summarize the state of the art of metabolite detection and identification, and illustrate these techniques with four case study applications: (i) comparing metabolite composition within and between species, (ii) assessing spatio-temporal metabolic changes during plant development, (iii) mining characteristic metabolites of plants in different ecological environments and upon exposure to various stresses, and (iv) assessing the performance of metabolomics as a means of functional gene identification , metabolic pathway elucidation, and metabolomics-assisted breeding through analyzing plant populations with diverse genetic variations. In addition, we highlight the prominent contributions of joint analyses of plant metabolomics and other omics datasets, including those from genomics, transcriptomics, proteomics, epigenomics, phenomics, microbiomes, and ion-omics studies. Finally, we discuss future directions and challenges exploiting metabolomics-centered approaches in understanding plant metabolic diversity.

The Chromosome-level genome of Aesculus wilsonii provides new insights into terpenoid biosynthesis and Aesculus evolution

Aesculus L. (buckeye and horse chestnut) are woody plant species with important horticultural and medicinal values. Aesculus seeds are widely used as biomedicine and cosmetic ingredients due to their saponins. We report a chromosomal-scale genome of Aesculus wilsonii. Sequences amounting to a total of 579.01 Mb were assembled into 20 chromosomes. More than half of the genome (54.46%) were annotated as repetitive sequences, and 46,914 protein-coding genes were predicted. In addition to the widespread gamma event with core eudicots, a unique whole-genome duplication (WGD) event (17.69 Mya) occurred in Aesculus after buckeye differentiated from longan. Due to WGD events and tandem duplications, the related synthetic genes of triterpene saponins unique to Aesculus increased significantly. Combined with transcriptome characterization, the study preliminarily resolved the biosynthetic pathway of triterpenoid saponins like aescin in A. wilsonii genome. Analyses of the resequencing of 104 buckeye accessions revealed clear relationship between the geographic distribution and genetic differentiation of buckeye trees in China. We found that the buckeye species found in southern Shaanxi is A. wilsonii rather than A. chinensis. Population dynamics analysis further suggests that the population size and evolution of existing buckeye species have been influenced by climate fluctuations during the Pleistocene and recent domestication events. The genome of A. wilsonii and population genomics of Aesculus provide a resource for future research on Hippocastanaceae. These findings will contribute to the utilization and diversity protection of Aesculus.

Plant metabolic gene clusters in the multi-omics era

Secondary metabolism in plants gives rise to a vast array of small-molecule natural products. The discovery of operon-like gene clusters in plants has provided a new perspective on the evolution of specialized metabolism and the opportunity to rapidly advance the metabolic engineering of natural product production. Here, we review historical aspects of the study of plant metabolic gene clusters as well as general strategies for identifying plant metabolic gene clusters in the multi-omics era. We also emphasize the exploration of their natural variation and evolution, as well as new strategies for the prospecting of plant metabolic gene clusters and a deeper understanding of how their structure influences their function.

Integrative omics approaches for biosynthetic pathway discovery in plants

Covering: up to 2022With the emergence of large amounts of omics data, computational approaches for the identification of plant natural product biosynthetic pathways and their genetic regulation have become increasingly important. While genomes provide clues regarding functional associations between genes based on gene clustering, metabolome mining provides a foundational technology to chart natural product structural diversity in plants, and transcriptomics has been successfully used to identify new members of their biosynthetic pathways based on coexpression. Thus far, most approaches utilizing transcriptomics and metabolomics have been targeted towards specific pathways and use one type of omics data at a time. Recent technological advances now provide new opportunities for integration of multiple omics types and untargeted pathway discovery. Here, we review advances in plant biosynthetic pathway discovery using genomics, transcriptomics, and metabolomics, as well as recent efforts towards omics integration. We highlight how transcriptomics and metabolomics provide complementary information to link genes to metabolites, by associating temporal and spatial gene expression levels with metabolite abundance levels across samples, and by matching mass-spectral features to enzyme families. Furthermore, we suggest that elucidation of gene regulatory networks using time-series data may prove useful for efforts to unwire the complexities of biosynthetic pathway components based on regulatory interactions and events.

Plant biosynthetic gene clusters in the context of metabolic evolution

Covering: up to 2022Plants produce a wide range of structurally and biosynthetically diverse natural products to interact with their environment. These specialised metabolites typically evolve in limited taxonomic groups presumably in response to specific selective pressures. With the increasing availability of sequencing data, it has become apparent that in many cases the genes encoding biosynthetic enzymes for specialised metabolic pathways are not randomly distributed on the genome. Instead they are physically linked in structures such as arrays, pairs and clusters. The exact function of these clusters is debated. In this review we take a broad view of gene arrangement in plant specialised metabolism, examining types of structures and variation. We discuss the evolution of biosynthetic gene clusters in the wider context of metabolism, populations and epigenetics. Finally, we synthesise our observations to propose a new hypothesis for biosynthetic gene cluster formation in plants.

Convergent gene clusters underpin hyperforin biosynthesis in St John's wort

The meroterpenoid hyperforin is responsible for the antidepressant activity of St John's wort extracts, but the genes controlling its biosynthesis are unknown. Using genome mining and biochemical work, we characterize two biosynthetic gene clusters (BGCs) that encode the first three steps in the biosynthesis of hyperforin precursors. The findings of syntenic and phylogenetic analyses reveal the parallel assembly of the two BGCs. The syntenous BGC in Mesua ferrea indicates that the first cluster was assembled before the divergence of the Hypericaceae and Calophyllaceae families. The assembly of the second cluster is the result of a coalescence of genomic fragments after a major duplication event. The differences between the two BGCs - in terms of gene expression, response to methyl jasmonate, substrate specificity and subcellular localization of key enzymes - suggest that the presence of the two clusters could serve to generate separate pools of precursors. The parallel assembly of two BGCs with similar compositions in a single plant species is uncommon, and our work provides insights into how and when these gene clusters form. Our discovery helps to advance our understanding of the evolution of plant specialized metabolism and its genomic organization. Additionally, our results offer a foundation from which hyperforin biosynthesis can be more fully understood, and which can be used in future metabolic engineering applications.

Genomic Analysis Based on Chromosome-Level Genome Assembly Reveals an Expansion of Terpene Biosynthesis of Azadirachta indica

Azadirachta indica (neem), an evergreen tree of the Meliaceae family, is a source of the potent biopesticide azadirachtin. The lack of a chromosome-level assembly impedes an in-depth understanding of its genome architecture and the comparative genomic analysis of A. indica. Here, a high-quality genome assembly of A. indica was constructed using a combination of data from Illumina, PacBio, and Hi-C technology, which is the first chromosome-scale genome assembly of A. indica. Based on the length of our assembly, the genome size of A. indica is estimated to be 281 Mb anchored to 14 chromosomes (contig N50 = 6 Mb and scaffold N50 = 19 Mb). The genome assembly contained 115 Mb repetitive elements and 25,767 protein-coding genes. Evolutional analysis revealed that A. indica didn't experience any whole-genome duplication (WGD) event after the core eudicot γ event, but some genes and genome segment might likely experienced recent duplications. The secondary metabolite clusters, TPS genes, and CYP genes were also identified. Comparative genomic analysis revealed that most of the A. indica-specific TPS genes and CYP genes were located on the terpene-related clusters on chromosome 13. It is suggested that chromosome 13 may play an important role in the specific terpene biosynthesis of A. indica. The gene duplication events may be responsible for the terpene biosynthesis expansion in A. indica. The genomic dataset and genomic analysis created for A. indica will shed light on terpene biosynthesis in A. indica and facilitate comparative genomic research of the family Meliaceae.

Integrated Metabolo-transcriptomics Reveals the Defense Response of Homogentisic Acid in Wheat against Puccinia striiformis f. sp. tritici

Stripe rust is a widespread and harmful wheat disease caused by Puccinia striiformis f. sp. tritici (Pst) worldwide. Targeted metabolome and transcriptomics analyses of CYR23 infected leaves were performed to identify the differential metabolites and differentially expressed genes related to wheat disease resistance. We observed upregulation of 33 metabolites involved in the primary and secondary metabolism, especially for homogentisic acid (HGA), p-coumaroylagmatine, and saccharopine. These three metabolites were mainly involved in the phenylpropanoid metabolic pathway, hydroxycinnamic acid amides pathway, and saccharopine pathway. Combined with transcriptome data on non-compatible interaction, the synthesis-related genes of these three differential metabolites were all upregulated significantly. The gene regulatory network involved in response to Pst infection was constructed, which revealed that several transcription factor families including WRKYs, MYBs, and bZIPs were identified as potentially hubs in wheat resistance response against Pst. An in vitro test showed that HGA effectively inhibited the germination of stripe rust fungus urediniospores and reduced the occurrence of wheat stripe rust. The results of gene silencing and overexpression of HGA synthesis-related gene 4-hydroxyphenylpyruvate dioxygenase proved that HGA was involved in wheat disease resistance. These results provided a further understanding of the disease resistance of wheat and indicated that HGA can be developed as a potential agent against Pst.

Unlocking plant metabolic diversity: A (pan)-genomic view

Plants produce a remarkable diversity of structurally and functionally diverse natural chemicals that serve as adaptive compounds throughout their life cycles. However, unlocking this metabolic diversity is significantly impeded by the size, complexity, and abundant repetitive elements of typical plant genomes. As genome sequencing becomes routine, we anticipate that links between metabolic diversity and genetic variation will be strengthened. In addition, an ever-increasing number of plant genomes have revealed that biosynthetic gene clusters are not only a hallmark of microbes and fungi; gene clusters for various classes of compounds have also been found in plants, and many are associated with important agronomic traits. We present recent examples of plant metabolic diversification that have been discovered through the exploration and exploitation of various genomic and pan-genomic data. We also draw attention to the fundamental genomic and pan-genomic basis of plant chemodiversity and discuss challenges and future perspectives for investigating metabolic diversity in the coming pan-genomics era.

Probing the functions of friedelane-type triterpene cyclases from four celastrol-producing plants

Triterpenes are among the most diverse plant natural products, and their diversity is closely related to various triterpene skeletons catalyzed by different 2,3-oxidosqualene cyclases (OSCs). Celastrol, a friedelane-type triterpene with significant bioactivities, is specifically distributed in higher plants, such as Celastraceae species. Friedelin is an important precursor for the biosynthesis of celastrol, and it is synthesized through the cyclization of 2,3-oxidosqualene, with the highest number of rearrangements being catalyzed by friedelane-type triterpene cyclases. However, the molecular mechanisms underlying the catalysis of friedelin production by friedelane-type triterpene cyclases have not yet been fully elucidated. In this study, transcriptome data of four celastrol-producing plants from Celastraceae were used to identify a total of 21 putative OSCs. Through functional characterization, the friedelane-type triterpene cyclases were separately verified in the four plants. Analysis of the selection pressure showed that purifying selection acted on these OSCs, and the friedelane-type triterpene cyclases may undergo weaker selective restriction during evolution. Molecular docking and site-directed mutagenesis revealed that changes in some amino acids that are unique to friedelane-type triterpene cyclases may lead to variations in catalytic specificity or efficiency, thereby affecting the synthesis of friedelin. Our research explored the functional diversity of triterpene synthases from a multispecies perspective. It also provides some references for further research on the relative mechanisms of friedelin biosynthesis.

Gene cluster from plant to microbes: Their role in genome architecture, organism's development, specialized metabolism and drug discovery

Plants and microbes fulfil our daily requirements through different high-value chemicals, e.g., nutraceuticals, pharmaceuticals, cosmetics, and through varieties of fruits, crops, vegetables, and many more. Utmost care would therefore be taken for growth, development and sustainability of these important crops and medicinal plants and microbes. Homeobox genes and HOX clusters and their recently characterized expanded family members, including newly discovered homeobox, WOX gene from medicinal herb, Panax ginseng, significantly contributes in the growth and development of these organisms. On the other hand, secondary metabolites produced through secondary metabolism of plants and microbes are used as organisms defense as well as drugs/drug-like molecules for humans. Both the developmental HOX cluster and the biosynthetic gene-cluster (BGC) for secondary metabolites are organised in organisms genome. Genome mining and genomewide analysis of these clusters will definitely identify and characterize many more important molecules from unexplored plants and microbes and underexplored human microbiota and the evolution studies of these clusters will indicate their source of origin. Although genomics revolution now continues at a pace, till date only few hundred plant genome sequences are available. However, next-generation sequencing (NGS) technology now in market and may be applied even for plants with recalcitrant genomes, eventually may discover genomic potential towards production of secondary metabolites of diverse plants and micro-organisms present in the environment and microbiota. Additionally, the development of tools for genome mining e.g., antiSMASH, plantiSMASH, and more and more computational approaches that predicts hundreds of secondary metabolite BGCs will be discussed.

The confluence of big data and evolutionary genome mining for the discovery of natural products

This review covers literature between 2003-2021The development and application of genome mining tools has given rise to ever-growing genetic and chemical databases and propelled natural products research into the modern age of Big Data. Likewise, an explosion of evolutionary studies has unveiled genetic patterns of natural products biosynthesis and function that support Darwin's theory of natural selection and other theories of adaptation and diversification. In this review, we aim to highlight how Big Data and evolutionary thinking converge in the study of natural products, and how this has led to an emerging sub-discipline of evolutionary genome mining of natural products. First, we outline general principles to best utilize Big Data in natural products research, addressing key considerations needed to provide evolutionary context. We then highlight successful examples where Big Data and evolutionary analyses have been combined to provide bioinformatic resources and tools for the discovery of novel natural products and their biosynthetic enzymes. Rather than an exhaustive list of evolution-driven discoveries, we highlight examples where Big Data and evolutionary thinking have been embraced for the evolutionary genome mining of natural products. After reviewing the nascent history of this sub-discipline, we discuss the challenges and opportunities of genomic and metabolomic tools with evolutionary foundations and/or implications and provide a future outlook for this emerging and exciting field of natural product research.

Diverse triterpene skeletons are derived from the expansion and divergent evolution of 2,3-oxidosqualene cyclases in plants

Triterpenoids are one of the largest groups of secondary metabolites and exhibit diverse structures, which are derived from C30 skeletons that are biosynthesized via the isoprenoid pathway by cyclization of 2,3-oxidosqualene. Triterpenoids have a wide range of biological activities, and are used in functional foods, drugs, and as industrial materials. Due to the low content levels in their native plants and limited feasibility and efficiency of chemical synthesis, heterologous biosynthesis of triterpenoids is the most promising strategy. Herein, we classified 121 triterpene alcohols/ketones according to their conformation and ring numbers, among which 51 skeletons have been experimentally characterized as the products of oxidosqualene cyclases (OSCs). Interestingly, 24 skeletons that have not been reported from nature source were generated by OSCs in heterologous expression. Comprehensive evolutionary analysis of the identified 152 OSCs from 75 species in 25 plant orders show that several pentacyclic triterpene synthases repeatedly originated in multiple plant lineages. Comparative analysis of OSC catalytic reaction revealed that stabilization of intermediate cations, steric hindrance, and conformation of active center amino acid residues are primary factors affecting triterpene formation. Optimization of OSC could be achieved by changing of side-chain orientations of key residues. Recently, methods, such as rationally design of pathways, regulation of metabolic flow, compartmentalization engineering, etc., were introduced in improving chassis for the biosynthesis of triterpenoids. We expect that extensive study of natural variation of large number of OSCs and catalytical mechanism will provide basis for production of high level of triterpenoids by application of synthetic biology strategies.

BAHD acetyltransferase contributes to wound-induced biosynthesis of oleo-gum resin triterpenes in Boswellia

Triterpenes (30-carbon isoprene compounds) represent a large and highly diverse class of natural products that play various physiological functions in plants. The triterpene biosynthetic enzymes, particularly those catalyzing the late-stage regio-selective modifications are not well characterized. The bark of select Boswellia trees, e.g., B. serrata exudes specialized oleo-gum resin in response to wounding, which is enriched with boswellic acids (BAs), a unique class of C3α-epimeric pentacyclic triterpenes with medicinal properties. The bark possesses a network of resin secretory structures comprised of vertical and horizontal resin canals, and amount of BAs in bark increases considerably in response to wounding. To investigate BA biosynthetic enzymes, we conducted tissue-specific transcriptome profiling and identified a wound-responsive BAHD acetyltransferase (BsAT1) of B. serrata catalyzing the late-stage C3α-O-acetylation reactions in the BA biosynthetic pathway. BsAT1 catalyzed C3α-O-acetylation of αBA, βBA, and 11-keto-βBA in vitro and in planta assays to produce all the major C3α-O-acetyl-BAs (3-acetyl-αBA, 3-acetyl-βBA, and 3-acetyl-11-keto-βBA) found in B. serrata bark and oleo-gum resin. BsAT1 showed strict specificity for BA scaffold, whereas it did not acetylate the more common C3β-epimeric pentacyclic triterpenes. The analysis of steady-state kinetics using various BAs revealed distinct substrate affinity and catalytic efficiency. BsAT1 transcript expression coincides with increased levels of C3α-O-acetyl-BAs in bark in response to wounding, suggesting a role of BsAT1 in wound-induced biosynthesis of C3α-O-acetyl-BAs. Overall, the results provide new insights into the biosynthesis of principal chemical constituents of Boswellia oleo-gum resin.

Mining genomes to illuminate the specialized chemistry of life

All organisms produce specialized organic molecules, ranging from small volatile chemicals to large gene-encoded peptides, that have evolved to provide them with diverse cellular and ecological functions. As natural products, they are broadly applied in medicine, agriculture and nutrition. The rapid accumulation of genomic information has revealed that the metabolic capacity of virtually all organisms is vastly underappreciated. Pioneered mainly in bacteria and fungi, genome mining technologies are accelerating metabolite discovery. Recent efforts are now being expanded to all life forms, including protists, plants and animals, and new integrative omics technologies are enabling the increasingly effective mining of this molecular diversity.

Plant Metabolic Gene Clusters: Evolution, Organization, and Their Applications in Synthetic Biology

Plants are a remarkable source of high-value specialized metabolites having significant physiological and ecological functions. Genes responsible for synthesizing specialized metabolites are often clustered together for a coordinated expression, which is commonly observed in bacteria and filamentous fungi. Similar to prokaryotic gene clustering, plants do have gene clusters encoding enzymes involved in the biosynthesis of specialized metabolites. More than 20 gene clusters involved in the biosynthesis of diverse metabolites have been identified across the plant kingdom. Recent studies demonstrate that gene clusters are evolved through gene duplications and neofunctionalization of primary metabolic pathway genes. Often, these clusters are tightly regulated at nucleosome level. The prevalence of gene clusters related to specialized metabolites offers an attractive possibility of an untapped source of highly useful biomolecules. Accordingly, the identification and functional characterization of novel biosynthetic pathways in plants need to be worked out. In this review, we summarize insights into the evolution of gene clusters and discuss the organization and importance of specific gene clusters in the biosynthesis of specialized metabolites. Regulatory mechanisms which operate in some of the important gene clusters have also been briefly described. Finally, we highlight the importance of gene clusters to develop future metabolic engineering or synthetic biology strategies for the heterologous production of novel metabolites.

cblaster: a remote search tool for rapid identification and visualization of homologous gene clusters

Motivation

Genes involved in coordinated biological pathways, including metabolism, drug resistance and virulence, are often collocalized as gene clusters. Identifying homologous gene clusters aids in the study of their function and evolution, however, existing tools are limited to searching local sequence databases. Tools for remotely searching public databases are necessary to keep pace with the rapid growth of online genomic data.

Results

Here, we present cblaster, a Python-based tool to rapidly detect collocated genes in local and remote databases. cblaster is easy to use, offering both a command line and a user-friendly graphical user interface. It generates outputs that enable intuitive visualizations of large datasets and can be readily incorporated into larger bioinformatic pipelines. cblaster is a significant update to the comparative genomics toolbox.

Availability and implementation

cblaster source code and documentation is freely available from GitHub under the MIT license (github.com/gamcil/cblaster).

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Origin and Evolution of the Cannabinoid Oxidocyclase Gene Family

Cannabis is an ancient crop representing a rapidly increasing legal market, especially for medicinal purposes. Medicinal and psychoactive effects of Cannabis rely on specific terpenophenolic ligands named cannabinoids. Recent whole-genome sequencing efforts have uncovered variation in multiple genes encoding the final steps in cannabinoid biosynthesis. However, the origin, evolution, and phylogenetic relationships of these cannabinoid oxidocyclase genes remain unclear. To elucidate these aspects, we performed comparative genomic analyses of Cannabis, related genera within the Cannabaceae family, and selected outgroup species. Results show that cannabinoid oxidocyclase genes originated in the Cannabis lineage from within a larger gene expansion in the Cannabaceae family. Localization and divergence of oxidocyclase genes in the Cannabis genome revealed two main syntenic blocks, each comprising tandemly repeated cannabinoid oxidocyclase genes. By comparing these blocks with those in genomes from closely related species, we propose an evolutionary model for the origin, neofunctionalization, duplication, and diversification of cannabinoid oxidocycloase genes. Based on phylogenetic analyses, we propose a comprehensive classification of three main clades and seven subclades that are intended to aid unequivocal referencing and identification of cannabinoid oxidocyclase genes. Our data suggest that cannabinoid phenotype is primarily determined by the presence/absence of single-copy genes. Although wild populations of Cannabis are still unknown, increased sampling of landraces and wild/feral populations across its native geographic range is likely to uncover additional cannabinoid oxidocyclase sequence variants.

Unraveling the roles of plant specialized metabolites: using synthetic biology to design molecular biosensors

Plants are a rich source of specialized metabolites with a broad range of bioactivities and many applications in human daily life. Over the past decades significant progress has been made in identifying many such metabolites in different plant species and in elucidating their biosynthetic pathways. However, the biological roles of plant specialized metabolites remain elusive and proposed functions lack an identified underlying molecular mechanism. Understanding the roles of specialized metabolites frequently is hampered by their dynamic production and their specific spatiotemporal accumulation within plant tissues and organs throughout a plant's life cycle. In this review, we propose the employment of strategies from the field of Synthetic Biology to construct and optimize genetically encoded biosensors that can detect individual specialized metabolites in a standardized and high-throughput manner. This will help determine the precise localization of specialized metabolites at the tissue and single-cell levels. Such information will be useful in developing complete system-level models of specialized plant metabolism, which ultimately will demonstrate how the biosynthesis of specialized metabolites is integrated with the core processes of plant growth and development.

An independent evolutionary origin for insect deterrent cucurbitacins in Iberis amara

Pieris rapae and Phyllotreta nemorum are Brassicaceae specialists, but do not feed on Iberis amara spp. that contain cucurbitacins. The cucurbitacins are highly oxygenated triterpenoid, occurring widespread in cucurbitaceous species and in a few other plant families. Using de-novo assembled transcriptomics from I. amara, gene co-expression analysis and comparative genomics, we unraveled the evolutionary origin of the insect deterrent cucurbitacins in I. amara. Phylogenetic analysis of five oxidosqualene cyclases and heterologous expression allowed us to identify the first committed enzyme in cucurbitacin biosynthesis in I. amara, cucurbitadienol synthase (IaCPQ). In addition, two species-specific cytochrome P450s (CYP708A16 and CYP708A15) were identified that catalyse the unique C16 and C22 hydroxylation of the cucurbitadienol backbone, enzymatic steps that have not been reported before. Furthermore, the draft genome assembly of I. amara showed that the IaCPQ was localized to the same scaffold together with CYP708A15 but spanning over 100 kb, this contrasts with the highly organized cucurbitacin gene cluster in the cucurbits. These results reveal that cucurbitacin biosynthesis has evolved convergently via different biosynthetic routes in different families rather than through divergence from an ancestral pathway. This study thus provides new insight into the mechanism of recurrent evolution and diversification of a plant defensive chemical.

Polyploidy underlies co-option and diversification of biosynthetic triterpene pathways in the apple tribe

Whole-genome duplication (WGD) plays important roles in plant evolution and function, yet little is known about how WGD underlies metabolic diversification of natural products that bear significant medicinal properties, especially in nonmodel trees. Here, we reveal how WGD laid the foundation for co-option and differentiation of medicinally important ursane triterpene pathway duplicates, generating distinct chemotypes between species and between developmental stages in the apple tribe. After generating chromosome-level assemblies of a widely cultivated loquat variety and Gillenia trifoliata, we define differentially evolved, duplicated gene pathways and date the WGD in the apple tribe at 13.5 to 27.1 Mya, much more recent than previously thought. We then functionally characterize contrasting metabolic pathways responsible for major triterpene biosynthesis in G. trifoliata and loquat, which pre- and postdate the Maleae WGD, respectively. Our work mechanistically details the metabolic diversity that arose post-WGD and provides insights into the genomic basis of medicinal properties of loquat, which has been used in both traditional and modern medicines.

Expansion within the CYP71D subfamily drives the heterocyclization of tanshinones synthesis in Salvia miltiorrhiza

Tanshinones are the bioactive nor-diterpenoid constituents of the Chinese medicinal herb Danshen (Salvia miltiorrhiza). These groups of chemicals have the characteristic furan D-ring, which differentiates them from the phenolic abietane-type diterpenoids frequently found in the Lamiaceae family. However, how the 14,16-epoxy is formed has not been elucidated. Here, we report an improved genome assembly of Danshen using a highly homozygous genotype. We identify a cytochrome P450 (CYP71D) tandem gene array through gene expansion analysis. We show that CYP71D373 and CYP71D375 catalyze hydroxylation at carbon-16 (C16) and 14,16-ether (hetero)cyclization to form the D-ring, whereas CYP71D411 catalyzes upstream hydroxylation at C20. In addition, we discover a large biosynthetic gene cluster associated with tanshinone production. Collinearity analysis indicates a more specific origin of tanshinones in Salvia genus. It illustrates the evolutionary origin of abietane-type diterpenoids and those with a furan D-ring in Lamiaceae.

Nature's Chemists: The Discovery and Engineering of Phytochemical Biosynthesis

Plants produce a diverse array of natural products, many of which have high pharmaceutical value or therapeutic potential. However, these compounds often occur at low concentrations in uncultivated species. Producing phytochemicals in heterologous systems has the potential to address the bioavailability issues related to obtaining these molecules from their natural source. Plants are suitable heterologous systems for the production of valuable phytochemicals as they are autotrophic, derive energy and carbon from photosynthesis, and have similar cellular context to native producer plants. In this review we highlight the methods that are used to elucidate natural product biosynthetic pathways, including the approaches leading to proposing the sequence of enzymatic steps, selecting enzyme candidates and characterizing gene function. We will also discuss the advantages of using plant chasses as production platforms for high value phytochemicals. In addition, through this report we will assess the emerging metabolic engineering strategies that have been developed to enhance and optimize the production of natural and novel bioactive phytochemicals in heterologous plant systems.

Formation and diversification of a paradigm biosynthetic gene cluster in plants

Numerous examples of biosynthetic gene clusters (BGCs), including for compounds of agricultural and medicinal importance, have now been discovered in plant genomes. However, little is known about how these complex traits are assembled and diversified. Here, we examine a large number of variants within and between species for a paradigm BGC (the thalianol cluster), which has evolved recently in a common ancestor of the Arabidopsis genus. Comparisons at the species level reveal differences in BGC organization and involvement of auxiliary genes, resulting in production of species-specific triterpenes. Within species, the thalianol cluster is primarily fixed, showing a low frequency of deleterious haplotypes. We further identify chromosomal inversion as a molecular mechanism that may shuffle more distant genes into the cluster, so enabling cluster compaction. Antagonistic natural selection pressures are likely involved in shaping the occurrence and maintenance of this BGC. Our work sheds light on the birth, life and death of complex genetic and metabolic traits in plants.

Generation of a chromosome-scale genome assembly of the insect-repellent terpenoid-producing Lamiaceae species, Callicarpa americana

BACKGROUND

Plants exhibit wide chemical diversity due to the production of specialized metabolites that function as pollinator attractants, defensive compounds, and signaling molecules. Lamiaceae (mints) are known for their chemodiversity and have been cultivated for use as culinary herbs, as well as sources of insect repellents, health-promoting compounds, and fragrance.

FINDINGS

We report the chromosome-scale genome assembly of Callicarpa americana L. (American beautyberry), a species within the early-diverging Callicarpoideae clade of Lamiaceae, known for its metallic purple fruits and use as an insect repellent due to its production of terpenoids. Using long-read sequencing and Hi-C scaffolding, we generated a 506.1-Mb assembly spanning 17 pseudomolecules with N50 contig and N50 scaffold sizes of 7.5 and 29.0 Mb, respectively. In all, 32,164 genes were annotated, including 53 candidate terpene synthases and 47 putative clusters of specialized metabolite biosynthetic pathways. Our analyses revealed 3 putative whole-genome duplication events, which, together with local tandem duplications, contributed to gene family expansion of terpene synthases. Kolavenyl diphosphate is a gateway to many of the bioactive terpenoids in C. americana; experimental validation confirmed that CamTPS2 encodes kolavenyl diphosphate synthase. Syntenic analyses with Tectona grandis L. f. (teak), a member of the Tectonoideae clade of Lamiaceae known for exceptionally strong wood resistant to insects, revealed 963 collinear blocks and 21,297 C. americana syntelogs.

CONCLUSIONS

Access to the C. americana genome provides a road map for rapid discovery of genes encoding plant-derived agrichemicals and a key resource for understanding the evolution of chemical diversity in Lamiaceae.

Evolution of a plant gene cluster in Solanaceae and emergence of metabolic diversity

Plants produce phylogenetically and spatially restricted, as well as structurally diverse specialized metabolites via multistep metabolic pathways. Hallmarks of specialized metabolic evolution include enzymatic promiscuity and recruitment of primary metabolic enzymes and examples of genomic clustering of pathway genes. Solanaceae glandular trichomes produce defensive acylsugars, with sidechains that vary in length across the family. We describe a tomato gene cluster on chromosome 7 involved in medium chain acylsugar accumulation due to trichome specific acyl-CoA synthetase and enoyl-CoA hydratase genes. This cluster co-localizes with a tomato steroidal alkaloid gene cluster and is syntenic to a chromosome 12 region containing another acylsugar pathway gene. We reconstructed the evolutionary events leading to this gene cluster and found that its phylogenetic distribution correlates with medium chain acylsugar accumulation across the Solanaceae. This work reveals insights into the dynamics behind gene cluster evolution and cell-type specific metabolite diversity.

Co-Regulation of Clustered and Neo-Functionalized Genes in Plant-Specialized Metabolism

Current findings of neighboring genes involved in plant specialized metabolism provide the genomic signatures of metabolic evolution. Two such genomic features, namely, (i) metabolic gene cluster and (ii) neo-functionalization of tandem gene duplications, represent key factors corresponding to the creation of metabolic diversity of plant specialized metabolism. So far, several terpenoid and alkaloid biosynthetic genes have been characterized with gene clusters in some plants. On the other hand, some modification genes involved in flavonoid and glucosinolate biosynthesis were found to arise via gene neo-functionalization. Although the occurrence of both types of metabolic evolution are different, the neighboring genes are generally regulated by the same or related regulation factors. Therefore, the translation-based approaches associated with genomics, and transcriptomics are able to be employed for functional genomics focusing on plant secondary metabolism. Here, we present a survey of the current understanding of neighboring genes involved in plant secondary metabolism. Additionally, a genomic overview of neighboring genes of four model plants and transcriptional co-expression network neighboring genes to detect metabolic gene clusters in Arabidopsis is provided. Finally, the insights functional genomics have provided concerning the evolution and mechanistic regulation of both the formation and operation of metabolic neighboring clusters is discussed.