Changing plasma cytokine, chemokine and growth factor profiles upon differing malaria transmission intensities
Malaria Journal volume 18, Article number: 406 (2019)
Malaria epidemiological and immunological data suggest that parasite tolerance wanes in the absence of continuous exposure to the parasite, potentially enhancing pathogenesis. The expansion of control interventions and elimination campaigns raises the necessity to better understand the host factors leading to susceptibility or tolerance that are affected by rapid changes in malaria transmission intensity (MTI). Mediators of cellular immune responses are responsible for the symptoms and pathological alterations during disease and are expected to change rapidly upon malaria exposure or cessation.
The plasma concentrations of 30 cytokine, chemokine and growth factors in individuals of all ages from a malaria endemic area of southern Mozambique were compared between 2 years of different MTI: 2010 (lower, n = 234) and 2013 (higher, n = 143). The effect of the year on the correlations between cytokines, chemokines and growth factors and IgGs to Plasmodium falciparum (markers of exposure) was explored. The effects of age, sex, neighbourhood and parasitaemia on analyte levels and their interactions with year were also assessed.
An inverse correlation of several cellular immune mediators with malarial antibodies in 2013, and a lack of correlation or even a positive correlation in 2010 were observed. Most cytokines, chemokines and growth factors, regardless of their immune function, had higher concentrations in 2010 compared with 2013 in P. falciparum-infected and uninfected subjects. Age and neighbourhood showed an effect on analyte concentrations.
The results show a different regulation of the cellular immune response in 2010 vs 2013 which could be related to a loss of immune-tolerance after a decline in MTI in 2010 and previous years, and a rapid re-establishment of tolerance as a consequence of more continuous exposure as MTI began increasing in 2012. Cellular immune mediators warrant further investigation as possible surrogates of MTI-associated host susceptibility or tolerance.
The burden of malaria disease has experienced significant changes in endemic areas in the 21st century. Between 2000 and 2015, an expansion of malaria interventions helped to reduce malaria incidence by 37% globally, and by 42% in Africa. However, between 2015 and 2017 no significant progress was made, with an estimated 219 million cases and 435 000 related deaths in 2017 .
Geographical and temporal changes in malaria transmission intensity (MTI) can affect disease burden and Plasmodium parasite dynamics. Thus, the expansion of control interventions and elimination campaigns raises the necessity to better understand the host factors affected by rapid changes in MTI.
In endemic areas of Africa, naturally acquired immunity (NAI) to malaria is developed with age and exposure to Plasmodium falciparum infection. NAI is suggested to be comprised of two main components: (i) an anti-parasite component, resulting in control of parasite replication and parasite clearance, which takes years to be acquired and is never sterilizing [2, 3]; and (ii) an anti-disease component, consisting of the ability to tolerate parasites asymptomatically, which is acquired rapidly and can result in long periods without malaria symptoms in older individuals [4, 5].
Tolerance is a less understood phenomenon. From the immunological perspective, it is defined as any endogenous mechanism by which a potentially injurious immune response is prevented, suppressed, or shifted to a non-injurious response . In malaria, such tolerance developed by the host is suggested to be multi-factorial, including: (i) the neutralization of parasite toxins and other virulence factors; (ii) immuno-regulatory processes that reduce the damage triggered by excessive immune responses of the host; and (iii) cellular and systemic adaptive responses that limit the deleterious effects associated with stress imposed by pathogens and/or host immunity .
Epidemiological and immunological data suggest that anti-parasite immunity and tolerance wane in the absence of continuous exposure to the parasite , and changes in MTI likely affect anti-malarial immunity. In fact, the geographical distribution of malaria prevalence, morbidity and mortality depends directly upon MTI. In low MTI settings, exposed people are at a higher risk of severe disease. In high MTI settings, severe disease is limited to naïve individuals (visitors, infants), young children, and pregnant women, while the rest of adults tolerate the presence of parasites .
Temporally, immune adults who migrate to non-endemic areas are at significant risk of contracting malaria upon return to an endemic area, even increased susceptibility to severe malaria . Previous studies show that IgG responses to malaria-specific antigens are maintained to a large extent upon cessation of malaria exposure, suggesting a long-lasting anti-parasite immunity ; however, control of pro-inflammatory responses and tolerance to P. falciparum appeared to be reduced . Thus, a rapid decrease in exposure to P. falciparum would result in a lesser development of NAI in children, and a possible partial loss of previous immunity and tolerance in the older population if exposure was very low or discontinued, which may enhance malaria pathogenesis.
Cytokines, chemokines and growth factors mediate cellular responses and are responsible for the symptoms and pathological alterations during disease. The outcome of infection depends on the regulation of pro-inflammatory and anti-inflammatory responses, leading to protection or immunopathology . Thus, blood soluble mediators are expected to change rapidly upon malaria exposure or cessation, and to reflect changes in anti-disease immunity and tolerance associated with varying MTI. Some P. falciparum-specific cytokine responses have been associated with protection against clinical malaria, including interferon gamma (IFN-γ) [13,14,15,16,17], interleukin-10 (IL-10) and tumour necrosis factor (TNF) [18, 19]. Others, such as IL-6, have correlated with increased risk of clinical malaria [12, 20]. Serum cytokines like IL-5 and RANTES appear to be important in the pathogenesis of severe malaria [20,21,22,23].
In this study, changes in cellular immune mediator profiles in 2 years of different MTIs were assessed in individuals of all ages from a malaria endemic area of southern Mozambique. To address this, a comprehensive panel of 30 cytokines, chemokines and growth factors, several of them known to vary upon malaria infection and/or exposure [11, 24, 25], were measured in plasma samples collected in two cross-sectional surveys in 2010 and 2013 in the Manhiça District. Most analytes had higher concentrations in 2010 (lower MTI) compared to 2013 (higher MTI) in P. falciparum infected but also in uninfected subjects. This could be indicative of a loss of immune-tolerance after years of decline in MTI, and a re-establishment of the tolerance as a consequence of a more continuous exposure due to a rise in MTI starting before 2013. The possibility of cytokines as surrogates of loss of parasite tolerance upon changes in MTI deserves further investigation.
Study design, subjects and sample collection
This study was developed in the context of yearly cross-sectional surveys performed to monitor changes in malaria burden and MTI in the Manhiça District, Maputo Province, in southern Mozambique, which started in 2010 and are still on going. The characteristics of the study area have been described in detail elsewhere . The climate is subtropical and transmission of P. falciparum malaria is perennial and of moderate intensity, with two different seasons, a warm and rainy season from November to April, and a cool and dry season the rest of the year . The cross-sectionals were performed at the end of the peak of the transmission season (February–March) and lasted approximately 1 month. Data on monthly mean rain in the 5 months previous to the 2010 and 2013 surveys are shown in Additional file 1: Table S1 (CHIRPS: http://chg.geog.ucsb.edu/data/chirps/). The cumulative rainfall levels between the months of September 2009 and January 2010 was 310.76 mm, and doubled to 625.26 mm between the months of September 2012 and January 2013. Individuals of all ages were selected by random age-stratified sampling from the demographic surveillance system census and were invited to participate in the study. Weighted sampling was done according to the following age groups: Infants < 1 year; 1 ≤ 2 years; 2 ≤ 3 years; 3 ≤ 4 years; 4 ≤ 5 years; 5 ≤ 10 years; 10 ≤ 20 years; 20 ≤ 40 years; 40 ≤ 60 years; ≥ 60 years. Participants were recruited from six different neighbourhoods: Malavele, Manhiça and Maragra, historically considered of low MTI; and Ilha Josina, Palmeira and Taninga, of moderate MTI . A standardized questionnaire was filled-in with basic demographic information, use of malaria control tools and socio-economic status. Axillary temperature was measured and registered in the questionnaire. A blood smear was collected to determine Plasmodium parasitaemia according to standard, quality-controlled procedures .
Haemoglobin was measured using HemoCue portable devices (Ängelholm, Sweden). A blood aliquot was collected in an EDTA microtainer and plasma separated by centrifugation and cryopreserved at − 80 °C for immunological analyses. Blood was also collected onto filter paper for IgG quantification by ELISA , and for parasitaemia quantification by real time quantitative PCR (qPCR), as described elsewhere . Febrile infection was defined as the presence of asexual P. falciparum parasites in blood detected qPCR, together with fever or reported fever during the previous 24 h. Febrile infections with parasites also detected by microscopy were considered clinical malarias. Anemia was defined as haemoglobin < 11 g/dL. Participants presenting with parasitaemia, anaemia, fever or history of fever in the previous 24 h were treated according to standard procedures. All malaria infections were treated with the first-line anti-malarial treatment (Coartem), and anaemia cases received ferrous sulphate, according to national guidelines. Participants presenting signs/symptoms of severity were transferred to the Manhiça District Hospital.
Plasma samples analysed in this study were from 377 participants from the surveys performed in 2010 (N = 234, 96 infected and 138 non-infected) and 2013 (N = 143, 65 infected and 78 non-infected). Years 2010 and 2013 were defined as lower and higher MTI periods, respectively, based on the trends reported in Mozambique, with a decrease in malaria prevalence from 2007 to 2011, and an increase since 2012 [32, 33].
All plasmas from infected individuals (161 participants qPCR positive) and the plasmas from 216 non-infected participants (qPCR negative) were analysed separately to assess cytokine profiles in 2010 and 2013 during an infection and at baseline, respectively. The plasmas from non-infected individuals were randomly selected from the surveys, balancing between years and stratifying by specific age groups (1 ≤ 2, 2 ≤ 5, 5 ≤ 10, 20 ≤ 40 and ≥ 60), with the last two groups enriched to have a sample size powered to address older age effects.
Cytokine, chemokine and growth factor multiplex bead array assay
The Cytokine Human Magnetic 30-Plex Panel from Life Technologies™ was used to measure the concentrations (pg/mL) of the following cytokines, chemokines and growth factors in plasma: epidermal growth factor (EGF), fibroblast growth factor (FGF), granulocyte colony-stimulating factor (G-CSF), granulocyte-macrophage colony-stimulating factor (GM-CSF), hepatocyte growth factor (HGF), vascular endothelial growth factor (VEGF), tumour necrosis factor (TNF), interferon (IFN)-α, IFN-γ, interleukin (IL)-1RA, IL-1β, IL-2, IL-2R, IL-4, IL-5, IL-6, IL-7, IL-8, IL-10, IL-12(p40/p70), IL-13, IL-15, IL-17, IFN-γ induced protein (IP-10), monocyte chemoattractant protein (MCP-1), monokine induced by IFN-γ (MIG), macrophage inflammatory protein (MIP)-1α, MIP-1β and regulated on activation normal T cell expressed and secreted (RANTES) and eotaxin. This panel has been used in previous studies showing several analytes varying upon malaria infection and/or exposure [11, 24, 25].
Twenty-five microlitres of all plasmas were tested in single replicates applying a modification of the manufacturer’s protocol that implies using half the volume of each reagent except for the washing buffer ; this modification was previously tested, showing no difference in assay performance compared to the original protocol. Each plate included 16 serial dilutions (two fold) of a standard sample provided by the vendor with known concentrations of each analyte, two blank controls and three positive controls of high, medium and low concentrations in duplicate prepared from a reference sample for quality assurance/quality control purposes. Samples from infected and non-infected individuals were assayed in two separate batches. In each batch, plates were balanced across the two cross-sectionals and age groups. Samples were acquired on a Luminex® 100/200 instrument and analysed in xPONENT® software 3.1. The concentration of each analyte was obtained by interpolating the median fluorescent intensity (MFI) (after blank MFI subtraction) to a 5-parameter logistic regression curve automatically calculated by xPONENT® software. Any value below the lower limit of detection (mean of blanks + 2 standard deviations) was assigned half the expected concentration at the low limit of quantification for that analyte.
Enzyme-linked immunosorbent assay
Antibody data obtained in a separate study of seroconversion rates performed in the context of the Manhiça cross-sectional surveys were analysed in relation to the cytokine data. Briefly, 2.5 mm disks of dried blood (≈ 1.5 µL of blood) were cut from filter paper spots and incubated overnight at room temperature (RT) with 150 mL of PBS/0.05% Tween 20/0.01% sodium azide (w/v) .
Reconstituted sera were stored at − 20 °C until use. The ELISA assay was performed as previously described . High-binding 96-well titer plates (Immulon 4HBX, Thermo Scientific, Inc.) were coated with merozoite surface protein-1 (MSP-142) (3D7 strain) or apical membrane protein-1 (AMA-1) (FVO strain) produced at the Walter Reed Army Institute of Research (MD, USA), at 0.5 mg/mL of carbonate buffer. PBS with 0.05% Tween (PBS-T) was used to wash plates between incubations. Plates were blocked with 1% skimmed milk powder (Sigma-Aldrich, Inc.) in PBS-T (blocking buffer) for 3 h at RT. Reconstituted antibodies were transferred to the ELISA plates (1/1000 and 1/2000 final dilution for MSP-142- and AMA-1-coated plates, respectively) and incubated overnight at 4 °C. HRP-conjugated rabbit anti-human IgG (Dako/Agilent Technologies, Inc.) was applied at 1/5000 in PBS-T and incubated for 3 h at RT. OPD development substrate was applied and incubated in the dark at RT for 20–25 min, stopping development with the addition of 2 M H2SO4. Plates were read on Bio-Tek ELx 50 plate reader using KC Junior software package (version 1.10, Biotek Instruments Inc.). Normalized optical density (OD) was calculated as the mean background-adjusted OD divided by the mean of a hyperimmune plasma. The same standard dilution was used for normalizing each plate in the study. Test samples were assayed in duplicate and included in the analysis if the coefficient of variation was less than 50% for all values greater than 0.1 OD.
Plasmas from malaria-infected and non-infected volunteers were analysed separately. The studied population was categorized into 6 age groups (≤ 2 years, > 2 ≤ 5 years, > 5 ≤ 10 years, > 10 ≤ 20 years, > 20 ≤ 60, and ≥ 60 years) according to commonly observed immunological patterns. Demographic continuous variables were analysed using the non-parametric Wilcoxon rank-sum test. Comparisons between groups for categorical variables were done using Fisher’s exact test. Concentrations of cytokines, chemokines, growth factors, parasite densities and IgG levels (OD) were log10 transformed for further analysis.
Comparison of cellular analyte concentrations between two or more groups was performed through Wilcoxon rank-sum test or Kruskall–wallis test, respectively. Comparisons of IgG levels between years were performed by Wilcoxon rank sum test with continuity correction for all age groups together or into each age category separately. Trends of IgG levels along age groups were assessed by Spearman correlations (p-trend). The correlations of IgG levels and analyte concentrations were assessed with univariable linear models, with IgG levels as independent variable, separate for each analyte. The effects of year, age group, neighbourhood, sex and parasite density (predictors) on analyte levels (outcome) were assessed through univariable and multivariable (adjusted models) linear regressions, separate for each analyte. Interaction tests were performed to determine if there was an interaction of age group, neighbourhood, sex, parasite density and antibody levels with year on the analyte levels. When interactions were statistically significant, a stratified analysis was performed.
All p-values were considered statistically significant when < 0.05. P-values were adjusted for multiple testing to control the false discovery rate using the Benjamini–Hochberg approach in each study endpoint separately (effect of year in infected subjects, effect of year in uninfected subjects, comparison of infected and uninfected subjects, effect of age, neighbourhood, sex, parasitaemia and AMA-1 and MSP-1 IgG levels) except for the comparison of antibody levels between years into age categories for which the Holm method was used. In all endpoints evaluating cytokine associations of this exploratory study, the Benjamini–Hochberg approach was used because it has more power than family-wise-error rate methods, such as Bonferroni or Holm, and allows identifying as many significant associations as possible while incurring a relatively low proportion of false positives . For the AMA-1 and MSP-1 antibody data a family-wise error rate method was used to be more conservative and because only two markers were analysed. All data collected were analysed using the R software version 3.2.4 (2016-03-10) . The ggplot2 package  was used to perform boxplot graphs and scatter plots. The compare Groups package was used for computing p-trend values . The fmsb package  was used to draw the radar charts with the radarchart function. The reshape and dyplr packages were used for data manipulation purposes  and the ReporteRs package for exporting data tables in Microsoft Word documents .
Malaria prevalence and antibody levels in the 2010 and 2013 surveys
A total of 981 and 980 individuals participated in the malaria cross-sectional surveys in 2010 and 2013, respectively. In 2010, 975 (99.5%) individuals were tested by microscopy and 970 (99%) by qPCR; while 944 (96.3%) and 808 (82.4%) individuals were tested by microscopy and qPCR, respectively, in 2013. Weighted prevalence of P. falciparum infection in 2010 and 2013 were 1.1% and 5.2% by microscopy, and 12.6% and 23.9% by qPCR, respectively; and the numbers of febrile infections were 15 and 43, being 3 and 18 clinical malarias, respectively (Table 1). When stratifying by neighbourhood, Palmeira, Ilha Josina and Taninga showed the highest increases in malaria prevalence (qPCR) from 2010 to 2013, with Palmeira and Ilha Josina presenting the highest increases in clinical cases.
Levels of IgG against AMA-1 and MSP-142 were higher in 2013 compared to 2010 (p < 0.001 for both IgGs) mirroring the MTI trends between 2010 and 2013, and increased with age within each cross-sectional (p trend < 0.001 in all four cases), reflecting the continuous exposure to P. falciparum infection (Fig. 1a). When stratifying by P. falciparum infection during the survey, similar trends were observed in the infected and uninfected populations with overall higher IgG levels in the infected subjects (Fig. 1b).
Description of the subset of participants included in this study
Among infected participants, there were no differences in age median or distribution of age groups between years, while in 2013 the non-infected individuals were younger than in 2010, with more subjects below 10 years old and less above 20 years old (Table 2). Participants’ distribution between neighbourhoods was also different between years in both groups (Table 2). Among the infected subjects, there were no differences in parasitaemia by qPCR between years (Table 2). In total, there were 11 cases of clinical malaria, 8 (12.3%) in 2013 and 3 (3.1%) in 2010 (p = 0.017).
Cytokine, chemokine and growth factor profiles differed between 2010 and 2013
The main aim of the study was to evaluate differences in cellular analyte concentrations between years of different MTI in P. falciparum-infected and uninfected individuals separately, as differences by MTI were expected to be altered by infection. Overall, 25 analytes out of 30 had higher levels in 2010 compared to 2013 in infected and/or uninfected volunteers, except for RANTES that showed an opposite pattern (Fig. 2, Additional file 2 and Additional file 3). The analytes that increased in 2010 were: the pro-inflammatory cytokines IL-1β, the IL-1 inhibitor IL-1RA, TNF and IL-6; the anti-inflammatory cytokine IL-10; the TH2 cytokine IL-13; the regulatory cytokine IL-7; the chemokines IL-8, IP-10, MCP-1, MIP-1α and MIP-1β; and the growth factors EGF, G-CSF, GM-CSF, HGF and VEGF. The regulatory or TH1-related cytokine IL-15, the TH1 cytokine receptor IL-2R, the anti-inflammatory cytokine IL-17 and the growth factor FGF were higher in 2010 only in the infected individuals (Fig. 2 and Additional file 2). The pro-inflammatory cytokine IFN-α, the TH1 cytokine IL-12 and the TH2 IL-4 were higher in 2010 only in the non-infected subjects (Fig. 2 and Additional file 3).
Levels between infected and non-infected individuals were also compared. As expected, P. falciparum infection was associated with higher levels of several cytokines, chemokines and growth factors in both years (Fig. 3): IFN-γ, IL-12, IL-17, IL-10, IL-2, IL-2R, IL-15, MIG, EGF, FGF, G-CSF and HGF. However, some markers were only affected in 2013: IL-1β, IL-6, IL-8, MIP-1α, MIP-1β and VEGF; and IL-5 only in 2010. There were also some analytes that were lower in the infected vs the uninfected individuals: IL-13, IL-7, eotaxin, IP-10 and RANTES in both years; IL-4 in 2010; and IL-1RA in 2013. IFN-α and MCP1 showed opposite patterns depending on the year. Additional file 1: Table S2 shows the effect of the infection on analyte levels and the interaction with year.
Correlations of cytokines, chemokines and growth factors with antibodies differ between 2010 and 2013
The correlations of cellular analyte concentrations with malaria antibodies as markers of exposure were explored. When considering data from 2010 and 2013 together (Table 3), the overall trend showed an inverse correlation of several cytokines and chemokines with antibodies to AMA-1 and MSP-1 in the univariable and multivariable models. In the univariable models, IL-13, IL-1RA, IP-10 and IL-2R negatively correlated with both antibodies, IL-12 was negatively correlated with IgG to AMA-1, and VEGF negatively correlated with IgG to MSP-1. The correlations of AMA-1 and MSP-1 IgG levels with IFN-α, AMA-1 IgG with VEGF and MSP-1 IgG with MCP-1 were different in 2010 compared to 2013 (Table 3, Additional file 4).
Year-dependent effects of age, sex and neighbourhood on analyte concentrations
Effect of age
The effect of age on the cytokine, chemokine and growth factor responses was assessed, with particular interest in the children and the elderly who have a developing and a senescent immune system, respectively. Age had a statistically significant effect on some analyte concentrations in infected and uninfected individuals, with different patterns depending on the analyte and infection status (Figs. 4 and 5). Levels of IL-12, IL-2, IL-2R, IL-15, FGF and HGF decreased with age in both groups. IFN-γ, IL-17, IL-13, IL-4, IL-5, MIG and G-CSF only decreased in the infected subjects, with some of the analytes slightly increasing in the elderly. Although not statistically significant, other analytes also showed a decreasing trend with respect to age in the infected group like IFN-α, IL-10, IL-15 and GM-CSF. The most marked drop in concentrations was always between age 2 and 10 years old. RANTES levels increased with age in both groups. Eotaxin and IP-10 profiles showed a U-shape with higher values in the groups < 2 years and > 60 years old. A trend of U-shape was also observed for some other analytes like the TH2 cytokines IL-4, IL-5, IL-13 and the G-CSF in the infected group. Overall, the children and the elderly showed opposite patterns for some cytokines but similar for others, reflecting the intrinsic characteristics of the immune system in these two age groups and suggesting a process of immunosenescence in the older one.
The interaction of age with year on analyte concentrations was assessed, and statistically significant interactions were found for eotaxin and VEGF in infected individuals, and IL-6, MCP-1, EGF, G-CSF, HGF and VEGF in uninfected individuals (Additional file 1: Table S3 and Additional file 5), however, after correcting by multiple testing significance disappeared (Additional file 1: Table S3).
Effect of sex
Sex did not have an effect on analyte concentrations in either infected or uninfected subjects (Additional file 1: Table S4). However, an interaction of sex with year was found for RANTES and GM-CSF concentrations in infected individuals (p-values = 0.002 and 0.041, respectively), although the significance was only maintained for RANTES after adjusting for multiple comparisons. RANTES presented higher levels in males compared to females in 2010, but the opposite was observed in 2013 (Additional file 6). Sex and year interactions were found in uninfected individuals for TNF, IFN-γ and IL-4, although they were not significant after adjusting for multiple testing (Additional file 1: Table S4 and Additional file 6).
Effect of neighbourhood
The neighbourhood had an effect on analyte concentrations only in the uninfected subjects (Fig. 6), with lower levels of several of them (IL-1β, IL-1RA, IL-6, IFN-α, IL-8, MCP-1, MIP-1α, MIP-1β and VEGF) in Palmeira and Taninga, and a trend of higher levels in Ilha Josina and Maragra; and the opposite trend for RANTES. No effect was observed in infected individuals (Additional file 7). However, the neighbourhood only interacted with year for IL-5 and FGF in the infected individuals, and the interaction only remained significant for FGF after adjusting by multiple testing (Additional file 1: Table S5).
Effect of parasitaemia on analyte concentrations
Among P. falciparum-infected subjects, parasite density had an effect on different analytes (Table 4). IL-10, IL-2R, IL-12 and MCP-1 increased with higher parasitaemia, and IFN-γ, IL-13, IL-5, IL-7 and eotaxin decreased with increasing parasitaemia (Table 4). The effect only remained significant for IL-10, IL-13 and IL-2R when adjusting for multiple testing. The effect of parasitaemia on IL-10 and GM-CSF was different by year (Table 4 and Additional file 8). However, after adjusting for multiple testing, interaction only remained statistically significant for IL-10, showing a stronger correlation with parasitaemia in 2013 compared to 2010.
Adjusted effect of year on cytokine, chemokine and growth factor concentrations
The effect of the year on analyte concentrations was also assessed in regression models adjusting by age and neighbourhood. Results were similar to the unadjusted analyses, with lower levels of almost the same analytes in 2013 compared to 2010 in infected and uninfected volunteers, and RANTES showing an opposite pattern (Table 5). Analytes that were lower in 2013 in infected and non-infected individuals were: the pro-inflammatory cytokines IL-1β, IL-1RA, TNF, IL-6; the TH1 cytokine IL-2R; the anti-inflammatory cytokine IL-10; the TH2 cytokine IL-13; the regulatory cytokine IL-7; the chemokines IL-8, IP-10, MCP1, MIP-1α and MIP-1β; and the growth factors EGF, G-CSF, GM-CSF, HGF and VEGF. The regulatory or TH1-related cytokine IL-15 and the growth factor FGF were lower in 2013 only in the infected subjects. In the non-infected individuals, there were some additional analytes also with lower levels in 2013: the pro-inflammatory cytokine IFN-α, the TH1 cytokine IL-12, and the TH2 cytokine IL-4.
In spite of significant reductions in the burden of malaria over the last 10 years , scarce monitoring has been done on the impact on malaria immunity. Mozambique is one of the countries with the highest malaria burden in the world, although there are increasing efforts to move towards elimination in southern provinces . Malaria burden decreased in all provinces from a malaria incidence rate over 6 million in 2007 to 3 million in 2011, but an increase was observed again in 2012 reaching near 6 million in 2014 [32, 33, 45]. In this study, the association of shifting epidemiological patterns with the systemic cellular immune profiles, with or without a current infection, was assessed through quantification of several plasma cytokines, chemokines and growth factors.
Previous studies have found that the balance between pro-inflammatory and anti-inflammatory cytokines determines host protection and injury [22, 46, 47], and the same for growth factors and chemokines [48, 49]. Given the importance of these proteins in the immune response and the consequent control of immunopathology, it was hypothesized that they could be importantly affected by changes in MTI. Results from this study show that in both infected and uninfected subjects most of the analytes were at higher concentrations in 2010 than 2013, suggesting the possibility of a blunted cytokine response with higher MTI, which could be associated to host tolerance. Albeit weakly, analytes were negatively associated with antibody markers of malaria exposure, further suggesting a blunted response with higher exposure. For some analytes, the correlations with antibodies differed by year, indicating a different regulation of the cellular response. The higher cytokine, chemokine and growth factor levels in 2010 could be related to a loss of immune tolerance after a decline in MTI, and their lower levels in 2013 could be indicative of a rapid re-establishment of tolerance as a consequence of more continuous exposure. Thus, these cellular immune mediators could be candidate surrogates of MTI and the associated host susceptibility and tolerance deserving further investigation.
Recent studies suggest that repeated exposure to P. falciparum leads to the establishment of tolerance , which may be associated with the loss and/or altered function of several immune cell types, including γδ T cells , αβ T cells [50,51,52], B cells , and myeloid cells  that show less proliferation and cytokine production . In addition, malaria exposure also induces changes in the innate immune response [55, 56]. A recent study identifies individuals who are primed to respond favorably to P. falciparum infection by controlling inflammatory symptoms (disease tolerance) and parasitaemia . Accordingly, the lower levels of cytokines, chemokines and growth factors observed with increasing MTI could be reflecting a more controlled cellular response attributed to some degree of immuno-tolerance. In agreement with these results, a recent study in children with clinical malaria from Ghana reported a decrease of several cytokines with increasing MTI .
Another study showed diminished cytokine response with age in children, in part probably due to increased exposure, which could also be indicative of tolerogenic mechanisms . A decrease in MTI is usually followed by a delayed acquisition of immunity to clinical and severe disease; data from a health facility in one of the higher endemic areas of Manhiça confirm that after years of sustained decrease in malaria incidence, there is an increase in the mean age of clinical malaria events and severe forms of disease . Accordingly, a study on infected pregnant women from Manhiça between 2003 and 2012, showed a decrease of antibodies against P. falciparum and an increase of malaria adverse consequences after the decline of MTI . Interestingly, in the present study, the differing cytokine, chemokine and growth factor levels detected in the non-infected individuals between 2010 and 2013 suggests that the impact of MTI could be beyond the response to P. falciparum infection at the time of the survey, suggesting that the rate of previous exposures could shape the basal immunological status. This difference in the basal immune system agrees with the observation by Tran et al.  of a marked different transcriptomic and cellular profiles between malaria-protected and malaria-susceptible children prior to Plasmodium infection. This study also suggests that previous exposure to malaria, and possibly other pathogens or commensals, shape the basal immune system.
In the same line, a previous study by our group in semi-immune African adults, migrants and European travellers, showed that naïve adults had stronger cytokine responses upon infection than semi-immune adults; and that migrants, in the absence of continuous exposure, presented higher concentrations of cytokines and chemokines than semi-immune individuals, but lower than individuals with a first infection . There, increased levels of IL-2, IFN-γ, IL-8 and IL-5 were associated with loss of exposure to P. falciparum . Similarly, in the present study, higher levels of IL-8 were observed in 2010. In addition, increased levels of IL-6 and IL-10, and decreased levels of RANTES were detected in 2010; these three cytokines have been associated with severe malaria [22, 47, 62]. Several studies indicate a strong positive role of RANTES against Plasmodium species susceptibility [63, 64] and malaria severity , being down-regulated in severe malaria compared to uncomplicated malaria . This immunomodulatory role of RANTES reducing the pathogenesis of malaria would agree with this chemokine being at lower levels in the lower MTI period in which a loss of immune-tolerance is hypothesized.
Despite the lower cellular immune mediators observed in 2013, infected individuals showed higher levels of several cytokines and chemokines compared to the uninfected ones, differences that were not observed in 2010; among them, the pro-inflammatory cytokines IFN-α, IL-1β and IL-6, the chemokines IL-8, MCP-1, MIP-1α and MIP-1β and the growth factor VEGF. This may be related to the lower basal levels of cellular immune mediators in 2013 compared to 2010. In 2010, basal levels of cellular mediators were so high that infections might have not been able to induce further increases; whereas in 2013, cytokine, chemokine and growth factor basal levels were low and infections increased them significantly, though still below the basal levels from 2010. This observation supports the hypothesis of a more controlled, non-harmful cellular immune response upon re-infection in 2013.
Despite the overall decrease in cellular immune mediators in 2013, a higher prevalence of clinical malaria (12%) was registered compared to 2010 (3%) among the infected individuals, similar to the prevalence in the surveys (18/162 vs 3/106). However, the total number of clinical malaria cases was low, and clinical malaria depends not only on tolerogenic responses but also on anti-parasite immunity. During lower MTI periods there is a delay in the acquisition of immunity in younger children and a reduction in premunition in older children and adults; thus, in the context of a sudden increase of MTI, anti-parasite immunity may be acquired slower than tolerogenic responses. It is likely that a higher number of clinical cases would have been observed in 2013 had some level of tolerance not been developed.
As expected, among the infected individuals the level of parasitaemia affected several analytes, some of them similar to previous reports [11, 66]. Interestingly, a different effect of parasitaemia on IL-10 and GM-CSF was found depending on the year, reinforcing the idea of different cellular immune responses at different MTIs. A different effect of parasitaemia on cytokine levels depending on MTI was also recently reported , with several cytokines correlating with parasitaemia only in the lower MTI area.
Age-related differences in immunity are suggested to also explain the different susceptibility to malaria disease in children vs adults  independently of cumulative exposure [8, 68]. Previous studies showed that adult newcomers into a hyperendemic area developed NAI relatively quickly while their children remained susceptible, suggesting that age is relevant for NAI development. NAI against malaria may also be diminished at older ages because of changes in the immune system, known as immunosenescence, which contribute to make elderly more susceptible to infections, cancer and autoimmunity. For example, there is a decline in the output of regulatory T cells after the age of 50 [69, 70] and altered cytokine levels have been observed [71, 72]. How immunosenescence may affect the acquired protection against malaria is still unknown, and the impact of changing MTI could be different. Interestingly, results from this study show that different cellular immune mediators present diverse age patterns, and some were different in different MTI periods while others not. For example, the TH1 cytokines IL-2, IL-12 and IL-2R and the TH1 related IL-15 decreased continuously with age; while eotaxin, IP-10 and the TH2 cytokines IL-13, IL-4 and IL-5 presented a U-shape, i.e. higher levels in children and the elderly and lower levels in young adults. Previous studies have shown that some cytokines increase in the elderly , while others seem to decrease [70, 73], but none of the cytokines were the same observed in this study. A U-shape suggests similarities between the children and elderly immune systems. High levels of some analytes in the age 1-2 years old group could be related to a lower immunity due to previous lack of exposure to malaria, but also to the intrinsic different characteristics of the immune system in early life. High levels of some analytes in the elderly could be related to the process of immunosenescence that may cause loss of anti-disease immunity, therefore behaving similar to a first infection or like a less exposed population presenting higher cellular responses upon re-challenge .
Previous studies have also shown that the severity of malaria infection differs between males and females , with men developing more severe parasitaemia and pathology than women . In addition, it has been described that in general women produce more intense humoral and cell-mediated immune responses than males [77, 78]. However, no effect of sex was observed on the analyte concentrations in either infected or uninfected subjects, although a sex interaction was found with year for TNF, IFN-γ, IL-4, RANTES and GM-CSF, the impact of which is difficult to interpret. MTI is also expected to present spatial heterogeneities due to social and natural factors  but here the neighbourhood did not have an effect on the association between year and cytokine, chemokine or growth factor levels.
There are some limitations in this study that could have affected the results. First of all, the different storage times for plasma samples collected in 2010 vs samples collected in 2013, which may have differently affected the cytokine concentrations. However, storage time tends to decrease the concentration of most of the cytokines . Therefore, if there was any general storage effect, it would have been in the opposite direction of what was observed. Secondly, total cytokines in plasma samples were measured instead of cytokines produced by peripheral blood mononuclear cells (PBMCs) upon in vitro stimulations with malaria antigens, which limits the interpretation of malaria-specific responses. However, PBMCs were not available. Third, many comparisons were performed and despite adjusting for multiple testing, some of the significant differences may be false positives. Nevertheless, many consistent and biologically feasible associations of cellular mediators with MTI were found. Fourth, the HIV serostatus of the study participants was not known, and the study area has a high HIV prevalence  which may be affecting the results. Finally, it is important to keep in mind that the average rainfall before 2013 survey doubled that before 2010 survey, therefore exposure to infections other than malaria or hydration of study subjects may also vary between years, affecting the systemic homeostatic profile of the inhabitants and their basal immune activation and analyte concentrations, thus it may have influenced or confounded the results. However, the higher P. falciparum IgG levels in 2013 compared to 2010 across all age groups and negative correlations with cellular immune mediators suggests that P. falciparum exposure could be driving the changes described. Future studies should have a longitudinal design and comprise different MTI periods to address the role of cellular immune mediators as surrogates of MTI and tolerance/susceptibility to malaria. In addition, controlling for other factors influencing the basal immune system would also be important. Understanding the mechanisms of longevity of immunity and tolerance may help in developing approaches to promote long-lasting anti-disease immunity, even after sustained interrupted exposure, which could be used in areas of low MTI and in elimination campaigns.
Cytokine, chemokine and growth factor profiles varied between years of different MTIs in P. falciparum infected and non-infected individuals, which could be related to a loss of immune tolerance when MTI declines and a rapid re-establishment of immune tolerance after an increase in P. falciparum exposure. This finding warrants further investigations on plasma cellular immune mediators as potential surrogate candidates of MTI and the associated host susceptibility and tolerance.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
malaria transmission intensity
naturally acquired immunity
tumour necrosis factor
regulated upon activation, normal T cell expressed, and secreted
real-time quantitative PCR
enzyme-linked immunosorbent assay
merozoite surface protein
apical membrane antigen
phosphate buffered saline
Tween-phosphate buffered saline
epidermal growth factor
fibroblast growth factor
granulocyte colony-stimulating factor
granulocyte-macrophage colony-stimulating factor
hepatocyte growth factor
vascular endothelial growth factor
IFN-γ induced protein
monocyte chemoattractant protein
monokine induced by IFN-γ
macrophage inflammatory protein
- TH :
T helper lymphocyte
peripheral blood mononuclear cells
WHO. World malaria report 2018. Geneva: World Health Organization; 2018. https://www.who.int/malaria/publications/world-malaria-report-2018/en/. Accessed 1 Apr 2019
Marsh K, Kinyanjui S. Immune effector mechanisms in malaria. Parasite Immunol. 2006;28:51–60.
Langhorne J, Ndungu FM, Sponaas AM, Marsh K. Immunity to malaria: more questions than answers. Nat Immunol. 2008;9:725–32.
Portugal S, Moebius J, Skinner J, Doumbo S, Doumtabe D, Kone Y, et al. Exposure-dependent control of malaria-induced inflammation in children. PLoS Pathog. 2014;10:e1004079.
Jagannathan P, Kim CC, Greenhouse B, Nankya F, Bowen K, Eccles-James I, et al. Loss and dysfunction of Vdelta2(+) gammadelta T cells are associated with clinical tolerance to malaria. Sci Transl Med. 2014;6:251ra117.
Riley EM, Wahl S, Perkins DJ, Schofield L. Regulating immunity to malaria. Parasite Immunol. 2006;28:35–49.
Galatas B, Bassat Q, Mayor A. Malaria parasites in the asymptomatic: looking for the hay in the haystack. Trends Parasitol. 2016;32:296–308.
Doolan DL, Dobano C, Baird JK. Acquired immunity to malaria. Clin Microbiol Rev. 2009;22:13–36.
Farnert A, Wyss K, Dashti S, Naucler P. Duration of residency in a non-endemic area and risk of severe malaria in African immigrants. Clin Microbiol Infect. 2015;21:494–501.
Moncunill G, Mayor A, Jimenez A, Nhabomba A, Casas-Vila N, Puyol L, et al. High antibody responses against Plasmodium falciparum in immigrants after extended periods of interrupted exposure to malaria. PLoS ONE. 2013;8:e73624.
Moncunill G, Mayor A, Bardaji A, Puyol L, Nhabomba A, Barrios D, et al. Cytokine profiling in immigrants with clinical malaria after extended periods of interrupted exposure to Plasmodium falciparum. PLoS ONE. 2013;8:e73360.
Rovira-Vallbona E, Moncunill G, Bassat Q, Aguilar R, Machevo S, Puyol L, et al. Low antibodies against Plasmodium falciparum and imbalanced pro-inflammatory cytokines are associated with severe malaria in Mozambican children: a case-control study. Malar J. 2012;11:181.
Reece WH, Pinder M, Gothard PK, Milligan P, Bojang K, Doherty T, et al. A CD4(+) T-cell immune response to a conserved epitope in the circumsporozoite protein correlates with protection from natural Plasmodium falciparum infection and disease. Nat Med. 2004;10:406–10.
Luty AJ, Lell B, Schmidt-Ott R, Lehman LG, Luckner D, Greve B, et al. Interferon-gamma responses are associated with resistance to reinfection with Plasmodium falciparum in young African children. J Infect Dis. 1999;179:980–8.
John CC, Moormann AM, Sumba PO, Ofulla AV, Pregibon DC, Kazura JW. Gamma interferon responses to Plasmodium falciparum liver-stage antigen 1 and thrombospondin-related adhesive protein and their relationship to age, transmission intensity, and protection against malaria. Infect Immun. 2004;72:5135–42.
Todryk SM, Bejon P, Mwangi T, Plebanski M, Urban B, Marsh K, et al. Correlation of memory T cell responses against TRAP with protection from clinical malaria, and CD4 CD25 high T cells with susceptibility in Kenyans. PLoS ONE. 2008;3:e2027.
Lyke KE, Daou M, Diarra I, Kone A, Kouriba B, Thera MA, et al. Cell-mediated immunity elicited by the blood stage malaria vaccine apical membrane antigen 1 in Malian adults: results of a Phase I randomized trial. Vaccine. 2009;27:2171–6.
Kurtis JD, Lanar DE, Opollo M, Duffy PE. Interleukin-10 responses to liver-stage antigen 1 predict human resistance to Plasmodium falciparum. Infect Immun. 1999;67:3424–9.
Migot-Nabias F, Deloron P, Ringwald P, Dubois B, Mayombo J, Minh TN, et al. Immune response to Plasmodium falciparum liver stage antigen-1: geographical variations within Central Africa and their relationship with protection from clinical malaria. Trans R Soc Trop Med Hyg. 2000;94:557–62.
Robinson LJ, D’Ombrain MC, Stanisic DI, Taraika J, Bernard N, Richards JS, et al. Cellular tumor necrosis factor, gamma interferon, and interleukin-6 responses as correlates of immunity and risk of clinical Plasmodium falciparum malaria in children from Papua New Guinea. Infect Immun. 2009;77:3033–43.
John CC, Opika-Opoka R, Byarugaba J, Idro R, Boivin MJ. Low levels of RANTES are associated with mortality in children with cerebral malaria. J Infect Dis. 2006;194:837–45.
Prakash D, Fesel C, Jain R, Cazenave PA, 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.
Lourembam SD, Sawian CE, Baruah S. Dysregulation of cytokines expression in complicated falciparum malaria with increased TGF-beta and IFN-gamma and decreased IL-2 and IL-12. Cytokine. 2013;64:503–8.
Moncunill G, Mayor A, Jimenez A, Nhabomba A, Puyol L, Manaca MN, et al. Cytokine and antibody responses to Plasmodium falciparum in naive individuals during a first malaria episode: effect of age and malaria exposure. PLoS ONE. 2013;8:e55756.
Quelhas D, Puyol L, Quinto L, Nhampossa T, Serra-Casas E, Macete E, et al. Intermittent preventive treatment with sulfadoxine-pyrimethamine does not modify plasma cytokines and chemokines or intracellular cytokine responses to Plasmodium falciparum in Mozambican children. BMC Immunol. 2012;13:5.
Aranda C, Aponte JJ, Saute F, Casimiro S, Pinto J, Sousa C, et al. Entomological characteristics of malaria transmission in Manhica, a rural area in southern Mozambique. J Med Entomol. 2005;42:180–6.
Alonso PL, Sacarlal J, Aponte JJ, Leach A, Macete E, Milman J, et al. Efficacy of the RTS, S/AS02A vaccine against Plasmodium falciparum infection and disease in young African children: randomised controlled trial. Lancet. 2004;364:1411–20.
WHO. Malaria surveillance, monitoring and evaluation: a reference manual. Geneva: World Health Organization; 2018.
Swysen C, Vekemans J, Bruls M, Oyakhirome S, Drakeley C, Kremsner P, et al. Development of standardized laboratory methods and quality processes for a phase III study of the RTS, S/AS01 candidate malaria vaccine. Malar J. 2011;10:223.
Stewart L, Gosling R, Griffin J, Gesase S, Campo J, Hashim R, et al. Rapid assessment of malaria transmission using age-specific sero-conversion rates. PLoS ONE. 2009;4:e6083.
Mayor A, Serra-Casas E, Bardaji A, Sanz S, Puyol L, Cistero P, et al. Sub-microscopic infections and long-term recrudescence of Plasmodium falciparum in Mozambican pregnant women. Malar J. 2009;8:9.
United Nations Development Programme. Human Development Index; 2014. http://hdr.undp.org/en/content/human-development-index-hdi. Accessed 28 Mar 2016.
Mozambique Ministry of Health, National Malaria Control Programme. Annual Report; 2014.
Ochola LA, Ayieko C, Kisia L, Magak NG, Shabani E, Ouma C, et al. Changes in antigen-specific cytokine and chemokine responses to Plasmodium falciparum antigens in a highland area of Kenya after a prolonged absence of malaria exposure. Infect Immun. 2014;82:3775–82.
Corran PH, Cook J, Lynch C, Leendertse H, Manjurano A, Griffin J, et al. Dried blood spots as a source of anti-malarial antibodies for epidemiological studies. Malar J. 2008;7:195.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc. 1995;57:289–300.
R: A language and environment for statistical computing. 2015. Available from http://www.r-project.org/. Accessed 1 Apr 2019
Wickham H. ggplot2: elegant graphics for data analysis. New York: Springer; 2009.
Subirana I, Sanz H, Vila J. Building bivariate tables: the compare groups. Package for R. J Stat Softw. 2014. https://doi.org/10.18637/jss.v057.i12.
Nakazawa, M. R package "fmsb": Functions for medical statistics book with some demographic data. 2018. p.1-51
Wickham H. Reshaping data with the reshape package. J Stat Softw. 2007;21:1–20.
Gohel D. ReporteRs: microsoft word and PowerPoint documents generation. Version 0.8.8; 2017. https://davidgohel.github.io/ReporteRs/index.html. Accessed 1 Apr 2019
Desai M, ter Kuile FO, Nosten F, McGready R, Asamoa K, Brabin B, et al. Epidemiology and burden of malaria in pregnancy. Lancet Infect Dis. 2007;7:93–104.
Aide P, Candrinho B, Galatas B, Munguambe K, Guinovart C, Luis F, et al. Setting the scene and generating evidence for malaria elimination in Southern Mozambique. Malar J. 2019;18:190.
Moonasar D, Maharaj R, Kunene S, Candrinho B, Saute F, Ntshalintshali N, et al. Towards malaria elimination in the MOSASWA (Mozambique, South Africa and Swaziland) region. Malar J. 2016;15:419.
Day NP, Hien TT, Schollaardt T, Loc PP, Chuong LV, Chau TT, et al. The prognostic and pathophysiologic role of pro- and antiinflammatory cytokines in severe malaria. J Infect Dis. 1999;180:1288–97.
Lyke KE, Burges R, Cissoko Y, Sangare L, Dao M, Diarra I, et al. Serum levels of the proinflammatory cytokines interleukin-1 beta (IL-1beta), IL-6, IL-8, IL-10, tumor necrosis factor alpha, and IL-12(p70) in Malian children with severe Plasmodium falciparum malaria and matched uncomplicated malaria or healthy controls. Infect Immun. 2004;72:5630–7.
Wilson NO, Bythwood T, Solomon W, Jolly P, Yatich N, Jiang Y, et al. Elevated levels of IL-10 and G-CSF associated with asymptomatic malaria in pregnant women. Infect Dis Obstet Gynecol. 2010;2010:e317430.
Ioannidis LJ, Nie CQ, Hansen DS. The role of chemokines in severe malaria: more than meets the eye. Parasitology. 2014;141:602–13.
Illingworth J, Butler NS, Roetynck S, Mwacharo J, Pierce SK, Bejon P, et al. Chronic exposure to Plasmodium falciparum is associated with phenotypic evidence of B and T cell exhaustion. J Immunol. 2013;190:1038–47.
Hafalla JC, Claser C, Couper KN, Grau GE, Renia L, de Souza JB, et al. The CTLA-4 and PD-1/PD-L1 inhibitory pathways independently regulate host resistance to Plasmodium-induced acute immune pathology. PLoS Pathog. 2012;8:e1002504.
Horne-Debets JM, Faleiro R, Karunarathne DS, Liu XQ, Lineburg KE, Poh CM, et al. PD-1 dependent exhaustion of CD8+ T cells drives chronic malaria. Cell Rep. 2013;5:1204–13.
Weiss GE, Crompton PD, Li S, Walsh LA, Moir S, Traore B, et al. Atypical memory B cells are greatly expanded in individuals living in a malaria-endemic area. J Immunol. 2009;183:2176–82.
Pinzon-Charry A, Woodberry T, Kienzle V, McPhun V, Minigo G, Lampah DA, et al. Apoptosis and dysfunction of blood dendritic cells in patients with falciparum and vivax malaria. J Exp Med. 2013;210:1635–46.
Natama HM, Moncunill G, Rovira-Vallbona E, Sanz H, Sorgho H, Aguilar R, et al. Modulation of innate immune responses at birth by prenatal malaria exposure and association with malaria risk during the first year of life. BMC Med. 2018;16:198.
Schrum JE, Crabtree JN, Dobbs KR, Kiritsy MC, Reed GW, Gazzinelli RT, et al. Cutting edge: Plasmodium falciparum induces trained innate immunity. J Immunol. 2018;200:1243–8.
Tran TM, Guha R, Portugal S, Skinner J, Ongoiba A, Bhardwaj J, et al. A molecular signature in blood reveals a role for p53 in regulating malaria-induced inflammation. Immunity. 2019;51(750–765):e10.
Ademolue TW, Aniweh Y, Kusi KA, Awandare GA. Patterns of inflammatory responses and parasite tolerance vary with malaria transmission intensity. Malar J. 2017;16:145.
Farrington L, Vance H, Rek J, Prahl M, Jagannathan P, Katureebe A, et al. Both inflammatory and regulatory cytokine responses to malaria are blunted with increasing age in highly exposed children. Malar J. 2017;16:499.
Galatas B, Guinovart C, Bassat Q, Aponte JJ, Nhamussua L, Macete E, Saute F, et al. A prospective cohort study to assess the micro-epidemiology of Plasmodium falciparum clinical malaria in Ilha Josina Machel (Manhica, Mozambique). Malar J. 2016;15:444.
Ndam NT, Mbuba E, Gonzalez R, Cistero P, Kariuki S, Sevene E, et al. Resisting and tolerating P. falciparum in pregnancy under different malaria transmission intensities. BMC Med. 2017;15:130.
Were T, Hittner JB, Ouma C, Otieno RO, Orago AS, Ong’echa JM, et al. Suppression of RANTES in children with Plasmodium falciparum malaria. Haematologica. 2006;91:1396–9.
Bostrom S, Giusti P, Arama C, Persson JO, Dara V, Traore B, et al. Changes in the levels of cytokines, chemokines and malaria-specific antibodies in response to Plasmodium falciparum infection in children living in sympatry in Mali. Malar J. 2012;11:109.
Reutersward P, Bergstrom S, Orikiiriza J, Lindquist E, Andersson Svahn H, Ayoglu B, et al. Levels of human proteins in plasma associated with acute paediatric malaria. Malar J. 2018;17:426.
Bujarbaruah D, Kalita MP, Baruah V, Basumatary TK, Hazarika S, Begum RH, et al. RANTES levels as a determinant of falciparum malaria severity or recovery. Parasite Immunol. 2017;39:e12452.
Rodrigues-da-Silva RN, Lima-Junior Jda C, Fonseca Bde P, Antas PR, Baldez A, Storer FL, Santos F, et al. Alterations in cytokines and haematological parameters during the acute and convalescent phases of Plasmodium falciparum and Plasmodium vivax infections. Mem Inst Oswaldo Cruz. 2014;109:154–62.
Baird JK, Masbar S, Basri H, Tirtokusumo S, Subianto B, Hoffman SL. Age-dependent susceptibility to severe disease with primary exposure to Plasmodium falciparum. J Infect Dis. 1998;178:592–5.
Baird JK. Host age as a determinant of naturally acquired immunity to Plasmodium falciparum. Parasitol Today. 1995;11:105–11.
Lang PO, Govind S, Aspinall R. Reversing T cell immunosenescence: why, who, and how. Age (Dordr). 2013;35:609–20.
Chou JP, Effros RB. T cell replicative senescence in human aging. Curr Pharm Des. 2013;19:1680–98.
Ventura MT, Casciaro M, Gangemi S, Buquicchio R. Immunosenescence in aging: between immune cells depletion and cytokines up-regulation. Clin Mol Allergy. 2017;15:21.
Ter Horst R, Jaeger M, Smeekens SP, Oosting M, Swertz MA, Li Y, et al. Host and environmental factors influencing individual human cytokine responses. Cell. 2016;167:1111–24.
Huang H, Patel DD, Manton KG. The immune system in aging: roles of cytokines, T cells and NK cells. Front Biosci. 2005;10:192–215.
Kollmann TR, Levy O, Montgomery RR, Goriely S. Innate immune function by Toll-like receptors: distinct responses in newborns and the elderly. Immunity. 2012;37:771–83.
Klein SL. Hormonal and immunological mechanisms mediating sex differences in parasite infection. Parasite Immunol. 2004;26:247–64.
Landgraf B, Kollaritsch H, Wiedermann G, Wernsdorfer WH. Parasite density of Plasmodium falciparum malaria in Ghanaian schoolchildren: evidence for influence of sex hormones? Trans R Soc Trop Med Hyg. 1994;88:73–4.
Roberts CW, Walker W, Alexander J. Sex-associated hormones and immunity to protozoan parasites. Clin Microbiol Rev. 2001;14:476–88.
Hernandez-Bello R, Nava-Castro K, Muniz-Hernandez S, Nava-Luna P, Trejo-Sanchez I, Tiempos-Guzman N, et al. Beyond the reproductive effect of sex steroids: their role during immunity to helminth parasite infections. Mini Rev Med Chem. 2012;12:1071–80.
Zacarias OP, Andersson M. Spatial and temporal patterns of malaria incidence in Mozambique. Malar J. 2011;10:189.
de Jager W, Bourcier K, Rijkers GT, Prakken BJ, Seyfert-Margolis V. Prerequisites for cytokine measurements in clinical trials with multiplex immunoassays. BMC Immunol. 2009;10:52.
Gonzalez R, Munguambe K, Aponte J, Bavo C, Nhalungo D, Macete E, et al. High HIV prevalence in a southern semi-rural area of Mozambique: a community-based survey. HIV Med. 2012;13:581–8.
We thank all the participants of the study. We thank the CISM field workers, phlebotomists, microscopists and other lab workers. Augusto Nhabomba, Chenjerai Jairoce, Berta Grau and Alfons Jiménez for assistance with management of samples and serological analysis, and Evelina Angov and Sheetij Dutta for the P. falciparum proteins. Laura Puyol, Diana Barrios and Helder Bulo for laboratory logistics support.
The study received funds from MINECO (SAF2008-00743, RYC-2008-02631 to C.D.), the Spanish Agency for International Cooperation and Development (AECID), the Fundación Ramon Areces, the Agència de Gestió d’Ajuts Universitaris i de Recerca (AGAUR 2010FI_B 00168 to J.J.C.), the Department of Health of the Catalan Government (SLT006/17/00109 to G.M.) and PATH Malaria Vaccine Initiative. ISGlobal is a member of the CERCA Programme, Generalitat de Catalunya.
Ethics approval and consent to participate
Written informed consent was obtained from parents or guardians of any child participant before sample collection. Approval for the protocols was obtained from the Hospital Clínic of Barcelona Ethics Review Committee and the National Mozambican Ethics Review Committee. All methods were performed in accordance with the relevant guidelines and regulations.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Monthly mean rainfall in the study area in the 5 months previous to the 2010 and 2013 surveys. Table S2. Effect of infection on analyte concentrations and the interactions with year. Table S3. Interaction of age with year on analyte levels in infected and uninfected individuals. Table S4. Effect of sex on analyte concentrations and the interactions with year in infected and uninfected individuals. Table S5. Interaction of neighborhood with year on analyte levels in infected and uninfected individuals.
Differences in cellular immune mediator concentrations between 2010 (lower MTI) and 2013 (higher MTI) in P. falciparum-infected subjects. Box plots representing the median and interquartile range of each analyte concentration (log10 pg/mL) in infected subjects recruited in 2010 and in 2013. Levels between both years were compared by Wilcoxon rank-sum test and p-values were adjusted for multiple testing by Benjamini–Hochberg approach.
Differences in cellular immune mediator concentrations between 2010 (lower MTI) and 2013 (higher MTI) in uninfected subjects. Box plots representing the median and interquartile range of each analyte concentration (log10 pg/mL) in uninfected subjects recruited in 2010 and in 2013. Levels between both years were compared by Wilcoxon rank-sum test and p-values were adjusted for multiple testing by the Benjamini–Hochberg approach.
Correlations of VEGF, MCP-1 and IFN-α with P. falciparum antibody levels stratified by year. Data are presented in scatter plots with trend lines, with analytes concentrations in the Y axes and antibody levels in the X axes. R2 and p-values were obtained through linear regression models. Only analytes that had a statistically significant interaction with year for the correlation with antibodies are shown.
Cellular immune mediator concentrations in 2010 (low MTI) and 2013 (high MTI) stratified by age group. Box plots representing the median and interquartile range of analytes concentrations (log10 pg/mL) in P. falciparum infected (a) and uninfected (b) subjects. Only analytes in which age and year had a significant p-value for the interaction test (before correcting for multiple testing) are shown.
Cellular immune mediator concentrations in 2010 (low MTI) and 2013 (high MTI) stratified by sex. Box plots representing the median and interquartile range of analyte concentrations (log10 pg/mL) in P. falciparum infected (a) and uninfected (b) subjects. Only analytes in which sex and year had a significant p-value for the interaction test (before correcting for multiple testing) are shown.
Differences in cellular immune mediator concentrations between areas in P. falciparum infected subjects. Box plots representing the median and interquartile range of each analyte concentration (log10 pg/mL) in infected subjects stratified by neighborhood. Levels between areas have been compared by Kruskal–Wallis test.
Effect of parasitema on IL-10 and GM-CSF concentrations stratified by year. Scatter plots with trend line representing the distribution of analytes concentration by parasitemia stratified by year. Only analytes in which parasitemia and year had a significant p-value for the interaction test (before correcting for multiple testing) are shown.
About this article
Cite this article
Aguilar, R., Campo, J.J., Chicuecue, S. et al. Changing plasma cytokine, chemokine and growth factor profiles upon differing malaria transmission intensities. Malar J 18, 406 (2019). https://doi.org/10.1186/s12936-019-3038-x