Analysis of small nucleolar RNAs reveals unique genetic features in malaria parasites

Starvation induces changes in abundance and small RNA cargo of extracellular vesicles released from Plasmodium falciparum infected red blood cells

The lethal malaria parasite Plasmodium falciparum needs to constantly respond and adapt to changes within the human host in order to survive and transmit. One such change is composed of nutritional limitation, which is augmented with increased parasite loads and intimately linked to severe disease development. Extracellular vesicles released from infected red blood cells have been proposed as important mediators of disease pathogenesis and intercellular communication but whether important for the parasite response to nutritional availability is unknown. Therefore, we investigated the abundance and small RNA cargo of extracellular vesicles released upon short-term nutritional starvation of P. falciparum in vitro cultures. We show that primarily ring-stage parasite cultures respond to glucose and amino acid deprivation with an increased release of extracellular vesicles. Small RNA sequencing of these extracellular vesicles further revealed human miRNAs and parasitic tRNA fragments as the main constituent biotypes. Short-term starvations led to alterations in the transcriptomic profile, most notably in terms of the over-represented biotypes. These data suggest a potential role for extracellular vesicles released from P. falciparum infected red blood cells in the response to nutritional perturbations, their potential as prognostic biomarkers and point towards an evolutionary conserved role among protozoan parasites.

Ribosome Specialization in Protozoa Parasites

Ribosomes, in general, are viewed as constitutive macromolecular machines where protein synthesis takes place; however, this view has been recently challenged, supporting the hypothesis of ribosome specialization and opening a completely new field of research. Recent studies have demonstrated that ribosomes are heterogenous in their nature and can provide another layer of gene expression control by regulating translation. Heterogeneities in ribosomal RNA and ribosomal proteins that compose them favor the selective translation of different sub-pools of mRNAs and functional specialization. In recent years, the heterogeneity and specialization of ribosomes have been widely reported in different eukaryotic study models; however, few reports on this topic have been made on protozoa and even less on protozoa parasites of medical importance. This review analyzes heterogeneities of ribosomes in protozoa parasites highlighting the specialization in their functions and their importance in parasitism, in the transition between stages in their life cycle, in the change of host and in response to environmental conditions.

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

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

Structural Analysis of the Plasmodial Proteins ZNHIT3 and NUFIP1 Provides Insights into the Selectivity of a Conserved Interaction

Malaria is a widespread and lethal disease caused by the Plasmodium parasites that can infect human beings through Anopheles mosquitoes. For that reason, the biology of Plasmodium needs to be studied to develop antimalarial treatments. By determining the three-dimensional structures of macromolecules, structural biology helps to understand the function of proteins and can reveal how interactions occur between biological partners. Here, we studied the ZNHIT3 and NUFIP1 proteins from Plasmodium falciparum, two proteins tightly linked to the ribosome biology. Due to their important functions in post-translational modifications of ribosomal RNAs and in ribophagy, these proteins participate in the survival of cells. In this study, we solved the three-dimensional structure of a thermally stable and species-dependent complex between fragments of these proteins. Our results were compared to the AlphaFold predictions, which motivated the study of the free ZNHIT3 fragment that binds NUFIP1. We showed that the latter fragment multimerized in vitro but also had the inner ability to change its conformation to escape the solvent exposition of key hydrophobic residues involved in the interaction with NUFIP1. Our data could open the gate to selective drug discovery processes involving these two proteins.

'2TM proteins': an antigenically diverse superfamily with variable functions and export pathways

Malaria is a disease that affects millions of people annually. An intracellular habitat and lack of protein synthesizing machinery in erythrocytes pose numerous difficulties for survival of the human pathogen Plasmodium falciparum. The parasite refurbishes the infected red blood cell (iRBC) by synthesis and export of several proteins in an attempt to suffice its metabolic needs and evade the host immune response. Immune evasion is largely mediated by surface display of highly polymorphic protein families known as variable surface antigens. These include the two trans-membrane (2TM) superfamily constituted by multicopy repetitive interspersed family (RIFINs), subtelomeric variable open reading frame (STEVORs) and Plasmodium falciparum Maurer's cleft two trans-membrane proteins present only in P. falciparum and some simian infecting Plasmodium species. Their hypervariable region flanked by 2TM domains exposed on the iRBC surface is believed to generate antigenic diversity. Though historically named "2TM superfamily," several A-type RIFINs and some STEVORs assume one trans-membrane topology. RIFINs and STEVORs share varied functions in different parasite life cycle stages like rosetting, alteration of iRBC rigidity and immune evasion. Additionally, a member of the STEVOR family has been implicated in merozoite invasion. Differential expression of these families in laboratory strains and clinical isolates propose them to be important for host cell survival and defense. The role of RIFINs in modulation of host immune response and presence of protective antibodies against these surface exposed molecules in patient sera highlights them as attractive targets of antimalarial therapies and vaccines. 2TM proteins are Plasmodium export elements positive, and several of these are exported to the infected erythrocyte surface after exiting through the classical secretory pathway within parasites. Cleaved and modified proteins are trafficked after packaging in vesicles to reach Maurer's clefts, while information regarding delivery to the iRBC surface is sparse. Expression and export timing of the RIFIN and Plasmodium falciparum erythrocyte membrane protein1 families correspond to each other. Here, we have compiled and comprehended detailed information regarding orthologues, domain architecture, surface topology, functions and trafficking of members of the "2TM superfamily." Considering the large repertoire of proteins included in the 2TM superfamily and recent advances defining their function in malaria biology, a surge in research carried out on this important protein superfamily is likely.

Genome-wide analysis of small nucleolar RNAs of Leishmania major reveals a rich repertoire of RNAs involved in modification and processing of rRNA

Trypanosomatids are protozoan parasites and the causative agent of infamous infectious diseases. These organisms regulate their gene expression mainly at the post-transcriptional level and possess characteristic RNA processing mechanisms. In this study, we analyzed the complete repertoire of Leishmania major small nucleolar (snoRNA) RNAs by performing RNA-seq analysis on RNAs that were affinity-purified using the C/D snoRNA core protein, SNU13, and the H/ACA core protein, NHP2. This study revealed a large collection of C/D and H/ACA snoRNAs, organized in gene clusters generally containing both snoRNA types. Abundant snoRNAs were identified and predicted to guide trypanosome-specific rRNA cleavages. The repertoire of snoRNAs was compared to that of the closely related Trypanosoma brucei, and 80% of both C/D and H/ACA molecules were found to have functional homologues. The comparative analyses elucidated how snoRNAs evolved to generate molecules with analogous functions in both species. Interestingly, H/ACA RNAs have great flexibility in their ability to guide modifications, and several of the RNA species can guide more than one modification, compensating for the presence of single hairpin H/ACA snoRNA in these organisms. Placing the predicted modifications on the rRNA secondary structure revealed hypermodification regions mostly in domains which are modified in other eukaryotes, in addition to trypanosome-specific modifications.

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.

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.

Directional gene expression and antisense transcripts in sexual and asexual stages of Plasmodium falciparum

BACKGROUND

It has been shown that nearly a quarter of the initial predicted gene models in the Plasmodium falciparum genome contain errors. Although there have been efforts to obtain complete cDNA sequences to correct the errors, the coverage of cDNA sequences on the predicted genes is still incomplete, and many gene models for those expressed in sexual or mosquito stages have not been validated. Antisense transcripts have widely been reported in P. falciparum; however, the extent and pattern of antisense transcripts in different developmental stages remain largely unknown.

RESULTS

We have sequenced seven bidirectional libraries from ring, early and late trophozoite, schizont, gametocyte II, gametocyte V, and ookinete, and four strand-specific libraries from late trophozoite, schizont, gametocyte II, and gametocyte V of the 3D7 parasites. Alignment of the cDNA sequences to the 3D7 reference genome revealed stage-specific antisense transcripts and novel intron-exon splicing junctions. Sequencing of strand-specific cDNA libraries suggested that more genes are expressed in one direction in gametocyte than in schizont. Alternatively spliced genes, antisense transcripts, and stage-specific expressed genes were also characterized.

CONCLUSIONS

It is necessary to continue to sequence cDNA from different developmental stages, particularly those of non-erythrocytic stages. The presence of antisense transcripts in some gametocyte and ookinete genes suggests that these antisense RNA may play an important role in gene expression regulation and parasite development. Future gene expression studies should make use of directional cDNA libraries. Antisense transcripts may partly explain the observed discrepancy between levels of mRNA and protein expression.

3' processing in protists

Molecular biologists have traditionally focused on the very small corner of eukaryotic evolution that includes yeast and animals; even plants have been neglected. In this article, we describe the scant information that is available concerning RNA processing in the other four major eukaryotic groups, especially pathogenic protists. We focus mainly on polyadenylation and nuclear processing of stable RNAs. These processes have--where examined--been shown to be conserved, but there are many novel details. We also briefly mention other processing reactions such as splicing.

New insights into the blood-stage transcriptome of Plasmodium falciparum using RNA-Seq

Recent advances in high-throughput sequencing present a new opportunity to deeply probe an organism's transcriptome. In this study, we used Illumina-based massively parallel sequencing to gain new insight into the transcriptome (RNA-Seq) of the human malaria parasite, Plasmodium falciparum. Using data collected at seven time points during the intraerythrocytic developmental cycle, we (i) detect novel gene transcripts; (ii) correct hundreds of gene models; (iii) propose alternative splicing events; and (iv) predict 5' and 3' untranslated regions. Approximately 70% of the unique sequencing reads map to previously annotated protein-coding genes. The RNA-Seq results greatly improve existing annotation of the P. falciparum genome with over 10% of gene models modified. Our data confirm 75% of predicted splice sites and identify 202 new splice sites, including 84 previously uncharacterized alternative splicing events. We also discovered 107 novel transcripts and expression of 38 pseudogenes, with many demonstrating differential expression across the developmental time series. Our RNA-Seq results correlate well with DNA microarray analysis performed in parallel on the same samples, and provide improved resolution over the microarray-based method. These data reveal new features of the P. falciparum transcriptional landscape and significantly advance our understanding of the parasite's red blood cell-stage transcriptome.

Assessing functional annotation transfers with inter-species conserved coexpression: application to Plasmodium falciparum

BACKGROUND

Plasmodium falciparum is the main causative agent of malaria. Of the 5 484 predicted genes of P. falciparum, about 57% do not have sufficient sequence similarity to characterized genes in other species to warrant functional assignments. Non-homology methods are thus needed to obtain functional clues for these uncharacterized genes. Gene expression data have been widely used in the recent years to help functional annotation in an intra-species way via the so-called Guilt By Association (GBA) principle.

RESULTS

We propose a new method that uses gene expression data to assess inter-species annotation transfers. Our approach starts from a set of likely orthologs between a reference species (here S. cerevisiae and D. melanogaster) and a query species (P. falciparum). It aims at identifying clusters of coexpressed genes in the query species whose coexpression has been conserved in the reference species. These conserved clusters of coexpressed genes are then used to assess annotation transfers between genes with low sequence similarity, enabling reliable transfers of annotations from the reference to the query species. The approach was used with transcriptomic data sets of P. falciparum, S. cerevisiae and D. melanogaster, and enabled us to propose with high confidence new/refined annotations for several dozens hypothetical/putative P. falciparum genes. Notably, we revised the annotation of genes involved in ribosomal proteins and ribosome biogenesis and assembly, thus highlighting several potential drug targets.

CONCLUSIONS

Our approach uses both sequence similarity and gene expression data to help inter-species gene annotation transfers. Experiments show that this strategy improves the accuracy achieved when using solely sequence similarity and outperforms the accuracy of the GBA approach. In addition, our experiments with P. falciparum show that it can infer a function for numerous hypothetical genes.