Exploring the genetic diversity pattern of PvEBP/DBP2: A promising candidate for an effective Plasmodium vivax vaccine.

Malaria remains a public health challenge. Since many control strategies have proven ineffective in eradicating this disease, new strategies are required, among which the design of a multivalent vaccine stands out. However, the effectiveness of this strategy has been hindered, among other reasons, by the genetic diversity observed in parasite antigens. In Plasmodium vivax, the Erythrocyte Binding Protein (PvEBP, also known as DBP2) is an alternate ligand to Duffy Binding Protein (DBP); given its structural resemblance to DBP, EBP/DBP2 is proposed as a promising antigen for inclusion in vaccine design. However, the extent of genetic diversity within the locus encoding this protein has not been comprehensively assessed. Thus, this study aimed to characterize the genetic diversity of the locus encoding the P. vivax EBP/DBP2 protein and to determine the evolutionary mechanisms modulating this diversity. Several intrapopulation genetic variation parameters were estimated using 36 gene sequences of PvEBP/DBP2 from Colombian P. vivax clinical isolates and 186 sequences available in databases. The study then evaluated the worldwide genetic structure and the evolutionary forces that may influence the observed patterns of genetic variation. It was found that the PvEBP/DBP2 gene exhibits one of the lowest levels of genetic diversity compared to other vaccine-candidate antigens. Four major haplotypes were shared worldwide. Analysis of the protein's 3D structure and epitope prediction identified five regions with potential antigenic properties. The results suggest that the PvEBP/DBP2 protein possesses ideal characteristics to be considered when designing a multivalent effective antimalarial vaccine against P. vivax.

Evaluating the genetic diversity of the Plasmodium vivax siap2 locus: A promising candidate for an effective malaria vaccine?

Malaria is the deadliest parasitic disease in the world. Traditional control measures have become less effective; hence, there is a need to explore alternative strategies, such as antimalarial vaccines. However, designing an anti-Plasmodium vivax vaccine is considered a challenge due to the complex parasite biology and the antigens' high genetic diversity. Recently, the sporozoite invasion-associated protein 2 (SIAP2) has been suggested as a potential antigen to be considered in vaccine design due to its significance during hepatocyte invasion. However, its use may be limited by the incomplete understanding of gene/protein diversity. Here, the genetic diversity of pvsiap2 using P. vivax DNA samples from Colombia was assessed. Through PCR amplification and sequencing, we compared the Colombian sequences with available worldwide sequences, revealing that pvsiap2 displays low genetic diversity. Molecular evolutionary analyses showed that pvsiap2 appears to be influenced by directional selection. Moreover, the haplotypes found differ by a few mutational steps and several of them were shared between different geographical areas. On the other hand, several conserved regions within PvSIAP2 were predicted as potential B-cell or T-cell epitopes. Considering these characteristics and its role in hepatocyte invasion, the PvSIAP2 protein emerges as a promising antigen to be considered in a multi-antigen-multi-stage (multivalent) fully effective vaccine against P. vivax malaria.

Immunoglobulin heavy constant gamma gene evolution is modulated by both the divergent and birth-and-death evolutionary models.

Immunoglobulin G (IgG) is one of the five antibody classes produced in mammals as part of the humoral responses accountable for protecting the organisms from infection. Its antibody heavy chain constant region is encoded by the Ig heavy-chain gamma gene (IGHG). In humans, there are four IGHG genes which encode the four subclasses, each with a specialized effector function. Although four subclasses of IgG proteins have also been reported in macaques, this does not appear to be the rule for all primates. In Platyrrhini, IgG has been stated to be encoded by a single-copy gene. To date, it remains unknown how the IGHG has expanded or contracted in the primate order; consequently, we have analyzed data from 38 primate genome sequences to identify IGHG genes and describe the evolution of IGHG genes in primate order. IGHG belongs to a multigene family that evolves by the birth-death evolutionary model in primates. Whereas Strepsirrhini and Platyrrhini have a single-copy gene, in Catarrhini, it has expanded to several paralogs in their genomes; some deleted and others pseudogenized. Furthermore, episodic positive selection may have promoted a species-specific IgG effector function. We propose that IgG evolved to reach an optimal number of copies per genome to adapt their humoral immune responses to different environmental conditions. This study has implications for biomedical trials using non-human primates.

IgG subclasses in New World Monkeys: an issue for debate?

Immunoglobulin G (IgG) is an essential antibody in adaptive immunity; a differential expansion of the gene encoding the Fc region (IGHG) of this antibody has been observed in mammals. Like humans, animal biomedical models, such as mice and macaques, have four functional genes encoding 4 IgG subclasses; however, the data for New World monkeys (NWM) seems contentious. Some publications argue for the existence of a single-copy gene for IgG Fc; however, a recent paper has suggested the presence of IgG subclasses in some NWM species. Here, we evaluated the genetic distances and phylogenetic relationships in NWM to assess the presence of IgG subclasses using the sequences of IGHG genes from 13 NWM species recovered from genomic data and lab PCR and cloning-based procedures available in GenBank. The results show that several sequences do not cluster into the expected taxon, probably due to cross-contamination during laboratory procedures, and consequently, they appear to be wrongly assigned. Additionally, several sequences reported as subclasses were shown to be 100% identical in the CH domains. The data presented here suggests that there is not enough evidence to establish the presence of IgG subclasses in NWM.

Genetic diversity of SARS-CoV-2 in South America: demographic history and structuration signals.

In 2020, the emergence of SARS-CoV-2 caused a global public health crisis with significant mortality rates and a large socioeconomic burden. The rapid spread of this new virus has led to the appearance of new variants, making the characterization and monitoring of genetic diversity necessary to understand the population dynamics and evolution of the virus. Here, a population-genetics-based study was performed starting with South American genome sequences available in the GISAID database to investigate the genetic diversity of SARS-CoV-2 on this continent and the evolutionary mechanisms that modulate it.

Two 20-Residue-Long Peptides Derived from Plasmodium vivax Merozoite Surface Protein 10 EGF-Like Domains Are Involved in Binding to Human Reticulocytes.

Plasmodium parasites' invasion of their target cells is a complex, multi-step process involving many protein-protein interactions. Little is known about how complex the interaction with target cells is in Plasmodium vivax and few surface molecules related to reticulocytes' adhesion have been described to date. Natural selection, functional and structural analysis were carried out on the previously described vaccine candidate P. vivax merozoite surface protein 10 (PvMSP10) for evaluating its role during initial contact with target cells. It has been shown here that the recombinant carboxyl terminal region (rPvMSP10-C) bound to adult human reticulocytes but not to normocytes, as validated by two different protein-cell interaction assays. Particularly interesting was the fact that two 20-residue-long regions (388DKEECRCRANYMPDDSVDYF407 and 415KDCSKENGNCDVNAECSIDK434) were able to inhibit rPvMSP10-C binding to reticulocytes and rosette formation using enriched target cells. These peptides were derived from PvMSP10 epidermal growth factor (EGF)-like domains (precisely, from a well-defined electrostatic zone) and consisted of regions having the potential of being B- or T-cell epitopes. These findings provide evidence, for the first time, about the fragments governing PvMSP10 binding to its target cells, thus highlighting the importance of studying them for inclusion in a P. vivax antimalarial vaccine.

Plasmodium vivax Cell Traversal Protein for Ookinetes and Sporozoites (CelTOS) Functionally Restricted Regions Are Involved in Specific Host-Pathogen Interactions.

Following the injection of Plasmodium sporozoites by a female Anopheles mosquito into the dermis, they become engaged on a long journey to hepatic tissue where they must migrate through different types of cell to become established in parasitophorous vacuoles in hepatocytes. Studies have shown that proteins such as cell traversal protein for Plasmodium ookinetes and sporozoites (CelTOS) play a crucial role in cell-traversal ability. Although CelTOS has been extensively studied in various species and included in pre-clinical assays it remains unknown which P. vivax CelTOS (PvCelTOS) regions are key in its interaction with traversed or target cells (Kupffer or hepatocytes) and what type of pressure, association and polymorphism these important regions could have to improve their candidacy as important vaccine antigens. This work has described producing a recombinant PvCelTOS which was recognized by ~30% P. vivax-infected individuals, thereby confirming its ability for inducing a natural immune response. PvCelTOS' genetic diversity in Colombia and its ability to interact with HeLa (traversal cell) and/or HepG2 cell (target cell) external membrane have been assessed. One region in the PvCelTOS amino-terminal region and another in its C-terminus were seen to be participating in host-pathogen interactions. These regions had important functional constraint signals (ω < 0.3 and several sites under negative selection) and were able to inhibit specific rPvCelTOS binding to HeLa cells. This led to suggesting that sequences between aa 41-60 (40833) and 141-160 (40838) represent promising candidates for an anti-P. vivax subunit-based vaccine.

Igh locus structure and evolution in Platyrrhines: new insights from a genomic perspective.

Non-human primates have been used as animal models because of their phylogenetic closeness to humans. However, the genetic differences between humans and non-human primates must be considered to select the appropriate animal models. Recently, New World monkeys (Platyrrhines) have generated a higher interest in biomedical research, especially in assessing vaccine safety and immunogenicity. Given the continued and renewed interest in Platyrrhines as biomedical models, it is a necessary to have a better and more complete understanding of their immune system and its implications for research. Immunoglobulins (Ig) are the main proteins that mediate humoral immunity. These proteins have evolved as part of an adaptive immune response system derived from ancient vertebrates. There are at least four Ig classes in Prosimians, whereas five have been reported in Catarrhines. Information on the structure and evolution of the loci containing immunoglobulin heavy chain constant genes (Igh) in Platyrrhines, however, is limited. Here, Igh loci were characterized in 10 Platyrrhines using the available whole genome sequences. Human and Macaca Igh loci were also assessed to compare them with their Platyrrhines counterparts. Differences in Igh locus structure were observed between Platyrrhines and Catarrhines. Noteworthy changes occur in the γ gene, which encodes a key Ig involved in organism defense that would favor protection after vaccination. The remarkable differences between the immunoglobulin proteins of Platyrrhines and Catarrhines warrant a cautionary message to biomedical researchers.

Using next-generation sequencing for characterising HLA-DRB1 and DQB1 loci in a cohort of Colombian women.

The Colombian population is characterised by a high genetic diversity, secondary to the ethnic mixture arising from colonisation. Unfortunately, few reports are available regarding HLA-DRB1 and DQB1 diversity in Colombia to date. HLA-DRB1 and DQB1 diversity was identified in this study using next-generating sequencing (NGS) on a cohort of Colombian women. Cervical samples taken from 276 women were used for typing DRB1 and DQB1 loci by Illumina MiSeq. Allele and haplotype frequencies were calculated using an expectation-maximisation algorithm. Hardy-Weinberg Equilibrium and linkage disequilibrium (LD) between loci were evaluated. Forty-seven DRB1 alleles and 14 DQB1 alleles were identified. DRB1*04:07:01G and DQB1*03:02:01G alleles occurred most frequently in the target population. Significant LD was found in 44 out of the 144 identified haplotypes, within which DRB1*04:07:01G-DQB1*03:02:01G occurred most frequently (6.56%). The alleles and haplotypes found with NGS agreed with that found in previous reports involving lower resolution for the Colombian population, and greater genetic variability was found, especially concerning DRB1. Comparing allele and haplotype frequency distribution in the target population to that of other populations denoted HLA system intra- and inter-population diversity.

Behavior and abundance of Anopheles darlingi in communities living in the Colombian Amazon riverside.

In the past few years, relative frequencies of malaria parasite species in communities living in the Colombian Amazon riverside have changed, being Plasmodium vivax (61.4%) and Plasmodium malariae (43.8%) the most frequent. Given this epidemiological scenario, it is important to determine the species of anophelines involved in these parasites' transmission. This study was carried out in June 2016 in two indigenous communities living close to the tributaries of the Amazon River using protected human bait. The results of this study showed a total abundance of 1,085 mosquitos, of which 99.2% corresponded to Anopheles darlingi. Additionally, only two anopheline species were found, showing low diversity in the study areas. Molecular confirmation of some individuals was then followed by evolutionary analysis by using the COI gene. Nested PCR was used for identifying the three Plasmodium species circulating in the study areas. Of the two species collected in this study, 21.0% of the An. darlingi mosquitoes were infected with P. malariae, 21.9% with P. vivax and 10.3% with Plasmodium falciparum. It exhibited exophilic and exophagic behavior in both study areas, having marked differences regarding its abundance in each community (Tipisca first sampling 49.4%, Tipisca second sampling 39.6% and Doce de Octubre 10.9%). Interestingly, An. mattogrossensis infected by P. vivax was found for the first time in Colombia (in 50% of the four females collected). Analysis of An. darlingi COI gene diversity indicated a single population maintaining a high gene flow between the study areas. The An. darlingi behavior pattern found in both communities represents a risk factor for the region's inhabitants living/working near these sites. This highlights the need for vector control efforts such as the use of personal repellents and insecticides for use on cattle, which must be made available in order to reduce this Anopheline's abundance.

On the Evolution and Function of Plasmodium vivax Reticulocyte Binding Surface Antigen (pvrbsa).

The RBSA protein is encoded by a gene described in Plasmodium species having tropism for reticulocytes. Since this protein is antigenic in natural infections and can bind to target cells, it has been proposed as a potential candidate for an anti-Plasmodium vivax vaccine. However, genetic diversity (a challenge which must be overcome for ensuring fully effective vaccine design) has not been described at this locus. Likewise, the minimum regions mediating specific parasite-host interaction have not been determined. This is why the rbsa gene's evolutionary history is being here described, as well as the P. vivax rbsa (pvrbsa) genetic diversity and the specific regions mediating parasite adhesion to reticulocytes. Unlike what has previously been reported, rbsa was also present in several parasite species belonging to the monkey-malaria clade; paralogs were also found in Plasmodium parasites invading reticulocytes. The pvrbsa locus had less diversity than other merozoite surface proteins where natural selection and recombination were the main evolutionary forces involved in causing the observed polymorphism. The N-terminal end (PvRBSA-A) was conserved and under functional constraint; consequently, it was expressed as recombinant protein for binding assays. This protein fragment bound to reticulocytes whilst the C-terminus, included in recombinant PvRBSA-B (which was not under functional constraint), did not. Interestingly, two PvRBSA-A-derived peptides were able to inhibit protein binding to reticulocytes. Specific conserved and functionally important peptides within PvRBSA-A could thus be considered when designing a fully-effective vaccine against P. vivax.

Simultaneous detection of Plasmodium vivax dhfr, dhps, mdr1 and crt-o resistance-associated mutations in the Colombian Amazonian region.

BACKGROUND: Malaria continues being a public health problem worldwide. Plasmodium vivax is the species causing the largest number of cases of malaria in Asia and South America. Due to the lack of a completely effective anti-malarial vaccine, controlling this disease has been based on transmission vector management, rapid diagnosis and suitable treatment. However, parasite resistance to anti-malarial drugs has become a major yet-to-be-overcome challenge. This study was thus aimed at determining pvmdr1, pvdhfr, pvdhps and pvcrt-o gene mutations and haplotypes from field samples obtained from an endemic area in the Colombian Amazonian region. METHODS: Fifty samples of parasite DNA infected by a single P. vivax strain from symptomatic patients from the Amazonas department in Colombia were analysed by PCR and the pvdhfr, pvdhps, pvmdr1 and pvcrt-o genes were sequenced. Diversity estimators were calculated from the sequences and the haplotypes circulating in the Colombian Amazonian region were obtained. CONCLUSION: pvdhfr, pvdhps, pvmdr1 and pvcrt-o genes in the Colombian Amazonian region are characterized by low genetic diversity. Some resistance-associated mutations were found circulating in this population. New variants are also being reported. A selective sweep signal was located in pvdhfr and pvmdr1 genes, suggesting that these mutations (or some of them) could be providing an adaptive advantage.

Identifying Potential Plasmodium vivax Sporozoite Stage Vaccine Candidates: An Analysis of Genetic Diversity and Natural Selection.

Parasite antigen genetic diversity represents a great obstacle when designing a vaccine against malaria caused by Plasmodium vivax. Selecting vaccine candidate antigens has been focused on those fulfilling a role in invasion and which are conserved, thus avoiding specific-allele immune responses. Most antigens described to date belong to the blood stage, thereby blocking parasite development within red blood cells, whilst studying antigens from other stages has been quite restricted. Antigens from different parasite stages are required for developing a completely effective vaccine; thus, pre-erythrocyte stage antigens able to block the first line of infection becoming established should also be taken into account. However, few antigens from this stage have been studied to date. Several P. falciparum sporozoite antigens are involved in invasion. Since 77% of genes are orthologous amongst Plasmodium parasites, P. vivax sporozoite antigen orthologs to those of P. falciparum might be present in its genome. Although these genes might have high genetic diversity, conserved functionally-relevant regions (ideal for vaccine development) could be predicted by comparing genetic diversity patterns and evolutionary rates. This study was thus aimed at searching for putative P. vivax sporozoite genes so as to analyse their genetic diversity for determining their potential as vaccine candidates. Several DNA sequence polymorphism estimators were computed at each locus. The evolutionary force (drift, selection and recombination) drawing the genetic diversity pattern observed was also determined by using tests based on polymorphism frequency spectrum as well as the type of intra- and inter-species substitutions. Likewise, recombination was assessed both indirectly and directly. The results showed that sporozoite genes were more conserved than merozoite genes evaluated to date. Putative domains implied in cell traversal, gliding motility and hepatocyte interaction had a negative selection signal, being conserved amongst different species in the genus. PvP52, PvP36, PvSPATR, PvPLP1, PvMCP1, PvTLP, PvCelTOS, and PvMB2 antigens or functionally restricted regions within them would thus seem promising vaccine candidates and could be used when designing a pre-erythrocyte and/or multi-stage vaccine against P. vivax to avoid allele-specific immune responses that could reduce vaccine efficacy.

PvGAMA reticulocyte binding activity: predicting conserved functional regions by natural selection analysis.

BACKGROUND: Adhesin proteins are used by Plasmodium parasites to bind and invade target cells. Hence, characterising molecules that participate in reticulocyte interaction is key to understanding the molecular basis of Plasmodium vivax invasion. This study focused on predicting functionally restricted regions of the P. vivax GPI-anchored micronemal antigen (PvGAMA) and characterising their reticulocyte binding activity. RESULTS: The pvgama gene was initially found in P. vivax VCG-I strain schizonts. According to the genetic diversity analysis, PvGAMA displayed a size polymorphism very common for antigenic P. vivax proteins. Two regions along the antigen sequence were highly conserved among species, having a negative natural selection signal. Interestingly, these regions revealed a functional role regarding preferential target cell adhesion. CONCLUSIONS: To our knowledge, this study describes PvGAMA reticulocyte binding properties for the first time. Conserved functional regions were predicted according to natural selection analysis and their binding ability was confirmed. These findings support the notion that PvGAMA may have an important role in P. vivax merozoite adhesion to its target cells.

Evidence of functional divergence in MSP7 paralogous proteins: a molecular-evolutionary and phylogenetic analysis.

BACKGROUND: The merozoite surface protein 7 (MSP7) is a Plasmodium protein which is involved in parasite invasion; the gene encoding it belongs to a multigene family. It has been proposed that MSP7 paralogues seem to be functionally redundant; however, recent experiments have suggested that they could have different roles. RESULTS: The msp7 multigene family has been described in newly available Plasmodium genomes; phylogenetic relationships were established in 12 species by using different molecular evolutionary approaches for assessing functional divergence amongst MSP7 members. Gene expansion and contraction rule msp7 family evolution; however, some members could have had concerted evolution. Molecular evolutionary analysis showed that relaxed and/or intensified selection modulated Plasmodium msp7 paralogous evolution. Furthermore, episodic diversifying selection and changes in evolutionary rates suggested that some paralogous proteins have diverged functionally. CONCLUSIONS: Even though msp7 has mainly evolved in line with a birth-and-death evolutionary model, gene conversion has taken place between some paralogous genes allowing them to maintain their functional redundancy. On the other hand, the evolutionary rate of some MSP7 paralogs has become altered, as well as undergoing relaxed or intensified (positive) selection, suggesting functional divergence. This could mean that some MSP7s can form different parasite protein complexes and/or recognise different host receptors during parasite invasion. These results highlight the importance of this gene family in the Plasmodium genus.

Size polymorphism and low sequence diversity in the locus encoding the Plasmodium vivax rhoptry neck protein 4 (PvRON4) in Colombian isolates.

BACKGROUND: Designing a vaccine against Plasmodium vivax has focused on selecting antigens involved in invasion mechanisms that must have domains with low polymorphism for avoiding allele-specific immune responses. The rhoptry neck protein 4 (RON4) forms part of the tight junction, which is essential in the invasion of hepatocytes and/or erythrocytes; however, little is known about this locus' genetic diversity. METHODS: DNA sequences from 73 Colombian clinical isolates from pvron4 gene were analysed for characterizing their genetic diversity; pvron4 haplotype number and distribution, as well as the evolutionary forces determining diversity pattern, were assessed by population genetics and molecular evolutionary approaches. RESULTS: ron4 has low genetic diversity in P. vivax at sequence level; however, a variable amount of tandem repeats at the N-terminal region leads to extensive size polymorphism. This region seems to be exposed to the immune system. The central region has a putative esterase/lipase domain which, like the protein's C-terminal fragment, is highly conserved at intra- and inter-species level. Both regions are under purifying selection. CONCLUSIONS: pvron4 is the locus having the lowest genetic diversity described to date for P. vivax. The repeat regions in the N-terminal region could be associated with immune evasion mechanisms while the central region and the C-terminal region seem to be under functional or structural constraint. Bearing such results in mind, the PvRON4 central and/or C-terminal portions represent promising candidates when designing a subunit-based vaccine as they are aimed at avoiding an allele-specific immune response, which might limit vaccine efficacy.

Inferring natural selection signals in Plasmodium vivax-encoded proteins having a potential role in merozoite invasion.

Detecting natural selection signals in Plasmodium parasites antigens might be used for identifying potential new vaccine candidates. Fifty-nine Plasmodium vivax-Sal-I genes encoding proteins having a potential role in invasion were used as query for identifying them in recent P. vivax strain genome sequences and two closely-related Plasmodium species. Several measures of DNA sequence variation were then calculated and selection signatures were detected by using different approaches. Our results may be used for determining which genes expressed during P. vivax merozoite stage could be prioritised for further population genetics or functional studies for designing a P. vivax vaccine which would avoid allele-specific immune responses.

Heterogeneous genetic diversity pattern in Plasmodium vivax genes encoding merozoite surface proteins (MSP) -7E, -7F and -7L.

BACKGROUND: The msp-7 gene has become differentially expanded in the Plasmodium genus; Plasmodium vivax has the highest copy number of this gene, several of which encode antigenic proteins in merozoites. METHODS: DNA sequences from thirty-six Colombian clinical isolates from P. vivax (pv) msp-7E, -7F and -7L genes were analysed for characterizing and studying the genetic diversity of these pvmsp-7 members which are expressed during the intra-erythrocyte stage; natural selection signals producing the variation pattern so observed were evaluated. RESULTS: The pvmsp-7E gene was highly polymorphic compared to pvmsp-7F and pvmsp-7L which were seen to have limited genetic diversity; pvmsp-7E polymorphism was seen to have been maintained by different types of positive selection. Even though these copies seemed to be species-specific duplications, a search in the Plasmodium cynomolgi genome (P. vivax sister taxon) showed that both species shared the whole msp-7 repertoire. This led to exploring the long-term effect of natural selection by comparing the orthologous sequences which led to finding signatures for lineage-specific positive selection. CONCLUSIONS: The results confirmed that the P. vivax msp-7 family has a heterogeneous genetic diversity pattern; some members are highly conserved whilst others are highly diverse. The results suggested that the 3'-end of these genes encode MSP-7 proteins' functional region whilst the central region of pvmsp-7E has evolved rapidly. The lineage-specific positive selection signals found suggested that mutations occurring in msp-7s genes during host switch may have succeeded in adapting the ancestral P. vivax parasite population to humans.

Low genetic diversity in the locus encoding the Plasmodium vivax P41 protein in Colombia's parasite population.

BACKGROUND: The development of malaria vaccine has been hindered by the allele-specific responses produced by some parasite antigens' high genetic diversity. Such antigen genetic diversity must thus be evaluated when designing a completely effective vaccine. Plasmodium falciparum P12, P38 and P41 proteins have red blood cell binding regions in the s48/45 domains and are located on merozoite surface, P41 forming a heteroduplex with P12. These three genes have been identified in Plasmodium vivax and share similar characteristics with their orthologues in Plasmodium falciparum. Plasmodium vivax pv12 and pv38 have low genetic diversity but pv41 polymorphism has not been described. METHODS: The present study was aimed at evaluating the P. vivax p41 (pv41) gene's polymorphism. DNA sequences from Colombian clinical isolates from pv41 gene were analysed for characterising and studying the genetic diversity and the evolutionary forces that produced the variation pattern so observed. RESULTS: Similarly to other members of the 6-Cys family, pv41 had low genetic polymorphism. pv41 3'-end displayed the highest nucleotide diversity value; several substitutions found there were under positive selection. Negatively selected codons at inter-species level were identified in the s48/45 domains; p41 would thus seem to have functional/structural constraints due to the presence of these domains. CONCLUSIONS: In spite of the functional constraints of Pv41 s48/45 domains, immune system pressure seems to have allowed non-synonymous substitutions to become fixed within them as an adaptation mechanism; including Pv41 s48/45 domains in a vaccine should thus be carefully evaluated due to these domains containing some allele variants.

Low genetic diversity and functional constraint in loci encoding Plasmodium vivax P12 and P38 proteins in the Colombian population.

BACKGROUND: Plasmodium vivax is one of the five species causing malaria in human beings, affecting around 391 million people annually. The development of an anti-malarial vaccine has been proposed as an alternative for controlling this disease. However, its development has been hampered by allele-specific responses produced by the high genetic diversity shown by some parasite antigens. Evaluating these antigens' genetic diversity is thus essential when designing a completely effective vaccine. METHODS: The gene sequences of Plasmodium vivax p12 (pv12) and p38 (pv38), obtained from field isolates in Colombia, were used for evaluating haplotype polymorphism and distribution by population genetics analysis. The evolutionary forces generating the variation pattern so observed were also determined. RESULTS: Both pv12 and pv38 were shown to have low genetic diversity. The neutral model for pv12 could not be discarded, whilst polymorphism in pv38 was maintained by balanced selection restricted to the gene's 5' region. Both encoded proteins seemed to have functional/structural constraints due to the presence of s48/45 domains, which were seen to be highly conserved. CONCLUSIONS: Due to the role that malaria parasite P12 and P38 proteins seem to play during invasion in Plasmodium species, added to the Pv12 and Pv38 antigenic characteristics and the low genetic diversity observed, these proteins might be good candidates to be evaluated in the design of a multistage/multi-antigen vaccine.