How does long-term malaria exposure shape baseline immune states in healthy adults?

Four Kenyan sites, one within-country comparison

The paper asks whether long-term residence in a high-malaria-transmission environment leaves a measurable imprint on the baseline immune state of adults who are not currently sick. To reduce the confounding that comes with comparing Africans to Europeans, the authors recruited all 60 participants inside Kenya, from four sites chosen for different Plasmodium falciparum transmission intensity: Ngerenya (low to negligible), Nairobi (minimal to no transmission), Junju (intermediate) and Ahero (high transmission). Fifteen adults came from each site, all clinically healthy and all with at least 10 years of residence in their region. The cohort was predominantly male (68%) with a median age of 26 years.

Exposure was not defined by a single test. The authors combined regional transmission history with anti-P. falciparum IgG responses and serological profiles against eight additional common viral and parasitic pathogens, including Epstein-Barr virus, cytomegalovirus, HSV-1, adenovirus, Ascaris lumbricoides, Schistosoma mansoni, Strongyloides stercoralis and Entamoeba histolytica. This matters for reading every later result: the comparison is between low and high cumulative pathogen exposure, with malaria transmission as the organising gradient rather than an isolated variable.

Antibody profiles separate the exposure environments

Site-level heatmaps showed broader increases across multiple pathogen-specific IgG responses in the higher-malaria-transmission sites, Junju and Ahero, relative to Nairobi and Ngerenya. Principal component analysis of the IgG profiles showed partial grouping of lower- and higher-exposure sites, with some overlap between categories. P. falciparum-specific IgG differed across regions, with the largest contrast between Ahero and Nairobi (FDR = 4.37 × 10⁻⁷).

That overlap is part of the finding, not a flaw to skip. Participants were then analysed as low-exposure (n = 30) and high-exposure (n = 30) groups, and the serology shows the grouping is a gradient of cumulative exposure rather than two cleanly separated populations. The authors also note that the relatively low helminth seroprevalence in the cohort is consistent with reduced helminth transmission in Kenya following sustained mass drug administration, so the measured pathogen panel does not capture every possible exposure difference between sites.

658 genes shift with higher pathogen exposure

The core molecular measurement was bulk RNA sequencing of peripheral blood mononuclear cells (PBMCs), the mixed white-cell fraction of blood, combined with a 65-plex plasma immunoassay for circulating immune mediators. After filtering, 14,576 genes entered differential expression analysis using edgeR, with a Benjamini-Hochberg false discovery rate threshold below 0.05 and library preparation batch included as a covariate. The primary model compared high versus low pathogen exposure and identified 658 differentially expressed genes.

Age and sex were not in the primary model: sex distribution did not differ between exposure groups, and age was treated as a proxy for cumulative exposure. This is a deliberate choice with a cost. Sensitivity analyses adding age, sex, or both reduced the number of significant genes but largely preserved the core exposure-associated signature, which is why the authors kept the unadjusted model as primary. Readers should treat the 658-gene count as the size of the unadjusted signal, not as a fixed effect size that survives every adjustment.

Adaptive activation rises while inflammatory programmes fall

The direction of the shift is the paper's most distinctive result. Higher exposure was associated with increased expression of pathways related to B-cell activation, antigen presentation, immune regulation and post-transcriptional RNA processing. In parallel, inflammatory, chemotactic and vascular gene programmes were reduced. The authors describe this as adaptive immune activation combined with restrained inflammatory transcription, rather than generalised immune activation.

Weighted gene co-expression network analysis (WGCNA) grouped genes into modules and tested their association with exposure status. Three modules illustrate the pattern: GM17 centred on B-cell activation and receptor signalling around CD19 and including B-cell hub genes BANK1 and MS4A1; GM20 capturing RNA processing and splicing machinery around HNRNPA1; and GM30 enriched for inflammatory mediators centred on IL1B and IL6. The module structure is the paper's evidence that these are coordinated programmes rather than scattered gene changes, and it is a statistical association with exposure group, not a demonstration that exposure caused the change.

Circulating mediators rise despite lower inflammatory transcription

A reader could reasonably expect that reduced inflammatory gene expression would come with lower circulating inflammatory proteins. The paper reports the opposite pattern: circulating immune mediators were modestly elevated despite the reduced inflammatory gene programmes. The authors call this a partial dissociation between systemic mediator levels and cellular transcription.

This is a measurement distinction worth keeping separate. The gene expression readout comes from PBMC transcripts, while the mediator readout comes from plasma protein concentrations in pg/mL measured on the Luminex platform. They answer different questions about different compartments and different molecular layers, so one cannot be substituted for the other. The paper's interpretation is that a dampened inflammatory transcriptional baseline in circulating cells can coexist with modestly raised plasma mediators, which is why the authors frame the state as recalibration rather than simple suppression.

An integrated score separates the two exposure groups

To test whether the separate readouts form a coherent state, the authors standardised the selected gene modules and computationally inferred immune cell populations to Z-scores and combined them with LASSO logistic regression, using low versus high exposure as the binary outcome. The resulting integrated immune state score differed between groups with a two-sided Wilcoxon rank-sum P value of 1.08 × 10⁻¹¹.

Two cautions belong next to that number. First, the score was built and evaluated on the same 60 participants, 30 per group, so it demonstrates separation within this cohort rather than validated predictive performance in a new population. Second, the immune cell proportions are computationally inferred from bulk RNA sequencing using deconvolution methods (ENIGMA with the LM22 signature and an independent ABIS replication), not measured by flow cytometry. The score is best read as evidence that the transcriptional, cellular and mediator features move together, which the paper supports with a Spearman correlation matrix among integrated features.

Malaria sits inside broader endemic exposure

The study design was built around malaria transmission, but the authors are explicit that the gradient also reflects wider ecological variation in exposure to co-endemic pathogens. Because malaria transmission formed the basis of site selection and was the most pronounced exposure gradient, they argue malaria is likely a key contributor to the observed immune variation, operating within cumulative exposure to multiple endemic pathogens. The paper therefore describes the result as environmental immune adaptation rather than a malaria-specific effect.

This scope condition changes how the finding should be used. The 658 genes, the module directions and the integrated score describe adults in these four Kenyan settings; they do not establish which pathogen drove which pathway, and they do not show that malaria alone produces this state. The authors connect the pattern to prior reports of adaptive immune expansion in high-endemic settings and of dampened inflammatory monocyte phenotypes in malaria-exposed adults, but those comparisons are context for interpretation, not additional evidence generated here.

A baseline state that may shape later responses

The paper's forward-looking claim is conditional: this distinct baseline immune state may influence responses to subsequent infections and vaccination. That is a hypothesis about consequences, not a measured outcome. Nothing in the study tests how these participants responded to an infection or a vaccine, so the result should not be read as evidence that high-exposure adults will respond worse or better.

The most useful takeaway is a comparison rather than a single number. Among these 60 healthy adults, higher cumulative pathogen exposure tracked with upregulated B-cell and antigen-presentation programmes and downregulated inflammatory, chemotactic and vascular transcription, while plasma mediators moved modestly in the opposite direction from the cellular inflammatory signal. Anyone assessing immune baseline data from endemic settings should therefore expect adaptive activation and inflammatory restraint to appear together, and should not treat a low inflammatory transcript profile as proof of a low-mediator systemic state.

Sources used in this answer