Subtelomeric assembly of a multi-gene pathway for antimicrobial defense compounds in cereals

Complete biosynthesis of the potent vaccine adjuvant QS-21

QS-21 is a potent vaccine adjuvant currently sourced by extraction from the Chilean soapbark tree. It is a key component of human vaccines for shingles, malaria, coronavirus disease 2019 and others under development. The structure of QS-21 consists of a glycosylated triterpene scaffold coupled to a complex glycosylated 18-carbon acyl chain that is critical for immunostimulant activity. We previously identified the early pathway steps needed to make the triterpene glycoside scaffold; however, the biosynthetic route to the acyl chain, which is needed for stimulation of T cell proliferation, was unknown. Here, we report the biogenic origin of the acyl chain, characterize the series of enzymes required for its synthesis and addition and reconstitute the entire 20-step pathway in tobacco, thereby demonstrating the production of QS-21 in a heterologous expression system. This advance opens up unprecedented opportunities for bioengineering of vaccine adjuvants, investigating structure-activity relationships and understanding the mechanisms by which these compounds promote the human immune response.

A widely targeted metabolite modificomics strategy for modified metabolites identification in tomato

The structural and functional diversity of plant metabolites is largely created via chemical modification of a basic backbone. However, metabolite modifications in plants have still not been thoroughly investigated by metabolomics approaches. In this study, a widely targeted metabolite modificomics (WTMM) strategy was developed based on ultra-high performance liquid chromatography-quadrupole-linear ion trap (UHPLC-Q-Trap) and UHPLC-Q-Exactive-Orbitrap (UHPLC-QE-Orbitrap), which greatly improved the detection sensitivity and the efficiency of identification of modified metabolites. A metabolite modificomics study was carried out using tomato as a model, and over 34,000 signals with MS2 information were obtained from approximately 232 neutral loss transitions. Unbiased metabolite profiling was also performed by utilizing high-resolution mass spectrometry data to annotate a total of 2,118 metabolites with 125 modification types; of these, 165 modified metabolites were identified in this study. Next, the WTMM database was used to assess diseased tomato tissues and 29 biomarkers were analyzed. In summary, the WTMM strategy is not only capable of large-scale detection and quantitative analysis of plant-modified metabolites in plants, but also can be used for plant biomarker development.

Progress in Identification of UDP-Glycosyltransferases for Ginsenoside Biosynthesis

Ginsenosides, the primary pharmacologically active constituents of the Panax genus, have demonstrated a variety of medicinal properties, including anticardiovascular disease, cytotoxic, antiaging, and antidiabetes effects. However, the low concentration of ginsenosides in plants and the challenges associated with their extraction impede the advancement and application of ginsenosides. Heterologous biosynthesis represents a promising strategy for the targeted production of these natural active compounds. As representative triterpenoids, the biosynthetic pathway of the aglycone skeletons of ginsenosides has been successfully decoded. While the sugar moiety is vital for the structural diversity and pharmacological activity of ginsenosides, the mining of uridine diphosphate-dependent glycosyltransferases (UGTs) involved in ginsenoside biosynthesis has attracted a lot of attention and made great progress in recent years. In this paper, we summarize the identification and functional study of UGTs responsible for ginsenoside synthesis in both plants, such as Panax ginseng and Gynostemma pentaphyllum, and microorganisms including Bacillus subtilis and Saccharomyces cerevisiae. The UGT-related microbial cell factories for large-scale ginsenoside production are also mentioned. Additionally, we delve into strategies for UGT mining, particularly potential rapid screening or identification methods, providing insights and prospects. This review provides insights into the study of other unknown glycosyltransferases as candidate genetic elements for the heterologous biosynthesis of rare ginsenosides.

Subtelomeric plasticity contributes to gene family expansion in the human parasitic flatworm Schistosoma mansoni

BACKGROUND

The genomic region that lies between the telomere and chromosome body, termed the subtelomere, is heterochromatic, repeat-rich, and frequently undergoes rearrangement. Within this region, large-scale structural changes enable gene diversification, and, as such, large multicopy gene families are often found at the subtelomere. In some parasites, genes associated with proliferation, invasion, and survival are often found in these regions, where they benefit from the subtelomere's highly plastic, rapidly changing nature. The increasing availability of complete (or near complete) parasite genomes provides an opportunity to investigate these typically poorly defined and overlooked genomic regions and potentially reveal relevant gene families necessary for the parasite's lifestyle.

RESULTS

Using the latest chromosome-scale genome assembly and hallmark repeat richness observed at chromosome termini, we have identified and characterised the subtelomeres of Schistosoma mansoni, a metazoan parasitic flatworm that infects over 250 million people worldwide. Approximately 12% of the S. mansoni genome is classified as subtelomeric, and, in line with other organisms, we find these regions to be gene-poor but rich in transposable elements. We find that S. mansoni subtelomeres have undergone extensive interchromosomal recombination and that these sites disproportionately contribute to the 2.3% of the genome derived from segmental duplications. This recombination has led to the expansion of subtelomeric gene clusters containing 103 genes, including the immunomodulatory annexins and other gene families with unknown roles. The largest of these is a 49-copy plexin domain-containing protein cluster, exclusively expressed in the tegument-the tissue located at the host-parasite physical interface-of intramolluscan life stages.

CONCLUSIONS

We propose that subtelomeric regions act as a genomic playground for trial-and-error of gene duplication and subsequent divergence. Owing to the importance of subtelomeric genes in other parasites, gene families implicated in this subtelomeric expansion within S. mansoni warrant further characterisation for a potential role in parasitism.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Recent advances in triterpenoid pathway elucidation and engineering

Triterpenoids are among the most assorted class of specialized metabolites found in all the taxa of living organisms. Triterpenoids are the leading active ingredients sourced from plant species and are utilized in pharmaceutical and cosmetic industries. The triterpenoid precursor 2,3-oxidosqualene, which is biosynthesized via the mevalonate (MVA) pathway is structurally diversified by the oxidosqualene cyclases (OSCs) and other scaffold-decorating enzymes such as cytochrome P450 monooxygenases (P450s), UDP-glycosyltransferases (UGTs) and acyltransferases (ATs). A majority of the bioactive triterpenoids are harvested from the native hosts using the traditional methods of extraction and occasionally semi-synthesized. These methods of supply are time-consuming and do not often align with sustainability goals. Recent advancements in metabolic engineering and synthetic biology have shown prospects for the green routes of triterpenoid pathway reconstruction in heterologous hosts such as Escherichia coli, Saccharomyces cerevisiae and Nicotiana benthamiana, which appear to be quite promising and might lead to the development of alternative source of triterpenoids. The present review describes the biotechnological strategies used to elucidate complex biosynthetic pathways and to understand their regulation and also discusses how the advances in triterpenoid pathway engineering might aid in the scale-up of triterpenoid production in engineered hosts.

Knockout of CAFFEOYL-COA 3-O-METHYLTRANSFERASE 6/6L enhances the S/G ratio of lignin monomers and disease resistance in Nicotiana tabacum

Background

Nicotiana tabacum is an important economic crop, which is widely planted in the world. Lignin is very important for maintaining the physiological and stress-resistant functions of tobacco. However, higher lignin content will produce lignin gas, which is not conducive to the formation of tobacco quality. To date, how to precisely fine-tune lignin content or composition remains unclear.

Results

Here, we annotated and screened 14 CCoAOMTs in Nicotiana tabacum and obtained homozygous double mutants of CCoAOMT6 and CCoAOMT6L through CRSIPR/Cas9 technology. The phenotype showed that the double mutants have better growth than the wild type whereas the S/G ratio increased and the total sugar decreased. Resistance against the pathogen test and the extract inhibition test showed that the transgenic tobacco has stronger resistance to tobacco bacterial wilt and brown spot disease, which are infected by Ralstonia solanacearum and Alternaria alternata, respectively. The combined analysis of metabolome and transcriptome in the leaves and roots suggested that the changes of phenylpropane and terpene metabolism are mainly responsible for these phenotypes. Furthermore, the molecular docking indicated that the upregulated metabolites, such as soyasaponin Bb, improve the disease resistance due to highly stable binding with tyrosyl-tRNA synthetase targets in Ralstonia solanacearum and Alternaria alternata.

Conclusions

CAFFEOYL-COA 3-O-METHYLTRANSFERASE 6/6L can regulate the S/G ratio of lignin monomers and may affect tobacco bacterial wilt and brown spot disease resistance by disturbing phenylpropane and terpene metabolism in leaves and roots of Nicotiana tabacum, such as soyasaponin Bb.

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

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

Pangenome-based trajectories of intracellular gene transfers in Poaceae unveil high cumulation in Triticeae

Intracellular gene transfers (IGTs) between the nucleus and organelles, including plastids and mitochondria, constantly reshape the nuclear genome during evolution. Despite the substantial contribution of IGTs to genome variation, the dynamic trajectories of IGTs at the pangenomic level remain elusive. Here, we developed an approach, IGTminer, that maps the evolutionary trajectories of IGTs using collinearity and gene reannotation across multiple genome assemblies. We applied IGTminer to create a nuclear organellar gene (NOG) map across 67 genomes covering 15 Poaceae species, including important crops. The resulting NOGs were verified by experiments and sequencing data sets. Our analysis revealed that most NOGs were recently transferred and lineage specific and that Triticeae species tended to have more NOGs than other Poaceae species. Wheat (Triticum aestivum) had a higher retention rate of NOGs than maize (Zea mays) and rice (Oryza sativa), and the retained NOGs were likely involved in photosynthesis and translation pathways. Large numbers of NOG clusters were aggregated in hexaploid wheat during 2 rounds of polyploidization, contributing to the genetic diversity among modern wheat accessions. We implemented an interactive web server to facilitate the exploration of NOGs in Poaceae. In summary, this study provides resources and insights into the roles of IGTs in shaping interspecies and intraspecies genome variation and driving plant genome evolution.

Subtelomeric 5-enolpyruvylshikimate-3-phosphate synthase copy number variation confers glyphosate resistance in Eleusine indica

Genomic structural variation (SV) has profound effects on organismal evolution; often serving as a source of novel genetic variation. Gene copy number variation (CNV), one type of SV, has repeatedly been associated with adaptive evolution in eukaryotes, especially with environmental stress. Resistance to the widely used herbicide, glyphosate, has evolved through target-site CNV in many weedy plant species, including the economically important grass, Eleusine indica (goosegrass); however, the origin and mechanism of these CNVs remain elusive in many weed species due to limited genetic and genomic resources. To study this CNV in goosegrass, we present high-quality reference genomes for glyphosate-susceptible and -resistant goosegrass lines and fine-assembles of the duplication of glyphosate's target site gene 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS). We reveal a unique rearrangement of EPSPS involving chromosome subtelomeres. This discovery adds to the limited knowledge of the importance of subtelomeres as genetic variation generators and provides another unique example for herbicide resistance evolution.

Triticeae crop genome biology: an endless frontier

Triticeae, the wheatgrass tribe, includes several major cereal crops and their wild relatives. Major crops within the Triticeae are wheat, barley, rye, and oat, which are important for human consumption, animal feed, and rangeland protection. Species within this tribe are known for their large genomes and complex genetic histories. Powered by recent advances in sequencing technology, researchers worldwide have made progress in elucidating the genomes of Triticeae crops. In addition to assemblies of high-quality reference genomes, pan-genome studies have just started to capture the genomic diversities of these species, shedding light on our understanding of the genetic basis of domestication and environmental adaptation of Triticeae crops. In this review, we focus on recent signs of progress in genome sequencing, pan-genome analyses, and resequencing analysis of Triticeae crops. We also propose future research avenues in Triticeae crop genomes, including identifying genome structure variations, the association of genomic regions with desired traits, mining functions of the non-coding area, introgression of high-quality genes from wild Triticeae resources, genome editing, and integration of genomic resources.

Genome-wide comparative analysis of the valine glutamine motif containing genes in four Ipomoea species

BACKGROUND

Genes with valine glutamine (VQ) motifs play an essential role in plant growth, development, and resistance to biotic and abiotic stresses. However, little information on the VQ genes in sweetpotato and other Ipomoea species is available.

RESULTS

This study identified 55, 58, 50 and 47 VQ genes from sweetpotato (I. batatas), I.triflida, I. triloba and I. nil, respectively. The phylogenetic analysis revealed that the VQ genes formed eight clades (I-VII), and the members in the same group exhibited similar exon-intron structure and conserved motifs distribution. The distribution of the VQ genes among the chromosomes of Ipomoea species was disproportional, with no VQ genes mapped on a few of each species' chromosomes. Duplication analysis suggested that segmental duplication significantly contributes to their expansion in sweetpotato, I.trifida, and I.triloba, while the segmental and tandem duplication contributions were comparable in I.nil. Cis-regulatory elements involved in stress responses, such as W-box, TGACG-motif, CGTCA-motif, ABRE, ARE, MBS, TCA-elements, LTR, and WUN-motif, were detected in the promoter regions of the VQ genes. A total of 30 orthologous groups were detected by syntenic analysis of the VQ genes. Based on the analysis of RNA-seq datasets, it was found that the VQ genes are expressed distinctly among different tissues and hormone or stress treatments. A total of 40 sweetpotato differentially expressed genes (DEGs) refer to biotic (sweetpotato stem nematodes and Ceratocystis fimbriata pathogen infection) or abiotic (cold, salt and drought) stress treatments were detected. Moreover, IbVQ8, IbVQ25 and IbVQ44 responded to the five stress treatments and were selected for quantitative reverse-transcription polymerase chain reaction (qRT-PCR) analysis, and the results were consistent with the transcriptome analysis.

CONCLUSIONS

Our study may provide new insights into the evolution of VQ genes in the four Ipomoea genomes and contribute to the future molecular breeding of sweetpotatoes.

Elucidation of the pathway for biosynthesis of saponin adjuvants from the soapbark tree

The Chilean soapbark tree (Quillaja saponaria) produces soap-like molecules called QS saponins that are important vaccine adjuvants. These highly valuable compounds are sourced by extraction from the bark, and their biosynthetic pathway is unknown. Here, we sequenced the Q. saponaria genome. Through genome mining and combinatorial expression in tobacco, we identified 16 pathway enzymes that together enable the production of advanced QS pathway intermediates that represent a bridgehead for adjuvant bioengineering. We further identified the enzymes needed to make QS-7, a saponin with excellent therapeutic properties and low toxicity that is present in low abundance in Q. saponaria bark extract. Our results enable the production of Q. saponaria vaccine adjuvants in tobacco and open the way for new routes to access and engineer natural and new-to-nature immunostimulants.

Population structure analysis and genome-wide association study of a hexaploid oat landrace and cultivar collection

INTRODUCTION

Oat (Avena sativa L.) is an important cereal crop grown worldwide for grain and forage, owing to its high adaptability to diverse environments. However, the genetic and genomics research of oat is lagging behind that of other staple cereal crops.

METHODS

In this study, a collection of 288 oat lines originating worldwide was evaluated using 2,213 single nucleotide polymorphism (SNP) markers obtained from an oat iSelect 6K-beadchip array to study its genetic diversity, population structure, and linkage disequilibrium (LD) as well as the genotype-phenotype association for hullessness and lemma color.

RESULTS

The average gene diversity and polymorphic information content (PIC) were 0.324 and 0.262, respectively. The first three principal components (PCs) accounted for 30.33% of the genetic variation, indicating that the population structure of this panel of oat lines was stronger than that reported in most previous studies. In addition, accessions could be classified into two subpopulations using a Bayesian clustering approach, and the clustering pattern of accessions was closely associated with their region of origin. Additionally, evaluation of LD decay using 2,143 mapped markers revealed that the intrachromosomal whole-genome LD decayed rapidly to a critical r² value of 0.156 for marker pairs separated by a genetic distance of 1.41 cM. Genome-wide association study (GWAS) detected six significant associations with the hullessness trait. Four of these six markers were located on the Mrg21 linkage group between 194.0 and 205.7 cM, while the other two significant markers mapped to Mrg05 and Mrg09. Three significant SNPs, showing strong association with lemma color, were located on linkage groups Mrg17, Mrg18, and Mrg20.

DISCUSSION

Our results discerned relevant patterns of genetic diversity, population structure, and LD among members of a worldwide collection of oat landraces and cultivars proposed to be 'typical' of the Qinghai-Tibetan Plateau. These results have important implications for further studies on association mapping and practical breeding in high-altitude oat.

Comparative analysis of the MYB gene family in seven Ipomoea species

The MYB transcription factors regulate plant growth, development, and defense responses. However, information about the MYB gene family in Ipomoea species is rare. Herein, we performed a comprehensive genome-wide comparative analysis of this gene family among seven Ipomoea species, sweet potato (I. batatas), I. trifida, I. triloba, I. nil, I. purpurea, I. cairica, and I. aquatic, and identified 296, 430, 411, 291, 226, 281, and 277 MYB genes, respectively. The identified MYB genes were classified into five types: 1R-MYB (MYB-related), 2R-MYB (R2R3-MYB), 3R-MYB (R1R2R3-MYB), 4R-MYB, and 5R-MYB, and the MYB-related or R2R3-MYB type was the most abundant MYB genes in the seven species. The Ipomoea MYB genes were classed into distinct subgroups based on the phylogenetic topology and the classification of the MYB superfamily in Arabidopsis. Analysis of gene structure and protein motifs revealed that members within the same phylogenetic group presented similar exon/intron and motif organization. The identified MYB genes were unevenly mapped on the chromosomes of each Ipomoea species. Duplication analysis indicated that segmental and tandem duplications contribute to expanding the Ipomoea MYB genes. Non-synonymous substitution (Ka) to synonymous substitution (Ks) [Ka/Ks] analysis showed that the duplicated Ipomoea MYB genes are mainly under purifying selection. Numerous cis-regulatory elements related to stress responses were detected in the MYB promoters. Six sweet potato transcriptome datasets referring to abiotic and biotic stresses were analyzed, and MYB different expression genes' (DEGs') responses to stress treatments were detected. Moreover, 10 sweet potato MYB DEGs were selected for qRT-PCR analysis. The results revealed that four responded to biotic stress (stem nematodes and Ceratocystis fimbriata pathogen infection) and six responded to the biotic stress (cold, drought, and salt). The results may provide new insights into the evolution of MYB genes in the Ipomoea genome and contribute to the future molecular breeding of sweet potatoes.

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.

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.

Complex scaffold remodeling in plant triterpene biosynthesis

Triterpenes with complex scaffold modifications are widespread in the plant kingdom. Limonoids are an exemplary family that are responsible for the bitter taste in citrus (e.g., limonin) and the active constituents of neem oil, a widely used bioinsecticide (e.g., azadirachtin). Despite the commercial value of limonoids, a complete biosynthetic route has not been described. We report the discovery of 22 enzymes, including a pair of neofunctionalized sterol isomerases, that catalyze 12 distinct reactions in the total biosynthesis of kihadalactone A and azadirone, products that bear the signature limonoid furan. These results enable access to valuable limonoids and provide a template for discovery and reconstitution of triterpene biosynthetic pathways in plants that require multiple skeletal rearrangements and oxidations.

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

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

Emerging mechanistic insights into the regulation of specialized metabolism in plants

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

The Vinca minor genome highlights conserved evolutionary traits in monoterpene indole alkaloid synthesis

Vinca minor, also known as the lesser periwinkle, is a well-known species from the Apocynaceae, native to central and southern Europe. This plant synthesizes monoterpene indole alkaloids, which are a class of specialized metabolites displaying a wide range of bioactive- and pharmacologically important properties. Within the almost 50 monoterpene indole alkaloids it produces, V. minor mainly accumulates vincamine, which is commercially used as a nootropic. Using a combination of Oxford Nanopore Technologies long read- and Illumina short-read sequencing, a 679,098 Mb V. minor genome was assembled into 296 scaffolds with an N50 scaffold length of 6 Mb, and encoding 29,624 genes. These genes were functionally annotated and used in a comparative genomic analysis to establish gene families and to investigate gene family expansion and contraction across the phylogenetic tree. Furthermore, homology-based monoterpene indole alkaloid gene predictions together with a metabolic analysis across 4 different V. minor tissue types guided the identification of candidate monoterpene indole alkaloid genes. These candidates were finally used to identify monoterpene indole alkaloid gene clusters, which combined with synteny analysis allowed for the discovery of a functionally validated vincadifformine-16-hydroxylase, reinforcing the potential of this dataset for monoterpene indole alkaloids gene discovery. It is expected that access to these resources will facilitate the elucidation of unknown monoterpene indole alkaloid biosynthetic routes with the potential of transferring these pathways to heterologous expression systems for large-scale monoterpene indole alkaloid production.

Oat chromosome and genome evolution defined by widespread terminal intergenomic translocations in polyploids

Structural chromosome rearrangements involving translocations, fusions and fissions lead to evolutionary variation between species and potentially reproductive isolation and variation in gene expression. While the wheats (Triticeae, Poaceae) and oats (Aveneae) all maintain a basic chromosome number of x=7, genomes of oats show frequent intergenomic translocations, in contrast to wheats where these translocations are relatively rare. We aimed to show genome structural diversity and genome relationships in tetraploid, hexaploid and octoploid Avena species and amphiploids, establishing patterns of intergenomic translocations across different oat taxa using fluorescence in situ hybridization (FISH) with four well-characterized repetitive DNA sequences: pAs120, AF226603, Ast-R171 and Ast-T116. In A. agadiriana (2n=4x=28), the selected probes hybridized to all chromosomes indicating that this species originated from one (autotetraploid) or closely related ancestors with the same genomes. Hexaploid amphiploids were confirmed as having the genomic composition AACCDD, while octoploid amphiploids showed three different genome compositions: AACCCCDD, AAAACCDD or AABBCCDD. The A, B, C, and D genomes of oats differ significantly in their involvement in non-centromeric, intercalary translocations. There was a predominance of distal intergenomic translocations from the C- into the D-genome chromosomes. Translocations from A- to C-, or D- to C-genome chromosomes were less frequent, proving that at least some of the translocations in oat polyploids are non-reciprocal. Rare translocations from A- to D-, D- to A- and C- to B-genome chromosomes were also visualized. The fundamental research has implications for exploiting genomic biodiversity in oat breeding through introgression from wild species potentially with contrasting chromosomal structures and hence deleterious segmental duplications or large deletions in amphiploid parental lines.

Emerging trends in genomic and epigenomic regulation of plant specialised metabolism

Regulation of specialised metabolism genes is multilayered and complex, influenced by an array of genomic, epigenetic and epigenomic mechanisms. Here, we review the most recent knowledge in this field, drawing from discoveries in several plant species. Our aim is to improve understanding of how plant genome structure and function influence specialised metabolism. We also highlight key areas for future exploration. Gene regulatory mechanisms influencing specialised metabolism include gene duplication and neo-functionalization, conservation of operon-like clusters of specialised metabolism genes, local chromatin modifications, and the organisation of higher order chromatin structures within the nucleus. Genomic and epigenomic research to-date in the discipline have focused on a relatively small number of plant species, primarily at whole organ or tissue level. This is largely due to the technical demands of the experimental methods needed. However, a high degree of cell-type specificity of function exists in specialised metabolism, driven by similarly specific gene regulation. In this review we focus on the genomic characteristics of genes that are found in different types of clusters within the genome. We propose that acquisition of cell-resolution epigenomic datasets in emerging models, such as the glandular trichomes of Cannabis sativa, will yield important advances. Data such as chromatin accessibility and histone modification profiles can pinpoint which regulatory sequences are active in individual cell types and at specific times in development. These could provide fundamental biological insight as well as novel targets for genetic engineering and crop improvement.

Plant metabolic gene clusters in the multi-omics era

Secondary metabolism in plants gives rise to a vast array of small-molecule natural products. The discovery of operon-like gene clusters in plants has provided a new perspective on the evolution of specialized metabolism and the opportunity to rapidly advance the metabolic engineering of natural product production. Here, we review historical aspects of the study of plant metabolic gene clusters as well as general strategies for identifying plant metabolic gene clusters in the multi-omics era. We also emphasize the exploration of their natural variation and evolution, as well as new strategies for the prospecting of plant metabolic gene clusters and a deeper understanding of how their structure influences their function.

Integrative omics approaches for biosynthetic pathway discovery in plants

Covering: up to 2022With the emergence of large amounts of omics data, computational approaches for the identification of plant natural product biosynthetic pathways and their genetic regulation have become increasingly important. While genomes provide clues regarding functional associations between genes based on gene clustering, metabolome mining provides a foundational technology to chart natural product structural diversity in plants, and transcriptomics has been successfully used to identify new members of their biosynthetic pathways based on coexpression. Thus far, most approaches utilizing transcriptomics and metabolomics have been targeted towards specific pathways and use one type of omics data at a time. Recent technological advances now provide new opportunities for integration of multiple omics types and untargeted pathway discovery. Here, we review advances in plant biosynthetic pathway discovery using genomics, transcriptomics, and metabolomics, as well as recent efforts towards omics integration. We highlight how transcriptomics and metabolomics provide complementary information to link genes to metabolites, by associating temporal and spatial gene expression levels with metabolite abundance levels across samples, and by matching mass-spectral features to enzyme families. Furthermore, we suggest that elucidation of gene regulatory networks using time-series data may prove useful for efforts to unwire the complexities of biosynthetic pathway components based on regulatory interactions and events.

Architecture, Chromatin and Gene Organization of Toxoplasma gondii Subtelomeres

Subtelomeres (ST) are chromosome regions that separate telomeres from euchromatin and play relevant roles in various biological processes of the cell. While their functions are conserved, ST structure and genetic compositions are unique to each species. This study aims to identify and characterize the subtelomeric regions of the 13 Toxoplasma gondii chromosomes of the Me49 strain. Here, STs were defined at chromosome ends based on poor gene density. The length of STs ranges from 8.1 to 232.4 kbp, with a gene density of 0.049 genes/kbp, lower than the Me49 genome (0.15 kbp). Chromatin organization showed that H3K9me3, H2A.X, and H3.3 are highly enriched near telomeres and the 5′ end of silenced genes, decaying in intensity towards euchromatin. H3K4me3 and H2A.Z/H2B.Z are shown to be enriched in the 5′ end of the ST genes. Satellite DNA was detected in almost all STs, mainly the sat350 family and a novel satellite named sat240. Beyond the STs, only short dispersed fragments of sat240 and sat350 were found. Within STs, there were 12 functional annotated genes, 59 with unknown functions (Hypothetical proteins), 15 from multigene FamB, and 13 from multigene family FamC. Some genes presented low interstrain synteny associated with the presence of satellite DNA. Orthologues of FamB and FamC were also detected in Neospora caninum and Hammondia hammondi. A re-analysis of previous transcriptomic data indicated that ST gene expression is strongly linked to the adaptation to different situations such as extracellular passage (evolve and resequencing study) and changes in metabolism (lack of acetyl-CoA cofactor). In conclusion, the ST region of the T. gondii chromosomes was defined, the STs genes were determined, and it was possible to associate them with high interstrain plasticity and a role in the adaptability of T. gondii to environmental changes.

GGDB: A Grameneae genome alignment database of homologous genes hierarchically related to evolutionary events

The genomes of Gramineae plants have been preferentially sequenced owing to their economic value. These genomes are often quite complex, for example harboring many duplicated genes, and are the main source of genetic innovation and often the result of recurrent polyploidization. Deciphering these complex genome structures and linking duplicated genes to specific polyploidization events are important for understanding the biology and evolution of plants. However, efforts have been hampered by the complexity of analyzing these genomes. Here, we analyzed 29 well-assembled and up-to-date Gramineae genome sequences by hierarchically relating duplicated genes in collinear regions to specific polyploidization or speciation events. We separated duplicated genes produced by each event, established lists of paralogous and orthologous genes, and ultimately constructed an online database, GGDB (http://www.grassgenome.com/). Homologous gene lists from each plant and between plants can be displayed, searched, and downloaded from the database. Interactive comparison tools are deployed to demonstrate homology among user-selected plants and to draw genome-scale or local alignment figures and gene-based phylogenetic trees corrected by exploiting gene collinearity. Using these tools and figures, users can easily detect structural changes in genomes and explore the effects of paleo-polyploidy on crop genome structure and function. The GGDB will provide a useful platform for improving our understanding of genome changes and functional innovation in Gramineae plants.

Development of SSR Markers Linked to Stress Responsive Genes along Tomato Chromosome 3 (Solanum lycopersicum L.)

This study aimed to develop novel SSR markers in tomato. Several BAC clones along chromosome 3 in tomato were selected based on their content. The criteria was the availability of genes, either directly or indirectly related to stress response (drought, salinity, and heat) in tomato. A total of 20 novel in silico SSR markers were developed and 96 important nearby genes were identified. The identified nearby genes represent different tomato genes involved in plant growth and development and biotic and abiotic stress tolerance. The developed SSR markers were assessed using tomato landraces. A total of 29 determinate and semi-determinate local tomato landraces collected from diverse environments were utilized. A total of 33 alleles with mean of 1.65 alleles per locus were scored, showing 100% polymorphic patterns, with a mean of 0.18 polymorphism information content (PIC) values. The mean of observed and expected heterozygosity were 0.19 and 0.24, respectively. The mean value of the Jaccard similarity index was used for clustering the landraces. The developed microsatellite markers showed potential to assess genetic variability among tomato landraces. The genetic distance information reported in this study can be used by breeders in future genetic improvement of tomato for tolerance against diverse stresses.

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.

Functions and biosynthesis of plant signaling metabolites mediating plant-microbe interactions

Covering: 2015-2022Plants and microbes have coevolved since their appearance, and their interactions, to some extent, define plant health. A reasonable fraction of small molecules plants produced are involved in mediating plant-microbe interactions, yet their functions and biosynthesis remain fragmented. The identification of these compounds and their biosynthetic genes will open up avenues for plant fitness improvement by manipulating metabolite-mediated plant-microbe interactions. Herein, we integrate the current knowledge on their chemical structures, bioactivities, and biosynthesis with the view of providing a high-level overview on their biosynthetic origins and evolutionary trajectory, and pinpointing the yet unknown and key enzymatic steps in diverse biosynthetic pathways. We further discuss the theoretical basis and prospects for directing plant signaling metabolite biosynthesis for microbe-aided plant health improvement in the future.

Convergent gene clusters underpin hyperforin biosynthesis in St John's wort

The meroterpenoid hyperforin is responsible for the antidepressant activity of St John's wort extracts, but the genes controlling its biosynthesis are unknown. Using genome mining and biochemical work, we characterize two biosynthetic gene clusters (BGCs) that encode the first three steps in the biosynthesis of hyperforin precursors. The findings of syntenic and phylogenetic analyses reveal the parallel assembly of the two BGCs. The syntenous BGC in Mesua ferrea indicates that the first cluster was assembled before the divergence of the Hypericaceae and Calophyllaceae families. The assembly of the second cluster is the result of a coalescence of genomic fragments after a major duplication event. The differences between the two BGCs - in terms of gene expression, response to methyl jasmonate, substrate specificity and subcellular localization of key enzymes - suggest that the presence of the two clusters could serve to generate separate pools of precursors. The parallel assembly of two BGCs with similar compositions in a single plant species is uncommon, and our work provides insights into how and when these gene clusters form. Our discovery helps to advance our understanding of the evolution of plant specialized metabolism and its genomic organization. Additionally, our results offer a foundation from which hyperforin biosynthesis can be more fully understood, and which can be used in future metabolic engineering applications.

The genome of the recretohalophyte Limonium bicolor provides insights into salt gland development and salinity adaptation during terrestrial evolution

Halophytes have evolved specialized strategies to cope with high salinity. The extreme halophyte sea lavender (Limonium bicolor) lacks trichomes but possesses salt glands on its epidermis that can excrete harmful ions, such as sodium, to avoid salt damage. Here, we report a high-quality, 2.92-Gb, chromosome-scale L. bicolor genome assembly based on a combination of Illumina short reads, single-molecule, real-time long reads, chromosome conformation capture (Hi-C) data, and Bionano genome maps, greatly enriching the genomic information on recretohalophytes with multicellular salt glands. Although the L. bicolor genome contains genes that show similarity to trichome fate genes from Arabidopsis thaliana, it lacks homologs of the decision fate genes GLABRA3, ENHANCER OF GLABRA3, GLABRA2, TRANSPARENT TESTA GLABRA2, and SIAMESE, providing a molecular explanation for the absence of trichomes in this species. We identified key genes (LbHLH and LbTTG1) controlling salt gland development among classical trichome homologous genes and confirmed their roles by showing that their mutations markedly disrupted salt gland initiation, salt secretion, and salt tolerance, thus offering genetic support for the long-standing hypothesis that salt glands and trichomes may share a common origin. In addition, a whole-genome duplication event occurred in the L. bicolor genome after its divergence from Tartary buckwheat and may have contributed to its adaptation to high salinity. The L. bicolor genome resource and genetic evidence reported in this study provide profound insights into plant salt tolerance mechanisms that may facilitate the engineering of salt-tolerant crops.

Beyond a reference genome: pangenomes and population genomics of underutilized and orphan crops for future food and nutrition security

Underutilized crops are, by definition, under-researched compared to staple crops yet come with traits that may be especially important given climate change and the need to feed a globally increasing population. These crops are often stress-tolerant, and this combined with unique and beneficial nutritional profiles. Whilst progress is being made by generating reference genome sequences, in this Tansley Review, we show how this is only the very first step. We advocate that going 'beyond a reference genome' should be a priority, as it is only at this stage one can identify the specific genes and the adaptive alleles that underpin the valuable traits. We sum up how population genomic and pangenomic approaches have led to the identification of stress- and disease-tolerant alleles in staple crops and compare this to the small number of examples from underutilized crops. We also demonstrate how previously underutilized crops have benefitted from genomic advances and that many breeding targets in underutilized crops are often well studied in staple crops. This cross-crop population-level resequencing could lead to an understanding of the genetic basis of adaptive traits in underutilized crops. This level of investment may be crucial for fully understanding the value of these crops before they are lost.

Pathogen-induced biosynthetic pathways encode defense-related molecules in bread wheat

Wheat is a widely grown food crop that suffers major yield losses due to attack by pests and pathogens. A better understanding of biotic stress responses in wheat is thus of major importance. The recently assembled bread wheat genome coupled with extensive transcriptomic resources provides unprecedented new opportunities to investigate responses to pathogen challenge. Here, we analyze gene coexpression networks to identify modules showing consistent induction in response to pathogen exposure. Within the top pathogen-induced modules, we identify multiple clusters of physically adjacent genes that correspond to six pathogen-induced biosynthetic pathways that share a common regulatory network. Functional analysis reveals that these pathways, all of which are encoded by biosynthetic gene clusters, produce various different classes of compounds—namely, flavonoids, diterpenes, and triterpenes, including the defense-related compound ellarinacin. Through comparative genomics, we also identify associations with the known rice phytoalexins momilactones, as well as with a defense-related gene cluster in the grass model plant Brachypodium distachyon. Our results significantly advance the understanding of chemical defenses in wheat and open up avenues for enhancing disease resistance in this agriculturally important crop. They also exemplify the power of transcriptional networks to discover the biosynthesis of chemical defenses in plants with large, complex genomes.

Computational prediction of plant metabolic pathways

Uncovering genes encoding enzymes responsible for the biosynthesis of diverse plant metabolites is essential for metabolic engineering and production of plant metabolite-derived medicine. With the availability of multi-omics data for an ever-increasing number of plant species and the development of computational approaches, the metabolic pathways of many important plant compounds can be predicted, complementing a more traditional genetic and/or biochemical approach. Here, we summarize recent progress in predicting plant metabolic pathways using genome, transcriptome, proteome, interactome, and/or metabolome data, and the utility of integrating these data with machine learning to further improve metabolic pathway predictions.

A conserved mechanism affecting hydride shifting and deprotonation in the synthesis of hopane triterpenes as compositions of wax in oat

Hopanoids are a group of biologically important triterpene scaffolds found in nature, but the discovery of hopane-type triterpene synthases in plants has not been reported. We discovered two types of triterpene synthases synthesizing hopanoid skeletons from monocot and dicot plants and elucidated a mechanism involving the deprotonation at different sites by site-directed mutagenesis experiments and the quantum mechanics and molecular mechanics calculation. Our results provide a genetic element for synthesizing biologically active hopane-type triterpenoids and serve as a foundation for studying the molecular mechanisms of methyl and hydride transfer in the triterpene cyclization mechanism.

Triterpenoids are biologically active metabolites synthesized from a common linear precursor catalyzed by 2,3-oxidosqualene cyclases (OSCs) to form diverse triterpenoid skeletons. OSCs corresponding to many discovered triterpene alcohols in nature have not been functionally and mechanistically characterized due to the diversity of chemical structures and complexity of the cyclization mechanism. We carried out a genome-wide investigation of OSCs from Avena strigosa and discovered two triterpene synthases, namely, AsHS1 and AsHS2, using a Nicotiana benthamiana expression system. These synthases produce hopenol B and hop-17(21)-en-3β-ol, which are components of surface wax in oat panicles and sheathes, respectively. We demonstrated that substitutions of two to three amino acid residues in AsHS1 with corresponding residues from AsHS2 allowed it to be completely converted into a hop-17(21)-en-3β-ol synthase. AsHS2 mutants with a substitution at site 410 could synthesize hopenol B alone or mixed with a side product isomotiol. The combined quantum mechanics and molecular mechanics calculation demonstrated that the side chain size of the residue at site 410 regulated the relative orientations between the hopyl C22 cation and Phe257, leading to a difference in deprotonation positions through providing or not providing cation–π interaction between the aromatic ring of F257 and the carbocation intermediate. A similar mechanism could be applied to a hopenol B synthase from a dicotyledonous plant Aquilegia. This study provided mechanistic insight into triterpenoid synthesis and discovered key amino acid residues acting on hydride transfer and a deprotonation site to differentiate between hopane-type scaffolds in diverse plant species.

Chromosome-level assembly of the Neolamarckia cadamba genome provides insights into the evolution of cadambine biosynthesis

Neolamarckia cadamba (Roxb.), a close relative of Coffea canephora and Ophiorrhiza pumila, is an important traditional medicine in Southeast Asia. Three major glycosidic monoterpenoid indole alkaloids (MIAs), cadambine and its derivatives 3β-isodihydrocadambine and 3β-dihydrocadambine, accumulate in the bark and leaves, and exhibit antimalarial, antiproliferative, antioxidant, anticancer and anti-inflammatory activities. Here, we report a chromosome-scale N. cadamba genome, with 744.5 Mb assembled into 22 pseudochromosomes with contig N50 and scaffold N50 of 824.14 Kb and 29.20 Mb, respectively. Comparative genomic analysis of N. cadamba with Co. canephora revealed that N. cadamba underwent a relatively recent whole-genome duplication (WGD) event after diverging from Co. canephora, which contributed to the evolution of the MIA biosynthetic pathway. We determined the key intermediates of the cadambine biosynthetic pathway and further showed that NcSTR1 catalyzed the synthesis of strictosidine in N. cadamba. A new component, epoxystrictosidine (C27H34N2O10, m/z 547.2285), was identified in the cadambine biosynthetic pathway. Combining genome-wide association study (GWAS), population analysis, multi-omics analysis and metabolic gene cluster prediction, this study will shed light on the evolution of MIA biosynthetic pathway genes. This N. cadamba reference sequence will accelerate the understanding of the evolutionary history of specific metabolic pathways and facilitate the development of tools for enhancing bioactive productivity by metabolic engineering in microbes or by molecular breeding in plants.

Genome-wide comparative analysis of the nucleotide-binding site-encoding genes in four Ipomoea species

The nucleotide-binding site (NBS)-encoding gene is a major type of resistance (R) gene, and its diverse evolutionary patterns were analyzed in different angiosperm lineages. Until now, no comparative studies have been done on the NBS encoding genes in Ipomoea species. In this study, various numbers of NBS-encoding genes were identified across the whole genome of sweet potato (Ipomoea batatas) (#889), Ipomoea trifida (#554), Ipomoea triloba (#571), and Ipomoea nil (#757). Gene analysis showed that the CN-type and N-type were more common than the other types of NBS-encoding genes. The phylogenetic analysis revealed that the NBS-encoding genes formed three monophyletic clades: CNL, TNL, and RNL, which were distinguished by amino acid motifs. The distribution of the NBS-encoding genes among the chromosomes was non-random and uneven; 83.13, 76.71, 90.37, and 86.39% of the genes occurred in clusters in sweet potato, I. trifida, I. triloba, and I. nil, respectively. The duplication pattern analysis reveals the presence of higher segmentally duplicated genes in sweet potatoes than tandemly duplicated ones. The opposite trend was found for the other three species. A total of 201 NBS-encoding orthologous genes were found to form synteny gene pairs between any two of the four Ipomea species, suggesting that each of the synteny gene pairs was derived from a common ancestor. The gene expression patterns were acquired by analyzing using the published datasets. To explore the candidate resistant genes in sweet potato, transcriptome analysis has been carried out using two resistant (JK20 and JK274) and susceptible cultivars (Tengfei and Santiandao) of sweet potato for stem nematodes and Ceratocystis fimbriata pathogen, respectively. A total of 11 differentially expressed genes (DEGs) were found in Tengfei and JK20 for stem nematodes and 19 DEGs in Santiandao and JK274 for C. fimbriata. Moreover, six DEGs were further selected for quantitative reverse-transcription polymerase chain reaction (qRT-PCR) analysis, and the results were consistent with the transcriptome analysis. The results may provide new insights into the evolution of NBS-encoding genes in the Ipomoea genome and contribute to the future molecular breeding of sweet potatoes.

Characterization of a Highly Selective 2″-O-Galactosyltransferase from Trollius chinensis and Structure-Guided Engineering for Improving UDP-Glucose Selectivity

A novel 2″-O-galactosyltransferase TcOGT4 was discovered from Trollius chinensis. TcOGT4 could regio-specifically catalyze 2″-O-galactosylation of 17 flavone 8-C-β-d-glucosides and shows high preference for UDP-Gal. Molecular docking indicated that Pro361 may play a key role in sugar donor selectivity, and the P361W mutant exhibited significantly enhanced selectivity toward UDP-Glc. A total of 21 products including 17 new compounds were obtained, and 5 of them showed potent COX-2 inhibitory activities. TcOGT4 is the first reported 2″-O-galactosyltransferase for flavone C-glycosides, and could be a powerful biocatalyst to synthesize bioactive flavone glycosides.

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.