Over 100 Million Years of Enzyme Evolution Underpinning the Production of Morphine in the Papaveraceae Family of Flowering Plants

Phylogeography and genetic structure of Papaver bracteatum populations in Iran based on genotyping-by-sequencing (GBS)

Papaver bracteatum, known for its high thebaine content and absence of morphine, has emerged as a promising alternative to opium poppy for codeine production. In this study, our objective was to create a diverse panel representing the natural variation of this species in Iran. To achieve this, we employed genotyping-by-sequencing to obtain genome-wide distributed single-nucleotide polymorphisms (SNPs) for phylogeographic analysis, population structure assessment, and evaluation of genetic diversity within P. bracteatum populations. A total of 244 P. bracteatum individuals from 13 distinct populations formed seven genetic groups, along with one highly admixed population. We observed a clear split between the populations inhabiting the Alborz Mts. in the east and Zagros Mts. in the west. In between these mountain ranges, the population of Kachal Mangan exhibited a high degree of genetic admixture between both genetic groups. At or after the end of the last glacial maximum, when climate conditions rapidly changed, all P. bracteatum populations experienced a strong demographic bottleneck reducing the already small effective population sizes further before they increased to their recent strengths. Our results suggest that the ongoing climate change together with human pressure on the species' habitats and limited seed-dispersal ability are potential factors contributing today to rising genetic isolation of P. bracteatum populations. Our results provide genetic data that can be used for conservation measures to safeguard the species' genetic diversity as a resource for future breeding approaches in this medicinally important species.

A cornucopia of diversity-Ranunculales as a model lineage

The Ranunculales are a hyperdiverse lineage in many aspects of their phenotype, including growth habit, floral and leaf morphology, reproductive mode, and specialized metabolism. Many Ranunculales species, such as opium poppy and goldenseal, have a high medicinal value. In addition, the order includes a large number of commercially important ornamental plants, such as columbines and larkspurs. The phylogenetic position of the order with respect to monocots and core eudicots and the diversity within this lineage make the Ranunculales an excellent group for studying evolutionary processes by comparative studies. Lately, the phylogeny of Ranunculales was revised, and genetic and genomic resources were developed for many species, allowing comparative analyses at the molecular scale. Here, we review the literature on the resources for genetic manipulation and genome sequencing, the recent phylogeny reconstruction of this order, and its fossil record. Further, we explain their habitat range and delve into the diversity in their floral morphology, focusing on perianth organ identity, floral symmetry, occurrences of spurs and nectaries, sexual and pollination systems, and fruit and dehiscence types. The Ranunculales order offers a wealth of opportunities for scientific exploration across various disciplines and scales, to gain novel insights into plant biology for researchers and plant enthusiasts alike.

Divergent proteomic profiles of opium poppy cultivars

We examined the proteomic profiles of three registered opium poppy cultivars (Papaver somniferum L.) with varying alkaloid contents. The study was conducted on both the stem and capsule organs. A high number of differentially expressed proteins (DEPs) were identified between the cultivars and the organs. We analyzed DEPs for their contribution in GO terms and KEGG pathways. The upregulated DEPs were significantly enriched in photosynthesis and translation for morphine-rich and noscapine-rich cultivars, respectively. The data indicated that photosynthesis is crucial for benzylisoquinoline alkaloid (BIA) biosynthesis, but different processes are also effective in morphine and noscapine biosynthesis, which occur at different branches in the biosynthetic pathway. The proteomics profiles revealed that energy demand is more effective in morphine biosynthesis, while translational control plays a leading role in noscapine biosynthesis. This study represents the first report demonstrating organ-based and cultivar-based protein expression differences in mature poppy plants.

From genomics to metabolomics: Deciphering sanguinarine biosynthesis in Dicranostigma leptopodum

Dicranostigma leptopodum (Maxim) Fedde (DLF) is a renowned medicinal plant in China, known to be rich in alkaloids. However, the unavailability of a reference genome has impeded investigation into its plant metabolism and genetic breeding potential. Here we present a high-quality chromosomal-level genome assembly for DLF, derived using a combination of Nanopore long-read sequencing, Illumina short-read sequencing and Hi-C technologies. Our assembly genome spans a size of 621.81 Mb with an impressive contig N50 of 93.04 Mb. We show that the species-specific whole-genome duplication (WGD) of DLF and Papaver somniferum corresponded to two rounds of WGDs of Papaver setigerum. Furthermore, we integrated comprehensive homology searching, gene family analyses and construction of a gene-to-metabolite network. These efforts led to the discovery of co-expressed transcription factors, including NAC and bZIP, alongside sanguinarine (SAN) pathway genes CYP719 (CFS and SPS). Notably, we identified P6H as a promising gene for enhancing SAN production. By providing the first reference genome for Dicranostigma, our study confirms the genomic underpinning of SAN biosynthesis and establishes a foundation for advancing functional genomic research on Papaveraceae species. Our findings underscore the pivotal role of high-quality genome assemblies in elucidating genetic variations underlying the evolutionary origin of secondary metabolites.

The CYP80A and CYP80G Are Involved in the Biosynthesis of Benzylisoquinoline Alkaloids in the Sacred Lotus (Nelumbo nucifera)

Bisbenzylisoquinoline and aporphine alkaloids are the two main pharmacological compounds in the ancient sacred lotus (Nelumbo nucifera). The biosynthesis of bisbenzylisoquinoline and aporphine alkaloids has attracted extensive attention because bisbenzylisoquinoline alkaloids have been reported as potential therapeutic agents for COVID-19. Our study showed that NnCYP80A can catalyze C-O coupling in both (R)-N-methylcoclaurine and (S)-N-methylcoclaurine to produce bisbenzylisoquinoline alkaloids with three different linkages. In addition, NnCYP80G catalyzed C-C coupling in aporphine alkaloids with extensive substrate selectivity, specifically using (R)-N-methylcoclaurine, (S)-N-methylcoclaurine, coclaurine and reticuline as substrates, but the synthesis of C-ring alkaloids without hydroxyl groups in the lotus remains to be elucidated. The key residues of NnCYP80G were also studied using the 3D structure of the protein predicted using Alphafold 2, and six key amino acids (G39, G69, A211, P288, R425 and C427) were identified. The R425A mutation significantly decreased the catalysis of (R)-N-methylcoclaurine and coclaurine inactivation, which might play important role in the biosynthesis of alkaloids with new configurations.

Exploring the Evolvability of Plant Specialized Metabolism: Uniqueness Out Of Uniformity and Uniqueness Behind Uniformity

The huge structural diversity exhibited by plant specialized metabolites has primarily been considered to result from the catalytic specificity of their biosynthetic enzymes. Accordingly, enzyme gene multiplication and functional differentiation through spontaneous mutations have been established as the molecular mechanisms that drive metabolic evolution. Nevertheless, how plants have assembled and maintained such metabolic enzyme genes and the typical clusters that are observed in plant genomes, as well as why identical specialized metabolites often exist in phylogenetically remote lineages, is currently only poorly explained by a concept known as convergent evolution. Here, we compile recent knowledge on the co-presence of metabolic modules that are common in the plant kingdom but have evolved under specific historical and contextual constraints defined by the physicochemical properties of each plant specialized metabolite and the genetic presets of the biosynthetic genes. Furthermore, we discuss a common manner to generate uncommon metabolites (uniqueness out of uniformity) and an uncommon manner to generate common metabolites (uniqueness behind uniformity). This review describes the emerging aspects of the evolvability of plant specialized metabolism that underlie the vast structural diversity of plant specialized metabolites in nature.

Evolutionary analysis of conserved non-coding elements subsequent to whole-genome duplication in opium poppy

Whole-genome duplication (WGD) leads to the duplication of both coding and non-coding sequences within an organism's genome, providing an abundant supply of genetic material that can drive evolution, ultimately contributing to plant adaptation and speciation. Although non-coding sequences contain numerous regulatory elements, they have been understudied compared to coding sequences. In order to address this gap, we explored the evolutionary patterns of regulatory sequences, coding sequences and transcriptomes using conserved non-coding elements (CNEs) as regulatory element proxies following the recent WGD event in opium poppy (Papaver somniferum). Our results showed similar evolutionary patterns in subgenomes of regulatory and coding sequences. Specifically, the biased or unbiased retention of coding sequences reflected the same pattern as retention levels in regulatory sequences. Further, the divergence of gene expression patterns mediated by regulatory element variations occurred at a more rapid pace than that of gene coding sequences. However, gene losses were purportedly dependent on relaxed selection pressure in coding sequences. Specifically, the rapid evolution of tissue-specific benzylisoquinoline alkaloid production in P. somniferum was associated with regulatory element changes. The origin of a novel stem-specific ACR, which utilized ancestral cis-elements as templates, is likely to be linked to the evolutionary trajectory behind the transition of the PSMT1-CYP719A21 cluster from high levels of expression solely in P. rhoeas root tissue to its elevated expression in P. somniferum stem tissue. Our findings demonstrate that rapid regulatory element evolution can contribute to the emergence of new phenotypes and provide valuable insights into the high evolvability of regulatory elements.

Genome-wide analysis of bHLH gene family in Coptis chinensis provides insights into the regulatory role in benzylisoquinoline alkaloid biosynthesis

Coptis chinensis Franch is a perennial species with high medical value. The rhizome of C. chinensis is a traditional Chinese medicine widely used for more than 2000 years in China. Its principal active ingredients are benzylisoquinoline alkaloids (BIAs). The basic helix-loop-helix (bHLH) transcription factors play an important regulatory role in the biosynthesis of plant secondary metabolites. However, the bHLH genes in C. chinensis have not been described, and little is known about their roles in alkaloid biosynthesis. In this study, a total of 143 CcbHLH genes (CcbHLHs) were identified and unevenly distributed on nine chromosomes. Phylogenetic analysis divided the 143 CcbHLH proteins into 26 subfamilies by comparison with Arabidopsis thaliana bHLH proteins. The majority CcbHLHs in each subgroup had similar gene structures and conserved motifs. Furthermore, the physicochemical properties, conserved motif, intron/exon composition, and cis-acting elements of CcbHLHs were analyzed. Transcriptome analysis revealed that 30 CcbHLHs were significantly expressed in the rhizomes of C. chinensis. Co-expression analysis revealed that 11 CcbHLHs were highly positively correlated with contents of various alkaloids of C. chinensis. Moreover, yeast one-hybrid experiments verified that CcbHLH001 and CcbHLH0002 could interact with the promoters of berberine biosynthesis pathway genes CcBBE and CcCAS, suggesting their regulatory roles in BIA biosynthesis. This study provides comprehensive insights into the bHLH gene family in C. chinensis and will support in-depth functional characterization of CcbHLHs involved in the regulation of protoberberine-type alkaloid biosynthesis.

Isolation and Characterization of the Genes Involved in the Berberine Synthesis Pathway in Asian Blue Cohosh, Caulophyllum robustum

Caulophyllum robustum, commonly named Asian blue cohosh, is a perennial herb in the family Berberidaceae. It has traditionally been used for folk medicine in China. We isolated berberine from the leaves, stem, roots, and fruits of C. robustum, and this is the first report on berberine in this species. Transcriptome analysis was conducted for the characterization of berberine biosynthesis genes in C. robustum, in which, all the genes for berberine biosynthesis were identified. From 40,094 transcripts, using gene ontology (GO) analysis, 26,750 transcripts were assigned their functions in the categories of biological process, molecular function, and cellular component. In the analysis of genes expressed in different tissues, the numbers of genes in the categories of intrinsic component of membrane and transferase activity were up-regulated in leaves versus stem. The berberine synthesis genes in C. robustum were characterized by phylogenetic analysis with corresponding genes from other berberine-producing species. The co-existence of genes from different plant families in the deepest branch subclade implies that the differentiation of berberine synthesis genes occurred early in the evolution of berberine-producing plants. Furthermore, the copy number increment of the berberine synthesis genes was detected at the species level.

California poppy (Eschscholzia californica), the Papaveraceae golden girl model organism for evodevo and specialized metabolism

California poppy or golden poppy (Eschscholzia californica) is the iconic state flower of California, with native ranges from Northern California to Southwestern Mexico. It grows well as an ornamental plant in Mediterranean climates, but it might be invasive in many parts of the world. California poppy was also highly prized by Native Americans for its medicinal value, mainly due to its various specialized metabolites, especially benzylisoquinoline alkaloids (BIAs). As a member of the Ranunculales, the sister lineage of core eudicots it occupies an interesting phylogenetic position. California poppy has a short-lived life cycle but can be maintained as a perennial. It has a comparatively simple floral and vegetative morphology. Several genetic resources, including options for genetic manipulation and a draft genome sequence have been established already with many more to come. Efficient cell and tissue culture protocols are established to study secondary metabolite biosynthesis and its regulation. Here, we review the use of California poppy as a model organism for plant genetics, with particular emphasis on the evolution of development and BIA biosynthesis. In the future, California poppy may serve as a model organism to combine two formerly separated lines of research: the regulation of morphogenesis and the regulation of secondary metabolism. This can provide insights into how these two integral aspects of plant biology interact with each other.

New Low Morphine Opium Poppy Genotype Obtained by TILLING Approach

The opium poppy's ability to produce various alkaloids is both useful and problematic. Breeding of new varieties with varying alkaloid content is therefore an important task. In this paper, the breeding technology of new low morphine poppy genotypes, based on a combination of a TILLING approach and single-molecule real-time NGS sequencing, is presented. Verification of the mutants in the TILLING population was obtained using RT-PCR and HPLC methods. Only three of the single-copy genes of the morphine pathway among the eleven genes were used for the identification of mutant genotypes. Point mutations were obtained only in one gene (CNMT) while an insertion was obtained in the other (SalAT). Only a few expected transition SNPs from G:C to A:T were obtained. In the low morphine mutant genotype, the production of morphine was decreased to 0.1% from 1.4% in the original variety. A comprehensive description of the breeding process, a basic characterization of the main alkaloid content, and a gene expression profile for the main alkaloid-producing genes is provided. Difficulties with the TILLING approach are also described and discussed.

Elucidation of the (R)-enantiospecific benzylisoquinoline alkaloid biosynthetic pathways in sacred lotus (Nelumbo nucifera)

Benzylisoquinoline alkaloids (BIAs) are a structurally diverse group of plant specialized metabolites found mainly in members of the order Ranunculales, including opium poppy (Papaver somniferum), for which BIA biosynthetic pathways leading to the critical drugs morphine, noscapine, and sanguinarine have been elucidated. Sacred lotus (Nelumbo nucifera), in the order Proteales, accumulates medicinal BIAs in the proaporphine, aporphine, and bisbenzylisoquinoline structural subgroups with a prevalence of R enantiomers, opposed to the dominant S configuration occurring in the Ranunculales. Nevertheless, distinctive BIA biosynthetic routes in sacred lotus have not been explored. In planta labeling experiments and in vitro assays with recombinant enzymes and plant protein extracts showed that dopamine and 4-hydroxyphenylacetaldehyde derived from L-tyrosine serve as precursors for the formation of (R,S)-norcoclaurine in sacred lotus, whereas only (R)-norcoclaurine byproducts are favored in the plant by action of R-enantiospecific methyltransferases and cytochrome P450 oxidoreductases (CYPs). Enzymes responsible for the R-enantiospecific formation of proaporphine (NnCYP80Q1) and bisbenzylisoquinoline (NnCYP80Q2) scaffolds, and a methylenedioxy bridge introduction on aporphine substrates (NnCYP719A22) were identified, whereas additional aspects of the biosynthetic pathways leading to the distinctive alkaloid profile are discussed. This work expands the availability of molecular tools that can be deployed in synthetic biology platforms for the production of high-value alkaloids.

Mining Therapeutic Efficacy from Treasure Chest of Biodiversity and Chemodiversity: Pharmacophylogeny of Ranunculales Medicinal Plants

Ranunculales, comprising of 7 families that are rich in medicinal species frequently utilized by traditional medicine and ethnomedicine, represents a treasure chest of biodiversity and chemodiversity. The phylogenetically related species often have similar chemical profile, which makes them often possess similar therapeutic spectrum. This has been validated by both ethnomedicinal experiences and pharmacological investigations. This paper summarizes molecular phylogeny, chemical constituents, and therapeutic applications of Ranunculales, i.e., a pharmacophylogeny study of this representative medicinal order. The phytochemistry/metabolome, ethnomedicine and bioactivity/pharmacology data are incorporated within the phylogenetic framework of Ranunculales. The most studied compounds of this order include benzylisoquinoline alkaloid, flavonoid, terpenoid, saponin and lignan, etc. Bisbenzylisoquinoline alkaloids are especially abundant in Berberidaceae and Menispermaceae. The most frequent ethnomedicinal uses are arthritis, heat-clearing and detoxification, carbuncle-abscess and sore-toxin. The most studied bioactivities are anticancer/cytotoxic, antimicrobial, and anti-inflammatory activities, etc. The pharmacophylogeny analysis, integrated with both traditional and modern medicinal uses, agrees with the molecular phylogeny based on chloroplast and nuclear DNA sequences, in which Ranunculales is divided into Ranunculaceae, Berberidaceae, Menispermaceae, Lardizabalaceae, Circaeasteraceae, Papaveraceae, and Eupteleaceae families. Chemical constituents and therapeutic efficacy of each taxonomic group are reviewed and the underlying connection between phylogeny, chemodiversity and clinical uses is revealed, which facilitate the conservation and sustainable utilization of Ranunculales pharmaceutical resources, as well as developing novel plant-based pharmacotherapy.

The genome of Corydalis reveals the evolution of benzylisoquinoline alkaloid biosynthesis in Ranunculales

Species belonging to the order Ranunculales have attracted much attention because of their phylogenetic position as a sister group to all other eudicot lineages and their ability to produce unique yet diverse benzylisoquinoline alkaloids (BIAs). The Papaveraceae family in Ranunculales is often used as a model system for studying BIA biosynthesis. Here, we report the chromosome-level genome assembly of Corydalis tomentella, a species of Fumarioideae, one of the two subfamilies of Papaveraceae. Based on comparisons of sequenced Ranunculalean species, we present clear evidence of a shared whole-genome duplication (WGD) event that has occurred before the divergence of Ranunculales but after its divergence from other eudicot lineages. The C. tomentella genome enabled us to integrate isotopic labeling and comparative genomics to reconstruct the BIA biosynthetic pathway for both sanguinarine biosynthesis shared by papaveraceous species and the cavidine biosynthesis that is specific to Corydalis. Also, our comparative analysis revealed that gene duplications, especially tandem gene duplications, underlie the diversification of BIA biosynthetic pathways in Ranunculales. In particular, tandemly duplicated berberine bridge enzyme-like genes appear to be involved in cavidine biosynthesis. In conclusion, our study of the C. tomentella genome provides important insights into the occurrence of WGDs during the early evolution of eudicots, as well as into the evolution of BIA biosynthesis in Ranunculales.

A functionally conserved STORR gene fusion in Papaver species that diverged 16.8 million years ago

The STORR gene fusion event is considered essential for the evolution of the promorphinan/morphinan subclass of benzylisoquinoline alkaloids (BIAs) in opium poppy as the resulting bi-modular protein performs the isomerization of (S)- to (R)-reticuline essential for their biosynthesis. Here, we show that of the 12 Papaver species analysed those containing the STORR gene fusion also contain promorphinans/morphinans with one important exception. P. californicum encodes a functionally conserved STORR but does not produce promorphinans/morphinans. We also show that the gene fusion event occurred only once, between 16.8-24.1 million years ago before the separation of P. californicum from other Clade 2 Papaver species. The most abundant BIA in P. californicum is (R)-glaucine, a member of the aporphine subclass of BIAs, raising the possibility that STORR, once evolved, contributes to the biosynthesis of more than just the promorphinan/morphinan subclass of BIAs in the Papaveraceae.

Endogenous opiates and behavior: 2020

This paper is the forty-third consecutive installment of the annual anthological review of research concerning the endogenous opioid system, summarizing articles published during 2020 that studied the behavioral effects of molecular, pharmacological and genetic manipulation of opioid peptides and receptors as well as effects of opioid/opiate agonists and antagonists. The review is subdivided into the following specific topics: molecular-biochemical effects and neurochemical localization studies of endogenous opioids and their receptors (1), the roles of these opioid peptides and receptors in pain and analgesia in animals (2) and humans (3), opioid-sensitive and opioid-insensitive effects of nonopioid analgesics (4), opioid peptide and receptor involvement in tolerance and dependence (5), stress and social status (6), learning and memory (7), eating and drinking (8), drug abuse and alcohol (9), sexual activity and hormones, pregnancy, development and endocrinology (10), mental illness and mood (11), seizures and neurologic disorders (12), electrical-related activity and neurophysiology (13), general activity and locomotion (14), gastrointestinal, renal and hepatic functions (15), cardiovascular responses (16), respiration and thermoregulation (17), and immunological responses (18).

Pervasive misannotation of microexons that are evolutionarily conserved and crucial for gene function in plants

It is challenging to identify the smallest microexons (≤15-nt) due to their small size. Consequently, these microexons are often misannotated or missed entirely during genome annotation. Here, we develop a pipeline to accurately identify 2,398 small microexons in 10 diverse plant species using 990 RNA-seq datasets, and most of them have not been annotated in the reference genomes. Analysis reveals that microexons tend to have increased detained flanking introns that require post-transcriptional splicing after polyadenylation. Examination of 45 conserved microexon clusters demonstrates that microexons and associated gene structures can be traced back to the origin of land plants. Based on these clusters, we develop an algorithm to genome-wide model coding microexons in 132 plants and find that microexons provide a strong phylogenetic signal for plant organismal relationships. Microexon modeling reveals diverse evolutionary trajectories, involving microexon gain and loss and alternative splicing. Our work provides a comprehensive view of microexons in plants.

The small size (≤15-nt) of micorexons poses difficulties for genome annotation and identification using standard RNA sequence mapping approaches. Here, the authors develop computational pipelines to discover and predict microexons in plants and reveal diverse evolutionary trajectories via genomewide microexon modeling.

Distribution of Therapeutic Efficacy of Ranunculales Plants Used by Ethnic Minorities on the Phylogenetic Tree of Chinese Species

The medicinal properties of plants can be evolutionarily predicted by phylogeny-based methods, which, however, have not been used to explore the regularity of therapeutic effects of Chinese plants utilized by ethnic minorities. This study aims at exploring the distribution law of therapeutic efficacy of Ranunculales plants on the phylogenetic tree of Chinese species. We collected therapeutic efficacy data of 551 ethnomedicinal species belonging to five species-rich families of Ranunculales; these therapeutic data were divided into 15 categories according to the impacted tissues and organs. The phylogenetic tree of angiosperm species was used to analyze the phylogenetic signals of ethnomedicinal plants by calculating the net relatedness index (NRI) and nearest taxon index (NTI) in R language. The NRI results revealed a clustered structure for eight medicinal categories (poisoning/intoxication, circulatory, gastrointestinal, eyesight, oral, pediatric, skin, and urinary disorders) and overdispersion for the remaining seven (neurological, general, hepatobiliary, musculoskeletal, otolaryngologic, reproductive, and respiratory disorders), while the NTI metric identified the clustered structure for all. Statistically, NRI and NTI values were significant in 5 and 11 categories, respectively. It was found that Mahonia eurybracteata has therapeutic effects on all categories. iTOL was used to visualize the distribution of treatment efficacy on species phylogenetic trees. By figuring out the distribution of therapeutic effects of Ranunculales medicinal plants, the importance of phylogenetic methods in finding potential medicinal resources is highlighted; NRI, NTI, and similar indices can be calculated to help find taxonomic groups with medicinal efficacy based on the phylogenetic tree of flora in different geographic regions.

Dissection of full-length transcriptome and metabolome of Dichocarpum (Ranunculaceae): implications in evolution of specialized metabolism of Ranunculales medicinal plants

Several main families of Ranunculales are rich in alkaloids and other medicinal compounds; many species of these families are used in traditional and folk medicine. Dichocarpum is a representative medicinal genus of Ranunculaceae, but the genetic basis of its metabolic phenotype has not been investigated, which hinders its sustainable conservation and utilization. We use the third-generation high-throughput sequencing and metabolomic techniques to decipher the full-length transcriptomes and metabolomes of five Dichocarpum species endemic in China, and 71,598 non-redundant full-length transcripts were obtained, many of which are involved in defense, stress response and immunity, especially those participating in the biosynthesis of specialized metabolites such as benzylisoquinoline alkaloids (BIAs). Twenty-seven orthologs extracted from trancriptome datasets were concatenated to reconstruct the phylogenetic tree, which was verified by the clustering analysis based on the metabolomic profile and agreed with the Pearson correlation between gene expression patterns of Dichocarpum species. The phylogenomic analysis of phytometabolite biosynthesis genes, e.g., (S)-norcoclaurine synthase, methyltransferases, cytochrome p450 monooxygenases, berberine bridge enzyme and (S)-tetrahydroprotoberberine oxidase, revealed the evolutionary trajectories leading to the chemodiversity, especially that of protoberberine type, aporphine type and bis-BIA abundant in Dichocarpum and related genera. The biosynthesis pathways of these BIAs are proposed based on full-length transcriptomes and metabolomes of Dichocarpum. Within Ranunculales, the gene duplications are common, and a unique whole genome duplication is possible in Dichocarpum. The extensive correlations between metabolite content and gene expression support the co-evolution of various genes essential for the production of different specialized metabolites. Our study provides insights into the transcriptomic and metabolomic landscapes of Dichocarpum, which will assist further studies on genomics and application of Ranunculales plants.

Three chromosome-scale Papaver genomes reveal punctuated patchwork evolution of the morphinan and noscapine biosynthesis pathway

For millions of years, plants evolve plenty of structurally diverse secondary metabolites (SM) to support their sessile lifestyles through continuous biochemical pathway innovation. While new genes commonly drive the evolution of plant SM pathway, how a full biosynthetic pathway evolves remains poorly understood. The evolution of pathway involves recruiting new genes along the reaction cascade forwardly, backwardly, or in a patchwork manner. With three chromosome-scale Papaver genome assemblies, we here reveal whole-genome duplications (WGDs) apparently accelerate chromosomal rearrangements with a nonrandom distribution towards SM optimization. A burst of structural variants involving fusions, translocations and duplications within 7.7 million years have assembled nine genes into the benzylisoquinoline alkaloids gene cluster, following a punctuated patchwork model. Biosynthetic gene copies and their total expression matter to morphinan production. Our results demonstrate how new genes have been recruited from a WGD-induced repertoire of unregulated enzymes with promiscuous reactivities to innovate efficient metabolic pathways with spatiotemporal constraint.

Structural studies of codeinone reductase reveal novel insights into aldo-keto reductase function in benzylisoquinoline alkaloid biosynthesis

Benzylisoquinoline alkaloids (BIAs) are a class of specialized metabolites with a diverse range of chemical structures and physiological effects. Codeine and morphine are two closely related BIAs with particularly useful analgesic properties. The aldo-keto reductase (AKR) codeinone reductase (COR) catalyzes the final and penultimate steps in the biosynthesis of codeine and morphine, respectively, in opium poppy (Papaver somniferum). However, the structural determinants that mediate substrate recognition and catalysis are not well defined. Here, we describe the crystal structure of apo-COR determined to a resolution of 2.4 Å by molecular replacement using chalcone reductase as a search model. Structural comparisons of COR to closely related plant AKRs and more distantly related homologues reveal a novel conformation in the β1α1 loop adjacent to the BIA-binding pocket. The proximity of this loop to several highly conserved active-site residues and the expected location of the nicotinamide ring of the NADP(H) cofactor suggest a model for BIA recognition that implies roles for several key residues. Using site-directed mutagenesis, we show that substitutions at Met-28 and His-120 of COR lead to changes in AKR activity for the major and minor substrates codeinone and neopinone, respectively. Our findings provide a framework for understanding the molecular basis of substrate recognition in COR and the closely related 1,2-dehydroreticuline reductase responsible for the second half of a stereochemical inversion that initiates the morphine biosynthesis pathway.

Analysis of the Coptis chinensis genome reveals the diversification of protoberberine-type alkaloids

Chinese goldthread (Coptis chinensis Franch.), a member of the Ranunculales, represents an important early-diverging eudicot lineage with diverse medicinal applications. Here, we present a high-quality chromosome-scale genome assembly and annotation of C. chinensis. Phylogenetic and comparative genomic analyses reveal the phylogenetic placement of this species and identify a single round of ancient whole-genome duplication (WGD) shared by the Ranunculaceae. We characterize genes involved in the biosynthesis of protoberberine-type alkaloids in C. chinensis. In particular, local genomic tandem duplications contribute to member amplification of a Ranunculales clade-specific gene family of the cytochrome P450 (CYP) 719. The functional versatility of a key CYP719 gene that encodes the (S)-canadine synthase enzyme involved in the berberine biosynthesis pathway may play critical roles in the diversification of other berberine-related alkaloids in C. chinensis. Our study provides insights into the genomic landscape of early-diverging eudicots and provides a valuable model genome for genetic and applied studies of Ranunculales.

Synthesis of Morphinans through Anodic Aryl-Aryl Coupling

The morphinans are an important class of structurally fascinating and physiologically important natural products as exemplified by the famous opium alkaloids of the morphine family. Although this class of secondary metabolites from the juice of the opium poppy capsule was already used for medicinal purposes thousands of years ago, chemical modifications are still being applied to the core structure today in order to achieve the most specific effect on the various receptor subtypes possible with the fewest possible side effects. The unusual architecture of the morphinan core has also proven to be a highly challenging target for total synthesis. This review highlights electrosynthetic approaches towards natural and semisynthetic morphinan alkaloids. The historical progress in applying anodic aryl-aryl couplings to the construction of the morphinan framework is described in chronological order while particular benefits and challenges concerning the electrochemical transformations are grouped together, including the influence of substitution patterns, protecting groups, and reaction conditions.

Comparative analysis using the draft genome sequence of California poppy (Eschscholzia californica) for exploring the candidate genes involved in benzylisoquinoline alkaloid biosynthesis

Genome characterization of California poppy (Eschscholzia californica cv. “Hitoezaki”), which produces pharmaceutically important benzylisoquinoline alkaloids (BIAs), was carried out using the draft genome sequence. The numbers of tRNA and rRNA genes were close to those of the other plant species tested, whereas the frequency of repetitive sequences was distinct from those species. Comparison of the predicted genes with those of Amborella trichopoda, Nelumbo nucifera, Solanum lycopersicum, and Arabidopsis thaliana, and analyses of gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway indicated that the enzyme genes involved in BIA biosynthesis were highly enriched in the California poppy genome. Further comparative analysis using the genome information of Papaver somniferum and Aquilegia coerulea, both BIA-producing plants, revealed that many genes encoding BIA biosynthetic enzymes, transcription factors, transporters, and candidate proteins, possibly related to BIA biosynthesis, were specifically distributed in these plant species.

PR10/Bet v1-like Proteins as Novel Contributors to Plant Biochemical Diversity

Pathogenesis-related (PR) proteins constitute a broad class of plant proteins with analogues found throughout nature from bacteria to higher eukaryotes. PR proteins were first noted in plants as part of the hypersensitive response, but have since been assigned an array of biological roles. The PR10/Bet v1-like proteins are a subset of PR proteins characterized by an ability to bind a wide range of lipophilic ligands, uniquely positioning them as contributors to specialized biosynthetic pathways. PR10/Bet v1-like proteins participate in the production of plant alkaloids and phenolics including flavonoids, both as general binding proteins and in special cases as catalysts. Owing initially to the perceived allergenic properties of PR10/Bet v1-like proteins, many were studied at the structural level to elucidate the basis for ligand binding. These studies provided a foundation for more recent efforts to understand higher-level structural order and how PR10/Bet v1-like proteins catalyse key reactions in plant pathways. Synthetic biology aimed at reconstituting plant-specialized metabolism in microorganisms uses knowledge of these proteins to fine-tune performance in new systems.

The scaffold-forming steps of plant alkaloid biosynthesis

Alkaloids from plants are characterised by structural diversity and bioactivity, and maintain a privileged position in both modern and traditional medicines. In recent years, there have been significant advances in elucidating the biosynthetic origins of plant alkaloids. In this review, I will describe the progress made in determining the metabolic origins of the so-called true alkaloids, specialised metabolites derived from amino acids containing a nitrogen heterocycle. By identifying key biosynthetic steps that feature in the majority of pathways, I highlight the key roles played by modifications to primary metabolism, iminium reactivity and spontaneous reactions in the molecular and evolutionary origins of these pathways.