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.

Entamoeba histolytica and pathogenesis: A calcium connection.

Calcium signaling plays a key role in many essential processes in almost all eukaryotic systems. It is believed that it may also be an important signaling system of the protist parasite Entamoeba histolytica. Motility, adhesion, cytolysis, and phagocytosis/trogocytosis are important steps in invasion and pathogenesis of E. histolytica, and Ca2+ signaling is thought to be associated with these processes leading to tissue invasion. There are a large number of Ca2+-binding proteins (CaBPs) in E. histolytica, and a number of these proteins appear to be associated with different steps in pathogenesis. The genome encodes 27 EF-hand-containing CaBPs in addition to a number of other Ca2+-binding domain/motif-containing proteins, which suggest intricate calcium signaling network in this parasite. Unlike other eukaryotes, a typical calmodulin-like protein has not been seen in E. histolytica. Though none of the CaBPs display sequence similarity with a typical calmodulin, extensive structural similarity has been seen in spite of lack of significant functional overlap with that of typical calmodulins. One of the unique features observed in E. histolytica is the identification of CaBPs (EhCaBP1, EhCaBP3) that have the ability to directly bind actin and modulate actin dynamics. Direct interaction of CaBPs with actin has not been seen in any other system. Pseudopod formation and phagocytosis are some of the processes that require actin dynamics, and some of the amoebic CaBPs (EhC2Pk, EhCaBP1, EhCaBP3, EhCaBP5) participate in this process. None of these E. histolytica CaBPs have any homolog in organisms other than different species of Entamoeba, suggesting a novel Ca2+ signaling pathway that has evolved in this genus.

Tissue-specific isoform expression of GNE gene in human tissues.

Mutations in the sialic acid biosynthesis enzyme GNE lead to a late-onset, debilitating neuromuscular disorder, GNE myopathy, characterized by progressive skeletal muscle weakness. The mechanisms responsible for skeletal muscle specificity, late-onset, and disease progression are unknown. Our main aim is to understand the reason for skeletal muscle-specific phenotype. To answer this question, we have analyzed the expression profile of the GNE gene and its multiple mRNA variants in different human tissues. A combinatorial approach encompassing bioinformatics tools and molecular biology techniques was used. NCBI, Ensembl, and GTEx were used for data mining. The expression analysis of GNE and its variants was performed with cDNA tissue panel using PCR and targeted RNA-seq. Among nine different GNE isoforms reported in this study, transcript variants 1, X1, and X2 were not tissue specific. Transcript variants 1, 6, X1, and X2, were found in skeletal muscles suggesting their possible role in GNE myopathy. In the current study, we present new data about GNE expression patterns in human tissues. Our results suggest that there may be a link between tissue-specific pathology and isoform pattern in skeletal muscles, which could provide clues for the development of new treatment strategies for GNE myopathy.

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.

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.

Tissue specific expression of sialic acid metabolic pathway: role in GNE myopathy.

GNE myopathy is an adult-onset degenerative muscle disease that leads to extreme disability in patients. Biallelic mutations in the rate-limiting enzyme UDP-N-acetylglucosamine-2-epimerase/N-acetylmannosamine-kinase (GNE) of sialic acid (SA) biosynthetic pathway, was shown to be the cause of this disease. Other genetic disorders with muscle pathology where defects in glycosylation are known. It is yet not clear why a defect in SA biosynthesis and glycosylation affect muscle cells selectively even though they are ubiquitously present in all tissues. Here we have comprehensively examined the complete SA metabolic pathway involving biosynthesis, sialylation, salvage, and catabolism. To understand the reason for tissue-specific phenotype caused by mutations in genes of this pathway, we analysed the expression of different SA pathway genes in various tissues, during the muscle tissue development and in muscle tissues from GNE myopathy patients (p.Met743Thr) using publicly available databases. We have also analysed gene co-expression networks with GNE in different tissues as well as gene interactions that are unique to muscle tissues only. The results do show a few muscle specific interactions involving ANLN, MYO16 and PRAMEF25 that could be involved in specific phenotype. Overall, our results suggest that SA biosynthetic and catabolic genes are expressed at a very low level in skeletal muscles that also display a unique gene interaction network.

Myg1 exonuclease couples the nuclear and mitochondrial translational programs through RNA processing.

Semi-autonomous functioning of mitochondria in eukaryotic cell necessitates coordination with nucleus. Several RNA species fine-tune mitochondrial processes by synchronizing with the nuclear program, however the involved components remain enigmatic. In this study, we identify a widely conserved dually localized protein Myg1, and establish its role as a 3'-5' RNA exonuclease. We employ mouse melanoma cells, and knockout of the Myg1 ortholog in Saccharomyces cerevisiae with complementation using human Myg1 to decipher the conserved role of Myg1 in selective RNA processing. Localization of Myg1 to nucleolus and mitochondrial matrix was studied through imaging and confirmed by sub-cellular fractionation studies. We developed Silexoseqencing, a methodology to map the RNAse trail at single-nucleotide resolution, and identified in situ cleavage by Myg1 on specific transcripts in the two organelles. In nucleolus, Myg1 processes pre-ribosomal RNA involved in ribosome assembly and alters cytoplasmic translation. In mitochondrial matrix, Myg1 processes 3'-termini of the mito-ribosomal and messenger RNAs and controls translation of mitochondrial proteins. We provide a molecular link to the possible involvement of Myg1 in chronic depigmenting disorder vitiligo. Our study identifies a key component involved in regulating spatially segregated organellar RNA processing and establishes the evolutionarily conserved ribonuclease as a coordinator of nucleo-mitochondrial crosstalk.

Novel regulatory roles of PtdIns(4,5)P2 generating enzyme EhPIPKI in actin dynamics and phagocytosis of Entamoeba histolytica.

Motility and phagocytosis are the two important processes that are intricately linked to survival and virulence potential of the protist parasite Entamoeba histolytica. These processes primarily rely on actin-dependent pathways, and regulation of these pathways is critical for understanding the pathology of E. histolytica. Generally, phosphoinositides dynamics have not been explored in amoebic actin dynamics and particularly during phagocytosis in E. histolytica. We have explored the roles of PtdIns(4,5)P2 as well as the enzyme that produces this metabolite, EhPIPKI during phagocytosis. Immunofluorescence and live cell images showed enrichment of EhPIPKI in different stages of phagocytosis from initiation till the cups progressed towards closure. However, the enzyme was absent after phagosomes are pinched off from the membrane. Overexpression of a dominant negative mutant revealed a reduction in the formation of phagocytic cups and inhibition in the rate of engulfment of erythrocytes. Moreover, EhPIPKI binds directly to F and G-actin unlike PIPKs from other organisms. PtdIns(4,5)P2 , the product of the enzyme, also followed a similar distribution pattern during phagocytosis as determined by a GFP-tagged PH-domain from PLCδ, which specifically binds PtdIns(4,5)P2 in trophozoites. In summary, EhPIPKI regulates initiation of phagocytosis by regulating actin dynamics.

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.

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.

Calcium-binding protein EhCaBP3 is recruited to the phagocytic complex of Entamoeba histolytica by interacting with Arp2/3 complex subunit 2.

Phagocytosis is involved in invasive disease of the parasite Entamoeba histolytica. Upon binding of red blood cells, there is a sequential recruitment of EhC2PK, EhCaBP1, EhAK1, and Arp2/3 complex during the initiation phase. In addition, EhCaBP3 is also recruited to the site and, along with myosin 1B, is thought to be involved in progression of phagocytic cups from initiation to phagosome formation. However, it is not clear how EhCaBP3 gets recruited to the rest of the phagocytic machinery. Here, we show that EhARPC2, a subunit of Arp2/3 complex, interacts with EhCaBP3 in a Ca2+ -dependent manner both in vivo and in vitro. Imaging and pull down experiments suggest that interaction with EhARPC2 is required for the closure of cups and formation of phagosomes. Moreover, downregulation of EhARPC2 prevents localisation of EhCaBP3 to phagocytic cups, suggesting that EhCaBP3 is part of EhC2PK-EhCaBP1-EhAK1-Arp2/3 complex (EhARPC1) pathway. In conclusion, these results suggest that the EhCaBP3-EhARPC2 interaction helps to recruit EhCaBP3 along with myosin 1B to the phagocytic machinery that plays an indispensable role in E. histolytica phagocytosis.

EhRho1 regulates phagocytosis by modulating actin dynamics through EhFormin1 and EhProfilin1 in Entamoeba histolytica.

The protist parasite Entamoeba histolytica causes amoebiasis, a major public health problem in developing countries and a major cause of morbidity and mortality. Invasive infection in amoebiasis mostly affects intestinal epithelial cell lining but can also involve other organs, such as liver, lungs, or brain. Phagocytosis is an essential mode of nutrition in amoeba and has often been associated with virulence behaviour of E. histolytica. E. histolytica possesses a highly dynamic and actin-rich cytoskeleton that is thought to be involved in many processes, such as motility, pseudopod formation, and pathogenesis. Rho GTPases are known to be key regulators of the actin cytoskeleton and consequently influence the shape and movement of cells. Our study is mainly focused to understand the role of EhRho1 in the phagocytosis process of E. histolytica. EhRho1 got enriched in the phagocytic cups along with EhActin and remains attached with phagosomal membrane. However, there was no direct binding of EhRho1 with G- or F-actin, though binding was observed with the actin nucleating proteins EhFormin1 and EhProfilin1. Overexpression of dominant negative mutant or lowering the expression by antisense RNA of EhRho1 in trophozoites caused delocalisation of EhFormin1 and EhProfilin1 from phagocytic cups, which results in impairment of phagocytic process and decrease in F-actin content. The overall results show that EhRho1 regulates phagocytosis by modulating actin dynamics through recruitment of EhFormin1 and EhProfilin1 at the phagocytosis nucleation site in E. histolytica.

EhRho1 regulates plasma membrane blebbing through PI3 kinase in Entamoeba histolytica.

The protozoan parasite Entamoeba histolytica causes amoebiasis, a major public health problem in developing countries. Motility of E. histolytica is important for its pathogenesis. Blebbing is an essential process contributing to cellular motility in many systems. In mammalian cells, formation of plasma membrane blebs is regulated by Rho-GTPases through its effectors, such as Rho kinase, mDia1, and acto-myosin proteins. In this study, we have illuminated the role of EhRho1 in bleb formation and motility of E. histolytica. EhRho1 was found at the site of bleb formation in plasma membrane of trophozoites. Overexpression of mutant EhRho1 defective for Guanosine triphosphate (GTP)-binding or down-regulating EhRho1 by antisense RNA resulted in reduced blebbing and motility. Moreover, serum-starvation reduced blebbing that was restored on serum-replenishment. Lysophosphatidic acid treatment induced bleb formation, whereas wortmannin inhibited the process. In all these cases, concentration of GTP-EhRho1 (active) and Phosphatidylinositol 4,5-bisphosphate (PIP2) inversely correlated with the level of plasma membrane blebbing. Our study suggests the role of EhRho1 in blebbing and bleb-based motility through PI3 kinase pathway in E. histolytica.

The Entamoeba histolytica, Arp2/3 Complex Is Recruited to Phagocytic Cups through an Atypical Kinase EhAK1.

The parasite Entamoeba histolytica is the etiological agent of amoebiasis and phagocytosis plays a key role in virulence of this organism. Signaling pathways involved in activation of cytoskeletal dynamics required for phagocytosis remain to be elucidated. Phagocytosis is initiated with sequential recruitment of EhC2PK, EhCaBP1, EhCaBP3 and an atypical kinase EhAK1 after particle attachment. Here we show that EhARPC1, an essential subunit of the actin branching complex Arp 2/3 is recruited to the phagocytic initiation sites by EhAK1. Imaging, expression knockdown of different molecules and pull down experiments suggest that EhARPC1 interacts with EhAK1 and that it is required during initiation of phagocytosis and phagosome formation. Moreover, recruitment of EhARPC2 at the phagocytosis initiation by EhAK1 is also observed, indicating that the Arp 2/3 complex is recruited. In conclusion, these results suggests a novel mechanism of recruitment of Arp 2/3 complex during phagocytosis in E. histolytica.

SINE polymorphism reveals distinct strains of Entamoeba histolytica from North India.

Entamoeba histolytica is an intestinal parasite causing significant morbidity and mortality in the developing world. More tools are needed to understand the epidemiology and molecular pathogenesis of amebiasis. Virulence pattern of E. histolytica could be linked with the genotype of a strain. Several loci showing insertion polymorphism of retrotransposable short interspersed nuclear elements EhSINE1 and EhSINE2 have been reported among laboratory strains of E. histolytica. The present study was undertaken to validate this observation in clinical isolates from north India. Our results indicate that the Indian samples show a different propensity of SINE retention or loss at two of the polymorphic loci compared with non-Indian samples. Statistical analysis of different loci revealed Locus 17 of EhSINE1as a potential geographical marker for distinguishing Indian isolates from non Indian isolates.

A novel alpha kinase EhAK1 phosphorylates actin and regulates phagocytosis in Entamoeba histolytica.

Phagocytosis plays a key role in nutrient uptake and virulence of the protist parasite Entamoeba histolytica. Phagosomes have been characterized by proteomics, and their maturation in the cells has been studied. However, there is so far not much understanding about initiation of phagocytosis and formation of phagosomes at the molecular level. Our group has been studying initiation of phagocytosis and formation of phagosomes in E. histolytica, and have described some of the molecules that play key roles in the process. Here we show the involvement of EhAK1, an alpha kinase and a SH3 domain containing protein in the pathway that leads to formation of phagosomes using red blood cell as ligand particle. A number of approaches, such as proteomics, biochemical, confocal imaging using specific antibodies or GFP tagged molecules, expression down regulation by antisense RNA, over expression of wild type and mutant proteins, were used to understand the role of EhAK1 in phagocytosis. EhAK1 was found in the phagocytic cups during the progression of cups, until closure of phagosomes, but not in the phagosomes themselves. It is recruited to the phagosomes through interaction with the calcium binding protein EhCaBP1. A reduction in phagocytosis was observed when EhAK1 was down regulated by antisense RNA, or by over expression of the kinase dead mutant. G-actin was identified as one of the major substrates of EhAK1. Phosphorylated actin preferentially accumulated at the phagocytic cups and over expression of a phosphorylation defective actin led to defects in phagocytosis. In conclusion, we describe an important component of the pathway that is initiated on attachment of red blood cells to E. histolytica cells. The main function of EhAK1 is to couple signalling events initiated after accumulation of EhC2PK to actin dynamics.

The Calmodulin-like calcium binding protein EhCaBP3 of Entamoeba histolytica regulates phagocytosis and is involved in actin dynamics.

Phagocytosis is required for proliferation and pathogenesis of Entamoeba histolytica and erythrophagocytosis is considered to be a marker of invasive amoebiasis. Ca²âº has been found to play a central role in the process of phagocytosis. However, the molecular mechanisms and the signalling mediated by Ca²âº still remain largely unknown. Here we show that Calmodulin-like calcium binding protein EhCaBP3 of E. histolytica is directly involved in disease pathomechanism by its capacity to participate in cytoskeleton dynamics and scission machinery during erythrophagocytosis. Using imaging techniques EhCaBP3 was found in phagocytic cups and newly formed phagosomes along with actin and myosin IB. In vitro studies confirmed that EhCaBP3 directly binds actin, and affected both its polymerization and bundling activity. Moreover, it also binds myosin 1B in the presence of Ca²âº. In cells where EhCaBP3 expression was down regulated by antisense RNA, the level of RBC uptake was reduced, myosin IB was found to be absent at the site of pseudopod cup closure and the time taken for phagocytosis increased, suggesting that EhCaBP3 along with myosin 1B mediate the closure of phagocytic cups. Experiments with EhCaBP3 mutant defective in Ca²âº-binding showed that Ca²âº binding is required for phagosome formation. Liposome binding assay revealed that EhCaBP3 recruitment and enrichment to membrane is independent of any cellular protein as it binds directly to phosphatidylserine. Taken together, our results suggest a novel pathway mediating phagocytosis in E. histolytica, and an unusual mechanism of modulation of cytoskeleton dynamics by two calcium binding proteins, EhCaBP1 and EhCaBP3 with mostly non-overlapping functions.