Strand-specific RNA-Seq reveals widespread and developmentally regulated transcription of natural antisense transcripts in Plasmodium falciparum

Time to switch gears: how long noncoding RNAs function as epigenetic regulators in Apicomplexan parasites

Long noncoding RNAs (lncRNA) are emerging as important regulators of gene expression in eukaryotes. In recent years, a large repertoire of lncRNA were discovered in Apicomplexan parasites and were implicated in several mechanisms of gene expression, including marking genes for activation, contributing to the formation of subnuclear compartments and organization, regulating the deposition of epigenetic modifications, influencing chromatin and chromosomal structure and manipulating host gene expression. Here, we aim to update recent knowledge on the role of lncRNAs as regulators in Apicomplexan parasites and highlight the possible molecular mechanisms by which they function. We hope that some of the hypotheses raised here will contribute to further investigation and lead to new mechanistic insight and better understanding of the role of lncRNA in parasite's biology.

Positive-unlabeled learning identifies vaccine candidate antigens in the malaria parasite Plasmodium falciparum

Malaria vaccine development is hampered by extensive antigenic variation and complex life stages of Plasmodium species. Vaccine development has focused on a small number of antigens, many of which were identified without utilizing systematic genome-level approaches. In this study, we implement a machine learning-based reverse vaccinology approach to predict potential new malaria vaccine candidate antigens. We assemble and analyze P. falciparum proteomic, structural, functional, immunological, genomic, and transcriptomic data, and use positive-unlabeled learning to predict potential antigens based on the properties of known antigens and remaining proteins. We prioritize candidate antigens based on model performance on reference antigens with different genetic diversity and quantify the protein properties that contribute most to identifying top candidates. Candidate antigens are characterized by gene essentiality, gene ontology, and gene expression in different life stages to inform future vaccine development. This approach provides a framework for identifying and prioritizing candidate vaccine antigens for a broad range of pathogens.

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

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

Comparative spatial proteomics of Plasmodium-infected erythrocytes

Plasmodium parasites contribute to one of the highest global infectious disease burdens. To achieve this success, the parasite has evolved a range of specialized subcellular compartments to extensively remodel the host cell for its survival. The information to fully understand these compartments is likely hidden in the so far poorly characterized Plasmodium species spatial proteome. To address this question, we determined the steady-state subcellular location of more than 12,000 parasite proteins across five different species by extensive subcellular fractionation of erythrocytes infected by Plasmodium falciparum, Plasmodium knowlesi, Plasmodium yoelii, Plasmodium berghei, and Plasmodium chabaudi. This comparison of the pan-species spatial proteomes and their expression patterns indicates increasing species-specific proteins associated with the more external compartments, supporting host adaptations and post-transcriptional regulation. The spatial proteome offers comprehensive insight into the different human, simian, and rodent Plasmodium species, establishing a powerful resource for understanding species-specific host adaptation processes in the parasite.

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.

Exosome-derived long noncoding RNAs: Mediators of host-Plasmodium parasite communication

Overcoming challenges associated with malaria eradication proves to be a formidable task due to the complicated life cycle exhibited by the malaria parasite and the lack of safe and enduring vaccines against malaria. Investigating the interplay between Plasmodium parasites and their mammalian hosts is crucial for the development of novel vaccines. Long noncoding RNAs (lncRNAs) derived from Plasmodium parasites or host cells have emerged as potential signaling molecules involved in the trafficking of proteins, RNA (mRNAs, miRNAs, and ncRNAs), and DNA. These lncRNAs facilitate the interaction between hosts and parasites, impacting normal physiology or pathology in malaria-infected individuals. Moreover, they possess the capacity to regulate immune responses and associated signaling pathways, thus potentially influencing chromatin organization, epigenetic modifications, mRNA processing, splicing, and translation. However, the functional role of exosomal lncRNAs in malaria remains poorly understood. This review offers a comprehensive analysis of lncRNA and exosomal lncRNA profiles during malaria infection. It presents an overview of recent progress in elucidating the involvement of exosomal lncRNAs in host-parasite interactions. Additionally, potential exosomal lncRNAs linked to the domains of innate and adaptive immunity in the context of malaria are proposed. These findings may contribute to the discovery of new diagnostic and therapeutic strategies for malaria. Furthermore, the need for additional research was highlighted that aimed to elucidate the mechanisms underlying lncRNA transportation into host cells and their targeting of specific genes to regulate the host's immune response. This knowledge gap presents an opportunity for future investigations, offering innovative approaches to enhance malarial control. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Small Molecule-RNA Interactions RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA in Disease and Development > RNA in Disease.

Decoding the impact of nuclear organization on antigenic variation in parasites

Antigenic variation as a strategy to evade the host adaptive immune response has evolved in divergent pathogens. Antigenic variation involves restricted, and often mutually exclusive, expression of dominant antigens and a periodic switch in antigen expression during infection. In eukaryotes, nuclear compartmentalization, including three-dimensional folding of the genome and physical separation of proteins in compartments or condensates, regulates mutually exclusive gene expression and chromosomal translocations. In this Review, we discuss the impact of nuclear organization on antigenic variation in the protozoan pathogens Trypanosoma brucei and Plasmodium falciparum. In particular, we highlight the relevance of nuclear organization in both mutually exclusive antigen expression and genome stability, which underlie antigenic variation.

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.

Discovery of RUF6 ncRNA-interacting proteins involved in P. falciparum immune evasion

Non-coding RNAs (ncRNAs) are emerging regulators of immune evasion and transmission of Plasmodium falciparum RUF6 is an ncRNA gene family that is transcribed by RNA polymerase III but actively regulates the Pol II-transcribed var virulence gene family. Understanding how RUF6 ncRNA connects to downstream effectors is lacking. We developed an RNA-directed proteomic discovery (ChIRP-MS) protocol to identify in vivo RUF6 ncRNA-protein interactions. The RUF6 ncRNA interactome was purified with biotinylated antisense oligonucleotides. Quantitative label-free mass spectrometry identified several unique proteins linked to gene transcription including RNA Pol II subunits, nucleosome assembly proteins, and a homologue of DEAD box helicase 5 (DDX5). Affinity purification of Pf-DDX5 identified proteins originally found by our RUF6-ChIRP protocol, validating the technique's robustness for identifying ncRNA interactomes in P. falciparum Inducible displacement of nuclear Pf-DDX5 resulted in significant down-regulation of the active var gene. Our work identifies a RUF6 ncRNA-protein complex that interacts with RNA Pol II to sustain the var gene expression, including a helicase that may resolve G-quadruplex secondary structures in var genes to facilitate transcriptional activation and progression.

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

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

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.

Functional Characterization of the Thrombospondin-Related Paralogous Proteins Rhoptry Discharge Factors 1 and 2 Unveils Phenotypic Plasticity in Toxoplasma gondii Rhoptry Exocytosis

To gain access to the intracellular cytoplasmic niche essential for their growth and replication, apicomplexan parasites such as Toxoplasma gondii rely on the timely secretion of two types of apical organelles named micronemes and rhoptries. Rhoptry proteins are key to host cell invasion and remodeling, however, the molecular mechanisms underlying the tight control of rhoptry discharge are poorly understood. Here, we report the identification and functional characterization of two novel T. gondii thrombospondin-related proteins implicated in rhoptry exocytosis. The two proteins, already annotated as MIC15 and MIC14, were renamed rhoptry discharge factor 1 (RDF1) and rhoptry discharge factor 2 (RDF2) and found to be exclusive of the Coccidia class of apicomplexan parasites. Furthermore, they were shown to have a paralogous relationship and share a C-terminal transmembrane domain followed by a short cytoplasmic tail. Immunofluorescence analysis of T. gondii tachyzoites revealed that RDF1 presents a diffuse punctate localization not reminiscent of any know subcellular compartment, whereas RDF2 was not detected. Using a conditional knockdown approach, we demonstrated that RDF1 loss caused a marked growth defect. The lack of the protein did not affect parasite gliding motility, host cell attachment, replication and egress, whereas invasion was dramatically reduced. Notably, while RDF1 depletion did not result in altered microneme exocytosis, rhoptry discharge was found to be heavily impaired. Interestingly, rhoptry secretion was reversed by spontaneous upregulation of the RDF2 gene in knockdown parasites grown under constant RDF1 repression. Collectively, our results identify RDF1 and RDF2 as additional key players in the pathway controlling rhoptry discharge. Furthermore, this study unveils a new example of compensatory mechanism contributing to phenotypic plasticity in T. gondii.

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.

Precise gene models using long-read sequencing reveal a unique poly(A) signal in Giardia lamblia

During pre-mRNA processing, the poly(A) signal is recognized by a protein complex that ensures precise cleavage and polyadenylation of the nascent transcript. The location of this cleavage event establishes the length and sequence of the 3' UTR of an mRNA, thus determining much of its post-transcriptional fate. Using long-read sequencing, we characterize the polyadenylation signal and related sequences surrounding Giardia lamblia cleavage sites for over 2600 genes. We find that G. lamblia uses an AGURAA poly(A) signal, which differs from the mammalian AAUAAA. We also describe how G. lamblia lacks common auxiliary elements found in other eukaryotes, along with the proteins that recognize them. Further, we identify 133 genes with evidence of alternative polyadenylation. These results suggest that despite pared-down cleavage and polyadenylation machinery, 3' end formation still appears to be an important regulatory step for gene expression in G. lamblia.

RNA Sequencing Reveals Widespread Transcription of Natural Antisense RNAs in Entamoeba Species

Entamoeba is a genus of Amoebozoa that includes the intestine-colonizing pathogenic species Entamoeba histolytica. To understand the basis of gene regulation in E. histolytica from an evolutionary perspective, we have profiled the transcriptomes of its closely related species E. dispar, E. moshkovskii and E. invadens. Genome-wide identification of transcription start sites (TSS) and polyadenylation sites (PAS) revealed the similarities and differences of their gene regulatory sequences. In particular, we found the widespread initiation of antisense transcription from within the gene coding sequences is a common feature among all Entamoeba species. Interestingly, we observed the enrichment of antisense transcription in genes involved in several processes that are common to species infecting the human intestine, e.g., the metabolism of phospholipids. These results suggest a potentially conserved and compact gene regulatory system in Entamoeba.

Resetting var Gene Transcription in Plasmodium falciparum

One of the key mechanisms contributing to the virulence of Plasmodium falciparum is its ability to undergo antigenic switching among antigenically distinct variants of the PfEMP1 adhesive proteins, encoded by the var gene family. To avoid premature exposure of its antigenic repertoire, the parasite transcribes its var genes in a mutually exclusive manner, and switch expression at a very slow rate. This process is epigenetically regulated and it relies on "epigenetic memory," which imprints the single active var gene to remain active for multiple replication cycles. Erasing this epigenetic memory in parasites grown in culture resembles parasites, which egress from the liver. It could therefore be of interest for investigating var switching patterns at the onset of malaria infections. In addition, this procedure could be used for creating heterogeneity of var expression among parasite populations. The methodology described here for resetting of var gene expression is based on promoter titration, also known as molecular sponging.

A cascade of transcriptional repression determines sexual commitment and development in Plasmodium falciparum

Gametocytogenesis, the process by which malaria parasites produce sexual forms that can infect mosquitoes, is essential for the transmission of malaria. A transcriptional switch of the pfap2-g gene triggers sexual commitment, but how the complex multi-step process is precisely programed remains largely unknown. Here, by systematic functional screening of a panel of ApiAP2 transcription factors, we identify six new ApiAP2 members associated with gametocytogenesis in Plasmodium falciparum. Among these, PfAP2-G5 (PF3D7_1139300) was found to be indispensable for gametocytogenesis. This factor suppresses the transcriptional activity of the pfap2-g gene via binding to both the upstream region and exonic gene body, the latter is linked to the maintenance of local heterochromatin structure, thereby preventing initiation of sexual commitment. Removal of this repressive effect through pfap2-g5 knockout disrupts the asexual replication cycle and promotes sexual commitment accompanied by upregulation of pfap2-g expression. However, the gametocytes produced fail to mature fully. Further analyses show that PfAP2-G5 is essential for gametocyte maturation, and causes the down-regulation of pfap2-g and a set of early gametocyte genes activated by PfAP2-G prior to gametocyte development. Collectively, our findings reveal a regulation cascade of gametocyte production in malaria parasites, and provide a new target for transmission blocking interventions.

Non-coding Natural Antisense Transcripts: Analysis and Application

Non-coding natural antisense transcripts (ncNATs) are regulatory RNA sequences that are transcribed in the opposite direction to protein-coding or non-coding transcripts. These transcripts are implicated in a broad variety of biological and pathological processes, including tumorigenesis and oncogenic progression. With this complex field still in its infancy, annotations, expression profiling and functional characterisations of ncNATs are far less comprehensive than those for protein-coding genes, pointing out substantial gaps in the analysis and characterisation of these regulatory transcripts. In this review, we discuss ncNATs from an analysis perspective, in particular regarding the use of high-throughput sequencing strategies, such as RNA-sequencing, and summarize the unique challenges of investigating the antisense transcriptome. Finally, we elaborate on their potential as biomarkers and future targets for treatment, focusing on cancer.

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.

Identification of long noncoding natural antisense transcripts (lncNATs) correlated with drought stress response in wild rice (Oryza nivara)

Background

Wild rice, including Oryza nivara and Oryza rufipogon, which are considered as the ancestors of Asian cultivated rice (Oryza sativa), possess high genetic diversity and serve as a crucial resource for breeding novel cultivars of cultivated rice. Although rice domestication related traits, such as seed shattering and plant architecture, have been intensively studied at the phenotypic and genomic levels, further investigation is needed to understand the molecular basis of phenotypic differences between cultivated and wild rice. Drought stress is one of the most severe abiotic stresses affecting rice growth and production. Adaptation to drought stress involves a cascade of genes and regulatory factors that form complex networks. O. nivara inhabits swampy areas with a seasonally dry climate, which is an ideal material to discover drought tolerance alleles. Long noncoding natural antisense transcripts (lncNATs), a class of long noncoding RNAs (lncRNAs), regulate the corresponding sense transcripts and play an important role in plant growth and development. However, the contribution of lncNATs to drought stress response in wild rice remains largely unknown.

Results

Here, we conducted strand-specific RNA sequencing (ssRNA-seq) analysis of Nipponbare (O. sativa) and two O. nivara accessions (BJ89 and BJ278) to determine the role of lncNATs in drought stress response in wild rice. A total of 1246 lncRNAs were identified, including 1091 coding–noncoding NAT pairs, of which 50 were expressed only in Nipponbare, and 77 were expressed only in BJ89 and/or BJ278. Of the 1091 coding–noncoding NAT pairs, 240 were differentially expressed between control and drought stress conditions. Among these 240 NAT pairs, 12 were detected only in Nipponbare, and 187 were detected uniquely in O. nivara. Furthermore, 10 of the 240 coding–noncoding NAT pairs were correlated with genes enriched in stress responsive GO terms; among these, nine pairs were uniquely found in O. nivara, and one pair was shared between O. nivara and Nipponbare.

Conclusion

We identified lncNATs associated with drought stress response in cultivated rice and O. nivara. These results will improve our understanding of the function of lncNATs in drought tolerance and accelerate rice breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07754-4.

Gene regulation of the avian malaria parasite Plasmodium relictum, during the different stages within the mosquito vector

The malaria parasite Plasmodium relictum is one of the most widespread species of avian malaria. As in the case of its human counterparts, bird Plasmodium undergoes a complex life cycle infecting two hosts: the arthropod vector and the vertebrate host. In this study, we examined transcriptomes of P. relictum (SGS1) during crucial timepoints within its vector, Culex pipiens quinquefasciatus. Differential gene-expression analyses identified genes linked to the parasites life-stages at: i) a few minutes after the blood meal is ingested, ii) during peak oocyst production phase, iii) during peak sporozoite phase and iv) during the late-stages of the infection. A large amount of genes coding for functions linked to host-immune invasion and multifunctional genes was active throughout the infection cycle. One gene associated with a conserved Plasmodium membrane protein with unknown function was upregulated throughout the parasite development in the vector, suggesting an important role in the successful completion of the sporogonic cycle. Gene expression analysis further identified genes, with unknown functions to be significantly differentially expressed during the infection in the vector as well as upregulation of reticulocyte-binding proteins, which raises the possibility of the multifunctionality of these RBPs. We establish the existence of highly stage-specific pathways being overexpressed during the infection. This first study of gene-expression of a non-human Plasmodium species in its vector provides a comprehensive insight into the molecular mechanisms of the common avian malaria parasite P. relictum and provides essential information on the evolutionary diversity in gene regulation of the Plasmodium's vector stages.

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.

Direct Nanopore Sequencing of mRNA Reveals Landscape of Transcript Isoforms in Apicomplexan Parasites

Alternative splicing is a widespread phenomenon in metazoans by which single genes are able to produce multiple isoforms of the gene product. However, this has been poorly characterized in apicomplexans, a major phylum of some of the most important global parasites. Efforts have been hampered by atypical transcriptomic features, such as the high AU content of Plasmodium RNA, but also the limitations of short-read sequencing in deciphering complex splicing events. In this study, we utilized the long read direct RNA sequencing platform developed by Oxford Nanopore Technologies to survey the alternative splicing landscape of Toxoplasma gondii and Plasmodium falciparum. We find that while native RNA sequencing has a reduced throughput, it allows us to obtain full-length or nearly full-length transcripts with comparable quantification to Illumina sequencing. By comparing these data with available gene models, we find widespread alternative splicing, particularly intron retention, in these parasites. Most of these transcripts contain premature stop codons, suggesting that in these parasites, alternative splicing represents a pathway to transcriptomic diversity, rather than expanding proteomic diversity. Moreover, alternative splicing rates are comparable between parasites, suggesting a shared splicing machinery, despite notable transcriptomic differences between the parasites. This study highlights a strategy in using long-read sequencing to understand splicing events at the whole-transcript level and has implications in the future interpretation of transcriptome sequencing studies.

IMPORTANCE We have used a novel nanopore sequencing technology to directly analyze parasite transcriptomes. The very long reads of this technology reveal the full-length genes of the parasites that cause malaria and toxoplasmosis. Gene transcripts must be processed in a process called splicing before they can be translated to protein. Our analysis reveals that these parasites very frequently only partially process their gene products, in a manner that departs dramatically from their human hosts.

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.

Analysis of Long Non-Coding RNA in Cryptosporidium parvum Reveals Significant Stage-Specific Antisense Transcription

Cryptosporidium is a protist parasite that has been identified as the second leading cause of moderate to severe diarrhea in children younger than two and a significant cause of mortality worldwide. Cryptosporidium has a complex, obligate, intracellular but extra cytoplasmic lifecycle in a single host. How genes are regulated in this parasite remains largely unknown. Long non-coding RNAs (lncRNAs) play critical regulatory roles, including gene expression across a broad range of organisms. Cryptosporidium lncRNAs have been reported to enter the host cell nucleus and affect the host response. However, no systematic study of lncRNAs in Cryptosporidium has been conducted to identify additional lncRNAs. In this study, we analyzed a C. parvum in vitro strand-specific RNA-seq developmental time series covering both asexual and sexual stages to identify lncRNAs associated with parasite development. In total, we identified 396 novel lncRNAs, mostly antisense, with 86% being differentially expressed. Surprisingly, nearly 10% of annotated mRNAs have an antisense transcript. lncRNAs occur most often at the 3' end of their corresponding sense mRNA. Putative lncRNA regulatory regions were identified and many appear to encode bidirectional promoters. A positive correlation between lncRNA and upstream mRNA expression was observed. Evolutionary conservation and expression of lncRNA candidates was observed between C. parvum, C. hominis and C. baileyi. Ten C. parvum protein-encoding genes with antisense transcripts have P. falciparum orthologs that also have antisense transcripts. Three C. parvum lncRNAs with exceptional properties (e.g., intron splicing) were experimentally validated using RT-PCR and RT-qPCR. This initial characterization of the C. parvum non-coding transcriptome facilitates further investigations into the roles of lncRNAs in parasite development and host-pathogen interactions.

Genomic determinants for initiation and length of natural antisense transcripts in Entamoeba histolytica

Natural antisense transcripts (NAT) have been reported in prokaryotes and eukaryotes. While the functions of most reported NATs remain unknown, their potentials in regulating the transcription of their counterparts have been speculated. Entamoeba histolytica, which is a unicellular eukaryotic parasite, has a compact protein-coding genome with very short intronic and intergenic regions. The regulatory mechanisms of gene expression in this compact genome are under-described. In this study, by genome-wide mapping of RNA-Seq data in the genome of E. histolytica, we show that a substantial fraction of its protein-coding genes (28%) has significant transcription on their opposite strand (i.e. NAT). Intriguingly, we found the location of transcription start sites or polyadenylation sites of NAT are determined by the specific motifs encoded on the opposite strand of the gene coding sequences, thereby providing a compact regulatory system for gene transcription. Moreover, we demonstrated that NATs are globally up-regulated under various environmental conditions including temperature stress and pathogenicity. While NATs do not appear to be consequences of spurious transcription, they may play a role in regulating gene expression in E. histolytica, a hypothesis which needs to be tested.

Analysis of Plasmodium vivax schizont transcriptomes from field isolates reveals heterogeneity of expression of genes involved in host-parasite interactions

Plasmodium vivax gene regulation remains difficult to study due to the lack of a robust in vitro culture method, low parasite densities in peripheral circulation and asynchronous parasite development. We adapted an RNA-seq protocol “DAFT-seq” to sequence the transcriptome of four P. vivax field isolates that were cultured for a short period ex vivo before using a density gradient for schizont enrichment. Transcription was detected from 78% of the PvP01 reference genome, despite being schizont-enriched samples. This extensive data was used to define thousands of 5′ and 3′ untranslated regions, some of which overlapped with neighbouring transcripts, and to improve the gene models of 352 genes, including identifying 20 novel gene transcripts. This dataset has also significantly increased the known amount of heterogeneity between P. vivax schizont transcriptomes from individual patients. The majority of genes found to be differentially expressed between the isolates lack Plasmodium falciparum homologs and are predicted to be involved in host-parasite interactions, with an enrichment in reticulocyte binding proteins, merozoite surface proteins and exported proteins with unknown function. An improved understanding of the diversity within P. vivax transcriptomes will be essential for the prioritisation of novel vaccine targets.

Characterization of natural antisense transcripts arisen from the locus encoding Toxoplasma gondii ubiquitin-like protease

Natural antisense transcripts (NATs) are non-protein coding RNAs that could play an important role in regulating the expression of their counterpart protein encoding sense transcript. Although NATs are widespread in most eukaryotic genomes, very little is known about their functions. This study focuses on gaining a better understanding of the function of NATs in Toxoplasma gondii, a pathogenic unicellular eukaryote. Previously, we characterized the gene encoding the first committed enzyme in sumoylation, named ubiquitin-like protease 1 (TgUlp1), and showed that the expression of TgUlp1 is vital to the life cycle of T. gondii. Interestingly, the locus of TgUlp1 also transcribes a NAT species. Using a dual luciferase assay, we identified the promoter of TgUlp1 NAT to be located within the 3'-region of its counterpart coding sequence. While TgUlp1 mRNA level was detected at a lower level throughout the life cycle of T. gondii, its NAT level was upregulated when the parasite converts from actively replicating tachyzoite form to slowly growing bradyzoite form. To investigate the effect of TgUlp1 NAT on the expression of its counterpart mRNA, we used a reporter system bearing TgUlp1 mRNA sequences and showed that the single-stranded TgUlp1 NAT and its in vitro RNase III processed products have the ability to lower the expression of the reporter system. Using a transgenic Dicer-knockout (TgDicer-KO) strain, we showed that TgDicer is required for the function of TgUlp1 NAT in vivo. The findings strongly suggest that the RNA interference pathway is necessary for the function of TgUlp1 NAT.

Genome wide distribution of G-quadruplexes and their impact on gene expression in malaria parasites

Mechanisms of transcriptional control in malaria parasites are still not fully understood. The positioning patterns of G-quadruplex (G4) DNA motifs in the parasite's AT-rich genome, especially within the var gene family which encodes virulence factors, and in the vicinity of recombination hotspots, points towards a possible regulatory role of G4 in gene expression and genome stability. Here, we carried out the most comprehensive genome-wide survey, to date, of G4s in the Plasmodium falciparum genome using G4Hunter, which identifies G4 forming sequences (G4FS) considering their G-richness and G-skewness. We show an enrichment of G4FS in nucleosome-depleted regions and in the first exon of var genes, a pattern that is conserved within the closely related Laverania Plasmodium parasites. Under G4-stabilizing conditions, i.e., following treatment with pyridostatin (a high affinity G4 ligand), we show that a bona fide G4 found in the non-coding strand of var promoters modulates reporter gene expression. Furthermore, transcriptional profiling of pyridostatin-treated parasites, shows large scale perturbations, with deregulation affecting for instance the ApiAP2 family of transcription factors and genes involved in ribosome biogenesis. Overall, our study highlights G4s as important DNA secondary structures with a role in Plasmodium gene expression regulation, sub-telomeric recombination and var gene biology.

Actin-related protein Arp4 regulates euchromatic gene expression and development through H2A.Z deposition in blood-stage Plasmodium falciparum

BACKGROUND

Malaria caused by Plasmodium spp. is still a major threat to public health globally. The various approaches to developing new antimalarial agents rely on the understanding of the complex regulatory mechanisms of dynamic gene expression in the life-cycle of these malaria parasites. The nuclear members of the evolutionarily conserved actin-related protein nuclear (ARP) superfamily are the major components of nucleosome remodelling complexes. In the human malaria parasite Plasmodium falciparum, bioinformatics analysis has predicted three ARP orthologues: PfArp1, PfArp4 and PfArp6. However, little is known about the biological functions of putative PfArp4. In this study, we aimed to investigate the function and the underlying mechanisms of PfArp4 gene regulation.

METHODS

A conditional gene knockdown approach was adopted by incorporating the glucosamine-inducible glmS ribozyme sequence into the 3' UTR of the PfArp4 and PfArp6 genes. The transgenic parasites PfArp4-Ty1-Ribo, PfArp6-Ty1-Ribo and pL6-PfArp4-Ty1::PfArp6-HA were generated by the CRISPR-Cas9 technique. The knockdown effect in the transgenic parasite was measured by growth curve assay and western blot (WB) analysis. The direct interaction between PfArp4 and PfArp6 was validated by co-IFA and co-IP assays. The euchromatic gene expression mediated through H2A.Z (histone H2A variant) deposition and H3K9ac modification at promoters and regulated by PfArp4, was determined by RNA-seq and ChIP-seq.

RESULTS

The inducible knockdown of PfArp4 inhibited blood-stage development of P. falciparum. PfArp4 and PfArp6 were colocalized in the nucleus of P. falciparum parasites. PfArp4 gene knockdown altered the global transcriptome. PfArp4 protein colocalized with the histone variant H2A.Z and euchromatic marker H3K9ac in intergenic regions. The inducible downregulation of PfArp4 resulted in the depletion of H2A.Z and lower H3K9ac levels at the upstream regions of eukaryotic genes, thereby repressing the transcriptional abundance of H2A.Z-dependent genes.

CONCLUSIONS

Our findings suggest that PfArp4 regulates the cell cycle by controlling H2A.Z deposition and affecting centromere function, contributing to the understanding the complex epigenetic regulation of gene expression and the development of P. falciparum.

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.

Differential gene expression of Plasmodium homocircumflexum (lineage pCOLL4) across two experimentally infected passerine bird species

Plasmodium parasites are present in a wide range of host species, some of which tend to be more susceptible than others, potentially as an outcome of evolved tolerance or resistance. Common starlings seem to cope with malaria infection while common crossbills are more susceptible to the same infections. That raises the question if the parasites rely on the same molecular mechanisms regardless of host species or do Plasmodium parasites change gene-expressions in accordance to the environment different hosts might provide? We used RNA-sequencing from starlings and crossbills, experimentally infected with Plasmodium homocircumflexum (lineage pCOLL4). The assembled transcriptome contained a total of 26,733 contigs. Parasite expression patterns differed between bird species. Parasites had higher expression of cell-invasion genes when infecting crossbills compared to starlings whereas in starlings genes related to apoptosis or/and oxidative stress showed higher expression levels. This article reveals how a Plasmodium parasite might adjust its expression and gene function depending on the host species infected.

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.

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.

Transcriptome analysis of Plasmodium berghei during exo-erythrocytic development

Background

The complex life cycle of malaria parasites requires well-orchestrated stage specific gene expression. In the vertebrate host the parasites grow and multiply by schizogony in two different environments: within erythrocytes and within hepatocytes. Whereas erythrocytic parasites are well-studied in this respect, relatively little is known about the exo-erythrocytic stages.

Methods

In an attempt to fill this gap, genome wide RNA-seq analyses of various exo-erythrocytic stages of Plasmodium berghei including sporozoites, samples from a time-course of liver stage development and detached cells were performed. These latter contain infectious merozoites and represent the final step in exo-erythrocytic development.

Results

The analysis represents the complete transcriptome of the entire life cycle of P. berghei parasites with temporal detailed analysis of the liver stage allowing comparison of gene expression across the progression of the life cycle. These RNA-seq data from different developmental stages were used to cluster genes with similar expression profiles, in order to infer their functions. A comparison with published data from other parasite stages confirmed stage-specific gene expression and revealed numerous genes that are expressed differentially in blood and exo-erythrocytic stages. One of the most exo-erythrocytic stage-specific genes was PBANKA_1003900, which has previously been annotated as a “gametocyte specific protein”. The promoter of this gene drove high GFP expression in exo-erythrocytic stages, confirming its expression profile seen by RNA-seq.

Conclusions

The comparative analysis of the genome wide mRNA expression profiles of erythrocytic and different exo-erythrocytic stages could be used to improve the understanding of gene regulation in Plasmodium parasites and can be used to model exo-erythrocytic stage metabolic networks toward the identification of differences in metabolic processes during schizogony in erythrocytes and hepatocytes.

Dissecting the Gene Expression, Localization, Membrane Topology, and Function of the Plasmodium falciparum STEVOR Protein Family

During its intraerythrocytic development, the malaria parasite Plasmodium falciparum exposes variant surface antigens (VSAs) on infected erythrocytes to establish and maintain an infection. One family of small VSAs is the polymorphic STEVOR proteins, which are marked for export to the host cell surface through their PEXEL signal peptide. Interestingly, some STEVORs have also been reported to localize to the parasite plasma membrane and apical organelles, pointing toward a putative function in host cell egress or invasion. Using deep RNA sequencing analysis, we characterized P. falciparumstevor gene expression across the intraerythrocytic development cycle, including free merozoites, in detail and used the resulting stevor expression profiles for hierarchical clustering. We found that most stevor genes show biphasic expression oscillation, with maximum expression during trophozoite stages and a second peak in late schizonts. We selected four STEVOR variants, confirmed the expected export of these proteins to the host cell membrane, and tracked them to a secondary location, either to the parasite plasma membrane or the secretory organelles of merozoites in late schizont stages. We investigated the function of a particular STEVOR that showed rhoptry localization and demonstrated its role at the parasite-host interface during host cell invasion by specific antisera and targeted gene disruption. Experimentally determined membrane topology of this STEVOR revealed a single transmembrane domain exposing the semiconserved as well as variable protein regions to the cell surface.IMPORTANCE Malaria claims about half a million lives each year. Plasmodium falciparum, the causative agent of the most severe form of the disease, uses proteins that are translocated to the surface of infected erythrocytes for immune evasion. To circumvent the detection of these gene products by the immune system, the parasite evolved a complex strategy that includes gene duplications and elaborate sequence polymorphism. STEVORs are one family of these variant surface antigens and are encoded by about 40 genes. Using deep RNA sequencing of blood-stage parasites, including free merozoites, we first established stevor expression of the cultured isolate and compared it with published transcriptomes. We reveal a biphasic expression of most stevor genes and confirm this for individual STEVORs at the protein level. The membrane topology of a rhoptry-associated variant was experimentally elucidated and linked to host cell invasion, underlining the importance of this multifunctional protein family for parasite proliferation.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Disruption of the RNA exosome reveals the hidden face of the malaria parasite transcriptome

Antisense transcription emerges as a key regulator of important biological processes in the human malaria parasite Plasmodium falciparum. RNA-processing factors, however, remain poorly characterized in this pathogen. Here, we purified the multiprotein RNA exosome complex of malaria parasites by affinity chromatography, using HA-tagged PfRrp4 and PfDis3 as the ligands. Seven distinct core exosome subunits (PfRrp41, PfMtr3, PfRrp42, PfRrp45, PfRrp4, PfRrp40, PfCsl4) and two exoribonuclease proteins PfRrp6 and PfDis3 are identified by mass spectrometry. Western blot analysis detects Dis3 and Rrp4 predominantly in the cytoplasmic fraction during asexual blood stage development. An inducible gene knock out of the PfDis3 subunit reveals the upregulation of structural and coding RNA, but the vast majority belongs to antisense RNA. Furthermore, we detect numerous types of cryptic unstable transcripts (CUTs) linked to virulence gene families including antisense RNA in the rif gene family. Our work highlights the limitations of steady-state RNA analysis to predict transcriptional activity and link the RNA surveillance machinery directly with post-transcriptional control and gene expression in malaria parasites.

Characterization of mRNA polyadenylation in the apicomplexa

Messenger RNA polyadenylation is a universal aspect of gene expression in eukaryotes. In well-established model organisms, this process is mediated by a conserved complex of 15–20 subunits. To better understand this process in apicomplexans, a group of unicellular parasites that causes serious disease in humans and livestock, a computational and high throughput sequencing study of the polyadenylation complex and poly(A) sites in several species was conducted. BLAST-based searches for orthologs of the human polyadenylation complex yielded clear matches to only two—poly(A) polymerase and CPSF73—of the 19 proteins used as queries in this analysis. As the human subunits that recognize the AAUAAA polyadenylation signal (PAS) were not immediately obvious, a computational analysis of sequences adjacent to experimentally-determined apicomplexan poly(A) sites was conducted. The results of this study showed that there exists in apicomplexans an A-rich region that corresponds in position to the AAUAAA PAS. The set of experimentally-determined sites in one species, Sarcocystis neurona, was further analyzed to evaluate the extent and significance of alternative poly(A) site choice in this organism. The results showed that almost 80% of S. neurona genes possess more than one poly(A) site, and that more than 780 sites showed differential usage in the two developmental stages–extracellular merozoites and intracellular schizonts–studied. These sites affected more than 450 genes, and included a disproportionate number of genes that encode membrane transporters and ribosomal proteins. Taken together, these results reveal that apicomplexan species seem to possess a poly(A) signal analogous to AAUAAA even though genes that may encode obvious counterparts of the AAUAAA-recognizing proteins are absent in these organisms. They also indicate that, as is the case in other eukaryotes, alternative polyadenylation is a widespread phenomenon in S. neurona that has the potential to impact growth and development.

Widespread antisense transcription of Populus genome under drought

Antisense transcription is widespread in many genomes and plays important regulatory roles in gene expression. The objective of our study was to investigate the extent and functional relevance of antisense transcription in forest trees. We employed Populus, a model tree species, to probe the antisense transcriptional response of tree genome under drought, through stranded RNA-seq analysis. We detected nearly 48% of annotated Populus gene loci with antisense transcripts and 44% of them with co-transcription from both DNA strands. Global distribution of reads pattern across annotated gene regions uncovered that antisense transcription was enriched in untranslated regions while sense reads were predominantly mapped in coding exons. We further detected 1185 drought-responsive sense and antisense gene loci and identified a strong positive correlation between the expression of antisense and sense transcripts. Additionally, we assessed the antisense expression in introns and found a strong correlation between intronic expression and exonic expression, confirming antisense transcription of introns contributes to transcriptional activity of Populus genome under drought. Finally, we functionally characterized drought-responsive sense-antisense transcript pairs through gene ontology analysis and discovered that functional groups including transcription factors and histones were concordantly regulated at both sense and antisense transcriptional level. Overall, our study demonstrated the extensive occurrence of antisense transcripts of Populus genes under drought and provided insights into genome structure, regulation pattern and functional significance of drought-responsive antisense genes in forest trees. Datasets generated in this study serve as a foundation for future genetic analysis to improve our understanding of gene regulation by antisense transcription.

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

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

Plasmodium vivax Biology: Insights Provided by Genomics, Transcriptomics and Proteomics

During the last decade, the vast omics field has revolutionized biological research, especially the genomics, transcriptomics and proteomics branches, as technological tools become available to the field researcher and allow difficult question-driven studies to be addressed. Parasitology has greatly benefited from next generation sequencing (NGS) projects, which have resulted in a broadened comprehension of basic parasite molecular biology, ecology and epidemiology. Malariology is one example where application of this technology has greatly contributed to a better understanding of Plasmodium spp. biology and host-parasite interactions. Among the several parasite species that cause human malaria, the neglected Plasmodium vivax presents great research challenges, as in vitro culturing is not yet feasible and functional assays are heavily limited. Therefore, there are gaps in our P. vivax biology knowledge that affect decisions for control policies aiming to eradicate vivax malaria in the near future. In this review, we provide a snapshot of key discoveries already achieved in P. vivax sequencing projects, focusing on developments, hurdles, and limitations currently faced by the research community, as well as perspectives on future vivax malaria research.

Identification and characterization of L1-specific endo-siRNAs essential for early embryonic development in pig

Small noncoding RNAs (sncRNAs) play important roles in RNA interference (RNAi). In addition to microRNA (miRNA) and Piwi-interacting RNA (piRNA), one key member of sncRNAs group is endogenous small interfering RNA (endo-siRNA). Some studies do show the role of endo-siRNAs in Dicer and/or Ago mutants, however, the biological functions of specific endo-siRNAs remains mostly unanswered. In the study, we have performed a comparative analysis of endo-siRNAs present in porcine sperms, oocytes and zygotes, identified by deep sequencing and bioinformatics analysis. Further, we observe a large amount of endo-siRNAs specific binding on ORF2 and 3′ UTR of porcine L1 (L1-siRNAs). And, 9 L1-siRNAs generated from a dsRNA formed between L1 transcript and a newly identified an antisense noncoding RNA was characterized. We show the L1-siRNAs regulate early embryonic development by inhibiting the activity of L1 retrotransposition. This work can contribute to understanding the functional role of abundant endo-siRNAs in embryonic development.

An ApiAP2 member regulates expression of clonally variant genes of the human malaria parasite Plasmodium falciparum

Variegated surface antigen expression is key to chronic infection and pathogenesis of the human malaria parasite Plasmodium falciparum. This protozoan parasite expresses distinct surface molecules that are encoded by clonally variant gene families such as var, rif and stevor. The molecular mechanisms governing activation of individual members remain ill-defined. To investigate the molecular events of the initial transcriptional activation process we focused on a member of the apicomplexan ApiAP2 transcription factor family predicted to bind to the 5′ upstream regions of the var gene family, AP2-exp (PF3D7_1466400). Viable AP2-exp mutant parasites rely on expressing no less than a short truncated protein including the N-terminal AP2 DNA-binding domain. RNA-seq analysis in mutant parasites revealed transcriptional changes in a subset of exported proteins encoded by clonally variant gene families. Upregulation of RIFINs and STEVORs was validated at the protein levels. In addition, morphological alterations were observed on the surface of the host cells infected by the mutants. This work points to a complex regulatory network of clonally variant gene families in which transcription of a subset of members is regulated by the same transcription factor. In addition, we highlight the importance of the non-DNA binding AP2 domain in functional gene regulation.

Genome-wide identification of novel intergenic enhancer-like elements: implications in the regulation of transcription in Plasmodium falciparum

Background

The molecular mechanisms of transcriptional regulation are poorly understood in Plasmodium falciparum. In addition, most of the genes in Plasmodium falciparum are transcriptionally poised and only a handful of cis-regulatory elements are known to operate in transcriptional regulation. Here, we employed an epigenetic signature based approach to identify significance of previously uncharacterised intergenic regions enriched with histone modification marks leading to discovery of enhancer-like elements.

Results

We found that enhancer-like elements are significantly enriched with H3K4me1, generate unique non-coding bi-directional RNAs and majority of them can function as cis-regulators. Furthermore, functional enhancer reporter assay demonstrates that the enhancer-like elements regulate transcription of target genes in Plasmodium falciparum. Our study also suggests that the Plasmodium genome segregates functionally related genes into discrete housekeeping and pathogenicity/virulence clusters, presumably for robust transcriptional control of virulence/pathogenicity genes.

Conclusions

This report contributes to the understanding of parasite regulatory genomics by identification of enhancer-like elements, defining their epigenetic and transcriptional features and provides a resource of functional cis-regulatory elements that may give insights into the virulence/pathogenicity of Plasmodium falciparum.

Electronic supplementary material

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

Characterization of P. vivax blood stage transcriptomes from field isolates reveals similarities among infections and complex gene isoforms

Our understanding of the structure and regulation of Plasmodium vivax genes is limited by our inability to grow the parasites in long-term in vitro cultures. Most P. vivax studies must therefore rely on patient samples, which typically display a low proportion of parasites and asynchronous parasites. Here, we present stranded RNA-seq data generated directly from a small volume of blood from three Cambodian vivax malaria patients collected before treatment. Our analyses show surprising similarities of the parasite gene expression patterns across infections, despite extensive variations in parasite stage proportion. These similarities contrast with the unique gene expression patterns observed in sporozoites isolated from salivary glands of infected Colombian mosquitoes. Our analyses also indicate that more than 10% of P. vivax genes encode multiple, often undescribed, protein-coding sequences, potentially increasing the diversity of proteins synthesized by blood stage parasites. These data also greatly improve the annotations of P. vivax gene untranslated regions, providing an important resource for future studies of specific genes.

A long and abundant non-coding RNA in Lactobacillus salivarius

Lactobacillus salivarius, found in the intestinal microbiota of humans and animals, is studied as an example of the sub-dominant intestinal commensals that may impart benefits upon their host. Strains typically harbour at least one megaplasmid that encodes functions contributing to contingency metabolism and environmental adaptation. RNA sequencing (RNA-seq)transcriptomic analysis of L. salivarius strain UCC118 identified the presence of a novel unusually abundant long non-coding RNA (lncRNA) encoded by the megaplasmid, and which represented more than 75 % of the total RNA-seq reads after depletion of rRNA species. The expression level of this 520 nt lncRNA in L. salivarius UCC118 exceeded that of the 16S rRNA, it accumulated during growth, was very stable over time and was also expressed during intestinal transit in a mouse. This lncRNA sequence is specific to the L. salivarius species; however, among 45 L. salivarius genomes analysed, not all (only 34) harboured the sequence for the lncRNA. This lncRNA was produced in 27 tested L. salivarius strains, but at strain-specific expression levels. High-level lncRNA expression correlated with high megaplasmid copy number. Transcriptome analysis of a deletion mutant lacking this lncRNA identified altered expression levels of genes in a number of pathways, but a definitive function of this new lncRNA was not identified. This lncRNA presents distinctive and unique properties, and suggests potential basic and applied scientific developments of this phenomenon.