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

Transcriptome sequencing of medical herb Salvia Rosmarinus (Rosemary) revealed the phenylpropanoid biosynthesis pathway genes and their phylogenetic relationships

BACKGROUND

The Salvia rosmarinus spenn. (rosemary) is considered an economically important ornamental and medicinal plant and is widely utilized in culinary and for treating several diseases. However, the procedure behind synthesizing secondary metabolites-based bioactive compounds at the molecular level in S. rosmarinus is not explored completely.

METHODS AND RESULTS

We performed transcriptomic sequencing of the pooled sample from leaf and stem tissues on the Illumina HiSeqTM X10 platform. The transcriptomics analysis led to the generation of 29,523,608 raw reads, followed by data pre-processing which generated 23,208,592 clean reads, and de novo assembly of S. rosmarinus obtained 166,849 unigenes. Among them, nearly 75.1% of unigenes i.e., 28,757 were interpreted against a non-redundant protein database. The gene ontology-based annotation classified them into 3 main categories and 55 sub-categories, and clusters of orthologous genes annotation categorized them into 23 functional categories. The Kyoto Encyclopedia of Genes and Genomes database-based pathway analysis confirmed the involvement of 13,402 unigenes in 183 biochemical pathways, among these unigenes, 1,186 are involved in the 17 secondary metabolite production pathways. Several key enzymes involved in producing aromatic amino acids and phenylpropanoids were identified from the transcriptome database. Among the identified 48 families of transcription factors from coding unigenes, bHLH, MYB, WRKYs, NAC, C2H2, C3H, and ERF are involved in flavonoids and other secondary metabolites biosynthesis.

CONCLUSION

The phylogenetic analysis revealed the evolutionary relationship between the phenylpropanoid pathway genes of rosemary with other members of Lamiaceae. Our work reveals a new molecular mechanism behind the biosynthesis of phenylpropanoids and their regulation in rosemary plants.

Screening and Genomic Analysis of Alkaloid-Producing Endophytic Fungus Fusarium solani Strain MC503 from Macleaya cordata

The extensive harvesting of Macleaya cordata, as a biomedicinal plant and a wild source of quaternary benzo[c]phenanthridine alkaloids, has led to a rapid decline in its population. An alternative approach to the production of these bioactive compounds, which are known for their diverse pharmacological effects, is needed. Production of these compounds using alkaloid-producing endophytic fungi is a promising potential approach. In this research, we isolated an alkaloid-producing endophytic fungus, strain MC503, from the roots of Macleaya cordata. Genomic analysis was conducted to elucidate its metabolic pathways and identify the potential genes responsible for alkaloid biosynthesis. High-performance liquid chromatography (HPLC) and liquid chromatography-mass spectrometry (LC-MS) analyses revealed the presence and quantified the content of sanguinarine (536.87 μg/L) and chelerythrine (393.31 μg/L) in the fungal fermentation extract. Based on our analysis of the morphological and micromorphological characteristics and the ITS region of the nuclear ribosomal DNA of the alkaloid-producing endophyte, it was identified as Fusarium solani strain MC503. To the best of our knowledge, there is no existing report on Fusarium solani from Macleaya cordata or other medicinal plants that produce sanguinarine and chelerythrine simultaneously. These findings provide valuable insights into the capability of Fusarium solani to carry out isoquinoline alkaloid biosynthesis and lay the foundation for further exploration of its potential applications in pharmaceuticals.

Chromosome-scale reference genome of broccoli (Brassica oleracea var. italica Plenck) provides insights into glucosinolate biosynthesis

Broccoli (Brassica oleracea var. italica Plenck) is an important vegetable crop, as it is rich in health-beneficial glucosinolates (GSLs). However, the genetic basis of the GSL diversity in Brassicaceae remains unclear. Here we report a chromosome-level genome assembly of broccoli generated using PacBio HiFi reads and Hi-C technology. The final genome assembly is 613.79 Mb in size, with a contig N50 of 14.70 Mb. The GSL profile and content analysis of different B. oleracea varieties, combined with a phylogenetic tree analysis, sequence alignment, and the construction of a 3D model of the methylthioalkylmalate synthase 1 (MAM1) protein, revealed that the gene copy number and amino acid sequence variation both contributed to the diversity of GSL biosynthesis in B. oleracea. The overexpression of BoMAM1 (BolI0108790) in broccoli resulted in high accumulation and a high ratio of C4-GSLs, demonstrating that BoMAM1 is the key enzyme in C4-GSL biosynthesis. These results provide valuable insights for future genetic studies and nutritive component applications of Brassica crops.

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.

Phylogenomic analyses sheds new light on the phylogeny and diversification of Corydalis DC. in Himalaya-Hengduan Mountains and adjacent regions

The Himalaya-Hengduan Mountains (HHM), a renowned biodiversity hotspot of the world, harbors the most extensive habitats for alpine plants with extraordinary high levels of endemism. Although the general evolution pattern has been elucidated, the underlying processes driving spectacular radiations in many species-rich groups remain elusive. Corydalis DC. is widely distributed throughout the Northern Hemisphere containing more than 500 species, with high diversity in HHM and adjacent regions. Using 95 plastid genes, 3,258,640 nuclear single nucleotide polymorphisms (SNPs) and eight single-copy nuclear genes (SCNs) generated from genome skimming data, we reconstructed a robust time-calibrated phylogeny of Corydalis comprising more than 100 species that represented all subgenera and most sections. Molecular dating indicated that all main clades of Corydalis began to diverge in the Eocene, with the majority of extant species in HHM emerged from a diversification burst after the middle Miocene. Global pattern of mean divergence times indicated that species distributed in HHM were considerably younger than those in other regions, particularly for the two most species-rich clades (V and VI) of Corydalis. The early divergence and the recent diversification of Corydalis were most likely promoted by the continuous orogenesis and climate change associated with the uplift of the Qinghai-Tibetan Plateau (QTP). Our study demonstrates the effectivity of phylogenomic analyses with genome skimming data on the phylogeny of species-rich taxa, and sheds lights on how the uplift of QTP has triggered the evolutionary radiations of large plant genera in HHM and adjacent regions.

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.

The genome of Stephania japonica provides insights into the biosynthesis of cepharanthine

Stephania japonica is an early-diverging eudicotyledon plant with high levels of cepharanthine, proven to be effective in curing coronavirus infections. Here, we report a high-quality S. japonica genome. The genome size is 688.52 Mb, and 97.37% sequences anchor to 11 chromosomes. The genome comprises 67.46% repetitive sequences and 21,036 genes. It is closely related to two Ranunculaceae species, which diverged from their common ancestor 55.90-71.02 million years ago (Mya) with a whole-genome duplication 85.59-96.75 Mya. We further reconstruct ancestral karyotype of Ranunculales. Several cepharanthine biosynthesis genes are identified and verified by western blot. Two genes (Sja03G0243 and Sja03G0241) exhibit catalytic activity as shown by liquid chromatography-mass spectrometry. Then, cepharanthine biosynthesis genes, transcription factors, and CYP450 family genes are used to construct a comprehensive network. Finally, we construct an early-diverging eudicotyledonous genome resources (EEGR) database. As the first genome of the Menispermaceae family to be released, this study provides rich resources for genomic studies.

Cepharanthine analogs mining and genomes of Stephania accelerate anti-coronavirus drug discovery

Cepharanthine is a secondary metabolite isolated from Stephania. It has been reported that it has anti-conronaviruses activities including severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Here, we assemble three Stephania genomes (S. japonica, S. yunnanensis, and S. cepharantha), propose the cepharanthine biosynthetic pathway, and assess the antiviral potential of compounds involved in the pathway. Among the three genomes, S. japonica has a near telomere-to-telomere assembly with one remaining gap, and S. cepharantha and S. yunnanensis have chromosome-level assemblies. Following by biosynthetic gene mining and metabolomics analysis, we identify seven cepharanthine analogs that have broad-spectrum anti-coronavirus activities, including SARS-CoV-2, Guangxi pangolin-CoV (GX_P2V), swine acute diarrhoea syndrome coronavirus (SADS-CoV), and porcine epidemic diarrhea virus (PEDV). We also show that two other genera, Nelumbo and Thalictrum, can produce cepharanthine analogs, and thus have the potential for antiviral compound discovery. Results generated from this study could accelerate broad-spectrum anti-coronavirus drug discovery.

RNA editing events and expression profiles of mitochondrial protein-coding genes in the endemic and endangered medicinal plant, Corydalis saxicola

Corydalis saxicola, an endangered medicinal plant endemic to karst habitats, is widely used in Traditional Chinese Medicine to treat hepatitis, abdominal pain, bleeding hemorrhoids and other conditions. However, to date, the mitochondrial (mt) genome of C. saxicola has not been reported, which limits our understanding of the genetic and biological mechanisms of C. saxicola. Here, the mt genome of C. saxicola was assembled by combining the Nanopore and Illumina reads. The mt genome of C. saxicola is represented by a circular chromosome which is 587,939 bp in length, with an overall GC content of 46.50%. 40 unique protein-coding genes (PCGs), 22 tRNA genes and three rRNA genes were identified. Codon usage of the PCGs was investigated and 167 simple sequence repeats were identified. Twelve homologous fragments were identified between the mt and ct genomes of C. saxicola, accounting for 1.04% of the entire mt genome. Phylogenetic examination of the mt genomes of C. saxicola and 30 other taxa provided an understanding of their evolutionary relationships. We also predicted 779 RNA editing sites in 40 C. saxicola mt PCGs and successfully validated 506 (65%) of these using PCR amplification and Sanger sequencing. In addition, we transcriptionally profiled 24 core mt PCGs in C. saxicola roots treated with different concentrations of CaCl2, as well as in other organs. These investigations will be useful for effective utilization and molecular breeding, and will also provide a reference for further studies of the genus Corydalis.

Genome sequencing and functional analysis of a multipurpose medicinal herb Tinospora cordifolia (Giloy)

Tinospora cordifolia (Willd.) Hook.f. & Thomson, also known as Giloy, is among the most important medicinal plants that have numerous therapeutic applications in human health due to the production of a diverse array of secondary metabolites. To gain genomic insights into the medicinal properties of T. cordifolia, the genome sequencing was carried out using 10× Genomics linked read and Nanopore long-read technologies. The draft genome assembly of T. cordifolia was comprised of 1.01 Gbp, which is the genome sequenced from the plant family Menispermaceae. We also performed the genome size estimation for T. cordifolia, which was found to be 1.13 Gbp. The deep sequencing of transcriptome from the leaf tissue was also performed. The genome and transcriptome assemblies were used to construct the gene set, resulting in 17,245 coding gene sequences. Further, the phylogenetic position of T. cordifolia was also positioned as basal eudicot by constructing a genome-wide phylogenetic tree using multiple species. Further, a comprehensive comparative evolutionary analysis of gene families contraction/expansion and multiple signatures of adaptive evolution was performed. The genes involved in benzyl iso-quinoline alkaloid, terpenoid, lignin and flavonoid biosynthesis pathways were found with signatures of adaptive evolution. These evolutionary adaptations in genes provide genomic insights into the presence of diverse medicinal properties of this plant. The genes involved in the common symbiosis signalling pathway associated with endosymbiosis (Arbuscular Mycorrhiza) were found to be adaptively evolved. The genes involved in adventitious root formation, peroxisome biogenesis, biosynthesis of phytohormones, and tolerance against abiotic and biotic stresses were also found to be adaptively evolved in T. cordifolia.

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.

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.

(±)-Yanhusuomide A, a pair of ornithine-fused benzylisoquinoline enantiomers from Corydalis yanhusuo

(±)-Yanhusuomide A (1), a novel enantiomeric pair of ornithine-fused benzylisoquinoline, were characterized from the dried tubers of Corydalis yanhusuo, along with a biogenetically related intermediate oblongine (2). Yanhusuomide A features an unprecedented skeleton based on a benzylisoquinoline coupled with an ornithine derivative to form a rare 5,6-dihydro-4H-pyrido[3,4,5-de]quinazoline motif. Plausible biosynthetic pathway of 1 was proposed, and (±)-yanhusuomide A (1) presented potential inhibitory bioactivity against acetylcholinesterase (AChE) with IC₅₀ = 14.07 ± 2.38 μM. The simulation of molecular docking displayed that 1 generated strong interaction with Asp-74 and Trp-86 residues of AChE through attractive charge of the quaternary nitrogen.

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.

Integration of high-throughput omics technologies in medicinal plant research: The new era of natural drug discovery

Medicinal plants are natural sources to unravel novel bioactive compounds to satisfy human pharmacological potentials. The world's demand for herbal medicines is increasing year by year; however, large-scale production of medicinal plants and their derivatives is still limited. The rapid development of modern technology has stimulated multi-omics research in medicinal plants, leading to a series of breakthroughs on key genes, metabolites, enzymes involved in biosynthesis and regulation of active compounds. Here, we summarize the latest research progress on the molecular intricacy of medicinal plants, including the comparison of genomics to demonstrate variation and evolution among species, the application of transcriptomics, proteomics and metabolomics to explore dynamic changes of molecular compounds, and the utilization of potential resources for natural drug discovery. These multi-omics research provide the theoretical basis for environmental adaptation of medicinal plants and allow us to understand the chemical diversity and composition of bioactive compounds. Many medicinal herbs' phytochemical constituents and their potential health benefits are not fully explored. Given their large diversity and global distribution as well as the impacts of growth duration and environmental factors on bioactive phytochemicals in medicinal plants, it is crucial to emphasize the research needs of using multi-omics technologies to address basic and applied problems in medicinal plants to aid in developing new and improved medicinal plant resources and discovering novel medicinal ingredients.

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.

New Insights Into the Backbone Phylogeny and Character Evolution of Corydalis (Papaveraceae) Based on Plastome Data

A robust backbone phylogeny is fundamental for developing a stable classification and is instructive for further research. However, it was still not available for Corydalis DC., a species-rich (> 500 species), ecologically and medically important, but taxonomically notoriously difficult genus. Here, we constructed backbone phylogeny and estimated the divergence of Corydalis based on the plastome data from 39 Corydalis species (32 newly sequenced), which represent ca. 80% of sections and series across this genus. Our phylogenetic analyses recovered six fully supported main clades (I-VI) and provided full support for the majority of lineages within Corydalis. Section Archaeocapnos was unexpectedly turned out to be sister to the rest of the subg. Corydalis s. l. (clades IV-VI), thus treating as a distinct clade (clade III) to render all the main clades monophyletic. Additionally, some unusual plastome structural rearrangements were constantly detected within Corydalis and were proven to be lineage-specific in this study, which, in turn, provided further support to our phylogeny. A segment containing five genes (trnV-UAC-rbcL) in the plastome's LSC region was either normally located downstream of the ndhC gene in clade I species or translocated downstream of the atpH gene in clade II species or translocated to downstream of the trnK-UUU gene in clade III-VI species. The unique large inversion (ca. 50 kb) in the plastome LSC region of clade III species, representing an intermediate stage of the above translocation in clades IV-VI, firmly supported clade III as a distinct and early diverged clade within this large lineage (clades III-VI). Our phylogeny contradicted substantially with the morphology-based taxonomy, rejected the treatment of tuberous species as an independent evolutionary group, and proved that some commonly used diagnostic characters (e.g., root and rhizome) were results of convergent evolution, suggestive of unreliability in Corydalis. We dated the origin of crown Corydalis to the early Eocene (crown age 49.08 Ma) and revealed possible explosive radiation around 25 Ma, coinciding with the drastic uplift of the Qinghai-Tibetan Plateau in Oligocene and Miocene. This study provided the most reliable and robust backbone phylogeny of Corydalis to date and shed some new insights on the evolution of Corydalis.