Differences in genetic population structures of Plasmodium falciparum isolates from patients along Thai-Myanmar border with severe or uncomplicated malaria
- Pannapa Susomboon1, 2,
- Moritoshi Iwagami1,
- Noppadon Tangpukdee2,
- Srivicha Krusood2,
- Sornchai Looareesuwan^2 and
- Shigeyuki Kano1Email author
© Susomboon et al; licensee BioMed Central Ltd. 2008
Received: 24 April 2008
Accepted: 21 October 2008
Published: 21 October 2008
There have been many reports on the population genetic structures of Plasmodium falciparum from different endemic regions, but few studies have examined the characteristics of isolates from patients with different clinical outcomes. The population genetic structures of P. falciparum isolates from patients with either severe or uncomplicated malaria were examined.
Twelve microsatellite DNA loci from P. falciparum were used to assess the population genetic structures of 50 isolates (i.e., 25 isolates from patients with severe malaria and 25 from patients with uncomplicated malaria) collected in the Thai-Myanmar border area between 2002 and 2005.
Genetic diversity and effective population sizes were greater in the uncomplicated malaria group than in the severe malaria group. Evidence of genetic bottlenecks was not observed in either group. Strong linkage disequilibrium was observed in the uncomplicated malaria group. The groups demonstrated significant genetic differentiation (P < 0.05), and allele frequencies for 3 of the 12 microsatellite loci differed significantly between the two groups.
These findings suggest that the genetic structure of P. falciparum populations in patients with severe malaria differs from that in patients with uncomplicated malaria. The microsatellite loci used in this study were presumably unrelated to antigenic features of the parasites, but, these findings suggest that some loci may influence the clinical outcome of malaria.
Plasmodium falciparum is the most pathogenic of the protozoan parasites that cause human malaria. According to the World Health Organization, more than 1 million people die from malaria each year . The clinical manifestations of malaria are quite pleomorphic, ranging from mild or asymptomatic parasitemia to potentially fatal conditions such as coma and multi-organ failure. Although the molecular basis of severe malaria has been well studied in recent years , determinants of the clinical outcomes of malaria remain unknown. Several factors, including host and parasite genetic characteristics, are thought to contribute to the clinical outcome of malaria.
Several studies have reported evidence of a relationship between parasite genotype and clinical outcome [2–8]. Ariey et al found that certain alleles in the polymorphic microsatellite loci of clinical isolates from French Guyana were associated with disease severity ; however, Ferreira et al found no evidence to suggest that the parasite genotypes of the isolates from malaria patients in Vietnam were associated with disease severity . These results demonstrate the difficulties associated with evaluating the role of parasite genetic factors in malaria pathogenesis.
The present study examines the relationship between parasite genetic factors and clinical outcomes. To this end, the genetic population structures of parasites isolated from patients with either uncomplicated or severe malaria were compared using highly polymorphic microsatellite DNA loci from P. falciparum.
A total of 50 P. falciparum isolates were collected from patients who contracted malaria along the Thai-Myanmar border and who had been admitted to the Bangkok Hospital for Tropical Diseases, Faculty of Tropical Medicine, Mahidol University between 2002 and 2005. Inclusion criteria of the study are: age ≥ 15 years old, both male and female, positive for P. falciparum malaria infection with microscopic confirmation showing acute manifestation, informed consent to take part in the study obtained either from the patient or the legal guardian, and signs of severe malaria according to the World Health Organization criteria  (for patients in severe malaria group). Exclusion criteria are: pregnant or lactating women, pre-treatment with any anti-malarial drugs within the past two weeks, evidence of mixed infection on admission, and the evidence of any contra-indication of blood sampling collection. Consequently, the 50 isolates were divided into two clinical groups (25 isolates from patients with severe malaria and 25 from patients with uncomplicated malaria) (Additional File 1, 2). Hyperparasitemia (> 250,000 parasite/μL), peripheral schizontemia, blood urea nitrogen (> 60 mg/dL) and elevation of serum enzymes (AST, ALT) were noted for patients in the severe group. Clinical manifestations such as nausea/vomiting and jaundice were also taken into account of the clinical severity of the patients.
The geographic and temporal distributions of the presumed area of infection were similar among both groups. Blood samples and clinical data were collected from patients after obtaining informed consent. Blood samples were frozen at -80°C until use.
This study was approved by the Ethics Committee of the Faculty of Tropical Medicine at Mahidol University in Thailand, and complied with the ethical guidelines for epidemiological studies set forth by the Japanese Ministry of Education, Culture, Sports, Science, and Technology and the Japanese Ministry of Health, Labour, and Welfare.
Preparation of parasite genomic DNA and genotyping analyses
Genomic DNA was extracted from the parasites in 2 mL of frozen blood sample using previously described methods .
Twelve microsatellite DNA loci were amplified by semi-nested PCR. The loci were as follows: TA1 (Chromosome 6), TA40 (Chromosome 10), Poly α (Chromosome 4), TA60 (Chromosome 13), ARAII (Chromosome 11), 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 [11, 12]. Sizes of fluorescence-labeled 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.
Expected heterozygosity (H) was calculated for each locus based on the allele frequencies of the 12 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 [13–15]. The infinite-alleles model (IAM) and the stepwise mutation model (SMM) were used to estimate Ne.
Each population was examined for evidence of a recent genetic bottleneck (i.e., a severe decrease in population size). In non-bottlenecked populations approaching mutation-drift equilibrium, the H value for each locus was calculated based on the number of alleles, and the sample size was equal to the observed Hardy-Weinberg equilibrium heterozygosity (He) [16, 17]. After a bottleneck event, the number of alleles and the expected heterozygosity of a population are predicted to decrease; however, rare alleles are purged more quickly, reducing the number of alleles without altering the extent of heterozygosity [16, 17]. Furthermore, the mode shift in allele frequency distribution for the presence of rare alleles was examined. The BOTTLENECK program (version 1.2.02)  was used to search for evidence of heterozygosity excess and mode-shift. 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) [19, 20]. This analysis was performed using the LIAN 3.5 Web interface . 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. To enable comparison of different data sets, the results were standardized by the number of loci [19, 22].
The extent of population subdivision was estimated using Weir and Cockerham's theta estimator for determining F statistics (F ST ) . F ST were calculated using the FSTAT program (version 22.214.171.124: available at http://www2.unil.ch/popgen/softwares/fstat.htm)  and tested for significant difference from 0, based on 1,000 random permutations of the data set.
Fisher's exact test was used to calculate differences in allele frequencies for each locus in the uncomplicated and severe malaria groups (SAS for Windows, version 9.1.3; SAS Institute Inc., Cary, NC, USA). A p-value of less than 0.05 was considered statistically significant.
Allele frequencies for each locus are shown in Additional File 3. Of the 50 isolates examined, 45 (i.e., 90%) represented single-genotype infections involving all 12 loci. The remaining five isolates (i.e., 10%) represented multiple-genotype infections involving one to four loci. Three of the five isolates were obtained from the severe malaria group, while the remaining two isolates were obtained from the uncomplicated malaria group. Data from multiple loci were excluded from the analyses. We excluded data from loci that originated from multi-genotype infection.
Sample size, mean number of alleles, and expected heterozygosity in the two study populations.
No. of isolates
No. of alleles
Expected heterozygosity (H)
mean ± SE
mean ± SE
mean ± SE
24.50 ± 0.15
5.17 ± 0.66
0.60 ± 0.07
24.67 ± 0.19
6.58 ± 0.66
0.72 ± 0.05
Effective population size
Effective size (Ne) values for the two study populations.
2,356 (1,012, 5,366)
4,120 (1,770, 9,385)
3,944 (1,695, 8984)
8,890 (3,820, 20,251)
Observed versus expected heteroxygosity in the study populations.
No. of loci
P < 0.05
Multilocus linkage disequilibrium
Multilocus linkage disequilibrium in the two study populations.
No. of loci
No. of isolates
I A S
Genetic differentiation and allele frequency-distribution patterns
The genetic structures of two parasite populations collected from patients with either severe or uncomplicated malaria who contracted malaria, from the Thai-Myanmar border region were assessed. This area is known to be endemic for multidrug-resistant malaria [25–27]; transmission is unstable and hypoendemic . However, morbidity and mortality in the area have gradually decreased in recent years, with a 2006 parasite incidence of 0.49 per 1,000 people .
In the present study, only 10% of the isolates were obtained from multi-genotype infections involving one to four loci. This is significantly lower than a previous estimate (i.e., 40%) reported for the same region in 1997 and 1998 . The frequency of multiple-genotype infection by malaria parasites is known to increase as malaria transmission becomes more prevalent. For example, more than 50% of isolates are obtained from multiple-genotype infection in high-transmission areas (e.g., Uganda, Congo, and Zimbabwe), while fewer than 10% are obtained from multiple-genotype infection in low-transmission areas (e.g., Colombia and Brazil) . Therefore, the present study indicates that malaria transmission has decreased in the Thai-Myanmar border region in recent years.
The genetic structure of the P. falciparum population in the severe and uncomplicated malaria groups differed in regards to several genetic indexes. Specifically, H values in the severe population were lower than those in the uncomplicated population. H values determined for both the severe and uncomplicated groups were higher than previously reported values for the same area . Ne values for the severe population were also lower than those in the uncomplicated population, and our values were generally higher than previously reported values for this area . Discrepancies between our study and previous reports may have resulted from differences in the length of sampling periods. In the present study, samples were collected between 2002 and 2005, while other investigators collected samples between December 1997 and January 1998. Thus, the longer sampling period used in our study may have resulted in higher H and N e values in the P. falciparum populations.
Significant genetic differentiation was observed between the severe and uncomplicated malaria groups (p < 0.05), even though the groups were collected from the same endemic area. This may reflect differences in transmission intensity between the two groups.
Allele frequency distributions for 3 (i.e., Pfg377, TA109, and TA42) of 12 loci differed significantly between the uncomplicated and severe malaria groups (P < 0.05). Pfg377 is located within the telomeric region of chromosome 12, while TA109 is located within a hypothetical protein on the centromeric region of chromosome 6, and TA42 is located within a hypothetical protein on the telomeric region of chromosome 5.
Differences in allele frequencies between the severe and uncomplicated malaria groups are mainly attributed to variations in minor alleles among the respective groups. The sizes and frequencies of the major alleles were similar in the two groups, while the numbers of minor alleles in the uncomplicated group were higher than those in the severe group.
In a previous study of the Western Amazon, Martha et al found that allele frequency distribution within the TA42 locus differs between symptomatic and non-symptomatic malaria patients . In light of these results, the authors suggested that specific alleles within the TA42 locus are dominant in non-symptomatic malaria patients; however, evidence of such phenomena was not observed in either the severe or uncomplicated populations.
Other possibilities that are not linked with the parasite genetic background and that could influence the population structure, such as age of the patient, or number of days between onset of the symptom and admission to the hospital, were excluded by the statistical analysis of the two groups (Additional File 1, 2).
Many population geneticists have accepted a model of P. falciparum infection that depicts low-transmission regions as having high levels of monoclonal infection, inbreeding, and rare recombination; while high-transmission regions are depicted as having frequent mosquito inoculation, multiple-genotype infections, frequent outbreeding, and extensive recombination . According to this model, genetic diversity increases in high-transmission regions and decreases in low-transmission regions. The results described in the present study appear to link genetic diversity with transmission levels in both the severe and uncomplicated populations. Luxemburger et al found that only 5% of malaria cases in the Thai-Myanmar border region in 1992 were severe: the remainder was uncomplicated .
The genetic indexes examined in this study were generally consistent with transmission levels; however, the level of linkage disequilibrium in the severe malaria group was lower than that in the uncomplicated malaria group. For reasons that remain unclear, significant linkage disequilibrium was evident in the uncomplicated malaria group but not in the severe group. Further studies are required to explore this finding in greater detail.
Twelve highly polymorphic microsatellite loci were examined and the genetic structure of P. falciparum populations in patients with severe malaria was demonstrated to differ from that in patients with uncomplicated malaria. The microsatellite loci used in this study were presumably unrelated to the antigenic features of the parasites; however, it was suggested that the loci might somehow influence the clinical outcome of malaria.
The authors would like to thank the malaria patients for their willingness to participate in this study. We gratefully acknowledge the staff of Bangkok Hospital for Tropical Diseases, Faculty of Tropical Medicine, Mahidol University, Thailand for their help, as well as Dr. T. Mizoue from the Department of Epidemiology and International Health at the International Medical Center of Japan for assistance with statistical analyses. This study was supported by Grants-in-Aid for Scientific Research (B16406012 and B19406013) awarded by the Japanese Ministry of Education, Culture, Sports, Science, and Technology, and by Health and Labour Sciences Research Grants for Research on Emerging and Re-emerging Infectious Diseases awarded by the Japanese Ministry of Health, Labour, and Welfare (H18-SHINKOU-IPPAN-009). This study was also supported by the Tokyo Biochemical Research Foundation.
- World Health Organization: Introduction. World Malaria Report 2005. 2005, Geneva, World Health Organization, 1-3.View ArticleGoogle Scholar
- Chen Q, Schlichtherle M, Wahlgren M: Molecular aspects of severe malaria. Clin Microbiol Rev. 2000, 13: 439-450. 10.1128/CMR.13.3.439-450.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- James SP, Nicol WD, Shute PG: A study of induced malignant tertian malaria. Proc R Soc Med. 1932, 25: 1153-1181.PubMed CentralPubMedGoogle Scholar
- James SP, Nicol WD, Shute PG: Clinical and parasitological observations on induced malaria. Proc R Soc Lond. 1936, B 29: 879-894.Google Scholar
- Gupta S, Hill AVS, Kwiatkowski D, Greenwood AM, Greenwood BM, Day KP: Parasite virulence and disease patterns in Plasmodium falciparum malaria. Proc Natl Acad Sci USA. 1994, 91: 3715-3719. 10.1073/pnas.91.9.3715.PubMed CentralView ArticlePubMedGoogle Scholar
- Ariey F, Hommel D, Le Scanf C, Duchemin JB, Peneau C, Hulin A, Sarthou JL, Reynes JM, Fandeur T, Mercereau-Puijalon O: Association of severe malaria with a specific Plasmodium falciparum genotype in French Guiana. J Infect Dis. 2001, 184: 237-241. 10.1086/322012.View ArticlePubMedGoogle Scholar
- Martha RC, Tada MS, Ferreira RG, Silva LH, Wunderlich G: Microsatellite characterization of Plasmodium falciparum from symptomatic and non-symptomatic infections from the Western Amazon reveals the existence of non-symptomatic infection-associated genotypes. Mem Inst Oswaldo Cruz. 2007, 102: 293-298. 10.1590/S0074-02762007005000044.View ArticlePubMedGoogle Scholar
- Ferreira MU, Nair S, Hyunh TV, Kawamoto F, Anderson TJ: Microsatellite characterization of Plasmodium falciparum from cerebral and uncomplicated malaria patients in southern Vietnam. J Clin Microbiol. 2002, 40: 1854-1857. 10.1128/JCM.40.5.1854-1857.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- World Health Organization: Severe falciparum malaria. Trans R Soc Trop Med Hyg. 2000, 94 (Suppl 1): 1-90. 10.1016/S0035-9203(00)90300-6.View ArticleGoogle Scholar
- Sambrook J, Russell DW: Molecular Cloning: A Laboratory Manual. 2001, New York: Cold Harbor Laboratory Press, 3Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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-1353. 10.1126/science.286.5443.1351.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Nei M, Kumar S: Molecular Evolution and Phylogenetics. 2000, New York: Oxford University PressGoogle Scholar
- 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 CentralPubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Maynard-Smith J, Smith NH, O'Rourke M, Spratt BG: How clonal are bacteria?. Proc Natl Acad Sci USA. 1993, 90: 4384-43848. 10.1073/pnas.90.10.4384.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Haubold B, Hudson RR: LIAN 3.0: detecting linkage disequilibrium in multilocus data. Bioinformatics. 2000, 16: 847-848. 10.1093/bioinformatics/16.9.847.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.View ArticleGoogle Scholar
- Goudet J: FSTAT (version 1.2): a computer program to calculate F statistics. J Hered. 1995, 86: 485-486.Google Scholar
- Congpuong K, Na Bangchang K, Mungthin M, Bualombai P, Wernsdorfer WH: Molecular epidemiology of drug resistance markers of Plasmodium falciparum malaria in Thailand. Trop Med Int Health. 2005, 10: 717-722. 10.1111/j.1365-3156.2005.01450.x.View ArticlePubMedGoogle Scholar
- Nosten F, ter Kuile F, Chongsuphajaisiddhi T, Luxemburger C, Webtser HK, Edstein M, Phaipun L, Thew KL, White NJ: Mefloquine-resistant falciparum malaria on the Thai-Burmese border. Lancet. 1991, 337: 1140-1143. 10.1016/0140-6736(91)92798-7.View ArticlePubMedGoogle Scholar
- Nelson AL, Purfield A, McDaniel P, Uthaimongkol N, Buathong N, Sriwichai S, Miller RS, Wongsrichanalai C, Meshnick SR: pfmdr1 Genotyping and in vivo mefloquine resistance on the Thai-Myanmar border. Am J Trop Med Hyg. 2005, 72: 586-592.PubMedGoogle Scholar
- Luxemburger C, Thwai KL, White NJ, Webster HK, Kyle DE: The epidemiology of malaria in a Karen population on the western border of Thailand. Trans R Soc Trop Med Hyg. 1996, 90: 105-111. 10.1016/S0035-9203(96)90102-9.View ArticleGoogle Scholar
- Ministry of Public Health of Thailand 2003. [http://www.rbm.who.int/wmr2005/profiles/thailand.pdf]
- Luxemburger C, Ricci F, Nosten F, Raimond D, Bathet S, White NJ: The epidemiology of severe malaria in an area of low transmission in Thailand. Trans R Soc Trop Med Hyg. 1997, 91: 256-262. 10.1016/S0035-9203(97)90066-3.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.