Next Article in Journal
NOV/CCN3 Promotes Cell Migration and Invasion in Intrahepatic Cholangiocarcinoma via miR-92a-3p
Previous Article in Journal
Genotype-Phenotype Correlations in PMM2-CDG
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Diversity of Plasmodium vivax Cysteine-Rich Protective Antigen (PvCyRPA) in Field Isolates from Five Different Areas of the Brazilian Amazon

by
Lana Bitencourt Chaves
1,
Glaucia de Oliveira Guimarães
1,
Daiana de Souza Perce-da-Silva
2,
Dalma Maria Banic
3,
Paulo Renato Rivas Totino
4,
Ricardo Luiz Dantas Machado
5,
Rodrigo Nunes Rodrigues-da-Silva
6,
Lilian Rose Pratt-Riccio
4,
Cláudio Tadeu Daniel-Ribeiro
4 and
Josué da Costa Lima-Junior
1,*
1
Laboratory of Immunoparasitology, Oswaldo Cruz Institute, Fiocruz, Rio de Janeiro 21040-900, Brazil
2
Laboratory of Basic and Applied Immunology, Arthur Sá Earp Neto University Center, Petrópolis 25680-120, Brazil
3
Laboratory of Clinical Immunology, Oswaldo Cruz Institute, Fiocruz, Rio de Janeiro 21040-900, Brazil
4
Laboratory of Malaria Research, Oswaldo Cruz Institute, Fiocruz, Rio de Janeiro 21040-900, Brazil
5
Department of Microbiology and Parasitology, Biomedical Institute, Fluminense Federal University, Rio de Janeiro 24210-130, Brazil
6
Laboratory of Monoclonal Antibodies Technology, Institute of Immunobiology Technology, Fiocruz, Rio de Janeiro 21040-900, Brazil
*
Author to whom correspondence should be addressed.
Genes 2021, 12(11), 1657; https://doi.org/10.3390/genes12111657
Submission received: 27 August 2021 / Revised: 27 September 2021 / Accepted: 27 September 2021 / Published: 21 October 2021
(This article belongs to the Section Human Genomics and Genetic Diseases)

Abstract

:
The Plasmodium vivax Cysteine-Rich Protective Antigen (PvCyRPA) has an important role in erythrocyte invasion and has been considered a target for vivax malaria vaccine development. Nonetheless, its genetic diversity remains uncharted in Brazilian malaria-endemic areas. Therefore, we investigated the pvcyrpa genetic polymorphism in 98 field isolates from the Brazilian Amazon and its impact on the antigenicity of predicted B-cell epitopes. Genetic diversity parameters, population genetic analysis, neutrality test and the median-joining network were analyzed, and the potential amino acid polymorphism participation in B-cell epitopes was investigated. One synonymous and 26 non-synonymous substitutions defined fifty haplotypes. The nucleotide diversity and Tajima’s D values varied across the coding gene. The exon-1 sequence had greater diversity than those of exon-2. Concerning the prediction analysis, seven sequences were predicted as linear B cell epitopes, the majority contained in conformational epitopes. Moreover, important amino acid polymorphism was detected in regions predicted to contain residues participating in B-cell epitopes. Our data suggest that the pvcyrpa gene presents a moderate polymorphism in the studied isolates and such polymorphisms alter amino acid sequences contained in potential B cell epitopes, an important observation considering the antigen potentiality as a vaccine candidate to cover distinct P. vivax endemic areas worldwide.

Graphical Abstract

1. Introduction

Malaria remains an important public health problem in several countries of tropical and subtropical regions of the world. In 2019, the disease caused an estimated 229 million clinical cases and around 409,000 deaths worldwide [1]. Among the Plasmodium species causing malaria in humans, Plasmodium vivax is the most widely distributed and prevalent outside of Africa [2]. In Brazil, endemic regions are restricted to the Legal Amazon, a region that currently accounts for the majority (>99%) of the countrywide malaria burden [3] and where P. vivax is predominant, with approximately 90% of the reported cases [4]. Several exclusive features of P. vivax biology, including the dormant liver stage, make it more resistant than other Plasmodium species to malaria elimination [5]. Thus, P. vivax presents a difficult obstacle to malaria elimination in endemic countries [6]. Therefore, it is very important to develop new methods and intervention strategies to block or reduce this transmission.
The complex life cycle of the Plasmodium includes an erythrocytic phase that is responsible for the clinical symptoms of malaria [7]. In this phase, P. vivax preferentially invades reticulocytes [8] in a process that occurs by sequential multiple molecule interactions, with each step mediated by antigens belonging to different protein families present on the merozoite surface and its apical organelles (i.e., micronemes and rhoptries) [9], which interact with a series of specific receptors on the erythrocyte surface to complete the invasion process [10]. The Cysteine-Rich Protective Antigen (CyRPA) is localized in the micronemes and is involved in the invasion process of merozoites into erythrocytes [11]. In Plasmodium falciparum, studies with the PfCyRPA protein have already demonstrated its inhibitory role in the invasion of merozoites both in vitro and in vivo, suggesting its potential as a candidate for asexual blood phase vaccine [12,13,14]. However, data demonstrating the potential of P. vivax CyRPA (PvCyRPA) as a vaccine candidate are still scarce and conflicting. França and collaborators demonstrated that antibodies against PvCyRPA are strongly related to protection. Interestingly, the protective effect of antibodies directed against PvCyRPA was higher than other proteins classically described as vaccine candidates, such as MSP-1, -3, -9 and AMA-1 [15]. On the other hand, in vitro studies of Ndegwa et al. (2021) showed that polyclonal antibodies raised against full-length PvCyRPA did not affect P. knowlesi growth [16].
A vaccine able to produce antibodies that effectively prevent the invasion process after the release of merozoites into the bloodstream may decrease parasite burden, disease symptoms and, indirectly, malaria transmission [17]. However, extensive allelic polymorphism in erythrocyte invasion pathways is known to limit the action of neutralizing antibodies against merozoite candidate vaccine antigens [18]. Malaria parasites have abundant genetic polymorphisms, much of which have evolved to escape host immune responses and thus present a major obstacle to the development of an effective malaria vaccine [19,20]. The genetic diversity and population structure of P. vivax for each candidate antigen is an important priority to the understanding of the malaria transmission dynamics [21]. In this scenario, many studies have been proposed to investigate the global diversity of leading vaccine antigens [22], and only one was recently addressed to the PvCyRPA [23], which does not include Brazilian malaria-endemic areas. Therefore, to understand the potential of PvCyRPA in vaccine development, we proposed to identify pvcyrpa gene in clinical isolates from different regions of the Brazilian Amazon and to study the potential impacts of the genetic diversity in predicted epitopes through bioinformatics tools.

2. Materials and Methods

2.1. Study Areas and Blood Sample Collection

Most cases of malaria in Brazil are concentrated in the Amazon Region, an endemic area for the disease [24]. For that reason, the study was carried out in five different regions of the Brazilian Amazon, with a set of 98 P. vivax-infected individuals previously described by our group: 31 from Cruzeiro do Sul, 17 from Mâncio Lima, 4 from Guajará, 37 from Manaus and 9 individuals from Oiapoque [25].
All 98 P. vivax participants were enrolled according to the following criteria: sought medical assistance for clinical malaria symptoms, presented uncomplicated malaria symptoms, were >18 years of age, and had a positive P. vivax malaria diagnosis. Pregnant women and P. falciparum-infected individuals were excluded from the study. Thin and thick blood smears were examined for the identification of the malaria parasite by a technician experienced in malaria diagnosis from the Brazilian Malaria Health Services. Thick blood smears from all subjects were stained with Giemsa, and a total of 200 microscopic fields were examined under a 1000-fold magnification. Thin blood smears of the positive samples were examined for species identification. To increase the sensitivity of parasite detection, molecular analyses using specific primers for genus (Plasmodium sp.) and species (P. falciparum and P. vivax) were performed in all the samples as previously described [26]. Donors positive for P. vivax and/or P. falciparum at the time of blood collection were subsequently treated by the chemotherapeutic regimen recommended by the Brazilian Ministry of Health.

2.2. Ethical Considerations

The study protocol was approved by the Research Ethics Committee of each locality, which included obtaining the following patients’ written consents for research use of their blood samples: Cruzeiro do Sul, Mâncio Lima and Guajará were reviewed and approved by the Fundação Oswaldo Cruz Research Ethics Committee, CEP-Fiocruz CAAE 46084015.1.0000.5248. In addition, the protocol of other blood sample collection was approved by the Research Ethics Committee of each locality: Manaus (CEP-Fiocruz): 346–613; Oiapoque (Hospital Municipal do Oiapoque/AP): 68980-000.

2.3. Genomic DNA Extraction

The DNA from 98 blood samples was previously extracted using the QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions and, then, stored at −20 °C until amplification.

2.4. Design of Pvcyrpa Specific Primers

The specific primers of pvcyrpa gene (1101 bp) were designed using standard gene sequences of P. vivax Salvador-1 (Sal-1) strain from GenBank NCBI Reference Sequence: XM_001615090.1 (Gene ID: PVX_090240). All oligonucleotides were designed and checked for specificity by using the Primer-BLAST tool provided by the National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/tools/primer-blast/, accessed on 13 August 2019) and the design quality of the oligonucleotides was evaluated by OligoAnalyzer v.3.1 (https://www.idtDNA.com/calc/analyzer/, accessed on 13 August 2019) to avoid homodimers and heterodimers (Table 1). The specific primers were chemically synthesized to perform PCR and DNA sequencing. The pvcyrpa gene has a structure consisting of two exons separated by a small well-conserved intron located on chromosome five, encoding for a microneme protein [22,27,28,29]. Consequently, from the extracted genomic DNA, the two exons of the pvcyrpa gene were amplified separately using two different primer sets, resulting in two gene lengths.

2.5. PCR Amplification of Pvcyrpa Gene

All the pvcyrpa genes reported in this study were amplified by conventional PCR using the two pairs of primers designed and described above. PCR reactions of the pvcyrpa gene were carried out in 25 μL volume that included 3 μL of DNA, 10 pmol/μL of each primer and the Master Mix kit (Promega, Madison, WI, USA) containing Taq DNA polymerase, PCR buffer and 10 nmol of each deoxynucleotide triphosphate (dNTP, Promega, Madison, WI, USA). The conventional PCR reaction was carried out using a GeneAmp PCR system 9700 (Applied Biosystems, Foster City, CA, USA) and the amplification conditions were as follows: for the first exon, one step at 95 °C for 2 min, 30 cycles at 95 °C for 1 min, 56 °C for 1 min and 72 °C for 1 min, and the last step at 72 °C for 1 min. For the second exon, the temperatures and the number of cycles remained the same, except for the annealing temperature, which was 59 °C for 1 min. The PCR conditions for the amplification of PvCyRPA F1/R1 generated the first pvcyrpa gene fragment (600 bp) and the PCR conditions for amplification of PvCyRPA F2/R2 yielded the second pvcyrpa gene fragment (466 bp). In all reactions, two negative controls (one without DNA and the other with DNA extracted from in vitro culture of P. falciparum PSS1 strain) and positive control (P. vivax-infected sample) were used. To confirm the presence of DNA from the in vitro culture of P. falciparum and that the lack of amplification was due to the specificity of the primers for pvcyrpa, we performed the amplification of the P. falciparum p126 gene fragment and electrophoresis as previously described [30]. After PCR, ten μL of amplified products were size-fractionated by electrophoresis within 2% agarose gel (Sigma Aldrich, Missouri, USA) in 1× TAE buffer (0.04 M TRIS-acetate, 1 mM EDTA) in the presence of 1× GelRed nucleic acid stain (Biotium, Fremont, CA, USA). PCR products were visualized by ultraviolet (UV) illumination. The sizing of products was performed using a GeneRuler 100 bp Plus DNA Ladder (Thermo Scientific, Waltham, MA, USA). Then, amplicons were purified using the GE Healthcare Lifesciences kit following the manufacturer’s instructions. Afterward, 5–50 ng of DNA was used per sequencing reaction employing the Sanger method, using forward and reverse primers.

2.6. DNA Sequencing and Polymorphism Analysis

The specificity of the assay was confirmed by sequencing the PCR products from all positive samples using a Big Dye Terminator Sequencing Kit (Applied Biosystems, Foster City, CA, USA) following the manufacturer’s instructions. The DNA sequencing was carried out on the ABI 3730xl DNA analyzer (Applied Biosystems, Foster City, CA, USA) with the support of Fiocruz Genomic Platform and all the results were analyzed using DNASTAR’s sequence alignment software [31]. Moreover, the sequences were also analyzed in BioEdit sequence alignment editor to better-visualized SNP positions, employing ClustalW multiple sequence alignment and the Sal-1 strain as a reference sequence.
To analyze hypothetical PvCyRPA protein derivatives from P. vivax genome data available and Mexico isolates already available, multiple sequence alignments were conducted with Clustal Omega using the MegAlign Pro 15 (Lasergene DNASTAR) program. The following sequences were used: India VII—GenBank: gb|KMZ81773.1; North Korean—GenBank: gb|KNA00954.1; Mauritania I—Genbank: gb|KMZ94334.1; Brazil I—GenBank: gb|KMZ87926.1; Sanger Institute: SCO 66052.1; Sanger Institute: SCO71483.1; Sanger Institute: SGX76259.1; Mexico-Southern Mexican [23]. Additionally, our Brazilian Amazon isolates were compared to P01 strain (PVP01_0532400), a new reference genome for P. vivax from an Indonesian clinical isolate [32], and the circular map of protein alignment was generated using the software GenVision v15 (Lasergene DNASTAR).
Multiple alignment Clustal Omega, distance matrix, and the phylogenetic tree were conducted using the MegAlign Pro 15 (Lasergene DNASTAR) program and the circular map of protein alignment was generated using the software GenVision v15 (Lasergene DNASTAR).

2.7. Genetic Analysis of the Coding Gene

Genetic diversity of pvcyrpa sequences was analyzed using the DnaSP v6 software [33] to estimate within-population diversity based on the genetic diversity parameters as the number of segregation sites (S), the number of haplotypes (h), haplotype diversity (Hd), and nucleotide diversity (π).
Natural selection in pvcyrpa was assessed by the Tajima’s D and Z-test. To test the neutral theory of evolution, Tajima’s D values [34] were calculated using the total number of mutations also estimated with DnaSP v.6 software. This test informs about the selection and demographic forces acting on a population. Positive values might be suggestive of positive or balancing selection. This force maintains alleles at balanced frequencies. On the other hand, negative values suggest purifying selection or recent population expansion [34]. The Z-test method was performed with MEGA7 v.6.0; the rates of non-synonymous (dN) to synonymous (dS) substitutions (dN/dS) (1000 bootstrapping replicates) were estimated with Nei and Gojobor’s method [35] and with the Jukes and Cantor correction, in which p < 0.05 was considered significant.
Haplotype data also were generated using DnaSP v.6 and the haplotype network was constructed using PopArt v.1.7 with the median-joining algorithm [36] to explore the parasite relationships based on the pvcyrpa gene. Mutational steps represent the connections between haplotypes, and empty squares show the non-sampled or extinct haplotypes. The color of the circles represents the geographic origins of each haplotype, while the size of the circle represents the frequency of each haplotype.

2.8. Prediction of Linear B-Cell Epitopes

The prediction of linear B-cell epitopes was carried out using the Ellipro algorithm, as well as confirmed by the overlap between the predictions of at least two more algorithms (BCPred, BepiPred, ABCpred and Emini). This software takes a single sequence in FASTA format input and each amino acid receives a prediction score profile of known antigens and incorporates propensity scale methods based on hydrophilicity and secondary structure prediction. For each input sequence, the server outputs a prediction score. The positions of the linear B-cell epitopes are predicted to be located at the residues with the highest scores. In addition, the software ElliPro predicts linear and discontinuous antibody epitopes based on a protein antigen’s 3D structure and accepts two types of input data: protein sequence or structure (PDB format) [37]. This server associates each predicted epitope with a score, defined as a PI (Protrusion Index) value averaged over epitope residues. In the method, residues with larger scores are associated with greater solvent accessibility.

3. Results

3.1. Molecular Characterization of the Pvcyrpa Gene in the Studied Regions

To identify the gene encoding the PvCyRPA in isolates from Brazilian endemic areas, 98 blood samples from P. vivax-infected individuals living in the cities of Cruzeiro do Sul, Mâncio Lima, Guajará, Manaus e Oiapoque had the DNA extracted and subjected to molecular diagnosis by conventional PCR.
The pvcyrpa sequence encodes 1101 bp with the two exons sequence and so it was divided into two regions (exon-1 and exon-2). The primer combinations (PvCyRPA_F1/R1 and PvCyRPA_F2/R2) designed to cover the length of the targeted pvcyrpa gene resulted in amplification of two fragments in 100% of samples (Figure 1). PvCyRPA F1/R1 primer combination amplified a fragment of 600 bp whereas PvCyRPA F2/R2 primer combination amplified a fragment of 466 bp. Additionally, P. falciparum specimens from in vitro culture were tested for quality assurance, resulting in negative PCR amplification of the pvcyrpa gene (Figure 1B). Therefore, the 98 samples from individuals infected with P. vivax amplified by PCR were subjected to sequencing reactions to screen the possible single nucleotide polymorphisms of the pvcyrpa gene. All amplified fragments were sequenced and aligned for sequence analysis.

3.2. Genetic Diversity of the Pvcyrpa Gene

Compared to reference sequence Sal-1 (Gene ID: PVX_090240), 27 polymorphic sites were observed, of which one was synonymous and 26 were non-synonymous substitutions. In addition, the isolates presented a frequent deletion of a single codon in position E267 (Table 2). Of these non-synonymous substitutions, two amino acids positions showed two variant alleles each—Q142 (Q142K and Q142R) and D145 (D145G and D145N)—that sometimes were present together in one population. Cruzeiro do Sul and Mâncio Lima presented the two variants existing in position Q142. As well as Cruzeiro do Sul, Mâncio Lima, Manaus and Oiapoque also presented the two variants in position D145, except Guajará. The Cruzeiro do Sul and Mâncio Lima regions shared all the 26 non-synonymous substitutions. Overall, R122K (N = 80%; 82%), K131E (N = 77%; 79%), D149G (N = 62%; 63%), A154D (N = 60%; 61%) and E159D (N = 66%; 67%) SNPs were the most frequent in Brazilian Amazon isolates, while E86 (N = 91%; 93%), Q142 (N = 95%; 97%), A187 (N = 92%; 94%) and V287 (N = 91%; 93%) presented high frequency of similarity to the wild type sequence. N (%): Frequency of substitution in each locality (Table 2 and Figure 2).

3.3. Population Genetic Analysis

In P. vivax isolates from the Brazilian Amazon, the genetic diversity was heterogeneously distributed between the regions coding the two exons, with higher values for exon-1 comparing to exon-2 among localities (Table 3). Using the entire coding gene, it was shown that Mâncio Lima isolates had the highest nucleotide diversity (π) (0.01272 ± 0.00064), while Manaus had the lowest (0.01116 ± 0.00030). Oiapoque had the highest haplotype diversity (Hd) (1.000 ± 0.052), while Guajará had the lowest (0.833 ± 0.222). Some similar differences were detected when analyzing each exon separately. Concerning exon-1, the highest nucleotide diversity was also observed in the Mâncio Lima group (0.01517 ± 0.00095) among all five populations, while Manaus sequences displayed the lowest nucleotide diversity (0.01285 ± 0.00050). Moreover, parasites from Oiapoque presented the highest estimate of haplotype diversity (0.917 ± 0.092), whereas parasites from Guajará showed the lowest Hd (0.667 ± 0.204). In comparison to exon-2, the highest nucleotide diversity was observed in the Mâncio Lima group (0.00938 ± 0.00077). In contrast, Guajará sequences displayed the lowest nucleotide diversity (0.00792 ± 0.00420) in exon-2. Isolates from Oiapoque presented the highest estimate of haplotype diversity (Hd) (0.889 ± 0.091), whereas parasites from Guajará showed the lowest Hd (0.500 ± 0.265).
We performed the Tajima’s D and Z-test to determine whether natural selection was affecting the pvcyrpa gene. The Tajima’s D (TjD) test was performed to assess if there is selective pressure on the pvcyrpa gene and it was calculated using the total number of mutations. In coding regions, an excess of non-synonymous relative to synonymous changes suggests a clear signal of positive selection. The TjD test showed significant positive values for entire coding in Cruzeiro do Sul (2.24953, * p < 0.05) and Manaus (2.75595, ** p < 0.01). Similarly, for exon-1 in Cruzeiro do Sul (2.17555, * p < 0.05), Guajará (2.24818, * p < 0.05) and Manaus (2.60865, ** p < 0.01) showed significant positives values, while the other locations showed no significant positive values. In contrast, the TjD test presented a significant positive value at exon-2 only in Manaus (2.42821, * p < 0.05). The Z test of selection shows that the exon-1 sequence had positive values in parasites of different origins. These values were higher in all parasites except Guajará (Table 3).

3.4. Haplotype Network Analysis

The median-joining haplotype network constructed by PopArt 1.7 using the 98 sequences produced a total of 50 haplotypes with some of which consisting of more than one sequence from the Brazilian Amazon (Figure 3). The haplotype network to explore the parasite relationships based on the pvcyrpa gene and comprising mutations at 34 segregating sites. All 50 haplotypes were found closely related. The haplotypes Hap_1 and Hap_11 had high frequency and shared parasites from all five localities. The haplotype Hap_4 (Cruzeiro do Sul, Mâncio Lima and Guajará) and Hap_8 (Cruzeiro do Sul, Mâncio Lima and Manaus) both shared parasites from three localities. Moreover, the haplotypes Hap_2 (Cruzeiro do Sul and Mâncio Lima) and Hap_9 (Cruzeiro do Sul and Manaus), Hap_14 (Cruzeiro do Sul and Mâncio Lima), Hap_23 (Manaus and Oiapoque) shared sequences from 2 localities. The other haplotypes had sequences from only one location.

3.5. Comparison of Amino Acid Variations in PvCyRPA among Genome Sequences Available Worldwide

The PvCyRPA amino acid substitutions identified in genome sequences worldwide, including those from Brazilian Amazon, are resumed in Table 4. As observed in the protein sequence alignments, the PvCyRPA coding gene had an excess of non-synonymous mutations, which were more frequent in exon-1 than in exon-2. In addition, we subsequently aligned the protein sequence of these mutant field isolates with other hypothetical CyRPA proteins derivative from P. vivax genome data available in the GenBank database and also aligned with the isolate Mexico-Southern Mexican [23]. Among a total of 31 in PvCyRPA protein observed in P. vivax sequences worldwide, 26 amino acid substitutions are also present in our isolates. Interestingly, only our isolates showed a new substitution at K150R, being found in the localities of Cruzeiro do Sul, Mâncio Lima, and Guajará. In addition, we can observe a high genetic variability among the isolates of each locality. Curiously, Cruzeiro do Sul and Mâncio Lima present the two variants existing in position Q142 when compared to the genome sequences. Additionally, Cruzeiro do Sul, Mâncio Lima, Manaus and Oiapoque also presented the two variants in position D145, except Guajará. Likewise, some amino acid substitutions of the PvCyRPA protein were found but were rare in only some genomes: L180H in India VII; I63T and Y361H in North Korea; R125T and Q147K in SCO 66052.1 (Sanger Institute) (Table 4). Furthermore, we also compared the consensus sequences of our Brazilian Amazon isolates with reference sequence Sal-1 (PVX_090240) and P01 strain (PVP01_0532400) and generated the circular map of protein alignment using the software GenVision v15 (Lasergene DNASTAR) (Figure 4A). As expected, a significantly high degree of identity was observed across the sequences analyzed, maintaining the mutations found in relation to reference sequence Sal-1. The analysis showed a high identity among our isolates and P01 strain, despite the deletion of 4 amino acids present in the P01 sequence at positions 13–16 (FLFS). According to pairwise distance, the percent identity ranged from 94.5% (P01 vs. GJ) to 99.7% (CZS vs. OIA) (Figure 4A,B).

3.6. Polymorphisms and Potential B-Cell Epitopes

We performed in silico prediction for the identification of B cell epitopes present in the PvCyRPA protein using Sal-1 reference and then, seven amino acid sequences were predicted as linear epitopes, (Table 5). All epitopes were initially predicted by the Ellipro algorithm and were confirmed by the overlap between the predictions of at least two more algorithms (BCPred, BepiPred, ABCpred and Emini). The sequences varied from 9 to 19 amino acids and five sequences were inserted in conformational epitopes. The protein appears to have epitopes mainly in the central region and the C terminal region, while the N region terminal does not contain antigenic sequences. Of note, amino acid polymorphisms were detected in regions predicted to contain residues participating in B-cell epitopes (Table 6).

4. Discussion

The invasion of the red blood cell by Plasmodium merozoites is essential for parasite survival and proliferation. The merozoites have therefore evolved multiple pathways, using various antigenic proteins which aid in the invasion process. Among the merozoite’s invasive proteins are Cysteine-Rich Protective Antigen (CyRPA), which seems to be essential for the parasite’s life cycle during the invasion of erythrocytes and a ligand for reticulocyte invasion [38]. The discovery of the antigen has revamped hope in the search for an effective malaria blood-stage vaccine of P. vivax. However, one of the major obstacles to malaria vaccine development is still the low efficiency of proteins used as immunogens in inducing protection, which, in part, can be explained by genetic polymorphisms [39]. It is important to understand the mechanisms of genetic recombination and sequence variation that represent the repertoire of polymorphic malarial surface antigens and that may help in designing vaccines [29,40]. The genetic diversity of these proteins in hyperendemic areas has been described as a limiting factor for the rapid acquisition of protective immunity and, consequently for the development of an effective vaccine. Furthermore, the antigenic polymorphism of P. vivax vaccine candidates has been little discussed in unstable transmission areas such as the Brazilian endemic regions [41]. Thus, considering that the epidemiology of malaria in Brazil presents unstable transmission and the knowledge about the genetic polymorphism of pvcyrpa remains unknown, we aimed to identify the pvcyrpa gene in isolates from different regions of the Brazilian Amazon and to study the potential impacts of the genetic diversity in potential B-cell epitopes.
The identification and analysis of the genetic diversity of the pvcyrpa gene in isolates from different geographic regions of the Brazilian Amazon have not been previously studied. Considering the distance among the studied localities and the possible existence of a gene flow of Plasmodium vivax genome among the studied populations, associated with migration of people, could promote the gene flow of the parasite [22] and impact the parasite transmission and dispersion [42,43]. Our first results showed that the pvcyrpa gene has high genetic variability in relation to reference sequence Sal-1, presenting 27 polymorphic sites of which one was synonymous and 26 non-synonymous substitutions throughout the sequence. Among these non-synonymous substitutions, two amino acids positions—Q142 (Q142K and Q142R) and D145 (D145G and D145N)—presented one or two variants in our study areas. Overall, R122K (N = 80%; 82%), K131E (N = 77%; 79%), D149G (N = 62%; 63%), A154D (N = 60%; 61%) and E159D (N = 66%; 67%) mutations were the most frequent in our Brazilian Amazon isolates. The analysis of the pvcyrpa gene from the Brazilian Amazon showed that mutations have contributed to generating nucleotide and haplotype diversity. The similarity in the genetic diversity pattern suggests that similar evolutionary forces act on pvcyrpa parasites and that the structural and/or functional properties are consistent. To evade the immune response, genes encoding for antigenic proteins accumulate non-synonymous mutations, which leads to an increase in diversity. In this study, the pvcyrpa gene presented non-synonymous mutation accumulation in parasites of different regions, mainly in exon-1. Significant positive values of Tajima’s D indicate balancing selection and population bottlenecks, while negative values suggest the presence of purifying selection or population expansion [34]. Exon-1 had significant positive values in Cruzeiro do Sul, Guajará and Manaus for the Tajima’s Test (TjD) as well as in Manaus at exon-2. The results suggest that polymorphism at pvcyrpa exon-1 is generated by mutation and recombination, and is probably maintained by positive balancing selection pressure, which might represent an evolutionary advantage to the parasite. Exon-1 codes for highly variant domains exposed on the surface of infected red blood cells (RBCs), while exon-2 codes for the more conserved segment [29]. Furthermore, the level of genetic diversity in blood-stage antigens seems to be associated with the degree of exposure to the immune system [29,44].
P. vivax biological and genetic characteristics, host immunity, and local vectors may contribute to their different patterns of demographic expansion [45]. Some discrete P. vivax lineages can remain stable across time in one of the areas with the highest malaria transmission in the Americas. Relapses can account for some clonal persistence because P. vivax strains are repeatedly reintroduced in the population as hypnozoites reactivate [46]. Maybe this context can be to explain why only Mâncio Lima and Oiapoque showed no significant positive values of Tajima’s D. However, genomic epidemiology approaches can help better to reveal the complex distribution of this parasite in the Brazilian Amazon, as well as the relationships with the worldwide genetic diversity.
Moreover, it was possible to identify 50 different haplotypes of pvcyrpa gene among the 98 P. vivax field isolates from the regions that were analyzed. The haplotype network explores the parasite relationships based on the pvcyrpa gene and comprising mutations at 34 segregating sites, and confirmed the extensive genetic diversity observed in pvcyrpa sequences. All 50 haplotypes were found to be closely related, with some of which consisting of more than one sequence from the Brazilian Amazon. Regarding the pvcyrpa sequences, we observed that haplotypes Hap_1 and Hap_11 had high frequency and shared parasites from all five localities. These findings suggest a global distribution of parasites containing similar pvcyrpa genotypes. Additionally, to compare our findings with the PvCyRPA sequences around the world, we observed that it presents a similar genetic profile among the complete genomes of P. vivax available on the GenBank Database. Among a total of 31 amino acid substitutions of PvCyRPA protein observed in P. vivax sequences worldwide, 26 amino acid substitutions are also present in our isolates.
To develop an effective malaria vaccine that can work in different regions of the world, it is important to include alleles that can induce the immune response and cover the antigenic diversity of P. vivax population. Consequently, the existence of the same haplotypes in different malaria-endemic areas and similar genetic profiles worldwide in their results will be important for the rationale of malaria vaccine designs. Moreover, as the immune system could act as selective pressure and the PvCyRPA is emerging as an alternative antigen in vaccine development, we also evaluated the impact of non-synonymous polymorphisms in relation to predicted B-cell epitope sequences.
Amino acid variation was present at peptide regions potentially participating in B-cell epitopes, which supports the idea that this molecule is under selective immune pressure. The seven B-cell potential epitopes have been identified in the PvCyRPA protein, most of which are contained in conformational epitopes, which corroborates its potential as an antibody target. A characteristic of malaria blood-stage antigens is their participation in merozoite invasion and immune evasion. Immunogenicity studies and molecular modeling are essential to determine the importance of PvCyRPA as a vaccine candidate. Targeting molecules important for the Plasmodium life cycle might be limited by their antigenic polymorphism or low immunogenicity. Molecular studies provide information about the dynamics of vaccine antigen polymorphisms that can be used to make informed decisions about which parasite alleles to include in vaccine formulations, and to evaluate accurately the efficacy of vaccines tested in malaria-endemic areas [21]. Thus, an effective antigen vaccine should include alleles that induce host immune responses that are sufficiently broad to cover the existing antigenic diversity. Nevertheless, because of the higher genetic diversity of P. vivax compared to P. falciparum, generating a broad cross-reactive immune response against highly polymorphic asexual stage antigens faces even greater challenges [47].

5. Conclusions

In summary, the present study explored the genetic polymorphism of PvCyRPA in field isolates from distinct endemic areas in Brazil, showing a moderate sequence variation, which could influence the potential B-cell epitopes and, consequently, antibody recognition. Despite the observed amino acid changes in the studied population and sequences worldwide, the potential antibody targets did not seem to be significantly affected. However, due to the paucity of information on PvCyRPA genetic diversity and its potential as a vaccine candidate, more studies are necessary to confirm the impact of PvCyRPA polymorphism in naturally acquired immune response and/or vaccine development.

Author Contributions

Conceptualization, L.B.C. and J.d.C.L.-J.; Methodology, L.B.C. and G.d.O.G.; software, L.B.C. and J.d.C.L.-J.; validation, L.B.C.; formal analysis, L.B.C. and J.d.C.L.-J.; investigation, L.B.C.; resources, J.d.C.L.-J., D.M.B., C.T.D.-R. and R.L.D.M.; Visualization, R.N.R.-d.-S.; writing—original draft preparation, L.B.C. and J.d.C.L.-J.; writing—review and editing, D.d.S.P.-d.-S., D.M.B., P.R.R.T., C.T.D.-R., R.L.D.M., L.R.P.-R., C.T.D.-R. and J.d.C.L.-J.; supervision, J.d.C.L.-J.; project administration, J.d.C.L.-J.; funding acquisition, J.d.C.L.-J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by FAPERJ (E26/202.854/2019—JCNE and E26/211.112/2019—Apoio a Grupos Emergentes de Pesquisa), Inovabio grant program (BIO-004-FIO-17-2-11—Bio-Manguinhos/Fiocruz) and FIOCRUZ (PAEF II, IOC-023-FIO-18-2-47). JCLJ and CTDR receive a Productivity Research Fellowship CNPq; JCLJ and PRRT are Jovem Cientista do Nosso Estado and CTDR is Cientista do Nosso Estado—FAPERJ. The Laboratório de Pesquisa em Malária (LPM, IOC, Fiocruz) is a National Institutes of Science & Technology (INCT) Associated Laboratory. The LPM is also an Associated Laboratory of the Neuroinflammation Network (FAPERJ). The APC was funded by Fiocruz.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Fundação Oswaldo Cruz Research Ethics Committee, CEP-Fiocruz CAAE 46084015.1.0000.5248/2016.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data supporting the conclusions of this article are provided within the article.

Acknowledgments

We are grateful to all individuals who participated in this study for their cooperation and generous donation of blood, which made this study possible. We are also indebted to the sequencing platform of the Program for Technological Development of Health Products (PDTIS-Fiocruz) for use of their facilities.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. WHO. World Malaria Report 2020; World Health Organization: Geneva, Switzerland, 2020. [Google Scholar]
  2. Dalrymple, U.; Mappin, B.; Gething, P.W. Malaria mapping: Understanding the global endemicity of falciparum and vivax malaria. BMC Med. 2015, 13, 140. [Google Scholar] [CrossRef] [Green Version]
  3. da Silva-Nunes, M.; Moreno, M.; Conn, J.E.; Gamboa, D.; Abeles, S.; Vinetz, J.M.; Ferreira, M.U. Amazonian malaria: Asymptomatic human reservoirs, diagnostic challenges, environmentally driven changes in mosquito vector populations, and the mandate for sustainable control strategies. Acta Trop. 2012, 121, 281–291. [Google Scholar] [CrossRef] [Green Version]
  4. Saúde SdVeSMd. Boletim Epidemiológico Malária; Secretaria de Vigilância em Saúde | Ministério da Saúde: Brasilia, Brazil, 2021. [Google Scholar]
  5. Hulden, L.; Hulden, L. Activation of the hypnozoite: A part of Plasmodium vivax life cycle and survival. Malar. J. 2011, 10, 90. [Google Scholar] [CrossRef] [Green Version]
  6. Cotter, C.; Sturrock, H.J.; Hsiang, M.S.; Liu, J.; Phillips, A.A.; Hwang, J.; Gueye, C.S.; Fullman, N.; Gosling, R.D.; Feachem, R.G. The changing epidemiology of malaria elimination: New strategies for new challenges. Lancet 2013, 382, 900–911. [Google Scholar] [CrossRef]
  7. Miller, L.H.; Baruch, D.I.; Marsh, K.; Doumbo, O.K. The pathogenic basis of malaria. Nature 2002, 415, 673–679. [Google Scholar] [CrossRef] [PubMed]
  8. Mons, B.; Collins, W.E.; Skinner, J.C.; van der Star, W.; Croon, J.J.; van der Kaay, H.J. Plasmodium vivax: In vitro growth and reinvasion in red blood cells of Aotus nancymai. Exp. Parasitol. 1988, 66, 183–188. [Google Scholar] [CrossRef]
  9. Kanjee, U.; Rangel, G.W.; Clark, M.A.; Duraisingh, M.T. Molecular and cellular interactions defining the tropism of Plasmodium vivax for reticulocytes. Curr. Opin. Microbiol. 2018, 46, 109–115. [Google Scholar] [CrossRef]
  10. Mitchell, G.H.; Bannister, L.H. Malaria parasite invasion: Interactions with the red cell membrane. Crit. Rev. Oncol./Hematol. 1988, 8, 225–310. [Google Scholar] [CrossRef]
  11. Favuzza, P.; Guffart, E.; Tamborrini, M.; Scherer, B.; Dreyer, A.M.; Rufer, A.C.; Erny, J.; Hoernschemeyer, J.; Thoma, R.; Schmid, G.; et al. Structure of the malaria vaccine candidate antigen CyRPA and its complex with a parasite invasion inhibitory antibody. eLife 2017, 6, e20383. [Google Scholar] [CrossRef] [Green Version]
  12. Tamborrini, M.; Hauser, J.; Schafer, A.; Amacker, M.; Favuzza, P.; Kyungtak, K.; Fleury, S.; Pluschke, G. Vaccination with virosomally formulated recombinant CyRPA elicits protective antibodies against Plasmodium falciparum parasites in preclinical in vitro and in vivo models. NPJ Vaccines 2020, 5, 9. [Google Scholar] [CrossRef] [Green Version]
  13. Favuzza, P.; Blaser, S.; Dreyer, A.M.; Riccio, G.; Tamborrini, M.; Thoma, R.; Matile, H.; Pluschke, G. Generation of Plasmodium falciparum parasite-inhibitory antibodies by immunization with recombinantly-expressed CyRPA. Malar. J. 2016, 15, 161. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Dreyer, A.M.; Matile, H.; Papastogiannidis, P.; Kamber, J.; Favuzza, P.; Voss, T.S.; Wittlin, S.; Pluschke, G. Passive immunoprotection of Plasmodium falciparum-infected mice designates the CyRPA as candidate malaria vaccine antigen. J. Immunol. 2012, 188, 6225–6237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Franca, C.T.; White, M.T.; He, W.Q.; Hostetler, J.B.; Brewster, J.; Frato, G.; Malhotra, I.; Gruszczyk, J.; Huon, C.; Lin, E.; et al. Identification of highly-protective combinations of Plasmodium vivax recombinant proteins for vaccine development. eLife 2017, 6, e28673. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Ndegwa, D.N.; Kundu, P.; Hostetler, J.B.; Marin-Menendez, A.; Sanderson, T.; Mwikali, K.; Verzier, L.H.; Coyle, R.; Adjalley, S.; Rayner, J.C. Using Plasmodium knowlesi as a model for screening Plasmodium vivax blood-stage malaria vaccine targets reveals new candidates. PLoS Pathog. 2021, 17, e1008864. [Google Scholar] [CrossRef] [PubMed]
  17. Beeson, J.G.; Drew, D.R.; Boyle, M.J.; Feng, G.; Fowkes, F.J.; Richards, J.S. Merozoite surface proteins in red blood cell invasion, immunity and vaccines against malaria. FEMS Microbiol. Rev. 2016, 40, 343–372. [Google Scholar] [CrossRef] [Green Version]
  18. Crompton, P.D.; Pierce, S.K.; Miller, L.H. Advances and challenges in malaria vaccine development. J. Clin. Investig. 2010, 120, 4168–4178. [Google Scholar] [CrossRef]
  19. Barry, A.E.; Arnott, A. Strategies for designing and monitoring malaria vaccines targeting diverse antigens. Front. Immunol. 2014, 5, 359. [Google Scholar] [CrossRef] [Green Version]
  20. The malERA Consultative Group on Vaccines, A research agenda for malaria eradication: Vaccines. PLoS Med. 2011, 8, e1000398. [CrossRef]
  21. Takala, S.L.; Plowe, C.V. Genetic diversity and malaria vaccine design, testing and efficacy: Preventing and overcoming ‘vaccine resistant malaria’. Parasite Immunol. 2009, 31, 560–573. [Google Scholar] [CrossRef] [Green Version]
  22. de Souza-Neiras, W.C.; de Melo, L.M.; Machado, R.L. The genetic diversity of Plasmodium vivax—A review. Mem. Inst. Oswaldo Cruz 2007, 102, 245–254. [Google Scholar] [CrossRef] [Green Version]
  23. Gonzalez-Ceron, L.; Cebrian-Carmona, J.; Mesa-Valle, C.M.; Garcia-Maroto, F.; Santillan-Valenzuela, F.; Garrido-Cardenas, J.A. Plasmodium vivax Cysteine-Rich Protective Antigen Polymorphism at Exon-1 Shows Recombination and Signatures of Balancing Selection. Genes 2020, 12, 29. [Google Scholar] [CrossRef] [PubMed]
  24. Oliveira-Ferreira, J.; Lacerda, M.V.; Brasil, P.; Ladislau, J.L.; Tauil, P.L.; Daniel-Ribeiro, C.T. Malaria in Brazil: An overview. Malar. J. 2010, 9, 115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Chaves, L.B.; Perce-da-Silva, D.S.; Totino, P.R.R.; Riccio, E.K.P.; Baptista, B.O.; de Souza, A.B.L.; Rodrigues-da-Silva, R.N.; Machado, R.L.D.; de Souza, R.M.; Daniel-Ribeiro, C.T.; et al. Plasmodium vivax ookinete surface protein (Pvs25) is highly conserved among field isolates from five different regions of the Brazilian Amazon. Infect. Genet. Evol. J. Mol. Epidemiol. Evol. Genet. Infect. Dis. 2019, 73, 287–294. [Google Scholar] [CrossRef] [PubMed]
  26. Snounou, G.; Viriyakosol, S.; Zhu, X.P.; Jarra, W.; Pinheiro, L.; do Rosario, V.E.; Thaithong, S.; Brown, K.N. High sensitivity of detection of human malaria parasites by the use of nested polymerase chain reaction. Mol. Biochem. Parasitol. 1993, 61, 315–320. [Google Scholar] [CrossRef]
  27. Garrido-Cardenas, J.A.; Gonzalez-Ceron, L.; Manzano-Agugliaro, F.; Mesa-Valle, C. Plasmodium genomics: An approach for learning about and ending human malaria. Parasitol. Res. 2019, 118, 1–27. [Google Scholar] [CrossRef]
  28. Singh, V.; Gupta, P.; Pande, V. Revisiting the multigene families: Plasmodium var and vir genes. J. Vector Borne Dis. 2014, 51, 75–81. [Google Scholar]
  29. Ferreira, M.U.; da Silva Nunes, M.; Wunderlich, G. Antigenic diversity and immune evasion by malaria parasites. Clin. Diagn. Lab. Immunol. 2004, 11, 987–995. [Google Scholar] [CrossRef] [Green Version]
  30. Pratt-Riccio, L.R.; Sallenave-Sales, S.; de Oliveira-Ferreira, J.; da Silva, B.T.; Guimaraes, M.L.; Santos, F.; de Simone, T.S.; Morgado, M.G.; de Simone, S.G.; Ferreira-Da-Cruz Mde, F.; et al. Evaluation of the genetic polymorphism of Plasmodium falciparum P126 protein (SERA or SERP) and its influence on naturally acquired specific antibody responses in malaria-infected individuals living in the Brazilian Amazon. Malar. J. 2008, 7, 144. [Google Scholar] [CrossRef] [Green Version]
  31. Burland, T.G. DNASTAR’s Lasergene sequence analysis software. Methods Mol. Biol. 2000, 132, 71–91. [Google Scholar] [CrossRef] [PubMed]
  32. Auburn, S.; Bohme, U.; Steinbiss, S.; Trimarsanto, H.; Hostetler, J.; Sanders, M.; Gao, Q.; Nosten, F.; Newbold, C.I.; Berriman, M.; et al. A new Plasmodium vivax reference sequence with improved assembly of the subtelomeres reveals an abundance of pir genes. Welcome Open Res. 2016, 1, 4. [Google Scholar] [CrossRef]
  33. Rozas, J.; Ferrer-Mata, A.; Sanchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sanchez-Gracia, A. DnaSP 6: DNA Sequence Polymorphism Analysis of Large Data Sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef] [PubMed]
  34. Tajima, F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 1989, 123, 585–595. [Google Scholar] [CrossRef] [PubMed]
  35. Nei, M.; Gojobori, T. Simple methods for estimating the numbers of synonymous and non-synonymous nucleotide substitutions. Mol. Biol. Evol. 1986, 3, 418–426. [Google Scholar] [CrossRef]
  36. Bandelt, H.J.; Forster, P.; Rohl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 1999, 16, 37–48. [Google Scholar] [CrossRef]
  37. Ponomarenko, J.; Bui, H.H.; Li, W.; Fusseder, N.; Bourne, P.E.; Sette, A.; Peters, B. ElliPro: A new structure-based tool for the prediction of antibody epitopes. BMC Bioinform. 2008, 9, 514. [Google Scholar] [CrossRef] [Green Version]
  38. Knuepfer, E.; Wright, K.E.; Kumar Prajapati, S.; Rawlinson, T.A.; Mohring, F.; Koch, M.; Lyth, O.R.; Howell, S.A.; Villasis, E.; Snijders, A.P.; et al. Divergent roles for the RH5 complex components, CyRPA and RIPR in human-infective malaria parasites. PLoS Pathog. 2019, 15, e1007809. [Google Scholar] [CrossRef]
  39. Miller, L.H.; Roberts, T.; Shahabuddin, M.; McCutchan, T.F. Analysis of sequence diversity in the Plasmodium falciparum merozoite surface protein-1 (MSP-1). Mol. Biochem. Parasitol. 1993, 59, 1–14. [Google Scholar] [CrossRef]
  40. Arnott, A.; Barry, A.E.; Reeder, J.C. Understanding the population genetics of Plasmodium vivax is essential for malaria control and elimination. Malar. J. 2012, 11, 14. [Google Scholar] [CrossRef] [Green Version]
  41. Ferreira, M.U.; Karunaweera, N.D.; da Silva-Nunes, M.; da Silva, N.S.; Wirth, D.F.; Hartl, D.L. Population structure and transmission dynamics of Plasmodium vivax in rural Amazonia. J. Infect. Dis. 2007, 195, 1218–1226. [Google Scholar] [CrossRef] [Green Version]
  42. Ferreira, M.U.; Castro, M.C. Challenges for malaria elimination in Brazil. Malar. J. 2016, 15, 284. [Google Scholar] [CrossRef] [PubMed]
  43. Souza, P.F.; Xavier, D.R.; Suarez Mutis, M.C.; da Mota, J.C.; Peiter, P.C.; de Matos, V.P.; Magalhaes, M.; Barcellos, C. Spatial spread of malaria and economic frontier expansion in the Brazilian Amazon. PLoS ONE 2019, 14, e0217615. [Google Scholar] [CrossRef]
  44. Renia, L.; Goh, Y.S. Malaria Parasites: The Great Escape. Front. Immunol. 2016, 7, 463. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Rougeron, V.; Elguero, E.; Arnathau, C.; Acuna Hidalgo, B.; Durand, P.; Houze, S.; Berry, A.; Zakeri, S.; Haque, R.; Shafiul Alam, M.; et al. Human Plasmodium vivax diversity, population structure and evolutionary origin. PLoS Negl. Trop. Dis. 2020, 14, e0008072. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. de Oliveira, T.C.; Corder, R.M.; Early, A.; Rodrigues, P.T.; Ladeia-Andrade, S.; Alves, J.M.P.; Neafsey, D.E.; Ferreira, M.U. Population genomics reveals the expansion of highly inbred Plasmodium vivax lineages in the main malaria hotspot of Brazil. PLoS Negl. Trop. Dis. 2020, 14, e0008808. [Google Scholar] [CrossRef] [PubMed]
  47. Neafsey, D.E.; Galinsky, K.; Jiang, R.H.; Young, L.; Sykes, S.M.; Saif, S.; Gujja, S.; Goldberg, J.M.; Young, S.; Zeng, Q.; et al. The malaria parasite Plasmodium vivax exhibits greater genetic diversity than Plasmodium falciparum. Nat. Genet. 2012, 44, 1046–1050. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Molecular detection of the pvcyrpa gene by PCR. (A) Schematic representation of pvcyrpa gene, showing the fragments (exons 1 and 2) that were successfully amplified by two different primers sets (PvCyRPA F1/R1 and PvCyRPA F2/R2). (B) PCR amplification of the pvcyrpa gene. The figure shows agarose gel of pvcyrpa fragments amplified using PvCyRPA_F1/R1 primer and PvCyRPA_F2/R2 primer sets, respectively: Lane 1 and 9:100 bp molecular marker; Lane 2: negative control (water); Lane 3: P. falciparum in vitro culture (amplification with pvcyrpa primers); Lane 4: P. falciparum in vitro culture (amplification with p126 primers); Lane 5 and 7: PCR positive control (P. vivax infected sample); Lanes 6 and 8: P. vivax samples.
Figure 1. Molecular detection of the pvcyrpa gene by PCR. (A) Schematic representation of pvcyrpa gene, showing the fragments (exons 1 and 2) that were successfully amplified by two different primers sets (PvCyRPA F1/R1 and PvCyRPA F2/R2). (B) PCR amplification of the pvcyrpa gene. The figure shows agarose gel of pvcyrpa fragments amplified using PvCyRPA_F1/R1 primer and PvCyRPA_F2/R2 primer sets, respectively: Lane 1 and 9:100 bp molecular marker; Lane 2: negative control (water); Lane 3: P. falciparum in vitro culture (amplification with pvcyrpa primers); Lane 4: P. falciparum in vitro culture (amplification with p126 primers); Lane 5 and 7: PCR positive control (P. vivax infected sample); Lanes 6 and 8: P. vivax samples.
Genes 12 01657 g001
Figure 2. Frequency of substitutions in the PvCyRPA of P. vivax isolates from Brazilian Amazon. The two colors (black and gray) represent the two alleles within a population: reference allele and variant allele, respectively. The identification code represents the amino acid in the Sal-1 reference sequence (first character), followed by the position of this residue (number), and the replacing amino acid observed at the same position (last character).
Figure 2. Frequency of substitutions in the PvCyRPA of P. vivax isolates from Brazilian Amazon. The two colors (black and gray) represent the two alleles within a population: reference allele and variant allele, respectively. The identification code represents the amino acid in the Sal-1 reference sequence (first character), followed by the position of this residue (number), and the replacing amino acid observed at the same position (last character).
Genes 12 01657 g002
Figure 3. Median-joining network of pvcyrpa haplotypes. Each circle represents a unique haplotype and the color of the circles represents the geographic origins of each haplotype, while the size of the circle represents the frequency of each haplotype. Lines separating haplotypes represent mutational steps.
Figure 3. Median-joining network of pvcyrpa haplotypes. Each circle represents a unique haplotype and the color of the circles represents the geographic origins of each haplotype, while the size of the circle represents the frequency of each haplotype. Lines separating haplotypes represent mutational steps.
Genes 12 01657 g003
Figure 4. Homology analysis of PvCyRPA protein across isolates. (A) Circular alignment of amino acid sequences of PvCyRPA protein with reference sequence Sal-1 and P01 strain. The outer circle shows the amino acid scale. Green and gray bars on the second circle show the percent matching among all sequences used in the analysis. Inner circles show the sequence alignment in which each amino acid was represented by a different color. The legend presents the ordination to outer to inner rings: Scale; Percent Matching; Consensus OIA; Consensus MAO; Consensus GJ; Consensus ML; Consensus CZS; P01 strain; Sal-1 strain. (B) The pairwise distance among all sequences studied.
Figure 4. Homology analysis of PvCyRPA protein across isolates. (A) Circular alignment of amino acid sequences of PvCyRPA protein with reference sequence Sal-1 and P01 strain. The outer circle shows the amino acid scale. Green and gray bars on the second circle show the percent matching among all sequences used in the analysis. Inner circles show the sequence alignment in which each amino acid was represented by a different color. The legend presents the ordination to outer to inner rings: Scale; Percent Matching; Consensus OIA; Consensus MAO; Consensus GJ; Consensus ML; Consensus CZS; P01 strain; Sal-1 strain. (B) The pairwise distance among all sequences studied.
Genes 12 01657 g004
Table 1. Polymerase chain reaction (PCR) primers used for the amplification of the pvcyrpa gene.
Table 1. Polymerase chain reaction (PCR) primers used for the amplification of the pvcyrpa gene.
PCR PrimersDirectionSequence (5′–3′)Gene Length (bp)
Exon-1PvCyRPA_F1Forward primerTGCAATTTTCCTCTTCTTCTTCC600
PvCyRPA_R1Reverse primerCCCCATGTTCTTCCCTTGTCTT
Exon-2PvCyRPA_F2Forward primerCCTGATAAACTACGAAGAGCTCCA466
PvCyRPA_R2Reverse primerCCTCGTATAGTAAAGCGTGTGT
Specific primers used for conventional PCR and sequencing of the pvcyrpa gene.
Table 2. Polymorphisms in the PvCyRPA of P. vivax isolates from the Brazilian Amazon.
Table 2. Polymorphisms in the PvCyRPA of P. vivax isolates from the Brazilian Amazon.
Sal-1 aSubstitutions bIsolates
CZS (31) N (%)ML (17) N (%)GJ (4) N (%)MAO (37) N (%)OIA (9) N (%)Total (98) N (%)
G69G69A12 (39%)7 (41%)2 (50%)19 (51%)3 (33%)43 (44%)
E86E86Q4 (13%)1 (6%)002 (22%)7 (7%)
T90T90K14 (45%)10 (59%)2 (50%)19 (51%)2 (22%)47 (48%)
K93K93E14 (45%)9 (53%)2 (50%)19 (51%)3 (33%)47 (48%)
H95H95N14 (45%)9 (53%)2 (50%)19 (51%)3 (33%)47 (48%)
R122R122K24 (77%)15 (88%)4 (100%)30 (81%)7 (78%)80 (82%)
E126E126K14 (45%)7 (41%)2 (50%)17 (46%)3 (33%)43 (44%)
G127G127E14 (45%)7 (41%)2 (50%)17 (46%)2 (22%)42 (43%)
I129I129S9 (29%)6 (35%)2 (50%)11 (30%)3 (33%)31 (32%)
K131K131E23 (74%)13 (76%)4 (100%)30 (81%)7 (78%)77 (79%)
Q142Q142K17 (55%)8 (47%)2 (50%)15 (41%)6 (67%)48 (49%)
Q142R1 (3%)2 (12%)0003 (3%)
D145D145G9 (29%)6 (35%)2 (50%)11 (30%)3 (33%)31 (32%)
D145N8 (26%)2 (12%)06 (16%)3 (33%)19 (19%)
D149D149G23 (74%)13 (76%)4 (100%)16 (43%)6 (67%)62 (63%)
K150K150R6 (19%)5 (29%)2 (50%)0013 (13%)
A154A154D20 (65%)11 (65%)4 (100%)19 (51%)6 (67%)60 (61%)
E159E159D23 (74%)14 (82%)4 (100%)19 (51%)6 (67%)66 (67%)
N170N170D12 (39%)8 (47%)2 (50%)7 (19%)029 (30%)
K185K185E10 (32%)6 (35%)2 (50%)0 (0%)1 (11%)19 (19%)
A187A187V2 (6%)2 (12%)02 (5%)06 (6%)
V220V220I2 (6%)2 (12%)013 (35%)2 (22%)19 (19%)
D232D232N9 (29%)3 (18%)3 (75%)7 (19%)3 (33%)25 (26%)
D259D259E12 (39%)8 (47%)1 (25%)18 (49%)2 (22%)41 (42%)
V260V260E12 (39%)8 (47%)1 (25%)18 (49%)2 (22%)41 (42%)
T261T261P12 (39%)8 (47%)1 (25%)18 (49%)2 (22%)41 (42%)
G264G264S14 (45%)8 (47%)1 (25%)18 (49%)2 (22%)43 (44%)
E267***14 (45%)8 (47%)1 (25%)18 (49%)2 (22%)43 (44%)
V287V287I4 (13%)2 (12%)01 (3%)07 (7%)
a Reference Sequence Salvador-1; b Substitutions: The first letter represents the amino acid in the Sal-1 reference sequence and the last in bold the replacing amino acid; Study areas: CZS, Cruzeiro do Sul; ML, Mâncio Lima; GJ, Guajará; MAO, Manaus; OIA, Oiapoque. The number of samples per locality is indicates in parentheses; ***: Deletion of amino acid in position E267; N (%): Frequency of substitution in each locality.
Table 3. Diversity parameters and natural selection for the pvcyrpa gene in P. vivax isolates from the Brazilian Amazon.
Table 3. Diversity parameters and natural selection for the pvcyrpa gene in P. vivax isolates from the Brazilian Amazon.
DiversityEntire CodingExon-1Exon-2
Cruzeiro do Sul (n = 31)
No. of segregating sites (S)322210
No. of haplotypes (h)201111
Haplotype diversity (Hd)0.955 ± 0.0220.875 ± 0.0310.867 ± 0.038
Nucleotide diversity (π)0.01248 ± 0.000370.01489 ± 0.000500.00921 ± 0.00048
Tajima’s test (D)2.24953 *2.17555 *1.96672 ns1
Tajima’s test (Z)0.7741.3300.884
Mâncio Lima (n = 17)
No. of segregating sites (S)322210
No. of haplotypes (h)1185
Haplotype diversity (Hd)0.926 ± 0.0450.846 ± 0.0620.750 ± 0.088
Nucleotide diversity (π)0.01272 ± 0.000640.01517 ± 0.000950.00938 ± 0.00077
Tajima’s test (D)1.63472 ns1.59064 ns1.47204 ns
Tajima’s test (Z)1.0141.3910.743
Guajará (n = 4)
No. of segregating sites (S)20137
No. of haplotypes (h)322
Haplotype diversity (Hd)0.833 ± 0.2220.667 ± 0.2040.500 ± 0.265
Nucleotide diversity (π)0.01168 ± 0.003380.01444 ± 0.004420.00792 ± 0.00420
Tajima’s test (D)1.18178 ns2.24818 *−0.81734 ns
Tajima’s test (Z)0.5020.6860.440
Manaus (n = 37)
No. of segregating sites (S)27189
No. of haplotypes (h)19116
Haplotype diversity (Hd)0.943 ± 0.0210.886 ± 0.0310.818 ± 0.024
Nucleotide diversity (π)0.01116 ± 0.000300.01285 ± 0.000500.00887 ± 0.00034
Tajima’s test (D)2.75595 **2.60865 **2.42821 *
Tajima’s test (Z)0.9081.2580.685
Oiapoque (n = 9)
No. of segregating sites (S)281810
No. of haplotypes (h)976
Haplotype diversity (Hd)1.000 ± 0.0520.917 ± 0.0920.889 ± 0.091
Nucleotide diversity (π)0.01136 ± 0.001130.01361 ± 0.001710.00830 ± 0.00233
Tajima’s test (D)0.74327 ns1.13901 ns−0.01607 ns
Tajima’s test (Z)0.9321.1560.402
InDels were not included in the analysis; ns not significant (p > 0.10); ns1 not significant (0.10 > p > 0.05); * p < 0.05; ** p < 0.01.
Table 4. Comparison of amino acid variations of the PvCyRPA protein between the Brazilian Amazon and genome sequences available worldwide.
Table 4. Comparison of amino acid variations of the PvCyRPA protein between the Brazilian Amazon and genome sequences available worldwide.
IsolatesExon-1Exon-2
Codon Number/Amino Acid Residue
636986909395122125126127129131142145147149150154159170180185187220232 260261264287361
Sal-1IGETKHRREGIKQDQDKAENLKAVDDVTGVY
Brazilian AmazonCruzeiro do SulG/AE/QT/KK/EH/NR/KE/KG/EI/SK/EQ/R/KD/G/ND/GK/RA/DE/DN/DK/EA/VV/ID/ND/EV/ET/PG/SV/I
Mâncio LimaG/AE/QT/KK/EH/NR/KE/KG/EI/SK/EQ/R/KD/G/ND/GK/RA/DE/DN/DK/EA/VV/ID/ND/EV/ET/PG/SV/I
GuajaráG/AT/KK/EH/NR/KE/KG/EI/SK/EQ/KD/GD/GK/RA/DE/DN/DK/ED/ND/EV/ET/PG/S
ManausG/AT/KK/EH/NR/KE/KG/EI/SK/EQ/KD/G/ND/GA/DE/DN/DA/VID/ND/EV/ET/PG/SV/I
OiapoqueG/AE/QT/KK/EH/NR/KE/KG/EI/SK/EQ/KD/G/ND/GA/DE/DK/EV/ID/ND/EV/ET/PG/S
GenBankgb|KMZ81773.1 IndiaVIIKGHN
gb|KNA00954.1 North KoreanTKKKEERGDDVIH
gb|KMZ94334.1 Mauritania IK SEKGGNEEPS
gb|KMZ87926.1 Brazil IAKENKKEEKNGDDVN
Sanger InstituteSCO 66052.1ANNKTKEETKGDDDVII
SCO 71483.1 KKERGDVN
SGX 76259.1AQKENKKEEKDEPEEP
MexicoSouthern MexicanAENKENDDDEEPS
The amino acid variants of the PvCyRPA protein were compared to the reference Sal-1 reference sequence (PVX_090240). • Indicates identical amino acid residues compared to the Sal-1 strain. Codons from 69 to 187 and from 220 to 363 correspond to exon-1 and exon-2, respectively.
Table 5. B-cell epitope mapping of PvCyRPA (Salvador-1 Strain).
Table 5. B-cell epitope mapping of PvCyRPA (Salvador-1 Strain).
SequenceStartEndLenghtConformational EpitopeOverlapped Prediction
ElliproBCPredABCpredBepiPredEmini
PvCyRPA-B1819111YesXXXX-
PvCyRPA-B211912911NoX-XXX
PvCyRPA-B313415118NoXXXXX
PvCyRPA-B418119212YesX--XX
PvCyRPA-B52412499YesXX-X-
PvCyRPA-B625727216YesXX--X
PvCyRPA-B731233019YesXXX--
Full-length protein was analyzed for B-cell linear epitopes by the Ellipro algorithm, following confirmation by the overlap between the predictions of the algorithms BCPred, ABCpred, BepiPred, and Emini. (X): confirmed; (-): no confirmed.
Table 6. Predicted B cell epitopes of PvCyRPA and related polymorphisms.
Table 6. Predicted B cell epitopes of PvCyRPA and related polymorphisms.
EpitopeSequenceEpitope
PvCyRPA(I81–L91)Sal-1INSTWETQTTL
E86QINSTWQTQTTL
T90KINSTWETQTKL
T90NINSTWETQTNL
PvCyRPA(Y119–I129)Sal-1YKQRSKREGTI
R122KYKQKSKREGTI
R125TYKQRSKTEGTI
E126KYKQRSKRKGTI
G127EYKQRSKREGTI
I129SYKQRSKREGTS
PvCyRPA(N134–151)Sal-1NSVTGTIYQKEDVQIDKE
Q142RNSVTGTIYRKEDVQIDKE
Q142KNSVTGTIYKKEDVQIDKE
Q142TNSVTGTIYTKEDVQIDKE
D145GNSVTGTIYQKEDVQIDKE
D145NNSVTGTIYQKENVQIDKE
Q147KNSVTGTIYQKEDVKIDKE
D149GNSVTGTIYQKEDVQIGKE
K150RNSVTGTIYQKEDVQIDRE
PvCyRPA(S181–F192)Sal-1SYEYKTANKDNF
K185ESYEYETANKDNF
A187VSYEYKTVNKDNF
A187PSYEYKTPNKDNF
PvCyRPA(R241–R249)Sal-1RISTNNTAR
PvCyRPA(T257–C272)Sal-1TLDVTNEGKKEYKFKC
D259ETLEVTNEGKKEYKFKC
V260ETLDETNEGKKEYKFKC
T261PTLDVPNEGKKEYKFKC
G264STLDVTNESKKEYKFKC
PvCyRPA(T312–G330)Sal-1TEQNAIVVKPKVQNDDLNG
Sequences containing polymorphism inserted on predicted epitopes and marked in red letters.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bitencourt Chaves, L.; Guimarães, G.d.O.; Perce-da-Silva, D.d.S.; Banic, D.M.; Totino, P.R.R.; Machado, R.L.D.; Rodrigues-da-Silva, R.N.; Pratt-Riccio, L.R.; Daniel-Ribeiro, C.T.; Lima-Junior, J.d.C. Genetic Diversity of Plasmodium vivax Cysteine-Rich Protective Antigen (PvCyRPA) in Field Isolates from Five Different Areas of the Brazilian Amazon. Genes 2021, 12, 1657. https://doi.org/10.3390/genes12111657

AMA Style

Bitencourt Chaves L, Guimarães GdO, Perce-da-Silva DdS, Banic DM, Totino PRR, Machado RLD, Rodrigues-da-Silva RN, Pratt-Riccio LR, Daniel-Ribeiro CT, Lima-Junior JdC. Genetic Diversity of Plasmodium vivax Cysteine-Rich Protective Antigen (PvCyRPA) in Field Isolates from Five Different Areas of the Brazilian Amazon. Genes. 2021; 12(11):1657. https://doi.org/10.3390/genes12111657

Chicago/Turabian Style

Bitencourt Chaves, Lana, Glaucia de Oliveira Guimarães, Daiana de Souza Perce-da-Silva, Dalma Maria Banic, Paulo Renato Rivas Totino, Ricardo Luiz Dantas Machado, Rodrigo Nunes Rodrigues-da-Silva, Lilian Rose Pratt-Riccio, Cláudio Tadeu Daniel-Ribeiro, and Josué da Costa Lima-Junior. 2021. "Genetic Diversity of Plasmodium vivax Cysteine-Rich Protective Antigen (PvCyRPA) in Field Isolates from Five Different Areas of the Brazilian Amazon" Genes 12, no. 11: 1657. https://doi.org/10.3390/genes12111657

APA Style

Bitencourt Chaves, L., Guimarães, G. d. O., Perce-da-Silva, D. d. S., Banic, D. M., Totino, P. R. R., Machado, R. L. D., Rodrigues-da-Silva, R. N., Pratt-Riccio, L. R., Daniel-Ribeiro, C. T., & Lima-Junior, J. d. C. (2021). Genetic Diversity of Plasmodium vivax Cysteine-Rich Protective Antigen (PvCyRPA) in Field Isolates from Five Different Areas of the Brazilian Amazon. Genes, 12(11), 1657. https://doi.org/10.3390/genes12111657

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop