Abstract
Introduction:
Extensive Plasmodium falciparum genetic diversity plays a role in immune evasion, and antibody responses can be strain-specific or broadly reactive depending on the epitope. Controlled human malaria infection (CHMI) allows investigation of immune responses to variant parasite proteins after a single infection with a known strain.
Methods:
We designed a novel diversity-reflecting peptide microarray containing 638,817 unique peptides representing 22,655 variants of 227 proteins from 23 P. falciparum genome assemblies and 379 field isolates. Using this array, we probed sera from 38 malaria naïve adults before and 28 days after CHMI with one of two genetically distinct P. falciparum strains, NF54 (n = 21) or 7G8 (n = 17). We examined fold-increase in antibody response (intensity) and cross-reactivity to protein variants (breadth). ABCPred was used to predict linear epitopes for all 227 proteins. We used MEME to identify enriched motifs in regions of high intensity or breadth, which were presumed to be potential epitopes.
Results:
While the two CHMI groups had similar intensity of responses to all proteins on the array, 20 proteins on the array had differential breadth of responses and participants infected with 7G8 strain had a higher breadth of responses to 17 of them. Of 543 ABCPred-predicted epitopes, 66 overlapped with MEME-identified epitopes, six of which were highly cross-reactive with >95% of peptide variants serorecognized by at least one CHMI group.
Discussion:
Overall, we found most antibody responses to be comparable after infection with the NF54 strain or 7G8 strain, but we saw notable differences for ~10% of proteins on the array. While many MEME-identified epitopes from highly cross-reactive proteins were asparagine rich, an epitope from PF3D7_1033200 (ETRAMP10.2) was not. Highly cross-reactive responses to ETRAMP10.2 could be further characterized and ETRAMP10.2 could be considered for inclusion in a next generation vaccine.
1 Introduction
Malaria continues to cause substantial morbidity and mortality worldwide, with over 263 million cases and 597,000 deaths in 2023 attributable to malaria, despite decades of research and prevention efforts (). RTS,S/AS01 is now recommended by the World Health Organization (WHO) for widespread use in countries with moderate to high transmission of Plasmodium falciparum (). But the four-dose regimen poses a logistical challenge, and supply is limited (, ). Also, RTS,S-associated efficacy is modest and wanes over time (, ). Another vaccine, R21/Matrix M, is also recommended by the WHO, which may alleviate supply chain issues (). However, the R21/Matrix M is also given as a four-dose regimen, and the long-term efficacy and required booster frequency remain unclear (, ). Additional and longer lasting malaria vaccines are needed to achieve WHO goals for malaria vaccines, including >90% efficacy lasting >12 months (). Two major impediments to a highly efficacious malaria vaccine are the complex P. falciparum life cycle and its extensive genetic diversity.
Vaccine candidates that target antigens from each phase of the P. falciparum life cycle are currently under development and may be candidates for incorporation into multi-stage vaccines (). But thousands of additional uncharacterized proteins exist with uncertain timing of expression and propensity to elicit protective immunity. Several challenges exist to identify novel candidates for P. falciparum vaccines. Commonly used assays, including ELISA, generally assess only a few antigens at a time, but the P. falciparum genome encodes over 5,000 proteins (). Additionally, methods such as ELISA and protein microarrays usually examine responses to whole proteins or large segments of proteins. Some malaria proteins have repetitive regions that elicit high antibody levels that correlate more with exposure to malaria than protection from infection or disease, and these immunodominant repetitive regions will drive the measured signal in an ELISA (). Focusing on the most immunodominant region, which is effectively represented when evaluating one signal per antigen, can hinder detection of signals corresponding to informative antibody binding to less immunogenic epitopes. One example is the circumsporozoite protein (CSP), which has both an immunodominant central repeat region that induces high seroreactivity, and additional epitopes, including the binding sites of the CIS43 and L9 monoclonal antibodies (mAbs), with potent neutralizing activity (, ). Studying humoral responses to overlapping peptides that represent the full length of a protein permits differentiation of responses to multiple epitopes, including subdominant epitopes.
The extensive genetic diversity that is a hallmark of P. falciparum surface proteins facilitates host immune evasion (–). Determining the immune response to single antigenic variants, as occurs with most antibody-based assays, represents only partial information that does not consider locally circulating parasite strains or potential cross-reactivity between strains. Previous studies have used protein microarrays to study natural and vaccine-induced antibody responses to conveniently available laboratory strains but have not systematically assessed antigenic diversity in field isolates (–). Allele-specific vaccine efficacy demonstrated in other studies underscores the need to rationally design vaccines effective against antigenically diverse malaria parasites in endemic areas (–). Depending on the antigen and epitopes studied, antibody responses to P. falciparum can be broadly reactive, allele-specific, or a combination of both (i.e., broad reactivity and allele-specific functional activity). Ideally, a malaria vaccine would be broadly reactive to variants of essential parasite epitopes or would target conserved regions. Better knowledge of which proteins elicit broadly reactive versus allele-specific antibody responses could, therefore, inform development of a broadly protective malaria vaccine.
To address the challenges of finding subdominant epitopes and measuring cross-reactivity, we use high-density peptide microarrays to query diverse sequences from different strains/isolates and identify putative epitopes with amino acid-level precision. With our method, individual proteins from reference strains and field isolates are represented by 16-amino acid-long peptides that span the protein with a high degree of overlap (e.g., 15 amino acid overlap). In the current study, we measured the antibody responses elicited by single malaria infections in participants who were experimentally infected with malaria, i.e., received a controlled human malaria infection (CHMI). To disentangle strain-specific vs. cross-reactive antibody responses at the epitope level, we used high-density peptide microarrays to measure antibody responses to 227 P. falciparum proteins inclusive of antigenic variants in two groups of participants who received CHMI with different P. falciparum strains: NF54, thought to be of West African origin, and 7G8, a Brazilian clone (). To our knowledge, these results represent the most comprehensive study to date of strain-specific and cross-reactive responses after a single episode of malaria with a known P. falciparum strain.
2 Materials and methods
2.1 Study participants
Serum samples originated from two CHMI dose optimization studies, NCT01546389 and NCT02780154 (, ). Included participants received either NF54 (n=21) or 7G8 (n=17) parasites, developed blood-stage malaria infection, received standard of care malaria treatment, and consented to future use. We probed sera from baseline (Day 1) and post-CHMI (Day 29) from these 38 participants on our diversity-reflecting peptide microarray (Table 1, Supplementary Table S1). The Institutional Review Board at the University of Maryland, Baltimore, approved the research protocols for both clinical trials, and appropriate informed consent for future use of specimens was obtained from study participants. The studies from which the samples were collected were conducted in accordance with the Declaration of Helsinki principles.
Table 1
| Number of participants | Route1 | Dose (sporozoites) | Strain | Method of Diagnosis2 | Treatment3 |
|---|---|---|---|---|---|
| 3 | ID | 10,000 | NF54 | TBS | CQ |
| 3 | ID | 10,000 | NF54 | TBS | CQ |
| 4 | ID | 10,000 | NF54 | TBS | CQ |
| 2 | ID | 10,000 | NF54 | TBS | CQ |
| 2 | ID | 50,000 | NF54 | TBS | CQ |
| 2 | ID | 50,000 | NF54 | TBS | CQ |
| 3 | DVI | 800 | 7G8 | uPCR | A/P |
| 4 | DVI | 1,600 | 7G8 | uPCR | A/P |
| 8 | DVI | 3,200 | 7G8 | uPCR | A/P |
| 2 | DVI | 4,800 | 7G8 | uPCR | A/P |
| 5 | DVI | 3,200 | NF54 | uPCR | A/P |
Details of controlled human malaria infection procedures for participants included in this study.
1ID, intradermal injection, DVI, direct venous injection.
2TBS, thick blood smear; uPCR, ultrasensitive polymerase chain reaction.
3CQ = 600 mg chloroquine base at time 0, 300 mg chloroquine base at hours 6, 24, and 48, A/P = atovaquone/proguanil 1000/400 mg once daily for 3 days.
ID group = Lyke KE et al. Optimizing intradermal administration of cryopreserved Plasmodium falciparum sporozoites in controlled human malaria infection. Am J Trop Med Hyg. 2015;93(6):1274-1284.
DVI group = Laurens et al. Dose-dependent infectivity of aseptic, purified, cryopreserved Plasmodium falciparum 7G8 sporozoites in malaria-naïve adults. J Infect Dis. 2019;220:1962-6.
2.2 Diversity-reflecting peptide microarray design
We designed a diversity-reflecting peptide array that included peptides mapping to a total of 22,655 protein variants encoded by 227 genes. The encoded proteins included known surface-exposed proteins, proteins previously observed to bind antibodies following recent malaria exposure, gametocyte and sexual stage proteins, proteins involved in metabolic pathways, and proteins included on previous peptide microarrays, with considerable overlap between these groupings. The variants originated from 23 de novo genome assemblies from East and West Africa and Southeast Asia (Ifeonu, Moser et al., in preparation), and alleles reconstructed from whole genome sequencing data generated in-house from 379 samples from Mali, Malawi, Brazil, Cambodia, Myanmar, Thailand, Laos, and Guinea, as well as reconstructed alleles from publicly available sequences from similar locations (Supplementary Table S2) (, ). Each protein was represented by 16-amino acid-long peptides that span the length of the protein, with a 15 amino acid overlap between consecutive peptides. A total of 638,817 unique peptides were included. We also included an additional 10,175 peptides comprised of random amino acids adapted from a collection previously used to distinguish antibody responses to influenza vaccines in triplicate (). Details of in situ peptide synthesis, quality control, serum probing, and slide imaging can be found in the Supplementary Materials.
2.3 Quality control and identification of outliers
We visualized distributions of raw peptide log2-transformed fluorescence intensities across all samples using boxplots and cumulative distribution functions. We inspected the data for globally outlying samples using hierarchical clustering, multidimensional scaling, and principal component analysis.
2.4 Reference protein sequence alignments and peptide filtering
Proteins from the P. falciparum NF54 strain were used as reference for peptide mapping. For each protein variant, we mapped peptides against the reference protein sequence using ssearch36 (Version 36.3.8g). We excluded any peptides for which the Smith-Waterman alignment score was less than three standard deviations below the mean, the alignment length was two or more amino acids shorter than the full-length peptide, or the alignment had more than one gap in the reference or the query. Peptides and protein variants that did not have a protein mapping were excluded from the final analysis. After filtering, 520,015 (81.4%) of 638,817 peptides corresponding to 227 proteins had one or more primary reference protein mapping, including >99% of peptides that originated from NF54 or 7G8.
2.5 Data normalization
We executed median normalization of study data using intensities collected for 10,175 randomer peptide sequences included on the arrays. For each sample, the sample-specific median for the 10,175 randomer log2 fluorescence intensities was subtracted from the global median, which was the median of all sample-specific medians, to create a scaling factor for each sample. Thus, each sample specific randomer median was centered to the global randomer median FI. The full data set was then normalized by subtracting the sample-specific scaling factor from each non-randomer peptide log2 intensity. We tested other methods of data normalization to choose the least biased method - additional details of the other methods of normalization and scaling procedures are in the Supplementary Materials.
2.6 Sliding window smoothing procedure
Raw peptide array fluorescence intensities were log2 transformed. As consecutive peptides overlapped by about 94%, we expected intensities to be correlated as they may contain the same epitope. To increase the data signal-to-noise ratio, we applied a sliding window-based average smoothing procedure (). Each peptide was represented by its midpoint (end position minus beginning position divided by two). Each amino acid position along the sequence was represented by a smoothed log2 fluorescence intensity equal to the average log2 fluorescence intensity for all peptides with a midpoint four amino acids before or after the given position. For each participant-sample and protein variant, we applied the sliding window smoothing procedure to calculate the average fluorescence intensity corresponding to each amino acid position on the reference protein sequence.
2.7 Serorecognition
We defined serorecognition thresholds for each peptide as the mean plus 2.5 times the standard deviation of the log2 fluorescence intensity at Day 1 from all participants (regardless of CHMI group). We classified a peptide as serorecognized by an individual participant if its log2 fluorescence intensity was greater than the serorecognition threshold for that peptide. We then applied a sliding window smoothing procedure so that at a given amino acid position along each protein, the number of serorecognized peptides with mid-points within four amino acids before or after that position were summed. We then calculated the sum of the smoothed log2 fluorescence intensities (SSI) across each protein variant for each sample and used it to define the serorecognition threshold for a protein variant (mean plus 2.5 times the SSI standard deviation for all baseline samples). For each post-CHMI sample, we defined a protein variant as serorecognized if the Day 29 SSI exceeded the serorecognition threshold for that protein variant. On the cohort level (CHMI group), we defined each feature (peptide and protein variant) as serorecognized by a CHMI group if the log2 fluorescence intensity exceeded the serorecognition threshold for >15% participants in that CHMI group.
To describe cross-reactivity (breadth of antibody response), for each protein, we captured the mean number of serorecognized peptides for each CHMI group at each amino acid position, and the maximum mean number of serorecognized peptide variants (MMNS) and its corresponding location. The MMNS and its location represents the epitope with the highest cross-reactivity (breadth of antibody response) for each protein and CHMI group.
2.8 Seroreactivity
Using the sliding window procedure, for each CHMI group we calculated the mean log2 fold change between D1 and D29 for each amino acid along the reference protein sequence for each protein variant. For each protein at each position, we averaged the log2 fold changes for all included protein variants, and calculated each CHMI group mean (mean of the mean of log2 fold changes, referring to the mean of the protein variants and the mean of the individuals in each cohort). We then summed all the position-specific mean of the mean log2 fold changes for each protein to represent mean (amongst protein variants and for each cohort) area under the curve for each protein (MF-AUC). We ascertained the amino acid positions and values for maximum mean of the mean fold changes (MMMF) for each protein by CHMI group. The MMMF and its location represents the epitope with the highest antibody binding (antibody intensity) for each protein and CHMI group.
2.9 Identification and comparison of serorecognized and seroreactive protein variants between the two CHMI groups
We assessed differences in the numbers of peptides and protein variants serorecognized after exposure to NF54 vs. 7G8 strain CHMI on the peptide and protein variant level using a Fisher’s exact test. To evaluate seroreactivity differences between the two CHMI groups, we implemented an ANOVA model in the limma R package fitted to MF-AUC including a fixed effect for CHMI group (NF54 and 7G8 strains) (). We used a contrast to assess statistical significance of the mean MF-AUC difference between NF54 and 7G8 (H0: mean NF54 – mean 7G8 = 0, H1: mean NF54 - mean 7G8 ≠ 0, on the log2 scale).
We adjusted for multiple comparisons using the Benjamini-Hochberg procedure (as implemented in the R p.adjust function). We considered any responses with a false discovery rate (FDR) adjusted p-value < 0.05 to be proteins with differential serorecognition or seroreactivity.
2.10 Identification and comparison of epitopes
We used MEME (Version 5.0.5) to find dominant motifs that could represent epitopes in regions that overlapped with areas of high serorecognition or seroreactivity (MMNS or MMMF), which we then termed “serologically identified epitopes” (). For each protein, at positions of interest (MMNS or MMMF), we used the peptide sequence per protein variant with the highest mean fluorescence intensity value at Day 29 across participants separately for each CHMI cohort as input for identifying motifs. We set the number of motifs detected per sequence to 0 or 1 occurrence (-mod zoops option) and the minimum and maximum motif width to 8 and 16, respectively (-minw and -maxw and options, respectively). The background sequence collection used for the null model included the unique set of all peptide sequences that did not overlap with either MMNS or MMMF regions. Motif discovery was analyzed separately for each CHMI group (NF54 and 7G8 strains). We identified shared motifs/epitopes between CHMI groups using the TOMTOM (Version 5.0.5) motif comparison software ().
2.11 Assessment of serorecognition and seroreactivity for a set of known epitopes
We assessed serorecognition and seroreactivity for a set of literature-curated epitopes (Supplementary Table S3) based on a literature search and computationally predicted B-cell epitopes produced by the ABCPred prediction algorithm (Supplementary Table S4) (). For ABCPred, we used a score cut off of >0.9 to define likely B-cell epitopes. We assessed the overlap between literature-curated or ABCPred-predicted and serologically identified epitopes using the TOMTOM motif comparison tool ().
2.12 Summarization of epitope results across individual participants
Average epitope MMNS and MMMF for individual participants were visualized for literature-curated epitopes and serologically identified epitopes using heatmaps utilizing (1) a data driven approach (hierarchical clustering of epitopes and participants) and a (2) pre-ordered approach (participants sorted by CHMI group) to identify similar or divergent serorecognition/seroreactivity clustering profiles between CHMI groups.
2.13 Software
Data was analyzed using R (Version 3.6.0, 26APR2019) and R Bioconductor packages on the Ubuntu (Version 13.04) operating system.
3 Results
3.1 Quality control
The two CHMI groups (NF54 and 7G8) did not cluster into mutually exclusive groups by Euclidean distances, principal component analysis (PCA) or 1-Spearman correlation distance, and there were no strong outliers (Supplementary Figures S1, S2). Samples were run in four batches, and CHMI groups and time points were well balanced between batches. Samples did not cluster by batch in Euclidean distances, PCA or 1-Spearman correlation distance plots (Supplementary Figure S3). Overall, inspection of sample-specific intensities across the 638,817 peptides showed that NF54 and 7G8 CHMI groups had similar pre-CHMI median intensities, but the 7G8 CHMI group tended to have higher median intensities post-CHMI and higher fold changes relative to pre-CHMI compared to the NF54 CHMI group (aggregated by protocol and group in Figure 1A). Therefore, we explored different methods of normalization to mitigate bias, including spatial and background correction, blank spot normalization, and randomer normalization (Figures 1B-D). Randomer normalization reduced systematic differences between the groups optimally, so we used randomer normalized data for our data analysis. Boxplots of individual samples and timepoints after randomer normalization can be visualized in Supplementary Figure S4.
Figure 1
3.2 Serorecognition and seroreactivity
First, we assessed the breadth of responses based on serorecognition to all peptides on the array (across all proteins and variants) for the two CHMI groups. Each peptide was classified according to the strains and geographic locations in which they were originally identified, and peptides could belong in multiple groups (Table 2, first column). For every comparison that we evaluated, we observed higher serorecognition (breadth of responses) in the NF54 CHMI group compared to the 7G8 CHMI group (odds ratios [OR] 1.68-2.84, p <0.0001 for all, Table 2).
Table 2
| Strain/Geographical location | Total number of peptides | Number of serorecognized peptides (7G8 CHMI group) | Number of serorecognized peptides (NF54 CHMI group) | Odds ratio (95% CI) [NF54 vs. 7G8] | P-value |
|---|---|---|---|---|---|
| NF54 (W. Africa) | 278,643 | 5,979 | 10,719 | 1.82 (1.77,1.88) | <0.0001* |
| 7G8 (Brazil) | 258,175 | 5,658 | 10,205 | 1.84 (1.78,1.9) | <0.0001* |
| NF54-specific | 20,279 | 344 | 751 | 2.24 (1.96,2.55) | <0.0001* |
| 7G8-specific | 7,832 | 267 | 712 | 2.84 (2.45,3.29) | <0.0001* |
| Mali | 421,891 | 10,840 | 24,359 | 2.32 (2.27,2.38) | <0.0001* |
| Malawi | 458,701 | 12,271 | 27,473 | 2.32 (2.27,2.37) | <0.0001* |
| Brazil | 286,233 | 6,428 | 12,131 | 1.93 (1.87,1.99) | <0.0001* |
| Cambodia | 333,259 | 8,190 | 16,074 | 2.01 (1.96,2.07) | <0.0001* |
| Myanmar | 297,953 | 6,958 | 13,514 | 1.99 (1.93,2.05) | <0.0001* |
| Thailand | 275,717 | 6,317 | 12,054 | 1.95 (1.89,2.01) | <0.0001* |
| Laos | 259,996 | 5,778 | 10,383 | 1.83 (1.77,1.89) | <0.0001* |
| Guinea | 252,627 | 5,340 | 9,448 | 1.8 (1.74,1.86) | <0.0001* |
| Conserved across >1 locations | 348,633 | 8,678 | 18,473 | 2.19 (2.14,2.25) | <0.0001* |
| Conserved across all locations | 225,394 | 4,722 | 8,028 | 1.73 (1.66,1.79) | <0.0001* |
| Conserved across all strains and locations | 208,126 | 4,229 | 6,991 | 1.68 (1.61,1.74) | <0.0001* |
| Any location | 629,618 | 17,483 | 39,871 | 2.37 (2.32,2.41) | <0.0001* |
| Any strain or location | 638,817 | 17,574 | 40,001 | 2.36 (2.32,2.41) | <0.0001* |
Serorecognition (breadth of responses) of peptides from reference strains and field isolates.
For every comparison that we evaluated, serorecognition of peptides was higher in the NF54 CHMI group compared to the 7G8 CHMI group. * indicates P-value < 0.05.
Next, we investigated serorecognition on a protein variant level, i.e., the numbers of protein variants per protein that were serorecognized by each CHMI group. Overall, we observed higher protein variant serorecognition (breadth of responses) by the 7G8 CHMI group than the NF54 CHMI group (OR 0.42, p<0.001 for protein variants from any strain or location, Table 3). This finding was statistically significant for the comparisons with more protein variants, i.e., most of the country specific comparisons (ORs 0.39-0.47). Of the 22,656 protein variants from 227 proteins on the array, only 364 and 347 protein variants were unique to the NF54 strain and the 7G8 strain, respectively. So, the majority of serorecognition of protein variants seen was related to cross-reactivity. No variant of the 227 proteins probed on the array was conserved across all strains and locations.
Table 3
| Strain/Geographical location | Total number of protein variants | Number of serorecognized protein variants (7G8 CHMI group) | Number of serorecognized protein variants (NF54 CHMI group) | Odds ratio (95% CI) [NF54 vs. 7G8] | P-value |
|---|---|---|---|---|---|
| NF54-specific (W. Africa) | 364 | 22 | 12 | 0.53 (0.24,1.14) | 0.1127 |
| 7G8-specific (Brazil) | 347 | 22 | 11 | 0.48 (0.21,1.06) | 0.073 |
| Mali | 7,373 | 572 | 275 | 0.46 (0.4,0.53) | <0.0001* |
| Malawi | 9,813 | 805 | 357 | 0.42 (0.37,0.48) | <0.0001* |
| Brazil | 1,067 | 73 | 35 | 0.46 (0.3,0.71) | 0.0002* |
| Cambodia | 3,620 | 321 | 157 | 0.47 (0.38,0.57) | <0.0001* |
| Myanmar | 1,447 | 124 | 51 | 0.39 (0.27,0.55) | <0.0001* |
| Thailand | 1,073 | 100 | 46 | 0.44 (0.3,0.63) | <0.0001* |
| Laos | 400 | 29 | 12 | 0.4 (0.18,0.81) | 0.0095 |
| Guinea | 252 | 17 | 10 | 0.57 (0.23,1.35) | 0.2348 |
| Conserved across >1 locations | 1,735 | 131 | 86 | 0.64 (0.48,0.85) | 0.002* |
| Conserved across all locations | 78 | 3 | 1 | 0.33 (0.01,4.17) | 0.6201 |
| Conserved across all strains and locations | 0 | 0 | 0 | 0 (0,0) | >0.9999 |
| Any location | 21,872 | 1,822 | 799 | 0.42 (0.38,0.45) | <0.0001* |
| Any strain or location | 22,656 | 1,871 | 827 | 0.42 (0.39,0.46) | <0.0001* |
Serorecognition (breadth of responses) of protein variants from reference strains and field isolates.
Serorecognition of protein variants was higher in the 7G8 CHMI group, except for protein variants specific to NF54, specific to 7G8, specific to Guinea, or conserved across all locations, which were similar in the two CHMI groups. * indicates P-value < 0.05.
The 7G8 CHMI group serorecognized more protein variants than the NF54 CHMI group for 17 proteins and the NF54 CHMI group serorecognized more protein variants for three proteins based on a Fisher’s Exact Test (Figure 2). Of the 17 proteins with more serorecognized variants by the 7G8 CHMI group, the 7G8 CHMI group serorecognized 8–230 more protein variants than the NF54 CHMI group (7-100% of total variants). For the three proteins with more serorecognized variants by the NF54 CHMI group, the NF54 CHMI group serorecognized 126–187 more protein variants than the 7G8 CHMI group (90 – 97% of total variants).
Figure 2
When comparing antibody levels elicited by CHMI (MF-AUC, a measure of aggregate signal across each protein), both CHMI groups showed similar seroreactivity to all 227 proteins on the array (Supplementary Table S5).
3.3 Identification and comparison of serologically identified epitopes
For each CHMI group, along each protein sequence, we identified one epitope that had the maximum mean number of serorecognized peptide variants (MMNS), implying the largest breadth of responses along that protein, and one epitope that had the maximum mean of the mean fold changes (MMMF), implying the highest intensity responses along that protein. We termed these “serologically identified epitopes.” For each CHMI group, any peptide sequences mapping to the MMNS (breadth) or MMMF (intensity) were used as input for MEME to find dominant motifs in each of the 227 proteins. Using MEME, we found 16 and 14 serologically identified epitopes for the 7G8 and NF54 CHMI groups, respectively, based on breadth of responses (MMNS) (Supplementary Table S6). Based on TOMTOM analysis, seven of these serologically identified epitopes overlapped between the two CHMI groups with one 7G8 epitope overlapping with two NF54 epitopes (Supplementary Table S6). Heatmaps of individual responses segregated by CHMI group are displayed in Figure 3. While some of the overlapping epitopes clustered by Euclidean distances, some did not. Responses did not cluster by CHMI group and responses were overall heterogeneous between participants. Three of the overlapping serologically identified epitopes mapped to single proteins for both groups: PF3D7_1021700, PF3D7_1033200, PF3D7_1325900. The others mapped to 3–67 proteins for each group. Many of the serologically identified epitopes found are highly repetitive and contain multiple asparagine residues.
Figure 3
Using the second metric, MMMF, which represents the intensity of responses, we found nine and two serologically identified epitopes for the 7G8 and NF54 CHMI groups, respectively (Supplementary Table S7). One epitope was shared among the two CHMI groups based on TOMTOM analysis, which mapped to 65 proteins for the 7G8 serologically identified epitope and 1 protein (PF3D7_0110600) for the NF54 serologically identified epitope. Heatmaps of individual responses segregated by CHMI group are shown in Figure 4. We did not see clustering of the overlapping epitope by Euclidean distances, nor did we see clustering of responses by CHMI group. Serologically identified epitopes found using maximum mean seroreactivity were repetitive and contained many asparagine residues. For all methods of ascertaining serologically identified epitopes, some serologically identified epitopes were found in multiple proteins (up to 70 proteins). These were highly repetitive epitopes with more than 50% of the epitope comprised of asparagine residues.
Figure 4
3.4 Comparison of literature-curated epitopes and ABCPred-predicted epitopes to serologically identified epitopes
We used the TOMTOM motif comparison tool to find similarities between literature-curated epitopes and serologically identified epitopes. We found one 7G8 serologically identified epitope (EPNQANKE) based on intensity of responses (MMMF) that had similarities to two literature-curated epitopes, both in merozoite surface protein 2 (MSP-2, PF3D7_0206800): EPNQANKE in the variable region and ECTDGNKE in the conserved C-terminal region (Supplementary Table S8).
When comparing ABCPred-predicted and serologically identified epitopes with TOMTOM, 65 of the 543 total ABCPred-predicted epitopes from 19 proteins overlapped with serologically identified epitopes from at least one CHMI group based on MNS or MMMF (Supplementary Table S9). Six of these 19 proteins with overlapping epitopes had highly cross-reactive protein variants with ≥99% of protein variants serorecognized by at least one CHMI group: PF3D7_0305300, PF3D7_0726400, PF3D7_0828100, PF3D7_1033200, PF3D7_1205500, PF3D7_1447800 (Table 4, Figure 5). Five were CHMI strain specific: PF3D7_0305300, PF3D7_0726400, PF3D7_0828100, PF3D7_1033200, and PF3D7_1447800, i.e., one CHMI group serorecognized all variants, but the other group serorecognized few or none. One protein, PF3D7_1205500, was highly cross-reactive with 100% of protein variants serorecognized by both CHMI groups.
Table 4
| Protein | Name | Total # of variants on the array | # of protein variants serorecognized by 7G8 CHMI group | # of protein variants serorecognized by NF54 CHMI group | Maximum proportion of variants serorecognized | Lifecycle stage | GO terms |
|---|---|---|---|---|---|---|---|
| PF3D7_0110600 | phosphatidylinositol-4-phosphate 5-kinase | 104 | 3 | 0 | 3% | Sporozoite | calcium ion binding, ATP binding, phosphatidylinositol phosphate kinase activity, small GTPase binding |
| PF3D7_0206800 | MSP2 | 176 | 12 | 0 | 7% | Merozoite | cell adhesion, nucleus, plasma membrane, integral component of plasma membrane, anchored component of plasma membrane |
| PF3D7_0305300 | conserved Plasmodium membrane protein, unknown function | 192 | 5 | 192 | 100% | Unavailable | cell surface, membrane, integral component of membrane |
| PF3D7_0310300 | phosphoglycerate mutase, putative | 170 | 101 | 0 | 59% | Gametocyte | nucleus, cytoplasm, mRNA binding, intramolecular transferase activity |
| PF3D7_0420000 | Zinc finger protein, putative | 201 | 0 | 0 | 0% | Sporozoite, Trophozoite, Gametocyte | nucleus, cytoplasm, membrane, zinc ion binding |
| PF3D7_0726400 | conserved Plasmodium membrane protein, unknown function | 200 | 198 | 0 | 99% | Sporozoite, Gametocyte | DNA repair, integral component of membrane, nuclease activity, 5’-flap endonuclease activity, 5’-3’ exodeoxyribonuclease activity, flap endonuclease activity |
| PF3D7_0731800 | alpha/beta hydrolase, putative | 128 | 0 | 0 | 0% | Unavailable | Membrane, food vacuole, lipase activity |
| PF3D7_0828100 | conserved Plasmodium protein, unknown function | 197 | 197 | 0 | 100% | Unavailable | No data available |
| PF3D7_0909000 | conserved Plasmodium protein, unknown function | 128 | 13 | 0 | 10% | Gametocyte | No data available |
| PF3D7_0930500 | diacylglycerol kinase, putative (DGK1) | 219 | 22 | 0 | 10% | Sporozoite | Signal transduction, integral component of membrane, kinase activity |
| PF3D7_1021100 | conserved Plasmodium protein, unknown function | 239 | 0 | 3 | 1% | Unavailable | No data available |
| PF3D7_1021700 | conserved Plasmodium membrane protein, unknown function | 121 | 3 | 0 | 2% | Sporozoite | Integral component of membrane |
| PF3D7_1033200 | early transcribed membrane protein 10.2 | 50 | 50 | 0 | 100% | Trophozoite | Nucleus, symbiont-containing vacuole membrane, Maurer’s cleft |
| PF3D7_1205500 | zinc finger protein | 248 | 248 | 248 | 100% | Merozoite, Gametocyte | Protein ubiquitination, nucleus, cytoplasm, protein binding |
| PF3D7_1229100 | Multidrug resistance-associated protein 2 (MRP2) | 274 | 2 | 0 | 1% | Sporozoite | Schizogony, response to drug, xenobiotic transport, transmembrane transport, integral component of plasma membrane, ATP binding |
| PF3D7_1324600 | conserved Plasmodium protein, unknown function | 110 | 96 | 11 | 87% | Unavailable | Integral component of membrane |
| PF3D7_1325900 | conserved Plasmodium protein, unknown function | 198 | 137 | 0 | 69% | Sporozoite | No data available |
| PF3D7_1343300 | conserved plasmodium protein | 99 | 0 | 0 | 0% | Trophozoite, Gametocyte | Nucleus |
| PF3D7_1447800 | calponin homology domain-containing protein | 230 | 230 | 0 | 100% | Sporozoite, Gametocyte | Protein binding |
Proteins containing ABCPred predicted epitopes that overlapped with serologically identified epitopes from at least one CHMI group.
Figure 5
4 Discussion
Overall, we observed key differences in antibody responses to peptides and protein variants between two groups of malaria-naïve participants infected with either the NF54 or 7G8 strain of P. falciparum via CHMI. Because the participants were malaria-naïve and the strain with which they were infected was known, the cross-reactive responses seen were driven solely by the single malaria episode. By contrast, in studies of participants from malaria endemic countries, breadth of responses is driven by both past exposures and cross-reactivity, and is, therefore, more difficult to interpret in the context of identifying responses that could drive the broad protection needed for a highly effective malaria vaccine. Our analysis of peptide-level antibody responses in malaria-naïve participants will be informative to subsequent studies in malaria endemic populations. For 20 of the 227 proteins, one strain group showed greater breadth (cross-reactivity) of antibody responses than the other, with the 7G8 CHMI group having higher serorecognition than the NF54 CHMI group for 17 of these 20 proteins. Aggregate signals across each protein, representing the amount of antibody response, did not differ between CHMI groups. This combination of findings is interesting in the context of the differences in study design. Most (76%) of the NF54 CHMI group samples were obtained from a trial where malaria was diagnosed by thick blood smear (TBS), but all of the 7G8 CHMI group samples were obtained from a trial where malaria was diagnosed by ultrasensitive PCR (usPCR, Table 1). Because usPCR becomes positive during malaria infections 1–3 days prior to TBS, and treatment is initiated upon diagnosis by TBS or uPCR, the 7G8 CHMI group participants had presumably shorter exposure to blood stage malaria. The shorter exposure time of the 7G8 CHMI group might therefore be expected to result in lower antibody responses; however, we saw more cross-reactive antibody responses in the 7G8 CHMI group and similar intensity of antibody responses for both CHMI groups. Although it is possible that NF54 has become attenuated over time (), in a small head to head comparison study of 7G8 and NF54, infectivity was similar for the two strains (). Even when considering these alternative explanations, our findings suggest that protein variants from some P. falciparum strains elicit more cross-reactive antibody responses than others.
We identified protein areas of high antibody binding (either breadth or intensity) and used MEME to identify dominant motifs within these areas. The serologically identified epitopes identified by MEME were highly repetitive and contained many asparagine residues. Repetitive proteins are known to be immunogenic and cross-reactive, but it is uncertain whether high antibody responses to these repeat regions are protective in individuals living in malaria endemic areas (, 46). Raghavan et al. showed enrichment of repetitive elements in antibody targets when examining serum from children and adults in high and moderate transmission settings, and higher breadth of responses to repetitive peptides in children in high transmission settings compared to moderate transmission settings (). The moderate and high transmission settings differed both in number of exposures and time between exposure and sample collection, so the authors postulated that the responses to repetitive peptides could be exposure-dependent and/or short-lived. The Raghavan et al. study used the phage-immunoprecipitation sequencing (PhIP-Seq) technology, which results in longer peptides (up to 90 amino acids in length) and can allow for some capture of conformational epitopes but loses the resolution of epitope mapping possible with peptide arrays (, ). Additionally, PhIP-Seq uses serial enrichment and amplification, which can result in potential amplification of non-specific phages and decreased specificity (); whereas the peptide arrays used in the current study are fabricated onto slides such that each antibody-bound peptide has a known sequence. Our study shows that antibodies against repetitive motifs can develop after a single exposure. Follow-up testing in CHMI studies could further elucidate the longevity of responses.
A serologically identified epitope in our study, EPNQANKE, matched two literature-curated epitopes in merozoite surface protein-2 (MSP-2, PF3D7_0206800), EPNQANKE and ECTDGNKE. Prior studies of the association between anti-MSP-2 antibodies and protection from clinical malaria have been conflicting, perhaps due to the use of different methods, including assessment of different IgG subclasses and antigens (–). While EPNQANKE lies within the variable region of MSP-2, ECTDGNKE is in the conserved C-terminal region (55). A study of a mouse monoclonal antibody against the conserved C-terminal region of MSP-2 identified a minimal binding epitope (NKENCGAA) that shares the asparagine-lysine-glutamate residues present in both epitopes, evoking the idea that a single cross-reactive antibody could bind both epitopes identified in the current study (56). The paradigm of cross-reactive epitopes within the same protein is exemplified by CSP. CIS43 and L9, two CSP monoclonal antibodies in clinical development for primary malaria prophylaxis, display cross reactivity with the major (NANP) repetitive motif of the central repeat region (immunodominant) and their primary target (immunoprotective) (, 57). Further work isolating and characterizing monoclonal antibodies targeting MSP-2 could determine whether similar cross-reactivity could be important for protective responses.
Serologically identified epitopes from one or both CHMI groups overlapped with epitopes predicted by ABCPred in 19 proteins. Breadth of responses were high for six of these proteins, PF3D7_0305300 (transporter, putative); PF3D7_0726400 (conserved Plasmodium membrane protein, unknown function); PF3D7_0828100 (conserved Plasmodium protein, unknown function); PF3D7_1033200 (early transcribed membrane protein 10.2, ETRAMP10.2); PF3D7_1205500 (zinc finger protein, putative); and PF3D7_1447800 (calponin homology domain-containing protein, putative), with ≥99% of protein variants serorecognized by one or both CHMI groups. The 7G8 CHMI group had more cross-reactive responses, serorecognizing ≥99% of the variants of four of the six proteins compared to the NF54 CHMI group, which serorecognized none or very few. This result shows that strain choice in the design of malaria vaccines could have significant impact on cross strain protection. Many of the serologically identified epitopes that overlapped with ABCPred predicted epitopes were asparagine rich, so may not be good vaccine candidates (). But a cross-reactive epitope from PF3D7_1033200 (ETRAMP10.2, YPTTGPNPNTHGPPS) was not asparagine-rich and was unique to this protein. ETRAMP10.2 is expressed in the parasitophorous vacuole membrane during the blood stage and may be involved in red blood cell binding (58, 59). The 7G8 CHMI group serorecognized all 50 variants of ETRAMP10.2 on the array. A portion of the ETRAMP10.2 epitope (NPNTHGPP) was a dominant motif identified by MEME based on serorecognition for the NF54 CHMI group. However, NF54 group antibody responses to the whole ETRAMP10.2 protein were below the serorecognition threshold, which suggests that examining whole protein serorecognition to define antigenicity can be less sensitive than studying individual epitopes. This may have implications for which variant is the best choice should ETRAMP10.2 become a vaccine candidate.
These findings highlight the importance of considering strain-specific responses when designing novel malaria vaccine targets. For several proteins on the array, antibodies elicited by the 7G8 CHMI were more cross-reactive with the non-CHMI strain protein variants that than antibodies elicited by the NF54 CHMI, highlighting that the strain used for CHMI had a notable impact on responses. Because of the size and complexity of the P. falciparum genome, it is not possible with current technologies to study humoral responses to all P. falciparum proteins encoded in the genome with comprehensive representation of variants with one assay. Our assay was designed to study a broad range of proteins from different lifecycle stages that were predominantly known surface-exposed proteins. However, we recognize that our results are limited to the proteins selected and we likely missed additional proteins with differential responses between the two groups due to the design of the peptide array. Our results are limited to serological responses and, therefore, do not characterize the contributions of cell mediated immunity. While we recognize that serological responses are only a part of the complex human immune response to malaria, antibody responses have been shown to play an important role in protection from symptomatic malaria and in studies of RTS,S- and R21-induced protection (, , , 60–64). Our participants were malaria-naïve adults at the time of CHMI, as is the standard, so findings may not be generalizable to endemic countries where first malaria exposures occur early in childhood. But, to our knowledge, this is the first study to compare CHMI with the NF54 and 7G8 strains to identify differences in humoral responses to a vast array of diverse proteins. Because these findings describe responses that develop after a single malaria episode and do not provide information on future protection from malaria, further work is needed to determine their clinical relevance. However, we were able to discriminate important epitopes with overlap in known and predicted epitopes. Our approach has the potential to find important strain-specific and cross-reactive responses highly relevant to malaria vaccine development.
Statements
Data availability statement
The datasets generated for this study can be found in the Figshare repository, https://doi.org/10.6084/m9.figshare.28278923.
Ethics statement
The studies involving humans were approved by University of Maryland School of Medicine IRB. The studies were conducted in accordance with the local legislation and institutional requirements and with the Declaration of Helsinki principles. The participants provided their written informed consent to participate in this study.
Author contributions
DF: Conceptualization, Validation, Writing – review & editing, Methodology, Writing – original draft. TJ: Formal analysis, Writing – review & editing, Methodology, Visualization, Conceptualization, Validation. KL: Methodology, Writing – review & editing, Resources, Conceptualization. ML: Methodology, Resources, Writing – review & editing, Conceptualization. JS: Resources, Conceptualization, Writing – review & editing, Methodology. ES: Writing – review & editing, Resources. AO: Writing – review & editing, Resources. OI: Resources, Writing – review & editing. TH: Resources, Writing – review & editing. KM: Writing – review & editing, Resources. CG: Formal analysis, Writing – review & editing, Validation. JG: Formal analysis, Visualization, Conceptualization, Validation, Methodology, Writing – review & editing. SH: Conceptualization, Resources, Writing – review & editing. JP: Methodology, Writing – review & editing, Investigation. RP: Writing – review & editing, Methodology, Investigation. JT: Investigation, Methodology, Writing – review & editing, Resources. GD: Conceptualization, Writing – review & editing. ST: Writing – review & editing. MT: Writing – review & editing, Conceptualization, Methodology. AB: Methodology, Conceptualization, Validation, Writing – review & editing, Writing – original draft.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. This work was supported by NIAID under the Vaccine and Treatment Evaluation Unit (VTEU) contract number HHSN272201300022I (PI: Karen Kotloff, Task Order PI: AB). DF-K was supported by the Burroughs Wellcome Fund/American Society of Tropical Medicine and Hygiene Postdoctoral Fellowship in Tropical Infectious Diseases, and NIAID under K23AI155838. AB was supported by an NIAID K23AI125720 and a Passano Foundation Clinician-Investigator Award. National Institutes of Health (NIH) awards U19 AI110820, R01 AI141900, and a “Dean’s Challenge Award” (an intramural program from the University of Maryland School of Medicine), provided funds to sequence the PfSPZ strains and P. falciparum clinical isolates, and support for KM, TH, and JS. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.
Acknowledgments
We would like to thank the study participants for their time and effort. We also thank our colleagues who helped define the literature-curated epitopes: Albert Zhou, Kyeisha Laurence, Biraj Shrestha, and Noah Ventimiglia. We would like to thank the Division of Microbiology and Infectious Diseases at the National Institute of Allergy and Infectious Diseases (NIAID) for their support and advice.
Conflict of interest
Authors JP, RP, JT were employed by the company Nimble Therapeutics, Inc, which is now a part of AbbVie, Inc. GD is currently an employee of AstraZeneca but was employed by the NIAID at the time of this work. Author SH was employed by the company Sanaria, Inc.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2025.1641280/full#supplementary-material
References
1
World Health Organization. World malaria report 2024: addressing inequity in the global malaria response. World Health Organization, Geneva (2024). Available online at: https://www.who.int/teams/global-malaria-programme/reports/world-malaria-report-2024 (Accessed May 23, 2025).
2
World Health Organization. WHO recommends groundbreaking malaria vaccine for children at risk. World Health Organization, Geneva (2021).
3
SamarasekeraU. Cautious optimism for malaria vaccine roll-out. Lancet. (2021) 398:1394. doi: 10.1016/S0140-6736(21)02237-6
4
WillyardC. The slow roll-out of the world’s first malaria vaccine. Nature. (2022) 612:S48–S9. doi: 10.1038/d41586-022-04343-7
5
RTSS Clinical Trials Partnership. Efficacy and safety of the RTS,S/AS01 malaria vaccine during 18 months after vaccination: a phase 3 randomized, controlled trial in children and young infants at 11 African sites. PloS Med. (2014) 11:e1001685. doi: 10.1371/journal.pmed.1001685
6
RTSS Clinical Trials Partnership. Efficacy and safety of RTS,S/AS01 malaria vaccine with or without a booster dose in infants and children in Africa: final results of a phase 3, individually randomised, controlled trial. Lancet. (2015) 386:31–45. doi: 10.1016/S0140-6736(15)60721-8
7
World Health Organization. WHO recommends R21/Matrix-M vaccine for malaria prevention in updated advice on immunization(2023). Available online at: https://www.who.int/news/item/02-10-2023-who-recommends-r21-matrix-m-vaccine-for-malaria-prevention-in-updated-advice-on-immunization (Accessed October 16).
8
DatooMSNatamaMHSomeATraoreORouambaTBellamyDet al. Efficacy of a low-dose candidate malaria vaccine, R21 in adjuvant Matrix-M, with seasonal administration to children in Burkina Faso: a randomised controlled trial. Lancet. (2021) 397:1809–18. doi: 10.1016/S0140-6736(21)00943-0
9
DatooMSNatamaHMSomeABellamyDTraoreORouambaTet al. Efficacy and immunogenicity of R21/Matrix-M vaccine against clinical malaria after 2 years’ follow-up in children in Burkina Faso: a phase 1/2b randomised controlled trial. Lancet Infect Dis. (2022) 22:1728–36. doi: 10.1016/S1473-3099(22)00442-X
10
World Health Organization. Malaria vaccines: preferred product characteristics and clinical development considerations. Geneva (2022).
11
HeppnerDGJr.KesterKEOckenhouseCFTornieporthNOforiOLyonJAet al. Towards an RTS,S-based, multi-stage, multi-antigen vaccine against falciparum malaria: progress at the Walter Reed Army Institute of Research. Vaccine. (2005) 23:2243–50. doi: 10.1016/j.vaccine.2005.01.142
12
GardnerMJHallNFungEWhiteOBerrimanMHymanRWet al. Genome sequence of the human malaria parasite Plasmodium falciparum. Nature. (2002) 419:498–511. doi: 10.1038/nature01097
13
HoffmanSLOsterCNPloweCVWoollettGRBeierJCChulayJDet al. Naturally acquired antibodies to sporozoites do not prevent malaria: vaccine development implications. Science. (1987) 237:639–42. doi: 10.1126/science.3299709
14
KisaluNKIdrisAHWeidleCFlores-GarciaYFlynnBJSackBKet al. A human monoclonal antibody prevents malaria infection by targeting a new site of vulnerability on the parasite. Nat Med. (2018) 24:408–16. doi: 10.1038/nm.4512
15
TanJSackBKOyenDZenklusenIPiccoliLBarbieriSet al. A public antibody lineage that potently inhibits malaria infection through dual binding to the circumsporozoite protein. Nat Med. (2018) 24:401–7. doi: 10.1038/nm.4513
16
FlanaganKLWilsonKLPlebanskiM. Polymorphism in liver-stage malaria vaccine candidate proteins: immune evasion and implications for vaccine design. Expert Rev Vaccines. (2016) 15:389–99. doi: 10.1586/14760584.2016.1125785
17
ReederJCBrownGV. Antigenic variation and immune evasion in Plasmodium falciparum malaria. Immunol Cell Biol. (1996) 74:546–54. doi: 10.1038/icb.1996.88
18
FerreiraMUda Silva NunesMWunderlichG. Antigenic diversity and immune evasion by malaria parasites. Clin Diagn Lab Immunol. (2004) 11:987–95. doi: 10.1128/CDLI.11.6.987-995.2004
19
DoolanDLMuYUnalBSundareshSHirstSValdezCet al. Profiling humoral immune responses to P. falciparum infection with protein microarrays. Proteomics. (2008) 8:4680–94. doi: 10.1002/pmic.200800194
20
TrieuAKayalaMABurkCMolinaDMFreilichDARichieTLet al. Sterile protective immunity to malaria is associated with a panel of novel P. falciparum antigens. Mol Cell Proteomics. (2011) 10:M111 007948. doi: 10.1074/mcp.M111.007948
21
FelgnerPLRoestenbergMLiangLHungCJainAPabloJet al. Pre-erythrocytic antibody profiles induced by controlled human malaria infections in healthy volunteers under chloroquine prophylaxis. Sci Rep. (2013) 3:3549. doi: 10.1038/srep03549
22
DentAENakajimaRLiangLBaumEMoormannAMSumbaPOet al. Plasmodium falciparum protein microarray antibody profiles correlate with protection from symptomatic malaria in Kenya. J Infect Dis. (2015) 212:1429–38. doi: 10.1093/infdis/jiv224
23
GandhiKTheraMACoulibalyDTraoreKGuindoABOuattaraAet al. Variation in the circumsporozoite protein of Plasmodium falciparum: vaccine development implications. PloS One. (2014) 9:e101783. doi: 10.1371/journal.pone.0101783
24
OuattaraANiangalyAAdamsMCoulibalyDKoneAKTraoreKet al. Epitope-based sieve analysis of Plasmodium falciparum sequences from a FMP2.1/AS02(A) vaccine trial is consistent with differential vaccine efficacy against immunologically relevant AMA1 variants. Vaccine. (2020) 38:5700–6. doi: 10.1016/j.vaccine.2020.06.035
25
NeafseyDEJuraskaMBedfordTBenkeserDValimCGriggsAet al. Genetic diversity and protective efficacy of the RTS,S/AS01 malaria vaccine. N Engl J Med. (2015) 373:2025–37. doi: 10.1056/NEJMoa1505819
26
OuattaraADwivediAAdamsMNiangalyALaurensMBNyuntMMet al. An in silico analysis of malaria pre-erythrocytic-stage antigens interpreting worldwide genetic data to suggest vaccine candidate variants and epitopes. Microorganisms. (2022) 10:1090. doi: 10.3390/microorganisms10061090
27
BaileyJABerryAATravassosMAOuattaraABoudovaSDotseyEYet al. Microarray analyses reveal strain-specific antibody responses to Plasmodium falciparum apical membrane antigen 1 variants following natural infection and vaccination. Sci Rep. (2020) 10:3952. doi: 10.1038/s41598-020-60551-z
28
BaileyJAPabloJNiangalyATravassosMAOuattaraACoulibalyDet al. Seroreactivity to a large panel of field-derived Plasmodium falciparum apical membrane antigen 1 and merozoite surface protein 1 variants reflects seasonal and lifetime acquired responses to malaria. Am J Trop Med Hyg. (2015) 92:9–12. doi: 10.4269/ajtmh.14-0140
29
TakalaSLPloweCV. Genetic diversity and malaria vaccine design, testing and efficacy: preventing and overcoming ‘vaccine resistant malaria’. Parasite Immunol. (2009) 31:560–73. doi: 10.1111/j.1365-3024.2009.01138.x
30
TakalaSLCoulibalyDTheraMABatchelorAHCummingsMPEscalanteAAet al. Extreme polymorphism in a vaccine antigen and risk of clinical malaria: implications for vaccine development. Sci Transl Med. (2009) 1:2ra5. doi: 10.1126/scitranslmed.3000257
31
DuttaSLeeSYBatchelorAHLanarDE. Structural basis of antigenic escape of a malaria vaccine candidate. Proc Natl Acad Sci U.S.A. (2007) 104:12488–93. doi: 10.1073/pnas.0701464104
32
OuattaraATakala-HarrisonSTheraMACoulibalyDNiangalyASayeRet al. Molecular basis of allele-specific efficacy of a blood-stage malaria vaccine: vaccine development implications. J Infect Dis. (2013) 207:511–9. doi: 10.1093/infdis/jis709
33
TheraMADoumboOKCoulibalyDLaurensMBOuattaraAKoneAKet al. A field trial to assess a blood-stage malaria vaccine. N Engl J Med. (2011) 365:1004–13. doi: 10.1056/NEJMoa1008115
34
MoserKADrabekEFDwivediAStuckeEMCrabtreeJDaraAet al. Strains used in whole organism Plasmodium falciparum vaccine trials differ in genome structure, sequence, and immunogenic potential. Genome Med. (2020) 12:6. doi: 10.1186/s13073-019-0708-9
35
LykeKELaurensMBStraussKAdamsMBillingsleyPFJamesEet al. Optimizing intradermal administration of cryopreserved plasmodium falciparum sporozoites in controlled human malaria infection. Am J Trop Med Hyg. (2015) 93:1274–84. doi: 10.4269/ajtmh.15-0341
36
LaurensMBBerryAATravassosMAStraussKAdamsMShresthaBet al. Dose-dependent infectivity of aseptic, purified, cryopreserved plasmodium falciparum 7G8 sporozoites in malaria-naive adults. J Infect Dis. (2019) 220:1962–6. doi: 10.1093/infdis/jiz410
37
AgrawalSMoserKAMortonLCummingsMPPariharADwivediAet al. Association of a novel mutation in the plasmodium falciparum chloroquine resistance transporter with decreased piperaquine sensitivity. J Infect Dis. (2017) 216:468–76. doi: 10.1093/infdis/jix334
38
LegutkiJBJohnstonSA. Immunosignatures can predict vaccine efficacy. Proc Natl Acad Sci U.S.A. (2013) 110:18614–9. doi: 10.1073/pnas.1309390110
39
ImholteGCSauteraudRKorberBBailerRTTurkETShenXet al. A computational framework for the analysis of peptide microarray antibody binding data with application to HIV vaccine profiling. J Immunol Methods. (2013) 395:1–13. doi: 10.1016/j.jim.2013.06.001
40
SmythGK. Limma: linear models for microarray data. In: Bioinformatics and computational biology solutions using R and Bioconductor. Springer, New York, NY (2005). p. 397–420.
41
BaileyTLBodenMBuskeFAFrithMGrantCEClementiLet al. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. (2009) 37:W202–8. doi: 10.1093/nar/gkp335
42
GuptaSStamatoyannopoulosJABaileyTLNobleWS. Quantifying similarity between motifs. Genome Biol. (2007) 8:R24. doi: 10.1186/gb-2007-8-2-r24
43
CostaJGFaccendiniPLSfercoSJLagierCMMarciparIS. Evaluation and comparison of the ability of online available prediction programs to predict true linear B-cell epitopes. Protein Pept Lett. (2013) 20:724–30. doi: 10.2174/0929866511320060011
44
Friedman-KlabanoffDJLaurensMBBerryAATravassosMAAdamsMStraussKAet al. The controlled human malaria infection experience at the university of maryland. Am J Trop Med Hyg. (2019) 100:556–65. doi: 10.4269/ajtmh.18-0476
45
DaviesHMNofalSDMcLaughlinEJOsborneAR. Repetitive sequences in malaria parasite proteins. FEMS Microbiol Rev. (2017) 41:923–40. doi: 10.1093/femsre/fux046
46
ArdeshirFHowardRFViriyakosolSAradOReeseRT. Cross-reactive asparagine-rich determinants shared between several blood-stage antigens of Plasmodium falciparum and the circumsporozoite protein. Mol Biochem Parasitol. (1990) 40:113–28. doi: 10.1016/0166-6851(90)90085-z
47
RaghavanMKalantarKLDuarteETeyssierNTakahashiSKungAFet al. Antibodies to repeat-containing antigens in Plasmodium falciparum are exposure-dependent and short-lived in children in natural malaria infections. Elife. (2023) 12:e81401. doi: 10.7554/eLife.81401
48
HuangZGunarathneSMSLiuWZhouYJiangYLiSet al. PhIP-Seq: methods, applications and challenges. Front Bioinform. (2024) 4:1424202. doi: 10.3389/fbinf.2024.1424202
49
PolleySDConwayDJCavanaghDRMcBrideJSLoweBSWilliamsTNet al. High levels of serum antibodies to merozoite surface protein 2 of Plasmodium falciparum are associated with reduced risk of clinical malaria in coastal Kenya. Vaccine. (2006) 24:4233–46. doi: 10.1016/j.vaccine.2005.06.030
50
al-YamanFGentonBAndersRTaraikaJGinnyMMellorSet al. Assessment of the role of the humoral response to Plasmodium falciparum MSP2 compared to RESA and SPf66 in protecting Papua New Guinean children from clinical malaria. Parasite Immunol. (1995) 17:493–501. doi: 10.1111/j.1365-3024.1995.tb00920.x
51
MetzgerWGOkenuDMCavanaghDRRobinsonJVBojangKAWeissHAet al. Serum IgG3 to the Plasmodium falciparum merozoite surface protein 2 is strongly associated with a reduced prospective risk of malaria. Parasite Immunol. (2003) 25:307–12. doi: 10.1046/j.1365-3024.2003.00636.x
52
TaylorRRAllenSJGreenwoodBMRileyEM. IgG3 antibodies to Plasmodium falciparum merozoite surface protein 2 (MSP2): increasing prevalence with age and association with clinical immunity to malaria. Am J Trop Med Hyg. (1998) 58:406–13. doi: 10.4269/ajtmh.1998.58.406
53
OsierFHFeganGPolleySDMurungiLVerraFTettehKKet al. Breadth and magnitude of antibody responses to multiple Plasmodium falciparum merozoite antigens are associated with protection from clinical malaria. Infect Immun. (2008) 76:2240–8. doi: 10.1128/IAI.01585-07
54
FowkesFJRichardsJSSimpsonJABeesonJG. The relationship between anti-merozoite antibodies and incidence of Plasmodium falciparum malaria: A systematic review and meta-analysis. PLoS Med. (2010) 7:e1000218. doi: 10.1371/journal.pmed.1000218
55
FerreiraMUHartlDL. Plasmodium falciparum: worldwide sequence diversity and evolution of the malaria vaccine candidate merozoite surface protein-2 (MSP-2). Exp Parasitol. (2007) 115:32–40. doi: 10.1016/j.exppara.2006.05.003
56
SeowJMoralesRAMacRaildCAKrishnarjunaBMcGowanSDingjanTet al. Structure and characterisation of a key epitope in the conserved C-terminal domain of the malaria vaccine candidate MSP2. J Mol Biol. (2017) 429:836–46. doi: 10.1016/j.jmb.2017.02.003
57
WangLTPereiraLSFlores-GarciaYO’ConnorJFlynnBJSchonAet al. A potent anti-malarial human monoclonal antibody targets circumsporozoite protein minor repeats and neutralizes sporozoites in the liver. Immunity. (2020) 53:733–44.e8. doi: 10.1016/j.immuni.2020.08.014
58
SpielmannTFergusenDJBeckHP. etramps, a new Plasmodium falciparum gene family coding for developmentally regulated and highly charged membrane proteins located at the parasite-host cell interface. Mol Biol Cell. (2003) 14:1529–44. doi: 10.1091/mbc.e02-04-0240
59
GarciaJCurtidorHObando-MartinezAZVizcainoCPintoMMartinezNLet al. Synthetic peptides from conserved regions of the Plasmodium falciparum early transcribed membrane and ring exported proteins bind specifically to red blood cell proteins. Vaccine. (2009) 27:6877–86. doi: 10.1016/j.vaccine.2009.09.009
60
CromptonPDKayalaMATraoreBKayentaoKOngoibaAWeissGEet al. A prospective analysis of the Ab response to Plasmodium falciparum before and after a malaria season by protein microarray. Proc Natl Acad Sci U.S.A. (2010) 107:6958–63. doi: 10.1073/pnas.1001323107
61
DatooMSDickoATintoHOuedraogoJBHamalubaMOlotuAet al. Safety and efficacy of malaria vaccine candidate R21/Matrix-M in African children: a multicentre, double-blind, randomised, phase 3 trial. Lancet. (2024) 403:533–44. doi: 10.1016/S0140-6736(23)02511-4
62
DobanoCSanzHSorghoHDosooDMpinaMUbillosIet al. Concentration and avidity of antibodies to different circumsporozoite epitopes correlate with RTS,S/AS01E malaria vaccine efficacy. Nat Commun. (2019) 10:2174. doi: 10.1038/s41467-019-10195-z
63
WhiteMTVerityRGriffinJTAsanteKPOwusu-AgyeiSGreenwoodBet al. Immunogenicity of the RTS,S/AS01 malaria vaccine and implications for duration of vaccine efficacy: secondary analysis of data from a phase 3 randomised controlled trial. Lancet Infect Dis. (2015) 15:1450–8. doi: 10.1016/S1473-3099(15)00239-X
64
KesterKECummingsJFOfori-AnyinamOOckenhouseCFKrzychUMorisPet al. Randomized, double-blind, phase 2a trial of falciparum malaria vaccines RTS,S/AS01B and RTS,S/AS02A in malaria-naive adults: safety, efficacy, and immunologic associates of protection. J Infect Dis. (2009) 200:337–46. doi: 10.1086/600120
Summary
Keywords
malaria, Plasmodium falciparum, controlled human malaria infection, peptide microarrays, humoral immunity to malaria, epitope
Citation
Friedman-Klabanoff DJ, Jensen TL, Lyke KE, Laurens MB, Silva JC, Stucke EM, Ouattara A, Ifeonu OO, Hodges T, Moser KA, Gelber CE, Goll JB, Hoffman SL, Patel JJ, Pinapati RS, Tan JC, Deye GA, Takala-Harrison S, Travassos MA and Berry AA (2025) Controlled human malaria infection with NF54 and 7G8 strains elicit differential antibody responses to Plasmodium falciparum peptides. Front. Immunol. 16:1641280. doi: 10.3389/fimmu.2025.1641280
Received
04 June 2025
Accepted
12 August 2025
Published
15 September 2025
Volume
16 - 2025
Edited by
Xinming Tang, Chinese Academy of Agricultural Sciences (CAAS), China
Reviewed by
Kwadwo Asamoah Kusi, University of Ghana, Ghana
Selorme Adukpo, University of Ghana, Ghana
Updates
Copyright
© 2025 Friedman-Klabanoff, Jensen, Lyke, Laurens, Silva, Stucke, Ouattara, Ifeonu, Hodges, Moser, Gelber, Goll, Hoffman, Patel, Pinapati, Tan, Deye, Takala-Harrison, Travassos and Berry.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Andrea A. Berry, aberry@som.umaryland.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.