Complexity of Infection and Genetic Diversity in Cambodian Plasmodium vivax

Single-cell analyses of polyclonal Plasmodium vivax infections and their consequences on parasite transmission

Most Plasmodium vivax infections contain genetically distinct parasites, but the consequences of this polyclonality on the development of asexual parasites, their sexual differentiation, and their transmission remain unknown. We describe infections of Saimiri monkeys with two strains of P. vivax and the analyses of 117,350 parasites characterized by single cell RNA sequencing and individually genotyped. In our model, consecutive inoculations fail to establish polyclonal infections. By contrast, simultaneous inoculations of two strains lead to sustained polyclonal infections, although without detectable differences in parasite regulation or sexual commitment. Analyses of sporozoites dissected from mosquitoes fed on coinfected monkeys show that all genotypes are successfully transmitted to mosquitoes. However, after sporozoite inoculation, not all genotypes contribute to the subsequent blood infections, highlighting an important bottleneck during pre-erythrocytic development. Overall, these studies provide new insights on the mechanisms regulating the establishment of polyclonal P. vivax infections and their consequences for disease transmission.

Genomic complexity of parasites and vectors challenges malaria control in Southeast Asia

Southeast Asia is a uniquely complex region of malaria transmission that maintains an astounding level of species diversity among potential malaria vectors and also generates drug-resistant and quickly diverging populations of malaria parasites. All five human malaria species circulate in Southeast Asia with over 50 Anopheles species that vary in their ability to transmit these pathogens. The intricate relationships of these parasites and vectors are not well-understood. Human activity in Southeast Asian countries has created an increasingly fragmented landscape, bringing humans and mosquitoes into more frequent contact, sustaining malaria transmission in a region where few control tools are effective. Genomic shifts at the species, population, and individual level in parasites and vectors introduce variation that has produced drug- and insecticide resistance. The goal of this review is to highlight genomic studies of Southeast Asian malaria parasites and vectors that demonstrate how diversity in these organisms presents unique challenges and opportunities for global malaria control and eradication efforts.

Genetic diversity of Plasmodium vivax isolates from pregnant women in the Western Brazilian Amazon: a prospective cohort study

Background

Each year, 92 million pregnant women are at risk of contracting malaria during pregnancy, with the underestimation of the mortality and morbidity burden associated with Plasmodium vivax. During pregnancy, P. vivax infection is associated with low birth weight, maternal anaemia, premature delivery, and stillbirth. In the State of Acre (Brazil), high transmission leaves pregnant women at greater risk of contracting malaria and having a greater number of recurrences. The study of genetic diversity and the association of haplotypes with adverse pregnancy effects is of great importance for the control of the disease. Here we investigate the genetic diversity of P. vivax parasites infecting pregnant women across their pregnancies.

Methods

P. vivax DNA was extracted from 330 samples from 177 women followed during pregnancy, collected in the State of Acre, Brazil. All samples were negative for Plasmodium falciparum DNA. Sequence data for the Pvmsp1 gene was analysed alongside data from six microsatellite (MS) markers. Allelic frequencies, haplotype frequencies, expected heterozygosity (HE) were calculated. Whole genome sequencing (WGS) was conducted on four samples from pregnant women and phylogenetic analysis performed with other samples from South American regions.

Findings

Initially, the pregnant women were stratified into two groups-1 recurrence and 2 or more recurrences-in which no differences were observed in clinical gestational outcomes or in placental histological changes between the two groups. Then we evaluated the parasites genetically. An average of 18.5 distinct alleles were found at each of the MS loci, and the HE calculated for each marker indicates a high genetic diversity occurring within the population. There was a high percentage of polyclonal infections (61.7%, 108/175), and one haplotype (H1) occurred frequently (20%), with only 9 of the haplotypes appearing in more than one patient.

Interpretation

Most pregnant women had polyclonal infections that could be the result of relapses and/or re-infections. The high percentage of H1 parasites, along with the low frequency of many other haplotypes are suggestive of a clonal expansion. Phylogenetic analysis shows that P. vivax population within pregnant women clustered with other Brazilian samples in the region.

Funding

FAPESP and CNPq - Brazil.

Molecular identification of vivax malaria relapse patients in the Yunnan Province based on homology analysis of the Plasmodium vivax circumsporozoite protein gene

More than 85% of the malaria burden in the Yunnan Province is caused by imported vivax malaria, and Yunnan is also where the majority of vivax malaria patients are diagnosed in China. Timely removal of the infection sources of Plasmodium vivax and its breeding environment remains the key to eliminating the secondary transmission of imported malaria. To that end, blood samples were collected from cases diagnosed and revalidated as single species infection with P. vivax in the Yunnan Province from 2013 to 2020. Specifically, samples from vivax malaria patients with suspected relapses episodes were subjected to PCR amplification, product sequencing, and analysis of the P. vivax circumsporozoite protein (pvcsp) gene. In total, 77 suspected relapse patients were identified out of 2484 cases infected with P. vivax, with a total of 81 recurrent episodes. A total of 156 CDS (coding DNA sequence) chains were obtained through PCR amplification and sequencing of the pvcsp gene from 159 blood samples, 121 of which can be matched to the paired sequences of 59 vivax malaria patients with both primary attack and recurrent experience. Of the 59 pairs of pvcsp gene sequences, every one of 31 pairs showed only one haplotype and no variant sites (VS), meaning every two paired sequence was completely homologous. Every one of the remaining 28 paired sequences had two haplotypes but no length polymorphism, indicating that the paired sequences was "weakly heterologous" with no fragment insertions (or deletions). All 59 vivax malaria patients with recurrences were caused by the activation of P. vivax hypnozoites originated from the same population as the primary infection. The paired analysis of the similarity between high variant genes allowed the identification of relapse episodes caused by P. vivax homologous hypnozoites and also demonstrated pvcsp gene as one of the candidate molecular markers for tracing infection origin.

Genomic approaches for monitoring transmission dynamics of malaria: A case for malaria molecular surveillance in Sub-Saharan Africa

Transmission dynamics is an important indicator for malaria control and elimination. As we move closer to eliminating malaria in Sub-Saharan Africa (sSA), transmission indices with higher resolution (genomic approaches) will complement our current measurements of transmission. Most of the present programmatic knowledge of malaria transmission patterns are derived from assessments of epidemiologic and clinical data, such as case counts, parasitological estimates of parasite prevalence, and Entomological Inoculation Rates (EIR). However, to eliminate malaria from endemic areas, we need to track changes in the parasite population and how they will impact transmission. This is made possible through the evolving field of genomics and genetics, as well as the development of tools for more in-depth studies on the diversity of parasites and the complexity of infections, among other topics. If malaria elimination is to be achieved globally, country-specific elimination activities should be supported by parasite genomic data from regularly collected blood samples for diagnosis, surveillance and possibly from other programmatic interventions. This presents a unique opportunity to track the spread of malaria parasites and shed additional light on intervention efficacy. In this review, various genetic techniques are highlighted along with their significance for an enhanced understanding of transmission patterns in distinct topological settings throughout Sub-Saharan Africa. The importance of these methods and their limitations in malaria surveillance to guide control and elimination strategies, are explored.

A maximum-likelihood method to estimate haplotype frequencies and prevalence alongside multiplicity of infection from SNP data

The introduction of genomic methods facilitated standardized molecular disease surveillance. For instance, SNP barcodes in Plasmodium vivax and Plasmodium falciparum malaria allows the characterization of haplotypes, their frequencies and prevalence to reveal temporal and spatial transmission patterns. A confounding factor is the presence of multiple genetically distinct pathogen variants within the same infection, known as multiplicity of infection (MOI). Disregarding ambiguous information, as usually done in ad-hoc approaches, leads to less confident and biased estimates. We introduce a statistical framework to obtain maximum-likelihood estimates (MLE) of haplotype frequencies and prevalence alongside MOI from malaria SNP data, i.e., multiple biallelic marker loci. The number of model parameters increases geometrically with the number of genetic markers considered and no closed-form solution exists for the MLE. Therefore, the MLE needs to be derived numerically. We use the Expectation-Maximization (EM) algorithm to derive the maximum-likelihood estimates, an efficient and easy-to-implement algorithm that yields a numerically stable solution. We also derive expressions for haplotype prevalence based on either all or just the unambiguous genetic information and compare both approaches. The latter corresponds to a biased ad-hoc estimate of prevalence. We assess the performance of our estimator by systematic numerical simulations assuming realistic sample sizes and various scenarios of transmission intensity. For reasonable sample sizes, and number of loci, the method has little bias. As an example, we apply the method to a dataset from Cameroon on sulfadoxine-pyrimethamine resistance in P. falciparum malaria. The method is not confined to malaria and can be applied to any infectious disease with similar transmission behavior. An easy-to-use implementation of the method as an R-script is provided.

Novel highly-multiplexed AmpliSeq targeted assay for Plasmodium vivax genetic surveillance use cases at multiple geographical scales

Although the power of genetic surveillance tools has been acknowledged widely, there is an urgent need in malaria endemic countries for feasible and cost-effective tools to implement in national malaria control programs (NMCPs) that can generate evidence to guide malaria control and elimination strategies, especially in the case of Plasmodium vivax. Several genetic surveillance applications (‘use cases’) have been identified to align research, technology development, and public health efforts, requiring different types of molecular markers. Here we present a new highly-multiplexed deep sequencing assay (Pv AmpliSeq). The assay targets the 33-SNP vivaxGEN-geo panel for country-level classification, and a newly designed 42-SNP within-country barcode for analysis of parasite dynamics in Vietnam and 11 putative drug resistance genes in a highly multiplexed NGS protocol with easy workflow, applicable for many different genetic surveillance use cases. The Pv AmpliSeq assay was validated using: 1) isolates from travelers and migrants in Belgium, and 2) routine collections of the national malaria control program at sentinel sites in Vietnam. The assay targets 229 amplicons and achieved a high depth of coverage (mean 595.7 ± 481) and high accuracy (mean error-rate of 0.013 ± 0.007). P. vivax parasites could be characterized from dried blood spots with a minimum of 5 parasites/µL and 10% of minority-clones. The assay achieved good spatial specificity for between-country prediction of origin using the 33-SNP vivaxGEN-geo panel that targets rare alleles specific for certain countries and regions. A high resolution for within-country diversity in Vietnam was achieved using the designed 42-SNP within-country barcode that targets common alleles (median MAF 0.34, range 0.01-0.49. Many variants were detected in (putative) drug resistance genes, with different predominant haplotypes in the pvmdr1 and pvcrt genes in different provinces in Vietnam. The capacity of the assay for high resolution identity-by-descent (IBD) analysis was demonstrated and identified a high rate of shared ancestry within Gia Lai Province in the Central Highlands of Vietnam, as well as between the coastal province of Binh Thuan and Lam Dong. Our approach performed well in geographically differentiating isolates at multiple spatial scales, detecting variants in putative resistance genes, and can be easily adjusted to suit the needs in other settings in a country or region. We prioritize making this tool available to researchers and NMCPs in endemic countries to increase ownership and ensure data usage for decision-making and malaria policy.

Dynamic Classification of Plasmodium vivax Malaria Recurrence: An Application of Classifying Unknown Cause of Failure in Competing Risks

A standard competing risks set-up requires both time to event and cause of failure to be fully observable for all subjects. However, in application, the cause of failure may not always be observable, thus impeding the risk assessment. In some extreme cases, none of the causes of failure is observable. In the case of a recurrent episode of Plasmodium vivax malaria following treatment, the patient may have suffered a relapse from a previous infection or acquired a new infection from a mosquito bite. In this case, the time to relapse cannot be modeled when a competing risk, a new infection, is present. The efficacy of a treatment for preventing relapse from a previous infection may be underestimated when the true cause of infection cannot be classified. In this paper, we developed a novel method for classifying the latent cause of failure under a competing risks set-up, which uses not only time to event information but also transition likelihoods between covariates at the baseline and at the time of event occurrence. Our classifier shows superior performance under various scenarios in simulation experiments. The method was applied to Plasmodium vivax infection data to classify recurrent infections of malaria.

Whole genome sequencing of Plasmodium vivax isolates reveals frequent sequence and structural polymorphisms in erythrocyte binding genes

Plasmodium vivax malaria is much less common in Africa than the rest of the world because the parasite relies primarily on the Duffy antigen/chemokine receptor (DARC) to invade human erythrocytes, and the majority of Africans are Duffy negative. Recently, there has been a dramatic increase in the reporting of P. vivax cases in Africa, with a high number of them being in Duffy negative individuals, potentially indicating P. vivax has evolved an alternative invasion mechanism that can overcome Duffy negativity. Here, we analyzed single nucleotide polymorphism (SNP) and copy number variation (CNV) in Whole Genome Sequence (WGS) data from 44 P. vivax samples isolated from symptomatic malaria patients in southwestern Ethiopia, where both Duffy positive and Duffy negative individuals are found. A total of 123,711 SNPs were detected, of which 22.7% were nonsynonymous and 77.3% were synonymous mutations. The largest number of SNPs were detected on chromosomes 9 (24,007 SNPs; 19.4% of total) and 10 (16,852 SNPs, 13.6% of total). There were particularly high levels of polymorphism in erythrocyte binding gene candidates including merozoite surface protein 1 (MSP1) and merozoite surface protein 3 (MSP3.5, MSP3.85 and MSP3.9). Two genes, MAEBL and MSP3.8 related to immunogenicity and erythrocyte binding function were detected with significant signals of positive selection. Variation in gene copy number was also concentrated in genes involved in host-parasite interactions, including the expansion of the Duffy binding protein gene (PvDBP) on chromosome 6 and MSP3.11 on chromosome 10. Based on the phylogeny constructed from the whole genome sequences, the expansion of these genes was an independent process among the P. vivax lineages in Ethiopia. We further inferred transmission patterns of P. vivax infections among study sites and showed various levels of gene flow at a small geographical scale. The genomic features of P. vivax provided baseline data for future comparison with those in Duffy-negative individuals and allowed us to develop a panel of informative Single Nucleotide Polymorphic markers diagnostic at a micro-geographical scale.

Plasmodium vivax is the most geographically widespread parasite species that causes malaria in humans. Although it occurs in Africa as a member of a mix of Plasmodium species, P. vivax is dominant in other parts of the world outside of Africa (e.g., Brazil). It was previously thought that most African populations were immune to P. vivax infections due to the absence of Duffy antigen chemokine receptor (DARC) gene expression required for erythrocyte invasion. However, several recent reports have indicated the emergence and potential spread of P. vivax across human populations in Africa. Compared to Southeast Asia and South America where P. vivax is highly endemic, data on polymorphisms in erythrocyte binding gene candidates of P. vivax from Africa is limited. Filling this knowlege gap is critical for identifying functional genes in erythrocyte invasion, biomarkers for tracking the P. vivax isolates from Africa, as well as potential gene targets for vaccine development. This paper examined the level of genetic polymorphisms in a panel of 43 potential erythrocyte binding protein genes based on whole genome sequences and described transmission patterns of P. vivax infections from different study sites in Ethiopia based on the genetic variants. Our analyses showed that chromosomes 9 and 10 of the P. vivax genomes isolated in Ethiopia had the most high-quality genetic polymorphisms. Among all erythrocyte binding protein gene candidates, the merozoite surface proteins 1 and merozoite surface protein 3 showed high levels of polymorphism. MAEBL and MSP3.8 related to immunogenicity and erythrocyte binding function were detected with significant signals of positive selection. The expansion of the Duffy binding protein and merozoite surface protein 3 gene copies was an independent process among the P. vivax lineages in Ethiopia. Various levels of gene flow were observed even at a smaller geographical scale. Our study provided baseline data for future comparison with P. vivax in Duffy negative individuals and help develop a panel of genetic markers that are informative at a micro-geographical scale.

Analysis of Plasmodium vivax schizont transcriptomes from field isolates reveals heterogeneity of expression of genes involved in host-parasite interactions

Plasmodium vivax gene regulation remains difficult to study due to the lack of a robust in vitro culture method, low parasite densities in peripheral circulation and asynchronous parasite development. We adapted an RNA-seq protocol “DAFT-seq” to sequence the transcriptome of four P. vivax field isolates that were cultured for a short period ex vivo before using a density gradient for schizont enrichment. Transcription was detected from 78% of the PvP01 reference genome, despite being schizont-enriched samples. This extensive data was used to define thousands of 5′ and 3′ untranslated regions, some of which overlapped with neighbouring transcripts, and to improve the gene models of 352 genes, including identifying 20 novel gene transcripts. This dataset has also significantly increased the known amount of heterogeneity between P. vivax schizont transcriptomes from individual patients. The majority of genes found to be differentially expressed between the isolates lack Plasmodium falciparum homologs and are predicted to be involved in host-parasite interactions, with an enrichment in reticulocyte binding proteins, merozoite surface proteins and exported proteins with unknown function. An improved understanding of the diversity within P. vivax transcriptomes will be essential for the prioritisation of novel vaccine targets.

Dynamics of Plasmodium vivax populations in border areas of the Greater Mekong sub-region during malaria elimination

Background

Countries within the Greater Mekong Sub-region (GMS) of Southeast Asia have committed to eliminating malaria by 2030. Although the malaria situation has greatly improved, malaria transmission remains at international border regions. In some areas, Plasmodium vivax has become the predominant parasite. To gain a better understanding of transmission dynamics, knowledge on the changes of P. vivax populations after the scale-up of control interventions will guide more effective targeted control efforts.

Methods

This study investigated genetic diversity and population structures in 206 P. vivax clinical samples collected at two time points in two international border areas: the China-Myanmar border (CMB) (n = 50 in 2004 and n = 52 in 2016) and Thailand-Myanmar border (TMB) (n = 50 in 2012 and n = 54 in 2015). Parasites were genotyped using 10 microsatellite markers.

Results

Despite intensified control efforts, genetic diversity remained high (H_E = 0.66–0.86) and was not significantly different among the four populations (P > 0.05). Specifically, H_E slightly decreased from 0.76 in 2004 to 0.66 in 2016 at the CMB and increased from 0.80 in 2012 to 0.86 in 2015 at the TMB. The proportions of polyclonal infections varied significantly among the four populations (P < 0.05), and showed substantial decreases from 48.0% in 2004 to 23.7 at the CMB and from 40.0% in 2012 to 30.7% in 2015 at the TMB, with corresponding decreases in the multiplicity of infection. Consistent with the continuous decline of malaria incidence in the GMS over time, there were also increases in multilocus linkage disequilibrium, suggesting more fragmented and increasingly inbred parasite populations. There were considerable genetic differentiation and sub-division among the four tested populations. Temporal genetic differentiation was observed at each site (F_ST = 0.081 at the CMB and F_ST = 0.133 at the TMB). Various degrees of clustering were evident between the older parasite samples collected in 2004 at the CMB and the 2016 CMB and 2012 TMB populations, suggesting some of these parasites had shared ancestry. In contrast, the 2015 TMB population was genetically distinctive, which may reflect a process of population replacement. Whereas the effective population size (N_e) at the CMB showed a decrease from 4979 in 2004 to 3052 in 2016 with the infinite allele model, the N_e at the TMB experienced an increase from 6289 to 10,259.

Conclusions

With enhanced control efforts on malaria, P. vivax at the TMB and CMB showed considerable spatial and temporal differentiation, but the presence of large P. vivax reservoirs still sustained genetic diversity and transmission. These findings provide new insights into P. vivax transmission dynamics and population structure in these border areas of the GMS. Coordinated and integrated control efforts on both sides of international borders are essential to reach the goal of regional malaria elimination.

Malaria Molecular Epidemiology: An Evolutionary Genetics Perspective

Malaria is a vector-borne disease that involves multiple parasite species in a variety of ecological settings. However, the parasite species causing the disease, the prevalence of subclinical infections, the emergence of drug resistance, the scale-up of interventions, and the ecological factors affecting malaria transmission, among others, are aspects that vary across areas where malaria is endemic. Such complexities have propelled the study of parasite genetic diversity patterns in the context of epidemiologic investigations. Importantly, molecular studies indicate that the time and spatial distribution of malaria cases reflect epidemiologic processes that cannot be fully understood without characterizing the evolutionary forces shaping parasite population genetic patterns. Although broad in scope, this review in the Microbiology Spectrum Curated Collection: Advances in Molecular Epidemiology highlights the need for understanding population genetic concepts when interpreting parasite molecular data. First, we discuss malaria complexity in terms of the parasite species involved. Second, we describe how molecular data are changing our understanding of malaria incidence and infectiousness. Third, we compare different approaches to generate parasite genetic information in the context of epidemiologically relevant questions related to malaria control. Finally, we describe a few Plasmodium genomic studies as evidence of how these approaches will provide new insights into the malaria disease dynamics. *This article is part of a curated collection.

Efficient Transmission of Mixed Plasmodium falciparum/vivax Infections From Humans to Mosquitoes

BACKGROUND

In Southeast Asia, people are often coinfected with different species of malaria (Plasmodium falciparum [Pf] and Plasmodium vivax [Pv]) as well as with multiple clones of the same species. Whether particular species or clones within mixed infections are more readily transmitted to mosquitoes remains unknown.

METHODS

Laboratory-reared Anopheles dirus were fed on blood from 119 Pf-infected Cambodian adults, with 5950 dissected to evaluate for transmitted infection. Among 12 persons who infected mosquitoes, polymerase chain reaction and amplicon deep sequencing were used to track species and clone-specific transmission to mosquitoes.

RESULTS

Seven of 12 persons that infected mosquitoes harbored mixed Pf/Pv infection. Among these 7 persons, all transmitted Pv with 2 transmitting both Pf and Pv, leading to Pf/Pv coinfection in 21% of infected mosquitoes. Up to 4 clones of each species were detected within persons. Shifts in clone frequency were detected during transmission. However, in general, all parasite clones in humans were transmitted to mosquitoes, with individual mosquitoes frequently carrying multiple transmitted clones.

CONCLUSIONS

Malaria diversity in human hosts was maintained in the parasite populations recovered from mosquitoes fed on their blood. However, in persons with mixed Pf/Pv malaria, Pv appears to be transmitted more readily, in association with more prevalent patent gametocytemia.

Recrudescence, Reinfection, or Relapse? A More Rigorous Framework to Assess Chloroquine Efficacy for Plasmodium vivax Malaria

Background

Plasmodium vivax resistance to chloroquine (CQ) has been reported worldwide, although the World Health Organization clinical drug efficacy studies protocol does not permit classification of patient outcomes.

Methods

We enrolled 40 patients with P. vivax malaria in northeastern Cambodia, where >17% treatment failures were previously reported. Patients were treated with CQ (30 mg/kg) and followed for 2 months, with frequent clinical examination and capillary blood sample collection for microscopy, molecular parasite detection and genotyping, and drug concentration measurements. Reinfections were prevented by relocating patients to a transmission-free area.

Results

P. vivax parasites were eliminated in all patients by day 3. Genomic analyses revealed that all clones in polyclonal infections were cleared at the same rate, indicating their equal susceptibility to CQ. CQ blood concentrations were below the therapeutic level in all recurrent infections (24 of 40 patients), which were efficiently cleared by a second course of CQ treatment. Genotyping (128 SNPs barcode) and sequences of entire parasite genome (Whole-Genome Sequencing, Illumina) indicated that two thirds (6 of 8) of the recurrent parasites resulted from heterologous relapses whose 50% are from by sibling/recombinant clones.

Conclusions

No evidence of CQ resistance was observed. Our data suggest that P. vivax antimalarial drug resistance is likely overestimated and that the current guidelines for clinical drug studies of P. vivax malaria need to be revised.

Humoral immunity prevents clinical malaria during Plasmodium relapses without eliminating gametocytes

Plasmodium relapses are attributed to the activation of dormant liver-stage parasites and are responsible for a significant number of recurring malaria blood-stage infections. While characteristic of human infections caused by P. vivax and P. ovale, their relative contribution to malaria disease burden and transmission remains poorly understood. This is largely because it is difficult to identify ‘bona fide’ relapse infections due to ongoing transmission in most endemic areas. Here, we use the P. cynomolgi–rhesus macaque model of relapsing malaria to demonstrate that clinical immunity can form after a single sporozoite-initiated blood-stage infection and prevent illness during relapses and homologous reinfections. By integrating data from whole blood RNA-sequencing, flow cytometry, P. cynomolgi-specific ELISAs, and opsonic phagocytosis assays, we demonstrate that this immunity is associated with a rapid recall response by memory B cells that expand and produce anti-parasite IgG1 that can mediate parasite clearance of relapsing parasites. The reduction in parasitemia during relapses was mirrored by a reduction in the total number of circulating gametocytes, but importantly, the cumulative proportion of gametocytes increased during relapses. Overall, this study reveals that P. cynomolgi relapse infections can be clinically silent in macaques due to rapid memory B cell responses that help to clear asexual-stage parasites but still carry gametocytes.

Plasmodium vivax contributes significantly to global malaria morbidity and remains a major obstacle for malaria elimination due to its ability to form dormant stages in the liver. These forms can become activated to cause relapsing blood-stage infections. Relapses remain poorly understood because it is difficult to verify whether P. vivax blood-stage infections in patients are due to new infections or relapses in most cases. Here, we use a nonhuman primate model of Plasmodium vivax malaria in concert with state-of-the-art immunological and molecular techniques to assess pathogenesis, host responses, and circulating gametocyte levels during relapses. We found that relapses were clinically silent compared to initial infections, and they were associated with a robust memory B cell response. This response resulted in the production of antibodies that were able to mediate clearance of asexual parasites. Despite this rapid immune protection, the sexual-stage gametocytes continued to circulate. Our study provides mechanistic insights into the host-parasite interface during Plasmodium relapse infections and demonstrates that clinically silent relapses can harbor gametocytes that may be infectious to mosquitoes.

Human Migration and the Spread of the Nematode Parasite Wuchereria bancrofti

The human disease lymphatic filariasis causes the debilitating effects of elephantiasis and hydrocele. Lymphatic filariasis currently affects the lives of 90 million people in 52 countries. There are three nematodes that cause lymphatic filariasis, Brugia malayi, Brugia timori, and Wuchereria bancrofti, but 90% of all cases of lymphatic filariasis are caused solely by W. bancrofti (Wb). Here we use population genomics to reconstruct the probable route and timing of migration of Wb strains that currently infect Africa, Haiti, and Papua New Guinea (PNG). We used selective whole genome amplification to sequence 42 whole genomes of single Wb worms from populations in Haiti, Mali, Kenya, and PNG. Our results are consistent with a hypothesis of an Island Southeast Asia or East Asian origin of Wb. Our demographic models support divergence times that correlate with the migration of human populations. We hypothesize that PNG was infected at two separate times, first by the Melanesians and later by the migrating Austronesians. The migrating Austronesians also likely introduced Wb to Madagascar where later migrations spread it to continental Africa. From Africa, Wb spread to the New World during the transatlantic slave trade. Genome scans identified 17 genes that were highly differentiated among Wb populations. Among these are genes associated with human immune suppression, insecticide sensitivity, and proposed drug targets. Identifying the distribution of genetic diversity in Wb populations and selection forces acting on the genome will build a foundation to test future hypotheses and help predict response to current eradication efforts.

Plasmodium vivax transcriptomes reveal stage-specific chloroquine response and differential regulation of male and female gametocytes

Studies of Plasmodium vivax gene expression are complicated by the lack of in vitro culture system and the difficulties associated with studying clinical infections that often contain multiple clones and a mixture of parasite stages. Here, we characterize the transcriptomes of P. vivax parasites from 26 malaria patients. We show that most parasite mRNAs derive from trophozoites and that the asynchronicity of P. vivax infections is therefore unlikely to confound gene expression studies. Analyses of gametocyte genes reveal two distinct clusters of co-regulated genes, suggesting that male and female gametocytes are independently regulated. Finally, we analyze gene expression changes induced by chloroquine and show that this antimalarial drug efficiently eliminates most P. vivax parasite stages but, in contrast to P. falciparum, does not affect trophozoites.

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.

Genetic diversity in two Plasmodium vivax protein ligands for reticulocyte invasion

The interaction between Plasmodium vivax Duffy binding protein (PvDBP) and Duffy antigen receptor for chemokines (DARC) has been described as critical for the invasion of human reticulocytes, although increasing reports of P. vivax infections in Duffy-negative individuals questions its unique role. To investigate the genetic diversity of the two main protein ligands for reticulocyte invasion, PvDBP and P. vivax Erythrocyte Binding Protein (PvEBP), we analyzed 458 isolates collected in Cambodia and Madagascar from individuals genotyped as Duffy-positive. First, we observed a high proportion of isolates with multiple copies PvEBP from Madagascar (56%) where Duffy negative and positive individuals coexist compared to Cambodia (19%) where Duffy-negative population is virtually absent. Whether the gene amplification observed is responsible for alternate invasion pathways remains to be tested. Second, we found that the PvEBP gene was less diverse than PvDBP gene (12 vs. 33 alleles) but provided evidence for an excess of nonsynonymous mutations with the complete absence of synonymous mutations. This finding reveals that PvEBP is under strong diversifying selection, and confirms the importance of this protein ligand in the invasion process of the human reticulocytes and as a target of acquired immunity. These observations highlight how genomic changes in parasite ligands improve the fitness of P. vivax isolates in the face of immune pressure and receptor polymorphisms.

Large and finite sample properties of a maximum-likelihood estimator for multiplicity of infection

Reliable measures of transmission intensities can be incorporated into metrics for monitoring disease-control interventions. Genetic (molecular) measures like multiplicity of infection (MOI) have several advantages compared with traditional measures, e.g., R0. Here, we investigate the properties of a maximum-likelihood approach to estimate MOI and pathogen-lineage frequencies. By verifying regulatory conditions, we prove asymptotical unbiasedness, consistency and efficiency of the estimator. Finite sample properties concerning bias and variance are evaluated over a comprehensive parameter range by a systematic simulation study. Moreover, the estimator's sensitivity to model violations is studied. The estimator performs well for realistic sample sizes and parameter ranges. In particular, the lineage-frequency estimates are almost unbiased independently of sample size. The MOI estimate's bias vanishes with increasing sample size, but might be substantial if sample size is too small. The estimator's variance matrix agrees well with the Cramér-Rao lower bound, even for small sample size. The numerical and analytical results of this study can be used for study design. This is exemplified by a malaria data set from Venezuela. It is shown how the results can be used to determine the necessary sample size to achieve certain performance goals. An implementation of the likelihood method and a simulation algorithm for study design, implemented as an R script, is available as S1 File alongside a documentation (S2 File) and example data (S3 File).

Plasmodium vivax Biology: Insights Provided by Genomics, Transcriptomics and Proteomics

During the last decade, the vast omics field has revolutionized biological research, especially the genomics, transcriptomics and proteomics branches, as technological tools become available to the field researcher and allow difficult question-driven studies to be addressed. Parasitology has greatly benefited from next generation sequencing (NGS) projects, which have resulted in a broadened comprehension of basic parasite molecular biology, ecology and epidemiology. Malariology is one example where application of this technology has greatly contributed to a better understanding of Plasmodium spp. biology and host-parasite interactions. Among the several parasite species that cause human malaria, the neglected Plasmodium vivax presents great research challenges, as in vitro culturing is not yet feasible and functional assays are heavily limited. Therefore, there are gaps in our P. vivax biology knowledge that affect decisions for control policies aiming to eradicate vivax malaria in the near future. In this review, we provide a snapshot of key discoveries already achieved in P. vivax sequencing projects, focusing on developments, hurdles, and limitations currently faced by the research community, as well as perspectives on future vivax malaria research.

Genomic Analyses Reveal the Common Occurrence and Complexity of Plasmodium vivax Relapses in Cambodia

Plasmodium vivax parasites have a unique dormant stage that can cause relapses weeks or months after the initial infection. These dormant parasites are among the main challenges of vivax malaria control as they constitute a reservoir that is difficult to eliminate. Since field studies are confounded by reinfections and possible recrudescence of drug-resistant parasites, most analyses of P. vivax relapses have focused on travelers returning from regions of malaria endemicity. However, it is not clear whether these individuals accurately recapitulate the relapse patterns of repeatedly infected individuals residing in areas of endemicity. Here, we present analyses of vivax malaria patients enrolled in a tightly controlled field study in Cambodia. After antimalarial drug treatment was administered, we relocated 20 individuals to a nontransmission area and followed them for 60 days, with blood collection performed every second day. Our analyses reveal that 60% of the patients relapsed during the monitoring period. Using whole-genome sequencing and high-throughput genotyping, we showed that relapses in Cambodia are often polyclonal and that the relapsing parasites harbor various degrees of relatedness to the parasites present in the initial infection. Our analyses also showed that clone populations differed dynamically, with new clones emerging during the course of the relapsing infections. Overall, our study data show that it is possible to investigate the patterns, dynamics, and diversity of P. vivax relapses of individuals living in a region of malaria endemicity and reveal that P. vivax relapses are much more pervasive and complex than previously considered. (This study has been registered at ClinicalTrials.gov under registration no. NCT02118090)IMPORTANCEP. vivax parasites can remain dormant in the liver and relapse weeks or months after the initial infection, greatly complicating malaria control and elimination efforts. The few investigations of this dormant stage have relied on travelers and military personnel returning from areas of malaria endemicity. However, it is not clear whether these individuals, exposed to a limited number of infections, accurately represent the patterns of relapses of individuals living in areas of endemicity, who are repeatedly infected by P. vivax parasites. Our study combined tightly controlled fieldwork with comprehensive genomic analyses, and our report provides a first opportunity to investigate the patterns, dynamics, and diversity of P. vivax relapses directly with individuals living in areas of endemicity.

Frequent Spread of Plasmodium vivax Malaria Maintains High Genetic Diversity at the Myanmar-China Border, Without Distance and Landscape Barriers

Background

In Myanmar, civil unrest and the establishment of internally displaced person (IDP) settlements along the Myanmar-China border have impacted malaria transmission.

Methods

Microsatellite markers were used to examine source-sink dynamics for Plasmodium vivax between IDP settlements and surrounding villages in the border region. Genotypic structure and diversity were compared across the 3 years following the establishment of IDP settlements, to infer demographic history. We investigated whether human migration and landscape heterogeneity contributed to P. vivax transmission.

Results

P. vivax from IDP settlements and local communities consistently exhibited high genetic diversity within populations but low polyclonality within individuals. No apparent genetic structure was observed among populations and years. P. vivax genotypes in China were similar to those in Myanmar, and parasite introduction was unidirectional. Landscape factors, including distance, elevation, and land cover, do not appear to impede parasite gene flow.

Conclusions

The admixture of P. vivax genotypes suggested that parasite gene flow via human movement contributes to the spread of malaria both locally in Myanmar and across the international border. Our genetic findings highlight the presence of large P. vivax gene reservoirs that can sustain transmission. Thus, it is important to reinforce and improve existing control efforts along border areas.

Nationwide genetic surveillance of Plasmodium vivax in Papua New Guinea reveals heterogeneous transmission dynamics and routes of migration amongst subdivided populations

The Asia Pacific Leaders in Malaria Alliance (APLMA) have committed to eliminate malaria from the region by 2030. Papua New Guinea (PNG) has the highest malaria burden in the Asia-Pacific region but with the intensification of control efforts since 2005, transmission has been dramatically reduced and Plasmodium vivax is now the dominant malaria infection in some parts of the country. To gain a better understanding of the transmission dynamics and migration patterns of P. vivax in PNG, here we investigate population structure in eight geographically and ecologically distinct regions of the country. A total of 219 P. vivax isolates (16-30 per population) were successfully haplotyped using 10 microsatellite markers. A wide range of genetic diversity (H_e=0.37-0.87, R_s=3.60-7.58) and significant multilocus linkage disequilibrium (LD) was observed in six of the eight populations (I_A^S=0.08-0.15 p-value<0.05) reflecting a spectrum of transmission intensities across the country. Genetic differentiation between regions was evident (Jost's D=0.07-0.72), with increasing divergence of populations with geographic distance. Overall, P. vivax isolates clustered into three major genetic populations subdividing the Mainland lowland and coastal regions, the Islands and the Highlands. P. vivax gene flow follows major human migration routes, and there was higher gene flow amongst Mainland parasite populations than among Island populations. The Central Province (samples collected in villages close to the capital city, Port Moresby), acts as a sink for imported infections from the three major endemic areas. These insights into P. vivax transmission dynamics and population networks will inform targeted strategies to contain malaria infections and to prevent the spread of drug resistance in PNG.

Determination of multiple-clone infection at allelic dimorphism site of Plasmodium vivax merozoite surface protein-1 in the Republic of Korea by pyrosequencing assay

Allelic diversity leading to multiple gene polymorphisms of vivax malaria parasites has been shown to greatly contribute to antigenic variation and drug resistance, increasing the potential for multiple-clone infections within the host. Therefore, to identify multiple-clone infections and the predominant haplotype of Plasmodium vivax in a South Korean population, P. vivax merozoite surface protein-1 (PvMSP-1) was analyzed by pyrosequencing. Pyrosequencing of 156 vivax malaria-infected samples yielded 97 (62.18%) output pyrograms showing two main types of peak patterns of the dimorphic allele for threonine and alanine (T1476A). Most of the samples evaluated (88.66%) carried multiple-clone infections (wild- and mutant-types), whereas 11.34% of the same population carried only the mutant-type (1476A). In addition, each allele showed a high frequency of guanine (G) base substitution at both the first and third positions (86.07% and 81.13%, respectively) of the nucleotide combinations. Pyrosequencing of the PvMSP-1 42-kDa fragment revealed a heterogeneous parasite population, with the mutant-type dominant compared to the wild-type. Understanding the genetic diversity and multiple-clone infection rates may lead to improvements in vivax malaria prevention and strategic control plans. Further studies are needed to improve the efficacy of the pyrosequencing assay with large sample sizes and additional nucleotide positions.

Development of amplicon deep sequencing markers and data analysis pipeline for genotyping multi-clonal malaria infections

Background

Amplicon deep sequencing permits sensitive detection of minority clones and improves discriminatory power for genotyping multi-clone Plasmodium falciparum infections. New amplicon sequencing and data analysis protocols are needed for genotyping in epidemiological studies and drug efficacy trials of P. falciparum.

Methods

Targeted sequencing of molecular marker csp and novel marker cpmp was conducted in duplicate on mixtures of parasite culture strains and 37 field samples. A protocol allowing to multiplex up to 384 samples in a single sequencing run was applied. Software “HaplotypR” was developed for data analysis.

Results

Cpmp was highly diverse (H_e = 0.96) in contrast to csp (H_e = 0.57). Minority clones were robustly detected if their frequency was >1%. False haplotype calls owing to sequencing errors were observed below that threshold.

Conclusions

To reliably detect haplotypes at very low frequencies, experiments are best performed in duplicate and should aim for coverage of >10′000 reads/amplicon. When compared to length polymorphic marker msp2, highly multiplexed amplicon sequencing displayed greater sensitivity in detecting minority clones.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4260-y) contains supplementary material, which is available to authorized users.

Malaria Epidemiology at the Clone Level

Genotyping to distinguish between parasite clones is nowadays a standard in many molecular epidemiological studies of malaria. It has become crucial in drug trials and to follow individual clones in epidemiological studies, and to understand how drug resistance emerges and spreads. Here, we review the applications of the increasingly available genotyping tools and whole-genome sequencing data, and argue for a better integration of population genetics findings into malaria-control strategies.

Characterization of P. vivax blood stage transcriptomes from field isolates reveals similarities among infections and complex gene isoforms

Our understanding of the structure and regulation of Plasmodium vivax genes is limited by our inability to grow the parasites in long-term in vitro cultures. Most P. vivax studies must therefore rely on patient samples, which typically display a low proportion of parasites and asynchronous parasites. Here, we present stranded RNA-seq data generated directly from a small volume of blood from three Cambodian vivax malaria patients collected before treatment. Our analyses show surprising similarities of the parasite gene expression patterns across infections, despite extensive variations in parasite stage proportion. These similarities contrast with the unique gene expression patterns observed in sporozoites isolated from salivary glands of infected Colombian mosquitoes. Our analyses also indicate that more than 10% of P. vivax genes encode multiple, often undescribed, protein-coding sequences, potentially increasing the diversity of proteins synthesized by blood stage parasites. These data also greatly improve the annotations of P. vivax gene untranslated regions, providing an important resource for future studies of specific genes.

Genome-wide scans for the identification of Plasmodium vivax genes under positive selection

BACKGROUND

The current trend of Plasmodium vivax cases imported from Southeast Asia into China has sharply increased recently, especially from the China-Myanmar border (CMB) area. High recombination rates of P. vivax populations associated with varied transmission intensity might cause distinct local selective pressures. The information on the genetic variability of P. vivax in this area is scant. Hence, this study assessed the genetic diversity of P. vivax genome sequence in CMB area and aimed to provide information on the positive selection of new gene loci.

RESULTS

This study reports a genome-wide survey of P. vivax in CMB area, using blood samples from local patients to identify population-specific selective processes. The result showed that considerable genetic diversity and mean pair-wise divergence among the sequenced P. vivax isolates were higher in some important gene families. Using the standardized integrated haplotype score (|iHS|) for all SNPs in chromosomal regions with SNPs above the top 1% distribution, it was observed that the top score locus involved 356 genes and most of them are associated with red blood cell invasion and immune evasion. The XP-EHH test was also applied and some important genes associated with anti-malarial drug resistance were observed in high positive scores list. This result suggests that P. vivax in CMB area is facing more pressure to survive than any other region and this has led to the strong positive selection of genes that are associated with host-parasite interactions.

CONCLUSIONS

This study suggests that greater genetic diversity in P. vivax from CMB area and positive selection signals in invasion and drug resistance genes are consistent with the history of drug use during malaria elimination programme in CMB area. Furthermore, this result also demonstrates that haplotype-based detecting selection can assist the genome-wide methods to identify the determinants of P. vivax diversity.

THE REAL McCOIL: A method for the concurrent estimation of the complexity of infection and SNP allele frequency for malaria parasites

As many malaria-endemic countries move towards elimination of Plasmodium falciparum, the most virulent human malaria parasite, effective tools for monitoring malaria epidemiology are urgent priorities. P. falciparum population genetic approaches offer promising tools for understanding transmission and spread of the disease, but a high prevalence of multi-clone or polygenomic infections can render estimation of even the most basic parameters, such as allele frequencies, challenging. A previous method, COIL, was developed to estimate complexity of infection (COI) from single nucleotide polymorphism (SNP) data, but relies on monogenomic infections to estimate allele frequencies or requires external allele frequency data which may not available. Estimates limited to monogenomic infections may not be representative, however, and when the average COI is high, they can be difficult or impossible to obtain. Therefore, we developed THE REAL McCOIL, Turning HEterozygous SNP data into Robust Estimates of ALelle frequency, via Markov chain Monte Carlo, and Complexity Of Infection using Likelihood, to incorporate polygenomic samples and simultaneously estimate allele frequency and COI. This approach was tested via simulations then applied to SNP data from cross-sectional surveys performed in three Ugandan sites with varying malaria transmission. We show that THE REAL McCOIL consistently outperforms COIL on simulated data, particularly when most infections are polygenomic. Using field data we show that, unlike with COIL, we can distinguish epidemiologically relevant differences in COI between and within these sites. Surprisingly, for example, we estimated high average COI in a peri-urban subregion with lower transmission intensity, suggesting that many of these cases were imported from surrounding regions with higher transmission intensity. THE REAL McCOIL therefore provides a robust tool for understanding the molecular epidemiology of malaria across transmission settings.

Selective sweep suggests transcriptional regulation may underlie Plasmodium vivax resilience to malaria control measures in Cambodia

Cambodia, in which both Plasmodium vivax and Plasmodium falciparum are endemic, has been the focus of numerous malaria-control interventions, resulting in a marked decline in overall malaria incidence. Despite this decline, the number of P vivax cases has actually increased. To understand better the factors underlying this resilience, we compared the genetic responses of the two species to recent selective pressures. We sequenced and studied the genomes of 70 P vivax and 80 P falciparum isolates collected between 2009 and 2013. We found that although P falciparum has undergone population fracturing, the coendemic P vivax population has grown undisrupted, resulting in a larger effective population size, no discernable population structure, and frequent multiclonal infections. Signatures of selection suggest recent, species-specific evolutionary differences. Particularly, in contrast to P falciparum, P vivax transcription factors, chromatin modifiers, and histone deacetylases have undergone strong directional selection, including a particularly strong selective sweep at an AP2 transcription factor. Together, our findings point to different population-level adaptive mechanisms used by P vivax and P falciparum parasites. Although population substructuring in P falciparum has resulted in clonal outgrowths of resistant parasites, P vivax may use a nuanced transcriptional regulatory approach to population maintenance, enabling it to preserve a larger, more diverse population better suited to facing selective threats. We conclude that transcriptional control may underlie P vivax's resilience to malaria control measures. Novel strategies to target such processes are likely required to eradicate P vivax and achieve malaria elimination.

Dissecting malaria biology and epidemiology using population genetics and genomics

Molecular approaches have an increasingly recognized utility in surveillance of malaria parasite populations, not only in defining prevalence and incidence with higher sensitivity than traditional methods, but also in monitoring local and regional parasite transmission patterns. In this review, we provide an overview of population genetic and genomic studies of human-infecting Plasmodium species, highlighting recent advances in the field. In accordance with the renewed impetus for malaria eradication, many studies are now using genetic and genomic epidemiology to support local evidence-based intervention strategies. Microsatellite genotyping remains a popular approach for both Plasmodium falciparum and Plasmodium vivax. However, with the increasing availability of whole genome sequencing data enabling effective single nucleotide polymorphism-based panels tailored to a given study question and setting, this approach is gaining popularity. The availability of new reference genomes for Plasmodium malariae and Plasmodium ovale should see a surge in similar molecular studies on these currently neglected species. Genomic studies are revealing new insights into important adaptive mechanisms of the parasite including antimalarial drug resistance. The advent of new methodologies such as selective whole genome amplification for dealing with extensive human DNA in low density field isolates should see genome-wide approaches becoming routine for parasite surveillance once the economic costs outweigh the current cost benefits of targeted approaches.