In silico analysis of EST and genomic sequences allowed the prediction of cis-regulatory elements for Entamoeba histolytica mRNA polyadenylation

RNA Sequencing Reveals Widespread Transcription of Natural Antisense RNAs in Entamoeba Species

Entamoeba is a genus of Amoebozoa that includes the intestine-colonizing pathogenic species Entamoeba histolytica. To understand the basis of gene regulation in E. histolytica from an evolutionary perspective, we have profiled the transcriptomes of its closely related species E. dispar, E. moshkovskii and E. invadens. Genome-wide identification of transcription start sites (TSS) and polyadenylation sites (PAS) revealed the similarities and differences of their gene regulatory sequences. In particular, we found the widespread initiation of antisense transcription from within the gene coding sequences is a common feature among all Entamoeba species. Interestingly, we observed the enrichment of antisense transcription in genes involved in several processes that are common to species infecting the human intestine, e.g., the metabolism of phospholipids. These results suggest a potentially conserved and compact gene regulatory system in Entamoeba.

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

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

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

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

Genomic determinants for initiation and length of natural antisense transcripts in Entamoeba histolytica

Natural antisense transcripts (NAT) have been reported in prokaryotes and eukaryotes. While the functions of most reported NATs remain unknown, their potentials in regulating the transcription of their counterparts have been speculated. Entamoeba histolytica, which is a unicellular eukaryotic parasite, has a compact protein-coding genome with very short intronic and intergenic regions. The regulatory mechanisms of gene expression in this compact genome are under-described. In this study, by genome-wide mapping of RNA-Seq data in the genome of E. histolytica, we show that a substantial fraction of its protein-coding genes (28%) has significant transcription on their opposite strand (i.e. NAT). Intriguingly, we found the location of transcription start sites or polyadenylation sites of NAT are determined by the specific motifs encoded on the opposite strand of the gene coding sequences, thereby providing a compact regulatory system for gene transcription. Moreover, we demonstrated that NATs are globally up-regulated under various environmental conditions including temperature stress and pathogenicity. While NATs do not appear to be consequences of spurious transcription, they may play a role in regulating gene expression in E. histolytica, a hypothesis which needs to be tested.

mRNA Polyadenylation Machineries in Intestinal Protozoan Parasites

In humans, mRNA polyadenylation involves the participation of about 20 factors in four main complexes that recognize specific RNA sequences. Notably, CFIm25, CPSF73, and PAP have essential roles for poly(A) site selection, mRNA cleavage, and adenosine residues polymerization. Besides the relevance of polyadenylation for gene expression, information is scarce in intestinal protozoan parasites that threaten human health. To better understand polyadenylation in Entamoeba histolytica, Giardia lamblia, and Cryptosporidium parvum, which represent leading causes of diarrhea worldwide, genomes were screened for orthologs of human factors. Results showed that Entamoeba histolytica and C. parvum have 16 and 12 proteins out of the 19 human proteins used as queries, respectively, while G. lamblia seems to have the smallest polyadenylation machinery with only six factors. Remarkably, CPSF30, CPSF73, CstF77, PABP2, and PAP, which were found in all parasites, could represent the core polyadenylation machinery. Multiple genes were detected for several proteins in Entamoeba, while gene redundancy is lower in Giardia and Cryptosporidium. Congruently with their relevance in the polyadenylation process, CPSF73 and PAP are present in all parasites, and CFIm25 is only missing in Giardia. They conserve the functional domains and predicted folding of human proteins, suggesting they may have the same roles in polyadenylation.

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.

Life and Death of mRNA Molecules in Entamoeba histolytica

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

Silencing the cleavage factor CFIm25 as a new strategy to control Entamoeba histolytica parasite

The 25 kDa subunit of the Clevage Factor Im (CFIm25) is an essential factor for messenger RNA polyadenylation in human cells. Therefore, here we investigated whether the homologous protein of Entamoeba histolytica, the protozoan responsible for human amoebiasis, might be considered as a biochemical target for parasite control. Trophozoites were cultured with bacterial double-stranded RNA molecules targeting the EhCFIm25 gene, and inhibition of mRNA and protein expression was confirmed by RT-PCR and Western blot assays, respectively. EhCFIm25 silencing was associated with a significant acceleration of cell proliferation and cell death. Moreover, trophozoites appeared as larger and multinucleated cells. These morphological changes were accompanied by a reduced mobility, and erythrophagocytosis was significantly diminished. Lastly, the knockdown of EhCFIm25 affected the poly(A) site selection in two reporter genes and revealed that EhCFIm25 stimulates the utilization of downstream poly(A) sites in E. histolytica mRNA. Overall, our data confirm that targeting the polyadenylation process represents an interesting strategy for controlling parasites, including E. histolytica. To our best knowledge, the present study is the first to have revealed the relevance of the cleavage factor CFIm25 as a biochemical target in parasites.

Untranslated regions of mRNA and their role in regulation of gene expression in protozoan parasites

Protozoan parasites are one of the oldest living entities in this world that throughout their existence have shown excellent resilience to the odds of survival and have adapted beautifully to ever changing rigors of the environment. In view of the dynamic environment encountered by them throughout their life cycle, and in establishing pathogenesis, it is unsurprising that modulation of gene expression plays a fundamental role in their survival. In higher eukaryotes, untranslated regions (UTRs) of transcripts are one of the crucial regulators of gene expression (influencing mRNA stability and translation efficiency). Parasitic protozoan genome studies have led to the characterization (in silico, in vitro and in vivo) of a large number of their genes. Comparison of higher eukaryotic UTRs with parasitic protozoan UTRs reveals the existence of several similar and dissimilar facets of the UTRs. This review focuses on the elements of UTRs of medically important protozoan parasites and their regulatory role in gene expression. Such information may be useful to researchers in designing gene targeting strategies linked with perturbation of host-parasite relationships leading to control of specific parasites.

Regulation of gene expression in the protozoan parasite Entamoeba invadens: identification of core promoter elements and promoters with stage-specific expression patterns

Developmental switching between life-cycle stages is a common feature among many pathogenic organisms. Entamoeba histolytica is an important human pathogen and is a leading parasitic cause of death globally. During its life cycle, Entamoeba converts between cysts (essential for disease transmission) and trophozoites (responsible for tissue invasion). Despite being central to its biology, the triggers that are involved in the developmental pathways of this parasite are not well understood. In order to define the transcriptional network associated with stage conversion we used Entamoeba invadens which serves as a model system for Entamoeba developmental biology, and performed RNA sequencing at different developmental time points. In this study RNA-Seq data was utilised to define basal transcriptional control elements as well as to identify promoters which regulate stage-specific gene expression patterns. We discovered that the 5' and 3' untranslated regions of E. invadens genes are short, a median of 20 nucleotides (nt) and 26 nt respectively. Bioinformatics analysis of DNA sequences proximate to the start and stop codons identified two conserved motifs: (i) E. invadens Core Promoter Motif - GAAC-Like (EiCPM-GL) (GAACTACAAA), and (ii) E. invadens 3'-U-Rich Motif (Ei3'-URM) (TTTGTT) in the 5' and 3' flanking regions, respectively. Electrophoretic mobility shift assays demonstrated that both motifs specifically bind nuclear protein(s) from E. invadens trophozoites. Additionally, we identified select genes with stage-specific expression patterns and analysed the ability of each gene promoter to drive a luciferase reporter gene during the developmental cycle. This approach confirmed three trophozoite-specific, four encystation-specific and two excystation-specific promoters. This work lays the framework for use of stage-specific promoters to express proteins of interest in a particular life-cycle stage, adding to the molecular toolbox for genetic manipulation of E. invadens and allowing further dissection of factors controlling Entamoeba developmental biology.

RNA polyadenylation sites on the genomes of microorganisms, animals, and plants

Pre–messenger RNA (mRNA) 3′-end cleavage and subsequent polyadenylation strongly regulate gene expression. In comparison with the upstream or downstream motifs, relatively little is known about the feature differences of polyadenylation [poly(A)] sites among major kingdoms. We suspect that the precise poly(A) sites are very selective, and we therefore mapped mRNA poly(A) sites on complete and nearly complete genomes using mRNA sequences available in the National Center for Biotechnology Information (NCBI) Nucleotide database. In this paper, we describe the mRNA nucleotide [i.e., the poly(A) tail attachment position] that is directly in attachment with the poly(A) tail and the pre-mRNA nucleotide [i.e., the poly(A) tail starting position] that corresponds to the first adenosine of the poly(A) tail in the 29 most-mapped species (2 fungi, 2 protists, 18 animals, and 7 plants). The most representative pre-mRNA dinucleotides covering these two positions were UA, CA, and GA in 17, 10, and 2 of the species, respectively. The pre-mRNA nucleotide at the poly(A) tail starting position was typically an adenosine [i.e., A-type poly(A) sites], sometimes a uridine, and occasionally a cytidine or guanosine. The order was U>C>G at the attachment position but A>>U>C≥G at the starting position. However, in comparison with the mRNA nucleotide composition (base composition), the poly(A) tail attachment position selected C over U in plants and both C and G over U in animals, in both A-type and non-A-type poly(A) sites. Animals, dicot plants, and monocot plants had clear differences in C/G ratios at the poly(A) tail attachment position of the non-A-type poly(A) sites. This study of poly(A) site evolution indicated that the two positions within poly(A) sites had distinct nucleotide compositions and were different among kingdoms.

The 25 kDa subunit of cleavage factor Im Is a RNA-binding protein that interacts with the poly(A) polymerase in Entamoeba histolytica

In eukaryotes, polyadenylation of pre-mRNA 3´ end is essential for mRNA export, stability and translation. Taking advantage of the knowledge of genomic sequences of Entamoeba histolytica, the protozoan responsible for human amoebiasis, we previously reported the putative polyadenylation machinery of this parasite. Here, we focused on the predicted protein that has the molecular features of the 25 kDa subunit of the Cleavage Factor Im (CFIm25) from other organisms, including the Nudix (nucleoside diphosphate linked to another moiety X) domain, as well as the RNA binding domain and the PAP/PAB interacting region. The recombinant EhCFIm25 protein (rEhCFIm25) was expressed in bacteria and used to generate specific antibodies in rabbit. Subcellular localization assays showed the presence of the endogenous protein in nuclear and cytoplasmic fractions. In RNA electrophoretic mobility shift assays, rEhCFIm25 was able to form specific RNA-protein complexes with the EhPgp5 mRNA 3´ UTR used as probe. In addition, Pull-Down and LC/ESI-MS/MS tandem mass spectrometry assays evidenced that the putative EhCFIm25 was able to interact with the poly(A) polymerase (EhPAP) that is responsible for the synthesis of the poly(A) tail in other eukaryotic cells. By Far-Western experiments, we confirmed the interaction between the putative EhCFIm25 and EhPAP in E. histolytica. Taken altogether, our results showed that the putative EhCFIm25 is a conserved RNA binding protein that interacts with the poly(A) polymerase, another member of the pre-mRNA 3´ end processing machinery in this protozoan parasite.

Quantification of stochastic noise of splicing and polyadenylation in Entamoeba histolytica

Alternative splicing and polyadenylation were observed pervasively in eukaryotic messenger RNAs. These alternative isoforms could either be consequences of physiological regulation or stochastic noise of RNA processing. To quantify the extent of stochastic noise in splicing and polyadenylation, we analyzed the alternative usage of splicing and polyadenylation sites in Entamoeba histolytica using RNA-Seq. First, we identified a large number of rarely spliced alternative junctions and then showed that the occurrence of these alternative splicing events is correlated with splicing site sequence, occurrence of constitutive splicing events and messenger RNA abundance. Our results implied the majority of these alternative splicing events are likely to be stochastic error of splicing machineries, and we estimated the corresponding error rates. Second, we observed extensive microheterogeneity of polyadenylation cleavage sites, and the extent of such microheterogeneity is correlated with the occurrence of constitutive cleavage events, suggesting most of such microheterogeneity is likely to be stochastic. Overall, we only observed a small fraction of alternative splicing and polyadenylation isoforms that are unlikely to be solely stochastic, implying the functional relevance of alternative splicing and polyadenylation in E. histolytica is limited. Lastly, we revised the gene models and annotated their 3′UTR in AmoebaDB, providing valuable resources to the community.

3' processing in protists

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

Serum-dependent selective expression of EhTMKB1-9, a member of Entamoeba histolytica B1 family of transmembrane kinases

Entamoeba histolytica transmembrane kinases (EhTMKs) can be grouped into six distinct families on the basis of motifs and sequences. Analysis of the E. histolytica genome revealed the presence of 35 EhTMKB1 members on the basis of sequence identity (≥95%). Only six homologs were full length containing an extracellular domain, a transmembrane segment and an intracellular kinase domain. Reverse transcription followed by polymerase chain reaction (RT-PCR) of the kinase domain was used to generate a library of expressed sequences. Sequencing of randomly picked clones from this library revealed that about 95% of the clones were identical with a single member, EhTMKB1-9, in proliferating cells. On serum starvation, the relative number of EhTMKB1-9 derived sequences decreased with concomitant increase in the sequences derived from another member, EhTMKB1-18. The change in their relative expression was quantified by real time PCR. Northern analysis and RNase protection assay were used to study the temporal nature of EhTMKB1-9 expression after serum replenishment of starved cells. The results showed that the expression of EhTMKB1-9 was sinusoidal. Specific transcriptional induction of EhTMKB1-9 upon serum replenishment was further confirmed by reporter gene (luciferase) expression and the upstream sequence responsible for serum responsiveness was identified. EhTMKB1-9 is one of the first examples of an inducible gene in Entamoeba. The protein encoded by this member was functionally characterized. The recombinant kinase domain of EhTMKB1-9 displayed protein kinase activity. It is likely to have dual specificity as judged from its sensitivity to different kinase inhibitors. Immuno-localization showed EhTMKB1-9 to be a surface protein which decreased on serum starvation and got relocalized on serum replenishment. Cell lines expressing either EhTMKB1-9 without kinase domain, or EhTMKB1-9 antisense RNA, showed decreased cellular proliferation and target cell killing. Our results suggest that E. histolytica TMKs of B1 family are functional kinases likely to be involved in serum response and cellular proliferation.

The presence of a vast array of putative transmembrane kinase genes suggests an extensive network of signaling systems in E. histolytica, particularly the ability to perceive signals from the extracellular environment and transduce these intracellularly. However, it has been very difficult to work with these molecules due to the presence of a large number of homologs. It is also not clear if these molecules are indeed protein kinases, as no kinase activity has yet been shown associated with these molecules. In this report, we show that EhTMKB1-9 is a protein kinase and it is one of the early serum-induced genes. It is a predominant EhTMKB1 molecule that is expressed in proliferating cells and its expression is modulated by serum. Cells containing a reduced level of EhTMKB1-9 or high level of a mutant protein result in decreased proliferation, target cell killing and adherence. The results presented in this report suggest that EhTMKB1-9 is an important signaling molecule likely to be involved in E. histolytica proliferation and virulence. We have also identified a serum starvation induced response where expression of EhTMKB1-18 was found to be induced.

Regulation of gene expression in protozoa parasites

Infections with protozoa parasites are associated with high burdens of morbidity and mortality across the developing world. Despite extensive efforts to control the transmission of these parasites, the spread of populations resistant to drugs and the lack of effective vaccines against them contribute to their persistence as major public health problems. Parasites should perform a strict control on the expression of genes involved in their pathogenicity, differentiation, immune evasion, or drug resistance, and the comprehension of the mechanisms implicated in that control could help to develop novel therapeutic strategies. However, until now these mechanisms are poorly understood in protozoa. Recent investigations into gene expression in protozoa parasites suggest that they possess many of the canonical machineries employed by higher eukaryotes for the control of gene expression at transcriptional, posttranscriptional, and epigenetic levels, but they also contain exclusive mechanisms. Here, we review the current understanding about the regulation of gene expression in Plasmodium sp., Trypanosomatids, Entamoeba histolytica and Trichomonas vaginalis.