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

Four human Plasmodium species quantification using droplet digital PCR.

Droplet digital polymerase chain reaction (ddPCR) is a partial PCR based on water-oil emulsion droplet technology. It is a highly sensitive method for detecting and delineating minor alleles from complex backgrounds and provides absolute quantification of DNA targets. The ddPCR technology has been applied for detection of many pathogens. Here the sensitive assay utilizing ddPCR for detection and quantification of Plasmodium species was investigated. The assay was developed for two levels of detection, genus specific for all Plasmodium species and for specific Plasmodium species detection. The ddPCR assay was developed based on primers and probes specific to the Plasmodium genus 18S rRNA gene. Using ddPCR for ultra-sensitive P. falciparum assessment, the lower level of detection from concentrated DNA obtained from a high volume (1 mL) blood sample was 11 parasites/mL. For species identification, in particular for samples with mixed infections, a duplex reaction was developed for detection and quantification P. falciparum/ P. vivax and P. malariae/ P. ovale. Amplification of each Plasmodium species in the duplex reaction showed equal sensitivity to singleplex single species detection. The duplex ddPCR assay had higher sensitivity to identify minor species in 32 subpatent parasitaemia samples from Cambodia, and performed better than real-time PCR. The ddPCR assay shows high sensitivity to assess very low parasitaemia of all human Plasmodium species. This provides a useful research tool for studying the role of the asymptomatic parasite reservoir for transmission in regions aiming for malaria elimination

Four human <i>Plasmodium</i> species quantification using droplet digital PCR

<div><p>Droplet digital polymerase chain reaction (ddPCR) is a partial PCR based on water-oil emulsion droplet technology. It is a highly sensitive method for detecting and delineating minor alleles from complex backgrounds and provides absolute quantification of DNA targets. The ddPCR technology has been applied for detection of many pathogens. Here the sensitive assay utilizing ddPCR for detection and quantification of <i>Plasmodium</i> species was investigated. The assay was developed for two levels of detection, genus specific for all <i>Plasmodium</i> species and for specific <i>Plasmodium</i> species detection. The ddPCR assay was developed based on primers and probes specific to the <i>Plasmodium</i> genus <i>18S rRNA</i> gene. Using ddPCR for ultra-sensitive <i>P</i>. <i>falciparum</i> assessment, the lower level of detection from concentrated DNA obtained from a high volume (1 mL) blood sample was 11 parasites/mL. For species identification, in particular for samples with mixed infections, a duplex reaction was developed for detection and quantification <i>P</i>. <i>falciparum</i>/ <i>P</i>. <i>vivax</i> and <i>P</i>. <i>malariae</i>/ <i>P</i>. <i>ovale</i>. Amplification of each <i>Plasmodium</i> species in the duplex reaction showed equal sensitivity to singleplex single species detection. The duplex ddPCR assay had higher sensitivity to identify minor species in 32 subpatent parasitaemia samples from Cambodia, and performed better than real-time PCR. The ddPCR assay shows high sensitivity to assess very low parasitaemia of all human <i>Plasmodium</i> species. This provides a useful research tool for studying the role of the asymptomatic parasite reservoir for transmission in regions aiming for malaria elimination.</p></div

An online mapping database of molecular markers of drug resistance in Plasmodium falciparum: the ACT Partner Drug Molecular Surveyor

Abstract Background Prior to this project, only a handful of online visualizations existed for exploring the published literature on molecular markers of antimalarial drug resistance, and none specifically for the markers associated with Plasmodium falciparum resistance to the partner drugs in artemisinin-based combination therapy (ACT). Molecular information is collected in studies with different designs, using a variety of molecular methodologies and data analysis strategies, making it difficult to compare across studies. The purpose of this project was to develop a free online tool, which visualizes the widely published data on molecular markers of antimalarial drug resistance, starting with the two genes pfcrt and pfmdr-1, associated with resistance to the three most common partner drugs; amodiaquine, lumefantrine and mefloquine. Methods A literature review was conducted, and a standardized method was used to extract data from publications, and critical decisions on visualization were made. A global geospatial database was developed of specific pfmdr1 and pfcrt single nucleotide polymorphisms and pfmdr1 copy number variation. An informatics framework was developed that allowed flexibility in development of the tool over time and efficient adaptation to different source data. Results The database discussed in this paper has pfmdr1 and pfcrt marker prevalence information, from 579 geographic sites in 76 different countries, including results from over 86,000 samples from 456 articles published January 2001–May 2017. The ACT Partner Drugs Molecular Surveyor was launched by the WorldWide Antimalarial Resistance Network (WWARN) in March 2015 and it has attracted over 3000 unique visitors since then. Presented here is a demonstration of how the Surveyor database can be explored to monitor local, temporal changes in the prevalence of molecular markers. Here publications up to May 2017 were included, however the online ACT partner drug Molecular Surveyor is continuously updated with new data and relevant markers. Conclusions The WWARN ACT Partner Drugs Molecular Surveyor summarizes data on resistance markers in the pfmdr1 and pfcrt genes. The database is fully accessible, providing users with a rich resource to explore and analyze, and thus utilize a centralized, standardized database for different purposes. This open-source software framework can be adapted to other data, as demonstrated by the subsequent launch of the Artemisinin Molecular Surveyor and the Vivax Surveyor

Genetic population of Plasmodium knowlesi during pre-malaria elimination in Thailand

Background Thailand is committed to eliminating malaria by 2024. From 2013 to 2020, the total number of malaria cases have decreased, from 37,741 to 4474 (an 88.1% reduction). However, infections with Plasmodium knowlesi, a monkey malarial pathogen that can also infect humans, have been increasingly observed. This study focused on the molecular analysis of P. knowlesi parasites causing malaria in Thailand. Methods Under Thailand’s integrated Drug Efficacy Surveillance (iDES), which includes drug-resistance monitoring as part of routine case-based surveillance and responses, specimens were collected from malaria patients (n = 966) between 2018 and 2020. Thirty-one mono P. knowlesi infections (3.1%), most of which were from eastern and southern Thailand, were observed and confirmed by nested PCR assay and DNA sequencing. To evaluate whether these pathogens were from different lineages, cluster analysis based on seven microsatellite genotyping markers and the merozoite surface protein 1 (pkmsp1) gene was carried out. The P. knowlesi pyrimethamine resistance gene dihydrofolate reductase (pkdhfr) was sequenced and homology modelling was constructed. Results The results of analysing the seven microsatellite markers and pkmsp1 sequence demonstrated that P. knowlesi parasites from eastern Thailand were of the same lineage as those isolated in Cambodia, while the parasites causing malaria in southern Thailand were the same lineage as those isolated from Malaysia. The sequencing results for the pkdhfr genes indicated the presence of two mutations, Arg34Leu and a deletion at position 105. On analysis with homology modelling, the two mutations were not associated with anti-malarial drug resistance. Conclusions This report compared the genetic populations of P. knowlesi parasites in Thailand from 2018 to 2020 and have shown similar lineages as those isolated in Cambodia and Malaysia of P. knowlesi infection in Thailand and demonstrated that the P. knowlesi parasites were of the same lineages as those isolated in Cambodia and Malaysia. The parasites were also shown to be sensitive to pyrimethamine

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

Nanopore sequencing of drug-resistance-associated genes in malaria parasites, Plasmodium falciparum

Here, we report the application of a portable sequencer, MinION, for genotyping the malaria parasite Plasmodium falciparum. In the present study, an amplicon mixture of nine representative genes causing resistance to anti-malaria drugs is diagnosed. First, we developed the procedure for four laboratory strains (3D7, Dd2, 7G8, and K1), and then applied the developed procedure to ten clinical samples. We sequenced and re-sequenced the samples using the obsolete flow cell R7.3 and the most recent flow cell R9.4. Although the average base-call accuracy of the MinION sequencer was 74.3%, performing >50 reads at a given position improves the accuracy of the SNP call, yielding a precision and recall rate of 0.92 and 0.8, respectively, with flow cell R7.3. These numbers increased significantly with flow cell R9.4, in which the precision and recall are 1 and 0.97, respectively. Based on the SNP information, the drug resistance status in ten clinical samples was inferred. We also analyzed K13 gene mutations from 54 additional clinical samples as a proof of concept. We found that a novel amino-acid changing variation is dominant in this area. In addition, we performed a small population-based analysis using 3 and 5 cases (K13) and 10 and 5 cases (PfCRT) from Thailand and Vietnam, respectively. We identified distinct genotypes from the respective regions. This approach will change the standard methodology for the sequencing diagnosis of malaria parasites, especially in developing countries

Evolution of multidrug resistance in Plasmodium falciparum: A longitudinal study of genetic resistance markers in the Greater Mekong subregion

Increasing resistance in Plasmodium falciparum to artemisinins and their artemisinin combination therapy (ACT) partner drugs jeopardizes effective antimalarial treatment. Resistance is worst in the Greater Mekong subregion. Monitoring genetic markers of resistance can help to guide antimalarial therapy. Markers of resistance to artemisinins (PfKelch mutations), mefloquine (amplification of P. falciparum multidrug resistance-1 [PfMDR1]), and piperaquine (PfPlasmepsin2/3 amplification and specific P. falciparum chloroquine resistance transporter [PfCRT] mutations) were assessed in 6,722 P. falciparum samples from Vietnam, Lao People’s Democratic Republic (PDR), Cambodia, Thailand, and Myanmar between 2007 and 2019. Against a high background prevalence of PfKelch mutations, PfMDR1 and PfPlasmepsin2/3 amplification closely followed regional drug pressures over time. PfPlasmepsin2/3 amplification preceded piperaquine resistance-associated PfCRT mutations in Cambodia and reached a peak prevalence of 23/28 (82%) in 2015. This declined to 57/156 (38%) after first-line treatment was changed from dihydroartemisinin-piperaquine to artesunate-mefloquine (ASMQ) between 2014 and 2017. The frequency of PfMDR1 amplification increased from 0/293 (0%) between 2012 and 2017 to 12/156 (8%) in 2019. Amplification of PfMDR1 and PfPlasmepsin2/3 in the same parasites was extremely rare (4/6,722 [0.06%]) and was dispersed over time. The mechanisms conferring mefloquine and piperaquine resistance may be counterbalancing. This supports the development of ASMQ plus piperaquine as a triple artemisinin combination therapy