Plant genomes remain highly fragmented and are often characterized by hundreds to thousands of assembly gaps. Here, we report chromosome-level reference and phased genome assembly of Ophiorrhiza pumila, a camptothecin-producing medicinal plant, through an ordered multi-scaffolding and experimental validation approach. With 21 assembly gaps and a contig N50 of 18.49 Mb, Ophiorrhiza genome is one of the most complete plant genomes assembled to date. We also report 273 nitrogen-containing metabolites, including diverse monoterpene indole alkaloids (MIAs). A comparative genomics approach identifies strictosidine biogenesis as the origin of MIA evolution. The emergence of strictosidine biosynthesis-catalyzing enzymes precede downstream enzymes’ evolution post γ whole-genome triplication, which occurred approximately 110 Mya in O. pumila, and before the whole-genome duplication in Camptotheca acuminata identified here. Combining comparative genome analysis, multi-omics analysis, and metabolic gene-cluster analysis, we propose a working model for MIA evolution, and a pangenome for MIA biosynthesis, which will help in establishing a sustainable supply of camptothecin.

No abstract

Aconitum carmichaelii is an important medicinal herb used widely in China, Japan, India, Korea, and other Asian countries. While extensive research on the characterization of metabolic extracts of A. carmichaelii has shown accumulation of numerous bioactive metabolites including aconitine and aconitine-type diterpene alkaloids, its biosynthetic pathway remains largely unknown. Biosynthesis of these secondary metabolites is tightly controlled and mostly occurs in a tissue-specific manner; therefore, transcriptome analysis across multiple tissues is an attractive method to identify the molecular components involved for further functional characterization. In order to understand the biosynthesis of secondary metabolites, Illumina-based deep transcriptome profiling and analysis was performed for four tissues (flower, bud, leaf, and root) of A. carmichaelii, resulting in 5.5 Gbps clean RNA-seq reads assembled into 128,183 unigenes. Unigenes annotated as possible rate-determining steps of an aconitine-type biosynthetic pathway were highly expressed in the root, in accordance with previous reports describing the root as the accumulation site for these metabolites. We also identified 21 unigenes annotated as cytochrome P450s and highly expressed in roots, which represent candidate unigenes involved in the diversification of secondary metabolites. Comparative transcriptome analysis of A. carmichaelii with A. heterophyllum identified 20,232 orthogroups, representing 30,633 unigenes of A. carmichaelii, gene ontology enrichment analysis of which revealed essential biological process together with a secondary metabolic process to be highly enriched. Unigenes identified in this study are strong candidates for aconitine-type diterpene alkaloid biosynthesis, and will serve as useful resources for further validation studies.

Secondary metabolites play a key role in coordinating ecology and defense strategies of plants. Diversity of these metabolites arises by conjugation of core structures with diverse chemical moieties, such as sugars in glycosylation. Active pools of phytohormones, including those involved in plant stress response, are also regulated by glycosylation. While much is known about the enzymes involved in glycosylation, we know little about their regulation or coordination with other processes. We characterized the flavonoid pathway transcription factor TRANSPARENT TESTA8 (TT8) in Arabidopsis (Arabidopsis thaliana) using an integrative omics strategy. This approach provides a systems-level understanding of the cellular machinery that is used to generate metabolite diversity by glycosylation. Metabolomics analysis of TT8 loss-of-function and inducible overexpression lines showed that TT8 coordinates glycosylation of not only flavonoids, but also nucleotides, thus implicating TT8 in regulating pools of activated nucleotide sugars. Transcriptome and promoter network analyses revealed that the TT8 regulome included sugar transporters, proteins involved in sugar binding and sequestration, and a number of carbohydrate-active enzymes. Importantly, TT8 affects stress response, along with brassinosteroid and jasmonic acid biosynthesis, by directly binding to the promoters of key genes of these processes. This combined effect on metabolite glycosylation and stress hormones by TT8 inducible overexpression led to significant increase in tolerance toward multiple abiotic and biotic stresses. Conversely, loss of TT8 leads to increased sensitivity to these stresses. Thus, the transcription factor TT8 is an integrator of secondary metabolism and stress response. These findings provide novel approaches to improve broad-spectrum stress tolerance.