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 ======================================================================================================================================================================================================================== * Jonathan J. Juliano * David J. Giesbrecht * Alfred Simkin * Abebe A. Fola * Beatus M. Lyimo * Dativa Pereus * Catherine Bakari * Rashid A. Madebe * Misago D. Seth * Celine I. Mandara * Zachary R. Popkin-Hall * Ramadhan Moshi * Ruth B. Mbwambo * Karamoko Niaré * Bronwyn MacInnis * Filbert Francis * Daniel Mbwambo * Issa Garimo * Frank Chacky * Sijenunu Aaron * Abdallah Lusasi * Fabrizio Molteni * Ritha J. A. Njau * Samwel Lazaro * Ally Mohamed * Jeffrey A. Bailey * Deus S. Ishengoma ## Abstract **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. **Evidence before this study** We did a literature search via PubMed for research articles published from January 2014 to October 2023 using the search term “Africa” and “Artemisinin resistance” linked to “R561H” or “A675V” or “R622I”, returning 32 studies. The published literature shows the emergence and establishment of these three validated *Plasmodium falciparum* kelch13 (K13) mutations associated with artemisinin partial resistance (ART-R) in Africa. Large molecular studies of 675V in Uganda and 622I in Ethiopia have defined the regional spread of these mutations. However, limited data is available from recent studies about the spread and origins of the 561H mutation in the Great Lakes region of East Africa. In particular, detailed studies of the regions of Tanzania that border Rwanda have not been carried out since the mutation was detected in Rwanda. These data are needed for malaria control programs to define and implement strategies for controlling the spread of ART-R in Africa, a potential global public health disaster and the potential obstacle to the ongoing elimination strategies. **Added value of this study** This analysis reports the first large-scale analysis of antimalarial resistance in Tanzania, with a focus on the regions bordering Rwanda since the 561H mutation reached high frequency in the area. Using 6,278 *P. falciparum* positive samples sequenced using molecular inversion probes (MIPs), we show that the mutation has become frequent in the districts of Kagera bordering Rwanda. Importantly, we provide evidence for the separate emergence of a different extended haplotype around 561H in Tanzania. This is the first evidence that multiple independent emergences of the 561H ART-R have occurred in Africa, as was seen within the last two decades in Southeast Asia. **Implications of all the available evidence** These findings highlight that, similar to 622I and 675V in other parts of Africa, we can expect the 561H mutation to continue to spread in the region. In addition, it highlights that we need to be watchful for new origins of mutations beyond the spread of existing resistant parasite lineages. ART-R appears to now be well established in multiple areas in Eastern Africa. Intensive control in these regions to prevent spread and monitoring for partner drug resistance emergence in affected areas will be critical for preventing further reversal of malaria control efforts in the region and support progress to the elimination targets by 2023. ## INTRODUCTION Resistance to drugs used for treatment and prevention of malaria poses one of the greatest threats to global control and is of grave concern in Africa where the vast majority of cases and deaths occur.1 Historically, the emergence and spread of chloroquine and sulfadoxine-pyrimethamine (SP) resistance resulted in the collapse of effective treatment of malaria with significant increases in morbidity and mortality.2 Antimalarial combination therapies, consisting of an artemisinin derivative and a partner drug (artemisinin-based combination therapy - ACTs), are currently the predominant therapeutic options for treating uncomplicated falciparum malaria. Given gains in malaria control have already plateaued and are reversing in some countries, the emergence of artemisinin partial resistance (ART-R) in Africa could be a global public health disaster if partner drug resistance emerges in concert resulting in frank ACT failure.3 Over a decade has passed since the emergence of ART-R in SouthEast Asian *P. falciparum* populations, leading to decreased drug efficacy.4 Clinical ART-R was first demonstrated in the late 2000s in studies conducted in Western Cambodia. The emergence of ART-R in Western Cambodia set the stage for the eventual failure of ACTs, as resistance to the partner drugs also emerged, indicated by increasing treatment failures, parasite clearance times and partner drug IC50s over a very short timeframe.3 Early in the emergence of ACT resistance, many areas of SEA experienced more than 50% failure in patients treated with ACTs due to the combined effect of artemisinin and partner drug resistance.3 Mutations in *Plasmodium falciparum* Kelch13 (K13) (PF3D7_1343700) are the key mediator of ART-R. These mutations were originally identified through drug pressure experiments and validated in the field and by genetic engineering studies.5 They are thought to alter ubiquitination patterns and help parasites to resist accumulation of polyubiquitinated proteins.5 A number of key mutations in K13 propeller domain are now validated markers of ART-R. Unfortunately, validated K13 mutations have now been found extensively in Eastern Africa, particularly in the Horn of Africa (R622I), Uganda (C469Y, A675V), Rwanda (R561H), and Tanzania (R561H).6–10 The regional emergence and spread of these mutations appears to be developing as the norm rather than the exception. In the Horn of Africa, the validated 622I mutation was first reported in 2014.11 Since that time, these mutations have spread across Eritrea and Ethiopia.7,8 In Uganda, longitudinal molecular surveillance at 16 sites has painted a clear picture of spread within Uganda of the 469Y and 675V variants.6 It is important to note that all three of these mutations have been associated with prolonged parasite clearance, day 3 positivity by microscopy, increasing prevalence over time, and are in regions that now meet WHO criteria for ART-R.7,8 The emergence and spread of the 561H mutation in the Great Lakes Region of East Africa is not yet as clear when compared to other mutations. Originally described in Rwanda in samples from 2014 and 2015, this mutation appears to have emerged within the country over the past decade.10,12 In 2015, 7.4% of samples collected in Masaka harbored the 561H mutation.10 By 2018, 561H prevalence had increased to 19.6% in Masaka and 22% in Rukara during a therapeutic efficacy study.13 In this study, 50% of isolates with delayed clearance (day 3 positive parasitemia) carried 561H. On this basis, Masaka met the WHO criteria for endemic ART-R as defined by >5% of patients carrying K13 resistance-confirmed mutation for which they had persistent parasitemia by microscopy on day 3.14 Genome sequencing of isolates from 2015 also confirmed a single haplotype of 561H in Rwanda that was not of Asian origin, suggesting *de novo* mutation within Africa.10 Importantly, K13 mutations should never be studied in isolation as it requires partner drug resistance to lead to clinical failures of ACTs. Further, mutations to drugs no longer used for therapy, but that remain in use for chemoprevention (e.g. SP), are also important to characterize for malaria control programs. Molecular surveillance of antimalarial resistance should therefore rely on platforms that can broadly detect different resistance mutations. Highly multiplex amplicon deep sequencing is one approach that has shown promise. Another approach is molecular inversion probes (MIPs), which have now been used extensively to characterize drug resistance and population structure in parasites in Africa.8,9,15 The ability to create and combine different highly multiplexed panels for antimalarial resistance mutations, copy number variation, gene deletions and other genome wide polymorphisms to study complexity of infection, parasite relatedness or population structure makes the MIP platform a highly flexible and cost-effective means of conducting malaria molecular surveillance (MMS). The Molecular Surveillance of Malaria in Tanzania (MSMT) project was developed to provide nationwide longitudinal surveillance of parasite populations to understand key aspects of parasite biology that may impact malaria control and interventions.16 The emergence of 561H is a prime concern for Tanzania given its proximity to Rwanda and previous studies documenting isolated cases across the country, with two cases in the Chato district and one on the Eastern coast.9,17 Here we describe the initial assessment of 6,278 successfully genotyped samples collected across 7,782 malaria positive individuals in the first year of the project using high-throughput MIP analysis. The goal was to assess the status of antimalarial resistance in Tanzania, with a focus on the border with Rwanda, to understand the distribution of the 561H mutation, partner drug resistance and resistance to chemoprevention drugs. This was further supported by whole genome sequencing of specific isolates to understand flanking haplotypes and gain insight into the origins of 561H mutations in Tanzania. ## MATERIALS AND METHODS ### Study Design and Participants We analyzed dried blood spot (DBS) samples from patients with a positive malaria rapid diagnostic test drawn from cross-sectional surveys that involved 100 health facilities in 10 regions (n=7,148) and asymptomatic individuals in community studies (n=634) in three additional regions (**Figure 1**). Details of the study sites’ selection and overall sampling have been provided elsewhere.16 The three community surveys were conducted in regions of Tanzania which were involved in previous studies undertaken by the National Institute for Medical Research (NIMR).9,18 Four regions (Kagera, Mara, Tabora, and Kigoma) were deemed high priority areas for antimalarial resistance surveillance due to their proximity to Rwanda and to their location in the Lake Zone of the country with higher transmission (**Figure 1**). Informed consent was obtained for each patient and de-identified DBS samples were processed at NIMR in Tanzania, Brown University, and University of North Carolina (USA) according to IRB requirements of the Tanzanian Medical Research Coordinating Committee (MRCC) of NIMR. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/11/30/2023.11.07.23298207/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2023/11/30/2023.11.07.23298207/F1) Figure 1. Tanzania malaria prevalence and study design. Panel A depicts a map of Tanzania showing malaria prevalence in children aged 0-59 months in 2017 as a color gradient (National Malaria Control Program data). Health facilities where sampling occurred for this study are shown as green dots, while community cross sectional collection sites are shown as yellow squares. Panel B shows the study design flowchart showing the number of samples collected at health facilities and cross sectional surveys and the success of genotyping. Priority for genotyping and analysis was given to regions in the northwest. ### Molecular Inversion Probe Analysis DNA was extracted from the DBS using a Chelex-Tween protocol and MIP captures and sequencing were conducted as previously described.15 This study used a panel specific for drug resistance polymorphism detection and a panel to look at genomic diversity.15 Approximately two to three thousand samples were run together on each NextSeq 500 run, and sample libraries lacking sufficient read depth were rebalanced and resequenced. For samples from the high priority regions (**Figure 1**), an additional MIP capture and high-depth sequencing was conducted. Resulting data was analyzed using MIPtools software with freebayes variant calling ([https://github.com/bailey-lab/MIPTools](https://github.com/bailey-lab/MIPTools)).15 Controls for each MIP capture and sequencing included DNA from 3D7 and 7G8 as well as no template and no probe controls. Variant calling was conducted as previously described.15 We kept samples that had at least one haplotype that mapped in the expected locations of the genome for any of our drug resistance MIPs. Antimalarial resistance prevalence was calculated for all variants with a UMI count of 3 or greater and if heterozygous with the alternate allele having 1 UMI or greater using a Python script and maps were created using the sf package. Analysis of haplotypes involves only samples where complete genotypes across the involved loci are available. Inheritance by descent (IBD) analysis of parasite relatedness among ART-R parasitemias was done as previously described.8 ### Whole Genome Sequencing Whole genome sequencing of selective whole genome amplification (sWGA) products was attempted for 23 pure 561H, 5 mixed R561H, and 45 wildtype parasite infections based on MIP genotyping. sWGA was performed in triplicate for each sample using a previously published protocol and pooled.19 The pooled sWGA product was sheared using a LE220R-plus Covaris Sonicator and libraries prepared using dual indexing with the Kappa Hyper Prep Kit (Roche, Indianapolis, IN). Pooled libraries were sequenced on a NovaSeq6000 using 2 x 150 bp chemistry at the University of North Carolina (UNC) High Throughput Sequencing Facility. We also downloaded publicly available WGS data (n=25) from *P. falciparum* isolates collected in 2014/15 in Rwanda.10 Whole genome sequencing data was analyzed using GATK4 following previously published methods ([https://github.com/Karaniare/Optimized\_GATK4_pipeline](https://github.com/Karaniare/Optimized_GATK4_pipeline)).20 Briefly, reads were mapped to the 3D7 reference genome using bwa mem, variants were called using GATK4 and SNPs and indels were filtered using variant quality score recalibration. SNPs and indels passing filters were visualized in R 4.2.1 using the gt package. In order to detect patterns of selection signals between Southeast Asian (SEA), Rwanda and Tanzania haplotypes, we did extended-haplotype homozygosity (EHH) statistics focusing 561H drug resistance SNP using filtered biallelic SNPs and with low missingness data from a VCF file. All associated EHH calculations were carried out using the R-package rehh (version 2.0.4). Rwanda and SEA data was downloaded from Pf7K. ## RESULTS We successfully genotyped 6,278 of the 7,782 (80.7%) samples attempted (**Figure 1B**). Sequencing across the K13 gene revealed three WHO validated ART-R mutations, 561H, 622I, and 675V. The 561H mutations were predominantly found in Kagera, the northwest most region bordering Rwanda and Uganda. (**Figure 2A**, **Table 1**). The overall prevalence in Kagera was 7.7% (50/649) and most of these mutant parasites were found in the west near the Rwandan border in the districts of Karagwe at 22.8% (31/136), Kyerwa at 14.4%% (17/118) and pfdhfr and pfdhps mutations were common in the area studied (Table 2). Notably, we identified Ngara, at 1.4% (2/144) (**Figure 2B**). Beyond Kagera, parasites with the 561H mutation were also detected at lower prevalence in three other regions, including Tabora at 0.5% (2/438), Manyara at 0.5% (1/179), and Njombe at 0.4% (1/279), as far as 800 km from the Rwanda border. The K13 622I mutation was found in a single isolate from Njombe in Southwestern Tanzania (**Table 1**). The 675V mutation was found in one isolate from Kagera and one isolate from Tabora. Other K13 propeller domain mutations not known to be associated with drug resistance were detected sporadically (**Table S1**). Outside the K13 propeller domain, we found several polymorphisms including 189T in 18.3% (397/2,168) of samples country-wide (**Table S1**). Other mutations have putatively been associated with partial artemisinin resistance outside of K13, including *Plasmodium falciparum* ATPase 6 (ATPase6) and *P. falciparum* AP-2mu (AP2mu) (**Table S1**). The ATPase6 623E mutation was found in four Kagera isolates and seven Mara isolates (0.6% and 2.3%, respectively), and none of these isolates contained parasites with mutations at K13 561H. *P. falciparum* coronin was not genotyped. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/11/30/2023.11.07.23298207/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2023/11/30/2023.11.07.23298207/F2) Figure 2. Antimalarial resistance polymorphisms in Tanzania. Panels A and B show K13 561H mutations in Tanzania by region (A) and by district in Kagera (B). In Panel A, white areas represent no sampling, gray areas represent sampled areas where no K13 561H was found, and colored areas are shaded by the frequency of samples with the 56H mutation. The red box surrounds the area plotted in Panel B, which is a map of Kagera showing K13 561H prevalence in each district sampled. DHFR 164L (Panel C) and DHPS 581G (Panel D) prevalence in sampled regions are shown in color gradients. Regions with no sampling are shown in white and regions where a polymorphism was not detected is shown in grey. View this table: [Table 1.](http://medrxiv.org/content/early/2023/11/30/2023.11.07.23298207/T1) Table 1. Prevalence of antimalarial resistance mutations by region. Prevalence reflects the alternate allele. Extended haplotypes around antimalarial genes have been used to study the origin and spread of mutations.10 With positive directional selection of resistance by drug pressure, large genomic sections are spread in the population, a selective sweep, which are eventually broken down through recombination. We achieved sufficient depth and quality on 29 of 63 Tanzanian isolates from 2021 to assess the variation surrounding K13 and used publicly available WGS data from the isolates collected in 2014/15 in Rwanda.10 We identify a shared haplotype between the older Rwandan isolates and the contemporary Tanzanian isolates (n=4 from Tanzania; haplotype RW/TZ1), suggesting cross-border spread resulting from a single origin event. However, a second extended haplotype (n=5; TZ2) is also seen within the Tanzanian parasites (**Figure 3**). RW/TZ1 and TZ2 differ at the closest flanking single nucleotide polymorphisms (within 1kb of the 561H mutation). The TZ2 haplotype also does not appear to be of SEA origin (**Figure S1**). Extended haplotype homozygosity analysis (EHH) shows extended haplotypes for RW/TZ1and TZ2 relative to wild type parasites, consistent with positive selection for both 561H haplotypes (**Figure 4**). IBD analysis showed that 561H isolates from Kagera were closely related in multiple small clusters (**Figure S2**). Of the more distant isolates, the two from Tabora were related to Kagera isolates (IBD≥0.25), with the other two being more distant (**Figure S3**). ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/11/30/2023.11.07.23298207/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2023/11/30/2023.11.07.23298207/F3) Figure 3. Extended flanking haplotype plot around K13 561H mutations. Pure mutant 561H are shown in maroon. Rwandan isolates are in the blue box, while Tanzanian isolates are shown in the red boxes. Grey regions represent loci where a call could not be made. ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/11/30/2023.11.07.23298207/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2023/11/30/2023.11.07.23298207/F4) Figure 4. Extended Haplotype Heterozygosity (EHH) analysis of 561H haplotypes. Panel A shows the extended haplotype for the RW/TZ1 561H haplotype (red-derived) compared to wildtype parasites (blue-ancestral). Panel B shows the extended haplotype of the TZ2 561H haplotype. To investigate associations with mutations previously found to occur in drug resistant parasites in Southeast Asia in combination with K13 mutations, we used the MIP panel genotypes and quantified mutations occurring in the same isolates. None of *P. falciparum* chloroquine resistance transporter (CRT) 326S, CRT 356T, *P.*multidrug resistance protein 2 (MDR2) 484I, *Plasmodium falciparum* Protein phosphatase (PPH) 1157L, or putative phosphoinositide-binding protein (PIB7) 1484F were associated with K13 561H mutant isolates. However, ferredoxin (FD) 193Y polymorphisms were found in 13 K13 561H mutant isolates in Kagera (**Table S1**). Known partner drug resistance mutations were also assessed (**Table 1**). In Kagera, 39 samples (39/586) had mutations at CRT 76T, and 31 of these isolates were from Ngara bordering Rwanda and Burundi (ASAQ is the most widely used ACT in Burundi). The multidrug resistance transporter 1 (MDR1) N86 allele associated with lumefantrine resistance was found in 99.9% of samples (3828/3830). The MDR1 NFD haplotype (N86, 184F, D1246) was found in 52.6% of samples (1420/2698). Markers for SP resistance were also genotyped, revealing high frequency of many mutations and the emergence of dihydrofolate reductase (DHFR) 164L in Kagera (**Table 1**, **Figure 2C**). The folate synthesis gene triple mutations DHFR 51I, 59R and 108N were found to be near fixation with the IRN haplotype at 92.5% (2,893/3,128). Notably, the DHFR 164L mutation was found at 15.2% (80/526) in Kagera region, accounting for 80% of the IRNL haplotypes found in the country (100/2,524 nationally) (**Figure 2C**). In the dihydropteroate synthase (DHPS) gene, two markers for elevated SP resistance were found to be above WHO guideline thresholds in some regions. DHPS 540E was found at 93.06% (3,485/3,745) in Tanzania with little geographic variation (**Table 1**), while DHPS 581G prevalence ranged from 0.4% (1/239) in Mtwara to 38.5% (25/65) in Tanga (**Figure 2D**). ## DISCUSSION Drug resistance in malaria parasites has emerged multiple times to every frontline antimalarial, each time with severe consequences for affected populations.21 The emerging threat of ART-R in Africa has the potential to be a global public health disaster. Thus, understanding the emergence and spread of these mutations is critical in making plans to contain this threat. While it is clear that K13 622I and 675V are spreading in the Horn of Africa and Uganda, respectively, the pattern of origin and spread of the 561H mutation have not been fully defined. Here we demonstrate that the mutation has risen to high frequencies in Kagera (7.7%), a region bordering Rwanda, which had little evidence of resistance in a 2017 study that used MIPs to genotype samples from Kagera in a national schoolkids survey (data not published). This trend would be consistent with outward spread from the Rwanda-Tanzania border, and also supports the idea that there has been and will continue to be regional spread of the mutation, especially given the shared extended haplotype around some of the Tanzanian isolates and the older reported Rwandan isolates. More concerning is the evidence of a potential unique origin of the 561H mutation in the region, suggesting that ART-R may follow patterns similar to those seen in Southeast Asia with multiple independent origins of the same mutation.22 Haplotype analysis suggests multiple origins, but does not prove the mechanism behind them.6 If multiple origins occur, containment efforts are more difficult, due to the need to closely monitor for new haplotypes and the spread of specific haplotypes, which may not become extensive without partner drug resistance. However, the presence of multiple haplotypes may be due to gene conversion or recombination events adjacent to the mutation. Several mutations in the K13 propeller domain that are now associated with delayed parasite clearance are present in East Africa.6–8 The most concerning of these, 469Y, 561H, 622I and 675V, have shown clinical and *in vitro* validation of ART-R. Determining the exact origin of these mutations with certainty is difficult due to the patchy nature of the surveillance data, but one (561H) of the three appears to have originated on the Rwanda-Tanzania border. 561H was first detected in Rwanda in 2015 as part of therapeutic efficacy studies and later in DRC and Tanzania.9,10 To date, genotyping efforts to describe the status of K13 mutations in Tanzania have been limited and do not capture the risk posed by ART-R parasites across the country. The high rate of human movement across the border with Rwanda, where people have close historical ties and One Stop Border Posts (OSPB) that routinely permit thousands of individuals and hundreds of vehicles to cross daily, makes this a primary concern for malaria control in Tanzania.23 While there have been sporadic reports of validated ART-R polymorphisms in the past, systematic longitudinal surveillance of these mutations is necessary to identify areas at risk, as has been done in Uganda.6,9,17 The speed at which data is generated and reported to control programs is critical. The MSMT project will provide these nationwide data for Tanzania given the identified threat shown here of 561H in Kagera, by allowing yearly sampling and in-country sequencing and analysis of data as the project develops. Partner drug resistance surveillance showed important patterns. An intriguing pattern of CRT 76T mutations was found to be clustered near the Burundi border. A possible explanation is that the use of ASAQ in Burundi is contributing to continued selection for chloroquine resistance and that these resistant parasites are being imported into Kagera. Given artemether-lumefantrine remains the first line antimalarial in Tanzania, it is not surprising that MDR1 N86 was near fixation and the NFD haplotype was very common, occurring in over 50% of samples. Markers of elevated SP resistance are increasing in Tanzania, which has potential impacts on the use of SP as a chemopreventive antimalarial. The geographic distribution of the DHPS 581G mutations in Tanzania has historically been limited to Tanga and the Lake districts.24,25 We confirm that this remains the case with little evidence of expansion into other regions of Tanzania. The reasons for the lack of spread of this mutation in Tanzania remains unclear and warrants further evaluation, especially given the evidence that spread is occurring in the neighboring Democratic Republic of the Congo.26 Previous work has suggested that multiple origins of parasites bearing 581G have occurred in East Africa, suggesting local origins and limited spread.27 However, that work was based on five microsatellites and may not take into account the complex stepwise evolutionary history of mutations in DHPS. Importantly, a new focus of the quadruple mutation in DHFR (51I, 59R, 108N, 164L) was identified in Kagera. In the laboratory, this form of the enzyme binds pyrimethamine 600 times less tightly than the wild type and about seven times less than the triple mutation (51I, 59R, 108N).28 This results in parasites being resistant to therapy *in vitro* at levels higher than what can be reached *in vivo*. The DHFR/DHPS sextuple mutation (DHFR 51I, 59R, 108N / DHPS 437G, 540E, 581G), has been associated with the compromise of intermittent preventive treatment in pregnancy (IPTp) and was found in 11.2% (323/2,894) of samples nationally.29 Given the regional high prevalence of DHPS 581G and the emergence of DHFR 164L in Northwest Tanzania, additional studies are warranted to evaluate the benefits of SP based chemoprevention in this region of Tanzania. A major advantage of this study is the large geographic range and large number of samples available for analysis. However, this study design also engenders some limitations. The large number of samples resulted in overall low levels of coverage for the sequencing for most regions (other than the four high priority regions that received extra depth). The samples collected at the three cross-sectional sites performed poorly with our approach, likely due to low parasitemia. The overall coverage in many regions may directly affect the ability to find minority drug resistant variants in the data and we may have filtered out minor variants below 1%. The detection of K13 561H at prevalence above 22% in Karagwe district of Kagera, Tanzania, must raise alarms thatART-R is emerging and that ACTs efficacy is under threat. The independent origin of a new K13 561H extended haplotype raises concern that multiple origins of ART-R mutations will occur in Africa, complicating control. Evolution of partner drug resistance must be monitored carefully and therapeutic efficacy studies are urgently needed to understand the susceptibility of currently circulating K13 mutations in the region. ## Authors contribution DSI, JJJ and JB formulated the original idea. DSI, CIM, RAM, MDS and CB conducted the field surveys. DG, AF, AS, BL, ZPH, KN and AL performed data analysis under the guidance of DSI, JJJ and JB. JJJ, AS, AF, DG, BL, ZPH, CB, DP, MDS, JB, and DSI wrote and edited the manuscript. All authors contributed to the article and approved the submitted version. ## Supporting information Supplemental Material [[supplements/298207_file02.docx]](pending:yes) Table S1 [[supplements/298207_file03.xlsx]](pending:yes) ## Funding This work was supported, in whole, by the Bill & Melinda Gates Foundation [grant number 002202]. Under the grant conditions of the Foundation, a Creative Commons Attribution 4.0 Generic License has already been assigned to the Author Accepted Manuscript version that might arise from this submission. This study was also funded by the National Institute for Allergy and Infectious Diseases (R01AI156267 to JAB, DSI and JJJ; K24AI134990 to JJJ). ## Data Availability All sequencing data has been submitted to SRA (# pending). Other metadata is available from the corresponding author upon reasonable request. Code for analysis is available at: [https://github.com/bailey-lab/MSMT\_2021\_DR\_analyses](https://github.com/bailey-lab/MSMT_2021_DR_analyses). ## Acknowledgements The authors wish to thank participants and parents/guardians of all children who took part in the surveillance. We acknowledge the contribution of the following project staff and other colleagues who participated in data collection and/or laboratory processing of samples; Raymond Kitengeso, Ezekiel Malecela, Muhidin Kassim, Athanas Mhina, August Nyaki, Juma Tupa, Anangisye Malabeja, Emmanuel Kessy, George Gesase, Tumaini Kamna, Grace Kanyankole, Oswald Osca, Richard Makono, Ildephonce Mathias, Godbless Msaki, Rashid Mtumba, Gasper Lugela, Gineson Nkya, Daniel Chale, Richard Malisa, Sawaya Msangi, Ally Idrisa, Francis Chambo, Kusa Mchaina, Raymond Kitengeso, Neema Barua, Christian Msokame, Rogers Msangi, Salome Simba, Hatibu Athumani, Mwanaidi Mtui, Rehema Mtibusa, Jumaa Akida, Ambele Yatinga, Tilaus Gustav, Oksana Kharabora and Claudia Gaither. The finance, administrative and logistic support team at NIMR: Christopher Masaka, Millen Meena, Beatrice Mwampeta, Gracia Sanga, Neema Manumbu, Halfan Mwanga, Arison Ekoni, Twalipo Mponzi, Pendael Nasary, Denis Byakuzana, Alfred Sezary, Emmanuel Mnzava, John Samwel, Daud Mjema, Seth Nguhu, Thomas Semdoe, Sadiki Yusuph, Alex Mwakibinga, Rodrick Ulomi and Andrea Kimboi. Management of the National Institute for Medical Research, National Malaria Control Program and President’s Office-Regional Administration and Local Government (Regional administrative secretaries of the 13 regions, and district officials, staff from all 100 HFs and Community Health Workers from the 3 community cross sectional regions. Technical and logistics support from the Bill and Melinda Gates Foundation team is highly appreciated. Permission to publish the manuscript was sought and obtained from the Director General of NIMR. ## Footnotes * + co-first author * * co-senior authors * We have updated the analysis to include additional antimalarial resistance mutations. * Received November 7, 2023. * Revision received November 30, 2023. * Accepted November 30, 2023. * © 2023, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## REFERENCES 1. World Health Organization. World malaria report 2022. World Health Organization, 2022. 2. Trape JF. The public health impact of chloroquine resistance in Africa. Am J Trop Med Hyg 2001; 64: 12–7. [Abstract](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo3OiI2NC8xLzEyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTEvMzAvMjAyMy4xMS4wNy4yMzI5ODIwNy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 3. Spring MD, Lin JT, Manning JE, et al. Dihydroartemisinin-piperaquine failure associated with a triple mutant including kelch13 C580Y in Cambodia: an observational cohort study. Lancet Infect Dis 2015; 15: 683–91. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(15)70049-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25877962&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F30%2F2023.11.07.23298207.atom) 4. Programme GM. Report on antimalarial drug efficacy, resistance and response: 10 years of surveillance (2010-2019). 2020; published online Nov 19. [https://www.who.int/publications/i/item/9789240012813](https://www.who.int/publications/i/item/9789240012813) (accessed Oct 30, 2023). 5. Ward KE, Fidock DA, Bridgford JL. Plasmodium falciparum resistance to artemisinin-based combination therapies. Curr Opin Microbiol 2022; 69: 102193. 6. Conrad MD, Asua V, Garg S, et al. Evolution of Partial Resistance to Artemisinins in Malaria Parasites in Uganda. N Engl J Med 2023; 389: 722–32. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2211803)&link_type=DOI) 7. Mihreteab S, Platon L, Berhane A, et al. Increasing Prevalence of Artemisinin-Resistant HRP2-Negative Malaria in Eritrea. N Engl J Med 2023; 389: 1191–202. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2210956&link_type=DOI) 8. Fola AA, Feleke SM, Mohammed H, et al. Plasmodium falciparum resistant to artemisinin and diagnostics have emerged in Ethiopia. Nat Microbiol 2023; 8: 1911–9. 9. Moser KA, Madebe RA, Aydemir O, et al. Describing the current status of Plasmodium falciparum population structure and drug resistance within mainland Tanzania using molecular inversion probes. Mol Ecol 2020; published online Oct 26. DOI:10.1111/mec.15706. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/mec.15706&link_type=DOI) 10. Uwimana A, Legrand E, Stokes BH, et al. Emergence and clonal expansion of in vitro artemisinin-resistant Plasmodium falciparum kelch13 R561H mutant parasites in Rwanda. Nat Med 2020; published online Aug 3. DOI:10.1038/s41591-020-1005-2. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-1005-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32747827&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F30%2F2023.11.07.23298207.atom) 11. Bayih AG, Getnet G, Alemu A, Getie S, Mohon AN, Pillai DR. A Unique Plasmodium falciparum K13 Gene Mutation in Northwest Ethiopia. Am J Trop Med Hyg 2016; 94: 132– 5. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiI5NC8xLzEzMiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzExLzMwLzIwMjMuMTEuMDcuMjMyOTgyMDcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 12. Kirby R, Giesbrecht D, Karema C, et al. Examining the Early Distribution of the Artemisinin-Resistant Plasmodium falciparum kelch13 R561H Mutation in Areas of Higher Transmission in Rwanda. Open Forum Infect Dis 2023; 10: ofad149. 13. Uwimana A, Umulisa N, Venkatesan M, et al. Association of Plasmodium falciparum kelch13 R561H genotypes with delayed parasite clearance in Rwanda: an open-label, single-arm, multicentre, therapeutic efficacy study. Lancet Infect Dis 2021; 21: 1120–8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(21)00142-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33864801&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F30%2F2023.11.07.23298207.atom) 14. WHO | Status report on artemisinin and ACT resistance (April 2017). World Health Organization, 2017 [https://www.who.int/malaria/publications/atoz/artemisinin-resistance-april2017/en/](https://www.who.int/malaria/publications/atoz/artemisinin-resistance-april2017/en/) (accessed Feb 3, 2020). 15. Verity R, Aydemir O, Brazeau NF, et al. The impact of antimalarial resistance on the genetic structure of Plasmodium falciparum in the DRC. Nat Commun 2020; 11: 2107. 16. Rogier E, Battle N, Bakari C, et al. Plasmodium falciparum pfhrp2 and pfhrp3 gene deletions among patients enrolled at 100 health facilities throughout Tanzania: February to July 2021. bioRxiv. 2023; published online Aug 5. DOI:10.1101/2023.07.29.23293322. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1101/2023.07.29.23293322&link_type=DOI) 17. Bwire GM, Ngasala B, Mikomangwa WP, Kilonzi M, Kamuhabwa AAR. Detection of mutations associated with artemisinin resistance at k13-propeller gene and a near complete return of chloroquine susceptible falciparum malaria in Southeast of Tanzania. Sci Rep 2020; 10: 1–7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-020-59121-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F30%2F2023.11.07.23298207.atom) 18. Ishengoma DS, Mmbando BP, Mandara CI, et al. Trends of Plasmodium falciparum prevalence in two communities of Muheza district North-eastern Tanzania: correlation between parasite prevalence, malaria interventions and rainfall in the context of re-emergence of malaria after two decades of progressively declining transmission. Malar J 2018; 17: 252. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-018-2395-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29976204&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F30%2F2023.11.07.23298207.atom) 19. Oyola SO, Ariani CV, Hamilton WL, et al. Whole genome sequencing of Plasmodium falciparum from dried blood spots using selective whole genome amplification. Malar J 2016; 15. DOI:10.1186/s12936-016-1641-7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-016-1641-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27998271&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F30%2F2023.11.07.23298207.atom) 20. Niaré K, Greenhouse B, Bailey JA. An optimized GATK4 pipeline for Plasmodium falciparum whole genome sequencing variant calling and analysis. Malar J 2023; 22: 207. 21. Blasco B, Leroy D, Fidock DA. Antimalarial drug resistance: linking Plasmodium falciparum parasite biology to the clinic. Nat Med 2017; 23: 917–28. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm.4381&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28777791&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F30%2F2023.11.07.23298207.atom) 22. Takala-Harrison S, Jacob CG, Arze C, et al. Independent emergence of artemisinin resistance mutations among Plasmodium falciparum in Southeast Asia. J Infect Dis 2015; 211: 670–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiu491&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25180241&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F30%2F2023.11.07.23298207.atom) 23. Katalyeba PO. Rural - Rural Movements in East Africa: a Study of Rwandan and Southern Ugandans in Karagwe District, Tanzania. 2002. [http://researcharchive.vuw.ac.nz/handle/10063/461](http://researcharchive.vuw.ac.nz/handle/10063/461) (accessed Oct 28, 2023). 24. Gesase S, Gosling RD, Hashim R, et al. High resistance of Plasmodium falciparum to sulphadoxine/pyrimethamine in northern Tanzania and the emergence of dhps resistance mutation at Codon 581. PLoS One 2009; 4: e4569. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0004569&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19238219&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F30%2F2023.11.07.23298207.atom) 25. Ngondi JM, Ishengoma DS, Doctor SM, et al. Surveillance for sulfadoxine-pyrimethamine resistant malaria parasites in the Lake and Southern Zones, Tanzania, using pooling and next-generation sequencing. Malar J 2017; 16: 236. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-017-1886-9&link_type=DOI) 26. Deutsch-Feldman M, Aydemir O, Carrel M, et al. The changing landscape of Plasmodium falciparum drug resistance in the Democratic Republic of Congo. BMC Infect Dis 2019; 19: 872. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12879-019-4523-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F30%2F2023.11.07.23298207.atom) 27. Taylor SM, Antonia AL, Harrington WE, et al. Independent lineages of highly sulfadoxine-resistant Plasmodium falciparum haplotypes, eastern Africa. Emerg Infect Dis 2014; 20: 1140–8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2007.131720&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24960247&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F30%2F2023.11.07.23298207.atom) 28. Hyde JE. Antifolate resistance in Africa and the 164-dollar question. Trans R Soc Trop Med Hyg 2008; 102. DOI:10.1016/j.trstmh.2008.01.003. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.trstmh.2008.01.003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18316104&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F30%2F2023.11.07.23298207.atom) 29. Chico RM, Cano J, Ariti C, et al. Influence of malaria transmission intensity and the 581G mutation on the efficacy of intermittent preventive treatment in pregnancy: systematic review and meta-analysis. Trop Med Int Health 2015; 20: 1621–33. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F30%2F2023.11.07.23298207.atom)