Insights into the invasion biology of Plasmodium vivax

Population-based genomic study of Plasmodium vivax malaria in seven Brazilian states and across South America

Background

Brazil is a unique and understudied setting for malaria, with complex foci of transmission associated with human and environmental conditions. An understanding of the population genomic diversity of P. vivax parasites across Brazil can support malaria control strategies.

Methods

Through whole genome sequencing of P. vivax isolates across 7 Brazilian states, we use population genomic approaches to compare genetic diversity within country (n = 123), continent (6 countries, n = 315) and globally (26 countries, n = 885).

Findings

We confirm that South American isolates are distinct, have more ancestral populations than the other global regions, with differentiating mutations in genes under selective pressure linked to antimalarial drugs (pvmdr1, pvdhfr-ts) and mosquito vectors (pvcrmp3, pvP45/48, pvP47). We demonstrate Brazil as a distinct parasite population, with signals of selection including ABC transporter (PvABCI3) and PHIST exported proteins.

Interpretation

Brazil has a complex population structure, with evidence of P. simium infections and Amazonian parasites separating into multiple clusters. Overall, our work provides the first Brazil-wide analysis of P. vivax population structure and identifies important mutations, which can inform future research and control measures.

Funding

AI is funded by an MRC LiD PhD studentship. TGC is funded by the Medical Research Council (Grant no. MR/M01360X/1, MR/N010469/1, MR/R025576/1, MR/R020973/1 and MR/X005895/1). SC is funded by Medical Research Council UK grants (MR/M01360X/1, MR/R025576/1, MR/R020973/1 and MR/X005895/1) and Bloomsbury SET (ref. CCF17-7779). FN is funded by The Shloklo Malaria Research Unit - part of the Mahidol Oxford Research Unit, supported by the Wellcome Trust (Grant no. 220211). ARSB is funded by São Paulo Research Foundation - FAPESP (Grant no. 2002/09546-1). RLDM is funded by Brazilian National Council for Scientific and Technological Development - CNPq (Grant no. 302353/2003-8 and 471605/2011-5); CRFM is funded by FAPESP (Grant no. 2020/06747-4) and CNPq (Grant no. 302917/2019-5 and 408636/2018-1); JGD is funded by FAPESP fellowships (2016/13465-0 and 2019/12068-5) and CNPq (Grant no. 409216/2018-6).

Plasmodium malariae Detected by Microscopy in the International Bordering Area of Mizoram, a Northeastern State of India

Northeastern states of India share international borders with Myanmar, China, Bangladesh, and Bhutan, contributing 7.45% of the overall malaria cases in the country. Mizoram accounts for the highest malaria burden in the northeastern states, with perennial transmission in the hilly and deep-forested areas. Plasmodium falciparum (93%) is the most prevalent human Plasmodium species, followed by P. vivax; however, information on P. ovale and P. malariae is negligible. Rapid diagnostic tests (RDTs) are the most preferred malaria diagnostic tool followed by microscopy in this high malaria-endemic region. The present epidemiological study was carried out in July and August 2019 to assess the malaria burden in and around the Chawngte primary health center, Lawngtlai District of Mizoram, using RDTs and microscopy as diagnostic tools. World Health Organization-certified level I microscopists examined the blood smears. Diagnosis using RDTs resulted in 151 malaria cases (P. falciparum: 136; P. vivax: 15) out of 948 screened fever cases. However, blood smear examination detected 179 cases (P. falciparum: 154; P. vivax: 17; mixed P. falciparum + P. vivax infection: 3; P. malariae: 5). Analysis revealed that the risk of malaria infection was higher in the ≥5-year-old subjects than in the under-5 age group. The mean parasite density of P. malariae (1455.00/μL blood) was the lowest; cf. with P. falciparum: 12,275.08/μL blood. Surveillance at the point-of-care level using microscopy was able to detect all the four human Plasmodium species and their mixed infections, including P. malariae, which were missed with RDTs. Thus, the quality of microscopy along with trained manpower should be strengthened to diagnose all human malaria parasite species (particularly P. malariae and P. ovale) until the molecular tools are deployed at the field level to achieve malaria elimination by 2030.

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.

Self-assembling functional programmable protein array for studying protein-protein interactions in malaria parasites

Background

Plasmodium vivax is the most widespread malarial species, causing significant morbidity worldwide. Knowledge is limited regarding the molecular mechanism of invasion due to the lack of a continuous in vitro culture system for these species. Since protein–protein and host–cell interactions play an essential role in the microorganism’s invasion and replication, elucidating protein function during invasion is critical when developing more effective control methods. Nucleic acid programmable protein array (NAPPA) has thus become a suitable technology for studying protein–protein and host–protein interactions since producing proteins through the in vitro transcription/translation (IVTT) method overcomes most of the drawbacks encountered to date, such as heterologous protein production, stability and purification.

Results

Twenty P. vivax proteins on merozoite surface or in secretory organelles were selected and successfully cloned using gateway technology. Most constructs were displayed in the array expressed in situ, using the IVTT method. The Pv12 protein was used as bait for evaluating array functionality and co-expressed with P. vivax cDNA display in the array. It was found that Pv12 interacted with Pv41 (as previously described), as well as PvMSP1_42kDa, PvRBP1a, PvMSP8 and PvRAP1.

Conclusions

NAPPA is a high-performance technique enabling co-expression of bait and query in situ, thereby enabling interactions to be analysed rapidly and reproducibly. It offers a fresh alternative for studying protein–protein and ligand–receptor interactions regarding a parasite which is difficult to cultivate (i.e. P. vivax).

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2414-2) contains supplementary material, which is available to authorized users.

A large scale Plasmodium vivax- Saimiri boliviensis trophozoite-schizont transition proteome

Plasmodium vivax is a complex protozoan parasite with over 6,500 genes and stage-specific differential expression. Much of the unique biology of this pathogen remains unknown, including how it modifies and restructures the host reticulocyte. Using a recently published P. vivax reference genome, we report the proteome from two biological replicates of infected Saimiri boliviensis host reticulocytes undergoing transition from the late trophozoite to early schizont stages. Using five database search engines, we identified a total of 2000 P. vivax and 3487 S. boliviensis proteins, making this the most comprehensive P. vivax proteome to date. PlasmoDB GO-term enrichment analysis of proteins identified at least twice by a search engine highlighted core metabolic processes and molecular functions such as glycolysis, translation and protein folding, cell components such as ribosomes, proteasomes and the Golgi apparatus, and a number of vesicle and trafficking related clusters. Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 enriched functional annotation clusters of S. boliviensis proteins highlighted vesicle and trafficking-related clusters, elements of the cytoskeleton, oxidative processes and response to oxidative stress, macromolecular complexes such as the proteasome and ribosome, metabolism, translation, and cell death. Host and parasite proteins potentially involved in cell adhesion were also identified. Over 25% of the P. vivax proteins have no functional annotation; this group includes 45 VIR members of the large PIR family. A number of host and pathogen proteins contained highly oxidized or nitrated residues, extending prior trophozoite-enriched stage observations from S. boliviensis infections, and supporting the possibility of oxidative stress in relation to the disease. This proteome significantly expands the size and complexity of the known P. vivax and Saimiri host iRBC proteomes, and provides in-depth data that will be valuable for ongoing research on this parasite’s biology and pathogenesis.

Quantitative Proteomics Analysis of Plasmodium vivax Induced Alterations in Human Serum during the Acute and Convalescent Phases of Infection

The radial distribution of Plasmodium vivax malaria burden has evoked enormous concern among the global research community. In this study, we have investigated the serum proteome alterations in non-severe vivax malaria patients before and during patient recuperation starting from the early febrile to the defervescence and convalescent stages of the infection. We have also performed an extensive quantitative proteomics analysis to compare the serum proteome profiles of vivax malaria patients with low (LPVM) and moderately-high (MPVM) parasitemia with healthy community controls. Interestingly, some of the serum proteins such as Serum amyloid A, Apolipoprotein A1, C-reactive protein, Titin and Haptoglobin, were found to be sequentially altered with respect to increased parasite counts. Analysis of a longitudinal cohort of malaria patients indicated reversible alterations in serum levels of some proteins such as Haptoglobin, Apolipoprotein E, Apolipoprotein A1, Carbonic anhydrase 1, and Hemoglobin subunit alpha upon treatment; however, the levels of a few other proteins did not return to the baseline even during the convalescent phase of the infection. Here we present the first comprehensive serum proteomics analysis of vivax malaria patients with different levels of parasitemia and during the acute and convalescent phases of the infection.

What Is Known about the Immune Response Induced by Plasmodium vivax Malaria Vaccine Candidates?

Malaria caused by Plasmodium vivax continues being one of the most important infectious diseases around the world; P. vivax is the second most prevalent species and has the greatest geographic distribution. Developing an effective antimalarial vaccine is considered a relevant control strategy in the search for means of preventing the disease. Studying parasite-expressed proteins, which are essential in host cell invasion, has led to identifying the regions recognized by individuals who are naturally exposed to infection. Furthermore, immunogenicity studies have revealed that such regions can trigger a robust immune response that can inhibit sporozoite (hepatic stage) or merozoite (erythrocyte stage) invasion of a host cell and induce protection. This review provides a synthesis of the most important studies to date concerning the antigenicity and immunogenicity of both synthetic peptide and recombinant protein candidates for a vaccine against malaria produced by P. vivax.

The Aotus nancymaae erythrocyte proteome and its importance for biomedical research

The Aotus nancymaae species has been of great importance in researching the biology and pathogenesis of malaria, particularly for studying Plasmodium molecules for including them in effective vaccines against such microorganism. In spite of the forgoing, there has been no report to date describing the biology of parasite target cells in primates or their biomedical importance. This study was thus designed to analyse A. nancymaae erythrocyte protein composition using MS data collected during a previous study aimed at characterising the Plasmodium vivax proteome and published in the pertinent literature. Most peptides identified were similar to those belonging to 1189 Homo sapiens molecules; >95% of them had orthologues in New World primates. GO terms revealed a correlation between categories having the greatest amount of proteins and vital cell function. Integral membrane molecules were also identified which could be possible receptors facilitating interaction with Plasmodium species. The A. nancymaae erythrocyte proteome is described here for the first time, as a starting point for more in-depth/extensive studies. The data reported represents a source of invaluable information for laboratories interested in carrying out basic and applied biomedical investigation studies which involve using this primate.

SIGNIFICANCE

An understanding of the proteomics characteristics of A. nancymaae erythrocytes represents a fascinating area for research regarding the study of the pathogenesis of malaria since these are the main target for Plasmodium invasion. However, and even though Aotus is one of the non-human primate models considered most appropriate for biomedical research, knowledge of its proteome, particularly its erythrocytes, remains unknown. According to the above and bearing in mind the lack of information about the A. nancymaae species genome and transcriptome, this study involved a search for primate proteins for comparing their MS/MS spectra with the available information for Homo sapiens. The great similarity found between the primate's molecules and those for humans supported the use of the monkeys or their cells for continuing assays involved in studying malaria. Integral membrane receptors used by Plasmodium for invading cells were also found; this required timely characterisation for evaluating their therapeutic role. The list of erythrocyte protein composition reported here represents a useful source of basic knowledge for advancing biomedical investigation in this field.

Plasmodium vivax trophozoite-stage proteomes

Plasmodium vivax is the causative infectious agent of 80-300 million annual cases of malaria. Many aspects of this parasite's biology remain unknown. To further elucidate the interaction of P. vivax with its Saimiri boliviensis host, we obtained detailed proteomes of infected red blood cells, representing the trophozoite-enriched stage of development. Data from two of three biological replicate proteomes, emphasized here, were analyzed using five search engines, which enhanced identifications and resulted in the most comprehensive P. vivax proteomes to date, with 1375 P. vivax and 3209 S. boliviensis identified proteins. Ribosome subunit proteins were noted for both P. vivax and S. boliviensis, consistent with P. vivax's known reticulocyte host-cell specificity. A majority of the host and pathogen proteins identified belong to specific functional categories, and several parasite gene families, while 33% of the P. vivax proteins have no reported function. Hemoglobin was significantly oxidized in both proteomes, and additional protein oxidation and nitration was detected in one of the two proteomes. Detailed analyses of these post-translational modifications are presented. The proteins identified here significantly expand the known P. vivax proteome and complexity of available host protein functionality underlying the host-parasite interactive biology, and reveal unsuspected oxidative modifications that may impact protein function.

BIOLOGICAL SIGNIFICANCE

Plasmodium vivax malaria is a serious neglected disease, causing an estimated 80 to 300 million cases annually in 95 countries. Infection can result in significant morbidity and possible death. P. vivax, unlike the much better-studied Plasmodium falciparum species, cannot be grown in long-term culture, has a dormant form in the liver called the hypnozoite stage, has a reticulocyte host-cell preference in the blood, and creates caveolae vesicle complexes at the surface of the infected reticulocyte membranes. Studies of stage-specific P. vivax expressed proteomes have been limited in scope and focused mainly on pathogen proteins, thus limiting understanding of the biology of this pathogen and its host interactions. Here three P. vivax proteomes are reported from biological replicates based on purified trophozoite-infected reticulocytes from different Saimiri boliviensis infections (the main non-human primate experimental model for P. vivax biology and pathogenesis). An in-depth analysis of two of the proteomes using 2D LC/MS/MS and multiple search engines identified 1375 pathogen proteins and 3209 host proteins. Numerous functional categories of both host and pathogen proteins were identified, including several known P. vivax protein family members (e.g., PHIST, eTRAMP and VIR), and 33% of protein identifications were classified as hypothetical. Ribosome subunit proteins were noted for both P. vivax and S. boliviensis, consistent with this parasite species' known reticulocyte host-cell specificity. In two biological replicates analyzed for post-translational modifications, hemoglobin was extensively oxidized, and various other proteins were also oxidized or nitrated in one of the two replicates. The cause of such protein modification remains to be determined but could include oxidized heme and oxygen radicals released from the infected red blood cell's parasite-induced acidic digestive vacuoles. In any case, the data suggests the presence of distinct infection-specific conditions whereby both the pathogen and host infected red blood cell proteins may be subject to significant oxidative stress.

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

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