Skip to main content

Effects of transmission reduction by insecticide-treated bed nets (ITNs) on parasite genetics population structure: I. The genetic diversity of Plasmodium falciparum parasites by microsatellite markers in western Kenya

Abstract

Background

Insecticide-treated bed nets (ITNs) reduce malaria transmission and are an important prevention tool. However, there are still information gaps on how the reduction in malaria transmission by ITNs affects parasite genetics population structure. This study examined the relationship between transmission reduction from ITN use and the population genetic diversity of Plasmodium falciparum in an area of high ITN coverage in western Kenya.

Methods

Parasite genetic diversity was assessed by scoring eight single copy neutral multilocus microsatellite (MS) markers in samples collected from P. falciparum- infected children (< five years) before introduction of ITNs (1996, baseline, n = 69) and five years after intervention (2001, follow-up, n = 74).

Results

There were no significant changes in overall high mixed infections and unbiased expected heterozygosity between baseline (%MA = 94% and He = 0.75) and follow up (%MA = 95% and He = 0.79) years. However, locus specific analysis detected significant differences for some individual loci between the two time points. Pfg377 loci, a gametocyte-specific MS marker showed significant increase in mixed infections and He in the follow up survey (%MA = 53% and He = 0.57) compared to the baseline (%MA = 30% and He = 0.29). An opposite trend was observed in the erythrocyte binding protein (EBP) MS marker. There was moderate genetic differentiation at the Pfg377 and TAA60 loci (FST = 0.117 and 0.137 respectively) between the baseline and post-ITN parasite populations. Further analysis revealed linkage disequilibrium (LD) of the microsatellites in the baseline (14 significant pair-wise tests and IS A = 0.016) that was broken in the follow up parasite population (6 significant pairs and IS A = 0.0003). The locus specific change in He, the moderate population differentiation and break in LD between the baseline and follow up years suggest an underlying change in population sub-structure despite the stability in the overall genetic diversity and multiple infection levels.

Conclusions

The results from this study suggest that although P. falciparum population maintained an overall stability in genetic diversity after five years of high ITN coverage, there was significant locus specific change associated with gametocytes, marking these for further investigation.

Background

Malaria continues to be a major global public health burden, causing 250 million clinical cases and over 1 million deaths each year. Sub-Saharan Africa accounts for 90% of these cases [1]. To combat malaria, insecticide-treated bed nets (ITNs) have emerged as an efficacious and cost-effective malaria prevention tool. Several previous trials conducted in areas of different malaria transmission patterns have demonstrated that ITNs reduce Plasmodium falciparum malaria transmission by 70-90%. Most importantly, these trials have provided substantial evidence that use of ITNs significantly reduces all-cause mortality and malaria morbidity in children less than five years of age [2]. Additionally, ITNs have been associated with significant reduction in the adverse effects of malaria during pregnancy [3]. The remarkable effectiveness of ITNs has led to an up-scaling of their use in malaria endemic regions in conjunction with other control and prevention measures [4]. Recently, the World Health Organization (WHO) reported that in countries where ITNs have been effectively scaled up, substantial reductions in malaria cases and deaths have occurred [1].

Insecticide-treated bed nets work by killing mosquitoes on contact and also by repelling or deterring the vectors from entering houses, thereby reducing malaria transmission [5]. Thus, the use of ITNs or increased distribution of ITNs not only affects the mosquito populations but also changes the dynamics of parasite dispersion in both human hosts and mosquito vectors, which could in turn modify vector-parasite-host interactions, ultimately affecting parasite populations. Several studies have shown that significant suppression of mosquito populations, changes in species distribution and vector behaviour, and changes in population genetic structure and susceptibility of mosquitoes to insecticides are associated with community-based ITNs intervention [6–9]. However, there are still information gaps on how the reduction in malaria transmission by ITNs affects parasite population genetic structure although there were a few earlier studies that reported no change in the proportion of multiple infections after transmission reduction by use of ITNs and curtains [10, 11].

Human Plasmodium parasites undergo asexual multiplication in the human host and obligate sexual reproduction in the mosquito vector, each stage shaping the parasite population genetic structure. Although the asexual multiplication by haploid parasites in humans is clonal, polymorphism can arise from insertion/deletion of tandem repeats through slippage in the parasite DNA sequences or natural mutations from various pressures in the host-parasite relationship [12–14]. On the other hand, transmission of malaria parasites from human to mosquito, which is solely accomplished by a small number of infective male and female gametocytes generated in humans, creates an opportunity for the generation of parasite diversity and the emergence of novel genetic traits [15]. The parasite sexual reproduction stage in mosquitoes allows for recombination and re-assortment of genetic material between genomes of gametes to form diploid zygotes during the development of oocysts. The degree of inbreeding or outcrossing in mosquitoes influences the number of clones that are infective to human [16]. Numerous factors including host responses can indirectly influence competitive advantages or suppression of specific parasite clones. However, the level of transmission intensity has a direct effect on the number of infected hosts and number of parasite clones per infected individual, which affects parasite population genetic structure in different endemic settings [17, 18]. Therefore, it is important to evaluate whether, and how, the transmission reduction (by use of ITNs or other methods), particularly in high transmission areas, affects the parasite population including the extent of multiple infections, genetic diversity, genes involved in transmission, drug resistance and polymorphism of vaccine candidate genes.

The study of parasite population structure explores the extent of genetic diversity, allele frequency, genotype distribution and degree of genetic admixture among other measures using statistical methods [19]. Common statistical measurements include expected heterozygosity (He) to test genetic variation, linkage disequilibrium (LD) to assess association of alleles between loci, and fixation index (FST) to evaluate population differentiation [20–22]. Natural Plasmodium parasite populations display extensive genetic variability within species at different geographic locations and different transmission intensity levels, with no predominant overall structure. Some studies of P. falciparum population structure report that areas with intense malaria transmission have higher He and higher rates of outcrossing and recombination which breaks LD, resulting in a more panmictic population structure [17, 23, 24]. Such settings allow faster emergence of novel genotypes reflected as multiple infections. The reverse is expected where transmission is lower with consequent lower He, stronger LD, and higher degree of selfing, resulting in a more clonal parasite population structure [17, 25, 26]. Yet other studies in areas with high malaria transmission have observed strong LD and non-random distribution of specific genotypes, implying inbreeding may be more extensive than expected even in areas with perennial transmission [27, 28]. Although the conflicting results generated from different geographic regions could be partially due to the differences in genetic markers used, methods for estimation of allele frequencies or sampling of parasites at different life cycle stages, they underscore the need to study the relationship between the transmission intensity and the P. falciparum population genetic structure in same locality where changes in transmission intensity can be monitored. Interventions which impose reduction in transmission, such as ITNs at high coverage in malaria-holoendemic areas, provide a field experimental system for research on these questions. The information from such studies is also useful in designing molecular surveillance systems for ITNs and for other adjunct control programmes [29, 30].

This study is part of a two-phase large-scale community-based trial conducted in western Kenya and designed to investigate the impact of ITNs on malaria morbidity and all cause mortality. The overall goal of these parasite population genetics studies is to assess the effects of transmission reduction by ITNs on the population genetic structure of P. falciparum parasites for a sustained period. The current study employed eight single copy multilocus neutral microsatellite markers to study the genetic diversity of P. falciparum using blood stage parasites collected from children less than five years old in the same area prior to and five years after the introduction of ITNs. Genetic diversity of the parasites between the baseline and post-ITNs was assessed by quantifying the extent of multiple infections, allele frequencies, He, LD, and genetic differentiation.

Methods

Study area and study samples

This study was part of a two-phase community-based ITN trial conducted in Asembo area (Bondo District) of western Kenya, where malaria is holoendemic. The design and characteristics of the ITN trial have been detailed elsewhere [31–33]. Briefly, biannual population censuses and annual cross-sectional surveys were conducted in 60 villages in the 200 km2 trial area during the rainy season of March to May between 1996 and 2001 to determine the effects of ITNs on malaria morbidity and all cause mortality in children below five years of age. During each cross sectional survey, blood samples were collected from children and parasitological, clinical, and demographic information were documented. Entomological monitoring of Anopheles density and sporozoite infection rates were conducted regularly throughout both phases of the trial. Approximately 98% of malaria infections were due to P. falciparum in the trial area. Before the ITN trial, entomological inoculation rate (EIR) was reported at 61.3 infective bites per person per year [34] and prevalence of parasitaemia was about 70% in children aged less than five years [31]. After the introduction of ITNs, it was estimated that ITNs reduced transmission by 90% at the early stage of the two-phase trials [6]. The EIR and prevalence of parasitaemia in May 2001 were 1.3 infective bites per person per year and about 34% in children less than five years old, respectively [33]. For the study presented here, the samples collected in the 1996 survey, just prior to ITN introduction for baseline measurements, and the samples from the 2001 follow-up survey, five years post-ITN intervention were utilized. Microscopically confirmed malaria positive blood samples were randomly selected in a subset of villages from the 1996 survey and were further matched by the villages in 2001 survey. Calculation of sample size was based on a hypothesized significant difference in overall heterozygosity (He) between baseline and post-ITN parasite populations using a confidence level at 95% and margin of error at 5% [35]. Assuming a He of 0.7 observed in high transmission areas [17] for baseline and a conservative change in He to 0.5 observed in areas with medium transmission intensity [17] for ITN post intervention, a sample size of approximately 63 (+ 10-15%) was deemed adequate allowing for failure in laboratory testing. In total, a sampling frame of 69 samples from baseline and 74 samples from the post-ITN survey was achieved for this study. Parasite genomic DNA extraction from the stored blood samples was by the QIAamp DNA Mini Kit (Qiagen, CA, USA). Extracted DNA was stored at -20C until use. The study protocol was approved by the Ethical Review Committee of the Kenya Medical Research Institute, Nairobi, Kenya, the Institutional Review Board of Michigan State University, and the Institutional Review Board of Centers for Disease Control (CDC) Atlanta, Georgia.

Microsatellite markers and genotyping

Eight single copy microsatellite markers (MS) were used for genotyping as listed in Table 1. Broadly, the microsatellites included: 1) five putatively neutral MS (Poly-α, PfPK2, ADL, TAA60, and TAA 109), 2) one MS (Pfg377) linked to the protein gene exclusively expressed during maturation of gametocytes transmission stage of parasites, and 3) two MS (EBP and P195) linked to the genes of asexual stage antigens under possible natural immune selection [36]. Among these MS markers, Poly-α, PfPK2 and Pfg377 are in the coding region [37]. The selected neutral MS markers were described earlier and have been used previously to study changes in genetic structure of Plasmodium parasites [17, 27, 37]. All amplifications were carried out using single reaction PCR with thermocycling conditions described elsewhere [38]. Fluorescent-labelled primers incorporated with either HEX (green) or FAM (Blue) dyes were used and the PCR products read on ABI (Applied Biosystems 3100) capillary sequencer. GeneMapper software (ABI) was used to automate measurement of microsatellite base-pair length and quantify peak height. Allele identity per locus was obtained directly after allocation of all peaks above 200 fluorescent units. Multiple alleles were quantified based on a method described previously [17] with identification of minor alleles set at peak heights of ≥30% of the predominant allele. Where amplification failed for any of the microsatellites, data was reported as missing and not used for haplotype definition.

Table 1 List of microsatellites and PCR primers used for microsatellites amplification

Data analysis

Since the ITN intervention is one of long-term goals for malaria control programmes, the sampling strategy in the current study aimed at testing changes in genetic diversity of P. falciparum after an extended use of ITN while minimizing any spatial effects on the population. Potential biases were envisioned in the data analyses as allele frequencies might change over time in finite populations [39]. However, the prediction in this study was that the ITN-mediated transmission reduction would precipitate changes in the parasite population from panmictic (higher parasite diversity, expected for high transmission intensity) to a more clonal structure (lower parasite diversity, expected for low transmission intensity). Changes in genetic diversity of the parasites between the baseline and post-ITNs were, therefore, assessed by quantifying and comparing multiple infections, allele frequencies, He, LD, and FST.

Initial microsatellite data checking and data conversion was done using Excel Microsatellite Tool Kit, an add-in programme used to format raw microsatellite data in Microsoft Excel for consequent use in different genetics softwares [40]. Analysis of multiple infections was based on both predominant and minor alleles. Outcome measures for this analysis were the proportion of infections with more than one allele, the mean allele count for individual microsatellite loci, the overall proportion of infections with at least two alleles and the overall mean of the highest number of allele count detected by any of the microsatellites. The difference in proportions of infections with at least two alleles and the difference in the mean allele count between the baseline and post-intervention parasite populations were assessed using Pearson's chi-square and Wilcoxon tests respectively. Conversely, only the predominant allele defined by the highest peak in each electropherogram per locus was used to analyze the allele frequency and allele richness, He, LD, and FST. There are possible biases from using predominant allele techniques to determine allele frequency in multiple infections [41, 42]. However, the method applied here was shown to be appropriate in previous studies using some of the microsatellite targets selected in this study [17, 37, 43, 44].

The infinite allele model which is more appropriate for analysis of the complex patterns observed in microsatellite loci in P. falciparum was used for genetic analysis [45]. Allele number, frequency and richness per locus in each parasite population were obtained by FSTAT2 [46]. The measure for allele number, frequency and richness was to illustrate the composition and distribution of alleles in the population. Unbiased He at each locus was calculated as He = [n/(n-1)][1-Σni = 1p2 i ] where n is the number of isolates sampled and p i is the frequency of the i th allele while sampling variance for He (Vs(He)) was calculated as Vs(He) = 2(n-1)/n3[2(n-2)][Σpi3-(Σpi2)2] [17, 22]. The difference in single locus heterozygosities between the two parasite populations was tested with standard error (SE) of the sampling variance by the method described earlier [22] using the z absolute values to obtain the p-values.

The prediction of the ITN-mediated decrease in transmission intensity and consequent reduction in genetic diversity was a possible increase in LD in the post-ITN parasite population. LD measures the degree of association between or among gene loci under the null hypothesis of no association [47]. As such,

individual pair-wise LDs for baseline and post-ITN parasite populations were first obtained respectively by the Fisher's exact test adapted for haploid data using ARLEQUINS 3.11 programme [48]. Multilocus LDs were further assessed from the overall index of association (IS A ) using LIAN programme [49] for both baseline and post-ITNs respectively. The multilocus LD test measures non-random association among all loci. The test compares the variance of differences at LD with the variance expected in LD derived from 10,000 simulated data sets, H0: VD = Ve. A significant LD was when the observed variance (VD) was greater than expected (Ve). Index of association was expressed as IS A = (VD/Ve-1)(r-1) where r is the number of loci tested. A 95% confidence limit was determined by Monte Carlo simulation.

In order to further assess genetic diversity before and after ITNs intervention, genetic differentiation was tested using the FST estimator [50] implemented by FSTAT programme [46]. FST is a comparison of the sum of genetic variability within and between populations based on allele frequency differences in populations. Interpretation of the FST values at each locus was based on three categories defined earlier as no differentiation (0), low genetic differentiation (0 > 0.05), moderate differentiation (0.05-0.15) and great differentiation (0.15-0.25) [51].

When multiple tests were conducted, significant levels of p-value for the comparisons were adjusted using Bonferroni's correction [52].

Results

Results in this study comprise outcomes from the analysis of the eight microsatellite loci on 69 and 74 P. falciparum positive samples from baseline and post-ITN surveys, respectively. Locus P195 had an overall amplification rate of 90% for the baseline parasite population. All other loci had amplification rates of 94-100%.

Multiple infections

The extent of multiple infections was assessed based on the proportion of multiple alleles and mean allele counts (Table 2). Overall, more than 90% of samples from both baseline and post-ITN surveys had at least two or more alleles detected by any of the microsatellites targets. There was no significant difference in the overall proportions of multiple alleles between the baseline and post-ITN parasite populations. However, analysis of individual loci revealed that the Pfg377 marker had a significantly lower proportion of multiple alleles in the baseline (30%) compared to post-ITN parasite population (53%). All other loci showed similar proportions of multiple alleles in the two parasite populations.

Table 2 Comparison of proportion of multiple alleles and mean allele counts in the baseline and post-ITN parasite populations

Results from the mean allele counts showed that PfPK2 had significantly higher allele counts in the baseline compared to that in the post-ITN survey. In contrast, Pfg377 had a relatively higher mean allele count (1.62 ± 0.09) in the post-ITN survey than in the baseline parasite population (1.36 ± 0.07) although it was not a statistically significant difference (p = 0.0076) after Bonferroni correction. There was no significant difference in the mean allele counts for all other microsatellite markers between the two parasite populations. Overall, mean allele count was similar in the baseline parasite population (3.4 ± 0.15) compared to the post-ITN parasite population (3.1 ± 0.12) (Table 2).

Allele frequency and heterozygosity

The allele size, number, frequency and richness for each MS shown in Figure 1 and Table 3 illustrated the composition and distribution of alleles. There were no differences between the allele number and richness as sample numbers used in the baseline and post-ITN surveys did not differ significantly for each locus. The allele number per locus ranged from a minimum of 4 (Pfg377) to maximum of 18 (Poly-α) in the baseline parasite population and from 3 (Pfg377) to 17 (Poly-α) in the post-ITN parasite population. This suggests marked variation between loci in the baseline and post-ITN surveys and was further tested by He. Overall, He was high and similar in the baseline (0.75 ± 0.072) and post-ITN parasite populations (0.79 ± 0.038). However, He of individual loci showed extensive range in gene diversity (Table 3). Poly-α locus had the highest level of He while Pfg377 locus showed the lowest in both parasite populations. Notably, Pfg377 locus had significantly lower He in the baseline parasite population compared to the He in the post-ITN survey (p = 0.044).

Figure 1
figure 1

Comparison of allele size and composition (base pairs, X-axis) and frequency distribution (Y-axis) of the eight individual microsatellites in the baseline (blue) and post-ITN (red) surveys. The standardized Y-axis scale was used to depict the proportion of different alleles in each locus.

Table 3 Comparisons of genetic diversity between the baseline and post-ITN parasite populations

Linkage disequilibrium and genetic differentiation

Pair-wise LD for each individual MS locus for the baseline and post-ITN parasite populations was assessed. Of the 28 possible pair-wise tests LD was highly significant in a total of 14 pairs in the baseline parasite population while only 6 pairs were significant in the post-ITN population (p = 0.0018) as shown in Table 4. Specifically, in the baseline survey MS P195 locus (Chr9) and TAA60 (Chr9) had strong LD with loci located on different chromosomes including Poly-α (Chr4), Pfg377 (Chr12), PfPK2 (Chr12), EBP (Chr13) and ADL (Chr13). This association was however, reduced in the post-ITN population where 6 pairs showed significant LD. Only LD between P195 and ADL, and between TAA60 and EBP loci were maintained while the remaining four loci pairs were new associations in the post-ITN survey.

Table 4 Estimates of pair-wise linkage disequilibrium (LD) in the baseline and post-ITN parasite populations

The multilocus analysis showed significant LD in the baseline population with a variance in difference (VD) of 1.21 and variance expected (Ve) of 1.10 (p = 0.01) compared to non-significant LD in the post-ITN population where VD and Ve were 1.183 and 1.180 respectively (p = 0.57) as shown in Table 5. Consequently the index of association (IS A ) was significant in the baseline survey (0.016) compared to the post-ITN survey (0.0003). The multilocus LD results coupled with the significant pair-wise LD observed in individual microsatellites suggest that the existing non-random association between the MS loci in the baseline was broken in the post-ITN parasite population.

Table 5 Estimates of multilocus linkage disequilibrium (LD) for baseline and post-ITN parasite populations

Results of the overall and locus specific genetic differentiation between the baseline and post-ITN surveys are shown in Table 6. Overall FST was low at 0.027. In the single locus FST, both Pfg377 and TAA60 markers showed moderate genetic differentiation (FST = 0.117 and 0.137) respectively while all other remaining six loci had little differentiation. The results suggest that a moderately significant genetic variability at Pfg377 and TAA60 arose from differences between the two parasite populations while the majority of genetic variability remains within each population. Conversely, low genetic differentiation in the other six markers suggested much of the genetic variability resulted from within each parasite population.

Table 6 Genetic differentiation index (FST) between baseline and post-ITN parasite populations

Discussion

The effect of five years of high coverage with ITNs on the genetic diversity of P. falciparum parasites was examined in this study. The overall proportion of mixed infections and heterozygosity were high at over 90% and 0.75 respectively both before and after ITN use with no significant reduction in these two parameters as well as the overall mean allele counts. This indicates an extensive multiplicity of circulating parasites in the area in spite of a dramatic decline in EIR post-ITN intervention. The results from this study are consistent with those of an earlier study conducted in areas with EIRs ranging from 0.4 to 31.8 in western Kenya, which also recorded over 80% mixed infections in both low and high malaria transmission areas [44]. This suggests the presence of a steady mix of circulating Plasmodium parasites in western Kenya despite reduction in EIRs.

The stable overall genetic diversity after dramatic reduction in transmission intensity observed in the current study was unexpected by the initial prediction. The counter-intuitive results suggest that other factors may be involved in offsetting the effect of transmission reduction on parasite genetic diversity and/or stabilization of the overall genetic diversity of malaria parasite. Indeed, several previous studies suggest that genetic diversity of malaria parasite measured by different markers could be shaped directly or indirectly by multiple factors such as seasonality, geographic scale, migration, disease severity, and host age and immunity [53–55] in addition to transmission intensity per se and natural selection. To minimize variation in host age, seasonality and geographic scale between baseline and post-ITN surveys, the current study sampled children less than five years of age, during similar transmission seasons and matched by villages in the two surveys. In addition, previous studies conducted in the ITN trial area showed that transmission reduction by use of ITNs changed humoral immunity in children and reduced childhood malaria morbidity and infant mortality resulting in overall decreased anti-malarial treatment [33, 56, 57]. The clinical and immunological outcomes after ITN intervention in the study area could potentially counteract the effect of transmission reduction on parasite genetic diversity and/or sustain overall higher genetic diversity although what mechanisms govern such a process within hosts is unknown. It is also possible that gene flow due to migration of mosquitoes and humans from surrounding non-ITN trial areas might contribute to the overall unchanged genetic diversity. However, the current study was not able to quantify the gene flow as the original ITN trial was not designed to include surrounding non-ITN areas after five years post-ITNs for comparison. Considering the parasite diversity could be influenced by multiple factors listed above, detection of change in parasite diversity within five years time window in the current study might not be sufficient. Currently, further studies on parasite population genetics are ongoing, which includes analysis of samples from approximately a decade later in the same ITN trial area and surrounding areas as well. Taken together, the unchanged overall genetic diversity observed in this study suggests a strong resilience of malaria parasite in response to dramatic transmission reduction after five years of sustained ITN use and possible involvement of other factors in stabilizing the overall parasite genetic diversity.

While the overall stability in the parasite genetic diversity show the transmission reduction by ITNs had insignificant impact on parasite population, locus specific changes suggest there were some differences in the parasite population sub-structure. For example, PfPK2 microsatellite marker showed a decrease in the mean allele counts in the post-ITN survey, while Pfg377 microsatellite locus showed a significant increase in the proportion of infections with more than one allele. There was also a decrease in genetic diversity (He) in the EBP marker, but an increase at the Pfg377 locus. EBP MS locus flanks the functionally important erythrocyte binding protein gene expressed in the asexual stage of the life cycle and the gene may be under selection by human immune response [58]. It is possible that the decrease in He for EBP MS marker observed in the post-ITN survey could reflect an indirect effect of ITNs on parasite genetic diversity but this will need further investigation. Likewise, PfPK2 MS which showed a decrease in mean allele counts flanks a putative protein kinase gene expressed in young trophozoite although the exact function is still not clear [59]. On the contrary, the MS located in the coding region of Pfg377 antigen gene specific for 'gametocyte-producing' parasites [60] showed an increased diversity in the post-ITN parasite population. Because there was either decrease or no significant changes in genetic diversity in other MS, the increase in genetic diversity for Pfg377 locus in the post-ITN parasites most likely reflected selection rather than genetic drift. Interestingly, the gametocyte carriage in the ITN trial area was significantly lower in the baseline survey (proportion 17% and mean density 12.4/ul) compared to the five years post-ITN survey (proportion 23%, mean density 41.2/ul) (CDC unpublished data). Taken together, this suggests that there is a possible relationship between the increased genetic diversity of Pfg377 and an increase in gametocyte carriage. The increased gametocytaemia and genetic diversity of Pfg377 locus could be an adaptive mechanism for transmission reduction to enhance the potential for parasite transmission to mosquitoes to maintain the life cycle for survival. This hypothesis requires testing to assess whether this gene has been a target of selection.

Consistent with the genetic diversity data described above, overall genetic differentiation between the baseline and post-ITN parasite populations was low, mainly arising from variations in the Pfg377 and TAA60 microsatellite markers. The differentiation observed in this study for Pfg377 and TAA60 were higher between the baseline and post-ITN surveys than that observed between three geographically different areas in western Kenya [44]. Yet FST was much lower in our study at other remaining loci examined in the same study in western Kenya [44]. The lower FST estimates at other loci observed in this study are expected since the samples were from the same area for baseline and post-ITN surveys. It is possible that the differentiation at TAA60 could represent random temporal effect/drift on the parasite population independent of transmission reduction. However, the differentiation observed at Pfg377 most likely resulted from the decreased transmission intensity by the use of ITNs rather than mere temporal effect since Pfg377 locus showed consistent increases in multiple infection and He after ITN intervention.

The inter-relationship among LD, transmission intensity and genetic diversity of malaria parasites is complex and is still far from conclusive. The stronger LD observed in the baseline survey in our study area is consistent with the trend observed in previous studies conducted in the Democratic Republic of Congo, Zimbabwe and western Kenya lowland areas where malaria transmission is intense [17, 27, 28, 44], suggesting the occurrence of high inbreeding in P. falciparum even in areas with intense and perennial transmission. It would be expected that decreasing transmission intensity by use of ITNs increases LD level based on a generalized assertion of higher LD in low transmission areas [17]. However, the results from five years post-ITN intervention in this study were surprising and interesting. After ITN intervention the pair-wise LDs were broken in 65% of physically unlinked loci (Table 4) and the multilocus LD was also not significant compared to the baseline survey (Table 5). The LD result from post-ITN parasite population could suggest that the overall parasite population became more panmictic after bed net intervention, which is contrary to earlier prediction of a more clonal structure after an ITN mediated transmission reduction. However, it is also possible that the unexpected decrease in LD in the post-ITN parasite population is partially masked by the increase in genetic diversity of Pfg377, the 'gametocyte specific' MS marker, but this will require further investigation.

Conclusion

This study suggests that although the parasite population maintained an overall stability after ITN use, there were locus specific changes in the P. falciparum parasites contributing to the observed differentiation between the two parasite populations. Of note, the data on Pfg377 locus showed an increase in diversity ecologically associated with reduction in transmission intensity. Further studies are necessary to evaluate the usefulness of this marker and other gametocyte-specific gene markers as molecular tools for monitoring how changes in transmission reflect gametocyte population dynamics [61]. It is also important to monitor genetic structure of P. falciparum for extended periods, and in different geographic areas and in changing ITN coverage.

References

  1. WHO: World Malaria Report. WHO Library Cataloguing-in-Publication Data. 2008, 25-

    Google Scholar 

  2. Lengeler C: Insecticide-treated bed nets and curtains for preventing malaria. Cochrane Database Syst Rev. 2004, CD000363-

    Google Scholar 

  3. Gamble C, Ekwaru JP, ter Kuile FO: Insecticide-treated nets for preventing malaria in pregnancy. Cochrane Database of Systematic Reviews. 2006, 19 (2): CD003755-

    Google Scholar 

  4. Hill J, Lines J, Rowland M, David HM: Insecticide-Treated Nets. Adv Parasitol. 2006, Academic Press, 61: 77-128. 10.1016/S0065-308X(05)61003-2.

    Google Scholar 

  5. Curtis CF, Jana-Kara B, Maxwell CA: Insecticide treated nets: impact on vector populations and relevance of initial intensity of transmission and pyrethroid resistance. J Vector Borne Dis. 2003, 40: 1-8.

    CAS  PubMed  Google Scholar 

  6. Gimnig JE, Kolczak MS, Hightower AW, Vulule JM, Schoute E, Kamau L, Phillips-Howard PA, ter Kuile FO, Nahlen BL, Hawley WA: Effect of permethrin-treated bed nets on the spatial distribution of malaria vectors in western Kenya. Am J Trop Med Hyg. 2003, 68: 115-120.

    PubMed  Google Scholar 

  7. Hossain MI, Curtis CF: Permethrin-impregnated bednets: behavioural and killing effects on mosquitoes. Med Vet Entomol. 1989, 4: 367-376. 10.1111/j.1365-2915.1989.tb00243.x.

    Article  Google Scholar 

  8. Vulule JM, Beach RF, Atieli FK, Roberts JM, Mount DL, Mwangi RW: Reduced susceptibility of Anopheles gambiae to permethrin associated with the use of permethrin-impregnated bednets and curtains in Kenya. Med Vet Entomol. 1994, 8: 71-75. 10.1111/j.1365-2915.1994.tb00389.x.

    Article  CAS  PubMed  Google Scholar 

  9. Lindblade KA, Gimnig JE, Kamau L, Hawley WA, Odhiambo F, Olang G, Ter Kuile FO, Vulule JM, Slutsker L: Impact of sustained use of insecticide-treated bednets on malaria vector species distribution and culicine mosquitoes. J Med Entomol. 2006, 43: 428-432. 10.1603/0022-2585(2006)043[0428:IOSUOI]2.0.CO;2.

    Article  CAS  PubMed  Google Scholar 

  10. Bolad A, Nebie I, Esposito F, Berzins K: The use of impregnated curtains does not affect antibody responses against Plasmodium falciparum and complexity of infecting parasite populations in children from Burkina Faso. Acta Trop. 2004, 90: 237-247. 10.1016/j.actatropica.2003.12.004.

    Article  CAS  PubMed  Google Scholar 

  11. Smith T, Felger I, Fraser-Hurt N, Beck HP: Effect of insecticide-treated bed nets on the dynamics of multiple Plasmodium falciparum infections. Trans R Soc Trop Med Hyg. 1999, 93 (Suppl 1): 53-57. 10.1016/S0035-9203(99)90328-0.

    Article  PubMed  Google Scholar 

  12. Carvalho TG, Menard R: Manipulating the Plasmodium genome. Curr Issues Mol Biol. 2005, 7: 39-55.

    CAS  PubMed  Google Scholar 

  13. Rich SM, Ayala FJ: Population structure and recent evolution of Plasmodium falciparum. Proc Natl Acad Sci USA. 2000, 97: 6994-7001. 10.1073/pnas.97.13.6994.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Ambroise-Thomas P: [Malaria, genetics and molecular biology: myths, hopes and realities]. Annales pharmaceutiques françaises. 2001, 59: 93-104.

    CAS  PubMed  Google Scholar 

  15. Read AF, Anwar M, Shutler D, Nee S: Sex allocation and population structure in malaria and related parasitic protozoa. Proc Biol Sci. 1995, 260: 359-363. 10.1098/rspb.1995.0105.

    Article  CAS  PubMed  Google Scholar 

  16. Day KP, Koella JC, Nee S, Gupta S, Read AF: Population genetics and dynamics of Plasmodium falciparum: an ecological view. Parasitology. 1992, 104 (Suppl): S35-52. 10.1017/S0031182000075235.

    Article  PubMed  Google Scholar 

  17. Anderson TJ, Haubold B, Williams JT, Estrada-Franco JG, Richardson L, Mollinedo R, Bockarie M, Mokili J, Mharakurwa S, French N, Whitworth J, Velez ID, Brockman AH, Nosten F, Ferreira MU, Day KP: Microsatellite markers reveal a spectrum of population structures in the malaria parasite Plasmodium falciparum. Mol Biol Evol. 2000, 17: 1467-1482.

    Article  CAS  PubMed  Google Scholar 

  18. Wargo AR, de Roode JC, Huijben S, Drew DR, Read AF: Transmission stage investment of malaria parasites in response to in-host competition. Proc Biol Sci. 2007, 274: 2629-2638. 10.1098/rspb.2007.0873.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Tibayrenc M: Genetic epidemiology of parasitic protozoa and other infectious agents: the need for an integrated approach. Int J Parasitol. 1998, 28: 81-104. 10.1016/S0020-7519(97)00180-X.

    Google Scholar 

  20. Nei M: Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics. 1978, 89: 583-590.

    PubMed Central  CAS  PubMed  Google Scholar 

  21. Cockerham CC, Weir BS: Estimation of inbreeding parameters in stratified populations. Ann Hum Genet. 1986, 50: 271-281. 10.1111/j.1469-1809.1986.tb01048.x.

    Article  CAS  PubMed  Google Scholar 

  22. Nei M: Molecular evolutionary genetics, Genetic variation within species. 1987, New York: Columbia University Press

    Google Scholar 

  23. Gauthier C, Tibayrenc M: Population structure of malaria parasites: The driving epidemiological forces. Acta Trop. 2005, 94: 241-250.

    Article  CAS  PubMed  Google Scholar 

  24. Mu J, Awadalla P, Duan J, McGee KM, Joy DA, McVean GA, Su XZ: Recombination hotspots and population structure in Plasmodium falciparum. PLoS Biol. 2005, 3: e335-10.1371/journal.pbio.0030335.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Babiker HA, Walliker D: Current views on the population structure of Plasmodium falciparum: Implications for control. Parasitol Today. 1997, 13: 262-267. 10.1016/S0169-4758(97)01075-2.

    Article  CAS  PubMed  Google Scholar 

  26. Urdaneta L, Lal A, Barnabe C, Oury B, Goldman I, Ayala FJ, Tibayrenc M: Evidence for clonal propagation in natural isolates of Plasmodium falciparum from Venezuela. Proc Natl Acad Sci USA. 2001, 98: 6725-6729. 10.1073/pnas.111144998.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Durand P, Michalakis Y, Cestier S, Oury B, Leclerc MC, Tibayrenc M, Renaud F: Significant linkage disequilibrium and high genetic diversity in a population of Plasmodium falciparum from an area (Republic of the Congo) highly endemic for malaria. Am J Trop Med Hyg. 2003, 68: 345-349.

    CAS  PubMed  Google Scholar 

  28. Razakandrainibe FG, Durand P, Koella JC, De Meeus T, Rousset F, Ayala FJ, Renaud F: "Clonal" population structure of the malaria agent Plasmodium falciparum in high-infection regions. Proc Natl Acad Sci USA. 2005, 102: 17388-17393. 10.1073/pnas.0508871102.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Escalante AA, Smith DL, Kim Y: The dynamics of mutations associated with anti-malarial drug resistance in Plasmodium falciparum. Trends Parasitol. 2009, 25: 557-563. 10.1016/j.pt.2009.09.008.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Takala SL, Plowe CV: Genetic diversity and malaria vaccine design, testing and efficacy: preventing and overcoming 'vaccine resistant malaria'. Parasite Immunol. 2009, 31: 560-573. 10.1111/j.1365-3024.2009.01138.x.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. ter Kuile FO, Terlouw DJ, Kariuki SK, Phillips-Howard PA, Mirel LB, Hawley WA, Friedman JF, Shi YP, Kolczak MS, Lal AA, Vulule JM, Nahlen BL: Impact of permethrin-treated bed nets on malaria, anemia, and growth in infants in an area of intense perennial malaria transmission in western Kenya. Am J Trop Med Hyg. 2003, 68: 68-77.

    PubMed  Google Scholar 

  32. Phillips-Howard PA, Nahlen BL, Kolczak MS, Hightower AW, ter Kuile FO, Alaii JA, Gimnig JE, Arudo J, Vulule JM, Odhacha A, Kachur SP, Schoute E, Rosen DH, Sexton JD, Oloo AJ, Hawley WA: Efficacy of permethrin-treated bed nets in the prevention of mortality in young children in an area of high perennial malaria transmission in western Kenya. Am J Trop Med Hyg. 2003, 68: 23-29.

    PubMed  Google Scholar 

  33. Lindblade KA, Eisele TP, Gimnig JE, Alaii JA, Odhiambo F, ter Kuile FO, Hawley WA, Wannemuehler KA, Phillips-Howard PA, Rosen DH, Nahlen BL, Terlouw DJ, Adazu K, Vulule JM, Slutsker L: Sustainability of reductions in malaria transmission and infant mortality in western Kenya with use of insecticide-treated bednets: 4 to 6 years of follow-up. JAMA. 2004, 291: 2571-2580. 10.1001/jama.291.21.2571.

    Article  CAS  PubMed  Google Scholar 

  34. Singer LM, Mirel LB, ter Kuile FO, Branch OH, Vulule JM, Kolczak MS, Hawley WA, Kariuki SK, Kaslow DC, Lanar DE, Lal AA: The effects of varying exposure to malaria transmission on development of antimalarial antibody responses in preschool children. XVI. Asembo Bay Cohort Project. J Infect Dis. 2003, 187: 1756-1764. 10.1086/375241.

    Article  PubMed  Google Scholar 

  35. Kirkwood BR: Essentials of medical statistics. 2001, Blackwell Science

    Google Scholar 

  36. Su X, Wellems TE: Toward a high-resolution Plasmodium falciparum linkage map: polymorphic markers from hundreds of simple sequence repeats. Genomics. 1996, 33: 430-444. 10.1006/geno.1996.0218.

    Article  CAS  PubMed  Google Scholar 

  37. Anderson TJ, Su XZ, Bockarie M, Lagog M, Day KP: Twelve microsatellite markers for characterization of Plasmodium falciparum from finger-prick blood samples. Parasitology. 1999, 119 (Pt 2): 113-125. 10.1017/S0031182099004552.

    Article  CAS  PubMed  Google Scholar 

  38. Nair S, Williams JT, Brockman A, Paiphun L, Mayxay M, Newton PN, Guthmann JP, Smithuis FM, Hien TT, White NJ, Nosten F, Anderson TJ: A selective sweep driven by pyrimethamine treatment in southeast asian malaria parasites. Mol Biol Evol. 2003, 20: 1526-1536. 10.1093/molbev/msg162.

    Article  CAS  PubMed  Google Scholar 

  39. Heath DD, Busch C, Kelly J, Atagi DY: Temporal change in genetic structure and effective population size in steelhead trout (Oncorhynchus mykiss). Mol Ecol. 2002, 11: 197-214. 10.1046/j.1365-294X.2002.01434.x.

    Article  CAS  PubMed  Google Scholar 

  40. Park SDE: Trypanotolerance in West African Cattle and the Population Genetic Effects of Selection. PhD Thesis. 2001, Internet Download

    Google Scholar 

  41. Holzmuller P, Herder S, Cuny G, De Meeus T: From clonal to sexual: a step in T. congolense evolution?. Trends Parasitol. 2010, 26: 56-60. 10.1016/j.pt.2009.11.006.

    Article  PubMed  Google Scholar 

  42. Hastings IM, Smith TA: MalHaploFreq: a computer programme for estimating malaria haplotype frequencies from blood samples. Malar J. 2008, 7: 130-10.1186/1475-2875-7-130.

    Article  PubMed Central  PubMed  Google Scholar 

  43. McCollum AM, Poe AC, Hamel M, Huber C, Zhou Z, Shi YP, Ouma P, Vulule J, Bloland P, Slutsker L, Barnwell JW, Udhayakumar V, Escalante AA: Antifolate resistance in Plasmodium falciparum: multiple origins and identification of novel dhfr alleles. J Infect Dis. 2006, 194: 189-197. 10.1086/504687.

    Article  CAS  PubMed  Google Scholar 

  44. Zhong D, Afrane Y, Githeko A, Yang Z, Cui L, Menge DM, Temu EA, Yan G: Plasmodium falciparum genetic diversity in western Kenya highlands. Am J Trop Med Hyg. 2007, 77: 1043-1050.

    PubMed  Google Scholar 

  45. Anderson TJ, Su XZ, Roddam A, Day KP: Complex mutations in a high proportion of microsatellite loci from the protozoan parasite Plasmodium falciparum. Mol Ecol. 2000, 9: 1599-1608. 10.1046/j.1365-294x.2000.01057.x.

    Article  CAS  PubMed  Google Scholar 

  46. Goudet J: FSTAT (version 1.2): a computer program to calculate F-statistics. Journal of heredity. 1995, 86: 485-

    Google Scholar 

  47. Slatkin M: Linkage disequilibrium in growing and stable populations. Genetics. 1994, 137: 331-336.

    PubMed Central  CAS  PubMed  Google Scholar 

  48. Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online. 2005, 1: 47-50.

    PubMed Central  CAS  Google Scholar 

  49. Haubold B, Hudson RR: LIAN 3.0: detecting linkage disequilibrium in multilocus data. Linkage Analysis. Bioinformatics. 2000, 9: 847-848. 10.1093/bioinformatics/16.9.847.

    Article  Google Scholar 

  50. Weir BS: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-10.2307/2408641.

    Article  Google Scholar 

  51. Balloux F, Lugon-Moulin N: The estimation of population differentiation with microsatellite markers. Mol Ecol. 2002, 11: 155-165. 10.1046/j.0962-1083.2001.01436.x.

    Article  PubMed  Google Scholar 

  52. Rice W: Analyzing Tables of Statistical Tests. Evolution. 1989, 43: 223-225. 10.2307/2409177.

    Article  Google Scholar 

  53. Babiker HA, Abdel-Muhsin AM, Ranford-Cartwright LC, Satti G, Walliker D: Characteristics of Plasmodium falciparum parasites that survive the lengthy dry season in eastern Sudan where malaria transmission is markedly seasonal. Am J Trop Med Hyg. 1998, 59: 582-590.

    CAS  PubMed  Google Scholar 

  54. Susomboon P, Iwagami M, Tangpukdee N, Krusood S, Looareesuwan S, Kano S: Differences in genetic population structures of Plasmodium falciparum isolates from patients along Thai-Myanmar border with severe or uncomplicated malaria. Malar J. 2008, 7: 212-10.1186/1475-2875-7-212.

    Article  PubMed Central  PubMed  Google Scholar 

  55. Weedall GD, Conway DJ: Detecting signatures of balancing selection to identify targets of anti-parasite immunity. Trends Parasitol. 26: 363-369. 10.1016/j.pt.2010.04.002.

  56. Kariuki SK, Lal AA, Terlouw DJ, ter Kuile FO, Ong'echa JM, Phillips-Howard PA, Orago AS, Kolczak MS, Hawley WA, Nahlen BL, Shi YP: Effects of permethrin-treated bed nets on immunity to malaria in western Kenya II. Antibody responses in young children in an area of intense malaria transmission. Am J Trop Med Hyg. 2003, 68: 108-114.

    PubMed  Google Scholar 

  57. ter Kuile FO, Terlouw DJ, Phillips-Howard PA, Hawley WA, Friedman JF, Kolczak MS, Kariuki SK, Shi YP, Kwena AM, Vulule JM, Nahlen BL: Impact of permethrin-treated bed nets on malaria and all-cause morbidity in young children in an area of intense perennial malaria transmission in western Kenya: cross-sectional survey. Am J Trop Med Hyg. 2003, 68: 100-107.

    PubMed  Google Scholar 

  58. Baum J, Thomas AW, Conway DJ: Evidence for diversifying selection on erythrocyte-binding antigens of Plasmodium falciparum and P. vivax. Genetics. 2003, 163: 1327-1336.

    PubMed Central  CAS  PubMed  Google Scholar 

  59. Zhao Y, Kappes B, Yang J, Franklin RM: Molecular cloning, stage-specific expression and cellular distribution of a putative protein kinase from Plasmodium falciparum. Eur J Biochem. 1992, 207: 305-313. 10.1111/j.1432-1033.1992.tb17051.x.

    Article  CAS  PubMed  Google Scholar 

  60. Alano P, Read D, Bruce M, Aikawa M, Kaido T, Tegoshi T, Bhatti S, Smith DK, Luo C, Hansra S, Carter R, Elliott JF: COS cell expression cloning of Pfg377, a Plasmodium falciparum gametocyte antigen associated with osmiophilic bodies. Mol Biochem Parasitol. 1995, 74: 143-156. 10.1016/0166-6851(95)02491-3.

    Article  CAS  PubMed  Google Scholar 

  61. Alano P: Molecular approaches to monitor parasite genetic complexity in the transmission of Plasmodium falciparum malaria. Parassitologia. 2005, 47: 199-203.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors express their gratitude to the children and families who participated in the ITN trial. We also thank Dr Paula Marcet of Malaria Branch, Division of Parasitic Diseases, CDC Atlanta for her tremendous assistance in genetic data analysis and Dr Ananias Escalante of Arizona State University for providing advices on genetic data analysis and valuable suggestions on the manuscript. We thank the Director, Kenya Medical Research Institute for permission to publish this paper. This study was supported by the Multilateral Initiative on Malaria grant # A40046 through the WHO Special Programme for Research and Training in Tropical Diseases and was partially supported by U.S. National Science Foundation, Ecology of Infectious Diseases grant # EF-0723770.

Disclaimer: The opinions and assertions contained in this manuscript are those of authors and do not necessarily represent the view of the Centers for Disease Control and Prevention.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ya Ping Shi.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

WG carried out genotyping work and genetic data analysis, and wrote the manuscript. SK was responsible for sample processing and laboratory diagnosis for the ITN trial, and participated in the design of this study. WH, FTK, DT, PPH, JG and KL implemented and conducted the Kenya ITN trial including collection of samples and epidemiological data. WH, BN and KL were PIs for the two phase ITN trial. JG, EW, MH and LS participated in the design of this study. SC and JW assisted in sample size calculation, sample selection and data analysis. YPS was the PI and responsible for the design of this study, participated in data analysis, and wrote manuscript. All authors contributed to data interpretation, read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Gatei, W., Kariuki, S., Hawley, W. et al. Effects of transmission reduction by insecticide-treated bed nets (ITNs) on parasite genetics population structure: I. The genetic diversity of Plasmodium falciparum parasites by microsatellite markers in western Kenya. Malar J 9, 353 (2010). https://doi.org/10.1186/1475-2875-9-353

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-2875-9-353

Keywords