- Open Access
PCR/RFLP-Based Analysis of Genetically Distinct Plasmodium vivax Population of Pvmsp-3α and Pvmsp-3β genes in Pakistan
Malaria Journalvolume 13, Article number: 355 (2014)
Plasmodium vivax is one of the widespread human malarial parasites accounting for 75% of malaria epidemics. However, there is no baseline information about the status and nature of genetic variation of Plasmodium species circulating in various parts of Pakistan. The present study was aimed at observing the molecular epidemiology and genetic variation of Plasmodium vivax by analysing its merozoite surface protein-3α (msp-3α) and merozoite surface protein-3β (msp-3β) genes, by using suballele, species-specific, combined nested PCR/RFLP detection techniques.
A total of 230 blood samples from suspected subjects tested slide positive for vivax malaria were collected from Punjab, Sindh, Khyber Pakhtunkhwa, and Balochistan during the period May 2012 to December 2013. Combined nested PCR/RFLP technique was conducted using Pvmsp-3α and Pvmsp-3β genetic markers to detect extent of genetic variation in clinical isolates of P. vivax in the studied areas of Pakistan.
By PCR, P. vivax, 202/230 (87.82%), was found to be widely distributed in the studied areas. PCR/RFLP analysis showed a high range of allelic variations for both msp-3α and msp-3β genetic markers of P. vivax, i.e., 21 alleles for msp-3α and 19 for msp-3β. Statistically a significant difference (p ≤ 0.05) was observed in the genetic diversity of the suballelic variants of msp-3α and msp-3β genes of P. vivax.
It is concluded that P. vivax populations are highly polymorphic and diverse allelic variants of Pvmsp-3α and Pvmsp-3β are present in Pakistan.
Malaria is the major threat to public health and economic development in many nations . It kills more than a million people a year, and approximately 40% of the world’s populations live in malarious countries . Plasmodium vivax is the most widespread species of human malaria parasites in the world and is endemic in many countries of Asia, Central and South America, the Middle East, and parts of Africa, with an estimated burden of 70-80 million cases annually [3, 4]. In Pakistan P. vivax is the common malarial species, contributing to 70% of the malaria burden . Although P. vivax is highly prevalent, it has received scant scientific attention in Pakistan, thus, a paucity of baseline data exists on various aspects of P. vivax, such as population structure and drug resistance patterns.
Genetic diversity is defined as the total number of heritable characteristics in the genetic make-up of the species . Studies have shown that Plasmodium species show diversity in several parameters, including morphology, biochemistry, relapse patterns, symptoms, course and duration of infection, immunological responses, drug resistance, and transmissibility by anopheline vectors. However, the extent of genetic diversity observed depends on the transmission rates, immune pressure and natural selection of the parasite in an area . Usually in endemic areas, where transmission rates are high, extensively diverse variants are in circulation, posing a threat to public health, since such variants have possibly acquired increased virulence and resistance to drugs in order to fit and survive .
Studies of the population structure of malaria parasites are important for understanding the evolution of parasite virulence and the role of parasite diversity in malaria transmission, and for designing control tools, including vaccines, as well as evaluating the impact of malaria control measures [9, 10].
Pvmsp-3α and Pvmsp-3β have been used as markers in population genetic studies worldwide [11–14]. The best approach for detecting genetic diversity is to analyse more than one marker gene, because by doing so the probability of different clones sharing the same genotype accidentally, is substantially reduced . Accordingly, it is suggested that analysis of both Pvmsp-3α and Pvmsp3β genes of P. vivax enables greater ability in identifying parasite haplotypes and detecting mixed strain infections [11, 12]. Limited data regarding the genetic diversity of Pvmsp-3α and Pvmsp-3β is available from one endemic area of Khyber Pakhtunkhwa province of Pakistan , it is therefore, hypothesized that new allelic variants may be present in other provinces of the country. This study will seek to elaborate on the genetic polymorphism of P. vivax malaria isolates and compare them with those reported from other parts of the world.
This study was conducted with the approval of the Medical Superintendents of source hospitals and the ethics committee of Kohat University of Science and Technology (KUST) Kohat Pakistan. A total of 230 blood samples were collected, after the informed consent from symptomatic patients, with the assessment of malaria control laboratories in the district headquarters hospitals of the four provinces: Punjab, Sindh, Khyber Pakhtunkhwa, and Balochistan during May 2012 to December 2013. Out of these 230 samples, 58 were collected from Punjab, 54 from Sindh, 68 from Khyber Pakhtunkhwa, and 50 from Balochistan. All the samples were screened for the detection of malarial parasites by microscope and PCR.
DNA of blood samples were extracted with the help of DNA extraction kit (GF1, Vivantus, USA).
PCR amplification for Plasmodium vivax identification
The nested PCR amplification procedures used by Snounou et al. and Singh et al.[17, 18] were applied for amplification and detection of the 18S rRNA type A gene size 1.2 Kb, using species-specific primers rVIV1 5′-CGCTTCTAGCTTAATCCACATAACTGATAC-3′ and rVIV2 5′-ACTTCCAAGCCGAAGCAAAGAAAGTCCTTA-3′ for the detection of P. vivax.
PCR of Pvmsp-3α genes
By using PCR technique, the allelic variants of P. vivax merozoite surface protein-3α were detected by standard protocol of Bruce et al.. Primary PCR of Pvmsp-3α was performed in 20 ml using 3 ml of DNA, using primers P1 (5′-CAGCAGACACCATTTAAGG-3′) and P2 (5′-CCGTTTGTTGATTAGTTGC-3′), while nested reactions were done in 20 μl with primers N1 (5′-GACCAGTGTGATACCATTAACC-3′) and N2 (5′-ATACTGGTTCTTCGTCTTCA GG-3′) using 3 μl of the primary PCR product. Taq DNA polymerase (2.6 units) was used and the following cyclic conditions were performed: 95°C for 3 min, 94°C for 30 sec, 58°C for 30 sec, 68°C for 2.5 min (2, 30 times) and 72°C for 5 min. Two per cent agarose gel was stained with ethidium bromide in order to visualize the PCR separated product, under UV illumination. Gel documentation system (Clearver Scientific, USA) was used for the photography.
RFLP analysis of Pvmsp-3α genes PCR products
For the RFLP analysis of Pvmsp-3α genes, Alu1 restriction enzyme (Fermentas, USA) was used. The total volume of the reaction medium was kept 20 μl having the following reagents Alu1 buffer (1.2 μl), PCR water (9.8 μl), Alu1 enzyme (1.0 μl) and product of msp- 3α gene (8.0 μl).
PCR of Pvmsp-3β genes
The allelic variants of P. vivax merozoite’s surface protein-3 were detected by PCR analysis, using standard protocol of Yang et al. . For Pvmsp-3β, primary PCR was performed in 20 μl using 3 μl of DNA using primers P1 (5′-GTATTCTTCGCAACACTC-3′) and P2 (5′-CTTCTGATGTTATTTCCAG-3′). Nested reactions were done in 20 ml with primers N1 (5′-CGAGGGGCGAAATTGTAAACC-3′) and N2 (5′-GCTGCTTCTTTTGCAAAGG-3′) using 1 μl of the primary PCR product as the template. Taq DNA polymerase (2.6 units) (Vivantis, USA) was used and the following cyclic conditions were performed: 94°C for 20 sec, 54°C for 30 sec, 68°C for 2.5 min, 35 cycles. Agarose was used to visualize the PCR product under UV illumination. Gel documentation system (Clearver Scientific, USA) was used for the photography.
RFLP analysis of Pvmsp-3β PCR products
The RFLP analysis of Pvmsp- 3β gene was carried out, using Pest 1 restriction enzyme (Fermentas, USA) described by Yang et al.. The total volume of the reaction medium was kept as 20 μl. The master mix for RFLP analysis was composed of Pest 1 Buffer 1.2 μl, PCR water 9.8 μl, Pest 1 enzyme 1.0 μl and product of msp-3β gene 8.0 μl.
The data were statistically analysed by using SPSS (Version 20) with application of Kruskal Wallis test. The p value less than 0.05 was assumed to be statistically significant.
Out of the total 230 samples, 202 samples were found positive for P. vivax both by microscopy and PCR. The distribution of the positive samples in the four provinces of Pakistan remained: 58/68 (85.29%) for Khyber Pakhtunkhwa, 53/58 (91.37%) for Punjab, 47/54 (87.03%) for Sindh and 44/50 (88%) for Balochistan (Figure 1).
Analysis of pvmsp- 3α gene by PCR /RFLP
A sum of 130/202 samples of blood were successfully amplified for pvmsp-3α gene and were identified for different genotypes of the aforementioned gene. Three major allelic variants of different sizes, i.e., 2.2 kb (Type-A), 2.0 kb (Type-B), 1.1 kb (Type-C) were detected. Besides these, an unusual band size of about 0.75 kb, named as Type-D and mixed genotypes designated as M were also observed (Figure 2). Out of the 130 amplified samples, the frequency of the different types of alleles remained 64/130 (49.29%) for Type-A, 27/130 (20.76%) for Type-B, 27/130 (20.76%) for Type-C, 9/130 (6.92 %) for Type-D, and 3/130 (2.30 %) for mixed. Overall 12 suballelic variants (A1-A12) for Type A, 5 (B1-B5) for Type B and 4 (C1-C4) for Type-C were observed. The alleles A3, A9, B1, B2, and C2 were frequently observed in all provinces of Pakistan (Table 1).
Analysis of pvmsp-3β gene by PCR/RFLP
In the current study 72/202 samples of blood had been successfully amplified for pvmsp-3β gene with a high frequency of polymorphism among their different allelic variants. The amplified DNA fragments for these alleles were classified according to their size variation. Based on variation in size of the fragments, the amplified samples could be differentiated into three allele sizes: for Type A (2.0-2.5 kb), Type B (1.6-2.0 kb) and Type C (1.5 kb) (Figure 3).
In all these amplified samples, a couple of samples reflect more than one band of different sizes. These multiple bands show mixed infection, which revealed 2.77% of the allelic variants. Out of the total 72 amplified samples, 16 (22.22%) isolates for Type-A, 32 (44.44%) for Type-B and 22 (30.55%) for Type-C were observed. Among these isolates Type-B was found to have highest per cent frequency and Type-A was found the lowest in per cent frequency (Table 2).
In this study pvmsp-3β gene was analysed using Pst 1 restriction enzyme and 7 suballelic variants, (A1-A7) for Type-A, 10 (B1-B10) for Type-B and 6 (C1-C6) for Type-C were observed. The suballeles B7-B10 were not found in the present study (Table 2). The different samples resolved through pst 1 restriction enzyme, show a high polymorphism in size and distribution among samples (Figure 3).
Pakistan is endemic for malaria and there are reports showing the distribution of various allelic variants of P. vivax. So, this study was aimed with the main objectives of identifying the new variants and to observe the epidemiology of level of pvmsp-3α and pvmsp-3β polymorphism of P. vivax in Pakistan. In the current study four distinct sizes of PCR products for Pvmsp-3α were detected. Among these, three different major allelic variants of different sizes, i.e., 2.2 kb (Type-A), 2.0 kb (Type-B), 1.1 kb (Type-C) along with unusual band size of almost 0.75 kb or 300 bp name as Type-D and mixed genotypes designated as M, were observed in the study areas, while in Thailand , Afghanistan  and the Federally Administered Tribal Areas (FATA) of Pakistan , only three allele sizes (A, B and C) were obtained for this gene.
The frequency of the different types of alleles remained 64/130 (49.29%) for Type-A, 27/130 (20.76%) for Type-B, 27/130 (20.76%) for Type-C, 9/130 (6.92%) for Type-D, and 3/130 (2.30%) for mixed, comparing with the study of Khatoon et al. where 82% (41/50) isolates having allele sizes 1.9 kb (Type-A), 6% (3/50) with the 1.5 kb fragment (Type-B), 8% (4/50) with the 1.2 kb fragment (Type-C), 2% (1/50) with the 0.3 kb fragment (Type-D) and one isolate (2%) exhibited mixed-strain infection .
In this study restriction digestion of the Pvmsp-3α PCR product with Alu I yielded 12 suballelic variants types: (A1-A12) for Type-A, 5 (B1-B5) for Type-B and 4 (C1-C4) for Type-C were observed. The alleles A3 11 (8.46%) and A9 13 (10%) for Type-A, allele B1 and B2 11 (8.46%) for Type-B and C2 15 (11.53%) for Type-C were frequently observed in all provinces of Pakistan, while in the similar study carried out by Khatoon et al., 12 different alleles designated as A1-A7, B1, B2, C1, C2, and D, allele A3 being the most abundant 12/50 (24%) .
The PCR-RFLP patterns for P. vivax Pvmsp-3β gene produced three categories Type-A (2.0-2.5 kb), Type-B (1.6-2.0 kb) and Type-C (1.5 kb). Similarly, in northwest Thailand, three distinct allele types (A, B and C (~0.65 kb)) were found . In contrast, studies carried out in district Bannu by Khatoon et al. found two types of allele A (1.7-2.2 kb) and B (1.4-1.5 kb) of size polymorphisms for P. vivax Pvmsp-3β gene. Compared with P. vivax isolates from previous work conducted in Asia, the results of the current study further confirmed the existence of small geographic differentiation among P. vivax populations. An earlier study in western Thailand showed that the B type of pvmsp-3β is more abundant (60.4%) than other types, whereas in Chinese Bengbu and Guangxi samples, both A and B types were similarly prevalent . According to Zhong et al., the Type-A allele was the most abundant in all four parasite populations (>57% more abundant than other types), i.e., Anhui, Hainan, Yunnan, and Myanmar . In the current study, Type-B was found to have the highest per cent frequency and Type-A was found to have lowest per cent frequency.
Data from PCR-RFLP of the Pvmsp-3α and Pvmsp-3β loci showed that 5.02% of the Pakistani P. vivax isolates exhibited mixed-strain infections, which is comparable to sites in Bengbu (China) (5.6%) and is lower than observed in Thailand (20.5%)  and FATA of Pakistan (30%), which directly shares a border with Afghanistan . In northern Iran  and Hongshuihe (China), no mixed genotypes were detected . Apart from differences in transmission intensities in different areas, the observed variations in parasite types could be attributed to factors such as sampling bias, host immune selective pressure on particular types and/or spatiotemporal changes in the availability of different mosquito species that can transmit specific parasite types in a particular area over different times or seasons, with areas that share the same mosquito types tending to have similar parasite types .
This study indicates that P.vivax populations in four provinces of Pakistan: Punjab, Sindh, Khyber Pakhtunkhwa, and Balochistan, are highly diverse. Such heterogeneity might cause differences in parasite virulence, transmissibility and responses to chemotherapy, with important implications for malaria control measures in this populous region.
Carlton JM, Adams JH, Silva JC, Bidwell SL, Lorenzi H, Caler E, Crabtree J, Angiuoli SV, Merino EF, Amedeo P, Cheng Q, Coulson RM, Crabb BS, Del Portillo HA, Essien K, Feldblyum TV, Fernandez-Becerra C, Gilson PR, Gueye AH, Guo X, Kang’a S, Kooij TW, Korsinczky M, Meyer EV, Nene V, Paulsen I, White O, Ralph SA, Ren Q, Sargeant TJ: Comparative genomics of the neglected human malaria parasite Plasmodium vivax. Nature. 2008, 455: 757-763. 10.1038/nature07327.
Tsuboi T, Takeo S, Arumugam TU, Otsuki H, Torii M: The wheat germ cell-free protein synthesis system: a key tool for novel malaria vaccine candidate discovery. Acta Trop. 2010, 114: 171-176. 10.1016/j.actatropica.2009.10.024.
Mendis K, Sina BJ, Marchesini P, Carter R: The neglected burden of Plasmodium vivax malaria. Am J Trop Med Hyg. 2001, 64: 97-106.
Guerra CA, Snow RW, Hay SI: Mapping the global extent of malaria in 2005. Trends Parasitol. 2006, 22: 353-358. 10.1016/j.pt.2006.06.006.
WHO: World Malaria Report 2011. 2011, Geneva: World Health Organization,http://www.who.int/malaria/world_malaria_report_2011/en/,
Genetic diversity [Internet]: 2010, Available from: en.wikipedia.org/ wiki /Genetic_diversity
Escalante AA, Rojas A, Udhayakumar V, Lal AA: Assessing the effect of natural selection in malaria parasites. Trends Parasitol. 2004, 20: 388-395. 10.1016/j.pt.2004.06.002.
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.
Thakur A, Alam MT, Sharma YD: Genetic diversity in the C-terminal 42 kDa region of merozoite surface protein-1 of Plasmodium vivax (PvMSP-1(42)) among Indian isolates. Acta Trop. 2008, 108: 58-63. 10.1016/j.actatropica.2008.08.011.
Mzilahowa T, McCall PJ, Hastings IM: “Sexual” population structure and genetics of the malaria agent P. falciparum. PLoS ONE. 2007, 2: e613-10.1371/journal.pone.0000613.
Bruce MC, Galinski MR, Barnwell JW, Snounou G, Day KP: Polymorphism at the merozoite surface protein-3alpha locus of Plasmodium vivax: global and local diversity. Am J Trop Med Hyg. 1999, 61: 518-525.
Yang Z, Miao J, Huang Y, Li X, Putaporntip C, Jongwutiwes S, Gao Q, Udomsangpetch R, Sattabongkot J, Cui L: Genetic structures of geographically distinct Plasmodium vivax populations assessed by PCR/RFLP analysis of the merozoite surface protein 3beta gene. Acta Trop. 2006, 100: 205-212. 10.1016/j.actatropica.2006.10.011.
Cui L, Mascorro CN, Fan Q, Rzomp KA, Khuntirat B, Zhou G, Chen H, Yan G, Sattabongkot J: Genetic diversity and multiple infections of Plasmodium vivax malaria in Western Thailand. Am J Trop Med Hyg. 2003, 68: 613-619.
Mueller I, Kaiok J, Reeder JC, Cortes A: The population structure of Plasmodium falciparum and Plasmodium vivax during an epidemic of malaria in the Eastern Highlands of Papua New Guinea. Am J Trop Med Hyg. 2002, 67: 459-464.
Schoepflin S, Valsangiacomo F, Lin E, Kiniboro B, Mueller I, Felger I: Comparison of Plasmodium falciparum allelic frequency distribution in different endemic settings by high-resolution genotyping. Malar J. 2009, 8: 250-10.1186/1475-2875-8-250.
Khatoon L, Baliraine FN, Bonizzoni M, Kalik SA, Yan G: Genetic structure of Plasmodium vivax and Plasmodium falciparum in the Bannu district of Pakistan. Malar J. 2010, 9: 112-10.1186/1475-2875-9-112.
Snounou G, Viriyakosol S, Zhu XP, Jarra W, Pinheiro L, Do Rosario VE, Thaithong S, Brown KN: High sensitivity of detection of human malaria parasites by the use of nested polymerase chain reaction. Mol Biochem Parasitol. 1993, 61: 315-320. 10.1016/0166-6851(93)90077-B.
Singh B, Bobogare A, Cox-Singh J, Snounou G, Abdullah MS, Rahman HA: A genus- and species-specific nested polymerase chain reaction malaria detection assay for epidemiologic studies. Am J Trop Med Hyg. 1999, 60: 687-692.
Zakeri S, Safi N, Afsharpad M, Butt W, Ghasemi F, Mehrizi AA, Atta H, Zamani G, Djadid ND: Genetic structure of Plasmodium vivax isolates from two malaria endemic areas in Afghanistan. Acta Trop. 2009, 113: 12-19.
Zakeri S, Raeisi A, Afsharpad M, Kakar Q, Ghasemi F, Atta H, Zamani G, Memon MS, Salehi M, Djadid ND: Molecular characterization of Plasmodium vivax clinical isolates in Pakistan and Iran using pvmsp-1, pvmsp-3alpha and pvcsp genes as molecular markers. Parasitol Int. 2010, 59: 15-21. 10.1016/j.parint.2009.06.006.
Zhong D, Bonizzoni M, Zhou G, Wang G, Chen B, Vardo-Zalik A, Cui L, Yan G, Zheng B: Genetic diversity of Plasmodium vivax malaria in China and Myanmar. Infect Genet Evol. 2011, 11: 1419-1425. 10.1016/j.meegid.2011.05.009.
We thank staff at the malaria clinics and malaria laboratories in different hospitals for assistance in sample collection, and thank the patients for participating in the study.
The authors declare that they have no competing interests.
SNK and SK designed the study and prepared the manuscript. AK and MAK performed the laboratory work. SA, IA, SA, analysed and interpreted the data. JK and AHS review the MS critically. All authors read and approved the final version of the manuscript.
Shahid Niaz Khan, Asif Khan contributed equally to this work.