Genetic diversity and dynamics of plasmodium falciparum and P. vivax populations in multiply infected children with asymptomatic malaria infections in Papua New Guinea

Measuring changes in Plasmodium falciparum census population size in response to sequential malaria control interventions

Here we introduce a new endpoint ″census population size″ to evaluate the epidemiology and control of Plasmodium falciparum infections, where the parasite, rather than the infected human host, is the unit of measurement. To calculate census population size, we rely on a definition of parasite variation known as multiplicity of infection (MOI var ), based on the hyper-diversity of the var multigene family. We present a Bayesian approach to estimate MOI var from sequencing and counting the number of unique DBLα tags (or DBLα types) of var genes, and derive from it census population size by summation of MOI var in the human population. We track changes in this parasite population size and structure through sequential malaria interventions by indoor residual spraying (IRS) and seasonal malaria chemoprevention (SMC) from 2012 to 2017 in an area of high-seasonal malaria transmission in northern Ghana. Following IRS, which reduced transmission intensity by > 90% and decreased parasite prevalence by ~40-50%, significant reductions in var diversity, MOI var , and population size were observed in ~2,000 humans across all ages. These changes, consistent with the loss of diverse parasite genomes, were short lived and 32-months after IRS was discontinued and SMC was introduced, var diversity and population size rebounded in all age groups except for the younger children (1-5 years) targeted by SMC. Despite major perturbations from IRS and SMC interventions, the parasite population remained very large and retained the var population genetic characteristics of a high-transmission system (high var diversity; low var repertoire similarity) demonstrating the resilience of P. falciparum to short-term interventions in high-burden countries of sub-Saharan Africa.

From multiplicity of infection to force of infection for sparsely sampled Plasmodium falciparum populations at high transmission

High multiplicity of infection or MOI, the number of genetically distinct parasite strains co-infecting a single human host, characterizes infectious diseases including falciparum malaria at high transmission. It accompanies high asymptomatic Plasmodium falciparum prevalence despite high exposure, creating a large transmission reservoir challenging intervention. High MOI and asymptomatic prevalence are enabled by immune evasion of the parasite achieved via vast antigenic diversity. Force of infection or FOI, the number of new infections acquired by an individual host over a given time interval, is the dynamic sister quantity of MOI, and a key epidemiological parameter for monitoring the impact of antimalarial interventions and assessing vaccine or drug efficacy in clinical trials. FOI remains difficult, expensive, and labor-intensive to accurately measure, especially in high-transmission regions, whether directly via cohort studies or indirectly via the fitting of epidemiological models to repeated cross-sectional surveys. We propose here the application of queuing theory to obtain FOI on the basis of MOI, in the form of either a two-moment approximation method or Little's law. We illustrate these methods with MOI estimates obtained under sparse sampling schemes with the recently proposed " v a r coding" method, based on sequences of the v a r multigene family encoding for the major variant surface antigen of the blood stage of malaria infection. The methods are evaluated with simulation output from a stochastic agent-based model, and are applied to an interrupted time-series study from Bongo District in northern Ghana before and immediately after a three-round transient indoor residual spraying (IRS) intervention. We incorporate into the sampling of the simulation output, limitations representative of those encountered in the collection of field data, including under-sampling of v a r genes, missing data, and usage of antimalarial drug treatment. We address these limitations in MOI estimates with a Bayesian framework and an imputation bootstrap approach. We demonstrate that both proposed methods give good and consistent FOI estimates across various simulated scenarios. Their application to the field surveys shows a pronounced reduction in annual FOI during intervention, of more than 70%. The proposed approach should be applicable to the many geographical locations where cohort or cross-sectional studies with regular and frequent sampling are lacking but single-time-point surveys under sparse sampling schemes are available, and for MOI estimates obtained in different ways. They should also be relevant to other pathogens of humans, wildlife and livestock whose immune evasion strategies are based on large antigenic variation resulting in high multiplicity of infection.

Dynamics of parasite growth in genetically diverse Plasmodium falciparum isolates

Multiple parasite lineages with different proliferation rates or fitness may coexist within a clinical malaria isolate, resulting in complex growth interactions and variations in phenotype. To elucidate the dynamics of parasite growth in multiclonal isolates, we measured growth rates (GRs) of three Plasmodium falciparum Cambodian isolates, including IPC_3445 (MRA-1236), IPC_5202 (MRA-1240), IPC_6403 (MRA-1285), and parasite lineages previously cloned from each of these isolates by limiting dilution. Following synchronization, in vitro cultures of each parasite line were maintained over four consecutive asexual cycles (192 h), with thin smears prepared at each 48-h cycle to estimate GR and fold change in parasitemia (FCP). Cell cycle time (CCT), the duration it takes for ring-stage parasites to develop into mature schizonts, was measured by monitoring the development of 0-3-h post-invasion rings for up to 52 h post-incubation. Laboratory lines 3D7 (MRA-102) and Dd2 (MRA-150) were used as controls. Significant differences in GR, FCP, and CCT were observed between parasite isolates and clonal lineages from each isolate. The parasite lines studied here have well-defined growth phenotypes and will facilitate basic malaria research and development of novel malaria interventions. These lines are available to malaria researchers through the MR4 collection of NIAID's BEI Resources Program.

Genetic diversity of Plasmodium vivax populations from the China-Myanmar border identified by genotyping merozoite surface protein markers

BACKGROUND

Parasite diversity and population structure influence malaria control measures. Malaria transmission at international borders affects indigenous residents and migrants, defying management efforts and resulting in malaria re-introduction. Here we aimed to determine the extent and distribution of genetic variations in Plasmodium vivax populations and the complexity of infections along the China-Myanmar border.

METHODS

We collected clinical P. vivax samples from local and migrant malaria patients from Laiza and Myitsone, Kachin State, Myanmar, respectively. We characterized the polymorphisms in two P. vivax merozoite surface protein markers, Pvmsp-3α and Pvmsp-3β, by PCR-restriction fragment length polymorphism (PCR-RFLP) analysis. We sought to determine whether these genetic markers could differentiate these two neighboring parasite populations.

RESULTS

PCR revealed three major size variants for Pvmsp-3α and four for Pvmsp-3β among the 370 and 378 samples, respectively. PCR-RFLP resolved 26 fragment-size alleles by digesting Pvmsp-3α with Alu I and Hha I and 28 alleles by digesting Pvmsp-3β with Pst I. PCR-RFLP analysis of Pvmsp-3α found that infections in migrant laborers from Myitsone bore more alleles than did infections in residents of Laiza, while such difference was not evident from genotyping Pvmsp-3β. Infections originating from these two places contained distinct but overlapping subpopulations of P. vivax. Infections from Myitsone had a higher multiplicity of infection as judged by the size of the Pvmsp-3α amplicons and alleles after Alu I/Hha I digestion.

CONCLUSIONS

Migrant laborers from Myitsone and indigenous residents from Laiza harbored overlapping but genetically distinct P. vivax parasite populations. The results suggested a more diverse P. vivax population in Myitsone than in the border town of Laiza. PCR-RFLP of Pvmsp-3α offers a convenient method to determine the complexity of P. vivax infections and differentiate parasite populations.

Reciprocal positive effects on parasitemia between coinfecting haemosporidian parasites in house sparrows

BACKGROUND

Hosts are often simultaneously infected with several parasite species. These co-infections can lead to within-host interactions of parasites, including mutualism and competition, which may affect both virulence and transmission. Birds are frequently co-infected with different haemosporidian parasites, but very little is known about if and how these parasites interact in natural host populations and what consequences there are for the infected hosts. We therefore set out to study Plasmodium and Haemoproteus parasites in house sparrows Passer domesticus with naturally acquired infections using a protocol where the parasitemia (infection intensity) is quantified by qPCR separately for the two parasites. We analysed infection status (presence/absence of the parasite) and parasitemia of parasites in the blood of both adult and juvenile house sparrows repeatedly over the season.

RESULTS

Haemoproteus passeris and Plasmodium relictum were the two dominating parasite species, found in 99% of the analyzed Sanger sequences. All birds were infected with both Plasmodium and Haemoproteus parasites during the study period. Seasonality explained infection status for both parasites in the adults: H. passeris was completely absent in the winter while P. relictum was present all year round. Among adults infected with H. passeris there was a positive effect of P. relictum parasitemia on H. passeris parasitemia and likewise among adults infected with P. relictum there was a positive effect of H. passeris parasitemia on P. relictum parasitemia. No such associations on parasitemia were seen in juvenile house sparrows.

CONCLUSIONS

The reciprocal positive relationships in parasitemia between P. relictum and H. passeris in adult house sparrows suggests either mutualistic interactions between these frequently occurring parasites or that there is variation in immune responses among house sparrow individuals, hence some individuals suppress the parasitemia of both parasites whereas other individuals suppress neither. Our detailed screening of haemosporidian parasites over the season shows that co-infections are very frequent in both juvenile and adult house sparrows, and since co-infections often have stronger negative effects on host fitness than the single infection, it is imperative to use screening systems with the ability to detect multiple parasites in ecological studies of host-parasite interactions.

Global diversity and balancing selection of 23 leading Plasmodium falciparum candidate vaccine antigens

Investigation of the diversity of malaria parasite antigens can help prioritize and validate them as vaccine candidates and identify the most common variants for inclusion in vaccine formulations. Studies of vaccine candidates of the most virulent human malaria parasite, Plasmodium falciparum, have focused on a handful of well-known antigens, while several others have never been studied. Here we examine the global diversity and population structure of leading vaccine candidate antigens of P. falciparum using the MalariaGEN Pf3K (version 5.1) resource, comprising more than 2600 genomes from 15 malaria endemic countries. A stringent variant calling pipeline was used to extract high quality antigen gene 'haplotypes' from the global dataset and a new R-package named VaxPack was used to streamline population genetic analyses. In addition, a newly developed algorithm that enables spatial averaging of selection pressure on 3D protein structures was applied to the dataset. We analysed the genes encoding 23 leading and novel candidate malaria vaccine antigens including csp, trap, eba175, ama1, rh5, and CelTOS. Our analysis shows that current malaria vaccine formulations are based on rare haplotypes and thus may have limited efficacy against natural parasite populations. High levels of diversity with evidence of balancing selection was detected for most of the erythrocytic and pre-erythrocytic antigens. Measures of natural selection were then mapped to 3D protein structures to predict targets of functional antibodies. For some antigens, geographical variation in the intensity and distribution of these signals on the 3D structure suggests adaptation to different human host or mosquito vector populations. This study provides an essential framework for the diversity of P. falciparum antigens to be considered in the design of the next generation of malaria vaccines.

'Bouncing Back' From Subclinical Malaria: Inflammation and Erythrocytosis After Resolution of P. falciparum Infection in Gambian Children

Recent malaria is associated with an increased risk of systemic bacterial infection. The aetiology of this association is unclear but malaria-related haemolysis may be one contributory factor. To characterise the physiological consequences of persistent and recently resolved malaria infections and associated haemolysis, 1650 healthy Gambian children aged 8-15 years were screened for P. falciparum infection (by 18sRNA PCR) and/or anaemia (by haematocrit) at the end of the annual malaria transmission season (t₁). P. falciparum-infected children and children with moderate or severe anaemia (haemoglobin concentration < 11g/dl) were age matched to healthy, uninfected, non-anaemic controls and screened again 2 months later (t₂). Persistently infected children (PCR positive at t₁ and t₂) had stable parasite burdens and did not differ significantly haematologically or in terms of proinflammatory markers from healthy, uninfected children. However, among persistently infected children, IL-10 concentrations were positively correlated with parasite density suggesting a tolerogenic response to persistent infection. By contrast, children who naturally resolved their infections (positive at t₁ and negative at t₂) exhibited mild erythrocytosis and concentrations of pro-inflammatory markers were raised compared to other groups of children. These findings shed light on a 'resetting' and potential overshoot of the homeostatic haematological response following resolution of malaria infection. Interestingly, the majority of parameters tested were highly heterogeneous in uninfected children, suggesting that some may be harbouring cryptic malaria or other infections.

Genetic diversity of Plasmodium vivax merozoite surface protein-3 alpha and beta from diverse geographic areas of Thailand

Malaria is parasitic disease cause by Plasmodium infection. In Thailand, co-infections of Plasmodium vivax and P. falciparum are commonly found. P. vivax infection has been increasing in the past decade. The objective of this study was to investigate the genetic diversity patterns of Plasmodium vivax merozoite surface protein 3 (PvMSP-3) genes in total of 450 isolates collected from Thai-neighboring border during two different periods (2009-20 14 and 2015 -2016) using polymerase chain reaction (PCR) - restriction fragment length polymorphism (RFLP) method. Three major types of PvMSP-3α (A, B, and C) and PvMSP-3β (A, B, and C) were detected based on PCR products size. Forty five and 23 of PvMSP-3α and, 41 and 30 alleles of PvMSP-3β genes from the first period and second period, respectively, with difference frequencies of samples were distinguished. The results strongly indicate genetic diversity patterns of PvMSP-3 in the second period especially samples from Thai-Myanmar border. These two polymorphic genes could be used as a molecular epidemiologic marker for genotyping P. vivax isolate in Thailand.

Allelic variation of msp-3α gene in Plasmodium vivax isolates and its correlation with the severity of disease in vivax malaria

Malaria is a global socio-economic burden of which Plasmodium vivax contributes for about 70-80 million cases on an annual basis worldwide and 60-65% cases in India. Diversity observed in highly polymorphic Merozoite Surface Protein-3α (msp-3α) encoded by MSP-3 gene family, has been used efficiently for genotyping of P. vivax infection. This study aims to correlate the severity of clinical symptoms with parasite load, genotype of P. vivax and multiplicity of infection. Based on clinical symptoms classification, 31 (67.9%) out of 46 cases were found to be severe while 15 (32.6%) were non-severe and correlation of the severity of vivax infection with parasite load was not observed. Analysis of msp3-α allele genotype showed that out of 31 severe cases, 19 (61.2%) were single-clone infection cases whereas 12 (38.7%) were multi-clone infections. Similarly, out of 15 non-severe cases, 9 (60%) were single clone and 6 (40%) were multi-clone infections indicating the absence of a correlation between the multiplicity of infection and disease severity. Allele frequency observed was 65.9%, 23.4%, 23.4%, and 28.2% for allele A, B, C and D, respectively. An important finding was the greater distribution of allele D than alleles B and C, which has been reported as a rare allele otherwise. Further, of 13 cases with allele D, 76.9% (10/13) cases were severe. This study showed the absence of a correlation between the severity of clinical symptoms with parasite load and multiplicity of infection but at the same time drives a possibility of severe vivax malarial symptoms to have an association with the persistence of allele D in the population. This upon exploration can lead to the development of a target in detection of severe cases of malaria.

Estimation of parasite age and synchrony status in Plasmodium falciparum infections

Human malaria parasites have complex but poorly understood population dynamics inside their human host. In some but not all infections, parasites progress synchronously through the 48 h lifecycle following erythrocyte invasion, such that at any one time there is a limited spread of parasites at a particular time (hours) post-invasion. Patients presenting with older parasites, and with asynchronous infections, have been reported to have higher risks of fatal outcomes, associated with higher parasite biomass and multiplication rates respectively. However, practical tools to assess synchrony and estimate parasite age post-invasion in patient samples are lacking. We have developed a novel method based on three genes differentially expressed over the parasite intra-erythrocytic lifecycle, and applied it to samples from patients with uncomplicated malaria attending two health clinics in Ghana. We found that most patients presented with synchronous infections, and with parasites within 12 h of erythrocyte invasion. Finally we investigated if clinical features such as fever and parasite density could act as predictors of parasite age and synchrony. The new method is a simple and practicable approach to study parasite dynamics in naturally-infected patients, and is a significant improvement on the subjective microscopical methods for parasite staging in vivo, aiding patient management.

Molecular surveillance over 14 years confirms reduction of Plasmodium vivax and falciparum transmission after implementation of Artemisinin-based combination therapy in Papua, Indonesia

Genetic epidemiology can provide important insights into parasite transmission that can inform public health interventions. The current study compared long-term changes in the genetic diversity and structure of co-endemic Plasmodium falciparum and P. vivax populations. The study was conducted in Papua Indonesia, where high-grade chloroquine resistance in P. falciparum and P. vivax led to a universal policy of Artemisinin-based Combination Therapy (ACT) in 2006. Microsatellite typing and population genetic analyses were undertaken on available isolates collected between 2004 and 2017 from patients with uncomplicated malaria (n = 666 P. falciparum and n = 615 P. vivax). The proportion of polyclonal P. falciparum infections fell from 28% (38/135) before policy change (2004-2006) to 18% (22/125) at the end of the study (2015-2017); p<0.001. Over the same period, polyclonal P. vivax infections fell from 67% (80/119) to 35% (33/93); p<0.001. P. falciparum strains persisted for up to 9 years compared to 3 months for P. vivax, reflecting higher rates of outbreeding in the latter. Sub-structure was observed in the P. falciparum population, but not in P. vivax, confirming different patterns of outbreeding. The P. falciparum population exhibited 4 subpopulations that changed in frequency over time. Notably, a sharp rise was observed in the frequency of a minor subpopulation (K2) in the late post-ACT period, accounting for 100% of infections in late 2016-2017. The results confirm epidemiological evidence of reduced P. falciparum and P. vivax transmission over time. The smaller change in P. vivax population structure is consistent with greater outbreeding associated with relapsing infections and highlights the need for radical cure to reduce recurrent infections. The study emphasizes the challenge in disrupting P. vivax transmission and demonstrates the potential of molecular data to inform on the impact of public health interventions.

Clinical Implications of Asymptomatic Plasmodium falciparum Infections in Malawi

Background

Asymptomatic Plasmodium falciparum infections are common in Malawi; however, the implications of these infections for the burden of malaria illness are unknown. Whether asymptomatic infections eventually progress to malaria illness, persist without causing symptoms, or clear spontaneously remains undetermined. We identified asymptomatic infections and evaluated the associations between persistent asymptomatic infections and malaria illness.

Methods

Children and adults (N = 120) who presented at a health facility with uncomplicated malaria were followed monthly for 2 years. During follow-up visits, participants with malaria symptoms were tested and, if positive, treated. Samples from all visits were tested for parasites using both microscopy and polymerase chain reaction, and all malaria infections underwent genotyping. Cox frailty models were used to estimate the temporal association between asymptomatic infections and malaria illness episodes. Mixed models were used to estimate the odds of clinical symptoms associated with new versus persistent infections.

Results

Participants had a median follow-up time of 720 days. Asymptomatic infections were detected during 23% of visits. Persistent asymptomatic infections were associated with decreased risk of malaria illness in all ages (hazard ratio 0.50, P < .001). When asymptomatic infections preceded malaria illness, newly-acquired infections were detected at 92% of subsequent clinical episodes, independent of presence of persistent infections. Malaria illness among children was more likely due to newly-acquired infections (odds ratio, 1.4; 95% confidence interval, 1.3-1.5) than to persistent infections.

Conclusions

Asymptomatic P. falciparum infections are associated with decreased incidence of malaria illness, but do not protect against disease when new infection occurs.

Dynamics of Afebrile Plasmodium falciparum Infections in Mozambican Men

Background

Afebrile Plasmodium falciparum infections usually remain undetected and untreated in the community and could potentially contribute to sustaining local malaria transmission in areas aiming for malaria elimination.

Methods

Thirty-two men with afebrile P. falciparum infections detected with rapid diagnostic test (RDTs) were followed for 28 days. Kaplan-Meier estimates were computed to estimate probability of parasite positivity and of reducing parasitemia by half of its initial level by day 28. Trends of parasite densities quantified by microscopy and real-time quantitative polymerase chain reaction (qPCR) were assessed using Poisson regression models, and the microscopy-to-qPCR positivity ratio was calculated at each time point. Three survival distributions (Gompertz, Weibull, and gamma) were used to evaluate their strength of fit to the data and to predict the median lifetime of infection.

Results

The cumulative probability of parasite qPCR positivity by day 28 was 81% (95% confidence interval [CI], 60.2–91.6). Geometric mean parasitemia at recruitment was 516.1 parasites/μL and fell to <100 parasites/μL by day 3, reaching 56.7 parasites/μL on day 28 (P < .001). The ratio of P. falciparum–positive samples by microscopy to qPCR decreased from 0.9 to 0.52 from recruitment to day 28. The best model fit to the data was obtained assuming a Gompertz distribution.

Conclusions

Afebrile P. falciparum infections detectable by RDT in semi-immune adults fall and stabilize at low-density levels during the first 4 days after detection, suggesting a rapid decline of potential transmissibility in this hidden parasite reservoir.

Clincial trials registration

NCT02698748

The Challenge of Quantifying Synchrony in Malaria Parasites

Malaria infection is often accompanied by periodic fevers, triggered by synchronous cycles of parasite replication within the host. The degree of synchrony in parasite development influences the efficacy of drugs and immune defenses and is therefore relevant to host health and infectiousness. Synchrony is thought to vary over the course of infection and across different host-parasite genotype or species combinations, but the evolutionary significance - if any - of this diversity remains elusive. Standardized methods are lacking, but the most common metric for quantifying synchrony is the percentage of parasites in a particular developmental stage. We use a heuristic model to show that this metric is often unacceptably biased. Methodological challenges must be addressed to characterize diverse patterns of synchrony and their consequences for disease severity and spread.

The Antibody Response to Plasmodium falciparum: Cues for Vaccine Design and the Discovery of Receptor-Based Antibodies

Plasmodium falciparum remains a serious public health problem and a continuous challenge for the immune system due to the complexity and diversity of the pathogen. Recent advances from several laboratories in the characterization of the antibody response to the parasite have led to the identification of critical targets for protection and revealed a new mechanism of diversification based on the insertion of host receptors into immunoglobulin genes, leading to the production of receptor-based antibodies. These advances have opened new possibilities for vaccine design and passive antibody therapies to provide sterilizing immunity and control blood-stage parasites.

Intra-host dynamics of co-infecting parasite genotypes in asymptomatic malaria patients

Malaria-infected individuals often harbor mixtures of genetically distinct parasite genotypes. We studied intra-host dynamics of parasite genotypes co-infecting asymptomatic adults in an area of intense malaria transmission in Chikhwawa, Malawi. Serial blood samples (5 ml) were collected over seven consecutive days from 25 adults with asymptomatic Plasmodium falciparum malaria and analyzed to determine whether a single peripheral blood sample accurately captures within-host parasite diversity. Blood samples from three of the participants were also analyzed by limiting dilution cloning and SNP genotyping of the parasite clones isolated to examine both the number and relatedness of co-infecting parasite haplotypes. We observed rapid turnover of co-infecting parasite genotypes in 88% of the individuals sampled (n = 22) such that the genetic composition of parasites infecting these individuals changed dramatically over the course of seven days of follow up. Nineteen of the 25 individuals sampled (76%) carried multiple parasite genotypes at baseline. Analysis of serial blood samples from three of the individuals revealed that they harbored 6, 12 and 17 distinct parasite haplotypes respectively. Approximately 70% of parasite haplotypes recovered from the three extensively sampled individuals were unrelated (proportion of shared alleles <83.3%) and were deemed to have primarily arisen from superinfection (inoculation of unrelated parasite haplotypes through multiple mosquito bites). The rest were related at the half-sib level or greater and were deemed to have been inoculated into individual human hosts via parasite co-transmission from single mosquito bites. These findings add further to the growing weight of evidence indicating that a single blood sample poorly captures within-host parasite diversity and underscore the importance of repeated blood sampling to accurately capture within-host parasite ecology. Our data also demonstrate a more pronounced role for parasite co-transmission in generating within-host parasite diversity in high transmission settings than previously assumed. Taken together, these findings have important implications for understanding the evolution of drug resistance, malaria transmission, parasite virulence, allocation of gametocyte sex ratios and acquisition of malaria immunity.

Networks of genetic similarity reveal non-neutral processes shape strain structure in Plasmodium falciparum

Pathogens compete for hosts through patterns of cross-protection conferred by immune responses to antigens. In Plasmodium falciparum malaria, the var multigene family encoding for the major blood-stage antigen PfEMP1 has evolved enormous genetic diversity through ectopic recombination and mutation. With 50-60 var genes per genome, it is unclear whether immune selection can act as a dominant force in structuring var repertoires of local populations. The combinatorial complexity of the var system remains beyond the reach of existing strain theory and previous evidence for non-random structure cannot demonstrate immune selection without comparison with neutral models. We develop two neutral models that encompass malaria epidemiology but exclude competitive interactions between parasites. These models, combined with networks of genetic similarity, reveal non-neutral strain structure in both simulated systems and an extensively sampled population in Ghana. The unique population structure we identify underlies the large transmission reservoir characteristic of highly endemic regions in Africa.

The persistence and oscillations of submicroscopic Plasmodium falciparum and Plasmodium vivax infections over time in Vietnam: an open cohort study

BACKGROUND

A substantial proportion of Plasmodium species infections are asymptomatic with densities too low to be detectable with standard diagnostic techniques. The importance of such asymptomatic plasmodium infections in malaria transmission is probably related to their duration and density. To explore the duration of asymptomatic plasmodium infections and changes in parasite densities over time, a cohort of participants who were infected with Plasmodium parasites was observed over a 2-year follow-up period.

METHODS

In this open cohort study, inhabitants of four villages in Vietnam were invited to participate in baseline and subsequent 3-monthly surveys up to 24 months, which included the collection of venous blood samples. Samples were batch-screened using ultra-sensitive (u)PCR (lower limit of detection of 22 parasites per mL). Participants found to be infected by uPCR during any of these surveys were invited to join a prospective cohort and provide monthly blood samples. We estimated the persistence of Plasmodium falciparum and Plasmodium vivax infections and changes in parasite densities over a study period of 24 months.

FINDINGS

Between Dec 1, 2013, and Jan 8, 2016, 356 villagers participated in between one and 22 surveys. These study participants underwent 4248 uPCR evaluations (11·9 tests per participant). 1874 (32%) of 4248 uPCR tests indicated a plasmodium infection; 679 (36%) of 1874 tests were P falciparum monoinfections, 507 (27%) were P vivax monoinfections, 463 (25%) were co-infections with P falciparum and P vivax, and 225 (12%) were indeterminate species of Plasmodium. The median duration of P falciparum infection was 2 months (IQR 1-3); after accounting for censoring, participants had a 20% chance of having parasitaemia for 4 months or longer. The median duration of P vivax infection was 6 months (3-9), and participants had a 59% chance of having parasitaemia for 4 months or longer. The parasite densities of persistent infections oscillated; following ultralow-density infections, high-density infections developed frequently.

INTERPRETATION

Persistent largely asymptomatic P vivax and P falciparum infections are common in this area of low seasonal malaria transmission. Infections with low-density parasitaemias can develop into much higher density infections at a later time, which are likely to sustain malaria endemicity.

FUNDING

The Wellcome Trust, Bill & Melinda Gates Foundation.

Modeling the genetic relatedness of Plasmodium falciparum parasites following meiotic recombination and cotransmission

Unlike in most pathogens, multiple-strain (polygenomic) infections of P. falciparum are frequently composed of genetic siblings. These genetic siblings are the result of sexual reproduction and can coinfect the same host when cotransmitted by the same mosquito. The degree with which coinfecting strains are related varies among infections and populations. Because sexual recombination occurs within the mosquito, the relatedness of cotransmitted strains could depend on transmission dynamics, but little is actually known of the factors that influence the relatedness of cotransmitted strains. Part of the uncertainty stems from an incomplete understanding of how within-host and within-vector dynamics affect cotransmission. Cotransmission is difficult to examine experimentally but can be explored using a computational model. We developed a malaria transmission model that simulates sexual reproduction in order to understand what determines the relatedness of cotransmitted strains. This study highlights how the relatedness of cotransmitted strains depends on both within-host and within-vector dynamics including the complexity of infection. We also used our transmission model to analyze the genetic relatedness of polygenomic infections following a series of multiple transmission events and examined the effects of superinfection. Understanding the factors that influence the relatedness of cotransmitted strains could lead to a better understanding of the population-genetic correlates of transmission and therefore be important for public health.

Genomic studies of P. falciparum reveal that multi-strain infections can include genetically related strains. P. falciparum must reproduce sexually in the mosquito vector. One consequence of sexual reproduction is that parasites cotransmitted by the same mosquito are related to one another. The degree of genetic relatedness of these parasites can be as great as that of full-siblings. However, our understanding of the cotransmission process is incomplete, and little is known of the role of cotransmission in influencing population genomic processes. To help bridge this gap, we developed a simulation model to determine which of the steps involved in transmission have the greatest impact on the relatedness of parasites cotransmitted by a mosquito vector. The primary goal of this study is to characterize the outcomes of cotransmission following single or multiple transmission events. Our model yields new insights into the cotransmission process, which we believe will be useful for understanding the results from more complicated population models and epidemiological conditions. Such an understanding is important for the use of population genomics to inform public health decisions as well as for understanding of parasite evolution.

Higher Complexity of Infection and Genetic Diversity of Plasmodium vivax Than Plasmodium falciparum Across All Malaria Transmission Zones of Papua New Guinea

Plasmodium falciparum and Plasmodium vivax have varying transmission dynamics that are informed by molecular epidemiology. This study aimed to determine the complexity of infection and genetic diversity of P. vivax and P. falciparum throughout Papua New Guinea (PNG) to evaluate transmission dynamics across the country. In 2008-2009, a nationwide malaria indicator survey collected 8,936 samples from all 16 endemic provinces of PNG. Of these, 892 positive P. vivax samples were genotyped at PvMS16 and PvmspF3, and 758 positive P. falciparum samples were genotyped at Pfmsp2. The data were analyzed for multiplicity of infection (MOI) and genetic diversity. Overall, P. vivax had higher polyclonality (71%) and mean MOI (2.32) than P. falciparum (20%, 1.39). These measures were significantly associated with prevalence for P. falciparum but not for P. vivax. The genetic diversity of P. vivax (PvMS16: expected heterozygosity = 0.95, 0.85-0.98; PvMsp1F3: 0.78, 0.66-0.89) was higher and less variable than that of P. falciparum (Pfmsp2: 0.89, 0.65-0.97). Significant associations of MOI with allelic richness (rho = 0.69, P = 0.009) and expected heterozygosity (rho = 0.87, P < 0.001) were observed for P. falciparum. Conversely, genetic diversity was not correlated with polyclonality nor mean MOI for P. vivax. The results demonstrate higher complexity of infection and genetic diversity of P. vivax across the country. Although P. falciparum shows a strong association of these parameters with prevalence, a lack of association was observed for P. vivax and is consistent with higher potential for outcrossing of this species.

Antibody-independent mechanisms regulate the establishment of chronic Plasmodium infection

Malaria is caused by parasites of the genus Plasmodium. All human-infecting Plasmodium species can establish long-lasting chronic infections^1-5, creating an infectious reservoir to sustain transmission^1,6. It is widely accepted that the maintenance of chronic infection involves evasion of adaptive immunity by antigenic variation⁷. However, genes involved in this process have been identified in only two of five human-infecting species: Plasmodium falciparum and Plasmodium knowlesi. Furthermore, little is understood about the early events in the establishment of chronic infection in these species. Using a rodent model we demonstrate that from the infecting population, only a minority of parasites, expressing one of several clusters of virulence-associated pir genes, establishes a chronic infection. This process occurs in different species of parasites and in different hosts. Establishment of chronicity is independent of adaptive immunity and therefore different from the mechanism proposed for maintenance of chronic P. falciparum infections^7-9. Furthermore, we show that the proportions of parasites expressing different types of pir genes regulate the time taken to establish a chronic infection. Because pir genes are common to most, if not all, species of Plasmodium¹⁰, this process may be a common way of regulating the establishment of chronic infections.

Numerical Distributions of Parasite Densities During Asymptomatic Malaria

BACKGROUND

Asymptomatic parasitemia is common even in areas of low seasonal malaria transmission, but the true proportion of the population infected has not been estimated previously because of the limited sensitivity of available detection methods.

METHODS

Cross-sectional malaria surveys were conducted in areas of low seasonal transmission along the border between eastern Myanmar and northwestern Thailand and in western Cambodia. DNA was quantitated by an ultrasensitive polymerase chain reaction (uPCR) assay (limit of accurate detection, 22 parasites/mL) to characterize parasite density distributions for Plasmodium falciparum and Plasmodium vivax, and the proportions of undetected infections were imputed.

RESULTS

The prevalence of asymptomatic malaria as determined by uPCR was 27.5% (1303 of 4740 people tested). Both P. vivax and P. falciparum density distributions were unimodal and log normal, with modal values well within the quantifiable range. The estimated proportions of all parasitemic individuals identified by uPCR were >70% among individuals infected with P. falciparum and >85% among those infected with P. vivax. Overall, 83% of infections were predicted to be P. vivax infections, 13% were predicted to be P. falciparum infections, and 4% were predicted to be mixed infections. Geometric mean parasite densities were similar; 5601 P. vivax parasites/mL and 5158 P. falciparum parasites/mL.

CONCLUSIONS

This uPCR method identified most infected individuals in malaria-endemic areas. Malaria parasitemia persists in humans at levels that optimize the probability of generating transmissible gametocyte densities without causing illness.

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.

Molecular genetic analysis of Plasmodium vivax isolates from Eastern and Central Sudan using pvcsp and pvmsp-3α genes as molecular markers

In Sudan, Plasmodium vivax accounts for approximately 5-10% of malaria cases. This study was carried out to determine the genetic diversity of P. vivax population from Sudan by analyzing the polymorphism of P. vivax csp (pvcsp) and pvmsp-3α genes. Blood samples (n=76) were taken from suspected malaria cases from 2012-2013 in three health centers of Eastern and Central Sudan. Parasite detection was performed by microscopy and molecular techniques, and genotyping of both genes was performed by PCR-RFLP followed by DNA sequence for only pvcsp gene (n=30). Based on microscopy analysis, 76 (%100) patients were infected with P. vivax, whereas nested-PCR results showed that 86.8% (n=66), 3.9% (n=3), and 3.9% (n=3) of tested samples had P. vivax as well as Plasmodium falciparum mono- and mixed infections, respectively. Four out of 76 samples had no results in molecular diagnosis. All sequenced samples were found to be of VK210 (100%) genotype with six distinct amino acid haplotypes, and 210A (66.7%) was the most prevalent haplotype. The Sudanese isolates displayed variations in the peptide repeat motifs (PRMs) ranging from 17 to 19 with GDRADGQPA (PRM1), GDRAAGQPA (PRM2) and DDRAAGQPA (PRM3). Also, 54 polymorphic sites with 56 mutations were found in repeat and post-repeat regions of the pvcsp and the overall nucleotide diversity (π) was 0.02149±0.00539. A negative value of dN-dS (-0.0344) was found that suggested a significant purifying selection of Sudanese pvcsp, (Z test, P<0.05). Regarding pvmsp-3α, three types were detected: types A (94.6%, 52/55), type C (3.6%, 2/55), and type B (1.8%, 1/55). No multiclonal infections were detected, and RFLP analysis identified 13 (Hha I, A1-A11, B1, and C1) and 16 (Alu I, A1-A14, B1, and C1) distinct allelic forms. In conclusion, genetic investigation among Sudanese P. vivax isolates indicated that this antigen showed limited antigenic diversity.

Within-host competition does not select for virulence in malaria parasites; studies with Plasmodium yoelii

In endemic areas with high transmission intensities, malaria infections are very often composed of multiple genetically distinct strains of malaria parasites. It has been hypothesised that this leads to intra-host competition, in which parasite strains compete for resources such as space and nutrients. This competition may have repercussions for the host, the parasite, and the vector in terms of disease severity, vector fitness, and parasite transmission potential and fitness. It has also been argued that within-host competition could lead to selection for more virulent parasites. Here we use the rodent malaria parasite Plasmodium yoelii to assess the consequences of mixed strain infections on disease severity and parasite fitness. Three isogenic strains with dramatically different growth rates (and hence virulence) were maintained in mice in single infections or in mixed strain infections with a genetically distinct strain. We compared the virulence (defined as harm to the mammalian host) of mixed strain infections with that of single infections, and assessed whether competition impacted on parasite fitness, assessed by transmission potential. We found that mixed infections were associated with a higher degree of disease severity and a prolonged infection time. In the mixed infections, the strain with the slower growth rate was often responsible for the competitive exclusion of the faster growing strain, presumably through host immune-mediated mechanisms. Importantly, and in contrast to previous work conducted with Plasmodium chabaudi, we found no correlation between parasite virulence and transmission potential to mosquitoes, suggesting that within-host competition would not drive the evolution of parasite virulence in P. yoelii.

The duration of Plasmodium falciparum infections

Plasmodium vivax and Plasmodium ovale are often considered the malaria parasites best adapted to long-term survival in the human host because of their latent exo-erythrocytic forms. The prevailing opinion until the middle of the last century was that the maximum duration of Plasmodium falciparum infections was less than two years. Case reports and series investigating blood donors following accidental malaria infection of blood transfusion recipients and other sporadic malaria cases in non-endemic countries have shown clearly that asymptomatic P. falciparum infections may persist for up to a decade or longer (maximum confirmed 13 years). Current policies in malaria-free countries of excluding blood donors who have lived in malarious areas are justified. Vigilance for longer than three years after declaring elimination in an area may be needed.

Differential PfEMP1 expression is associated with cerebral malaria pathology

Plasmodium falciparum is unique among human malarias in its ability to sequester in post-capillary venules of host organs. The main variant antigens implicated are the P. falciparum erythrocyte membrane protein 1 (PfEMP1), which can be divided into three major groups (A–C). Our study was a unique examination of sequestered populations of parasites for genetic background and expression of PfEMP1 groups. We collected post-mortem tissue from twenty paediatric hosts with pathologically different forms of cerebral malaria (CM1 and CM2) and parasitaemic controls (PC) to directly examine sequestered populations of parasites in the brain, heart and gut. Use of two different techniques to investigate this question produced divergent results. By quantitative PCR, group A var genes were upregulated in all three organs of CM2 and PC cases. In contrast, in CM1 infections displaying high levels of sequestration but negligible vascular pathology, there was high expression of group B var. Cloning and sequencing of var transcript tags from the same samples indicated a uniformly low expression of group A-like var. Generally, within an organ sample, 1–2 sequences were expressed at dominant levels. 23% of var tags were detected in multiple patients despite the P. falciparum infections being genetically distinct, and two tags were observed in up to seven hosts each with high expression in the brains of 3–4 patients. This study is a novel examination of the sequestered parasites responsible for fatal cerebral malaria and describes expression patterns of the major cytoadherence ligand in three organ-derived populations and three pathological states.

One of the most severe forms of malarial disease is cerebral malaria, which disproportionally affects young children. In this disease, the parasite places proteins on the red blood cell surface, providing a “smokescreen” by which they evade host immunity and hide in organ blood vessels, blocking them and causing tissue damage. It is impossible to study parasites in the organs during life and autopsy studies on children with malaria are exceedingly rare. In Malawi, we examined parasites from the brain, heart and intestine of twenty cases of fatal malaria including controls with low numbers of malaria parasites but another identified cause of death. We found little difference in the category of proteins the parasites used in controls and cerebral malaria, although a small number of specific proteins were detected in multiple infections. In an alternative form of malaria in which the brain is heavily infected but shows no evidence of damage, we found a different set of proteins at high proportion. However, as these children were typically older and most were infected with HIV, we could not determine which of these factors was most important. Interactions between host and parasite have the potential to influence disease outcomes.

The Plasmodium vivax merozoite surface protein 3β sequence reveals contrasting parasite populations in southern and northwestern Thailand

BACKGROUND

Malaria control efforts have a significant impact on the epidemiology and parasite population dynamics. In countries aiming for malaria elimination, malaria transmission may be restricted to limited transmission hot spots, where parasite populations may be isolated from each other and experience different selection forces. Here we aim to examine the Plasmodium vivax population divergence in geographically isolated transmission zones in Thailand.

METHODOLOGY

We employed the P. vivax merozoite surface protein 3β (PvMSP3β) as a molecular marker for characterizing P. vivax populations based on the extensive diversity of this gene in Southeast Asian parasite populations. To examine two parasite populations with different transmission levels in Thailand, we obtained 45 P. vivax isolates from Tak Province, northwestern Thailand, where the annual parasite incidence (API) was more than 2%, and 28 isolates from Yala and Narathiwat Provinces, southern Thailand, where the API was less than 0.02%. We sequenced the PvMSP3β gene and examined its genetic diversity and molecular evolution between the parasite populations.

PRINCIPAL FINDINGS

Of 58 isolates containing single PvMSP3β alleles, 31 sequence types were identified. The overall haplotype diversity was 0.77 ± 0.06 and nucleotide diversity 0.0877±0.0054. The northwestern vivax malaria population exhibited extensive haplotype diversity (HD) of PvMSP3β (HD=1.0). In contrast, the southern parasite population displayed a single PvMSP3β allele (HD=0), suggesting a clonal population expansion. This result revealed that the extent of allelic diversity in P. vivax populations in Thailand varies among endemic areas.

CONCLUSION

Malaria parasite populations in a given region may vary significantly in genetic diversity, which may be the result of control and influenced by the magnitude of malaria transmission intensity. This is an issue that should be taken into account for the implementation of P. vivax control measures such as drug policy and vaccine development.

Understanding the relationship between Plasmodium falciparum growth rate and multiplicity of infection

Natural infections with Plasmodium falciparum are often composed of multiple concurrent genetically distinct parasite clones. Such multiclonal infections are more common in areas of high transmission, and the frequency of multiclonal infection also varies with age. A number of studies have suggested that multiclonal infection predicts the risk of subsequent clinical malaria. The multiplicity of infection is determined by the rate of new infections, the number of clones inoculated at each mosquito bite, and the duration of infections. Here, we used a mathematical modeling approach to understand how variation in the growth rate of blood-stage parasites affects the observed multiplicity of infection (MOI), as well as the relationship between the MOI and the risk of subsequent malaria. We then analyzed data from a study of multiclonal infection and malaria in an malaria-endemic area in Tanzania and show that the proportion of multiclonal infections varies with age and that the observed relationship between multiclonal infection and subsequent clinical events can be explained by a reduction in blood-stage parasite growth with age in this population.

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.

Strategies for designing and monitoring malaria vaccines targeting diverse antigens

After more than 50 years of intensive research and development, only one malaria vaccine candidate, “RTS,S,” has progressed to Phase 3 clinical trials. Despite only partial efficacy, this candidate is now forecast to become the first licensed malaria vaccine. Hence, more efficacious second-generation malaria vaccines that can significantly reduce transmission are urgently needed. This review will focus on a major obstacle hindering development of effective malaria vaccines: parasite antigenic diversity. Despite extensive genetic diversity in leading candidate antigens, vaccines have been and continue to be formulated using recombinant antigens representing only one or two strains. These vaccine strains represent only a small fraction of the diversity circulating in natural parasite populations, leading to escape of non-vaccine strains and challenging investigators’ abilities to measure strain-specific efficacy in vaccine trials. Novel strategies are needed to overcome antigenic diversity in order for vaccine development to succeed. Many studies have now cataloged the global diversity of leading Plasmodium falciparum and Plasmodium vivax vaccine antigens. In this review, we describe how population genetic approaches can be applied to this rich data source to predict the alleles that best represent antigenic diversity, polymorphisms that contribute to it, and to identify key polymorphisms associated with antigenic escape. We also suggest an approach to summarize the known global diversity of a given antigen to predict antigenic diversity, how to select variants that best represent the strains circulating in natural parasite populations and how to investigate the strain-specific efficacy of vaccine trials. Use of these strategies in the design and monitoring of vaccine trials will not only shed light on the contribution of genetic diversity to the antigenic diversity of malaria, but will also maximize the potential of future malaria vaccine candidates.

Distribution, diversity and drivers of blood-borne parasite co-infections in Alaskan bird populations

Avian species are commonly infected by multiple parasites, however few studies have investigated the environmental determinants of the prevalence of co-infection over a large scale. Here we believe that we report the first, detailed ecological study of the prevalence, diversity and co-infections of four avian blood-borne parasite genera: Plasmodium spp., Haemoproteus spp., Leucocytozoon spp. and Trypanosoma spp. We collected blood samples from 47 resident and migratory bird species across a latitudinal gradient in Alaska. From the patterns observed at collection sites, random forest models were used to provide evidence of associations between bioclimatic conditions and the prevalence of parasite co-infection distribution. Molecular screening revealed a higher prevalence of haematozoa (53%) in Alaska than previously reported. Leucocytozoons had the highest diversity, prevalence and prevalence of co-infection. Leucocytozoon prevalence (35%) positively correlated with Trypanosoma prevalence (11%), negatively correlated with Haemoproteus prevalence (14%) and had no correlation with Plasmodium prevalence (7%). We found temperature, precipitation and tree cover to be the primary environmental drivers that show a relationship with the prevalence of co-infection. The results provide insight into the impacts of bioclimatic drivers on parasite ecology and intra-host interactions, and have implications for the study of infectious diseases in rapidly changing environments.

Genetic shifting: a novel approach for controlling vector-borne diseases

Rendering populations of vectors of diseases incapable of transmitting pathogens through genetic methods has long been a goal of vector geneticists. We outline a method to achieve this goal that does not involve the introduction of any new genetic variants to the target population. Rather we propose that shifting the frequencies of naturally occurring alleles that confer refractoriness to transmission can reduce transmission below a sustainable level. The program employs methods successfully used in plant and animal breeding. Because no artificially constructed genetically modified organisms (GMOs) are introduced into the environment, the method is minimally controversial. We use Aedes aegypti and dengue virus (DENV) for illustrative purposes but point out that the proposed program is generally applicable to vector-borne disease control.

Polymorphic patterns of the merozoite surface protein-3β in Korean isolates of Plasmodium vivax

BACKGROUND

The merozoite surface protein-3β of Plasmodium vivax (PvMSP-3β) is one of the candidate antigens for blood stage malaria vaccine development. The polymorphisms in PvMSP-3β have been reported in certain P. vivax isolates. However, the diversity of PvMSP-3β throughout its global distribution has not been well understood. In this study, the genetic diversity and the effects of natural selection in PvMSP-3β among P. vivax Korean isolates were analysed.

METHODS

Blood samples were collected from 95 patients with vivax malaria in Korea. The region flanking full-length PvMSP-3β was amplified by polymerase chain reaction and cloned into a TA cloning vector. The PvMSP-3β sequence of each isolate was determined and the polymorphic characteristics and effects of natural selection were analysed using the DNASTAR, MEGA4, and DnaSP programs.

RESULTS

Five different subtypes of PvMSP-3β were identified based on single nucleotide polymorphisms (SNPs), insertions, and deletions. Although a high level of sequence diversity was observed in the PvMSP-3β gene, the coiled-coil tertiary structure of the PvMSP-3β protein was well conserved in all of the sequences. The PvMSP-3β of Korean isolates is under natural selection. DNA polymerase slippage and intragenic recombination likely contributed to PvMSP-3β diversity in Korean P. vivax isolates.

CONCLUSIONS

The PvMSP-3β of Korean P. vivax isolates displayed polymorphisms, with SNPs, insertions and deletions scattered throughout of the gene. These results of parasite heterogeneity are relevant to the development of a PvMSP-3β based vaccine against P. vivax and the implementation of malaria control programmes in Korea.

Genotyping Plasmodium vivax isolates from the 2011 outbreak in Greece

BACKGROUND

Plasmodium vivax malaria was common in Greece until the 1950s with epidemics involving thousands of cases every year. Greece was declared free of malaria by the World Health Organization in 1974. From 1974 to 2010, an average of 39 cases per year were reported, which were mainly imported. However, in 2009 and 2010 six and one autochthonous cases were reported culminating with a total of 40 autochthonous cases reported in 2011, of which 34 originated from a single region: Laconia of Southern Peloponnese. In this study the genotypic complexity of the P. vivax infections from the outbreak in Greece during 2011 is described, to elucidate the possible origin and spread of the disease.

METHODS

Three polymorphic markers of P. vivax were used; Pvmsp-3α and the microsatellites m1501 and m3502 on P. vivax isolates sampled from individuals diagnosed in Greece. Thirty-nine isolates were available for this study (20 autochthonous and 19 imported), mostly from Evrotas municipality in Laconia region, in southern Greece, (n = 29), with the remaining representing sporadic cases originating from other areas of Greece.

RESULTS

Genotyping the Evrotas samples revealed seven different haplotypes where the majority of the P. vivax infections expressed two particular Pvmsp-3α-m1501-m3502 haplotypes, A10-128-151 (n = 14) and A10-121-142 (n = 7). These haplotypes appeared throughout the period in autochthonous and imported cases, indicating continuous transmission. In contrast, the P. vivax autochthonous cases from other parts of Greece were largely comprised of unique haplotypes, indicating limited transmission in these other areas.

CONCLUSIONS

The results indicate that several P. vivax strains were imported into various areas of Greece in 2011, thereby increasing the risk of re-introduction of malaria. In the region of Evrotas ongoing transmission occurred exemplifying that further control measures are urgently needed in this region of southern Europe. In circumstances where medical or travel history is scarce, methods of molecular epidemiology may prove highly useful for the correct classification of the cases.

Epistatic Interactions between apolipoprotein E and hemoglobin S Genes in regulation of malaria parasitemia

Apolipoprotein E is a monomeric protein secreted by the liver and responsible for the transport of plasma cholesterol and triglycerides. The APOE gene encodes 3 isoforms Ɛ4, Ɛ3 and Ɛ2 with APOE Ɛ4 associated with higher plasma cholesterol levels and increased pathogenesis in several infectious diseases (HIV, HSV). Given that cholesterol is an important nutrient for malaria parasites, we examined whether APOE Ɛ4 was a risk factor for Plasmodium infection, in terms of prevalence or parasite density. A cross sectional survey was performed in 508 children aged 1 to 12 years in Gabon during the wet season. Children were screened for Plasmodium spp. infection, APOE and hemoglobin S (HbS) polymorphisms. Median parasite densities were significantly higher in APOE Ɛ4 children for Plasmodium spp. densities compared to non-APOE Ɛ4 children. When stratified for HbS polymorphisms, median Plasmodium spp. densities were significantly higher in HbAA children if they had an APOE Ɛ4 allele compared to those without an APOE Ɛ4 allele. When considering non-APOE Ɛ4 children, there was no quantitative reduction of Plasmodium spp. parasite densities for HbAS compared to HbAA phenotypes. No influence of APOE Ɛ4 on successful Plasmodium liver cell invasion was detected by multiplicity of infection. These results show that the APOE Ɛ4 allele is associated with higher median malaria parasite densities in children likely due to the importance of cholesterol availability to parasite growth and replication. Results suggest an epistatic interaction between APOE and HbS genes such that sickle cell trait only had an effect on parasite density in APOE Ɛ4 children. This suggests a linked pathway of regulation of parasite density involving expression of these genes. These findings have significance for understanding host determinants of regulation of malaria parasite density, the design of clinical trials as well as studies of co-infection with Plasmodium and other pathogens.

Merozoite surface protein-3 alpha as a genetic marker for epidemiologic studies in Plasmodium vivax: a cautionary note

Background

Plasmodium vivax is the most widespread of the human malaria parasites in terms of geography, and is thought to present unique challenges to local efforts aimed at control and elimination. Parasite molecular markers can provide much needed data on P. vivax populations, but few such markers have been critically evaluated. One marker that has seen extensive use is the gene encoding merozoite surface protein 3-alpha (MSP-3α), a blood-stage antigen known to be highly variable among P. vivax isolates. Here, a sample of complete msp-3α gene sequences is analysed in order to assess its utility as a molecular marker for epidemiologic investigations.

Methods

Amplification, cloning and sequencing of additional P. vivax isolates from different geographic locations, including a set of Venezuelan field isolates (n = 10), yielded a sample of 48 complete msp-3α coding sequences. Characterization of standard population genetic measures of diversity, phylogenetic analysis, and tests for recombination were performed. This allowed comparisons to patterns inferred from the in silico simulation of a polymerase chain reaction restriction fragment length polymorphism (PCR-RFLP) protocol used widely.

Results

The larger sample of MSP-3α diversity revealed incongruence between the observed levels of nucleotide polymorphism, which were high in all populations, and the pattern of PCR-RFLP haplotype diversity. Indeed, PCR-RFLP haplotypes were not informative of a population’s genetic diversity and identical haplotypes could be produced from analogous bands in the commonly used protocol. Evidence of frequent and variable insertion-deletion mutations and recurrent recombination between MSP-3α haplotypes complicated the inference of genetic diversity patterns and reduced the phylogenetic signal.

Conclusions

The genetic diversity of P. vivax msp-3α involves intragenic recombination events. Whereas the high genetic diversity of msp-3α makes it a promising marker for some epidemiological applications, the ability of msp-3α PCR-RFLP analysis to accurately track parasites is limited. Local studies of the circulating alleles are needed before implementing PCR-RFLP approaches. Furthermore, evidence from the global sample analysed here suggests such msp-3α PCR-RFLP methods are not suitable for broad geographic studies or tracking parasite populations for an extended period of time.

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.

High levels of genetic diversity of Plasmodium falciparum populations in Papua New Guinea despite variable infection prevalence

High levels of genetic diversity in Plasmodium falciparum populations are an obstacle to malaria control. Here, we investigate the relationship between local variation in malaria epidemiology and parasite genetic diversity in Papua New Guinea (PNG). Cross-sectional malaria surveys were performed in 14 villages spanning four distinct malaria-endemic areas on the north coast, including one area that was sampled during the dry season. High-resolution msp2 genotyping of 2,147 blood samples identified 761 P. falciparum infections containing a total of 1,392 clones whose genotypes were used to measure genetic diversity. Considerable variability in infection prevalence and mean multiplicity of infection was observed at all of the study sites, with the area sampled during the dry season showing particularly striking local variability. Genetic diversity was strongly associated with multiplicity of infection but not with infection prevalence. In highly endemic areas, differences in infection prevalence may not translate into a decrease in parasite population diversity.

The dynamics of natural Plasmodium falciparum infections

BACKGROUND

Natural immunity to Plasmodium falciparum has been widely studied, but its effects on parasite dynamics are poorly understood. Acquisition and clearance rates of untreated infections are key elements of the dynamics of malaria, but estimating these parameters is challenging because of frequent super-infection and imperfect detectability of parasites. Consequently, information on effects of host immune status or age on infection dynamics is fragmentary.

METHODS

An age-stratified cohort of 347 individuals from Northern Ghana was sampled six times at 2 month intervals. High-throughput capillary electrophoresis was used to genotype the msp-2 locus of all P. falciparum infections detected by PCR. Force of infection (FOI) and duration were estimated for each age group using an immigration-death model that allows for imperfect detection of circulating parasites.

RESULTS

Allowing for imperfect detection substantially increased estimates of FOI and duration. Effects of naturally acquired immunity on the FOI and duration would be reflected in age dependence in these indices, but in our cohort data FOI tended to increase with age in children. Persistence of individual parasite clones was characteristic of all age-groups. Duration peaked in 5-9 year old children (average duration 319 days, 95% confidence interval 318;320).

CONCLUSIONS

The main age-dependence is on parasite densities, with only small age-variations in the FOI and persistence of infections. This supports the hypothesis that acquired immunity controls transmission mainly by limiting blood-stage parasite densities rather than changing rates of acquisition or clearance of infections.

Estimating the numbers of malaria infections in blood samples using high-resolution genotyping data

People living in endemic areas often habour several malaria infections at once. High-resolution genotyping can distinguish between infections by detecting the presence of different alleles at a polymorphic locus. However the number of infections may not be accurately counted since parasites from multiple infections may carry the same allele. We use simulation to determine the circumstances under which the number of observed genotypes are likely to be substantially less than the number of infections present and investigate the performance of two methods for estimating the numbers of infections from high-resolution genotyping data. The simulations suggest that the problem is not substantial in most datasets: the disparity between the mean numbers of infections and of observed genotypes was small when there was 20 or more alleles, 20 or more blood samples, a mean number of infections of 6 or less and where the frequency of the most common allele was no greater than 20%. The issue of multiple infections carrying the same allele is unlikely to be a major component of the errors in PCR-based genotyping. Simulations also showed that, with heterogeneity in allele frequencies, the observed frequencies are not a good approximation of the true allele frequencies. The first method that we proposed to estimate the numbers of infections assumes that they are a good approximation and hence did poorly in the presence of heterogeneity. In contrast, the second method by Li et al estimates both the numbers of infections and the true allele frequencies simultaneously and produced accurate estimates of the mean number of infections.

Broadly reactive antibodies specific for Plasmodium falciparum MSP-119 are associated with the protection of naturally exposed children against infection

Background

The 19 kDa C-terminal region of Plasmodium falciparum Merozoite Surface Protein-1 is a known target of naturally acquired humoral immunity and a malaria vaccine candidate. MSP-1₁₉ has four predominant haplotypes resulting in amino acid changes labelled EKNG, QKNG, QTSR and ETSR. IgG antibodies directed against all four variants have been detected, but it is not known if these variant specific antibodies are associated with haplotype-specific protection from infection.

Methods

Blood samples from 201 healthy Kenyan adults and children who participated in a 12-week treatment time-to-infection study were evaluated. Venous blood drawn at baseline (week 0) was examined for functional and serologic antibodies to MSP-1₁₉ and MSP-1₄₂ variants. MSP-1₁₉ haplotypes were detected by a multiplex PCR assay at baseline and weekly throughout the study. Generalized linear models controlling for age, baseline MSP-1₁₉ haplotype and parasite density were used to determine the relationship between infecting P. falciparum MSP-1₁₉ haplotype and variant-specific antibodies.

Results

A total of 964 infections resulting in 1,533 MSP-1₁₉ haplotypes detected were examined. The most common haplotypes were EKNG and QKNG, followed by ETSR and QTSR. Children had higher parasite densities, greater complexity of infection (>1 haplotype), and more frequent changes in haplotypes over time compared to adults. Infecting MSP-1₁₉ haplotype at baseline (week 0) had no influence on haplotypes detected over the subsequent 11 weeks among children or adults. Children but not adults with MSP-1₁₉ and some MSP-1₄₂ variant antibodies detected by serology at baseline had delayed time-to-infection. There was no significant association of variant-specific serology or functional antibodies at baseline with infecting haplotype at baseline or during 11 weeks of follow up among children or adults.

Conclusions

Variant transcending IgG antibodies to MSP-1₁₉ are associated with protection from infection in children, but not adults. These data suggest that inclusion of more than one MSP-1₁₉ variant may not be required in a malaria blood stage vaccine.

High degree of Plasmodium vivax diversity in the Peruvian Amazon demonstrated by tandem repeat polymorphism analysis

Molecular tools to distinguish strains of Plasmodium vivax are important for studying the epidemiology of malaria transmission. Two sets of markers-tandem repeat (TR) polymorphisms and MSP3α-were used to study Plasmodium vivax in patients in the Peruvian Amazon region of Iquitos. Of 110 patients, 90 distinct haplotypes were distinguished using 9 TR markers. An MSP3α polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) using HhaI and AluI revealed 8 and 9 profiles, respectively, and 36 profiles when analyzed in combination. Combining TR and PCR-RFLP markers, 101 distinct molecular profiles were distinguished among these 110 patients. Nine TR markers arrayed along a 100 kB stretch of a P. vivax chromosome containing the gene for circumsporozoite protein showed non-linear linkage disequilibrium (I(SA) = 0.03, P = 0.001). These findings demonstrate the potential use of TR markers for molecular epidemiology studies.

Close kinship within multiple-genotype malaria parasite infections

Malaria infections containing multiple parasite genotypes are ubiquitous in nature, and play a central role in models of recombination, intra-host dynamics, virulence, sex ratio, immunity and drug resistance evolution in Plasmodium. While these multiple infections (MIs) are often assumed to result from superinfection (bites from multiple infected mosquitoes), we know remarkably little about their composition or generation. We isolated 336 parasite clones from eight patients from Malawi (high transmission) and six from Thailand (low transmission) by dilution cloning. These were genotyped using 384 single-nucleotide polymorphisms, revealing 22 independent haplotypes in Malawi (2–6 per MI) and 15 in Thailand (2–5 per MI). Surprisingly, all six patients from Thailand and six of eight from Malawi contained related haplotypes, and haplotypes were more similar within- than between-infections. These results argue against a simple superinfection model. Instead, the observed kinship patterns may be explained by inoculation of multiple related haploid sporozoites from single mosquito bites, by immune suppression of parasite subpopulations within infections, and serial transmission of related parasites between people. That relatedness is maintained in endemic areas in the face of repeated bites from infected mosquitoes has profound implications for understanding malaria transmission, immunity and intra-host dynamics of co-infecting parasite genotypes.

Quantification of Plasmodium falciparum malaria from complex infections in the Peruvian Amazon using quantitative PCR of the merozoite surface protein 1, block 2 (PfMSP1-B2): in vitro dynamics reveal density-dependent interactions

The majority of Plasmodium falciparum field isolates are defined as complex infections because they contain multiple genetically distinct clones. Studying interactions between clones in complex infections in vivo and in vitro could elucidate important phenomena in malaria infection, transmission and treatment. Using quantitative PCR (qPCR) of the P. falciparum merozoite surface protein 1, block 2 (PfMSP1-B2), we provide a sensitive and efficient genotyping method. This is important for epidemiological studies because it makes it possible to study genotype-specific growth dynamics. We compared 3 PfMSP1-B2 genotyping methods by analysing 79 field isolates from the Peruvian Amazon. In vivo observations from other studies using these techniques led to the hypothesis that clones within complex infections interact. By co-culturing clones with different PfMSP1-B2 genotypes, and measuring parasitaemia using qPCR, we found that suppression of clonal expansion was a factor of the collective density of all clones present in a culture. PfMSP1-B2 qPCR enabled us to find in vitro evidence for parasite-parasite interactions and could facilitate future investigations of growth trends in naturally occurring complex infections.

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.

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.

Genetic diversity in the merozoite surface protein 1 and 2 genes of Plasmodium falciparum from the Artibonite Valley of Haiti

Describing genetic diversity of the Plasmodium falciparum parasite provides important information about the local epidemiology of malaria. In this study, we examined the genetic diversity of P. falciparum isolates from the Artibonite Valley in Haiti using the allelic families of merozoite surface protein 1 and 2 genes (msp-1 and msp-2). The majority of study subjects infected with P. falciparum had a single parasite genotype (56% for msp-1 and 69% for msp-2: n=79); 9 distinct msp-1 genotypes were identified by size differences on agarose gels. K1 was the most polymorphic allelic family with 5 genotypes (amplicons from 100 to 300 base pairs [bp]); RO33 was the least polymorphic, with a single genotype (120-bp). Although both msp-2 alleles (3D7/IC1, FC27) had similar number of genotypes (n=4), 3D7/IC1 was more frequent (85% vs. 26%). All samples were screened for the presence of the K76T mutation on the P. falciparum chloroquine resistance transporter (pfcrt) gene with 10 of 79 samples positive. Of the 2 (out of 10) samples from individuals follow-up for 21 days, P. falciparum parasites were present through day 7 after treatment with chloroquine. No parasites were found on day 21. Our results suggest that the level of genetic diversity is low in this area of Haiti, which is consistent with an area of low transmission.