New insights into the Plasmodium vivax transcriptome using RNA-Seq

De novo assembly of plasmodium interspersed repeat (pir) genes from Plasmodium vivax RNAseq data suggests geographic conservation of sub-family transcription

BACKGROUND

The plasmodium interspersed repeats (pir) multigene family is found across malaria parasite genomes, first discovered in the human-infecting species Plasmodium vivax, where they were initially named the virs. Their function remains unknown, although studies have suggested a role in virulence of the asexual blood stages. Sub-families of the P. vivax pir/virs have been identified, and are found in isolates from across the world, however their transcription at different localities and in different stages of the life cycle have not been quantified. Multiple transcriptomic studies of the parasite have been conducted, but many map the pir reads to existing reference genomes (as part of standard bioinformatic practice), which may miss members of the multigene family due to its inherent variability. This obscures our understanding of how the pir sub-families in P. vivax may be contributing to human/vector infection.

RESULTS

To overcome the issue of hidden pir diversity from utilising a reference genome, we employed de novo transcriptome assembly to construct the pir 'reference' of different parasite isolates from published and novel RNAseq datasets. For this purpose, a pipeline was written in Nextflow, and first tested on data from the rodent-infecting P. c. chabaudi parasite to ascertain its efficacy on a sample with a full, genome-based set of pir gene sequences. The pipeline assembled hundreds of pirs from the studies included. By performing BLAST sequence identity comparisons with reference genome pirs (including P. vivax and related species) we found a clustered network of transcripts which corresponded well with prior sub-family annotations, albeit requiring some updated nomenclature. Mapping the RNAseq datasets to the de novo transcriptome references revealed that the transcription of these updated pir gene sub-families is generally consistent across the different geographical regions. From this transcriptional quantification, a time course of mosquito bloodmeals (after feeding on an infected patient) highlighted the first evidence of ookinete stage pir transcription in a human-infective malaria parasite.

CONCLUSIONS

De novo transcriptome assembly is a valuable tool for understanding highly variable multigene families from Plasmodium spp., and with pipeline software these can be applied more easily and at scale. Despite a global distribution, P. vivax has a conserved pir sub-family structure-both in terms of genome copy number and transcription. We suggest that this indicates important roles of the distinct sub-families, or a genetic mechanism maintaining their preservation. Furthermore, a burst of pir transcription in the mosquito stages of development is the first glint of ookinete pir expression for a human-infective malaria parasite, suggesting a role for the gene family at a new stage of the lifecycle.

A tale of two parasites: a glimpse into the RNA methylome of patient-derived Plasmodium falciparum and Plasmodium vivax isolates

BACKGROUND

Understanding the molecular mechanisms of the malarial parasites in hosts is crucial for developing effective treatments. Epitranscriptomic research on pathogens has unveiled the significance of RNA methylation in gene regulation and pathogenesis. This is the first report investigating methylation signatures and alternative splicing events using Nanopore Direct RNA Sequencing to single-base resolution in Plasmodium falciparum and Plasmodium vivax clinical isolates with hepatic dysfunction complications.

METHODS

Direct RNA Sequencing using Nanopore from clinical isolates of P. falciparum and P. vivax showing hepatic dysfunction manifestation was performed. Subsequently, transcriptome reconstruction using FLAIR and transcript classification using SQANTI3, followed by methylation detection using CHEUI and m6Anet to identify N6-methyladenosine (m6A) and 5-methylcytosine (m5C) methylation signatures, was done. The alternative splicing events from both the datasets were documented.

RESULTS

The reference genome of Plasmodium reports > 5000 genes out of which ~ 50% was identified as expressed in the two sequenced isolates, including novel isoforms and intergenic transcripts, highlighting extensive transcriptome diversity. The distinct RNA methylation profiles of m6A and m5C from the expressed transcripts were observed in sense, Natural Antisense Transcripts (NATs) and intergenic categories hinting at species-specific regulatory mechanisms. Dual modification events were observed in a significant number of transcripts in both the parasites. Modified transcripts originating from apicoplast and mitochondrial genomes have also been detected. These modifications are unevenly present in the annotated regions of the mRNA, potentially influencing mRNA export and translation. Several splicing events were observed, with alternative 3' and 5' end splicing predominating in the datasets suggesting differences in translational kinetics and possible protein characteristics in these disease conditions.

CONCLUSION

The data shows the presence of modified sense, NATs and alternatively spliced transcripts. These phenomena together suggest the presence of multiple regulatory layers which decides the post-translational proteome of the parasites in particular disease conditions. Studies like these will help to decipher the post-translational environments of malaria parasites in vivo and elucidate their inherent proteome plasticity, thus allowing the conceptualization of novel strategies for interventions.

Liver stage P. falciparum antigens highly targeted by CD4+ T cells in malaria-exposed Ugandan children

T cell responses against liver stage Plasmodium help protect against reinfection, but the antigens and epitopes targeted by these T cells are largely unknown. This knowledge gap has impeded mechanistic studies to identify the effector functions most critical for protection. We performed a bioinformatic analysis of gene expression datasets to identify plasmodial genes that are highly and selectively expressed during liver stage infection and predict epitopes within them likely to bind MHC-II molecules prevalent in Uganda. We then tested their recognition by malaria-exposed Ugandan children. In over two-thirds of the children, we detected a peripheral blood CD4+ T cell response to one or more antigens. The most highly targeted antigen, LISP1, contained several epitopes, including one that was promiscuously presented and recognized by most participants. These novel liver stage P. falciparum epitopes should be explored as potential vaccine targets and will facilitate the development of tools to interrogate antimalarial immunity at the single-cell level and inform future vaccine development efforts.

Population genomics and transcriptomics of Plasmodium falciparum in Cambodia and Vietnam uncover key components of the artemisinin resistance genetic background

The emergence of Plasmodium falciparum parasites resistant to artemisinins compromises the efficacy of Artemisinin Combination Therapies (ACTs), the global first-line malaria treatment. Artemisinin resistance is a complex genetic trait in which nonsynonymous SNPs in PfK13 cooperate with other genetic variations. Here, we present population genomic/transcriptomic analyses of P. falciparum collected from patients with uncomplicated malaria in Cambodia and Vietnam between 2018 and 2020. Besides the PfK13 SNPs, several polymorphisms, including nonsynonymous SNPs (N1131I and N821K) in PfRad5 and an intronic SNP in PfWD11 (WD40 repeat-containing protein on chromosome 11), appear to be associated with artemisinin resistance, possibly as new markers. There is also a defined set of genes whose steady-state levels of mRNA and/or splice variants or antisense transcripts correlate with artemisinin resistance at the base level. In vivo transcriptional responses to artemisinins indicate the resistant parasite's capacity to decelerate its intraerythrocytic developmental cycle (IDC), which can contribute to the resistant phenotype. During this response, PfRAD5 and PfWD11 upregulate their respective alternatively/aberrantly spliced isoforms, suggesting their contribution to the protective response to artemisinins. PfRAD5 and PfWD11 appear under selective pressure in the Greater Mekong Sub-region over the last decade, suggesting their role in the genetic background of the artemisinin resistance.

MYST regulates DNA repair and forms a NuA4-like complex in the malaria parasite Plasmodium falciparum

Histone lysine acetyltransferase MYST-associated NuA4 complex is conserved from yeast to humans and plays key roles in cell cycle regulation, gene transcription, and DNA replication/repair. Here, we identified a Plasmodium falciparum MYST-associated complex, PfNuA4, which contains 11 of the 13 conserved NuA4 subunits. Reciprocal pulldowns using PfEAF2, a shared component between the NuA4 and SWR1 complexes, not only confirmed the PfNuA4 complex but also identified the PfSWR1 complex, a histone remodeling complex, although their identities are low compared to the homologs in yeast or humans. Notably, both H2A.Z/H2B.Z were associated with the PfSWR1 complex, indicating that this complex is involved in the deposition of H2A.Z/H2B.Z, the variant histone pair that is enriched in the activated promoters. Overexpression of PfMYST resulted in earlier expression of genes involved in cell cycle regulation, DNA replication, and merozoite invasion, and upregulation of the genes related to antigenic variation and DNA repair. Consistently, PfMYST overexpression led to high basal phosphorylated PfH2A (γ-PfH2A), the mark of DNA double-strand breaks, and conferred protection against genotoxic agent methyl methanesulfonate (MMS), X-rays, and artemisinin, the first-line antimalarial drug. In contrast, the knockdown of PfMYST caused a delayed parasite recovery upon MMS treatment. MMS induced the gradual disappearance of PfMYST in the cytoplasm and concomitant accumulation of PfMYST in the nucleus, suggesting cytoplasm-nucleus shuttling of PfMYST. Meanwhile, PfMYST colocalized with the γ-PfH2A, indicating PfMYST was recruited to the DNA damage sites. Collectively, PfMYST plays critical roles in cell cycle regulation, gene transcription, and DNA replication/DNA repair in this low-branching parasitic protist.IMPORTANCEUnderstanding gene regulation and DNA repair in malaria parasites is critical for identifying targets for antimalarials. This study found PfNuA4, a PfMYST-associated, histone modifier complex, and PfSWR1, a chromatin remodeling complex in malaria parasite Plasmodium falciparum. These complexes are divergent due to the low identities compared to their homologs from yeast and humans. Furthermore, overexpression of PfMYST resulted in substantial transcriptomic changes, indicating that PfMYST is involved in regulating the cell cycle, antigenic variation, and DNA replication/repair. Consistently, PfMYST was found to protect against DNA damage caused by the genotoxic agent methyl methanesulfonate, X-rays, and artemisinin, the first-line antimalarial drug. Additionally, DNA damage led to the relocation of cytoplasmic PfMYST to the nucleus and colocalization of PfMYST with γ-PfH2A, the mark of DNA damage. In summary, this study demonstrated that the PfMYST complex has critical functions in regulating cell cycle, antigenic variation, and DNA replication/DNA repair in P. falciparum.

Characterisation of Plasmodium vivax lactate dehydrogenase dynamics in P. vivax infections

Plasmodium vivax lactate dehydrogenase (PvLDH) is an essential enzyme in the glycolytic pathway of P. vivax. It is widely used as a diagnostic biomarker and a measure of total-body parasite biomass in vivax malaria. However, the dynamics of PvLDH remains poorly understood. Here, we developed mathematical models that capture parasite and matrix PvLDH dynamics in ex vivo culture and the human host. We estimated key biological parameters characterising in vivo PvLDH dynamics based on longitudinal data of parasitemia and PvLDH concentration collected from P. vivax-infected humans, with the estimates informed by the ex vivo data as prior knowledge in a Bayesian hierarchical framework. We found that the in vivo accumulation rate of intraerythrocytic PvLDH peaks at 10-20 h post-invasion (late ring stage) with a median estimate of intraerythrocytic PvLDH mass at the end of the life cycle to be 9.4 × 10-3ng. We also found that the median estimate of in vivo PvLDH half-life was approximately 21.9 h. Our findings provide a foundation with which to advance our quantitative understanding of P. vivax biology and will facilitate the improvement of PvLDH-based diagnostic tools.

Comparative transcriptomics reveal differential gene expression among Plasmodium vivax geographical isolates and implications on erythrocyte invasion mechanisms

The documentation of Plasmodium vivax malaria across Africa especially in regions where Duffy negatives are dominant suggests possibly alternative erythrocyte invasion mechanisms. While the transcriptomes of the Southeast Asian and South American P. vivax are well documented, the gene expression profile of P. vivax in Africa is unclear. In this study, we examined the expression of 4,404 gene transcripts belong to 12 functional groups and 43 erythrocyte binding gene candidates in Ethiopian isolates and compared them with the Cambodian and Brazilian P. vivax transcriptomes. Overall, there were 10-26% differences in the gene expression profile amongst geographical isolates, with the Ethiopian and Cambodian P. vivax being most similar. Majority of the gene transcripts involved in protein transportation, housekeeping, and host interaction were highly transcribed in the Ethiopian isolates. Members of the reticulocyte binding protein PvRBP2a and PvRBP3 expressed six-fold higher than Duffy binding protein PvDBP1 and 60-fold higher than PvEBP/DBP2 in the Ethiopian isolates. Other genes including PvMSP3.8, PvMSP3.9, PvTRAG2, PvTRAG14, and PvTRAG22 also showed relatively high expression. Differential expression patterns were observed among geographical isolates, e.g., PvDBP1 and PvEBP/DBP2 were highly expressed in the Cambodian but not the Brazilian and Ethiopian isolates, whereas PvRBP2a and PvRBP2b showed higher expression in the Ethiopian and Cambodian than the Brazilian isolates. Compared to Pvs25, gametocyte genes including PvAP2-G, PvGAP (female gametocytes), and Pvs47 (male gametocytes) were highly expressed across geographical samples.

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.

Unusually Divergent Ubiquitin Genes and Proteins in Plasmodium Species

Ubiquitin is an extraordinarily highly conserved 76 amino acid protein encoded by three different types of gene, where the primary translation products are fusions either of ubiquitin with one of two ribosomal proteins (RPs) or of multiple ubiquitin monomers from head to tail. Here, we investigate the evolution of ubiquitin genes in mammalian malaria parasites (Plasmodium species). The ubiquitin encoded by the RPS27a fusion gene is highly divergent, as previously found in a variety of protists. However, we also find that two other forms of divergent ubiquitin sequence, each previously thought to be extremely rare, have arisen recently during the divergence of Plasmodium subgenera. On two occasions, in two distinct lineages, the ubiquitin encoded by the RPL40 fusion gene has rapidly diverged. In addition, in one of these lineages, the polyubiquitin genes have undergone a single codon insertion, previously considered a unique feature of Rhizaria. There has been disagreement whether the multiple ubiquitin coding repeats within a genome exhibit concerted evolution or undergo a birth-and-death process; the Plasmodium ubiquitin genes show clear signs of concerted evolution, including the spread of this codon insertion to multiple repeats within the polyubiquitin gene.

Comparative transcriptomics reveal differential gene expression in Plasmodium vivax geographical isolates and implications on erythrocyte invasion mechanisms

Plasmodium vivax uses Duffy binding protein (PvDBP1) to bind to the Duffy Antigen-Chemokine Receptor (DARC) to invade human erythrocytes. Individuals who lack DARC expression (Duffy-negative) are thought to be resistance to P. vivax. In recent years, P. vivax malaria is becoming more prevalent in Africa with a portion of these cases detected in Duffy-negatives. Apart from DBP1, members of the reticulocyte binding protein (RBP) and tryptophan-rich antigen (TRAg) families may also play a role in erythrocyte invasion. While the transcriptomes of the Southeast Asian and South American P. vivax are well documented, the gene expression profile of P. vivax in Africa and more specifically the expression level of several erythrocyte binding gene candidates as compared to DBP1 are largely unknown. This paper characterized the first P. vivax transcriptome in Africa and compared with those from the Southeast Asian and South American isolates. The expression of 4,404 gene transcripts belong to 12 functional groups including 43 specific erythrocyte binding gene candidates were examined. Overall, there were 10-26% differences in the gene expression profile amongst the geographical isolates, with the Ethiopian and Cambodian P. vivax being most similar. Majority of the gene transcripts involved in protein transportation, housekeeping, and host interaction were highly transcribed in the Ethiopian P. vivax. Erythrocyte binding genes including PvRBP2a and PvRBP3 expressed six-fold higher than PvDBP1and 60-fold higher than PvEBP/DBP2. Other genes including PvRBP1a, PvMSP3.8, PvMSP3.9, PvTRAG2, PvTRAG14, and PvTRAG22 also showed relatively high expression. Differential expression was observed among geographical isolates, e.g., PvDBP1 and PvEBP/DBP2 were highly expressed in the Cambodian but not the Brazilian and Ethiopian isolates, whereas PvRBP2a and PvRBP2b showed higher expression in the Ethiopian and Cambodian than the Brazilian isolates. Compared to Pvs25, the standard biomarker for detecting female gametocytes, PvAP2-G (PVP01_1440800), GAP (PVP01_1403000), and Pvs47 (PVP01_1208000) were highly expressed across geographical samples. These findings provide an important baseline for future comparisons of P. vivax transcriptomes from Duffy-negative infections and highlight potential biomarkers for improved gametocyte detection.

Evaluation of transmission-blocking potential of Pv22 using clinical Plasmodium vivax infections and transgenic Plasmodium berghei

Antigens expressed during the sexual development of malaria parasites are transmission-blocking vaccine (TBV) targets. Pb22, a protein expressed and localized to the plasma membrane of gametes and ookinetes in Plasmodium berghei, is an excellent TBV candidate. Here, we evaluated the TB potential of the Plasmodium vivax ortholog Pv22 using a transgenic P. berghei parasite line and P. vivax clinical isolates. The full-length recombinant Pv22 (rPv22) protein was produced and used to immunize mice and rabbits to obtain antibodies. We generated a transgenic P. berghei line (TrPv22Pb) by inserting the pv22 gene into the pb22 locus and showed that Pv22 expression completely rescued the defects in male gametogenesis of the pb22 deletion parasite. Since Pv22 in the transgenic parasite showed similar expression and localization patterns to Pb22, we used the TrPv22Pb parasite as a surrogate to evaluate the TB potential of Pv22. In mosquito feeding assays, mosquitoes feeding on rPv22-immunized mice infected with TrPv22Pb parasites showed a 49.3-53.3 % reduction in the oocyst density compared to the control group. In vitro assays showed that the rPv22 immune sera significantly inhibited exflagellation and ookinete formation of the TrPv22Pb parasites. In a direct membrane feeding assay using three clinical P. vivax isolates, the rabbit anti-rPv22 antibodies also significantly decreased the oocyst density by 53.7, 30.2, and 26.2 %, respectively. This study demonstrated the feasibility of using transgenic P. berghei parasites expressing P. vivax antigens as a potential tool to evaluate TBV candidates. However, the much weaker TB activity of Pv22 obtained from two complementary assays suggest that Pv22 may not be a promising TBV candidate for P. vivax.

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.

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.

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.

Comparative analysis of the kinomes of Plasmodium falciparum, Plasmodium vivax and their host Homo sapiens

Background

Novel antimalarials should be effective across all species of malaria parasites that infect humans, especially the two species that bear the most impact, Plasmodium falciparum and Plasmodium vivax. Protein kinases encoded by pathogens, as well as host kinases required for survival of intracellular pathogens, carry considerable potential as targets for antimalarial intervention (Adderley et al. Trends Parasitol 37:508–524, 2021;  Wei et al. Cell Rep Med 2:100423, 2021). To date, no comprehensive P. vivax kinome assembly has been conducted; and the P. falciparum kinome, first assembled in 2004, requires an update. The present study, aimed to fill these gaps, utilises a recently published structurally-validated multiple sequence alignment (MSA) of the human kinome (Modi et al. Sci Rep 9:19790, 2019). This MSA is used as a scaffold to assist the alignment of all protein kinase sequences from P. falciparum and P. vivax, and (where possible) their assignment to specific kinase groups/families.

Results

We were able to assign six P. falciparum previously classified as OPK or ‘orphans’ (i.e. with no clear phylogenetic relation to any of the established ePK groups) to one of the aforementioned ePK groups. Direct phylogenetic comparison established that despite an overall high level of similarity between the P. falciparum and P. vivax kinomes, which will help in selecting targets for intervention, there are differences that may underlie the biological specificities of these species. Furthermore, we highlight a number of Plasmodium kinases that have a surprisingly high level of similarity with their human counterparts and therefore not well suited as targets for drug discovery.

Conclusions

Direct comparison of the kinomes of Homo sapiens, P. falciparum and P. vivax sheds additional light on the previously documented divergence of many P. falciparum and P. vivax kinases from those of their human host. We provide the first direct kinome comparison between the phylogenetically distinct species of P. falciparum and P. vivax, illustrating the key similarities and differences which must be considered in the context of kinase-directed antimalarial drug discovery, and discuss the divergences and similarities between the human and Plasmodium kinomes to inform future searches for selective antimalarial intervention.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08457-0.

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.

A unique GCN5 histone acetyltransferase complex controls erythrocyte invasion and virulence in the malaria parasite Plasmodium falciparum

The histone acetyltransferase GCN5-associated SAGA complex is evolutionarily conserved from yeast to human and functions as a general transcription co-activator in global gene regulation. In this study, we identified a divergent GCN5 complex in Plasmodium falciparum, which contains two plant homeodomain (PHD) proteins (PfPHD1 and PfPHD2) and a plant apetela2 (AP2)-domain transcription factor (PfAP2-LT). To dissect the functions of the PfGCN5 complex, we generated parasite lines with either the bromodomain in PfGCN5 or the PHD domain in PfPHD1 deleted. The two deletion mutants closely phenocopied each other, exhibiting significantly reduced merozoite invasion of erythrocytes and elevated sexual conversion. These domain deletions caused dramatic decreases not only in histone H3K9 acetylation but also in H3K4 trimethylation, indicating synergistic crosstalk between the two euchromatin marks. Domain deletion in either PfGCN5 or PfPHD1 profoundly disturbed the global transcription pattern, causing altered expression of more than 60% of the genes. At the schizont stage, these domain deletions were linked to specific down-regulation of merozoite genes involved in erythrocyte invasion, many of which contain the AP2-LT binding motif and are also regulated by AP2-I and BDP1, suggesting targeted recruitment of the PfGCN5 complex to the invasion genes by these specific factors. Conversely, at the ring stage, PfGCN5 or PfPHD1 domain deletions disrupted the mutually exclusive expression pattern of the entire var gene family, which encodes the virulent factor PfEMP1. Correlation analysis between the chromatin state and alteration of gene expression demonstrated that up- and down-regulated genes in these mutants are highly correlated with the silent and active chromatin states in the wild-type parasite, respectively. Collectively, the PfGCN5 complex represents a novel HAT complex with a unique subunit composition including an AP2 transcription factor, which signifies a new paradigm for targeting the co-activator complex to regulate general and parasite-specific cellular processes in this low-branching parasitic protist.

Epigenetic regulation of gene expression plays essential roles in orchestrating the general and parasite-specific cellular pathways in the malaria parasite Plasmodium falciparum. To better understand the epigenetic mechanisms in this parasite, we characterized the histone acetyltransferase GCN5-mediated transcription regulation during intraerythrocytic development of the parasite. Using tandem affinity purification and proteomic characterization, we identified that the PfGCN5-associated complex contains nine core components, including two PHD domain proteins (PfPHD1 and PfPHD2) and an AP2-domain transcription factor, which is divergent from the canonical GCN5 complexes evolutionarily conserved from yeast to human. To understand the functions of the PfGCN5 complex, we performed domain deletions in two subunits of this complex, PfGCN5 and PfPHD1. We found that the two deletion mutants displayed very similar growth phenotypes, including significantly reduced merozoite invasion rates and elevated sexual conversion. These two mutants were associated with dramatic decreases in histone H3K9 acetylation and H3K4 trimethylation, which led to global changes in chromatin states and gene expression. Consistent with the phenotypes, genes significantly affected by the PfGCN5 and PfPHD1 gene disruption include those participating in parasite-specific pathways such as invasion, virulence, and sexual development. In conclusion, this study presents a new model of the PfGCN5 complex for targeting the co-activator complex to regulate general and parasite-specific cellular processes in this low-branching parasitic protist.

Identification of a protein unique to the genus Plasmodium that contains a WD40 repeat domain and extensive low-complexity sequence

Proteins containing WD40 domains play important roles in the formation of multiprotein complexes. Little is known about WD40 proteins in the malaria parasite. This report contains the initial description of a WD40 protein that is unique to the genus Plasmodium and possibly closely related genera. The N-terminal portion of this protein consists of seven WD40 repeats that are highly conserved in all Plasmodium species. Following the N-terminal region is a central region that is conserved within the major Plasmodium clades, such as parasites of great apes, monkeys, rodents, and birds, but partially conserved across all Plasmodium species. This central region contains extensive low-complexity sequence and is predicted to have a disordered structure. Proteins with disordered structure generally function in molecular interactions. The C-terminal region is semi-conserved across all Plasmodium species and has no notable features. This WD40 repeat protein likely functions in some aspect of parasite biology that is unique to Plasmodium and this uniqueness makes the protein a possible target for therapeutic intervention.

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.

In silico identification of novel open reading frames in Plasmodium falciparum oocyte and salivary gland sporozoites using proteogenomics framework

Background

Plasmodium falciparum causes the deadliest form of malaria, which remains one of the most prevalent infectious diseases. Unfortunately, the only licensed vaccine showed limited protection and resistance to anti-malarial drug is increasing, which can be largely attributed to the biological complexity of the parasite’s life cycle. The progression from one developmental stage to another in P. falciparum involves drastic changes in gene expressions, where its infectivity to human hosts varies greatly depending on the stage. Approaches to identify candidate genes that are responsible for the development of infectivity to human hosts typically involve differential gene expression analysis between stages. However, the detection may be limited to annotated proteins and open reading frames (ORFs) predicted using restrictive criteria.

Methods

The above problem is particularly relevant for P. falciparum; whose genome annotation is relatively incomplete given its clinical significance. In this work, systems proteogenomics approach was used to address this challenge, as it allows computational detection of unannotated, novel Open Reading Frames (nORFs), which are neglected by conventional analyses. Two pairs of transcriptome/proteome were obtained from a previous study where one was collected in the mosquito-infectious oocyst sporozoite stage, and the other in the salivary gland sporozoite stage with human infectivity. They were then re-analysed using the proteogenomics framework to identify nORFs in each stage.

Results

Translational products of nORFs that map to antisense, intergenic, intronic, 3′ UTR and 5′ UTR regions, as well as alternative reading frames of canonical proteins were detected. Some of these nORFs also showed differential expression between the two life cycle stages studied. Their regulatory roles were explored through further bioinformatics analyses including the expression regulation on the parent reference genes, in silico structure prediction, and gene ontology term enrichment analysis.

Conclusion

The identification of nORFs in P. falciparum sporozoites highlights the biological complexity of the parasite. Although the analyses are solely computational, these results provide a starting point for further experimental validation of the existence and functional roles of these nORFs,

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.

Humic-acid-driven escape from eye parasites revealed by RNA-seq and target-specific metabarcoding

Background

Next generation sequencing (NGS) technologies are extensively used to dissect the molecular mechanisms of host-parasite interactions in human pathogens. However, ecological studies have yet to fully exploit the power of NGS as a rich source for formulating and testing new hypotheses.

Methods

We studied Eurasian perch (Perca fluviatilis) and its eye parasite (Trematoda, Diplostomidae) communities in 14 lakes that differed in humic content in order to explore host-parasite-environment interactions. We hypothesised that high humic content along with low pH would decrease the abundance of the intermediate hosts (gastropods), thus limiting the occurrence of diplostomid parasites in humic lakes. This hypothesis was initially invoked by whole eye RNA-seq data analysis and subsequently tested using PCR-based detection and a novel targeted metabarcoding approach.

Results

Whole eye transcriptome results revealed overexpression of immune-related genes and the presence of eye parasite sequences in RNA-seq data obtained from perch living in clear-water lakes. Both PCR-based and targeted-metabarcoding approach showed that perch from humic lakes were completely free from diplostomid parasites, while the prevalence of eye flukes in clear-water lakes that contain low amounts of humic substances was close to 100%, with the majority of NGS reads assigned to Tylodelphys clavata.

Conclusions

High intraspecific diversity of T. clavata indicates that massively parallel sequencing of naturally pooled samples represents an efficient and powerful strategy for shedding light on cryptic diversity of eye parasites. Our results demonstrate that perch populations in clear-water lakes experience contrasting eye parasite pressure compared to those from humic lakes, which is reflected by prevalent differences in the expression of immune-related genes in the eye. This study highlights the utility of NGS to discover novel host-parasite-environment interactions and provide unprecedented power to characterize the molecular diversity of cryptic parasites.

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.

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.

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.

TDR Targets 6: driving drug discovery for human pathogens through intensive chemogenomic data integration

The volume of biological, chemical and functional data deposited in the public domain is growing rapidly, thanks to next generation sequencing and highly-automated screening technologies. These datasets represent invaluable resources for drug discovery, particularly for less studied neglected disease pathogens. To leverage these datasets, smart and intensive data integration is required to guide computational inferences across diverse organisms. The TDR Targets chemogenomics resource integrates genomic data from human pathogens and model organisms along with information on bioactive compounds and their annotated activities. This report highlights the latest updates on the available data and functionality in TDR Targets 6. Based on chemogenomic network models providing links between inhibitors and targets, the database now incorporates network-driven target prioritizations, and novel visualizations of network subgraphs displaying chemical- and target-similarity neighborhoods along with associated target-compound bioactivity links. Available data can be browsed and queried through a new user interface, that allow users to perform prioritizations of protein targets and chemical inhibitors. As such, TDR Targets now facilitates the investigation of drug repurposing against pathogen targets, which can potentially help in identifying candidate targets for bioactive compounds with previously unknown targets. TDR Targets is available at https://tdrtargets.org.

The transcriptome of the rumen ciliate Entodinium caudatum reveals some of its metabolic features

Background

Rumen ciliates play important roles in rumen function by digesting and fermenting feed and shaping the rumen microbiome. However, they remain poorly understood due to the lack of definitive direct evidence without influence by prokaryotes (including symbionts) in co-cultures or the rumen. In this study, we used RNA-Seq to characterize the transcriptome of Entodinium caudatum, the most predominant and representative rumen ciliate species.

Results

Of a large number of transcripts, > 12,000 were annotated to the curated genes in the NR, UniProt, and GO databases. Numerous CAZymes (including lysozyme and chitinase) and peptidases were represented in the transcriptome. This study revealed the ability of E. caudatum to depolymerize starch, hemicellulose, pectin, and the polysaccharides of the bacterial and fungal cell wall, and to degrade proteins. Many signaling pathways, including the ones that have been shown to function in E. caudatum, were represented by many transcripts. The transcriptome also revealed the expression of the genes involved in symbiosis, detoxification of reactive oxygen species, and the electron-transport chain. Overall, the transcriptomic evidence is consistent with some of the previous premises about E. caudatum. However, the identification of specific genes, such as those encoding lysozyme, peptidases, and other enzymes unique to rumen ciliates might be targeted to develop specific and effective inhibitors to improve nitrogen utilization efficiency by controlling the activity and growth of rumen ciliates. The transcriptomic data will also help the assembly and annotation in future genomic sequencing of E. caudatum.

Conclusion

As the first transcriptome of a single species of rumen ciliates ever sequenced, it provides direct evidence for the substrate spectrum, fermentation pathways, ability to respond to various biotic and abiotic stimuli, and other physiological and ecological features of E. caudatum. The presence and expression of the genes involved in the lysis and degradation of microbial cells highlight the dependence of E. caudatum on engulfment of other rumen microbes for its survival and growth. These genes may be explored in future research to develop targeted control of Entodinium species in the rumen. The transcriptome can also facilitate future genomic studies of E. caudatum and other related rumen ciliates.

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.

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.

Vivax Sporozoite Consortium, Muller I, Jex AR, Kappe SHI, Mikolajczak SA, Sattabongkot J, Patrapuvich R, Lindner S, Flannery EL, Koepfli C, Ansell B, Lerch A, Emery-Corbin SJ, Charnaud S, Smith J, Merrienne N, Swearingen KE, Moritz RL, Petter M, Duffy MF, Chuenchob V. 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.

Alternative Splicing in Apicomplexan Parasites

Alternative splicing is a widespread, essential, and complex component of gene regulation. Apicomplexan parasites have long been recognized to produce alternatively spliced transcripts for some genes and can produce multiple protein products that are essential for parasite growth.

Alternative splicing is a widespread, essential, and complex component of gene regulation. Apicomplexan parasites have long been recognized to produce alternatively spliced transcripts for some genes and can produce multiple protein products that are essential for parasite growth. Recent approaches are now providing more wide-ranging surveys of the extent of alternative splicing; some indicate that alternative splicing is less widespread than in other model eukaryotes, whereas others suggest levels comparable to those of previously studied groups. In many cases, apicomplexan alternative splicing events appear not to generate multiple alternative proteins but instead produce aberrant or noncoding transcripts. Nonetheless, appropriate regulation of alternative splicing is clearly essential in Plasmodium and Toxoplasma parasites, suggesting a biological role for at least some of the alternative splicing observed. Several studies have now disrupted conserved regulators of alternative splicing and demonstrated lethal effects in apicomplexans. This minireview discusses methods to accurately determine the extent of alternative splicing in Apicomplexa and discuss potential biological roles for this conserved process in a phylum of parasites with compact genomes.

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.

Plasmodium vivax readiness to transmit: implication for malaria eradication

Background

The lack of a continuous long-term in vitro culture system for Plasmodium vivax severely limits our knowledge of pathophysiology of the most widespread malaria parasite. To gain direct understanding of P. vivax human infections, we used Next Generation Sequencing data mining to unravel parasite in vivo expression profiles for P. vivax, and P. falciparum as comparison.

Results

We performed cloud and local computing to extract parasite transcriptomes from publicly available raw data of human blood samples. We developed a Poisson Modelling (PM) method to confidently identify parasite derived transcripts in mixed RNAseq signals of infected host tissues. We successfully retrieved and reconstructed parasite transcriptomes from infected patient blood as early as the first blood stage cycle; and the same methodology did not recover any significant signal from controls. Surprisingly, these first generation blood parasites already show strong signature of transmission, which indicates the commitment from asexual-to-sexual stages. Further, we place the results within the context of P. vivax’s complex life cycle, by developing mathematical models for P. vivax and P. falciparum and using sensitivity analysis assess the relative epidemiological impact of possible early stage transmission.

Conclusion

The study uncovers the earliest onset of P. vivax blood pathogenesis and highlights the challenges of P. vivax eradication programs.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0669-4) contains supplementary material, which is available to authorized users.

Highlights on the Application of Genomics and Bioinformatics in the Fight Against Infectious Diseases: Challenges and Opportunities in Africa

Genomics and bioinformatics are increasingly contributing to our understanding of infectious diseases caused by bacterial pathogens such as Mycobacterium tuberculosis and parasites such as Plasmodium falciparum. This ranges from investigations of disease outbreaks and pathogenesis, host and pathogen genomic variation, and host immune evasion mechanisms to identification of potential diagnostic markers and vaccine targets. High throughput genomics data generated from pathogens and animal models can be combined with host genomics and patients’ health records to give advice on treatment options as well as potential drug and vaccine interactions. However, despite accounting for the highest burden of infectious diseases, Africa has the lowest research output on infectious disease genomics. Here we review the contributions of genomics and bioinformatics to the management of infectious diseases of serious public health concern in Africa including tuberculosis (TB), dengue fever, malaria and filariasis. Furthermore, we discuss how genomics and bioinformatics can be applied to identify drug and vaccine targets. We conclude by identifying challenges to genomics research in Africa and highlighting how these can be overcome where possible.

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.

A synthetic promoter for multi-stage expression to probe complementary functions of Plasmodium adhesins

Gene expression of malaria parasites is mediated by the apicomplexan Apetala2 (ApiAP2) transcription factor family. Different ApiAP2s control gene expression at distinct stages in the complex life cycle of the parasite, ensuring timely expression of stage-specific genes. ApiAP2s recognize short cis-regulatory elements that are enriched in the upstream/promoter region of their target genes. This should, in principle, allow the generation of 'synthetic' promoters that drive gene expression at desired stages of the Plasmodium life cycle. Here we test this concept by combining cis-regulatory elements of two genes expressed successively within the mosquito part of the life cycle. Our tailored 'synthetic' promoters, named Spooki 1.0 and Spooki 2.0, activate gene expression in early and late mosquito stages, as shown by the expression of a fluorescent reporter. We used these promoters to address the specific functionality of two related adhesins that are exclusively expressed either during the early or late mosquito stage. By modifying the expression profile of both adhesins in absence of their counterpart we were able to test for complementary functions in gliding and invasion. We discuss the possible advantages and drawbacks of our approach.This article has an associated First Person interview with the first author of the paper.

In silico analysis of putative dormancy genes in Plasmodium vivax

Plasmodium vivax is the most widely spread species causing human malaria. The control of malaria caused by P. vivax has been largely hampered by its ability to develop a dormant liver stage that can generate a new blood infection at different periods of time. Unfortunately, the mechanisms of dormancy in P. vivax have not been thoroughly elucidated to date. In this study, the putative dormancy genes were analyzed to select genes with less genetic variability to maintain the function of relapsing. Expression data concerning these genes were searched to support the selection. Protein interactions among selected gene products were identified based on known and predicted protein-protein interaction using String database. Potentially interacting proteins (n = 15) were used to propose a mechanism involved in dormancy based on the differential vesicular transport due to the iron available in the hepatocyte.

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.

Whole blood transcriptome changes following controlled human malaria infection in malaria pre-exposed volunteers correlate with parasite prepatent period

Malaria continues to be one of mankind’s most devastating diseases despite the many and varied efforts to combat it. Indispensable for malaria elimination and eventual eradication is the development of effective vaccines. Controlled human malaria infection (CHMI) is an invaluable tool for vaccine efficacy assessment and investigation of early immunological and molecular responses against Plasmodium falciparum infection. Here, we investigated gene expression changes following CHMI using RNA-Seq. Peripheral blood samples were collected in Bagamoyo, Tanzania, from ten adults who were injected intradermally (ID) with 2.5x10⁴ aseptic, purified, cryopreserved P. falciparum sporozoites (Sanaria® PfSPZ Challenge). A total of 2,758 genes were identified as differentially expressed following CHMI. Transcriptional changes were most pronounced on day 5 after inoculation, during the clinically silent liver phase. A secondary analysis, grouping the volunteers according to their prepatent period duration, identified 265 genes whose expression levels were linked to time of blood stage parasitemia detection. Gene modules associated with these 265 genes were linked to regulation of transcription, cell cycle, phosphatidylinositol signaling and erythrocyte development. Our study showed that in malaria pre-exposed volunteers, parasite prepatent period in each individual is linked to magnitude and timing of early gene expression changes after ID CHMI.

Long non-coding RNAs as molecular players in plant defense against pathogens

Long non-coding RNAs (lncRNAs) has significant role in of gene expression and silencing pathways for several biological processes in eukaryotes. lncRNAs has been reported as key player in remodeling chromatin and genome architecture, RNA stabilization and transcription regulation, including enhancer-associated activity. Host lncRNAs are reckoned as compulsory elements of plant defense. In response to pathogen attack, plants protect themselves with the help of lncRNAs -dependent immune systems in which lncRNAs regulate pathogen-associated molecular patterns (PAMPs) and other effectors. Role of lncRNAs in plant microbe interaction has been studied extensively but regulations of several lncRNAs still need extensive research. In this study we discussed and provide as overview the topical advancements and findings relevant to pathogen attack and plant defense mediated by lncRNAs. It is hoped that lncRNAs would be exploited as a mainstream player to achieve food security by tackling different plant diseases.

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.

Advancing Research Models and Technologies to Overcome Biological Barriers to Plasmodium vivax Control

Malaria prevalence has declined in the past 10 years, especially outside of sub-Saharan Africa. However, the proportion of cases due to Plasmodium vivax is increasing, accounting for up to 90-100% of the malaria burden in endemic regions. Nonetheless, investments in malaria research and control still prioritize Plasmodium falciparum while largely neglecting P. vivax. Specific biological features of P. vivax, particularly invasion of reticulocytes, occurrence of dormant liver forms of the parasite, and the potential for transmission of sexual-stage parasites prior to onset of clinical illness, promote its persistence and hinder development of research tools and interventions. This review discusses recent advances in P. vivax research, current knowledge of its unique biology, and proposes priorities for P. vivax research and control efforts.