Nascent RNA sequencing reveals mechanisms of gene regulation in the human malaria parasite Plasmodium falciparum

Charting new territory: The Plasmodium falciparum tRNA modification landscape

Ribonucleoside modifications comprising the epitranscriptome are present in all organisms and all forms of RNA, including mRNA, rRNA and tRNA, the three major RNA components of the translational machinery. Of these, tRNA is the most heavily modified and the tRNA epitranscriptome has the greatest diversity of modifications. In addition to their roles in tRNA biogenesis, quality control, structure, cleavage, and codon recognition, tRNA modifications have been shown to regulate gene expression post-transcriptionally in prokaryotes and eukaryotes, including humans. However, studies investigating the impact of tRNA modifications on gene expression in the malaria parasite Plasmodium falciparum are currently scarce. Current evidence shows that the parasite has a limited capacity for transcriptional control, which points to a heavier reliance on strategies for posttranscriptional regulation such as tRNA epitranscriptome reprogramming. This review addresses the known functions of tRNA modifications in the biology of P. falciparum while highlighting the potential therapeutic opportunities and the value of using P. falciparum as a model organism for addressing several open questions related to the tRNA epitranscriptome.

The DEAD-box RNA helicase PfDOZI imposes opposing actions on RNA metabolism in Plasmodium falciparum

In malaria parasites, the regulation of mRNA translation, storage and degradation during development and life-stage transitions remains largely unknown. Here, we functionally characterized the DEAD-box RNA helicase PfDOZI in P. falciparum. Disruption of pfdozi enhanced asexual proliferation but reduced sexual commitment and impaired gametocyte development. By quantitative transcriptomics, we show that PfDOZI is involved in the regulation of invasion-related genes and sexual stage-specific genes during different developmental stages. PfDOZI predominantly participates in processing body-like mRNPs in schizonts but germ cell granule-like mRNPs in gametocytes to impose opposing actions of degradation and protection on different mRNA targets. We further show the formation of stress granule-like mRNPs during nutritional deprivation, highlighting an essential role of PfDOZI-associated mRNPs in stress response. We demonstrate that PfDOZI participates in distinct mRNPs to maintain mRNA homeostasis in response to life-stage transition and environmental changes by differentially executing post-transcriptional regulation on the target mRNAs.

Differential Gene Expression of Malaria Parasite in Response to Red Blood Cell-Specific Glycolytic Intermediate 2,3-Diphosphoglycerate (2,3-DPG)

Innovative strategies to control malaria are urgently needed. Exploring the interplay between Plasmodium sp. parasites and host red blood cells (RBCs) offers opportunities for novel antimalarial interventions. Pyruvate kinase deficiency (PKD), characterized by heightened 2,3-diphosphoglycerate (2,3-DPG) concentration, has been associated with protection against malaria. Elevated levels of 2,3-DPG, a specific mammalian metabolite, may hinder glycolysis, prompting us to hypothesize its potential contribution to PKD-mediated protection. We investigated the impact of the extracellular supplementation of 2,3-DPG on the Plasmodium falciparum intraerythrocytic developmental cycle in vitro. The results showed an inhibition of parasite growth, resulting from significantly fewer progeny from 2,3-DPG-treated parasites. We analyzed differential gene expression and the transcriptomic profile of P. falciparum trophozoites, from in vitro cultures subjected or not subjected to the action of 2,3-DPG, using Nanopore Sequencing Technology. The presence of 2,3-DPG in the culture medium was associated with the significant differential expression of 71 genes, mostly associated with the GO terms nucleic acid binding, transcription or monoatomic anion channel. Further, several genes related to cell cycle control were downregulated in treated parasites. These findings suggest that the presence of this RBC-specific glycolytic metabolite impacts the expression of genes transcribed during the parasite trophozoite stage and the number of merozoites released from individual schizonts, which supports the potential role of 2,3-DPG in the mechanism of protection against malaria by PKD.

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.

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.

Plasmodium falciparum Development from Gametocyte to Oocyst: Insight from Functional Studies

Malaria elimination may never succeed without the implementation of transmission-blocking strategies. The transmission of Plasmodium spp. parasites from the human host to the mosquito vector depends on circulating gametocytes in the peripheral blood of the vertebrate host. Once ingested by the mosquito during blood meals, these sexual forms undergo a series of radical morphological and metabolic changes to survive and progress from the gut to the salivary glands, where they will be waiting to be injected into the vertebrate host. The design of effective transmission-blocking strategies requires a thorough understanding of all the mechanisms that drive the development of gametocytes, gametes, sexual reproduction, and subsequent differentiation within the mosquito. The drastic changes in Plasmodium falciparum shape and function throughout its life cycle rely on the tight regulation of stage-specific gene expression. This review outlines the mechanisms involved in Plasmodium falciparum sexual stage development in both the human and mosquito vector, and zygote to oocyst differentiation. Functional studies unravel mechanisms employed by P. falciparum to orchestrate the expression of stage-specific functional products required to succeed in its complex life cycle, thus providing us with potential targets for developing new therapeutics. These mechanisms are based on studies conducted with various Plasmodium species, including predominantly P. falciparum and the rodent malaria parasites P. berghei. However, the great potential of epigenetics, genomics, transcriptomics, proteomics, and functional genetic studies to improve the understanding of malaria as a disease remains partly untapped because of limitations in studies using human malaria parasites and field isolates.

Histone modification analysis reveals common regulators of gene expression in liver and blood stage merozoites of Plasmodium parasites

Gene expression in malaria parasites is subject to various layers of regulation, including histone post-translational modifications (PTMs). Gene regulatory mechanisms have been extensively studied during the main developmental stages of Plasmodium parasites inside erythrocytes, from the ring stage following invasion to the schizont stage leading up to egress. However, gene regulation in merozoites that mediate the transition from one host cell to the next is an understudied area of parasite biology. Here, we sought to characterize gene expression and the corresponding histone PTM landscape during this stage of the parasite lifecycle through RNA-seq and ChIP-seq on P. falciparum blood stage schizonts, merozoites, and rings, as well as P. berghei liver stage merozoites. In both hepatic and erythrocytic merozoites, we identified a subset of genes with a unique histone PTM profile characterized by a region of H3K4me3 depletion in their promoter. These genes were upregulated in hepatic and erythrocytic merozoites and rings, had roles in protein export, translation, and host cell remodeling, and shared a DNA motif. These results indicate that similar regulatory mechanisms may underlie merozoite formation in the liver and blood stages. We also observed that H3K4me2 was deposited in gene bodies of gene families encoding variant surface antigens in erythrocytic merozoites, which may facilitate switching of gene expression between different members of these families. Finally, H3K18me and H2K27me were uncoupled from gene expression and were enriched around the centromeres in erythrocytic schizonts and merozoites, suggesting potential roles in the maintenance of chromosomal organization during schizogony. Together, our results demonstrate that extensive changes in gene expression and histone landscape occur during the schizont-to-ring transition to facilitate productive erythrocyte infection. The dynamic remodeling of the transcriptional program in hepatic and erythrocytic merozoites makes this stage attractive as a target for novel anti-malarial drugs that may have activity against both the liver and blood stages.

Plasmodium falciparum gametocytes display global chromatin remodelling during sexual differentiation

BACKGROUND

The protozoan malaria parasite Plasmodium falciparum has a complex life cycle during which it needs to differentiate into multiple morphologically distinct life forms. A key process for transmission of the disease is the development of male and female gametocytes in the human blood, yet the mechanisms determining sexual dimorphism in these haploid, genetically identical sexual precursor cells remain largely unknown. To understand the epigenetic program underlying the differentiation of male and female gametocytes, we separated the two sexual forms by flow cytometry and performed RNAseq as well as comprehensive ChIPseq profiling of several histone variants and modifications.

RESULTS

We show that in female gametocytes the chromatin landscape is globally remodelled with respect to genome-wide patterns and combinatorial usage of histone variants and histone modifications. We identified sex specific differences in heterochromatin distribution, implicating exported proteins and ncRNAs in sex determination. Specifically in female gametocytes, the histone variants H2A.Z/H2B.Z were highly enriched in H3K9me3-associated heterochromatin. H3K27ac occupancy correlated with stage-specific gene expression, but in contrast to asexual parasites this was unlinked to H3K4me3 co-occupancy at promoters in female gametocytes.

CONCLUSIONS

Collectively, we defined novel combinatorial chromatin states differentially organising the genome in gametocytes and asexual parasites and unravelled fundamental, sex-specific differences in the epigenetic code. Our chromatin maps represent an important resource for future understanding of the mechanisms driving sexual differentiation in P. falciparum.

PfSRPK1 Regulates Asexual Blood Stage Schizogony and Is Essential for Male Gamete Formation

Serine/arginine-rich protein kinases (SRPKs) are cell cycle-regulated serine/threonine protein kinases and are important regulators of splicing factors. In this study, we functionally characterize SRPK1 of the human malaria parasite Plasmodium falciparum. P. falciparum SRPK1 (PfSRPK1) was expressed in asexual blood-stage and sexual-stage gametocytes. Pfsrpk1^- parasites formed asexual schizonts that generated far fewer merozoites than wild-type parasites, causing reduced replication rates. Pfsrpk1^- parasites also showed a severe defect in the differentiation of male gametes, causing a complete block in parasite transmission to mosquitoes. RNA sequencing (RNA-seq) analysis of wild-type PfNF54 and Pfsrpk1^- stage V gametocytes suggested a role for PfSRPK1 in regulating transcript splicing and transcript abundance of genes coding for (i) microtubule/cilium morphogenesis-related proteins, (ii) proteins involved in cyclic nucleotide metabolic processes, (iii) proteins involved in signaling such as PfMAP2, (iv) lipid metabolism enzymes, (v) proteins of osmophilic bodies, and (vi) crystalloid components. Our study reveals an essential role for PfSRPK1 in parasite cell morphogenesis and suggests this kinase as a target to prevent malaria transmission from humans to mosquitoes. IMPORTANCE Plasmodium sexual stages represent a critical bottleneck in the parasite life cycle. Gametocytes taken up in an infectious blood meal by female anopheline mosquito get activated to form gametes and fuse to form short-lived zygotes, which transform into ookinetes to infect mosquitoes. In the present study, we demonstrate that PfSRPK1 is important for merozoite formation and critical for male gametogenesis and is involved in transcript homeostasis for numerous parasite genes. Targeting PfSRPK1 and its downstream pathways may reduce parasite replication and help achieve effective malaria transmission-blocking strategies.

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 log2 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.

Epigenetic and Epitranscriptomic Gene Regulation in Plasmodium falciparum and How We Can Use It against Malaria

Malaria, caused by Plasmodium parasites, is still one of the biggest global health challenges. P. falciparum is the deadliest species to humans. In this review, we discuss how this parasite develops and adapts to the complex and heterogenous environments of its two hosts thanks to varied chromatin-associated and epigenetic mechanisms. First, one small family of transcription factors, the ApiAP2 proteins, functions as master regulators of spatio-temporal patterns of gene expression through the parasite life cycle. In addition, chromatin plasticity determines variable parasite cell phenotypes that link to parasite growth, virulence and transmission, enabling parasite adaptation within host conditions. In recent years, epitranscriptomics is emerging as a new regulatory layer of gene expression. We present evidence of the variety of tRNA and mRNA modifications that are being characterized in Plasmodium spp., and the dynamic changes in their abundance during parasite development and cell fate. We end up outlining that new biological systems, like the mosquito model, to decipher the unknowns about epigenetic mechanisms in vivo; and novel methodologies, to study the function of RNA modifications; are needed to discover the Achilles heel of the parasite. With this new knowledge, future strategies manipulating the epigenetics and epitranscriptomic machinery of the parasite have the potential of providing new weapons against malaria.

RNAseq of Infected Erythrocyte Surface Antigen-Encoding Genes

Massive parallel sequencing technology has greatly increased the breadth and depth of transcriptomic data that can be captured from P. falciparum samples. This has revolutionized in vitro studies but uptake has been slower in the analysis of clinical samples. The principal barriers are the removal of contaminating white blood cells in a malaria endemic setting and preservation of the RNA. We provide here detailed methods for the collection of purified infected erythrocytes and the preservation and extraction of RNA. We also provide methods for assessing and addressing contaminating RNA from erythroid cells, and a protocol for RNAseq library preparation optimized to maximize yield from low amounts of parasite mRNA. Finally, we provide some examples of RNAseq library characteristics that may fail quality control for other species but are in fact satisfactory for P. falciparum RNAseq.

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.

CRISPR/Cas9-engineered inducible gametocyte producer lines as a valuable tool for Plasmodium falciparum malaria transmission research

The malaria parasite Plasmodium falciparum replicates inside erythrocytes in the blood of infected humans. During each replication cycle, a small proportion of parasites commits to sexual development and differentiates into gametocytes, which are essential for parasite transmission via the mosquito vector. Detailed molecular investigation of gametocyte biology and transmission has been hampered by difficulties in generating large numbers of these highly specialised cells. Here, we engineer P. falciparum NF54 inducible gametocyte producer (iGP) lines for the routine mass production of synchronous gametocytes via conditional overexpression of the sexual commitment factor GDV1. NF54/iGP lines consistently achieve sexual commitment rates of 75% and produce viable gametocytes that are transmissible by mosquitoes. We also demonstrate that further genetic engineering of NF54/iGP parasites is a valuable tool for the targeted exploration of gametocyte biology. In summary, we believe the iGP approach developed here will greatly expedite basic and applied malaria transmission stage research.

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.

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.

Identification of long regulatory elements in the genome of Plasmodium falciparum and other eukaryotes

Long regulatory elements (LREs), such as CpG islands, polydA:dT tracts or AU-rich elements, are thought to play key roles in gene regulation but, as opposed to conventional binding sites of transcription factors, few methods have been proposed to formally and automatically characterize them. We present here a computational approach named DExTER (Domain Exploration To Explain gene Regulation) dedicated to the identification of candidate LREs (cLREs) and apply it to the analysis of the genomes of P. falciparum and other eukaryotes. Our analyses show that all tested genomes contain several cLREs that are somewhat conserved along evolution, and that gene expression can be predicted with surprising accuracy on the basis of these long regions only. Regulation by cLREs exhibits very different behaviours depending on species and conditions. In P. falciparum and other Apicomplexan organisms as well as in Dictyostelium discoideum, the process appears highly dynamic, with different cLREs involved at different phases of the life cycle. For multicellular organisms, the same cLREs are involved in all tissues, but a dynamic behavior is observed along embryonic development stages. In P. falciparum, whose genome is known to be strongly depleted of transcription factors, cLREs are predictive of expression with an accuracy above 70%, and our analyses show that they are associated with both transcriptional and post-transcriptional regulation signals. Moreover, we assessed the biological relevance of one LRE discovered by DExTER in P. falciparum using an in vivo reporter assay. The source code (python) of DExTER is available at https://gite.lirmm.fr/menichelli/DExTER.

From Genes to Transcripts, a Tightly Regulated Journey in Plasmodium

Over the past decade, we have witnessed significant progresses in understanding gene regulation in Apicomplexa including the human malaria parasite, Plasmodium falciparum. This parasite possesses the ability to convert in multiple stages in various hosts, cell types, and environments. Recent findings indicate that P. falciparum is talented at using efficient and complementary molecular mechanisms to ensure a tight control of gene expression at each stage of its life cycle. Here, we review the current understanding on the contribution of the epigenome, atypical transcription factors, and chromatin organization to regulate stage conversion in P. falciparum. The adjustment of these regulatory mechanisms occurring during the progression of the life cycle will be extensively discussed.

Differential expression of predisposing HLA-DQ2.5 alleles in DR5/DR7 celiac disease patients affects the pathological immune response to gluten

The DR5-DQ7/DR7-DQ2 genotype is very frequent among patients affected by celiac disease (CD), in Europe. This genotype, associated to high risk of CD, carries the HLA-DQA1*05 and HLA-DQB1*02 predisposing alleles, in trans configuration. The alleles encode the DQ2.5 heterodimer responsible of gluten peptide presentation on the surface of antigen-presenting cells (APCs), and consequent pathogenic CD4⁺ T cell activation. We demonstrated that DR5/DR7 APCs induce an anti-gluten CD4⁺ T cell response, of comparable intensity to that observed with APCs carrying DR1/DR3 genotype, which risk alleles are in cis configuration. In addition, we showed that DR5/DR7 APCs from celiac patients stimulated an effector CD4⁺ T cell response higher with respect to that induced by DR5/DR7 APCs from healthy subjects. To explain these findings, we assessed the DQ2.5 RNA and protein quantity. We showed that the expression of DQA1*05 and DQB1*02 risk alleles is much higher than the expression of non-CD-associated alleles, in agreement with the previous results obtained with DR1/DR3 genotype. The differential expression of transcripts influences the quantity of DQα1*05 and DQβ1*02 chains and, as consequence, the cell surface density of DQ2.5 heterodimers. Moreover, both RNA and proteins, are more abundant in APCs from celiac patients than controls. Finally, to unravel the mechanism regulating the expression of predisposing DQA1*05 and DQB1*02 alleles, we quantified the new synthetized RNA and found that the differential expression is explained by their transcription rate. Our results confirmed that the strength of antigen-specific CD4⁺ T cell response is mainly determined by the amount of gluten in the diet and provided a new possible approach for a personalized diagnosis and for risk stratification.

Dynamic Chromatin Structure and Epigenetics Control the Fate of Malaria Parasites

Multiple hosts and various life cycle stages prompt the human malaria parasite, Plasmodium falciparum, to acquire sophisticated molecular mechanisms to ensure its survival, spread, and transmission to its next host. To face these environmental challenges, increasing evidence suggests that the parasite has developed complex and complementary layers of regulatory mechanisms controlling gene expression. Here, we discuss the recent developments in the discovery of molecular components that contribute to cell replication and differentiation and highlight the major contributions of epigenetics, transcription factors, and nuclear architecture in controlling gene regulation and life cycle progression in Plasmodium spp.

Rrp6 Regulates Heterochromatic Gene Silencing via ncRNA RUF6 Decay in Malaria Parasites

Malaria remains a major public health and economic burden. The heterochromatin environment controls the silencing of genes associated with the fate of malaria parasites. Previous studies have demonstrated that a group of GC-rich ncRNAs (RUF6) is associated with the mutually exclusive expression of var genes, but the underlying mechanisms remain elusive. Here, through a series of genetic manipulation and genome-wide multiomics analysis, we have identified the plasmodial orthologue of RNA exosome-associated Rrp6 as an upstream regulator of RUF6 expression and revealed that the dysregulation of RUF6 upon Rrp6 knockdown triggered local chromatin alteration, thereby activating most heterochromatic genes via direct interaction of RUF6 and distal gene loci. This finding not only uncovered the in-depth mechanism of RUF6-mediated regulation of heterochromatic genes but also identified Rrp6 as a novel regulator of gene expression in human malaria parasites, which provides a new target for developing intervention strategies against malaria.

The heterochromatin environment plays a central role in silencing genes associated with the malaria parasite’s development, survival in the host, and transmission to the mosquito vector. However, the underlying mechanism regulating the dynamic chromatin structure is not understood yet. Here, we have uncovered that Plasmodium falciparum Rrp6, an orthologue of eukaryotic RNA exosome-associated RNase, controls the silencing of heterochromatic genes. PfRrp6 knockdown disrupted the singular expression of the GC-rich ncRNA RUF6 family, a known critical regulator of virulence gene expression, through the stabilization of the nascent transcripts. Mechanistic investigation showed that the accumulation of the multiple RUF6 ncRNAs triggered local chromatin remodeling in situ, which activated their adjacent var genes. Strikingly, chromatin isolation by RNA purification analysis (ChIRP-seq) revealed that a remarkable RUF6 ncRNA had interacted with distal heterochromatin regions directly and stimulated a global derepression effect on heterochromatic genes, including all variant gene families and the sexual commitment-associated regulator ap2-g gene. Collectively, Rrp6 appears to conduct the epigenetic surveillance of heterochromatic gene expression through controlling RUF6 levels, thereby securing antigenic variation and sexual commitment of malaria parasites during the infection of the host.

Plasmodium falciparum translational machinery condones polyadenosine repeats

Plasmodium falciparum is a causative agent of human malaria. Sixty percent of mRNAs from its extremely AT-rich (81%) genome harbor long polyadenosine (polyA) runs within their ORFs, distinguishing the parasite from its hosts and other sequenced organisms. Recent studies indicate polyA runs cause ribosome stalling and frameshifting, triggering mRNA surveillance pathways and attenuating protein synthesis. Here, we show that P. falciparum is an exception to this rule. We demonstrate that both endogenous genes and reporter sequences containing long polyA runs are efficiently and accurately translated in P. falciparum cells. We show that polyA runs do not elicit any response from No Go Decay (NGD) or result in the production of frameshifted proteins. This is in stark contrast to what we observe in human cells or T. thermophila, an organism with similar AT-content. Finally, using stalling reporters we show that Plasmodium cells evolved not to have a fully functional NGD pathway.

The second life of Plasmodium in the mosquito host: gene regulation on the move

Malaria parasites face dynamically changing environments and strong selective constraints within human and mosquito hosts. To survive such hostile and shifting conditions, Plasmodium switches transcriptional programs during development and has evolved mechanisms to adjust its phenotype through heterogeneous patterns of gene expression. In vitro studies on culture-adapted isolates have served to set the link between chromatin structure and functional gene expression. Yet, experimental evidence is limited to certain stages of the parasite in the vertebrate, i.e. blood, while the precise mechanisms underlying the dynamic regulatory landscapes during development and in the adaptation to within-host conditions remain poorly understood. In this review, we discuss available data on transcriptional and epigenetic regulation in Plasmodium mosquito stages in the context of sporogonic development and phenotypic variation, including both bet-hedging and environmentally triggered direct transcriptional responses. With this, we advocate the mosquito offers an in vivo biological model to investigate the regulatory networks, transcription factors and chromatin-modifying enzymes and their modes of interaction with regulatory sequences, which might be responsible for the plasticity of the Plasmodium genome that dictates stage- and cell type-specific blueprints of gene expression.

Functional genomics of simian malaria parasites and host-parasite interactions

Two simian malaria parasite species, Plasmodium knowlesi and Plasmodium cynomolgi, cause zoonotic infections in Southeast Asia, and they have therefore gained recognition among scientists and public health officials. Notwithstanding, these species and others including Plasmodium coatneyi have served for decades as sources of knowledge on the biology, genetics and evolution of Plasmodium, and the diverse ramifications and outcomes of malaria in their monkey hosts. Experimental analysis of these species can help to fill gaps in knowledge beyond what may be possible studying the human malaria parasites or rodent parasite species. The genome sequences for these simian malaria parasite species were reported during the last decade, and functional genomics research has since been pursued. Here research on the functional genomics analysis involving these species is summarized and their importance is stressed, particularly for understanding host–parasite interactions, and potentially testing novel interventions. Importantly, while Plasmodium falciparum and Plasmodium vivax can be studied in small New World monkeys, the simian malaria parasites can be studied more effectively in the larger Old World monkey macaque hosts, which are more closely related to humans. In addition to ex vivo analyses, experimental scenarios can include passage through Anopheline mosquito hosts and longitudinal infections in monkeys to study acute and chronic infections, as well as relapses, all in the context of the in vivo host environment. Such experiments provide opportunities for understanding functional genomic elements that govern host–parasite interactions, immunity and pathogenesis in-depth, addressing hypotheses not possible from in vitro cultures or cross-sectional clinical studies with humans.

Screening the Medicines for Malaria Venture Pathogen Box for invasion and egress inhibitors of the blood stage of Plasmodium falciparum reveals several inhibitory compounds

With emerging resistance to frontline treatments, it is vital that new drugs are identified to target Plasmodium falciparum. One of the most critical processes during parasites asexual lifecycle is the invasion and subsequent egress of red blood cells (RBCs). Many unique parasite ligands, receptors and enzymes are employed during egress and invasion that are essential for parasite proliferation and survival, therefore making these processes druggable targets. To identify potential inhibitors of egress and invasion, we screened the Medicines for Malaria Venture Pathogen Box, a 400 compound library against neglected tropical diseases, including 125 with antimalarial activity. For this screen, we utilised transgenic parasites expressing a bioluminescent reporter, nanoluciferase (Nluc), to measure inhibition of parasite egress and invasion in the presence of the Pathogen Box compounds. At a concentration of 2 µM, we found 15 compounds that inhibited parasite egress by >40% and 24 invasion-specific compounds that inhibited invasion by >90%. We further characterised 11 of these inhibitors through cell-based assays and live cell microscopy, and found two compounds that inhibited merozoite maturation in schizonts, one compound that inhibited merozoite egress, one compound that directly inhibited parasite invasion and one compound that slowed down invasion and arrested ring formation. The remaining compounds were general growth inhibitors that acted during the egress and invasion phase of the cell cycle. We found the sulfonylpiperazine, MMV020291, to be the most invasion-specific inhibitor, blocking successful merozoite internalisation within human RBCs and having no substantial effect on other stages of the cell cycle. This has significant implications for the possible development of an invasion-specific inhibitor as an antimalarial in a combination based therapy, in addition to being a useful tool for studying the biology of the invading parasite.

The cryptic unstable transcripts are associated with developmentally regulated gene expression in blood-stage Plasmodium falciparum

The tight gene expression regulation controls the development and pathogenesis of human malaria parasite Plasmodium falciparum throughout the complex life cycle. Recent studies have revealed the pervasive nascent transcripts in the genome of P. falciparum, suggesting the existence of a hidden transcriptome involved in the dynamic gene expression. However, the landscape and related biological functions of nascent non-coding RNAs (ns-ncRNAs) are still poorly explored. Here we profiled the transcription dynamics of nascent RNAs by rRNA-depleted and stranded RNA sequencing over the course of 48-h intraerythrocytic developmental cycle (IDC). We identified the genome-wide sources of a total of 2252 ns-ncRNAs, mostly originating from intergenic and untranslated regions of annotated genes. By integrating the nascent RNA abundances with ATAC-seq and ChIP-seq analysis, we uncovered the euchromatic microenvironment surrounding the ns-ncRNA loci, and revealed a positive correlation between ns-ncRNAs and corresponding mRNA abundances. Finally, by gene knock-down strategy, we showed that the cooperation of RNA exosome catalytic subunit PfDis3 and PfMtr4 cofactor played a major role in ns-ncRNAs degradation. Collectively, this study contributes to understanding of the potential roles of short-lived nascent ncRNAs in regulating gene expression in malaria parasites.

Adaptation of Translational Machinery in Malaria Parasites to Accommodate Translation of Poly-Adenosine Stretches Throughout Its Life Cycle

Malaria is caused by unicellular apicomplexan parasites of the genus Plasmodium, which includes the major human parasite Plasmodium falciparum. The complex cycle of the malaria parasite in both mosquito and human hosts has been studied extensively. There is tight control of gene expression in each developmental stage, and at every level of gene synthesis: from RNA transcription, to its subsequent translation, and finally post-translational modifications of the resulting protein. Whole-genome sequencing of P. falciparum has laid the foundation for significant biological advances by revealing surprising genomic information. The P. falciparum genome is extremely AT-rich (∼80%), with a substantial portion of genes encoding intragenic polyadenosine (polyA) tracks being expressed throughout the entire parasite life cycle. In most eukaryotes, intragenic polyA runs act as negative regulators of gene expression. Recent studies have shown that translation of mRNAs containing 12 or more consecutive adenosines results in ribosomal stalling and frameshifting; activating mRNA surveillance mechanisms. In contrast, P. falciparum translational machinery can efficiently and accurately translate polyA tracks without activating mRNA surveillance pathways. This unique feature of P. falciparum raises interesting questions: (1) How is P. falciparum able to efficiently and correctly translate polyA track transcripts, and (2) What are the specifics of the translational machinery and mRNA surveillance mechanisms that separate P. falciparum from other organisms? In this review, we analyze possible evolutionary shifts in P. falciparum protein synthesis machinery that allow efficient translation of an AU rich-transcriptome. We focus on physiological and structural differences of P. falciparum stage specific ribosomes, ribosome-associated proteins, and changes in mRNA surveillance mechanisms throughout the complete parasite life cycle, with an emphasis on the mosquito and liver stages.

Hierarchical transcriptional control regulates Plasmodium falciparum sexual differentiation

BACKGROUND

Malaria pathogenesis relies on sexual gametocyte forms of the malaria parasite to be transmitted between the infected human and the mosquito host but the molecular mechanisms controlling gametocytogenesis remains poorly understood. Here we provide a high-resolution transcriptome of Plasmodium falciparum as it commits to and develops through gametocytogenesis.

RESULTS

The gametocyte-associated transcriptome is significantly different from that of the asexual parasites, with dynamic gene expression shifts characterizing early, intermediate and late-stage gametocyte development and results in differential timing for sex-specific transcripts. The transcriptional dynamics suggest strict transcriptional control during gametocytogenesis in P. falciparum, which we propose is mediated by putative regulators including epigenetic mechanisms (driving active repression of proliferation-associated processes) and a cascade-like expression of ApiAP2 transcription factors.

CONCLUSIONS

The gametocyte transcriptome serves as the blueprint for sexual differentiation and will be a rich resource for future functional studies on this critical stage of Plasmodium development, as the intraerythrocytic transcriptome has been for our understanding of the asexual cycle.

TRIBE Uncovers the Role of Dis3 in Shaping the Dynamic Transcriptome in Malaria Parasites

Identification of RNA targets of RNA-binding proteins (RBPs) is essential for complete understanding of their biological functions. However, it is still a challenge to identify the biologically relevant targets of RBPs through in vitro strategies of RIP-seq, HITS-CLIP, or GoldCLIP due to the potentially high background and complicated manipulation. In malaria parasites, RIP-seq and gene disruption are the few tools available currently for identification of RBP targets. Here, we have adopted the TRIBE (Targets of RNA binding proteins identified by editing) system to in vivo identify the RNA targets of PfDis3, a key exoribonuclease subunit of RNA exosome in Plasmodium falciparum. We generated a transgenic parasite line of PfDis3-ADARcd, which catalyzes an adenosine (A)-to-inosine (I) conversion at the potential interacting sites of PfDis3-targeting RNAs. Most of PfDis3 target genes contain one edit site. The majority of the edit sites detected by PfDis3-TRIBE locate in exons and spread across the entire coding regions. The nucleotides adjacent to the edit sites contain ∼75% of A + T. PfDis3-TRIBE target genes are biases toward higher RIP enrichment, suggesting that PfDis3-TRIBE preferentially detects stronger PfDis3 RIP targets. Collectively, PfDis3-TRIBE is a favorable tool to identify in vivo target genes of RBP with high efficiency and reproducibility. Additionally, the PfDis3-targeting genes are involved in stage-related biological processes during the blood-stage development. Thus PfDis3 appears to shape the dynamic transcriptional transcriptome of malaria parasites through post-transcriptional degradation of a variety of unwanted transcripts from both strands in the asexual blood stage.

The role of epigenetics and chromatin structure in transcriptional regulation in malaria parasites

Due to the unique selective pressures and extreme changes faced by the human malaria parasite Plasmodium falciparum throughout its life cycle, the parasite has evolved distinct features to alter its gene expression patterns. Along with classical gene regulation by transcription factors (TFs), of which only one family, the AP2 TFs, has been described in the parasite genome, a large body of evidence points toward chromatin structure and epigenetic factors mediating the changes in gene expression associated with parasite life cycle stages. These attributes may be critically important for immune evasion, host cell invasion and development of the parasite in its two hosts, the human and the Anopheles vector. Thus, the factors involved in the maintenance and regulation of chromatin and epigenetic features represent potential targets for antimalarial drugs. In this review, we discuss the mechanisms in P. falciparum that regulate chromatin structure, nucleosome landscape, the 3-dimensional structure of the genome and additional distinctive features created by parasite-specific genes and gene families. We review conserved traits of chromatin in eukaryotes in order to highlight what is unique in the parasite.

Predicting gene expression in the human malaria parasite Plasmodium falciparum using histone modification, nucleosome positioning, and 3D localization features

Empirical evidence suggests that the malaria parasite Plasmodium falciparum employs a broad range of mechanisms to regulate gene transcription throughout the organism's complex life cycle. To better understand this regulatory machinery, we assembled a rich collection of genomic and epigenomic data sets, including information about transcription factor (TF) binding motifs, patterns of covalent histone modifications, nucleosome occupancy, GC content, and global 3D genome architecture. We used these data to train machine learning models to discriminate between high-expression and low-expression genes, focusing on three distinct stages of the red blood cell phase of the Plasmodium life cycle. Our results highlight the importance of histone modifications and 3D chromatin architecture in Plasmodium transcriptional regulation and suggest that AP2 transcription factors may play a limited regulatory role, perhaps operating in conjunction with epigenetic factors.

Ivermectin Impairs the Development of Sexual and Asexual Stages of Plasmodium falciparum In Vitro

Ivermectin is the drug of choice for many parasitic infections, with more than one billion doses being distributed in onchocerciasis programs. The drug has been put into focus recently by the malaria community because of its potential to kill blood-sucking mosquitoes, thereby reducing malaria transmission. However, the activity of ivermectin against the malaria parasite itself has been only partly investigated. This study aimed to investigate the in vitro activity of ivermectin against asexual and sexual stages of Plasmodium falciparum Both asexual and late-stage gametocytes were incubated with ivermectin and control drugs in vitro The growth-inhibiting effects were assessed for asexual stages of different Plasmodium falciparum laboratory strains and culture-adapted clinical isolates using the histidine-rich protein 2 enzyme-linked immunosorbent assay technique. The effect against stage IV/V gametocytes was evaluated based on ATP quantification. Ivermectin showed activities at nanomolar concentrations against asexual stages (50% inhibitory concentration of ∼100 nM) and stage IV/V gametocytes (500 nM) of P. falciparum Stage-specific assays suggested that ivermectin arrests the parasite cycle at the trophozoite stage. Ivermectin might add a feature to its "wonder drug" properties with activity against asexual stages of the malaria parasite Plasmodium falciparum The observed activities might be difficult to reach with current regimens but will be more relevant with future high-dose regimens under investigation. Further studies should be performed to confirm these results in vitro and in vivo.

Clinical expression and antigenic profiles of a Plasmodium vivax vaccine candidate: merozoite surface protein 7 (PvMSP-7)

Background

Vivax malaria is the predominant form of malaria outside Africa, affecting about 14 million people worldwide, with about 2.5 billion people exposed. Development of a Plasmodium vivax vaccine is a priority, and merozoite surface protein 7 (MSP-7) has been proposed as a plausible candidate. The P. vivax genome contains 12 MSP-7 genes, which contribute to erythrocyte invasion during blood-stage infection. Previous analysis of MSP-7 sequence diversity suggested that not all paralogs are functionally equivalent. To explore MSP-7 functional diversity, and to identify the best vaccine candidate within the family, MSP-7 expression and antigenicity during bloodstream infections were examined directly from clinical isolates.

Methods

Merozoite surface protein 7 gene expression was profiled using RNA-seq data from blood samples isolated from ten human patients with vivax malaria. Differential expression analysis and co-expression cluster analysis were used to relate PvMSP-7 expression to genetic markers of life cycle stage. Plasma from vivax malaria patients was also assayed using a custom peptide microarray to measure antibody responses against the coding regions of 12 MSP-7 paralogs.

Results

Ten patients presented diverse transcriptional profiles that comprised four patient groups. Two MSP-7 paralogs, 7A and 7F, were expressed abundantly in all patients, while other MSP-7 genes were uniformly rare (e.g. 7J). MSP-7H and 7I were significantly more abundant in patient group 4 only, (two patients having experienced longer patency), and were co-expressed with a schizont-stage marker, while negatively associated with liver-stage and gametocyte-stage markers. Screening infections with a PvMSP-7 peptide array identified 13 linear B-cell epitopes in five MSP-7 paralogs that were recognized by plasma from all patients.

Conclusions

These results show that MSP-7 family members vary in expression profile during blood infections; MSP-7A and 7F are expressed throughout the intraerythrocytic development cycle, while expression of other paralogs is focused on the schizont. This may reflect developmental regulation, and potentially functional differentiation, within the gene family. The frequency of B-cell epitopes among paralogs also varies, with MSP-7A and 7L consistently the most immunogenic. Thus, MSP-7 paralogs cannot be assumed to have equal potential as vaccines. This analysis of clinical infections indicates that the most abundant and immunogenic paralog is MSP-7A.

Electronic supplementary material

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

Parasite-Related Genetic and Epigenetic Aspects and Host Factors Influencing Plasmodium falciparum Invasion of Erythrocytes

Malaria, a disease caused by Plasmodium parasites, is widespread throughout tropical and sub-tropical regions worldwide; it mostly affects children and pregnant woman. Eradication has stalled despite effective prevention measures and medication being available for this disease; this has mainly been due to the parasite's resistance to medical treatment and the mosquito vector's resistance to insecticides. Tackling such resistance involves using renewed approaches and techniques for accruing a deep understanding of the parasite's biology, and developing new drugs and vaccines. Studying the parasite's invasion of erythrocytes should shed light on its ability to switch between invasion phenotypes related to the expression of gene sets encoding proteins acting as ligands during target cell invasion, thereby conferring mechanisms for evading a particular host's immune response and adapting to changes in target cell surface receptors. This review considers some factors influencing the expression of such phenotypes, such as Plasmodium's genetic, transcriptional and epigenetic characteristics, and explores some host-related aspects which could affect parasite phenotypes, aiming at integrating knowledge regarding this topic and the possible relationship between the parasite's biology and host factors playing a role in erythrocyte invasion.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

Studies of Plasmodium vivax gene expression are complicated by the lack of in vitro culture system and the difficulties associated with studying clinical infections that often contain multiple clones and a mixture of parasite stages. Here, we characterize the transcriptomes of P. vivax parasites from 26 malaria patients. We show that most parasite mRNAs derive from trophozoites and that the asynchronicity of P. vivax infections is therefore unlikely to confound gene expression studies. Analyses of gametocyte genes reveal two distinct clusters of co-regulated genes, suggesting that male and female gametocytes are independently regulated. Finally, we analyze gene expression changes induced by chloroquine and show that this antimalarial drug efficiently eliminates most P. vivax parasite stages but, in contrast to P. falciparum, does not affect trophozoites.

Schizont transcriptome variation among clinical isolates and laboratory-adapted clones of the malaria parasite Plasmodium falciparum

BACKGROUND

Malaria parasites are genetically polymorphic and phenotypically plastic. In studying transcriptome variation among parasites from different infections, it is challenging to overcome potentially confounding technical and biological variation between samples. We investigate variation in the major human parasite Plasmodium falciparum, generating RNA-seq data on multiple independent replicate sample preparations of merozoite-containing intra-erythrocytic schizonts from a panel of clinical isolates and from long-term laboratory-adapted clones, with a goal of robustly identifying differentially expressed genes.

RESULTS

Analysis of biological sample replicates shows that increased numbers improve the true discovery rate of differentially expressed genes, and that six independent replicates of each parasite line allowed identification of most differences that could be detected with larger numbers. For highly expressed genes, focusing on the top quartile at schizont stages, there was more power to detect differences. Comparing cultured clinical isolates and laboratory-adapted clones, genes more highly expressed in the laboratory-adapted clones include those encoding an AP2 transcription factor (PF3D7_0420300), a ubiquitin-binding protein and two putative methyl transferases. In contrast, higher expression in clinical isolates was seen for the merozoite surface protein gene dblmsp2, proposed to be a marker of schizonts forming merozoites committed to sexual differentiation. Variable expression was extremely strongly, but not exclusively, associated with genes known to be targeted by Heterochromatin Protein 1. Clinical isolates show variable expression of several known merozoite invasion ligands, as well as other genes for which new RT-qPCR assays validate the quantitation and allow characterisation in samples with more limited material. Expression levels of these genes vary among schizont preparations of different clinical isolates in the first ex vivo cycle in patient erythrocytes, but mean levels are similar to those in continuously cultured clinical isolates.

CONCLUSIONS

Analysis of multiple biological sample replicates greatly improves identification of genes variably expressed between different cultured parasite lines. Clinical isolates recently established in culture show differences from long-term adapted clones in transcript levels of particular genes, and are suitable for analyses requiring biological replicates to understand parasite phenotypes and variable expression likely to be relevant in nature.

tRNA epitranscriptomics and biased codon are linked to proteome expression in Plasmodium falciparum

Among components of the translational machinery, ribonucleoside modifications on tRNAs are emerging as critical regulators of cell physiology and stress response. Here, we demonstrate highly coordinated behavior of the repertoire of tRNA modifications of Plasmodium falciparum throughout the intra-erythrocytic developmental cycle (IDC). We observed both a synchronized increase in 22 of 28 modifications from ring to trophozoite stage, consistent with tRNA maturation during translational up-regulation, and asynchronous changes in six modifications. Quantitative analysis of ~2,100 proteins across the IDC revealed that up- and down-regulated proteins in late but not early stages have a marked codon bias that directly correlates with parallel changes in tRNA modifications and enhanced translational efficiency. We thus propose a model in which tRNA modifications modulate the abundance of stage-specific proteins by enhancing translation efficiency of codon-biased transcripts for critical genes. These findings reveal novel epitranscriptomic and translational control mechanisms in the development and pathogenesis of Plasmodium parasites.

Fluorescence-based assay for accurate measurement of transcriptional activity in trypanosomatid parasites

Trypanosomatid parasites are causative agents of neglected human diseases. Their lineage diverged early from the common eukaryotic ancestor, and they evolved singular mechanisms of gene expression that are crucial for their survival. Studies on unusual and essential molecular pathways lead to new drug targets. In this respect, assays to analyze transcriptional activity will provide useful information to identify essential and specific factors. However, the current methods are laborious and do not provide global and accurate measures. For this purpose, a previously reported radiolabeling in vitro nascent mRNA methodology was used to establish an alternative fluorescent-based assay that is able to precisely quantify nascent mRNA using both flow cytometry and a high-content image system. The method allowed accurate and global measurements in Trypanosoma brucei, a representative species of trypanosomatid parasites. We obtained data demonstrating that approximately 70% of parasites from a population under normal growth conditions displayed mRNA transcriptional activity, whilst the treatment with α-amanitin (75 µg/ml) inhibited the polymerase II activity. The adaptation of the method also allowed the analyses of the transcriptional activity during the cell cycle. Therefore, the methodology described herein contributes to obtaining precise measurements of transcriptional rates using multiparametric analysis. This alternative method can facilitate investigations of genetic and biochemical processes in trypanosome parasites and consequently provide additional information related to new treatment or prophylaxis strategies involving these important human parasites.

Genome-wide real-time in vivo transcriptional dynamics during Plasmodium falciparum blood-stage development

Genome-wide analysis of transcription in the malaria parasite Plasmodium falciparum has revealed robust variation in steady-state mRNA abundance throughout the 48-h intraerythrocytic developmental cycle (IDC), suggesting that this process is highly dynamic and tightly regulated. Here, we utilize rapid 4-thiouracil (4-TU) incorporation via pyrimidine salvage to specifically label, capture, and quantify newly-synthesized RNA transcripts at every hour throughout the IDC. This high-resolution global analysis of the transcriptome captures the timing and rate of transcription for each newly synthesized mRNA in vivo, revealing active transcription throughout all IDC stages. Using a statistical model to predict the mRNA dynamics contributing to the total mRNA abundance at each timepoint, we find varying degrees of transcription and stabilization for each mRNA corresponding to developmental transitions. Finally, our results provide new insight into co-regulation of mRNAs throughout the IDC through regulatory DNA sequence motifs, thereby expanding our understanding of P. falciparum mRNA dynamics.

Transcriptomic analysis often doesn’t differentiate between newly synthesized and stabilized mRNAs. Using rapid 4-thiouracil incorporation, Painter et al. here define genome-wide active transcription throughout Plasmodium blood-stage developmental stages and identify associated regulatory DNA sequence motifs.

Regulation of Sexual Commitment and Gametocytogenesis in Malaria Parasites

Sexual differentiation of malaria parasites from the asexual blood stage into gametocytes is an essential part of the life cycle, as gametocytes are the form that is taken up by the mosquito host. Because of the essentiality of this process for transmission to the mosquito, gametocytogenesis is an extremely attractive target for therapeutic interventions. The subject of this review is the considerable progress that has been made in recent years in elucidating the molecular mechanisms governing this important differentiation process. In particular, a number of critical transcription factors and epigenetic regulators have emerged as crucial elements in the regulation of commitment. The identification of these factors has allowed us to understand better than ever before the events occurring prior to and during commitment to sexual development and offers potential for new therapeutic interventions.

Changes in genome organization of parasite-specific gene families during the Plasmodium transmission stages

The development of malaria parasites throughout their various life cycle stages is coordinated by changes in gene expression. We previously showed that the three-dimensional organization of the Plasmodium falciparum genome is strongly associated with gene expression during its replication cycle inside red blood cells. Here, we analyze genome organization in the P. falciparum and P. vivax transmission stages. Major changes occur in the localization and interactions of genes involved in pathogenesis and immune evasion, host cell invasion, sexual differentiation, and master regulation of gene expression. Furthermore, we observe reorganization of subtelomeric heterochromatin around genes involved in host cell remodeling. Depletion of heterochromatin protein 1 (PfHP1) resulted in loss of interactions between virulence genes, confirming that PfHP1 is essential for maintenance of the repressive center. Our results suggest that the three-dimensional genome structure of human malaria parasites is strongly connected with transcriptional activity of specific gene families throughout the life cycle.

Characterization of the Catalytic Subunits of the RNA Exosome-like Complex in Plasmodium falciparum

The eukaryotic ribonucleic acid (RNA) exosome is a versatile multiribonuclease complex that mediates the processing, surveillance, and degradation of virtually all classes of RNA in both the nucleus and cytoplasm. The complex, composed of 10 to 11 subunits, has been widely described in many organisms. Bioinformatic analyses revealed that there may be also an exosome‐like complex in Plasmodium falciparum, a parasite of great importance in public health, with eight predicted subunits having high sequence similarity to their counterparts in yeast and human. In this work, the putative RNA catalytic components, designated as PfRrp4, PfRrp41, PfDis3, and PfRrp6, were identified and systematically analyzed. Quantitative polymerase chain reaction (QPCR) analyses suggested that all of them were transcribed steadily throughout the asexual stage. The expression of these proteins was determined by Western blot, and their localization narrowed to the cytoplasm of the parasite by indirect immunofluorescence. The recombinant proteins of PfRrp41, PfDis3, and PfRrp6 exhibited catalytic activity for single‐stranded RNA (ssRNA), whereas PfRrp4 showed no processing activity of both ssRNA and dsRNA. The identification of these putative components of the RNA exosome complex opens up new perspectives for a deep understanding of RNA metabolism in the malarial parasite P. falciparum.

Chromatin Accessibility-Based Characterization of the Gene Regulatory Network Underlying Plasmodium falciparum Blood-Stage Development

Underlying the development of malaria parasites within erythrocytes and the resulting pathogenicity is a hardwired program that secures proper timing of gene transcription and production of functionally relevant proteins. How stage-specific gene expression is orchestrated in vivo remains unclear. Here, using the assay for transposase accessible chromatin sequencing (ATAC-seq), we identified ∼4,000 regulatory regions in P. falciparum intraerythrocytic stages. The vast majority of these sites are located within 2 kb upstream of transcribed genes and their chromatin accessibility pattern correlates positively with abundance of the respective mRNA transcript. Importantly, these regions are sufficient to drive stage-specific reporter gene expression and DNA motifs enriched in stage-specific sets of regulatory regions interact with members of the P. falciparum AP2 transcription factor family. Collectively, this study provides initial insights into the in vivo gene regulatory network of P. falciparum intraerythrocytic stages and should serve as a valuable resource for future studies.

Single-cell RNA-seq reveals hidden transcriptional variation in malaria parasites

Single-cell RNA-sequencing is revolutionising our understanding of seemingly homogeneous cell populations but has not yet been widely applied to single-celled organisms. Transcriptional variation in unicellular malaria parasites from the Plasmodium genus is associated with critical phenotypes including red blood cell invasion and immune evasion, yet transcriptional variation at an individual parasite level has not been examined in depth. Here, we describe the adaptation of a single-cell RNA-sequencing (scRNA-seq) protocol to deconvolute transcriptional variation for more than 500 individual parasites of both rodent and human malaria comprising asexual and sexual life-cycle stages. We uncover previously hidden discrete transcriptional signatures during the pathogenic part of the life cycle, suggesting that expression over development is not as continuous as commonly thought. In transmission stages, we find novel, sex-specific roles for differential expression of contingency gene families that are usually associated with immune evasion and pathogenesis.