- Open Access
Genetic polymorphism and natural selection of circumsporozoite protein in Myanmar Plasmodium vivax
Malaria Journal volume 19, Article number: 303 (2020)
Circumsporozoite surface protein (CSP) of malaria parasites has been recognized as one of the leading vaccine candidates. Clinical trials of vaccines for vivax malaria incorporating Plasmodium vivax CSP (PvCSP) have demonstrated their effectiveness in preventing malaria, at least in part. However, genetic diversity of pvcsp in the natural population remains a major concern.
A total of 171 blood samples collected from patients infected with Plasmodium vivax in Myanmar were analysed in this study. The pvcsp was amplified by polymerase chain reaction, followed by cloning and sequencing. Polymorphic characteristics and natural selection of pvcsp population in Myanmar were analysed using DNASTAR, MEGA6 and DnaSP programs. The polymorphic pattern and natural selection of publicly accessible global pvcsp sequences were also comparatively analysed.
Myanmar pvcsp sequences were divided into two subtypes VK210 and VK247 comprising 143 and 28 sequences, respectively. The VK210 subtypes showed higher levels of genetic diversity and polymorphism than the VK247 subtypes. The N-terminal non-repeat region of pvcsp displayed limited genetic variations in the global population. Different patterns of octapeptide insertion (ANKKAEDA in VK210 and ANKKAGDA in VK247) and tetrapeptide repeat motif (GGNA) were identified in the C-terminal region of global pvcsp population. Meanwhile, the central repeat region (CRR) of Myanmar and global pvcsp, both in VK210 and VK247 variants, was highly polymorphic. The high level of genetic diversity in the CRR has been attributed to the different numbers, types and combinations of peptide repeat motifs (PRMs). Interestingly, 27 and 5 novel PRMs were found in Myanmar VK210 and VK247 variants, respectively.
Comparative analysis of the global pvcsp population suggests a complex genetic profile of pvcsp in the global population. These results widen understanding of the genetic make-up of pvcsp in the global P. vivax population and provide valuable information for the development of a vaccine based on PvCSP.
Malaria caused by Plasmodium species has threatened human health since ancient time . The mortality and morbidity of malaria has greatly decreased in recent years; however, an estimated 219 million cases and 435,000 deaths due to malaria were reported globally in 2017 . Among five human malaria parasites, Plasmodium vivax is the most prevalent species outside of Africa. More than 7.5 million cases of malaria were caused by the parasite, which accounted for 56% of total malaria cases in South East Asia . Plasmodium vivax has been linked to benign malarial infection due to its mild clinical manifestations compared with Plasmodium falciparum; however, concerns that P. vivax also causes serious clinical illnesses and even death have increased [3, 4]. The emergence and spread of drug-resistant strains also increases the burden of the parasite [5, 6]. Furthermore, the liver stages in the life cycle of P. vivax, known as hypnozoites, can remain dormant and survive for weeks to months in liver cells, and be reactivated . This unique biological characteristic is a challenge in the control and elimination of this parasite. Therefore, the development of an effective vaccine is imperative for effective control and elimination of the parasite.
Circumsporozoite protein (CSP) is the most abundantly expressed protein on the surface of sporozoites. CSP plays multiple and crucial roles in the development, migration and invasion of sporozoites into hepatocytes . Thus, this protein has been studied extensively as one of the most promising vaccine candidates. The gene encoding CSP consists of three distinct regions: a conserved non-repeat N-terminal region, a highly polymorphic central repeat region (CRR), and a conserved non-repeat C-terminal region (Additional file 1: Figure S1). The CRR of P. vivax CSP (PvCSP) consists of peptide repeat motifs (PRMs). Based on the composition of one of the three major PRMs in the CRR, pvcsp is classified into three allelic variants: VK210, VK247 and P. vivax-like. The two most abundant alleles include VK210 and VK247, which carry nonapeptide repeat motifs GDRA(D/A)GQPA and ANGAG(N/D)QPG, respectively . Meanwhile, the P. vivax-like variant contains APGANQ(E/G)GAA motifs in the CRR .
Although no licensed vaccine against vivax malaria is currently available, a few notable advances in the development of a potential vaccine have been reported. A vivax malaria subunit vaccine named VMP001 has been clinically tested. This vaccine uses a chimeric recombinant protein containing repeat sequences of two major alleles of pvcsp, including VK210 and VK247 . The phase 1 trial with VMP001 showed a significant delay in parasitaemia, even though vaccination did not induce fully sterilizing protective immunity . However, the geographical variation in the genetic make-up of pvcsp population may influence the effectiveness of PvCSP-based vaccine, which may be a barrier for the development of a universal vaccine. Therefore, a comprehensive analysis of genetic diversity and structure of the global pvcsp population is necessary. In this study, genetic polymorphism and natural selection of Myanmar pvcsp has been examined. A comparative analysis of global pvcsp was also conducted to obtain a deeper insight into the genetic nature of pvcsp in the global P. vivax population.
One hundred seventy-one blood samples from P. vivax-infected Myanmar patients were used in this study. The blood samples were collected from the patients in field surveys in towns and villages in Naung Cho, Pyin Oo Lwin, and Tha Beik Kyin in Upper Myanmar during 2013 to 2015 (Fig. 1). The age range of patients was 13–62 years, with median age of 28.4 years. Initial screening for malaria infection was done by microscopic examination of thin and thick blood smears. Finger-prick blood samples were taken from the P. vivax-infected patients before drug treatment and spotted on Whatman 3 MM filter paper (GE Healthcare, Maidstone, UK) for confirmation by species-specific polymerase chain reaction (PCR) targeting 18S ribosomal RNA (rRNA) gene [13, 14]. Informed consent was obtained from all of the patients before blood collection. The study protocol was reviewed and approved by either the Ethics committee of the Ministry of Health, Myanmar (97/Ethics 2015) and the Biomedical Research Ethics Review Board of Inha University School of Medicine, Republic of Korea (INHA 15-013).
Genomic DNA extraction and amplification of pvcsp
Genomic DNA was extracted from dried blood spots using QIAamp DNA Blood Kit (Qiagen, Hilden, Germany) following the manufacturer’s protocol. Amplification of pvcsp gene was performed with primers by using nested PCR. The primers for first round PCR were 5ʹ-ATGTAGATCTGTCCAAGGCCATAAA-3ʹ and 5ʹ-AATTGAATAATGCTAGGACTAACAATATG-3ʹ. The thermal cycling parameters for primary PCR were as follows: one cycle of initial denaturation at 95 °C for 5 min, 25 cycles of 94 °C for 1 min, annealing at 58 °C for 2 min and extension at 72 °C for 2 min, followed by a final extension at 72 °C for 5 min. The nested PCR was performed with primers, 5ʹ-GCAGAACCAAAAAATCCACGTGAAAATAAG-3ʹ and 5ʹ-CCAACGGTAGCTCTAACTTTATCTAGGTAT-3ʹ, and similar amplification condition except the annealing temperature was 68 °C. Ex Taq DNA polymerase (Takara, Otsu, Japan) with proof-reading activity was used in all PCR steps to minimize the nucleotide mismatching during the amplification. Each PCR product was analysed by electrophoresis on 2% agarose gel. The resulting PCR product was extracted from the gel and was cloned into T&A cloning vector (Real Biotech Corporation, Banqiao City, Taiwan). Each ligation mixture was transformed into Escherichia coli DH5α competent cells. To identify positive clones with appropriate insert, colony PCR with nested PCR primers was performed. The nucleotide sequences of the cloned pvcsp were analysed by the Sanger method with M13 forward and reverse primers. Plasmids from at least two independent clones from each transformation mixture were analysed to confirm the sequence accuracy. These nucleotide sequences analysed in this study were deposited at GenBank under the accession numbers MN821829–MN821999.
Analyses of genetic diversity and natural selection in Myanmar pvcsp
The nucleotide and deduced amino acid sequences of Myanmar pvcsp were analysed using Editseq and Seqman in the DNASTAR package (DNASTAR, Madison, WI, USA). Two major variants of pvcsp sequences, Salvador I (Sal I; GU339059) and Papua New Guinea (PNG; M69059), were used as reference sequences to analyse Myanmar pvcsp sequences. The values of number of segregating sites (S), haplotypes (H), haplotype diversity (Hd), nucleotide diversity (π), and average number of pair-wise nucleotide differences within a population (K) were calculated with DnaSP ver. 5.10.00 . The rate of synonymous (dS) and nonsynonymous (dN) substitutions were calculated and compared using the Z-test (P < 0.05) with MEGA6 program  using Nei and Gojobori’s method  with the Juke and Cantor (JC) correction of 1000 bootstrap replications. Based on the acquired vales, the dN–dS value was calculated. Positive value for dN–dS imply to positive natural selection whereas negative value correspond to negative or purifying natural selection . Tajima’s D test, Fu and Li’s D and F statistics analysis was performed using DnaSP ver. 5.10.00  to evaluate the neutral theory of natural selection. An excess of high-frequency variation is consistent with balancing selection and is indicated by a positive Tajima’s D and/or Fu and Li’s D [19, 20]. A negative value of Tajima’s D and/or Fu and Li’s D indicates an excess of rare alleles, which may result from a recent selective sweep or purifying selection.
Genetic diversity and natural of pvcsp among the global P. vivax population
The genetic diversity of pvcsp among the global P. vivax population was also analysed. The pvcsp sequences deposited in public database were used in this study; Cambodia (n = 41), India (n = 79), Iran (n = 50), South Korea (n = 39), Brazil (n = 41), Mexico (n = 19), Colombia (n = 25), Sudan (n = 30), and Vanuatu (n = 21) (Additional file 2: Table S1). Genetic polymorphism and test of neutrality were examined for each pvcsp population with DnaSP ver 5.10.00  and MEGA6  as described above. To analyse the polymorphic patterns of the N- and C-terminal non-repeat regions in global pvcsp, a logo plot was constructed for each population using the WebLogo program (https://weblogo.berkeley.edu/logo.cgi).
Amplification and sequence analysis of Myanmar pvcsp
A total of 171 Myanmar pvcsp genes were successfully amplified form the genomic DNA samples used in this study. The size of the amplified pvcsp genes ranged from 0.5 to 1.3 kb. Sequence analysis of the amplified pvcsp genes revealed only two variants of pvcsp, VK210 and VK247, in Myanmar pvcsp, but not the P. vivax-like variant. The VK210 variants were prevalent (n = 143, 83.6%) and the frequency of VK247 variants was 16.4% (n = 28). No mixed infection with the two different variants was detected.
Genetic diversity of the N-terminal non-repeat region of Myanmar pvcsp
The N-terminal non-repeat region of Myanmar pvcsp showed a limited range of genetic diversity. Alignment of the deduced amino acid sequences of Myanmar pvcsp revealed 7 distinct haplotypes of VK210 variants and 7 haplotypes of VK247 variants (Fig. 2). The N-terminal non-repeat region of Myanmar VK210 variants was highly conserved. Compared with Sal I (GU339059) sequence, only few amino acid substitutions was found in the latter portion of the N-terminal non-repeat region. Haplotype 1, which was identical to Sal I sequence, was predominant (n = 108, 75.5%). An alanine insertion at the end of the RI conserved motif (KLKQP) was identified in haplotypes 2, 3 and 4. The N86I was found in haplotype 3. Two dimorphic (K91R and P95S) and one trimorphic (K93E/N) amino acid changes were observed in the RI motif of haplotypes 4, 5, 6, and 7 (Fig. 2a). Changes in the amino acid sequences were also identified in VK247 variants of Myanmar pvcsp. Haplotype 1, which shared the same sequence with the PNG (M69059) sequence, was the most prevalent, accounting for 57.1% of 28 Myanmar VK247 variants. All amino acid changes identified in Myanmar VK247 variants were all dimorphic (E96G/A, G99R, and N100D). Furthermore, eight amino acids (97DGAGNQPG104) were not detected in haplotypes 6 and 7 (Fig. 2b).
Genetic polymorphisms of the N-terminal non-repeat region in global pvcsp
Overall genetic polymorphisms of the N-terminal non-repeat region in the global pvcsp population were analysed. A comparative analysis of the region revealed that the region is relatively well-conserved in the global pvcsp. Alanine insertion at the end of the RI in the VK210 variants was the major variation identified in the global pvcsp, but the prevalence of the insertion differed by geographically (Fig. 3a). The pvcsp from Sudan showed 100% alanine insertion followed by pvcsp from Cambodia (67.7%), Vanuatu (47.6%), Myanmar (21.7%), Brazil (9.8%), and India (7.6%). No alanine insertion was found in the VK210 sequences identified in Iran, South Korea, and Mexico. Indian VK210 variants showed the highest genetic diversity with amino acid substitutions at 10 positions including E83K/G, K85T/Q, N86K/Y/I, P87A, R88G, N90I, L92V, K93N, Q94H/P and P95G, even though their frequencies were low, less than 15%. Meanwhile, A82T, N86I, N86S, K91R, P95S, and K93E/N showed uneven geographic distribution in the global VK210 variants with very low frequencies. The N-terminal non-repeat region of global VK247 variants were also well-conserved, although low frequencies of uneven amino acid changes were identified (Fig. 3b). The most remarkable variation involved N100D, which was observed in VK247 variants from Colombia (100%), Mexico (100%), Iran (27.3%), and Myanmar (17.9%). The amino acid changes such as E96G/A, A99E, G100R, and N101D showed uneven geographic distribution with low frequencies.
Polymorphic pattern of the CRR in Myanmar pvcsp
The CRR of Myanmar pvcsp showed extreme diversity in both VK210 and VK247 variants. A total of 118 and 23 haplotypes were identified in VK210 and VK247 variants, respectively. As expected, the greatest diversity of Myanmar pvcsp CRR was mainly attributed to differences in numbers, types, and arrangements of PRMs in each haplotype. A total of 47 different types of PRMs have been identified in the CRR of Myanmar VK210 variants (Fig. 4). Among these PRMs, two major types GDRADGQPA and GDRAAGQPA were the dominant ones for VK210 variants. Twenty-seven novel PRMs including GDRVAGQPA, GDRAHGQPA, GDRADGKPA, GDRADRQPA, GDGAGGQAA, EDRAAGQPA, GDKAAGQPA, GDRAAGLPA, GDRADVQPA, GDRADGQPV, GDRADGRPA, GDRADGLPA, GDRADGQPT, GDRAARQPA, GDRAAGRPA, GDGAGGQPA, GDRAAGQSA, SDRADGRPA, GDRAAGQPT, GDRAYGQPA, SDRAAGQPA, RDRADGQPA, GDRASGQPA, GGRADGRPA, GDRADQQPA, GDRADGPPA and GNGADGQPA, which were not previously reported, were found in Myanmar VK210 variants. Interestingly, 109 haplotypes out of 118 VK210 variants were terminated with GNGAGGQAA motif. The number of PRMs consisting CRR of Myanmar VK210 haplotypes varied from 1 to 29. Sequences with 18 PRMs were the most prevalent accounting for 18.9% of 143 Myanmar VK210 sequences (Fig. 5).
Compared with global VK210 variants, the CRR of Myanmar VK210 variants showed a high level of length polymorphisms. The VK210 variants from other countries analysed in this study showed only a few (2 to 5) variations in length of polymorphisms in CRR; however, the CRR of Myanmar VK210 variants had 24 different length polymorphisms consistent with various types and different compositions of PRMs. The overall genetic diversity in CRR of Myanmar VK247 variants was much lower than that of Myanmar VK210 variants. A total of 23 VK247 haplotypes, each CRR comprising different numbers and combinations of 8 types of PRMs, were identified (Fig. 6). Five of 8 PRMs, ANGAGNQSG, AYGAGNQPG, VNGAGNQPG, ANGVGNQPG and AYGAGNQPG, identified in Myanmar VK247 variants were novel ones that have not reported previously. The CRR of Myanmar VK247 variants carried different numbers of PRMs ranging from 1 to 22, and the CRR with 2 PRMs was the most prevalent (14.3%) (Fig. 7). Similar to VK210 variants, the Myanmar VK247 variants also displayed higher levels of diversity than the VK247 variants identified from other geographical origins. Sixteen different size polymorphisms of CRR were identified in Myanmar VK247 variants, whereas less size variations (1 to 5) were found in CRR of all other countries analysed in this study.
Genetic diversity in the C-terminal non-repeat region of Myanmar and global pvcsp
Sequence analysis of the C-terminal non-repeat region of Myanmar VK210 variants revealed 27 distinct haplotypes (Fig. 8a). These sequence diversities were attributed to differences in the arrangement of ANKKAEDA octapeptide insertion and GGNA tetrapeptide repeat motifs with different amino acid substitutions throughout the region. The sequence of haplotype 22 was identical with Sal I (GU339059), and accounted for 9.8% of all the VK210 sequences. The octapeptide insertions were observed in haplotypes 1 to 16. All the inserted sequences in these 16 haplotypes were ANKKAEDA except for haplotypes 8 and 13, which contained ANKKAENA and ANKEAENA, respectively. The C-terminal non-repeat region of Myanmar VK247 variants showed a lower level of genetic diversity than that of Myanmar VK210 variants (Fig. 8b). A total of 10 haplotypes were identified, and haplotype 1, which carried a sequence identical to that of the reference sequence of PNG (M69059), was the most prevalent haplotype with a frequency of 46.4%. The most noteworthy polymorphic characteristics identified in the C-terminal non-repeat region of Myanmar VK247 variants were the deletions of GGQAAGGNAANKKAGDAG in haplotype 7 and ANKKAGDAG in haplotypes 8, 9 and 10. Analysis of sequence polymorphisms in the C-terminal non-repeat region of the global pvcsp suggested a high level of genetic diversity. Among VK210 variants, the frequency of ANKKAEDA insertion differed by country (Fig. 9a). All sequences from Iran and South Korea VK210 variants contained an octapeptide insertion, but the frequencies in VK210 variants from Sudan, India, Mexico and Cambodia were 96.7, 89.9, 63.6 and 9.7%, respectively. Interestingly, no insertion of the octapeptide was identified in VK210 variants from Brazil and Vanuatu. The numbers of GGNA motifs in the C-terminal non-repeat region of global VK210 variants also differed by country (Fig. 9b). The number of repeated GGNA motifs found in global VK210 variants ranged from 0 to 6. Similar to VK210 variants, the presence of ANKKAGDAG octapeptide insertion and the number of GGNA motifs also differed in the C-terminal non-repeat region of global VK247 variants. The frequency of ANKKAGDAG octapeptide insertion was high in VK247 variants isolated from Cambodia, Iran, Mexico, and Colombia, but was low in Myanmar VK247 variants (Fig. 9c). The number of GGNA motifs found in the C-terminal non-repeat region of global VK247 variants ranged from 0 to 3 (Fig. 9d). The number of GGNA motifs in global VK247 also differed by country.
Nucleotide diversity and natural selection in the N- and C-terminal non-repeat regions of Myanmar pvcsp
Since the CRR of Myanmar pvcsp sequences showed a high degree of length polymorphisms, the nucleotide diversity and genetic differentiation of the N- and C-terminal non-repeat regions were analysed separately by omitting the CRR. In the N-terminal non-repeat region of Myanmar VK210 variants, the average number of nucleotide differences (K), overall haplotype diversity (Hd), and nucleotide diversity (π) were 0.098, 0.096 ± 0.034, and 0.0024 ± 0.0009, respectively (Table 1). The estimated dN–dS value in the N-terminal non-repeat regions was 0.0008. These results suggested that the N-terminal non-repeat region of Myanmar VK210 variants was under positive natural selection. Meanwhile, the average number of nucleotide differences (K), the overall haplotype diversity (Hd), and the nucleotide diversity (π) for the C-terminal non-repeat region of Myanmar VK210 variants were 0.209, 0.186 ± 0.044, and 0.0035 ± 0.0009, respectively (Table 1). The dN–dS value was –0.0035. These findings indicated that the C-terminal region was influenced by negative natural selection. The Tajima’s D test was also performed to further elucidate the effect of natural selection on the N- and C-terminal non-repeat regions in Myanmar VK210 variants. Tajima’s D values for the N- and C-terminal non-repeat regions were –1.9359 (P < 0.05) and –2.2251 (P < 0.01), respectively (Table 1). The Fu and Li’s D and F values of these regions also showed negative values. In the Myanmar VK247 variants, the average number of nucleotide differences (K), overall haplotype diversity (Hd) and nucleotide diversity (π) of the N-terminal non-repeat region were 0.405, 0.405 ± 0.094 and 0.0090 ± 0.0021, respectively (Table 1). These values for the C-terminal non-repeat region were 0.495, 0.331 ± 0.114 and 0.0067 ± 0.0027, respectively. The dN–dS values of the N- and C-terminal non-repeat regions of Myanmar VK247 variants were 0.0117 and –0.0150, respectively (Table 1). The Tajima’s D values were –0.4445 (P > 0.1) and –1.9719 (P < 0.05) for the N- and C-terminal non-repeat regions, respectively (Table 1). The Fu and Li’s D and F values for both regions were all negative.
Nucleotide diversity and natural selection in the N- and C-terminal non-repeat regions of global VK210 variants
To further examine the nucleotide diversity and natural selection in the global pvcsp population, the nucleotide diversity of the N- and C-terminal non-repeat regions of global pvcsp was analysed. For the N-terminal non-repeat region of VK210 variants, nucleotide diversity and pattern of natural selection were differed by country (Table 2). VK210 variants from India showed the highest nucleotide diversity; the values of K, Hd, and π were 1.972, 0.661 ± 0.063, and 0.0470 ± 0.0074, respectively. The N-terminal non-repeat region of Cambodia VK210 variants also showed relatively high nucleotide diversity comparable to Myanmar VK210 variants. Substantial nucleotide diversity was found in the VK210 variants from South Korea and Brazil. However, VK210 variants from Iran, Mexico, Sudan, and Vanuatu were genetically well-conserved. The N-terminal non-repeat region of VK210 variants showed different patterns of natural selection by country. The dN–dS value was positive for VK210 variants from Myanmar, Cambodia, and South Korea, which indicated positive natural selection may occur in the region. Meanwhile, the values for VK210 variants from India and Brazil were negative, suggesting negative selection. The values of Tajima’s D for all VK210 variants derived from Myanmar, Cambodia, India, South Korea, and Brazil were negative, indicating that they were under purifying selection. The C-terminal non-repeat region of global VK210 variants also revealed nucleotide diversity and pattern of natural selection (Table 2). Indian VK210 variants showed the greatest nucleotide diversity with K, Hd, and π values of were 0.276, 0.123 ± 0.051, and 0.0043 ± 0.0020, respectively. Meanwhile, no nucleotide diversity was detected in VK210 variants from Cambodia, Iran, Brazil, Mexico, and Vanuatu. Similar to the N-terminal non-repeat region, the C-terminal non-repeat region of global VK210 was also under the effect of purifying selection based on the negative values of Tajima’s D. The values of Fu and Li’s D and Fu and Li’s F were also negative for the C-terminal region of VK210 variants from Myanmar, India, South Korea, and Sudan.
Nucleotide diversity and natural selection in the N-terminal and C-terminal non-repeat regions of global VK247 variants
The nucleotide diversity and natural selection in the global VK247 variants were analysed (Table 3). Analysis of the N-terminal non-repeat region of VK247 variants revealed the greatest nucleotide diversity in VK247 variants from Iran, with the values of K, Hd, and π of 1.309, 0.436 ± 0.133, and 0.0190 ± 0.0058, respectively. VK247 variants from Cambodia also revealed nucleotide diversity in the region. However, no nucleotide diversity was found in VK247 variants from Mexico and Colombia. The N-terminal non-repeat region of VK247 variants from Iran showed negative dN–dS (–0.0543) and positive Tajima’s D (0.9518), Fu and Li’s D (1.1271), and Fu and Li’s F (1.2185). However, Myanmar and Cambodia VK247 variants revealed positive values for dN–dS and negative values for Tajima’s D, Fu and Li’s D, and Fu and Li’s F. The C-terminal non-repeat region of VK267 variants from Cambodia and Colombia also showed nucleotide diversity, which was lower than that of Myanmar VK247 variants. Similar to Myanmar VK247 variants, the C-terminal non-repeat region of Cambodia VK247 variants showed a negative values of dN–dS (–0.0172), Tajima’s D (–1.4009), Fu and Li’s D (–1.5866), and Fu and Li’s F (–1.7190). However, the dN–dS of Colombia variants was estimated to be positive (0.0017), although the values of Tajima’s D, Fu and Li’s D, and Fu and Li’s F were negative.
PvCSP is one of the leading candidates for vivax malaria vaccine. Several recent PvCSP-derived vaccines were designed as multivalent formulations or chimeric molecules in an attempt to induce protective immunity against P. vivax [11, 21, 22]. However, the impact of natural genetic variations in the global pvcsp population on vaccine efficacy remains unclear. In this study, genetic polymorphism and natural selection in Myanmar pvcsp and global pvcsp populations were comprehensively analysed. Among the three pvcsp allelic variants, only two allelic variants, VK210 and VK247, were identified in P. vivax isolates from Myanmar analysed in this study. The VK210 was the dominant allele occurring in 83.6% of Myanmar pvcsp population, which is consistent with previous studies reporting that VK210 is the most common allelic variant in the pvcsp populations from Iran, India, China, Brazil, Thailand, Bangladesh, and Azerbaijan [23,24,25,26,27,28,29,30,31]. Meanwhile, VK247 was more prevalent in certain regions of Colombia . It has been suggested that the distribution of Anopheles mosquito species may affect the prevalence of VK210 and VK247 variants in endemic areas [32, 33]. In Myanmar, diverse species of Anopheles are distributed throughout the country, and at least 10 species including Anopheles dirus, Anopheles minimus, and Anopheles aconitus may transmit malaria . Differences in infectivity of VK210 and VK247 variants in certain Anopheles species is not clear; however, the differences in adaptability of each allelic variant to mosquito vectors affects the prevalence of VK210 in Myanmar. The predominance of VK210 may also be associated with its genetic diversity. VK247 showed a lower level of genetic diversity than VK210 in the Myanmar pvcsp population, suggesting different polymorphisms may increase the adaptability of VK210 to different mosquito vectors.
The N-terminal non-repeat region of Myanmar pvcsp was relatively well-conserved. Only a few amino acid changes were identified in the latter portion of the N-terminal non-repeat region of Myanmar pvcsp. The RI motif (KLKQP), a cell adhesive motif exposed by proteolytic cleavage after the interaction between sporozoite and hepatocyte , was well conserved in both Myanmar VK210 and VK247 variants. Global VK210 variants also displayed limited polymorphism in the N-terminal non-repeat region. An alanine insertion at the end of the RI was the major variation identified in global pvcsp, although the prevalence of the insertion differed by geographically. Amino acid changes were also identified in global VK210 variants, but they showed uneven geographic distribution with low frequencies. The N-terminal non-repeat region of global VK247 variants was also relatively well-conserved, although uneven and minor amino acid changes were also identified. It has been shown that the invasion of malaria parasite is inhibited by antibodies directed against RI motif . The highly conserved genetic characteristic of the N-terminal region of global pvcsp suggest that this region including RI motif may represent an attractive candidate for formulation of a PvCSP-based vaccine.
As expected, high level of genetic diversity was identified in the CRR of Myanmar pvcsp. Both VK210 and VK247 variants showed great genetic diversity in CRR resulting in 118 and 23 different haplotypes for VK210 and VK247, respectively. Insertions and deletions of nonapeptide sequences probably resulted from either sexual recombination during meiosis or intrahelical strand slippage during mitotic DNA replication , which may generate novel haplotypes in Myanmar pvcsp. Interestingly, 47 different types of PRMs including 27 novel types were identified in the CRR of Myanmar VK210 variants. In the case of VK247 variants, a total of 8 distinct PRMs were identified, in which 5 were novel ones. These findings also suggest point mutations as one of the major factors increasing the genetic complexity of Myanmar pvcsp. Similar patterns of high genetic polymorphisms were also identified in the CRR of global pvcsp population; however, the overall diversity due to different types and repeats of PRMs was greater in the pvcsp from Myanmar compared with other countries. Especially, the repeated numbers of PRMs varied widely in both Myanmar VK210 and VK247 variants compared with those from other countries.
The C-terminal regions of Myanmar and global pvcsp also showed genetic polymorphisms. Differences in ANKKAEDA octapeptide insertion and GGNA tetrapeptide repeat motifs and amino acid substitutions throughout the region contributed to genetic diversity in Myanmar and global VK210 variants. Despite the unclear biological functions of ANKKAEDA octapeptide and GGNA tetrapeptide in VK210 variants, they may also be generated via intrahelical recombination, and parasites carrying the ANKKAEDA octapeptide show a high degree of delayed infections . Similar to VK210 variants, ANKKAGDAG octapeptide and GGNA motifs were important factors contributing to genetic diversity in Myanmar and global VK247 variants. Although the frequency of octapeptide insertion and number of GGNA motifs in global pvcsp, both in VK210 and VK247 variants, were differed by country, no clear geographic clustering was identified. Considering the limited genetic information of the C-terminal region in global pvcsp in current, a further study employing larger numbers of pvcsp sequences in diverse geographical origins is necessary to understand the genetic nature and evolutionary aspects of the C-terminal region in global pvcsp population.
Analysis of the two non-repeat regions in global pvcsp suggests that these regions are likely to be under natural selection, which may maintain or generate genetic diversity in the global pvcsp population. The dN–dS values for Myanmar and Cambodia VK210 and VK247 variants were positive, implying that balancing selection might act in these regions. However, the values were 0 or negative for both VK210 and VK247 variants from other countries except the N-terminal region of South Korea VK210 and the C-terminal region of Colombia VK247. Tajima’s D and Fu and Li’s D and F values also revealed complex patterns of natural selection that were unique to pvcsp from each country. While pvcsp populations from a few countries showed no genetic diversity involving the N- and C-terminal regions, a few pvcsp populations showed negative Tajima’s D values, suggesting purifying natural selection. The overall trend indicated that the global pvcsp population decreased the genetic diversity in the two terminal non-repeat regions. Due to size polymorphisms induced by different numbers and arrangements of PRMs in CRR, a direct analysis of natural selection in CRR was not possible in this study. However, evidence supporting natural selection in the CRR have been reported [36,37,38]. These findings collectively suggest that complex natural selection phenomena may act on the global pvcsp.
Despite the remarkable reduction in malaria transmission in Myanmar during the last decades [2, 38], the high levels of genetic diversity among malaria parasites appears to persist in the country. As demonstrated in this study, the overall genetic polymorphisms of Myanmar pvcsp were greater than those reported from other countries analysed. Similar patterns of high genetic diversity in several polymorphic genetic markers including apical membrane antigen-1, merozoite surface protein-1 and -2, and CSP of Myanmar P. falciparum population have also been reported [39,40,41,42]. Despite the lack of clear insight into the high genetic diversity underlying the population decline, asymptomatic carriers may act as fundamental reservoirs contributing to malaria transmission, which provides adequate population size to maintain or generate genetic diversity of malaria parasites in Myanmar . This study also had a limitation. The Myanmar P. vivax isolates analysed in this study were collected in restricted areas of Myanmar and, therefore, not fully representative of the nation-wide genetic diversity and population structure of Myanmar pvcsp. A further comprehensive study involving a larger number of P. vivax isolates collected from different regions in Myanmar is necessary.
PvCSP is a leading candidate for vivax vaccine; however, the impact of natural genetic variations in the global pvcsp population on the efficacy of PvCSP-based vaccine remains unclear. This study highlights an additional level of complexity in the application of this antigen for the development of a vivax vaccine. Global pvcsp population revealed polymorphic characters, especially in the CRR. The size polymorphisms and the appearance of novel PRMs as well as new CRR arrays primarily contribute to the genetic diversity of global pvcsp. Natural selection may also be a major force affecting the genetic diversity of global pvcsp. The results of this study provide not only an insight into the genetic nature of global pvcsp but also valuable information for the development of a universal vaccine based on PvCSP. Continuous monitoring of genetic diversity of global pvcsp population is necessary to elucidate the polymorphic nature and evolutionary aspect of pvcsp in global P. vivax population.
Availability of data and materials
The data supporting the conclusions of this article are provided within the article and its additional files. The original datasets analysed in this current study are available from the corresponding author upon request. The nucleotide sequences reported in this study have been deposited in the GenBank database under the Accession Numbers MN821829–MN821999.
Central repeat region
Circumsporozoite surface protein
Rate of non-synonymous mutations
Rate of synonymous mutations
Number of haplotypes
- K :
Average number of nucleotide differences
Polymerase chain reaction
Peptide repeat motif
Plasmodium vivax circumsporozoite surface protein
18S ribosomal RNA
Number of segregating sites
Observed average pairwise nucleotide diversity
Cox FE. History of human parasitology. Clin Microbiol Rev. 2002;15:595–612.
WHO. World Malaria report 2018. Geneva, World Health Organization, 2018.
Rogerson SJ, Carter R. Severe vivax malaria: newly recognised or rediscovered. PLoS Med. 2008;5:e136.
Anstey NM, Russell B, Yeo TW, Price RN. The pathophysiology of vivax malaria. Trends Parasitol. 2009;25:220–7.
Baird JK. Chloroquine resistance in Plasmodium vivax. Antimicrob Agents Chemother. 2004;48:4075–83.
Price RN, Tjitra E, Guerra CA, Yeung S, White NJ, Anstey NM. Vivax malaria: neglected and not benign. Am J Trop Med Hyg. 2007;77:79–87.
Miller LH, Baruch DI, Marsh K, Doumbo OK. The pathogenic basis of malaria. Nature. 2002;415:673–9.
Coppi A, Natarajan R, Pradel G, Bennett BL, James ER, Roggero MA, et al. The malaria circumsporozoite protein has two functional domains, each with distinct roles as sporozoites journey from mosquito to mammalian host. J Exp Med. 2011;208:341–56.
Rosenberg R, Wirtz RA, Lanar DE, Sattabongkot J, Hall T, Waters AP, et al. Circumsporozoite protein heterogeneity in the human malaria parasite Plasmodium vivax. Science. 1989;245:973–6.
Qari SH, Shi YP, Goldman IF, Udhayakumar V, Alpers MP, Collins WE, et al. Identification of Plasmodium vivax-like human malaria parasite. Lancet. 1993;341:780–3.
Yadava A, Sattabongkot J, Washington MA, Ware LA, Majam V, Zheng H, et al. A novel chimeric Plasmodium vivax circumsporozoite protein induces biologically functional antibodies that recognize both VK210 and VK247 sporozoites. Infect Immun. 2007;75:1177–85.
Bennett JW, Yadava A, Tosh D, Sattabongkot J, Komisar J, Ware LA, et al. Phase 1/2a trial of Plasmodium vivax malaria vaccine candidate VMP001/AS01B in malaria-naive adults: safety, immunogenicity, and efficacy. PLoS Negl Trop Dis. 2016;10:e0004423.
Snounou G, Viriyakosol S, Jarra W, Thaithong S, Brown KN. Identification of the four human malaria parasite species in field samples by the polymerase chain reaction and detection of a high prevalence of mixed infections. Mol Biochem Parasitol. 1993;58:283–92.
Kang JM, Cho PY, Moe M, Lee J, Jun H, Lee HW, et al. Comparison of the diagnostic performance of microscopic examination with nested polymerase chain reaction for optimum malaria diagnosis in Upper Myanmar. Malar J. 2017;16:119.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Nei M, Gojobori T. Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol. 1986;3:418–26.
Escalante AA, Cornejo OE, Rojas A, Udhayakumar V, Lal AA. Assessing the effect of natural selection in malaria parasites. Trends Parasitol. 2004;20:388–95.
Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.
Fu YX, Li WH. Statistical tests of neutrality of mutations. Genetics. 1993;133:693–709.
Gimenez AM, Lima LC, Francoso KS, Denapoli PMA, Panatieri R, Bargieri DY, et al. Vaccine containing the three allelic variants of the Plasmodium vivax circumsporozoite antigen induces protection in mice after challenge with a transgenic rodent malaria parasite. Front Immunol. 2017;8:1275.
Teixeira LH, Tararam CA, Lasaro MO, Camacho AG, Ersching J, Leal MT, et al. Immunogenicity of a prime-boost vaccine containing the circumsporozoite proteins of Plasmodium vivax in rodents. Infect Immun. 2014;82:793–807.
Shabani SH, Zakeri S, Mehrizi AA, Mortazavi Y, Djadid ND. Population genetics structure of Plasmodium vivax circumsporozoite protein during the elimination process in low and unstable malaria transmission areas, southeast of Iran. Acta Trop. 2016;160:23–34.
Kim JR, Imwong M, Nandy A, Chotivanich K, Nontprasert A, Tonomsing N, et al. Genetic diversity of Plasmodium vivax in Kolkata, India. Malar J. 2006;5:71.
Machado RL, Povoa MM. Distribution of Plasmodium vivax variants (VK210, VK247 and P. vivax-like) in three endemic areas of the Amazon region of Brazil and their correlation with chloroquine treatment. Trans R Soc Trop Med Hyg. 2000;94:377–81.
Cui L, Mascorro CN, Fan Q, Rzomp KA, Khuntirat B, Zhou G, et al. Genetic diversity and multiple infections of Plasmodium vivax malaria in Western Thailand. Am J Trop Med Hyg. 2003;68:613–9.
Leclerc MC, Menegon M, Cligny A, Noyer JL, Mammadov S, Aliyev N, et al. Genetic diversity of Plasmodium vivax isolates from Azerbaijan. Malar J. 2004;3:40.
Pratt-Riccio LR, Baptista BO, Torres VR, Bianco-Junior C, Perce-Da-Silva DS, Riccio EKP, et al. Chloroquine and mefloquine resistance profiles are not related to the circumsporozoite protein (CSP) VK210 subtypes in field isolates of Plasmodium vivax from Manaus. Brazilian Amazon. Mem Inst Oswaldo Cruz. 2019;114:e190054.
Kibria MG, Elahi R, Mohon AN, Khan WA, Haque R, Alam MS. Genetic diversity of Plasmodium vivax in clinical isolates from Bangladesh. Malar J. 2015;14:267.
Li YC, Wang GZ, Meng F, Zeng W, He CH, Hu XM, et al. Genetic diversity of Plasmodium vivax population before elimination of malaria in Hainan Province China. Malar J. 2015;14:78.
Kaur H, Sehgal R, Kumar A, Sehgal A, Bharti PK, Bansal D, et al. Exploration of genetic diversity of Plasmodium vivax circumsporozoite protein (Pvcsp) and Plasmodium vivax sexual stage antigen (Pvs25) among North Indian isolates. Malar J. 2019;18:308.
Gonzalez JM, Hurtado S, Arevalo-Herrera M, Herrera S. Variants of the Plasmodium vivax circumsporozoite protein (VK210 and VK247) in Colombian isolates. Mem Inst Oswaldo Cruz. 2001;96:709–12.
Rodriguez MH, Gonzalez-Ceron L, Hernandez JE, Nettel JA, Villarreal C, Kain KC, et al. Different prevalences of Plasmodium vivax phenotypes VK210 and VK247 associated with the distribution of Anopheles albimanus and Anopheles pseudopunctipennis in Mexico. Am J Trop Med Hyg. 2000;62:122–7.
Ool TT, Storch V, Becker N. Review of the anopheline mosquitoes of Myanmar. J Vector Ecol. 2004;29:21–40.
Ying P, Shakibaei M, Patankar MS, Clavijo P, Beavis RC, Clark GF, et al. The malaria circumsporozoite protein: interaction of the conserved regions I and II-plus with heparin-like oligosaccharides in heparan sulfate. Exp Parasitol. 1997;85:168–82.
Hughes AL. The evolution of amino acid repeat arrays in Plasmodium and other organisms. J Mol Evol. 2004;59:528–35.
Patil A, Orjuela-Sanchez P, da Silva-Nunes M, Ferreira MU. Evolutionary dynamics of the immunodominant repeats of the Plasmodium vivax malaria-vaccine candidate circumsporozoite protein (CSP). Infect Genet Evol. 2010;10:298–303.
Parobek CM, Bailey JA, Hathaway NJ, Socheat D, Rogers WO, Juliano JJ. Differing patterns of selection and geospatial genetic diversity within two leading Plasmodium vivax candidate vaccine antigens. PLoS Negl Trop Dis. 2014;8:e2796.
WHO. Malaria in the Greater Mekong Subregion: regional and country profiles. Geneva, World Health Organization, 2010.
Kang JM, Lee J, Moe M, Jun H, Le HG, Kim TI, et al. Population genetic structure and natural selection of Plasmodium falciparum apical membrane antigen-1 in Myanmar isolates. Malar J. 2018;17:71.
Le HG, Kang JM, Moe M, Jun H, Thai TL, Lee J, et al. Genetic polymorphism and natural selection of circumsporozoite surface protein in Plasmodium falciparum field isolates from Myanmar. Malar J. 2018;17:361.
Le HG, Kang JM, Jun H, Lee J, Thai TL, Myint MK, et al. Changing pattern of the genetic diversities of Plasmodium falciparum merozoite surface protein-1 and merozoite surface protein-2 in Myanmar isolates. Malar J. 2019;18:241.
We thank the staffs in the Department of Medical Research Pyin Oo Lwin Branch and the health professionals in Naung Cho and Pyin Oo Lwin townships for their contribution and technical support during blood collection.
This research was supported by the National Research Foundation of Korea (NRF) Grants funded by the Korean Government (NRF-2019K1A3A9A01000005).
Ethics approval and consent to participate
This study was approved by the Ethics Review Committee, Department of Medical Research, Myanmar (97/Ethics 2015) and by the Ethical Review Committee of Inha University School of Medicine, Korea (INHA 15-013). Informed written consent and permission were obtained from each individual.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Schematic structure of pvcsp. The gene is separated into three regions; an N-terminal non-repeat region, a central repeat region (CRR), and a C-terminal non-repeat region. The CRR consists of two major repeat peptide motifs (PRMs), termed VK210 and VK247.
Global pvcsp sequences analysed in this study.
About this article
Cite this article
Võ, T.C., Lê, H.G., Kang, JM. et al. Genetic polymorphism and natural selection of circumsporozoite protein in Myanmar Plasmodium vivax. Malar J 19, 303 (2020). https://doi.org/10.1186/s12936-020-03366-7
- Plasmodium vivax
- Circumsporozoite protein
- Genetic polymorphism
- Natural selection