The Rauvolfia tetraphylla genome suggests multiple distinct biosynthetic routes for yohimbane monoterpene indole alkaloids.

Monoterpene indole alkaloids (MIAs) are a structurally diverse family of specialized metabolites mainly produced in Gentianales to cope with environmental challenges. Due to their pharmacological properties, the biosynthetic modalities of several MIA types have been elucidated but not that of the yohimbanes. Here, we combine metabolomics, proteomics, transcriptomics and genome sequencing of Rauvolfia tetraphylla with machine learning to discover the unexpected multiple actors of this natural product synthesis. We identify a medium chain dehydrogenase/reductase (MDR) that produces a mixture of four diastereomers of yohimbanes including the well-known yohimbine and rauwolscine. In addition to this multifunctional yohimbane synthase (YOS), an MDR synthesizing mainly heteroyohimbanes and the short chain dehydrogenase vitrosamine synthase also display a yohimbane synthase side activity. Lastly, we establish that the combination of geissoschizine synthase with at least three other MDRs also produces a yohimbane mixture thus shedding light on the complex mechanisms evolved for the synthesis of these plant bioactives.

Identification of a second 16-hydroxytabersonine-O-methyltransferase suggests an evolutionary relationship between alkaloid and flavonoid metabolisms in Catharanthus roseus.

The medicinal plant Catharanthus roseus biosynthesizes many important drugs for human health, including the anticancer monoterpene indole alkaloids (MIAs) vinblastine and vincristine. Over the past decades, the continuous increase in pharmaceutical demand has prompted several research groups to characterize MIA biosynthetic pathways for considering future metabolic engineering processes of supply. In line with previous work suggesting that diversification can potentially occur at various steps along the vindoline branch, we were here interested in investigating the involvement of distinct isoforms of tabersonine-16-O-methyltransferase (16OMT) which plays a pivotal role in the MIA biosynthetic pathway. By combining homology searches based on the previously characterized 16OMT1, phylogenetic analyses, functional assays in yeast, and biochemical and in planta characterizations, we identified a second isoform of 16OMT, referred to as 16OMT2. 16OMT2 appears to be a multifunctional enzyme working on both MIA and flavonoid substrates, suggesting that a constrained evolution of the enzyme for accommodating the MIA substrate has probably occurred to favor the apparition of 16OMT2 from an ancestral specific flavonoid-O-methyltransferase. Since 16OMT1 and 16OMT2 displays a high sequence identity and similar kinetic parameters for 16-hydroxytabersonine, we postulate that 16OMT1 may result from a later 16OMT2 gene duplication accompanied by a continuous neofunctionalization leading to an almost complete loss of flavonoid O-methyltransferase activity. Overall, these results participate in increasing our knowledge on the evolutionary processes that have likely led to enzyme co-optation for MIA synthesis.

The Vinca minor genome highlights conserved evolutionary traits in monoterpene indole alkaloid synthesis.

Vinca minor, also known as the lesser periwinkle, is a well-known species from the Apocynaceae, native to central and southern Europe. This plant synthesizes monoterpene indole alkaloids, which are a class of specialized metabolites displaying a wide range of bioactive- and pharmacologically important properties. Within the almost 50 monoterpene indole alkaloids it produces, V. minor mainly accumulates vincamine, which is commercially used as a nootropic. Using a combination of Oxford Nanopore Technologies long read- and Illumina short-read sequencing, a 679,098 Mb V. minor genome was assembled into 296 scaffolds with an N50 scaffold length of 6 Mb, and encoding 29,624 genes. These genes were functionally annotated and used in a comparative genomic analysis to establish gene families and to investigate gene family expansion and contraction across the phylogenetic tree. Furthermore, homology-based monoterpene indole alkaloid gene predictions together with a metabolic analysis across 4 different V. minor tissue types guided the identification of candidate monoterpene indole alkaloid genes. These candidates were finally used to identify monoterpene indole alkaloid gene clusters, which combined with synteny analysis allowed for the discovery of a functionally validated vincadifformine-16-hydroxylase, reinforcing the potential of this dataset for monoterpene indole alkaloids gene discovery. It is expected that access to these resources will facilitate the elucidation of unknown monoterpene indole alkaloid biosynthetic routes with the potential of transferring these pathways to heterologous expression systems for large-scale monoterpene indole alkaloid production.

Genome Assembly of the Medicinal Plant Voacanga thouarsii.

The Apocynaceae tree Voacanga thouarsii, native to southern Africa and Madagascar, produces monoterpene indole alkaloids (MIA), which are specialized metabolites with a wide range of bioactive properties. Voacanga species mainly accumulates tabersonine in seeds making these species valuable medicinal plants currently used for industrial MIA production. Despite their importance, the MIA biosynthesis in Voacanga species remains poorly studied. Here, we report the first genome assembly and annotation of a Voacanga species. The combined assembly of Oxford Nanopore Technologies long-reads and Illumina short-reads resulted in 3,406 scaffolds with a total length of 1,354.26 Mb and an N50 of 3.04 Mb. A total of 33,300 protein-coding genes were predicted and functionally annotated. These genes were then used to establish gene families and to investigate gene family expansion and contraction across the phylogenetic tree. A transposable element (TE) analysis showed the highest proportion of TE in Voacanga thouarsii compared with all other MIA-producing plants. In a nutshell, this first reference genome of V. thouarsii will thus contribute to strengthen future comparative and evolutionary studies in MIA-producing plants leading to a better understanding of MIA pathway evolution. This will also allow the potential identification of new MIA biosynthetic genes for metabolic engineering purposes.

A microbial supply chain for production of the anti-cancer drug vinblastine.

Monoterpene indole alkaloids (MIAs) are a diverse family of complex plant secondary metabolites with many medicinal properties, including the essential anti-cancer therapeutics vinblastine and vincristine1. As MIAs are difficult to chemically synthesize, the world's supply chain for vinblastine relies on low-yielding extraction and purification of the precursors vindoline and catharanthine from the plant Catharanthus roseus, which is then followed by simple in vitro chemical coupling and reduction to form vinblastine at an industrial scale2,3. Here, we demonstrate the de novo microbial biosynthesis of vindoline and catharanthine using a highly engineered yeast, and in vitro chemical coupling to vinblastine. The study showcases a very long biosynthetic pathway refactored into a microbial cell factory, including 30 enzymatic steps beyond the yeast native metabolites geranyl pyrophosphate and tryptophan to catharanthine and vindoline. In total, 56 genetic edits were performed, including expression of 34 heterologous genes from plants, as well as deletions, knock-downs and overexpression of ten yeast genes to improve precursor supplies towards de novo production of catharanthine and vindoline, from which semisynthesis to vinblastine occurs. As the vinblastine pathway is one of the longest MIA biosynthetic pathways, this study positions yeast as a scalable platform to produce more than 3,000 natural MIAs and a virtually infinite number of new-to-nature analogues.

RNA-seq Analysis of Monoterpene Indole Alkaloid Biosynthetic Pathway Elucidation in Catharanthus roseus.

Increased affordability and availability of high-throughput next-generation sequencing (NGS) technologies have resulted in an explosion of available RNA-seq data, igniting a variety of data-mining methodologies, valuable for plant-specialized biosynthetic pathway discovery. When combined with traditional homology-based annotations, these methods can facilitate short-listing candidate genes for downstream functional validation screenings. Genes related to common pathways often display homogenous expression patterns across different tissue types and experimental conditions. Here, we describe bioinformatic protocols for exploiting such coexpression to shortlist candidate genes of the well-described monoterpene indole alkaloid (MIA) pathway of Catharanthus roseus. These methods aim to inspire researchers to utilize this publicly available RNA-seq treasure trove to guide their own endeavors in the characterization of missing steps in plant metabolic pathways.

Predicting Monoterpene Indole Alkaloid-Related Genes from Expression Data with Artificial Neural Networks.

Elucidation of biological pathways leading to specialized metabolites remains a complex task. It is however a mandatory step to allow bioproduction into heterologous hosts. Many steps have already been identified using conventional approaches, enlarging the space of known possible chemical steps. In the recent past years, identification of missing steps has been fueled by the generation of genomic and transcriptomic data for nonmodel species. The analysis of gene expression profiles has revealed that in many cases, genes encoding enzymes involved in the same biosynthetic pathways are coexpressed across different tissue types and environmental conditions. Hence, coexpressed studies, either in the form of differential gene expression, gene coexpression network, or unsupervised clustering methods, have helped deciphering missing steps to complete knowledge on biosynthetic pathways. Already identified biosynthetic steps can be used as baits to capture the remaining unknown steps. The present protocol shows how supervised machine learning in the form of artificial neural networks (ANNs) can efficiently classify genes as specialized metabolism related or not according to their expression levels. Using Catharanthus roseus as an example, we show that ANN trained on a minimal set of bait genes results in many true positives (correctly predicted genes) while keeping false positives low (containing possible candidate genes).

Tonoplast and Peroxisome Targeting of γ-tocopherol N-methyltransferase Homologs Involved in the Synthesis of Monoterpene Indole Alkaloids.

Many plant species from the Apocynaceae, Loganiaceae and Rubiaceae families evolved a specialized metabolism leading to the synthesis of a broad palette of monoterpene indole alkaloids (MIAs). These compounds are believed to constitute a cornerstone of the plant chemical arsenal but above all several MIAs display pharmacological properties that have been exploited for decades by humans to treat various diseases. It is established that MIAs are produced in planta due to complex biosynthetic pathways engaging a multitude of specialized enzymes but also a complex tissue and subcellular organization. In this context, N-methyltransferases (NMTs) represent an important family of enzymes indispensable for MIA biosynthesis but their characterization has always remained challenging. In particular, little is known about the subcellular localization of NMTs in MIA-producing plants. Here, we performed an extensive analysis on the subcellular localization of NMTs from four distinct medicinal plants but also experimentally validated that two putative NMTs from Catharanthus roseus exhibit NMT activity. Apart from providing unprecedented data regarding the targeting of these enzymes in planta, our results point out an additional layer of complexity to the subcellular organization of the MIA biosynthetic pathway by introducing tonoplast and peroxisome as new actors of the final steps of MIA biosynthesis.

Chromosome-scale genomes throw light on plant drug biosynthesis.

Recent chromosome-scale genomes from prominent medicinal plants have provided unprecedented insight into the architecture and evolution of some prominent metabolic pathways. These new technologies also facilitate the identification of plant drug biosynthetic genes and would likely accelerate the development of new bioengineering procedures to secure the supply of important plant-derived pharmaceuticals.

An updated version of the Madagascar periwinkle genome.

The Madagascar periwinkle, Catharanthus roseus, belongs to the Apocynaceae family. This medicinal plant, endemic to Madagascar, produces many important drugs including the monoterpene indole alkaloids (MIA) vincristine and vinblastine used to treat cancer worldwide. Here, we provide a new version of the C. roseus genome sequence obtained through the combination of Oxford Nanopore Technologies long-reads and Illumina short-reads. This more contiguous assembly consists of 173 scaffolds with a total length of 581.128 Mb and an N50 of 12.241 Mb. Using publicly available RNAseq data, 21,061 protein coding genes were predicted and functionally annotated. A total of 42.87% of the genome was annotated as transposable elements, most of them being long-terminal repeats. Together with the increasing access to MIA-producing plant genomes, this updated version should ease evolutionary studies leading to a better understanding of MIA biosynthetic pathway evolution.

Improved virus-induced gene silencing allows discovery of a serpentine synthase gene in Catharanthus roseus.

Specialized metabolites are chemically complex small molecules with a myriad of biological functions. To investigate plant-specialized metabolite biosynthesis more effectively, we developed an improved method for virus-induced gene silencing (VIGS). We designed a plasmid that incorporates fragments of both the target gene and knockdown marker gene (phytoene desaturase, PDS), which identifies tissues that have been successfully silenced in planta. To demonstrate the utility of this method, we used the terpenoid indole alkaloid (TIA) pathway in Madagascar periwinkle (Catharanthus roseus) as a model system. Catharanthus roseus is a medicinal plant well known for producing many bioactive compounds, such as vinblastine and vincristine. Our VIGS method enabled the discovery of a previously unknown biosynthetic enzyme, serpentine synthase (SS). This enzyme is a cytochrome P450 (CYP) that produces the ß-carboline alkaloids serpentine and alstonine, compounds with strong blue autofluorescence and potential pharmacological activity. The discovery of this enzyme highlights the complexity of TIA biosynthesis and demonstrates the utility of this improved VIGS method for discovering unidentified metabolic enzymes in plants.

Computational biotechnology guides elucidation of the biosynthesis of the plant anticancer drug camptothecin.

Camptothecin is a clinically important monoterpene indole alkaloid (MIAs) used for treating various cancers. Currently, the production of this biopharmaceutical hinges on its extraction from camptothecin-producing plants, leading to high market prices and supply bottlenecks. While synthetic biology combined with metabolic approaches could represent an attractive alternative approach to manufacturing, it requires firstly a complete biosynthetic pathway elucidation, which is, unfortunately, severely missing in species naturally accumulating camptothecin. This knowledge gap can be attributed to the lack of high-quality genomic resources of medicinal plant species. In such a perspective, Yamazaki and colleagues produced the first described and experimentally validated chromosome-level plant genome assembly of Ophiorrhiza pumila, a prominent source plant of camptothecin for the pharmaceutical industry. More specifically, they have developed a method incorporating Illumina reads, PacBio single-molecule reads, optical mapping and Hi-C sequencing, followed by the experimental validation of contig orientation within scaffolds, using fluorescence in situ hybridization (FISH) analysis. This relevant strategy resulted in the most contiguous and complete de novo plant reference genome described to date, which can streamline the sequencing of new plant genomes. Further mining approaches, including integrative omics analysis, phylogenetics, gene cluster evaluation and comparative genomics were successfully used to puzzle out the evolutionary origins of MIA metabolism and revealed a short-list of high confidence MIA biosynthetic genes for functional validation.

Alternative splicing creates a pseudo-strictosidine ß-d-glucosidase modulating alkaloid synthesis in Catharanthus roseus.

Deglycosylation is a key step in the activation of specialized metabolites involved in plant defense mechanisms. This reaction is notably catalyzed by ß-glucosidases of the glycosyl hydrolase 1 (GH1) family such as strictosidine ß-d-glucosidase (SGD) from Catharanthus roseus. SGD catalyzes the deglycosylation of strictosidine, forming a highly reactive aglycone involved in the synthesis of cytotoxic monoterpene indole alkaloids (MIAs) and in the crosslinking of aggressor proteins. By exploring C. roseus transcriptomic resources, we identified an alternative splicing event of the SGD gene leading to the formation of a shorter isoform of this enzyme (shSGD) that lacks the last 71-residues and whose transcript ratio with SGD ranges from 1.7% up to 42.8%, depending on organs and conditions. Whereas it completely lacks ß-glucosidase activity, shSGD interacts with SGD and causes the disruption of SGD multimers. Such disorganization drastically inhibits SGD activity and impacts downstream MIA synthesis. In addition, shSGD disrupts the metabolic channeling of downstream biosynthetic steps by hampering the recruitment of tetrahydroalstonine synthase in cell nuclei. shSGD thus corresponds to a pseudo-enzyme acting as a regulator of MIA biosynthesis. These data shed light on a peculiar control mechanism of ß-glucosidase multimerization, an organization common to many defensive GH1 members.

Puzzling Out the Colchicine Biosynthetic Pathway.

Colchicine is among the oldest plant natural products (NPs) still used for treating a broad spectrum of human diseases including gout and other articular inflammation disorders. This molecule is synthesized by several herbaceous species related to the Liliaceae family, but in very low quantities in whole plants. As for many pharmaceutical compounds from plants, the production of colchicine still depends on the natural resource from which it is extracted. From the past decade, metabolic engineering has progressively become a credible alternative for the cost-effective large-scale production of several valuable NPs. In the same vein, Nett and colleagues recently reported an unprecedented advance in the field for colchicine. By using a combination of transcriptomics, metabolomics and pathway reconstitution, Sattely's group deciphered a near-complete biosynthetic pathway to colchicine without prior knowledge of biosynthetic genes. Besides constituting a benchmark for the elucidation of natural product biosynthetic pathways, it opens unprecedented perspectives regarding metabolic engineering of colchicine biosynthesis.

Identifying Genes Involved in alkaloid Biosynthesis in Vinca minor Through Transcriptomics and Gene Co-Expression Analysis.

The lesser periwinkle Vinca minor accumulates numerous monoterpene indole alkaloids (MIAs) including the vasodilator vincamine. While the biosynthetic pathway of MIAs has been largely elucidated in other Apocynaceae such as Catharanthus roseus, the counterpart in V. minor remains mostly unknown, especially for reactions leading to MIAs specific to this plant. As a consequence, we generated a comprehensive V. minor transcriptome elaborated from eight distinct samples including roots, old and young leaves exposed to low or high light exposure conditions. This optimized resource exhibits an improved completeness compared to already published ones. Through homology-based searches using C. roseus genes as bait, we predicted candidate genes for all common steps of the MIA pathway as illustrated by the cloning of a tabersonine/vincadifformine 16-O-methyltransferase (Vm16OMT) isoform. The functional validation of this enzyme revealed its capacity of methylating 16-hydroxylated derivatives of tabersonine, vincadifformine and lochnericine with a Km 0.94 ± 0.06 µM for 16-hydroxytabersonine. Furthermore, by combining expression of fusions with yellow fluorescent proteins and interaction assays, we established that Vm16OMT is located in the cytosol and forms homodimers. Finally, a gene co-expression network was performed to identify candidate genes of the missing V. minor biosynthetic steps to guide MIA pathway elucidation.

Transcriptomics of the Rooibos (Aspalathus linearis) Species Complex.

Rooibos (Aspalathus linearis), widely known as a herbal tea, is endemic to the Cape Floristic Region of South Africa (SA). It produces a wide range of phenolic compounds that have been associated with diverse health promoting properties of the plant. The species comprises several growth forms that differ in their morphology and biochemical composition, only one of which is cultivated and used commercially. Here, we established methodologies for non-invasive transcriptome research of wild-growing South African plant species, including (1) harvesting and transport of plant material suitable for RNA sequencing; (2) inexpensive, high-throughput biochemical sample screening; (3) extraction of high-quality RNA from recalcitrant, polysaccharide- and polyphenol rich plant material; and (4) biocomputational analysis of Illumina sequencing data, together with the evaluation of programs for transcriptome assembly (Trinity, IDBA-Trans, SOAPdenovo-Trans, CLC), protein prediction, as well as functional and taxonomic transcript annotation. In the process, we established a biochemically characterized sample pool from 44 distinct rooibos ecotypes (1-5 harvests) and generated four in-depth annotated transcriptomes (each comprising on average ≈86,000 transcripts) from rooibos plants that represent distinct growth forms and differ in their biochemical profiles. These resources will serve future rooibos research and plant breeding endeavours.

Visualization of Aspalathin in Rooibos (Aspalathus linearis) Plant and Herbal Tea Extracts Using Thin-Layer Chromatography.

Aspalathin, the main polyphenol of rooibos (Aspalathus linearis), is associated with diverse health promoting properties of the tea. During fermentation, aspalathin is oxidized and concentrations are significantly reduced. Standardized methods for quality control of rooibos products do not investigate aspalathin, since current techniques of aspalathin detection require expensive equipment and expertise. Here, we describe a simple and fast thin-layer chromatography (TLC) method that can reproducibly visualize aspalathin in rooibos herbal tea and plant extracts at a limit of detection (LOD) equal to 178.7 ng and a limit of quantification (LOQ) equal to 541.6 ng. Aspalathin is a rare compound, so far only found in A. linearis and its (rare) sister species A. pendula. Therefore, aspalathin could serve as a marker compound for authentication and quality control of rooibos products, and the described TLC method represents a cost-effective approach for high-throughput screening of plant and herbal tea extracts.

Novel phages of healthy skin metaviromes from South Africa.

Recent skin metagenomic studies have investigated the harbored viral diversity and its possible influence on healthy skin microbial populations, and tried to establish global patterns of skin-phage evolution. However, the detail associated with the phages that potentially play a role in skin health has not been investigated. While skin metagenome and -metavirome studies have indicated that the skin virome is highly site specific and shows marked interpersonal variation, they have not assessed the presence/absence of individual phages. Here, we took a semi-culture independent approach (metaviromic) to better understand the composition of phage communities on skin from South African study participants. Our data set adds over 130 new phage species of the skin to existing databases. We demonstrated that identical phages were present on different individuals and in different body sites, and we conducted a detailed analysis of the structural organization of these phages. We further found that a bacteriophage related to the Staphylococcus capitis phage Stb20 may be a common skin commensal virus potentially regulating its host and its activities on the skin.