Self-circularizing 5'-ETS RNAs accumulate along with unprocessed pre ribosomal RNAs in growth-stressed Entamoeba histolytica

In Silico Identification and Characterization of circRNAs as Potential Virulence-Related miRNA/siRNA Sponges from Entamoeba histolytica and Encystment-Related circRNAs from Entamoeba invadens

Ubiquitous eukaryotic non-coding circular RNAs regulate transcription and translation. We have reported full-length intronic circular RNAs (flicRNAs) in Entamoeba histolytica with esterified 3′ss and 5′ss. Their 5′ss GU-rich elements are essential for their biogenesis and their suggested role in transcription regulation. Here, we explored whether exonic, exonic-intronic, and intergenic circular RNAs are also part of the E. histolytica and E. invadens ncRNA RNAome and investigated their possible functions. Available RNA-Seq libraries were analyzed with the CIRI-full software in search of circular exonic RNAs (circRNAs). The robustness of the analyses was validated using synthetic decoy sequences with bona fide back splice junctions. Differentially expressed (DE) circRNAs, between the virulent HM1:IMSS and the nonvirulent Rahman E. histolytica strains, were identified, and their miRNA sponging potential was analyzed using the intaRNA software. Respectively, 188 and 605 reverse overlapped circRNAs from E. invadens and E. histolytica were identified. The sequence composition of the circRNAs was mostly exonic although different to human circRNAs in other attributes. 416 circRNAs from E. histolytica were virulent-specific and 267 were nonvirulent-specific. Out of the common circRNAs, 32 were DE between strains. Finally, we predicted that 8 of the DE circRNAs could function as sponges of the bioinformatically reported miRNAs in E. histolytica, whose functions are still unknown. Our results extend the E. histolytica RNAome and allow us to devise a hypothesis to test circRNAs/miRNAs/siRNAs interactions in determining the virulent/nonvirulent phenotypes and to explore other regulatory mechanisms during amoebic encystment.

An extended catalogue of ncRNAs in Streptomyces coelicolor reporting abundant tmRNA, RNase-P RNA and RNA fragments derived from pre-ribosomal RNA leader sequences

Streptomyces coelicolor is a model organism for studying streptomycetes. This genus possesses relevant medical and economical roles, because it produces many biologically active metabolites of pharmaceutical interest, including the majority of commercialized antibiotics. In this bioinformatic study, the transcriptome of S. coelicolor has been analyzed to identify novel RNA species and quantify the expression of both annotated and novel transcripts in solid and liquid growth medium cultures at different times. The major characteristics disclosed in this study are: (i) the diffuse antisense transcription; (ii) the great abundance of transfer-messenger RNAs (tmRNA); (iii) the abundance of rnpB transcripts, paramount for the RNase-P complex; and (iv) the presence of abundant fragments derived from pre-ribosomal RNA leader sequences of unknown biological function. Overall, this study extends the catalogue of ncRNAs in S. coelicolor and suggests an important role of non-coding transcription in the regulation of biologically active molecule production.

Transcriptomic analysis of ribosome biogenesis and pre-rRNA processing during growth stress in Entamoeba histolytica

Ribosome biogenesis, a multi-step process involving transcription, modification, folding and processing of rRNA, is the major consumer of cellular energy. It involves sequential assembly of ribosomal proteins (RP)s via more than 200 ribogenesis factors. Unlike model organisms where transcription of rRNA and RP genes slows down during stress, in Entamoeba histolytica, pre-rRNA synthesis continues, and unprocessed pre-rRNA accumulates. Northern hybridization from different spacer regions depicted the accumulation of unprocessed intermediates during stress. To gain insight into the vast repertoire of ribosome biogenesis factors and understand the major components playing role during stress we computationally identified ribosome biogenesis factors in E. histolytica. Of the ∼279 Saccharomyces cerevisiae proteins, we could only find 188 proteins in E. histolytica. Some of the proteins missing in E. histolytica were also missing in humans. A number of proteins represented by multiple genes in S. cerevisiae had a single copy in E. histolytica. Interestingly E. histolytica lacked mitochondrial ribosome biogenesis factors and had far less RNase components compared to S. cerevisiae. Transcriptomic studies revealed the differential regulation of ribosomal factors both in serum starved and RRP6 down-regulation conditions. These included the NEP1 and TSR3 proteins that chemically modify 18S-rRNA. Pre-rRNA precursors accumulate upon downregulation of the latter proteins in S. cerevisiae and humans. These data reveal the major factors that regulate pre-rRNA processing during stress in E. histolytica and provide the first complete repertoire of ribosome biogenesis factors in this early-branching protist.

Advances in Entamoeba histolytica Biology Through Transcriptomic Analysis

A large number of transcriptome-level studies in Entamoeba histolytica, the protozoan parasite that causes amoebiasis, have investigated gene expression patterns to help understand the pathology and biology of the organism. They have compared virulent and avirulent strains in lab culture and after tissue invasion, cells grown under different stress conditions, response to anti-amoebic drug treatments, and gene expression changes during the process of encystation. These studies have revealed interesting molecules/pathways that will help increase our mechanistic understanding of differentially expressed genes during growth perturbations and tissue invasion. Some of the important insights obtained from transcriptome studies include the observations that regulation of carbohydrate metabolism may be an important determinant for tissue invasion, while the novel up-regulated genes during encystation include phospholipase D, and meiotic genes, suggesting the possibility of meiosis during the process. Classification of genes according to expression levels showed that amongst the highly transcribed genes in cultured E. histolytica trophozoites were some virulence factors, raising the question of the role of these factors in normal parasite growth. Promoter motifs associated with differential gene expression and regulation were identified. Some of these motifs associated with high gene expression were located downstream of start codon, and were required for efficient transcription. The listing of E. histolytica genes according to transcript expression levels will help us determine the scale of post-transcriptional regulation, and the possible roles of predicted promoter motifs. The small RNA transcriptome is a valuable resource for detailed structural and functional analysis of these molecules and their regulatory roles. These studies provide new drug targets and enhance our understanding of gene regulation in E. histolytica.

Transcriptomic analysis reveals novel downstream regulatory motifs and highly transcribed virulence factor genes of Entamoeba histolytica

Background

Promoter motifs in Entamoeba histolytica were earlier analysed using microarray data with lower dynamic range of gene expression. Additionally, previous transcriptomic studies did not provide information on the nature of highly transcribed genes, and downstream promoter motifs important for gene expression. To address these issues we generated RNA-Seq data and identified the high and low expressing genes, especially with respect to virulence potential. We analysed sequences both upstream and downstream of start site for important motifs.

Results

We used RNA-Seq data to classify genes according to expression levels, which ranged six orders of magnitude. Data were validated by reporter gene expression. Virulence-related genes (except AIG1) were amongst the highly expressed, while some kinases and BspA family genes were poorly expressed. We looked for conserved motifs in sequences upstream and downstream of the initiation codon. Following enrichment by AME we found seven motifs significantly enriched in high expression- and three in low expression-classes. Two of these motifs (M4 and M6) were located downstream of AUG, were exclusively enriched in high expression class, and were mostly found in ribosomal protein, and translation-related genes. Motif deletion resulted in drastic down regulation of reporter gene expression, showing functional relevance. Distribution of core promoter motifs (TATA, GAAC, and Inr) in all genes revealed that genes with downstream motifs were not preferentially associated with TATA-less promoters. We looked at gene expression changes in cells subjected to growth stress by serum starvation, and experimentally validated the data. Genes showing maximum up regulation belonged to the low or medium expression class, and included genes in signalling pathways, lipid metabolism, DNA repair, Myb transcription factors, BspA, and heat shock. Genes showing maximum down regulation belonged to the high or medium expression class. They included genes for signalling factors, actin, Ariel family, and ribosome biogenesis factors.

Conclusion

Our analysis has added important new information about the E. histolytica transcriptome. We report for the first time two downstream motifs required for gene expression, which could be used for over expression of E. histolytica genes. Most of the virulence-related genes in this parasite are highly expressed in culture.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5570-z) contains supplementary material, which is available to authorized users.

Stress-induced nuclear depletion of Entamoeba histolytica 3'-5' exoribonuclease EhRrp6 and its role in growth and erythrophagocytosis

The 3'-5' exoribonuclease Rrp6 is a key enzyme in RNA homeostasis involved in processing and degradation of many stable RNA precursors, aberrant transcripts, and noncoding RNAs. We previously have shown that in the protozoan parasite Entamoeba histolytica, the 5'-external transcribed spacer fragment of pre-rRNA accumulates under serum starvation-induced growth stress. This fragment is a known target of degradation by Rrp6. Here, we computationally and biochemically characterized EhRrp6 and found that it contains the catalytically important EXO and HRDC domains and exhibits exoribonuclease activity with both unstructured and structured RNA substrates, which required the conserved DEDD-Y catalytic-site residues. It lacked the N-terminal PMC2NT domain for binding of the cofactor Rrp47, but could functionally complement the growth defect of a yeast rrp6 mutant. Of note, no Rrp47 homologue was detected in E. histolytica Immunolocalization studies revealed that EhRrp6 is present both in the nucleus and cytosol of normal E. histolytica cells. However, growth stress induced its complete loss from the nuclei, reversed by proteasome inhibitors. EhRrp6-depleted E. histolytica cells were severely growth restricted, and EhRrp6 overexpression protected the cells against stress, suggesting that EhRrp6 functions as a stress sensor. Importantly EhRrp6 depletion reduced erythrophagocytosis, an important virulence determinant of E. histolytica This reduction was due to a specific decrease in transcript levels of some phagocytosis-related genes (Ehcabp3 and Ehrho1), whereas expression of other genes (Ehcabp1, Ehcabp6, Ehc2pk, and Eharp2/3) was unaffected. This is the first report of the role of Rrp6 in cell growth and stress responses in a protozoan parasite.

Postsplicing-Derived Full-Length Intron Circles in the Protozoan Parasite Entamoeba histolytica

Noncoding circular RNAs are widespread in the tree of life. Particularly, intron-containing circular RNAs which apparently upregulate their parental gene expression. Entamoeba histolytica, the causative agent of dysentery and liver abscesses in humans, codes for several noncoding RNAs, including circular ribosomal RNAs, but no intron containing circular RNAs have been described to date. Divergent RT-PCR and diverse molecular approaches, allowed us to detect bona fide full-length intronic circular RNA (flicRNA) molecules. Self-splicing reactions, RNA polymerase II inhibition with Actinomycin D, and second step of splicing-inhibition with boric acid showed that the production of flicRX13 (one of the flicRNAs found in this work, and our test model) depends on mRNA synthesis and pre-mRNA processing instead of self-splicing. To explore the cues and factors involved in flicRX13 biogenesis in vivo, splicing assays were carried out in amoeba transformants where splicing factors and Dbr1 (intron lariat debranching enzyme 1) were silenced or overexpressed, or where Rabx13 wild-type and mutant 5'ss (splice site) and branch site minigene constructs were overexpressed. Whereas SF1 (splicing factor 1) is not involved, the U2 auxiliary splicing factor, Dbr1, and the GU-rich 5'ss are involved in postsplicing flicRX13 biogenesis, probably by Dbr1 stalling, in a similar fashion to the formation of ciRNAs (circular intronic RNAs), but with distinctive 5'-3'ss ligation points. Different from the reported functions of ciRNAs, the 5'ss GU-rich element of flicRX13 possibly interacts with transcription machinery to silence its own gene in cis. Furthermore, introns of E. histolytica virulence-related genes are also processed as flicRNAs.

Application of SHAPE reveals in vivo RNA folding under normal and growth-stressed conditions in the human parasite Entamoeba histolytica

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) is a versatile sequence independent method to probe RNA structure in vivo and in vitro. It has so far been tried mainly with model organisms. We show that cells of Entamoeba histolytica, a protozoan parasite of humans are hyper-sensitive to the in vivo SHAPE reagent, NAI, and show rapid loss of viability and RNA integrity. We optimized treatment conditions with 5.8S rRNA and Eh_U3 snoRNA to obtain NAI-modification while retaining RNA integrity. The modification patterns were highly reproducible. The in vivo folding was different from in vitro and correlated well with known interactions of 5.8S rRNA with proteins in vivo. The Eh_U3 snoRNA also showed many differences in its in vivo versus in vitro folding, which correlated with conserved interactions of this RNA with 18S rRNA and 5'-ETS. Further, Eh_U3 snoRNA obtained from serum-starved cells showed an open 3'-hinge structure, indicating disruption of 5'-ETS interaction. This could contribute to the observed slow processing of pre-rRNA in starved cells. Our work shows the applicability of SHAPE to study in vivo RNA folding in a parasite and will encourage the use of this reagent for RNA structure analysis in other such organisms.

Identification of EhTIF-IA: The putative E. histolytica orthologue of the human ribosomal RNA transcription initiation factor-IA

Initiation of rDNA transcription requires the assembly of a specific multi-protein complex at the rDNA promoter containing the RNA Pol I with auxiliary factors. One of these factors is known as Rrn3P in yeast and Transcription Initiation Factor IA (TIF-IA) in mammals. Rrn3p/TIF-IA serves as a bridge between RNA Pol I and the pre-initiation complex at the promoter. It is phosphorylated at multiple sites and is involved in regulation of rDNA transcription in a growth-dependent manner. In the early branching parasitic protist Entamoeba histolytica, the rRNA genes are present exclusively on circular extra chromosomal plasmids. The protein factors involved in regulation of rDNA transcription in E. histolytica are not known. We have identified the E. histolytica equivalent of TIF-1A (EhTIF-IA) by homology search within the database and was further cloned and expressed. Immuno-localization studies showed that EhTIF-IA co-localized partially with fibrillarin in the peripherally localized nucleolus. EhTIF-IA was shown to interact with the RNA Pol I-specific subunit RPA12 both in vivo and in vitro. Mass spectroscopy data identified RNA Pol I-specific subunits and other nucleolar proteins to be the interacting partners of EhTIF-IA. Our study demonstrates for the first time a conserved putative RNA Pol I transcription factor TIF-IA in E. histolytica.

Post-transcriptional regulation of ribosomal protein genes during serum starvation in Entamoeba histolytica

Ribosome synthesis involves all three RNA polymerases which are co-ordinately regulated to produce equimolar amounts of rRNAs and ribosomal proteins (RPs). Unlike model organisms where transcription of rRNA and RP genes slows down during stress, in E. histolytica rDNA transcription continues but pre-rRNA processing slows down and unprocessed pre-rRNA accumulates during serum starvation. To investigate the regulation of RP genes under stress we measured transcription of six selected RP genes from the small- and large-ribosomal subunits (RPS6, RPS3, RPS19, RPL5, RPL26, RPL30) representing the early-, mid-, and late-stages of ribosomal assembly. Transcripts of these genes persisted in growth-stressed cells. Expression of luciferase reporter under the control of two RP genes (RPS19 and RPL30) was studied during serum starvation and upon serum replenishment. Although luciferase transcript levels remained unchanged during starvation, luciferase activity steadily declined to 7.8% and 15% of control cells, respectively. After serum replenishment the activity increased to normal levels, suggesting post-transcriptional regulation of these genes. Mutations in the sequence -2 to -9 upstream of AUG in the RPL30 gene resulted in the phenotype expected of post-transcriptional regulation. Transcription of luciferase reporter was unaffected in this mutant, and luciferase activity did not decline during serum starvation, showing that this sequence is required to repress translation of RPL30 mRNA, and mutations in this region relieve repression. Our data show that during serum starvation E. histolytica blocks ribosome biogenesis post-transcriptionally by inhibiting pre-rRNA processing on the one hand, and the translation of RP mRNAs on the other.

Ribosomal RNA and protein transcripts persist in the cysts of Entamoeba invadens

In most organisms rDNA transcription ceases under conditions of growth stress. However, we have earlier shown that pre-rRNA accumulates during encystation in Entamoeba invadens. We labeled newly-synthesized rRNA during encystation, with [methyl-(3)H] methionine in the presence of chitinase to enable uptake of isotope. Incorporation rate reduced after 24h, and then increased to reach levels comparable with normal cells. The label was rapidly chased to the ribosomal pellet in dividing cells, while at late stages of encystation the ratio of counts going to the pellet dropped 3-fold. The transcript levels of selected ribosomal protein genes also went down initially but went up again at later stages of encystation. This suggested that rRNA and ribosomal protein transcription may be coordinately regulated. Our data shows that encysting E. invadens cells accumulate transcripts of both the RNA and protein components of the ribosome, which may ensure rapid synthesis of new ribosomes when growth resumes.

Analysis of U3 snoRNA and small subunit processome components in the parasitic protist Entamoeba histolytica

In the early branching parasitic protist Entamoeba histolytica, pre-rRNA synthesis continues when cells are subjected to growth stress, but processing slows down and unprocessed pre-rRNA accumulates. To gain insight into the regulatory mechanisms leading to accumulation, it is necessary to define the pre-rRNA processing machinery in E. histolytica. We searched the E. histolytica genome sequence for homologs of the SSU processome, which contains the U3snoRNA, and 72 proteins in yeast. We could identify 57 of the proteins with high confidence. Of the rest, 6 were absent in human, and 4 were non-essential in yeast. The remaining 5 were absent in other parasite genomes as well. Analysis of U3snoRNA showed that the E. histolytica U3snoRNA adopted the same conserved secondary structure as seen in yeast and human. The predicted structure was verified by chemical modification followed by primer extension (SHAPE). Further we showed that the predicted interactions of Eh_U3snoRNA boxes A and A' with pre-18S rRNA were highly conserved both in position and sequence. The predicted interactions of 5'-hinge and 3'-hinge sequences of Eh_U3 snoRNA with the 5'-ETS sequences were conserved in position but not in sequence. Transcription of selected genes of SSU processome was tested by northern analysis, and transcripts of predicted sizes were obtained. During serum starvation, when unprocessed pre-RNA accumulated, the transcript levels of some of these genes declined. This is the first report on pre-rRNA processing machinery in E. histolytica, and shows that the components are well conserved with respect to yeast and human.

Homologous Recombination Occurs in Entamoeba and Is Enhanced during Growth Stress and Stage Conversion

Homologous recombination (HR) has not been demonstrated in the parasitic protists Entamoeba histolytica or Entamoeba invadens, as no convenient method is available to measure it. However, HR must exist to ensure genome integrity, and possible genetic exchange, especially during stage conversion from trophozoite to cyst. Here we show the up regulation of mitotic and meiotic HR genes in Entamoeba during serum starvation, and encystation. To directly demonstrate HR we use a simple PCR-based method involving inverted repeats, which gives a reliable read out, as the recombination junctions can be determined by sequencing the amplicons. Using this read out, we demonstrate enhanced HR under growth stress in E. histolytica, and during encystation in E. invadens. We also demonstrate recombination between chromosomal inverted repeats. This is the first experimental demonstration of HR in Entamoeba and will help future investigations into this process, and to explore the possibility of meiosis in Entamoeba.

Long non-coding RNAs coordinate cellular responses to stress

Following the initial discovery of the heat shock RNA omega (hsrω) gene of Drosophila melanogaster to be non-coding (nc) and also inducible by cell stress, other stress-inducible long non-coding RNAs (lncRNA) have been described in diverse organisms. In view of the rapid sequence divergence of lncRNAs, present knowledge of stress trasncriptome is limited and fragmented. Several known stress-related lncRNAs, associated with specific nuclear speckled domains or nucleolus, provide structural base for sequestering diverse RNA-processing/regulatory proteins. Others have roles in transcriptional or translational inhibition during stress or in signaling pathways; functions of several other lncRNAs are not yet known. Most stress-related lncRNAs act primarily by modulating activity of the proteins to which they bind or by sequestering specific sets of proteins away from the active pool. A common emerging theme is that a given lncRNA targets one or more protein/s with key role/s in the cascade of events triggered by the stress and therefore has a widespread integrative effect. Since proteins associate with RNA through short sequence motifs, the overall base sequence of functionally similar ncRNAs is often not conserved except for specific motifs. The rapid evolvability of ncRNA sequences provides elegant modules for adaptability to changing environment as binding of one or the other protein to ncRNA can alter its structure and functions in distinct ways. Thus the stress-related lncRNAs act as hubs in the cellular networks to coordinate activities of the members within and between different networks to maintain cellular homeostasis for survival or to trigger cell death.