Skip to main content

Genetic diversity of Plasmodium vivax population in Anhui province of China

Abstract

Background

Although the numbers of malaria cases in China have been declining in recent years, outbreaks of Plasmodium vivax malaria were still being reported in rural areas south of the Yellow River. To better understand the transmission dynamics of P. vivax parasites in China, the extent of genetic diversity of P. vivax populations circulating in Bozhou of Anhui province of China were investigated using three polymorphic genetic markers: merozoite surface proteins 1 and 3α (pvmsp-1 and pvmsp-3α) and circumsporozoite protein (pvcsp).

Methods

Forty-five P. vivax clinical isolates from Bouzhou of Anhui province were collected from 2009 to 2010 and were analysed using PCR/RFLP or DNA sequencing.

Results

Seven and six distinct allelic variants were identified using PCR/RFLP analysis of pvmsp-3α with Hha I and Alu I, respectively. DNA sequence analysis of pvmsp-1 (variable block 5) revealed that there were Sal-I and recombinant types but not Belem type, and seven distinct allelic variants in pvmsp-1 were detected, with recombinant subtype 2 (R2) being predominant (66.7%). All the isolates carried pvcsp with VK210 type but not VK247 or P. vivax-like types in the samples. Sequence analysis of pvcsp gene revealed 12 distinct allelic variants, with VK210-1 being predominant (41.5%).

Conclusions

The present data indicate that there is some degree of genetic diversity among P. vivax populations in Anhui province of China. The genetic data obtained may assist in the surveillance of P. vivax infection in endemic areas or in tracking potential future disease outbreak.

Background

Malaria kills more than a million people a year, and approximately 40% of the world’s populations live in malarious countries [1]. 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 [2, 3]. In China, the malaria eradication campaign initiated in the 1990s, using vector control and drug treatment of febrile individuals had been highly successful. However, episodes of P. vivax malaria infection were still being reported, and vivax malaria re-emerged in many counties of central China such as provinces of Anhui, Shandong, Hubei, Henan, and Jiangsu in recent years [4]. Approximately 37% to 68% of malaria cases reported during 2003 to 2007 were from central China, where the predominant vector mosquito was Anopheles sinensis. Most importantly, the numbers of vivax malaria cases in Anhui province of China were the highest in the country during those years [5].

A countrywide malaria elimination policy was launched by the Ministry of Health of P.R. China in 2010, with the goal to eliminate malaria by 2015 in a majority of the regions with the exception of the border region in Yunnan province, and to completely eliminate malaria from China by 2020. To achieve this ambitious goal, studying population structure, genetic diversity, and transmission of Plasmodium parasites in endemic areas will provide important information for disease control and management, including baseline data essential for monitoring drug resistance and for predicting the origin and spread of parasite variants within and between populations, and the performance of malaria vaccines under development [6, 7]. Several genetic markers including microsatellites and the genes encoding circumsporozoite protein (csp), merozoite surface protein 1 (msp-1), and merozoite surface proteins 3α and 3β (msp-3α and msp-3β), have been used to study P. vivax population diversity [811]. However, so far the population structures of P. vivax are less well studied and understood and only a limited number of studies on the population genetics of P. vivax isolates in China have been conducted [7]. In this study, the population diversity of clinical P. vivax isolates from Anhui province of central China was evaluated using the above three polymorphic markers, and the molecular diversity data may provide important information for vivax malaria surveillance and for tracking parasites in future outbreaks in this area.

Methods

Study area

The study was conducted in Bozhou city of Anhui province, which borders six provinces: the Jiangsu province in the east, the Zhejiang province in the southeast, the Jiangxi province in the south, the Hubei province in the west, the Henan province in the northwest, and the Shandong province in the north (Figure 1), located in an area between latitudes 32°51′ and 35°05′, and longitudes 115°53′ and 116°49′. The population of this city was estimated to be 485,000 in 2010. The area has a distinct seasonal climate characterized by humid summer from June to September and dry winter from October to December, with temperature ranging from 14°C to 17°C and rainfall ranging from 500 mm to 750 mm per year [12]. Bozhou city was a vivax endemic area with seasonal transmission peaking from June to September. Climatic and ecological conditions of this area make the environment favorable for mosquito breeding, and the main malaria vector is An. sinensis and the secondary one is Anopheles anthropophagus.

Figure 1
figure 1

Map of Anhui province in central China showing the location of Bozhou city where P. vivax isolates were collected. Scale bar = 100 kilometer.

Patients and sample collection

This study was reviewed and approved by the Ethics Committee of Sun Yat-sen University, China (IRB 2009A0101). Blood samples were collected from patients with symptomatic P. vivax malaria admitted to the local hospitals for antimalarial drug treatment in the city from June 2009 to October 2010. Inclusion criteria were patients infected with P. vivax, but not other human malaria species, as confirmed by peripheral smear examination. Preparation and staining of the blood slides were done following the procedures of WHO outlined in basic malaria microscopy. From patients who had the fever-chill cycle and visited a designated hospital for malaria diagnosis and treatment in Bozhou area, forty-five P. vivax-infected patients were diagnosed by Giemsa-stained thin and thick blood smears, in which only two out of the 45 P. vivax infected patients had a travel history outside of Anhui province, and none of the patients had a travel history abroad. After informed consent from all adults or legal guardians of children, 1.0 ml of whole blood sample was collected in EDTA and stored at -20°C until DNA extraction.

DNA extraction

Plasmodium vivax genomic DNA was extracted from 200 μl of each infected blood sample using DNA blood kit following the manufacturer’s instructions (TaKaRa, Japan) with minor modifications. The DNA was dissolved in TE buffer (10 mM Tris–HCl, pH 8.0, 0.1 M EDTA) and stored at -20°C until use. The quality of total DNA was analysed by running 5 μL of each DNA sample on a 1.0% agarose gel stained with ethidium bromide and visualized with UV illumination.

PCR amplification of pvmsp-3α, pvmsp-1, and pvcsp genes

To amplify the pvmsp-3α, pvmsp-1, and pvcsp genes, a nested PCR amplification method was used following previously reported protocols with some minor modifications [10, 13, 14]. Oligonucleotide primers and cycling conditions are listed in Additional file 1. All amplification reactions were carried out in a total volume of 25 μL containing 17.0 μl of dH2O, 1.0 μl of each primer (10 pM), 1.5 μl of MgCl2 (25 mM), 0.5 μl of Taq polymerase (5 U/μl), 0.5 μl of dNTP mixture (10 mM each), and 2.5 μl of 10 × PCR buffer, following the manufacturer’s instructions (BBI, Canada). Primary amplification reactions were initiated with the addition of 1.0 μL of template genomic DNA prepared from the blood samples, and 1.0 μL of the primary reaction amplification was used as template in the secondary amplification reactions. The amplified PCR products were resolved on 1.0% agarose gel, and the sizes of the PCR products were determined using a 100 or 500 bp DNA ladder (NewProbe, China).

PCR/RFLP analysis of pvmsp-3α gene

For RFLP analysis of pvmsp-3α gene, the PCR products were digested individually with restriction enzymes Hha I and/or Alu I in 20 μl reaction volumes at 37°C for 4 h as previously described [10]. Briefly, all digestion reactions were carried out in the presence of 10.0 μl PCR product, 7.0 μl of dH2O, 1.0 μl of enzymes Hha I or Alu I (5 U/μl), and 2.0 μl of buffer according to the manufacturer’s instructions (TaKaRa, Japan). After electrophoresis on 2.5% agarose gel, the enzyme-digested fragments were visualized under UV illumination. The sizes of the digested fragments were estimated using a 100 bp ladder of molecular weight markers. The results were recorded and analysed on a Gel Doc XR image analyzer using Quantity One software (Tanon, China).

Sequence analysis of pvmsp-1 and pvcsp genes

The nested PCR products of pvmsp-1 and pvcsp genes were directly sequenced in both directions using an ABI PRISM3730 DNA sequencer by Sangon Biotech (Shanghai, China). Nucleotide and amino acid sequences were aligned and compared using CLUSTAL W of BioEdit 7.0 program with the following published sequences: VK210 (accession no. M28746) and VK247 (accession no. M28745) of pvcsp and Sal-I (accession no. M75674) and Belem (accession no. M60807) of pvmsp-1, respectively. Sequence relationship trees of the pvmsp-1 and pvcsp genes from the Anhui isolates and published sequences from isolates of different geographic locations of Asia were constructed using neighbour-joining method implemented in MEGA 4 program [15]. Bootstrap proportions were used to assess the robustness of the tree with 1,000 bootstrap replications.

Results

PCR/RFLP analysis of the pvmsp-3α gene

The pvmsp-3α gene was successfully amplified from all the 45 P. vivax isolates examined. All of the PCR products from P. vivax isolates were the same size of approximately 1.9 kb (type A); no type B (~1.5 kb) or type C (~1.1 kb) was detected. To further characterize variation in the gene, the PCR products were digested with restriction enzymes Hha I and Alu I, and 7 (PH1 to PH7) and 6 (PA1 to PA6) distinct variants were detected, respectively, based on the restriction patterns of the PCR products (Figures 2A, B). In the Hha I digestion, allelic variants PH1, PH2, and PH3 were the most common patterns, with frequencies of 42.2%, 17.8%, and 17.8%, respectively (Figures 2C), representing 77.8% of the total samples tested. Similarly, 82.2% of the Alu I digested samples were allele variants of PA1, PA2, and PA5, representing 40.0%, 20.0%, and 22.2%, respectively (Figures 2C).

Figure 2
figure 2

PCR/RPLF genotypes and allelic frequencies of 45 P. vivax Anhui isolates based on pvmsp-3α . (A) The amplification products of pvmsp-3α observed in Anhui P. vivax isolates were digested by Hha I, and seven variants were obtained (PH1 to PH7). (B) The amplification products of pvmsp-3α observed in Anhui P. vivax isolates were digested by Alu I, and six variants were obtained (PA1 to PA6). The lane with the molecular weight marker (100 bp ladder) is labeled as M. (C) The allelic frequencies of pvmsp-3α after Hha I or Alu I restriction from Anhui P. vivax isolates. PH: The allelic type of pvmsp-3α was based on the fragment sizes of Hha I-digested PCR products; PA: The allelic type of pvmsp-3α was based on the fragment sizes of Alu I -digested PCR products.

DNA sequence polymorphism and clustering of pvmsp-1

All 45 isolates from Bozhou were successfully amplified for pvmsp-1 gene. The size of PCR products varied from 420 to 520 bp for the isolates examined, with the 420 bp product being the most common fragment (80.2%). Sequence analyses of the PCR products showed that the isolates could be divided into two types: Sal-I (S) type (24.4%) or a recombination (R) type (75.6%) that was characterized by a combination of a Sal-I-like sequence at the 5′-end and a Belem-like sequence at the 3′-end, including a polymorphic poly Q segment [8, 16]; no Belem type was found in these isolates. In total, seven distinct variants were found in the isolates, with five variant subtypes for the Sal-I type and two variant subtypes for the R-type, and with 13 and 10 poly Qs, respectively (Figure 3A). Of the seven sequence variants, R2-type was the most prevalent, accounting for 66.7% of all isolates examined. Sal-I type 5 variants were evenly distributed among the isolates, with frequencies of 6.7%, 4.4%, 4.4%, 4.4%, 4.4%, and 6.7% for Sal-I type S1, S2, S3, S4, S5, and R1 types, respectively (Figure 3B).

Figure 3
figure 3

Alignment of amino acid sequences of 7 pvmsp-1 distinct allelic variants and frequencies of pvmsp-1 allelic variants identified from 45 P. vivax Anhui isolates. (A) Sequences were constructed to the published amino acid sequences of Belem (M60807) and Sal-I (M75674). Dots and dashes represent identical residues and deletions, respectively. Amino acid changes resulting from nucleotide substitutions are shown in bold. S: Sal-I type; R: Recombination type. (B) Allelic frequencies of 7 pvmsp-1 allelic variants obtained from P. vivax Anhui isolates. S: Sal-I type; R: Recombination type of pvmsp-1.

A dendrogram of the pvmsp-1 variable block 5 sequences was constructed using nucleotide sequences from the Anhui isolates and 57 published pvmsp-1 sequences obtained from the isolates of different geographic locations of China (seven from Anhui province, seven from Fujian province, eight from Hainan province, six from Liaoning province, and five from Yunnan province of China) and other Asian countries (three from Afghanistan, two from Bangladesh, three from India, three from Iran, two from Myanmar, one from North Korea, three from South Korea, four from Papua New Guinea, three from Thailand, and two including the Sal-I strain and Belem strain from Brazil). This analysis grouped the P. vivax Anhui population into two major clusters (Sal-I type and R-type) with different subtypes in each group (Figure 4). The Sal-I type strains clustered with isolates from Hainan, Fujian, and Anhui of China as well as those from Afghanistan, Myanmar, Bangladesh, India, Thailand, and Papua New Guinea, whereas the R-type parasites clustered with parasites from Liaoning of China (northern China), North Korea, Iran, Afghanistan, Myanmar, and Papua New Guinea. Some of the P. vivax isolates examined in this study showed 100% identity with strains from other regions of China or even from different Asian countries, including Anhui province [AY465398 with S3 (Sal-I type)], Liaoning province [JQ606831 with S3, and JQ606834 with R2], Fujian province [AY538672 with R2], Hainan province [AY465377 with S2], Yunnan province [AY465391 with R2], South Korea [HQ171940 with S3, and HQ171937 with R2], North Korea [AF153032 with R2], India [AY229866 with S1], and Iran [AY502161 with R1]. Sal-I type S4 and S5 allele variants were new alleles identified in this study.

Figure 4
figure 4

Dendrogram of the pvmsp-1 gene based on the nucleotide acid sequences of 7 allelic variants from 45 P. vivax Anhui isolates and from 52 published sequences around the world in public databases. The tree was constructed using neighbor-joining method implemented in MEGA4 program. The geographical origin of the 52 pvmsp-1 published sequences is as follows: Anhui province (AY465395, AY 465398, AY465400, AY465401, AY465402, AY465405, and AY465406), Liaoning province (JQ319388, JQ582422, JQ390468, JQ606831, JQ606833, and JQ606834), Fujian province (AY538667, AY538668, AY538669, AY538670, AY538671, AY538672, and AY538673), Hainan province (AY465381, AY465377, AY465384, AY465378, AY465383, AY465385, AY465380, and AY465387), and Yunnan province (AY465389, AY465391, AY465392, AY465393, and AY465386) of China, Afghanistan (FJ600704 and FJ600714), Bangladesh (AF435617 and AF435618), India (AY229866 and FJ490829), Iran (FJ619327 and AF502161), Myanmar (EU048259 and EU048266), North Korea (AF216677), Papua New Guinea (DQ376102, DQ376105, and DQ376106), South Korea (HQ171940, JQ317282, and HQ171937), Thailand (GQ890955 and GQ890961), Belem (M60807), and Sal-I (M75674). The length of the line (bottom) is proportional to the genetic differences (%). Numbers on the branches indicate bootstrap proportions (1000 replicates). Only bootstrap values above 50% are displayed on the tree. S: Sal-I type; R: Recombination type.

DNA sequence polymorphism and clustering of pvcsp

To analyse polymorphisms in the pvcsp gene, the central repeat and its flanking regions of the gene from the samples examined were amplified and sequenced. All of the 45 isolates were successfully amplified for the pvcsp gene with PCR products having sizes ranging from 540 to 720 bp, with a 600 bp fragment being the predominant product among these isolates (58.5%). All of the DNA sequences belonged to VK210 type with 12 distinct variants; no VK247 and P. vivax-like types were detected within these isolates (Figure 5A) [1719]. All variants started with the same pre-repeat sequence (KLKQP Region). The isolates displayed variations in the central peptide repeat motifs GDRA (A/D) GQPA with alternations of the repeating units, ranging from 7 to 17 repeats among these isolates. Of the 12 variants, the most prevalent sequence variant was VK210-1 subtype (42.2%, 19/45), followed by subtypes VK210-6 (15.6%, 7/45) and VK210-7 (13.3%, 6/45). The remaining nine variants of pvcsp gene were evenly distributed among the isolates (Figure 5B).

Figure 5
figure 5

Sequence alignments of amino acid sequences of 12 pvcsp distinct allelic variants and frequencies of pvcsp allelic variants identified from 45 P. vivax Anhui isolates. (A) Sequences were aligned to the amino acid sequences of the reference strain VK210 (M28746). Dots and dashes represent identical residues and deletions, respectively. Amino acid changes resulting from nucleotide substitutions are shown in bold. (B)  Allelic frequencies of 12 pvcsp allelic variants obtained from P. vivax Anhui isolates.

A dendrogram of the pvcsp sequences was constructed using neighbor-joining method based on the nucleotide sequences from Anhui isolates in this study and 22 published pvcsp sequences from the isolates of different geographic locations of Asia (seven from China, one from Indonesia, two from India, two from Iran, two from Myanmar, two from North Korea, two from Papua New Guinea, two from the Solomon Islands, and two from South Korea) (Figure 6). The sequences clustered into two distinct groups for the VK210 and VK247 type isolates. And the analysis clearly showed that all subtypes of pvcsp from Anhui isolates clustered with the VK210 type in the tree. Some of the P. vivax Anhui isolates in this study showed 100% identity with other strains from China and Asian countries previously published, such as China’s Tibet (FJ601725 with VK210-6), North Korea (M20670 with VK210-6), and South Korea (DQ859770 with VK210-10); whereas the remaining VK210 subtypes were new alleles identified in this study.

Figure 6
figure 6

Dendrogram of the pvcsp gene constructed (neighbor-joining method using by MEGA4 program) based on the nucleotide acid sequences from 12 allelic variants of 45 P. vivax Anhui isolates and from 22 published sequences around the world. The geographical origin of the 22 pvcsp published sequences is as follows: China’s Tibet (FJ601725, FJ601755, FJ601744, and FJ601732), China (U08977, U08978, and U08979), India (JQ308526 and FJ491139), Indonesia (JQ308527), Iran (AY632261 and AY632311), Myanmar (EU048256 and AF316584), North Korea (M20670 and EU401928), Papua New Guinea (EU031819 and EU031828), the Solomon Islands (U08983 and U08982), South Korea (DQ859768 and AF236846), VK210 (M28746), and VK247 (M28745). The length of the line (bottom) is proportional to the genetic differences (%). Numbers on the branches indicate bootstrap proportions (1000 replicates). Only bootstrap values above 50% are displayed on the tree.

Haplotypes from combined pvmsp-3α, pvmsp-1, and pvcsp allelic variants

By analysing variants of the three genetic markers (pvmsp-3α, pvmsp-1, and pvcsp) in combination, a total of 16 P. vivax haplotypes were identified in the 45 Anhui isolates with haplotype A of PH1 and PA1 of pvmsp-3α, R-2 of pvmsp-1, and VK210-1 of pvcsp being the most common haplotypes (37.8%), followed by haplotype D of PH2 and PA2 of pvmsp-3α, R-2 of pvmsp-1, and VK210-6 of pvcsp (Table 1). These haplotypes also contained the most common allelic variants of the individual markers. PH1 and PA1 of pvmsp-3α, R-2 of pvmsp-1, and VK210-1 of pvcsp were detected with high frequency in their corresponding genetic markers: PH1 and PA1 of pvmsp-3α (42.2% and 40.0%, respectively), R-2 of pvmsp-1 (66.7%), and VK210-1 of pvcsp (42.2%).

Table 1 Combined haplotypes and frequencies of pvmsp-3α , pvmsp-1 , and pvcsp from 45 vivax isolates of Anhui province

Discussion

Malaria is one of the major parasitic diseases having a wide distribution in China with the prevalence gradually decreasing from south to north. Compared with southern (Hainan province) and southwestern (Yunnan province) regions, central China (including Anhui province) has a low endemicity of malaria, with P. vivax being the only species after the year 2000. Anhui province is the most seriously affected area in central China, with the largest number of malaria cases reported during the years 2005–2010. Different parasite genotypes are circulating in each endemic area of malaria-hypoendemic regions of Southeast Asia, and that geographic isolation may exist [7]. Additionally, with increasing movement of human populations, the P. vivax populations and transmission pattern in the regions are also changing constantly [20]. However, the genetic diversity of P. vivax parasites circulating in the area of central China is poorly understood. To better understand the P. vivax populations and transmission dynamics in Anhui of central China, the extent of genetic diversity of P. vivax parasites circulating in Bozhou of Anhui province were investigated using three polymorphic genetic marker pvmsp-3α, pvmsp-1, and pvcsp genes, and their relationships with parasites from other regions of China and Asian countries were analysed in this study. The data in the present study indicate that there was some degree of genetic diversity of P. vivax in Anhui province, and there appeared to be signs of genetic exchange and genetic recombination among vivax parasites in the region. Although a relatively small number of samples (45) could be obtained, the patients came from many parts of the city. Additionally, the diverse genetic backgrounds from the three genetic markers suggested that the parasites were not from outbreaks of 1–2 parasite strains and should be able to represent the genetic backgrounds of the parasite populations in the region.

It has been proposed that pvmsp-3α is a reliable molecular marker for analysis of P. vivax populations [11]. Based on the length of PCR products, four different allelic types of pvmsp-3α have been characterized around the world: type A (~1.9 Kb), type B (~1.5 Kb), type C (~1.1 Kb), and type D (~0.3 Kb) [6, 21]. In the present study, only type A (~1.9 Kb) was identified among the 45 Anhui isolates tested. The observation was different from a previous report from several locations in China [7], in which type A, B, and C were detected in P. vivax isolates from Suzhou city of Anhui province, Sanya city of Hainan province, and Dehong prefecture of Yunnan province. The lack of other pvmsp-3α allele types in Bozhou city of Anhui province could simply be due to the relatively small sample size of this study. The results in some degree are consistent with the reports of type A being the most prevalent type around the world, with a frequency 70% to 100% (average ~80%) in many regions of China and countries in Asia and South America [7, 2127]. Seven patterns were detected in the pvmsp-3α gene after digestion of the PCR products with HhaI, with the most frequent allele variant being PH1 subtype. The results were consistent with those reported from other regions of the world, including those from Myanmar (n = 14) [7], Iran (n = 6) [14], Colombia (n = 9) [22], French Guiana (n = 11) [23], Thailand (n = 6) [28], Venezuela (n = 9) [29], and Peru (n = 17) [30]. However, outside Anhui province, the allele PH3, the second most common one in the present study, was the most prevalent in other areas of China including Hainan and Yunnan provinces, and also in Myanmar [7]. All of the Hha I allelic types found in the present study have been reported in other regions of the world, including Myanmar (PH1, PH3, PH4, PH5, and PH6) [7], India Chennai (PH1, PH4, PH5, and PH6) [11, 21], Colombia (PH1, PH2, PH3, PH4, PH5, and PH6) [22], Brazil (PH3 and PH4) [26], and Thailand (PH1, PH3, PH4, PH6, and PH7) [28]. In China, five of the seven alleles were reported previously, e.g., PH1, PH3, PH4, PH5, and PH6 from Suzhou city of Anhui province; PH3, PH4, PH5, and PH6 from Hainan province; and PH1, PH3, and PH4 from Yunnan province [7].

Six patterns were detected in the pvmsp-3α gene after digestion of the PCR products with AluI in the present study, with the most predominated PA1 pattern in Bozhou city. This observation has not been reported in other areas of the world. Some allele types of Alu I digestion were comparable to those reported in other areas of the world, such as Colombia (PA2, PA4, PA5, and PA6), India (PA1, PA2, PA4, PA5, and PA6), and Iran (PA1 and PA2), suggesting that these allele types of P. vivax may have a global distribution [14, 21, 22]. The allele variant PA3 was a new allele identified in this study. The analyses from pvmsp-3α suggest that P. vivax populations in Bozhou city shared the majority of allelic variants with other parts of China and the world; however, some new allelic types were emerging.

To obtain more information of genetic polymorphism of the vivax malaria samples, the DNA sequences encoding the pvmsp-1 and pvcsp genes were analysed. Three allele-types, Belem, Sal-I, and recombination types of the variable block 5 have been described from P. vivax isolates worldwide [8, 16]. Genetic studies of P. vivax malaria parasites using the pvmsp-1 marker have been reported from many countries, including India, Myanmar, Pakistan, and Korea [15, 25, 31, 32]. In the present study, only Sal-I and R-types were detected among the Anhui isolates with R-type being the predominant allele types (75.6%). It has been reported that the R-type result from intragenic recombination of Belem and Sal-I types in mosquito vector [33, 34]. The results indicate that inter-allelic recombination can be an important source of pvmsp-1 genetic diversity among Anhui province locations. The high level of inter-allelic recombinant types in this isolate was in concordance with previous reports in Pakistan [25], India [32], Colombia [32], Thailand [35], and Brazil [36]. Although the Belem type is the predominant allelic type in other P. vivax malaria endemic regions in the world, including northern Iran [14], south-eastern and south-western parts of Afghanistan [37], and Azerbaijan [38], Belem type was not found among the Anhui isolates in this study. The absence of Belem type is intriguing, considering the presence of the recombinant types, which could be related to population changes of mosquito species or sampling bias. It has been reported that the distributions of different P. vivax populations in southern Mexico were largely determined by their infectivity to two species of anopheline vectors in the areas [39]. A similar situation could occur in China too; changes in mosquito strains or species could favor a subpopulation of vivax parasites, leading to changes in parasite population structure.

For pvcsp gene, three sequence types of VK210, VK247, and P. vivax-like have been detected based on a central repeat domain that varies in sequences and number of repeat units [1719, 40, 41]. VK210 being the predominant type has been observed in Myanmar (98.3%) [15], Iran (69.3%) [24, 42], Thailand (77%) [28], Afghanistan (86.6%) [37], Pakistan (95.7%) [37], Azerbaijan (100%) [38], Brazil (86%) [40], Guyana (92%) [41], and India (99.3%) [43], although VK247 was also reported to be the predominant type previously [4446]. The VK247 and P. vivax-like types were not detected in the Anhui isolates in the present study. It has been proposed that the geographical differences in distribution of different parasite types may be related to the susceptibility of different species of Anopheline vectors to the infections by different malaria parasites in endemic areas [47, 48]. The main malaria vectors in Anhui province are An. sinensis and An. anthropophagus, but the density of An. anthropophagus is very low in recent years, and the re-emergence of P. vivax malaria was mainly caused by An. sinensis[49]. Therefore, the prevalence of VK210 type may be correlated with the increase of An. sinensis malaria vector in the region. A total of 12 VK210-1 subtypes were identified in 45 blood samples from Anhui province, two subtypes of which, VK210-6 and VK210-10, were also detected in other endemic areas, including China’s Tibet, North Korea, and South Korea through phylogenetic and Blast analysis.

In summary, P. vivax populations in Bozhou city of Anhui province have some degree of genetic diversity in three highly polymorphic markers, with 16 unique haplotypes observed among 45 samples (~36%). However, the present results also showed that P. vivax populations in Bozhou city of Anhui province had relative even major genotypes, pvmsp-3α type A (100%), pvmsp-1 R2 (75.6%), and pvcsp VK210 subtype (100%), suggesting its relative uniform parasite populations. It is interesting that pvmsp-1 recombinant type was detected, which may be related to the recent changes in vector population (increase in An. sinensis mosquito) in the region, or may have been simply sampled some subpopulations from recent outbreaks. The data in this study showed shared genotypes with those from Korea and other countries; there is a possibility that these parasites have the same origin and were disseminated by travelers. Although limited samples were included in this study, and the data may not represent the whole picture of genome diversity of the parasite populations in the region, the results suggest parasite populations with relatively diverse genetic backgrounds in the region. The genetic diversity information obtained may provide valuable information for tracking and monitoring future malaria infections in this region.

References

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

    Article  CAS  PubMed  Google Scholar 

  2. Mendis K, Sina BJ, Marchesini P, Carter R: The neglected burden of Plasmodium vivax malaria. Am J Trop Med Hyg. 2001, 64: 97-106.

    CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  PubMed  Google Scholar 

  4. Zhang W, Wang L, Fang L, Ma J, Xu Y, Jiang J, Hui F, Wang J, Liang S, Yang H, Cao W: Spatial analysis of malaria in Anhui province, China. Malar J. 2008, 7: 206-10.1186/1475-2875-7-206.

    Article  PubMed Central  PubMed  Google Scholar 

  5. Zhu G, Xia H, Zhou H, Li J, Lu F, Liu Y, Cao J, Gao Q, Sattabongkot J: Susceptibility of Anopheles sinensis to Plasmodium vivax in malarial outbreak areas of central China. Parasit Vectors. 2013, 6: 176-10.1186/1756-3305-6-176.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Cui L, Escalante AA, Imwong M, Snounou G: The genetic diversity of Plasmodium vivax populations. Trends Parasitol. 2003, 19: 220-226. 10.1016/S1471-4922(03)00085-0.

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Premawansa S, Snewin VA, Khoury E, Mendis KN, David PH: Plasmodium vivax: recombination between potential allelic types of the merozoite surface protein MSP-1 in parasites isolated from patient. Exp Parasitol. 1993, 76: 192-199. 10.1006/expr.1993.1022.

    Article  CAS  PubMed  Google Scholar 

  9. Qari SH, Collins WE, Lobel HO, Taylor F, Lal AA: A study of polymorphism in the circumsporozoite protein of huaman malaria parasites. Am J Trop Med Hyg. 1994, 501: 45-51.

    Google Scholar 

  10. Bruce MC, Galinski MR, Barnwell JW, Snounou G, Day KP: Polymorphism at the merozoite surface protein-3 alpha locus of Plasmodium vivax: global and local diversity. Am J Trop Med Hyg. 1999, 614: 518-525.

    Google Scholar 

  11. Zakeri S, Barjesteh H, Djadid ND: Merozoite surface protein-3 is a reliable marker for population genetic analysis of Plasmodium vivax. Malar J. 2006, 5: 53-10.1186/1475-2875-5-53.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Ding S, Ye R, Zhang D, Sun X, Zhou H, McCutchan TF, Pan W: Anti-folate combination therapies and their effect on the development of drug resistance in Plasmodium vivax. Sci Rep. 2013, 3: 1008-

    PubMed Central  PubMed  Google Scholar 

  13. Imwong M, Pukrittayakamee S, Gruner AC, Renia L, Letourneur F, Looareesuwan S, White NJ, Snounou G: Practical PCR genotyping protocols for Plasmodium vivax using Pvcs and Pvmsp1. Malar J. 2005, 4: 20-10.1186/1475-2875-4-20.

    Article  PubMed Central  PubMed  Google Scholar 

  14. Zakeri S, Mehrizi AA, Mamaghani S, Noorizadeh S, Snounou G, Djadid ND: Population structure analysis of Plasmodium vivax in areas of Iran with different malaria endemicity. Am J Trop Med Hyg. 2006, 74: 394-400.

    CAS  PubMed  Google Scholar 

  15. Moon SU, Lee HW, Kim JY, Na BK, Cho SH, Lin K, Sohn WM, Kim TS: High frequency of genetic diversity of Plasmodium vivax field isolates in Myanmar. Acta Trop. 2009, 109: 30-36. 10.1016/j.actatropica.2008.09.006.

    Article  CAS  PubMed  Google Scholar 

  16. Del Portillo HA, Longacre S, Khouri E, David PH: Primary structure of the merozoite surface antigen 1 of Plasmodium vivax reveals sequences conserved between different Plasmodium species. Proc Natl Acad Sci USA. 1991, 88: 4030-4034. 10.1073/pnas.88.9.4030.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Rosenberg R, Wirtz RA, Lanar DE, Sattabongkot J, Hall T, Waters AP, Prasittisuk C: Circumsporozoite protein heterogeneity in the human malaria parasite Plasmodium vivax. Science. 1989, 245: 973-976. 10.1126/science.2672336.

    Article  CAS  PubMed  Google Scholar 

  18. Qari SH, Shi YP, Goldman IF, Udhayakumar V, Alpers MP, Collins WE, Lal AA: Identification of Plasmodium vivax-like human malaria parasite. Lancet. 1993, 341: 780-783. 10.1016/0140-6736(93)90559-Y.

    Article  CAS  PubMed  Google Scholar 

  19. Qari SH, Shi YP, Povoa MM, Alpers MP, Deloron P, Murphy GS, Harjosuwarno S, Lal AA: Global occurrence of Plasmodium vivax-like human malaria parasite. J Infect Dis. 1993, 168: 1485-1489. 10.1093/infdis/168.6.1485.

    Article  CAS  PubMed  Google Scholar 

  20. Pindolia DK, Garcia AJ, Huang Z, Smith DL, Alegana VA, Noor AM, Snow RW, Tatem AJ: The demographics of human and malaria movement and migration patterns in East Africa. Malar J. 2013, 12: 397-10.1186/1475-2875-12-397.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Prajapati SK, Joshi H, Valecha N: Plasmodium vivax merozoite surface protein-3α: a high-resolution marker for genetic diversity studies. J Vector Borne Dis. 2010, 47: 85-90.

    PubMed  Google Scholar 

  22. Cristiano FA, Perez MA, Nicholls RS, Guerra AP: Polymorphism in the Plasmodium vivax msp 3α gene in field samples from Tierralta, Colombia. Mem Inst Oswaldo Cruz. 2008, 103: 493-496. 10.1590/S0074-02762008000500015.

    Article  CAS  PubMed  Google Scholar 

  23. Veron V, Legrand E, Yrinesi J, Volney B, Simon S, Carme B: Genetic diversity of msp3α and msp1_b5 markers of Plasmodium vivax in French Guiana. Malar J. 2009, 8: 40-10.1186/1475-2875-8-40.

    Article  PubMed Central  PubMed  Google Scholar 

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

    Article  PubMed Central  PubMed  Google Scholar 

  25. 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-3α and pvcsp genes as molecular markers. Parasitol Int. 2010, 591: 15-21.

    Article  Google Scholar 

  26. Ribeiro RS, Ladeira L, Rezende AM, Fontes CJF, Carvalho LH, Alves de Brito CF: Analysis of the genetic variability of PvMSP-3α among Plasmodium vivax in Brazilian field isolates. Mem Inst Oswaldo Cruz. 2011, 106: 27-33.

    Article  CAS  PubMed  Google Scholar 

  27. Rungsihirunrat K, Chaijaroenkul W, Siripoon N, Seugorn A, Na-Bangchang K: Geneotyping of polymorphic markar (MSP3α and MSP3β) genes of Plasmodium vivax field isolates from malaria endemic of Thailand. Trop Med Int Health. 2011, 16: 794-801. 10.1111/j.1365-3156.2011.02771.x.

    Article  PubMed  Google Scholar 

  28. Cui LW, Mascorro CN, Fan Q, Rzomp KA, Khuntirat B, Zhou G, Chen H, Yan GY, Sattabongkot J: Genetic diversity and multiple infections of Plasmodium vivax malaria in western Thailand. Am J Trop Med Hyg. 2003, 685: 613-619.

    Google Scholar 

  29. Ord R, Polley S, Tami A, Sutherland CJ: High sequence diversity and evidence of balancing selection in the Pvmsp3α gene of Plasmodium vivax in the Venezuelan Amazon. Mol Biochem Parasitol. 2005, 144: 86-93. 10.1016/j.molbiopara.2005.08.005.

    Article  CAS  PubMed  Google Scholar 

  30. Sutton PL, Neyra V, Hernandez JN, Branch OH: Plasmodium falciparum and Plasmodium vivax infections in the Peruvian Amazon: propagation of complex, multiple allele-type infections without super-infection. Am J Trop Med Hyg. 2009, 816: 950-960.

    Article  Google Scholar 

  31. Lim CS, Kim SH, Kwon SI, Song JK, Lee KN: Analysis of Plasmodium vivax merozoite surface protein-1 gene sequence from resurgent Korean isolates. Am J Trop Med Hyg. 2000, 62: 261-265.

    CAS  PubMed  Google Scholar 

  32. Maestre A, Sunil S, Ahmad G, Mohmmed A, Echeverri M, Corredor M, Blair S, Chauhan VS, Malhotra P: Inter-allelic recombination in the Plasmodium vivax merozoite surface protein 1 gene among Indian and Colombian isolates. Malar J. 2004, 3: 4-10.1186/1475-2875-3-4.

    Article  PubMed Central  PubMed  Google Scholar 

  33. Gutierrez A, Vicini J, Patarroyo ME, Murillo LA, Patarroyo MA: Plasmodium vivax: polymorphism in the merozoite surface protein 1 gene from wild Colombian isolates. Exp Parasitol. 2000, 95: 215-219. 10.1006/expr.2000.4534.

    Article  CAS  PubMed  Google Scholar 

  34. Putaporntip C, Jongwutiwes S, Sakihama N, Ferreira M, Kho W, Kaneko A, Kanbara H, Hattori T, Tanabe K: Mosaic organization and heterogeneity in frequency of allelic recombination of the Plasmodium vivax merozoite surface protein-1 locus. Proc Natl Acad Sci USA. 2002, 99: 16348-16353. 10.1073/pnas.252348999.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Putaporntip C, Jongwutiwes S, Tanabe K, Thaithong S: Interallelic recombination in the merozoite surface protein (MSP-1) gene of Plasmodium vivax from Thai isolates. Mol Biochem Parasitol. 1997, 84: 49-56. 10.1016/S0166-6851(96)02786-7.

    Article  CAS  PubMed  Google Scholar 

  36. Santos-Ciminera PD, Alecrim MG, Roberts DR, Quinnan GVJ: Molecular epidemiology of Plasmodium vivax in the state of Amazonas, Brazil. Acta Trop. 2007, 102: 38-46. 10.1016/j.actatropica.2007.02.013.

    Article  CAS  PubMed  Google Scholar 

  37. 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. 2010, 11: 12-19.

    Article  Google Scholar 

  38. Leclerc MC, Menegon M, Cligny A, Noyer JL, Mammadov S, Aliyev N, Gasimov E, Majori G, Severini C: Genetic diversity of Plasmodium vivax isolates from Azerbaijan. Malar J. 2004, 3: 40-10.1186/1475-2875-3-40.

    Article  PubMed Central  PubMed  Google Scholar 

  39. Joy DA, Gonzalez-Ceron L, Carlton JM, Gueye A, Fay M, McCutchan TF, Su XZ: Local adaptation and vector-mediated population structure in Plasmodium vivax malaria. Mol Biol Evol. 2008, 25: 1245-1252. 10.1093/molbev/msn073.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Machado RL, Povoa MM: Distribution of Plasmodium vivax variants (VK210, VK247 and P. vivax-like) in three endemic areas of the Amazon region of Brazil and their correlation with chloroquine treatment. Trans R Soc Trop Med Hyg. 2000, 94: 377-381. 10.1016/S0035-9203(00)90110-X.

    Article  CAS  PubMed  Google Scholar 

  41. Bonilla JA, Validum L, Cumming R, Palmer CJ: Genetic diversity of Plasmodium vivax Pvcsp and Pvmsp1 in Guyana, South America. Am J Trop Med Hyg. 2006, 75: 830-835.

    CAS  PubMed  Google Scholar 

  42. Zakeri S, Mehrizi AA, Mamaghani S, Djadid ND, Snounou G: Circumsporozoite protein gene diversity amony temperate and tropical Plasmodium vivax isolates from Iran. Trop Med Int Health. 2006, 11: 729-737. 10.1111/j.1365-3156.2006.01613.x.

    Article  CAS  PubMed  Google Scholar 

  43. Kim JR, Imwong M, Nandy A, Chotivanich K, Nontprasert A, Tonomsing N, Maji A, Addy M, Day NPJ, White NJ, Pukrittayakamee S: Genetic diversity of Plasmodium vivax in Kolkata, India. Malar J. 2006, 5: 71-10.1186/1475-2875-5-71.

    Article  PubMed Central  PubMed  Google Scholar 

  44. Burket TR, Wirtz RA, Paru R, Garner P, Alpers MP: The population dynamics in mosquitoes and humans of two Plasmodium vivax polymorphs distinguished by different circumsporozoite protein repeat regions. Am J Trop Med Hyg. 1992, 47: 778-786.

    Google Scholar 

  45. Kain KC, Brown AE, Webster HK, Wirtz RA, Keystone JS, Rodrigues MH, Kinahan J, Rowland M, Lanar DE: Circumsporozoite genotyping of global isolates of Plasmodium vivax from dried blood specimens. J Clin Microbiol. 1992, 30: 1863-1866.

    PubMed Central  CAS  PubMed  Google Scholar 

  46. Gonzalez JM, Hurtado S, Arevalo-Herrera M, Herrera S: Variants of the Plasmodium vivax circumsporozoite protein (VK210 and VK247) in Colombian isolates. Mem Inst Oswaldo Cruz. 2001, 96: 709-712. 10.1590/S0074-02762001000500023.

    Article  CAS  PubMed  Google Scholar 

  47. Gonzalez-Ceron L, Rodriguez MH, Nettel JC, Villarreal C, Kain KC, Hernandez JE: Differential susceptibilities of Anopheles albimanus and Anopheles psedopunctipennis to infections with coindigenous Plasmodium vivax variants VK210 and VK247 in southern Mexico. Infect Immun. 1999, 67: 410-420.

    PubMed Central  CAS  PubMed  Google Scholar 

  48. Rodriguez MH, Gonzalez-Ceron L, Hernandez JE, Nettel JA, Villarreal C, Kain KC, Wirtz RA: Different prevalences of Plasmodium vivax phenotypes VK210 and VK247 associated with the distribution of Anopheles albimanus and Anopheles pseudopunctipennis in Mexico. Am J Trop Med Hyg. 2000, 62: 122-127.

    CAS  PubMed  Google Scholar 

  49. Zhou S, Zhang S, Wang J, Zhang X, Huang F, Li W, Xu X, Zhang H: Spatial correlation between malaria cases and water-bodies in Anopheles sinensis dominated areas of Huang-Huai plain, China. Parasit Vectors. 2012, 5: 106-10.1186/1756-3305-5-106.

    Article  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

Research reported in this publication was supported in part by the Fogarty International Center of the National Institutes of Health under Award Number R01TW008151, the Divisions of Intramural Research at the National Institute of Allergy and Infectious Diseases, National Institutes of Health, and the “111 Project” sponsored by the State Bureau of Foreign Experts and Ministry of Education of China (no. B12003). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. We thank Cindy Clark at NIH Library Writing Center for editing.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Fangli Lu.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

XZS and FL designed and supervised the study, analysed the data, and wrote and revised the manuscript. BH and SH carried out the experiments, and analysed the data, and wrote the manuscript. HG, YX, FX, XH, YY, and SW participated in field work and preliminary data analysis. All authors read and approved the final manuscript.

Bo Huang, Shiguang Huang contributed equally to this work.

Electronic supplementary material

12936_2013_3079_MOESM1_ESM.doc

Additional file 1: Sequences of the primers and cycling conditions used to amplify the pvmsp-3α , pvmsp-1 , and pvcsp genes of P. vivax isolates from Anhui province. The data provided oligonucleotide primers and cycling conditions of pvmsp-3α, pvmsp-1, and pvcsp genes. (DOC 40 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Huang, B., Huang, S., Su, Xz. et al. Genetic diversity of Plasmodium vivax population in Anhui province of China. Malar J 13, 13 (2014). https://doi.org/10.1186/1475-2875-13-13

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-2875-13-13

Keywords