Genetic diversity among Plasmodium falciparum field isolates in Pakistan measured with PCR genotyping of the merozoite surface protein 1 and 2
- Najia K Ghanchi†1, 2,
- Andreas Mårtensson†3, 4,
- Johan Ursing3,
- Sana Jafri1,
- Sándor Bereczky3, 5,
- Rabia Hussain1 and
- Mohammad A Beg1Email author
© Ghanchi et al; licensee BioMed Central Ltd. 2010
Received: 9 July 2009
Accepted: 1 January 2010
Published: 1 January 2010
The genetic diversity of Plasmodium falciparum has been extensively studied in various parts of the world. However, limited data are available from Pakistan. This study aimed to establish molecular characterization of P. falciparum field isolates in Pakistan measured with two highly polymorphic genetic markers, i.e. the merozoite surface protein 1 (msp-1)and 2 (msp-2).
Between October 2005 and October 2007, 244 blood samples from patients with symptomatic blood-slide confirmed P. falciparum mono-infections attending the Aga Khan University Hospital, Karachi, or its collection units located in Sindh and Baluchistan provinces, Pakistan were collected. The genetic diversity of P. falciparum was analysed by length polymorphism following gel electrophoresis of DNA products from nested polymerase chain reactions (PCR) targeting block 2 of msp-1 and block 3 of msp-2, including their respective allelic families KI, MAD 20, RO33, and FC27, 3D7/IC.
A total of 238/244 (98%) patients had a positive PCR outcome in at least one genetic marker; the remaining six were excluded from analysis. A majority of patients had monoclonal infections. Only 56/231 (24%) and 51/236 (22%) carried multiple P. falciparum genotypes in msp-1 and msp-2, respectively. The estimated total number of genotypes was 25 msp-1 (12 KI; 8 MAD20; 5 RO33) and 33 msp-2 (14 FC27; 19 3D7/IC).
This is the first report on molecular characterization of P. falciparum field isolates in Pakistan with regards to multiplicity of infection. The genetic diversity and allelic distribution found in this study is similar to previous reports from India and Southeast Asian countries with low malaria endemicity.
Plasmodium falciparum is an important public health problem in Pakistan causing at least 500,000 cases of malaria annually . An increasing number of cases are being reported every year and the relative frequency of P. falciparum amongst slide positive malaria cases has increased from 45% in 1995 to 68% in 2006 [2–4]. In line with this, the National Malaria Control Programme has reported a six-fold increase in P. falciparum during the last decade. Despite this malaria has received scarce scientific attention, particularity with regards to molecular characterization of the local parasite population.
The introduction of PCR genotyping in malaria research has paved the way for major improvements in the understanding of parasite biology. For example, a number of highly diverse genetic markers of P. falciparum have been identified and extensively studied as potential vaccine candidates [5, 6]. The most commonly used markers for genotyping of P. falciparum are the surface antigens merozoite surface protein 1 (msp-1), merozoite surface protein 2 (msp-2) and the glutamate-rich protein (glurp) [7, 8]. These three genetic markers are unlinked, i.e. located on different chromosomes, single copy genes with extensive polymorphism, both with regards to sequence and size, which is mostly generated by intragenic repeats that are variable in copy number and length of the repeat unit . These features make them attractive candidates for studies where identification and enumeration of genetically distinct P. falciparum parasite sub-populations are of interest. As such they have proven to be useful tools both in molecular epidemiology studies in different epidemiological settings as well as to distinguish treatment failures from new infections in anti-malarial drug trials [10, 11].
The genetic diversity among the P. falciparum population is an important indicator of the malaria transmission intensity in an area [12, 13]. A high endemic area is generally characterized by extensive parasite diversity and infected humans often carry multiple genotypes. Conversely, the parasite population in a low transmission area has a limited genetic diversity and the majority of infections are monoclonal [14–17].
Plasmodium falciparum field isolates have been characterized in both Iran and India, but not previously in Pakistan, using the above mentioned molecular markers. This study, therefore, aimed to bridge this gap by molecular characterization of P. falciparum field isolates in Pakistan measured using two highly polymorphic markers, i.e. msp-1 and msp-2.
Study setting, participants, and ethics
Between October 2005 and October 2007, patients of all ages presenting with microscopy confirmed asexual P. falciparum mono-infections, regardless of parasite density and disease severity, at the Aga Khan University Hospital, Karachi or its established chain of health centres/collection units located in Sindh and Baluchistan provinces, Pakistan were enrolled in this study.
The study was conducted in accordance with the latest version of the Declaration of Helsinki and Good Clinical Practice . Informed consent was obtained from all participants or in case of children from their parents/legal guardians. The study was approved by the ethical review committee of Aga University Hospital, Karachi, Pakistan.
Blood collection and microscopy
An intravenous blood sample of 2 ml was collected in an EDTA tube in accordance with routine clinical practice for all patients referred for laboratory investigation of malaria infection. Initially, a blood smear was microscopically analysed using Leishman staining for screening of presence of malaria parasites . In case of a positive screening slide a confirmatory thick and thin Giemsa-stained blood film was analysed for species identification and parasite density. For patients with confirmed P. falciparum mono-infection asexual parasites were counted against 200 white blood cells on the thick film. Parasite density was quantified (parasites/μl) by assuming an average of 8,000 leucocytes per μl blood . All blood slides were examined by experienced microscopists at the clinical laboratory, Aga Khan University. Quality control was done for 10% of the slides by an independent microscopist blinded to the initial report.
For all patients enrolled in the study the remaining blood volume collected in the EDTA tube was transferred to cryovials and kept frozen at -80°C until used for DNA extraction. A brief epidemiological/demographic history was also collected from each participant using a structured questionnaire.
DNA extraction and molecular analyses
Genomic DNA was extracted from a total of 200 μl whole blood per patient using the Qiagen DNA extraction kit (QIAGEN, USA) according to manufacturer's instructions. Nested PCR genotyping was performed both for the variable block 2 regions of msp-1 and block 3 of msp-2, considered to be the two most informative genetic markers for assessment of multiplicity of P. falciparum infection .
The initial amplification was followed by individual nested PCR reactions using family specific primers for msp-1 (KI, MAD20 and R033), and msp-2 (FC27 and 3D7/IC), respectively, based on previously described standard protocols . Positive and negative controls were systematically incorporated in each PCR run. The msp-1 and msp-2 PCR products were loaded on 3% and 2% agarose gels, respectively, stained with ethidium-bromide, separated by electrophoresis and visualized under UV trans-illumination (GelDoc®, Biorad, Hercules, USA).
Analyses of number of genotypes and size polymorphism were digitalized using Quantity One® software (Biorad, Hercules, USA). For assessment of the overall number of genotypes present within the P. falciparum population and their respective prevalence, individual genotypes were defined arbitrarily by binning 20 base-pair (bp) intervals together. The median genotype for each family of the respective genetic markers was identified. The absolute size of the identified median band +/- 10 bp formed the initial bin. Thereafter, each 20 bp interval below and above the median band were defined as representing a distinct genotype.
Study design and statistical analyses
This was a descriptive/exploratory study, which precludes a power calculation of sample size. Data were entered in Microsoft Excel and exported to SPSS 15.0 software for analyses. Arithmetical means and medians, where applicable, were calculated for all continuous baseline demographic variables, except for asexual parasite density (geometric mean). Proportions of binary outcomes were compared using non-parametric tests. Spearman's rank correlation coefficients were calculated to assess association between multiplicity of infection and parasite densities. Statistical significance was defined as a P value ≤ 0.05.
Baseline demographic data
A total of 244 patients with microscopy confirmed P. falciparum mono-infection were enrolled. Six (2.5%) were excluded from the analysis due to negative PCR outcome in both the msp-1 and msp-2 marker. The remaining 238 patients all presented with symptomatic malaria infection, 164 (69%) with uncomplicated and the remaining 74 (31%) with complicated/severe disease.
The geographical distribution of the 238 patients included in the analysis was: 176 (74%) from Karachi, 21 (9%) from Sukkur, 31 (13%) from rural areas of Sindh province and the remaining 10 (4%) from Baluchistan. The asexual P. falciparum parasitaemia ranged from 50 to 500,000 parasites/μl with a geometric mean of 8,412 parasites/μl.
Plasmodium falciparum allelic diversity
Number of Plasmodium falciparum genotypes and base pair ranges observed in msp-1 and msp-2, as well as their respective allelic families
Base pair range
No. of genotypes
Total no. of genotypes
Overall multiplicity of infection
Association between mean multiplicity of infection, parasite densities and age
Multiplicity of Plasmodium falciparum infection as assessed with the msp-1 and msp-2 marker in different geographic locations of Southern Pakistan
Number of PCR fragments
Multiplicity of infection*
Multiplicity of infection*
Plasmodium falciparum parasite density stratified by number of msp-1 and msp-2 genotypes per infection.
Number of PCR fragments*
Plasmodium falciparum parasite density and mean multiplicity of infection in msp-1 and msp-2 stratified by age
This is the first report on genetic diversity of P. falciparum field isolates from Pakistan and as such, it bridges an important gap in the understanding of the molecular characteristics of the parasite population in South Asia.
The results reveal a relatively limited genetic diversity and multiplicity of infection. These findings are compatible with the general belief that Pakistan, despite potential substantial regional variations and recent increases in P. falciparum incidence, represents a low malaria transmission setting [13, 14, 16, 22, 23]. Moreover, the results are in line with several previous reports from India, Iran and Southeast Asian countries with low malaria transmission . In contrast, the results of this study do not corroborate recent finding from Iran where Zakeri et al reported a relatively high multiplicity of infection in both the msp-1 and msp-2 marker from a presumably low endemic area .
However, it should be emphasized that a majority of patients (74%) presented with symptomatic malaria infection at a tertiary level of care, i.e. the Aga Khan University Hospital. Consequently, some of these patients may have received anti-malarial treatment prior to enrolment in the study. As treatment is likely to reduce the number of genotypes in an infected individual, the P. falciparum identified in these patients may not be representative for the Pakistani parasite population. Moreover, reports of genetic diversity from low endemic areas generally study symptomatic P. falciparum infections, since asymptomatic carriers rarely are observed in setting where semi-immunity cannot be acquired . This is of importance since data indicate that symptomatic infections generally appear to harbour a lower multiplicity of infection as compared to asymptomatic children residing in high transmission areas in Africa . Lastly, the inability to detect all parasite sub-populations present in an individual with one single blood sample, either due to inborn limitations of the PCR technique to detect minority clones or parasite population dynamics, may also be considered as a factor potentially underestimating the multiplicity of infection retrieved in our study . Thus, the results of this study may reflect a significant underestimation of the genetic diversity of the P. falciparum population in Pakistan.
The reported study represents a first attempt to establish molecular techniques for malaria research at the Aga Khan University, Karachi. Future studies will be designed to geographically expand the collection of P. falciparum field isolates to include all four provinces of Pakistan. During the next phase, analyses of genetic markers related to anti-malarial drug resistance will also be included to gain a more thorough picture of P. falciparum molecular epidemiology in Pakistan.
The genetic diversity and allelic distribution found in this study is similar to previous reports from India and Southeast Asian countries with low malaria endemicity.
We thank Malaria Research and Reference Reagent Resource Centre (MR4), ATCC® Manassas, VA, USA for providing us positive controls (3D7, KI, Dd2 genomic DNA). This study was financially supported by the University Research Council grant (URC # 051012 P&M) by Aga Khan University (MAB), the Swedish South Asian Studies Network (SASNET) (AM), the Swedish institute (AM, JU and NKG) and Higher Education Commission, Government of Pakistan (NKG) (Ref # 042-130358-Ls2305).
- Yasinzai MI, Kakarsulemankhel JK: Prevalence of human malaria infection in bordering areas of East Balochistan, adjoining with Punjab: Loralai and Musakhel. J Pak Med Assoc. 2009, 59: 132-135.PubMedGoogle Scholar
- Durrani AB, Durrani IU, Abbas N, Jabeen M: Epidemiology of cerebral malaria and its mortality. J Pak Med Assoc. 1997, 47: 213-215.PubMedGoogle Scholar
- Rab MA, Freeman TW, Durrani N, de Poerck D, Rowland MW: Resistance of Plasmodium falciparum malaria to chloroquine is widespread in eastern Afghanistan. Ann Trop Med Parasitol. 2001, 95: 41-46. 10.1080/00034980020035906.View ArticlePubMedGoogle Scholar
- Hozhabri S, Akhtar S, Rahbar MH, Luby SP: Prevalence of plasmodium slide positivity among the children treated for malaria, Jhangara, Sindh. J Pak Med Assoc. 2000, 50: 401-405.PubMedGoogle Scholar
- Diggs CL, Ballou WR, Miller LH: The major merozoite surface protein as a malaria vaccine target. Parasitol Today. 1993, 9: 300-302. 10.1016/0169-4758(93)90130-8.View ArticlePubMedGoogle Scholar
- Greenhouse B, Myrick A, Dokomajilar C, Woo JM, Carlson EJ, Rosenthal PJ, Dorsey G: Validation of microsatellite markers for use in genotyping polyclonal Plasmodium falciparum infections. Am J Trop Med Hyg. 2006, 75: 836-842.PubMed CentralPubMedGoogle Scholar
- Snounou G, Beck HP: The use of PCR genotyping in the assessment of recrudescence or reinfection after antimalarial drug treatment. Parasitol Today. 1998, 14: 462-467. 10.1016/S0169-4758(98)01340-4.View ArticlePubMedGoogle Scholar
- Smythe JA, Coppel RL, Day KP, Martin RK, Oduola AM, Kemp DJ, Anders RF: Structural diversity in the Plasmodium falciparum merozoite surface antigen 2. Proc Natl Acad Sci USA. 1991, 88: 1751-1755. 10.1073/pnas.88.5.1751.PubMed CentralView ArticlePubMedGoogle Scholar
- Farnert A, Arez AP, Babiker HA, Beck HP, Benito A, Bjorkman A, Bruce MC, Conway DJ, Day KP, Henning L, Mercereau-Puijalon O, Ranford-Cartwright LC, Rubio JM, Snounou G, Walliker D, Zwetyenga J, Do Rosario VE: Genotyping of Plasmodium falciparum infections by PCR: a comparative multicentre study. Trans R Soc Trop Med Hyg. 2001, 95: 225-232. 10.1016/S0035-9203(01)90175-0.View ArticlePubMedGoogle Scholar
- Cattamanchi A, Kyabayinze D, Hubbard A, Rosenthal PJ, Dorsey G: Distinguishing recrudescence from reinfection in a longitudinal antimalarial drug efficacy study: comparison of results based on genotyping of msp-1, msp-2, and glurp. Am J Trop Med Hyg. 2003, 68: 133-139.PubMedGoogle Scholar
- Collins WJ, Greenhouse B, Rosenthal PJ, Dorsey G: The use of genotyping in antimalarial clinical trials: a systematic review of published studies from 1995-2005. Malar J. 2006, 5: 122-10.1186/1475-2875-5-122.PubMed CentralView ArticlePubMedGoogle Scholar
- Babiker HA, Charlwood JD, Smith T, Walliker D: Gene flow and cross-mating in Plasmodium falciparum in households in a Tanzanian village. Parasitology. 1995, 111: 433-442. 10.1017/S0031182000065938.View ArticlePubMedGoogle Scholar
- Paul RE, Hackford I, Brockman A, Muller-Graf C, Price R, Luxemburger C, White NJ, Nosten F, Day KP: Transmission intensity and Plasmodium falciparum diversity on the northwestern border of Thailand. Am J Trop Med Hyg. 1998, 58: 195-203.PubMedGoogle Scholar
- Haddad D, Snounou G, Mattei D, Enamorado IG, Figueroa J, Stahl S, Berzins K: Limited genetic diversity of Plasmodium falciparum in field isolates from Honduras. Am J Trop Med Hyg. 1999, 60: 30-34.PubMedGoogle Scholar
- Babiker HA, Lines J, Hill WG, Walliker D: Population structure of Plasmodium falciparum in villages with different malaria endemicity in east Africa. Am J Trop Med Hyg. 1997, 56: 141-147.PubMedGoogle Scholar
- Gomez D, Chaparro J, Rubiano C, Rojas MO, Wasserman M: Genetic diversity of Plasmodium falciparum field samples from an isolated Colombian village. Am J Trop Med Hyg. 2002, 67: 611-616.PubMedGoogle Scholar
- Peyerl-Hoffmann G, Jelinek T, Kilian A, Kabagambe G, Metzger WG, von Sonnenburg F: Genetic diversity of Plasmodium falciparum and its relationship to parasite density in an area with different malaria endemicities in West Uganda. Trop Med Int Health. 2001, 6: 607-613. 10.1046/j.1365-3156.2001.00761.x.View ArticlePubMedGoogle Scholar
- Williams JR: The Declaration of Helsinki and public health. Bull World Health Organ. 2008, 86: 650-652. 10.2471/BLT.08.050955.PubMed CentralView ArticlePubMedGoogle Scholar
- Warhurst DC, Williams JE: ACP Broadsheet no 148. July 1996. Laboratory diagnosis of malaria. J Clin Pathol. 1996, 49: 533-538. 10.1136/jcp.49.7.533.PubMed CentralView ArticlePubMedGoogle Scholar
- Moody A: Rapid diagnostic tests for malaria parasites. Clin Microbiol Rev. 2002, 15: 66-78. 10.1128/CMR.15.1.66-78.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Zakeri S, Bereczky S, Naimi P, Pedro Gil J, Djadid ND, Farnert A, Snounou G, Bjorkman A: Multiple genotypes of the merozoite surface proteins 1 and 2 in Plasmodium falciparum infections in a hypoendemic area in Iran. Trop Med Int Health. 2005, 10: 1060-1064. 10.1111/j.1365-3156.2005.01477.x.View ArticlePubMedGoogle Scholar
- Joshi H, Valecha N, Verma A, Kaul A, Mallick PK, Shalini S, Prajapati SK, Sharma SK, Dev V, Biswas S, Nanda N, Malhotra MS, Subbarao SK, Dash AP: Genetic structure of Plasmodium falciparum field isolates in eastern and north-eastern India. Malar J. 2007, 6: 60-10.1186/1475-2875-6-60.PubMed CentralView ArticlePubMedGoogle Scholar
- Ariey F, Chalvet W, Hommel D, Peneau C, Hulin A, Mercereau-Puijalon O, Duchemin JB, Sarthou JL, Reynes JM, Fandeur T: Plasmodium falciparum parasites in French Guiana: limited genetic diversity and high selfing rate. Am J Trop Med Hyg. 1999, 61: 978-985.PubMedGoogle Scholar
- Heidari A, Keshavarz H, Rokni MB, Jelinek T: Genetic diversity in merozoite surface protein (MSP)-1 and MSP-2 genes of Plasmodium falciparum in a major endemic region of Iran. Korean J Parasitol. 2007, 45: 59-63. 10.3347/kjp.2007.45.1.59.PubMed CentralView ArticlePubMedGoogle Scholar
- Beck HP, Felger I, Huber W, Steiger S, Smith T, Weiss N, Alonso P, Tanner M: Analysis of multiple Plasmodium falciparum infections in Tanzanian children during the phase III trial of the malaria vaccine SPf66. J Infect Dis. 1997, 175: 921-926. 10.1086/513991.View ArticlePubMedGoogle Scholar
- Farnert A, Rooth I, Svensson , Snounou G, Bjorkman A: Complexity of Plasmodium falciparum infections is consistent over time and protects against clinical disease in Tanzanian children. J Infect Dis. 1999, 179: 989-995. 10.1086/314652.View ArticlePubMedGoogle Scholar
- Martensson A, Ngasala B, Ursing J, Isabel Veiga M, Wiklund L, Membi C, Montgomery SM, Premji Z, Farnert A, Bjorkman A: Influence of consecutive-day blood sampling on polymerase chain reaction-adjusted parasitological cure rates in an antimalarial-drug trial conducted in Tanzania. J Infect Dis. 2007, 195: 597-601. 10.1086/510910.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.