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

Systems and synthetic biology for plant natural product pathway elucidation

Plants are one of the major reservoirs of medicinal compounds, serving as a cornerstone of both traditional and modern medicine. However, despite their importance, the complex biosynthetic pathways of many plant-derived compounds remain only partially understood, hindering their full potential in therapeutic applications. This review paper summarizes the advances in systems and synthetic biology utilized in the characterization and engineering of plant metabolic pathways. We discuss various strategies such as (1) co-expression analysis, (2) gene cluster identification, (3) metabolite profiling, (4) deep learning approaches, (5) genome-wide association studies, and (6) protein complex identification. Through case studies on several biosynthesis pathways, we highlight how these methods are applied to unravel complex pathways and enhance the production of important natural products. Finally, we discuss future directions in the context of metabolic engineering, including metabolon engineering, AI integration, and sustainable production strategies, underscoring the potential for cheaper and greener production of plant natural products.

Diversified quantity, gene structure, and expression profile of OPR gene family of A. annua

Artemisia annua, the source of artemisinin production, is a traditional herb used for treating malaria for thousand years. The genetic background is of high heterozygosity and traits (plant height, biomass, artemisinin content, etc.) are diverse across different germplasms. Unraveling the key genes associated with growth and secondary metabolism is essential for the efficient production of artemisinin. The 12-oxo-phytodienoic acid reductase (OPR) genes, crucial for plant growth and development and stress resistance, remain unexplored in A. annua. In this study, nine OPR genes (named as AaOPR1 to AaOPR9) were identified in A. annua, including two pairs of genes formed from recent tandem duplications. The number of OPRs varied among different haplotype genomes, and each OPR gene exhibiting distinct expression pattern. Moreover, the OPR family displayed evolutionarily activity with significant variations in numbers and gene structures observed across different plant species. Widespread gene duplication of OPRs, observed in the majority of analyzed plant genomes, brought evolutionary potential. DBR2, a member of AaOPRs involved in artemisinin biosynthesis, had two copies (AaOPR1/DBR2.1 and AaOPR2/DBR2.2) with different expression patterns, one of which was a recently generated copy with a significant 7-amino acids truncation. Heterologous protein expression and functional characterization of the two copies of DBR2 yielded multiple isomers with identical molecular weights but different arrangements, indicating neofunctionalization of the newly generated copy. The polymorphism within the OPR gene family merely scratches the surface of the genetic diversity in A. annua, and further investigation of genetic features is needed for the screening of elite germplasm resources.

The origins and spread of the opium poppy (Papaver somniferum L.) revealed by genomics and seed morphometrics

The opium poppy (Papaver somniferum L.) is one of the most important plants in human history. It is the main source of opiates used as analgesic medicines or psychotropic drugs, the latter related to addiction problems, illegal trafficking and geopolitical issues. Poppyseed is also used in cooking. The prehistoric origins, domestication and cultivation spread of the opium poppy remain unresolved. Traditionally, Papaver setigerum has been considered the wild ancestor with early cultivation presumed to have occurred in the Western Mediterranean region, where setigerum is autochthonous. Other theories suggest that somniferum may have been introduced by Southwest Asian early farmers as a weed. To investigate these hypotheses, we analysed 190 accessions from 15 Papaver species using genotype-by-sequencing and geometric morphometric (GMM) techniques. Our analysis revealed that setigerum is the only taxa genetically close to somniferum and can be better described as a subspecies. The domesticated plants are, however, distinct from setigerum. Additionally, GMM analysis of seeds also revealed morphological differences between setigerum and somniferum. Some phenotypically wild setigerum accessions exhibited intermediate genetic features, suggesting introgression events. Two major populations were found in somniferum and, to some extent, these correspond to differences in seed form. These two populations may reflect recent attempts to breed varieties rich in opiates, as opposed to varieties used for poppyseed production. This study supports the idea that opium poppy cultivation began in the Western Mediterranean, with setigerum as the wild progenitor, although some wild varieties are likely to be feral forms, which can confound domestication studies.This article is part of the theme issue 'Unravelling domestication: multi-disciplinary perspectives on human and non-human relationships in the past, present and future'.

Application and Progress of Genomics in Deciphering the Genetic Regulation Mechanisms of Plant Secondary Metabolites

This review aims to systematically dissect the genetic regulatory mechanisms of plant secondary metabolites in the era of genomics, while comprehensively summarizing the progress and potential impact of genomics in plant secondary metabolism research. By integrating methodologies such as high-throughput sequencing, structural genomics, comparative genomics, and functional genomics, we elucidate the principles underlying plant secondary metabolism and identify functional genes. The application of these technologies has deepened our understanding of secondary metabolic pathways and driven advancements in plant molecular genetics and genomics. The development of genomics has enabled scientists to gain profound insights into the biosynthetic pathways of secondary metabolites in plants such as ginseng (Panax ginseng) and grapevine (Vitis vinifera), while offering novel possibilities for precise regulation of these pathways. Despite remarkable progress in studying the genetic regulation of plant secondary metabolites, significant challenges persist. Future research must focus on integrating multi-omics data, developing advanced bioinformatics tools, and exploring effective genetic improvement strategies to fully harness the medicinal potential of plants and enhance their capacity to synthesize secondary metabolites.

SOI: robust identification of orthologous synteny with the Orthology Index and broad applications in evolutionary genomics

With the explosive growth of whole-genome datasets, accurate detection of orthologous synteny has become crucial for reconstructing evolutionary history. However, current methods for identifying orthologous synteny face great limitations, particularly in scaling with varied polyploidy histories and accurately removing out-paralogous synteny. In this study, we developed a scalable and robust approach, based on the Orthology Index (OI), to effectively identify orthologous synteny. Our evaluation across a large-scale empirical dataset with diverse polyploidization events demonstrated the high reliability and robustness of the OI method. Simulation-based benchmarks further validated the accuracy of our method, showing its superior performance against existing methods across a wide range of scenarios. Additionally, we explored its broad applications in reconstructing the evolutionary histories of plant genomes, including the inference of polyploidy, identification of reticulation, and phylogenomics. In conclusion, OI offers a robust, interpretable, and scalable approach for identifying orthologous synteny, facilitating more accurate and efficient analyses in plant evolutionary genomics.

The chromosome-scale genomes of two Tinospora species reveal differential regulation of the MEP pathway in terpenoid biosynthesis

BACKGROUND

The relationship between gene family expansion and the resulting changes in plant phenotypes has shown remarkable complexity during the evolution. The gene family expansion has contributed to the diversity in plant phenotypes, specifically metabolites through neo-functionalization and sub-functionalization. However, the negative regulatory effects associated with the gene family expansion remain poorly understood.

RESULTS

Here, we present the chromosome-scale genomes of Tinospora crispa and Tinospora sinensis. Comparative genomic analyses demonstrated conserved chromosomal evolution within the Menispermaceae family. KEGG analysis revealed a significant enrichment of genes related to terpenoid biosynthesis in T. sinensis. However, T. crispa exhibited a higher abundance of terpenoids compared to T. sinensis. Detailed analysis revealed the expansion of genes encoding 1-hydroxy-2-methyl 2-(E)-butenyl 4-diphosphate synthase (HDS), a key enzyme in the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway of terpenoid biosynthesis in T. sinensis. TsiHDS4 retained the ancestral function of converting methylerythritol cyclic diphosphate (MEcPP) to (E)-4-hydroxy-3-methylbut-2-enyl diphosphate (HMBPP). However, the noncanonical CDS-derived small peptide TsiHDS5 was shown to interact with TsiHDS4, inhibiting its catalytic activity. This interaction reduced the levels of HMBPP and isopentenyl pyrophosphate (IPP), which represent key substrates for downstream terpenoid biosynthesis.

CONCLUSIONS

These findings offer clues to decipher the variations in the MEP pathway of terpenoid biosynthesis between T. crispa and T. sinensis and form a basis for further detailed research on the negative regulation of expanded genes.

Tracing the Origin and Evolution of the Fungal Mycophenolic Acid Biosynthesis Pathway

Like bacteria and plants, fungi produce a remarkable diversity of small molecules with potent activities for human health known as natural products or secondary metabolites. One such example is mycophenolic acid, a powerful immunosuppressant drug that is administered daily to millions of transplant recipients worldwide. Production of mycophenolic acid is restricted to a very limited number of filamentous fungi, and little is known about its biosynthetic modalities. It is therefore a particular challenge to improve our knowledge of the biosynthesis of this valuable natural compound, as this would contribute to a better understanding of the specialized metabolism of fungi and could also lead to the identification of new fungal producers for the supply of immunosuppressants. Here, we were interested in deciphering the origin and evolution of the fungal mycophenolic acid biosynthetic pathway. Large-scale analyses of fungal genomic resources led us to identify several new species that harbor a gene cluster for mycophenolic acid biosynthesis. Phylogenomic analysis suggests that the mycophenolic acid biosynthetic gene cluster originated early in a common ancestor of the fungal family Aspergillaceae but was repeatedly lost and it is now present in a narrow but diverse set of filamentous fungi. Moreover, a comparison of the inosine 5'-monophosphate dehydrogenase protein sequences that are the target of the mycophenolic acid drug as well as analysis of mycophenolic acid production and susceptibility suggest that all mycophenolic acid fungal producers are resistant to this toxic compound, but that this resistance is likely to be based on different molecular mechanisms. Our study provides new insight into the evolution of the biosynthesis of the important secondary metabolite mycophenolic acid in fungi.

Evolution and diversification of the momilactone biosynthetic gene cluster in the genus Oryza

Plants are master chemists and collectively are able to produce hundreds of thousands of different organic compounds. The genes underlying the biosynthesis of many specialized metabolites are organized in biosynthetic gene clusters (BGCs), which is hypothesized to ensure their faithful coinheritance and to facilitate their coordinated expression. In rice (Oryza sativa), momilactones are diterpenoids that act in plant defence and various organismic interactions. Many of the genes essential for momilactone biosynthesis are grouped in a BGC. We applied comparative genomics of diploid and allotetraploid Oryza species to reconstruct the species-specific architecture, evolutionary trajectory, and sub-functionalisation of the momilactone biosynthetic gene cluster (MBGC) in the Oryza genus. Our data show that the evolution of the MBGC is marked by lineage-specific rearrangements and gene copy number variation, as well as by occasional cluster loss. We identified a distinct cluster architecture in Oryza coarctata, which represents the first instance of an alternative architecture of the MBGC in Oryza and strengthens the idea of a common origin of the cluster in Oryza and the distantly related genus Echinochloa. Our research illustrates the evolutionary and functional dynamics of a biosynthetic gene cluster within a plant genus.

Subfunctionalization and epigenetic regulation of a biosynthetic gene cluster in Solanaceae

Biosynthetic gene clusters (BGCs) are sets of often heterologous genes that are genetically and functionally linked. Among eukaryotes, BGCs are most common in plants and fungi and ensure the coexpression of the different enzymes coordinating the biosynthesis of specialized metabolites. Here, we report the identification of a withanolide BGC in Physalis grisea (ground-cherry), a member of the nightshade family (Solanaceae). A combination of transcriptomic, epigenomic, and metabolic analyses revealed that, following a duplication event, this BGC evolved two tissue-specifically expressed subclusters, containing several pairs of paralogs that contribute to related but distinct biochemical processes; this subfunctionalization is tightly associated with epigenetic features and the local chromatin environment. The two subclusters appear strictly isolated from each other at the structural chromatin level, each forming a highly self-interacting chromatin domain with tissue-dependent levels of condensation. This correlates with gene expression in either above- or below-ground tissue, thus spatially separating the production of different withanolide compounds. By comparative phylogenomics, we show that the withanolide BGC most likely evolved before the diversification of the Solanaceae family and underwent lineage-specific diversifications and losses. The tissue-specific subfunctionalization is common to species of the Physalideae tribe but distinct from other, independent duplication events outside of this clade. In sum, our study reports on an instance of an epigenetically modulated subfunctionalization within a BGC and sheds light on the biosynthesis of withanolides, a highly diverse group of steroidal triterpenoids important in plant defense and amenable to pharmaceutical applications due to their anti-inflammatory, antibiotic, and anticancer properties.

Genomic and cell-specific regulation of benzylisoquinoline alkaloid biosynthesis in opium poppy

Opium poppy is a crop of great commercial value as a source of several opium alkaloids for the pharmaceutical industries including morphine, codeine, thebaine, noscapine, and papaverine. Most enzymes involved in benzylisoquinoline alkaloid (BIA) biosynthesis in opium poppy have been functionally characterized, and opium poppy currently serves as a model system to study BIA metabolism in plants. BIA biosynthesis in opium poppy involves two biosynthetic gene clusters associated respectively with the morphine and noscapine branches. Recent reports have shown that genes in the same cluster are co-expressed, suggesting they might also be co-regulated. However, the transcriptional regulation of opium poppy BIA biosynthesis is not well studied. Opium poppy BIA biosynthesis involves three cell types associated with the phloem system: companion cells, sieve elements, and laticifers. The transcripts and enzymes associated with BIA biosynthesis are distributed across cell types, requiring the translocation of key enzymes and pathway intermediates between cell types. Together, these suggest that the regulation of BIA biosynthesis in opium poppy is multilayered and complex, involving biochemical, genomic, and physiological mechanisms. In this review, we highlight recent advances in genome sequencing and single cell and spatial transcriptomics with a focus on how these efforts can improve our understanding of the genomic and cell-specific regulation of BIA biosynthesis. Such knowledge is vital for opium poppy genetic improvement and metabolic engineering efforts targeting the modulation of alkaloid yield and composition.

The ABC transporter SmABCG1 mediates tanshinones export from the peridermic cells of Salvia miltiorrhiza root

Plants have mechanisms to transport secondary metabolites from where they are biosynthesized to the sites where they function, or to sites such as the vacuole for detoxification. However, current research has mainly focused on metabolite biosynthesis and regulation, and little is known about their transport. Tanshinone, a class diterpenoid with medicinal properties, is biosynthesized in the periderm of Salvia miltiorrhiza roots. Here, we discovered that tanshinone can be transported out of peridermal cells and secreted into the soil environment and that the ABC transporter SmABCG1 is involved in the efflux of tanshinone ⅡA and tanshinone Ⅰ. The SmABCG1 gene is adjacent to the diterpene biosynthesis gene cluster in the S. miltiorrhiza genome. The temporal-spatial expression pattern of SmABCG1 is consistent with tanshinone accumulation profiles. SmABCG1 is located on the plasma membrane and preferentially accumulates in the peridermal cells of S. miltiorrhiza roots. Heterologous expression in Xenopus laevis oocytes demonstrated that SmABCG1 can export tanshinone ⅡA and tanshinone Ⅰ. CRISPR/Cas9-mediated mutagenesis of SmABCG1 in S. miltiorrhiza hairy roots resulted in a significant decrease in tanshinone contents in both hairy roots and the culture medium, whereas overexpression of this gene resulted in increased tanshinone contents. CYP76AH3 transcript levels increased in hairy roots overexpressing SmABCG1 and decreased in knockout lines, suggesting that SmABCG1 may affect the expression of CYP76AH3, indirectly regulating tanshinone biosynthesis. Finally, tanshinone ⅡA showed cytotoxicity to Arabidopsis roots. These findings offer new perspectives on plant diterpenoid transport and provide a new genetic tool for metabolic engineering and synthetic biology research.

Glycyrrhiza, a commonly used medicinal herb: Review of species classification, pharmacology, active ingredient biosynthesis, and synthetic biology

BACKGROUND

Licorice is extensively and globally utilized as a medicinal herb and is one of the traditional Chinese herbal medicines with valuable pharmacological effects. Its therapeutic components primarily reside within its roots and rhizomes, classifying it as a tonifying herb. As more active ingredients in licorice are unearthed and characterized, licorice germplasm resources are gaining more and more recognition. However, due to the excessive exploitation of wild licorice resources, the degrading germplasm reserves fail to meet the requirements of chemical extraction and clinical application.

AIM OF REVIEW

This article presents a comprehensive review of the classification and phylogenetic relationships of species in genus Glycyrrhiza, types of active components and their pharmacological activities, licorice omics, biosynthetic pathways of active compounds in licorice, and metabolic engineering. It aims to offer a unique and comprehensive perspective on Glycyrrhiza, integrating knowledge from diverse fields to offer a comprehensive understanding of this genus. It will serve as a valuable resource and provide a solid foundation for future research and development in the molecular breeding and synthetic biology fields of Glycyrrhiza.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Licorice has an abundance of active constituents, primarily triterpenoids, flavonoids, and polysaccharides. Modern pharmacological research unveiled its multifaceted effects encompassing anti-inflammatory, analgesic, anticancer, antiviral, antioxidant, and hepatoprotective activities. Many resources of Glycyrrhiza species remain largely untapped, and multiomic studies of the Glycyrrhiza lineage are expected to facilitate new discoveries in the fields of medicine and human health. Therefore, strategies for breeding high-yield licorice plants and developing effective biosynthesis methods for bioactive compounds will provide valuable insights into resource conservation and drug development. Metabolic engineering and microorganism-based green production provide alternative strategies to improve the production efficiency of natural products.

Structural diversity, evolutionary origin, and metabolic engineering of plant specialized benzylisoquinoline alkaloids

Covering: up to June 2024Benzylisoquinoline alkaloids (BIAs) represent a diverse class of plant specialized metabolites derived from L-tyrosine, exhibiting significant pharmacological properties such as anti-microbial, anti-spasmodic, anti-cancer, cardiovascular protection, and analgesic effects. The industrial production of valuable BIAs relies on extraction from plants; however, challenges concerning their low concentration and efficiency hinder drug development. Hence, alternative approaches, including biosynthesis and chemoenzymatic synthesis, have been explored. Model species like Papaver somniferum and Coptis japonica have played a key role in unraveling the biosynthetic pathways of BIAs; however, many aspects, particularly modified steps like oxidation and methylation, remain unclear. Critical enzymes, e.g., CYP450s and methyltransferases, play a substantial role in BIA backbone formation and modification, which is essential for understanding the origin and adaptive evolution of these plant specialized metabolites. This review comprehensively analyzes the structural diversity of reported BIAs and their distribution in plant lineages. In addition, the progress in understanding biosynthesis, evolution, and catalytic mechanisms underlying BIA biosynthesis is summarized. Finally, we discuss the progress and challenges in metabolic engineering, providing valuable insights into BIA drug development and the sustainable utilization of BIA-producing plants.

Two ubiquitous aldo-keto reductases in the genus Papaver support a patchwork model for morphine pathway evolution

The evolution of morphinan alkaloid biosynthesis in plants of the genus Papaver includes permutation of several processes including gene duplication, fusion, neofunctionalization, and deletion resulting in the present chemotaxonomy. A critical gene fusion event resulting in the key bifunctional enzyme reticuline epimerase (REPI), which catalyzes the stereochemical inversion of (S)-reticuline, was suggested to precede neofunctionalization of downstream enzymes leading to morphine biosynthesis in opium poppy (Papaver somniferum). The ancestrally related aldo-keto reductases 1,2-dehydroreticuline reductase (DRR), which occurs in some species as a component of REPI, and codeinone reductase (COR) catalyze the second and penultimate steps, respectively, in the pathway converting (S)-reticuline to morphine. Orthologs for each enzyme isolated from the transcriptomes of 12 Papaver species were shown to catalyze their respective reactions in species that capture states of the metabolic pathway prior to key evolutionary events, including the gene fusion event leading to REPI, thus suggesting a patchwork model for pathway evolution. Analysis of the structure and substrate preferences of DRR orthologs in comparison with COR orthologs revealed structure-function relationships underpinning the functional latency of DRR and COR orthologs in the genus Papaver, thus providing insights into the molecular events leading to the evolution of the pathway.

Herbgenomics meets Papaveraceae: a promising -omics perspective on medicinal plant research

Herbal medicines were widely used in ancient and modern societies as remedies for human ailments. Notably, the Papaveraceae family includes well-known species, such as Papaver somniferum and Chelidonium majus, which possess medicinal properties due to their latex content. Latex-bearing plants are a rich source of diverse bioactive compounds, with applications ranging from narcotics to analgesics and relaxants. With the advent of high-throughput technologies and advancements in sequencing tools, an opportunity exists to bridge the knowledge gap between the genetic information of herbs and the regulatory networks underlying their medicinal activities. This emerging discipline, known as herbgenomics, combines genomic information with other -omics studies to unravel the genetic foundations, including essential gene functions and secondary metabolite biosynthesis pathways. Furthermore, exploring the genomes of various medicinal plants enables the utilization of modern genetic manipulation techniques, such as Clustered Regularly-Interspaced Short Palindromic Repeats (CRISPR/Cas9) or RNA interference. This technological revolution has facilitated systematic studies of model herbs, targeted breeding of medicinal plants, the establishment of gene banks and the adoption of synthetic biology approaches. In this article, we provide a comprehensive overview of the recent advances in genomic, transcriptomic, proteomic and metabolomic research on species within the Papaveraceae family. Additionally, it briefly explores the potential applications and key opportunities offered by the -omics perspective in the pharmaceutical industry and the agrobiotechnology field.

The centromere landscapes of four karyotypically diverse Papaver species provide insights into chromosome evolution and speciation

Understanding the roles played by centromeres in chromosome evolution and speciation is complicated by the fact that centromeres comprise large arrays of tandemly repeated satellite DNA, which hinders high-quality assembly. Here, we used long-read sequencing to generate nearly complete genome assemblies for four karyotypically diverse Papaver species, P. setigerum (2n = 44), P. somniferum (2n = 22), P. rhoeas (2n = 14), and P. bracteatum (2n = 14), collectively representing 45 gapless centromeres. We identified four centromere satellite (cenSat) families and experimentally validated two representatives. For the two allopolyploid genomes (P. somniferum and P. setigerum), we characterized the subgenomic distribution of each satellite and identified a "homogenizing" phase of centromere evolution in the aftermath of hybridization. An interspecies comparison of the peri-centromeric regions further revealed extensive centromere-mediated chromosome rearrangements. Taking these results together, we propose a model for studying cenSat competition after hybridization and shed further light on the complex role of the centromere in speciation.

A reference-grade genome of the xerophyte Ammopiptanthus mongolicus sheds light on its evolution history in legumes and drought-tolerance mechanisms

Plants that grow in extreme environments represent unique sources of stress-resistance genes and mechanisms. Ammopiptanthus mongolicus (Leguminosae) is a xerophytic evergreen broadleaf shrub native to semi-arid and desert regions; however, its drought-tolerance mechanisms remain poorly understood. Here, we report the assembly of a reference-grade genome for A. mongolicus, describe its evolutionary history within the legume family, and examine its drought-tolerance mechanisms. The assembled genome is 843.07 Mb in length, with 98.7% of the sequences successfully anchored to the nine chromosomes of A. mongolicus. The genome is predicted to contain 47 611 protein-coding genes, and 70.71% of the genome is composed of repetitive sequences; these are dominated by transposable elements, particularly long-terminal-repeat retrotransposons. Evolutionary analyses revealed two whole-genome duplication (WGD) events at 130 and 58 million years ago (mya) that are shared by the genus Ammopiptanthus and other legumes, but no species-specific WGDs were found within this genus. Ancestral genome reconstruction revealed that the A. mongolicus genome has undergone fewer rearrangements than other genomes in the legume family, confirming its status as a "relict plant". Transcriptomic analyses demonstrated that genes involved in cuticular wax biosynthesis and transport are highly expressed, both under normal conditions and in response to polyethylene glycol-induced dehydration. Significant induction of genes related to ethylene biosynthesis and signaling was also observed in leaves under dehydration stress, suggesting that enhanced ethylene response and formation of thick waxy cuticles are two major mechanisms of drought tolerance in A. mongolicus. Ectopic expression of AmERF2, an ethylene response factor unique to A. mongolicus, can markedly increase the drought tolerance of transgenic Arabidopsis thaliana plants, demonstrating the potential for application of A. mongolicus genes in crop improvement.

Chromosome genome assembly of the Camphora longepaniculata (Gamble) with PacBio and Hi-C sequencing data

Introduction

Camphora longepaniculata, a crucial commercial crop and a fundamental component of traditional Chinese medicine, is renowned for its abundant production of volatile terpenoids. However, the lack of available genomic information has hindered pertinent research efforts in the past.

Methods

To bridge this gap, the present study aimed to use PacBio HiFi, short-read, and highthroughput chromosome conformation capture sequencing to construct a chromosome-level assembly of the C. longepaniculata genome.

Results and discussion

With twelve chromosomes accounting for 99.82% (766.69 Mb) of the final genome assembly, which covered 768.10 Mb, it was very complete. Remarkably, the assembly's contig and scaffold N50 values are exceptional as well-41.12 and 63.78 Mb, respectively-highlighting its excellent quality and intact structure. Furthermore, a total of 39,173 protein-coding genes were predicted, with 38,766 (98.96%) of them being functionally annotated. The completeness of the genome was confirmed by the Benchmarking Universal Single-Copy Ortholog evaluation, which revealed 99.01% of highly conserved plant genes. As the first comprehensive assembly of the C. longepaniculata genome, it provides a crucial starting point for deciphering the complex pathways involved in terpenoid production. Furthermore, this excellent genome serves as a vital resource for upcoming research on the breeding and genetics of C. longepaniculata.

Genomic incongruence accompanies the evolution of flower symmetry in Eudicots: a case study in the poppy family (Papaveraceae, Ranunculales)

Reconstructing evolutionary trajectories and transitions that have shaped floral diversity relies heavily on the phylogenetic framework on which traits are modelled. In this study, we focus on the angiosperm order Ranunculales, sister to all other eudicots, to unravel higher-level relationships, especially those tied to evolutionary transitions in flower symmetry within the family Papaveraceae. This family presents an astonishing array of floral diversity, with actinomorphic, disymmetric (two perpendicular symmetry axes), and zygomorphic flowers. We generated nuclear and plastid datasets using the Angiosperms353 universal probe set for target capture sequencing (of 353 single-copy nuclear ortholog genes), together with publicly available transcriptome and plastome data mined from open-access online repositories. We relied on the fossil record of the order Ranunculales to date our phylogenies and to establish a timeline of events. Our phylogenomic workflow shows that nuclear-plastid incongruence accompanies topological uncertainties in Ranunculales. A cocktail of incomplete lineage sorting, post-hybridization introgression, and extinction following rapid speciation most likely explain the observed knots in the topology. These knots coincide with major floral symmetry transitions and thus obscure the order of evolutionary events.

Pindel-TD: A Tandem Duplication Detector Based on A Pattern Growth Approach

Tandem duplication (TD) is a major type of structural variations (SVs) that plays an important role in novel gene formation and human diseases. However, TDs are often missed or incorrectly classified as insertions by most modern SV detection methods due to the lack of specialized operation on TD-related mutational signals. Herein, we developed a TD detection module for the Pindel tool, referred to as Pindel-TD, based on a TD-specific pattern growth approach. Pindel-TD is capable of detecting TDs with a wide size range at single nucleotide resolution. Using simulated and real read data from HG002, we demonstrated that Pindel-TD outperforms other leading methods in terms of precision, recall, F1-score, and robustness. Furthermore, by applying Pindel-TD to data generated from the K562 cancer cell line, we identified a TD located at the seventh exon of SAGE1, providing an explanation for its high expression. Pindel-TD is available for non-commercial use at https://github.com/xjtu-omics/pindel.

Lineage-Specific CYP80 Expansion and Benzylisoquinoline Alkaloid Diversity in Early-Diverging Eudicots

Menispermaceae species, as early-diverging eudicots, can synthesize valuable benzylisoquinoline alkaloids (BIAs) like bisbenzylisoquinoline alkaloids (bisBIAs) and sinomenines with a wide range of structural diversity. However, the evolutionary mechanisms responsible for their chemo-diversity are not well understood. Here, a chromosome-level genome assembly of Menispermum dauricum is presented and demonstrated the occurrence of two whole genome duplication (WGD) events that are shared by Ranunculales and specific to Menispermum, providing a model for understanding chromosomal evolution in early-diverging eudicots. The biosynthetic pathway for diverse BIAs in M. dauricum is reconstructed by analyzing the transcriptome and metabolome. Additionally, five catalytic enzymes - one norcoclaurine synthase (NCS) and four cytochrome P450 monooxygenases (CYP450s) - from M. dauricum are responsible for the formation of the skeleton, hydroxylated modification, and C-O/C-C phenol coupling of BIAs. Notably, a novel leaf-specific MdCYP80G10 enzyme that catalyzes C2'-C4a phenol coupling of (S)-reticuline into sinoacutine, the enantiomer of morphinan compounds, with predictable stereospecificity is discovered. Moreover, it is found that Menispermum-specific CYP80 gene expansion, as well as tissue-specific expression, has driven BIA diversity in Menispermaceae as compared to other Ranunculales species. This study sheds light on WGD occurrences in early-diverging eudicots and the evolution of diverse BIA biosynthesis.

NextDenovo: an efficient error correction and accurate assembly tool for noisy long reads

Long-read sequencing data, particularly those derived from the Oxford Nanopore sequencing platform, tend to exhibit high error rates. Here, we present NextDenovo, an efficient error correction and assembly tool for noisy long reads, which achieves a high level of accuracy in genome assembly. We apply NextDenovo to assemble 35 diverse human genomes from around the world using Nanopore long-read data. These genomes allow us to identify the landscape of segmental duplication and gene copy number variation in modern human populations. The use of NextDenovo should pave the way for population-scale long-read assembly using Nanopore long-read data.

A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum

BACKGROUND

Rhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism.

RESULTS

We assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening.

CONCLUSION

Overall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.

Genomic data and ecological niche modeling reveal an unusually slow rate of molecular evolution in the Cretaceous Eupteleaceae

Living fossils are evidence of long-term sustained ecological success. However, whether living fossils have little molecular changes remains poorly known, particularly in plants. Here, we have introduced a novel method that integrates phylogenomic, comparative genomic, and ecological niche modeling analyses to investigate the rate of molecular evolution of Eupteleaceae, a Cretaceous relict angiosperm family endemic to East Asia. We assembled a high-quality chromosome-level nuclear genome, and the chloroplast and mitochondrial genomes of a member of Eupteleaceae (Euptelea pleiosperma). Our results show that Eupteleaceae is most basal in Ranunculales, the earliest-diverging order in eudicots, and shares an ancient whole-genome duplication event with the other Ranunculales. We document that Eupteleaceae has the slowest rate of molecular changes in the observed angiosperms. The unusually low rate of molecular evolution of Eupteleaceae across all three independent inherited genomes and genes within each of the three genomes is in association with its conserved genome architecture, ancestral woody habit, and conserved niche requirements. Our findings reveal the evolution and adaptation of living fossil plants through large-scale environmental change and also provide new insights into early eudicot diversification.

The RanOmics group, Becker A, Bachelier JB, Carrive L, Conde e Silva N, Damerval C, Del Rio C, Deveaux Y, Di Stilio VS, Gong Y, Jabbour F, Kramer EM, Nadot S, Pabón-Mora N, Wang W. 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.

Molecular insights into self-incompatibility systems: From evolution to breeding

Plants have evolved diverse self-incompatibility (SI) systems for outcrossing. Since Darwin's time, considerable progress has been made toward elucidating this unrivaled reproductive innovation. Recent advances in interdisciplinary studies and applications of biotechnology have given rise to major breakthroughs in understanding the molecular pathways that lead to SI, particularly the strikingly different SI mechanisms that operate in Solanaceae, Papaveraceae, Brassicaceae, and Primulaceae. These best-understood SI systems, together with discoveries in other "nonmodel" SI taxa such as Poaceae, suggest a complex evolutionary trajectory of SI, with multiple independent origins and frequent and irreversible losses. Extensive exploration of self-/nonself-discrimination signaling cascades has revealed a comprehensive catalog of male and female identity genes and modifier factors that control SI. These findings also enable the characterization, validation, and manipulation of SI-related factors for crop improvement, helping to address the challenges associated with development of inbred lines. Here, we review current knowledge about the evolution of SI systems, summarize key achievements in the molecular basis of pollen‒pistil interactions, discuss potential prospects for breeding of SI crops, and raise several unresolved questions that require further investigation.

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.

Identification and Validation of Magnolol Biosynthesis Genes in Magnolia officinalis

Bacterial infections pose a significant risk to human health. Magnolol, derived from Magnolia officinalis, exhibits potent antibacterial properties. Synthetic biology offers a promising approach to manufacture such natural compounds. However, the plant-based biosynthesis of magnolol remains obscure, and the lack of identification of critical genes hampers its synthetic production. In this study, we have proposed a one-step conversion of magnolol from chavicol using laccase. After leveraging 20 transcriptomes from diverse parts of M. officinalis, transcripts were assembled, enriching genome annotation. Upon integrating this dataset with current genomic information, we could identify 30 laccase enzymes. From two potential gene clusters associated with magnolol production, highly expressed genes were subjected to functional analysis. In vitro experiments confirmed MoLAC14 as a pivotal enzyme in magnolol synthesis. Improvements in the thermal stability of MoLAC14 were achieved through selective mutations, where E345P, G377P, H347F, E346C, and E346F notably enhanced stability. By conducting alanine scanning, the essential residues in MoLAC14 were identified, and the L532A mutation further boosted magnolol production to an unprecedented level of 148.83 mg/L. Our findings not only elucidated the key enzymes for chavicol to magnolol conversion, but also laid the groundwork for synthetic biology-driven magnolol production, thereby providing valuable insights into M. officinalis biology and comparative plant science.

Intergrative metabolomic and transcriptomic analyses reveal the potential regulatory mechanism of unique dihydroxy fatty acid biosynthesis in the seeds of an industrial oilseed crop Orychophragmus violaceus

BACKGROUND

Orychophragmus violaceus is a potentially important industrial oilseed crop due to the two 24-carbon dihydroxy fatty acids (diOH-FA) that was newly identified from its seed oil via a 'discontinuous elongation' process. Although many research efforts have focused on the diOH-FA biosynthesis mechanism and identified the potential co-expressed diacylglycerol acyltranferase (DGAT) gene associated with triacylglycerol (TAG)-polyestolides biosynthesis, the dynamics of metabolic changes during seed development of O. violaceus as well as its associated regulatory network changes are poorly understood.

RESULTS

In this study, by combining metabolome and transcriptome analysis, we identified that 1,003 metabolites and 22,479 genes were active across four stages of seed development, which were further divided into three main clusters based on the patterns of metabolite accumulation and/or gene expression. Among which, cluster2 was mostly related to diOH-FA biosynthesis pathway. We thus further constructed transcription factor (TF)-structural genes regulatory map for the genes associated with the flavonoids, fatty acids and diOH-FA biosynthesis pathway in this cluster. In particular, several TF families such as bHLH, B3, HD-ZIP, MYB were found to potentially regulate the metabolism associated with the diOH-FA pathway. Among which, multiple candidate TFs with promising potential for increasing the diOH-FA content were identified, and we further traced the evolutionary history of these key genes among species of Brassicaceae.

CONCLUSION

Taken together, our study provides new insight into the gene resources and potential relevant regulatory mechanisms of diOH-FA biosynthesis uniquely in seeds of O. violaceus, which will help to promote the downstream breeding efforts of this potential oilseed crop and advance the bio-lubricant industry.

Strategies on biosynthesis and production of bioactive compounds in medicinal plants

Medicinal plants are a valuable source of essential medicines and herbal products for healthcare and disease therapy. Compared with chemical synthesis and extraction, the biosynthesis of natural products is a very promising alternative for the successful conservation of medicinal plants, and its rapid development will greatly facilitate the conservation and sustainable utilization of medicinal plants. Here, we summarize the advances in strategies and methods concerning the biosynthesis and production of natural products of medicinal plants. The strategies and methods mainly include genetic engineering, plant cell culture engineering, metabolic engineering, and synthetic biology based on multiple "OMICS" technologies, with paradigms for the biosynthesis of terpenoids and alkaloids. We also highlight the biosynthetic approaches and discuss progress in the production of some valuable natural products, exemplifying compounds such as vindoline (alkaloid), artemisinin and paclitaxel (terpenoids), to illustrate the power of biotechnology in medicinal plants.

Structural Characterization and Molecular Dynamics Study of the REPI Fusion Protein from Papaver somniferum L

REPI is a pivotal point enzyme in plant benzylisoquinoline alkaloid metabolism as it promotes the evolution of the biosynthetic branch of morphinan alkaloids. Experimental studies of its activity led to the identification of two modules (DRS and DRR) that catalyze two sequential steps of the epimerization of (S)- to (R)-reticuline. Recently, special attention has been paid to its genetic characterization and evolutionary history, but no structural analyses of the REPI protein have been conducted to date. We present here a computational structural characterization of REPI with heme and NADP cofactors in the apo state and in three complexes with substrate (S)-reticuline in DRS and intermediate 1,2-dehydroreticuline in DRS and in DRR. Since no experimental structure exists for REPI, we used its AlphaFold model as a scaffold to build up these four systems, which were submitted to all-atom molecular dynamics (MD) simulations. A comparison of MD results for the four systems revealed key dynamic changes associated with cofactor and ligand binding and provided a dynamic picture of the evolution of their structures and interactions. We also explored the possible dynamic occurrence of tunnels and electrostatic highways potentially involved in alternative mechanisms for channeling the intermediate from DRS to DRR.

Genome-wide expansion and reorganization during grass evolution: from 30 Mb chromosomes in rice and Brachypodium to 550 Mb in Avena

BACKGROUND

The BOP (Bambusoideae, Oryzoideae, and Pooideae) clade of the Poaceae has a common ancestor, with similarities to the genomes of rice, Oryza sativa (2n = 24; genome size 389 Mb) and Brachypodium, Brachypodium distachyon (2n = 10; 271 Mb). We exploit chromosome-scale genome assemblies to show the nature of genomic expansion, structural variation, and chromosomal rearrangements from rice and Brachypodium, to diploids in the tribe Aveneae (e.g., Avena longiglumis, 2n = 2x = 14; 3,961 Mb assembled to 3,850 Mb in chromosomes).

RESULTS

Most of the Avena chromosome arms show relatively uniform expansion over the 10-fold to 15-fold genome-size increase. Apart from non-coding sequence diversification and accumulation around the centromeres, blocks of genes are not interspersed with blocks of repeats, even in subterminal regions. As in the tribe Triticeae, blocks of conserved synteny are seen between the analyzed species with chromosome fusion, fission, and nesting (insertion) events showing deep evolutionary conservation of chromosome structure during genomic expansion. Unexpectedly, the terminal gene-rich chromosomal segments (representing about 50 Mb) show translocations between chromosomes during speciation, with homogenization of genome-specific repetitive elements within the tribe Aveneae. Newly-formed intergenomic translocations of similar extent are found in the hexaploid A. sativa.

CONCLUSIONS

The study provides insight into evolutionary mechanisms and speciation in the BOP clade, which is valuable for measurement of biodiversity, development of a clade-wide pangenome, and exploitation of genomic diversity through breeding programs in Poaceae.

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.

Subcellular compartmentalization in the biosynthesis and engineering of plant natural products

Plant natural products (PNPs) are specialized metabolites with diverse bioactivities. They are extensively used in the pharmaceutical, cosmeceutical and food industries. PNPs are synthesized in plant cells by enzymes that are distributed in different subcellular compartments with unique microenvironments, such as ions, co-factors and substrates. Plant metabolic engineering is an emerging and promising approach for the sustainable production of PNPs, for which the knowledge of the subcellular compartmentalization of their biosynthesis is instrumental. In this review we describe the state of the art on the role of subcellular compartments in the biosynthesis of major types of PNPs, including terpenoids, phenylpropanoids, alkaloids and glucosinolates, and highlight the efforts to target biosynthetic pathways to subcellular compartments in plants. In addition, we will discuss the challenges and strategies in the field of plant synthetic biology and subcellular engineering. We expect that newly developed methods and tools, together with the knowledge gained from the microbial chassis, will greatly advance plant metabolic engineering.

Subgenome phasing for complex allopolyploidy: case-based benchmarking and recommendations

Accurate subgenome phasing is crucial for understanding the origin, evolution and adaptive potential of polyploid genomes. SubPhaser and WGDI software are two common methodologies for subgenome phasing in allopolyploids, particularly in scenarios lacking known diploid progenitors. Triggered by a recent debate over the subgenomic origins of the cultivated octoploid strawberry, we examined four well-documented complex allopolyploidy cases as benchmarks, to evaluate and compare the accuracy of the two software. Our analysis demonstrates that the subgenomic structure phased by both software is in line with prior research, effectively tracing complex allopolyploid evolutionary trajectories despite the limitations of each software. Furthermore, using these validated methodologies, we revisited the controversial issue regarding the progenitors of the octoploid strawberry. The results of both methodologies reaffirm Fragaria vesca and Fragaria iinumae as progenitors of the octoploid strawberry. Finally, we propose recommendations for enhancing the accuracy of subgenome phasing in future studies, recognizing the potential of integrated tools for advanced complex allopolyploidy research and offering a new roadmap for robust subgenome-based phylogenetic analysis.

Characterization of the horse chestnut genome reveals the evolution of aescin and aesculin biosynthesis

Horse chestnut (Aesculus chinensis) is an important medicinal tree that contains various bioactive compounds, such as aescin, barrigenol-type triterpenoid saponins (BAT), and aesculin, a glycosylated coumarin. Herein, we report a 470.02 Mb genome assembly and characterize an Aesculus-specific whole-genome duplication event, which leads to the formation and duplication of two triterpenoid biosynthesis-related gene clusters (BGCs). We also show that AcOCS6, AcCYP716A278, AcCYP716A275, and AcCSL1 genes within these two BGCs along with a seed-specific expressed AcBAHD6 are responsible for the formation of aescin. Furthermore, we identify seven Aesculus-originated coumarin glycoside biosynthetic genes and achieve the de novo synthesis of aesculin in E. coli. Collinearity analysis shows that the collinear BGC segments can be traced back to early-diverging angiosperms, and the essential gene-encoding enzymes necessary for BAT biosynthesis are recruited before the splitting of Aesculus, Acer, and Xanthoceras. These findings provide insight on the evolution of gene clusters associated with medicinal tree metabolites.

Unlocking secrets of nature's chemists: Potential of CRISPR/Cas-based tools in plant metabolic engineering for customized nutraceutical and medicinal profiles

Plant species have evolved diverse metabolic pathways to effectively respond to internal and external signals throughout their life cycle, allowing adaptation to their sessile and phototropic nature. These pathways selectively activate specific metabolic processes, producing plant secondary metabolites (PSMs) governed by genetic and environmental factors. Humans have utilized PSM-enriched plant sources for millennia in medicine and nutraceuticals. Recent technological advances have significantly contributed to discovering metabolic pathways and related genes involved in the biosynthesis of specific PSM in different food crops and medicinal plants. Consequently, there is a growing demand for plant materials rich in nutrients and bioactive compounds, marketed as "superfoods". To meet the industrial demand for superfoods and therapeutic PSMs, modern methods such as system biology, omics, synthetic biology, and genome editing (GE) play a crucial role in identifying the molecular players, limiting steps, and regulatory circuitry involved in PSM production. Among these methods, clustered regularly interspaced short palindromic repeats-CRISPR associated protein (CRISPR/Cas) is the most widely used system for plant GE due to its simple design, flexibility, precision, and multiplexing capabilities. Utilizing the CRISPR-based toolbox for metabolic engineering (ME) offers an ideal solution for developing plants with tailored preventive (nutraceuticals) and curative (therapeutic) metabolic profiles in an ecofriendly way. This review discusses recent advances in understanding the multifactorial regulation of metabolic pathways, the application of CRISPR-based tools for plant ME, and the potential research areas for enhancing plant metabolic profiles.

Functional divergence of CYP76AKs shapes the chemodiversity of abietane-type diterpenoids in genus Salvia

The genus Salvia L. (Lamiaceae) comprises myriad distinct medicinal herbs, with terpenoids as one of their major active chemical groups. Abietane-type diterpenoids (ATDs), such as tanshinones and carnosic acids, are specific to Salvia and exhibit taxonomic chemical diversity among lineages. To elucidate how ATD chemical diversity evolved, we carried out large-scale metabolic and phylogenetic analyses of 71 Salvia species, combined with enzyme function, ancestral sequence and chemical trait reconstruction, and comparative genomics experiments. This integrated approach showed that the lineage-wide ATD diversities in Salvia were induced by differences in the oxidation of the terpenoid skeleton at C-20, which was caused by the functional divergence of the cytochrome P450 subfamily CYP76AK. These findings present a unique pattern of chemical diversity in plants that was shaped by the loss of enzyme activity and associated catalytic pathways.

Revealing evolution of tropane alkaloid biosynthesis by analyzing two genomes in the Solanaceae family

Tropane alkaloids (TAs) are widely distributed in the Solanaceae, while some important medicinal tropane alkaloids (mTAs), such as hyoscyamine and scopolamine, are restricted to certain species/tribes in this family. Little is known about the genomic basis and evolution of TAs biosynthesis and specialization in the Solanaceae. Here, we present chromosome-level genomes of two representative mTAs-producing species: Atropa belladonna and Datura stramonium. Our results reveal that the two species employ a conserved biosynthetic pathway to produce mTAs despite being distantly related within the nightshade family. A conserved gene cluster combined with gene duplication underlies the wide distribution of TAs in this family. We also provide evidence that branching genes leading to mTAs likely have evolved in early ancestral Solanaceae species but have been lost in most of the lineages, with A. belladonna and D. stramonium being exceptions. Furthermore, we identify a cytochrome P450 that modifies hyoscyamine into norhyoscyamine. Our results provide a genomic basis for evolutionary insights into the biosynthesis of TAs in the Solanaceae and will be useful for biotechnological production of mTAs via synthetic biology approaches.

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.

Uncovering a miltiradiene biosynthetic gene cluster in the Lamiaceae reveals a dynamic evolutionary trajectory

The spatial organization of genes within plant genomes can drive evolution of specialized metabolic pathways. Terpenoids are important specialized metabolites in plants with diverse adaptive functions that enable environmental interactions. Here, we report the genome assemblies of Prunella vulgaris, Plectranthus barbatus, and Leonotis leonurus. We investigate the origin and subsequent evolution of a diterpenoid biosynthetic gene cluster (BGC) together with other seven species within the Lamiaceae (mint) family. Based on core genes found in the BGCs of all species examined across the Lamiaceae, we predict a simplified version of this cluster evolved in an early Lamiaceae ancestor. The current composition of the extant BGCs highlights the dynamic nature of its evolution. We elucidate the terpene backbones generated by the Callicarpa americana BGC enzymes, including miltiradiene and the terpene (+)-kaurene, and show oxidization activities of BGC cytochrome P450s. Our work reveals the fluid nature of BGC assembly and the importance of genome structure in contributing to the origin of metabolites.

Emerging mechanistic insights into the regulation of specialized metabolism in plants

Plants biosynthesize a broad range of natural products through specialized and species-specific metabolic pathways that are fuelled by core metabolism, together forming a metabolic network. Specialized metabolites have important roles in development and adaptation to external cues, and they also have invaluable pharmacological properties. A growing body of evidence has highlighted the impact of translational, transcriptional, epigenetic and chromatin-based regulation and evolution of specialized metabolism genes and metabolic networks. Here we review the forefront of this research field and extrapolate to medicinal plants that synthetize rare molecules. We also discuss how this new knowledge could help in improving strategies to produce useful plant-derived pharmaceuticals.

The genome of the king protea, Protea cynaroides

The king protea (Protea cynaroides), an early-diverging eudicot, is the most iconic species from the Megadiverse Cape Floristic Region, and the national flower of South Africa. Perhaps best known for its iconic flower head, Protea is a key genus for the South African horticulture industry and cut-flower market. Ecologically, the genus and the family Proteaceae are important models for radiation and adaptation, particularly to soils with limited phosphorus bio-availability. Here, we present a high-quality chromosome-scale assembly of the P. cynaroides genome as the first representative of the fynbos biome. We reveal an ancestral whole-genome duplication event that occurred in the Proteaceae around the late Cretaceous that preceded the divergence of all crown groups within the family and its extant diversity in all Southern continents. The relatively stable genome structure of P. cynaroides is invaluable for comparative studies and for unveiling paleopolyploidy in other groups, such as the distantly related sister group Ranunculales. Comparative genomics in sequenced genomes of the Proteales shows loss of key arbuscular mycorrhizal symbiosis genes likely ancestral to the family, and possibly the order. The P. cynaroides genome empowers new research in plant diversification, horticulture and adaptation, particularly to nutrient-poor soils.

Haplotype-resolved Chinese male genome assembly based on high-fidelity sequencing

The advantages of both the length and accuracy of high-fidelity (HiFi) reads enable chromosome-scale haplotype-resolved genome assembly. In this study, we sequenced a cell line named HJ, established from a Chinese Han male individual by using HiFi and Hi-C. We assembled two high-quality haplotypes of the HJ genome (haplotype 1 (H1): 3.1 Gb, haplotype 2 (H2): 2.9 Gb). The continuity (H1: contig N50 = 28.2 Mb, H2: contig N50 = 25.9 Mb) and completeness (BUSCO: H1 = 94.9%, H2 = 93.5%) are substantially better than those of other Chinese genomes, for example, HX1, NH1.0, and YH2.0. By comparing HJ genome with GRCh38, we reported the mutation landscape of HJ and found that 176 and 213 N-gaps were filled in H1 and H2, respectively. In addition, we detected 12.9 Mb and 13.4 Mb novel sequences containing 246 and 135 protein-coding genes in H1 and H2, respectively. Our results demonstrate the advantages of HiFi reads in haplotype-resolved genome assembly and provide two high-quality haplotypes of a potential Chinese genome as a reference for the Chinese Han population.

Characterization of benzylisoquinoline alkaloid methyltransferases in Liriodendron chinense provides insights into the phylogenic basis of angiosperm alkaloid diversity

Benzylisoquinoline alkaloids (BIAs) are a class of plant secondary metabolites with great pharmacological value. Their biosynthetic pathways have been extensively elucidated in the species from the Ranunculales order, such as poppy and Coptis japonica, in which methylation events play central roles and are directly responsible for BIA chemodiversity. Here, we combined BIA quantitative profiling and transcriptomic analyses to identify novel BIA methyltransferases (MTs) from Liriodendron chinense, a basal angiosperm plant. We identified an N-methyltransferase (LcNMT1) and two O-methyltransferases (LcOMT1 and LcOMT3), and characterized their biochemical functions in vitro. LcNMT1 methylates (S)-coclaurine to produce mono- and dimethylated products. Mutagenesis experiments revealed that a single-residue alteration is sufficient to change its substrate selectivity. LcOMT1 methylates (S)-norcoclaurine at the C6 site and LcOMT3 methylates (S)-coclaurine at the C7 site, respectively. Two key residues of LcOMT3, A115 and T301, are identified as important contributors to its catalytic activity. Compared with Ranunculales-derived NMTs, Magnoliales-derived NMTs were less abundant and had narrower substrate specificity, indicating that NMT expansion has contributed substantially to BIA chemodiversity in angiosperms, particularly in Ranunculales species. In summary, we not only characterized three novel enzymes that could be useful in the biosynthetic production of valuable BIAs but also shed light on the molecular origin of BIAs during angiosperm evolution.

Natural products of medicinal plants: biosynthesis and bioengineering in post-genomic era

Globally, medicinal plant natural products (PNPs) are a major source of substances used in traditional and modern medicine. As we human race face the tremendous public health challenge posed by emerging infectious diseases, antibiotic resistance and surging drug prices etc., harnessing the healing power of medicinal plants gifted from mother nature is more urgent than ever in helping us survive future challenge in a sustainable way. PNP research efforts in the pre-genomic era focus on discovering bioactive molecules with pharmaceutical activities, and identifying individual genes responsible for biosynthesis. Critically, systemic biological, multi- and inter-disciplinary approaches integrating and interrogating all accessible data from genomics, metabolomics, structural biology, and chemical informatics are necessary to accelerate the full characterization of biosynthetic and regulatory circuitry for producing PNPs in medicinal plants. In this review, we attempt to provide a brief update on the current research of PNPs in medicinal plants by focusing on how different state-of-the-art biotechnologies facilitate their discovery, the molecular basis of their biosynthesis, as well as synthetic biology. Finally, we humbly provide a foresight of the research trend for understanding the biology of medicinal plants in the coming decades.

Lateral transfers lead to the birth of momilactone biosynthetic gene clusters in grass

Momilactone A, an important plant labdane-related diterpenoid, functions as a phytoalexin against pathogens and an allelochemical against neighboring plants. The genes involved in the biosynthesis of momilactone A are found in clusters, i.e., momilactone A biosynthetic gene clusters (MABGCs), in the rice and barnyardgrass genomes. In addition, we know little about the origin and evolution of MABGCs. Here, we integrated results from comprehensive phylogeny and comparative genomic analyses of the core genes of MABGC-like clusters and MABGCs in 40 monocot plant genomes, providing convincing evidence for the birth and evolution of MABGCs in grass species. The MABGCs found in the PACMAD clade of the core grass lineage (including Panicoideae and Chloridoideae) originated from a MABGC-like cluster in Triticeae (BOP clade) via lateral gene transfer (LGT) and followed by recruitment of MAS1/2 and CYP76L1 genes. The MABGCs in Oryzoideae originated from PACMAD through another LGT event and lost CYP76L1 afterwards. The Oryza MABGC and another Oryza diterpenoid cluster c2BGC are two distinct clusters, with the latter originating from gene duplication and relocation within Oryzoideae. Further comparison of the expression patterns of the MABGC genes between rice and barnyardgrass in response to pathogen infection and allelopathy provides novel insights into the functional innovation of MABGCs in plants. Our results demonstrate LGT-mediated origination of MABGCs in grass and shed lights into the evolutionary innovation and optimization of plant biosynthetic pathways.