The malaria parasite Plasmodium vivax exhibits greater genetic diversity than Plasmodium falciparum

Multiple Origins of Mutations in the mdr1 Gene--A Putative Marker of Chloroquine Resistance in P. vivax

Background

Chloroquine combined with primaquine has been the recommended antimalarial treatment of Plasmodium vivax malaria infections for six decades but the efficacy of this treatment regimen is threatened by chloroquine resistance (CQR). Single nucleotide polymorphisms (SNPs) in the multidrug resistance gene, Pvmdr1 are putative determinants of CQR but the extent of their emergence at population level remains to be explored.

Objective

In this study we describe the prevalence of SNPs in the Pvmdr1 among samples collected in seven P. vivax endemic countries and we looked for molecular evidence of drug selection by characterising polymorphism at microsatellite (MS) loci flanking the Pvmdr1 gene.

Methods

We examined the prevalence of SNPs in the Pvmdr1 gene among 267 samples collected from Pakistan, Afghanistan, Sri Lanka, Nepal, Sudan, São Tomé and Ecuador. We measured and diversity in four microsatellite (MS) markers flanking the Pvmdr1 gene to look evidence of selection on mutant alleles.

Results

SNP polymorphism in the Pvmdr1 gene was largely confined to codons T958M, Y976F and F1076L. Only 2.4% of samples were wildtype at all three codons (TYF, n = 5), 13.3% (n = 28) of the samples were single mutant MYF, 63.0% of samples (n = 133) were double mutant MYL, and 21.3% (n = 45) were triple mutant MFL. Clear geographic differences in the prevalence of these Pvmdr mutation combinations were observed.

Significant linkage disequilibrium (LD) between Pvmdr1 and MS alleles was found in populations sampled in Ecuador, Nepal and Sri Lanka, while significant LD between Pvmdr1 and the combined 4 MS locus haplotype was only seen in Ecuador and Sri Lanka. When combining the 5 loci, high level diversity, measured as expected heterozygosity (H_e), was seen in the complete sample set (H_e = 0.99), while H_e estimates for individual loci ranged from 0.00–0.93. Although Pvmdr1 haplotypes were not consistently associated with specific flanking MS alleles, there was significant differentiation between geographic sites which could indicate directional selection through local drug pressure.

Conclusions

Our observations suggest that Pvmdr1 mutations emerged independently on multiple occasions even within the same population. In Sri Lanka population analysis at multiple sites showed evidence of local selection and geographical dispersal of Pvmdr1 mutations between sites.

Chloroquine combined with primaquine has been the recommended antimalarial treatment for Plasmodium vivax malaria infections for sixty years but the efficacy of this treatment regimen is threatened by chloroquine resistance. In this study we describe the prevalence of mutations in the P. vivax gene, Pvmdr1 among samples collected in seven endemic countries. The mutations are thought to be associated with chloroquine resistance and here we looked for evidence of drug selection by characterising polymorphism in DNA repeat regions (microsatellite (MS) loci) flanking the Pvmdr1 gene. Mutations in the Pvmdr1 gene were mainly identified at codons T958M, Y976F and F1076L. Just 2.4% of samples were wildtype at all three codons, while 63% were single mutants (MYF). Clear geographic differences in the prevalence of these Pvmdr mutation combinations were observed. At the flanking MS loci, we found high levels of diversity, and significant differentiation between geographic sites. This pattern of variation could indicate directional selection through local drug pressure. In summary, our observations suggest that Pvmdr1 mutations and thus, chloroquine resistance has emerged independently on multiple occasions even within the same population.

African origin of the malaria parasite Plasmodium vivax

Plasmodium vivax is the leading cause of human malaria in Asia and Latin America but is absent from most of central Africa due to the near fixation of a mutation that inhibits the expression of its receptor, the Duffy antigen, on human erythrocytes. The emergence of this protective allele is not understood because P. vivax is believed to have originated in Asia. Here we show, using a non-invasive approach, that wild chimpanzees and gorillas throughout central Africa are endemically infected with parasites that are closely related to human P. vivax. Sequence analyses reveal that ape parasites lack host specificity and are much more diverse than human parasites, which form a monophyletic lineage within the ape parasite radiation. These findings indicate that human P. vivax is of African origin and likely selected for the Duffy-negative mutation. All extant human P. vivax parasites are derived from a single ancestor that escaped out of Africa.

Population genomic structure and adaptation in the zoonotic malaria parasite Plasmodium knowlesi

Malaria cases caused by the zoonotic parasite Plasmodium knowlesi are being increasingly reported throughout Southeast Asia and in travelers returning from the region. To test for evidence of signatures of selection or unusual population structure in this parasite, we surveyed genome sequence diversity in 48 clinical isolates recently sampled from Malaysian Borneo and in five lines maintained in laboratory rhesus macaques after isolation in the 1960s from Peninsular Malaysia and the Philippines. Overall genomewide nucleotide diversity (π = 6.03 × 10(-3)) was much higher than has been seen in worldwide samples of either of the major endemic malaria parasite species Plasmodium falciparum and Plasmodium vivax. A remarkable substructure is revealed within P. knowlesi, consisting of two major sympatric clusters of the clinical isolates and a third cluster comprising the laboratory isolates. There was deep differentiation between the two clusters of clinical isolates [mean genomewide fixation index (FST) = 0.21, with 9,293 SNPs having fixed differences of FST = 1.0]. This differentiation showed marked heterogeneity across the genome, with mean FST values of different chromosomes ranging from 0.08 to 0.34 and with further significant variation across regions within several chromosomes. Analysis of the largest cluster (cluster 1, 38 isolates) indicated long-term population growth, with negatively skewed allele frequency distributions (genomewide average Tajima's D = -1.35). Against this background there was evidence of balancing selection on particular genes, including the circumsporozoite protein (csp) gene, which had the top Tajima's D value (1.57), and scans of haplotype homozygosity implicate several genomic regions as being under recent positive selection.

The genomes of three stocks comprising the most widely utilized live sporozoite Theileria parva vaccine exhibit very different degrees and patterns of sequence divergence

BACKGROUND

There are no commercially available vaccines against human protozoan parasitic diseases, despite the success of vaccination-induced long-term protection against infectious diseases. East Coast fever, caused by the protist Theileria parva, kills one million cattle each year in sub-Saharan Africa, and contributes significantly to hunger and poverty in the region. A highly effective, live, multi-isolate vaccine against T. parva exists, but its component isolates have not been characterized. Here we sequence and compare the three component T. parva stocks within this vaccine, the Muguga Cocktail, namely Muguga, Kiambu5 and Serengeti-transformed, aiming to identify genomic features that contribute to vaccine efficacy.

RESULTS

We find that Serengeti-transformed, originally isolated from the wildlife carrier, the African Cape buffalo, is remarkably and unexpectedly similar to the Muguga isolate. The 420 detectable non-synonymous SNPs were distributed among only 53 genes, primarily subtelomeric antigens and antigenic families. The Kiambu5 isolate is considerably more divergent, with close to 40,000 SNPs relative to Muguga, including >8,500 non-synonymous mutations distributed among >1,700 (42.5 %) of the predicted genes. These genetic markers of the component stocks can be used to characterize the composition of new batches of the Muguga Cocktail.

CONCLUSIONS

Differences among these three isolates, while extensive, represent only a small proportion of the genetic variation in the entire species. Given the efficacy of the Muguga Cocktail in inducing long-lasting protection against infections in the field, our results suggest that whole-organism vaccines against parasitic diseases can be highly efficacious despite considerable genome-wide differences relative to the isolates against which they protect.

Identification of Immunodominant B-cell Epitope Regions of Reticulocyte Binding Proteins in Plasmodium vivax by Protein Microarray Based Immunoscreening

Plasmodium falciparum can invade all stages of red blood cells, while Plasmodium vivax can invade only reticulocytes. Although many P. vivax proteins have been discovered, their functions are largely unknown. Among them, P. vivax reticulocyte binding proteins (PvRBP1 and PvRBP2) recognize and bind to reticulocytes. Both proteins possess a C-terminal hydrophobic transmembrane domain, which drives adhesion to reticulocytes. PvRBP1 and PvRBP2 are large (> 326 kDa), which hinders identification of the functional domains. In this study, the complete genome information of the P. vivax RBP family was thoroughly analyzed using a prediction server with bioinformatics data to predict B-cell epitope domains. Eleven pvrbp family genes that included 2 pseudogenes and 9 full or partial length genes were selected and used to express recombinant proteins in a wheat germ cell-free system. The expressed proteins were used to evaluate the humoral immune response with vivax malaria patients and healthy individual serum samples by protein microarray. The recombinant fragments of 9 PvRBP proteins were successfully expressed; the soluble proteins ranged in molecular weight from 16 to 34 kDa. Evaluation of the humoral immune response to each recombinant PvRBP protein indicated a high antigenicity, with 38-88% sensitivity and 100% specificity. Of them, N-terminal parts of PvRBP2c (PVX_090325-1) and PvRBP2 like partial A (PVX_090330-1) elicited high antigenicity. In addition, the PvRBP2-like homologue B (PVX_116930) fragment was newly identified as high antigenicity and may be exploited as a potential antigenic candidate among the PvRBP family. The functional activity of the PvRBP family on merozoite invasion remains unknown.

Assessment of an automated capillary system for Plasmodium vivax microsatellite genotyping

BACKGROUND

Several platforms have been used to generate the primary data for microsatellite analysis of malaria parasite genotypes. Each has relative advantages but share a limitation of being time- and cost-intensive. A commercially available automated capillary gel cartridge system was assessed in the microsatellite analysis of Plasmodium vivax diversity in the Peruvian Amazon.

METHODS

The reproducibility and accuracy of a commercially-available automated capillary system, QIAxcel, was assessed using a sequenced PCR product of 227 base pairs. This product was measured 42 times, then 27 P. vivax samples from Peruvian Amazon subjects were analyzed with this instrument using five informative microsatellites. Results from the QIAxcel system were compared with a Sanger-type sequencing machine, the ABI PRISM(®) 3100 Genetic Analyzer.

RESULTS

Significant differences were seen between the sequenced amplicons and the results from the QIAxcel instrument. Different runs, plates and cartridges yielded significantly different results. Additionally, allele size decreased with each run by 0.045, or 1 bp, every three plates. QIAxcel and ABI PRISM systems differed in giving different values than those obtained by ABI PRISM, and too many (i.e. inaccurate) alleles per locus were also seen with the automated instrument.

CONCLUSIONS

While P. vivax diversity could generally be estimated using an automated capillary gel cartridge system, the data demonstrate that this system is not sufficiently precise for reliably identifying parasite strains via microsatellite analysis. This conclusion reached after systematic analysis was due both to inadequate precision and poor reproducibility in measuring PCR product size.

Parasites and human evolution

Our understanding of human evolutionary and population history can be advanced by ecological and evolutionary studies of our parasites. Many parasites flourish only in the presence of very specific human behaviors and in specific habitats, are wholly dependent on us, and have evolved with us for thousands or millions of years. Therefore, by asking when and how we first acquired those parasites, under which environmental and cultural conditions we are the most susceptible, and how the parasites have evolved and adapted to us and we in response to them, we can gain considerable insight into our own evolutionary history. As examples, the tapeworm life cycle is dependent on our consumption of meat, the divergence of body and head lice may have been subsequent to the development of clothing, and malaria hyperendemicity may be associated with agriculture. Thus, the evolutionary and population histories of these parasites are likely intertwined with critical aspects of human biology and culture. Here I review the mechanics of these and multiple other parasite proxies for human evolutionary history and discuss how they currently complement our fossil, archeological, molecular, linguistic, historical, and ethnographic records. I also highlight potential future applications of this promising model for the field of evolutionary anthropology.

Next-Generation Sequencing of Plasmodium vivax Patient Samples Shows Evidence of Direct Evolution in Drug-Resistance Genes

Understanding the mechanisms of drug resistance in Plasmodium vivax, the parasite that causes the most widespread form of human malaria, is complicated by the lack of a suitable long-term cell culture system for this parasite. In contrast to P. falciparum, which can be more readily manipulated in the laboratory, insights about parasite biology need to be inferred from human studies. Here we analyze the genomes of parasites within 10 human P. vivax infections from the Peruvian Amazon. Using next-generation sequencing we show that some P. vivax infections analyzed from the region are likely polyclonal. Despite their polyclonality we observe limited parasite genetic diversity by showing that three or fewer haplotypes comprise 94% of the examined genomes, suggesting the recent introduction of parasites into this geographic region. In contrast we find more than three haplotypes in putative drug-resistance genes, including the gene encoding dihydrofolate reductase-thymidylate synthase and the P. vivax multidrug resistance associated transporter, suggesting that resistance mutations have arisen independently. Additionally, several drug-resistance genes are located in genomic regions with evidence of increased copy number. Our data suggest that whole genome sequencing of malaria parasites from patients may provide more insight about the evolution of drug resistance than genetic linkage or association studies, especially in geographical regions with limited parasite genetic diversity.

Population Genetics, Evolutionary Genomics, and Genome-Wide Studies of Malaria: A View Across the International Centers of Excellence for Malaria Research

The study of the three protagonists in malaria-the Plasmodium parasite, the Anopheles mosquito, and the human host-is key to developing methods to control and eventually eliminate the disease. Genomic technologies, including the recent development of next-generation sequencing, enable interrogation of this triangle to an unprecedented level of scrutiny, and promise exciting progress toward real-time epidemiology studies and the study of evolutionary adaptation. We discuss the use of genomics by the International Centers of Excellence for Malaria Research, a network of field sites and laboratories in malaria-endemic countries that undertake cutting-edge research, training, and technology transfer in malarious countries of the world.

Plasmodium falciparum population structure in Sudan post artemisinin-based combination therapy

Over the past decade, Sudan has stepped up malaria control backed by WHO, and this has resulted in significant reduction in parasite rate, malaria morbidity and mortality. The present study analyzed Plasmodium falciparum parasites in four geographical separated areas, to examine whether the success in malaria control following the use of artemisinin-based combination therapy (ACT) has disrupted the population structure and evolution of the parasite. We examined 319 P. falciparum isolates collected between October 2009 and October 2012 in four different areas in Sudan (Jazira [central Sudan], Southern Darfur [western Sudan], Upper Nile [southern Sudan] and Kasala [eastern Sudan]). Twelve microsatellites were analyzed for allelic diversity, multi-locus haplotype and inter-population differentiation. Level of diversity was compared to that detected for three of the above microsatellites among P. falciparum parasites in central and eastern Sudan in 1999, prior to introduction of ACT. Diversity at each locus (unbiased heterozygosity [H]) was high in all areas (Jazira, H=0.67), (Southern Darfur, H=0.71), (Upper Nile, H=0.71), and (Kasala, H=0.63). Microsatellites were distributed widely and private alleles, detected in a single population, were rare. The extent of diversity in the above sites was similar to that seen, in 1999, in central (Khartoum, H=0.73) and eastern Sudan (Gedaref, H=0.75). Significant Linkage disequilibrium (LD) was observed between the microsatellites in all populations. Pairwise FST analysis revealed that parasites in the four areas could be considered as one population. However, the parasites in Sudan clustered away from parasites in West Africa and the Arabian Peninsula. Despite marked reduction in malaria risk in Sudan, the extent of diversity and parasite genetic structure are indicative of a large population size. Further considerable reduction in transmission would be needed before fragmented sub-population can be seen. In addition, the large divergence of P. falciparum in Sudan from West Africa and Arabian Peninsula populations may result from differential evolutionary pressures acting at the population level, which shall be considered in eradication plans.

Plasmodium vivax Diversity and Population Structure across Four Continents

Plasmodium vivax is the geographically most widespread human malaria parasite. To analyze patterns of microsatellite diversity and population structure across countries of different transmission intensity, genotyping data from 11 microsatellite markers was either generated or compiled from 841 isolates from four continents collected in 1999–2008. Diversity was highest in South-East Asia (mean allelic richness 10.0–12.8), intermediate in the South Pacific (8.1–9.9) Madagascar and Sudan (7.9–8.4), and lowest in South America and Central Asia (5.5–7.2). A reduced panel of only 3 markers was sufficient to identify approx. 90% of all haplotypes in South Pacific, African and SE-Asian populations, but only 60–80% in Latin American populations, suggesting that typing of 2–6 markers, depending on the level of endemicity, is sufficient for epidemiological studies. Clustering analysis showed distinct clusters in Peru and Brazil, but little sub-structuring was observed within Africa, SE-Asia or the South Pacific. Isolates from Uzbekistan were exceptional, as a near-clonal parasite population was observed that was clearly separated from all other populations (F_ST>0.2). Outside Central Asia F_ST values were highest (0.11–0.16) between South American and all other populations, and lowest (0.04–0.07) between populations from South-East Asia and the South Pacific. These comparisons between P. vivax populations from four continents indicated that not only transmission intensity, but also geographical isolation affect diversity and population structure. However, the high effective population size results in slow changes of these parameters. This persistency must be taken into account when assessing the impact of control programs on the genetic structure of parasite populations.

Plasmodium vivax is the predominant malaria parasite in Latin America, Asia and the South Pacific. Different factors are expected to shape diversity and population structure across continents, e.g. transmission intensity which is much lower in South America as compared to Southeast-Asia and the South Pacific, or geographical isolation of P. vivax populations in the South Pacific. We have compiled data from 841 isolates from South and Central America, Africa, Central Asia, Southeast-Asia and the South Pacific typed with a panel of 11 microsatellite markers. Diversity was highest in Southeast-Asia, where transmission is intermediate-high and migration of infected hosts is high, and lowest in South America and Central Asia where malaria transmission is low and focal. Reducing the panel of microsatellites showed that 2–6 markers are sufficient for genotyping for most drug trials and epidemiological studies, as these markers can identify >90% of all haplotypes. Parasites clustered according to continental origin, with high population differentiation between South American and Central Asian populations and the other populations, and lowest differences between Southeast-Asia and the South Pacific. Current attempts to reduce malaria transmission might change this pattern, but only after transmission is reduced for an extended period of time.

Vector sequence contamination of the Plasmodium vivax sequence database in PlasmoDB and In silico correction of 26 parasite sequences

We found a 47 aa protein sequence that occurs 17 times in the Plasmodium vivax nucleotide database published on PlasmoDB. Coding sequence analysis showed multiple restriction enzyme sites within the 141 bp nucleotide sequence, and a His6 tag attached to the 3' end, suggesting cloning vector origins. Sequences with vector contamination were submitted to NCBI, and BLASTN was used to cross-examine whole-genome shotgun contigs (WGS) from four recently deposited P. vivax whole genome sequencing projects. There are at least 26 genes listed in the PlasmoDB database that incorporate this cloning vector sequence into their predicted provisional protein products.

The evolutionary origins of Southeast Asian Ovalocytosis

Southeast Asian Ovalocytosis (SAO) is a common red blood cell disorder that is maintained as a balanced polymorphism in human populations. In individuals heterozygous for the SAO-causing mutation there are minimal detrimental effects and well-documented protection from severe malaria caused by Plasmodium vivax and Plasmodium falciparum; however, the SAO-causing mutation is fully lethal in utero when homozygous. The present-day high frequency of SAO in Island Southeast Asia indicates the trait is maintained by strong heterozygote advantage. Our study elucidates the evolutionary origin of SAO by characterizing DNA sequence variation in a 9.5 kilobase region surrounding the causal mutation in the SLC4A1 gene. We find substantial haplotype diversity among SAO chromosomes and estimate the age of the trait to be approximately 10,005 years (95% CI: 4930-23,200 years). This date is far older than any other human malaria-resistance trait examined previously in Southeast Asia, and considerably pre-dates the widespread adoption of agriculture associated with the spread of speakers of Austronesian languages some 4000 years ago. Using a genealogy-based method we find no evidence of historical positive selection acting on SAO (s=0.0, 95% CI: 0.0-0.03), in sharp contrast to the strong present-day selection coefficient (e.g., 0.09) estimated from the frequency of this recessively lethal trait. This discrepancy may be due to a recent increase in malaria-driven selection pressure following the spread of agriculture, with SAO targeted as a standing variant by positive selection in malarial populations.

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.

The utility of genomic data for Plasmodium vivax population surveillance

Genetic polymorphisms identified from genomic sequencing can be used to track changes in parasite populations through time. Such tracking is particularly informative when applying control strategies and evaluating their effectiveness. Using genomic approaches may also enable improved ability to categorise populations and to stratify them according to the likely effectiveness of intervention. Clinical applications of genomic approaches also allow relapses to be classified according to reinfection or recrudescence. These tools can be used not only to assess the effectiveness of malaria interventions but also to appraise the strategies for malaria elimination.

Uncovering the transmission dynamics of Plasmodium vivax using population genetics

Population genetic analysis of malaria parasites has the power to reveal key insights into malaria epidemiology and transmission dynamics with the potential to deliver tools to support control and elimination efforts. Analyses of parasite genetic diversity have suggested that Plasmodium vivax populations are more genetically diverse and less structured than those of Plasmodium falciparum indicating that P. vivax may be a more ancient parasite of humans and/or less susceptible to population bottlenecks, as well as more efficient at disseminating its genes. These population genetic insights into P. vivax transmission dynamics provide an explanation for its relative resilience to control efforts. Here, we describe current knowledge on P. vivax population genetic structure, its relevance to understanding transmission patterns and relapse and how this information can inform malaria control and elimination programmes.

Plasmodium vivax populations are more genetically diverse and less structured than sympatric Plasmodium falciparum populations

Introduction

The human malaria parasite, Plasmodium vivax, is proving more difficult to control and eliminate than Plasmodium falciparum in areas of co-transmission. Comparisons of the genetic structure of sympatric parasite populations may provide insight into the mechanisms underlying the resilience of P. vivax and can help guide malaria control programs.

Methodology/Principle findings

P. vivax isolates representing the parasite populations of four areas on the north coast of Papua New Guinea (PNG) were genotyped using microsatellite markers and compared with previously published microsatellite data from sympatric P. falciparum isolates. The genetic diversity of P. vivax (H_e = 0.83–0.85) was higher than that of P. falciparum (H_e = 0.64–0.77) in all four populations. Moderate levels of genetic differentiation were found between P. falciparum populations, even over relatively short distances (less than 50 km), with 21–28% private alleles and clear geospatial genetic clustering. Conversely, very low population differentiation was found between P. vivax catchments, with less than 5% private alleles and no genetic clustering observed. In addition, the effective population size of P. vivax (30353; 13043–69142) was larger than that of P. falciparum (18871; 8109–42986).

Conclusions/Significance

Despite comparably high prevalence, P. vivax had higher diversity and a panmictic population structure compared to sympatric P. falciparum populations, which were fragmented into subpopulations. The results suggest that in comparison to P. falciparum, P. vivax has had a long-term large effective population size, consistent with more intense and stable transmission, and limited impact of past control and elimination efforts. This underlines suggestions that more intensive and sustained interventions will be needed to control and eventually eliminate P. vivax. This research clearly demonstrates how population genetic analyses can reveal deeper insight into transmission patterns than traditional surveillance methods.

The neglected human malaria parasite Plasmodium vivax is responsible for a large proportion of the global malaria burden. Efforts to control malaria have revealed that P. vivax is more resilient than the other major human malaria parasite, Plasmodium falciparum. This study utilised population genetics to compare patterns of P. vivax and P. falciparum transmission in Papua New Guinea, a region where infection rates of the two species are similar. The results demonstrated that P. vivax populations are more genetically diverse than those of P. falciparum suggestive of a parasite population that is more resilient to environmental challenges, undergoing higher levels of interbreeding locally and between distant parasite populations. Unique characteristics of P. vivax such as relapse, which allows different strains from past infections to produce subsequent infections, may provide more opportunities for the exchange and dissemination of genetic material. The contrasting patterns observed for the two species may be the result of a differential impact of past elimination attempts and indicate that more rigorous interventions will be needed in efforts to control and eventually eliminate P. vivax.

Plasmodium vivax: N-terminal diversity in the blood stage SERA genes from Indian isolates

Worldwide malaria risk due to Plasmodium vivax makes development of vaccine against P. vivax, a high priority. Serine Repeat Antigen of P. vivax (PvSERA) is a multigene family of blood stage proteins with 12 homologues. Sequence diversity studies are important for understanding them as potential vaccine candidates. No information on N-terminal diversity of these genes is available in literature. In this paper, we evaluate the genetic polymorphism of N-terminal regions of the highly expressed member PvSERA4 and PvSERA5 genes from Indian field isolates. Our results show that PvSERA4 has deletions and insertions in Glutamine rich tetrameric repeat units contributing to its diversity. PvSERA5 also exhibits high genetic diversity with non-synonymous substitutions leading to identification of novel haplotypes from India. Our first report helps in elucidating the allelic variants of PvSERA genes in this region and contributes to evaluating their efficacy as vaccine candidates.

Plasmodium vivax and Plasmodium falciparum at the crossroads of exchange among islands in Vanuatu: implications for malaria elimination strategies

Understanding the transmission and movement of Plasmodium parasites is crucial for malaria elimination and prevention of resurgence. Located at the limit of malaria transmission in the Pacific, Vanuatu is an ideal candidate for elimination programs due to low endemicity and the isolated nature of its island setting. We analyzed the variation in the merozoite surface protein 1 (msp1) and the circumsporozoite protein (csp) of P. falciparum and P. vivax populations to examine the patterns of gene flow and population structures among seven sites on five islands in Vanuatu. Genetic diversity was in general higher in P. vivax than P. falciparum from the same site. In P. vivax, high genetic diversity was likely maintained by greater extent of gene flow among sites and among islands. Consistent with the different patterns of gene flow, the proportion of genetic variance found among islands was substantially higher in P. falciparum (28.81–31.23%) than in P. vivax (-0.53–3.99%). Our data suggest that the current island-by-island malaria elimination strategy in Vanuatu, while adequate for P. falciparum elimination, might need to be complemented with more centrally integrated measures to control P. vivax movement across islands.

Development of a single nucleotide polymorphism barcode to genotype Plasmodium vivax infections

Plasmodium vivax, one of the five species of Plasmodium parasites that cause human malaria, is responsible for 25–40% of malaria cases worldwide. Malaria global elimination efforts will benefit from accurate and effective genotyping tools that will provide insight into the population genetics and diversity of this parasite. The recent sequencing of P. vivax isolates from South America, Africa, and Asia presents a new opportunity by uncovering thousands of novel single nucleotide polymorphisms (SNPs). Genotyping a selection of these SNPs provides a robust, low-cost method of identifying parasite infections through their unique genetic signature or barcode. Based on our experience in generating a SNP barcode for P. falciparum using High Resolution Melting (HRM), we have developed a similar tool for P. vivax. We selected globally polymorphic SNPs from available P. vivax genome sequence data that were located in putatively selectively neutral sites (i.e., intergenic, intronic, or 4-fold degenerate coding). From these candidate SNPs we defined a barcode consisting of 42 SNPs. We analyzed the performance of the 42-SNP barcode on 87 P. vivax clinical samples from parasite populations in South America (Brazil, French Guiana), Africa (Ethiopia) and Asia (Sri Lanka). We found that the P. vivax barcode is robust, as it requires only a small quantity of DNA (limit of detection 0.3 ng/μl) to yield reproducible genotype calls, and detects polymorphic genotypes with high sensitivity. The markers are informative across all clinical samples evaluated (average minor allele frequency > 0.1). Population genetic and statistical analyses show the barcode captures high degrees of population diversity and differentiates geographically distinct populations. Our 42-SNP barcode provides a robust, informative, and standardized genetic marker set that accurately identifies a genomic signature for P. vivax infections.

Plasmodium vivax malaria is a major global public health problem, with nearly 2.5 billion people at risk for infection and approximately 132–391 million clinical infections annually. It has a wide geographical range, with a high disease burden in Asia, Central and South America, the Middle East, Oceania, and East Africa. Advances in sequencing technology and sample processing have made it possible to characterize the genetic diversity of P. vivax populations. This genetic variation provides a means to identify parasites by unique genetic signatures, or “barcodes.” We developed such a genetic barcode for P. vivax, composed of 42 robust and informative variants. Here we report its development and validation based on 87 clinical samples identified by microscopy to contain P. vivax from geographically diverse parasite populations from South America (Brazil, French Guiana), Africa (Ethiopia) and Asia (Sri Lanka). We show that the SNP barcode provides a genotyping tool that can be performed at low cost, providing a means to uniquely identify parasite infections and distinguish geographic origins, and that barcode data may offer new insights into P. vivax population structure and diversity.

Comparative analysis of field-isolate and monkey-adapted Plasmodium vivax genomes

Significant insights into the biology of Plasmodium vivax have been gained from the ability to successfully adapt human infections to non-human primates. P. vivax strains grown in monkeys serve as a renewable source of parasites for in vitro and ex vivo experimental studies and functional assays, or for studying in vivo the relapse characteristics, mosquito species compatibilities, drug susceptibility profiles or immune responses towards potential vaccine candidates. Despite the importance of these studies, little is known as to how adaptation to a different host species may influence the genome of P. vivax. In addition, it is unclear whether these monkey-adapted strains consist of a single clonal population of parasites or if they retain the multiclonal complexity commonly observed in field isolates. Here we compare the genome sequences of seven P. vivax strains adapted to New World monkeys with those of six human clinical isolates collected directly in the field. We show that the adaptation of P. vivax parasites to monkey hosts, and their subsequent propagation, did not result in significant modifications of their genome sequence and that these monkey-adapted strains recapitulate the genomic diversity of field isolates. Our analyses also reveal that these strains are not always genetically homogeneous and should be analyzed cautiously. Overall, our study provides a framework to better leverage this important research material and fully utilize this resource for improving our understanding of P. vivax biology.

In this study we compare the genome sequences of Plasmodium vivax collected directly from patients with those of parasites propagated in laboratory monkeys. We show that the adaptation and continuous propagation of Plasmodium vivax in monkeys does not induce systematic changes in the genome and, therefore, that these parasites constitute an unbiased resource for studying this important pathogen. Our analyses also reveal that some monkey-adapted Plasmodium vivax strains are not genetically homogenous and retain multiple genetically different parasites present in the original patient infection. Overall, our study confirms the utility of monkey-adapted Plasmodium vivax strains for malaria research but also shows that this resource should be analyzed cautiously as different samples of the same strain might provide different biological material.

In silico analysis for structure, function and T-cell epitopes of a hypothetical conserved (HP-C) protein coded by PVX_092425 in Plasmodium vivax

OBJECTIVE

Plasmodium spp. merozoite glycosylphosphatidylinositol-anchored proteins (GPI-APs) considered as protective immunogen in novel vaccines against malaria. To analyze the structure and function of a hypothetical conserved (HP-C) GPI-AP coded by gene PVX_092425 from Plasmodium vivax, and find its potential T-cell epitopes for further vivax malaria vaccine study.

METHODS

The structure, function and T-cell epitopes of the HP-C protein named Pvx_092425 were analyzed and predicted by online and offline bioinformatics software.

RESULTS

The bioinformatics data showed that the Pvx_092425 is an 830 amino acid (AA) long polypeptide encoded by five exons gene PVX_092425.It contains a pectin lyase-like superfamily, an AA repeats region, a cys-rich region and a transmembrane domain (TM) in C-terminal region. The alignment analysis drew it has a unique AA repeats region among Plasmodium spp. It was located in the cytoplasm, secretory system or cellular nucleus of P. vivax merozoite. For the sequence, the fragment of I823-V829 inserts in the interior side of the membrane, and M1--A812 belongs to the cytoplasmic tail. It has seven protein-protein binding sites. The peptides with the best predicted binding affinities were human leucocyte antigen (HLA) HLA-A*0203, HLA-DRB1*0101 and HLA- DRB1*0701.Among these predicted peptides, 582FLWDKALFD590 epitope interacted with HLA-DRB1*0101 allele showed best binding affinity compared to others. Structural analysis explained that the epitope fits well into the epitope-binding groove of HLA-DRB1*0101.

CONCLUSIONS

It proposes that the Pvx_092425 plays a key role during erythrocyte stage and generates information that is useful for development of blood-stage vaccine to block the merozoites invasion.

The biology of Plasmodium vivax explored through genomics

Malaria is a mosquito-borne disease caused by the Plasmodium parasite. Of the four Plasmodium species that routinely cause human malaria, Plasmodium vivax is the most widespread species outside Africa, causing ∼18.9 million cases in 2012. P. vivax cannot be cultured continuously in vitro, which severely hampers research in nonendemic and endemic countries alike. Consequently, whole-genome sequencing has become an effective means to interrogate the biology of the P. vivax parasite. Our comparative genomic analysis of five P. vivax reference genomes and several whole-genome sequences of the closely related monkey malaria species P. cynomolgi has revealed an extraordinary level of genetic diversity and enabled characterization of novel multigene families and important single-copy genes. The generation of whole-genome sequences from multiple clinical isolates is also driving forward knowledge concerning the biology and evolution of the species. Understanding the biology of P. vivax is crucial to develop potential antimalarial drugs and vaccines and to achieve the goal of eliminating malaria.

Paths to a malaria vaccine illuminated by parasite genomics

•

Discovery of vaccine candidate antigens by parasite genome sequence analyses.

•

Genetic crosses, linkage group selection, and functional studies on parasites.

•

Characterizing developmental and epigenetic variation alongside allelic polymorphism.

•

Selection by naturally acquired immune responses helps to focus vaccine design.

More human death and disease is caused by malaria parasites than by all other eukaryotic pathogens combined. As early as the sequencing of the first human genome, malaria parasite genomics was prioritized to fuel the discovery of vaccine candidate antigens. This stimulated increased research on malaria, generating new understanding of the cellular and molecular mechanisms of infection and immunity. This review of recent developments illustrates how new approaches in parasite genomics, and increasingly large amounts of data from population studies, are helping to identify antigens that are promising lead targets. Although these results have been encouraging, effective discovery and characterization need to be coupled with more innovation and funding to translate findings into newly designed vaccine products for clinical trials.

Genome-wide patterns of genetic polymorphism and signatures of selection in Plasmodium vivax

Plasmodium vivax is the most prevalent human malaria parasite outside of Africa. Yet, studies aimed to identify genes with signatures consistent with natural selection are rare. Here, we present a comparative analysis of the pattern of genetic variation of five sequenced isolates of P. vivax and its divergence with two closely related species, Plasmodium cynomolgi and Plasmodium knowlesi, using a set of orthologous genes. In contrast to Plasmodium falciparum, the parasite that causes the most lethal form of human malaria, we did not find significant constraints on the evolution of synonymous sites genome wide in P. vivax. The comparative analysis of polymorphism and divergence across loci allowed us to identify 87 genes with patterns consistent with positive selection, including genes involved in the “exportome” of P. vivax, which are potentially involved in evasion of the host immune system. Nevertheless, we have found a pattern of polymorphism genome wide that is consistent with a significant amount of constraint on the replacement changes and prevalent negative selection. Our analyses also show that silent polymorphism tends to be larger toward the ends of the chromosomes, where many genes involved in antigenicity are located, suggesting that natural selection acts not only by shaping the patterns of variation within the genes but it also affects genome organization.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/1475-2875-13-495) contains supplementary material, which is available to authorized users.

The impact of genomics on population genetics of parasitic diseases

Parasites, defined as eukaryotic microbes and parasitic worms that cause global diseases of human and veterinary importance, span many lineages in the eukaryotic Tree of Life. Historically challenging to study due to their complicated life-cycles and association with impoverished settings, their inherent complexities are now being elucidated by genome sequencing. Over the course of the last decade, projects in large sequencing centers, and increasingly frequently in individual research labs, have sequenced dozens of parasite reference genomes and field isolates from patient populations. This 'tsunami' of genomic data is answering questions about parasite genetic diversity, signatures of evolution orchestrated through anti-parasitic drug and host immune pressure, and the characteristics of populations. This brief review focuses on the state of the art of parasitic protist genomics, how the peculiar genomes of parasites are driving creative methods for their sequencing, and the impact that next-generation sequencing is having on our understanding of parasite population genomics and control of the diseases they cause.

A comprehensive evaluation of rodent malaria parasite genomes and gene expression

BACKGROUND

Rodent malaria parasites (RMP) are used extensively as models of human malaria. Draft RMP genomes have been published for Plasmodium yoelii, P. berghei ANKA (PbA) and P. chabaudi AS (PcAS). Although availability of these genomes made a significant impact on recent malaria research, these genomes were highly fragmented and were annotated with little manual curation. The fragmented nature of the genomes has hampered genome wide analysis of Plasmodium gene regulation and function.

RESULTS

We have greatly improved the genome assemblies of PbA and PcAS, newly sequenced the virulent parasite P. yoelii YM genome, sequenced additional RMP isolates/lines and have characterized genotypic diversity within RMP species. We have produced RNA-seq data and utilised it to improve gene-model prediction and to provide quantitative, genome-wide, data on gene expression. Comparison of the RMP genomes with the genome of the human malaria parasite P. falciparum and RNA-seq mapping permitted gene annotation at base-pair resolution. Full-length chromosomal annotation permitted a comprehensive classification of all subtelomeric multigene families including the 'Plasmodium interspersed repeat genes' (pir). Phylogenetic classification of the pir family, combined with pir expression patterns, indicates functional diversification within this family.

CONCLUSIONS

Complete RMP genomes, RNA-seq and genotypic diversity data are excellent and important resources for gene-function and post-genomic analyses and to better interrogate Plasmodium biology. Genotypic diversity between P. chabaudi isolates makes this species an excellent parasite to study genotype-phenotype relationships. The improved classification of multigene families will enhance studies on the role of (variant) exported proteins in virulence and immune evasion/modulation.

The distribution of incubation and relapse times in experimental human infections with the malaria parasite Plasmodium vivax

BACKGROUND

The distributions of incubation and relapse periods are key components of infectious disease models for the malaria parasite Plasmodium vivax; however, detailed distributions based upon experimental data are lacking.

METHODS

Using a range of historical, experimental mosquito-transmitted human infections, Bayesian estimation with non-informative priors was used to determine parametric distributions that can be readily implemented for the incubation period and time-to-first relapse in P. vivax infections, including global subregions by parasite source. These analyses were complemented with a pooled analysis of observational human infection data with infections that included malaria chemoprophylaxis and long-latencies. The epidemiological impact of these distributional assumptions was explored using stochastic epidemic simulations at a fixed reproductive number while varying the underlying distribution of incubation periods.

RESULTS

Using the Deviance Information Criteria to compare parameterizations, experimental incubation periods are most closely modeled with a shifted log-logistic distribution; a log-logistic mixture is the best fit for incubations in observational studies. The mixture Gompertz distribution was the best fit for experimental times-to-relapse among the tested parameterizations, with some variation by geographic subregions. Simulations suggest underlying distributional assumptions have critically important impacts on both the time-scale and total case counts within epidemics.

CONCLUSIONS

These results suggest that the exponential and gamma distributions commonly used for modeling incubation periods and relapse times inadequately capture the complexity in the distributions of event times in P. vivax malaria infections. In future models, log-logistic and Gompertz distributions should be utilized for general incubation periods and relapse times respectively, and region-specific distributions should be considered to accurately model and predict the epidemiology of this important human pathogen.

Global and local genetic diversity at two microsatellite loci in Plasmodium vivax parasites from Asia, Africa and South America

Background

Even though Plasmodium vivax has the widest worldwide distribution of the human malaria species and imposes a serious impact on global public health, the investigation of genetic diversity in this species has been limited in comparison to Plasmodium falciparum. Markers of genetic diversity are vital to the evaluation of drug and vaccine efficacy, tracking of P. vivax outbreaks, and assessing geographical differentiation between parasite populations.

Methods

The genetic diversity of eight P. vivax populations (n = 543) was investigated by using two microsatellites (MS), m1501 and m3502, chosen because of their seven and eight base-pair (bp) repeat lengths, respectively. These were compared with published data of the same loci from six other P. vivax populations.

Results

In total, 1,440 P. vivax samples from 14 countries on three continents were compared. There was highest heterozygosity within Asian populations, where expected heterozygosity (H_e) was 0.92-0.98, and alleles with a high repeat number were more common. Pairwise F_ST revealed significant differentiation between most P. vivax populations, with the highest divergence found between Asian and South American populations, yet the majority of the diversity (~89%) was found to exist within rather than between populations.

Conclusions

The MS markers used were informative in both global and local P. vivax population comparisons and their seven and eight bp repeat length facilitated population comparison using data from independent studies. A complex spatial pattern of MS polymorphisms among global P. vivax populations was observed which has potential utility in future epidemiological studies of the P. vivax parasite.

Electronic supplementary material

The online version of this article (doi:10.1186/1475-2875-13-392) contains supplementary material, which is available to authorized users.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/1475-2875-13-388) contains supplementary material, which is available to authorized users.

Genetic diversity and natural selection of three blood-stage 6-Cys proteins in Plasmodium vivax populations from the China-Myanmar endemic border

Pv12, Pv38 and Pv41, the three 6-Cys family proteins which are expressed in the blood-stage of vivax malaria, might be involved in merozoite invasion activity and thus be potential vaccine candidate antigens of Plasmodium vivax. However, little information is available concerning the genetic diversity and natural selection of these three proteins. In the present study, we analyzed the amino acid sequences of P. vivax blood-stage 6-Cys family proteins in comparison with the homologue proteins of Plasmodium cynomolgi strain B using bioinformatic methods. We also investigated genetic polymorphisms and natural selection of these three genes in P. vivax populations from the China-Myanmar endemic border. The three P. vivax blood-stage 6-Cys proteins were shown to possess a signal peptide at the N-terminus, containing two s48/45 domains, and Pv12 and Pv38 have a GPI-anchor motif at the C-terminus. Then, 22, 21 and 29 haplotypes of pv12, pv38 and pv41 were identified out of 45, 38 and 40 isolates, respectively. The dN/dS values for Domain II of pv38 and pv41 were 3.33880 and 5.99829, respectively, suggesting positive balancing selection for these regions. Meanwhile, the C-terminus of pv41 showed high nucleotide diversity, and Tajima's D test suggested that this fragment could be under positive balancing selection. Overall, our results have significant implications, providing a genetic basis for blood-stage malaria vaccine development based on these three 6-Cys proteins.

Genetic diversity of Plasmodium vivax over time and space: a community-based study in rural Amazonia

To examine how community-level genetic diversity of the malaria parasite Plasmodium vivax varies across time and space, we investigated the dynamics of parasite polymorphisms during the early phases of occupation of a frontier settlement in the Amazon Basin of Brazil. Microsatellite characterization of 84 isolates of P. vivax sampled over 3 years revealed a moderate-to-high genetic diversity (mean expected heterozygosity, 0.699), with a large proportion (78.5%) of multiple-clone infections (MCI), but also a strong multilocus linkage disequilibrium (LD) consistent with rare outcrossing. Little temporal and no spatial clustering was observed in the distribution of parasite haplotypes. A single microsatellite haplotype was shared by 3 parasites collected during an outbreak; all other 81 haplotypes were recovered only once. The lowest parasite diversity, with the smallest proportion of MCI and the strongest LD, was observed at the time of the outbreak, providing a clear example of epidemic population structure in a human pathogen. Population genetic parameters returned to pre-outbreak values during last 2 years of study, despite the concomitant decline in malaria incidence. We suggest that parasite genotyping can be useful for tracking the spread of new parasite strains associated with outbreaks in areas approaching malaria elimination.

Plasmodium falciparum and Plasmodium vivax specific lactate dehydrogenase: genetic polymorphism study from Indian isolates

Control and eradication of malaria is hindered by the acquisition of drug resistance by Plasmodium species. This has necessitated a persistent search for novel drugs and more efficient targets. Plasmodium species specific lactate dehydrogenase is one of the potential therapeutic and diagnostic targets, because of its indispensable role in endoerythrocytic stage of the parasite. A target molecule that is highly conserved in the parasite population can be more effectively used in diagnostics and therapeutics, hence, in the present study polymorphism in PfLDH (Plasmodiumfalciparum specific LDH) and PvLDH (Plasmodiumvivax specific LDH) genes was analyzed using PCR-single strand confirmation polymorphism (PCR-SSCP) and sequencing. Forty-six P. falciparum and thirty-five P. vivax samples were screened from different states of India. Our findings have revealed presence of a single PfLDH genotype and six PvLDH genotypes among the studied samples. Interestingly, along with synonymous substitutions, nonsynonymous substitutions were reported to be present for the first time in the PvLDH genotypes. Further, through amino acid sequence alignment and homology modeling studies we observed that the catalytic residues were conserved in all PvLDH genotypes and the nonsynonymous substitutions have not altered the enzyme structure significantly. Evolutionary genetics studies have confirmed that PfLDH and PvLDH loci are under strong purifying selection. Phylogenetic analysis of the pLDH gene sequences revealed that P. falciparum compared to P. vivax, has recent origin. The study therefore supports PfLDH and PvLDH as suitable therapeutic and diagnostic targets as well as phylogenetic markers to understand the genealogy of malaria species.

Distinct patterns of diversity, population structure and evolution in the AMA1 genes of sympatric Plasmodium falciparum and Plasmodium vivax populations of Papua New Guinea from an area of similarly high transmission

Background

As Plasmodium falciparum and Plasmodium vivax co-exist in most malaria-endemic regions outside sub-Saharan Africa, malaria control strategies in these areas must target both species in order to succeed. Population genetic analyses can predict the effectiveness of interventions including vaccines, by providing insight into patterns of diversity and evolution. The aim of this study was to investigate the population genetics of leading malaria vaccine candidate AMA1 in sympatric P. falciparum and P. vivax populations of Papua New Guinea (PNG), an area of similarly high prevalence (Pf = 22.3 to 38.8%, Pv = 15.3 to 31.8%).

Methods

A total of 72 Pfama1 and 102 Pvama1 sequences were collected from two distinct areas, Madang and Wosera, on the highly endemic PNG north coast.

Results

Despite a greater number of polymorphic sites in the AMA1 genes of P. falciparum (Madang = 52; Wosera = 56) compared to P. vivax (Madang = 36, Wosera = 34), the number of AMA1 haplotypes, haplotype diversity (Hd) and recombination (R) was far lower for P. falciparum (Madang = 12, Wosera = 20; Hd ≤0.92, R ≤45.8) than for P. vivax (Madang = 50, Wosera = 38; Hd = 0.99, R = ≤70.9). Balancing selection was detected only within domain I of AMA1 for P. vivax, and in both domains I and III for P. falciparum.

Conclusions

Higher diversity in the genes encoding P. vivax AMA1 than in P. falciparum AMA1 in this highly endemic area has important implications for development of AMA1-based vaccines in PNG and beyond. These results also suggest a smaller effective population size of P. falciparum compared to P. vivax, a finding that warrants further investigation. Differing patterns of selection on the AMA1 genes indicate that critical antigenic sites may differ between the species, highlighting the need for independent investigations of these two leading vaccine candidates.

A high resolution case study of a patient with recurrent Plasmodium vivax infections shows that relapses were caused by meiotic siblings

Plasmodium vivax infects a hundred million people annually and endangers 40% of the world's population. Unlike Plasmodium falciparum, P. vivax parasites can persist as a dormant stage in the liver, known as the hypnozoite, and these dormant forms can cause malaria relapses months or years after the initial mosquito bite. Here we analyze whole genome sequencing data from parasites in the blood of a patient who experienced consecutive P. vivax relapses over 33 months in a non-endemic country. By analyzing patterns of identity, read coverage, and the presence or absence of minor alleles in the initial polyclonal and subsequent monoclonal infections, we show that the parasites in the three infections are likely meiotic siblings. We infer that these siblings are descended from a single tetrad-like form that developed in the infecting mosquito midgut shortly after fertilization. In this natural cross we find the recombination rate for P. vivax to be 10 kb per centimorgan and we further observe areas of disequilibrium surrounding major drug resistance genes. Our data provide new strategies for studying multiclonal infections, which are common in all types of infectious diseases, and for distinguishing P. vivax relapses from reinfections in malaria endemic regions. This work provides a theoretical foundation for studies that aim to determine if new or existing drugs can provide a radical cure of P. vivax malaria.

Plasmodium vivax is capable of remaining dormant in the human liver for months to years after an initial infection, creating an asymptomatic human reservoir. This unique aspect of parasite biology makes eliminating P. vivax distinctly different from P. falciparum elimination, and yet very little is known about this dormant parasite stage. Lack of knowledge about the dormant liver stage prevents the creation of new drugs and public health interventions directed at P. vivax. In order to better understand this particular parasite stage, we used whole genome sequencing to analyze three sequential P. vivax infections, two of which could be definitively categorized as having arisen from dormant liver stages. Our whole genome sequencing data demonstrates that dormant liver stage parasites are closely related yet not, as had previously been postulated, identical. These data highlight the need for a new paradigm to investigate P. vivax dormant liver stages in order to design the next generation of P. vivax drugs and effective global health interventions.

The origin and diversification of the merozoite surface protein 3 (msp3) multi-gene family in Plasmodium vivax and related parasites

The genus Plasmodium is a diversified group of parasites with more than 200 known species that includes those causing malaria in humans. These parasites use numerous proteins in a complex process that allows them to invade the red blood cells of their vertebrate hosts. Many of those proteins are part of multi-gene families; one of which is the merozoite surface protein-3 (msp3) family. The msp3 multi-gene family is considered important in the two main human parasites, Plasmodium vivax and Plasmodium falciparum, as its paralogs are simultaneously expressed in the blood stage (merozoite) and are immunogenic. There are large differences among Plasmodium species in the number of paralogs in this family. Such differences have been previously explained, in part, as adaptations that allow the different Plasmodium species to invade their hosts. To investigate this, we characterized the array containing msp3 genes among several Plasmodium species, including P. falciparum and P. vivax. We first found no evidence indicating that the msp3 family of P. falciparum was homologous to that of P. vivax. Subsequently, by focusing on the diverse clade of nonhuman primate parasites to which P. vivax is closely related, where homology was evident, we found no evidence indicating that the interspecies variation in the number of paralogs was an adaptation related to changes in host range or host switches. Overall, we hypothesize that the evolution of the msp3 family in P. vivax is consistent with a model of multi-allelic diversifying selection where the paralogs may have functionally redundant roles in terms of increasing antigenic diversity. Thus, we suggest that the expressed MSP3 proteins could serve as "decoys", via antigenic diversity, during the critical process of invading the host red blood cells.

Identification of a new export signal in Plasmodium yoelii: identification of a new exportome

Development of the erythrocytic malaria parasite requires targeting of parasite proteins into multiple compartments located within and beyond the parasite confine. Beyond the PEXEL/VTS pathway and its characterized players, increasing amount of evidence has highlighted the existence of proteins exported using alternative export-signal(s)/pathway(s); hence, the exportomes currently predicted are incomplete. The nature of these exported proteins which could have a prominent role in most of the Plasmodium species remains elusive. Using P.  yoelii variant proteins, we identified a signal associated to lipophilic region that mediates export of P.  yoelii proteins. This non-PEXEL signal termed PLASMED is defined by semi-conserved residues and possibly a secondary structure. In vivo characterization of exported-proteins indicated that PLASMED is a bona fide export-signal that allowed us to identify an unseen P.  yoelii exportome. The repertoire of the newly predicted exported proteins opens up perspectives for unravelling the remodelling of the host-cell by the parasite, against which new therapies could be elaborated.

Dynamics of clonal diversity in natural infections of the malaria parasite Plasmodium mexicanum in its free-ranging lizard host

Within mixed-genotype infections of malaria parasites (Plasmodium), the number of genetic clones present is associated with variation in important life history traits of the infection, including virulence. Although the number of clones present is important, how the proportion of those clones varies over time is poorly known. Clonal proportions of the lizard malaria parasite, Plasmodium mexicanum, were assessed in naturally infected free-ranging lizards followed in a mark-recapture program over as long as two warm seasons, the typical life span of the lizard. Clonal proportions were determined by amplifying two microsatellite markers, a method previously verified for accuracy. Most blood samples had been stored for over a decade, so a verification test determined that these samples had not degraded. Although the environment experienced by the parasite (its host) varies over the seasons and transmission occurs over the entire warm season, 68% of infections were stable over time, harboring a single clone (37% of infections) or multiple clones changing only 1-12% maximum comparing any two samples (31% of infections). The maximum change seen in any infection (comparing any two sample periods) was only 30%. A new clone entered three infections (only once successfully), and a clone was lost in only three infections. These results mirror those seen for a previous study of experimentally induced infections that showed little change in relative proportions over time. The results of this study, the first look at how clonal proportions vary over time for any malaria parasite of a nonhuman vertebrate host for natural infections, were surprising because experimental studies show clones of P. mexicanum appear to interact, yet relative proportions of clones typically remain constant over time.

Identification of pathogen genomic variants through an integrated pipeline

BACKGROUND

Whole-genome sequencing represents a powerful experimental tool for pathogen research. We present methods for the analysis of small eukaryotic genomes, including a streamlined system (called Platypus) for finding single nucleotide and copy number variants as well as recombination events.

RESULTS

We have validated our pipeline using four sets of Plasmodium falciparum drug resistant data containing 26 clones from 3D7 and Dd2 background strains, identifying an average of 11 single nucleotide variants per clone. We also identify 8 copy number variants with contributions to resistance, and report for the first time that all analyzed amplification events are in tandem.

CONCLUSIONS

The Platypus pipeline provides malaria researchers with a powerful tool to analyze short read sequencing data. It provides an accurate way to detect SNVs using known software packages, and a novel methodology for detection of CNVs, though it does not currently support detection of small indels. We have validated that the pipeline detects known SNVs in a variety of samples while filtering out spurious data. We bundle the methods into a freely available package.

Analysis of the causes of spawning of large-scale, severe malarial epidemics and their rapid total extinction in western Provence, historically a highly endemic region of France (1745-1850)

BACKGROUND

The two main puzzles of this study are the onset and then sudden stopping of severe epidemics in western Provence (a highly malaria-endemic region of Mediterranean France) without any deliberate counter-measures and in the absence of significant population flux.

METHODS

Malaria epidemics during the period from 1745 to 1850 were analysed against temperature and rainfall records and several other potentially relevant factors.

RESULTS

Statistical analyses indicated that relatively high temperatures in early spring and in September/October, rainfall during the previous winter (principally December) and even from November to September and epidemics during the previous year could have played a decisive role in the emergence of these epidemics. Moreover, the epidemics were most likely not driven by other parameters (e.g., social, cultural, agricultural and geographical). Until 1776, very severe malarial epidemics affected large areas, whereas after this date, they were rarer and generally milder for local people and were due to canal digging activities. In the latter period, decreased rainfall in December, and more extreme and variable temperatures were observed. It is known that rainfall anomalies and temperature fluctuations may be detrimental to vector and parasite development.

CONCLUSION

This study showed the particular characteristics of malaria in historical Provence. Contrary to the situation in most other Mediterranean areas, Plasmodium falciparum was most likely not involved (during the years with epidemics, the mean temperature during the months of July and August, among other factors, did not play a role) and the population had no protective mutation. The main parasite species was Plasmodium vivax, which was responsible for very severe diseases, but contrary to in northern Europe, it is likely that transmission occurred only during the period where outdoor sporogony was possible, and P. vivax sporogony was always feasible, even during colder summers. Possible key elements in the understanding of the course of malaria epidemics include changes in the virulence of P. vivax strains, the refractoriness of anophelines and/or the degree or efficiency of acquired immunity. This study could open new lines of investigation into the comprehension of the conditions of disappearance/emergence of severe malaria epidemics in highly endemic areas.

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.

The Sri Lankan paradox: high genetic diversity in Plasmodium vivax populations despite decreasing levels of malaria transmission

Here we examined whether the recent dramatic decline in malaria transmission in Sri Lanka led to a major bottleneck in the local Plasmodium vivax population, with a substantial decrease in the effective population size. To this end, we typed 14 highly polymorphic microsatellite markers in 185 P. vivax patient isolates collected from 13 districts in Sri Lanka over a period of 5 years (2003-2007). Overall, we found a high degree of polymorphism, with 184 unique haplotypes (12-46 alleles per locus) and average genetic diversity (expected heterozygosity) of 0·8744. Almost 69% (n = 127) isolates had multiple-clone infections (MCI). Significant spatial and temporal differentiation (F ST = 0·04-0·25; P⩽0·0009) between populations was observed. The effective population size was relatively high but showed a decline from 2003-4 to 2006-7 periods (estimated as 45 661 to 22 896 or 10 513 to 7057, depending on the underlying model used). We used three approaches - namely, mode-shift in allele frequency distribution, detection of heterozygote excess and the M-ratio statistics - to test for evidence of a recent population bottleneck but only the low values of M-ratio statistics (ranging between 0·15-0·33, mean 0·26) were suggestive of such a bottleneck. The persistence of high genetic diversity and high proportion of MCI, with little change in effective population size, despite the collapse in demographic population size of P. vivax in Sri Lanka indicates the importance of maintaining stringent control and surveillance measures to prevent resurgence.

Operational research needs toward malaria elimination in China

Owing to the implementation of a national malaria elimination programme from 2010 to 2020, we performed a systematic review to assess research challenges in the People's Republic of China (P.R. China) and define research priorities in the next few years. A systematic search was conducted for articles published from January 2000 to December 2012 in international journals from PubMed and Chinese journals from the China National Knowledge Infrastructure (CNKI). In total, 2532 articles from CNKI and 308 articles from PubMed published between 2010 and 2012 related to malaria after unrelated references and review or comment were further excluded, and a set of research gaps have been identified that could hinder progress toward malaria elimination in P.R. China. For example, there is a lack of sensitive and specific tests for the diagnosis of malaria cases with low parasitemia, and there is a need for surveillance tools that can evaluate the epidemic status for guiding the elimination strategy. Hence, we argue that malaria elimination will be accelerated in P.R. China through the development of new tests, such as detection of parasite or drug resistance, monitoring glucose-6-phosphate dehydrogenase (G6PD) deficiency, active malaria screening methods, and understanding the effects of the environment and climate variation on vector distribution.

A call to arms: on refining Plasmodium vivax microsatellite marker panels for comparing global diversity

BACKGROUND

Microsatellite (MS) markers have become an important tool for studying the population diversity, evolutionary history and multiplicity of infection (MOI) of malaria parasite infections. MS are typically selected on the basis of being highly polymorphic. However, it is known that the polymorphic potential (mutability) of each marker can vary as much as two orders of magnitude, which radically changes how diversity is represented in the genome from one marker to the next. Over the past decade, approximately 240 Plasmodium vivax MS have been published, comprising nine major panels of markers. Inconsistent usage of each panel has resulted in a surfeit of descriptive genetic diversity data that are largely incomparable between populations. The objective of this study was to statistically evaluate the quality of individual MS markers in order to validate a refined panel of markers that will provide a balanced picture of P. vivax population diversity.

METHODS

All previously published data, including genetic diversity indices, MS parameters, and population parameters, were assembled from 18 different global studies into a flat file to facilitate statistical analysis and modelling using JMP® Genomics 6.0 (SAS Institute Inc, Cary, NC, USA). Statistical modeling was employed to down-select markers with extreme variation among the mean number of alleles, expected heterozygosity, maximum repeat length and/or chromosomal location of the repeat. Individual MS were analysed by step-down whole model linear regression and standard least squares fit models, both stratified by annual parasite incidence to identify MS markers with values significantly different from the mean.

RESULTS

Of the 42 MS under evaluation in this study, 18 (nine high priority) were identified as ideal candidates for measuring population diversity between global regions, while five (two high priority) additional markers were identified as candidates for MOI studies.

CONCLUSIONS

MS diversity was found to be a function of endemicity and motif structure. Evaluation of individual MS permitted the assembly of a refined panel of markers that can be reliably utilized in the field to compare population structures between global regions.

De novo assembly of a field isolate genome reveals novel Plasmodium vivax erythrocyte invasion genes

Recent sequencing of Plasmodium vivax field isolates and monkey-adapted strains enabled characterization of SNPs throughout the genome. These analyses relied on mapping short reads onto the P. vivax reference genome that was generated using DNA from the monkey-adapted strain Salvador I. Any genomic locus deleted in this strain would be lacking in the reference genome sequence and missed in previous analyses. Here, we report de novo assembly of a P. vivax field isolate genome. Out of 2,857 assembled contigs, we identify 362 contigs, each containing more than 5 kb of contiguous DNA sequences absent from the reference genome sequence. These novel P. vivax DNA sequences account for 3.8 million nucleotides and contain 792 predicted genes. Most of these contigs contain members of multigene families and likely originate from telomeric regions. Interestingly, we identify two contigs containing predicted protein coding genes similar to known Plasmodium red blood cell invasion proteins. One gene encodes the reticulocyte-binding protein gene orthologous to P. cynomolgi RBP2e and P. knowlesi NBPXb. The second gene harbors all the hallmarks of a Plasmodium erythrocyte-binding protein, including conserved Duffy-binding like and C-terminus cysteine-rich domains. Phylogenetic analysis shows that this novel gene clusters separately from all known Plasmodium Duffy-binding protein genes. Additional analyses showing that this gene is present in most P. vivax genomes and transcribed in blood-stage parasites suggest that P. vivax red blood cell invasion mechanisms may be more complex than currently understood. The strategy employed here complements previous genomic analyses and takes full advantage of next-generation sequencing data to provide a comprehensive characterization of genetic variations in this important malaria parasite. Further analyses of the novel protein coding genes discovered through de novo assembly have the potential to identify genes that influence key aspects of P. vivax biology, including alternative mechanisms of human erythrocyte invasion.

Plasmodium vivax is responsible for most malaria cases outside Africa, but is poorly understood, as the parasite is difficult to study in vitro. Genome sequencing studies offer a novel and exciting opportunity to better understand this parasite but, so far, have directly mapped reads onto the reference genome sequence generated from a single P. vivax strain. Here, we use sequence data generated from a field isolate to reconstruct long DNA sequences without relying on the reference genome. Our analyses reveal many P. vivax DNA sequences that are absent from the reference genome and contain 792 predicted genes. One of these novel genes encodes a predicted protein similar to known Plasmodium proteins involved in red blood cell invasion. This new gene is present in all P. vivax strains sequenced so far, except for the strain used to generate the reference genome, and is transcribed in blood-stage parasites. Overall, our analyses show that the catalogue of P. vivax genes was incomplete and that potentially important genes have been missed. We notably identified one putative invasion gene that seems functional and could dramatically change our understanding of the mechanisms determining red blood cell invasion by this important malaria parasite.