-
PDF
- Split View
-
Views
-
Cite
Cite
Vanessa C. Nicolete, Sarah Frischmann, Susana Barbosa, Christopher L. King, Marcelo U. Ferreira, Naturally Acquired Binding-Inhibitory Antibodies to Plasmodium vivax Duffy Binding Protein and Clinical Immunity to Malaria in Rural Amazonians, The Journal of Infectious Diseases, Volume 214, Issue 10, 15 November 2016, Pages 1539–1546, https://doi.org/10.1093/infdis/jiw407
- Share Icon Share
Antibodies to the cysteine-rich domain II of Plasmodium vivax Duffy binding protein (PvDBP) can inhibit binding of this parasite ligand to its receptor on red blood cells, the Duffy antigen/receptor for chemokines. These binding-inhibitory antibodies (BIAbs) also inhibit P. vivax invasion of reticulocytes in vitro.
To investigate whether naturally acquired anti-PvDBP antibodies are associated with reduced risk of clinical malaria in a population exposed to low levels of P. vivax transmission, we measured total levels of immunoglobulin G antibodies to 5 PvDBP variants and used a functional in vitro assay to quantify their binding-inhibitory activity in a cohort of 466 rural Amazonians followed up for up to 37 months.
No association between total immunoglobulin G antibody responses to any PvDBP variant and risk of symptomatic, laboratory-confirmed vivax malaria was observed in this cohort. However, a Cox proportional hazards model, adjusted for age, sex, and genotype for the Duffy antigen/receptor for chemokines, showed a >40% decrease in the prospective risk of clinical vivax malaria in subjects with the strongest BIAb responses (upper and middle terciles). High BIAb responses were mostly PvDBP variant transcending and stable over time.
Strong naturally acquired BIAb responses are associated with a reduced risk of clinical P. vivax malaria in rural Amazonians.
Clinical manifestations of Plasmodium vivax malaria are caused by asexual parasites multiplying inside red blood cells (RBCs). Children exposed to intense transmission experience repeated P. vivax malaria attacks during their first years of life, whereas adolescents and adults may harbor low-grade infections but usually remain free of symptoms once parasitized [1]. Naturally acquired immunity does not prevent infection but limits parasite multiplication and reduces the proportion of infections that progress to disease [2]. Clinical immunity also develops under the conditions of substantially lower P. vivax transmission that prevail in the Amazon. Although malaria episodes are experienced by children and adults, an increased proportion of infections remain asymptomatic in subjects with >5–8 years of continuous exposure to infection [3–5].
Several receptor-ligand interactions are required for P. vivax invasion of RBCs. A key step involves the cysteine-rich domain II of the P. vivax Duffy binding protein (PvDBP) and its cognate receptor on RBCs, the Duffy antigen/receptor for chemokines (DARC). PvDBP binding to DARC is needed to form an irreversible junction between invading P. vivax merozoites and their host cells, immature RBCs known as reticulocytes [6]. RBCs lacking DARC are typically refractory to P. vivax infection [7].
Strong naturally acquired antibody responses to PvDBP are associated with reduced prospective risk of high-density parasitemia in Papua New Guinean children exposed to intense P. vivax transmission [8]. Children with high levels of binding-inhibitory antibodies (BIAbs), which are able to inhibit PvDBP binding to DARC in vitro, are protected from subsequent P. vivax infections detected with microscopy [9]. Moreover, human anti-PvDBP antibodies purified from the serum of these children, as well as anti- Duffy binding protein antibodies elicited in immunized rabbits, can inhibit P. vivax invasion of reticulocytes in vitro [10].
To investigate whether naturally acquired anti-PvDBP antibodies reduce the risk of clinical malaria in a population exposed to low-level P. vivax transmission, we measured levels of immunoglobulin (Ig) G antibodies to 5 PvDBP variants in a cohort of rural Amazonians and used a functional assay to quantify their binding-inhibitory activity. We show that high levels of BIAbs, but not of total anti-PvDBP IgG antibodies, persist for several months, are variant transcending, and confer protection from subsequent clinical vivax malaria in this population.
METHODS
Cross-sectional Surveys
Between March 2010 and May 2013, 7 cross-sectional surveys were carried out in farming settlements scattered in the equatorial rain forest of Remansinho area, northwestern Brazil (Supplementary Figure 1). All inhabitants older than 3 months were invited to contribute 5-mL venous blood samples, irrespective of any clinical symptoms, for malaria diagnosis, antibody assays, and Duffy blood group genotyping [11]. Study protocols were approved by the National Human Research Ethics Committee of the Ministry of Health of Brazil (approval No. 551/2010). Written informed consent was obtained from all study participants or their parents or guardians.
Malaria Diagnosis and Surveillance
Conventional microscopy of Giemsa-stained thick smears and quantitative polymerase chain reaction (qPCR) were used to diagnose malaria in cross-sectional surveys [11]. At least 100 fields were examined for malaria parasites by 2 experienced microscopists, under ×1000 magnification, before slides were declared negative. Species-specific primers were used to amplify a 100–base pair fragment of the 18S ribosomal RNA genes of P. vivax and Plasmodium. falciparum by real-time qPCR, with a detection threshold of 2 parasites/microliter of blood [12]. To estimate the incidence of clinical malaria over the study period, trained local health workers carried out weekly house-to-house visits from March 2010 to October 2013. Thick smears were systematically obtained from subjects with fever, chills, headache, or any other malaria-related symptoms and or signs at the time of blood collection or those reporting these symptoms within the past 48 hours. Slide-positive subjects were treated according to the current malaria therapy guidelines of the Ministry of Health of Brazil [13].
DARC Genotyping
TaqMan assays [14] were used to genotype 2 DARC polymorphisms: the T-33C substitution in the RBC-specific GATA1 transcription factor binding motif (rs2814778), which suppresses DARC expression on RBC surface (FY*BES allele), and the G125A polymorphism (rs12075), which defines the FY*B (wild type) and FY*A (mutated) alleles.
Recombinant Antigens
We measured IgG antibodies to the following: (1) 5 variants (Sal I, AH, P, C, and O [8]) of the domain II (residues 194–521) of PvDBP [15]; (2) the C-terminal, 19-kDa region of P. vivax merozoite surface protein (MSP) 1 (PvMSP-119), Belém strain, kindly provided by Anthony Stowers (National Institute of Allergy and Infectious Diseases, National Institutes of Health); (3) block II (nucleotides 1246–2058) of P. vivax MSP-3α (PvMSP-3α) from the Belém strain [16]; (4) the N-terminal region (residues 34–193) of P. vivax MSP-9 (PvMSP-9) from the Belém strain [17], kindly provided by Mary Galinski (Emory University); and (5) the F2 region of P. falciparum erythrocyte-binding antigen (PfEBA) 175, kindly provided by Chetan Chitnis (International Centre for Genetic Engineering and Biotechnology).
Bead-Based Multiplex Assay for Antibodies
Spectrally unique carboxylated xMAP microspheres (Luminex) coupled to recombinant proteins were used to detect antigen-specific total IgG antibodies [18]. Briefly, 1000 coupled beads per well for each antigen were incubated with patients' plasma at a 1:200 dilution in 96-well microplates (MultiScreenBV MSBVN1210; Millipore). R-phycoerythrin–conjugated, anti-human IgG (Jackson Immunoresearch) was used to detect antibody binding. Serially diluted plasma samples from a pool of infected adults from Papua New Guinea were used as positive controls, starting at 1:50 dilution, and plasma samples from 20 North Americans who had not been exposed to malaria served as negative controls. (We did not use samples from Amazonians as negative controls because we would be unable to rule out past exposure to malaria.) We read fluorescence from 75 beads per antigen on Bio-Plex 200 equipment with Bioplex Manger 6.1 software (Bio-Rad). The mean fluorescence intensity (MFI) obtained with negative controls for each antigen, plus 3 standard deviations, was used as the cutoff value to define a positive response. Quantitative results were expressed as reactivity indices, calculated as the sample MFI divided by the cutoff MFI value; an index >1 was considered positive.
Microplate-Based Binding Inhibition Assay
Functional assays were performed in duplicate to test plasma samples for the presence of BIAbs [19]. Plasma samples were preincubated at 1:20 dilution with the recombinant Sal I variant of PvDBP (0.011 µg/mL; kindly provided by Niraj Tolia, Washington University School of Medicine) and transferred to Immulon 2B microplates (ThermoScientific) coated (at 12.5 ng per well) with the N-terminal region of the human DARC protein ligated to the Fc region of human IgG (nDARC-Fc). Bound PvDBP was detected by anti-PvDBP antibodies raised in rabbits (1:5000 dilution), followed by a horseradish peroxidase–conjugated goat, anti-rabbit IgG (Millipore) at 1:5000 dilution. As a positive control we used a human monoclonal antibody against PvDBP, Sal-I variant; plasma samples from unexposed North Americans served as negative controls. Net sample absorbance values (measured at 450 nm) were obtained after subtracting background absorbance readings from wells with neither nDARC-Fc nor plasma. The percentage of binding inhibition was calculated as 1 − (net sample absorbance/net control absorbance) × 100. To assess the variant specificity of BIAbs, we tested serially diluted plasma samples (from 1:20 to 1:640) from selected study participants and nonexposed controls for their ability to inhibit the binding of 3 PvDBP recombinant proteins (variants Sal I, P, and C) to solid-phase nDARC-Fc. The BIAb end-point titer was defined as the reciprocal of the highest plasma dilution at which a blocking-inhibitory activity ≥90% could still be detected.
pvdbp Gene Sequencing
To determine whether PvDBP variants in Remansinho differed from those expressed as recombinant antigens, we sequenced the domain II (nucleotide residues 870–1545) of the pvdbp gene from 46 local P. vivax isolates collected between 2010 and 2011. The 676–base pair target sequence was amplified with PCR [20] and sequenced on an ABI Prism 3100 Genetic Analyzer (Applied Biosystems). Nucleotide sequences were deposited into the GenBank database (accession Nos. KP036999–KP037006).
Statistical Analysis
All analyses were performed with R statistical software (version 3.3.0). Proportions were compared using standard χ2 tests or Mantel-Haenszel χ2 tests for linear trend, while continuous variables were compared using the nonparametric Mann–Whitney U test. Correlations were evaluated using nonparametric Spearman correlation. Statistical significance was defined at the 5% level; 95% confidence intervals (CIs) and interquartile ranges were estimated whenever appropriate.
Survival analysis was used to compare the time to the first clinical vivax malaria episode in subjects with varying levels of specific antibodies. The outcome was clinical vivax malaria, defined as a laboratory-confirmed infection with P. vivax, regardless of parasite density, in a subject with fever, headache, or any other malaria-related symptoms or and signs at enrollment or reporting these symptoms in the past 48 hours. At each cross-sectional survey, study participants who were free of malaria infection (negative qPCR and microscopy) had their antibody levels measured. They were followed up until the next antibody status measurement. The unit of analysis was the time (in days) between the first (baseline) and the next antibody measurement, between the baseline antibody measurement and the date when subjects left the study or between the baseline antibody measurement and the next vivax malaria episode, whatever came first; each study participant contributed up to 7 observation periods.
We used mixed-effects Cox proportional hazards models to compare hazard ratios (HRs) for the time to the first clinical vivax malaria episode across terciles of antibody levels while adjusting for subjects' age and sex and for DARC genotype; the clustering of repeated observations within individuals at different time points was modeled as a random effect [21]. DARC-negative subjects were excluded from survival analysis because they are at negligible risk of vivax malaria [7]; all DARC-positive subjects (regardless of their genotype) were included. A similar approach was used to analyze HRs for the time to the first clinical vivax malaria episode according to levels of anti-PvDBP BIAbs. For this purpose, subjects were classified into terciles for the percentage of binding-inhibitory activity against the Sal I variant of PvDBP.
RESULTS
Naturally Acquired IgG Antibodies to PvDBP

Prevalence of immunoglobulin G antibodies to 5 Plasmodium vivax Duffy binding protein (PvDBP) variants (Sal I, AH, C, P, and O) in rural Amazonians according to age. Study participants had 1–7 consecutive samples tested for antibodies; this figure shows antibody prevalence data at enrollment for 451 participants with known age. Age groups roughly correspond to quintiles, with 87–91 subjects in each group. Extended Mantel-Haenszel χ2 tests for linear trend showed a significant increase in antibody prevalence with increasing age for all PvDBP variants (P < .001).
Sequence Diversity and Antibody Recognition of PvDBP Variants
None of the 8 PvDBP variants characterized in Remansinho was identical to any PvDBP variant expressed as recombinant antigen (Table 1); pairwise differences between local variants and antigens ranged from 1 (AH vs BR7) to 8 amino acid residues (P vs BR1 and BR2). Interestingly, most study subjects recognized PvDBP variants to which they are most likely not exposed. Moreover, levels of IgG antibodies to different PvDBP recombinant proteins at study subjects' enrollment were highly correlated to each other (P < .001 for all; Supplementary Table 3), suggesting that sequence diversity in PvDBP domain II of recombinant antigens had relatively little impact on specific IgG antibody measurements.
Polymorphic Amino Acid Residues in PvDBP Variants Expressed as Recombinant Antigens and Those Characterized in 46 P. vivax Isolates From the Study Site (Remansinho, Brazil)
PvDBP Variant . | Amino Acid Residue (Position According to Sal I Sequence)a . | Proportion in Remansinho, % . | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
308 . | 333 . | 371 . | 375 . | 384 . | 385 . | 386 . | 390 . | 417 . | 424 . | 437 . | 447 . | 503 . | ||
Sal I | R | L | K | N | D | E | K | R | N | L | W | S | I | 0 |
AH | S | … | E | … | G | … | Q | … | K | I | R | … | K | 0 |
C | S | … | … | … | G | … | Q | … | K | I | R | … | K | 0 |
O | S | … | … | … | G | … | … | H | … | I | … | K | … | 0 |
P | S | F | … | D | G | K | N | H | K | I | R | … | K | 0 |
BR1 | … | … | E | … | G | … | … | … | K | I | R | … | … | 34.8 |
BR2 | … | … | … | … | G | … | … | H | … | … | … | … | K | 23.9 |
BR3 | … | … | … | … | G | K | N | H | … | … | R | … | K | 15.2 |
BR4 | S | … | E | … | G | K | N | … | K | I | R | … | K | 10.9 |
BR5 | S | … | … | … | G | K | N | H | … | I | … | … | K | 8.7 |
BR6 | … | … | … | … | G | … | … | … | K | I | R | … | … | 2.2 |
BR7 | S | … | E | … | G | … | … | … | K | I | R | … | K | 2.2 |
BR8 | … | … | … | … | G | K | … | H | … | … | R | … | … | 2.2 |
PvDBP Variant . | Amino Acid Residue (Position According to Sal I Sequence)a . | Proportion in Remansinho, % . | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
308 . | 333 . | 371 . | 375 . | 384 . | 385 . | 386 . | 390 . | 417 . | 424 . | 437 . | 447 . | 503 . | ||
Sal I | R | L | K | N | D | E | K | R | N | L | W | S | I | 0 |
AH | S | … | E | … | G | … | Q | … | K | I | R | … | K | 0 |
C | S | … | … | … | G | … | Q | … | K | I | R | … | K | 0 |
O | S | … | … | … | G | … | … | H | … | I | … | K | … | 0 |
P | S | F | … | D | G | K | N | H | K | I | R | … | K | 0 |
BR1 | … | … | E | … | G | … | … | … | K | I | R | … | … | 34.8 |
BR2 | … | … | … | … | G | … | … | H | … | … | … | … | K | 23.9 |
BR3 | … | … | … | … | G | K | N | H | … | … | R | … | K | 15.2 |
BR4 | S | … | E | … | G | K | N | … | K | I | R | … | K | 10.9 |
BR5 | S | … | … | … | G | K | N | H | … | I | … | … | K | 8.7 |
BR6 | … | … | … | … | G | … | … | … | K | I | R | … | … | 2.2 |
BR7 | S | … | E | … | G | … | … | … | K | I | R | … | K | 2.2 |
BR8 | … | … | … | … | G | K | … | H | … | … | R | … | … | 2.2 |
Abbreviation: PvDBP, Plasmodium vivax Duffy binding protein.
a Ellipses indicate identity with the Sal I sequence. Sal I, AH, C, O, and P are PvDBP haplotypes expressed as recombinant antigens for antibody detection in this study; BR1–BR8 are the PvDBP haplotypes found in the P. vivax population of the study site.
Polymorphic Amino Acid Residues in PvDBP Variants Expressed as Recombinant Antigens and Those Characterized in 46 P. vivax Isolates From the Study Site (Remansinho, Brazil)
PvDBP Variant . | Amino Acid Residue (Position According to Sal I Sequence)a . | Proportion in Remansinho, % . | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
308 . | 333 . | 371 . | 375 . | 384 . | 385 . | 386 . | 390 . | 417 . | 424 . | 437 . | 447 . | 503 . | ||
Sal I | R | L | K | N | D | E | K | R | N | L | W | S | I | 0 |
AH | S | … | E | … | G | … | Q | … | K | I | R | … | K | 0 |
C | S | … | … | … | G | … | Q | … | K | I | R | … | K | 0 |
O | S | … | … | … | G | … | … | H | … | I | … | K | … | 0 |
P | S | F | … | D | G | K | N | H | K | I | R | … | K | 0 |
BR1 | … | … | E | … | G | … | … | … | K | I | R | … | … | 34.8 |
BR2 | … | … | … | … | G | … | … | H | … | … | … | … | K | 23.9 |
BR3 | … | … | … | … | G | K | N | H | … | … | R | … | K | 15.2 |
BR4 | S | … | E | … | G | K | N | … | K | I | R | … | K | 10.9 |
BR5 | S | … | … | … | G | K | N | H | … | I | … | … | K | 8.7 |
BR6 | … | … | … | … | G | … | … | … | K | I | R | … | … | 2.2 |
BR7 | S | … | E | … | G | … | … | … | K | I | R | … | K | 2.2 |
BR8 | … | … | … | … | G | K | … | H | … | … | R | … | … | 2.2 |
PvDBP Variant . | Amino Acid Residue (Position According to Sal I Sequence)a . | Proportion in Remansinho, % . | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
308 . | 333 . | 371 . | 375 . | 384 . | 385 . | 386 . | 390 . | 417 . | 424 . | 437 . | 447 . | 503 . | ||
Sal I | R | L | K | N | D | E | K | R | N | L | W | S | I | 0 |
AH | S | … | E | … | G | … | Q | … | K | I | R | … | K | 0 |
C | S | … | … | … | G | … | Q | … | K | I | R | … | K | 0 |
O | S | … | … | … | G | … | … | H | … | I | … | K | … | 0 |
P | S | F | … | D | G | K | N | H | K | I | R | … | K | 0 |
BR1 | … | … | E | … | G | … | … | … | K | I | R | … | … | 34.8 |
BR2 | … | … | … | … | G | … | … | H | … | … | … | … | K | 23.9 |
BR3 | … | … | … | … | G | K | N | H | … | … | R | … | K | 15.2 |
BR4 | S | … | E | … | G | K | N | … | K | I | R | … | K | 10.9 |
BR5 | S | … | … | … | G | K | N | H | … | I | … | … | K | 8.7 |
BR6 | … | … | … | … | G | … | … | … | K | I | R | … | … | 2.2 |
BR7 | S | … | E | … | G | … | … | … | K | I | R | … | K | 2.2 |
BR8 | … | … | … | … | G | K | … | H | … | … | R | … | … | 2.2 |
Abbreviation: PvDBP, Plasmodium vivax Duffy binding protein.
a Ellipses indicate identity with the Sal I sequence. Sal I, AH, C, O, and P are PvDBP haplotypes expressed as recombinant antigens for antibody detection in this study; BR1–BR8 are the PvDBP haplotypes found in the P. vivax population of the study site.
Antibody Responses to Other Malaria Antigens
IgG response rates were high for PvMSP-119 (67.5%) but substantially lower for PvMSP-3α and PvMSP-9 (Supplementary Table 1). Only 32.6% of the samples had IgG antibodies to PfEBA-175, a member of the Duffy binding–like erythrocyte-binding protein (erythrocyte-binding antigen) family of P. falciparum with no known orthologue with significant sequence similarity in P. vivax. Levels of IgG antibodies to PvMSP-119, PvMSP-9, and PfEBA-175, but not to PvMSP-3α, correlated weakly but positively with the subjects' age (Supplementary Table 2).
Malaria Prevalence and Antibody Levels Over Time

Temporal variation in Plasmodium vivax prevalence and levels of immunoglobulin (Ig) G antibodies to P. vivax Duffy binding protein (PvDBP) Sal I and the C-terminal, 19-kDa region of P. vivax merozoite surface protein 1 (PvMSP-119) in rural Amazonians. Top panel, Prevalence of P. vivax parasitemia detected with microscopy and/or quantitative polymerase chain reaction (qPCR) in 7 consecutive cross-sectional surveys in March–May 2010 (baseline), May–July 2010, March–April 2011, October–November 2011, April–May 2012, October–November 2012, and April–May 2013 (solid line); dashed lines represent 95% confidence intervals lines. Middle and bottom panels, Median reactivity indices of IgG antibodies to PvDBP Sal I (middle) and PvMSP-119 (bottom), detected with a multiplex assay during the same cross-sectional surveys (solid lines); dashed lines represent interquartile ranges. Data for other PvDBP variants and merozoite surface protein antigens are shown in Supplementary Figure 2.
Antibody Levels and Prospective Risk of Clinical Vivax Malaria
An average of 8.4 clinical vivax malaria episodes per 100 person-months at risk were recorded in the study population over the study period, with a peak of 39.3 cases per 100 person-months in October 2011 (Supplementary Figure 4). A total of 529 clinical vivax malaria episodes were recorded over the whole study period. Survival analysis showed no delay in time to the first clinical vivax malaria episode in the upper and middle terciles of IgG response to any PvDBP variant or MSP, compared with the lower tercile of antibody response (Supplementary Figure 5). We compared HRs for the time to the first clinical vivax malaria episode across terciles of specific antibody levels, using Cox proportional hazards models adjusted for age, sex, and DARC genotype, and found no significant association between levels of specific IgG antibodies and prospective risk of clinical vivax malaria (Table 2).
Association Between Levels of Naturally Acquired Antibodies to PvDBP and PvMSP Variants and Prospective Risk of Clinical Vivax Malaria in Rural Amazonians
Antigen . | Variant . | Middle vs Lower Tercile of Reactivity Indicesa . | Upper vs Lower Tercile of Reactivity Indicesa . | ||
---|---|---|---|---|---|
HR (95% CI)b . | P Value . | HR (95% CI)b . | P Value . | ||
PvDBP | Sal I | 1.05 (.73–1.51) | .78 | 1.14 (.78–1.67) | .49 |
AH | 0.88 (.62–1.25) | .47 | 0.88 (.60–1.27) | .49 | |
C | 0.91 (.64–1.29) | .60 | 0.97 (.67–1.40) | .89 | |
O | 0.83 (.59–1.16) | .28 | 0.89 (.61–1.26) | .48 | |
P | 0.81 (.56–1.15) | .25 | 0.93 (.64–1.34) | .69 | |
PvMSP-119 | Belém | 0.82 (.58–1.16) | .26 | 1.20 (.85–1.68) | .31 |
PvMSP-3α | Belém | 1.11 (.76–1.61) | .60 | 1.21 (.83–1.77) | .33 |
PvMSP-9 | Belém | 0.86 (.60–1.24) | .42 | 1.01 (.72–1.43) | .93 |
Antigen . | Variant . | Middle vs Lower Tercile of Reactivity Indicesa . | Upper vs Lower Tercile of Reactivity Indicesa . | ||
---|---|---|---|---|---|
HR (95% CI)b . | P Value . | HR (95% CI)b . | P Value . | ||
PvDBP | Sal I | 1.05 (.73–1.51) | .78 | 1.14 (.78–1.67) | .49 |
AH | 0.88 (.62–1.25) | .47 | 0.88 (.60–1.27) | .49 | |
C | 0.91 (.64–1.29) | .60 | 0.97 (.67–1.40) | .89 | |
O | 0.83 (.59–1.16) | .28 | 0.89 (.61–1.26) | .48 | |
P | 0.81 (.56–1.15) | .25 | 0.93 (.64–1.34) | .69 | |
PvMSP-119 | Belém | 0.82 (.58–1.16) | .26 | 1.20 (.85–1.68) | .31 |
PvMSP-3α | Belém | 1.11 (.76–1.61) | .60 | 1.21 (.83–1.77) | .33 |
PvMSP-9 | Belém | 0.86 (.60–1.24) | .42 | 1.01 (.72–1.43) | .93 |
Abbreviations: CI, confidence interval; HR, hazard ratio; PvDBP, Plasmodium vivax Duffy binding protein; P. vivax PvMSP, merozoite surface protein.
a Antibody levels were stratified into terciles of reactivity indices for analysis.
b HRs were obtained with Cox proportional hazards models adjusted for age, sex, and genotype for the Duffy antigen/receptor for chemokines.
Association Between Levels of Naturally Acquired Antibodies to PvDBP and PvMSP Variants and Prospective Risk of Clinical Vivax Malaria in Rural Amazonians
Antigen . | Variant . | Middle vs Lower Tercile of Reactivity Indicesa . | Upper vs Lower Tercile of Reactivity Indicesa . | ||
---|---|---|---|---|---|
HR (95% CI)b . | P Value . | HR (95% CI)b . | P Value . | ||
PvDBP | Sal I | 1.05 (.73–1.51) | .78 | 1.14 (.78–1.67) | .49 |
AH | 0.88 (.62–1.25) | .47 | 0.88 (.60–1.27) | .49 | |
C | 0.91 (.64–1.29) | .60 | 0.97 (.67–1.40) | .89 | |
O | 0.83 (.59–1.16) | .28 | 0.89 (.61–1.26) | .48 | |
P | 0.81 (.56–1.15) | .25 | 0.93 (.64–1.34) | .69 | |
PvMSP-119 | Belém | 0.82 (.58–1.16) | .26 | 1.20 (.85–1.68) | .31 |
PvMSP-3α | Belém | 1.11 (.76–1.61) | .60 | 1.21 (.83–1.77) | .33 |
PvMSP-9 | Belém | 0.86 (.60–1.24) | .42 | 1.01 (.72–1.43) | .93 |
Antigen . | Variant . | Middle vs Lower Tercile of Reactivity Indicesa . | Upper vs Lower Tercile of Reactivity Indicesa . | ||
---|---|---|---|---|---|
HR (95% CI)b . | P Value . | HR (95% CI)b . | P Value . | ||
PvDBP | Sal I | 1.05 (.73–1.51) | .78 | 1.14 (.78–1.67) | .49 |
AH | 0.88 (.62–1.25) | .47 | 0.88 (.60–1.27) | .49 | |
C | 0.91 (.64–1.29) | .60 | 0.97 (.67–1.40) | .89 | |
O | 0.83 (.59–1.16) | .28 | 0.89 (.61–1.26) | .48 | |
P | 0.81 (.56–1.15) | .25 | 0.93 (.64–1.34) | .69 | |
PvMSP-119 | Belém | 0.82 (.58–1.16) | .26 | 1.20 (.85–1.68) | .31 |
PvMSP-3α | Belém | 1.11 (.76–1.61) | .60 | 1.21 (.83–1.77) | .33 |
PvMSP-9 | Belém | 0.86 (.60–1.24) | .42 | 1.01 (.72–1.43) | .93 |
Abbreviations: CI, confidence interval; HR, hazard ratio; PvDBP, Plasmodium vivax Duffy binding protein; P. vivax PvMSP, merozoite surface protein.
a Antibody levels were stratified into terciles of reactivity indices for analysis.
b HRs were obtained with Cox proportional hazards models adjusted for age, sex, and genotype for the Duffy antigen/receptor for chemokines.
Naturally Acquired Anti-PvDBP BIAbs

Naturally acquired binding-inhibitory antibodies (BIAbs) to Plasmodium vivax Duffy binding protein (PvDBP) in rural Amazonians. BIAb activity was expressed as the percentage of inhibition of PvDBP (Sal I variant) binding to solid-phase Duffy antigen/receptor for chemokines (DARC) in vitro. A, Frequency distribution of the percentage of binding inhibition of anti-PvDBP antibodies at enrollment in 263 rural Amazonians; quite similar results were obtained with the complete data set (n = 572 samples). B, Kaplan–Meier curve showing the proportion of rural Amazonians who remained free of microscopically confirmed clinical vivax malaria during follow-up, according to their levels of anti-PvDBP BIAbs. Study participants (n = 229) who were DARC positive and free of P. vivax infection at the time of antibody measurement were grouped into terciles of BIAb activity (upper, ≥50%; middle, 6%–49%; lower, <6%) and followed up for up to 37 months. The unit of analysis was the period of observation (in days) between the baseline and the next BIAb measurement, the first vivax malaria episode, or the date when subjects left the study, whichever came first. Each subject contributed between 1 and 6 periods of observation, with 456 BIAb measurements (n = 152 in each BIAb tercile) and 126 events (33, 39, and 54 in the upper, middle, and lower BIAb terciles, respectively) analyzed.
Variant Specificity of BIAbs

Pairwise comparisons of binding-inhibitory antibody (BIAb) end-point titers to Plasmodium vivax Duffy binding protein variants in rural Amazonians. Top panel, Sal I variant (columns) versus P variant (rows). Middle panel, Sal I variant (columns) versus C variant (rows). Bottom panel, P variant (columns) versus C variant (rows). Black boxes indicate comparisons with a ≤2-fold difference in end-point titer between variants (“variant-transcending” BIAbs); gray boxes, comparisons with a 4-fold difference in end-point titer (“moderately variant-specific” BIAbs); and white boxes, comparisons with a >4-fold difference in end-point titer (“variant-specific” BIAbs). Numbers of samples with each end-point titer are indicated within boxes; results for 164 plasma samples from 77 subjects (1–7 samples per study participant) are analyzed.
BIAbs and Prospective Risk of Clinical Vivax Malaria
We next tested whether strong BIAb responses protected from clinical malaria. To this end, we grouped study participants into terciles of BIAb response. Survival analysis showed a delay in time to the first clinical vivax malaria episode among subjects in the upper and middle terciles of BIAb response, compared with their counterparts in the lower tercile (Figure 3B). A Cox proportional hazards model adjusted for age, sex, and DARC genotype showed a 41%–42% reduction in the prospective risk of clinical malaria among subjects in the upper and middle terciles of BIAb response at baseline, relative to those in the lower tercile (Table 3). Age (in years) was another independent predictor of reduced risk of clinical vivax malaria in the fully adjusted Cox model (HR, 0.98; 95% CI, .97–1.00; P = .005), suggesting that age-related factors other than BIAbs are also associated with protection. Female sex (HR, 1.45; 95% CI, 1.01–2.08; P = .04) was independently associated with an increased risk of malaria, after adjustment for BIAb response, but we have no clear-cut explanation for this finding. Interestingly, the Cox model revealed no significant interaction between DARC genotype and BIAb response (P values between .27 and .98) in their association with the outcome, despite the recent finding of greater in vitro binding-inhibitory activity of anti-PvDBP antibodies with Fya RBCs, compared with Fyb RBCs [23].
Association Between Binding-Inhibitory Activity of Naturally Acquired Antibodies to PvDBP and Prospective Risk of Clinical Vivax Malaria in Rural Amazonians
Cox Model . | Middle vs Lower Tercile of BIAb Responsesa . | Upper vs Lower Tercile of BIAb Responsesa . | ||
---|---|---|---|---|
HR (95% CI)b . | P Value . | HR (95% CI)b . | P Value . | |
Unadjusted | 0.59 (.39–.90) | .01 | 0.51 (.33–.79) | .003 |
Adjusted for age | 0.62 (.41–.94) | .02 | 0.59 (.38–.93) | .02 |
Adjusted for sex | 0.57 (.43–.97) | .04 | 0.48 (.37–.86) | .008 |
Adjusted for DARC genotype | 0.56 (.37–.86) | .008 | 0.54 (.34–.84) | .007 |
Adjusted for all variables above | 0.58 (.38–.89) | .01 | 0.59 (.37–.93) | .02 |
Cox Model . | Middle vs Lower Tercile of BIAb Responsesa . | Upper vs Lower Tercile of BIAb Responsesa . | ||
---|---|---|---|---|
HR (95% CI)b . | P Value . | HR (95% CI)b . | P Value . | |
Unadjusted | 0.59 (.39–.90) | .01 | 0.51 (.33–.79) | .003 |
Adjusted for age | 0.62 (.41–.94) | .02 | 0.59 (.38–.93) | .02 |
Adjusted for sex | 0.57 (.43–.97) | .04 | 0.48 (.37–.86) | .008 |
Adjusted for DARC genotype | 0.56 (.37–.86) | .008 | 0.54 (.34–.84) | .007 |
Adjusted for all variables above | 0.58 (.38–.89) | .01 | 0.59 (.37–.93) | .02 |
Abbreviations: BIAb, binding-inhibitory antibody; CI, confidence interval; DARC, Duffy antigen/receptor for chemokines; HR, hazard ratio; PvDBP, Plasmodium vivax Duffy binding protein.
a Antibody responses were stratified into terciles of the percentage of binding inhibition: upper (≥50%), middle (6%–49%), and lower (<6%).
b HRs were obtained with Cox proportional hazards models.
Association Between Binding-Inhibitory Activity of Naturally Acquired Antibodies to PvDBP and Prospective Risk of Clinical Vivax Malaria in Rural Amazonians
Cox Model . | Middle vs Lower Tercile of BIAb Responsesa . | Upper vs Lower Tercile of BIAb Responsesa . | ||
---|---|---|---|---|
HR (95% CI)b . | P Value . | HR (95% CI)b . | P Value . | |
Unadjusted | 0.59 (.39–.90) | .01 | 0.51 (.33–.79) | .003 |
Adjusted for age | 0.62 (.41–.94) | .02 | 0.59 (.38–.93) | .02 |
Adjusted for sex | 0.57 (.43–.97) | .04 | 0.48 (.37–.86) | .008 |
Adjusted for DARC genotype | 0.56 (.37–.86) | .008 | 0.54 (.34–.84) | .007 |
Adjusted for all variables above | 0.58 (.38–.89) | .01 | 0.59 (.37–.93) | .02 |
Cox Model . | Middle vs Lower Tercile of BIAb Responsesa . | Upper vs Lower Tercile of BIAb Responsesa . | ||
---|---|---|---|---|
HR (95% CI)b . | P Value . | HR (95% CI)b . | P Value . | |
Unadjusted | 0.59 (.39–.90) | .01 | 0.51 (.33–.79) | .003 |
Adjusted for age | 0.62 (.41–.94) | .02 | 0.59 (.38–.93) | .02 |
Adjusted for sex | 0.57 (.43–.97) | .04 | 0.48 (.37–.86) | .008 |
Adjusted for DARC genotype | 0.56 (.37–.86) | .008 | 0.54 (.34–.84) | .007 |
Adjusted for all variables above | 0.58 (.38–.89) | .01 | 0.59 (.37–.93) | .02 |
Abbreviations: BIAb, binding-inhibitory antibody; CI, confidence interval; DARC, Duffy antigen/receptor for chemokines; HR, hazard ratio; PvDBP, Plasmodium vivax Duffy binding protein.
a Antibody responses were stratified into terciles of the percentage of binding inhibition: upper (≥50%), middle (6%–49%), and lower (<6%).
b HRs were obtained with Cox proportional hazards models.
DISCUSSION
This cohort study provides the first evidence that naturally acquired anti-PvDBP BIAb responses are associated with reduced risk of clinical vivax malaria. Of note, high BIAb responses developed under conditions of low malaria endemicity that are typical of rural Amazonian communities and, once acquired, were predominantly variant transcending and long lasting. Moreover, the protective effect of anti-PvDBP BIAbs was not modulated by hosts' DARC genotype. These results further support PvDBP as the most promising vaccine target for P. vivax blood stages.
Only 8.7% of children aged 5–14 years and exposed to very intense malaria transmission in Papua New Guinea have been found to acquire BIAbs with >90% inhibitory activity [9]. Interestingly, these children with high levels of BIAbs at baseline had a reduced risk of P. vivax malaria diagnosed by means of conventional microscopy over the next 25 weeks, compared with those with <50% inhibitory activity, although no protection was observed against P. vivax infections diagnosed with a more sensitive molecular method. Moreover, these children, once infected, had comparatively reduced P. vivax densities [9]. These results indicate that high levels of BIAbs in hyperendemic settings are associated with reduced parasite growth, often resulting in low-density P. vivax infections that are missed by conventional microscopy, rather than fully sterilizing immunity. This previous study was unable to examine the association of high levels of BIAbs with protection against clinical immunity because by age 5 years almost all the children in this population were immune to P. vivax illness.
Our data underline the need for functional assays to detect anti-PvDBP antibodies mediating protective immunity, rather than antibodies that are markers of increased exposure to infection [24]. We found no association between high total IgG responses to PvDBP and clinical immunity to malaria in rural Amazonians, suggesting that our conventional serology measures both protective and nonprotective antibodies. Whether BIAbs target mostly polymorphic or conserved epitopes of PvDBP domain II remains to be determined, but our titration experiments with different PvDBP variants suggest that most strongly inhibitory antibodies are variant transcending. To better analyze the fine specificity of naturally acquired BIAbs, we are currently generating a panel of human monoclonal antibodies, from rural Amazonians whose anti-PvDBP antibodies have high and long-lasting binding-inhibitory activity, for use in epitope mapping and further functional characterization of these potentially protective responses.
Given our inability to distinguish protective from nonprotective antibodies by conventional serology, the absence of significant associations between IgG responses to PvMSP-119, PvMSP-3α, and PvMSP-9 and protection from clinical vivax malaria is not entirely surprising. A single cohort study of young Papua New Guinean children aged 1–3 years has shown an association between high levels of antibodies to both PvMSP-3α and PvMSP-9 and reduced incidence of clinical vivax malaria over 16 months of follow-up [25]. Of note, no association between antibody responses to PvMSP-119 and reduced risk of P. vivax infection [8, 26, 27] or high-density parasitemia [8, 27] has been found in cohort studies so far. In contrast, the association between anti-PfMSP-119 antibodies and protection from P. falciparum malaria has been well demonstrated in Africa and Papua New Guinea [28].
We conclude that high titers of long-lasting BIAbs develop and can confer broadly specific clinical immunity to P. vivax under conditions of low malaria transmission, with clear implications for the development of PvDBP-based vaccines. As shown elsewhere [29], repeated exposure to the parasite is not necessarily required to maintain strong BIAb responses, once they develop. Moreover, BIAbs seem to protect against several PvDBP variants, putatively because they target conserved, rather than variant-specific epitopes [24]. Therefore, a PvDBP vaccine that elicits antibody responses with these properties is very likely to be successful in preventing P. vivax infection and clinical disease.
Notes
Acknowledgments. We thank all inhabitants in Remansinho for their enthusiastic participation in this study; Amanda B. Gozze, Pablo S. Fontoura, Kézia Katiani G. Scopel, Nathália F. Lima, Carlos E. Cavasini, Raquel M. Gonçalves, Maria José Menezes, Rosely S. Malafronte, Camilla L. Batista, Ariel M. Silber, Cristiana F. Alves de Brito, Mônica da Silva-Nunes, Carla Roberta O. Carvalho, and Mauro R. Tucci for the help in field work and clinical care of patients; Cleide F. Nunes and Eusueli Arraes da Silva for microscopic diagnosis of malaria; Márcio C. Santana, Andrecresa N. Duarte, and Francisco Naildo C. Leitão for overall logistic support; and Lenore Carias for laboratory support.
Disclaimer. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Financial support. This work was supported by the National Institute of Allergy and Infectious Diseases, National Institutes of Health (International Centers of Excellence in Malaria Research program; research grants U19 AI089681 to Joseph M. Vinetz, U19 AI089686 to James W. Kazura, and AI064478 to C. L. K.); the Veterans Administration Research Service (C. L. K.), the Fundação de Amparo à Pesquisa do Estado de São Paulo, Brazil (grants 2009/52729-9 to M. U. F. and 2013/23770-6 to S. B.), and the Conselho Nacional de Desenvolvimento Científico e Tecnológico of Brazil (scholarship to V. C. N. and senior researcher scholarship to M. U. F.).
Potential conflicts of interest. All authors: No reported conflicts. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed.
References
Author notes
C. L. K. and M. U. F. contributed equally to this work.