Alterations in Plasmodium falciparum genetic structure two years after increased malaria control efforts in western Kenya

Exploring how space, time, and sampling impact our ability to measure genetic structure across Plasmodium falciparum populations

A primary use of malaria parasite genomics is identifying highly related infections to quantify epidemiological, spatial, or temporal factors associated with patterns of transmission. For example, spatial clustering of highly related parasites can indicate foci of transmission and temporal differences in relatedness can serve as evidence for changes in transmission over time. However, for infections in settings of moderate to high endemicity, understanding patterns of relatedness is compromised by complex infections, overall high forces of infection, and a highly diverse parasite population. It is not clear how much these factors limit the utility of using genomic data to better understand transmission in these settings. In particular, further investigation is required to determine which patterns of relatedness we expect to see with high quality, densely sampled genomic data in a high transmission setting and how these observations change under different study designs, missingness, and biases in sample collection. Here we investigate two identity-by-state measures of relatedness and apply them to amplicon deep sequencing data collected as part of a longitudinal cohort in Western Kenya that has previously been analysed to identify individual-factors associated with sharing parasites with infected mosquitoes. With these data we use permutation tests, to evaluate several hypotheses about spatiotemporal patterns of relatedness compared to a null distribution. We observe evidence of temporal structure, but not of fine-scale spatial structure in the cohort data. To explore factors associated with the lack of spatial structure in these data, we construct a series of simplified simulation scenarios using an agent based model calibrated to entomological, epidemiological and genomic data from this cohort study to investigate whether the lack of spatial structure observed in the cohort could be due to inherent power limitations of this analytical method. We further investigate how our hypothesis testing behaves under different sampling schemes, levels of completely random and systematic missingness, and different transmission intensities.

Genomic miscellany and allelic frequencies of Plasmodium falciparum msp-1, msp-2 and glurp in parasite isolates

Introduction

The genomic miscellany of malaria parasites can help inform the intensity of malaria transmission and identify potential deficiencies in malaria control programs. This study was aimed at investigating the genomic miscellany, allele frequencies, and MOI of P. falciparum infection.

Methods

A total of 85 P. falciparum confirmed isolates out of 100 were included in this study that were collected from P. falciparum patients aged 4 months to 60 years in nine districts of Khyber Pakhtunkhwa Province. Parasite DNA was extracted from 200µL whole blood samples using the Qiagen DNA extraction kit following the manufacturer’s instructions. The polymorphic regions of msp-1, msp-2 and glurp loci were genotyped using nested PCR followed by gel electrophoresis for amplified fragments identification and subsequent data analysis.

Results

Out of 85 P. falciparum infections detected, 30 were msp-1 and 32 were msp-2 alleles specific. Successful amplification occurred in 88.23% (75/85) isolates for msp-1, 78.9% (67/85) for msp-2 and 70% (60/85) for glurp gene. In msp-1, the K1 allelic family was predominantly prevalent as 66.66% (50/75), followed by RO33 and MAD20. The frequency of samples with single infection having only K1, MAD20 and RO33 were 21.34% (16/75), 8% (6/75), and 10.67% (8/75), respectively. In msp-2, both the FC27 and 3D7 allelic families revealed almost the same frequencies as 70.14% (47/67) and 67.16% (45/67), respectively. Nine glurp RII region alleles were identified in 60 isolates. The overall mean multiplicity of infection for msp genes was 1.6 with 1.8 for msp-1 and 1.4 for msp-2, while for glurp the MOI was 1.03. There was no significant association between multiplicity of infection and age groups (Spearman’s rank coefficient = 0.050; P = 0.6) while MOI and parasite density correlated for only msp-2 allelic marker.

Conclusions

The study showed high genetic diversity and allelic frequency with multiple clones of msp-1, msp-2 and glurp in P. falciparum isolates in Khyber Pakhtunkhwa, Pakistan. In the present study the genotype data may provide valuable information essential for monitoring the impact of malaria eradication efforts in this region.

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.

The duration of chemoprophylaxis against malaria after treatment with artesunate-amodiaquine and artemether-lumefantrine and the effects of pfmdr1 86Y and pfcrt 76T: a meta-analysis of individual patient data

BACKGROUND

The majority of Plasmodium falciparum malaria cases in Africa are treated with the artemisinin combination therapies artemether-lumefantrine (AL) and artesunate-amodiaquine (AS-AQ), with amodiaquine being also widely used as part of seasonal malaria chemoprevention programs combined with sulfadoxine-pyrimethamine. While artemisinin derivatives have a short half-life, lumefantrine and amodiaquine may give rise to differing durations of post-treatment prophylaxis, an important additional benefit to patients in higher transmission areas.

METHODS

We analyzed individual patient data from 8 clinical trials of AL versus AS-AQ in 12 sites in Africa (n = 4214 individuals). The time to PCR-confirmed reinfection after treatment was used to estimate the duration of post-treatment protection, accounting for variation in transmission intensity between settings using hidden semi-Markov models. Accelerated failure-time models were used to identify potential effects of covariates on the time to reinfection. The estimated duration of chemoprophylaxis was then used in a mathematical model of malaria transmission to determine the potential public health impact of each drug when used for first-line treatment.

RESULTS

We estimated a mean duration of post-treatment protection of 13.0 days (95% CI 10.7-15.7) for AL and 15.2 days (95% CI 12.8-18.4) for AS-AQ overall. However, the duration varied significantly between trial sites, from 8.7-18.6 days for AL and 10.2-18.7 days for AS-AQ. Significant predictors of time to reinfection in multivariable models were transmission intensity, age, drug, and parasite genotype. Where wild type pfmdr1 and pfcrt parasite genotypes predominated (<=20% 86Y and 76T mutants, respectively), AS-AQ provided ~ 2-fold longer protection than AL. Conversely, at a higher prevalence of 86Y and 76T mutant parasites (> 80%), AL provided up to 1.5-fold longer protection than AS-AQ. Our simulations found that these differences in the duration of protection could alter population-level clinical incidence of malaria by up to 14% in under-5-year-old children when the drugs were used as first-line treatments in areas with high, seasonal transmission.

CONCLUSION

Choosing a first-line treatment which provides optimal post-treatment prophylaxis given the local prevalence of resistance-associated markers could make a significant contribution to reducing malaria morbidity.

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.

Emerging implications of policies on malaria treatment: genetic changes in the Pfmdr-1 gene affecting susceptibility to artemether-lumefantrine and artesunate-amodiaquine in Africa

Artemether–lumefantrine (AL) and artesunate–amodiaquine (AS-AQ) are the most commonly used artemisinin-based combination therapies (ACT) for treatment of Plasmodium falciparum in Africa. Both treatments remain efficacious, but single nucleotide polymorphisms (SNPs) in the Plasmodium falciparum multidrug resistance 1 (Pfmdr1) gene may compromise sensitivity. AL and AS-AQ exert opposing selective pressures: parasites with genotype 86Y, Y184 and 1246Y are partially resistant to AS-AQ treatment, while N86, 184 F and D1246 are favoured by AL treatment. Through a systematic review, we identified 397 surveys measuring the prevalence of Pfmdr1 polymorphisms at positions 86 184 or 1246 in 30 countries in Africa. Temporal trends in SNP frequencies after introduction of AL or AS-AQ as first-line treatment were analysed in 32 locations, and selection coefficients estimated. We examined associations between antimalarial policies, consumption, transmission intensity and rate of SNP selection. 1246Y frequency decreased on average more rapidly in locations where national policy recommended AL (median selection coefficient(s) of −0.083), compared with policies of AS-AQ or both AL and AS-AQ (median s=−0.035 and 0.021, p<0.001 respectively). 86Y frequency declined markedly after ACT policy introduction, with a borderline significant trend for a more rapid decline in countries with AL policies (p=0.055). However, these trends could also be explained by a difference in initial SNP frequencies at the time of ACT introduction. There were non-significant trends for faster selection of N86 and D1246 in areas with higher AL consumption and no trend with transmission intensity. Recorded consumption of AS-AQ was low in the locations and times Pfmdr1 data were collected. SNP trends in countries with AL policies suggest a broad increase in sensitivity of parasites to AS-AQ, by 7–10 years after AL introduction. Observed rates of selection have implications for planning strategies to cycle drugs or use multiple first-line therapies to maintain drug efficacy.

Genetic Diversity of the Plasmodium falciparum Glutamate-Rich Protein R2 Region Before and Twelve Years after Introduction of Artemisinin Combination Therapies among Febrile Children in Nigeria

The genetic diversity of glutamate-rich protein (GLURP) R2 region in Plasmodium falciparum isolates collected before and 12 years after the introduction of artemisinin combination treatment of malaria in Osogbo, Osun State, Nigeria, was compared in this study. Blood samples were collected on filter paper in 2004 and 2015 from febrile children from ages 1-12 years. The R2 region of the GLURP gene was genotyped using nested polymerase chain reaction and by nucleotide sequencing. In all, 12 GLURP alleles were observed in a total of 199 samples collected in the two study years. The multiplicity of infection (MOI) marginally increased over the two study years; however, the differences were statistically insignificant (2004 samples MOI = 1.23 versus 2015 samples MOI = 1.47). Some alleles were stable in their prevalence, whereas two GLURP alleles, VIII and XI, showed considerable variability between both years. This variability was replicated when GLURP sequences from other regions were compared with ours. The expected heterozygosity (He) values (He = 0.87) were identical for the two groups. High variability in the rearrangement of the amino acid repeat units in the R2 region were observed, with the amino acid repeat sequence DKNEKGQHEIVEVEEILPE more prevalent in both years, compared with the two other repeat sequences observed in the study. The parasite population characterized in this study displayed extensive genetic diversity. The detailed genetic profile of the GLURP R2 region has the potential to help guide further epidemiological studies aimed toward the rational design of novel chemotherapies that are antagonistic toward malaria.

Impacts of Antimalarial Drugs on Plasmodium falciparum Drug Resistance Markers, Western Kenya, 2003-2015

Antimalarial drug resistance has threatened global malaria control since chloroquine (CQ)-resistant Plasmodium falciparum emerged in Asia in the 1950s. Understanding the impacts of changing antimalarial drug policy on resistance is critical for resistance management. Plasmodium falciparum isolates were collected from 2003 to 2015 in western Kenya and analyzed for genetic markers associated with resistance to CQ (Pfcrt), sulfadoxine-pyrimethamine (SP) (Pfdhfr/Pfdhps), and artemether-lumefantrine (AL) (PfKelch13/Pfmdr1) antimalarials. In addition, household antimalarial drug use surveys were administered. Pfcrt 76T prevalence decreased from 76% to 6% from 2003 to 2015. Pfdhfr/Pfdhps quintuple mutants decreased from 70% in 2003 to 14% in 2008, but increased to near fixation by 2015. SP "super resistant" alleles Pfdhps 581G and 613S/T were not detected in the 2015 samples that were assessed. The Pfmdr1 N86-184F-D1246 haplotype associated with decreased lumefantrine susceptibility increased significantly from 4% in 2005 to 51% in 2015. No PfKelch13 mutations that have been previously associated with artemisinin resistance were detected in the study populations. The increase in Pfdhfr/Pfdhps quintuple mutants that associates with SP resistance may have resulted from the increased usage of SP for intermittent preventative therapy in pregnancy (IPTp) and for malaria treatment in the community. Prevalent Pfdhfr/Pfdhps mutations call for careful monitoring of SP resistance and effectiveness of the current IPTp program in Kenya. In addition, the commonly occurring Pfmdr1 N86-184F-D1246 haplotype associated with increased lumefantrine tolerance calls for surveillance of AL efficacy in Kenya, as well as consideration for a rotating artemisinin-combination therapy regimen.

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

The effect of local variation in malaria transmission on the prevalence of sulfadoxine-pyrimethamine resistant haplotypes and selective sweep characteristics in Malawi

BACKGROUND

Persistence of sulfadoxine-pyrimethamine (SP) resistance has been described in an urban setting in Malawi where malaria transmission is relatively low. Higher malaria transmission is associated with greater genetic diversity and more frequent genetic recombination, which could lead to a more rapid re-emergence of SP-sensitive parasites, as well as more rapid degradation of selective sweeps. In this study, the impact of local variation in malaria transmission on the prevalence of SP-resistant haplotypes and selective sweep characteristics was investigated at an urban site with low parasite prevalence and two rural sites with moderate and high parasite prevalence.

METHODS

Samples from three sites with different parasite prevalence were genotyped for resistance markers within pfdhfr-ts and pfdhps and at microsatellites flanking these genes. Expected heterozygosity (He) was estimated to evaluate genetic diversity.

RESULTS

No difference in the prevalence of highly resistant DHFR 51I/59R/108N and DHPS 437G/540E was found between sites. Small differences in He flanking pfdhfr-ts and pfdhps were seen between rural-moderate and the other sites, as well as some shared haplotypes between the rural-high and urban-low sites.

CONCLUSIONS

The results do not show an effect of local variation in malaria transmission, as inferred from parasite prevalence, on SP-resistant haplotype prevalence.

The emerging threat of artemisinin resistance in malaria: focus on artemether-lumefantrine

The development of artemisinin resistance in the Greater Mekong Subregion poses a significant threat to malaria elimination. Artemisinin-based combination therapies including artemether-lumefantrine (AL) are recommended by WHO as first-line treatment for uncomplicated Plasmodium falciparum malaria. This article provides a comprehensive review of the existing and latest data as a basis for interpretation of observed variability in parasite sensitivity to AL over the last 5 years. Clinical efficacy and preclinical data from a range of endemic countries are summarized, including potential molecular markers of resistance. Overall, AL remains effective in the treatment of uncomplicated P. falciparum malaria in most regions. Establishing validated molecular markers for resistance and strict efficacy monitoring will reinforce timely updates of treatment policies.

Population genetics analysis during the elimination process of Plasmodium falciparum in Djibouti

Background

Case management of imported malaria within the context of malaria pre-elimination is increasingly considered to be relevant because of the risk of resurgence. The assessment of malaria importation would provide key data i) to select countries with propitious conditions for pre-elimination phase and ii) to predict its feasibility. Recently, a sero-prevalence study in Djibouti indicated low malaria prevalence, which is propitious for the implementation of pre-elimination, but data on the extent of malaria importation remain unknown.

Methods

Djiboutian plasmodial populations were analysed over an eleven-year period (1998, 1999, 2002 and 2009). The risk of malaria importation was indirectly assessed by using plasmodial population parameters. Based on 5 microsatellite markers, expected heterozygosity (H.e.), multiplicity of infection, pairwise Fst index, multiple correspondence analysis and individual genetic relationship were determined. The prevalence of single nucleotide polymorphisms associated with pyrimethamine resistance was also determined.

Results

Data indicated a significant decline in genetic diversity (0.51, 0.59, 0.51 and 0 in 1998, 1999, 2002 and 2009, respectively) over the study period, which is inconsistent with the level of malaria importation described in a previous study. This suggested that Djiboutian malaria situation may have benefited from the decline of malaria prevalence that occurred in neighbouring countries, in particular in Ethiopia. The high Fst indices derived from plasmodial populations from one study period to another (0.12 between 1999 and 2002, and 0.43 between 2002 and 2009) suggested a random sampling of parasites, probably imported from neighbouring countries, leading to oligo-clonal expansion of few different strains during each transmission season. Nevertheless, similar genotypes observed during the study period suggested recurrent migrations and imported malaria.

Conclusion

In the present study, the extent of genetic diversity was used to assess the risk of malaria importation in the low malaria transmission setting of Djibouti. The molecular approach highlights i) the evolution of Djiboutian plasmodial population profiles that are consistent and compatible with Djiboutian pre-elimination goals and ii) the necessity to implement the monitoring of plasmodial populations and interventions at the regional scale in the Horn of Africa to ensure higher efficiency of malaria control and elimination.