The complex genetic architecture of the metabolome

OsLSC6 Regulates Leaf Sheath Color and Cold Tolerance in Rice Revealed by Metabolite Genome Wide Association Study

Plant metabolites including anthocyanins play an important role in the growth of plants, as well as in regulating biotic and abiotic stress responses to the environment. Here we report comprehensive profiling of 3315 metabolites and a further metabolic-based genome-wide association study (mGWAS) based on 292,485 SNPs obtained from 311 rice accessions, including 160 wild and 151 cultivars. We identified hundreds of common variants affecting a large number of secondary metabolites with large effects at high throughput. Finally, we identified a novel gene namely OsLSC6 (Oryza sativa leaf sheath color 6), which encoded a UDP 3-O-glucosyltransferase and involved in the anthocyanin biosynthesis of Cyanidin-3-Galc (sd1825) responsible for leaf sheath color, and resulted in significant different accumulation of sd1825 between wild (purple) and cultivars (green). The results of knockout transgenic experiments showed that OsLSC6 regulated the biosynthesis and accumulation of sd1825, controlled the purple leaf sheath. Our further research revealed that OsLSC6 also confers resistance to cold stress during the seedling stage in rice. And we identified that a SNP in OsLSC6 was responsible for the leaf sheath color and chilling tolerance, supporting the importance of OsLSC6 in plant adaption. Our study could not only demonstrate that OsLSC6 is a vital regulator during anthocyanin biosynthesis and abiotic stress responses, but also provide a powerful complementary tool based on metabolites-to-genes analysis by mGWAS for functional gene identification andpromising candidate in future rice breeding and improvement.

Genome-wide association studies identify loci controlling specialized seed metabolites in Arabidopsis

Plants synthesize specialized metabolites to facilitate environmental and ecological interactions. During evolution, plants diversified in their potential to synthesize these metabolites. Quantitative differences in metabolite levels of natural Arabidopsis (Arabidopsis thaliana) accessions can be employed to unravel the genetic basis for metabolic traits using genome-wide association studies (GWAS). Here, we performed metabolic GWAS on seeds of a panel of 315 A. thaliana natural accessions, including the reference genotypes C24 and Col-0, for polar and semi-polar seed metabolites using untargeted ultra-performance liquid chromatography-mass spectrometry. As a complementary approach, we performed quantitative trait locus (QTL) mapping of near-isogenic introgression lines between C24 and Col-0 for specific seed specialized metabolites. Besides common QTL between seeds and leaves, GWAS revealed seed-specific QTL for specialized metabolites, indicating differences in the genetic architecture of seeds and leaves. In seeds, aliphatic methylsulfinylalkyl and methylthioalkyl glucosinolates associated with the ALKENYL HYDROXYALKYL PRODUCING loci (GS-ALK and GS-OHP) on chromosome 4 containing alkenyl hydroxyalkyl producing 2 (AOP2) and 3 (AOP3) or with the GS-ELONG locus on chromosome 5 containing methylthioalkyl malate synthase (MAM1) and MAM3. We detected two unknown sulfur-containing compounds that were also mapped to these loci. In GWAS, some of the annotated flavonoids (kaempferol 3-O-rhamnoside-7-O-rhamnoside, quercetin 3-O-rhamnoside-7-O-rhamnoside) were mapped to transparent testa 7 (AT5G07990), encoding a cytochrome P450 75B1 monooxygenase. Three additional mass signals corresponding to quercetin-containing flavonols were mapped to UGT78D2 (AT5G17050). The association of the loci and associating metabolic features were functionally verified in knockdown mutant lines. By performing GWAS and QTL mapping, we were able to leverage variation of natural populations and parental lines to study seed specialized metabolism. The GWAS data set generated here is a high-quality resource that can be investigated in further studies.

Genetic control of abiotic stress-related specialized metabolites in sunflower

BACKGROUND

Abiotic stresses in plants include all the environmental conditions that significantly reduce yields, like drought and heat. One of the most significant effects they exert at the cellular level is the accumulation of reactive oxygen species, which cause extensive damage. Plants possess two mechanisms to counter these molecules, i.e. detoxifying enzymes and non-enzymatic antioxidants, which include many classes of specialized metabolites. Sunflower, the fourth global oilseed, is considered moderately drought resistant. Abiotic stress tolerance in this crop has been studied using many approaches, but the control of specialized metabolites in this context remains poorly understood. Here, we performed the first genome-wide association study using abiotic stress-related specialized metabolites as molecular phenotypes in sunflower. After analyzing leaf specialized metabolites of 450 hybrids using liquid chromatography-mass spectrometry, we selected a subset of these compounds based on their association with previously known abiotic stress-related quantitative trait loci. Eventually, we characterized these molecules and their associated genes.

RESULTS

We putatively annotated 30 compounds which co-localized with abiotic stress-related quantitative trait loci and which were associated to seven most likely candidate genes. A large proportion of these compounds were potential antioxidants, which was in agreement with the role of specialized metabolites in abiotic stresses. The seven associated most likely candidate genes, instead, mainly belonged to cytochromes P450 and glycosyltransferases, two large superfamilies which catalyze greatly diverse reactions and create a wide variety of chemical modifications. This was consistent with the high plasticity of specialized metabolism in plants.

CONCLUSIONS

This is the first characterization of the genetic control of abiotic stress-related specialized metabolites in sunflower. By providing hints concerning the importance of antioxidant molecules in this biological context, and by highlighting some of the potential molecular mechanisms underlying their biosynthesis, it could pave the way for novel applications in breeding. Although further analyses will be required to better understand this topic, studying how antioxidants contribute to the tolerance to abiotic stresses in sunflower appears as a promising area of research.

Unlocking secrets of nature's chemists: Potential of CRISPR/Cas-based tools in plant metabolic engineering for customized nutraceutical and medicinal profiles

Plant species have evolved diverse metabolic pathways to effectively respond to internal and external signals throughout their life cycle, allowing adaptation to their sessile and phototropic nature. These pathways selectively activate specific metabolic processes, producing plant secondary metabolites (PSMs) governed by genetic and environmental factors. Humans have utilized PSM-enriched plant sources for millennia in medicine and nutraceuticals. Recent technological advances have significantly contributed to discovering metabolic pathways and related genes involved in the biosynthesis of specific PSM in different food crops and medicinal plants. Consequently, there is a growing demand for plant materials rich in nutrients and bioactive compounds, marketed as "superfoods". To meet the industrial demand for superfoods and therapeutic PSMs, modern methods such as system biology, omics, synthetic biology, and genome editing (GE) play a crucial role in identifying the molecular players, limiting steps, and regulatory circuitry involved in PSM production. Among these methods, clustered regularly interspaced short palindromic repeats-CRISPR associated protein (CRISPR/Cas) is the most widely used system for plant GE due to its simple design, flexibility, precision, and multiplexing capabilities. Utilizing the CRISPR-based toolbox for metabolic engineering (ME) offers an ideal solution for developing plants with tailored preventive (nutraceuticals) and curative (therapeutic) metabolic profiles in an ecofriendly way. This review discusses recent advances in understanding the multifactorial regulation of metabolic pathways, the application of CRISPR-based tools for plant ME, and the potential research areas for enhancing plant metabolic profiles.

Metabolome-Based Genome-Wide Association Study of Duck Meat Leads to Novel Genetic and Biochemical Insights

Meat is among the most consumed foods worldwide and has a unique flavor and high nutrient density in the human diet. However, the genetic and biochemical bases of meat nutrition and flavor are poorly understood. Here, 3431 metabolites and 702 volatiles in 423 skeletal muscle samples are profiled from a gradient consanguinity segregating population generated by Pekin duck × Liancheng duck crosses using metabolomic approaches. The authors identified 2862 metabolome-based genome-wide association studies (mGWAS) signals and 48 candidate genes potentially modulating metabolite and volatile levels, 79.2% of which are regulated by cis-regulatory elements. The level of plasmalogen is significantly associated with TMEM189 encoding plasmanylethanolamine desaturase 1. The levels of 2-pyrrolidone and glycerophospholipids are regulated by the gene expression of AOX1 and ACBD5, which further affects the levels of volatiles, 2-pyrrolidone and decanal, respectively. Genetic variations in GADL1 and CARNMT2 determine the levels of 49 metabolites including L-carnosine and anserine. This study provides novel insights into the genetic and biochemical basis of skeletal muscle metabolism and constitutes a valuable resource for the precise improvement of meat nutrition and flavor.

Limitations and advantages of using metabolite-based genome-wide association studies: focus on fruit quality traits

In the last decades, linkage mapping has help in the location of metabolite quantitative trait loci (QTL) in many species; however, this approach shows some limitations. Recently, thanks to the most recent advanced in high-throughput genotyping technologies like next-generation sequencing, metabolite genome-wide association study (mGWAS) has been proposed a powerful tool to identify the genetic variants in polygenic agrinomic traits. Fruit flavor is a complex interaction of aroma volatiles and taste being sugar and acid ratio key parameter for flavor acceptance. Here, we review recent progress of mGWAS in pinpoint gene polymorphisms related to flavor-related metabolites in fruits. Despite clear successes in discovering novel genes or regions associated with metabolite accumulation affecting sensory attributes in fruits, GWAS incurs in several limitations summarized in this review. In addition, in our own work, we performed mGWAS on 194 Citrus grandis accessions to investigate the genetic control of individual primary and lipid metabolites in ripe fruit. We have identified a total of 667 associations for 14 primary metabolites including amino acids, sugars, and organic acids, and 768 associations corresponding to 47 lipids. Furthermore, candidate genes related to important metabolites related to fruit quality such as sugars, organic acids and lipids were discovered.

Natural variation of respiration-related traits in plants

Plant respiration is one of the greatest global metabolic fluxes, but rates of respiration vary massively both within different cell types as well as between different individuals and different species. Whilst this is well known, few studies have detailed population-level variation of respiration until recently. The last 20 years have seen a renaissance in studies of natural variance. In this review, we describe how experimental breeding populations and collections of large populations of accessions can be used to determine the genetic architecture of plant traits. We further detail how these approaches have been used to study the rate of respiration per se as well as traits that are intimately associated with respiration. The review highlights specific breakthroughs in these areas but also concludes that the approach should be more widely adopted in the study of respiration per se as opposed to the more frequently studied respiration-related traits.

Plant and microbial sciences as key drivers in the development of metabolomics research

This year marks the 25th anniversary of the coinage of the term metabolome [S. G. Oliver et al., Trends Biotech. 16, 373-378 (1998)]. As the field rapidly advances, it is important to take stock of the progress which has been made to best inform the disciplines future. While a medical-centric perspective on metabolomics has recently been published [M. Giera et al., Cell Metab. 34, 21-34 (2022)], this largely ignores the pioneering contributions made by the plant and microbial science communities. In this perspective, we provide a contemporary overview of all fields in which metabolomics is employed with particular emphasis on both methodological and application breakthroughs made in plant and microbial sciences that have shaped this evolving research discipline from the very early days of its establishment. This will not cover all types of metabolomics assays currently employed but will focus mainly on those utilizing mass spectrometry-based measurements since they are currently by far the most prominent. Having established the historical context of metabolomics, we will address the key challenges currently facing metabolomics and offer potential approaches by which these can be faced. Most salient among these is the fact that the vast majority of mass features are as yet not annotated with high confidence; what we may refer to as definitive identification. We discuss the potential of both standard compound libraries and artificial intelligence technologies to address this challenge and the use of natural variance-based approaches such as genome-wide association studies in attempt to assign specific functions to the myriad of structurally similar and complex specialized metabolites. We conclude by stating our contention that as these challenges are epic and that they will need far greater cooperative efforts from biologists, chemists, and computer scientists with an interest in all kingdoms of life than have been made to date. Ultimately, a better linkage of metabolome and genome data will likely also be needed particularly considering the Earth BioGenome Project.

Genome-wide association studies of five free amino acid levels in rice

Rice (Oryza sativa L.) is one of the important staple foods for human consumption and livestock use. As a complex quality trait, free amino acid (FAA) content in rice is of nutritional importance. To dissect the genetic mechanism of FAA level, five amino acids' (Val, Leu, Ile, Arg, and Trp) content and 4,325,832 high-quality SNPs of 448 rice accessions were used to conduct genome-wide association studies (GWAS) with nine different methods. Of these methods, one single-locus method (GEMMA), seven multi-locus methods (mrMLM, pLARmEB, FASTmrEMMA, pKWmEB, FASTmrMLM, ISIS EM-BLASSO, and FarmCPU), and the recent released 3VmrMLM were adopted for methodological comparison of quantitative trait nucleotide (QTN) detection and identification of stable quantitative trait nucleotide loci (QTLs). As a result, 987 QTNs were identified by eight multi-locus GWAS methods; FASTmrEMMA detected the most QTNs (245), followed by 3VmrMLM (160), and GEMMA detected the least QTNs (0). Among 88 stable QTLs identified by the above methods, 3VmrMLM has some advantages, such as the most common QTNs, the highest LOD score, and the highest proportion of all detected stable QTLs. Around these stable QTLs, candidate genes were found in the GO classification to be involved in the primary metabolic process, biosynthetic process, and catalytic activity, and shown in KEGG analysis to have participated in metabolic pathways, biosynthesis of amino acids, and tryptophan metabolism. Natural variations of candidate genes resulting in the content alteration of five FAAs were identified in this association panel. In addition, 95 QTN-by-environment interactions (QEIs) of five FAA levels were detected by 3VmrMLM only. GO classification showed that the candidate genes got involved in the primary metabolic process, transport, and catalytic activity. Candidate genes of QEIs played important roles in valine, leucine, and isoleucine degradation (QEI_09_03978551 and candidate gene LOC_Os09g07830 in the Leu dataset), tryptophan metabolism (QEI_01_00617184 and candidate gene LOC_Os01g02020 in the Trp dataset), and glutathione metabolism (QEI_12_09153839 and candidate gene LOC_Os12g16200 in the Arg dataset) pathways through KEGG analysis. As an alternative of the multi-locus GWAS method, these findings suggested that the application of 3VmrMLM may provide new insights into better understanding FAA accumulation and facilitate the molecular breeding of rice with high FAA level.

Metabolomics-Driven Mining of Metabolite Resources: Applications and Prospects for Improving Vegetable Crops

Vegetable crops possess a prominent nutri-metabolite pool that not only contributes to the crop performance in the fields, but also offers nutritional security for humans. In the pursuit of identifying, quantifying and functionally characterizing the cellular metabolome pool, biomolecule separation technologies, data acquisition platforms, chemical libraries, bioinformatics tools, databases and visualization techniques have come to play significant role. High-throughput metabolomics unravels structurally diverse nutrition-rich metabolites and their entangled interactions in vegetable plants. It has helped to link identified phytometabolites with unique phenotypic traits, nutri-functional characters, defense mechanisms and crop productivity. In this study, we explore mining diverse metabolites, localizing cellular metabolic pathways, classifying functional biomolecules and establishing linkages between metabolic fluxes and genomic regulations, using comprehensive metabolomics deciphers of the plant’s performance in the environment. We discuss exemplary reports covering the implications of metabolomics, addressing metabolic changes in vegetable plants during crop domestication, stage-dependent growth, fruit development, nutri-metabolic capabilities, climatic impacts, plant-microbe-pest interactions and anthropogenic activities. Efforts leading to identify biomarker metabolites, candidate proteins and the genes responsible for plant health, defense mechanisms and nutri-rich crop produce are documented. With the insights on metabolite-QTL (mQTL) driven genetic architecture, molecular breeding in vegetable crops can be revolutionized for developing better nutritional capabilities, improved tolerance against diseases/pests and enhanced climate resilience in plants.

Multi-omics analyses of 398 foxtail millet accessions reveal genomic regions associated with domestication, metabolite traits, and anti-inflammatory effects

Foxtail millet (Setaria italica), which was domesticated from the wild species green foxtail (Setaria viridis), is a rich source of phytonutrients for humans. To evaluate how breeding changed the metabolome of foxtail millet grains, we generated and analyzed the datasets encompassing the genomes, transcriptomes, metabolomes, and anti-inflammatory indices from 398 foxtail millet accessions. We identified hundreds of common variants that influence numerous secondary metabolites. We observed tremendous differences in natural variations of the metabolites and their underlying genetic architectures between distinct sub-groups of foxtail millet. Furthermore, we found that the selection of the gene alleles associated with yellow grains led to altered profiles of metabolites such as carotenoids and endogenous phytohormones. Using CRISPR-mediated genome editing we validated the function of PHYTOENE SYNTHASE 1 (PSY1) gene in affecting millet grain color and quality. Interestingly, our in vitro cell inflammation assays showed that 83 metabolites in millet grains have anti-inflammatory effects. Taken together, our multi-omics study illustrates how the breeding history of foxtail millet has shaped its metabolite profile. The datasets we generated in this study also provide important resources for further understanding how millet grain quality is affected by different metabolites, laying the foundations for future millet genetic research and metabolome-assisted improvement.

Genome-Wide Association Study With Growth-Related Traits and Secondary Metabolite Contents in Red- and White-Heart Chinese Fir

Chinese fir [Cunninghamia lanceolata (Lamb.) Hook] is an important evergreen coniferous tree species that is widely distributed in many southern provinces of China and has important economic value. The Chinese fir accounts for 1/4 and 1/3 of the total artificial forest area and stock volume, respectively. Red-heart Chinese fir is popular in the market because of its high density and red heartwood. The long-growth cycle hindered the breeding process of Chinese fir, while molecular marker-assisted breeding could accelerate it. However, Chinese fir, a perennial conifer species, has a large genome, which has not yet been published. In this study, the growth-related traits and secondary metabolite contents of red- and white-heart Chinese fir were measured and found to be different between them. There are extremely significant differences among growth-related traits (p < 0.001), but secondary metabolite contents have different correlations due to differences in chemical structure. Moreover, genotype effect analysis of the substantially correlated single nucleotide polymorphisms (SNPs) revealed that most of the loci related to each growth-related traits were different from each other, indicating a type specificity of the genes regulated different growth-related traits. Furthermore, among the loci related to secondary metabolite contents, nine loci associated with multiple metabolite phenotypes such as Marker21022_4, Marker21022_172, Marker24559_31, Marker27425_37, Marker20748_85, Marker18841_115, Marker18841_198, Marker65846_146, and Marker21486_163, suggesting the presence of pleiotropic genes. This study identified the potential SNP markers associated with secondary metabolites in Chinese fir, thus setting the basis for molecular marker-assisted selection.

Generalizable approaches for genomic prediction of metabolites in plants

Plant metabolites are important traits for plant breeders seeking to improve nutrition and agronomic performance yet integrating selection for metabolomic traits can be limited by phenotyping expense and degree of genetic characterization, especially of uncommon metabolites. As such, developing generalizable genomic selection methods based on biochemical pathway biology for metabolites that are transferable across plant populations would benefit plant breeding programs. We tested genomic prediction accuracy for >600 metabolites measured by gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS) in oat (Avena sativa L.) seed. Using a discovery germplasm panel, we conducted metabolite genome-wide association study (mGWAS) and selected loci to use in multikernel models that encompassed metabolome-wide mGWAS results or mGWAS from specific metabolite structures or biosynthetic pathways. Metabolite kernels developed from LC-MS metabolites in the discovery panel improved prediction accuracy of LC-MS metabolite traits in the validation panel consisting of more advanced breeding lines. No approach, however, improved prediction accuracy for GC-MS metabolites. We ranked model performance by metabolite and found that metabolites with similar polarity had consistent rankings of models. Overall, testing biological rationales for developing kernels for genomic prediction across populations contributes to developing frameworks for plant breeding for metabolite traits.

A Perspective on Plant Phenomics: Coupling Deep Learning and Near-Infrared Spectroscopy

The trait-based approach in plant ecology aims at understanding and classifying the diversity of ecological strategies by comparing plant morphology and physiology across organisms. The major drawback of the approach is that the time and financial cost of measuring the traits on many individuals and environments can be prohibitive. We show that combining near-infrared spectroscopy (NIRS) with deep learning resolves this limitation by quickly, non-destructively, and accurately measuring a suite of traits, including plant morphology, chemistry, and metabolism. Such an approach also allows to position plants within the well-known CSR triangle that depicts the diversity of plant ecological strategies. The processing of NIRS through deep learning identifies the effect of growth conditions on trait values, an issue that plagues traditional statistical approaches. Together, the coupling of NIRS and deep learning is a promising high-throughput approach to capture a range of ecological information on plant diversity and functioning and can accelerate the creation of extensive trait databases.

Computational prediction of plant metabolic pathways

Uncovering genes encoding enzymes responsible for the biosynthesis of diverse plant metabolites is essential for metabolic engineering and production of plant metabolite-derived medicine. With the availability of multi-omics data for an ever-increasing number of plant species and the development of computational approaches, the metabolic pathways of many important plant compounds can be predicted, complementing a more traditional genetic and/or biochemical approach. Here, we summarize recent progress in predicting plant metabolic pathways using genome, transcriptome, proteome, interactome, and/or metabolome data, and the utility of integrating these data with machine learning to further improve metabolic pathway predictions.

Methylated Cytochrome P450 and the Solute Carrier Family of Genes Correlate With Perturbations in Bile Acid Metabolism in Parkinson’s Disease

Parkinson’s disease (PD) is second most prevalent neurodegenerative disorder following Alzheimer’s disease. Parkinson’s disease is hypothesized to be caused by a multifaceted interplay between genetic and environmental factors. Herein, and for the first time, we describe the integration of metabolomics and epigenetics (genome-wide DNA methylation; epimetabolomics) to profile the frontal lobe from people who died from PD and compared them with age-, and sex-matched controls. We identified 48 metabolites to be at significantly different concentrations (FDR q < 0.05), 4,313 differentially methylated sites [5’-C-phosphate-G-3’ (CpGs)] (FDR q < 0.05) and increased DNA methylation age in the primary motor cortex of people who died from PD. We identified Primary bile acid biosynthesis as the major biochemical pathway to be perturbed in the frontal lobe of PD sufferers, and the metabolite taurine (p-value = 5.91E-06) as being positively correlated with CpG cg14286187 (SLC25A27; CYP39A1) (FDR q = 0.002), highlighting previously unreported biochemical changes associated with PD pathogenesis. In this novel multi-omics study, we identify regulatory mechanisms which we believe warrant future translational investigation and central biomarkers of PD which require further validation in more accessible biomatrices.

Metabolite Diversity and Metabolic Genome-Wide Marker Association Studies (Mgwas) for Health Benefiting Nutritional Traits in Pearl Millet Grains

As efforts are made to increase food security, millets are gaining increasing importance due to their excellent nutritional credentials. Among the millets, pearl millet is the predominant species possessing several health benefiting nutritional traits in its grain that are helpful in mitigating chronic illnesses such as type−2 diabetes and obesity. In this paper, we conducted metabolomic fingerprinting of 197 pearl millet inbred lines drawn randomly from within the world collection of pearl millet germplasm and report the extent of genetic variation for health benefitting metabolites in these genotypes. Metabolites were extracted from seeds and assessed using flow infusion high-resolution mass spectrometry (FIE-HRMS). Metabolite features (m/z), whose levels significantly differed among the germplasm inbred lines, were identified by ANOVA corrected for FDR and subjected to functional pathway analysis. A number of health-benefiting metabolites linked to dietary starch, antioxidants, vitamins, and lipid metabolism-related compounds were identified. Metabolic genome-wide association analysis (mGWAS) performed using the 396 m/z as phenotypic traits and the 76 K SNP as genotypic variants identified a total of 897 SNPs associated with health benefiting nutritional metabolite at the -log p-value ≤ 4.0. From these associations, 738 probable candidate genes were predicted to have an important role in starch, antioxidants, vitamins, and lipid metabolism. The mGWAS analysis focused on genes involved in starch branching (α-amylase, β-amylase), vitamin-K reductase, UDP-glucuronosyl, and UDP-glucosyl transferase (UGTs), L-ascorbate oxidase, and isoflavone 2′-monooxygenase genes, which are known to be linked to increases in human health benefiting metabolites. We demonstrate how metabolomic, genomic, and statistical approaches can be utilized to pinpoint genetic variations and their functions linked to key nutritional properties in pearl millet, which in turn can be bred into millets and other cereals crops using plant breeding methods.

Genome-wide association studies: assessing trait characteristics in model and crop plants

GWAS involves testing genetic variants across the genomes of many individuals of a population to identify genotype–phenotype association. It was initially developed and has proven highly successful in human disease genetics. In plants genome-wide association studies (GWAS) initially focused on single feature polymorphism and recombination and linkage disequilibrium but has now been embraced by a plethora of different disciplines with several thousand studies being published in model and crop species within the last decade or so. Here we will provide a comprehensive review of these studies providing cases studies on biotic resistance, abiotic tolerance, yield associated traits, and metabolic composition. We also detail current strategies of candidate gene validation as well as the functional study of haplotypes. Furthermore, we provide a critical evaluation of the GWAS strategy and its alternatives as well as future perspectives that are emerging with the emergence of pan-genomic datasets.

Prediction of complex phenotypes using the Drosophila melanogaster metabolome

Understanding the genotype-phenotype map and how variation at different levels of biological organization is associated are central topics in modern biology. Fast developments in sequencing technologies and other molecular omic tools enable researchers to obtain detailed information on variation at DNA level and on intermediate endophenotypes, such as RNA, proteins and metabolites. This can facilitate our understanding of the link between genotypes and molecular and functional organismal phenotypes. Here, we use the Drosophila melanogaster Genetic Reference Panel and nuclear magnetic resonance (NMR) metabolomics to investigate the ability of the metabolome to predict organismal phenotypes. We performed NMR metabolomics on four replicate pools of male flies from each of 170 different isogenic lines. Our results show that metabolite profiles are variable among the investigated lines and that this variation is highly heritable. Second, we identify genes associated with metabolome variation. Third, using the metabolome gave better prediction accuracies than genomic information for four of five quantitative traits analyzed. Our comprehensive characterization of population-scale diversity of metabolomes and its genetic basis illustrates that metabolites have large potential as predictors of organismal phenotypes. This finding is of great importance, e.g., in human medicine, evolutionary biology and animal and plant breeding.

Understanding Omics Driven Plant Improvement and de novo Crop Domestication: Some Examples

In the current era, one of biggest challenges is to shorten the breeding cycle for rapid generation of a new crop variety having high yield capacity, disease resistance, high nutrient content, etc. Advances in the "-omics" technology have revolutionized the discovery of genes and bio-molecules with remarkable precision, resulting in significant development of plant-focused metabolic databases and resources. Metabolomics has been widely used in several model plants and crop species to examine metabolic drift and changes in metabolic composition during various developmental stages and in response to stimuli. Over the last few decades, these efforts have resulted in a significantly improved understanding of the metabolic pathways of plants through identification of several unknown intermediates. This has assisted in developing several new metabolically engineered important crops with desirable agronomic traits, and has facilitated the de novo domestication of new crops for sustainable agriculture and food security. In this review, we discuss how "omics" technologies, particularly metabolomics, has enhanced our understanding of important traits and allowed speedy domestication of novel crop plants.

The genetics underlying metabolic signatures in a brown rice diversity panel and their vital role in human nutrition

Brown rice (Oryza sativa) possesses various nutritionally dense bioactive phytochemicals exhibiting a wide range of antioxidant, anti-cancer, and anti-diabetic properties known to promote various human health benefits. However, despite the wide claims made about the importance of brown rice for human nutrition the underlying metabolic diversity has not been systematically explored. Non-targeted metabolite profiling of developing and mature seeds of a diverse genetic panel of 320 rice cultivars allowed quantification of 117 metabolites. The metabolite genome-wide association study (mGWAS) detected genetic variants influencing diverse metabolic targets in developing and mature seeds. We further interlinked genetic variants on chromosome 7 (6.06-6.43 Mb region) with complex epistatic genetic interactions impacting multi-dimensional nutritional targets, including complex carbohydrate starch quality, the glycemic index, antioxidant catechin, and rice grain color. Through this nutrigenomics approach rare gene bank accessions possessing genetic variants in bHLH and IPT5 genes were identified through haplotype enrichment. These variants were associated with a low glycemic index, higher catechin levels, elevated total flavonoid contents, and heightened antioxidant activity in the whole grain with elevated anti-cancer properties being confirmed in cancer cell lines. This multi-disciplinary nutrigenomics approach thus allowed us to discover the genetic basis of human health-conferring diversity in the metabolome of brown rice.

Translating insights from the seed metabolome into improved prediction for lipid-composition traits in oat (Avena sativa L.)

Oat (Avena sativa L.) seed is a rich resource of beneficial lipids, soluble fiber, protein, and antioxidants, and is considered a healthful food for humans. Little is known regarding the genetic controllers of variation for these compounds in oat seed. We characterized natural variation in the mature seed metabolome using untargeted metabolomics on 367 diverse lines and leveraged this information to improve prediction for seed quality traits. We used a latent factor approach to define unobserved variables that may drive covariance among metabolites. One hundred latent factors were identified, of which 21% were enriched for compounds associated with lipid metabolism. Through a combination of whole-genome regression and association mapping, we show that latent factors that generate covariance for many metabolites tend to have a complex genetic architecture. Nonetheless, we recovered significant associations for 23% of the latent factors. These associations were used to inform a multi-kernel genomic prediction model, which was used to predict seed lipid and protein traits in two independent studies. Predictions for 8 of the 12 traits were significantly improved compared to genomic best linear unbiased prediction when this prediction model was informed using associations from lipid-enriched factors. This study provides new insights into variation in the oat seed metabolome and provides genomic resources for breeders to improve selection for health-promoting seed quality traits. More broadly, we outline an approach to distill high-dimensional "omics" data to a set of biologically meaningful variables and translate inferences on these data into improved breeding decisions.

Genome-wide association study on metabolite accumulation in a wild barley NAM population reveals natural variation in sugar metabolism

Metabolites play a key role in plants as they are routing plant developmental processes and are involved in biotic and abiotic stress responses. Their analysis can offer important information on the underlying processes. Regarding plant breeding, metabolite concentrations can be used as biomarkers instead of or in addition to genetic markers to predict important phenotypic traits (metabolic prediction). In this study, we applied a genome-wide association study (GWAS) in a wild barley nested association mapping (NAM) population to identify metabolic quantitative trait loci (mQTL). A set of approximately 130 metabolites, measured at early and late sampling dates, was analysed. For four metabolites from the early and six metabolites from the late sampling date significant mQTL (grouped as 19 mQTL for the early and 25 mQTL for the late sampling date) were found. Interestingly, all of those metabolites could be classified as sugars. Sugars are known to be involved in signalling, plant growth and plant development. Sugar-related genes, encoding mainly sugar transporters, have been identified as candidate genes for most of the mQTL. Moreover, several of them co-localized with known flowering time genes like Ppd-H1, HvELF3, Vrn-H1, Vrn-H2 and Vrn-H3, hinting on the known role of sugars in flowering. Furthermore, numerous disease resistance-related genes were detected, pointing to the signalling function of sugars in plant resistance. An mQTL on chromosome 1H in the region of 13 Mbp to 20 Mbp stood out, that alone explained up to 65% of the phenotypic variation of a single metabolite. Analysis of family-specific effects within the diverse NAM population showed the available natural genetic variation regarding sugar metabolites due to different wild alleles. The study represents a step towards a better understanding of the genetic components of metabolite accumulation, especially sugars, thereby linking them to biological functions in barley.

New insights into the role of MADS-box transcription factor gene CmANR1 on root and shoot development in chrysanthemum (Chrysanthemum morifolium)

BACKGROUND

MADS-box transcription factors (TFs) are the key regulators of multiple developmental processes in plants; among them, a chrysanthemum MADS-box TF CmANR1 has been isolated and described as functioning in root development in response to high nitrate concentration signals. However, how CmANR1 affects root and shoot development remains unclear.

RESULTS

We report that CmANR1 plays a positive role in root system development in chrysanthemum throughout the developmental stages of in vitro tissue cultures. Metabolomics combined with transcriptomics assays show that CmANR1 promotes robust root system development by facilitating nitrate assimilation, and influencing the metabolic pathways of amino acid, glycolysis, and the tricarboxylic acid cycle (TCA) cycle. Also, we found that the expression levels of TFs associated with the nitrate signaling pathways, such as AGL8, AGL21, and LBD29, are significantly up-regulated in CmANR1-transgenic plants relative to the wild-type (WT) control; by contrast, the expression levels of RHD3-LIKE, LBD37, and GATA23 were significantly down-regulated. These results suggest that these nitrate signaling associated TFs are involved in CmANR1-modulated control of root development. In addition, CmANR1 also acts as a positive regulator to control shoot growth and development.

CONCLUSIONS

These findings provide potential mechanisms of MADS-box TF CmANR1 modulation of root and shoot development, which occurs by regulating a series of nitrate signaling associated TFs, and influencing the metabolic pathways of amino acid and glycolysis, as well as TCA cycle and nitrate assimilation.

Allelic composition of carotenoid metabolic genes in 13 founders influences carotenoid composition in juice sac tissues of fruits among Japanese citrus breeding population

To enrich carotenoids, especially β-cryptoxanthin, in juice sac tissues of fruits via molecular breeding in citrus, allele mining was utilized to dissect allelic variation of carotenoid metabolic genes and identify an optimum allele on the target loci characterized by expression quantitative trait (eQTL) analysis. SNPs of target carotenoid metabolic genes in 13 founders of the Japanese citrus breeding population were explored using the SureSelect target enrichment method. An independent allele was determined based on the presence or absence of reliable SNPs, using trio analysis to confirm inheritability between parent and offspring. Among the 13 founders, there were 7 PSY alleles, 7 HYb alleles, 11 ZEP alleles, 5 NCED alleles, and 4 alleles for the eQTL that control the transcription levels of PDS and ZDS among the ancestral species, indicating that some founders acquired those alleles from them. The carotenoid composition data of 263 breeding pedigrees in juice sac tissues revealed that the phenotypic variance of carotenoid composition was similar to that in the 13 founders, whereas the mean of total carotenoid content increased. This increase in total carotenoid content correlated with the increase in either or both β-cryptoxanthin and violaxanthin in juice sac tissues. Bayesian statistical analysis between allelic composition of target genes and carotenoid composition in 263 breeding pedigrees indicated that PSY-a and ZEP-e alleles at PSY and ZEP loci had strong positive effects on increasing the total carotenoid content, including β-cryptoxanthin and violaxanthin, in juice sac tissues. Moreover, the pyramiding of these alleles also increased the β-cryptoxanthin content. Interestingly, the offset interaction between the alleles with increasing and decreasing effects on carotenoid content and the epistatic interaction among carotenoid metabolic genes were observed and these interactions complexed carotenoid profiles in breeding population. These results revealed that allele composition would highly influence the carotenoid composition in citrus fruits. The allelic genotype information for the examined carotenoid metabolic genes in major citrus varieties and the trio-tagged SNPs to discriminate the optimum alleles (PSY-a and ZEP-e) from the rest would promise citrus breeders carotenoid enrichment in fruit via molecular breeding.

Genomic Heritability and Genome-Wide Association Studies of Plasma Metabolites in Crossbred Beef Cattle

Metabolites, substrates or products of metabolic processes, are involved in many biological functions, such as energy metabolism, signaling, stimulatory and inhibitory effects on enzymes and immunological defense. Metabolomic phenotypes are influenced by combination of genetic and environmental effects allowing for metabolome-genome-wide association studies (mGWAS) as a powerful tool to investigate the relationship between these phenotypes and genetic variants. The objectives of this study were to estimate genomic heritability and perform mGWAS and in silico functional enrichment analyses for a set of plasma metabolites in Canadian crossbred beef cattle. Thirty-three plasma metabolites and 45,266 single nucleotide polymorphisms (SNPs) were available for 475 animals. Genomic heritability for all metabolites was estimated using genomic best linear unbiased prediction (GBLUP) including genomic breed composition as covariates in the model. A single-step GBLUP implemented in BLUPF90 programs was used to determine SNP P values and the proportion of genetic variance explained by SNP windows containing 10 consecutive SNPs. The top 10 SNP windows that explained the largest genetic variation for each metabolite were identified and mapped to detect corresponding candidate genes. Functional enrichment analyses were performed on metabolites and their candidate genes using the Ingenuity Pathway Analysis software. Eleven metabolites showed low to moderate heritability that ranged from 0.09 ± 0.15 to 0.36 ± 0.15, while heritability estimates for 22 metabolites were zero or negligible. This result indicates that while variations in 11 metabolites were due to genetic variants, the majority are largely influenced by environment. Three significant SNP associations were detected for betaine (rs109862186), L-alanine (rs81117935), and L-lactic acid (rs42009425) based on Bonferroni correction for multiple testing (family wise error rate <0.05). The SNP rs81117935 was found to be located within the Catenin Alpha 2 gene (CTNNA2) showing a possible association with the regulation of L-alanine concentration. Other candidate genes were identified based on additive genetic variance explained by SNP windows of 10 consecutive SNPs. The observed heritability estimates and the candidate genes and networks identified in this study will serve as baseline information for research into the utilization of plasma metabolites for genetic improvement of crossbred beef cattle.

Dissection of the domestication-shaped genetic architecture of lettuce primary metabolism

Lettuce (Lactuca sativa L.) is an important vegetable crop species worldwide. The primary metabolism of this species is essential for its growth, development and reproduction as well as providing a considerable direct source of energy and nutrition for humans. Here, through investigating 77 primary metabolites in 189 accessions including all major horticultural types and wild lettuce L. serriola we showed that the metabolites in L. serriola were different from those in cultivated lettuce. The findings were consistent with the demographic model of lettuce and supported a single domestication event for this species. Selection signals among these metabolic traits were detected. Specifically, galactinol, malate, quinate and threonate were significantly affected by the domestication process and cultivar differentiation of lettuce. Galactinol and raffinose might have been selected during stem lettuce cultivation as an adaption to the local environments in China. Furthermore, we identified 154 loci significantly associated with the level of 51 primary metabolites. Three genes (LG8749721, LG8763094 and LG5482522) responsible for the levels of galactinol, raffinose, quinate and chlorogenic acid were further dissected, which may have been the target of domestication and/or affected by local adaptation. Additionally, our findings strongly suggest that human selection resulted in reduced quinate and chlorogenic acid levels in cultivated lettuce. Our study thus provides beneficial genetic resources for lettuce quality improvement and sheds light on the domestication and evolution of this important leafy green.

Metabolomics as an Emerging Tool in the Search for Astrobiologically Relevant Biomarkers

It is now routinely possible to sequence and recover microbial genomes from environmental samples. To the degree it is feasible to assign transcriptional and translational functions to these genomes, it should be possible, in principle, to largely understand the complete molecular inputs and outputs of a microbial community. However, gene-based tools alone are presently insufficient to describe the full suite of chemical reactions and small molecules that compose a living cell. Metabolomic tools have developed quickly and now enable rapid detection and identification of small molecules within biological and environmental samples. The convergence of these technologies will soon facilitate the detection of novel enzymatic activities, novel organisms, and potentially extraterrestrial life-forms on solar system bodies. This review explores the methodological problems and scientific opportunities facing researchers who hope to apply metabolomic methods in astrobiology-related fields, and how present challenges might be overcome.

Assessment of metabolic variability and diversity present in leaf, peel and pulp tissue of diploid and triploid Musa spp

Banana (Musa spp.) plants produce many health promoting compounds in leaf, peel and pulp. For a robust metabolic analysis of these tissues, leaf at five developmental stages were compared to assess suitable sampling practices. Results confirmed that the common sampling practise of leaf 3 is applicable for metabolic comparisons. The developed work flow was applied to analyse the metabolite diversity present in 18 different Musa varieties, providing baseline levels of metabolites in leaf, peel and pulp tissue. Correlation analysis was then used to ascertain whether similar trends can be detected in the three plant tissues of the diversity panel. The genome group displayed a dominant role in the composition of the metabolome in all three tissues. This led to the conclusion that a correlation between tissues was only possible within a genome group as the different parental backgrounds caused too great a variation in the metabolomes. It also suggests the metabolome could be used to monitor the interaction/hybridisation of genomes during breeding programmes.

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.

An efficient multi-locus mixed model framework for the detection of small and linked QTLs in F2

In the genetic system that regulates complex traits, metabolites, gene expression levels, RNA editing levels and DNA methylation, a series of small and linked genes exist. To date, however, little is known about how to design an efficient framework for the detection of these kinds of genes. In this article, we propose a genome-wide composite interval mapping (GCIM) in F₂. First, controlling polygenic background via selecting markers in the genome scanning of linkage analysis was replaced by estimating polygenic variance in a genome-wide association study. This can control large, middle and minor polygenic backgrounds in genome scanning. Then, additive and dominant effects for each putative quantitative trait locus (QTL) were separately scanned so that a negative logarithm P-value curve against genome position could be separately obtained for each kind of effect. In each curve, all the peaks were identified as potential QTLs. Thus, almost all the small-effect and linked QTLs are included in a multi-locus model. Finally, adaptive least absolute shrinkage and selection operator (adaptive lasso) was used to estimate all the effects in the multi-locus model, and all the nonzero effects were further identified by likelihood ratio test for true QTL identification. This method was used to reanalyze four rice traits. Among 25 known genes detected in this study, 16 small-effect genes were identified only by GCIM. To further demonstrate GCIM, a series of Monte Carlo simulation experiments was performed. As a result, GCIM is demonstrated to be more powerful than the widely used methods for the detection of closely linked and small-effect QTLs.

Genetic and metabolomic architecture of variation in diet restriction-mediated lifespan extension in Drosophila

In most organisms, dietary restriction (DR) increases lifespan. However, several studies have found that genotypes within the same species vary widely in how they respond to DR. To explore the mechanisms underlying this variation, we exposed 178 inbred Drosophila melanogaster lines to a DR or ad libitum (AL) diet, and measured a panel of 105 metabolites under both diets. Twenty four out of 105 metabolites were associated with the magnitude of the lifespan response. These included proteinogenic amino acids and metabolites involved in α-ketoglutarate (α-KG)/glutamine metabolism. We confirm the role of α-KG/glutamine synthesis pathways in the DR response through genetic manipulations. We used covariance network analysis to investigate diet-dependent interactions between metabolites, identifying the essential amino acids threonine and arginine as “hub” metabolites in the DR response. Finally, we employ a novel metabolic and genetic bipartite network analysis to reveal multiple genes that influence DR lifespan response, some of which have not previously been implicated in DR regulation. One of these is CCHa2R, a gene that encodes a neuropeptide receptor that influences satiety response and insulin signaling. Across the lines, variation in an intronic single nucleotide variant of CCHa2R correlated with variation in levels of five metabolites, all of which in turn were correlated with DR lifespan response. Inhibition of adult CCHa2R expression extended DR lifespan of flies, confirming the role of CCHa2R in lifespan response. These results provide support for the power of combined genomic and metabolomic analysis to identify key pathways underlying variation in this complex quantitative trait.

Dietary restriction extends lifespan across most organisms in which it has been tested. However, several studies have now demonstrated that this effect can vary dramatically across different genotypes within a population. Within a population, dietary restriction might be beneficial for some, yet detrimental for others. Here, we measure the metabolome of 178 genetically characterized fly strains on fully fed and restricted diets. The fly strains vary widely in their lifespan response to dietary restriction. We then use information about each strain’s genome and metabolome (a measure of small molecules circulating in flies) to pinpoint cellular pathways that govern this variation in response. We identify a novel pathway involving the gene CCHa2R, which encodes a neuropeptide receptor that has not previously been implicated in dietary restriction or age-related signaling pathways. This study demonstrates the power of leveraging systems biology and network biology methods to understand how and why different individuals vary in their response to health and lifespan-extending interventions.

PANOMICS meets germplasm

Genotyping-by-sequencing has enabled approaches for genomic selection to improve yield, stress resistance and nutritional value. More and more resource studies are emerging providing 1000 and more genotypes and millions of SNPs for one species covering a hitherto inaccessible intraspecific genetic variation. The larger the databases are growing, the better statistical approaches for genomic selection will be available. However, there are clear limitations on the statistical but also on the biological part. Intraspecific genetic variation is able to explain a high proportion of the phenotypes, but a large part of phenotypic plasticity also stems from environmentally driven transcriptional, post-transcriptional, translational, post-translational, epigenetic and metabolic regulation. Moreover, regulation of the same gene can have different phenotypic outputs in different environments. Consequently, to explain and understand environment-dependent phenotypic plasticity based on the available genotype variation we have to integrate the analysis of further molecular levels reflecting the complete information flow from the gene to metabolism to phenotype. Interestingly, metabolomics platforms are already more cost-effective than NGS platforms and are decisive for the prediction of nutritional value or stress resistance. Here, we propose three fundamental pillars for future breeding strategies in the framework of Green Systems Biology: (i) combining genome selection with environment-dependent PANOMICS analysis and deep learning to improve prediction accuracy for marker-dependent trait performance; (ii) PANOMICS resolution at subtissue, cellular and subcellular level provides information about fundamental functions of selected markers; (iii) combining PANOMICS with genome editing and speed breeding tools to accelerate and enhance large-scale functional validation of trait-specific precision breeding.

PANOMICS meets germplasm

Genotyping-by-sequencing has enabled approaches for genomic selection to improve yield, stress resistance and nutritional value. More and more resource studies are emerging providing 1000 and more genotypes and millions of SNPs for one species covering a hitherto inaccessible intraspecific genetic variation. The larger the databases are growing, the better statistical approaches for genomic selection will be available. However, there are clear limitations on the statistical but also on the biological part. Intraspecific genetic variation is able to explain a high proportion of the phenotypes, but a large part of phenotypic plasticity also stems from environmentally driven transcriptional, post-transcriptional, translational, post-translational, epigenetic and metabolic regulation. Moreover, regulation of the same gene can have different phenotypic outputs in different environments. Consequently, to explain and understand environment-dependent phenotypic plasticity based on the available genotype variation we have to integrate the analysis of further molecular levels reflecting the complete information flow from the gene to metabolism to phenotype. Interestingly, metabolomics platforms are already more cost-effective than NGS platforms and are decisive for the prediction of nutritional value or stress resistance. Here, we propose three fundamental pillars for future breeding strategies in the framework of Green Systems Biology: (i) combining genome selection with environment-dependent PANOMICS analysis and deep learning to improve prediction accuracy for marker-dependent trait performance; (ii) PANOMICS resolution at subtissue, cellular and subcellular level provides information about fundamental functions of selected markers; (iii) combining PANOMICS with genome editing and speed breeding tools to accelerate and enhance large-scale functional validation of trait-specific precision breeding.

Can metabolic prediction be an alternative to genomic prediction in barley?

Like other crop species, barley, the fourth most important crop worldwide, suffers from the genetic bottleneck effect, where further improvements in performance through classical breeding methods become difficult. Therefore, indirect selection methods are of great interest. Here, genomic prediction (GP) based on 33,005 SNP markers and, alternatively, metabolic prediction (MP) based on 128 metabolites with sampling at two different time points in one year, were applied to predict multi-year agronomic traits in the nested association mapping (NAM) population HEB-25. We found prediction abilities of up to 0.93 for plant height with SNP markers and of up to 0.61 for flowering time with metabolites. Interestingly, prediction abilities in GP increased after reducing the number of incorporated SNP markers. The estimated effects of GP and MP were highly concordant, indicating MP as an interesting alternative to GP, being able to reflect a stable genotype-specific metabolite profile. In MP, sampling at an early developmental stage outperformed sampling at a later stage. The results confirm the value of GP for future breeding. With MP, an interesting alternative was also applied successfully. However, based on our results, usage of MP alone cannot be recommended in barley. Nevertheless, MP can assist in unravelling physiological pathways for the expression of agronomically important traits.

Pathogen Genetic Control of Transcriptome Variation in the Arabidopsis thaliana - Botrytis cinerea Pathosystem

In plant-pathogen relations, disease symptoms arise from the interaction of the host and pathogen genomes. Host-pathogen functional gene interactions are well described, whereas little is known about how the pathogen genetic variation modulates both organisms' transcriptomes. To model and generate hypotheses on a generalist pathogen control of gene expression regulation, we used the Arabidopsis thaliana-Botrytis cinerea pathosystem and the genetic diversity of a collection of 96 B. cinerea isolates. We performed expression-based genome-wide association (eGWA) for each of 23,947 measurable transcripts in Arabidopsis (host), and 9267 measurable transcripts in B. cinerea (pathogen). Unlike other eGWA studies, we detected a relative absence of locally acting expression quantitative trait loci (cis-eQTL), partly caused by structural variants and allelic heterogeneity hindering their identification. This study identified several distantly acting trans-eQTL linked to eQTL hotspots dispersed across Botrytis genome that altered only Botrytis transcripts, only Arabidopsis transcripts, or transcripts from both species. Gene membership in the trans-eQTL hotspots suggests links between gene expression regulation and both known and novel virulence mechanisms in this pathosystem. Genes annotated to these hotspots provide potential targets for blocking manipulation of the host response by this ubiquitous generalist necrotrophic pathogen.

Identifying genetic variants underlying phenotypic variation in plants without complete genomes

Structural variants and presence/absence polymorphisms are common in plant genomes, yet they are routinely overlooked in genome-wide association studies (GWAS). Here, we expand the type of genetic variants detected in GWAS to include major deletions, insertions and rearrangements. We first use raw sequencing data directly to derive short sequences, k-mers, that mark a broad range of polymorphisms independently of a reference genome. We then link k-mers associated with phenotypes to specific genomic regions. Using this approach, we reanalyzed 2,000 traits in Arabidopsis thaliana, tomato and maize populations. Associations identified with k-mers recapitulate those found with SNPs, but with stronger statistical support. Importantly, we discovered new associations with structural variants and with regions missing from reference genomes. Our results demonstrate the power of performing GWAS before linking sequence reads to specific genomic regions, which allows the detection of a wider range of genetic variants responsible for phenotypic variation.

Using networks to identify and interpret natural variation

Studies on natural variation and network biology inherently work to summarize vast amounts of information and data. The combination of these two areas of study while creating datasets of immense complexity is critical to their mutual progress. Networks are necessary as a way to work to reduce the dimensionality inherent in natural variation with 100 s to 1000 s of genotypes. Correspondingly natural variation is essential for testing how networks may or may not be shared across individuals or species. Advances in this area of cross-fertilization including using networks directly as phenotypes and the use of networks to help in prioritizing candidate gene validation efforts. Interesting new observations on frequent presence-absence variation in gene content and adaptation is beginning to highlight the potential for natural variation in network presence-absence. This review attempts to delve into these new insights.

Metabolite database for root, tuber, and banana crops to facilitate modern breeding in understudied crops

Roots, tubers, and bananas (RTB) are vital staples for food security in the world's poorest nations. A major constraint to current RTB breeding programmes is limited knowledge on the available diversity due to lack of efficient germplasm characterization and structure. In recent years large-scale efforts have begun to elucidate the genetic and phenotypic diversity of germplasm collections and populations and, yet, biochemical measurements have often been overlooked despite metabolite composition being directly associated with agronomic and consumer traits. Here we present a compound database and concentration range for metabolites detected in the major RTB crops: banana (Musa spp.), cassava (Manihot esculenta), potato (Solanum tuberosum), sweet potato (Ipomoea batatas), and yam (Dioscorea spp.), following metabolomics-based diversity screening of global collections held within the CGIAR institutes. The dataset including 711 chemical features provides a valuable resource regarding the comparative biochemical composition of each RTB crop and highlights the potential diversity available for incorporation into crop improvement programmes. Particularly, the tropical crops cassava, sweet potato and banana displayed more complex compositional metabolite profiles with representations of up to 22 chemical classes (unknowns excluded) than that of potato, for which only metabolites from 10 chemical classes were detected. Additionally, over 20% of biochemical signatures remained unidentified for every crop analyzed. Integration of metabolomics with the on-going genomic and phenotypic studies will enhance 'omics-wide associations of molecular signatures with agronomic and consumer traits via easily quantifiable biochemical markers to aid gene discovery and functional characterization.

The evolution of metabolism: How to test evolutionary hypotheses at the genomic level

The origin of primordial metabolism and its expansion to form the metabolic networks extant today represent excellent systems to study the impact of natural selection and the potential adaptive role of novel compounds. Here we present the current hypotheses made on the origin of life and ancestral metabolism and present the theories and mechanisms by which the large chemical diversity of plants might have emerged along evolution. In particular, we provide a survey of statistical methods that can be used to detect signatures of selection at the gene and population level, and discuss potential and limits of these methods for investigating patterns of molecular adaptation in plant metabolism.

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.

Systems genetics of the Drosophila metabolome

How effects of DNA sequence variants are transmitted through intermediate endophenotypes to modulate organismal traits remains a central question in quantitative genetics. This problem can be addressed through a systems approach in a population in which genetic polymorphisms, gene expression traits, metabolites, and complex phenotypes can be evaluated on the same genotypes. Here, we focused on the metabolome, which represents the most proximal link between genetic variation and organismal phenotype, and quantified metabolite levels in 40 lines of the Drosophila melanogaster Genetic Reference Panel. We identified sex-specific modules of genetically correlated metabolites and constructed networks that integrate DNA sequence variation and variation in gene expression with variation in metabolites and organismal traits, including starvation stress resistance and male aggression. Finally, we asked to what extent SNPs and metabolites can predict trait phenotypes and generated trait- and sex-specific prediction models that provide novel insights about the metabolomic underpinnings of complex phenotypes.

The Structure and Function of Major Plant Metabolite Modifications

Plants produce a myriad of structurally and functionally diverse metabolites that play many different roles in plant growth and development and in plant response to continually changing environmental conditions as well as abiotic and biotic stresses. This metabolic diversity is, to a large extent, due to chemical modification of the basic skeletons of metabolites. Here, we review the major known plant metabolite modifications and summarize the progress that has been achieved and the challenges we are facing in the field. We focus on discussing both technical and functional aspects in studying the influences that various modifications have on biosynthesis, degradation, transport, and storage of metabolites, as well as their bioactivity and toxicity. Finally, we discuss some emerging insights into the evolution of metabolic pathways and metabolite functionality.

Integration of Metabolomic and Other Omics Data in Population-Based Study Designs: An Epidemiological Perspective

It is not controversial that study design considerations and challenges must be addressed when investigating the linkage between single omic measurements and human phenotypes. It follows that such considerations are just as critical, if not more so, in the context of multi-omic studies. In this review, we discuss (1) epidemiologic principles of study design, including selection of biospecimen source(s) and the implications of the timing of sample collection, in the context of a multi-omic investigation, and (2) the strengths and limitations of various techniques of data integration across multi-omic data types that may arise in population-based studies utilizing metabolomic data.

Metabolome-Scale Genome-Wide Association Studies Reveal Chemical Diversity and Genetic Control of Maize Specialized Metabolites

Cultivated maize (Zea mays) has retained much of the genetic diversity of its wild ancestors. Here, we performed nontargeted liquid chromatography-mass spectrometry metabolomics to analyze the metabolomes of the 282 maize inbred lines in the Goodman Diversity Panel. This analysis identified a bimodal distribution of foliar metabolites. Although 15% of the detected mass features were present in >90% of the inbred lines, the majority were found in <50% of the samples. Whereas leaf bases and tips were differentiated by flavonoid abundance, maize varieties (stiff-stalk, nonstiff-stalk, tropical, sweet maize, and popcorn) showed differential accumulation of benzoxazinoid metabolites. Genome-wide association studies (GWAS), performed for 3,991 mass features from the leaf tips and leaf bases, showed that 90% have multiple significantly associated loci scattered across the genome. Several quantitative trait locus hotspots in the maize genome regulate the abundance of multiple, often structurally related mass features. The utility of maize metabolite GWAS was demonstrated by confirming known benzoxazinoid biosynthesis genes, as well as by mapping isomeric variation in the accumulation of phenylpropanoid hydroxycitric acid esters to a single linkage block in a citrate synthase-like gene. Similar to gene expression databases, this metabolomic GWAS data set constitutes an important public resource for linking maize metabolites with biosynthetic and regulatory genes.

Nonlinear phenotypic variation uncovers the emergence of heterosis in Arabidopsis thaliana

Heterosis describes the phenotypic superiority of hybrids over their parents in traits related to agronomic performance and fitness. Understanding and predicting nonadditive inheritance such as heterosis is crucial for evolutionary biology as well as for plant and animal breeding. However, the physiological bases of heterosis remain debated. Moreover, empirical data in various species have shown that diverse genetic and molecular mechanisms are likely to explain heterosis, making it difficult to predict its emergence and amplitude from parental genotypes alone. In this study, we examined a model of physiological dominance initially proposed by Sewall Wright to explain the nonadditive inheritance of traits like metabolic fluxes at the cellular level. We evaluated Wright's model for two fitness-related traits at the whole-plant level, growth rate and fruit number, using 450 hybrids derived from crosses among natural accessions of A. thaliana. We found that allometric relationships between traits constrain phenotypic variation in a nonlinear and similar manner in hybrids and accessions. These allometric relationships behave predictably, explaining up to 75% of heterosis amplitude, while genetic distance among parents at best explains 7%. Thus, our findings are consistent with Wright's model of physiological dominance and suggest that the emergence of heterosis on plant performance is an intrinsic property of nonlinear relationships between traits. Furthermore, our study highlights the potential of a geometric approach of phenotypic relationships for predicting heterosis of major components of crop productivity and yield.

Identification and characterization of metabolite quantitative trait loci in tomato leaves and comparison with those reported for fruits and seeds

INTRODUCTION

To date, most studies of natural variation and metabolite quantitative trait loci (mQTL) in tomato have focused on fruit metabolism, leaving aside the identification of genomic regions involved in the regulation of leaf metabolism.

OBJECTIVE

This study was conducted to identify leaf mQTL in tomato and to assess the association of leaf metabolites and physiological traits with the metabolite levels from other tissues.

METHODS

The analysis of components of leaf metabolism was performed by phenotypying 76 tomato ILs with chromosome segments of the wild species Solanum pennellii in the genetic background of a cultivated tomato (S. lycopersicum) variety M82. The plants were cultivated in two different environments in independent years and samples were harvested from mature leaves of non-flowering plants at the middle of the light period. The non-targeted metabolite profiling was obtained by gas chromatography time-of-flight mass spectrometry (GC-TOF-MS). With the data set obtained in this study and already published metabolomics data from seed and fruit, we performed QTL mapping, heritability and correlation analyses.

RESULTS

Changes in metabolite contents were evident in the ILs that are potentially important with respect to stress responses and plant physiology. By analyzing the obtained data, we identified 42 positive and 76 negative mQTL involved in carbon and nitrogen metabolism.

CONCLUSIONS

Overall, these findings allowed the identification of S. lycopersicum genome regions involved in the regulation of leaf primary carbon and nitrogen metabolism, as well as the association of leaf metabolites with metabolites from seeds and fruits.

Metabolite profiling and genome-wide association studies reveal response mechanisms of phosphorus deficiency in maize seedling

Inorganic phosphorus (Pi) is an essential element in numerous metabolic reactions and signaling pathways, but the molecular details of these pathways remain largely unknown. In this study, metabolite profiles of maize (Zea mays L.) leaves and roots were compared between six low-Pi-sensitive lines and six low-Pi-tolerant lines under Pi-sufficient and Pi-deficient conditions to identify pathways and genes associated with the low-Pi stress response. Results showed that under Pi deprivation the concentrations of nucleic acids, organic acids and sugars were increased, but that the concentrations of phosphorylated metabolites, certain amino acids, lipid metabolites and nitrogenous compounds were decreased. The levels of secondary metabolites involved in plant immune reactions, including benzoxazinoids and flavonoids, were significantly different in plants grown under Pi-deficient conditions. Among them, the 11 most stable metabolites showed significant differences under low- and normal-Pi conditions based on the coefficient of variation (CV). Isoleucine and alanine were the most stable metabolites for the identification of Pi-sensitive and Pi-resistant maize inbred lines. With the significant correlation between morphological traits and metabolites, five low-Pi-responding consensus genes associated with morphological traits and simultaneously involved in metabolic pathways were mined by combining metabolites profiles and genome-wide association study (GWAS). The consensus genes induced by Pi deficiency in maize seedlings were also validated by reverse-transcription quantitative polymerase chain reaction (RT-qPCR). Moreover, these genes were further validated in a recombinant inbred line (RIL) population, in which the glucose-6-phosphate-1-epimerase encoding gene mediated yield and correlated traits to phosphorus availability. Together, our results provide a framework for understanding the metabolic processes underlying Pi-deficient responses and give multiple insights into improving the efficiency of Pi use in maize.

Genomic architecture of complex traits in loblolly pine

Dissecting the genetic and genomic architecture of complex traits is essential to understand the forces maintaining the variation in phenotypic traits of ecological and economical importance. Whole-genome resequencing data were used to generate high-resolution polymorphic single nucleotide polymorphism (SNP) markers and genotype individuals from common gardens across the loblolly pine (Pinus taeda) natural range. Genome-wide associations were tested with a large phenotypic dataset comprising 409 variables including morphological traits (height, diameter, carbon isotope discrimination, pitch canker resistance), and molecular traits such as metabolites and expression of xylem development genes. Our study identified 2335 new SNP × trait associations for the species, with many SNPs located in physical clusters in the genome of the species; and the genomic location of hotspots for metabolic × genotype associations. We found a highly polygenic basis of quantitative inheritance, with significant differences in number, effects size, genomic location and frequency of alleles contributing to variation in phenotypes in the different traits. While mutation-selection balance might be shaping the genetic variation in metabolic traits, balancing selection is more likely to shape the variation in expression of xylem development genes. Our work contributes to the study of complex traits in nonmodel plant species by identifying associations at a whole-genome level.

Metabolic GWAS-based dissection of genetic bases underlying the diversity of plant metabolism

Plants have served as sources providing humans with metabolites for food and nutrition, biomaterials for living, and treatment for pain and disease. Plants produce a huge array of metabolites, with an immense diversity at both the population and individual levels. Dissection of the genetic bases for metabolic diversity has attracted increasing research attention. The concept of genome-wide association study (GWAS) was extended to studies on the diversity of plant metabolome that benefitted from the development of mass-spectrometry-based analytical systems and genome sequencing technologies. Metabolic genome-wide association study (mGWAS) is one of the most powerful tools for global identification of genetic determinants for diversity of plant metabolism. Recently, mGWAS has been performed for various species with continuous improvements, providing deeper insights into the genetic bases of metabolic diversity. In this review, we discuss fully the achievements to date and remaining challenges that are associated with both mGWAS and mGWAS-based multi-dimensional analysis. We begin with a summary of GWAS and its development based on statistical methods and populations. As variation in targeted traits is essential for GWAS, we review metabolic diversity and its rise at both the population and individual levels. Subsequently, the application of mGWAS for plants and its corresponding achievements are fully discussed. We address the current knowledge on mGWAS-based multi-dimensional analysis and emerging insights into the diversity of metabolism.