Geospatial joint modeling of vector and parasite serology to microstratify malaria transmission.

The World Health Organization identifies a strong surveillance system for malaria and its mosquito vector as an essential pillar of the malaria elimination agenda. Anopheles salivary antibodies are emerging biomarkers of exposure to mosquito bites that potentially overcome sensitivity and logistical constraints of traditional entomological surveys. Using samples collected by a village health volunteer network in 104 villages in Southeast Myanmar during routine surveillance, the present study employs a Bayesian geostatistical modeling framework, incorporating climatic and environmental variables together with Anopheles salivary antigen serology, to generate spatially continuous predictive maps of Anopheles biting exposure. Our maps quantify fine-scale spatial and temporal heterogeneity in Anopheles salivary antibody seroprevalence (ranging from 9 to 99%) that serves as a proxy of exposure to Anopheles bites and advances current static maps of only Anopheles occurrence. We also developed an innovative framework to perform surveillance of malaria transmission. By incorporating antibodies against the vector and the transmissible form of malaria (sporozoite) in a joint Bayesian geostatistical model, we predict several foci of ongoing transmission. In our study, we demonstrate that antibodies specific for Anopheles salivary and sporozoite antigens are a logistically feasible metric with which to quantify and characterize heterogeneity in exposure to vector bites and malaria transmission. These approaches could readily be scaled up into existing village health volunteer surveillance networks to identify foci of residual malaria transmission, which could be targeted with supplementary interventions to accelerate progress toward elimination.

Sulfonylpiperazine compounds prevent Plasmodium falciparum invasion of red blood cells through interference with actin-1/profilin dynamics.

With emerging resistance to frontline treatments, it is vital that new antimalarial drugs are identified to target Plasmodium falciparum. We have recently described a compound, MMV020291, as a specific inhibitor of red blood cell (RBC) invasion, and have generated analogues with improved potency. Here, we generated resistance to MMV020291 and performed whole genome sequencing of 3 MMV020291-resistant populations. This revealed 3 nonsynonymous single nucleotide polymorphisms in 2 genes; 2 in profilin (N154Y, K124N) and a third one in actin-1 (M356L). Using CRISPR-Cas9, we engineered these mutations into wild-type parasites, which rendered them resistant to MMV020291. We demonstrate that MMV020291 reduces actin polymerisation that is required by the merozoite stage parasites to invade RBCs. Additionally, the series inhibits the actin-1-dependent process of apicoplast segregation, leading to a delayed death phenotype. In vitro cosedimentation experiments using recombinant P. falciparum proteins indicate that potent MMV020291 analogues disrupt the formation of filamentous actin in the presence of profilin. Altogether, this study identifies the first compound series interfering with the actin-1/profilin interaction in P. falciparum and paves the way for future antimalarial development against the highly dynamic process of actin polymerisation.

Plasmodium GPI-anchored micronemal antigen is essential for parasite transmission through the mosquito host.

Plasmodium parasites, the eukaryotic pathogens that cause malaria, feature three distinct invasive forms tailored to the host environment they must navigate and invade for life cycle progression. One conserved feature of these invasive forms is the micronemes, apically oriented secretory organelles involved in egress, motility, adhesion, and invasion. Here we investigate the role of GPI-anchored micronemal antigen (GAMA), which shows a micronemal localization in all zoite forms of the rodent-infecting species Plasmodium berghei. ∆GAMA parasites are severely defective for invasion of the mosquito midgut. Once formed, oocysts develop normally, however, sporozoites are unable to egress and exhibit defective motility. Epitope-tagging of GAMA revealed tight temporal expression late during sporogony and showed that GAMA is shed during sporozoite gliding motility in a similar manner to circumsporozoite protein. Complementation of P. berghei knockout parasites with full-length P. falciparum GAMA partially restored infectivity to mosquitoes, indicating conservation of function across Plasmodium species. A suite of parasites with GAMA expressed under the promoters of CTRP, CAP380, and TRAP, further confirmed the involvement of GAMA in midgut infection, motility, and vertebrate infection. These data show GAMA's involvement in sporozoite motility, egress, and invasion, implicating GAMA as a regulator of microneme function.

Population-level genome-wide STR discovery and validation for population structure and genetic diversity assessment of Plasmodium species.

Short tandem repeats (STRs) are highly informative genetic markers that have been used extensively in population genetics analysis. They are an important source of genetic diversity and can also have functional impact. Despite the availability of bioinformatic methods that permit large-scale genome-wide genotyping of STRs from whole genome sequencing data, they have not previously been applied to sequencing data from large collections of malaria parasite field samples. Here, we have genotyped STRs using HipSTR in more than 3,000 Plasmodium falciparum and 174 Plasmodium vivax published whole-genome sequence data from samples collected across the globe. High levels of noise and variability in the resultant callset necessitated the development of a novel method for quality control of STR genotype calls. A set of high-quality STR loci (6,768 from P. falciparum and 3,496 from P. vivax) were used to study Plasmodium genetic diversity, population structures and genomic signatures of selection and these were compared to genome-wide single nucleotide polymorphism (SNP) genotyping data. In addition, the genome-wide information about genetic variation and other characteristics of STRs in P. falciparum and P. vivax have been available in an interactive web-based R Shiny application PlasmoSTR (https://github.com/bahlolab/PlasmoSTR).

Whole-genome analysis of Malawian Plasmodium falciparum isolates identifies possible targets of allele-specific immunity to clinical malaria.

Individuals acquire immunity to clinical malaria after repeated Plasmodium falciparum infections. Immunity to disease is thought to reflect the acquisition of a repertoire of responses to multiple alleles in diverse parasite antigens. In previous studies, we identified polymorphic sites within individual antigens that are associated with parasite immune evasion by examining antigen allele dynamics in individuals followed longitudinally. Here we expand this approach by analyzing genome-wide polymorphisms using whole genome sequence data from 140 parasite isolates representing malaria cases from a longitudinal study in Malawi and identify 25 genes that encode possible targets of naturally acquired immunity that should be validated immunologically and further characterized for their potential as vaccine candidates.

VIVID: A Web Application for Variant Interpretation and Visualization in Multi-dimensional Analyses.

Large-scale comparative genomics- and population genetic studies generate enormous amounts of polymorphism data in the form of DNA variants. Ultimately, the goal of many of these studies is to associate genetic variants to phenotypes or fitness. We introduce VIVID, an interactive, user-friendly web application that integrates a wide range of approaches for encoding genotypic to phenotypic information in any organism or disease, from an individual or population, in three-dimensional (3D) space. It allows mutation mapping and annotation, calculation of interactions and conservation scores, prediction of harmful effects, analysis of diversity and selection, and 3D visualization of genotypic information encoded in Variant Call Format on AlphaFold2 protein models. VIVID enables the rapid assessment of genes of interest in the study of adaptive evolution and the genetic load, and it helps prioritizing targets for experimental validation. We demonstrate the utility of VIVID by exploring the evolutionary genetics of the parasitic protist Plasmodium falciparum, revealing geographic variation in the signature of balancing selection in potential targets of functional antibodies.

Global diversity and balancing selection of 23 leading Plasmodium falciparum candidate vaccine antigens.

Investigation of the diversity of malaria parasite antigens can help prioritize and validate them as vaccine candidates and identify the most common variants for inclusion in vaccine formulations. Studies of vaccine candidates of the most virulent human malaria parasite, Plasmodium falciparum, have focused on a handful of well-known antigens, while several others have never been studied. Here we examine the global diversity and population structure of leading vaccine candidate antigens of P. falciparum using the MalariaGEN Pf3K (version 5.1) resource, comprising more than 2600 genomes from 15 malaria endemic countries. A stringent variant calling pipeline was used to extract high quality antigen gene 'haplotypes' from the global dataset and a new R-package named VaxPack was used to streamline population genetic analyses. In addition, a newly developed algorithm that enables spatial averaging of selection pressure on 3D protein structures was applied to the dataset. We analysed the genes encoding 23 leading and novel candidate malaria vaccine antigens including csp, trap, eba175, ama1, rh5, and CelTOS. Our analysis shows that current malaria vaccine formulations are based on rare haplotypes and thus may have limited efficacy against natural parasite populations. High levels of diversity with evidence of balancing selection was detected for most of the erythrocytic and pre-erythrocytic antigens. Measures of natural selection were then mapped to 3D protein structures to predict targets of functional antibodies. For some antigens, geographical variation in the intensity and distribution of these signals on the 3D structure suggests adaptation to different human host or mosquito vector populations. This study provides an essential framework for the diversity of P. falciparum antigens to be considered in the design of the next generation of malaria vaccines.

Emergence of artemisinin-resistant Plasmodium falciparum with kelch13 C580Y mutations on the island of New Guinea.

The rapid and aggressive spread of artemisinin-resistant Plasmodium falciparum carrying the C580Y mutation in the kelch13 gene is a growing threat to malaria elimination in Southeast Asia, but there is no evidence of their spread to other regions. We conducted cross-sectional surveys in 2016 and 2017 at two clinics in Wewak, Papua New Guinea (PNG) where we identified three infections caused by C580Y mutants among 239 genotyped clinical samples. One of these mutants exhibited the highest survival rate (6.8%) among all parasites surveyed in ring-stage survival assays (RSA) for artemisinin. Analyses of kelch13 flanking regions, and comparisons of deep sequencing data from 389 clinical samples from PNG, Indonesian Papua and Western Cambodia, suggested an independent origin of the Wewak C580Y mutation, showing that the mutants possess several distinctive genetic features. Identity by descent (IBD) showed that multiple portions of the mutants' genomes share a common origin with parasites found in Indonesian Papua, comprising several mutations within genes previously associated with drug resistance, such as mdr1, ferredoxin, atg18 and pnp. These findings suggest that a P. falciparum lineage circulating on the island of New Guinea has gradually acquired a complex ensemble of variants, including kelch13 C580Y, which have affected the parasites' drug sensitivity. This worrying development reinforces the need for increased surveillance of the evolving parasite populations on the island, to contain the spread of resistance.

Community-based molecular and serological surveillance of subclinical malaria in Myanmar.

BACKGROUND: In the Greater Mekong Subregion (GMS), current malaria surveillance strategies rely on a network of village health volunteers (VHVs) reporting the results of rapid diagnostic tests (RDTs), known to miss many asymptomatic infections. Integration of more sensitive diagnostic molecular and serological measures into the VHV network may improve surveillance of residual malaria transmission in hard-to-reach areas in the region and inform targeted interventions and elimination responses. However, data on residual malaria transmission that would be captured by these measures in the VHV-led testing and treatment surveillance network in the GMS is unknown. METHODS: A total of 114 VHVs were trained to collect dried blood spots from villagers undergoing routine RDTs as part of VHV-led active and passive case detection from April 2015 to June 2016. Samples were subjected to molecular testing (quantitative polymerase chain reaction [qPCR]) to determine Plasmodium falciparum and P. vivax infection and serological testing (against P. falciparum and P. vivax antigens) to determine exposure to P. falciparum and P. vivax. RESULTS: Over 15 months, 114 VHVs performed 32,194 RDTs and collected samples for molecular (n = 13,157) and serological (n = 14,128) testing. The prevalence of molecular-detectable P. falciparum and P. vivax infection was 3.2% compared to the 0.16% prevalence of Plasmodium spp. by RDT, highlighting the large burden of infections undetected by standard surveillance. Peaks in anti-P. falciparum, but not P. vivax, merozoite IgG seroprevalence coincided with seasonal P. falciparum transmission peaks, even in those with no molecularly detectable parasites. At the individual level, antibody seropositivity was associated with reduced odds of contemporaneous P. falciparum (OR for PfCSP 0.51 [95%CI 0.35, 0.76], p = 0.001, PfAMA1 0.70 [95%CI 0.52, 0.93], p = 0.01, and PfMSP2 0.81 [95%CI 0.61, 1.08], p = 0.15), but not P. vivax infection (OR PvAMA1 1.02 [95%CI 0.73, 1.43], p = 0.89) indicating a potential role of immunity in protection against molecular-detectable P. falciparum parasitaemia. CONCLUSIONS: We demonstrated that integration and implementation of sample collection for molecular and serological surveillance into networks of VHV servicing hard-to-reach populations in the GMS is feasible, can capture significant levels of ongoing undetected seasonal malaria transmission and has the potential to supplement current routine RDT testing. Improving malaria surveillance by advancing the integration of molecular and serological techniques, through centralised testing approaches or novel point-of-contact tests, will advance progress, and tracking, towards malaria elimination goals in the GMS.

Evaluation of the effectiveness of topical repellent distributed by village health volunteer networks against Plasmodium spp. infection in Myanmar: A stepped-wedge cluster randomised trial.

BACKGROUND: The World Health Organization has yet to endorse deployment of topical repellents for malaria prevention as part of public health campaigns. We aimed to quantify the effectiveness of repellent distributed by the village health volunteer (VHV) network in the Greater Mekong Subregion (GMS) in reducing malaria in order to advance regional malaria elimination. METHODS AND FINDINGS: Between April 2015 and June 2016, a 15-month stepped-wedge cluster randomised trial was conducted in 116 villages in Myanmar (stepped monthly in blocks) to test the effectiveness of 12% N,N-diethylbenzamide w/w cream distributed by VHVs, on Plasmodium spp. infection. The median age of participants was 18 years, approximately half were female, and the majority were either village residents (46%) or forest dwellers (40%). No adverse events were reported during the study. Generalised linear mixed modelling estimated the effect of repellent on infection detected by rapid diagnostic test (RDT) (primary outcome) and polymerase chain reaction (PCR) (secondary outcome). Overall Plasmodium infection detected by RDT was low (0.16%; 50/32,194), but infection detected by PCR was higher (3%; 419/13,157). There was no significant protection against RDT-detectable infection (adjusted odds ratio [AOR] = 0.25, 95% CI 0.004-15.2, p = 0.512). In Plasmodium-species-specific analyses, repellent protected against PCR-detectable P. falciparum (adjusted relative risk ratio [ARRR] = 0.67, 95% CI 0.47-0.95, p = 0.026), but not P. vivax infection (ARRR = 1.41, 95% CI 0.80-2.47, p = 0.233). Repellent effects were similar when delayed effects were modelled, across risk groups, and regardless of village-level and temporal heterogeneity in malaria prevalence. The incremental cost-effectiveness ratio was US$256 per PCR-detectable infection averted. Study limitations were a lower than expected Plasmodium spp. infection rate and potential geographic dilution of the intervention. CONCLUSIONS: In this study, we observed apparent protection against new infections associated with the large-scale distribution of repellent by VHVs. Incorporation of repellent into national strategies, particularly in areas where bed nets are less effective, may contribute to the interruption of malaria transmission. Further studies are warranted across different transmission settings and populations, from the GMS and beyond, to inform WHO public health policy on the deployment of topical repellents for malaria prevention. TRIAL REGISTRATION: Australian and New Zealand Clinical Trials Registry (ACTRN12616001434482).

Monitoring Plasmodium falciparum and Plasmodium vivax using microsatellite markers indicates limited changes in population structure after substantial transmission decline in Papua New Guinea.

Monitoring the genetic structure of pathogen populations may be an economical and sensitive approach to quantify the impact of control on transmission dynamics, highlighting the need for a better understanding of changes in population genetic parameters as transmission declines. Here we describe the first population genetic analysis of two major human malaria parasites, Plasmodium falciparum (Pf) and Plasmodium vivax (Pv), following nationwide distribution of long-lasting insecticide-treated nets (LLINs) in Papua New Guinea (PNG). Parasite isolates from pre- (2005-2006) and post-LLIN (2010-2014) were genotyped using microsatellite markers. Despite parasite prevalence declining substantially (East Sepik Province: Pf = 54.9%-8.5%, Pv = 35.7%-5.6%, Madang Province: Pf = 38.0%-9.0%, Pv: 31.8%-19.7%), genetically diverse and intermixing parasite populations remained. Pf diversity declined modestly post-LLIN relative to pre-LLIN (East Sepik: Rs  = 7.1-6.4, HE  = 0.77-0.71; Madang: Rs  = 8.2-6.1, HE  = 0.79-0.71). Unexpectedly, population structure present in pre-LLIN populations was lost post-LLIN, suggesting that more frequent human movement between provinces may have contributed to higher gene flow. Pv prevalence initially declined but increased again in one province, yet diversity remained high throughout the study period (East Sepik: Rs  = 11.4-9.3, HE  = 0.83-0.80; Madang: Rs  = 12.2-14.5, HE  = 0.85-0.88). Although genetic differentiation values increased between provinces over time, no significant population structure was observed at any time point. For both species, a decline in multiple infections and increasing clonal transmission and significant multilocus linkage disequilibrium post-LLIN were positive indicators of impact on the parasite population using microsatellite markers. These parameters may be useful adjuncts to traditional epidemiological tools in the early stages of transmission reduction.

SNP barcodes provide higher resolution than microsatellite markers to measure Plasmodium vivax population genetics.

BACKGROUND: Genomic surveillance of malaria parasite populations has the potential to inform control strategies and to monitor the impact of interventions. Barcodes comprising large numbers of single nucleotide polymorphism (SNP) markers are accurate and efficient genotyping tools, however may need to be tailored to specific malaria transmission settings, since 'universal' barcodes can lack resolution at the local scale. A SNP barcode was developed that captures the diversity and structure of Plasmodium vivax populations of Papua New Guinea (PNG) for research and surveillance. METHODS: Using 20 high-quality P. vivax genome sequences from PNG, a total of 178 evenly spaced neutral SNPs were selected for development of an amplicon sequencing assay combining a series of multiplex PCRs and sequencing on the Illumina MiSeq platform. For initial testing, 20 SNPs were amplified in a small number of mono- and polyclonal P. vivax infections. The full barcode was then validated by genotyping and population genetic analyses of 94 P. vivax isolates collected between 2012 and 2014 from four distinct catchment areas on the highly endemic north coast of PNG. Diversity and population structure determined from the SNP barcode data was then benchmarked against that of ten microsatellite markers used in previous population genetics studies. RESULTS: From a total of 28,934,460 reads generated from the MiSeq Illumina run, 87% mapped to the PvSalI reference genome with deep coverage (median = 563, range 56-7586) per locus across genotyped samples. Of 178 SNPs assayed, 146 produced high-quality genotypes (minimum coverage = 56X) in more than 85% of P. vivax isolates. No amplification bias was introduced due to either polyclonal infection or whole genome amplification (WGA) of samples before genotyping. Compared to the microsatellite panels, the SNP barcode revealed greater variability in genetic diversity between populations and geographical population structure. The SNP barcode also enabled assignment of genotypes according to their geographic origins with a significant association between genetic distance and geographic distance at the sub-provincial level. CONCLUSIONS: High-throughput SNP barcoding can be used to map variation of malaria transmission dynamics at sub-national resolution. The low cost per sample and genotyping strategy makes the transfer of this technology to field settings highly feasible.

The epidemiology of Plasmodium falciparum and Plasmodium vivax in East Sepik Province, Papua New Guinea, pre- and post-implementation of national malaria control efforts.

BACKGROUND: In the past decade, national malaria control efforts in Papua New Guinea (PNG) have received renewed support, facilitating nationwide distribution of free long-lasting insecticidal nets (LLINs), as well as improvements in access to parasite-confirmed diagnosis and effective artemisinin-combination therapy in 2011-2012. METHODS: To study the effects of these intensified control efforts on the epidemiology and transmission of Plasmodium falciparum and Plasmodium vivax infections and investigate risk factors at the individual and household level, two cross-sectional surveys were conducted in the East Sepik Province of PNG; one in 2005, before the scale-up of national campaigns and one in late 2012-early 2013, after 2 rounds of LLIN distribution (2008 and 2011-2012). Differences between studies were investigated using Chi square (χ2), Fischer's exact tests and Student's t-test. Multivariable logistic regression models were built to investigate factors associated with infection at the individual and household level. RESULTS: The prevalence of P. falciparum and P. vivax in surveyed communities decreased from 55% (2005) to 9% (2013) and 36% to 6%, respectively. The mean multiplicity of infection (MOI) decreased from 1.8 to 1.6 for P. falciparum (p = 0.08) and from 2.2 to 1.4 for P. vivax (p < 0.001). Alongside these reductions, a shift towards a more uniform distribution of infections and illness across age groups was observed but there was greater heterogeneity across the study area and within the study villages. Microscopy positive infections and clinical cases in the household were associated with high rate infection households (> 50% of household members with Plasmodium infection). CONCLUSION: After the scale-up of malaria control interventions in PNG between 2008 and 2012, there was a substantial reduction in P. falciparum and P. vivax infection rates in the studies villages in East Sepik Province. Understanding the extent of local heterogeneity in malaria transmission and the driving factors is critical to identify and implement targeted control strategies to ensure the ongoing success of malaria control in PNG and inform the development of tools required to achieve elimination. In household-based interventions, diagnostics with a sensitivity similar to (expert) microscopy could be used to identify and target high rate households.

Re-emergence of yaws after single mass azithromycin treatment followed by targeted treatment: a longitudinal study.

BACKGROUND: Yaws is a substantial cause of chronic disfiguring ulcers in children in at least 14 countries in the tropics. WHO's newly adopted strategy for yaws eradication uses a single round of mass azithromycin treatment followed by targeted treatment programmes, and data from pilot studies have shown a short-term significant reduction of yaws. We assessed the long-term efficacy of the WHO strategy for yaws eradication. METHODS: Between April 15, 2013, and Oct 24, 2016, we did a longitudinal study on a Papua New Guinea island (Lihir; 16 092 population) in which yaws was endemic. In the initial study, the participants were followed for 12 months; in this extended follow-up study, clinical, serological, and PCR surveys were continued every 6 months for 42 months. We used genotyping and travel history to identify importation events. Active yaws confirmed by PCR specific for Treponema pallidum was the primary outcome indicator. The study is registered with ClinicalTrials.gov, number NCT01955252. FINDINGS: Mass azithromycin treatment (coverage rate of 84%) followed by targeted treatment programmes reduced the prevalence of active yaws from 1·8% to a minimum of 0·1% at 18 months (difference from baseline -1·7%, 95% CI, -1·9 to -1·4; p<0·0001), but the infection began to re-emerge after 24 months with a significant increase to 0·4% at 42 months (difference from 18 months 0·3%, 95% CI 0·1 to 0·4; p<0·0001). At each timepoint after baseline, more than 70% of the total community burden of yaws was found in individuals who had not had the mass treatment or as new infections in non-travelling residents. At months 36 and 42, five cases of active yaws, all from the same village, showed clinical failure following azithromycin treatment, with PCR-detected mutations in the 23S ribosomal RNA genes conferring resistance to azithromycin. A sustained decrease in the prevalence of high-titre latent yaws from 13·7% to <1·5% in asymptomatic children aged 1-5 years old and of genetic diversity of yaws strains from 0·139 to less than 0·046 between months 24 and 42 indicated a reduction in transmission of infection. INTERPRETATION: The implementation of the WHO strategy did not, in the long-term, achieve elimination in a high-endemic community mainly due to the individuals who were absent at the time of mass treatment in whom yaws reactivated; repeated mass treatment might be necessary to eliminate yaws. To our knowledge, this is the first report of the emergence of azithromycin-resistant T p pertenue and spread within one village. Communities' surveillance should be strengthened to detect any possible treatment failure and biological markers of resistance. FUNDING: ISDIN laboratories, Newcrest Mining Limited, and US Public Health Service National Institutes of Health.

Structurally conserved erythrocyte-binding domain in Plasmodium provides a versatile scaffold for alternate receptor engagement.

Understanding how malaria parasites gain entry into human red blood cells is essential for developing strategies to stop blood stage infection. Plasmodium vivax preferentially invades reticulocytes, which are immature red blood cells. The organism has two erythrocyte-binding protein families: namely, the Duffy-binding protein (PvDBP) and the reticulocyte-binding protein (PvRBP) families. Several members of the PvRBP family bind reticulocytes, specifically suggesting a role in mediating host cell selectivity of P. vivax. Here, we present, to our knowledge, the first high-resolution crystal structure of an erythrocyte-binding domain from PvRBP2a, solved at 2.12 Å resolution. The monomeric molecule consists of 10 α-helices and one short ß-hairpin, and, although the structural fold is similar to that of PfRh5--the essential invasion ligand in Plasmodium falciparum--its surface properties are distinct and provide a possible mechanism for recognition of alternate receptors. Sequence alignments of the crystallized fragment of PvRBP2a with other PvRBPs highlight the conserved placement of disulfide bonds. PvRBP2a binds mature red blood cells through recognition of an erythrocyte receptor that is neuraminidase- and chymotrypsin-resistant but trypsin-sensitive. By examining the patterns of sequence diversity within field isolates, we have identified and mapped polymorphic residues to the PvRBP2a structure. Using mutagenesis, we have also defined the critical residues required for erythrocyte binding. Characterization of the structural features that govern functional erythrocyte binding for the PvRBP family provides a framework for generating new tools that block P. vivax blood stage infection.

Antibodies to Intercellular Adhesion Molecule 1-Binding Plasmodium falciparum Erythrocyte Membrane Protein 1-DBLß Are Biomarkers of Protective Immunity to Malaria in a Cohort of Young Children from Papua New Guinea.

Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) mediates parasite sequestration to the cerebral microvasculature via binding of DBLß domains to intercellular adhesion molecule 1 (ICAM1) and is associated with severe cerebral malaria. In a cohort of 187 young children from Papua New Guinea (PNG), we examined baseline levels of antibody to the ICAM1-binding PfEMP1 domain, DBLß3PF11_0521, in comparison to four control antigens, including NTS-DBLα and CIDR1 domains from another group A variant and a group B/C variant. Antibody levels for the group A antigens were strongly associated with age and exposure. Antibody responses to DBLß3PF11_0521 were associated with a 37% reduced risk of high-density clinical malaria in the follow-up period (adjusted incidence risk ratio [aIRR] = 0.63 [95% confidence interval {CI}, 0.45 to 0.88; P = 0.007]) and a 25% reduction in risk of low-density clinical malaria (aIRR = 0.75 [95% CI, 0.55 to 1.01; P = 0.06]), while there was no such association for other variants. Children who experienced severe malaria also had significantly lower levels of antibody to DBLß3PF11_0521 and the other group A domains than those that experienced nonsevere malaria. Furthermore, a subset of PNG DBLß sequences had ICAM1-binding motifs, formed a distinct phylogenetic cluster, and were similar to sequences from other areas of endemicity. PfEMP1 variants associated with these DBLß domains were enriched for DC4 and DC13 head structures implicated in endothelial protein C receptor (EPCR) binding and severe malaria, suggesting conservation of dual binding specificities. These results provide further support for the development of specific classes of PfEMP1 as vaccine candidates and as biomarkers for protective immunity against clinical P. falciparum malaria.

Sustained Malaria Control Over an 8-Year Period in Papua New Guinea: The Challenge of Low-Density Asymptomatic Plasmodium Infections.

Background: The scale-up of effective malaria control in the last decade has resulted in a substantial decline in the incidence of clinical malaria in many countries. The effects on the proportions of asymptomatic and submicroscopic infections and on transmission potential are yet poorly understood. Methods: In Papua New Guinea, vector control has been intensified since 2008, and improved diagnosis and treatment was introduced in 2012. Cross-sectional surveys were conducted in Madang Province in 2006 (with 1280 survey participants), 2010 (with 2117 participants), and 2014 (with 2516 participants). Infections were quantified by highly sensitive quantitative polymerase chain reaction (PCR) analysis, and gametocytes were quantified by reverse-transcription qPCR analysis. Results: Plasmodium falciparum prevalence determined by qPCR decreased from 42% in 2006 to 9% in 2014. The P. vivax prevalence decreased from 42% in 2006 to 13% in 2010 but then increased to 20% in 2014. Parasite densities decreased 5-fold from 2006 to 2010; 72% of P. falciparum and 87% of P. vivax infections were submicroscopic in 2014. Gametocyte density and positivity correlated closely with parasitemia, and population gametocyte prevalence decreased 3-fold for P. falciparum and 29% for P. vivax from 2010 to 2014. Conclusions: Sustained control has resulted in reduced malaria transmission potential, but an increasing proportion of gametocyte carriers are asymptomatic and submicroscopic and represent a challenge to malaria control.

Phylogeography of var gene repertoires reveals fine-scale geospatial clustering of Plasmodium falciparum populations in a highly endemic area.

Plasmodium falciparum malaria is a major global health problem that is being targeted for progressive elimination. Knowledge of local disease transmission patterns in endemic countries is critical to these elimination efforts. To investigate fine-scale patterns of malaria transmission, we have compared repertoires of rapidly evolving var genes in a highly endemic area. A total of 3680 high-quality DBLα-sequences were obtained from 68 P. falciparum isolates from ten villages spread over two distinct catchment areas on the north coast of Papua New Guinea (PNG). Modelling of the extent of var gene diversity in the two parasite populations predicts more than twice as many var gene alleles circulating within each catchment (Mugil = 906; Wosera = 1094) than previously recognized in PNG (Amele = 369). In addition, there were limited levels of var gene sharing between populations, consistent with local parasite population structure. Phylogeographic analyses demonstrate that while neutrally evolving microsatellite markers identified population structure only at the catchment level, var gene repertoires reveal further fine-scale geospatial clustering of parasite isolates. The clustering of parasite isolates by village in Mugil, but not in Wosera was consistent with the physical and cultural isolation of the human populations in the two catchments. The study highlights the microheterogeneity of P. falciparum transmission in highly endemic areas and demonstrates the potential of var genes as markers of local patterns of parasite population structure.

Limited antigenic diversity of Plasmodium falciparum apical membrane antigen 1 supports the development of effective multi-allele vaccines.

BACKGROUND: Polymorphism in antigens is a common mechanism for immune evasion used by many important pathogens, and presents major challenges in vaccine development. In malaria, many key immune targets and vaccine candidates show substantial polymorphism. However, knowledge on antigenic diversity of key antigens, the impact of polymorphism on potential vaccine escape, and how sequence polymorphism relates to antigenic differences is very limited, yet crucial for vaccine development. Plasmodium falciparum apical membrane antigen 1 (AMA1) is an important target of naturally-acquired antibodies in malaria immunity and a leading vaccine candidate. However, AMA1 has extensive allelic diversity with more than 60 polymorphic amino acid residues and more than 200 haplotypes in a single population. Therefore, AMA1 serves as an excellent model to assess antigenic diversity in malaria vaccine antigens and the feasibility of multi-allele vaccine approaches. While most previous research has focused on sequence diversity and antibody responses in laboratory animals, little has been done on the cross-reactivity of human antibodies. METHODS: We aimed to determine the extent of antigenic diversity of AMA1, defined by reactivity with human antibodies, and to aid the identification of specific alleles for potential inclusion in a multi-allele vaccine. We developed an approach using a multiple-antigen-competition enzyme-linked immunosorbent assay (ELISA) to examine cross-reactivity of naturally-acquired antibodies in Papua New Guinea and Kenya, and related this to differences in AMA1 sequence. RESULTS: We found that adults had greater cross-reactivity of antibodies than children, although the patterns of cross-reactivity to alleles were the same. Patterns of antibody cross-reactivity were very similar between populations (Papua New Guinea and Kenya), and over time. Further, our results show that antigenic diversity of AMA1 alleles is surprisingly restricted, despite extensive sequence polymorphism. Our findings suggest that a combination of three different alleles, if selected appropriately, may be sufficient to cover the majority of antigenic diversity in polymorphic AMA1 antigens. Antigenic properties were not strongly related to existing haplotype groupings based on sequence analysis. CONCLUSIONS: Antigenic diversity of AMA1 is limited and a vaccine including a small number of alleles might be sufficient for coverage against naturally-circulating strains, supporting a multi-allele approach for developing polymorphic antigens as malaria vaccines.

Distinct patterns of diversity, population structure and evolution in the AMA1 genes of sympatric Plasmodium falciparum and Plasmodium vivax populations of Papua New Guinea from an area of similarly high transmission.

BACKGROUND: As Plasmodium falciparum and Plasmodium vivax co-exist in most malaria-endemic regions outside sub-Saharan Africa, malaria control strategies in these areas must target both species in order to succeed. Population genetic analyses can predict the effectiveness of interventions including vaccines, by providing insight into patterns of diversity and evolution. The aim of this study was to investigate the population genetics of leading malaria vaccine candidate AMA1 in sympatric P. falciparum and P. vivax populations of Papua New Guinea (PNG), an area of similarly high prevalence (Pf = 22.3 to 38.8%, Pv = 15.3 to 31.8%). METHODS: A total of 72 Pfama1 and 102 Pvama1 sequences were collected from two distinct areas, Madang and Wosera, on the highly endemic PNG north coast. RESULTS: Despite a greater number of polymorphic sites in the AMA1 genes of P. falciparum (Madang = 52; Wosera = 56) compared to P. vivax (Madang = 36, Wosera = 34), the number of AMA1 haplotypes, haplotype diversity (Hd) and recombination (R) was far lower for P. falciparum (Madang = 12, Wosera = 20; Hd ≤0.92, R ≤45.8) than for P. vivax (Madang = 50, Wosera = 38; Hd = 0.99, R = ≤70.9). Balancing selection was detected only within domain I of AMA1 for P. vivax, and in both domains I and III for P. falciparum. CONCLUSIONS: Higher diversity in the genes encoding P. vivax AMA1 than in P. falciparum AMA1 in this highly endemic area has important implications for development of AMA1-based vaccines in PNG and beyond. These results also suggest a smaller effective population size of P. falciparum compared to P. vivax, a finding that warrants further investigation. Differing patterns of selection on the AMA1 genes indicate that critical antigenic sites may differ between the species, highlighting the need for independent investigations of these two leading vaccine candidates.