# Erythropoietin Levels Increase during Cerebral Malaria and Correlate with Heme, Interleukin-10 and Tumor Necrosis Factor-Alpha in India

E. Dalko, N. Tchitchek, Laurent Pays, F. Herbert, P. A. Cazenave, B.

Ravindran, S. Sharma, Serge Nataf, B. Das, S. Pied

► To cite this version:

E. Dalko, N. Tchitchek, Laurent Pays, F. Herbert, P. A. Cazenave, et al., Erythropoietin Levels Increase during Cerebral Malaria and Correlate with Heme, Interleukin-10 and Tumor Necrosis Factor-Alpha in India. PLoS ONE, 2016, 11 (7), pp.e0158420. 10.1371/journal.pone.0158420. hal-01850300

HAL Id: hal-01850300

https://hal.science/hal-01850300

Submitted on 27 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Erythropoietin Levels Increase during Cerebral Malaria and Correlate with Heme, Interleukin-10 and Tumor Necrosis Factor-Alpha in India

Esther Dalko¹ᵃ, Nicolas Tchitchek², Laurent Pays³, Fabien Herbert¹, Pierre-André Cazenave¹ᵇ, Balachandran Ravindran⁴, Shobhona Sharma⁵, Serge Nataf³, Bidyut Das⁶, Sylviane Pied¹*

1 Centre for Infection and Immunity of Lille, INSERM U1019, CNRS UMR 8204, Université Lille Nord de France, Institut Pasteur de Lille, Lille 59019, France, 2 CEA, DSV/iMETI, Immunology of viral infections and autoimmune diseases research unit, UMR1184, IDMIT infrastructure, Fontenay-aux-Roses, France, 3 Lyon 1 University, CarMeN Laboratory, INSERM U-1060, INRA USC-1235, 69921, Oullins, France; Banque de Tissus et de Cellules des Hospices Civils de Lyon, Hôpital Edouard Herriot, Lyon, France, 4 Institute of Life Sciences, Bhubaneswar, Odisha, India, 5 Department of Biological Sciences, Tata Institute of Fundamental Research, Mumbai, Maharashtra 400005, India, 6 SCB Medical College, Cuttack, Odisha 753007, India

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Dalko E, Tchitchek N, Pays L, Herbert F, Cazenave P-A, Ravindran B, et al. (2016) Erythropoietin Levels Increase during Cerebral Malaria and Correlate with Heme, Interleukin-10 and Tumor Necrosis Factor-Alpha in India. PLoS ONE 11 (7): e0158420. doi:10.1371/journal.pone.0158420

**Editor:** Luzia Helena Carvalho, Centro de Pesquisa Rene Rachou/Fundação Oswaldo Cruz (Fiocruz-Minas), BRAZIL

**Received:** December 18, 2015

**Accepted:** June 15, 2016

**Published:** July 21, 2016

**Copyright:** © 2016 Dalko et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information file.

**Funding:** This work was supported by the Indo-French Centre for the Promotion of Advanced Research (IFCPAR), the Associated International Laboratory Systems Immunology and Genetics of Infectious Diseases (SIGID), the Department of Biotechnology from the Ministry of Science and Technology of India (DBT), and intramural funds from Tata Institute of Fundamental Research. The project was supported by the Raman-Charpak award from

## Abstract

Cerebral malaria (CM) caused by *Plasmodium falciparum* parasites often leads to the death of infected patients or to persisting neurological sequelae despite anti-parasitic treatments. Erythropoietin (EPO) was recently suggested as a potential adjunctive treatment for CM. However diverging results were obtained in patients from Sub-Saharan countries infected with *P. falciparum*. In this study, we measured EPO levels in the plasma of well-defined groups of *P. falciparum*-infected patients, from the state of Odisha in India, with mild malaria (MM), CM, or severe non-CM (NCM). EPO levels were then correlated with biological parameters, including parasite biomass, heme, tumor necrosis factor (TNF)-α, interleukin (IL)-10, interferon gamma-induced protein (IP)-10, and monocyte chemoattractant protein (MCP)-1 plasma concentrations by Spearman's rank and multiple correlation analyses. We found a significant increase in EPO levels with malaria severity degree, and more specifically during fatal CM. In addition, EPO levels were also found correlated positively with heme, TNF-α, IL-10, IP-10 and MCP-1 during CM. We also found a significant multivariate correlation between EPO, TNF-α, IL-10, IP-10 MCP-1 and heme, suggesting an association of EPO with a network of immune factors in CM patients. The contradictory levels of circulating EPO reported in CM patients in India when compared to Africa highlights the need for the optimization of adjunctive treatments according to the targeted population.

<footnote>ᵃa Current address: Département des Sciences Biologiques, Université du Québec à Montréal, Montréal (QC) H2X 1Y4, Canada</footnote>

<footnote>ᵇb Current address: Immunologie—Immunopathologie—Immunothérapie, UPMC/CNRS UMR 7211, Paris 75651, France</footnote>

<footnote>* sylviane.pied@pasteur-lille.fr</footnote>
the IFCPAR and by the AAP n10 award from the
College Doctoral Lille Nord de France.

**Competing Interests:** The authors have declared that no competing interests exist.

**Introduction**

Malaria eradication is a worldwide public health priority. Despite tremendous efforts, children and adults still die from this plague every year due to severe disease manifestations. Cerebral malaria (CM) is the deadliest complication, and as reported by the World Health Organization it is most of the time caused by *Plasmodium falciparum* infections [1]. The pathophysiology of CM is complex and far from being completely understood. It involves brain inflammation triggered by the sequestration of parasitized erythrocytes, cell death and infiltrating immune cells in the central nervous system, as well as the systemic production of pro-inflammatory cytokines [2]. In 18% of the cases, patients die despite receiving anti-parasitic treatment, and surviving patients often suffer from neurological sequelae [2,3]. Thus, there is an important need in finding appropriate adjunctive treatments against CM.

In this perspective, administration of erythropoietin (EPO) has recently raised interest in the treatment of neurodegenerative diseases, including CM. EPO, first described as a hemato-poietic hormone, is principally produced by peritubular capillary endothelial cells of the kidney in response to hypoxia. Upon binding to its homodimeric receptor, EPO-R, EPO stimulates the expression of transcription factors such as Bcl-X<sub>L</sub> that further stimulates the maturation of erythroid progenitors [4]. Interestingly, the *in vitro* analysis of human brain cells have shown that astrocytes also constitutively produce low levels of EPO, which vary in response to cyto-kines [5]. In addition, neurons, glial cells, and endothelial cells, express EPO-R, suggesting a plausible localized impact of this hormone on the central nervous system [5,6]. EPO may exert extra-hematopoietic functions in the brain, diminishing the production of pro-inflammatory cytokines and preventing neuronal cell death, which appear to limit brain damage and neuro-logical sequelae [7-9].

Several studies have suggested that EPO plays a pivotal role in CM by exerting protective effects within the central nervous system. Elevated EPO levels were found in plasma of children with severe *P. falciparum* malaria in Mozambique, including severe malarial anemia, but not within the CM group [10]. Furthermore, high levels of endogenous plasma EPO were associ-ated with reduced risk of developing acute neurological sequelae after CM in Kenya [11]. EPO was also shown to efficiently prevent experimental CM during *P. berghei* ANKA infection of C57BL/6 and CBA/J mice [12-14]. It was likely that the beneficial impact of recombinant human EPO (rhEPO) on CM was concurrent to decreased production of the pro-inflammatory cytokines interferon (IFN)-γ and tumor necrosis factor (TNF)-α, and decreased recruitment of T lymphocytes within the brain [12-14]. These results were the basis for conducting a safety clinical trial in Mali, using 1,500 IU/kg rhEPO as adjunctive treatment during three days after admission of the patients, and combined to quinine for CM [15]. However, a very recent study published by Shabani et al., reported elevated plasma EPO levels in CM cases in children with prolonged coma and increased mortality [16]. In addition, considerable uncertainty still remains regarding the use of rhEPO as a neuroprotective strategy for the treatment of other neurodegenerative disease such as acute stroke. In fact, rhEPO given at the admission of the patients was first associated with smaller lesion size and a better functional recovery of the patients compared to the placebo [17], but a later clinical trial with a higher number of patients has reported deleterious outcomes for the survival of acute stroke patients [18].

In light of these discrepancies and the absence of studies concerning populations with differ-
ent age and environmental factors, we have investigated EPO levels in groups of *P. falciparum*-
infected patients with mild malaria (MM), severe non-cerebral malaria (NCM), or CM from
the state of Odisha, India. In 2010, approximately 85% of the Indian population lived in
malaria-endemic areas, representing 27% of malaria cases outside Africa [19]. Unlike in Sub-
Saharan Africa, mainly adults develop severe malaria in India, with very heterogeneous clinical
pictures that are usually not associated with anemia [20–22]. We addressed the relationship between EPO levels and the clinical severity in *P. falciparum*-infected adults. In addition, we examined the association of EPO levels with biological parameters associated to the lysis of red blood cells [hemoglobin, *P. falciparum* histidine-rich protein-2 (PfHRP-2), hemopexin and total heme levels], and to pro- and anti-inflammatory cytokine associated to heme levels and involved in CM pathophysiology.

## Materials and Methods

### Ethics

This study was conducted according to the guidelines set out in the Declaration of Helsinki. The study was approved by the Institutional Human Ethics Committee of SCB Medical College, Cuttack, India, and the Institutional Review Board of all three collaborating institutes: (i) the Institute of Life Sciences, Bhubaneswar, India, (ii) the Tata Institute of Fundamental Research, Mumbai, India, and (iii) Institut Pasteur de Lille, Lille, France. The National Health Office Ethics committee in India approved the study design. All blood samples were collected after obtaining written consent from participants or, in the case of comatose patients, accompanying relatives.

### Study site and participants

This study is a retrospective analysis of 130 patients admitted at the Department of Medicine of the Sriram Chandra Bhanj Medical College & Hospital (SCB Medical College), Cuttack, Odisha in India, between 2008 and 2011. We have already analyzed these subjects in other studies [20,21]. All patients of age $\ge$ 15 years, with a history of fever and with evidence of altered sensorium, jaundice, oliguria, respiratory distress, shock, and/or bleeding diathesis were screened for *Plasmodium* infection by Giemsa blood smears and an immune chromatography test (SD Bio Standard Diagnostics, India). As described earlier, patients who were positive by immune chromatography test but negative by blood smears were subjected to nested PCR [23]. Only *P. falciparum*-infected patients were included; all *P. vivax* co-infections were excluded from the study.

Following confirmation of *P. falciparum* infection, patients were categorized into distinct clinical groups; MM patients only developed fever, and severe malaria (SM) patients were defined as those having at least one laboratory or clinical feature of severe malaria reported by the WHO [24]. Further categorization of SM patients into two groups, CM and NCM, was based on distinct clinical features. CM were defined as patients with a Glasgow Come Score (GCS) $\le$ 10 after exclusion of other causes of encephalopathy such as encephalitis, meningitis and metabolic encephalopathy by biochemical investigations in the cerebrospinal fluid, and NCM patients were SM patients without cerebral involvement. All severe anemia cases, defined by the WHO as hemoglobin levels < 7 g/dL for adults [1], were excluded from the study. Endemic controls (EC) of identical ethnicity and coming from a similar geographical background were enrolled from among patients’ relatives (n = 25), and did not report a history of clinical malaria in the five years preceding our study.

### Blood collection and assessment of biological parameters

Peripheral venous blood (5 mL) was collected on day 0, before treatment, by vein puncture in sterile EDTA tubes. Plasma was obtained by centrifuging blood samples at 4500g for 15 min, and was stored at -80°C until further use. Complete blood count was performed with 2 mL of blood collected in EDTA and using an automated cell counter (Sysmex XT-2000 Hematology Analyzer, Japan). All biological parameters described earlier were measured in peripheral
venous blood including blood sugar, electrolytes, creatinine, transaminases and lactate levels [20,21]. Those with evidence of bacterial (urinary tract infection, sepsis) and viral (hepatitis A, E, B, C) infections were excluded from the study, as well as patients affected with metabolic diseases, including diabetes type 2.

## Quantification of plasma levels of *P. falciparum* histidine rich protein-2

PfHRP-2 levels were measured in plasma by ELISA, and as described previously [25]. Briefly, an anti-PfHRP-2 IgM (MPFM-55A; Santa Cruz Biotechnology) was used as the capture antibody; an anti-PfHRP-2 IgG1 was used as a detection antibody (MPFG-55; Santa Cruz Biotechnology), and an HRP-coupled monoclonal anti-IgG1 antibody (SB77e; Abcam) was used as the secondary antibody for detection. A standard curve was added in duplicate to each plate, and it was obtained through serial dilutions of a plasma sample from a previously described cohort with known parasitemia [26]. Samples out of the range of detection were more or less diluted. The lower limit of detection was 0.5 arbitrary units (AU).

## Determination of EPO, total heme, hemopexin and heme oxygenase (HO)-1 levels in plasma

Plasma levels of EPO, hemopexin and HO-1 were measured by ELISA using, respectively, the Human EPO Platinum ELISA (eBioscience), the Human Hemopexin ELISA kit (Abcam), and the Human HO-1 ELISA kit (Enzo Life Sciences), and according to the manufacturer's instructions. Total heme levels were measured in plasma using the QuantiChrom Heme Assay kit (BioAssay Systems). The colorimetric reaction was measured at 405 nm.

## Assay of cytokines

TNF-α, interleukin (IL)-10, interferon gamma-induced protein (IP)-10, and monocyte chemoattractant protein (MCP)-1 levels were measured in plasma using the MILLIPLEX MAP multiplex assay kit (Millipore), as described previously [20].

## Statistical analysis

Statistical analyses of erythropoietin levels (EPO) presented in Fig 1 have been performed using GraphPad Prism Software (version 7) and R [27]. Data have been tested for variance equality using the Levene test, and for non-normality using D'Agostino-Pearson omnibus normality tests (p-values < 0.0001). The Kruskal-Wallis test showed that the distribution of ranks was significantly different in some groups (p-value < 0.0001). Sample groups were then compared using the non-parametric Mann-Whitney test and significant differences were defined as comparisons with a p-value below 0.05.

The correlation between EPO plasma levels and biological parameters were determined using the Spearman coefficient of correlation. Multiple correlations analyses have been calculated using linear fitting models. Computations of the linear models have been performed using the 'lm' function of the native 'stat' package in R [27]. For each regression model used in this paper, the estimated value, standard error, t-value and p-value are provided for each coefficient. Moreover for each regression model, the Spearman coefficient of correlation between data and fitted data is provided with its confident interval and the associated p-value. A p-value < 0.05 was considered as significant.

The GeneMANIA Cytoscape plugin (v3.4.0) was used to generate 2 distinct co-expression networks from the clinical, biological and cytokine parameters obtained in severe NCM and CM patients [28]. The severe NCM network was composed of 532 interactions and the CM
<figure><img src="image_4.png" /><figcaption>**Fig 1. Erythropoietin levels increase during cerebral malaria and with severity of the infection. (A)**</figcaption></figure>

Erythropoietin (EPO) levels were measured in the plasma of endemic controls (EC), mild malaria (MM) and severe malaria (SM) patients, and (B) in the SM sub-groups of non-cerebral malaria (NCM) and cerebral malaria (CM). (C) EPO levels were compared between CM patients who survived and those who did not. Significant differences with the EC group or between selected groups are indicated as the following: \*p ≤ 0.05; \*\*p ≤ 0.01; \*\*\*p ≤ 0.001.

doi:10.1371/journal.pone.0158420.g001

network of 628 interactions. A supervised clustering was then performed to identify the top ten parameters that co-regulated with EPO levels. For presentation purposes, nodes with the highest level of connectivity were assigned a different shape using Cytoscape (v3.2.1), an open-source network visualization software [29].

## Results

### Baseline characteristics

A total of 130 patients were included in this study: 28 MM, 45 NCM and 57 CM patients. The strong male bias is in agreement with studies previously reported in India [22,30]. There were no differences in patients' age, and hemoglobin levels between MM, NCM and CM groups (Table 1). PfHRP-2 levels were shown to be good indicators of total parasite biomass [31], and they were significantly higher in the plasma of NCM and CM compared to MM patients (Table 1).

### EPO levels increase during CM

To investigate if EPO levels were associated with a "protective phenotype" against CM in Odisha, EPO levels were determined in the plasma of *P. falciparum*-infected patients and EC. To

<table><caption>**Table 1. Baseline characteristics.**</caption><thead><tr><th></th><th>MM</th><th>NCM</th><th>CM</th><th>All</th></tr></thead><tbody><tr><td>No.</td><td>28</td><td>45</td><td>57</td><td>130</td></tr><tr><td>Age (year) #</td><td>30 (24–35)</td><td>28 (23–41)</td><td>30 (20–43)</td><td>30 (22–41)</td></tr><tr><td>% Male</td><td>60.7</td><td>77.8</td><td>86.0</td><td>77.7</td></tr><tr><td>Hb (g/dL) #</td><td>10.6 (9.9–12.3)</td><td>10.0 (9.0–12.0)</td><td>10.0 (8.6–11.3)</td><td>10.2 (9.0–12.0)</td></tr><tr><td>PfHRP-2 (AU) #</td><td>11.9 (1.1–25.3)</td><td>22.0 (8.3–38.5)*</td><td>29.7 (21.2–47.8)***</td><td>23.2 (9.7–40.8)</td></tr></tbody></table>

\# Values are medians (interquartile ranges).

\* MM vs. NCM: *p*-value < 0.05.

\*\* MM vs. CM: *p*-value < 0.001.

doi:10.1371/journal.pone.0158420.t001
eliminate a possible bias resulting from low hemoglobin levels, patients with severe anemia (Hb < 7 g/dL) were excluded from the study [1]. Surprisingly, we found that EPO levels were significantly higher in the plasma of SM patients compared to EC, and MM (Fig 1A). More specifically, only CM patients had significantly higher plasma EPO compared to MM patients (Fig 1B, left panel). EPO levels were not significantly different between MM and NCM patients, as well as between CM and NCM. In our groups of patients, only those developing neurological outcomes died, and we show a 6.7-fold increase in the median levels of EPO among the group of CM patients who died (n = 16) compared to CM survivors (n = 41) (Fig 1B, right panel); hemoglobin levels were not significantly different between these two groups (p-value = 0.3559). In addition, EPO levels were significantly higher in plasma of CM non-survivors compared to NCM and MM patients (Fig 1B, right panel), suggesting that high levels of EPO are associated with lethality in Odisha. It is to note that age, sex, PfHRP-2 levels and hemopexin levels were independently associated with mortality.

## EPO levels correlated with biological parameters

During blood stage malaria, EPO production by the kidney increases to counteract the loss of red blood cells. Considering that in our groups of malaria-infected patients hemoglobin levels correlated significantly and negatively with PfHRP-2 plasma levels (R = -0.217; *p*-value = 0.0013), we investigated correlations between EPO and hemoglobin, PfHRP-2, total heme and hemopexin plasma levels to further characterize the conditions associated to enhanced EPO levels. When considering all *P. falciparum*-infected patients, EPO correlated significantly with hemoglobin and PfHRP-2 plasma levels, with a negative and a positive factor, respectively. In addition, EPO plasma levels correlated with hemoglobin for MM and CM patients, but not for the NCM ones. A correlation between EPO and PfHRP-2 levels was found for both NCM and CM (Table 2). We have previously reported that total heme and hemopexin levels are good indicators for free heme levels [21]. Hemopexin levels correlated negatively with EPO when analyzing all *P. falciparum*-infected patients, as well as both SM groups, the NCM and the CM group. EPO correlated positively with total heme only for CM patients (Table 2). As shown by multiple correlation analyses, EPO levels correlated significantly to both sets of data: [hemoglobin + PfHRP-2 + hemopexin + total heme] (R = 0.50; CI (0.359–0.618); *p*-value = 0.0001), and [hemopexin + total heme] (R = 0.45; CI (0.302–0.577); *p*-value = 0.0002).

<table><caption>**Table 2.** Correlation of erythropoietin plasma levels, with biological parameters.</caption><thead><tr><th></th><th>MM</th><th>NCM</th><th>CM</th><th>All</th></tr></thead><tbody><tr><td><strong>Hemoglobin</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td>-0.42</td><td>-0.25</td><td>-0.35</td><td>-0.34</td></tr><tr><td><em>p</em>-value</td><td>0.0421</td><td>0.1070</td><td>0.0118</td><td>0.0002</td></tr><tr><td><strong>PfHRP-2</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td>0.32</td><td>0.38</td><td>0.27</td><td>0.41</td></tr><tr><td><em>p</em>-value</td><td>0.0922</td><td>0.0127</td><td>0.0419</td><td>&lt; 0.0001</td></tr><tr><td><strong>Hemopexin</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td>-0.22</td><td>-0.41</td><td>-0.41</td><td>-0.44</td></tr><tr><td><em>p</em>-value</td><td>0.2534</td><td>0.0069</td><td>0.0024</td><td>&lt; 0.0001</td></tr><tr><td><strong>Total heme</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td>-0.04</td><td>0.03</td><td>0.30</td><td>0.17</td></tr><tr><td><em>p</em>-value</td><td>0.8500</td><td>0.8532</td><td>0.0209</td><td>0.0526</td></tr></tbody></table>

R: Spearman coefficient of correlation.

doi:10.1371/journal.pone.0158420.t002
<table><caption>Table 3. Multiple correlation analyses of [EPO] with [hemoglobin + PfHRP-2 + hemopexin + total heme].</caption><thead><tr><th></th><th>Estimate</th><th>Std. Error</th><th>t-value</th><th>Pr (&gt;|t|)</th></tr></thead><tbody><tr><td>(Intercept)</td><td>456.3511</td><td>126.2689</td><td>3.614</td><td>0.0005</td></tr><tr><td>Hemoglobin</td><td>-25.6799</td><td>11.6281</td><td>-2.208</td><td>0.0292</td></tr><tr><td>PfHRP-2</td><td>0.8461</td><td>0.7130</td><td>1.187</td><td>0.2378</td></tr><tr><td>Hemopexin</td><td>-0.3550</td><td>0.1464</td><td>-2.425</td><td>0.0169</td></tr><tr><td>Total heme</td><td>0.0103</td><td>0.3874</td><td>0.027</td><td>0.9788</td></tr></tbody></table>

The estimated coefficient (Estimate) corresponds to the fitted variable value, the standard error (Std. Error) indicates the variability of the estimated coefficient, the t-value (t-value) measures if the estimated coefficient for this variable is meaningful for the model, and the variable p-value (Pr (>|t|)) indicates the probability that the variable is not relevant in the model. R = 0.50; CI (0.359–0.618); p-value = 0.0001.

doi:10.1371/journal.pone.0158420.t003

<table><caption>Table 4. Multiple correlation analyses of [EPO] with [hemopexin + total heme].</caption><thead><tr><th></th><th>Estimate</th><th>Std. Error</th><th>t-value</th><th>Pr (&gt;|t|)</th></tr></thead><tbody><tr><td>(Intercept)</td><td>227.5893</td><td>41.9219</td><td>5.429</td><td>2.77e-07</td></tr><tr><td>Hemopexin</td><td>-0.4708</td><td>0.1268</td><td>-3.713</td><td>0.0003</td></tr><tr><td>Total heme</td><td>0.1330</td><td>0.3660</td><td>0.363</td><td>0.7169</td></tr></tbody></table>

The estimated coefficient (Estimate) corresponds to the fitted variable value, the standard error (Std. Error) indicates the variability of the estimated coefficient, the t-value measures whether or not the estimated coefficient for this variable is meaningful for the model, and the variable p-value (Pr (>|t|)) indicates the probability that the variable is not relevant in the model. R = 0.45 CI (0.302–0.577); p-value = 0.0001.

doi:10.1371/journal.pone.0158420.t004

value = 0.0001) (Tables 3 and 4). According to our analyses, mainly hemoglobin and hemopexin levels are likely to associate with EPO levels as suggested by their low Pr (>|t|), and with a negative factor of the estimated coefficient.

Given our findings on the heme/hemopexin balance that influences cytokine profiles during severe malaria [32], we first investigated a possible correlation between EPO, TNF-α, IL-10, IP-10 and MCP-1 plasma levels. These cytokines have been shown to contribute to the pathophysiology of human CM in India [20,22]. All four cytokines were positively correlated to EPO levels during malaria, as well as to CM cases. EPO levels were not correlated to IP-10 and MCP-1 levels during MM (Table 5). Interestingly, the set of data [EPO + hemopexin + total heme] correlated significantly to [TNF-α + IL-10 + IP-10 + MCP-1] during infection, and more specifically for the CM, and NCM groups, but not for MM patients (Tables 6–9). During *P. falciparum* infections, IL-10 and MCP-1 associated positively to [EPO + hemopexin + total heme] as suggested by their respective estimated coefficient 0.0616 and 0.0602 that associated to a low variable p-value (Table 9). IP-10 was also shown to contribute to the statistical model during CM, but associated to [EPO + hemopexin + total heme] with a negative factor as indicated by the estimated coefficient of -0.0198 (Table 8).

Then, to complement our analyses, we used GeneMANIA, a plugin of the Cytoscape software that allows co-expression or co-regulation links to be identified between measurable objects [28]. This bioinformatics tool, classically used to generate co-expression networks from microarray data, was adapted here to the correlation analysis of quantifiable clinical, biological parameters and cytokine production correlating with disease severity in our *P. falciparum* infected Odisha patients. In particular, we sought to identify cytokines and bio-clinical parameters that correlated with EPO levels during CM but did not correlated with EPO levels in NCM malaria patients. In both groups, we found that the 10 parameters that most tightly correlated with EPO levels included IL-10 and MCP-1 (Fig 2), which confirms our multiple correlation
<table><caption>Table 5. Correlation of EPO plasma levels, with cytokines.</caption><thead><tr><th></th><th>MM</th><th>NCM</th><th>CM</th><th>All</th></tr></thead><tbody><tr><td><strong>TNF-α</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td>0.42</td><td>0.13</td><td>0.42</td><td>0.32</td></tr><tr><td><em>p</em>-value</td><td>0.0355</td><td>0.4046</td><td>0.0021</td><td>0.0005</td></tr><tr><td><strong>IL-10</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td>0.40</td><td>0.57</td><td>0.51</td><td>0.52</td></tr><tr><td><em>p</em>-value</td><td>0.0479</td><td>0.0003</td><td>0.0002</td><td>&lt; 0.0001</td></tr><tr><td><strong>IP-10</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td>0.10</td><td>0.64</td><td>0.35</td><td>0.49</td></tr><tr><td><em>p</em>-value</td><td>0.6117</td><td>&lt; 0.0001</td><td>0.0111</td><td>&lt; 0.0001</td></tr><tr><td><strong>MCP-1</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>R</td><td>0.05</td><td>0.58</td><td>0.50</td><td>0.50</td></tr><tr><td><em>p</em>-value</td><td>0.8042</td><td>0.0002</td><td>0.0002</td><td>&lt; 0.0001</td></tr></tbody></table>

R: correlation coefficient.

doi:10.1371/journal.pone.0158420.t005

<table><caption>Table 6. Multiple correlation analyses of [EPO + hemopexin + total heme] with [TNF-α + IL-10 + IP-10 + MCP-1] during MM.</caption><thead><tr><th></th><th>Estimate</th><th>Std. Error</th><th>t-value</th><th>Pr (&gt;|t|)</th></tr></thead><tbody><tr><td>(Intercept)</td><td>345.1181</td><td>60.3095</td><td>5.722</td><td>1.11e-05</td></tr><tr><td>TNF-α</td><td>0.5480</td><td>0.1463</td><td>3.747</td><td>0.0012</td></tr><tr><td>IL-10</td><td>0.0788</td><td>0.0641</td><td>1.229</td><td>0.2326</td></tr><tr><td>IP-10</td><td>0.0307</td><td>0.0164</td><td>1.875</td><td>0.0748</td></tr><tr><td>MCP-1</td><td>-0.2493</td><td>0.1644</td><td>-1.516</td><td>0.1444</td></tr></tbody></table>

The estimated coefficient (Estimate) corresponds to the fitted variable value, the standard error (Std. Error) indicates the variability of the estimated coefficient, the *t*-value measures whether or not the estimated coefficient for this variable is meaningful for the model, and the variable *p*-value (Pr (>|t|)) indicates the probability that the variable is not relevant in the model. R = 0.33 CI (-0.049–0.626); *p*-value = 0.1054.

doi:10.1371/journal.pone.0158420.t006

analyses. However, it was striking to observe that EPO levels closely correlated with HO-1 levels in NCM patients (Fig 2A), while it correlated with heme during CM (Fig 2B).

## Discussion

In the light of divergences reported in the literature on circulating levels of EPO during human CM, we investigated EPO levels in a population of patients with distinct disease manifestations,

<table><caption>Table 7. Multiple correlation analyses of [EPO + hemopexin + total heme] with [TNF-α + IL-10 + IP-10 + MCP-1] during NCM.</caption><thead><tr><th></th><th>Estimate</th><th>Std. Error</th><th>t-value</th><th>Pr (&gt;|t|)</th></tr></thead><tbody><tr><td>(Intercept)</td><td>203.0414</td><td>82.0068</td><td>2.476</td><td>0.0180</td></tr><tr><td>TNF-α</td><td>-0.0392</td><td>0.0977</td><td>-0.401</td><td>0.6910</td></tr><tr><td>IL-10</td><td>0.0532</td><td>0.0355</td><td>1.500</td><td>0.1420</td></tr><tr><td>IP-10</td><td>0.0132</td><td>0.0126</td><td>1.050</td><td>0.3010</td></tr><tr><td>MCP-1</td><td>0.0364</td><td>0.0378</td><td>0.962</td><td>0.3420</td></tr></tbody></table>

The estimated coefficient (Estimate) corresponds to the fitted variable value, the standard error (Std. Error) indicates the variability of the estimated coefficient, the *t*-value measures whether or not the estimated coefficient for this variable is meaningful for the model, and the variable *p*-value (Pr (>|t|)) indicates the probability that the variable is not relevant in the model. R = 0.50 CI (0.242–0.691); *p*-value = 0.0008.

doi:10.1371/journal.pone.0158420.t007
<table><caption>**Table 8.** Multiple correlation analyses of [EPO + hemopexin + total heme] with [TNF-α + IL-10 + IP-10 + MCP-1] during CM.</caption><thead><tr><th></th><th>Estimate</th><th>Std. Error</th><th>t-value</th><th>Pr (&gt;|t|)</th></tr></thead><tbody><tr><td>(Intercept)</td><td>330.9000</td><td>60.8000</td><td>5.442</td><td>1.59e-06</td></tr><tr><td>TNF-α</td><td>-0.0001</td><td>0.1243</td><td>-0.001</td><td>0.9991</td></tr><tr><td>IL-10</td><td>0.0856</td><td>0.0296</td><td>2.894</td><td>0.0056</td></tr><tr><td>IP-10</td><td>-0.0198</td><td>0.0097</td><td>-2.038</td><td>0.0468</td></tr><tr><td>MCP-1</td><td>0.1017</td><td>0.0318</td><td>3.200</td><td>0.0024</td></tr></tbody></table>

The estimated coefficient (Estimate) corresponds to the fitted variable value, the standard error (Std. Error) indicates the variability of the estimated coefficient, the *t*-value measures whether or not the estimated coefficient for this variable is meaningful for the model, and the variable *p*-value (Pr (>|t|)) indicates the probability that the variable is not relevant in the model. R = 0.49 CI (0.264–0.665); *p*-value = 0.0002.

doi:10.1371/journal.pone.0158420.t008

<table><caption>**Table 9.** Multiple correlation analyses of [EPO + hemopexin + total heme] with [TNF-α + IL-10 + IP-10 + MCP-1] for all *P. falciparum*-infected patients.</caption><thead><tr><th></th><th>Estimate</th><th>Std. Error</th><th>t-value</th><th>Pr (&gt;|t|)</th></tr></thead><tbody><tr><td>(Intercept)</td><td>320.3575</td><td>41.2283</td><td>7.770</td><td>3.2e-12</td></tr><tr><td>TNF-α</td><td>0.0278</td><td>0.0691</td><td>0.403</td><td>0.6879</td></tr><tr><td>IL-10</td><td>0.0616</td><td>0.0224</td><td>2.750</td><td>0.0069</td></tr><tr><td>IP-10</td><td>-0.0067</td><td>0.0071</td><td>-0.941</td><td>0.3487</td></tr><tr><td>MCP-1</td><td>0.0602</td><td>0.0252</td><td>2.390</td><td>0.0185</td></tr></tbody></table>

The estimated coefficient (Estimate) corresponds to the fitted variable value, the standard error (Std. Error) indicates the variability of the estimated coefficient, the *t*-value measures whether or not the estimated coefficient for this variable is meaningful for the model, and the variable *p*-value (Pr (>|t|)) indicates the probability that the variable is not relevant in the model. R = 0.28 CI (0.049–0.466); *p*-value = 0.0019.

doi:10.1371/journal.pone.0158420.t009

epidemiology, environmental parameters and genetic background than those described earlier in the literature [10,11,16]. We found an increase in EPO levels in the plasma of *P. falciparum*-infected patients from Odisha with severe malaria, and particularly with fatal CM. It is noteworthy that a positive correlation of EPO levels with prolonged coma and survival was also reported earlier in Uganda, but only in patients with less than 8 g/dL hemoglobin [16]. As the destruction of red blood cells and the inhibition of erythropoiesis are known to trigger the production of EPO in response to hypoxia, it is conceivable that increased EPO levels observed in patients from Uganda were attributable to severe anemia and not to CM [33]. Moreover, our NCM and CM patients did not exhibit the same response to EPO, and we observed a correlation between HO-1 and EPO levels during NCM but not CM. We have previously reported that increased MIP-1α levels discriminate NCM from CM in Odisha [20], while MIP-1α was shown by others to inhibit the formation of erythroid progenitors [34]. Therefore, a bystander effect of MIP-1α on erythropoiesis may partly explain the lack of correlation between EPO and hemoglobin levels during NCM, reinforcing the idea that NCM and CM are two different disease phenotypes governed by distinct pathological mechanisms.

Previous studies suggest that EPO production increases with age, and within the normal range of hemoglobin [16,35]. Within our cohort, EPO levels were not affected by age, neither in control nor in infected patients in Odisha (S1 Fig). As most of the reported studies were conducted in children, further analysis in younger Odissi individuals would be necessary to better explain the contradictory results between the different populations studied. One can note that the discrepancy due to age may also reflect the level of clinical immunity, which is known to develop with time under constant exposition to the parasite according to the endemicity of *P*
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /><figcaption>**Fig 2.** Sub-networks showing the ten parameters most related to plasma level of EPO during malaria.
Sub-networks obtained for (A) the NCM group and the (B) CM group. Nodes with the highest number of connections with other parameters are represented by hexagons. For each parameter, the size of the node is associated with its level of correlation with EPO levels. Interactions between parameters are represented by edges (lines), thicker lines representing a higher confidence level for a given interaction.</figcaption></figure>

doi:10.1371/journal.pone.0158420.g002
*falciparum* [36]. While Kenya and Mozambique are hyperendemic for *P. falciparum*, Odisha is mesoendemic with a transmission once a year peaking during the rainy season.

Other than age, distinct environmental determinants such as the genetic background and hematological parameters may also influence the pathophysiology of CM, and consequently the role of EPO in the outcome of malaria. In fact, molecules released during parasite multiplication within red blood cells, such as heme, may contribute to increased EPO levels during malaria [32]. We have previously described a rise in heme levels with malaria severity and, particularly during CM, in the same group of patients [21]. In the same line, Spearman's and multivariate analyses have shown that EPO, total heme and hemopexin levels were significantly correlated during CM. Given our findings, we propose that the release of these endogenous molecules up-regulates EPO levels. Nevertheless, the correlation observed between HO-1 and EPO levels during NCM, but not CM, suggests a heme-dependent amplification loop of EPO that may be linked to the degradation of heme by HO-1. In addition, the significant multivariate correlation of EPO, hemopexin and total heme with TNF-α, IL-10, IP-10 and MCP-1 suggests the implication of EPO in the complex network of factors of the immune system influenced by heme during severe malaria. Free heme is more likely to induce a pro-inflammatory response during acute malarial infection, and, similarly to EPO, we have previously reported a positive correlation between heme and TNF-α, IL-10, IP-10 and MCP-1 during *P. falciparum* malaria in India [21]. Unlike IL-10 that is protective against CM, TNF-α, IP-10 and MCP-1 were shown to be important mediators in the induction of rodent and human CM [20,37–39]. However, it is unclear whether high levels of EPO are contributing to, or are insufficient in preventing inflammation associated to CM. In fact, EPO is likely to affect the production of cytokines by immune cells but the overall response is still not clear. For instance, EPO mediates the pro-inflammatory response of macrophages and dendritic cells stimulated with lipopolysaccharide *in vitro*, as suggested by their increased production of IL-12 and decreased production of IL-10 [40,41]. However, Nairz et al. have reported a decreased production of TNF-α by lipopolysaccharide-stimulated primary macrophages that were pretreated with EPO [42]. In the line with these observations, exogenous EPO was shown to decrease the expression of pro-inflammatory cytokines in the brain of C57BL/6 and CBA/J mice infected with *P. berghei* ANKA [12,13].

Rodent models are essential for understanding the pathophysiology of malaria. However, they are insufficient as proof of principle for the development of human therapeutic targets. Our analysis suggests that, in contrast to what was reported earlier from Sub-Saharan Africa, EPO levels increased in CM patients from Odisha. These observations suggest differences in the pathogenesis of malaria between Sub-Saharan African and South Asian populations and highlight the need to better evaluate malaria-related parameters in different geographical and epidemiological contexts for the optimization of anti-malarial and adjunctive treatments.

**Supporting Information**

S1 Fig. Spearman rank's correlation analyzes of plasma EPO levels with age in different clinical groups of malaria-infected patients.

(PDF)

**Acknowledgments**

We thank Dr. Tatiana Scorza (Département des Sciences Biologiques, Université du Québec à Montréal) for kindly reviewing the manuscript.
## Author Contributions

Conceived and designed the experiments: ED SP. Performed the experiments: ED. Analyzed the data: ED NT SP SS PAC BD LP SN FH. Contributed reagents/materials/analysis tools: BD. Wrote the paper: ED NT SP SN. Approved and corrected the manuscript: PAC BD SS BR SN FH.

## References

1. World Health Organization. Management of severe malaria [Internet]. 2012. Available: http://apps.who.int/iris/bitstream/10665/79317/1/9789241548526_eng.pdf. Accessed November 2014.

2. Idro R, Jenkins NE, Newton CRJC. Pathogenesis, clinical features, and neurological outcome of cerebral malaria. Lancet Neurol. 2005; 4: 827–840. PMID: 16297841

3. Idro R, Kakooza-Mwesige A, Balyejjussa S, Mirembe G, Mugasha C, Tugumisirize J, et al. Severe neurological sequelae and behaviour problems after cerebral malaria in Ugandan children. BMC Res Notes. 2010; 3: 104. doi: 10.1186/1756-0500-3-104 PMID: 20398391

4. Dolznig H, Habermann B, Stangl K, Deiner EM, Moriggl R, Beug H, et al. Apoptosis protection by the Epo target Bcl-XL allows factor-independent differentiation of primary erythroblasts. Curr Biol. 2002; 12: 1076–1085. PMID: 12121614

5. Nagai A, Nakagawa E, Choi HB, Hatori K, Kobayashi S, Kim SU. Erythropoietin and erythropoietin receptors in human CNS neurons, astrocytes, microglia, and oligodendrocytes grown in culture. J Neuropathol Exp Neurol. 2001; 60: 386–392. PMID: 11305874

6. Medana IM, Day NP, Hien TT, White NJ, Turner GD. Erythropoietin and its receptors in the brainstem of adults with fatal falciparum malaria. Malar J. 2009; 8: 261. doi: 10.1186/1475-2875-8-261 PMID: 19930602

7. Casals-Pascual C, Idro R, Picot S, Roberts DJ, Newton CRJC. Can erythropoietin be used to prevent brain damage in cerebral malaria? Trends in Parasitology. 2009; 25: 30–36. doi: 10.1016/j.pt.2008.10.002 PMID: 19008152

8. Bond WS, Rex TS. Evidence that erythropoietin modulates neuroinflammation through differential action on neurons, astrocytes, and microglia. Frontiers in Immunology. 2014; 5: 523. doi: 10.3389/fimmu.2014.00523 PMID: 25374571

9. Subirós N, Del Barco DG, Coro-Antich RM. Erythropoietin: still on the neuroprotection road. Ther Adv Neurol Disord. 2012; 5: 161–73. doi: 10.1177/1756285611434926 PMID: 22590480

10. Díez-Padrisa N, Aguilar R, Machevo S, Morais L, Nharnossa T, O'Callaghan-Gordo C, et al. Erythropoietin levels are not independently associated with malaria-attributable severe disease in Mozambican children. PLoS One. 2011; 6: e24090. doi: 10.1371/journal.pone.0024090 PMID: 21912616

11. Casals-Pascual C, Idro R, Gicheru N, Gwer S, Kitsao B, Gitau E, et al. High levels of erythropoietin are associated with protection against neurological sequelae in African children with cerebral malaria. Proc Natl Acad Sci U S A. 2008; 105: 2634–2639. doi: 10.1073/pnas.0709715105 PMID: 18263734

12. Kaiser K, Texier A, Ferrandiz J, Buguet A, Meiller A, Latour C, et al. Recombinant human erythropoietin prevents the death of mice during cerebral malaria. J Infect Dis. 2006; 193: 987–995. PMID: 16518761

13. Wiese L, Hempel C, Penkowa M, Kirkby N, Kurtzhals JAL. Recombinant human erythropoietin increases survival and reduces neuronal apoptosis in a murine model of cerebral malaria. Malar J. 2008; 7: 3. doi: 10.1186/1475-2875-7-3 PMID: 18179698

14. Wei X, Li Y, Sun X, Zhu X, Feng Y, Liu J, et al. Erythropoietin protects against murine cerebral malaria through actions on host cellular immunity. Infect Immun. 2014; 82: 165–173. doi: 10.1128/IAI.00929-13 PMID: 24126529

15. Picot S, Bienvenu A-L, Konate S, Sissoko S, Barry A, Diarra E, et al. Safety of epoietin beta-quinine drug combination in children with cerebral malaria in Mali. Malar J. 2009; 8: 169. doi: 10.1186/1475-2875-8-169 PMID: 19630971

16. Shabani E, Opoka RO, Idro R, Schmidt R, Park GS, Bangirana P, et al. High plasma erythropoietin levels are associated with prolonged coma duration and increased mortality in children with cerebral malaria. Clin Infect Dis. 2015; 60: 27–35. doi: 10.1093/cid/ciu735 PMID: 25228706

17. Ehrenreich H, Hasselblatt M, Dembowski C, Cepek L, Lewczuk P, Stiefel M, et al. Erythropoietin therapy for acute stroke is both safe and beneficial. Mol Med. 2002; 8: 495–505. PMID: 12435860

18. Ehrenreich H, Weissenborn K, Prange H, Schneider D, Weimar C, Wartenberg K, et al. Recombinant human erythropoietin in the treatment of acute ischemic stroke. Stroke. 2009; 40: e647–656. doi: 10.1161/STROKEAHA.109.564872 PMID: 19834012
19. World Health Organization. World Malaria Report. 2010. Available: http://www.who.int/malaria/world_malaria_report_2010/worldmalariareport2010.pdf?ua=1. Accessed November 2014.

20. Herbert F, Tchitchek N, Bansal D, Jacques J, Pathak S, Bécavin C, et al. Evidence of IL-17, IP-10, and IL-10 involvement in multiple-organ dysfunction and IL-17 pathway in acute renal failure associated to *Plasmodium falciparum* malaria. J Transl Med. 2015; 13: 369. doi: 10.1186/s12967-015-0731-6 PMID: 26602091

21. Dalko E, Das B, Herbert F, Fesel C, Pathak S, Tripathy R, et al. Multifaceted role of heme during severe *Plasmodium falciparum* infections in India. Infect Immun. 2015; 83: 3793-3799. doi: 10.1128/IAI.00531-15 PMID: 26169278

22. Prakash D, Fesel C, Jain R, Cazenave P-A, Mishra GC, Pied S. Clusters of cytokines determine malaria severity in *Plasmodium falciparum*-infected patients from endemic areas of Central India. J Infect Dis. 2006; 194: 198-207. PMID: 16779726

23. Panda AK, Panda SK, Sahu AN, Tripathy R, Ravindran B, Das BK. Association of ABO blood group with severe *falciparum* malaria in adults: case control study and meta-analysis. Malar J. 2011; 10: 309. doi: 10.1186/1475-2875-10-309 PMID: 22011404

24. World Health Organization. Guidelines for the treatment of malaria. 2010. Available: http://whqlibdoc.who.int/publications/2010/9789241547925_eng.pdf?ua=1. Accessed November 2014.

25. Noedl H, Bronnert J, Yingyuen K, AttImayr B, Kollaritsch H, Fukuda M. Simple histidine-rich protein 2 double-site sandwich enzyme-linked immunosorbent assay for use in malaria drug sensitivity testing. Antimicrob Agents Chemother. 2005; 49: 3575-3577. PMID: 16048989

26. Bansal D, Herbert F, Lim P, Deshpande P, Bécavin C, Guiyedi V, et al. IgG autoantibody to brain beta tubulin III associated with cytokine cluster-II discriminate cerebral malaria in central India. PLoS One. 2009; 4: e8245. doi: 10.1371/journal.pone.0008245 PMID: 20011600

27. R Development Core Team (2008). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN: 3-900051-07-0. Available: http://www.R-project.org.

28. Montojo J, Zuberi K, Rodriguez H, Kazi F, Wright G, Donaldson SL, et al. GeneMANIA cytoscape plugin: Fast gene function predictions on the desktop. Bioinformatics. 2010; 26: 2927-2928. doi: 10.1093/bioinformatics/btq562 PMID: 20926419

29. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: A software Environment for integrated models of biomolecular interaction networks. Genome Res. 2003; 13: 2498-2504. PMID: 14597658

30. Pathak S, Rege M, Gogtay NJ, Aigal U, Sharma SK, Valecha N, et al. Age-dependent sex bias in clinical malarial disease in hypoendemic regions. PLoS One. 2012; 7: e35592. doi: 10.1371/journal.pone.0035592 PMID: 22558172

31. Dondorp AM, Desakorn V, Pongtavornpinyo W, Sahassananda D, Silamut K, Chotivanich K, et al. Estimation of the total parasite biomass in acute *falciparum* malaria from plasma PfHRP2. PLoS Med. 2005; 2: e0788-0797.

32. Dalko E, Gaudreault V, Sanchez-Dardon J, Moreau R, Scorza T. Preconditioning with hemin decreases *Plasmodium chabaudi adami* parasitemia and inhibits erythropoiesis in BALB/c mice. PLoS One. 2013; 8: e54744. doi: 10.1371/journal.pone.0054744 PMID: 23358441

33. Fandrey J. Oxygen-dependent and tissue-specific regulation of erythropoietin gene expression. Am J Physiol Regul Integr Comp Physiol. 2004; 286: R977-R988. PMID: 15142852

34. Su S, Mukaida N, Wang J, Zhang Y, Takami A, Nakao S, et al. Inhibition of immature erythroid progenitor cell proliferation by macrophage inflammatory protein-1 alpha by interacting mainly with a C-C chemokine receptor, CCR1. Blood. 1997; 90: 605-11. PMID: 9226160

35. Eckardt KU, Ratcliffe PJ, Tan CC, Bauer C, Kurtz A. Age-dependent expression of the erythropoietin gene in rat liver and kidneys. J Clin Invest. 1992; 89: 753-760. PMID: 1541670

36. Doolan DL, Dobaño C, Baird JK. Acquired immunity to malaria. Clin Microbiol Rev. 2009; 22: 13-36. doi: 10.1128/CMR.00025-08 PMID: 19136431

37. Kossodo S, Monso C, Juillard P, Velu T, Goldman M, Grau GE. Interleukin-10 modulates susceptibility in experimental cerebral malaria. Immunology. 1997; 91: 536-540. PMID: 9378491

38. Grau GE, Piguet PF, Vassalli P, Lambert PH. Tumor-necrosis factor and other cytokines in cerebral malaria: experimental and clinical data. Immunol Rev. 1989; 112: 49-70. PMID: 2575074

39. Rudin W, Eugster HP, Bordmann G, Bonato J, Müller M, Yamage M, et al. Resistance to cerebral malaria in tumor necrosis factor-alpha/beta-deficient mice is associated with a reduction of intercellular adhesion molecule-1 up-regulation and T helper type 1 response. Am J Pathol. 1997; 150: 257-66. PMID: 9006341
40. Lifshitz L, Tabak G, Gassmann M, Mittelman M, Neumann D. Macrophages as novel target cells for erythropoietin. Haematologica. 2010; 95: 1823-1831. doi: 10.3324/haematol.2010.025015 PMID: 20595096

41. Lifshitz L, Prutchi-Sagiv S, Avneon M, Gassmann M, Mittelman M, Neumann D. Non-erythroid activities of erythropoietin: Functional effects on murine dendritic cells. Mol Immunol. 2009; 46: 713-721. doi: 10.1016/j.molimm.2008.10.004 PMID: 19027957

42. Nairz M, Schroll A, Moschen AR, Sonnweber T, Theurl M, Theurl I, et al. Erythropoietin contrastingly affects bacterial infection and experimental colitis by inhibiting nuclear factor-κB-inducible immune pathways. Immunity. 2011; 34: 61-74. doi: 10.1016/j.immuni.2011.01.002 PMID: 21256055