Demonstrating the reliability of in vivo metabolomics based chemical grouping: towards best practice

While grouping/read-across is widely used to fill data gaps, chemical registration dossiers are often rejected due to weak category justifications based on structural similarity only. Metabolomics provides a route to robust chemical categories via evidence of shared molecular effects across source and target substances. To gain international acceptance, this approach must demonstrate high reliability, and best-practice guidance is required. The MetAbolomics ring Trial for CHemical groupING (MATCHING), comprising six industrial, government and academic ring-trial partners, evaluated inter-laboratory reproducibility and worked towards best-practice. An independent team selected eight substances (WY-14643, 4-chloro-3-nitroaniline, 17α-methyl-testosterone, trenbolone, aniline, dichlorprop-p, 2-chloroaniline, fenofibrate); ring-trial partners were blinded to their identities and modes-of-action. Plasma samples were derived from 28-day rat tests (two doses per substance), aliquoted, and distributed to partners. Each partner applied their preferred liquid chromatography-mass spectrometry (LC-MS) metabolomics workflows to acquire, process, quality assess, statistically analyze and report their grouping results to the European Chemicals Agency, to ensure the blinding conditions of the ring trial. Five of six partners, whose metabolomics datasets passed quality control, correctly identified the grouping of eight test substances into three categories, for both male and female rats. Strikingly, this was achieved even though a range of metabolomics approaches were used. Through assessing intrastudy quality-control samples, the sixth partner observed high technical variation and was unable to group the substances. By comparing workflows, we conclude that some heterogeneity in metabolomics methods is not detrimental to consistent grouping, and that assessing data quality prior to grouping is essential. We recommend development of international guidance for quality-control acceptance criteria. This study demonstrates the reliability of metabolomics for chemical grouping and works towards best-practice.

Lipidome characterisation and sex-specific differences in type 1 and type 2 diabetes mellitus

BACKGROUND

In this study, we evaluated the lipidome alterations caused by type 1 diabetes (T1D) and type 2 diabetes (T2D), by determining lipids significantly associated with diabetes overall and in both sexes, and lipids associated with the glycaemic state.

METHODS

An untargeted lipidomic analysis was performed to measure the lipid profiles of 360 subjects (91 T1D, 91 T2D, 74 with prediabetes and 104 controls (CT)) without cardiovascular and/or chronic kidney disease. Ultra-high performance liquid chromatography-electrospray ionization mass spectrometry (UHPLC-ESI-MS) was conducted in two ion modes (positive and negative). We used multiple linear regression models to (1) assess the association between each lipid feature and each condition, (2) determine sex-specific differences related to diabetes, and (3) identify lipids associated with the glycaemic state by considering the prediabetes stage. The models were adjusted by sex, age, hypertension, dyslipidaemia, body mass index, glucose, smoking, systolic blood pressure, triglycerides, HDL cholesterol, LDL cholesterol, alternate Mediterranean diet score (aMED) and estimated glomerular filtration rate (eGFR); diabetes duration and glycated haemoglobin (HbA1c) were also included in the comparison between T1D and T2D.

RESULTS

A total of 54 unique lipid subspecies from 15 unique lipid classes were annotated. Lysophosphatidylcholines (LPC) and ceramides (Cer) showed opposite effects in subjects with T1D and subjects with T2D, LPCs being mainly up-regulated in T1D and down-regulated in T2D, and Cer being up-regulated in T2D and down-regulated in T1D. Also, Phosphatidylcholines were clearly down-regulated in subjects with T1D. Regarding sex-specific differences, ceramides and phosphatidylcholines exhibited important diabetes-associated differences due to sex. Concerning the glycaemic state, we found a gradual increase of a panel of 1-deoxyceramides from normoglycemia to prediabetes to T2D.

CONCLUSIONS

Our findings revealed an extensive disruption of lipid metabolism in both T1D and T2D. Additionally, we found sex-specific lipidome changes associated with diabetes, and lipids associated with the glycaemic state that can be linked to previously described molecular mechanisms in diabetes.

Simultaneously discovering the fate and biochemical effects of pharmaceuticals through untargeted metabolomics

Untargeted metabolomics is an established approach in toxicology for characterising endogenous metabolic responses to xenobiotic exposure. Detecting the xenobiotic and its biotransformation products as part of the metabolomics analysis provides an opportunity to simultaneously gain deep insights into its fate and metabolism, and to associate the internal relative dose directly with endogenous metabolic responses. This integration of untargeted exposure and response measurements into a single assay has yet to be fully demonstrated. Here we assemble a workflow to discover and analyse pharmaceutical-related measurements from routine untargeted UHPLC-MS metabolomics datasets, derived from in vivo (rat plasma and cardiac tissue, and human plasma) and in vitro (human cardiomyocytes) studies that were principally designed to investigate endogenous metabolic responses to drug exposure. Our findings clearly demonstrate how untargeted metabolomics can discover extensive biotransformation maps, temporally-changing relative systemic exposure, and direct associations of endogenous biochemical responses to the internal dose.

The Precision Toxicology Initiative

The goal of PrecisionTox is to overcome conceptual barriers to replacing traditional mammalian chemical safety testing by accelerating the discovery of evolutionarily conserved toxicity pathways that are shared by descent among humans and more distantly related animals. An international consortium is systematically testing the toxicological effects of a diverse set of chemicals on a suite of five model species comprising fruit flies, nematodes, water fleas, and embryos of clawed frogs and zebrafish along with human cell lines. Multiple forms of omics and comparative toxicology data are integrated to map the evolutionary origins of biomolecular interactions, which are predictive of adverse health effects, to major branches of the animal phylogeny. These conserved elements of adverse outcome pathways (AOPs) and their biomarkers are expect to provide mechanistic insight useful for regulating groups of chemicals based on their shared modes of action. PrecisionTox also aims to quantify risk variation within populations by recognizing susceptibility as a heritable trait that varies with genetic diversity. This initiative incorporates legal experts and collaborates with risk managers to address specific needs within European chemicals legislation, including the uptake of new approach methodologies (NAMs) for setting precise regulatory limits on toxic chemicals.

The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2022 update

Galaxy is a mature, browser accessible workbench for scientific computing. It enables scientists to share, analyze and visualize their own data, with minimal technical impediments. A thriving global community continues to use, maintain and contribute to the project, with support from multiple national infrastructure providers that enable freely accessible analysis and training services. The Galaxy Training Network supports free, self-directed, virtual training with >230 integrated tutorials. Project engagement metrics have continued to grow over the last 2 years, including source code contributions, publications, software packages wrapped as tools, registered users and their daily analysis jobs, and new independent specialized servers. Key Galaxy technical developments include an improved user interface for launching large-scale analyses with many files, interactive tools for exploratory data analysis, and a complete suite of machine learning tools. Important scientific developments enabled by Galaxy include Vertebrate Genome Project (VGP) assembly workflows and global SARS-CoV-2 collaborations.

Suitability of Dried Blood Spots for Accelerating Veterinary Biobank Collections and Identifying Metabolomics Biomarkers With Minimal Resources

Biomarker discovery using biobank samples collected from veterinary clinics would deliver insights into the diverse population of pets and accelerate diagnostic development. The acquisition, preparation, processing, and storage of biofluid samples in sufficient volumes and at a quality suitable for later analysis with most suitable discovery methods remain challenging. Metabolomics analysis is a valuable approach to detect health/disease phenotypes. Pre-processing changes during preparation of plasma/serum samples may induce variability that may be overcome using dried blood spots (DBSs). We report a proof of principle study by metabolite fingerprinting applying UHPLC-MS of plasma and DBSs acquired from healthy adult dogs and cats (age range 1–9 years), representing each of 4 dog breeds (Labrador retriever, Beagle, Petit Basset Griffon Vendeen, and Norfolk terrier) and the British domestic shorthair cat (n = 10 per group). Blood samples (20 and 40 μL) for DBSs were loaded onto filter paper, air-dried at room temperature (3 h), and sealed and stored (4°C for ~72 h) prior to storage at −80°C. Plasma from the same blood draw (250 μL) was prepared and stored at −80°C within 1 h of sampling. Metabolite fingerprinting of the DBSs and plasma produced similar numbers of metabolite features that had similar abilities to discriminate between biological classes and correctly assign blinded samples. These provide evidence that DBSs, sampled in a manner amenable to application in in-clinic/in-field processing, are a suitable sample for biomarker discovery using UHPLC-MS metabolomics. Further, given appropriate owner consent, the volumes tested (20–40 μL) make the acquisition of remnant blood from blood samples drawn for other reasons available for biobanking and other research activities. Together, this makes possible large-scale biobanking of veterinary samples, gaining sufficient material sooner and enabling quicker identification of biomarkers of interest.

Integrating in vitro metabolomics with a 96-well high-throughput screening platform

INTRODUCTION

High-throughput screening (HTS) is emerging as an approach to support decision-making in chemical safety assessments. In parallel, in vitro metabolomics is a promising approach that can help accelerate the transition from animal models to high-throughput cell-based models in toxicity testing.

OBJECTIVE

In this study we establish and evaluate a high-throughput metabolomics workflow that is compatible with a 96-well HTS platform employing 50,000 hepatocytes of HepaRG per well.

METHODS

Low biomass cell samples were extracted for metabolomics analyses using a newly established semi-automated protocol, and the intracellular metabolites were analysed using a high-resolution spectral-stitching nanoelectrospray direct infusion mass spectrometry (nESI-DIMS) method that was modified for low sample biomass.

RESULTS

The method was assessed with respect to sensitivity and repeatability of the entire workflow from cell culturing and sampling to measurement of the metabolic phenotype, demonstrating sufficient sensitivity (> 3000 features in hepatocyte extracts) and intra- and inter-plate repeatability for polar nESI-DIMS assays (median relative standard deviation < 30%). The assays were employed for a proof-of-principle toxicological study with a model toxicant, cadmium chloride, revealing changes in the metabolome across five sampling times in the 48-h exposure period. To allow the option for lipidomics analyses, the solvent system was extended by establishing separate extraction methods for polar metabolites and lipids.

CONCLUSIONS

Experimental, analytical and informatics workflows reported here met pre-defined criteria in terms of sensitivity, repeatability and ability to detect metabolome changes induced by a toxicant and are ready for application in metabolomics-driven toxicity testing to complement HTS assays.

Automated Sample Preparation and Data Collection Workflow for High-Throughput In Vitro Metabolomics

Regulatory bodies have started to recognise the value of in vitro screening and metabolomics as two types of new approach methodologies (NAMs) for chemical risk assessments, yet few high-throughput in vitro toxicometabolomics studies have been reported. A significant challenge is to implement automated sample preparation of the low biomass samples typically used for in vitro screening. Building on previous work, we have developed, characterised and demonstrated an automated sample preparation and analysis workflow for in vitro metabolomics of HepaRG cells in 96-well microplates using a Biomek i7 Hybrid Workstation (Beckman Coulter) and Orbitrap Elite (Thermo Scientific) high-resolution nanoelectrospray direct infusion mass spectrometry (nESI-DIMS), across polar metabolites and lipids. The experimental conditions evaluated included the day of metabolite extraction, order of extraction of samples in 96-well microplates, position of the 96-well microplate on the instrument's deck and well location within a microplate. By using the median relative standard deviation (mRSD (%)) of spectral features, we have demonstrated good repeatability of the workflow (final mRSD < 30%) with a low percentage of features outside the threshold applied for statistical analysis. To improve the quality of the automated workflow further, small method modifications were made and then applied to a large cohort study (4860 sample infusions across three nESI-DIMS assays), which confirmed very high repeatability of the whole workflow from cell culturing to metabolite measurements, whilst providing a significant improvement in sample throughput. It is envisioned that the automated in vitro metabolomics workflow will help to advance the application of metabolomics (as a part of NAMs) in chemical safety, primarily as an approach for high throughput screening and prioritisation.

An Extensive Metabolomics Workflow to Discover Cardiotoxin-Induced Molecular Perturbations in Microtissues

Discovering modes of action and predictive biomarkers of drug-induced structural cardiotoxicity offers the potential to improve cardiac safety assessment of lead compounds and enhance preclinical to clinical translation during drug development. Cardiac microtissues are a promising, physiologically relevant, in vitro model, each composed of ca. 500 cells. While untargeted metabolomics is capable of generating hypotheses on toxicological modes of action and discovering metabolic biomarkers, applying this technology to low-biomass microtissues in suspension is experimentally challenging. Thus, we first evaluated a filtration-based approach for harvesting microtissues and assessed the sensitivity and reproducibility of nanoelectrospray direct infusion mass spectrometry (nESI-DIMS) measurements of intracellular extracts, revealing samples consisting of 28 pooled microtissues, harvested by filtration, are suitable for profiling the intracellular metabolome and lipidome. Subsequently, an extensive workflow combining nESI-DIMS untargeted metabolomics and lipidomics of intracellular extracts with ultra-high performance liquid chromatography-mass spectrometry (UHPLC-MS/MS) analysis of spent culture medium, to profile the metabolic footprint and quantify drug exposure concentrations, was implemented. Using the synthetic drug and model cardiotoxin sunitinib, time-resolved metabolic and lipid perturbations in cardiac microtissues were investigated, providing valuable data for generating hypotheses on toxicological modes of action and identifying putative biomarkers such as disruption of purine metabolism and perturbation of polyunsaturated fatty acid levels.

ISA API: An open platform for interoperable life science experimental metadata

BACKGROUND

The Investigation/Study/Assay (ISA) Metadata Framework is an established and widely used set of open source community specifications and software tools for enabling discovery, exchange, and publication of metadata from experiments in the life sciences. The original ISA software suite provided a set of user-facing Java tools for creating and manipulating the information structured in ISA-Tab-a now widely used tabular format. To make the ISA framework more accessible to machines and enable programmatic manipulation of experiment metadata, the JSON serialization ISA-JSON was developed.

RESULTS

In this work, we present the ISA API, a Python library for the creation, editing, parsing, and validating of ISA-Tab and ISA-JSON formats by using a common data model engineered as Python object classes. We describe the ISA API feature set, early adopters, and its growing user community.

CONCLUSIONS

The ISA API provides users with rich programmatic metadata-handling functionality to support automation, a common interface, and an interoperable medium between the 2 ISA formats, as well as with other life science data formats required for depositing data in public databases.

Progress towards an OECD reporting framework for transcriptomics and metabolomics in regulatory toxicology

Omics methodologies are widely used in toxicological research to understand modes and mechanisms of toxicity. Increasingly, these methodologies are being applied to questions of regulatory interest such as molecular point-of-departure derivation and chemical grouping/read-across. Despite its value, widespread regulatory acceptance of omics data has not yet occurred. Barriers to the routine application of omics data in regulatory decision making have been: 1) lack of transparency for data processing methods used to convert raw data into an interpretable list of observations; and 2) lack of standardization in reporting to ensure that omics data, associated metadata and the methodologies used to generate results are available for review by stakeholders, including regulators. Thus, in 2017, the Organisation for Economic Co-operation and Development (OECD) Extended Advisory Group on Molecular Screening and Toxicogenomics (EAGMST) launched a project to develop guidance for the reporting of omics data aimed at fostering further regulatory use. Here, we report on the ongoing development of the first formal reporting framework describing the processing and analysis of both transcriptomic and metabolomic data for regulatory toxicology. We introduce the modular structure, content, harmonization and strategy for trialling this reporting framework prior to its publication by the OECD.

Acoustic Mist Ionization Mass Spectrometry for Ultrahigh-Throughput Metabolomics Screening

Incorporating safety data early in the drug discovery pipeline is key to reducing costly lead candidate failures. For a single drug development project, we estimate that several thousand samples per day require screening (<10 s per acquisition). While chromatography-based metabolomics has proven value at predicting toxicity from metabolic biomarker profiles, it lacks sufficiently high sample throughput. Acoustic mist ionization mass spectrometry (AMI-MS) is an atmospheric pressure ionization approach that can measure metabolites directly from 384-well plates with unparalleled speed. We sought to implement a signal processing and data analysis workflow to produce high-quality AMI-MS metabolomics data and to demonstrate its application to drug safety screening. An existing direct infusion mass spectrometry workflow was adapted, extended, optimized, and tested, utilizing three AMI-MS data sets acquired from technical and biological replicates of metabolite standards and HepG2 cell lysates and a toxicity study. Driven by criteria to minimize variance and maximize feature counts, an algorithm to extract the pulsed scan data was designed; parameters for signal-to-noise-ratio, replicate filter, sample filter, missing value filter, and RSD filter were all optimized; normalization and batch correction strategies were adapted; and cell phenotype filtering was implemented to exclude high cytotoxicity samples. The workflow was demonstrated using a highly replicated HepG2 toxicity data set, comprising 2772 samples from exposures to 16 drugs across 9 concentrations and generated in under 5 h, revealing metabolic phenotypes and individual metabolite changes that characterize specific modes of action. This AMI-MS workflow opens the door to ultrahigh-throughput metabolomics screening, increasing the rate of sample analysis by approximately 2 orders of magnitude.

Struct: an R/bioconductor-based framework for standardised metabolomics data analysis and beyond

Summary

Implementing and combining methods from a diverse range of R/Bioconductor packages into ‘omics’ data analysis workflows represents a significant challenge in terms of standardization, readability and reproducibility. Here, we present an R/Bioconductor package, named struct (Statistics in R using Class-based Templates), which defines a suite of class-based templates that allows users to develop and implement highly standardized and readable statistical analysis workflows. Struct integrates with the STATistics Ontology to ensure consistent reporting and maximizes semantic interoperability. We also present a toolbox, named structToolbox, which includes an extensive set of commonly used data analysis methods that have been implemented using struct. This toolbox can be used to build data-analysis workflows for metabolomics and other omics technologies.

Availability and implementation

struct and structToolbox are implemented in R, and are freely available from Bioconductor (http://bioconductor.org/packages/struct and http://bioconductor.org/packages/structToolbox), including documentation and vignettes. Source code is available and maintained at https://github.com/computational-metabolomics.

Characterization of Monophasic Solvent-Based Tissue Extractions for the Detection of Polar Metabolites and Lipids Applying Ultrahigh-Performance Liquid Chromatography-Mass Spectrometry Clinical Metabolic Phenotyping Assays

Metabolic phenotyping of tissues uses metabolomics and lipidomics to measure the relative polar and nonpolar (lipid) metabolite levels in biological samples. This approach aims to understand disease biochemistry and identify biochemical markers of disease. Sample preparation methods must be reproducible, sensitive (high metabolite and lipid yield), and ideally rapid. We evaluated three biphasic methods for polar and nonpolar compound extraction (chloroform/methanol/water, dichloromethane/methanol/water, and methyl tert-butyl ether [MTBE]/methanol/water), a monophasic method for polar compound extraction (acetonitrile/methanol/water), and a monophasic method for nonpolar compound extraction (isopropanol/water). All methods were applied to mammalian heart, kidney, and liver tissues. Polar extracts were analyzed by hydrophilic interaction chromatography (HILIC) ultrahigh-performance liquid chromatography-mass spectrometry (UHPLC-MS) and nonpolar extracts by C₁₈ reversed-phase UHPLC-MS. Method reproducibility and yield were assessed using multiple annotated endogenous compounds (putatively and MS/MS annotated). Monophasic methods had the highest yield and high reproducibility for both polar (positive ion: median relative standard deviation (RSD) < 18%; negative ion: median RSD < 28%) and nonpolar (positive and negative ion: median RSD < 15%) extractions for heart, kidneys, and liver. The polar monophasic method extracted higher levels of lipid than biphasic polar extractions, and these lipids caused minimal detection suppression for other compounds during HILIC UHPLC-MS. The nonpolar monophasic method had similar or greater detection responses of all detected lipid classes compared to biphasic methods (including increased phosphatidylinositol, phosphatidylserine, and cardiolipin responses). Monophasic methods are quicker and simpler than biphasic methods and are therefore most suited for future automation.

Metabolic characterisation of disturbances in the APOC3/triglyceride-rich lipoprotein pathway through sample-based recall by genotype

INTRODUCTION

High plasma triacylglyceride levels are known to be associated with increased risk of atherosclerotic cardiovascular disease. Apolipoprotein C-III (apoC-III) is a key regulator of plasma triacylglyceride levels and is associated with hypertriglyceridemia via a number of pathways. There is consistent evidence for an association of cardiovascular events with blood apoC-III level, with support from human genetic studies of APOC3 variants. As such, apoC-III has been recognised as a potential therapeutic target for patients with severe hypertriglyceridaemia with one of the most promising apoC-III-targeting drugs, volanesorsen, having recently progressed through Phase III trials.

OBJECTIVES

To exploit a rare loss of function variant in APOC3 (rs138326449) to characterise the potential long-term treatment effects of apoC-III targeting interventions on the metabolome.

METHODS

In a recall-by-genotype study, 115 plasma samples were analysed by UHPLC-MS to acquire non-targeted metabolomics data. The study included samples from 57 adolescents and 33 adults. Overall, 12 985 metabolic features were tested for an association with APOC3 genotype.

RESULTS

161 uniquely annotated metabolites were found to be associated with rs138326449(APOC3). The highest proportion of associated metabolites belonged to the acyl-acyl glycerophospholipid and triacylglyceride metabolite classes. In addition to the anticipated (on-target) reduction of metabolites in the triacylglyceride and related classes, carriers of the rare variant exhibited previously unreported increases in levels of a number of metabolites from the acyl-alkyl glycerophospholipid class.

CONCLUSION

Overall, our results suggest that therapies targeting apoC-III may potentially achieve a broad shift in lipid profile that favours better metabolic health.

Assessment of human plasma and urine sample preparation for reproducible and high-throughput UHPLC-MS clinical metabolic phenotyping

In this study we assess multiple sample preparation methods for UHPLC-MS metabolic phenotyping analysis of human urine and plasma. All methods are discussed in terms of metabolite and lipid coverage and reproducibility.

The metaRbolomics Toolbox in Bioconductor and beyond

Metabolomics aims to measure and characterise the complex composition of metabolites in a biological system. Metabolomics studies involve sophisticated analytical techniques such as mass spectrometry and nuclear magnetic resonance spectroscopy, and generate large amounts of high-dimensional and complex experimental data. Open source processing and analysis tools are of major interest in light of innovative, open and reproducible science. The scientific community has developed a wide range of open source software, providing freely available advanced processing and analysis approaches. The programming and statistics environment R has emerged as one of the most popular environments to process and analyse Metabolomics datasets. A major benefit of such an environment is the possibility of connecting different tools into more complex workflows. Combining reusable data processing R scripts with the experimental data thus allows for open, reproducible research. This review provides an extensive overview of existing packages in R for different steps in a typical computational metabolomics workflow, including data processing, biostatistics, metabolite annotation and identification, and biochemical network and pathway analysis. Multifunctional workflows, possible user interfaces and integration into workflow management systems are also reviewed. In total, this review summarises more than two hundred metabolomics specific packages primarily available on CRAN, Bioconductor and GitHub.

Use cases, best practice and reporting standards for metabolomics in regulatory toxicology

Metabolomics is a widely used technology in academic research, yet its application to regulatory science has been limited. The most commonly cited barrier to its translation is lack of performance and reporting standards. The MEtabolomics standaRds Initiative in Toxicology (MERIT) project brings together international experts from multiple sectors to address this need. Here, we identify the most relevant applications for metabolomics in regulatory toxicology and develop best practice guidelines, performance and reporting standards for acquiring and analysing untargeted metabolomics and targeted metabolite data. We recommend that these guidelines are evaluated and implemented for several regulatory use cases.

Multiple metabolic pathways are predictive of ricin intoxication in a rat model

INTRODUCTION

Exposure to ricin can be lethal and treatments that are under development have short windows of opportunity for administration after exposure. It is therefore essential to achieve early detection of ricin exposure to provide the best prognosis for exposed individuals. Ricin toxin can be detected in clinical samples via several antibody-based techniques, but the efficacy of these can be limited due to the rapid processing and cellular uptake of toxin in the body and subsequent low blood ricin concentrations. Other diagnostic tools that perform, in an orthogonal manner, are therefore desirable.

OBJECTIVES

To determine time-dependent metabolic changes in Sprague-Dawley rats following intravenous exposure to ricin.

METHODS

Sprague-Dawley rats were intravenously exposed to ricin and multiple blood samples were collected from each animal for up to 48 h following exposure in two independent studies. Plasma samples were analysed applying HILIC and C₁₈ reversed phase UHPLC-MS assays followed by univariate and multivariate analysis.

RESULTS

In Sprague-Dawley rats we have demonstrated that metabolic changes measured in blood can distinguish between rats exposed intravenously to ricin and controls prior to the onset of behavioral signs of intoxication after 24 h. A total of 37 metabolites were significantly altered following exposure to ricin when compared to controls. The arginine/proline, bile acid and triacylglyceride metabolic pathways were highlighted as being important with two triacylglycerides at 8 h post exposure giving an AUROC score of 0.94. At 16 h and 24 h the AUROC score increased to 0.98 and 1.0 with the number of metabolites in the panel increasing to 5 and 7, respectively.

CONCLUSIONS

These data demonstrate that metabolites may be a useful tool to diagnose and detect ricin exposure, thus increasing the effectiveness of supportive therapy and future ricin-specific medical treatments.

PhenoMeNal: processing and analysis of metabolomics data in the cloud

BACKGROUND

Metabolomics is the comprehensive study of a multitude of small molecules to gain insight into an organism's metabolism. The research field is dynamic and expanding with applications across biomedical, biotechnological, and many other applied biological domains. Its computationally intensive nature has driven requirements for open data formats, data repositories, and data analysis tools. However, the rapid progress has resulted in a mosaic of independent, and sometimes incompatible, analysis methods that are difficult to connect into a useful and complete data analysis solution.

FINDINGS

PhenoMeNal (Phenome and Metabolome aNalysis) is an advanced and complete solution to set up Infrastructure-as-a-Service (IaaS) that brings workflow-oriented, interoperable metabolomics data analysis platforms into the cloud. PhenoMeNal seamlessly integrates a wide array of existing open-source tools that are tested and packaged as Docker containers through the project's continuous integration process and deployed based on a kubernetes orchestration framework. It also provides a number of standardized, automated, and published analysis workflows in the user interfaces Galaxy, Jupyter, Luigi, and Pachyderm.

CONCLUSIONS

PhenoMeNal constitutes a keystone solution in cloud e-infrastructures available for metabolomics. PhenoMeNal is a unique and complete solution for setting up cloud e-infrastructures through easy-to-use web interfaces that can be scaled to any custom public and private cloud environment. By harmonizing and automating software installation and configuration and through ready-to-use scientific workflow user interfaces, PhenoMeNal has succeeded in providing scientists with workflow-driven, reproducible, and shareable metabolomics data analysis platforms that are interfaced through standard data formats, representative datasets, versioned, and have been tested for reproducibility and interoperability. The elastic implementation of PhenoMeNal further allows easy adaptation of the infrastructure to other application areas and 'omics research domains.

Pulmonary toxicity of inhaled nano-sized cerium oxide aerosols in Sprague-Dawley rats

Cerium oxide nanoparticles (CeO₂NPs), used in some diesel fuel additives to improve fuel combustion efficiency and exhaust filter operation, have been detected in ambient air and concerns have been raised about their potential human health impact. The majority of CeO₂NP inhalation studies undertaken to date have used aerosol particles of larger sizes than the evidence suggests are emitted from vehicles using such fuel additives. Hence, the objective of this study was to investigate the effects of inhaled CeO₂NP aerosols of a more environmentally relevant size, utilizing a combination of methods, including untargeted multi-omics to enable the broadest possible survey of molecular responses and synchrotron X-ray spectroscopy to investigate cerium speciation. Male Sprague-Dawley rats were exposed by nose-only inhalation to aerosolized CeO₂NPs (mass concentration 1.8 mg/m³, aerosol count median diameter 40 nm) for 3 h/d for 4 d/week, for 1 or 2 weeks and sacrificed at 3 and 7 d post-exposure. Markers of inflammation changed significantly in a dose- and time-dependent manner, which, combined with results from lung histopathology and gene expression analyses suggest an inflammatory response greater than that seen in studies using micron-sized ceria aerosols. Lipidomics of lung tissue revealed changes to minor lipid species, implying specific rather than general cellular effects. Cerium speciation analysis indicated a change in Ce³⁺/Ce⁴⁺ ratio within lung tissue. Collectively, these results in conjunction with earlier studies emphasize the importance of aerosol particle size on toxicity determination. Furthermore, the limited effect resolution within 7 d suggested the possibility of longer-term effects.

Applying 'omics technologies in chemicals risk assessment: Report of an ECETOC workshop

Prevailing knowledge gaps in linking specific molecular changes to apical outcomes and methodological uncertainties in the generation, storage, processing, and interpretation of 'omics data limit the application of 'omics technologies in regulatory toxicology. Against this background, the European Centre for Ecotoxicology and Toxicology of Chemicals (ECETOC) convened a workshop Applying 'omics technologies in chemicals risk assessment that is reported herein. Ahead of the workshop, multi-expert teams drafted frameworks on best practices for (i) a Good-Laboratory Practice-like context for collecting, storing and curating 'omics data; (ii) the processing of 'omics data; and (iii) weight-of-evidence approaches for integrating 'omics data. The workshop participants confirmed the relevance of these Frameworks to facilitate the regulatory applicability and use of 'omics data, and the workshop discussions provided input for their further elaboration. Additionally, the key objective (iv) to establish approaches to connect 'omics perturbations to phenotypic alterations was addressed. Generally, it was considered promising to strive to link gene expression changes and pathway perturbations to the phenotype by mapping them to specific adverse outcome pathways. While further work is necessary before gene expression changes can be used to establish safe levels of substance exposure, the ECETOC workshop provided important incentives towards achieving this goal.

The future of metabolomics in ELIXIR

Metabolomics, the youngest of the major omics technologies, is supported by an active community of researchers and infrastructure developers across Europe. To coordinate and focus efforts around infrastructure building for metabolomics within Europe, a workshop on the "Future of metabolomics in ELIXIR" was organised at Frankfurt Airport in Germany. This one-day strategic workshop involved representatives of ELIXIR Nodes, members of the PhenoMeNal consortium developing an e-infrastructure that supports workflow-based metabolomics analysis pipelines, and experts from the international metabolomics community. The workshop established metabolite identification as the critical area, where a maximal impact of computational metabolomics and data management on other fields could be achieved. In particular, the existing four ELIXIR Use Cases, where the metabolomics community - both industry and academia - would benefit most, and which could be exhaustively mapped onto the current five ELIXIR Platforms were discussed. This opinion article is a call for support for a new ELIXIR metabolomics Use Case, which aligns with and complements the existing and planned ELIXIR Platforms and Use Cases.

mzML2ISA & nmrML2ISA: generating enriched ISA-Tab metadata files from metabolomics XML data

SUMMARY

Submission to the MetaboLights repository for metabolomics data currently places the burden of reporting instrument and acquisition parameters in ISA-Tab format on users, who have to do it manually, a process that is time consuming and prone to user input error. Since the large majority of these parameters are embedded in instrument raw data files, an opportunity exists to capture this metadata more accurately. Here we report a set of Python packages that can automatically generate ISA-Tab metadata file stubs from raw XML metabolomics data files. The parsing packages are separated into mzML2ISA (encompassing mzML and imzML formats) and nmrML2ISA (nmrML format only). Overall, the use of mzML2ISA & nmrML2ISA reduces the time needed to capture metadata substantially (capturing 90% of metadata on assay and sample levels), is much less prone to user input errors, improves compliance with minimum information reporting guidelines and facilitates more finely grained data exploration and querying of datasets.

AVAILABILITY AND IMPLEMENTATION

mzML2ISA & nmrML2ISA are available under version 3 of the GNU General Public Licence at https://github.com/ISA-tools. Documentation is available from http://2isa.readthedocs.io/en/latest/.

CONTACT

reza.salek@ebi.ac.uk or isatools@googlegroups.com.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Protein Corona Modulates Uptake and Toxicity of Nanoceria via Clathrin-Mediated Endocytosis

Particles present in diesel exhaust have been proposed as a significant contributor to the development of acute and chronic lung diseases, including respiratory infection and allergic asthma. Nanoceria (CeO2 nanoparticles) are used to increase fuel efficiency in internal combustion engines, are present in exhaust fumes, and could affect cells of the airway. Components from the environment such as biologically derived proteins, carbohydrates, and lipids can form a dynamic layer, commonly referred to as the "protein corona" which alters cellular nanoparticle interactions and internalization. Using confocal reflectance microscopy, we quantified nanoceria uptake by lung-derived cells in the presence and absence of a serum-derived protein corona. Employing mass spectrometry, we identified components of the protein corona, and demonstrated that the interaction between transferrin in the protein corona and the transferrin receptor is involved in mediating the cellular entry of nanoceria via clathrin-mediated endocytosis. Furthermore, under these conditions nanoceria does not affect cell growth, viability, or metabolism, even at high concentration. Alternatively, despite the antioxidant capacity of nanoceria, in serum-free conditions these nanoparticles induce plasma membrane disruption and cause changes in cellular metabolism. Thus, our results identify a specific receptor-mediated mechanism for nanoceria entry, and provide significant insight into the potential for nanoparticle-dependent toxicity.

Galaxy-M: a Galaxy workflow for processing and analyzing direct infusion and liquid chromatography mass spectrometry-based metabolomics data

BACKGROUND

Metabolomics is increasingly recognized as an invaluable tool in the biological, medical and environmental sciences yet lags behind the methodological maturity of other omics fields. To achieve its full potential, including the integration of multiple omics modalities, the accessibility, standardization and reproducibility of computational metabolomics tools must be improved significantly.

RESULTS

Here we present our end-to-end mass spectrometry metabolomics workflow in the widely used platform, Galaxy. Named Galaxy-M, our workflow has been developed for both direct infusion mass spectrometry (DIMS) and liquid chromatography mass spectrometry (LC-MS) metabolomics. The range of tools presented spans from processing of raw data, e.g. peak picking and alignment, through data cleansing, e.g. missing value imputation, to preparation for statistical analysis, e.g. normalization and scaling, and principal components analysis (PCA) with associated statistical evaluation. We demonstrate the ease of using these Galaxy workflows via the analysis of DIMS and LC-MS datasets, and provide PCA scores and associated statistics to help other users to ensure that they can accurately repeat the processing and analysis of these two datasets. Galaxy and data are all provided pre-installed in a virtual machine (VM) that can be downloaded from the GigaDB repository. Additionally, source code, executables and installation instructions are available from GitHub.

CONCLUSIONS

The Galaxy platform has enabled us to produce an easily accessible and reproducible computational metabolomics workflow. More tools could be added by the community to expand its functionality. We recommend that Galaxy-M workflow files are included within the supplementary information of publications, enabling metabolomics studies to achieve greater reproducibility.

Non-targeted UHPLC-MS metabolomic data processing methods: a comparative investigation of normalisation, missing value imputation, transformation and scaling

INTRODUCTION

The generic metabolomics data processing workflow is constructed with a serial set of processes including peak picking, quality assurance, normalisation, missing value imputation, transformation and scaling. The combination of these processes should present the experimental data in an appropriate structure so to identify the biological changes in a valid and robust manner.

OBJECTIVES

Currently, different researchers apply different data processing methods and no assessment of the permutations applied to UHPLC-MS datasets has been published. Here we wish to define the most appropriate data processing workflow.

METHODS

We assess the influence of normalisation, missing value imputation, transformation and scaling methods on univariate and multivariate analysis of UHPLC-MS datasets acquired for different mammalian samples.

RESULTS

Our studies have shown that once data are filtered, missing values are not correlated with m/z, retention time or response. Following an exhaustive evaluation, we recommend PQN normalisation with no missing value imputation and no transformation or scaling for univariate analysis. For PCA we recommend applying PQN normalisation with Random Forest missing value imputation, glog transformation and no scaling method. For PLS-DA we recommend PQN normalisation, KNN as the missing value imputation method, generalised logarithm transformation and no scaling. These recommendations are based on searching for the biologically important metabolite features independent of their measured abundance.

CONCLUSION

The appropriate choice of normalisation, missing value imputation, transformation and scaling methods differs depending on the data analysis method and the choice of method is essential to maximise the biological derivations from UHPLC-MS datasets.

The effect of perhexiline on myocardial protection during coronary artery surgery: a two-centre, randomized, double-blind, placebo-controlled trial

OBJECTIVES

Perhexiline is thought to modulate metabolism by inhibiting mitochondrial carnitine palmitoyltransferase-1, reducing fatty acid uptake and increasing carbohydrate utilization. This study assessed whether preoperative perhexiline improves markers of myocardial protection in patients undergoing coronary artery bypass graft surgery and analysed its effect on the myocardial metabolome.

METHODS

In a prospective, randomized, double-blind, placebo-controlled trial, patients at two centres were randomized to receive either oral perhexiline or placebo for at least 5 days prior to surgery. The primary outcome was a low cardiac output episode in the first 6 h. All pre-specified analyses were conducted according to the intention-to-treat principle with a statistical power of 90% to detect a relative risk of 0.5 and a conventional one-sided α-value of 0.025. A subset of pre-ischaemic left ventricular biopsies was analysed using mass spectrometry-based metabolomics.

RESULTS

Over a 3-year period, 286 patients were randomized, received the intervention and were included in the analysis. The incidence rate of a low cardiac output episode in the perhexiline arm was 36.7% (51/139) vs 34.7% (51/147) in the control arm [odds ratio (OR) 0.92, 95% confidence interval (CI) 0.56–1.50, P = 0.74]. Perhexiline was associated with a reduction in the cardiac index at 6 h [difference in means 0.19, 95% CI 0.07–0.31, P = 0.001] and an increase in inotropic support in the first 12 h (OR 0.55, 95% CI 0.34–0.89, P = 0.015). There were no significant differences in myocardial injury with troponin-T or electrocardiogram, reoperation, renal dysfunction or length of stay. No difference in the preischaemic left ventricular metabolism was identified between groups on metabolomics analysis.

CONCLUSIONS

Preoperative perhexiline does not improve myocardial protection in patients undergoing coronary surgery and in fact reduced perioperative cardiac output, increasing the need for inotropic support. Perhexiline has no significant effect on the mass spectrometry-visible polar myocardial metabolome in vivo in humans, supporting the suggestion that it acts via a pathway that is independent of myocardial carnitine palmitoyltransferase inhibition and may explain the lack of clinical benefit observed following surgery.

ClinicalTrials.Gov ID

NCT00845364.

HAMMER: automated operation of mass frontier to construct in silico mass spectral fragmentation libraries

SUMMARY

Experimental MS(n) mass spectral libraries currently do not adequately cover chemical space. This limits the robust annotation of metabolites in metabolomics studies of complex biological samples. In silico fragmentation libraries would improve the identification of compounds from experimental multistage fragmentation data when experimental reference data are unavailable. Here, we present a freely available software package to automatically control Mass Frontier software to construct in silico mass spectral libraries and to perform spectral matching. Based on two case studies, we have demonstrated that high-throughput automation of Mass Frontier allows researchers to generate in silico mass spectral libraries in an automated and high-throughput fashion with little or no human intervention required.

AVAILABILITY AND IMPLEMENTATION

Documentation, examples, results and source code are available at http://www.biosciences-labs.bham.ac.uk/viant/hammer/.

A stable-isotope mass spectrometry-based metabolic footprinting approach to analyze exudates from phytoplankton

Phytoplankton exudates play an important role in pelagic ecology and biogeochemical cycles of elements. Exuded compounds fuel the microbial food web and often encompass bioactive secondary metabolites like sex pheromones, allelochemicals, antibiotics, or feeding attractants that mediate biological interactions. Despite this importance, little is known about the bioactive compounds present in phytoplankton exudates. We report a stable-isotope metabolic footprinting method to characterise exudates from aquatic autotrophs. Exudates from (13)C-enriched alga were concentrated by solid phase extraction and analysed by high-resolution Fourier transform ion cyclotron resonance mass spectrometry. We used the harmful algal bloom forming dinoflagellate Alexandrium tamarense to prove the method. An algorithm was developed to automatically pinpoint just those metabolites with highly (13)C-enriched isotope signatures, allowing us to discover algal exudates from the complex seawater background. The stable-isotope pattern (SIP) of the detected metabolites then allowed for more accurate assignment to an empirical formula, a critical first step in their identification. This automated workflow provides an effective way to explore the chemical nature of the solutes exuded from phytoplankton cells and will facilitate the discovery of novel dissolved bioactive compounds.

Characterization of isotopic abundance measurements in high resolution FT-ICR and Orbitrap mass spectra for improved confidence of metabolite identification

Currently there is limited information available on the accuracy and precision of relative isotopic abundance (RIA) measurements using high-resolution direct-infusion mass spectrometry (HR DIMS), and it is unclear if this information can benefit automated peak annotation in metabolomics. Here we characterize the accuracy of RIA measurements on the Thermo LTQ FT Ultra (resolution of 100,000-750,000) and LTQ Orbitrap (R = 100,000) mass spectrometers. This first involved reoptimizing the SIM-stitching method (Southam, A. D. Anal. Chem. 2007, 79, 4595-4602) for the LTQ FT Ultra, which achieved a ca. 3-fold sensitivity increase compared to the original method while maintaining a root-mean-squared mass error of 0.16 ppm. Using this method, we show the quality of RIA measurements is highly dependent on signal-to-noise ratio (SNR), with RIA accuracy increasing with higher SNR. Furthermore, a negative offset between the theoretical and empirically calculated numbers of carbon atoms was observed for both mass spectrometers. Increasing the resolution of the LTQ FT Ultra lowered both the sensitivity and the quality of RIA measurements. Overall, although the errors in the empirically calculated number of carbons can be large (e.g., 10 carbons), we demonstrate that RIA measurements do improve automated peak annotation, increasing the number of single empirical formula assignments by >3-fold compared to using accurate mass alone.

Discriminating between different acute chemical toxicities via changes in the daphnid metabolome

Currently, there is widespread interest in exploiting "omics" approaches to screen the toxicity of chemicals, potentially enabling their rapid categorization into classes of defined mode of action (MOA). Direct infusion Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS) metabolomics provides a sensitive and nontargeted analysis of potentially a thousand endogenous metabolites. Our previous work has shown that mass spectra can be recorded from whole-organism homogenate or hemolymph of single adult Daphnia magna. Here we develop multivariate models and discover perturbations to specific metabolic pathways that can discriminate between the acute toxicities of four chemicals to D. magna using FT-ICR MS metabolomics. We focus on model toxicants (cadmium, fenvalerate, dinitrophenol, and propranolol) with different MOAs. First, we confirmed that a toxicant-induced metabolic effect could be determined for each chemical in both the hemolymph and the whole-organism metabolome, with between 9 and 660 mass spectral peaks changing intensities significantly, dependent upon toxicant and sample type. Subsequently, supervised multivariate models were built that discriminated significantly all four acute metabolic toxicities, yielding mean classification error rates (across all classes) of 3.9 and 6.9% for whole-organism homogenates and hemolymph, respectively. Following extensive peak annotation, we discovered toxicant-specific perturbations to putatively identified metabolic pathways, including propranolol-induced disruption of fatty acid metabolism and eicosanoid biosynthesis and fenvalerate-induced disruption of amino sugar metabolism. We conclude that the metabolic profiles of whole-daphnid homogenates are more discriminatory for toxicant action than hemolymph. Furthermore, our findings highlight the capability of metabolomics to discover early-event metabolic responses that can discriminate between the acute toxicities of chemicals.