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

A virus associated with the zoonotic pathogen Plasmodium knowlesi causing human malaria is a member of a diverse and unclassified viral taxon

Apicomplexa are single-celled eukaryotes that can infect humans and include the mosquito-borne parasite Plasmodium, the cause of malaria. Increasing rates of drug resistance in human-only Plasmodium species are reducing the efficacy of control efforts and antimalarial treatments. There are also rising cases of P. knowlesi, the only zoonotic Plasmodium species that causes severe disease and death in humans. Thus, there is a need to develop additional innovative strategies to combat malaria. Viruses that infect non-Plasmodium spp. Disease-causing protozoa have been shown to affect pathogen life cycle and disease outcomes. However, only one virus (Matryoshka RNA virus 1) has been identified in Plasmodium, and none have been identified in zoonotic Plasmodium species. The rapid expansion of the known RNA virosphere using structure- and artificial intelligence-based methods suggests that this dearth is due to the divergent nature of RNA viruses that infect protozoa. We leveraged these newly uncovered data sets to explore the virome of human-infecting Plasmodium species collected in Sabah, east (Borneo) Malaysia. We identified a highly divergent RNA virus in two human-infecting P. knowlesi isolates that is related to the unclassified group 'ormycoviruses'. By characterising fifteen additional ormycoviruses identified in the transcriptomes of arthropods we show that this group of viruses exhibits a complex ecology at the arthropod-mammal interface. Through the application of artificial intelligence methods, we then demonstrate that the ormycoviruses are part of a diverse and unclassified viral taxon. This is the first observation of an RNA virus in a zoonotic Plasmodium species. By linking small-scale experimental data to large-scale virus discovery advances, we characterise the diversity and genomic architecture of an unclassified viral taxon. This approach should be used to further explore the virome of disease-causing Apicomplexa and better understand how protozoa-infecting viruses may affect parasite fitness, pathobiology, and treatment outcomes.

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

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

Gene expression analyses reveal differences in children's response to malaria according to their age

In Bandiagara, Mali, children experience on average two clinical malaria episodes per year. However, even in the same transmission area, the number of uncomplicated symptomatic infections, and their parasitemia, can vary dramatically among children. We simultaneously characterize host and parasite gene expression profiles from 136 Malian children with symptomatic falciparum malaria and examine differences in the relative proportion of immune cells and parasite stages, as well as in gene expression, associated with infection and or patient characteristics. Parasitemia explains much of the variation in host and parasite gene expression, and infections with higher parasitemia display proportionally more neutrophils and fewer T cells, suggesting parasitemia-dependent neutrophil recruitment and/or T cell extravasation to secondary lymphoid organs. The child's age also strongly correlates with variations in gene expression: Plasmodium falciparum genes associated with age suggest that older children carry more male gametocytes, while variations in host gene expression indicate a stronger innate response in younger children and stronger adaptive response in older children. These analyses highlight the variability in host responses and parasite regulation during P. falciparum symptomatic infections and emphasize the importance of considering the children's age when studying and treating malaria infections.

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

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

In-vitro susceptibility and ex-vivo evaluation of macrocyclic lactone endectocides sub-lethal concentrations against Plasmodium vivax oocyst development in Anopheles arabiensis

BACKGROUND

Asymptomatic malaria transmission has become a public health concern across malaria-endemic Africa including Ethiopia. Specifically, Plasmodium vivax is more efficient at transmitting earlier in the infection and at lower densities than Plasmodium falciparum. Consequently, a greater proportion of individuals infected with P. vivax can transmit without detectable gametocytaemia. Mass treatment of livestock with macrocyclic lactones (MLs), e.g., ivermectin and doramectin, was suggested as a complementary malaria vector tool because of their insecticidal effects. However, the effects of MLs on P. vivax in Anopheles arabiensis has not yet been fully explored. Hence, comparative in-vitro susceptibility and ex-vivo studies were conducted to evaluate the effects of ivermectin, doramectin and moxidectin sub-lethal concentrations on P. vivax oocyst development in An. arabiensis.

METHODS

The 7-day sub-lethal concentrations of 25% (LC25) and 5% (LC5) were determined from in-vitro susceptibility tests on female An. arabiensis in Hemotek® membrane feeding assay. Next, an ex-vivo study was conducted using P. vivax gametocytes infected patient's blood spiked with the LC25 and LC5 of the MLs. At 7-days post-feeding, each mosquito was dissected under a dissection stereo microscope, stained with 0.5% (w/v) mercurochrome solution, and examined for the presence of P. vivax oocysts. Statistical analysis was based on a generalized mixed model with binomially distributed error terms.

RESULTS

A 7-day lethal concentration of 25% (LC25, in ng/mL) of 7.1 (95% CI: [6.3;8.0]), 20.0 (95%CI:[17.8;22.5]) and 794.3 (95%CI:[716.4;1516.3]) were obtained for ivermectin, doramectin and moxidectin, respectively. Similarly, a lethal concentration of 5% (LC5, in ng/mL) of 0.6 (95% CI: [0.5;0.7]), 1.8 (95% CI:[1.6;2.0]) and 53.7 (95% CI:[ 48.4;102.5]) were obtained respectively for ivermectin, doramectin and moxidectin. The oocyst prevalence in treatment and control groups did not differ significantly (p > 0.05) from each other. Therefore, no direct effect of ML endectocides on P. vivax infection in An. arabiensis mosquitoes was observed at the sub-lethal concentration (LC25 and LC5).

CONCLUSIONS

The effects of ivermectin and doramectin on malaria parasite is more likely via indirect effects, particularly by reducing the vectors lifespan and causing mortality before completing the parasite's sporogony cycle or reducing their vector capacity as it affects the locomotor activity of the mosquito.

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.

Gene expression analyses reveal differences in children's response to malaria according to their age

In Bandiagara, Mali, children experience on average two clinical malaria episodes per season. However, even in the same transmission area, the number of uncomplicated symptomatic infections, and their parasitemia, vary dramatically among children. To examine the factors contributing to these variations, we simultaneously characterized the host and parasite gene expression profiles from 136 children with symptomatic falciparum malaria and analyzed the expression of 9,205 human and 2,484 Plasmodium genes. We used gene expression deconvolution to estimate the relative proportion of immune cells and parasite stages in each sample and to adjust the differential gene expression analyses. Parasitemia explained much of the variation in both host and parasite gene expression and revealed that infections with higher parasitemia had more neutrophils and fewer T cells, suggesting parasitemia-dependent neutrophil recruitment and/or T cell extravasation to secondary lymphoid organs. The child's age was also strongly correlated with gene expression variations. Plasmodium falciparum genes associated with age suggested that older children carried more male gametocytes, while host genes associated with age indicated a stronger innate response (through TLR and NLR signaling) in younger children and stronger adaptive immunity (through TCR and BCR signaling) in older children. These analyses highlight the variability in host responses and parasite regulation during P. falciparum symptomatic infections and emphasize the importance of considering the children's age when studying and treating malaria infections.

Gene expression analyses reveal differences in children's response to malaria according to their age

In Bandiagara, Mali, children experience on average two clinical malaria episodes per season. However, even in the same transmission area, the number of uncomplicated symptomatic infections, and their parasitemia, vary dramatically among children. To examine the factors contributing to these variations, we simultaneously characterized the host and parasite gene expression profiles from 136 children with symptomatic falciparum malaria and analyzed the expression of 9,205 human and 2,484 Plasmodium genes. We used gene expression deconvolution to estimate the relative proportion of immune cells and parasite stages in each sample and to adjust the differential gene expression analyses. Parasitemia explained much of the variation in both host and parasite gene expression and revealed that infections with higher parasitemia had more neutrophils and fewer T cells, suggesting parasitemia-dependent neutrophil recruitment and/or T cell extravasation to secondary lymphoid organs. The child's age was also strongly correlated with gene expression variations. Plasmodium falciparum genes associated with age suggested that older children carried more male gametocytes, while host genes associated with age indicated a stronger innate response (through TLR and NLR signaling) in younger children and stronger adaptive immunity (through TCR and BCR signaling) in older children. These analyses highlight the variability in host responses and parasite regulation during P. falciparum symptomatic infections and emphasize the importance of considering the children's age when studying and treating malaria infections.

Plasmodium vivax transcriptomics: What is new?

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

Detection of Plasmodium simium gametocytes in non-human primates from the Brazilian Atlantic Forest

BACKGROUND

Plasmodium species of non-human primates (NHP) are of great interest because they can naturally infect humans. Plasmodium simium, a parasite restricted to the Brazilian Atlantic Forest, was recently shown to cause a zoonotic outbreak in the state of Rio de Janeiro. The potential of NHP to act as reservoirs of Plasmodium infection presents a challenge for malaria elimination, as NHP will contribute to the persistence of the parasite. The aim of the current study was to identify and quantify gametocytes in NHP naturally-infected by P. simium.

METHODS

Whole blood samples from 35 NHP were used in quantitative reverse transcription PCR (RT-qPCR) assays targeting 18S rRNA, Pss25 and Pss48/45 malaria parasite transcripts. Absolute quantification was performed in positive samples for 18S rRNA and Pss25 targets. Linear regression was used to compare the quantification cycle (Cq) and the Spearman's rank correlation coefficient was used to assess the correlation between the copy numbers of 18S rRNA and Pss25 transcripts. The number of gametocytes/µL was calculated by applying a conversion factor of 4.17 Pss25 transcript copies per gametocyte.

RESULTS

Overall, 87.5% of the 26 samples, previously diagnosed as P. simium, were positive for 18S rRNA transcript amplification, of which 13 samples (62%) were positive for Pss25 transcript amplification and 7 samples (54%) were also positive for Pss48/45 transcript. A strong positive correlation was identified between the Cq of the 18S rRNA and Pss25 and between the Pss25 and Pss48/45 transcripts. The 18S rRNA and Pss25 transcripts had an average of 1665.88 and 3.07 copies/µL, respectively. A positive correlation was observed between the copy number of Pss25 and 18S rRNA transcripts. Almost all gametocyte carriers exhibited low numbers of gametocytes (< 1/µL), with only one howler monkey having 5.8 gametocytes/µL.

CONCLUSIONS

For the first time, a molecular detection of P. simium gametocytes in the blood of naturally-infected brown howler monkeys (Alouatta guariba clamitans) was reported here, providing evidence that they are likely to be infectious and transmit P. simium infection, and, therefore, may act as a reservoir of malaria infection for humans in the Brazilian Atlantic Forest.

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.

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.

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

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

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.

Gene Polymorphisms Among Plasmodium vivax Geographical Isolates and the Potential as New Biomarkers for Gametocyte Detection

The unique biological features of Plasmodium vivax not only make it difficult to control but also to eliminate. For the transmission of the malaria parasite from infected human to the vector, gametocytes play a major role. The transmission potential of a malarial infection is inferred based on microscopic detection of gametocytes and molecular screening of genes in the female gametocytes. Microscopy-based detection methods could grossly underestimate the reservoirs of infection as gametocytes may occur as submicroscopic or as micro- or macro-gametocytes. The identification of genes that are highly expressed and polymorphic in male and female gametocytes is critical for monitoring changes not only in their relative proportions but also the composition of gametocyte clones contributing to transmission over time. Recent transcriptomic study revealed two distinct clusters of highly correlated genes expressed in the P. vivax gametocytes, indicating that the male and female terminal gametocytogeneses are independently regulated. However, the detective power of these genes is unclear. In this study, we compared genetic variations of 15 and 11 genes expressed, respectively, in the female and male gametocytes among P. vivax isolates from Southeast Asia, Africa, and South America. Further, we constructed phylogenetic trees to determine the resolution power and clustering patterns of gametocyte clones. As expected, Pvs25 (PVP01_0616100) and Pvs16 (PVP01_0305600) expressed in the female gametocytes were highly conserved in all geographical isolates. In contrast, genes including 6-cysteine protein Pvs230 (PVP01_0415800) and upregulated in late gametocytes ULG8 (PVP01_1452800) expressed in the female gametocytes, as well as two CPW-WPC family proteins (PVP01_1215900 and PVP01_1320100) expressed in the male gametocytes indicated considerably high nucleotide and haplotype diversity among isolates. Parasite samples expressed in male and female gametocyte genes were observed in separate phylogenetic clusters and likely represented distinct gametocyte clones. Compared to Pvs25, Pvs230 (PVP01_0415800) and a CPW-WPC family protein (PVP01_0904300) showed higher expression in a subset of Ethiopian P. vivax samples. Thus, Pvs230, ULG8, and CPW-WPC family proteins including PVP01_0904300, PVP01_1215900, and PVP01_1320100 could potentially be used as novel biomarkers for detecting both sexes of P. vivax gametocytes in low-density infections and estimating transmission reservoirs.

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.

Plasmodium vivax gametocytes and transmission

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

High Proportion of Genome-Wide Homology and Increased Pretreatment pvcrt Levels in Plasmodium vivax Late Recurrences: a Chloroquine Therapeutic Efficacy Study

Chloroquine (CQ) is the first-line treatment for Plasmodium vivax malaria in most countries where malaria is endemic. Monitoring P. vivax CQ resistance (CQR) is critical but remains challenged by the difficulty to distinguish real treatment failure from reinfection or liver relapse. The therapeutic efficacy of CQ against uncomplicated P. vivax malaria was evaluated in Gia Lai Province, Vietnam. Sixty-seven patients were enrolled and followed for 42 days using microscopy and quantitative PCR. Adequate clinical and parasitological response (ACPR) was 100% (66/66) on day 28 but 75.4% (49/65) on day 42. Eighteen recurrences (27.7%) were detected, with a median time to recurrence of 42 days (interquartile range [IQR], 35 to 42) and blood CQ concentration of <100 ng/ml. Primary infections leading to recurrence occurred in younger individuals (median age for ACPR = 25 years [IQR, 20 to 28]; recurrences = 18 [16 to 21]; P = 0.002) had a longer parasite clearance time (PCT for ACPR = 47.5 h [IQR, 36.2 to 59.8 h]; recurrences = 54.2 [48.4 to 62.0]; P = 0.035) and higher pvcrt gene expression (median relative expression ratio for ACPR = 0.09 [IQR, 0.05 to 0.22]; recurrences = 0.20 [0.15 to 0.56]; P = 0.002), but showed no differences in ex vivo CQ sensitivity. Parasite genotyping by microsatellites, single nucleotide polymorphism (SNP) barcoding, and whole-genome sequencing (WGS) identified a majority of homologous recurrences, with 80% (8/10) showing >98% identity by descent to paired day 0 samples. This study shows that CQ remained largely efficacious to treat P. vivax in Gia Lai; i.e., recurrences occurred late (>day 28) and in the presence of low blood CQ concentrations. However, the combination of both WGS and gene expression analysis (pvcrt) data with clinical data (PCT) allowed us to identify potential emergence of low-grade CQR, which should be closely monitored. (This study has been registered at ClinicalTrials.gov under identifier NCT02610686.).

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.

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.

The host targeting effect of chloroquine in malaria

Due to the rapid onset and spread of the COVID-19 pandemic, the treatment of COVID-19 patients by hydroxychloroquine alone or in combination with other drugs has captured a great deal of attention and triggered considerable debate. Historically, the worldwide use of quinoline based-drugs has led to a spectacular reduction in death from malaria. Unfortunately, scientists have been forced to seek alternative drugs to treat malaria due to the emergence of chloroquine-resistant parasites in the 1960s. The repurposing of hydroxychloroquine against viral infections, various types of cancer and autoimmune diseases has been ongoing for more than 70 years, with no clear understanding of its mechanism of action (MOA). Here, we closely examine the MOA of this old but influential drug in and beyond malaria. Better insights into how chloroquine targets the host's cellular and immune responses may help to develop applications against to new pathogens and diseases, and perhaps even restore the clinical utility of chloroquine against malaria.

Host and Parasite Transcriptomic Changes upon Successive Plasmodium falciparum Infections in Early Childhood

We show that dual RNA-seq from patient blood samples allows characterization of host/parasite interactions during malaria infections and can provide a solid framework to study the acquisition of antimalarial immunity, as well as the adaptations of P. falciparum to malaria-experienced hosts.

Children are highly susceptible to clinical malaria, and in regions where malaria is endemic, their immune systems must face successive encounters with Plasmodium falciparum parasites before they develop immunity, first against severe disease and later against uncomplicated malaria. Understanding cellular and molecular interactions between host and parasites during an infection could provide insights into the processes underlying this gradual acquisition of immunity, as well as to how parasites adapt to infect hosts that are successively more malaria experienced. Here, we describe methods to analyze the host and parasite gene expression profiles generated simultaneously from blood samples collected from five consecutive symptomatic P. falciparum infections in three Malian children. We show that the data generated enable statistical assessment of the proportions of (i) each white blood cell subset and (ii) the parasite developmental stages, as well as investigations of host-parasite gene coexpression. We also use the sequences generated to analyze allelic variations in transcribed regions and determine the complexity of each infection. While limited by the modest sample size, our analyses suggest that host gene expression profiles primarily clustered by individual, while the parasite gene expression profiles seemed to differentiate early from late infections. Overall, this study provides a solid framework to examine the mechanisms underlying acquisition of immunity to malaria infections using whole-blood transcriptome sequencing (RNA-seq).

IMPORTANCE We show that dual RNA-seq from patient blood samples allows characterization of host/parasite interactions during malaria infections and can provide a solid framework to study the acquisition of antimalarial immunity, as well as the adaptations of P. falciparum to malaria-experienced hosts.

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.

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.

Immunity against sexual stage Plasmodium falciparum and Plasmodium vivax parasites

The efficient spread of malaria from infected humans to mosquitoes is a major challenge for malaria elimination initiatives. Gametocytes are the only Plasmodium life stage infectious to mosquitoes. Here, we summarize evidence for naturally acquired anti-gametocyte immunity and the current state of transmission blocking vaccines (TBV). Although gametocytes are intra-erythrocytic when present in infected humans, developing Plasmodium falciparum gametocytes may express proteins on the surface of red blood cells that elicit immune responses in naturally exposed individuals. This immune response may reduce the burden of circulating gametocytes. For both P. falciparum and Plasmodium vivax, there is a solid evidence that antibodies against antigens present on the gametocyte surface, when co-ingested with gametocytes, can influence transmission to mosquitoes. Transmission reducing immunity, reducing the burden of infection in mosquitoes, is a well-acknowledged but poorly quantified phenomenon that forms the basis for the development of TBV. Transmission enhancing immunity, increasing the likelihood or intensity of transmission to mosquitoes, is more speculative in nature but is convincingly demonstrated for P. vivax. With the increased interest in malaria elimination, TBV and monoclonal antibodies have moved to the center stage of malaria vaccine development. Methodologies to prioritize and evaluate products are urgently needed.

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

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