parasitic infections induce a variety of biochemical responses in the host , including modulation of energy metabolism , inflammation , anemia , and nutrient deficiencies . the effective control and eventual elimination of these invading organisms is generally mediated by both cellular and humoral immune mechanisms , and includes different t cell - effector mechanisms , such as those reflected by the increase in th1 , th17 , th2 subtype cd4 cells , as well as regulatory t - cells . the th1/th2 paradigm originally postulated by mosmann and coffman in 1986 described a framework for understanding immunological control mechanisms,(20 ) which were simplified into two main mechanisms : th1 or  cell - mediated immunity and th2 or intracellular parasites ( e.g. , bacteria and protozoa ) are effectively fought by a th1 response and the promotion of direct cytotoxic activity , which , in turn , can lead to tissue damage in the host . helminth infections , in contrast , are controlled by th2-based immune mechanisms , which involve secretion of antibodies , promotion of eosinophilia ( il-5),(21 ) production of ige and igg1 ( il-4 ) , mastocytosis , and mediation of helminth expulsion . however , this original th1/th2 dogma is somewhat oversimplified and it is necessary to consider emerging evidence of other immunological processes . for example , th17 subtype cd4 cells , which are reported to control microparasites effectively ( e.g. , bacteria and fungi ) also promote certain autoimmune mechanisms , whereas t - regulatory cells ( tregs ) are activated in an exclusive fashion to th17 cells and fine - tune the immune - response of the host . the immune outcome and the associated cytokine profile manifested by the host in response to parasitic infection are dependent upon the parasitic species , intensity of infection , the first organ of antigen exposure and also the cytokine milieu already present within the host.(26 ) in addition to immunologic changes , parasitic infections cause direct cellular events such as altered energy metabolism . metabolic profiling strategies combining spectroscopic fingerprinting with mathematical modeling have previously been applied to characterize the global response of the host at the metabolic level and have been discussed in the light of specific measured inflammatory components in particular disease conditions , such as breast cancer(27 ) and colitis,(28 ) but thus far metabolic signatures of pathology have not been systematically modeled with immunological profiles . several single - parasite infection models have been investigated in rodents , using h nuclear magnetic resonance ( nmr ) spectroscopy to generate biochemical signatures of various biological fluids and tissues . these studies have defined panels of metabolic biomarkers reflective of the various components of the biological response to specific parasites . in addition to biochemical changes induced as a direct consequence of tissue pathology , some components of these metabolic responses are likely to be directly linked to immunologic events , in particular the inflammatory response . for example , a trypanosoma brucei brucei infection in the mouse has been associated with increased circulating levels of acetylated glycoproteins,(29 ) which are normally secreted in the liver in response to pro - inflammatory cytokines . conversely , increased plasma levels of lactate and depleted glucose have also been documented to be a feature of a t. b. brucei infection in the mouse and imply increased glycolytic activity in the parasite independent of the inflammation.(29 ) here , we assess covariance between immunological parameters as defined by the pattern of multiple cytokines , and the metabolic phenotypes , as defined by the h nmr spectral profiles , in biofluids obtained from four parasiterodent models . we propose a new strategy for correlating the plasma concentrations of several selected th1/th2 cytokines with the biochemical signature of the different parasiterodent models , including protozoan ( th1-inducing ) and helminth ( th2-inducing ) infections and elucidate model - dependent and model - independent associations of inflammatory measures with metabolic phenotypes with a view to highlighting the component of the metabolic signature , which is associated with , or modulated by , parasite - specific host inflammatory responses . the levels of selected plasma cytokines were coanalyzed with pre - existing spectroscopic data from four different parasiterodent models . all animal experiments were carried out at the swiss tropical institute ( basel , switzerland ) and complied with swiss local and national regulations on animal welfare ( permission no . three separate batches of nmri female mice were purchased from rcc ( itingen , switzerland ) . animals were kept under environmentally controlled conditions ( light / dark cycle , 12/12 h ; temperature , 22 c ; relative humidity , 6070% ) with unrestricted access to community tap water and standard rodent diet ( nafag ; gossau , switzerland ) . sample collection in the plasmodium bergheimouse , t. b. bruceimouse and schistosoma mansonimouse models have been described elsewhere . , in brief , urine samples were collected into a petri dish by gently rubing the abdomen until a minimum volume of 20 l was obtained . whole blood was collected from the tail tip of the mice using sodium heparin - coated hematocrit capillaries . the plasma was separated from the cellular components by centrifuging the whole blood sample at 4,000 g for 4 min . both urine and plasma samples were transferred into eppendorf tubes on dry ice immediately after collection prior to storage at 40 c . additionally , a fasciola hepaticarat model was included , which allowed subsequent comparison between results obtained from mice and rats . the experimental details from each of the parasiterodent models are summarized in table 1 and a reference provided to indicate the rationale for selection of the infection intensities and sampling time points for each model . for the cytokine analysis , 5 infected and 5 noninfected control animals were analyzed for the t. b. bruceimouse model ( 28 days postinfection ) , and 5 infected and 4 noninfected control animals for the f. hepaticarat model ( 22 days postinfection ) . in the case of s. mansonimouse ( 53 days postinfection ) , 7 infected and 5 noninfected animals were used due to the observed variation in metabolic profiles in the infected group . since the duration of the p. berghei model was only 4 days , an age - matched control group was deemed unnecessary , particularly in view of the 3r rules ( reduce , replace , and refine ) , and hence the postinfection samples group ( day 4 ) was directly compared with the preinfection samples group . animals : ( c ) , non - infected animals ; ( t ) , infected animals ; and ( p ) , pre - infection ;  increased upon infection ,  decreased upon infection . altered cytokines : top 4 cytokines that co - map with the infection in the pca analysis . * urine and plasma samples were analyzed using standard one - dimensional ( 1d ) h nmr spectroscopy at 600 mhz with suppression of the water signal , according to a published protocol.(30 ) cytokine measurements were made using the  meso scale discovery ( msd)-based multi cytokine approach.(34 ) the commercial mouse - specific 96-well multispot plate of th1/th2 cytokines ( tnf- , mkc , il-5 , il-4 , il-2 , il-1 , il-12 , il-10 , and ifn- ) was used for the murine models and the corresponding rat demonstration kit ( ifn- , il-1 , il-4 , il-5 , il-13 , kc / gro , and tnf- ) was selected for the f. hepaticarat model . the concentrations of cytokines were quantified using an 8-point calibration curve constructed from a plot of the signal intensity for a series of known concentrations of the multiplex standard provided by the kit manufacturer . the data were analyzed using msd workbench software and the curves fitted using the built - in 4pl fit with 1/y weighting ( y = b2 + ( ( b1  b2)/(1 + ( x / b3 ) ) ) . validation of the msd platform and comparison with luminex , a different multiplex technology has been conducted previously.(35 ) each plasma sample was measured in duplicate and an aliquot of 20 l placed into separate wells . the measurement procedure was performed according to the manufacturer s protocol with the exception that the incubation time used was 4 h. plasma samples from all selected animals were obtained after nondestructive h nmr spectral acquisition ( e.g. , a 1:1 dilution in saline , containing 90% deuterated water , was implemented ) ; hence , the incubation time was extended from 2 to 4 h in order to compensate for the dilution factor and to increase the sensitivity of the assay . a msd sector imager ( sector imager 6000 ; gaithersburg , md ) was utilized to read the plates . principal component analysis ( pca ) was applied to each individual parasitehost data set and also to the combined murine data for s. mansoni , p. berghei , and t. b. brucei infections . since pca is an unsupervised method and makes no assumption as to the origin or class of samples , it allows the major sources of variance in a data set to be defined without incorporating inherent bias . application to a complex matrix , such as those represented here by the nmr spectral data or the cytokine data ( x - matrix ) , achieves a reduction in the dimensionality of the data matrix , identifies outliers and facilitates visualization of the data distribution . to maximize the variance in the nmr data related to the response to the infection , projection to latent structure discriminant analysis ( pls - da ) was conducted on both the nmr data and the cytokine data separately , using a binary code to indicate infection or noninfection as the y - matrix . the information on class membership allows optimal extraction of the particular observations ( e.g. , metabolites ) , which are responsible for the differentiation of groups and can maximize separation in the pls - da model . however , pls - da can also include variation that is not response class or y - related . by incorporation of an orthogonal filter ( o - pls ) , systematic variation unrelated to infection status can be removed and thereby facilitates the interpretation of infection - related effects.(38 ) here , an in - house adapted o - pls script was used(39 ) to perform a cross - matrix correlation between h nmr - derived spectral information ( x ) and identity matched single cytokine levels ( y ) in order to find connections between metabolic and immune metrics . the establishment of infection in the 4 parasiterodent models employed was verified by hematologic parameters and histological examination . all parasitic infections induced clear pathology in all infected animals as described in table 1 , with the exception of the s. mansonimouse model , where the establishment of infection could not be confirmed for 2 out of 7 infected mice . images of hematoxylin and eosin stained selected sections of host organs are provided for each of the four infection models in figure 1 , showing pathological events in spleen tissue sections from p. berghei ( figure 1a ) and s. mansoni - infected mice ( figure 1c ) and in liver sections from f. hepatica - infected rats ( figure 1b ) , and t. b. brucei - infected mice ( figure 1d ) . pathological changes were also observed in the kidney of t. b. brucei - infected mice,(30 ) which showed interstitial nephritis and f. hepatica - infected rats demonstrated splenic follicular hyperplasia with inflammatory infiltrates found in the kidney ( unpublished data ) ( table 1 ) . ( a ) reactive follicular hyperplasia detected in a mouse spleen 4 days after infection with p. berghei . red pulp congestion was further indicated by abundant hemozoin in cordal macrophages ; ( b ) liver tissue sections from a f. hepatica - infected rat showing typical cavity - like necrotic tunnels that were demarked by a foreign body reaction ( 100 ) ; ( c ) lymphofollicular hyperplasia was detected in the spleen of s. mansoni - infected mice and indicated by prominent germinal centers including tangible body macrophages ; ( d ) hepatic hemopoiesis after infection with t. b. brucei . the parasite - infected animals were easily differentiated from noninfected controls using pca and o - pls - da models and the disease - related biological markers were extracted from urinary and plasma metabolite profiles . metabolic characterization of the t. b. bruceimouse,(30 ) and p. bergheimouse(3 ) models has been conducted and described previously . in brief , the infection with both protozoan parasites induced depleted plasma glucose concentrations and subsequently increased lactate levels , reflecting a higher glycolytic activity postinfection . increased levels of o - acetyl glycoproteins were found in the plasma of t. b. brucei , but not p. berghei - infected mice . the urinary manifestation of both protozoan species shared some similarities , including an increase in 2-oxoisovalerate levels , but also elicited disease specific markers , such as increased pipecolic acid ( p. berghei ) and increased 4-hydroxyphenyl acetic acid ( t. b. brucei ) . f. hepatica induced systematic changes in the plasma metabolite profile as illustrated in the o - pls - da coefficient plot ( figure 2a ) . f. hepatica infection caused a decrease in plasma levels of tmao , and increased relative levels of the o- and n - acetyl glycoproteins ( figure 2a ) , with a concomitant depletion of urinary creatinine and taurine . the scores from both the o - pls - da and the pca models constructed from the plasma h nmr spectral data of f. hepatica - infected rats on day 22 postinfection , when overt pathology was established , demonstrated discrete clustering of h nmr plasma spectra relating to control and infected rats ( figure 2b ) ( o - pls - da : ry = 0.90 , qy = 0.57 ; pca : rcumx = 0.96 , qx = 0.40 ) . one o - pls component and one orthogonal component were used for the o - pls - da model , whereas for the pca analysis three principal components ( pcs ) were used to construct the model . differentiation of plasma samples obtained from a pairwise comparison of plasma h nmr data of f. hepatica - infected vs. noninfected control rats at day 22 postinfection . ( a ) o - pls - da coefficient plot ; ( b ) pca of h nmr plasma data ; ( c ) pca scores plot ; and ( d ) pca coefficients plot based on the relative levels of selected cytokines assessed on day 22 postinfection in f. hepatica - infected rats ( orange ) and noninfected control rats ( black ) . global pca and pls - da analyses on the combined spectral data sets derived from the three murine models showed that the samples from the protozoan infections were biochemically distinct from samples from s. mansoni - infected mice , which formed a separate cluster along the axis of the third pc ( figure 3a and supporting information ) . the pls - da model was constructed using three pcs ( pls - da : rx = 0.40 , ry = 0.71 , qy = 0.33 ) . from the scores and the corresponding loadings , it was apparent that the biochemical composition of the control samples relating to each parasitehost model differed according to batch , particularly with respect to t. b. brucei . these observed batch differences are unsurprising given that an outbred mouse strain was used as the host . furthermore , age - dependent variation might contribute to the differences within the control groups ; all control animals were approximately the same age prior to the infection , but due to the different length of time required to establish a patent infection , the time of sampling ranged from 4 days postinfection in the shortest model ( p. berghei ) to 53 days postinfection in the longest ( s. mansoni ) . however , despite the interbatch variation , it was obvious that the s. mansoni infection drove the metabolite profiles in a different direction to the two protozoan infections and that the batch differences could be largely eliminated by use of o - pls modeling ( data not shown ) . global pls - da scores ( a ) and loadings line plots for discriminatory components 2 and 3 ( b ) from plasma h nmr spectral data showing differentiation of the samples from noninfected control animals ( black ) and the 3 parasite - mouse models ; s. mansoni ( red ) , t. b. brucei ( green ) , and p. berghei ( blue ) . the protozoan infections are clearly differentiated from the noninfected and s. mansoni - infected groups along the first pc , while s. mansoni is differentiated from all other groups in the second pc . ( c ) pca scores plot and ( d ) loadings plots derived from relative cytokine levels from matched animals from each infection model . the recovery of urine and plasma metabolic candidate biomarkers ( 1d standard pulse program ) for s. mansoni via o - pls - da comprised higher levels of urinary phenylacetylglycine and 3-ureidopropanoate , and lower relative levels of 2-oxoisocaproate , hippurate and 2-oxoiadipate , while the plasma profile was characterized by lower concentrations of acetate , glucose and citrate , but relatively higher lipoprotein fractions ( table 1 ) . the spectral pca scores plot of plasma samples from s. mansoni - infected and noninfected mice showed three separate groupings ( figure 4a ) , whereby in addition to clear separation of infected and noninfected animals , a third cluster was evident , consisting of two individual mice that underwent initial infection but where infection failed to establish , as evidenced by zero worm counts upon dissection of these mice . the pls - da model was built using two components ( pls - da : rx = 0.32 , ry = 0.83 , qy = 0.20 ) . ( a ) pls - da scores plot showing plasma metabolite - based differentiation of s. mansoni - infected ( red ) , noninfected control animals ( black ) and mice , which did not establish infection ( possibly demonstrating early expulsion ) as assessed upon dissection ( blue ) , and ( b ) loadings line plot for the first and second pcs . ( c ) pca scores plot of the relative plasma cytokine concentrations , and ( d ) main corresponding loadings showing similar infection - related clustering as the spectral data ; ( e ) bar chart showing plasma cytokine concentrations for s. mansoni - infected mice . values are mean  standard deviation . significance determined by student s t - test ( 2-tailed , unpaired ) p = * < 0.05 , * * < 0.01 , * * * < 0.001 ( since there were only two animals in which infection failed to establish , statistical significance was not determined for the intermediate group ) . note : cytokine concentrations have been scaled to allow visual comparison across cytokines . il-4 scaled 1:2 , il-5 scaled 1:10 , il-12 scaled 1:50 and mkc scaled 1:200 . * separate pls - da models were constructed comparing infected and noninfected control animals for each parasiterodent model , or in the case of p. berghei , pre- and postinfection status . the top 4 discriminating variables ( i.e. , cytokines ) were extracted . for p. berghei infection , il-12 , mkc , ifn- , and il-4 were the 4 most discriminatory cytokines between the pre- and postinfection time points , whereas the t. b. brucei - infected mice were discriminated from the controls at day 28 postinfection mainly by il-1 and il-5 . the main discriminatory cytokines for each pairwise comparison of infection and noninfection status investigated are listed in table 1 . the cytokine distribution of rats infected with f. hepatica was also clearly differentiated from that of controls as shown in the pca scores and loadings plots ( figures 2c , d ) , whereby separation was mainly influenced by increased il-5 levels in infected animals ( rcumx = 0.99 , qx = 0.71 ) . three pcs were used to build the model . similarly to the analysis of the metabolic data , pca analysis of cytokine data of the s. mansoni - infection model indicated the presence of three different groups , one of which corresponded to the two animals in the infected group in which infection failed to establish ( figures 4c , d ) ( rcumx = 0.67 , qx = 0.25 ) . samples from these animals were clearly separated from both uninfected controls and the infected mice , and the mapping position from these anomalous mice was mainly influenced by the plasma level of mkc , which was relatively higher in these two mice than in either the infected or the noninfected control animals ( figures 4c , d ) . to compare plasma levels of the cytokines across the different parasitic infections in the same host species ( i.e. , nmri mouse ) , pca analyses were carried out on the integrated murine data set using two pcs ( rcumx = 0.67 , qx = 0.18 ) . the distribution of the samples in a multivariate space was compared for the three murine infection models based on the cytokine composition ( figure 3c , d ) , as determined via the pca scores and loadings plots ( s. mansoni on day 53 postinfection , t. b. brucei on day 28 postinfection , p. berghei on day 4 postinfection ) . the two protozoan infections demonstrated a similar trend in mapping position in the plots and deviated from the control group in the first pc , but t. b. brucei - infected mice had a wider spread compared to the animals infected with p. berghei , indicating greater variation in response to infection . in contrast , the trematode s. mansoni induced an orthogonal shift in the second pc direction compared to the two protozoan infections , with a marked discrimination from the controls also along the second pc axis . higher levels of il-4 and il-5 dominated the separation between s. mansoni - infected animals from both protozoan infected and noninfected control mice . increased concentrations of il-12 , tnf- , il-10 , ifn- , and il-2 were the main contributors to the discrimination of the mice infected with protozoa . the matrices for both the nmr - characterized metabolite profiles and the cytokine data were statistically related using o - pls analysis in order to obtain information on potentially causal events at the interface of the metabolic and immune systems . the model was constructed using one pls component and one orthogonal component ( o - pls - da : ry = 0.99 , qy = 0.77 ) . there were several statistically significant correlations between the relative cytokine levels and the corresponding plasma spectra of 5 mice preinfection and at days 2 and 4 postinfection in mice infected with p. berghei ( figure 5 ) . a positive correlation was found between ifn- and plasma concentrations of n - acetyl glycoprotein fragments , lactate , creatine , and glycine , whereas ifn- was anticorrelated with glucose . the same correlation script was applied to urine and plasma spectra from each of the parasitehost models . although trends in covariation between metabolites and cytokines were established , the correlations between the relative cytokine levels and metabolites were not significant , which may be due to the relatively small group sizes . o - pls correlation analysis of plasma obtained from p. berghei - infected mice ( n = 5 ) before infection and on days 2 and 4 postinfection correlated with the individual matched relative levels of ifn-. upward - oriented peaks represent positive correlation of metabolites with ifn- concentration and vice versa . each of the assessed parasiterodent models gave rise to a distinct pathological phenotype on the basis of both the metabolic and immunologic profiles . the metabolic fingerprint of each infection was composed of disease - specific features such as relatively increased urinary levels of pipecolic acid due to a p. berghei infection , or increased 4-hydroxyphenylacetic acid as response to a t. b. brucei infection ( table 1 ) . on the other hand , several biomarkers were indicative of more than one infection , suggesting some commonly affected metabolic pathways across infections . interestingly , many of the gut microbial - derived urinary metabolites seem to undergo similar directional changes in response to parasite infection , such as depleted concentrations of hippurate ( t. b. brucei , s. mansoni , and f. hepatica ) , increased trimethylamine ( t. b. brucei and s. mansoni ) and increased phenylacetylglycine ( p. berghei and s. mansoni ) . the association of hippurate , phenylacetylglycine , and the oxidized form of trimethylamine ( tmao ) with gut microbial composition and/or activity has been underlined by studies on germ - free rats and mice with a reduced microbial composition . the involvement of the gut microbiota has been discussed in more detail previously with regard to parasite infections . a common metabolic feature between t. b. brucei and f. hepatica was the elevation in levels of acetylated glycoproteins in plasma ( both n- and o - acetyl glycoproteins ) . acute phase proteins , such as transferrin , -1-antitrypsin or haptoglobin , are secreted by hepatocytes in response to tissue damage and act as inflammatory mediators , and thus , the increased concentrations of plasma acetylated glycoproteins are likely to reflect an inflammatory response . further analysis of s. mansoni - infected , noninfected control mice and the intermediate which group of 2 mice , which underwent experimental infection but which failed to establish adult worms , showed a comapping of the branched chain amino acids ( bcaas ) with the established infection group , whereas the phenylacetylglycine , tyrosine and trimethylamine contributed most to the separation of the  intermediate  group , which may suggest an involvement of the gut microbial environment in susceptibility and resistance to a helminth infection and merits further investigation . comparison of the metabolic profiles of plasma across the three mouse models showed a spatial separation of the helminth infection from the two protozoan infections along the third pc . the corresponding loadings revealed relatively higher levels of lactate , acetate and creatine in the protozoan infections , whereas the ( ch2)n - lipid signals were relatively lower . the differentiation of f. hepatica - infected from noninfected control rats was dominated by increased plasma il-13 , il-5 , and il-1 levels ( table 1 ) . il-13 , characteristically a th2 cytokine , has been shown to directly inhibit pro - inflammatory cytokines in macrophages and plays an important role in minimizing the hepatic inflammation caused by the migrating larvae of f. hepatica . the cytokine response to the extracellular protozoan parasite t. b. brucei consisted of mixed th - subset mechanisms reflecting both cell - mediated and humoral immune mechanisms . the main differential markers of infection with p. berghei , il-12 , and ifn- , are typical th1 cytokines , whereas s. mansoni induced a th2-dominated cytokine pattern ( e.g. , il-4 and il-5 ) . as observed with the metabolic profiles , the s. mansoni - infected mice were further divided into two subgroups corresponding to 5 mice where infection was established as confirmed by worm counts in the liver and mesenteric veins(46 ) and 2 mice where infection failed to establish ( intermediate group ) . the driving cytokine component in the intermediate group was mkc , whereas the main descriptors for established infection were the th2 cytokines il-4 , il-5 , and il-10 . although , the most effective response to a worm infection is believed to be th2-based , s. mansoni seems more efficiently controlled by an initial strong cell - mediated ( th1 ) mechanism , which transfers to a th2-mediated response with progressing course of disease in order to counteract the chronic inflammation caused by the eggs in the liver and other tissues . the mouse functional analogue of human il-8 , mkc , has been shown to play a major role in acute inflammation by recruiting and activating neutrophils to the sites of inflammation , which raises the question if those functions would be crucial to the expulsion of s. mansoni in particular . direct comparison of the plasma cytokine profiles from s. mansoni with those from both protozoan infections in a pca analysis ( figure 3c ) showed clustering of all corresponding control animals and separation of the three infected groups , wherein the 2 groups of mice infected by the protozoan parasites were differentiated from the controls in the first pc and from the trematode infection in the second component . however , the two protozoan infections could also be differentiated from each other with respect to the infection - induced alterations in host cytokine levels . it is widely established that helminth infections induce a strong th2-mediated immune response and here the s. mansoni - induced cytokine expression at a chronic stage time point was consistent with th2-dominated immune mechanisms . the th1 subset of cytokines involved in the response to t. b. brucei and p. berghei , on the other hand , was stronger , as suggested by increased il-12 and ifn- , which are ranked higher in the list of the discriminatory cytokines in the protozoan infections , particularly with respect to p. berghei . many protozoan parasites have intracellular stages where they are cleared by cell - mediated immune mechanisms induced by the th1 cytokine subset . even though t. b. brucei is an extracellular parasite like s. mansoni , its effect on the host immune system is closer to that of p. berghei , which manifests both intra- and extracellular stages . african trypanosome species trigger a strong cellular immune response ( th1-mediated ) in infected mice via their variant surface glycoproteins ( vsgs ) , which cover the parasite coat , and for which constant remodeling contributes to successful evasion from the humoral immunity.(50 ) the relatively lower plasma levels of triglycerides ( vldl / ldl ) and phosphatidylcholine observed following a t. b. brucei infection , with subsequent increase of choline in the plasma of t. b. brucei - infected mice may be associated with the constant remodeling of the glycosylphosphatidylinositol ( gpi ) anchor of t. b. brucei , in order to enable the rapidly changing vsg composition . vsgs also induce secretion of tnf- , which initiates b - cell proliferation and has a direct cytotoxic effect on t. b. brucei.(51 ) to obtain an overview of the global similarities / dissimilarities between the different infection - induced signatures of the 3 parasite - mouse models , comparison of the plasma ( figure 3 ) and urine samples ( data not shown ) was made . both the metabolic and the cytokine data yielded a similar distribution of samples ( figure 3 ) , in that samples of the protozoan infections were differentiated both from the control and from the helminth infection samples . a tight clustering of the cytokine data from all control samples was noted in the scores plot generated from the cytokine data . however , with respect to the metabolic data , batch differences between the control samples in both the urine ( data not shown ) and the plasma data sets were evident . this is unsurprising given that the nmri mouse is an outbred strain and the h nmr metabolic profiles are composed of thousands of signals , which contribute to the model and which all have the potential to exhibit batch variation , whereas the cytokine panel represents a relatively coarse profile consisting of only 9 variables . temporal and batch differences in h nmr spectra of biofluids are well documented and result from the high sensitivity of the technology . however , despite the observed batch - to - batch variation between the 3 parasitemouse models , the parasite - induced alteration in the metabolic signature was stronger and the presence of these batch differences did not impede detection of clear infection - induced responses with further differentiation of the protozoan and helminth infections . analysis of each parasiterodent model uncovered covariation between a range of metabolic and immune features , which may potentially be mechanistically related . in the correlation assessment of plasma h nmr spectral data and cytokine concentrations in the p. berghei data set , direct statistical connections with biological coherence were evident , such as the positive correlations of ifn- , with n - acetyl glycoproteins , whereby both are important pro - inflammatory signaling molecules . on the other hand , the positive correlation of ifn- with lactate and subsequent negative correlation with glucose deriving from the p. berghei - induced increase in the glycolytic activity may simply reflect the fact that these changes dominated the cytokine and metabolic profiles , respectively , and thereby the statistical correlation may result from covariation without necessarily inferring a direct mechanistic link . however , as all other assessed cytokines did not show any significant correlation with either glucose or lactate , a more direct involvement of the immune system and the components of the glycolysis may still be implicated . to explore causality , with respect to the observed covariation , further studies are required . focusing on the direct correlation of metabolic and immune components , this approach will ultimately allow the extraction of directly related events at the interface of metabolic and immune system and provide a mechanistic tool for probing pathology . the approach employed in the current analysis , based on the parallel and integrated assessment of the metabolite and cytokine profiles from experimental infections in rodent models , enabled the extraction of the component of the metabolic signature directly and indirectly associated with immune mechanisms . a deeper understanding of the association between metabolic and immune events in these well - controlled laboratory models could open new avenues in understanding immunoregulatory mechanisms and their functions acting at a global metabolic level across a wide range of parasitic infections and diseases .
parasitic infections cause a myriad of responses in their mammalian hosts , on immune as well as on metabolic level . a multiplex panel of cytokines and metabolites derived from four parasite - rodent models , namely , plasmodium bergheimouse , trypanosoma brucei bruceimouse , schistosoma mansonimouse , and fasciola hepaticarat were statistically coanalyzed . 
 1h nmr spectroscopy and multivariate statistical analysis were used to characterize the urine and plasma metabolite profiles in infected and noninfected animals . 
 each parasite generated a unique metabolic signature in the host . 
 plasma cytokine concentrations were obtained using the  meso scale discovery multi cytokine assay platform . 
 multivariate data integration methods were subsequently used to elucidate the component of the metabolic signature which is associated with inflammation and to determine specific metabolic correlates with parasite - induced changes in plasma cytokine levels . 
 for example , the relative levels of acetyl glycoproteins extracted from the plasma metabolite profile in the p. berghei - infected mice were statistically correlated with ifn- , whereas the same cytokine was anticorrelated with glucose levels . 
 both the metabolic and the cytokine data showed a similar spatial distribution in principal component analysis scores plots constructed for the combined murine data , with samples from all infected animals clustering according to the parasite species and whereby the protozoan infections ( p. berghei and t. b. brucei ) grouped separately from the helminth infection ( s. mansoni ) . for s. mansoni , 
 the main infection - responsive cytokines were il-4 and il-5 , which covaried with lactate , choline , and d-3-hydroxybutyrate . 
 this study demonstrates that the inherently differential immune response to single- and multicellular parasites not only manifests in the cytokine expression , but also consequently imprints on the metabolic signature , and calls for in - depth analysis to further explore direct links between immune features and biochemical pathways .