- Open Access
Differences in genetic population structures of Plasmodium falciparum isolates from patients along Thai-Myanmar border with severe or uncomplicated malaria
Malaria Journalvolume 7, Article number: 212 (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 188.8.131.52: 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.
The genetic diversity of each population was assessed by determining the number of alleles per locus in each population and by calculating the expected H values (Table 1). The mean numbers of alleles ± SE in the severe and uncomplicated malaria populations were 5.17 ± 0.66 and 6.58 ± 0.66, respectively. The H values ± SE for the severe and uncomplicated malaria populations were 0.60 ± 0.07 and 0.71 ± 0.05, respectively. These results indicated that the genetic diversity in the uncomplicated population was slightly greater than that in the severe population, although the differences were not statistically significant.
Effective population size
Ne values were calculated from the mean expected heterozygosity and mutation rates of P. falciparum microsatellite loci using the infinite-allele model (IAM) and the stepwise mutation model (SMM) (Table 2) [14, 15]. The sizes of the severe and uncomplicated malaria populations were 2,356 and 3,944, respectively, based on IAM, and 4,120 and 8,890, respectively, based on SMM. These findings indicated that the severe malaria population was genetically less divergent than the uncomplicated malaria population.
Evidence of genetic bottleneck was assessed based on heterozygosity excess and patterns of allele frequency distribution (i.e., mode-shift) [16, 17]. Table 3 shows the number of loci corresponding to H excess and deficiency. Statistically high levels of H excess compared with H deficiency were not observed in either study population when IAM and SMM were applied to the analyses, indicating the absence of genetic bottleneck events. The mode-shift indicator test revealed a normal L-shaped allele frequency distribution in the two study populations (Table 3), further demonstrating the absence of genetic bottleneck events.
Multilocus linkage disequilibrium
I A Svalues were calculated for 12 loci from 22 isolates in the severe malaria group and 23 isolates in the uncomplicated malaria group. I A Svalues were equal to 0.010 and 0.127 for the severe and uncomplicated groups, respectively (Table 4). Significant linkage disequilibrium was observed in the uncomplicated malaria population (p < 0.0001).
Genetic differentiation and allele frequency-distribution patterns
Genetic differentiation between the two populations was indicated by an F ST value of 0.014. This value was found to be significantly different from 0 (p < 0.05). The distribution of alleles for each locus in the severe malaria group was compared with that in the uncomplicated group. Distribution patterns for 3 (i.e., Pfg377, TA109, and TA42) of the 12 loci were significantly different (p < 0.05) between the two populations (Figure 1). Specifically, distribution patterns in the uncomplicated malaria population were more variable than those in the severe malaria population.
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.
World Health Organization: Introduction. World Malaria Report 2005. 2005, Geneva, World Health Organization, 1-3.
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.
James SP, Nicol WD, Shute PG: A study of induced malignant tertian malaria. Proc R Soc Med. 1932, 25: 1153-1181.
James SP, Nicol WD, Shute PG: Clinical and parasitological observations on induced malaria. Proc R Soc Lond. 1936, B 29: 879-894.
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.
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.
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.
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.
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.
Sambrook J, Russell DW: Molecular Cloning: A Laboratory Manual. 2001, New York: Cold Harbor Laboratory Press, 3
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.
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.
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.
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.
Nei M, Kumar S: Molecular Evolution and Phylogenetics. 2000, New York: Oxford University Press
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.
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.
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.
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.
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.
Haubold B, Hudson RR: LIAN 3.0: detecting linkage disequilibrium in multilocus data. Bioinformatics. 2000, 16: 847-848. 10.1093/bioinformatics/16.9.847.
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.
Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.
Goudet J: FSTAT (version 1.2): a computer program to calculate F statistics. J Hered. 1995, 86: 485-486.
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.
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.
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.
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.
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.
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.
The authors declare that they have no competing interests.
PS and MI carried out the molecular genetic studies, performed the population genetic analysis and drafted the manuscript. NT, SKr and SL performed the clinical management and collected the patients' blood samples as well as helping with the writing of the manuscript. SKa participated in the design of the study, acquisition of funding, coordination and writing of the manuscript.
Electronic supplementary material
Additional file 1: Clinical data of 25 severe malaria patients at time of admission. RBC, red blood cell; Hb, hemoglobin; Ht, hematocrit; WBC, white blood cell; Plt, platelet; BUN, blood urea nitrogen; Cr, creatinine; Alb, albumin; AST, aspartate aminotransferase; ALT, alanine aminotransferase; PCT, parasite clearance time; FCT, fever clearance time. (XLS 42 KB)
Additional file 2: Clinical data of 25 uncomplicated malaria patients at time of admission. RBC, red blood cell; Hb, hemoglobin; Ht, hematocrit; WBC, white blood cell; Plt, platelet; BUN, blood urea nitrogen; Cr, creatinine; Alb, albumin; AST, aspartate aminotransferase; ALT, alanine aminotransferase; PCT, parasite clearance time; FCT, fever clearance time; NA, not available. (XLS 39 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.