Mari van Reenen

Contribution towards a metabolite profile of the detoxification of benzoic acid through glycine conjugation: an intervention study

Benzoic acid is widely used as a preservative in food products and is detoxified in humans through glycine conjugation. Different viewpoints prevail on the physiological significance of the glycine conjugation reaction and concerns have been raised on potential public health consequences following uncontrolled benzoic acid ingestion. We performed a metabolomics study which used commercial benzoic acid containing flavored water as vehicle for designed interventions, and report here on the controlled consumption of the benzoic acid by 21 cases across 6 time points for a total of 126 time points. Metabolomics data from urinary samples analyzed by nuclear magnetic resonance spectroscopy were generated in a time-dependent cross-over study. We used ANOVA-simultaneous component analysis (ASCA), repeated measures analysis of variance (RM-ANOVA) and unfolded principal component analysis (unfolded PCA) to supplement conventional statistical methods to uncover fully the metabolic perturbations due to the xenobiotic intervention, encapsulated in the metabolomics tensor (three-dimensional matrices having cases, spectral areas and time as axes). Identification of the biologically important metabolites by the novel combination of statistical methods proved the power of this approach for metabolomics studies having complex data structures in general. The study disclosed a high degree of inter-individual variation in detoxification of the xenobiotic and revealed metabolic information, indicating that detoxification of benzoic acid through glycine conjugation to hippuric acid does not indicate glycine depletion, but is supplemented by ample glycine regeneration. The observations lend support to the view of maintenance of glycine homeostasis during detoxification. The study indicates also that time-dependent metabolomics investigations, using designed interventions, provide a way of interpreting the variation induced by the different factors of a designed experiment–an approach with potential to advance significantly our understanding of normal and pathophysiological perturbations of endogenous or exogenous origin.

Variable selection for binary classification using error rate p-values applied to metabolomics data

BackgroundMetabolomics datasets are often high-dimensional though only a limited number of variables are expected to be informative given a specific research question. The important task of selecting informative variables can therefore become complex. In this paper we look at discriminating between two groups. Two tasks need to be performed: (i) finding variables which differ between the two groups; and (ii) determining how the selected variables can be used to classify new subjects. We introduce an approach using minimum classification error rates as test statistics to find discriminatory and therefore informative variables. The thresholds resulting in the minimum error rates can be used to classify new subjects. This approach transforms error rates into p-values and is referred to as ERp.ResultsWe show that non-parametric hypothesis testing, based on minimum classification error rates as test statistics, can find statistically significantly shifted variables. The discriminatory ability of variables becomes more apparent when error rates are evaluated based on their corresponding p-values, as relatively high error rates can still be statistically significant. ERp can handle unequal and small group sizes, as well as account for the cost of misclassification. ERp retains (if known) or reveals (if unknown) the shift direction, aiding in biological interpretation. The threshold resulting in the minimum error rate can immediately be used to classify new subjects.We use NMR generated metabolomics data to illustrate how ERp is able to discriminate subjects diagnosed with Mycobacterium tuberculosis infected meningitis from a control group. The list of discriminatory variables produced by ERp contains all biologically relevant variables with appropriate shift directions discussed in the original paper from which this data is taken.ConclusionsERp performs variable selection and classification, is non-parametric and aids biological interpretation while handling unequal group sizes and misclassification costs. All this is achieved by a single approach which is easy to perform and interpret. ERp has the potential to address many other characteristics of metabolomics data. Future research aims to extend ERp to account for a large proportion of observations below the detection limit, as well as expand on interactions between variables.

Metabolic risks at birth of neonates exposed in utero to HIV-antiretroviral therapy relative to unexposed neonates: an NMR metabolomics study of cord blood

IntroductionAntiretroviral therapy (ART) for HIV-infected pregnant women is highly effective in preventing mother-to-child transmission (PMTCT) of the virus, but deleterious metabolic and mitochondrial observations in infants born to HIV-infected women treated with ART during pregnancy are periodically reported.ObjectivesThis study addresses the concern of HIV-ART-induced metabolic perturbations through a metabolomics study of cord blood collected during transitional neonatal hypoglycaemia following birth from newborns either exposed or unexposed to fetal HIV-ART.MethodsProton magnetic resonance spectra from cord blood of 11 in utero HIV-ART-exposed and 14 unexposed newborns, as well as serum from 8 control infants, generated 114 spectral bins which were used to identify significant metabolites by means of univariate and multivariate statistical analyses.ResultsThe metabolite profiles differed significantly between that from the unexposed newborns and that from infants—interpreted to characterize the state of transitional neonatal hypoglycaemia (low glucose and high lactic acid and ketone bodies). Quantitative analysis of potential ATP generation showed no meaningful difference in the global metabolite profiles of HIV-ART-exposed and unexposed neonates, but Volcano plot analysis, affirmed by odds ratios, indicated that exposure to HIV-ART affected the plasma 3-hydroxybutyric acid and hypoxanthine concentrations.ConclusionsThe metabolite profile for transitional neonatal hypoglycaemia indicated that HIV-ART did not compromise the exposed neonates to the energy stress of allostasis experienced at birth. Increased hypoxanthine and 3-hydroxybutyric acid indicates metabolic stress at birth in some of the newborns exposed to HIV-ART and raises a concern about unrecognized prolonged allostasis with potential neurological consequences for these infants.

Metabolomics variable selection and classification in the presence of observations below the detection limit using an extension of ERp

BackgroundERp is a variable selection and classification method for metabolomics data. ERp uses minimized classification error rates, based on data from a control and experimental group, to test the null hypothesis of no difference between the distributions of variables over the two groups. If the associated p-values are significant they indicate discriminatory variables (i.e. informative metabolites). The p-values are calculated assuming a common continuous strictly increasing cumulative distribution under the null hypothesis. This assumption is violated when zero-valued observations can occur with positive probability, a characteristic of GC-MS metabolomics data, disqualifying ERp in this context. This paper extends ERp to address two sources of zero-valued observations: (i) zeros reflecting the complete absence of a metabolite from a sample (true zeros); and (ii) zeros reflecting a measurement below the detection limit. This is achieved by allowing the null cumulative distribution function to take the form of a mixture between a jump at zero and a continuous strictly increasing function. The extended ERp approach is referred to as XERp.ResultsXERp is no longer non-parametric, but its null distributions depend only on one parameter, the true proportion of zeros. Under the null hypothesis this parameter can be estimated by the proportion of zeros in the available data. XERp is shown to perform well with regard to bias and power. To demonstrate the utility of XERp, it is applied to GC-MS data from a metabolomics study on tuberculosis meningitis in infants and children. We find that XERp is able to provide an informative shortlist of discriminatory variables, while attaining satisfactory classification accuracy for new subjects in a leave-one-out cross-validation context.ConclusionXERp takes into account the distributional structure of data with a probability mass at zero without requiring any knowledge of the detection limit of the metabolomics platform. XERp is able to identify variables that discriminate between two groups by simultaneously extracting information from the difference in the proportion of zeros and shifts in the distributions of the non-zero observations. XERp uses simple rules to classify new subjects and a weight pair to adjust for unequal sample sizes or sensitivity and specificity requirements.

Metabolic risks of neonates at birth following in utero exposure to HIV-ART: the amino acid profile of cord blood

IntroductionUntargeted metabolomics of cord blood indicated that antiretroviral therapy to HIV-infected mothers (HIV-ART) did not compromise the exposed neonates with regard to the stress of neonatal hypoglycaemia at birth. However, identified biomarkers reflected stress in their energy metabolism, raising concern over developmental risks in some newborns exposed to ART.ObjectivesThis study addresses the concern over HIV-ART-induced metabolic perturbations by expanding the metabolomics study to the amino acid profiles in cord blood collected at birth from newborns either exposed or unexposed to HIV-ART in utero.MethodsAmino acid profiles derived from liquid chromatographic triple quadruple spectra of cord blood from neonates exposed and unexposed to HIV-ART (cohort 1) were investigated using a metabolomics approach. Amino acid data, generated by ultra performance liquid chromatography–tandem mass spectrometry from similar cases (cohort 2), were included for comparison.ResultsMultivariate and supporting statistics indicated differentiation between the exposed and unexposed neonates in both cohorts, caused by a general decrease or downregulation of amino acid concentrations in the cord blood samples from the exposed cases. Specifically, significant upregulation of aspartic acid in both cohorts and downregulation of arginine, and of threonine, tryptophan and lysine in cohorts 1 and 2, respectively, were observed.ConclusionsThe benefits of ART for HIV-infected pregnant women are well established. However, the amino acid profile of cord blood, obtained from the two independent cohorts, adds to observed metabolic risks of in utero HIV-ART-exposed newborns. These risks could potentially have adverse consequences for the future health of some exposed infants.

Trends in malaria case management following changes in the treatment policy to artemisinin combination therapy at the Mbakong Health Centre, Cameroon 2006–2012: a retrospective study

National malaria treatment policies are devised to guide health professionals and to facilitate diagnosis and case management. Following the recommendations of the WHO, Cameroon changed its malaria treatment policy from monotherapy to artemisinin-based combination therapy (ACT) as the first-line treatment for uncomplicated malaria. We report an investigation into trends of case management following this change in policy. Data was collected retrospectively, through consultation and perusal of laboratory and prescription registers of the Mbakong Health Centre. Analysis of data was done using SPSS and SAS Statistics. Data presented herein demonstrate that from 2006 to 2012, a total of 2484 (58.7%) of the total prescriptions included an anti-malarial, 1989 (47.0%) included an antibiotic and 1935(45.7%) included an antipyretic. The anti-malarials prescribed were Anti-malaria combination therapy (ACT) - 1216 (47.6%), quinine 1044 (40.8%) or SP 296 (11.6%). Of the 1216 patients prescribed an ACT, 441(36.3%) had a positive malaria parasite confirmation, 746 (61.3%) were negative for plasmodium. Overall, 29 patients (2.4%) were treated either with an ACT without any test performed. Quinine intake was recorded in 566 (54.2%) patients positive for plasmodium. ACT prescription increased from 23% in 2007 to between 44 and 45% in 2008-2009. During this period there was a corresponding drop in the prescription of quinine from 38% in 2007 to 13% in 2009 (r=-0.43, p>0.05). Sulphadoxine-Pyrimethamine (SP) was restrictively prescribed to women of childbearing age (97.0%) after 2008. Antibiotics prescription dropped from 53.7% to 39.3% from 2010 to 2012. The odds of being prescribed an antibiotic was significantly higher in patients with a malaria negative result compared to malaria positive patients (OR=6.12, CI 4.74-7.91, p<0.00001). Overall, there is an over treatment of malaria, thus departing from the WHO guidelines of appropriate treatment. Although there is an overall increase in the prescription of ACT, less prescription of quinine and a noticeable restrain from prescription of SP to febrile cases, the old practice was still rampant. There is need for healthcare workers to adhere to guidelines in order to enhance the rational use of drugs to achieve appropriate treatment of uncomplicated malaria according to WHO guidelines.