A global view of the nonprotein-coding transcriptome in Plasmodium 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.

The role of long noncoding RNAs in malaria parasites

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Emerging biology of noncoding RNAs in malaria parasites

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

Non-coding RNAs in malaria infection

Malaria is one of the most severe infectious diseases affecting humans and it is caused by protozoan pathogens of the species Plasmodium (spp.). The malaria parasite Plasmodium is characterized by a complex, multistage life cycle that requires tight gene regulation which allows for host invasion and defense against host immune responses. Unfortunately, the mechanisms regulating gene expression during Plasmodium infection remain largely elusive, though several lines of evidence implicate a major involvement of non-coding RNAs (ncRNAs). The ncRNAs have been found to play a key role in regulating transcriptional and post-transcriptional events in a broad range of organisms including Plasmodium. In Plasmodium ncRNAs have been shown to regulate key events in the multistage life cycle and virulence ability. Here we review recent progress involving ncRNAs (microRNAs, long non-coding RNAs, and circular RNAs) and their role as regulators of gene expression during Plasmodium infection in human hosts with focus on the possibility of using these molecules as biomarkers for monitoring disease status. We also discuss the surprising function of ncRNAs in mediating the complex interplay between parasite and human host and future perspectives of the field. This article is categorized under: RNA in Disease and Development > RNA in Disease.

Phage Display Screening for Alba Superfamily Proteins from the Human Malaria Parasite, Plasmodium falciparum Reveals a High Level of Association with Protein Modification Pathways and Hints at New Drug Targets

PURPOSE

A 2016 study estimated that over 3 billion people are currently at risk of contracting malaria. Although a wide variety of medications are available to treat malaria, the parasites have started to exhibit resistance to many commonly used therapeutics necessitating a push for new investigations to identify novel drug targets.

METHODS

In this study, nucleic acid-binding Alba superfamily proteins of the human malaria parasite, Plasmodium falciparum were investigated to identify interacting protein motifs. A high-throughput molecular screening technique, phage display, coupled with next-generation sequencing was applied to assess large data sets.

RESULTS

Four P. falciparum Alba proteins were used for screening which appear to have distinct roles in parasite biology based on the results of this work. The majority of the peptide motifs identified from phage display were involved in post-translational modification pathways, thus suggesting that parasite-specific gene regulatory mechanisms are involved which could serve as drug targets for novel therapeutics.

CONCLUSION

This study found 18 peptide motifs which potentially have strong interactions with one or more of the Alba superfamily proteins from P. falciparum. Considering the large fraction of post-translational modification-related peptide motifs identified from this work, one or more of the protein modification pathways could serve as a good target for malaria treatment.

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

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

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.

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.

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.

Utilization of Small RNA Genes to Distinguish Vibrio cholerae Biotypes via Multiplex Polymerase Chain Reaction

The diarrheal disease "cholera" is caused by Vibrio cholerae, and is primarily confined to endemic regions, mostly in Africa and Asia. It is punctuated by outbreaks and creates severe challenges to public health. The disease-causing strains are most-often members of serogroups O1 and O139. PCR-based methods allow rapid diagnosis of these pathogens, including the identification of their biotypes. However, this necessitates the selection of specific target sequences to differentiate even the closely related biotypes of V. cholerae. Oligonucleotides for selective amplification of small RNA (sRNA) genes that are specific to these V. cholerae subtypes were designed. The resulting multiplex PCR assay was validated using V. cholerae cultures (i.e., 19 V. cholerae and 22 non-V. cholerae isolates) and spiked stool samples. The validation using V. cholerae cultures and spiked stool suspensions revealed detection limits of 10-100 pg DNA per reaction and 1.5 cells/mL suspension, respectively. The multiplex PCR assay that targets sRNA genes for amplification enables the sensitive and specific detection, as well as the differentiation of V. cholerae-O1 classical, O1 El Tor, and O139 biotypes. Most importantly, the assay enables fast and cheaper diagnosis compared with classic culture-based methods.

Post-Genomic Approaches to Understanding Malaria Parasite Biology: Linking Genes to Biological Functions

Plasmodium species are evolutionarily distant from model eukaryotes, and as a consequence they exhibit many non-canonical cellular processes. In the post-genomic era, functional "omics" disciplines (transcriptomics, proteomics, and metabolomics) have accelerated our understanding of unique aspects of the biology of malaria parasites. Functional "omics" tools, in combination with genetic manipulations, have offered new opportunities to investigate the function of previously uncharacterized genes. Knowledge of basic parasite biology is fundamental to understanding drug modes of action, mechanisms of drug resistance, and relevance of vaccine candidates. This Perspective highlights recent "omics"-based discoveries in basic biology and gene function of the most virulent human malaria parasite, Plasmodium falciparum.

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.

Unraveling the 3D genome of human malaria parasites

The chromosomes within the eukaryotic cell nucleus are highly dynamic and adopt complex hierarchical structures. Understanding how this three-dimensional (3D) nuclear architectureaffects gene regulation, cell cycle progression and disease pathogenesis are important biological questions in development and disease. Recently, many genome-wide technologies including chromosome conformation capture (3C) and 3C-based methodologies (4C, 5C, and Hi-C) have been developed to investigate 3D chromatin structure. In this review, we introduce 3D genome methodologies, with a focus on their application for understanding the nuclear architecture of the human malaria parasite, Plasmodium falciparum. An increasing amount of evidence now suggests that gene regulation in the parasite is largely regulated by epigenetic mechanisms and nuclear reorganization. Here, we explore the 3D genome architecture of P. falciparum, including local and global chromatin structure. In addition, molecular components important for maintaining 3D chromatin organization including architectural proteins and long non-coding RNAs are discussed. Collectively, these studies contribute to our understanding of how the plasticity of 3D genome architecture regulates gene expression and cell cycle progression in this deadly parasite.

Independent regulation of Plasmodium falciparum rif gene promoters

All Plasmodium species express variant antigens which may mediate immune escape in the vertebrate host. In Plasmodium falciparum, the rif gene family encodes variant antigens which are partly exposed on the infected red blood cell surface and may function as virulence factors. Not all rif genes are expressed at the same time and it is unclear what controls rif gene expression. In this work, we addressed global rif transcription using plasmid vectors with two drug resistance markers, one controlled by a rif 5′ upstream region and the second by a constitutively active promoter. After spontaneous integration into the genome of one construct, we observed that the resistance marker controlled by the rif 5′ upstream region was expressed dependent on the applied drug pressure. Then, the global transcription of rif genes in these transfectants was compared in the presence or absence of drugs. The relative transcript quantities of all rif loci did not change profoundly between strains grown with or without drug. We conclude that either there is no crosstalk between rif loci or that the elusive system of allelic exclusion of rif gene transcription is not controlled by their 5′ upstream region alone.

Pervasive Transcription of Mitochondrial, Plastid, and Nucleomorph Genomes across Diverse Plastid-Bearing Species

Organelle genomes exhibit remarkable diversity in content, structure, and size, and in their modes of gene expression, which are governed by both organelle- and nuclear-encoded machinery. Next generation sequencing (NGS) has generated unprecedented amounts of genomic and transcriptomic data, which can be used to investigate organelle genome transcription. However, most of the available eukaryotic RNA-sequencing (RNA-seq) data are used to study nuclear transcription only, even though large numbers of organelle-derived reads can typically be mined from these experiments. Here, we use publicly available RNA-seq data to assess organelle genome transcription in 59 diverse plastid-bearing species. Our RNA mapping analyses unraveled pervasive (full or near-full) transcription of mitochondrial, plastid, and nucleomorph genomes. In all cases, 85% or more of the organelle genome was recovered from the RNA data, including noncoding (intergenic and intronic) regions. These results reinforce the idea that organelles transcribe all or nearly all of their genomic material and are dependent on post-transcriptional processing of polycistronic transcripts. We explore the possibility that transcribed intergenic regions are producing functional noncoding RNAs, and that organelle genome noncoding content might provide raw material for generating regulatory RNAs.

RNomic identification and evaluation of npcTB_6715, a non-protein-coding RNA gene as a potential biomarker for the detection of Mycobacterium tuberculosis

Technological advances in RNA biology greatly improved transcriptome profiling during the last two decades. Besides the discovery of many small RNAs (sRNA) that are involved in the physiological and pathophysiological regulation of various cellular circuits, it becomes evident that the corresponding RNA genes might also serve as potential biomarkers to monitor the progression of disease and treatment. sRNA gene candidate npcTB_6715 was previously identified via experimental RNomic (unpublished data), and we report its application as potential biomarker for the detection of Mycobacterium tuberculosis (MTB) in patient samples. For proof of principle, we developed a multiplex PCR assay and report its validation with 500 clinical cultures, positive for Mycobacteria. The analysis revealed 98.9% sensitivity, 96.1% specificity, positive and negative predictive values of 98.6% and 96.8%, respectively. These results underscore the diagnostic value of the sRNA gene as diagnostic marker for the specific detection of MTB in clinical samples. Its successful application and the general ease of PCR‐based detection compared to standard bacterial culture techniques might be the first step towards ‘point‐of‐care’ diagnostics of Mycobacteria. To the best of our knowledge, this is the first time for the design of diagnostic applications based on sRNA genes, in Mycobacteria.

Epigenetic landscapes underlining global patterns of gene expression in the human malaria parasite, Plasmodium falciparum

The dynamic chromatin landscape displaying combinatorial complexity of the epigenome impacts gene expression that underlies many events of differentiation and cell cycle progression. In the past few years, epigenetic mechanisms have emerged as important processes involved in the tight gene regulation in malaria parasites, Plasmodium spp. Focusing predominantly on Plasmodium falciparum, the species associated with the most severe form of the disease, many advances have been made in our understanding of the interaction between transcriptional regulation and epigenetic mechanisms as the pivotal processes in regulating life cycle progression, host parasite interactions and parasite adaptation to the host environment. This review focuses on the epigenome and its effect on transcriptional regulation in P. falciparum, highlighting its unique, evolutionary diverse features.

The Role of Chromatin Structure in Gene Regulation of the Human Malaria Parasite

The human malaria parasite, Plasmodium falciparum, depends on a coordinated regulation of gene expression for development and propagation within the human host. Recent developments suggest that gene regulation in the parasite is largely controlled by epigenetic mechanisms. Here, we discuss recent advancements contributing to our understanding of the mechanisms controlling gene regulation in the parasite, including nucleosome landscape, histone modifications, and nuclear architecture. In addition, various processes involved in regulation of parasite-specific genes and gene families are examined. Finally, we address the use of epigenetic processes as targets for novel antimalarial therapies. Collectively, these topics highlight the unique biology of P. falciparum, and contribute to our understanding of mechanisms regulating gene expression in this deadly parasite.

Trans-acting GC-rich non-coding RNA at var expression site modulates gene counting in malaria parasite

Monoallelic expression of the var multigene family enables immune evasion of the malaria parasite Plasmodium falciparum in its human host. At a given time only a single member of the 60-member var gene family is expressed at a discrete perinuclear region called the 'var expression site'. However, the mechanism of var gene counting remains ill-defined. We hypothesize that activation factors associating specifically with the expression site play a key role in this process. Here, we investigate the role of a GC-rich non-coding RNA (ncRNA) gene family composed of 15 highly homologous members. GC-rich genes are positioned adjacent to var genes in chromosome-central gene clusters but are absent near subtelomeric var genes. Fluorescence in situ hybridization demonstrates that GC-rich ncRNA localizes to the perinuclear expression site of central and subtelomeric var genes in trans. Importantly, overexpression of distinct GC-rich ncRNA members disrupts the gene counting process at the single cell level and results in activation of a specific subset of var genes in distinct clones. We identify the first trans-acting factor targeted to the elusive perinuclear var expression site and open up new avenues to investigate ncRNA function in antigenic variation of malaria and other protozoan pathogens.

Translational regulation in blood stages of the malaria parasite Plasmodium spp.: systems-wide studies pave the way

The malaria parasite Plasmodium spp. varies the expression profile of its genes depending on the host it resides in and its developmental stage. Virtually all messenger RNA (mRNA) is expressed in a monocistronic manner, with transcriptional activation regulated at the epigenetic level and by specialized transcription factors. Furthermore, recent systems-wide studies have identified distinct mechanisms of post-transcriptional and translational control at various points of the parasite lifecycle. Taken together, it is evident that 'just-in-time' transcription and translation strategies coexist and coordinate protein expression during Plasmodium development, some of which we review here. In particular, we discuss global and specific mechanisms that control protein translation in blood stages of the human malaria parasite Plasmodium falciparum, once a cytoplasmic mRNA has been generated, and its crosstalk with mRNA decay and storage. We also focus on the widespread translational delay observed during the 48-hour blood stage lifecycle of P. falciparum-for over 30% of transcribed genes, including virulence factors required to invade erythrocytes-and its regulation by cis-elements in the mRNA, RNA-processing enzymes and RNA-binding proteins; the first-characterized amongst these are the DNA- and RNA-binding Alba proteins. More generally, we conclude that translational regulation is an emerging research field in malaria parasites and propose that its elucidation will not only shed light on the complex developmental program of this parasite, but may also reveal mechanisms contributing to drug resistance and define new targets for malaria intervention strategies. WIREs RNA 2016, 7:772-792. doi: 10.1002/wrna.1365 For further resources related to this article, please visit the WIREs website.

Landscape and Dynamics of Transcription Initiation in the Malaria Parasite Plasmodium falciparum

A comprehensive map of transcription start sites (TSSs) across the highly AT-rich genome of P. falciparum would aid progress toward deciphering the molecular mechanisms that underlie the timely regulation of gene expression in this malaria parasite. Using high-throughput sequencing technologies, we generated a comprehensive atlas of transcription initiation events at single-nucleotide resolution during the parasite intra-erythrocytic developmental cycle. This detailed analysis of TSS usage enabled us to define architectural features of plasmodial promoters. We demonstrate that TSS selection and strength are constrained by local nucleotide composition. Furthermore, we provide evidence for coordinate and stage-specific TSS usage from distinct sites within the same transcription unit, thereby producing transcript isoforms, a subset of which are developmentally regulated. This work offers a framework for further investigations into the interactions between genomic sequences and regulatory factors governing the complex transcriptional program of this major human pathogen.

Analysis of nucleosome positioning landscapes enables gene discovery in the human malaria parasite Plasmodium falciparum

Background

Plasmodium falciparum, the deadliest malaria-causing parasite, has an extremely AT-rich (80.7 %) genome. Because of high AT-content, sequence-based annotation of genes and functional elements remains challenging. In order to better understand the regulatory network controlling gene expression in the parasite, a more complete genome annotation as well as analysis tools adapted for AT-rich genomes are needed. Recent studies on genome-wide nucleosome positioning in eukaryotes have shown that nucleosome landscapes exhibit regular characteristic patterns at the 5’- and 3’-end of protein and non-protein coding genes. In addition, nucleosome depleted regions can be found near transcription start sites. These unique nucleosome landscape patterns may be exploited for the identification of novel genes. In this paper, we propose a computational approach to discover novel putative genes based exclusively on nucleosome positioning data in the AT-rich genome of P. falciparum.

Results

Using binary classifiers trained on nucleosome landscapes at the gene boundaries from two independent nucleosome positioning data sets, we were able to detect a total of 231 regions containing putative genes in the genome of Plasmodium falciparum, of which 67 highly confident genes were found in both data sets. Eighty-eight of these 231 newly predicted genes exhibited transcription signal in RNA-Seq data, indicative of active transcription. In addition, 20 out of 21 selected gene candidates were further validated by RT-PCR, and 28 out of the 231 genes showed significant matches using BLASTN against an expressed sequence tag (EST) database. Furthermore, 108 (47 %) out of the 231 putative novel genes overlapped with previously identified but unannotated long non-coding RNAs. Collectively, these results provide experimental validation for 163 predicted genes (70.6 %). Finally, 73 out of 231 genes were found to be potentially translated based on their signal in polysome-associated RNA-Seq representing transcripts that are actively being translated.

Conclusion

Our results clearly indicate that nucleosome positioning data contains sufficient information for novel gene discovery. As distinct nucleosome landscapes around genes are found in many other eukaryotic organisms, this methodology could be used to characterize the transcriptome of any organism, especially when coupled with other DNA-based gene finding and experimental methods (e.g., RNA-Seq).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2214-9) contains supplementary material, which is available to authorized users.

Challenges of drug-resistant malaria

Over the past six decades, the drug resistance of Plasmodium falciparum has become an issue of utmost concern. Despite the remarkable progress that has been made in recent years in reducing the mortality rate to about 30% with the scaling-up of vector control, introduction of artemisinin-based combination therapies and other malaria control strategies, the confirmation of artemisinin resistance on the Cambodia–Thailand border threatened all the previous success. This review addresses the global scenario of antimalarial resistance and factors associated with it, with the main emphasis on futuristic approaches like nanotechnology and stem cell therapy that may impede resistant malaria, along with novel medications which are preparing to enter the global antimalarial market. These novel studies are likely to escalate over the coming years and will hopefully help to reduce the burden of malaria.

Alternative processing as evolutionary mechanism for the origin of novel nonprotein coding RNAs

The evolution of new genes can ensue through either gene duplication and the neofunctionalization of one of the copies or the formation of a de novo gene from hitherto nonfunctional, neutrally evolving intergenic or intronic genomic sequences. Only very rarely are entire genes created de novo. Mostly, nonfunctional sequences are coopted as novel parts of existing genes, such as in the process of exonization whereby introns become exons through changes in splicing. Here, we report a case in which a novel nonprotein coding RNA evolved by intron-sequence recruitment into its structure. cDNAs derived from rat brain small RNAs, revealed a novel small nucleolar RNA (snoRNA) originating from one of the Snord115 copies in the rat Prader–Willi syndrome locus. We suggest that a single-point substitution in the Snord115 region led to the expression of a longer snoRNA variant, designated as L-Snord115. Cell culture and footprinting experiments confirmed that a single nucleotide substitution at Snord115 position 67 destabilized the kink-turn motif within the canonical snoRNA, while distal intronic sequences provided an alternate D-box region. The exapted sequence displays putative base pairing to 28S rRNA and mRNA targets.

Noncoding RNAs as emerging regulators of Plasmodium falciparum virulence gene expression

The eukaryotic unicellular pathogen Plasmodium falciparum tightly regulates gene expression, both during development and in adaptation to dynamic host environments. This regulation is evident in the mutually exclusive expression of members of clonally variant virulence multigene families. While epigenetic regulators have been selectively identified at active or repressed virulence genes, their specific recruitment remains a mystery. In recent years, noncoding RNAs (ncRNAs) have emerged as lynchpins of eukaryotic gene regulation; by binding to epigenetic regulators, they provide target specificity to otherwise non-specific enzyme complexes. Not surprisingly, there is great interest in understanding the role of ncRNA in P. falciparum, in particular, their contribution to the mutually exclusive expression of virulence genes. The current repertoire of P. falciparum ncRNAs includes, but is not limited to, subtelomeric ncRNAs, virulence gene-associated ncRNAs and natural antisense RNA transcripts. Continued improvement in high-throughput sequencing methods is sure to expand this repertoire. Here, we summarize recent advances in P. falciparum ncRNA biology, with an emphasis on ncRNA-mediated epigenetic modes of gene regulation.

Dataset of natural antisense transcripts in P. vivax clinical isolates derived using custom designed strand-specific microarray

Natural antisense transcripts (NATs) have been detected in many organisms and shown to regulate gene expression. Similarly, NATs have also been observed in malaria parasites with most studies focused on Plasmodium falciparum. There were no reports on the presence of NATs in Plasmodium vivax, which has also been shown to cause severe malaria like P. falciparum, until a recent study published by us. To identify in vivo prevalence of antisense transcripts in P. vivax clinical isolates, we performed whole genome expression profiling using a custom designed strand-specific microarray that contains probes for both sense and antisense strands. Here we describe the experimental methods and analysis of the microarray data available in Gene Expression Omnibus (GEO) under GSE45165. Our data provides a resource for exploring the presence of antisense transcripts in P. vivax isolated from patients showing varying clinical symptoms. Related information about the description and interpretation of the data can be found in a recent publication by Boopathi and colleagues in Infection, Genetics and Evolution 2013.

Deep profiling of the novel intermediate-size noncoding RNAs in intraerythrocytic Plasmodium falciparum

Intermediate-size noncoding RNAs (is-ncRNAs) have been shown to play important regulatory roles in the development of several eukaryotic organisms. However, they have not been thoroughly explored in Plasmodium falciparum, which is the most virulent malaria parasite infecting human being. By using Illumina/Solexa paired-end sequencing of an is-ncRNA-specific library, we performed a systematic identification of novel is-ncRNAs in intraerythrocytic P. falciparum, strain 3D7. A total of 1,198 novel is-ncRNA candidates, including antisense, intergenic, and intronic is-ncRNAs, were identified. Bioinformatics analyses showed that the intergenic is-ncRNAs were the least conserved among different Plasmodium species, and antisense is-ncRNAs were more conserved than their sense counterparts. Twenty-two novel snoRNAs were identified, and eight potential novel classes of P. falciparum is-ncRNAs were revealed by clustering analysis. The expression of randomly selected novel is-ncRNAs was confirmed by RT-PCR and northern blotting assays. An obvious different expressional profile of the novel is-ncRNA between the early and late intraerythrocytic developmental stages of the parasite was observed. The expression levels of the antisense RNAs correlated with those of their cis-encoded sense RNA counterparts, suggesting that these is-ncRNAs are involved in the regulation of gene expression of the parasite. In conclusion, we accomplished a deep profiling analysis of novel is-ncRNAs in P. falciparum, analysed the conservation and structural features of these novel is-ncRNAs, and revealed their differential expression patterns during the development of the parasite. These findings provide important information for further functional characterisation of novel is-ncRNAs during the development of P. falciparum.

Natural antisense transcripts in Plasmodium falciparum isolates from patients with complicated malaria

Mechanisms regulating gene expression in malaria parasites are not well understood. Little is known about how the parasite regulates its gene expression during transition from one developmental stage to another and in response to various environmental conditions. Parasites in a diseased host face environments which differ from the static, well adapted in vitro conditions. Parasites thus need to adapt quickly and effectively to these conditions by establishing transcriptional states which are best suited for better survival. With the discovery of natural antisense transcripts (NATs) in this parasite and considering the various proposed mechanisms by which NATs might regulate gene expression, it has been speculated that these might be playing a critical role in gene regulation. We report here the diversity of NATs in this parasite, using isolates taken directly from patients with differing clinical symptoms caused by malaria infection. Using a custom designed strand specific whole genome microarray, a total of 797 NATs targeted against annotated loci have been detected. Out of these, 545 NATs are unique to this study. The majority of NATs were positively correlated with the expression pattern of the sense transcript. However, 96 genes showed a change in sense/antisense ratio on comparison between uncomplicated and complicated disease conditions. The antisense transcripts map to a broad range of biochemical/metabolic pathways, especially pathways pertaining to the central carbon metabolism and stress related pathways. Our data strongly suggests that a large group of NATs detected here are unannotated transcription units antisense to annotated gene models. The results reveal a previously unknown set of NATs that prevails in this parasite, their differential regulation in disease conditions and mapping to functionally well annotated genes. The results detailed here call for studies to deduce the possible mechanism of action of NATs, which would further help in understanding the in vivo pathological adaptations of these parasites.

Genome-wide identification and functional annotation of Plasmodium falciparum long noncoding RNAs from RNA-seq data

The life cycle of Plasmodium falciparum is very complex, with an erythrocytic stage that involves the invasion of red blood cells and the survival and growth of the parasite within the host. Over the past several decades, numbers of studies have shown that proteins exported by P. falciparum to the surface of infected red blood cells play a critical role in recognition and interaction with host receptors and are thus essential for the completion of the life cycle of P. falciparum. However, little is known about long noncoding RNAs (lncRNAs). In this study, we designed a computational pipeline to identify new lncRNAs of P. falciparum from published RNA-seq data and analyzed their sequences and expression features. As a result, 164 novel lncRNAs were found. The sequences and expression features of P. falciparum lncRNAs were similar to those of humans and mice: there was a lack of sequence conservation, low expression levels, and high expression coefficient of variance and co-expression with nearby coding sequences in the genome. Next, a coding/noncoding gene co-expression network for P. falciparum was constructed to further annotate the functions of novel and known lncRNAs. In total, the functions of 69 lncRNAs, including 44 novel lncRNAs, were annotated. The main functions of the lncRNAs included metabolic processes, biosynthetic processes, regulation of biological processes, establishment of localization, catabolic processes, cellular component organization, and interspecies interactions between organisms. Our results will provide clues to further the investigation of interactions between human hosts and parasites and the mechanisms of P. falciparum infection.

Parasite epigenetics and immune evasion: lessons from budding yeast

The remarkable ability of many parasites to evade host immunity is the key to their success and pervasiveness. The immune evasion is directly linked to the silencing of the members of extended families of genes that encode for major parasite antigens. At any time only one of these genes is active. Infrequent switches to other members of the gene family help the parasites elude the immune system and cause prolonged maladies. For most pathogens, the detailed mechanisms of gene silencing and switching are poorly understood. On the other hand, studies in the budding yeast Saccharomyces cerevisiae have revealed similar mechanisms of gene repression and switching and have provided significant insights into the molecular basis of these phenomena. This information is becoming increasingly relevant to the genetics of the parasites. Here we summarize recent advances in parasite epigenetics and emphasize the similarities between S. cerevisiae and pathogens such as Plasmodium, Trypanosoma, Candida, and Pneumocystis. We also outline current challenges in the control and the treatment of the diseases caused by these parasites and link them to epigenetics and the wealth of knowledge acquired from budding yeast.

Biases in small RNA deep sequencing data

High-throughput RNA sequencing (RNA-seq) is considered a powerful tool for novel gene discovery and fine-tuned transcriptional profiling. The digital nature of RNA-seq is also believed to simplify meta-analysis and to reduce background noise associated with hybridization-based approaches. The development of multiplex sequencing enables efficient and economic parallel analysis of gene expression. In addition, RNA-seq is of particular value when low RNA expression or modest changes between samples are monitored. However, recent data uncovered severe bias in the sequencing of small non-protein coding RNA (small RNA-seq or sRNA-seq), such that the expression levels of some RNAs appeared to be artificially enhanced and others diminished or even undetectable. The use of different adapters and barcodes during ligation as well as complex RNA structures and modifications drastically influence cDNA synthesis efficacies and exemplify sources of bias in deep sequencing. In addition, variable specific RNA G/C-content is associated with unequal polymerase chain reaction amplification efficiencies. Given the central importance of RNA-seq to molecular biology and personalized medicine, we review recent findings that challenge small non-protein coding RNA-seq data and suggest approaches and precautions to overcome or minimize bias.

Diversity of mitochondrial genome structure in the phylum Apicomplexa

Mitochondria are ubiquitous organelles in all eukaryotes that are essential for a range of cellular processes and cellular signaling. Nearly all mitochondria have their own DNA or mitochondrial (mt) genome, which varies considerably in size, structure and organization. The phylum Apicomplexa includes a variety of unicellular eukaryotes, some of which are parasites of clinical or economic importance. Recent studies have demonstrated that apicomplexan mt genomes, which include the smallest 6 kb genome of the malaria parasites, exhibit remarkably diverse structures. Apicomplexan parasites are interesting model organisms in order to understand the evolution of mt genomes. This review summarizes the structure of apicomplexan mt genomes and highlights the unique features and the evolution of the mt genome.

microRNAs in parasites and parasite infection

miRNAs, a subclass of small regulatory RNAs, are present from ancient unicellular protozoans to parasitic helminths and parasitic arthropods. The miRNA-silencing mechanism appears, however, to be absent in a number of protozoan parasites. Protozoan miRNAs and components of their silencing machinery possess features different from other eukaryotes, providing some clues on the evolution of the RNA-induced silencing machinery. miRNA functions possibly associate with neoblast biology, development, physiology, infection and immunity of parasites. Parasite infection can alter host miRNA expression that can favor both parasite clearance and infection. miRNA pathways are, thus, a potential target for the therapeutic control of parasitic diseases.

Impact of chromosome ends on the biology and virulence of Plasmodium falciparum

In recent years, many studies have focused on heterochromatin located at chromosome ends, which plays an important role in regulating gene expression in many organisms ranging from yeast to humans. Similarly, in the protozoan Plasmodium falciparum, which is the most virulent human malaria parasite, the heterochromatin present in telomeres and subtelomeric regions exerts a silencing effect on the virulence gene families located therein. Studies addressing P. falciparum chromosome ends have demonstrated that these regions participate in other functions, such as the formation of the T-loop structure, the replication of telomeric regions, the regulation of telomere length and the formation of telomeric heterochromatin. In addition, telomeres are involved in anchoring chromosome ends to the nuclear periphery, thereby playing an important role in nuclear architecture and gene expression regulation. Here, we review the current understanding of chromosome ends, the proteins that bind to these regions and their impact on the biology and virulence of P. falciparum.

Computational prediction and validation of C/D, H/ACA and Eh_U3 snoRNAs of Entamoeba histolytica

BACKGROUND

Small nucleolar RNAs are a highly conserved group of small RNAs found in eukaryotic cells. Genes encoding these RNAs are diversely located throughout the genome. They are functionally conserved, performing post transcriptional modification (methylation and pseudouridylation) of rRNA and other nuclear RNAs. They belong to two major categories: the C/D box and H/ACA box containing snoRNAs. U3 snoRNA is an exceptional member of C/D box snoRNAs and is involved in early processing of pre-rRNA. An antisense sequence is present in each snoRNA which guides the modification or processing of target RNA. However, some snoRNAs lack this sequence and often they are called orphan snoRNAs.

RESULTS

We have searched snoRNAs of Entamoeba histolytica from the genome sequence using computational programmes (snoscan and snoSeeker) and we obtained 99 snoRNAs (C/D and H/ACA box snoRNAs) along with 5 copies of Eh_U3 snoRNAs. These are located diversely in the genome, mostly in intergenic regions, while some are found in ORFs of protein coding genes, intron and UTRs. The computationally predicted snoRNAs were validated by RT-PCR and northern blotting. The expected sizes were in agreement with the observed sizes for all C/D box snoRNAs tested, while for some of the H/ACA box there was indication of processing to generate shorter products.

CONCLUSION

Our results showed the presence of snoRNAs in E. histolytica, an early branching eukaryote, and the structural features of E. histolytica snoRNAs were well conserved when compared with yeast and human snoRNAs. This study will help in understanding the evolution of these conserved RNAs in diverse phylogenetic groups.

The fragmented mitochondrial ribosomal RNAs of Plasmodium falciparum

Background

The mitochondrial genome in the human malaria parasite Plasmodium falciparum is most unusual. Over half the genome is composed of the genes for three classic mitochondrial proteins: cytochrome oxidase subunits I and III and apocytochrome b. The remainder encodes numerous small RNAs, ranging in size from 23 to 190 nt. Previous analysis revealed that some of these transcripts have significant sequence identity with highly conserved regions of large and small subunit rRNAs, and can form the expected secondary structures. However, these rRNA fragments are not encoded in linear order; instead, they are intermixed with one another and the protein coding genes, and are coded on both strands of the genome. This unorthodox arrangement hindered the identification of transcripts corresponding to other regions of rRNA that are highly conserved and/or are known to participate directly in protein synthesis.

Principal Findings

The identification of 14 additional small mitochondrial transcripts from P. falcipaurm and the assignment of 27 small RNAs (12 SSU RNAs totaling 804 nt, 15 LSU RNAs totaling 1233 nt) to specific regions of rRNA are supported by multiple lines of evidence. The regions now represented are highly similar to those of the small but contiguous mitochondrial rRNAs of Caenorhabditis elegans. The P. falciparum rRNA fragments cluster on the interfaces of the two ribosomal subunits in the three-dimensional structure of the ribosome.

Significance

All of the rRNA fragments are now presumed to have been identified with experimental methods, and nearly all of these have been mapped onto the SSU and LSU rRNAs. Conversely, all regions of the rRNAs that are known to be directly associated with protein synthesis have been identified in the P. falciparum mitochondrial genome and RNA transcripts. The fragmentation of the rRNA in the P. falciparum mitochondrion is the most extreme example of any rRNA fragmentation discovered.

Two long non-coding RNAs generated from subtelomeric regions accumulate in a novel perinuclear compartment in Plasmodium falciparum

Chromosome ends have been implicated in the default silencing of clonally variant gene families in the human malaria parasite Plasmodium falciparum. These chromosome regions are organized into heterochromatin, as defined by the presence of a repressive histone H3 lysine 9 trimethylated marker and heterochromatin protein 1. Here, we show that the non-coding subtelomeric region adjacent to virulence genes forms facultative heterochromatin in a cell cycle-dependent manner. We demonstrate that telomere-associated repeat elements (TAREs) and telomeres are transcribed as long non-coding RNAs (lncRNAs) during schizogony. Northern blot assays revealed two classes of lncRNAs: a ~4-kb transcript composed of telomere sequences and a TARE-3 element, and a >6-kb transcript composed of 21-bp repeats from TARE-6. These lncRNAs are transcribed by RNA polymerase II as single-stranded molecules. RNA-FISH analysis showed that these lncRNAs form several nuclear foci during the schizont stage, whereas in the ring stage, they are located in a single perinuclear compartment that does not co-localize with any known nuclear subcompartment. Furthermore, the TARE-6 lncRNA is predicted to form a stable and repetitive hairpin structure that is able to bind histones. Consequently, the characterization of the molecular interactions of these lncRNAs with nuclear proteins may reveal novel modes of gene regulation and nuclear function in P. falciparum.

New Agilent platform DNA microarrays for transcriptome analysis of Plasmodium falciparum and Plasmodium berghei for the malaria research community

Background

DNA microarrays have been a valuable tool in malaria research for over a decade but remain in limited use in part due their relatively high cost, poor availability, and technical difficulty. With the aim of alleviating some of these factors next-generation DNA microarrays for genome-wide transcriptome analysis for both Plasmodium falciparum and Plasmodium berghei using the Agilent 8x15K platform were designed.

Methods

Probe design was adapted from previously published methods and based on the most current transcript predictions available at the time for P. falciparum or P. berghei. Array performance and transcriptome analysis was determined using dye-coupled, aminoallyl-labelled cDNA and streamlined methods for hybridization, washing, and array analysis were developed.

Results

The new array design marks a notable improvement in the number of transcripts covered and average number of probes per transcript. Array performance was excellent across a wide range of transcript abundance, with low inter-array and inter-probe variability for relative abundance measurements and it recapitulated previously observed transcriptional patterns. Additionally, improvements in sensitivity permitted a 20-fold reduction in necessary starting RNA amounts, further reducing experimental costs and widening the range of application.

Conclusions

DNA microarrays utilizing the Agilent 8x15K platform for genome-wide transcript analysis in P. falciparum and P. berghei mark an improvement in coverage and sensitivity, increased availability to the research community, and simplification of the experimental methods.

Genomics and integrated systems biology in Plasmodium falciparum: a path to malaria control and eradication

The first draft of the human malaria parasite's genome was released in 2002. Since then, the malaria scientific community has witnessed a steady embrace of new and powerful functional genomic studies. Over the years, these approaches have slowly revolutionized malaria research and enabled the comprehensive, unbiased investigation of various aspects of the parasite's biology. These genome-wide analyses delivered a refined annotation of the parasite's genome, delivered a better knowledge of its RNA, proteins and metabolite derivatives, and fostered the discovery of new vaccine and drug targets. Despite the positive impacts of these genomic studies, most research and investment still focus on protein targets, drugs and vaccine candidates that were known before the publication of the parasite genome sequence. However, recent access to next-generation sequencing technologies, along with an increased number of genome-wide applications, is expanding the impact of the parasite genome on biomedical research, contributing to a paradigm shift in research activities that may possibly lead to new optimized diagnosis and treatments. This review provides an update of Plasmodium falciparum genome sequences and an overview of the rapid development of genomics and system biology applications that have an immense potential of creating powerful tools for a successful malaria eradication campaign.

Harnessing genomics and genome biology to understand malaria biology

Malaria is an important human disease and is the target of a global eradication campaign. New technological and informatics advancements in population genomics are being leveraged to identify genetic loci under selection in the malaria parasite and to find variants that are associated with key clinical phenotypes, such as drug resistance. This article provides a timely Review of how population-genetics-based strategies are being applied to Plasmodium falciparum both to identify genetic loci as key targets of interventions and to develop monitoring and surveillance tools that are crucial for the successful elimination and eradication of malaria.