Genomic clustering of cyanogenic glucoside biosynthetic genes aids their identification in Lotus japonicus and suggests the repeated evolution of this chemical defence pathway

Identification of a Unique Genomic Region in Sweet Chestnut (Castanea sativa Mill.) That Controls Resistance to Asian Chestnut Gall Wasp Dryocosmus kuriphilus Yasumatsu

The Asian chestnut gall wasp (ACGW) (Hymenoptera Dryocosmus kuriphilus Yasumatsu) is a severe pest of sweet chestnut (Castanea sativa Mill.) with a strong impact on growth and nut production. A comparative field trial in Central Italy, including provenances from Spain, Italy, and Greece, was screened for ACGW infestation over consecutive years. The Greek provenance Hortiatis expressed a high proportion of immune plants and was used to perform a genome-wide association study based on DNA pool sequencing (Pool-GWAS) by comparing two DNA pools from 25 susceptible and 25 resistant plants. DNA pools were sequenced with 50X coverage depth. Sequence reads were aligned to a C. mollissima reference genome and the pools were compared to identify SNPs associated with resistance. Twenty-one significant SNPs were identified and highlighted a small genomic region on pseudochromosome 3 (Chr 3), containing 12 candidate genes of three gene families: Cytochrome P450, UDP-glycosyltransferase, and Rac-like GTP-binding protein. Functional analyses revealed a putative metabolic gene cluster related to saccharide biosynthesis in the genomic regions associated with resistance that could be involved in the production of a toxic metabolite against parasites. The comparison with previous genetic studies confirmed the involvement of Chr 3 in the control of resistance to ACGW.

Lineage-Specific CYP80 Expansion and Benzylisoquinoline Alkaloid Diversity in Early-Diverging Eudicots

Menispermaceae species, as early-diverging eudicots, can synthesize valuable benzylisoquinoline alkaloids (BIAs) like bisbenzylisoquinoline alkaloids (bisBIAs) and sinomenines with a wide range of structural diversity. However, the evolutionary mechanisms responsible for their chemo-diversity are not well understood. Here, a chromosome-level genome assembly of Menispermum dauricum is presented and demonstrated the occurrence of two whole genome duplication (WGD) events that are shared by Ranunculales and specific to Menispermum, providing a model for understanding chromosomal evolution in early-diverging eudicots. The biosynthetic pathway for diverse BIAs in M. dauricum is reconstructed by analyzing the transcriptome and metabolome. Additionally, five catalytic enzymes - one norcoclaurine synthase (NCS) and four cytochrome P450 monooxygenases (CYP450s) - from M. dauricum are responsible for the formation of the skeleton, hydroxylated modification, and C-O/C-C phenol coupling of BIAs. Notably, a novel leaf-specific MdCYP80G10 enzyme that catalyzes C2'-C4a phenol coupling of (S)-reticuline into sinoacutine, the enantiomer of morphinan compounds, with predictable stereospecificity is discovered. Moreover, it is found that Menispermum-specific CYP80 gene expansion, as well as tissue-specific expression, has driven BIA diversity in Menispermaceae as compared to other Ranunculales species. This study sheds light on WGD occurrences in early-diverging eudicots and the evolution of diverse BIA biosynthesis.

Metabolism of phenolic compounds catalyzed by Tomato CYP736A61

Cytochrome P450s (CYPs) regulate plant growth and stress responses by producing diverse primary and secondary metabolites. However, the function of many plant CYPs remains unknown because, despite their structural similarity, predicting the enzymatic activity of CYPs is difficult. In this study, one member of the CYP736A subfamily (CYP736A61) from tomatoes was isolated and characterized its enzymatic functions. CYP736A61 was successfully expressed in Escherichia coli through co-expression with molecular chaperones. The purified CYP736A61 showed hydroxylation activity toward 7-ethoxycoumarin, producing 7-hydroxycoumarin or 3-hydroxy 7-ethoxycoumarin. Further substrate screening revealed that dihydrochalcone and stilbene derivates (resveratrol and polydatin) are the substrates of CYP736A61. CYP736A61 also mediated the hydroxylation of resveratrol and polydatin, albeit with low activity. Importantly, CYP736A61 mediated the cleavage of resveratrol and polydatin as well as pinostilbene and pterostilbene. Interestingly, CY736A61 also converted phloretin to naringenin chalcone. These results suggest that CYP736A61 is a novel CYP enzyme with stilbene cleavage activity.

Plant Cyanogenic-Derived Metabolites and Herbivore Counter-Defences

The release of cyanide from cyanogenic precursors is the central core of the plant defences based on the cyanogenesis process. Although cyanide is formed as a coproduct of some metabolic routes, its production is mostly due to the degradation of cyanohydrins originating from cyanogenic glycosides in cyanogenic plants and the 4-OH-ICN route in Brassicaceae. Cyanohydrins are then hydrolysed in a reversible reaction generating cyanide, being both, cyanohydrins and cyanide, toxic compounds with potential defensive properties against pests and pathogens. Based on the production of cyanogenic-derived molecules in response to the damage caused by herbivore infestation, in this review, we compile the actual knowledge of plant cyanogenic events in the plant-pest context. Besides the defensive potential, the mode of action, and the targets of the cyanogenic compounds to combat phytophagous insects and acari, special attention has been paid to arthropod responses and the strategies to overcome the impact of cyanogenesis. Physiological and behavioural adaptations, as well as cyanide detoxification by β-cyanoalanine synthases, rhodaneses, and cyanases are common ways of phytophagous arthropods defences against the cyanide produced by plants. Much experimental work is needed to further understand the complexities and specificities of the defence-counter-defence system to be applied in breeding programs.

Nitrile biosynthesis in nature: how and why?

Covering: up to the end of 2023Natural nitriles comprise a small set of secondary metabolites which however show intriguing chemical and functional diversity. Various patterns of nitrile biosynthesis can be seen in animals, plants, and microorganisms with the characteristics of both evolutionary divergence and convergence. These specialized compounds play important roles in nitrogen metabolism, chemical defense against herbivores, predators and pathogens, and inter- and/or intraspecies communications. Here we review the naturally occurring nitrile-forming pathways from a biochemical perspective and discuss the biological and ecological functions conferred by diversified nitrile biosyntheses in different organisms. Elucidation of the mechanisms and evolutionary trajectories of nitrile biosynthesis underpins better understandings of nitrile-related biology, chemistry, and ecology and will ultimately benefit the development of desirable nitrile-forming biocatalysts for practical applications.

Disentangling hydroxynitrile glucoside biosynthesis in a barley (Hordeum vulgare) metabolon provides access to elite malting barleys for ethyl carbamate-free whisky production

Barley produces several specialized metabolites, including five α-, β-, and γ-hydroxynitrile glucosides (HNGs). In malting barley, presence of the α-HNG epiheterodendrin gives rise to undesired formation of ethyl carbamate in the beverage production, especially after distilling. Metabolite-GWAS identified QTLs and underlying gene candidates possibly involved in the control of the relative and absolute content of HNGs, including an undescribed MATE transporter. By screening 325 genetically diverse barley accessions, we discovered three H. vulgare ssp. spontaneum (wild barley) lines with drastic changes in the relative ratios of the five HNGs. Knock-out (KO)-lines, isolated from the barley FIND-IT resource and each lacking one of the functional HNG biosynthetic genes (CYP79A12, CYP71C103, CYP71C113, CYP71U5, UGT85F22 and UGT85F23) showed unprecedented changes in HNG ratios enabling assignment of specific and mutually dependent catalytic functions to the biosynthetic enzymes involved. The highly similar relative ratios between the five HNGs found across wild and domesticated barley accessions indicate assembly of the HNG biosynthetic enzymes in a metabolon, the functional output of which was reconfigured in the absence of a single protein component. The absence or altered ratios of the five HNGs in the KO-lines did not change susceptibility to the fungal phytopathogen Pyrenophora teres causing net blotch. The study provides a deeper understanding of the organization of HNG biosynthesis in barley and identifies a novel, single gene HNG-0 line in an elite spring barley background for direct use in breeding of malting barley, eliminating HNGs as a source of ethyl carbamate formation in whisky production.

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.

The Putative GATA Transcription Factor SbGATA22 as a Novel Regulator of Dhurrin Biosynthesis

Cyanogenic glucosides are specialized metabolites produced by over 3000 species of higher plants from more than 130 families. The deployment of cyanogenic glucosides is influenced by biotic and abiotic factors in addition to being developmentally regulated, consistent with their roles in plant defense and stress mitigation. Despite their ubiquity, very little is known regarding the molecular mechanisms that regulate their biosynthesis. The biosynthetic pathway of dhurrin, the cyanogenic glucoside found in the important cereal crop sorghum (Sorghum bicolor (L.) Moench), was described over 20 years ago, and yet no direct regulator of the biosynthetic genes has been identified. To isolate regulatory proteins that bind to the promoter region of the key dhurrin biosynthetic gene of sorghum, SbCYP79A1, yeast one-hybrid screens were performed. A bait fragment containing 1204 base pairs of the SbCYP79A1 5' regulatory region was cloned upstream of a reporter gene and introduced into Saccharomyces cerevisiae. Subsequently, the yeast was transformed with library cDNA representing RNA from two different sorghum developmental stages. From these screens, we identified SbGATA22, an LLM domain B-GATA transcription factor that binds to the putative GATA transcription factor binding motifs in the SbCYP79A1 promoter region. Transient assays in Nicotiana benthamiana show that SbGATA22 localizes to the nucleus. The expression of SbGATA22, in comparison with SbCYP79A1 expression and dhurrin concentration, was analyzed over 14 days of sorghum development and in response to nitrogen application, as these conditions are known to affect dhurrin levels. Collectively, these findings suggest that SbGATA22 may act as a negative regulator of SbCYP79A1 expression and provide a preliminary insight into the molecular regulation of dhurrin biosynthesis in sorghum.

The telomere-to-telomere (T2T) genome of Peucedanum praeruptorum Dunn provides insights into the genome evolution and coumarin biosynthesis

BACKGROUND

Traditional Chinese medicine has used Peucedanum praeruptorum Dunn (Apiaceae) for a long time. Various coumarins, including the significant constituents praeruptorin (A-E), are the active constituents in the dried roots of P. praeruptorum. Previous transcriptomic and metabolomic studies have attempted to elucidate the distribution and biosynthetic network of these medicinal-valuable compounds. However, the lack of a high-quality reference genome impedes an in-depth understanding of genetic traits and thus the development of better breeding strategies.

RESULTS

A telomere-to-telomere (T2T) genome was assembled for P. praeruptorum by combining PacBio HiFi, ONT ultra-long, and Hi-C data. The final genome assembly was approximately 1.798 Gb, assigned to 11 chromosomes with genome completeness >98%. Comparative genomic analysis suggested that P. praeruptorum experienced 2 whole-genome duplication events. By the transcriptomic and metabolomic analysis of the coumarin metabolic pathway, we presented coumarins' spatial and temporal distribution and the expression patterns of critical genes for its biosynthesis. Notably, the COSY and cytochrome P450 genes showed tandem duplications on several chromosomes, which may be responsible for the high accumulation of coumarins.

CONCLUSIONS

A T2T genome for P. praeruptorum was obtained, providing molecular insights into the chromosomal distribution of the coumarin biosynthetic genes. This high-quality genome is an essential resource for designing engineering strategies for improving the production of these valuable compounds.

Biosynthesis, herbivore induction, and defensive role of phenylacetaldoxime glucoside

Aldoximes are well-known metabolic precursors for plant defense compounds such as cyanogenic glycosides, glucosinolates, and volatile nitriles. They are also defenses themselves produced in response to herbivory; however, it is unclear whether aldoximes can be stored over a longer term as defense compounds and how plants protect themselves against the potential autotoxic effects of aldoximes. Here, we show that the Neotropical myrmecophyte tococa (Tococa quadrialata, recently renamed Miconia microphysca) accumulates phenylacetaldoxime glucoside (PAOx-Glc) in response to leaf herbivory. Sequence comparison, transcriptomic analysis, and heterologous expression revealed that 2 cytochrome P450 enzymes, CYP79A206 and CYP79A207, and the UDP-glucosyltransferase UGT85A123 are involved in the formation of PAOx-Glc in tococa. Another P450, CYP71E76, was shown to convert PAOx to the volatile defense compound benzyl cyanide. The formation of PAOx-Glc and PAOx in leaves is a very local response to herbivory but does not appear to be regulated by jasmonic acid signaling. In contrast to PAOx, which was only detectable during herbivory, PAOx-Glc levels remained high for at least 3 d after insect feeding. This, together with the fact that gut protein extracts of 3 insect herbivore species exhibited hydrolytic activity toward PAOx-Glc, suggests that the glucoside is a stable storage form of a defense compound that may provide rapid protection against future herbivory. Moreover, the finding that herbivory or pathogen elicitor treatment also led to the accumulation of PAOx-Glc in 3 other phylogenetically distant plant species suggests that the formation and storage of aldoxime glucosides may represent a widespread plant defense response.

Exploring the Evolvability of Plant Specialized Metabolism: Uniqueness Out Of Uniformity and Uniqueness Behind Uniformity

The huge structural diversity exhibited by plant specialized metabolites has primarily been considered to result from the catalytic specificity of their biosynthetic enzymes. Accordingly, enzyme gene multiplication and functional differentiation through spontaneous mutations have been established as the molecular mechanisms that drive metabolic evolution. Nevertheless, how plants have assembled and maintained such metabolic enzyme genes and the typical clusters that are observed in plant genomes, as well as why identical specialized metabolites often exist in phylogenetically remote lineages, is currently only poorly explained by a concept known as convergent evolution. Here, we compile recent knowledge on the co-presence of metabolic modules that are common in the plant kingdom but have evolved under specific historical and contextual constraints defined by the physicochemical properties of each plant specialized metabolite and the genetic presets of the biosynthetic genes. Furthermore, we discuss a common manner to generate uncommon metabolites (uniqueness out of uniformity) and an uncommon manner to generate common metabolites (uniqueness behind uniformity). This review describes the emerging aspects of the evolvability of plant specialized metabolism that underlie the vast structural diversity of plant specialized metabolites in nature.

Metabolic GWAS-based dissection of genetic basis underlying nutrient quality variation and domestication of cassava storage root

BACKGROUND

Metabolites play critical roles in regulating nutritional qualities of plants, thereby influencing their consumption and human health. However, the genetic basis underlying the metabolite-based nutrient quality and domestication of root and tuber crops remain largely unknown.

RESULTS

We report a comprehensive study combining metabolic and phenotypic genome-wide association studies to dissect the genetic basis of metabolites in the storage root (SR) of cassava. We quantify 2,980 metabolic features in 299 cultivated cassava accessions. We detect 18,218 significant marker-metabolite associations via metabolic genome-wide association mapping and identify 12 candidate genes responsible for the levels of metabolites that are of potential nutritional importance. Me3GT, MeMYB4, and UGT85K4/UGT85K5, which are involved in flavone, anthocyanin, and cyanogenic glucoside metabolism, respectively, are functionally validated through in vitro enzyme assays and in vivo gene silencing analyses. We identify a cluster of cyanogenic glucoside biosynthesis genes, among which CYP79D1, CYP71E7b, and UGT85K5 are highly co-expressed and their allelic combination contributes to low linamarin content. We find MeMYB4 is responsible for variations in cyanidin 3-O-glucoside and delphinidin 3-O-rutinoside contents, thus controlling SR endothelium color. We find human selection affects quercetin 3-O-glucoside content and SR weight per plant. The candidate gene MeFLS1 is subject to selection during cassava domestication, leading to decreased quercetin 3-O-glucoside content and thus increased SR weight per plant.

CONCLUSIONS

These findings reveal the genetic basis of cassava SR metabolome variation, establish a linkage between metabolites and agronomic traits, and offer useful resources for genetically improving the nutrition of cassava and other root crops.

Resequencing of global Lotus corniculatus accessions reveals population distribution and genetic loci, associated with cyanogenic glycosides accumulation and growth traits

BACKGROUND

Lotus corniculatus is a widely distributed perennial legume whose great adaptability to different environments and resistance to barrenness make it an excellent forage and ecological restoration plant. However, its molecular genetics and genomic relationships among populations are yet to be uncovered.

RESULT

Here we report on a genomic variation map from worldwide 272 L. corniculatus accessions by genome resequencing. Our analysis suggests that L. corniculatus accessions have high genetic diversity and could be further divided into three subgroups, with the genetic diversity centers were located in Transcaucasia. Several candidate genes and SNP site associated with CNglcs content and growth traits were identified by genome-wide associated study (GWAS). A non-synonymous in LjMTR was responsible for the decreased expression of CNglcs synthesis genes and LjZCD was verified to positively regulate CNglcs synthesis gene CYP79D3. The LjZCB and an SNP in LjZCA promoter were confirmed to be involved in plant growth.

CONCLUSION

This study provided a large number of genomic resources and described genetic relationship and population structure among different accessions. Moreover, we attempt to provide insights into the molecular studies and breeding of CNglcs and growth traits in L. corniculatus.

Identification and characterization of cytochrome P450 CYP77A59 of loquat (Rhaphiolepis bibas) responsible for biosynthesis of phenylacetonitrile, a floral nitrile compound

MAIN CONCLUSION

Cytochrome P450 CYP77A59 is responsible for the biosynthesis of phenylacetonitrile in loquat flowers. Flowers of some plants emit volatile nitrile compounds, but the biosynthesis of these compounds is unclear. Loquat (Rhaphiolepis bibas) flowers emit characteristic N-containing volatiles, such as phenylacetonitrile (PAN), (E/Z)-phenylacetaldoxime (PAOx), and (2-nitroethyl)benzene (NEB). These volatiles likely play a defense role against pathogens and insects. PAN and NEB are commonly biosynthesized from L-phenylalanine via (E/Z)-PAOx. Two cytochrome P450s-CYP79D80 and "promiscuous fatty acid ω-hydroxylase" CYP94A90, which catalyze the formation of (E/Z)-PAOx from L-phenylalanine and NEB from (E/Z)-PAOx, respectively-are involved in NEB biosynthesis. However, the enzymes catalyzing the formation of PAN from (E/Z)-PAOx in loquat have not been identified. In this study, we aimed to identify candidate cytochrome P450s catalyzing PAN formation in loquat flowers. Yeast whole-cell biocatalyst assays showed that among nine candidate cytochrome P450s, CYP77A58 and CYP77A59 produced PAN from (E/Z)-PAOx. CYP77As catalyzed the dehydration of aldoximes, which is atypical of cytochrome P450; the reaction was NADPH-dependent, with an optimum temperature and pH of 40 °C and 8.0, respectively. CYP77As acted on (E/Z)-PAOx, (E/Z)-4-hydroxyphenylacetaldoxime, and (E/Z)-indole-3-acetaldoxime. Previously characterized CYP77As are known to hydroxylate fatty acids; loquat CYP77As did not act on tested fatty acids. We observed higher expression of CYP77A59 in flowers than in buds; expression of CYP77A58 was remarkably reduced in the flowers. Because the flowers, but not buds, emit PAN, CYP77A59 is likely responsible for the biosynthesis of PAN in loquat flowers. This study will help us understand the biosynthesis of floral nitrile compounds.

CRISPR-Cas9-mediated knockout of CYP79D1 and CYP79D2 in cassava attenuates toxic cyanogen production

Cassava (Manihot esculenta) is a starchy root crop that supports over a billion people in tropical and subtropical regions of the world. This staple, however, produces the neurotoxin cyanide and requires processing for safe consumption. Excessive consumption of insufficiently processed cassava, in combination with protein-poor diets, can have neurodegenerative impacts. This problem is further exacerbated by drought conditions which increase this toxin in the plant. To reduce cyanide levels in cassava, we used CRISPR-mediated mutagenesis to disrupt the cytochrome P450 genes CYP79D1 and CYP79D2 whose protein products catalyze the first step in cyanogenic glucoside biosynthesis. Knockout of both genes eliminated cyanide in leaves and storage roots of cassava accession 60444; the West African, farmer-preferred cultivar TME 419; and the improved variety TMS 91/02324. Although knockout of CYP79D2 alone resulted in significant reduction of cyanide, mutagenesis of CYP79D1 did not, indicating these paralogs have diverged in their function. The congruence of results across accessions indicates that our approach could readily be extended to other preferred or improved cultivars. This work demonstrates cassava genome editing for enhanced food safety and reduced processing burden, against the backdrop of a changing climate.

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

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

Recruitment of distinct UDP-glycosyltransferase families demonstrates dynamic evolution of chemical defense within Eucalyptus L'Hér

The economic and ecologically important genus Eucalyptus is rich in structurally diverse specialized metabolites. While some specialized metabolite classes are highly prevalent across the genus, the cyanogenic glucoside prunasin is only produced by c. 3% of species. To investigate the evolutionary mechanisms behind prunasin biosynthesis in Eucalyptus, we compared de novo assembled transcriptomes, together with online resources between cyanogenic and acyanogenic species. Identified genes were characterized in vivo and in vitro. Pathway characterization of cyanogenic Eucalyptus camphora and Eucalyptus yarraensis showed for the first time that the final glucosylation step from mandelonitrile to prunasin is catalyzed by a novel UDP-glucosyltransferase UGT87. This step is typically catalyzed by a member of the UGT85 family, including in Eucalyptus cladocalyx. The upstream conversion of phenylalanine to mandelonitrile is catalyzed by three cytochrome P450 (CYP) enzymes from the CYP79, CYP706, and CYP71 families, as previously shown. Analysis of acyanogenic Eucalyptus species revealed the loss of different ortholog prunasin biosynthetic genes. The recruitment of UGTs from different families for prunasin biosynthesis in Eucalyptus demonstrates important pathway heterogeneities and unprecedented dynamic pathway evolution of chemical defense within a single genus. Overall, this study provides relevant insights into the tremendous adaptability of these long-lived trees.

Emerging trends in genomic and epigenomic regulation of plant specialised metabolism

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

Biosynthesis of iridoid sex pheromones in aphids

Plants, animals, and microbes produce a plethora of natural products that are important for defense and communication. Most of these compounds show a phylogenetically restricted occurrence, but, in rare instances, the same natural product is biosynthesized by organisms in two different kingdoms. The monoterpene-derived iridoids, for example, have been found in more than 50 plant families but are also observed in several insect orders. The discovery of the aphid iridoid pathway, one of the longest and most chemically complex insect-derived natural product biosynthetic pathways reported to date, highlights the mechanisms underlying the convergent evolution of metabolic enzymes in insects and plants, including the recruitment of different enzyme classes to catalyze the same chemical processes.

Iridoid monoterpenes, widely distributed in plants and insects, have many ecological functions. While the biosynthesis of iridoids has been extensively studied in plants, little is known about how insects synthesize these natural products. Here, we elucidated the biosynthesis of the iridoids cis-trans-nepetalactol and cis-trans-nepetalactone in the pea aphid Acyrthosiphon pisum (Harris), where they act as sex pheromones. The exclusive production of iridoids in hind legs of sexual female aphids allowed us to identify iridoid genes by searching for genes specifically expressed in this tissue. Biochemical characterization of candidate enzymes revealed that the iridoid pathway in aphids proceeds through the same sequence of intermediates as described for plants. The six identified aphid enzymes are unrelated to their counterparts in plants, conclusively demonstrating an independent evolution of the entire iridoid pathway in plants and insects. In contrast to the plant pathway, at least three of the aphid iridoid enzymes are likely membrane bound. We demonstrated that a lipid environment facilitates the cyclization of a reactive enol intermediate to the iridoid cyclopentanoid-pyran scaffold in vitro, suggesting that membranes are an essential component of the aphid iridoid pathway. Altogether, our discovery of this complex insect metabolic pathway establishes the genetic and biochemical basis for the formation of iridoid sex pheromones in aphids, and this discovery also serves as a foundation for understanding the convergent evolution of complex metabolic pathways between kingdoms.

Plant metabolic gene clusters in the multi-omics era

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

Integrative omics approaches for biosynthetic pathway discovery in plants

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

Investigating the Metabolic Mesosulfuron-Methyl Resistance in Aegilops tauschii Coss. By Transcriptome Sequencing Combined with the Reference Genome

Aegilops tauschii Coss. is a malignant weed in wheat fields in China, its herbicide resistance has been threatening crop production. This study identified one mesosulfuron-methyl-resistant(R) population, JJMHN2018-05 (R), without target resistance mutations. To fully understand the resistance mechanism, non-target site resistance was investigated by using transcriptome sequencing combined with a reference genome. Results showed that the cytochrome P450 monooxygenase (P450) inhibitor malathion significantly increased the mesosulfuron-methyl sensitivity in R plants, and greater herbicide-induced glutathione S-transferase (GST) activity was also confirmed. Liquid chromatography with tandem mass spectrometry analysis further supported the enhanced mesosulfuron-methyl metabolism in R plants. Gene expression data analysis and qRT-PCR validation indicated that eight P450s, six GSTs, two glycosyltransferases (GTs), four peroxidases, and one aldo-keto reductase (AKRs) stably upregulated in R plants. This research demonstrates that the P450s and GSTs involved in enhanced mesosulfuron-methyl metabolism contribute to mesosulfuron-methyl resistance in A. tauschii and identifies potential contributors from metabolic enzyme families.

Hydroxynitrile lyase defends Arabidopsis against Tetranychus urticae

Plant-pest interactions involve multifaceted processes encompassing a complex crosstalk of pathways, molecules, and regulators aimed at overcoming defenses developed by each interacting organism. Among plant defensive compounds against phytophagous arthropods, cyanide-derived products are toxic molecules that directly target pest physiology. Here, we identified the Arabidopsis (Arabidopsis thaliana) gene encoding hydroxynitrile lyase (AtHNL, At5g10300) as one gene induced in response to spider mite (Tetranychus urticae) infestation. AtHNL catalyzes the reversible interconversion between cyanohydrins and derived carbonyl compounds with free cyanide. AtHNL loss- and gain-of-function Arabidopsis plants showed that specific activity of AtHNL using mandelonitrile as substrate was higher in the overexpressing lines than in wild-type (WT) and mutant lines. Concomitantly, mandelonitrile accumulated at higher levels in mutant lines than in WT plants and was significantly reduced in the AtHNL overexpressing lines. After mite infestation, mandelonitrile content increased in WT and overexpressing plants but not in mutant lines, while hydrogen cyanide (HCN) accumulated in the three infested Arabidopsis genotypes. Feeding bioassays demonstrated that the AtHNL gene participated in Arabidopsis defense against T. urticae. The reduced leaf damage detected in the AtHNL overexpressing lines reflected the mite's reduced ability to feed on leaves, which consequently restricted mite fecundity. In turn, mites upregulated TuCAS1 encoding β-cyanoalanine synthase to avoid the respiratory damage produced by HCN. This detoxification effect was functionally demonstrated by reduced mite fecundity observed when dsRNA-TuCAS-treated mites fed on WT plants and hnl1 mutant lines. These findings add more players in the Arabidopsis-T. urticae interplay to overcome mutual defenses.

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-wide transcriptome analysis and characterization of the cytochrome P450 flavonoid biosynthesis genes in pigeon pea (Cajanus cajan)

226 CcCYP450 genes were identified at the genomic level and were classified into 45 clades based on phylogenetic analysis. CcCYP75B165 gene was found that might play important roles in the biosynthesis of flavonoids in pigeon pea, and was significantly induced by methyl jasmonate (MeJA). The cytochrome P450 mono-oxygenase (CYP450) superfamily plays a key role in the flavonoid biosynthesis pathway and resists different kinds of stresses. Several CYP450 genes have been identified to be involved in the biosynthesis of crop protection agents. However, the CcCYP450 genes from pigeon pea have not been identified. Here, 226 CcCYP450 genes were identified at the genomic level by analysing the gene structure, distribution on chromosomes, gene duplication, and conserved motifs and were classified into 45 clades based on phylogenetic analysis. RNA-seq analysis revealed clear details of CcCYP450 genes that varied with time of MeJA (methyl jasmonate) induction. Among them, six CcCYP450 subfamily genes were found that might play important roles in the biosynthesis of flavonoids in pigeon pea. The overexpression of CcCYP75B165 in pigeon pea significantly induced the accumulation of genistin and downregulated the contents of cajaninstilbene acid, apigenin, isovitexin, and genistein and the expression of flavonoid synthase genes. This study provides theoretical guidance and plant genetic resources for cultivating new pigeon pea varieties with high flavonoid contents and exploring the molecular mechanisms of the biosynthesis of flavonoids under MeJA treatment.

MeJA-responsive bHLH transcription factor LjbHLH7 regulates cyanogenic glucoside biosynthesis in Lotus japonicus

Cyanogenic glucosides (CNglcs) play an important role in plant defense response; however, the mechanism of regulation of CNglc synthesis by the external environment and endogenous hormones is largely unclear. In this study, we found that jasmonates (JAs) promoted the synthesis of CNglcs by activating the expression of CNglc biosynthesis genes in Lotus japonicus. Several differentially expressed basic helix-loop-helix (bHLH) family genes related to the synthesis of CNglcs were identified by RNA-seq. LjbHLH7 can directly activate the expression of CYP79D3 gene, the first step of CNglc synthesis, by binding to the G-box sequence of its promoter. Transgenic plants overexpressing LjbHLH7 exhibited higher relative CNglc content and enhanced insect resistance compared with the wild type. Furthermore, the transcriptional activity of LjbHLH7 was suppressed by the interaction with the L. japonicus JASMONATE-ZIM DOMAIN protein LjJAZ4. Based on these results, we propose that LjbHLH7 acts as an activator and LjJAZ4 acts as a repressor of JA-induced regulation of CNglc biosynthesis in L. japonicus.

Cyanogenesis in cassava and its molecular manipulation for crop improvement

While cassava is one of the most important staple crops worldwide, it has received the least investment per capita consumption of any of the major global crops. This is in part due to cassava being a crop of subsistence farmers that is grown in countries with limited resources for crop improvement. While its starchy roots are rich in calories, they are poor in protein and other essential nutrients. In addition, they contain potentially toxic levels of cyanogenic glycosides which must be reduced to safe levels before consumption. Furthermore, cyanogens compromise the shelf life of harvested roots due to cyanide-induced inhibition of mitochondrial respiration, and associated production of reactive oxygen species that accelerate root deterioration. Over the past two decades, the genetic, biochemical, and developmental factors that control cyanogen synthesis, transport, storage, and turnover have largely been elucidated. It is now apparent that cyanogens contribute substantially to whole-plant nitrogen metabolism and protein synthesis in roots. The essential role of cyanogens in root nitrogen metabolism, however, has confounded efforts to create acyanogenic varieties. This review proposes alternative molecular approaches that integrate accelerated cyanogen turnover with nitrogen reassimilation into root protein that may offer a solution to creating a safer, more nutritious cassava crop.

A flavin-dependent monooxygenase produces nitrogenous tomato aroma volatiles using cysteine as a nitrogen source

Aroma is an important factor in consumer perception and acceptance of fresh tomatoes and involves a cocktail of several dozen compounds. Tomato fruits produce uncommon nitrogen-containing volatiles derived mainly from the amino acids leucine and phenylalanine. These volatiles have strong positive correlations with consumer liking. We show that an enzyme active in ripening tomatoes is responsible for the production of all nitrogenous volatiles in tomato fruit, at the expense of substrates derived from cysteine and volatile aldehydes. This discovery defines a cysteine-dependent route to nitrogenous volatiles in plants, prompting a reconsideration of the impact of sulfur metabolism on tomato flavor and quality.

Tomato (Solanum lycopersicum) produces a wide range of volatile chemicals during fruit ripening, generating a distinct aroma and contributing to the overall flavor. Among these volatiles are several aromatic and aliphatic nitrogen-containing compounds for which the biosynthetic pathways are not known. While nitrogenous volatiles are abundant in tomato fruit, their content in fruits of the closely related species of the tomato clade is highly variable. For example, the green-fruited species Solanum pennellii are nearly devoid, while the red-fruited species S. lycopersicum and Solanum pimpinellifolium accumulate high amounts. Using an introgression population derived from S. pennellii, we identified a locus essential for the production of all the detectable nitrogenous volatiles in tomato fruit. Silencing of the underlying gene (SlTNH1;Solyc12g013690) in transgenic plants abolished production of aliphatic and aromatic nitrogenous volatiles in ripe fruit, and metabolomic analysis of these fruit revealed the accumulation of 2-isobutyl-tetrahydrothiazolidine-4-carboxylic acid, a known conjugate of cysteine and 3-methylbutanal. Biosynthetic incorporation of stable isotope-labeled precursors into 2-isobutylthiazole and 2-phenylacetonitrile confirmed that cysteine provides the nitrogen atom for all nitrogenous volatiles in tomato fruit. Nicotiana benthamiana plants expressing SlTNH1 readily transformed synthetic 2-substituted tetrahydrothiazolidine-4-carboxylic acid substrates into a mixture of the corresponding 2-substituted oxime, nitro, and nitrile volatiles. Distinct from other known flavin-dependent monooxygenase enzymes in plants, this tetrahydrothiazolidine-4-carboxylic acid N-hydroxylase catalyzes sequential hydroxylations. Elucidation of this pathway is a major step forward in understanding and ultimately improving tomato flavor quality.

22R- but not 22S-hydroxycholesterol is recruited for diosgenin biosynthesis

Diosgenin is an important compound in the pharmaceutical industry and it is biosynthesized in several eudicot and monocot species, herein represented by fenugreek (a eudicot), and Dioscorea zingiberensis (a monocot). Formation of diosgenin can be achieved by the early C22,16-oxidations of cholesterol followed by a late C26-oxidation. This study reveals that, in both fenugreek and D. zingiberensis, the early C22,16-oxygenase(s) shows strict 22R-stereospecificity for hydroxylation of the substrates. Evidence against the recently proposed intermediacy of 16S,22S-dihydroxycholesterol in diosgenin biosynthesis was also found. Moreover, in contrast to the eudicot fenugreek, which utilizes a single multifunctional cytochrome P450 (TfCYP90B50) to perform the early C22,16-oxidations, the monocot D. zingiberensis has evolved two separate cytochrome P450 enzymes, with DzCYP90B71 being specific for the 22R-oxidation and DzCYP90G6 for the C16-oxidation. We suggest that the DzCYP90B71/DzCYP90G6 pair represent more broadly conserved catalysts for diosgenin biosynthesis in monocots.

Cyanogenesis in the Sorghum Genus: From Genotype to Phenotype

Domestication has resulted in a loss of genetic diversity in our major food crops, leading to susceptibility to biotic and abiotic stresses linked with climate change. Crop wild relatives (CWR) may provide a source of novel genes potentially important for re-gaining climate resilience. Sorghum bicolor is an important cereal crop with wild relatives that are endemic to Australia. Sorghum bicolor is cyanogenic, but the cyanogenic status of wild Sorghum species is not well known. In this study, leaves of wild species endemic in Australia are screened for the presence of the cyanogenic glucoside dhurrin. The direct measurement of dhurrin content and the potential for dhurrin-derived HCN release (HCNp) showed that all the tested Australian wild species were essentially phenotypically acyanogenic. The unexpected low dhurrin content may reflect the variable and generally nutrient-poor environments in which they are growing in nature. Genome sequencing of six CWR and PCR amplification of the CYP79A1 gene from additional species showed that a high conservation of key amino acids is required for correct protein function and dhurrin synthesis, pointing to the transcriptional regulation of the cyanogenic phenotype in wild sorghum as previously shown in elite sorghum.

A monocot-specific hydroxycinnamoylputrescine gene cluster contributes to immunity and cell death in rice

Phenolamides (PAs), a diverse group of specialized metabolites, including hydroxycinnamoylputrescine (HP), hydroxycinnamoylagmatine, and hydroxycinnamoyltryptamine, are important in plant resistance to biotic stress. However, the genes involved in the biosynthesis and modulation of PAs have not been fully elucidated. This study identified an HP biosynthetic gene cluster in rice (Oryza sativa) comprising one gene (OsODC) encoding a decarboxylase and two tandem-duplicated genes (OsPHT3 and OsPHT4) encoding putrescine hydroxycinnamoyl acyltransferases coexpressed in different tissues. OsODC catalyzes the conversion of ornithine to putrescine, which is used in HP biosynthesis involving OsPHT3 and OsPHT4. OsPHT3 or OsPHT4 overexpression causes HP accumulation and cell death and putrescine hydroxycinnamoyl acyltransferases (PHT) activity-dependent resistance against the fungal pathogen Magnaporthe oryzae. OsODC overexpression plants also confer enhanced resistance to M. oryzae. Notably, the basic leucine zipper transcription factor APIP5, a negative regulator of cell death, directly binds to the OsPHT4 promoter, repressing its transcription. Moreover, APIP5 suppression induces OsPHT4 expression and HP accumulation. Comparative genomic analysis revealed that the HP biosynthetic gene cluster is conserved in monocots. These results characterized a previously unidentified monocot-specific gene cluster that is involved in HP biosynthesis and contributes to defense and cell death in rice.

Gene cluster from plant to microbes: Their role in genome architecture, organism's development, specialized metabolism and drug discovery

Plants and microbes fulfil our daily requirements through different high-value chemicals, e.g., nutraceuticals, pharmaceuticals, cosmetics, and through varieties of fruits, crops, vegetables, and many more. Utmost care would therefore be taken for growth, development and sustainability of these important crops and medicinal plants and microbes. Homeobox genes and HOX clusters and their recently characterized expanded family members, including newly discovered homeobox, WOX gene from medicinal herb, Panax ginseng, significantly contributes in the growth and development of these organisms. On the other hand, secondary metabolites produced through secondary metabolism of plants and microbes are used as organisms defense as well as drugs/drug-like molecules for humans. Both the developmental HOX cluster and the biosynthetic gene-cluster (BGC) for secondary metabolites are organised in organisms genome. Genome mining and genomewide analysis of these clusters will definitely identify and characterize many more important molecules from unexplored plants and microbes and underexplored human microbiota and the evolution studies of these clusters will indicate their source of origin. Although genomics revolution now continues at a pace, till date only few hundred plant genome sequences are available. However, next-generation sequencing (NGS) technology now in market and may be applied even for plants with recalcitrant genomes, eventually may discover genomic potential towards production of secondary metabolites of diverse plants and micro-organisms present in the environment and microbiota. Additionally, the development of tools for genome mining e.g., antiSMASH, plantiSMASH, and more and more computational approaches that predicts hundreds of secondary metabolite BGCs will be discussed.

Regulation of dhurrin pathway gene expression during Sorghum bicolor development

Developmental and organ-specific expression of genes in dhurrin biosynthesis, bio-activation, and recycling offers dynamic metabolic responses optimizing growth and defence responses in Sorghum. Plant defence models evaluate the costs and benefits of resource investments at different stages in the life cycle. Poor understanding of the molecular regulation of defence deployment and remobilization hampers accuracy of the predictions. Cyanogenic glucosides, such as dhurrin are phytoanticipins that release hydrogen cyanide upon bio-activation. In this study, RNA-seq was used to investigate the expression of genes involved in the biosynthesis, bio-activation and recycling of dhurrin in Sorghum bicolor. Genes involved in dhurrin biosynthesis were highly expressed in all young developing vegetative tissues (leaves, leaf sheath, roots, stems), tiller buds and imbibing seeds and showed gene specific peaks of expression in leaves during diel cycles. Genes involved in dhurrin bio-activation were expressed early in organ development with organ-specific expression patterns. Genes involved in recycling were expressed at similar levels in the different organ during development, although post-floral initiation when nutrients are remobilized for grain filling, expression of GSTL1 decreased > tenfold in leaves and NITB2 increased > tenfold in stems. Results are consistent with the establishment of a pre-emptive defence in young tissues and regulated recycling related to organ senescence and increased demand for nitrogen during grain filling. This detailed characterization of the transcriptional regulation of dhurrin biosynthesis, bioactivation and remobilization genes during organ and plant development will aid elucidation of gene regulatory networks and signalling pathways that modulate gene expression and dhurrin levels. In-depth knowledge of dhurrin metabolism could improve the yield, nitrogen use efficiency and stress resilience of Sorghum.

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

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

Plant cytochrome P450 plasticity and evolution

The superfamily of cytochrome P450 (CYP) enzymes plays key roles in plant evolution and metabolic diversification. This review provides a status on the CYP landscape within green algae and land plants. The 11 conserved CYP clans known from vascular plants are all present in green algae and several green algae-specific clans are recognized. Clan 71, 72, and 85 remain the largest CYP clans and include many taxa-specific CYP (sub)families reflecting emergence of linage-specific pathways. Molecular features and dynamics of CYP plasticity and evolution are discussed and exemplified by selected biosynthetic pathways. High substrate promiscuity is commonly observed for CYPs from large families, favoring retention of gene duplicates and neofunctionalization, thus seeding acquisition of new functions. Elucidation of biosynthetic pathways producing metabolites with sporadic distribution across plant phylogeny reveals multiple examples of convergent evolution where CYPs have been independently recruited from the same or different CYP families, to adapt to similar environmental challenges or ecological niches. Sometimes only a single or a few mutations are required for functional interconversion. A compilation of functionally characterized plant CYPs is provided online through the Plant P450 Database (erda.dk/public/vgrid/PlantP450/).

Genome-wide association study of cyanogenic glycosides, proline, sugars, and pigments in Eucalyptus cladocalyx after 18 consecutive dry summers

Natural variation of cyanogenic glycosides, soluble sugars, proline, and nondestructive optical sensing of pigments (chlorophyll, flavonols, and anthocyanins) was examined in ex situ natural populations of Eucalyptus cladocalyx F. Muell. grown under dry environmental conditions in the southern Atacama Desert, Chile. After 18 consecutive dry seasons, considerable plant-to-plant phenotypic variation for all the traits was observed in the field. For example, leaf hydrogen cyanide (HCN) concentrations varied from 0 (two acyanogenic individuals) to 1.54 mg cyanide g^-1 DW. Subsequent genome-wide association study revealed associations with several genes with a known function in plants. HCN content was associated robustly with genes encoding Cytochrome P450 proteins, and with genes involved in the detoxification mechanism of HCN in cells (β-cyanoalanine synthase and cyanoalanine nitrilase). Another important finding was that sugars, proline, and pigment content were linked to genes involved in transport, biosynthesis, and/or catabolism. Estimates of genomic heritability (based on haplotypes) ranged between 0.46 and 0.84 (HCN and proline content, respectively). Proline and soluble sugars had the highest predictive ability of genomic prediction models (PA = 0.65 and PA = 0.71, respectively). PA values for HCN content and flavonols were relatively moderate, with estimates ranging from 0.44 to 0.50. These findings provide new understanding on the genetic architecture of cyanogenic capacity, and other key complex traits in cyanogenic E. cladocalyx.

Comparison of glucosinolate diversity in the crucifer tribe Cardamineae and the remaining order Brassicales highlights repetitive evolutionary loss and gain of biosynthetic steps

We review glucosinolate (GSL) diversity and analyze phylogeny in the crucifer tribe Cardamineae as well as selected species from Brassicaceae (tribe Brassiceae) and Resedaceae. Some GSLs occur widely, while there is a scattered distribution of many less common GSLs, tentatively sorted into three classes: ancient, intermediate and more recently evolved. The number of conclusively identified GSLs in the tribe (53 GSLs) constitute 60% of all GSLs known with certainty from any plant (89 GSLs) and apparently unique GSLs in the tribe constitute 10 of those GSLs conclusively identified (19%). Intraspecific, qualitative GSL polymorphism is known from at least four species in the tribe. The most ancient GSL biosynthesis in Brassicales probably involved biosynthesis from Phe, Val, Leu, Ile and possibly Trp, and hydroxylation at the β-position. From a broad comparison of families in Brassicales and tribes in Brassicaceae, we estimate that a common ancestor of the tribe Cardamineae and the family Brassicaceae exhibited GSL biosynthesis from Phe, Val, Ile, Leu, possibly Tyr, Trp and homoPhe (ancient GSLs), as well as homologs of Met and possibly homoIle (intermediate age GSLs). From the comparison of phylogeny and GSL diversity, we also suggest that hydroxylation and subsequent methylation of indole GSLs and usual modifications of Met-derived GSLs (formation of sulfinyls, sulfonyls and alkenyls) occur due to conserved biochemical mechanisms and was present in a common ancestor of the family. Apparent loss of homologs of Met as biosynthetic precursors was deduced in the entire genus Barbarea and was frequent in Cardamine (e.g. C. pratensis, C. diphylla, C. concatenata, possibly C. amara). The loss was often associated with appearance of significant levels of unique or rare GSLs as well as recapitulation of ancient types of GSLs. Biosynthetic traits interpreted as de novo evolution included hydroxylation at rare positions, acylation at the thioglucose and use of dihomoIle and possibly homoIle as biosynthetic precursors. Biochemical aspects of the deduced evolution are discussed and testable hypotheses proposed. Biosyntheses from Val, Leu, Ile, Phe, Trp, homoPhe and homologs of Met are increasingly well understood, while GSL biosynthesis from mono- and dihomoIle is poorly understood. Overall, interpretation of known diversity suggests that evolution of GSL biosynthesis often seems to recapitulate ancient biosynthesis. In contrast, unprecedented GSL biosynthetic innovation seems to be rare.

Influence of domestication on specialized metabolic pathways in fruit crops

During the process of plant domestication, the selection and traditional breeding for desired characters such as flavor, juiciness and nutritional value of fruits, probably have resulted in gain or loss of specialized metabolites contributing to these traits. Their appearance in fruits is likely due to the acquisition of novel and specialized metabolic pathways and their regulation, driven by systematic molecular evolutionary events facilitated by traditional breeding. Plants change their armory of specialized metabolism to adapt and survive in diverse ecosystems. This may occur through molecular evolutionary events, such as single nucleotide polymorphism, gene duplication and transposition, leading to convergent or divergent evolution of biosynthetic pathways producing such specialized metabolites. Breeding and selection for improved specific and desired traits (fruit size, color, taste, flavor, etc.) in fruit crops through conventional breeding approaches may further alter content and profile of specialized metabolites. Biosynthetic routes of these metabolites have been studied in various plants. Here, we explore the influence of plant domestication and breeding processes on the selection of biosynthetic pathways of favorable specialized metabolites in fruit crops. An orderly clustered arrangement of genes associated with their production is observed in many fruit crops. We further analyzed selection-based acquisition of specialized metabolic pathways comparing first the metabolic profiles and genes involved in their biosynthesis, followed by the genomic organization of such genes between wild and domesticated horticultural crops. Domestication of crop plants favored the acquisition and retention of metabolic pathways that enhanced the fruit value while eliminated those which produced toxic or unfavorable metabolites. Interestingly, unintentional reorganization of complex metabolic pathways by selection and traditional breeding processes has endowed us with flavorful, juicy and nutritionally rich fruits.

Large-scale genome-wide association study, using historical data, identifies conserved genetic architecture of cyanogenic glucoside content in cassava (Manihot esculenta Crantz) root

Manihot esculenta (cassava) is a root crop originating from South America that is a major staple in the tropics, including in marginal environments. This study focused on South American and African germplasm and investigated the genetic architecture of hydrogen cyanide (HCN), a major component of root quality. HCN, representing total cyanogenic glucosides, is a plant defense component against herbivory but is also toxic for human consumption. We genotyped 3354 landraces and modern breeding lines originating from 26 Brazilian states and 1389 individuals were phenotypically characterized across multi-year trials for HCN. All plant material was subjected to high-density genotyping using genotyping by sequencing. We performed genome-wide association mapping to characterize the genetic architecture and gene mapping of HCN. Field experiments revealed strong broad- and narrow-sense trait heritability (0.82 and 0.41, respectively). Two major loci were identified, encoding for an ATPase and a MATE protein, and contributing up to 7 and 30% of the HCN concentration in roots, respectively. We developed diagnostic markers for breeding applications, validated trait architecture consistency in African germplasm and investigated further evidence for the domestication of sweet and bitter cassava. Fine genomic characterization revealed: (i) the major role played by vacuolar transporters in regulating HCN content; (ii) the co-domestication of sweet and bitter cassava major alleles are dependent upon geographical zone; and (iii) the major loci allele for high HCN in M. esculenta Crantz seems to originate from its ancestor, M. esculenta subsp. flabellifolia. Taken together, these findings expand our insights into cyanogenic glucosides in cassava roots and its glycosylated derivatives in plants.

Selection of a subspecies-specific diterpene gene cluster implicated in rice disease resistance

Diterpenoids are the major group of antimicrobial phytoalexins in rice^1,2. Here, we report the discovery of a rice diterpenoid gene cluster on chromosome 7 (DGC7) encoding the entire biosynthetic pathway to 5,10-diketo-casbene, a member of the monocyclic casbene-derived diterpenoids. We revealed that DGC7 is regulated directly by JMJ705 through methyl jasmonate-mediated epigenetic control³. Functional characterization of pathway genes revealed OsCYP71Z21 to encode a casbene C10 oxidase, sought after for the biosynthesis of an array of medicinally important diterpenoids. We further show that DGC7 arose relatively recently in the Oryza genus, and that it was partly formed in Oryza rufipogon and positively selected for in japonica during domestication. Casbene-synthesizing enzymes that are functionally equivalent to OsTPS28 are present in several species of Euphorbiaceae but gene tree analysis shows that these and other casbene-modifying enzymes have evolved independently. As such, combining casbene-modifying enzymes from these different families of plants may prove effective in producing a diverse array of bioactive diterpenoid natural products.

A flavin-dependent monooxygenase catalyzes the initial step in cyanogenic glycoside synthesis in ferns

Cyanogenic glycosides form part of a binary plant defense system that, upon catabolism, detonates a toxic hydrogen cyanide bomb. In seed plants, the initial step of cyanogenic glycoside biosynthesis-the conversion of an amino acid to the corresponding aldoxime-is catalyzed by a cytochrome P450 from the CYP79 family. An evolutionary conundrum arises, as no CYP79s have been identified in ferns, despite cyanogenic glycoside occurrence in several fern species. Here, we report that a flavin-dependent monooxygenase (fern oxime synthase; FOS1), catalyzes the first step of cyanogenic glycoside biosynthesis in two fern species (Phlebodium aureum and Pteridium aquilinum), demonstrating convergent evolution of biosynthesis across the plant kingdom. The FOS1 sequence from the two species is near identical (98%), despite diversifying 140 MYA. Recombinant FOS1 was isolated as a catalytic active dimer, and in planta, catalyzes formation of an N-hydroxylated primary amino acid; a class of metabolite not previously observed in plants.

Biosynthesis of cyanogenic glucosides in Phaseolus lunatus and the evolution of oxime-based defenses

Lima bean, Phaseolus lunatus, is a crop legume that produces the cyanogenic glucosides linamarin and lotaustralin. In the legumes Lotus japonicus and Trifolium repens, the biosynthesis of these two α-hydroxynitrile glucosides involves cytochrome P450 enzymes of the CYP79 and CYP736 families and a UDP-glucosyltransferase. Here, we identify CYP79D71 as the first enzyme of the pathway in P. lunatus, producing oximes from valine and isoleucine. A second CYP79 family member, CYP79D72, was shown to catalyze the formation of leucine-derived oximes, which act as volatile defense compounds in Phaseolus spp. The organization of the biosynthetic genes for cyanogenic glucosides in a gene cluster aided their identification in L. japonicus. In the available genome sequence of P. vulgaris, the gene orthologous to CYP79D71 is adjacent to a member of the CYP83 family. Although P. vulgaris is not cyanogenic, it does produce oximes as volatile defense compounds. We cloned the genes encoding two CYP83s (CYP83E46 and CYP83E47) and a UDP-glucosyltransferase (UGT85K31) from P. lunatus, and these genes combined form a complete biosynthetic pathway for linamarin and lotaustralin in Lima bean. Within the genus Phaseolus, the occurrence of linamarin and lotaustralin as functional chemical defense compounds appears restricted to species belonging to the closely related Polystachios and Lunatus groups. A preexisting ability to produce volatile oximes and nitriles likely facilitated evolution of cyanogenesis within the Phaseolus genus.

Drivers of metabolic diversification: how dynamic genomic neighbourhoods generate new biosynthetic pathways in the Brassicaceae

Plants produce an array of specialized metabolites with important ecological functions. The mechanisms underpinning the evolution of new biosynthetic pathways are not well-understood. Here, we exploit available genome sequence resources to investigate triterpene biosynthesis across the Brassicaceae. Oxidosqualene cyclases (OSCs) catalyze the first committed step in triterpene biosynthesis. Systematic analysis of 13 sequenced Brassicaceae genomes was performed to identify all OSC genes. The genome neighbourhoods (GNs) around a total of 163 OSC genes were investigated to identify Pfam domains significantly enriched in these regions. All-vs-all comparisons of OSC neighbourhoods and phylogenomic analysis were used to investigate the sequence similarity and evolutionary relationships of the numerous candidate triterpene biosynthetic gene clusters (BGCs) observed. Functional analysis of three representative BGCs was carried out and their triterpene pathway products were elucidated. Our results indicate that plant genomes are remarkably plastic, and that dynamic GNs generate new biosynthetic pathways in different Brassicaceae lineages by shuffling the genes encoding a core palette of triterpene-diversifying enzymes, presumably in response to strong environmental selection pressure. These results illuminate a genomic basis for diversification of plant-specialized metabolism through natural combinatorics of enzyme families, which can be mimicked using synthetic biology to engineer diverse bioactive molecules.

Evolution of a plant gene cluster in Solanaceae and emergence of metabolic diversity

Plants produce phylogenetically and spatially restricted, as well as structurally diverse specialized metabolites via multistep metabolic pathways. Hallmarks of specialized metabolic evolution include enzymatic promiscuity and recruitment of primary metabolic enzymes and examples of genomic clustering of pathway genes. Solanaceae glandular trichomes produce defensive acylsugars, with sidechains that vary in length across the family. We describe a tomato gene cluster on chromosome 7 involved in medium chain acylsugar accumulation due to trichome specific acyl-CoA synthetase and enoyl-CoA hydratase genes. This cluster co-localizes with a tomato steroidal alkaloid gene cluster and is syntenic to a chromosome 12 region containing another acylsugar pathway gene. We reconstructed the evolutionary events leading to this gene cluster and found that its phylogenetic distribution correlates with medium chain acylsugar accumulation across the Solanaceae. This work reveals insights into the dynamics behind gene cluster evolution and cell-type specific metabolite diversity.

Ab initio GO-based mining for non-tandem-duplicated functional clusters in three model plant diploid genomes

A functional Non-Tandem Duplicated Cluster (FNTDC) is a group of non-tandem-duplicated genes that are located closer than expected by mere chance and have a role in the same biological function. The identification of secondary-compounds–related FNTDC has gained increased interest in recent years, but little ab-initio attempts aiming to the identification of FNTDCs covering all biological functions, including primary metabolism compounds, have been carried out. We report an extensive FNTDC dataset accompanied by a detailed assessment on parameters used for genome scanning and their impact on FNTDC detection. We propose 70% identity and 70% alignment coverage as intermediate settings to exclude tandem duplicated genes and a dynamic scanning window of 24 genes. These settings were applied to rice, arabidopsis and grapevine genomes to call for FNTDCs. Besides the best-known secondary metabolism clusters, we identified many FNTDCs associated to primary metabolism ranging from macromolecules synthesis/editing, TOR signalling, ubiquitination, proton and electron transfer complexes. Using the intermediate FNTDC setting parameters (at P-value 1e^-6), 130, 70 and 140 candidate FNTDCs were called in rice, arabidopsis and grapevine, respectively, and 20 to 30% of GO tags associated to called FNTDC were common among the 3 genomes. The datasets developed along with this work provide a rich framework for pinpointing candidate FNTDCs reflecting all GO-BP tags covering both primary and secondary metabolism with large macromolecular complexes/metabolons as the most represented FNTDCs. Noteworthy, several FNTDCs are tagged with GOs referring to organelle-targeted multi-enzyme complex, a finding that suggest the migration of endosymbiont gene chunks towards nuclei could be at the basis of these class of candidate FNTDCs. Most FNTDC appear to have evolved prior of genome duplication events. More than one-third of genes interspersed/adjacent to called FNTDCs lacked any functional annotation; however, their co-localization may provide hints towards a candidate biological role.

Genomics-enabled analysis of specialized metabolism in bioenergy crops: current progress and challenges

Plants produce a staggering diversity of specialized small molecule metabolites that play vital roles in mediating environmental interactions and stress adaptation. This chemical diversity derives from dynamic biosynthetic pathway networks that are often species-specific and operate under tight spatiotemporal and environmental control. A growing divide between demand and environmental challenges in food and bioenergy crop production has intensified research on these complex metabolite networks and their contribution to crop fitness. High-throughput omics technologies provide access to ever-increasing data resources for investigating plant metabolism. However, the efficiency of using such system-wide data to decode the gene and enzyme functions controlling specialized metabolism has remained limited; due largely to the recalcitrance of many plants to genetic approaches and the lack of 'user-friendly' biochemical tools for studying the diverse enzyme classes involved in specialized metabolism. With emphasis on terpenoid metabolism in the bioenergy crop switchgrass as an example, this review aims to illustrate current advances and challenges in the application of DNA synthesis and synthetic biology tools for accelerating the functional discovery of genes, enzymes and pathways in plant specialized metabolism. These technologies have accelerated knowledge development on the biosynthesis and physiological roles of diverse metabolite networks across many ecologically and economically important plant species and can provide resources for application to precision breeding and natural product metabolic engineering.

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

Current findings of neighboring genes involved in plant specialized metabolism provide the genomic signatures of metabolic evolution. Two such genomic features, namely, (i) metabolic gene cluster and (ii) neo-functionalization of tandem gene duplications, represent key factors corresponding to the creation of metabolic diversity of plant specialized metabolism. So far, several terpenoid and alkaloid biosynthetic genes have been characterized with gene clusters in some plants. On the other hand, some modification genes involved in flavonoid and glucosinolate biosynthesis were found to arise via gene neo-functionalization. Although the occurrence of both types of metabolic evolution are different, the neighboring genes are generally regulated by the same or related regulation factors. Therefore, the translation-based approaches associated with genomics, and transcriptomics are able to be employed for functional genomics focusing on plant secondary metabolism. Here, we present a survey of the current understanding of neighboring genes involved in plant secondary metabolism. Additionally, a genomic overview of neighboring genes of four model plants and transcriptional co-expression network neighboring genes to detect metabolic gene clusters in Arabidopsis is provided. Finally, the insights functional genomics have provided concerning the evolution and mechanistic regulation of both the formation and operation of metabolic neighboring clusters is discussed.

The Sweet Side of Plant-Specialized Metabolism

Glycosylation plays a major role in the structural diversification of plant natural products. It influences the properties of molecules by modifying the reactivity and solubility of the corresponding aglycones, so influencing cellular localization and bioactivity. Glycosylation of plant natural products is usually carried out by uridine diphosphate(UDP)-dependent glycosyltransferases (UGTs) belonging to the carbohydrate-active enzyme glycosyltransferase 1 (GT1) family. These enzymes transfer sugars from UDP-activated sugar moieties to small hydrophobic acceptor molecules. Plant GT1s generally show high specificity for their sugar donors and recognize a single UDP sugar as their substrate. In contrast, they are generally promiscuous with regard to acceptors, making them attractive biotechnological tools for small molecule glycodiversification. Although microbial hosts have traditionally been used for heterologous engineering of plant-derived glycosides, transient plant expression technology offers a potentially disruptive platform for rapid characterization of new plant glycosyltransferases and biosynthesis of complex glycosides.

Extractive silylation method for high throughput GC analysis of flaxseed cyanogenic glycosides

The reported methods for the estimation of CG are indirect, long and tedious uses solvent extraction which results in lower recoveries due to several co-eluted components. The present study aims to develop and validate a high throughput method for the analysis of cyanogenic glycosides (CG) in flaxseed using extractive silylation. The experimental procedure comprised of preparation of trimethylsilyl (TMS) derivatives of CG, their quantitation through gas chromatography with flame ionization detector (GC-FID) and further characterization by LC-MS/MS. Different validation parameters determined in the experiment include the relative standard deviation both inter-day and intraday less than 5%, recovery in the range of 79.9-112.7%, limit of detection 4.72-6.43 µg/mL and limit of quantitation 14.31-19.50 µg/mL. Combinations of silylation reagent were screened in a central composite experimental design in order to examine their effect on the extraction recovery. Finally, the developed method was applied successfully to quantify CG in various flaxseed cultivars. Advantages of the extractive silylation are simple preparation, short reaction times and the combination of extraction and silylation in one step which indicate that the method has the potential to sensitively and accurately determine CG where large numbers of samples are being routinely analysed.