Molecular Markers of Sulfadoxine-Pyrimethamine Resistance in Samples from Children with Uncomplicated Plasmodium falciparum at Three Sites in Angola in 2019.

Sulfadoxine-pyrimethamine (SP) is used for prevention of malaria in pregnant women in Angola. We sequenced the Plasmodium falciparum dihydrofolate reductase (pfdhfr) and dihydropteroate synthase (pfdhps) genes, implicated in SP resistance, in samples collected during a 2019 study of artemisinin-based combination therapy efficacy in Benguela, Lunda Sul, and Zaire provinces. A total of 90 day 0 and day of failure samples were individually sequenced, while 508 day 0 samples from participants without recurrent parasitemia were pooled after DNA extraction into 61 pools. The N51I, C59R, and S108N pfdhfr mutations and A437G pfdhps mutations were present at high proportions in all provinces (weighted allele frequencies, 62% to 100%). The K540E pfdhps mutation was present at lower proportions (10% to 14%). The A581G pfdhps mutation was only observed in Zaire, at a 4.6% estimated prevalence. The I431V and A613S mutations were also only observed in Zaire, at a prevalence of 2.8% to 2.9%. The most common (27% to 66%) reconstructed haplotype in all three provinces was the canonical quadruple pfdhfr pfdhps mutant. The canonical quintuple mutant was absent in Lunda Sul and Benguela and present in 7.9% of samples in Zaire. A single canonical sextuple (2.6%) mutant was observed in Zaire Province. Proportions of the pfdhps K540E and A581G mutations were well below the World Health Organization thresholds for meaningful SP resistance (prevalence of 95% for K540E and 10% for A581G). Samples from therapeutic efficacy studies represent a convenient source of samples for monitoring SP resistance markers.

Portable and cost-effective genetic detection and characterization of Plasmodium falciparum hrp2 using the MinION sequencer.

The prevalence of Plasmodium falciparum hrp2 (pfhrp2)-deleted parasites threatens the efficacy of the most used and sensitive malaria rapid diagnostic tests and highlights the need for continued surveillance for this gene deletion. While PCR methods are adequate for determining pfhrp2 presence or absence, they offer a limited view of its genetic diversity. Here, we present a portable sequencing method using the MinION. Pfhrp2 amplicons were generated from individual samples, barcoded, and pooled for sequencing. To overcome potential crosstalk between barcodes, we implemented a coverage-based threshold for pfhrp2 deletion confirmation. Amino acid repeat types were then counted and visualized with custom Python scripts following de novo assembly. We evaluated this assay using well-characterized reference strains and 152 field isolates with and without pfhrp2 deletions, of which 38 were also sequenced on the PacBio platform to provide a standard for comparison. Of 152 field samples, 93 surpassed the positivity threshold, and of those samples, 62/93 had a dominant pfhrp2 repeat type. PacBio-sequenced samples with a dominant repeat-type profile from the MinION sequencing data matched the PacBio profile. This field-deployable assay can be used alone for surveilling pfhrp2 diversity or as a sequencing-based addition to the World Health Organization's existing deletion surveillance protocol.

Evaluation of a parasite-density based pooled targeted amplicon deep sequencing (TADS) method for molecular surveillance of Plasmodium falciparum drug resistance genes in Haiti.

Sequencing large numbers of individual samples is often needed for countrywide antimalarial drug resistance surveillance. Pooling DNA from several individual samples is an alternative cost and time saving approach for providing allele frequency (AF) estimates at a population level. Using 100 individual patient DNA samples of dried blood spots from a 2017 nationwide drug resistance surveillance study in Haiti, we compared codon coverage of drug resistance-conferring mutations in four Plasmodium falciparum genes (crt, dhps, dhfr, and mdr1), for the same deep sequenced samples run individually and pooled. Samples with similar real-time PCR cycle threshold (Ct) values (+/- 1.0 Ct value) were combined with ten samples per pool. The sequencing success for samples in pools were higher at a lower parasite density than the individual samples sequence method. The median codon coverage for drug resistance-associated mutations in all four genes were greater than 3-fold higher in the pooled samples than in individual samples. The overall codon coverage distribution for pooled samples was wider than the individual samples. The sample pools with < 40 parasites/µL blood showed more discordance in AF calls for dhfr and mdr1 between the individual and pooled samples. This discordance in AF estimation may be due to low amounts of parasite DNA, which could lead to variable PCR amplification efficiencies. Grouping samples with an estimated ≥ 40 parasites/µL blood prior to pooling and deep sequencing yielded the expected population level AF. Pooling DNA samples based on estimates of > 40 parasites/µL prior to deep sequencing can be used for rapid genotyping of a large number of samples for these four genes and possibly other drug resistant markers in population-based studies. As Haiti is a low malaria transmission country with very few mixed infections and continued chloroquine sensitivity, the pooled sequencing approach can be used for routine national molecular surveillance of resistant parasites.

Plasmodium falciparum kelch 13 Mutations, 9 Countries in Africa, 2014-2018.

The spread of drug resistance to antimalarial treatments poses a serious public health risk globally. To combat this risk, molecular surveillance of drug resistance is imperative. We report the prevalence of mutations in the Plasmodium falciparum kelch 13 propeller domain associated with partial artemisinin resistance, which we determined by using Sanger sequencing samples from patients enrolled in therapeutic efficacy studies from 9 sub-Saharan countries during 2014-2018. Of the 2,865 samples successfully sequenced before treatment (day of enrollment) and on the day of treatment failure, 29 (1.0%) samples contained 11 unique nonsynonymous mutations and 83 (2.9%) samples contained 27 unique synonymous mutations. Two samples from Kenya contained the S522C mutation, which has been associated with delayed parasite clearance; however, no samples contained validated or candidate artemisinin-resistance mutations.

Targeted deep amplicon sequencing of antimalarial resistance markers in Plasmodium falciparum isolates from Cameroon.

BACKGROUND: Recent studies showed the first emergence of the R561H artemisinin-associated resistance marker in Africa, which highlights the importance of continued molecular surveillance to assess the selection and spread of this and other drug resistance markers in the region. METHOD: In this study, we used targeted amplicon deep sequencing of 116 isolates collected in two areas of Cameroon to genotype the major drug resistance genes, k13, crt, mdr1, dhfr, and dhps, and the cytochrome b gene (cytb) in Plasmodium falciparum. RESULTS: No confirmed or associated artemisinin resistance markers were observed in Pfk13. In comparison, both major and minor alleles associated with drug resistance were found in Pfcrt, Pfmdr1, Pfdhfr, and Pfdhps. Notably, a high frequency of other nonsynonymous mutations was observed across all the genes, except for Pfcytb, suggesting continued selection pressure. CONCLUSIONS: The results from this study supported the continued use of artemisinin-based combination therapy and administration of sulfadoxine-pyrimethamine for intermittent preventive therapy in pregnant women, and for seasonal chemoprevention in these study sites in Cameroon.

Continued Low Efficacy of Artemether-Lumefantrine in Angola in 2019.

Biennial therapeutic efficacy monitoring is a crucial activity for ensuring the efficacy of currently used artemisinin-based combination therapy in Angola. Children with acute uncomplicated Plasmodium falciparum infection in sentinel sites in the Benguela, Zaire, and Lunda Sul Provinces were treated with artemether-lumefantrine (AL) or artesunate-amodiaquine (ASAQ) and monitored for 28 days to assess clinical and parasitological responses. Molecular correction was performed using seven microsatellite markers. Samples from treatment failures were genotyped for the pfk13, pfcrt, and pfmdr1 genes. Day 3 clearance rates were ≥95% in all arms. Uncorrected day 28 Kaplan-Meier efficacy estimates ranged from 84.2 to 90.1% for the AL arms and 84.7 to 100% for the ASAQ arms. Corrected day 28 estimates were 87.6% (95% confidence interval [CI], 81 to 95%) for the AL arm in Lunda Sul, 92.2% (95% CI, 87 to 98%) for AL in Zaire, 95.6% (95% CI, 91 to 100%) for ASAQ in Zaire, 98.4% (95% CI, 96 to 100%) for AL in Benguela, and 100% for ASAQ in Benguela and Lunda Sul. All 103 analyzed samples had wild-type pfk13 sequences. The 76T pfcrt allele was found in most (92%; 11/12) ASAQ late-failure samples but in only 16% (4/25) of AL failure samples. The N86 pfmdr1 allele was found in 97% (34/35) of treatment failures. The AL efficacy in Lunda Sul was below the 90% World Health Organization threshold, the third time in four rounds that this threshold was crossed for an AL arm in Angola. In contrast, the observed ASAQ efficacy has not been below 95% to date in Angola, including this latest round.

Evaluation of an ensemble-based distance statistic for clustering MLST datasets using epidemiologically defined clusters of cyclosporiasis.

Outbreaks of cyclosporiasis, a food-borne illness caused by the coccidian parasite Cyclospora cayetanensis have increased in the USA in recent years, with approximately 2300 laboratory-confirmed cases reported in 2018. Genotyping tools are needed to inform epidemiological investigations, yet genotyping Cyclospora has proven challenging due to its sexual reproductive cycle which produces complex infections characterized by high genetic heterogeneity. We used targeted amplicon deep sequencing and a recently described ensemble-based distance statistic that accommodates heterogeneous (mixed) genotypes and specimens with partial genotyping data, to genotype and cluster 648 C. cayetanensis samples submitted to CDC in 2018. The performance of the ensemble was assessed by comparing ensemble-identified genetic clusters to analogous clusters identified independently based on common food exposures. Using these epidemiologic clusters as a gold standard, the ensemble facilitated genetic clustering with 93.8% sensitivity and 99.7% specificity. Hence, we anticipate that this procedure will greatly complement epidemiologic investigations of cyclosporiasis.

Malaria Risk and Prevention in Asian Migrants to Angola.

The number of Asian migrants working in sub-Saharan developing countries like Angola has been increasing. Their malaria risk, prevention, and care-seeking practices have not been characterized. A cross-sectional survey was conducted in 733 Chinese and Southeast Asian migrants in Angola. Respondents were interviewed and provided blood samples. Samples were analyzed to detect Plasmodium antigen and characterize host anti-Plasmodium response. Positive samples were genotyped using the pfs47 marker. Most respondents (72%; 95% CI: 68-75) reported using bed nets, but less than 1% reported using chemoprophylaxis. Depending on the assay, 1-4% of respondents had evidence of active malaria infection. By contrast, 55% (95% CI: 52-59) were seropositive for Plasmodium antibodies. Most infections were Plasmodium falciparum, but infection and/or exposure to Plasmodium vivax and Plasmodium malariae was also detected. Seroprevalence by time in Angola showed most exposure occurred locally. One respondent had sufficiently high parasitemia for pfs47 genotyping, which showed that the infection was likely locally acquired despite recent travel to home country. Asian migrants to Angola are at substantial risk of malaria. Employers should consider enhanced malaria prevention programs, including chemoprophylaxis; embassies should encourage prevention practices. Angolan healthcare workers should be aware of high malaria exposure in Asian migrants.

Efficacy and safety of artesunate-amodiaquine and artemether-lumefantrine and prevalence of molecular markers associated with resistance, Guinea: an open-label two-arm randomised controlled trial.

BACKGROUND: Anti-malarial resistance is a threat to recent gains in malaria control. This study aimed to assess the efficacy and safety of artesunate-amodiaquine (ASAQ) and artemether-lumefantrine (AL) in the management of uncomplicated malaria and to measure the prevalence of molecular markers of resistance of Plasmodium falciparum in sentinel sites in Maferinyah and Labé Health Districts in Guinea in 2016. METHODS: This was a two-arm randomised controlled trial of the efficacy of AL and ASAQ among children aged 6-59 months with uncomplicated Plasmodium falciparum malaria in two sites. Children were followed for 28 days to assess clinical and parasitological response. The primary outcome was the Kaplan-Meier estimate of Day 28 (D28) efficacy after correction by microsatellite-genotyping. Pre-treatment (D0) and day of failure samples were assayed for molecular markers of resistance in the pfk13 and pfmdr1 genes. RESULTS: A total of 421 participants were included with 211 participants in the Maferinyah site and 210 in Labé. No early treatment failure was observed in any study arms. However, 22 (5.3%) participants developed a late treatment failure (8 in the ASAQ arm and 14 in the AL arm), which were further classified as 2 recrudescences and 20 reinfections. The Kaplan-Meier estimate of the corrected efficacy at D28 was 100% for both AL and ASAQ in Maferinyah site and 99% (95% Confidence Interval: 97.2-100%) for ASAQ and 99% (97.1-100%) for AL in Labé. The majority of successfully analysed D0 (98%, 380/389) and all day of failure (100%, 22/22) samples were wild type for pfk13. All 9 observed pfk13 mutations were polymorphisms not associated with artemisinin resistance. The NFD haplotype was the predominant haplotype in both D0 (197/362, 54%) and day of failure samples (11/18, 61%) successfully analysed for pfmdr1. CONCLUSION: This study observed high efficacy and safety of both ASAQ and AL in Guinea, providing evidence for their continued use to treat uncomplicated malaria. Continued monitoring of ACT efficacy and safety and molecular makers of resistance in Guinea is important to detect emergence of parasite resistance and to inform evidence-based malaria treatment policies.

Targeted deep amplicon sequencing of kelch 13 and cytochrome b in Plasmodium falciparum isolates from an endemic African country using the Malaria Resistance Surveillance (MaRS) protocol.

BACKGROUND: Routine molecular surveillance for imported drug-resistant malaria parasites to the USA and European Union is an important public health activity. The obtained molecular data are used to help keep chemoprophylaxis and treatment guidelines up to date for persons traveling to malaria endemic countries. Recent advances in next-generation sequencing (NGS) technologies provide a new and effective way of tracking malaria drug-resistant parasites. METHODS: As part of a technology transfer arrangement between the CDC Malaria Branch and the Istituto Superiore di Sanità (ISS), Rome, Italy, the recently described Malaria Resistance Surveillance (MaRS) protocol was used to genotype 148 Plasmodium falciparum isolates from Eritrea for kelch 13 (k13) and cytochrome b (cytb) genes, molecular markers associated with resistance to artemisinin (ART) and atovaquone/proguanil (AP), respectively. RESULTS: Spanning the full-length k13 gene, seven non-synonymous single nucleotide polymorphisms (SNPs) were found (K189N, K189T, E208K, D281V, E401Q, R622I and T535M), of which none have been associated with artemisinin resistance. No mutations were found in cytochrome b. CONCLUSION: All patients successfully genotyped carried parasites susceptible to ART and AP treatment. Future studies between CDC Malaria Branch and ISS are planned to expand the MaRS system, including data sharing, in an effort to maintain up to date treatment guidelines for travelers to malaria endemic countries.

Evolution and Genetic Diversity of the k13 Gene Associated with Artemisinin Delayed Parasite Clearance in Plasmodium falciparum.

Mutations in the Plasmodium falciparumk13 (Pfk13) gene are linked to delayed parasite clearance in response to artemisinin-based combination therapies (ACTs) in Southeast Asia. To explore the evolutionary rate and constraints acting on this gene, k13 orthologs from species sharing a recent common ancestor with P. falciparum and Plasmodium vivax were analyzed. These comparative studies were followed by genetic polymorphism analyses within P. falciparum using 982 complete Pfk13 sequences from public databases and new data obtained by next-generation sequencing from African and Haitian isolates. Although k13 orthologs evolve at heterogeneous rates, the gene was conserved across the genus, with only synonymous substitutions being found at residues where mutations linked to the delayed parasite clearance phenotype have been reported. This suggests that those residues were under constraint from undergoing nonsynonymous changes during evolution of the genus. No fixed nonsynonymous differences were found between Pfk13 and its orthologs in closely related species found in African apes. This indicates that all nonsynonymous substitutions currently found in Pfk13 are younger than the time of divergence between P. falciparum and its closely related species. At the population level, no mutations linked to delayed parasite clearance were found in our samples from Africa and Haiti. However, there is a high number of single Pfk13 mutations segregating in P. falciparum populations, and two predominant alleles are distributed worldwide. This pattern is discussed in terms of how changes in the efficacy of natural selection, affected by population expansion, may have allowed for the emergence of mutations tolerant to ACTs.

Using the Plasmodium mitochondrial genome for classifying mixed-species infections and inferring the geographical origin of P. falciparum parasites imported to the U.S.

The ability to identify mixed-species infections and track the origin of Plasmodium parasites can further enhance the development of treatment and prevention recommendations as well as outbreak investigations. Here, we explore the utility of using the full Plasmodium mitochondrial genome to classify Plasmodium species, detect mixed infections, and infer the geographical origin of imported P. falciparum parasites to the United States (U.S.). Using the recently developed standardized, high-throughput Malaria Resistance Surveillance (MaRS) protocol, the full Plasmodium mitochondrial genomes of 265 malaria cases imported to the U.S. from 2014-2017 were sequenced and analyzed. P. falciparum infections were found in 94.7% (251/265) of samples. Five percent (14/265) of samples were identified as mixed- Plasmodium species or non-P. falciparum, including P. vivax, P. malariae, P. ovale curtisi, and P. ovale wallikeri. P. falciparum mitochondrial haplotypes analysis revealed greater than eighteen percent of samples to have at least two P. falciparum mitochondrial genome haplotypes, indicating either heteroplasmy or multi-clonal infections. Maximum-likelihood phylogenies of 912 P. falciparum mitochondrial genomes with known country origin were used to infer the geographical origin of thirteen samples from persons with unknown travel histories as: Africa (country unspecified) (n = 10), Ghana (n = 1), Southeast Asia (n = 1), and the Philippines (n = 1). We demonstrate the utility and current limitations of using the Plasmodium mitochondrial genome to classify samples with mixed-infections and infer the geographical origin of imported P. falciparum malaria cases to the U.S. with unknown travel history.

Next-Generation Sequencing and Bioinformatics Protocol for Malaria Drug Resistance Marker Surveillance.

The recent advances in next-generation sequencing technologies provide a new and effective way of tracking malaria drug-resistant parasites. To take advantage of this technology, an end-to-end Illumina targeted amplicon deep sequencing (TADS) and bioinformatics pipeline for molecular surveillance of drug resistance in P. falciparum, called malaria resistance surveillance (MaRS), was developed. TADS relies on PCR enriching genomic regions, specifically target genes of interest, prior to deep sequencing. MaRS enables researchers to simultaneously collect data on allele frequencies of multiple full-length P. falciparum drug resistance genes (crt, mdr1, k13, dhfr, dhps, and the cytochrome b gene), as well as the mitochondrial genome. Information is captured at the individual patient level for both known and potential new single nucleotide polymorphisms associated with drug resistance. The MaRS pipeline was validated using 245 imported malaria cases that were reported to the Centers for Disease Control and Prevention (CDC). The chloroquine resistance crt CVIET genotype (mutations underlined) was observed in 42% of samples, the highly pyrimethamine-resistant dhpsIRN triple mutant in 92% of samples, and the sulfadoxine resistance dhps mutation SGEAA in 26% of samples. The mdr1 NFSND genotype was found in 40% of samples. With the exception of two cases imported from Cambodia, no artemisinin resistance k13 alleles were identified, and 99% of patients carried parasites susceptible to atovaquone-proguanil. Our goal is to implement MaRS at the CDC for routine surveillance of imported malaria cases in the United States and to aid in the adoption of this system at participating state public health laboratories, as well as by global partners.

Prevalence of molecular markers of artemisinin and lumefantrine resistance among patients with uncomplicated Plasmodium falciparum malaria in three provinces in Angola, 2015.

BACKGROUND: Artemisinin-based combination therapy is the first-line anti-malarial treatment for uncomplicated Plasmodium falciparum infection in Angola. To date, the prevalence of polymorphisms in the pfk13 gene, associated with artemisinin resistance, and pfmdr1, associated with lumefantrine resistance, have not been systematically studied in Angola. METHODS: DNA was isolated from pretreatment and late treatment failure dried blood spots collected during the 2015 round of therapeutic efficacy studies in Benguela, Lunda Sul, and Zaire Provinces in Angola. The pfk13 propeller domain and pfmdr1 gene were sequenced and analysed for polymorphisms. Pfmdr1 copy number variation was assessed using a real-time PCR method. The association between pfmdr1 and pfk13 mutations and treatment failure was investigated. RESULTS: The majority of pretreatment (99%, 466/469) and all late treatment failure (100%, 50/50) samples were wild type for pfk13. Three of the pretreatment samples (1%) carried the A578S mutation commonly observed in Africa and not associated with artemisinin resistance. All 543 pretreatment and day of late treatment failure samples successfully analysed for pfmdr1 copy number variation carried one copy of pfmdr1. The NYD haplotype was the predominant pfmdr1 haplotype, present in 63% (308/491) of pretreatment samples, followed by NFD, which was present in 32% (157/491) of pretreatment samples. The pfmdr1 N86 allele was overrepresented in day of late treatment failure samples from participants receiving artemether-lumefantrine (p value 0.03). CONCLUSIONS: The pretreatment parasites in patients participating in therapeutic efficacy studies in 2015 in Angola's three sentinel sites showed genetic evidence of susceptibility to artemisinins, consistent with clinical outcome data showing greater than 99% day 3 clearance rates. The lack of increased pfmdr1 copy number is consistent with previous reports from sub-Saharan Africa. Although pfmdr1 NYD and NFD haplotypes were overrepresented in artemether-lumefantrine late treatment failure samples, their role as markers of resistance was unclear given that these haplotypes were also present in the majority of successfully treated patients in the artemether-lumefantrine treatment arms.

Molecular Profile of Malaria Drug Resistance Markers of Plasmodium falciparum in Suriname.

In Suriname, an artesunate monotherapy therapeutic efficacy trial was recently conducted to evaluate partial artemisinin resistance emerging in Plasmodium falciparum We genotyped the PfK13 propeller domain of P. falciparum in 40 samples as well as other mutations proposed to be associated with artemisinin-resistant mutants. We did not find any mutations previously associated with artemisinin resistance in Southeast Asia, but we found fixed resistance mutations for chloroquine (CQ) and sulfadoxine-pyrimethamine. Additionally, the PfCRT C350R mutation, associated with reversal of CQ resistance and piperaquine-selective pressure, was present in 62% of the samples. Our results from neutral microsatellite data also confirmed a high parasite gene flow in the Guiana Shield. Although recruiting participants for therapeutic efficacy studies is challenging in areas where malaria endemicity is very low due to the low number of malaria cases reported, conducting these studies along with molecular surveillance remains essential for the monitoring of artemisinin-resistant alleles and for the characterization of the population structure of P. falciparum in areas targeted for malaria elimination.

Efficacy of artemether-lumefantrine, artesunate-amodiaquine, and dihydroartemisinin-piperaquine for treatment of uncomplicated Plasmodium falciparum malaria in Angola, 2015.

BACKGROUND: Recent anti-malarial resistance monitoring in Angola has shown efficacy of artemether-lumefantrine (AL) in certain sites approaching the key 90% lower limit of efficacy recommended for artemisinin-based combination therapy. In addition, a controversial case of malaria unresponsive to artemisinins was reported in a patient infected in Lunda Sul Province in 2013. METHODS: During January-June 2015, investigators monitored the clinical and parasitological response of children with uncomplicated Plasmodium falciparum infection treated with AL, artesunate-amodiaquine (ASAQ), or dihydroartemisinin-piperaquine (DP). The study comprised two treatment arms in each of three provinces: Benguela (AL, ASAQ), Zaire (AL, DP), and Lunda Sul (ASAQ, DP). Samples from treatment failures were analysed for molecular markers of resistance for artemisinin (K13) and lumefantrine (pfmdr1). RESULTS: A total of 467 children reached a study endpoint. Fifty-four treatment failures were observed: four early treatment failures, 40 re-infections and ten recrudescences. Excluding re-infections, the 28-day microsatellite-corrected efficacy was 96.3% (95% CI 91-100) for AL in Benguela, 99.9% (95-100) for ASAQ in Benguela, 88.1% (81-95) for AL in Zaire, and 100% for ASAQ in Lunda Sul. For DP, the 42-day corrected efficacy was 98.8% (96-100) in Zaire and 100% in Lunda Sul. All treatment failures were wild type for K13, but all AL treatment failures had pfmdr1 haplotypes associated with decreased lumefantrine susceptibility. CONCLUSIONS: No evidence was found to corroborate the specific allegation of artemisinin resistance in Lunda Sul. The efficacy below 90% of AL in Zaire matches findings from 2013 from the same site. Further monitoring, particularly including measurement of lumefantrine blood levels, is recommended.