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

Plasmodium vivax spleen-dependent protein 1 and its role in extracellular vesicles-mediated intrasplenic infections

Recent studies indicate that human spleen contains over 95% of the total parasite biomass during chronic asymptomatic infections caused by Plasmodium vivax. Previous studies have demonstrated that extracellular vesicles (EVs) secreted from infected reticulocytes facilitate binding to human spleen fibroblasts (hSFs) and identified parasite genes whose expression was dependent on an intact spleen. Here, we characterize the P. vivax spleen-dependent hypothetical gene (PVX_114580). Using CRISPR/Cas9, PVX_114580 was integrated into P. falciparum 3D7 genome and expressed during asexual stages. Immunofluorescence analysis demonstrated that the protein, which we named P. vivax Spleen-Dependent Protein 1 (PvSDP1), was located at the surface of infected red blood cells in the transgenic line and this localization was later confirmed in natural infections. Plasma-derived EVs from P. vivax-infected individuals (PvEVs) significantly increased cytoadherence of 3D7_PvSDP1 transgenic line to hSFs and this binding was inhibited by anti-PvSDP1 antibodies. Single-cell RNAseq of PvEVs-treated hSFs revealed increased expression of adhesion-related genes. These findings demonstrate the importance of parasite spleen-dependent genes and EVs from natural infections in the formation of intrasplenic niches in P. vivax, a major challenge for malaria elimination.

A single cell atlas of sexual development in Plasmodium falciparum

The developmental decision made by malaria parasites to become sexual underlies all malaria transmission. Here, we describe a rich atlas of short- and long-read single-cell transcriptomes of over 37,000 Plasmodium falciparum cells across intraerythrocytic asexual and sexual development. We used the atlas to explore transcriptional modules and exon usage along sexual development and expanded it to include malaria parasites collected from four Malian individuals naturally infected with multiple P. falciparum strains. We investigated genotypic and transcriptional heterogeneity within and among these wild strains at the single-cell level, finding differential expression between different strains even within the same host. These data are a key addition to the Malaria Cell Atlas interactive data resource, enabling a deeper understanding of the biology and diversity of transmission stages.

Understanding the complex formation of falstatin; an endogenous macromolecular inhibitor of falcipains

Proteolytic activity constitutes a fundamental process essential for the survival of the malaria parasite and is thus highly regulated. Falstatin, a protease inhibitor of Plasmodium falciparum, tightly regulates the activity of cysteine hemoglobinases, falcipain-2 and 3 (FP2, FP3), by inhibiting FP2 through a single surface exposed loop. However, the multimeric nature of falstatin and its interaction with FP2 remained unexplored. Here we report that the N-terminal falstatin region is highly disordered, and needs chaperone activity (heat-shock protein 70, HSP70) for its folding. Protein-protein interaction assays showed a significant interaction between falstatin and HSP70. Further, characterization of the falstatin multimer through a series of biophysical techniques identified the formation of a falstatin decamer, which was extremely thermostable. Computational analysis of the falstatin decamer showed the presence of five falstatin dimers, with each dimer aligned in a head-to-tail orientation. Further, the falstatin C-terminal region was revealed to be primarily involved in the oligomerization process. Stoichiometric analysis of the FP2-falstatin multimer showed the formation of a heterooligomeric complex in a 1:1 ratio, with the participation of ten subunits of each protein. Taken together, our results report a novel protease-inhibitor complex and strengthens our understanding of the regulatory mechanisms of major plasmodium hemoglobinases.

Identifying the best PCR enzyme for library amplification in NGS

Background. PCR amplification is a necessary step in many next-generation sequencing (NGS) library preparation methods [1, 2]. Whilst many PCR enzymes are developed to amplify single targets efficiently, accurately and with specificity, few are developed to meet the challenges imposed by NGS PCR, namely unbiased amplification of a wide range of different sizes and GC content. As a result PCR amplification during NGS library prep often results in bias toward GC neutral and smaller fragments. As NGS has matured, optimized NGS library prep kits and polymerase formulations have emerged and in this study we have tested a wide selection of available enzymes for both short-read Illumina library preparation and long fragment amplification ahead of long-read sequencing.We tested over 20 different hi-fidelity PCR enzymes/NGS amplification mixes on a range of Illumina library templates of varying GC content and composition, and find that both yield and genome coverage uniformity characteristics of the commercially available enzymes varied dramatically. Three enzymes Quantabio RepliQa Hifi Toughmix, Watchmaker Library Amplification Hot Start Master Mix (2X) 'Equinox' and Takara Ex Premier were found to give a consistent performance, over all genomes, that mirrored closely that observed for PCR-free datasets. We also test a range of enzymes for long-read sequencing by amplifying size fractionated S. cerevisiae DNA of average size 21.6 and 13.4 kb, respectively.The enzymes of choice for short-read (Illumina) library fragment amplification are Quantabio RepliQa Hifi Toughmix, Watchmaker Library Amplification Hot Start Master Mix (2X) 'Equinox' and Takara Ex Premier, with RepliQa also being the best performing enzyme from the enzymes tested for long fragment amplification prior to long-read sequencing.

Reading the epitranscriptome of the human malaria parasite

Epigenetic machinery has emerged as a central player in gene regulation and chromatin organization in Plasmodium spp. Epigenetic modifications on histones and their role in antigenic variation in P. falciparum are widely studied. Recent discoveries on nucleic acid methylome are exciting and provide a new dimension to the apicomplexan protozoan parasite's gene regulatory process. Reports have confirmed that N6-methyl adenosine (m6A) methylation plays a crucial role in the translational plasticity of the human malaria parasite during its development in RBC. The YTH domain (YT521-B Homology) protein in P. falciparum binds to m6A epitranscriptome modifications on the mRNA and regulates protein translation. The binding of the PfYTH domain protein to the m6A-modified mRNA is mediated through a binding pocket formed by aromatic amino acids. The P. falciparum genome encodes two members of YTH domain proteins, i.e., YTH1 and YTH2, and both have distinct roles in dictating the epitranscriptome in human malaria parasites. This review highlights recent advancements in the functions and mechanisms of YTH domain protein's role in translational plasticity in the various developmental stages of the parasite.

The genetic landscape of origins of replication in P. falciparum

Various origin mapping approaches have enabled genome-wide identification of origins of replication (ORI) in model organisms, but only a few studies have focused on divergent organisms. By employing three complementary approaches we provide a high-resolution map of ORIs in Plasmodium falciparum, the deadliest human malaria parasite. We profiled the distribution of origin of recognition complex (ORC) binding sites by ChIP-seq of two PfORC subunits and mapped active ORIs using NFS and SNS-seq. We show that ORIs lack sequence specificity but are not randomly distributed, and group in clusters. Licensing is biased towards regions of higher GC content and associated with G-quadruplex forming sequences (G4FS). While strong transcription likely enhances firing, active origins are depleted from transcription start sites. Instead, most accumulate in transcriptionally active gene bodies. Single molecule analysis of nanopore reads containing multiple initiation events, which could have only come from individual nuclei, showed a relationship between the replication fork pace and the distance to the nearest origin. While some similarities were drawn with the canonic eukaryote model, the distribution of ORIs in P. falciparum is likely shaped by unique genomic features such as extreme AT-richness-a product of evolutionary pressure imposed by the parasitic lifestyle.

Beyond the MEP Pathway: A novel kinase required for prenol utilization by malaria parasites

A proposed treatment for malaria is a combination of fosmidomycin and clindamycin. Both compounds inhibit the methylerythritol 4-phosphate (MEP) pathway, the parasitic source of farnesyl and geranylgeranyl pyrophosphate (FPP and GGPP, respectively). Both FPP and GGPP are crucial for the biosynthesis of several essential metabolites such as ubiquinone and dolichol, as well as for protein prenylation. Dietary prenols, such as farnesol (FOH) and geranylgeraniol (GGOH), can rescue parasites from MEP inhibitors, suggesting the existence of a missing pathway for prenol salvage via phosphorylation. In this study, we identified a gene in the genome of P. falciparum, encoding a transmembrane prenol kinase (PolK) involved in the salvage of FOH and GGOH. The enzyme was expressed in Saccharomyces cerevisiae, and its FOH/GGOH kinase activities were experimentally validated. Furthermore, conditional knockout parasites (Δ-PolK) were created to investigate the biological importance of the FOH/GGOH salvage pathway. Δ-PolK parasites were viable but displayed increased susceptibility to fosmidomycin. Their sensitivity to MEP inhibitors could not be rescued by adding prenols. Additionally, Δ-PolK parasites lost their capability to utilize prenols for protein prenylation. Experiments using culture medium supplemented with whole/delipidated human plasma in transgenic parasites revealed that human plasma has components that can diminish the effectiveness of fosmidomycin. Mass spectrometry tests indicated that both bovine supplements used in culture and human plasma contain GGOH. These findings suggest that the FOH/GGOH salvage pathway might offer an alternate source of isoprenoids for malaria parasites when de novo biosynthesis is inhibited. This study also identifies a novel kind of enzyme related to isoprenoid metabolism.

Promising antimalarial hits from phenotypic screens: a review of recently-described multi-stage actives and their modes of action

Over the last two decades, global malaria cases caused by Plasmodium falciparum have declined due to the implementation of effective treatments and the use of insecticides. However, the COVID-19 pandemic caused major disruption in the timely delivery of medical goods and diverted public health resources, impairing malaria control. The emergence of resistance to all existing frontline antimalarials underpins an urgent need for new antimalarials with novel mechanisms of action. Furthermore, the need to reduce malaria transmission and/or prevent malaria infection has shifted the focus of antimalarial research towards the discovery of compounds that act beyond the symptomatic blood stage and also impact other parasite life cycle stages. Phenotypic screening has been responsible for the majority of new antimalarial lead compounds discovered over the past 10 years. This review describes recently reported novel antimalarial hits that target multiple parasite stages and were discovered by phenotypic screening during the COVID-19 pandemic. Their modes of action and targets in blood stage parasites are also discussed.

Effects of the G-quadruplex-binding drugs quarfloxin and CX-5461 on the malaria parasite Plasmodium falciparum

Plasmodium falciparum is the deadliest causative agent of human malaria. This parasite has historically developed resistance to most drugs, including the current frontline treatments, so new therapeutic targets are needed. Our previous work on guanine quadruplexes (G4s) in the parasite's DNA and RNA has highlighted their influence on parasite biology, and revealed G4 stabilising compounds as promising candidates for repositioning. In particular, quarfloxin, a former anticancer agent, kills blood-stage parasites at all developmental stages, with fast rates of kill and nanomolar potency. Here we explored the molecular mechanism of quarfloxin and its related derivative CX-5461. In vitro, both compounds bound to P. falciparum-encoded G4 sequences. In cellulo, quarfloxin was more potent than CX-5461, and could prevent establishment of blood-stage malaria in vivo in a murine model. CX-5461 showed clear DNA damaging activity, as reported in human cells, while quarfloxin caused weaker signatures of DNA damage. Both compounds caused transcriptional dysregulation in the parasite, but the affected genes were largely different, again suggesting different modes of action. Therefore, CX-5461 may act primarily as a DNA damaging agent in both Plasmodium parasites and mammalian cells, whereas the complete antimalarial mode of action of quarfloxin may be parasite-specific and remains somewhat elusive.

Plasmodium falciparum GCN5 plays a key role in regulating artemisinin resistance-related stress responses

Plasmodium falciparum causes the most severe malaria and is exposed to various environmental and physiological stresses in the human host. Given that GCN5 plays a critical role in regulating stress responses in model organisms, we aimed to elucidate PfGCN5's function in stress responses in P. falciparum. The protein level of PfGCN5 was substantially induced under three stress conditions [heat shock, low glucose starvation, and dihydroartemisinin, the active metabolite of artemisinin (ART)]. With a TetR-DOZI conditional knockdown (KD) system, we successfully down-regulated PfGCN5 to ~50% and found that KD parasites became more sensitive to all three stress conditions. Transcriptomic analysis via RNA-seq identified ~1,000 up- and down-regulated genes in the wild-type (WT) and KD parasites under these stress conditions. Importantly, DHA induced transcriptional alteration of many genes involved in many aspects of stress responses, which were heavily shared among the altered genes under heat shock and low glucose conditions, including ART-resistance-related genes such as K13 and coronin. Based on the expression pattern between WT and KD parasites under three stress conditions, ~300-400 genes were identified to be involved in PfGCN5-dependent, general, and stress-condition-specific responses with high levels of overlaps among three stress conditions. Notably, using ring-stage survival assay, we found that KD or inhibition of PfGCN5 could sensitize the ART-resistant parasites to the DHA treatment. All these indicate that PfGCN5 is pivotal in regulating general and ART-resistance-related stress responses in malaria parasites, implicating PfGCN5 as a potential target for malaria intervention.

Extracellular vesicles could be a putative posttranscriptional regulatory mechanism that shapes intracellular RNA levels in Plasmodium falciparum

Plasmodium falciparum secretes extracellular vesicles (PfEVs) that contain parasite-derived RNA. However, the significance of the secreted RNA remains unexplored. Here, we compare secreted and intracellular RNA from asexual cultures of six P. falciparum lines. We find that secretion of RNA via extracellular vesicles is not only periodic throughout the asexual intraerythrocytic developmental cycle but is also highly conserved across P. falciparum isolates. We further demonstrate that the phases of RNA secreted via extracellular vesicles are discernibly shifted compared to those of the intracellular RNA within the secreting whole parasite. Finally, transcripts of genes with no known function during the asexual intraerythrocytic developmental cycle are enriched in PfEVs compared to the whole parasite. We conclude that the secretion of extracellular vesicles could be a putative posttranscriptional RNA regulation mechanism that is part of or synergise the classic RNA decay processes to maintain intracellular RNA levels in P. falciparum.

Epigenetic Regulation and Chromatin Remodeling in Malaria Parasites

Plasmodium falciparum, the human malaria parasite, infects two hosts and various cell types, inducing distinct morphological and physiological changes in the parasite in response to different environmental conditions. These variations required the parasite to adapt and develop elaborate molecular mechanisms to ensure its spread and transmission. Recent findings have significantly improved our understanding of the regulation of gene expression in P. falciparum. Here, we provide an up-to-date overview of technologies used to highlight the transcriptomic adjustments occurring in the parasite throughout its life cycle. We also emphasize the complementary and complex epigenetic mechanisms regulating gene expression in malaria parasites. This review concludes with an outlook on the chromatin architecture, the remodeling systems, and how this 3D genome organization is critical in various biological processes.

The global transcriptome of Plasmodium falciparum mid-stage gametocytes (stages II-IV) appears largely conserved and gametocyte-specific gene expression patterns vary in clinical isolates

Our overall understanding of the developmental biology of malaria parasites has been greatly enhanced by recent advances in transcriptomic analysis. However, most of these investigations rely on laboratory strains (LS) that were adapted into in vitro culture many years ago, and the transcriptomes of clinical isolates (CI) circulating in human populations have not been assessed. In this study, RNA-seq was used to compare the global transcriptome of mid-stage gametocytes derived from three short-term cultured CI, with gametocytes derived from the NF54 reference laboratory strain. The core transcriptome appeared to be consistent between CI- and LS-derived gametocyte preparations, but some important differences were also observed. A majority of gametocyte-specific genes (43/53) appear to have relatively higher expression in CI-derived gametocytes than in LS-derived gametocytes, but a K-means clustering analysis showed that genes involved in flagellum- and microtubule-based processes (movement/motility) were more abundant in both groups, albeit with some differences between them. In addition, gametocytes from one CI described as CI group II gametocytes (CI:GGII) showed gene expression variation in the form of reduced gametocyte-specific gene expression compared to the other two CI-derived gametocytes (CI gametocyte group I, CI:GGI), although the mixed developmental stages used in our study is a potential confounder, only partially mitigated by the inclusion of multiple replicates for each CI. Overall, our study suggests that there may be subtle differences in the gene expression profiles of mid-stage gametocytes from CI relative to the NF54 reference strain of Plasmodium falciparum. Thus, it is necessary to deploy gametocyte-producing clinical parasite isolates to fully understand the diversity of gene expression strategies that may occur during the sequestered development of parasite sexual stages. IMPORTANCE Maturing gametocytes of Plasmodium falciparum are known to sequester away from peripheral circulation into the bone marrow until they are mature. Blocking gametocyte sequestration can prevent malaria transmission from humans to mosquitoes, but most studies aim to understand gametocyte development utilizing long-term adapted laboratory lines instead of clinical isolates. This is a particular issue for our understanding of the sexual stages, which are known to decrease rapidly during adaptation to long-term culture, meaning that many LS are unable to produce transmissible gametocytes. Using RNA-seq, we investigated the global transcriptome of mid-stage gametocytes derived from three clinical isolates and a reference strain (NF54). This identified important differences in gene expression profiles between immature gametocytes of CI and the NF54 reference strain of P. falciparum, suggesting increased investment in gametocytogenesis in clinical isolates. Our transcriptomic data highlight the use of clinical isolates in studying the morphological, cellular features and molecular biology of gametocytes.

A novel Modulator of Ring Stage Translation (MRST) gene alters artemisinin sensitivity in Plasmodium falciparum

The implementation of artemisinin (ART) combination therapies (ACTs) has greatly decreased deaths caused by Plasmodium falciparum malaria, but increasing ACT resistance in Southeast Asia and Africa could reverse this progress. Parasite population genetic studies have identified numerous genes, single-nucleotide polymorphisms (SNPs), and transcriptional signatures associated with altered artemisinin activity with SNPs in the Kelch13 (K13) gene being the most well-characterized artemisinin resistance marker. However, there is an increasing evidence that resistance to artemisinin in P. falciparum is not related only to K13 SNPs, prompting the need to characterize other novel genes that can alter ART responses in P. falciparum. In our previous analyses of P. falciparum piggyBac mutants, several genes of unknown function exhibited increased sensitivity to artemisinin that was similar to a mutant of K13. Further analysis of these genes and their gene co-expression networks indicated that the ART sensitivity cluster was functionally linked to DNA replication and repair, stress responses, and maintenance of homeostatic nuclear activity. In this study, we have characterized PF3D7_1136600, another member of the ART sensitivity cluster. Previously annotated as a conserved Plasmodium gene of unknown function, we now provide putative annotation of this gene as a Modulator of Ring Stage Translation (MRST). Our findings reveal that the mutagenesis of MRST affects gene expression of multiple translation-associated pathways during the early ring stage of asexual development via putative ribosome assembly and maturation activity, suggesting an essential role of MRST in protein biosynthesis and another novel mechanism of altering the parasite's ART drug response.IMPORTANCEPlasmodium falciparum malaria killed more than 600,000 people in 2021, though ACTs have been critical in reducing malaria mortality as a first-line treatment for infection. However, ACT resistance in Southeast Asia and emerging resistance in Africa are detrimental to this progress. Mutations to Kelch13 (K13) have been identified to confer increased artemisinin tolerance in field isolates, however, genes other than K13 are implicated in altering how the parasite responds to artemisinin prompts additional analysis. Therefore, in this study we have characterized a P. falciparum mutant clone with altered sensitivity to artemisinin and identified a novel gene (PF3D7_1136600) that is associated with alterations to parasite translational metabolism during critical timepoints for artemisinin drug response. Many genes of the P. falciparum genome remain unannotated, posing a challenge for drug-gene characterizations in the parasite. Therefore, through this study, we have putatively annotated PF3D7_1136600 as a novel MRST gene and have identified a potential link between MRST and parasite stress response mechanisms.

Novel insights into the role of long non-coding RNA in the human malaria parasite, Plasmodium falciparum

The complex life cycle of Plasmodium falciparum requires coordinated gene expression regulation to allow host cell invasion, transmission, and immune evasion. Increasing evidence now suggests a major role for epigenetic mechanisms in gene expression in the parasite. In eukaryotes, many lncRNAs have been identified to be pivotal regulators of genome structure and gene expression. To investigate the regulatory roles of lncRNAs in P. falciparum we explore the intergenic lncRNA distribution in nuclear and cytoplasmic subcellular locations. Using nascent RNA expression profiles, we identify a total of 1768 lncRNAs, of which 718 (~41%) are novels in P. falciparum. The subcellular localization and stage-specific expression of several putative lncRNAs are validated using RNA-FISH. Additionally, the genome-wide occupancy of several candidate nuclear lncRNAs is explored using ChIRP. The results reveal that lncRNA occupancy sites are focal and sequence-specific with a particular enrichment for several parasite-specific gene families, including those involved in pathogenesis and sexual differentiation. Genomic and phenotypic analysis of one specific lncRNA demonstrate its importance in sexual differentiation and reproduction. Our findings bring a new level of insight into the role of lncRNAs in pathogenicity, gene regulation and sexual differentiation, opening new avenues for targeted therapeutic strategies against the deadly malaria parasite.

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.

The role of long noncoding RNAs in malaria parasites

The human malaria parasites, including Plasmodium falciparum, persist as a major cause of global morbidity and mortality. The recent stalling of progress toward malaria elimination substantiates a need for novel interventions. Controlled gene expression is central to the parasite's numerous life cycle transformations and adaptation. With few specific transcription factors (TFs) identified, crucial roles for chromatin states and epigenetics in parasite transcription have become evident. Although many chromatin-modifying enzymes are known, less is known about which factors mediate their impacts on transcriptional variation. Like those of higher eukaryotes, long noncoding RNAs (lncRNAs) have recently been shown to have integral roles in parasite gene regulation. This review aims to summarize recent developments and key findings on the role of lncRNAs in P. falciparum.

Roles of the RON3 C-terminal fragment in erythrocyte invasion and blood-stage parasite proliferation in Plasmodium falciparum

Plasmodium species cause malaria, and in the instance of Plasmodium falciparum is responsible for a societal burden of over 600,000 deaths annually. The symptoms and pathology of malaria are due to intraerythocytic parasites. Erythrocyte invasion is mediated by the parasite merozoite stage, and is accompanied by the formation of a parasitophorous vacuolar membrane (PVM), within which the parasite develops. The merozoite apical rhoptry organelle contains various proteins that contribute to erythrocyte attachment and invasion. RON3, a rhoptry bulb membrane protein, undergoes protein processing and is discharged into the PVM during invasion. RON3-deficient parasites fail to develop beyond the intraerythrocytic ring stage, and protein export into erythrocytes by the Plasmodium translocon of exported proteins (PTEX) apparatus is abrogated, as well as glucose uptake into parasites. It is known that truncated N- and C-terminal RON3 fragments are present in rhoptries, but it is unclear which RON3 fragments contribute to protein export by PTEX and glucose uptake through the PVM. To investigate and distinguish the roles of the RON3 C-terminal fragment at distinct developmental stages, we used a C-terminus tag for conditional and post-translational control. We demonstrated that RON3 is essential for blood-stage parasite survival, and knockdown of RON3 C-terminal fragment expression from the early schizont stage induces a defect in erythrocyte invasion and the subsequent development of ring stage parasites. Protein processing of full-length RON3 was partially inhibited in the schizont stage, and the RON3 C-terminal fragment was abolished in subsequent ring-stage parasites compared to the RON3 N-terminal fragment. Protein export and glucose uptake were abrogated specifically in the late ring stage. Plasmodial surface anion channel (PSAC) activity was partially retained, facilitating small molecule traffic across the erythrocyte membrane. The knockdown of the RON3 C-terminal fragment after erythrocyte invasion did not alter parasite growth. These data suggest that the RON3 C-terminal fragment participates in erythrocyte invasion and serves an essential role in the progression of ring-stage parasite growth by the establishment of the nutrient-permeable channel in the PVM, accompanying the transport of ring-stage parasite protein from the plasma membrane to the PVM.

A type II protein arginine methyltransferase regulates merozoite invasion in Plasmodium falciparum

Protein arginine methyltransferases (PRMTs) regulate many important cellular processes, such as transcription and RNA processing in model organisms but their functions in human malaria parasites are not elucidated. Here, we characterize PfPRMT5 in Plasmodium falciparum, which catalyzes symmetric dimethylation of histone H3 at R2 (H3R2me2s) and R8, and histone H4 at R3 in vitro. PfPRMT5 disruption results in asexual stage growth defects primarily due to lower invasion efficiency of the merozoites. Transcriptomic analysis reveals down-regulation of many transcripts related to invasion upon PfPRMT5 disruption, in agreement with H3R2me2s being an active chromatin mark. Genome-wide chromatin profiling detects extensive H3R2me2s marking of genes of different cellular processes, including invasion-related genes in wildtype parasites and PfPRMT5 disruption leads to the depletion of H3R2me2s. Interactome studies identify the association of PfPRMT5 with invasion-related transcriptional regulators such as AP2-I, BDP1, and GCN5. Furthermore, PfPRMT5 is associated with the RNA splicing machinery, and PfPRMT5 disruption caused substantial anomalies in RNA splicing events, including those for invasion-related genes. In summary, PfPRMT5 is critical for regulating parasite invasion and RNA splicing in this early-branching eukaryote.

Generation of a mutator parasite to drive resistome discovery in Plasmodium falciparum

In vitro evolution of drug resistance is a powerful approach for identifying antimalarial targets, however, key obstacles to eliciting resistance are the parasite inoculum size and mutation rate. Here we sought to increase parasite genetic diversity to potentiate resistance selections by editing catalytic residues of Plasmodium falciparum DNA polymerase δ. Mutation accumulation assays reveal a ~5-8 fold elevation in the mutation rate, with an increase of 13-28 fold in drug-pressured lines. Upon challenge with the spiroindolone PfATP4-inhibitor KAE609, high-level resistance is obtained more rapidly and at lower inocula than wild-type parasites. Selections also yield mutants with resistance to an "irresistible" compound, MMV665794 that failed to yield resistance with other strains. We validate mutations in a previously uncharacterised gene, PF3D7_1359900, which we term quinoxaline resistance protein (QRP1), as causal for resistance to MMV665794 and a panel of quinoxaline analogues. The increased genetic repertoire available to this "mutator" parasite can be leveraged to drive P. falciparum resistome discovery.

Long-Read Genome Assembly and Gene Model Annotations for the Rodent Malaria Parasite Plasmodium yoelii 17XNL

Malaria causes over 600 thousand fatalities each year, with most cases attributed to the human-infectious Plasmodium falciparum species. Many rodent-infectious Plasmodium species, like Plasmodium berghei and Plasmodium yoelii, have been used as model species that can expedite studies of this pathogen. P. yoelii is an especially good model for investigating the mosquito and liver stages of parasite development because key attributes closely resemble those of P. falciparum. Because of its importance, in 2002 the 17XNL strain of P. yoelii was the first rodent malaria parasite to be sequenced. While a breakthrough effort, the assembly consisted of >5000 contiguous sequences that adversely impacted the annotated gene models. While other rodent malaria parasite genomes have been sequenced and annotated since then, including the related P. yoelii 17X strain, the 17XNL strain has not. As a result, genomic data for 17X has become the de facto reference genome for the 17XNL strain while leaving open questions surrounding possible differences between the 17XNL and 17X genomes. In this work, we present a high-quality genome assembly for P. yoelii 17XNL using PacBio DNA sequencing. In addition, we use Nanopore and Illumina RNA sequencing of mixed blood stages to create complete gene models that include coding sequences, alternate isoforms, and UTR designations. A comparison of the 17X and this new 17XNL assembly revealed biologically meaningful differences between the strains due to the presence of coding sequence variants. Taken together, our work provides a new genomic framework for studies with this commonly used rodent malaria model species.

Cripowellins Pause Plasmodium falciparum Intraerythrocytic Development at the Ring Stage

Cripowellins from Crinum erubescens are known pesticidal and have potent antiplasmodial activity. To gain mechanistic insights to this class of natural products, studies to determine the timing of action of cripowellins within the asexual intraerythrocytic cycle of Plasmodium falciparum were performed and led to the observation that this class of natural products induced reversible cytostasis in the ring stage within the first 24 h of treatment. The transcriptional program necessary for P. falciparum to progress through the asexual intraerythrocytic life cycle is well characterized. Whole transcriptome abundance analysis showed that cripowellin B “pauses” the transcriptional program necessary to progress through the intraerythrocytic life cycle coinciding with the lack of morphological progression of drug treated parasites. In addition, cripowellin B-treated parasites re-enter transcriptional progression after treatment was removed. This study highlights the use of cripowellins as chemical probes to reveal new aspects of cell cycle progression of the asexual ring stage of P. falciparum which could be leveraged for the generation of future antimalarial therapeutics.

Regulators of male and female sexual development are critical for the transmission of a malaria parasite

Malaria transmission to mosquitoes requires a developmental switch in asexually dividing blood-stage parasites to sexual reproduction. In Plasmodium berghei, the transcription factor AP2-G is required and sufficient for this switch, but how a particular sex is determined in a haploid parasite remains unknown. Using a global screen of barcoded mutants, we here identify genes essential for the formation of either male or female sexual forms and validate their importance for transmission. High-resolution single-cell transcriptomics of ten mutant parasites portrays the developmental bifurcation and reveals a regulatory cascade of putative gene functions in the determination and subsequent differentiation of each sex. A male-determining gene with a LOTUS/OST-HTH domain as well as the protein interactors of a female-determining zinc-finger protein indicate that germ-granule-like ribonucleoprotein complexes complement transcriptional processes in the regulation of both male and female development of a malaria parasite.

Protein KIC5 is a novel regulator of artemisinin stress response in the malaria parasite Plasmodium falciparum

Artemisinin combination therapies (ACTs) have led to a significant decrease in Plasmodium falciparum malaria mortality. This progress is now threatened by emerging artemisinin resistance (ART-R) linked originally in SE Asia to polymorphisms in the Kelch propeller protein (K13) and more recently to several other seemingly unrelated genetic mutations. To better understand the parasite response to ART, we are characterizing a P. falciparum mutant with altered sensitivity to ART that was created via piggyBac transposon mutagenesis. The transposon inserted near the putative transcription start site of a gene defined as a "Plasmodium-conserved gene of unknown function," now functionally linked to K13 as the Kelch13 Interacting Candidate 5 protein (KIC5). Phenotype analysis of the KIC5 mutant during intraerythrocytic asexual development identified transcriptional changes associated with DNA stress response and altered mitochondrial metabolism, linking dysregulation of the KIC5 gene to the parasite's ability to respond to ART exposure. Through characterization of the KIC5 transcriptome, we hypothesize that this gene may be essential under ART exposure to manage gene expression of the wild-type stress response at early ring stage, thereby providing a better understanding of the parasite's processes that can alter ART sensitivity.

Long-Read Genome Assembly and Gene Model Annotations for the Rodent Malaria Parasite Plasmodium yoelii 17XNL

Malaria causes over 200 million infections and over 600 thousand fatalities each year, with most cases attributed to a human-infectious Plasmodium species, Plasmodium falciparum . Many rodent-infectious Plasmodium species, like Plasmodium berghei, Plasmodium chabaudi , and Plasmodium yoelii , have been used as genetically tractable model species that can expedite studies of this pathogen. In particular, P. yoelii is an especially good model for investigating the mosquito and liver stages of parasite development because key attributes closely resemble those of P. falciparum . Because of its importance to malaria research, in 2002 the 17XNL strain of P. yoelii was the first rodent malaria parasite to be sequenced. While sequencing and assembling this genome was a breakthrough effort, the final assembly consisted of >5000 contiguous sequences that impacted the creation of annotated gene models. While other important rodent malaria parasite genomes have been sequenced and annotated since then, including the related P. yoelii 17X strain, the 17XNL strain has not. As a result, genomic data for 17X has become the de facto reference genome for the 17XNL strain while leaving open questions surrounding possible differences between the 17XNL and 17X genomes. In this work, we present a high-quality genome assembly for P. yoelii 17XNL using HiFi PacBio long-read DNA sequencing. In addition, we use Nanopore long-read direct RNA-seq and Illumina short-read sequencing of mixed blood stages to create complete gene models that include not only coding sequences but also alternate transcript isoforms, and 5' and 3' UTR designations. A comparison of the 17X and this new 17XNL assembly revealed biologically meaningful differences between the strains due to the presence of coding sequence variants. Taken together, our work provides a new genomic and gene expression framework for studies with this commonly used rodent malaria model species.

A manually curated annotation characterises genomic features of P. falciparum lncRNAs

BACKGROUND

Important regulation occurs at the level of transcription in Plasmodium falciparum and growing evidence suggests that these apicomplexan parasites have complex regulatory networks. Recent studies implicate long noncoding RNAs (lncRNAs) as transcriptional regulators in P. falciparum. However, due to limited research and the lack of necessary experimental tools, our understanding of their role in the malaria-causing parasite remains largely unelucidated. In this work, we address one of these limitations, the lack of an updated and improved lncRNA annotation in P. falciparum.

RESULTS

We generated long-read RNA sequencing data and integrated information extracted and curated from multiple sources to manually annotate lncRNAs. We identified 1119 novel lncRNAs and validated and refined 1250 existing annotations. Utilising the collated datasets, we generated evidence-based ranking scores for each annotation and characterised the distinct genomic contexts and features of P. falciparum lncRNAs. Certain features indicated subsets with potential biological significance such as 25 lncRNAs containing multiple introns, 335 lncRNAs lacking mutations in piggyBac mutagenic studies and lncRNAs associated with specific biologic processes including two new types of lncRNAs found proximal to var genes.

CONCLUSIONS

The insights and the annotation presented in this study will serve as valuable tools for researchers seeking to understand the role of lncRNAs in parasite biology through both bioinformatics and experimental approaches.

Transcriptomic complexity of the human malaria parasite Plasmodium falciparum revealed by long-read sequencing

The Plasmodium falciparum human malaria parasite genome is incompletely annotated and does not accurately represent the transcriptomic diversity of this species. To address this need, we performed long-read transcriptomic sequencing. 5' capped mRNA was enriched from samples of total and nuclear-fractionated RNA from intra-erythrocytic stages and converted to cDNA library. The cDNA libraries were sequenced on PacBio and Nanopore long-read platforms. 12,495 novel isoforms were annotated from the data. Alternative 5' and 3' ends represent the majority of isoform events among the novel isoforms, with retained introns being the next most common event. The majority of alternative 5' ends correspond to genomic regions with features similar to those of the reference transcript 5' ends. However, a minority of alternative 5' ends showed markedly different features, including locations within protein-coding regions. Alternative 3' ends showed similar features to the reference transcript 3' ends, notably adenine-rich termination signals. Distinguishing features of retained introns could not be observed, except for a tendency towards shorter length and greater GC content compared with spliced introns. Expression of antisense and retained intron isoforms was detected at different intra-erythrocytic stages, suggesting developmental regulation of these isoform events. To gain insights into the possible functions of the novel isoforms, their protein-coding potential was assessed. Variants of P. falciparum proteins and novel proteins encoded by alternative open reading frames suggest that P. falciparum has a greater proteomic repertoire than the current annotation. We provide a catalog of annotated transcripts and encoded alternative proteins to support further studies on gene and protein regulation of this pathogen.

Inhibitors of ApiAP2 protein DNA binding exhibit multistage activity against Plasmodium parasites

Plasmodium parasites are reliant on the Apicomplexan AP2 (ApiAP2) transcription factor family to regulate gene expression programs. AP2 DNA binding domains have no homologs in the human or mosquito host genomes, making them potential antimalarial drug targets. Using an in-silico screen to dock thousands of small molecules into the crystal structure of the AP2-EXP (Pf3D7_1466400) AP2 domain (PDB:3IGM), we identified putative AP2-EXP interacting compounds. Four compounds were found to block DNA binding by AP2-EXP and at least one additional ApiAP2 protein. Our top ApiAP2 competitor compound perturbs the transcriptome of P. falciparum trophozoites and results in a decrease in abundance of log₂ fold change > 2 for 50% (46/93) of AP2-EXP target genes. Additionally, two ApiAP2 competitor compounds have multi-stage anti-Plasmodium activity against blood and mosquito stage parasites. In summary, we describe a novel set of antimalarial compounds that interact with AP2 DNA binding domains. These compounds may be used for future chemical genetic interrogation of ApiAP2 proteins or serve as starting points for a new class of antimalarial therapeutics.

Plasmodium parasites are the causative agent of malaria, which resulted in over 600,000 deaths in 2021. Due to resistance arising for every antimalarial therapeutic deployed to date, new drug targets and druggable pathways must be explored. To address this concern, we used a molecular docking screen to predict competitors of DNA binding by the parasite specific family of Apicomplexan AP2 (ApiAP2) transcription factor proteins for testing in vitro and in vivo. We find that ApiAP2 competing compounds have antimalarial activity consistent with the disruption of gene regulation. This work will further our understanding of both the biological role and targetability of parasite transcriptional regulation.

Genomic insights into host and parasite interactions during intracellular infection by Toxoplasma gondii

To gain insights into the molecular interactions of an intracellular pathogen and its host cell, we studied the gene expression and chromatin states of human fibroblasts infected with the Apicomplexan parasite Toxoplasma gondii. We show a striking activation of host cell genes that regulate a number of cellular processes, some of which are protective of the host cell, others likely to be advantageous to the pathogen. The simultaneous capture of host and parasite genomic information allowed us to gain insights into the regulation of the T. gondii genome. We show how chromatin accessibility and transcriptional profiling together permit novel annotation of the parasite's genome, including more accurate mapping of known genes and the identification of new genes and cis-regulatory elements. Motif analysis reveals not only the known T. gondii AP2 transcription factor-binding site but also a previously-undiscovered candidate TATA box-containing motif at one-quarter of promoters. By inferring the transcription factor and upstream cell signaling responses involved in the host cell, we can use genomic information to gain insights into T. gondii's perturbation of host cell physiology. Our resulting model builds on previously-described human host cell signalling responses to T. gondii infection, linked to induction of specific transcription factors, some of which appear to be solely protective of the host cell, others of which appear to be co-opted by the pathogen to enhance its own survival.

Comparative single-cell transcriptional atlases of Babesia species reveal conserved and species-specific expression profiles

Babesia is a genus of apicomplexan parasites that infect red blood cells in vertebrate hosts. Pathology occurs during rapid replication cycles in the asexual blood stage of infection. Current knowledge of Babesia replication cycle progression and regulation is limited and relies mostly on comparative studies with related parasites. Due to limitations in synchronizing Babesia parasites, fine-scale time-course transcriptomic resources are not readily available. Single-cell transcriptomics provides a powerful unbiased alternative for profiling asynchronous cell populations. Here, we applied single-cell RNA sequencing to 3 Babesia species (B. divergens, B. bovis, and B. bigemina). We used analytical approaches and algorithms to map the replication cycle and construct pseudo-synchronized time-course gene expression profiles. We identify clusters of co-expressed genes showing "just-in-time" expression profiles, with gradually cascading peaks throughout asexual development. Moreover, clustering analysis of reconstructed gene curves reveals coordinated timing of peak expression in epigenetic markers and transcription factors. Using a regularized Gaussian graphical model, we reconstructed co-expression networks and identified conserved and species-specific nodes. Motif analysis of a co-expression interactome of AP2 transcription factors identified specific motifs previously reported to play a role in DNA replication in Plasmodium species. Finally, we present an interactive web application to visualize and interactively explore the datasets.

Efficient generation of mNeonGreen Plasmodium falciparum reporter lines enables quantitative fitness analysis

CRISPR editing has enabled the rapid creation of fluorescent Plasmodium transgenic lines, facilitating a deeper understanding of parasite biology. The impact of genetic perturbations such as gene disruption or the introduction of drug resistance alleles on parasite fitness is typically quantified in competitive growth assays between the query line and a wild type reference. Although fluorescent reporter lines offer a facile and frequently used method to measure relative growth, this approach is limited by the strain background of the existing reporter, which may not match the growth characteristics of the query strains, particularly if these are slower-growing field isolates. Here, we demonstrate an efficient CRISPR-based approach to generate fluorescently labelled parasite lines using mNeonGreen derived from the LanYFP protein in Branchiostoma lanceolatum, which is one of the brightest monomeric green fluorescent proteins identified. Using a positive-selection approach by insertion of an in-frame blasticidin S deaminase marker, we generated a Dd2 reporter line expressing mNeonGreen under the control of the pfpare (P. falciparum Prodrug Activation and Resistance Esterase) locus. We selected the pfpare locus as an integration site because it is highly conserved across P. falciparum strains, expressed throughout the intraerythrocytic cycle, not essential, and offers the potential for negative selection to further enrich for integrants. The mNeonGreen@pare line demonstrates strong fluorescence with a negligible fitness defect. In addition, the construct developed can serve as a tool to fluorescently tag other P. falciparum strains for in vitro experimentation.

Emerging biology of noncoding RNAs in malaria parasites

In eukaryotic organisms, noncoding RNAs (ncRNAs) have been implicated as important regulators of multifaceted biological processes, including transcriptional, posttranscriptional, and epigenetic regulation of gene expression. In recent years, it is becoming clear that protozoan parasites encode diverse ncRNA transcripts; however, little is known about their cellular functions. Recent advances in high-throughput “omic” studies identified many novel long ncRNAs (lncRNAs) in apicomplexan parasites, some of which undergo splicing, polyadenylation, and encode small proteins. To date, only a few of them are characterized, leaving a big gap in our understanding regarding their origin, mode of action, and functions in parasite biology. In this review, we focus on lncRNAs of the human malaria parasite Plasmodium falciparum and highlight their cellular functions and possible mechanisms of action.

Evaluation of the Effect of Gene Duplication by Genome Editing on Drug Resistance in Plasmodium falciparum

The emergence and spread of drug-resistant Plasmodium falciparum have compromised antimalarial efficacy and threatened the global malaria elimination campaign using artemisinin combination therapies. The impacts of amino acid substitutions in antimalarial drug resistance-associated genes on drug susceptibility have been investigated; however, the effects of amplification of those genes remain unexplored due to the lack of robust genetic approaches. Here, we generated transgenic P. falciparum parasites with an additional copy of a drug resistance-associated gene using the highly efficient CRISPR/Cas9 system and investigated their drug response. Insertion of a drug resistance-associated gene expression cassette in the genome resulted in a roughly twofold increase of mRNA levels of the target gene mdr1, which encodes multidrug resistance protein 1. The gene duplication event contributed to resistance to mefloquine, lumefantrine, and dihydroartemisinin, while the duplication of a genomic region encoding plasmepsin 2 and plasmepsin 3 did not affect resistance to antimalarial drugs, including piperaquine. We further demonstrated that mdr1 mRNA expression levels are strongly associated with mefloquine resistance in several field-derived P. falciparum lines with various genetic backgrounds. This study provides compelling evidence that mdr1 could serve as a molecular marker for the surveillance of mefloquine-resistant parasites. Long DNA integration into parasite genomes using the CRISPR/Cas9 system provides a useful tool for the evaluation of the effect of copy number variation on drug response.

A Leak-Free Inducible CRISPRi/a System for Gene Functional Studies in Plasmodium falciparum

By fusing catalytically dead Cas9 (dCas9) to active domains of histone deacetylase (Sir2a) or acetyltransferase (GCN5), this CRISPR interference/activation (CRISPRi/a) system allows gene regulation at the transcriptional level without causing permanent changes in the parasite genome. However, the constitutive expression of dCas9 poses a challenge for studying essential genes, which may lead to adaptive changes in the parasite, masking the true phenotypes. Here, we developed a leak-free inducible CRISPRi/a system by integrating the DiCre/loxP regulon to allow the expression of dCas9-GCN5/-Sir2a upon transient induction with rapamycin, which allows convenient transcriptional regulation of a gene of interest by introducing a guide RNA targeting its transcription start region. Using eight genes that are either silent or expressed from low to high levels during asexual erythrocytic development, we evaluated the robustness and versatility of this system in the asexual parasites. For most genes analyzed, this inducible CRISPRi/a system led to 1.5- to 3-fold up-or downregulation of the target genes at the mRNA level. Alteration in the expression of PfK13 and PfMYST resulted in altered sensitivities to artemisinin. For autophagy-related protein 18, an essential gene related to artemisinin resistance, a >2-fold up- or downregulation was obtained by inducible CRISPRi/a, leading to growth retardation. For the master regulator of gametocytogenesis, PfAP2-G, a >10-fold increase of the PfAP2-G transcripts was obtained by CRISPRa, resulting in >4-fold higher gametocytemia in the induced parasites. Additionally, inducible CRISPRi/a could also regulate gene expression in gametocytes. This inducible epigenetic regulation system offers a fast way of studying gene functions in Plasmodium falciparum. IMPORTANCE Understanding the fundamental biology of malaria parasites through functional genetic/genomic studies is critical for identifying novel targets for antimalarial development. Conditional knockout/knockdown systems are required to study essential genes in the haploid blood stages of the parasite. In this study, we developed an inducible CRISPRi/a system via the integration of DiCre/loxP. We evaluated the robustness and versatility of this system by activating or repressing eight selected genes and achieved up- and downregulation of the targeted genes located in both the euchromatin and heterochromatin regions. This system offers the malaria research community another tool for functional genetic studies.

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.

The transcriptome from asexual to sexual in vitro development of Cystoisospora suis (Apicomplexa: Coccidia)

The apicomplexan parasite Cystoisospora suis is an enteropathogen of suckling piglets with woldwide distribution. As with all coccidian parasites, its lifecycle is characterized by asexual multiplication followed by sexual development with two morphologically distinct cell types that presumably fuse to form a zygote from which the oocyst arises. However, knowledge of the sexual development of C. suis is still limited. To complement previous in vitro studies, we analysed transcriptional profiles at three different time points of development (corresponding to asexual, immature and mature sexual stages) in vitro via RNASeq. Overall, transcription of genes encoding proteins with important roles in gametes biology, oocyst wall biosynthesis, DNA replication and axonema formation as well as proteins with important roles in merozoite biology was identified. A homologue of an oocyst wall tyrosine rich protein of Toxoplasma gondii was expressed in macrogametes and oocysts of C. suis. We evaluated inhibition of sexual development in a host-free culture for C. suis by antiserum specific to this protein to evaluate whether it could be exploited as a candidate for control strategies against C. suis. Based on these data, targets can be defined for future strategies to interrupt parasite transmission during sexual development.

Functional characterization of 5' UTR cis-acting sequence elements that modulate translational efficiency in Plasmodium falciparum and humans

Background

The eukaryotic parasite Plasmodium falciparum causes millions of malarial infections annually while drug resistance to common anti-malarials is further confounding eradication efforts. Translation is an attractive therapeutic target that will benefit from a deeper mechanistic understanding. As the rate limiting step of translation, initiation is a primary driver of translational efficiency. It is a complex process regulated by both cis and trans acting factors, providing numerous potential targets. Relative to model organisms and humans, P. falciparum mRNAs feature unusual 5′ untranslated regions suggesting cis-acting sequence complexity in this parasite may act to tune levels of protein synthesis through their effects on translational efficiency.

Methods

Here, in vitro translation is deployed to compare the role of cis-acting regulatory sequences in P. falciparum and humans. Using parasite mRNAs with high or low translational efficiency, the presence, position, and termination status of upstream “AUG”s, in addition to the base composition of the 5′ untranslated regions, were characterized.

Results

The density of upstream “AUG”s differed significantly among the most and least efficiently translated genes in P. falciparum, as did the average “GC” content of the 5′ untranslated regions. Using exemplars from highly translated and poorly translated mRNAs, multiple putative upstream elements were interrogated for impact on translational efficiency. Upstream “AUG”s were found to repress translation to varying degrees, depending on their position and context, while combinations of upstream “AUG”s had non-additive effects. The base composition of the 5′ untranslated regions also impacted translation, but to a lesser degree. Surprisingly, the effects of cis-acting sequences were remarkably conserved between P. falciparum and humans.

Conclusions

While translational regulation is inherently complex, this work contributes toward a more comprehensive understanding of parasite and human translational regulation by examining the impact of discrete cis-acting features, acting alone or in context.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-04024-2.

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

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

G-quadruplex RNA motifs influence gene expression in the malaria parasite Plasmodium falciparum

G-quadruplexes are non-helical secondary structures that can fold in vivo in both DNA and RNA. In human cells, they can influence replication, transcription and telomere maintenance in DNA, or translation, transcript processing and stability of RNA. We have previously showed that G-quadruplexes are detectable in the DNA of the malaria parasite Plasmodium falciparum, despite a very highly A/T-biased genome with unusually few guanine-rich sequences. Here, we show that RNA G-quadruplexes can also form in P. falciparum RNA, using rG4-seq for transcriptome-wide structure-specific RNA probing. Many of the motifs, detected here via the rG4seeker pipeline, have non-canonical forms and would not be predicted by standard in silico algorithms. However, in vitro biophysical assays verified formation of non-canonical motifs. The G-quadruplexes in the P. falciparum transcriptome are frequently clustered in certain genes and associated with regions encoding low-complexity peptide repeats. They are overrepresented in particular classes of genes, notably those that encode PfEMP1 virulence factors, stress response genes and DNA binding proteins. In vitro translation experiments and in vivo measures of translation efficiency showed that G-quadruplexes can influence the translation of P. falciparum mRNAs. Thus, the G-quadruplex is a novel player in post-transcriptional regulation of gene expression in this major human pathogen.

Identifying transcript 5' capped ends in Plasmodium falciparum

Background

The genome of the human malaria parasite Plasmodium falciparum is poorly annotated, in particular, the 5' capped ends of its mRNA transcripts. New approaches are needed to fully catalog P. falciparum transcripts for understanding gene function and regulation in this organism.

Methods

We developed a transcriptomic method based on next-generation sequencing of complementary DNA (cDNA) enriched for full-length fragments using eIF4E, a 5' cap-binding protein, and an unenriched control. DNA sequencing adapter was added after enrichment of full-length cDNA using two different ligation protocols. From the mapped sequence reads, enrichment scores were calculated for all transcribed nucleotides and used to calculate P-values of 5' capped nucleotide enrichment. Sensitivity and accuracy were increased by combining P-values from replicate experiments. Data were obtained for P. falciparum ring, trophozoite and schizont stages of intra-erythrocytic development.

Results

5' capped nucleotide signals were mapped to 17,961 non-overlapping P. falciparum genomic intervals. Analysis of the dominant 5' capped nucleotide in these genomic intervals revealed the presence of two groups with distinctive epigenetic features and sequence patterns. A total of 4,512 transcripts were annotated as 5' capped based on the correspondence of 5' end with 5' capped nucleotide annotated from full-length cDNA data.

Discussion

The presence of two groups of 5' capped nucleotides suggests that alternative mechanisms may exist for producing 5' capped transcript ends in P. falciparum. The 5' capped transcripts that are antisense, outside of, or partially overlapping coding regions may be important regulators of gene function in P. falciparum.

Plasmodium falciparum transcription in different clinical presentations of malaria associates with circulation time of infected erythrocytes

Following Plasmodium falciparum infection, individuals can remain asymptomatic, present with mild fever in uncomplicated malaria cases, or show one or more severe malaria symptoms. Several studies have investigated associations between parasite transcription and clinical severity, but no broad conclusions have yet been drawn. Here, we apply a series of bioinformatic approaches based on P. falciparum's tightly regulated transcriptional pattern during its ~48-hour intraerythrocytic developmental cycle (IDC) to publicly available transcriptomes of parasites obtained from malaria cases of differing clinical severity across multiple studies. Our analysis shows that within each IDC, the circulation time of infected erythrocytes without sequestering to endothelial cells decreases with increasing parasitaemia or disease severity. Accordingly, we find that the size of circulating infected erythrocytes is inversely related to parasite density and disease severity. We propose that enhanced adhesiveness of infected erythrocytes leads to a rapid increase in parasite burden, promoting higher parasitaemia and increased disease severity.

The RNA structurome in the asexual blood stages of malaria pathogen plasmodium falciparum

Plasmodium falciparum is a deadly human pathogen responsible for the devastating disease called malaria. In this study, we measured the differential accumulation of RNA secondary structures in coding and non-coding transcripts from the asexual developmental cycle in P. falciparum in human red blood cells. Our comprehensive analysis that combined high-throughput nuclease mapping of RNA structures by duplex RNA-seq, SHAPE-directed RNA structure validation, immunoaffinity purification and characterization of antisense RNAs collectively measured differentially base-paired RNA regions throughout the parasite's asexual RBC cycle. Our mapping data not only aligned to a diverse pool of RNAs with known structures but also enabled us to identify new structural RNA regions in the malaria genome. On average, approximately 71% of the genes with secondary structures are found to be protein coding mRNAs. The mapping pattern of these base-paired RNAs corresponded to all regions of mRNAs, including the 5' UTR, CDS and 3' UTR as well as the start and stop codons. Histone family genes which are known to form secondary structures in their mRNAs and transcripts from genes which are important for transcriptional and post-transcriptional control, such as the unique plant-like transcription factor family, ApiAP2, DNA-/RNA-binding protein, Alba3 and proteins important for RBC invasion and malaria cytoadherence also showed strong accumulation of duplex RNA reads in various asexual stages in P. falciparum. Intriguingly, our study determined stage-specific, dynamic relationships between mRNA structural contents and translation efficiency in P. falciparum asexual blood stages, suggesting an essential role of RNA structural changes in malaria gene expression programs.

Malaria in the 'Omics Era'

Genomics has revolutionised the study of the biology of parasitic diseases. The first Eukaryotic parasite to have its genome sequenced was the malaria parasite Plasmodium falciparum. Since then, Plasmodium genomics has continued to lead the way in the study of the genome biology of parasites, both in breadth—the number of Plasmodium species’ genomes sequenced—and in depth—massive-scale genome re-sequencing of several key species. Here, we review some of the insights into the biology, evolution and population genetics of Plasmodium gained from genome sequencing, and look at potential new avenues in the future genome-scale study of its biology.

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

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

Plasmodium Reproduction, Cell Size, and Transcription: How to Cope With Increasing DNA Content?

Plasmodium, the unicellular parasite that causes malaria, evolved a highly unusual mode of reproduction. During its complex life cycle, invasive or transmissive stages alternate with proliferating stages, where a single parasite can produce tens of thousands of progeny. In the clinically relevant blood stage of infection, the parasite replicates its genome up to thirty times and forms a multinucleated cell before daughter cells are assembled. Thus, within a single cell cycle, Plasmodium develops from a haploid to a polypoid cell, harboring multiple copies of its genome. Polyploidy creates several biological challenges, such as imbalances in genome output, and cells can respond to this by changing their size and/or alter the production of RNA species and protein to achieve expression homeostasis. However, the effects and possible adaptations of Plasmodium to the massively increasing DNA content are unknown. Here, we revisit and embed current Plasmodium literature in the context of polyploidy and propose potential mechanisms of the parasite to cope with the increasing gene dosage.

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.

Full-Length Transcriptome Analysis of Plasmodium falciparum by Single-Molecule Long-Read Sequencing

Malaria, an infectious disease caused by Plasmodium parasites, still accounts for amounts of deaths annually in last decades. Despite the significance of Plasmodium falciparum as a model organism of malaria parasites, our understanding of gene expression of this parasite remains largely elusive since lots of progress on its genome and transcriptome are based on assembly with short sequencing reads. Herein, we report the new version of transcriptome dataset containing all full-length transcripts over the whole asexual blood stages by adopting a full-length sequencing approach with optimized experimental conditions of cDNA library preparation. We have identified a total of 393 alternative splicing (AS) events, 3,623 long non-coding RNAs (lncRNAs), 1,555 alternative polyadenylation (APA) events, 57 transcription factors (TF), 1,721 fusion transcripts in P. falciparum. Furthermore, the shotgun proteome was performed to validate the full-length transcriptome of P. falciparum. More importantly, integration of full-length transcriptomic and proteomic data identified 160 novel small proteins in lncRNA regions. Collectively, this full-length transcriptome dataset with high quality and accuracy and the shotgun proteome analyses shed light on the complex gene expression in malaria parasites and provide a valuable resource for related functional and mechanistic researches on P. falciparum genes.

Single-cell RNA sequencing reveals developmental heterogeneity among Plasmodium berghei sporozoites

In the malaria-causing parasite's life cycle, Plasmodium sporozoites must travel from the midgut of a mosquito to the salivary glands before they can infect a mammalian host. However, only a fraction of sporozoites complete the journey. Since salivary gland invasion is required for transmission of sporozoites, insights at the molecular level can contribute to strategies for malaria prevention. Recent advances in single-cell RNA sequencing provide an opportunity to assess sporozoite heterogeneity at a resolution unattainable by bulk RNA sequencing methods. In this study, we use a droplet-based single-cell RNA sequencing workflow to analyze the transcriptomes of over 8000 Plasmodium berghei sporozoites derived from the midguts and salivary glands of Anopheles stephensi mosquitoes. The detection of known marker genes confirms the successful capture and sequencing of samples composed of a mixed population of sporozoites. Using data integration, clustering, and trajectory analyses, we reveal differences in gene expression profiles of individual sporozoites, and identify both annotated and unannotated markers associated with sporozoite development. Our work highlights the utility of a high-throughput workflow for the transcriptomic profiling of Plasmodium sporozoites, and provides new insights into gene usage during the parasite's development in the mosquito.

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,