SciELO - Scientific Electronic Library Online

 
vol.108 issue3New structures for goat corrals to control peridomestic populations of Triatoma infestans (Hemiptera: Reduviidae) in the Gran Chaco of ArgentinaSerodiagnosis of human neurocysticercosis using antigenic components of Taenia solium metacestodes derived from the unbound fraction from jacalin affinity chromatography author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

Share


Memórias do Instituto Oswaldo Cruz

Print version ISSN 0074-0276

Mem. Inst. Oswaldo Cruz vol.108 no.3 Rio de Janeiro May 2013

http://dx.doi.org/10.1590/S0074-02762013000300015 

Articles

Sequence polymorphisms in Pvs48/45 and Pvs47 gametocyte and gamete surface proteins in Plasmodium vivax isolated in Korea

Mi Kyung Woo1 

Kyeong Ah Kim1 

JuYeon Kim1 

Jun Seo Oh2 

Eun Taek Han3 

Seong Soo A An4 

Chae Seung Lim1  + 

1Korea University Guro Hospital, Department of Laboratory Medicine

2Korea University, College of Medicine, Graduate School of Medicine, Laboratory of Cellular Oncology, Seoul, Republic of Korea

3Kangwon National University, College of Medicine, Department of Parasitology, Chuncheon, Republic of Korea

4Gachon University, Gachon BioNano Research Institute, Department of BioNano Technology, Gyeonggi Do, Republic of Korea

ABSTRACT

Nucleotide sequence analyses of the Pvs48/45 and Pvs47 genes were conducted in 46 malaria patients from the Republic of Korea (ROK) (n = 40) and returning travellers from India (n = 3) and Indonesia (n = 3). The domain structures, which were based on cysteine residue position and secondary protein structure, were similar between Plasmodium vivax (Pvs48/45 and Pvs47) and Plasmodium falciparum (Pfs48/45 and Pfs47). In comparison to the Sal-1 reference strain (Pvs48/45, PVX_083235 and Pvs47, PVX_083240), Korean isolates revealed seven polymorphisms (E35K, H211N, K250N, D335Y, A376T, I380T and K418R) in Pvs48/45. These isolates could be divided into five haplotypes with the two major types having frequencies of 47.5% and 20%, respectivelfy. In Pvs47, 10 polymorphisms (F22L, F24L, K27E, D31N, V230I, M233I, E240D, I262T, I273M and A373V) were found and they could be divided into four haplotypes with one major type having a frequency of 75%. The Pvs48/45 isolates from India showed a unique amino acid substitution site (K26R). Compared to the Sal-1 and ROK isolates, the Pvs47 isolates from travellers returning from India and Indonesia had amino acid substitutions (S57T and I262K). The current data may contribute to the development of the malaria transmission-blocking vaccine in future clinical trials.

Key words: Plasmodium vivax; Pvs48/45; Pvs47; polymorphism; malaria vaccine; transmission blocking vaccine; Republic of Korea

Malaria is a highly infectious disease and has the highest worldwide mortality rate. According to the 2011 World Malaria Report, the number of suspected malaria cases reached 216 million with 655,000 malaria-related deaths in 2010. The prevalence of Plasmodium vivax infection in 2009 was 71% in South America and 49% in both Southern Asia and the Western Pacific. In the Republic of Korea (ROK), P. vivax is the only indigenous strain of malaria (Mendis et al. 2001, WHO 2011).

The life cycle of the malaria parasite is divided into asexual and sexual stages. The malaria parasite can be spread from mosquitoes to humans and can reside and multiply in the liver and red blood cells (RBCs). The asexually transmitted parasites in RBCs can then induce malaria symptoms.

To reduce communal malaria infections, a malaria transmission-blocking vaccine (TBV) has been investigated using surface proteins of male and female gametocytes, zygotes and ookinetes from Plasmodium falcipa-rum (Roeffen et al. 1995, 2001, van Schaijk et al. 2006, Outchkourov et al. 2007, 2008). Pfs48/45 is expressed on the surface of male and female gametocytes/gametes and is a necessary factor for male fertility (van Dijk et al. 2001). The Pfs48/45 protein can induce antibodies against the surface proteins of gametocytes and the monoclonal antibodies against Pfs48/45 can block the transmission of the parasite to mosquitoes (Carter et al. 2000, Carter 2001, Stowers & Carter 2001). With additional supporting results from other investigators, Pfs48/45 has been accepted as an effective TBV candidate against P. falciparum (Roeffen et al. 1995, van Dijk et al. 2001, Drakeley et al. 2006, Pradel 2007). In addition, Pfs47, which is expressed in female gametocytes and gametes, has also demonstrated an ability to induce antibody responses and to serve as a TBV candidate against P. falciparum (van Schaijk et al. 2006, Pradel 2007).

Because P. vivax is the second most common infectious Plasmodium species, developing a TBV using the surface proteins of gametocytes has been proposed as a possible public health defence mechanism. Pvs48/45 of P. vivax, which is a homolog of Pfs48/45 of P. falciparum, is also a TBV candidate (Galinski & Barnwell 2008). However, TBV research on Pvs48/45 is in the early stages due to difficulties in producing recombinant Pvs48/45 protein. Pvs47 of P. vivax has not yet been investigated as a TBV candidate (Galinski & Barnwell 2008). Hence, additional studies using the Pvs48/45 and Pvs47 proteins of P. vivax should be performed in the development of a P. vivax TBV.

P. vivax malaria is an endemic disease in the ROK. The incidence of P. vivax decreased sharply with the launch of the WHO malaria control project resulting in an almost complete eradication of P. vivax. Unexpectedly, P. vivax re-emerged in 1993 and malaria cases climbed to 4,141 by 2000 (Yeom et al. 2005). Fortunately, the number of patients with malaria decreased to 841 in 2011. Malaria transmission in the ROK occurs mostly in the regions near the Delimited Militarised Zone (DMZ). Hence, this study was conducted in Kimpo and Pajoo, which are high-risk malaria areas, located in the northwest region of the ROK near the DMZ. Korean P. vivax infection presents with mild, but characteristic malaria symptoms. From June-October, the majority of malaria transmission occurs in areas near the DMZ (Yeom et al. 2005). Vector control and chemoprophylaxis have been performed in the ROK near the DMZ, but the malaria outbreaks are not controlled in North Korea. To control the P. vivax malaria infection rate on the Korean peninsula, use of the malaria TBV against mosquito hosts in risky areas would be an ideal and long-term way to eliminate malaria in both regions. Therefore, we analysed the polymorphisms in the Pvs48/45 and Pvs47 proteins from the clinical P. vivax isolates from the ROK. The results from these Korean isolates may provide valuable baseline data for the development of a TBV and future clinical trials.

SUBJECTS, MATERIALS AND METHODS

This study was conducted between 2006-2010 at the Korea University Hospital at Guro, ROK. Whole blood samples were collected into ethylenediaminetetraacetic acid tubes from febrile patients for malaria diagnosis. The confirmed samples, which were infected with indigenous Korean P. vivax according to Giemsa-stain microscopic analysis, were included in this study. The patients had no history of travel to a foreign country within three years prior to the study. The patients were 22-49 years old (mean = 37.6 years) with a male/female ratio of 5.5:1. The following samples were collected: seven in 2006, 16 in 2007, eight in 2008 and nine in 2010. Each patient with malaria received treatments of oral chloroquine (800 mg of hydroxychloroquine sulphate) followed by 400 mg after 6 h, 24 h and 48 h for a total of 2,000 mg. After completing chloroquine therapy, oral primaquine was prescribed at 15 mg/day for 14 days. All of the patients responded well to the standard chloroquine therapy. The parasitaemia level was indirectly calculated by assessing the number of parasites per 200 white blood cells (WBCs) in the blood film and using the results of the WBC counts from an automated haematology analyser (Cell-Dyn 4000; Abbott Diagnostics, Santa Clara, CA, USA). The mean ± standard deviation parasitaemia level was 4,236 (± 13,210)/µL with a range of 46-31,000/µL. These patients were residents of or had a history of visiting Kimpo (n = 26) and Pajoo (n = 14). Informed consent was obtained before collection of the blood samples. Six isolates from travellers returning from Indonesia and India were compared against the Korean sequences. The blood samples were stored at -80ºC. The current protocol was approved by the institutional review board of the Korea University Guro Hospital.

The P. vivax DNA was extracted from whole blood using an AquaPure Genomic DNA Kit (Bio-Rad Laboratories, USA) according to the manufacturer's instructions. The Pvs48/45 and Pvs47 P. vivax genes were amplified by polymerase chain reaction (PCR) using the indicated primers. For amplification of the Pvs48/45 gene, the forward and reverse primers were 5′-CACGCAGATTGATGTCGC-3′ (CM000453, 289234-289251) and 5′-CAAACAGGTGCACGCGC-3′ (CM000453, 290823-290839), respectively. For amplification of the Pvs47 gene, the forward and reverse primers were 5′-CACACCACCGCAAACAGG-3′ (CM000453, 286169-286186) and 5′-GTGCACATTCCGCGGTTG-3′ (CM000453, 287676-287693), respectively. Primers were used at a final concentration of 0.1 µM in 25 µL of reaction mixture (10 mM Tris-HCl, pH 8.3; 50 mM KCl; 1.5 mM MgCl2; 0.2 mM deoxyribonucleotide triphosphate) with 1 µL of extracted DNA from whole blood and 2.5 units of ExTaq polymerase (TaKaRa Bio Inc, Japan). The PCR mixtures were cycled 35 times using a TaKaRa PCR Thermal Cycler (TaKaRa Bio Inc, Japan) with the following conditions: each cycle included 1 min of denaturation at 95ºC, 1 min of annealing at 67ºC for Pvs48/45 or 65°C for Pvs47 and 1.5 min of extension at 72ºC. Amplified DNA fragments of Pvs48/45 and Pvs47 were size fractionated by electrophoresis using a 1.5% agarose gel containing 0.5 mg/mL ethidium bromide. Amplified PCR products were gel purified using the QIAGEN DNA purification kit (QIAGEN, Germany) following the manufacturer's instructions. The sequences were determined using the BigDye Terminator v3.1 Cycle Sequencing Kit and ABI 3730XL automatic DNA sequencer. The sequencing primers for the Pvs48/45 gene were 5′-CGGGTGCTTTAAAAATGGAA-3′ (pvs48R2; CM000453, 290138-290157) and 5′-TGTAAAATCTG-CGGACGTGA-3′ (pvs48F2; CM000453, 290162-290181). The sequencing primers for the Pvs47 gene were 5′-GAAGAAAGGGGAGGACCAAG-3′ (pvs47F2; PVX_083240, 950-969) and 5′-TTTTCCGT-TTTCTGCGTCTT-3′ (pvs47R2; PVX_083240, 1336-1317). The 46 sequences for the Pvs48/45 and Pvs47 genes were aligned with the reference sequence in addition to the Pvs48 and Pvs47 genes from the Sal-1 strain (PVX_083235 and PVX_083240, respectively) using MultAlin software (Corpet 1988), which is a progressive multiple-sequence alignment algorithm. Corrections were carried out manually. The newly obtained sequences were submitted to the GenBank database (accessions JN582344-JN582353 and JQ435579-JQ435660). A phylogenetic tree was constructed based on the maximum-likelihood method from the MEGA 5.0 package (Tamura et al. 2011) and the reliability of each branch was assessed using 500 bootstrap replications. The rate of nucleotide substitutions in pairwise comparisons between haplotypes was computed by the number of synonymous substitutions per synonymous sites and nonsynonymous substitutions per nonsynonymous sites (Nei & Gojobori 1986) using the MEGA 5.0 package (Tamura et al. 2011).

To predict each domain in Pvs48/45 and Pvs47, the Pvs48/45 (PVX_083235) and Pvs47 (PVX_083240) gene sequences were aligned with the Pfs48/45 (XM_001350145) and Pfs47 (XM_001350146) gene sequences of P. falciparum, respectively, using the MultAlin computer program and analysing secondary structure using Professor_King's program (aber.ac.uk/~phiwww/prof/) (Ouali & King 2000). Three-dimensional (3D) structure was predicted using the template-based protein structure internet modelling program, RaptorX (Källberg et al. 2012).

RESULTS

Prediction of cysteine-rich domains (CMs) in Pvs48/45 and Pvs47 - A 56% amino acid sequence homology existed between Pvs48/45 (PVX_083235) and Pfs48/45 (PF13_0247) and 42% homology existed between Pvs47 (PVX_083240) and Pfs47 (PF13_0248). The conserved cysteine residues were also found in Pvs48/45 and Pvs47 (Fig. 1). The predicted domains for comparing the positions of the cysteine residues were found at amino acid positions 45-182 (CM1), 183-292 (CM2) and 296-435 (CM3) in Pvs48/45 and 26-175 (CM1), 176-277 (CM2) and 279-417 (CM3) in Pvs47 (van Dijk et al. 2001, Gerloff et al. 2005) (Fig. 1). Pvs48/45 and Pvs47 contained conserved cysteines that formed the disulfide bond in a position similar to that of Pfs48/45 and Pfs47, respectively. Furthermore, in the secondary structures of Pvs48/45, Pvs47, Pfs48/45 and Pfs47, the locations of β-sheets were similar between P. falciparum and P. vivax (Fig. 1). Therefore, the predicted domain and epitope in Pvs48/45 and Pvs47 were similar to those reported for Pfs48/45 and Pfs47 (Gerloff et al. 2005).

Fig. 1 locations of the cysteine-rich domain (CM)1, CM2 and CM3 domains of Pvs48/45 (A) and Pvs47 (B) in the Plasmodium vivax Sal-1 strain. The locations of the above-mentioned domains were predicted via comparison against the Pfs48/45 and Pfs47 sequences of Plasmodium falciparum, respectively, and comparison of the positions of the cysteine residues and β-sheets (van Dijk et al. 2001, Gerloff et al. 2005). A: three domains in Pfs48/45 were recognised using different transmission blocking monoclonal antibodies and their targeted epitopes (I, IIb, III and V) were described (Roeffen et al. 2001); B: the Pvs47 epitope was not yet determined; ►◄: range of CM; dark gray: cysteine; light gray: β-sheet. 

When polymorphisms from Pvs48/45 were used to obtain secondary structure information using RaptorX, a template-based protein structure modelling program (Källberg et al. 2012) positions H211N, K250N and K418R was found to belong to a loop/irregular structure without disorder, which may provide a new epitope for antibody generation (Figs 2, 3A-C). In contrast, positions E35K, D335Y, A376T and I380T were found to belong to a beta-sheet region, which caused these residues to be buried in the inner core of the protein, thereby resulting in no accessibility for antibody binding.

Fig. 2 secondary structure prediction of Pvs48/45 protein. Red colour represented the loop/irregular structures and blue and yellow referred to alpha helical and beta-sheet conformations, respectively. Asterisks indicated the mutations. 

Fig. 3 three-dimensional structure prediction of Pvs48/45 (A-C) and Pvs47 (D-F) proteins. Red coloured atoms represented the mutations in Pvs48/45 (A-C) and blue referred to the mutations in Pvs47 (D-F), respectively. 

For Pvs47, polymorphisms of F22L, F24L, D31N and A373V in Figs 3D-F, 4 may fall into the loop regions, whereas mutations V230I, M233I and E240D may participate in forming alpha helix and mutations of I262T, I273M and K27E mutations may result in beta-sheet conformations.

Fig. 4 secondary structure prediction of Pvs47 protein. Red colour represented the loop/irregular structures and blue and yellow referred to alpha helical and beta-sheet conformations, respectively. Asterisks indicated the mutations. 

Polymorphisms in Pvs48/45 and Pvs47 - From the Pvs48/45 nucleotide and amino acid sequence analyses, one synonymous and eight nonsynonymous substitutions were found in comparison to those in the Sal-1 strain (Table I) in the present study.

TABLE I Polymorphisms of Pvs48/45 in Plasmodium vivax isolates from Republic of Korea (ROK), India and Indonesia in comparison with Sal-1 strain 

  Pre-CM a
CM2
CM3
 
Position of nucleotides substitution (amino acids) 46 77 (26) 103 (35) 631 (211) 750 (250) 1003 (335) 1126 (376) 1139 (380) 1253 (418) ROK isolates n (%)
Sal-1 b C AAA K GAG E CAT H AAA K GAT D GCT A ATT I AAG K -
KPVS48-2010-05 (JN582348) C AAA K AAG K AAT N AAC N TAT Y ACT T ATT I AGG R 19 (47.5)
KPVS48-2010-34 (JN582349) C AAA K AAG K AAT N AAC N TAT Y ACT T ACT T AGG R 8 (20)
KPVS48-2010-06 (JN582350) C AAA K GAG E AAT N AAC N GAT D GCT A ACT T AAG K 6 (15)
KPVS48-2008-21 (JN582352) C AAA K GAG E AAT N AAC N TAT Y ACT T ATT I AGG R 5 (12.5)
KPVS48-2007-46 (JN582353) C AAA K GAG E CAT H AAC N TAT Y ACT T ATT I AGG R 2 (5)
India-1 (JQ435655) C AAA K GAG E AAT N AAC N TAT Y ACT T ATT I AGG R -
India-2 (JQ435656) C AGA R GAG E AAT N AAC N TAT Y ACT T ATT I AGG R -
India-3 (JQ435657) C AGA R GAG E AAT N AAC N TAT Y ACT T ATT I AGG R -
Indonesia-1 (JQ435658) T AAA K AAG K AAT N AAC N TAT Y ACT T ATT I AGG R -
Indonesia-2 (JQ435659) C AAA K GAG E AAT N AAC N TAT Y ACT T ATT I AGG R -
Indonesia-3 (JQ435660) C AAA K GAG E AAT N AAC N TAT Y ACT T ATT I AGG R -

a: this region was located between signal peptide and cysteine-rich domain (CM)1; b: reference sequence of Pvs48/45 of P. vivax (Sal-1 strain, PVX_083235). Nucleotides and amino acids polymorphisms were shaded at corresponding positions.

In the Korean isolates, seven nonsynonymous substitutions (E35K, H211N, K250N, D335Y, A376T, I380T and K418R) and five haplotypes were found (Table I). The two major types, KPVS48-2010-5 (KNNYTIR) and KPVS48-2010-34 (KNNYTTR), had frequencies of 47.5% and 20%, respectively (Table I). The two minor types, KPVS48-2010-6 and KPVS48-2008-21, had frequencies of 15% and 12.5%, respectively (Table I). Another haplotype, KPVS48-2007-46, was found in two Korean isolates with a frequency of 5%. In comparison to SAL-1 (Table I), six amino acid substitutions at the K26R, H211N, K250N, D335Y, A376T and K418R positions (RNNYTR), were found in the isolates of patients returning from India. Indonesian isolates showed one synonymous substitution at 46T>C and six nonsynonymous substitutions (KNNYTR) at the respective positions. The K26R substitution was found only in the isolates of the patients returning from India and the I380T substitution was found only from Korean isolates.

For the nucleotide and amino acid sequence analysis of Pvs47, two synonymous (78A>G and 1116C>T) and 11 nonsynonymous substitutions (F22L, F24L, K27E, D31N, S57T, V230I, M233I, E240D, I262T/K, I273M, A373V) were found in comparison to the Sal-1 strain.

In the Korean P. vivax isolates, one synonymous (78A>G) and 10 nonsynonymous substitutions (F22L, F24L, K27E, D31N, V230I, M233I, E240D, I262T, I273M, A373V) were found in comparison to the Sal-1 strain. The Korean isolates revealed four haplotypes. The major haplotype (75% frequency) was KPVS47-2010-15 with LLEIITV substitutions (Table II). The second major haplotype (20% frequency) was KPVS47-2010-3 with LLEIIT substitutions. The two minor haplotypes, KPVS47-2007-24 with LLENIITV substitutions and KPVS47-2008-20 with LLEDTMV substitutions, were found only in a single Korean isolate. In the isolates from the patients returning from India, one synonymous (78A>G) and six nonsynonymous substitutions were found at the following positions: F22L, F24L, K27E, S57T, M233I, I262T. Among the amino acid substitution sites, the S57T of CM1 was found only in the isolates of the returning patients from India. From the isolates of the patients returning from Indonesia, one synonymous (78A>G) and seven nonsynonymous substitutions were found at F22L, F24L, K27E, V230I, M233I, I262K and A373V positions in comparison with SAL-1 (Table II). Compared with ROK and Sal-1 isolates, the Indian and Indonesian isolates had specific amino acid substitutions at S57T and I262K in CM1 and CM2, respectively.

TABLE II Polymorphisms of Pvs47 in Plasmodium vivax isolates from Republic of Korea (ROK), India and Indonesia in comparison with Sal-1 strain 

  Pre-CM a
  CM1
 
Position of nucleotides substitution (amino acids) 66 (22) 70 (24)   78 79 (27) 91 (31) 170 (57)  


Sal-1 b TTC F TTT F   G AAG K GAC D AGC S  
KPVS47-2010-15 (JN582344) TTA L CTT L   A GAG E GAC D AGC S  
KPVS47-2010-03 (JN582345) TTA L CTT L   A GAG E GAC D AGC S  
KPVS47-2007-24 (JN582346) TTA L CTT L   A GAG E AAC N AGC S  
KPVS47-2008-20 (JN582346) TTA L CTT L   A GAG E GAC D AGC S  
India-1 (JQ435615) TTA L CTT L   A GAG E GAC D ACC T  
India-2 (JQ435616) TTA L CTT L   A GAG E GAC D ACC T  
India-3 (JQ435617) TTA L CTT L   A GAG E GAC D ACC T  
Indonesia-1 (JQ435618) TTA L CTT L   A GAG E GAC D AGC S  
Indonesia-2 (JQ435619) TTA L CTT L   A GAG E GAC D AGC S  
Indonesia-3 (JQ435620) TTA L CTT L   A GAG E GAC D AGC S  

  CM2
CM3
688 (230) 699 (233) 720 (240) 785 (262) 819 (273) 1116   1118 (373) ROK isolates n (%)
Sal-1 b GTT V ATG M GAG E ATA I ATA I T   GCA A -
KPVS47-2010-15 (JN582344) ATT I ATT I GAG E ACA T ATA I T   GTA V 30 (75)
KPVS47-2010-03 (JN582345) ATT I ATT I GAG E ACA T ATA I T   GCA A 8 (20)
KPVS47-2007-24 (JN582346) ATT I ATT I GAG E ACA T ATA I T   GTA V 1 (2.5)
KPVS47-2008-20 (JN582346) GTT V ATG M GAT D ACA T ATG M T   GTA V 1 (2.5)
India-1 (JQ435615) GTT V ATT I GAG E AAA K ATA I T   GCA A -
India-2 (JQ435616) GTT V ATT I GAG E AAA K ATA I T   GCA A -
India-3 (JQ435617) GTT V ATT I GAG E AAA K ATA I T   GCA A -
Indonesia-1 (JQ435618) ATT I ATT I GAG E AAA K ATA I C   GTA V -
Indonesia-2 (JQ435619) ATT I ATT I GAG E AAA K ATA I T   GCA A -
Indonesia-3 (JQ435620) ATT I ATT I GAG E AAA K ATA I T   GCA A -

a: this region was located between signal peptide and cysteine-rich domain (CM)1; b: reference sequence of Pvs47 of P. vivax (Sal-1 strain, PVX_083240). Nucleotides and amino acids polymorphisms were shaded at corresponding positions.

The full-length sequences of the 47 isolates of Pvs48/45 and Pvs47 were determined to assess genetic diversity as compared to Sal-1, which is the worldwide P. vivax reference sequence. The nucleotide diversity per site was 0.00147 in Pvs48/45 and 0.00109 in Pvs47 with an average number of nucleotide differences of 1.987 and 1.424, respectively. The synonymous and nonsynonymous nucleotide diversity per site was determined to be 0.00186 and 0.00029 in Pvs48/45, respectively and 0.00014 and 0.00132 in Pvs47, respectively.

Phylogenetic analysis showed that the Pvs48/45 and Pvs47 haplotypes from the Asian isolates appeared to be scattered without clusters and not influenced by their geographic origins (Fig. 5).

Fig. 5: maximum-likelihood tree of Plasmodium vivax populations from Korea and world isolates inferred from Pvs48/45 (A) and Pvs47 (B) based on pairwise comparison. Values of the bootstrap confidence level are given as per cent branch searches out of 500 replications. The distance corresponding to five changes per 10,000 positions is indicated on the lower left. 

DISCUSSION

Pvs48/45 and Pvs47 proteins belong to cysteine-rich protein families that contain three domains, including A types (CM1 and CM3) and a B type (CM2), with conserved cysteine residue motifs (Gerloff et al. 2005, Pradel 2007).

Four epitopes, including epitope V (CM1), epitope IIb (CM2), epitope III (CM2) and epitope I (CM3), have been constructed in P. falciparum Pfs48/45 (Gerloff et al. 2005). Based on the results from the membrane-feeding assay, the transmission of parasites is blocked by monoclonal antibodies against epitope I and the effectiveness of blocking transmission by monoclonal antibodies against epitopes IIb, III and V is much weaker (Roeffen et al. 1995, 2001). Other studies have suggested that epitope I (CM3) and epitope III (CM2) may be crucial domains for inducing effective TB antibodies (Outchkourov et al. 2007, 2008). Hence, the baseline data for polymorphisms in CM2 and CM3 would be useful for developing a Pfs48/45-based TBV. The present results predicted the locations of CM1, CM2 and CM3 in P. vivax Pvs48/45 in comparison to the Sal-1 strain (PVX_083235) (Fig. 1).

Haplotypes with specific substitutions in the CM domains are shown. Nucleotide and amino acid substitutions in Pvs48/45 from Korean isolates were located at Pre-CM (E35K), CM2 (H211N and K250N) and CM3 (D335Y, A376T, I380T and K418R) as shown in Table I. The substitutions in the isolates from India and Indonesia were also similar to the ROK-type isolates, except for 46T>C, K26R and E35K in PreCM, as well as I380T in CM3. In the Pfs48/45 (P. falciparum) polymorphism study by Anthony et al. (2007) five nonsynonymous substitutions were found in CM2 (2 substitutions) and CM3 (3 substitutions). Our results and those from the study by Anthony et al. (2007) revealed no substitution in the CM1 of Pvs48/45 (Table I). As mentioned above, CM2 and CM3 could be highly efficient as a TBV target antigen. Hence, the substitutions in CM2 and CM3 of Pvs48/45 and Pfs48/45 could have occurred as immune protection against the host. Therefore, our polymorphism data and the above epitope studies suggest that CM2 and CM3 of Pvs48/45 could be crucial domain sets for developing an efficient TBV.

Pvs47 also had three domains (CM1, CM2 and CM3), but the results of the epitope study for both Pvs47 and Pfs47 are not yet available. As shown in Table II, two synonymous and 11 nonsynonymous polymorphisms were found. The following substitutions of Pvs47 were found in the Korean isolates: F22L and F24L in Pre-CM, K27E and D31N in CM1, V230I, M233I, E240D, I262T and I273M in CM2 and A373V in CM3. Among the isolates from India and Indonesia, three substitutions in CM1 (A>G, K27E and S57T) and four substitutions in CM2 (V230I, M233I, I262K and 1666C>T) were found. However, V230I was found only in the isolates from Indonesia. In the study by (Anthony et al. 2007), diverse substitutions were also found in all domains. Therefore, the substitutions in Pvs47 and Pfs47 occurred more often than in Pvs48/45. To test if these substitutions lead to immune protection, an immune response study with recombinant Pvs47 is necessary.

From the 3D protein structure prediction of Pvs48/45, the H211N, K250N and K418R polymorphisms in the loop regions might serve as potential sites for targeted vaccine generation. However, the potential influence of the buried residues on the global protein structure cannot be ruled out, especially if the protein becomes unfolded, thereby exposing the buried residues.

For Pvs47, mutations found in the loop regions, F22L, F24L, D31N, A373V, may directly influence epitope recognition and residues in the alpha helical or beta-sheet conformations may not be directly involved when creating a potent TBV.

The effectiveness of TBV vaccines is dependent on the concentration of the induced antibodies and the quantity of the fertilisation-competent gametocytes. The cysteine residues in the central and C-terminal modules of Pvs48/45 are crucial for proper presentation of the TB epitopes (Outchkourov et al. 2007). However, amino acid substitutions of cysteine residues in Pvs48/45 and Pvs47 were not found in the present study. Therefore, the present results suggested that a TBV based on Pvs48/45 and Pvs47 may induce a TB antibody against antigens with various polymorphisms. Importantly, P. vivax isolates in this study were collected between June 2006-October 2010. Yearly fluctuations in isolate types were not significant.

We expect that the polymorphism analyses of Pvs48/45 and Pvs47 of P. vivax from the ROK, India and Indonesia may have significant contributions towards the development of an efficient TBV against Pvs48/45 and Pvs47 proteins in future clinical trials.

REFERENCES

Anthony TG, Polley SD, Vogler AP, Conway DJ 2007. Evidence of non-neutral polymorphism in Plasmodium falciparum gamete surface protein genes Pfs47 and Pfs48/45. Mol Biochem Parasitol 156: 117-123. [ Links ]

Carter R 2001. Transmission blocking malaria vaccines. Vaccine 19: 2309-2314. [ Links ]

Carter R, Mendis KN, Miller LH, Molineaux L, Saul A 2000. Malaria transmission-blocking vaccines - how can their development be supported? Nat Med 6: 241-244. [ Links ]

Corpet F 1988. Multiple sequence alignment with hierarchical clustering. Nucleic Acids Res 16: 10881-10890. [ Links ]

Drakeley CJ, Bousema JT, Akim NI, Teelen K, Roeffen W, Lensen AH, Bolmer M, Eling W, Sauerwein RW 2006. Transmission-reducing immunity is inversely related to age in Plasmodium falciparum gametocyte carriers. Parasite Immunol 28: 185-190. [ Links ]

Galinski MR, Barnwell JW 2008. Plasmodium vivax: who cares? Malar J 7 (Suppl. 1): S9. [ Links ]

Gerloff DL, Creasey A, Maslau S, Carter R 2005. Structural models for the protein family characterized by gamete surface protein Pfs230 of Plasmodium falciparum. Proc Natl Acad Sci USA 102: 13598-13603. [ Links ]

Källberg M, Wang H, Wang S, Peng J, Wang Z, Lu H, Xu J 2012. Template-based protein structure modeling using the RaptorX web server. Nat Protoc 7: 1511-1522. [ Links ]

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

Nei M, Gojobori T 1986. Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol 5: 418-426. [ Links ]

Ouali M, King RD 2000. Cascaded multiple classifiers for secondary structure prediction. Protein Sci 9: 1162-1176. [ Links ]

Outchkourov N, Vermunt A, Jansen J, Kaan A, Roeffen W, Teelen K, Lasonder E, Braks A, van de Vegte-Bolmer M, Qiu LY, Sauerwein R, Stunnenberg HG 2007. Epitope analysis of the malaria surface antigen pfs48/45 identifies a subdomain that elicits transmission blocking antibodies. J Biol Chem 282: 17148-17156. [ Links ]

Outchkourov NS, Roeffen W, Kaan A, Jansen J, Luty A, Schuiffel D, van Gemert GJ, van de Vegte-Bolmer M, Sauerwein RW, Stunnenberg HG 2008. Correctly folded Pfs48/45 protein of Plasmodium falciparum elicits malaria transmission-blocking immunity in mice. Proc Natl Acad Sci USA 105: 4301-4305. [ Links ]

Pradel G 2007. Proteins of the malaria parasite sexual stages: expression, function and potential for transmission blocking strategies . Parasitology 134: 1911-1929. [ Links ]

Roeffen W, Lensen T, Mulder B, Teelen K, Sauerwein R, Van Druten J, Eling W, Meuwissen JH, Beckers PJ 1995. A comparison of transmission-blocking activity with reactivity in a Plasmodium falciparum 48/45-kD molecule-specific competition enzyme-linked immunosorbent assay. Am J Trop Med Hyg 52: 60-65. [ Links ]

Roeffen W, Teelen K, van As J, vd Vegte-Bolmer M, Eling W, Sauerwein R 2001. Plasmodium falciparum: production and characterization of rat monoclonal antibodies specific for the sexual-stage Pfs48/45 antigen. Exp Parasitol 97: 45-49. [ Links ]

Stowers A, Carter R 2001. Current developments in malaria transmission-blocking vaccines. Expert Opin Biol Ther 1: 619-628. [ Links ]

Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S 2011. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance and maximum parsimony methods. Mol Biol Evol 28: 2731-2739. [ Links ]

van Dijk MR, Janse CJ, Thompson J, Waters AP, Braks JA, Dodemont HJ, Stunnenberg HG, van Gemert GJ, Sauerwein RW, Eling W 2001. A central role for P48/45 in malaria parasite male gamete fertility. Cell 104: 153-164. [ Links ]

van Schaijk BC, van Dijk MR, van de Vegte-Bolmer M, van Gemert GJ, van Dooren MW, Eksi S, Roeffen WF, Janse CJ, Waters AP, Sauerwein RW 2006. Pfs47, paralog of the male fertility factor Pfs48/45, is a female specific surface protein in Plasmodium falciparum. Mol Biochem Parasitol 149: 216-222. [ Links ]

WHO - World Health Organization 2011. WHO global malaria programme. World malaria report 2011. Available from: who.int/malaria/world_malaria_report_2011/9789241564403_eng.pdf. [ Links ]

Yeom JS, Ryu SH, Oh S, Lee WJ, Kim TS, Kim KH, Kim YA, Ahn SY, Cha JE, Park JW 2005. Status of Plasmodium vivax malaria in the republic of Korea during 2001-2003. Am J Trop Med Hyg 73: 604-608. [ Links ]

Korea University, R & D Grant (2011)

Received: August 09, 2012; Accepted: October 09, 2012

+ Corresponding author: malarim@korea.ac.kr