Gene and genome duplications in the evolution of chemodiversity: perspectives from studies of Lamiaceae

Identification of the ERF gene family of Mangifera indica and the defense response of MiERF4 to Xanthomonas campestris pv. mangiferaeindicae

An important regulatory role for ethylene-responsive transcription factors (ERFs) is in plant growth and development, stress response, and hormone signaling. However, AP2/ERF family genes in mango have not been systematically studied. In this study, a total of 113 AP2/ERF family genes were identified from the mango genome and phylogenetically classified into five subfamilies: AP2 (28 genes), DREB (42 genes), ERF (33 genes), RAV (6 genes), and Soloist (4 genes). Of these, the ERF family, in conjunction with Arabidopsis and rice, forms a phylogenetic tree divided into seven groups, five of which have MiERF members. Analysis of gene structure and cis-elements showed that each MiERF gene contains only one AP2 structural domain, and that MiERF genes contain a large number of cis-elements associated with hormone signaling and stress response. Collinearity tests revealed a high degree of homology between MiERFs and CsERFs. Tissue-specific and stress-responsive expression profiling revealed that MiERF genes are primarily involved in the regulation of reproductive growth and are differentially and positively expressed in response to external hormones and pathogenic bacteria. Physiological results from a gain-of-function analysis of MiERF4 transiently overexpressed in tobacco and mango showed that transient expression of MiERF4 resulted in decreased colony count and callose deposition, as well as varying degrees of response to hormonal signals such as ETH, JA, and SA. Thus, MiERF4 may be involved in the JA/ETH signaling pathway to enhance plant defense against pathogenic bacteria. This study provides a basis for further research on the function and regulation of MiERF genes and lays a foundation for the selection of disease-resistant genes in mango.

Unlocking plant bioactive pathways: omics data harnessing and machine learning assisting

Plant bioactives hold immense potential in the medicine and food industry. The recent advancements in omics applied in deciphering specialized metabolic pathways underscore the importance of high-quality genome releases and the wealth of data in metabolomics and transcriptomics. While harnessing data, whether integrated or standalone, has proven successful in unveiling plant natural product (PNP) biosynthetic pathways, the democratization of machine learning in biology opens exciting new opportunities for enhancing the exploration of these pathways. This review highlights the recent breakthroughs in disrupting plant-specialized biosynthetic pathways through the utilization of omics data harnessing and machine learning techniques.

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.

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

Multiomics analyses of two Leonurus species illuminate leonurine biosynthesis and its evolution

The Lamiaceae family is renowned for its terpenoid-based medicinal components, but Leonurus, which has traditional medicinal uses, stands out for its alkaloid-rich composition. Leonurine, the principal active compound found in Leonurus, has demonstrated promising effects in reducing blood lipids and treating strokes. However, the biosynthetic pathway of leonurine remains largely unexplored. Here, we present the chromosome-level genome sequence assemblies of Leonurus japonicus, known for its high leonurine production, and Leonurus sibiricus, characterized by very limited leonurine production. By integrating genomics, RNA sequencing, metabolomics, and enzyme activity assay data, we constructed the leonurine biosynthesis pathway and identified the arginine decarboxylase (ADC), uridine diphosphate glucosyltransferase (UGT), and serine carboxypeptidase-like (SCPL) acyltransferase enzymes that catalyze key reactions in this pathway. Further analyses revealed that the UGT-SCPL gene cluster evolved by gene duplication in the ancestor of Leonurus and neofunctionalization of SCPL in L. japonicus, which contributed to the accumulation of leonurine specifically in L. japonicus. Collectively, our comprehensive study illuminates leonurine biosynthesis and its evolution in Leonurus.

Medicinal plants enter the single-cell multi-omics era

Elucidating biosynthetic pathways of plant specialized metabolites is a tricky but essential task for the biotechnological production of plant drugs. In a new report, Li et al. used a single-cell multi-omics approach to provide an integrative view of the architecture and regulation of anticancer alkaloid routes in Madagascar periwinkle.

Integrated evolutionary pattern analyses reveal multiple origins of steroidal saponins in plants

Steroidal saponins are a class of specialized metabolites essential for plant's response to biotic and abiotic stresses. They are also important raw materials for the industrial production of steroid drugs. Steroidal saponins are present in some monocots, such as Dioscorea and Paris, but their distribution, origin, and evolution in plants remain poorly understood. By reconstructing the evolutionary history of the steroidal saponin-associated module (SSAM) in plants, we reveal that the steroidal saponin pathway has its origin in Asparagus and Dioscorea. Through evaluating the distribution and evolutionary pattern of steroidal saponins in angiosperms, we further show that steroidal saponins originated multiple times in angiosperms, and exist in early diverged lineages of certain monocot lineages including Asparagales, Dioscoreales, and Liliales. In these lineages, steroidal saponins are synthesized through the high copy and/or high expression mechanisms of key genes in SSAM. Together with shifts in gene evolutionary rates and amino acid usage, these molecular mechanisms shape the current distribution and diversity of steroidal saponins in plants. Consequently, our results provide new insights into the distribution, diversity and evolutionary history of steroidal saponins in plants, and enhance our understanding of plants' resistance to abiotic and biotic stresses. Additionally, fundamental understanding of the steroidal saponin biosynthesis will facilitate their industrial production and pharmacological applications.

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.

Differential RNA-Seq Analysis Predicts Genes Related to Terpene Tailoring in Caryopteris × clandonensis

Enzymatic terpene functionalization is an essential part of plant secondary metabolite diversity. Within this, multiple terpene-modifying enzymes are required to enable the chemical diversity of volatile compounds essential in plant communication and defense. This work sheds light on the differentially transcribed genes within Caryopteris × clandonensis that are capable of functionalizing cyclic terpene scaffolds, which are the product of terpene cyclase action. The available genomic reference was subjected to further improvements to provide a comprehensive basis, where the number of contigs was minimized. RNA-Seq data of six cultivars, Dark Knight, Grand Bleu, Good as Gold, Hint of Gold, Pink Perfection, and Sunny Blue, were mapped on the reference, and their distinct transcription profile investigated. Within this data resource, we detected interesting variations and additionally genes with high and low transcript abundancies in leaves of Caryopteris × clandonensis related to terpene functionalization. As previously described, different cultivars vary in their modification of monoterpenes, especially limonene, resulting in different limonene-derived molecules. This study focuses on predicting the cytochrome p450 enzymes underlying this varied transcription pattern between investigated samples. Thus, making them a reasonable explanation for terpenoid differences between these plants. Furthermore, these data provide the basis for functional assays and the verification of putative enzyme activities.

Single-cell multi-omics in the medicinal plant Catharanthus roseus

Advances in omics technologies now permit the generation of highly contiguous genome assemblies, detection of transcripts and metabolites at the level of single cells and high-resolution determination of gene regulatory features. Here, using a complementary, multi-omics approach, we interrogated the monoterpene indole alkaloid (MIA) biosynthetic pathway in Catharanthus roseus, a source of leading anticancer drugs. We identified clusters of genes involved in MIA biosynthesis on the eight C. roseus chromosomes and extensive gene duplication of MIA pathway genes. Clustering was not limited to the linear genome, and through chromatin interaction data, MIA pathway genes were present within the same topologically associated domain, permitting the identification of a secologanin transporter. Single-cell RNA-sequencing revealed sequential cell-type-specific partitioning of the leaf MIA biosynthetic pathway that, when coupled with a single-cell metabolomics approach, permitted the identification of a reductase that yields the bis-indole alkaloid anhydrovinblastine. We also revealed cell-type-specific expression in the root MIA pathway.

Comparative genomics reveals the diversification of triterpenoid biosynthesis and origin of ocotillol-type triterpenes in Panax

Gene duplication is assumed to be the major force driving the evolution of metabolite biosynthesis in plants. Freed from functional burdens, duplicated genes can mutate toward novelties until fixed due to selective fitness. However, the extent to which this mechanism has driven the diversification of metabolite biosynthesis remains to be tested. Here we performed comparative genomics analysis and functional characterization to evaluate the impact of gene duplication on the evolution of triterpenoid biosynthesis using Panax species as models. We found that whole-genome duplications (WGDs) occurred independently in Araliaceae and Apiaceae lineages. Comparative genomics revealed the evolutionary trajectories of triterpenoid biosynthesis in plants, which was mainly promoted by WGDs and tandem duplication. Lanosterol synthase (LAS) was likely derived from a tandem duplicate of cycloartenol synthase that predated the emergence of Nymphaeales. Under episodic diversifying selection, the LAS gene duplicates produced by γ whole-genome triplication have given rise to triterpene biosynthesis in core eudicots through neofunctionalization. Moreover, functional characterization revealed that oxidosqualene cyclases (OSCs) responsible for synthesizing dammarane-type triterpenes in Panax species were also capable of producing ocotillol-type triterpenes. Genomic and biochemical evidence suggested that Panax genes encoding the above OSCs originated from the specialization of one OSC gene duplicate produced from a recent WGD shared by Araliaceae (Pg-β). Our results reveal the crucial role of gene duplication in diversification of triterpenoid biosynthesis in plants and provide insight into the origin of ocotillol-type triterpenes in Panax species.

A chromosome-level genome assembly reveals that a bipartite gene cluster formed via an inverted duplication controls monoterpenoid biosynthesis in Schizonepeta tenuifolia

Biosynthetic gene clusters (BGCs) are regions of a genome where genes involved in a biosynthetic pathway are in proximity. The origin and evolution of plant BGCs as well as their role in specialized metabolism remain largely unclear. In this study, we have assembled a chromosome-scale genome of Japanese catnip (Schizonepeta tenuifolia) and discovered a BGC that contains multiple copies of genes involved in four adjacent steps in the biosynthesis of p-menthane monoterpenoids. This BGC has an unprecedented bipartite structure, with mirrored biosynthetic regions separated by 260 kilobases. This bipartite BGC includes identical copies of a gene encoding an old yellow enzyme, a type of flavin-dependent reductase. In vitro assays and virus-induced gene silencing revealed that this gene encodes the missing isopiperitenone reductase. This enzyme evolved from a completely different enzyme family to isopiperitenone reductase from closely related Mentha spp., indicating convergent evolution of this pathway step. Phylogenomic analysis revealed that this bipartite BGC has emerged uniquely in the S. tenuifolia lineage and through insertion of pathway genes into a region rich in monoterpene synthases. The cluster gained its bipartite structure via an inverted duplication. The discovered bipartite BGC for p-menthane biosynthesis in S. tenuifolia has similarities to the recently described duplicated p-menthane biosynthesis gene pairs in the Mentha longifolia genome, providing an example of the convergent evolution of gene order. This work expands our understanding of plant BGCs with respect to both form and evolution, and highlights the power of BGCs for gene discovery in plant biosynthetic pathways.

A chromosome-level genome assembly reveals that tandem-duplicated CYP706V oxidase genes control oridonin biosynthesis in the shoot apex of Isodon rubescens

The ent-kaurenoids (e.g., oridonin and enmein) from the Isodon genus (Lamiaceae) are one class of diterpenoids with rich structural diversity and intriguing pharmaceutical activity. In contrast to the well-established gibberellin pathway, oxidative modifications diversifying the ent-kaurene skeleton in Isodon have remained undetermined for half a century. Here we report a chromosome-level genome assembly of I. rubescens, a well-recognized oridonin producer long favored by Asian people as a traditional herb with antitumor effects. The shoot apex was confirmed to be the actual region actively producing ent-kaurene diterpenoids. Through comparative genomics and phylogenetic analyses, we discovered a cluster of tandem-duplicated CYP706V oxygenase-encoding genes located on an ancient genomic block widely distributed in eudicots, whereas almost exclusively emerged in Isodon plants. In the shoot apex, IrCYP706V2 and IrCYP706V7 oxidized the ent-kaurene core in the initial stage of oridonin biosynthesis. Loss of CYP706Vs in other Lamiaceae plants offered an explanation for the specific kaurenoid production in Isodon plants. Moreover, we found that the Isodon genomes encode multiple diterpenoid synthases that are potentially involved in generating diterpenoid diversity. These findings provided new insights into the evolution of the lineage-specific diterpenoid pathway and laid a foundation for improving production of bioactive ent-kaurene-type diterpenoids by molecular breeding and synthetic biology approaches.

The genomes of medicinal skullcaps reveal the polyphyletic origins of clerodane diterpene biosynthesis in the family Lamiaceae

The presence of anticancer clerodane diterpenoids is a chemotaxonomic marker for the traditional Chinese medicinal plant Scutellaria barbata, although the molecular mechanisms behind clerodane biosynthesis are unknown. Here, we report a high-quality assembly of the 414.98 Mb genome of S. barbata into 13 pseudochromosomes. Using phylogenomic and biochemical data, we mapped the plastidial metabolism of kaurene (gibberellins), abietane, and clerodane diterpenes in three species of the family Lamiaceae (Scutellaria barbata, Scutellaria baicalensis, and Salvia splendens), facilitating the identification of genes involved in the biosynthesis of the clerodanes, kolavenol, and isokolavenol. We show that clerodane biosynthesis evolved through recruitment and neofunctionalization of genes from gibberellin and abietane metabolism. Despite the assumed monophyletic origin of clerodane biosynthesis, which is widespread in species of the Lamiaceae, our data show distinct evolutionary lineages and suggest polyphyletic origins of clerodane biosynthesis in the family Lamiaceae. Our study not only provides significant insights into the evolution of clerodane biosynthetic pathways in the mint family, Lamiaceae, but also will facilitate the production of anticancer clerodanes through future metabolic engineering efforts.

Broad diversity in monoterpene-sesquiterpene balance across wild sunflowers: Implications of leaf and floral volatiles for biotic interactions

PREMISE

As plant lineages diversify across environmental gradients, species are predicted to encounter divergent biotic pressures. This study investigated the evolution of volatile secondary metabolism across species of Helianthus.

METHODS

Leaves and petals of 40 species of wild Helianthus were analyzed via gas chromatography-mass spectrometry to determine volatile secondary metabolite profiles.

RESULTS

Across all species, 500 compounds were identified; 40% were sesquiterpenes, 18% monoterpenes, 3% diterpenes, 4% fatty acid derivatives, and 35% other compounds such as phenolics and small organic molecules. Qualitatively, annuals and species from more arid western climates had leaf compositions with a higher proportion of total monoterpenes, while erect perennials and species from more mesic eastern habitats contained a higher proportion of total sesquiterpenes. Among species, mass-based leaf monoterpene and sesquiterpene abundance were identified as largely orthogonal axes of variation by principal component analysis. Profiles for leaves were not strongly correlated with those of petals.

CONCLUSIONS

Volatile metabolites were highly diverse among wild Helianthus, indicating the value of this genus as a model system and rich genetic resource. The independence of leaf and petal volatile profiles indicates a low level of phenotypic integration between vegetative and reproductive structures, implying vegetative defense and reproductive defense or pollinator attraction functions mediated by terpene profiles in these two organs can evolve without major trade-offs. The major biosynthetic pathways for the major terpenes in wild Helianthus are already well described, providing a road map to deeper inquiry into the drivers of this diversity.

Chromosome-level genome assembly of bunching onion illuminates genome evolution and flavor formation in Allium crops

The Allium genus is cultivated globally as vegetables, condiments, or medicinal plants and is characterized by large genomes and strong pungency. However, the genome evolution and genomic basis underlying their unique flavor formation remain poorly understood. Herein, we report an 11.27-Gb chromosome-scale genome assembly for bunching onion (A. fistulosum). The uneven bursts of long-terminal repeats contribute to diversity in genome constituents, and dispersed duplication events largely account for gene expansion in Allium genomes. The extensive duplication and differentiation of alliinase and lachrymatory factor synthase manifest as important evolutionary events during flavor formation in Allium crops. Furthermore, differential selective preference for flavor-related genes likely lead to the variations in isoalliin content in bunching onions. Moreover, we reveal that China is the origin and domestication center for bunching onions. Our findings provide insights into Allium genome evolution, flavor formation and domestication history and enable future genome-assisted breeding of important traits in these crops.

De novo genome assembly of the medicinal plant Gentiana macrophylla provides insights into the genomic evolution and biosynthesis of iridoids

Gentiana macrophylla is a perennial herb in the Gentianaceae family, whose dried roots are used in traditional Chinese medicine. Here, we assembled a chromosome-level genome of G. macrophylla using a combination of Nanopore, Illumina, and Hi-C scaffolding approaches. The final genome size was ~1.79 Gb (contig N50 = 720.804 kb), and 98.89% of the genome sequences were anchored on 13 pseudochromosomes (scaffold N50 = 122.73 Mb). The genome contained 55,337 protein-coding genes, and 73.47% of the assemblies were repetitive sequences. Genome evolution analysis indicated that G. macrophylla underwent two rounds of whole-genome duplication after the core eudicot γ genome triplication event. We further identified candidate genes related to the biosynthesis of iridoids, and the corresponding gene families mostly expanded in G. macrophylla. In addition, we found that root-specific genes are enriched in pathways involved in defense responses, which may greatly improve the biological adaptability of G. macrophylla. Phylogenomic analyses showed a sister relationship of asterids and rosids, and all Gentianales species formed a monophyletic group. Our study contributes to the understanding of genome evolution and active component biosynthesis in G. macrophylla and provides important genomic resource for the genetic improvement and breeding of G. macrophylla.

Identification and bioinformatic analysis of the CaCesA/Csls family members and the expression of the CaCslD1 in the flower buds of CMS/Rf system in pepper

The cellulose synthase gene superfamily contains cellulose synthase (CesA) and cellulose synthase-like (Csl) gene families, which synthesize cellulose and hemicellulose in plant cell walls and play a crucial role in plant growth and development. However, the CesA/Csl gene family has not been reported in pepper. Therefore, the genome-wide research of the CaCesA/CaCsl gene family was conducted in pepper. In this study, a total of 39 CaCesA/CaCsls genes (10 CesAs genes and 29 Csls genes) were identified in pepper and unevenly distributed on 11 chromosomes. These CaCesA/Csls were divided into seven subfamilies (CesAs, CslAs, CslBs, CslCs, CslDs, CslEs, CslGs), and most of CaCesA/Csls genes are closely related to AtCesA/Csls genes. The cis-acting elements in the promoters of CaCesA/Csls genes are mainly related to hormone response and stress response. There are ten collinear gene pairs between the CesA/Csls gene family of pepper and Arabidopsis, and four fragment duplication gene pairs of the CaCesA/Csls genes were discovered. RNA-seq analysis shows that the majority of CaCesA/Csls are expressed in a variety of plant tissues, indicating that most CaCesA/Csls gene expression patterns are not organ-specific, and CaCslD1/D4 have the highest expression in anthers, followed by petal, ovary, and F9. RNA-seq analysis shows that most CaCesA/Csls are responsive to five hormones (IAA, GA₃, ABA, SA, and MeJA). The tissue-specific expression analysis of the CaCslD1 gene shows that the CaCslD1 gene is expressed specifically in flowers. In the flower buds IV of cytoplasmic male sterility (CMS) and its restoration of fertility (Rf) system, CaCslD1 reach the highest expression respectively. However, the relative expression level of CaCslD1 in the fertile accessions is extremely significantly higher than in the sterile accessions. This study shows an overall understanding of the CaCesA/Csls gene family and provides a new insight for understanding the function of CaCslD1 in pollen development and exploring the fertility restoration of CMS in pepper.

The sage genome provides insight into the evolutionary dynamics of diterpene biosynthesis gene cluster in plants

The widely cultivated medicinal and ornamental plant sage (Salvia officinalis L.) is an evergreen shrub of the Lamiaceae family, native to the Mediterranean. We assembled a high-quality sage genome of 480 Mb on seven chromosomes, and identified a biosynthetic gene cluster (BGC) encoding two pairs of diterpene synthases (diTPSs) that, together with the cytochromes P450 (CYPs) genes located inside and outside the cluster, form two expression cascades responsible for the shoot and root diterpenoids, respectively, thus extending BGC functionality from co-regulation to orchestrating metabolite production in different organs. Phylogenomic analysis indicates that the Salvia clades diverged in the early Miocene. In East Asia, most Salvia species are herbaceous and accumulate diterpenoids in storage roots. Notably, in Chinese sage S. miltiorrhiza, the diterpene BGC has contracted and the shoot cascade has been lost. Our data provide genomic insights of micro-evolution of growth type-associated patterning of specialized metabolite production in plants.

Plant biosynthetic gene clusters in the context of metabolic evolution

Covering: up to 2022Plants produce a wide range of structurally and biosynthetically diverse natural products to interact with their environment. These specialised metabolites typically evolve in limited taxonomic groups presumably in response to specific selective pressures. With the increasing availability of sequencing data, it has become apparent that in many cases the genes encoding biosynthetic enzymes for specialised metabolic pathways are not randomly distributed on the genome. Instead they are physically linked in structures such as arrays, pairs and clusters. The exact function of these clusters is debated. In this review we take a broad view of gene arrangement in plant specialised metabolism, examining types of structures and variation. We discuss the evolution of biosynthetic gene clusters in the wider context of metabolism, populations and epigenetics. Finally, we synthesise our observations to propose a new hypothesis for biosynthetic gene cluster formation in plants.

Identification of Abietane-Type Diterpenoids and Phenolic Acids Biosynthesis Genes in Salvia apiana Jepson Through Full-Length Transcriptomic and Metabolomic Profiling

Salvia apiana (S. apiana) Jepson is a medicinal plant that is frequently used by the Chumash Indians in southern California as a diaphoretic, calmative, diuretic, or antimicrobial agent. Abietane-type diterpenoids (ATDs) and phenolic acids (PAs) are the main bioactive ingredients in S. apiana. However, few studies have looked into the biosynthesis of ATDs and PAs in S. apiana. In this study, using metabolic profiling focused on the ATDs and PAs in the roots and leaves of S. apiana, we found a distinctive metabolic feature with all-around accumulation of ATDs, but absence of salvianolic acid B. To identify the candidate genes involved in these biosynthesis pathways, full-length transcriptome was performed by PacBio single-molecule real-time (SMRT) sequencing. A total of 50 and 40 unigenes were predicted to be involved in ATDs and PAs biosynthesis, respectively. Further transcriptional profile using Illumina HiSeq sequencing showed that the transcriptional variations of these pathways were consistent with the accumulation patterns of corresponding metabolites. A plant kingdom-wide phylogenetic analysis of cytochromes (CYPs) identified two CYP76AK and two CYP76AH subfamily genes that might contribute for the specific ATDs biosynthesis in S. apiana. We also noticed that the clade VII laccase gene family was significantly expanded in Salvia miltiorrhiza compared with that of S. apiana, indicating their involvements in the formation of salvianolic acid B. In conclusion, our results will enable the further understanding of ATDs and PAs biosynthesis in S. apiana and Salvia genus.

High Diversity of Type I Polyketide Genes in Bacidia rubella as Revealed by the Comparative Analysis of 23 Lichen Genomes

Fungi involved in lichen symbioses produce a large array of secondary metabolites that are often diagnostic in the taxonomic delimitation of lichens. The most common lichen secondary metabolites—polyketides—are synthesized by polyketide synthases, particularly by Type I PKS (TI-PKS). Here, we present a comparative genomic analysis of the TI-PKS gene content of 23 lichen-forming fungal genomes from Ascomycota, including the de novo sequenced genome of Bacidia rubella. Firstly, we identify a putative atranorin cluster in B. rubella. Secondly, we provide an overview of TI-PKS gene diversity in lichen-forming fungi, and the most comprehensive Type I PKS phylogeny of lichen-forming fungi to date, including 624 sequences. We reveal a high number of biosynthetic gene clusters and examine their domain composition in the context of previously characterized genes, confirming that PKS genes outnumber known secondary substances. Moreover, two novel groups of reducing PKSs were identified. Although many PKSs remain without functional assignments, our findings highlight that genes from lichen-forming fungi represent an untapped source of novel polyketide compounds.

Unlocking plant metabolic diversity: A (pan)-genomic view

Plants produce a remarkable diversity of structurally and functionally diverse natural chemicals that serve as adaptive compounds throughout their life cycles. However, unlocking this metabolic diversity is significantly impeded by the size, complexity, and abundant repetitive elements of typical plant genomes. As genome sequencing becomes routine, we anticipate that links between metabolic diversity and genetic variation will be strengthened. In addition, an ever-increasing number of plant genomes have revealed that biosynthetic gene clusters are not only a hallmark of microbes and fungi; gene clusters for various classes of compounds have also been found in plants, and many are associated with important agronomic traits. We present recent examples of plant metabolic diversification that have been discovered through the exploration and exploitation of various genomic and pan-genomic data. We also draw attention to the fundamental genomic and pan-genomic basis of plant chemodiversity and discuss challenges and future perspectives for investigating metabolic diversity in the coming pan-genomics era.

Tonoplast and Peroxisome Targeting of γ-tocopherol N-methyltransferase Homologs Involved in the Synthesis of Monoterpene Indole Alkaloids

Many plant species from the Apocynaceae, Loganiaceae and Rubiaceae families evolved a specialized metabolism leading to the synthesis of a broad palette of monoterpene indole alkaloids (MIAs). These compounds are believed to constitute a cornerstone of the plant chemical arsenal but above all several MIAs display pharmacological properties that have been exploited for decades by humans to treat various diseases. It is established that MIAs are produced in planta due to complex biosynthetic pathways engaging a multitude of specialized enzymes but also a complex tissue and subcellular organization. In this context, N-methyltransferases (NMTs) represent an important family of enzymes indispensable for MIA biosynthesis but their characterization has always remained challenging. In particular, little is known about the subcellular localization of NMTs in MIA-producing plants. Here, we performed an extensive analysis on the subcellular localization of NMTs from four distinct medicinal plants but also experimentally validated that two putative NMTs from Catharanthus roseus exhibit NMT activity. Apart from providing unprecedented data regarding the targeting of these enzymes in planta, our results point out an additional layer of complexity to the subcellular organization of the MIA biosynthetic pathway by introducing tonoplast and peroxisome as new actors of the final steps of MIA biosynthesis.

Metabolomics and complementary techniques to investigate the plant phytochemical cosmos

Covering: up to 2021Plants and their associated microbial communities are known to produce millions of metabolites, a majority of which are still not characterized and are speculated to possess novel bioactive properties. In addition to their role in plant physiology, these metabolites are also relevant as existing and next-generation medicine candidates. Elucidation of the plant metabolite diversity is thus valuable for the successful exploitation of natural resources for humankind. Herein, we present a comprehensive review on recent metabolomics approaches to illuminate molecular networks in plants, including chemical isolation and enzymatic production as well as the modern metabolomics approaches such as stable isotope labeling, ultrahigh-resolution mass spectrometry, metabolome imaging (spatial metabolomics), single-cell analysis, cheminformatics, and computational mass spectrometry. Mass spectrometry-based strategies to characterize plant metabolomes through metabolite identification and annotation are described in detail. We also highlight the use of phytochemical genomics to mine genes associated with specialized metabolites' biosynthesis. Understanding the metabolic diversity through biotechnological advances is fundamental to elucidate the functions of the plant-derived specialized metabolome.

Flavonoids and saponins in plant rhizospheres: roles, dynamics, and the potential for agriculture

Plants are in constant interaction with a myriad of soil microorganisms in the rhizosphere, an area of soil in close contact with plant roots. Recent research has highlighted the importance of plant-specialized metabolites (PSMs) in shaping and modulating the rhizosphere microbiota; however, the molecular mechanisms underlying the establishment and function of the microbiota mostly remain unaddressed. Flavonoids and saponins are a group of PSMs whose biosynthetic pathways have largely been revealed. Although these PSMs are abundantly secreted into the rhizosphere and exert various functions, the secretion mechanisms have not been clarified. This review summarizes the roles of flavonoids and saponins in the rhizosphere with a special focus on interactions between plants and the rhizosphere microbiota. Furthermore, this review introduces recent advancements in the dynamics of these metabolites in the rhizosphere and indicates potential applications of PSMs for crop production and discusses perspectives in this emerging research field.

Plant Metabolic Gene Clusters: Evolution, Organization, and Their Applications in Synthetic Biology

Plants are a remarkable source of high-value specialized metabolites having significant physiological and ecological functions. Genes responsible for synthesizing specialized metabolites are often clustered together for a coordinated expression, which is commonly observed in bacteria and filamentous fungi. Similar to prokaryotic gene clustering, plants do have gene clusters encoding enzymes involved in the biosynthesis of specialized metabolites. More than 20 gene clusters involved in the biosynthesis of diverse metabolites have been identified across the plant kingdom. Recent studies demonstrate that gene clusters are evolved through gene duplications and neofunctionalization of primary metabolic pathway genes. Often, these clusters are tightly regulated at nucleosome level. The prevalence of gene clusters related to specialized metabolites offers an attractive possibility of an untapped source of highly useful biomolecules. Accordingly, the identification and functional characterization of novel biosynthetic pathways in plants need to be worked out. In this review, we summarize insights into the evolution of gene clusters and discuss the organization and importance of specific gene clusters in the biosynthesis of specialized metabolites. Regulatory mechanisms which operate in some of the important gene clusters have also been briefly described. Finally, we highlight the importance of gene clusters to develop future metabolic engineering or synthetic biology strategies for the heterologous production of novel metabolites.

Linking a Gene Cluster to Atranorin, a Major Cortical Substance of Lichens, through Genetic Dereplication and Heterologous Expression

The depside and depsidone series compounds of polyketide origin accumulate in the cortical or medullary layers of lichen thalli. Despite the taxonomic and ecological significance of lichen chemistry and its pharmaceutical potentials, there has been no single piece of genetic evidence linking biosynthetic genes to lichen substances. Thus, we systematically analyzed lichen polyketide synthases (PKSs) for categorization and identification of the biosynthetic gene cluster (BGC) involved in depside/depsidone production. Our in-depth analysis of the interspecies PKS diversity in the genus Cladonia and a related Antarctic lichen, Stereocaulon alpinum, identified 45 BGC families, linking lichen PKSs to 15 previously characterized PKSs in nonlichenized fungi. Among these, we identified highly syntenic BGCs found exclusively in lichens producing atranorin (a depside). Heterologous expression of the putative atranorin PKS gene (coined atr1) yielded 4-O-demethylbarbatic acid, found in many lichens as a precursor compound, indicating an intermolecular cross-linking activity of Atr1 for depside formation. Subsequent introductions of tailoring enzymes into the heterologous host yielded atranorin, one of the most common cortical substances of macrolichens. Phylogenetic analysis of fungal PKS revealed that the Atr1 is in a novel PKS clade that included two conserved lichen-specific PKS families likely involved in biosynthesis of depsides and depsidones. Here, we provide a comprehensive catalog of PKS families of the genus Cladonia and functionally characterize a biosynthetic gene cluster from lichens, establishing a cornerstone for studying the genetics and chemical evolution of diverse lichen substances. IMPORTANCE Lichens play significant roles in ecosystem function and comprise about 20% of all known fungi. Polyketide-derived natural products accumulate in the cortical and medullary layers of lichen thalli, some of which play key roles in protection from biotic and abiotic stresses (e.g., herbivore attacks and UV irradiation). To date, however, no single lichen product has been linked to respective biosynthetic genes with genetic evidence. Here, we identified a gene cluster family responsible for biosynthesis of atranorin, a cortical substance found in diverse lichen species, by categorizing lichen polyketide synthase and reconstructing the atranorin biosynthetic pathway in a heterologous host. This study will help elucidate lichen secondary metabolism, harnessing the lichen's chemical diversity, hitherto obscured due to limited genetic information on lichens.

Genome sequencing of four culinary herbs reveals terpenoid genes underlying chemodiversity in the Nepetoideae

Species within the mint family, Lamiaceae, are widely used for their culinary, cultural, and medicinal properties due to production of a wide variety of specialized metabolites, especially terpenoids. To further our understanding of genome diversity in the Lamiaceae and to provide a resource for mining biochemical pathways, we generated high-quality genome assemblies of four economically important culinary herbs, namely, sweet basil (Ocimum basilicum L.), sweet marjoram (Origanum majorana L.), oregano (Origanum vulgare L.), and rosemary (Rosmarinus officinalis L.), and characterized their terpenoid diversity through metabolite profiling and genomic analyses. A total 25 monoterpenes and 11 sesquiterpenes were identified in leaf tissue from the 4 species. Genes encoding enzymes responsible for the biosynthesis of precursors for mono- and sesqui-terpene synthases were identified in all four species. Across all 4 species, a total of 235 terpene synthases were identified, ranging from 27 in O. majorana to 137 in the tetraploid O. basilicum. This study provides valuable resources for further investigation of the genetic basis of chemodiversity in these important culinary herbs.

A New Way of Rice Breeding: Polyploid Rice Breeding

Polyploid rice, first discovered by Japanese scientist Eiiti Nakamori in 1933, has a history of nearly 90 years. In the following years, polyploid rice studies have mainly focused on innovations in breeding theory, induction technology and the creation of new germplasm, the analysis of agronomic traits and nutritional components, the study of gametophyte development and reproduction characteristics, DNA methylation modification and gene expression regulation, distant hybridization and utilization among subspecies, species and genomes. In recent years, PMeS lines and neo-tetraploid rice lines with stable high seed setting rate characteristics have been successively selected, breaking through the bottleneck of low seed setting rate of polyploid rice. Following, a series of theoretical and applied studies on high seed setting rate tetraploid rice were carried out. This has pushed research on polyploid rice to a new stage, opening new prospects for polyploid rice breeding.

A route to de novo domestication of wild allotetraploid rice

Cultivated rice varieties are all diploid, and polyploidization of rice has long been desired because of its advantages in genome buffering, vigorousness, and environmental robustness. However, a workable route remains elusive. Here, we describe a practical strategy, namely de novo domestication of wild allotetraploid rice. By screening allotetraploid wild rice inventory, we identified one genotype of Oryza alta (CCDD), polyploid rice 1 (PPR1), and established two important resources for its de novo domestication: (1) an efficient tissue culture, transformation, and genome editing system and (2) a high-quality genome assembly discriminated into two subgenomes of 12 chromosomes apiece. With these resources, we show that six agronomically important traits could be rapidly improved by editing O. alta homologs of the genes controlling these traits in diploid rice. Our results demonstrate the possibility that de novo domesticated allotetraploid rice can be developed into a new staple cereal to strengthen world food security.

Chromosome-level genome assembly of Ophiorrhiza pumila reveals the evolution of camptothecin biosynthesis

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

Gene Coexpression Network Reveals Insights into the Origin and Evolution of a Theanine-Associated Regulatory Module in Non-Camellia and Camellia Species

Theanine (thea) is one of the most important plant-derived characteristic secondary metabolites and a major healthcare product because of its beneficial biological activities, such as being an antianxiety agent, promoting memory, and lowering blood pressure. Thea mostly accumulates in Camellia plants and is especially rich in Camellia sinensis (tea plant). Although some functional genes (e.g., TS, GOGAT, and GS) attributed to thea accumulation have been separately well explored in tea plants, the evolution of a regulatory module (highly interacting gene group) related to thea metabolism remains to be elaborated. Herein, a thea-associated regulatory module (TARM) was mined by using a comprehensive analysis of a weighted gene coexpression network in Camellia and non-Camellia species. Comparative genomic analysis of 84 green plant species revealed that TARM originated from the ancestor of green plants (algae) and that TARM genes were recruited from different evolutionary nodes with the most gene duplication events at the early stage. Among the TARM genes, two core transcription factors of NAC080 and LBD38 were deduced, which may play a crucial role in regulating the biosynthesis of thea. Our findings provide the first insights into the origin and evolution of TARM and indicate a promising paradigm for identifying vital regulatory genes involved in thea metabolism.

Identifying Genes Involved in alkaloid Biosynthesis in Vinca minor Through Transcriptomics and Gene Co-Expression Analysis

The lesser periwinkle Vinca minor accumulates numerous monoterpene indole alkaloids (MIAs) including the vasodilator vincamine. While the biosynthetic pathway of MIAs has been largely elucidated in other Apocynaceae such as Catharanthus roseus, the counterpart in V. minor remains mostly unknown, especially for reactions leading to MIAs specific to this plant. As a consequence, we generated a comprehensive V. minor transcriptome elaborated from eight distinct samples including roots, old and young leaves exposed to low or high light exposure conditions. This optimized resource exhibits an improved completeness compared to already published ones. Through homology-based searches using C. roseus genes as bait, we predicted candidate genes for all common steps of the MIA pathway as illustrated by the cloning of a tabersonine/vincadifformine 16-O-methyltransferase (Vm16OMT) isoform. The functional validation of this enzyme revealed its capacity of methylating 16-hydroxylated derivatives of tabersonine, vincadifformine and lochnericine with a Km 0.94 ± 0.06 µM for 16-hydroxytabersonine. Furthermore, by combining expression of fusions with yellow fluorescent proteins and interaction assays, we established that Vm16OMT is located in the cytosol and forms homodimers. Finally, a gene co-expression network was performed to identify candidate genes of the missing V. minor biosynthetic steps to guide MIA pathway elucidation.

The scaffold-forming steps of plant alkaloid biosynthesis

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