Population genetic analysis of the Plasmodium falciparum 6-cys protein Pf38 in Papua New Guinea reveals domain-specific balancing selection

<p>Abstract</p> <p>Background</p> <p>The <it>Plasmodium falciparum </it>merozoite surface protein Pf38 is targeted by antibodies of malaria immune adults and has been shown to be under balancing (immune) selection in a Gambian parasite population, indicating potential as a malaria vaccine candidate. This study explores the population genetics of <it>Pf</it>38 in Papua New Guinea, to determine the extent and geographic distribution of diversity and to measure selective pressure along the length of the gene.</p> <p>Methods</p> <p>Using samples collected during community-based cross-sectional surveys in the Mugil and Wosera regions, the <it>Pf38 </it>genes of 59 <it>P. falciparum </it>isolates were amplified and sequenced. These sequences, along with previously sequenced Gambian and laboratory isolates, were then subjected to an array of population genetic analyses, examining polymorphisms, haplotype diversity and balancing selection. In addition to whole-gene analysis, the two 6-cys domains were considered separately, to investigate domain specific polymorphism and selection.</p> <p>Results</p> <p>Nineteen polymorphic sites were identified in the <it>Pf </it>38 gene. Of these, 13 were found in the Gambia, 10 in Mugil and 8 in Wosera. Notably, the majority of common polymorphisms were confined to domain I. Although only moderate levels of nucleotide diversity were observed, the haplotype diversity was high in all populations, suggesting extensive recombination. Analyses of the full-length sequence provided only modest evidence for balancing selection. However, there was a strong contrast between domain I, which showed strong evidence for positive balancing selection, and domain II which was neutral. Analyses of the geographic distribution of Pf38 haplotypes showed that four haplotypes accounted for the majority of sequences found world-wide, but there were many more haplotypes unique to the African than the PNG populations.</p> <p>Conclusion</p> <p>This study confirmed previous findings that <it>Pf38 </it>is a polymorphic gene under balancing selection. However, analysing polymorphism and selection across the length of the gene painted a considerably different picture. Domain I is highly polymorphic and the target of significant balancing selection. In contrast, domain II is relatively conserved and does not show evidence of immune selective pressure. The findings have implications for future population genetic studies on vaccine candidates, showing that the biological context must also be considered as a framework for analysis.</p

Correction: Population Genomics of the Immune Evasion (var) Genes of Plasmodium falciparum

Var genes encode the major surface antigen (PfEMP1) of the blood stages of the human malaria parasite Plasmodium falciparum. Differential expression of up to 60 diverse var genes in each parasite genome underlies immune evasion. We compared the diversity of the DBLalpha domain of var genes sampled from 30 parasite isolates from a malaria endemic area of Papua New Guinea (PNG) and 59 from widespread geographic origins (global). Overall, we obtained over 8,000 quality-controlled DBLalpha sequences. Within our sampling frame, the global population had a total of 895 distinct DBLalpha "types" and negligible overlap among repertoires. This indicated that var gene diversity on a global scale is so immense that many genomes would need to be sequenced to capture its true extent. In contrast, we found a much lower diversity in PNG of 185 DBLalpha types, with an average of approximately 7% overlap among repertoires. While we identify marked geographic structuring, nearly 40% of types identified in PNG were also found in samples from different countries showing a cosmopolitan distribution for much of the diversity. We also present evidence to suggest that recombination plays a key role in maintaining the unprecedented levels of polymorphism found in these immune evasion genes. This population genomic framework provides a cost effective molecular epidemiological tool to rapidly explore the geographic diversity of var genes

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

BackgroundPolymorphism 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.MethodsWe 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.ResultsWe 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.ConclusionsAntigenic 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

Blood-Stage Parasitaemia and Age Determine Plasmodium falciparum and P. vivax Gametocytaemia in Papua New Guinea

A better understanding of human-to-mosquito transmission is crucial to control malaria. In order to assess factors associated with gametocyte carriage, 2083 samples were collected in a cross-sectional survey in Papua New Guinea. Plasmodium species were detected by light microscopy and qPCR and gametocytes by detection of pfs25 and pvs25 mRNA transcripts by reverse-transcriptase PCR (qRT-PCR). The parasite prevalence by PCR was 18.5% for Plasmodium falciparum and 13.0% for P. vivax. 52.5% of all infections were submicroscopic. Gametocytes were detected in 60% of P. falciparum-positive and 51% of P. vivax-positive samples. Each 10-fold increase in parasite density led to a 1.8-fold and 3.3-fold increase in the odds of carrying P. falciparum and P. vivax gametocytes. Thus the proportion of gametocyte positive and gametocyte densities was highest in young children carrying high asexual parasite densities and in symptomatic individuals. Dilution series of gametocytes allowed absolute quantification of gametocyte densities by qRT-PCR and showed that pvs25 expression is 10-20 fold lower than pfs25 expression. Between 2006 and 2010 parasite prevalence in the study site has decreased by half. 90% of the remaining infections were asymptomatic and likely constitute an important reservoir of transmission. However, mean gametocyte densities were low (approx. 1-2 gametocyte/muL) and it remains to be determined to what extent low-density gametocyte positive individuals are infective to mosquitos

Multilocus haplotypes reveal variable levels of diversity and population structure of Plasmodium falciparum in Papua New Guinea, a region of intense perennial transmission

<p>Abstract</p> <p>Background</p> <p>The South West Pacific nation of Papua New Guinea has intense year round transmission of <it>Plasmodium falciparum </it>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 <it>P. falciparum </it>in the region. This study describes the population genetics of <it>P. falciparum </it>in thirteen villages spread over four distinct catchment areas of Papua New Guinea.</p> <p>Methods</p> <p>Ten microsatellite loci were genotyped in 318 <it>P. falciparum </it>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.</p> <p>Results</p> <p>Although 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 (<it>F_ST</it>) 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.</p> <p>Conclusion</p> <p>The observed population genetics of <it>P. falciparum </it>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 <it>P. falciparum </it>has important implications for malaria control strategies and warrants further fine scale sampling throughout Papua New Guinea.</p

Population Genomics of the Immune Evasion (var) Genes of Plasmodium falciparum

Var genes encode the major surface antigen (PfEMP1) of the blood stages of the human malaria parasite Plasmodium falciparum. Differential expression of up to 60 diverse var genes in each parasite genome underlies immune evasion. We compared the diversity of the DBLα domain of var genes sampled from 30 parasite isolates from a malaria endemic area of Papua New Guinea (PNG) and 59 from widespread geographic origins (global). Overall, we obtained over 8,000 quality-controlled DBLα sequences. Within our sampling frame, the global population had a total of 895 distinct DBLα “types” and negligible overlap among repertoires. This indicated that var gene diversity on a global scale is so immense that many genomes would need to be sequenced to capture its true extent. In contrast, we found a much lower diversity in PNG of 185 DBLα types, with an average of approximately 7% overlap among repertoires. While we identify marked geographic structuring, nearly 40% of types identified in PNG were also found in samples from different countries showing a cosmopolitan distribution for much of the diversity. We also present evidence to suggest that recombination plays a key role in maintaining the unprecedented levels of polymorphism found in these immune evasion genes. This population genomic framework provides a cost effective molecular epidemiological tool to rapidly explore the geographic diversity of var genes

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

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  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

Increasingly inbred and fragmented populations of Plasmodium vivax associated with the eastward decline in malaria transmission across the Southwest Pacific

The human malaria parasite Plasmodium vivax is more resistant to malaria control strategies than Plasmodium falciparum, and maintains high genetic diversity even when transmission is low. To investigate whether declining P. vivax transmission leads to increasing population structure that would facilitate elimination, we genotyped samples from across the Southwest Pacific region, which experiences an eastward decline in malaria transmission, as well as samples from two time points at one site (Tetere, Solomon Islands) during intensified malaria control. Analysis of 887 P. vivax microsatellite haplotypes from hyperendemic Papua New Guinea (PNG, n = 443), meso-hyperendemic Solomon Islands (n = 420), and hypoendemic Vanuatu (n = 24) revealed increasing population structure and multilocus linkage disequilibrium yet a modest decline in diversity as transmission decreases over space and time. In Solomon Islands, which has had sustained control efforts for 20 years, and Vanuatu, which has experienced sustained low transmission for many years, significant population structure was observed at different spatial scales. We conclude that control efforts will eventually impact P. vivax population structure and with sustained pressure, populations may eventually fragment into a limited number of clustered foci that could be targeted for elimination