Genome scale metabolic network modelling for metabolic profile predictions.

Metabolic profiling (metabolomics) aims at measuring small molecules (metabolites) in complex samples like blood or urine for human health studies. While biomarker-based assessment often relies on a single molecule, metabolic profiling combines several metabolites to create a more complex and more specific fingerprint of the disease. However, in contrast to genomics, there is no unique metabolomics setup able to measure the entire metabolome. This challenge leads to tedious and resource consuming preliminary studies to be able to design the right metabolomics experiment. In that context, computer assisted metabolic profiling can be of strong added value to design metabolomics studies more quickly and efficiently. We propose a constraint-based modelling approach which predicts in silico profiles of metabolites that are more likely to be differentially abundant under a given metabolic perturbation (e.g. due to a genetic disease), using flux simulation. In genome-scale metabolic networks, the fluxes of exchange reactions, also known as the flow of metabolites through their external transport reactions, can be simulated and compared between control and disease conditions in order to calculate changes in metabolite import and export. These import/export flux differences would be expected to induce changes in circulating biofluid levels of those metabolites, which can then be interpreted as potential biomarkers or metabolites of interest. In this study, we present SAMBA (SAMpling Biomarker Analysis), an approach which simulates fluxes in exchange reactions following a metabolic perturbation using random sampling, compares the simulated flux distributions between the baseline and modulated conditions, and ranks predicted differentially exchanged metabolites as potential biomarkers for the perturbation. We show that there is a good fit between simulated metabolic exchange profiles and experimental differential metabolites detected in plasma, such as patient data from the disease database OMIM, and metabolic trait-SNP associations found in mGWAS studies. These biomarker recommendations can provide insight into the underlying mechanism or metabolic pathway perturbation lying behind observed metabolite differential abundances, and suggest new metabolites as potential avenues for further experimental analyses.

Suggesting disease associations for overlooked metabolites using literature from metabolic neighbors.

In human health research, metabolic signatures extracted from metabolomics data have a strong added value for stratifying patients and identifying biomarkers. Nevertheless, one of the main challenges is to interpret and relate these lists of discriminant metabolites to pathological mechanisms. This task requires experts to combine their knowledge with information extracted from databases and the scientific literature. However, we show that most compounds (>99%) in the PubChem database lack annotated literature. This dearth of available information can have a direct impact on the interpretation of metabolic signatures, which is often restricted to a subset of significant metabolites. To suggest potential pathological phenotypes related to overlooked metabolites that lack annotated literature, we extend the "guilt-by-association" principle to literature information by using a Bayesian framework. The underlying assumption is that the literature associated with the metabolic neighbors of a compound can provide valuable insights, or an a priori, into its biomedical context. The metabolic neighborhood of a compound can be defined from a metabolic network and correspond to metabolites to which it is connected through biochemical reactions. With the proposed approach, we suggest more than 35,000 associations between 1,047 overlooked metabolites and 3,288 diseases (or disease families). All these newly inferred associations are freely available on the FORUM ftp server (see information at https://github.com/eMetaboHUB/Forum-LiteraturePropagation).

Pathway analysis in metabolomics: Recommendations for the use of over-representation analysis.

Over-representation analysis (ORA) is one of the commonest pathway analysis approaches used for the functional interpretation of metabolomics datasets. Despite the widespread use of ORA in metabolomics, the community lacks guidelines detailing its best-practice use. Many factors have a pronounced impact on the results, but to date their effects have received little systematic attention. Using five publicly available datasets, we demonstrated that changes in parameters such as the background set, differential metabolite selection methods, and pathway database used can result in profoundly different ORA results. The use of a non-assay-specific background set, for example, resulted in large numbers of false-positive pathways. Pathway database choice, evaluated using three of the most popular metabolic pathway databases (KEGG, Reactome, and BioCyc), led to vastly different results in both the number and function of significantly enriched pathways. Factors that are specific to metabolomics data, such as the reliability of compound identification and the chemical bias of different analytical platforms also impacted ORA results. Simulated metabolite misidentification rates as low as 4% resulted in both gain of false-positive pathways and loss of truly significant pathways across all datasets. Our results have several practical implications for ORA users, as well as those using alternative pathway analysis methods. We offer a set of recommendations for the use of ORA in metabolomics, alongside a set of minimal reporting guidelines, as a first step towards the standardisation of pathway analysis in metabolomics.

FORUM: building a Knowledge Graph from public databases and scientific literature to extract associations between chemicals and diseases.

MOTIVATION: Metabolomics studies aim at reporting a metabolic signature (list of metabolites) related to a particular experimental condition. These signatures are instrumental in the identification of biomarkers or classification of individuals, however their biological and physiological interpretation remains a challenge. To support this task, we introduce FORUM: a Knowledge Graph (KG) providing a semantic representation of relations between chemicals and biomedical concepts, built from a federation of life science databases and scientific literature repositories. RESULTS: The use of a Semantic Web framework on biological data allows us to apply ontological-based reasoning to infer new relations between entities. We show that these new relations provide different levels of abstraction and could open the path to new hypotheses. We estimate the statistical relevance of each extracted relation, explicit or inferred, using an enrichment analysis, and instantiate them as new knowledge in the KG to support results interpretation/further inquiries. AVAILABILITY AND IMPLEMENTATION: A web interface to browse and download the extracted relations, as well as a SPARQL endpoint to directly probe the whole FORUM KG, are available at https://forum-webapp.semantic-metabolomics.fr. The code needed to reproduce the triplestore is available at https://github.com/eMetaboHUB/Forum-DiseasesChem. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Multiplatform metabolomics for an integrative exploration of metabolic syndrome in older men.

BACKGROUND: Metabolic syndrome (MetS), a cluster of factors associated with risks of developing cardiovascular diseases, is a public health concern because of its growing prevalence. Considering the combination of concomitant components, their development and severity, MetS phenotypes are largely heterogeneous, inducing disparity in diagnosis. METHODS: A case/control study was designed within the NuAge longitudinal cohort on aging. From a 3-year follow-up of 123 stable individuals, we present a deep phenotyping approach based on a multiplatform metabolomics and lipidomics untargeted strategy to better characterize metabolic perturbations in MetS and define a comprehensive MetS signature stable over time in older men. FINDINGS: We characterize significant changes associated with MetS, involving modulations of 476 metabolites and lipids, and representing 16% of the detected serum metabolome/lipidome. These results revealed a systemic alteration of metabolism, involving various metabolic pathways (urea cycle, amino-acid, sphingo- and glycerophospholipid, and sugar metabolisms…) not only intrinsically interrelated, but also reflecting environmental factors (nutrition, microbiota, physical activity…). INTERPRETATION: These findings allowed identifying a comprehensive MetS signature, reduced to 26 metabolites for future translation into clinical applications for better diagnosing MetS. FUNDING: The NuAge Study was supported by a research grant from the Canadian Institutes of Health Research (CIHR; MOP-62842). The actual NuAge Database and Biobank, containing data and biologic samples of 1,753 NuAge participants (from the initial 1,793 participants), are supported by the Fonds de recherche du Québec (FRQ; 2020-VICO-279753), the Quebec Network for Research on Aging, a thematic network funded by the Fonds de Recherche du Québec - Santé (FRQS) and by the Merck-Frost Chair funded by La Fondation de l'Université de Sherbrooke. All metabolomics and lipidomics analyses were funded and performed within the metaboHUB French infrastructure (ANR-INBS-0010). All authors had full access to the full data in the study and accept responsibility to submit for publication.

Improving lipid mapping in Genome Scale Metabolic Networks using ontologies.

INTRODUCTION: To interpret metabolomic and lipidomic profiles, it is necessary to identify the metabolic reactions that connect the measured molecules. This can be achieved by putting them in the context of genome-scale metabolic network reconstructions. However, mapping experimentally measured molecules onto metabolic networks is challenging due to differences in identifiers and level of annotation between data and metabolic networks, especially for lipids. OBJECTIVES: To help linking lipids from lipidomics datasets with lipids in metabolic networks, we developed a new matching method based on the ChEBI ontology. The implementation is freely available as a python library and in MetExplore webserver. METHODS: Our matching method is more flexible than an exact identifier-based correspondence since it allows establishing a link between molecules even if a different level of precision is provided in the dataset and in the metabolic network. For instance, it can associate a generic class of lipids present in the network with the molecular species detailed in the lipidomics dataset. This mapping is based on the computation of a distance between molecules in ChEBI ontology. RESULTS: We applied our method to a chemical library (968 lipids) and an experimental dataset (32 modulated lipids) and showed that using ontology-based mapping improves and facilitates the link with genome scale metabolic networks. Beyond network mapping, the results provide ways for improvements in terms of network curation and lipidomics data annotation. CONCLUSION: This new method being generic, it can be applied to any metabolomics data and therefore improve our comprehension of metabolic modulations.

MetaboRank: network-based recommendation system to interpret and enrich metabolomics results.

Motivation: Metabolomics has shown great potential to improve the understanding of complex diseases, potentially leading to therapeutic target identification. However, no single analytical method allows monitoring all metabolites in a sample, resulting in incomplete metabolic fingerprints. This incompleteness constitutes a stumbling block to interpretation, raising the need for methods that can enrich those fingerprints. We propose MetaboRank, a new solution inspired by social network recommendation systems for the identification of metabolites potentially related to a metabolic fingerprint. Results: MetaboRank method had been used to enrich metabolomics data obtained on cerebrospinal fluid samples from patients suffering from hepatic encephalopathy (HE). MetaboRank successfully recommended metabolites not present in the original fingerprint. The quality of recommendations was evaluated by using literature automatic search, in order to check that recommended metabolites could be related to the disease. Complementary mass spectrometry experiments and raw data analysis were performed to confirm these suggestions. In particular, MetaboRank recommended the overlooked α-ketoglutaramate as a metabolite which should be added to the metabolic fingerprint of HE, thus suggesting that metabolic fingerprints enhancement can provide new insight on complex diseases. Availability and implementation: Method is implemented in the MetExplore server and is available at www.metexplore.fr. A tutorial is available at https://metexplore.toulouse.inra.fr/com/tutorials/MetaboRank/2017-MetaboRank.pdf. Supplementary information: Supplementary data are available at Bioinformatics online.

Large-Scale Modeling Approach Reveals Functional Metabolic Shifts during Hepatic Differentiation.

Being able to explore the metabolism of broad metabolizing cells is of critical importance in many research fields. This article presents an original modeling solution combining metabolic network and omics data to identify modulated metabolic pathways and changes in metabolic functions occurring during differentiation of a human hepatic cell line (HepaRG). Our results confirm the activation of hepato-specific functionalities and newly evidence modulation of other metabolic pathways, which could not be evidenced from transcriptomic data alone. Our method takes advantage of the network structure to detect changes in metabolic pathways that do not have gene annotations and exploits flux analyses techniques to identify activated metabolic functions. Compared to the usual cell-specific metabolic network reconstruction approaches, it limits false predictions by considering several possible network configurations to represent one phenotype rather than one arbitrarily selected network. Our approach significantly enhances the comprehensive and functional assessment of cell metabolism, opening further perspectives to investigate metabolic shifts occurring within various biological contexts.

MetExplore: collaborative edition and exploration of metabolic networks.

Metabolism of an organism is composed of hundreds to thousands of interconnected biochemical reactions responding to environmental or genetic constraints. This metabolic network provides a rich knowledge to contextualize omics data and to elaborate hypotheses on metabolic modulations. Nevertheless, performing this kind of integrative analysis is challenging for end users with not sufficiently advanced computer skills since it requires the use of various tools and web servers. MetExplore offers an all-in-one online solution composed of interactive tools for metabolic network curation, network exploration and omics data analysis. In particular, it is possible to curate and annotate metabolic networks in a collaborative environment. The network exploration is also facilitated in MetExplore by a system of interactive tables connected to a powerful network visualization module. Finally, the contextualization of metabolic elements in the network and the calculation of over-representation statistics make it possible to interpret any kind of omics data. MetExplore is a sustainable project maintained since 2010 freely available at https://metexplore.toulouse.inra.fr/metexplore2/.

MetExploreViz: web component for interactive metabolic network visualization.

SUMMARY: MetExploreViz is an open source web component that can be easily embedded in any web site. It provides features dedicated to the visualization of metabolic networks and pathways and thus offers a flexible solution to analyse omics data in a biochemical context. AVAILABILITY AND IMPLEMENTATION: Documentation and link to GIT code repository (GPL 3.0 license) are available at this URL: http://metexplore.toulouse.inra.fr/metexploreViz/doc/.

A Computational Solution to Automatically Map Metabolite Libraries in the Context of Genome Scale Metabolic Networks.

This article describes a generic programmatic method for mapping chemical compound libraries on organism-specific metabolic networks from various databases (KEGG, BioCyc) and flat file formats (SBML and Matlab files). We show how this pipeline was successfully applied to decipher the coverage of chemical libraries set up by two metabolomics facilities MetaboHub (French National infrastructure for metabolomics and fluxomics) and Glasgow Polyomics (GP) on the metabolic networks available in the MetExplore web server. The present generic protocol is designed to formalize and reduce the volume of information transfer between the library and the network database. Matching of metabolites between libraries and metabolic networks is based on InChIs or InChIKeys and therefore requires that these identifiers are specified in both libraries and networks. In addition to providing covering statistics, this pipeline also allows the visualization of mapping results in the context of metabolic networks. In order to achieve this goal, we tackled issues on programmatic interaction between two servers, improvement of metabolite annotation in metabolic networks and automatic loading of a mapping in genome scale metabolic network analysis tool MetExplore. It is important to note that this mapping can also be performed on a single or a selection of organisms of interest and is thus not limited to large facilities.

Dynamic Metabolic Disruption in Rats Perinatally Exposed to Low Doses of Bisphenol-A.

Along with the well-established effects on fertility and fecundity, perinatal exposure to endocrine disrupting chemicals, and notably to xeno-estrogens, is strongly suspected of modulating general metabolism. The metabolism of a perinatally exposed individual may be durably altered leading to a higher susceptibility of developing metabolic disorders such as obesity and diabetes; however, experimental designs involving the long term study of these dynamic changes in the metabolome raise novel challenges. 1H-NMR-based metabolomics was applied to study the effects of bisphenol-A (BPA, 0; 0.25; 2.5, 25 and 250 µg/kg BW/day) in rats exposed perinatally. Serum and liver samples of exposed animals were analyzed on days 21, 50, 90, 140 and 200 in order to explore whether maternal exposure to BPA alters metabolism. Partial Least Squares-Discriminant Analysis (PLS-DA) was independently applied to each time point, demonstrating a significant pair-wise discrimination for liver as well as serum samples at all time-points, and highlighting unequivocal metabolic shifts in rats perinatally exposed to BPA, including those exposed to lower doses. In BPA exposed animals, metabolism of glucose, lactate and fatty acids was modified over time. To further explore dynamic variation, ANOVA-Simultaneous Component Analysis (A-SCA) was used to separate data into blocks corresponding to the different sources of variation (Time, Dose and Time*Dose interaction). A-SCA enabled the demonstration of a dynamic, time/age dependent shift of serum metabolome throughout the rats' lifetimes. Variables responsible for the discrimination between groups clearly indicate that BPA modulates energy metabolism, and suggest alterations of neurotransmitter signaling, the latter finding being compatible with the neurodevelopmental effect of this xenoestrogen. In conclusion, long lasting metabolic effects of BPA could be characterized over 200 days, despite physiological (and thus metabolic) changes connected with sexual maturation and aging.

TrypanoCyc: a community-led biochemical pathways database for Trypanosoma brucei.

The metabolic network of a cell represents the catabolic and anabolic reactions that interconvert small molecules (metabolites) through the activity of enzymes, transporters and non-catalyzed chemical reactions. Our understanding of individual metabolic networks is increasing as we learn more about the enzymes that are active in particular cells under particular conditions and as technologies advance to allow detailed measurements of the cellular metabolome. Metabolic network databases are of increasing importance in allowing us to contextualise data sets emerging from transcriptomic, proteomic and metabolomic experiments. Here we present a dynamic database, TrypanoCyc (http://www.metexplore.fr/trypanocyc/), which describes the generic and condition-specific metabolic network of Trypanosoma brucei, a parasitic protozoan responsible for human and animal African trypanosomiasis. In addition to enabling navigation through the BioCyc-based TrypanoCyc interface, we have also implemented a network-based representation of the information through MetExplore, yielding a novel environment in which to visualise the metabolism of this important parasite.

Exposure to pesticides and risk of childhood cancer: a meta-analysis of recent epidemiological studies.

OBJECTIVES: The authors performed a meta-analysis of case-control and cohort studies to clarify the possible relationship between exposure to pesticides and childhood cancers. METHODS: Two cohort and 38 case-control studies were selected for the first meta-analysis. After evaluating homogeneity among studies using the Cochran Q test, the authors calculated a pooled meta-OR stratified on each cancer site. The authors then constructed a list of variables believed to play an important role in explaining the relation between parental exposure to pesticide and childhood cancer, and performed a series of meta-analyses. The authors also performed a distinct meta-analysis for three cohort studies with RR data. RESULTS: Meta-analysis of the three cohort studies did not show any positive links between parental pesticide exposure and childhood cancer incidence. However, the meta-analysis of the 40 studies with OR values showed that the risk of lymphoma and leukaemia increased significantly in exposed children when their mother was exposed during the prenatal period (OR=1.53; 95% CI 1.22 to 1.91 and OR=1.48; 95% CI 1.26 to 1.75). The risk of brain cancer was correlated with paternal exposure either before or after birth (OR=1.49; 95% CI 1.23 to 1.79 and OR=1.66; 95% CI 1.11 to 2.49). The OR of leukaemia and lymphoma was higher when the mother was exposed to pesticides (through household use or professional exposure). Conversely, the incidence of brain cancer was influenced by the father's exposure (occupational activity or use of household or garden pesticides). CONCLUSION: Despite some limitations in this study, the incidence of childhood cancer does appear to be associated with parental exposure during the prenatal period.

MetExplore: a web server to link metabolomic experiments and genome-scale metabolic networks.

High-throughput metabolomic experiments aim at identifying and ultimately quantifying all metabolites present in biological systems. The metabolites are interconnected through metabolic reactions, generally grouped into metabolic pathways. Classical metabolic maps provide a relational context to help interpret metabolomics experiments and a wide range of tools have been developed to help place metabolites within metabolic pathways. However, the representation of metabolites within separate disconnected pathways overlooks most of the connectivity of the metabolome. By definition, reference pathways cannot integrate novel pathways nor show relationships between metabolites that may be linked by common neighbours without being considered as joint members of a classical biochemical pathway. MetExplore is a web server that offers the possibility to link metabolites identified in untargeted metabolomics experiments within the context of genome-scale reconstructed metabolic networks. The analysis pipeline comprises mapping metabolomics data onto the specific metabolic network of an organism, then applying graph-based methods and advanced visualization tools to enhance data analysis. The MetExplore web server is freely accessible at http://metexplore.toulouse.inra.fr.

High-throughput single nucleotide polymorphism genotyping in wheat (Triticum spp.).

Over the past few years, considerable progress has been made in high-throughput single nucleotide polymorphism (SNP) genotyping technologies, largely through the investment of the human genetics community. These technologies are well adapted to diploid species. For plant breeding purposes, it is important to determine whether these genotyping methods are adapted to polyploidy, as most major crops are former or recent polyploids. To address this problem, we tested the capacity of the multiplex technology SNPlex with a set of 47 wheat SNPs to genotype DNAs of 1314 lines that were organized in four 384-well plates. These lines represented different taxa of tetra- and hexaploid Triticum species and their wild diploid relatives. We observed 40 markers which gave less than 20% missing data. Different methods, based on either Sanger sequencing or the MassARRAY genotyping technology, were then used to validate the genotypes obtained by SNPlex for 11 markers. The concordance of the genotypes obtained by SNPlex with the results obtained by the different validation methods was 96%, except for one discarded marker. Furthermore, a mapping study on six markers showed the expected genetic positions previously described. To conclude, this study showed that high-throughput genotyping technologies developed for diploid species can be used successfully in polyploids, although there is a need for manual reading. For the first time in wheat species, a core of 39 SNPs is available that can serve as the basis for the development of a complete SNPlex set of 48 markers.

A worldwide bread wheat core collection arrayed in a 384-well plate.

Bread wheat (Triticum aestivum), one of the world's major crops, is genetically very diverse. In order to select a representative sample of the worldwide wheat diversity, 3,942 accessions originating from 73 countries were analysed with a set of 38 genomic simple sequence repeat (SSR) markers. The number of alleles at each locus ranged from 7 to 45 with an average of 23.9 alleles per locus. The 908 alleles detected were used together with passport data to select increasingly large sub-samples that maximised both the number of observed alleles at SSR loci and the number of geographical origins. A final core of 372 accessions (372CC) was selected with this M strategy. All the different geographical areas and more than 98% of the allelic diversity at the 38 polymorphic loci were represented in this core. The method used to build the core was validated, by using a second set of independent markers [44 expressed sequence tag (EST)-SSR markers] on a larger sample of 744 accessions: 96.74% of the alleles observed at these loci had already been captured in the 372CC. So maximizing the diversity with a first set of markers also maximised the diversity at a second independent set of locus. To relate the genetic structure of wheat germplasm to its geographical origins, the two sets of markers were used to compute a dissimilarity matrix between geographical groups. Current worldwide wheat diversity is clearly divided according to wheat's European and Asian origins, whereas the diversity within each geographical group might be the result of the combined effects of adaptation of an initial germplasm to different environmental conditions and specific breeding practices. Seeds from each accession of the 372CC were multiplied and are now available to the scientific community. The genomic DNA of the 372CC, which can be entirely contained in a 384-deep-well storage plate, will be a useful tool for future studies of wheat genetic diversity.