Search

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.

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.

For the protist parasite Entamoeba histolytica, endocytic processes, such as phagocytosis, are essential for its survival in the human gut. The actin cytoskeleton is involved in the formation of pseudopods and phagosomal vesicles by incorporating a number of actin-binding and modulating proteins along with actin in a temporal manner. The actin dynamics, which comprises polymerization, branching, and depolymerization is very tightly regulated and takes place directionally at the sites of initiation of phagocytosis. Formin and profilin are two actin-binding proteins that are known to regulate actin cytoskeleton dynamics and thereby, endocytic processes. In this article, we report the participation of formin and profilin in E. histolytica phagocytosis and propose that these two proteins interact with each other and their sequential recruitment at the site is required for the successful completion of phagocytosis. The evidence is based on detailed microscopic, live imaging, interaction studies, and expression downregulation. The cells downregulated for expression of formin show absence of profilin at the site of phagocytosis, whereas downregulation of profilin does not affect formin localization.

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.

During amoebiasis, colonization of the gut by Entamoeba histolytica can lead to alterations of the host microbiota. In this study we have compared the gut microbiota of patients of amoebiasis with healthy controls using 16S rRNA gene variable regions, (V1-V3) and (V3-V5) of the bacterial genome. From this 16S rRNA gene amplicon data one paired end and two single end data sets were selected and compared by the number of OTUs obtained, sequence count and diversity analysis. Our results showed that the V1-V3- paired end data set gave the maximum number of OTUs in comparison to the two single end data sets studied. The amoebiasis patient samples showed a significant drop in richness in the alpha diversity measurements and lower intra group similarity compared to the healthy controls. Bacteria of genus Prevotella, Sutterella and Collinsella were more abundant in healthy controls whereas Escherichia, Klebsiella and Ruminococcus were more abundant in the E. histolytica positive patients. All the healthy controls harboured bacteria belonging to Faecalibacterium, Prevotella, Ruminococcus, Subdoligranulum and Escherichia genera while all the E. histolytica positive patients samples contained genus Enterobacter.

GTPases are molecular switches, which regulate a variety of cellular processes such as cell polarity, gene transcription, microtubule dynamics, cell-cycle etc. In this paper, we characterize a Ca2+-binding protein from Entamoeba histolytica (EhCaBP6) as a novel GTPase. We locate the active site for GTP hydrolysis within the C-terminal domain of EhCaBP6, although it requires full length protein for its complete range of activity. Using NMR studies, we observe that GTP binding induces conformational change in EhCaBP6. The identification of this novel and unusual Ca2+-dependent GTPase is important to elucidate the unconventional cell cycle of E. histolytica.

Ca2+-binding proteins are present in almost all living organisms and different types display different levels of binding affinities for the cation. Here, we applied two new scoring schemes enabling the user to manipulate the binding affinities of such proteins. We specifically designed a unique EF-hand loop capable of binding calcium with high affinity by altering five residues of the loop based on the scoring scheme. We worked on the N-terminal domain of Entamoeba histolytica calcium-binding protein1 (NtEhCaBP1), and used site-directed mutagenesis to incorporate the designed loop sequence into the second EF hand motif of this protein. The binding isotherms calculated using ITC calorimetry showed a ∼500-fold greater association constant (Ka) for the mutant. The crystal structure of the mutant was also determined, and displayed more compact Ca2+-coordination spheres in both of its EF loops than did the structure of the wildtype protein, consistent with the greater calcium-binding affinities of the mutant. The NtEhCaBP1 mutant was also shown to form a hexamer rather than just a trimer, and this hexamer formation was attributed to the position of the last helix of the mutant having been changed as a result of the strong calcium coordination. Further dynamic correlation analysis revealed that the mutation in the second EF loop changed the entire residue network of the monomer, resulting in a stronger coordination of Ca2+ even in other EF hand loop.

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.

Entamoeba histolytica is an intestinal protist parasite that causes amoebiasis, a major source of morbidity and mortality in developing countries. Phosphoinositides are involved in signalling systems that have a role in invasion and pathogenesis of this parasite. Phosphatidylinositol 4-phosphate 5-kinase (PIP5K) catalyses the generation of phosphatidylinositol(4,5)bisphosphate (PtdIns(4,5)P2 ), a key species of phosphoinositide that regulates various cellular processes. However, phosphatidylinositol phosphate kinase (PIPK) family of enzymes have not been characterized in E. histolytica. Here, we report the identification and characterization of type I PIPK (EhPIPKI) of E. histolytica. Computational analysis revealed homologs of type I and III PIPK family in E. histolytica and the absence of type II PIPK. In spite of low overall sequence identity, the kinase domain was found to be highly conserved. Interestingly, a unique insertion of a tandem repeat motif was observed in EhPIPKI distinguishing it from existing PIPKs of other organisms. Substrate profiling showed that EhPIPKI could phosphorylate at third and fifth hydroxyl positions of phosphatidylinositols, though the predominant substrate was phosphatidylinositol 4-phosphate (PtdIns(4)P). Furthermore, EhPIPKI underwent intracellular cleavage close to the amino-terminal, generating two distinct fragments Nter-EhPIPKI (27p) and Cter-EhPIPKI (47p). Immunofluorescence and cellular fractionation revealed that the full-length EhPIPKI and the Cter-EhPIPKI containing carboxyl-terminal activation loop were present in the plasma membrane while the Nter-EhPIPKI was observed in the cytosolic region. In conclusion, E. histolytica has a single EhPIPKI gene that displays novel properties of post-translational processing, the presence of a repeat domain and substrate specificity not observed in any PIPK enzyme so far.

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.

SARS-CoV-2 is a member of the Coronavirus family which recently originated from the Wuhan province of China and spread very rapidly through the world infecting more than 4 million people. In the past, other Coronaviruses have also been found to cause human infection, but not as widespread as COVID-19. Since Coronavirus sequences constantly change due to mutation and recombination, it is important to understand the pattern of changes and likely path the virus can take in the future. In this study, we have used the Shewhart control chart to identify and analyze hypervariable (hotspots) and hypovariable (coldspots) regions of the virus. Our analysis shows that SARS-CoV-2 has changed in a few regions of the genome. Analysis of SARS-CoV-1 and MERS sequences suggests that over time, mutations start accumulating in different regions and most likely SARS-CoV-2 may also follow a similar path. The results suggest a possible emergence of modified viruses over some time. Electronic supplementary material The online version of this article (10.1007/s12038-020-00131-5) contains supplementary material, which is available to authorized users.

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.

The actin cytoskeleton, a dynamic network of actin filaments and associated F-actin–binding proteins, is fundamentally important in eukaryotes. α-Actinins are major F-actin bundlers that are inhibited by Ca(2+) in nonmuscle cells. Here we report the mechanism of Ca(2+)-mediated regulation of Entamoeba histolytica α-actinin-2 (EhActn2) with features expected for the common ancestor of Entamoeba and higher eukaryotic α-actinins. Crystal structures of Ca(2+)-free and Ca(2+)-bound EhActn2 reveal a calmodulin-like domain (CaMD) uniquely inserted within the rod domain. Integrative studies reveal an exceptionally high affinity of the EhActn2 CaMD for Ca(2+), binding of which can only be regulated in the presence of physiological concentrations of Mg(2+). Ca(2+) binding triggers an increase in protein multidomain rigidity, reducing conformational flexibility of F-actin–binding domains via interdomain cross-talk and consequently inhibiting F-actin bundling. In vivo studies uncover that EhActn2 plays an important role in phagocytic cup formation and might constitute a new drug target for amoebic dysentery.

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.

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.

GTPases are molecular switches, which regulate a variety of cellular processes such as cell polarity, gene transcription, microtubule dynamics, cell-cycle etc. In this paper, we characterize a Ca

Chronic myelocytic leukemia cell line K562 undergoes differentiation by phorbol esters to megakaryocytes and we have used this system to understand miRNA processing leading to isomiR generation. PMA treatment significantly altered the production of miRNA in K562 cells. Expression of 24.4% of miRNAs were found to be stimulated whereas expression of 10% miRNAs were inhibited by PMA treatment. Our results suggest that miRNA precursors are processed into isomiRs in a deterministic manner. The relative levels of different isomiRs of a miRNA remained mainly unchanged even after PMA treatment irrespective of overall changes in expression (either up-regulation or down-regulation). However, not all miRNAs behave in the same way, about 7% showed a variation of isomiR profiles after PMA treatment. Most of the later class of miRNAs were found to be oncogenic miRNAs. Further, it was also found that number of isomiRs was independent of abundance of a miRNA. Functional importance of different isomiRs was demonstrated using three different isomiRs of miR-22. Our results showed that different isomiRs could inhibit expression of targets genes with different efficiencies. Our study suggests that the heterogeneity of a miRNA population generated during processing is in general regulated and that variation in the generation of an isomiR can be a functionally important regulatory feature.

Selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) is a versatile sequence independent method to probe RNA structure in vivo and in vitro. It has so far been tried mainly with model organisms. We show that cells of Entamoeba histolytica, a protozoan parasite of humans are hyper-sensitive to the in vivo SHAPE reagent, NAI, and show rapid loss of viability and RNA integrity. We optimized treatment conditions with 5.8S rRNA and Eh_U3 snoRNA to obtain NAI-modification while retaining RNA integrity. The modification patterns were highly reproducible. The in vivo folding was different from in vitro and correlated well with known interactions of 5.8S rRNA with proteins in vivo. The Eh_U3 snoRNA also showed many differences in its in vivo versus in vitro folding, which correlated with conserved interactions of this RNA with 18S rRNA and 5′-ETS. Further, Eh_U3 snoRNA obtained from serum-starved cells showed an open 3′-hinge structure, indicating disruption of 5′-ETS interaction. This could contribute to the observed slow processing of pre-rRNA in starved cells. Our work shows the applicability of SHAPE to study in vivo RNA folding in a parasite and will encourage the use of this reagent for RNA structure analysis in other such organisms.

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.

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.

We report almost complete sequence specific 1H, 13C and 15N NMR assignments of an unusual Ca2+-binding protein from Entamoeba histolytica (EhCaBP6) in its apo form as a prelude to its structural and functional characterization.

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.

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)P

The highly conserved proteins of the 14-3-3 family are universal adaptors known to regulate an enormous range of cellular processes in eukaryotes. However, their biological functions remain largely uncharacterized in pathogenic protists comprising of several 14-3-3 protein isoforms. In this study, we report the role of 14-3-3 in coordinating cytoskeletal dynamics during phagocytosis in a professional phagocytic protist Entamoeba histolytica, the etiological agent of human amebiasis. There are three isoforms of 14-3-3 protein in amoeba and here we have investigated Eh14-3-3 Protein 3 (EhP3). Live and fixed cell imaging studies revealed the presence of this protein throughout the parasite phagocytosis process, with high rate of accumulation at the phagocytic cups and closed phagosomes. Conditional suppression of EhP3 expression caused significant defects in phagocytosis accompanied by extensive diminution of F-actin at the site of cup formation. Downregulated cells also exhibited defective recruitment of an F-actin stabilizing protein, EhCoactosin at the phagocytic cups. In addition, mass spectrometry based analysis further revealed a large group of EhP3-associated proteins, many of these proteins are known to regulate cytoskeletal architecture in E histolytica. The dynamics of these proteins may also be controlled by EhP3. Taken together, our findings strongly suggest that EhP3 is a novel and a key regulatory element of actin dynamics and phagocytosis in E. histolytica., Author summary Phagocytosis of host cells is central to pathogenesis of protist parasite Entamoeba histolytica, the etiological agent of human amebiasis. It is a complex and multistep process that requires dynamic remodelling of the actin cytoskeleton by a large number of scaffolding, signaling and actin-binding proteins (ABPs). Although several parasite ligands such as EhC2PK, EhCaBP1, EhCaBP3, EhAK1, Arp2/3 complex and EhCoactosin that participate in the phagocytic machinery have been identified, the mechanistic insights to their regulation process remain largely elusive. We have in this study identified and characterized the important role of scaffolding protein EhP3 in modulating cytoskeletal dynamics and regulating phagocytosis in E. histolytica. Expression knockdown, imaging and interaction studies suggest that EhP3 function as an adaptor molecule that controls the localization of an F-actin stabilizing protein EhCoactosin and thus the dynamics of F-actin rearrangement during phagocytosis. EhP3 also interact with other actin dynamics regulating proteins that may in coordination regulate cytoskeletal dynamics and thereby phagocytosis in Entamoeba.

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. Electronic supplementary material The online version of this article (10.1186/s12864-019-5570-z) contains supplementary material, which is available to authorized users.

The protist parasite Entamoeba histolytica causes amoebiasis, a major public health problem in developing countries. Only a small fraction of patients infected with the parasite display invasive disease involving colon or extra intestinal tissues such as liver. E. histolytica exists as two distinct forms, cysts, the infective form, and trophozoites, that are responsible for disease pathology. The latter multiply in the large intestine occasionally causing disease. The large intestine in humans is populated by a number of different bacterial communities and amoebic cells grow in their midst using some as food material. Several studies have shown relationship between bacteria and E. histolytica growth and virulence. However, an understanding of this relationship in human gut environment is not clear. We have investigated the possibility that there may be specific interaction of amoeba with different bacteria present in the gut environment by using a metagenomic pipe line. This was done by incubating bacteria isolated from human fecal material with E. histolytica and then identifying the bacterial population isolated from amoebic cells using a rRNA based metagenomic approach. Our results show that the parasite prefers a few bacterial species. One of these species is Lactobacillus ruminus which has never shown to be associated with E. histolytica.