Exploration of genetic diversity of Plasmodium vivax circumsporozoite protein (Pvcsp) and Plasmodium vivax sexual stage antigen (Pvs25) among North Indian isolates

Exploring genetic polymorphisms among Plasmodium vivax isolates from the Thai-Myanmar borders using circumsporozoite protein (pvcsp) and ookinete surface protein (pvs25) encoding genes

Plasmodium vivax malaria cases remain high along the Thai-Myanmar and Thai-Cambodia borders. Plasmodium vivax circumsporozoite protein (pvcsp) and Plasmodium vivax ookinete surface protein (pvs25) genes are promising molecular markers of the genetic diversity of P. vivax. This study investigated the genetic diversity of pvcsp and pvs25 in P. vivax isolates collected from the Thai-Myanmar border. The DNA samples were amplified, and the genotypes were analyzed by PCR-RFLP and DNA sequencing. Pvcsp genotypes, VK210, VK247, and mixed types, were found in 203 (91.9%), 15 (6.8%), and 3 (1.3%) of the isolates, respectively. Twenty-four allelic variants were observed, of which a high prevalence of VK210E and VK247E were reported. Two pvcsp variants, VK210C and VK210M showed significantly higher parasite density (46,234 (1154-144,000) vs. 25,606 (1373-68,878), respectively). The genetic diversity of pvcsp along the Thai-Myanmar border during 2002-2015 showed dynamic changes with both positive and negative selection. The frequency and distribution of pvcsp pattern might be changed over time and might be other factors contributing to gene selection. Three amino acid substitutions of pvs25, i.e., E97Q, I130T, and Q131K, were investigated with frequencies of 10 (4.5%), 221 (100%), and 204 (92.3%) isolates, respectively. There was no association between parasite density and pvs25 polymorphisms. The frequency of pvs25 polymorphism was similar to that previously reported, with the absence of random mutation. In conclusion, the genetic variation of pvcsp was changed over times whereas the genetic diversity of pvs25 was limited; these variations would be helpful for further vaccine development against P. vivax malaria.

Genetic Diversity of Plasmodium vivax Field Isolates from the Thai–Myanmar Border during the Period of 2006–2016

High levels of genetic variants of Plasmodium vivax have previously been reported in Thailand. Circumsporozoite surface protein (CSP), merozoite surface protein (MSP), and microsatellite markers were used to determine the genetic polymorphisms of P. vivax. This study aimed to investigate the molecular epidemiology of P. vivax populations at the Thai–Myanmar border by genotyping the PvCSP, PvMSP-3α, and PvMSP-3β genes. Four hundred and forty P. vivax clinical isolates were collected from the Mae Sot and Sai Yok districts from 2006–2007 and 2014–2016. Polymerase chain reaction with restriction fragment length polymorphism (RFLP) was used to investigate the genetic polymorphisms of the target genes. Based on PCR band size variations, 14 different PvCSP alleles were identified: eight for VK210 and six for VK247. The VK210 genotype was the dominant variant during both sample collection periods. Based on PCR genotyping, three distinct types (A, B, and C) for both PvMSP-3α and PvMSP-3β were observed. Following RFLP, 28 and 14 allelic variants of PvMSP-3α and 36 and 20 allelic variants of PvMSP-3β with varying frequencies were identified during the first and second periods, respectively. High genetic variants of PvMSP-3 and PvCSP were found in the study area. PvMSP-3β exhibited a higher level of genetic diversity and multiple-genotype infection versus PvMSP-3α.

Identification of immunodominant epitopes in allelic variants VK210 and VK247 of Plasmodium Vivax Circumsporozoite immunogen

Plasmodium vivax-induced malaria is among the leading causes of morbidity and mortality in sub-tropical and tropical regions and infect 2.85 billion people globally. The continual rise and propagation of resistance against anti-malarial drugs is a prerequisite to develop a potent vaccine candidate for Plasmodium vivax (P. vivax). Circumsporozoite protein (CSP) is an important immunogen of malaria parasite that has the conserved CSP structure as an immune dominant B-cell epitope. In current study, we focused on designing multi-epitope vaccines (MEVs) using various immunoinformatics tools against Pakistani based allelic variants VK210 and VK247 of P. vivax CSP (PvCSP) gene. Antigenicity, allergic potential and physicochemical parameters of both PvCSP variants were assessed for the designed MEVs and they were within acceptable range suitable for post experimental investigations. The three-dimensional structures of both MEVs have been predicted ab initio, optimized, and validated by using different online servers. The both MEVs candidates were stable and free from aggregation-prone regions. The stability of both MEVs had been improved by a disulfide engineering approach. To estimate the binding energy and stability of the MEVs, molecular docking simulation and binding free energy calculations with TLR-4 immune receptor have been conducted. The docking score of PvCSP210 and PvCSP247 for TLR-4 was -6.34 kJ/mol and - 2.3 kJ/mol, respectively. For PvCSP210-TLR4 system, mean RMSD was 4.96 Å while PvCSP247-TLR4 system, average RMSD was 4.49 Å. The binding free energy of PvCSP210-TLR4 complex and PvCSP247-TLR4 complex was -50.49/-117.15 kcal/mol (MMGBSA/MMPSA) and -52.94/-96.26 kcal/mol (MMGBSA/MMPSA), respectively. The expression of both MEVs produced in Escherichia coli K12 expression system by in silico cloning was significant. Immune simulation revealed that the proposed MEVs induce strong humoral and cellular immunological responses, in addition to significant production of interleukins and cytokines. In conclusions, we believed that the MEVs proposed in current research, using combine approach of immunoinformatics, structural biology and biophysical approaches, could induce protective and effective immune responses against P. vivax and the experimental validation of our findings could contribute to the development of potential malaria vaccine.

Genetic characterization of Plasmodium vivax isolates from Pakistan using circumsporozoite protein (pvcsp) and merozoite surface protein-1 (pvmsp-1) genes as genetic markers

Background

Plasmodium vivax contributes to over 70% malaria burden in Pakistan, but limited data exists on various aspects including genetic diversity of the parasite as compared to other parts of the world. Since the information about the genetic diversity of P. vivax assists to understand the population dynamics of the parasite, the current study was designed to understand population divergence of P. vivax in Pakistan using circumsporozoite protein (pvcsp) and merozoite surface protein-1 (pvmsp-1) genes as molecular markers.

Methods

The PCR for pvcsp and pvmsp-1 genes was carried out for 150 P. vivax isolates, followed by DNA sequencing of 35 and 30, respectively. Genetic diversity and polymorphism were analysed using ChromasPro, ClustalW, MEGA7, DnaSP v.5 and WebLogo programs.

Results

The PCR for pvcsp and pvmsp-1 genes was carried out for 150 P. vivax isolates and resulting the PCR products of 1100 bp for pvcsp and ~ 400 bp for pvmsp-1 genes, respectively. In the central-repeat region (CRR) of pvcsp gene, sequences comprised of four variable repeats of PRMs, out of which GDRADGQPA (PRM1), GDRAAGQPA (PRM2) were more extensively dispersed among the P. vivax isolates. Partial sequences (~ 400 bp) of block 2 of pvmsp-1 gene depicted high level of diversity.

Conclusion

The results revealed the polymorphism and genetic diversity especially at the CRR of pvcsp and block 2 of pvmsp-1 genes, respectively. The base-line data presented here warrants future studies to investigate more into the genetic diversity of P. vivax with large sample size from across the country for better understanding of population dynamics of P. vivax that will help to control malaria at individual and community level.

Balancing selection and high genetic diversity of Plasmodium vivax circumsporozoite central region in parasites from Brazilian Amazon and Rio de Janeiro Atlantic Forest

Circumsporozoite protein (CSP) is the primary pre-erythrocytic vaccine target in Plasmodium species. Knowledge about their genetic diversity can help predict vaccine efficacy and the spread of novel parasite variants. Thus, we investigated pvcsp gene polymorphisms in 219 isolates (136 from Brazilian Amazon [BA], 71 from Rio de Janeiro Atlantic Forest [AF], and 12 from non-Brazilian countries [NB]). Forty-eight polymorphic sites were detected, 46 in the central repeat region (CR), and two in the C-terminal region. Also, the CR presents InDels and a variable number of repeats. All samples correspond to the VK210 variant, and 24 VK210 subtypes based on CR. Nucleotide diversity (π = 0.0135) generated a significant number of haplotypes (168) with low genetic differentiation between the Brazilian regions (Fst = 0.208). The haplotype network revealed similar distances among the BA and AF regions. The linkage disequilibrium indicates that recombination does not seem to be acting in diversity, reinforcing natural selection's role in accelerating adaptive evolution. The high diversity (low Fst) and polymorphism frequencies could be indicators of balancing selection. Although malaria in BA and AF have distinct vector species and different host immune pressures, consistent genetic signature was found in two regions. The immunodominant B-cell epitope mapped in the CR varies from seven to 19 repeats. The CR T-cell epitope is conserved only in 39 samples. Concerning to C-terminal region, the Th2R epitope presented nonsynonymous SNP only in 6% of Brazilian samples, and the Th3R epitope remained conserved in all studied regions. We conclude that, although the uneven distribution of alleles may jeopardize the deployment of vaccines directed to a specific variable locus, a unique vaccine formulation could protect populations in all Brazilian regions.

Genetic polymorphism of vir genes of Plasmodium vivax in Myanmar

The Plasmodium vivax variant proteins encoded by vir genes are highly polymorphic antigens and are considered as one of key proteins of P. vivax for host immune evasion via antigenic variations. Because genetic diversity of these antigens is a critical hurdle in the development of an effective vaccine, understanding the genetic nature of the vir genes in natural population is important. In this study, we selected four vir genes (vir 4, vir 12, vir 21, and vir 27) previously used for genetic analysis in several studies and evaluated the genetic polymorphisms and phylogenetic relationship of these 4 vir genes in Myanmar P. vivax population. Taken all genetic diversity values, the vir 12 (S = 168, H = 17, Hd = 0.854, Tajima's D value = 2.91524) was the most genetically diverse gene and the vir 4 (S = 9, H = 4, Hd = 0.744, Tajima's D value = -0.49151) was the most conserved gene. All phylogenetic trees showed two clades, and vir 4 and 12 haplotypes from Myanmar were clustered in a distinct clade with those from India and Republic of Korea. These results confirmed the pattern of high genetic polymorphism of vir genes and provided information on vir gene for further functional research and studies focused toward the practical use of vir genes.

Genetic polymorphism and natural selection of circumsporozoite protein in Myanmar Plasmodium vivax

Background

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.

Methods

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.

Results

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.

Conclusion

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.