Functional characterization of spliceosomal introns and identification of U2, U4, and U5 snRNAs in the deep-branching eukaryote Entamoeba histolytica

Proteomics approaches to understand cell biology and virulence of Entamoeba histolytica protozoan parasite

Entamoeba histolytica is the primitive eukaryotic parasite responsible of human amoebiasis, a disease characterized by bloody intestinal diarrhea and invasive extraintestinal illness. The knowledge of the complete genome sequence of virulent E. histolytica and related non-pathogenic species allowed the development of novel genome-wide methodological approaches including protein expression profiling and cellular proteomics in the so called post-genomic era. Proteomics studies have greatly increased our understanding of the cell biology of this ancient parasite. This review summarizes the current works concerning proteomics studies on cell biology, life cycle, virulence and pathogenesis, novel therapies, and protein expression regulation mechanisms in E. histolytica parasite. Also, we discuss the use of proteomics data for the development of novel therapies, the identification of potential disease biomarkers and differential diagnosis between species. SIGNIFICANCE: Entamoeba histolytica is the unicellular protozoan parasite responsible of human amoebiasis, a serious disease with worldwide distribution characterized by bloody intestinal diarrhea and invasive extraintestinal illness including peritonitis and liver, pulmonary and brain abscesses. The post-genomic era allowed the development of proteomic studies including protein expression profiling and cellular proteomics. These proteomics studies have greatly increased our understanding on cell biology, life cycle (cyst-trophozoite conversion), virulence, pathogenesis, novel therapies, and protein expression regulation mechanisms in E. histolytica. Importantly, proteomics has revealed the identity of proteins related to novel therapies, and the identification of potential disease biomarkers and proteins with use in diagnosis between species. Hopefully in the coming years, and through the use of more sophisticated omics tools, including deep proteomics, a more complete set of proteins involved in the aforementioned cellular processes can be obtained to understand the biology of this ancient eukaryote.

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.

The Entamoeba histolytica Syf1 Homolog Is Involved in the Splicing of AG-Dependent and AG-Independent Transcripts

Syf1 is a tetratricopeptide repeat (TPR) protein implicated in transcription elongation, spliceosome conformation, mRNA nuclear-cytoplasmic export and transcription-coupled DNA repair. Recently, we identified the spliceosomal components of the human parasite Entamoeba histolytica, among them is EhSyf. Molecular predictions confirmed that EhSyf contains 15 type 1 TPR tandem α-antiparallel array motifs. Amoeba transformants carrying plasmids overexpressing HA-tagged or EhSyf silencing plasmids were established to monitor the impact of EhSyf on the splicing of several test Entamoeba transcripts. EhSyf Entamoeba transformants efficiently silenced or overexpressed the proteins in the nucleus. The overexpression or absence of EhSyf notably enhanced or blocked splicing of transcripts irrespective of the strength of their 3′ splice site. Finally, the absence of EhSyf negatively affected the transcription of an intron-less transcript. Altogether our data suggest that EhSyf is a bona fide Syf1 ortholog involved in transcription and splicing.

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.

The Tudor Staphylococcal Nuclease Protein of Entamoeba histolytica Participates in Transcription Regulation and Stress Response

Entamoeba histolytica is the protozoa parasite responsible of human amoebiasis, disease that causes from 40,000 to 100,000 deaths annually worldwide. However, few are known about the expression regulation of molecules involved in its pathogenicity. Transcription of some virulence-related genes is positively controlled by the cis-regulatory element named URE1. Previously we identified the transcription factor that binds to URE1, which displayed a nuclear and cytoplasmic localization. This protein belongs to the Tudor Staphyococcal nuclease (TSN) family, which in other systems participates in virtually all pathways of gene expression, suggesting that this amoebic transcription factor (EhTSN; former EhURE1BP) could also play multiple functions in E. histolytica. The aim of this study was to identify the possible cellular events where EhTSN is involved. Here, we found that EhTSN in nucleus is located in euchromatin and close to, but not into, heterochromatin. We also showed the association of EhTSN with proteins involved in transcription and that the knockdown of EhTSN provokes a diminishing in the mRNA level of the EhRabB gene, which in its promoter region contains the URE1 motif, confirming that EhTSN participates in transcription regulation. In cytoplasm, this protein was found linked to the membrane of small vesicles and to plasma membrane. Through pull-down assays and mass spectrometry we identity thirty two candidate proteins to interact with EhTSN. These proteins participate in transcription, metabolism, signaling, and stress response, among other cellular processes. Interaction of EhTSN with some candidate proteins involved in metabolism, and signaling was validated by co-immunoprecipitation or co-localization. Finally we showed the co-localization of EhTSN and HSP70 in putative stress granules during heat shock and that the knockdown of EhTSN increases the cell death during heat shock treatment, reinforcing the hypothesis that EhTSN has a role during stress response. All data support the proposal that EhTSN is a multifunctional protein of E. histolytica.

Proteomic analysis of Entamoeba histolytica in vivo assembled pre-mRNA splicing complexes

The genome of the human intestinal parasite Entamoeba histolytica contains nearly 3000 introns and bioinformatic predictions indicate that major and minor spliceosomes occur in Entamoeba. However, except for the U2-, U4-, U5- and U6 snRNAs, no other splicing factor has been cloned and characterized. Here, we HA-tagged cloned the snRNP component U1A and assessed its expression and nuclear localization. Because the snRNP-free U1A form interacts with polyadenylate-binding protein, HA-U1A immunoprecipitates could identify early and late splicing complexes. Avoiding Entamoeba's endonucleases and ensuring the precipitation of RNA-binding proteins, parasite cultures were UV cross-linked prior to nuclear fraction immunoprecipitations with HA antibodies, and precipitates were subjected to tandem mass spectrometry (MS/MS) analyses. To discriminate their nuclear roles (chromatin-, co-transcriptional-, splicing-related), MS/MS analyses were carried out with proteins eluted with MS2-GST-sepharose from nuclear extracts of an MS2 aptamer-tagged Rabx13 intron amoeba transformant. Thus, we probed thirty-six Entamoeba proteins corresponding to 32 cognate splicing-specific factors, including 13 DExH/D helicases required for all stages of splicing, and 12 different splicing-related helicases were identified also. Furthermore 50 additional proteins, possibly involved in co-transcriptional processes were identified, revealing the complexity of co-transcriptional splicing in Entamoeba. Some of these later factors were not previously found in splicing complex analyses.

BIOLOGICAL SIGNIFICANCE

Numerous facts about the splicing of the nearly 3000 introns of the Entamoeba genome have not been unraveled, particularly the splicing factors and their activities. Considering that many of such introns are located in metabolic genes, the knowledge of the splicing cues has the potential to be used to attack or control the parasite. We have found numerous new splicing-related factors which could have therapeutic benefit. We also detected all the DExH/A RNA helicases involved in splicing and splicing proofreading control. Still, Entamoeba is very inefficient in splicing fidelity, thus we may have found a possible model system to study these processes.

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

Background

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

Results

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

Conclusions

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

A phylogenetic study of Drosophila splicing assembly chaperone RNP-4F associated U4-/U6-snRNA secondary structure

The rnp-4f gene in Drosophila melanogaster encodes nuclear protein RNP-4F. This encoded protein is represented by homologs in other eukaryotic species, where it has been shown to function as an intron splicing assembly factor. Here, RNP-4F is believed to initially bind to a recognition sequence on U6-snRNA, serving as a chaperone to facilitate its association with U4-snRNA by intermolecular hydrogen bonding. RNA conformations are a key factor in spliceosome function, so that elucidation of changing secondary structures for interacting snRNAs is a subject of considerable interest and importance. Among the five snRNAs which participate in removal of spliceosomal introns, there is a growing consensus that U6-snRNA is the most structurally dynamic and may constitute the catalytic core. Previous studies by others have generated potential secondary structures for free U4- and U6-snRNAs, including the Y-shaped U4-/U6-snRNA model. These models were based on study of RNAs from relatively few species, and the popular Y-shaped model remains to be systematically re-examined with reference to the many new sequences generated by recent genomic sequencing projects. We have utilized a comparative phylogenetic approach on 60 diverse eukaryotic species, which resulted in a revised and improved U4-/U6-snRNA secondary structure. This general model is supported by observation of abundant compensatory base mutations in every stem, and incorporates more of the nucleotides into base-paired associations than in previous models, thus being more energetically stable. We have extensively sampled the eukaryotic phylogenetic tree to its deepest roots, but did not find genes potentially encoding either U4- or U6-snRNA in the Giardia and Trichomonas data-bases. Our results support the hypothesis that nuclear introns in these most deeply rooted eukaryotes may represent evolutionary intermediates, sharing characteristics of both group II and spliceosomal introns. An unexpected result of this study was discovery of a potential competitive binding site for Drosophila splicing assembly factor RNP-4F to a 5’-UTR regulatory region within its own premRNA, which may play a role in negative feedback control.

Targets of the Entamoeba histolytica transcription factor URE3-BP

The Entamoeba histolytica transcription factor Upstream Regulatory Element 3-Binding Protein (URE3-BP) is a calcium-responsive regulator of two E. histolytica virulence genes, hgl5 and fdx1. URE3-BP was previously identified by a yeast one-hybrid screen of E. histolytica proteins capable of binding to the sequence TATTCTATT (Upstream Regulatory Element 3 (URE3)) in the promoter regions of hgl5 and fdx1. In this work, precise definition of the consensus URE3 element was performed by electrophoretic mobility shift assays (EMSA) using base-substituted oligonucleotides, and the consensus motif validated using episomal reporter constructs. Transcriptome profiling of a strain induced to produce a dominant-positive URE3-BP was then used to identify additional genes regulated by URE3-BP. Fifty modulated transcripts were identified, and of these the EMSA defined motif T[atg]T[tc][cg]T[at][tgc][tg] was found in over half of the promoters (54% p<0.0001). Fifteen of the URE3-BP regulated genes were potential membrane proteins, suggesting that one function of URE3-BP is to remodel the surface of E. histolytica in response to a calcium signal. Induction of URE3-BP leads to an increase in tranwell migration, suggesting a possible role in the regulation of cellular motility.

Spliceosomal snRNAs in the unicellular eukaryote Trichomonas vaginalis are structurally conserved but lack a 5'-cap structure

Few genes in the divergent eukaryote Trichomonas vaginalis have introns, despite the unusually large gene repertoire of this human-infective parasite. These introns are characterized by extended conserved regulatory motifs at the 5' and 3' boundaries, a feature shared with another divergent eukaryote, Giardia lamblia, but not with metazoan introns. This unusual characteristic of T. vaginalis introns led us to examine spliceosomal small nuclear RNAs (snRNAs) predicted to mediate splicing reactions via interaction with intron motifs. Here we identify T. vaginalis U1, U2, U4, U5, and U6 snRNAs, present predictions of their secondary structures, and provide evidence for interaction between the U2/U6 snRNA complex and a T. vaginalis intron. Structural models predict that T. vaginalis snRNAs contain conserved sequences and motifs similar to those found in other examined eukaryotes. These data indicate that mechanisms of intron recognition as well as coordination of the two catalytic steps of splicing have been conserved throughout eukaryotic evolution. Unexpectedly, we found that T. vaginalis spliceosomal snRNAs lack the 5' trimethylguanosine cap typical of snRNAs and appear to possess unmodified 5' ends. Despite the lack of a cap structure, U1, U2, U4, and U5 genes are transcribed by RNA polymerase II, whereas the U6 gene is transcribed by RNA polymerase III.

Computational screen for spliceosomal RNA genes aids in defining the phylogenetic distribution of major and minor spliceosomal components

The RNA molecules of the spliceosome are critical for specificity and catalysis during splicing of eukaryotic pre-mRNA. In order to examine the evolution and phylogenetic distribution of these RNAs, we analyzed 149 eukaryotic genomes representing a broad range of phylogenetic groups. RNAs were predicted using high-sensitivity local alignment methods and profile HMMs in combination with covariance models. The results provide the most comprehensive view so far of the phylogenetic distribution of spliceosomal RNAs. RNAs were predicted in many phylogenetic groups where these RNA were not previously reported. Examples are RNAs of the major (U2-type) spliceosome in all fungal lineages, in lower metazoa and many protozoa. We also identified the minor (U12-type) spliceosomal U11 and U6atac RNAs in Acanthamoeba castellanii, where U12 spliceosomal RNA as well as minor introns were reported recently. In addition, minor-spliceosome-specific RNAs were identified in a number of phylogenetic groups where previously such RNAs were not observed, including the nematode Trichinella spiralis, the slime mold Physarum polycephalum and the fungal lineages Zygomycota and Chytridiomycota. The detailed map of the distribution of the U12-type RNA genes supports an early origin of the minor spliceosome and points to a number of occasions during evolution where it was lost.

Implications of the new eukaryotic systematics for parasitologists

An accurate understanding of evolutionary relationships is central in biology. For parasitologists, understanding the relationships among eukaryotic organisms allows the prediction of virulence mechanisms, reconstruction of metabolic pathways, identification of potential drug targets, elucidation of parasite-specific cellular processes and understanding of interactions with the host or vector. Here we consider the impact of major recent revisions of eukaryotic systematics and taxonomy on parasitology. The previous, ladder-like model placed some protists as early diverging, with the remaining eukaryotes "progressing" towards a "crown radiation" of animals, plants, Fungi and some additional protistan lineages. This model has been robustly disproven. The new model is based on vastly increased amounts of molecular sequence data, integration with morphological information and the rigorous application of phylogenetic methods to those data. It now divides eukaryotes into six major supergroups; the relationships between those groups and the order of branching remain unknown. This new eukaryotic phylogeny emphasizes that organisms including Giardia, Trypanosoma and Trichomonas are not primitive, but instead highly evolved and specialised for their specific environments. The wealth of newly available comparative genomic data has also allowed the reconstruction of ancient suites of characteristics and mapping of character evolution in diverse parasites. For example, the last common eukaryotic ancestor was apparently complex, suggesting that lineage-specific adaptations and secondary losses have been important in the evolution of protistan parasites. Referring to the best evidence-based models for eukaryotic evolution will allow parasitologists to make more accurate and reliable inferences about pathogens that cause significant morbidity and mortality.

Functional and evolutionary analysis of alternatively spliced genes is consistent with an early eukaryotic origin of alternative splicing

Background

Alternative splicing has been reported in various eukaryotic groups including plants, apicomplexans, diatoms, amoebae, animals and fungi. However, whether widespread alternative splicing has evolved independently in the different eukaryotic groups or was inherited from their last common ancestor, and may therefore predate multicellularity, is still unknown. To better understand the origin and evolution of alternative splicing and its usage in diverse organisms, we studied alternative splicing in 12 eukaryotic species, comparing rates of alternative splicing across genes of different functional classes, cellular locations, intron/exon structures and evolutionary origins.

Results

For each species, we find that genes from most functional categories are alternatively spliced. Ancient genes (shared between animals, fungi and plants) show high levels of alternative splicing. Genes with products expressed in the nucleus or plasma membrane are generally more alternatively spliced while those expressed in extracellular location show less alternative splicing. We find a clear correspondence between incidence of alternative splicing and intron number per gene both within and between genomes. In general, we find several similarities in patterns of alternative splicing across these diverse eukaryotes.

Conclusion

Along with previous studies indicating intron-rich genes with weak intron boundary consensus and complex spliceosomes in ancestral organisms, our results suggest that at least a simple form of alternative splicing may already have been present in the unicellular ancestor of plants, fungi and animals. A role for alternative splicing in the evolution of multicellularity then would largely have arisen by co-opting the preexisting process.