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

Search for new phenomena in final states with an energetic jet and large missing transverse momentum in pp collisions at √ s = 8 TeV with the ATLAS detector

Results of a search for new phenomena in final states with an energetic jet and large missing transverse momentum are reported. The search uses 20.3 fb−1 of √ s = 8 TeV data collected in 2012 with the ATLAS detector at the LHC. Events are required to have at least one jet with pT > 120 GeV and no leptons. Nine signal regions are considered with increasing missing transverse momentum requirements between Emiss T > 150 GeV and Emiss T > 700 GeV. Good agreement is observed between the number of events in data and Standard Model expectations. The results are translated into exclusion limits on models with either large extra spatial dimensions, pair production of weakly interacting dark matter candidates, or production of very light gravitinos in a gauge-mediated supersymmetric model. In addition, limits on the production of an invisibly decaying Higgs-like boson leading to similar topologies in the final state are presente

Genetics Meets Metabolomics: A Genome-Wide Association Study of Metabolite Profiles in Human Serum

The rapidly evolving field of metabolomics aims at a comprehensive measurement of ideally all endogenous metabolites in a cell or body fluid. It thereby provides a functional readout of the physiological state of the human body. Genetic variants that associate with changes in the homeostasis of key lipids, carbohydrates, or amino acids are not only expected to display much larger effect sizes due to their direct involvement in metabolite conversion modification, but should also provide access to the biochemical context of such variations, in particular when enzyme coding genes are concerned. To test this hypothesis, we conducted what is, to the best of our knowledge, the first GWA study with metabolomics based on the quantitative measurement of 363 metabolites in serum of 284 male participants of the KORA study. We found associations of frequent single nucleotide polymorphisms (SNPs) with considerable differences in the metabolic homeostasis of the human body, explaining up to 12% of the observed variance. Using ratios of certain metabolite concentrations as a proxy for enzymatic activity, up to 28% of the variance can be explained (p-values 10−16 to 10−21). We identified four genetic variants in genes coding for enzymes (FADS1, LIPC, SCAD, MCAD) where the corresponding metabolic phenotype (metabotype) clearly matches the biochemical pathways in which these enzymes are active. Our results suggest that common genetic polymorphisms induce major differentiations in the metabolic make-up of the human population. This may lead to a novel approach to personalized health care based on a combination of genotyping and metabolic characterization. These genetically determined metabotypes may subscribe the risk for a certain medical phenotype, the response to a given drug treatment, or the reaction to a nutritional intervention or environmental challenge

Metabolomics in Early Alzheimer's Disease: Identification of Altered Plasma Sphingolipidome Using Shotgun Lipidomics

The development of plasma biomarkers could facilitate early detection, risk assessment and therapeutic monitoring in Alzheimer's disease (AD). Alterations in ceramides and sphingomyelins have been postulated to play a role in amyloidogensis and inflammatory stress related neuronal apoptosis; however few studies have conducted a comprehensive analysis of the sphingolipidome in AD plasma using analytical platforms with accuracy, sensitivity and reproducibility.We prospectively analyzed plasma from 26 AD patients (mean MMSE 21) and 26 cognitively normal controls in a non-targeted approach using multi-dimensional mass spectrometry-based shotgun lipidomics to determine the levels of over 800 molecular species of lipids. These data were then correlated with diagnosis, apolipoprotein E4 genotype and cognitive performance. Plasma levels of species of sphingolipids were significantly altered in AD. Of the 33 sphingomyelin species tested, 8 molecular species, particularly those containing long aliphatic chains such as 22 and 24 carbon atoms, were significantly lower (p<0.05) in AD compared to controls. Levels of 2 ceramide species (N16:0 and N21:0) were significantly higher in AD (p<0.05) with a similar, but weaker, trend for 5 other species. Ratios of ceramide to sphingomyelin species containing identical fatty acyl chains differed significantly between AD patients and controls. MMSE scores were correlated with altered mass levels of both N20:2 SM and OH-N25:0 ceramides (p<0.004) though lipid abnormalities were observed in mild and moderate AD. Within AD subjects, there were also genotype specific differences.In this prospective study, we used a sensitive multimodality platform to identify and characterize an essentially uniform but opposite pattern of disruption in sphingomyelin and ceramide mass levels in AD plasma. Given the role of brain sphingolipids in neuronal function, our findings provide new insights into the AD sphingolipidome and the potential use of metabolomic signatures as peripheral biomarkers

An Effective Assessment of Simvastatin-Induced Toxicity with NMR-Based Metabonomics Approach

BACKGROUND: Simvastatin, which is used to control elevated cholesterol levels, is one of the most widely prescribed drugs. However, a daily excessive dose can induce drug-toxicity, especially in muscle and liver. Current markers for toxicity reflect mostly the late stages of tissue damage; thus, more efficient methods of toxicity evaluation are desired. METHODOLOGY/PRINCIPAL FINDINGS: As a new way to evaluate toxicity, we performed NMR-based metabonomics analysis of urine samples. Compared to conventional markers, such as AST, ALT, and CK, the urine metabolic profile provided clearer distinction between the pre- and post-treatment groups treated with toxic levels of simvastatin. Through multivariate statistical analysis, we identified marker metabolites associated with the toxicity. Importantly, we observed that the treatment group could be further categorized into two subgroups based on the NMR profiles: weak toxicity (WT) and high toxicity (HT). The distinction between these two groups was confirmed by the enzyme values and histopathological exams. Time-dependent studies showed that the toxicity at 10 days could be reliably predicted from the metabolic profiles at 6 days. CONCLUSIONS/SIGNIFICANCE: This metabonomics approach may provide a non-invasive and effective way to evaluate the simvastatin-induced toxicity in a manner that can complement current measures. The approach is expected to find broader application in other drug-induced toxicity assessments

Enteric Microbiome Metabolites Correlate with Response to Simvastatin Treatment

Although statins are widely prescribed medications, there remains considerable variability in therapeutic response. Genetics can explain only part of this variability. Metabolomics is a global biochemical approach that provides powerful tools for mapping pathways implicated in disease and in response to treatment. Metabolomics captures net interactions between genome, microbiome and the environment. In this study, we used a targeted GC-MS metabolomics platform to measure a panel of metabolites within cholesterol synthesis, dietary sterol absorption, and bile acid formation to determine metabolite signatures that may predict variation in statin LDL-C lowering efficacy. Measurements were performed in two subsets of the total study population in the Cholesterol and Pharmacogenetics (CAP) study: Full Range of Response (FR), and Good and Poor Responders (GPR) were 100 individuals randomly selected from across the entire range of LDL-C responses in CAP. GPR were 48 individuals, 24 each from the top and bottom 10% of the LDL-C response distribution matched for body mass index, race, and gender. We identified three secondary, bacterial-derived bile acids that contribute to predicting the magnitude of statin-induced LDL-C lowering in good responders. Bile acids and statins share transporters in the liver and intestine; we observed that increased plasma concentration of simvastatin positively correlates with higher levels of several secondary bile acids. Genetic analysis of these subjects identified associations between levels of seven bile acids and a single nucleotide polymorphism (SNP), rs4149056, in the gene encoding the organic anion transporter SLCO1B1. These findings, along with recently published results that the gut microbiome plays an important role in cardiovascular disease, indicate that interactions between genome, gut microbiome and environmental influences should be considered in the study and management of cardiovascular disease. Metabolic profiles could provide valuable information about treatment outcomes and could contribute to a more personalized approach to therapy

A Novel Conserved Isoform of the Ubiquitin Ligase UFD2a/UBE4B Is Expressed Exclusively in Mature Striated Muscle Cells

Yeast Ufd2p was the first identified E4 multiubiquitin chain assembly factor. Its vertebrate homologues later referred to as UFD2a, UBE4B or E4B were also shown to have E3 ubiquitin ligase activity. UFD2a function in the brain has been well established in vivo, and in vitro studies have shown that its activity is essential for proper condensation and segregation of chromosomes during mitosis. Here we show that 2 alternative splice forms of UFD2a, UFD2a-7 and -7/7a, are expressed sequentially during myoblast differentiation of C2C12 cell cultures and during cardiotoxin-induced regeneration of skeletal muscle in mice. UFD2a-7 contains an alternate exon 7, and UFD2a-7/7a, the larger of the 2 isoforms, contains an additional novel exon 7a. Analysis of protein or mRNA expression in mice and zebrafish revealed that a similar pattern of isoform switching occurs during developmental myogenesis of cardiac and skeletal muscle. In vertebrates (humans, rodents, zebrafish), UFD2a-7/7a is expressed only in mature striated muscle. This unique tissue specificity is further validated by the conserved presence of 2 muscle-specific splicing regulatory motifs located in the 3′ introns of exons 7 and 7a. UFD2a interacts with VCP/p97, an AAA-type ATPase implicated in processes whose functions appear to be regulated, in part, through their interaction with one or more of 15 previously identified cofactors. UFD2a-7/7a did not interact with VCP/p97 in yeast 2-hybrid experiments, which may allow the ATPase to bind cofactors that facilitate its muscle-specific functions. We conclude that the regulated expression of these UFD2a isoforms most likely imparts divergent functions that are important for myogenisis

Epistasis between COMT and MTHFR in Maternal-Fetal Dyads Increases Risk for Preeclampsia

Preeclampsia is a leading cause of perinatal morbidity and mortality. This disorder is thought to be multifactorial in origin, with multiple genes, environmental and social factors, contributing to disease. One proposed mechanism is placental hypoxia-driven imbalances in angiogenic and anti-angiogenic factors, causing endothelial cell dysfunction. Catechol-O-methyltransferase (Comt)-deficient pregnant mice have a preeclampsia phenotype that is reversed by exogenous 2-methoxyestradiol (2-ME), an estrogen metabolite generated by COMT. 2-ME inhibits Hypoxia Inducible Factor 1α, a transcription factor mediating hypoxic responses. COMT has been shown to interact with methylenetetrahydrofolate reductase (MTHFR), which modulates the availability of S-adenosylmethionine (SAM), a COMT cofactor. Variations in MTHFR have been associated with preeclampsia. By accounting for allelic variation in both genes, the role of COMT has been clarified. COMT allelic variation is linked to enzyme activity and four single nucleotide polymorphisms (SNPs) (rs6269, rs4633, rs4680, and rs4818) form haplotypes that characterize COMT activity. We tested for association between COMT haplotypes and the MTHFR 677 C→T polymorphism and preeclampsia risk in 1103 Chilean maternal-fetal dyads. The maternal ACCG COMT haplotype was associated with reduced risk for preeclampsia (P = 0.004), and that risk increased linearly from low to high activity haplotypes (P = 0.003). In fetal samples, we found that the fetal ATCA COMT haplotype and the fetal MTHFR minor “T” allele interact to increase preeclampsia risk (p = 0.022). We found a higher than expected number of patients with preeclampsia with both the fetal risk alleles alone (P = 0.052) and the fetal risk alleles in combination with a maternal balancing allele (P<0.001). This non-random distribution was not observed in controls (P = 0.341 and P = 0.219, respectively). Our findings demonstrate a role for both maternal and fetal COMT in preeclampsia and highlight the importance of including allelic variation in MTHFR

Quantitative metabolomics based on gas chromatography mass spectrometry: status and perspectives

Metabolomics involves the unbiased quantitative and qualitative analysis of the complete set of metabolites present in cells, body fluids and tissues (the metabolome). By analyzing differences between metabolomes using biostatistics (multivariate data analysis; pattern recognition), metabolites relevant to a specific phenotypic characteristic can be identified. However, the reliability of the analytical data is a prerequisite for correct biological interpretation in metabolomics analysis. In this review the challenges in quantitative metabolomics analysis with regards to analytical as well as data preprocessing steps are discussed. Recommendations are given on how to optimize and validate comprehensive silylation-based methods from sample extraction and derivatization up to data preprocessing and how to perform quality control during metabolomics studies. The current state of method validation and data preprocessing methods used in published literature are discussed and a perspective on the future research necessary to obtain accurate quantitative data from comprehensive GC-MS data is provided

Measurement of the t¯tZ and t¯tW cross sections in proton-proton collisions at √s=13 TeV with the ATLAS detector

A measurement of the associated production of a top-quark pair (t¯t) with a vector boson (W, Z) in proton-proton collisions at a center-of-mass energy of 13 TeV is presented, using 36.1  fb−1 of integrated luminosity collected by the ATLAS detector at the Large Hadron Collider. Events are selected in channels with two same- or opposite-sign leptons (electrons or muons), three leptons or four leptons, and each channel is further divided into multiple regions to maximize the sensitivity of the measurement. The t¯tZ and t¯tW production cross sections are simultaneously measured using a combined fit to all regions. The best-fit values of the production cross sections are σt¯tZ=0.95±0.08stat±0.10syst pb and σt¯tW=0.87±0.13stat±0.14syst pb in agreement with the Standard Model predictions. The measurement of the t¯tZ cross section is used to set constraints on effective field theory operators which modify the t¯tZ vertex