Metabolic Gene Clusters in Eukaryotes

Comprehensive mutant chemotyping reveals embedding of a lineage-specific biosynthetic gene cluster in wider plant metabolism

Plants produce diverse specialized metabolites with important ecological functions. It has recently become apparent that the genes for many of these pathways are not dispersed in plant genomes, but rather are arranged like beads on a string in biosynthetic gene clusters (BGCs). Pathways encoded by BGCs are as a rule dedicated linear pathways that do not form parts of wider metabolic networks. In contrast, the genes for the biosynthesis of widely distributed more ancestral metabolites such as carotenoids and anthocyanins are not clustered. Little is known about how these more recently evolved clustered pathways interact with general plant metabolism. We recently characterized a 12-gene BGC for the biosynthesis of the antimicrobial defense compound avenacin A-1, a triterpene glycoside produced by oats. Avenacin A-1 is acylated with the fluorophore N-methyl anthranilate and confers bright blue fluorescence of oat root tips under ultraviolet light. Here, we exploit a suite of >100 avenacin-deficient mutants identified by screening for reduced root fluorescence to identify genes required for the function of this paradigm BGC. Using a combination of mutant chemotyping, biochemical and molecular analysis, and genome resequencing, we identify two nonclustered genes (Sad4 and Pal2) encoding enzymes that synthesize the donors required for avenacin glycosylation and acylation (recruited from the phenylpropanoid and tryptophan pathways). Our finding of these Cluster Auxiliary Enzymes (CAEs) provides insights into the interplay between general plant metabolism and a newly evolved lineage-specific BGC.

Overexpression of the RNA-binding protein NrdA affects global gene expression and secondary metabolism in Aspergillus species

RNA-binding protein Nrd1 plays a role in RNA polymerase II transcription termination. In this study, we showed that the orthologous NrdA is important in global mRNA expression and secondary metabolism in Aspergillus species. We constructed an nrdA conditional expression strain using the Tet-On system in Aspergillus luchuenesis mut. kawachii. Downregulation of nrdA caused a severe growth defect, indicating that NrdA is essential for the proliferation of A. kawachii. Parallel RNA-sequencing and RNA immunoprecipitation-sequencing analysis identified potential NrdA-interacting transcripts, corresponding to 32% of the predicted protein-coding genes of A. kawachii. Subsequent gene ontology analysis suggested that overexpression of NrdA affects the production of secondary metabolites. To clarify this, we constructed Aspergillus nidulans, Aspergillus fumigatus, and Aspergillus oryzae strains overexpressing NrdA in the early developmental stage. Overexpression of NrdA reduced the production of sterigmatocystin and penicillin in A. nidulans, as well as that of helvolic acid and pyripyropene A in A. fumigatus. Moreover, it increased the production of kojic acid and reduced the production of penicillin in A. oryzae. These effects were accompanied by almost consistent changes in the mRNA levels of relevant genes. Collectively, these results suggest that NrdA is the essential RNA-binding protein, which plays a significant role in global gene expression and secondary metabolism in Aspergillus species.IMPORTANCENrd1, a component of the Nrd1-Nab3-Sen1 complex, is an essential RNA-binding protein involved in transcriptional termination in yeast. However, its role in filamentous fungi has not been studied. In this study, we characterized an orthologous NrdA in the Aspergillus species, identified potential NrdA-interacting mRNA, and investigated the effect of overexpression of NrdA on mRNA expression in Aspergillus luchuensis mut. kawachii. The results indicated that NrdA controls global gene expression involved in versatile metabolic pathways, including the secondary metabolic process, at least in the early developmental stage. We demonstrated that NrdA overexpression significantly affected the production of secondary metabolites in Aspergillus nidulans, Aspergillus oryzae, and Aspergillus fumigatus. Our findings are of importance to the fungal research community because the secondary metabolism is an industrially and clinically important aspect for the Aspergillus species.

The palmitoyl-CoA ligase Fum16 is part of a Fusarium verticillioides fumonisin subcluster involved in self-protection

Fusarium verticillioides produces the mycotoxin fumonisin B1 (FB1), which disrupts sphingolipid biosynthesis by inhibiting ceramide synthase and affects the health of plants, animals, and humans. The means by which F. verticillioides protects itself from its own mycotoxin are not completely understood. Some fumonisin (FUM) cluster genes do not contribute to the biosynthesis of the compound, but their function has remained enigmatic. Recently, we showed that FUM17, FUM18, and FUM19 encode two ceramide synthases and an ATP-binding cassette transporter, respectively, which play a role in antagonizing the toxicity mediated by FB1. In the present work, we uncovered functions of two adjacent genes, FUM15 and FUM16. Using homologous and heterologous expression systems, in F. verticillioides and Saccharomyces cerevisiae, respectively, we provide evidence that both contribute to protection against FB1. Our data indicate a potential role for the P450 monooxygenase Fum15 in the modification and detoxification of FB1 since the deletion and overexpression of the respective gene affected extracellular FB1 levels in both hosts. Furthermore, relative quantification of ceramide intermediates and an in vitro enzyme assay revealed that Fum16 is a functional palmitoyl-CoA ligase. It co-localizes together with the ceramide synthase Fum18 to the endoplasmic reticulum, where they contribute to sphingolipid biosynthesis. Thereby, FUM15-19 constitute a subcluster within the FUM biosynthetic gene cluster dedicated to the fungal self-protection against FB1.IMPORTANCEThe study identifies a five-gene FUM subcluster (FUM15-19) in Fusarium verticillioides involved in self-protection against FB1. FUM16 (palmitoyl-CoA ligase), FUM17, and FUM18 (ceramide synthases) enzymatically supplement ceramide biosynthesis, while FUM19 (ATP-binding cassette transporter) acts as a repressor of the FUM cluster. The evolutionary conservation of FUM15 (P450 monooxygenase) in Fusarium and Aspergillus FUM clusters is discussed, and its effect on extracellular FB1 levels in both native (F. verticillioides) and heterologous (Saccharomyces cerevisiae) hosts is highlighted. These findings enhance our understanding of mycotoxin self-protection mechanisms and could inform strategies for predicting biological activity of unknown secondary metabolites, managing mycotoxin contamination, and developing resistant crop cultivars.

Molecular fingerprints of cell size sensing and mating type differentiation in pennate diatoms

A unique cell size-sensing mechanism is at the heart of the life cycle of diatoms. During population growth, cell size decreases until a sexual size threshold (SST) is reached, below which cells become sexually competent. In most pennate diatoms, the two mating types undergo biochemical and behavioral differentiation below the SST, although the molecular pathways underlying their size-dependent maturation remain unknown. Here, we developed a method to shorten the generation time of Cylindrotheca closterium through single-cell microsurgery, enabling the transcriptomic comparison of genetically identical large and undifferentiated cells with small, sexually competent cells for six different genotypes. We identified 21 genes upregulated in small cells regardless of their mating type, revealing how cells undergo specific transcriptional reprogramming when passing the SST. Furthermore, we revealed a size-regulated gene cluster with three mating type-specific genes susceptible to sex-inducing pheromones. In addition, comparative transcriptomics confirmed the shared mating type specificity of Mating-type Related Minus 2 homologs in three pennate diatoms, suggesting them to be part of a conserved partner recognition mechanism. This study sheds light on how diatoms acquire sexual competence in a strictly size-dependent manner, revealing a complex machinery underlying size-dependent maturation, mating behavior, and heterothally in pennate diatoms.

The giant genome of lily provides insights into the hybridization of cultivated lilies

Lilies are economically important monocots known for their ornamental flowers, bulbs, and large genomes. The absence of their genomic information has impeded evolutionary studies and genome-based breeding efforts. Here, we present reference genomes for Lilium sargentiae (lily, 35.66 Gb) and Gloriosa superba (flame lily, 5.09 Gb). The giant lily genome is shaped by recent long terminal repeat retroelements. Phylogenetic analysis reveals diverse, independent origins of lily cultivars. Gene families involved in sucrose and starch metabolism are significantly expanded in the lily genome. Key homologs of XTH22, SOC1, and AP1/FUL-like genes regulate the development, bud growth transition, and floral bud growth transition of lily bulbs. Colchicine biosynthetic gene clusters are identified in G. superba but are absent in L. sargentiae, highlighting independent colchicine evolution in Colchicaceae. These genomic insights enhance understanding of Liliales evolution, providing a foundation for future breeding and molecular research.

Genomic and cell-specific regulation of benzylisoquinoline alkaloid biosynthesis in opium poppy

Opium poppy is a crop of great commercial value as a source of several opium alkaloids for the pharmaceutical industries including morphine, codeine, thebaine, noscapine, and papaverine. Most enzymes involved in benzylisoquinoline alkaloid (BIA) biosynthesis in opium poppy have been functionally characterized, and opium poppy currently serves as a model system to study BIA metabolism in plants. BIA biosynthesis in opium poppy involves two biosynthetic gene clusters associated respectively with the morphine and noscapine branches. Recent reports have shown that genes in the same cluster are co-expressed, suggesting they might also be co-regulated. However, the transcriptional regulation of opium poppy BIA biosynthesis is not well studied. Opium poppy BIA biosynthesis involves three cell types associated with the phloem system: companion cells, sieve elements, and laticifers. The transcripts and enzymes associated with BIA biosynthesis are distributed across cell types, requiring the translocation of key enzymes and pathway intermediates between cell types. Together, these suggest that the regulation of BIA biosynthesis in opium poppy is multilayered and complex, involving biochemical, genomic, and physiological mechanisms. In this review, we highlight recent advances in genome sequencing and single cell and spatial transcriptomics with a focus on how these efforts can improve our understanding of the genomic and cell-specific regulation of BIA biosynthesis. Such knowledge is vital for opium poppy genetic improvement and metabolic engineering efforts targeting the modulation of alkaloid yield and composition.

Golden Gate Assembly of Transcriptional Unit Libraries into a Rearrangeable Gene Cluster

Both regulatory sequences and genome organization contribute to the production of diverse transcript isoforms, which can influence how genes, or sets of genes, are expressed. An efficient, modular approach is needed to generate the combinatorial complexity required to empirically test many combinations of different regulatory sequences and different gene orders. Golden Gate assembly provides such a tool for seamless one-pot cleavage and ligation, by using type IIS restriction enzymes, which cleave outside of their recognition site. In addition to reducing the number of steps, this one-pot reaction can improve correct assemblies by the continued cleavage of self-ligation products that retain the recognition site. Switching the specific restriction enzyme used between steps allows for modular assembly of several units. A protocol to perform modular assemblies with two type IIS restriction enzymes, namely BsaI-v2-HF and BsmBI-v2, is described here. This protocol includes a description for generating destination vectors that add loxPsym sites between transcriptional units, allowing for diversification of gene order, orientation, and spacing.

Post-genomic illumination of paclitaxel biosynthesis

Paclitaxel rapidly became one of the most effective anticancer drugs. However, the production of paclitaxel is hindered by substantial challenges, particularly considering the significant quantities of drug required and the inherently low concentration of paclitaxel and its intermediates in plants. Paclitaxel is currently produced in a so-called semi-synthesis in which baccatin III is extracted from Taxus species and chemically converted to paclitaxel. Despite the fact that many of the intermediates of paclitaxel biosynthesis are yet to be experimentally determined, a set of recent papers-facilitated by the sequencing and assembly of three Taxus genomes-has uncovered the minimal gene sets for both baccatin III and paclitaxel biosynthesis. Here we summarize the key milestones towards our understanding of paclitaxel biosynthesis and highlight recent advancements made possible by genome-level analysis of potential key genes involved. We argue that these studies will ultimately pave the way towards the elucidation of the entire paclitaxel biosynthetic pathway and facilitate the industrial production of paclitaxel via synthetic biology approaches. However, several major challenges lie ahead before we can fully tap into the amazing curative potential that taxanes provide.

Current and future perspectives for enhancing our understanding of the evolution of plant metabolism

The special issue 'The evolution of plant metabolism' has brought together original research, reviews and opinions that cover various aspects from the full breath of plant metabolism including its interaction with the environment including other species. Here, we briefly summarize these efforts and attempts to extract a consensus opinion of the best manner in which to tackle this subject both now and in the future. This article is part of the theme issue 'The evolution of plant metabolism'.

Understanding metabolic diversification in plants: branchpoints in the evolution of specialized metabolism

Plants are chemical engineers par excellence. Collectively they make a vast array of structurally diverse specialized metabolites. The raw materials for building new pathways (genes encoding biosynthetic enzymes) are commonly recruited directly or indirectly from primary metabolism. Little is known about how new metabolic pathways and networks evolve in plants, or what key nodes contribute to branches that lead to the biosynthesis of diverse chemicals. Here we review the molecular mechanisms underlying the generation of biosynthetic branchpoints. We also consider examples in which new metabolites are formed through the joining of precursor molecules arising from different biosynthetic routes, a scenario that greatly increases both the diversity and complexity of specialized metabolism. Given the emerging importance of metabolic gene clustering in helping to identify new enzymes and pathways, we further cover the significance of biosynthetic gene clusters in relation to metabolic networks and dedicated biosynthetic pathways. In conclusion, an improved understanding of the branchpoints between metabolic pathways will be key in order to be able to predict and illustrate the complex structure of metabolic networks and to better understand the plasticity of plant metabolism. This article is part of the theme issue 'The evolution of plant metabolism'.

Review on computer-assisted biosynthetic capacities elucidation to assess metabolic interactions and communication within microbial communities

Microbial communities thrive through interactions and communication, which are challenging to study as most microorganisms are not cultivable. To address this challenge, researchers focus on the extracellular space where communication events occur. Exometabolomics and interactome analysis provide insights into the molecules involved in communication and the dynamics of their interactions. Advances in sequencing technologies and computational methods enable the reconstruction of taxonomic and functional profiles of microbial communities using high-throughput multi-omics data. Network-based approaches, including community flux balance analysis, aim to model molecular interactions within and between communities. Despite these advances, challenges remain in computer-assisted biosynthetic capacities elucidation, requiring continued innovation and collaboration among diverse scientists. This review provides insights into the current state and future directions of computer-assisted biosynthetic capacities elucidation in studying microbial communities.

Variation in a Poaceae-conserved fatty acid metabolic gene cluster controls rice yield by regulating male fertility

A wide variety of metabolic gene clusters exist in eukaryotic genomes, but fatty acid metabolic gene clusters have not been discovered. Here, combining with metabolic and phenotypic genome-wide association studies, we identify a major locus containing a six-gene fatty acid metabolic gene cluster on chromosome 3 (FGC3) that controls the cutin monomer hydroxymonoacylglycerols (HMGs) contents and rice yield, possibly through variation in the transcription of FGC3 members. We show that HMGs are sequentially synthesized in the endoplasmic reticulum by OsFAR2, OsKCS11, OsGPAT6, OsCYP704B2 and subsequently transported to the apoplast by OsABCG22 and OsLTPL82. Mutation of FGC3 members reduces HMGs, leading to defective male reproductive development and a significant decrease in yield. OsMADS6 and OsMADS17 directly regulate FGC3 and thus influence male reproduction and yield. FGC3 is conserved in Poaceae and likely formed prior to the divergence of Pharus latifolius. The eukaryotic fatty acid and plant primary metabolic gene cluster we identified show a significant impact on the origin and evolution of Poaceae and has potential for application in hybrid crop breeding.

Hordedane diterpenoid phytoalexins restrict Fusarium graminearum infection but enhance Bipolaris sorokiniana colonization of barley roots

Plant immunity is a multilayered process that includes recognition of patterns or effectors from pathogens to elicit defense responses. These include the induction of a cocktail of defense metabolites that typically restrict pathogen virulence. Here, we investigate the interaction between barley roots and the fungal pathogens Bipolaris sorokiniana (Bs) and Fusarium graminearum (Fg) at the metabolite level. We identify hordedanes, a previously undescribed set of labdane-related diterpenoids with antimicrobial properties, as critical players in these interactions. Infection of barley roots by Bs and Fg elicits hordedane synthesis from a 600-kb gene cluster. Heterologous reconstruction of the biosynthesis pathway in yeast and Nicotiana benthamiana produced several hordedanes, including one of the most functionally decorated products 19-β-hydroxy-hordetrienoic acid (19-OH-HTA). Barley mutants in the diterpene synthase genes of this cluster are unable to produce hordedanes but, unexpectedly, show reduced Bs colonization. By contrast, colonization by Fusarium graminearum, another fungal pathogen of barley and wheat, is 4-fold higher in the mutants completely lacking hordedanes. Accordingly, 19-OH-HTA enhances both germination and growth of Bs, whereas it inhibits other pathogenic fungi, including Fg. Analysis of microscopy and transcriptomics data suggest that hordedanes delay the necrotrophic phase of Bs. Taken together, these results show that adapted pathogens such as Bs can subvert plant metabolic defenses to facilitate root colonization.

Photosynthetic capacity and assimilate transport of the lower canopy influence maize yield under high planting density

Photosynthesis is a major trait of interest for the development of high-yield crop plants. However, little is known about the effects of high-density planting on photosynthetic responses at the whole-canopy level. Using the high-yielding maize (Zea mays L.) cultivars "LY66," "MC670," and "JK968," we conducted a 2-yr field experiment to assess ear development in addition to leaf characteristics and photosynthetic parameters in each canopy layer at 4 planting densities. Increased planting density promoted high grain yield and population-scale biomass accumulation despite reduced per-plant productivity. MC670 had the strongest adaptability to high-density planting conditions. A physiological analysis showed that increased planting density primarily led to decreases in the single-leaf area above the ear for LY66 and MC670 and below the ear for JK968. Furthermore, high planting density decreased chlorophyll content and the photosynthetic rate due to decreased canopy transmission, leading to severe decreases in single-plant biomass accumulation in the lower canopy. Moreover, increased planting density improved presilking biomass transfer, especially in the lower canopy. The yield showed significant positive relationships with photosynthesis and biomass in the lower canopy, demonstrating the important contributions of these leaves to grain yield under dense planting conditions. Increased planting density led to retarded ear development as a consequence of reduced glucose and fructose contents in the ears, indicating reductions in sugar transport that were associated with limited sink organ development, reduced kernel number, and yield loss. Overall, these findings highlighted the photosynthetic capacities of the lower canopy as promising targets for improving maize yield under dense planting conditions.

Anthocyanin gene enrichment in the distal region of cotton chromosome A07: mechanisms of reproductive organ coloration

Introduction

The biosynthesis of secondary metabolites like anthocyanins is often governed by metabolic gene clusters (MGCs) in the plant ancestral genome. However, the existence of gene clusters specifically regulating anthocyanin accumulation in certain organs is not well understood.

Methods and results

In this study, we identify MGCs linked to the coloration of cotton reproductive organs, such as petals, spots, and fibers. Through genetic analysis and map-based cloning, we pinpointed key genes on chromosome A07, such as PCC/GhTT19, which is involved in anthocyanin transport, and GbBM and GhTT2-3A, which are associated with the regulation of anthocyanin and proanthocyanidin biosynthesis. Our results demonstrate the coordinated control of anthocyanin and proanthocyanidin pathways, highlighting the evolutionary significance of MGCs in plant adaptation. The conservation of these clusters in cotton chromosome A07 across species underscores their importance in reproductive development and color variation. Our study sheds light on the complex biosynthesis and transport mechanisms for plant pigments, emphasizing the role of transcription factors and transport proteins in pigment accumulation.

Discussion

This research offers insights into the genetic basis of color variation in cotton reproductive organs and the potential of MGCs to enhance our comprehension of plant secondary metabolism.

Noncontiguous operon atlas for the Staphylococcus aureus genome

Bacteria synchronize the expression of genes with related functions by organizing genes into operons so that they are cotranscribed together in a single polycistronic messenger RNA. However, some cellular processes may benefit if the simultaneous production of the operon proteins coincides with the inhibition of the expression of an antagonist gene. To coordinate such situations, bacteria have evolved noncontiguous operons (NcOs), a subtype of operons that contain one or more genes that are transcribed in the opposite direction to the other operon genes. This structure results in overlapping transcripts whose expression is mutually repressed. The presence of NcOs cannot be predicted computationally and their identification requires a detailed knowledge of the bacterial transcriptome. In this study, we used direct RNA sequencing methodology to determine the NcOs map in the Staphylococcus aureus genome. We detected the presence of 18 NcOs in the genome of S. aureus and four in the genome of the lysogenic prophage 80α. The identified NcOs comprise genes involved in energy metabolism, metal acquisition and transport, toxin-antitoxin systems, and control of the phage life cycle. Using the menaquinone operon as a proof of concept, we show that disarrangement of the NcO architecture results in a reduction of bacterial fitness due to an increase in menaquinone levels and a decrease in the rate of oxygen consumption. Our study demonstrates the significance of NcO structures in bacterial physiology and emphasizes the importance of combining operon maps with transcriptomic data to uncover previously unnoticed functional relationships between neighbouring genes.

Full-length transcriptome analysis of Ophioglossum vulgatum: effects of experimentally identified chloroplast gene clusters on expression and evolutionary patterns

Genes with similar or related functions in chloroplasts are often arranged in close proximity, forming clusters on chromosomes. These clusters are transcribed coordinated to facilitate the expression of genes with specific function. Our previous study revealed a significant negative correlation between the chloroplast gene expression level of the rare medicinal fern Ophioglossum vulgatum and its evolutionary rates as well as selection pressure. Therefore, in this study, we employed a combination of SMRT and Illumina sequencing technology to analyze the full-length transcriptome sequencing of O. vulgatum for the first time. In particular, we experimentally identified gene clusters based on transcriptome data and investigated the effects of chloroplast gene clustering on expression and evolutionary patterns. The results revealed that the total sequenced data volume of the full-length transcriptome of O. vulgatum amounted to 71,950,652,163 bp, and 110 chloroplast genes received transcript coverage. Nine different types of gene clusters were experimentally identified in their transcripts. The chloroplast cluster genes may cause a decrease in non-synonymous substitution rate and selection pressure, as well as a reduction in transversion rate, transition rate, and their ratio. While expression levels of chloroplast cluster genes in leaf, sporangium, and stem would be relatively elevated. The Mann-Whitney U test indicated statistically significant in the selection pressure, sporangia and leaves groups (P < 0.05). We have contributed novel full-length transcriptome data resources for ferns, presenting new evidence on the effects of chloroplast gene clustering on expression land evolutionary patterns, and offering new theoretical support for transgenic research through gene clustering.

Wheat stripe rust resistance locus YR63 is a hot spot for evolution of defence genes - a pangenome discovery

BACKGROUND

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), poses a threat to global wheat production. Deployment of widely effective resistance genes underpins management of this ongoing threat. This study focused on the mapping of stripe rust resistance gene YR63 from a Portuguese hexaploid wheat landrace AUS27955 of the Watkins Collection.

RESULTS

YR63 exhibits resistance to a broad spectrum of Pst races from Australia, Africa, Asia, Europe, Middle East and South America. It was mapped to the short arm of chromosome 7B, between two single nucleotide polymorphic (SNP) markers sunCS_YR63 and sunCS_67, positioned at 0.8 and 3.7 Mb, respectively, in the Chinese Spring genome assembly v2.1. We characterised YR63 locus using an integrated approach engaging targeted genotyping-by-sequencing (tGBS), mutagenesis, resistance gene enrichment and sequencing (MutRenSeq), RNA sequencing (RNASeq) and comparative genomic analysis with tetraploid (Zavitan and Svevo) and hexaploid (Chinese Spring) wheat genome references and 10+ hexaploid wheat genomes. YR63 is positioned at a hot spot enriched with multiple nucleotide-binding and leucine rich repeat (NLR) and kinase domain encoding genes, known widely for defence against pests and diseases in plants and animals. Detection of YR63 within these gene clusters is not possible through short-read sequencing due to high homology between members. However, using the sequence of a NLR member we were successful in detecting a closely linked SNP marker for YR63 and validated on a panel of Australian bread wheat, durum and triticale cultivars.

CONCLUSIONS

This study highlights YR63 as a valuable source for resistance against Pst in Australia and elsewhere. The closely linked SNP marker will facilitate rapid introgression of YR63 into elite cultivars through marker-assisted selection. The bottleneck of this study reinforces the necessity for a long-read sequencing such as PacBio or Oxford Nanopore based techniques for accurate detection of the underlying resistance gene when it is part of a large gene cluster.

Elucidation of the melitidin biosynthesis pathway in pummelo

Specialized plant metabolism is a rich resource of compounds for drug discovery. The acylated flavonoid glycoside melitidin is being developed as an anti-cholesterol statin drug candidate, but its biosynthetic route in plants has not yet been fully characterized. Here, we describe the gene discovery and functional characterization of a new flavonoid gene cluster (UDP-glucuronosyltransferases (CgUGTs), 1,2 rhamnosyltransferase (Cg1,2RhaT), acyltransferases (CgATs)) that is responsible for melitidin biosynthesis in pummelo (Citrus grandis (L.) Osbeck). Population variation analysis indicated that the tailoring of acyltransferases, specific for bitter substrates, mainly determine the natural abundance of melitidin. Moreover, 3-hydroxy-3-methylglutaryl-CoA reductase enzyme inhibition assays showed that the product from this metabolic gene cluster, melitidin, may be an effective anti-cholesterol statin drug candidate. Co-expression of these clustered genes in Nicotiana benthamiana resulted in the formation of melitidin, demonstrating the potential for metabolic engineering of melitidin in a heterologous plant system. This study establishes a biosynthetic pathway for melitidin, which provides genetic resources for the breeding and genetic improvement of pummelo aimed at fortifying the content of biologically active metabolites.

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.

Functional Crosstalk between Discrete Indole Terpenoid Gene Clusters in Tolypocladium album

Indole terpenoids make up a large group of secondary metabolites that display an enticing array of bioactivities. While indole diterpene (IDT) and rarely indole sesquiterpene (IST) pathways have been found individually in filamentous fungi, here we show that both cluster types are encoded within the genome of Tolypocladium album. Through heterologous reconstruction, we demonstrate the SES cluster encodes for IST biosynthesis and can tailor IDT substrates produced by the TER cluster.

Characterization of the horse chestnut genome reveals the evolution of aescin and aesculin biosynthesis

Horse chestnut (Aesculus chinensis) is an important medicinal tree that contains various bioactive compounds, such as aescin, barrigenol-type triterpenoid saponins (BAT), and aesculin, a glycosylated coumarin. Herein, we report a 470.02 Mb genome assembly and characterize an Aesculus-specific whole-genome duplication event, which leads to the formation and duplication of two triterpenoid biosynthesis-related gene clusters (BGCs). We also show that AcOCS6, AcCYP716A278, AcCYP716A275, and AcCSL1 genes within these two BGCs along with a seed-specific expressed AcBAHD6 are responsible for the formation of aescin. Furthermore, we identify seven Aesculus-originated coumarin glycoside biosynthetic genes and achieve the de novo synthesis of aesculin in E. coli. Collinearity analysis shows that the collinear BGC segments can be traced back to early-diverging angiosperms, and the essential gene-encoding enzymes necessary for BAT biosynthesis are recruited before the splitting of Aesculus, Acer, and Xanthoceras. These findings provide insight on the evolution of gene clusters associated with medicinal tree metabolites.

Unraveling the evolutionary dynamics of the TPS gene family in land plants

Terpenes and terpenoids are key natural compounds for plant defense, development, and composition of plant oil. The synthesis and accumulation of a myriad of volatile terpenoid compounds in these plants may dramatically alter the quality and flavor of the oils, which provide great commercial utilization value for oil-producing plants. Terpene synthases (TPSs) are important enzymes responsible for terpenic diversity. Investigating the differentiation of the TPS gene family could provide valuable theoretical support for the genetic improvement of oil-producing plants. While the origin and function of TPS genes have been extensively studied, the exact origin of the initial gene fusion event - it occurred in plants or microbes - remains uncertain. Furthermore, a comprehensive exploration of the TPS gene differentiation is still pending. Here, phylogenetic analysis revealed that the fusion of the TPS gene likely occurred in the ancestor of land plants, following the acquisition of individual C- and N- terminal domains. Potential mutual transfer of TPS genes was observed among microbes and plants. Gene synteny analysis disclosed a differential divergence pattern between TPS-c and TPS-e/f subfamilies involved in primary metabolism and those (TPS-a/b/d/g/h subfamilies) crucial for secondary metabolites. Biosynthetic gene clusters (BGCs) analysis suggested a correlation between lineage divergence and potential natural selection in structuring terpene diversities. This study provides fresh perspectives on the origin and evolution of the TPS gene family.

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.

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.

A fungal sesquiterpene biosynthesis gene cluster critical for mutualist-pathogen transition in Colletotrichum tofieldiae

Plant-associated fungi show diverse lifestyles from pathogenic to mutualistic to the host; however, the principles and mechanisms through which they shift the lifestyles require elucidation. The root fungus Colletotrichum tofieldiae (Ct) promotes Arabidopsis thaliana growth under phosphate limiting conditions. Here we describe a Ct strain, designated Ct3, that severely inhibits plant growth. Ct3 pathogenesis occurs through activation of host abscisic acid pathways via a fungal secondary metabolism gene cluster related to the biosynthesis of sesquiterpene metabolites, including botrydial. Cluster activation during root infection suppresses host nutrient uptake-related genes and changes mineral contents, suggesting a role in manipulating host nutrition state. Conversely, disruption or environmental suppression of the cluster renders Ct3 beneficial for plant growth, in a manner dependent on host phosphate starvation response regulators. Our findings indicate that a fungal metabolism cluster provides a means by which infectious fungi modulate lifestyles along the parasitic-mutualistic continuum in fluctuating environments.

3D organization of regulatory elements for transcriptional regulation in Arabidopsis

BACKGROUND

Although spatial organization of compartments and topologically associating domains at large scale is relatively well studied, the spatial organization of regulatory elements at fine scale is poorly understood in plants.

RESULTS

Here we perform high-resolution chromatin interaction analysis using paired-end tag sequencing approach. We map chromatin interactions tethered with RNA polymerase II and associated with heterochromatic, transcriptionally active, and Polycomb-repressive histone modifications in Arabidopsis. Analysis of the regulatory repertoire shows that distal active cis-regulatory elements are linked to their target genes through long-range chromatin interactions with increased expression of the target genes, while poised cis-regulatory elements are linked to their target genes through long-range chromatin interactions with depressed expression of the target genes. Furthermore, we demonstrate that transcription factor MYC2 is critical for chromatin spatial organization, and propose that MYC2 occupancy and MYC2-mediated chromatin interactions coordinately facilitate transcription within the framework of 3D chromatin architecture. Analysis of functionally related gene-defined chromatin connectivity networks reveals that genes implicated in flowering-time control are functionally compartmentalized into separate subdomains via their spatial activity in the leaf or shoot apical meristem, linking active mark- or Polycomb-repressive mark-associated chromatin conformation to coordinated gene expression.

CONCLUSION

The results reveal that the regulation of gene transcription in Arabidopsis is not only by linear juxtaposition, but also by long-range chromatin interactions. Our study uncovers the fine scale genome organization of Arabidopsis and the potential roles of such organization in orchestrating transcription and development.

A Comprehensive Self-Resistance Gene Database for Natural-Product Discovery with an Application to Marine Bacterial Genome Mining

In the world of microorganisms, the biosynthesis of natural products in secondary metabolism and the self-resistance of the host always occur together and complement each other. Identifying resistance genes from biosynthetic gene clusters (BGCs) helps us understand the self-defense mechanism and predict the biological activity of natural products synthesized by microorganisms. However, a comprehensive database of resistance genes is still lacking, which hinders natural product annotation studies in large-scale genome mining. In this study, we compiled a resistance gene database (RGDB) by scanning the four available databases: CARD, MIBiG, NCBIAMR, and UniProt. Every resistance gene in the database was annotated with resistance mechanisms and possibly involved chemical compounds, using manual annotation and transformation from the resource databases. The RGDB was applied to analyze resistance genes in 7432 BGCs in 1390 genomes from a marine microbiome project. Our calculation showed that the RGDB successfully identified resistance genes for more than half of the BGCs, suggesting that the database helps prioritize BGCs that produce biologically active natural products.

Conserved H3K27me3-associated chromatin looping mediates physical interactions of gene clusters in plants

Higher-order chromatin organization is essential for transcriptional regulation, genome stability maintenance, and other genome functions. Increasing evidence has revealed significant differences in 3D chromatin organization between plants and animals. However, the extent, pattern, and rules of chromatin organization in plants are still unclear. In this study, we systematically identified and characterized long-range chromatin loops in the Arabidopsis 3D genome. We identified hundreds of long-range cis chromatin loops and found their anchor regions are closely associated with H3K27me3 epigenetic modifications. Furthermore, we demonstrated that these chromatin loops are dependent on Polycomb group (PcG) proteins, suggesting that the Polycomb repressive complex 2 (PRC2) complex is essential for establishing and maintaining these novel loops. Although most of these PcG-medicated chromatin loops are stable, many of these loops are tissue-specific or dynamically regulated by different treatments. Interestingly, tandemly arrayed gene clusters and metabolic gene clusters are enriched in anchor regions. Long-range H3K27me3-marked chromatin interactions are associated with the coregulation of specific gene clusters. Finally, we also identified H3K27me3-associated chromatin loops associated with gene clusters in Oryza sativa and Glycine max, indicating that these long-range chromatin loops are conserved in plants. Our results provide novel insights into genome evolution and transcriptional coregulation in plants.

Endophytic fungi: hidden treasure chest of antimicrobial metabolites interrelationship of endophytes and metabolites

Endophytic fungi comprise host-associated fungal communities which thrive within the tissues of host plants and produce a diverse range of secondary metabolites with various bioactive attributes. The metabolites such as phenols, polyketides, saponins, alkaloids help to mitigate biotic and abiotic stresses, fight against pathogen attacks and enhance the plant immune system. We present an overview of the association of endophytic fungal communities with a plant host and discuss molecular mechanisms induced during their symbiotic interaction. The overview focuses on the secondary metabolites (especially those of terpenoid nature) secreted by endophytic fungi and their respective function. The recent advancement in multi-omics approaches paved the way for identification of these metabolites and their characterization via comparative analysis of extensive omics datasets. This study also elaborates on the role of diverse endophytic fungi associated with key agricultural crops and hence important for sustainability of agriculture.

Regulatory architecture of housekeeping genes is driven by promoter assemblies

Genes that are key to cell identity are generally regulated by cell-type-specific enhancer elements bound by transcription factors, some of which facilitate looping to distant gene promoters. In contrast, genes that encode housekeeping functions, whose regulation is essential for normal cell metabolism and growth, generally lack interactions with distal enhancers. We find that Ronin (Thap11) assembles multiple promoters of housekeeping and metabolic genes to regulate gene expression. This behavior is analogous to how enhancers are brought together with promoters to regulate cell identity genes. Thus, Ronin-dependent promoter assemblies provide a mechanism to explain why housekeeping genes can forgo distal enhancer elements and why Ronin is important for cellular metabolism and growth control. We propose that clustering of regulatory elements is a mechanism common to cell identity and housekeeping genes but is accomplished by different factors binding distinct control elements to establish enhancer-promoter or promoter-promoter interactions, respectively.

On the evolution of natural product biosynthesis

Natural products are the raw material for drug discovery programmes. Bioactive natural products are used extensively in medicine and agriculture and have found utility as antibiotics, immunosuppressives, anti-cancer drugs and anthelminthics. Remarkably, the natural role and what mechanisms drive evolution of these molecules is relatively poorly understood. The exponential increase in genome and chemical data in recent years, coupled with technical advances in bioinformatics and genetics have enabled progress to be made in understanding the evolution of biosynthetic gene clusters and the products of their enzymatic machinery. Here we discuss the diversity of natural products, incorporating the mechanisms that govern evolution of metabolic pathways and how this can be applied to biosynthetic gene clusters. We build on the nomenclature of natural products in terms of primary, integrated, secondary and specialised metabolism and place this within an ecology-evolutionary-developmental biology framework. This eco-evo-devo framework we believe will help to clarify the nature and use of the term specialised metabolites in the future.

Vacuolar MATE/DTX protein-mediated cucurbitacin C transport is co-regulated with bitterness biosynthesis in cucumber

Membrane-localized transporters constitute important components for specialized metabolism in plants. However, due to the vast array of specialized metabolites produced by plants, and the large families of transporter genes, knowledge about the intracellular and intercellular transport of plant metabolites is still in its infancy. Cucurbitacins are bitter and defensive triterpenoids produced mainly in the cucurbits. Using a comparative genomics and multi-omics approach, a MATE gene (CsMATE1), physically clustered with cucurbitacin C (CuC) biosynthetic genes, was identified and functionally shown to sequester CuC in cucumber leaf mesophyll cells. Notably, the CuC transport process is strictly co-regulated with CuC biosynthesis. CsMATE1 clustering with bitterness biosynthesis genes may provide benefits and a basis for this feedback regulation on CuC sequestration and biosynthesis. Identification of transport systems for plant-specialized metabolites can accelerate the metabolic engineering of high-value-added compounds by simplifying their purification process.

Functional Clustering of Metabolically Related Genes Is Conserved across Dikarya

Transcriptional regulation is vital for organismal survival, with many layers and mechanisms collaborating to balance gene expression. One layer of this regulation is genome organization, specifically the clustering of functionally related, co-expressed genes along the chromosomes. Spatial organization allows for position effects to stabilize RNA expression and balance transcription, which can be advantageous for a number of reasons, including reductions in stochastic influences between the gene products. The organization of co-regulated gene families into functional clusters occurs extensively in Ascomycota fungi. However, this is less characterized within the related Basidiomycota fungi despite the many uses and applications for the species within this clade. This review will provide insight into the prevalence, purpose, and significance of the clustering of functionally related genes across Dikarya, including foundational studies from Ascomycetes and the current state of our understanding throughout representative Basidiomycete species.

Plant (di)terpenoid evolution: from pigments to hormones and beyond

Covering: up to 2014-2022.Diterpenoid biosynthesis in plants builds on the necessary production of (E,E,E)-geranylgeranyl diphosphate (GGPP) for photosynthetic pigment production, with diterpenoid biosynthesis arising very early in land plant evolution, enabling stockpiling of the extensive arsenal of (di)terpenoid natural products currently observed in this kingdom. This review will build upon that previously published in the Annual Review of Plant Biology, with a stronger focus on enzyme structure-function relationships, as well as additional insights into the evolution of (di)terpenoid metabolism since generated.

Analysis of Arabidopsis non-reference accessions reveals high diversity of metabolic gene clusters and discovers new candidate cluster members

Metabolic gene clusters (MGCs) are groups of genes involved in a common biosynthetic pathway. They are frequently formed in dynamic chromosomal regions, which may lead to intraspecies variation and cause phenotypic diversity. We examined copy number variations (CNVs) in four Arabidopsis thaliana MGCs in over one thousand accessions with experimental and bioinformatic approaches. Tirucalladienol and marneral gene clusters showed little variation, and the latter was fixed in the population. Thalianol and especially arabidiol/baruol gene clusters displayed substantial diversity. The compact version of the thalianol gene cluster was predominant and more conserved than the noncontiguous version. In the arabidiol/baruol cluster, we found a large genomic insertion containing divergent duplicates of the CYP705A2 and BARS1 genes. The BARS1 paralog, which we named BARS2, encoded a novel oxidosqualene synthase. The expression of the entire arabidiol/baruol gene cluster was altered in the accessions with the duplication. Moreover, they presented different root growth dynamics and were associated with warmer climates compared to the reference-like accessions. In the entire genome, paired genes encoding terpene synthases and cytochrome P450 oxidases were more variable than their nonpaired counterparts. Our study highlights the role of dynamically evolving MGCs in plant adaptation and phenotypic diversity.

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.

Whole genome sequence characterization of Aspergillus terreus ATCC 20541 and genome comparison of the fungi A. terreus

Aspergillus terreus is well-known for lovastatin and itaconic acid production with biomedical and commercial importance. The mechanisms of metabolite formation have been extensively studied to improve their yield through genetic engineering. However, the combined repertoire of carbohydrate-active enzymes (CAZymes), cytochrome P450s (CYP) enzymes, and secondary metabolites (SMs) in the different A. terreus strains has not been well studied yet, especially with respect to the presence of biosynthetic gene clusters (BGCs). Here we present a 30 Mb whole genome sequence of A. terreus ATCC 20541 in which we predicted 10,410 protein-coding genes. We compared the CAZymes, CYPs enzyme, and SMs across eleven A. terreus strains, and the results indicate that all strains have rich pectin degradation enzyme and CYP52 families. The lovastatin BGC of lovI was linked with lovF in A. terreus ATCC 20541, and the phenomenon was not found in the other strains. A. terreus ATCC 20541 lacked a non-ribosomal peptide synthetase (AnaPS) participating in acetylaszonalenin production, which was a conserved protein in the ten other strains. Our results present a comprehensive analysis of CAZymes, CYPs enzyme, and SM diversities in A. terreus strains and will facilitate further research in the function of BGCs associated with valuable SMs.

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.

Transperons: RNA operons as effectors of coordinated gene expression in eukaryotes

Coordinated gene expression allows spatiotemporal control of cellular processes and is achieved by the cotranscription/translation of functionally related genes/proteins. Prokaryotes evolved polycistronic messages (operons) to confer expression from a single promoter to efficiently cotranslate proteins functioning on the same pathway. Yet, despite having far greater diversity (e.g., gene number, distribution, modes of expression), eukaryotic cells employ individual promoters and monocistronic messages. Although gene expression is modular, it does not account for how eukaryotes achieve coordinated localized translation. The RNA operon theory states that mRNAs derived from different chromosomes assemble into ribonucleoprotein particles (RNPs) that act as functional operons to generate protein cohorts upon cotranslation. Work in yeast has now validated this theory and shown that intergenic associations and noncanonical histone functions create pathway-specific RNA operons (transperons) that regulate cell physiology. Herein the involvement of chromatin organization in transperon formation and programmed gene coexpression is discussed.

Predicting Cloned Disease Resistance Gene Homologs (CDRHs) in Radish, Underutilised Oilseeds, and Wild Brassicaceae Species

Brassicaceae crops, including Brassica, Camelina and Raphanus species, are among the most economically important crops globally; however, their production is affected by several diseases. To predict cloned disease resistance (R) gene homologs (CDRHs), we used the protein sequences of 49 cloned R genes against fungal and bacterial diseases in Brassicaceae species. In this study, using 20 Brassicaceae genomes (17 wild and 3 domesticated species), 3172 resistance gene analogs (RGAs) (2062 nucleotide binding-site leucine-rich repeats (NLRs), 497 receptor-like protein kinases (RLKs) and 613 receptor-like proteins (RLPs)) were identified. CDRH clusters were also observed in Arabis alpina, Camelina sativa and Cardamine hirsuta with assigned chromosomes, consisting of 62 homogeneous (38 NLR, 17 RLK and 7 RLP clusters) and 10 heterogeneous RGA clusters. This study highlights the prevalence of CDRHs in the wild relatives of the Brassicaceae family, which may lay the foundation for rapid identification of functional genes and genomics-assisted breeding to develop improved disease-resistant Brassicaceae crop cultivars.

The Chromosome-level genome of Aesculus wilsonii provides new insights into terpenoid biosynthesis and Aesculus evolution

Aesculus L. (buckeye and horse chestnut) are woody plant species with important horticultural and medicinal values. Aesculus seeds are widely used as biomedicine and cosmetic ingredients due to their saponins. We report a chromosomal-scale genome of Aesculus wilsonii. Sequences amounting to a total of 579.01 Mb were assembled into 20 chromosomes. More than half of the genome (54.46%) were annotated as repetitive sequences, and 46,914 protein-coding genes were predicted. In addition to the widespread gamma event with core eudicots, a unique whole-genome duplication (WGD) event (17.69 Mya) occurred in Aesculus after buckeye differentiated from longan. Due to WGD events and tandem duplications, the related synthetic genes of triterpene saponins unique to Aesculus increased significantly. Combined with transcriptome characterization, the study preliminarily resolved the biosynthetic pathway of triterpenoid saponins like aescin in A. wilsonii genome. Analyses of the resequencing of 104 buckeye accessions revealed clear relationship between the geographic distribution and genetic differentiation of buckeye trees in China. We found that the buckeye species found in southern Shaanxi is A. wilsonii rather than A. chinensis. Population dynamics analysis further suggests that the population size and evolution of existing buckeye species have been influenced by climate fluctuations during the Pleistocene and recent domestication events. The genome of A. wilsonii and population genomics of Aesculus provide a resource for future research on Hippocastanaceae. These findings will contribute to the utilization and diversity protection of Aesculus.

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.

Core circadian clock and light signaling genes brought into genetic linkage across the green lineage

The circadian clock is conserved at both the level of transcriptional networks as well as core genes in plants, ensuring that biological processes are phased to the correct time of day. In the model plant Arabidopsis (Arabidopsis thaliana), the core circadian SHAQKYF-type-MYB (sMYB) genes CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) and REVEILLE (RVE4) show genetic linkage with PSEUDO-RESPONSE REGULATOR 9 (PRR9) and PRR7, respectively. Leveraging chromosome-resolved plant genomes and syntenic ortholog analysis enabled tracing this genetic linkage back to Amborella trichopoda, a sister lineage to the angiosperm, and identifying an additional evolutionarily conserved genetic linkage in light signaling genes. The LHY/CCA1-PRR5/9, RVE4/8-PRR3/7, and PIF3-PHYA genetic linkages emerged in the bryophyte lineage and progressively moved within several genes of each other across an array of angiosperm families representing distinct whole-genome duplication and fractionation events. Soybean (Glycine max) maintained all but two genetic linkages, and expression analysis revealed the PIF3-PHYA linkage overlapping with the E4 maturity group locus was the only pair to robustly cycle with an evening phase, in contrast to the sMYB-PRR morning and midday phase. While most monocots maintain the genetic linkages, they have been lost in the economically important grasses (Poaceae), such as maize (Zea mays), where the genes have been fractionated to separate chromosomes and presence/absence variation results in the segregation of PRR7 paralogs across heterotic groups. The environmental robustness model is put forward, suggesting that evolutionarily conserved genetic linkages ensure superior microhabitat pollinator synchrony, while wide-hybrids or unlinking the genes, as seen in the grasses, result in heterosis, adaptation, and colonization of new ecological niches.

Genome sequencing of Colletotrichum gloeosporioides ESO026 reveals plausible pathway of HupA

BACKGROUND

Colletotrichum gloeosporioides ES026, isolated as an endophytic fungal strain, was found to produce the important medicinal compound HuperzineA (HupA). In a genetic context, ES026 showed potential in elucidating the biosynthetic pathway of HupA.

METHODS AND RESULTS

The ES026 strain was sequenced using de-novo Illumina sequencing methods in this study. Assembling the cleaned data resulted in 58,594,804bp, consisting of 404 scaffolds. The G + C mol % content of this genome was 52.53%. The genome progressive-alignment with other 4 Colletotrichum strains revealed that ES026 showed closer relation with 030206, SMCG1#C and Nara gc5. More than 60 putative biosynthetic clusters were predicted with the fungal version antiSMASH4.0 program. More than 33 types I polyketide-related biosynthetic gene clusters were distributed, containing PKS and PKS-NRPS (polyketide-nonribosomal peptides) hybrid gene clusters. Another 8 NRPS biosynthetic gene clusters were distributed among the genome of ES026. The prenyltransferases, probably involved in aromatic prenyl-compounds and terpenoid biosynthesis, were analyzed using bioinformatics tools like MEGA.

CONCLUSION

We predicted a new possible biosynthetic pathway for the HupA from the pipecolic acid, based on the published HupA biosynthesis proposed pathway, the biosynthesis and pipecolic acid-derived compounds. We hypothesize that a hybrid PKS-NRPS mega-enzyme was probably involved in the biosynthesis of HupA with the pipecolic acid, the building block of rapamycin, as a HupA precursor. The rapamycin is produced from a polyketide biosynthesis pathway, and the domain incorporating the pipecolic acid is studied.

Lateral transfers lead to the birth of momilactone biosynthetic gene clusters in grass

Momilactone A, an important plant labdane-related diterpenoid, functions as a phytoalexin against pathogens and an allelochemical against neighboring plants. The genes involved in the biosynthesis of momilactone A are found in clusters, i.e., momilactone A biosynthetic gene clusters (MABGCs), in the rice and barnyardgrass genomes. In addition, we know little about the origin and evolution of MABGCs. Here, we integrated results from comprehensive phylogeny and comparative genomic analyses of the core genes of MABGC-like clusters and MABGCs in 40 monocot plant genomes, providing convincing evidence for the birth and evolution of MABGCs in grass species. The MABGCs found in the PACMAD clade of the core grass lineage (including Panicoideae and Chloridoideae) originated from a MABGC-like cluster in Triticeae (BOP clade) via lateral gene transfer (LGT) and followed by recruitment of MAS1/2 and CYP76L1 genes. The MABGCs in Oryzoideae originated from PACMAD through another LGT event and lost CYP76L1 afterwards. The Oryza MABGC and another Oryza diterpenoid cluster c2BGC are two distinct clusters, with the latter originating from gene duplication and relocation within Oryzoideae. Further comparison of the expression patterns of the MABGC genes between rice and barnyardgrass in response to pathogen infection and allelopathy provides novel insights into the functional innovation of MABGCs in plants. Our results demonstrate LGT-mediated origination of MABGCs in grass and shed lights into the evolutionary innovation and optimization of plant biosynthetic pathways.

A gene cluster in Ginkgo biloba encodes unique multifunctional cytochrome P450s that initiate ginkgolide biosynthesis

The ginkgo tree (Ginkgo biloba) is considered a living fossil due to its 200 million year's history under morphological stasis. Its resilience is partly attributed to its unique set of specialized metabolites, in particular, ginkgolides and bilobalide, which are chemically complex terpene trilactones. Here, we use a gene cluster-guided mining approach in combination with co-expression analysis to reveal the primary steps in ginkgolide biosynthesis. We show that five multifunctional cytochrome P450s with atypical catalytic activities generate the tert-butyl group and one of the lactone rings, characteristic of all G. biloba trilactone terpenoids. The reactions include scarless C-C bond cleavage as well as carbon skeleton rearrangement (NIH shift) occurring on a previously unsuspected intermediate. The cytochrome P450s belong to CYP families that diversifies in pre-seed plants and gymnosperms, but are not preserved in angiosperms. Our work uncovers the early ginkgolide pathway and offers a glance into the biosynthesis of terpenoids of the Mesozoic Era.

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.

Genome doubling enabled the expansion of yeast vesicle traffic pathways

Vesicle budding and fusion in eukaryotes depend on a suite of protein types, such as Arfs, Rabs, coats and SNAREs. Distinct paralogs of these proteins act at distinct intracellular locations, suggesting a link between gene duplication and the expansion of vesicle traffic pathways. Genome doubling, a common source of paralogous genes in fungi, provides an ideal setting in which to explore this link. Here we trace the fates of paralog doublets derived from the 100-Ma-old hybridization event that gave rise to the whole genome duplication clade of budding yeast. We find that paralog doublets involved in specific vesicle traffic functions and pathways are convergently retained across the entire clade. Vesicle coats and adaptors involved in secretory and early-endocytic pathways are retained as doublets, at rates several-fold higher than expected by chance. Proteins involved in later endocytic steps and intra-Golgi traffic, including the entire set of multi-subunit and coiled-coil tethers, have reverted to singletons. These patterns demonstrate that selection has acted to expand and diversify the yeast vesicle traffic apparatus, across species and time.