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

Microsatellites reveal high polymorphism and high potential for use in anti-malarial efficacy studies in areas with different transmission intensities in mainland Tanzania

BACKGROUND

Tanzania is currently implementing therapeutic efficacy studies (TES) in areas of varying malaria transmission intensities as per the World Health Organization (WHO) recommendations. In TES, distinguishing reinfection from recrudescence is critical for the determination of anti-malarial efficacy. Recently, the WHO recommended genotyping polymorphic coding genes, merozoite surface proteins 1 and 2 (msp1 and msp2), and replacing the glutamate-rich protein (glurp) gene with one of the highly polymorphic microsatellites in Plasmodium falciparum to adjust the efficacy of antimalarials in TES. This study assessed the polymorphisms of six neutral microsatellite markers and their potential use in TES, which is routinely performed in Tanzania.

METHODS

Plasmodium falciparum samples were obtained from four TES sentinel sites, Kibaha (Pwani), Mkuzi (Tanga), Mlimba (Morogoro) and Ujiji (Kigoma), between April and September 2016. Parasite genomic DNA was extracted from dried blood spots on filter papers using commercial kits. Genotyping was done using six microsatellites (Poly-α, PfPK2, TA1, C3M69, C2M34 and M2490) by capillary method, and the data were analysed to determine the extent of their polymorphisms and genetic diversity at the four sites.

RESULTS

Overall, 83 (88.3%) of the 94 samples were successfully genotyped (with positive results for ≥ 50.0% of the markers), and > 50.0% of the samples (range = 47.6-59.1%) were polyclonal, with a mean multiplicity of infection (MOI) ranging from 1.68 to 1.88 among the four sites. There was high genetic diversity but limited variability among the four sites based on mean allelic richness (RS = 7.48, range = 7.27-8.03, for an adjusted minimum sample size of 18 per site) and mean expected heterozygosity (He = 0.83, range = 0.80-0.85). Cluster analysis of haplotypes using STRUCTURE, principal component analysis, and pairwise genetic differentiation (FST) did not reveal population structure or clustering of parasites according to geographic origin. Of the six markers, Poly-α was the most polymorphic, followed by C2M34, TA1 and C3M69, while M2490 was the least polymorphic.

CONCLUSION

Microsatellite genotyping revealed high polyclonality and genetic diversity but no significant population structure. Poly-α, C2M34, TA1 and C3M69 were the most polymorphic markers, and Poly-α alone or with any of the other three markers could be adopted for use in TES in Tanzania.

Measurably recombining malaria parasites

Genomic epidemiology has guided research and policy for various viral pathogens and there has been a parallel effort towards using genomic epidemiology to combat diseases that are caused by eukaryotic pathogens, such as the malaria parasite. However, the central concept of viral genomic epidemiology, namely that of measurably mutating pathogens, does not apply easily to sexually recombining parasites. Here we introduce the related but different concept of measurably recombining malaria parasites to promote convergence around a unifying theoretical framework for malaria genomic epidemiology. Akin to viral phylodynamics, we anticipate that an inferential framework developed around recombination will help guide practical research and thus realize the full public health potential of genomic epidemiology for malaria parasites and other sexually recombining pathogens.

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

Distinguishing gene flow between malaria parasite populations

Measuring gene flow between malaria parasite populations in different geographic locations can provide strategic information for malaria control interventions. Multiple important questions pertaining to the design of such studies remain unanswered, limiting efforts to operationalize genomic surveillance tools for routine public health use. This report examines the use of population-level summaries of genetic divergence (FST) and relatedness (identity-by-descent) to distinguish levels of gene flow between malaria populations, focused on field-relevant questions about data size, sampling, and interpretability of observations from genomic surveillance studies. To do this, we use P. falciparum whole genome sequence data and simulated sequence data approximating malaria populations evolving under different current and historical epidemiological conditions. We employ mobile-phone associated mobility data to estimate parasite migration rates over different spatial scales and use this to inform our analysis. This analysis underscores the complementary nature of divergence- and relatedness-based metrics for distinguishing gene flow over different temporal and spatial scales and characterizes the data requirements for using these metrics in different contexts. Our results have implications for the design and implementation of malaria genomic surveillance studies.

Genetic diversity and population structure of the human malaria parasite Plasmodium falciparum surface protein Pfs47 in isolates from the lowlands in Western Kenya

Plasmodium falciparum parasites have evolved genetic adaptations to overcome immune responses mounted by diverse Anopheles vectors hindering malaria control efforts. Plasmodium falciparum surface protein Pfs47 is critical in the parasite’s survival by manipulating the vector’s immune system hence a promising target for blocking transmission in the mosquito. This study aimed to examine the genetic diversity, haplotype distribution, and population structure of Pfs47 and its implications on malaria infections in endemic lowlands in Western Kenya. Cross-sectional mass blood screening was conducted in malaria endemic regions in the lowlands of Western Kenya: Homa Bay, Kombewa, and Chulaimbo. Dried blood spots and slide smears were simultaneously collected in 2018 and 2019. DNA was extracted using Chelex method from microscopic Plasmodium falciparum positive samples and used to genotype Pfs47 using polymerase chain reaction (PCR) and DNA sequencing. Thirteen observed haplotypes of the Pfs47 gene were circulating in Western Kenya. Population-wise, haplotype diversity ranged from 0.69 to 0.77 and the nucleotide diversity 0.10 to 0.12 across all sites. All the study sites displayed negative Tajima’s D values although not significant. However, the negative and significant Fu’s Fs statistical values were observed across all the study sites, suggesting population expansion or positive selection. Overall genetic differentiation index was not significant (F_ST = -0.00891, P > 0.05) among parasite populations. All Nm values revealed a considerable gene flow in these populations. These results could have important implications for the persistence of high levels of malaria transmission and should be considered when designing potential targeted control interventions.

Genetic diversity and population structure of Plasmodium falciparum in Nigeria: insights from microsatellite loci analysis

BACKGROUND

Malaria remains a public health burden especially in Nigeria. To develop new malaria control and elimination strategies or refine existing ones, understanding parasite population diversity and transmission patterns is crucial.

METHODS

In this study, characterization of the parasite diversity and structure of Plasmodium falciparum isolates from 633 dried blood spot samples in Nigeria was carried out using 12 microsatellite loci of P. falciparum. These microsatellite loci were amplified via semi-nested polymerase chain reaction (PCR) and fragments were analysed using population genetic tools.

RESULTS

Estimates of parasite genetic diversity, such as mean number of different alleles (13.52), effective alleles (7.13), allelic richness (11.15) and expected heterozygosity (0.804), were high. Overall linkage disequilibrium was weak (0.006, P < 0.001). Parasite population structure was low (Fst: 0.008-0.105, AMOVA: 0.039).

CONCLUSION

The high level of parasite genetic diversity and low population structuring in this study suggests that parasite populations circulating in Nigeria are homogenous. However, higher resolution methods, such as the 24 SNP barcode and whole genome sequencing, may capture more specific parasite genetic signatures circulating in the country. The results obtained can be used as a baseline for parasite genetic diversity and structure, aiding in the formulation of appropriate therapeutic and control strategies in Nigeria.

Population structure and diversity of Plasmodium falciparum in children with asymptomatic malaria living in different ecological zones of Ghana

BACKGROUND

Genetic diversity in Plasmodium falciparum populations can be used to describe the resilience and spatial distribution of the parasite in the midst of intensified intervention efforts. This study used microsatellite analysis to evaluate the genetic diversity and population dynamics of P. falciparum parasites circulating in three ecological zones of Ghana.

METHODS

A total of 1168 afebrile children aged between 3 to 13 years were recruited from five (5) Primary schools in 3 different ecological zones (Sahel (Tamale and Kumbungu), Forest (Konongo) and Coastal (Ada and Dodowa)) of Ghana. Asymptomatic malaria parasite carriage was determined using microscopy and PCR, whilst fragment analysis of 6 microsatellite loci was used to determine the diversity and population structure of P. falciparum parasites.

RESULTS

Out of the 1168 samples examined, 16.1 and 39.5% tested positive for P. falciparum by microscopy and nested PCR respectively. The genetic diversity of parasites in the 3 ecological zones was generally high, with an average heterozygosity (He) of 0.804, 0.787 and 0.608 the rainy (peak) season for the Sahel, Forest and Coastal zones respectively. The mean He for the dry (off-peak) season were 0.562, 0.693 and 0.610 for the Sahel, Forest and Coastal zones respectively. Parasites from the Forest zone were more closely related to those from the Sahel than from the Coastal zone, despite the Coastal zone being closer in physical distance to the Forest zone. The fixation indexes among study sites ranged from 0.049 to 0.112 during the rainy season and 0.112 to 0.348 during the dry season.

CONCLUSION

A large asymptomatic parasite reservoir was found in the school children during both rainy and dry seasons, especially those in the Forest and Sahel savannah zones where parasites were also found to be related compared to those from the Coastal zone. Further studies are recommended to understand why despite the roll out of several malaria interventions in Ghana, high transmission still persist.

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: R_s  = 7.1-6.4, H_E  = 0.77-0.71; Madang: R_s  = 8.2-6.1, H_E  = 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: R_s  = 11.4-9.3, H_E  = 0.83-0.80; Madang: R_s  = 12.2-14.5, H_E  = 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.

Allelic variation of msp-3α gene in Plasmodium vivax isolates and its correlation with the severity of disease in vivax malaria

Malaria is a global socio-economic burden of which Plasmodium vivax contributes for about 70-80 million cases on an annual basis worldwide and 60-65% cases in India. Diversity observed in highly polymorphic Merozoite Surface Protein-3α (msp-3α) encoded by MSP-3 gene family, has been used efficiently for genotyping of P. vivax infection. This study aims to correlate the severity of clinical symptoms with parasite load, genotype of P. vivax and multiplicity of infection. Based on clinical symptoms classification, 31 (67.9%) out of 46 cases were found to be severe while 15 (32.6%) were non-severe and correlation of the severity of vivax infection with parasite load was not observed. Analysis of msp3-α allele genotype showed that out of 31 severe cases, 19 (61.2%) were single-clone infection cases whereas 12 (38.7%) were multi-clone infections. Similarly, out of 15 non-severe cases, 9 (60%) were single clone and 6 (40%) were multi-clone infections indicating the absence of a correlation between the multiplicity of infection and disease severity. Allele frequency observed was 65.9%, 23.4%, 23.4%, and 28.2% for allele A, B, C and D, respectively. An important finding was the greater distribution of allele D than alleles B and C, which has been reported as a rare allele otherwise. Further, of 13 cases with allele D, 76.9% (10/13) cases were severe. This study showed the absence of a correlation between the severity of clinical symptoms with parasite load and multiplicity of infection but at the same time drives a possibility of severe vivax malarial symptoms to have an association with the persistence of allele D in the population. This upon exploration can lead to the development of a target in detection of severe cases of malaria.

Confirmation of the absence of local transmission and geographic assignment of imported falciparum malaria cases to China using microsatellite panel

Background

Current methods to classify local and imported malaria infections depend primarily on patient travel history, which can have limited accuracy. Genotyping has been investigated as a complementary approach to track the spread of malaria and identify the origin of imported infections.

Methods

An extended panel of 26 microsatellites (16 new microsatellites) for Plasmodium falciparum was evaluated in 602 imported infections from 26 sub-Saharan African countries to the Jiangsu Province of People’s Republic of China. The potential of the 26 microsatellite markers to assign imported parasites to their geographic origin was assessed using a Bayesian method with Markov Chain Monte Carlo (MCMC) as implemented in the program Smoothed and Continuous Assignments (SCAT) with a modification to incorporate haploid genotype data.

Results

The newly designed microsatellites were polymorphic and are not in linkage disequilibrium with the existing microsatellites, supporting previous findings of high rate of recombination in sub-Saharan Africa. Consistent with epidemiology inferred from patients’ travel history, no evidence for local transmission was found; nearly all genetically related infections were identified in people who travelled to the same country near the same time. The smoothing assignment method assigned imported cases to their likely geographic origin with an accuracy (Angola: 59%; Nigeria: 51%; Equatorial Guinea: 40%) higher than would be achieved at random, reaching statistical significance for Angola and Equatorial Guinea.

Conclusions

Genotyping using an extended microsatellite panel is valuable for malaria case classification and programme evaluation in an elimination setting. A Bayesian method for assigning geographic origin of mammals based on genetic data was adapted for malaria and showed potential for identification of the origin of imported infections.

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 (χ²), 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.

Diversity and Multiplicity of P. falciparum infections among asymptomatic school children in Mbita, Western Kenya

Multiplicity of infection (MOI) and genetic diversity of P. falciparum infections are important surrogate indicators for assessing malaria transmission intensity in different regions of endemicity. Determination of MOI and diversity of P. falciparum among asymptomatic carriers will enhance our understanding of parasite biology and transmission to mosquito vectors. This study examined the MOI and genetic diversity of P. falciparum parasite populations circulating in Mbita, a region characterized as one of the malaria hotspots in Kenya. The genetic diversity and multiplicity of P. falciparum infections in 95 asymptomatic school children (age 5–15 yrs.) residing in Mbita, western Kenya were assessed using 10 polymorphic microsatellite markers. An average of 79.69% (Range: 54.84–95.74%) of the isolates analysed in this study were polyclonal infections as detected in at least one locus. A high mean MOI of 3.39 (Range: 2.24–4.72) and expected heterozygosity (He) of 0.81 (Range: 0.57–0.95) was reported in the study population. The analysed samples were extensively polyclonal infections leading to circulation of highly genetically diverse parasite populations in the study area. These findings correlated with the expectations of high malaria transmission intensity despite scaling up malaria interventions in the area thereby indicating the need for a robust malaria interventions particularly against asymptomatic carriers in order to attain elimination in the region.

Malaria Molecular Epidemiology: An Evolutionary Genetics Perspective

Malaria is a vector-borne disease that involves multiple parasite species in a variety of ecological settings. However, the parasite species causing the disease, the prevalence of subclinical infections, the emergence of drug resistance, the scale-up of interventions, and the ecological factors affecting malaria transmission, among others, are aspects that vary across areas where malaria is endemic. Such complexities have propelled the study of parasite genetic diversity patterns in the context of epidemiologic investigations. Importantly, molecular studies indicate that the time and spatial distribution of malaria cases reflect epidemiologic processes that cannot be fully understood without characterizing the evolutionary forces shaping parasite population genetic patterns. Although broad in scope, this review in the Microbiology Spectrum Curated Collection: Advances in Molecular Epidemiology highlights the need for understanding population genetic concepts when interpreting parasite molecular data. First, we discuss malaria complexity in terms of the parasite species involved. Second, we describe how molecular data are changing our understanding of malaria incidence and infectiousness. Third, we compare different approaches to generate parasite genetic information in the context of epidemiologically relevant questions related to malaria control. Finally, we describe a few Plasmodium genomic studies as evidence of how these approaches will provide new insights into the malaria disease dynamics. *This article is part of a curated collection.

Genetic diversity and population structure of Plasmodium falciparum in Kenyan-Ugandan border areas

Kenya has, in the last decade, made tremendous progress in the fight against malaria. Nevertheless, continued surveillance of the genetic diversity and population structure of Plasmodium falciparum is required to refine malaria control and to adapt and improve elimination strategies. Twelve neutral microsatellite loci were genotyped in 201 P. falciparum isolates obtained from the Kenyan-Ugandan border (Busia) and from two inland malaria-endemic sites situated in western (Nyando) and coastal (Msambweni) Kenya. Analyses were done to assess the genetic diversity (allelic richness and expected heterozygosity, [H_e ]), multilocus linkage disequilibrium ( I S A ) and population structure. A similarly high degree of genetic diversity was observed among the three parasite populations surveyed (mean H_e  = 0.76; P > 0.05). Except in Msambweni, random association of microsatellite loci was observed, indicating high parasite out-breeding. Low to moderate genetic structure (F_ST   = 0.022-0.076; P < 0.0001) was observed with only 5% variance in allele frequencies observed among the populations. This study shows that the genetic diversity of P. falciparum populations at the Kenyan-Ugandan border is comparable to the parasite populations from inland Kenya. In addition, high genetic diversity, panmixia and weak population structure in this study highlight the fitness of Kenyan P. falciparum populations to successfully withstand malaria control interventions.

Genetic Diversity of Plasmodium vivax Causing Epidemic Malaria in the Republic of Korea

Plasmodium vivax is more challenging to control and eliminate than P. falciparum due to its more asymptomatic infections with low parasite densities making diagnosis more difficult, in addition to its unique biological characteristics. The potential re-introduction of incidence cases, either through borders or via human migrations, is another major hurdle to sustained control and elimination. The Republic of Korea has experienced re-emergence of vivax malaria in 1993 but is one of the 32 malaria-eliminating countries to-date. Despite achieving successful nationwide control and elimination of vivax malaria, the evolutionary characteristics of vivax malaria isolates in the Republic of Korea have not been fully understood. In this review, we present an overview of the genetic variability of such isolates to increase understanding of the epidemiology, diversity, and dynamics of vivax populations in the Republic of Korea.

Lack of significant recovery of chloroquine sensitivity in Plasmodium falciparum parasites following discontinuance of chloroquine use in Papua New Guinea

Background

Chloroquine treatment for Plasmodium falciparum has been discontinued in almost all endemic regions due to the spread of resistant isolates. Reversal of chloroquine susceptibility after chloroquine discontinuation has been reported in dozens of endemic regions. However, this phenomenon has been mostly observed in Africa and is not well documented in other malaria endemic regions. To investigate this, an ex vivo study on susceptibility to chloroquine and lumefantrine was conducted during 2016–2018 in Wewak, Papua New Guinea where chloroquine had been removed from the official malaria treatment regimen in 2010. Genotyping of pfcrt and pfmdr1 was also performed.

Results

In total, 368 patients were enrolled in this study. Average IC₅₀ values for chloroquine were 106.6, 80.5, and 87.6 nM in 2016, 2017, and 2018, respectively. These values were not significantly changed from those obtained in 2002/2003 (108 nM). The majority of parasites harboured a pfcrt K76T the mutation responsible for chloroquine resistance. However, a significant upward trend was observed in the frequency of the K76 (wild) allele from 2.3% in 2016 to 11.7% in 2018 (P = 0.008; Cochran–Armitage trend test).

Conclusions

Eight years of chloroquine withdrawal has not induced a significant recovery of susceptibility in Papua New Guinea. However, an increasing tendency of parasites harbouring chloroquine-susceptible K76 suggests a possibility of resurgence of chloroquine susceptibility in the future.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2585-x) contains supplementary material, which is available to authorized users.

Genetic diversity of the msp-1, msp-2, and glurp genes of Plasmodium falciparum isolates in Northwest Ethiopia

Background

Determination of the genetic diversity of malaria parasites can inform the intensity of transmission and identify potential deficiencies in malaria control programmes. This study was conducted to characterize the genetic diversity and allele frequencies of Plasmodium falciparum in Northwest Ethiopia along the Eritrea and Sudan border.

Methods

A total of 90 isolates from patients presenting to the local health centre with uncomplicated P. falciparum were collected from October 2014 to January 2015. DNA was extracted and the polymorphic regions of the msp-1, msp-2 and glurp loci were genotyped by nested polymerase chain reactions followed by gel electrophoresis for fragment analysis.

Results

Allelic variation in msp-1, msp-2 and glurp were identified in 90 blood samples. A total of 34 msp alleles (12 for msp-1 and 22 for msp-2) were detected. For msp-1 97.8% (88/90), msp-2 82.2% (74/90) and glurp 46.7% (42/90) were detected. In msp-1, MAD20 was the predominant allelic family detected in 47.7% (42/88) of the isolates followed by RO33 and K1. For msp-2, the frequency of FC27 and IC/3D7 were 77% (57/74) and 76% (56/74), respectively. Nine glurp RII region genotypes were identified. Seventy percent of isolates had multiple genotypes and the overall mean multiplicity of infection was 2.6 (95% CI 2.25–2.97). The heterozygosity index was 0.82, 0.62 and 0.20 for msp-1, msp-2 and glurp, respectively. There was no significant association between multiplicity of infection and age or parasite density.

Conclusions

There was a high degree of genetic diversity with multiple clones in P. falciparum isolates from Northwest Ethiopia suggesting that there is a need for improved malaria control efforts in this region.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2540-x) contains supplementary material, which is available to authorized users.

Rapid selection of sulphadoxine-resistant Plasmodium falciparum and its effect on within-population genetic diversity in Papua New Guinea

The ability of the human malarial parasite Plasmodium falciparum to adapt to environmental changes depends considerably on its ability to maintain within-population genetic variation. Strong selection, consequent to widespread antimalarial drug usage, occasionally elicits a rapid expansion of drug-resistant isolates, which can act as founders. To investigate whether this phenomenon induces a loss of within-population genetic variation, we performed a population genetic analysis on 302 P. falciparum cases detected during two cross-sectional surveys in 2002/2003, just after the official introduction of sulphadoxine/pyrimethamine as a first-line treatment, and again in 2010/2011, in highly endemic areas in Papua New Guinea. We found that a single-origin sulphadoxine-resistant parasite isolate rapidly increased from 0% in 2002/2003 to 54% in 2010 and 84% in 2011. However, a considerable number of pairs exhibited random associations among 10 neutral microsatellite markers located in various chromosomes, suggesting that outcrossing effectively reduced non-random associations, albeit at a low average multiplicity of infection (1.35–1.52). Within-population genetic diversity was maintained throughout the study period. This indicates that the parasites maintained within-population variation, even after a clonal expansion of drug-resistant parasites. Outcrossing played a role in the preservation of within-population genetic diversity despite low levels of multiplicity of infection.

Signatures of competition and strain structure within the major blood-stage antigen of Plasmodium falciparum in a local community in Ghana

The concept of niche partitioning has received considerable theoretical attention at the interface of ecology and evolution of infectious diseases. Strain theory postulates that pathogen populations can be structured into distinct nonoverlapping strains by frequency-dependent selection in response to intraspecific competition for host immune space. The malaria parasite Plasmodium falciparum presents an opportunity to investigate this phenomenon in nature, under conditions of high recombination rate and extensive antigenic diversity. The parasite's major blood-stage antigen, Pf EMP1, is encoded by the hyperdiverse var genes. With a dataset that includes thousands of var DBLα sequence types sampled from asymptomatic cases within an area of high endemicity in Ghana, we address how var diversity is distributed within isolates and compare this to the distribution of microsatellite allelic diversity within isolates to test whether antigenic and neutral regions of the genome are structured differently. With respect to var DBLα sequence types, we find that on average isolates exhibit significantly lower overlap than expected randomly, but that there also exists frequent pairs of isolates that are highly related. Furthermore, the linkage network of var DBLα sequence types reveals a pattern of nonrandom modularity unique to these antigenic genes, and we find that modules of highly linked DBLα types are not explainable by neutral forces related to var recombination constraints, microsatellite diversity, sampling location, host age, or multiplicity of infection. These findings of reduced overlap and modularity among the var antigenic genes are consistent with a role for immune selection as proposed by strain theory. Identifying the evolutionary and ecological dynamics that are responsible for the nonrandom structure in P. falciparum antigenic diversity is important for designing effective intervention in endemic areas.

Nationwide genetic surveillance of Plasmodium vivax in Papua New Guinea reveals heterogeneous transmission dynamics and routes of migration amongst subdivided populations

The Asia Pacific Leaders in Malaria Alliance (APLMA) have committed to eliminate malaria from the region by 2030. Papua New Guinea (PNG) has the highest malaria burden in the Asia-Pacific region but with the intensification of control efforts since 2005, transmission has been dramatically reduced and Plasmodium vivax is now the dominant malaria infection in some parts of the country. To gain a better understanding of the transmission dynamics and migration patterns of P. vivax in PNG, here we investigate population structure in eight geographically and ecologically distinct regions of the country. A total of 219 P. vivax isolates (16-30 per population) were successfully haplotyped using 10 microsatellite markers. A wide range of genetic diversity (H_e=0.37-0.87, R_s=3.60-7.58) and significant multilocus linkage disequilibrium (LD) was observed in six of the eight populations (I_A^S=0.08-0.15 p-value<0.05) reflecting a spectrum of transmission intensities across the country. Genetic differentiation between regions was evident (Jost's D=0.07-0.72), with increasing divergence of populations with geographic distance. Overall, P. vivax isolates clustered into three major genetic populations subdividing the Mainland lowland and coastal regions, the Islands and the Highlands. P. vivax gene flow follows major human migration routes, and there was higher gene flow amongst Mainland parasite populations than among Island populations. The Central Province (samples collected in villages close to the capital city, Port Moresby), acts as a sink for imported infections from the three major endemic areas. These insights into P. vivax transmission dynamics and population networks will inform targeted strategies to contain malaria infections and to prevent the spread of drug resistance in PNG.

Cost-effectiveness of artemisinin-naphthoquine versus artemether-lumefantrine for the treatment of uncomplicated malaria in Papua New Guinean children

BACKGROUND

A recent randomized trial showed that artemisinin-naphthoquine (AN) was non-inferior to artemether-lumefantrine (AL) for falciparum malaria and superior for vivax malaria in young Papua New Guinean children. The aim of this study was to compare the cost-effectiveness of these two regimens.

METHODS

An incremental cost-effectiveness analysis was performed using data from 231 children with Plasmodium falciparum and/or Plasmodium vivax infections in an open-label, randomized, parallel-group trial. Recruited children were randomized 1:1 to receive once daily AN for 3 days with water or twice daily AL for 3 days given with fat. World Health Organisation (WHO) definitions were used to determine clinical/parasitological outcomes. The cost of transport between the home and clinic, plus direct health-care costs, served as a basis for determining each regimen's incremental cost per incremental treatment success relative to AL by Day 42 and its cost per life year saved.

RESULTS

In the usual care setting, AN was more effective for the treatment of uncomplicated malaria in children aged 0.5-5.9 years. AL and AN were equally efficacious for the treatment of falciparum malaria, however AN had increased anti-malarial treatment costs per patient of $10.46, compared with AL. AN was the most effective regimen for treatment of vivax malaria, but had increased treatment costs of $14.83 per treatment success compared with AL.

CONCLUSIONS

Whilst AN has superior overall efficacy for the treatment of uncomplicated malaria in PNG children, AL was the less costly regimen. An indicative extrapolation estimated the cost per life year saved by using AN instead of AL to treat uncomplicated malaria to be $12,165 for girls and $12,469 for boys (discounted), which means AN may not be cost-effective and affordable for PNG at current cost. However, AN may become acceptable should it become WHO prequalified and/or should donated/subsidized drug supply become available.

Oduro A, Koram KA, Day KP. Lack of Geospatial Population Structure Yet Significant Linkage Disequilibrium in the Reservoir of Plasmodium falciparum in Bongo District, Ghana

Malaria control in West Africa is impeded by the large reservoir of chronic asymptomatic Plasmodium falciparum infections in the human population. This study aimed to assess the extent of diversity in the P. falciparum reservoir in Bongo District (BD), Ghana, at the end of the dry season, the lowest point in malaria transmission over the course of the year. Analysis of the variation in 12 microsatellite loci was completed for 200 P. falciparum isolates collected from a cross-sectional survey of residents of all ages from two catchment areas in BD. Analysis of the multilocus haplotypes showed high levels of genetic diversity (He = 0.74), no population differentiation yet significant linkage disequilibrium (LD) (ISA = 0.0127, P = 0.006) in BD. Multilocus LD was significant between and within catchment areas even though every haplotype in the population was unique and the majority of individuals (84.0%) harbored multiple-clone infections. The linkage structure among multilocus haplotypes was not associated with sampling location. These data provide the first study with deep sampling of the P. falciparum reservoir in an area of seasonal malaria transmission in West Africa. The co-occurrence of high multiplicity of infection (multiple-clone infections) with significant multilocus LD is surprising given the likelihood of high recombination rates in BD. The results suggest that the linkage structure among multilocus haplotypes has not been shaped by geographic separation of parasite populations. Furthermore, the observed LD levels provide a baseline population genetic metric with putatively neutral markers to evaluate the effects of seasonality and malaria control efforts in BD.

Hemoglobin Levels and the Risk of Malaria in Papua New Guinean Infants: A Nested Cohort Study

Studies are available that assess the risk of malaria in accordance to the body's iron store and the systematic iron supplementation of preschool children. However, only a few studies evaluated the temporal association between hemoglobin and malaria and their results are opposing. A total of 1,650 3-month-old Papua New Guinean infants were enrolled in this study and followed-up for 12 months. The risk of malaria was assessed in all children every 3 months and with each episode of fever. The incidence of clinical malaria between 3 and 15 months of age was 249 cases per 1,000 infants per year. After adjustment for potential confounding factors, a decrease of 1 g/dL of hemoglobin was associated with a nonsignificant increase of 11% for risk of malaria infection (hazard ratio, 1.11, 95% confidence interval; CI, 0.99-1.25, P = 0.076). Only children with severe anemia (hemoglobin < 8.0 g/dL) at baseline were at higher risk of malaria infection (hazard ratio, 1.72, 95% CI, 1.08-2.76, P = 0.023) during the follow-up year compared with the control group (Hemoglobin > 10.0 g/dL). This association was not statistically significant if only clinical malaria episodes were taken into account (hazard ratio, 1.42, 95% CI, 0.77-2.61, P = 0.26). Our study suggests that infants with lower hemoglobin levels are not protected against malaria infection. Further research that examines the risk of malaria in relation to both hemoglobin and iron store levels would be important to better understand this complex interaction.

Transmission dynamics of co-endemic Plasmodium vivax and P. falciparum in Ethiopia and prevalence of antimalarial resistant genotypes

Ethiopia is one of the few African countries where Plasmodium vivax is co-endemic with P. falciparum. Malaria transmission is seasonal and transmission intensity varies mainly by landscape and climate. Although the recent emergence of drug resistant parasites presents a major issue to malaria control in Ethiopia, little is known about the transmission pathways of parasite species and prevalence of resistant markers. This study used microsatellites to determine population diversity and gene flow patterns of P. falciparum (N = 226) and P. vivax (N = 205), as well as prevalence of drug resistant markers to infer the impact of gene flow and existing malaria treatment regimes. Plasmodium falciparum indicated a higher rate of polyclonal infections than P. vivax. Both species revealed moderate genetic diversity and similar population structure. Populations in the northern highlands were closely related to the eastern Rift Valley, but slightly distinct from the southern basin area. Gene flow via human migrations between the northern and eastern populations were frequent and mostly bidirectional. Landscape genetic analyses indicated that environmental heterogeneity and geographical distance did not constrain parasite gene flow. This may partly explain similar patterns of resistant marker prevalence. In P. falciparum, a high prevalence of mutant alleles was detected in codons related to chloroquine (pfcrt and pfmdr1) and sulfadoxine-pyrimethamine (pfdhps and pfdhfr) resistance. Over 60% of the samples showed pfmdr1 duplications. Nevertheless, no mutation was detected in pfK13 that relates to artemisinin resistance. In P. vivax, while sequences of pvcrt-o were highly conserved and less than 5% of the samples showed pvmdr duplications, over 50% of the samples had pvmdr1 976F mutation. It remains to be tested if this mutation relates to chloroquine resistance. Monitoring the extent of malaria spread and markers of drug resistance is imperative to inform policy for evidence-based antimalarial choice and interventions. To effectively reduce malaria burden in Ethiopia, control efforts should focus on seasonal migrant populations.

Higher Complexity of Infection and Genetic Diversity of Plasmodium vivax Than Plasmodium falciparum Across All Malaria Transmission Zones of Papua New Guinea

Plasmodium falciparum and Plasmodium vivax have varying transmission dynamics that are informed by molecular epidemiology. This study aimed to determine the complexity of infection and genetic diversity of P. vivax and P. falciparum throughout Papua New Guinea (PNG) to evaluate transmission dynamics across the country. In 2008-2009, a nationwide malaria indicator survey collected 8,936 samples from all 16 endemic provinces of PNG. Of these, 892 positive P. vivax samples were genotyped at PvMS16 and PvmspF3, and 758 positive P. falciparum samples were genotyped at Pfmsp2. The data were analyzed for multiplicity of infection (MOI) and genetic diversity. Overall, P. vivax had higher polyclonality (71%) and mean MOI (2.32) than P. falciparum (20%, 1.39). These measures were significantly associated with prevalence for P. falciparum but not for P. vivax. The genetic diversity of P. vivax (PvMS16: expected heterozygosity = 0.95, 0.85-0.98; PvMsp1F3: 0.78, 0.66-0.89) was higher and less variable than that of P. falciparum (Pfmsp2: 0.89, 0.65-0.97). Significant associations of MOI with allelic richness (rho = 0.69, P = 0.009) and expected heterozygosity (rho = 0.87, P < 0.001) were observed for P. falciparum. Conversely, genetic diversity was not correlated with polyclonality nor mean MOI for P. vivax. The results demonstrate higher complexity of infection and genetic diversity of P. vivax across the country. Although P. falciparum shows a strong association of these parameters with prevalence, a lack of association was observed for P. vivax and is consistent with higher potential for outcrossing of this species.

Genetic clustering and polymorphism of the merozoite surface protein-3 of Plasmodium knowlesi clinical isolates from Peninsular Malaysia

Background

The simian malaria parasite Plasmodium knowlesi has been reported to cause significant numbers of human infection in South East Asia. Its merozoite surface protein-3 (MSP3) is a protein that belongs to a multi-gene family of proteins first found in Plasmodium falciparum. Several studies have evaluated the potential of P. falciparum MSP3 as a potential vaccine candidate. However, to date no detailed studies have been carried out on P. knowlesi MSP3 gene (pkmsp3). The present study investigates the genetic diversity, and haplotypes groups of pkmsp3 in P. knowlesi clinical samples from Peninsular Malaysia.

Methods

Blood samples were collected from P. knowlesi malaria patients within a period of 4 years (2008–2012). The pkmsp3 gene of the isolates was amplified via PCR, and subsequently cloned and sequenced. The full length pkmsp3 sequence was divided into Domain A and Domain B. Natural selection, genetic diversity, and haplotypes of pkmsp3 were analysed using MEGA6 and DnaSP ver. 5.10.00 programmes.

Results

From 23 samples, 48 pkmsp3 sequences were successfully obtained. At the nucleotide level, 101 synonymous and 238 non-synonymous mutations were observed. Tests of neutrality were not significant for the full length, Domain A or Domain B sequences. However, the dN/dS ratio of Domain B indicates purifying selection for this domain. Analysis of the deduced amino acid sequences revealed 42 different haplotypes. Neighbour Joining phylogenetic tree and haplotype network analyses revealed that the haplotypes clustered into two distinct groups.

Conclusions

A moderate level of genetic diversity was observed in the pkmsp3 and only the C-terminal region (Domain B) appeared to be under purifying selection. The separation of the pkmsp3 into two haplotype groups provides further evidence of the existence of two distinct P. knowlesi types or lineages. Future studies should investigate the diversity of pkmsp3 among P. knowlesi isolates in North Borneo, where large numbers of human knowlesi malaria infection still occur.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1935-1) contains supplementary material, which is available to authorized users.

Dissecting malaria biology and epidemiology using population genetics and genomics

Molecular approaches have an increasingly recognized utility in surveillance of malaria parasite populations, not only in defining prevalence and incidence with higher sensitivity than traditional methods, but also in monitoring local and regional parasite transmission patterns. In this review, we provide an overview of population genetic and genomic studies of human-infecting Plasmodium species, highlighting recent advances in the field. In accordance with the renewed impetus for malaria eradication, many studies are now using genetic and genomic epidemiology to support local evidence-based intervention strategies. Microsatellite genotyping remains a popular approach for both Plasmodium falciparum and Plasmodium vivax. However, with the increasing availability of whole genome sequencing data enabling effective single nucleotide polymorphism-based panels tailored to a given study question and setting, this approach is gaining popularity. The availability of new reference genomes for Plasmodium malariae and Plasmodium ovale should see a surge in similar molecular studies on these currently neglected species. Genomic studies are revealing new insights into important adaptive mechanisms of the parasite including antimalarial drug resistance. The advent of new methodologies such as selective whole genome amplification for dealing with extensive human DNA in low density field isolates should see genome-wide approaches becoming routine for parasite surveillance once the economic costs outweigh the current cost benefits of targeted approaches.

Genetic Diversity and Population Structure of Plasmodium falciparum in Lake Victoria Islands, A Region of Intense Transmission

Understanding the genetic structure and transmission dynamics of Plasmodium falciparum parasites in malaria-endemic regions is crucial before the implementation of interventions. Located in a high-transmission region of western Kenya where P. falciparum is the predominant species, the Lake Victoria islands are ideal for feasibility of malaria elimination studies. We analyzed genetic variation in eight microsatellite loci to examine parasite population structure and gene flow patterns across five sites. High levels of genetic diversity were measured throughout the region (mean heterozygosity index = 0.84). The overall fixation index value between the sites was 0.044, indicating that approximately 5% of the overall allelic variation is due to differences between the populations. Based on these results, we concluded that parasite population structure in the studied islands is shaped by human migration patterns that maintain extensive parasite gene flow between the sites. Consequently, any malaria elimination and interventions strategies in the study area will have to be carried out broadly on all four islands and adjoining mainland region.

Genetically Determined Response to Artemisinin Treatment in Western Kenyan Plasmodium falciparum Parasites

Genetically determined artemisinin resistance in Plasmodium falciparum has been described in Southeast Asia. The relevance of recently described Kelch 13-propeller mutations for artemisinin resistance in Sub-Saharan Africa parasites is still unknown. Southeast Asia parasites have low genetic diversity compared to Sub-Saharan Africa, where parasites are highly genetically diverse. This study attempted to elucidate whether genetics provides a basis for discovering molecular markers in response to artemisinin drug treatment in P. falciparum in Kenya. The genetic diversity of parasites collected pre- and post- introduction of artemisinin combination therapy (ACT) in western Kenya was determined. A panel of 12 microsatellites and 91 single nucleotide polymorphisms (SNPs) distributed across the P. falciparum genome were genotyped. Parasite clearance rates were obtained for the post-ACT parasites. The 12 microsatellites were highly polymorphic with post-ACT parasites being significantly more diverse compared to pre-ACT (p < 0.0001). The median clearance half-life was 2.55 hours for the post-ACT parasites. Based on SNP analysis, 15 of 90 post-ACT parasites were single-clone infections. Analysis revealed 3 SNPs that might have some causal association with parasite clearance rates. Further, genetic analysis using Bayesian tree revealed parasites with similar clearance phenotypes were more closely genetically related. With further studies, SNPs described here and genetically determined response to artemisinin treatment might be useful in tracking artemisinin resistance in Kenya.

Genetic diversity of Plasmodium falciparum in human malaria cases in Mali

Background

In Mali, Plasmodium falciparum malaria is highly endemic and remains stable despite the implementation of various malaria control measures. Understanding P. falciparum population structure variations across the country could provide new insights to guide malaria control programmes. In this study, P. falciparum genetic diversity and population structure in regions of varying patterns of malaria transmission in Mali were analysed.

Methods

A total of 648 blood isolates adsorbed onto filter papers during population surveillance surveys (December 2012–March 2013, October 2013) in four distinct sites of Mali were screened for the presence of P. falciparum via quantitative PCR (qPCR). Multiple loci variable number of tandem repeats analysis (MLVA) using eight microsatellite markers was then performed on positive qPCR samples. Complete genotypes were then analysed for genetic diversity, genetic differentiation and linkage disequilibrium.

Results

Of 156 qPCR-positive samples, complete genotyping of 112 samples was achieved. The parasite populations displayed high genetic diversity (mean He = 0.77), which was consistent with a high level of malaria transmission in Mali. Genetic differentiation was low (F_ST < 0.02), even between sites located approximately 900 km apart, thereby illustrating marked gene flux amongst parasite populations. The lack of linkage disequilibrium further revealed an absence of local clonal expansion, which was corroborated by the genotype relationship results. In contrast to the stable genetic diversity level observed throughout the country, mean multiplicity of infection increased from north to south (from 1.4 to 2.06) and paralleled malaria transmission levels observed locally.

Conclusions

In Mali, the high level of genetic diversity and the pronounced gene flux amongst P. falciparum populations may represent an obstacle to control malaria. Indeed, results suggest that parasite populations are polymorphic enough to adapt to their host and to counteract interventions, such as anti-malarial vaccination. Additionally, the panmictic parasite population structure imply that resistance traits may disseminate freely from one area to another, making control measures performed at a local level ineffective.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1397-0) contains supplementary material, which is available to authorized users.

Microsatellite genotyping and genome-wide single nucleotide polymorphism-based indices of Plasmodium falciparum diversity within clinical infections

Background

In regions where malaria is endemic, individuals are often infected with multiple distinct parasite genotypes, a situation that may impact on evolution of parasite virulence and drug resistance. Most approaches to studying genotypic diversity have involved analysis of a modest number of polymorphic loci, although whole genome sequencing enables a broader characterisation of samples.

Methods

PCR-based microsatellite typing of a panel of ten loci was performed on Plasmodium falciparum in 95 clinical isolates from a highly endemic area in the Republic of Guinea, to characterize within-isolate genetic diversity. Separately, single nucleotide polymorphism (SNP) data from genome-wide short-read sequences of the same samples were used to derive within-isolate fixation indices (F_ws), an inverse measure of diversity within each isolate compared to overall local genetic diversity. The latter indices were compared with the microsatellite results, and also with indices derived by randomly sampling modest numbers of SNPs.

Results

As expected, the number of microsatellite loci with more than one allele in each isolate was highly significantly inversely correlated with the genome-wide F_ws fixation index (r = −0.88, P < 0.001). However, the microsatellite analysis revealed that most isolates contained mixed genotypes, even those that had no detectable genome sequence heterogeneity. Random sampling of different numbers of SNPs showed that an F_ws index derived from ten or more SNPs with minor allele frequencies of >10 % had high correlation (r > 0.90) with the index derived using all SNPs.

Conclusions

Different types of data give highly correlated indices of within-infection diversity, although PCR-based analysis detects low-level minority genotypes not apparent in bulk sequence analysis. When whole-genome data are not obtainable, quantitative assay of ten or more SNPs can yield a reasonably accurate estimate of the within-infection fixation index (F_ws).

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1324-4) contains supplementary material, which is available to authorized users.

Genetic Diversity of Plasmodium falciparum Populations in Malaria Declining Areas of Sabah, East Malaysia

Malaysia has a national goal to eliminate malaria by 2020. Understanding the genetic diversity of malaria parasites in residual transmission foci can provide invaluable information which may inform the intervention strategies used to reach elimination targets. This study was conducted to determine the genetic diversity level of P. falciparum isolates in malaria residual foci areas of Sabah. Malaria active case detection was conducted in Kalabakan and Kota Marudu. All individuals in the study sites were screened for malaria infection by rapid diagnostic test. Blood from P. falciparum-infected individuals were collected on filter paper prior to DNA extraction. Genotyping was performed using merozoite surface protein-1 (MSP-1), merozoite surface protein-2 (MSP-2), glutamate rich protein (GLURP) and 10 neutral microsatellite loci markers. The size of alleles, multiplicity of infection (MOI), mean number of alleles (Na), expected heterozygosity (He), linkage disequilibrium (LD) and genetic differentiation (FST) were determined. In Kalabakan, the MSP-1 and MSP-2 alleles were predominantly K1 and FC27 family types, respectively. The GLURP genotype VI (751-800 bp) was predominant. The MOI for MSP-1 and MSP-2 were 1.65 and 1.20, respectively. The Na per microsatellite locus was 1.70. The He values for MSP-1, MSP-2, GLURP and neutral microsatellites were 0.17, 0.37, 0.70 and 0.33, respectively. In Kota Marudu, the MSP-1 and MSP-2 alleles were predominantly MAD20 and 3D7 family types, respectively. The GLURP genotype IV (651-700 bp) was predominant. The MOI for both MSP-1 and MSP-2 was 1.05. The Na per microsatellite locus was 3.60. The He values for MSP-1, MSP-2, GLURP and neutral microsatellites were 0.24, 0.25, 0.69 and 0.30, respectively. A significant LD was observed in Kalabakan (0.495, p<0.01) and Kota Marudu P. falciparum populations (0.601, p<0.01). High genetic differentiation between Kalabakan and Kota Marudu P. falciparum populations was observed (FST = 0.532). The genetic data from the present study highlighted the limited diversity and contrasting genetic pattern of P. falciparum populations in the malaria declining areas of Sabah.

Genetic variability and population structure of Plasmodium falciparum parasite populations from different malaria ecological regions of Kenya

Transmission intensity, movement of human and vector hosts, biogeographical features, and malaria control measures are some of the important factors that determine Plasmodium falciparum parasite genetic variability and population structure. Kenya has different malaria ecologies which might require different disease intervention methods. Refined parasite population genetic studies are critical for informing malaria control and elimination strategies. This study describes the genetic diversity and population structure of P. falciparum parasites from the different malaria ecological zones in Kenya. Twelve multi-locus microsatellite (MS) loci previously described were genotyped in 225 P. falciparum isolates collected between 2012 and 2013 from five sites; three in lowland endemic regions (Kisumu, Kombewa, and Malindi) and two in highland, epidemic regions (Kisii and Kericho). Parasites from the lowland endemic and highland epidemic regions of western Kenya had high genetic diversity compared to coastal lowland endemic region of Kenya [Malindi]. The Kenyan parasites had a mean genetic differentiation index (FST) of 0.072 (p=0.011). The multi-locus genetic analysis of the 12 MS revealed all the parasites had unique haplotypes. Significant linkage disequilibrium (LD) was observed in all the five parasite populations. Kisumu had the most significant index of association values (0.16; p<0.0001) whereas Kisii had the least significant index of association values (0.03; p<0.0001). Our data suggest high genetic diversity in Kenyan parasite population with the exception of parasite from Malindi where malaria has been on the decline. The presence of significant LD suggests that there is occurrence of inbreeding in the parasite population. Parasite populations from Kisii showed the strongest evidence for epidemic population structure whereas the rest of the regions showed panmixia. Defining the genetic diversity of the parasites in different ecological regions of Kenya after introduction of the artemether-lumefantrine is important in refining the spread of drug resistant strains and malaria transmission for more effective control and eventual elimination of malaria in Kenya.

Accuracy of an HRP-2/panLDH rapid diagnostic test to detect peripheral and placental Plasmodium falciparum infection in Papua New Guinean women with anaemia or suspected malaria

Background

The diagnosis of malaria during pregnancy is complicated by placental sequestration, asymptomatic infection, and low-density peripheral parasitaemia. Where intermittent preventive treatment (IPT) with sulfadoxine-pyrimethamine is threatened by drug resistance, or is inappropriate due to low transmission, intermittent screening and treatment (ISTp) with rapid diagnostic tests for malaria (RDT) could be a valuable alternative. Therefore, the accuracy of RDTs to detect peripheral and placental infection was assessed in a declining transmission setting in Papua New Guinea (PNG).

Methods

The performance of a combination RDT detecting histidine-rich protein-2 (HRP-2) and Plasmodium lactate dehydrogenase (pLDH), and light microscopy (LM), to diagnose peripheral Plasmodium falciparum and Plasmodium vivax infections during pregnancy, were assessed using quantitative real-time PCR (qPCR) as the reference standard. Participants in a malaria prevention trial in PNG with a haemoglobin ≤90 g/L, or symptoms suggestive of malaria, were tested. Ability of RDT and LM to detect active placental infection on histology was evaluated in some participants.

Results

Among 876 women, 1162 RDTs were undertaken (anaemia: 854 [73.5 %], suspected malaria: 308 [26.5 %]). qPCR detected peripheral infection during 190 RDT episodes (165 P. falciparum, 19 P. vivax, 6 mixed infections). Overall, RDT detected peripheral P. falciparum infection with 45.6 % sensitivity (95 % CI 38.0–53.4), a specificity of 96.4 % (95.0–97.4), a positive predictive value of 68.4 % (59.1–76.8), and a negative predictive value of 91.1 % (89.2–92.8). RDT performance to detect P. falciparum was inferior to LM, more so amongst anaemic women (18.6 vs 45.3 % sensitivity, Liddell’s exact test, P < 0.001) compared to symptomatic women (72.9 vs 82.4 % sensitivity, P = 0.077). RDT and LM missed 88.0 % (22/25) and 76.0 % (19/25) of P. vivax infections, respectively. In a subset of women tested at delivery and who had placental histology (n = 158) active placental infection was present in 19.6 %: all three peripheral blood infection detection methods (RDT, LM, qPCR) missed >50 % of these infections.

Conclusions

In PNG, HRP-2/pLDH RDTs may be useful to diagnose peripheral P. falciparum infections in symptomatic pregnant women. However, they are not sufficiently sensitive for use in intermittent screening amongst asymptomatic (anaemic) women. These findings have implications for the management of malaria in pregnancy. The adverse impact of infections undetected by RDT or LM on pregnancy outcomes needs further evaluation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0927-5) contains supplementary material, which is available to authorized users.

Variation in Complexity of Infection and Transmission Stability between Neighbouring Populations of Plasmodium vivax in Southern Ethiopia

Background

P. vivax is an important public health burden in Ethiopia, accounting for almost half of all malaria cases. Owing to heterogeneous transmission across the country, a stronger evidence base on local transmission dynamics is needed to optimise allocation of resources and improve malaria interventions.

Methodology and Principal Findings

In a pilot evaluation of local level P. vivax molecular surveillance in southern Ethiopia, the diversity and population structure of isolates collected between May and November 2013 were investigated. Blood samples were collected from microscopy positive P. vivax patients recruited to clinical and cross-sectional surveys from four sites: Arbaminch, Halaba, Badawacho and Hawassa. Parasite genotyping was undertaken at nine tandem repeat markers. Eight loci were successfully genotyped in 197 samples (between 36 and 59 per site). Heterogeneity was observed in parasite diversity and structure amongst the sites. Badawacho displayed evidence of unstable transmission, with clusters of identical clonal infections. Linkage disequilibrium in Badawacho was higher (I_AS = 0.32, P = 0.010) than in the other populations (I_AS range = 0.01–0.02) and declined markedly after adjusting for identical infections (I_AS = 0.06, P = 0.010). Other than Badawacho (H_E = 0.70), population diversity was equivalently high across the sites (H_E = 0.83). Polyclonal infections were more frequent in Hawassa (67%) than the other populations (range: 8–44%). Despite the variable diversity, differentiation between the sites was low (F_ST range: 5 x 10⁻³–0.03).

Conclusions

Marked variation in parasite population structure likely reflects differing local transmission dynamics. Parasite genotyping in these heterogeneous settings has potential to provide important complementary information with which to optimise malaria control interventions.

Genetic Diversity of Plasmodium falciparum in Haiti: Insights from Microsatellite Markers

Hispaniola, comprising Haiti and the Dominican Republic, has been identified as a candidate for malaria elimination. However, incomplete surveillance data in Haiti hamper efforts to assess the impact of ongoing malaria control interventions. Characteristics of the genetic diversity of Plasmodium falciparum populations can be used to assess parasite transmission, which is information vital to evaluating malaria elimination efforts. Here we characterize the genetic diversity of P. falciparum samples collected from patients at seven sites in Haiti using 12 microsatellite markers previously employed in population genetic analyses of global P. falciparum populations. We measured multiplicity of infections, level of genetic diversity, degree of population geographic substructure, and linkage disequilibrium (defined as non-random association of alleles from different loci). For low transmission populations like Haiti, we expect to see few multiple infections, low levels of genetic diversity, high degree of population structure, and high linkage disequilibrium. In Haiti, we found low levels of multiple infections (12.9%), moderate to high levels of genetic diversity (mean number of alleles per locus = 4.9, heterozygosity = 0.61), low levels of population structure (highest pairwise F_st = 0.09 and no clustering in principal components analysis), and moderate linkage disequilibrium (ISA = 0.05, P<0.0001). In addition, population bottleneck analysis revealed no evidence for a reduction in the P. falciparum population size in Haiti. We conclude that the high level of genetic diversity and lack of evidence for a population bottleneck may suggest that Haiti’s P. falciparum population has been stable and discuss the implications of our results for understanding the impact of malaria control interventions. We also discuss the relevance of parasite population history and other host and vector factors when assessing transmission intensity from genetic diversity data.

Strategies for understanding and reducing the Plasmodium vivax and Plasmodium ovale hypnozoite reservoir in Papua New Guinean children: a randomised placebo-controlled trial and mathematical model

BACKGROUND

The undetectable hypnozoite reservoir for relapsing Plasmodium vivax and P. ovale malarias presents a major challenge for malaria control and elimination in endemic countries. This study aims to directly determine the contribution of relapses to the burden of P. vivax and P. ovale infection, illness, and transmission in Papua New Guinean children.

METHODS AND FINDINGS

From 17 August 2009 to 20 May 2010, 524 children aged 5-10 y from East Sepik Province in Papua New Guinea (PNG) participated in a randomised double-blind placebo-controlled trial of blood- plus liver-stage drugs (chloroquine [CQ], 3 d; artemether-lumefantrine [AL], 3 d; and primaquine [PQ], 20 d, 10 mg/kg total dose) (261 children) or blood-stage drugs only (CQ, 3 d; AL, 3 d; and placebo [PL], 20 d) (263 children). Participants, study staff, and investigators were blinded to the treatment allocation. Twenty children were excluded during the treatment phase (PQ arm: 14, PL arm: 6), and 504 were followed actively for 9 mo. During the follow-up time, 18 children (PQ arm: 7, PL arm: 11) were lost to follow-up. Main primary and secondary outcome measures were time to first P. vivax infection (by qPCR), time to first clinical episode, force of infection, gametocyte positivity, and time to first P. ovale infection (by PCR). A basic stochastic transmission model was developed to estimate the potential effect of mass drug administration (MDA) for the prevention of recurrent P. vivax infections. Targeting hypnozoites through PQ treatment reduced the risk of having at least one qPCR-detectable P. vivax or P. ovale infection during 8 mo of follow-up (P. vivax: PQ arm 0.63/y versus PL arm 2.62/y, HR = 0.18 [95% CI 0.14, 0.25], p < 0.001; P. ovale: 0.06 versus 0.14, HR = 0.31 [95% CI 0.13, 0.77], p = 0.011) and the risk of having at least one clinical P. vivax episode (HR = 0.25 [95% CI 0.11, 0.61], p = 0.002). PQ also reduced the molecular force of P. vivax blood-stage infection in the first 3 mo of follow-up (PQ arm 1.90/y versus PL arm 7.75/y, incidence rate ratio [IRR] = 0.21 [95% CI 0.15, 0.28], p < 0.001). Children who received PQ were less likely to carry P. vivax gametocytes (IRR = 0.27 [95% CI 0.19, 0.38], p < 0.001). PQ had a comparable effect irrespective of the presence of P. vivax blood-stage infection at the time of treatment (p = 0.14). Modelling revealed that mass screening and treatment with highly sensitive quantitative real-time PCR, or MDA with blood-stage treatment alone, would have only a transient effect on P. vivax transmission levels, while MDA that includes liver-stage treatment is predicted to be a highly effective strategy for P. vivax elimination. The inclusion of a directly observed 20-d treatment regime maximises the efficiency of hypnozoite clearance but limits the generalisability of results to real-world MDA programmes.

CONCLUSIONS

These results suggest that relapses cause approximately four of every five P. vivax infections and at least three of every five P. ovale infections in PNG children and are important in sustaining transmission. MDA campaigns combining blood- and liver-stage treatment are predicted to be a highly efficacious intervention for reducing P. vivax and P. ovale transmission.

TRIAL REGISTRATION

ClinicalTrials.gov NCT02143934.

Malaria Molecular Epidemiology: Lessons from the International Centers of Excellence for Malaria Research Network

Molecular epidemiology leverages genetic information to study the risk factors that affect the frequency and distribution of malaria cases. This article describes molecular epidemiologic investigations currently being carried out by the International Centers of Excellence for Malaria Research (ICEMR) network in a variety of malaria-endemic settings. First, we discuss various novel approaches to understand malaria incidence and gametocytemia, focusing on Plasmodium falciparum and Plasmodium vivax. Second, we describe and compare different parasite genotyping methods commonly used in malaria epidemiology and population genetics. Finally, we discuss potential applications of molecular epidemiological tools and methods toward malaria control and elimination efforts.

Evaluation of the population structure and genetic diversity of Plasmodium falciparum in southern China

Background

Yunnan and Hainan provinces are the two major endemic regions for Plasmodiumfalciparum malaria in China. However, few studies have investigated the characteristics of this parasite. Therefore, this study aimed to evaluate the genetic diversity and population structure of P. falciparum to predict the geographic origin of falciparum malaria.

Methods

Thirteen highly polymorphic microsatellite loci were studied to estimate the genetic diversity and population structure of 425 P. falciparum isolates obtained from blood samples collected from Yunnan and Hainan provinces of South China. The isolates were analysed for genetic diversity, linkage disequilibrium, and population structure. The parasite populations were clustered into two subgroups (i.e., Yunnan and Hainan) and a classification algorithm was used to identify molecular markers for classifying the P. falciparum populations.

Results

All 13 microsatellite loci were highly polymorphic, with the number of alleles per locus varying from 5 to 20. The mean expected heterozygosity (He) in Yunnan and Hainan was 0.766 ± 0.036 and 0.677 ± 0.039, respectively, revealing a moderate high level of genetic diversity. Significant linkage disequilibrium was found for some regions of Yunnan (Lazan county and Xishuangbanna region) and Hainan (Dongfang city and Sanya city) province. According to the classification algorithm, a combination of three microsatellites could be used as a discriminatory marker to identify the origin of P. falciparum isolates.

Conclusions

The results on the genetic structure of P. falciparum populations from South China provide a basis for developing a genetic marker-based tool to trace the source of the parasite infections and consequently improve malaria control and elimination strategies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0786-0) contains supplementary material, which is available to authorized users.

Mitochondrial population genomic analyses reveal population structure and demography of Indian Plasmodium falciparum

Inference on the genetic diversity of Plasmodium falciparum populations could help in better management of malaria. A very recent study with mitochondrial (mt) genomes in global P. falciparum had revealed interesting evolutionary genetic patterns of Indian isolates in comparison to global ones. However, no population genetic study using the whole mt genome sequences of P. falciparum isolates collected in the entire distribution range in India has yet been performed. We herewith have analyzed 85 whole mt genomes (48 already published and 37 entirely new) sampled from eight differentially endemic Indian locations to estimate genetic diversity and infer population structure and historical demography of Indian P. falciparum. We found 19 novel Indian-specific Single Nucleotide Polymorphisms (SNPs) and 22 novel haplotypes segregating in Indian P. falciparum. Accordingly, high haplotype and nucleotide diversities were detected in Indian P. falciparum in comparison to many other global isolates. Indian P. falciparum populations were found to be moderately sub-structured with four different genetic clusters. Interestingly, group of local populations aggregate to form each cluster; while samples from Jharkhand and Odisha formed a single cluster, P. falciparum isolates from Asom formed an independent one. Similarly, Surat, Bilaspur and Betul formed a single cluster and Goa and Mangalore formed another. Interestingly, P. falciparum isolates from the two later populations were significantly genetically differentiated from isolates collected in other six Indian locations. Signature of historical population expansion was evident in five population samples, and the onset of expansion event was found to be very similar to African P. falciparum. In agreement with the previous finding, the estimated Time to Most Recent Common Ancestor (TMRCA) and the effective population size were high in Indian P. falciparum. All these genetic features of Indian P. falciparum with high mt genome diversity are somehow similar to Africa, but quite different from other Asian population samples.

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/μL) and it remains to be determined to what extent low-density gametocyte positive individuals are infective to mosquitos.

Contrasting Transmission Dynamics of Co-endemic Plasmodium vivax and P. falciparum: Implications for Malaria Control and Elimination

BACKGROUND

Outside of Africa, P. falciparum and P. vivax usually coexist. In such co-endemic regions, successful malaria control programs have a greater impact on reducing falciparum malaria, resulting in P. vivax becoming the predominant species of infection. Adding to the challenges of elimination, the dormant liver stage complicates efforts to monitor the impact of ongoing interventions against P. vivax. We investigated molecular approaches to inform the respective transmission dynamics of P. falciparum and P. vivax and how these could help to prioritize public health interventions.

METHODOLOGY/PRINCIPAL FINDINGS

Genotype data generated at 8 and 9 microsatellite loci were analysed in 168 P. falciparum and 166 P. vivax isolates, respectively, from four co-endemic sites in Indonesia (Bangka, Kalimantan, Sumba and West Timor). Measures of diversity, linkage disequilibrium (LD) and population structure were used to gauge the transmission dynamics of each species in each setting. Marked differences were observed in the diversity and population structure of P. vivax versus P. falciparum. In Bangka, Kalimantan and Timor, P. falciparum diversity was low, and LD patterns were consistent with unstable, epidemic transmission, amenable to targeted intervention. In contrast, P. vivax diversity was higher and transmission appeared more stable. Population differentiation was lower in P. vivax versus P. falciparum, suggesting that the hypnozoite reservoir might play an important role in sustaining local transmission and facilitating the spread of P. vivax infections in different endemic settings. P. vivax polyclonality varied with local endemicity, demonstrating potential utility in informing on transmission intensity in this species.

CONCLUSIONS/SIGNIFICANCE

Molecular approaches can provide important information on malaria transmission that is not readily available from traditional epidemiological measures. Elucidation of the transmission dynamics circulating in a given setting will have a major role in prioritising malaria control strategies, particularly against the relatively neglected non-falciparum species.

Uncovering the transmission dynamics of Plasmodium vivax using population genetics

Population genetic analysis of malaria parasites has the power to reveal key insights into malaria epidemiology and transmission dynamics with the potential to deliver tools to support control and elimination efforts. Analyses of parasite genetic diversity have suggested that Plasmodium vivax populations are more genetically diverse and less structured than those of Plasmodium falciparum indicating that P. vivax may be a more ancient parasite of humans and/or less susceptible to population bottlenecks, as well as more efficient at disseminating its genes. These population genetic insights into P. vivax transmission dynamics provide an explanation for its relative resilience to control efforts. Here, we describe current knowledge on P. vivax population genetic structure, its relevance to understanding transmission patterns and relapse and how this information can inform malaria control and elimination programmes.

Plasmodium vivax populations are more genetically diverse and less structured than sympatric Plasmodium falciparum populations

Introduction

The 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 (H_e = 0.83–0.85) was higher than that of P. falciparum (H_e = 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/Significance

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

The neglected human malaria parasite Plasmodium vivax is responsible for a large proportion of the global malaria burden. Efforts to control malaria have revealed that P. vivax is more resilient than the other major human malaria parasite, Plasmodium falciparum. This study utilised population genetics to compare patterns of P. vivax and P. falciparum transmission in Papua New Guinea, a region where infection rates of the two species are similar. The results demonstrated that P. vivax populations are more genetically diverse than those of P. falciparum suggestive of a parasite population that is more resilient to environmental challenges, undergoing higher levels of interbreeding locally and between distant parasite populations. Unique characteristics of P. vivax such as relapse, which allows different strains from past infections to produce subsequent infections, may provide more opportunities for the exchange and dissemination of genetic material. The contrasting patterns observed for the two species may be the result of a differential impact of past elimination attempts and indicate that more rigorous interventions will be needed in efforts to control and eventually eliminate P. vivax.

Sulphadoxine-pyrimethamine plus azithromycin for the prevention of low birthweight in Papua New Guinea: a randomised controlled trial

BACKGROUND

Intermittent preventive treatment in pregnancy has not been evaluated outside of Africa. Low birthweight (LBW, <2,500 g) is common in Papua New Guinea (PNG) and contributing factors include malaria and reproductive tract infections.

METHODS

From November 2009 to February 2013, we conducted a parallel group, randomised controlled trial in pregnant women (≤ 26 gestational weeks) in PNG. Sulphadoxine-pyrimethamine (1,500/75 mg) plus azithromycin (1 g twice daily for 2 days) (SPAZ) monthly from second trimester (intervention) was compared against sulphadoxine-pyrimethamine and chloroquine (450 to 600 mg, daily for three days) (SPCQ) given once, followed by SPCQ placebo (control). Women were assigned to treatment (1:1) using a randomisation sequence with block sizes of 32. Participants were blinded to assignments. The primary outcome was LBW. Analysis was by intention-to-treat.

RESULTS

Of 2,793 women randomised, 2,021 (72.4%) were included in the primary outcome analysis (SPCQ: 1,008; SPAZ: 1,013). The prevalence of LBW was 15.1% (305/2,021). SPAZ reduced LBW (risk ratio [RR]: 0.74, 95% CI: 0.60-0.91, P = 0.005; absolute risk reduction (ARR): 4.5%, 95% CI: 1.4-7.6; number needed to treat: 22), and preterm delivery (0.62, 95% CI: 0.43-0.89, P = 0.010), and increased mean birthweight (41.9 g, 95% CI: 0.2-83.6, P = 0.049). SPAZ reduced maternal parasitaemia (RR: 0.57, 95% CI: 0.35-0.95, P = 0.029) and active placental malaria (0.68, 95% CI: 0.47-0.98, P = 0.037), and reduced carriage of gonorrhoea (0.66, 95% CI: 0.44-0.99, P = 0.041) at second visit. There were no treatment-related serious adverse events (SAEs), and the number of SAEs (intervention 13.1% [181/1,378], control 12.7% [174/1,374], P = 0.712) and AEs (intervention 10.5% [144/1,378], control 10.8% [149/1,374], P = 0.737) was similar. A major limitation of the study was the high loss to follow-up for birthweight.

CONCLUSIONS

SPAZ was efficacious and safe in reducing LBW, possibly acting through multiple mechanisms including the effect on malaria and on sexually transmitted infections. The efficacy of SPAZ in the presence of resistant parasites and the contribution of AZ to bacterial antibiotic resistance require further study. The ability of SPAZ to improve pregnancy outcomes warrants further evaluation.

TRIAL REGISTRATION

ClinicalTrials.gov NCT01136850 (06 April 2010).

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.

Genetic Evidence of Importation of Drug-Resistant Plasmodium falciparum to Guatemala from the Democratic Republic of the Congo

Molecular markers and population genetics were effective tracking tools.

Imported malaria threatens control and elimination efforts in countries that have low rates of transmission. In 2010, an outbreak of Plasmodium falciparum malaria was reported among United Nations peacekeeping soldiers from Guatemala who had recently returned from the Democratic Republic of the Congo (DRC). Epidemiologic evidence suggested that the soldiers were infected in the DRC, but local transmission could not be ruled out in all cases. We used population genetic analyses of neutral microsatellites to determine the outbreak source. Genetic relatedness was compared among parasites found in samples from the soldiers and parasite populations collected in the DRC and Guatemala; parasites identified in the soldiers were more closely related to those from the DRC. A phylogenetic clustering analysis confirms this identification with >99.9% confidence. Thus, results support the hypothesis that the soldiers likely imported malaria from the DRC. This study demonstrates the utility of molecular genotyping in outbreak investigations.

Microsatellite markers reveal low levels of population sub-structuring of Plasmodium falciparum in southwestern Nigeria

BACKGROUND

Genetic diversity studies provide evidence of Plasmodium falciparum differentiation that could affect fitness and adaptation to drugs and target antigens for vaccine development. This study describes the genetic structure of P. falciparum populations in urban and rural sites from southwestern Nigeria.

METHODOLOGY

Ten neutral microsatellite loci were genotyped in 196 P. falciparum infections from three localities: Aramoko-Ekiti, a rural community; Lekki, an urban location and Badagry, a peri-urban border settlement. Analysis was performed on the genetic diversity, linkage disequilibrium, population structure and inter-population differentiation.

RESULTS

Allelic diversity values were similar across all populations, with mean expected heterozygosity (HE) values between 0.65 and 0.79. No matching multilocus haplotypes were found and analysis of multilocus LD showed no significant index of association. Genetic differentiation between populations was low (ΦPT = 0.017).

CONCLUSION

The absence of detectable population structure of P. falciparum in southwestern Nigeria is evident in the lack of significant differentiation between populations separated by about 200 km. This implies that a fairly uniform malaria control strategy may be effective over a wide geographic range in this highly endemic region. However, more wide-scale survey across the country will be required to inform malaria control in this large and densely populated endemic region.

The geography of malaria genetics in the Democratic Republic of Congo: A complex and fragmented landscape

Understanding how malaria parasites move between populations is important, particularly given the potential for malaria to be reintroduced into areas where it was previously eliminated. We examine the distribution of malaria genetics across seven sites within the Democratic Republic of Congo (DRC) and two nearby countries, Ghana and Kenya, in order to understand how the relatedness of malaria parasites varies across space, and whether there are barriers to the flow of malaria parasites within the DRC or across borders. Parasite DNA was retrieved from dried blood spots from 7 Demographic and Health Survey sample clusters in the DRC. Malaria genetic characteristics of parasites from Ghana and Kenya were also obtained. For each of 9 geographic sites (7 DRC, 1 Ghana and 1 Kenya), a pair-wise RST statistic was calculated, indicating the genetic distance between malaria parasites found in those locations. Mapping genetics across the spatial extent of the study area indicates a complex genetic landscape, where relatedness between two proximal sites may be relatively high (RST > 0.64) or low (RST < 0.05), and where distal sites also exhibit both high and low genetic similarity. Mantel's tests suggest that malaria genetics differ as geographic distances increase. Principal Coordinate Analysis suggests that genetically related samples are not co-located. Barrier analysis reveals no significant barriers to gene flow between locations. Malaria genetics in the DRC have a complex and fragmented landscape. Limited exchange of genes across space is reflected in greater genetic distance between malaria parasites isolated at greater geographic distances. There is, however, evidence for close genetic ties between distally located sample locations, indicating that movement of malaria parasites and flow of genes is being driven by factors other than distance decay. This research demonstrates the contributions that spatial disease ecology and landscape genetics can make to understanding the evolutionary dynamics of infectious diseases.

Controlled human malaria infection of Tanzanians by intradermal injection of aseptic, purified, cryopreserved Plasmodium falciparum sporozoites

Controlled human malaria infection (CHMI) by mosquito bite has been used to assess anti-malaria interventions in > 1,500 volunteers since development of methods for infecting mosquitoes by feeding on Plasmodium falciparum (Pf) gametocyte cultures. Such CHMIs have never been used in Africa. Aseptic, purified, cryopreserved Pf sporozoites, PfSPZ Challenge, were used to infect Dutch volunteers by intradermal injection. We conducted a double-blind, placebo-controlled trial to assess safety and infectivity of PfSPZ Challenge in adult male Tanzanians. Volunteers were injected intradermally with 10,000 (N = 12) or 25,000 (N = 12) PfSPZ or normal saline (N = 6). PfSPZ Challenge was well tolerated and safe. Eleven of 12 and 10 of 11 subjects, who received 10,000 and 25,000 PfSPZ respectively, developed parasitemia. In 10,000 versus 25,000 PfSPZ groups geometric mean days from injection to Pf positivity by thick blood film was 15.4 versus 13.5 (P = 0.023). Alpha-thalassemia heterozygosity had no apparent effect on infectivity. PfSPZ Challenge was safe, well tolerated, and infectious.