Co-Regulation of Clustered and Neo-Functionalized Genes in Plant-Specialized Metabolism

Comparative genomics points to tandem duplications of SAD gene clusters as drivers of increased α-linolenic (ω-3) content in S. hispanica seeds

Salvia hispanica L. (chia) is a source of abundant ω-3 polyunsaturated fatty acids (ω-3-PUFAs) that are highly beneficial to human health. The genomic basis for this accrued ω-3-PUFA content in this emerging crop was investigated through the assembly and comparative analysis of a chromosome-level reference genome for S. hispanica. The highly contiguous 321.5-Mbp genome assembly covering all six chromosomes enabled the identification of 32,922 protein-coding genes. Two whole-genome duplications (WGD) events were identified in the S. hispanica lineage. However, these WGD events could not be linked to the high α-linolenic acid (ALA, ω-3) accumulation in S. hispanica seeds based on phylogenomics. Instead, our analysis supports the hypothesis that evolutionary expansion through tandem duplications of specific lipid gene families, particularly the stearoyl-acyl carrier protein desaturase (ShSAD) gene family, is the main driver of the abundance of ω-3-PUFAs in S. hispanica seeds. The insights gained from the genomic analysis of S. hispanica will help establish a molecular breeding target that can be leveraged through genome editing techniques to increase ω-3 content in oil crops.

Key FAD2, FAD3, and SAD Genes Involved in the Fatty Acid Synthesis in Flax Identified Based on Genomic and Transcriptomic Data

FAD (fatty acid desaturase) and SAD (stearoyl-ACP desaturase) genes play key roles in the synthesis of fatty acids (FA) and determination of oil composition in flax (Linum usitatissimum L.). We searched for FAD and SAD genes in the most widely used flax genome of the variety CDC Bethune and three available long-read assembled flax genomes-YY5, 3896, and Atlant. We identified fifteen FAD2, six FAD3, and four SAD genes. Of all the identified genes, 24 were present in duplicated pairs. In most cases, two genes from a pair differed by a significant number of gene-specific SNPs (single nucleotide polymorphisms) or even InDels (insertions/deletions), except for FAD2a-1 and FAD2a-2, where only seven SNPs distinguished these genes. Errors were detected in the FAD2a-1, FAD2a-2, FAD3c-1, and FAD3d-2 sequences in the CDC Bethune genome assembly but not in the long-read genome assemblies. Expression analysis of the available transcriptomic data for different flax organs/tissues revealed that FAD2a-1, FAD2a-2, FAD3a, FAD3b, SAD3-1, and SAD3-2 were specifically expressed in embryos/seeds/capsules and could play a crucial role in the synthesis of FA in flax seeds. In contrast, FAD2b-1, FAD2b-2, SAD2-1, and SAD2-2 were highly expressed in all analyzed organs/tissues and could be involved in FA synthesis in whole flax plants. FAD2c-2, FAD2d-1, FAD3c-1, FAD3c-2, FAD3d-1, FAD3d-2, SAD3-1, and SAD3-2 showed differential expression under stress conditions-Fusarium oxysporum infection and drought. The obtained results are essential for research on molecular mechanisms of fatty acid synthesis, FAD and SAD editing, and marker-assisted and genomic selection for breeding flax varieties with a determined fatty acid composition of oil.

Species-specific 'specialized' genomic region provides the new insights into the functional genomics characterizing metabolic polymorphisms in plants

Functional genomics approaches with comparative omics analyses of wild-accessions and cultivars/wild species, as well as comparative genomic analyses in plant species focusing on gene clusters, have successfully detected key metabolic polymorphisms in plant specialized metabolism. In recent decades, (i) intra-species specific metabolic polymorphisms, (ii) new functionalization of tandem duplicated genes, and (iii) metabolic gene clusters were found as the main factors creating metabolic diversity of specialized metabolites in plants. However, given findings aware us that the identification of genes in plant specialized metabolism requires strategic approaches depending on the target metabolic pathways. The increasing availability of plant genome sequences and transcriptome data has facilitated inter-specific comparative analyses, including genomic analysis and gene co-expression network analysis. Here, we introduce functional genomics approaches with the integration of inter-/intra-species comparative metabolomics, their key roles in providing genomic signatures of metabolic evolution, and discuss future prospects of functional genomics on plant specialized metabolism.

Identification of P450 Candidates Associated with the Biosynthesis of Physalin-Class Compounds in Physalis angulata

The Physalis genus has long been used as traditional medicine in the treatment of various diseases. Physalins, the characteristic class of compounds in this genus, are major bioactive constituents. To date, the biogenesis of physalins remains largely unknown, except for the recently established knowledge that 24-methyldesmosterol is a precursor of physalin. To identify the genes encoding P450s that are putatively involved in converting 24-methyldesmosterol to physalins, a total of 306 P450-encoding unigenes were retrieved from our recently constructed P. angulata transcriptome. Extensive phylogenetic analysis proposed 21 P450s that might participate in physalin biosynthesis. To validate the candidates, we developed a virus-induced gene silencing (VIGS) system for P. angulata, and four P450 candidates were selected for the VIGS experiments. The reduction in the transcripts of the four P450 candidates by VIGS all led to decreased levels of physalin-class compounds in the P. angulata leaves. Thus, this study provides a number of P450 candidates that are likely associated with the biosynthesis of physalin-class compounds, forming a strong basis to reveal the unknown physalin biosynthetic pathway in the future.

Post-Cyclization Skeletal Rearrangements in Plant Triterpenoid Biosynthesis by a Pair of Branchpoint Isomerases

Triterpenoids possess potent biological activities, but their polycyclic skeletons are challenging to synthesize. The skeletal diversity of triterpenoids in plants is generated by oxidosqualene cyclases based on epoxide-triggered cationic rearrangement cascades. Normally, triterpenoid skeletons then remain unaltered during subsequent tailoring steps. In contrast, the highly modified triterpenoids found in Sapindales plants imply the existence of post-cyclization skeletal rearrangement enzymes that have not yet been found. We report here a biosynthetic pathway in Sapindales plants for the modification of already cyclized tirucallane triterpenoids, controlling the pathway bifurcation between different plant triterpenoid classes. Using a combination of bioinformatics, heterologous expression in plants and chemical analyses, we identified a cytochrome P450 monooxygenase and two isomerases which harness the epoxidation-rearrangement biosynthetic logic of triterpene cyclizations for modifying the tirucallane scaffold. The two isomerases share the same epoxide substrate made by the cytochrome P450 monooxygenase CYP88A154, but generate two different rearrangement products, one containing a cyclopropane ring. Our findings reveal a process for skeletal rearrangements of triterpenoids in nature that expands their scaffold diversity after the initial cyclization. In addition, the enzymes described here are crucial for the biotechnological production of limonoid, quassinoid, apoprotolimonoid, and glabretane triterpenoids.

Identification of a second 16-hydroxytabersonine-O-methyltransferase suggests an evolutionary relationship between alkaloid and flavonoid metabolisms in Catharanthus roseus

The medicinal plant Catharanthus roseus biosynthesizes many important drugs for human health, including the anticancer monoterpene indole alkaloids (MIAs) vinblastine and vincristine. Over the past decades, the continuous increase in pharmaceutical demand has prompted several research groups to characterize MIA biosynthetic pathways for considering future metabolic engineering processes of supply. In line with previous work suggesting that diversification can potentially occur at various steps along the vindoline branch, we were here interested in investigating the involvement of distinct isoforms of tabersonine-16-O-methyltransferase (16OMT) which plays a pivotal role in the MIA biosynthetic pathway. By combining homology searches based on the previously characterized 16OMT1, phylogenetic analyses, functional assays in yeast, and biochemical and in planta characterizations, we identified a second isoform of 16OMT, referred to as 16OMT2. 16OMT2 appears to be a multifunctional enzyme working on both MIA and flavonoid substrates, suggesting that a constrained evolution of the enzyme for accommodating the MIA substrate has probably occurred to favor the apparition of 16OMT2 from an ancestral specific flavonoid-O-methyltransferase. Since 16OMT1 and 16OMT2 displays a high sequence identity and similar kinetic parameters for 16-hydroxytabersonine, we postulate that 16OMT1 may result from a later 16OMT2 gene duplication accompanied by a continuous neofunctionalization leading to an almost complete loss of flavonoid O-methyltransferase activity. Overall, these results participate in increasing our knowledge on the evolutionary processes that have likely led to enzyme co-optation for MIA synthesis.

Flavan-3-ol Galloylation-Related Functional Gene Cluster and the Functional Diversification of SCPL Paralogs in Camellia sp

The high accumulation of galloylated flavan-3-ols in Camellia sp. is a noteworthy phenomenon. We identified a flavan-3-ol galloylation-related functional gene cluster in tannin-rich plant Camellia sp., which included UGT84A22 and SCPL-AT gene clusters. We investigated the possible correlation between the accumulation of metabolites and the expression of SCPL-ATs and UGT84A22. The results revealed that C. sinensis, C. ptilophylla, and C. oleifera accumulated galloylated cis-flavan-3-ols (EGCG), galloylated trans-flavan-3-ols (GCG), and hydrolyzed tannins, respectively; however, C. nitidissima did not accumulate any galloylated compounds. C. nitidissima exhibited no expression of SCPL-AT or UGT84A22, whereas the other three species of Camellia exhibited various expression patterns. This indicated that the functions of the paralogs of SCPL-AT vary. Enzymatic analysis revealed that SCPL5 was neofunctionalized as a noncatalytic chaperone paralog, a type of chaerone-like protein, associating with flavan-3-ol galloylation; moreover, CsSCPL4 was subfunctionalized in association with the galloylation of cis- and trans-flavan-3-ols. In C. nitidissima, an SCPL4 homolog was noted with mutations in two cysteine residues forming a disulfide bond, which suggested that this homolog was defunctionalized. The findings of this study improve our understanding of the functional diversification of SCPL paralogs in Camellia sp.

Pathogen Resistance Depending on Jacalin-Dirigent Chimeric Proteins Is Common among Poaceae but Absent in the Dicot Arabidopsis as Evidenced by Analysis of Homologous Single-Domain Proteins

MonocotJRLs are Poaceae-specific two-domain proteins that consist of a jacalin-related lectin (JRL) and a dirigent (DIR) domain which participate in multiple developmental processes, including disease resistance. For OsJAC1, a monocotJRL from rice, it has been confirmed that constitutive expression in transgenic rice or barley plants facilitates broad-spectrum disease resistance. In this process, both domains of OsJAC1 act cooperatively, as evidenced from experiments with artificially separated JRL- or DIR-domain-containing proteins. Interestingly, these chimeric proteins did not evolve in dicotyledonous plants. Instead, proteins with a single JRL domain, multiple JRL domains or JRL domains fused to domains other than DIR domains are present. In this study, we wanted to test if the cooperative function of JRL and DIR proteins leading to pathogen resistance was conserved in the dicotyledonous plant Arabidopsis thaliana. In Arabidopsis, we identified 50 JRL and 24 DIR proteins, respectively, from which seven single-domain JRL and two single-domain DIR candidates were selected. A single-cell transient gene expression assay in barley revealed that specific combinations of the Arabidopsis JRL and DIR candidates reduced the penetration success of barley powdery mildew. Strikingly, one of these pairs, AtJAX1 and AtDIR19, is encoded by genes located next to each other on chromosome one. However, when using natural variation and analyzing Arabidopsis ecotypes that express full-length or truncated versions of AtJAX1, the presence/absence of the full-length AtJAX1 protein could not be correlated with resistance to the powdery mildew fungus Golovinomyces orontii. Furthermore, an analysis of the additional JRL and DIR candidates in a bi-fluorescence complementation assay in Nicotiana benthamiana revealed no direct interaction of these JRL/DIR pairs. Since transgenic Arabidopsis plants expressing OsJAC1-GFP also did not show increased resistance to G. orontii, it was concluded that the resistance mediated by the synergistic activities of DIR and JRL proteins is specific for members of the Poaceae, at least regarding the resistance against powdery mildew. Arabidopsis lacks the essential components of the DIR-JRL-dependent resistance pathway.

Chromosome-scale assembly of the Moringa oleifera Lam. genome uncovers polyploid history and evolution of secondary metabolism pathways through tandem duplication

The African Orphan Crops Consortium (AOCC) selected the highly nutritious, fast growing and drought tolerant tree crop moringa (Moringa oleifera Lam.) as one of the first of 101 plant species to have its genome sequenced and a first draft assembly was published in 2019. Given the extensive uses and culture of moringa, often referred to as the multipurpose tree, we generated a significantly improved new version of the genome based on long-read sequencing into 14 pseudochromosomes equivalent to n = 14 haploid chromosomes. We leveraged this nearly complete version of the moringa genome to investigate main drivers of gene family and genome evolution that may be at the origin of relevant biological innovations including agronomical favorable traits. Our results reveal that moringa has not undergone any additional whole-genome duplication (WGD) or polyploidy event beyond the gamma WGD shared by all core eudicots. Moringa duplicates retained following that ancient gamma events are also enriched for functions commonly considered as dosage balance sensitive. Furthermore, tandem duplications seem to have played a prominent role in the evolution of specific secondary metabolism pathways including those involved in the biosynthesis of bioactive glucosinolate, flavonoid, and alkaloid compounds as well as of defense response pathways and might, at least partially, explain the outstanding phenotypic plasticity attributed to this species. This study provides a genetic roadmap to guide future breeding programs in moringa, especially those aimed at improving secondary metabolism related traits.

Genomic insights into the evolution of plant chemical defense

Plant trait evolution can be impacted by common mechanisms of genome evolution, including whole-genome and small-scale duplication, rearrangement, and selective pressures. With the increasing accessibility of genome sequencing for non-model species, comparative studies of trait evolution among closely related or divergent lineages have supported investigations into plant chemical defense. Plant defensive compounds include major chemical classes, such as terpenoids, alkaloids, and phenolics, and are used in primary and secondary plant functions. These include the promotion of plant health, facilitation of pollination, defense against pathogens, and responses to a rapidly changing climate. We discuss mechanisms of genome evolution and use examples from recent studies to impress a stronger understanding of the link between genotype and phenotype as it relates to the evolution of plant chemical defense. We conclude with considerations for how to leverage genomics, transcriptomics, metabolomics, and functional assays for studying the emergence and evolution of chemical defense systems.

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.

Using Gene Expression to Study Specialized Metabolism-A Practical Guide

Plants produce a vast array of chemical compounds that we use as medicines and flavors, but these compounds' biosynthetic pathways are still poorly understood. This paucity precludes us from modifying, improving, and mass-producing these specialized metabolites in suitable bioreactors. Many of the specialized metabolites are expressed in a narrow range of organs, tissues, and cell types, suggesting a tight regulation of the responsible biosynthetic pathways. Fortunately, with unprecedented ease of generating gene expression data and with >200,000 publicly available RNA sequencing samples, we are now able to study the expression of genes from hundreds of plant species. This review demonstrates how gene expression can elucidate the biosynthetic pathways by mining organ-specific genes, gene expression clusters, and applying various types of co-expression analyses. To empower biologists to perform these analyses, we showcase these analyses using recently published, user-friendly tools. Finally, we analyze the performance of co-expression networks and show that they are a valuable addition to elucidating multiple the biosynthetic pathways of specialized metabolism.

Terpenes and Terpenoids in Plants: Interactions with Environment and Insects

The interactions of plants with environment and insects are bi-directional and dynamic. Consequently, a myriad of mechanisms has evolved to engage organisms in different types of interactions. These interactions can be mediated by allelochemicals known as volatile organic compounds (VOCs) which include volatile terpenes (VTs). The emission of VTs provides a way for plants to communicate with the environment, including neighboring plants, beneficiaries (e.g., pollinators, seed dispersers), predators, parasitoids, and herbivores, by sending enticing or deterring signals. Understanding terpenoid distribution, biogenesis, and function provides an opportunity for the design and implementation of effective and efficient environmental calamity and pest management strategies. This review provides an overview of plant-environment and plant-insect interactions in the context of terpenes and terpenoids as important chemical mediators of these abiotic and biotic interactions.