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.     

Metabolic profiling of strawberry (Fragaria x ananassa Duch.) during fruit development and maturation

Strawberry (Fragaria × ananassa Duch), a fruit of economic and nutritional importance, is also a model species for fleshy fruits and genomics in Rosaceae. Strawberry fruit quality at different harvest stages is a function of the fruit's metabolite content, which results from physiological changes during fruit growth and ripening. In order to investigate strawberry fruit development, untargeted (GC-MS) and targeted (HPLC) metabolic profiling analyses were conducted. Principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA) were employed to explore the non-polar and polar metabolite profiles from fruit samples at seven developmental stages. Different cluster patterns and a broad range of metabolites that exerted influence on cluster formation of metabolite profiles were observed. Significant changes in metabolite levels were found in both fruits turning red and fruits over-ripening in comparison with red-ripening fruits. The levels of free amino acids decreased gradually before the red-ripening stage, but increased significantly in the over-ripening stage. Metabolite correlation and network analysis revealed the interdependencies of individual metabolites and metabolic pathways. Activities of several metabolic pathways, including ester biosynthesis, the tricarboxylic acid cycle, the shikimate pathway, and amino acid metabolism, shifted during fruit growth and ripening. These results not only confirmed published metabolic data but also revealed new insights into strawberry fruit composition and metabolite changes, thus demonstrating the value of metabolomics as a functional genomics tool in characterizing the mechanism of fruit quality formation, a key developmental stage in most economically important fruit crops.     

The potential of metabolite profiling as a selection tool for genotype discrimination in Populus

Differences between wild-type Populus tremulaxalba and two transgenic lines with modified lignin monomer composition, were interrogated using metabolic profiling. Analysis of metabolite abundance data by GC-MS, coupled with principal components analysis (PCA), successfully differentiated between lines that had distinct phenotypes, whether samples were taken from the cambial zone or non-lignifying suspension tissue cultures. Interestingly, the GC-MS analysis detected relatively few phenolic metabolites in cambial extracts, although a single metabolite associated with the differentiation between lines was directly related to the phenylpropanoid pathway or other down-stream aspects of lignin biosynthesis. In fact, carbohydrates, which have only an indirect relationship with the modified lignin monomer composition, featured strongly in the line-differentiating aspects of the statistical analysis. Traditional HPLC analysis was employed to verify the GC-MS data. These findings demonstrate that metabolic traits can be dissected reliably and accurately by metabolomic analyses, enabling the discrimination of individual genotypes of the same tree species that exhibit marked differences in industrially relevant wood traits. Furthermore, this validates the potential of using metabolite profiling techniques for marker generation in the context of plant/tree breeding for industrial applications.     

1H NMR metabolomics identification of markers of hypoxia-induced metabolic shifts in a breast cancer model system

Hypoxia can promote invasive behavior in cancer cells and alters the response to therapeutic intervention as a result of changes in the expression many genes, including genes involved in intermediary metabolism. Although metabolomics technologies are capable of simultaneously measuring a wide range of metabolites in an untargeted manner, these methods have been relatively under utilized in the study of cancer cell responses to hypoxia. Thus, (1)H NMR metabolomics was used to examine the effects of hypoxia in the MDA-MB-231 human breast cancer cell line, both in vitro and in vivo. Cell cultures were compared with respect to their metabolic responses during growth under either hypoxic (1% O(2)) or normoxic conditions. Orthogonal partial least squares discriminant analysis (OPLS-DA) was used to identify a set of metabolites that were responsive to hypoxia. Via intracardiac administration, MDA-MB-231 cells were also used to generate widespread metastatic disease in immuno-compromised mice. Serum metabolite analysis was conducted to compare animals with and without a large tumor burden. Intriguingly, using a cross-plot of the OPLS loadings, both the in vitro and in vivo samples yielded a subset of metabolites that were significantly altered by hypoxia. These included primarily energy metabolites and amino acids, indicative of known alterations in energy metabolism, and possibly protein synthesis or catabolism. The results suggest that the metabolite pattern identified might prove useful as a marker for intra-tumoral hypoxia.     

Organic acid profile of isovaleric acidemia: a comprehensive metabolomics approach

Isovaleric acidemia (IVA, MIM 248600) can be a severe and potentially life-threatening disease in affected neonates, but with a positive prognosis on treatment for some phenotypes. This study presents the first application of metabolomics to evaluate the metabolite profiles derived from urine samples of untreated and treated IVA patients as well as of obligate heterozygotes. All IVA patients carried the same homozygous c.367 G > A nucleotide change in exon 4 of the IVD gene but manifested phenotypic diversity. Concurrent class analysis (CONCA) was used to compare all the metabolites from the original complete data set obtained from the three case and two control groups used in this investigation. This application of CONCA has not been reported previously, and is used here to compare four different modes of scaling of all metabolites. The variables important in discrimination from the CONCA thus enabled the recognition of different metabolic patterns encapsulated within the data sets that would not have been revealed by using only one mode of scaling. Application of multivariate and univariate analyses disclosed 11 important metabolites that distinguished untreated IVA from controls. These included well-established diagnostic biomarkers of IVA, endogenous detoxification markers, and 3-hydroxycaproic acid, an indicator of ketosis, but not reported previously for this disease. Nine metabolites were identified that reflected the effect of treatment of IVA. They included detoxification products and indicators related to the high carbohydrate and low protein diet which formed the hallmark of the treatment. This investigation also provides the first comparative metabolite profile for heterozygotes of this inherited metabolic disorder. The detection of informative metabolites in even very low concentrations in all three experimental groups highlights the potential advantage of the holistic mode of analysis of inherited metabolic diseases in a metabolomics investigation.     

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.     

A functional genomics strategy that uses metabolome data to reveal the phenotype of silent mutations

A large proportion of the 6,000 genes present in the genome of Saccharomyces cerevisiae, and of those sequenced in other organisms, encode proteins of unknown function. Many of these genes are "silent, " that is, they show no overt phenotype, in terms of growth rate or other fluxes, when they are deleted from the genome. We demonstrate how the intracellular concentrations of metabolites can reveal phenotypes for proteins active in metabolic regulation. Quantification of the change of several metabolite concentrations relative to the concentration change of one selected metabolite can reveal the site of action, in the metabolic network, of a silent gene. In the same way, comprehensive analyses of metabolite concentrations in mutants, providing "metabolic snapshots," can reveal functions when snapshots from strains deleted for unstudied genes are compared to those deleted for known genes. This approach to functional analysis, using comparative metabolomics, we call FANCY-an abbreviation for functional analysis by co-responses in yeast.     

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.     

Comparative metabolomics charts the impact of genotype-dependent methionine accumulation in Arabidopsis thaliana

Methionine (Met) is an essential amino acid for all organisms. In plants, Met also functions as a precursor of plant hormones, polyamines, and defense metabolites. The regulatory mechanism of Met biosynthesis is highly complex and, despite its great importance, remains unclear. To investigate how accumulation of Met influences metabolism as a whole in Arabidopsis, three methionine over-accumulation (mto) mutants were examined using a gas chromatography–mass spectrometry-based metabolomics approach. Multivariate statistical analyses of the three mto mutants (mto1, mto2, and mto3) revealed distinct metabolomic phenotypes. Orthogonal projection to latent structures–discriminant analysis highlighted discriminative metabolites contributing to the separation of each mutant and the corresponding control samples. Though Met accumulation in mto1 had no dramatic effect on other metabolic pathways except for the aspartate family, metabolite profiles of mto2 and mto3 indicated that several extensive pathways were affected in addition to over-accumulation of Met. The pronounced changes in metabolic pathways in both mto2 and mto3 were associated with polyamines. The findings suggest that our metabolomics approach not only can reveal the impact of Met over-accumulation on metabolism, but also may provide clues to identify crucial pathways for regulation of metabolism in plants.     

Water-use efficiency and whole-plant performance of nine tropical tree species at two sites with contrasting water availability in Panama

Across their natural distributions, tropical tree species are regularly exposed to seasonal droughts of varying intensities. Their ability to tolerate drought stress plays a vital role in determining growth and mortality rates, as well as shaping the functional composition of tropical forests. In order to assess the ability of species to acclimate to contrasting levels of drought stress, physiological and structural traits involved in drought adaptation—wood C isotope discrimination (δ¹³C), wood specific gravity, and wood C content—of 2-year-old saplings of nine tropical tree species were evaluated in common garden experiments at two study sites in Panama with contrasting seasonality. We assessed co-variation in wood traits with relative growth rates (RGR_BD), aboveground biomass, and basal diameter and the plasticity of wood traits across study sites. Overall, species responded to lower water availability by increasing intrinsic water-use efficiency, i.e., less negative wood δ¹³C, but did not exhibit a uniform, directional response for wood specific gravity or wood C content. Trait plasticity for all wood traits was independent of RGR_BD and tree size. We found that the adaptive value of intrinsic water-use efficiency varied with water availability. Intrinsic water-use efficiency increased with decreasing RGR_BD at the more seasonal site, facilitating higher survival of slower growing species. Conversely, intrinsic water-use efficiency increased with tree size at the less seasonal site, which conferred a competitive advantage to larger individuals at the cost of greater susceptibility to drought-induced mortality. Our results illustrate that acclimation to water availability has negligible impacts on tree growth over short periods, but eventually could favor slow-growing species with conservative water-use strategies in tropical regions experiencing increasingly frequent and severe droughts.     

Mining the Unknown: A Systems Approach to Metabolite Identification Combining Genetic and Metabolic Information

Recent genome-wide association studies (GWAS) with metabolomics data linked genetic variation in the human genome to differences in individual metabolite levels. A strong relevance of this metabolic individuality for biomedical and pharmaceutical research has been reported. However, a considerable amount of the molecules currently quantified by modern metabolomics techniques are chemically unidentified. The identification of these “unknown metabolites” is still a demanding and intricate task, limiting their usability as functional markers of metabolic processes. As a consequence, previous GWAS largely ignored unknown metabolites as metabolic traits for the analysis. Here we present a systems-level approach that combines genome-wide association analysis and Gaussian graphical modeling with metabolomics to predict the identity of the unknown metabolites. We apply our method to original data of 517 metabolic traits, of which 225 are unknowns, and genotyping information on 655,658 genetic variants, measured in 1,768 human blood samples. We report previously undescribed genotype–metabotype associations for six distinct gene loci (SLC22A2, COMT, CYP3A5, CYP2C18, GBA3, UGT3A1) and one locus not related to any known gene (rs12413935). Overlaying the inferred genetic associations, metabolic networks, and knowledge-based pathway information, we derive testable hypotheses on the biochemical identities of 106 unknown metabolites. As a proof of principle, we experimentally confirm nine concrete predictions. We demonstrate the benefit of our method for the functional interpretation of previous metabolomics biomarker studies on liver detoxification, hypertension, and insulin resistance. Our approach is generic in nature and can be directly transferred to metabolomics data from different experimental platforms.     

Mass spectrometry methods in metabolomics

The review deals with metabolomics, a new and rapidly growing area directed to the comprehensive analysis of metabolites of biological objects. Metabolites are characterized by various physical and chemical properties, traditionally studied by methods of analytical chemistry focused on certain groups of chemical substances. However, current progress in mass spectrometry has led to formation of rather unified methods, such as metabolic fingerprinting and metabolomic profiling, which allow defining thousands of metabolites in one biological sample and therefore draw “a modern portrait of metabolomics.” This review describes basic characteristics of these methods, ways of metabolite separation, and analysis of metabolites by mass spectrometry. The examples shown in this review, allow to estimate these methods and to compare their advantages and disadvantages. Besides that, we consider the methods, which are of the most frequent use in metabolomics; these include the methods for data processing and the required resources, such as software for mass spectra processing and metabolite search database. In the conclusion, general suggestions for successful metabolomic experiments are given.     

Genetic Influences on Metabolite Levels: A Comparison across Metabolomic Platforms

Metabolomic profiling is a powerful approach to characterize human metabolism and help understand common disease risk. Although multiple high-throughput technologies have been developed to assay the human metabolome, no technique is capable of capturing the entire human metabolism. Large-scale metabolomics data are being generated in multiple cohorts, but the datasets are typically profiled using different metabolomics platforms. Here, we compared analyses across two of the most frequently used metabolomic platforms, Biocrates and Metabolon, with the aim of assessing how complimentary metabolite profiles are across platforms. We profiled serum samples from 1,001 twins using both targeted (Biocrates, n = 160 metabolites) and non-targeted (Metabolon, n = 488 metabolites) mass spectrometry platforms. We compared metabolite distributions and performed genome-wide association analyses to identify shared genetic influences on metabolites across platforms. Comparison of 43 metabolites named for the same compound on both platforms indicated strong positive correlations, with few exceptions. Genome-wide association scans with high-throughput metabolic profiles were performed for each dataset and identified genetic variants at 7 loci associated with 16 unique metabolites on both platforms. The 16 metabolites showed consistent genetic associations and appear to be robustly measured across platforms. These included both metabolites named for the same compound across platforms as well as unique metabolites, of which 2 (nonanoylcarnitine (C9) [Biocrates]/Unknown metabolite X-13431 [Metabolon] and PC aa C28:1 [Biocrates]/1-stearoylglycerol [Metabolon]) are likely to represent the same or related biochemical entities. The results demonstrate the complementary nature of both platforms, and can be informative for future studies of comparative and integrative metabolomics analyses in samples profiled on different platforms.     

Consistency analysis of metabolic correlation networks

Background  

Metabolic correlation networks are derived from the covariance of metabolites in replicates of metabolomics experiments. They constitute an interesting intermediate between topology (i.e. the system's architecture defined by the set of reactions between metabolites) and dynamics (i.e. the metabolic concentrations observed as fluctuations around steady-state values in the metabolic network).     

Identifying metabolic alterations associated with coral growth anomalies using 1H NMR metabolomics

Coral growth anomalies (GAs) are tumor-like protrusions that are detrimental to coral health, affecting both the coral skeleton and soft tissues. These lesions are increasingly found throughout the tropics and are commonly associated with high human population density, yet little is known about the molecular pathology of the disease. Here, we investigate the metabolic impacts of GAs through ¹H nuclear magnetic resonance (NMR) metabolomics in Porites compressa tissues from a site of high disease prevalence (Coconut Island, Hawaii). We putatively identified 18 metabolites (8.1% of total annotated features) through complementary ¹H and ¹H–¹³C heteronuclear single quantum correlation NMR data that increase confidence in pathway analyses and may bolster future coral metabolite annotation efforts. Extract yield was elevated in both GA and unaffected (normal tissue from a diseased colony) compared to reference (normal tissue from GA-free colony) samples, potentially indicating elevated metabolic activity in GA-impacted colonies. Relatively high variation in metabolomic profiles among coral samples of the same treatment (i.e., inter-colony variation) confounded data interpretation, however, analyses of paired GA and unaffected samples identified 73 features that differed between these respective metabolome types. These features were largely annotated as unknowns, but 1-methylnicotinamide and trigonelline were found to be elevated in GA samples, while betaine, glycine, and histamine were lower in GA samples. Pathway analyses indicate decreased choline oxidation in GA samples, making this a pathway of interest for future targeted studies. Collectively, our results provide unique insights into GA pathophysiology by showing these lesions alter both the absolute and relative metabolism of affected colonies and by identifying features (metabolites and unknowns) and metabolic pathways of interest in GA pathophysiology going forward.

     

Comprehensive analysis of the metabolome of Pseudomonas putida S12 grown on different carbon sources

Metabolomics is an emerging, powerful, functional genomics technology that involves the comparative non-targeted analysis of the complete set of metabolites in an organism. We have set-up a robust quantitative metabolomics platform that allows the analysis of 'snapshot' metabolomes. In this study, we have applied this platform for the comprehensive analysis of the metabolite composition of Pseudomonas putida S12 grown on four different carbon sources, i.e. fructose, glucose, gluconate and succinate. This paper focuses on the microbial aspects of analyzing comprehensive metabolomes, and demonstrates that metabolomes can be analyzed reliably. The technical (i.e. sample work-up and analytical) reproducibility was on average 10%, while the biological reproducibility was approximately 40%. Moreover, the energy charge values of the microbial samples generated were determined, and indicated that no biotic or abiotic changes had occurred during sample work-up and analysis. In general, the metabolites present and their concentrations were very similar after growth on the different carbon sources. However, specific metabolites showed large differences in concentration, especially the intermediates involved in the degradation of the carbon sources studied. Principal component discriminant analysis was applied to identify metabolites that are specific for, i.e. not necessarily the metabolites that show those largest differences in concentration, cells grown on either of these four carbon sources. For selected enzymatic reactions, i.e. the glucose-6-phosphate isomerase, triosephosphate isomerase and phosphoglyceromutase reactions, the apparent equilibrium constants (K(app)) were calculated. In several instances a carbon source-dependent deviation between the apparent equilibrium constant (K(app)) and the thermodynamic equilibrium constant (K(eq)) was observed, hinting towards a potential point of metabolic regulation or towards bottlenecks in biosynthesis routes. For glucose-6-phosphate isomerase and phosphoglyceromutase, the K(app) was larger than K(eq), and the results suggested that the specific enzymatic activities of these two enzymes were too low to reach the thermodynamic equilibrium in growing cells. In contrast, with triosephosphate isomerase the K(app) was smaller than K(eq), and the results suggested that this enzyme is kinetically controlled.     

Metabolite profiling of CHO cells with different growth characteristics

Mammalian cell cultures are the predominant system for the production of recombinant proteins requiring post-translational modifications. As protein yields are a function of growth performance (among others), and performance varies greatly between culture medium (e.g., different growth rates and peak cell densities), an understanding of the biological mechanisms underpinning this variability would facilitate rational medium and process optimization, increasing product yields, and reducing costs. We employed a metabolomics approach to analyze differences in metabolite concentrations of CHO cells cultivated in three different media exhibiting different growth rates and maximum viable cell densities. Analysis of intra- and extracellular metabolite concentrations over the course of the cultures using a combination of HPLC and GC-MS, readily detected medium specific and time dependent changes. Using multivariate data analysis, we identified a range of metabolites correlating with growth rate, illustrating how metabolomics can be used to relate gross phenotypic changes to the fine details of cellular metabolism.