Small RNA pyrosequencing in the protozoan parasite Entamoeba histolytica reveals strain-specific small RNAs that target virulence genes
© Zhang et al.; licensee BioMed Central Ltd. 2013
Received: 17 June 2012
Accepted: 2 January 2013
Published: 25 January 2013
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.
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.
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.
Keywords5′-Polyphosphate small RNA Argonaute High-throughput pyrosequencing Parasite
RNA interference (RNAi) is a well-conserved gene regulatory pathway found in most eukaryotes [1, 2]. Many important biological functions are controlled by RNAi such as developmental regulation , genome protection against viruses and transposons , and DNA elimination . Small RNA (sRNA) molecules, usually 20 – 30nt, are the key elements for RNAi. Guided by their associated protein complexes, they base-pair to the targeted transcripts or genomic loci to trigger gene silencing at either the transcriptional or post-transcriptional level . In recent years, high throughput sequencing has facilitated the identification of diverse species of small RNAs in different organisms .
Several protozoan parasites such as Trypanosoma brucei, Toxoplasma gondii, Giardia lamblia, Trichomonas vaginalis, and Entamoeba histolytica contain key genes of the RNAi pathway in their genomes [8, 9]. The functions of RNAi in parasite biology include retrotransposon control in T. brucei , gene regulation in E. histolytica [11, 12], and control of antigenic variation in G. lamblia [13, 14]. In T. gondii and T. vaginalis the studies have been limited to endogenous small RNA sequencing with no functional studies yet reported [15, 16].
E. histolytica causes dysentery and liver abscesses in humans and affects 500 million people worldwide . The study of this important human parasite has been hampered by lack of standard molecular genetic tools due to the polyploid nature of the E. histolytica genome . Recently, several RNAi-based gene knockdown approaches -- dsRNA/siRNA [19, 20], short-hairpin RNA  and a transcriptional gene silencing approach in the G3 parasite strain  -- have been established in this organism. We have shown that E. histolytica has a 27nt small RNA population, which has 5′-polyphosphate (5′-polyP) and 3′-OH termini and associates with EhAGO2-2 . Additionally, we have demonstrated that gene silencing in the E. histolytica G3 strain is mediated through a siRNA pathway .
The E. histolytica genome encodes three Argonaute proteins (EHI_125650, EHI_186850, and EHI_177170) of which EhAGO2-2 (EHI_125650) is highly expressed and associates with 27nt small RNAs . In this report, we immunoprecipitated small RNAs bound to EhAGO2-2 and sequenced them using a high throughput pyrosequencing approach, which generated over 360,000 small RNA reads. Analysis of these endogenous small RNAs revealed that their peak-length is 27nt and that there was a strong G bias at the 5′-nucleotide. Genome analysis showed that small RNAs largely mapped to annotated protein coding genes (overall mapping to ~4% of protein coding genes in the genome), which can be categorized into three groups: (I) genes with only antisense small RNAs; (II) genes with both antisense and sense small RNAs; and (III) genes with only sense small RNAs. Biochemical analysis revealed that both sense and antisense small RNAs (that map to group I or II genes) have 5′-polyP termini. Strand specific RT-PCR showed both sense and antisense transcripts can be detected for group II gene loci, suggesting that bi-directional transcripts are available to serve as templates for small RNA generation. Identification of small RNAs that map to introns as well as exon-exon junctions indicated that E. histolytica RNA-dependent RNA polymerase (RdRP) could use both nascent and mature transcript as a template for generating small RNAs. Whole genome microarray data indicated that most protein coding genes with large numbers of antisense small RNAs are not expressed, indicating that antisense small RNAs likely play a role in gene silencing. We further sequenced a size-fractionated small RNA library from the non-virulent E. histolytica Rahman strain and identified strain-specific patterns of antisense small RNAs mapping to a virulence gene in a manner consistent with its regulation by small RNAs. In summary, we have characterized the endogenous small RNA populations from a virulent and non-virulent strain of E. histolytica using a high-throughput small RNA pyrosequencing approach. Our analysis indicates that the small RNA pathway likely regulates expression of strain-specific genes in this parasite system, including some virulence determinants.
High-throughput pyrophosphate sequencing of EhAGO2-2-bound small RNAs
The algorithm for small RNA processing for the EhAGO2-2 IP small RNA library from E. histolytica HM-1:IMSS
5′-P independent cloning (5′CIP+PNK method)
HM-1:IMSS AGO2-2 IP library
Total reads by Pyrosequencing
Primer trimming, size-limiting (15–40nt)
Unique sequences (percentage)
209,513 unique sequences (158,904 cloned only once; 75.8%) (50,609 cloned more than once; 24.2%)
Scan for tRNA, rRNA
tRNA: 1,640 sequences rRNA: 3,456 sequences
Scan for SINE/LINE, EhERE elements
Map to E. histolytica genome
Map to E. histolytica predicted ORFs
Small RNAs are largely 27nt size with a striking 5′-G bias
We analyzed the size distribution and nucleotide composition of small RNAs and found that all reads (excluding reads from tRNAs, rRNAs and repetitive elements) peaked sharply at 27nt (Figure 1B), consistent with the size of small RNAs previously noted to associate with EhAGO2-2 . When nucleotide frequency was plotted at each position for the 27nt population we noticed a striking 5′-G bias (Figure 1C). The nucleotide frequency for the 26nt and 28nt small RNAs also shows a 5′-G propensity, but this was not the case for 17nt small RNAs indicating that the smaller sequences are likely degradation products (Additional file 1: Figure S3). Given that the E. histolytica genome is very AT-rich (approximately 75% AT) , the 5′-G bias in small RNAs is remarkable when compared to all the remaining plotted positions for small RNAs. A 5′-G bias in small RNAs has been reported in two other organisms, i.e. C. elegans and Ascaris, where 22G populations with 5′-polyP termini are defined as secondary siRNAs and thought to be generated by RdRP [29, 30]. Thus, we reason that the 5′-G biased 27nt small RNAs are likely RdRP-related and function as silencing siRNAs in E. histolytica.
Small RNA distribution patterns in the genome
The E. histolytica genome was first published in 2005  and a second version including new assemblies and reannotation was released in 2010 . The current genome assembly is still in scaffold stage, which contains 1,496 supercontigs and 8,201 genes. To gain an overview of small RNA distribution in the genome, we mapped small RNAs using the Bowtie alignment tool . We identified 216 supercontigs that had ≥100 small RNAs, and which accounted for 52% of the genome size (10.8 Mb/20.8 Mb) and which contained 96% of the total small RNAs. Specifically, the small RNAs were highly enriched on 19 supercontigs that although only 0.4 Mb in total size contained ~50% of all sequenced small RNAs (Additional file 2: Table S1). As the structure of the genome is unknown at present, we do not know if the 19 supercontigs that were enriched in small RNAs belong to centromeric or telomeric regions. To analyze the overall small RNA density distribution on the genome, we scanned the genome using a 500 bp window, and counted numbers of small RNAs in each window. We defined a “hot spot” as containing ≥100 small RNAs per window; using these parameters, there were 784 “hot spot” windows out of 41,600 genomic windows. The graphic views of “hot spot” distribution further revealed that most small RNAs arose either from some large clusters (≥ 3 genes and up to 14 genes that have abundant small RNAs) (Additional file 1: Figure S4A) or from isolated peaks (Additional file 1: Figure S4B). There are regions of the genome that had only a few mapped small RNAs (Additional file 1: Figure S4C).
We analyzed the protein coding genes to which small RNAs mapped, and found that many genes had only a few small RNAs that mapped to them and hence are likely artifacts. Thus, we tested our dataset with different cutoffs for the number of small RNAs mapping to a gene; we used four measures (no cutoff, ≥10, ≥25 or ≥50 small RNAs mapping to a gene). For each cutoff, we identified the number of protein coding genes in each category (antisense only small RNAs; sense and antisense small RNAs; and sense only small RNAs) and plotted the microarray expression value for these genes (see later section on "A global assessment of genes potentially regulated by small RNAs in E. histolytica"). For the two least stringent criteria (no cutoff or ≥10 small RNAs mapping to a gene), we observed no significant difference in the microarray expression value for the three categories of protein coding genes (Additional file 1: Figure S5A and Additional file 1: Figure S5B). However, when we used either the ≥25 or ≥50 small RNA cutoff, we identified significantly lower expression among genes with antisense or sense/antisense small RNAs compared to genes with sense small RNAs (Additional file 1: Figure S5C and Additional file 1: Figure S5D). The total number of protein coding genes using either cutoff (≥25 and ≥50 small RNAs) was relatively similar (420 and 358 genes, respectively). To be as stringent as possible, we decided to use a cutoff of ≥50 small RNAs mapping to a gene for further analysis. Overall, 358 protein coding genes (~4% of the genome) had ≥50 small RNAs that mapped to them. These protein coding genes could be categorized into three groups: (I) 226 genes with only antisense small RNAs; (II) 45 genes with both antisense and sense small RNAs; and (III) 87 genes with only sense small RNAs. Most genes in group I and II are annotated as hypothetical proteins (206 out of 271). However, a few gene families were represented including AIG1 family proteins (28 genes), beta-amylase (8 genes), deoxyuridine 5′-triphosphate nucleotidohydrolase domain proteins (6 genes), DNA polymerase (5 genes), and C2 domain proteins (2 genes).
In order to determine whether protein coding genes with small RNAs are in proximity to each other, we characterized the patterns of genes to which small RNAs map. A cluster is defined as ≥3 contiguous genes (with ≥50 small RNAs mapping to each gene). A pair is defined as 2 contiguous genes (with ≥50 small RNAs mapping to each gene) that are ≤1000 bp apart. There are a total of 358 protein coding genes that have ≥50 small RNAs and of these the majority are in clusters (167 genes) or in pairs (46 genes). These clustered/paired genes were largely in group I and II categories. We next looked at transcript orientation in the paired genes. The 46 protein coding genes in pairs had 6 divergent pairs (defined as tail to tail), 7 convergent pairs (defined as head to head), and 10 tandem pairs (defined as two adjacent genes in the same orientation). A total of 145 protein coding genes are unpaired with 63 genes in group I, 8 genes in group II, and 74 genes in group III.
When we characterized clustered protein coding genes (in group of ≥3) with small RNAs, we identified 26 clusters, which ranged in size from a cluster of 3 genes to clusters as large as 14 genes (Additional file 2: Table S2). A substantial number (19 out of 26) of clusters are in previously identified regions of D1, D2 and D4 genome duplication, which are segmental duplications with each segment flanked on both ends by inverted repeats such as IR/EhERE1/EhLINEs . Clusters that are not in regions of D1-D4 genome duplications are still flanked by repetitive elements at either one or both ends (Additional file 2: Table S2). Thus, clustered genes (that have large numbers of small RNAs) are more likely to be associated with repetitive elements.
In order to determine whether small RNAs are enriched on paired or clustered protein coding genes or in the intergenic DNA regions, we calculated the small RNA density on these paired/clustered genes as well as the intergenic regions between genes. This was calculated as small RNA/bp. We identified that the density of small RNAs mapping to intergenic regions was significantly lower compared to small RNA density mapping to paired/clustered genes (Additional file 1: Figure S6)(mean value 0.12 vs 0.54; p-value < 2.2e-16). This indicates that small RNA synthesis is most likely templated using a given gene rather than a long template covering several genes. For intergenic regions that had high small RNA density, we found these small RNAs are often in discrete sections or adjacent to predicted genes. Thus, we postulate that this may be due to small RNAs mapping to an unannotated gene or UTRs.
In summary, our analysis suggests that the small RNA targeted protein coding genes tend to be in pairs or clusters, and that clusters of genes with small RNAs are more often associated with repetitive elements. Small RNA density on paired/clustered genes versus intergenic regions implies that it is unlikely for either DNA or a long transcript covering several genes to be used as a template, but rather that transcript derived from each gene is the most likely template.
Small RNA distribution patterns within protein coding genes
In contrast, when small RNA distribution was plotted for group III genes (87 genes with only sense small RNAs), it became apparent that sense small RNAs were enriched towards the 3′-end of each gene (81/87 or 93%) (Figure 2C). We noticed that genes with only sense small RNAs mainly code for four highly expressed gene families (i.e. 40S ribosomal protein S16 (4 genes); Gal/GalNAc lectin (6 genes); three protein kinase families (44 genes); and hypothetical proteins (4 genes)). Since our cloning method will capture all types of small RNAs (5′-cap, 5′-PPP, 5′-P, 5′-OH), degradation products with 5′-OH species would also be cloned; this may suggest that at least some of the sense small RNAs might be mRNA degradation products. However, we found no sense small RNAs to other highly expressed genes, indicating some specificity for sense small RNAs to these specific loci. Since the mapping of sense small RNAs to many genes often had abrupt boundaries within/near 3′ end of genes and sometimes extended into intergenic regions (Additional file 1: Figure S7), we felt that some mapping artifacts could be due to poor genome quality in these regions. As the E. histolytica genome is not complete at present, we cannot rule out whether or not these sense small RNAs are derived from some other loci.
Although sense small RNAs have been identified in almost all cloning libraries, these types of small RNAs have largely been ignored since it is hard to evaluate their identity as true endogenous small RNAs versus non-specific hydrolysis fragments. In both C. elegans and Ascaris suum, where 5′-polyP small RNAs were found, sense small RNAs were identified but were not characterized [30, 33]. For the E. histolytica group III genes that have only sense small RNAs mapped to them we have not made any assumptions and have not further characterized the structure, derivation, or function of this category of sense small RNAs.
Antisense and sense small RNAs (for group I and II genes) have 5′-polyphosphate termini
The secondary 5′-polyP small RNAs in nematodes (both C. elegans and Ascaris suum) are thought to be generated by an amplified gene silencing mechanism, likely through RdRP [30, 33]. We have previously demonstrated that in E. histolytica both endogenous and exogenously triggered antisense small RNAs have 5′-polyP termini [11, 12]. Since our analysis to date had focused on the antisense small RNAs, we wanted to determine whether the small RNAs that map sense to genes (in the group II category) also have 5′-polyP termini. For this purpose, we analyzed small RNAs that mapped sense to the EHI_020920 locus. This gene is highly enriched for both antisense and sense small RNAs with small RNA coverage extending into the apparent 5′-upstream region (Figure 3A). Northern blot analysis detected signal at 27nt with probes that detect antisense small RNAs (probes 1, 2 and 3) and with probes that detect sense small RNAs (probes 4 and 5) (Figure 3B). All the blots were performed using the same membrane, and therefore the relative band intensity reflects the abundance of small RNAs. We observed a correlation of 5'-enriched distribution of small RNAs with the intensity of probe 1 > probe 2 > probe 3. We then performed a Terminator exonuclease assay and a 5′-end Capping assay on the total RNA sample. The signal for sense small RNAs (probes 4 and 5) were resistant to Terminator treatment and shifted up after capping assay, indicating that these small RNAs that map sense to the EHI_020920 locus have 5′-polyP termini (Figure 3C). Probe 1, which would detect antisense small RNAs had the expected biochemical features consistent with 5′-polyP termini. The control (a “spiked” synthetic 21-mer RNA species, which was labeled to have a 5′-monophosphate structure) was, as expected, degraded by Terminator enzyme and unaffected by treatment with capping enzyme.
We further examined a second example: the locus (DS571502:1400–2800 bp) that contains both antisense and sense small RNAs to EHI_130480 and EHI_130490 as well as one potential unannotated gene (Additional file 1: Figure S8). Probes for detecting antisense (probe 7) and sense (probe 6) small RNAs were chosen for EHI_130480. Northern blot analysis showed signal for both probes at ~27nt; the sense small RNA was resistant to Terminator exonuclease and thus has 5′-polyphosphate termini. Thus, for the two loci tested (both of which had abundant antisense and sense small RNAs) we determined that the sense small RNAs also have 5′-polyP termini.
The biochemical features of sense small RNAs (in group II genes) having 5′-polyP termini was unexpected. Analysis for pairing between antisense and sense small RNAs showed no enriched pairs for these gene loci (data not shown). In a typical RNAi pathway, dsRNA is chopped into a siRNA duplex, where enriched pairing between antisense and sense small RNAs can be found [34, 35]. The lack of pairing between antisense and sense small RNAs in group II genes and the 5′-polyP termini for both antisense and sense small RNAs indicates that these small RNAs could be individually processed from bidirectional transcripts at these loci. To determine if this may have occurred, we performed strand-specific RT-PCR for EHI_020920, EHI_130480 and EHI_130490. At all three loci, both sense and antisense transcripts can be detected, albeit with antisense transcripts at much lower level than the sense transcript (Figure 3D, Additional file 1: Figure S8C). Overall transcript levels as assayed by RT-PCR correlated with the abundance of antisense and sense small RNAs mapped to these loci. We further applied the strand-specific RT-PCR assay to two additional protein coding genes (EHI_192590 and EHI_148800), which had both antisense and sense small RNAs and identified transcripts in both directions for these loci (data not shown). Thus, we conclude that although bidirectional transcription has not been previously documented in E. histolytica, both antisense and sense small RNAs are likely templated from bidirectional transcripts and generated by an RdRP-dependent mechanism resulting in small RNAs with 5′-polyphosphate termini. Natural antisense transcripts have been shown to be a major source of siRNA generation in Drosophila melanogaster  and in plants . In the parasite, Giardia lamblia, bidirectional transcription produces abundant sterile antisense transcripts . Further studies to characterize the extent of natural antisense transcripts in E. histolytica and the mechanism for generation of 5′-polyP small RNAs are needed. Our demonstration that small RNAs to group II genes (genes with both antisense and sense small RNAs) have 5′-polyP termini and strand-specific RT-PCR detecting transcripts in both directions indicates that these sense small RNAs are likely generated from antisense transcript, and indicates that some portion of “sense” small RNAs in our dataset although called “sense”, are truly antisense to the transcript derived from the opposite strand of the examined gene.
Small RNAs are derived from both unspliced and spliced transcripts
In order to identify whether small RNAs map to spliced or unspliced transcripts, we downloaded both genomic and mRNA sequences for all E. histolytica protein coding genes with at least one predicted intron. Small RNAs mapping to introns are defined as those that map to the genomic gene sequence but not the mRNA sequence. Small RNAs mapping to exon-exon junctions are those that map to mRNA sequence but not the genomic gene sequence. We found a total of 52 small RNA reads that spanned exon-exon junctions (51 mapped in the antisense orientation and 1 mapped in the sense orientation) and 1,187 small RNA reads that mapped to predicted introns (908 mapped in the antisense orientation and 279 mapped in the sense orientation) (Additional file 2: Table S3).
We further examined the protein coding genes with at least 50 small RNAs and which also have at least one intron (Additional file 3: Table S8). For this list, we first checked for potential “false introns” caused by genome sequence error, and excluded EHI_018150, EHI_180820 and EHI_137120 as they have many Ns in their intron sequences. We then checked the predicted intronic sequence in the remaining genes for an in-frame stop codon or frame disruption, as this strongly suggests that the intron is correctly predicted (Additional file 3: Table S8). Lastly, we checked for paralogs within these genes. With these criteria, we examined the three categories of genes to which small RNAs mapped to identify small RNAs that mapped to exons, introns and exon-exon junctions (Additional file 3: Table S8). For the group I genes, 18 unique genes had small RNAs that mapped to predicted introns and among them, 4 unique genes also had small RNAs that mapped to exon-exon junctions. For the group II genes, 4 unique genes had antisense small RNAs that mapped to introns and 3 of these also had sense small RNAs that mapped to introns. None of the genes in this category had small RNAs that mapped to exon-exon junctions. For the group III genes, 4 unique genes had small RNAs mapped to predicted introns, but none mapped to exon-exon junctions. Overall, we made a number of observations: (i) a greater number of small RNAs mapped to exons than to the introns for all three groups of genes; (ii) all but one intron-containing gene in groups I and II had small RNAs that mapped to introns; and (iii) only a limited number of genes (4 unique genes) had small RNAs that mapped to exon-exon junctions. These data indicate that both spliced and unspliced transcripts are capable of being used as templates to produce small RNAs in E. histolytica. As an example, the mapping of small RNAs to exons, introns and exon-exon junctions are shown for EHI_135940 and EHI_197360 genes (Additional file 1: Figure S9). Further calculations of the small RNA density revealed four-fold greater density of small RNAs in exons than in introns [EHI_135940, exon (0.77) vs. intron (0.19); EHI_197360, exon (0.39) vs. intron (0.10)] (Additional file 2: Table S4). The difference could suggest that spliced transcripts are preferred as templates to unspliced transcripts, or alternatively may simply be a reflection of the ratio of spliced and unspliced transcripts available in the cell.
In C. elegans, EGO-1, an RdRP, is critical for C. elegans germline development and is responsible for producing 5′-polyP antisense small RNAs from mRNA-derived loci. Small RNA sequencing has shown that small RNAs often span exon-exon junctions and rarely map to introns, indicating EGO-1 uses processed mRNA as a template . RdRP could theoretically template on genomic DNA, nascent transcripts, or processed mRNAs. Small RNAs that map to exon-exon junctions provide evidence that a spliced mRNA template is used to generate these small RNAs, whereas small RNAs that map to introns indicate that non-spliced templates can also be used to generate small RNAs. Based on the observations that there are small RNA free genomic regions between genes with antisense small RNAs, and that many more small RNAs map to exons than introns, we conclude that the E. histolytica RNAi machinery prefers mature transcript as a template for generating small RNAs. However, the machinery in E. histolytica also seems capable of using unspliced transcripts as template, although at reduced levels. Whether this is indicative of the inherent preference of the E. histolytica machinery or due instead to the low abundance of unspliced mRNA is not clear at present.
Small RNAs that map to tRNAs, rRNAs and retrotransposon elements
In order to identify small RNAs that map to the tRNAs, rRNAs and retrotransposon elements, we followed the outline in Additional file 1: Figure S1. E. histolytica has uniquely organized tRNA genes that are in multiple tandem-array units, likely arranged at subtelomeric regions and spaced by tandem repeats of AT-rich sequences . The E. histolytica rRNA genes reside on an extrachromosomal circular plasmid and two rRNA transcription units are organized as inverted repeats . We mapped the small RNA reads to the tRNA repeat units and the rRNA plasmid. We found that nearly all small RNA reads were in the sense orientation to the coded tRNAs and rRNAs. Additionally, we plotted the size distribution and nucleotide frequency for rRNA reads and noted that small RNAs in these categories did not peak at 27nt (Additional file 1: Figure S10), and did not have a 5′-G enrichment (Additional file 1: Figure S11) indicating that these small RNAs are most likely degradation products resulting from these highly expressed structural RNAs. However, we are aware of some recent reports showing that small RNAs could originate from tRNAs and snoRNAs [42, 43]. Upon a closer examination of the tRNA reads, we noticed a slight peak at 27nt with some degree of 5′-G enrichment for this 27nt population only (Additional file 1: Figure S10 and Additional file 1: Figure S11). Whether or not this indicates that these are functional small RNAs in Entamoeba needs further study.
Transposons and repetitive DNA are abundant in Entamoeba and hundreds of copies of the long interspersed nuclear elements (LINEs) and short interspersed nuclear elements (SINEs) can be found in the E. histolytica genome [44, 45]. Our small RNA dataset contains 5% of reads that mapped to LINE and SINE elements. Analysis of the lengths of these small RNAs showed two peaks (one at 27nt and the other at 17nt) (Additional file 1: Figure S10). When nucleotide composition is plotted, the 5′-G propensity is apparent for the 27nt peak, but not for the 17nt peak (Additional file 1: Figure S12). Thus, the 27nt small RNA population that maps to LINE/SINE elements had features similar to those that map to coding regions and are likely not artifacts. As an example, we mapped small RNAs to the EhRLE5 sequence, which has been categorized in the EhLINE1 family . The small RNAs are scattered along the whole region on both strands and cover the whole EhRLE unit, with a slight increase in small RNAs near each end (Additional file 1: Figure S13). We got a positive signal from Northern blot analysis using several probes to retrotransposon elements although the size by Northern blot analysis was slightly higher (~32nt) than the cloned small RNA. This indicates that small RNAs could derive from these retrotransposon elements.
The mapping of small RNAs to D1-D4 repetitive segments showed a large number of small RNA reads on D1, D2 and D4 segments but not on D3 segments (Additional file 2: Table S2). Annotated protein coding genes in these duplication regions appear to be covered by large numbers of antisense small RNAs, forming a large cluster (Additional file 1: Figure S14 and Additional file 2: Table S2). Thus, although the overall numbers of small RNAs that associate with EhAGO2-2 and map to repeat or retrotransposon elements is low, they may play a functional role in controlling genome stability as has been shown in other systems [47, 48]. Alternatively, small RNAs may have a role in controlling these retrotransposon elements but may do so by associating with the two other Argonaute proteins in E. histolytica.
A global assessment of genes potentially regulated by small RNAs in E. histolytica
For the genes in group III (≥50 mapped sense small RNAs), 87 genes met our criteria with 64 represented on the microarray. Plotting the number of mapped reads for each gene as a function of its normalized microarray gene expression value showed no direct link between numbers of sense small RNAs to a gene and the expression level of that gene (Figures 4C and 4D). We reasoned that for genes with high expression value, the sense small RNAs may represent degradation products. However, for genes with very low mRNA expression (15 of the 87 on the microarray), we have no good explanation on how these sense small RNAs were generated. Whether this represents an artifact of errors in genome assembly or some other factor is not clear at present. Interestingly, none of these genes are regulated under other conditions tested (data not shown).
In summary, the strongest correlation between gene expression and small RNA abundance was for genes that had either abundant antisense small RNAs or abundant sense and antisense small RNAs. Although simply correlative at present, the inverse correlation between antisense small RNA abundance and gene expression suggests that antisense small RNAs mediate target gene silencing in E. histolytica and are potentially involved in gene regulation under various conditions.
Small RNA sequencing from the nonvirulent E. histolytica Rahman strain
It has been observed that pathogenicity varies greatly among different E. histolytica strains. E. histolytica HM-1:IMSS is a virulent strain, while E. histolytica Rahman is a nonvirulent strain . We have previously demonstrated that gene expression profiles are substantially different between these two strains and have used the distinct strain-specific expression profiles to identify virulence genes . Given the data above, and to explore whether small RNAs play a role in strain-specific and/or virulence gene regulation, we constructed a small RNA library from trophozoites of the E. histolytica Rahman strain.
We visualized the small RNA populations in E. histolytica Rahman by separating total RNA on a 12% denaturing polyacrylamide gel followed by Sybr gold staining and visualized an abundant 27nt small RNA population (data not shown). A small RNA library was constructed from size-fractionated RNA (15-30nt) using a 5′-P independent cloning approach and limited pyrophosphate sequencing was performed generating 151,656 reads. For the purpose of mapping, we used E. histolytica HM-1:IMSS genome as a reference genome rather than the current E. histolytica Rahman assembly, based on the following facts: (a) the current E. histolytica Rahman genome assembly (Amoeba DB, http://amoebadb.org) is in a preliminary stage containing 17,378 small contigs and is unannotated, (b) there is a high similarity between these two strains and one previous study has estimated that only 5 out of a sample of 1,817 (0.3%) genes were identified as highly or significantly divergent , and (c) Affymetrix platform microarrays found no difference in overall hybridization efficiency levels compared to HM-1:IMSS (i.e. the signal scale factor was similar between two arrays and thus both RNA samples hybridized efficiently to the array chip), indicating a high level of sequence identity for the protein coding genes . We realize that sequence differences between the E. histolytica HM-1:IMSS and Rahman strains may cause us to lose some data. However, the advantages of being able to map to an annotated genome and thus determine how many small RNAs map to protein coding genes (and whether they map sense or antisense) and to intergenic regions were significant enough that we proceeded with the data generated by aligning the E. histolytica Rahman small RNA library to the E. histolytica HM-1:IMSS genome sequence.
The algorithm for small RNA processing for the E. histolytica Rahman size-fractionated small RNA library
5′-P independent cloning (5′CIP+PNK method)
Rahman total RNA, size fractionation (15-30nt) library
Total reads by Pyrosequencing
Primer trimming, size-limiting (15–40nt)
Unique sequences (percentage)
98,414 unique sequences (82,780 cloned only once; 84.1%) (15,634 cloned more than once; 15.9%)
Scan for tRNA, rRNA
tRNA: 517 sequences rRNA:1,753 sequences
Scan for SINE/LINE, EhERE elements
Map to E. histolytica genome
Map to E. histolytica predicted ORFs
To be confident that our mapping using the HM-1:IMSS genome represented the overall picture for the Rahman strain, we also aligned these small RNA reads to the current Rahman assembly and compared the number of mapped reads. Aligning the Rahman small RNA dataset to the E. histolytica Rahman assembly rather than the HM-1:IMSS genome only increased the number of mapped reads by 7% (7,232 reads). Overall, this indicates that mapping of small RNAs from Rahman to the HM-1:IMSS genome is representative of the overall picture and the greater information gained by mapping the small RNAs to an annotated assembly outweighed the negative effects of using a genome sequence from a different strain.
Genes with strain-specific expression patterns and roles in virulence have small RNAs that map to them
We cannot directly compare the HM-1:IMSS and Rahman libraries for read frequencies and patterns of small RNA mapping, as they were generated differently and were sequenced to different depths. However, common features in both libraries (peak at 27nt, 5′-G bias, and mapping antisense to genes) suggested that EhAGO2-2 bound species are represented in the Rahman library. Thus, this overlap of small RNA coverage between the two libraries allowed us to look for genes, which might be regulated strain specifically on the basis of small RNA abundance. An important caveat to the analysis is that while the presence of small RNAs is meaningful in either strain, the absence of small RNAs in the Rahman library is less meaningful and could be due to the limited sequencing depth or use of a size-selected small RNA library.
Using the same criteria as used for HM-1:IMSS (genes with ≥50 small RNAs that map to them) the small RNA library generated from E. histolytica Rahman identified 223 genes with small RNAs (175 genes with antisense small RNAs; 29 genes with both antisense and sense small RNAs; and 19 genes with only sense small RNAs). When the genes from these three categories were compared between the Rahman and HM-1:IMSS strains, we found significant overlap for genes with antisense and sense/antisense small RNAs. A total of 90 genes with only antisense small RNAs were common between the two strains (out of a total of 226 genes for HM-1:IMSS and 175 genes for Rahman); 16 genes with both antisense and sense small RNAs were common between the two strains (out of a total of 45 genes for HM-1:IMSS and 29 genes for Rahman). However, for the genes with sense only small RNAs, no overlap was identified (out of a total of 87 genes for HM-1:IMSS and 19 genes for Rahman). These data further support the idea that antisense small RNAs are likely playing roles in conferring strain-specific gene expression profiles, whereas the small RNAs that map sense only to genes are likely to be random degradation products.
In order to compare gene expression patterns for genes with strain-specific small RNAs, we used previously published microarray data for the two E. histolytica strains and a two-fold difference in normalized expression and a p-value <0.05 as cutoff . Of the 160 genes that have lower expression in HM-1:IMSS than in the Rahman strain, 11 genes could be potentially regulated by antisense or antisense/sense small RNAs in HM-1:IMSS (these genes had small RNAs in the HM-1:IMSS strain but none in the Rahman strain). Of the 130 genes that have higher expression in the HM-1:IMSS strain compared to the Rahman strain, 19 genes could be potentially regulated by antisense or antisense/sense small RNAs in Rahman (these genes had small RNAs in the Rahman strain, but none in the HM-1:IMSS strain).
The E. histolytica serine, threonine, isoleucine, rich protein (EhSTIRP) was previously identified as a virulence determinant as it has high expression in E. histolytica virulent strains and no expression in nonvirulent strains and genetically proven to have a role in virulence . There were 519 small RNAs that mapped antisense to this gene from our Rahman small RNA dataset compared to virtually no antisense small RNAs from the HM-1:IMSS small RNA dataset (Figure 5D). The distribution of the antisense small RNAs to EhSTIRP gene showed a clear enrichment at the 5′-end. We further confirmed the presence of EhSTIRP antisense small RNAs in the Rahman strain (and absence in the HM-1:IMSS strain) by Northern blot analysis (Figure 5D) using total RNA samples from two strains. Lastly, we attempted to overexpress a Myc-tagged EhSTIRP1 (EHI_025700), which resulted in successful EhSTIRP overexpression in the HM-1:IMSS strain but not in the Rahman strain even at very high levels of drug selection (data not shown). The lack of ability to overexpress the EhSTIRP gene in the E. histolytica Rahman strain that has abundant antisense small RNAs that map to it raises the intriguing possibility that the antisense small RNAs control the silencing of this virulence gene in the E. histolytica Rahman strain. Future studies will be needed to address this question.
In this study, we performed pyrosequencing of E. histolytica endogenous small RNAs that associate with the EhAGO2-2 protein in the virulent E. histolytica HM-1:IMSS strain and size selected small RNAs in the non-virulent E. histolytica Rahman strain. Characterization of these small RNAs showed that they are distinct in size (27nt), have preference for 5′-G, and have 5′-polyphosphate termini. Genome mapping revealed that these 27G-RNAs are mainly derived from coding genes with a much smaller population coming from retrotransposons. Comparison of the 27nt small RNA population from E. histolytica strains with varying virulence and expression profiles demonstrated an inverse correlation between antisense small RNA and gene expression levels, hinting that antisense 27G-small RNAs may be involved in the regulation of strain-specific genes, including known virulence determinants, such as EhSTIRP.
We have previously demonstrated that 27nt antisense small RNA in E. histolytica have 5′-polyP termini . In this study, we further characterized small RNAs that associate with EhAGO2-2, including those small RNAs that map both antisense and sense to genes. For all the small RNAs that were detectable by Northern blot analysis, we were able to show that they have 5′-polyP termini, indicating these abundant small RNA species in E. histolytica are reminiscent of secondary small RNA from C. elegans and Ascaris [29, 30]. In worms, the biogenesis and stability of 5′-polyphosphate small RNAs are dependent on RdRPs and WAGOs respectively, with these small RNAs being a component of the siRNA amplification pathway . Signal amplification is controlled by using primary trigger siRNAs to instigate secondary siRNAs through RdRP for the enforced silencing, but limiting secondary siRNAs from doing further signal amplification . Although small RNAs that mapped sense to coding regions were found in C. elegans and Ascaris suum 5'-monoP independent libraries, their existence and functionality were not confirmed; instead they were generally treated as non-specific degradation products [30, 33]. We have confirmed that in E. histolytica small RNAs that map sense to genes (which also have antisense small RNAs) are detected in Northern blot analyses as discrete band and bear the same 5′-polyP termini. Strand-specific RT-PCRs detected transcripts in both directions for these loci (albeit at a lower level for the antisense transcript), implying that RdRP-based small RNA generation could occur for both sense and antisense transcripts in this parasite. The E. histolytica genome encodes one full RdRP gene (EHI_139420) and two genes with partial RdRP domains (EHI_179800 and EHI_086260). The functions of the RdRP genes in E. histolytica RNAi pathway are still elusive and need further investigation.
Due to the fact that gene knockout is not feasible in E. histolytica we have been unable to dissect how the parasite RNAi components affect the levels of these small RNAs. The comparison of 5′-polyP small RNAs among the three organisms in which they have been described (C. elegans, Ascaris and E. histolytica) show several differences: (a) 5′-polyP small RNA size in E. histolytica is 27nt, whereas in C. elegans and Ascaris these RNAs are 22nt; (b) the distribution pattern of antisense small RNAs to the targeted gene loci is enriched at the 5′-end for E. histolytica and Ascaris whereas in C. elegans there is enrichment at the 3′-end of transcripts; (c) localization of EhAGO2-2 and bound 27nt small RNAs are mostly localized to the parasite nucleus , whereas C. elegans 22G-RNAs can associate with several different WAGOs and have both perinuclear  and nuclear localization [59–61]; and (d) C. elegans strongly prefers spliced transcripts as RdRP template for generating 5′-polyP small RNAs whereas both mature and nascent transcripts appear to function as templates in E. histolytica. Future studies aimed at elucidating these different mechanisms are needed.
Our previous limited Sanger sequencing has shown that small RNAs in E. histolytica largely mapped to the coding genes . Our pyrosequencing data further confirmed this mapping, which is in contrast to other parasitic systems T. gondii, G. intestinalis and T. brucei where the small RNAs are 5′ single phosphate and mostly derived from repetitive elements, retrotransposons [10, 15, 16]. The genome of E. histolytica contains hundreds of copies of LINE and SINE elements, with SINE elements actively transcribed and LINE1 transcript detected by Northern blot analysis [55, 62]. How retrotransposons are controlled in ameba is not known. As noted in our small RNA dataset analysis, one way that the RNAi pathway could function in this parasite is to target the unwanted transposon insertion events where segments are flanked on both ends by inverted repeats (IR/EhERE1/EhLINEs) as many small RNAs map to these regions. Whether RNAi directly targets retrotransposons in ameba is still an open question. Further characterization of the other two Argonaute proteins (EhAGO2-1 and EhAGO2-3) will provide a complete picture of small RNA populations in ameba and their functions in retrotransposon silencing.
The small RNA sequencing from different strains (HM-1:IMSS and Rahman) clearly indicated that expression of a subset of genes, including the virulence factor EhSTIRP, appears to be controlled by small RNAs in a strain-specific manner. In order to determine direct effects of the small RNA repertoire on parasite pathogenesis, the next step will be to perform functional studies to demonstrate direct roles for these small RNAs in regulating strain-specific virulence gene expression.
In summary, we present two pyrosequencing small RNA datasets from the parasite E. histolytica: one an EhAGO2-2 IP library from the virulent E. histolytica HM-1:IMSS strain and the other a size selected small RNA library from the non-virulent E. histolytica Rahman strain. Our analysis identified a number of new findings: (1) amebic 27nt small RNAs have 5′ G preference; (2) antisense small RNA targeted genes are in pairs or clustered and notably most of clusters are from segmental duplications D1, D2 and D4; (3) characterization of group II gene loci shows that both sense and antisense small RNAs have 5′-polyphosphate termini; (4) small RNAs mapping to introns and exon-exon junctions were found indicating that both spliced and unspliced mRNA can serve as the templates for small RNA production; (5) few small RNAs are found in intergenic regions between paired/clustered genes indicating that RNA transcript from each gene was used as template; (6) small RNAs targeting retrotransponsons have similar features to the small RNA targeting the mRNAs, but are not highly abundant; and (7) antisense small RNAs may contribute to differential gene expression between virulent and nonvirulent amebic strains including the known virulence gene EhSTIRP. Thus, the two small RNA datasets in this study will provide important data for the community to study small RNA-mediated gene regulation in this important human pathogen.
Parasite culture and RNA preparation
Entamoeba histolytica trophozoites (HM-1:IMSS and Rahman) were grown under standard conditions as previously published . A transfectant cell line expressing N-terminal Myc tagged EhAGO2-2 was previously described  and was maintained at 24μg/ml G418. For isolation of RNA that immunoprecipitated (IP) with EhAGO2-2, anti-Myc antibody (Pierce) was incubated with parasite lysate (2 hours, 4ºC), washed twice with 1x IP solution and pelleted. RNA was then isolated using mirVANA kit (Ambion). Small RNA enriched material from E. histolytica Rahman strain was prepared according to the mirVANA kit protocol.
Library construction and sequencing
Small RNA cloning was based on the previously published protocol of the 5′-phosphate independent cloning method (CIP + PNK method) . For the EhAGO2-2 IP RNA library, RNA was extracted from the IP material, directly ligated to the 3′-adapter oligonucleotide, size fractionated and treated with CIP and PNK, and then ligated to the 5´-adapter oligonucleotide. For the Rahman small RNA library, 100 μg of small RNA enriched RNA was size-fractionated on a 12% TBE urea polyacrylamide gel, the 15-30nt fraction was excised, followed by 3′-adapter ligation, CIP + PNK treatment, and 5′-adapter ligation as above.
For both libraries, the final ligated RNA with both 5′- and 3´-adapters was converted to single-stranded cDNA using Superscript II reverse transcriptase (Invitrogen, CA, USA). The cDNA was PCR-amplified using 454 Primers (A and B) for 20 cycles, and resolved on a 4% low melting point agarose gel. The band at 100 bp was excised and purified, and then further heat-denatured and purified from a 6% PAGE-urea gel. After a second round of purification, the recovered DNA was resuspended in Qiagen elution buffer, and pipelined into the 454 sequencing procedure.
Small RNA read processing and mapping
All sequencing reads were processed by first removing the linker from both ends; the resulting sequences were analyzed with Unix tools and unique sequences selected. The unique reads were mapped against the E. histolytica HM-1:IMSS genome, release 1.3, (http://amoebadb.org/amoeba/) using the program Bowtie (http://bowtie-bio.sourceforge.net/index.shtml) with parameters set as -v 1, -k 5. Mapped reads were visualized with the genome browser IGV (http://www.broadinstitute.org/igv/). To identify small RNAs that map to exon-exon junctions, TopHat (2.0.0 release 4/09/2012, http://tophat.cbcb.umd.edu) was used with the following parameters: -report-secondary-alignments -G Ehistolytica_AmoebaDB-1.3.gtf -i 20. For the scaffold view of mapped small RNAs, histograms were generated in R (http://cran.r-project.org/), using a window size of 500 bp to divide the scaffolds. For tRNA and rRNA analysis, we downloaded all tRNA array sequences from NCBI based on previously published analysis . For repetitive element analysis, we made a custom dataset using all SINEs/LINEs/EhEREs coordinates recently deposited to AmoebaDB (personal communication, Omar Harb), then aligned all sequences with Bowtie using the parameters -v 2, -k 5. All sequences from both small RNA libraries have been deposited at AmoebaDB (http://amoebadb.org/amoeba/). The small RNA sequence data from this study have been submitted to the NCBI Gene Expression Omnibus (GEO) under accession number GSE43668.
Nucleotide composition of small RNAs and distribution of small RNAs on gene loci
For each subset of small RNA populations the small RNA read sequences were extracted. The nucleotide composition at each position was counted using the R package ShortRead, (http://bioconductor.org/ packages/release/bioc/html/ShortRead.html) and frequency bar plots generated from these data. For small RNA distribution on gene loci, only genes with ≥ 50 small RNAs mapping were considered. Each small RNA sequence was assigned a position value based on the position of the starting nucleotide along the gene (normalized by gene length; each gene has a total length of 1). Histograms for these values were plotted in R.
Small RNA Northern blot analysis
Small RNA Northern blot analysis was done as previously published . Briefly, 20 μg - 50 μg of small RNA enriched samples were separated on a denaturing 12% polyacrylamide gel and transferred to a membrane. This was probed with end-labeled 32P-labeled oligonucleotides in perfectHyb buffer (Sigma) at 37°C and washed using low (2X SSC, 0.1% SDS at 37°C for 15 min) and medium (1X SSC, 0.1% SDS at 37°C for 15 min) stringency conditions. All probes used are listed in Additional file 2: Table S6.
Small RNA molecular analyses
For enzymatic analyses of RNA material, assays were done as previously published . Briefly, either 10 μg or 50 μg small RNA enriched RNA sample was spiked with a control sample (a synthetic 21-mer RNA with the 5′-end labeled with 32P). For the Terminator assay, the sample mixture was treated with Terminator enzyme (Epicentre), following the provided protocol from the manufacturer. For the capping assay, the ScriptCap m7G capping system (Epicentre) was used with the alternate cap zero capping protocol. After enzymatic treatment, samples were phenol:chloroform extracted and resolved on a 12% polyacrylamide gel. Northern blot analysis was performed using a radiolabeled probe to detect the small RNA of interest.
Strand specific RT-PCR analyses
We used SuperScript III first-strand synthesis kit (Invitrogen) for strand-specific cDNA synthesis and PCR analysis. E. histolytica HM-1:IMSS total RNA was treated with DNase I, and purified strand-specific primer for each gene (or oligo dT primer) was added to 0.5 μg total RNA reaction and heated to 65°C for 5 min. The temperature was lowered to 55°C, and prewarmed cDNA synthesis mix was added to the reaction and incubated at 55°C for 50 min. The reaction was terminated at 85°C for 5 min, chilled on ice, and 1 μl of RNase H was added to each tube and incubated for 20 min at 37°C before proceeding to PCR. In each primer reaction, both + RT and –RT reactions were performed, and the final cDNA volume was 20 μl. PCR was performed using 1 μl cDNA for a 30 μl PCR reaction for 33 cycles (94°C 15''; 55°C 30''; 72°C 1 min). Half of the volume of PCR reaction was loaded on a gel for visualization. Primers used are listed in Additional file 2: Table S7.
We thank all members of the Singh lab for helpful comments and suggestions. We also thank staff of the University of Liverpool Centre for Genomic Research for the generation of sequence data. The work was supported by NIH grants AI085178 and AI053724 to US. NH was supported by a Wolfson Merit Award from the Royal Society of Great Britain.
- Agrawal N, Dasaradhi PV, Mohmmed A, Malhotra P, Bhatnagar RK, Mukherjee SK: RNA interference: biology, mechanism, and applications. Microbiol Mol Biol Rev. 2003, 67 (4): 657-685. 10.1128/MMBR.67.4.657-685.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Cerutti H, Casas-Mollano JA: On the origin and functions of RNA-mediated silencing: from protists to man. Curr Genet. 2006, 50 (2): 81-99. 10.1007/s00294-006-0078-x.PubMed CentralView ArticlePubMedGoogle Scholar
- Ambros V: MicroRNAs and developmental timing. Curr Opin Genet Dev. 2011, 21 (4): 511-517. 10.1016/j.gde.2011.04.003.PubMed CentralView ArticlePubMedGoogle Scholar
- Aliyari R, Ding SW: RNA-based viral immunity initiated by the Dicer family of host immune receptors. Immunol Rev. 2009, 227 (1): 176-188. 10.1111/j.1600-065X.2008.00722.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Schoeberl UE, Mochizuki K: Keeping the soma free of transposons: programmed DNA elimination in ciliates. J Biol Chem. 2011, 286 (43): 37045-37052. 10.1074/jbc.R111.276964.PubMed CentralView ArticlePubMedGoogle Scholar
- Ghildiyal M, Zamore PD: Small silencing RNAs: an expanding universe. Nat Rev. 2009, 10 (2): 94-108. 10.1038/nrg2504.View ArticleGoogle Scholar
- Han H, Nutiu R, Moffat J, Blencowe BJ: SnapShot: High-throughput sequencing applications. Cell. 2011, 146 (6): 1044-10.1016/j.cell.2011.09.002. 1044 e1041-1042View ArticlePubMedGoogle Scholar
- Ullu E, Tschudi C, Chakraborty T: RNA interference in protozoan parasites. Cell Microbiol. 2004, 6 (6): 509-519. 10.1111/j.1462-5822.2004.00399.x.View ArticlePubMedGoogle Scholar
- Kolev NG, Tschudi C, Ullu E: RNA interference in protozoan parasites: achievements and challenges. Eukaryot Cell. 2011, 10 (9): 1156-1163. 10.1128/EC.05114-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Patrick KL, Shi H, Kolev NG, Ersfeld K, Tschudi C, Ullu E: Distinct and overlapping roles for two Dicer-like proteins in the RNA interference pathways of the ancient eukaryote Trypanosoma brucei. Proc Natl Acad Sci U S A. 2009, 106 (42): 17933-17938. 10.1073/pnas.0907766106.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang H, Ehrenkaufer GM, Pompey JM, Hackney JA, Singh U: Small RNAs with 5′-polyphosphate termini associate with a Piwi-related protein and regulate gene expression in the single-celled eukaryote Entamoeba histolytica. PLoS Pathogens. 2008, 4 (11): e1000219-10.1371/journal.ppat.1000219.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang H, Alramini H, Tran V, Singh U: Nucleus-localized antisense small RNAs with 5′-polyphosphate termini regulate long term transcriptional gene silencing in Entamoeba histolytica G3 strain. J Biol Chem. 2011, 286 (52): 44467-44479. 10.1074/jbc.M111.278184.PubMed CentralView ArticlePubMedGoogle Scholar
- Saraiya AA, Li W, Wang CC: A microRNA derived from an apparent canonical biogenesis pathway regulates variant surface protein gene expression in Giardia lamblia. RNA (New York, NY). 2011, 17 (12): 2152-2164. 10.1261/rna.028118.111.View ArticleGoogle Scholar
- Prucca CG, Slavin I, Quiroga R, Elias EV, Rivero FD, Saura A, Carranza PG, Lujan HD: Antigenic variation in Giardia lamblia is regulated by RNA interference. Nature. 2008, 456 (7223): 750-754. 10.1038/nature07585.View ArticlePubMedGoogle Scholar
- Braun L, Cannella D, Ortet P, Barakat M, Sautel CF, Kieffer S, Garin J, Bastien O, Voinnet O, Hakimi MA: A complex small RNA repertoire is generated by a plant/fungal-like machinery and effected by a metazoan-like Argonaute in the single-cell human parasite Toxoplasma gondii. PLoS Pathogens. 2010, 6 (5): e1000920-10.1371/journal.ppat.1000920.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen XS, Collins LJ, Biggs PJ, Penny D: High throughput genome-wide survey of small RNAs from the parasitic protists Giardia intestinalis and Trichomonas vaginalis. Genome Biol Evol. 2009, 1: 165-175.View ArticlePubMedGoogle Scholar
- Stanley SL: Amoebiasis. Lancet. 2003, 361 (9362): 1025-1034. 10.1016/S0140-6736(03)12830-9.View ArticlePubMedGoogle Scholar
- Zhang H, Pompey JM, Singh U: RNA interference in Entamoeba histolytica: implications for parasite biology and gene silencing. Future Microbiol. 2011, 6 (1): 103-117. 10.2217/fmb.10.154.PubMed CentralView ArticlePubMedGoogle Scholar
- Kaur G, Lohia A: Inhibition of gene expression with double strand RNA interference in Entamoeba histolytica. Biochem Biophys Res Commun. 2004, 320 (4): 1118-1122. 10.1016/j.bbrc.2004.06.064.View ArticlePubMedGoogle Scholar
- Vayssie L, Vargas M, Weber C, Guillen N: Double-stranded RNA mediates homology-dependent gene silencing of gamma-tubulin in the human parasite Entamoeba histolytica. Mol Biochem Parasitol. 2004, 138 (1): 21-28. 10.1016/j.molbiopara.2004.07.005.View ArticlePubMedGoogle Scholar
- Linford AS, Moreno H, Good KR, Zhang H, Singh U, Petri WA: Short hairpin RNA-mediated knockdown of protein expression in Entamoeba histolytica. BMC Microbiol. 2009, 9: 38-10.1186/1471-2180-9-38.PubMed CentralView ArticlePubMedGoogle Scholar
- Bracha R, Nuchamowitz Y, Anbar M, Mirelman D: Transcriptional silencing of multiple genes in trophozoites of Entamoeba histolytica. PLoS Pathogens. 2006, 2 (5): e48-10.1371/journal.ppat.0020048.PubMed CentralView ArticlePubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10 (3): R25-10.1186/gb-2009-10-3-r25.PubMed CentralView ArticlePubMedGoogle Scholar
- Djikeng A, Shi H, Tschudi C, Ullu E: RNA interference in Trypanosoma brucei: cloning of small interfering RNAs provides evidence for retroposon-derived 24-26-nucleotide RNAs. RNA (New York, NY). 2001, 7 (11): 1522-1530.Google Scholar
- Tschudi C, Shi H, Franklin JB, Ullu E: Small interfering RNA-producing loci in the ancient parasitic eukaryote Trypanosoma brucei. BMC Genomics. 2012, 13: 427-10.1186/1471-2164-13-427.PubMed CentralView ArticlePubMedGoogle Scholar
- Aravin AA, Lagos-Quintana M, Yalcin A, Zavolan M, Marks D, Snyder B, Gaasterland T, Meyer J, Tuschl T: The small RNA profile during Drosophila melanogaster development. Dev Cell. 2003, 5 (2): 337-350. 10.1016/S1534-5807(03)00228-4.View ArticlePubMedGoogle Scholar
- Djupedal I, Kos-Braun IC, Mosher RA, Soderholm N, Simmer F, Hardcastle TJ, Fender A, Heidrich N, Kagansky A, Bayne E: Analysis of small RNA in fission yeast; centromeric siRNAs are potentially generated through a structured RNA. EMBO J. 2009, 28 (24): 3832-3844. 10.1038/emboj.2009.351.PubMed CentralView ArticlePubMedGoogle Scholar
- Clark CG, Alsmark UC, Tazreiter M, Saito-Nakano Y, Ali V, Marion S, Weber C, Mukherjee C, Bruchhaus I, Tannich E: Structure and content of the Entamoeba histolytica genome. Adv Parasitol. 2007, 65: 51-190.View ArticlePubMedGoogle Scholar
- Pak J, Fire A: Distinct populations of primary and secondary effectors during RNAi in C. elegans. Science. 2007, 315 (5809): 241-244. 10.1126/science.1132839.View ArticlePubMedGoogle Scholar
- Wang J, Czech B, Crunk A, Wallace A, Mitreva M, Hannon GJ, Davis RE: Deep small RNA sequencing from the nematode Ascaris reveals conservation, functional diversification, and novel developmental profiles. Genome Res. 2011, 21 (9): 1462-1477. 10.1101/gr.121426.111.PubMed CentralView ArticlePubMedGoogle Scholar
- Loftus B, Anderson I, Davies R, Alsmark UCM, Samuelson J, Amedeo P, Roncaglia P, Berriman M, Hirt RP, Mann BJ: The genome of the protist parasite Entamoeba histolytica. Nature. 2005, 433 (7028): 865-868. 10.1038/nature03291.View ArticlePubMedGoogle Scholar
- Lorenzi HA, Puiu D, Miller JR, Brinkac LM, Amedeo P, Hall N, Caler EV: New assembly, reannotation and analysis of the Entamoeba histolytica genome reveal new genomic features and protein content information. PLoS Negl Trop Dis. 2010, 4 (6): e716-10.1371/journal.pntd.0000716.PubMed CentralView ArticlePubMedGoogle Scholar
- Gu W, Shirayama M, Conte D, Gu W, Shirayama M, Conte D, Vasale J, Batista PJ, Claycomb JM, Moresco JJ, Youngman EM, Keys J, Stoltz MJ: Distinct argonaute-mediated 22 G-RNA pathways direct genome surveillance in the C. elegans germline. Mol Cell. 2009, 36 (2): 231-244. 10.1016/j.molcel.2009.09.020.PubMed CentralView ArticlePubMedGoogle Scholar
- Elbashir SM, Harborth J, Lendeckel W, Yalcin A, Weber K, Tuschl T: Duplexes of 21-nucleotide RNAs mediate RNA interference in cultured mammalian cells. Nature. 2001, 411 (6836): 494-498. 10.1038/35078107.View ArticlePubMedGoogle Scholar
- Gan J, Tropea JE, Austin BP, Court DL, Waugh DS, Ji X: Structural insight into the mechanism of double-stranded RNA processing by ribonuclease III. Cell. 2006, 124 (2): 355-366. 10.1016/j.cell.2005.11.034.View ArticlePubMedGoogle Scholar
- Okamura K, Balla S, Martin R, Liu N, Lai EC: Two distinct mechanisms generate endogenous siRNAs from bidirectional transcription in Drosophila melanogaster. Nat Struct Mol Biol. 2008, 15 (9): 998-View ArticlePubMedGoogle Scholar
- Borsani O, Zhu J, Verslues PE, Sunkar R, Zhu JK: Endogenous siRNAs derived from a pair of natural cis-antisense transcripts regulate salt tolerance in Arabidopsis. Cell. 2005, 123 (7): 1279-1291. 10.1016/j.cell.2005.11.035.PubMed CentralView ArticlePubMedGoogle Scholar
- Teodorovic S, Walls CD, Elmendorf HG: Bidirectional transcription is an inherent feature of Giardia lamblia promoters and contributes to an abundance of sterile antisense transcripts throughout the genome. Nucleic Acids Res. 2007, 35 (8): 2544-2553. 10.1093/nar/gkm105.PubMed CentralView ArticlePubMedGoogle Scholar
- Maniar JM, Fire AZ: EGO-1, a C. elegans RdRP, modulates gene expression via production of mRNA-templated short antisense RNAs. Current Biol: CB. 2011, 21 (6): 449-459. 10.1016/j.cub.2011.02.019.PubMed CentralView ArticleGoogle Scholar
- Clark CG, Ali IK, Zaki M, Loftus BJ, Hall N: Unique organisation of tRNA genes in Entamoeba histolytica. Mol Biochem Parasitol. 2006, 146 (1): 24-29. 10.1016/j.molbiopara.2005.10.013.View ArticlePubMedGoogle Scholar
- Bhattacharya S, Som I, Bhattacharya A: The ribosomal DNA plasmids of entamoeba. Parasitol Today. 1998, 14 (5): 181-185. 10.1016/S0169-4758(98)01222-8.View ArticlePubMedGoogle Scholar
- Ender C, Krek A, Friedlander MR, Beitzinger M, Weinmann L, Chen W, Pfeffer S, Rajewsky N, Meister G: A human snoRNA with microRNA-like functions. Mol Cell. 2008, 32 (4): 519-528. 10.1016/j.molcel.2008.10.017.View ArticlePubMedGoogle Scholar
- Cole C, Sobala A, Lu C, Thatcher SR, Bowman A, Brown JW, Green PJ, Barton GJ, Hutvagner G: Filtering of deep sequencing data reveals the existence of abundant Dicer-dependent small RNAs derived from tRNAs. RNA (New York, NY). 2009, 15 (12): 2147-2160. 10.1261/rna.1738409.View ArticleGoogle Scholar
- Bakre AA, Rawal K, Ramaswamy R, Bhattacharya A, Bhattacharya S: The LINEs and SINEs of Entamoeba histolytica: comparative analysis and genomic distribution. Exp Parasitol. 2005, 110 (3): 207-213. 10.1016/j.exppara.2005.02.009.View ArticlePubMedGoogle Scholar
- Lorenzi H, Thiagarajan M, Haas B, Wortman J, Hall N, Caler E: Genome wide survey, discovery and evolution of repetitive elements in three Entamoeba species. BMC Genom. 2008, 9: 595-10.1186/1471-2164-9-595.View ArticleGoogle Scholar
- Sharma R, Bagchi A, Bhattacharya A, Bhattacharya S: Characterization of a retrotransposon-like element from Entamoeba histolytica. Mol Biochem Parasitol. 2001, 116 (1): 45-53. 10.1016/S0166-6851(01)00300-0.View ArticlePubMedGoogle Scholar
- Girard A, Hannon GJ: Conserved themes in small-RNA-mediated transposon control. Trends Cell Biol. 2008, 18 (3): 136-148. 10.1016/j.tcb.2008.01.004.PubMed CentralView ArticlePubMedGoogle Scholar
- Ketting RF: The many faces of RNAi. Dev Cell. 2011, 20 (2): 148-161. 10.1016/j.devcel.2011.01.012.View ArticlePubMedGoogle Scholar
- Ehrenkaufer GM, Haque R, Hackney JA, Eichinger DJ, Singh U: Identification of developmentally regulated genes in Entamoeba histolytica: insights into mechanisms of stage conversion in a protozoan parasite. Cell Microbiol. 2007, 9 (6): 1426-1444. 10.1111/j.1462-5822.2006.00882.x.View ArticlePubMedGoogle Scholar
- Ehrenkaufer GM, Eichinger DJ, Singh U: Trichostatin A effects on gene expression in the protozoan parasite Entamoeba histolytica. BMC Genom. 2007, 8: 216-10.1186/1471-2164-8-216.View ArticleGoogle Scholar
- Ali IK, Ehrenkaufer GM, Hackney JA, Singh U: Growth of the protozoan parasite Entamoeba histolytica in 5-azacytidine has limited effects on parasite gene expression. BMC Genom. 2007, 8: 7-10.1186/1471-2164-8-7.View ArticleGoogle Scholar
- Hackney JA, Ehrenkaufer GM, Singh U: Identification of putative transcriptional regulatory networks in Entamoeba histolytica using Bayesian inference. Nucleic Acids Res. 2007, 35 (7): 2141-2152. 10.1093/nar/gkm028.PubMed CentralView ArticlePubMedGoogle Scholar
- Gilchrist CA, Houpt E, Trapaidze N, Fei Z, Crasta O, Asgharpour A, Evans C, Martino-Catt S, Baba DJ, Stroup S: Impact of intestinal colonization and invasion on the Entamoeba histolytica transcriptome. Mol Biochem Parasitol. 2006, 147 (2): 163-176. 10.1016/j.molbiopara.2006.02.007.View ArticlePubMedGoogle Scholar
- Dvorak JA, Kobayashi S, Nozaki T, Takeuchi T, Matsubara C: Induction of permeability changes and death of vertebrate cells is modulated by the virulence of Entamoeba spp. isolates. Parasitol Int. 2003, 52 (2): 169-173. 10.1016/S1383-5769(02)00090-9.View ArticlePubMedGoogle Scholar
- MacFarlane RC, Singh U: Identification of differentially expressed genes in virulent and nonvirulent Entamoeba species: potential implications for amebic pathogenesis. Infect Immun. 2006, 74 (1): 340-351. 10.1128/IAI.74.1.340-351.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Shah PH, MacFarlane RC, Bhattacharya D, Matese JC, Demeter J, Stroup SE, Singh U: Comparative genomic hybridizations of Entamoeba strains reveal unique genetic fingerprints that correlate with virulence. Eukaryot Cell. 2005, 4 (3): 504-515. 10.1128/EC.4.3.504-515.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- MacFarlane RC, Singh U: Identification of an Entamoeba histolytica serine-, threonine-, and isoleucine-rich protein with roles in adhesion and cytotoxicity. Eukaryot Cell. 2007, 6 (11): 2139-2146. 10.1128/EC.00174-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Pak J, Maniar JM, Mello CC, Fire A: Protection from feed-forward amplification in an amplified RNAi mechanism. Cell. 2012, 151 (4): 885-899. 10.1016/j.cell.2012.10.022.PubMed CentralView ArticlePubMedGoogle Scholar
- Burton NO, Burkhart KB, Kennedy S: Nuclear RNAi maintains heritable gene silencing in Caenorhabditis elegans. Proc Natl Acad Sci U S A. 2011, 108 (49): 19683-19688. 10.1073/pnas.1113310108.PubMed CentralView ArticlePubMedGoogle Scholar
- Buckley BA, Burkhart KB, Gu SG, Spracklin G, Kershner A, Fritz H, Kimble J, Fire A, Kennedy S: A nuclear Argonaute promotes multigenerational epigenetic inheritance and germline immortality. Nature. 2012, 489 (7416): 447-451. 10.1038/nature11352.PubMed CentralView ArticlePubMedGoogle Scholar
- Guang S, Bochner AF, Pavelec DM, Burkhart KB, Harding S, Lachowiec J, Kennedy S: An Argonaute transports siRNAs from the cytoplasm to the nucleus. Science. 2008, 321 (5888): 537-541. 10.1126/science.1157647.PubMed CentralView ArticlePubMedGoogle Scholar
- Huntley DM, Pandis I, Butcher SA, Ackers JP: Bioinformatic analysis of Entamoeba histolytica SINE1 elements. BMC Genom. 2010, 11: 321-10.1186/1471-2164-11-321.View ArticleGoogle Scholar