Transcriptome profiling reveals functional variation in Plasmodium falciparum parasites from controlled human malaria infection studies

Zombie malaria parasites

Natural infections are distinct from those of laboratory-or zombie-strains.

A novel computational pipeline for var gene expression augments the discovery of changes in the Plasmodium falciparum transcriptome during transition from in vivo to short-term in vitro culture

The pathogenesis of severe Plasmodium falciparum malaria involves cytoadhesive microvascular sequestration of infected erythrocytes, mediated by P. falciparum erythrocyte membrane protein 1 (PfEMP1). PfEMP1 variants are encoded by the highly polymorphic family of var genes, the sequences of which are largely unknown in clinical samples. Previously, we published new approaches for var gene profiling and classification of predicted binding phenotypes in clinical P. falciparum isolates (Wichers et al., 2021), which represented a major technical advance. Building on this, we report here a novel method for var gene assembly and multidimensional quantification from RNA-sequencing that outperforms the earlier approach of Wichers et al., 2021, on both laboratory and clinical isolates across a combination of metrics. Importantly, the tool can interrogate the var transcriptome in context with the rest of the transcriptome and can be applied to enhance our understanding of the role of var genes in malaria pathogenesis. We applied this new method to investigate changes in var gene expression through early transition of parasite isolates to in vitro culture, using paired sets of ex vivo samples from our previous study, cultured for up to three generations. In parallel, changes in non-polymorphic core gene expression were investigated. Modest but unpredictable var gene switching and convergence towards var2csa were observed in culture, along with differential expression of 19% of the core transcriptome between paired ex vivo and generation 1 samples. Our results cast doubt on the validity of the common practice of using short-term cultured parasites to make inferences about in vivo phenotype and behaviour.

Comparative Transcriptomics in Atlantic Salmon Head Kidney and SHK-1 Cell Line Exposed to the Sea Louse Cr-Cathepsin

The development of vaccines against sea lice in salmon farming is complex, expensive, and takes several years for commercial availability. Recently, transcriptome studies in sea louse have provided valuable information for identifying relevant molecules with potential use for fish vaccines. However, the bottleneck is the in vivo testing of recombinant protein candidates, the dosage, and the polyvalent formulation strategies. This study explored a cell-based approach to prospect antigens as candidate vaccines against sea lice by comparison with immunized fish. Herein, SHK-1 cells and Atlantic salmon head kidney tissue were exposed to the antigen cathepsin identified from the sea louse Caligus rogercresseyi. The cathepsin protein was cloned and recombinantly expressed in Escherichia coli, and then SHK-1 cell lines were stimulated with 100 ng/mL cathepsin recombinant for 24 h. In addition, Atlantic salmons were vaccinated with 30 ug/mL recombinant protein, and head kidney samples were then collected 30 days post-immunization. SHK-1 cells and salmon head kidney exposed to cathepsin were analyzed by Illumina RNA sequencing. The statistical comparisons showed differences in the transcriptomic profiles between SHK-1 cells and the salmon head kidney. However, 24.15% of the differentially expressed genes were shared. Moreover, putative gene regulation through lncRNAs revealed tissue-specific transcription patterns. The top 50 up and downregulated lncRNAs were highly correlated with genes involved in immune response, iron homeostasis, pro-inflammatory cytokines, and apoptosis. Also, highly enriched pathways related to the immune system and signal transduction were shared between both tissues. These findings highlight a novel approach to evaluating candidate antigens for sea lice vaccine development, improving the antigens screening in the SHK-1 cell line model.

Controlled Human Infection Models To Accelerate Vaccine Development

The timelines for developing vaccines against infectious diseases are lengthy, and often vaccines that reach the stage of large phase 3 field trials fail to provide the desired level of protective efficacy. The application of controlled human challenge models of infection and disease at the appropriate stages of development could accelerate development of candidate vaccines and, in fact, has done so successfully in some limited cases. Human challenge models could potentially be used to gather critical information on pathogenesis, inform strain selection for vaccines, explore cross-protective immunity, identify immune correlates of protection and mechanisms of protection induced by infection or evoked by candidate vaccines, guide decisions on appropriate trial endpoints, and evaluate vaccine efficacy. We prepared this report to motivate fellow scientists to exploit the potential capacity of controlled human challenge experiments to advance vaccine development. In this review, we considered available challenge models for 17 infectious diseases in the context of the public health importance of each disease, the diversity and pathogenesis of the causative organisms, the vaccine candidates under development, and each model's capacity to evaluate them and identify correlates of protective immunity. Our broad assessment indicated that human challenge models have not yet reached their full potential to support the development of vaccines against infectious diseases. On the basis of our review, however, we believe that describing an ideal challenge model is possible, as is further developing existing and future challenge models.

Pervasive sequence-level variation in the transcriptome of Plasmodium falciparum

Single-nucleotide variations (SNVs) in RNA, arising from co- and post-transcriptional phenomena including transcription errors and RNA-editing, are well studied in a range of organisms. In the malaria parasite Plasmodium falciparum, stage-specific and non-specific gene-expression variations accompany the parasite's array of developmental and morphological phenotypes over the course of its complex life cycle. However, the extent, rate and effect of sequence-level variation in the parasite's transcriptome are unknown. Here, we report the presence of pervasive, non-specific SNVs in the P. falciparum transcriptome. SNV rates for a gene were correlated to gene length (r[Formula: see text]0.65-0.7) but not to the AT-content of that gene. Global SNV rates for the P. falciparum lines we used, and for publicly available P. vivax and P. falciparum clinical isolate datasets, were of the order of 10-3 per base, ∼10× higher than rates we calculated for bacterial datasets. These variations may reflect an intrinsic transcriptional error rate in the parasite, and RNA editing may be responsible for a subset of them. This seemingly characteristic property of the parasite may have implications for clinical outcomes and the basic biology and evolution of P. falciparum and parasite biology more broadly. We anticipate that our study will prompt further investigations into the exact sources, consequences and possible adaptive roles of these SNVs.

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.

Dramatic transcriptomic differences in Macaca mulatta and Macaca fascicularis with Plasmodium knowlesi infections

Plasmodium knowlesi, a model malaria parasite, is responsible for a significant portion of zoonotic malaria cases in Southeast Asia and must be controlled to avoid disease severity and fatalities. However, little is known about the host-parasite interactions and molecular mechanisms in play during the course of P. knowlesi malaria infections, which also may be relevant across Plasmodium species. Here we contrast P. knowlesi sporozoite-initiated infections in Macaca mulatta and Macaca fascicularis using whole blood RNA-sequencing and transcriptomic analysis. These macaque hosts are evolutionarily close, yet malaria-naïve M. mulatta will succumb to blood-stage infection without treatment, whereas malaria-naïve M. fascicularis controls parasitemia without treatment. This comparative analysis reveals transcriptomic differences as early as the liver phase of infection, in the form of signaling pathways that are activated in M. fascicularis, but not M. mulatta. Additionally, while most immune responses are initially similar during the acute stage of the blood infection, significant differences arise subsequently. The observed differences point to prolonged inflammation and anti-inflammatory effects of IL10 in M. mulatta, while M. fascicularis undergoes a transcriptional makeover towards cell proliferation, consistent with its recovery. Together, these findings suggest that timely detection of P. knowlesi in M. fascicularis, coupled with control of inflammation while initiating the replenishment of key cell populations, helps contain the infection. Overall, this study points to specific genes and pathways that could be investigated as a basis for new drug targets that support recovery from acute malaria.

Expression Patterns of Plasmodium falciparum Clonally Variant Genes at the Onset of a Blood Infection in Malaria-Naive Humans

Clonally variant genes (CVGs) play fundamental roles in the adaptation of Plasmodium falciparum to fluctuating conditions of the human host. However, their expression patterns under the natural conditions of the blood circulation have been characterized in detail for only a few specific gene families. Here, we provide a detailed characterization of the complete P. falciparum transcriptome across the full intraerythrocytic development cycle (IDC) at the onset of a blood infection in malaria-naive human volunteers. We found that the vast majority of transcriptional differences between parasites obtained from the volunteers and the parental parasite line maintained in culture occurred in CVGs. In particular, we observed a major increase in the transcript levels of most genes of the pfmc-2tm and gbp families and of specific genes of other families, such as phist, hyp10, rif, or stevor, in addition to previously reported changes in var and clag3 gene expression. Increased transcript levels of individual pfmc-2tm, rif, and stevor genes involved activation in small subsets of parasites. Large transcriptional differences correlated with changes in the distribution of heterochromatin, confirming their epigenetic nature. Furthermore, the similar expression of several CVGs between parasites collected at different time points along the blood infection suggests that the epigenetic memory for multiple CVG families is lost during transmission stages, resulting in a reset of their transcriptional state. Finally, the CVG expression patterns observed in a volunteer likely infected by a single sporozoite suggest that new epigenetic patterns are established during liver stages.

The transcriptome of circulating sexually committed Plasmodium falciparum ring stage parasites forecasts malaria transmission potential

Malaria is spread by the transmission of sexual stage parasites, called gametocytes. However, with Plasmodium falciparum, gametocytes can only be detected in peripheral blood when they are mature and transmissible to a mosquito, which complicates control efforts. Here, we identify the set of genes overexpressed in patient blood samples with high levels of gametocyte-committed ring stage parasites. Expression of all 18 genes is regulated by transcription factor AP2-G, which is required for gametocytogenesis. We select three genes, not expressed in mature gametocytes, to develop as biomarkers. All three biomarkers we validate in vitro using 6 different parasite lines and develop an algorithm that predicts gametocyte production in ex vivo samples and volunteer infection studies. The biomarkers are also sensitive enough to monitor gametocyte production in asymptomatic P. falciparum carriers allowing early detection and treatment of infectious reservoirs, as well as the in vivo analysis of factors that modulate sexual conversion.

Malaria gametocytes are sexual-stage parasites transmitted from mammalian host’s blood back to their insect vector. Here, Prajapati et al. identify gametocyte-committed ring-stage biomarkers allowing to forecast malaria transmission potential.

A comprehensive RNA handling and transcriptomics guide for high-throughput processing of Plasmodium blood-stage samples

Background

Sequencing technology advancements opened new opportunities to use transcriptomics for studying malaria pathology and epidemiology. Even though in recent years the study of whole parasite transcriptome proved to be essential in understanding parasite biology there is no compiled up-to-date reference protocol for the efficient generation of transcriptome data from growing number of samples. Here, a comprehensive methodology on how to preserve, extract, amplify, and sequence full-length mRNA transcripts from Plasmodium-infected blood samples is presented that can be fully streamlined for high-throughput studies.

Results

The utility of various commercially available RNA-preserving reagents in a range of storage conditions was evaluated. Similarly, several RNA extraction protocols were compared and the one most suitable method for the extraction of high-quality total RNA from low-parasitaemia and low-volume blood samples was established. Furthermore, the criteria needed to evaluate the quality and integrity of Plasmodium RNA in the presence of human RNA was updated. Optimization of SMART-seq2 amplification method to better suit AT-rich Plasmodium falciparum RNA samples allowed us to generate high-quality transcriptomes from as little as 10 ng of total RNA and a lower parasitaemia limit of 0.05%. Finally, a modified method for depletion of unwanted human haemoglobin transcripts using in vitro CRISPR-Cas9 treatment was designed, thus improving parasite transcriptome coverage in low parasitaemia samples. To prove the functionality of the pipeline for both laboratory and field strains, the highest  2-hour resolution RNA-seq transcriptome for P. falciparum 3D7 intraerythrocytic life cycle available to  date was generated, and the entire protocol was applied to create the largest transcriptome data from Southeast Asian field isolates.

Conclusions

Overall, the presented methodology is an inclusive pipeline for generation of good quality transcriptomic data from a diverse range of Plasmodium-infected blood samples with varying parasitaemia and RNA inputs. The flexibility of this pipeline to be adapted to robotic handling will facilitate both small and large-scale future transcriptomic studies in the field of malaria.

Refining the transcriptome of the human malaria parasite Plasmodium falciparum using amplification-free RNA-seq

Background

Plasmodium parasites undergo several major developmental transitions during their complex lifecycle, which are enabled by precisely ordered gene expression programs. Transcriptomes from the 48-h blood stages of the major human malaria parasite Plasmodium falciparum have been described using cDNA microarrays and RNA-seq, but these assays have not always performed well within non-coding regions, where the AT-content is often 90–95%.

Results

We developed a directional, amplification-free RNA-seq protocol (DAFT-seq) to reduce bias against AT-rich cDNA, which we have applied to three strains of P. falciparum (3D7, HB3 and IT). While strain-specific differences were detected, overall there is strong conservation between the transcriptional profiles. For the 3D7 reference strain, transcription was detected from 89% of the genome, with over 78% of the genome transcribed into mRNAs. We also find that transcription from bidirectional promoters frequently results in non-coding, antisense transcripts. These datasets allowed us to refine the 5′ and 3′ untranslated regions (UTRs), which can be variable, long (> 1000 nt), and often overlap those of adjacent transcripts.

Conclusions

The approaches applied in this study allow a refined description of the transcriptional landscape of P. falciparum and demonstrate that very little of the densely packed P. falciparum genome is inactive or redundant. By capturing the 5′ and 3′ ends of mRNAs, we reveal both constant and dynamic use of transcriptional start sites across the intraerythrocytic developmental cycle that will be useful in guiding the definition of regulatory regions for use in future experimental gene expression studies.

Identification of a conserved var gene in different Plasmodium falciparum strains

BACKGROUND

The multicopy var gene family of Plasmodium falciparum is of crucial importance for pathogenesis and antigenic variation. So far only var2csa, the var gene responsible for placental malaria, was found to be highly conserved among all P. falciparum strains. Here, a new conserved 3D7 var gene (PF3D7_0617400) is identified in several field isolates.

METHODS

DNA sequencing, transcriptional analysis, Cluster of Differentiation (CD) 36-receptor binding, indirect immunofluorescence with PF3D7_0617400-antibodies and quantification of surface reactivity against semi-immune sera were used to characterize an NF54 clone and a Gabonese field isolate clone (MOA C3) transcribing the gene. A population of 714 whole genome sequenced parasites was analysed to characterize the conservation of the locus in African and Asian isolates. The genetic diversity of two var2csa fragments was compared with the genetic diversity of 57 microsatellites fragments in field isolates.

RESULTS

PFGA01_060022400 was identified in a Gabonese parasite isolate (MOA) from a chronic infection and found to be 99% identical with PF3D7_0617400 of the 3D7 genome strain. Transcriptional analysis and immunofluorescence showed expression of the gene in an NF54 and a MOA clone but CD36 binding assays and surface reactivity to semi-immune sera differed markedly in the two clones. Long-read Pacific bioscience whole genome sequencing showed that PFGA01_060022400 is located in the internal cluster of chromosome 6. The full length PFGA01_060022400 was detected in 36 of 714 P. falciparum isolates and 500 bp fragments were identified in more than 100 isolates. var2csa was in parts highly conserved (H_e = 0) but in other parts as variable (H_e = 0.86) as the 57 microsatellites markers (H_e = 0.8).

CONCLUSIONS

Individual var gene sequences exhibit conservation in the global parasite population suggesting that purifying selection may limit overall genetic diversity of some var genes. Notably, field and laboratory isolates expressing the same var gene exhibit markedly different phenotypes.

Human unconventional T cells in Plasmodium falciparum infection

Malaria is an old scourge of humankind and has a large negative impact on the economic development of affected communities. Recent success in malaria control and reduction of mortality seems to have stalled emphasizing that our current intervention tools need to be complemented by malaria vaccines. Different populations of unconventional T cells such as mucosal-associated invariant T (MAIT) cells, invariant natural killer T (iNKT) cells and γδ T cells are gaining attention in the field of malaria immunology. Significant advances in our basic understanding of unconventional T cell biology in rodent malaria models have been made, however, their roles in humans during malaria are less clear. Unconventional T cells are abundant in skin, gut and liver tissues, and long-lasting expansions and functional alterations were observed upon malaria infection in malaria naïve and malaria pre-exposed volunteers. Here, we review the current understanding of involvement of unconventional T cells in anti-Plasmodium falciparum immunity and highlight potential future research avenues.