Skip to main content

Genetic diversity and population structure of Plasmodium falciparum in the Philippines

Abstract

Background

In the Philippines, malaria morbidity and mortality have decreased since the 1990s by effective malaria control. Several epidemiological surveys have been performed in the country, but the characteristics of the Plasmodium falciparum populations are not yet fully understood. In this study, the genetic structure of P. falciparum populations in the Philippines was examined.

Methods

Population genetic analyses based on polymorphisms of 10 microsatellite loci of the parasite were conducted on 92 isolates from three provinces (Kalinga, Palawan, and Davao del Norte) with different malaria endemicity.

Results

The levels of genetic diversity and the effective population sizes of P. falciparum in the Philippines were similar to those reported in the mainland of Southeast Asia or South America. In the low malaria transmission area (Kalinga), there was a low level of genetic diversity and a strong linkage disequilibrium (LD) when the single-clone haplotype (SCH) was used in the multilocus LD analysis, while in the high malaria transmission areas (Palawan and Davao del Norte), there was a high level of genetic diversity and a weak LD when SCH was used in the multilocus LD analysis. On the other hand, when the unique haplotypes were used in the multilocus LD analysis, no significant LD was observed in the Kalinga and the Palawan populations. The Kalinga and the Palawan populations were, therefore, estimated to have an epidemic population structure. The three populations were moderately differentiated from each other.

Conclusion

In each area, the level of genetic diversity correlates with the local malaria endemicity. These findings confirm that population genetic analyses using microsatellite loci are a useful tool for evaluating malaria endemicity.

Background

Malaria is still one of the major public health problems in the Philippines, although the morbidity and the mortality have decreased in the last decade [1, 2]. The Philippine archipelago comprises 7,109 islands in the western Pacific Ocean. Out of the 79 provinces, malaria is endemic to 65 provinces [3]. Epidemiological data suggest that malaria endemicity across the Philippines is quite different from one area to another. Therefore, it is likely that the parasite populations in each endemic area differ genetically.

Understanding the genetic structure of malaria parasites is essential to predict how fast phenotypes of interest, such as novel antigenic variants or drug resistance, originate and spread in populations [4, 5]. Anderson et al reported a spectrum of population structures in Plasmodium falciparum: in high transmission areas, such as sub-Saharan Africa, weak linkage disequilibrium, high genetic diversity and low levels of genetic differentiation between populations were observed, while in low transmission areas, such as the Brazilian Amazon, significant linkage disequilibrium, low genetic diversity and high levels of genetic differentiation between populations were observed [6].

Population genetic analyses of P. falciparum in Southeast Asia have been reported from Thailand and Malaysian Borneo [6, 7]. In the Philippines, several epidemiological surveys and some molecular biological studies on drug resistant malaria and antigenic molecules have been reported [1, 2, 8–11]. However, basic genetic analyses of the parasite populations in the country using multilocus neutral markers have not been reported so far.

The objective of the present study was to estimate the genetic structure of P. falciparum populations in the Philippines using 10 highly polymorphic microsatellite loci and to discuss the correlation between the levels of genetic diversity and the level of malaria endemicity in this country. To investigate P. falciparum populations in the Philippines, the population genetic analyses were conducted using a method consistent with those used to assess the genetic structure of the parasite populations in other countries, so as to enable a comparison of the present data with others [6, 7, 12]. The authors also considered the implications for the emergence of anti-malarial drug resistance, based on the data of this study and data of the other populations that were previously reported from other countries [6, 7, 12].

Methods

Collection sites and epidemiological data

The study areas, sample sizes (n), and years of the patient blood collection were as follows: Kalinga province (Northern part of Luzon Island), n = 22, collected by field surveys in 2003 through 2005, Palawan province (Palawan Island), n = 40, collected by field surveys in 2003 through 2006, and at Davao Regional Hospital in Tagum City, Davao del Norte (Southeastern part of Mindanao Island), n = 30, collected in 1999 through 2001 (Figure 1). These three study areas are sentinel provinces in the Philippines, where malaria is endemic perennially, but do not represent all the P. falciparum endemic foci of the country. The field isolates were collected in Tabuk, Kalinga province, where slide positivity rate (SPR) and annual parasite incidence (API; per 1,000 people/year) were 4.3 and 3.6 respectively (2004). Other isolates were collected in Puerto Princesa City, Palawan province where SPR and API were 18.5 and 16.3 respectively (2004). The other isolates were collected in Davao del Norte province, which had a 4.5–15.2, SPR and 2.3–7.8, API (1999–2001). APIs were collected from Annual Reports of the Provincial Health Offices of Kalinga, Palawan and Davao del Norte. Based on this epidemiological information, Kalinga can be classified as hypoendemic, Palawan as mesoendemic, and Davao del Norte as hypoendemic to mesoendemic. Information pertaining to entomological inoculation rates was not available to this study.

Figure 1
figure 1

Study areas in the Philippines. Levels of genetic differentiation (FST values) for each pairwise population comparison are shown with numbers and arrows. All values are significantly different from 0 (* P < 0.001).

Informed consent was obtained from all the patients prior to the blood collection. This study was given ethical approval by the University of the Philippines, Manila, and performed according to the ethical guidelines for epidemiological studies provided by the Ministry of Education, Culture, Sports, Science and Technology and the Ministry of Health, Labour and Welfare of Japan.

DNA extraction, PCR amplification

The patient blood samples were frozen in liquid nitrogen or preserved on filter papers (IsoCode Stix, Schleicher & Schuell co., Germany) and kept at room temperature until examined. The parasite DNA was extracted from frozen whole blood samples by phenol-chloroform extraction after proteinase K digestion or from the blood-spot samples on the filter papers, as per manufacturer's instructions [13].

Genotyping

Ten microsatellite DNA loci were amplified by semi-nested PCR. The loci were as follows: TA1 (chromosome 6), TA40 (chromosome 10), Poly α (chromosome 4), Pfg377 (chromosome 12), PfPK2 (chromosome 12), TA109 (chromosome 6), TA87 (chromosome 6), TA81 (chromosome 5), TA42 (chromosome 5) and 2490 (chromosome 10). The PCR primer sets and amplification conditions were consistent with the protocol of Anderson et al using a modified TA40 primer set [14, 15]. Sizes of fluorescence-labelled PCR products were measured on an Applied Biosystems Prism Genetic Analyzer 310 using Gene Scan version 3.1.2 with a 500 ROX size standard (ABI, CA, USA).

Different-sized PCR products amplified using the same primer set were considered to be individual alleles within a locus, as size variation among isolates is consistent with the repeat number in a microsatellite locus [14]. The electropherogram shows peak profiles for the microsatellite loci, based on fluorescence intensity of the labelled PCR products in this analysis. Multiple alleles per locus were scored if minor peaks were taller than at least one-third the height of the predominant allele for each locus. Multiple-genotype infections (MGIs) were defined as those in which at least one of the 10 loci contained more than one allele [6].

Data analysis

Expected heterozygosity (H) was calculated for each locus based on the allele frequencies of the 10 examined microsatellite loci. H values were calculated using H = [n/(n - 1)] [1 - ∑pi2], where n corresponds to the number of isolates examined and pi is the frequency of the i th allele.

Effective population size (Ne) was estimated based on H and the microsatellite mutation rate (μ = 1.59 × 10-4; 95% confidence interval: 6.98 × 10-5, 3.7 × 10-4) for P. falciparum [16–18]. The infinite-alleles model (IAM) and stepwise mutation model (SMM) were used to estimate Ne. For the IAM, the formula Neμ = H/4(1 - H) was used, whereas for the SMM, the formula Neμ = 1/8{[1/(1 - H)]2 - 1} was used [17, 18].

Each population was examined for evidence of a recent genetic bottleneck (ie, a severe decrease in population size). Heterozygosity excess was used for evidence of genetic bottlenecks [19, 20]. In a population that underwent a bottleneck, heterozygosity excess is observed, in which the value of He (the observed Hardy-Weinberg equilibrium heterozygosity) is greater than that of H (expected heterozygosity based on the number of alleles and the sample size). The mode shift in allele frequency distribution for the presence of rare alleles was also examined. A normal L-shaped distribution indicates evidence of a non-bottlenecked population, whereas a shifted mode of distribution indicates evidence of a bottlenecked population. The BOTTLENECK program, version 1.2.02, was used to search for evidence of heterozygosity excess and mode-shift [19, 21]. The Sign Test and the Wilcoxon Signed-Rank Test were used to evaluate statistical significance.

Multilocus linkage disequilibrium was assessed using the standardized index of association (IAS) [22, 23]. This analysis was performed using the LIAN 3.5 Web interface [24]. IAS was calculated using the formula IAS = (VD/Ve - 1)/(l - 1) with permutation testing of the null hypothesis of complete linkage equilibrium (IAS = 0), where VD is the observed mismatch variance, Ve is the expected mismatch variance, and l is the number of examined loci. Significances of the observed IAS values were calculated by Monte-Carlo simulation, using 10,000 random permutations of the data. This statistic is a variation of the method proposed by Maynard-Smith et al [22]. The results were standardized by the number of loci, to enable a comparison of different data sets [6, 22]. This test was applied to the data sets from each population in two ways. First, the mixed-clone infections were excluded so that only the single-clone infections were analysed, giving absolute confidence in the haplotype profile. Second, any multilocus genotype found in more than one isolate was only counted once in the analysis, ie unique haplotypes only, reducing the sample size slightly and thereby removing the possible effect of recent epidemic expansion of particular clones [6]. Since the LIAN cannot analyse isolates possessing missing data, the 57 isolates obtaining the allelic data of 10 loci were used.

The extent of population subdivision between the populations in the Philippines was estimated using Weir and Cockerham's Ø estimator for determining F statistics (FST) [25]. FST were calculated using the FSTAT program version 2.9.3.2 and tested for significant difference from 0, based on 1,000 random permutations of the data set [26].

Results

Allele frequencies of each locus in the three populations are shown in Additional File 1. Of the total of 92 isolates examined, allelic data of the microsatellite loci were obtained from the 74 isolates (80.4%)(16 from Kalinga, 28 from Palawan and 30 from Davao del Norte), while in the remaining 18 isolates (19.6%), the allelic data were not available (6 from Kalinga and 12 from Palawan) because of failures in acquiring PCR products during the genotyping. Multiple-genotype infection (MGI) was observed in some loci in three isolates from the Palawan population, and two isolates from the Davao del Norte population. No MGI was observed from the Kalinga population. Data from loci that originated from MGI were excluded from the analyses.

Genetic diversity

The genetic diversity of each population was assessed by determining the number of alleles per locus in each population and by calculating the expected heterozygosity (H) (Table 1). The mean numbers of alleles ± SE of the Kalinga, Palawan and Davao del Norte populations were 2.40 ± 0.37, 4.80 ± 0.83 and 4.50 ± 0.87, respectively. The mean of H ± SE of the Kalinga, Palawan and Davao del Norte populations were 0.39 ± 0.10, 0.60 ± 0.09 and 0.51 ± 0.08, respectively.

Table 1 Mean numbers of alleles and expected heterozygosity (H) in the three P. falciparum populations

Effective population size

Ne values of the three populations were calculated from the mean expected heterozygosity and mutation rates of P. falciparum microsatellite loci using the IAM and the SMM (Table 2). The sizes of the Kalinga, Palawan and Davao del Norte populations were 997, 2388 and 1650, respectively, based on the IAM, and 1313, 4202 and 2515 respectively, based on the SMM.

Table 2 Effective population sizes (Ne) of the three P. falciparum populations

Genetic bottleneck

Evidence of genetic bottlenecks was assessed based on heterozygosity (H) excess and patterns of allele frequency distribution (ie, mode-shift) [19, 20]. Table 3 shows the number of loci that correspond with H excess and deficiency. Statistically high levels of H excess compared with H deficiency (P < 0.05) were observed in the Kalinga and the Palawan populations when the IAM was applied to the analyses, indicating evidence of genetic bottleneck events. In contrast, there was no significant difference between the levels of H excess and H deficiency in the Davao del Norte population when the IAM and the SMM were applied to the analyses, indicating the absence of a genetic bottleneck event.

Table 3 Observed versus expected heterozygosity in the three P. falciparum populations

In the Davao del Norte population, the mode-shift indicator test showed a normal L-shaped distribution, indicating evidence of a non-bottlenecked population. This test could not be applied to the other two populations due to small sample sizes (<30) [21].

Multilocus linkage disequilibrium

IAS values were calculated for the three populations, with permutation testing of the null hypothesis of IAS = 0 (equilibrium of multilocus frequencies) (Table 4). IAS values of the three populations were highly variable. When the single-clone haplotype was used in the analysis, the IAS values ranged from 0.043 to 0.104, whereas when the unique haplotypes were used in the analysis, the IAS values ranged from 0.003 to 0.029. Significant linkage disequilibrium was observed in all three populations (P < 0.05) when the single-clone haplotype was examined using the 10 loci in the analyses. However, no significant linkage disequilibrium was found in the Kalinga and Palawan populations when the unique haplotypes were examined in the analyses.

Table 4 Multilocus linkage disequilibrium in the three P. falciparum Populations

Genetic differentiation

Levels of genetic differentiation between each pair of the three different populations were indicated by FST values using the Weir and Cockerham estimator [25, 26]. The values for each of the pairwise population comparisons ranged from 0.096 to 0.144. All values were significantly different from 0 (P < 0.001) (Figure 1).

Discussion

The failure to obtain allelic data of the 18 isolates might be due to low parasite density and/or misdiagnosis of other malaria species such as P. vivax by microscopic observations in the field during the surveys in Kalinga and Palawan. The samples of Davao del Norte were collected at a hospital where more accurate microscopic diagnosis was feasible. A possibility that differences in sample size between the Kalinga population and the other two populations might have affected the results of the analyses cannot be ruled out, and this could be the limitation of the study. However, some of the data analyses such as expected heterozygosity (H), effective population size (Ne) and multilocus linkage disequilibrium (IAS), were corrected for sample size, which enabled an appropriate comparison of different data sets.

Levels of genetic diversity of the three P. falciparum populations in the Philippines were different from each other. The Kalinga population showed a relatively low number of alleles and low levels of expected heterozygosity (H), being similar to those of South American populations of P. falciparum (H: 0.30–0.40) [6, 12]. In contrast, the Palawan and the Davao del Norte populations showed intermediate levels of H, being similar to those of mainland Southeast Asian (H: 0.51), Papua New Guinean (0.62–0.65), and Malaysian Borneo (0.46–0.63) populations of P. falciparum [6, 7].

The levels of genetic diversity seemed to be correlated with the levels of malaria endemicity. The API in Kalinga province was lower than that of the other 2 provinces. According to the report from the Center for Health Development-CAR in Kalinga province, the mean of API ± SD was 4.8 ± 1.9 (1996–2005). According to the report from the Kilusan Ligtas Malaria in Palawan province, the mean of API ± SD was 15.6 ± 4.4 (1999–2005). This difference of malaria endemicity between the two provinces probably affected the levels of genetic diversity of the two populations of P. falciparum.

Levels of genetic differentiations (FST) between the populations of P. falciparum in the Philippines (0.096–0.144) were moderate, and were similar to those between the mainland Southeast Asian and Papua New Guinean populations of P. falciparum (0.121–0.140), and lower than those between Sabah and Sarawak populations of P. falciparum (0.217–376) in Malaysian Borneo, as previously reported [6, 7].

In the present study, possibilities of genetic bottlenecks with the Kalinga and the Palawan populations of P. falciparum were found. However, the significant heterozygosity excess (evidence of genetic bottlenecks) was observed only when the IAM was used for the analysis. This result is not robust because modes of microsatellite mutation theoretically fit the SMM. Luikart and Cornuet recommended that only the SMM should be used for microsatellite data to test for genetic bottlenecks, because the IAM may detect heterozygosity excess even in non-bottlenecked populations [20].

Significant multilocus linkage disequilibrium was obtained from all three populations of P. falciparum in the Philippines, when a single clone was used in the analysis [22–24]. The Kalinga and the Palawan populations showed decreased linkage disequilibrium with no levels of significance when the unique haplotypes were included in the data set. This is because the identical haplotype was frequently found from two or more isolates (patients) within the same population in Kalinga and Palawan. In particular, in the Kalinga population, 57% of the haplotypes (four out of seven haplotypes) were shared among two or more isolates (patients). In fact, such isolates were collected from the same site and on the same day. Maynard-Smith et al proposed simple methods to distinguish between "clonal" population structure and "epidemic" population structure of microbial pathogens [22]. They demonstrated that epidemic population structure could be identified by treating the multiply-represented haplotype as single individuals and then re-examining the linkage disequilibrium. If a decreased level of linkage disequilibrium is observed in the second analysis, the population is estimated to have an epidemic population structure. The results indicated that the Kalinga and Palawan populations had an epidemic population structure.

On the contrary, the Davao del Norte population persisted with significant linkage disequilibrium (P < 0.05) even when the multiply-represented haplotype was treated as single individuals. This result indicated that a recombination rate was estimated as low in the Davao del Norte population, and that the population has possibly not experienced the recent epidemic expansion of a particular haplotype.

Generally, levels of malaria transmission have a negative correlation with levels of linkage disequilibrium [6, 7, 12]. In areas where malaria transmission is high, levels of linkage disequilibrium are likely to decrease because one malaria patient can be infected with more than two genetically different clones of the parasite, and a mosquito can frequently acquire more than two genetically different clones in the same blood meal. Then, meiotic recombination of the parasite genomes will occur in the mosquito midgut and decrease the levels of linkage disequilibrium. In contrast, in areas where malaria transmission is low, levels of linkage disequilibrium are likely to increase because, generally, the frequency of multiple-genotype infection (MGI) and the number of clones per patient are lower than those in high transmission areas, and the chances of acquiring more than two genetically different clones with one blood feeding by a mosquito will also be lower than those in high transmission areas.

In the present study, MGI in a patient was observed in the Palawan (three of 28 isolates: 10.7%) and the Davao del Norte (two of 30 isolates: 6.7%) populations, but not in the Kalinga population. MGI is important to the maintenance of genetic diversity of the P. falciparum population. Without multiple infections, this parasite has no chance of encountering genetically different clones in a mosquito midgut. In the present study, a clear correlation between the levels of genetic diversity (H) and the frequencies of MGI was found. The H values of the three populations were elevated as follows: Kalinga (0.39) < Davao del Norte (0.51) < Palawan (0.60), and the frequencies of MGI were also elevated as follows: Kalinga (0.0%) < Davao del Norte (6.7%) < Palawan (10.7%). In contrast, a negative correlation between the levels of linkage disequilibrium and the frequencies of MGI of the parasite was found. When single clones were used in the data set for the linkage disequilibrium analysis (Table 4), the levels of linkage disequilibrium were elevated as follows: Palawan < Davao del Norte < Kalinga, while the frequencies of MGI were elevated as follows: Kalinga (0.0%) < Davao del Norte (6.7%) < Palawan (10.7%). These results are consistent with the previously reported population genetics trends of P. falciparum [6, 7, 12].

However, caution is needed for detecting MGI by microsatellite typing, because it is likely to underestimate the prevalence of MGI of P. falciparum. In Brazil, over 40% of local isolates were MGI detected by a single-copy antigen-coding gene (merozoite surface protein-1: msp-1), however less than 20% of local isolates were MGI detected by microsatellite typing [6, 27]. In fact, nested PCR genotyping of block 2 of the msp-1 gene and block 3 of msp-2 [28] in the present samples showed that more than 20% of the isolates of each of the three populations proved to be MGI (Additional File 2).

The levels of genetic differentiation between the three populations in the Philippines were moderate, and lower than those reported in Malaysian Borneo where the parasite populations were highly fragmented from each other [7]. High levels of genetic differentiation of the populations in Malaysian Borneo are considered to be a result of fragmentation of the population structure owing to an effective malaria control and low migration rate of people in the endemic areas. Therefore it is expected that an elevation of FST values between the three populations in the Philippines will be detected in the near future, with the onset of effective malaria control.

However, the migration rate of people to the endemic areas in the Philippines may be different from that in Malaysian Borneo. In Malaysian Borneo, malaria endemic areas are typically in remote, isolated mountainous villages. In contrast, in the Philippines, malaria is not only endemic in remote areas (e.g. forests, mountainous villages) but also in relatively big cities or ports (e.g. Puerto Princesa City, Palawan province). For example, in Palawan, malaria is highly endemic in villages near seashores where fishermen live, and they sail from island to island, freely catching and selling fish.

Although mechanisms for the emergence of drug resistance have not been fully understood, levels of genetic diversity and characteristics of population structure of P. falciparum would be the key factors to estimate a potential of the emergence of drug resistance. Ariey et al demonstrated that a metapopulation concept would be applicable to explain the emergence of drug resistant malaria [29, 30]. They showed that the P. falciparum populations in Southeast Asia and South America, where chloroquine resistance emerged, were consistent with a typical metapopulation structure, ie, parasite populations were differentiated in several subpopulations but not completely isolated from each other. The levels of genetic diversity and the characteristics of population structure of P. falciparum in the Philippines were similar to those reported in mainland Southeast Asia or South America where chloroquine resistant malaria emerged [6]. Indeed, Chen et al indicated that some chloroquine resistant mutations seemed to have emerged in the Philippines [10]. Moreover, a chloroquine resistant malaria determined by an in vitro chloroquine susceptibility test has been reported from Kalinga, Palawan and Davao del Norte [8, 9, 31]. The results in this study suggest that the populations in the Philippines have the potential to make the drug-resistant mutation(s) spread into different endemic areas in this country.

Conclusion

The genetic diversity and population structure of P. falciparum in the Philippines were examined using 10 highly polymorphic microsatellite loci of the parasite. The levels of the genetic diversity and the frequencies of MGI correlated with the local malaria endemicity. On the other hand, the levels of linkage disequilibrium correlated negatively with the local malaria endemicity. Genetic population analyses of P. falciparum using the microsatellite loci provide insight not only into understanding the basic biology of this organism, but also the epidemiology and control of the disease.

References

  1. Tongol-Rivera P: Milestones in the history of malaria research and its control in the Philippines. Malaria in Asia. Edited by: Kano S, Tongol-Rivera P. 2005, The Federation of Asian Parasitologists, 135-166.

    Google Scholar 

  2. Tongol-Rivera P, Villacorte EA, Escueta AS, Hatabu T, Iwagami M, Kawazu SI, Kano S: A study of chloroquine resistance of Plasmodium falciparum using the In-vitro sensitivity test and polymerase chain reaction (PCR). Acta Medica Philippina. 2005, 39: 7-10.

    Google Scholar 

  3. World Health Organization: Philippines country profile. World Malaria Report. 2005, Geneva, Switzerland: WHO/UNICEF/RBM, [http://rbm.who.int/wmr2005/profiles/philippines.pdf]

    Google Scholar 

  4. 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 

  5. Zilversmit M, Hartl DL: Evolutionary history and population genetics of human malaria parasites. Molecular approaches to malaria. Edited by: Sherman IW. 2005, Washington, DC: American Society for Microbiology Press, 95-109.

    Chapter  Google Scholar 

  6. 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 

  7. Anthony TG, Conway DJ, Cox-Singh J, Matusop A, Ratnam S, Shamsul S, Singh B: Fragmented population structure of Plasmodium falciparum in a region of declining endemicity. J Infect Dis. 2005, 191: 1558-1564. 10.1086/429338.

    Article  PubMed  Google Scholar 

  8. Hatabu T, Kawazu S, Suzuki J, Valenzuela RF, Villacorte EA, Suzuki M, Rivera PT, Kano S: In vitro susceptibility of Plasmodium falciparum isolates to chloroquine and mefloquine in southeastern Mindanao Island, the Philippines. Southeast Asian J Trop Med Public Health. 2003, 34: 546-551.

    CAS  PubMed  Google Scholar 

  9. Iwagami M, Taguchi N, Angluben RU, Escueta AS, Villacorte EA, Rivera PT, Kano S, Kawazu SI: Evaluation of blood preservation methods in the performance of the WHO in vitro micro-test for Plasmodium falciparum in the field. Tropical Medicine and Health. 2007, 35: 337-341. 10.2149/tmh.35.337.

    Article  Google Scholar 

  10. Chen N, Wilson DW, Pasay C, Bell D, Martin LB, Kyle D, Cheng Q: Origin and dissemination of chloroquine-resistant Plasmodium falciparum with mutant pfcrt alleles in the Philippines. Antimicrob Agents Chemother. 2005, 49: 2102-2105. 10.1128/AAC.49.5.2102-2105.2005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Sakihama N, Nakamura M, Palanca AA, Argubano RA, Realon EP, Larracas AL, Espina RL, Tanabe K: Allelic diversity in the merozoite surface protein 1 gene of Plasmodium falciparum on Palawan Island, the Philippines. Parasitol Int. 2007, 56: 185-194. 10.1016/j.parint.2007.01.011.

    Article  CAS  PubMed  Google Scholar 

  12. Machado RL, Povoa MM, Calvosa VS, Ferreira MU, Rossit AR, dos Santos EJ, Conway DJ: Genetic structure of Plasmodium falciparum populations in the Brazilian Amazon region. J Infect Dis. 2004, 190: 1547-1555. 10.1086/424601.

    Article  CAS  PubMed  Google Scholar 

  13. Sambrook J, Russell DW: Molecular Cloning: A Laboratory Manual. 2001, New York: Cold Harbor Laboratory Press, 3

    Google Scholar 

  14. 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: 113-125. 10.1017/S0031182099004552.

    Article  CAS  PubMed  Google Scholar 

  15. Su XZ, 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 

  16. Su X, Ferdig MT, Huang Y, Huynh CQ, Liu A, You J, Wootton JC, Wellems TE: A genetic map and recombination parameters of the human malaria parasite Plasmodium falciparum. Science. 1999, 286: 1351-1253. 10.1126/science.286.5443.1351.

    Article  CAS  PubMed  Google Scholar 

  17. Schug MD, Mackay TF, Aquadro CF: Low mutation rates of microsatellite loci in Drosophila melanogaster. Nat Genet. 1997, 15: 99-102. 10.1038/ng0197-99.

    Article  CAS  PubMed  Google Scholar 

  18. Nei M, Kumar S: Molecular Evolution and Phylogenetics. 2000, New York: Oxford University Press

    Google Scholar 

  19. Cornuet JM, Luikart G: Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics. 1996, 144: 2001-2014.

    PubMed Central  CAS  PubMed  Google Scholar 

  20. Luikart G, Cornuet JM: Empirical Evaluation of a test for identifying recently bottlenecked populations from allele frequency data. Conservation Biology. 1998, 12: 228-237. 10.1046/j.1523-1739.1998.96388.x. [http://www3.interscience.wiley.com/cgi-bin/fulltext/120829627/HTMLSTART]

    Article  Google Scholar 

  21. Piry S, Luikart G, Cornuet JM: BOTTLENECK: a computer program for detecting recent reductions in the effective population size using allele frequency data. J Hered. 1999, 90: 502-503. 10.1093/jhered/90.4.502. [http://www.montpellier.inra.fr/URLB/bottleneck/bottleneck.html]

    Article  Google Scholar 

  22. Maynard-Smith J, Smith NH, O'Rourke M, Spratt BG: How clonal are bacteria?. Proc Natl Acad Sci USA. 1993, 90: 4384-4388. 10.1073/pnas.90.10.4384.

    Article  Google Scholar 

  23. Hudson RR: Analytical results concerning linkage disequilibrium in models with genetic transformation and recombination. J Evol Biol. 1994, 7: 535-548. 10.1046/j.1420-9101.1994.7050535.x.

    Article  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  Google Scholar 

  26. Goudet J: FSTAT (version 1.2): a computer program to calculate F statistics. J Hered. 1995, 86: 485-486. [http://www2.unil.ch/popgen/softwares/fstat.htm]

    Google Scholar 

  27. Ferreira MU, Karunaweera ND, da Silva-Nunes M, da Silva NS, Wirth DF, Hartl DL: Population structure and transmission dynamics of Plasmodium vivax in rural Amazonia. J Infect Dis. 2007, 195: 1218-1226. 10.1086/512685.

    Article  PubMed  Google Scholar 

  28. Snounou G, Zhu X, Siripoon N, Jarra W, Thaithong S, Brown KN, Viriyakosol S: Biased distribution of msp1 and msp2 allelic variants in Plasmodium falciparum populations in Thailand. Trans R Soc Trop Med Hyg. 1999, 93: 369-374. 10.1016/S0035-9203(99)90120-7. Erratum in: Snounou G: Trans R Soc Trop Med Hyg 2000, 94:65

    Article  CAS  PubMed  Google Scholar 

  29. Ariey F, Duchemin JB, Robert V: Metapopulation concepts applied to falciparum malaria and their impacts on the emergence and spread of chloroquine resistance. Infect Genet Evol. 2003, 2: 185-192. 10.1016/S1567-1348(02)00099-0.

    Article  CAS  PubMed  Google Scholar 

  30. Hanski I, Simberloff D: The metapopulation approach, its history, conceptual domain, and application to conservation. Metapopulation Biology, Ecology, Genetics, and Evolution. Edited by: Hanski I, Gilpin ME. 1997, San Diego: Academic Press

    Google Scholar 

  31. Hatabu H, Iwagami M, Kawazu S, Taguchi N, Escueta DA, Villacorte EA, Tongol-Rivera P, Kano S: Association of molecular markers in Plasmodium falciparum crt and mdr1 with in vitro chloroquine resistance: A Philippines study. Parasitol Int. 2009, 58: 166-170. 10.1016/j.parint.2009.01.010.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Authors would like to thank the following people for kind assistance in obtaining samples and epidemiological data: Dr. Nao Taguchi, Gunma University School of Health Sciences, Gunma, Japan: Mr. Ray U. Angluben, Mr. Darius de la Cruz, Mr. Conrado E. Mendoza Jr, Kilusan Ligtas Malaria, Puerto Princesa City, Palawan province, the Philippines: Dr. Gavino, Provincial Health Officer, Kalinga province, the Philippines: Ms. Norma B. Bunaton, Ms. Linda Esquivel, Center for Health Development-CAR, Tabuk City, Kalinga province, the Philippines: the Dr. Romulo Busuego, Director of Davao Regional Hospital, Tagum City, Dr. Alfredo Lacerona, Municipal Health Officer, Kapalong, Davao del Norte province, the Philippines. The authors also thank Pilipinas Shell Foundation Inc., Palawan, for supporting the field survey in Palawan province. This work was supported by a Grant-in-Aid for Scientific Research (B) (16406012, 19406013) from the Ministry of Education, Culture, Sports, Science and Technology of Japan.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Shigeyuki Kano.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MI carried out the molecular genetic studies, performed the population genetic analyses and drafted the manuscript. PTR, EAV, ADE, THat and SIK collected the patients' blood samples as well as helping with the writing of the manuscript. THay and KT helped with the population genetic analyses and helped with the writing of the manuscript. SK participated in the design of the study, acquisition of funding, coordination and writing of the manuscript.

Electronic supplementary material

12936_2009_845_MOESM1_ESM.xls

Additional File 1: Allele frequencies and number of isolates (n) of 10 microsatellite loci of the three Plasmodium falciparum populations in the Philippines. Allele names refer to the base pair sizes of the product of polymerase chain reaction. (XLS 46 KB)

12936_2009_845_MOESM2_ESM.xls

Additional File 2: Percentages of multiple-genotype infection (MGI) of msp1 and msp2 of the three Plasmodium falciparum populations in the Philippines. No. indicates the number of isolates for each measure.msp1: block 2 of merozoite surface protein-1 gene, msp2: block 3 of merozoite surface protein-2 gene, msp1 or msp2 represents the percentage of isolates that showed MGI in either msp1 or msp2. (XLS 42 KB)

Authors’ original submitted files for images

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

Authors’ original file for figure 1

Authors’ original file for figure 2

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

Iwagami, M., Rivera, P.T., Villacorte, E.A. et al. Genetic diversity and population structure of Plasmodium falciparum in the Philippines. Malar J 8, 96 (2009). https://doi.org/10.1186/1475-2875-8-96

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-2875-8-96

Keywords