Plasmodium vivax Biology: Insights Provided by Genomics, Transcriptomics and Proteomics
- 1Laboratory of Tropical Diseases, Department of Genetics, Evolution, Microbiology and Immunology, University of Campinas - UNICAMP, Campinas, Brazil
- 2Laboratory of Regulation of Gene Expression, Instituto Carlos Chagas, Curitiba, Brazil
- 3Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
During the last decade, the vast omics field has revolutionized biological research, especially the genomics, transcriptomics and proteomics branches, as technological tools become available to the field researcher and allow difficult question-driven studies to be addressed. Parasitology has greatly benefited from next generation sequencing (NGS) projects, which have resulted in a broadened comprehension of basic parasite molecular biology, ecology and epidemiology. Malariology is one example where application of this technology has greatly contributed to a better understanding of Plasmodium spp. biology and host-parasite interactions. Among the several parasite species that cause human malaria, the neglected Plasmodium vivax presents great research challenges, as in vitro culturing is not yet feasible and functional assays are heavily limited. Therefore, there are gaps in our P. vivax biology knowledge that affect decisions for control policies aiming to eradicate vivax malaria in the near future. In this review, we provide a snapshot of key discoveries already achieved in P. vivax sequencing projects, focusing on developments, hurdles, and limitations currently faced by the research community, as well as perspectives on future vivax malaria research.
The last 10 years brought breakthroughs in biological research with the emergence of and improved accessibility to sequencing technology (Margulies et al., 2005). Nowadays, the increasingly available, user-friendly and less costly next generation sequencing (NGS) (Metzker, 2010) platforms are major technological contributions. Genomic exploration of organisms in greater detail through whole genome sequencing (WGS) is a reality. Additional valuable insights are coming from whole transcriptome (WTS) and proteome sequencing projects within the wider “omics” field, which start to shed light on multiple aspects of biology. Data integration within the “omics,” quantitative, functional and regulatory analysis - systems biology - is revolutionizing our understanding of the mechanisms of life. All these interconnected genomic, transcriptomic and proteomic platforms bring a new world of possibilities for hypothesis-driven research.
Genome sequencing projects revealed Plasmodium spp. genetic diversity and contributed to a better knowledge of parasite biology and host-parasite interaction (Hemingway et al., 2016). A significant number of fully sequenced and high quality Plasmodium spp. reference genomes are available (Carlton et al., 2008a; Dharia et al., 2010; Menard et al., 2010, 2013; Westenberger et al., 2010; Bright et al., 2012; Chan et al., 2012; Neafsey et al., 2012; Flannery et al., 2015; Winter et al., 2015; Hupalo et al., 2016; Pearson et al., 2016), and sequence curation efforts should continue. Today, the malaria research community ventures more and more into transcriptomics (Bozdech et al., 2008; Hoo et al., 2016; Zhu et al., 2016) and proteomics (Ray et al., 2016, 2017), areas to capitalize on for understanding Plasmodium spp. biology, especially the dynamics of RNA and protein expression and regulation through its complex multi-staged life cycle, in host and vector interaction contexts, and under different environmental selective pressures. Hence, the expertise provided by in-depth sequencing projects with clear data integration for understanding metabolic pathways in a systematic way is welcome by the parasite research community. Not only does it provide important pieces of information to understand different immune evasion and host invasion strategies, it is also a way to monitor and find new means to combat the rapidly increasing transmission of drug resistant parasites, and to identify molecular targets as starting points for effective vaccine development.
Plasmodium vivax malaria research has historically faced numerous technical adversities and has been largely neglected. This situation has led to a general lack of knowledge of P. vivax biology, and consequently impaired our capacity for making the best decisions on transmission control measures, and in the long run, for vivax malaria eradication. Currently vivax malaria is acknowledged as a disease that should no longer be neglected as it has been shown to result in considerable morbidity and mortality (Alexandre et al., 2010; Andrade et al., 2010; Lacerda et al., 2012a; Quispe et al., 2014; Rodriguez-Morales et al., 2015; Siqueira et al., 2015).
In this review, we show the main achievements on P. vivax biology accomplished based on the published genome, transcriptome and proteome sequencing projects. In particular, (1) the genome-wide comparative studies showing evolutionary relationships between parasites of the same genus, (2) the broad genetic diversity landscape within the P. vivax populations reported as to understand specific selection pressures (environmental and host/vector related) acting presently on parasite populations, (3) the more recent expression profile datasets and regulation mechanisms emerging from sensitive high-throughput WTS (RNA-seq) of P. vivax in different stages, and (4) the attempts to identify parasite metabolic pathways and antigens as possible diagnosis biomarkers through mass spectrometry (MS) based proteomics analysis of parasites and human host profiling from vivax malaria patient samples. Also, we present the main challenges encountered, remaining gaps and possible research avenues being explored and developed by the vivax malaria community.
Vivax Malaria: An Overview
Human malaria infections can be caused by five different Plasmodium species. P. falciparum is considered the deadliest parasite, causing the most severe clinical outcomes, whereas P. vivax is the most geographically spread within densely populated regions, thus accentuating the socio-economic burden caused by the disease (Gething et al., 2012; WHO, 2015). Recently, vivax malaria has re-emerged in regions formerly considered malaria free (Severini et al., 2004; Kim et al., 2009; Bitoh et al., 2011). Worldwide, about 2.85 billion people have been estimated to be at risk of infection by P. vivax (Price et al., 2007; Guerra et al., 2010; Battle et al., 2012; Gething et al., 2012). Following the decline of P. falciparum infections, P. vivax is now the dominant malaria species in several endemic regions (Coura et al., 2006; Gething et al., 2012; Hussain et al., 2013; WHO, 2015), where reports show a higher incidence, especially in young children (Marsh et al., 1995; Williams et al., 1997; Price et al., 2007; Genton et al., 2008; Tjitra et al., 2008). Several clinical complications that were normally associated with P. falciparum infections have been reported for vivax malaria (Kochar et al., 2005, 2009a; Hutchinson and Lindsay, 2006; Baird, 2007; Barcus et al., 2007; Alexandre et al., 2010; Rahimi et al., 2014). The most observed include anemia (Haldar and Mohandas, 2009; Quintero et al., 2011), haemolytic, coagulation disorders, jaundice (Sharma et al., 1993; Erhart et al., 2004; Lacerda et al., 2004; Saharan et al., 2009) and acute respiratory distress syndrome (Anstey et al., 2002, 2007; Suratt and Parsons, 2006; Tan et al., 2008; Lacerda et al., 2012b; Lanca et al., 2012), followed by nephropathology (Chung et al., 2008), porphyria (Kochar et al., 2009b), rhabdomyolysis (Siqueira et al., 2010), splenic rupture (de Lacerda et al., 2007; Gupta, 2010), and cerebral malaria (Lampah et al., 2011; Tanwar et al., 2011). Vivax malaria during pregnancy causing spontaneous abortions, premature and low weight new-borns (McGready et al., 2004; Poespoprodjo et al., 2008) is another major health concern, challenging the pre-established view of P. vivax as a “benign” parasite (Mendis et al., 2001; Baird, 2007; Anstey et al., 2009; Mueller et al., 2009; Gething et al., 2012; Naing et al., 2014).
According to the World Health Organization (WHO) guidelines (WHO, 2015), the first line P. vivax chemotherapy is chloroquine (CQ) plus primaquine (PQ), the only approved drug targeting the latent parasite form (Beutler et al., 2007). In high transmission areas presenting cases of drug resistance, an artemisinin based combination therapy (ACT) is recommended (WHO, 2010). The constant increase and spread of anti-malarial drug resistance remains of great concern (Baird, 2004; de Santana Filho et al., 2007; Suwanarusk et al., 2007; Poespoprodjo et al., 2008; Russell et al., 2008; Tjitra et al., 2008; Price et al., 2009, 2014).
P. vivax had an evolutionary path distinct from P. falciparum, being more closely related to P. cynomolgi, a sister taxon that infects Asian macaque monkeys (Duval et al., 2010; Liu et al., 2014; Luo et al., 2015). Probably as a consequence of such a unique evolutionary path, P. vivax shows unique biological features (Baird, 2007; Mueller et al., 2009; Gething et al., 2012) that distinguish it from P. falciparum (Figure 1):
(1) Preference for invading reticulocytes (RTs) (Field and Shute, 1956), increasing their deformability, size, fragility and permeability (Kitchen, 1938; Suwanarusk et al., 2004; Handayani et al., 2009; Desai, 2014), which would allow P. vivax to evade the host immune system (del Portillo et al., 2004; Suwanarusk et al., 2004; Handayani et al., 2009) and maintain its characteristic low biomass parasitemia (Mueller et al., 2009).
(2) Earlier production of sexual stages during the infection (Boyd and Kitchen, 1937), with a characteristic spherical shape seen in peripheral blood circulation even before the beginning of clinical symptoms, which might function as a reservoir to promote successful transmission to the mosquitoes (Boyd and Kitchen, 1937; Bousema and Drakeley, 2011).
(3) Formation of hypnozoites, which remain in the liver in a latent state (Krotoski et al., 1982; Baird et al., 2007), and for which the reactivation mechanisms are still unknown (Mueller et al., 2009).
Figure 1. Plasmodium vivax and Plasmodium falciparum life cycle comparison. (A) Anopheles blood meal infection: Malaria infection occurs once a female Anopheles spp. mosquito inoculates Plasmodium sporozoites into the skin of the human host during its feeding (Kiszewski et al., 2004). The sporozoites eventually reach the bloodstream. (B) Pre-erythrocytic stage infection: The sporozoites are transported through the vascular system to the liver, where they migrate across Kupffer or endothelial cells and enter hepatocytes. When a hepatocyte is found and within a period of 10–12 days, sporozoites from both Plasmodium species form a parasitophorous vacuole membrane and differentiate into schizonts. The schizogony process involves thousands of mitotic replications giving rise to a high number of merozoites packed into merosomes, which are then released into the bloodstream. P. vivax sporozoites can also differentiate into dormant long lasting liver forms called hypnozoites that upon activation start schizogony with consequent release of merozoites into the vasculature, and consequently causing clinical relapses (Krotoski, 1985). (C) Asexual erythrocytic stage: During this 48 h stage, P. vivax merozoites have a preference for reticulocyte invasion, whereas P. falciparum invade normocytes. Upon invasion, the parasite promotes several alterations in these blood cells, enlarging and deforming them (Suwanarusk et al., 2004) and promoting the formation of caveola-vesicle like complexes (CVC) and cytoplasmic cleft structures (Barnwell et al., 1990), resulting in an appropriate environment for asexual schizogony in the bloodstream. The species-specific set of surface proteins produced by the parasites influences the proportion of different parasite stages observed in the patient's peripheral blood and are suggested to be linked to the different pyrogenic capacity, biomass and parasitemia levels presented by vivax or falciparum malaria patients. RBC, Red Blood Cell. (D) Intra-erythrocytic gametocyte development: A proportion of asexual parasites are known to undergo gametocytogenesis into round-shaped sexual micro- and macrogametocytes, quite early (4 days) after P. vivax infection compared to P. falciparum (15 days post-infection) (Pelle et al., 2015), before detection of any clinical symptoms. (E) Mosquito Stage: Gametocytes circulating in the bloodstream can be taken up by the next Anopheles blood meal and engage in the sexual cycle with release of male and female mature gametes, fertilization (zygote) and formation of motile ookinetes that cross the mosquito midgut epithelium and further mature into oocysts. The oocyst, a new asexual sporogonic replicative entity, generates and releases thousands of sporozoites. These sporozoites migrate and invade the salivary glands of the mosquito vector and can then be transmitted to a new human host, completing the complex life cycle of the parasite (Mueller et al., 2009).
P. vivax research has been greatly restricted by lack of a reliable and reproducible in vitro system for long term culture (culturing efforts reviewed by Udomsangpetch et al., 2008; Noulin et al., 2013). This limitations are reflected both in the quantity and quality of functional assays performed in ex vivo short term cultures (Noulin et al., 2013). Furthermore, the alternative use of monkey-based studies with adapted strains of P. vivax (Beeson and Crabb, 2007; Panichakul et al., 2007) is of no easy access. Nevertheless, the available molecular tools (reviewed in Escalante et al., 2015) emerge as a valuable option to better estimate the prevalence and incidence of vivax malaria, its dynamics and differential contributions between host and vectors.
Plasmodium Vivax Genomics
Genome Sequencing Breakthroughs and Projects
In 2008, Carlton and colleagues achieved the complete generation, assembly and analysis of the first P. vivax WGS (Carlton, 2003; Feng et al., 2003; Carlton et al., 2008a), 6 years after the publication of P. falciparum reference genome (Gardner et al., 2002; Figure 2). The sequenced P. vivax Salvador-1 (Sal-1) strain came from a patient isolate from La Paz, close to El Salvador, and was passed to human volunteers, Aotus owl monkeys, and subsequently adapted to growth in Saimiri squirrel monkeys by mosquito and blood infection during the 1970s (Collins et al., 1972). Only later, it was possible to extract enough genomic DNA (gDNA) for whole genome shotgun sequencing (WGSS) and de novo assembly (Carlton et al., 2008a). Some technical hurdles were overcome, notably that the genome GC content of the parasite and laboratory host were similar and little information was available for the squirrel monkey gDNA sequence to exclude possible contamination (Carlton, 2003).
Figure 2. Time line of the major Plasmodium spp. sequencing projects. Genome, transcriptome, and proteome sequencing projects in the last four decades have contributed to our current parasite biology knowledge. P. falciparum (light red) and P. vivax (light blue) sequencing projects are shown, emphasizing the delay in time of around 5–6 years for the latter species. The first P. falciparum reference genome became available much earlier in 2002 (circled in light gray), then the corresponding P. vivax Sal-1 reference genome in 2008 (circled in light gray). WGS, whole genome sequencing; IDC, intraerythrocytic developmental cycle; MudPIT, multidimensional protein identification technology; MS, mass spectrometry; 1D LC-MS/MS, one-dimension liquid chromatography–mass spectrometry; RNA-seq, RNA sequencing; ChIP-seq, chromatin immunoprecipitation (ChIP) DNA sequencing.
The constraints initially faced by these studies to acquire sufficient gDNA of acceptable quality for WGS were alleviated by the use of NGS platforms. These allowed direct sequencing of the first patient ex vivo isolate in 2010 (Dharia et al., 2010) without the need for in vitro propagation (Figure 2) or use the of more traditional and laborious cloning and expression approaches (del Portillo et al., 1988). Since then, sequencing technology sensitivity has increased while required genetic material amounts have decreased, leading to lower costs. As a result, NGS is now a portable benchtop technology much closer to endemic areas. Parasite DNA can be enriched in low parasitemia blood samples and human contamination eliminated (Auburn et al., 2013). To achieve an effective preparation of P. vivax field isolates for WGS, a technique combining a CF11-based human lymphocyte filtration and the short-term ex vivo culture for schizont maturation was optimized and is currently applied with success in several labs set in endemic areas. An alternative approach to preventing host contamination and enriching samples in P. vivax DNA is the direct sequencing of parasite gDNA by hybrid selection (Melnikov et al., 2011). This sequencing methodology has allowed the first worldwide characterization of P. vivax isolates (Dharia et al., 2010; Menard et al., 2013).
Recently, several patient isolates from Peru, Madagascar, Malagasy, Cambodia, and a Belem monkey adapted strain were sequenced coupled with the resequencing of the Sal-1 strain (Carlton et al., 2008a) using whole genome capture (WGC), WGS and NGS techniques (Dharia et al., 2010; Bright et al., 2012; Chan et al., 2012; Neafsey et al., 2012; Winter et al., 2015; Figure 2). In 2016, major WGS analysis studies were published, including one sequencing 195 P. vivax genomes including 182 new high depth and quality sequences from clinical isolates from 11 different countries worldwide (Hupalo et al., 2016) and 13 already published sequences from both clinical isolates and monkey-adapted laboratory lines. In addition, 70 Cambodian P. vivax isolates, further compared to 80 P. falciparum isolates from the same region (Parobek et al., 2016) were sequenced, considerably augmenting the genomic data available to the malaria community (Luo et al., 2015; Hupalo et al., 2016; Parobek et al., 2016; Figure 2). Expectations are that these data will bring great advances in population genetics at a genome scale and that comparative genomic analysis (CGA) will function as a powerful method for parasite evolution hypothesis generation and testing (Feng et al., 2003).
Recently, a newly assembled and annotated reference genome named P. vivax P01 was published (Auburn et al., 2016a; Figure 2). The data was generated by using high-depth Illumina® sequencing from a Papua Indonesia patient isolate (Auburn et al., 2016a). For assembly, manual curation and comparison between data sets, the draft assemblies for the P. vivax isolates C01 (China) and T01 (Thailand) were used (Auburn et al., 2016a). The higher degree of assembly of these three genomes (~11 times) largely exceeds the Sal-1 reference genome (Carlton et al., 2008b; Auburn et al., 2016a). Not only was there an enhancement on scaffolds assembly with reduction from 2,500 unassembled scaffolds in Sal-1 to 226 in P01 reference genomes, as an improvement of the subtelomeric regions scaffold assembly. Therefore, annotation was also improved with the attribution of functions to 58% against the 38% in Sal-1 in a total of 4,465 core genes, defined as genes 1:1 orthologous between P. vivax P01 and P. falciparum 3D7 strains (Auburn et al., 2016a; Figure 3 and Table 1). The improvements in assembly and annotation quality in this sequencing project contributed to the generation of this very important new resource to study vivax malaria.
Figure 3. Plasmodium spp. morphology and genome architecture. Illustration showing the principal characteristics of Apicomplexan Plasmodium spp. parasite morphology, with its 3 parasite genomes, nuclear (ncDNA), mitochondrial (mtDNA) and apicoplast (apDNA).
Genome Architecture and Evolution
The biggest landmark achieved with the WGS of Plasmodium spp. was data acquisition to perform alignments and construct synteny maps with improved annotation, i.e., identifying conserved order of genetic loci along chromosomes (Ureta-Vidal et al., 2003). Through CGA one can now understand how the key forces that shape the evolutionary processes of organisms (mutation, selection, and genetic drift) came into play, by assuming that the genomes analyzed share a common ancestor (Frazer et al., 2003). For instance, CGA has allowed us to understand the integration of the apicoplast into apicomplexan species (Kohler et al., 1997) and to identify genes and their target signals transferred from the second endosymbiotic event into Plasmodium (Foth et al., 2003). Furthermore, CGA based on the P. falciparum and P. yoelii yoelii WGSs was the first such comparison between tropical pathogens of the same genus (Carlton et al., 2002; Feng et al., 2003), suggesting that chromosomal reshuffling might have been involved in speciation inside the Plasmodium genus.
Three different genomes, one nuclear, one mitochondrial and one apicoplast, characterize P. vivax (Figure 3). The ~26.8 megabases (Mb) nuclear genome is distributed among 14 chromosomes and has 42.3% GC average content, although presenting an isochore structure with high GC content in most internal chromosome regions, interspersed by high AT stretches mainly in subtelomeric regions (Carlton, 2003; Carlton et al., 2008a; Taylor et al., 2013; Table 1). The significance of Plasmodium spp. isochore structure remains unknown (Eyre-Walker and Hurst, 2001), however GC-rich genes evolve faster than AT-rich subtelomeric genes (Carlton, 2003; Carlton et al., 2008a). The codon usage in P. vivax is balanced (Nc~54.2), with less biased use of the 61 codons (Carlton et al., 2008b; Yadav and Swati, 2012; Cornejo et al., 2015). In particular, gene location in telomeric, centromeric or chromosome interstitial regions was recently shown to present different codon biased compositions (Cornejo et al., 2015). Genes showing the highest codon usage bias encode housekeeping proteins, the most highly expressed. Interestingly, genes in this category also include the Pf/Pv-fam family putatively responsible for antigenic variation and associated with early gametocyte formation (Cornejo et al., 2015). These patterns have led to speculation regarding the strong influence of the mutation rate on local nucleotide (nt) composition. This may indicate dramatic variance in gene expression potential and profiles, and that genes encoded in such regions could intervene in the specific pathophysiology of P. vivax malaria (McCutchan et al., 1984; Carlton et al., 2008b; Cornejo et al., 2015). As previously observed for other organisms, P. vivax DNA enriched in CpG motifs stimulates the Toll-like Receptor 9 (TLR9) increasing the host inflammatory response (Parroche et al., 2007), previously thought to be mainly caused by hemozoin toxicity alone. Given that P. vivax has a higher GC content, a greater TLR9 activation upon parasite DNA uptake by the host has been hypothesized (Anstey et al., 2009), potentially contributing to the increased pyrogenicity of this parasite (Karunaweera et al., 1992; Hemmer et al., 2006).
With reservations for sequence misassembly, it was estimated that the P. vivax parasite has around 5,400 genes (Carlton et al., 2008a; Table 1), of which a large majority was found to be orthologous between several Plasmodium spp. (CGA on human parasites P. vivax and P. falciparum, primate parasites P. cynomolgi and P. knowlesi and the rodent parasites P. yoelii. yoelii, P. chabaudi, and P. berghei; Carlton et al., 2008a). Orthologous genes identified had conserved genome positions with no genome breakage hot spots (Carlton et al., 2008a). Specifically, the more closely related primate species P. cynomolgi and P. knowlesi (Glazko and Nei, 2003) showed a higher degree of conservation throughout the 14 chromosomes (Cornejo et al., 2015) with exception for multigenic regions, where frequent microsyntenic breaks were identified (Carlton et al., 2008b). Calculation of substitution rates at synonymous (dS) vs. non-synonymous (dN) sites has allowed estimation of the relative importance of selection and genetic drift (dN/dS) throughout P. vivax evolution. Although a great difference between average dS and dN values was observed in ~3,300 high-confidence P. vivax/P. knowlesi orthologous gene pairs, these values correlate within and between syntenic regions of chromosomes of the two species. Carlton and colleagues have suggested heterogeneous mutation rates across the genome (Carlton, 2003). In Plasmodium, genes encoding membrane-anchored, transmembrane, cell adhesion, exported proteins or extracellular proteins with signal peptide motifs, most of them in direct contact with the host defense system, evolve relatively faster than the ones encoding housekeeping proteins (Carlton et al., 2008a). Selective constraint analysis between the highly conserved P. vivax/P. cynomolgi genomes reported a minority of genes (<2%) as positively selected and consequently evolving quicker than others with heavily constrained functions (Carlton et al., 2008a). Therefore, it appears that evolutionary rates have been greatly influenced by host-parasite interactions leading to a high degree of variation between different gene classes in Plasmodium (Carlton et al., 2008a; Joy et al., 2008).
It has been suggested that P. vivax had an African origin (Koepfli et al., 2015; Winter et al., 2015) and subsequently spread through the world (Culleton et al., 2011; Liu et al., 2014). The recently analyzed P. vivax samples were gathered in two main groups, Old and New World (Hupalo et al., 2016). Today the diversity between Old World P. vivax and P. falciparum samples is significantly higher. This is indicative of a large effective population size in P. vivax that shows signs of population specific natural selection, where P. vivax is constantly adapting to human host and mosquito vector regional differences, and more recently, to antimalarial drug pressures (Neafsey et al., 2012; Winter et al., 2015; Hupalo et al., 2016; Parobek et al., 2016). These results were also seen at the local level for a high number of P. vivax samples from the Asian-Pacific region (Brazeau et al., 2016; Parobek et al., 2016; Pearson et al., 2016) and, in a smaller scale, in South American samples (Winter et al., 2015). The non-uniform distribution of the single nucleotide polymorphisms (SNPs) in P. vivax, clearly enriched at subtelomeric regions, is being interpreted as a result of local high recombination rates. This leads to the hypothesis that selection not only acts on genes, but also has been shaping the P. vivax genomic architecture. Natural selection can thus act faster on these regions where the generation of antigenic variation occurs and also have a direct influence on genome architecture and consequently, impact P. vivax adaptation (Cornejo et al., 2015). Further, genome-wide analysis of P. vivax genetic diversity points at interesting evolutionary implications. For instance, the P. vivax genome is generally under very strong constraint with negative selection, and only very few genes are being positively selected (Cornejo et al., 2015; Parobek et al., 2016). Such results suggest a longer evolutionary history of primate infection (Liu et al., 2014; Luo et al., 2015), clearly closer to other monkey-infecting Plasmodium spp. (Carlton et al., 2013) than to P. falciparum, making it more adapted to survive and replicate in primate hosts.
These results have supported the use of monkey malaria Plasmodium spp. as models for the adaptation of P. vivax parasite populations (Chan et al., 2015), Moreover, amplification in monkey hosts of patient isolates from several endemic areas (Brazil, North Korea, India, and Mauritania) (Galinski and Barnwell, 2012) was crucial to collect enough high-quality gDNA for WGS (Neafsey et al., 2012; Carlton et al., 2013). However, analysis of these sequences should take into account that the parasites isolated from vivax malaria patients were passed through multiple infection cycles and were adapted to grow in Saimiri boliviensis monkeys (Carlton et al., 2008b). Host switching and adaptation to a new in vivo immune environment could affect mutation rates, and consequently the degree of P. vivax sequence variation, expression profiles and multiclonal complexity of parasite populations (Hester et al., 2013). Nevertheless, the study of Chan and colleagues proved that research using P. vivax monkey-adapted strains resulted in useful data, however data interpretation should be done with caution, given that these strains are not always genetically homogeneous (Chan et al., 2015). Overall, WGS reveals specific differences within human malaria Plasmodium spp. which should be taken into account in future interventions (Carlton et al., 2008a; Pain et al., 2008; Frech and Chen, 2011; Tachibana et al., 2012; Hester et al., 2013; Cornejo et al., 2015; Winter et al., 2015; Hupalo et al., 2016; Parobek et al., 2016; Pearson et al., 2016).
Plasmodium Vivax Transcriptomics
Whole Transcriptome Sequencing Projects
The first great effort into this parasite's transcriptome profile immediately followed the publication of P. vivax Sal-1 reference genome (Figure 2). By using a customized microarray platform of their own design, Bozdech et al. directly accessed P. vivax mRNA levels throughout 48 h intraerythrocytic developmental cycle (IDC) of three distinct isolates (Bozdech et al., 2008). Expression data analysis identified stage-specific differential expression of certain groups of genes. These genes were predicted to encode proteins with a role in parasitic development, virulence capacity and/or host-parasite interaction. In the following years, several microarray platform studies were published (Bozdech et al., 2008; Westenberger et al., 2010; Boopathi et al., 2013, 2014; Figure 2). Supporting the previous findings of Bozdech et al. on transcript levels at different lifecycle stages, Westenberger and colleagues described dramatic changes in messenger RNA (mRNA) co-expression of various genes predicted to be involved in developmental processes, suggesting that they could modulate parasitic progress through its different stages (Westenberger et al., 2010). In the 5′ region of co-expressed genes conserved motifs across Plasmodium spp. were identified as possible sites for regulatory protein binding with important roles in stage specific transcriptional regulation. Several multigene families and genes predicted to encode exported proteins displayed synchronized transcription in P. vivax sporozoites, but no difference (different gene set or mRNA levels) was observed between parasites showing capacity for hepatocyte invasion or hypnozoite development (Westenberger et al., 2010). Boopathi et al. reported a positive correlation of Natural Antisense Transcripts (NATs) with sense transcripts level, indicating that differing sense/antisense transcript ratios are involved in differential regulation of gene expression in diverse clinical conditions (Bozdech et al., 2008; Westenberger et al., 2010; Boopathi et al., 2013, 2014).
Very recently, RNA-seq was successfully applied to re-sequence two P. vivax isolates (Bozdech et al., 2008) throughout their IDC (Zhu et al., 2016; Figure 2). RNA-seq is more sensitive than gene expression microarrays as it detects expression nuances. In addition, RNA-seq includes all transcripts and not only the ones previously characterized present on microarray slides. The IDC transcriptome map produced by the high-resolution Illumina® HiSeq platform allows a better comprehension of the regulation behind the parasite gene expression profiles for specific biological functions. Strand-specific RNA-seq was just published for three isolates from Cambodian vivax malaria patient before anti-malarial treatment, and transcripts were assembled de novo revealing homogenous parasite gene expression profile regardless of the proportion of different stages (Kim et al., 2017). These results contrast with the gene expression pattern reported before (Bozdech et al., 2008; Zhu et al., 2016) and with data for sporozoites isolated from salivary glands of infected Colombian mosquitos (Kim et al., 2017). The authors recognize the fact that the patient isolates have multi-staged parasites, where one asexual stage might be transcriptionally more active than the others and, irrespective of its proportion, will deliver most transcripts, homogenizing the final expression patterns observed (Kim et al., 2017).
Differential Expression and Regulation
Plasmodium parasites do not have the typical eukaryotic transcriptional apparatus, presenting a poor repertoire of transcription associated proteins (TAP) somewhat similar among P. falciparum, P. vivax, and P. knowlesi, but with a significant range of putative regulatory sequences in the genome. This raises the idea of a dynamic, complex and non-standard gene expression regulation system (Carlton et al., 2008a; Boopathi et al., 2013, 2014; Adjalley et al., 2016). From the microarray and RNA-seq studies, some of the P. vivax transcriptome architecture characteristics have emerged. Throughout IDC of P. vivax (and other Plasmodium spp.), the overall timing of gene expression profiles does not differ significantly between strains and isolates. Hence, transcriptional differences or variations observed between species might reflect distinct evolution paths (Bozdech et al., 2008). Species-specific gene expression patterns in members of multigene families may mirror separate environmental challenges and multiple and varied host-parasite interactions faced by the parasites (Bozdech et al., 2008; Westenberger et al., 2010; Boopathi et al., 2013; Zhu et al., 2016). Considerable changes were found on orthologous genes (30% divergence) between P. falciparum (Bozdech et al., 2003; Le Roch et al., 2003) and P. vivax (Bozdech et al., 2008; Westenberger et al., 2010; Zhu et al., 2016). Expression of these genes peaks at different moments and possibly contributes to a specific timing reflecting species-specific biological functions. Among these genes are those related to host-parasite interactions, red blood cells (RBC) invasion, early intraerythrocytic development and antigenic variation. P. vivax antigen-presenting gene families (e.g., vir and trag genes) are activated mainly at the schizont-ring stage transition, indicative of different transcriptional control between the two species and, once more, underscoring the marked biological distinctions between the two most important human malaria parasites (Bozdech et al., 2003, 2008). In addition, contrary to what was seen for P. falciparum (Bozdech et al., 2003), P. vivax seems to have an enhanced capacity to adapt to its host and modify its virulence factors, as shown by its extensive intra-isolate gene expression diversity (Bozdech et al., 2008). RNA-seq data analysis revealed a transcriptional hotspot of vir genes on chromosome 2, suggesting them as important for immune evasion mechanisms. Structurally, P. vivax genes present unusually long 5′ untranslated regions (UTRs) and multiple transcription start sites (TSSs) (Zhu et al., 2016; Kim et al., 2017). Supporting previous observations, Hoo et al. published the first integrated transcriptome analysis for six Plasmodium spp. showing that conserved syntenic orthologs present early parasitic stage transcriptional divergence, where expression patterns follow the corresponding mammalian hosts (Hoo et al., 2016) and changes in important putative transcriptional regulators might explain the observed transcriptional diversity. Conversely, orthologs presenting similar transcriptional profiles across all Plasmodium spp. might be responsible for conserved and crucial functions (Hoo et al., 2016).
The importance of clusters of TSSs, the significance of long 5′ vs. shorter 3′ UTRs, alternative splicing events and the profusion of non-coding RNAs (ncRNAs), in a context of genome-wide different GC content, is starting to be characterized (Zhu et al., 2016; Kim et al., 2017). Different sequence motifs have already been associated with stage specific expression regulation (Westenberger et al., 2010), and recently Zhu and colleagues observed a substantial variation of the 5′ UTR owing to differential selection of TSS, although these phenomena do not seem to be involved in regulation of the asexual parasite cycle. P. falciparum transcriptome analyses reveal uniform usage of several clusters of TTSs, but in a differentiated way (Figure 4A; Watanabe et al., 2002; Lenhard et al., 2012; Adjalley et al., 2016). Such sites are spaced along the genome, indicative of dispersed patterns of transcription both in coding and intergenic regions (Watanabe et al., 2002; Lenhard et al., 2012; Adjalley et al., 2016). With a broad bidirectional promoter sequence possibly controlled by multiple regulatory elements, also characteristic of other eukaryotic species, the architecture and sequence properties of Plasmodium spp. chromosomes can play an important role in transcription initiation and regulation by means of chromatin remodeling (Watanabe et al., 2002; Lenhard et al., 2012; Adjalley et al., 2016). Moreover, as for P. falciparum, P. vivax shows a cycling transcription pattern throughout the IDC (Bozdech et al., 2003; Le Roch et al., 2003), with stage-specific regulation of transcription initiation events, well correlated with gene expression levels (Figure 4B; Siegel et al., 2014; Adjalley et al., 2016).
Figure 4. Plasmodium spp. gene expression and regulation dynamics. (A) Scheme illustrating the average observed nucleotide content and nucleosome occupancy surrounding the transcriptional start site blocks (black round tip arrow) controlling gene expression, with a preferential pyrimidine-purine (T/A) di-nucleotide usage (gray line) and a drop-in guanine-cytosine (GC) content (red line) around the starting site. (B) Generic representation showing an amplification of the suggested subtelomeric (blue strip on the chromosome drawing) transcriptional start site clusters (round tip arrows) dynamics for differential gene expression during parasitic intraerythrocytic developmental cycle (~48 h). The graph displays arbitrary subtelomeric gene expression profiles under different TSS clusters control, during ring (green, peak at 6 h), trophozoite (red, peak at 24 h) and schizont (blue, peak at 36 h) stages. IDC, intraerythrocytic developmental cycle; TSS, transcriptional start site.
Alternative splicing as a mechanism of producing transcriptome variability was reported for P. falciparum with a role in parasite sexual development (Iriko et al., 2009; Otto et al., 2010; Sorber et al., 2011), contrasting with the scarcity reported until now for P. vivax (Zhu et al., 2016). However, the few splicing events identified are linked to a late schizont stage, suggesting its significance at the gene function level (Zhu et al., 2016). Recently, Kim and colleagues described a high percentage of P. vivax genes encoding multiple, often uncharacterized, protein-coding sequences (Kim et al., 2017). In addition to long 5′ UTRs, many of the isoforms significantly differ in their 5′- or 3′-UTRs, possibly as a result of differential exon splicing. P. vivax may use such protein isoforms in transcription and/or translation regulation mechanisms. Importantly, examples of such cases support the idea that even for well-described genes conferring drug resistance, novel isoforms can be revealed through transcriptomic data analysis, which could aid in decoding molecular mechanisms responsible for antimalarial drug resistance (Hoo et al., 2016; Kim et al., 2017).
All P. vivax transcriptomic data allowing for ncRNAs identification reveal the presence of these RNAs at high levels, notably in the antisense direction (Zhu et al., 2016; Kim et al., 2017). In 2013, the NATs identified in P. vivax directly isolated from infected patients with different malaria complications by genome-wide transcriptome analysis using customized high density tiling microarray sequencing technology (Boopathi et al., 2014), suggested a possible role in differential regulation of P. vivax gene expression. The involvement of distinct types of transcriptional machinery and different mechanisms, such as transcriptional run-through, bidirectional and antisense-specific promoters, produces a varied landscape of differential regulation observed for P. vivax from malaria patients in diverse clinical conditions (Boopathi et al., 2013). Surprisingly, the extensively investigated micro RNAs (miRNAs) class of ncRNAs appears to be absent in some Apicomplexan species including Plasmodium. However, an increasing number of studies have described a role of host miRNAs in host-parasite interactions (reviewed by Judice et al., 2016), as host miRNA expression can change after parasite infection and, within this context, suggesting a crucial role for the host immune system signaling.
The incorporation of RNA sequencing data will be extremely important to further work on annotated genes and continue the curation and further characterization of putative genes sets that are being identified in new P. vivax isolates under WGS projects. On the other hand, transcript data will indicate the significance of SNPs and increase of copy number variation (CNV) during the hypothesized recent population expansion in P. vivax (Mu et al., 2005; Taylor et al., 2013; Cornejo et al., 2015; Parobek et al., 2016). Furthermore, genome-wide transcriptome analysis already contributes to a better understanding of the dynamics of P. vivax gene expression by annotating new coding and non-coding transcripts, determining absolute levels of some transcripts, and characterizing the biological importance of UTRs, TSSs and alternative splicing sites (Zhu et al., 2016). Most importantly, it gives an understanding of the role of the transcriptional regulatory machinery and fine-tuned mechanisms that promote the large and diverse population upkeep of P. vivax, ability to stand up to selective sweeps, including those pressures caused by malaria control measures (Parobek et al., 2016). Even though only few studies into P. vivax transcriptional regulation have been published (Bozdech et al., 2008; Hoo et al., 2016; Zhu et al., 2016), some results suggest that specific biological processes that have enabled the parasite to avoid the traditional malaria control measures might be under tight transcriptional control. For example, early gametogenesis (Bousema and Drakeley, 2011) along with the asexual cycle seems to be achieved by transcriptional repression of an AP2 transcription factor (Yuda et al., 2015; Hoo et al., 2016). Another example is the biological process underlying hypnozoite latent stage formation and activation, suggested to be under epigenetic regulation. Studies using simian hepatocyte cultures report an acceleration of hypnozoite activation when histone methylation is inhibited. Thus, methylation of histones promotes suppression of transcription, and ultimately, maintain hypnozoite dormancy (Barnwell and Galinski, 2014; Dembele et al., 2014). Other example of P. vivax transcriptional regulation are related to mechanisms of drug resistance emergence, has shown for differential expression of CRT protein, rather than pvcrt gene coding sequence alteration, mediating CQ resistance (Fernandez-Becerra et al., 2009; Melo et al., 2014; Pava et al., 2015).
Plasmodium Vivax Proteomics and Metabolomics
Parasite and Patient Serum Proteome Analysis
In spite of all advances mentioned above, the functions of half of the predicted proteins in P. vivax remain unknown. Unsurprisingly, P. vivax and P. falciparum share a similar metabolic potential with key housekeeping pathways and functions and a range of putative membrane transporters, including the important apicoplast metabolism, essential for the parasite (Carlton et al., 2008a). The apicoplast proteome plays an important role in several major parasite metabolic processes, including some described for P. falciparum (Ralph et al., 2004) including isopentenyl diphosphate and iron sulfur cluster assembly, type II fatty acid synthesis and haem synthesis pathways, the latter subdivided between the apicoplast and mitochondria. The localization of several other central pathways is being clarified. For example, it was confirmed that the glyoxalase pathway occurs in the apicoplast, although in P. vivax, thiamine pyrophosphate biosynthesis takes place in the cytosol (Baird, 2010).
Although P. falciparum studies gave us some indications on the principal features of the P. vivax proteome, parasite-specific biology probably involves expression of different gene families. Thus, alternative invasion and host escape pathways still remain largely unknown. Stage-specific analysis of P. vivax IDC from natural isolates is extremely challenging due to low parasitemias, characteristically asynchronous parasite populations, and frequent polyclonal infections. In order to overcome the low parasitemia data challenge, some studies have pooled samples from different patients. However, such procedures can potentially increase the assortment of proteins detected, render their identification and association with disease pathogenesis more difficult, and affect the reproducibility of vivax proteomic data and follow up validation studies. Proteome studies are expected to contribute greatly by disclosing proteins expressed in clinical isolates and reveal new unique parasite pathways involved in malaria pathophysiology, leading to identification of new targets for drug and vaccine development.
Efforts are being made to identify potential vivax malaria biomarkers to understand host responses. The earliest proteome study, published in 2009 (Acharya et al., 2009), identified 16 proteins from a single patient infected with P. vivax blood-stage parasites (Figure 2). In 2011, the same authors were the first to examine the P. vivax proteome directly from a pool of clinical isolates, identifying a total of 153 P. vivax proteins, most of them with unknown homology to expressed genes (Acharya et al., 2011). That same year Roobsoong et al. (2011) published the P. vivax schizont proteome from multi-patient schizont cultures and enriched samples, where they reported 316 proteins of which 36 were unique to this parasite, and half of these were unique hypothetical proteins. Such proteins could be P. vivax-specific targets for diagnosis and treatment, as the four novel P. vivax antigens described (Roobsoong et al., 2011). More recently, 238 trophozoite proteins were identified in P. vivax strain VCG-1 together with 485 from the Aotus host (Moreno-Perez et al., 2014). The in-depth analysis of two P. vivax Sal-1 proteomes (S. boliviensis monkey iRBCs at the trophozoite stage) allowed the identification of 1375 parasite and 3209 host proteins and their post-translational modifications, mostly N-terminal acetylation (Anderson et al., 2015).
High-throughput screening has been used to study P. vivax immune proteomes, leading to identification of 44 antigens from a total of 152 protein putative candidates, and characterization and confirmation of rhoptry-associated membrane antigens (PvRAMA) as relevant serological markers of recent exposure to infections (Lu et al., 2014). Furthermore, differences between P. vivax and P. falciparum antigenic genes were revealed. Most of these highly immunoreactive proteins are hypothetical, but there have been suggestions of their importance for P. vivax invasion and evasion mechanisms (Figure 5). For instance, reticulocyte binding protein 2 and RAMA maybe involved in selectivity for invasion of RTs and/or exported protein 1 and 2, histidine-rich knob protein homolog and aspartic protease PM5 for host immune evasion (Lu et al., 2014).
Figure 5. P. vivax variant proteins. Slices of the P. vivax infected erythrocyte (Pv-iE) are shown with known variant proteins (VIR, PvTRAg, DBP, RBP, PvHIST, MSP) involved in P. vivax merozoite invasion (A), transport structures mediating the export of parasite ligands on the surface that promote changes in the host erythrocyte (B) and proteins associated with host evasion mechanisms, e.g., the rosetting phenotype (C). RBC, red blood cell; iE, infected erythrocyte; VIR, variable interspersed repeat multigene family proteins; PvTRAg, P. vivax tryptophan-rich antigen family proteins; PvHIST, P. vivax Plasmodia Helical Interspersed SubTelomeric multigene family proteins; DBP, Duffy binding protein; RBP, reticulocyte binding proteins; MSP, Merozoite surface proteins; Pexel, Plasmodium export element protein motif.
Humoral immunity of vivax malaria patients has also been interrogated. Chen and co-authors cloned and expressed 89 proteins that were screened against sera from vivax malaria patients using protein arrays to show 18 highly immunogenic of which 7 had been already characterized as vaccine candidates and the remaining were uncharacterized (Chen et al., 2010). Later, Ray and colleagues (Ray et al., 2012a,b) identified proteins such as apolipoprotein A1 and E, serum amyloid A and P, haptoglobin, ceruloplasmin, and hemopexin as differentially expressed in uncomplicated malaria patients by serum proteome analysis. Collectively, these results indicate that several physiologically important host pathways have been modulated during parasitic infection including vitamin D metabolism, hemostasis and the coagulation cascade, acute phase response and interleukin signaling, and the complement pathway (Ray et al., 2012a).
Host-Parasite Metabolic Pathways
To uncover pathway modulation occurring during uncomplicated (Ray et al., 2012a,b) vs. severe vivax malaria (Ray et al., 2016), researchers are going further into clinicopathological analysis and proteomic profiling. Recently, human serum proteome comparative studies were published (Ray et al., 2016) showing different levels of parasitemia for severe and non-severe infections, as well as during the acute and convalescent phases of the infection (Ray et al., 2017). The involvement of the acute phase response was confirmed, including acute phase reactants or proteins, cytokine signaling, oxidative stress and anti-oxidative pathways, muscle contraction and cytoskeletal regulation, lipid metabolism and transport, complement cascades, and several coagulation and hemostasis associated proteins in malaria pathophysiology (Ray et al., 2016, 2017). However, alterations in the blood coagulation cascade, as reported for severe falciparum malaria, were not identified for severe vivax malaria (Ray et al., 2016). Moreover, these proteomic studies allowed the identification of prospective new host markers, such as the differentially expressed proteins superoxide dismutase, ceruloplasmin, vitronectin, titin, and nebulin. Furthermore, they confirmed the already investigated apolipoprotein A1 and E, serum amyloid A and haptoglobin markers. Those are important for future improved diagnosis, discrimination of other infections, as well as disease progression and shift toward severe clinical manifestations, response to new therapies and outcome prediction (Ray et al., 2016). Analysis of a recent longitudinal cohort of vivax malaria patients revealed that some serum levels of proteins such as haptoglobin, apolipoprotein E, apolipoprotein A1, carbonic anhydrase 1, and hemoglobin subunit alpha, revert to baseline levels under treatment, but others did not show the same behavior in the convalescent phase of the infection (Ray et al., 2017), which might reflect the effect on the host of different parasitemias present. Serum proteome results also prompt the malaria community for a clearer definition of severity parameters for P. vivax malaria and show a marked difference between the pathogenesis of infection caused by different Plasmodium parasites (Ray et al., 2016).
The panorama emerging from all malaria proteome studies is that a great fraction of proteins identified until now are uncharacterized and of unknown functions, alongside with several housekeeping proteins. However, other proteins of major metabolic pathways for parasite survival show us that we are in a good position to better understand immune evasion and host cell invasion mechanisms within the pathophysiology of vivax malaria. These pathways comprise metabolism (glycolysis, hemoglobin digestion, nucleic acid synthesis) and cellular invasion (binding proteins, protein synthesis, modification, and degradation) (Figure 5A). Intracellular transport, translocation and presentation of variable antigen proteins (Figure 5B), which promote erythrocyte modification (Figure 5C), have been carefully analyzed, in particular proteins directly related to drug resistance. All chaperonin complexes, some involved in heat shock, oxidative stress and other counteractive responses, have been identified (Acharya et al., 2009, 2011; Chen et al., 2010; Roobsoong et al., 2011; Ray et al., 2012a,b, 2016; Lu et al., 2014; Moreno-Perez et al., 2014; Anderson et al., 2015).
Although there was evidence of transcription for genes encoding almost all P. vivax proteins identified, no correlation was seen between levels of mRNA transcripts and the corresponding protein expression (Bozdech et al., 2008; Westenberger et al., 2010; Bautista et al., 2014; Moreno-Perez et al., 2014; Anderson et al., 2015; Parobek et al., 2016). Integrating genomic and transcriptomic information is paving the way for systems biology approaches. For instance, chokepoint analysis has the potential to uncover enzymes predicted to have only one substrate or product, thus pushing the way to new drug target discoveries. Lists of compounds that could target each chokepoint are being generated using bioinformatics tools from “omics” data analysis.
Analysis of the hypnozoites proteome faces great experimental hurdles. A list of candidate genes includes all homologs of other dormancy genes previously discovered in other organisms (Carlton et al., 2008a). An in vitro system for culturing liver stages of this parasite (Mazier et al., 1984; Hollingdale et al., 1985; Sattabongkot et al., 2006; Cui et al., 2009) and in vivo experimental models are currently under development and will be of importance to approach hypnozoite proteomics. Taken together, blood proteome studies shed light on vivax malaria pathogenesis and pave the way for future integrated multi-omics investigations. However, such an approach has to be interpreted in light of host-parasite interaction, where some reported alterations might be host or parasite related and/or the result of cumulative effects between both parties. Ex vivo specific functional assays may provide further insights.
Plasmodium vivax Genetic and Transcriptomic Diversity, Host-Parasite Dynamics, and Disease Control
The first P. vivax genetic diversity reports were based on a few loci. Researchers relied on the amplification of a small set of polymorphic antigens to assess the genetic diversity of P. vivax. In general, the circumsporozoite protein (CSP), the merozoite surface proteins (MSP), the apical membrane antigen (AMA-1) (Mueller et al., 2002; Cui et al., 2003; Kim et al., 2006; Aresha et al., 2008; Moon et al., 2009; Zakeri et al., 2010), Duffy binding protein (DBP) (Cole-Tobian and King, 2003) and mitochondrial DNA were analyzed (Jongwutiwes et al., 2005; Cornejo and Escalante, 2006) and brought forward the polyclonal nature of P. vivax natural infections. Feng et al. published the first comparative population genomics study between P. falciparum and P. vivax, revealing shared macroscale conserved genetic patterns, but with P. vivax populations showing a highly polymorphic genome for which the degree of diversity depends on the gene context. Thus, coding regions are more conserved than intergenic ones, and central chromosomal regions more than (sub)telomeric, indicative of a purifying selection. This implies an ancient evolutionary history shaped by distinct selection pressures (Feng et al., 2003; Chang et al., 2012). Genes encoding transcription factors (type AP2) and ABC transporters associated with multidrug resistance showed higher dN/dS rate (Dharia et al., 2010; Parobek et al., 2016). The genetic variability seen for P. vivax transcription factors (Hoo et al., 2016) could indicate that malaria parasites evolve by changing and/or expanding their transcription factor DNA binding domains in response to drug administration and host immune challenges (Dharia et al., 2010; Parobek et al., 2016).
Several studies have highlighted the great degree of diversity, which is shared between P. vivax samples of different geographic locations at the population level, strongly suggesting that recent evolutionary selective pressures are currently acting upon some genes at a few loci, mainly those related with drug resistance (Mu et al., 2005; Imwong et al., 2007a; Gunawardena et al., 2010; Neafsey et al., 2012; Lin et al., 2013). Compared to P. falciparum natural isolates from similar geographical regions (Cheeseman et al., 2016), results exposed a greater genetic diversity in P. vivax natural isolates, both for SNP and microsatellites (Karunaweera et al., 2008; Neafsey et al., 2012; Orjuela-Sanchez et al., 2013; Taylor et al., 2013; Barry et al., 2015; Winter et al., 2015; Brazeau et al., 2016; Friedrich et al., 2016; Hupalo et al., 2016; Pearson et al., 2016). While microsatellite marker analysis has revealed the genetic diversity of P. vivax infections, they may provide over- or underestimates. This can be attributed to their unstable nature in comparison to SNPs, especially because of the smaller number and/or less informative use of these markers. Therefore, P. vivax high density tiling microarrays (Dharia et al., 2010), WGC sequencing (Bright et al., 2012) and WGS (Bright et al., 2012; Chan et al., 2012; Neafsey et al., 2012; Hester et al., 2013; Winter et al., 2015; Hupalo et al., 2016; Pearson et al., 2016) are being used to find informative SNPs. Whole genome analysis of genetic diversity between five P. vivax natural isolates (Neafsey et al., 2012) and the Sal-1 strain (Carlton et al., 2008b) has revealed shared gene diversity, but also significant numbers of unique SNPs for each isolate (Hester et al., 2013; Cornejo et al., 2015; Hupalo et al., 2016). Genomic and transcriptomic sequencing of several patient P. vivax isolates resulted in drug resistance phenotype identification (Bright et al., 2012; Mideo et al., 2013; Lin et al., 2015), emergence detection, transmission rates follow-up, as well as the identification (also through proteomics) of rapidly evolving antigens (Boyd and Kitchen, 1937; Cornejo et al., 2015; Winter et al., 2015; Hupalo et al., 2016; Parobek et al., 2016; Pearson et al., 2016). Genomics data suggest that evolutionary pressures acted upon and currently shape several P. vivax loci associated with invasion of RBC (Figures 5A,B), host and vector immune evasion mechanisms (Figures 5B,C) and antifolate drug resistance, among others (Imwong et al., 2003; Korsinczky et al., 2004; Dharia et al., 2010; Flannery et al., 2015; Winter et al., 2015; Hupalo et al., 2016; Pearson et al., 2016), which could support target discovery (Carlton et al., 2013; Cornejo et al., 2015; Hupalo et al., 2016; Pearson et al., 2016). For instance, WGS of two patient isolates from Madagascar allowed the identification of recently arisen SNPs present in the DBP region II (Chan et al., 2012). The detection of pvdbp gene duplication (Menard et al., 2013; Hupalo et al., 2016; Pearson et al., 2016) suggests that this gene is under strong positive selection through influence of the human host, and is rapidly evolving in response to the Duffy- allele present in most of sub-Sharan African population, where P. vivax is believed to have originated (Miller et al., 1976; Culleton et al., 2011; Liu et al., 2014). Surprisingly, similar duplications at high frequency have been recently seen in other geographic locations where almost all individuals are Duffy+ (Howes et al., 2011; Pearson et al., 2016). More recently, Auburn and colleagues reported amplification breakpoints for the multidrug resistance 1 gene (pvmdr1) (Auburn et al., 2016b). Thus, P. vivax is considered within an evolutionary neutral model, possibly the outcome of repeated gene duplications (population size expansion) of the ancestral lineage during the course of evolution, or that several sub-populations are present (Parobek et al., 2016). Multi-species WTS data suggests that variability in non-coding sequences of genes but other features such as transcription factors, upstream bidirectional promoter regions, post-transcriptional control including epigenetic regulation, chromatin remodeling events or ncRNAs may have a greater impact on transcription modulation across the Plasmodium spp. (Hoo et al., 2016).
With few exceptions verified in small-scale studies (Winter et al., 2015), P. vivax natural infections are almost always polyclonal (Parobek et al., 2016). However, some haplotype reconstruction studies have reported that 2–4 strains account for the majority of P. vivax gDNA, set against a situation of several different strains equally abundant in one patient (Chan et al., 2012; Neafsey et al., 2012; Winter et al., 2015; Friedrich et al., 2016; Pearson et al., 2016). For instance, deep genome sequencing of more than 200 clinical Asian-Pacific samples confirmed the complex genetic structure of P. vivax populations showing variations both in the number of dominant clones in individual infection, but also in their degree of relatedness and inbreeding (Pearson et al., 2016). Still, this remarkable polyclonality can have several origins: parasite infection could originate from a single meiosis event, from multiple infections with unrelated parasites as a consequence of several mosquito bites, or from a combination of a relapse and a new infection event, all of them having implications on parasite population diversity (Lin et al., 2015). Microsatellite marker identification and WGS have allowed the characterization of relapse infections as the result of heterologous hypnozoite activation (Imwong et al., 2007b; Carlton et al., 2008a; Restrepo et al., 2011; de Araujo et al., 2012; Bright et al., 2014; Lin et al., 2015). Furthermore, hypnozoites present hidden in the hosts function as a reservoir and a dynamic and continuous source and flow of genetic diversity in the present P. vivax population (White, 2011; Neafsey et al., 2012; Menard et al., 2013; White et al., 2016).
Even though WGS allows determination of informative SNPs, it lacks sensitivity to detect low frequency CNVs present within a polyclonal infection context. The tailored and more affordable P. vivax SNPs barcode project (Baniecki et al., 2015), transcriptome profiling microarrays (Boopathi et al., 2016) and other breakpoint-specific PCR approaches (Auburn and Barry, 2017) can be used as important tools to not only discriminate P. vivax infections and their origin, but also to deliver information of its population dynamics of transmission (Friedrich et al., 2016). Additionally, they are helping to characterize and discern between recrudescence, relapses (Lin et al., 2015), de novo infections (Friedrich et al., 2016) and monitoring drug resistance emergence (Auburn et al., 2016b; Brazeau et al., 2016). As for P. falciparum (Cheeseman et al., 2016), this greater genetic/expression diversity may provide P. vivax with a wide range of alternative paths for successful host erythrocyte invasion (Figures 5A,B) and immune system evasion (Figure 5C), thus enhancing its capacity of adaptation and survival (Cornejo et al., 2015). In addition, the production of an effective vivax malaria vaccine is hampered by the existence of such huge diversity of antigens concentrated in invasion and immune evasion related genes (Neafsey et al., 2012; Winter et al., 2015; Hupalo et al., 2016; Pearson et al., 2016) in combination with the absence of allelic exclusion and clonal interference mechanisms (Fernandez-Becerra et al., 2005).
Plasmodium vivax Omics Outlook: Achievements, Gaps, and Future Prospects
During the last 10 years, several P. vivax sequencing projects within the three big omics branches, genomics, transcriptomics and proteomics, have greatly contributed with insights into the parasite's biology (Table 2). The lack of a P. vivax in vitro long-term culturing system is the root of all limitations for its study. Nevertheless, we are living in times of rapid development, expansion and application of high-throughput and highly sensitive technology tools that were absurdly expensive only a few years ago. Exciting progress is being made in single-cell (Nair et al., 2014) and single-molecule sequencing, bringing prospects of complex problem-solving capacity that were previously unthinkable. The development of methodologies such as CRISPR/Cas9 system for genome editing opens new possibilities to understand and further explore gene expression and regulation of malaria parasites. The application of such tools to in vitro cell or animal experimental models may help us explore further some aspects of P. vivax biology and the pathophysiology of disease (Table 3). Specifically, the development, establishment and validation of in vitro cell culture and animal experimental models (as the published Mazier et al., 1984; Hollingdale et al., 1985; Sattabongkot et al., 2006; Cui et al., 2009; Mikolajczak et al., 2015; Mehlotra et al., 2017 reports) to support P. vivax growth would overcome the current number of parasites available for projects in all omics branches. Clearly, this would positively impact on our understanding of several molecular mechanisms taking place in the parasite during cell invasion, development, stage transition, host immune system evasion and environmental stress responses (Table 3).
The extensive and detailed omics datasets that are continuously made available by several research groups, under which technical and experimental progress has and is being achieved brings forward an old problem: how to analyze and integrate large and complex information within and between datasets in a standard and easy-to-share platform, accessible by all. Several initiatives, such as the comparative databases PlasmoDB (Bahl et al., 2003), support malaria research communities by offering computational biology training courses. Such platforms should be continuously developed with the most recent powerful bioinformatics tools and information combined with suitable mathematical modeling and epidemiologic analysis. Mainstream bioinformatics institutes around the world are providing expert support to the Plasmodium research community, critically needed given the extraordinary challenges posed by the complex nature of the P. vivax genome and population genetics. Worldwide, major efforts must be taken by the entire community in order to establish some criteria for dealing with data collection, analysis and annotation, so that it can be easily accessed and fully explored in order to design adequate interventions to ever changing conditions, focusing on the control, elimination and prevention of malaria.
CB and LA contributed to literature review and writing of the manuscript. CB and AK were responsible for the original figure concept and design. LA, PS, and FC helped with critical revision of the manuscript for intellectual content. All authors read and approved the final manuscript.
This work was supported by Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP 2012/16525-2), the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), and the Swedish Research Council (2014-4234). CB was supported by a FAPESP fellowship (FAPESP 2013/20509-5). FC is a CNPq research fellow and AK was supported by a CNPq fellowship (150030/2017-7).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
ACT, Artemisinin based Combination Therapy; AMA, Apical Membrane Antigen; apDNA, apicoplast DNA; bp, base pairs; CGA, Comparative Genome Analysis; ChIP-seq, chromatin immunoprecipitation (ChIP) DNA sequencing; CQ, Chloroquine; CNV, Copy Number Variation; CSP, Circumsporozoite Protein; CVC, caveola-vesicle like complexes; DBP, Duffy Binding Protein; dN, non-synonymous substitution rate; dS, synonymous substitution rate; ETRAMP, Early Transcribed Membrane Protein; gDNA, genomic DNA; IDC, Intraerythrocytic Development Cycle; iRBC, infected Red Blood Cell; Kb, kilobases; Mb, megabases; miRNA, micro RNA; mRNA, messenger RNA; MS, Mass Spectrometry; MSP, Merozoite Surface Protein; MudPIT, multidimensional protein identification technology; mtDNA, mitochondrial DNA; NAT, Natural Antisense Transcript; ncDNA, nuclear DNA; ncRNA, non-coding RNA; NGS, Next Generation Sequensing; nt, nucleotide; Pexel, Plasmodium export element protein motif; PQ, Primaquine; PST-A, lysophospholipase; PvHIST, P. vivax Plasmodia Helical Interspersed SubTelomeric multigene family proteins (also named as Pf-fam-b) and RAD protein (Pv-fam-e); PvTRAg, P. vivax tryptophan-rich antigen family proteins, also named as Pv-fam-a and tryptophan-rich antigens; PvSTP-1, P. vivax STP1 protein; RBC, Red Blood Cell; RBP, Reticulocyte Binding Protein; RAMA, Rhoptry-Associated Membrane Antigen; RNA-seq, RNA sequencing; RT, Reticulocytes; SNP, Single Nucleotide Polymorphism; TAP, Transcriptional Associated Protein; TLR, Toll-like Receptor; TSS, Transcription Starting Site; UTR, Untranslated Regions; VIR, Variable Interspersed Repeat; WGC, Whole Genome Capture; WGS, Whole Genome Sequencing; WGSS, Whole Genome Shotgun Sequencing; WHO, World Health Organization; WTS, Whole Transcriptome Sequencing; 1D LC-MS/MS, one-dimension liquid chromatography–mass spectrometry.
Acharya, P., Pallavi, R., Chandran, S., Chakravarti, H., Middha, S., Acharya, J., et al. (2009). A glimpse into the clinical proteome of human malaria parasites Plasmodium falciparum and Plasmodium vivax. Proteomics Clin. Appl. 3, 1314–1325. doi: 10.1002/prca.200900090
Acharya, P., Pallavi, R., Chandran, S., Dandavate, V., Sayeed, S. K., Rochani, A., et al. (2011). Clinical proteomics of the neglected human malarial parasite Plasmodium vivax. PLoS ONE 6:e26623. doi: 10.1371/journal.pone.0026623
Adjalley, S. H., Chabbert, C. D., Klaus, B., Pelechano, V., and Steinmetz, L. M. (2016). Landscape and dynamics of transcription initiation in the malaria parasite Plasmodium falciparum. Cell Rep. 14, 2463–2475. doi: 10.1016/j.celrep.2016.02.025
Alexandre, M. A., Ferreira, C. O., Siqueira, A. M., Magalhaes, B. L., Mourao, M. P., Lacerda, M. V., et al. (2010). Severe Plasmodium vivax malaria, Brazilian Amazon. Emerging Infect. Dis. 16, 1611–1614. doi: 10.3201/eid1610.100685
Anderson, D. C., Lapp, S. A., Akinyi, S., Meyer, E. V., Barnwell, J. W., Korir-Morrison, C., et al. (2015). Plasmodium vivax trophozoite-stage proteomes. J. Proteomics 115, 157–176. doi: 10.1016/j.jprot.2014.12.010
Andrade, B. B., Reis-Filho, A., Souza-Neto, S. M., Clarencio, J., Camargo, L. M., Barral, A., et al. (2010). Severe Plasmodium vivax malaria exhibits marked inflammatory imbalance. Malar. J. 9:13. doi: 10.1186/1475-2875-9-13
Anstey, N. M., Handojo, T., Pain, M. C., Kenangalem, E., Tjitra, E., Price, R. N., et al. (2007). Lung injury in vivax malaria: pathophysiological evidence for pulmonary vascular sequestration and posttreatment alveolar-capillary inflammation. J. Infect. Dis. 195, 589–596. doi: 10.1086/510756
Anstey, N. M., Jacups, S. P., Cain, T., Pearson, T., Ziesing, P. J., Fisher, D. A., et al. (2002). Pulmonary manifestations of uncomplicated falciparum and vivax malaria: cough, small airways obstruction, impaired gas transfer, and increased pulmonary phagocytic activity. J. Infect. Dis. 185, 1326–1334. doi: 10.1086/339885
Aresha, M., Sanath, M., Deepika, F., Renu, W., Anura, B., Chanditha, H., et al. (2008). Genotyping of Plasmodium vivax infections in Sri Lanka using Pvmsp-3alpha and Pvcs genes as markers: a preliminary report. Trop. Biomed. 25, 100–106. Available online at: http://www.msptm.org/files/100_-_106_Aresha_Manamperi.pdf
Auburn, S., Bohme, U., Steinbiss, S., Trimarsanto, H., Hostetler, J., Sanders, M., et al. (2016a). A new Plasmodium vivax reference sequence with improved assembly of the subtelomeres reveals an abundance of pir genes. Wellcome Open Res 1, 4. doi: 10.12688/wellcomeopenres.9876.1
Auburn, S., Marfurt, J., Maslen, G., Campino, S., Ruano Rubio, V., Manske, M., et al. (2013). Effective preparation of Plasmodium vivax field isolates for high-throughput whole genome sequencing. PLoS ONE 8:e53160. doi: 10.1371/journal.pone.0053160
Auburn, S., Serre, D., Pearson, R. D., Amato, R., Sriprawat, K., To, S., et al. (2016b). Genomic analysis reveals a common breakpoint in amplifications of the Plasmodium vivax multidrug resistance 1 locus in Thailand. J. Infect. Dis. 214, 1235–1242. doi: 10.1093/infdis/jiw323
Bahl, A., Brunk, B., Crabtree, J., Fraunholz, M. J., Gajria, B., Grant, G. R., et al. (2003). PlasmoDB: the Plasmodium genome resource. A database integrating experimental and computational data. Nucleic Acids Res 31, 212–215. doi: 10.1093/nar/gkg081
Baniecki, M. L., Faust, A. L., Schaffner, S. F., Park, D. J., Galinsky, K., Daniels, R. F., et al. (2015). Development of a single nucleotide polymorphism barcode to genotype Plasmodium vivax infections. PLoS Negl. Trop. Dis. 9:e0003539. doi: 10.1371/journal.pntd.0003539
Barcus, M. J., Basri, H., Picarima, H., Manyakori, C., Sekartuti, Elyazar, I., et al. (2007). Demographic risk factors for severe and fatal vivax and falciparum malaria among hospital admissions in northeastern Indonesian Papua. Am. J. Trop. Med. Hyg. 77, 984–991. doi: 10.4269/ajtmh.2007.77.984
Barnwell, J. W., Ingravallo, P., Galinski, M. R., Matsumoto, Y., and Aikawa, M. (1990). Plasmodium vivax: malarial proteins associated with the membrane-bound caveola-vesicle complexes and cytoplasmic cleft structures of infected erythrocytes. Exp. Parasitol. 70, 85–99. doi: 10.1016/0014-4894(90)90088-T
Barry, A. E., Waltmann, A., Koepfli, C., Barnadas, C., and Mueller, I. (2015). Uncovering the transmission dynamics of Plasmodium vivax using population genetics. Pathog. Glob. Health 109, 142–152. doi: 10.1179/2047773215Y.0000000012
Battle, K. E., Gething, P. W., Elyazar, I. R., Moyes, C. L., Sinka, M. E., Howes, R. E., et al. (2012). The global public health significance of Plasmodium vivax. Adv. Parasitol. 80, 1–111. doi: 10.1016/B978-0-12-397900-1.00001-3
Bautista, J. M., Marin-Garcia, P., Diez, A., Azcarate, I. G., and Puyet, A. (2014). Malaria proteomics: insights into the parasite-host interactions in the pathogenic space. J. Proteomics 97, 107–125. doi: 10.1016/j.jprot.2013.10.011
Beutler, E., Duparc, S., and Group, G. P. D. W. (2007). Glucose-6-phosphate dehydrogenase deficiency and antimalarial drug development. Am. J. Trop. Med. Hyg. 77, 779–789. doi: 10.4269/ajtmh.2007.77.779
Bitoh, T., Fueda, K., Ohmae, H., Watanabe, M., and Ishikawa, H. (2011). Risk analysis of the re-emergence of Plasmodium vivax malaria in Japan using a stochastic transmission model. Environ. Health Prev. Med. 16, 171–177. doi: 10.1007/s12199-010-0184-8
Boopathi, P. A., Subudhi, A. K., Garg, S., Middha, S., Acharya, J., Pakalapati, D., et al. (2013). Revealing natural antisense transcripts from Plasmodium vivax isolates: evidence of genome regulation in complicated malaria. Infect. Genet. Evol. 20, 428–443. doi: 10.1016/j.meegid.2013.09.026
Boopathi, P. A., Subudhi, A. K., Garg, S., Middha, S., Acharya, J., Pakalapati, D., et al. (2014). Dataset of natural antisense transcripts in P. vivax clinical isolates derived using custom designed strand-specific microarray. Genom Data 2, 199–201. doi: 10.1016/j.gdata.2014.06.024
Boopathi, P. A., Subudhi, A. K., Middha, S., Acharya, J., Mugasimangalam, R. C., Kochar, S. K., et al. (2016). Design, construction and validation of a Plasmodium vivax microarray for the transcriptome profiling of clinical isolates. Acta Trop. 164, 438–447. doi: 10.1016/j.actatropica.2016.10.002
Bousema, T., and Drakeley, C. (2011). Epidemiology and infectivity of Plasmodium falciparum and Plasmodium vivax gametocytes in relation to malaria control and elimination. Clin. Microbiol. Rev. 24, 377–410. doi: 10.1128/CMR.00051-10
Boyd, M. F., and Kitchen, S. F. (1937). On the infectiousness of patients infected with Plasmodium vivax and Plasmodium falciparum. Am. J. Trop. Med. Hyg. 17, 253–262. doi: 10.4269/ajtmh.1937.s1-17.253
Bozdech, Z., Llinas, M., Pulliam, B. L., Wong, E. D., Zhu, J., and DeRisi, J. L. (2003). The transcriptome of the intraerythrocytic developmental cycle of Plasmodium falciparum. PLoS Biol. 1:E5. doi: 10.1371/journal.pbio.0000005
Bozdech, Z., Mok, S., Hu, G., Imwong, M., Jaidee, A., Russell, B., et al. (2008). The transcriptome of Plasmodium vivax reveals divergence and diversity of transcriptional regulation in malaria parasites. Proc. Natl. Acad. Sci. U.S.A. 105, 16290–16295. doi: 10.1073/pnas.0807404105
Brazeau, N. F., Hathaway, N., Parobek, C. M., Lin, J. T., Bailey, J. A., Lon, C., et al. (2016). Longitudinal pooled deep sequencing of the Plasmodium vivax K12 kelch gene in cambodia reveals a lack of selection by artemisinin. Am. J. Trop. Med. Hyg. 95, 1409–1412. doi: 10.4269/ajtmh.16-0566
Bright, A. T., Manary, M. J., Tewhey, R., Arango, E. M., Wang, T., Schork, N. J., et al. (2014). A high resolution case study of a patient with recurrent Plasmodium vivax infections shows that relapses were caused by meiotic siblings. PLoS Negl. Trop. Dis. 8:e2882. doi: 10.1371/journal.pntd.0002882
Bright, A. T., Tewhey, R., Abeles, S., Chuquiyauri, R., Llanos-Cuentas, A., Ferreira, M. U., et al. (2012). Whole genome sequencing analysis of Plasmodium vivax using whole genome capture. BMC Genomics 13:262. doi: 10.1186/1471-2164-13-262
Carlton, J. M., Adams, J. H., Silva, J. C., Bidwell, S. L., Lorenzi, H., Caler, E., et al. (2008a). Comparative genomics of the neglected human malaria parasite Plasmodium vivax. Nature 455, 757–763. doi: 10.1038/nature07327
Carlton, J. M., Angiuoli, S. V., Suh, B. B., Kooij, T. W., Pertea, M., Silva, J. C., et al. (2002). Genome sequence and comparative analysis of the model rodent malaria parasite Plasmodium yoelii yoelii. Nature 419, 512–519. doi: 10.1038/nature01099
Chan, E. R., Barnwell, J. W., Zimmerman, P. A., and Serre, D. (2015). Comparative analysis of field-isolate and monkey-adapted Plasmodium vivax genomes. PLoS Negl. Trop. Dis. 9:e0003566. doi: 10.1371/journal.pntd.0003566
Chan, E. R., Menard, D., David, P. H., Ratsimbasoa, A., Kim, S., Chim, P., et al. (2012). Whole genome sequencing of field isolates provides robust characterization of genetic diversity in Plasmodium vivax. PLoS Negl. Trop. Dis. 6:e1811. doi: 10.1371/journal.pntd.0001811
Chang, H. H., Park, D. J., Galinsky, K. J., Schaffner, S. F., Ndiaye, D., Ndir, O., et al. (2012). Genomic sequencing of Plasmodium falciparum malaria parasites from Senegal reveals the demographic history of the population. Mol. Biol. Evol. 29, 3427–3439. doi: 10.1093/molbev/mss161
Cheeseman, I. H., Miller, B., Tan, J. C., Tan, A., Nair, S., Nkhoma, S. C., et al. (2016). Population structure shapes copy number variation in malaria parasites. Mol. Biol. Evol. 33, 603–620. doi: 10.1093/molbev/msv282
Chen, J. H., Jung, J. W., Wang, Y., Ha, K. S., Lu, F., Lim, C. S., et al. (2010). Immunoproteomics profiling of blood stage Plasmodium vivax infection by high-throughput screening assays. J. Proteome Res. 9, 6479–6489. doi: 10.1021/pr100705g
Chung, B. H., Lee, S. W., Lee, S. E., Hwang, T. J., and Shin, H. S. (2008). Predictors of Plasmodium vivax malaria-induced nephropathy in young Korean men. Nephron Clin. Pract. 110, c172–c177. doi: 10.1159/000167023
Collins, W. E., Contacos, P. G., Krotoski, W. A., and Howard, W. A. (1972). Transmission of four Central American strains of Plasmodium vivax from monkey to man. J. Parasitol. 58, 332–335. doi: 10.2307/3278097
Cornejo, O. E., Fisher, D., and Escalante, A. A. (2015). Genome-wide patterns of genetic polymorphism and signatures of selection in Plasmodium vivax. Genome Biol. Evol. 7, 106–119. doi: 10.1093/gbe/evu267
Coura, J. R., Suarez-Mutis, M., and Ladeia-Andrade, S. (2006). A new challenge for malaria control in Brazil: asymptomatic Plasmodium infection–a review. Mem. Inst. Oswaldo Cruz 101, 229–237. doi: 10.1590/S0074-02762006000300001
Cui, L., Mascorro, C. N., Fan, Q., Rzomp, K. A., Khuntirat, B., Zhou, G., et al. (2003). Genetic diversity and multiple infections of Plasmodium vivax malaria in Western Thailand. Am. J. Trop. Med. Hyg. 68, 613–619. doi: 10.4269/ajtmh.2003.68.613
Cui, L., Trongnipatt, N., Sattabongkot, J., and Udomsangpetch, R. (2009). Culture of exoerythrocytic stages of the malaria parasites Plasmodium falciparum and Plasmodium vivax. Methods Mol. Biol. 470, 263–273. doi: 10.1007/978-1-59745-204-5_18
Culleton, R., Coban, C., Zeyrek, F. Y., Cravo, P., Kaneko, A., Randrianarivelojosia, M., et al. (2011). The origins of African Plasmodium vivax; insights from mitochondrial genome sequencing. PLoS ONE 6:e29137. doi: 10.1371/journal.pone.0029137
de Araujo, F. C., de Rezende, A. M., Fontes, C. J., Carvalho, L. H., and Alves de Brito, C. F. (2012). Multiple-clone activation of hypnozoites is the leading cause of relapse in Plasmodium vivax infection. PLoS ONE 7:e49871. doi: 10.1371/journal.pone.0049871
de Lacerda, M. V., Zackiewicz, C., Alecrim, W. D., and Alecrim, M. (2007). The neglected Plasmodium vivax: are researchers from endemic areas really concerned about new treatment options? Rev. Soc. Bras. Med. Trop. 40, 489–490. doi: 10.1590/S0037-86822007000400026
del Portillo, H. A., Gysin, J., Mattei, D. M., Khouri, E., Udagama, P. V., Mendis, K. N., et al. (1988). Plasmodium vivax: cloning and expression of a major blood-stage surface antigen. Exp. Parasitol. 67, 346–353. doi: 10.1016/0014-4894(88)90081-1
del Portillo, H. A., Lanzer, M., Rodriguez-Malaga, S., Zavala, F., and Fernandez-Becerra, C. (2004). Variant genes and the spleen in Plasmodium vivax malaria. Int. J. Parasitol. 34, 1547–1554. doi: 10.1016/j.ijpara.2004.10.012
Dembele, L., Franetich, J. F., Lorthiois, A., Gego, A., Zeeman, A. M., Kocken, C. H., et al. (2014). Persistence and activation of malaria hypnozoites in long-term primary hepatocyte cultures. Nat. Med. 20, 307–312. doi: 10.1038/nm.3461
de Santana Filho, F. S., Arcanjo, A. R., Chehuan, Y. M., Costa, M. R., Martinez-Espinosa, F. E., Vieira, J. L., et al. (2007). Chloroquine-resistant Plasmodium vivax, Brazilian Amazon. Emerging Infect. Dis. 13, 1125–1126. doi: 10.3201/eid1307.061386
Dharia, N. V., Bright, A. T., Westenberger, S. J., Barnes, S. W., Batalov, S., Kuhen, K., et al. (2010). Whole-genome sequencing and microarray analysis of ex vivo Plasmodium vivax reveal selective pressure on putative drug resistance genes. Proc. Natl. Acad. Sci. U.S.A. 107, 20045–20050. doi: 10.1073/pnas.1003776107
Duval, L., Fourment, M., Nerrienet, E., Rousset, D., Sadeuh, S. A., Goodman, S. M., et al. (2010). African apes as reservoirs of Plasmodium falciparum and the origin and diversification of the Laverania subgenus. Proc. Natl. Acad. Sci. U.S.A. 107, 10561–10566. doi: 10.1073/pnas.1005435107
Erhart, L. M., Yingyuen, K., Chuanak, N., Buathong, N., Laoboonchai, A., Miller, R. S., et al. (2004). Hematologic and clinical indices of malaria in a semi-immune population of western Thailand. Am. J. Trop. Med. Hyg. 70, 8–14. doi: 10.4269/ajtmh.2004.70.8
Escalante, A. A., Ferreira, M. U., Vinetz, J. M., Volkman, S. K., Cui, L., Gamboa, D., et al. (2015). Malaria molecular epidemiology: lessons from the international centers of excellence for malaria research network. Am. J. Trop. Med. Hyg. 93, 79–86. doi: 10.4269/ajtmh.15-0005
Feng, X., Carlton, J. M., Joy, D. A., Mu, J., Furuya, T., Suh, B. B., et al. (2003). Single-nucleotide polymorphisms and genome diversity in Plasmodium vivax. Proc. Natl. Acad. Sci. U.S.A. 100, 8502–8507. doi: 10.1073/pnas.1232502100
Fernandez-Becerra, C., Pein, O., de Oliveira, T. R., Yamamoto, M. M., Cassola, A. C., Rocha, C., et al. (2005). Variant proteins of Plasmodium vivax are not clonally expressed in natural infections. Mol. Microbiol. 58, 648–658. doi: 10.1111/j.1365-2958.2005.04850.x
Fernandez-Becerra, C., Pinazo, M. J., Gonzalez, A., Alonso, P. L., del Portillo, H. A., and Gascon, J. (2009). Increased expression levels of the pvcrt-o and pvmdr1 genes in a patient with severe Plasmodium vivax malaria. Malar. J. 8:55. doi: 10.1186/1475-2875-8-55
Field, J., and Shute, P. (1956). “Plasmodium vivax,” in The Microscopuic Diagnosis of Human Malaria. II. A Morphological Study of the Erythrocytic Parasites, eds Institute for Medical Research. Studies from the Institute for Medical Research, Federated Malay States (Kuala Lumpur: Govt. Press), 251.
Flannery, E. L., Wang, T., Akbari, A., Corey, V. C., Gunawan, F., Bright, A. T., et al. (2015). Next-generation sequencing of Plasmodium vivax patient samples shows evidence of direct evolution in drug-resistance genes. ACS Infect Dis 1, 367–379. doi: 10.1021/acsinfecdis.5b00049
Foth, B. J., Ralph, S. A., Tonkin, C. J., Struck, N. S., Fraunholz, M., Roos, D. S., et al. (2003). Dissecting apicoplast targeting in the malaria parasite Plasmodium falciparum. Science 299, 705–708. doi: 10.1126/science.1078599
Frazer, K. A., Elnitski, L., Church, D. M., Dubchak, I., and Hardison, R. C. (2003). Cross-species sequence comparisons: a review of methods and available resources. Genome Res. 13, 1–12. doi: 10.1101/gr.222003
Frech, C., and Chen, N. (2011). Genome comparison of human and non-human malaria parasites reveals species subset-specific genes potentially linked to human disease. PLoS Comput. Biol. 7:e1002320. doi: 10.1371/journal.pcbi.1002320
Friedrich, L. R., Popovici, J., Kim, S., Dysoley, L., Zimmerman, P. A., Menard, D., et al. (2016). Complexity of infection and genetic diversity in cambodian Plasmodium vivax. PLoS Negl. Trop. Dis. 10:e0004526. doi: 10.1371/journal.pntd.0004526
Galinski, M. R., and Barnwell, J. W. (2012). “Chapter 5 - Nonhuman Primate Models for Human Malaria Research A2 - Abee, Christian R,” in Nonhuman Primates in Biomedical Research, 2nd Edn., eds K. Mansfield, S. Tardif and T. Morris (Boston, MA: Academic Press), 299–323.
Gardner, M. J., Hall, N., Fung, E., White, O., Berriman, M., Hyman, R. W., et al. (2002). Genome sequence of the human malaria parasite Plasmodium falciparum. Nature 419, 498–511. doi: 10.1038/nature01097
Genton, B., D'Acremont, V., Rare, L., Baea, K., Reeder, J. C., Alpers, M. P., et al. (2008). Plasmodium vivax and mixed infections are associated with severe malaria in children: a prospective cohort study from Papua New Guinea. PLoS Med. 5:e127. doi: 10.1371/journal.pmed.0050127
Gething, P. W., Elyazar, I. R., Moyes, C. L., Smith, D. L., Battle, K. E., Guerra, C. A., et al. (2012). A long neglected world malaria map: Plasmodium vivax endemicity in 2010. PLoS Negl. Trop. Dis. 6:e1814. doi: 10.1371/journal.pntd.0001814
Guerra, C. A., Howes, R. E., Patil, A. P., Gething, P. W., Van Boeckel, T. P., Temperley, W. H., et al. (2010). The international limits and population at risk of Plasmodium vivax transmission in 2009. PLoS Negl. Trop. Dis. 4:e774. doi: 10.1371/journal.pntd.0000774
Gunawardena, S., Karunaweera, N. D., Ferreira, M. U., Phone-Kyaw, M., Pollack, R. J., Alifrangis, M., et al. (2010). Geographic structure of Plasmodium vivax: microsatellite analysis of parasite populations from Sri Lanka, Myanmar, and Ethiopia. Am. J. Trop. Med. Hyg. 82, 235–242. doi: 10.4269/ajtmh.2010.09-0588
Gupta, B. K. (2010). Reply to comment: left upper quadrant abdominal pain in malaria: suspect pathological splenic rupture first. Trans. R. Soc. Trop. Med. Hyg. 104, 690. doi: 10.1016/j.trstmh.2010.07.001
Handayani, S., Chiu, D. T., Tjitra, E., Kuo, J. S., Lampah, D., Kenangalem, E., et al. (2009). High deformability of Plasmodium vivax-infected red blood cells under microfluidic conditions. J. Infect. Dis. 199, 445–450. doi: 10.1086/596048
Hemingway, J., Shretta, R., Wells, T. N., Bell, D., Djimde, A. A., Achee, N., et al. (2016). Tools and strategies for malaria control and elimination: what do we need to achieve a grand convergence in malaria? PLoS Biol. 14:e1002380. doi: 10.1371/journal.pbio.1002380
Hemmer, C. J., Holst, F. G., Kern, P., Chiwakata, C. B., Dietrich, M., and Reisinger, E. C. (2006). Stronger host response per parasitized erythrocyte in Plasmodium vivax or ovale than in Plasmodium falciparum malaria. Trop. Med. Int. Health 11, 817–823. doi: 10.1111/j.1365-3156.2006.01635.x
Hester, J., Chan, E. R., Menard, D., Mercereau-Puijalon, O., Barnwell, J., Zimmerman, P. A., et al. (2013). De novo assembly of a field isolate genome reveals novel Plasmodium vivax erythrocyte invasion genes. PLoS Negl. Trop. Dis. 7:e2569. doi: 10.1371/journal.pntd.0002569
Hollingdale, M. R., Collins, W. E., Campbell, C. C., and Schwartz, A. L. (1985). In vitro culture of two populations (dividing and nondividing) of exoerythrocytic parasites of Plasmodium vivax. Am. J. Trop. Med. Hyg. 34, 216–222. doi: 10.4269/ajtmh.1985.34.216
Hoo, R., Zhu, L., Amaladoss, A., Mok, S., Natalang, O., Lapp, S. A., et al. (2016). Integrated analysis of the Plasmodium species transcriptome. EBioMedicine 7, 255–266. doi: 10.1016/j.ebiom.2016.04.011
Hupalo, D. N., Luo, Z., Melnikov, A., Sutton, P. L., Rogov, P., Escalante, A., et al. (2016). Population genomics studies identify signatures of global dispersal and drug resistance in Plasmodium vivax. Nat. Genet. 48, 953–958. doi: 10.1038/ng.3588
Hussain, M. M., Sohail, M., Abhishek, K., and Raziuddin, M. (2013). Investigation on Plasmodium falciparum and Plasmodium vivax infection influencing host haematological factors in tribal dominant and malaria endemic population of Jharkhand. Saudi J. Biol. Sci. 20, 195–203. doi: 10.1016/j.sjbs.2013.01.003
Imwong, M., Nair, S., Pukrittayakamee, S., Sudimack, D., Williams, J. T., Mayxay, M., et al. (2007a). Contrasting genetic structure in Plasmodium vivax populations from Asia and South America. Int. J. Parasitol. 37, 1013–1022. doi: 10.1016/j.ijpara.2007.02.010
Imwong, M., Pukrittayakamee, S., Renia, L., Letourneur, F., Charlieu, J. P., Leartsakulpanich, U., et al. (2003). Novel point mutations in the dihydrofolate reductase gene of Plasmodium vivax: evidence for sequential selection by drug pressure. Antimicrob. Agents Chemother. 47, 1514–1521. doi: 10.1128/AAC.47.5.1514-1521.2003
Imwong, M., Snounou, G., Pukrittayakamee, S., Tanomsing, N., Kim, J. R., Nandy, A., et al. (2007b). Relapses of Plasmodium vivax infection usually result from activation of heterologous hypnozoites. J. Infect. Dis. 195, 927–933. doi: 10.1086/512241
Iriko, H., Jin, L., Kaneko, O., Takeo, S., Han, E. T., Tachibana, M., et al. (2009). A small-scale systematic analysis of alternative splicing in Plasmodium falciparum. Parasitol. Int. 58, 196–199. doi: 10.1016/j.parint.2009.02.002
Jongwutiwes, S., Putaporntip, C., Iwasaki, T., Ferreira, M. U., Kanbara, H., and Hughes, A. L. (2005). Mitochondrial genome sequences support ancient population expansion in Plasmodium vivax. Mol. Biol. Evol. 22, 1733–1739. doi: 10.1093/molbev/msi168
Joy, D. A., Gonzalez-Ceron, L., Carlton, J. M., Gueye, A., Fay, M., McCutchan, T. F., et al. (2008). Local adaptation and vector-mediated population structure in Plasmodium vivax malaria. Mol. Biol. Evol. 25, 1245–1252. doi: 10.1093/molbev/msn073
Judice, C. C., Bourgard, C., Kayano, A. C., Albrecht, L., and Costa, F. T. (2016). MicroRNAs in the host-apicomplexan parasites interactions: a review of immunopathological aspects. Front. Cell. Infect. Microbiol. 6:5. doi: 10.3389/fcimb.2016.00005
Karunaweera, N. D., Ferreira, M. U., Munasinghe, A., Barnwell, J. W., Collins, W. E., King, C. L., et al. (2008). Extensive microsatellite diversity in the human malaria parasite Plasmodium vivax. Gene 410, 105–112. doi: 10.1016/j.gene.2007.11.022
Karunaweera, N. D., Grau, G. E., Gamage, P., Carter, R., and Mendis, K. N. (1992). Dynamics of fever and serum levels of tumor necrosis factor are closely associated during clinical paroxysms in Plasmodium vivax malaria. Proc. Natl. Acad. Sci. U.S.A. 89, 3200–3203. doi: 10.1073/pnas.89.8.3200
Kim, A., Popovici, J., Vantaux, A., Samreth, R., Bin, S., Kim, S., et al. (2017). Characterization of P. vivax blood stage transcriptomes from field isolates reveals similarities among infections and complex gene isoforms. Sci. Rep. 7:7761. doi: 10.1038/s41598-017-07275-9
Kim, H. C., Pacha, L. A., Lee, W. J., Lee, J. K., Gaydos, J. C., Sames, W. J., et al. (2009). Malaria in the Republic of Korea, 1993-2007. Variables related to re-emergence and persistence of Plasmodium vivax among Korean populations and U.S. forces in Korea. Mil. Med. 174, 762–769. doi: 10.7205/MILMED-D-01-6208
Kiszewski, A., Mellinger, A., Spielman, A., Malaney, P., Sachs, S. E., and Sachs, J. (2004). A global index representing the stability of malaria transmission. Am. J. Trop. Med. Hyg. 70, 486–498. doi: 10.4269/ajtmh.2004.70.486
Kochar, D. K., Das, A., Kochar, S. K., Saxena, V., Sirohi, P., Garg, S., et al. (2009a). Severe Plasmodium vivax malaria: a report on serial cases from Bikaner in northwestern India. Am. J. Trop. Med. Hyg. 80, 194–198. doi: 10.4269/ajtmh.2009.80.194
Kochar, S. K., Mahajan, M., Gupta, R. P., Middha, S., Acharya, J., Kochar, A., et al. (2009b). Acute attack of AIP (acute intermittent porphyria) with severe vivax malaria associated with convulsions: a case report. J. Vector Borne Dis. 46, 307–309.
Koepfli, C., Rodrigues, P. T., Antao, T., Orjuela-Sanchez, P., Van den Eede, P., Gamboa, D., et al. (2015). Plasmodium vivax diversity and population structure across four continents. PLoS Negl. Trop. Dis. 9:e0003872. doi: 10.1371/journal.pntd.0003872
Korsinczky, M., Fischer, K., Chen, N., Baker, J., Rieckmann, K., and Cheng, Q. (2004). Sulfadoxine resistance in Plasmodium vivax is associated with a specific amino acid in dihydropteroate synthase at the putative sulfadoxine-binding site. Antimicrob. Agents Chemother. 48, 2214–2222. doi: 10.1128/AAC.48.6.2214-2222.2004
Krotoski, W. A., Collins, W. E., Bray, R. S., Garnham, P. C., Cogswell, F. B., Gwadz, R. W., et al. (1982). Demonstration of hypnozoites in sporozoite-transmitted Plasmodium vivax infection. Am. J. Trop. Med. Hyg. 31, 1291–1293. doi: 10.4269/ajtmh.1982.31.1291
Lacerda, M. V., Alexandre, M. A., Santos, P. D., Arcanjo, A. R., Alecrim, W. D., and Alecrim, M. G. (2004). Idiopathic thrombocytopenic purpura due to vivax malaria in the Brazilian Amazon. Acta Trop. 90, 187–190. doi: 10.1016/j.actatropica.2003.12.001
Lacerda, M. V., Fragoso, S. C., Alecrim, M. G., Alexandre, M. A., Magalhaes, B. M., Siqueira, A. M., et al. (2012a). Postmortem characterization of patients with clinical diagnosis of Plasmodium vivax malaria: to what extent does this parasite kill? Clin. Infect. Dis. 55, e67–e74. doi: 10.1093/cid/cis615
Lacerda, M. V., Mourao, M. P., Alexandre, M. A., Siqueira, A. M., Magalhaes, B. M., Martinez-Espinosa, F. E., et al. (2012b). Understanding the clinical spectrum of complicated Plasmodium vivax malaria: a systematic review on the contributions of the Brazilian literature. Malar. J. 11:12. doi: 10.1186/1475-2875-11-12
Lampah, D. A., Yeo, T. W., Hardianto, S. O., Tjitra, E., Kenangalem, E., Sugiarto, P., et al. (2011). Coma associated with microscopy-diagnosed Plasmodium vivax: a prospective study in Papua, Indonesia. PLoS Negl. Trop. Dis. 5:e1032. doi: 10.1371/journal.pntd.0001032
Lanca, E. F., Magalhaes, B. M., Vitor-Silva, S., Siqueira, A. M., Benzecry, S. G., Alexandre, M. A., et al. (2012). Risk factors and characterization of Plasmodium vivax-associated admissions to pediatric intensive care units in the Brazilian Amazon. PLoS ONE 7:e35406. doi: 10.1371/journal.pone.0035406
Le Roch, K. G., Zhou, Y., Blair, P. L., Grainger, M., Moch, J. K., Haynes, J. D., et al. (2003). Discovery of gene function by expression profiling of the malaria parasite life cycle. Science 301, 1503–1508. doi: 10.1126/science.1087025
Lin, J. T., Hathaway, N. J., Saunders, D. L., Lon, C., Balasubramanian, S., Kharabora, O., et al. (2015). Using amplicon deep sequencing to detect genetic signatures of Plasmodium vivax relapse. J. Infect. Dis. 212, 999–1008. doi: 10.1093/infdis/jiv142
Lin, J. T., Patel, J. C., Kharabora, O., Sattabongkot, J., Muth, S., Ubalee, R., et al. (2013). Plasmodium vivax isolates from Cambodia and Thailand show high genetic complexity and distinct patterns of P. vivax multidrug resistance gene 1 (pvmdr1) polymorphisms. Am. J. Trop. Med. Hyg. 88, 1116–1123. doi: 10.4269/ajtmh.12-0701
Lu, F., Li, J., Wang, B., Cheng, Y., Kong, D. H., Cui, L., et al. (2014). Profiling the humoral immune responses to Plasmodium vivax infection and identification of candidate immunogenic rhoptry-associated membrane antigen (RAMA). J. Proteomics 102, 66–82. doi: 10.1016/j.jprot.2014.02.029
Margulies, M., Egholm, M., Altman, W. E., Attiya, S., Bader, J. S., Bemben, L. A., et al. (2005). Genome sequencing in microfabricated high-density picolitre reactors. Nature 437, 376–380. doi: 10.1038/nature03959
Marsh, K., Forster, D., Waruiru, C., Mwangi, I., Winstanley, M., Marsh, V., et al. (1995). Indicators of life-threatening malaria in African children. N. Engl. J. Med. 332, 1399–1404. doi: 10.1056/NEJM199505253322102
McGready, R., Davison, B. B., Stepniewska, K., Cho, T., Shee, H., Brockman, A., et al. (2004). The effects of Plasmodium falciparum and P. vivax infections on placental histopathology in an area of low malaria transmission. Am. J. Trop. Med. Hyg. 70, 398–407. doi: 10.4269/ajtmh.2004.70.398
Mehlotra, R. K., Blankenship, D., Howes, R. E., Rakotomanga, T. A., Ramiranirina, B., Ramboarina, S., et al. (2017). Long-term in vitro culture of Plasmodium vivax isolates from Madagascar maintained in Saimiri boliviensis blood. Malar. J. 16, 442. doi: 10.1186/s12936-017-2090-7
Melnikov, A., Galinsky, K., Rogov, P., Fennell, T., Van Tyne, D., Russ, C., et al. (2011). Hybrid selection for sequencing pathogen genomes from clinical samples. Genome Biol. 12:R73. doi: 10.1186/gb-2011-12-8-r73
Melo, G. C., Monteiro, W. M., Siqueira, A. M., Silva, S. R., Magalhaes, B. M., Alencar, A. C., et al. (2014). Expression levels of pvcrt-o and pvmdr-1 are associated with chloroquine resistance and severe Plasmodium vivax malaria in patients of the Brazilian Amazon. PLoS ONE 9:e105922. doi: 10.1371/journal.pone.0105922
Menard, D., Barnadas, C., Bouchier, C., Henry-Halldin, C., Gray, L. R., Ratsimbasoa, A., et al. (2010). Plasmodium vivax clinical malaria is commonly observed in Duffy-negative Malagasy people. Proc. Natl. Acad. Sci. U.S.A. 107, 5967–5971. doi: 10.1073/pnas.0912496107
Menard, D., Chan, E. R., Benedet, C., Ratsimbasoa, A., Kim, S., Chim, P., et al. (2013). Whole genome sequencing of field isolates reveals a common duplication of the Duffy binding protein gene in Malagasy Plasmodium vivax strains. PLoS Negl. Trop. Dis. 7:e2489. doi: 10.1371/journal.pntd.0002489
Mideo, N., Kennedy, D. A., Carlton, J. M., Bailey, J. A., Juliano, J. J., and Read, A. F. (2013). Ahead of the curve: next generation estimators of drug resistance in malaria infections. Trends Parasitol. 29, 321–328. doi: 10.1016/j.pt.2013.05.004
Mikolajczak, S. A., Vaughan, A. M., Kangwanrangsan, N., Roobsoong, W., Fishbaugher, M., Yimamnuaychok, N., et al. (2015). Plasmodium vivax liver stage development and hypnozoite persistence in human liver-chimeric mice. Cell Host Microbe 17, 526–535. doi: 10.1016/j.chom.2015.02.011
Moon, S. U., Lee, H. W., Kim, J. Y., Na, B. K., Cho, S. H., Lin, K., et al. (2009). High frequency of genetic diversity of Plasmodium vivax field isolates in Myanmar. Acta Trop. 109, 30–36. doi: 10.1016/j.actatropica.2008.09.006
Moreno-Perez, D. A., Degano, R., Ibarrola, N., Muro, A., and Patarroyo, M. A. (2014). Determining the Plasmodium vivax VCG-1 strain blood stage proteome. J. Proteomics 113C, 268–280. doi: 10.1016/j.jprot.2014.10.003
Mu, J., Joy, D. A., Duan, J., Huang, Y., Carlton, J., Walker, J., et al. (2005). Host switch leads to emergence of Plasmodium vivax malaria in humans. Mol. Biol. Evol. 22, 1686–1693. doi: 10.1093/molbev/msi160
Mueller, I., Galinski, M. R., Baird, J. K., Carlton, J. M., Kochar, D. K., Alonso, P. L., et al. (2009). Key gaps in the knowledge of Plasmodium vivax, a neglected human malaria parasite. Lancet Infect. Dis. 9, 555–566. doi: 10.1016/S1473-3099(09)70177-X
Mueller, I., Kaiok, J., Reeder, J. C., and Cortes, A. (2002). The population structure of Plasmodium falciparum and Plasmodium vivax during an epidemic of malaria in the Eastern Highlands of Papua New Guinea. Am. J. Trop. Med. Hyg. 67, 459–464. doi: 10.4269/ajtmh.2002.67.459
Naing, C., Whittaker, M. A., Nyunt Wai, V., and Mak, J. W. (2014). Is Plasmodium vivax malaria a severe malaria? A systematic review and meta-analysis. PLoS Negl. Trop. Dis. 8:e3071. doi: 10.1371/journal.pntd.0003071
Nair, S., Nkhoma, S. C., Serre, D., Zimmerman, P. A., Gorena, K., Daniel, B. J., et al. (2014). Single-cell genomics for dissection of complex malaria infections. Genome Res. 24, 1028–1038. doi: 10.1101/gr.168286.113
Neafsey, D. E., Galinsky, K., Jiang, R. H., Young, L., Sykes, S. M., Saif, S., et al. (2012). The malaria parasite Plasmodium vivax exhibits greater genetic diversity than Plasmodium falciparum. Nat. Genet. 44, 1046–1050. doi: 10.1038/ng.2373
Noulin, F., Borlon, C., Van Den Abbeele, J., D'Alessandro, U., and Erhart, A. (2013). 1912-2012: a century of research on Plasmodium vivax in vitro culture. Trends Parasitol. 29, 286–294. doi: 10.1016/j.pt.2013.03.012
Orjuela-Sanchez, P., Sa, J. M., Brandi, M. C., Rodrigues, P. T., Bastos, M. S., Amaratunga, C., et al. (2013). Higher microsatellite diversity in Plasmodium vivax than in sympatric Plasmodium falciparum populations in Pursat, Western Cambodia. Exp. Parasitol. 134, 318–326. doi: 10.1016/j.exppara.2013.03.029
Otto, T. D., Wilinski, D., Assefa, S., Keane, T. M., Sarry, L. R., Bohme, U., et al. (2010). New insights into the blood-stage transcriptome of Plasmodium falciparum using RNA-Seq. Mol. Microbiol. 76, 12–24. doi: 10.1111/j.1365-2958.2009.07026.x
Pain, A., Bohme, U., Berry, A. E., Mungall, K., Finn, R. D., Jackson, A. P., et al. (2008). The genome of the simian and human malaria parasite Plasmodium knowlesi. Nature 455, 799–803. doi: 10.1038/nature07306
Panichakul, T., Sattabongkot, J., Chotivanich, K., Sirichaisinthop, J., Cui, L., and Udomsangpetch, R. (2007). Production of erythropoietic cells in vitro for continuous culture of Plasmodium vivax. Int. J. Parasitol. 37, 1551–1557. doi: 10.1016/j.ijpara.2007.05.009
Parobek, C. M., Lin, J. T., Saunders, D. L., Barnett, E. J., Lon, C., Lanteri, C. A., et al. (2016). Selective sweep suggests transcriptional regulation may underlie Plasmodium vivax resilience to malaria control measures in Cambodia. Proc. Natl. Acad. Sci. U.S.A. 113, E8096–E8105. doi: 10.1073/pnas.1608828113
Parroche, P., Lauw, F. N., Goutagny, N., Latz, E., Monks, B. G., Visintin, A., et al. (2007). Malaria hemozoin is immunologically inert but radically enhances innate responses by presenting malaria DNA to Toll-like receptor 9. Proc. Natl. Acad. Sci. U.S.A. 104, 1919–1924. doi: 10.1073/pnas.0608745104
Pava, Z., Handayuni, I., Wirjanata, G., To, S., Trianty, L., Noviyanti, R., et al. (2015). Expression of Plasmodium vivax crt-o is related to parasite stage but not ex vivo chloroquine susceptibility. Antimicrob. Agents Chemother. 60, 361–367. doi: 10.1128/AAC.02207-15
Pearson, R. D., Amato, R., Auburn, S., Miotto, O., Almagro-Garcia, J., Amaratunga, C., et al. (2016). Genomic analysis of local variation and recent evolution in Plasmodium vivax. Nat. Genet. 48, 959–964. doi: 10.1038/ng.3599
Pelle, K. G., Oh, K., Buchholz, K., Narasimhan, V., Joice, R., Milner, D. A., et al. (2015). Transcriptional profiling defines dynamics of parasite tissue sequestration during malaria infection. Genome Med. 7, 19. doi: 10.1186/s13073-015-0133-7
Poespoprodjo, J. R., Fobia, W., Kenangalem, E., Lampah, D. A., Warikar, N., Seal, A., et al. (2008). Adverse pregnancy outcomes in an area where multidrug-resistant plasmodium vivax and Plasmodium falciparum infections are endemic. Clin. Infect. Dis. 46, 1374–1381. doi: 10.1086/586743
Price, R. N., Douglas, N. M., and Anstey, N. M. (2009). New developments in Plasmodium vivax malaria: severe disease and the rise of chloroquine resistance. Curr. Opin. Infect. Dis. 22, 430–435. doi: 10.1097/QCO.0b013e32832f14c1
Price, R. N., von Seidlein, L., Valecha, N., Nosten, F., Baird, J. K., and White, N. J. (2014). Global extent of chloroquine-resistant Plasmodium vivax: a systematic review and meta-analysis. Lancet Infect. Dis. 14, 982–991. doi: 10.1016/S1473-3099(14)70855-2
Quintero, J. P., Siqueira, A. M., Tobon, A., Blair, S., Moreno, A., Arevalo-Herrera, M., et al. (2011). Malaria-related anaemia: a Latin American perspective. Mem. Inst. Oswaldo Cruz 106(Suppl. 1), 91–104. doi: 10.1590/S0074-02762011000900012
Quispe, A. M., Pozo, E., Guerrero, E., Durand, S., Baldeviano, G. C., Edgel, K. A., et al. (2014). Plasmodium vivax hospitalizations in a monoendemic malaria region: severe vivax malaria? Am. J. Trop. Med. Hyg. 91, 11–17. doi: 10.4269/ajtmh.12-0610
Rahimi, B. A., Thakkinstian, A., White, N. J., Sirivichayakul, C., Dondorp, A. M., and Chokejindachai, W. (2014). Severe vivax malaria: a systematic review and meta-analysis of clinical studies since 1900. Malar. J. 13:481. doi: 10.1186/1475-2875-13-481
Ralph, S. A., van Dooren, G. G., Waller, R. F., Crawford, M. J., Fraunholz, M. J., Foth, B. J., et al. (2004). Tropical infectious diseases: metabolic maps and functions of the Plasmodium falciparum apicoplast. Nat. Rev. Microbiol. 2, 203–216. doi: 10.1038/nrmicro843
Ray, S., Kamath, K. S., Srivastava, R., Raghu, D., Gollapalli, K., Jain, R., et al. (2012a). Serum proteome analysis of vivax malaria: an insight into the disease pathogenesis and host immune response. J. Proteomics 75, 3063–3080. doi: 10.1016/j.jprot.2011.10.018
Ray, S., Patel, S. K., Venkatesh, A., Bhave, A., Kumar, V., Singh, V., et al. (2016). Clinicopathological analysis and multipronged quantitative proteomics reveal oxidative stress and cytoskeletal proteins as possible markers for severe vivax malaria. Sci. Rep. 6:24557. doi: 10.1038/srep24557
Ray, S., Patel, S. K., Venkatesh, A., Chatterjee, G., Ansari, N. N., Gogtay, N. J., et al. (2017). Quantitative proteomics analysis of Plasmodium vivax induced alterations in human serum during the acute and convalescent phases of infection. Sci. Rep. 7:4400. doi: 10.1038/s41598-017-04447-5
Ray, S., Renu, D., Srivastava, R., Gollapalli, K., Taur, S., Jhaveri, T., et al. (2012b). Proteomic investigation of falciparum and vivax malaria for identification of surrogate protein markers. PLoS ONE 7:e41751. doi: 10.1371/journal.pone.0041751
Restrepo, E., Imwong, M., Rojas, W., Carmona-Fonseca, J., and Maestre, A. (2011). High genetic polymorphism of relapsing P. vivax isolates in northwest Colombia. Acta Trop. 119, 23–29. doi: 10.1016/j.actatropica.2011.03.012
Rodriguez-Morales, B.-M. A., Alarcón-Olave, C., and Calvo-Betancourt, L. S. (2015). “Plasmodium vivax Malaria in Latin America. Neglected Tropical Diseases - Latin America and the Caribbean,” in Neglected Tropical Diseases, eds C. Franco-Paredes and J. Santos-Preciado (Vienna: Springer), 89–111.
Roobsoong, W., Roytrakul, S., Sattabongkot, J., Li, J., Udomsangpetch, R., and Cui, L. (2011). Determination of the Plasmodium vivax schizont stage proteome. J. Proteomics 74, 1701–1710. doi: 10.1016/j.jprot.2011.03.035
Russell, B., Chalfein, F., Prasetyorini, B., Kenangalem, E., Piera, K., Suwanarusk, R., et al. (2008). Determinants of in vitro drug susceptibility testing of Plasmodium vivax. Antimicrob. Agents Chemother. 52, 1040–1045. doi: 10.1128/AAC.01334-07
Sattabongkot, J., Yimamnuaychoke, N., Leelaudomlipi, S., Rasameesoraj, M., Jenwithisuk, R., Coleman, R. E., et al. (2006). Establishment of a human hepatocyte line that supports in vitro development of the exo-erythrocytic stages of the malaria parasites Plasmodium falciparum and P. vivax. Am. J. Trop. Med. Hyg. 74, 708–715. doi: 10.4269/ajtmh.2006.74.708
Severini, C., Menegon, M., Di Luca, M., Abdullaev, I., Majori, G., Razakov, S. A., et al. (2004). Risk of Plasmodium vivax malaria reintroduction in Uzbekistan: genetic characterization of parasites and status of potential malaria vectors in the Surkhandarya region. Trans. R. Soc. Trop. Med. Hyg. 98, 585–592. doi: 10.1016/j.trstmh.2004.01.003
Siegel, T. N., Hon, C. C., Zhang, Q., Lopez-Rubio, J. J., Scheidig-Benatar, C., Martins, R. M., et al. (2014). Strand-specific RNA-Seq reveals widespread and developmentally regulated transcription of natural antisense transcripts in Plasmodium falciparum. BMC Genomics 15:150. doi: 10.1186/1471-2164-15-150
Siqueira, A. M., Lacerda, M. V., Magalhaes, B. M., Mourao, M. P., Melo, G. C., Alexandre, M. A., et al. (2015). Characterization of Plasmodium vivax-associated admissions to reference hospitals in Brazil and India. BMC Med. 13:57. doi: 10.1186/s12916-015-0302-y
Siqueira, A. M., Alexandre, M. A., Mourao, M. P., Santos, V. S., Nagahashi-Marie, S. K., Alecrim, M. G., et al. (2010). Severe rhabdomyolysis caused by Plasmodium vivax malaria in the Brazilian Amazon. Am. J. Trop. Med. Hyg. 83, 271–273. doi: 10.4269/ajtmh.2010.10-0027
Sorber, K., Dimon, M. T., and DeRisi, J. L. (2011). RNA-Seq analysis of splicing in Plasmodium falciparum uncovers new splice junctions, alternative splicing and splicing of antisense transcripts. Nucleic Acids Res. 39, 3820–3835. doi: 10.1093/nar/gkq1223
Suwanarusk, R., Cooke, B. M., Dondorp, A. M., Silamut, K., Sattabongkot, J., White, N. J., et al. (2004). The deformability of red blood cells parasitized by Plasmodium falciparum and P. vivax. J. Infect. Dis. 189, 190–194. doi: 10.1086/380468
Suwanarusk, R., Russell, B., Chavchich, M., Chalfein, F., Kenangalem, E., Kosaisavee, V., et al. (2007). Chloroquine resistant Plasmodium vivax: in vitro characterisation and association with molecular polymorphisms. PLoS ONE 2:e1089. doi: 10.1371/journal.pone.0001089
Tachibana, M., Sato, C., Otsuki, H., Sattabongkot, J., Kaneko, O., Torii, M., et al. (2012). Plasmodium vivax gametocyte protein Pvs230 is a transmission-blocking vaccine candidate. Vaccine 30, 1807–1812. doi: 10.1016/j.vaccine.2012.01.003
Tan, L. K., Yacoub, S., Scott, S., Bhagani, S., and Jacobs, M. (2008). Acute lung injury and other serious complications of Plasmodium vivax malaria. Lancet Infect. Dis. 8, 449–454. doi: 10.1016/S1473-3099(08)70153-1
Tanwar, G. S., Khatri, P. C., Sengar, G. S., Kochar, A., Kochar, S. K., Middha, S., et al. (2011). Clinical profiles of 13 children with Plasmodium vivax cerebral malaria. Ann. Trop. Paediatr. 31, 351–356. doi: 10.1179/1465328111Y.0000000040
Taylor, J. E., Pacheco, M. A., Bacon, D. J., Beg, M. A., Machado, R. L., Fairhurst, R. M., et al. (2013). The evolutionary history of Plasmodium vivax as inferred from mitochondrial genomes: parasite genetic diversity in the Americas. Mol. Biol. Evol. 30, 2050–2064. doi: 10.1093/molbev/mst104
Tjitra, E., Anstey, N. M., Sugiarto, P., Warikar, N., Kenangalem, E., Karyana, M., et al. (2008). Multidrug-resistant Plasmodium vivax associated with severe and fatal malaria: a prospective study in Papua, Indonesia. PLoS Med. 5:e128. doi: 10.1371/journal.pmed.0050128
Watanabe, J., Sasaki, M., Suzuki, Y., and Sugano, S. (2002). Analysis of transcriptomes of human malaria parasite Plasmodium falciparum using full-length enriched library: identification of novel genes and diverse transcription start sites of messenger RNAs. Gene 291, 105–113. doi: 10.1016/S0378-1119(02)00552-8
Westenberger, S. J., McClean, C. M., Chattopadhyay, R., Dharia, N. V., Carlton, J. M., Barnwell, J. W., et al. (2010). A systems-based analysis of Plasmodium vivax lifecycle transcription from human to mosquito. PLoS Negl. Trop. Dis. 4:e653. doi: 10.1371/journal.pntd.0000653
White, M. T., Shirreff, G., Karl, S., Ghani, A. C., and Mueller, I. (2016). Variation in relapse frequency and the transmission potential of Plasmodium vivax malaria. Proc. Biol. Sci. 283:20160048. doi: 10.1098/rspb.2016.0048
Williams, T. N., Maitland, K., Phelps, L., Bennett, S., Peto, T. E., Viji, J., et al. (1997). Plasmodium vivax: a cause of malnutrition in young children. QJM 90, 751–757. doi: 10.1093/qjmed/90.12.751
Winter, D. J., Pacheco, M. A., Vallejo, A. F., Schwartz, R. S., Arevalo-Herrera, M., Herrera, S., et al. (2015). Whole genome sequencing of field isolates reveals extensive genetic diversity in Plasmodium vivax from Colombia. PLoS Negl. Trop. Dis. 9:e0004252. doi: 10.1371/journal.pntd.0004252
Yuda, M., Iwanaga, S., Kaneko, I., and Kato, T. (2015). Global transcriptional repression: an initial and essential step for Plasmodium sexual development. Proc. Natl. Acad. Sci. U.S.A. 112, 12824–12829. doi: 10.1073/pnas.1504389112
Zakeri, S., Raeisi, A., Afsharpad, M., Kakar, Q., Ghasemi, F., Atta, H., et al. (2010). Molecular characterization of Plasmodium vivax clinical isolates in Pakistan and Iran using pvmsp-1, pvmsp-3alpha and pvcsp genes as molecular markers. Parasitol. Int. 59, 15–21. doi: 10.1016/j.parint.2009.06.006
Keywords: malaria, Plasmodium vivax, genome, transcriptome, proteome
Citation: Bourgard C, Albrecht L, Kayano ACAV, Sunnerhagen P and Costa FTM (2018) Plasmodium vivax Biology: Insights Provided by Genomics, Transcriptomics and Proteomics. Front. Cell. Infect. Microbiol. 8:34. doi: 10.3389/fcimb.2018.00034
Received: 29 November 2017; Accepted: 25 January 2018;
Published: 08 February 2018.
Edited by:Nathan W. Schmidt, University of Louisville, United States
Reviewed by:João Marcelo Pereira Alves, University of São Paulo, Brazil
Tracey Lamb, University of Utah, United States
Copyright © 2018 Bourgard, Albrecht, Kayano, Sunnerhagen and Costa. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fabio T. M. Costa, firstname.lastname@example.org