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

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

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

Origin of the human malaria parasite Plasmodium vivax

The geographic origin of Plasmodium vivax, a leading cause of human malaria, has been the subject of much speculation. Here we review the evolutionary history of P. vivax and P. vivax-like parasites in humans and non-human primates on three continents, providing overwhelming evidence for an African origin. This conclusion is consistent with recent reports showing that Duffy-negative humans in Africa are, in fact, susceptible to P. vivax, with parasites invading Duffy-antigen-expressing erythroid precursors. Thus, the African origin of P. vivax not only explains the distribution of the Duffy-negative genotype but also provides new insight into the history and status of P. vivax malaria in Africa and efforts geared toward its eradication.

Ancient Plasmodium genomes shed light on the history of human malaria

Malaria-causing protozoa of the genus Plasmodium have exerted one of the strongest selective pressures on the human genome, and resistance alleles provide biomolecular footprints that outline the historical reach of these species1. Nevertheless, debate persists over when and how malaria parasites emerged as human pathogens and spread around the globe1,2. To address these questions, we generated high-coverage ancient mitochondrial and nuclear genome-wide data from P. falciparum, P. vivax and P. malariae from 16 countries spanning around 5,500 years of human history. We identified P. vivax and P. falciparum across geographically disparate regions of Eurasia from as early as the fourth and first millennia BCE, respectively; for P. vivax, this evidence pre-dates textual references by several millennia3. Genomic analysis supports distinct disease histories for P. falciparum and P. vivax in the Americas: similarities between now-eliminated European and peri-contact South American strains indicate that European colonizers were the source of American P. vivax, whereas the trans-Atlantic slave trade probably introduced P. falciparum into the Americas. Our data underscore the role of cross-cultural contacts in the dissemination of malaria, laying the biomolecular foundation for future palaeo-epidemiological research into the impact of Plasmodium parasites on human history. Finally, our unexpected discovery of P. falciparum in the high-altitude Himalayas provides a rare case study in which individual mobility can be inferred from infection status, adding to our knowledge of cross-cultural connectivity in the region nearly three millennia ago.

High-Throughput Genotyping of Plasmodium vivax in the Peruvian Amazon via Molecular Inversion Probes

Plasmodium vivax transmission occurs throughout the tropics and is an emerging threat in areas of Plasmodium falciparum decline, causing relapse infections that complicate treatment and control. Targeted sequencing for P. falciparum has been widely deployed to detect population structure and the geographic spread of antimalarial and diagnostic resistance. However, there are fewer such tools for P. vivax . Leveraging global variation data, we designed four molecular inversion probe (MIP) genotyping panels targeting geographically differentiating SNPs, neutral SNPs, putative antimalarial resistance genes, and vaccine candidate genes. We deployed these MIP panels on 866 infections from the Peruvian Amazon and identified transmission networks with clonality (IBD>0.99), copy number variation in Pvdbp and multiple Pvrbps , fixation of putative antimalarial resistance, and balancing selection in 13 vaccine candidate genes. Our MIP panels are the broadest genotyping panel currently available and are poised for successful deployment in other regions of P. vivax transmission.

Plasmodium vivax tryptophan-rich antigen reduces type I collagen secretion via the NF-κBp65 pathway in splenic fibroblasts

BACKGROUND

The spleen plays a critical role in the immune response against malaria parasite infection, where splenic fibroblasts (SFs) are abundantly present and contribute to immune function by secreting type I collagen (collagen I). The protein family is characterized by Plasmodium vivax tryptophan-rich antigens (PvTRAgs), comprising 40 members. PvTRAg23 has been reported to bind to human SFs (HSFs) and affect collagen I levels. Given the role of type I collagen in splenic immune function, it is important to investigate the functions of the other members within the PvTRAg protein family.

METHODS

Protein structural prediction was conducted utilizing bioinformatics analysis tools and software. A total of 23 PvTRAgs were successfully expressed and purified using an Escherichia coli prokaryotic expression system, and the purified proteins were used for co-culture with HSFs. The collagen I levels and collagen-related signaling pathway protein levels were detected by immunoblotting, and the relative expression levels of inflammatory factors were determined by quantitative real-time PCR.

RESULTS

In silico analysis showed that P. vivax has 40 genes encoding the TRAg family. The C-terminal region of all PvTRAgs is characterized by the presence of a domain rich in tryptophan residues. A total of 23 recombinant PvTRAgs were successfully expressed and purified. Only five PvTRAgs (PvTRAg5, PvTRAg16, PvTRAg23, PvTRAg30, and PvTRAg32) mediated the activation of the NF-κBp65 signaling pathway, which resulted in the production of inflammatory molecules and ultimately a significant reduction in collagen I levels in HSFs.

CONCLUSIONS

Our research contributes to the expansion of knowledge regarding the functional role of PvTRAgs, while it also enhances our understanding of the immune evasion mechanisms utilized by parasites.

Cellular dormancy: A widespread phenomenon that perpetuates infectious diseases

Under adverse environmental conditions, microorganisms are able to enter a state of cellular dormancy which consists of cell cycle arrest and interruption of multiplication. This process ensures their perpetuation in the infected host organism and enables the spread of disease. Throughout biological evolution, dormancy allowed microorganisms to persist in a harsh niche until favorable conditions for their reactivation were re-established. Here, we propose to discuss the dormancy of bacteria and protozoa pathogens focusing on the potential mechanisms and components associated with dormancy.

Population genomic evidence of structured and connected Plasmodium vivax populations under host selection in Latin America

Pathogen genomic epidemiology has the potential to provide a deep understanding of population dynamics, facilitating strategic planning of interventions, monitoring their impact, and enabling timely responses, and thereby supporting control and elimination efforts of parasitic tropical diseases. Plasmodium vivax, responsible for most malaria cases outside Africa, shows high genetic diversity at the population level, driven by factors like sub-patent infections, a hidden reservoir of hypnozoites, and early transmission to mosquitoes. While Latin America has made significant progress in controlling Plasmodium falciparum, it faces challenges with residual P. vivax. To characterize genetic diversity and population structure and dynamics, we have analyzed the largest collection of P. vivax genomes to date, including 1474 high-quality genomes from 31 countries across Asia, Africa, Oceania, and America. While P. vivax shows high genetic diversity globally, Latin American isolates form a distinctive population, which is further divided into sub-populations and occasional clonal pockets. Genetic diversity within the continent was associated with the intensity of transmission. Population differentiation exists between Central America and the North Coast of South America, vs. the Amazon Basin, with significant gene flow within the Amazon Basin, but limited connectivity between the Northwest Coast and the Amazon Basin. Shared genomic regions in these parasite populations indicate adaptive evolution, particularly in genes related to DNA replication, RNA processing, invasion, and motility - crucial for the parasite's survival in diverse environments. Understanding these population-level adaptations is crucial for effective control efforts, offering insights into potential mechanisms behind drug resistance, immune evasion, and transmission dynamics.

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.

The Genome of Plasmodium gonderi: Insights into the Evolution of Human Malaria Parasites

Plasmodium species causing malaria in humans are not monophyletic, sharing common ancestors with nonhuman primate parasites. Plasmodium gonderi is one of the few known Plasmodium species infecting African old-world monkeys that are not found in apes. This study reports a de novo assembled P. gonderi genome with complete chromosomes. The P. gonderi genome shares codon usage, syntenic blocks, and other characteristics with the human parasites Plasmodium ovale s.l. and Plasmodium malariae, also of African origin, and the human parasite Plasmodium vivax and species found in nonhuman primates from Southeast Asia. Using phylogenetically aware methods, newly identified syntenic blocks were found enriched with conserved metabolic genes. Regions outside those blocks harbored genes encoding proteins involved in the vertebrate host-Plasmodium relationship undergoing faster evolution. Such genome architecture may have facilitated colonizing vertebrate hosts. Phylogenomic analyses estimated the common ancestor between P. vivax and an African ape parasite P. vivax-like, within the Asian nonhuman primates parasites clade. Time estimates incorporating P. gonderi placed the P. vivax and P. vivax-like common ancestor in the late Pleistocene, a time of active migration of hominids between Africa and Asia. Thus, phylogenomic and time-tree analyses are consistent with an Asian origin for P. vivax and an introduction of P. vivax-like into Africa. Unlike other studies, time estimates for the clade with Plasmodium falciparum, the most lethal human malaria parasite, coincide with their host species radiation, African hominids. Overall, the newly assembled genome presented here has the quality to support comparative genomic investigations in Plasmodium.

Vivax malaria: a possible stumbling block for malaria elimination in India

Plasmodium vivax is geographically the most widely dispersed human malaria parasite species. It has shown resilience and a great deal of adaptability. Genomic studies suggest that P. vivax originated from Asia or Africa and moved to the rest of the world. Although P. vivax is evolutionarily an older species than Plasmodium falciparum, its biology, transmission, pathology, and control still require better elucidation. P. vivax poses problems for malaria elimination because of the ability of a single primary infection to produce multiple relapses over months and years. P. vivax malaria elimination program needs early diagnosis, and prompt and complete radical treatment, which is challenging, to simultaneously exterminate the circulating parasites and dormant hypnozoites lodged in the hepatocytes of the host liver. As prompt surveillance and effective treatments are rolled out, preventing primaquine toxicity in the patients having glucose-6-phosphate dehydrogenase (G6PD) deficiency should be a priority for the vivax elimination program. This review sheds light on the burden of P. vivax, changing epidemiological patterns, the hurdles in elimination efforts, and the essential tools needed not just in India but globally. These tools encompass innovative treatments for eliminating dormant parasites, coping with evolving drug resistance, and the development of potential vaccines against the parasite.

Host kinase regulation of Plasmodium vivax dormant and replicating liver stages

Upon transmission to the liver, Plasmodium vivax parasites form replicating schizonts, which continue to initiate blood-stage infection, or dormant hypnozoites that reactivate weeks to months after initial infection. P. vivax phenotypes in the field vary significantly, including the ratio of schizonts to hypnozoites formed and the frequency and timing of relapse. Evidence suggests that both parasite genetics and environmental factors underly this heterogeneity. We previously demonstrated that data on the effect of a panel of kinase inhibitors with overlapping targets on Plasmodium liver stage infection, in combination with a computational approach called kinase regression (KiR), can be used to uncover novel host regulators of infection. Here, we applied KiR to evaluate the extent to which P. vivax liver-stage parasites are susceptible to changes in host kinase activity. We identified a role for a subset of host kinases in regulating schizont and hypnozoite infection and schizont size and characterized overlap as well as variability in host phosphosignaling dependencies between parasite forms and across multiple patient isolates. Striking, our data point to variability in host dependencies across P. vivax isolates, suggesting one possible origin of the heterogeneity observed across P. vivax in the field.

A new Plasmodium vivax reference genome for South American isolates

BACKGROUND

Plasmodium vivax is the second most important cause of human malaria worldwide, and accounts for the majority of malaria cases in South America. A high-quality reference genome exists for Papua Indonesia (PvP01) and Thailand (PvW1), but is lacking for South America. A reference genome specifically for South America would be beneficial though, as P. vivax is a genetically diverse parasite with geographical clustering.

RESULTS

This study presents a new high-quality assembly of a South American P. vivax isolate, referred to as PvPAM (P. vivax Peruvian AMazon). The genome was obtained from a low input patient sample from the Peruvian Amazon and sequenced using PacBio technology, resulting in a highly complete assembly with 6497 functional genes. Telomeric ends were present in 17 out of 28 chromosomal ends, and additional (sub)telomeric regions are present in 12 unassigned contigs. A comparison of multigene families between PvPAM and the PvP01 genome revealed remarkable variation in vir genes, and the presence of merozoite surface proteins (MSP) 3.6 and 3.7. Three dhfr and dhps drug resistance associated mutations are present in PvPAM, similar to those found in other Peruvian isolates. Mapping of publicly available South American whole genome sequencing (WGS) data to PvPAM resulted in significantly fewer variants and truncated reads compared to the use of PvP01 or PvW1 as reference genomes. To minimize the number of core genome variants in non-South American samples, PvW1 is most suited for Southeast Asian isolates, both PvPAM and PvW1 are suited for South Asian isolates, and PvPAM is recommended for African isolates. Interestingly, non-South American samples still contained the least subtelomeric variants when mapped to PvPAM, indicating high quality of the PvPAM subtelomeric regions.

CONCLUSIONS

Our findings show that the PvPAM reference genome more accurately represents South American P. vivax isolates in comparison to PvP01 and PvW1. In addition, PvPAM has a high level of completeness, and contains a similar number of annotated genes as PvP01 or PvW1. The PvPAM genome therefore will be a valuable resource to improve future genomic analyses on P. vivax isolates from the South American continent.

Malaria Genomics, Vaccine Development, and Microbiome

Recent advances in malaria genetics and genomics have transformed many aspects of malaria research in areas of molecular evolution, epidemiology, transmission, host-parasite interaction, drug resistance, pathogenicity, and vaccine development. Here, in addition to introducing some background information on malaria parasite biology, parasite genetics/genomics, and genotyping methods, we discuss some applications of genetic and genomic approaches in vaccine development and in studying interactions with microbiota. Genetic and genomic data can be used to search for novel vaccine targets, design an effective vaccine strategy, identify protective antigens in a whole-organism vaccine, and evaluate the efficacy of a vaccine. Microbiota has been shown to influence disease outcomes and vaccine efficacy; studying the effects of microbiota in pathogenicity and immunity may provide information for disease control. Malaria genetics and genomics will continue to contribute greatly to many fields of malaria research.

How can the complex epidemiology of malaria in India impact its elimination?

Malaria is a human health hazard in the tropical and subtropical zones of the globe and is poised to be eliminated by the year 2030. Despite a decrease in incidence in the past two decades, many endemic countries, including India, report cases regularly. The epidemiology of malaria in India is unique owing to several features of the Plasmodium parasites, Anopheles vectors, ecoepidemiological situations conducive to disease transmission, and susceptible humans living in rural and forested areas. Limitations in public health reach, and poor health-seeking behaviour of vulnerable populations living in hard-to-reach areas, add to the problem. We bring all of these factors together in a comprehensive framework and opine that, in spite of complexities, targeted elimination of malaria in India is achievable with planned programmatic approaches.

A case of atypical vivax malaria with a global review of reports on myriads of morpho-variations in parasitized red blood cells

Plasmodium vivax, one of the major species associated with human malaria, continues to be a major public health problem in many parts of the world. Numerous studies related to vivax malaria have described quantitative haematological findings (level of haemoglobin, thrombocytopaenia, haematocrit values), but diverse morphological changes of parasite forms within infected red blood cells (iRBCs) have been mentioned only in few studies. Here we report a case of a 13-year-old boy who presented with fever, significant low platelet counts and hypovolaemia that created a diagnostic dilemma. Detection of microgametocytes by microscopic examinations, further confirmed by multiplex nested PCR assays and response to anti-malarials, helped to make the diagnosis. We present an atypical case of vivax malaria with a review of morpho-variations of iRBCs and have summarized the characteristics that aid in creating increased awareness among laboratory health professionals and public health workers.

Retrospective analysis of Plasmodium vivax genomes from a pre-elimination China inland population in the 2010s

Introduction

In malaria-free countries, imported cases are challenging because interconnections with neighboring countries with higher transmission rates increase the risk of parasite reintroduction. Establishing a genetic database for rapidly identifying malaria importation or reintroduction is crucial in addressing these challenges. This study aimed to examine genomic epidemiology during the pre-elimination stage by retrospectively reporting whole-genome sequence variation of 10 Plasmodium vivax isolates from inland China.

Methods

The samples were collected during the last few inland outbreaks from 2011 to 2012 when China implemented a malaria control plan. After next-generation sequencing, we completed a genetic analysis of the population, explored the geographic specificity of the samples, and examined clustering of selection pressures. We also scanned genes for signals of positive selection.

Results

China's inland populations were highly structured compared to the surrounding area, with a single potential ancestor. Additionally, we identified genes under selection and evaluated the selection pressure on drug-resistance genes. In the inland population, positive selection was detected in some critical gene families, including sera, msp3, and vir. Meanwhile, we identified selection signatures in drug resistance, such as ugt, krs1, and crt, and noticed that the ratio of wild-type dhps and dhfr-ts increased after China banned sulfadoxine-pyrimethamine (SP) for decades.

Discussion

Our data provides an opportunity to investigate the molecular epidemiology of pre-elimination inland malaria populations, which exhibited lower selection pressure on invasion and immune evasion genes than neighbouring areas, but increased drug resistance in low transmission settings. Our results revealed that the inland population was severely fragmented with low relatedness among infections, despite a higher incidence of multiclonal infections, suggesting that superinfection or co-transmission events are rare in low-endemic circumstances. We identified selective signatures of resistance and found that the proportion of susceptible isolates fluctuated in response to the prohibition of specific drugs. This finding is consistent with the alterations in medication strategies during the malaria elimination campaign in inland China. Such findings could provide a genetic basis for future population studies, assessing changes in other pre-elimination countries.

A genetically modified Plasmodium berghei parasite as a surrogate for whole-sporozoite vaccination against P. vivax malaria

Two malaria parasite species, Plasmodium falciparum (Pf) and P. vivax (Pv) are responsible for most of the disease burden caused by malaria. Vaccine development against this disease has focused mainly on Pf. Whole-sporozoite (WSp) vaccination, targeting pre-erythrocytic (PE) parasite stages, is a promising strategy for immunization against malaria and several PfWSp-based vaccine candidates are currently undergoing clinical evaluation. In contrast, no WSp candidates have been developed for Pv, mainly due to constraints in the production of Pv sporozoites in the laboratory. Recently, we developed a novel approach for WSp vaccination against Pf based on the use of transgenic rodent P. berghei (Pb) sporozoites expressing immunogens of this human-infective parasite. We showed that this platform can be used to deliver PE Pf antigens, eliciting both targeted humoral responses and cross-species cellular immune responses against Pf. Here we explored this WSp platform for the delivery of Pv antigens. As the Pv circumsporozoite protein (CSP) is a leading vaccine candidate antigen, we generated a transgenic Pb parasite, PbviVac, that, in addition to its endogenous PbCSP, expresses PvCSP under the control of a strictly PE promoter. Immunofluorescence microscopy analyses confirmed that both the PbCSP and the PvCSP antigens are expressed in PbviVac sporozoites and liver stages and that PbviVac sporozoite infectivity of hepatic cells is similar to that of its wild-type Pb counterpart. Immunization of mice with PbviVac sporozoites elicits the production of anti-PvCSP antibodies that efficiently recognize and bind to Pv sporozoites. Our results warrant further development and evaluation of PbviVac as a surrogate for WSp vaccination against Pv malaria.

Long read single cell RNA sequencing reveals the isoform diversity of Plasmodium vivax transcripts

Plasmodium vivax infections often consist of heterogenous populations of parasites at different developmental stages and with distinct transcriptional profiles, which complicates gene expression analyses. The advent of single cell RNA sequencing (scRNA-seq) enabled disentangling this complexity and has provided robust and stage-specific characterization of Plasmodium gene expression. However, scRNA-seq information is typically derived from the end of each mRNA molecule (usually the 3'-end) and therefore fails to capture the diversity in transcript isoforms documented in bulk RNA-seq data. Here, we describe the sequencing of scRNA-seq libraries using Pacific Biosciences (PacBio) chemistry to characterize full-length Plasmodium vivax transcripts from single cell parasites. Our results show that many P. vivax genes are transcribed into multiple isoforms, primarily through variations in untranslated region (UTR) length or splicing, and that the expression of many isoforms is developmentally regulated. Our findings demonstrate that long read sequencing can be used to characterize mRNA molecules at the single cell level and provides an additional resource to better understand the regulation of gene expression throughout the Plasmodium life cycle.

The Impact, Emerging Needs, and New Research Questions Arising from 12 Years of the Center for the Study of Complex Malaria in India

The Center for the Study of Complex Malaria in India (CSCMi) was launched in 2010 with the overall goal of addressing major gaps in our understanding of "complex malaria" in India through projects on the epidemiology, transmission, and pathogenesis of the disease. The Center was mandated to adopt an integrated approach to malaria research, including building capacity, developing infrastructure, and nurturing future malaria leaders while conducting relevant and impactful studies to assist India as it moves from control to elimination. Here, we will outline some of the interactions and impacts the Center has had with malaria policy and control counterparts in India, as well as describe emerging needs and new research questions that have become apparent over the past 12 years.

Population genomics in neglected malaria parasites

Malaria elimination includes neglected human malaria parasites Plasmodium vivax, Plasmodium ovale spp., and Plasmodium malariae. Biological features such as association with low-density infection and the formation of hypnozoites responsible for relapse make their elimination challenging. Studies on these parasites rely primarily on clinical samples due to the lack of long-term culture techniques. With improved methods to enrich parasite DNA from clinical samples, whole-genome sequencing of the neglected malaria parasites has gained increasing popularity. Population genomics of more than 2200 P. vivax global isolates has improved our knowledge of parasite biology and host-parasite interactions, identified vaccine targets and potential drug resistance markers, and provided a new way to track parasite migration and introduction and monitor the evolutionary response of local populations to elimination efforts. Here, we review advances in population genomics for neglected malaria parasites, discuss how the rich genomic information is being used to understand parasite biology and epidemiology, and explore opportunities for the applications of malaria genomic data in malaria elimination practice.

Plasmodium vivax Duffy Binding Protein-Based Vaccine: a Distant Dream

The neglected but highly prevalent Plasmodium vivax in South-east Asia and South America poses a great challenge, with regards to long-term in-vitro culturing and heavily limited functional assays. Such visible challenges as well as narrowed progress in development of experimental research tools hinders development of new drugs and vaccines. The leading vaccine candidate antigen Plasmodium vivax Duffy Binding Protein (PvDBP), is essential for reticulocyte invasion by binding to its cognate receptor, the Duffy Antigen Receptor for Chemokines (DARC), on the host’s reticulocyte surface. Despite its highly polymorphic nature, the amino-terminal cysteine-rich region II of PvDBP (PvDBPII) has been considered as an attractive target for vaccine-mediated immunity and has successfully completed the clinical trial Phase 1. Although this molecule is an attractive vaccine candidate against vivax malaria, there is still a question on its viability due to recent findings, suggesting that there are still some aspects which needs to be looked into further. The highly polymorphic nature of PvDBPII and strain-specific immunity due to PvDBPII allelic variation in Bc epitopes may complicate vaccine efficacy. Emergence of various blood-stage antigens, such as PvRBP, PvEBP and supposedly many more might stand in the way of attaining full protection from PvDBPII. As a result, there is an urgent need to assess and re-assess various caveats connected to PvDBP, which might help in designing a long-term promising vaccine for P. vivax malaria. This review mainly deals with a bunch of rising concerns for validation of DBPII as a vaccine candidate antigen for P. vivax malaria.

Comparative genomic analysis of Babesia duncani responsible for human babesiosis

Background

Human babesiosis, caused by parasites of the genus Babesia, is an emerging and re-emerging tick-borne disease that is mainly transmitted by tick bites and infected blood transfusion. Babesia duncani has caused majority of human babesiosis in Canada; however, limited data are available to correlate its genomic information and biological features.

Results

We generated a B. duncani reference genome using Oxford Nanopore Technology (ONT) and Illumina sequencing technology and uncovered its biological features and phylogenetic relationship with other Apicomplexa parasites. Phylogenetic analyses revealed that B. duncani form a clade distinct from B. microti, Babesia spp. infective to bovine and ovine species, and Theileria spp. infective to bovines. We identified the largest species-specific gene family that could be applied as diagnostic markers for this pathogen. In addition, two gene families show signals of significant expansion and several genes that present signatures of positive selection in B. duncani, suggesting their possible roles in the capability of this parasite to infect humans or tick vectors.

Conclusions

Using ONT sequencing and Illumina sequencing technologies, we provide the first B. duncani reference genome and confirm that B. duncani forms a phylogenetically distinct clade from other Piroplasm parasites. Comparative genomic analyses show that two gene families are significantly expanded in B. duncani and may play important roles in host cell invasion and virulence of B. duncani. Our study provides basic information for further exploring B. duncani features, such as host-parasite and tick-parasite interactions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01361-9.

Whole Genome Sequencing Contributions and Challenges in Disease Reduction Focused on Malaria

Malaria elimination remains an important goal that requires the adoption of sophisticated science and management strategies in the era of the COVID-19 pandemic. The advent of next generation sequencing (NGS) is making whole genome sequencing (WGS) a standard today in the field of life sciences, as PCR genotyping and targeted sequencing provide insufficient information compared to the whole genome. Thus, adapting WGS approaches to malaria parasites is pertinent to studying the epidemiology of the disease, as different regions are at different phases in their malaria elimination agenda. Therefore, this review highlights the applications of WGS in disease management, challenges of WGS in controlling malaria parasites, and in furtherance, provides the roles of WGS in pursuit of malaria reduction and elimination. WGS has invaluable impacts in malaria research and has helped countries to reach elimination phase rapidly by providing required information needed to thwart transmission, pathology, and drug resistance. However, to eliminate malaria in sub-Saharan Africa (SSA), with high malaria transmission, we recommend that WGS machines should be readily available and affordable in the region.

Primate malarias as a model for cross-species parasite transmission

Parasites regularly switch into new host species, representing a disease burden and conservation risk to the hosts. The distribution of these parasites also gives insight into characteristics of ecological networks and genetic mechanisms of host-parasite interactions. Some parasites are shared across many species, whereas others tend to be restricted to hosts from a single species. Understanding the mechanisms producing this distribution of host specificity can enable more effective interventions and potentially identify genetic targets for vaccines or therapies. As ecological connections between human and local animal populations increase, the risk to human and wildlife health from novel parasites also increases. Which of these parasites will fizzle out and which have the potential to become widespread in humans? We consider the case of primate malarias, caused by Plasmodium parasites, to investigate the interacting ecological and evolutionary mechanisms that put human and nonhuman primates at risk for infection. Plasmodium host switching from nonhuman primates to humans led to ancient introductions of the most common malaria-causing agents in humans today, and new parasite switching is a growing threat, especially in Asia and South America. Based on a wild host-Plasmodium occurrence database, we highlight geographic areas of concern and potential areas to target further sampling. We also discuss methodological developments that will facilitate clinical and field-based interventions to improve human and wildlife health based on this eco-evolutionary perspective.

Spatiotemporal Changes in Plasmodium vivax msp142 Haplotypes in Southern Mexico: From the Control to the Pre-Elimination Phase

For 20 years, Plasmodium vivax has been the only prevalent malaria species in Mexico, and cases have declined significantly and continuously. Spatiotemporal genetic studies can be helpful for understanding parasite dynamics and developing strategies to weaken malaria transmission, thus facilitating the elimination of the parasite. The aim of the current contribution was to analyze P. vivax-infected blood samples from patients in southern Mexico during the control (1993–2007) and pre-elimination phases (2008–2011). Nucleotide and haplotype changes in the pvmsp142 fragment were evaluated over time. The majority of multiple genotype infections occurred in the 1990s, when the 198 single nucleotide sequences exhibited 57 segregating sites, 64 mutations, and 17 haplotypes. Nucleotide and genetic diversity parameters showed subtle fluctuations from across time, in contrast to the reduced haplotype diversity and the increase in the R2 index and Tajima’s D value from 2008 to 2011. The haplotype network consisted of four haplogroups, the geographical distribution of which varied slightly over time. Haplogroup-specific B-cell epitopes were predicted. Since only high-frequency and divergent haplotypes persisted, there was a contraction of the parasite population. Given that 84% of haplotypes were exclusive to Mesoamerica, P. vivax flow is likely circumscribed to this region, representing important information for parasite surveillance.

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.

Progress in understanding the phylogeny of the Plasmodium vivax lineage

There has been some controversy about the evolutionary origin of Plasmodium vivax, particularly whether it is of Asian or African origin. Recently, a new malaria species which closely related to ape P. vivax was found in chimpanzees, in addition, the host switches of P. vivax from ape to human was confirmed. These findings support the African origin of P. vivax. Previous phylogenetic analyses have shown the position of P. vivax within the Asian primate malaria parasite clade. This suggested an Asian origin of P. vivax. Recent analyses using massive gene data, however, positioned P. vivax after the branching of the African Old World monkey parasite P. gonderi, and before the branching of the common ancestor of Asian primate malaria parasites. This position is consistent with an African origin of P. vivax. We here review the history of phylogenetic analyses on P. vivax, validate previous analyses, and finally present a definitive analysis using currently available data that indicate a tree in which P. vivax is positioned at the base of the Asian primate malaria parasite clade, and thus that is consistent with an African origin of P. vivax.

Genetic Diversity of Plasmodium vivax Cysteine-Rich Protective Antigen (PvCyRPA) in Field Isolates from Five Different Areas of the Brazilian Amazon

The Plasmodium vivax Cysteine-Rich Protective Antigen (PvCyRPA) has an important role in erythrocyte invasion and has been considered a target for vivax malaria vaccine development. Nonetheless, its genetic diversity remains uncharted in Brazilian malaria-endemic areas. Therefore, we investigated the pvcyrpa genetic polymorphism in 98 field isolates from the Brazilian Amazon and its impact on the antigenicity of predicted B-cell epitopes. Genetic diversity parameters, population genetic analysis, neutrality test and the median-joining network were analyzed, and the potential amino acid polymorphism participation in B-cell epitopes was investigated. One synonymous and 26 non-synonymous substitutions defined fifty haplotypes. The nucleotide diversity and Tajima’s D values varied across the coding gene. The exon-1 sequence had greater diversity than those of exon-2. Concerning the prediction analysis, seven sequences were predicted as linear B cell epitopes, the majority contained in conformational epitopes. Moreover, important amino acid polymorphism was detected in regions predicted to contain residues participating in B-cell epitopes. Our data suggest that the pvcyrpa gene presents a moderate polymorphism in the studied isolates and such polymorphisms alter amino acid sequences contained in potential B cell epitopes, an important observation considering the antigen potentiality as a vaccine candidate to cover distinct P. vivax endemic areas worldwide.

Genetic diversity in two leading Plasmodium vivax malaria vaccine candidates AMA1 and MSP119 at three sites in India

Plasmodium vivax, a major contributor to the malaria burden in India, has the broadest geographic distribution and shows higher genetic diversity than P. falciparum. Here, we investigated the genetic diversity of two leading P. vivax vaccine candidate antigens, at three geographically diverse malaria-endemic regions in India. Pvama1 and Pvmsp119 partial coding sequences were generated from one hundred P. vivax isolates in India (Chennai n = 28, Nadiad n = 50 and Rourkela n = 22) and ~1100 published sequences from Asia, South America, North America, and Oceania regions included. These data were used to assess the genetic diversity and potential for vaccine candidacy of both antigens on a global scale. A total of 44 single nucleotide polymorphism (SNPs) were identified among 100 Indian Pvama1 sequences, including 10 synonymous and 34 nonsynonymous mutations. Nucleotide diversity was higher in Rourkela and Nadiad as compared to Chennai. Nucleotide diversity measures showed a strong balancing selection in Indian and global population for domain I of Pvama1, which suggests that it is a dominant target of the protective immune response. In contrast, the Pvmsp119 region showed highly conserved sequences in India and across the Oceania, South America, North America and Asia, demonstrating low genetic diversity in the global population when compared to Pvama1. Results suggest the possibility of including Pvmsp119 in a multivalent vaccine formulation against P. vivax infections. However, the high genetic diversity seen in Pvama1 would be more challenging for vaccine development.

Advances and opportunities in malaria population genomics

Almost 20 years have passed since the first reference genome assemblies were published for Plasmodium falciparum, the deadliest malaria parasite, and Anopheles gambiae, the most important mosquito vector of malaria in sub-Saharan Africa. Reference genomes now exist for all human malaria parasites and nearly half of the ~40 important vectors around the world. As a foundation for genetic diversity studies, these reference genomes have helped advance our understanding of basic disease biology and drug and insecticide resistance, and have informed vaccine development efforts. Population genomic data are increasingly being used to guide our understanding of malaria epidemiology, for example by assessing connectivity between populations and the efficacy of parasite and vector interventions. The potential value of these applications to malaria control strategies, together with the increasing diversity of genomic data types and contexts in which data are being generated, raise both opportunities and challenges in the field. This Review discusses advances in malaria genomics and explores how population genomic data could be harnessed to further support global disease control efforts.

Variation in selective constraints along the Plasmodium life cycle

Plasmodium parasites, the cause of malaria, have a complex life cycle, infecting alternatively vertebrate hosts and female Anopheles mosquitoes and undergoing intra- and extra-cellular development in several organs of these hosts. Most of the ~5000 protein-coding genes present in Plasmodium genomes are only expressed at specific life stages, and different genes might therefore be subject to different selective pressures depending on the biological activity of the parasite and its microenvironment at this point in development. Here, we estimate the selective constraints on the protein-coding sequences of all annotated genes of rodent and primate Plasmodium parasites and, using data from scRNA-seq experiments spanning many developmental stages, analyze their variation with regard to when these genes are expressed in the parasite life cycle. Our study reveals extensive variation in selective constraints throughout the parasites' development and highlights stages that are evolving more rapidly than others. These findings provide novel insights into the biology of these parasites and could provide important information to develop better treatment strategies or vaccines against these medically-important organisms.

Sequence Analysis of Pvama-1 among Plasmodium Vivax Isolates in Sistan-Baluchistan

Background

Apical Membrane antigen 1 (AMA-1) is an important membrane protein that presents in all Plasmodium species and participates in critical phases in the attraction of cells. In human, it is one of the most immunodominant antigens with a protective immune response simulation role Apical Membrane antigen 1 (AMA-1) is an important membrane protein which presents in all Plasmodium species and is located on the surface of merozoite and sporozoites that participates in critical phases in attraction of human red blood cells by merozoites and hepatocytes by sporozoites, so in human, it is one of the most immunodominant antigens with a protective immune response simulation role. Since extra information is necessary to lighten of AMA-1 scope, we equaled genetic variation in P.vivax AMA-1 from 40 Iranian isolates with those reported from the other malarious countries.

Methods

Blood samples were collected from 40 patients' positive of P.vivax, and genomic DNA was extracted from the blood. The nucleotide sequence for 446 amino acid (AA) residues (42-488 of PvAMA-1) of AMA-1 gene was amplified via PCR and then sequenced.

Result

A total of 24 different haplotypes were recognized between samples. No new haplotype was determined in this research that was reported previously in other regions of Iran and the world. We detected 37-point mutations at the nucleotide level in their sequences and showed 43 amino acid variations, at 37 positions in which 6 sites demonstrate trimorphic polymorphism, and the others were dimorphic.

Conclusion

Sequence analysis of the major haplotype showed 95% similarity with P.vivax Sal-1 AMA-1 gene and high level of allelic diversity at the domain I of PvAMA-1 among P. vivax isolates of Iran. Because PvAMA-1 is noticeable as vaccine candidate antigen, these documents provide valuable information for the development of malaria vaccine.

Ape Origins of Human Malaria

African apes harbor at least twelve Plasmodium species, some of which have been a source of human infection. It is now well established that Plasmodium falciparum emerged following the transmission of a gorilla parasite, perhaps within the last 10,000 years, while Plasmodium vivax emerged earlier from a parasite lineage that infected humans and apes in Africa before the Duffy-negative mutation eliminated the parasite from humans there. Compared to their ape relatives, both human parasites have greatly reduced genetic diversity and an excess of nonsynonymous mutations, consistent with severe genetic bottlenecks followed by rapid population expansion. A putative new Plasmodium species widespread in chimpanzees, gorillas, and bonobos places the origin of Plasmodium malariae in Africa. Here, we review what is known about the origins and evolutionary history of all human-infective Plasmodium species, the time and circumstances of their emergence, and the diversity, host specificity, and zoonotic potential of their ape counterparts.

Malaria in the 'Omics Era'

Genomics has revolutionised the study of the biology of parasitic diseases. The first Eukaryotic parasite to have its genome sequenced was the malaria parasite Plasmodium falciparum. Since then, Plasmodium genomics has continued to lead the way in the study of the genome biology of parasites, both in breadth—the number of Plasmodium species’ genomes sequenced—and in depth—massive-scale genome re-sequencing of several key species. Here, we review some of the insights into the biology, evolution and population genetics of Plasmodium gained from genome sequencing, and look at potential new avenues in the future genome-scale study of its biology.

Population genomic evidence of Plasmodium vivax Southeast Asian origin

Plasmodium vivax is the most common and widespread human malaria parasite. It was recently proposed that P. vivax originates from sub-Saharan Africa based on the circulation of its closest genetic relatives (P. vivax-like) among African great apes. However, the limited number of genetic markers and samples investigated questions the robustness of this hypothesis. Here, we extensively characterized the genomic variations of 447 human P. vivax strains and 19 ape P. vivax-like strains collected worldwide. Phylogenetic relationships between human and ape Plasmodium strains revealed that P. vivax is a sister clade of P. vivax-like, not included within the radiation of P. vivax-like By investigating various aspects of P. vivax genetic variation, we identified several notable geographical patterns in summary statistics in function of the increasing geographic distance from Southeast Asia, suggesting that P. vivax may have derived from a single area in Asia through serial founder effects.

Next Generation Sequencing and Bioinformatics Methodologies for Infectious Disease Research and Public Health: Approaches, Applications, and Considerations for Development of Laboratory Capacity

Next generation sequencing (NGS) combined with bioinformatics has successfully been used in a vast array of analyses for infectious disease research of public health relevance. For instance, NGS and bioinformatics approaches have been used to identify outbreak origins, track transmissions, investigate epidemic dynamics, determine etiological agents of a disease, and discover novel human pathogens. However, implementation of high-quality NGS and bioinformatics in research and public health laboratories can be challenging. These challenges mainly include the choice of the sequencing platform and the sequencing approach, the choice of bioinformatics methodologies, access to the appropriate computation and information technology infrastructure, and recruiting and retaining personnel with the specialized skills and experience in this field. In this review, we summarize the most common NGS and bioinformatics workflows in the context of infectious disease genomic surveillance and pathogen discovery, and highlight the main challenges and considerations for setting up an NGS and bioinformatics-focused infectious disease research public health laboratory. We describe the most commonly used sequencing platforms and review their strengths and weaknesses. We review sequencing approaches that have been used for various pathogens and study questions, as well as the most common difficulties associated with these approaches that should be considered when implementing in a public health or research setting. In addition, we provide a review of some common bioinformatics tools and procedures used for pathogen discovery and genome assembly, along with the most common challenges and solutions. Finally, we summarize the bioinformatics of advanced viral, bacterial, and parasite pathogen characterization, including types of study questions that can be answered when utilizing NGS and bioinformatics.

Progress towards the development of a P. vivax vaccine

INTRODUCTION

Plasmodium vivax causes significant public health problems in endemic regions. A vaccine to prevent disease is critical, considering the rapid spread of drug-resistant parasite strains, and the development of hypnozoites in the liver with potential for relapse. A minimally effective vaccine should prevent disease and transmission while an ideal vaccine provides sterile immunity.

AREAS COVERED

Despite decades of research, the complex life cycle, technical challenges and a lack of funding have hampered progress of P. vivax vaccine development. Here, we review the progress of potential P. vivax vaccine candidates from different stages of the parasite life cycle. We also highlight the challenges and important strategies for rational vaccine design. These factors can significantly increase immune effector mechanisms and improve the protective efficacy of these candidates in clinical trials to generate sustained protection over longer periods of time.

EXPERT OPINION

A vaccine that presents functionally-conserved epitopes from multiple antigens from various stages of the parasite life cycle is key to induce broadly neutralizing strain-transcending protective immunity to effectively disrupt parasite development and transmission.

Human Plasmodium vivax mosquito experimental transmission

Plasmodium vivax bench research greatly lags behind Plasmodium falciparum because of an inability to culture in vitro. A century ago, intentionally inducing a malaria infection was a strategy commonly used to cure late-stage syphilis. These controlled human malaria infections were used with expertise and persisted to the end of World War II. While controlled malaria liver-stage infection has been achieved for both P. vivax and P. falciparum, controlled human transmission to mosquitoes falls short for both species. In this issue of the JCI, Collins et al. present groundbreaking work that establishes a system to transmit P. vivax gametocytes from humans to mosquitoes. The authors injected a unique human isolate of P. vivax that reached high gametocyte density within weeks. This study provides a technical advance that will facilitate the study and eradication of the human parasite P. vivax.

Zoonotic Transmissions and Host Switches of Malaria Parasites

Malaria is a deadly disease that affects the health of hundreds of millions of people annually. There are five Plasmodium parasite species that can naturally infect humans, including Plasmodium falciparum, Plasmodium vivax, Plasmodium malariae, Plasmodium ovale and Plasmodium knowlesi. Some of the parasites can also infect various non-human primates. Parasites mainly infecting monkeys such as Plasmodium cynomolgi (in fact P. knowlesi was considered as a parasite of monkeys for years) can also be transmitted to human hosts. Recently, many new Plasmodium species were discovered in African apes, and it is possible that some of the parasites can be transmitted to humans in the future. Here, we searched PubMed and the internet via Google and selected articles concerning zoonotic transmission and evolution of selected malaria parasite species. We reviewed the current advances in the relevant topics emphasizing on transmissions of malaria parasites between humans and non-human primates. We also briefly discuss the transmissions of some avian malaria parasites between wild birds and domestic fowls. Zoonotic malaria transmissions are widespread, which poses a threat to public health. More studies on parasite species identification in non-human primates, transmission, and evolution are needed to reduce or prevent transmission of malaria parasites from non-human primates to humans.

Antigenicity and adhesiveness of a Plasmodium vivax VIR-E protein from Brazilian isolates

BACKGROUND

Plasmodium vivax, the major cause of malaria in Latin America, has a large subtelomeric multigene family called vir. In the P. vivax genome, about 20% of its sequences are vir genes. Vir antigens are grouped in subfamilies according to their sequence similarities and have been shown to have distinct roles and subcellular locations. However, little is known about vir subfamilies, especially when comes to their functions.

OBJECTIVE

To evaluate the diversity, antigenicity, and adhesiveness of Plasmodium vivax VIR-E.

METHODS

Vir-E genes were amplified from six P. vivax isolates from Manaus, North of Brazil. The presence of naturally acquired antibodies to recombinant PvBrVIR-E and PvAMA-1 was evaluated by ELISA. Binding capacity of recombinant PvBrVIR-E was assessed by adhesion assay to CHO-ICAM1 cells.

FINDINGS

Despite vir-E sequence diversity, among those identified sequences, a representative one was chosen to be expressed as recombinant protein. The presence of IgM or IgG antibodies to PvBrVIR-E was detected in 23.75% of the study population while the presence of IgG antibodies to PvAMA-1 antigen was 66.25% in the same population. PvBrVIR-E was adhesive to CHO-ICAM1.

MAIN CONCLUSIONS

PvBrVIR-E was antigenic and adhesive to CHO-ICAM1.

Monitoring Plasmodium falciparum and Plasmodium vivax using microsatellite markers indicates limited changes in population structure after substantial transmission decline in Papua New Guinea

Monitoring the genetic structure of pathogen populations may be an economical and sensitive approach to quantify the impact of control on transmission dynamics, highlighting the need for a better understanding of changes in population genetic parameters as transmission declines. Here we describe the first population genetic analysis of two major human malaria parasites, Plasmodium falciparum (Pf) and Plasmodium vivax (Pv), following nationwide distribution of long-lasting insecticide-treated nets (LLINs) in Papua New Guinea (PNG). Parasite isolates from pre- (2005-2006) and post-LLIN (2010-2014) were genotyped using microsatellite markers. Despite parasite prevalence declining substantially (East Sepik Province: Pf = 54.9%-8.5%, Pv = 35.7%-5.6%, Madang Province: Pf = 38.0%-9.0%, Pv: 31.8%-19.7%), genetically diverse and intermixing parasite populations remained. Pf diversity declined modestly post-LLIN relative to pre-LLIN (East Sepik: R_s  = 7.1-6.4, H_E  = 0.77-0.71; Madang: R_s  = 8.2-6.1, H_E  = 0.79-0.71). Unexpectedly, population structure present in pre-LLIN populations was lost post-LLIN, suggesting that more frequent human movement between provinces may have contributed to higher gene flow. Pv prevalence initially declined but increased again in one province, yet diversity remained high throughout the study period (East Sepik: R_s  = 11.4-9.3, H_E  = 0.83-0.80; Madang: R_s  = 12.2-14.5, H_E  = 0.85-0.88). Although genetic differentiation values increased between provinces over time, no significant population structure was observed at any time point. For both species, a decline in multiple infections and increasing clonal transmission and significant multilocus linkage disequilibrium post-LLIN were positive indicators of impact on the parasite population using microsatellite markers. These parameters may be useful adjuncts to traditional epidemiological tools in the early stages of transmission reduction.

Balancing selection and high genetic diversity of Plasmodium vivax circumsporozoite central region in parasites from Brazilian Amazon and Rio de Janeiro Atlantic Forest

Circumsporozoite protein (CSP) is the primary pre-erythrocytic vaccine target in Plasmodium species. Knowledge about their genetic diversity can help predict vaccine efficacy and the spread of novel parasite variants. Thus, we investigated pvcsp gene polymorphisms in 219 isolates (136 from Brazilian Amazon [BA], 71 from Rio de Janeiro Atlantic Forest [AF], and 12 from non-Brazilian countries [NB]). Forty-eight polymorphic sites were detected, 46 in the central repeat region (CR), and two in the C-terminal region. Also, the CR presents InDels and a variable number of repeats. All samples correspond to the VK210 variant, and 24 VK210 subtypes based on CR. Nucleotide diversity (π = 0.0135) generated a significant number of haplotypes (168) with low genetic differentiation between the Brazilian regions (Fst = 0.208). The haplotype network revealed similar distances among the BA and AF regions. The linkage disequilibrium indicates that recombination does not seem to be acting in diversity, reinforcing natural selection's role in accelerating adaptive evolution. The high diversity (low Fst) and polymorphism frequencies could be indicators of balancing selection. Although malaria in BA and AF have distinct vector species and different host immune pressures, consistent genetic signature was found in two regions. The immunodominant B-cell epitope mapped in the CR varies from seven to 19 repeats. The CR T-cell epitope is conserved only in 39 samples. Concerning to C-terminal region, the Th2R epitope presented nonsynonymous SNP only in 6% of Brazilian samples, and the Th3R epitope remained conserved in all studied regions. We conclude that, although the uneven distribution of alleles may jeopardize the deployment of vaccines directed to a specific variable locus, a unique vaccine formulation could protect populations in all Brazilian regions.

Exploration of Plasmodium vivax merozoite surface proteins 1 and 7 genetic diversity in Brazilian Amazon and Rio de Janeiro Atlantic Forest

Plasmodium vivax merozoite surface proteins (PvMSP) 1 and 7 are considered vaccine targets. Genetic diversity knowledge is crucial to assess their potential as immunogens and to provide insights about population structure in different epidemiological contexts. Here, we investigate the variability of pvmsp-1₄₂, pvmsp-7E, and pvmsp-7F genes in 227 samples from the Brazilian Amazon (BA) and Rio de Janeiro Atlantic Forest (AF). pvmsp-1₄₂ has 63 polymorphisms - 57 nonsynonymous - generating a nucleotide diversity of π = 0.009 in AF, and π = 0.018 in BA. In pvmsp-7E, 134 polymorphisms - 103 nonsynonymous - generate the nucleotide diversity of π = 0.027 in AF, and π = 0.042 in BA. The pvmsp-7F has only two SNPs - A610G and A1054T -, with nucleotide diversity of π = 0.0004 in AF, and π = 0.0007 in BA. The haplotype diversity of pvmsp-1₄₂, pvmsp-7E, and pvmsp-7F genes is 0.997, 1.00, and 0.649, respectively. None of the pvmsp-1₄₂ or pvmsp-7E sequences are identical to the Salvador 1 strain's sequence. Conversely, most of pvmsp-7F sequences (94/48%) are identical to Sal-1. We evaluated eight B-cell epitopes in pvmsp-7E, four of them showed higher nucleotide diversity compared to pvmsp-7E's epitopes. Positive selection was detected in pvmsp-1₄₂, pvmsp-7E central region, and pvmsp-7F with Tajima's D. In pvmsp-7E, the significant nucleotide and haplotype diversities with low genetic differentiation, could be indicative of balancing selection. The genetic differentiation of pvmsp-1₄₂ (0.315) and pvmsp-7F (0.354) genes between AF and BA regions is significant, which is not the case for pvmsp-7E (0.193). We conclude that pvmsp-1₄₂ and pvmsp-7E have great genetic diversity even in AF region, an enclosure area with deficient transmission levels of P. vivax zoonotic malaria. In both Brazilian regions, pvmsp-1₁₉, pvmsp-7E, and pvmsp-7F are conserved, most likely due to their roles in parasite survival, and could be considered potential targets for a "blood-stage vaccine".

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.

Population genomics identifies a distinct Plasmodium vivax population on the China-Myanmar border of Southeast Asia

Plasmodium vivax has become the predominant malaria parasite and a major challenge for malaria elimination in the Greater Mekong Subregion (GMS). Yet, our knowledge about the evolution of P. vivax populations in the GMS is fragmental. We performed whole genome sequencing on 23 P. vivax samples from the China-Myanmar border (CMB) and used 21 high-coverage samples to compare to over 200 samples from the rest of the GMS. Using genome-wide single nucleotide polymorphisms (SNPs), we analyzed population differentiation, genetic structure, migration and potential selection using an array of methods. The CMB parasites displayed a higher proportion of monoclonal infections, and 52% shared over 90% of their genomes in identity-by-descent segments with at least one other sample from the CMB, suggesting preferential expansion of certain parasite strains in this region, likely resulting from the P. vivax outbreaks occurring during this study period. Principal component, admixture, fixation index and phylogenetic analyses all identified that parasites from the CMB were genetically distinct from parasites from eastern parts of the GMS (Cambodia, Laos, Vietnam, and Thailand), whereas the eastern GMS parasite populations were largely undifferentiated. Such a genetic differentiation pattern of the P. vivax populations from the GMS parasite was largely explainable through geographic distance. Using the genome-wide SNPs, we narrowed down to a set of 36 SNPs for differentiating parasites from different areas of the GMS. Genome-wide scans to determine selection in the genome with two statistical methods identified genes potentially under drug selection, including genes associated with antifolate resistance and genes linked to chloroquine resistance in Plasmodium falciparum.

Kerteszia cruzii and extra-Amazonian malaria in Brazil: Challenges due to climate change in the Atlantic Forest

Kerteszia cruzii is a sylvatic mosquito and the primary vector of Plasmodium spp., which can cause malaria in humans in areas outside the Amazon River basin in Brazil. Anthropic changes in the natural environments are the major drivers of massive deforestation and local climate change, with serious impacts on the dynamics of mosquito communities and on the risk of acquiring malaria. Considering the lack of information on the dynamics of malaria transmission in areas across the Atlantic Forest biome, where Ke. cruzii is the dominant vector, and the impact of climate drivers of malaria, the present study aimed to: (i) investigate the occurrence and survival rate of Ke. cruzii based on the distinct vegetation profiles found in areas across the coastal region of the Brazilian Atlantic Forest biome; (ii) estimate the extrinsic incubation period (EIP) and survival rates of P. vivax and P. falciparum parasites in Ke. cruzii under current and future scenarios. The potential distribution of Plasmodium spp. was estimated using simulation analyses under distinct scenarios of average temperature increases from 1 °C to 3.7 °C. Our results showed that two conditions are necessary to explain the occurrence and survival of Ke. cruzii: warm temperature and presence of the Atlantic Forest biome. Moreover, both Plasmodium species showed a tendency to decrease their EIP and increase their estimated survival rates in a scenario of higher temperature. Our findings support that the high-risk malaria areas may include the southern region of the distribution range of the Atlantic Forest biome in the coming years. Despite its limitations and assumptions, the present study provides robust evidence of areas with potential to be impacted by malaria incidence in a future scenario. These areas should be monitored in the next decades regarding the occurrence of the mosquito vector and the potential for malaria persistence and increased occurrence.

Potential and Limitations of Cross-Protective Vaccine against Malaria by Blood-Stage Naturally Attenuated Parasite

Human malaria vaccine trials have revealed vaccine efficacy but improvement is still needed. In this study, we aimed to re-evaluate vaccination with blood-stage naturally attenuated parasites, as a whole-organism vaccine model against cross-strain and cross-species malaria, to establish a better vaccination strategy. C57BL/6 mice controlled blood-stage Plasmodium yoelii 17XNL (PyNL) within 1 month of infection, while mice with a variety of immunodeficiencies demonstrated different susceptibilities to PyNL, including succumbing to hyperparasitemia. However, after recovery, survivors had complete protection against a challenge with the lethal strain PyL. Unlike cross-strain protection, PyNL-recovered mice failed to induce sterile immunity against Plasmodium berghei ANKA, although prolonged survival was observed in some vaccinated mice. Splenomegaly is a typical characteristic of malaria; the splenic structure became reorganized to prioritize extra-medullary hematopoiesis and to eliminate parasites. We also found that the peritoneal lymph node was enlarged, containing activated/memory phenotype cells that did not confer protection against PyL challenge. Hemozoins remained in the spleen several months after PyNL infection. Generation of an attenuated human blood-stage parasite expressing proteins from multiple species of malaria would greatly improve anti-malaria vaccination.

Selective whole genome amplification of Plasmodium malariae DNA from clinical samples reveals insights into population structure

The genomic diversity of Plasmodium malariae malaria parasites is understudied, partly because infected individuals tend to present with low parasite densities, leading to difficulties in obtaining sufficient parasite DNA for genome analysis. Selective whole genome amplification (SWGA) increases the relative levels of pathogen DNA in a clinical sample, but has not been adapted for P. malariae parasites. Here we design customized SWGA primers which successfully amplify P. malariae DNA extracted directly from unprocessed clinical blood samples obtained from patients with P. malariae-mono-infections from six countries, and further test the efficacy of SWGA on mixed infections with other Plasmodium spp. SWGA enables the successful whole genome sequencing of samples with low parasite density (i.e. one sample with a parasitaemia of 0.0064% resulted in 44% of the genome covered by ≥ 5 reads), leading to an average 14-fold increase in genome coverage when compared to unamplified samples. We identify a total of 868,476 genome-wide SNPs, of which 194,709 are unique across 18 high-quality isolates. After exclusion of the hypervariable subtelomeric regions, a high-quality core subset of 29,899 unique SNPs is defined. Population genetic analysis suggests that P. malariae parasites display clear geographical separation by continent. Further, SWGA successfully amplifies genetic regions of interest such as orthologs of P. falciparum drug resistance-associated loci (Pfdhfr, Pfdhps, Pfcrt, Pfk13 and Pfmdr1), and several non-synonymous SNPs were detected in these genes. In conclusion, we have established a robust SWGA approach that can assist whole genome sequencing of P. malariae, and thereby facilitate the implementation of much-needed large-scale multi-population genomic studies of this neglected malaria parasite. As demonstrated in other Plasmodia, such genetic diversity studies can provide insights into the biology underlying the disease and inform malaria surveillance and control measures.