Multiple Novel Mutations in Plasmodium falciparum Chloroquine Resistance Transporter Gene during Implementation of Artemisinin Combination Therapy in Thailand

Distribution patterns of molecular markers of antimalarial drug resistance in Plasmodium falciparum isolates on the Thai-Myanmar border during the periods of 1993-1998 and 2002-2008

BACKGROUND

Polymorphisms of Plasmodium falciparum chloroquine resistance transporter (pfcrt), Plasmodium falciparum multi-drug resistance 1 (pfmdr1) and Plasmodium falciparum kelch 13-propeller (pfk13) genes are accepted as valid molecular markers of quinoline antimalarials and artemisinins. This study investigated the distribution patterns of these genes in P. falciparum isolates from the areas along the Thai-Myanmar border during the two different periods of antimalarial usage in Thailand.

RESULTS

Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) were used to detect pfcrt mutations at codons 76, 220, 271, 326, 356, and 371 as well as pfmdr1 mutation at codon 86. The prevalence of pfcrt mutations was markedly high (96.4-99.7%) in samples collected during both periods. The proportions of mutant genotypes (number of mutant/total isolate) at codons 76, 220, 271, 326, 356 and 371 in the isolates collected during 1993-1998 (period 1) compared with 2002-2008 (period 2) were 97.9% (137/140) vs. 97.1% (401/413), 97.9% (140/143) vs. 98.8% (171/173), 97.2% (139/143) vs. 97.1% (333/343), 98.6% (140/142) vs. 99.7% (385/386), 96.4% (134/139) vs. 98.2% (378/385) and 97.8% (136/139) vs. 98.9% (375/379), respectively. Most isolates carried pfmdr1 wild-type at codon 86, with a significant difference in proportions genotypes (number of wild type/total sample) in samples collected during period 1 [92.9% (130/140)] compared with period 2 [96.9% (379/391)]. Investigation of pfmdr1 copy number was performed by real-time PCR. The proportions of isolates carried 1, 2, 3 and 4 or more than 4 copies of pfmdr1 (number of isolates carried correspondent copy number/total isolate) were significantly different between the two sample collecting periods (65.7% (90/137) vs. 87.8% (390/444), 18.2% (25/137) vs. 6.3%(28/444), 5.1% (7/137) vs. 1.4% (6/444) and 11.0% (15/137) vs. 4.5% (20/444), for period 1 vs. period 2, respectively). No pfk13 mutation was detected by nested PCR and nucleotide sequencing in all samples with successful analysis (n = 68).

CONCLUSIONS

The persistence of pfcrt mutations and pfmdr1 wild-types at codon 86, along with gene amplification in P. falciparum, contributes to the continued resistance of chloroquine and mefloquine in P. falciparum isolates in the study area. Regular surveillance of antimalarial drug resistance in P. falciparum, incorporating relevant molecular markers and treatment efficacy assessments, should be conducted.

Analysis of sequence diversity in Plasmodium falciparum glutamic acid-rich protein (PfGARP), an asexual blood stage vaccine candidate

Glutamic acid-rich protein of Plasmodium falciparum (PfGARP) binds to erythrocyte band 3 and may enhance cytoadherence of infected erythrocytes. Naturally acquired anti-PfGARP antibodies could confer protection against high parasitemia and severe symptoms. While whole genome sequencing analysis has suggested high conservation in this locus, little is known about repeat polymorphism in this vaccine candidate antigen. Direct sequencing was performed from the PCR-amplified complete PfGARP gene of 80 clinical isolates from four malaria endemic provinces in Thailand and an isolate from a Guinean patient. Publicly available complete coding sequences of this locus were included for comparative analysis. Six complex repeat (RI-RVI) and two homopolymeric glutamic acid repeat (E1 and E2) domains were identified in PfGARP. The erythrocyte band 3-binding ligand in domain RIV and the epitope for mAB7899 antibody eliciting in vitro parasite killing property were perfectly conserved across isolates. Repeat lengths in domains RIII and E1-RVI-E2 seemed to be correlated with parasite density of the patients. Sequence variation in PfGARP exhibited genetic differentiation across most endemic areas of Thailand. Phylogenetic tree inferred from this locus has shown that most Thai isolates formed closely related lineages, suggesting local expansion/contractions of repeat-encoding regions. Positive selection was observed in non-repeat region preceding domain RII which corresponded to a helper T cell epitope predicted to be recognized by a common HLA class II among Thai population. Predicted linear B cell epitopes were identified in both repeat and non-repeat domains. Besides length variation in some repeat domains, sequence conservation in non-repeat regions and almost all predicted immunogenic epitopes have suggested that PfGARP-derived vaccine may largely elicit strain-transcending immunity.

Comparative analysis of targeted next-generation sequencing for Plasmodium falciparum drug resistance markers

Well-defined molecular resistance markers are available for a range of antimalarial drugs, and molecular surveillance is increasingly important for monitoring antimalarial drug resistance. Different genotyping platforms are available, but these have not been compared in detail. We compared Targeted Amplicon Deep sequencing (TADs) using Ion Torrent PGM with Illumina MiSeq for the typing of antimalarial drug resistance genes. We developed and validated protocols to type the molecular resistance markers pfcrt, pfdhfr, pfdhps, pfmdr1, pfkelch, and pfcytochrome b, in Plasmodium falciparum for the Ion Torrent PGM and Illumina MiSeq sequencing platforms. With P. falciparum 3D7 and K1 as reference strains, whole blood samples (N = 20) and blood spots from Rapid Diagnostic Test (RDT) samples (N = 5) from patients with uncomplicated falciparum malaria from Ubon Ratchathani were assessed on both platforms and compared for coverage (average reads per amplicon), sequencing accuracy, variant accuracy, false positive rate, false negative rate, and alternative allele detection, with conventional Sanger sequencing as the reference method for SNP calling. Both whole blood and RDT samples could be successfully sequenced using the Ion Torrent PGM and Illumina MiSeq platforms. Coverage of reads per amplicon was higher with Illumina MiSeq (28,886 reads) than with Ion Torrent PGM (1754 reads). In laboratory generated artificial mixed infections, the two platforms could detect the minor allele down to 1% density at 500X coverage. SNPs calls from both platforms were in complete agreement with conventional Sanger sequencing. The methods can be multiplexed with up to 96 samples per run, which reduces cost by 86% compared to conventional Sanger sequencing. Both platforms, using the developed TAD protocols, provide an accurate method for molecular surveillance of drug resistance markers in P. falciparum, but Illumina MiSeq provides higher coverage than Ion Torrent PGM.

Plasmodium falciparum phenotypic and genotypic resistance profile during the emergence of Piperaquine resistance in Northeastern Thailand

Malaria remains a public health problem in Thailand, especially along its borders where highly mobile populations can contribute to persistent transmission. This study aimed to determine resistant genotypes and phenotypes of 112 Plasmodium falciparum isolates from patients along the Thai-Cambodia border during 2013-2015. The majority of parasites harbored a pfmdr1-Y184F mutation. A single pfmdr1 copy number had CVIET haplotype of amino acids 72-76 of pfcrt and no pfcytb mutations. All isolates had a single pfk13 point mutation (R539T, R539I, or C580Y), and increased % survival in the ring-stage survival assay (except for R539I). Multiple copies of pfpm2 and pfcrt-F145I were detected in 2014 (12.8%) and increased to 30.4% in 2015. Parasites containing either multiple pfpm2 copies with and without pfcrt-F145I or a single pfpm2 copy with pfcrt-F145I exhibited elevated IC₉₀ values of piperaquine. Collectively, the emergence of these resistance patterns in Thailand near Cambodia border mirrored the reports of dihydroartemisinin-piperaquine treatment failures in the adjacent province of Cambodia, Oddar Meanchey, suggesting a migration of parasites across the border. As malaria elimination efforts ramp up in Southeast Asia, host nations militaries and other groups in border regions need to coordinate the proposed interventions.

Polymorphism Analysis of pfmdr1 and pfcrt from Plasmodium falciparum Isolates in Northwestern Nigeria Revealed the Major Markers Associated with Antimalarial Resistance

Suspicion of failure in the effectiveness of artemisinin-based combination therapies (currently the first-line treatment of malaria, worldwide) is leading to the unofficial use of alternative antimalarials, including chloroquine and sulfadoxine/pyrimethamine, across northern Nigeria. To facilitate evidence-based resistance management, antimalarial resistance mutations were investigated in Plasmodium falciparum multidrug resistance-1 (pfmdr1) and chloroquine resistance transporter (pfcrt), in isolates from Kano, northwestern Nigeria. Out of the 88 samples genotyped for pfmdr1N86Y mutation using PCR/restriction fragment length polymorphism, one sample contained the 86Y mutation (86Yfrequency = 1.14%). The analysis of 610 bp fragments of pfmdr1 from 16 isolates revealed two polymorphic sites and low haplotype diversity (Hd = 0.492), with only 86 Y mutations in one isolate, and 184 F replacements in five isolates (184Ffrequency = 31.25%). The analysis of 267 bp fragments of pfcrt isolates revealed high polymorphism (Hd = 0.719), with six haplotypes and seven non-synonymous polymorphic sites. Eleven isolates (61.11%) were chloroquine-resistant, CQR (C72V73I74E75T76 haplotype), two of which had an additional mutation, D57E. An additional sequence was CQR, but of the C72V73M74E75T76 haplotype, while the rest of the sequences (33.33%) were chloroquine susceptible (C72V73M74N75K76 haplotype). The findings of these well characterized resistance markers should be considered when designing resistance management strategies in the northwestern Nigeria.

Transmission of Artemisinin-Resistant Malaria Parasites to Mosquitoes under Antimalarial Drug Pressure

Resistance to artemisinin-based combination therapy (ACT) in the Plasmodium falciparum parasite is threatening to reverse recent gains in reducing global deaths from malaria. While resistance manifests as delayed parasite clearance in patients, the phenotype can only spread geographically via the sexual stages and mosquito transmission. In addition to their asexual killing properties, artemisinin and its derivatives sterilize sexual male gametocytes. Whether resistant parasites overcome this sterilizing effect has not, however, been fully tested.

Resistance to artemisinin-based combination therapy (ACT) in the Plasmodium falciparum parasite is threatening to reverse recent gains in reducing global deaths from malaria. While resistance manifests as delayed parasite clearance in patients, the phenotype can only spread geographically via the sexual stages and mosquito transmission. In addition to their asexual killing properties, artemisinin and its derivatives sterilize sexual male gametocytes. Whether resistant parasites overcome this sterilizing effect has not, however, been fully tested. Here, we analyzed P. falciparum clinical isolates from the Greater Mekong Subregion, each demonstrating delayed clinical clearance and known resistance-associated polymorphisms in the Kelch13 (PfK13^var) gene. As well as demonstrating reduced asexual sensitivity to drug, certain PfK13^var isolates demonstrated a marked reduction in sensitivity to artemisinin in an in vitro male gamete formation assay. Importantly, this same reduction in sensitivity was observed when the most resistant isolate was tested directly in mosquito feeds. These results indicate that, under artemisinin drug pressure, while sensitive parasites are blocked, resistant parasites continue transmission. This selective advantage for resistance transmission could favor acquisition of additional host-specificity or polymorphisms affecting partner drug sensitivity in mixed infections. Favored resistance transmission under ACT coverage could have profound implications for the spread of multidrug-resistant malaria beyond Southeast Asia.