Complete biosynthesis of QS-21 in engineered yeast.

QS-21 is a potent vaccine adjuvant and remains the only saponin-based adjuvant that has been clinically approved for use in humans1,2. However, owing to the complex structure of QS-21, its availability is limited. Today, the supply depends on laborious extraction from the Chilean soapbark tree or on low-yielding total chemical synthesis3,4. Here we demonstrate the complete biosynthesis of QS-21 and its precursors, as well as structural derivatives, in engineered yeast strains. The successful biosynthesis in yeast requires fine-tuning of the host's native pathway fluxes, as well as the functional and balanced expression of 38 heterologous enzymes. The required biosynthetic pathway spans seven enzyme families-a terpene synthase, P450s, nucleotide sugar synthases, glycosyltransferases, a coenzyme A ligase, acyl transferases and polyketide synthases-from six organisms, and mimics in yeast the subcellular compartmentalization of plants from the endoplasmic reticulum membrane to the cytosol. Finally, by taking advantage of the promiscuity of certain pathway enzymes, we produced structural analogues of QS-21 using this biosynthetic platform. This microbial production scheme will allow for the future establishment of a structure-activity relationship, and will thus enable the rational design of potent vaccine adjuvants.

Metabolic Gene Clusters in Eukaryotes.

In bacteria, more than half of the genes in the genome are organized in operons. In contrast, in eukaryotes, functionally related genes are usually dispersed across the genome. There are, however, numerous examples of functional clusters of nonhomologous genes for metabolic pathways in fungi and plants. Despite superficial similarities with operons (physical clustering, coordinate regulation), these clusters have not usually originated by horizontal gene transfer from bacteria, and (unlike operons) the genes are typically transcribed separately rather than as a single polycistronic message. This clustering phenomenon raises intriguing questions about the origins of clustered metabolic pathways in eukaryotes and the significance of clustering for pathway function. Here we review metabolic gene clusters from fungi and plants, highlight commonalities and differences, and consider how these clusters form and are regulated. We also identify opportunities for future research in the areas of large-scale genomics, synthetic biology, and experimental evolution.

Complete biosynthesis of the potent vaccine adjuvant QS-21.

QS-21 is a potent vaccine adjuvant currently sourced by extraction from the Chilean soapbark tree. It is a key component of human vaccines for shingles, malaria, coronavirus disease 2019 and others under development. The structure of QS-21 consists of a glycosylated triterpene scaffold coupled to a complex glycosylated 18-carbon acyl chain that is critical for immunostimulant activity. We previously identified the early pathway steps needed to make the triterpene glycoside scaffold; however, the biosynthetic route to the acyl chain, which is needed for stimulation of T cell proliferation, was unknown. Here, we report the biogenic origin of the acyl chain, characterize the series of enzymes required for its synthesis and addition and reconstitute the entire 20-step pathway in tobacco, thereby demonstrating the production of QS-21 in a heterologous expression system. This advance opens up unprecedented opportunities for bioengineering of vaccine adjuvants, investigating structure-activity relationships and understanding the mechanisms by which these compounds promote the human immune response.

Identification and functional validation of super-enhancers in Arabidopsis thaliana.

Super-enhancers (SEs) are exceptionally large enhancers and are recognized to play prominent roles in cell identity in mammalian species. We surveyed the genomic regions containing large clusters of accessible chromatin regions (ACRs) marked by deoxyribonuclease (DNase) I hypersensitivity in Arabidopsis thaliana. We identified a set of 749 putative SEs, which have a minimum length of 1.5 kilobases and represent the top 2.5% of the largest ACR clusters. We demonstrate that the genomic regions associating with these SEs were more sensitive to DNase I than other nonpromoter ACRs. The SEs were preferentially associated with topologically associating domains. Furthermore, the SEs and their predicted cognate genes were frequently associated with organ development and tissue identity in A. thaliana. Therefore, the A. thaliana SEs and their cognate genes mirror the functional characteristics of those reported in mammalian species. We developed CRISPR/Cas-mediated deletion lines of a 3,578-bp SE associated with the thalianol biosynthetic gene cluster (BGC). Small deletions (131-157 bp) within the SE resulted in distinct phenotypic changes and transcriptional repression of all five thalianol genes. In addition, T-DNA insertions in the SE region resulted in transcriptional alteration of all five thalianol genes. Thus, this SE appears to play a central role in coordinating the operon-like expression pattern of the thalianol BGC.

Pathogen-induced biosynthetic pathways encode defense-related molecules in bread wheat.

Wheat is a widely grown food crop that suffers major yield losses due to attack by pests and pathogens. A better understanding of biotic stress responses in wheat is thus of major importance. The recently assembled bread wheat genome coupled with extensive transcriptomic resources provides unprecedented new opportunities to investigate responses to pathogen challenge. Here, we analyze gene coexpression networks to identify modules showing consistent induction in response to pathogen exposure. Within the top pathogen-induced modules, we identify multiple clusters of physically adjacent genes that correspond to six pathogen-induced biosynthetic pathways that share a common regulatory network. Functional analysis reveals that these pathways, all of which are encoded by biosynthetic gene clusters, produce various different classes of compounds­namely, flavonoids, diterpenes, and triterpenes, including the defense-related compound ellarinacin. Through comparative genomics, we also identify associations with the known rice phytoalexins momilactones, as well as with a defense-related gene cluster in the grass model plant Brachypodium distachyon. Our results significantly advance the understanding of chemical defenses in wheat and open up avenues for enhancing disease resistance in this agriculturally important crop. They also exemplify the power of transcriptional networks to discover the biosynthesis of chemical defenses in plants with large, complex genomes.

Biosynthesis of saponin defensive compounds in sea cucumbers.

Soft-bodied slow-moving sea creatures such as sea stars and sea cucumbers lack an adaptive immune system and have instead evolved the ability to make specialized protective chemicals (glycosylated steroids and triterpenes) as part of their innate immune system. This raises the intriguing question of how these biosynthetic pathways have evolved. Sea star saponins are steroidal, while those of the sea cucumber are triterpenoid. Sterol biosynthesis in animals involves cyclization of 2,3-oxidosqualene to lanosterol by the oxidosqualene cyclase (OSC) enzyme lanosterol synthase (LSS). Here we show that sea cucumbers lack LSS and instead have two divergent OSCs that produce triterpene saponins and that are likely to have evolved from an ancestral LSS by gene duplication and neofunctionalization. We further show that sea cucumbers make alternate sterols that confer protection against self-poisoning by their own saponins. Collectively, these events have enabled sea cucumbers to evolve the ability to produce saponins and saponin-resistant sterols concomitantly.

Vacuolar MATE/DTX protein-mediated cucurbitacin C transport is co-regulated with bitterness biosynthesis in cucumber.

Membrane-localized transporters constitute important components for specialized metabolism in plants. However, due to the vast array of specialized metabolites produced by plants, and the large families of transporter genes, knowledge about the intracellular and intercellular transport of plant metabolites is still in its infancy. Cucurbitacins are bitter and defensive triterpenoids produced mainly in the cucurbits. Using a comparative genomics and multi-omics approach, a MATE gene (CsMATE1), physically clustered with cucurbitacin C (CuC) biosynthetic genes, was identified and functionally shown to sequester CuC in cucumber leaf mesophyll cells. Notably, the CuC transport process is strictly co-regulated with CuC biosynthesis. CsMATE1 clustering with bitterness biosynthesis genes may provide benefits and a basis for this feedback regulation on CuC sequestration and biosynthesis. Identification of transport systems for plant-specialized metabolites can accelerate the metabolic engineering of high-value-added compounds by simplifying their purification process.

OeBAS and CYP716C67 catalyze the biosynthesis of health-beneficial triterpenoids in olive (Olea europaea) fruits.

The bioactive properties of olive (Olea europaea) fruits and olive oil are largely attributed to terpenoid compounds, including diverse triterpenoids such as oleanolic, maslinic and ursolic acids, erythrodiol, and uvaol. They have applications in the agri-food, cosmetics, and pharmaceutical industries. Some key steps involved in the biosynthesis of these compounds are still unknown. Genome mining, biochemical analysis, and trait association studies have been used to identify major gene candidates controlling triterpenoid content of olive fruits. Here, we identify and functionally characterize an oxidosqualene cyclase (OeBAS) required for the production of the major triterpene scaffold ß-amyrin, the precursor of erythrodiol, oleanolic and maslinic acids, and a cytochrome P450 (CYP716C67) that mediates 2α oxidation of the oleanane- and ursane-type triterpene scaffolds to produce maslinic and corosolic acids, respectively. To confirm the enzymatic functions of the entire pathway, we have reconstituted the olive biosynthetic pathway for oleanane- and ursane-type triterpenoids in the heterologous host, Nicotiana benthamiana. Finally, we have identified genetic markers associated with oleanolic and maslinic acid fruit content on the chromosomes carrying the OeBAS and CYP716C67 genes. Our results shed light on the biosynthesis of olive triterpenoids and provide new gene targets for germplasm screening and breeding for high triterpenoid content.

Active and repressed biosynthetic gene clusters have spatially distinct chromosome states.

While colocalization within a bacterial operon enables coexpression of the constituent genes, the mechanistic logic of clustering of nonhomologous monocistronic genes in eukaryotes is not immediately obvious. Biosynthetic gene clusters that encode pathways for specialized metabolites are an exception to the classical eukaryote rule of random gene location and provide paradigmatic exemplars with which to understand eukaryotic cluster dynamics and regulation. Here, using 3C, Hi-C, and Capture Hi-C (CHi-C) organ-specific chromosome conformation capture techniques along with high-resolution microscopy, we investigate how chromosome topology relates to transcriptional activity of clustered biosynthetic pathway genes in Arabidopsis thaliana Our analyses reveal that biosynthetic gene clusters are embedded in local hot spots of 3D contacts that segregate cluster regions from the surrounding chromosome environment. The spatial conformation of these cluster-associated domains differs between transcriptionally active and silenced clusters. We further show that silenced clusters associate with heterochromatic chromosomal domains toward the periphery of the nucleus, while transcriptionally active clusters relocate away from the nuclear periphery. Examination of chromosome structure at unrelated clusters in maize, rice, and tomato indicates that integration of clustered pathway genes into distinct topological domains is a common feature in plant genomes. Our results shed light on the potential mechanisms that constrain coexpression within clusters of nonhomologous eukaryotic genes and suggest that gene clustering in the one-dimensional chromosome is accompanied by compartmentalization of the 3D chromosome.

Identification of key enzymes responsible for protolimonoid biosynthesis in plants: Opening the door to azadirachtin production.

Limonoids are natural products made by plants belonging to the Meliaceae (Mahogany) and Rutaceae (Citrus) families. They are well known for their insecticidal activity, contribution to bitterness in citrus fruits, and potential pharmaceutical properties. The best known limonoid insecticide is azadirachtin, produced by the neem tree (Azadirachta indica). Despite intensive investigation of limonoids over the last half century, the route of limonoid biosynthesis remains unknown. Limonoids are classified as tetranortriterpenes because the prototypical 26-carbon limonoid scaffold is postulated to be formed from a 30-carbon triterpene scaffold by loss of 4 carbons with associated furan ring formation, by an as yet unknown mechanism. Here we have mined genome and transcriptome sequence resources for 3 diverse limonoid-producing species (A. indica, Melia azedarach, and Citrus sinensis) to elucidate the early steps in limonoid biosynthesis. We identify an oxidosqualene cyclase able to produce the potential 30-carbon triterpene scaffold precursor tirucalla-7,24-dien-3ß-ol from each of the 3 species. We further identify coexpressed cytochrome P450 enzymes from M. azedarach (MaCYP71CD2 and MaCYP71BQ5) and C. sinensis (CsCYP71CD1 and CsCYP71BQ4) that are capable of 3 oxidations of tirucalla-7,24-dien-3ß-ol, resulting in spontaneous hemiacetal ring formation and the production of the protolimonoid melianol. Our work reports the characterization of protolimonoid biosynthetic enzymes from different plant species and supports the notion of pathway conservation between both plant families. It further paves the way for engineering crop plants with enhanced insect resistance and producing high-value limonoids for pharmaceutical and other applications by expression in heterologous hosts.

A noncanonical vacuolar sugar transferase required for biosynthesis of antimicrobial defense compounds in oat.

Plants produce an array of natural products with important ecological functions. These compounds are often decorated with oligosaccharide groups that influence bioactivity, but the biosynthesis of such sugar chains is not well understood. Triterpene glycosides (saponins) are a large family of plant natural products that determine important agronomic traits, as exemplified by avenacins, antimicrobial defense compounds produced by oats. Avenacins have a branched trisaccharide moiety consisting of l-arabinose linked to 2 d-glucose molecules that is critical for antifungal activity. Plant natural product glycosylation is usually performed by uridine diphosphate-dependent glycosyltransferases (UGTs). We previously characterized the arabinosyltransferase that initiates the avenacin sugar chain; however, the enzymes that add the 2 remaining d-glucose molecules have remained elusive. Here we characterize the enzymes that catalyze these last 2 glucosylation steps. AsUGT91G16 is a classical cytosolic UGT that adds a 1,2-linked d-glucose molecule to l-arabinose. Unexpectedly, the enzyme that adds the final 1,4-linked d-glucose (AsTG1) is not a UGT, but rather a sugar transferase belonging to Glycosyl Hydrolase family 1 (GH1). Unlike classical UGTs, AsTG1 is vacuolar. Analysis of oat mutants reveals that AsTG1 corresponds to Sad3, a previously uncharacterized locus shown by mutation to be required for avenacin biosynthesis. AsTG1 and AsUGT91G16 form part of the avenacin biosynthetic gene cluster. Our demonstration that a vacuolar transglucosidase family member plays a critical role in triterpene biosynthesis highlights the importance of considering other classes of carbohydrate-active enzymes in addition to UGTs as candidates when elucidating pathways for the biosynthesis of glycosylated natural products in plants.

Modulation of Arabidopsis root growth by specialized triterpenes.

Plant roots are specialized belowground organs that spatiotemporally shape their development in function of varying soil conditions. This root plasticity relies on intricate molecular networks driven by phytohormones, such as auxin and jasmonate (JA). Loss-of-function of the NOVEL INTERACTOR OF JAZ (NINJA), a core component of the JA signaling pathway, leads to enhanced triterpene biosynthesis, in particular of the thalianol gene cluster, in Arabidopsis thaliana roots. We have investigated the biological role of thalianol and its derivatives by focusing on Thalianol Synthase (THAS) and Thalianol Acyltransferase 2 (THAA2), two thalianol cluster genes that are upregulated in the roots of ninja mutant plants. THAS and THAA2 activity was investigated in yeast, and metabolite and phenotype profiling of thas and thaa2 loss-of-function plants was carried out. THAA2 was shown to be responsible for the acetylation of thalianol and its derivatives, both in yeast and in planta. In addition, THAS and THAA2 activity was shown to modulate root development. Our results indicate that the thalianol pathway is not only controlled by phytohormonal cues, but also may modulate phytohormonal action itself, thereby affecting root development and interaction with the environment.

Analysis of Two New Arabinosyltransferases Belonging to the Carbohydrate-Active Enzyme (CAZY) Glycosyl Transferase Family1 Provides Insights into Disease Resistance and Sugar Donor Specificity.

Glycosylation of small molecules is critical for numerous biological processes in plants, including hormone homeostasis, neutralization of xenobiotics, and synthesis and storage of specialized metabolites. Glycosylation of plant natural products is usually performed by uridine diphosphate-dependent glycosyltransferases (UGTs). Triterpene glycosides (saponins) are a large family of plant natural products that determine important agronomic traits such as disease resistance and flavor and have numerous pharmaceutical applications. Most characterized plant natural product UGTs are glucosyltransferases, and little is known about enzymes that add other sugars. Here we report the discovery and characterization of AsAAT1 (UGT99D1), which is required for biosynthesis of the antifungal saponin avenacin A-1 in oat (Avena strigosa). This enzyme adds l-Ara to the triterpene scaffold at the C-3 position, a modification critical for disease resistance. The only previously reported plant natural product arabinosyltransferase is a flavonoid arabinosyltransferase from Arabidopsis (Arabidopsis thaliana). We show that AsAAT1 has high specificity for UDP-ß-l-arabinopyranose, identify two amino acids required for sugar donor specificity, and through targeted mutagenesis convert AsAAT1 into a glucosyltransferase. We further identify a second arabinosyltransferase potentially implicated in the biosynthesis of saponins that determine bitterness in soybean (Glycine max). Our investigations suggest independent evolution of UDP-Ara sugar donor specificity in arabinosyltransferases in monocots and eudicots.

Multigenome analysis implicates miniature inverted-repeat transposable elements (MITEs) in metabolic diversification in eudicots.

Plants produce a plethora of natural products, including many drugs. It has recently emerged that the genes encoding different natural product pathways may be organized as biosynthetic gene clusters in plant genomes, with >30 examples reported so far. Despite superficial similarities with microbes, these clusters have not arisen by horizontal gene transfer, but rather by gene duplication, neofunctionalization, and relocation via unknown mechanisms. Previously we reported that two Arabidopsis thaliana biosynthetic gene clusters are located in regions of the genome that are significantly enriched in transposable elements (TEs). Other plant biosynthetic gene clusters also harbor abundant TEs. TEs can mediate genomic rearrangement by providing homologous sequences that enable illegitimate recombination and gene relocation. Thus, TE-mediated recombination may contribute to plant biosynthetic gene cluster formation. TEs may also facilitate establishment of regulons. However, a systematic analysis of the TEs associated with plant biosynthetic gene clusters has not been carried out. Here we investigate the TEs associated with clustered terpene biosynthetic genes in multiple plant genomes and find evidence to suggest a role for miniature inverted-repeat transposable elements in cluster formation in eudicots. Through investigation of the newly sequenced Amborella trichopoda, Aquilegia coerulea, and Kalanchoe fedtschenkoi genomes, we further show that the "block" mechanism of founding of biosynthetic gene clusters through duplication and diversification of pairs of terpene synthase and cytochrome P450 genes that is prevalent in the eudicots arose around 90-130 million years ago, after the appearance of the basal eudicots and before the emergence of the superrosid clade.

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

Plants produce an array of specialized metabolites with important ecological functions. The mechanisms underpinning the evolution of new biosynthetic pathways are not well-understood. Here, we exploit available genome sequence resources to investigate triterpene biosynthesis across the Brassicaceae. Oxidosqualene cyclases (OSCs) catalyze the first committed step in triterpene biosynthesis. Systematic analysis of 13 sequenced Brassicaceae genomes was performed to identify all OSC genes. The genome neighbourhoods (GNs) around a total of 163 OSC genes were investigated to identify Pfam domains significantly enriched in these regions. All-vs-all comparisons of OSC neighbourhoods and phylogenomic analysis were used to investigate the sequence similarity and evolutionary relationships of the numerous candidate triterpene biosynthetic gene clusters (BGCs) observed. Functional analysis of three representative BGCs was carried out and their triterpene pathway products were elucidated. Our results indicate that plant genomes are remarkably plastic, and that dynamic GNs generate new biosynthetic pathways in different Brassicaceae lineages by shuffling the genes encoding a core palette of triterpene-diversifying enzymes, presumably in response to strong environmental selection pressure. These results illuminate a genomic basis for diversification of plant-specialized metabolism through natural combinatorics of enzyme families, which can be mimicked using synthetic biology to engineer diverse bioactive molecules.

Unearthing a sesterterpene biosynthetic repertoire in the Brassicaceae through genome mining reveals convergent evolution.

Sesterterpenoids are a rare terpene class harboring untapped chemodiversity and bioactivities. Their structural diversity originates primarily from the scaffold-generating sesterterpene synthases (STSs). In fungi, all six known STSs are bifunctional, containing C-terminal trans-prenyltransferase (PT) and N-terminal terpene synthase (TPS) domains. In plants, two colocalized PT and TPS gene pairs from Arabidopsis thaliana were recently reported to synthesize sesterterpenes. However, the landscape of PT and TPS genes in plant genomes is unclear. Here, using a customized algorithm for systematically searching plant genomes, we reveal a suite of physically colocalized pairs of PT and TPS genes for the biosynthesis of a large sesterterpene repertoire in the wider Brassicaceae. Transient expression of seven TPSs from A. thaliana, Capsella rubella, and Brassica oleracea in Nicotiana benthamiana yielded fungal-type sesterterpenes with tri-, tetra-, and pentacyclic scaffolds, and notably (-)-ent-quiannulatene, an enantiomer of the fungal metabolite (+)-quiannulatene. Protein and structural modeling analysis identified an amino acid site implicated in structural diversification. Mutation of this site in one STS (AtTPS19) resulted in premature termination of carbocation intermediates and accumulation of bi-, tri-, and tetracyclic sesterterpenes, revealing the cyclization path for the pentacyclic sesterterpene (-)-retigeranin B. These structural and mechanistic insights, together with phylogenetic analysis, suggest convergent evolution of plant and fungal STSs, and also indicate that the colocalized PT-TPS gene pairs in the Brassicaceae may have originated from a common ancestral gene pair present before speciation. Our findings further provide opportunities for rapid discovery and production of sesterterpenes through metabolic and protein engineering.

The protosteryl and dammarenyl cation dichotomy in polycyclic triterpene biosynthesis revisited: has this 'rule' finally been broken?

Covering: 1948 up to the end of 2018. The triterpene alcohols represent an important and diverse class of natural products. This diversity is believed to originate from the differential enzymatically controlled cyclisation of 2,3-oxidosqualene. It is now a well-established presumption that all naturally occurring tetra- and penta-cyclic triterpene alcohols can be rationalised by the resolution of one of two intermediary tetracyclic cations, termed the protosteryl and dammarenyl cations. Here, a discussion of typical key triterpene structures and their proposed derivation from either of these progenitors is followed by comparison with a recently reported novel pentacyclic triterpene orysatinol which appears to correspond to an unprecedented divergence from this dichotomous protosteryl/dammarenyl view of triterpene biogenesis. Not only does this discovery widen the potential scope of triterpene scaffolds that could exist in nature, it could call into question the reliability of stereochemical assignments of some existing triterpene structures that are supported by only limited spectroscopic evidence. The discovery of orysatinol provides direct experimental evidence to support considering more flexibility in the stereochemical interpretation of the biogenic isoprene rule.

Light-Induced Artemisinin Biosynthesis Is Regulated by the bZIP Transcription Factor AaHY5 in Artemisia annua.

Artemisinin, the frontline drug against malaria, is a sesquiterpenoid extracted from Artemisia annua. Light has been proposed to play an important role in the activation of artemisinin biosynthesis. Here, we report the basic leucine zipper transcription factor (TF) AaHY5 as a key regulator of light-induced biosynthesis of artemisinin. We show that AaHY5 transcription overlaps with that of artemisinin biosynthesis genes in response to light and in A. annua tissues. Analysis of AaHY5 overexpression and RNAi-suppression lines suggests that AaHY5 is a positive regulator of the expression of artemisinin biosynthesis genes and accumulation of artemisinin. We show that AaHY5 complements the hy5 mutant in Arabidopsis thaliana. Our data further suggest that AaHY5 interacts with AaCOP1, the ubiquitin E3 ligase CONSTITUTIVE PHOTOMORPHOGENIC1 in A. annua. In yeast one-hybrid and transient expression assays, we demonstrate that AaHY5 acts via the TF GLANDULAR TRICHOME-SPECIFIC WRKY 1 (AaGSW1) in artemisinin regulation. In summary, we present a novel regulator of artemisinin gene expression and propose a model in which AaHY5 indirectly controls artemisinin production in response to changing light conditions.

Towards take-all control: a C-21ß oxidase required for acylation of triterpene defence compounds in oat.

Oats produce avenacins, antifungal triterpenes that are synthesized in the roots and provide protection against take-all and other soilborne diseases. Avenacins are acylated at the carbon-21 position of the triterpene scaffold, a modification critical for antifungal activity. We have previously characterized several steps in the avenacin pathway, including those required for acylation. However, transfer of the acyl group to the scaffold requires the C-21ß position to be oxidized first, by an as yet uncharacterized enzyme. We mined oat transcriptome data to identify candidate cytochrome P450 enzymes that may catalyse C-21ß oxidation. Candidates were screened for activity by transient expression in Nicotiana benthamiana. We identified a cytochrome P450 enzyme AsCYP72A475 as a triterpene C-21ß hydroxylase, and showed that expression of this enzyme together with early pathway steps yields C-21ß oxidized avenacin intermediates. We further demonstrate that AsCYP72A475 is synonymous with Sad6, a previously uncharacterized locus required for avenacin biosynthesis. sad6 mutants are compromised in avenacin acylation and have enhanced disease susceptibility. The discovery of AsCYP72A475 represents an important advance in the understanding of triterpene biosynthesis and paves the way for engineering the avenacin pathway into wheat and other cereals for control of take-all and other diseases.

plantiSMASH: automated identification, annotation and expression analysis of plant biosynthetic gene clusters.

Plant specialized metabolites are chemically highly diverse, play key roles in host-microbe interactions, have important nutritional value in crops and are frequently applied as medicines. It has recently become clear that plant biosynthetic pathway-encoding genes are sometimes densely clustered in specific genomic loci: biosynthetic gene clusters (BGCs). Here, we introduce plantiSMASH, a versatile online analysis platform that automates the identification of candidate plant BGCs. Moreover, it allows integration of transcriptomic data to prioritize candidate BGCs based on the coexpression patterns of predicted biosynthetic enzyme-coding genes, and facilitates comparative genomic analysis to study the evolutionary conservation of each cluster. Applied on 48 high-quality plant genomes, plantiSMASH identifies a rich diversity of candidate plant BGCs. These results will guide further experimental exploration of the nature and dynamics of gene clustering in plant metabolism. Moreover, spurred by the continuing decrease in costs of plant genome sequencing, they will allow genome mining technologies to be applied to plant natural product discovery. The plantiSMASH web server, precalculated results and source code are freely available from http://plantismash.secondarymetabolites.org.