The LINEs and SINEs of Entamoeba histolytica: comparative analysis and genomic distribution

A Small Genome amidst the Giants: Evidence of Genome Reduction in a Small Tubulinid Free-Living Amoeba

This study investigates the genomic characteristics of Echinamoeba silvestris, a small-sized amoeba within the Tubulinea clade of the Amoebozoa supergroup. Despite Tubulinea's significance in various fields, genomic data for this clade have been scarce. E. silvestris presents the smallest free-living amoeba genome within Tubulinea and Amoebozoa to date. Comparative analysis reveals intriguing parallels with parasitic lineages in terms of genome size and predicted gene numbers, emphasizing the need to understand the consequences of reduced genomes in free-living amoebae. Functional categorization of predicted genes in E. silvestris shows similar percentages of ortholog groups to other amoebae in various categories, but a distinctive feature is the extensive gene contraction in orphan (ORFan) genes and those involved in biological processes. Notably, among the few genes that underwent expansion, none are related to cellular components, suggesting adaptive processes that streamline biological processes and cellular components for efficiency and energy conservation. Additionally, our investigation into noncoding and repetitive elements sheds light on the evolution of genome size in amoebae, with E. silvestris distinguished by low percentage of repetitive elements. Furthermore, the analysis reveals that E. silvestris has the lowest mean number of introns per gene among the species studied, providing further support for its observed compact genome. Overall, this research underscores the diversity within Tubulinea, highlights knowledge gaps in Amoebozoa genomics, and positions E. silvestris as a valuable addition to genomic data sets, prompting further exploration of complexities in Amoebozoa diversity and genome evolution.

A Small Genome Amidst the Giants: Evidence of Genome Reduction in a Small Tubulinid Free-Living Amoeba

This study investigates the genomic characteristics of Echinamoeba silvestris , a small-sized amoeba within the Tubulinea clade of the Amoebozoa supergroup. Despite Tubulinea's significance in various fields, genomic data for this clade have been scarce . E. silvestris presents the smallest free-living amoeba genome within Tubulinea and Amoebozoa to date. Comparative analysis reveals intriguing parallels with parasitic lineages in terms of genome size and predicted gene numbers, emphasizing the need to understand the consequences of reduced genomes in free-living amoebae. Functional categorization of predicted genes in E. silvestris shows similar percentages of ortholog groups to other amoebae in various categories, but a distinctive feature is the extensive gene contraction in orphan (ORFan) genes and those involved in biological processes. Notably, among the few genes that underwent expansion, none are related to cellular components, suggesting adaptive processes that streamline biological processes and cellular components for efficiency and energy conservation. The investigation delves into genomic structural evidence, including gene content and repetitive elements, illuminating the distinctive genomic traits of E. silvestris and providing reinforcement for its compact genome size. Overall, this research underscores the diversity within Tubulinea, highlights knowledge gaps in Amoebozoa genomics, and positions E. silvestris as a valuable addition to genomic datasets, prompting further exploration of complexities in Amoebozoa diversity and genome evolution.

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.

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

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

Identification of oligo-adenylated small RNAs in the parasite Entamoeba and a potential role for small RNA control

Background

The RNA interference (RNAi) pathway is a gene regulation mechanism that utilizes small RNA (sRNA) and Argonaute (Ago) proteins to silence target genes. Our previous work identified a functional RNAi pathway in the protozoan parasite Entamoeba histolytica, including abundant 27 nt antisense sRNA populations which associate with EhAgo2–2 protein. However, there is lack of understanding about the sRNAs that are bound to two other EhAgos (EhAgo2–1 and 2–3), and the mechanism of sRNA regulation itself is unclear in this parasite. Therefore, identification of the entire pool of sRNA species and their sub-populations that associate with each individual EhAgo protein would be a major step forward.

Results

In the present study, we sequenced sRNA libraries from both total RNAs and EhAgo bound RNAs. We identified a new population of 31 nt sRNAs that results from the addition of a non-templated 3–4 adenosine nucleotides at the 3′-end of the 27 nt sRNAs, indicating a non-templated RNA-tailing event in the parasite. The relative abundance of these two sRNA populations is linked to the efficacy of gene silencing for the target gene when parasites are transfected with an RNAi-trigger construct, indicating that non-templated sRNA-tailing likely play a role in sRNA regulation in this parasite. We found that both sRNA populations (27 nt and 31 nt) are present in the related parasite Entamoeba invadens, and are unchanged during the development. In sequencing the sRNAs associating with the three EhAgo proteins, we observed that despite distinct cellular localization, all three EhAgo sRNA libraries contain 27 nt sRNAs with 5′-polyphosphate (5′-polyP) structure and share a largely overlapping sRNA repertoire. In addition, our data showed that a fraction of 31 nt sRNAs associate with EhAgo2–2 but not with its mutant protein (C-terminal deletion), nor other two EhAgos, indicating a specific EhAgo site may be required for sRNA modification process in the parasite.

Conclusion

We identified a new population of sRNA with non-templated oligo-adenylation modification, which is the first such observation amongst single celled protozoan parasites. Our sRNA sequencing libraries provide the first comprehensive sRNA dataset for all three Entamoeba Ago proteins, which can serve as a useful database for the amoeba community.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07275-6.

Whole genome sequencing of Entamoeba nuttalli reveals mammalian host-related molecular signatures and a novel octapeptide-repeat surface protein

The enteric protozoa Entamoeba histolytica is the causative agent of amebiasis, which is one of the most common parasitic diseases in developed and developing countries. Entamoeba nuttalli is the genetically closest species to E. histolytica in current phylogenetic analyses of Entamoeba species, and is prevalent in wild macaques. Therefore, E. nuttalli may be a key organism in which to investigate molecules required for infection of human or non-human primates. To explore the molecular signatures of host-parasite interactions, we conducted de novo assembly of the E. nuttalli genome, utilizing self-correction of PacBio long reads and polishing corrected reads using Illumina short reads, followed by comparative genomic analysis with two other mammalian and a reptilian Entamoeba species. The final draft assembly of E. nuttalli included 395 contigs with a total length of approximately 23 Mb, and 9,647 predicted genes, of which 6,940 were conserved with E. histolytica. In addition, we found an E. histolytica-specific repeat known as ERE2 in the E. nuttalli genome. GO-term enrichment analysis of mammalian host-related molecules indicated diversification of transmembrane proteins, including AIG1 family and BspA-like proteins that may be involved in the host-parasite interaction. Furthermore, we identified an E. nuttalli-specific protein that contained 42 repeats of an octapeptide ([G,E]KPTDTPS). This protein was shown to be localized on the cell surface using immunofluorescence. Since many repeat-containing proteins in parasites play important roles in interactions with host cells, this unique octapeptide repeat-containing protein may be involved in colonization of E. nuttalli in the intestine of macaques. Overall, our draft assembly provides a valuable resource for studying Entamoeba evolution and host-parasite selection.

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.

The trend in distribution of Q223R mutation of leptin receptor gene in amoebic liver abscess patients from North India: a prospective study

Host genetic susceptibility is an important risk factor in infectious diseases. We explored the distribution of Q223R mutation in leptin receptor gene of amoebic liver abscess (ALA) patients of North India. A total of 55 ALA samples along with 102 controls were subjected to PCR-RFLP analysis. The frequency of allele “G” (coding for arginine) was in general high in Indian population irrespective of the disease. Our results of Fisher exact test shows that heterozygous mutant (QQ versus QR, P = 0.049) and homozygous mutant (QQ versus RR, P = 0.004) were significantly associated with amoebic liver abscess when compared with homozygous wild (QQ).

Stress granule formation in Entamoeba histolytica: cross-talk between EhMLBP, EhRLE3 reverse transcriptase and polyubiquitinated proteins

The Entamoeba histolytica-methylated LINE-binding protein (EhMLBP) binds to methylated repetitive DNA and is a positive regulator of a reverse transcriptase of a long interspersed nucleotide element (LINE). This protein protects trophozoites against heat shock by reducing protein aggregation. The presence of EhMLBP and polyubiquitinated proteins in heat shock-induced protein aggregates raised the question whether these proteins interact. This assumption was confirmed by co-immunoprecipitation experiments: ubiquitinated proteins were detected in the perinuclear region of non-stressed E. histolytica trophozoites, whereas ubiquitinated proteins were detected in the perinuclear region and colocalized with EhMLBP in cytoplasmic granules in heat-shocked trophozoites. We also observed that overexpression of the reverse transcriptase of EhRLE3 induced the upregulation of EhMLBP expression and the formation of these EhMLBP-containing granules. Since (i) these EhMLBP-containing granules in the cytoplasm of heat-shocked E. histolytica trophozoites also contain polyubiquitinated proteins and poly(A)(+) mRNA and (ii) their formation is promoted by sodium arsenate, puromycin, and pateamine A and is inhibited by cycloheximide, we propose that these cytoplasmic EhMLBP-containing granules are stress granules. Our data also suggest that the formation of these granules is dependent upon EhMLBP and LINE.

Homologous Recombination Occurs in Entamoeba and Is Enhanced during Growth Stress and Stage Conversion

Homologous recombination (HR) has not been demonstrated in the parasitic protists Entamoeba histolytica or Entamoeba invadens, as no convenient method is available to measure it. However, HR must exist to ensure genome integrity, and possible genetic exchange, especially during stage conversion from trophozoite to cyst. Here we show the up regulation of mitotic and meiotic HR genes in Entamoeba during serum starvation, and encystation. To directly demonstrate HR we use a simple PCR-based method involving inverted repeats, which gives a reliable read out, as the recombination junctions can be determined by sequencing the amplicons. Using this read out, we demonstrate enhanced HR under growth stress in E. histolytica, and during encystation in E. invadens. We also demonstrate recombination between chromosomal inverted repeats. This is the first experimental demonstration of HR in Entamoeba and will help future investigations into this process, and to explore the possibility of meiosis in Entamoeba.

Genomic distribution of SINEs in Entamoeba histolytica strains: implication for genotyping

BACKGROUND

The major clinical manifestations of Entamoeba histolytica infection include amebic colitis and liver abscess. However the majority of infections remain asymptomatic. Earlier reports have shown that some E. histolytica isolates are more virulent than others, suggesting that virulence may be linked to genotype. Here we have looked at the genomic distribution of the retrotransposable short interspersed nuclear elements EhSINE1 and EhSINE2. Due to their mobile nature, some EhSINE copies may occupy different genomic locations among isolates of E. histolytica possibly affecting adjacent gene expression; this variability in location can be exploited to differentiate strains.

RESULTS

We have looked for EhSINE1- and EhSINE2-occupied loci in the genome sequence of Entamoeba histolytica HM-1:IMSS and searched for homologous loci in other strains to determine the insertion status of these elements. A total of 393 EhSINE1 and 119 EhSINE2 loci were analyzed in the available sequenced strains (Rahman, DS4-868, HM1:CA, KU48, KU50, KU27 and MS96-3382. Seventeen loci (13 EhSINE1 and 4 EhSINE2) were identified where a EhSINE1/EhSINE2 sequence was missing from the corresponding locus of other strains. Most of these loci were unoccupied in more than one strain. Some of the loci were analyzed experimentally for SINE occupancy using DNA from strain Rahman. These data helped to correctly assemble the nucleotide sequence at three loci in Rahman. SINE occupancy was also checked at these three loci in 7 other axenically cultivated E. histolytica strains and 16 clinical isolates. Each locus gave a single, specific amplicon with the primer sets used, making this a suitable method for strain typing. Based on presence/absence of SINE and amplification with locus-specific primers, the 23 strains could be divided into eleven genotypes. The results obtained by our method correlated with the data from other typing methods. We also report a bioinformatic analysis of EhSINE2 copies.

CONCLUSIONS

Our results reveal several loci with extensive polymorphism of SINE occupancy among different strains of E. histolytica and prove the principle that the genomic distribution of SINEs is a valid method for typing of E. histolytica strains.

Small RNA pyrosequencing in the protozoan parasite Entamoeba histolytica reveals strain-specific small RNAs that target virulence genes

BACKGROUND

Small RNA mediated gene silencing is a well-conserved regulatory pathway. In the parasite Entamoeba histolytica an endogenous RNAi pathway exists, however, the depth and diversity of the small RNA population remains unknown.

RESULTS

To characterize the small RNA population that associates with E. histolytica Argonaute-2 (EhAGO2-2), we immunoprecipitated small RNAs that associate with it and performed one full pyrosequencing run. Data analysis revealed new features of the 27nt small RNAs including the 5'-G predominance, distinct small RNA distribution patterns on protein coding genes, small RNAs mapping to both introns and exon-exon junctions, and small RNA targeted genes that are clustered particularly in sections of genome duplication. Characterization of genomic loci to which both sense and antisense small RNAs mapped showed that both sets of small RNAs have 5'-polyphosphate termini; strand-specific RT-PCR detected transcripts in both directions at these loci suggesting that both transcripts may serve as template for small RNA generation. In order to determine whether small RNA abundance patterns account for strain-specific gene expression profiles of E. histolytica virulent and non-virulent strains, we sequenced small RNAs from a non-virulent strain and found that small RNAs mapped to genes in a manner consistent with their regulation of strain-specific virulence genes.

CONCLUSIONS

We provided a full spectrum analysis for E. histolytica AGO2-2 associated 27nt small RNAs. Additionally, comparative analysis of small RNA populations from virulent and non-virulent amebic strains indicates that small RNA populations may regulate virulence genes.

Recombinant SINEs are formed at high frequency during induced retrotransposition in vivo

Non-long terminal repeat Retrotransposons are referred to as long interspersed nuclear elements (LINEs) and their non-autonomous partners are short interspersed nuclear elements (SINEs). It is believed that an active SINE copy, upon retrotransposition, generates near identical copies of itself, which subsequently accumulate mutations resulting in sequence polymorphism. Here we show that when a retrotransposition-competent cell line of the parasitic protist Entamoeba histolytica, transfected with a marked SINE copy, is induced to retrotranspose, >20% of the newly retrotransposed copies are neither identical to the marked SINE nor to the mobilized resident SINEs. Rather they are recombinants of resident SINEs and the marked SINE. They are a consequence of retrotransposition and not DNA recombination, as they are absent in cells not expressing the retrotransposition functions. This high-frequency recombination provides a new explanation for the existence of mosaic SINEs, which may impact on genetic analysis of SINE lineages, and measurement of phylogenetic distances.

Distribution of MGEs and their insertion sites in the Macaca mulatta genome

Mobile genetic elements (MGEs) are fragments of DNA that can move around within the genome through retrotransposition. These are responsible for various important events such as gene inactivation, transduction, regulation of gene expression and genome expansion. The present work involves the identification and study of the distribution of Alu and L1 retrotransposons in the genome of Macaca mulatta, an extensively used organism in biomedical studies. We also make comparisons with MGE distributions in other primate genomes and study the physicochemical properties of the local DNA structure around the transposon insertion site using ELAN. The present work also includes computational testing of the pre-insertion loci in order to detect unique features based on DNA structure, thermodynamic considerations and protein interaction measures. Although there is significant sequence divergence between the elements of M. mulatta and H. sapiens, their genome wide distribution is very similar; comparing the distribution of L1's in all available X chromosome sequences suggests a common mechanism behind the spread of MGE's in primate genomes.

Differential distribution of a SINE element in the Entamoeba histolytica and Entamoeba dispar genomes: role of the LINE-encoded endonuclease

BACKGROUND

Entamoeba histolytica and Entamoeba dispar are closely related protistan parasites but while E. histolytica can be invasive, E. dispar is completely non pathogenic. Transposable elements constitute a significant portion of the genome in these species; there being three families of LINEs and SINEs. These elements can profoundly influence the expression of neighboring genes. Thus their genomic location can have important phenotypic consequences. A genome-wide comparison of the location of these elements in the E. histolytica and E. dispar genomes has not been carried out. It is also not known whether the retrotransposition machinery works similarly in both species. The present study was undertaken to address these issues.

RESULTS

Here we extracted all genomic occurrences of full-length copies of EhSINE1 in the E. histolytica genome and matched them with the homologous regions in E. dispar, and vice versa, wherever it was possible to establish synteny. We found that only about 20% of syntenic sites were occupied by SINE1 in both species. We checked whether the different genomic location in the two species was due to differences in the activity of the LINE-encoded endonuclease which is required for nicking the target site. We found that the endonucleases of both species were essentially very similar, both in their kinetic properties and in their substrate sequence specificity. Hence the differential distribution of SINEs in these species is not likely to be influenced by the endonuclease. Further we found that the physical properties of the DNA sequences adjoining the insertion sites were similar in both species.

CONCLUSIONS

Our data shows that the basic retrotransposition machinery is conserved in these sibling species. SINEs may indeed have occupied all of the insertion sites in the genome of the common ancestor of E. histolytica and E. dispar but these may have been subsequently lost from some locations. Alternatively, SINE expansion took place after the divergence of the two species. The absence of SINE1 in 80% of syntenic loci could affect the phenotype of the two species, including their pathogenic properties, which needs to be explored.

Bioinformatic analysis of Entamoeba histolytica SINE1 elements

BACKGROUND

Invasive amoebiasis, caused by infection with the human parasite Entamoeba histolytica remains a major cause of morbidity and mortality in some less-developed countries. Genetically E. histolytica exhibits a number of unusual features including having approximately 20% of its genome comprised of repetitive elements. These include a number of families of SINEs - non-autonomous elements which can, however, move with the help of partner LINEs. In many eukaryotes SINE mobility has had a profound effect on gene expression; in this study we concentrated on one such element - EhSINE1, looking in particular for evidence of recent transposition.

RESULTS

EhSINE1s were detected in the newly reassembled E. histolytica genome by searching with a Hidden Markov Model developed to encapsulate the key features of this element; 393 were detected. Examination of their sequences revealed that some had an internal structure showing one to four 26-27 nt repeats. Members of the different classes differ in a number of ways and in particular those with two internal repeats show the properties expected of fairly recently transposed SINEs - they are the most homogeneous in length and sequence, they have the longest (i.e. the least decayed) target site duplications and are the most likely to show evidence (in a cDNA library) of active transcription. Furthermore we were able to identify 15 EhSINE1s (6 pairs and one triplet) which appeared to be identical or very nearly so but inserted into different sites in the genome; these provide good evidence that if mobility has now ceased it has only done so very recently.

CONCLUSIONS

Of the many families of repetitive elements present in the genome of E. histolytica we have examined in detail just one - EhSINE1. We have shown that there is evidence for waves of transposition at different points in the past and no evidence that mobility has entirely ceased. There are many aspects of the biology of this parasite which are not understood, in particular why it is pathogenic while the closely related species E. dispar is not, the great genetic diversity found amongst patient isolates and the fact, which may be related, that only a small proportion of those infected develop clinical invasive amoebiasis. Mobile genetic elements, with their ability to alter gene expression may well be important in unravelling these puzzles.

The biology and evolution of transposable elements in parasites

Transposable elements (TEs) are dynamic elements that can reshape host genomes by generating rearrangements with the potential to create or disrupt genes, to shuffle existing genes, and to modulate their patterns of expression. In the genomes of parasites that infect mammals several TEs have been identified that probably have been maintained throughout evolution due to their contribution to gene function and regulation of gene expression. This review addresses how TEs are organized, how they colonize the genomes of mammalian parasites, the functional role these elements play in parasite biology, and the interactions between these elements and the parasite genome.

Identification of four Entamoeba histolytica organellar DNA polymerases of the family B and cellular localization of the Ehodp1 gene and EhODP1 protein

We report the identification of a family of four active genes (Ehodp1, Ehodp2, Ehodp3, and Ehodp4) encoding putative DNA polymerases in Entamoeba histolytica, the protozoan parasite responsible of human amoebiasis. The four Ehodp genes show similarity to DNA polymerases encoded in fungi and plant mitochondrial plasmids. EhODP polypeptides conserve the 3′-5′ exonuclease II and 5′-3′ polymerization domains, and they have the I, II, and III conserved boxes that characterize them as DNA polymerases of family B. Furthermore, we found in EhODP polymerases two novel A and B boxes, present also in DNA polymerases encoded in fungi mitochondrial plasmids. By in situ PCR, Ehodp1 gene was located in nuclei and in DNA-containing cytoplasmic structures. Additionally, using polyclonal antibodies against a recombinant rEhODP1-168 polypeptide, and confocal microscopy, EhODP1 was located in cytoplasmic DNA-containing structures.

Genome wide survey, discovery and evolution of repetitive elements in three Entamoeba species

Background

Identification and mapping of repetitive elements is a key step for accurate gene prediction and overall structural annotation of genomes. During the assembly and annotation of three highly repetitive amoeba genomes, Entamoeba histolytica, Entamoeba dispar, and Entamoeba invadens, we performed comparative sequence analysis to identify and map all class I and class II transposable elements in their sequences.

Results

Here, we report the identification of two novel Entamoeba-specific repeats: ERE1 and ERE2; ERE1 is spread across the three genomes and associated with different repeats in a species-specific manner, while ERE2 is unique to E. histolytica. We also report the identification of two novel subfamilies of LINE and SINE retrotransposons in E. dispar and provide evidence for how the different LINE and SINE subfamilies evolved in these species. Additionally, we found a putative transposase-coding gene in E. histolytica and E. dispar related to the mariner transposon Hydargos from E. invadens. The distribution of transposable elements in these genomes is markedly skewed with a tendency of forming clusters. More than 70% of the three genomes have a repeat density below their corresponding average value indicating that transposable elements are not evenly distributed. We show that repeats and repeat-clusters are found at syntenic break points between E. histolytica and E. dispar and hence, could work as recombination hot spots promoting genome rearrangements.

Conclusion

The mapping of all transposable elements found in these parasites shows that repeat coverage is up to three times higher than previously reported. LINE, ERE1 and mariner elements were present in the common ancestor to the three Entamoeba species while ERE2 was likely acquired by E. histolytica after its separation from E. dispar. We demonstrate that E. histolytica and E. dispar share their entire repertoire of LINE and SINE retrotransposons and that Eh_SINE3/Ed_SINE1 originated as a chimeric SINE from Eh/Ed_SINE2 and Eh_SINE1/Ed_SINE3. Our work shows that transposable elements are organized in clusters, frequently found at syntenic break points providing insights into their contribution to chromosome instability and therefore, to genomic variation and speciation in these parasites.

DNA transposons and the evolution of eukaryotic genomes

Transposable elements are mobile genetic units that exhibit broad diversity in their structure and transposition mechanisms. Transposable elements occupy a large fraction of many eukaryotic genomes and their movement and accumulation represent a major force shaping the genes and genomes of almost all organisms. This review focuses on DNA-mediated or class 2 transposons and emphasizes how this class of elements is distinguished from other types of mobile elements in terms of their structure, amplification dynamics, and genomic effect. We provide an up-to-date outlook on the diversity and taxonomic distribution of all major types of DNA transposons in eukaryotes, including Helitrons and Mavericks. We discuss some of the evolutionary forces that influence their maintenance and diversification in various genomic environments. Finally, we highlight how the distinctive biological features of DNA transposons have contributed to shape genome architecture and led to the emergence of genetic innovations in different eukaryotic lineages.

Structure and content of the Entamoeba histolytica genome

The intestinal parasite Entamoeba histolytica is one of the first protists for which a draft genome sequence has been published. Although the genome is still incomplete, it is unlikely that many genes are missing from the list of those already identified. In this chapter we summarise the features of the genome as they are currently understood and provide previously unpublished analyses of many of the genes.

Antisense RNA and RNAi in protozoan parasites: working hard or hardly working?

The complex life cycles of many protozoan parasites require the ability to respond to environmental and developmental cues through regulated gene expression. Traditionally, parasitologists have investigated these mechanisms by identifying and characterizing proteins that are necessary for the regulated expression of the genetic material. Although often successful, it is clear that protein-mediated gene regulation is only part of a complex story in which RNA itself is endowed with regulatory functions. Herein, we review both the known and potential regulatory roles of two types of RNA pathways within protozoan parasites: the RNA interference pathway and natural antisense transcripts. A better understanding of the native role of these pathways will not only enhance our understanding of the biology of these organisms but also aid in the development of more robust tools for reverse genetic analysis in this post-genomic era.

Growth of the protozoan parasite Entamoeba histolytica in 5-azacytidine has limited effects on parasite gene expression

Background

In higher eukaryotes DNA methylation regulates important biological functions including silencing of gene expression and protection from adverse effects of retrotransposons. In the protozoan parasite Entamoeba histolytica, a DNA methyltransferase has been identified and treatment with 5-azacytidine (5-AzaC), a potent inhibitor of DNA methyltransferase, has been reported to attenuate parasite virulence. However, the overall extent of DNA methylation and its subsequent effects on global gene expression in this parasite are currently unknown.

Results

In order to identify the genome-wide effects of DNA methylation in E. histolytica, we used a short oligonucleotide microarray representing 9,435 genes (~95% of all annotated amebic genes) and compared the expression profile of E. histolytica HM-1:IMSS parasites with those treated with 23 μM 5-AzaC for up to one week. Overall, 2.1% of genes tested were transcriptionally modulated under these conditions. 68 genes were upregulated and 131 genes down regulated (2-fold change; p-value < 0.05). Sodium-bisulfite treatment and sequencing of genes indicated that there were at least two subsets of genes with genomic DNA methylation in E. histolytica: (i) genes that were endogenously silenced by genomic DNA methylation and for which 5-AzaC treatment induced transcriptional de-repression, and (ii) genes that have genomic DNA methylation, but which were not endogenously silenced by the methylation. We identified among the genes down regulated by 5-AzaC treatment a cysteine proteinase (2.m00545) and lysozyme (52.m00148) both of which have known roles in amebic pathogenesis. Decreased expression of these genes in the 5-AzaC treated E. histolytica may account in part for the parasites reduced cytolytic abilities.

Conclusion

This work represents the first genome-wide analysis of DNA-methylation in Entamoeba histolytica and indicates that DNA methylation has relatively limited effects on gene expression in this parasite.

Molecular methods for diagnosis of Entamoeba histolytica in a clinical setting: an overview

The range of clinical outcomes following Entamoeba histolytica infection is likely to be influenced by the different strains of the parasite already existing in our population. There is a need for developing faster, reliable and reproducible methods for identifying the different strains of E. histolytica. This would have a major impact on the subsequent course of treatment given to patients. In the post-genomic era, different loci of the Entamoeba genome have been targeted for developing suitable probes and genetic markers. This review highlights the development made in this direction and the possibility of using these methods for routine testing of this parasite in clinical samples.

SINE elements of Entamoeba dispar

Entamoeba histolytica and E. dispar are closely related protozoan parasites; the former causes clinical amoebiasis in humans while the latter appears to be non-pathogenic. The molecular biology of E. histolytica shows a number of unusual features, one of which is the abundance of polyadenylated but apparently untranslatable mRNAs produced; many of these are the product of at least three families of SINEs (EhSINE1-3). Here we show that the genome of E. dispar contains numerous copies of a SINE element (EdSINE1) whose 5'- and 3'-ends are very similar to those of EhSINE1 but with a much less similar middle portion. Twelve out of 18 copies examined had target site duplications. In none out of six cases examined was there a SINE element in the homologous region of the E. histolytica genome but a single copy of EdSINE1 is present in E. histolytica where it is identified as EhSINE3.

Identification of insertion hot spots for non-LTR retrotransposons: computational and biochemical application to Entamoeba histolytica

The genome of the human pathogen Entamoeba histolytica contains non-long terminal repeat (LTR) retrotransposons, the EhLINEs and EhSINEs, which lack targeted insertion. We investigated the importance of local DNA structure, and sequence preference of the element-encoded endonuclease (EN) in selecting target sites for retrotransposon insertion. Pre-insertion loci were tested computationally to detect unique features based on DNA structure, thermodynamic considerations and protein interaction measures. Target sites could readily be distinguished from other genomic sites based on these criteria. The contribution of the EhLINE1-encoded EN in target site selection was investigated biochemically. The sequence-specificity of the EN was tested in vitro with a variety of mutated substrates. It was possible to assign a consensus sequence, 5′-GCATT-3′, which was efficiently nicked between A-T and T-T. The upstream G residue enhanced EN activity, possibly serving to limit retrotransposition in the A+T-rich E.histolytica genome. Mutated substrates with poor EN activity showed structural differences compared with normal substrates. Analysis of retrotransposon insertion sites from a variety of organisms showed that, in general, regions of favorable DNA structure were recognized for retrotransposition. A combination of favorable DNA structure and preferred EN nicking sequence in the vicinity of this structure may determine the genomic hotspots for retrotransposition.

Stress by heat shock induces massive down regulation of genes and allows differential allelic expression of the Gal/GalNAc lectin in Entamoeba histolytica

Gene expression analysis by microarray assay revealed that when exposed to stress, Entamoeba histolytica exhibits a specific heat shock response, together with a dramatic overall reduction in gene transcription as well as differential allelic expression of key genes participating in virulence, such as the galactose/N-acetylgalactosamine (Gal/GalNAc) lectin.

Epigenetic silencing of gene expression in Entamoeba histolytica

Transcriptional silencing of an amebapore (ap-a) gene occurred in Entamoeba histolytica following the transfection of plasmids containing a DNA segment (473 bp) homologous to the 5' upstream region of the gene. This segment contains the promoter region of the ap-a gene, a T-rich stretch, followed by a truncated SINE1 (short interspersed element) that is transcribed from the opposite strand. The downstream silencing of the ap-a gene did not occur with plasmids containing the entire SINE1 sequence or lacking the entire SINE1 sequences including the T-rich stretch. Such plasmids promoted the overexpression of the ap-a gene. The transcription of the SINE element required both the T-rich stretch as well as sequences from the 5' end of SINE. RNA extracts from gene-silenced cultures showed small amounts of short (approximately 140 nt), single-stranded molecules with homology to SINE1 transcripts but no siRNA. Chromatin immunoprecipitation (ChIP) analysis of silenced G3 trophozoites with an antibody against methylated K4 of histone H3 revealed a demethylation of K4 at the domain of the ap-a gene indicating transcriptional inactivation. These results suggest the involvement of the SINE1 element in triggering the gene silencing and the role of histone modification in its epigenetic maintenance. The avirulent phenotype of the silenced trophozoites was demonstrated in various assays and the results suggest they may have a potential use for vaccination.

Involvement of a short interspersed element in epigenetic transcriptional silencing of the amoebapore gene in Entamoeba histolytica

Transcriptional silencing of an amoebapore (ap-a) gene occurred in Entamoeba histolytica following the transfection of plasmids containing a DNA segment (473 bp) homologous to the 5' upstream region of the gene (R. Bracha, Y. Nuchamowitz, and D. Mirelman, Eukaryot. Cell 2:295-305, 2003). This segment contains the promoter region of the ap-a gene, a T-rich stretch, followed by a truncated SINE1 (short interspersed element 1) that is transcribed from the antisense strand. Transfection of plasmids containing truncated SINE1 sequences which lack their 3' regulatory elements upstream of the ap-a gene was essential for the downstream silencing of the ap-a gene while transfection with plasmids containing the entire SINE1 sequence or without the T-rich stretch promoted the overexpression of the ap-a gene. Both the T-rich stretch and sequences of the 5' SINE1 were essential for the transcription of SINE1. RNA extracts from gene-silenced cultures showed small amounts of short (approximately 140-nucleotide), single-stranded molecules with homology to SINE1 but no short interfering RNA. Chromatin immunoprecipitation analysis with an antibody against methylated K4 of histone H3 showed a demethylation of K4 at the domain of the ap-a gene, indicating transcriptional inactivation. These results suggest the involvement of SINE1 in triggering the gene silencing and the role of histone modification in its epigenetic maintenance.

Identification of differentially expressed genes in virulent and nonvirulent Entamoeba species: potential implications for amebic pathogenesis

Entamoeba histolytica is a protozoan parasite that causes colitis and liver abscesses. Several Entamoeba species and strains with differing levels of virulence have been identified. E. histolytica HM-1:IMSS is a virulent strain, E. histolytica Rahman is a nonvirulent strain, and Entamoeba dispar is a nonvirulent species. We used an E. histolytica DNA microarray consisting of 2,110 genes to assess the transcriptional differences between these species/strains with the goal of identifying genes whose expression correlated with a virulence phenotype. We found 415 genes expressed at lower levels in E. dispar and 32 genes with lower expression in E. histolytica Rahman than in E. histolytica HM-1:IMSS. Overall, 29 genes had decreased expression in both the nonvirulent species/strains than the virulent E. histolytica HM-1:IMSS. Interestingly, a number of genes with potential roles in stress response and virulence had decreased expression in either one or both nonvirulent Entamoeba species/strains. These included genes encoding Fe hydrogenase (9.m00419), peroxiredoxin (176.m00112), type A flavoprotein (6.m00467), lysozyme (6.m00454), sphingomyelinase C (29.m00231), and a hypothetical protein with homology to both a Plasmodium sporozoite threonine-asparagine-rich protein (STARP) and a streptococcal hemagglutinin (238.m00054). The function of these genes in Entamoeba and their specific roles in parasite virulence need to be determined. We also found that a number of the non-long-terminal-repeat retrotransposons (EhLINEs and EhSINEs), which have been shown to modulate gene expression and genomic evolution, had lower expression in the nonvirulent species/strains than in E. histolytica HM-1:IMSS. Our results, identifying expression profiles and patterns indicative of a virulence phenotype, may be useful in characterizing the transcriptional framework of virulence.

Comparative genomics of Dictyostelium discoideum and Entamoeba histolytica

Amoebozoa represent one of the earliest branches from the last common ancestor of all eukaryotes and contain some of the most dangerous human pathogens. Two amoebozoan genomes -- from the model organism Dictyostelium discoideum and the human pathogen Entamoeba histolytica -- have been published this year. Owing to their high A+T content, both genomes were difficult to sequence. In addition to nine amoebozoan expressed sequence tag projects, efforts are underway for comparative sequencing of four additional Entamoeba species. The completed genome sequences of D. discoideum and E. histolytica revealed unusual telomere structures, a high percentage of repetitive elements and a remarkably high gene content that is close to the one of Drosophila melanogaster. Finally, both organisms are brilliant examples of the influence of the lifestyle of an organism on its genome.