The widespread distribution and relapsing nature of Plasmodium vivax infection present major challenges for malaria elimination. To characterise the genetic diversity of this parasite within individual infections and across the population, we performed deep genome sequencing of >200 clinical samples collected across the Asia-Pacific region, and analysed data on >300,000 SNPs and 9 regions of the genome with large copy number variations. Individual infections showed complex patterns of genetic structure, with variation not only in the number of dominant clones but also in their level of relatedness and inbreeding. At the population level, we observed strong signals of recent evolutionary selection both in known drug resistance genes and at novel loci, and these varied markedly between geographical locations. These findings reveal a dynamic landscape of local evolutionary adaptation in P. vivax populations, and provide a foundation for genomic surveillance to guide effective strategies for control and elimination.

The EBAs appear important targets of acquired protective immunity. These findings support their further development as vaccine candidates.

BackgroundThe South West Pacific nation of Papua New Guinea has intense year round transmission of Plasmodium falciparum on the coast and in the low-lying inland areas. Local heterogeneity in the epidemiology of malaria suggests that parasites from multiple locations will need to be surveyed to define the population biology of P. falciparum in the region. This study describes the population genetics of P. falciparum in thirteen villages spread over four distinct catchment areas of Papua New Guinea.MethodsTen microsatellite loci were genotyped in 318 P. falciparum isolates from the parasite populations of two inland catchment areas, namely Wosera (number of villages (n) = 7) and Utu (n = 1) and; and two coastal catchments, Malala (n = 3) and Mugil (n = 3). Analysis of the resultant multilocus haplotypes was done at different spatial scales (2-336 km) to define the genetic diversity (allelic richness and expected heterozygosity), linkage disequilibrium and population structure throughout the study area.ResultsAlthough genetic diversity was high in all parasite populations, it was also variable with a lower allelic richness and expected heterozygosity for inland populations compared to those from the more accessible coast. This variability was not correlated with two proxy measures of transmission intensity, the infection prevalence and the proportion multiple infections. Random associations among the microsatellite loci were observed in all four catchments showing that a substantial degree of out-crossing occurs in the region. Moderate to very high levels of population structure were found but the amount of genetic differentiation (FST) did not correlate with geographic distance suggesting that parasite populations are fragmented. Population structure was also identified between villages within the Malala area, with the haplotypes of one parasite population clustering with the neighbouring catchment of Mugil.ConclusionThe observed population genetics of P. falciparum in this region is likely to be a consequence of the high transmission intensity combined with the isolation of human and vector populations, especially those located inland and migration of parasites via human movement into coastal populations. The variable genetic diversity and population structure of P. falciparum has important implications for malaria control strategies and warrants further fine scale sampling throughout Papua New Guinea.

IntroductionThe human malaria parasite, Plasmodium vivax, is proving more difficult to control and eliminate than Plasmodium falciparum in areas of co-transmission. Comparisons of the genetic structure of sympatric parasite populations may provide insight into the mechanisms underlying the resilience of P. vivax and can help guide malaria control programs.Methodology/Principle findings P. vivax isolates representing the parasite populations of four areas on the north coast of Papua New Guinea (PNG) were genotyped using microsatellite markers and compared with previously published microsatellite data from sympatric P. falciparum isolates. The genetic diversity of P. vivax (He = 0.83–0.85) was higher than that of P. falciparum (He = 0.64–0.77) in all four populations. Moderate levels of genetic differentiation were found between P. falciparum populations, even over relatively short distances (less than 50 km), with 21–28% private alleles and clear geospatial genetic clustering. Conversely, very low population differentiation was found between P. vivax catchments, with less than 5% private alleles and no genetic clustering observed. In addition, the effective population size of P. vivax (30353; 13043–69142) was larger than that of P. falciparum (18871; 8109–42986).Conclusions/SignificanceDespite comparably high prevalence, P. vivax had higher diversity and a panmictic population structure compared to sympatric P. falciparum populations, which were fragmented into subpopulations. The results suggest that in comparison to P. falciparum, P. vivax has had a long-term large effective population size, consistent with more intense and stable transmission, and limited impact of past control and elimination efforts. This underlines suggestions that more intensive and sustained interventions will be needed to control and eventually eliminate P. vivax. This research clearly demonstrates how population genetic analyses can reveal deeper insight into transmission patterns than traditional surveillance methods.

Background: The protection afforded by human erythrocyte polymorphisms against the malaria parasite, Plasmodium falciparum, has been proposed to be due to reduced ability of the parasite to invade or develop in erythrocytes. If this were the case, variable levels of parasitaemia and rates of seroconversion to infected-erythrocyte variant surface antigens (VSA) should be seen in different host genotypes.

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.

Abs targeting blood-stage Ags of Plasmodium falciparum are important in acquired immunity to malaria, but major targets remain unclear. The P. falciparum reticulocyte-binding homologs (PfRh) are key ligands used by merozoites during invasion of erythrocytes. PfRh2a and PfRh2b are functionally important members of this family and may be targets of protective immunity, but their potential role in human immunity has not been examined. We expressed eight recombinant proteins covering the entire PfRh2 common region, as well as PfRh2a- and PfRh2b-specific regions. Abs were measured among a cohort of 206 Papua New Guinean children who were followed prospectively for 6 mo for reinfection and malaria. At baseline, Abs were associated with increasing age and active infection. High levels of IgG to all PfRh2 protein constructs were strongly associated with protection from symptomatic malaria and high-density parasitemia. The predominant IgG subclasses were IgG1 and IgG3, with little IgG2 and IgG4 detected. To further understand the significance of PfRh2 as an immune target, we analyzed PfRh2 sequences and found that polymorphisms are concentrated in an N-terminal region of the protein and seem to be under diversifying selection, suggesting immune pressure. Cluster analysis arranged the sequences into two main groups, suggesting that many of the haplotypes identified may be antigenically similar. These findings provide evidence suggesting that PfRh2 is an important target of protective immunity in humans and that Abs act by controlling blood-stage parasitemia and support its potential for vaccine development.

No abstract