Simplivariate Models: Ideas and First Examples

One of the new expanding areas in functional genomics is metabolomics: measuring the metabolome of an organism. Data being generated in metabolomics studies are very diverse in nature depending on the design underlying the experiment. Traditionally, variation in measurements is conceptually broken down in systematic variation and noise where the latter contains, e.g. technical variation. There is increasing evidence that this distinction does not hold (or is too simple) for metabolomics data. A more useful distinction is in terms of informative and non-informative variation where informative relates to the problem being studied. In most common methods for analyzing metabolomics (or any other high-dimensional x-omics) data this distinction is ignored thereby severely hampering the results of the analysis. This leads to poorly interpretable models and may even obscure the relevant biological information. We developed a framework from first data analysis principles by explicitly formulating the problem of analyzing metabolomics data in terms of informative and non-informative parts. This framework allows for flexible interactions with the biologists involved in formulating prior knowledge of underlying structures. The basic idea is that the informative parts of the complex metabolomics data are approximated by simple components with a biological meaning, e.g. in terms of metabolic pathways or their regulation. Hence, we termed the framework ‘simplivariate models’ which constitutes a new way of looking at metabolomics data. The framework is given in its full generality and exemplified with two methods, IDR analysis and plaid modeling, that fit into the framework. Using this strategy of ‘divide and conquer’, we show that meaningful simplivariate models can be obtained using a real-life microbial metabolomics data set. For instance, one of the simple components contained all the measured intermediates of the Krebs cycle of E. coli. Moreover, these simplivariate models were able to uncover regulatory mechanisms present in the phenylalanine biosynthesis route of E. coli

Plasma Metabolomic Profiles Reflective of Glucose Homeostasis in Non-Diabetic and Type 2 Diabetic Obese African-American Women

Insulin resistance progressing to type 2 diabetes mellitus (T2DM) is marked by a broad perturbation of macronutrient intermediary metabolism. Understanding the biochemical networks that underlie metabolic homeostasis and how they associate with insulin action will help unravel diabetes etiology and should foster discovery of new biomarkers of disease risk and severity. We examined differences in plasma concentrations of >350 metabolites in fasted obese T2DM vs. obese non-diabetic African-American women, and utilized principal components analysis to identify 158 metabolite components that strongly correlated with fasting HbA1c over a broad range of the latter (r = −0.631; p<0.0001). In addition to many unidentified small molecules, specific metabolites that were increased significantly in T2DM subjects included certain amino acids and their derivatives (i.e., leucine, 2-ketoisocaproate, valine, cystine, histidine), 2-hydroxybutanoate, long-chain fatty acids, and carbohydrate derivatives. Leucine and valine concentrations rose with increasing HbA1c, and significantly correlated with plasma acetylcarnitine concentrations. It is hypothesized that this reflects a close link between abnormalities in glucose homeostasis, amino acid catabolism, and efficiency of fuel combustion in the tricarboxylic acid (TCA) cycle. It is speculated that a mechanism for potential TCA cycle inefficiency concurrent with insulin resistance is “anaplerotic stress” emanating from reduced amino acid-derived carbon flux to TCA cycle intermediates, which if coupled to perturbation in cataplerosis would lead to net reduction in TCA cycle capacity relative to fuel delivery

Evaluation of metabolomics profiles of grain from maize hybrids derived from near-isogenic GM positive and negative segregant inbreds demonstrates that observed differences cannot be attributed unequivocally to the GM trait

Introduction: Past studies on plant metabolomes have highlighted the influence of growing environments and varietal differences in variation of levels of metabolites yet there remains continued interest in evaluating the effect of genetic modification (GM). Objectives: Here we test the hypothesis that metabolomics differences in grain from maize hybrids derived from a series of GM (NK603, herbicide tolerance) inbreds and corresponding negative segregants can arise from residual genetic variation associated with backcrossing and that the effect of insertion of the GM trait is negligible. Methods: Four NK603-positive and negative segregant inbred males were crossed with two different females (testers). The resultant hybrids, as well as conventional comparator hybrids, were then grown at three replicated field sites in Illinois, Minnesota, and Nebraska during the 2013 season. Metabolomics data acquisition using gas chromatography–time of flight-mass spectrometry (GC–TOF-MS) allowed the measurement of 367 unique metabolite features in harvested grain, of which 153 were identified with small molecule standards. Multivariate analyses of these data included multi-block principal component analysis and ANOVA-simultaneous component analysis. Univariate analyses of all 153 identified metabolites was conducted based on significance testing (α = 0.05), effect size evaluation (assessing magnitudes of differences), and variance component analysis. Results: Results demonstrated that the largest effects on metabolomic variation were associated with different growing locations and the female tester. They further demonstrated that differences observed between GM and non-GM comparators, even in stringent tests utilizing near-isogenic positive and negative segregants, can simply reflect minor genomic differences associated with conventional back-crossing practices. Conclusion: The effect of GM on metabolomics variation was determined to be negligible and supports that there is no scientific rationale for prioritizing GM as a source of variation.</p

Phase 1b randomized, double-blind study of namilumab, an anti-granulocyte macrophage colony-stimulating factor monoclonal antibody, in mild-to-moderate rheumatoid arthritis

Change from baseline in swollen (a) and tender (b) joint counts. *Error bars show upper SE for placebo and lower SE for namilumab. SE standard error, SJC swollen joint count, TJC tender joint count. (PDF 1292 kb

Between Metabolite Relationships: an essential aspect of metabolic change

Not only the levels of individual metabolites, but also the relations between the levels of different metabolites may indicate (experimentally induced) changes in a biological system. Component analysis methods in current ‘standard’ use for metabolomics, such as Principal Component Analysis (PCA), do not focus on changes in these relations. We therefore propose the concept of ‘Between Metabolite Relationships’ (BMRs): common changes in the covariance (or correlation) between all metabolites in an organism. Such structural changes may indicate metabolic change brought about by experimental manipulation but which are lost with standard data analysis methods. These BMRs can be analysed by the INdividual Differences SCALing (INDSCAL) method. First the BMR quantification is described and subsequently the INDSCAL method. Finally, two studies illustrate the power and the applicability of BMRs in metabolomics. The first study is about the induced plant response of cabbage to herbivory, of which BMRs are a considerable part. In the second study—a human nutritional intervention study of green tea extract—standard data analysis tools did not reveal any metabolic change, although the BMRs were considerably affected. The presented results show that BMRs can be easily implemented in a wide variety of metabolomic studies. They provide a new source of information to describe biological systems in a way that fits flawlessly into the next generation of systems biology questions, dealing with personalized responses

Development of a kinetic metabolic model: application to Catharanthus roseus hairy root

A kinetic metabolic model describing Catharanthus roseus hairy root growth and nutrition was developed. The metabolic network includes glycolysis, pentose-phosphate pathway, TCA cycle and the catabolic reactions leading to cell building blocks such as amino acids, organic acids, organic phosphates, lipids and structural hexoses. The central primary metabolic network was taken at pseudo-steady state and metabolic flux analysis technique allowed reducing from 31 metabolic fluxes to 20 independent pathways. Hairy root specific growth rate was described as a function of intracellular concentration in cell building blocks. Intracellular transport and accumulation kinetics for major nutrients were included. The model uses intracellular nutrients as well as energy shuttles to describe metabolic regulation. Model calibration was performed using experimental data obtained from batch and medium exchange liquid cultures of C. roseus hairy root using a minimal medium in Petri dish. The model is efficient in estimating the growth rate

Covering Chemical Diversity of Genetically-Modified Tomatoes Using Metabolomics for Objective Substantial Equivalence Assessment

As metabolomics can provide a biochemical snapshot of an organism's phenotype it is a promising approach for charting the unintended effects of genetic modification. A critical obstacle for this application is the inherently limited metabolomic coverage of any single analytical platform. We propose using multiple analytical platforms for the direct acquisition of an interpretable data set of estimable chemical diversity. As an example, we report an application of our multi-platform approach that assesses the substantial equivalence of tomatoes over-expressing the taste-modifying protein miraculin. In combination, the chosen platforms detected compounds that represent 86% of the estimated chemical diversity of the metabolites listed in the LycoCyc database. Following a proof-of-safety approach, we show that % had an acceptable range of variation while simultaneously indicating a reproducible transformation-related metabolic signature. We conclude that multi-platform metabolomics is an approach that is both sensitive and robust and that it constitutes a good starting point for characterizing genetically modified organisms

Oral microbe-host interactions: influence of β-glucans on gene expression of inflammatory cytokines and metabolome profile

Background: The aim of this study was to evaluate the effects of β-glucan on the expression of inflammatory mediators and metabolomic profile of oral cells [keratinocytes (OBA-9) and fibroblasts (HGF-1) in a dual-chamber model] infected by Aggregatibacter actinomycetemcomitans. The periodontopathogen was applied and allowed to cross the top layer of cells (OBA-9) to reach the bottom layer of cells (HGF-1) and induce the synthesis of immune factors and cytokines in the host cells. β-glucan (10 μg/mL or 20 μg/mL) were added, and the transcriptional factors and metabolites produced were quantified in the remaining cell layers and supernatant. Results: The relative expression of interleukin (IL)-1-α and IL-18 genes in HGF-1 decreased with 10 μg/mL or 20 μg/mL of β-glucan, where as the expression of PTGS-2 decreased only with 10 μg/mL. The expression of IL-1-α increased with 20 μg/mL and that of IL-18 increased with 10 μg/mL in OBA-9; the expression of BCL 2, EP 300, and PTGS-2 decreased with the higher dose of β-glucan. The production of the metabolite 4-aminobutyric acid presented lower concentrations under 20 μg/mL, whereas the concentrations of 2-deoxytetronic acid NIST and oxalic acid decreased at both concentrations used. Acetophenone, benzoic acid, and pinitol presented reduced concentrations only when treated with 10 μg/mL of β-glucan. Conclusions: Treatment with β-glucans positively modulated the immune response and production of metabolites