Multi-omics analysis to characterize molecular adaptation of Entamoeba histolytica during serum stress.

Entamoeba histolytica is responsible for dysentery and extraintestinal disease in humans. To establish successful infection, it must generate adaptive response against stress due to host defense mechanisms. We have developed a robust proteomics workflow by combining miniaturized sample preparation, low flow-rate chromatography, and ultra-high sensitivity mass spectrometry, achieving increased proteome coverage, and further integrated proteomics and RNA-seq data to decipher regulation at translational and transcriptional levels. Label-free quantitative proteomics led to identification of 2344 proteins, an improvement over the maximum number identified in E. histolytica proteomic studies. In serum-starved cells, 127 proteins were differentially abundant and were associated with functions including antioxidant activity, cytoskeleton, translation, catalysis, and transport. The virulence factor, Gal/GalNAc-inhibitable lectin subunits, was significantly altered. Integration of transcriptomic and proteomic data revealed that only 30% genes were coordinately regulated at both transcriptional and translational levels. Some highly expressed transcripts did not change in protein abundance. Conversely, genes with no transcriptional change showed enhanced protein abundance, indicating post-transcriptional regulation. This multi-omics approach enables more refined gene expression analysis to understand the adaptive response of E. histolytica during growth stress.

The non-LTR retrotransposons of Entamoeba histolytica: genomic organization and biology.

Genome sequence analysis of Entamoeba species revealed various classes of transposable elements. While E. histolytica and E. dispar are rich in non-long terminal repeat (LTR) retrotransposons, E. invadens contains predominantly DNA transposons. Non-LTR retrotransposons of E. histolytica constitute three families of long interspersed nuclear elements (LINEs), and their short, nonautonomous partners, SINEs. They occupy ~ 11% of the genome. The EhLINE1/EhSINE1 family is the most abundant and best studied. EhLINE1 is 4.8 kb, with two ORFs that encode functions needed for retrotransposition. ORF1 codes for the nucleic acid-binding protein, and ORF2 has domains for reverse transcriptase (RT) and endonuclease (EN). Most copies of EhLINEs lack complete ORFs. ORF1p is expressed constitutively, but ORF2p is not detected. Retrotransposition could be demonstrated upon ectopic over expression of ORF2p, showing that retrotransposition machinery is functional. The newly retrotransposed sequences showed a high degree of recombination. In transcriptomic analysis, RNA-Seq reads were mapped to individual EhLINE1 copies. Although full-length copies were transcribed, no full-length 4.8 kb transcripts were seen. Rather, sense transcripts mapped to ORF1, RT and EN domains. Intriguingly, there was strong antisense transcription almost exclusively from the RT domain. These unique features of EhLINE1 could serve to attenuate retrotransposition in E. histolytica.

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.

Cytosine DNA methylation at promoter of non LTR retrotransposon and heat shock protein gene (HSP70) of Entamoeba histolytica and lack of correlation with transcription status.

Non LTR retrotransposons (EhLINEs and EhSINEs) occupy 11% of the Entamoeba histolytica genome. Since promoter DNA methylation at cytosines has been correlated with transcriptional silencing of transposable elements in model organisms we checked whether this was the case in EhLINE1. We located promoter activity in a 841bp fragment at 5'-end of this element by luciferase reporter assay. From RNAseq and RT-PCR analyses we selected a transcriptionally active and silent copy to study cytosine DNA methylation of the promoter region by bisulfite sequencing. None of the cytosines were methylated in either copy. Further, we looked at methylation status of a few selected cytosines in all 5'-intact EhLINE1 copies by single nucleotide incorporation opposite cytosine in bisulfite-treated DNA, where dGTP would be incorporated if the cytosine was methylated. Again we did not find evidence of cytosine methylation, indicating that expression status of this element was not correlated with promoter DNA methylation. To test for any role of cytosine methylation in transcriptional regulation of the E. histolytica Hsp70 gene in which the promoter is fully methylated under normal growth conditions, we checked methylation status and found that the promoter remained fully methylated during heat-shock as well, although transcription was greatly enhanced by heat-shock, showing that cytosine methylation is not a repressive mark for EhHsp70. Our data present direct evidence that promoter methylation, a common mode of transposon silencing, is unlikely to be involved in transcriptional regulation of EhLINE1, and reinforce the conclusion that promoter DNA methylation may not be a major contributor to transcriptional regulation in E. histolytica.

Functionally conserved RNA-binding and protein-protein interaction properties of LINE-ORF1p in an ancient clade of non-LTR retrotransposons of Entamoeba histolytica.

Retrotransposons are mobile genetic elements found in most organisms. Their origin and evolution is not very well understood. Retrotransposons that lack long terminal repeats (non-LTR) have been classified based on their reverse transcriptase (RT) and endonuclease sequences into groups, of which R2 is the most ancient. Its members contain a single open reading frame (ORF) while there are two ORFs in the other groups, of which ORF2 contains the RT and endonuclease sequences. It is thought that ORF1 was added later to the single-ORF-containing elements, and codes for a protein with nucleic acid binding activity. We have examined the non-LTR retrotransposons in Entamoeba histolytica, an early-branching parasitic protist, which belongs to the R2 group. However, unlike other members of R2, E. histolytica contains two ORFs. Here we show that EhLINE1-ORF1p is functionally related to the ORF1p found in the non-R2 groups. Its N-terminal region has RNA-binding activity and its C-terminal has a coiled coil domain which participates in protein-protein interaction. It lacks sequence-specificity of RNA-binding and binds to EhLINE1-RNA fragment and ribosomal RNA with comparable affinities. Our study suggests that ORF1p could have evolved independently to maintain functional conservation.

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.

Autophosphorylation of Ser428 of EhC2PK plays a critical role in regulating erythrophagocytosis in the parasite Entamoeba histolytica.

The protozoan parasite Entamoeba histolytica can invade both intestinal and extra intestinal tissues resulting in amoebiasis. During the process of invasion E. histolytica ingests red blood and host cells using phagocytic processes. Though phagocytosis is considered to be a key virulence determinant, the mechanism is not very well understood in E. histolytica. We have recently demonstrated that a novel C2 domain-containing protein kinase, EhC2PK is involved in the initiation of erythrophagocytosis. Because cells overexpressing the kinase-dead mutant of EhC2PK displayed a reduction in erythrophagocytosis, it appears that kinase activity is necessary for initiation. Biochemical analysis showed that EhC2PK is an unusual Mn(2+)-dependent serine kinase. It has a trans-autophosphorylated site at Ser(428) as revealed by mass spectrometric and biochemical analysis. The autophosphorylation defective mutants (S428A, KDΔC) showed a reduction in auto and substrate phosphorylation. Time kinetics of in vitro kinase activity suggested two phases, an initial short slow phase followed by a rapid phase for wild type protein, whereas mutations in the autophosphorylation sites that cause defect (S428A) or conferred phosphomimetic property (S428E) displayed no distinct phases, suggesting that autophosphorylation may be controlling kinase activity through an autocatalytic mechanism. A reduction and delay in erythrophagocytosis was observed in E. histolytica cells overexpressing S428A and KDΔC proteins. These results indicate that enrichment of EhC2PK at the site of phagocytosis enhances the rate of trans-autophosphorylation, thereby increasing kinase activity and regulating the initiation of erythrophagocytosis in E. histolytica.