Data standards can boost metabolomics research, and if there is a will, there is a way.

Thousands of articles using metabolomics approaches are published every year. With the increasing amounts of data being produced, mere description of investigations as text in manuscripts is not sufficient to enable re-use anymore: the underlying data needs to be published together with the findings in the literature to maximise the benefit from public and private expenditure and to take advantage of an enormous opportunity to improve scientific reproducibility in metabolomics and cognate disciplines. Reporting recommendations in metabolomics started to emerge about a decade ago and were mostly concerned with inventories of the information that had to be reported in the literature for consistency. In recent years, metabolomics data standards have developed extensively, to include the primary research data, derived results and the experimental description and importantly the metadata in a machine-readable way. This includes vendor independent data standards such as mzML for mass spectrometry and nmrML for NMR raw data that have both enabled the development of advanced data processing algorithms by the scientific community. Standards such as ISA-Tab cover essential metadata, including the experimental design, the applied protocols, association between samples, data files and the experimental factors for further statistical analysis. Altogether, they pave the way for both reproducible research and data reuse, including meta-analyses. Further incentives to prepare standards compliant data sets include new opportunities to publish data sets, but also require a little "arm twisting" in the author guidelines of scientific journals to submit the data sets to public repositories such as the NIH Metabolomics Workbench or MetaboLights at EMBL-EBI. In the present article, we look at standards for data sharing, investigate their impact in metabolomics and give suggestions to improve their adoption.

1H NMR-based profiling reveals differential immune-metabolic networks during influenza virus infection in obese mice.

Obese individuals are at greater risk for death from influenza virus infection. Paralleling human evidence, obese mice are also more susceptible to influenza infection mortality. However, the underlying mechanisms driving greater influenza severity in the obese remain unclear. Metabolic profiling has been utilized in infectious disease models to enhance prognostic or diagnostic methods, and to gain insight into disease pathogenesis by providing a more global picture of dynamic infection responses. Herein, metabolic profiling was used to develop a deeper understanding of the complex processes contributing to impaired influenza protection in obese mice and to facilitate generation of new explanatory hypotheses. Diet-induced obese and lean mice were infected with influenza A/Puerto Rico/8/34. 1H nuclear magnetic resonance-based metabolic profiling of urine, feces, lung, liver, mesenteric white adipose tissue, bronchoalveolar lavage fluid and serum revealed distinct metabolic signatures in infected obese mice, including perturbations in nucleotide, vitamin, ketone body, amino acid, carbohydrate, choline and lipid metabolic pathways. Further, metabolic data was integrated with immune analyses to obtain a more comprehensive understanding of potential immune-metabolic interactions. Of interest, uncovered metabolic signatures in urine and feces allowed for discrimination of infection status in both lean and obese mice at an early influenza time point, which holds prognostic and diagnostic implications for this methodology. These results confirm that obesity causes distinct metabolic perturbations during influenza infection and provide a basis for generation of new hypotheses and use of this methodology in detection of putative biomarkers and metabolic patterns to predict influenza infection outcome.

Modelling short time series in metabolomics: a functional data analysis approach.

Metabolomics is the study of the complement of small molecule metabolites in cells, biofluids and tissues. Many metabolomic experiments are designed to compare changes observed over time under two or more experimental conditions (e.g. a control and drug-treated group), thus producing time course data. Models from traditional time series analysis are often unsuitable because, by design, only very few time points are available and there are a high number of missing values. We propose a functional data analysis approach for modelling short time series arising in metabolomic studies which overcomes these obstacles. Our model assumes that each observed time series is a smooth random curve, and we propose a statistical approach for inferring this curve from repeated measurements taken on the experimental units. A test statistic for detecting differences between temporal profiles associated with two experimental conditions is then presented. The methodology has been applied to NMR spectroscopy data collected in a pre-clinical toxicology study.

Time-resolved metabolic footprinting for nonlinear modeling of bacterial substrate utilization.

Untargeted profiling of small-molecule metabolites from microbial culture supernatants (metabolic footprinting) has great potential as a phenotyping tool. We used time-resolved metabolic footprinting to compare one Escherichia coli and three Pseudomonas aeruginosa strains growing on complex media and show that considering metabolite changes over the whole course of growth provides much more information than analyses based on data from a single time point. Most strikingly, there was pronounced selectivity in metabolite uptake, even when the bacteria were growing apparently exponentially, with certain groups of metabolites not taken up until others had been entirely depleted from the medium. In addition, metabolite excretion showed some complex patterns. Fitting nonlinear equations (four-parameter sigmoids) to individual metabolite data allowed us to model these changes for metabolite uptake and visualize them by back-projecting the curve-fit parameters onto the original growth curves. These "uptake window" plots clearly demonstrated strain differences, with the uptake of some compounds being reversed in order between different strains. Comparison of an undefined rich medium with a defined complex medium designed to mimic cystic fibrosis sputum showed many differences, both qualitative and quantitative, with a greater proportion of excreted to utilized metabolites in the defined medium. Extending the strain comparison to a more closely related set of isolates showed that it was possible to discriminate two species of the Burkholderia cepacia complex based on uptake dynamics alone. We believe time-resolved metabolic footprinting could be a valuable tool for many questions in bacteriology, including isolate comparisons, phenotyping deletion mutants, and as a functional complement to taxonomic classifications.

Human metabolic phenotype diversity and its association with diet and blood pressure.

Metabolic phenotypes are the products of interactions among a variety of factors-dietary, other lifestyle/environmental, gut microbial and genetic. We use a large-scale exploratory analytical approach to investigate metabolic phenotype variation across and within four human populations, based on 1H NMR spectroscopy. Metabolites discriminating across populations are then linked to data for individuals on blood pressure, a major risk factor for coronary heart disease and stroke (leading causes of mortality worldwide). We analyse spectra from two 24-hour urine specimens for each of 4,630 participants from the INTERMAP epidemiological study, involving 17 population samples aged 40-59 in China, Japan, UK and USA. We show that urinary metabolite excretion patterns for East Asian and western population samples, with contrasting diets, diet-related major risk factors, and coronary heart disease/stroke rates, are significantly differentiated (P < 10(-16)), as are Chinese/Japanese metabolic phenotypes, and subgroups with differences in dietary vegetable/animal protein and blood pressure. Among discriminatory metabolites, we quantify four and show association (P < 0.05 to P < 0.0001) of mean 24-hour urinary formate excretion with blood pressure in multiple regression analyses for individuals. Mean 24-hour urinary excretion of alanine (direct) and hippurate (inverse), reflecting diet and gut microbial activities, are also associated with blood pressure of individuals. Metabolic phenotyping applied to high-quality epidemiological data offers the potential to develop an area of aetiopathogenetic knowledge involving discovery of novel biomarkers related to cardiovascular disease risk.

Comparative metabonomics of differential hydrazine toxicity in the rat and mouse.

Interspecies variation between rats and mice has been studied for hydrazine toxicity using a novel metabonomics approach. Hydrazine hydrochloride was administered to male Sprague-Dawley rats (30 mg/kg, n = 10 and 90 mg/kg, n = 10) and male B6C3F mice (100 mg/kg, n = 8 and 250 mg/kg, n = 8) by oral gavage. In each species, the high dose was selected to produce the major histopathologic effect, hepatocellular lipid accumulation. Urine samples were collected at sequential time points up to 168 h post dose and analyzed by 1H NMR spectroscopy. The metabolites of hydrazine, namely diacetyl hydrazine and 1,4,5,6-tetrahydro-6-oxo-3-pyridazine carboxylic acid (THOPC), were detected in both the rat and mouse urine samples. Monoacetyl hydrazine was detected only in urine samples from the rat and its absence in the urine of the mouse was attributed to a higher activity of N-acetyl transferases in the mouse compared with the rat. Differential metabolic effects observed between the two species included elevated urinary beta-alanine, 3-D-hydroxybutyrate, citrulline, N-acetylcitrulline, and reduced trimethylamine-N-oxide excretion unique to the rat. Metabolic principal component (PC) trajectories highlighted the greater degree of toxic response in the rat. A data scaling method, scaled to maximum aligned and reduced trajectories (SMART) analysis, was used to remove the differences between the metabolic starting positions of the rat and mouse and varying magnitudes of effect, to facilitate comparison of the response geometries between the rat and mouse. Mice followed "biphasic" open PC trajectories, with incomplete recovery 7 days after dosing, whereas rats followed closed "hairpin" time profiles, indicating functional reversibility. The greater magnitude of metabolic effects observed in the rat was supported by the more pronounced effect on liver pathology in the rat when compared with the mouse.