Historical shifts in Brazilian P. falciparum population structure and drug resistance alleles

Malaria Molecular Surveillance in the Peruvian Amazon with a Novel Highly Multiplexed Plasmodium falciparum AmpliSeq Assay

Molecular surveillance for malaria has great potential to support national malaria control programs (NMCPs). To bridge the gap between research and implementation, several applications (use cases) have been identified to align research, technology development, and public health efforts. For implementation at NMCPs, there is an urgent need for feasible and cost-effective tools. We designed a new highly multiplexed deep sequencing assay (Pf AmpliSeq), which is compatible with benchtop sequencers, that allows high-accuracy sequencing with higher coverage and lower cost than whole-genome sequencing (WGS), targeting genomic regions of interest. The novelty of the assay is its high number of targets multiplexed into one easy workflow, combining population genetic markers with 13 nearly full-length resistance genes, which is applicable for many different use cases. We provide the first proof of principle for hrp2 and hrp3 deletion detection using amplicon sequencing. Initial sequence data processing can be performed automatically, and subsequent variant analysis requires minimal bioinformatic skills using any tabulated data analysis program. The assay was validated using a retrospective sample collection (n = 254) from the Peruvian Amazon between 2003 and 2018. By combining phenotypic markers and a within-country 28-single-nucleotide-polymorphism (SNP) barcode, we were able to distinguish different lineages with multiple resistance haplotypes (in dhfr, dhps, crt and mdr1) and hrp2 and hrp3 deletions, which have been increasing in recent years. We found no evidence to suggest the emergence of artemisinin (ART) resistance in Peru. These findings indicate a parasite population that is under drug pressure but is susceptible to current antimalarials and demonstrate the added value of a highly multiplexed molecular tool to inform malaria strategies and surveillance systems. IMPORTANCE While the power of next-generation sequencing technologies to inform and guide malaria control programs has become broadly recognized, the integration of genomic data for operational incorporation into malaria surveillance remains a challenge in most countries where malaria is endemic. The main obstacles include limited infrastructure, limited access to high-throughput sequencing facilities, and the need for local capacity to run an in-country analysis of genomes at a large-enough scale to be informative for surveillance. In addition, there is a lack of standardized laboratory protocols and automated analysis pipelines to generate reproducible and timely results useful for relevant stakeholders. With our standardized laboratory and bioinformatic workflow, malaria genetic surveillance data can be readily generated by surveillance researchers and malaria control programs in countries of endemicity, increasing ownership and ensuring timely results for informed decision- and policy-making.

Genetic diversity and genotype multiplicity of Plasmodium falciparum infection in patients with uncomplicated malaria in Chewaka district, Ethiopia

Background

Genetic diversity in Plasmodium falciparum poses a major threat to malaria control and elimination interventions. Characterization of the genetic diversity of P. falciparum strains can be used to assess intensity of parasite transmission and identify potential deficiencies in malaria control programmes, which provides vital information to evaluating malaria elimination efforts. This study investigated the P. falciparum genetic diversity and genotype multiplicity of infection in parasite isolates from cases with uncomplicated P. falciparum malaria in Southwest Ethiopia.

Methods

A total of 80 P. falciparum microscopy and qPCR positive blood samples were collected from study participants aged 6 months to 60 years, who visited the health facilities during study evaluating the efficacy of artemether-lumefantrine from September–December, 2017. Polymorphic regions of the msp-1 and msp-2 were genotyped by nested polymerase chain reactions (nPCR) followed by gel electrophoresis for fragment analysis.

Results

Of 80 qPCR-positive samples analysed for polymorphisms on msp-1 and msp-2 genes, the efficiency of msp-1 and msp-2 gene amplification reactions with family-specific primers were 95% and 98.8%, respectively. Allelic variation of 90% (72/80) for msp-1 and 86.2% (69/80) for msp-2 were observed. K1 was the predominant msp-1 allelic family detected in 20.8% (15/72) of the samples followed by MAD20 and RO33. Within msp-2, allelic family FC27 showed a higher frequency (26.1%) compared to IC/3D7 (15.9%). Ten different alleles were observed in msp-1 with 6 alleles for K1, 3 alleles for MAD20 and 1 allele for RO33. In msp-2, 19 individual alleles were detected with 10 alleles for FC27 and 9 alleles for 3D7. Eighty percent (80%) of isolates had multiple genotypes and the overall mean multiplicity of infection was 3.2 (95% CI 2.87–3.46). The heterozygosity indices were 0.43 and 0.85 for msp-1 and msp-2, respectively. There was no significant association between multiplicity of infection and age or parasite density.

Conclusions

The study revealed high levels of genetic diversity and mixed-strain infections of P. falciparum populations in Chewaka district, Ethiopia, suggesting that both endemicity level and malaria transmission remain high and that strengthened control efforts are needed in Ethiopia.

Novel Quinolinyl-pyrrolo[3,4-d]pyrimidine-2,5-dione Derivatives Against Chloroquine-resistant Plasmodium falciparum

INTRODUCTION

In this work DHPMs were combined with the quinoline nucleus to obtain new quinolinyl-pyrrolo[3,4-d]pyrimidine-2,5-dione compounds with improved antiplasmodial activity as well as decreased cytotoxicity. Nineteen quinolinyl-pyrrolo[3,4-d]pyrimidine-2,5-dione derivatives connected by a linker group to quinolone ring moieties with different substituents were synthesized and assayed against P. falciparum.

MATERIALS AND METHODS

Nineteen quinolinyl-pyrrolo[3,4-d]pyrimidine-2,5-dione derivatives connected by a linker group to quinoline ring moieties with different substituents were synthesized and assayed against chloroquine-resistant Plasmodium falciparum, along with the reference drug chloroquine. Among these compounds, the derivatives with two methylene carbon spacers showed the best activity accompanied by low cytotoxicity.

RESULTS

The derivative without substituents on the aromatic ring (2a) and the derivative with a chlorine group at position 4 (2d) provided the best results, with IC50 = 1.15 µM and 1.5 µM, respectively.

CONCLUSION

Compared to the parent drugs, these compounds presented marked decreases in cytotoxicity, with MDL50 values over 1,000 µM and selectivity indexes of >869.5 and >666.6, respectively. The quinolinyl-pyrrolo[3,4-d]pyrimidine-2,5-dione framework appears to be promising for further studies as an antimalarial for overcoming the burden of resistance in P. falciparum.

Diverse origin of Plasmodium falciparum in northwest Ecuador

BACKGROUND

Ecuador plans to eliminate malaria by 2020, and the country has already seen a decrease in the number of cases from more than 100,000 in 2000 to only 618 in 2015. Around 30% of malaria infections in Ecuador are caused by Plasmodium falciparum. Most malaria population genetics studies performed in Latin America, especially in the Pacific Coast, indicate a high clonality and a clear structure of P. falciparum populations. It was shown that an outbreak of P. falciparum in northwest Ecuador was the result of a clonal expansion of parasites circulating at low levels in the country or re-invading Ecuador from neighbouring territories. However, general characteristics of P. falciparum circulating in the northwest coast of Ecuador have not been determined. The main goal of this study was to genetically characterize the population structure of P. falciparum in coastal Ecuadorian localities bordering with Colombia.

METHODS

Molecular investigation of 41 samples collected from 2013 to 2016 in San Lorenzo County, northwest Ecuador was performed using seven neutral microsatellite markers.

RESULTS

The genetic population structure of P. falciparum in northwest Ecuador is clearly defined as three different genetic groups previously reported in Ecuador, Peru and Colombia.

CONCLUSIONS

The limited number of P. falciparum clonal types that are circulating in northwest Ecuador, are related to ancestral parasite clonal lineages reported in the Pacific Coast. These parasites could be a product of migration from neighbouring regions or residual clonal types circulating in the country in low proportions. Studies of the genetic characterization of P. falciparum in eliminating areas help determine the possible origin of parasites in order to create strategies to prevent the entrance of new lineages and achieve local elimination of malaria.

Genome-wide diversity and differentiation in New World populations of the human malaria parasite Plasmodium vivax

Background

The Americas were the last continent colonized by humans carrying malaria parasites. Plasmodium falciparum from the New World shows very little genetic diversity and greater linkage disequilibrium, compared with its African counterparts, and is clearly subdivided into local, highly divergent populations. However, limited available data have revealed extensive genetic diversity in American populations of another major human malaria parasite, P. vivax.

Methods

We used an improved sample preparation strategy and next-generation sequencing to characterize 9 high-quality P. vivax genome sequences from northwestern Brazil. These new data were compared with publicly available sequences from recently sampled clinical P. vivax isolates from Brazil (BRA, total n = 11 sequences), Peru (PER, n = 23), Colombia (COL, n = 31), and Mexico (MEX, n = 19).

Principal findings/Conclusions

We found that New World populations of P. vivax are as diverse (nucleotide diversity π between 5.2 × 10⁻⁴ and 6.2 × 10⁻⁴) as P. vivax populations from Southeast Asia, where malaria transmission is substantially more intense. They display several non-synonymous nucleotide substitutions (some of them previously undescribed) in genes known or suspected to be involved in antimalarial drug resistance, such as dhfr, dhps, mdr1, mrp1, and mrp-2, but not in the chloroquine resistance transporter ortholog (crt-o) gene. Moreover, P. vivax in the Americas is much less geographically substructured than local P. falciparum populations, with relatively little between-population genome-wide differentiation (pairwise F_ST values ranging between 0.025 and 0.092). Finally, P. vivax populations show a rapid decline in linkage disequilibrium with increasing distance between pairs of polymorphic sites, consistent with very frequent outcrossing. We hypothesize that the high diversity of present-day P. vivax lineages in the Americas originated from successive migratory waves and subsequent admixture between parasite lineages from geographically diverse sites. Further genome-wide analyses are required to test the demographic scenario suggested by our data.

Plasmodium vivax is the most common human malaria parasite in the Americas, but how and when this species arrived in the New World remains unclear. Here we describe high-quality whole-genome sequence data for nine P. vivax isolates from Brazil, a country that accounts for 37% of the malaria burden in this continent, and compare these data with additional publicly available P. vivax genomes from Brazil, Peru, Colombia, and Mexico. P. vivax populations from the New World were found to be as diverse as their counterparts from areas with substantially higher malaria transmission, such as Southeast Asia, and to carry several non-synonymous substitutions in candidate drug-resistance genes. Moreover, genome-wide patterns of linkage disequilibrium between pairs of polymorphic sites are consistent with very frequent outcrossing in these populations. Interestingly, local P. vivax is more polymorphic, with less between-population differentiation, than sympatric populations of P. falciparum, possibly as a result of different demographic histories of these two species in the Americas. We hypothesize that local P. vivax lineages originated from successive migratory waves and subsequent admixture between parasites from geographically diverse sites.

Molecular Profile of Malaria Drug Resistance Markers of Plasmodium falciparum in Suriname

In Suriname, an artesunate monotherapy therapeutic efficacy trial was recently conducted to evaluate partial artemisinin resistance emerging in Plasmodium falciparum We genotyped the PfK13 propeller domain of P. falciparum in 40 samples as well as other mutations proposed to be associated with artemisinin-resistant mutants. We did not find any mutations previously associated with artemisinin resistance in Southeast Asia, but we found fixed resistance mutations for chloroquine (CQ) and sulfadoxine-pyrimethamine. Additionally, the PfCRT C350R mutation, associated with reversal of CQ resistance and piperaquine-selective pressure, was present in 62% of the samples. Our results from neutral microsatellite data also confirmed a high parasite gene flow in the Guiana Shield. Although recruiting participants for therapeutic efficacy studies is challenging in areas where malaria endemicity is very low due to the low number of malaria cases reported, conducting these studies along with molecular surveillance remains essential for the monitoring of artemisinin-resistant alleles and for the characterization of the population structure of P. falciparum in areas targeted for malaria elimination.

Plasmodium falciparum Drug-Resistant Haplotypes and Population Structure in Postearthquake Haiti, 2010

Chloroquine (CQ) remains the first-line treatment of malaria in Haiti. Given the challenges of conducting in vivo drug efficacy trials in low-endemic settings like Haiti, molecular surveillance for drug resistance markers is a reasonable approach for detecting resistant parasites. In this study, 349 blood spots were collected from suspected malaria cases in areas in and around Port-au-Prince from March to July 2010. Among them, 121 samples that were Plasmodium falciparum positive by polymerase chain reaction were genotyped for drug-resistant pfcrt, pfdhfr, pfdhps, and pfmdr1 alleles. Among the 108 samples that were successfully sequenced for CQ resistant markers in pfcrt, 107 were wild type (CVMNK), whereas one sample carried a CQ-resistant allele (CVIET). Neutral microsatellite genotyping revealed that the CQ-resistant isolate was distinct from all other samples in this study. Furthermore, the remaining parasite specimens appeared to be genetically distinct from other reported Central and South American populations.

Genetic Diversity of Plasmodium falciparum Populations in Malaria Declining Areas of Sabah, East Malaysia

Malaysia has a national goal to eliminate malaria by 2020. Understanding the genetic diversity of malaria parasites in residual transmission foci can provide invaluable information which may inform the intervention strategies used to reach elimination targets. This study was conducted to determine the genetic diversity level of P. falciparum isolates in malaria residual foci areas of Sabah. Malaria active case detection was conducted in Kalabakan and Kota Marudu. All individuals in the study sites were screened for malaria infection by rapid diagnostic test. Blood from P. falciparum-infected individuals were collected on filter paper prior to DNA extraction. Genotyping was performed using merozoite surface protein-1 (MSP-1), merozoite surface protein-2 (MSP-2), glutamate rich protein (GLURP) and 10 neutral microsatellite loci markers. The size of alleles, multiplicity of infection (MOI), mean number of alleles (Na), expected heterozygosity (He), linkage disequilibrium (LD) and genetic differentiation (FST) were determined. In Kalabakan, the MSP-1 and MSP-2 alleles were predominantly K1 and FC27 family types, respectively. The GLURP genotype VI (751-800 bp) was predominant. The MOI for MSP-1 and MSP-2 were 1.65 and 1.20, respectively. The Na per microsatellite locus was 1.70. The He values for MSP-1, MSP-2, GLURP and neutral microsatellites were 0.17, 0.37, 0.70 and 0.33, respectively. In Kota Marudu, the MSP-1 and MSP-2 alleles were predominantly MAD20 and 3D7 family types, respectively. The GLURP genotype IV (651-700 bp) was predominant. The MOI for both MSP-1 and MSP-2 was 1.05. The Na per microsatellite locus was 3.60. The He values for MSP-1, MSP-2, GLURP and neutral microsatellites were 0.24, 0.25, 0.69 and 0.30, respectively. A significant LD was observed in Kalabakan (0.495, p<0.01) and Kota Marudu P. falciparum populations (0.601, p<0.01). High genetic differentiation between Kalabakan and Kota Marudu P. falciparum populations was observed (FST = 0.532). The genetic data from the present study highlighted the limited diversity and contrasting genetic pattern of P. falciparum populations in the malaria declining areas of Sabah.

Clonal population expansion in an outbreak of Plasmodium falciparum on the northwest coast of Ecuador

Background

Determining the source of malaria outbreaks in Ecuador and identifying remaining transmission foci will help in malaria elimination efforts. In this study, the genetic signatures of Plasmodium falciparum isolates, obtained from an outbreak that occurred in northwest Ecuador from 2012 to 2013, were characterized.

Methods

Molecular investigation of the outbreak was performed using neutral microsatellites, drug resistance markers and pfhrp2 and pfhrp3 genotyping.

Results

A majority of parasite isolates (31/32) from this outbreak were of a single clonal type that matched a clonal lineage previously described on the northern coast of Peru and a historical isolate from Ecuador. All but one isolate carried a chloroquine-resistant pfcrt genotype and sulfadoxine- and pyrimethamine-sensitive pfdhps and pfdhfr genotypes. Pfmdr1 mutations were identified in codons 184 and 1042. In addition, most samples (97 %) showed presence of pfhrp2 gene.

Conclusions

This study indicates that parasites from a single clonal lineage largely contributed to this outbreak and this lineage was found to be genetically related to a lineage previously reported in the Peruvian coast and historical Ecuadorian parasites.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-1019-2) contains supplementary material, which is available to authorized users.

Molecular Investigation into a Malaria Outbreak in Cusco, Peru: Plasmodium falciparum BV1 Lineage is Linked to a Second Outbreak in Recent Times

In November 2013, a Plasmodium falciparum malaria outbreak of 11 cases occurred in Cusco, southern Peru, where falciparum malaria had not been reported since 1946. Although initial microscopic diagnosis reported only Plasmodium vivax infection in each of the specimens, subsequent examination by the national reference laboratory confirmed P. falciparum infection in all samples. Molecular typing of four available isolates revealed identity as the B-variant (BV1) strain that was responsible for a malaria outbreak in Tumbes, northern Peru, between 2010 and 2012. The P. falciparum BV1 strain is multidrug resistant, can escape detection by PfHRP2-based rapid diagnostic tests, and has contributed to two malaria outbreaks in Peru. This investigation highlights the importance of accurate species diagnosis given the potential for P. falciparum to be reintroduced to regions where it may have been absent. Similar molecular epidemiological investigations can track the probable source(s) of outbreak parasite strains for malaria surveillance and control purposes.

Genetic Diversity of Plasmodium falciparum in Haiti: Insights from Microsatellite Markers

Hispaniola, comprising Haiti and the Dominican Republic, has been identified as a candidate for malaria elimination. However, incomplete surveillance data in Haiti hamper efforts to assess the impact of ongoing malaria control interventions. Characteristics of the genetic diversity of Plasmodium falciparum populations can be used to assess parasite transmission, which is information vital to evaluating malaria elimination efforts. Here we characterize the genetic diversity of P. falciparum samples collected from patients at seven sites in Haiti using 12 microsatellite markers previously employed in population genetic analyses of global P. falciparum populations. We measured multiplicity of infections, level of genetic diversity, degree of population geographic substructure, and linkage disequilibrium (defined as non-random association of alleles from different loci). For low transmission populations like Haiti, we expect to see few multiple infections, low levels of genetic diversity, high degree of population structure, and high linkage disequilibrium. In Haiti, we found low levels of multiple infections (12.9%), moderate to high levels of genetic diversity (mean number of alleles per locus = 4.9, heterozygosity = 0.61), low levels of population structure (highest pairwise F_st = 0.09 and no clustering in principal components analysis), and moderate linkage disequilibrium (ISA = 0.05, P<0.0001). In addition, population bottleneck analysis revealed no evidence for a reduction in the P. falciparum population size in Haiti. We conclude that the high level of genetic diversity and lack of evidence for a population bottleneck may suggest that Haiti’s P. falciparum population has been stable and discuss the implications of our results for understanding the impact of malaria control interventions. We also discuss the relevance of parasite population history and other host and vector factors when assessing transmission intensity from genetic diversity data.

Longitudinal analysis of Plasmodium falciparum genetic variation in Turbo, Colombia: implications for malaria control and elimination

Background

Malaria programmes estimate changes in prevalence to evaluate their efficacy. In this study, parasite genetic data was used to explore how the demography of the parasite population can inform about the processes driving variation in prevalence. In particular, how changes in treatment and population movement have affected malaria prevalence in an area with seasonal malaria.

Methods

Samples of Plasmodium falciparum collected over 8 years from a population in Turbo, Colombia were genotyped at nine microsatellite loci and three drug-resistance loci. These data were analysed using several population genetic methods to detect changes in parasite genetic diversity and population structure. In addition, a coalescent-based method was used to estimate substitution rates at the microsatellite loci.

Results

The estimated mean microsatellite substitution rates varied between 5.35 × 10⁻³ and 3.77 × 10⁻² substitutions/locus/month. Cluster analysis identified six distinct parasite clusters, five of which persisted for the full duration of the study. However, the frequencies of the clusters varied significantly between years, consistent with a small effective population size.

Conclusions

Malaria control programmes can detect re-introductions and changes in transmission using rapidly evolving microsatellite loci. In this population, the steadily decreasing diversity and the relatively constant effective population size suggest that an increase in malaria prevalence from 2004 to 2007 was primarily driven by local rather than imported cases.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0887-9) contains supplementary material, which is available to authorized users.

Serologically defined variations in malaria endemicity in Pará state, Brazil

BACKGROUND

Measurement of malaria endemicity is typically based on vector or parasite measures. A complementary approach is the detection of parasite specific IgG antibodies. We determined the antibody levels and seroconversion rates to both P. vivax and P. falciparum merozoite antigens in individuals living in areas of varying P. vivax endemicity in Pará state, Brazilian Amazon region.

METHODOLOGY/PRINCIPAL FINDINGS

The prevalence of antibodies to recombinant antigens from P. vivax and P. falciparum was determined in 1,330 individuals. Cross sectional surveys were conducted in the north of Brazil in Anajás, Belém, Goianésia do Pará, Jacareacanga, Itaituba, Trairão, all in the Pará state, and Sucuriju, a free-malaria site in the neighboring state Amapá. Seroprevalence to any P. vivax antigens (MSP1 or AMA-1) was 52.5%, whereas 24.7% of the individuals were seropositive to any P. falciparum antigens (MSP1 or AMA-1). For P. vivax antigens, the seroconversion rates (SCR) ranged from 0.005 (Sucuriju) to 0.201 (Goianésia do Pará), and are strongly correlated to the corresponding Annual Parasite Index (API). We detected two sites with distinct characteristics: Goianésia do Pará where seroprevalence curve does not change with age, and Sucuriju where seroprevalence curve is better described by a model with two SCRs compatible with a decrease in force of infection occurred 14 years ago (from 0.069 to 0.005). For P. falciparum antigens, current SCR estimates varied from 0.002 (Belém) to 0.018 (Goianésia do Pará). We also detected a putative decrease in disease transmission occurred ∼29 years ago in Anajás, Goianésia do Pará, Itaituba, Jacareacanga, and Trairão.

CONCLUSIONS

We observed heterogeneity of serological indices across study sites with different endemicity levels and temporal changes in the force of infection in some of the sites. Our study provides further evidence that serology can be used to measure and monitor transmission of both major species of malaria parasite.

Conservation efforts and malaria in the Brazilian Amazon

We respond to Valle and Clark, who assert that "conservation efforts may increase malaria burden in the Brazilian Amazon," because the relationship between forest cover and malaria incidence was stronger than the effect of the deforestation rate. We contend that their conclusion is flawed because of limitations in their methodology that we discuss in detail. Most important are the exclusion of one-half the original data without a discussion of selection bias, the lack of model adjustment for either population growth or migration, and the crude classifications of land cover and protected areas that lead to aggregation bias. Of greater significance, we stress the need for caution in the interpretation of data that could have profound effects on regional land use decisions.

Genetic structure of Plasmodium falciparum populations across the Honduras-Nicaragua border

Background

The Caribbean coast of Central America remains an area of malaria transmission caused by Plasmodium falciparum despite the fact that morbidity has been reduced in recent years. Parasite populations in that region show interesting characteristics such as chloroquine susceptibility and low mortality rates. Genetic structure and diversity of P. falciparum populations in the Honduras-Nicaragua border were analysed in this study.

Methods

Seven neutral microsatellite loci were analysed in 110 P. falciparum isolates from endemic areas of Honduras (n = 77) and Nicaragua (n = 33), mostly from the border region called the Moskitia. Several analyses concerning the genetic diversity, linkage disequilibrium, population structure, molecular variance, and haplotype clustering were conducted.

Results

There was a low level of genetic diversity in P. falciparum populations from Honduras and Nicaragua. Expected heterozigosity (H_e) results were similarly low for both populations. A moderate differentiation was revealed by the F_ST index between both populations, and two putative clusters were defined through a structure analysis. The main cluster grouped most of samples from Honduras and Nicaragua, while the second cluster was smaller and included all the samples from the Siuna community in Nicaragua. This result could partially explain the stronger linkage disequilibrium (LD) in the parasite population from that country. These findings are congruent with the decreasing rates of malaria endemicity in Central America.

Population genetics analysis during the elimination process of Plasmodium falciparum in Djibouti

Background

Case management of imported malaria within the context of malaria pre-elimination is increasingly considered to be relevant because of the risk of resurgence. The assessment of malaria importation would provide key data i) to select countries with propitious conditions for pre-elimination phase and ii) to predict its feasibility. Recently, a sero-prevalence study in Djibouti indicated low malaria prevalence, which is propitious for the implementation of pre-elimination, but data on the extent of malaria importation remain unknown.

Methods

Djiboutian plasmodial populations were analysed over an eleven-year period (1998, 1999, 2002 and 2009). The risk of malaria importation was indirectly assessed by using plasmodial population parameters. Based on 5 microsatellite markers, expected heterozygosity (H.e.), multiplicity of infection, pairwise Fst index, multiple correspondence analysis and individual genetic relationship were determined. The prevalence of single nucleotide polymorphisms associated with pyrimethamine resistance was also determined.

Results

Data indicated a significant decline in genetic diversity (0.51, 0.59, 0.51 and 0 in 1998, 1999, 2002 and 2009, respectively) over the study period, which is inconsistent with the level of malaria importation described in a previous study. This suggested that Djiboutian malaria situation may have benefited from the decline of malaria prevalence that occurred in neighbouring countries, in particular in Ethiopia. The high Fst indices derived from plasmodial populations from one study period to another (0.12 between 1999 and 2002, and 0.43 between 2002 and 2009) suggested a random sampling of parasites, probably imported from neighbouring countries, leading to oligo-clonal expansion of few different strains during each transmission season. Nevertheless, similar genotypes observed during the study period suggested recurrent migrations and imported malaria.

Conclusion

In the present study, the extent of genetic diversity was used to assess the risk of malaria importation in the low malaria transmission setting of Djibouti. The molecular approach highlights i) the evolution of Djiboutian plasmodial population profiles that are consistent and compatible with Djiboutian pre-elimination goals and ii) the necessity to implement the monitoring of plasmodial populations and interventions at the regional scale in the Horn of Africa to ensure higher efficiency of malaria control and elimination.