Pharmacokinetics of single low dose primaquine in Ugandan and Congolese children with falciparum malaria

Background: There are no pharmacokinetic data of single low dose primaquine (SLDPQ) as transmission blocking in African children with acute Plasmodium falciparum and glucose-6-phosphate dehydrogenase deficiency (G6PDd). Methods: Primaquine pharmacokinetics of age-dosed SLDPQ (shown previously to be gametocytocidal with similar tolerability as placebo) were characterised in falciparum-infected Ugandan and Congolese children aged 6 months to 11 years, treated on admission with standard 3-day dihydroartemisinin-piperaquine or artemether-lumefantrine plus SLDPQ: 6 m–<1 y: 1.25 mg, 1–5 y: 2.5 mg, 6–9 y: 5 mg, 10–11 y: 7.5 mg. LC-MS/MS-measured plasma primaquine and carboxyprimaquine (baseline, 1, 1.5, 2, 4, 8, 12, 24 h) were analysed by noncompartmental analysis. Multivariable linear regression modelled associations between covariates, including cytochrome-P450 2D6 metaboliser status, and outcomes. Findings: 258 children (median age 5 [interquartile range (IQR) 3–7]) were sampled; 8 (3.1%) with early vomiting were excluded. Primaquine doses of 0.10–0.40 (median 0.21, IQR 0.16–0.25) mg base/kg resulted in primaquine maximum plasma concentrations (Cmax) of 2.3–447 (median 103.0, IQR 72.1–140.0) ng/mL between 1.0 and 8.0 (median 2) hours (Tmax) and median areas under the drug concentration curves (AUC0-last) 730.2 (6 m–<1 y, n = 12), 582.8 (1–5 y, n = 126), 871.1 (6–9 y, n = 80), and 931.0 (10–11 y, n = 32) ng∗h/mL. Median elimination half-live (T½) was 4.7 (IQR 3.8–5.6) hours. Primaquine clearance/kg peaked at 18 months, plateauing at 4 y. Increasing CYP2D6 metaboliser activity score [poor (3/250), intermediate (52/250), normal (150/250), ultrarapid (5/250), indeterminate (40/250)] and baseline haemoglobin were significantly associated with a lower primaquine AUC0-last,which increased with increasing mg/kg dose and age but was independent of the artemisinin treatment used. Interpretation: Age-dosed SLDPQ resulted in variable primaquine exposure that depended on bodyweight-adjusted dose, age, baseline haemoglobin and CYP2D6 metaboliser status, but not on dihydroartemisinin-piperaquine or artemether-lumefantrine. These data support age-dosed SLDPQ for transmission blocking in sub-Saharan Africa. Funding: This work was cofunded by the UK Medical Research Council, Wellcome Trust, and UK Aid through the Global Health Trials (grant reference MR/P006973/1). The funders had no role in the study design, execution, and analysis and decisions regarding publication

Measurement of gene amplifications related to drug resistance in Plasmodium falciparum using droplet digital PCR

Background Copy number variations (CNVs) of the Plasmodium falciparum multidrug resistance 1 (pfmdr1), P. falciparum plasmepsin2 (pfplasmepsin2) and P. falciparum GTP cyclohydrolase 1 (pfgch1) genes are associated with anti-malarial drug resistance in P. falciparum malaria. Droplet digital PCR (ddPCR) assays have been developed for accurate assessment of CNVs in several human genes. The aim of the present study was to develop and validate ddPCR assays for detection of the CNVs of P. falciparum genes associated with resistance to anti-malarial drugs. Methods A multiplex ddPCR assay was developed to detect the CNVs in the pfmdr1 and pfplasmepsin2 genes, while a duplex ddPCR assay was developed to detect CNV in the pfgch1 gene. The gene copy number (GCN) quantification limit, as well as the accuracy and precision of the ddPCR assays were determined and compared to conventional quantitative PCR (qPCR). In order to reduce the cost of testing, a multiplex ddPCR assay of two target genes, pfmdr1 and pfplasmepsin2, was validated. In addition, the CNVs of genes of field samples collected from Thailand from 2015 to 2019 (n = 84) were assessed by ddPCR and results were compared to qPCR as the reference assay. Results There were no significant differences between the GCN results obtained from uniplex and multiplex ddPCR assays for detection of CNVs in the pfmdr1 and pfplasmepsin2 genes (p = 0.363 and 0.330, respectively). Based on the obtained gene copy number quantification limit, the accuracy and percent relative standard deviation (%RSD) value of the multiplex ddPCR assay were 95% and 5%, respectively, for detection of the CNV of the pfmdr1 gene, and 91% and 5% for detection of the CNV of the pfplasmepsin2 gene. There was no significant difference in gene copy numbers assessed by uniplex or duplex ddPCR assays regarding CNV in the pfgch1 gene (p = 0.276). The accuracy and %RSD value of the duplex ddPCR assay were 95% and 4%, respectively, regarding pfgch1 GCN. In the P. falciparum field samples, pfmdr1 and pfplasmepsin2 GCNs were amplified in 15% and 27% of samples from Ubon Ratchathani, Thailand, while pfgch1 GCN was amplified in 50% of samples from Yala, Thailand. There was 100% agreement between the GCN results obtained from the ddPCR and qPCR assays (κ = 1.00). The results suggested that multiplex ddPCR assay is the optional assay for the accurate detection of gene copy number without requiring calibration standards, while the cost and required time are reduced. Based on the results of this study, criteria for GCN detection by ddPCR analysis were generated. Conclusions The developed ddPCR assays are simple, accurate, precise and cost-effective tools for detection of the CNVs in the pfmdr1, pfplasmepsin2 and pfgch1 genes of P. falciparum. The ddPCR assay is a useful additional tool for the surveillance of anti-malarial drug resistance

Utility of Plasmodium falciparum DNA from rapid diagnostic test kits for molecular analysis and whole genome amplification

Background Rapid diagnostic tests (RDTs) have become the most common diagnostic tool for detection of Plasmodium falciparum malaria, in particular in remote areas. RDT blood spots provide a source of parasite DNA for molecular analysis. In this study, the utility of RDTs for molecular analysis and the performance of different methods for whole genome amplification were investigated. Methods Positive P. falciparum RDTs were collected from Kayin, Myanmar from August 2014 to January 2016. The RDT samples were stored for 6 months, 9 months, 20 months, 21 months, and 32 months before DNA extraction and subsequent molecular analysis of P. falciparum kelch 13 (pfkelch13) mutations, P. falciparum multidrug resistance 1 (pfmdr1), and P. falciparum plasmepsin 2 (pfplasmepsin2) gene amplification. In addition, performance of four whole genome amplification (WGA) kits were compared, including REPLI-g®, MALBACTM, PicoPLEX®, and GenomePlex®, for which DNA quantity and quality were compared between original DNA and post-WGA products. Results The proportion of successful amplification of the different molecular markers was similar between blood spots analysed from RDTs stored for 6, 9, 20, 21, or 32 months. Successful amplification was dependent on the molecular markers fragment length (p value  Conclusions Blood spots from RDTs are a good source for molecular analysis of P. falciparum, even after storage up to 32 months. WGA of RDT-derived parasite DNA reliably increase DNA quantity with sufficient quality for molecular analysis of resistance markers.</p

Molecular characterization of Plasmodium falciparum antifolate resistance markers in Thailand between 2008 and 2016

Background Resistance to anti-malarials is a major threat to the control and elimination of malaria. Sulfadoxine–pyrimethamine (SP) anti-malarial treatment was used as a national policy for treatment of uncomplicated falciparum malaria in Thailand from 1973 to 1990. In order to determine whether withdrawal of this antifolate drug has led to restoration of SP sensitivity, the prevalence of genetic markers of SP resistance was assessed in historical Thai samples. Methods Plasmodium falciparum DNA was collected from the Thailand–Myanmar, Thailand–Malaysia and Thailand–Cambodia borders during 2008–2016 (N = 233). Semi-nested PCR and nucleotide sequencing were used to assess mutations in Plasmodium falciparum dihydrofolate reductase (pfdhfr), P. falciparum dihydropteroate synthase (pfdhps). Gene amplification of Plasmodium falcipaurm GTP cyclohydrolase-1 (pfgch1) was assessed by quantitative real-time PCR. The association between pfdhfr/pfdhps mutations and pfgch1 copy numbers were evaluated. Results Mutations in pfdhfr/pfdhsp and pfgch1 copy number fluctuated overtime through the study period. Altogether, 14 unique pfdhfr–pdfhps haplotypes collectively containing quadruple to octuple mutations were identified. High variation in pfdhfr–pfdhps haplotypes and a high proportion of pfgch1 multiple copy number (51% (73/146)) were observed on the Thailand–Myanmar border compared to other parts of Thailand. Overall, the prevalence of septuple mutations was observed for pfdhfr–pfdhps haplotypes. In particular, the prevalence of pfdhfr–pfdhps, septuple mutation was observed in the Thailand–Myanmar (50%, 73/146) and Thailand–Cambodia (65%, 26/40) border. In Thailand–Malaysia border, majority of the pfdhfr–pfdhps haplotypes transaction from quadruple (90%, 9/10) to quintuple (65%, 24/37) during 2008–2016. Within the pfdhfr–pfdhps haplotypes, during 2008–2013 the pfdhps A/S436F mutation was observed only in Thailand–Myanmar border (9%, 10/107), while it was not identified later. In general, significant correlation was observed between the prevalence of pfdhfr I164L (ϕ = 0.213, p-value = 0.001) or pfdhps K540E/N (ϕ = 0.399, p-value ≤ 0.001) mutations and pfgch1 gene amplification. Conclusions Despite withdrawal of SP as anti-malarial treatment for 17 years, the border regions of Thailand continue to display high prevalence of antifolate and anti-sulfonamide resistance markers in falciparum malaria. Significant association between pfgch1 amplification and pfdhfr (I164L) or pfdhps (K540E) resistance markers were observed, suggesting a compensatory mutation

Molecular characterization of Plasmodium falciparum antifolate resistance markers in Thailand between 2008 and 2016

Background Resistance to anti-malarials is a major threat to the control and elimination of malaria. Sulfadoxine–pyrimethamine (SP) anti-malarial treatment was used as a national policy for treatment of uncomplicated falciparum malaria in Thailand from 1973 to 1990. In order to determine whether withdrawal of this antifolate drug has led to restoration of SP sensitivity, the prevalence of genetic markers of SP resistance was assessed in historical Thai samples. Methods Plasmodium falciparum DNA was collected from the Thailand–Myanmar, Thailand–Malaysia and Thailand–Cambodia borders during 2008–2016 (N = 233). Semi-nested PCR and nucleotide sequencing were used to assess mutations in Plasmodium falciparum dihydrofolate reductase (pfdhfr), P. falciparum dihydropteroate synthase (pfdhps). Gene amplification of Plasmodium falcipaurm GTP cyclohydrolase-1 (pfgch1) was assessed by quantitative real-time PCR. The association between pfdhfr/pfdhps mutations and pfgch1 copy numbers were evaluated. Results Mutations in pfdhfr/pfdhsp and pfgch1 copy number fluctuated overtime through the study period. Altogether, 14 unique pfdhfr–pdfhps haplotypes collectively containing quadruple to octuple mutations were identified. High variation in pfdhfr–pfdhps haplotypes and a high proportion of pfgch1 multiple copy number (51% (73/146)) were observed on the Thailand–Myanmar border compared to other parts of Thailand. Overall, the prevalence of septuple mutations was observed for pfdhfr–pfdhps haplotypes. In particular, the prevalence of pfdhfr–pfdhps, septuple mutation was observed in the Thailand–Myanmar (50%, 73/146) and Thailand–Cambodia (65%, 26/40) border. In Thailand–Malaysia border, majority of the pfdhfr–pfdhps haplotypes transaction from quadruple (90%, 9/10) to quintuple (65%, 24/37) during 2008–2016. Within the pfdhfr–pfdhps haplotypes, during 2008–2013 the pfdhps A/S436F mutation was observed only in Thailand–Myanmar border (9%, 10/107), while it was not identified later. In general, significant correlation was observed between the prevalence of pfdhfr I164L (ϕ = 0.213, p-value = 0.001) or pfdhps K540E/N (ϕ = 0.399, p-value ≤ 0.001) mutations and pfgch1 gene amplification. Conclusions Despite withdrawal of SP as anti-malarial treatment for 17 years, the border regions of Thailand continue to display high prevalence of antifolate and anti-sulfonamide resistance markers in falciparum malaria. Significant association between pfgch1 amplification and pfdhfr (I164L) or pfdhps (K540E) resistance markers were observed, suggesting a compensatory mutation

Protective effect of Mediterranean-type glucose-6-phosphate dehydrogenase deficiency against Plasmodium vivax malaria

X-linked glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common human enzymopathy. The severe Mediterranean variant (G6PD Med) found across Europe and Asia is thought to confer protection against malaria, but its effect is unclear. We fitted a Bayesian statistical model to observed G6PD Med allele frequencies in 999 Pashtun patients presenting with acute Plasmodium vivax malaria and 1408 population controls. G6PD Med was associated with reductions in symptomatic P. vivax malaria incidence of 76% (95% credible interval [CI], 58–88) in hemizygous males and homozygous females combined and 55% (95% CI, 38–68) in heterozygous females. Unless there is very large population stratification within the Pashtun (confounding these results), the G6PD Med genotype confers a very large and gene-dose proportional protective effect against acute vivax malaria. The proportion of patients with vivax malaria at risk of haemolysis following 8-aminoquinoline radical cure is substantially overestimated by studies measuring G6PD deficiency prevalence in healthy subjects

Asymptomatic natural human infections with the Simian malaria parasites plasmodium cynomolgi and plasmodium knowlesi

Background In Southeast Asia, Plasmodium knowlesi, a parasite of long-tailed macaques (Macaca fascicularis), is an important cause of human malaria. Plasmodium cynomolgi also commonly infects these monkeys, but only one naturally acquired symptomatic human case has been reported previously. Methods Malariometric studies involving 5422 subjects (aged 6 months to 65 years) were conducted in 23 villages in Pailin and Battambang, western Cambodia. Parasite detection and genotyping was conducted on blood samples, using high-volume quantitative PCR (uPCR). Results Asymptomatic malaria parasite infections were detected in 1361 of 14732 samples (9.2%). Asymptomatic infections with nonhuman primate malaria parasites were found in 21 individuals living close to forested areas; P. cynomolgi was found in 11, P. knowlesi was found in 8, and P. vivax and P. cynomolgi were both found in 2. Only 2 subjects were female, and 14 were men aged 20–40 years. Geometric mean parasite densities were 3604 parasites/mL in P. cynomolgi infections and 52488 parasites/mL in P. knowlesi infections. All P. cynomolgi isolates had wild-type dihydrofolate reductase genes, in contrast to the very high prevalence of mutations in the human malaria parasites. Asymptomatic reappearance of P. cynomolgi occurred in 2 subjects 3 months after the first infection. Conclusions Asymptomatic naturally acquired P. cynomolgi and P. knowlesi infections can both occur in humans

The spread of artemisinin-resistant Plasmodium falciparum in the Greater Mekong Subregion: a molecular epidemiology observational study

BACKGROUND: Evidence suggests that the PfKelch13 mutations that confer artemisinin resistance in falciparum malaria have multiple independent origins across the Greater Mekong subregion, which has motivated a regional malaria elimination agenda. We aimed to use molecular genotyping to assess antimalarial drug resistance selection and spread in the Greater Mekong subregion. METHODS: In this observational study, we tested Plasmodium falciparum isolates from Myanmar, northeastern Thailand, southern Laos, and western Cambodia for PfKelch13 mutations and for Pfplasmepsin2 gene amplification (indicating piperaquine resistance). We collected blood spots from patients with microscopy or rapid test confirmed uncomplicated falciparum malaria. We used microsatellite genotyping to assess genetic relatedness. FINDINGS: As part of studies on the epidemiology of artemisinin-resistant malaria between Jan 1, 2008, and Dec 31, 2015, we collected 434 isolates. In 2014-15, a single long PfKelch13 C580Y haplotype (-50 to +31·5 kb) lineage, which emerged in western Cambodia in 2008, was detected in 65 of 88 isolates from northeastern Thailand, 86 of 111 isolates from southern Laos, and 14 of 14 isolates from western Cambodia, signifying a hard transnational selective sweep. Pfplasmepsin2 amplification occurred only within this lineage, and by 2015 these closely related parasites were found in ten of the 14 isolates from Cambodia and 15 of 15 isolates from northeastern Thailand. C580Y mutated parasites from Myanmar had a different genetic origin. INTERPRETATION: Our results suggest that the dominant artemisinin-resistant P falciparum C580Y lineage probably arose in western Cambodia and then spread to Thailand and Laos, outcompeting other parasites and acquiring piperaquine resistance. The emergence and spread of fit artemisinin-resistant P falciparum parasite lineages, which then acquire partner drug resistance across the Greater Mekong subregion, threatens regional malaria control and elimination goals. Elimination of falciparum malaria from this region should be accelerated while available antimalarial drugs still remain effective. FUNDING: The Wellcome Trust and the Bill and Melinda Gates Foundation

Molecular epidemiology of resistance to antimalarial drugs in the Greater Mekong subregion: an observational study

Background The Greater Mekong subregion is a recurrent source of antimalarial drug resistance in Plasmodium falciparum malaria. This study aimed to characterise the extent and spread of resistance across this entire region between 2007 and 2018. Methods P falciparum isolates from Myanmar, Thailand, Laos, and Cambodia were obtained from clinical trials and epidemiological studies done between Jan 1, 2007, and Dec 31, 2018, and were genotyped for molecular markers (pfkelch, pfcrt, pfplasmepsin2, and pfmdr1) of antimalarial drug resistance. Genetic relatedness was assessed using microsatellite and single nucleotide polymorphism typing of flanking sequences around target genes. Findings 10 632 isolates were genotyped. A single long pfkelch Cys580Tyr haplotype (from −50 kb to +31·5 kb) conferring artemisinin resistance (PfPailin) now dominates across the eastern Greater Mekong subregion. Piperaquine resistance associated with pfplasmepsin2 gene amplification and mutations in pfcrt downstream of the Lys76Thr chloroquine resistance locus has also developed. On the Thailand–Myanmar border a different pfkelch Cys580Tyr lineage rose to high frequencies before it was eliminated. Elsewhere in Myanmar the Cys580Tyr allele remains widespread at low allele frequencies. Meanwhile a single artemisinin-resistant pfkelch Phe446Ile haplotype has spread across Myanmar. Despite intense use of dihydroartemisinin–piperaquine in Kayin state, eastern Myanmar, both in treatment and mass drug administrations, no selection of piperaquine resistance markers was observed. pfmdr1 amplification, a marker of resistance to mefloquine, remains at low prevalence across the entire region. Interpretation Artemisinin resistance in P falciparum is now prevalent across the Greater Mekong subregion. In the eastern Greater Mekong subregion a multidrug resistant P falciparum lineage (PfPailin) dominates. In Myanmar a long pfkelch Phe446Ile haplotype has spread widely but, by contrast with the eastern Greater Mekong subregion, there is no indication of artemisinin combination therapy (ACT) partner drug resistance from genotyping known markers, and no evidence of spread of ACT resistant P falciparum from the east to the west. There is still a window of opportunity to prevent global spread of ACT resistance