Gene clustering in plant specialized metabolism

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.

Active and repressed biosynthetic gene clusters have spatially distinct chromosome states

While colocalization within a bacterial operon enables coexpression of the constituent genes, the mechanistic logic of clustering of nonhomologous monocistronic genes in eukaryotes is not immediately obvious. Biosynthetic gene clusters that encode pathways for specialized metabolites are an exception to the classical eukaryote rule of random gene location and provide paradigmatic exemplars with which to understand eukaryotic cluster dynamics and regulation. Here, using 3C, Hi-C, and Capture Hi-C (CHi-C) organ-specific chromosome conformation capture techniques along with high-resolution microscopy, we investigate how chromosome topology relates to transcriptional activity of clustered biosynthetic pathway genes in Arabidopsis thaliana Our analyses reveal that biosynthetic gene clusters are embedded in local hot spots of 3D contacts that segregate cluster regions from the surrounding chromosome environment. The spatial conformation of these cluster-associated domains differs between transcriptionally active and silenced clusters. We further show that silenced clusters associate with heterochromatic chromosomal domains toward the periphery of the nucleus, while transcriptionally active clusters relocate away from the nuclear periphery. Examination of chromosome structure at unrelated clusters in maize, rice, and tomato indicates that integration of clustered pathway genes into distinct topological domains is a common feature in plant genomes. Our results shed light on the potential mechanisms that constrain coexpression within clusters of nonhomologous eukaryotic genes and suggest that gene clustering in the one-dimensional chromosome is accompanied by compartmentalization of the 3D chromosome.

Genomic evidence for convergent evolution of gene clusters for momilactone biosynthesis in land plants

Momilactones are bioactive diterpenoids that contribute to plant defense against pathogens and allelopathic interactions between plants. Both cultivated and wild grass species of Oryza and Echinochloa crus-galli (barnyard grass) produce momilactones using a biosynthetic gene cluster (BGC) in their genomes. The bryophyte Calohypnum plumiforme (formerly Hypnum plumaeforme) also produces momilactones, and the bifunctional diterpene cyclase gene CpDTC1/HpDTC1, which is responsible for the production of the diterpene framework, has been characterized. To understand the molecular architecture of the momilactone biosynthetic genes in the moss genome and their evolutionary relationships with other momilactone-producing plants, we sequenced and annotated the C. plumiforme genome. The data revealed a 150-kb genomic region that contains two cytochrome P450 genes, the CpDTC1/HpDTC1 gene and the "dehydrogenase momilactone A synthase" gene tandemly arranged and inductively transcribed following stress exposure. The predicted enzymatic functions in yeast and recombinant assay and the successful pathway reconstitution in Nicotiana benthamiana suggest that it is a functional BGC responsible for momilactone production. Furthermore, in a survey of genomic sequences of a broad range of plant species, we found that momilactone BGC is limited to the two grasses (Oryza and Echinochloa) and C. plumiforme, with no synteny among these genomes. These results indicate that while the gene cluster in C. plumiforme is functionally similar to that in rice and barnyard grass, it is likely a product of convergent evolution. To the best of our knowledge, this report of a BGC for a specialized plant defense metabolite in bryophytes is unique.

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.

Gene clustering and copy number variation in alkaloid metabolic pathways of opium poppy

Genes in plant secondary metabolic pathways enable biosynthesis of a range of medically and industrially important compounds, and are often clustered on chromosomes. Here, we study genomic clustering in the benzylisoquinoline alkaloid (BIA) pathway in opium poppy (Papaver somniferum), exploring relationships between gene expression, copy number variation, and metabolite production. We use Hi-C to improve the existing draft genome assembly, yielding chromosome-scale scaffolds that include 35 previously unanchored BIA genes. We find that co-expression of BIA genes increases within clusters and identify candidates with unknown function based on clustering and covariation in expression and alkaloid production. Copy number variation in critical BIA genes correlates with stark differences in alkaloid production, linking noscapine production with an 11-gene deletion, and increased thebaine/decreased morphine production with deletion of a T6ODM cluster. Our results show that the opium poppy genome is still dynamically evolving in ways that contribute to medically and industrially important phenotypes.

Genome of Tripterygium wilfordii and identification of cytochrome P450 involved in triptolide biosynthesis

Triptolide is a trace natural product of Tripterygium wilfordii. It has antitumor activities, particularly against pancreatic cancer cells. Identification of genes and elucidation of the biosynthetic pathway leading to triptolide are the prerequisite for heterologous bioproduction. Here, we report a reference-grade genome of T. wilfordii with a contig N50 of 4.36 Mb. We show that copy numbers of triptolide biosynthetic pathway genes are impacted by a recent whole-genome triplication event. We further integrate genomic, transcriptomic, and metabolomic data to map a gene-to-metabolite network. This leads to the identification of a cytochrome P450 (CYP728B70) that can catalyze oxidation of a methyl to the acid moiety of dehydroabietic acid in triptolide biosynthesis. We think the genomic resource and the candidate genes reported here set the foundation to fully reveal triptolide biosynthetic pathway and consequently the heterologous bioproduction.

Tripterygium wilfordii is a medical plant that can produce antitumor activity compound triptolide. Here, the authors assemble its genome and identify a cytochrome P450 that can catalyze oxidation of a methyl to the acid moiety of dehydroabietic acid in triptolide biosynthesis by integrating multi-omics data.

A chromosome-anchored eggplant genome sequence reveals key events in Solanaceae evolution

With approximately 450 species, spiny Solanum species constitute the largest monophyletic group in the Solanaceae family, but a high-quality genome assembly from this group is presently missing. We obtained a chromosome-anchored genome assembly of eggplant (Solanum melongena), containing 34,916 genes, confirming that the diploid gene number in the Solanaceae is around 35,000. Comparative genomic studies with tomato (S. lycopersicum), potato (S. tuberosum) and pepper (Capsicum annuum) highlighted the rapid evolution of miRNA:mRNA regulatory pairs and R-type defense genes in the Solanaceae, and provided a genomic basis for the lack of steroidal glycoalkaloid compounds in the Capsicum genus. Using parsimony methods, we reconstructed the putative chromosomal complements of the key founders of the main Solanaceae clades and the rearrangements that led to the karyotypes of extant species and their ancestors. From 10% to 15% of the genes present in the four genomes were syntenic paralogs (ohnologs) generated by the pre-γ, γ and T paleopolyploidy events, and were enriched in transcription factors. Our data suggest that the basic gene network controlling fruit ripening is conserved in different Solanaceae clades, and that climacteric fruit ripening involves a differential regulation of relatively few components of this network, including CNR and ethylene biosynthetic genes.

plantiSMASH: automated identification, annotation and expression analysis of plant biosynthetic gene clusters

Plant specialized metabolites are chemically highly diverse, play key roles in host-microbe interactions, have important nutritional value in crops and are frequently applied as medicines. It has recently become clear that plant biosynthetic pathway-encoding genes are sometimes densely clustered in specific genomic loci: biosynthetic gene clusters (BGCs). Here, we introduce plantiSMASH, a versatile online analysis platform that automates the identification of candidate plant BGCs. Moreover, it allows integration of transcriptomic data to prioritize candidate BGCs based on the coexpression patterns of predicted biosynthetic enzyme-coding genes, and facilitates comparative genomic analysis to study the evolutionary conservation of each cluster. Applied on 48 high-quality plant genomes, plantiSMASH identifies a rich diversity of candidate plant BGCs. These results will guide further experimental exploration of the nature and dynamics of gene clustering in plant metabolism. Moreover, spurred by the continuing decrease in costs of plant genome sequencing, they will allow genome mining technologies to be applied to plant natural product discovery. The plantiSMASH web server, precalculated results and source code are freely available from http://plantismash.secondarymetabolites.org.

Natural variations in the MYB transcription factor MYB31 determine the evolution of extremely pungent peppers

Plants produce countless specialized metabolites crucial for their development and fitness, and many are useful bioactive compounds. Capsaicinoids are intriguing genus-specialized metabolites that confer a pungent flavor to Capsicum fruits, and they are widely applied in different areas. Among the five domesticated Capsicum species, Capsicum chinense has a high content of capsaicinoids, which results in an extremely hot flavor. However, the species-specific upregulation of capsaicinoid-biosynthetic genes (CBGs) and the evolution of extremely pungent peppers are not well understood. We conducted genetic and functional analyses demonstrating that the quantitative trait locus Capsaicinoid1 (Cap1), which is identical to Pun3 contributes to the level of pungency. The Cap1/Pun3 locus encodes the Solanaceae-specific MYB transcription factor MYB31. Capsicum species have evolved placenta-specific expression of MYB31, which directly activates expression of CBGs and results in genus-specialized metabolite production. The capsaicinoid content depends on MYB31 expression. Natural variations in the MYB31 promoter increase MYB31 expression in C. chinense via the binding of the placenta-specific expression of transcriptional activator WRKY9 and augmentation of CBG expression, which promotes capsaicinoid biosynthesis. Our findings provide insights into the evolution of extremely pungent C. chinense, which is due to natural variations in the master regulator, and offers targets for engineering or selecting flavor in Capsicum.

The Reference Genome Sequence of Scutellaria baicalensis Provides Insights into the Evolution of Wogonin Biosynthesis

Scutellaria baicalensis Georgi is important in Chinese traditional medicine where preparations of dried roots, "Huang Qin," are used for liver and lung complaints and as complementary cancer treatments. We report a high-quality reference genome sequence for S. baicalensis where 93% of the 408.14-Mb genome has been assembled into nine pseudochromosomes with a super-N50 of 33.2 Mb. Comparison of this sequence with those of closely related species in the order Lamiales, Sesamum indicum and Salvia splendens, revealed that a specialized metabolic pathway for the synthesis of 4'-deoxyflavone bioactives evolved in the genus Scutellaria. We found that the gene encoding a specific cinnamate coenzyme A ligase likely obtained its new function following recent mutations, and that four genes encoding enzymes in the 4'-deoxyflavone pathway are present as tandem repeats in the genome of S. baicalensis. Further analyses revealed that gene duplications, segmental duplication, gene amplification, and point mutations coupled to gene neo- and subfunctionalizations were involved in the evolution of 4'-deoxyflavone synthesis in the genus Scutellaria. Our study not only provides significant insight into the evolution of specific flavone biosynthetic pathways in the mint family, Lamiaceae, but also will facilitate the development of tools for enhancing bioactive productivity by metabolic engineering in microbes or by molecular breeding in plants. The reference genome of S. baicalensis is also useful for improving the genome assemblies for other members of the mint family and offers an important foundation for decoding the synthetic pathways of bioactive compounds in medicinal plants.

Cytochrome P450 enzymes: A driving force of plant diterpene diversity

In plant terpene biosynthesis, oxidation of the hydrocarbon backbone produced by terpene synthases is typically carried out by cytochrome P450 oxygenases (CYPs). The modifications introduced by CYPs include hydroxylations, sequential oxidations at one position and ring rearrangements and closures. These reactions significantly expand the structural diversity of terpenoids, but also provide anchoring points for further decorations by various transferases. In recent years, there has been a significant increase in reports of CYPs involved in plant terpene pathways. Plant diterpenes represent an important class of metabolites that includes hormones and a number of industrially relevant compounds such as pharmaceutical, aroma or food ingredients. In this review, we provide a comprehensive survey on CYPs reported to be involved in plant diterpene biosynthesis to date. A phylogenetic analysis showed that only few CYP clans are represented in diterpene biosynthesis, namely CYP71, CYP85 and CYP72. Remarkably few CYP families and subfamilies within those clans are involved, indicating specific expansion of these clades in plant diterpene biosynthesis. Nonetheless, the evolutionary trajectory of CYPs of specialized diterpene biosynthesis is diverse. Some are recently derived from gibberellin biosynthesis, while others have a more ancient history with recent expansions in specific plant families. Among diterpenoids, labdane-related diterpenoids represent a dominant class. The availability of CYPs from diverse plant species able to catalyze oxidations in specific regions of the labdane-related backbones provides opportunities for combinatorial biosynthesis to produce novel diterpene compounds that can be screened for biological activities of interest.

Charting Angiosperm Chemistry: Evolutionary Perspective on Specialized Metabolites Reflected in Chemical Property Space

Plants possess an outstanding chemical diversity of specialized metabolites developed to adapt to environmental niches and increase fitness. The observed diversity is hypothesized to result from various evolutionary mechanisms, such as the continuous branching off and extension of existing biosynthetic pathways or enhanced levels of catalytic promiscuity in certain enzymes. In this study, ChemGPS-NP has been employed to chart the distribution and diversity of physicochemical properties for selected types of specialized metabolites from the angiosperms. Utilizing these charts, it is analyzed how different properties of various types of specialized metabolites change in different plant groups, and the chemical diversity from the volume they occupy in chemical property space is evaluated. In this context, possible underlying evolutionary mechanisms are discussed, which could explain the observed distribution and behavior in chemical property space. Based on these studies, it is demonstrated that evolutionary processes in plant specialized metabolism and the resultant metabolic diversification are reflected in chemical property space.

Epigenetic modification, co-culture and genomic methods for natural product discovery

Fungi and bacteria are encountered in many habitats where they live in complex communities interacting with one another mainly by producing secondary metabolites, which are organic compounds that are not directly involved in the normal growth, development, or reproduction of the organism. These organisms appear as a promising source for the discovery of novel bioactive natural products that may find their application in medicine. However, the production of secondary metabolites by those organisms when cultured axenically is limited as only a subset of biosynthetic genes is expressed under standard laboratory conditions leading to the search of new methods for the activation of the silent genes including epigenetic modification and co-cultivation. Biosynthetic gene clusters which produce secondary metabolites are known to be present in a heterochromatin state in which the transcription of constitutive genes is usually regulated by epigenetic modification including DNA methylation and histone deacetylation. Therefore, small-molecule epigenetic modifiers which promote changes in the structure of chromatin could control the expression of silent genes and may be rationally employed for the discovery of novel bioactive compounds. Co-cultivation, which is also known as mixed-fermentation, usually implies two or more microorganisms in the same medium in which the resulting competition is known to enhance the production of constitutively present compounds and/or to lead to the induction of cryptic metabolites that were not detected in axenic cultures of the considered axenic microorganism. Genomic strategies could help to identify biosynthetic gene clusters in fungal genomes and link them to their products by the means of novel algorithms as well as integrative pan-genomic approaches. Despite that all these techniques are still in their infancy, they appear as promising sources for the discovery of new bioactive compounds. This chapter presents recent ecological techniques for the discovery of new secondary metabolites that might find application in medicine.

Genome sequence of Malania oleifera, a tree with great value for nervonic acid production

BACKGROUND

Malania oleifera, a member of the Olacaceae family, is an IUCN red listed tree, endemic and restricted to the Karst region of southwest China. This tree's seed is valued for its high content of precious fatty acids (especially nervonic acid). However, studies on its genetic makeup and fatty acid biogenesis are severely hampered by a lack of molecular and genetic tools.

FINDINGS

We generated 51 Gb and 135 Gb of raw DNA sequences, using Pacific Biosciences (PacBio) single-molecule real-time and 10× Genomics sequencing, respectively. A final genome assembly, with a scaffold N50 size of 4.65 Mb and a total length of 1.51 Gb, was obtained by primary assembly based on PacBio long reads plus scaffolding with 10× Genomics reads. Identified repeats constituted ∼82% of the genome, and 24,064 protein-coding genes were predicted with high support. The genome has low heterozygosity and shows no evidence for recent whole genome duplication. Metabolic pathway genes relating to the accumulation of long-chain fatty acid were identified and studied in detail.

CONCLUSIONS

Here, we provide the first genome assembly and gene annotation for M. oleifera. The availability of these resources will be of great importance for conservation biology and for the functional genomics of nervonic acid biosynthesis.

Finding and analysing the minimum set of driver nodes required to control multilayer networks

It is difficult to control multilayer networks in situations with real-world complexity. Here, we first define the multilayer control problem in terms of the minimum dominating set (MDS) controllability framework and mathematically demonstrate that simple formulas can be used to estimate the size of the minimum dominating set in multilayer (MDSM) complex networks. Second, we develop a new algorithm that efficiently identifies the MDSM in up to 6 layers, with several thousand nodes in each layer network. Interestingly, the findings reveal that the MDSM size for similar networks does not significantly differ from that required to control a single network. This result opens future directions for controlling, for example, multiple species by identifying a common set of enzymes or proteins for drug targeting. We apply our methods to 70 genome-wide metabolic networks across major plant lineages, unveiling some relationships between controllability in multilayer networks and metabolic functions at the genome scale.

A physical and genetic map of Cannabis sativa identifies extensive rearrangements at the THC/CBD acid synthase loci

Cannabis sativa is widely cultivated for medicinal, food, industrial, and recreational use, but much remains unknown regarding its genetics, including the molecular determinants of cannabinoid content. Here, we describe a combined physical and genetic map derived from a cross between the drug-type strain Purple Kush and the hemp variety "Finola." The map reveals that cannabinoid biosynthesis genes are generally unlinked but that aromatic prenyltransferase (AP), which produces the substrate for THCA and CBDA synthases (THCAS and CBDAS), is tightly linked to a known marker for total cannabinoid content. We further identify the gene encoding CBCA synthase (CBCAS) and characterize its catalytic activity, providing insight into how cannabinoid diversity arises in cannabis. THCAS and CBDAS (which determine the drug vs. hemp chemotype) are contained within large (>250 kb) retrotransposon-rich regions that are highly nonhomologous between drug- and hemp-type alleles and are furthermore embedded within ∼40 Mb of minimally recombining repetitive DNA. The chromosome structures are similar to those in grains such as wheat, with recombination focused in gene-rich, repeat-depleted regions near chromosome ends. The physical and genetic map should facilitate further dissection of genetic and molecular mechanisms in this commercially and medically important plant.

The Chrysanthemum nankingense Genome Provides Insights into the Evolution and Diversification of Chrysanthemum Flowers and Medicinal Traits

The Asteraceae (Compositae), a large plant family of approximately 24 000-35 000 species, accounts for ∼10% of all angiosperm species and contributes a lot to plant diversity. The most representative members of the Asteraceae are the economically important chrysanthemums (Chrysanthemum L.) that diversified through reticulate evolution. Biodiversity is typically created by multiple evolutionary mechanisms such as whole-genome duplication (WGD) or polyploidization and locally repetitive genome expansion. However, the lack of genomic data from chrysanthemum species has prevented an in-depth analysis of the evolutionary mechanisms involved in their diversification. Here, we used Oxford Nanopore long-read technology to sequence the diploid Chrysanthemum nankingense genome, which represents one of the progenitor genomes of domesticated chrysanthemums. Our analysis revealed that the evolution of the C. nankingense genome was driven by bursts of repetitive element expansion and WGD events including a recent WGD that distinguishes chrysanthemum from sunflower, which diverged from chrysanthemum approximately 38.8 million years ago. Variations of ornamental and medicinal traits in chrysanthemums are linked to the expansion of candidate gene families by duplication events including paralogous gene duplication. Collectively, our study of the assembled reference genome offers new knowledge and resources to dissect the history and pattern of evolution and diversification of chrysanthemum plants, and also to accelerate their breeding and improvement.

Enhancing the abiotic stress tolerance of plants: from chemical treatment to biotechnological approaches

Abiotic stresses affect crop plants and cause decreases in plant quality and productivity. Plants can overcome environmental stresses by activating molecular networks, including signal transduction, stress perception, metabolite production and expressions of specific stress-related genes. Recent research suggests that chemical priming is a promising field in crop stress management because plants can be primed by chemical agents to increase their tolerance to various environmental stresses. We present a concept to meet this objective and protect plants through priming of existing defense mechanisms avoiding manipulation of the genome. In addition, recent developments in plant molecular biology include the discovery of genes related to stress tolerance, including functional genes for protecting cells and regulatory genes for regulating stress responses. Therefore, enhancing abiotic stress tolerance using a transgenic approach to transfer these genes into plant genomes has attracted more investigations. Both chemical priming agents and genetic engineering can enhance regulatory and functional genes in plants and increase stress tolerance of plants. This review summarizes the latest findings of chemical priming agents and major achievements in molecular approaches that can potentially enhance the abiotic stress tolerance of plants.

A physical and genetic map of Cannabis sativa identifies extensive rearrangements at the THC/CBD acid synthase loci

Cannabis sativa is widely cultivated for medicinal, food, industrial, and recreational use, but much remains unknown regarding its genetics, including the molecular determinants of cannabinoid content. Here, we describe a combined physical and genetic map derived from a cross between the drug-type strain Purple Kush and the hemp variety “Finola.” The map reveals that cannabinoid biosynthesis genes are generally unlinked but that aromatic prenyltransferase (AP), which produces the substrate for THCA and CBDA synthases (THCAS and CBDAS), is tightly linked to a known marker for total cannabinoid content. We further identify the gene encoding CBCA synthase (CBCAS) and characterize its catalytic activity, providing insight into how cannabinoid diversity arises in cannabis. THCAS and CBDAS (which determine the drug vs. hemp chemotype) are contained within large (>250 kb) retrotransposon-rich regions that are highly nonhomologous between drug- and hemp-type alleles and are furthermore embedded within ∼40 Mb of minimally recombining repetitive DNA. The chromosome structures are similar to those in grains such as wheat, with recombination focused in gene-rich, repeat-depleted regions near chromosome ends. The physical and genetic map should facilitate further dissection of genetic and molecular mechanisms in this commercially and medically important plant.

Challenges and emergent solutions for LC-MS/MS based untargeted metabolomics in diseases

In the past decade, advances in liquid chromatography-mass spectrometry (LC-MS) have revolutionized untargeted metabolomics analyses. By mining metabolomes more deeply, researchers are now primed to uncover key metabolites and their associations with diseases. The employment of untargeted metabolomics has led to new biomarker discoveries and a better mechanistic understanding of diseases with applications in precision medicine. However, many major pertinent challenges remain. First, compound identification has been poor, and left an overwhelming number of unidentified peaks. Second, partial, incomplete metabolomes persist due to factors such as limitations in mass spectrometry data acquisition speeds, wide-range of metabolites concentrations, and cellular/tissue/temporal-specific expression changes that confound our understanding of metabolite perturbations. Third, to contextualize metabolites in pathways and biology is difficult because many metabolites partake in multiple pathways, have yet to be described species specificity, or possess unannotated or more-complex functions that are not easily characterized through metabolomics analyses. From a translational perspective, information related to novel metabolite biomarkers, metabolic pathways, and drug targets might be sparser than they should be. Thankfully, significant progress has been made and novel solutions are emerging, achieved through sustained academic and industrial community efforts in terms of hardware, computational, and experimental approaches. Given the rapidly growing utility of metabolomics, this review will offer new perspectives, increase awareness of the major challenges in LC-MS metabolomics that will significantly benefit the metabolomics community and also the broader the biomedical community metabolomics aspire to serve.

Bioinformatics tools for the identification of gene clusters that biosynthesize specialized metabolites

Specialized metabolites (also called natural products or secondary metabolites) derived from bacteria, fungi, marine organisms and plants constitute an important source of antibiotics, anti-cancer agents, insecticides, immunosuppressants and herbicides. Many specialized metabolites in bacteria and fungi are biosynthesized via metabolic pathways whose enzymes are encoded by clustered genes on a chromosome. Metabolic gene clusters comprise a group of physically co-localized genes that together encode enzymes for the biosynthesis of a specific metabolite. Although metabolic gene clusters are generally not known to occur outside of microbes, several plant metabolic gene clusters have been discovered in recent years. The discovery of novel metabolic pathways is being enabled by the increasing availability of high-quality genome sequencing coupled with the development of powerful computational toolkits to identify metabolic gene clusters. To provide a comprehensive overview of various bioinformatics methods for detecting gene clusters, we compare and contrast key aspects of algorithmic logic behind several computational tools, including 'NP.searcher', 'ClustScan', 'CLUSEAN', 'antiSMASH', 'SMURF', 'MIDDAS-M', 'ClusterFinder', 'CASSIS/SMIPS' and 'C-Hunter' among others. We also review additional tools such as 'NRPSpredictor' and 'SBSPKS' that can infer substrate specificity for previously identified gene clusters. The continual development of bioinformatics methods to predict gene clusters will help shed light on how organisms assemble multi-step metabolic pathways for adaptation to various ecological niches.

Specialized Plant Metabolism Characteristics and Impact on Target Molecule Biotechnological Production

Plant secondary metabolism evolved in the context of highly organized and differentiated cells and tissues, featuring massive chemical complexity operating under tight environmental, developmental and genetic control. Biotechnological demand for natural products has been continuously increasing because of their significant value and new applications, mainly as pharmaceuticals. Aseptic production systems of plant secondary metabolites have improved considerably, constituting an attractive tool for increased, stable and large-scale supply of valuable molecules. Surprisingly, to date, only a few examples including taxol, shikonin, berberine and artemisinin have emerged as success cases of commercial production using this strategy. The present review focuses on the main characteristics of plant specialized metabolism and their implications for current strategies used to produce secondary compounds in axenic cultivation systems. The search for consonance between plant secondary metabolism unique features and various in vitro culture systems, including cell, tissue, organ, and engineered cultures, as well as heterologous expression in microbial platforms, is discussed. Data to date strongly suggest that attaining full potential of these biotechnology production strategies requires being able to take advantage of plant specialized metabolism singularities for improved target molecule yields and for bypassing inherent difficulties in its rational manipulation.

Beyond pathways: genetic dissection of tocopherol content in maize kernels by combining linkage and association analyses

Although tocopherols play an important role in plants and animals, the genetic architecture of tocopherol content in maize kernels has remained largely unknown. In this study, linkage and association analyses were conducted to examine the genetic architecture of tocopherol content in maize kernels. Forty-one unique quantitative trait loci (QTLs) were identified by linkage mapping in six populations of recombinant inbred lines (RILs). In addition, 32 significant loci were detected via genome-wide association study (GWAS), 18 of which colocalized with the QTLs identified by linkage mapping. Fine mapping of a major QTL validated the accuracy of GWAS and QTL mapping results and suggested a role for nontocopherol pathway genes in the modulation of natural tocopherol variation. We provided genome-wide evidence that genes involved in fatty acid metabolism, chlorophyll metabolism and chloroplast function may affect natural variation in tocopherols. These findings were confirmed through mutant analysis of a particular gene from the fatty acid pathway. In addition, the favourable alleles for many of the significant SNPs/QTLs represented rare alleles in natural populations. Together, our results revealed many novel genes that are potentially involved in the variation of tocopherol content in maize kernels. Pyramiding of the favourable alleles of the newly elucidated genes and the well-known tocopherol pathway genes would greatly improve tocopherol content in maize.

A Highly Specific Genome-Wide Association Study Integrated with Transcriptome Data Reveals the Contribution of Copy Number Variations to Specialized Metabolites in Arabidopsis thaliana Accessions

Lineage-specific gene duplications contribute to a large variation in specialized metabolites among different plant species. There is also considerable variability in the specialized metabolites within a single plant species. However, it is unclear whether copy number variations (CNVs) derived from gene duplication events contribute to the diversity of specialized metabolites within species. We identified metabolome quantitative trait genes (mQTGs) associated with quantitative metabolite variations and examined the relationship between mQTGs and CNVs. We obtained 1,335 specialized metabolite signals from 53 worldwide A. thaliana accessions using liquid chromatography-quadrupole time-of-flight mass spectrometry. In this study, genes associated with specialized metabolites were inferred by either a generally authorized genome-wide association study (GWAS) approach or a novel analysis of the association between gene expression and metabolite accumulation. Genes qualified by both analyses are defined to be mQTGs. The integrated method enabled us to detect mQTGs with a low false positive rate (=5.71 × 10-4). We also identified 5,654 genes associated with 1,335 specialized metabolites. Of these genes, 4.4% were affected by CNVs, which was more than expected (χ2 test: P < 0.01). This result suggests that CNVs contribute to variations in specialized metabolites within a species. To assess the contribution of CNVs to adaptive evolution in A. thaliana, we examined the selective sweeps around the mQTGs. We observed that the mQTGs with CNVs tended to undergo selective sweeps. These observations imply that variations in specialized metabolites caused by CNVs contribute to the adaptive evolution of A. thaliana.

A eukaryotic nicotinate-inducible gene cluster: convergent evolution in fungi and bacteria

Nicotinate degradation has hitherto been elucidated only in bacteria. In the ascomycete Aspergillus nidulans, six loci, hxnS/AN9178 encoding the molybdenum cofactor-containing nicotinate hydroxylase, AN11197 encoding a Cys2/His2 zinc finger regulator HxnR, together with AN11196/hxnZ, AN11188/hxnY, AN11189/hxnP and AN9177/hxnT, are clustered and stringently co-induced by a nicotinate derivative and subject to nitrogen metabolite repression mediated by the GATA factor AreA. These genes are strictly co-regulated by HxnR. Within the hxnR gene, constitutive mutations map in two discrete regions. Aspergillus nidulans is capable of using nicotinate and its oxidation products 6-hydroxynicotinic acid and 2,5-dihydroxypyridine as sole nitrogen sources in an HxnR-dependent way. HxnS is highly similar to HxA, the canonical xanthine dehydrogenase (XDH), and has originated by gene duplication, preceding the origin of the Pezizomycotina. This cluster is conserved with some variations throughout the Aspergillaceae. Our results imply that a fungal pathway has arisen independently from bacterial ones. Significantly, the neo-functionalization of XDH into nicotinate hydroxylase has occurred independently from analogous events in bacteria. This work describes for the first time a gene cluster involved in nicotinate catabolism in a eukaryote and has relevance for the formation and evolution of co-regulated primary metabolic gene clusters and the microbial degradation of N-heterocyclic compounds.

High-quality assembly of the reference genome for scarlet sage, Salvia splendens, an economically important ornamental plant

Background

Salvia splendens Ker-Gawler, scarlet or tropical sage, is a tender herbaceous perennial widely introduced and seen in public gardens all over the world. With few molecular resources, breeding is still restricted to traditional phenotypic selection, and the genetic mechanisms underlying phenotypic variation remain unknown. Hence, a high-quality reference genome will be very valuable for marker-assisted breeding, genome editing, and molecular genetics.

Findings

We generated 66 Gb and 37 Gb of raw DNA sequences, respectively, from whole-genome sequencing of a largely homozygous scarlet sage inbred line using Pacific Biosciences (PacBio) single-molecule real-time and Illumina HiSeq sequencing platforms. The PacBio de novo assembly yielded a final genome with a scaffold N50 size of 3.12 Mb and a total length of 808 Mb. The repetitive sequences identified accounted for 57.52% of the genome sequence, and  54,008 protein-coding genes were predicted collectively with ab initio and homology-based gene prediction from the masked genome. The divergence time between S. splendens and Salvia miltiorrhiza was estimated at 28.21 million years ago (Mya). Moreover, 3,797 species-specific genes and 1,187 expanded gene families were identified for the scarlet sage genome.

Conclusions

We provide the first genome sequence and gene annotation for the scarlet sage. The availability of these resources will be of great importance for further breeding strategies, genome editing, and comparative genomics among related species.

The Vacuolar Transportome of Plant Specialized Metabolites

The plant vacuole is a cellular compartment that is essential to plant development and growth. Often plant vacuoles accumulate specialized metabolites, also called secondary metabolites, which constitute functionally and chemically diverse compounds that exert in planta many essential functions and improve the plant's fitness. These metabolites provide, for example, chemical defense against herbivorous and pathogens or chemical attractants (color and fragrance) to attract pollinators. The chemical composition of the vacuole is dynamic, and is altered during development and as a response to environmental changes. To some extent these alterations rely on vacuolar transporters, which import and export compounds into and out of the vacuole, respectively. During the past decade, significant progress was made in the identification and functional characterization of the transporters implicated in many aspects of plant specialized metabolism. Still, deciphering the molecular players underlying such processes remains a challenge for the future. In this review, we present a comprehensive summary of the most recent achievements in this field.

Annotation, classification, genomic organization and expression of the Vitis vinifera CYPome

Cytochromes P450 are enzymes that participate in a wide range of functions in plants, from hormonal signaling and biosynthesis of structural polymers, to defense or communication with other organisms. They represent one of the largest gene/protein families in the plant kingdom. The manual annotation of cytochrome P450 genes in the genome of Vitis vinifera PN40024 revealed 579 P450 sequences, including 279 complete genes. Most of the P450 sequences in grapevine genome are organized in physical clusters, resulting from tandem or segmental duplications. Although most of these clusters are small (2 to 35, median = 3), some P450 families, such as CYP76 and CYP82, underwent multiple duplications and form large clusters of homologous sequences. Analysis of gene expression revealed highly specific expression patterns, which are often the same within the genes in large physical clusters. Some of these genes are induced upon biotic stress, which points to their role in plant defense, whereas others are specifically activated during grape berry ripening and might be responsible for the production of berry-specific metabolites, such as aroma compounds. Our work provides an exhaustive and robust annotation including clear identification, structural organization, evolutionary dynamics and expression patterns for the grapevine cytochrome P450 families, paving the way to efficient functional characterization of genes involved in grapevine defense pathways and aroma biosynthesis.

Revealing the transcriptomic complexity of switchgrass by PacBio long-read sequencing

BACKGROUND

Switchgrass (Panicum virgatum L.) is an important bioenergy crop widely used for lignocellulosic research. While extensive transcriptomic analyses have been conducted on this species using short read-based sequencing techniques, very little has been reliably derived regarding alternatively spliced (AS) transcripts.

RESULTS

We present an analysis of transcriptomes of six switchgrass tissue types pooled together, sequenced using Pacific Biosciences (PacBio) single-molecular long-read technology. Our analysis identified 105,419 unique transcripts covering 43,570 known genes and 8795 previously unknown genes. 45,168 are novel transcripts of known genes. A total of 60,096 AS transcripts are identified, 45,628 being novel. We have also predicted 1549 transcripts of genes involved in cell wall construction and remodeling, 639 being novel transcripts of known cell wall genes. Most of the predicted transcripts are validated against Illumina-based short reads. Specifically, 96% of the splice junction sites in all the unique transcripts are validated by at least five Illumina reads. Comparisons between genes derived from our identified transcripts and the current genome annotation revealed that among the gene set predicted by both analyses, 16,640 have different exon-intron structures.

CONCLUSIONS

Overall, substantial amount of new information is derived from the PacBio RNA data regarding both the transcriptome and the genome of switchgrass.

Intrachromosomal colocalization strengthens co-expression, co-modification and evolutionary conservation of neighboring genes

Background

Gene order and location in chromosomes of species are non-random. Neighboring gene pairs tend to display some similarities, such as co-expression and co-modification. However, the contribution of linear proximity, spatial proximity, and evolutionary proximity to these similarities remain unclear, together with whether the presence of several types of proximity can strengthens the similarities.

Results

In this study, we investigated the properties of three kinds of colocalized gene pairs: intrachromosomal colocalized gene pairs, always-neighboring gene pairs, and evolutionary neighboring gene pairs. Our analysis showed that (1) Different types of colocalized genes differentially contribute to co-expression, co-modifications and conservation across species; (2) Intrachromosomal colocalization can strengthen co-expression and co-modification of neighboring gene pairs and their conservation across species; (3) The combination of the three kinds of colocalization can lead to the strongest co-modification and is most strongly conserved across species. (4) Colocalized gene pairs are indicative of phylogenetic relationships and whole genome duplications (WGDs).

Conclusions

These results provide valuable clues for future efforts to understand the characteristics of colocalized gene pairs and how the neighborhood affects their interactions.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4844-1) contains supplementary material, which is available to authorized users.

De novo genome assembly of Camptotheca acuminata, a natural source of the anti-cancer compound camptothecin

Camptotheca acuminata is 1 of a limited number of species that produce camptothecin, a pentacyclic quinoline alkaloid with anti-cancer activity due to its ability to inhibit DNA topoisomerase. While transcriptome studies have been performed previously with various camptothecin-producing species, no genome sequence for a camptothecin-producing species is available to date. We generated a high-quality de novo genome assembly for C. acuminata representing 403 174 860 bp on 1394 scaffolds with an N50 scaffold size of 1752 kbp. Quality assessments of the assembly revealed robust representation of the genome sequence including genic regions. Using a novel genome annotation method, we annotated 31 825 genes encoding 40 332 gene models. Based on sequence identity and orthology with validated genes from Catharanthus roseus as well as Pfam searches, we identified candidate orthologs for genes potentially involved in camptothecin biosynthesis. Extensive gene duplication including tandem duplication was widespread in the C. acuminata genome, with 2571 genes belonging to 997 tandem duplicated gene clusters. To our knowledge, this is the first genome sequence for a camptothecin-producing species, and access to the C. acuminata genome will permit not only discovery of genes encoding the camptothecin biosynthetic pathway but also reagents that can be used for heterologous expression of camptothecin and camptothecin analogs with novel pharmaceutical applications.

The effects of promoter variations of the N-Methylcanadine 1-Hydroxylase (CYP82Y1) gene on the noscapine production in opium poppy

Noscapine is an antitumor alkaloid produced in opium poppy (Papaver somniferum) and some members of the Papaveraceae family. It has been primarily used for its antitussive effects; more recently, its anticancer properties were shown. Herein, we detected an SSR embedded in the promoter region of the CYP82Y1 gene, which was found to be the first committed-step enzyme in the noscapine biosynthesis pathway, using the MISA program. Some collected ecotypes of P. somniferum were investigated for understanding of SSRs role in the regulation of gene expression and metabolite content. Quantitative PCR showed that a variation in the motif repeat number (either a decrease or increase) down-regulated the expression of the CYP82Y1 gene. Furthermore, the analysis of noscapine content suggested that a variation in the promoter region influence noscapine amount. Moreover, P. bracteatum was analyzed in both transcript and metabolite levels, and illustrated much less expression and metabolite level in comparison to P. somniferum. By exploiting the transcriptome data from the eight genera of the Papaveraceae family, we found that noscapine biosynthesis genes are present in P. bracteatum and are not shared in other genera of the Papaveraceae family. This results may explain production of a confined metabolite within a genus.

Genomic Identification and Analysis of Specialized Metabolite Biosynthetic Gene Clusters in Plants Using PlantiSMASH

Plants produce a vast diversity of specialized metabolites, which play important roles in the interactions with their microbiome, as well as with animals and other plants. Many such molecules have valuable biological activities that render them (potentially) useful as medicines, flavors and fragrances, nutritional ingredients, or cosmetics. Recently, plant scientists have discovered that the genes for many biosynthetic pathways for the production of such specialized metabolites are physically clustered on the chromosome within biosynthetic gene clusters (BGCs). The Plant Secondary Metabolite Analysis Shell (plantiSMASH) allows for the automated identification of such plant BGCs, facilitates comparison of BGCs across genomes, and helps users to predict the functional interactions of pairs of genes within and between BGCs based on coexpression analysis. In this chapter, we provide a detailed protocol on how to install and run plantiSMASH, and how to interpret its results to draw biological conclusions that are supported by the data.

The PhytoClust tool for metabolic gene clusters discovery in plant genomes

The existence of Metabolic Gene Clusters (MGCs) in plant genomes has recently raised increased interest. Thus far, MGCs were commonly identified for pathways of specialized metabolism, mostly those associated with terpene type products. For efficient identification of novel MGCs, computational approaches are essential. Here, we present PhytoClust; a tool for the detection of candidate MGCs in plant genomes. The algorithm employs a collection of enzyme families related to plant specialized metabolism, translated into hidden Markov models, to mine given genome sequences for physically co-localized metabolic enzymes. Our tool accurately identifies previously characterized plant MGCs. An exhaustive search of 31 plant genomes detected 1232 and 5531 putative gene cluster types and candidates, respectively. Clustering analysis of putative MGCs types by species reflected plant taxonomy. Furthermore, enrichment analysis revealed taxa- and species-specific enrichment of certain enzyme families in MGCs. When operating through our web-interface, PhytoClust users can mine a genome either based on a list of known cluster types or by defining new cluster rules. Moreover, for selected plant species, the output can be complemented by co-expression analysis. Altogether, we envisage PhytoClust to enhance novel MGCs discovery which will in turn impact the exploration of plant metabolism.

Tying the knot: occurrence and possible significance of gene fusions in plant metabolism and beyond

Gene fusions have recently attracted attention especially in the field of plant specialized metabolism. The occurrence of a gene fusion, in which originally separate gene products are combined into a single polypeptide, often corresponds to the functional association of individual components within a single metabolic pathway. Examples include gene fusions implicated in benzylisoquinoline alkaloid (BIA), terpenoid, and amino acid biosynthetic pathways, in which distinct domains within a fusion catalyze consecutive, yet independent reactions. Both genomic and transcriptional mechanisms result in the fusion of gene products, which can include partial or complete domain repeats and extensive domain shuffling as evident in the BIA biosynthetic enzyme norcoclaurine synthase. Artificial gene fusions are commonly deployed in attempts to engineer new or improved pathways in plants or microorganisms, based on the premise that fusions are advantageous. However, a survey of functionally characterized fusions in microbial systems shows that the functional impact of fused gene products is not straightforward. For example, whereas enzyme fusions might facilitate the metabolic channeling of unstable intermediates, this channeling can also occur between tightly associated independent enzymes. The frequent occurrence of both fused and unfused enzymes in plant and microbial metabolism adds additional complexity, in terms of both pathway functionality and evolution.

The Polyketide Components of Waxes and the Cer-cqu Gene Cluster Encoding a Novel Polyketide Synthase, the β-Diketone Synthase, DKS

The primary function of the outermost, lipophilic layer of plant aerial surfaces, called the cuticle, is preventing non-stomatal water loss. Its exterior surface is often decorated with wax crystals, imparting a blue–grey color. Identification of the barley Cer-c, -q and -u genes forming the 101 kb Cer-cqu gene cluster encoding a novel polyketide synthase—the β-diketone synthase (DKS), a lipase/carboxyl transferase, and a P450 hydroxylase, respectively, establishes a new, major pathway for the synthesis of plant waxes. The major product is a β-diketone (14,16-hentriacontane) aliphatic that forms long, thin crystalline tubes. A pathway branch leads to the formation of esterified alkan-2-ols.

Synthetic Botany

Plants are attractive platforms for synthetic biology and metabolic engineering. Plants' modular and plastic body plans, capacity for photosynthesis, extensive secondary metabolism, and agronomic systems for large-scale production make them ideal targets for genetic reprogramming. However, efforts in this area have been constrained by slow growth, long life cycles, the requirement for specialized facilities, a paucity of efficient tools for genetic manipulation, and the complexity of multicellularity. There is a need for better experimental and theoretical frameworks to understand the way genetic networks, cellular populations, and tissue-wide physical processes interact at different scales. We highlight new approaches to the DNA-based manipulation of plants and the use of advanced quantitative imaging techniques in simple plant models such as Marchantia polymorpha. These offer the prospects of improved understanding of plant dynamics and new approaches to rational engineering of plant traits.

The Genome of Medicinal Plant Macleaya cordata Provides New Insights into Benzylisoquinoline Alkaloids Metabolism

The overuse of antibiotics in animal agriculture and medicine has caused a series of potential threats to public health. Macleaya cordata is a medicinal plant species from the Papaveraceae family, providing a safe resource for the manufacture of antimicrobial feed additive for livestock. The active constituents from M. cordata are known to include benzylisoquinoline alkaloids (BIAs) such as sanguinarine (SAN) and chelerythrine (CHE), but their metabolic pathways have yet to be studied in this non-model plant. The active biosynthesis of SAN and CHE in M. cordata was first examined and confirmed by feeding ¹³C-labeled tyrosine. To gain further insights, we de novo sequenced the whole genome of M. cordata, the first to be sequenced from the Papaveraceae family. The M. cordata genome covering 378 Mb encodes 22,328 predicted protein-coding genes with 43.5% being transposable elements. As a member of basal eudicot, M. cordata genome lacks the paleohexaploidy event that occurred in almost all eudicots. From the genomics data, a complete set of 16 metabolic genes for SAN and CHE biosynthesis was retrieved, and 14 of their biochemical activities were validated. These genomics and metabolic data show the conserved BIA metabolic pathways in M. cordata and provide the knowledge foundation for future productions of SAN and CHE by crop improvement or microbial pathway reconstruction.

Integrated omics analysis of specialized metabolism in medicinal plants

Medicinal plants are a rich source of highly diverse specialized metabolites with important pharmacological properties. Until recently, plant biologists were limited in their ability to explore the biosynthetic pathways of these metabolites, mainly due to the scarcity of plant genomics resources. However, recent advances in high-throughput large-scale analytical methods have enabled plant biologists to discover biosynthetic pathways for important plant-based medicinal metabolites. The reduced cost of generating omics datasets and the development of computational tools for their analysis and integration have led to the elucidation of biosynthetic pathways of several bioactive metabolites of plant origin. These discoveries have inspired synthetic biology approaches to develop microbial systems to produce bioactive metabolites originating from plants, an alternative sustainable source of medicinally important chemicals. Since the demand for medicinal compounds are increasing with the world's population, understanding the complete biosynthesis of specialized metabolites becomes important to identify or develop reliable sources in the future. Here, we review the contributions of major omics approaches and their integration to our understanding of the biosynthetic pathways of bioactive metabolites. We briefly discuss different approaches for integrating omics datasets to extract biologically relevant knowledge and the application of omics datasets in the construction and reconstruction of metabolic models.

A bacterial artificial chromosome (BAC) genomic approach reveals partial clustering of the furanocoumarin pathway genes in parsnip

Furanocoumarins are specialized metabolites that are involved in the defense of plants against phytophagous insects. The molecular and functional characterization of the genes involved in their biosynthetic pathway is only partially complete. Many recent reports have described gene clusters responsible for the biosynthesis of specialized metabolites in plants. To investigate possible co-localization of the genes involved in the furanocoumarin pathway, we sequenced parsnip BAC clones spanning two different gene loci. We found that two genes previously identified in this pathway, CYP71AJ3 and CYP71AJ4, were located on the same BAC, whereas a third gene, PsPT1, belonged to a different BAC clone. Chromosome mapping using fluorescence in situ hybridization (FISH) indicated that PsPT1 and the CYP71AJ3-CYP71AJ4 clusters are located on two different chromosomes. Sequencing the BAC clone harboring PsPT1 led to the identification of a gene encoding an Fe(II) α-ketoglutarate-dependent dioxygenase (PsDIOX) situated in the neighborhood of PsPT1 and confirmed the occurrence of a second gene cluster involved in the furanocoumarin pathway. This enzyme metabolizes p-coumaroyl CoA, leading exclusively to the synthesis of umbelliferone, an important intermediate compound in furanocoumarin synthesis. This work provides an insight into the genomic organization of genes from the furanocoumarin biosynthesis pathway organized in more than one gene cluster. It also confirms that the screening of a genomic library and the sequencing of BAC clones represent a valuable tool to identify genes involved in biosynthetic pathways dedicated to specialized metabolite synthesis.

Functional Characterization of CYP716 Family P450 Enzymes in Triterpenoid Biosynthesis in Tomato

Triterpenoids are a group of structurally diverse specialized metabolites that frequently show useful bioactivities. These chemicals are biosynthesized from the common precursor 2,3-oxidosqualene in plants. The carbon skeletons produced by oxidosqualene cyclase (OSC) are usually modified by cytochrome P450 monooxygenases (P450s) and UDP-dependent glycosyltransferases. These biosynthetic enzymes contribute to the structural diversification of plant triterpenoids. Until now, many P450 enzymes have been characterized as triterpenoid oxidases. Among them, the CYP716 family P450 enzymes, which have been isolated from a wide range of plant families, seem to contribute to the triterpenoid structural diversification. Many CYP716 family P450 enzymes have been characterized as the multifunctional triterpene C-28 oxidases, which oxidize α-amyrin and β-amyrin to the widely distributed triterpenoids ursolic and oleanolic acids, respectively. Tomato (Solanum lycopersicum) is one of the most important solanaceous crops in the world. However, little information is known regarding its triterpenoid biosynthesis. To understand the mechanism of triterpenoid biosynthesis in tomato, we focused on the function of CYP716 family enzymes as triterpenoid oxidases. We isolated all six CYP716 family genes from the Micro-Tom cultivar of tomato, and functionally characterized them in the heterologous yeast expression system. The in vivo enzymatic assays showed that CYP716A44 and CYP716A46 exhibited the ordinary C-28 oxidation activity against α-amyrin and β-amyrin to produce ursolic and oleanolic acids, respectively. Interestingly, one CYP716E subfamily enzyme, CYP716E26, exhibited the previously unreported C-6β hydroxylation activity against β-amyrin to produce a rare bioactive triterpenoid, daturadiol (olean-12-ene-3β,6β-diol). To determine the roles of the CYP716 family genes in tomato triterpenoid biosynthesis, we analyzed the gene expression and triterpenoid accumulation patterns in different plant tissues by performing the quantitative real-time polymerase chain reaction (qPCR) and gas chromatography-mass spectrometry (GC-MS) analyses, respectively. High levels of the CYP716A44 gene expression and the accumulation of C-28-oxidized triterpenoids, ursolic acid, and oleanolic acid were observed in the roots, indicating a significant contribution of the CYP716A44 gene in the triterpenoid biosynthesis in tomato. Thus, our study partially elucidated the mechanism of triterpenoid biosynthesis in tomato, and identified CYP716E26 as a novel C-6β hydroxylase for its subsequent use in the combinatorial biosynthesis of bioactive triterpenoids.

Natural products from filamentous fungi and production by heterologous expression

Filamentous fungi represent an incredibly rich and rather overlooked reservoir of natural products, which often show potent bioactivity and find applications in different fields. Increasing the naturally low yields of bioactive metabolites within their host producers can be problematic, and yield improvement is further hampered by such fungi often being genetic intractable or having demanding culturing conditions. Additionally, total synthesis does not always represent a cost-effective approach for producing bioactive fungal-inspired metabolites, especially when pursuing assembly of compounds with complex chemistry. This review aims at providing insights into heterologous production of secondary metabolites from filamentous fungi, which has been established as a potent system for the biosynthesis of bioactive compounds. Numerous advantages are associated with this technique, such as the availability of tools that allow enhanced production yields and directing biosynthesis towards analogues of the naturally occurring metabolite. Furthermore, a choice of hosts is available for heterologous expression, going from model unicellular organisms to well-characterised filamentous fungi, which has also been shown to allow the study of biosynthesis of complex secondary metabolites. Looking to the future, fungi are likely to continue to play a substantial role as sources of new pharmaceuticals and agrochemicals-either as producers of novel natural products or indeed as platforms to generate new compounds through synthetic biology.

Draft genome assembly and annotation of Glycyrrhiza uralensis, a medicinal legume

Chinese liquorice/licorice (Glycyrrhiza uralensis) is a leguminous plant species whose roots and rhizomes have been widely used as a herbal medicine and natural sweetener. Whole-genome sequencing is essential for gene discovery studies and molecular breeding in liquorice. Here, we report a draft assembly of the approximately 379-Mb whole-genome sequence of strain 308-19 of G. uralensis; this assembly contains 34 445 predicted protein-coding genes. Comparative analyses suggested well-conserved genomic components and collinearity of gene loci (synteny) between the genome of liquorice and those of other legumes such as Medicago and chickpea. We observed that three genes involved in isoflavonoid biosynthesis, namely, 2-hydroxyisoflavanone synthase (CYP93C), 2,7,4'-trihydroxyisoflavanone 4'-O-methyltransferase/isoflavone 4'-O-methyltransferase (HI4OMT) and isoflavone-7-O-methyltransferase (7-IOMT) formed a cluster on the scaffold of the liquorice genome and showed conserved microsynteny with Medicago and chickpea. Based on the liquorice genome annotation, we predicted genes in the P450 and UDP-dependent glycosyltransferase (UGT) superfamilies, some of which are involved in triterpenoid saponin biosynthesis, and characterised their gene expression with the reference genome sequence. The genome sequencing and its annotations provide an essential resource for liquorice improvement through molecular breeding and the discovery of useful genes for engineering bioactive components through synthetic biology approaches.

The biosynthetic pathway of the nonsugar, high-intensity sweetener mogroside V from Siraitia grosvenorii

The consumption of sweeteners, natural as well as synthetic sugars, is implicated in an array of modern-day health problems. Therefore, natural nonsugar sweeteners are of increasing interest. We identify here the biosynthetic pathway of the sweet triterpenoid glycoside mogroside V, which has a sweetening strength of 250 times that of sucrose and is derived from mature fruit of luo-han-guo (Siraitia grosvenorii, monk fruit). A whole-genome sequencing of Siraitia, leading to a preliminary draft of the genome, was combined with an extensive transcriptomic analysis of developing fruit. A functional expression survey of nearly 200 candidate genes identified the members of the five enzyme families responsible for the synthesis of mogroside V: squalene epoxidases, triterpenoid synthases, epoxide hydrolases, cytochrome P450s, and UDP-glucosyltransferases. Protein modeling and docking studies corroborated the experimentally proven functional enzyme activities and indicated the order of the metabolic steps in the pathway. A comparison of the genomic organization and expression patterns of these Siraitia genes with the orthologs of other Cucurbitaceae implicates a strikingly coordinated expression of the pathway in the evolution of this species-specific and valuable metabolic pathway. The genomic organization of the pathway genes, syntenously preserved among the Cucurbitaceae, indicates, on the other hand, that gene clustering cannot account for this novel secondary metabolic pathway.