The ribosomal RNA transcription unit of Entamoeba invadens: Accumulation of unprocessed pre-rRNA and a long non coding RNA during encystation

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.

A landscape of gene regulation in the parasitic amoebozoa Entamoeba spp

Entamoeba are amoeboid extracellular parasites that represent an important group of organisms for which the regulatory networks must be examined to better understand how genes and functional processes are interrelated. In this work, we inferred the gene regulatory networks (GRNs) in four Entamoeba species, E. histolytica, E. dispar, E. nuttalli, and E. invadens, and the GRN topological properties and the corresponding biological functions were evaluated. From these analyses, we determined that transcription factors (TFs) of E. histolytica, E. dispar, and E. nuttalli are associated mainly with the LIM family, while the TFs in E. invadens are associated with the RRM_1 family. In addition, we identified that EHI_044890 regulates 121 genes in E. histolytica, EDI_297980 regulates 284 genes in E. dispar, ENU1_120230 regulates 195 genes in E. nuttalli, and EIN_249270 regulates 257 genes in E. invadens. Finally, we identified that three types of processes, Macromolecule metabolic process, Cellular macromolecule metabolic process, and Cellular nitrogen compound metabolic process, are the main biological processes for each network. The results described in this work can be used as a basis for the study of gene regulation in these organisms.

Comparative analysis of the ribosomal DNA repeat unit (rDNA) of Perna viridis (Linnaeus, 1758) and Perna canaliculus (Gmelin, 1791)

Perna viridis and P. canaliculus are economically and ecologically important species of shellfish. In this study, the complete ribosomal DNA (rDNA) unit sequences of these species were determined for the first time. The gene order, 18S rRNA-internal transcribed spacer (ITS) 1-5.8S rRNA-ITS2-28S rRNA-intergenic spacer (IGS), was similar to that observed in other eukaryotes. The lengths of the P. viridis and P. canaliculus rDNA sequences ranged from 8,432 to 8,616 bp and from 7,597 to 7,610 bp, respectively, this variability was mainly attributable to the IGS region. The putative transcription termination site and initiation site were confirmed. Perna viridis and P. canaliculus rDNA contained two (length: 93 and 40 bp) and one (length: 131 bp) repeat motifs, respectively. Individual intra-species differences mainly involved the copy number of repeat units. In P. viridis, three cytosine-guanine (CpG) sites with sizes of 440, 1,075 and 537 bp were found to cover nearly the entire IGS sequence, whereas in P. canaliculus, two CpG islands with sizes of 361 and 484 bp were identified. The phylogenetic trees constructed with maximum likelihood and neighbour-joining methods and based on ITS sequences were identical and included three major clusters. Species of the same genus were easily clustered together.

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.