Metabolite‐based genome‐wide association study enables dissection of the flavonoid decoration pathway of wheat kernels

The marriage of metabolomic approaches with genetic design has proven a powerful tool in dissecting diversity in the metabolome and has additionally enhanced our understanding of complex traits. That said, such studies have rarely been carried out in wheat. In this study, we detected 805 metabolites from wheat kernels and profiled their relative contents among 182 wheat accessions, conducting a metabolite‐based genome‐wide association study (mGWAS) utilizing 14 646 previously described polymorphic SNP markers. A total of 1098 mGWAS associations were detected with large effects, within which 26 candidate genes were tentatively designated for 42 loci. Enzymatic assay of two candidates indicated they could catalyse glucosylation and subsequent malonylation of various flavonoids and thereby the major flavonoid decoration pathway of wheat kernel was dissected. Moreover, numerous high‐confidence genes associated with metabolite contents have been provided, as well as more subdivided metabolite networks which are yet to be explored within our data. These combined efforts presented the first step towards realizing metabolomics‐associated breeding of wheat.     

Metabolomics analysis and metabolite‐agronomic trait associations using kernels of wheat (Triticum aestivum) recombinant inbred lines

Plants produce numerous metabolites that are important for their development and growth. However, the genetic architecture of the wheat metabolome has not been well studied. Here, utilizing a high‐density genetic map, we conducted a comprehensive metabolome study via widely targeted LC‐MS/MS to analyze the wheat kernel metabolism. We further combined agronomic traits and dissected the genetic relationship between metabolites and agronomic traits. In total, 1260 metabolic features were detected. Using linkage analysis, 1005 metabolic quantitative trait loci (mQTLs) were found distributed unevenly across the genome. Twenty‐four candidate genes were found to modulate the levels of different metabolites, of which two were functionally annotated by in vitro analysis to be involved in the synthesis and modification of flavonoids. Combining the correlation analysis of metabolite‐agronomic traits with the co‐localization of methylation quantitative trait locus (mQTL) and phenotypic QTL (pQTL), genetic relationships between the metabolites and agronomic traits were uncovered. For example, a candidate was identified using correlation and co‐localization analysis that may manage auxin accumulation, thereby affecting number of grains per spike (NGPS). Furthermore, metabolomics data were used to predict the performance of wheat agronomic traits, with metabolites being found that provide strong predictive power for NGPS and plant height. This study used metabolomics and association analysis to better understand the genetic basis of the wheat metabolism which will ultimately assist in wheat breeding.     

Trichome regulator SlMIXTA‐like directly manipulates primary metabolism in tomato fruit

Trichomes are storage compartments for specialized metabolites in many plant species. In trichome, plant primary metabolism is significantly changed, providing substrates for downstream secondary metabolism. However, little is known of how plants coordinate trichome formation and primary metabolism regulation. In this report, tomato (Solanum lycopersicum) trichome regulator SlMIXTA‐like is indicated as a metabolic regulation gene by mGWAS analysis. Overexpression of SlMIXTA‐like in tomato fruit enhances trichome formation. In addition, SlMIXTA‐like can directly bind to the promoter region of gene encoding 3‐deoxy‐7‐phosphoheptulonate synthase (SlDAHPS) to activate its expression. Induction of SlDAHPS expression enhances shikimate pathway activities and provides substrates for downstream secondary metabolism. Our data provide direct evidence that trichome regulator can directly manipulate primary metabolism, in which way plants can coordinate metabolic regulation and the formation of storage compartments for specialized metabolites. The newly identified SlMIXTA‐like can be used for future metabolic engineering.     

Metabolome‐genome‐wide association study dissects genetic architecture for generating natural variation in rice secondary metabolism

Plants produce structurally diverse secondary (specialized) metabolites to increase their fitness for survival under adverse environments. Several bioactive compounds for new drugs have been identified through screening of plant extracts. In this study, genome‐wide association studies (GWAS) were conducted to investigate the genetic architecture behind the natural variation of rice secondary metabolites. GWAS using the metabolome data of 175 rice accessions successfully identified 323 associations among 143 single nucleotide polymorphisms (SNPs) and 89 metabolites. The data analysis highlighted that levels of many metabolites are tightly associated with a small number of strong quantitative trait loci (QTLs). The tight association may be a mechanism generating strains with distinct metabolic composition through the crossing of two different strains. The results indicate that one plant species produces more diverse phytochemicals than previously expected, and plants still contain many useful compounds for human applications.     

Three‐dimensional genetic networks among seed oil‐related traits,metabolites and genes reveal the genetic foundations of oil synthesis in soybean

Although the biochemical and genetic basis of lipid metabolism is clear in Arabidopsis, there is limited information concerning the relevant genes in Glycine max (soybean). To address this issue, we constructed three‐dimensional genetic networks using six seed oil‐related traits, 52 lipid metabolism‐related metabolites and 54 294 SNPs in 286 soybean accessions in total. As a result, 284 and 279 candidate genes were found to be significantly associated with seed oil‐related traits and metabolites by phenotypic and metabolic genome‐wide association studies and multi‐omics analyses, respectively. Using minimax concave penalty (MCP) and smoothly clipped absolute deviation (SCAD) analyses, six seed oil‐related traits were found to be significantly related to 31 metabolites. Among the above candidate genes, 36 genes were found to be associated with oil synthesis (27 genes), amino acid synthesis (four genes) and the tricarboxylic acid (TCA) cycle (five genes), and four genes (GmFATB1a, GmPDAT, GmPLDα1 and GmDAGAT1) are already known to be related to oil synthesis. Using this information, 133 three‐dimensional genetic networks were constructed, 24 of which are known, e.g. pyruvate–GmPDAT–GmFATA2–oil content. Using these networks, GmPDAT, GmAGT and GmACP4 reveal the genetic relationships between pyruvate and the three major nutrients, and GmPDAT, GmZF351 and GmPgs1 reveal the genetic relationships between amino acids and seed oil content. In addition, GmCds1, along with average temperature in July and the rainfall from June to September, influence seed oil content across years. This study provides a new approach for the construction of three‐dimensional genetic networks and reveals new information for soybean seed oil improvement and the identification of gene function.     

Meta‐Omics‐ and Metabolic Modeling‐Assisted Deciphering of Human Microbiota Metabolism

The human microbiota is a complex community of commensal, symbiotic, and pathogenic microbes that play a crucial role in maintaining the homeostasis of human health. Such a homeostasis is maintained through the collective functioning of enzymatic genes responsible for the production of metabolites, enabling the interaction and signaling within microbiota as well as between microbes and the human host. Understanding microbial genes, their associated chemistries and functions would be valuable for engineering systemic metabolic pathways within the microbiota to manage human health and diseases. Given that there are many unknown gene metabolic functions and interactions, increasing efforts have been made to gain insights into the underlying functions of microbiota metabolism. This can be achieved through culture‐independent metagenomic approaches and metabolic modeling to simulate the microenvironment of human microbiota. In this article, the recent advances in metagenome mining and functional profiling for the discovery of the genetic and biochemical links in human microbiota metabolism as well as metabolic modeling for simulation and prediction of metabolic fluxes in the human microbiota are reviewed. This review provides useful insights into the understanding, reconstruction, and modulation of the human microbiota guided by the knowledge acquired from the basic understanding of the human microbiota metabolism.     

Genome‐scale gene/reaction essentiality and synthetic lethality analysis

Synthetic lethals are to pairs of non‐essential genes whose simultaneous deletion prohibits growth. One can extend the concept of synthetic lethality by considering gene groups of increasing size where only the simultaneous elimination of all genes is lethal, whereas individual gene deletions are not. We developed optimization‐based procedures for the exhaustive and targeted enumeration of multi‐gene (and by extension multi‐reaction) lethals for genome‐scale metabolic models. Specifically, these approaches are applied to iAF1260, the latest model of Escherichia coli, leading to the complete identification of all double and triple gene and reaction synthetic lethals as well as the targeted identification of quadruples and some higher‐order ones. Graph representations of these synthetic lethals reveal a variety of motifs ranging from hub‐like to highly connected subgraphs providing a birds‐eye view of the avenues available for redirecting metabolism and uncovering complex patterns of gene utilization and interdependence. The procedure also enables the use of falsely predicted synthetic lethals for metabolic model curation. By analyzing the functional classifications of the genes involved in synthetic lethals, we reveal surprising connections within and across clusters of orthologous group functional classifications.     

A Critical Review of Genome Editing and Synthetic Biology Applications in Metabolic Engineering of Microalgae and Cyanobacteria

Being the green gold of the future, microalgae and cyanobacteria have recently attracted considerable interest worldwide, for their metabolites such as lipids, protein, pigments, and bioactive compounds have immense potential for sustainable energy and pharmaceutical production capabilities. In the last decades, the efforts attended to enhance the usage of microalgae and cyanobacteria by genetic manipulation, synthetic and metabolic engineering. However, the development of photoautotrophic cell factories have rarely compared to the heterotrophic counterparts due to limited tools, bioinformatics, and multi‐omics database. Therefore, recent advances of their genome editing techniques by clustered regularly interspaced short palindromic repeats (CRISPR) technology, and potential applications of their metabolic engineering and regulation approaches are examined in this review. Moreover, the contemporary achievements of synthetic biology approaches of microalgae and cyanobacteria in carbon fixation and sequestration, lipid and triacylglycerol (TAG), and sustainable production of high value‐added chemicals, such as carotenoids and docosahexaenoic acid (DHA), have been also discussed. From recent genomic study to trends in metabolic regulation of microalgae and cyanobacteria and a comprehensive assessment of the current challenges and opportunities for microalgae and cyanobacteria is also conducted.     

Genetics of the human metabolome,what is next?

Increases in throughput and decreases in costs have facilitated large scale metabolomics studies, the simultaneous measurement of large numbers of biochemical components in biological samples. Initial large scale studies focused on biomarker discovery for disease or disease progression and helped to understand biochemical pathways underlying disease. The first population-based studies that combined metabolomics and genome wide association studies (mGWAS) have increased our understanding of the (genetic) regulation of biochemical conversions. Measurements of metabolites as intermediate phenotypes are a potentially very powerful approach to uncover how genetic variation affects disease susceptibility and progression. However, we still face many hurdles in the interpretation of mGWAS data. Due to the composite nature of many metabolites, single enzymes may affect the levels of multiple metabolites and, conversely, levels of single metabolites may be affected by multiple enzymes. Here, we will provide a global review of the current status of mGWAS. We will specifically discuss the application of prior biological knowledge present in databases to the interpretation of mGWAS results and discuss the potential of mathematical models. As the technology continuously improves to detect metabolites and to measure genetic variation, it is clear that comprehensive systems biology based approaches are required to further our insight in the association between genes, metabolites and disease. This article is part of a Special Issue entitled: From Genome to Function.     

Comparative analysis of metabolome of rice seeds at three developmental stages using a recombinant inbred line population

Plants are considered an important food and nutrition source for humans. Despite advances in plant seed metabolomics, knowledge about the genetic and molecular bases of rice seed metabolomes at different developmental stages is still limited. Here, using Zhenshan 97 (ZS97) and Minghui 63 (MH63), we performed a widely targeted metabolic profiling in seeds during grain filling, mature seeds and germinating seeds. The diversity between MH63 and ZS97 was characterized in terms of the content of metabolites and the metabolic shifting across developmental stages. Taking advantage of the ultra‐high‐density genetic map of a population of 210 recombinant inbred lines (RILs) derived from a cross between ZS97 and MH63, we identified 4681 putative metabolic quantitative trait loci (mQTLs) in seeds across the three stages. Further analysis of the mQTLs for the codetected metabolites across the three stages revealed that the genetic regulation of metabolite accumulation was closely related to developmental stage. Using in silico analyses, we characterized 35 candidate genes responsible for 30 structurally identified or annotated compounds, among which LOC_Os07g04970 and LOC_Os06g03990 were identified to be responsible for feruloylserotonin and l ‐asparagine content variation across populations, respectively. Metabolite?agronomic trait association and colocation between mQTLs and phenotypic quantitative trait loci (pQTLs) revealed the complexity of the metabolite?agronomic trait relationship and the corresponding genetic basis.     

Extreme‐phenotype genome‐wide association study (XP‐GWAS): a method for identifying trait‐associated variants by sequencing pools of individuals selected from a diversity panel

Although approaches for performing genome‐wide association studies (GWAS) are well developed, conventional GWAS requires high‐density genotyping of large numbers of individuals from a diversity panel. Here we report a method for performing GWAS that does not require genotyping of large numbers of individuals. Instead XP‐GWAS (extreme‐phenotype GWAS) relies on genotyping pools of individuals from a diversity panel that have extreme phenotypes. This analysis measures allele frequencies in the extreme pools, enabling discovery of associations between genetic variants and traits of interest. This method was evaluated in maize (Zea mays) using the well‐characterized kernel row number trait, which was selected to enable comparisons between the results of XP‐GWAS and conventional GWAS. An exome‐sequencing strategy was used to focus sequencing resources on genes and their flanking regions. A total of 0.94 million variants were identified and served as evaluation markers; comparisons among pools showed that 145 of these variants were statistically associated with the kernel row number phenotype. These trait‐associated variants were significantly enriched in regions identified by conventional GWAS. XP‐GWAS was able to resolve several linked QTL and detect trait‐associated variants within a single gene under a QTL peak. XP‐GWAS is expected to be particularly valuable for detecting genes or alleles responsible for quantitative variation in species for which extensive genotyping resources are not available, such as wild progenitors of crops, orphan crops, and other poorly characterized species such as those of ecological interest.     

Genome-scale models of bacterial metabolism: reconstruction and applications

Genome-scale metabolic models bridge the gap between genome-derived biochemical information and metabolic phenotypes in a principled manner, providing a solid interpretative framework for experimental data related to metabolic states, and enabling simple in silico experiments with whole-cell metabolism. Models have been reconstructed for almost 20 bacterial species, so far mainly through expert curation efforts integrating information from the literature with genome annotation. A wide variety of computational methods exploiting metabolic models have been developed and applied to bacteria, yielding valuable insights into bacterial metabolism and evolution, and providing a sound basis for computer-assisted design in metabolic engineering. Recent advances in computational systems biology and high-throughput experimental technologies pave the way for the systematic reconstruction of metabolic models from genomes of new species, and a corresponding expansion of the scope of their applications. In this review, we provide an introduction to the key ideas of metabolic modeling, survey the methods, and resources that enable model reconstruction and refinement, and chart applications to the investigation of global properties of metabolic systems, the interpretation of experimental results, and the re-engineering of their biochemical capabilities.     

Linkage disequilibrium clustering‐based approach for association mapping with tightly linked genomewide data

Genomewide association studies (GWAS) aim to identify genetic markers strongly associated with quantitative traits by utilizing linkage disequilibrium (LD) between candidate genes and markers. However, because of LD between nearby genetic markers, the standard GWAS approaches typically detect a number of correlated SNPs covering long genomic regions, making corrections for multiple testing overly conservative. Additionally, the high dimensionality of modern GWAS data poses considerable challenges for GWAS procedures such as permutation tests, which are computationally intensive. We propose a cluster‐based GWAS approach that first divides the genome into many large nonoverlapping windows and uses linkage disequilibrium network analysis in combination with principal component (PC) analysis as dimensional reduction tools to summarize the SNP data to independent PCs within clusters of loci connected by high LD. We then introduce single‐ and multilocus models that can efficiently conduct the association tests on such high‐dimensional data. The methods can be adapted to different model structures and used to analyse samples collected from the wild or from biparental F₂ populations, which are commonly used in ecological genetics mapping studies. We demonstrate the performance of our approaches with two publicly available data sets from a plant (Arabidopsis thaliana) and a fish (Pungitius pungitius), as well as with simulated data.