Evolution of Plasmodium falciparum drug resistance genes following artemisinin combination therapy in Sudan.

BACKGROUND: Malaria control efforts in Sudan rely heavily on case management. In 2004, health authorities adopted artemisinin-based combination therapies (ACTs) for the treatment of uncomplicated malaria. However, some recent surveys have reported ACT failure and a prevalent irrational malaria treatment practice. Here we examine whether the widespread use of ACT and failure to adhere to national guidelines have led to the evolution of drug resistance genes. METHODS: We genotyped known drug resistance markers (Pfcrt, Pfmdr-1, Pfdhfr, Pfdhps, Pfk13 propeller) and their flanking microsatellites among Plasmodium falciparum isolates obtained between 2009 and 2016 in different geographical regions in Sudan. Data were then compared with published findings pre-ACT (1992-2003). RESULTS: A high prevalence of Pfcrt76T, Pfmdr-1-86Y, Pfdhfr51I, Pfdhfr108N, Pfdhps37G was observed in all regions, while no Pfk13 mutations were detected. Compared with pre-ACT data, Pfcrt-76T and Pfmdr-1-86Y have decayed, while Pfdhfr-51I, Pfdhfr-108N and Pfdhps-437G strengthened. Haplotypes Pfcrt-CVIET, Pfmdr-1-NFSND/YFSND, Pfdhfr-ICNI and Pfdhps-SGKAA predominated in all sites. Microsatellites flanking drug resistance genes showed lower diversity than neutral ones, signifying high ACT pressure/selection. CONCLUSIONS: Evaluation of P. falciparum drug resistance genes in Sudan matches the drug deployment pattern. Regular monitoring of these genes, coupled with clinical response, should be considered to combat the spread of ACT resistance.

Associations between Season and Gametocyte Dynamics in Chronic Plasmodium falciparum Infections.

INTRODUCTION: In a markedly seasonal malaria setting, the transition from the transmission-free dry season to the transmission season depends on the resurgence of the mosquito population following the start of annual rains. The sudden onset of malaria outbreaks at the start of the transmission season suggests that parasites persist during the dry season and respond to either the reappearance of vectors, or correlated events, by increasing the production of transmission stages. Here, we investigate whether Plasmodium falciparum gametocyte density and the correlation between gametocyte density and parasite density show seasonal variation in chronic (largely asymptomatic) carriers in eastern Sudan. MATERIALS AND METHODS: We recruited and treated 123 malaria patients in the transmission season 2001. We then followed them monthly during four distinct consecutive epidemiological seasons: transmission season 1, transmission-free season, pre-clinical period, and transmission season 2. In samples collected from 25 participants who fulfilled the selection criteria of the current analysis, we used quantitative PCR (qPCR) and RT-qPCR to quantify parasite and gametocyte densities, respectively. RESULTS AND DISCUSSION: We observed a significant increase in gametocyte density and a significantly steeper positive correlation between gametocyte density and total parasite density during the pre-clinical period compared to the preceding transmission-free season. However, there was no corresponding increase in the density or prevalence of total parasites or gametocyte prevalence. The increase in gametocyte production during the pre-clinical period supports the hypothesis that P. falciparum may respond to environmental cues, such as mosquito biting, to modulate its transmission strategy. Thus, seasonal changes may be important to ignite transmission in unstable-malaria settings.

Comparative sequence analysis of domain I of Plasmodium falciparum apical membrane antigen 1 from Saudi Arabia and worldwide isolates.

The apical membrane antigen 1 of Plasmodium falciparum (PfAMA1) plays a crucial role in erythrocyte invasion and is a target of protective antibodies. Although domain I of PfAMA1 has been considered a promising vaccine component, extensive sequence diversity in this domain could compromise an effective vaccine design. To explore the extent of sequence diversity in domain I of PfAMA1, P. falciparum-infected blood samples from Saudi Arabia collected between 2007 and 2009 were analyzed and compared with those from worldwide parasite populations. Forty-six haplotypes and a novel codon change (M190V) were found among Saudi Arabian isolates. The haplotype diversity (0.948±0.004) and nucleotide diversity (0.0191±0.0008) were comparable to those from African hyperendemic countries. Positive selection in domain I of PfAMA1 among Saudi Arabian parasite population was observed because nonsynonymous nucleotide substitutions per nonsynonymous site (dN) significantly exceeded synonymous nucleotide substitutions per synonymous site (dS) and Tajima's D and its related statistics significantly deviated from neutrality in the positive direction. Despite a relatively low prevalence of malaria in Saudi Arabia, a minimum of 17 recombination events occurred in domain I. Genetic differentiation was significant between P. falciparum in Saudi Arabia and parasites from other geographic origins. Several shared or closely related haplotypes were found among parasites from different geographic areas, suggesting that vaccine derived from multiple shared epitopes could be effective across endemic countries.

Plasmodium falciparum population structure in Sudan post artemisinin-based combination therapy.

Over the past decade, Sudan has stepped up malaria control backed by WHO, and this has resulted in significant reduction in parasite rate, malaria morbidity and mortality. The present study analyzed Plasmodium falciparum parasites in four geographical separated areas, to examine whether the success in malaria control following the use of artemisinin-based combination therapy (ACT) has disrupted the population structure and evolution of the parasite. We examined 319 P. falciparum isolates collected between October 2009 and October 2012 in four different areas in Sudan (Jazira [central Sudan], Southern Darfur [western Sudan], Upper Nile [southern Sudan] and Kasala [eastern Sudan]). Twelve microsatellites were analyzed for allelic diversity, multi-locus haplotype and inter-population differentiation. Level of diversity was compared to that detected for three of the above microsatellites among P. falciparum parasites in central and eastern Sudan in 1999, prior to introduction of ACT. Diversity at each locus (unbiased heterozygosity [H]) was high in all areas (Jazira, H=0.67), (Southern Darfur, H=0.71), (Upper Nile, H=0.71), and (Kasala, H=0.63). Microsatellites were distributed widely and private alleles, detected in a single population, were rare. The extent of diversity in the above sites was similar to that seen, in 1999, in central (Khartoum, H=0.73) and eastern Sudan (Gedaref, H=0.75). Significant Linkage disequilibrium (LD) was observed between the microsatellites in all populations. Pairwise FST analysis revealed that parasites in the four areas could be considered as one population. However, the parasites in Sudan clustered away from parasites in West Africa and the Arabian Peninsula. Despite marked reduction in malaria risk in Sudan, the extent of diversity and parasite genetic structure are indicative of a large population size. Further considerable reduction in transmission would be needed before fragmented sub-population can be seen. In addition, the large divergence of P. falciparum in Sudan from West Africa and Arabian Peninsula populations may result from differential evolutionary pressures acting at the population level, which shall be considered in eradication plans.

Distinct haplotypes of dhfr and dhps among Plasmodium falciparum isolates in an area of high level of sulfadoxine-pyrimethamine (SP) resistance in eastern Sudan.

Typing of polymorphic microsatellites that are linked to drug resistance genes has shed light on the origin and pattern of spread of some anti-malarial drugs. Recent surveys revealed spread of a high-level pyrimethemine resistant lineage of Plasmodium falciparum, of Asian origin, across Africa. Here, we examined mutations in dihydrofolate reductase, dhfr [chromsosome 4], the dihydropteroate synthase, dhps [chromosome 8] associated with resistance to sulfadoxine-pyrimethamine (SP), and neighboring microsatellites among P. falciparum isolates in Asar village, eastern Sudan. This area lies at the fringes of malaria endemicity, where the remote P. falciparum parasites have some distinct genetic characteristics. Overall, 89% (84/94) of the examined isolates carried double mutations at dhfr (N51I and S108N), but the 59R and I164L mutations were not seen. Similarly, the majority, 43% (35/81) of the isolates carried double mutations at dhps (437G, 540E). Analysis of neighboring microsatellites revealed one major dhfr haplotype with mutations (51I, 108N) and one dhps haplotype with mutations (436S, 437G, 540E). These haplotypes differ from the major ones thought to drive resistance to SP across Africa. The resistant haplotypes of dhfr and dhps, in Asar, share some microsatellites with the wild genotypes suggesting that they were generated locally. Among isolates successfully examined, 40% shared identical haplotypes of the 2 loci, comprising a dominant resistant lineage. Undoubtedly, this lineage plays an important role in clinical failure to SP in this area.

Allelic dimorphism-associated restriction of recombination in Plasmodium falciparum msp1.

Allelic dimorphism is a characteristic feature of the Plasmodium falciparum msp1 gene encoding the merozoite surface protein 1, a strong malaria vaccine candidate. Meiotic recombination is a major mechanism for the generation of msp1 allelic diversity. Potential recombination sites have previously been mapped to specific regions within msp1 (a 5' 1-kb region and a 3' 0.4-kb region) with no evidence for recombination events in a central 3.5-kb region. However, evidence for the lack of recombination events is circumstantial and inconclusive because the number of msp1 sequences analysed is limited, and the frequency of recombination events has not been addressed previously in a high transmission area, where the frequency of meiotic recombination is expected to be high. In the present study, we have mapped potential allelic recombination sites in 34 full-length msp1 sequences, including 24 new sequences, from various geographic origins. We also investigated recombination events in blocks 6 to 16 by population genetic analysis of P. falciparum populations in Tanzania, where malaria transmission is intense. The results clearly provide no evidence of recombination events occurring between the two major msp1 allelic types, K1-type and Mad20-type, in the central region, but do show recombination events occurring throughout the entire gene within sequences of the Mad20-type. Thus, the present study indicates that allelic dimorphism of msp1 greatly affects inter-allelic recombination events, highlighting a unique feature of allelic diversity of P. falciparum msp1.

Increased density but not prevalence of gametocytes following drug treatment of Plasmodium falciparum.

We monitored post-treatment Plasmodium falciparum among patients treated with chloroquine (CQ) and sulfadoxine-pyrimethamine (SP; Fansidar in a village in eastern Sudan. Parasites were examined on day 0 (pre-treatment), day 7, day 14 and day 21 (post-treatment) during the transmission season. A further sample was taken 2 months later (day 80) at the start of the dry season. Asexual forms and gametocytes were detected by microscopy, and reverse transcriptase polymerase chain reaction (RT-PCR) was used to detect expression of gametocyte-specific proteins pfs 25 and pfg 377. Gametocyte carriage, as revealed by microscopy, increased significantly following CQ and SP treatment, reaching a maximum between days 7 and 14. When measured by RT-PCR, however, there was no significant difference in gametocyte rate between day 0 and days 7 or 14. RT-PCR gametocyte rates dropped dramatically by day 80 post treatment but were still 33% and 8% in the CQ- and SP-treated group at this time. Alleles associated with drug resistance of P. falciparum to chloroquine (the chloroquine resistance transporter, pfcrt, and multidrug resistance, pfmdr1) and to pyrimethamine (dihydrofolate reductase, dhfr) were seen at a high frequency at the beginning of treatment and increased further through time following both drug treatments. Infections with drug-resistant parasites tended to have higher gametocyte prevalence than drug-sensitive infections.

Impact of genetic complexity on longevity and gametocytogenesis of Plasmodium falciparum during the dry and transmission-free season of eastern Sudan.

Malaria in eastern Sudan is characterised by limited seasonal transmission, with the majority of the year remaining transmission-free. Some inhabitants who contract malaria during the transmission season retain long-lasting sub-patent infections, which probably initiate transmission the following year. Here we have monitored Plasmodium falciparum infection prevalence and gametocyte production during the dry season, and examined the impact of parasite genetic multiplicity on infection longevity. A cohort of 38 individuals who were infected with P. falciparum in November 2001 was monitored monthly by microscopy and PCR until December 2002. Reverse transcriptase polymerase chain reaction of the pfg377 gene was used to detect sub-patent gametocytes. In addition, all isolates were examined for msp-2 alleles and the mean number of parasite clones per infection was estimated. We found that a large proportion (40%) of the cohort retained gametocytes throughout the dry season. The majority of patients retained asexual infection for at least 7 months. Genetic multiplicity of P. falciparum significantly influenced longevity of asexual infection and its gametocyte production. Gametocytes from mixed genotype P. falciparum infections persisted three times longer than those from single genotype infections, suggesting that genetic diversity promotes persistence. These findings are discussed in the context of the parasite biology and malaria epidemiology in the study area.

Evolution of drug-resistance genes in Plasmodium falciparum in an area of seasonal malaria transmission in Eastern Sudan.

We investigated the evolution of drug-resistant Plasmodium falciparum in a village in eastern Sudan. The frequencies of alleles of 4 genes thought to be determinants of drug resistance were monitored from 1990 through 2001. Changes in frequencies of drug-resistance genes between wet and dry seasons were monitored from 1998 through 2000. Parasites were also typed for 3 putatively neutral microsatellite loci. No significant variation in frequencies was observed for the microsatellite loci over the whole study period or between seasons. However, genes involved in resistance to chloroquine showed consistent, significant increases in frequencies over time (rate of annual increase, 0.027/year for pfcrt and 0.018/year for pfmdr1). Genes involved in resistance to the second-line drug used in the area (Fansidar) remained at low frequencies between 1990 and 1993 but increased dramatically between 1998 and 2000, which is consistent with the advent of Fansidar usage during this period. For mutant alleles of the primary drug-resistance targets for chloroquine and pyrimethamine, higher frequencies were seen during the dry season than during the wet season. This cyclical fluctuation in drug-resistance genes most likely reflects seasonal variation in drug pressure and differences in the fitness of resistant and sensitive parasites.

Detection of mutations in the Plasmodium falciparum dihydrofolate reductase (dhfr) gene by dot-blot hybridization.

There is a need for a specific, sensitive, robust, and large-scale method for diagnosis of drug resistance genes in natural Plasmodium falciparum infections. Established polymerase chain reaction (PCR)-based methods may be compromised by the multiplicity of P. falciparum genotypes in natural infections. Here we adopt a dot-blot method to detect point mutations at nucleotide 323 (residue 108) in the P. falciparum dihydrofolate reductase (dhfr) gene using allele-specific oligonucleotide probes. Serine (Ser) or threonine (Thr) at this position are associated with sensitivity to pyrimethamine while asparagine (Asn) is associated with resistance. The method combines PCR amplification and hybridization of amplified products with radiolabeled allele-specific probes. This technique is specific and sensitive; it detects parasitemia of less than 100 parasites/microl of blood, and can identify a minority parasite genotype down to 1% in a mixture. Analysis of P. falciparum isolates from Sudan, of known response to pyrimethamine, has demonstrated the sensitivity and specificity of the method and its ability to detect multiple genotypes in single infections. Furthermore, it has confirmed the association between pyrimethamine responses and dhfr alleles. The method has been successfully extended for analysis of other point mutations in dhfr at residues 51 and 59, which are associated with a high level of pyrimethamine resistance.

Dynamics of gametocytes among Plasmodium falciparum clones in natural infections in an area of highly seasonal transmission.

The dynamics of gametocyte production in Plasmodium falciparum clones were studied in inhabitants of an area of highly seasonal malaria transmission in eastern Sudan. Reverse-transcriptase polymerase chain reaction was used to detect expression of 2 genes that encode gametocyte-specific proteins, pfs25 and pfg377, in parasites sampled from individuals throughout one year. Some patients who acquired infections during the wet season were found to harbor subpatent gametocytemia through the following dry season in the apparent absence of mosquito transmission. Genotyping of parasites in multiclonal infections showed considerable fluctuation of gametocyte production by individual clones. The gametocytes present at the end of the dry season provide the most probable source of the genetically complex cyclical malaria outbreaks following the rainy season in this region.

A marked seasonality of malaria transmission in two rural sites in eastern Sudan.

The ecology of Anopheles arabiensis and its relationship to malaria transmission was investigated in two villages in eastern Sudan. Seasonal malaria case incidence was compared with the number of vectors detected and with climatic variables. Following the end of the short rainy season in October the number of A. arabiensis detected dropped gradually until February when neither outdoor human bait trapping nor indoor spray catches revealed any mosquitoes. Vectors re-appeared in June as humidity rose with the onset of rain. Despite the apparent absence of the vector at the height of the long, hot dry season between February and May, sporadic asymptomatic malaria infections were detected in the two villages. The low endemicity of malaria in the area was reflected by the relatively low total September-December parasite and sporozoite rates (15 and 1.4%, respectively) measured in the villages. The entomological inoculation rate (EIR) was estimated to be around two to three infective bites per person per year, although heterogeneity in the transmission indices of malaria between the two villages was observed. The implications of these patterns of anopheline population dynamics for the epidemiology and control of malaria in eastern Sudan are considered.