Measuring protective efficacy and quantifying the impact of drug resistance: A novel malaria chemoprevention trial design and methodology.

BACKGROUND: Recently revised WHO guidelines on malaria chemoprevention have opened the door to more tailored implementation. Countries face choices on whether to replace old drugs, target additional age groups, and adapt delivery schedules according to local drug resistance levels and malaria transmission patterns. Regular routine assessment of protective efficacy of chemoprevention is key. Here, we apply a novel modelling approach to aid the design and analysis of chemoprevention trials and generate measures of protection that can be applied across a range of transmission settings. METHODS AND FINDINGS: We developed a model of genotype-specific drug protection, which accounts for underlying risk of infection and circulating genotypes. Using a Bayesian framework, we fitted the model to multiple simulated scenarios to explore variations in study design, setting, and participant characteristics. We find that a placebo or control group with no drug protection is valuable but not always feasible. An alternative approach is a single-arm trial with an extended follow-up (>42 days), which allows measurement of the underlying infection risk after drug protection wanes, as long as transmission is relatively constant. We show that the currently recommended 28-day follow-up in a single-arm trial results in low precision of estimated 30-day chemoprevention efficacy and low power in determining genotype differences of 12 days in the duration of protection (power = 1.4%). Extending follow-up to 42 days increased precision and power (71.5%) in settings with constant transmission over this time period. However, in settings of unstable transmission, protective efficacy in a single-arm trial was overestimated by 24.3% if recruitment occurred during increasing transmission and underestimated by 15.8% when recruitment occurred during declining transmission. Protective efficacy was estimated with greater precision in high transmission settings, and power to detect differences by resistance genotype was lower in scenarios where the resistant genotype was either rare or too common. CONCLUSIONS: These findings have important implications for the current guidelines on chemoprevention efficacy studies and will be valuable for informing where these studies should be optimally placed. The results underscore the need for a comparator group in seasonal settings and provide evidence that the extension of follow-up in single-arm trials improves the accuracy of measures of protective efficacy in settings with more stable transmission. Extension of follow-up may pose logistical challenges to trial feasibility and associated costs. However, these studies may not need to be repeated multiple times, as the estimates of drug protection against different genotypes can be applied to different settings by adjusting for transmission intensity and frequency of resistance.

Identification of the PfK13 mutations R561H and P441L in the Democratic Republic of Congo.

OBJECTIVES: Partial artemisinin resistance, mediated by Plasmodium falciparum K13 (PfK13) mutations, has been confirmed in certain areas of East Africa that are historically associated with high-level antimalarial resistance. The Democratic Republic of Congo (DRC) borders these areas in the East. This study aimed to determine the prevalence of resistance markers in six National Malaria Control Program surveillance sites; Boende, Kabondo, Kapolowe, Kimpese, Mikalayi, and Rutshuru. METHODS: The single nucleotide polymorphisms (SNPs) in P. falciparum genes PfK13, Pfdhfr, Pfdhps, Pfmdr1, and Pfcrt were assessed using targeted next-generation sequencing of isolates collected at enrollment in therapeutic efficacy studies. RESULTS: PfK13 SNPs were detected in two samples: in Kabondo (R561H) and in Rutshuru (P441L), both areas near Uganda and Rwanda. The Pfdhps ISGEGA haplotype, associated with reduced sulfadoxine-pyrimethamine chemoprevention efficacy, ranged from 0.8% in Mikalayi (central DRC) to 42.2% in Rutshuru (East DRC). CONCLUSIONS: R561H and P441L observed in eastern DRC are a concern, as they are associated with delayed artemisinin-based combination therapies-clearance and candidate marker of resistance, respectively. This is consistent with previous observations of shared drug resistance profiles in parasites of that region with bordering areas of Rwanda and Uganda. The likely circulation of parasites has important implications for the ongoing surveillance of partial artemisinin-resistant P. falciparum and for future efforts to mitigate its dispersal.

A snapshot of the prevalence of dihydropteroate synthase-431V mutation and other sulfadoxine-pyrimethamine resistance markers in Plasmodium falciparum isolates in Nigeria.

BACKGROUND: Malaria is a major public health issue with substantial risks among vulnerable populations. Currently, the World Health Organization (WHO) recommends SP-IPTp in the second and third trimesters. However, the efficacy of SP-IPTp is threatened by the emergence of sulfadoxine-pyrimethamine resistant malaria parasites due to single nucleotide polymorphisms in the Plasmodium falciparum dihydrofolate reductase and dihydropteroate synthetase genes. This study aimed to assess the current prevalence of Pfdhfr/Pfdhps mutations in P. falciparum isolates collected from individuals residing in Ile-Ife, Nigeria, and also present maps of the prevalence of Pfdhps 431V and 581G within Nigeria and surrounding countries. METHODS: Between October 2020 and April 2021, samples were collected as dried blood spots among 188 participants who showed malaria positivity with a histidine-rich-protein-based rapid diagnostic test (RDT). Nested PCR assays were used to confirm falciparum in the samples with RDT positivity, and to amplify fragments of the Pfdhfr/Pfdhps genes followed by targeted amplicon sequencing. Published data since 2007 on the prevalence of the Pfdhps genotypes in Nigeria and the neighbouring countries were used to produce maps to show the distribution of the mutant genotypes. RESULTS: Only 74 and 61 samples were successfully amplified for the Pfdhfr and Pfdhps genes, respectively. At codons resulting in N51I, C59R, and S108N, Pfdhfr carried mutant alleles of 97.3% (72/74), 97.3% (72/74) and 98.6% (73/74), respectively. The Pfdhps gene carried mutations at codons resulting in amino acid changes at 431-436-437-540-581-613; I431V [45.9%, (28/61)], A581G [31.1% (19/61)] and A613S [49.2% (30/61)]. Constructed haplotypes were mainly the triple Pfdhfr mutant 51I-59R-108N (95.9%), and the most common haplotypes observed for the Pfdhps gene were the ISGKAA (32.8%), ISGKGS (8.2%), VAGKAA (14.8%), VAGKAS (9.8%) and VAGKGS (14.8%). In the context of the previously published data, a high prevalence of 431V/581G mutations was found in the study population. It seems quite evident that the Pfdhps 431V, 581G and 613S often co-occur as Pfdhps-VAGKGS haplotype. CONCLUSION: This study showed that the prevalence of VAGKGS haplotype seems to be increasing in prevalence. If this is similar in effect to the emergence of 581G in East Africa, the efficacy of SP-IPTp in the presence of these novel Pfdhps mutants should be re-assessed.

Identification of Single-Nucleotide Polymorphisms in the Mitochondrial Genome and Kelch 13 Gene of Plasmodium falciparum in Different Geographical Populations.

The emergence of artemisinin-resistant Plasmodium falciparum parasites in Southeast Asia threatens malaria control and elimination. The interconnectedness of parasite populations may be essential to monitor the spread of resistance. Combining a published barcoding system of geographically restricted single-nucleotide polymorphisms (SNPs), mainly mitochondria of P. falciparum with SNPs in the K13 artemisinin resistance marker, could elucidate the parasite population structure and provide insight regarding the spread of drug resistance. We explored the diversity of mitochondrial SNPs (bp position 611-2825) and identified K13 SNPs from malaria patients in the districts of India (Ranchi), Tanzania (Korogwe), and Senegal (Podor, Richard Toll, Kaolack, and Ndoffane). DNA was amplified using a nested PCR and Sanger-sequenced. Overall, 199 K13 sequences (India: N = 92; Tanzania: N = 48; Senegal: N = 59) and 237 mitochondrial sequences (India: N = 93; Tanzania: N = 48; Senegal: N = 96) were generated. SNPs were identified by comparisons with reference genomes. We detected previously reported geographically restricted mitochondrial SNPs (T2175C and G1367A) as markers for parasites originating from the Indian subcontinent and several geographically unrestricted mitochondrial SNPs. Combining haplotypes with published P. falciparum mitochondrial genome data suggested possible regional differences within India. All three countries had G1692A, but Tanzanian and Senegalese SNPs were well-differentiated. Some mitochondrial SNPs are reported here for the first time. Four nonsynonymous K13 SNPs were detected: K189T (India, Tanzania, Senegal); A175T (Tanzania); and A174V and R255K (Senegal). This study supports the use of mitochondrial SNPs to determine the origin of the parasite and suggests that the P. falciparum populations studied were susceptible to artemisinin during sampling because all K13 SNPs observed were outside the propeller domain for artemisinin resistance.

Selective whole genome amplification of Plasmodium malariae DNA from clinical samples reveals insights into population structure.

The genomic diversity of Plasmodium malariae malaria parasites is understudied, partly because infected individuals tend to present with low parasite densities, leading to difficulties in obtaining sufficient parasite DNA for genome analysis. Selective whole genome amplification (SWGA) increases the relative levels of pathogen DNA in a clinical sample, but has not been adapted for P. malariae parasites. Here we design customized SWGA primers which successfully amplify P. malariae DNA extracted directly from unprocessed clinical blood samples obtained from patients with P. malariae-mono-infections from six countries, and further test the efficacy of SWGA on mixed infections with other Plasmodium spp. SWGA enables the successful whole genome sequencing of samples with low parasite density (i.e. one sample with a parasitaemia of 0.0064% resulted in 44% of the genome covered by ≥ 5 reads), leading to an average 14-fold increase in genome coverage when compared to unamplified samples. We identify a total of 868,476 genome-wide SNPs, of which 194,709 are unique across 18 high-quality isolates. After exclusion of the hypervariable subtelomeric regions, a high-quality core subset of 29,899 unique SNPs is defined. Population genetic analysis suggests that P. malariae parasites display clear geographical separation by continent. Further, SWGA successfully amplifies genetic regions of interest such as orthologs of P. falciparum drug resistance-associated loci (Pfdhfr, Pfdhps, Pfcrt, Pfk13 and Pfmdr1), and several non-synonymous SNPs were detected in these genes. In conclusion, we have established a robust SWGA approach that can assist whole genome sequencing of P. malariae, and thereby facilitate the implementation of much-needed large-scale multi-population genomic studies of this neglected malaria parasite. As demonstrated in other Plasmodia, such genetic diversity studies can provide insights into the biology underlying the disease and inform malaria surveillance and control measures.

A molecular barcode to inform the geographical origin and transmission dynamics of Plasmodium vivax malaria.

Although Plasmodium vivax parasites are the predominant cause of malaria outside of sub-Saharan Africa, they not always prioritised by elimination programmes. P. vivax is resilient and poses challenges through its ability to re-emerge from dormancy in the human liver. With observed growing drug-resistance and the increasing reports of life-threatening infections, new tools to inform elimination efforts are needed. In order to halt transmission, we need to better understand the dynamics of transmission, the movement of parasites, and the reservoirs of infection in order to design targeted interventions. The use of molecular genetics and epidemiology for tracking and studying malaria parasite populations has been applied successfully in P. falciparum species and here we sought to develop a molecular genetic tool for P. vivax. By assembling the largest set of P. vivax whole genome sequences (n = 433) spanning 17 countries, and applying a machine learning approach, we created a 71 SNP barcode with high predictive ability to identify geographic origin (91.4%). Further, due to the inclusion of markers for within population variability, the barcode may also distinguish local transmission networks. By using P. vivax data from a low-transmission setting in Malaysia, we demonstrate the potential ability to infer outbreak events. By characterising the barcoding SNP genotypes in P. vivax DNA sourced from UK travellers (n = 132) to ten malaria endemic countries predominantly not used in the barcode construction, we correctly predicted the geographic region of infection origin. Overall, the 71 SNP barcode outperforms previously published genotyping methods and when rolled-out within new portable platforms, is likely to be an invaluable tool for informing targeted interventions towards elimination of this resilient human malaria.

Markers of sulfadoxine-pyrimethamine resistance in Eastern Democratic Republic of Congo; implications for malaria chemoprevention.

BACKGROUND: Sulfadoxine-pyrimethamine (SP) is a cornerstone of malaria chemoprophylaxis and is considered for programmes in the Democratic Republic of Congo (DRC). However, SP efficacy is threatened by drug resistance, that is conferred by mutations in the dhfr and dhps genes. The World Health Organization has specified that intermittent preventive treatment for infants (IPTi) with SP should be implemented only if the prevalence of the dhps K540E mutation is under 50%. There are limited current data on the prevalence of resistance-conferring mutations available from Eastern DRC. The current study aimed to address this knowledge gap. METHODS: Dried blood-spot samples were collected from clinically suspected malaria patients [outpatient department (OPD)] and pregnant women attending antenatal care (ANC) in four sites in North and South Kivu, DRC. Quantitative PCR (qPCR) was performed on samples from individuals with positive and with negative rapid diagnostic test (RDT) results. Dhps K450E and A581G and dhfr I164L were assessed by nested PCR followed by allele-specific primer extension and detection by multiplex bead-based assays. RESULTS: Across populations, Plasmodium falciparum parasite prevalence was 47.9% (1160/2421) by RDT and 71.7 (1763/2421) by qPCR. Median parasite density measured by qPCR in RDT-negative qPCR-positive samples was very low with a median of 2.3 parasites/µL (IQR 0.5-25.2). Resistance genotyping was successfully performed in RDT-positive samples and RDT-negative/qPCR-positive samples with success rates of 86.2% (937/1086) and 55.5% (361/651), respectively. The presence of dhps K540E was high across sites (50.3-87.9%), with strong evidence for differences between sites (p < 0.001). Dhps A581G mutants were less prevalent (12.7-47.2%). The dhfr I164L mutation was found in one sample. CONCLUSIONS: The prevalence of the SP resistance marker dhps K540E exceeds 50% in all four study sites in North and South Kivu, DRC. K540E mutations regularly co-occurred with mutations in dhps A581G but not with the dhfr I164L mutation. The current results do not support implementation of IPTi with SP in the study area.

Deployment and utilization of next-generation sequencing of Plasmodium falciparum to guide anti-malarial drug policy decisions in sub-Saharan Africa: opportunities and challenges.

Parasite resistance against anti-malarial drugs is a major threat to the ongoing malaria control and elimination strategies. This is especially true since resistance to the currently recommended artemisinins and partner drugs has been confirmed in South East Asia (SEA) and new anti-malarial compounds are not expected to be available in the near future. Spread from SEA or independent emergence of artemisinin resistance in sub-Saharan Africa (SSA) could reverse the achievements in malaria control that have been attained in the past two decades and derail the ongoing elimination strategies. The current surveillance of clinical efficacy and resistance to anti-malarial drugs is based on efficacy trials to assess the clinical performance of anti-malarials, in vivo/ex vivo assessment of parasite susceptibility to anti-malarials and prevalence of known molecular markers of drug resistance. Whereas clinical efficacy trials are restricted by cost and the complex logistics of patient follow-up, molecular detection of genetic mutations associated with resistance or reduced susceptibility to anti-malarials is by contrast a simple and powerful tool for early detection and monitoring of the prevalence of resistant parasites at population level. This provides needed information before clinical failure emerges, allowing policy makers to anticipate problems and respond. The various methods previously used in detection of molecular markers of drug resistance share some limitations: low-throughput, and high costs per sample and demanding infrastructure. However, recent technological advances including next-generation sequencing (NGS) methodologies promise greatly increased throughput and reduced costs, essentially providing unprecedented potential to address different research and operational questions of relevance for drug policy. This review assesses the potential role of NGS to provide comprehensive information that could guide drug policies in malaria endemic countries and looks at the foreseeable challenges facing the establishment of NGS approaches for routine surveillance of parasite resistance to anti-malarials in SSA.

Effect of Plasmodium falciparum sulfadoxine-pyrimethamine resistance on the effectiveness of intermittent preventive therapy for malaria in pregnancy in Africa: a systematic review and meta-analysis.

BACKGROUND: Resistance of Plasmodium falciparum to sulfadoxine-pyrimethamine threatens the antimalarial effectiveness of intermittent preventive treatment during pregnancy (IPTp) in sub-Saharan Africa. We aimed to assess the associations between markers of sulfadoxine-pyrimethamine resistance in P falciparum and the effectiveness of sulfadoxine-pyrimethamine IPTp for malaria-associated outcomes. METHODS: For this systematic review and meta-analysis, we searched databases (from Jan 1, 1990 to March 1, 2018) for clinical studies (aggregated data) or surveys (individual participant data) that reported data on low birthweight (primary outcome) and malaria by sulfadoxine-pyrimethamine IPTp dose, and for studies that reported on molecular markers of sulfadoxine-pyrimethamine resistance. Studies that involved only HIV-infected women or combined interventions were excluded. We did a random-effects meta-analysis (clinical studies) or multivariate log-binomial regression (surveys) to obtain summarised dose-response data (relative risk reduction [RRR]) and multivariate meta-regression to explore the modifying effects of sulfadoxine-pyrimethamine resistance (as indicated by Ala437Gly, Lys540Glu, and Ala581Gly substitutions in the dhps gene). This study is registered with PROSPERO, number 42016035540. FINDINGS: Of 1097 records screened, 57 studies were included in the aggregated-data meta-analysis (including 59 457 births). The RRR for low birthweight declined with increasing prevalence of dhps Lys540Glu (ptrend=0·0060) but not Ala437Gly (ptrend=0·35). The RRR was 7% (95% CI 0 to 13) in areas of high resistance to sulfadoxine-pyrimethamine (Lys540Glu ≥90% in east and southern Africa; n=11), 21% (14 to 29) in moderate-resistance areas (Ala437Gly ≥90% [central and west Africa], or Lys540Glu ≥30% to <90% [east and southern Africa]; n=16), and 27% (21 to 33) in low-resistance areas (Ala437Gly <90% [central and west Africa], or Lys540Glu <30% [east and southern Africa]; n=30; ptrend=0·0054 [univariate], I2=69·5%). The overall RRR in all resistance strata was 21% (17 to 25). In the analysis of individual participant data from 13 surveys (42 394 births), sulfadoxine-pyrimethamine IPTp was associated with reduced prevalence of low birthweight in areas with a Lys540Glu prevalence of more than 90% and Ala581Gly prevalence of less than 10% (RRR 10% [7 to 12]), but not in those with an Ala581Gly prevalence of 10% or higher (pooled Ala581Gly prevalence 37% [range 29 to 46]; RRR 0·5% [-16 to 14]; 2326 births). INTERPRETATION: The effectiveness of sulfadoxine-pyrimethamine IPTp is reduced in areas with high resistance to sulfadoxine-pyrimethamine among P falciparum parasites, but remains associated with reductions in low birthweight even in areas where dhps Lys540Glu prevalence exceeds 90% but where the sextuple-mutant parasite (harbouring the additional dhps Ala581Gly mutation) is uncommon. Therapeutic alternatives to sulfadoxine-pyrimethamine IPTp are needed in areas where the prevalence of the sextuple-mutant parasite exceeds 37%. FUNDING: US Centers for Disease Control and Prevention, the Malaria in Pregnancy Consortium (funded through a grant from the Bill & Melinda Gates Foundation to the Liverpool School of Tropical Medicine), Worldwide Antimalarial Resistance Network, European and Developing Countries Clinical Trials Partnership.

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.

Global genetic diversity of var2csa in Plasmodium falciparum with implications for malaria in pregnancy and vaccine development.

Malaria infection during pregnancy, caused by the sequestering of Plasmodium falciparum parasites in the placenta, leads to high infant mortality and maternal morbidity. The parasite-placenta adherence mechanism is mediated by the VAR2CSA protein, a target for natural occurring immunity. Currently, vaccine development is based on its ID1-DBL2Xb domain however little is known about the global genetic diversity of the encoding var2csa gene, which could influence vaccine efficacy. In a comprehensive analysis of the var2csa gene in >2,000 P. falciparum field isolates across 23 countries, we found that var2csa is duplicated in high prevalence (>25%), African and Oceanian populations harbour a much higher diversity than other regions, and that insertions/deletions are abundant leading to an underestimation of the diversity of the locus. Further, ID1-DBL2Xb haplotypes associated with adverse birth outcomes are present globally, and African-specific haplotypes exist, which should be incorporated into vaccine design.

Molecular assays for antimalarial drug resistance surveillance: A target product profile.

Antimalarial drug resistance is a major constraint for malaria control and elimination efforts. Artemisinin-based combination therapy is now the mainstay for malaria treatment. However, delayed parasite clearance following treatment with artemisinin derivatives has now spread in the Greater Mekong Sub region and may emerge or spread to other malaria endemic regions. This spread is of great concern for malaria control programmes, as no alternatives to artemisinin-based combination therapies are expected to be available in the near future. There is a need to strengthen surveillance systems for early detection and response to the antimalarial drug resistance threat. Current surveillance is mainly done through therapeutic efficacy studies; however these studies are complex and both time- and resource-intensive. For multiple common antimalarials, parasite drug resistance has been correlated with specific genetic mutations, and the molecular markers associated with antimalarial drug resistance offer a simple and powerful tool to monitor the emergence and spread of resistant parasites. Different techniques to analyse molecular markers associated with antimalarial drug resistance are available, each with advantages and disadvantages. However, procedures are not adequately harmonized to facilitate comparisons between sites. Here we describe the target product profiles for tests to analyse molecular markers associated with antimalarial drug resistance, discuss how use of current techniques can be standardised, and identify the requirements for an ideal product that would allow malaria endemic countries to provide useful spatial and temporal information on the spread of resistance.

A reference genome and methylome for the Plasmodium knowlesi A1-H.1 line.

Plasmodium knowlesi, a common parasite of macaques, is recognised as a significant cause of human malaria in Malaysia. The P. knowlesi A1H1 line has been adapted to continuous culture in human erythrocytes, successfully providing an in vitro model to study the parasite. We have assembled a reference genome for the PkA1-H.1 line using PacBio long read combined with Illumina short read sequence data. Compared with the H-strain reference, the new reference has improved genome coverage and a novel description of methylation sites. The PkA1-H.1 reference will enhance the capabilities of the in vitro model to improve the understanding of P. knowlesi infection in humans.

Analysis of nuclear and organellar genomes of Plasmodium knowlesi in humans reveals ancient population structure and recent recombination among host-specific subpopulations.

The macaque parasite Plasmodium knowlesi is a significant concern in Malaysia where cases of human infection are increasing. Parasites infecting humans originate from genetically distinct subpopulations associated with the long-tailed (Macaca fascicularis (Mf)) or pig-tailed macaques (Macaca nemestrina (Mn)). We used a new high-quality reference genome to re-evaluate previously described subpopulations among human and macaque isolates from Malaysian-Borneo and Peninsular-Malaysia. Nuclear genomes were dimorphic, as expected, but new evidence of chromosomal-segment exchanges between subpopulations was found. A large segment on chromosome 8 originating from the Mn subpopulation and containing genes encoding proteins expressed in mosquito-borne parasite stages, was found in Mf genotypes. By contrast, non-recombining organelle genomes partitioned into 3 deeply branched lineages, unlinked with nuclear genomic dimorphism. Subpopulations which diverged in isolation have re-connected, possibly due to deforestation and disruption of wild macaque habitats. The resulting genomic mosaics reveal traits selected by host-vector-parasite interactions in a setting of ecological transition.

Mapping sulphadoxine-pyrimethamine-resistant Plasmodium falciparum malaria in infected humans and in parasite populations in Africa.

Intermittent preventive treatment (IPT) with sulphadoxine-pyrimethamine in vulnerable populations reduces malaria morbidity in Africa, but resistance mutations in the parasite dhps gene (combined with dhfr mutations) threaten its efficacy. We update a systematic review to map the prevalence of K540E and A581G mutations in 294 surveys of infected humans across Africa from 2004-present. Interpreting these data is complicated by multiclonal infections in humans, especially in high transmission areas. We extend statistical methods to estimate the frequency, i.e. the proportion of resistant clones in the parasite population at each location, and so standardise for varying transmission levels. Both K540E and A581G mutations increased in prevalence and frequency in 60% of areas after 2008, highlighting the need for ongoing surveillance. Resistance measures within countries were similar within 300 km, suggesting an appropriate spatial scale for surveillance. Spread of the mutations tended to accelerate once their prevalence exceeded 10% (prior to fixation). Frequencies of resistance in parasite populations are the same or lower than prevalence in humans, so more areas would be classified as likely to benefit from IPT if similar frequency thresholds were applied. We propose that the use of resistance frequencies as well as prevalence measures for policy decisions should be evaluated.

Genomic variation in Plasmodium vivax malaria reveals regions under selective pressure.

BACKGROUND: Although Plasmodium vivax contributes to almost half of all malaria cases outside Africa, it has been relatively neglected compared to the more deadly P. falciparum. It is known that P. vivax populations possess high genetic diversity, differing geographically potentially due to different vector species, host genetics and environmental factors. RESULTS: We analysed the high-quality genomic data for 46 P. vivax isolates spanning 10 countries across 4 continents. Using population genetic methods we identified hotspots of selection pressure, including the previously reported MRP1 and DHPS genes, both putative drug resistance loci. Extra copies and deletions in the promoter region of another drug resistance candidate, MDR1 gene, and duplications in the Duffy binding protein gene (PvDBP) potentially involved in erythrocyte invasion, were also identified. For surveillance applications, continental-informative markers were found in putative drug resistance loci, and we show that organellar polymorphisms could classify P. vivax populations across continents and differentiate between Plasmodia spp. CONCLUSIONS: This study has shown that genomic diversity that lies within and between P. vivax populations can be used to elucidate potential drug resistance and invasion mechanisms, as well as facilitate the molecular barcoding of the parasite for surveillance applications.

Genetic diversity of next generation antimalarial targets: A baseline for drug resistance surveillance programmes.

Drug resistance is a recurrent problem in the fight against malaria. Genetic and epidemiological surveillance of antimalarial resistant parasite alleles is crucial to guide drug therapies and clinical management. New antimalarial compounds are currently at various stages of clinical trials and regulatory evaluation. Using ∼2000 Plasmodium falciparum genome sequences, we investigated the genetic diversity of eleven gene-targets of promising antimalarial compounds and assessed their potential efficiency across malaria endemic regions. We determined if the loci are under selection prior to the introduction of new drugs and established a baseline of genetic variance, including potential resistant alleles, for future surveillance programmes.

Travel and the emergence of high-level drug resistance in Plasmodium falciparum in southwest Uganda: results from a population-based study.

BACKGROUND: The I164L mutation on the dhfr gene confers high level resistance to sulfadoxine-pyrimethamine (SP) but it is rare in Africa except in a cluster of reports where prevalence >10% in highland areas of southwest Uganda and eastern Rwanda. The occurrence of the dhfr I164L mutation was investigated in community surveys in this area and examined the relationship to migration. METHODS: A cross-sectional prevalence survey was undertaken in among villages within the catchment areas of two health facilities in a highland site (Kabale) and a highland fringe site (Rukungiri) in 2007. Sociodemographic details, including recent migration, were collected for each person included in the study. A total of 206 Plasmodium falciparum positive subjects were detected by rapid diagnostic test; 203 in Rukungiri and 3 in Kabale. Bloodspot samples were taken and were screened for dhfr I164L. RESULTS: Sequence analysis confirmed the presence of the I164L mutations in twelve P. falciparum positive samples giving an estimated prevalence of 8.6% in Rukungiri. Of the three parasite positive samples in Kabale, none had I164L mutations. Among the twelve I164L positives three were male, ages ranged from 5 to 90 years of age. None of those with the I164L mutation had travelled in the 8 weeks prior to the survey, although three were from households from which at least one household member had travelled during that period. Haplotypes were determined in non-mixed infections and showed the dhfr I164L mutation occurs in both as a N51I + S108N + I164L haplotype (n = 2) and N51I + C59R + S108N + I164L haplotype (n = 5). Genotyping of flanking microsatellite markers showed that the I164L occurred independently on the triple mutant (N51I, C59R + S108N) and double mutant (N51I + S108N) background. CONCLUSIONS: There is sustained local transmission of parasites with the dhfr I164L mutation in Rukungiri and no evidence to indicate its occurrence is associated with recent travel to highly resistant neighbouring areas. The emergence of a regional cluster of I164L in SW Uganda and Rwanda indicates that transmission of I164L is facilitated by strong drug pressure in low transmission areas potentially catalysed in those areas by travel and the importation of parasites from relatively higher transmission settings.

Unexpected selections of Plasmodium falciparum polymorphisms in previously treatment-naïve areas after monthly presumptive administration of three different anti-malarial drugs in Liberia 1976-78.

BACKGROUND: To assess the effect on malaria prevalence, village specific monthly administrations of pyrimethamine, chlorproguanil, chloroquine or placebo were given to children in four previously treatment-naïve Liberian villages, 1976-78. Plasmodium falciparum in vivo resistance developed to pyrimethamine only. Selection of molecular markers of P. falciparum resistance after 2 years of treatment are reported. METHODS: Blood samples were collected from 191 study children in a survey in 1978. Polymorphisms in pfcrt, pfmdr1, pfdhfr, pfdhps, pfmrp1 and pfnhe1 genes were determined using PCR-based methods. RESULTS: Pfcrt 72-76 CVIET was found in one chloroquine village sample, all remaining samples had pfcrt CVMNK. Pfmdr1 N86 prevalence was 100%. A pfmdr1 T1069ACT→ACG synonymous polymorphism was found in 30% of chloroquine village samples and 3% of other samples (P = 0.008). Variations in pfnhe1 block I were found in all except the chloroquine treated village (P < 0.001). Resistance associated pfdhfr 108N prevalence was 2% in the pyrimethamine village compared to 45-65% elsewhere, including the placebo village (P = 0.001). CONCLUSIONS: Chloroquine treatment possibly resulted in the development of pfcrt 72-76 CVIET. Selection of pfmdr1 T1069ACG and a pfnhe1 block 1 genotypes indicates that chloroquine treatment exerted a selective pressure on P. falciparum. Pyrimethamine resistance associated pfdhfr 108N was present prior to the introduction of any drug. Decreased pfdhfr 108N frequency concurrent with development of pyrimethamine resistance suggests a non-pfdhfr polymorphisms mediated resistance mechanism.