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

A new Plasmodium vivax reference genome for South American isolates

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Plasmodium vivax transcriptomics: What is new?

Malaria is the leading human parasitosis and is transmitted through the bite of anopheline mosquitoes infected with parasites of the genus Plasmodium spp. Among the seven species that cause malaria in humans, Plasmodium vivax is the most prevalent species in Latin America. In recent years, there have been an increasing number of reports of clinical complications caused by P. vivax infections, which were previously neglected and underestimated. P. vivax biology remains with large gaps. The emergence of next-generation sequencing technology has ensured a breakthrough in species knowledge. Coupled with this, the deposition of the P. vivax Sal-1 reference genome allowed an increase in transcriptomics projects by accessing messenger RNA. Thus, the regulation of differential gene expression according to the parasite life stage was verified, and several expressed genes were linked to different biological functions. Today, with the progress associated with RNA sequencing technologies, it is possible to detect nuances and obtain robust results. Discoveries provided by transcriptomic studies allow us to understand topics such as RNA expression and regulation and proteins and metabolic pathways involved during different stages of the parasite life cycle. The information obtained enables a better comprehension of immune system evasion mechanisms; invasion and adhesion strategies used by the parasite; as well as new vaccine targets, potential molecular markers, and others therapeutic targets. In this review, we provide new insights into P. vivax biology by summarizing recent findings in transcriptomic studies.

No Association between the Plasmodium vivax crt-o MS334 or In9pvcrt Polymorphisms and Chloroquine Failure in a Pre-Elimination Clinical Cohort from Malaysia with a Large Clonal Expansion

Increasing reports of resistance to a frontline malaria blood-stage treatment, chloroquine (CQ), raises concerns for the elimination of Plasmodium vivax. The absence of an effective molecular marker of CQ resistance in P. vivax greatly constrains surveillance of this emerging threat. A recent genetic cross between CQ sensitive (CQS) and CQ resistant (CQR) NIH-1993 strains of P. vivax linked a moderate CQR phenotype with two candidate markers in P. vivax CQ resistance transporter gene (pvcrt-o): MS334 and In9pvcrt. Longer TGAAGH motif lengths at MS334 were associated with CQ resistance, as were shorter motifs at the In9pvcrt locus. In this study, high-grade CQR clinical isolates of P. vivax from a low endemic setting in Malaysia were used to investigate the association between the MS334 and In9pvcrt variants and treatment efficacy. Among a total of 49 independent monoclonal P. vivax isolates assessed, high-quality MS334 and In9pvcrt sequences could be derived from 30 (61%) and 23 (47%), respectively. Five MS334 and six In9pvcrt alleles were observed, with allele frequencies ranging from 2 to 76% and 3 to 71%, respectively. None of the clinical isolates had the same variant as the NIH-1993 CQR strain, and none of the variants were associated with CQ treatment failure (all P > 0.05). Multi-locus genotypes (MLGs) at 9 neutral microsatellites revealed a predominant P. vivax strain (MLG6) accounting for 52% of Day 0 infections. The MLG6 strain comprised equal proportions of CQS and CQR infections. Our study reveals complexity in the genetic basis of CQ resistance in the Malaysian P. vivax pre-elimination setting and suggests that the proposed pvcrt-o MS334 and In9pvcrt markers are not reliable markers of CQ treatment efficacy in this setting. Further studies are needed in other endemic settings, applying hypothesis-free genome-wide approaches, and functional approaches to understand the biological impact of the TGAAGH repeats linked to CQ response in a cross are warranted to comprehend and track CQR P. vivax.

Reduced polymorphism of Plasmodium vivax early transcribed membrane protein (PvETRAMP) 11.2

BACKGROUND

ETRAMP11.2 (PVX_003565) is a well-characterized protein with antigenic potential. It is considered to be a serological marker for diagnostic tools, and it has been suggested as a potential vaccine candidate. Despite its immunological relevance, the polymorphism of the P. vivax ETRAMP11.2 gene (pvetramp11.2) remains undefined. The genetic variability of an antigen may limit the effectiveness of its application as a serological surveillance tool and in vaccine development and, therefore, the aim of this study was to investigate the genetic diversity of pvetramp11.2 in parasite populations from Amazonian regions and worldwide. We also evaluated amino acid polymorphism on predicted B-cell epitopes. The low variability of the sequence encoding PvETRAMP11.2 protein suggests that it would be a suitable marker in prospective serodiagnostic assays for surveillance strategies or in vaccine design against P. vivax malaria.

METHODS

The pvetramp11.2 of P. vivax isolates collected from Brazil (n = 68) and Peru (n = 36) were sequenced and analyzed to assess nucleotide polymorphisms, allele distributions, population differentiation, genetic diversity and signature of selection. In addition, sequences (n = 104) of seven populations from different geographical regions were retrieved from the PlasmoDB database and included in the analysis to study the worldwide allele distribution. Potential linear B-cell epitopes and their polymorphisms were also explored.

RESULTS

The multiple alignments of 208 pvetramp11.2 sequences revealed a low polymorphism and a marked geographical variation in allele diversity. Seven polymorphic sites and 11 alleles were identified. All of the alleles were detected in isolates from the Latin American region and five alleles were detected in isolates from the Southeast Asia/Papua New Guinea (SEA/PNG) region. Three alleles were shared by all Latin American populations (H1, H6 and H7). The H1 allele (reference allele from Salvador-1 strain), which was absent in the SEA/PNG populations, was the most represented allele in populations from Brazil (54%) and was also detected at high frequencies in populations from all other Latin America countries (range: 13.0% to 33.3%). The H2 allele was the major allele in SEA/PNG populations, but was poorly represented in Latin America populations (only in Brazil: 7.3%). Plasmodium vivax populations from Latin America showed a marked inter-population genetic differentiation (fixation index [Fst]) in contrast to SEA/PNG populations. Codon bias measures (effective number of codons [ENC] and Codon bias index [CBI]) indicated preferential use of synonymous codons, suggesting selective pressure at the translation level. Only three amino acid substitutions, located in the C-terminus, were detected. Linear B-cell epitope mapping predicted two epitopes in the Sal-1 PvETRAMP11.2 protein, one of which was fully conserved in all of the parasite populations analyzed.

CONCLUSIONS

We provide an overview of the allele distribution and genetic differentiation of ETRAMP11.2 antigen in P. vivax populations from different endemic areas of the world. The reduced polymorphism and the high degree of protein conservation supports the application of PvETRAMP11.2 protein as a reliable antigen for application in serological assays or vaccine design. Our findings provide useful information that can be used to inform future study designs.

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

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

Systems biology of malaria explored with nonhuman primates

"The Primate Malarias" book has been a uniquely important resource for multiple generations of scientists, since its debut in 1971, and remains pertinent to the present day. Indeed, nonhuman primates (NHPs) have been instrumental for major breakthroughs in basic and pre-clinical research on malaria for over 50 years. Research involving NHPs have provided critical insights and data that have been essential for malaria research on many parasite species, drugs, vaccines, pathogenesis, and transmission, leading to improved clinical care and advancing research goals for malaria control, elimination, and eradication. Whilst most malaria scientists over the decades have been studying Plasmodium falciparum, with NHP infections, in clinical studies with humans, or using in vitro culture or rodent model systems, others have been dedicated to advancing research on Plasmodium vivax, as well as on phylogenetically related simian species, including Plasmodium cynomolgi, Plasmodium coatneyi, and Plasmodium knowlesi. In-depth study of these four phylogenetically related species over the years has spawned the design of NHP longitudinal infection strategies for gathering information about ongoing infections, which can be related to human infections. These Plasmodium-NHP infection model systems are reviewed here, with emphasis on modern systems biological approaches to studying longitudinal infections, pathogenesis, immunity, and vaccines. Recent discoveries capitalizing on NHP longitudinal infections include an advanced understanding of chronic infections, relapses, anaemia, and immune memory. With quickly emerging new technological advances, more in-depth research and mechanistic discoveries can be anticipated on these and additional critical topics, including hypnozoite biology, antigenic variation, gametocyte transmission, bone marrow dysfunction, and loss of uninfected RBCs. New strategies and insights published by the Malaria Host-Pathogen Interaction Center (MaHPIC) are recapped here along with a vision that stresses the importance of educating future experts well trained in utilizing NHP infection model systems for the pursuit of innovative, effective interventions against malaria.

Transcriptome assembly and differential gene expression of the invasive avian malaria parasite Plasmodium relictum in Hawai'i

The malaria parasite Plasmodium relictum (lineage GRW4) was introduced less than a century ago to the native avifauna of Hawai'i, where it has since caused major declines of endemic bird populations. One of the native bird species that is frequently infected with GRW4 is the Hawai'i 'amakihi (Chlorodrepanis virens). To achieve a better understanding of the transcriptional activities of this virulent parasite, we performed a controlled challenge experiment of 15 'amakihi that were infected with GRW4. Blood samples containing malaria parasites were collected at two time points (intermediate and peak infection stages) from host individuals that were either experimentally infected by mosquitoes or inoculated with infected blood. We then used RNA sequencing to assemble a high-quality blood transcriptome of P. relictum GRW4, allowing us to quantify parasite expression levels inside individual birds. We found few significant differences (one to two transcripts) in GRW4 expression levels between host infection stages and between inoculation methods. However, 36 transcripts showed differential expression levels among all host individuals, indicating a potential presence of host-specific gene regulation across hosts. To reduce the extinction risk of the remaining native bird species in Hawai'i, genetic resources of the local Plasmodium lineage are needed to enable further molecular characterization of this parasite. Our newly built Hawaiian GRW4 transcriptome assembly, together with analyses of the parasite's transcriptional activities inside the blood of Hawai'i 'amakihi, can provide us with important knowledge on how to combat this deadly avian disease in the future.

Dual RNA Sequencing Meta-analysis in Plasmodium Infection Identifies Host-Parasite Interactions

Dual RNA sequencing (RNA-Seq) is the simultaneous transcriptomic analysis of interacting symbionts, for example, in malaria. Potential cross-species interactions identified by correlated gene expression might highlight interlinked signaling, metabolic, or gene regulatory pathways in addition to physically interacting proteins. Often, malaria studies address one of the interacting organisms-host or parasite-rendering the other "contamination." Here we perform a meta-analysis using such studies for cross-species expression analysis. We screened experiments for gene expression from host and Plasmodium. Out of 171 studies in Homo sapiens, Macaca mulatta, and Mus musculus, we identified 63 potential studies containing host and parasite data. While 16 studies (1,950 samples) explicitly performed dual RNA-Seq, 47 (1,398 samples) originally focused on one organism. We found 915 experimental replicates from 20 blood studies to be suitable for coexpression analysis and used orthologs for meta-analysis across different host-parasite systems. Centrality metrics from the derived gene expression networks correlated with gene essentiality in the parasites. We found indications of host immune response to elements of the Plasmodium protein degradation system, an antimalarial drug target. We identified well-studied immune responses in the host with our coexpression networks, as our approach recovers known broad processes interlinked between hosts and parasites in addition to individual host and parasite protein associations. The set of core interactions represents commonalities between human malaria and its model systems for prioritization in laboratory experiments. Our approach might also allow insights into the transferability of model systems for different pathways in malaria studies.IMPORTANCE Malaria still causes about 400,000 deaths a year and is one of the most studied infectious diseases. The disease is studied in mice and monkeys as lab models to derive potential therapeutic intervention in human malaria. Interactions between Plasmodium spp. and its hosts are either conserved across different host-parasite systems or idiosyncratic to those systems. Here we use correlation of gene expression from different RNA-Seq studies to infer common host-parasite interactions across human, mouse, and monkey studies. First, we find a set of very conserved interactors, worth further scrutiny in focused laboratory experiments. Second, this work might help assess to which extent experiments and knowledge on different pathways can be transferred from models to humans for potential therapy.

A suitable RNA preparation methodology for whole transcriptome shotgun sequencing harvested from Plasmodium vivax-infected patients

Plasmodium vivax is a world-threatening human malaria parasite, whose biology remains elusive. The unavailability of in vitro culture, and the difficulties in getting a high number of pure parasites makes RNA isolation in quantity and quality a challenge. Here, a methodological outline for RNA-seq from P. vivax isolates with low parasitemia is presented, combining parasite maturation and enrichment with efficient RNA extraction, yielding ~ 100 pg.µL⁻¹ of RNA, suitable for SMART-Seq Ultra-Low Input RNA library and Illumina sequencing. Unbiased coding transcriptome of ~ 4 M reads was achieved for four patient isolates with ~ 51% of transcripts mapped to the P. vivax P01 reference genome, presenting heterogeneous profiles of expression among individual isolates. Amongst the most transcribed genes in all isolates, a parasite-staged mixed repertoire of conserved parasite metabolic, membrane and exported proteins was observed. Still, a quarter of transcribed genes remain functionally uncharacterized. In parallel, a P. falciparum Brazilian isolate was also analyzed and 57% of its transcripts mapped against IT genome. Comparison of transcriptomes of the two species revealed a common trophozoite-staged expression profile, with several homologous genes being expressed. Collectively, these results will positively impact vivax research improving knowledge of P. vivax biology.

Alternative Invasion Mechanisms and Host Immune Response to Plasmodium vivax Malaria: Trends and Future Directions

Plasmodium vivax malaria is a neglected tropical disease, despite being more geographically widespread than any other form of malaria. The documentation of P. vivax infections in different parts of Africa where Duffy-negative individuals are predominant suggested that there are alternative pathways for P. vivax to invade human erythrocytes. Duffy-negative individuals may be just as fit as Duffy-positive individuals and are no longer resistant to P.vivax malaria. In this review, we describe the complexity of P. vivax malaria, characterize pathogenesis and candidate invasion genes of P. vivax, and host immune responses to P. vivax infections. We provide a comprehensive review on parasite ligands in several Plasmodium species that further justify candidate genes in P. vivax. We also summarize previous genomic and transcriptomic studies related to the identification of ligand and receptor proteins in P. vivax erythrocyte invasion. Finally, we identify topics that remain unclear and propose future studies that will greatly contribute to our knowledge of P. vivax.

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

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

Plasmodium vivax spleen-dependent genes encode antigens associated with cytoadhesion and clinical protection

In spite of low peripheral blood parasitemia, vivax malaria causes severe disease. This conundrum finds an explanation from reports suggesting that the spleen is a place for parasite sequestration. We performed a global transcriptional analysis of parasites that grew in the presence or absence of the spleen in a nonhuman primate model. We identified 67 spleen-dependent genes, including multigene variant families, and functionally demonstrated specific adherence to human spleen fibroblasts by a member of such families. Moreover, we further demonstrated that spleen-dependent Plasmodium vivax genes code for immunogenic proteins during natural infections. Our results indicate that this organ plays an important function in P. vivax malaria and call for deeper studies of the role of spleen in P. vivax infections.

Plasmodium vivax, the most widely distributed human malaria parasite, causes severe clinical syndromes despite low peripheral blood parasitemia. This conundrum is further complicated as cytoadherence in the microvasculature is still a matter of investigations. Previous reports in Plasmodium knowlesi, another parasite species shown to infect humans, demonstrated that variant genes involved in cytoadherence were dependent on the spleen for their expression. Hence, using a global transcriptional analysis of parasites obtained from spleen-intact and splenectomized monkeys, we identified 67 P. vivax genes whose expression was spleen dependent. To determine their role in cytoadherence, two Plasmodium falciparum transgenic lines expressing two variant proteins pertaining to VIR and Pv-FAM-D multigene families were used. Cytoadherence assays demonstrated specific binding to human spleen but not lung fibroblasts of the transgenic line expressing the VIR14 protein. To gain more insights, we expressed five P. vivax spleen-dependent genes as recombinant proteins, including members of three different multigene families (VIR, Pv-FAM-A, Pv-FAM-D), one membrane transporter (SECY), and one hypothetical protein (HYP1), and determined their immunogenicity and association with clinical protection in a prospective study of 383 children in Papua New Guinea. Results demonstrated that spleen-dependent antigens are immunogenic in natural infections and that antibodies to HYP1 are associated with clinical protection. These results suggest that the spleen plays a major role in expression of parasite proteins involved in cytoadherence and can reveal antigens associated with clinical protection, thus prompting a paradigm shift in P. vivax biology toward deeper studies of the spleen during infections.

Single-cell transcription analysis of Plasmodium vivax blood-stage parasites identifies stage- and species-specific profiles of expression

Plasmodium vivax and P. falciparum, the parasites responsible for most human malaria worldwide, exhibit striking biological differences, which have important clinical consequences. Unfortunately, P. vivax, unlike P. falciparum, cannot be cultivated continuously in vitro, which limits our understanding of its biology and, consequently, our ability to effectively control vivax malaria. Here, we describe single-cell gene expression profiles of 9,215 P. vivax parasites from bloodstream infections of Aotus and Saimiri monkeys. Our results show that transcription of most P. vivax genes occurs during short periods of the intraerythrocytic cycle and that this pattern of gene expression is conserved in other Plasmodium species. However, we also identify a strikingly high proportion of species-specific transcripts in late schizonts, possibly associated with the specificity of erythrocyte invasion. Our findings provide new and robust markers of blood-stage parasites, including some that are specific to the elusive P. vivax male gametocytes, and will be useful for analyzing gene expression data from laboratory and field samples.

Analysis of individual Plasmodium vivax parasites reveals the tight control of the expression of most genes during the intra-erythrocytic cycle and the differentiation of male and female gametocytes, and highlights differences between the development of P. vivax and P. falciparum.

Functional genomics of simian malaria parasites and host-parasite interactions

Two simian malaria parasite species, Plasmodium knowlesi and Plasmodium cynomolgi, cause zoonotic infections in Southeast Asia, and they have therefore gained recognition among scientists and public health officials. Notwithstanding, these species and others including Plasmodium coatneyi have served for decades as sources of knowledge on the biology, genetics and evolution of Plasmodium, and the diverse ramifications and outcomes of malaria in their monkey hosts. Experimental analysis of these species can help to fill gaps in knowledge beyond what may be possible studying the human malaria parasites or rodent parasite species. The genome sequences for these simian malaria parasite species were reported during the last decade, and functional genomics research has since been pursued. Here research on the functional genomics analysis involving these species is summarized and their importance is stressed, particularly for understanding host–parasite interactions, and potentially testing novel interventions. Importantly, while Plasmodium falciparum and Plasmodium vivax can be studied in small New World monkeys, the simian malaria parasites can be studied more effectively in the larger Old World monkey macaque hosts, which are more closely related to humans. In addition to ex vivo analyses, experimental scenarios can include passage through Anopheline mosquito hosts and longitudinal infections in monkeys to study acute and chronic infections, as well as relapses, all in the context of the in vivo host environment. Such experiments provide opportunities for understanding functional genomic elements that govern host–parasite interactions, immunity and pathogenesis in-depth, addressing hypotheses not possible from in vitro cultures or cross-sectional clinical studies with humans.

A glance of the blood stage transcriptome of a Southeast Asian Plasmodium ovale isolate

Plasmodium ovale accounts for a disproportionate number of travel-related malaria cases. This parasite is understudied since there is a reliance on clinical samples. We collected a P. ovale curtisi parasite isolate from a clinical case in western Thailand and performed RNA-seq analysis on the blood stage transcriptomes. Using both de novo assembly and alignment-based methods, we detected the transcripts for 6628 out of 7280 annotated genes. For those lacking evidence of expression, the vast majority belonged to the PIR and STP1 gene families. We identified new splicing patterns for over 2500 genes, and mapped at least one untranslated region for over half of all annotated genes. Our analysis also detected a notable presence of anti-sense transcripts for over 10% of P. ovale curtisi genes. This transcriptomic analysis provides new insights into the blood-stage biology of this neglected parasite.

Plasmodium vivax chloroquine resistance links to pvcrt transcription in a genetic cross

Mainstay treatment for Plasmodium vivax malaria has long relied on chloroquine (CQ) against blood-stage parasites plus primaquine against dormant liver-stage forms (hypnozoites), however drug resistance confronts this regimen and threatens malaria control programs. Understanding the basis of P. vivax chloroquine resistance (CQR) will inform drug discovery and malaria control. Here we investigate the genetics of P. vivax CQR by a cross of parasites differing in drug response. Gametocytogenesis, mosquito infection, and progeny production are performed with mixed parasite populations in nonhuman primates, as methods for P. vivax cloning and in vitro cultivation remain unavailable. Linkage mapping of progeny surviving >15 mg/kg CQ identifies a 76 kb region in chromosome 1 including pvcrt, an ortholog of the Plasmodium falciparum CQR transporter gene. Transcriptional analysis supports upregulated pvcrt expression as a mechanism of CQR.

Clinical expression and antigenic profiles of a Plasmodium vivax vaccine candidate: merozoite surface protein 7 (PvMSP-7)

Background

Vivax malaria is the predominant form of malaria outside Africa, affecting about 14 million people worldwide, with about 2.5 billion people exposed. Development of a Plasmodium vivax vaccine is a priority, and merozoite surface protein 7 (MSP-7) has been proposed as a plausible candidate. The P. vivax genome contains 12 MSP-7 genes, which contribute to erythrocyte invasion during blood-stage infection. Previous analysis of MSP-7 sequence diversity suggested that not all paralogs are functionally equivalent. To explore MSP-7 functional diversity, and to identify the best vaccine candidate within the family, MSP-7 expression and antigenicity during bloodstream infections were examined directly from clinical isolates.

Methods

Merozoite surface protein 7 gene expression was profiled using RNA-seq data from blood samples isolated from ten human patients with vivax malaria. Differential expression analysis and co-expression cluster analysis were used to relate PvMSP-7 expression to genetic markers of life cycle stage. Plasma from vivax malaria patients was also assayed using a custom peptide microarray to measure antibody responses against the coding regions of 12 MSP-7 paralogs.

Results

Ten patients presented diverse transcriptional profiles that comprised four patient groups. Two MSP-7 paralogs, 7A and 7F, were expressed abundantly in all patients, while other MSP-7 genes were uniformly rare (e.g. 7J). MSP-7H and 7I were significantly more abundant in patient group 4 only, (two patients having experienced longer patency), and were co-expressed with a schizont-stage marker, while negatively associated with liver-stage and gametocyte-stage markers. Screening infections with a PvMSP-7 peptide array identified 13 linear B-cell epitopes in five MSP-7 paralogs that were recognized by plasma from all patients.

Conclusions

These results show that MSP-7 family members vary in expression profile during blood infections; MSP-7A and 7F are expressed throughout the intraerythrocytic development cycle, while expression of other paralogs is focused on the schizont. This may reflect developmental regulation, and potentially functional differentiation, within the gene family. The frequency of B-cell epitopes among paralogs also varies, with MSP-7A and 7L consistently the most immunogenic. Thus, MSP-7 paralogs cannot be assumed to have equal potential as vaccines. This analysis of clinical infections indicates that the most abundant and immunogenic paralog is MSP-7A.

Electronic supplementary material

The online version of this article (10.1186/s12936-019-2826-7) contains supplementary material, which is available to authorized users.

Transcriptome and histone epigenome of Plasmodium vivax salivary-gland sporozoites point to tight regulatory control and mechanisms for liver-stage differentiation in relapsing malaria

Plasmodium vivax is the key obstacle to malaria elimination in Asia and Latin America, largely attributed to its ability to form resilient hypnozoites (sleeper cells) in the host liver that escape treatment and cause relapsing infections. The decision to form hypnozoites is made early in the liver infection and may already be set in sporozoites prior to invasion. To better understand these early stages of infection, we undertook a comprehensive transcriptomic and histone epigenetic characterization of P. vivax sporozoites. Through comparisons with recently published proteomic data for the P. vivax sporozoite, our study found that although highly transcribed, transcripts associated with functions needed for early infection of the vertebrate host are not detectable as proteins and may be regulated through translational repression. We identified differential transcription between the sporozoite and published transcriptomes of asexual blood stages and mixed versus hypnozoite-enriched liver stages. These comparisons point to multiple layers of transcriptional, post-transcriptional and post-translational control that appear active in sporozoites and to a lesser extent hypnozoites, but are largely absent in replicating liver schizonts or mixed blood stages. We also characterised histone epigenetic modifications in the P. vivax sporozoite and explored their role in regulating transcription. Collectively, these data support the hypothesis that the sporozoite is a tightly programmed stage to infect the human host and identify mechanisms for hypnozoite formation that may be further explored in liver stage models.

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

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

Unraveling the Plasmodium vivax sporozoite transcriptional journey from mosquito vector to human host

Malaria parasites transmitted by mosquito bite are remarkably efficient in establishing human infections. The infection process requires roughly 30 minutes and is highly complex as quiescent sporozoites injected with mosquito saliva must be rapidly activated in the skin, migrate through the body, and infect the liver. This process is poorly understood for Plasmodium vivax due to low infectivity in the in vitro models. To study this skin-to-liver-stage of malaria, we used quantitative bioassays coupled with transcriptomics to evaluate parasite changes linked with mammalian microenvironmental factors. Our in vitro phenotyping and RNA-seq analyses revealed key microenvironmental relationships with distinct biological functions. Most notable, preservation of sporozoite quiescence by exposure to insect-like factors coupled with strategic activation limits untimely activation of invasion-associated genes to dramatically increase hepatocyte invasion rates. We also report the first transcriptomic analysis of the P. vivax sporozoite interaction in salivary glands identifying 118 infection-related differentially-regulated Anopheles dirus genes. These results provide important new insights in malaria parasite biology and identify priority targets for antimalarial therapeutic interventions to block P. vivax infection.

Systems Biology-Based Investigation of Host-Plasmodium Interactions

Malaria is a serious, complex disease caused by parasites of the genus Plasmodium. Plasmodium parasites affect multiple tissues as they evade immune responses, replicate, sexually reproduce, and transmit between vertebrate and invertebrate hosts. The explosion of omics technologies has enabled large-scale collection of Plasmodium infection data, revealing systems-scale patterns, mechanisms of pathogenesis, and the ways that host and pathogen affect each other. Here, we provide an overview of recent efforts using systems biology approaches to study host-Plasmodium interactions and the biological themes that have emerged from these efforts. We discuss some of the challenges in using systems biology for this goal, key research efforts needed to address those issues, and promising future malaria applications of systems biology.

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

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