PCR/RFLP-Based Analysis of Genetically Distinct Plasmodium vivax Population of Pvmsp-3α and Pvmsp-3β genes in Pakistan

Background 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. Methods 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. Results 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. Conclusion It is concluded that P. vivax populations are highly polymorphic and diverse allelic variants of Pvmsp-3α and Pvmsp-3β are present in Pakistan.


Background
Malaria is the major threat to public health and economic development in many nations [1]. It kills more than a million people a year, and approximately 40% of the world's populations live in malarious countries [2]. 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 [5]. 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 [6]. 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 [7]. 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 [8].
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][12][13][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 [15]. 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 [16], 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.

Sample collection
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 purification
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′-CGCTTCTAGC TTAATCCACATAACTGATAC-3′ and rVIV2 5′-ACTTC CAAGCCGAAGCAAAGAAAGTCCTTA-3′ for the detection of P. vivax.

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. [12]. 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.

Statistical analysis
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.  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 Pst1 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 pst1 restriction enzyme, show a high polymorphism in size and distribution among samples ( Figure 3).

Discussion
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 [13], Afghanistan [19] and the Federally Administered Tribal Areas (FATA) of Pakistan [20], only three allele sizes (A, B and C) were obtained for this gene.   fragment (Type-D) and one isolate (2%) exhibited mixedstrain infection [16]. In this study restriction digestion of the Pvmsp-3α PCR product with AluI 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%) [16].
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 [12]. 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 [16]. 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 [12]. 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 [21]. 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%) [12] and FATA of Pakistan (30%), which directly shares a border with Afghanistan [20]. In northern Iran [20] and Hongshuihe (China), no mixed genotypes were detected [12]. 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 [20].

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