A strategy to detect metabolic changes induced by exposure to chemicals from large sets of condition-specific metabolic models computed with enumeration techniques.

BACKGROUND: The growing abundance of in vitro omics data, coupled with the necessity to reduce animal testing in the safety assessment of chemical compounds and even eliminate it in the evaluation of cosmetics, highlights the need for adequate computational methodologies. Data from omics technologies allow the exploration of a wide range of biological processes, therefore providing a better understanding of mechanisms of action (MoA) related to chemical exposure in biological systems. However, the analysis of these large datasets remains difficult due to the complexity of modulations spanning multiple biological processes. RESULTS: To address this, we propose a strategy to reduce information overload by computing, based on transcriptomics data, a comprehensive metabolic sub-network reflecting the metabolic impact of a chemical. The proposed strategy integrates transcriptomic data to a genome scale metabolic network through enumeration of condition-specific metabolic models hence translating transcriptomics data into reaction activity probabilities. Based on these results, a graph algorithm is applied to retrieve user readable sub-networks reflecting the possible metabolic MoA (mMoA) of chemicals. This strategy has been implemented as a three-step workflow. The first step consists in building cell condition-specific models reflecting the metabolic impact of each exposure condition while taking into account the diversity of possible optimal solutions with a partial enumeration algorithm. In a second step, we address the challenge of analyzing thousands of enumerated condition-specific networks by computing differentially activated reactions (DARs) between the two sets of enumerated possible condition-specific models. Finally, in the third step, DARs are grouped into clusters of functionally interconnected metabolic reactions, representing possible mMoA, using the distance-based clustering and subnetwork extraction method. The first part of the workflow was exemplified on eight molecules selected for their known human hepatotoxic outcomes associated with specific MoAs well described in the literature and for which we retrieved primary human hepatocytes transcriptomic data in Open TG-GATEs. Then, we further applied this strategy to more precisely model and visualize associated mMoA for two of these eight molecules (amiodarone and valproic acid). The approach proved to go beyond gene-based analysis by identifying mMoA when few genes are significantly differentially expressed (2 differentially expressed genes (DEGs) for amiodarone), bringing additional information from the network topology, or when very large number of genes were differentially expressed (5709 DEGs for valproic acid). In both cases, the results of our strategy well fitted evidence from the literature regarding known MoA. Beyond these confirmations, the workflow highlighted potential other unexplored mMoA. CONCLUSION: The proposed strategy allows toxicology experts to decipher which part of cellular metabolism is expected to be affected by the exposition to a given chemical. The approach originality resides in the combination of different metabolic modelling approaches (constraint based and graph modelling). The application to two model molecules shows the strong potential of the approach for interpretation and visual mining of complex omics in vitro data. The presented strategy is freely available as a python module ( https://pypi.org/project/manamodeller/ ) and jupyter notebooks ( https://github.com/LouisonF/MANA ).

Weight of evidence evaluation of the metabolism disrupting effects of triphenyl phosphate using an expert knowledge elicitation approach.

Identification of Endocrine-Disrupting Chemicals (EDCs) in a regulatory context requires a high level of evidence. However, lines of evidence (e.g. human, in vivo, in vitro or in silico) are heterogeneous and incomplete for quantifying evidence of the adverse effects and mechanisms involved. To date, for the regulatory appraisal of metabolism-disrupting chemicals (MDCs), no harmonised guidance to assess the weight of evidence has been developed at the EU or international level. To explore how to develop this, we applied a formal Expert Knowledge Elicitation (EKE) approach within the European GOLIATH project. EKE captures expert judgment in a quantitative manner and provides an estimate of uncertainty of the final opinion. As a proof of principle, we selected one suspected MDC -triphenyl phosphate (TPP) - based on its related adverse endpoints (obesity/adipogenicity) relevant to metabolic disruption and a putative Molecular Initiating Event (MIE): activation of peroxisome proliferator activated receptor gamma (PPARγ). We conducted a systematic literature review and assessed the quality of the lines of evidence with two independent groups of experts within GOLIATH, with the objective of categorising the metabolic disruption properties of TPP, by applying an EKE approach. Having followed the entire process separately, both groups arrived at the same conclusion, designating TPP as a "suspected MDC" with an overall quantitative agreement exceeding 85%, indicating robust reproducibility. The EKE method provides to be an important way to bring together scientists with diverse expertise and is recommended for future work in this area.

PathIntegrate: Multivariate modelling approaches for pathway-based multi-omics data integration.

As terabytes of multi-omics data are being generated, there is an ever-increasing need for methods facilitating the integration and interpretation of such data. Current multi-omics integration methods typically output lists, clusters, or subnetworks of molecules related to an outcome. Even with expert domain knowledge, discerning the biological processes involved is a time-consuming activity. Here we propose PathIntegrate, a method for integrating multi-omics datasets based on pathways, designed to exploit knowledge of biological systems and thus provide interpretable models for such studies. PathIntegrate employs single-sample pathway analysis to transform multi-omics datasets from the molecular to the pathway-level, and applies a predictive single-view or multi-view model to integrate the data. Model outputs include multi-omics pathways ranked by their contribution to the outcome prediction, the contribution of each omics layer, and the importance of each molecule in a pathway. Using semi-synthetic data we demonstrate the benefit of grouping molecules into pathways to detect signals in low signal-to-noise scenarios, as well as the ability of PathIntegrate to precisely identify important pathways at low effect sizes. Finally, using COPD and COVID-19 data we showcase how PathIntegrate enables convenient integration and interpretation of complex high-dimensional multi-omics datasets. PathIntegrate is available as an open-source Python package.

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.

PathIntegrate: Multivariate modelling approaches for pathway-based multi-omics data integration.

As terabytes of multi-omics data are being generated, there is an ever-increasing need for methods facilitating the integration and interpretation of such data. Current multi-omics integration methods typically output lists, clusters, or subnetworks of molecules related to an outcome. Even with expert domain knowledge, discerning the biological processes involved is a time-consuming activity. Here we propose PathIntegrate, a method for integrating multi-omics datasets based on pathways, designed to exploit knowledge of biological systems and thus provide interpretable models for such studies. PathIntegrate employs single-sample pathway analysis to transform multi-omics datasets from the molecular to the pathway-level, and applies a predictive single-view or multi-view model to integrate the data. Model outputs include multi-omics pathways ranked by their contribution to the outcome prediction, the contribution of each omics layer, and the importance of each molecule in a pathway. Using semi-synthetic data we demonstrate the benefit of grouping molecules into pathways to detect signals in low signal-to-noise scenarios, as well as the ability of PathIntegrate to precisely identify important pathways at low effect sizes. Finally, using COPD and COVID-19 data we showcase how PathIntegrate enables convenient integration and interpretation of complex high-dimensional multi-omics datasets. The PathIntegrate Python package is available at https://github.com/cwieder/PathIntegrate.

Challenges and perspectives for naming lipids in the context of lipidomics.

INTRODUCTION: Lipids are key compounds in the study of metabolism and are increasingly studied in biology projects. It is a very broad family that encompasses many compounds, and the name of the same compound may vary depending on the community where they are studied. OBJECTIVES: In addition, their structures are varied and complex, which complicates their analysis. Indeed, the structural resolution does not always allow a complete level of annotation so the actual compound analysed will vary from study to study and should be clearly stated. For all these reasons the identification and naming of lipids is complicated and very variable from one study to another, it needs to be harmonized. METHODS & RESULTS: In this position paper we will present and discuss the different way to name lipids (with chemoinformatic and semantic identifiers) and their importance to share lipidomic results. CONCLUSION: Homogenising this identification and adopting the same rules is essential to be able to share data within the community and to map data on functional networks.

Despite similar clinical features metabolomics reveals distinct signatures in insulin resistant and progressively obese minipigs

Obesity is a major contributor to the silent and progressive development of type 2 diabetes (T2D) whose prevention could be improved if individuals at risk were identified earlier. Our aim is to identify early phenotypes that precede T2D in diet-induced obese minipigs. We fed four groups of minipigs (n = 5–10) either normal-fat or high-fat high-sugar diet during 2, 4, or 6 months. Morphometric features were recorded, and metabolomics and clinical parameters were assessed on fasting plasma samples. Multivariate statistical analysis on 46 morphometrical and clinical parameters allowed to differentiate 4 distinct phenotypes: NFC (control group) and three others (HF2M, HF4M, HF6M) corresponding to the different stages of the obesity progression. Compared to NFC, we observed a rapid progression of body weight and fat mass (4-, 7-, and tenfold) in obese phenotypes. Insulin resistance (IR; 2.5-fold increase of HOMA-IR) and mild dyslipidemia (1.2- and twofold increase in total cholesterol and HDL) were already present in the HF2M and remained stable in HF4M and HF6M. Plasma metabolome revealed subtle changes of 23 metabolites among the obese groups, including a progressive switch in energy metabolism from amino acids to lipids, and a transient increase in de novo lipogenesis and TCA-related metabolites in HF2M. Low anti-oxidative capacities and anti-inflammatory response metabolites were found in the HF4M, and a perturbed hexose metabolism was observed in HF6M. Overall, we show that IR and progressively obese minipigs reveal phenotype-specific metabolomic signatures for which some of the identified metabolites could be considered as potential biomarkers of early progression to TD2. (AU)

Lipidomic Profiling of PFOA-Exposed Mouse Liver by Multi-Modal Mass Spectrometry Analysis.

Perfluorooctanoic acid (PFOA) is a synthetic perfluorinated chemical classified as a persistent organic pollutant. PFOA has been linked to many toxic effects, including liver injury. Many studies report that PFOA exposure alters serum and hepatic lipid metabolism. However, lipidomic pathways altered by PFOA exposure are largely unknown and only a few lipid classes, mostly triacylglycerol (TG), are usually considered in lipid analysis. Here, we performed a global lipidomic analysis on the liver of PFOA-exposed (high-dose and short-duration) and control mice by combining three mass spectrometry (MS) techniques: liquid chromatography with tandem mass spectrometry (LC-MS/MS), matrix-assisted laser desorption ionization mass spectrometry imaging (MALDI-MSI), and time-of-flight secondary ion mass spectrometry (TOF-SIMS). Among all hepatic lipids identified by LC-MS/MS analysis, more than 350 were statistically impacted (increased or decreased levels) after PFOA exposure, as confirmed by multi-variate data analysis. The levels of many lipid species from different lipid classes, most notably phosphatidylethanolamine (PE), phosphatidylcholine (PC), and TG, were significantly altered. Subsequent lipidomic analysis highlights the pathways significantly impacted by PFOA exposure, with the glycerophospholipid metabolism being the most impacted, and the changes in the lipidome network, which connects all the lipid species together. MALDI-MSI displays the heterogeneous distribution of the affected lipids and PFOA, revealing different areas of lipid expression linked to PFOA localization. TOF-SIMS localizes PFOA at the cellular level, supporting MALDI-MSI results. This multi-modal MS analysis unveils the lipidomic impact of PFOA in the mouse liver after high-dose and short-term exposure and opens new opportunities in toxicology.

Compared with Milk Protein, a Wheat and Pea Protein Blend Reduces High-Fat, High-Sucrose Induced Metabolic Dysregulations while Similarly Supporting Tissue Protein Anabolism in Rats.

BACKGROUND: Plant proteins (PPs) have been associated with better cardiovascular health than animal proteins (APs) in epidemiological studies. However, the underlying metabolic mechanisms remain mostly unknown. OBJECTIVES: Using a combination of cutting-edge isotopic methods, we aimed to better characterize the differences in protein and energy metabolisms induced by dietary protein sources (PP compared with AP) in a prudent or western dietary context. METHODS: Male Wistar rats (n = 44, 8 wk old) were fed for 4.5 mo with isoproteic diets differing in their protein isolate sources, either AP (100% milk) or PP (50%:50% pea: wheat) and being normal (NFS) or high (HFS) in sucrose (6% or 15% kcal) and saturated fat (7% or 20% kcal), respectively. We measured body weight and composition, hepatic enzyme activities and lipid content, and plasma metabolites. In the intestine, liver, adipose tissues, and skeletal muscles, we concomitantly assessed the extent of amino acid (AA) trafficking using a 15N natural abundance method, the rates of macronutrient routing to dispensable AA using a 13C natural abundance method, and the metabolic fluxes of protein synthesis (PS) and de novo lipogenesis using a 2H labeling method. Data were analyzed using ANOVA and Mixed models. RESULTS: At the whole-body level, PP limited HFS-induced insulin resistance (-27% in HOMA-IR between HFS groups, P < 0.05). In the liver, PP induced lower lipid content (-17%, P < 0.01) and de novo lipogenesis (-24%, P < 0.05). In the different tissues studied, PP induced higher AA transamination accompanied by higher routings of dietary carbohydrates and lipids toward dispensable AA synthesis by glycolysis and ß-oxidation, resulting in similar tissue PS and protein mass. CONCLUSIONS: In growing rats, compared with AP, a balanced blend of PP similarly supports protein anabolism while better limiting whole-body and tissue metabolic dysregulations through mechanisms related to their less optimal AA profile for direct channeling to PS.

Despite similar clinical features metabolomics reveals distinct signatures in insulin resistant and progressively obese minipigs.

Obesity is a major contributor to the silent and progressive development of type 2 diabetes (T2D) whose prevention could be improved if individuals at risk were identified earlier. Our aim is to identify early phenotypes that precede T2D in diet-induced obese minipigs. We fed four groups of minipigs (n = 5-10) either normal-fat or high-fat high-sugar diet during 2, 4, or 6 months. Morphometric features were recorded, and metabolomics and clinical parameters were assessed on fasting plasma samples. Multivariate statistical analysis on 46 morphometrical and clinical parameters allowed to differentiate 4 distinct phenotypes: NFC (control group) and three others (HF2M, HF4M, HF6M) corresponding to the different stages of the obesity progression. Compared to NFC, we observed a rapid progression of body weight and fat mass (4-, 7-, and tenfold) in obese phenotypes. Insulin resistance (IR; 2.5-fold increase of HOMA-IR) and mild dyslipidemia (1.2- and twofold increase in total cholesterol and HDL) were already present in the HF2M and remained stable in HF4M and HF6M. Plasma metabolome revealed subtle changes of 23 metabolites among the obese groups, including a progressive switch in energy metabolism from amino acids to lipids, and a transient increase in de novo lipogenesis and TCA-related metabolites in HF2M. Low anti-oxidative capacities and anti-inflammatory response metabolites were found in the HF4M, and a perturbed hexose metabolism was observed in HF6M. Overall, we show that IR and progressively obese minipigs reveal phenotype-specific metabolomic signatures for which some of the identified metabolites could be considered as potential biomarkers of early progression to TD2.

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.

Transgenerational metabolomic fingerprints in mice ancestrally exposed to the obesogen TBT.

BACKGROUND: Endocrine disrupting chemicals (EDCs) contribute to the etiology of metabolic disorders such as obesity, insulin resistance and hepatic dysfunction. Concern is growing about the consequences of perinatal EDC exposure on disease predisposition later in life. Metabolomics are promising approaches for studying long-term consequences of early life EDC exposure. These approaches allow for the identification and characterization of biomarkers of direct or ancestral exposures that could be diagnostic for individual susceptibility to disease and help to understand mechanisms through which EDCs act. OBJECTIVES: We sought to identify metabolomic fingerprints in mice ancestrally exposed to the model obesogen tributyltin (TBT), to assess whether metabolomics could discriminate potential trans-generational susceptibility to obesity and recognize metabolic pathways modulated by ancestral TBT exposure. METHODS: We used non-targeted 1H NMR metabolomic analyses of plasma and liver samples collected from male and female mice ancestrally exposed to TBT in two independent transgenerational experiments in which F3 and F4 males became obese when challenged with increased dietary fat. RESULTS: Metabolomics confirmed transgenerational obesogenic effects of environmentally relevant doses of TBT in F3 and F4 males, in two independent studies. Although females never became obese, their specific metabolomic fingerprint evidenced distinct transgenerational effects of TBT in female mice consistent with impaired capacity for liver biotransformation. DISCUSSION: This study is the first application of metabolomics to unveil the transgenerational effects of EDC exposure. Very early, significant changes in the plasma metabolome were observed in animals ancestrally exposed to TBT. These changes preceded the onset of obesogenic effects elicited by increased dietary fat in the TBT groups, and which ultimately resulted in significant changes in the liver metabolome. Development of metabolomic fingerprints could facilitate the identification of individuals carrying the signature of ancestral obesogen exposure that might increase their susceptibility to other risk factor such as increased dietary fat.

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.

Mitochondrial metabolism supports resistance to IDH mutant inhibitors in acute myeloid leukemia.

Mutations in IDH induce epigenetic and transcriptional reprogramming, differentiation bias, and susceptibility to mitochondrial inhibitors in cancer cells. Here, we first show that cell lines, PDXs, and patients with acute myeloid leukemia (AML) harboring an IDH mutation displayed an enhanced mitochondrial oxidative metabolism. Along with an increase in TCA cycle intermediates, this AML-specific metabolic behavior mechanistically occurred through the increase in electron transport chain complex I activity, mitochondrial respiration, and methylation-driven CEBPα-induced fatty acid ß-oxidation of IDH1 mutant cells. While IDH1 mutant inhibitor reduced 2-HG oncometabolite and CEBPα methylation, it failed to reverse FAO and OxPHOS. These mitochondrial activities were maintained through the inhibition of Akt and enhanced activation of peroxisome proliferator-activated receptor-γ coactivator-1 PGC1α upon IDH1 mutant inhibitor. Accordingly, OxPHOS inhibitors improved anti-AML efficacy of IDH mutant inhibitors in vivo. This work provides a scientific rationale for combinatory mitochondrial-targeted therapies to treat IDH mutant AML patients, especially those unresponsive to or relapsing from IDH mutant inhibitors.

DEXOM: Diversity-based enumeration of optimal context-specific metabolic networks.

The correct identification of metabolic activity in tissues or cells under different conditions can be extremely elusive due to mechanisms such as post-transcriptional modification of enzymes or different rates in protein degradation, making difficult to perform predictions on the basis of gene expression alone. Context-specific metabolic network reconstruction can overcome some of these limitations by leveraging the integration of multi-omics data into genome-scale metabolic networks (GSMN). Using the experimental information, context-specific models are reconstructed by extracting from the generic GSMN the sub-network most consistent with the data, subject to biochemical constraints. One advantage is that these context-specific models have more predictive power since they are tailored to the specific tissue, cell or condition, containing only the reactions predicted to be active in such context. However, an important limitation is that there are usually many different sub-networks that optimally fit the experimental data. This set of optimal networks represent alternative explanations of the possible metabolic state. Ignoring the set of possible solutions reduces the ability to obtain relevant information about the metabolism and may bias the interpretation of the true metabolic states. In this work we formalize the problem of enumerating optimal metabolic networks and we introduce DEXOM, an unified approach for diversity-based enumeration of context-specific metabolic networks. We developed different strategies for this purpose and we performed an exhaustive analysis using simulated and real data. In order to analyze the extent to which these results are biologically meaningful, we used the alternative solutions obtained with the different methods to measure: 1) the improvement of in silico predictions of essential genes in Saccharomyces cerevisiae using ensembles of metabolic network; and 2) the detection of alternative enriched pathways in different human cancer cell lines. We also provide DEXOM as an open-source library compatible with COBRA Toolbox 3.0, available at https://github.com/MetExplore/dexom.

An Integrated Analysis of miRNA and Gene Expression Changes in Response to an Obesogenic Diet to Explore the Impact of Transgenerational Supplementation with Omega 3 Fatty Acids.

Insulin resistance decreases the ability of insulin to inhibit hepatic gluconeogenesis, a key step in the development of metabolic syndrome. Metabolic alterations, fat accumulation, and fibrosis in the liver are closely related and contribute to the progression of comorbidities, such as hypertension, type 2 diabetes, or cancer. Omega 3 (n-3) polyunsaturated fatty acids, such as eicosapentaenoic acid (EPA), were identified as potent positive regulators of insulin sensitivity in vitro and in animal models. In the current study, we explored the effects of a transgenerational supplementation with EPA in mice exposed to an obesogenic diet on the regulation of microRNAs (miRNAs) and gene expression in the liver using high-throughput techniques. We implemented a comprehensive molecular systems biology approach, combining statistical tools, such as MicroRNA Master Regulator Analysis pipeline and Boolean modeling to integrate these biochemical processes. We demonstrated that EPA mediated molecular adaptations, leading to the inhibition of miR-34a-5p, a negative regulator of Irs2 as a master regulatory event leading to the inhibition of gluconeogenesis by insulin during the fasting-feeding transition. Omics data integration provided greater biological insight and a better understanding of the relationships between biological variables. Such an approach may be useful for deriving innovative data-driven hypotheses and for the discovery of molecular-biochemical mechanistic links.

Postprandial NMR-Based Metabolic Exchanges Reflect Impaired Phenotypic Flexibility across Splanchnic Organs in the Obese Yucatan Mini-Pig.

The postprandial period represents one of the most challenging phenomena in whole-body metabolism, and it can be used as a unique window to evaluate the phenotypic flexibility of an individual in response to a given meal, which can be done by measuring the resilience of the metabolome. However, this exploration of the metabolism has never been applied to the arteriovenous (AV) exploration of organs metabolism. Here, we applied an AV metabolomics strategy to evaluate the postprandial flexibility across the liver and the intestine of mini-pigs subjected to a high fat-high sucrose (HFHS) diet for 2 months. We identified for the first time a postprandial signature associated to the insulin resistance and obesity outcomes, and we showed that the splanchnic postprandial metabolome was considerably affected by the meal and the obesity condition. Most of the changes induced by obesity were observed in the exchanges across the liver, where the metabolism was reorganized to maintain whole body glucose homeostasis by routing glucose formed de novo from a large variety of substrates into glycogen. Furthermore, metabolites related to lipid handling and energy metabolism showed a blunted postprandial response in the obese animals across organs. Finally, some of our results reflect a loss of flexibility in response to the HFHS meal challenge in unsuspected metabolic pathways that must be further explored as potential new events involved in early obesity and the onset of insulin resistance.

An Optimized Dual Extraction Method for the Simultaneous and Accurate Analysis of Polar Metabolites and Lipids Carried out on Single Biological Samples.

The functional understanding of metabolic changes requires both a significant investigation into metabolic pathways, as enabled by global metabolomics and lipidomics approaches, and the comprehensive and accurate exploration of specific key pathways. To answer this pivotal challenge, we propose an optimized approach, which combines an efficient sample preparation, aiming to reduce the variability, with a biphasic extraction method, where both the aqueous and organic phases of the same sample are used for mass spectrometry analyses. We demonstrated that this double extraction protocol allows working with one single sample without decreasing the metabolome and lipidome coverage. It enables the targeted analysis of 40 polar metabolites and 82 lipids, together with the absolute quantification of 32 polar metabolites, providing comprehensive coverage and quantitative measurement of the metabolites involved in central carbon energy pathways. With this method, we evidenced modulations of several lipids, amino acids, and energy metabolites in HepaRG cells exposed to fenofibrate, a model hepatic toxicant, and metabolic modulator. This new protocol is particularly relevant for experiments involving limited amounts of biological material and for functional metabolic explorations and is thus of particular interest for studies aiming to decipher the effects and modes of action of metabolic disrupting compounds.

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.