Skip to main content
Advertisement
  • Loading metrics

Discovery of Infection Associated Metabolic Markers in Human African Trypanosomiasis

  • Sabrina D. Lamour,

    Affiliation Section of Biomolecular Medicine, Division of Computational and Systems Medicine, Department of Surgery and Cancer, Imperial College London, London, United Kingdom

  • Maria Gomez-Romero,

    Affiliation Section of Hepatology and Gastroenterology, Department of Medicine, Imperial College London, London, United Kingdom

  • Panagiotis A. Vorkas,

    Affiliation Section of Biomolecular Medicine, Division of Computational and Systems Medicine, Department of Surgery and Cancer, Imperial College London, London, United Kingdom

  • Vincent P. Alibu,

    Affiliation Section of Hepatology and Gastroenterology, Department of Medicine, Imperial College London, London, United Kingdom

  • Jasmina Saric,

    Affiliation Section of Biomolecular Medicine, Division of Computational and Systems Medicine, Department of Surgery and Cancer, Imperial College London, London, United Kingdom

  • Elaine Holmes ,

    ‡ These authors are joint senior authors on this work.

    Affiliation Section of Biomolecular Medicine, Division of Computational and Systems Medicine, Department of Surgery and Cancer, Imperial College London, London, United Kingdom

  • Jeremy M. Sternberg

    j.sternberg@abdn.ac.uk

    ‡ These authors are joint senior authors on this work.

    Affiliation Institute of Biological and Environmental Sciences, University Of Aberdeen, Aberdeen, United Kingdom

Abstract

Human African trypanosomiasis (HAT) remains a major neglected tropical disease in Sub-Saharan Africa. As clinical symptoms are usually non-specific, new diagnostic and prognostic markers are urgently needed to enhance the number of identified cases and optimise treatment. This is particularly important for disease caused by Trypanosoma brucei rhodesiense, where indirect immunodiagnostic approaches have to date been unsuccessful. We have conducted global metabolic profiling of plasma from T.b.rhodesiense HAT patients and endemic controls, using 1H nuclear magnetic resonance (NMR) spectroscopy and ultra-performance liquid chromatography, coupled with mass spectrometry (UPLC-MS) and identified differences in the lipid, amino acid and metabolite profiles. Altogether 16 significantly disease discriminatory metabolite markers were found using NMR, and a further 37 lipid markers via UPLC-MS. These included significantly higher levels of phenylalanine, formate, creatinine, N-acetylated glycoprotein and triglycerides in patients relative to controls. HAT patients also displayed lower concentrations of histidine, sphingomyelins, lysophosphatidylcholines, and several polyunsaturated phosphatidylcholines. While the disease metabolite profile was partially consistent with previous data published in experimental rodent infection, we also found unique lipid and amino acid profile markers highlighting subtle but important differences between the host response to trypanosome infections between animal models and natural human infections. Our results demonstrate the potential of metabolic profiling in the identification of novel diagnostic biomarkers and the elucidation of pathogenetic mechanisms in this disease.

Author Summary

Metabolic profiling of biofluids and tissues in disease and healthy individuals is a powerful approach to discover new markers for diagnosis. We have applied these techniques to the protozoan infection human African trypanosomiasis (HAT), otherwise known as sleeping sickness. The form of HAT endemic in East Africa, caused by Trypanosoma brucei rhodesiense, requires technically demanding direct microscopic diagnosis, as there are no indirect rapid diagnostic tests available. We studied the metabolite profiles in plasma from HAT patients and controls. Clear biochemical differences were discovered between control individuals and patients, including changes in the overall lipid composition and concentration of certain amino acids. These may have been caused by the inflammatory immune response to infection and the uptake of particular molecules by the parasites, although further research will be required for confirmation. We demonstrate that plasma metabolic profiles are characteristic for T. b. rhodesiense infection. While some of these changes are consistent with those observed in an experimental mouse infection model of HAT, many are unique to this clinical study and indicate the necessity of validating experimental animal study data in clinical disease studies. Our results also reveal biochemical changes in patients that will help us understand the development of disease.

Introduction

Human African Trypanosomiasis (HAT) is caused by infection with either of two subspecies of Trypanosoma brucei. Trypanosoma brucei (T. b.) gambiense causes chronic disease (that can last months or years) and is endemic in Western Africa while T. b. rhodiesiense causes a more acute illness in Eastern and Southern Africa, which is typically fatal within less than a year of infection if untreated [1,2]. Diagnosis is a critical challenge for treatment and control of this disease as patients display an array of non-specific inflammatory symptoms, often indistinguishable from other endemic illnesses such as malaria or enteric fever [1,3]. Infection with T. b. gambiense is routinely screened using the Card Agglutination Test for Trypanosomiasis (CATT) and also a new generation of lateral-flow rapid diagnostic tests are being deployed based on the host-response to commonly expressed variant surface glycoproteins (VSG) [4]. However no such serological approach has been successful for T. b. rhodesiense probably due to the greater diversity of VSG expression. New diagnostic techniques are urgently needed for this disease in order to meet the World Health Organization (WHO) target of elimination by 2020 [5,6].

In search for novel markers of disease, we have conducted global untargeted metabolic profiling of plasma from rhodesiense trypanosomiasis patients from Uganda and matched controls. Metabolic profiling (also termed metabonomics/metabolomics) can be used to characterise biochemical patterns associated with specific physiological and/or pathological states [7]. In the field of parasitic infections, this approach has shown capacity for characterising infection-induced metabolic changes in the host, within infected tissues and systemically (as measured by plasma or urine profiles). As yet, parasitic metabolic profiling studies have largely focused on in vitro assays [810] and experimental rodent models [1115]. There have been very few examples of identifying the metabolic signature of a specific parasitic infection in humans [16], largely because human profiles are confounded by strong genetic and environmental variation, and are often superimposed upon a background of other concurrent endemic infections. In mice, infection with T. b. brucei (a subspecies of T. brucei that is not infective to humans) resulted in augmented plasma levels of lactate, acetylglycoproteins and creatine, and reductions in phosphatidylcholine and lipoproteins, as well as elevated pro-inflammatory cytokines [11,17]. Increases in plasma pro-inflammatory cytokines have also been measured in clinical samples [18,19], although metabolic similarities between the mouse model and humans remain to be investigated. In this study we present the first characterisation the metabolic effects of T. brucei rhodesiense infection in humans, and highlight both similarities and differences to results from published mouse model infections with T.brucei brucei.

Plasma samples were analysed by proton nuclear magnetic resonance (1H NMR) spectroscopy and reversed-phase ultra-performance liquid chromatography, coupled to mass spectrometry (RP-UPLC-MS). 1H NMR is a highly reproducible method optimal for analysing complex biological mixtures, such as plasma, with minimal sample preparation [20]. Evidence from a range of experimental models have highlighted marked changes in plasma lipids following T. brucei ssp infection, such as hypertriglyceridaemia [11,2123]. Consequently, lipid profiles were separately characterised via UPLC-MS, to provide complementary information on the different lipid species altered in disease. We hypothesised that the metabolic phenotype of individuals with HAT rhodesiense infection would be distinct from that of uninfected individuals, yielding a range of discriminatory markers that may be of diagnostic value.

Methods

Ethics Statement

The human plasma samples in this study were collected under protocols approved by ethics committees in Uganda (UNCST) and UK (North of Scotland Research Ethics Committee), conforming to the principles of the Declaration of Helsinki. Ethical consent forms and information sheets were designed in English and translated into local languages. Informed consent was given as a signature or a thumb-print (as approved by the UNCST) after verbal explanation. For those aged under 18, affirmative consent as well as the consent of their legal guardian was obtained.

Patient Information and Sample Collection

A total of 46 HAT patients and 21 controls were recruited at Lwala Hospital, Kaberamaido District and Serere Health Centre, Serere District in Eastern Uganda, between November 2008 and March 2010. Controls were healthy individuals who were confirmed as non-infected with either trypanosomes or Plasmodium sp. having undergone parasitological assessment. Study sites, recruitment protocols, treatment regimens, disease progression characteristics and clinical examination methods have been published elsewhere [24]. Patients with intercurrent infections of malaria, filariasis or schistosomiasis were excluded. Plasma and cerebral spinal fluid (CSF) samples were collected from all patients prior to treatment as part of normal diagnostic and staging procedures. Staging was carried out in accordance with WHO criteria [25] defining late stage by the presence of parasites in the lumbar CSF and/or a CSF white blood cell count (WBC) > 5/μl. Aliquots of 1–2 ml of plasma were immediately frozen and then maintained in liquid nitrogen until transfer to the UK. After air-freight (24 hours on dry ice) samples were maintained at -80°C until analysis. A summary of the parasitological and demographic characteristics of the patients is presented in supplementary S1 Table.

1H NMR Acquisition and Processing

All sample preparation was performed in one batch in a randomised order. Prior to acquisition, plasma samples were diluted 1:1 with plasma buffer (0.142 M NaHPO4, 2 mM NaN3, 0.08% (volume/volume i.e. v/v) 3-(Trimethylsilyl)propionic-2,2,3,3-d4 acid sodium salt (TSP) solution, 20% (v/v) D2O; all SIGMA-Aldrich, Germany), and transferred into 5 mm NMR tubes. 1H NMR data were acquired (over 3 days) on a Bruker Avance 600 MHz (14.1 T) NMR Spectrometer with BB probe head and refrigerated SampleJet autosampler (Bruker, Germany), using Topsin (3.1) software (Bruker BioSpin, Germany). Samples were acquired at 300 K (~27°C) using a standard one-dimensional pulse with water suppression program for general metabolite screening [26]; sequence as follows: 2 second (s) relaxation delay (RD)– 90° pulse– 4 μs delay—90° pulse– 100 μs mixing time– 90° pulse– 2.73 s acquisition (ACQ) of free induction decay (FID) signal. Number of dummy scans (ds) was 4, and number of acquired FID’s (ns) was 64. To enhance the visualisation of low molecular weight molecules, a Carr-Purcell Meiboom-Gill (CPMG) spin-echo pulse spectrum was also acquired for each sample [26], with the following pulse sequence: 2 s RD– 90° pulse–(τ– 180°–τ)n-ACQ, where spin echo time τ = 300 μs, number of loops n = 128, total spin echo time 2nτ = 76.8 ms, ds = 4 and ns = 64. All FID’s were multiplied by an exponential function with a 0.3 Hz line broadening prior to Fourier transformation into spectral data (spectral width of 18 kHz collected over 131,072 data points). Spectra were processed using automatic phasing, baseline correction and calibrated to the TSP peak via an in-house algorithm in MATLAB (version R2013a, Mathworks Inc.) and Topspin. The residual water signal was removed prior to automatic spectral alignment [27], and probabilistic quotient normalisation [28], using an in-house MATLAB script. Additionally, artefactual peaks arising from EDTA buffer [29] present in blood collection tubes were removed from all plasma spectra for final analysis. The regions removed were 2.53–2.62, 2.7–2.735, 3.08–3.256, 3.60–3.66 ppm which may have resulted in the loss of some endogenous metabolic information. Spectral assignments of metabolite peaks were identified using Chenomx NMR suite profiler 8.1 software and known assignments from in-house NMR databases and literature [30,31].

Plasma Lipid Profiling UPLC-MS Acquisition and Processing

Due to sample limitations, additional aliquots of only a subset of plasma samples were used for lipid profiling (n = 30, see supplementary S1 Table). As for NMR, sample preparation for UPLC-MS was performed in one day in one batch, in a randomised manner. Samples underwent protein precipitation using isopropanol, as described by Sarafian et al, 2014 [32]. Briefly, 50 μl plasma was mixed 1:3 (v/v) with ice-cold isopropanol (SIGMA-Aldrich), incubated for 10 mins at room temperature, and then overnight at -20°C. Samples were then centrifuged at 20,817 x g (14,000 rpm) at 4°C for 20 mins and 100 μl of the supernatant (subsequently referred to as test samples) were transferred to MS vials (Waters Corp., UK) for analysis and placed in the auto-sampler (kept at 8°C). Additional aliquots from each plasma-isopropanol sample were pooled to generate a quality control (QC) master sample. The QC was injected before starting the run to optimise the injection volume and to condition the column, and then injected every 8 samples for assurance of a stable run.

Lipid profiling was performed on an Acquity UPLC system (Waters Ltd. Elstree, UK) coupled to a Q-TOF Premier mass spectrometer (Waters Ltd., Manchester, UK). Mass spectra for all samples were acquired continuously, firstly in positive and then negative ion electrospray (ESI+ and ESI-) modes, over the course of 3 days. The separation conditions for reversed-phase UPLC have been previously established [33,34]. Gradient elution was performed using a CSH C18 (1.7 μm, 2.1 x 100 mm) column (Waters Corporation, Milford, USA) kept at 55°C. The mobile phases consisted of 0.1% (v/v) formic acid and 10 mM ammonium formate in 60:40 (v/v) acetonitrile (ACN)/HPLC-grade water (mobile phase A) and 0.1% formic acid (v/v) and 10 mM ammonium formate in 90:10 (v/v) IPA/ACN (mobile phase B) at a flow rate of 0.4 ml/min (all solvents purchased from SIGMA-Aldrich except water which was from Fisher, Germany). The sample injection volume was 5 μl and 10 μl for ESI+ and ESI-, respectively.

The MS parameters were set as follows: capillary voltage, 3 kV (ESI+) and 2.5 kV (ESI-); sample cone voltage, 30 V (ESI+) and 25 V (ESI-); source temperature, 120°C; desolvation temperature, 400°C; desolvation gas flow, 800 L/h and cone gas flow, 25 L/h. For mass accuracy, a 0.2 ng/μl leucine encephalin (mass/charge ratio (m/z) 556.2771 in ESI+ and m/z 554.2615 in ESI−) solution at 20 μl/min was used as the lock mass. Data were collected in centroid mode with a scan range of 50–1200 m/z, with lock mass scans collected every 30 seconds and averaged over 3 scans to perform mass correction. Additionally, data-dependent acquisition (DDA) of the QC sample was performed for structural elucidation for each ionisation mode.

Chromatograms and spectra were displayed using MassLynx (version 4.1; Waters Corp.). The MassLynx package ‘Databridge’ was used to convert the data files into NetCDF format. The package XCMS (R software) and an in-house developed script was used for the pre-processing of the data (peak-picking, alignment, grouping, zero filling, and median-fold change normalisation) [35]. Integrals of metabolic features with a coefficient of variation (CV) greater than 30% in the QC samples (n = 6) were removed in Microsoft Excel (2013, Microsoft, USA) prior to analysis.

Statistical Analysis

Metabolic data from plasma 1H NMR spectra and lipid profiling UPLC-MS data were initially analysed in SIMCA software (V. 13.0, Umetrics) via Principal Component Analysis (PCA) [36], a non-supervised descriptive model used to investigate inherent similarities/differences across the data set and identify outliers. Analyses resulted in the removal of one NMR sample with an acquisition artefact thus bringing the final total to n = 66: 45 HAT and 21 controls (see S1 Table). PCA was additionally used to identify and remove contaminant/artefactual features, as described by Vorkas et al, 2015 [34], and to evaluate experimental reproducibility by assessing QC sample clustering for both ESI+ and ESI- runs. This was followed by Partial Least Squares Discriminant Analysis (PLS-DA), and Orthogonal PLS-DA (O-PLS-DA) [37]. Both are multivariate regression models used to assess the spectral differences that contribute towards classification of patients vs. controls, and to determine how well these differences can predict class separation. Contribution of within-class variation that is uncorrelated to the separation between classes is minimised in O-PLS-DA models, which were used for final extraction of discriminating metabolites. See supplementary protocols in supporting information (S1 File) for further details on data processing, analyses, and model validation.

For NMR data, metabolites integrals were then compared between HAT patients and controls via the two-tailed Welch’s (unequal variance) T-test in Microsoft Excel, followed by Benjamini-Hochberg multiple test correction [38]. Those where p<0.05 was observed were labelled as statistically significant discriminatory markers. For UPLS-MS data, threshold criteria for initial selection of markers were variables with a correlation value p(corr)[1]>0.5 and covariance p[1]>0.05, as visualized on the OPLS-DA S-plot, to extract the most discriminatory and robust features (see supporting protocols). This was followed by conducting a two-tailed Welch’s T-test with Benjamini-Hochberg correction, whereby the final statistical discriminatory criteria was also p<0.05 after multiple testing correction. Lipid species assignment of discriminatory features was based on the accurate mass-to-charge (m/z) ratio, searched against online MS database METLIN (https://metlin.scripps.edu) and previous published results generated in-house [34,39], as well as on the MS/MS information obtained in the DDA experiments of QC samples.

Integrals of these discriminatory variables were re-imported into SIMCA as new loadings for O-PLS-DA models, to evaluate the sensitivity and specificity of the combination of these metabolites for biomarker potential. O-PLS-DA models was performed twice, separately for NMR and UPLC-MS data. Firstly, models were generated using integrals of all significantly discriminatory markers, and secondly using only integrals of the top 5 discriminating markers, defined as those which displayed the greatest % difference between patients and controls. Calculations are further explained in supplementary protocols (see Supporting Information S1 File).

To explore links between potentially complementary metabolic datasets, NMR metabolite integrals were correlated with corresponding integral of discriminatory lipid profiling features from matching patient samples (n = 16), via Pearson correlation in MATLAB, followed by Benjamini-Hochberg correction. Similarly Pearson -based correlation analyses (with multiple test correction) was performed between metabolic (NMR and/or MS) data sets. Relationships between patient demographic data and metabolic datasets were assessed using Spearman correlation analysis, Mann-Whitney tests, and mixed linear models as appropriate, in JMP (v.10.0, SAS, USA).

Results

1H NMR Spectroscopy Shows Marked Differences in Plasma between Human African Trypanosomiasis Patients and Controls

Inherent differences between T. b. rhodesiense infected HAT patients and uninfected controls dominated overall dataset variability, demonstrated by the partial separation of these two populations along the first component in PCA (Fig 1A), indicative that the HAT infection produced a strong metabolic response. No significant differences between metabolic profiles based on demographic factors (gender, age, or diagnostic stage) were found. Metabolic disparities between patients and controls were further examined using O-PLS-DA (a supervised modelling method that maximises differences between classes whilst minimising within-class variation), which showed complete segregation between patients and controls, also along the first component (Fig 1B). The model had a Q2Y (model predictive capability parameter) value of 0.708, indicating that the model was robust. Further validation through permutation testing found that the likelihood of achieving the same discriminatory result by chance was zero (p = 0 for both R2Y and Q2Y). It is noteworthy that although the first two components of the PCA model only cumulatively account for 18.1% of the total variance (R2) in the data set, indicating variability in sample composition, the O-PLS-DA model indicates that the metabolic signature is robust.

thumbnail
Fig 1. Plasma in HAT patients display different 1H NMR metabolic profiles compared with controls.

PCA model (A) and O-PLS-DA model (B) score plots of plasma 1H NMR spectra across HAT patients and controls. Each circle represents a spectra from one sample, whereby patients are presented in red (n = 45) and controls in dark grey (n = 21). Abbreviations: R2X, model fit parameter for variation in spectral data; R2Y, model fit parameter for variation in classifier data (for O-PLS-DA); Q2, model predictive parameter for spectral data in PCA (Q2X) and for classifier data in O-PLS-DA (Q2Y). Individual component contribution of R2X are shown on the axes as percentage.

https://doi.org/10.1371/journal.pntd.0004200.g001

Examination of the corresponding O-PLS-DA loadings revealed that the lipid peaks accounted for the majority of the discrimination. In total, 34 endogenous metabolites were identified from NMR spectra though only 31 of these displayed distinct peaks that were used for integral calculation (signals from isoleucine, leucine and lysine were superimposed by overlapping lipid peaks), and compared between patients and controls (shown in S2 Table). Those metabolites that were significantly higher in disease are shown in Fig 2A, whilst metabolites significantly lower in disease are shown in Fig 2B. The majority of lipid moieties detected in plasma were significantly reduced in patients, including those containing unsaturated bonds (CH = CH and C = CCH2C = C), lysyl group in plasma albumin, and cholesterol (Fig 2B). In contrast, lipid glyceryl and ketene (CH2CO) groups were significantly higher in patients than in controls (Fig 2A). Several plasma amino acids were also significantly discriminatory between HAT and controls, including histidine and alanine, which were found in lower concentrations in patients, and phenylalanine, which was augmented. Other significant changes included an increase of creatinine, N-acetyl glycoprotein, formate and myo-inositol in HAT patients whilst citrate was higher in controls. No associations between discriminatory metabolites and patient demographic data (age, gender or diagnostic stage) were found though the sample size in this study limited the power of this analysis (this was particularly the case for staging as our cohort had relatively few early cases—see S1 Table).

thumbnail
Fig 2. Differences in plasma metabolites between HAT patients and controls detected by NMR.

Bar-charts show relative levels of plasma metabolites that were significantly altered between patients (shown in red, n = 45) and controls (shown in grey, n = 21), as measured by 1H NMR spectroscopy. (A) Metabolites higher in patients than controls. (B). Metabolites lower in patients than controls. Bars represent group mean average with standard error of the mean as error bars. Significance measured via Welch T-test with multiple test correction, shown as asterisks; * p<0.05, ** p<0.01, *** p<0.001. A.U., arbitrary units; NAG, N-acetyl glycoproteins.

https://doi.org/10.1371/journal.pntd.0004200.g002

Additional O-PLS-DA models were built, firstly using integrals of all 16 discriminatory metabolites, followed by models using only the top 5 metabolites. The latter were defined as those which displayed the biggest measured percentage difference between HAT and controls, representing metabolites of most practical diagnostic potential. These consisted of phenylalanine, formate and lipid moieties CH2CO, CH = CH, and glyceryl. In both O-PLS-DA models, results provided identical specificity scores for classifying patients and controls samples, namely 95.24% (with 95% confidence intervals [CI] of 86.1–100%). Model sensitivity was 88.89% (79.7–98.1% CI) using all NMR discriminatory markers, and 84.44% (73.9–95.0% CI) using just the top 5 metabolites (see S3 Table for confusion matrices).

Lipid Profiling UPLC-MS Reveals Compositional Differences in Plasma Lipid Species between HAT and Controls

Lipid plasma profiles were further examined using reversed-phase UPLC-MS. In total, 8771 metabolic features were detected in positive ionisation mode (ESI+), whilst 2543 features were observed in negative mode (ESI-). Following the exclusion of peaks with CV > 30% across the QC’s, as well as the removal of contaminant peaks present in either QC’s or blank controls, 6952 features in ESI+ and 1714 features in ESI- were retained for final analyses. PCA models showed tight clustering of QC samples compared with the remaining test samples (supplementary information S1 Fig), assuring the robustness of the MS experimental runs.

As noted with NMR, lipid profiling results showed an inherent separation between the two classes already in the first two components of PCA, observed in both ionisation modes (see S1 Fig). These differences were confirmed in the corresponding O-PLS-DA models (S1 Fig). Permutation testing confirmed model robustness (ESI+, p = 0.004 for R2Y and p = 0 for Q2Y and for ESI-, p = 0 for both R2Y and Q2Y). Metabolic features responsible for class discrimination (i.e. HAT patients vs. controls) were visualised in S-plots in SIMCA, with the cut-off values of p(corr)[1]>0.5 and covariance p[1]>0.05 were used as selection criteria (see Methods), shown in Fig 3A and 3B. Integrals of selected discriminatory features highlighted in the S-plots are shown in Fig 3C and 3D, demonstrating significant decreases in lysophosphatidylcholines 16:0 and 18:0, as well as increases in phosphatidylcholine 32:0 and 34:1.

thumbnail
Fig 3. Plasma in HAT patients display different lipid profiles compared with controls.

S-plots of O-PLS-DA model for plasma lipid profiling UPLC-MS features detected in ESI+ mode (A) and ESI- mode (B), whereby each circle represents one feature with a unique combination of m/z and retention time values. Discriminatory features selected surpassed p[1] and p(corr)[1] threshold criteria, highlighted in red boxes (see Methods). See S1 Fig for corresponding scores plots. Integrals of features highlighted in the S-plots are shown as bar-chart showing mean averages ± standard error of the mean, for both ESI+ (C) and ESI- (D). Patients levels are shown in red (n = 16) and controls in dark grey (n = 14). Significant differences, as measured via Welch T-test with multiple test correction, are labelled with asterisks, where *** p<0.001. Abbreviations: A.U., arbitrary units; LysoPC, lysophosphatidylcholine; PC, phosphatidylcholine.

https://doi.org/10.1371/journal.pntd.0004200.g003

Overall, lipid profiling analysis identified 37 unique features that were statistically different between HAT and controls, 12 of which were found to be significantly different in both of these modes (Table 1). Features included six lysophosphatidylcholines (LysoPCs 16:0, 18:0, 18:1, 18:2, 20:3 and 20:4), four phosphatidylcholines (PCs 32:0, O-34:1, 38:3 and 38:5) and two sphingomyelin species (SMs d40:1 and d41:1). Further information on UPLC-MS assignments of features can be found in S4 Table. The majority of lipid classes studied were glycerophospholipids (LysoPCs and PCs), followed by sphingolipids (mainly sphingomyelins) and triglycerides (see Table 1 and S4 Table). Free fatty acids were detected using our MS methodology but did not significantly differ between patients and controls.

thumbnail
Table 1. Discriminatory lipid species observed in plasma via UPLC-MS.

https://doi.org/10.1371/journal.pntd.0004200.t001

In total, six significantly altered lysoPCs measured in this study were found to be lower in patients than controls. Significantly lower intensities were also observed for eight out of nine sphingolipids, including seven sphingomyelins (SM) as well as nervonic ceramide (Cer d42:2). Only SM d34:1, the smallest of the analysed SM species, was found to be higher in patients than controls. Additionally, patients had augmented levels of six triglycerides (TG) species and depleted levels of cholesterol ester (CE) 18:2 in patients, in line with NMR findings where we detected higher glyceryl moieties and lower cholesterol in HAT.

In general, differences in lipid species concentrations between patients and controls were similar across the lipid subclasses, with the exception of the phosphatidylcholines. In contrast, our UPLC-MS results indicate that PCs that contained a high number of saturated bonds were significantly decreased in patients (e.g. PCs 38:5, 38:6, 40:5), whilst those with two or fewer unsaturated bonds were increased (e.g. PCs 32:0, 32:1, and 34:1, see Table 1). These results mirror the NMR findings, where lipid moieties with unsaturated bonds (CH = CH and C = CCH2C = C) were significantly lower in patients than controls (Fig 2).

Similar to the NMR data, additional O-PLS-DA models were generated using the integrals all 37 markers from both ionisation modes, followed by further models using only the top 5 features with the greatest % difference (as shown in Table 1). As most of the largest changes were observed in the positive mode, we selected the feature from this ionisation mode only. These included PC(34:3), PC(32:0), TG(51:2), PC(O-34:1) and SM(d34:1). Both O-PLS-DA models with either all markers or top 5 gave identifical sensitivity and specificity scores, namely 93.75% (81.9–100% CI) and 85.71% respectively (67.4–100% CI; see S5 Table for confusion matrices).

Potential links between UPLC-MS and NMR data was explored by correlating lipid metabolic features with NMR metabolite levels measured via NMR. Whilst multiple trends were observed, including links between lipid MS species and several NMR lipid moieties, 3-hydroxybutyrate (a ketone body generated from fatty acid break-down), and energy associated metabolites, e.g. creatinine and lactate), none of these associations remained statistically significant after Benjamini-Hochberg multiple test correction was performed, due to lack of statistical power. No significant links were observed between MS data and age and/or gender data (analyses could not be performed for diagnostic staging due to limited numbers of early stage within this subset of data—see S2 Table).

Discussion

The HAT infection was characterised by strong metabolic differences in plasma profiles, based on both NMR and UPLC-MS, which shows the potential of this approach in the development of novel diagnostic marker sets. The discriminatory profile was largely driven by variation in lipids and amino acids. Phenylalanine was the single most discriminatory metabolite as detected by NMR, with the largest percentage difference between HAT and controls in concentration. Raised levels of phenylalanine have been frequently linked with infection and inflammation [40]. Phenylalaninemia has previously been reported in patients with enteric fever [41], and with malaria [42] Evidence suggests that this phenomenon may be related to neurological effects seen in cerebral malaria [42], since high levels of phenylalanine are known to be neurotoxic [43]. Other differences in plasma amino acid concentrations included significant decreases in histidine and alanine in HAT. Histidine concentrations have been shown to be negatively correlated with inflammatory disease markers in humans [4446] and this would be consistent with clinical evidence of systemic inflammatory responses in HAT [3].

Creatinine concentration was also found to be significantly higher in HAT, potentially a sign of renal dysfunction which is one of the clinical features observed in sleeping sickness [47,48]. Additionally, concentrations of N-acetyl groups of plasma glycoproteins were also higher in patients. Most of these have been shown to belong to α1 acid glycoprotein [49,50], another inflammatory marker known to be associated with the acute-phase response. This result is similar to what has been found with experimental murine infection of T. b. brucei where plasma acetylated glycoproteins (dominated by O-acetyls) were found to be augmented [11].

Some of the most pronounced changes observed in the NMR spectra were the clear disparities in the plasma lipid composition between the two groups. In addition to their roles as energy reserves and major constituents of membranes, lipids also represent highly biologically active metabolites, with involvement (or function) in signalling and a range of inflammatory processes, e.g. prostaglandins and leukotrienes synthesis [51]. Our results showed that the majority of the concentrations of discriminatory lipids were significantly reduced in patients. These primarily included lipids containing unsaturated bonds (CH = CH and C = CCH2C = C), CH3 and CH2CH2CO moieties, and albumin lysyl groups, and high-density lipoprotein (HDL) cholesterol. Dyslipidaemia has been observed in a range of infections, measured primarily using lipid-specific biochemical serum assays, as opposed to comprehensive metabolic profiling, whereby significantly reduced HDL and total cholesterol have been reported in malaria [52], enteric fever [53] and others [54,55]. Decreased lipoprotein CH3 moieties were also observed in experimental murine T. b. brucei infection [11], confirming the need for deeper assessment into the variety of lipid species, to help dissect their different roles during infection.

Further investigation of the lipid profile with UPLC-MS indicated that HAT patients displayed significantly lower levels of polyunsaturated phosphatidylcholines (e.g. PC 40:6, 40:5, 38:5 and 38:3), although PC molecules with 2 or fewer unsaturated bonds were in contrast increased. This profile presents a more complex picture than was obtained in experimental mouse infections with T. b. brucei, where an overall decrease in plasma PC concentrations was measured by NMR [11].

A range of lysophophatidylcholines (e.g. LysoPC 16:0, 18:0, 18:2) intensities was also significantly lower in patients. Despite the fact that Trypanosoma brucei ssp have the capability to produce many of their lipids de novo, they are known to scavenge several lipid precursor molecules from the host and have been shown to require lipid uptake for optimal growth [56,57]. More specifically, these include fatty acids and lyso-phospholipids, which are assembled to make phospholipids, the most dominant lipid component of T. brucei ssp [58]. As fatty acyl chain composition in T. b. rhodesiense bloodstream forms have a high proportion of long-chain polyunsaturated fatty acids [59,60], it could be hypothesised that the lower plasma concentrations of polyunsaturated lipids in HAT patients may be due to increased uptake by the parasites.

Numerous sphingomyelin species were also found in lower concentrations in infection (e.g. SM d40:1, d41:1 and d42:1), accompanied by a significant increase of ceramide. The sphingolipid turnover is regulated by the balance between enzymatic degradation and synthesis, with ceramide being a major break-down product. Oxidative stress, inflammatory cytokines and infection can trigger sphingomyelinase enzymes to increase ceramide generation. Higher levels of ceramide, in turn, have been linked to increased cell autophagy and apoptosis (the latter particularly in endothelial cells leading to increased vascular permeability), additional pro-inflammatory cytokine and chemokine synthesis, and other metabolic disorders (e.g. insulin resistance and obesity) [61]. Thus, the changes in sphingolipids observed in HAT further highlights the widespread inflammatory effects on the host, associated with T. b. rhodesiense infection.

Another characteristic of HAT patients was hypertriglyceridaemia. Not only is this result is in agreement with a range of rabbit and non-human primate models of T. brucei spp infections [21,22], but is also consistent with marked increases in serum triglycerides in T. b. gambiense HAT [62]. However, the disease discriminatory potential of hypertriglyceridaemia (often combined with increases in circulating free fatty acids) will require further investigation as it has been shown to be a common feature across numerous infections in humans [5254].

In summary, the presented work is the first to apply a comprehensive metabolic profiling approach to the investigation of Trypanosoma brucei infection in humans and highlights the potential of metabonomic technology for developing a diagnostic platform for HAT. The models generated based on both NMR and UPLC-MS were robust, thus the utility of this technology in an exploratory capacity is evident. Similarities but importantly significant differences in discriminatory NMR metabolites were identified to previously published metabolic profiles in experimental rodent T. brucei spp. infections. We identified a range of discriminatory plasma metabolites in clinical HAT that collectively have diagnostic disease marker potential (up to 93.75% sensitivity and 95.24% specificity). While some of these markers have been associated with other infections or inflammatory processes, the effects of specific infection on a combined panel such as in this study has not previously been reported and may prove diagnostically effective in differentiating HAT from other endemics disease. However, it is obvious that neither NMR spectroscopy nor UPLC-MS are platforms that can be used at point-of-care in developing countries. Therefore the functionality of these platforms lies in characterisation of the metabolic consequences of a parasitic infection, with subsequent development of simpler clinical assays such as a biochemical dipstick based on the strongest differentiating biomarkers. Of particular interest from our results were differences in lipids, observed in both NMR and UPLC-MS, that dominated the variation in metabolic profiles between HAT and controls. Further investigation of inflammatory lipid mediators such prostaglandins and leukotrienes would also be of particular interest, to better understand pathology and provide key insights into links between lipid metabolism and immune responses in the host. Detailed lipid sub-species characterisation with the use of lipid reference standards and MS/MS are required to further isolate novel diagnostic candidates to ultimately aid in the detection and management of T. b. rhodiesiense HAT, moving closer towards the WHO goal of disease eradication.

Supporting Information

S1 Fig. Differences observed in plasma lipid-MS profiles between HAT patients and controls.

Multivariate model scores plots based on lipid MS features of plasma samples. PCA models based on first two component of all plasma test samples (as light grey circles, n = 30) and QC’s (as blue circles, n = 6) in positive ionisation (ESI+) mode (A) and negative (ESI-) mode (B). PCA models of the first two components of plasma test samples only, comparing HAT patients (red circles, n = 16) vs. controls (dark grey circles, n = 14), in ESI+ (C) and ESI- (D). O-PLS-DA models of test samples comparing HAT vs. controls (labelling and group sizes as described for PCA), in ESI+ (E) and ESI- (F). O-PLS-DA based on both predictive and orthogonal components for ESI+ but only on the first predictive component for ESI- (the orthogonal component for this model reduced the predictive ability of model and was thus omitted). Abbreviations: R2X, model fit parameter for variation in spectral data; R2Y, model fit parameter for variation in classifier data (for O-PLS-DA models); Q2, model predictive parameter (for spectral data Q2X in PCA/for classifier data Q2Y in O-PLS-DA), SD, standard deviation.

https://doi.org/10.1371/journal.pntd.0004200.s001

(TIF)

S1 Table. Patient and control cohort demographic information.

https://doi.org/10.1371/journal.pntd.0004200.s002

(PDF)

S2 Table. Plasma metabolite levels measured via 1H NMR spectroscopy.

https://doi.org/10.1371/journal.pntd.0004200.s003

(PDF)

S3 Table. Confusion matrices of HAT vs. control classification ability of NMR discriminatory markers.

https://doi.org/10.1371/journal.pntd.0004200.s004

(PDF)

S4 Table. List of Discriminatory Lipid Species detected via UPLC-MS in both ionisation modes.

https://doi.org/10.1371/journal.pntd.0004200.s005

(PDF)

S5 Table. Confusion matrices of HAT vs. control classification ability of NMR discriminatory markers.

https://doi.org/10.1371/journal.pntd.0004200.s006

(PDF)

S1 File. Supplementary protocols for statistical analyses.

https://doi.org/10.1371/journal.pntd.0004200.s007

(DOCX)

Author Contributions

Conceived and designed the experiments: JMS SDL MGR JS EH. Performed the experiments: SDL MGR. Analyzed the data: SDL PAV MGR JMS. Contributed reagents/materials/analysis tools: JMS VPA. Wrote the paper: SDL JMS EH MGR JS PAV. Contributed to data processing: SDL PAV.

References

  1. 1. Chappuis F, Loutan L, Simarro P, Lejon V, Buscher P. (2005) Options for field diagnosis of human african trypanosomiasis. Clin Microbiol Rev 18: 133–146. pmid:15653823
  2. 2. Fèvre E,M., Wissmann Bv, Welburn SC, Lutumba P. (2008) The burden of human african trypanosomiasis. PLoS Neglected Tropical Diseases 2: e333. pmid:19104653
  3. 3. Sternberg JM. (2004) Human african trypanosomiasis: Clinical presentation and immune response. Parasite Immunol 26: 469–476. pmid:15771682
  4. 4. Sternberg JM, Gierliński M, Biéler S, Ferguson MAJ, Ndung'u JM. (2014) Evaluation of the diagnostic accuracy of prototype rapid tests for human african trypanosomiasis. PLoS Negl Trop Dis 8: e3373. pmid:25521120
  5. 5. World Health Organization. (2012) Accelerating work to overcome the global impact of neglected tropical diseases: A roadmap for implementation.
  6. 6. Simarro PP, Franco JR, Diarra A, Ruiz Postigo JA, Jannin J. (2013) Diversity of human african trypanosomiasis epidemiological settings requires fine-tuning control strategies to facilitate disease elimination. Res Rep Trop Med 4: 1–6.
  7. 7. Lindon JC, Nicholson JK, Holmes E, Antti H, Bollard ME, et al. (2003) Contemporary issues in toxicology the role of metabonomics in toxicology and its evaluation by the COMET project. Toxicol Appl Pharmacol 187: 137–146. pmid:12662897
  8. 8. Teng R, Junankar PR, Bubb WA, Rae C, Mercier P, et al. (2009) Metabolite profiling of the intraerythrocytic malaria parasite plasmodium falciparum by 1H NMR spectroscopy. NMR Biomed 22: 292–302. pmid:19021153
  9. 9. Lamour SD, Choi B, Keun H, Mueller I, Saric J. (2012) Metabolic characterization of Leishmania major infection in activated and non-activated macrophages. J Proteome Res.
  10. 10. Vincent IM, Creek DJ, Burgess K, Woods DJ, Burchmore RJ, et al. (2012) Untargeted metabolomics reveals a lack of synergy between nifurtimox and eflornithine against trypanosoma brucei. PLoS Negl Trop Dis 6: e1618. pmid:22563508
  11. 11. Wang Y, Utzinger J, Saric J, Li JV, Burckhardt J, et al. (2008) Global metabolic responses of mice to trypanosoma brucei brucei infection. Proc Natl Acad Sci U S A 105: 6127–6132. pmid:18413599
  12. 12. Ghosh S, Sengupta A, Sharma S, Sonawat HM. (2013) Metabolic perturbations of kidney and spleen in murine cerebral malaria: 1H NMR-based metabolomic study. PLoS ONE 8: e73113. pmid:24039868
  13. 13. Girones N, Carbajosa S, Guerrero NA, Poveda C, Chillon-Marinas C, et al. (2014) Global metabolomic profiling of acute myocarditis caused by trypanosoma cruzi infection. PLoS Negl Trop Dis 8: e3337. pmid:25412247
  14. 14. Tritten L, Keiser J, Karwa T, Utzinger J, Holmes E, et al. (2014) Comparing systemic metabolic responses in mice to single or dual infection with plasmodium berghei and heligmosomoides bakeri. Mol BioSyst 10: 2358–2367. pmid:24960299
  15. 15. Lamour SD, Veselkov KA, Posma JM, Giraud E, Rogers ME, et al. (2015) Metabolic, immune, and gut microbial signals mount a systems response to leishmania major infection. J Proteome Res 14: 318–329. pmid:25369177
  16. 16. Balog CIA, Meissner A, Göraler S, Bladergroen MR, Vennervald BJ, et al. (2011) Metabonomic investigation of human schistosoma mansoni infection. Mol BioSyst: 1742–206X.
  17. 17. Saric J, Li JV, Swann JR, Utzinger J, Calvert G, et al. (2010) Integrated cytokine and metabolic analysis of pathological responses to parasite exposure in rodents. J Proteome Res 9: 2255–2264. pmid:20092362
  18. 18. Lejon V, Lardon J, Kenis G, Pinoges L, Legros D, et al. (2002) Interleukin (IL)-6, IL-8 and IL-10 in serum and CSF of trypanosoma brucei gambiense sleeping sickness patients before and after treatment. Transactions of The Royal Society of Tropical Medicine and Hygiene 96: 329–333. pmid:12174791
  19. 19. MacLean L, Chisi JE, Odiit M, Gibson WC, Ferris V, et al. (2004) Severity of human african trypanosomiasis in east africa is associated with geographic location, parasite genotype, and host inflammatory cytokine response profile. Infection and Immunity 72: 7040–7044. pmid:15557627
  20. 20. Dumas M, Maibaum EC, Teague C, Ueshima H, Zhou B, et al. (2006) Assessment of analytical reproducibility of 1H NMR spectroscopy based metabonomics for large-scale epidemiological research: The INTERMAP study. Anal Chem 78: 2199–2208. pmid:16579598
  21. 21. Rouzer CA, Cerami A. (1980) Hypertriglyceridemia associated with trypanosoma brucei brucei infection in rabbits: Role of defective triglyceride removal. Mol Biochem Parasitol 2: 31–38. pmid:7464858
  22. 22. Nakamura Y. (1998) Alterations of serum lipid, lipoprotein and inflammatory cytokine profiles of rabbits infected with trypanosoma brucei brucei. Vet Parasitol 80: 117–125. pmid:9870364
  23. 23. Gaithuma AK, Karanja SM, Ngotho M, Maathai RG, Kagira JM, et al. (2012) Lipid metabolism and other metabolic changes in vervet monkeys experimentally infected with trypanosoma brucei rhodesiense. J Med Primatol 41: 75–81. pmid:22070162
  24. 24. MacLean LM, Odiit M, Chisi JE, Kennedy PGE, Sternberg JM. (2010) Focus-specific clinical profiles in human african trypanosomiasis caused by Trypanosoma brucei rhodesiense. PLoS Negl Trop Dis 4: e906. pmid:21151878
  25. 25. World Health Organization. (1998) Control and surveillance of african trypanosomiasis: Report of a WHO expert committee. WHO Technical Report Series.
  26. 26. Beckonert O, Keun HC, Ebbels TM, Bundy J, Holmes E, et al. (2007) Metabolic profiling, metabolomic and metabonomic procedures for NMR spectroscopy of urine, plasma, serum and tissue extracts. Nat Protoc 2: 2692–2703. pmid:18007604
  27. 27. Veselkov KA, Lindon JC, Ebbels TM, Crockford D, Volynkin VV, et al. (2009) Recursive segment-wise peak alignment of biological (1)h NMR spectra for improved metabolic biomarker recovery. Anal Chem 81: 56–66. pmid:19049366
  28. 28. Dieterle F, Ross A, Schlotterbeck G, Senn H. (2006) Probabilistic quotient normalization as robust method to account for dilution of complex biological mixtures. application in 1H NMR metabonomics. Anal Chem 78: 4281–4290. pmid:16808434
  29. 29. Barton RH, Waterman D, Bonner FW, Holmes E, Clarke R, et al. (2010) The influence of EDTA and citrate anticoagulant addition to human plasma on information recovery from NMR-based metabolic profiling studies. Mol Biosyst 6: 215–224. pmid:20024083
  30. 30. Nicholson JK, Foxall PJD, Spraul M, Farrant RD, Lindon JC. (1995) 750 MHz 1H and 1H-13C NMR spectroscopy of human blood plasma. Anal Chem 67: 793–811. pmid:7762816
  31. 31. Liu M, Tang H, Nicholson JK, Lindon JC. (2002) Use of 1H NMR-determined diffusion coefficients to characterize lipoprotein fractions in human blood plasma. Magn Reson Chem 40: S83–S88.
  32. 32. Sarafian MH, Gaudin M, Lewis MR, Martin F, Holmes E, et al. (2014) Objective set of criteria for optimization of sample preparation procedures for ultra-high throughput untargeted blood plasma lipid profiling by ultra performance liquid Chromatography–Mass spectrometry. Anal Chem 86: 5766–5774. pmid:24820162
  33. 33. Isaac G, McDonald S, Astarita G. (2011) Lipid separation using UPLC with charged surface hybrid technology. Waters App note.
  34. 34. Vorkas PA, Isaac G, Anwar MA, Davies AH, Want EJ, et al. (2015) Untargeted UPLC-MS profiling pipeline to expand tissue metabolome coverage: Application to cardiovascular disease. Anal Chem 87: 4184–4193. pmid:25664760
  35. 35. Smith CA, Want EJ, O'Maille G, Abagyan R, Siuzdak G. (2006) XCMS:Processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal Chem 78: 779–787. pmid:16448051
  36. 36. Hotelling H. (1933) Analysis of a complex of statistical variables into principal components. J Educ Psychol 24: 417.
  37. 37. Bylesjö M, Rantalainen M, Cloarec O, Nicholson JK, Holmes E, et al. (2006) OPLS discriminant analysis: Combining the strengths of PLS-DA and SIMCA classification. J Chemometrics 20: 341–351.
  38. 38. Benjamini Y, Hochberg Y. (1995) Controlling the false discovery rate: A practical and powerful approach to multiple testing. Journal of the Royal Statistical Society.Series B (Methodological) 57: 289–300.
  39. 39. Vorkas PA, Shalhoub J, Isaac G, Want EJ, Nicholson JK, et al. (2015) Metabolic phenotyping of atherosclerotic plaques reveals latent associations between free cholesterol and ceramide metabolism in atherogenesis. J Proteome Res 14: 1389–1399. pmid:25565173
  40. 40. Wannemacher RW, Klainer AS, Dinterman RE, Beisel WR. (1976) The significance and mechanism of an increased serum phenylalanine-tyrosine ratio during infection. The American Journal of Clinical Nutrition 29: 997–1006. pmid:822705
  41. 41. Nasstrom E, Vu Thieu NT, Dongol S, Karkey A, Voong Vinh P, et al. (2014) Salmonella typhi and salmonella paratyphi A elaborate distinct systemic metabolite signatures during enteric fever. Elife 3:
  42. 42. Lopansri BK, Anstey NM, Stoddard GJ, Mwaikambo ED, Boutlis CS, et al. (2006) Elevated plasma phenylalanine in severe malaria and implications for pathophysiology of neurological complications. Infect Immun 74: 3355–3359. pmid:16714564
  43. 43. Kaufman S. (1989) An evaluation of the possible neurotoxicity of metabolites of phenylalanine. J Pediatr 114: 895–900. pmid:2654351
  44. 44. Enwonwu CO, Afolabi BM, Salako LO, Idigbe EO, Bashirelahi N. (2000) Increased plasma levels of histidine and histamine in falciparum malaria: Relevance to severity of infection. J Neural Transm 107: 1273–1287. pmid:11145003
  45. 45. Watanabe M, Suliman ME, Qureshi AR, Garcia-Lopez E, Bárány P, et al. (2008) Consequences of low plasma histidine in chronic kidney disease patients: Associations with inflammation, oxidative stress, and mortality. The American Journal of Clinical Nutrition 87: 1860–1866. pmid:18541578
  46. 46. Niu YC, Feng RN, Hou Y, Li K, Kang Z, et al. (2012) Histidine and arginine are associated with inflammation and oxidative stress in obese women. Br J Nutr 108: 57–61. pmid:21996294
  47. 47. Brun R, Blum J, Chappuis F, Burri C. (2010) Human african trypanosomiasis. The Lancet 375: 148–159.
  48. 48. Omeje JN, Anene BM. (2012) Comparative serum biochemical changes induced by experimental infection of T. brucei and T. congolense in pigs. Vet Parasitol 190: 368–374. pmid:22858639
  49. 49. Bell JD, Brown JCC, Nicholson JK, Sadler PJ. (1987) Assignment of resonances for ‘acute-phase’ glycoproteins in high resolution proton NMR spectra of human blood plasma. FEBS Lett 215: 311–315. pmid:2438159
  50. 50. Fischer K, Kettunen J, Würtz P, Haller T, Havulinna AS, et al. (2014) Biomarker profiling by nuclear magnetic resonance spectroscopy for the prediction of all-cause mortality: An observational study of 17,345 persons. PLoS Med 11: e1001606. pmid:24586121
  51. 51. Calder PC. (2001) Polyunsaturated fatty acids, inflammation, and immunity. Lipids 36: 1007–1024. pmid:11724453
  52. 52. Visser BJ, Wieten RW, Nagel IM, Grobusch MP. (2013) Serum lipids and lipoproteins in malaria-a systematic review and meta-analysis. Malaria J 12: 1–16.
  53. 53. Khosla SN, Goyle N, Seth RK. (1991) Lipid profile in enteric fever. J Assoc Physicians India 39: 260–262. pmid:1880095
  54. 54. Alvarez C, Ramos A. (1986) Lipids, lipoproteins, and apoproteins in serum during infection. Clin Chem 32: 142–145. pmid:3940695
  55. 55. Soares NM, Leal TF, Fiuza MC, Reis EA, Souza MA, et al. (2010) Plasma lipoproteins in visceral leishmaniasis and their effect on leishmania-infected macrophages. Parasite Immunol 32: 259–266. pmid:20398226
  56. 56. Black S, Vandeweerd V. (1989) Serum lipoproteins are required for multiplication of trypanosoma brucei brucei under axenic culture conditions. Mol Biochem Parasitol 37: 65–72. pmid:2615788
  57. 57. Coppens I, Levade T, Courtoy PJ. (1995) Host plasma low density lipoprotein particles as an essential source of lipids for the bloodstream forms of trypanosoma brucei. J Biol Chem 270: 5736–5741. pmid:7890701
  58. 58. Richmond GS, Gibellini F, Young SA, Major L, Denton H, et al. (2010) Lipidomic analysis of bloodstream and procyclic form trypanosoma brucei. Parasitology 137: 1357–1392. pmid:20602846
  59. 59. Dixon H, Williamson J. (1970) The lipid composition of blood and culture forms of trypanosoma lewisi and trypanosoma rhodesiense compared with that of their environment. Comp Biochem Physiol 33: 111–128. pmid:5440923
  60. 60. Smith TK, Bütikofer P. (2010) Lipid metabolism in trypanosoma brucei. Mol Biochem Parasitol 172: 66–79. pmid:20382188
  61. 61. Maceyka M, Spiegel S. (2014) Sphingolipid metabolites in inflammatory disease. Nature 510: 58–67. pmid:24899305
  62. 62. Huet G, Lemesre JL, Grard G, Noireau F, Boutignon F, et al. (1990) Serum lipid and lipoprotein abnormalities in human african trypanosomiasis. Trans R Soc Trop Med Hyg 84: 792–794. pmid:2096509