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

Beyond co-expression: pathway discovery for plant pharmaceuticals

Plant natural products have been an important source of medicinal molecules since ancient times. To gain access to the whole diversity of these molecules for pharmaceutical applications, it is important to understand their biosynthetic origins. Whilst co-expression is a reliable tool for identifying gene candidates, a variety of complementary methods can aid in screening or refining candidate selection. Here, we review recently employed plant biosynthetic pathway discovery approaches, and highlight future directions in the field.

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.

Bioinformatic Identification and Expression Analyses of the MAPK-MAP4K Gene Family Reveal a Putative Functional MAP4K10-MAP3K7/8-MAP2K1/11-MAPK3/6 Cascade in Wheat (Triticum aestivum L.)

The mitogen-activated protein kinase (MAPK) cascades act as crucial signaling modules that regulate plant growth and development, response to biotic/abiotic stresses, and plant immunity. MAP3Ks can be activated through MAP4K phosphorylation in non-plant systems, but this has not been reported in plants to date. Here, we identified a total of 234 putative TaMAPK family members in wheat (Triticum aestivum L.). They included 48 MAPKs, 17 MAP2Ks, 144 MAP3Ks, and 25 MAP4Ks. We conducted systematic analyses of the evolution, domain conservation, interaction networks, and expression profiles of these TaMAPK-TaMAP4K (representing TaMAPK, TaMAP2K, TaMAP3K, and TaMAP4K) kinase family members. The 234 TaMAPK-TaMAP4Ks are distributed on 21 chromosomes and one unknown linkage group (Un). Notably, 25 of these TaMAP4K family members possessed the conserved motifs of MAP4K genes, including glycine-rich motif, invariant lysine (K) motif, HRD motif, DFG motif, and signature motif. TaMAPK3 and 6, and TaMAP4K10/24 were shown to be strongly expressed not only throughout the growth and development stages but also in response to drought or heat stress. The bioinformatics analyses and qRT-PCR results suggested that wheat may activate the MAP4K10-MEKK7-MAP2K11-MAPK6 pathway to increase drought resistance in wheat, and the MAP4K10-MAP3K8-MAP2K1/11-MAPK3 pathway may be involved in plant growth. In general, our work identified members of the MAPK-MAP4K cascade in wheat and profiled their potential roles during their response to abiotic stresses and plant growth based on their expression pattern. The characterized cascades might be good candidates for future crop improvement and molecular breeding.

Chromosome-level genome assembly provides insights into the genome evolution and functional importance of the phenylpropanoid-flavonoid pathway in Thymus mongolicus

BACKGROUND

Thymus mongolicus (family Lamiaceae) is a Thyme subshrub with strong aroma and remarkable environmental adaptability. Limited genomic information limits the use of this plant.

RESULTS

Chromosome-level 605.2 Mb genome of T. mongolicus was generated, with 96.28% anchored to 12 pseudochromosomes. The repetitive sequences were dominant, accounting for 70.98%, and 32,593 protein-coding genes were predicted. Synteny analysis revealed that Lamiaceae species generally underwent two rounds of whole genome duplication; moreover, species-specific genome duplication was identified. A recent LTR retrotransposon burst and tandem duplication might play important roles in the formation of the Thymus genome. Using comparative genomic analysis, phylogenetic tree of seven Lamiaceae species was constructed, which revealed that Thyme plants evolved recently in the family. Under the phylogenetic framework, we performed functional enrichment analysis of the genes on nodes that contained the most gene duplication events (> 50% support) and of relevant significant expanded gene families. These genes were highly associated with environmental adaptation and biosynthesis of secondary metabolites. Combined transcriptome and metabolome analyses revealed that Peroxidases, Hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferases, and 4-coumarate-CoA ligases genes were the essential regulators of the phenylpropanoid-flavonoid pathway. Their catalytic products (e.g., apigenin, naringenin chalcone, and several apigenin-related compounds) might be responsible for the environmental tolerance and aromatic properties of T. mongolicus.

CONCLUSION

This study enhanced the understanding of the genomic evolution of T. mongolicus, enabling further exploration of its unique traits and applications, and contributed to the understanding of Lamiaceae genomics and evolutionary biology.

Genome-wide identification, stress- and hormone-responsive expression characteristics, and regulatory pattern analysis of Scutellaria baicalensis SbSPLs

SQUAMOSA PROMOTER BINDING PROTEIN-LIKEs (SPLs) encode plant-specific transcription factors that regulate plant growth and development, stress response, and metabolite accumulation. However, there is limited information on Scutellaria baicalensis SPLs. In this study, 14 SbSPLs were identified and divided into 8 groups based on phylogenetic relationships. SbSPLs in the same group had similar structures. Abscisic acid-responsive (ABRE) and MYB binding site (MBS) cis-acting elements were found in the promoters of 8 and 6 SbSPLs. Segmental duplications and transposable duplications were the main causes of SbSPL expansion. Expression analysis based on transcriptional profiling showed that SbSPL1, SbSPL10, and SbSPL13 were highly expressed in roots, stems, and flowers, respectively. Expression analysis based on quantitative real-time polymerase chain reaction (RT‒qPCR) showed that most SbSPLs responded to low temperature, drought, abscisic acid (ABA) and salicylic acid (SA), among which the expression levels of SbSPL7/9/10/12 were significantly upregulated in response to abiotic stress. These results indicate that SbSPLs are involved in the growth, development and stress response of S. baicalensis. In addition, 8 Sba-miR156/157 s were identified, and SbSPL1-5 was a potential target of Sba-miR156/157 s. The results of target gene prediction and coexpression analysis together indicated that SbSPLs may be involved in the regulation of L-phenylalanine (L-Phe), lignin and jasmonic acid (JA) biosynthesis. In summary, the identification and characterization of the SbSPL gene family lays the foundation for functional research and provides a reference for improved breeding of S. baicalensis stress resistance and quality traits.

Sequence-Structure Analysis Unlocking the Potential Functional Application of the Local 3D Motifs of Plant-Derived Diterpene Synthases

Plant-derived diterpene synthases (PdiTPSs) play a critical role in the formation of structurally and functionally diverse diterpenoids. However, the specificity or functional-related features of PdiTPSs are not well understood. For a more profound insight, we collected, constructed, and curated 199 functionally characterized PdiTPSs and their corresponding 3D structures. The complex correlations among their sequences, domains, structures, and corresponding products were comprehensively analyzed. Ultimately, our focus narrowed to the geometric arrangement of local structures. We found that local structural alignment can rapidly localize product-specific residues that have been validated by mutagenesis experiments. Based on the 3D motifs derived from the residues around the substrate, we successfully searched diterpene synthases (diTPSs) from the predicted terpene synthases and newly characterized PdiTPSs, suggesting that the identified 3D motifs can serve as distinctive signatures in diTPSs (I and II class). Local structural analysis revealed the PdiTPSs with more conserved amino acid residues show features unique to class I and class II, whereas those with fewer conserved amino acid residues typically exhibit product diversity and specificity. These results provide an attractive method for discovering novel or functionally equivalent enzymes and probing the product specificity in cases where enzyme characterization is limited.

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.

Novel Tri-Segmented Rhabdoviruses: A Data Mining Expedition Unveils the Cryptic Diversity of Cytorhabdoviruses

Cytorhabdoviruses (genus Cytorhabdovirus, family Rhabdoviridae) are plant-infecting viruses with enveloped, bacilliform virions. Established members of the genus Cytorhabdovirus have unsegmented single-stranded negative-sense RNA genomes (ca. 10-16 kb) which encode four to ten proteins. Here, by exploring large publicly available metatranscriptomics datasets, we report the identification and genomic characterization of 93 novel viruses with genetic and evolutionary cues of cytorhabdoviruses. Strikingly, five unprecedented viruses with tri-segmented genomes were also identified. This finding represents the first tri-segmented viruses in the family Rhabdoviridae, and they should be classified in a novel genus within this family for which we suggest the name "Trirhavirus". Interestingly, the nucleocapsid and polymerase were the only typical rhabdoviral proteins encoded by those tri-segmented viruses, whereas in three of them, a protein similar to the emaravirus (family Fimoviridae) silencing suppressor was found, while the other predicted proteins had no matches in any sequence databases. Genetic distance and evolutionary insights suggest that all these novel viruses may represent members of novel species. Phylogenetic analyses, of both novel and previously classified plant rhabdoviruses, provide compelling support for the division of the genus Cytorhabdovirus into three distinct genera. This proposed reclassification not only enhances our understanding of the evolutionary dynamics within this group of plant rhabdoviruses but also illuminates the remarkable genomic diversity they encompass. This study not only represents a significant expansion of the genomics of cytorhabdoviruses that will enable future research on the evolutionary peculiarity of this genus but also shows the plasticity in the rhabdovirus genome organization with the discovery of tri-segmented members with a unique evolutionary trajectory.

Gap-free genome assembly and CYP450 gene family analysis reveal the biosynthesis of anthocyanins in Scutellaria baicalensis

Scutellaria baicalensis Georgi, a member of the Lamiaceae family, is a widely utilized medicinal plant. The flavones extracted from S. baicalensis contribute to numerous health benefits, including anti-inflammatory, antiviral, and anti-tumor activities. However, the incomplete genome assembly hinders biological studies on S. baicalensis. This study presents the first telomere-to-telomere (T2T) gap-free genome assembly of S. baicalensis through the integration of Pacbio HiFi, Nanopore ultra-long and Hi-C technologies. A total of 384.59 Mb of genome size with a contig N50 of 42.44 Mb was obtained, and all sequences were anchored into nine pseudochromosomes without any gap or mismatch. In addition, we analysed the major cyanidin- and delphinidin-based anthocyanins involved in the determination of blue-purple flower using a widely-targeted metabolome approach. Based on the genome-wide identification of Cytochrome P450 (CYP450) gene family, three genes (SbFBH1, 2, and 5) encoding flavonoid 3'-hydroxylases (F3'Hs) and one gene (SbFBH7) encoding flavonoid 3'5'-hydroxylase (F3'5'H) were found to hydroxylate the B-ring of flavonoids. Our studies enrich the genomic information available for the Lamiaceae family and provide a toolkit for discovering CYP450 genes involved in the flavonoid decoration.

Tandem duplication and sub-functionalization of clerodane diterpene synthase originate the blooming of clerodane diterpenoids in Scutellaria barbata

Scutellaria barbata is a traditional Chinese herb medicine and a major source of bioactive clerodane diterpenoids. However, barely clerodanes have been isolated from the closely related S. baicalensis. Here we assembled a chromosome-level genome of S. barbata and identified three class II clerodane diterpene synthases (SbarKPS1, SbarKPS2 and SbaiKPS1) from these two organisms. Using in vitro and in vivo assays, SbarKPS1 was characterized as a monofunctional (-)-kolavenyl diphosphate synthases ((-)-KPS), while SbarKPS2 and SbaiKPS1 produced major neo-cleroda-4(18),13E-dienyl diphosphate with small amount of (-)-KPP. SbarKPS1 and SbarKPS2 shared a high protein sequence identity and formed a tandem gene pair, indicating tandem duplication and sub-functionalization probably led to the evolution of monofunctional (-)-KPS in S. barbata. Additionally, SbarKPS1 and SbarKPS2 were primarily expressed in the leaves and flowers of S. barbata, which was consistent with the distribution of major clerodane diterpenoids scutebarbatine A and B. In contrast, SbaiKPS1 was barely expressed in any tissue of S. baicalensis. We further explored the downstream class I diTPS by functional characterizing of SbarKSL3 and SbarKSL4. Unfortunately, no dephosphorylated product was detected in the coupled assays with SbarKSL3/KSL4 and four class II diTPSs (SbarKPS1, SbarKPS2, SbarCPS2 and SbarCPS4) when a phosphatase inhibitor cocktail was included. Co-expression of SbarKSL3/KSL4 with class II diTPSs in yeast cells did not increase the yield of the corresponding dephosphorylated products, either. Together, these findings elucidated the involvement of two class II diTPSs in clerodane biosynthesis in S. barbata, while the class I diTPS is likely not responsible for the subsequent dephosphorylation step.