Small RNAs with 5'-polyphosphate termini associate with a Piwi-related protein and regulate gene expression in the single-celled eukaryote Entamoeba histolytica

Work with me here: variations in genome content and emerging genetic tools in Entamoeba histolytica

Entamoeba histolytica is the causative agent of amoebiasis, a significant source of morbidity and mortality in developing nations. Despite this, E. histolytica is understudied, leading to few treatment options and a poor understanding of pathogenesis. Genetic tools have historically been limited. By applying modern approaches, it was recently revealed that the genome is aneuploid. Interestingly, gene expression levels do not correlate with ploidy, potentially highlighting the importance of RNAi in gene regulation. Characterization of the RNAi pathway has led to potent tools for targeted gene knockdown, and the advent of RNAi-based forward genetics. CRISPR/Cas tools for editing the endogenous genome are an exciting possibility on the horizon. We celebrate the gains that have made E. histolytica tractable and anticipate continued advances.

Avirulent UG10 Entamoeba histolytica mutant derived from HM-1:IMSS strain shows limited genome variability and aberrant 5-methyl cytosine genomic distribution

Entamoeba histolytica, an intestinal parasite of global significance, poses substantial health risks with its associated high morbidity and mortality rates. Despite the current repertoire of molecular tools for the study of gene function in, the regulatory mechanisms governing its pathogenicity remain largely unexplored. This knowledge gap underscores the need to elucidate key genetic determinants orchestrating cellular functions critical to its virulence. Previously, our group generated an avirulent strain, termed UG10, with the same genetic background as the HM1:IMSS strain. UG10 strain, despite showing normal expression levels of well-known virulence factors, was unable to perform in-vitro and in-vivo activities related to amoebic virulence. In this study, we aimed to uncover the genome-wide modifications that rendered the avirulent phenotype of the UG10 strain through whole-genome sequencing. As a complementary approach, we conducted Methylated DNA Immunoprecipitation coupled with sequencing (MeDIP-seq) analysis on both the highly virulent HM1:IMSS strain and the low-virulence UG10 strain to uncover the genome-wide methylation profile. These dual methodologies revealed two aspects of the UG10 avirulent strain. One is the random integration of fragments from the ribosomal gene cluster and tRNA genes, ranging from 120 to 400 bp; and secondly, a clear, enriched methylation profile in the coding and non-coding strand relative to the start codon sequence in genes encoding small GTPases, which is associated with the previously described avirulent phenotype. This study provides the foundation to explore other genetic and epigenetic regulatory circuitries in E. histolytica and novel targets to understand the pathogenic mechanism of this parasite.

Extracellular Vesicles and Their Impact on the Biology of Protozoan Parasites

Extracellular vesicles (EVs) are lipid-membrane-bound structures produced naturally by all cells and have a variety of functions. EVs act as vehicles for transporting important molecular signals from one cell to another. Several parasites have been shown to secrete EVs, and their biological functions have been extensively studied. EVs have been shown to facilitate communication with the host cells (such as modulation of the host's immune system or promoting attachment and invasion into the host cells) or for communication between parasitic cells (e.g., transferring drug-resistance genes or factors modulating stage conversion). It is clear that EVs play an important role in host-parasite interactions. In this review, we summarized the latest research on the EVs secreted by protozoan parasites and their role in host-parasite and parasite-parasite communications.

Stress Response in Entamoeba histolytica Is Associated with Robust Processing of tRNA to tRNA Halves

tRNA-derived fragments have been reported in many different organisms and have diverse cellular roles, such as regulating gene expression, inhibiting protein translation, silencing transposable elements, and modulating cell proliferation. In particular, tRNA halves, a class of tRNA fragments produced by the cleavage of tRNAs in the anti-codon loop, have been widely reported to accumulate under stress and regulate translation in cells. Here, we report the presence of tRNA-derived fragments in Entamoeba, with tRNA halves being the most abundant. We further established that tRNA halves accumulate in the parasites upon different stress stimuli such as oxidative stress, heat shock, and serum starvation. We also observed differential expression of tRNA halves during developmental changes of trophozoite-to-cyst conversion, with various tRNA halves accumulating during early encystation. In contrast to other systems, the stress response does not appear to be mediated by a few specific tRNA halves, as multiple tRNAs appear to be processed during the various stresses. Furthermore, we identified some tRNA-derived fragments associated with Entamoeba Argonaute proteins, EhAgo2-2 and EhAgo2-3, which have a preference for different tRNA-derived fragment species. Finally, we show that tRNA halves are packaged inside extracellular vesicles secreted by amoebas. The ubiquitous presence of tRNA-derived fragments, their association with the Argonaute proteins, and the accumulation of tRNA halves during multiple different stresses, including encystation, suggest a nuanced level of gene expression regulation mediated by different tRNA-derived fragments in Entamoeba. IMPORTANCE In the present study, we report for the first time the presence of tRNA-derived fragments in Entamoeba. tRNA-derived fragments were identified by bioinformatics analyses of small-RNA sequencing data sets from the parasites and also confirmed experimentally. We found that tRNA halves accumulated in parasites exposed to environmental stress or during the developmental process of encystation. We also found that shorter tRNA-derived fragments are bound to Entamoeba Argonaute proteins, indicating that they may have a potential role in the Argonaute-mediated RNA-interference pathway, which mediates robust gene silencing in Entamoeba. We noticed that in response to heat shock, the protein translation levels were elevated in the parasites. This effect was reversed in the presence of an analog of leucine, which also reduced the levels of the tRNA halves in the stressed cells. Our results suggest that tRNA-derived fragments in Entamoeba have a possible role in regulating gene expression during environmental stress.

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.

The non-LTR retrotransposons of Entamoeba histolytica: genomic organization and biology

Genome sequence analysis of Entamoeba species revealed various classes of transposable elements. While E. histolytica and E. dispar are rich in non-long terminal repeat (LTR) retrotransposons, E. invadens contains predominantly DNA transposons. Non-LTR retrotransposons of E. histolytica constitute three families of long interspersed nuclear elements (LINEs), and their short, nonautonomous partners, SINEs. They occupy ~ 11% of the genome. The EhLINE1/EhSINE1 family is the most abundant and best studied. EhLINE1 is 4.8 kb, with two ORFs that encode functions needed for retrotransposition. ORF1 codes for the nucleic acid-binding protein, and ORF2 has domains for reverse transcriptase (RT) and endonuclease (EN). Most copies of EhLINEs lack complete ORFs. ORF1p is expressed constitutively, but ORF2p is not detected. Retrotransposition could be demonstrated upon ectopic over expression of ORF2p, showing that retrotransposition machinery is functional. The newly retrotransposed sequences showed a high degree of recombination. In transcriptomic analysis, RNA-Seq reads were mapped to individual EhLINE1 copies. Although full-length copies were transcribed, no full-length 4.8 kb transcripts were seen. Rather, sense transcripts mapped to ORF1, RT and EN domains. Intriguingly, there was strong antisense transcription almost exclusively from the RT domain. These unique features of EhLINE1 could serve to attenuate retrotransposition in E. histolytica.

Establishment of quantitative RNAi-based forward genetics in Entamoeba histolytica and identification of genes required for growth

While Entamoeba histolytica remains a globally important pathogen, it is dramatically understudied. The tractability of E. histolytica has historically been limited, which is largely due to challenging features of its genome. To enable forward genetics, we constructed and validated the first genome-wide E. histolytica RNAi knockdown mutant library. This library allows for Illumina deep sequencing analysis for quantitative identification of mutants that are enriched or depleted after selection. We developed a novel analysis pipeline to precisely define and quantify gene fragments. We used the library to perform the first RNAi screen in E. histolytica and identified slow growth (SG) mutants. Among genes targeted in SG mutants, many had annotated functions consistent with roles in cellular growth or metabolic pathways. Some targeted genes were annotated as hypothetical or lacked annotated domains, supporting the power of forward genetics in uncovering functional information that cannot be gleaned from databases. While the localization of neither of the proteins targeted in SG1 nor SG2 mutants could be predicted by sequence analysis, we showed experimentally that SG1 localized to the cytoplasm and cell surface, while SG2 localized to the cytoplasm. Overexpression of SG1 led to increased growth, while expression of a truncation mutant did not lead to increased growth, and thus aided in defining functional domains in this protein. Finally, in addition to establishing forward genetics, we uncovered new details of the unusual E. histolytica RNAi pathway. These studies dramatically improve the tractability of E. histolytica and open up the possibility of applying genetics to improve understanding of this important pathogen.

RISC in Entamoeba histolytica: Identification of a Protein-Protein Interaction Network for the RNA Interference Pathway in a Deep-Branching Eukaryote

Entamoeba histolytica is a protozoan parasite that causes amebiasis in humans and is a major health concern in developing countries. Our previous work revealed a functional RNA interference (RNAi) pathway in Entamoeba. Several unusual features encompass the RNAi pathway in the parasite, including small RNAs (sRNAs) with a 5'-polyphosphate structure (identified to date only in Entamoeba and nematodes) and the conspicuous absence of a canonical Dicer enzyme. Currently, little is known about the Entamoeba RNA-induced silencing complex (RISC), which is critical in understanding how RNAi is achieved in the parasite. In this study, we examined the RISC of EhAgo2-2, the most highly expressed Argonaute protein in Entamoeba. We identified 43 protein components of EhAgo2-2 RISC with a broad range of functional activities. Two proteins with nucleosome assembly protein (NAP) domains, not previously observed in other RNAi systems, were identified as novel core members of amebic RISC. We further demonstrated the interaction of these NAPs with Ago using an in vitro recombinant system. Finally, we characterized the interaction network of five RISC components identified in this study to further elucidate the interactions of these RNAi pathway proteins. Our data suggest the presence of closely interacting protein groups within RISC and allowed us to build a map of protein-protein interactions in relation to Ago. Our work is the first to elucidate RISC components in Entamoeba and expands the current knowledge of RISC to a deep-branching single-celled eukaryote. IMPORTANCE Entamoeba histolytica is a leading parasitic cause of death in developing countries, and our efforts are focused on defining the molecular basis of RNA interference (RNAi) gene regulation in this parasite. The Entamoeba RNAi pathway effectively silences a subset of endogenous genes and has also been harnessed as a gene silencing tool to study gene function in this organism. However, little is known about the components of the Entamoeba RNA-induced silencing complex (RISC), which is critical in understanding how gene silencing is achieved in the parasite. This study characterizes, for the first time, the RISC components in Entamoeba and provides new insights in understanding the molecular regulatory mechanisms of RNAi in this parasite, including the demonstration of novel Ago protein-interacting partners. From an evolutionary point of view, our findings expand the current knowledge of RISC to a deep-branching single-celled eukaryote.

Development of a CRISPR/Cas9 system in Entamoeba histolytica: proof of concept

The protozoan parasite Entamoeba histolytica is an important human pathogen and a leading parasitic cause of death on a global scale. The lack of molecular tools for genome editing hinders the study of important biological functions of this parasite. Due to its versatility, the CRISPR (clustered regularly interspaced short palindromic repeat)-Cas9 system has been successfully used to induce site-specific genomic alterations, including in protozoan parasites. In this study, we optimised CRISPR-Cas9 for use as a genetic tool in E. histolytica. We chose a single plasmid approach containing both guide RNA (gRNA) and Cas9 nuclease expression cassettes. The amebic U6 promoter was used to drive the expression of the gRNA and its expression was confirmed by Northern blot analysis. Stable transfectant cell lines were obtained using a destabilising domain of dihydrofolate reductase fused to myc-tagged Cas9 (ddCas9). With this system, we were able to induce ddCas9 expression 16 h following treatment with the small molecule ligand trimethoprim (TMP). Stable cell lines expressing ddCas9 and Luc-gRNA or non-specific (NS)-gRNA were transiently transfected with a plasmid containing a mutated luciferase gene (pDeadLuc) targeted by Luc-gRNA and another plasmid with a truncated luciferase gene (pDonorLuc) to restore luciferase expression and consequent activity. We observed that luminescence signal increased for the cell line expressing Luc-gRNA, suggesting that homologous recombination was facilitated by Cas9 activity. This evidence is supported by the presence of chimeric DNA detected by PCR and confirmed by sequencing of the resulting repaired DNA obtained by homologous recombination. We believe this represents the first report of a CRISPR/Cas9 system use in Entamoeba and provides evidence that this genome editing approach can be useful for genetic studies in this early branching eukaryote.

Transcriptomic analysis of Entamoeba histolytica reveals domain-specific sense strand expression of LINE-encoded ORFs with massive antisense expression of RT domain

LINEs are retrotransposable elements found in diverse organisms. Their activity is kept in check by several mechanisms, including transcriptional silencing. Here we have analyzed the transcription status of LINE1 copies in the early-branching parasitic protist Entamoeba histolytica. Full-length EhLINE1 encodes ORF1, and ORF2 with reverse transcriptase (RT) and endonuclease (EN) domains. RNA-Seq analysis of EhLINE1 copies (both truncated and full-length) showed unique features. Firstly, although 20/41 transcribed copies were full-length, we failed to detect any full-length transcripts. Rather, sense-strand transcripts mapped to the functional domains- ORF1, RT and EN. Secondly, there was strong antisense transcription specifically from RT domain. No antisense transcripts were seen from ORF1. Antisense RT transcripts did not encode known functional peptides. They could possibly be involved in attenuating translation of RT domain, as we failed to detect ORF2p, whereas ORF1p was detectable. Lack of full-length transcripts and strong antisense RT expression may serve to limit EhLINE1 retrotransposition.

Small RNAs Are Implicated in Regulation of Gene and Transposable Element Expression in the Protist Trichomonas vaginalis

Trichomoniasis, caused by the protozoan Trichomonas vaginalis, is the most common nonviral sexually transmitted infection in humans. The millions of cases each year have sequelae that may include complications during pregnancy and increased risk of HIV infection.

Trichomonas vaginalis is the causative agent of trichomoniasis, the most prevalent nonviral sexually transmitted infection worldwide. Repetitive elements, including transposable elements (TEs) and virally derived repeats, comprise more than half of the ∼160-Mb T. vaginalis genome. An intriguing question is how the parasite controls its potentially lethal complement of mobile elements, which can disrupt transcription of protein-coding genes and genome functions. In this study, we generated high-throughput RNA sequencing (RNA-Seq) and small RNA-Seq data sets in triplicate for the T. vaginalis G3 reference strain and characterized the mRNA and small RNA populations and their mapping patterns along all six chromosomes. Mapping the RNA-Seq transcripts to the genome revealed that the majority of genes predicted within repetitive elements are not expressed. Interestingly, we identified a novel species of small RNA that maps bidirectionally along the chromosomes and is correlated with reduced protein-coding gene expression and reduced RNA-Seq coverage in repetitive elements. This novel small RNA family may play a regulatory role in gene and repetitive element expression. Our results identify a possible small RNA pathway mechanism by which the parasite regulates expression of genes and TEs and raise intriguing questions as to the role repeats may play in shaping T. vaginalis genome evolution and the diversity of small RNA pathways in general.

IMPORTANCE Trichomoniasis, caused by the protozoan Trichomonas vaginalis, is the most common nonviral sexually transmitted infection in humans. The millions of cases each year have sequelae that may include complications during pregnancy and increased risk of HIV infection. Given its evident success in this niche, it is paradoxical that T. vaginalis harbors in its genome thousands of transposable elements that have the potential to be extremely detrimental to normal genomic function. In many organisms, transposon expression is regulated by the activity of endogenously expressed short (∼21 to 35 nucleotides [nt]) small RNA molecules that effect gene silencing by targeting mRNAs for degradation or by recruiting epigenetic silencing machinery to locations in the genome. Our research has identified small RNA molecules correlated with reduced expression of T. vaginalis genes and transposons. This suggests that a small RNA pathway is a major contributor to gene expression patterns in the parasite and opens up new avenues for investigation into small RNA biogenesis, function, and diversity.

Possible role played by the SINE2 element in gene regulation, as demonstrated by differential processing and polyadenylation in avirulent strains of E. histolytica

Entamoeba histolytica represents a useful model in parasitic organisms due to its complex genomic organization and survival mechanisms. To counteract pathogenic organisms, it is necessary to characterize their molecular biology to design new strategies to combat them. In this report, we investigated a less-known genetic element, short interspersed nuclear element 2 (SINE2), that is present in this ameba and is highly transcribed and polyadenylated. In this study, we show that in two different nonvirulent strains of E. histolytica, SINE2 is differentially processed into two transcript fragments, that is, a full-length 560-nt fragment and a shorter 393-nt fragment bearing an approximately 18-nt polyadenylation tail. Sequence analysis of the SINE2 transcript showed that a Musashi-like protein may bind to it. Also, two putative Musashi-like sequences were identified on the transcript. Semiquantitative expression analysis of the two Musashi-like proteins identified in the E. histolytica genome (XP_648918 and XP_649094) showed that XP_64094 is overexpressed in the nonvirulent strains tested. The information available in the literature and the results presented in this report indicate that SINE2 may affect other genes, as observed with the epigenetic silencing of the G3 strain, by an antisense mechanism or via RNA-protein interactions that may ultimately be involved in the phenotype of nonvirulent strains of E. histolytica.

Identification of oligo-adenylated small RNAs in the parasite Entamoeba and a potential role for small RNA control

Background

The RNA interference (RNAi) pathway is a gene regulation mechanism that utilizes small RNA (sRNA) and Argonaute (Ago) proteins to silence target genes. Our previous work identified a functional RNAi pathway in the protozoan parasite Entamoeba histolytica, including abundant 27 nt antisense sRNA populations which associate with EhAgo2–2 protein. However, there is lack of understanding about the sRNAs that are bound to two other EhAgos (EhAgo2–1 and 2–3), and the mechanism of sRNA regulation itself is unclear in this parasite. Therefore, identification of the entire pool of sRNA species and their sub-populations that associate with each individual EhAgo protein would be a major step forward.

Results

In the present study, we sequenced sRNA libraries from both total RNAs and EhAgo bound RNAs. We identified a new population of 31 nt sRNAs that results from the addition of a non-templated 3–4 adenosine nucleotides at the 3′-end of the 27 nt sRNAs, indicating a non-templated RNA-tailing event in the parasite. The relative abundance of these two sRNA populations is linked to the efficacy of gene silencing for the target gene when parasites are transfected with an RNAi-trigger construct, indicating that non-templated sRNA-tailing likely play a role in sRNA regulation in this parasite. We found that both sRNA populations (27 nt and 31 nt) are present in the related parasite Entamoeba invadens, and are unchanged during the development. In sequencing the sRNAs associating with the three EhAgo proteins, we observed that despite distinct cellular localization, all three EhAgo sRNA libraries contain 27 nt sRNAs with 5′-polyphosphate (5′-polyP) structure and share a largely overlapping sRNA repertoire. In addition, our data showed that a fraction of 31 nt sRNAs associate with EhAgo2–2 but not with its mutant protein (C-terminal deletion), nor other two EhAgos, indicating a specific EhAgo site may be required for sRNA modification process in the parasite.

Conclusion

We identified a new population of sRNA with non-templated oligo-adenylation modification, which is the first such observation amongst single celled protozoan parasites. Our sRNA sequencing libraries provide the first comprehensive sRNA dataset for all three Entamoeba Ago proteins, which can serve as a useful database for the amoeba community.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07275-6.

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.

Characterization of Extracellular Vesicles from Entamoeba histolytica Identifies Roles in Intercellular Communication That Regulates Parasite Growth and Development

Extracellular vesicles (EVs) secreted by eukaryotic and prokaryotic cells to transport lipids, proteins, and nucleic acids to the external environment have important roles in cell-cell communication through cargo transfer. We identified and characterized EVs from Entamoeba histolytica, a protozoan parasite and a human pathogen. Conditioned medium from amebic parasites contained particles consistent with the expected size and morphology of EVs. Mass spectrometry was used to characterize the EV proteome and showed that it was enriched in common exosome marker proteins, including proteins associated with vesicle formation, cell signaling, and metabolism, as well as cytoskeletal proteins. Additionally, the EVs were found to selectively package small RNAs (sRNA), which were protected within the vesicles against RNase treatment. Sequencing analysis of the sRNA contained in EVs revealed that the majority were 27 nucleotides (nt) in size and represented a subset of the cellular antisense small RNA population that has previously been characterized in Entamoeba RNA interference (RNAi) pathway proteins, including Argonaute, were also present in amebic EVs. Interestingly, we found that the amebic EVs impacted intercellular communication between parasites and altered encystation efficiency. EVs isolated from encysting parasites promoted encystation in other parasites, whereas EVs from metabolically active trophozoites impeded encystation. Overall, the data reveal that Entamoeba secrete EVs that are similar in size and shape to previously characterized exosomes from other organisms and that these EVs contain a defined protein and small RNA cargo and have roles in intercellular communication among parasites and influence growth kinetics.

The NAD+ Responsive Transcription Factor ERM-BP Functions Downstream of Cellular Aggregation and Is an Early Regulator of Development and Heat Shock Response in Entamoeba

Entamoeba histolytica is a protozoan parasite and a major cause of dysentery and diarrheal disease in developing countries. Disease transmission from one host to another occurs via cysts which can survive in environmental extremes and are transmitted through contaminated food and water. Recent studies in our lab identified a novel transcription factor, Encystation Regulatory Motif- Binding Protein (ERM-BP), which is responsive to NAD+ and has an important role in encystation. The key residues important for ERM-BP function were demonstrated in vitro using recombinant protein. In this study we demonstrate the in vivo functional consequences of mutations in key domains and their impact on Entamoeba encystation. Our results show that mutations in the DNA binding domain (ERM-BP-DBM) and in the nicotinamidase domain (ERM-BP-C198A) lead to protein mis-localization in both trophozoites and cysts and significantly reduce encystation efficiency. Additionally, we showed that silencing of ERM-BP significantly decreased the size and number of multi-nucleated giant cells (MGC) that form during encystation, indicating that ERM-BP functions upstream of the cellular aggregation that precedes stage conversion. Dissection of epistatic interactions between ERM-BP and a second encystation-related transcription factor, NF-Y revealed that ERM-BP is upstream of NF-Y in controlling the developmental cascade and appears to be one of the earliest regulators of development identified to date in Entamoeba. We also demonstrated that ERM-BP is upregulated during heat stress in Entamoeba, another condition which increases intracellular NAD+ levels and that overexpression of ERM-BP makes E. histolytica and E. invadens parasites more resistant to heat stress. Overexpression of ERM-BP in E. histolytica also induced the formation of cyst-like quadrinucleated cells and formation of MGCs. Overall, our work has identified an important role of ERM-BP in Entamoeba stress response and links an NAD+-responsive transcription factor to both development and heat shock response. Characterization of stress and developmental cascades are important avenues to investigate for Entamoeba, an important human parasitic pathogen.

Functional Characterization of Entamoeba histolytica Argonaute Proteins Reveals a Repetitive DR-Rich Motif Region That Controls Nuclear Localization

The RNA interference (RNAi) pathway regulates gene expression in many eukaryotic organisms. Argonaute (Ago) proteins, together with bound small RNAs (sRNAs), are key effectors that mediate gene silencing function. However, there is limited knowledge of Ago proteins and their functions in nonmodel systems. In the protozoan parasite Entamoeba histolytica, RNAi is a robust means for stable gene silencing mediated via large populations of antisense sRNAs. Here, we report functional characterization of three Ago proteins in E. histolytica (EhAgo2-1, EhAgo2-2, and EhAgo2-3). Our data show that each EhAgo protein has a distinct subcellular localization and binds 27-nucleotide (nt) sRNAs and that the localization of EhAgo proteins is altered in response to stress conditions. Via mutagenesis analyses, we demonstrated that the Ago PAZ (Piwi/Argonaute/Zwille) domain in all three EhAgos is essential for sRNA binding. With mutation of the PAZ domain in EhAgo2-2, there was no effect on the nuclear localization of the protein but a strong phenotype and a growth defect. We further show that EhAgo2-2 contains an unusual repetitive DR-rich (aspartic acid, arginine-rich) motif region which functions as a nuclear localization signal (NLS) and is both necessary and sufficient to mediate nuclear localization. Overall, our data delineate the localization and sRNA binding features of the three E. histolytica Ago proteins and demonstrate that the PAZ domain is necessary for sRNA binding. The repetitive DR-rich motif region in EhAgo2-2 has not previously been defined in other systems, which adds to the novel observations that can be made when studies of the RNAi pathway are extended to nonmodel systems.IMPORTANCE The protozoan parasite Entamoeba histolytica, which causes amebiasis and affects over 50 million people worldwide, contains an important RNAi pathway for gene silencing. Gene silencing via the RNAi pathway is mediated by the Argonaute (Ago) proteins. However, we lack knowledge on Ago function(s) in this nonmodel system. In this paper, we discovered that three E. histolytica Ago proteins (EhAgo2-1, EhAgo2-2, and EhAgo2-3) all bind 27-nt small RNAs and have distinct subcellular localizations, which change in response to stress conditions. The EhAgos bind small RNA populations via their PAZ domains. An unusual repetitive DR-rich motif region is identified in EhAgo2-2 that functions as a nuclear localization signal. Our results show for the first time an active nuclear transport process of the EhAgo2-2 RNA-induced silencing complex (RISC) in this parasite. These data add to the novel observations that can be made when studies of the RNAi pathway are extended to nonmodel systems.

Trigger-induced RNAi gene silencing to identify pathogenicity factors of Entamoeba histolytica

Recently, Entamoeba histolytica clones derived from isolate HM-1:IMSS that differ in their pathogenicity were identified. Whereas some clones induce amoebic liver abscesses (ALAs) in animal models of amoebiasis, others provoke only minimal liver lesions. Based on transcriptome studies of pathogenic and nonpathogenic clones, differentially expressed genes associated with reduced or increased liver pathology can be identified. Here, to analyze the influence of these genes on ALA formation in more detail, an RNA interference-trigger mediated silencing approach was used. Using newly identified trigger sequences, the expression of 15 genes was silenced. The respective transfectants were analyzed for their ability to induce liver destruction in the murine model for the disease. Silencing of EHI_180390 (encoding an AIG1 protein) increased liver pathology induced by a nonpathogenic parent clone, whereas silencing of EHI_127670 (encoding a hypothetical protein) decreased the pathogenicity of an initially pathogenic parent clone. Additional phenotypical in vitro analyses of EHI_127670 silencing as well as overexpression transfectants indicated that this molecule has an influence on size, growth, and cysteine peptidase activity of E. histolytica. This work describes an example of how the sole operational method for effective gene silencing in E. histolytica can be used for comprehensive analyses of putative pathogenicity factors.-Matthiesen, J., Lender, C., Haferkorn, A., Fehling, H., Meyer, M., Matthies, T., Tannich, E., Roeder, T., Lotter, H., Bruchhaus, I. Trigger-induced RNAi gene silencing to identify pathogenicity factors of Entamoeba histolytica.

Life and Death of mRNA Molecules in Entamoeba histolytica

In eukaryotic cells, the life cycle of mRNA molecules is modulated in response to environmental signals and cell-cell communication in order to support cellular homeostasis. Capping, splicing and polyadenylation in the nucleus lead to the formation of transcripts that are suitable for translation in cytoplasm, until mRNA decay occurs in P-bodies. Although pre-mRNA processing and degradation mechanisms have usually been studied separately, they occur simultaneously and in a coordinated manner through protein-protein interactions, maintaining the integrity of gene expression. In the past few years, the availability of the genome sequence of Entamoeba histolytica, the protozoan parasite responsible for human amoebiasis, coupled to the development of the so-called "omics" technologies provided new opportunities for the study of mRNA processing and turnover in this pathogen. Here, we review the current knowledge about the molecular basis for splicing, 3' end formation and mRNA degradation in amoeba, which suggest the conservation of events related to mRNA life throughout evolution. We also present the functional characterization of some key proteins and describe some interactions that indicate the relevance of cooperative regulatory events for gene expression in this human parasite.

Structural and functional studies of a noncanonical Dicer from Entamoeba histolytica

RNaseIII proteins are dsRNA-specific endonucleases involved in many important biological processes, such as small RNA processing and maturation in eukaryotes. Various small RNAs have been identified in a protozoan parasite Entamoeba histolytica. EhRNaseIII is the only RNaseIII endonuclease domain (RIIID)-containing protein in E. histolytica. Here, we present three crystal structures that reveal several unique structural features of EhRNaseIII, especially the interactions between the two helixes (α1 and α7) flanking the RIIID core domain. Structure and sequence analysis indicate that EhRNaseIII is a noncanonical Dicer and it lacks a dsRBD in the C-terminal region (CTR). In vitro studies suggest that EhRNaseIII prefers to bind and cleave longer dsRNAs, generating products around 25 nucleotides in length. Truncation of the CTR or attaching the dsRBD of Aquifex aeolicus RNaseIII can enhance the binding and cleavage activities of EhRNaseIII. In combination with in vitro crosslinking assay, our results suggested that EhRNaseIII functions in a cooperative mode. We speculate that some partner proteins may exist in E. histolytica and regulates the activity of EhRNaseIII through interaction with its CTR. Our studies support that EhRNaseIII plays an important role in producing small RNAs in E. histolytica.

Recent advances in Entamoeba biology: RNA interference, drug discovery, and gut microbiome

In recent years, substantial progress has been made in understanding the molecular and cell biology of the human parasite Entamoeba histolytica, an important pathogen with significant global impact. This review outlines some recent advances in the Entamoeba field in the last five years, focusing on areas that have not recently been discussed in detail: (i) molecular mechanisms regulating parasite gene expression, (ii) new efforts at drug discovery using high-throughput drug screens, and (iii) the effect of gut microbiota on amoebiasis.

Serum stress responsive gene EhslncRNA of Entamoeba histolytica is a novel long noncoding RNA

Non coding RNAs are known to play important roles in regulating gene expression at the transcriptional and posttranscriptional levels in metazoans. There is very little information available about non coding RNAs in protists such as Entamoeba histolytica. Antisense and micro RNAs have been reported in E. histolytica, however no long non coding RNAs has been reported yet. Here, we report our findings on an in vitro serum stress-inducible gene EhslncRNA, a member of B1 transmembrane kinase family of E. histolytica. EhslncRNA encodes a transcript of 2.6 kb and sequence analysis revealed that there is no ORF >150 bp within this transcript. The transcript was found to be polyadenylated and mainly associated with monosomes in the cytoplasm under serum starvation. In normal proliferating cells this RNA is mainly present in the nucleus. The promoter element was mapped between 437 to 346 nucleotides upstream of transcriptional start site and has both positive and negative regulatory elements. Deletion of the negative element converted the promoter to serum inducible type. Oxygen and heat stress also increased expression levels of EhslncRNA. These observations suggest that EhslncRNA may be a long non coding RNA and likely to help cells withstand stressful conditions in the host.

Development of RNA Interference Trigger-Mediated Gene Silencing in Entamoeba invadens

Entamoeba histolytica, a protozoan parasite, is an important human pathogen and a leading parasitic cause of death. The organism has two life cycle stages, trophozoites, which are responsible for tissue invasion, and cysts, which are involved in pathogen transmission. Entamoeba invadens is the model system to study Entamoeba developmental biology, as high-grade regulated encystation and excystation are readily achievable. However, the lack of gene-silencing tools in E. invadens has limited the molecular studies that can be performed. Using the endogenous RNA interference (RNAi) pathway in Entamoeba, we developed an RNAi-based trigger gene-silencing approach inE. invadens We demonstrate that a gene's coding region that has abundant antisense small RNAs (sRNAs) can trigger silencing of a gene that is fused to it. The trigger fusion leads to the generation of abundant antisense sRNAs that map to the target gene, with silencing occurring independently of trigger location at the 5' or 3' end of a gene. Gene silencing is stably maintained during development, including encystation and excystation. We have used this approach to successfully silence two E. invadens genes: a putative rhomboid protease gene and a SHAQKY family Myb gene. The Myb gene is upregulated during oxidative stress and development, and its downregulation led, as predicted, to decreased viability under oxidative stress and decreased cyst formation. Thus, the RNAi trigger silencing method can be used to successfully investigate the molecular functions of genes inE. invadens Dissection of the molecular basis of Entamoeba stage conversion is now possible, representing an important technical advance for the system.

Technical advances in trigger-induced RNA interference gene silencing in the parasite Entamoeba histolytica

Entamoeba histolytica has a robust endogenous RNA interference (RNAi) pathway. There are abundant 27 nucleotide (nt) anti-sense small RNAs (AS sRNAs) that target genes for silencing and the genome encodes many genes involved in the RNAi pathway such as Argonaute proteins. Importantly, an E. histolytica gene with numerous AS sRNAs can function as a "trigger" to induce silencing of a gene that is fused to the trigger. Thus, the amebic RNAi pathway regulates gene expression relevant to amebic biology and has additionally been harnessed as a tool for genetic manipulation. In this study we have further improved the trigger-induced gene silencing method. We demonstrate that rather than using the full-length gene, a short portion of the coding region fused to a trigger is sufficient to induce silencing; the first 537 bp of the E. histolytica rhomboid gene (EhROM1) fused in-frame to the trigger was sufficient to silence EhROM1. We also demonstrated that the trigger method could silence two amebic genes concomitantly; fusion of the coding regions of EhROM1 and transcription factor, EhMyb, in-frame to a trigger gene resulted in both genes being silenced. Alternatively, two genes can be silenced sequentially: EhROM1-silenced parasites with no drug selection plasmid were transfected with trigger-EhMyb, resulting in parasites with both EhROM1 and EhMyb silenced. With all approaches tested, the trigger-mediated silencing was substantive and silencing was maintained despite loss of the G418 selectable marker. All gene silencing was associated with generation of AS sRNAs to the silenced gene. We tested the reversibility of the trigger system using inhibitors of histone modifications but found that the silencing was highly stable. This work represents a technical advance in the trigger gene silencing method in E. histolytica. Approaches that readily silence multiple genes add significantly to the genetic toolkit available to the ameba research community.

High Throughput Sequencing of Entamoeba 27nt Small RNA Population Reveals Role in Permanent Gene Silencing But No Effect on Regulating Gene Expression Changes during Stage Conversion, Oxidative, or Heat Shock Stress

The human parasite Entamoeba histolytica has an active RNA interference (RNAi) pathway with an extensive repertoire of 27nt small RNAs that silence genes. However the role of this pathway in regulating amebic biology remains unknown. In this study, we address whether silencing via 27nt small RNAs may be a mechanism for controlling gene expression changes during conversion between the trophozoite and cyst stages of the parasite. We sequenced small RNA libraries generated from trophozoites, early cysts, mature cysts, and excysting cells and mapped them to the E. invadens genome. Our results show that, as in E. histolytica, small RNAs in E. invadens are largely ~27nt in length, have an unusual 5'-polyphosphate structure and mediate gene silencing. However, when comparing the libraries from each developmental time-point we found few changes in the composition of the small RNA populations. Furthermore, genes targeted by small RNAs were permanently silenced with no changes in transcript abundance during development. Thus, the E. invadens 27nt small RNA population does not mediate gene expression changes during development. In order to assess the generalizability of our observations, we examined whether small RNAs may be regulating gene expression changes during stress response in E. histolytica. Comparison of the 27nt small RNA populations from E. histolytica trophozoites from basal conditions, or after heat shock or exposure to oxidative stress showed few differences. Similar to data in E. invadens development, genes targeted by small RNAs were consistently silenced and did not change expression under tested stress conditions. Thus, the biological roles of the 27nt small RNA population in Entamoeba remain elusive. However, as the first characterization of the RNAi pathway in E. invadens these data serve as a useful resource for the study of Entamoeba development and open the door to the development of RNAi-based gene silencing tools in E. invadens.

A Single RNaseIII Domain Protein from Entamoeba histolytica Has dsRNA Cleavage Activity and Can Help Mediate RNAi Gene Silencing in a Heterologous System

Dicer enzymes process double-stranded RNA (dsRNA) into small RNAs that target gene silencing through the RNA interference (RNAi) pathway. Dicer enzymes are complex, multi-domain RNaseIII proteins, however structural minimalism of this protein has recently emerged in parasitic and fungal systems. The most minimal Dicer, Saccharomyces castellii Dicer1, has a single RNaseIII domain and two double stranded RNA binding domains. In the protozoan parasite Entamoeba histolytica 27nt small RNAs are abundant and mediate silencing, yet no canonical Dicer enzyme has been identified. Although EhRNaseIII does not exhibit robust dsRNA cleavage in vitro, it can process dsRNA in the RNAi-negative background of Saccharomyces cerevisiae, and in conjunction with S. castellii Argonaute1 can partially reconstitute the RNAi pathway. Thus, although EhRNaseIII lacks the domain architecture of canonical or minimal Dicer enzymes, it has dsRNA processing activity that contributes to gene silencing via RNAi. Our data advance the understanding of small RNA biogenesis in Entamoeba as well as broaden the spectrum of non-canonical Dicer enzymes that contribute to the RNAi pathway.

Dimethylated H3K27 Is a Repressive Epigenetic Histone Mark in the Protist Entamoeba histolytica and Is Significantly Enriched in Genes Silenced via the RNAi Pathway

RNA interference (RNAi) is a fundamental biological process that plays a crucial role in regulation of gene expression in many organisms. Transcriptional gene silencing (TGS) is one of the important nuclear roles of RNAi. Our previous data show that Entamoeba histolytica has a robust RNAi pathway that links to TGS via Argonaute 2-2 (Ago2-2) associated 27-nucleotide small RNAs with 5'-polyphosphate termini. Here, we report the first repressive histone mark to be identified in E. histolytica, dimethylation of H3K27 (H3K27Me2), and demonstrate that it is enriched at genes that are silenced by RNAi-mediated TGS. An RNAi-silencing trigger can induce H3K27Me2 deposits at both episomal and chromosomal loci, mediating gene silencing. Our data support two phases of RNAi-mediated TGS: an active silencing phase where the RNAi trigger is present and both H3K27Me2 and Ago2-2 concurrently enrich at chromosomal loci; and an established silencing phase in which the RNAi trigger is removed, but gene silencing with H3K27Me2 enrichment persist independently of Ago2-2 deposition. Importantly, some genes display resistance to chromosomal silencing despite induction of functional small RNAs. In those situations, the RNAi-triggering plasmid that is maintained episomally gets partially silenced and has H3K27Me2 enrichment, but the chromosomal copy displays no repressive histone enrichment. Our data are consistent with a model in which H3K27Me2 is a repressive histone modification, which is strongly associated with transcriptional repression. This is the first example of an epigenetic histone modification that functions to mediate RNAi-mediated TGS in the deep-branching eukaryote E. histolytica.

RNAi pathway genes are resistant to small RNA mediated gene silencing in the protozoan parasite Entamoeba histolytica

The RNA interference pathway in the protist Entamoeba histolytica plays important roles in permanent gene silencing as well as in the regulation of virulence determinants. Recently, a novel RNA interference (RNAi)-based silencing technique was developed in this parasite that uses a gene endogenously silenced by small RNAs as a “trigger” to induce silencing of other genes that are fused to it. Fusion to a trigger gene induces the production of gene-specific antisense small RNAs, resulting in robust and permanent silencing of the cognate gene. This approach has silenced multiple genes including those involved in virulence and transcriptional regulation. We now demonstrate that all tested genes of the amebic RNAi pathway are unable to be silenced using the trigger approach, including Argonaute genes (Ago2-1, Ago2-2, and Ago2-3), RNaseIII, and RNA-dependent RNA polymerase (RdRP). In all situations (except for RdRP), fusion to a trigger successfully induces production of gene-specific antisense small RNAs to the cognate gene. These small RNAs are capable of silencing a target gene in trans, indicating that they are functional; despite this, however, they cannot silence the RNAi pathway genes. Interestingly, when a trigger is fused to RdRP, small RNA induction to RdRP does not occur, a unique phenotype hinting that either RdRP is highly resistant to being a target of small RNAs or that small RNA generation may be controlled by RdRP. The inability of the small RNA pathway to silence RNAi genes in E. histolytica, despite the generation of functional small RNAs to these loci suggest that epigenetic factors may protect certain genomic loci and thus determine susceptibility to small RNA mediated silencing.

Deciphering molecular mechanisms of mRNA metabolism in the deep-branching eukaryote Entamoeba histolytica

Although extraordinary rapid advance has been made in the knowledge of mechanisms regulating messenger RNA (mRNA) metabolism in mammals and yeast, little information is known in deep-branching eukaryotes. The complete genome sequence of Entamoeba histolytica, the protozoan parasite responsible for human amoebiasis, provided a lot of information for the identification and comparison of regulatory sequences and proteins potentially involved in mRNA synthesis, processing, and degradation. Here, we review the current knowledge of mRNA metabolism in this human pathogen. Several DNA motifs in promoter and nuclear factors involved in transcription, as well as conserved polyadenylation sequences in mRNA 3'-untranslated region and possible cleavage and polyadenylation factors, are described. In addition, we present recent data about proteins involved in mRNA decay with a special focus on the recently reported P-bodies in amoeba. Models for mechanisms of decapping and deadenylation-dependent pathways are discussed. We also review RNA-based gene silencing mechanisms and describe the DEAD/DExH box RNA helicases that are molecular players in all mRNA metabolism reactions. The functional characterization of selected proteins allows us to define a general framework to describe how mRNA synthesis, processing, and decay may occur in E. histolytica. Taken altogether, studies of mRNA metabolism in this single-celled eukaryotic model suggest the conservation of specific gene expression regulatory events through evolution.

Parallel evolution of genome structure and transcriptional landscape in the Epsilonproteobacteria

Background

Gene reshuffling, point mutations and horizontal gene transfer contribute to bacterial genome variation, but require the genome to rewire its transcriptional circuitry to ensure that inserted, mutated or reshuffled genes are transcribed at appropriate levels. The genomes of Epsilonproteobacteria display very low synteny, due to high levels of reshuffling and reorganisation of gene order, but still share a significant number of gene orthologs allowing comparison. Here we present the primary transcriptome of the pathogenic Epsilonproteobacterium Campylobacter jejuni, and have used this for comparative and predictive transcriptomics in the Epsilonproteobacteria.

Results

Differential RNA-sequencing using 454 sequencing technology was used to determine the primary transcriptome of C. jejuni NCTC 11168, which consists of 992 transcription start sites (TSS), which included 29 putative non-coding and stable RNAs, 266 intragenic (internal) TSS, and 206 antisense TSS. Several previously unknown features were identified in the C. jejuni transcriptional landscape, like leaderless mRNAs and potential leader peptides upstream of amino acid biosynthesis genes. A cross-species comparison of the primary transcriptomes of C. jejuni and the related Epsilonproteobacterium Helicobacter pylori highlighted a lack of conservation of operon organisation, position of intragenic and antisense promoters or leaderless mRNAs. Predictive comparisons using 40 other Epsilonproteobacterial genomes suggests that this lack of conservation of transcriptional features is common to all Epsilonproteobacterial genomes, and is associated with the absence of genome synteny in this subdivision of the Proteobacteria.

Conclusions

Both the genomes and transcriptomes of Epsilonproteobacteria are highly variable, both at the genome level by combining and division of multicistronic operons, but also on the gene level by generation or deletion of promoter sequences and 5′ untranslated regions. Regulatory features may have evolved after these species split from a common ancestor, with transcriptome rewiring compensating for changes introduced by genomic reshuffling and horizontal gene transfer.

Robust gene silencing mediated by antisense small RNAs in the pathogenic protist Entamoeba histolytica

RNA interference uses small RNAs (sRNA), which target genes for sequence-specific silencing. The parasite Entamoeba histolytica contains an abundant repertoire of 27 nt antisense (AS) sRNA with 5′-polyphosphate termini, but their roles in regulating gene expression have not been well established. We demonstrate that a gene-coding region to which large numbers of AS sRNAs map can serve as a ‘trigger’ and silence the gene fused to it. Silencing is mediated by generation of AS sRNAs with 5′-polyphosphate termini that have sequence specificity to the fused gene. The mechanism of silencing is independent of the placement of the trigger relative to the silenced gene but is dependent on the sRNA concentration to the trigger. Silencing requires transcription of the trigger-gene fusion and is maintained despite loss of the trigger plasmid. We used this approach to silence multiple amebic genes, including an E. histolytica Myb gene, which is upregulated during oxidative stress response. Silencing of the EhMyb gene decreased parasite viability under oxidative stress conditions. Thus, we have developed a new tool for genetic manipulation in E. histolytica with many advantages over currently available technologies. Additionally, these data shed mechanistic insights into a eukaryotic RNA interference pathway with many novel aspects.

The genome and transcriptome of the enteric parasite Entamoeba invadens, a model for encystation

Background

Several eukaryotic parasites form cysts that transmit infection. The process is found in diverse organisms such as Toxoplasma, Giardia, and nematodes. In Entamoeba histolytica this process cannot be induced in vitro, making it difficult to study. In Entamoeba invadens, stage conversion can be induced, but its utility as a model system to study developmental biology has been limited by a lack of genomic resources. We carried out genome and transcriptome sequencing of E. invadens to identify molecular processes involved in stage conversion.

Results

We report the sequencing and assembly of the E. invadens genome and use whole transcriptome sequencing to characterize changes in gene expression during encystation and excystation. The E. invadens genome is larger than that of E. histolytica, apparently largely due to expansion of intergenic regions; overall gene number and the machinery for gene regulation are conserved between the species. Over half the genes are regulated during the switch between morphological forms and a key signaling molecule, phospholipase D, appears to regulate encystation. We provide evidence for the occurrence of meiosis during encystation, suggesting that stage conversion may play a key role in recombination between strains.

Conclusions

Our analysis demonstrates that a number of core processes are common to encystation between distantly related parasites, including meiosis, lipid signaling and RNA modification. These data provide a foundation for understanding the developmental cascade in the important human pathogen E. histolytica and highlight conserved processes more widely relevant in enteric pathogens.

Identification and characterization of microRNAs from Entamoeba histolytica HM1-IMSS

BACKGROUND

Entamoeba histolytica is the causative agent of amebiasis, a disease that is a major source of morbidity and mortality in the developing world. MicroRNAs (miRNAs) are a large group of non-coding RNAs that play important roles in regulating gene expression and protein translation in animals. Genome-wide identification of miRNAs is a critical step to facilitating our understanding of genome organization, genome biology, evolution, and post-transcriptional regulation.

METHODOLOGY/PRINCIPAL FINDINGS

We sequenced a small RNA library prepared from a culture of trophozoites of Entamoeba histolytica Strain HM1-IMSS using a deep DNA sequencing approach. Deep sequencing yielded 16 million high-quality short sequence reads containing a total of 5 million non-redundant sequence reads. Based on a bioinformatics pipeline, we found that only 0.5% of these non-redundant small RNA reads were a perfect match with the drafted E. histolytica genome. We did not find miRNA homologs in plant or animal miRNAs. We discovered 199 new potential Entamoeba histolytica miRNAs. The expression and sequence of these Ehi-miRNAs were further validated through microarray by µParaflo Microfluidic Biochip Technology. Ten potential miRNAs were additionally confirmed by real time RT-PCR analysis. Prediction of target genes matched 32 known genes and 34 hypothetical genes.

CONCLUSIONS/SIGNIFICANCE

These results show that there is a number of regulatory miRNAs in Entamoeba histolytica. The collection of miRNAs in this parasite could be used as a new platform to study genomic structure, gene regulation and networks, development, and host-parasite interactions.

Comparative transcriptome analysis of small noncoding RNAs in different stages of Trypanosoma brucei

Trypanosoma brucei, a pathogen of human and domestic animals, is an early evolved parasitic protozoan with a complex life cycle. Most genes of this parasite are post-transcriptionally regulated. However, the mechanisms and the molecules involved remain largely unknown. We have deep-sequenced the small RNAs of two life stages of this parasite--the bloodstream form and the procyclic form. Our results show that the small RNAs of T. brucei could derive from multiple sources, including NATs (natural antisense transcripts), tRNAs, and rRNAs. Most of these small RNAs in the two stages were found to share uniform characteristics. However, our results demonstrate that their variety and expression show significant differences between different stages, indicating possible functional differentiation. Dicer-knockdown evidence further proved that some of the small interfering RNAs (siRNAs) could regulate the expression of genes. Based on the genome-wide analysis of the small RNAs in the two stages of T. brucei, our results not only provide evidence to study their differentiation but also shed light on questions regarding the origins and evolution of small RNA-based mechanisms in early eukaryotes.

Piwi-interacting RNAs protect DNA against loss during Oxytricha genome rearrangement

Genome duality in ciliated protozoa offers a unique system to showcase their epigenome as a model of inheritance. In Oxytricha, the somatic genome is responsible for vegetative growth, whereas the germline contributes DNA to the next sexual generation. Somatic nuclear development removes all transposons and other so-called "junk" DNA, which comprise ~95% of the germline. We demonstrate that Piwi-interacting small RNAs (piRNAs) from the maternal nucleus can specify genomic regions for retention in this process. Oxytricha piRNAs map primarily to the somatic genome, representing the ~5% of the germline that is retained. Furthermore, injection of synthetic piRNAs corresponding to normally deleted regions leads to their retention in later generations. Our findings highlight small RNAs as powerful transgenerational carriers of epigenetic information for genome programming.

Small RNA pyrosequencing in the protozoan parasite Entamoeba histolytica reveals strain-specific small RNAs that target virulence genes

BACKGROUND

Small RNA mediated gene silencing is a well-conserved regulatory pathway. In the parasite Entamoeba histolytica an endogenous RNAi pathway exists, however, the depth and diversity of the small RNA population remains unknown.

RESULTS

To characterize the small RNA population that associates with E. histolytica Argonaute-2 (EhAGO2-2), we immunoprecipitated small RNAs that associate with it and performed one full pyrosequencing run. Data analysis revealed new features of the 27nt small RNAs including the 5'-G predominance, distinct small RNA distribution patterns on protein coding genes, small RNAs mapping to both introns and exon-exon junctions, and small RNA targeted genes that are clustered particularly in sections of genome duplication. Characterization of genomic loci to which both sense and antisense small RNAs mapped showed that both sets of small RNAs have 5'-polyphosphate termini; strand-specific RT-PCR detected transcripts in both directions at these loci suggesting that both transcripts may serve as template for small RNA generation. In order to determine whether small RNA abundance patterns account for strain-specific gene expression profiles of E. histolytica virulent and non-virulent strains, we sequenced small RNAs from a non-virulent strain and found that small RNAs mapped to genes in a manner consistent with their regulation of strain-specific virulence genes.

CONCLUSIONS

We provided a full spectrum analysis for E. histolytica AGO2-2 associated 27nt small RNAs. Additionally, comparative analysis of small RNA populations from virulent and non-virulent amebic strains indicates that small RNA populations may regulate virulence genes.

mRNA decay proteins are targeted to poly(A)+ RNA and dsRNA-containing cytoplasmic foci that resemble P-bodies in Entamoeba histolytica

In higher eukaryotes, mRNA degradation and RNA-based gene silencing occur in cytoplasmic foci referred to as processing bodies (P-bodies). In protozoan parasites, the presence of P-bodies and their putative role in mRNA decay have yet to be comprehensively addressed. Identification of P-bodies might provide information on how mRNA degradation machineries evolved in lower eukaryotes. Here, we used immunofluorescence and confocal microscopy assays to investigate the cellular localization of mRNA degradation proteins in the human intestinal parasite Entamoeba histolytica and found evidence of the existence of P-bodies. Two mRNA decay factors, namely the EhXRN2 exoribonuclease and the EhDCP2 decapping enzyme, were localized in cytoplasmic foci in a pattern resembling P-body organization. Given that amoebic foci appear to be smaller and less rounded than those described in higher eukaryotes, we have named them "P-body-like structures". These foci contain additional mRNA degradation factors, including the EhCAF1 deadenylase and the EhAGO2-2 protein involved in RNA interference. Biochemical analysis revealed that EhCAF1 co-immunoprecipitated with EhXRN2 but not with EhDCP2 or EhAGO2-2, thus linking deadenylation to 5'-to-3' mRNA decay. The number of EhCAF1-containing foci significantly decreased after inhibition of transcription and translation with actinomycin D and cycloheximide, respectively. Furthermore, results of RNA-FISH assays showed that (i) EhCAF1 colocalized with poly(A)(+) RNA and (ii) during silencing of the Ehpc4 gene by RNA interference, EhAGO2-2 colocalized with small interfering RNAs in cytoplasmic foci. Our observation of decapping, deadenylation and RNA interference proteins within P-body-like foci suggests that these structures have been conserved after originating in the early evolution of eukaryotic lineages. To the best of our knowledge, this is the first study to report on the localization of mRNA decay proteins within P-body-like structures in E. histolytica. Our findings should open up opportunities for deciphering the mechanisms of mRNA degradation and RNA-based gene silencing in this deep-branching eukaryote.

Entamoeba histolytica: a snapshot of current research and methods for genetic analysis

Entamoeba histolytica represents one of the leading causes of parasitic death worldwide. Although identified as the causative agent of amebiasis since 1875, the molecular mechanisms by which the parasite causes disease are still not fully understood. Studying Entamoeba reveals insights into a eukaryotic cell that differs in many ways from better-studied model organisms. Thus, much can be learned from this protozoan parasite on evolution, cell biology, and RNA biology. In this review we discuss selected research highlights in Entamoeba research and focus on the development of molecular biological techniques to study this pathogen. We end by highlighting some of the many questions that remain to be answered in order to fully understand this important human pathogen.

Nucleus-localized antisense small RNAs with 5'-polyphosphate termini regulate long term transcriptional gene silencing in Entamoeba histolytica G3 strain

In the deep-branching eukaryotic parasite Entamoeba histolytica, transcriptional gene silencing (TGS) of the Amoebapore A gene (ap-a) in the G3 strain has been reported with subsequent development of this parasite strain for gene silencing. However, the mechanisms underlying this gene silencing approach are poorly understood. Here we report that antisense small RNAs (sRNAs) specific to the silenced ap-a gene can be identified in G3 parasites. Furthermore, when additional genes are silenced in the G3 strain, antisense sRNAs to the newly silenced genes can also be detected. Characterization of these sRNAs demonstrates that they are ~27 nucleotides in size, have 5'-polyphosphate termini, and persist even after removal of the silencing plasmid. Immunofluorescence analysis (IFA) and fluorescence in situ hybridization (FISH) show that both the Argonaute protein EhAGO2-2 and antisense sRNAs to the silenced genes are localized to the parasite nucleus. Furthermore, α-EhAGO2-2 immunoprecipitation confirmed the direct association of the antisense sRNAs with EhAGO2-2. Finally, chromatin immunoprecipitation (ChIP) assays demonstrate that the loci of the silenced genes are enriched for histone H3 and EhAGO2-2, indicating that both chromatin modification and the RNA-induced transcriptional silencing complex are involved in permanent gene silencing in G3 parasites. In conclusion, our data demonstrate that G3-based gene silencing in E. histolytica is mediated by an siRNA pathway, which utilizes antisense 5'-polyphosphate sRNAs. To our knowledge, this is the first study to show that 5'- polyphosphate antisense sRNAs can mediate TGS, and it is the first example of RNAi-mediated TGS in protozoan parasites.

A microRNA derived from an apparent canonical biogenesis pathway regulates variant surface protein gene expression in Giardia lamblia

We have previously shown that a snoRNA-derived microRNA, miR2, in Giardia lamblia potentially regulates the expression of 22 variant surface protein (VSP) genes. Here, we identified another miRNA, miR4, also capable of regulating the expression of several VSPs but derived from an unannotated open reading frame (ORF) rather than a snoRNA, suggesting a canonical miRNA biogenesis pathway in Giardia. miR4 represses expression of a reporter containing two miR4 antisense sequences at the 3' UTR without causing a corresponding decrease in the mRNA level. This repression requires the presence of the Giardia Argonaute protein (GlAgo) and is reversed by 2' O-methylated antisense oligo to miR4, suggesting an RNA-induced silencing complex (RISC)-mediated mechanism. Furthermore, in vivo and in vitro evidence suggested that the Giardia Dicer protein (GlDcr) is required for miR4 biogenesis. Coimmunoprecipitation of miR4 with GlAgo further verified miR4 as a miRNA. A total of 361 potential target sites for miR4 were bioinformatically identified in Giardia, out of which 69 (32.7%) were associated with VSP genes. miR4 reduces the expression of a reporter containing two copies of the target site from VSP (GL50803_36493) at the 3' UTR. Sixteen of the 69 VSP genes were further found to contain partially overlapping miR2 and miR4 targeting sites. Expression of a reporter carrying the two overlapping sites was inhibited by either miR2 or miR4, but the inhibition was neither synergistic nor additive, suggesting a complex mechanism of miRNA regulation of VSP expression and the presence of a rich miRNAome in Giardia.

Gene regulation in Giardia lambia involves a putative microRNA derived from a small nucleolar RNA

Two core microRNA (miRNA) pathway proteins, Dicer and Argonaute, are found in Giardia lamblia, a deeply branching parasitic protozoan. There are, however, no apparent homologues of Drosha or Exportin5 in the genome. Here, we report a 26 nucleotide (nt) RNA derived from a 106 nt Box C/D snoRNA, GlsR2. This small RNA, designated miR5, localizes to the 3' end of GlsR2 and has a 75 nt hairpin precursor. GlsR2 is processed by the Dicer from Giardia (GlDcr) and generated miR5. Immunoprecipitation of the Argonaute from Giardia (GlAgo) brought down miR5. When a Renilla Luciferase transcript with a 26 nt miR5 antisense sequence at the 3'-untranslated region (3' UTR) was introduced into Giardia trophozoites, Luciferase expression was reduced ∼25% when synthetic miR5 was also introduced. The Luciferase mRNA level remained, however, unchanged, suggesting translation repression by miR5. This inhibition was fully reversed by introducing also a 2'-O-methylated antisense inhibitor of miR5, suggesting that miR5 acts by interacting specifically with the antisense sequence in the mRNA. A partial antisense knock down of GlDcr or GlAgo in Giardia indicated that the former is needed for miR5 biogenesis whereas the latter is required for miR5-mediated translational repression. Potential targets for miR5 with canonical seed sequences were predicted bioinformatically near the stop codon of Giardia mRNAs. Four out of the 21 most likely targets were tested in the Luciferase reporter assay. miR5 was found to inhibit Luciferase expression (∼20%) of transcripts carrying these potential target sites, indicating that snoRNA-derived miRNA can regulate the expression of multiple genes in Giardia.

RNA interference in protozoan parasites: achievements and challenges

Protozoan parasites that profoundly affect mankind represent an exceptionally diverse group of organisms, including Plasmodium, Toxoplasma, Entamoeba, Giardia, trypanosomes, and Leishmania. Despite the overwhelming impact of these parasites, there remain many aspects to be discovered about mechanisms of pathogenesis and how these organisms survive in the host. Combined with the ever-increasing availability of sequenced genomes, RNA interference (RNAi), discovered a mere 13 years ago, has enormously facilitated the analysis of gene function, especially in organisms that are not amenable to classical genetic approaches. Here we review the current status of RNAi in studies of parasitic protozoa, with special emphasis on its use as a postgenomic tool.

RNA interference in Entamoeba histolytica: implications for parasite biology and gene silencing

Entamoeba histolytica is a major health threat to people in developing countries, where it causes invasive diarrhea and liver abscesses. The study of this important human pathogen has been hindered by a lack of tools for genetic manipulation. Recently, a number of genetic approaches based on variations of the RNAi method have been successfully developed and cloning of endogenous small-interfering RNAs from E. histolytica revealed an abundant population of small RNAs with an unusual 5´-polyphosphate structure. However, little is known about the implications of these findings to amebic biology or the mechanisms of gene silencing in this organism. In this article we review the literature relevant to RNAi in E. histolytica, discuss its implications for advances in gene silencing in this organism and outline potential future directions towards understanding the repertoire of RNAi and its impact on the biology of this deep-branching eukaryotic parasite.

Epigenetics in the unicellular parasite Entamoeba histolytica

Amoebiasis is a serious infectious disease that is caused by the unicellular parasite, Entamoeba histolytica. This parasite is mainly found in developing countries, and are named owing to its ability to destroy tissues. The molecular mechanisms that regulate the virulence of this parasite are not well understood. In recent years, an increasing interest in the epigenetic regulation of the parasite's virulence has emerged. In this article, an overview of our current knowledge about the role of DNA methylation, histone modifications and RNA-associated silencing in the biology of E. histolytica is provided. The relevance of some features of the parasite's unique epigenetic machinery to the development of new antiamoebic therapeutic molecules is discussed.

The emerging world of small silencing RNAs in protozoan parasites

A new RNA world has emerged in the past 10 years with the discovery of a plethora of 20- to 30-nucleotide long small RNAs that are involved in various gene silencing mechanisms. These small RNAs have considerably changed our view of the regulation of gene expression in eukaryotic organisms, with a major shift towards epigenetic and post-transcriptional mechanisms. In this article, we focus on the striking diversity of small silencing RNAs that have been identified in several protozoan parasites and their potential biological role.

Tissue destruction and invasion by Entamoeba histolytica

Entamoeba histolytica is the causative agent of amebiasis, a disease that is a major source of morbidity and mortality in the developing world. The potent cytotoxic activity of the parasite appears to underlie disease pathogenesis, although the mechanism is unknown. Recently, progress has been made in determining that the parasite activates apoptosis in target cells and some putative effectors have been identified. Recent studies have also begun to unravel the host genetic determinants that influence infection outcome. Thus, we are beginning to get a clearer picture of how this parasite manages to infect, invade and ultimately inflict devastating tissue destruction.

RNAi pathways in parasitic protists and worms

Tropical diseases caused by parasitic worms and protists are of major public health concern affecting millions of people worldwide. New therapeutic and diagnostic tools would be of great help in dealing with the public health and economic impact of these diseases. RNA interference (RNAi) pathways utilize small non-coding RNAs to regulate gene expression in a sequence-specific manner. In recent years, a wealth of data about the mechanisms and biological functions of RNAi pathways in distinct groups of eukaryotes has been described. Often, RNAi pathways have unique features that are restricted to groups of eukaryotes. The focus of this review will be on RNAi pathways in specific groups of parasitic eukaryotes that include Trypanosoma cruzi, Plasmodium and Schistosoma mansoni. These parasites are the causative agents of Chagas disease, Malaria, and Schistosomiasis, respectively, all of which are tropical diseases that would greatly benefit from the development of new diagnostic and therapeutic tools. In this context, we will describe specific features of RNAi pathways in each of these parasitic eukaryotic groups and discuss how they could be exploited for the treatment of tropical diseases.

Experimental design, preprocessing, normalization and differential expression analysis of small RNA sequencing experiments

Prior to the advent of new, deep sequencing methods, small RNA (sRNA) discovery was dependent on Sanger sequencing, which was time-consuming and limited knowledge to only the most abundant sRNA. The innovation of large-scale, next-generation sequencing has exponentially increased knowledge of the biology, diversity and abundance of sRNA populations. In this review, we discuss issues involved in the design of sRNA sequencing experiments, including choosing a sequencing platform, inherent biases that affect sRNA measurements and replication. We outline the steps involved in preprocessing sRNA sequencing data and review both the principles behind and the current options for normalization. Finally, we discuss differential expression analysis in the absence and presence of biological replicates. While our focus is on sRNA sequencing experiments, many of the principles discussed are applicable to the sequencing of other RNA populations.