Describing the current status of Plasmodium falciparum population structure and drug resistance within mainland Tanzania using molecular inversion probes

Prevalence of mutations associated with artemisinin partial resistance and sulfadoxine-pyrimethamine resistance in 13 regions in Tanzania in 2021: a cross-sectional survey

BACKGROUND

The emergence of the artemisinin partial resistance (ART-R) mutation in the Plasmodium falciparum kelch13 gene (k13), Arg561His, in Rwanda and the regional presence of polymorphisms affecting sulfadoxine-pyrimethamine have raised concern in neighbouring Tanzania. The goal of this study was to assess the status of antimalarial resistance in Tanzania, with a focus on the border with Rwanda, to understand the distribution of the Arg561His mutation, partner drug resistance, and resistance to chemoprevention drugs.

METHODS

In this cross-sectional survey, capillary dried blood spots were collected from malaria positive asymptomatic individuals in the community and symptomatic individuals in health facilities aged 6 months and older, in 13 regions of mainland Tanzania from Jan 31 to June 26, 2021. Exclusion criteria included residence of the areas other than the target sites, presenting to the health facility for care and treatment of conditions other than malaria, and not providing informed consent. Samples were assessed for antimalarial resistance polymorphisms and genetic relatedness using molecular inversion probes targeting P falciparum and short-read whole-genome sequencing. The primary outcome was the prevalence of molecular markers of antimalarial resistance at the region level, as well as at the district level in Kagera, a region in the northwest of the country at the border with Rwanda.

FINDINGS

6855 (88·1%) of 7782 capillary dried blood spot samples collected were successfully genotyped. The overall prevalence of k13 Arg561His in Kagera was 7·7% (90% CI 6·0-9·4; 50 of 649), with the highest prevalence in the districts near the Rwandan border (22·8% [31 of 136] in Karagwe, 14·4% [17 of 118]) in Kyerwa, and 1·4% [two of 144] in Ngara). k13 Arg561His was uncommon in the other regions. Haplotype analysis suggested that some of these parasites are related to isolates collected in Rwanda in 2015, supporting regional spread of Arg561His. However, a novel k13 Arg561His haplotype was observed, potentially indicating a second origin in the region. Other validated k13 resistance polymorphisms (one Arg622Ile and two Ala675Val isolates) were also identified. A region of prevalent dihydrofolate reductase Ile164Leu mutation, associated with sulfadoxine-pyrimethamine resistance, was also identified in Kagera (15·2% [12·6-17·8%]; 80 of 526). The mutant crt Lys76Thr mutation, associated with chloroquine and amodiaquine resistance, was uncommon, occurring only in 75 of 2861 genotyped isolates, whereases the wild-type mdr1 Asn86Tyr allele, associated with reduced sensitivity to lumefantrine, was found in 99·7% (3819 of 3830) of samples countrywide.

INTERPRETATION

These findings show that the k13 Arg561His mutation is common in northwest Tanzania and that multiple emergences of ART-R, similar as to what was seen in southeast Asia, have occurred. Mutations associated with high levels of sulfadoxine-pyrimethamine resistance are common. These results raise concerns about the long-term efficacy of artemisinin and chemoprevention antimalarials in the region. Understanding how multiple emergences interact with drivers of regional spread is essential for combating ART-R in Africa.

FUNDING

This study was funded by the Bill & Melinda Gates Foundation and the National Institutes of Health.

The emergence of artemisinin partial resistance in Africa: how do we respond?

Malaria remains one of the most important infectious diseases in the world, with the greatest burden in sub-Saharan Africa, primarily from Plasmodium falciparum infection. The treatment and control of malaria is challenged by resistance to most available drugs, but partial resistance to artemisinins (ART-R), the most important class for the treatment of malaria, was until recently confined to southeast Asia. This situation has changed, with the emergence of ART-R in multiple countries in eastern Africa. ART-R is mediated primarily by single point mutations in the P falciparum kelch13 protein, with several mutations present in African parasites that are now validated resistance mediators based on clinical and laboratory criteria. Major priorities at present are the expansion of genomic surveillance for ART-R mutations across the continent, more frequent testing of the efficacies of artemisinin-based regimens against uncomplicated and severe malaria in trials, more regular assessment of ex-vivo antimalarial drug susceptibilities, consideration of changes in treatment policy to deter the spread of ART-R, and accelerated development of new antimalarial regimens to overcome the impacts of ART-R. The emergence of ART-R in Africa is an urgent concern, and it is essential that we increase efforts to characterise its spread and mitigate its impact.

Phase I study on the pharmacokinetics of intravaginal, self-administered artesunate vaginal pessaries among women in Kenya

Cervical cancer remains a significant global health issue, especially in low- and middle-income countries (LMICs), where access to prevention and treatment is limited and women are at a higher risk of cervical cancer. Artesunate, a widely available drug used to treat malaria, has shown promise in treating human papillomavirus (HPV)-associated anogenital lesions including high-grade cervical precancer, in a recent Phase I studies in the United States. Data on the pharmacokinetics of artesunate following intravaginal use, and its implications on malaria resistance, are lacking.

Objectives

The primary objective of this study is to investigate the pharmacokinetics of Artesunate (AS) and its active metabolite, dihydroartemisinin (DHA) following intravaginal use at the dosing and frequency intended for cervical precancer treatment. A secondary objective is to assess safety among study participants.

Methods

We are conducting a single-arm, phase I trial with a sample size of 12 female volunteers. Participants will self-administer artesunate vaginal pessaries in the study clinic daily for 5 consecutive days. Participants will have their blood drawn prior to receiving the first dose of artesunate on day one of the study and then will receive 8 blood draws on study day five, prior to artesunate administration and at 15 minutes, 30 minutes, 1 hour, 2 hours, 4 hours, 6 hours, and 8 hours after pessary administration. Pharmacokinetic parameters of artesunate and DHA will be calculated by way of quantitative analysis of with determination of maximum concentration (Cmax), time to Cmax (Tmax), area under the plasma concentration versus time curve (AUC), apparent clearance, and elimination half-life (t1/2).

Plasmodium falciparum population dynamics in East Africa and genomic surveillance along the Kenya-Uganda border

East African countries accounted for ~ 10% of all malaria prevalence worldwide in 2022, with an estimated 23.8 million cases and > 53,000 deaths. Despite recent increases in malaria incidence, high-resolution genome-wide analyses of Plasmodium parasite populations are sparse in Kenya, Tanzania, and Uganda. The Kenyan-Ugandan border region is a particular concern, with Uganda confirming the emergence and spread of artemisinin resistant P. falciparum parasites. To establish genomic surveillance along the Kenyan-Ugandan border and analyse P. falciparum population dynamics within East Africa, we generated whole-genome sequencing (WGS) data for 38 parasites from Bungoma, Western Kenya. These sequences were integrated into a genomic analysis of available East African isolate data (n = 599) and revealed parasite subpopulations with distinct genetic structure and diverse ancestral origins. Ancestral admixture analysis of these subpopulations alongside isolates from across Africa (n = 365) suggested potential independent ancestral populations from other major African populations. Within isolates from Western Kenya, the prevalence of biomarkers associated with chloroquine resistance (e.g. Pfcrt K76T) were significantly reduced compared to wider East African populations and a single isolate contained the PfK13 V568I variant, potentially linked to reduced susceptibility to artemisinin. Overall, our work provides baseline WGS data and analysis for future malaria genomic surveillance in the region.

High-Throughput Genotyping of Plasmodium vivax in the Peruvian Amazon via Molecular Inversion Probes

Plasmodium vivax transmission occurs throughout the tropics and is an emerging threat in areas of Plasmodium falciparum decline, causing relapse infections that complicate treatment and control. Targeted sequencing for P. falciparum has been widely deployed to detect population structure and the geographic spread of antimalarial and diagnostic resistance. However, there are fewer such tools for P. vivax . Leveraging global variation data, we designed four molecular inversion probe (MIP) genotyping panels targeting geographically differentiating SNPs, neutral SNPs, putative antimalarial resistance genes, and vaccine candidate genes. We deployed these MIP panels on 866 infections from the Peruvian Amazon and identified transmission networks with clonality (IBD>0.99), copy number variation in Pvdbp and multiple Pvrbps , fixation of putative antimalarial resistance, and balancing selection in 13 vaccine candidate genes. Our MIP panels are the broadest genotyping panel currently available and are poised for successful deployment in other regions of P. vivax transmission.

Strong isolation by distance and evidence of population microstructure reflect ongoing Plasmodium falciparum transmission in Zanzibar

Background

The Zanzibar archipelago of Tanzania has become a low-transmission area for Plasmodium falciparum. Despite being considered an area of pre-elimination for years, achieving elimination has been difficult, likely due to a combination of imported infections from mainland Tanzania and continued local transmission.

Methods

To shed light on these sources of transmission, we applied highly multiplexed genotyping utilizing molecular inversion probes to characterize the genetic relatedness of 282 P. falciparum isolates collected across Zanzibar and in Bagamoyo district on the coastal mainland from 2016 to 2018.

Results

Overall, parasite populations on the coastal mainland and Zanzibar archipelago remain highly related. However, parasite isolates from Zanzibar exhibit population microstructure due to the rapid decay of parasite relatedness over very short distances. This, along with highly related pairs within shehias, suggests ongoing low-level local transmission. We also identified highly related parasites across shehias that reflect human mobility on the main island of Unguja and identified a cluster of highly related parasites, suggestive of an outbreak, in the Micheweni district on Pemba island. Parasites in asymptomatic infections demonstrated higher complexity of infection than those in symptomatic infections, but have similar core genomes.

Conclusions

Our data support importation as a main source of genetic diversity and contribution to the parasite population in Zanzibar, but they also show local outbreak clusters where targeted interventions are essential to block local transmission. These results highlight the need for preventive measures against imported malaria and enhanced control measures in areas that remain receptive to malaria reemergence due to susceptible hosts and competent vectors.

Funding

This research was funded by the National Institutes of Health, grants R01AI121558, R01AI137395, R01AI155730, F30AI143172, and K24AI134990. Funding was also contributed from the Swedish Research Council, Erling-Persson Family Foundation, and the Yang Fund. RV acknowledges funding from the MRC Centre for Global Infectious Disease Analysis (reference MR/R015600/1), jointly funded by the UK Medical Research Council (MRC) and the UK Foreign, Commonwealth & Development Office (FCDO), under the MRC/FCDO Concordat agreement and is also part of the EDCTP2 program supported by the European Union. RV also acknowledges funding by Community Jameel.

Therapeutic efficacy and safety of artemether-lumefantrine for uncomplicated Plasmodium falciparum malaria treatment in Metehara, Central-east Ethiopia

BACKGROUND

Malaria remains a major global health problem although there was a remarkable achievement between 2000 and 2015. Malaria drug resistance, along with several other factors, presents a significant challenge to malaria control and elimination efforts. Numerous countries in sub-Saharan Africa have documented the presence of confirmed or potential markers of partial resistance against artemisinin, the drug of choice for the treatment of uncomplicated Plasmodium falciparum malaria. The World Health Organization (WHO) recommends regular surveillance of artemisinin therapeutic efficacy to inform policy decisions.

METHODS

This study aimed to evaluate the therapeutic efficacy of artemether-lumefantrine (AL), which is the first-line treatment for uncomplicated P. falciparum malaria in Ethiopia since 2004. Using a single-arm prospective evaluation design, the study assessed the clinical and parasitological responses of patients with uncomplicated P. falciparum malaria in Metehara Health Centre, central-east Ethiopia. Out of 2332 malaria suspects (1187 males, 1145 females) screened, 80 (50 males, 30 females) were enrolled, followed up for 28 days, and 73 (44 males, 29 females) completed the follow up. The study was conducted and data was analysed by employing the per-protocol and Kaplan-Meier analyses following the WHO Malaria Therapeutic Efficacy Evaluation Guidelines 2009.

RESULTS

The results indicated rapid parasite clearance and resolution of clinical symptoms, with all patients achieving complete recovery from asexual parasitaemia and fever by day (D) 3. The prevalence of gametocytes decreased from 6.3% on D0 to 2.5% on D2, D3, D7, and ultimately achieving complete clearance afterward.

CONCLUSION

The overall cure rate for AL treatment was 100%, demonstrating its high efficacy in effectively eliminating malaria parasites in patients. No serious adverse events related to AL treatment were reported during the study, suggesting its safety and tolerability among the participants. These findings confirm that AL remains a highly efficacious treatment for uncomplicated P. falciparum malaria in the study site after 20 years of its introduction in Ethiopia.

Temporal genomic analysis of Plasmodium falciparum reveals increased prevalence of mutations associated with delayed clearance following treatment with artemisinin-lumefantrine in Choma District, Southern Province, Zambia

The emergence of antimalarial drug resistance is an impediment to malaria control and elimination in Africa. Analysis of temporal trends in molecular markers of resistance is critical to inform policy makers and guide malaria treatment guidelines. In a low and seasonal transmission region of southern Zambia, we successfully genotyped 85.5% (389/455) of Plasmodium falciparum samples collected between 2013-2018 from 8 spatially clustered health centres using molecular inversion probes (MIPs) targeting key drug resistance genes. Aside from one sample carrying K13 R622I, none of the isolates carried other World Health Organization-validated or candidate artemisinin partial resistance (ART-R) mutations in K13. However, 13% (CI, 9.6-17.2) of isolates had the AP2MU S160N mutation, which has been associated with delayed clearance following artemisinin combination therapy in Africa. This mutation increased in prevalence between 2015-2018 and bears a genomic signature of selection. During this time period, there was an increase in the MDR1 NFD haplotype that is associated with reduced susceptibility to lumefantrine. Sulfadoxine-pyrimethamine polymorphisms were near fixation. While validated ART-R mutations are rare, a mutation associated with slow parasite clearance in Africa appears to be under selection in southern Zambia.

Emergence, transmission dynamics and mechanisms of artemisinin partial resistance in malaria parasites in Africa

Malaria, mostly due to Plasmodium falciparum infection in Africa, remains one of the most important infectious diseases in the world. Standard treatment for uncomplicated P. falciparum malaria is artemisinin-based combination therapy (ACT), which includes a rapid-acting artemisinin derivative plus a longer-acting partner drug, and standard therapy for severe P. falciparum malaria is intravenous artesunate. The efficacy of artemisinins and ACT has been threatened by the emergence of artemisinin partial resistance in Southeast Asia, mediated principally by mutations in the P. falciparum Kelch 13 (K13) protein. High ACT treatment failure rates have occurred when resistance to partner drugs is also seen. Recently, artemisinin partial resistance has emerged in Rwanda, Uganda and the Horn of Africa, with independent emergences of different K13 mutants in each region. In this Review, we summarize our current knowledge of artemisinin partial resistance and focus on the emergence of resistance in Africa, including its epidemiology, transmission dynamics and mechanisms. At present, the clinical impact of emerging resistance in Africa is unclear and most available evidence suggests that the efficacies of leading ACTs remain excellent, but there is an urgent need to better appreciate the extent of the problem and its consequences for the treatment and control of malaria.

Persistent and multiclonal malaria parasite dynamics despite extended artemether-lumefantrine treatment in children

Standard diagnostics used in longitudinal antimalarial studies are unable to characterize the complexity of submicroscopic parasite dynamics, particularly in high transmission settings. We use molecular markers and amplicon sequencing to characterize post-treatment stage-specific malaria parasite dynamics during a 42 day randomized trial of 3- versus 5 day artemether-lumefantrine in 303 children with and without HIV (ClinicalTrials.gov number NCT03453840). The prevalence of parasite-derived 18S rRNA is >70% in children throughout follow-up, and the ring-stage marker SBP1 is detectable in over 15% of children on day 14 despite effective treatment. We find that the extended regimen significantly lowers the risk of recurrent ring-stage parasitemia compared to the standard 3 day regimen, and that higher day 7 lumefantrine concentrations decrease the probability of ring-stage parasites in the early post-treatment period. Longitudinal amplicon sequencing reveals remarkably dynamic patterns of multiclonal infections that include new and persistent clones in both the early post-treatment and later time periods. Our data indicate that post-treatment parasite dynamics are highly complex despite efficacious therapy, findings that will inform strategies to optimize regimens in the face of emerging partial artemisinin resistance in Africa.

Plasmodium falciparum pfhrp2 and pfhrp3 gene deletions among patients enrolled at 100 health facilities throughout Tanzania: February to July 2021

Plasmodium falciparum with the histidine rich protein 2 gene (pfhrp2) deleted from its genome can escape diagnosis by HRP2-based rapid diagnostic tests (HRP2-RDTs). The World Health Organization (WHO) recommends switching to a non-HRP2 RDT for P. falciparum clinical case diagnosis when pfhrp2 deletion prevalence causes ≥ 5% of RDTs to return false negative results. Tanzania is a country of heterogenous P. falciparum transmission, with some regions approaching elimination and others at varying levels of control. In concordance with the current recommended WHO pfhrp2 deletion surveillance strategy, 100 health facilities encompassing 10 regions of Tanzania enrolled malaria-suspected patients between February and July 2021. Of 7863 persons of all ages enrolled and providing RDT result and blood sample, 3777 (48.0%) were positive by the national RDT testing for Plasmodium lactate dehydrogenase (pLDH) and/or HRP2. A second RDT testing specifically for the P. falciparum LDH (Pf-pLDH) antigen found 95 persons (2.5% of all RDT positives) were positive, though negative by the national RDT for HRP2, and were selected for pfhrp2 and pfhrp3 (pfhrp2/3) genotyping. Multiplex antigen detection by laboratory bead assay found 135/7847 (1.7%) of all blood samples positive for Plasmodium antigens but very low or no HRP2, and these were selected for genotyping as well. Of the samples selected for genotyping based on RDT or laboratory multiplex result, 158 were P. falciparum DNA positive, and 140 had sufficient DNA to be genotyped for pfhrp2/3. Most of these (125/140) were found to be pfhrp2+/pfhrp3+, with smaller numbers deleted for only pfhrp2 (n = 9) or only pfhrp3 (n = 6). No dual pfhrp2/3 deleted parasites were observed. This survey found that parasites with these gene deletions are rare in Tanzania, and estimated that 0.24% (95% confidence interval: 0.08% to 0.39%) of false-negative HRP2-RDTs for symptomatic persons were due to pfhrp2 deletions in this 2021 Tanzania survey. These data provide evidence for HRP2-based diagnostics as currently accurate for P. falciparum diagnosis in Tanzania.

Microsatellites reveal high polymorphism and high potential for use in anti-malarial efficacy studies in areas with different transmission intensities in mainland Tanzania

BACKGROUND

Tanzania is currently implementing therapeutic efficacy studies (TES) in areas of varying malaria transmission intensities as per the World Health Organization (WHO) recommendations. In TES, distinguishing reinfection from recrudescence is critical for the determination of anti-malarial efficacy. Recently, the WHO recommended genotyping polymorphic coding genes, merozoite surface proteins 1 and 2 (msp1 and msp2), and replacing the glutamate-rich protein (glurp) gene with one of the highly polymorphic microsatellites in Plasmodium falciparum to adjust the efficacy of antimalarials in TES. This study assessed the polymorphisms of six neutral microsatellite markers and their potential use in TES, which is routinely performed in Tanzania.

METHODS

Plasmodium falciparum samples were obtained from four TES sentinel sites, Kibaha (Pwani), Mkuzi (Tanga), Mlimba (Morogoro) and Ujiji (Kigoma), between April and September 2016. Parasite genomic DNA was extracted from dried blood spots on filter papers using commercial kits. Genotyping was done using six microsatellites (Poly-α, PfPK2, TA1, C3M69, C2M34 and M2490) by capillary method, and the data were analysed to determine the extent of their polymorphisms and genetic diversity at the four sites.

RESULTS

Overall, 83 (88.3%) of the 94 samples were successfully genotyped (with positive results for ≥ 50.0% of the markers), and > 50.0% of the samples (range = 47.6-59.1%) were polyclonal, with a mean multiplicity of infection (MOI) ranging from 1.68 to 1.88 among the four sites. There was high genetic diversity but limited variability among the four sites based on mean allelic richness (RS = 7.48, range = 7.27-8.03, for an adjusted minimum sample size of 18 per site) and mean expected heterozygosity (He = 0.83, range = 0.80-0.85). Cluster analysis of haplotypes using STRUCTURE, principal component analysis, and pairwise genetic differentiation (FST) did not reveal population structure or clustering of parasites according to geographic origin. Of the six markers, Poly-α was the most polymorphic, followed by C2M34, TA1 and C3M69, while M2490 was the least polymorphic.

CONCLUSION

Microsatellite genotyping revealed high polyclonality and genetic diversity but no significant population structure. Poly-α, C2M34, TA1 and C3M69 were the most polymorphic markers, and Poly-α alone or with any of the other three markers could be adopted for use in TES in Tanzania.

Strong isolation by distance and evidence of population microstructure reflect ongoing Plasmodium falciparum transmission in Zanzibar

The Zanzibar archipelago of Tanzania has become a low-transmission area for Plasmodium falciparum. Despite being considered an area of pre-elimination for years, achieving elimination has been difficult, likely due to a combination of imported infections from mainland Tanzania, and continued local transmission. To shed light on these sources of transmission, we applied highly multiplexed genotyping utilizing molecular inversion probes to characterize the genetic relatedness of 282 P. falciparum isolates collected across Zanzibar and in Bagamoyo District on the coastal mainland from 2016-2018. Overall, parasite populations on the coastal mainland and Zanzibar archipelago remain highly related. However, parasite isolates from Zanzibar exhibit population microstructure due to rapid decay of parasite relatedness over very short distances. This, along with highly related pairs within shehias, suggests ongoing low level local transmission. We also identified highly related parasites across shehias that reflect human mobility on the main island of Unguja and identified a cluster of highly related parasites, suggestive of an outbreak, in the Micheweni district on Pemba island. Parasites in asymptomatic infections demonstrated higher complexity of infection than those in symptomatic infections, but have similar core genomes. Our data support importation as a main source of genetic diversity and contribution to the parasite population on Zanzibar, but they also show local outbreak clusters where targeted interventions are essential to block local transmission. These results highlight the need for preventive measures against imported malaria and enhanced control measures in areas that remain receptive for malaria reemergence due to susceptible hosts and competent vectors.

Using a mobile nanopore sequencing lab for end-to-end genomic surveillance of Plasmodium falciparum: A feasibility study

Genomic epidemiology holds promise for malaria control and elimination efforts, for example by informing on Plasmodium falciparum genetic diversity and prevalence of mutations conferring anti-malarial drug resistance. Limited sequencing infrastructure in many malaria-endemic areas prevents the rapid generation of genomic data. To address these issues, we developed and validated assays for P. falciparum nanopore sequencing in endemic sites using a mobile laboratory, targeting key antimalarial drug resistance markers and microhaplotypes. Using two multiplexed PCR reactions, we amplified six highly polymorphic microhaplotypes and ten drug resistance markers. We developed a bioinformatics workflow that allows genotyping of polyclonal malaria infections, including minority clones. We validated the panels on mock dried blood spot (DBS) and rapid diagnostic test (RDT) samples and archived DBS, demonstrating even, high read coverage across amplicons (range: 580x to 3,212x median coverage), high haplotype calling accuracy, and the ability to explore within-sample diversity of polyclonal infections. We field-tested the feasibility of rapid genotyping in Zanzibar in close collaboration with the local malaria elimination program using DBS and routinely collected RDTs as sample inputs. Our assay identified haplotypes known to confer resistance to known antimalarials in the dhfr, dhps and mdr1 genes, but no evidence of artemisinin partial resistance. Most infections (60%) were polyclonal, with high microhaplotype diversity (median HE = 0.94). In conclusion, our assays generated actionable data within a few days, and we identified current challenges for implementing nanopore sequencing in endemic countries to accelerate malaria control and elimination.

Emergence of Plasmodium falciparum strains with artemisinin partial resistance in East Africa and the Horn of Africa: is there a need to panic?

The emergence and spread of artemisinin partial resistance in East and Horn of Africa is alarming. However, artemisinin-based combination therapy (ACT) generally remains efficacious for the treatment of falciparum malaria. The emergence of partial artemisinin resistance does not currently meet the criteria to initiate change on treatment guidelines nor affect ACT routine procurement and distribution. It is high time for scientists and transitional researchers to be more critical and vigilant on further changes so that national programmes will be able to make informed decisions as well as remain alert and prepared for any change that may be required in the future.

Country wide surveillance reveals prevalent artemisinin partial resistance mutations with evidence for multiple origins and expansion of high level sulfadoxine-pyrimethamine resistance mutations in northwest Tanzania

Background

Emergence of artemisinin partial resistance (ART-R) in Plasmodium falciparum is a growing threat to the efficacy of artemisinin combination therapies (ACT) and the efforts for malaria elimination. The emergence of Plasmodium falciparum Kelch13 (K13) R561H in Rwanda raised concern about the impact in neighboring Tanzania. In addition, regional concern over resistance affecting sulfadoxine-pyrimethamine (SP), which is used for chemoprevention strategies, is high.

Methods

To enhance longitudinal monitoring, the Molecular Surveillance of Malaria in Tanzania (MSMT) project was launched in 2020 with the goal of assessing and mapping antimalarial resistance. Community and clinic samples were assessed for resistance polymorphisms using a molecular inversion probe platform.

Findings

Genotyping of 6,278 samples collected countrywide in 2021 revealed a focus of K13 561H mutants in northwestern Tanzania (Kagera) with prevalence of 7.7% (50/649). A small number of 561H mutants (about 1%) were found as far as 800 km away in Tabora, Manyara, and Njombe. Genomic analysis suggests some of these parasites are highly related to isolates collected in Rwanda in 2015, supporting regional spread of 561H. However, a novel haplotype was also observed, likely indicating a second origin in the region. Other validated resistance polymorphisms (622I and 675V) were also identified. A focus of high sulfadoxine-pyrimethamine drug resistance was also identified in Kagera with a prevalence of dihydrofolate reductase 164L of 15% (80/526).

Interpretation

These findings demonstrate the K13 561H mutation is entrenched in the region and that multiple origins of ART-R, similar as to what was seen in Southeast Asia, have occurred. Mutations associated with high levels of SP resistance are increasing. These results raise concerns about the long-term efficacy of artemisinin and chemoprevention antimalarials in the region.

Funding

This study was funded by the Bill and Melinda Gates Foundation and the National Institutes of Health.

Malaria surveillance reveals parasite relatedness, signatures of selection, and correlates of transmission across Senegal

We here analyze data from the first year of an ongoing nationwide program of genetic surveillance of Plasmodium falciparum parasites in Senegal. The analysis is based on 1097 samples collected at health facilities during passive malaria case detection in 2019; it provides a baseline for analyzing parasite genetic metrics as they vary over time and geographic space. The study's goal was to identify genetic metrics that were informative about transmission intensity and other aspects of transmission dynamics, focusing on measures of genetic relatedness between parasites. We found the best genetic proxy for local malaria incidence to be the proportion of polygenomic infections (those with multiple genetically distinct parasites), although this relationship broke down at low incidence. The proportion of related parasites was less correlated with incidence while local genetic diversity was uninformative. The type of relatedness could discriminate local transmission patterns: two nearby areas had similarly high fractions of relatives, but one was dominated by clones and the other by outcrossed relatives. Throughout Senegal, 58% of related parasites belonged to a single network of relatives, within which parasites were enriched for shared haplotypes at known and suspected drug resistance loci and at one novel locus, reflective of ongoing selection pressure.

Screening for K13-Propeller Mutations Associated with Artemisinin Resistance in Plasmodium falciparum in Yambio County (Western Equatoria State, South Sudan)

Artemisinin-combined treatments are the recommended first-line treatment of Plasmodium falciparum malaria, but they are being threatened by emerging artemisinin resistance. Mutations in pfk13 are the principal molecular marker for artemisinin resistance. This study characterizes the presence of mutations in pfk13 in P. falciparum in Western Equatoria State, South Sudan. We analyzed 468 samples from patients with symptomatic malaria and found 15 mutations (8 nonsynonymous and 7 synonymous). Each mutation appeared only once, and none were validated or candidate markers of artemisinin resistance. However, some mutations were in the same or following position of validated and candidate resistance markers, suggesting instability of the gene that could lead to resistance. The R561L nonsynonymous mutation was found in the same position as the R561H validated mutation. Moreover, the A578S mutation, which is widespread in Africa, was also reported in this study. We found a high diversity of other pfk13 mutations in low frequency. Therefore, routine molecular surveillance of resistance markers is highly recommended to promptly detect the emergence of resistance-related mutations and to limit their spread.

Plasmodium falciparum resistant to artemisinin and diagnostics have emerged in Ethiopia

Diagnosis and treatment of Plasmodium falciparum infections are required for effective malaria control and are pre-requisites for malaria elimination efforts; hence we need to monitor emergence, evolution and spread of drug- and diagnostics-resistant parasites. We deep sequenced key drug-resistance mutations and 1,832 SNPs in the parasite genomes of 609 malaria cases collected during a diagnostic-resistance surveillance study in Ethiopia. We found that 8.0% (95% CI 7.0-9.0) of malaria cases were caused by P. falciparum carrying the candidate artemisinin partial-resistance kelch13 (K13) 622I mutation, which was less common in diagnostic-resistant parasites mediated by histidine-rich proteins 2 and 3 (pfhrp2/3) deletions than in wild-type parasites (P = 0.03). Identity-by-descent analyses showed that K13 622I parasites were significantly more related to each other than to wild type (P < 0.001), consistent with recent expansion and spread of this mutation. Pfhrp2/3-deleted parasites were also highly related, with evidence of clonal transmissions at the district level. Of concern, 8.2% of K13 622I parasites also carried the pfhrp2/3 deletions. Close monitoring of the spread of combined drug- and diagnostic-resistant parasites is needed.

Molecular surveillance of Plasmodium falciparum drug-resistance markers in Vietnam using multiplex amplicon sequencing (2000-2016)

Emergence and spread of Plasmodium falciparum resistance to artemisinin-based combination therapies (ACT) is a major challenge for Greater Mekong Subregion countries in their goal to eliminate malaria by 2030. Tools to efficiently monitor drug resistance beyond resource-demanding therapeutic efficacy studies are necessary. A custom multiplex amplicon sequencing assay based on Illumina technology was designed to target the marker of partial resistance to artemisinin (K13), five candidate modulators of artemisinin resistance, the marker of resistance to chloroquine (crt), and four neutral microsatellite loci. The assay was used to genotype 635 P. falciparum-positive blood samples collected across seven provinces of Vietnam and one of Cambodia between 2000 and 2016. Markers of resistance to artemisinin partner-drugs piperaquine (copy number of plasmepsin-2) and mefloquine (copy number of multidrug-resistance 1) were determined by qPCR. Parasite population structure was further assessed using a 101-SNP barcode. Validated mutations of artemisinin partial resistance in K13 were found in 48.1% of samples, first detection was in 2000, and by 2015 prevalence overcame > 50% in Central Highlands and Binh Phuoc province. K13-C580Y variant became predominant country-wide, quickly replacing an outbreak of K13-I543T in Central Highlands. Mutations in candidate artemisinin resistance modulator genes paralleled the trends of K13 mutants, whereas resistance to piperaquine and mefloquine remained low (≈ 10%) by 2015-2016. Genomic tools applied to malaria surveillance generate comprehensive information on dynamics of drug resistance and population structure and reflect drug efficacy profiles from in vivo studies.

Temporal and spatial analysis of Plasmodium falciparum genomics reveals patterns of parasite connectivity in a low-transmission district in Southern Province, Zambia

BACKGROUND

Understanding temporal and spatial dynamics of malaria transmission will help to inform effective interventions and strategies in regions approaching elimination. Parasite genomics are increasingly used to monitor epidemiologic trends, including assessing residual transmission across seasons and importation of malaria into these regions.

METHODS

In a low and seasonal transmission setting of southern Zambia, a total of 441 Plasmodium falciparum samples collected from 8 neighbouring health centres between 2012 and 2018 were genotyped using molecular inversion probes (MIPs n = 1793) targeting a total of 1832 neutral and geographically informative SNPs distributed across the parasite genome. After filtering for quality and missingness, 302 samples and 1410 SNPs were retained and used for downstream population genomic analyses.

RESULTS

The analyses revealed most (67%, n = 202) infections harboured one clone (monogenomic) with some variation at local level suggesting low, but heterogenous malaria transmission. Relatedness identity-by-descent (IBD) analysis revealed variable distribution of IBD segments across the genome and 6% of pairs were highly-related (IBD ≥ 0.25). Some of the highly-related parasite populations persisted across multiple seasons, suggesting that persistence of malaria in this low-transmission region is fueled by parasites "seeding" across the dry season. For recent years, clusters of clonal parasites were identified that were dissimilar to the general parasite population, suggesting parasite populations were increasingly fragmented at small spatial scales due to intensified control efforts. Clustering analysis using PCA and t-SNE showed a lack of substantial parasite population structure.

CONCLUSION

Leveraging both genomic and epidemiological data provided comprehensive picture of fluctuations in parasite populations in this pre-elimination setting of southern Zambia over 7 years.

Multiplexed ddPCR-amplicon sequencing reveals isolated Plasmodium falciparum populations amenable to local elimination in Zanzibar, Tanzania

Zanzibar has made significant progress toward malaria elimination, but recent stagnation requires novel approaches. We developed a highly multiplexed droplet digital PCR (ddPCR)-based amplicon sequencing method targeting 35 microhaplotypes and drug-resistance loci, and successfully sequenced 290 samples from five districts covering both main islands. Here, we elucidate fine-scale Plasmodium falciparum population structure and infer relatedness and connectivity of infections using an identity-by-descent (IBD) approach. Despite high genetic diversity, we observe pronounced fine-scale spatial and temporal parasite genetic structure. Clusters of near-clonal infections on Pemba indicate persistent local transmission with limited parasite importation, presenting an opportunity for local elimination efforts. Furthermore, we observe an admixed parasite population on Unguja and detect a substantial fraction (2.9%) of significantly related infection pairs between Zanzibar and the mainland, suggesting recent importation. Our study provides a high-resolution view of parasite genetic structure across the Zanzibar archipelago and provides actionable insights for prioritizing malaria elimination efforts.

Performance of SNP barcodes to determine genetic diversity and population structure of Plasmodium falciparum in Africa

Panels of informative biallelic single nucleotide polymorphisms (SNPs) have been proposed to be an economical method to fast-track the population genetic analysis of Plasmodium falciparum in malaria-endemic areas. Whilst used successfully in low-transmission areas where infections are monoclonal and highly related, we present the first study to evaluate the performance of these 24- and 96-SNP molecular barcodes in African countries, characterised by moderate-to-high transmission, where multiclonal infections are prevalent. For SNP barcodes it is generally recommended that the SNPs chosen i) are biallelic, ii) have a minor allele frequency greater than 0.10, and iii) are independently segregating, to minimise bias in the analysis of genetic diversity and population structure. Further, to be standardised and used in many population genetic studies, these barcodes should maintain characteristics i) to iii) across various iv) geographies and v) time points. Using haplotypes generated from the MalariaGEN P. falciparum Community Project version six database, we investigated the ability of these two barcodes to fulfil these criteria in moderate-to-high transmission African populations in 25 sites across 10 countries. Predominantly clinical infections were analysed, with 52.3% found to be multiclonal, generating high proportions of mixed-allele calls (MACs) per isolate thereby impeding haplotype construction. Of the 24- and 96-SNPs, loci were removed if they were not biallelic and had low minor allele frequencies in all study populations, resulting in 20- and 75-SNP barcodes respectively for downstream population genetics analysis. Both SNP barcodes had low expected heterozygosity estimates in these African settings and consequently biased analyses of similarity. Both minor and major allele frequencies were temporally unstable. These SNP barcodes were also shown to identify weak genetic differentiation across large geographic distances based on Mantel Test and DAPC. These results demonstrate that these SNP barcodes are vulnerable to ascertainment bias and as such cannot be used as a standardised approach for malaria surveillance in moderate-to-high transmission areas in Africa, where the greatest genomic diversity of P. falciparum exists at local, regional and country levels.

Malaria surveillance reveals parasite relatedness, signatures of selection, and correlates of transmission across Senegal

Parasite genetic surveillance has the potential to play an important role in malaria control. We describe here an analysis of data from the first year of an ongoing, nationwide program of genetic surveillance of Plasmodium falciparum parasites in Senegal, intended to provide actionable information for malaria control efforts. Looking for a good proxy for local malaria incidence, we found that the best predictor was the proportion of polygenomic infections (those with multiple genetically distinct parasites), although that relationship broke down in very low incidence settings (r = 0.77 overall). The proportion of closely related parasites in a site was more weakly correlated ( r = -0.44) with incidence while the local genetic diversity was uninformative. Study of related parasites indicated their potential for discriminating local transmission patterns: two nearby study areas had similarly high fractions of relatives, but one area was dominated by clones and the other by outcrossed relatives. Throughout the country, 58% of related parasites proved to belong to a single network of relatives, within which parasites were enriched for shared haplotypes at known and suspected drug resistance loci as well as at one novel locus, reflective of ongoing selection pressure.

Design and implementation of multiplexed amplicon sequencing panels to serve genomic epidemiology of infectious disease: A malaria case study

Multiplexed PCR amplicon sequencing (AmpSeq) is an increasingly popular application for cost-effective monitoring of threatened species and managed wildlife populations, and shows strong potential for the genomic epidemiology of infectious disease. AmpSeq data from infectious microbes can inform disease control in multiple ways, such as by measuring drug resistance marker prevalence, distinguishing imported from local cases, and determining the effectiveness of therapeutics. We describe the design and comparative evaluation of two new AmpSeq assays for Plasmodium falciparum malaria parasites: a four-locus panel ("4CAST") composed of highly diverse antigens, and a 129-locus panel ("AMPLseq") composed of drug resistance markers, highly diverse loci for inferring relatedness, and a locus to detect Plasmodium vivax co-infection. We explore the performance of each panel in various public health use cases with in silico simulations as well as empirical experiments. The 4CAST panel appears highly suitable for evaluating the number of distinct parasite strains within samples (complexity of infection), showing strong performance across a wide range of parasitaemia levels without a DNA pre-amplification step. For relatedness inference, the larger AMPLseq panel performs similarly to two existing panels of comparable size, despite differences in the data and approach used for designing each panel. Finally, we describe an R package (paneljudge) that facilitates the design and comparative evaluation of genetic panels for relatedness estimation, and we provide general guidance on the design and implementation of AmpSeq panels for the genomic epidemiology of infectious disease.

Potential Opportunities and Challenges of Deploying Next Generation Sequencing and CRISPR-Cas Systems to Support Diagnostics and Surveillance Towards Malaria Control and Elimination in Africa

Recent developments in molecular biology and genomics have revolutionized biology and medicine mainly in the developed world. The application of next generation sequencing (NGS) and CRISPR-Cas tools is now poised to support endemic countries in the detection, monitoring and control of endemic diseases and future epidemics, as well as with emerging and re-emerging pathogens. Most low and middle income countries (LMICs) with the highest burden of infectious diseases still largely lack the capacity to generate and perform bioinformatic analysis of genomic data. These countries have also not deployed tools based on CRISPR-Cas technologies. For LMICs including Tanzania, it is critical to focus not only on the process of generation and analysis of data generated using such tools, but also on the utilization of the findings for policy and decision making. Here we discuss the promise and challenges of NGS and CRISPR-Cas in the context of malaria as Africa moves towards malaria elimination. These innovative tools are urgently needed to strengthen the current diagnostic and surveillance systems. We discuss ongoing efforts to deploy these tools for malaria detection and molecular surveillance highlighting potential opportunities presented by these innovative technologies as well as challenges in adopting them. Their deployment will also offer an opportunity to broadly build in-country capacity in pathogen genomics and bioinformatics, and to effectively engage with multiple stakeholders as well as policy makers, overcoming current workforce and infrastructure challenges. Overall, these ongoing initiatives will build the malaria molecular surveillance capacity of African researchers and their institutions, and allow them to generate genomics data and perform bioinformatics analysis in-country in order to provide critical information that will be used for real-time policy and decision-making to support malaria elimination on the continent.

Low Complexity of Infection Is Associated With Molecular Persistence of Plasmodium falciparum in Kenya and Tanzania

Background

Plasmodium falciparum resistance to artemisinin-based combination therapies (ACTs) is a threat to malaria elimination. ACT-resistance in Asia raises concerns for emergence of resistance in Africa. While most data show high efficacy of ACT regimens in Africa, there have been reports describing declining efficacy, as measured by both clinical failure and prolonged parasite clearance times.

Methods

Three hundred children aged 2-10 years with uncomplicated P. falciparum infection were enrolled in Kenya and Tanzania after receiving treatment with artemether-lumefantrine. Blood samples were taken at 0, 24, 48, and 72 h, and weekly thereafter until 28 days post-treatment. Parasite and host genetics were assessed, as well as clinical, behavioral, and environmental characteristics, and host anti-malarial serologic response.

Results

While there was a broad range of clearance rates at both sites, 85% and 96% of Kenyan and Tanzanian samples, respectively, were qPCR-positive but microscopy-negative at 72 h post-treatment. A greater complexity of infection (COI) was negatively associated with qPCR-detectable parasitemia at 72 h (OR: 0.70, 95% CI: 0.53-0.94), and a greater baseline parasitemia was marginally associated with qPCR-detectable parasitemia (1,000 parasites/uL change, OR: 1.02, 95% CI: 1.01-1.03). Demographic, serological, and host genotyping characteristics showed no association with qPCR-detectable parasitemia at 72 h. Parasite haplotype-specific clearance slopes were grouped around the mean with no association detected between specific haplotypes and slower clearance rates.

Conclusions

Identifying risk factors for slow clearing P. falciparum infections, such as COI, are essential for ongoing surveillance of ACT treatment failure in Kenya, Tanzania, and more broadly in sub-Saharan Africa.

Malaria protection due to sickle haemoglobin depends on parasite genotype

Host genetic factors can confer resistance against malaria¹, raising the question of whether this has led to evolutionary adaptation of parasite populations. Here we searched for association between candidate host and parasite genetic variants in 3,346 Gambian and Kenyan children with severe malaria caused by Plasmodium falciparum. We identified a strong association between sickle haemoglobin (HbS) in the host and three regions of the parasite genome, which is not explained by population structure or other covariates, and which is replicated in additional samples. The HbS-associated alleles include nonsynonymous variants in the gene for the acyl-CoA synthetase family member^2–4PfACS8 on chromosome 2, in a second region of chromosome 2, and in a region containing structural variation on chromosome 11. The alleles are in strong linkage disequilibrium and have frequencies that covary with the frequency of HbS across populations, in particular being much more common in Africa than other parts of the world. The estimated protective effect of HbS against severe malaria, as determined by comparison of cases with population controls, varies greatly according to the parasite genotype at these three loci. These findings open up a new avenue of enquiry into the biological and epidemiological significance of the HbS-associated polymorphisms in the parasite genome and the evolutionary forces that have led to their high frequency and strong linkage disequilibrium in African P. falciparum populations.

A strong association has been found between three regions of the Plasmodium falciparum genome and sickle haemoglobin in children with severe malaria, suggesting parasites have adapted to overcome natural host immunity.

Immune selection suppresses the emergence of drug resistance in malaria parasites but facilitates its spread

Although drug resistance in Plasmodium falciparum typically evolves in regions of low transmission, resistance spreads readily following introduction to regions with a heavier disease burden. This suggests that the origin and the spread of resistance are governed by different processes, and that high transmission intensity specifically impedes the origin. Factors associated with high transmission, such as highly immune hosts and competition within genetically diverse infections, are associated with suppression of resistant lineages within hosts. However, interactions between these factors have rarely been investigated and the specific relationship between adaptive immunity and selection for resistance has not been explored. Here, we developed a multiscale, agent-based model of Plasmodium parasites, hosts, and vectors to examine how host and parasite dynamics shape the evolution of resistance in populations with different transmission intensities. We found that selection for antigenic novelty (“immune selection”) suppressed the evolution of resistance in high transmission settings. We show that high levels of population immunity increased the strength of immune selection relative to selection for resistance. As a result, immune selection delayed the evolution of resistance in high transmission populations by allowing novel, sensitive lineages to remain in circulation at the expense of the spread of a resistant lineage.

In contrast, in low transmission settings, we observed that resistant strains were able to sweep to high population prevalence without interference. Additionally, we found that the relationship between immune selection and resistance changed when resistance was widespread. Once resistance was common enough to be found on many antigenic backgrounds, immune selection stably maintained resistant parasites in the population by allowing them to proliferate, even in untreated hosts, when resistance was linked to a novel epitope. Our results suggest that immune selection plays a role in the global pattern of resistance evolution.

Drug resistance in the malaria parasite, Plasmodium falciparum, presents an ongoing public health challenge, but aspects of its evolution are poorly understood. Although antimalarial resistance is common worldwide, it can typically be traced to just a handful of evolutionary origins. Counterintuitively, although Sub Saharan Africa bears 90% of the global malaria burden, resistance typically originates in regions where transmission intensity is low. In high transmission regions, infections are genetically diverse, and hosts have significant standing adaptive immunity, both of which are known to suppress the frequency of resistance within infections. However, interactions between immune-driven selection, transmission intensity, and resistance have not been investigated. Using a multiscale, agent-based model, we found that high transmission intensity slowed the evolution of resistance via its effect on host population immunity. High host immunity strengthened selection for antigenic novelty, interfering with selection for resistance and allowing sensitive lineages to suppress resistant lineages in untreated hosts. However, once resistance was common in the circulating parasite population, immune selection maintained it in the population at a high prevalence. Our findings provide a novel explanation for observations about the origin of resistance and suggest that adaptive immunity is a critical component of selection.

A review of the frequencies of Plasmodium falciparum Kelch 13 artemisinin resistance mutations in Africa

Artemisinin resistance (AR) emerged in South East Asia 13 years ago and the identification of the resistance conferring molecular marker, Plasmodium falciparum Kelch 13 (Pfk13), 7 years ago has provided an invaluable tool for monitoring AR in malaria endemic countries. Molecular Pfk13 surveillance revealed the resistance foci in the Greater Mekong Subregion, an independent emergence in Guyana, South America, and a low frequency of mutations in Africa. The recent identification of the R561H Pfk13 AR associated mutation in Tanzania, Uganda and in Rwanda, where it has been associated with delayed parasite clearance, should be a concern for the continent. In this review, we provide a summary of Pfk13 resistance associated propeller domain mutation frequencies across Africa from 2012 to 2020, to examine how many other countries have identified these mutations. Only four African countries reported a recent identification of the M476I, P553L, R561H, P574L, C580Y and A675V Pfk13 mutations at low frequencies and with no reports of clinical treatment failure, except for Rwanda. These mutations present a threat to malaria control across the continent, since the greatest burden of malaria remains in Africa. A rise in the frequency of these mutations and their spread would reverse the gains made in the reduction of malaria over the last 20 years, given the lack of new antimalarial treatments in the event artemisinin-based combination therapies fail. The review highlights the frequency of Pfk13 propeller domain mutations across Africa, providing an up-to-date perspective of Pfk13 mutations, and appeals for an urgent and concerted effort to monitoring antimalarial resistance markers in Africa and the efficacy of antimalarials by re-establishing sentinel surveillance systems.

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.

Describing the current status of Plasmodium falciparum population structure and drug resistance within mainland Tanzania using molecular inversion probes

High-throughput Plasmodium genomic data is increasingly useful in assessing prevalence of clinically important mutations and malaria transmission patterns. Understanding parasite diversity is important for identification of specific human or parasite populations that can be targeted by control programmes, and to monitor the spread of mutations associated with drug resistance. An up-to-date understanding of regional parasite population dynamics is also critical to monitor the impact of control efforts. However, this data is largely absent from high-burden nations in Africa, and to date, no such analysis has been conducted for malaria parasites in Tanzania countrywide. To this end, over 1,000 P. falciparum clinical isolates were collected in 2017 from 13 sites in seven administrative regions across Tanzania, and parasites were genotyped at 1,800 variable positions genome-wide using molecular inversion probes. Population structure was detectable among Tanzanian P. falciparum parasites, approximately separating parasites from the northern and southern districts and identifying genetically admixed populations in the north. Isolates from nearby districts were more likely to be genetically related compared to parasites sampled from more distant districts. Known drug resistance mutations were seen at increased frequency in northern districts (including two infections carrying pfk13-R561H), and additional variants with undetermined significance for antimalarial resistance also varied by geography. Malaria Indicator Survey (2017) data corresponded with genetic findings, including average region-level complexity-of-infection and malaria prevalence estimates. The parasite populations identified here provide important information on extant spatial patterns of genetic diversity of Tanzanian parasites, to which future surveys of genetic relatedness can be compared.