<title>Background</title> <p>Molecular mechanisms of invasion by <italic>Plasmodium vivax </italic>merozoites are mediated by the Duffy-binding-like (DBL) family of homologous erythrocyte binding protein (EBPs) located within the micronemes of merozoites that recognize specific receptors on red blood cell [##REF##1496004##1##]. <italic>PvDBP </italic>is a 140-kDa protein belonging to a family of erythrocyte binding proteins characterized by a functionally conserved cysteine-rich region. The similarity among DBL families of EBPs is the greatest within two adjacent cysteine-rich domains, designated as the amino cysteine-rich domain (N-cys) and carboxyl cysteine-rich domain (C-cys) [##REF##1849231##2##]. The amino cysteine-rich domain has been identified as the principle adhesion region binding to erythrocyte receptors, while the function of the carboxyl cysteine-rich domain remains unclear. This amino-cysteine-rich domain is present within the 330 amino acid region II, which has been shown to contain the binding motifs necessary for adherence to the Duffy antigen receptor for chemokines (DARC), required for erythrocyte invasion [##REF##8879225##3##]. Critical binding motifs in <italic>PvDBPII </italic>have been mapped to a 170 amino acid spanning region between cysteines 4 and 7 (amino acids 291–460) [##REF##10570199##4##]. Cysteine residues are conserved within the identified binding motif, whereas other amino acids are highly polymorphic, having a high ratio of nonsynonymous to synonymous mutations [##REF##7960140##5##].</p> <p>The polymorphic nature of the <italic>PvDBP</italic>, particularly its region II (<italic>DBPII</italic>) is a major impediment to the development of a broadly protective vaccine against vivax malaria. A study in Papua New Guinea showed 43 and 18 unique nonsynonymous mutations in 10 Madang [##REF##12195381##6##] and 24 Wosera vivax isolates [##REF##11163434##7##], respectively, whereas 19 unique nonsynonymous mutations were described in 10 Colombian vivax isolates [##REF##8813697##8##]. Often the mutations were found in corresponding locations in the protein sequence of isolates from different geographical areas. Amino acids substitutions at D384G, K386(N/Q), N417K, L424I, W437R, and I503K were found in both Papua New Guinea isolates and Colombia isolates [##REF##7960140##5##, ####REF##12195381##6##, ##REF##11163434##7##, ##REF##8813697##8####8813697##8##]. Such results demonstrate that many similar alleles are widely distributed among <italic>P. vivax </italic>from different geographical areas. Polymorphic patterns in the <italic>PvDBP </italic>cluster within the critical erythrocyte-binding segment imply the existence of a selection pressure on <italic>PvDBP </italic>[##REF##12672521##9##]. Under the influence of immune pressure, the distribution of existing antigenic polymorphisms affects the population structure in which newly selected mutants spread and may provide insight to the evolution and selection of parasite populations over time [##REF##9371796##10##]. Polymorphisms of <italic>PvDBPII </italic>require further exploration among vivax isolates from different geographic areas, particularly south-east Asia where a large proportion of vivax infections take place. Therefore, the gene polymorphisms within the region II domain of the <italic>P. vivax DBP </italic>among Thai <italic>P. vivax </italic>isolates were investigated.</p>
<title>Methods</title> <title><italic>P. vivax </italic>samples and DNA preparation</title> <p><italic>Plasmodium vivax </italic>isolates were obtained from malaria patients with informed consent at the outpatient clinic, Hospital for Tropical Diseases, Faculty of Tropical Medicine, Mahidol University, Bangkok, from May, 2002 to June, 2003 under a protocol approval by the Ethics Committee, Faculty of Tropical Medicine, Mahidol University. The confirmation of <italic>P. vivax </italic>infection was done by microscopic examination of the parasites in thick and thin blood smears. Approximately 1 ml of venous blood was collected from each individual in 0.5 M ethylene-diamine tetra-acetic acid (EDTA) tubes and kept at -20°C until use. Genomic DNA was extracted from 200 μl of blood using a High Pure PCR Template Preparation Kit (Roche, Mannhein, Germany) according to the manufacturer's protocol.</p> <title>Genotyping of <italic>P. vivax </italic>isolates at the <italic>PvDBPII </italic>locus</title> <p><italic>P. vivax DBPII </italic>gene was amplified by PCR using the oligonucleotide primers, DBPII-F; 5'-CACCACGATCTCTAGTGCTATTATA-3' and DBPII-R; 5'-TGTCACAACTTCCTGAGTATT-3' [##REF##10960177##11##]. PCR cycling conditions were 94°C for 3 min, 55°C for 1 min, and 72°C for 1 min for 35 cycles followed by a 10 min extension at 72°C. PCR products were stained with 0.1 μg/ml ethidium bromide and visualized by digital photography under ultraviolet light. The PCR products were run adjacently to 0.1 μg/lane of a 2-log DNA ladder (NewEngland BioLabs, Inc. MA, USA) as a standard DNA marker.</p> <title>Purification of PCR products, DNA cloning and DNA sequencing of <italic>PvDBPII </italic>gene</title> <p>PCR products were purified by QIAquick PCR purification Kit (Qiagen Inc., Valencia, CA) following the manufacturer's instructions. The PCR products were cloned directly into pCR-Blunt II-TOPO cloning vector using a Zero Blunt TOPO PCR cloning kit (Invitrogen, La Jolla, CA). For DNA sequencing of the species-specific amplicons, excess dNTPs and unincorporated primers were eliminated using the QIAquick PCR purification kit (Qiagen Inc., Valencia, CA). DNA sequencing was performed by fluorescence based methodologies using MegaBACE DNA Analysis Systems (Amersham Biosciences, Piscataway, NJ). The sequences have been submitted to the GenBank under the accession number <ext-link ext-link-type="gen" xlink:href="EF219451">EF219451</ext-link>, <ext-link ext-link-type="gen" xlink:href="EF368159">EF368159</ext-link>–<ext-link ext-link-type="gen" xlink:href="368180">368180</ext-link>, <ext-link ext-link-type="gen" xlink:href="EF379127">EF379127</ext-link>–<ext-link ext-link-type="gen" xlink:href="379135">379135</ext-link>. Due to the possibility of multiple genotypes within parasites taken from one individual, direct sequencing was not performed. Instead, at least five clones were sequenced from each corresponding PCR product and sequences of these had to be identical to meet our quality control standards.</p> <title>Analysis of <italic>PvDBPII </italic>gene sequences</title> <p>The alignment of complete sequences of 30 <italic>PvDBPII </italic>genes was analysed by Sequencher ver 4.2 software (Gene Code Corporation, Ann Arbor, MI) [##REF##1630306##12##,##REF##2513255##13##]. Sequences were aligned using CLUSTALX Multiple Sequence Alignment Program developed by the National Center for Biotechnology Information (NCBI, Bethesda, MD). The percent similarity was assessed using BLASTN (NCBI, Bethesda, MD) and BioEdit ver 7.0 software (Tom Hall Isis Therapeutics, Isis Pharmaceuticals, CA), DnaSP ver 4.0 [##REF##14668244##14##] and MEGA ver 4.0 Beta program [##REF##15260895##15##].</p> <title>Analysis of natural selection</title> <p>Evidence of positive natural selection was determined by comparing the rate of nonsynonymous and synonymous substitutions in region II of the <italic>P. vivax DBP </italic>gene. The rates of substitutions were estimated using the Nei and Gojobori's method [##REF##3444411##16##] with the Jukes and Cantor correction [##REF##10757745##17##,##REF##1346783##18##] as implemented in the MEGA program. Standard error was determined by 1000 bootstrap replications. Negative selection acting on most coding genes can be identified typically when the rate on nonsynonymous mutations is less than the rate of synonymous mutations (<italic>k</italic><sub>n </sub>&lt;<italic>K</italic><sub>s</sub>). However, when positive selection is acting on a gene, the rate of nonsynonymous mutations will exceed that of synonymous mutations (<italic>k</italic><sub>n </sub>&gt; <italic>K</italic><sub>s</sub>). Two other tests of neutrality were performed by Tajima's D test, Fu and Li's <italic>D- </italic>and <italic>F</italic>-tests on DnaSP 4 software using <italic>P. knowlesi </italic>DBP as an outgroup. Tajima's D test compares the estimation of nucleotide diversity calculated in two ways (θ calculated from the number of segregating sites and π calculated from average pairwise nucleotide diversity) in order to test for a departure from neutrality. Fu and Li's <italic>D </italic>and <italic>F </italic>test, departures from neutrality are identified as a deviation between estimates of θ (derived from the number of mutations in external branches of the phylogeny and from the total number of mutations giving the index <italic>D </italic>or from the average pairwise diversity π giving the index <italic>F</italic>) [##REF##8454210##19##].</p> <title>Phylogenetic analysis</title> <p>Phylogenetic analysis was used to investigate the associations of <italic>PvDBPII </italic>gene with sequences elucidated from different geographical regions. The gene tree was constructed using regions common to all available <italic>PvDBPII </italic>sequences. Forty five <italic>PvDBPII </italic>genes sequences found in GenBank were compiled and compared to Thai isolates including the sequences from a reference strain, Sal I, one individual sequence from Vietnam, Indonesia, Brazil and India as well as 2, 3, 3 and 31 sequences from Bangladesh, Colombia (COL), South Korea (SK) and Papua New Guinea (PNG), respectively (Table ##TAB##0##1##). A phylogenetic tree was derived from the aligned nucleotide sequences using the Neighbour-Joining (NJ) method with 1000 bootstrap replicates, the Tamura's three-parameter distance model as implemented in the MEGA ver 4.0 Beta program [##REF##9162471##20##].</p>
<title>Results</title> <title>PCR products and genetic polymorphisms of <italic>PvDBPII </italic>gene among Thai isolates</title> <p>The PCR amplification of the <italic>P. vivax DBPII </italic>gene from 30 Thai isolates yielded DNA fragments of approximately 990 base pairs in length (Figure ##FIG##0##1##). A 900 bp was selected after excluding primer regions and adjacent bases which did not give uniformly reliable sequences. DNA analysis in comparison to the <italic>PvDBP </italic>gene of Sal I reference strain revealed amino acid polymorphisms across the entire domain II of <italic>PvDBP </italic>protein (300 amino acids) among Thai isolates. Among the 30 polymorphic sites, 29 showed monomorphic mutation (changed into one amino acid-type) and only a position 386(N/Q) showed dimorphic mutation (changed into two amino acid-types) (Table ##TAB##1##2##). Among 300 amino acids within <italic>PvDBPII</italic>, 25 nonsynonymous (representing 83% of all isolates) and five synonymous mutations were identified (Table ##TAB##1##2##). High frequencies of variant amino acids (&gt;50%) were found and included L424I (86%), D384G (76%), W437R (63%), R390H (56%), and I503K (56%) residue polymorphisms.</p> <p>Based on the variation of nonsynonymous mutations and in comparison to Sal I the Thai haplotyes were designated as TA (TA1-TA25). TA1, TA5 and TA11 haplotypes were identified in two isolates while TA17 was identified in three isolates. The other 21 haplotypes were identified in single isolates (Table ##TAB##2##3##). The cluster analysis of <italic>PvDBPII </italic>amino acid sequences among 25 haplotypes based on Neighbour-Joining method was organized into five main groups (Figure ##FIG##1##2##).</p> <title>Evidence of natural selection</title> <p>To investigate whether natural selection contributed to generation of this diversity, the ratio of nonsynonymous (<italic>K</italic><sub>n</sub>) to synonymous substitutions (<italic>K</italic><sub>s</sub>) was compared between Thai and Sal I isolates. The number of nonsynonymous substitutions per nonsynonymous sites (<italic>K</italic><sub>n</sub>) (0.00956) exceeded the number of synonymous substitutions per synonymous sites (<italic>K</italic><sub>s</sub>) (0.00754). The <italic>K</italic><sub>n</sub>/<italic>K</italic><sub>s </sub>ratio was 1.270 predicting that positive selection may be occurring in region II of <italic>DBP </italic>favouring the fixation of amino acid replacement in certain areas of the protein. The statistic tests detecting departure from neutrality showed evidence that the observed polymorphism might be trending natural selection. Tajima's D value was -0.084549 with P &gt; 0.10. The Fu and Li, <italic>D </italic>and <italic>F </italic>values were -0.92031 with P &gt; 0.10 and -0.7654 with P &gt; 0.10, respectively.</p> <title>Phylogenetic tree of <italic>PvDBPII</italic></title> <p>Phylogenetic tree analysis of DNA sequences derived from region II of the <italic>P. vivax DBP </italic>constructs based on the Neighbour-Joining method using Tamura's three-parameter distance is shown in Figure ##FIG##2##3##. Variation of DNA sequence was categorized into nine distinct alleles groups. Thai isolates were widely distributed amongst groups 1, 3, 5 and 7 also containing isolates from different geographical regions. <italic>PvDBPII </italic>DNA sequence of Thai isolates was related to isolates sets from Papua New Guinea and India (group 1), Korea (group 3), Colombia (group 5) and a second Papua New Guinea set (group 7). The Sal I reference strain used in this study was included in group 2 composed of isolates from Colombia, Korea, Bangladesh, Indonesia, and Vietnam, but bootstrap analysis demonstrated that the Thai isolates were distinct from this group. Most isolates from Papua New Guinea formed a distinctive group comprised of 8 and 9.</p>
<title>Discussion</title> <p>The process invasion into erythrocyte is essential for survival of <italic>P. vivax </italic>parasite and requires <italic>DBP </italic>to bind the Duffy blood group antigen that act as erythrocyte receptor, Duffy antigen receptor for chemokines (DARC). <italic>Plasmodium vivax </italic>utilizes genetic diversity of individual functional proteins to evade the host immune system and resist many anti-malarial drugs [##REF##15478059##21##]. In this study, polymorphisms and natural selection was assessed and suggested positive natural selection in <italic>PvDBPII </italic>of Thai isolates. Although, non significant could be observed trending positive selection. As this is only a glimpse at the overall population, a wide range study including more subjects from different areas of Thailand would give a more detailed picture.</p> <p>Analysis of 300 amino acids of <italic>PvDBPII </italic>protein relative to the Sal I reference strain, resulted in 25 nonsynonymous mutations being found. In addition, 142 nonsynonymous mutations were reported in 76 Papua New Guinea isolates [##REF##17984360##22##]. The nonsynonymous mutations in <italic>PvDBPII </italic>might have an effect on parasite binding to erythrocyte receptors. Nonsynonymous mutations often times result in amino acid substitutions that alter charge and therefore host immune/antibody recognition [##REF##12195381##6##]. Thus, excessive stable polymorphism generation within the <italic>PvDBPII </italic>ligand domain from high rate of nonsynonymous polymorphisms may promote the parasite escape of host immune response [##REF##12672521##9##].</p> <p>The high frequency (&gt;50%) of L424I (86%), D384G (76%), W437R (63%), R390H (56%), and I503K (56%) residues were found in Thai isolates relative to Sal I <italic>PvDBP </italic>sequence. The finding is in contrast to previous studies showing R308S (67%), D384G (66%), and S447K (59%) in Papua New Guinea, D384G (59%) in Colombia and D384G (85%) in Brazil isolates [##REF##16384615##23##]. Of the following variants N417K, W437R and I503K, all having been shown to be involved in evasion to antibody neutralization, variant N417K was found in low frequency (36%) while W437R and I503K were high (63% and 56% respectively) in Thai isolates. These results were different from what has been currently reported in that the low frequencies of variants N417K (27.5%) and W437R (27.5%) and the high frequency of variant I503K (55%) found in Brazil and low frequencies of all three variants (N417K, W437R, I503K) found in either Colombia (47%, 18%, 12%, respectively) and/or Papua New Guinea (23%, 26%, 29%, respectively) (Table ##TAB##3##4##) [##REF##16384615##23##]. N417K, W437R and I503K have recently been shown to not be directly involved in erythrocyte binding, but possibly in inhibition of antibody binding to erythrocyte receptor [##REF##15498870##24##]. Moreover, the variant amino acid T404R was identified in <italic>PvDBPII </italic>domains of three Thai isolates (Table ##TAB##1##2##) which is similar to the recent data in Papua New Guinea [##REF##17984360##22##]. Other Thai variants, namely R308S, K371E, D384G, E385K, N386(N/Q), R390H, T404R, N417K, L424I, W437R, S447S and I503K, were also found in isolates of other regions, demonstrating commonality with a progenitor clone (Table ##TAB##3##4##). The other 18 polymorphic residues of <italic>PvDBPII </italic>ligand domain were recognized only among Thai isolates, indicative of divergence from a progenitor clone.</p> <p>The study showed that 25 <italic>PvDBPII </italic>haplotypes clustered into five main groups (dominant haplotypes) (Table ##TAB##2##3##) on the basis of nonsynonymous mutations is similar, but not entirely the same, with data from Papua New Guinea isolates whose 27 different <italic>PvDBPII </italic>haplotypes clustering into three dominant <italic>PvDBPII </italic>haplotypes. Potentially dominant <italic>PvDBPII </italic>haplotypes observed in these populations could be less immunogenic, thereby may help promoting parasite success and escaping the host immune response [##REF##17984360##22##].</p> <p>Phylogenetic analysis of <italic>PvDBPII </italic>suggested that Thai isolates fell into 6 different alleles groups. Among these, group 1 formed a group with a subset of Papua New Guinea isolates. Groups 4 and 6 were unique among Thai <italic>P. vivax </italic>isolates. However, the rest of the allele groups that Thai isolates fell into were more related to the isolates from Korea, India, and Colombia than from Papua New Guinea (Figure ##FIG##2##3##). The observations derived from the <italic>PvDBPII </italic>phylogeny suggested that single Indian isolate found in group 1 appeared to share a common ancestor close to a subset of Thai and Papua New Guinea isolates, demonstrating common ancestoral origins. Thai and Papua New Guinea isolate subsets in group 7 also may have shared common ancestor heritage. Recently, it was reported that travellers' malaria among foreigners at the Hospital for Tropical Diseases, Bangkok, Thailand including India and Papua New Guinea patients were though to have acquired their infections in their countries. Malaria importations from country to country can occur by either immigration or travel, and changing malaria attack rates in the countries of exposure are likely to influence the incidence of imported disease [##REF##16969060##25##]. Additionally, these results demonstrated that there is some difference in the nucleotide or amino acid variation in <italic>PvDBPII </italic>among isolates from Thai, Korea, India, Papua New Guinea and Colombia, but most importantly, unique variants from Thailand do exist and should be considered in future vaccine development of <italic>PvDBPII</italic>. It should be noted that the present highly diverse phylogenetic tree was constructed relying on polymorphisms within a single gene, demonstrating the high degree of topological variation that can be found using this technique, but also highlighting the problems that can be associated with gene sampling in phylogenetic studies [##REF##17567929##26##]. Use of a combination of several genes with different rates of evolutions has been suggested as an efficient way to overcome this difficulty to prove the existence of a unique tree relating these sequences [##REF##16387305##27##]. Finally, it must also be pointed out that in many situations a single-gene phylogeny may be interesting in itself. Awareness of the problems of orthology (genes in different species that derive from a common ancestor) assignment and tree reconstruction artifacts should be considered.</p>
<title>Conclusion</title> <p>The results indicated that <italic>PvDBPII </italic>gene among Thai isolates is genetically diverse. Nonsynonymous polymorphism in <italic>PvDBPII </italic>was the predominant result of which variant L424I showed the highest frequency (86%). Five dominant <italic>PvDBPII </italic>haplotypes could be clustered among Thai isolates. The high frequency of polymorphisms and the presence of distinct alleles within <italic>PvDBPII </italic>gene among <italic>P. vivax </italic>from different geographical areas could provide complication in malaria vaccine development, but should be considered important parts of the development process.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>The Duffy-binding protein II of <italic>Plasmodium vivax </italic>(<italic>PvDBPII</italic>) has been considered as an attractive target for vaccine-mediated immunity despite a possible highly polymorphic nature. Among seven <italic>PvDBP </italic>domains, domain II has been shown to exhibit a high rate of nonsynonymous polymorphism, which has been suggested to be a potential immune (antibody binding) evasion mechanism. This study aimed to determine the extent of genetic polymorphisms and positive natural selection at domain II of the <italic>PvDBP </italic>gene among a sampling of Thai <italic>P. vivax </italic>isolates.</p> <title>Methods</title> <p>The <italic>PvDBPII </italic>gene was PCR amplified and the patterns of polymorphisms were characterized from 30 Thai <italic>P. vivax </italic>isolates using DNA cloning and sequencing. Phylogenetic analysis of the sequences and positive selection were done using DnaSP ver 4.0 and MEGA ver 4.0 packages.</p> <title>Results</title> <p>This study demonstrated a high rate of nonsynonymous polymorphism. Using Sal I as the reference strain, a total of 30 point-mutations were observed in the <italic>PvDBPII </italic>gene among the set of Thai <italic>P. vivax </italic>isolates, of which 25 nonsynonymous and five synonymous were found. The highest frequency of polymorphism was found in five variant amino acids (residues D384G, R390H, L424I, W437R, I503K) with the variant L424I having the highest frequency. The difference between the rates of nonsynonymous and synonymous mutations estimated by the Nei and Gojobori's method suggested that <italic>PvDBPII </italic>antigen appears to be under selective pressure. Phylogenetic analysis of <italic>PvDBPII </italic>Thai <italic>P. vivax </italic>isolates to others found internationally demonstrated six distinct allele groups. Allele groups 4 and 6 were unique to Thailand.</p> <title>Conclusion</title> <p>Polymorphisms within <italic>PvDBPII </italic>indicated that Thai vivax malaria parasites are genetically diverse. Phylogenetic analysis of DNA sequences using the Neighbour-Joining method demonstrated that Thai isolates shared distinct alleles with <italic>P. vivax </italic>isolates from different geographical areas. The study reported here will be valuable for the development of <italic>PvDBPII</italic>-based malaria vaccine.</p>
<title>List of abbreviations used</title> <p><italic>PvDBPII</italic>: Domain II of <italic>Plasmodium vivax DBP </italic>gene; DARC: Duffy antigen receptor for chemokines</p> <title>Authors' contributions</title> <p>SRM and SK designed the study, SK, MMF, SRM, KES, DEL and TK, were responsible for the supervision of the work, KES, TK, DEL, and PG were responsible for the sequences and phylogenetic analysis, SK, MMF, SRM, KES, DEL, TK, and PG drafted the manuscript. All authors read and approved the final manuscript.</p>
<title>Acknowledgements</title> <p>We thank the patients for their kind participation in this study; staff of the Hospital for Tropical Diseases, Faculty of Tropical Medicine, Mahidol University for their assistance in blood collection and COL Douglas S. Walsh, US Army Medical Research Unit-Kenya for his valuable guidance and comments. This study was supported by the United States Army Medical Research and Materiel Command, Fort Detrick, MD and Faculty of Tropical Medicine, Mahidol University.</p> <p>Disclaimer</p> <p>The views of the authors are solely there own, and do not reflect official opinion of the United States Army Medical Research and Materiel Command or the Department of Defense.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p>Schematic diagram of <italic>PvDBPII </italic>gene showing the location of the primers.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p>Neighbour-Joining tree of <italic>PvDBPII </italic>amino acid sequence of 25 variants observed among 30 Thai isolates. The <italic>PvDBPII </italic>haplotype for each variant could be seen in Table 3.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p>Phylogenetic tree of <italic>PvDBPII </italic>obtaining from Thai <italic>P. vivax </italic>isolates and other isolates from different geographical malaria endemic areas, <italic>i.e</italic>. PNG (Papua New Guinea), India, COL (Colombia), Korea, Bangladesh, Brazil, Indonesia, Vietnam, and Sal I. Numbers at nodes indicate percentage support of 1000 bootstrap replicates (bootstrap support values below 10 are not presented).</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p><italic>PvDBPII </italic>sequences deposited in GenBank used in the study</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">species</td><td align="center">locality</td><td align="left">GenBank accession number</td></tr></thead><tbody><tr><td align="left"><italic>P. vivax</italic></td><td align="center">Papua New Guinea</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AF289480">AF289480</ext-link>–<ext-link ext-link-type="gen" xlink:href="289653">289653</ext-link>, <ext-link ext-link-type="gen" xlink:href="AF469522">AF469522</ext-link>–<ext-link ext-link-type="gen" xlink:href="469602">469602</ext-link>, <ext-link ext-link-type="gen" xlink:href="AY970848">AY970848</ext-link>–<ext-link ext-link-type="gen" xlink:href="970925">970925</ext-link>, <ext-link ext-link-type="gen" xlink:href="AF291096">AF291096</ext-link>, <ext-link ext-link-type="gen" xlink:href="AF695565">AF695565</ext-link>, <ext-link ext-link-type="gen" xlink:href="DQ156519">DQ156519</ext-link></td></tr><tr><td align="left"><italic>P. vivax</italic></td><td align="center">Colombia</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AY341907">AY341907</ext-link>, <ext-link ext-link-type="gen" xlink:href="DQ156513">DQ156513</ext-link>, <ext-link ext-link-type="gen" xlink:href="AY341899">AY341899</ext-link></td></tr><tr><td align="left"><italic>P. vivax</italic></td><td align="center">Korea</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AF215737">AF215737</ext-link>, <ext-link ext-link-type="gen" xlink:href="AF215738">AF215738</ext-link>, <ext-link ext-link-type="gen" xlink:href="DQ156523">DQ156523</ext-link></td></tr><tr><td align="left"><italic>P. vivax</italic></td><td align="center">India</td><td align="left"><ext-link ext-link-type="gen" xlink:href="DQ156514">DQ156514</ext-link></td></tr><tr><td align="left"><italic>P. vivax</italic></td><td align="center">Brazil</td><td align="left"><ext-link ext-link-type="gen" xlink:href="DQ153520">DQ153520</ext-link></td></tr><tr><td align="left"><italic>P. vivax</italic></td><td align="center">Vietnam</td><td align="left"><ext-link ext-link-type="gen" xlink:href="DQ156518">DQ156518</ext-link></td></tr><tr><td align="left"><italic>P. vivax</italic></td><td align="center">Indonesia</td><td align="left"><ext-link ext-link-type="gen" xlink:href="DQ156521">DQ156521</ext-link></td></tr><tr><td align="left"><italic>P. vivax</italic></td><td align="center">Sal I</td><td align="left"><ext-link ext-link-type="gen" xlink:href="M37514">M37514</ext-link>, <ext-link ext-link-type="gen" xlink:href="DQ156512">DQ156512</ext-link></td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Amino acid changes found in <italic>PvDBPII </italic>gene among 30 Thai isolates, comparing to <italic>DBP </italic>Sal I sequence.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center" colspan="5">Position of the amino acid</td><td/><td/><td/><td/><td/><td/><td/><td align="center">▼</td><td/><td align="center">▽</td><td align="center">▼</td></tr></thead><tbody><tr><td/><td align="center">268</td><td align="center">276</td><td align="center">281</td><td align="center">306</td><td align="center">308</td><td align="center">333</td><td align="center">351</td><td align="center">367</td><td align="center">371</td><td align="center">375</td><td align="center">378</td><td align="center">384</td><td align="center">385</td><td align="center">386</td><td align="center">390</td></tr><tr><td align="center">SalI</td><td align="center">CGT</td><td align="center">AAC</td><td align="center">GTT</td><td align="center">TTT</td><td align="center">AGG</td><td align="center">CTT</td><td align="center">AGT</td><td align="center">ATC</td><td align="center">AAA</td><td align="center">AAT</td><td align="center">CGC</td><td align="center">GAT</td><td align="center">GAA</td><td align="center">AAG</td><td align="center">CGT</td></tr><tr><td align="center">Thai</td><td align="center">AGT</td><td align="center">AAT</td><td align="center">GTA</td><td align="center">TTG</td><td align="center">AGT</td><td align="center">TTT</td><td align="center">TGT</td><td align="center">ACC</td><td align="center">GAA</td><td align="center">GAT</td><td align="center">CGT</td><td align="center">GGT</td><td align="center">AAA</td><td align="center">AAT/CAG</td><td align="center">CAT</td></tr><tr><td align="center">SalI</td><td align="center">R</td><td align="center">N</td><td align="center">V</td><td align="center">F</td><td align="center">R</td><td align="center">L</td><td align="center">S</td><td align="center">I</td><td align="center">K</td><td align="center">N</td><td align="center">R</td><td align="center">D</td><td align="center">E</td><td align="center">K</td><td align="center">R</td></tr><tr><td align="center">Thai</td><td align="center">S</td><td align="center">N</td><td align="center">V</td><td align="center">L</td><td align="center">S</td><td align="center">F</td><td align="center">C</td><td align="center">T</td><td align="center">E</td><td align="center">D</td><td align="center">R</td><td align="center">G</td><td align="center">K</td><td align="center">N/Q</td><td align="center">H</td></tr><tr><td align="center">a</td><td align="center">(1/30)</td><td align="center">(1/30)</td><td align="center">(1/30)</td><td align="center">(2/30)</td><td align="center">(8/30)</td><td align="center">(14/30)</td><td align="center">(1/30)</td><td align="center">(2/30)</td><td align="center">(6/30)</td><td align="center">(7/30)</td><td align="center">(30/30)</td><td align="center">(23/30)</td><td align="center">(14/30)</td><td align="center">(13/30)</td><td align="center">(17/30)</td></tr><tr><td align="center">b</td><td align="center">3</td><td align="center">3</td><td align="center">3</td><td align="center">6</td><td align="center">26</td><td align="center">46</td><td align="center">3</td><td align="center">6</td><td align="center">20</td><td align="center">23</td><td align="center">100</td><td align="center">76</td><td align="center">46</td><td align="center">43</td><td align="center">56</td></tr><tr><td colspan="16"><hr/></td></tr><tr><td align="center" colspan="5">Position of the amino acid</td><td align="center">▼</td><td/><td/><td align="center">▼</td><td/><td/><td/><td/><td align="center">▼</td><td/><td/></tr><tr><td colspan="16"><hr/></td></tr><tr><td/><td align="center">398</td><td align="center">404</td><td align="center">417</td><td align="center">419</td><td align="center">424</td><td align="center">433</td><td align="center">436</td><td align="center">437</td><td align="center">447</td><td align="center">464</td><td align="center">475</td><td align="center">486</td><td align="center">503</td><td align="center">507</td><td align="center">513</td></tr><tr><td align="center">Sal I</td><td align="center">TCT</td><td align="center">ACA</td><td align="center">AAT</td><td align="center">ATA</td><td align="center">TTA</td><td align="center">CAG</td><td align="center">AGA</td><td align="center">TGG</td><td align="center">TCA</td><td align="center">ATC</td><td align="center">CCA</td><td align="center">CAA</td><td align="center">ATA</td><td align="center">AAC</td><td align="center">ACG</td></tr><tr><td align="center">Thai</td><td align="center">ACT</td><td align="center">AGA</td><td align="center">AAA</td><td align="center">ATG</td><td align="center">ATA</td><td align="center">AAG</td><td align="center">ACA</td><td align="center">CGG</td><td align="center">TCC</td><td align="center">ATA</td><td align="center">GCA</td><td align="center">GAA</td><td align="center">AAA</td><td align="center">CAC</td><td align="center">AAG</td></tr><tr><td align="center">Sal I</td><td align="center">S</td><td align="center">T</td><td align="center">N</td><td align="center">I</td><td align="center">L</td><td align="center">Q</td><td align="center">R</td><td align="center">W</td><td align="center">S</td><td align="center">I</td><td align="center">P</td><td align="center">Q</td><td align="center">I</td><td align="center">N</td><td align="center">T</td></tr><tr><td align="center">Thai</td><td align="center">T</td><td align="center">R</td><td align="center">K</td><td align="center">M</td><td align="center">I</td><td align="center">K</td><td align="center">T</td><td align="center">R</td><td align="center">S</td><td align="center">I</td><td align="center">A</td><td align="center">E</td><td align="center">K</td><td align="center">H</td><td align="center">K</td></tr><tr><td align="center">a</td><td align="center">(3/30)</td><td align="center">(3/30)</td><td align="center">(11/30)</td><td align="center">(1/30)</td><td align="center">(26/30)</td><td align="center">(1/30)</td><td align="center">(1/30)</td><td align="center">(19/30)</td><td align="center">(1/30)</td><td align="center">(1/30)</td><td align="center">(3/30)</td><td align="center">(2/30)</td><td align="center">(17/30)</td><td align="center">(1/30)</td><td align="center">(1/30)</td></tr><tr><td align="center">b</td><td align="center">10</td><td align="center">10</td><td align="center">36</td><td align="center">3</td><td align="center">86</td><td align="center">3</td><td align="center">3</td><td align="center">63</td><td align="center">3</td><td align="center">3</td><td align="center">10</td><td align="center">6</td><td align="center">56</td><td align="center">3</td><td align="center">3</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p><italic>PvDBPII </italic>haplotypes observed among 30 Thai isolates.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">Haplotype name</td><td align="center">Amino acid haplotype*</td><td align="center">No. of isolate observed (%)</td></tr></thead><tbody><tr><td align="center">TA 1</td><td align="center">RFRFSIKDGKNHSTKIIQRRPQKNT</td><td align="center">2(6.6)</td></tr><tr><td align="center">TA 2</td><td align="center">RFSFSIKDGKNHSTKIIQRRPQKNT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA3</td><td align="center">RFSFSIKDGKKHSTKIIQRRPQKNT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 4</td><td align="center">RFRLSIKDGKNHSTKIIQRRPQKNT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 5</td><td align="center">RFRFSIENGKNRSTNIIQRRPQINT</td><td align="center">2(6.6)</td></tr><tr><td align="center">TA6</td><td align="center">RFRFSIENGKNRSTKIIQRRPQINT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 7</td><td align="center">RFSFSIKDGKNHSTNIIQRRPQINT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 8</td><td align="center">RFRLSIENGEKHSTKIIQRRPQKNT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 9</td><td align="center">RFRLSIENGEQRSTKIIQRRPQKNT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 10</td><td align="center">RFRLSIENGEKRSTKMIQRRPQINT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 11</td><td align="center">RLSLLSIKNGKNHTRNIIQRWPQINT</td><td align="center">2(6.6)</td></tr><tr><td align="center">TA 12</td><td align="center">SFSLSIKNGKNHTRNIQRWPEINT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 13</td><td align="center">RFRLSTKNGEKHSTNIIQRWPQKHT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 14</td><td align="center">RFRLSTKNGEKHSTNIIQRWPQKHT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 15</td><td align="center">RFSLSIKNGEKRSTNIIQRWPQKNT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 16</td><td align="center">RFRLSIKNGEKRSTNIIQRWPQINT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 17</td><td align="center">RFRFSIKNDEKRSTNIIQRRAQKNT</td><td align="center">3(10)</td></tr><tr><td align="center">TA 18</td><td align="center">RFRLSIKNDEKRSTNIIQRRPQKNT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA19</td><td align="center">RFRFSIKNDKKRSTNIIQRRPQKNT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 20</td><td align="center">RFRLSIKNDKKRSTKIIKTRPQINT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 21</td><td align="center">RFRLSIKNDKKRSTNIIQRWPQINT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 22</td><td align="center">RFSFSIKDGKNHSTNILQRRPQKNT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA23</td><td align="center">RFRLSIKNGEKHSTNILQRWPQKNK</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 24</td><td align="center">RFRLCIKNGEKHSTNILQRWPEINT</td><td align="center">1(3.3)</td></tr><tr><td align="center">TA 25</td><td align="center">RFRLSIKNGEKHSTNILQRWPQINT</td><td align="center">1(3.3)</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>Frequencies of the most common variant amino acids in <italic>PvDBPII</italic>, comparing to Sal I sequence (accession number <ext-link ext-link-type="gen" xlink:href="M37519">M37519</ext-link>)</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">#</td><td align="center">R308S</td><td align="center">D384G</td><td align="center">K386N</td><td align="center">K386Q</td><td align="center">N417K</td><td align="center">L424I</td><td align="center">W437R</td><td align="center">S447K</td><td align="center">I503K</td></tr></thead><tbody><tr><td align="center">Thai</td><td align="center">26</td><td align="center">76</td><td align="center">40</td><td align="center">3</td><td align="center">36</td><td align="center">86</td><td align="center">63</td><td align="center">0</td><td align="center">56</td></tr><tr><td align="center">COL*</td><td align="center">0</td><td align="center">59</td><td align="center">23</td><td align="center">0</td><td align="center">47</td><td align="center">47</td><td align="center">18</td><td align="center">0</td><td align="center">12</td></tr><tr><td align="center">PNG*</td><td align="center">67</td><td align="center">66</td><td align="center">8</td><td align="center">11</td><td align="center">23</td><td align="center">34</td><td align="center">26</td><td align="center">59</td><td align="center">29</td></tr><tr><td align="center">Brazil*</td><td align="center">12.5</td><td align="center">85</td><td align="center">12.5</td><td align="center">0</td><td align="center">27.5</td><td align="center">32.5</td><td align="center">27.5</td><td align="center">0</td><td align="center">55</td></tr></tbody></table></table-wrap>
<table-wrap-foot><p>a Showing the detection frequency of each amino acids variant on <italic>PvDBPII </italic>protein among 30 Thai isolates.</p><p>b Percentage of detection frequency calculated from value in panel "a".</p><p>▼ High prevalence (&gt;50%) of variant amino acids.</p><p>▽ Dimorphic mutation (changed into two amino acid-types)</p></table-wrap-foot> <table-wrap-foot><p>* Amino acid residues included in the haplotypes are as follows: 268, 306, 308, 333, 351, 367, 371, 375, 384, 385, 386, 390, 398, 404, 417, 419, 424, 433, 436, 437, 475, 486, 503, 507, 513.</p></table-wrap-foot> <table-wrap-foot><p>* The frequencies of the common variant amino acid data were referable according to Sousa et al. 2006.</p><p># The first letter represented the amino acid in that position in Sal I sequences, and the other letter represents the substituted amino acid.</p><p>COL: Colombia, PNG: Papua New Guinea</p></table-wrap-foot>
<graphic xlink:href="1475-2875-7-112-1"/> <graphic xlink:href="1475-2875-7-112-2"/> <graphic xlink:href="1475-2875-7-112-3"/>
[]
Malar J. 2008 Jun 26; 7:112