Extensive microsatellite diversity in the human malaria parasite Plasmodium vivax

Plasmodium vivax population structure and transmission dynamics in Sabah Malaysia

Despite significant progress in the control of malaria in Malaysia, the complex transmission dynamics of P. vivax continue to challenge national efforts to achieve elimination. To assess the impact of ongoing interventions on P. vivax transmission dynamics in Sabah, we genotyped 9 short tandem repeat markers in a total of 97 isolates (8 recurrences) from across Sabah, with a focus on two districts, Kota Marudu (KM, n = 24) and Kota Kinabalu (KK, n = 21), over a 2 year period. STRUCTURE analysis on the Sabah-wide dataset demonstrated multiple sub-populations. Significant differentiation (F_ST  = 0.243) was observed between KM and KK, located just 130 Km apart. Consistent with low endemic transmission, infection complexity was modest in both KM (mean MOI  = 1.38) and KK (mean MOI  = 1.19). However, population diversity remained moderate (H_E  = 0.583 in KM and H_E  = 0.667 in KK). Temporal trends revealed clonal expansions reflecting epidemic transmission dynamics. The haplotypes of these isolates declined in frequency over time, but persisted at low frequency throughout the study duration. A diverse array of low frequency isolates were detected in both KM and KK, some likely reflecting remnants of previous expansions. In accordance with clonal expansions, high levels of Linkage Disequilibrium (I_A^S >0.5 [P<0.0001] in KK and KM) declined sharply when identical haplotypes were represented once (I_A^S  = 0.07 [P = 0.0076] in KM, and I_A^S = -0.003 [P = 0.606] in KK). All 8 recurrences, likely to be relapses, were homologous to the prior infection. These recurrences may promote the persistence of parasite lineages, sustaining local diversity. In summary, Sabah's shrinking P. vivax population appears to have rendered this low endemic setting vulnerable to epidemic expansions. Migration may play an important role in the introduction of new parasite strains leading to epidemic expansions, with important implications for malaria elimination.

Primaquine for preventing relapse in people with Plasmodium vivax malaria treated with chloroquine

BACKGROUND

Plasmodium vivax infections are an important contributor to the malaria burden worldwide. The World Health Organization recommends a 14-day course of primaquine (0.25 mg/kg/day, giving an adult dose of 15 mg/day) to eradicate the liver stage of the parasite and prevent relapse of the disease. Many people find a 14-day primaquine regimen difficult to complete, and there is a potential risk of haemolytic anaemia in people with glucose-6-phosphate-dehydrogenase enzyme (G6PD) deficiency. This review evaluates primaquine in P. vivax, particularly alternatives to the standard 14-day course.

OBJECTIVES

To compare alternative primaquine regimens to the recommended 14-day regimen for preventing relapses (radical cure) in people with P. vivax malaria treated for blood stage infection with chloroquine. We also summarize trials comparing primaquine to no primaquine that led to the recommendation for the 14-day regimen.

SEARCH METHODS

We searched the Cochrane Infectious Diseases Group's Specialized Register, CENTRAL (The Cochrane Library), MEDLINE, EMBASE and LILACS up to 8 October 2013. We checked conference proceedings, trial registries and reference lists and contacted researchers and pharmaceutical companies for eligible studies.

SELECTION CRITERIA

Randomized controlled trials (RCTs) and quasi-RCTs comparing various primaquine dosing regimens with the standard primaquine regimen (15 mg/day for 14 days), or with no primaquine, in people with vivax malaria treated for blood stage infection with chloroquine.

DATA COLLECTION AND ANALYSIS

We independently assessed trial eligibility, trial quality, and extracted data. We calculated risk ratios (RR) with 95% confidence intervals (CI) for dichotomous data, and used the random-effects model in meta-analyses if there was significant heterogeneity. We assessed the overall quality of the evidence using the GRADE approach.

MAIN RESULTS

We included 15 trials (two cluster-RCTs) of 4377 adult and child participants. Most trials excluded people with G6PD deficiency. Trials compared various regimens of primaquine with the standard primaquine regimen, or with placebo or no treatment. All trials treated blood stage infection with chloroquine. Alternative primaquine regimens compared to 14-day primaquineRelapse rates were higher over six months with the five-day primaquine regimen than the standard 14-day regimen (RR 10.05, 95% CI 2.82 to 35.86; two trials, 186 participants, moderate quality evidence). Similarly, relapse over six months was higher with three days of primaquine than the standard 14-day regimen (RR 3.18, 95% CI 2.1 to 4.81; two trials, 262 participants, moderate quality evidence; six months follow-up); and with primaquine for seven days followed up over two months, compared to 14-day primaquine (RR 2.24, 95% CI 1.24 to 4.03; one trial, 126 participants, low quality evidence).Relapse with once-weekly supervised primaquine for eight weeks was little different over nine months follow-up compared to 14-day self-administered primaquine in one small study (RR 2.97, 95% CI 0.34 to 25.87; one trial, 129 participants, very low quality evidence). Primaquine regimens compared to no primaquineThe number of people that relapsed was similar between people given five days of primaquine or given placebo or no primaquine (four trials, 2213 participants, high quality evidence; follow-up six to 15 months); but lower with 14 days of primaquine (RR 0.6; 95% CI 0.48 to 0.75; ten trials, 1740 participants, high quality evidence; follow-up seven weeks to 15 months).No serious adverse events were reported. Treatment-limiting adverse events were rare and non-serious adverse events were mild and transient. Trial authors reported that people tolerated the drugs.We did not find trials comparing higher dose primaquine regimens (0.5 mg/kg/day or more) for five days or more with the 14-day regimen.

AUTHORS' CONCLUSIONS

The analysis confirms the current World Health Organization recommendation for 14-day primaquine (15 mg/day) to prevent relapse of vivax malaria. Shorter primaquine regimens at the same daily dose are associated with higher relapse rates. The comparative effects with weekly primaquine are promising, but require further trials to establish equivalence or non-inferiority compared to the 14-day regimen in high malaria transmission settings.

Neutral polymorphisms in putative housekeeping genes and tandem repeats unravels the population genetics and evolutionary history of Plasmodium vivax in India

The evolutionary history and age of Plasmodium vivax has been inferred as both recent and ancient by several studies, mainly using mitochondrial genome diversity. Here we address the age of P. vivax on the Indian subcontinent using selectively neutral housekeeping genes and tandem repeat loci. Analysis of ten housekeeping genes revealed a substantial number of SNPs (n = 75) from 100 P. vivax isolates collected from five geographical regions of India. Neutrality tests showed a majority of the housekeeping genes were selectively neutral, confirming the suitability of housekeeping genes for inferring the evolutionary history of P. vivax. In addition, a genetic differentiation test using housekeeping gene polymorphism data showed a lack of geographical structuring between the five regions of India. The coalescence analysis of the time to the most recent common ancestor estimate yielded an ancient TMRCA (232,228 to 303,030 years) and long-term population history (79,235 to 104,008) of extant P. vivax on the Indian subcontinent. Analysis of 18 tandem repeat loci polymorphisms showed substantial allelic diversity and heterozygosity per locus, and analysis of potential bottlenecks revealed the signature of a stable P. vivax population, further corroborating our ancient age estimates. For the first time we report a comparable evolutionary history of P. vivax inferred by nuclear genetic markers (putative housekeeping genes) to that inferred from mitochondrial genome diversity.

Microgeographical differences of Plasmodium vivax relapse and re-infection in the Peruvian Amazon

To determine the magnitude of Plasmodium vivax relapsing malaria in rural Amazonia, we carried out a study in four sites in northeastern Peru. Polymerase chain reaction-restriction fragment length polymorphism of PvMSP-3α and tandem repeat (TR) markers were compared for their ability to distinguish relapse versus reinfection. Of 1,507 subjects with P. vivax malaria, 354 developed > 1 episode during the study; 97 of 354 (27.5%) were defined as relapse using Pvmsp-3α alone. The addition of TR polymorphism analysis significantly reduced the number of definitively defined relapses to 26 of 354 (7.4%) (P < 0.05). Multivariate logistic regression modeling showed that the probability of having > 1 infection was associated with the following: subjects in Mazan (odds ratio [OR] = 2.56; 95% confidence interval [CI] 1.87, 3.51), 15-44 years of age (OR = 1.49; 95% CI 1.03, 2.15), traveling for job purposes (OR = 1.45; 95%CI 1.03, 2.06), and travel within past month (OR = 1.46; 95% CI 1.0, 2.14). The high discriminatory capacity of the molecular tools shown here is useful for understanding the micro-geography of malaria transmission.

A Large Plasmodium vivax Reservoir and Little Population Structure in the South Pacific

Introduction

The importance of Plasmodium vivax in malaria elimination is increasingly being recognized, yet little is known about its population size and population genetic structure in the South Pacific, an area that is the focus of intensified malaria control.

Methods

We have genotyped 13 microsatellite markers in 295 P. vivax isolates from four geographically distinct sites in Papua New Guinea (PNG) and one site from Solomon Islands, representing different transmission intensities.

Results

Diversity was very high with expected heterozygosity values ranging from 0.62 to 0.98 for the different markers. Effective population size was high (12′872 to 19′533 per site). In PNG population structuring was limited with moderate levels of genetic differentiation. F_ST values (adjusted for high diversity of markers) were 0.14–0.15. Slightly higher levels were observed between PNG populations and Solomon Islands (F_ST = 0.16).

Conclusions

Low levels of population structure despite geographical barriers to transmission are in sharp contrast to results from regions of low P. vivax endemicity. Prior to intensification of malaria control programs in the study area, parasite diversity and effective population size remained high.

The evolutionary history of Plasmodium vivax as inferred from mitochondrial genomes: parasite genetic diversity in the Americas

Plasmodium vivax is the most prevalent human malaria parasite in the Americas. Previous studies have contrasted the genetic diversity of parasite populations in the Americas with those in Asia and Oceania, concluding that New World populations exhibit low genetic diversity consistent with a recent introduction. Here we used an expanded sample of complete mitochondrial genome sequences to investigate the diversity of P. vivax in the Americas as well as in other continental populations. We show that the diversity of P. vivax in the Americas is comparable to that in Asia and Oceania, and we identify several divergent clades circulating in South America that may have resulted from independent introductions. In particular, we show that several haplotypes sampled in Venezuela and northeastern Brazil belong to a clade that diverged from the other P. vivax lineages at least 30,000 years ago, albeit not necessarily in the Americas. We propose that, unlike in Asia where human migration increases local genetic diversity, the combined effects of the geographical structure and the low incidence of vivax malaria in the Americas has resulted in patterns of low local but high regional genetic diversity. This could explain previous views that P. vivax in the Americas has low genetic diversity because these were based on studies carried out in limited areas. Further elucidation of the complex geographical pattern of P. vivax variation will be important both for diversity assessments of genes encoding candidate vaccine antigens and in the formulation of control and surveillance measures aimed at malaria elimination.

Higher microsatellite diversity in Plasmodium vivax than in sympatric Plasmodium falciparum populations in Pursat, Western Cambodia

Previous microsatellite analyses of sympatric populations of Plasmodium vivax and Plasmodium falciparum in Brazil revealed higher diversity in the former species. However, it remains unclear whether regional species-specific differences in prevalence and transmission levels might account for these findings. Here, we examine sympatric populations of P. vivax (n=87) and P. falciparum (n=164) parasites from Pursat province, Western Cambodia, where both species are similarly prevalent. Using 10 genome-wide microsatellites for P. falciparum and 13 for P. vivax, we found that the P. vivax population was more diverse than the sympatric P. falciparum population (average virtual heterozygosity [HE], 0.87 vs. 0.66, P=0.003), with more multiple-clone infections (89.6% vs. 47.6%) and larger mean number of alleles per marker (16.2 vs. 11.1, P=0.07). Both populations showed significant multi-locus linkage disequilibrium suggestive of a predominantly clonal mode of parasite reproduction. The higher microsatellite diversity found in P. vivax isolates, compared to sympatric P. falciparum isolates, does not necessarily result from local differences in transmission level and may reflect differences in population history between species or increased mutation rates in P. vivax.

Plasmodium vivax isolates from Cambodia and Thailand show high genetic complexity and distinct patterns of P. vivax multidrug resistance gene 1 (pvmdr1) polymorphisms

Plasmodium vivax accounts for an increasing fraction of malaria infections in Thailand and Cambodia. We compared P. vivax genetic complexity and antimalarial resistance patterns in the two countries. Use of a heteroduplex tracking assay targeting the merozoite surface protein 1 gene revealed that vivax infections in both countries are frequently polyclonal (84%), with parasites that are highly diverse (HE = 0.86) but closely related (GST = 0.18). Following a history of different drug policies in Thailand and Cambodia, distinct patterns of antimalarial resistance have emerged: most Cambodian isolates harbor the P. vivax multidrug resistance gene 1 (pvmdr1) 976F mutation associated with chloroquine resistance (89% versus 8%, P < 0.001), whereas Thai isolates more often display increased pvmdr1 copy number (39% versus 4%, P < 0.001). Finally, genotyping of paired isolates from individuals suspected of suffering relapse supports a complex scheme of relapse whereby recurrence of multiple identical variants is sometimes accompanied by the appearance of novel variants.

Microsatellite analysis of malaria parasites

Microsatellites have been increasingly used to investigate the population structure of malaria parasites, to map genetic loci contributing to phenotypes such as drug resistance and virulence in laboratory crosses and genome-wide association studies and to distinguish between treatment failures and new infections in clinical trials. Here, we provide optimized protocols for genotyping highly polymorphic microsatellites sampled from across the genomes of the human malaria parasites Plasmodium falciparum and P. vivax that have been extensively used in research laboratories worldwide.

Local population structure of Plasmodium: impact on malaria control and elimination

BACKGROUND

Regardless of the growing interest in detecting population structures in malarial parasites, there have been limited discussions on how to use this concept in control programmes. In such context, the effects of the parasite population structures will depend on interventions' spatial or temporal scales. This investigation explores the problem of identifying genetic markers, in this case microsatellites, to unveil Plasmodium genetic structures that could affect decisions in the context of elimination. The study was performed in a low-transmission area, which offers a good proxy to better understand problems associated with surveillance at the final stages of malaria elimination.

METHODS

Plasmodium vivax samples collected in Tumeremo, Venezuela, between March 2003 and November 2004 were analysed. Since Plasmodium falciparum also circulates in many low endemic areas, P. falciparum samples from the same locality and time period were included for comparison. Plasmodium vivax samples were assayed for an original set of 25 microsatellites and P. falciparum samples were assayed for 12 microsatellites.

RESULTS

Not all microsatellite loci assayed offered reliable local data. A complex temporal-cluster dynamics is found in both P. vivax and P. falciparum. Such dynamics affect the numbers and the type of microsatellites required for identifying individual parasites or parasite clusters when performing cross-sectional studies. The minimum number of microsatellites required to differentiate circulating P. vivax clusters differs from the minimum number of hyper-variable microsatellites required to distinguish individuals within these clusters. Regardless the extended number of microsatellites used in P. vivax, it was not possible to separate all individual infections.

CONCLUSIONS

Molecular surveillance has great potential; however, it requires preliminary local studies in order to properly interpret the emerging patterns in the context of elimination. Clonal expansions and clusters turnovers need to be taken into account when using molecular markers. Those affect the number and type of microsatellite markers, as well as, the expected genetic patterns in the context of operational investigations. By considering the local dynamics, elimination programmes could cost-effectively use molecular markers. However, population level studies need to consider the local limitations of a given set of loci in terms of providing epidemiologically relevant information.

Multiple-clone activation of hypnozoites is the leading cause of relapse in Plasmodium vivax infection

Background

Plasmodium vivax infection is characterized by a dormant hepatic stage, the hypnozoite that is activated at varying periods of time after clearance of the primary acute blood-stage, resulting in relapse. Differentiation between treatment failure and new infections requires characterization of initial infections, relapses, and clone multiplicity in vivax malaria infections.

Methodology/Principal Findings

Parasite DNA obtained from primary/relapse paired blood samples of 30 patients with P. vivax infection in Brazil was analyzed using 10 molecular markers (8 microsatellites and MSP-1 blocks 2 and 10). Cloning of PCR products and genotyping was used to identify low-frequency clones of parasites. We demonstrated a high frequency of multiple-clone infections in both primary and relapse infections. Few alleles were identified per locus, but the combination of these alleles produced many haplotypes. Consequently, the majority of parasites involved in relapse showed haplotypes that were distinct from those of primary infections. Plasmodium vivax relapse was characterized by temporal variations in the predominant parasite clones.

Conclusions/Significance

The high rate of low frequency alleles observed in both primary and relapse infections, along with temporal variation in the predominant alleles, might be the source of reported heterologous hypnozoite activation. Our findings complicate the concept of heterologous activation, suggesting the involvement of undetermined mechanisms based on host or environmental factors in the simultaneous activation of multiple clones of hypnozoites.

Phenotypic and genotypic characterisation of drug-resistant Plasmodium vivax

In this review we present recent developments in the analysis of Plasmodium vivax clinical trials and ex vivo drug-susceptibility assays, as well approaches currently being used to identify molecular markers of drug resistance. Clinical trials incorporating the measurement of in vivo drug concentrations and parasite clearance times are needed to detect early signs of resistance. Analysis of P. vivax growth dynamics ex vivo have defined the criteria for acceptable assay thresholds for drug susceptibility testing, and their subsequent interpretation. Genotyping and next-generation sequencing studies in P. vivax field isolates are set to transform our understanding of the molecular mechanisms of drug resistance.

Molecular analysis of reticulocyte binding protein-2 gene in Plasmodium vivax isolates from India

Background

Plasmodium vivax reticulocyte binding protein-2 (PvRBP-2) is a promising candidate for development of vaccine against parasite. DNA sequence polymorphism in pvrbp-2 which may hamper the vaccine development program has been identified in laboratory strains. Therefore, unraveling genetic polymorphism in pvrbp-2 from field isolates is a prerequisite for success in vaccine development. This study was designed with a primary aim to uncover genetic polymorphism in pvrbp-2 among P. vivax field isolates.

Results

Using virtual restriction mapping of pvrbp-2 sequences, two restriction enzymes (AluI and ApoI) were selected for the development of pvrbp-2 as a PCR-RFLP marker. Restriction fragment length polymorphism (RFLP) analysis revealed a high degree of genetic polymorphism in the pvrbp-2 gene among field isolates of P. vivax. ApoI-RFLP was found to be more efficient in identifying the extent of genetic polymorphism in pvrbp-2 compared to AluI-RFLP. Combined genotyping/haplotyping of RFLP pattern revealed a total of 36 distinct RFLP patterns among 83 P. vivax isolates analyzed. DNA sequence analysis also supports high degree of genetic polymorphism among field isolates of P. vivax. Pvrbp-2 PCR-RFLP method is able to distinguish multiple infection up to 16.86% and it revealed a low level of shared genetic pool between more than two populations.

Conclusion

The study suggests that pvrbp-2 is highly polymorphic genetic marker which can be used for population genetic analyses. RFLP analysis suggests presence of nearly similar proportion of Sal-1 and Belem alleles in Indian P. vivax populations. The larger extent of genetic polymorphism identified from limited samples advocates to screen genetic polymorphism in pvrbp-2 from malaria endemic geographical regions and countries for designing pvrbp-2 based anti-malarial control measures.

Evidence for a recent population bottleneck in an Apicomplexan parasite of caribou and reindeer, Besnoitia tarandi

The evolutionary history and epidemiology of parasites may be reflected in the extent and geographic distribution of their genetic variation. Among coccidian parasites, the population structure of only Toxoplasma gondii has been extensively examined. Intraspecific variation in other coccidia, for example, those assigned to the genus Besnoitia, remains poorly defined. Here, we characterize the extent of genetic variation among populations of Besnoitia tarandi, a parasite whose intermediate hosts include reindeer/caribou (Rangifer tarandus). Isolates from the Canadian Arctic and Finnish sub-Arctic were genotyped at six microsatellite loci, the first internal transcribed spacer region of nuclear rDNA, and the RNA polymerase β subunit (rpoB) encoded in the plastid genome. Remarkably, all isolates exhibited the same multilocus genotype, regardless of the isolate's geographic origin. This absolute monomorphism occurred despite the capacity of these loci to vary, as established by evident differentiation between B. tarandi and two other species of Besnoitia, and variation among four isolates of B. besnoiti. The surprising lack of genetic variation across the sampled range suggests that B. tarandi may have experienced a recent population bottleneck.

Plasmodium vivax populations revisited: mitochondrial genomes of temperate strains in Asia suggest ancient population expansion

BACKGROUND

Plasmodium vivax is the most widely distributed human malaria parasite outside of Africa, and its range extends well into the temperate zones. Previous studies provided evidence for vivax population differentiation, but temperate vivax parasites were not well represented in these analyses. Here we address this deficit by using complete mitochondrial (mt) genome sequences to elucidate the broad genetic diversity and population structure of P. vivax from temperate regions in East and Southeast Asia.

RESULTS

From the complete mtDNA sequences of 99 clinical samples collected in China, Myanmar and Korea, a total of 30 different haplotypes were identified from 26 polymorphic sites. Significant differentiation between different East and Southeast Asian parasite populations was observed except for the comparison between populations from Korea and southern China. Haplotype patterns and structure diversity analysis showed coexistence of two different groups in East Asia, which were genetically related to the Southeast Asian population and Myanmar population, respectively. The demographic history of P. vivax, examined using neutrality tests and mismatch distribution analyses, revealed population expansion events across the entire P. vivax range and the Myanmar population. Bayesian skyline analysis further supported the occurrence of ancient P. vivax population expansion.

CONCLUSIONS

This study provided further resolution of the population structure and evolution of P. vivax, especially in temperate/warm-temperate endemic areas of Asia. The results revealed divergence of the P. vivax populations in temperate regions of China and Korea from other populations. Multiple analyses confirmed ancient population expansion of this parasite. The extensive genetic diversity of the P. vivax populations is consistent with phenotypic plasticity of the parasites, which has implications for malaria control.

Molecular markers and genetic diversity of Plasmodium vivax

Enhanced understanding of the transmission dynamics and population genetics for Plasmodium vivax is crucial in predicting the emergence and spread of novel parasite phenotypes with major public health implications, such as new relapsing patterns, drug resistance and increased virulence. Suitable molecular markers are required for these population genetic studies. Here, we focus on two groups of molecular markers that are commonly used to analyse natural populations of P. vivax. We use markers under selective pressure, for instance, antigen-coding polymorphic genes, and markers that are not under strong natural selection, such as most minisatellite and microsatellite loci. First, we review data obtained using genes encoding for P. vivax antigens: circumsporozoite protein, merozoite surface proteins 1 and 3α, apical membrane antigen 1 and Duffy binding antigen. We next address neutral or nearly neutral molecular markers, especially microsatellite loci, providing a complete list of markers that have already been used in P. vivax populations studies. We also analyse the microsatellite loci identified in the P. vivax genome project. Finally, we discuss some practical uses for P. vivax genotyping, for example, detecting multiple-clone infections and tracking the geographic origin of isolates.

Individual Plasmodium vivax msp1 variants within polyclonal P. vivax infections display different propensities for relapse

Using a newly developed Plasmodium vivax merozoite surface protein 1 gene (Pvmsp1) heteroduplex tracking assay, we genotyped 107 P. vivax infections in individuals from Cambodia, 45 of whom developed recurrent parasitemia within 42 days. The majority of isolates were polyclonal, but recurrent parasitemias displayed fewer variants compared to initial parasitemias. Two Pvmsp1 gene variants occurred more frequently in the initial genotypes of those who developed recurrent parasitemia, representing the first time P. vivax variants associated with a higher risk of relapse have been described.

Understanding the population genetics of Plasmodium vivax is essential for malaria control and elimination

Traditionally, infection with Plasmodium vivax was thought to be benign and self-limiting, however, recent evidence has demonstrated that infection with P. vivax can also result in severe illness and death. Research into P. vivax has been relatively neglected and much remains unknown regarding the biology, pathogenesis and epidemiology of this parasite. One of the fundamental factors governing transmission and immunity is parasite diversity. An understanding of parasite population genetic structure is necessary to understand the epidemiology, diversity, distribution and dynamics of natural P. vivax populations. In addition, studying the population structure of genes under immune selection also enables investigation of the dynamic interplay between transmission and immunity, which is crucial for vaccine development. A lack of knowledge regarding the transmission and spread of P. vivax has been particularly highlighted in areas where malaria control and elimination programmes have made progress in reducing the burden of Plasmodium falciparum, yet P. vivax remains as a substantial obstacle. With malaria elimination back on the global agenda, mapping of global and local P. vivax population structure is essential prior to establishing goals for elimination and the roll-out of interventions. A detailed knowledge of the spatial distribution, transmission and clinical burden of P. vivax is required to act as a benchmark against which control targets can be set and measured. This paper presents an overview of what is known and what is yet to be fully understood regarding P. vivax population genetics, as well as the importance and application of P. vivax population genetics studies.

Plasmodium vivax lineages: geographical distribution, tandem repeat polymorphism, and phylogenetic relationship

Background

Multi-drug resistance and severe/complicated cases are the emerging phenotypes of vivax malaria, which may deteriorate current anti-malarial control measures. The emergence of these phenotypes could be associated with either of the two Plasmodium vivax lineages. The two lineages had been categorized as Old World and New World, based on geographical sub-division and genetic and phenotypical markers. This study revisited the lineage hypothesis of P. vivax by typing the distribution of lineages among global isolates and evaluated their genetic relatedness using a panel of new mini-satellite markers.

Methods

18S SSU rRNA S-type gene was amplified from 420 Plasmodium vivax field isolates collected from different geographical regions of India, Thailand and Colombia as well as four strains each of P. vivax originating from Nicaragua, Panama, Thailand (Pak Chang), and Vietnam (ONG). A mini-satellite marker panel was then developed to understand the population genetic parameters and tested on a sample subset of both lineages.

Results

18S SSU rRNA S-type gene typing revealed the distribution of both lineages (Old World and New World) in all geographical regions. However, distribution of Plasmodium vivax lineages was highly variable in every geographical region. The lack of geographical sub-division between lineages suggests that both lineages are globally distributed. Ten mini-satellites were scanned from the P. vivax genome sequence; these tandem repeats were located in eight of the chromosomes. Mini-satellites revealed substantial allelic diversity (7-21, AE = 14.6 ± 2.0) and heterozygosity (He = 0.697-0.924, AE = 0.857 ± 0.033) per locus. Mini-satellite comparison between the two lineages revealed high but similar pattern of genetic diversity, allele frequency, and high degree of allele sharing. A Neighbour-Joining phylogenetic tree derived from genetic distance data obtained from ten mini-satellites also placed both lineages together in every cluster.

Conclusions

The global lineage distribution, lack of genetic distance, similar pattern of genetic diversity, and allele sharing strongly suggested that both lineages are a single species and thus new emerging phenotypes associated with vivax malaria could not be clearly classified as belonging to a particular lineage on basis of their geographical origin.

Recent increase of genetic diversity in Plasmodium vivax population in the Republic of Korea

BACKGROUND

The reemergence of Plasmodium vivax in South Korea since 1993 represents a serious public health concern. Despite the importance in understanding genetic diversity for control strategies, however, studies remain inconclusive with the general premise that due to low rate of malaria transmission, there is generally low genetic diversity with very few strains involved. In this study, the genetic diversity and population structure of P. vivax in South Korea were explored by analysing microsatellite polymorphism.

METHODS

Sequences for 13 microsatellite loci distributed across the twelve chromosomes of P. vivax were obtained from 58 South Korean isolates collected during two sampling periods, namely 1997-2000 and 2007. The sequences were used for the analysis of expected heterozygosity and multilocus genotype diversity. Population structure was evaluated using STRUCTURE version 2.3.2. Linkage disequilibrium was also analysed to investigate the extent of outbreeding in the P. vivax population.

RESULTS

Mean expected heterozygosity significantly increased from 0.382 in 1997-2000 to 0.545 in 2007 (P < 0.05). The number of multilocus genotypes was 7 and 27; and genotype diversity was statistically significant (P < 0.01) at 0.661 and 0.995 in 1997-2000 and 2007, respectively. Analysis by STRUCTURE showed a more complex population structure in 2007 than in 1997-2000. Linkage disequilibrium between 13 microsatellites, although significant in both time points, was notably lower in 2007.

CONCLUSIONS

The present microsatellite analysis clearly showed recent increase of genetic diversity and recent relaxation of the strong population structure observed in 1997-2000. These results suggest that multiple genotypes not present previously recently migrated into South Korea, accompanied by substantial outbreeding between different genotypes.

High genetic polymorphism of relapsing P. vivax isolates in northwest Colombia

Genetic diversity of Plasmodium populations has been more extensively documented in Colombia for Plasmodium falciparum than for Plasmodium vivax. Recently, highly variable microsatellite markers have been described and used in population-level studies of genetic variation of P. vivax throughout the world. We applied this approach to understand the genetic structure of P. vivax populations and to identify recurrence-associated haplotypes. In this, three microsatellite markers of P. vivax were amplified and the combined size of the fragments was used to establish genotypes. Patients from an ongoing treatment efficacy trial who were kept either in endemic or non-endemic regions in the northwest of Colombia were included in the study. In total 58 paired clinical isolates, were amplified. A total of 54 haplotypes were observed among the two regions. Some haplotypes were exclusive to the endemic region where the highest degree of polymorphism was detected. In addition, we confirmed the different genotypes of recurrent-relapsing and primary infection isolates suggesting the activation of heterologous hypnozoite populations. We conclude that analysis of the three microsatellites is a valuable tool to establish the genetic characteristics of P. vivax populations in Colombia.

Characterization of treatment failure in efficacy trials of drugs against Plasmodium vivax by genotyping neutral and drug resistance-associated markers

Plasmodium vivax intervention trials customarily report uncorrected treatment failure rates. Application of recrudescence-reinfection genotyping and drug resistance single-nucleotide polymorphism typing to a 4-arm comparative efficacy trial illustrated that molecular approaches can assist in understanding the relative contributions of true drug resistance (recurrent with same genotype) and new infections to treatment failure. The PCR-corrected adequate clinical and parasitologic response may constitute an informative secondary endpoint in future P. vivax drug trials.

Genetic diversity of Plasmodium vivax Duffy Binding Protein II (PvDBPII) under unstable transmission and low intensity malaria in Sri Lanka

Elucidating the genetic diversity of the Duffy Binding Protein II (PvDBPII), a leading vaccine candidate for vivax malaria, in different geographical settings is vital. In Sri Lanka malaria transmission is unstable with low intensity. A relatively high level of allelic diversity, with 27 polymorphic nucleotides and 33 (aa) haplotypes was detected among the PvdbpII gene in 100 local Plasmodium vivax isolates collected from two hypoendemic areas, and from a non endemic area of the country. Mutations, recombination and balancing selection seem to maintain the observed local allelic diversity of PvdbpII. Lack of gene flow was evidenced by high Fst values between the two endemic study sites. Some of the aa polymorphisms may alter the binding and expression capacity of predicted T cell epitopes in PvDBPII. Of the 8 binding inhibitory linear B cell epitopes, 2 (H2 and M1) in the vicinity of the exact binding region of PvDBPII appeared to be highly conserved in Sri Lankan, Iran and Colombian isolates, while H3, M2, M3 and L3 neutralizing epitopes seem to be polymorphic globally, with H1 and L2 conserved in Colombian, South Korean and Iran isolates. In comparison to the reference Sal-1 strain, among 402 world-wide isolates (302 global and 100 local), 121 aa polymorphisms and 138 haplotypes were recorded of which 3 aa polymorphisms and 21 haplotypes seem to be unique to Sri Lanka. PvdbpII phylogeny suggests that local P. vivax parasites represent a sample of the global population. The ubiquitous presence of some PvDBPII aa haplotypes among both local and global P. vivax isolates may aid future vaccination strategies based on PvDBPII.

Consistent safety and infectivity in sporozoite challenge model of Plasmodium vivax in malaria-naive human volunteers

A safe and reproducible Plasmodium vivax infectious challenge method is required to evaluate the efficacy of malaria vaccine candidates. Seventeen healthy Duffy (+) and five Duffy (-) subjects were randomly allocated into three (A-C) groups and were exposed to the bites of 2-4 Anopheles albimanus mosquitoes infected with Plasmodium vivax derived from three donors. Duffy (-) subjects were included as controls for each group. Clinical manifestations of malaria and parasitemia were monitored beginning 7 days post-challenge. All Duffy (+) volunteers developed patent malaria infection within 16 days after challenge. Prepatent period determined by thick smear, was longer for Group A (median 14.5 d) than for Groups B and C (median 10 d/each). Infected volunteers recovered rapidly after treatment with no serious adverse events. The bite of as low as two P. vivax-infected mosquitoes provides safe and reliable infections in malaria-naive volunteers, suitable for assessing antimalarial and vaccine efficacy trials.

Plasmodium vivax sub-patent infections after radical treatment are common in Peruvian patients: results of a 1-year prospective cohort study

BACKGROUND

There is an increasing body of literature reporting treatment failure of the currently recommended radical treatment of Plasmodium vivax infections. As P. vivax is the main malaria species outside the African continent, emerging tolerance to its radical treatment regime could have major consequences in countries like Peru, where 80% of malaria cases are due to P. vivax. Here we describe the results of a 1-year longitudinal follow up of 51 confirmed P. vivax patients living around Iquitos, Peruvian Amazon, and treated according to the Peruvian national guidelines.

METHODOLOGY

Each month a blood sample for microscopy and later genotyping was systematically collected. Recent exposure to infection was estimated by detecting antibodies against the P. vivax circumsporozoite protein (CSP) and all PCR confirmed P. vivax infections were genotyped with 16 polymorphic microsatellites.

RESULTS

During a 1-year period, 84 recurrent infections, 22 positive also by microscopy, were identified, with a median survival time to first recurrent infection of 203 days. Most of them (71%) were asymptomatic; in 13 patients the infection persisted undetected by microscopy for several consecutive months. The genotype of mostly recurrent infections differed from that at day 0 while fewer differences were seen between the recurrent infections. The average expected heterozygosity was 0.56. There was strong linkage disequilibrium (I(A)(s) = 0.29, p<1.10(-4)) that remained also when analyzing only the unique haplotypes, suggesting common inbreeding.

CONCLUSION

In Peru, the P. vivax recurrent infections were common and displayed a high turnover of parasite genotypes compared to day 0. Plasmodium vivax patients, even when treated according to the national guidelines, may still represent an important parasite reservoir that can maintain transmission. Any elimination effort should consider such a hidden reservoir.

Limited polymorphism of the Plasmodium vivax merozoite surface protein 1 gene in isolates from Turkey

The 200-kD merozoite surface protein of Plasmodium vivax (PvMSP-1) is one of the leading vaccine candidates against P. vivax malaria. However, the gene encoding PvMSP-1 (pvmsp1) is highly polymorphic and is a major obstacle to effective vaccine development. To further understand polymorphism in pvmsp1, we obtained 30 full-length pvmsp1 sequences from southeastern Turkey. Comparative analysis of sequences from Turkey and other areas showed substantially limited polymorphism. Substitutions were found at 280 and 162 amino acid sites in samples from other regions and those from Turkey, respectively. Eight substitutions were unique to Turkey. In one of them, D/E at position 1706 in the C-terminal 19-kD region, the K/E change at 1709 was the only polymorphism previously known. Limited diversity was also observed in microsatellites. Data suggest a recent population bottleneck in Turkey that may have obscured a signature for balancing selection in the C-terminal 42-kD region, which was otherwise detectable in other areas.

Single-nucleotide polymorphism, linkage disequilibrium and geographic structure in the malaria parasite Plasmodium vivax: prospects for genome-wide association studies

Background

The ideal malaria parasite populations for initial mapping of genomic regions contributing to phenotypes such as drug resistance and virulence, through genome-wide association studies, are those with high genetic diversity, allowing for numerous informative markers, and rare meiotic recombination, allowing for strong linkage disequilibrium (LD) between markers and phenotype-determining loci. However, levels of genetic diversity and LD in field populations of the major human malaria parasite P. vivax remain little characterized.

Results

We examined single-nucleotide polymorphisms (SNPs) and LD patterns across a 100-kb chromosome segment of P. vivax in 238 field isolates from areas of low to moderate malaria endemicity in South America and Asia, where LD tends to be more extensive than in holoendemic populations, and in two monkey-adapted strains (Salvador-I, from El Salvador, and Belem, from Brazil). We found varying levels of SNP diversity and LD across populations, with the highest diversity and strongest LD in the area of lowest malaria transmission. We found several clusters of contiguous markers with rare meiotic recombination and characterized a relatively conserved haplotype structure among populations, suggesting the existence of recombination hotspots in the genome region analyzed. Both silent and nonsynonymous SNPs revealed substantial between-population differentiation, which accounted for ~40% of the overall genetic diversity observed. Although parasites clustered according to their continental origin, we found evidence for substructure within the Brazilian population of P. vivax. We also explored between-population differentiation patterns revealed by loci putatively affected by natural selection and found marked geographic variation in frequencies of nucleotide substitutions at the pvmdr-1 locus, putatively associated with drug resistance.

Conclusion

These findings support the feasibility of genome-wide association studies in carefully selected populations of P. vivax, using relatively low densities of markers, but underscore the risk of false positives caused by population structure at both local and regional levels.

See commentary: http://www.biomedcentral.com/1741-7007/8/90

Multilocus genotyping reveals high heterogeneity and strong local population structure of the Plasmodium vivax population in the Peruvian Amazon

Background

Peru is one of the Latin American countries with the highest malaria burden, mainly due to Plasmodium vivax infections. However, little is known about P. vivax transmission dynamics in the Peruvian Amazon, where most malaria cases occur. The genetic diversity and population structure of P. vivax isolates collected in different communities around Iquitos city, the capital of the Peruvian Amazon, was determined.

Methods

Plasmodium vivax population structure was determined by multilocus genotyping with 16 microsatellites on 159 P. vivax infected blood samples (mono-infections) collected in four sites around Iquitos city. The population characteristics were assessed only in samples with monoclonal infections (n = 94), and the genetic diversity was determined by calculating the expected heterozygosity and allelic richness. Both linkage disequilibrium and the genetic differentiation (θ) were estimated.

Results

The proportion of polyclonal infections varied substantially by site (11% - 70%), with the expected heterozygosity ranging between 0.44 and 0.69; no haplotypes were shared between the different populations. Linkage disequilibrium was present in all populations (I_A^S 0.14 - 0.61) but was higher in those with fewer polyclonal infections, suggesting inbreeding and a clonal population structure. Strong population differentiation (θ = 0.45) was found and the Bayesian inference cluster analysis identified six clusters based on distinctive allele frequencies.

Conclusion

The P. vivax populations circulating in the Peruvian Amazon basin are genetically diverse, strongly differentiated and they have a low effective recombination rate. These results are in line with the low and clustered pattern of malaria transmission observed in the region around Iquitos city.

Characterisation of pvmdr1 and pvdhfr genes associated with chemoresistance in Brazilian Plasmodium vivax isolates

Plasmodium vivax control is now being hampered by drug resistance. Orthologous Plasmodium falciparum genes linked to chloroquine or sulfadoxine-pyrimethamine chemoresistance have been identified in P. vivax parasites, but few studies have been performed. The goal of the present work is to characterise pvmdr1 and pvdhfr genes in parasite isolates from a Brazilian endemic area where no molecular investigation had been previously conducted. The pvmdr1 analysis revealed the existence of single (85.7%) and double (14.3%) mutant haplotypes, while the pvdhfr examination showed the presence of double (57.2%) and triple (42.8%) mutant haplotypes. The implications of these findings are discussed.

Microsatellite loci: determining the genetic variability of Plasmodium vivax

OBJECTIVE

To describe the genetic diversity of Plasmodium vivax isolates from different areas in the Brazilian Amazon using 11 polymorphic microsatellites and to evaluate the correlation between microsatellite variation and repeat array length.

METHODS

Microsatellites with variable repeat units and array lengths were selected using in silico search of the P. vivax genome. We designed primers and amplified the selected loci in DNA obtained from patients with P. vivax acute infections.

RESULTS

Positive correlation between repeat array length and microsatellite variation was detected independently of the size of repeat unit (di, tri, or tetranucleotide). We used these markers to describe the genetic variability of P. vivax isolates from four geographic regions of the Brazilian Amazon. Substantial variability was observed among P. vivax isolates within populations, concurrent with high levels of multiple-clone infections and high linkage disequilibrium. Overall, structured populations were observed with moderate to high genetic differentiation.

CONCLUSION

The markers studied are useful tools for assessing population structure of P. vivax, as demonstrated for Brazilian populations and for searching for evidence of recent selection events associated with different phenotypes, such as drug resistance.

High complexity of Plasmodium vivax infections in symptomatic patients from a rural community in central Vietnam detected by microsatellite genotyping

Fourteen published and three newly identified polymorphic microsatellites were used to genotype 69 Plasmodium vivax samples obtained from 39 patients detected over a period of two years who lived in a rural community of central Vietnam. All samples were polyclonal with an average expected heterozygosity of 0.86. Among the 39 patients, 16 experienced 1–5 recurrent episodes of P. vivax malaria, most of them (83%) with a different genotype profile compared with previous infections. The minimal set of microsatellites required for differentiating the genotype profiles of the recurrent infections compared with the full set of 17 microsatellites was explored. A combination of five markers was sufficient to identify all recurrent infections with an unrelated or different genotype profile compared with all previous episodes.

Geographic structure of Plasmodium vivax: microsatellite analysis of parasite populations from Sri Lanka, Myanmar, and Ethiopia

Genetic diversity and population structure of Plasmodium vivax parasites can predict the origin and spread of novel variants within a population enabling population specific malaria control measures. We analyzed the genetic diversity and population structure of 425 P. vivax isolates from Sri Lanka, Myanmar, and Ethiopia using 12 trinucleotide and tetranucleotide microsatellite markers. All three parasite populations were highly polymorphic with 3-44 alleles per locus. Approximately 65% were multiple-clone infections. Mean genetic diversity (H(E)) was 0.7517 in Ethiopia, 0.8450 in Myanmar, and 0.8610 in Sri Lanka. Significant linkage disequilibrium was maintained. Population structure showed two clusters (Asian and African) according to geography and ancestry. Strong clustering of outbreak isolates from Sri Lanka and Ethiopia was observed. Predictive power of ancestry using two-thirds of the isolates as a model identified 78.2% of isolates accurately as being African or Asian. Microsatellite analysis is a useful tool for mapping short-term outbreaks of malaria and for predicting ancestry.

From parasite genomes to one healthy world: Are we having fun yet?

In 1990, the Human Genome Sequencing Project was established. This laid the ground work for an explosion of sequence data that has since followed. As a result of this effort, the first complete genome of an animal, Caenorhabditis elegans was published in 1998. The sequence of Drosophila melanogaster was made available in March, 2000 and in the following year, working drafts of the human genome were generated with the completed sequence (92%) being released in 2003. Recent advancements and next-generation technologies have made sequencing common place and have infiltrated every aspect of biological research, including parasitology. To date, sequencing of 32 apicomplexa and 24 nematode genomes are either in progress or near completion, and over 600k nematode EST and 200k apicomplexa EST submissions fill the databases. However, the winds have shifted and efforts are now refocusing on how best to store, mine and apply these data to problem solving. Herein we tend not to summarize existing X-omics datasets or present new technological advances that promise future benefits. Rather, the information to follow condenses up-to-date-applications of existing technologies to problem solving as it relates to parasite research. Advancements in non-parasite systems are also presented with the proviso that applications to parasite research are in the making.

Large-scale genotyping and genetic mapping in Plasmodium parasites

The completion of many malaria parasite genomes provides great opportunities for genomewide characterization of gene expression and high-throughput genotyping. Substantial progress in malaria genomics and genotyping has been made recently, particularly the development of various microarray platforms for large-scale characterization of the Plasmodium falciparum genome. Microarray has been used for gene expression analysis, detection of single nucleotide polymorphism (SNP) and copy number variation (CNV), characterization of chromatin modifications, and other applications. Here we discuss some recent advances in genetic mapping and genomic studies of malaria parasites, focusing on the use of high-throughput arrays for the detection of SNP and CNV in the P. falciparum genome. Strategies for genetic mapping of malaria traits are also discussed.

Limited global diversity of the Plasmodium vivax merozoite surface protein 4 gene

Merozoite surface proteins (MSPs) of the malaria parasites are major candidates for vaccine development targeting asexual blood stages. However, the diverse antigenic repertoire of these antigens that induce strain-specific protective immunity in human is a major challenge for vaccine design and often determines the efficacy of a vaccine. Here we further assessed the genetic diversity of Plasmodium vivax MSP4 (PvMSP4) protein using 195 parasite samples collected mostly from Thailand, Indonesia and Brazil. Overall, PvMSP4 is highly conserved with only eight amino acid substitutions. The majority of the haplotype diversity was restricted to the two short tetrapeptide repeat arrays in exon 1 and 2, respectively. Selection and neutrality tests indicated that exon 1 and the entire coding region of PvMSP4 were under purifying selection. Despite the limited nucleotide polymorphism of PvMSP4, significant genetic differentiation among the three major parasite populations was detected. Moreover, microgeographical heterogeneity was also evident in the parasite populations from different endemic areas of Thailand.

The transcriptome of Plasmodium vivax reveals divergence and diversity of transcriptional regulation in malaria parasites

Plasmodium vivax causes over 100 million clinical infections each year. Primarily because of the lack of a suitable culture system, our understanding of the biology of this parasite lags significantly behind that of the more deadly species P. falciparum. Here, we present the complete transcriptional profile throughout the 48-h intraerythrocytic cycle of three distinct P. vivax isolates. This approach identifies strain specific patterns of expression for subsets of genes predicted to encode proteins associated with virulence and host pathogen interactions. Comparison to P. falciparum revealed significant differences in the expression of genes involved in crucial cellular functions that underpin the biological differences between the two parasite species. These data provide insights into the biology of P. vivax and constitute an important resource for the development of therapeutic approaches.

Polymorphism at the apical membrane antigen 1 locus reflects the world population history of Plasmodium vivax

Background

In malaria parasites (genus Plasmodium), ama-1 is a highly polymorphic locus encoding the Apical Membrane Protein-1, and there is evidence that the polymorphism at this locus is selectively maintained. We tested the hypothesis that polymorphism at the ama-1 locus reflects population history in Plasmodium vivax, which is believed to have originated in Southeast Asia and is widely geographically distributed. In particular, we tested for a signature of the introduction of P. vivax into the New World at the time of the European conquest and African slave trade and subsequent population expansion.

Results

One hundred and five ama-1 sequences were generated and analyzed from samples from six different Brazilian states and compared with database sequences from the Old World. Old World populations of P. vivax showed substantial evidence of population substructure, with high sequence divergence among localities at both synonymous and nonsynonymous sites, while Brazilian isolates showed reduced diversity and little population substructure.

Conclusion

These results show that genetic diversity in P. vivax AMA-1 reflects population history, with population substructure characterizing long-established Old World populations, whereas Brazilian populations show evidence of loss of diversity and recent population expansion.

Note

Nucleotide sequence data reported is this paper are available in the GenBank™ database under the accession numbers EF031154 – EF031216 and EF057446 – EF057487