- Open Access
Genetic diversity of transmission-blocking vaccine candidates Pvs25 and Pvs28 in Plasmodium vivax isolates from Yunnan Province, China
Parasites & Vectorsvolume 4, Article number: 224 (2011)
Transmission-blocking vaccines (TBVs) have been considered an important strategy for disrupting the malaria transmission cycle, especially for Plasmodium vivax malaria, which undergoes gametocytogenesis earlier during infection. Pvs25 and Pvs28 are transmission-blocking vaccine candidates for P. vivax malaria. Assessment of genetic diversity of the vaccine candidates will provide necessary information for predicting the performance of vaccines, which will guide us during the development of malaria vaccines.
We sequenced the coding regions of pvs25 and pvs28 from 30 P. vivax isolates from Yunnan Province, identifying five amino acid haplotypes of Pvs25 and seven amino acid haplotypes of Pvs28. Among a total of four mutant residues, the predominant haplotype of Pvs25 only had the I130T substitution. For Pvs28, a total of eight amino acid substitutions were identified. The predominant haplotype of Pvs28 had two substitution at positions 52 (M52L) and 140 (T140S) with 5-6 GSGGE/D tandem repeats at the end of fourth EGF-like domain. Most amino acid substitutions were common with previous reports from South Asian isolates. Although the nucleotide diversity of pvs28 (π = 0.0034 ± 0.0012) was significantly higher than pvs25 (π = 0.0013 ± 0.0009), it was still conserved when compared with the blood stage vaccine candidates.
Genetic analysis revealed limited genetic diversity of pvs25 and pvs28, suggesting antigenic diversity may not be a particular problem for Sal I based TBVs in most P. vivax-endemic areas of China.
Plasmodium vivax, the most widespread species of human malaria parasites, is responsible for the majority of malaria cases outside of Africa, and the most prevalent form of relapsing malaria. Although not as lethal as P. falciparum malaria, P. vivax has caused substantial morbidity for human populations residing in Asia and South America [1, 2]. Throughout the malaria control history, P. vivax malaria has displayed tremendous resilience to control efforts, which is in part due to two biological characteristics of this parasite. One is the formation of hypnozoites in the liver, which is responsible for relapses of the disease. The other is that P. vivax undergoes gametocytogenesis before manifestation of the disease symptoms, making transmission possible before treatment.
P. vivax is the predominant Plasmodium species in China, and in recent years, P. vivax cases accounted for more than 90% of all malaria cases . As The Ministry of Health of China has set the goal of malaria elimination by 2020 , interruption of vivax malaria transmission is a major challenge. Yunnan is one of the two provinces that have year-round local transmission of P. vivax and P. falciparum. It is located in southwest China and borders Myanmar to the west and Laos and Vietnam to the south . A total of 250,070 confirmed cases of P. vivax malaria and 44,465 cases of P. falciparum malaria have been reported between 1991 and 2006 in Yunnan . In recent years, we have witnessed a significant decline of malaria cases in Yunnan . A total of 32,566 confirmed cases of P. vivax malaria and 5,821 cases of P. falciparum malaria have been reported between 2006 and 2009 in Yunnan [3, 7–9]. Most of the malaria cases are now clustered in counties located at the China-Myanmar border area [3, 10]. Despite improvement in malaria situation in this region, malaria control along the borders is problematic, since reintroduction of malaria by human migration is still difficult to monitor. Thus, it is particularly important to develop effective strategies to control and eliminate malaria in this region.
Transmission-blocking vaccines (TBVs), which target the sexual stages of malaria parasites to prevent their further development within mosquitoes, are considered an important strategy for disrupting malaria transmission [11, 12]. To date, several sexual stage antigens have been characterized, which showed excellent transmission blocking activities. These include the pre-fertilization antigens P48/45 and post-fertilization antigens P25 and P28. P25 and P28, are specifically expressed on Plasmodium zygote and ookinete surface , and they are essential for the survival of ookinetes in the mosquito midgut, and subsequent penetration of the midgut epithelium and transformation into oocysts . Among the TBVs candidates, Pvs25 has received the most attention. Recombinant Pvs25 protein expressed in bacteria , yeasts [15, 16] and baculovirus  exhibits strong immunogenicity and immune sera display significant transmission-blocking activity on the development of sporozoites  and oocysts [15–17]. The immunogenicity trials of recombinant Pvs25 have been undertaken in rhesus monkeys  and human volunteers [19, 20]. A phase I clinical trial of Pvs25H vaccine in human volunteers has demonstrated that vaccine-induced antibodies have significant transmission-blocking activity . It has been reported that antisera to recombinant Pvs25 and Pvs28 based on the Sal-I strain recognized the corresponding molecules expressed by field-isolated parasites in Thailand, and that these antisera blocked transmission of field isolates .
Polymorphisms represent a major impediment in vaccine design. An effective vaccine must either be based on conserved regions or incorporate multiple allelic forms of the antigen. As antigenic polymorphisms limit the immunogenicity and immunoreactivity of vaccines, it is necessary to understand antigenic variations in order to generate effective vaccines against natural Plasmodium infections . Pvs25 and Pvs28 contained conserved structures, which are characterized by a secretory N-terminal signal sequence followed by four epidermal growth factor (EGF) domains and a glycosylphosphatidylinositol (GPI) anchor . To date, genetic diversity of pvs25 and pvs28 genes has been surveyed in several Asian countries, including South Korea , Indonesia , Iran , India [25, 26], Bangladesh , and Thailand . In this study, we aim to analyze the genetic diversity of these candidate TBVs genes in parasite isolates from Yunnan Province, China.
Yunnan Province is one of the most highly endemic regions in China with year-round transmission of P. vivax and P. falciparum. Since 2001, there has been a strong seasonal pattern of malaria incidence characterized by a peak of infection occurring from June to July and another peak from October to November . Counties with higher malaria incidence rates are clustered in the border area and the Yuanjiang River Basin [6, 27]. The predominant malaria vectors are Anopheles minimus and Anopheles dirus .
A total of 30 parasite samples were collected in 2004 from three regions in Yunnan. All volunteers in this study were symptomatic patients diagnosed with P. vivax malaria using Giemsa stained thin smear examination by malaria clinic staff. After informed consent/assent was obtained from either adults or parents or legal guardians of children, ~100 μl of finger-prick blood were collected from each patient on filter papers. The patients were Yunnan local residents, 13 to 50 years old, living in Baoshan (10 individuals), Ruili (10 individuals) and Longchuan (10 individuals) (showed in Figure 1). Patients were then treated according to the Ministry of Health drug policy. This study protocol was approved by the Biomedical Research Ethics Review Board at China Medical University.
Parasite DNA preparation and amplification
Three circles, 2.5 mm in diameter, were punched out from each dried blood spot and Plasmodium DNA was purified using the QIAamp DNA Mini kit (QIAGEN, Hilden, Germany) according to the manufacturer's protocol. DNA was eluted in 60 μl of H2O. pvs25 (GenBank: AF083502.1) and pvs28 (GenBank: AF083503.2) of the Salvador I (Sal I) strain were used as reference sequences. The full-length open reading frames (ORFs) of pvs25 and pvs28 were amplified using the following primer pairs: pvs25 F (5'-CACTTAGCCAAAATGAACTC-3') and pvs25 R (5'-AAAGGACAAGCAGGATGATA-3') for pvs25; pvs28 F (5'-CTACCACAGCTTGCTGTTCC-3') and pvs28 R (5'-TGACATCATGAAGAAGGCG-3') for pvs28. The PCR reaction contained 2 μl of 10 × KOD-Plus buffer, 2 μl of 2 mM dNTPs, 0.8 μl of 25 mM MgSO4, 1.0 μl of 10 μM of each primer, 0.2 units of KOD-Plus DNA polymerase (Toyobo, Osaka, Japan), and 0.5 μl genomic DNA template in a final volume of 20 μl. Amplification conditions were as follows: initial denaturing at 94°C for 2 min, 35 cycles of 94°C for 15 sec, 56°C (for pvs28) or 54°C (for pvs25) for 15 sec, and 68°C for 1 min, and a final extension at 68°C for 5 min. PCR products were cleaned up using Exonuclease I and Shrip alkaline phosphatase treatment according to the manufacturer's instructions (USB, CA, USA).
Pvs25 and pvs28 ORFs were sequenced using the ABI Prism® BigDye™ Terminator cycle sequencing kit (Applied Biosystems, Foster City, CA, USA) on an ABI PRISM® 310 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). Sequencing primers are: pvs25 F (5'-GAAACCCTAGGCAAAGCATG-3') and pvs25 R (5'- GGGACTTTGCCA ATAGCACA-3') for pvs25, and pvs28 F (5'-AACTGTGGAGACTACGCTG-3') and pvs28 R (5'-ATATTACAAGAGCACATGGTG-3') for pvs28.
The ORFs of pvs25 and pvs28 were aligned with the reference sequences using the CLUSTAL W program  with manual corrections. Nucleotide diversity (π) and its standard error (SE) were calculated with the Jukes and Cantor method using MEGA 4.0. The number of haplotypes (h), haplotypes diversity (Hd), segregating sites (S) were computed using the options available in DnaSP 5.0 .
Sequence polymorphisms of the pvs25 and pvs28 genes
Pvs25 and pvs28 genes were successfully amplified in 30 P. vivax samples obtained from Yunnan. The sequence of pvs25 (645 bp) and pvs28 (712 bp) were determined by direct sequencing of the PCR products. The 30 pvs25 sequences contained four polymorphic nucleotide sites (C103A, G289C, T389C, and C391A), which resulted in four amino acid substitutions (L35M, E97Q, I130T, and Q131K). These substitutions resulted in a total of five amino acid haplotypes, hitherto referred to as Pvs25-I, II, III, IV and V (Table 1). The predominant amino acid haplotype of Pvs25 only had the I130T substitution.
The 30 pvs28 sequences had nucleotide polymorphisms at 13 sites, which included four synonymous substitutions (T346C, T558C, C573T, and G579A) and eight nonsynonymous substitutions (G40T, A154C, T292A, G313A, T346G, C419G, G536A, and C666G), which resulted in eight amino acid substitutions (V14L, M52L, L98I, E105K, L116V, T140S, G179E, and I222M). A total of seven amino acid haplotypes were identified: Pvs28-I, II, III, IV, V, VI and VII (Table 2). The predominant amino acid haplotype of Pvs28 had two amino acid substitutions at positions 52 (M52L) and 140 (T140S). Compared with Sal-I strain, we found the majority of clones had six copies of GSGGE/D tandem repeats at the end of fourth EGF-like domain. Only one clone with Pvs28-I haplotype had five copies of repeat regions.
Compared with the worldwide isolates
Compared with previous reports of worldwide isolates, Pvs25-I haplotype showed 100% identity with isolates from Vietnam (GenBank: ABG29072.1), North Korea (GenBank: AAV33640.1) , Thailand (GenBank: AB091730.1)  and Mexico (GenBank: ABS70906.1-ABS70935.1) . The Pvs25-II haplotype was similar to some isolates from Bangladesh (GenBank: BAA94348.1-BAA94350.1)  and Thailand (GenBank: BAC66001.1) . The Pvs25-III haplotype was common with isolates from Iran (GenBank: ACJ54129.1, ACJ54128.1) , Thailand (GenBank: BAC66003.1) , Indonesia (GenBank: AAV33639.1) , India (GenBank: BAA94346.1)  and Bangladesh (GenBank: BAA94347.1) . Only Pvs25-IV was unique in China. Like in some South Korean isolates, position 35 was also polymorphic, but the substitution is L35M in the Yunnan isolate as compared to L35P in South Korean isolates (GenBank: ADJ18839.1) .
The Pvs28-I haplotype in Yunnan isolates was also common in isolates from Iran (GenBank: ACJ54136.1)  and Bangladesh (GenBank: BAA94352.1, BAA94353.1, BAA94371.1, BAA94362, BAA94363) . The Pvs28-II haplotype was similar with those from Thailand (GenBank: BAC66011.1)  and Bangladesh (GenBank: BAA94370.1, BAA94369.1) . Parasite isolates from Yunnan contained some region-specific haplotypes. Haplotype Pvs28-III with a new V14L substitution was found only in one isolate. The remaining haplotypes (Pvs28-IV, V and VI) were only identified in the Yunnan isolates. The G179E mutation was only found in South Korean (GenBank: ADJ18881.1)  and Chinese isolates.
The Pvs25 and Pvs28 amino acid substitutions identified from worldwide isolates including those from China are summarized in Tables 3 and 4. Among a total of 17 variant amino acids in Pvs25 from worldwide isolates, I130T was the most common amino acid substitution from American and Asian isolates. The Q87K substitution was common in American isolates such as those from Brazil  and Mexico , and Western Asian isolates such as those from Iran  and Turkey, but not detected in South and East Asian isolates. Amino acid residue E97Q was detected in some Asian isolates such as those from Iran , Bangladesh , Thailand , India [25, 26], Indonesia , South Korea , Central China  and Yunnan, but rare in American isolates. Compared with previous studies from Central China isolates, amino acid residue I130T was the major amino acid substitution in Yunnan, Zhejiang and Hubei isolates . The amino acid residue Q131K was not detected in Zhejiang and Hubei isolates , but it was very common in South Asian isolates such as those from Bangladesh , Thailand , India , and also in Yunnan isolates.
Among a total of 24 amino acid substitutions in Pvs28 from worldwide isolates, M52L and T140S were the most common substitutions from different geographic regions [15, 24–26, 31]. The amino acid residues L98I and L116V were very common in South Asia isolates such as those from Bangladesh , Thailand , India , and also in Yunnan isolates. The amino acid residue T65K was not identified in isolates from Yunnan and South Korea , but was common in those from Bangladesh , Thailand , India [25, 26] and Iran .
Genetic diversity of pvs25 and pvs28
There were a total of four polymorphic sites that generated five haplotypes in pvs25. There were a total of 12 polymorphic sites that yielded eight haplotypes in pvs28. Nucleotide diversity of pvs28 (π = 0.0034 ± 0.0012) was significantly higher than that of pvs25 (π = 0.0013 ± 0.0009) (p < 0.05) (Table 5).
In this study, we assessed the level of genetic diversity of pvs25 and pvs28 genes in P. vivax isolates from Yunnan Province, China. A total of four and eight amino acid substitutions have been identified in Pvs25 and Pvs28, respectively. Compared with P. vivax isolates from previous studies conducted in Asia, our results showed that the majority of the variant amino acids of Pvs25 and Pvs28 detected in Yunnan isolates were shared with those reported in isolates from South Asia (Bangladesh, Thailand and India) [15, 25]. Because of its location in southwest China, Yunnan isolates had mutant residues characteristic of South Asian isolates, such as mutant residues E97Q and Q131K in Pvs25, as well as L98I and L116V in Pvs28. It suggested that some parasite haplotypes may have different geographic distributions.
One of the major obstacles to the development of an effective malaria vaccine is the genetic polymorphism of many of the genes in natural parasite populations that otherwise would be promising vaccine candidates. A number of studies have shown that TBVs candidates have limited polymorphism compared to antigens expressed in asexual stage parasites [34, 35]. Our data showed that pvs25 (π = 0.0013) had a much lower level of diversity than P. vivax blood stage proteins DBP (π = 0.0086-0.0184)  and MSP1 (π = 0.1193-0.2055) . In addition, compared with P. falciparum sexual stage proteins pfs25 (π = 0.0035) , the lower genetic diversity of pvs25 suggested that genetic polymorphism of pvs25 was limited in Yunnan isolates. Although the nucleotide diversity of pvs28 was significantly higher (π = 0.0034) than that of pvs25, it was still conserved when compared with the blood stage vaccine candidates. Although direct sequencing of the PCR products in our study might underestimate the genetic diversity of these two genes, the effect should be minor since our study areas are malaria hypoendemic and mixed infections should not be prevalent. The limited diversity of sexual stage antigens such as P25 and P28 is perhaps attributed to the expression of these proteins in mosquito stages, which should avoid immune selection in the humans.
It is very important to evaluate the effect of genetic polymorphism of Pvs25 and Pvs28 on the efficacy of Sal-I based TBVs in different malaria epidemiological areas. The transmission-blocking assays conducted in Thailand clearly demonstrated that antisera to recombinant Pvs25 and Pvs28 based on the Sal-I strain of P. vivax recognized corresponding molecules expressed by field-isolated parasites in Thailand. Our data showed limited genetic diversity of pvs25 and pvs28 in the Yunnan isolates. The major amino acid haplotype of Pvs25 from this region was also shared with that from Thailand isolates, it provids a good prospect on implication of TBVs in China malaria areas.
Altogether, our study offers a first glimpse of the genetic diversity of two TBVs candidate antigens in Yunnan Province, China. Yet, this study only offered limited information of P. vivax parasites from Yunnan. Further studies on a larger sample scale and comparison analysis with Central China isolates will certainly improve our understanding of genetic polymorphisms of TBVs candidates in the malaria endemic areas of China.
TBVs are one of the important strategies for controlling malaria transmission. A key point in vaccine optimization is to understand the extent of genetic diversity of candidate antigens. The TBVs candidates Pvs25 and Pvs28 from Yunnan isolates showed limited genetic diversity. Further studies encompassing larger malaria endemic areas are needed to provide a thorough evaluation of situation in China.
World Health Organization: World malaria report 2009. Geneva. 2009
Mendis K, Sina BJ, Marchesini P, Carter R: The neglected burden of Plasmodium vivax malaria. Am J Trop Med Hyg. 2001, 64: 97-106.
Zhou SS, Wang Y, Xia ZG: Malaria situation in the People's Republic Of China in 2009. Zhongguo Ji Sheng Chong Xue Yu Ji Sheng Chong Bing Za Zhi. 2011, 29: 1-3.
Ministry of Health of the People's Republic of China: China to eliminate Malaria Action Plan (2010-2020). 2010
Cui L, Yan G, Sattabongkot J, Cao Y, Chen B, Chen X, Fan Q, Fang Q, Jongwutiwes S, Parker D, Sirichaisinthop J, Kyaw MP, Su XZ, Yang H, Yang Z, Wang B, Xu J, Zheng B, Zhong D, Zhou G: Malaria in the Greater Mekong Subregion: Heterogeneity and complexity. Acta Trop. 2011
Clements AC, Barnett AG, Cheng ZW, Snow RW, Zhou HN: Space-time variation of malaria incidence in Yunnan province, China. Malar J. 2009, 8: 180-10.1186/1475-2875-8-180.
Zhou SS, Wang Y, Tang LH: [Malaria situation in the People's Republic of China in 2006]. Zhongguo Ji Sheng Chong Xue Yu Ji Sheng Chong Bing Za Zhi. 2007, 25: 439-441.
Zhou SS, Wang Y, Fang W, Tang LH: [Malaria situation in the People's Republic Of China in 2007]. Zhongguo Ji Sheng Chong Xue Yu Ji Sheng Chong Bing Za Zhi. 2008, 26: 401-403.
Zhou SS, Wang Y, Fang W, Tang LH: [Malaria situation in the People's Republic of China in 2008]. Zhongguo Ji Sheng Chong Xue Yu Ji Sheng Chong Bing Za Zhi. 2009, 27: 457-455-456
Lin H, Lu L, Tian L, Zhou S, Wu H, Bi Y, Ho SC, Liu Q: Spatial and temporal distribution of falciparum malaria in China. Malar J. 2009, 8: 130-10.1186/1475-2875-8-130.
Kaslow DC: Transmission blocking vaccines. Malaria Vaccine Development. Edited by: Hofman SL. 1996, ASM books, 181-228.
Tsuboi T, Tachibana M, Kaneko O, Torii M: Transmission-blocking vaccine of vivax malaria. Parasitol Int. 2003, 52: 1-11. 10.1016/S1383-5769(02)00037-5.
Tomas AM, Margos G, Dimopoulos G, van Lin LH, de Koning-Ward TF, Sinha R, Lupetti P, Beetsma AL, Rodriguez MC, Karras M, Hager A, Mendoza J, Butcher GA, Kafatos F, Janse CJ, Waters AP, Sinden RE: P25 and P28 proteins of the malaria ookinete surface have multiple and partially redundant functions. EMBO J. 2001, 20: 3975-3983. 10.1093/emboj/20.15.3975.
Moon SU, Kim HH, Kim TS, Choi KM, Oh CM, Ahn YJ, Hwang SK, Sohn Y, Shin EH, Kim H, Lee HW: Blocking effect of a monoclonal antibody against recombinant Pvs25 on sporozoite development in Anopheles sinensis. Clin Vaccine Immunol. 2010, 17: 1183-1187. 10.1128/CVI.00101-10.
Sattabongkot J, Tsuboi T, Hisaeda H, Tachibana M, Suwanabun N, Rungruang T, Cao YM, Stowers AW, Sirichaisinthop J, Coleman RE, Torii M: Blocking of transmission to mosquitoes by antibody to Plasmodium vivax malaria vaccine candidates Pvs25 and Pvs28 despite antigenic polymorphism in field isolates. Am J Trop Med Hyg. 2003, 69: 536-541.
Miyata T, Harakuni T, Tsuboi T, Sattabongkot J, Kohama H, Tachibana M, Matsuzaki G, Torii M, Arakawa T: Plasmodium vivax ookinete surface protein Pvs25 linked to cholera toxin B subunit induces potent transmission-blocking immunity by intranasal as well as subcutaneous immunization. Infect Immun. 2010, 78: 3773-3782. 10.1128/IAI.00306-10.
Blagborough AM, Yoshida S, Sattabongkot J, Tsuboi T, Sinden RE: Intranasal and intramuscular immunization with Baculovirus Dual Expression System-based Pvs25 vaccine substantially blocks Plasmodium vivax transmission. Vaccine. 2010, 28: 6014-6020. 10.1016/j.vaccine.2010.06.100.
Saul A, Hensmann M, Sattabongkot J, Collins WE, Barnwell JW, Langermans JA, Wu Y, Long CA, Dubovsky F, Thomas AW: Immunogenicity in rhesus of the Plasmodium vivax mosquito stage antigen Pvs25H with Alhydrogel and Montanide ISA 720. Parasite Immunol. 2007, 29: 525-533. 10.1111/j.1365-3024.2007.00971.x.
Malkin EM, Durbin AP, Diemert DJ, Sattabongkot J, Wu Y, Miura K, Long CA, Lambert L, Miles AP, Wang J, Stowers A, Miller LH, Saul A: Phase 1 vaccine trial of Pvs25H: a transmission blocking vaccine for Plasmodium vivax malaria. Vaccine. 2005, 23: 3131-3138. 10.1016/j.vaccine.2004.12.019.
Wu Y, Ellis RD, Shaffer D, Fontes E, Malkin EM, Mahanty S, Fay MP, Narum D, Rausch K, Miles AP, Aebig J, Orcutt A, Muratova O, Song G, Lambert L, Zhu D, Miura K, Long C, Saul A, Miller LH, Durbin AP: Phase 1 trial of malaria transmission blocking vaccine candidates Pfs25 and Pvs25 formulated with montanide ISA 51. PLoS One. 2008, 3: e2636-10.1371/journal.pone.0002636.
Tsuboi T, Kaslow DC, Gozar MM, Tachibana M, Cao YM, Torii M: Sequence polymorphism in two novel Plasmodium vivax ookinete surface proteins, Pvs25 and Pvs28, that are malaria transmission-blocking vaccine candidates. Mol Med. 1998, 4: 772-782.
Han ET, Lee WJ, Sattabongkot J, Jang JW, Nam MH, An SS, Suh I, Lim CS: Sequence polymorphisms of Plasmodium vivax ookinete surface proteins (Pvs25 and Pvs28) from clinical isolates in Korea. Trop Med Int Health. 2010, 15: 1072-1076.
Escalante AA, Cornejo OE, Freeland DE, Poe AC, Durrego E, Collins WE, Lal AA: A monkey's tale: the origin of Plasmodium vivax as a human malaria parasite. Proc Natl Acad Sci USA. 2005, 102: 1980-1985. 10.1073/pnas.0409652102.
Zakeri S, Razavi S, Djadid ND: Genetic diversity of transmission blocking vaccine candidate (Pvs25 and Pvs28) antigen in Plasmodium vivax clinical isolates from Iran. Acta Trop. 2009, 109: 176-180. 10.1016/j.actatropica.2008.09.012.
Tsuboi T, Kaneko O, Cao YM, Tachibana M, Yoshihiro Y, Nagao T, Kanbara H, Torii M: A rapid genotyping method for the vivax malaria transmission-blocking vaccine candidates, Pvs25 and Pvs28. Parasitol Int. 2004, 53: 211-216. 10.1016/j.parint.2004.01.012.
Prajapati SK, Joshi H, Dua VK: Antigenic repertoire of Plasmodium vivax transmission-blocking vaccine candidates from the Indian subcontinent. Malar J. 2011, 10: 111-10.1186/1475-2875-10-111.
Hui FM, Xu B, Chen ZW, Cheng X, Liang L, Huang HB, Fang LQ, Yang H, Zhou HN, Yang HL, Zhou XN, Cao WC, Gong P: Spatio-temporal distribution of malaria in Yunnan Province, China. Am J Trop Med Hyg. 2009, 81: 503-509.
Cox FE: History of the discovery of the malaria parasites and their vectors. Parasit Vectors. 2010, 3: 5-10.1186/1756-3305-3-5.
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.
Gonzalez-Ceron L, Alvarado-Delgado A, Martinez-Barnetche J, Rodriguez MH, Ovilla-Munoz M, Perez F, Hernandez-Avila JE, Sandoval MA, Rodriguez Mdel C, Villarreal-Trevino C: Sequence variation of ookinete surface proteins Pvs25 and Pvs28 of Plasmodium vivax isolates from Southern Mexico and their association to local anophelines infectivity. Infect Genet Evol. 2010, 10: 645-654. 10.1016/j.meegid.2010.03.014.
Zheng L, Xu WM, Liu YJ, Yang YM, Cao YM: Transmission-blocking vaccine candidate of Plasmodium vivax Pvs25 is highly conservative among Chinese isolates. Zhongguo Ji Sheng Chong Xue Yu Ji Sheng Chong Bing Za Zhi. 2004, 22: 16-19.
Sattabongkot J, Tsuboi T, Zollner GE, Sirichaisinthop J, Cui L: Plasmodium vivax transmission: chances for control?. Trends Parasitol. 2004, 20: 192-198. 10.1016/j.pt.2004.02.001.
Escalante AA, Lal AA, Ayala FJ: Genetic polymorphism and natural selection in the malaria parasite Plasmodium falciparum. Genetics. 1998, 149: 189-202.
Barry AE, Schultz L, Buckee CO, Reeder JC: Contrasting population structures of the genes encoding ten leading vaccine-candidate antigens of the human malaria parasite, Plasmodium falciparum. PLoS One. 2009, 4: e8497-10.1371/journal.pone.0008497.
Cole-Tobian J, King CL: Diversity and natural selection in Plasmodium vivax Duffy binding protein gene. Mol Biochem Parasitol. 2003, 127: 121-132. 10.1016/S0166-6851(02)00327-4.
Putaporntip C, Jongwutiwes S, Sakihama N, Ferreira MU, Kho WG, 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.
Acknowledgements and Funding
We are grateful to Drs. Osamu Kaneko, Takafumi Tsuboi and Motomi Torii for the helpful suggestions. The work was supported by the National S & T Major Program (2008ZX10004-011), National Natural Science Foundation of China (30972774) and National Institutes of Health (U19AI089672).
The authors declare that they have no competing interests.
YC conceived of the study and helped to draft the manuscript. YY collected the blood spots specimens. HF carried out the studies, statistical analysis and drafted the manuscript. LZ, XZ, GW, YP, YL, YL participated in the molecular genetic studies and sequence alignment. LC helped with statistical analysis and critically revised the manuscript. All authors contributed to the writing of the manuscript and approved the submitted version of the manuscript.