Pooled Deep Sequencing of Drug Resistance Loci from Plasmodium falciparum Parasites across Ethiopia

Application of a new highly multiplexed amplicon sequencing tool to evaluate Plasmodium falciparum antimalarial resistance and relatedness in individual and pooled samples from Dschang, Cameroon

Background

Resistance to antimalarial drugs remains a major obstacle to malaria elimination. Multiplexed, targeted amplicon sequencing is being adopted for surveilling resistance and dissecting the genetics of complex malaria infections. Moreover, genotyping of parasites and detection of molecular markers drug resistance in resource-limited regions requires open-source protocols for processing samples, using accessible reagents, and rapid methods for processing numerous samples including pooled sequencing.

Methods

P lasmodium falciparum S treamlined M ultiplex A ntimalarial R esistance and R elatedness T esting ( Pf -SMARRT) is a PCR-based amplicon panel consisting of 15 amplicons targeting antimalarial resistance mutations and 9 amplicons targeting hypervariable regions. This assay uses oligonucleotide primers in two pools and a non-proprietary library and barcoding approach.

Results

We evaluated Pf -SMARRT using control mocked dried blood spots (DBS) at varying levels of parasitemia and a mixture of 3D7 and Dd2 strains at known frequencies, showing the ability to genotype at low parasite density and recall within-sample allele frequencies. We then piloted Pf -SMARRT to genotype 100 parasite isolates collected from uncomplicated malaria cases at three health facilities in Dschang, Western Cameroon. Antimalarial resistance genotyping showed high levels of sulfadoxine-pyrimethamine resistance mutations, including 31% prevalence of the DHPS A613S mutation. No K13 candidate or validated artemisinin partial resistance mutations were detected, but one low-level non-synonymous change was observed. Pf -SMARRT's hypervariable targets, used to assess complexity of infections and parasite diversity and relatedness, showed similar levels and patterns compared to molecular inversion probe (MIP) sequencing. While there was strong concordance of antimalarial resistance mutations between individual samples and pools, low-frequency variants in the pooled samples were often missed.

Conclusion

Overall, Pf -SMARRT is a robust tool for assessing parasite relatedness and antimalarial drug resistance markers from both individual and pooled samples. Control samples support that accurate genotyping as low as 1 parasite per microliter is routinely possible.

SCOPE STATEMENT 200

Malaria remains a critical global public health problem. Antimalarial drug resistance has repeatedly undermined control and the emergence of artemisinin partial resistance in Africa is the latest major challenge. Malaria molecular surveillance (MMS) has emerged as a powerful tool to monitor molecular markers of resistance and changes in the parasite population. Streamlined methods are needed that can be readily adopted in endemic countries. We developed P lasmodium falciparum S treamlined M ultiplex A ntimalarial R esistance and R elatedness T esting ( Pf -SMARRT), a multiplex amplicon deep sequencing approach that uses easily accessible products without proprietary steps and can be sequenced on any Illumina sequencer. We validated this tool using controls, including mocked dried blood spots, and then implemented it to evaluate resistance and parasite relatedness among 100 samples from Cameroon. The assay was able to reliably assess the within-sample allele frequency of antimalarial resistance markers and discriminate strains within and between individuals. We also evaluated a more cost-effective surveillance approach for antimalarial resistance polymorphisms using pooled samples. While within-pool frequencies of mutations were accurate in pools with higher numbers of samples, this resulted in the loss of the ability to detect variants uncommon in the pool. Overall Pf- SMARRT provides a new protocol for conducting MMS that is easily implementable in Africa.

High frequency of artemisinin partial resistance mutations in the great lake region revealed through rapid pooled deep sequencing

In Africa, the first Plasmodium falciparum Kelch13 (K13) artemisinin partial resistance mutation 561H was first detected and validated in Rwanda. Surveillance to better define the extent of the emergence in Rwanda and neighboring countries as other mutations arise in East Africa is critical. We employ a novel scheme of liquid blood drop preservation combined with pooled sequencing to provide a cost-effective rapid assessment of resistance mutation frequencies at multiple collection sites across Rwanda and neighboring countries. Malaria-positive samples (n=5,465) were collected from 39 health facilities in Rwanda, Uganda, Tanzania, and the Democratic Republic of the Congo (DRC) between May 2022 and March 2023 and sequenced in 199 pools. In Rwanda, K13 561H and 675V were detected in 90% and 65% of sites with an average frequency of 19.0% (0-54.5%) and 5.0% (0-35.5%), respectively. In Tanzania, 561H had high frequency in multiple sites while it was absent from the DRC although 675V was seen at low frequency. Conceringly candidate mutations were observed: 441L, 449A, and 469F co-occurred with validated mutations suggesting they are arising under the same pressures. Other resistance markers associated with artemether-lumefantrine are common: P. falciparum multidrug resistance protein 1 N86 at 98.0% and 184F at 47.0% (0-94.3%) and P. falciparum chloroquine resistance transporter 76T at 14.7% (0-58.6%). Additionally, sulfadoxine-pyrimethamine-associated mutations show high frequencies. Overall, K13 mutations are rapidly expanding in the region further endangering control efforts with the potential of engendering partner drug resistance.

Plasmodium falciparum transmission in the highlands of Ethiopia is driven by closely related and clonal parasites

Malaria cases are frequently recorded in the Ethiopian highlands even at altitudes above 2000 m. The epidemiology of malaria in the Ethiopian highlands, and, in particular, the role of importation by human migration from the highly endemic lowlands is not well understood. We sequenced 187 Plasmodium falciparum samples from two sites in the Ethiopian highlands, Gondar (n = 159) and Ziway (n = 28), using a multiplexed droplet digital PCR (ddPCR)-based amplicon sequencing method targeting 35 microhaplotypes and drug resistance loci. Here, we characterize the parasite population structure and genetic relatedness. We identify moderate parasite diversity (mean HE : 0.54) and low infection complexity (74.9% monoclonal). A significant percentage of infections share microhaplotypes, even across transmission seasons and sites, indicating persistent local transmission. We identify multiple clusters of clonal or near-clonal infections, highlighting high genetic relatedness. Only 6.3% of individuals diagnosed with P. falciparum reported recent travel. Yet, in clonal or near-clonal clusters, infections of travellers were frequently observed first in time, suggesting that parasites may have been imported and then transmitted locally. 31.1% of infections are pfhrp2-deleted and 84.4% pfhrp3-deleted, and 28.7% have pfhrp2/3 double deletions. Parasites with pfhrp2/3 deletions and wild-type parasites are genetically distinct. Mutations associated with resistance to sulphadoxine-pyrimethamine or suggested to reduce sensitivity to lumefantrine are observed at near-fixation. In conclusion, genomic data corroborate local transmission and the importance of intensified control in the Ethiopian highlands.

Molecular surveillance of Kelch 13 polymorphisms in Plasmodium falciparum isolates from Kenya and Ethiopia

BACKGROUND

Timely molecular surveillance of Plasmodium falciparum kelch 13 (k13) gene mutations is essential for monitoring the emergence and stemming the spread of artemisinin resistance. Widespread artemisinin resistance, as observed in Southeast Asia, would reverse significant gains that have been made against the malaria burden in Africa. The purpose of this study was to assess the prevalence of k13 polymorphisms in western Kenya and Ethiopia at sites representing varying transmission intensities between 2018 and 2022.

METHODS

Dried blood spot samples collected through ongoing passive surveillance and malaria epidemiological studies, respectively, were investigated. The k13 gene was genotyped in P. falciparum isolates with high parasitaemia: 775 isolates from four sites in western Kenya (Homa Bay, Kakamega, Kisii, and Kombewa) and 319 isolates from five sites across Ethiopia (Arjo, Awash, Gambella, Dire Dawa, and Semera). DNA sequence variation and neutrality were analysed within each study site where mutant alleles were detected.

RESULTS

Sixteen Kelch13 haplotypes were detected in this study. Prevalence of nonsynonymous k13 mutations was low in both western Kenya (25/783, 3.19%) and Ethiopia (5/319, 1.57%) across the study period. Two WHO-validated mutations were detected: A675V in three isolates from Kenya and R622I in four isolates from Ethiopia. Seventeen samples from Kenya carried synonymous mutations (2.17%). No synonymous mutations were detected in Ethiopia. Genetic variation analyses and tests of neutrality further suggest an excess of low frequency polymorphisms in each study site. Fu and Li's F test statistic in Semera was 0.48 (P > 0.05), suggesting potential population selection of R622I, which appeared at a relatively high frequency (3/22, 13.04%).

CONCLUSIONS

This study presents an updated report on the low frequency of k13 mutations in western Kenya and Ethiopia. The WHO-validated R622I mutation, which has previously only been reported along the north-west border of Ethiopia, appeared in four isolates collected from eastern Ethiopia. The rapid expansion of R622I across Ethiopia signals the need for enhanced monitoring of the spread of drug-resistant P. falciparum parasites in East Africa. Although ACT remains currently efficacious in the study areas, continued surveillance is necessary to detect early indicators of artemisinin partial resistance.

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.

Unraveling the Complexity of Imported Malaria Infections by Amplicon Deep Sequencing

Imported malaria and recurrent infections are becoming an emerging issue in many malaria non-endemic countries. This study aimed to determine the molecular patterns of the imported malaria infections and recurrence. Blood samples were collected from patients with imported malaria infections during 2016-2018 in Guangxi Zhuang Autonomous Region, China. Next-generation amplicon deep-sequencing approaches were used to assess parasite genetic diversity, multiplexity of infection, relapse, recrudescence, and antimalarial drug resistance. A total of 44 imported malaria cases were examined during the study, of which 35 (79.5%) had recurrent malaria infections within 1 year. The majority (91.4%) had one recurrent malaria episode, whereas two patients had two recurrences and one patient had three recurrences. A total of 19 recurrence patterns (the species responsible for primary and successive clinical episodes) were found in patients returning from malaria epidemic countries. Four parasite species were detected with a higher than usual proportion (46.2%) of non-falciparum infections or mixed-species infections. An increasing trend of recurrence infections and reduced drug treatment efficacy were observed among the cases of imported malaria. The high recurrence rate and complex patterns of imported malaria from Africa to non-endemic countries have the potential to initiate local transmission, thereby undermining efforts to eliminate locally acquired malaria. Our findings highlight the power of amplicon deep-sequencing applications in molecular epidemiological studies of the imported malaria recurrences.

Antimalarial Drug Resistance and Implications for the WHO Global Technical Strategy

PURPOSE OF REVIEW

Five years have passed since the World Health Organization released its Global Technical Strategy for Malaria (GTS). In that time, progress against malaria has plateaued. This review focuses on the implications of antimalarial drug resistance for the GTS and how interim progress in parasite genomics and antimalarial pharmacology offer a bulwark against it.

RECENT FINDINGS

For the first time, drug resistance-conferring genes have been identified and validated before their global expansion in malaria parasite populations. More efficient methods for their detection and elaboration have been developed, although low-density infections and polyclonality remain a nuisance to be solved. Clinical trials of alternative regimens for multidrug-resistant malaria have delivered promising results. New agents continue down the development pipeline, while a nascent infrastructure in sub-Saharan Africa for conducting phase I trials and trials of transmission-blocking agents has come to fruition after years of preparation.

SUMMARY

These and other developments can help inform the GTS as the world looks ahead to the next two decades of its implementation. To remain ahead of the threat that drug resistance poses, wider application of genomic-based surveillance and optimization of existing and forthcoming antimalarial drugs are essential.