Skip to main content

Temporal evolution of sulfadoxine-pyrimethamine resistance genotypes and genetic diversity in response to a decade of increased interventions against Plasmodium falciparum in northern Ghana

Abstract

Background

Anti-malarial drug resistance remains a key concern for the global fight against malaria. In Ghana sulfadoxine-pyrimethamine (SP) is used for intermittent preventive treatment of malaria in pregnancy and combined with amodiaquine for Seasonal Malaria Chemoprevention (SMC) during the high malaria season. Thus, surveillance of molecular markers of SP resistance is important to guide decision-making for these interventions in Ghana.

Methods

A total of 4469 samples from uncomplicated malaria patients collected from 2009 to 2018 was submitted to the Wellcome Trust Sanger Institute, UK for DNA sequencing using MiSeq. Genotypes were successfully translated into haplotypes in 2694 and 846 mono infections respectively for pfdhfr and pfdhps genes and the combined pfhdfr/pfdhps genes across all years.

Results

At the pfdhfr locus, a consistently high (> 60%) prevalence of parasites carrying triple mutants (IRNI) were detected from 2009 to 2018. Two double mutant haplotypes (NRNI and ICNI) were found, with haplotype NRNI having a much higher prevalence (average 13.8%) than ICNI (average 3.2%) across all years. Six pfdhps haplotypes were detected. Of these, prevalence of five fluctuated in a downward trend over time from 2009 to 2018, except a pfdhps double mutant (AGKAA), which increased consistently from 2.5% in 2009 to 78.2% in 2018. Across both genes, pfdhfr/pfdhps combined triple (NRNI + AAKAA) mutants were only detected in 2009, 2014, 2015 and 2018, prevalence of which fluctuated between 3.5 and 5.5%. The combined quadruple (IRNI + AAKAA) genotype increased in prevalence from 19.3% in 2009 to 87.5% in 2011 before fluctuating downwards to 19.6% in 2018 with an average prevalence of 37.4% within the nine years. Prevalence of parasites carrying the quintuple (IRNI + AGKAA or SGEAA) mutant haplotypes, which are highly refractory to SP increased over time from 14.0% in 2009 to 89.0% in 2016 before decreasing to 78.9 and 76.6% in 2017 and 2018 respectively. Though quintuple mutants are rising in prevalence in both malaria seasons, together these combined genotypes vary significantly within season but not between seasons.

Conclusions

Despite high prevalence of pfdhfr triple mutants and combined pfdhfr/pfdhps quadruple and quintuple mutants in this setting SP may still be efficacious. These findings are significant as they highlight the need to continuously monitor SP resistance, particularly using deep targeted sequencing to ascertain changing resistance patterns.

Background

Malaria disproportionately leads to childhood morbidity and mortality in sub-Saharan Africa (sSA). An estimated 24 million children were infected with Plasmodium falciparum in 2018 in sSA with Ghana being one of top ten countries in Africa with the highest absolute increases in cases of malaria in 2018 compared to the previous year [1]. Due to persistent high malaria transmission, the Ghana National Malaria Control Programme (NMCP) prioritized the northern regions for high-impact interventions, such as indoor residual spraying (IRS), seasonal malaria chemoprevention (SMC) among children under 5 years old [2], increased coverage of long-lasting insecticidal bed nets (LLINs) and intermittent preventive treatment in pregnancy (IPTp) [3, 4]. The World Health Organization (WHO) recommends use of SP plus amodiaquine (SP-AQ) for SMC in areas of high seasonal malaria transmission in the Sahel sub-region of Africa [5]. Therefore, despite the withdrawal of SP in 2005 as a first-line anti-malarial in Ghana [6], it remains a key component in interventions targeting vulnerable groups such as pregnant women and young children. However, the success of targeted interventions such as, IPTp and SMC, is largely dependent on the population prevalence of bifunctional dihydrofolate reductase-thymidylate synthase (dhfr) and dihydropteroate synthetase (dhps) SP resistance mutations [7,8,9].

Resistance of P. falciparum to SP is caused by well documented single nucleotide polymorphisms (SNPs) in the pfdhfr and pfdhps genes that encode enzymes in the folate metabolism pathway and are targeted by pyrimethamine and sulfadoxine, respectively [10, 11]. A change from wild-type Ser108 to Asn108 (S108N) in pfdhfr is associated with low level pyrimethamine tolerance both in vitro and in vivo [10,11,12]. This basal S108N amino acid substitution is characterized to have a tenfold increased risk of SP therapeutic failure [13]. The progressive accumulation of other mutations that result in altered amino acid substitutions including pfdhfr-C50R, pfdhfr-N51I, pfdhfr-C59R, and pfdhfr-I164L that can result in higher levels of resistance to pyrimethamine and diminish the efficacy of SP in vivo [14]. Parasites carrying multiple mutations in the pfdhfr gene have evolved independently in different populations throughout the malaria endemic world [15], and these haplotypes are often associated with higher levels of resistance as compared to single mutant genotypes. Typically, the most prevalent multi-locus genotype found in areas of high SP resistance in sSA is the triple pfdhfr mutant genotype (pfdhfr-51I/59R/108N/I164) [16, 17]. A quadruple pfdhfr mutant first detected in southeast Asia in the late 1980s, includes an additional I164L substitution [18]. Parasites carrying the quadruple haplotype are associated with high SP treatment failure [19]. Thus far, parasites carrying this allele have been only detected in parts of East Africa [16, 20], but not West Africa [11, 12].

Parasites carrying the triple pfdhr mutant (108N/51I/59R/I164) alleles and additional pfdhps double mutant (G437/540E) alleles have been associated with strong resistance to SP [21]. A further mutation at codon 581 of pfdhps results in a triple mutant allele (A437G/K540E/A581G), which in combination the pfdhfr triple mutant confers a high rate of SP treatment failure [22, 23]. The occurrence of the K540E mutation is a proxy marker for a variant quintuple mutant genotype containing both the pfdhps double mutant (A437G/K540E) and the pfdhfr triple mutant genotypes, which is also highly correlated with SP treatment failure in children [9, 21, 22].

The objectives of the present study were to determine the temporal trends in prevalence of multi-locus SP resistance markers in P. falciparum isolates collected from clinical sources between 2009 and 2018, a period spanning major interventions for reducing the burden of malaria and disrupting transmission including IRS, SMC and IPTp. This study also explores the seasonal risk of carrying SP resistant genotypes in northern Ghana.

Methods

Study site

The study was conducted in the Kassena-Nankana East Municipality and Kassena-Nankana West District, two adjoining administrative areas in the Upper East Region of northern Ghana, here referred to as KNDs. The KNDs cover a total area of about 1675km2 and lie between latitude 10.30′ and 11.10′ North and longitude 1.1′ West close to the Ghana-Burkina Faso border (Fig. 1). Average rainfall is estimated at about 1300 mm, but this has been consistently dwindling in recent years [24]. The KNDs population is estimated at 160,000 inhabitants with a population density of 91.5 per square kilometre [25]. The proportion of children under 5 years has been estimated at about 11.0% and there are approximately 4000 births per year. The KNDs population is under continuous demographic surveillance by the Navrongo Health and Socio-demographic Surveillance System. There is one main hospital, the Navrongo War Memorial Hospital (NWMH) that serves as a referral unit for the population, nine satellite clinics and several additional community-based health planning and services compounds strategically located to increase access to basic primary health care [26]. Malaria transmission is intense with seasonal fluctuations largely dependent on the rainfall pattern [24, 25]. The high transmission season coincides with the short rainy season from July to November, while the low season occurs during the dry months of December to May. The main malaria vector is Anopheles gambiae and the annual entomological inoculation rate (EIR) was estimated to vary between 1132 and 157 infective bites/person/year in the KNDs [27].

Fig. 1
figure1

Location of the Upper East Region in Ghana, The Kassena-Nankana Districts in the Upper East Region and Kassena-Nankana East and West Districts showing location of Health facilities from which samples were collected. CHPS; Community Health and Planning Services

Ethics statement

Scientific and ethical clearance was obtained from the Navrongo Health Research Centre Institutional Review Board (# NHRCIRB203). Informed consent was obtained from all individuals or their parents/guardians prior to enrollment into each study.

Study procedure

All-age patients with fever (or history of fever within the past 24hrs) reporting to the NWMH, and four clinics in the KNDs were screened for malaria from January to December of each year from 2009 to 2018, except 2012 when no survey was conducted. Plasmodium falciparum Histidine Rich Protein 2-based rapid diagnostic test (RDT)[CareStart™ malaria Pf (HRP2), Access Bio, NJ, USA] was used. Informed consent was documented for RDT positive patients and those willing to participate were enrolled. Four 50 µL dried blood spots (DBS) and malaria smears were prepared from a finger prick. Plasmodium falciparum density was determined by microscopy and asexual parasites scored against 200 white blood cells (WBCs) and smears were considered negative after examining 100 fields without detecting asexual or sexual parasites. Parasite density per 200 WBCs was converted to density per microlitre assuming 8000 WBCs.

Parasite DNA was extracted from DBS samples using the QIAamp DNA Midi kit (Qiagen, UK) as per the manufacturer's protocol.

Genotyping of Drug resistance markers

Samples between the 2009 and 2013 studies were submitted to the Malaria Genomic Epidemiology Network (MalariaGEN; http://www.malariagen.net/) at the Wellcome Sanger Institute, UK for whole genome sequencing as previously described [28]. Genotypes were assigned for variants with a read depth ≥ 5. Genotypes for the selected resistance genes were then extracted, translated and annotated into known SP resistance markers. For samples collected between 2014 and 2018, amplicon sequencing of SP resistance markers was performed by MalariaGEN using an amplicon-sequencing protocol (https://www.medrxiv.org/content/10.1101/2020.07.23.20159624v1 and https://www.malariagen.net/resource/29). In brief, locus-specific multiplexed primers were designed using MPprimer software [29] with some modifications. Primers were designed to generate amplicons each of between 190 and 250 bp and were assigned to one of three pools. A two-step PCR protocol was used for each pool to first amplify the regions of interest in the parasite genome, followed by a second PCR to incorporate sequencing and unique sample-level and primer-pool multiplexing adapters. Multiple samples were sequenced on a single MiSeq lane combining all 3 pools of the PCR amplicons. Sequenced samples were de-plexed using the unique multiplexing adapter IDs and the individual sample-level CRAM files were aligned to a modified amplicon P. falciparum reference genome. Genotyping of bi-allelic SNPs was undertaken using bcftools as well as custom scripts to determine genotypes which were translated and annotated into known drug resistance haplotypes (https://www.medrxiv.org/content/10.1101/2020.07.23.20159624v1 and https://www.malariagen.net/resource/29).

Within-host diversity estimation

Within-host diversity of (complexity) infections was determined using either the genome-wide Fws metric in the 2009 to 2013 sequence data or THE REAL McCOIL to estimate infection complexity in the 2014 to 2018 amplicon data. The McCOIL’s categorical method, which uses likelihood estimation to determine the number of distinct parasite genome-wide haplotypes (strains) within each sample was implemented [30].

Haplotype reconstruction

Owing to the high within host diversity, in the 2009 to 2013 collections, genotypes were translated into haplotypes for each isolate based on the major allele for pfdhfr and pfdhps loci. The amplicon data from 2014 to 2018 however, enabled the genotyping of pure or mixed (heterozygous) samples and the reconstruction of haplotypes for both pure and heterozygous samples, accounting for both major and minor clones in mixed infections. Only samples with Fws > 0.95 and MOI = 1 were used for construction of haplotypes across both genes to obtain pfdhfr/pfdhps combined genotypes.

Data and statistical analysis

To determine the temporal trends of multi-locus SP resistance markers, a year-on-year fluctuation in SP resistance markers was examined. To further explore the effect of malaria interventions, the data was grouped into three windows to reflect periods of major interventions-2009–2011(pre-IRS and SMC), 2013–2015 (IRS) and 2016–2018 (SMC). However, it is worth noting that IPTp and bednets were implemented across all study periods. The following analysis approaches were explored: a descriptive analysis of temporal trends in proportion of isolates with SP resistant marker haplotypes, mixed species infections and within-host diversity of infections. The χ2 test was used to compare categorical variables among groups. The non-parametric Kruskal–Wallis test was used for group comparisons as appropriate for the distribution. All statistical tests were two-tailed and statistical significance set at p < 0.05. Data were analysed using the open source statistical software R version 3.4.1.

Results

A total 4469 patients with uncomplicated malaria were enrolled into this study in the years 2009–2018. Of these patients, 2694 [(2009; n = 120); (2010; n = 184); (2011; n = 55); (2013; n = 33); (2014; n = 323); (2015; n = 178); (2016; n = 920); (2017; n = 560); (2018; n = 321)] were available for genetic analysis of pfdhfr and pfdhps loci. A total of 846 samples [(2009; n = 57); (2010; n = 81); (2011; n = 24); (2013; n = 22); (2014; n = 73); (2015; n = 70); (2016; n = 237); (2017; n = 175); (2018; n = 107)] passed as mono-infections (Fws > 0.95 or MOI = 1) and were included in genetic analysis across both genes (Table 1).

Table 1 Temporal trends in prevalence of SP resistant haplotypes from 2009 to 2018

Clinical and demographic characteristics

In the years 2009 to 2018, participants were aged between 3.7 to 21.5 years. Overall, there were more female patients (56.3%) than males in the study period 2009 to 2018. The results show a declining trend in parasite density from 2009 to 2018. Trends in monthly mean parasite densities from 2009 to 2018 show a similar unimodal peak during the high malaria season (rainy months of July to November) (Fig. 2 and Additional file 1: Table S1) (Additional file 2: Fig S1).

Fig. 2
figure2

Monthly trends in mean P. falciparum parasite density from 2009 to 2018

The median parasite density decreased fourfold from 38,280 (IQR; 62,320) in the 2009–2011 study period to 10,600 (IQR; 52,780) in the 2013–2015 study period but increased to 17,680 (IQR; 49,680) in the 2016–2018 period with a significant decline in the distribution of parasite density during peak transmission from the 2009–2011 study period and the two subsequent surveys (2009–2011 vs 2013–2015 (P < 0.001) or 2016–2018, P < 0.001) (Additional file 3: Fig. S2). Overall, the majority (> 97.8%) of infections were P. falciparum mono infections with 2.2% being co-infections with Plasmodium malariae as determined from amplicon sequencing data.

Within-host diversity

Within-host diversity in the earlier studies from 2009 to 2013 ranged from 0.003 to 0.864 (mean 0.254 ± 0.247 SD), estimated by the inbreeding co-efficient Fws. Overall, 79.5% (198/249) had Fws ≤ 0.5 in the 2009–2013 studies, indicating a high proportion of multi-genomic infections with high potential for outcrossing. Interestingly, overall, in the aggregated data, infections were significantly more complex in the low malaria transmission season (P = 0.022) (Additional file 4: Fig. S3A) but the year on year showed significantly higher diversity only in 2011 (P < 0.001) (Additional file 5: Fig. S4A). In the latter studies from 2014 to 2018 infection complexity ranged from one to three clones, except in 2015 when up to four clones were detected (Additional file 5: Fig. S4B). A higher proportion of monoclonal infections was observed in the 2016–2018 studies in both seasons.

Temporal trends in pfdhfr, pfdhps and the combined pfdhfr/pfdhps genotypes

The temporal trends in the pfdhfr, pfdhps haplotypes and pfdhfr/pfdhps combined haplotypes were examined on a year-on-year analysis of fluctuations in the prevalence of circulating haplotypes in these genes as shown in Fig. 3. At the Pfdhfr locus, year-on-year prevalence of triple (IRNI) Pfdhfr mutants associated with pyrimethamine resistance is on a general decline though still high (> 60%) (Fig. 3a). Two double mutant haplotypes (NRNI and ICNI) were detected, with haplotype NRNI having a higher year-on-year prevalence compared to haplotype ICNI, except in 2013 when their prevalence estimates were similar (Fig. 3a). The annual prevalence of NRNI fluctuated with time from 7.5 in 2009 to 25% in 2018 whilst prevalence of ICNI fluctuated over time from 3.3% in 2009 to 2.2% in 2018. In general, seasonality analysis showed no significant differences in the prevalence of pfdhfr mutants between the wet (high) and dry (low) seasons (Fig. Additional file 6: S5). However, the pfdhfr triple mutants increased in prevalence (> 75%) from 2015 to 2016 (when SMC was introduced) during the dry season and fluctuated sharply to about 50% and 59% in 2017 and 2018 respectively (Additional file 6: Fig. S5A). In the dry season comparing the double mutants (NRNI and ICNI), NRNI had about a 20-fold higher prevalence (≤ 20%) than the ICNI mutant (< 1%) between 2011 to 2015 when prevalence of NRNI decreased and that of ICNI increased slightly to narrow their prevalence difference in 2016, after which NRNI increased sharply to > 30% whilst ICNI decreased gradually to < 2% in 2018. On the other hand, during the wet season, whilst pfdhfr triple mutants declined similarly as observed in the dry season, albeit a gradual decrease after 2016, the double mutant (NRNI) fluctuated more rapidly at between 3 to 20% from 2009 to 2016 when prevalence of the two mutants equalized to about 2.5% (Additional file 6: Fig. S5B). Between 2016 to 2018, NRNI increased sharply to about 20% whilst ICNI decreased to below 1%. The I164L amino acid substitution associated with high SP treatment failure in some settings was not detected in the current study.

Fig. 3
figure3

Temporal trends in haplotype prevalence in SP resistance genes from 2009 to 2018. Panel A: pfdhfr double (ICNI or NRNI),triple (IRNI) and sensitive (NCSI), Panel B: pfdhps triple ([A/F/Y] AKAS or [A/S]GKAS or [S/AGEAA), double (C/S/A/GKAA), single (AAKAA) and fully sensitive (SAKAA). Panel C: pfdhr/pfhps combined haplotypes as defined in Table 1. Mutant amino acid substitutions are shown in bold

At the pfdhps locus, six different haplotypes were observed. Of these haplotypes 5 declined in prevalence over time from between 12.5 and 56% to less than 12.5% with the exception of the pfdhps double mutant (C/S/AGKAA), which fluctuated between 2009 and 2011 and increased consistently from about 9% in 2011 to 71% in 2016, decreased slightly to 69.5% and then increased to about 78% in 2018 (Table 1 and Fig. 3b). The results from seasonality analysis showed a similar trend with the same 5 haplotypes occurring in both seasons and showing a general decline from 2009 to 2018, particularly after 2014 with haplotype ((A/S)GKAS) having higher prevalence in dry season compared to the wet season. However, the double mutant (C/S/AGKAA), increased in prevalence sharply in both seasons rising from about 33% in 2014 to a peak of about 84% during the dry (low) season in 2016 and then decreasing to 69% in 2018. The double mutant also rose sharply from 46% in 2014 to peak at 85% 2016 and decrease to 82% in 2018 (Additional file 7: Fig. S6B).

Figure 3c shows the yearly fluctuations in the combined pfdhfr/pfdhps genotypes. Of the three haplotypes observed (i.e. triple, quadruple and quintuple), the prevalence of triple (NRNI + AAKAA) haplotype was generally low (0 to 6%) from 2009 to 2018. The quadruple (IRNI + AAKAA) mutants increased in frequency from 2009 (44%) to 2011 (88%) [pre-IRS] and decreased rapidly from 2011 to about 10% in 2016 before increasing sharply to about 21% in 2017 and 2018 (Table 1 and Fig. 3c). However, prevalence of parasites carrying quintuple mutations (IRNI + AGKAA or SGEAA), which are highly refractory to SP initially decreased from 2009 (53%) to 2011 (13%) and then increased rapidly over time to 77% in 2018 with a light decrease in prevalence between 2016 (89%) and 2017 (79%) (Fig. 3c). No sextuple mutants were detected from 2009 to 2018. The pfdhfr/pfdhps combined genotypes did not differ in prevalence between seasons (p > 0.05) but within season prevalence estimates were statistically significant (Wilcoxon, dry season; p < 0.001 & wet season; p = 0.006) (Additional file 8: Fig. S7A). Also, only the prevalence of quintuple haplotype was significantly different between pre-SMC (2009–2015) and post-SMC (2016–2018) (Additional file 8: Fig. S7B). However, the year-on-year dry season prevalence of the quadruple mutants increased sharply from 2009 (9.8%) to 2011 (65%) but rapidly decreased from 2013 (60%) to 2016 (10%) before increasing slightly to 19.8% in 2017 and falling back to 10% in 2018 (Additional file 9: Fig. S8A). The wet season dynamics showed marked fluctuation in the prevalence of quadruple mutants from slightly over 10% in 2009 to 22% in 2011, and about 18% in 2014 to 14% in 2015 before increasing from about 2.5 to 9% in 2017 and 2018 respectively. No quadruple mutants were detected in 2013 and 2016 during the wet season (Fig. S8B). On the one hand, during dry season, quintuple mutant prevalence fluctuated markedly from 10% (2009), 50% (2010) to 12% (2011), then increased gradually through 22% (2013), 28% (2014) to peak at about 48% prevalence in 2015 and 2016 before fluctuating down to 32% in 2017 and 38% in 2018. However, during the wet season quintuple mutants initially found at 5% prevalence were undetected consecutively in 2010 and 2011 but re-emerged in 2013 at 18% prevalence and increased consistently to peak at 46% in 2017 before dropping to a prevalence of 39% in 2018 (Additional file 9: Fig. S8B).

Discussion

The current study utilized data from Illumina whole genome sequencing and deep amplicon sequencing to determine SNP haplotypes for SP resistance markers following a stringent SNP calling criteria. This enabled reliable haplotype reconstruction in both single and multi-genome infections by accounting for minor clones and avoiding the ambiguity posed when using the standard methodology or individual genotyping assays on mixed infections as observed in this study [31]. A significant decline in the monthly mean prevalence of parasitaemia was observed between the 2009–2011 studies and the subsequent studies in 2013–2015 or 2016–2018 during the high malaria transmission season. The 2009 to 2011 studies were conducted prior to the deployment of strategic interventions, such as IRS and SMC in 2014 and 2016, respectively. Therefore, the marked decline in prevalence of parasitaemia in the post IRS study period 2016–2018 may be due to the direct impact of this intervention aimed at disrupting transmission. The combined effects of IRS at the beginning of the high transmission season and the deployment of SMC (i.e. 4 rounds) during high transmission may explain these observations. However, the observed significantly increased mean parasitaemia post-SMC (2013–2015 vs 2016–2018) corroborates our observed shift in malaria morbidity to higher age groups, which could be attributed to mass treatment campaigns such as SMC that protect young children from clinical malaria. This highlights unintended epidemiological consequences of such targeted mass treatment campaigns and suggest that older children may also benefit from these campaigns. The high within-host diversity observed in the 2009 to 2013 studies with over 70% of infections having Fws < 0.5 is an indication of transmission levels in the era in these settings that pre-date major interventions such as IRS aimed at transmission reduction. The present study highlights the impact of the combined effects of IRS, long-lasting bed nets and SMC, in driving down multi-clonal infections with about 40% of infections harbouring two clones per year from 2014 to 2018, and less than 5% of infections have 3–4 clones between 2014 and 2018. This is consistent with a drop in annual EIR estimates from previous > 250 infective bites/person/year [32] to about 50 bites/person/year with wide variation particularly during the high malaria season in these settings [27].

In the absence of an alternative drug to SP for IPTp and other intermittent preventive treatment programmes in infants such as SMC, monitoring of resistance markers in endemic populations is crucial. In characterizing SP-associated resistance loci in the current study settings, a persistently high prevalence of the triple-mutant (51I + 59R + 108 N + I164) pfdhfr haplotype from 2009 to 2018 was found. This haplotype, observed at > 60% prevalence year-on-year from 2009 to 2018, has been previously demonstrated to contribute significantly to SP treatment failure [22]. The fairly stable double mutant haplotypes (NRNI and ICNI) appear to differ in fitness with parasites carrying the NRNI haplotypes being more prevalent. No marked seasonality was found in distribution of these haplotypes year-on-year but the prevalence of NRNI was higher in dry (low) season suggesting the acquisition of the haplotype alleles by the infection reservoir during high transmission and subsequent seeding of infections in the dry season. Well-designed community studies are required to ascertain the seasonality of these haplotypes. Of note is the complete absence of the pfdhfr-164L allele, which has been reported in parts of East Africa [33], but thus far not in West Africa. The counterpart SP-resistance associated gene locus pfdhps also had several mutant haplotypes that have been previously characterized to confer high resistance phenotypes. Six different haplotypes were found in this gene, with two (AGEAA and AGKAS) containing markers of high-grade resistance with AGKAS increasing over time and within seasons in this setting. The haplotype AGEAA, which has the K540E mutation previously shown to cause increased SP resistance is observed at < 5% prevalence year-on-year.

The declining year-on-year fluctuations of pfdhfr and pfdhps haplotypes detected suggest a slow expansion of these haplotypes in this study settings. This may be driven by high recombination rates, particularly during the high malaria season. The higher complexity of infections observed in the current study attest to high transmission intensity in these settings, which is ripe for frequent recombination resulting in breakdown of these long haplotypes.

Also, temporal trends in prevalence of pfdhfr/pfdhps combined genotypes show a fluctuating trend of quintuple and quadruple mutants respectively until 2016 when SMC was deployed in our study communities. Since then, the yearly prevalence of both combined haplotypes began to fluctuate in opposing directions with a slight decrease in quadruple mutant prevalence and a slight increase in quintuple prevalence. It is interesting that the distribution of these combined genotypes was not significant between the wet and dry season but rather within each season (Additional file 8: Fig. S7A). This may be attributed to differences in transmissibility within each season. Only quintuple mutants were significantly different between pre-and post-SMC periods, suggesting a possible additional selection pressure imposed by the SMC campaign. This is supported by the observed increased year-on-year prevalence of quintuple mutants in the wet season (Additional file 9: Fig. S8). The overall survival of parasites carrying quintuple mutants all year round may be attributed to increased fitness in the presence of routine anti-malarials and from prophylaxis under SMC. Previous studies have suggested that childhood intervention campaigns such as SMC are likely to promote the spread of resistance in high transmission settings than other adult interventions, such as IPTp-SP [34]. The SMC drugs are given at the beginning of the rainy season (high transmission) in four rounds from July to October. These trends provide support for targeted interventions such as SMC to take place during the high transmission season in order to maximize the long acting effect of AQ partner drug in this intervention to kill these parasites that may have a higher fitness in the presence of SP. However, sentinel data from 2015 to 2017 showed PCR-corrected AS-AQ treatment efficacy of 98.2% in these settings [35], and there is also accumulating evidence of rising prevalence of AQ resistance markers in this population (data to be described elsewhere). Previous studies show that IPTp-SP may remain efficacious in the phase of circulating high grade resistance genotypes but additional mutations, such as 581G, 540E and 613S see the benefits of IPTp-SP begin to decline. Therefore, these mutations may have distinct impacts on different malaria interventions (eg SMC, IPTp among others) and the genotype–phenotype relations may not be perfect. Other factors, such as amplification of parasite GTP cyclohydrolase 1 (GCH1) gene that mediate de novo folate biosynthesis [36], host immunity and metabolism may limit these genotype–phenotype interactions.

Overall, barring any sampling bias due to the passive case detection approach in this study, these findings signal a real threat to obtaining the full benefit of key SP-based interventions, such as IPTp and SMC in this setting. However, it is worth noting that the prevalence levels obtained in the current study fall below the WHO limits of > 50% prevalence of 540E to signal failure of SP interventions at a particular setting [22].

Preventing malaria in vulnerable groups such as pregnant women and young children comes with enormous benefits, yet there is currently no effective alternative to SP for use in interventions targeting these groups. Hence, the need to continue to monitor these resistance markers in the general population and target populations to provide evidence-based guidance for the implementation of SP intervention programmes in Ghana. This study, however, is limited to the catchment area of the Navrongo War Memorial Hospital, which is only one out of the 10 NMCP sentinel sites for monitoring anti-malarial drug efficacy across three transmission zones in Ghana [37]. Therefore, there is the need to assess the prevalence of these molecular markers across all 10 sentinel sites with ecological variance across Ghana using targeted sequencing methods to make the data more valuable to the National Malaria Control Programme.

Conclusion

This study highlights high prevalence of parasites carrying pfdhfr triple mutants and rising prevalence of combined pfdhfr/pfdhps quintuple mutants in this setting. Thus, indicating the need to continuously monitor parasite response to SP across all three transmission zones in Ghana.

Availability of the data

Raw data is available on European Nucleotide Archive (ENA). All data generated is presented within the manuscript and its supplementary files.

Abbreviations

WHO:

World Health Organization

SP:

Sulfadoxine-pyrimethamine

RDT:

Rapid diagnostic test

NMCP:

National Malaria Control Programme

SMC:

Seasonal Malaria Chemoprevention

IPTp:

Intermittent Preventive Treatment in Pregnancy

LLINs:

Long-lasting Insecticidal Nets

IRS:

Indoor Residual Spraying

EIR:

Entomological Inoculation Rate

SNP:

Single nucleotide polymorphism

PCR:

Polymerase chain reaction

Pfdhfr:

Plasmodium falciparum dihydrofolate reductase

Pfdhps:

Plasmodium falciparum dihydropteroate synthetase

References

  1. 1.

    WHO. World malaria report. Geneva: World Health Organization; 2019.

    Google Scholar 

  2. 2.

    Nonvignon J, Aryeetey GC, Issah S, Ansah P, Malm KL, Ofosu W, et al. Cost-effectiveness of seasonal malaria chemoprevention in upper west region of Ghana. Malar J. 2016;15:367.

    Article  Google Scholar 

  3. 3.

    Steketee RW, Nahlen BL, Parise ME, Menendez C. The burden of malaria in pregnancy in malaria-endemic areas. Am J Trop Med Hyg. 2001;64:28–35.

    CAS  Article  Google Scholar 

  4. 4.

    Menéndez C, Bardají A, Sigauque B, Sanz S, Aponte JJ, Mabunda S, et al. Malaria prevention with IPTp during pregnancy reduces neonatal mortality. PLoS ONE. 2010;5:e9438.

    Article  Google Scholar 

  5. 5.

    WHO. Policy recommendation: Seasonal malaria chemoprevention (SMC) for Plasmodium falciparum malaria control in highly seasonal transmission areas of the Sahel sub-region in Africa. Geneva: World Health Organization; 2015.

    Google Scholar 

  6. 6.

    Duah NO, Matrevi SA, de Souza DK, Binnah DD, Tamakloe MM, Opoku VS, et al. Increased pfmdr1 gene copy number and the decline in pfcrt and pfmdr1 resistance alleles in Ghanaian Plasmodium falciparum isolates after the change of anti-malarial drug treatment policy. Malar J. 2013;12:377.

    Article  Google Scholar 

  7. 7.

    Iriemenam NC, Shah M, Gatei W, van Eijk AM, Ayisi J, Kariuki S, et al. Temporal trends of sulphadoxine-pyrimethamine (SP) drug-resistance molecular markers in Plasmodium falciparum parasites from pregnant women in western Kenya. Malar J. 2012;11:134.

    CAS  Article  Google Scholar 

  8. 8.

    Venkatesan M, Alifrangis M, Roper C, Plowe CV. Monitoring antifolate resistance in intermittent preventive therapy for malaria. Trends Parasitol. 2013;29:497–504.

    CAS  Article  Google Scholar 

  9. 9.

    Kublin JG, Dzinjalamala FK, Kamwendo DD, Malkin EM, Cortese JF, Martino LM, et al. Molecular markers for failure of sulfadoxine-pyrimethamine and chlorproguanil-dapsone treatment of Plasmodium falciparum malaria. J Infect Dis. 2002;185:380–8.

    CAS  Article  Google Scholar 

  10. 10.

    Roper C, Pearce R, Bredenkamp B, Gumede J, Drakeley C, Mosha F, et al. Antifolate antimalarial resistance in southeast Africa: a population-based analysis. Lancet. 2003;361:1174–81.

    Article  Google Scholar 

  11. 11.

    McCollum AM, Poe AC, Hamel M, Huber C, Zhou Z, Shi YP, et al. Antifolate resistance in Plasmodium falciparum: multiple origins and identification of novel dhfr alleles. J Infect Dis. 2006;194:189–97.

    CAS  Article  Google Scholar 

  12. 12.

    Hastings MD, Bates SJ, Blackstone EA, Monks SM, Mutabingwa TK, Sibley CH. Highly pyrimethamine-resistant alleles of dihydrofolate reductase in isolates of Plasmodium falciparum from Tanzania. Trans R Soc Trop Med Hyg. 2002;96:674–6.

    CAS  Article  Google Scholar 

  13. 13.

    Lozovsky ER, Chookajorn T, Brown KM, Imwong M, Shaw PJ, Kamchonwongpaisan S, et al. Stepwise acquisition of pyrimethamine resistance in the malaria parasite. Proc Natl Acad Sci USA. 2009;106:12025–30.

    CAS  Article  Google Scholar 

  14. 14.

    Mita T, Ohashi J, Venkatesan M, Marma ASP, Nakamura M, Plowe CV, et al. Ordered accumulation of mutations conferring resistance to sulfadoxine-pyrimethamine in the Plasmodium falciparum parasite. J Infect Dis. 2014;209:130–9.

    CAS  Article  Google Scholar 

  15. 15.

    Roper C, Pearce R, Nair S, Sharp B, Nosten F, Anderson T. Intercontinental spread of pyrimethamine-resistant malaria. Science. 2004;305:1124.

    CAS  Article  Google Scholar 

  16. 16.

    Naidoo I, Roper C. Mapping ‘partially resistant’, ‘fully resistant’, and ‘super resistant’ malaria. Trends Parasitol. 2013;29:505–15.

    Article  Google Scholar 

  17. 17.

    Tinto H, Guekoun L, Zongo I, Guiguemdé RT, D’Alessandro U, Ouédraogo JB. Chloroquine-resistance molecular markers (Pfcrt T76 and Pfmdr-1 Y86) and amodiaquine resistance in Burkina Faso. Trop Med Int Health. 2008;13:238–40.

    CAS  Article  Google Scholar 

  18. 18.

    Sirawaraporn W, Sathitkul T, Sirawaraporn R, Yuthavong Y, Santi DV. Antifolate-resistant mutants of Plasmodium falciparum dihydrofolate reductase. Proc Natl Acad Sci USA. 1997;94:1124–9.

    CAS  Article  Google Scholar 

  19. 19.

    Hyde JE. Antifolate resistance in Africa and the 164-dollar question. Trans R Soc Trop Med Hyg. 2008;102:301–3.

    Article  Google Scholar 

  20. 20.

    Maiga OM, Kayentao K, Traoré BT, Djimde A, Traoré B, Diallo M, et al. Superiority of 3 over 2 doses of intermittent preventive treatment with sulfadoxine-pyrimethamine for the prevention of malaria during pregnancy in Mali: a randomized controlled trial. Clin Infect Dis. 2011;53:215–23.

    CAS  Article  Google Scholar 

  21. 21.

    Omar SA, Adagu IS, Gump DW, Ndaru NP, Warhurst DC. Plasmodium falciparum in Kenya: high prevalence of drug-resistance-associated polymorphisms in hospital admissions with severe malaria in an epidemic area. Ann Trop Med Parasitol. 2001;95:661–9.

    CAS  Article  Google Scholar 

  22. 22.

    Naidoo I, Roper C. Drug resistance maps to guide intermittent preventive treatment of malaria in African infants. Parasitology. 2011;138:1469–79.

    Article  Google Scholar 

  23. 23.

    Matondo SI, Temba GS, Kavishe AA, Kauki JS, Kalinga A, van Zwetselaar M, et al. High levels of sulphadoxine-pyrimethamine resistance Pfdhfr-Pfdhps quintuple mutations: a cross sectional survey of six regions in Tanzania. Malar J. 2014;13:152.

    Article  Google Scholar 

  24. 24.

    Koram KA, Owusu-Agyei S, Fryauff DJ, Anto F, Atuguba F, Hodgson A, et al. Seasonal profiles of malaria infection, anaemia, and bed net use among age groups and communities in northern Ghana. Trop Med Int Health. 2003;8:793–802.

    Article  Google Scholar 

  25. 25.

    Oduro AR, Wak G, Azongo D, Debpuur C, Wontuo P, Kondayire F, et al. Profile of the Navrongo Health and Demographic Surveillance System. Int J Epidemiol. 2012;41:968–76.

    Article  Google Scholar 

  26. 26.

    Dalaba MA, Akweongo P, Savadogo G, Saronga H, Williams J, Sauerborn R, et al. Cost of maternal health services in selected primary care centres in Ghana: a step down allocation approach. BMC Health Serv Res. 2013;13:287.

    Article  Google Scholar 

  27. 27.

    Kasasa S, Asoala V, Gosoniu L, Anto F, Adjuik M, Tindana C, et al. Spatio-temporal malaria transmission patterns in Navrongo demographic surveillance site, northern Ghana. Malar J. 2013;12:63.

    Article  Google Scholar 

  28. 28.

    Manske M, Miotto O, Campino S, Auburn S, Almagro-Garcia J, Maslen G, et al. Analysis of Plasmodium falciparum diversity in natural infections by deep sequencing. Nature. 2012;487:375–9.

    CAS  Article  Google Scholar 

  29. 29.

    Shen Z, Qu W, Wang W, Lu Y, Wu Y, Li Z, et al. MPprimer: a program for reliable multiplex PCR primer design. BMC Bioinform. 2010;11:143.

    Article  Google Scholar 

  30. 30.

    Chang H-H, Worby CJ, Yeka A, Nankabirwa J, Kamya MR, Staedke SG, et al. The real McCoil: a method for the concurrent estimation of the complexity of infection and SNP allele frequency for malaria parasites. PLOS Comput Biol. 2017;13:e1005348.

    Article  Google Scholar 

  31. 31.

    Auburn S, Campino S, Miotto O, Djimde AA, Zongo I, Manske M, et al. Characterization of within-host Plasmodium falciparum diversity using next-generation sequence data. Snounou G, editor. PLoS One. 2012;7:e32891.

  32. 32.

    Appawu M, Owusu-Agyei S, Dadzie S, Asoala V, Anto F, Koram K, et al. Malaria transmission dynamics at a site in northern Ghana proposed for testing malaria vaccines. Trop Med Int Health. 2004;9:164–70.

    Article  Google Scholar 

  33. 33.

    Gesase S, Gosling RD, Hashim R, Ord R, Naidoo I, Madebe 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.

    Article  Google Scholar 

  34. 34.

    O’Meara WP, Smith DL, McKenzie FE. Potential impact of intermittent preventive treatment (IPT) on spread of drug-resistant malaria. PLoS Med. 2006;3:e141.

    Article  Google Scholar 

  35. 35.

    Harrington WE, Mutabingwa TK, Muehlenbachs A, Sorensen B, Bolla MC, Fried M, et al. Competitive facilitation of drug-resistant Plasmodium falciparum malaria parasites in pregnant women who receive preventive treatment. Proc Natl Acad Sci USA. 2009;106:9027–32.

    CAS  Article  Google Scholar 

  36. 36.

    Heinberg A, Kirkman L. The molecular basis of antifolate resistance in Plasmodium falciparum: looking beyond point mutations. Ann N Y Acad Sci. 2015;1342:10–8.

    CAS  Article  Google Scholar 

  37. 37.

    Abuaku B, Duah-Quashie NO, Quaye L, Matrevi SA, Quashie N, Gyasi A, et al. Therapeutic efficacy of artesunate–amodiaquine and artemether–lumefantrine combinations for uncomplicated malaria in 10 sentinel sites across Ghana: 2015–2017. Malar J. 2019;18:206.

    Article  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the Chiefs and people of the Kassena-Nankana Districts in the Upper East Region of northern Ghana, and the staff of the Navrongo Health Research Centre. We are extremely grateful to our study participants and their caregivers. Our gratitude goes to Mr Sebastian Nsoh who assessed the malaria blood smears, Mr Stephen Aberinga for supervising the field work over the years. The authors will also like to thank the staff of Wellcome Trust Sanger Institute Sample Management, Genotyping, Sequencing and Informatics teams for their support.

Genome sequencing was carried out at the Wellcome Sanger Institute and we thank the staff of the Wellcome Sanger Institute Sample Logistics, Sequencing, and Informatics facilities for their contribution.

Funding

This work was supported by MalariaGEN and the Parasites and Microbes program at the Wellcome Sanger Institute with funding from the Wellcome Trust (206194; 090770/Z/09/Z) and by the MRC Centre for Genomics and Global Health which is jointly funded by the Medical Research Council and the Department for International Development (DFID) (G0600718; M006212). The Resource Centre for Genomic Epidemiology of Malaria is supported by Wellcome (090770; 204911). Wellcome also provides core awards to the Wellcome Centre for Human Genetics (203141) and the Wellcome Sanger Institute (206194). Wellcome Trust (DELGEME grant 107740/z/15/z, Aspiring leader award).

Author information

Affiliations

Authors

Contributions

LA, KR and DK conceptualised and design the study. LA, VA, GA supervised the field data collection. LA analysed the data and drafted the manuscript, GA, SG, CJ, KR and DK edited the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Lucas N. Amenga-Etego.

Ethics declarations

Ethics approval and consent to participate

The study was conducted in accordance with the principles and Declaration of Helsinki (1964) as amended in 2015. Ethical approval was obtained from the Institutional Review Board of the Navrongo Health Research Centre (NHRCIRB203). Informed consent was obtained from all adult participants and assent/parental consent for children and minors obtained through their parents or Caregivers.

Consent for publication

All authors have consented for publication.

Competing interests

The authors declare that they have no competing interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Additional file 2: Figure S2.

Mean monthly P. falciparum parasitaemia from 2009 to 2018.

Additional file 3: Figure S1.

Distribution of parasite density by intervention period. 1; is 2009–2011 (pre-IRS) study, 2; is 2013–2015 study (IRS) and 3; is 2016–2018 study (SMC). High coverage of long-lasting insecticidal bed nets from 2013.

Additional file 4: Figure S3

Overall seasonal distribution of Plasmodium falciparum complexity of infections. Panel A: Genome-wide Fws metric from sequenced data (2009–2013). Panel B: Complexity of infection scored using COIL for amplicon data (2014–2018) studies. Study period 1-high season and 2-low season.

Additional file 5: Figure S4

Annual Seasonal distribution of Plasmodium falciparum complexity of infections from 2009–2018. Panel A: Genome-wide Fws metric from sequenced data (2009–2013). 1-low season and 2-high season. Panel B: Complexity of infection scored using COIL for amplicon data (2014–2018) studies.

Additional file 6: Figure S5

Seasonal trends of pfdhfr haplotypes from 2009 to 2018. Panel A: Low; dry season (low malaria transmission season); Panel B: High; Wet season (High malaria transmission season).

Additional file 7: Figure S6

Seasonal trends of pfdhps haplotypes from 2009 to 2018. Panel A: Low; dry season (low malaria transmission season); Panel B: High; Wet season (High malaria transmission season). Hap1-(A/F/Y/S)AKAS; Hap2 -(C/S/A)GKAA; Hap3 – AAKAA; Hap4 -(A/S)GKAS; Hap5-(S/A)GEAA; Hap6 -SAKAA

Additional file 8: Figure S7

Distribution of pfdhfr/pfdhps combined genotypes during high (wet) and low (dry) seasons and pre-and post-SMC periods.

Additional file 9: Figure S8

Temporal trends of pfdhfr/pfdhps combined genotypes during high (wet) and low (dry) seasons and pre-and post-SMC periods.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Amenga-Etego, L.N., Asoala, V., Agongo, G. et al. Temporal evolution of sulfadoxine-pyrimethamine resistance genotypes and genetic diversity in response to a decade of increased interventions against Plasmodium falciparum in northern Ghana. Malar J 20, 152 (2021). https://doi.org/10.1186/s12936-021-03693-3

Download citation

Keywords

  • Malaria
  • SP resistance
  • Genotypes
  • Haplotype
  • Pfdhfr
  • Pfdhps
  • Temporal