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

Malaria parasites require a divergent heme oxygenase for apicoplast gene expression and biogenesis

Malaria parasites have evolved unusual metabolic adaptations that specialize them for growth within heme-rich human erythrocytes. During blood-stage infection, Plasmodium falciparum parasites internalize and digest abundant host hemoglobin within the digestive vacuole. This massive catabolic process generates copious free heme, most of which is biomineralized into inert hemozoin. Parasites also express a divergent heme oxygenase (HO)-like protein (PfHO) that lacks key active-site residues and has lost canonical HO activity. The cellular role of this unusual protein that underpins its retention by parasites has been unknown. To unravel PfHO function, we first determined a 2.8 Å-resolution X-ray structure that revealed a highly α-helical fold indicative of distant HO homology. Localization studies unveiled PfHO targeting to the apicoplast organelle, where it is imported and undergoes N-terminal processing but retains most of the electropositive transit peptide. We observed that conditional knockdown of PfHO was lethal to parasites, which died from defective apicoplast biogenesis and impaired isoprenoid-precursor synthesis. Complementation and molecular-interaction studies revealed an essential role for the electropositive N-terminus of PfHO, which selectively associates with the apicoplast genome and enzymes involved in nucleic acid metabolism and gene expression. PfHO knockdown resulted in a specific deficiency in levels of apicoplast-encoded RNA but not DNA. These studies reveal an essential function for PfHO in apicoplast maintenance and suggest that Plasmodium repurposed the conserved HO scaffold from its canonical heme-degrading function in the ancestral chloroplast to fulfill a critical adaptive role in organelle gene expression.

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.

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 P. falciparum alternative histones Pf H2A.Z and Pf H2B.Z are dynamically acetylated and antagonized by PfSir2 histone deacetylases at heterochromatin boundaries

The Plasmodium falciparum alternative histones Pf H2A.Z and Pf H2B.Z are enriched in the same nucleosomes in intergenic euchromatin but depleted from heterochromatin. They occupy most promoters but are only dynamically associated with expression at var genes. In other organisms, acetylation of H2A.Z is important for its functions in gene expression and chromatin structure. Here, we show that acetylated Pf H2A.Z and Pf H2B.Z are dynamically associated with gene expression at promoters. In addition, acetylated Pf H2A.Z and Pf H2B.Z are antagonized by the sirtuin class III histone deacetylases (HDAC) PfSir2A and B at heterochromatin boundaries and encroach upon heterochromatin in parasites lacking PfSir2A or B. However, the majority of acetylated Pf H2A.Z and Pf H2B.Z are deacetylated by class I or II HDACs. Acetylated Pf H2A.Z and Pf H2B.Z are also dynamically associated with promoter activity of both canonical upstream var gene promoters and var gene introns. These findings suggest that both acetylated Pf H2A.Z and Pf H2B.Z play critical roles in gene expression and contribute to maintenance of chromatin structure at the boundaries of subtelomeric, facultative heterochromatin, critical for the variegated expression of genes that enable rapid adaptation to altered host environments.IMPORTANCEThe malaria parasite Plasmodium falciparum relies on variant expression of members of multi-gene families as a strategy for environmental adaptation to promote parasite survival and pathogenesis. These genes are located in transcriptionally silenced DNA regions. A limited number of these genes escape gene silencing, and switching between them confers variant fitness on parasite progeny. Here, we show that PfSir2 histone deacetylases antagonize DNA-interacting acetylated alternative histones at the boundaries between active and silent DNA. This finding implicates acetylated alternative histones in the mechanism regulating P. falciparum variant gene silencing and thus malaria pathogenesis. This work also revealed that acetylation of alternative histones at promoters is dynamically associated with promoter activity across the genome, implicating acetylation of alternative histones in gene regulation genome wide. Understanding mechanisms of gene regulation in P. falciparum may aid in the development of new therapeutic strategies for malaria, which killed 619,000 people in 2021.

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.

Cohesin contributes to transcriptional repression of stage-specific genes in the human malaria parasite

The complex life cycle of the human malaria parasite, Plasmodium falciparum, is driven by specific transcriptional programs, but it is unclear how most genes are activated or silenced at specific times. There is an association between transcription and spatial organization; however, the molecular mechanisms behind genome organization are unclear. While P. falciparum lacks key genome-organizing proteins found in metazoans, it has all core components of the cohesin complex. To investigate the role of cohesin in P. falciparum, we functionally characterize the cohesin subunit Structural Maintenance of Chromosomes protein 3 (SMC3). SMC3 knockdown during early stages of the intraerythrocytic developmental cycle (IDC) upregulates a subset of genes involved in erythrocyte egress and invasion, which are normally expressed at later stages. ChIP-seq analyses reveal that during the IDC, SMC3 enrichment at the promoter regions of these genes inversely correlates with gene expression and chromatin accessibility. These data suggest that SMC3 binding contributes to the repression of specific genes until their appropriate time of expression, revealing a new mode of stage-specific gene repression 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.

Coordination of apicoplast transcription in a malaria parasite by internal and host cues

The malaria parasite Plasmodium falciparum has a nonphotosynthetic plastid called the apicoplast, which contains its own genome. Regulatory mechanisms for apicoplast gene expression remain poorly understood, despite this organelle being crucial for the parasite life cycle. Here, we identify a nuclear-encoded apicoplast RNA polymerase σ subunit (sigma factor) which, along with the α subunit, appears to mediate apicoplast transcript accumulation. This has a periodicity reminiscent of parasite circadian or developmental control. Expression of the apicoplast subunit gene, apSig, together with apicoplast transcripts, increased in the presence of the blood circadian signaling hormone melatonin. Our data suggest that the host circadian rhythm is integrated with intrinsic parasite cues to coordinate apicoplast genome transcription. This evolutionarily conserved regulatory system might be a future target for malaria treatment.

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.

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.

Epitranscriptomics in parasitic protists: Role of RNA chemical modifications in posttranscriptional gene regulation

"Epitranscriptomics" is the new RNA code that represents an ensemble of posttranscriptional RNA chemical modifications, which can precisely coordinate gene expression and biological processes. There are several RNA base modifications, such as N6-methyladenosine (m6A), 5-methylcytosine (m5C), and pseudouridine (Ψ), etc. that play pivotal roles in fine-tuning gene expression in almost all eukaryotes and emerging evidences suggest that parasitic protists are no exception. In this review, we primarily focus on m6A, which is the most abundant epitranscriptomic mark and regulates numerous cellular processes, ranging from nuclear export, mRNA splicing, polyadenylation, stability, and translation. We highlight the universal features of spatiotemporal m6A RNA modifications in eukaryotic phylogeny, their homologs, and unique processes in 3 unicellular parasites-Plasmodium sp., Toxoplasma sp., and Trypanosoma sp. and some technological advances in this rapidly developing research area that can significantly improve our understandings of gene expression regulation in parasites.

The in vivo RNA structurome of the malaria parasite Plasmodium falciparum, a protozoan with an A/U-rich transcriptome

Plasmodium falciparum, a protozoan parasite and causative agent of human malaria, has one of the most A/T-biased genomes sequenced to date. This may give the genome and the transcriptome unusual structural features. Recent progress in sequencing techniques has made it possible to study the secondary structures of RNA molecules at the transcriptomic level. Thus, in this study we produced the in vivo RNA structurome of a protozoan parasite with a highly A/U-biased transcriptome. We showed that it is possible to probe the secondary structures of P. falciparum RNA molecules in vivo using two different chemical probes, and obtained structures for more than half of all transcripts in the transcriptome. These showed greater stability (lower free energy) than the same structures modelled in silico, and structural features appeared to influence translation efficiency and RNA decay. Finally, we compared the P. falciparum RNA structurome with the predicted RNA structurome of an A/U-balanced species, P. knowlesi, finding a bias towards lower overall transcript stability and more hairpins and multi-stem loops in P. falciparum. This unusual protozoan RNA structurome will provide a basis for similar studies in other protozoans and also in other unusual genomes.

5-methylcytosine modification by Plasmodium NSUN2 stabilizes mRNA and mediates the development of gametocytes

5-methylcytosine (m⁵C) is an important epitranscriptomic modification involved in messenger RNA (mRNA) stability and translation efficiency in various biological processes. However, it remains unclear if m⁵C modification contributes to the dynamic regulation of the transcriptome during the developmental cycles of Plasmodium parasites. Here, we characterize the landscape of m⁵C mRNA modifications at single nucleotide resolution in the asexual replication stages and gametocyte sexual stages of rodent (Plasmodium yoelii) and human (Plasmodium falciparum) malaria parasites. While different representations of m⁵C-modified mRNAs are associated with the different stages, the abundance of the m⁵C marker is strikingly enhanced in the transcriptomes of gametocytes. Our results show that m⁵C modifications confer stability to the Plasmodium transcripts and that a Plasmodium ortholog of NSUN2 is a major mRNA m⁵C methyltransferase in malaria parasites. Upon knockout of P. yoelii nsun2 (pynsun2), marked reductions of m⁵C modification were observed in a panel of gametocytogenesis-associated transcripts. These reductions correlated with impaired gametocyte production in the knockout rodent malaria parasites. Restoration of the nsun2 gene in the knockout parasites rescued the gametocyte production phenotype as well as m⁵C modification of the gametocytogenesis-associated transcripts. Together with the mRNA m⁵C profiles for two species of Plasmodium, our findings demonstrate a major role for NSUN2-mediated m⁵C modifications in mRNA transcript stability and sexual differentiation in malaria parasites.

Plasmodium falciparum GBP2 Is a Telomere-Associated Protein That Binds to G-Quadruplex DNA and RNA

In the early-diverging protozoan parasite Plasmodium, few telomere-binding proteins have been identified and several are unique. Plasmodium telomeres, like those of most eukaryotes, contain guanine-rich repeats that can form G-quadruplex structures. In model systems, quadruplex-binding drugs can disrupt telomere maintenance and some quadruplex-binding drugs are potent anti-plasmodial agents. Therefore, telomere-interacting and quadruplex-interacting proteins may offer new targets for anti-malarial therapy. Here, we report that P. falciparum GBP2 is such a protein. It was identified via 'Proteomics of Isolated Chromatin fragments', applied here for the first time in Plasmodium. In vitro, PfGBP2 binds specifically to G-rich telomere repeats in quadruplex form and it can also bind to G-rich RNA. In vivo, PfGBP2 partially colocalises with the known telomeric protein HP1 but is also found in the cytoplasm, probably due to its affinity for RNA. Consistently, its interactome includes numerous RNA-associated proteins. PfGBP2 is evidently a multifunctional DNA/RNA-binding factor in Plasmodium.

Contribution of Transcriptome to Elucidate the Biology of Plasmodium spp

In the present review, we discuss some of the new technologies that have been applied to elucidate how Plasmodium spp escape from the immune system and subvert the host physiology to orchestrate the regulation of its biological pathways. Our manuscript describes how techniques such as microarray approaches, RNA-Seq and single-cell RNA sequencing have contributed to the discovery of transcripts and changed the concept of gene expression regulation in closely related malaria parasite species. Moreover, the text highlights the contributions of high-throughput RNA sequencing for the current knowledge of malaria parasite biology, physiology, vaccine target and the revelation of new players in parasite signaling.

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.

Deep Learning Methods in Predicting Gene Expression Levels for the Malaria Parasite

Malaria is a mosquito-borne disease caused by single-celled blood parasites of the genus Plasmodium. The most severe cases of this disease are caused by the Plasmodium species, Falciparum. Once infected, a human host experiences symptoms of recurrent and intermittent fevers occurring over a time-frame of 48 hours, attributed to the synchronized developmental cycle of the parasite during the blood stage. To understand the regulated periodicity of Plasmodium falciparum transcription, this paper forecast and predict the P. falciparum gene transcription during its blood stage life cycle implementing a well-tuned recurrent neural network with gated recurrent units. Additionally, we also employ a spiking neural network to predict the expression levels of the P. falciparum gene. We provide results of this prediction on multiple genes including potential genes that express possible drug target enzymes. Our results show a high level of accuracy in being able to predict and forecast the expression levels of the different genes.

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.

The PfAP2-G2 transcription factor is a critical regulator of gametocyte maturation

Differentiation from asexual blood stages to mature sexual gametocytes is required for the transmission of malaria parasites. Here, we report that the ApiAP2 transcription factor, PfAP2-G2 (PF3D7_1408200) plays a critical role in the maturation of Plasmodium falciparum gametocytes. PfAP2-G2 binds to the promoters of a wide array of genes that are expressed at many stages of the parasite life cycle. Interestingly, we also find binding of PfAP2-G2 within the gene body of almost 3,000 genes, which strongly correlates with the location of H3K36me3 and several other histone modifications as well as Heterochromatin Protein 1 (HP1), suggesting that occupancy of PfAP2-G2 in gene bodies may serve as an alternative regulatory mechanism. Disruption of pfap2-g2 does not impact asexual development, but the majority of sexual parasites are unable to mature beyond stage III gametocytes. The absence of pfap2-g2 leads to overexpression of 28% of the genes bound by PfAP2-G2 and none of the PfAP2-G2 bound genes are downregulated, suggesting that it is a repressor. We also find that PfAP2-G2 interacts with chromatin remodeling proteins, a microrchidia (MORC) protein, and another ApiAP2 protein (PF3D7_1139300). Overall our data demonstrate that PfAP2-G2 establishes an essential gametocyte maturation program in association with other chromatin-related proteins.

Functional Characterization of the m6A-Dependent Translational Modulator PfYTH.2 in the Human Malaria Parasite

Posttranscriptional regulation of gene expression is central to the development and replication of the malaria parasite, Plasmodium falciparum, within its human host. The timely coordination of RNA maturation, homeostasis, and protein synthesis relies on the recruitment of specific RNA-binding proteins to their cognate target mRNAs. One possible mediator of such mRNA-protein interactions is the N⁶-methylation of adenosines (m⁶A), a prevalent mRNA modification of parasite mRNA transcripts. Here, we used RNA protein pulldowns, RNA modification mass spectrometry, and quantitative proteomics to identify two P. falciparum YTH domain proteins (PfYTH.1 and PfYTH.2) as m⁶A-binding proteins during parasite blood-stage development. Interaction proteomics revealed that PfYTH.2 associates with the translation machinery, including multiple subunits of the eukaryotic initiation factor 3 (eIF3) and poly(A)-binding proteins. Furthermore, knock sideways of PfYTH.2 coupled with ribosome profiling showed that this m⁶A reader is essential for parasite survival and is a repressor of mRNA translation. Together, these data reveal an important missing link in the m⁶A-mediated mechanism controlling mRNA translation in a unicellular eukaryotic pathogen.IMPORTANCE Infection with the unicellular eukaryotic pathogen Plasmodium falciparum causes malaria, a mosquito-borne disease affecting more than 200 million and killing 400,000 people each year. Underlying the asexual replication within human red blood cells is a tight regulatory network of gene expression and protein synthesis. A widespread mechanism of posttranscriptional gene regulation is the chemical modification of adenosines (m⁶A), through which the fate of individual mRNA transcripts can be changed. Here, we report on the protein machinery that "reads" this modification and "translates" it into a functional outcome. We provide mechanistic insight into one m⁶A reader protein and show that it interacts with the translational machinery and acts as a repressor of mRNA translation. This m⁶A-mediated phenotype has not been described in other eukaryotes as yet, and the functional characterization of the m⁶A interactome will ultimately open new avenues to combat the disease.

Plasmodium falciparum DDX17 is an RNA helicase crucial for parasite development

Malaria is one of the major global health concerns still prevailing in this 21st century. Even the effect of artemisinin combination therapies (ACT) have declined and causing more mortality across the globe. Therefore, it is important to understand the basic biology of malaria parasite in order to find novel drug targets. Helicases play important role in nucleic acid metabolism and are components of cellular machinery in various organisms. In this manuscript we have performed the biochemical characterization of homologue of DDX17 from Plasmodium falciparum (PfDDX17). Our results show that PfDDX17 is an active RNA helicase and uses mostly ATP for its function. The qRT-PCR experiment results suggest that PfDDX17 is highly expressed in the trophozoite stage and it is localised mainly in the cytoplasm and in infected RBC (iRBC) membrane mostly in the trophozoite stage. The dsRNA knockdown study suggests that PfDDX17 is important for cell cycle progression. These studies report the biochemical functions of PfDDX17 helicase and further augment the fundamental knowledge about helicase families of P. falciparum.

•

Biochemical characterization of homologue of DDX17 from Plasmodium falciparum (PfDDX17) is presented.

•

Results show that PfDDX17 is an active RNA helicase and uses mostly ATP for its function.

•

Results also suggest that PfDDX17 is highly expressed in the trophozoite stage.

•

dsRNA knockdown study revealed that PfDDX17 is important for cell cycle progression.

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.

Transcriptional heterogeneity and tightly regulated changes in gene expression during Plasmodium berghei sporozoite development

Despite the critical role of Plasmodium sporozoites in malaria transmission, we still know little about the mechanisms underlying their development in mosquitoes. Here, we use single-cell RNA sequencing to characterize the gene expression profiles of 16,038 Plasmodium berghei sporozoites isolated throughout their development from midgut oocysts to salivary glands, and from forced salivation experiments. Our results reveal a succession of tightly regulated changes in gene expression occurring during the maturation of sporozoites and highlight candidate genes that could play important roles in oocyst egress, sporozoite motility, and the mechanisms underlying the invasion of mosquito salivary glands and mammalian hepatocytes. In addition, the single-cell data reveal extensive transcriptional heterogeneity among parasites isolated from the same anatomical site, suggesting that Plasmodium development in mosquitoes is asynchronous and regulated by intrinsic as well as environmental factors. Finally, our analyses show a decrease in transcriptional activity preceding the translational repression observed in mature sporozoites and associated with their quiescent state in salivary glands, followed by a rapid reactivation of the transcriptional machinery immediately upon salivation.

Protein Modification Characteristics of the Malaria Parasite Plasmodium falciparum and the Infected Erythrocytes

Malaria elimination is still pending on the development of novel tools that rely on a deep understanding of parasite biology. Proteins of all living cells undergo myriad posttranslational modifications (PTMs) that are critical to multifarious life processes. An extensive proteome-wide dissection revealed a fine PTM map of most proteins in both Plasmodium falciparum, the causative agent of severe malaria, and the infected red blood cells. More than two-thirds of proteins of the parasite and its host cell underwent extensive and dynamic modification throughout the erythrocytic developmental stage. PTMs critically modulate the virulence factors involved in the host-parasite interaction and pathogenesis. Furthermore, P. falciparum stabilized the supporting proteins of erythrocyte origin by selective demodification. Collectively, our multiple omic analyses, apart from having furthered a deep understanding of the systems biology of P. falciparum and malaria pathogenesis, provide a valuable resource for mining new antimalarial targets.

Histone modifications associated with gene expression and genome accessibility are dynamically enriched at Plasmodium falciparum regulatory sequences

Background

The malaria parasite Plasmodium falciparum has an unusually euchromatic genome with poorly conserved positioning of nucleosomes in intergenic sequences and poorly understood mechanisms of gene regulation. Variant histones and histone modifications determine nucleosome stability and recruit trans factors, but their combinatorial contribution to gene regulation is unclear.

Results

Here, we show that the histone H3 acetylations H3K18ac and H3K27ac and the variant histone Pf H2A.Z are enriched together at regulatory sites upstream of genes. H3K18ac and H3K27ac together dynamically mark regulatory regions of genes expressed during the asexual life cycle. In contrast, H3K4me1 is depleted in intergenic sequence and dynamically depleted upstream of expressed genes. The temporal pattern of H3K27ac and H3K18ac enrichment indicates that they accumulate during S phase and mitosis and are retained at regulatory sequences until at least G1 phase and after cessation of expression of the cognate genes. We integrated our ChIPseq data with existing datasets to show that in schizont stages H3K18ac, H3K27ac and Pf H2A.Z colocalise with the transcription factor PfAP2-I and the bromodomain protein PfBDP1 and are enriched at stably positioned nucleosomes within regions of exposed DNA at active transcriptional start sites. Using transient transfections we showed that sequences enriched with colocalised H3K18ac, H3K27ac and Pf H2A.Z possess promoter activity in schizont stages, but no enhancer-like activity.

Conclusions

The dynamic H3 acetylations define P. falciparum regulatory sequences and contribute to gene activation. These findings expand the knowledge of the chromatin landscape that regulates gene expression in P. falciparum.

Definition of constitutive and stage-enriched promoters in the rodent malaria parasite, Plasmodium yoelii

Background

Well-defined promoters are essential elements for genetic studies in all organisms, and enable controlled expression of endogenous genes, transgene expression, and gene editing. Despite this, there is a paucity of defined promoters for the rodent-infectious malaria parasites. This is especially true for Plasmodium yoelii, which is often used to study the mosquito and liver stages of malarial infection, as well as host immune responses to infection.

Methods

Here six promoters were selected from across the parasite’s life cycle (clag-a, dynein heavy chain delta, lap4, trap, uis4, lisp2) that have been invoked in the literature as controlling their genes in a stage-specific manner. A minimal promoter length for the constitutive pybip promoter that confers strong expression levels was also determined, which is useful for expression of reporters and gene editing enzymes.

Results

Instead, it was observed that these promoters confer stage-enriched gene control, as some parasites also effectively use these promoters in other stages. Thus, when used alone, these promoters could complicate the interpretation of results obtained from promoter swaps, stage-targeted recombination, or gene editing experiments.

Conclusions

Together these data indicate that achieving stage-specific effects, such as gene editing, is likely best done using a two-component system with independent promoter activities overlapping only in the intended life cycle stage.

Decreased K13 Abundance Reduces Hemoglobin Catabolism and Proteotoxic Stress, Underpinning Artemisinin Resistance

Increased tolerance of Plasmodium falciparum to front-line artemisinin antimalarials (ARTs) is associated with mutations in Kelch13 (K13), although the precise role of K13 remains unclear. Here, we show that K13 mutations result in decreased expression of this protein, while mislocalization of K13 mimics resistance-conferring mutations, pinpointing partial loss of function of K13 as the relevant molecular event. K13-GFP is associated with ∼170 nm diameter doughnut-shaped structures at the parasite periphery, consistent with the location and dimensions of cytostomes. Moreover, the hemoglobin-peptide profile of ring-stage parasites is reduced when K13 is mislocalized. We developed a pulse-SILAC approach to quantify protein turnover and observe less disruption to protein turnover following ART exposure when K13 is mislocalized. Our findings suggest that K13 regulates digestive vacuole biogenesis and the uptake/degradation of hemoglobin and that ART resistance is mediated by a decrease in heme-dependent drug activation, less proteotoxicity, and increased survival of parasite ring stages.

PfKsgA1 functions as a transcription initiation factor and interacts with the N-terminal region of the mitochondrial RNA polymerase of Plasmodium falciparum

The small mitochondrial genome (mtDNA) of the malaria parasite is known to transcribe its genes polycistonically, although promoter element(s) have not yet been identified. An unusually large Plasmodium falciparum candidate mitochondrial phage-like RNA polymerase (PfmtRNAP) with an extended N-terminal region is encoded by the parasite nuclear genome. Using specific antibodies against the enzyme, we established that PfmtRNAP was targeted exclusively to the mitochondrion and interacted with mtDNA. Phylogenetic analysis showed that it is part of a separate apicomplexan clade. A search for PfmtRNAP-associated transcription initiation factors using sequence homology and in silico protein-protein interaction network analysis identified PfKsgA1. PfKsgA1 is a dual cytosol- and mitochondrion-targeted protein that functions as a small subunit rRNA dimethyltransferase in ribosome biogenesis. Chromatin immunoprecipitation showed that PfKsgA1 interacts with mtDNA, and in vivo crosslinking and pull-down experiments confirmed PfmtRNAP-PfKsgA1 interaction. The ability of PfKsgA1 to serve as a transcription initiation factor was demonstrated by complementation of yeast mitochondrial transcription factor Mtf1 function in Rpo41-driven in vitro transcription. Pull-down experiments using PfKsgA1 and PfmtRNAP domains indicated that the N-terminal region of PfmtRNAP interacts primarily with the PfKsgA1 C-terminal domain with some contacts being made with the linker and N-terminal domain of PfKsgA1. In the absence of full-length recombinant PfmtRNAP, solution structures of yeast mitochondrial RNA polymerase Rpo41 complexes with Mtf1 or PfKsgA1 were determined by small-angle X-ray scattering. Protein interaction interfaces thus identified matched with those reported earlier for Rpo41-Mtf1 interaction and overlaid with the PfmtRNAP-interfacing region identified experimentally for PfKsgA1. Our results indicate that in addition to a role in mitochondrial ribosome biogenesis, PfKsgA1 has an independent function as a transcription initiation factor for PfmtRNAP.

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

Background

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

Results

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

Conclusions

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

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.

Global Repertoire of Human Antibodies Against Plasmodium falciparum RIFINs, SURFINs, and STEVORs in a Malaria Exposed Population

Clinical immunity to malaria develops after repeated exposure to Plasmodium falciparum parasites. Broadly reactive antibodies against parasite antigens expressed on the surface of infected erythrocytes (variable surface antigens; VSAs) are candidates for anti-malaria therapeutics and vaccines. Among the VSAs, several RIFIN, STEVOR, and SURFIN family members have been demonstrated to be targets of naturally acquired immunity against malaria. For example, RIFIN family members are important ligands for opsonization of P. falciparum infected erythrocytes with specific immunoglobulins (IgG) acquiring broad protective reactivity. However, the global repertoire of human anti-VSAs IgG, its variation in children, and the key protective targets remain poorly understood. Here, we report wheat germ cell-free system-based production and serological profiling of a comprehensive library of A-RIFINs, B-RIFINs, STEVORs, and SURFINs derived from the P. falciparum 3D7 parasite strain. We observed that >98% of assayed proteins (n = 265) were immunogenic in malaria-exposed individuals in Uganda. The overall breadth of immune responses was significantly correlated with age but not with clinical malaria outcome among the study volunteers. However, children with high levels of antibodies to four RIFINs (PF3D7_0201000, PF3D7_1254500, PF3D7_1040600, PF3D7_1041100), STEVOR (PF3D7_0732000), and SURFIN 1.2 (PF3D7_0113600) had prospectively reduced the risk of developing febrile malaria, suggesting that the 5 antigens are important targets of protective immunity. Further studies on the significance of repeated exposure to malaria infection and maintenance of such high-level antibodies would contribute to a better understanding of susceptibility and naturally acquired immunity to malaria.

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.

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.

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.

An Endoplasmic Reticulum CREC Family Protein Regulates the Egress Proteolytic Cascade in Malaria Parasites

The endoplasmic reticulum (ER) is thought to play an essential role during egress of malaria parasites because the ER is assumed to be required for biogenesis and secretion of egress-related organelles. However, no proteins localized to the parasite ER have been shown to play a role in egress of malaria parasites. In this study, we generated conditional mutants of the Plasmodium falciparum endoplasmic reticulum-resident calcium-binding protein (PfERC), a member of the CREC family. Knockdown of the PfERC gene showed that this gene is essential for asexual growth of P. falciparum Analysis of the intraerythrocytic life cycle revealed that PfERC is essential for parasite egress but is not required for protein trafficking or calcium storage. We found that PfERC knockdown prevents the rupture of the parasitophorous vacuole membrane. This is because PfERC knockdown inhibited the proteolytic maturation of the subtilisin-like serine protease SUB1. Using double mutant parasites, we showed that PfERC is required for the proteolytic maturation of the essential aspartic protease plasmepsin X, which is required for SUB1 cleavage. Further, we showed that processing of substrates downstream of the proteolytic cascade is inhibited by PfERC knockdown. Thus, these data establish that the ER-resident CREC family protein PfERC is a key early regulator of the egress proteolytic cascade of malaria parasites.IMPORTANCE The divergent eukaryotic parasites that cause malaria grow and divide within a vacuole inside a host cell, which they have to break open once they finish cell division. The egress of daughter parasites requires the activation of a proteolytic cascade, and a subtilisin-like protease initiates a proteolytic cascade to break down the membranes blocking egress. It is assumed that the parasite endoplasmic reticulum plays a role in this process, but the proteins in this organelle required for egress remain unknown. We have identified an early ER-resident regulator essential for the maturation of the recently discovered aspartic protease in the egress proteolytic cascade, plasmepsin X, which is required for maturation of the subtilisin-like protease. Conditional loss of PfERC results in the formation of immature and inactive egress proteases that are unable to breakdown the vacuolar membrane barring release of daughter parasites.

Combining Stage Specificity and Metabolomic Profiling to Advance Antimalarial Drug Discovery

We report detailed susceptibility profiling of asexual blood stages of the malaria parasite Plasmodium falciparum to clinical and experimental antimalarials, combined with metabolomic fingerprinting. Results revealed a variety of stage-specific and metabolic profiles that differentiated the modes of action of clinical antimalarials including chloroquine, piperaquine, lumefantrine, and mefloquine, and identified late trophozoite-specific peak activity and stage-specific biphasic dose-responses for the mitochondrial inhibitors DSM265 and atovaquone. We also identified experimental antimalarials hitting previously unexplored druggable pathways as reflected by their unique stage specificity and/or metabolic profiles. These included several ring-active compounds, ones affecting hemoglobin catabolism through distinct pathways, and mitochondrial inhibitors with lower propensities for resistance than either DSM265 or atovaquone. This approach, also applicable to other microbes that undergo multiple differentiation steps, provides an effective tool to prioritize compounds for further development within the context of combination therapies.

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.

A network-based approach reveals novel invasion and Maurer's clefts-related proteins in Plasmodium falciparum

Malaria continues to be a major concern in developing countries despite continuous efforts to find a cure for the disease. Understanding the pathogenesis mechanism is necessary to identify more effective drug targets against malaria. Many years of experimental research have generated a large amount of data for the malarial parasite, Plasmodium falciparum. These data are useful to understand the importance of certain parasite proteins, but it often remains unclear how these proteins come together, interact with other proteins and carry out their function. Identification of all proteins involved in pathogenesis is an important step towards understanding the molecular mechanism of pathogenesis. In this study, dynamic stage-specific protein-protein interaction networks were created based on gene expression data during the parasite's intra-erythrocytic stages and static protein-protein interaction data. Using previously known proteins of a biological event as seed proteins, the random walk with restart (RWR) method was used on the dynamic protein-protein interaction networks to identify novel proteins related to that event. Two screening procedures namely, permutation test and GO enrichment test were performed to increase the reliability of the RWR predictions. The proposed method was first validated on Plasmodium falciparum proteins related to invasion, where it could reproduce the existing knowledge from a small set of seed proteins. It was then used to identify novel Maurer's clefts resident proteins, where it could identify 152 parasite proteins. We show that the current approach can annotate conserved proteins with unknown function. The predicted proteins can help build a mechanistic model for disease pathogenesis, which will be useful in identifying new drug targets.

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 molecular machinery of translational control in malaria parasites

Translational control regulates the levels of protein synthesized from its transcript and is key for the rapid adjustment of gene expression in response to environmental stimuli. The regulation of translation is of special importance for malaria parasites, which pass through a complex life cycle that includes various replication phases in the different organs of the human and mosquito hosts and a sexual reproduction phase in the mosquito midgut. In particular, the quiescent transmission stages rely on translational control to rapidly adapt to the new environment, once they switch over from the human to the mosquito and vice versa. Three control mechanisms are currently proposed in Plasmodium, (1) global regulation that acts on the translation initiation complex; (2) mRNA-specific regulation, involving cis control elements, mRNA-binding proteins and translational repressors; and (3) induced mRNA decay by the Ccr4-Not and the RNA exosome complex. The main molecules controlling translation are highly conserved in malaria parasites and an increasing number of studies shed light on the interwoven pathways leading to the up or downregulation of protein synthesis in the diverse plasmodial stages. We here highlight recent findings on translational control during life cycle progression of Plasmodium and discuss the molecules involved in regulating protein synthesis.

The Malaria Cell Atlas: Single parasite transcriptomes across the complete Plasmodium life cycle

Malaria parasites adopt a remarkable variety of morphological life stages as they transition through multiple mammalian host and mosquito vector environments. We profiled the single-cell transcriptomes of thousands of individual parasites, deriving the first high-resolution transcriptional atlas of the entire Plasmodium berghei life cycle. We then used our atlas to precisely define developmental stages of single cells from three different human malaria parasite species, including parasites isolated directly from infected individuals. The Malaria Cell Atlas provides both a comprehensive view of gene usage in a eukaryotic parasite and an open-access reference dataset for the study of malaria parasites.

Transcriptome-wide dynamics of extensive m6A mRNA methylation during Plasmodium falciparum blood-stage development

Malaria pathogenesis results from the asexual replication of Plasmodium falciparum within human red blood cells, which relies on a precisely timed cascade of gene expression over a 48-h life cycle. Although substantial post-transcriptional regulation of this hardwired program has been observed, it remains unclear how these processes are mediated on a transcriptome-wide level. To this end, we identified mRNA modifications in the P. falciparum transcriptome and performed a comprehensive characterization of N⁶-methyladenosine (m⁶A) over the course of blood-stage development. Using mass spectrometry and m⁶A RNA sequencing, we demonstrate that m⁶A is highly developmentally regulated, exceeding m⁶A levels known in any other eukaryote. We characterize a distinct m⁶A writer complex and show that knockdown of the putative m⁶A methyltransferase, PfMT-A70, by CRISPR interference leads to increased levels of transcripts that normally contain m⁶A. In accordance, we find an inverse correlation between m⁶A methylation and mRNA stability or translational efficiency. We further identify two putative m⁶A-binding YTH proteins that are likely to be involved in the regulation of these processes across the parasite's life cycle. Our data demonstrate unique features of an extensive m⁶A mRNA methylation programme in malaria parasites and reveal its crucial role in dynamically fine-tuning the transcriptional cascade of a unicellular eukaryote.

The endoplasmic reticulum chaperone PfGRP170 is essential for asexual development and is linked to stress response in malaria parasites

The vast majority of malaria mortality is attributed to one parasite species: Plasmodium falciparum. Asexual replication of the parasite within the red blood cell is responsible for the pathology of the disease. In Plasmodium, the endoplasmic reticulum (ER) is a central hub for protein folding and trafficking as well as stress response pathways. In this study, we tested the role of an uncharacterised ER protein, PfGRP170, in regulating these key functions by generating conditional mutants. Our data show that PfGRP170 localises to the ER and is essential for asexual growth, specifically required for proper development of schizonts. PfGRP170 is essential for surviving heat shock, suggesting a critical role in cellular stress response. The data demonstrate that PfGRP170 interacts with the Plasmodium orthologue of the ER chaperone, BiP. Finally, we found that loss of PfGRP170 function leads to the activation of the Plasmodium eIF2α kinase, PK4, suggesting a specific role for this protein in this parasite stress response pathway.

Plasmodium knowlesi as a model system for characterising Plasmodium vivax drug resistance candidate genes

Plasmodium vivax causes the majority of malaria outside Africa, but is poorly understood at a cellular level partly due to technical difficulties in maintaining it in in vitro culture conditions. In the past decades, drug resistant P. vivax parasites have emerged, mainly in Southeast Asia, but while some molecular markers of resistance have been identified, none have so far been confirmed experimentally, which limits interpretation of the markers, and hence our ability to monitor and control the spread of resistance. Some of these potential markers have been identified through P. vivax genome-wide population genetic analyses, which highlighted genes under recent evolutionary selection in Southeast Asia, where chloroquine resistance is most prevalent. These genes could be involved in drug resistance, but no experimental proof currently exists to support this hypothesis. In this study, we used Plasmodium knowlesi, the most closely related species to P. vivax that can be cultured in human erythrocytes, as a model system to express P. vivax genes and test for their role in drug resistance. We adopted a strategy of episomal expression, and were able to express fourteen P. vivax genes, including two allelic variants of several hypothetical resistance genes. Their expression level and localisation were assessed, confirming cellular locations conjectured from orthologous species, and suggesting locations for several previously unlocalised proteins, including an apical location for PVX_101445. These findings establish P. knowlesi as a suitable model for P. vivax protein expression. We performed chloroquine and mefloquine drug assays, finding no significant differences in drug sensitivity: these results could be due to technical issues, or could indicate that these genes are not actually involved in drug resistance, despite being under positive selection pressure in Southeast Asia. These data confirm that in vitro P. knowlesi is a useful tool for studying P. vivax biology. Its close evolutionary relationship to P. vivax, high transfection efficiency, and the availability of markers for colocalisation, all make it a powerful model system. Our study is the first of its kind using P. knowlesi to study unknown P. vivax proteins and investigate drug resistance mechanisms.

Progression of the canonical reference malaria parasite genome from 2002-2019

Here we describe the ways in which the sequence and annotation of the Plasmodium falciparum reference genome has changed since its publication in 2002. As the malaria species responsible for the most deaths worldwide, the richness of annotation and accuracy of the sequence are important resources for the P. falciparum research community as well as the basis for interpreting the genomes of subsequently sequenced species. At the time of publication in 2002 over 60% of predicted genes had unknown functions. As of March 2019, this number has been significantly decreased to 33%. The reduction is due to the inclusion of genes that were subsequently characterised experimentally and genes with significant similarity to others with known functions. In addition, the structural annotation of genes has been significantly refined; 27% of gene structures have been changed since 2002, comprising changes in exon-intron boundaries, addition or deletion of exons and the addition or deletion of genes. The sequence has also undergone significant improvements. In addition to the correction of a large number of single-base and insertion or deletion errors, a major miss-assembly between the subtelomeres of chromosome 7 and 8 has been corrected. As the number of sequenced isolates continues to grow rapidly, a single reference genome will not be an adequate basis for interpreting intra-species sequence diversity. We therefore describe in this publication a population reference genome of P. falciparum, called Pfref1. This reference will enable the community to map to regions that are not present in the current assembly. P. falciparum 3D7 will continue to be maintained, with ongoing curation ensuring continual improvements in annotation quality.

Progression of the canonical reference malaria parasite genome from 2002-2019

Here we describe the ways in which the sequence and annotation of the Plasmodium falciparum reference genome has changed since its publication in 2002. As the malaria species responsible for the most deaths worldwide, the richness of annotation and accuracy of the sequence are important resources for the P. falciparum research community as well as the basis for interpreting the genomes of subsequently sequenced species. At the time of publication in 2002 over 60% of predicted genes had unknown functions. As of March 2019, this number has been significantly decreased to 33%. The reduction is due to the inclusion of genes that were subsequently characterised experimentally and genes with significant similarity to others with known functions. In addition, the structural annotation of genes has been significantly refined; 27% of gene structures have been changed since 2002, comprising changes in exon-intron boundaries, addition or deletion of exons and the addition or deletion of genes. The sequence has also undergone significant improvements. In addition to the correction of a large number of single-base and insertion or deletion errors, a major miss-assembly between the subtelomeres of chromosome 7 and 8 has been corrected. As the number of sequenced isolates continues to grow rapidly, a single reference genome will not be an adequate basis for interpreting intra-species sequence diversity. We therefore describe in this publication a population reference genome of P. falciparum, called Pfref1. This reference will enable the community to map to regions that are not present in the current assembly. P. falciparum 3D7 will continue to be maintained, with ongoing curation ensuring continual improvements in annotation quality.

ApiAP2 Transcription Factors in Apicomplexan Parasites

Apicomplexan parasites are protozoan organisms that are characterised by complex life cycles and they include medically important species, such as the malaria parasite Plasmodium and the causative agents of toxoplasmosis (Toxoplasma gondii) and cryptosporidiosis (Cryptosporidium spp.). Apicomplexan parasites can infect one or more hosts, in which they differentiate into several morphologically and metabolically distinct life cycle stages. These developmental transitions rely on changes in gene expression. In the last few years, the important roles of different members of the ApiAP2 transcription factor family in regulating life cycle transitions and other aspects of parasite biology have become apparent. Here, we review recent progress in our understanding of the different members of the ApiAP2 transcription factor family in apicomplexan parasites.