The transcriptome of Plasmodium vivax reveals divergence and diversity of transcriptional regulation in malaria parasites

The biology and pathogenesis of vivax malaria

Plasmodium vivax contributes significantly to global malaria morbidity. Key advances include the discovery of pathways facilitating invasion by P. vivax merozoites of nascent reticulocytes, crucial for vaccine development. Humanized mouse models and hepatocyte culture systems have enhanced understanding of hypnozoite biology. The spleen has emerged as a major reservoir for asexual vivax parasites, replicating in an endosplenic life cycle, and contributing to recurrent and chronic infections, systemic inflammation, and anemia. Splenic accumulation of uninfected red cells is the predominant cause of anemia. Recurring and chronic infections cause progressive anemia, malnutrition, and death in young children in high-transmission regions. Endothelial activation likely contributes to vivax-associated organ dysfunction. The many recent advances in vivax pathobiology should help guide new approaches to prevention and management.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Sterile protection against P. vivax malaria by repeated blood stage infection in the Aotus monkey model

The malaria parasite Plasmodium vivax remains a major global public health challenge, and no vaccine is approved for use in humans. Here, we assessed whether P. vivax strain-transcendent immunity can be achieved by repeated infection in Aotus monkeys. Sterile immunity was achieved after two homologous infections, whereas subsequent heterologous challenge provided only partial protection. IgG levels based on P. vivax lysate ELISA and protein microarray increased with repeated infections and correlated with the level of homologous protection. Parasite transcriptional profiles provided no evidence of major antigenic switching upon homologous or heterologous challenge. However, we observed significant sequence diversity and transcriptional differences in the P. vivax core gene repertoire between the two strains used in the study, suggesting that partial protection upon heterologous challenge is due to molecular differences between strains rather than immune evasion by antigenic switching. Our study demonstrates that sterile immunity against P. vivax can be achieved by repeated homologous blood stage infection in Aotus monkeys, thus providing a benchmark to test the efficacy of candidate blood stage P. vivax malaria vaccines.

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

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

Functional characterization of Plasmodium vivax hexose transporter 1

Plasmodium vivax is the most widely distributed human malaria parasite. The eradication of vivax malaria remains challenging due to transmission of drug-resistant parasite and dormant liver form. Consequently, anti-malarial drugs with novel mechanisms of action are urgently demanded. Glucose uptake blocking strategy is suggested as a novel mode of action that leads to selective starvation in various species of malaria parasites. The role of hexose transporter 1 in Plasmodium species is glucose uptake, and its blocking strategies proved to successfully induce selective starvation. However, there is limited information on the glucose uptake properties via P. vivax hexose transporter 1 (PvHT1). Thus, we focused on the PvHT1 to precisely identify its properties of glucose uptake. The PvHT1 North Korean strain (PvHT1NK) expressed Xenopus laevis oocytes mediating the transport of [3H] deoxy-D-glucose (ddGlu) in an expression and incubation time-dependent manner without sodium dependency. Moreover, the PvHT1NK showed no exchange mode of glucose in efflux experiments and concentration-dependent results showed saturable kinetics following the Michaelis-Menten equation. Non-linear regression analysis revealed a Km value of 294.1 μM and a Vmax value of 1,060 pmol/oocyte/hr, and inhibition experiments showed a strong inhibitory effect by glucose, mannose, and ddGlu. Additionally, weak inhibition was observed with fructose and galactose. Comparison of amino acid sequence and tertiary structure between P. falciparum and P. vivax HT1 revealed a completely conserved residue in glucose binding pocket. This result supported that the glucose uptake properties are similar to P. falciparum, and PfHT1 inhibitor (compound 3361) works in P. vivax. These findings provide properties of glucose uptake via PvHT1NK for carbohydrate metabolism and support the approaches to vivax malaria drug development strategy targeting the PvHT1 for starving of the parasite.

Host kinase regulation of Plasmodium vivax dormant and replicating liver stages

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

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.

The parasite intraerythrocytic cycle and human circadian cycle are coupled during malaria infection

During infections with the malaria parasites Plasmodium vivax, patients exhibit rhythmic fevers every 48 h. These fever cycles correspond with the time the parasites take to traverse the intraerythrocytic cycle (IEC). In other Plasmodium species that infect either humans or mice, the IEC is likely guided by a parasite-intrinsic clock [Rijo-Ferreiraet al., Science 368, 746-753 (2020); Smith et al., Science 368, 754-759 (2020)], suggesting that intrinsic clock mechanisms may be a fundamental feature of malaria parasites. Moreover, because Plasmodium cycle times are multiples of 24 h, the IECs may be coordinated with the host circadian clock(s). Such coordination could explain the synchronization of the parasite population in the host and enable alignment of IEC and circadian cycle phases. We utilized an ex vivo culture of whole blood from patients infected with P. vivax to examine the dynamics of the host circadian transcriptome and the parasite IEC transcriptome. Transcriptome dynamics revealed that the phases of the host circadian cycle and the parasite IEC are correlated across multiple patients, showing that the cycles are phase coupled. In mouse model systems, host-parasite cycle coupling appears to provide a selective advantage for the parasite. Thus, understanding how host and parasite cycles are coupled in humans could enable antimalarial therapies that disrupt this coupling.

Using Serological Markers for the Surveillance of Plasmodium vivax Malaria: A Scoping Review

The utilisation of serological surveillance methods for malaria has the potential to identify individuals exposed to Plasmodium vivax, including asymptomatic carriers. However, the application of serosurveillance varies globally, including variations in methodology and transmission context. No systematic review exists describing the advantages and disadvantages of utilising serosurveillance in various settings. Collation and comparison of these results is a necessary first step to standardise and validate the use of serology for the surveillance of P. vivax in specific transmission contexts. A scoping review was performed of P. vivax serosurveillance applications globally. Ninety-four studies were found that met predefined inclusion and exclusion criteria. These studies were examined to determine the advantages and disadvantages of serosurveillance experienced in each study. If studies reported seroprevalence results, this information was also captured. Measurement of antibodies serves as a proxy by which individuals exposed to P. vivax may be indirectly identified, including those with asymptomatic infections, which may be missed by other technologies. Other thematic advantages identified included the ease and simplicity of serological assays compared to both microscopy and molecular diagnostics. Seroprevalence rates varied widely from 0-93%. Methodologies must be validated across various transmission contexts to ensure the applicability and comparability of results. Other thematic disadvantages identified included challenges with species cross-reactivity and determining changes in transmission patterns in both the short- and long-term. Serosurveillance requires further refinement to be fully realised as an actionable tool. Some work has begun in this area, but more is required.

Identification of novel Plasmodium vivax proteins associated with protection against clinical malaria

As progress towards malaria elimination continues, the challenge posed by the parasite species Plasmodium vivax has become more evident. In many regions co-endemic for P. vivax and Plasmodium falciparum, as transmission has declined the proportion of cases due to P. vivax has increased. Novel tools that directly target P. vivax are thus warranted for accelerated elimination. There is currently no advanced vaccine for P. vivax and only a limited number of potential candidates in the pipeline. In this study we aimed to identify promising P. vivax proteins that could be used as part of a subunit vaccination approach. We screened 342 P. vivax protein constructs for their ability to induce IgG antibody responses associated with protection from clinical disease in a cohort of children from Papua New Guinea. This approach has previously been used to successfully identify novel candidates. We were able to confirm previous results from our laboratory identifying the proteins reticulocyte binding protein 2b and StAR-related lipid transfer protein, as well as at least four novel candidates with similar levels of predicted protective efficacy. Assessment of these P. vivax proteins in further studies to confirm their potential and identify functional mechanisms of protection against clinical disease are warranted.

Short tandem repeat polymorphism in the promoter region of cyclophilin 19B drives its transcriptional upregulation and contributes to drug resistance in the malaria parasite Plasmodium falciparum

Resistance of the human malaria parasites, Plasmodium falciparum, to artemisinins is now fully established in Southeast Asia and is gradually emerging in Sub-Saharan Africa. Although nonsynonymous SNPs in the pfk13 Kelch-repeat propeller (KREP) domain are clearly associated with artemisinin resistance, their functional relevance requires cooperation with other genetic factors/alterations of the P. falciparum genome, collectively referred to as genetic background. Here we provide experimental evidence that P. falciparum cyclophilin 19B (PfCYP19B) may represent one putative factor in this genetic background, contributing to artemisinin resistance via its increased expression. We show that overexpression of PfCYP19B in vitro drives limited but significant resistance to not only artemisinin but also piperaquine, an important partner drug in artemisinin-based combination therapies. We showed that PfCYP19B acts as a negative regulator of the integrated stress response (ISR) pathway by modulating levels of phosphorylated eIF2α (eIF2α-P). Curiously, artemisinin and piperaquine affect eIF2α-P in an inverse direction that in both cases can be modulated by PfCYP19B towards resistance. Here we also provide evidence that the upregulation of PfCYP19B in the drug-resistant parasites appears to be maintained by a short tandem repeat (SRT) sequence polymorphism in the gene's promoter region. These results support a model that artemisinin (and other drugs) resistance mechanisms are complex genetic traits being contributed to by altered expression of multiple genes driven by genetic polymorphism at their promoter regions.

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.

Band 3–mediated Plasmodium vivax invasion is associated with transcriptional variation in PvTRAg genes

The Plasmodium vivax reticulocyte invasion process is still poorly understood, with only a few receptor-ligand interactions identified to date. Individuals with the Southeast Asian ovalocytosis (SAO) phenotype have a deletion in the band 3 protein on the surface of erythrocytes, and are reported to have a lower incidence of clinical P. vivax malaria. Based on this observation, band 3 has been put forward as a receptor for P. vivax invasion, although direct proof is still lacking. In this study, we combined functional ex vivo invasion assays and transcriptome sequencing to uncover a band 3–mediated invasion pathway in P. vivax and potential band 3 ligands. Invasion by P. vivax field isolates was 67%-71% lower in SAO reticulocytes compared with non-SAO reticulocytes. Reticulocyte invasion was decreased by 40% and 27%-31% when blocking with an anti-band 3 polyclonal antibody and a PvTRAg38 peptide, respectively. To identify new band 3 receptor candidates, we mRNA-sequenced schizont-stage isolates used in the invasion assays, and observed high transcriptional variability in multigene and invasion-related families. Transcriptomes of isolates with low or high dependency on band 3 for invasion were compared by differential expression analysis, which produced a list of band 3 ligand candidates with high representation of PvTRAg genes. Our ex vivo invasion assays have demonstrated that band 3 is a P. vivax invasion receptor and confirm previous in vitro studies showing binding between PvTRAg38 and band 3, although the lower and variable inhibition levels observed suggest the involvement of other ligands. By coupling transcriptomes and invasion phenotypes from the same isolates, we identified a list of band 3 ligand candidates, of which the overrepresented PvTRAg genes are the most promising for future research.

Determination of the Stage Composition of Plasmodium Infections from Bulk Gene Expression Data

Malaria symptoms are caused by the development of the parasites within the blood of an infected host. Bulk RNA sequencing (RNA-seq) of infected blood can reveal interactions between parasites and the host immune system during an infection, but because multiple developmental stages with distinct transcriptional profiles are concurrently present in infected blood, it is necessary to correct such analyses for differences in cell composition among samples. Gene expression deconvolution is a statistical approach that has been developed for inferring the cell composition of complex tissues characterized by bulk RNA-seq using gene expression profiles from reference cell types. Here, we describe the evaluation of a species-agnostic reference data set that can be used for efficient and accurate gene expression deconvolution of bulk RNA-seq data generated from any Plasmodium species and for correct gene expression analyses for biases caused by differences in stage composition among samples. IMPORTANCE Differences in cell type proportions among samples can introduce artifacts in gene expression analyses and mask genuine differences in gene regulation. Gene expression deconvolution allows estimation of the proportion of each cell type present in one sample directly from bulk RNA sequencing data, but this approach requires a reference data set with the signature profile of each cell type. Here, we evaluate the suitability of a rodent malaria parasite gene expression data set for estimating the proportions of each parasite developmental stage present in bulk RNA sequencing data generated from blood-stage infections with the human parasites Plasmodium falciparum and Plasmodium vivax. These analyses provide a species-agnostic approach for reliably estimating stage proportions in infected human blood and correcting subsequent gene expression analyses for these variations.

Single-cell RNA sequencing of Plasmodium vivax sporozoites reveals stage- and species-specific transcriptomic signatures

Background

Plasmodium vivax sporozoites reside in the salivary glands of a mosquito before infecting a human host and causing malaria. Previous transcriptome-wide studies in populations of these parasite forms were limited in their ability to elucidate cell-to-cell variation, thereby masking cellular states potentially important in understanding malaria transmission outcomes.

Methodology/Principal findings

In this study, we performed transcription profiling on 9,947 P. vivax sporozoites to assess the extent to which they differ at single-cell resolution. We show that sporozoites residing in the mosquito’s salivary glands exist in distinct developmental states, as defined by their transcriptomic signatures. Additionally, relative to P. falciparum, P. vivax displays overlapping and unique gene usage patterns, highlighting conserved and species-specific gene programs. Notably, distinguishing P. vivax from P. falciparum were a subset of P. vivax sporozoites expressing genes associated with translational regulation and repression. Finally, our comparison of single-cell transcriptomic data from P. vivax sporozoite and erythrocytic forms reveals gene usage patterns unique to sporozoites.

Conclusions/Significance

In defining the transcriptomic signatures of individual P. vivax sporozoites, our work provides new insights into the factors driving their developmental trajectory and lays the groundwork for a more comprehensive P. vivax cell atlas.

Plasmodium vivax is the second most common cause of malaria worldwide. It is particularly challenging for malaria elimination as it forms both active blood-stage infections, as well as asymptomatic liver-stage infections that can persist for extended periods of time. The activation of persister forms in the liver (hypnozoites) are responsible for relapsing infections occurring weeks or months following primary infection via a mosquito bite. How P. vivax persists in the liver remains a major gap in understanding of this organism. It has been hypothesized that there is pre-programming of the infectious sporozoite while it is in the salivary-glands that determines if the cell’s fate once in the liver is to progress towards immediate liver stage development or persist for long-periods as a hypnozoite. The aim of this study was to see if such differences were distinguishable at the transcript level in salivary-gland sporozoites. While we found significant variation amongst sporozoites, we did not find clear evidence that they are transcriptionally pre-programmed as has been suggested. Nevertheless, we highlight several intriguing patterns that appear to be P. vivax specific relative to non-relapsing species that cause malaria prompting further investigation.

Differential Trafficking and Expression of PIR Proteins in Acute and Chronic Plasmodium Infections

Plasmodium multigene families are thought to play important roles in the pathogenesis of malaria. Plasmodium interspersed repeat (pir) genes comprise the largest multigene family in many Plasmodium species. However, their expression pattern and localisation remain to be elucidated. Understanding protein subcellular localisation is fundamental to reveal the functional importance and cell-cell interactions of the PIR proteins. Here, we use the rodent malaria parasite, Plasmodium chabaudi chabaudi, as a model to investigate the localisation pattern of this gene family. We found that most PIR proteins are co-expressed in clusters during acute and chronic infection; members of the S7 clade are predominantly expressed during the acute-phase, whereas members of the L1 clade dominate the chronic-phase of infection. Using peptide antisera specific for S7 or L1 PIRS, we show that these PIRs have different localisations within the infected red blood cells. S7 PIRs are exported into the infected red blood cell cytoplasm where they are co-localised with parasite-induced host cell modifications termed Maurer’s clefts, whereas L1 PIRs are localised on or close to the parasitophorous vacuolar membrane. This localisation pattern changes following mosquito transmission and during progression from acute- to chronic-phase of infection. The presence of PIRs in Maurer’s clefts, as seen for Plasmodium falciparum RIFIN and STEVOR proteins, might suggest trafficking of the PIRs on the surface of the infected erythrocytes. However, neither S7 nor L1 PIR proteins detected by the peptide antisera are localised on the surface of infected red blood cells, suggesting that they are unlikely to be targets of surface variant-specific antibodies or to be directly involved in adhesion of infected red blood cells to host cells, as described for Plasmodium falciparum VAR proteins. The differences in subcellular localisation of the two major clades of Plasmodium chabaudi PIRs across the blood cycle, and the apparent lack of expression on the red cell surface strongly suggest that the function(s) of this gene family may differ from those of other multigene families of Plasmodium, such as the var genes of Plasmodium falciparum.

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.

Artemisinin resistance in the malaria parasite, Plasmodium falciparum, originates from its initial transcriptional response

The emergence and spread of artemisinin-resistant Plasmodium falciparum, first in the Greater Mekong Subregion (GMS), and now in East Africa, is a major threat to global malaria elimination ambitions. To investigate the artemisinin resistance mechanism, transcriptome analysis was conducted of 577 P. falciparum isolates collected in the GMS between 2016–2018. A specific artemisinin resistance-associated transcriptional profile was identified that involves a broad but discrete set of biological functions related to proteotoxic stress, host cytoplasm remodelling, and REDOX metabolism. The artemisinin resistance-associated transcriptional profile evolved from initial transcriptional responses of susceptible parasites to artemisinin. The genetic basis for this adapted response is likely to be complex.

Transcriptomic analysis of isolates from the malaria parasite (Plasmodium falciparum) in the Greater Mekong Subregion of Southeast Asia identifies gene expression patterns that are correlated with resistance to a common anti-malaria drug, artemisinin.

Contribution of Transcriptome to Elucidate the Biology of Plasmodium spp

In the present review, we discuss some of the new technologies that have been applied to elucidate how Plasmodium spp escape from the immune system and subvert the host physiology to orchestrate the regulation of its biological pathways. Our manuscript describes how techniques such as microarray approaches, RNA-Seq and single-cell RNA sequencing have contributed to the discovery of transcripts and changed the concept of gene expression regulation in closely related malaria parasite species. Moreover, the text highlights the contributions of high-throughput RNA sequencing for the current knowledge of malaria parasite biology, physiology, vaccine target and the revelation of new players in parasite signaling.

Decrypting the complexity of the human malaria parasite biology through systems biology approaches

The human malaria parasite, Plasmodium falciparum, is a unicellular protozoan responsible for over half a million deaths annually. With a complex life cycle alternating between human and invertebrate hosts, this apicomplexan is notoriously adept at evading host immune responses and developing resistance to all clinically administered treatments. Advances in omics-based technologies, increased sensitivity of sequencing platforms and enhanced CRISPR based gene editing tools, have given researchers access to more in-depth and untapped information about this enigmatic micro-organism, a feat thought to be infeasible in the past decade. Here we discuss some of the most important scientific achievements made over the past few years with a focus on novel technologies and platforms that set the stage for subsequent discoveries. We also describe some of the systems-based methods applied to uncover gaps of knowledge left through single-omics applications with the hope that we will soon be able to overcome the spread of this life-threatening disease.

Identification of the Recombinant Plasmodium vivax Surface-Related Antigen as a Possible Immune Evasion Factor Against Human Splenic Fibroblasts by Targeting ITGB1

Plasmodium vivax–infected erythrocytes can enter the spleen and evade spleen clearance to establish chronic infections. However, the mechanism underlying P. vivax immune evasion in the spleen is still unclear. Human splenic fibroblasts (HSF), also known as barrier cells, play an essential role in the immune function of spleen. A hypothesis holds that P. vivax—infected erythrocytes induce spleen structural remodeling to form barrier cells. Subsequently, these infected erythrocytes can selectively cytoadhere to these barrier cells to escape spleen clearance. In this work, we found that P. vivax surface-related antigen (PvSRA; PlasmoDB ID: PVX_084970), an exported protein on infected erythrocyte membrane, could bind with HSF. Considering the above hypothesis, we speculated that PvSRA might be involved in P. vivax immune evasion by changing HSF cell performance. To investigate this speculation, RNA sequencing, protein microarray, and bioinformatics analysis technologies were applied, and in vitro validations were further performed. The results showed that the recombinant PvSRA attracted HSF migration and interacted with HSF by targeting integrin β1 (ITGB1) along with changes in HSF cell performance, such as focal adhesion, extracellular matrix, actin cytoskeleton, and cell cycle. This study indicated that PvSRA might indeed participate in the immune evasion of P. vivax in the spleen by changing HSF function through PvSRA–ITGB1 axis.

Plasmodium vivax gametocytes and transmission

Malaria elimination means cessation of parasite transmission. At present, the declining malaria incidence in many countries has made elimination a feasible goal. Transmission control has thus been placed at the center of the national malaria control programs. The efficient transmission of Plasmodium vivax from humans to mosquitoes is a key factor that helps perpetuate malaria in endemic areas. A better understanding of transmission is crucial to the success of elimination efforts. Biological delineation of the parasite transmission process is important for identifying and prioritizing new targets of intervention. Identification of the infectious parasite reservoir in the community is key to devising an effective elimination strategy. Here we describe the fundamental characteristics of P. vivax gametocytes - the dynamics of their production, longevity, and the relationship with the total parasitemia - as well as recent advances in the molecular understanding of parasite sexual development. In relation to malaria elimination, factors influencing the human infectivity and the current evidence for a role of asymptomatic carriers in transmission are presented.

Serology for Plasmodium vivax surveillance: A novel approach to accelerate towards elimination

Plasmodium vivax is the most widespread causative agent of human malaria in the world. Despite the ongoing implementation of malaria control programs, the rate of case reduction has declined over the last 5 years. Hence, surveillance of malaria transmission should be in place to identify and monitor areas that require intensified malaria control interventions. Serological tools may offer additional insights into transmission intensity over parasite and entomological measures, especially as transmission levels decline. Antibodies can be detected in the host system for months to even years after parasite infections have been cleared from the blood, enabling malaria exposure history to be captured. Because the Plasmodium parasite expresses more than 5000 proteins, it is important to a) understand antibody longevity following infection and b) measure antibodies to more than one antigen in order to accurately inform on the exposure and/or immune status of populations. This review summarises current practices for surveillance of P. vivax malaria, the current state of research into serological exposure markers and their potential role for accelerating malaria elimination, and discusses further studies that need to be undertaken to see such technology implemented in malaria-endemic areas.

Variation in selective constraints along the Plasmodium life cycle

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

Shifts in gene expression variability in the blood-stage of Plasmodium relictum

Avian malaria is a common and widespread disease of birds caused by a diverse group of pathogens of the genera Plasmodium. We investigated the transcriptomal profiles of one of the most common species, Plasmodium relictum, lineage SGS1, at multiple timepoints during the blood stages of the infection under experimental settings. The parasite showed well separated overall transcriptome profiles between day 8 and 20 after the infection, shown by well separated PCA profiles. Moreover, gene expression becomes more heterogenous within the experimental group late in the infection, either due to adaptations to individual differences between the experimental hosts, or due to desynchronisation of the life-cycle of the parasite. Overall, this study shows how the avian malaria system can be used to study gene expression of the avian Plasmodium parasite under controlled experimental settings, thus allowing for future comparative analysis of gene responses of parasite with different life-history traits and host effects.

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.

Two 20-Residue-Long Peptides Derived from Plasmodium vivax Merozoite Surface Protein 10 EGF-Like Domains Are Involved in Binding to Human Reticulocytes

Plasmodium parasites’ invasion of their target cells is a complex, multi-step process involving many protein-protein interactions. Little is known about how complex the interaction with target cells is in Plasmodium vivax and few surface molecules related to reticulocytes’ adhesion have been described to date. Natural selection, functional and structural analysis were carried out on the previously described vaccine candidate P. vivax merozoite surface protein 10 (PvMSP10) for evaluating its role during initial contact with target cells. It has been shown here that the recombinant carboxyl terminal region (rPvMSP10-C) bound to adult human reticulocytes but not to normocytes, as validated by two different protein-cell interaction assays. Particularly interesting was the fact that two 20-residue-long regions (³⁸⁸DKEECRCRANYMPDDSVDYF⁴⁰⁷ and ⁴¹⁵KDCSKENGNCDVNAECSIDK⁴³⁴) were able to inhibit rPvMSP10-C binding to reticulocytes and rosette formation using enriched target cells. These peptides were derived from PvMSP10 epidermal growth factor (EGF)-like domains (precisely, from a well-defined electrostatic zone) and consisted of regions having the potential of being B- or T-cell epitopes. These findings provide evidence, for the first time, about the fragments governing PvMSP10 binding to its target cells, thus highlighting the importance of studying them for inclusion in a P. vivax antimalarial vaccine.

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.

Putative circumsporozoite protein (CSP) of Plasmodium vivax is considerably distinct from the well-known CSP and plays a role in the protein ubiquitination pathway

Amidst technical challenges which limit successful culture and genetic manipulation of P. vivax parasites, we used a computational approach to identify a critical target with evolutionary significance. The putative circumsporozoite protein on chromosome 13 of P. vivax (PvpuCSP)is distinct from the well-known vaccine candidate PfCSP. The aim of this study was to understand the role of PvpuCSP and its relatedness to the well-known CSP. The study revealed PvpuCSP as a membrane bound E3 ubiquitin ligase involved in ubiquitination. It has a species-specific tetra-peptide unit which is differentially repeated in various P. vivax strains. The PvpuCSP is different from CSP in terms of stage-specific expression and function. Since E3 ubiquitin ligases are known antimalarial drug targets targeting the proteasome pathway, PvpuCSP, with evolutionary connotation and a key role in orchestrating protein degradation in P. vivax, can be explored as a potential drug target.

•

PvpuCSP is predicted as E3 ubiquitin ligase, a part of ubiquitination pathway.

•

Tetra-peptide tandem repeat at C terminal of PvpuCSP is exclusive to P. vivax.

•

Moderately expressed during all parasitic stages in host and vector

•

Partially disordered protein with both structured domains and two distinct IDRs

•

A transmembrane protein with highly conserved functional domain across Apicomplexa

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.

Inferring Plasmodium vivax protein biology by using omics data

Deciphering Plasmodium vivax biology has long been a challenge for groups working on this parasite, mainly due to the complications involved in propagating it in vitro. However, adapting P. vivax strains in non-human primates and the arrival of high-performance analysis methods has led to increased knowledge regarding parasite protein composition and the ability of some molecules to trigger an immune response or participate in protein-protein interactions. This review describes the state of the art concerning proteomics-, immunomics- and interatomics-related P. vivax omic studies, discussing their potential use in developing disease control methods.

The Blood Stage Antigen RBP2-P1 of Plasmodium vivax Binds Reticulocytes and Is a Target of Naturally Acquired Immunity

The interactions between Plasmodium parasites and human erythrocytes are prime targets of blood stage malaria vaccine development. The reticulocyte binding protein 2-P1 (RBP2-P1) of Plasmodium vivax, a member of the reticulocyte binding protein family, has recently been shown to be highly antigenic in several settings endemic for malaria. Yet, its functional characteristics and the relevance of its antibody response in human malaria have not been examined. In this study, the potential function of RBP2-P1 as an invasion ligand of P. vivax was evaluated. The protein was found to be expressed in schizonts, be localized at the apical end of the merozoite, and preferentially bind reticulocytes over normocytes. Human antibodies to this protein also exhibit erythrocyte binding inhibition at physiologically relevant concentrations. Furthermore, RBP2-P1 antibodies are associated with lower parasitemia and tend to be higher in asymptomatic carriers than in patients. This study provides evidence supporting a role of RBP2-P1 as an invasion ligand and its consideration as a vaccine target.

Plasmodium vivax transcriptional profiling of low input cryopreserved isolates through the intraerythrocytic development cycle

Approximately one-third of the global population is at risk of Plasmodium vivax infection, and an estimated 7.51 million cases were reported in 2017. Although, P. vivax research is currently limited by the lack of a robust continuous in vitro culture system for this parasite, recent work optimizing short-term ex vivo culture of P. vivax from cryopreserved isolates has facilitated quantitative assays on synchronous parasites. Pairing this improved culture system with low-input Smart-seq2 RNAseq library preparation, we sought to determine whether transcriptional profiling of P. vivax would provide insight into the differential survival of parasites in different culture media. To this end we probed the transcriptional signature of three different ex vivo P. vivax samples in four different culture media using only 1000 cells for each time point taken during the course of the intraerythrocytic development cycle (IDC). Using this strategy, we achieved similar quality transcriptional data to previously reported P. vivax transcriptomes. We found little effect with varying culture media on parasite transcriptional signatures, identified many novel gametocyte-specific genes from transcriptomes of FACS-isolated gametocytes, and determined invasion ligand expression in schizonts in biological isolates and across the IDC. In total, these data demonstrate the feasibility and utility of P. vivax RNAseq-based transcriptomic studies using minimal biomass input to maximize experimental capacity.

Plasmodium vivax is the most prevalent malaria-causing parasite species outside of Sub-Saharan Africa and has many unique and poorly understood biological characteristics that make it particularly challenging to study and combat. Transcriptomic profiling of P. vivax under various conditions has the potential to unlock new experimental abilities and aid in elucidating biology and the development of clinical interventions. However, a lack of a robust in vitro culture system for this parasite has restricted transcriptomic studies to researchers with timely access to fresh human isolates from clinics, which often are in resource-poor settings, as well as nearby, well-equipped laboratories for sample processing. This study aimed to gain insight into the differential survival of P. vivax in various culture media from the parasites transcriptional signature in each media. By implementing low-input RNA library preparation strategies, this study obtains robust transcriptomic data at various parasite development stages and in different culture conditions from just 1000 FACS-purified, P. vivax-infected erythrocytes from viable cryopreserved patient isolates. With these data, we find culture media has little effect on transcriptional profile, we characterize invasion ligand expression across intraerythrocytic development and between clinical isolates, and we define the transcriptome of sexual, transmissible stages of the P. vivax parasite. These results highlight the establishment and utility of a powerful platform for studying the transcriptomic biology of this particularly challenging parasite.

From a basic to a functional approach for developing a blood stage vaccine against Plasmodium vivax

Introduction: Numerous challenges have hampered developing an anti-malarial vaccine against the most widespread malarial parasite worldwide: Plasmodium vivax. Despite the progress achieved in studying proteins in short-term in vitro culture or in experimental models, there is still no clear method for defining which antigens or their regions should be prioritized for including them in a vaccine.Areas covered: The methods used by research groups so far which have focused on the functional analysis of P. vivax blood stage antigens have been reviewed here. A logical strategy orientated toward resolving two of the most commonly occurring problems in designing vaccines against this species has thus been proposed (i.e. the search for candidates and evaluating/ascertaining their functional role in the invasion of such molecules).Expert commentary: Advances in knowledge regarding P. vivax biology have been extremely slow. Only two key receptor-ligand interactions concerning merozoite entry to reticulocytes have been reported during the last 20 years: PvDBP1-DARC and PvRBP2b-CD71. Despite increasing knowledge about the parasite's intimate preference for its host cells, it has yet to be determined which regions of the merozoite molecules characterized to date meet the requirement of inducing protective immune responses effectively blocking heterologous parasite entry to human cells.

Amplification of Duffy binding protein-encoding gene allows Plasmodium vivax to evade host anti-DBP humoral immunity

Antigenic variation, the capacity to produce a range of variable antigens, is a well-described strategy of Plasmodium and other parasites to evade host immunity. Here, we show that gene amplification is an additional evasion mechanism used by Plasmodium vivax to escape humoral immunity targeting PvDBP, the key ligand involved in reticulocyte invasion. PvDBP gene amplification leads to increased mRNA levels and protects P. vivax in vitro against invasion inhibitory human monoclonal antibodies targeting a conserved binding domain of DBP. Patient samples suggest that parasites with increased pvdbp copy number are able to infect individuals with naturally acquired antibodies highly blocking the binding of PvDBP to the Duffy receptor. These results show that gene copy number variation affect the parasite’s ability to evade anti-PvDBP humoral immunity.

Duffy binding protein (DBP) of Plasmodium vivax is important for invasion and is a potential vaccine candidate. Here, the authors show that PvDBP gene amplification protects P vivax in vitro against invasion inhibitory human monoclonal antibodies and is associated to infection of patients with PvDBP binding inhibitory antibodies.

Plasmodium vivax Reticulocyte Binding Proteins for invasion into reticulocytes

Plasmodium vivax is responsible for most of the malaria infections outside Africa and is currently the predominant malaria parasite in countries under elimination programs. P. vivax preferentially enters young red cells called reticulocytes. Advances in understanding the molecular and cellular mechanisms of entry are hampered by the inability to grow large numbers of P. vivax parasites in a long-term in vitro culture. Recent progress in understanding the biology of the P. vivax Reticulocyte Binding Protein (PvRBPs) family of invasion ligands has led to the identification of a new invasion pathway into reticulocytes, an understanding of their structural architecture and PvRBPs as targets of the protective immune response to P. vivax infection. This review summarises current knowledge on the role of reticulocytes in P. vivax infection, the function of the PvRBP family of proteins in generating an immune response in human populations, and the characterization of anti-PvRBP antibodies in blocking parasite invasion.

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.

Transcriptome analysis of Plasmodium berghei during exo-erythrocytic development

BACKGROUND

The complex life cycle of malaria parasites requires well-orchestrated stage specific gene expression. In the vertebrate host the parasites grow and multiply by schizogony in two different environments: within erythrocytes and within hepatocytes. Whereas erythrocytic parasites are well-studied in this respect, relatively little is known about the exo-erythrocytic stages.

METHODS

In an attempt to fill this gap, genome wide RNA-seq analyses of various exo-erythrocytic stages of Plasmodium berghei including sporozoites, samples from a time-course of liver stage development and detached cells were performed. These latter contain infectious merozoites and represent the final step in exo-erythrocytic development.

RESULTS

The analysis represents the complete transcriptome of the entire life cycle of P. berghei parasites with temporal detailed analysis of the liver stage allowing comparison of gene expression across the progression of the life cycle. These RNA-seq data from different developmental stages were used to cluster genes with similar expression profiles, in order to infer their functions. A comparison with published data from other parasite stages confirmed stage-specific gene expression and revealed numerous genes that are expressed differentially in blood and exo-erythrocytic stages. One of the most exo-erythrocytic stage-specific genes was PBANKA_1003900, which has previously been annotated as a "gametocyte specific protein". The promoter of this gene drove high GFP expression in exo-erythrocytic stages, confirming its expression profile seen by RNA-seq.

CONCLUSIONS

The comparative analysis of the genome wide mRNA expression profiles of erythrocytic and different exo-erythrocytic stages could be used to improve the understanding of gene regulation in Plasmodium parasites and can be used to model exo-erythrocytic stage metabolic networks toward the identification of differences in metabolic processes during schizogony in erythrocytes and hepatocytes.

Identification and Immunological Characterization of the Ligand Domain of Plasmodium vivax Reticulocyte Binding Protein 1a

Background

Erythrocyte invasion by malaria parasites is essential for blood-stage development. Consequently, parasite proteins critically involved in erythrocyte invasion, such as the Plasmodium vivax reticulocyte binding proteins (RBPs) that mediate preferential invasion of reticulocytes, are considered potential vaccine targets. Thus, targeting the RBPs could prevent blood-stage infection and disease. The RBPs are large, and little is known about their functional domains and whether individuals naturally exposed to P. vivax acquire binding-inhibitory antibodies to these critical binding regions. This study aims to functionally and immunologically characterize Plasmodium vivax RBP1a.

Methods

Recombinant proteins of overlapping fragments of RBP1a were used to determine binding specificity to erythrocytes and immunogenicity in laboratory animals. The naturally acquired antibody response to these proteins was evaluated using serum samples from individuals in regions of endemicity.

Results

The N-terminal extracellular region, RBP1157-650 (RBP1:F8), was determined to bind both reticulocytes and normocytes, with a preference for immature reticulocytes. Antibodies elicited against rRBP1:F8 blocked binding between RBP1:F8 and erythrocytes. Naturally acquired anti-RBP1 binding-inhibitory antibodies were detected in serum specimens from P. vivax-exposed individuals from Papua New Guinea and Brazil.

Conclusion

Recombinant RBP1:F8 binds human erythrocytes, elicits artificially induced functional blocking antibodies, and is a target of naturally acquired binding-inhibitory antibodies.

Revisiting gametocyte biology in malaria parasites

Gametocytes are the only form of the malaria parasite that is transmissible to the mosquito vector. They are present at low levels in blood circulation and significant knowledge gaps exist in their biology. Recent reductions in the global malaria burden have brought the possibility of elimination and eradication, with renewed focus on malaria transmission biology as a basis for interventions. This review discusses recent insights into gametocyte biology in the major human malaria parasite, Plasmodium falciparum and related species.

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.

Identification and characterization of Pv50, a novel Plasmodium vivax merozoite surface protein

Background

Plasmodium vivax contains approximately 5400 coding genes, more than 40% of which code for hypothetical proteins that have not been functionally characterized. In a previous preliminary screening using pooled serum samples, numerous hypothetical proteins were selected from among those that were highly transcribed in the schizont-stage of parasites, and highly antigenic P. vivax candidates including hypothetical proteins were identified. However, their immunological and functional activities in P. vivax remain unclear. From these candidates, we investigated a P. vivax 50-kDa protein (Pv50, PVX_087140) containing a highly conserved signal peptide that shows high transcription levels in blood-stage parasites.

Results

Recombinant Pv50 was expressed in a cell-free expression system and used for IgG prevalence analysis of patients with vivax malaria and healthy individuals. Immune responses were analyzed in immunized mice and mouse antibodies were used to detect the subcellular localization of the protein in blood-stage parasites by immunofluorescence assay. A protein array method was used to evaluate protein-protein interactions to predict protein functional activities during the invasion of parasites into erythrocytes. Recombinant Pv50 showed IgG prevalence in patient samples with a sensitivity of 42.9% and specificity of 93.8% compared to that in healthy individuals. The non-cytophilic antibodies IgG1 and IgG3 were the major components involved in the antibody response in Pv50-immunized mice. Pv50 localized on the surface of merozoites and a specific interaction between Pv50 and PvMSP1 was detected, suggesting that Pv50-PvMSP1 forms a heterodimeric complex in P. vivax.

Conclusions

Increased immune responses caused by native P. vivax parasites were detected, confirming its immunogenic effects. This study provides a method for detecting new malaria antigens, and Pv50 may be a vivax malaria vaccine candidate with PvMSP1.

Electronic supplementary material

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

Gene regulatory networks on transfer entropy (GRNTE): a novel approach to reconstruct gene regulatory interactions applied to a case study for the plant pathogen Phytophthora infestans

Background

The increasing amounts of genomics data have helped in the understanding of the molecular dynamics of complex systems such as plant and animal diseases. However, transcriptional regulation, although playing a central role in the decision-making process of cellular systems, is still poorly understood. In this study, we linked expression data with mathematical models to infer gene regulatory networks (GRN). We present a simple yet effective method to estimate transcription factors’ GRNs from transcriptional data.

Method

We defined interactions between pairs of genes (edges in the GRN) as the partial mutual information between these genes that takes into account time and possible lags in time from one gene in relation to another. We call this method Gene Regulatory Networks on Transfer Entropy (GRNTE) and it corresponds to Granger causality for Gaussian variables in an autoregressive model. To evaluate the reconstruction accuracy of our method, we generated several sub-networks from the GRN of the eukaryotic yeast model, Saccharomyces cerevisae. Then, we applied this method using experimental data of the plant pathogen Phytophthora infestans. We evaluated the transcriptional expression levels of 48 transcription factors of P. infestans during its interaction with one moderately resistant and one susceptible cultivar of yellow potato (Solanum tuberosum group Phureja), using RT-qPCR. With these data, we reconstructed the regulatory network of P. infestans during its interaction with these hosts.

Results

We first evaluated the performance of our method, based on the transfer entropy (GRNTE), on eukaryotic datasets from the GRNs of the yeast S. cerevisae. Results suggest that GRNTE is comparable with the state-of-the-art methods when the parameters for edge detection are properly tuned. In the case of P. infestans, most of the genes considered in this study, showed a significant change in expression from the onset of the interaction (0 h post inoculum - hpi) to the later time-points post inoculation. Hierarchical clustering of the expression data discriminated two distinct periods during the infection: from 12 to 36 hpi and from 48 to 72 hpi for both the moderately resistant and susceptible cultivars. These distinct periods could be associated with two phases of the life cycle of the pathogen when infecting the host plant: the biotrophic and necrotrophic phases.

Conclusions

Here we presented an algorithmic solution to the problem of network reconstruction in time series data. This analytical perspective makes use of the dynamic nature of time series data as it relates to intrinsically dynamic processes such as transcription regulation, were multiple elements of the cell (e.g., transcription factors) act simultaneously and change over time. We applied the algorithm to study the regulatory network of P. infestans during its interaction with two hosts which differ in their level of resistance to the pathogen. Although the gene expression analysis did not show differences between the two hosts, the results of the GRN analyses evidenced rewiring of the genes’ interactions according to the resistance level of the host. This suggests that different regulatory processes are activated in response to different environmental cues. Applications of our methodology showed that it could reliably predict where to place edges in the transcriptional networks and sub-networks. The experimental approach used here can help provide insights on the biological role of these interactions on complex processes such as pathogenicity. The code used is available at https://github.com/jccastrog/GRNTE under GNU general public license 3.0.

Electronic supplementary material

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

DNA aptamers for the recognition of HMGB1 from Plasmodium falciparum

Rapid Diagnostic Tests (RDTs) for malaria are restricted to a few biomarkers and antibody-mediated detection. However, the expression of commonly used biomarkers varies geographically and the sensibility of immunodetection can be affected by batch-to-batch differences or limited thermal stability. In this study we aimed to overcome these limitations by identifying a potential biomarker and by developing molecular sensors based on aptamer technology. Using gene expression databases, ribosome profiling analysis, and structural modeling, we find that the High Mobility Group Box 1 protein (HMGB1) of Plasmodium falciparum is highly expressed, structurally stable, and present along all blood-stages of P. falciparum infection. To develop biosensors, we used in vitro evolution techniques to produce DNA aptamers for the recombinantly expressed HMG-box, the conserved domain of HMGB1. An evolutionary approach for evaluating the dynamics of aptamer populations suggested three predominant aptamer motifs. Representatives of the aptamer families were tested for binding parameters to the HMG-box domain using microscale thermophoresis and rapid kinetics. Dissociation constants of the aptamers varied over two orders of magnitude between nano- and micromolar ranges while the aptamer-HMG-box interaction occurred in a few seconds. The specificity of aptamer binding to the HMG-box of P. falciparum compared to its human homolog depended on pH conditions. Altogether, our study proposes HMGB1 as a candidate biomarker and a set of sensing aptamers that can be further developed into rapid diagnostic tests for P. falciparum detection.

Tail Wags the Dog? Functional Gene Classes Driving Genome-Wide GC Content in Plasmodium spp

Plasmodium parasites are valuable models to understand how nucleotide composition affects mutation, diversification, and adaptation. No other observed eukaryotes have undergone such large changes in genomic Guanine-Cytosine (GC) content as seen in the genus Plasmodium (∼30% within 35-40 Myr). Although mutational biases are known to influence GC content in the human-infective Plasmodium vivax and Plasmodium falciparum; no study has addressed how different gene functional classes contribute to genus-wide compositional changes, or if Plasmodium GC content variation is driven by natural selection. Here, we tested the hypothesis that certain gene processes and functions drive variation in global GC content between Plasmodium species. We performed a large-scale comparative genomic analysis using the genomes and predicted genes of 17 Plasmodium species encompassing a wide genomic GC content range. Genic GC content was sorted and divided into ten equally sized quantiles that were then assessed for functional enrichment classes. In agreement that selection on gene classes may drive genomic GC content, trans-membrane proteins were enriched within extreme GC content quantiles (Q1 and Q10). Specifically, variant surface antigens, which primarily interact with vertebrate immune systems, showed skewed GC content distributions compared with other trans-membrane proteins. Although a definitive causation linking GC content, expression, and positive selection within variant surface antigens from Plasmodium vivax, Plasmodium berghei, and Plasmodium falciparum could not be established, we found that regardless of genomic nucleotide composition, genic GC content and expression were positively correlated during trophozoite stages. Overall, these data suggest that, alongside mutational biases, functional protein classes drive Plasmodium GC content change.

A fast and cost-effective microsampling protocol incorporating reduced animal usage for time-series transcriptomics in rodent malaria parasites

BACKGROUND

The transcriptional regulation that occurs in malaria parasites during the erythrocytic stages of infection can be studied in vivo with rodent malaria parasites propagated in mice. Time-series transcriptome profiling commonly involves the euthanasia of groups of mice at specific time points followed by the extraction of parasite RNA from whole blood samples. Current methodologies for parasite RNA extraction involve several steps and when multiple time points are profiled, these protocols are laborious, time-consuming, and require the euthanization of large cohorts of mice.

RESULTS

A simplified protocol has been designed for parasite RNA extraction from blood volumes as low as 20 μL (microsamples), serially bled from mice via tail snips and directly lysed with TRIzol reagent. Gene expression data derived from microsampling using RNA-seq were closely matched to those derived from larger volumes of leucocyte-depleted and saponin-treated blood obtained from euthanized mice with high reproducibility between biological replicates. Transcriptome profiling of microsamples taken at different time points during the intra-erythrocytic developmental cycle of the rodent malaria parasite Plasmodium vinckei revealed the transcriptional cascade commonly observed in malaria parasites.

CONCLUSIONS

Microsampling is a quick, robust and cost-efficient approach to sample collection for in vivo time-series transcriptomic studies in rodent malaria parasites.

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.

Forces acting on codon bias in malaria parasites

Malaria parasite genomes have a range of codon biases, with Plasmodium falciparum one of the most AT-biased genomes known. We examined the make up of synonymous coding sites and stop codons in the core genomes of representative malaria parasites, showing first that local DNA context influences codon bias similarly across P. falciparum, P. vivax and P. berghei, with suppression of CpG dinucleotides and enhancement of CpC dinucleotides, both within and aross codons. Intense asexual phase gene expression in P. falciparum and P. berghei is associated with increased A3:G3 bias but reduced T3:C3 bias at 2-fold sites, consistent with adaptation of codons to tRNA pools and avoidance of wobble tRNA interactions that potentially slow down translation. In highly expressed genes, the A3:G3 ratio can exceed 30-fold while the T3:C3 ratio can be less than 1, according to the encoded amino acid and subsequent base. Lysine codons (AAA/G) show distinctive behaviour with substantially reduced A3:G3 bias in highly expressed genes, perhaps because of selection against frameshifting when the AAA codon is followed by another adenine. Intense expression is also associated with a strong bias towards TAA stop codons (found in 94% and 89% of highly expressed P. falciparum and P. berghei genes respectively) and a proportional rise in the TAAA stop ‘tetranucleotide’. The presence of these expression-linked effects in the relatively AT-rich malaria parasite species adds weight to the suggestion that AT-richness in the Plasmodium genus might be a fitness adaptation. Potential explanations for the relative lack of codon bias in P. vivax include the distinct features of its lifecycle and its effective population size over evolutionary time.