Dear Editor, Rauvolfia tetraphylla(aka the Devil pepper)(Supplemental Figure 1)is a well-known medicinal plant that produces monoterpenoid indole alkaloids(MIAs).This MIA biosynthesis occurs in several or-gans,including leaves,stems,fruit,and roots,which accumulate the famous antiarrhythmic ajmaline(Kumar et al.2016a,2016b;Kumara et al.,2019).

The composition of bioactive polyphenols from grape canes, an important viticultural byproduct, was shown to be varietal-dependent; however, the influence of soil-related terroir factors remains unexplored. Using spatial metabolomics and correlation-based networks, we investigated how continuous changes in soil features and topography may impact the polyphenol composition in grape canes. Soil properties, topography, and grape cane extracts were analyzed at georeferenced points over 3 consecutive years, followed by UPLC-DAD-MS-based metabolomic analysis targeting 42 metabolites. Principal component analyses on intra-vintage metabolomic data presented a good reproducibility in relation to geographic coordinates. A correlation-driven approach was used to explore the combined influence of soil and topographic variables on metabolomic responses. As a result, a metabolic cluster including flavonoids was correlated with elevation and curvature. Spatial metabolomics driven by correlation-based networks represents a powerful approach to spatialize field-omics data and may serve as new field-phenotyping tool in precision agriculture.

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.

Monoterpenoid indole alkaloids (MIAs) represent a large class of plant natural products with marketed pharmaceutical activities against a wide range of indications, including cancer, malaria and hypertension. Halogenated MIAs have shown improved pharmaceutical properties; however, synthesis of new-to-nature halogenated MIAs remains a challenge. Here we demonstrate a platform for de novo biosynthesis of two MIAs, serpentine and alstonine, in baker's yeast Saccharomyces cerevisiae and deploy it to systematically explore the biocatalytic potential of refactored MIA pathways for the production of halogenated MIAs. From this, we demonstrate conversion of individual haloindole derivatives to a total of 19 different new-to-nature haloserpentine and haloalstonine analogs. Furthermore, by process optimization and heterologous expression of a modified halogenase in the microbial MIA platform, we document de novo halogenation and biosynthesis of chloroalstonine. Together, this study highlights a microbial platform for enzymatic exploration and production of complex natural and new-to-nature MIAs with therapeutic potential.

Catharanthus_roseus_v3.1_UT.cds (Predicted CDS). Catharanthus_roseus_v3.1_UT.gff (Genome annotation file (GFF)). Catharanthus_roseus_v3.1_UT.pep (Predicted proteins). Catharanthus_roseus_v3.1_UT.tr (Predicted transcripts). Cuello et al - F1000R - SuppMat.xlsx (Supplementary tables).

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.

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.

pPPO16, the first Ea-inducible promoter cloned from apple, can be a useful component of intragenic strategies to create fire blight resistant apple genotypes. Intragenesis is an important alternative to transgenesis to produce modified plants containing native DNA only. A key point to develop such a strategy is the availability of regulatory sequences controlling the expression of the gene of interest. With the aim of finding apple gene promoters either inducible by the fire blight pathogen Erwinia amylovora (Ea) or moderately constitutive, we focused on polyphenoloxidase genes (PPO). These genes encode oxidative enzymes involved in many physiological processes and have been previously shown to be upregulated during the Ea infection process. We found ten PPO and two PPO-like sequences in the apple genome and characterized the promoters of MdPPO16 (pPPO16) and MdKFDV02 PPO-like (pKFDV02) for their potential as Ea-inducible and low-constitutive regulatory sequences, respectively. Expression levels of reporter genes fused to these promoters and transiently or stably expressed in apple were quantified after various treatments. Unlike pKFDV02 which displayed a variable activity, pPPO16 allowed a fast and strong expression of transgenes in apple following Ea infection in a Type 3 Secretion System dependent manner. Altogether our results does not confirmed pKFDV02 as a constitutive and weak promoter whereas pPPO16, the first Ea-inducible promoter cloned from apple, can be a useful component of intragenic strategies to create fire blight resistant apple genotypes.

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.

DNA methylation is the most studied epigenetic mark in both plants and animals. The gold standard for assaying genome-wide DNA methylation at single-base resolution is whole-genome bisulfite sequencing (WGBS). Here, we describe an improved procedure for WGBS and original bioinformatic workflows applied to unravel tissue-specific variations of the methylome in relation to gene expression and accumulation of secondary metabolites in the medicinal plant Catharanthus roseus.

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.

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.

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.

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).

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.

Trihydroxycinnamoyl spermidines (THCSpd) are plant specialized metabolites with promising pharmacological activities as antifungals, antibacterial, antiviral, and antidepressant drugs. However, their characterization and potential pharmaceutical exploitation are greatly impaired by the sourcing of these compounds, restricted to the pollen of core Eudicot plant species. In this work, we developed a precursor-directed biosynthesis of THCSpd in yeast using a dual enzymatic system based on 4-coumarate-CoA ligases (4CL) and spermidine N-hydroxycinnamoyltransferases (SHT). The system relies on the yeast endogenous spermidine pool and only requires hydroxycinnamic acids as exogenous precursors. By exploring 4CL isoforms and SHT diversity among plants, we have driven the production of 8 natural THCSpd, using single or mixed hydroxycinnamic acid precursors. Substrate promiscuities of 4CL and SHT were genuinely exploited to produce 8 new-to-nature THCSpd from exotic hydroxycinnamic and dihydrohydroxycinnamic acids, together with 3 new-to-nature THCSpd containing halogenated hydroxycinnamoyl moieties. In this work, we established a versatile and modular biotechnological production platform allowing the tailor-made THCSpd synthesis, constituting pioneer metabolic engineering for access to these valuable natural products.

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.

Grape canes are viticulture byproducts representing a sustainable source of valuable bioactive polyphenols. However, varietal origin and postharvest treatment greatly influence the effective concentration of biomolecules, thus limiting industrial development. With the aim to develop grape cane extracts with high polyphenol contents, the selection of a specific grape variety combined with optimized postharvest treatment is the major determinant. A previously described postharvest treatment comprising cutting grapevine stalks of 0.5 cm length and storage at 15-20 degrees C over 2 weeks was applied on a selection of 44 grape varieties representative of the genetic diversity of the whole European collection and performing the screening of polyphenol contents. Varietal rankings according to major polyphenols (catechin, epicatechin, E-resveratrol, E-piceatannol, E-epsilon-viniferin, E-miyabenol C, ampelopsin A, E-vitisin B, hopeaphenol, and isophopeaphenol) were performed with and without the postharvest treatment. We observed that postharvest treatment greatly influenced the total polyphenol composition but also the ranking of polyphenol-rich varieties. This polyphenol screening of grape canes from a large collection of European varieties revealed the importance of postharvest treatment together with the selection of varieties to develop natural extracts based on grapevine wood biomass enriched with molecules with health benefits.

Today, the genus Scedosporium comprises at least ten species with four of them, Scedosporium apiospermum, Scedosporium boydii, Scedosporium aurantiacum and Scedosporium minutisporum capable of colonizing the lungs of patients with cystic fibrosis. Scedosporium dehoogii, which is also common in the soil, has never been reported as causing human pulmonary infections. Here we report the first genome sequence for S. dehoogii, an invaluable resource to understand the genetic bases of pathogenesis in the genus Scedosporium.