An Entamoeba histolytica LINE/SINE pair inserts at common target sites cleaved by the restriction enzyme-like LINE-encoded endonuclease

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.

Distinct Retrotransposon Evolution Profile in the Genome of Rabbit (Oryctolagus cuniculus)

Although the rabbit genome has already been annotated, it is mobilome remains largely unknown. Here, multiple pipelines were used to de novo mine and annotate the mobilome in rabbit. Four families and 19 subfamilies of LINE1s, two families and nine subfamilies of SINEs, and 12 ERV families were defined in rabbit based on sequence identity, structural organization, and phylogenetic tree. The analysis of insertion age and polymerase chain reaction suggests that a number of families are very young and may remain active, such as L1B, L1D, OcuSINEA, and OcuERV1. RepeatMasker annotation revealed a distinct transposable element landscape within the genome, with approximately two million copies of SINEs, representing the greatest proportion of the genome (19.61%), followed by LINEs (15.44%), and LTRs (4.11%), respectively, considerably different from most other mammal mobilomes except hedgehog and tree shrew, in which LINEs have the highest proportion. Furthermore, a very high rate of insertion polymorphisms (>85%) for the youngest subfamily (OcuSINEA1) was identified by polymerase chain reaction. The majority of retrotransposon insertions overlapped with protein-coding regions (>80%) and lncRNA (90%) genes. Genomic distribution bias was observed for retrotransposons, with those immediately upstream (-1 kb) and downstream (1 kb) of genes significantly depleted. Local GC content in 50-kb widows had significantly negative correlations with LINE (rs=-0.996) and LTR (rs=-0.829) insertions. The current study revealed a distinct mobilome landscape in rabbit, which will assist in the elucidation of the evolution of the genome of lagomorphs, and even other mammals.

Transcriptomic analysis of Entamoeba histolytica reveals domain-specific sense strand expression of LINE-encoded ORFs with massive antisense expression of RT domain

LINEs are retrotransposable elements found in diverse organisms. Their activity is kept in check by several mechanisms, including transcriptional silencing. Here we have analyzed the transcription status of LINE1 copies in the early-branching parasitic protist Entamoeba histolytica. Full-length EhLINE1 encodes ORF1, and ORF2 with reverse transcriptase (RT) and endonuclease (EN) domains. RNA-Seq analysis of EhLINE1 copies (both truncated and full-length) showed unique features. Firstly, although 20/41 transcribed copies were full-length, we failed to detect any full-length transcripts. Rather, sense-strand transcripts mapped to the functional domains- ORF1, RT and EN. Secondly, there was strong antisense transcription specifically from RT domain. No antisense transcripts were seen from ORF1. Antisense RT transcripts did not encode known functional peptides. They could possibly be involved in attenuating translation of RT domain, as we failed to detect ORF2p, whereas ORF1p was detectable. Lack of full-length transcripts and strong antisense RT expression may serve to limit EhLINE1 retrotransposition.

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

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

Site-specific non-LTR retrotransposons

Although most of non-long terminal repeat (non-LTR) retrotransposons are incorporated in the host genome almost randomly, some non-LTR retrotransposons are incorporated into specific sequences within a target site. On the basis of structural and phylogenetic features, non-LTR retrotransposons are classified into two large groups, restriction enzyme-like endonuclease (RLE)-encoding elements and apurinic/apyrimidinic endonuclease (APE)-encoding elements. All clades of RLE-encoding non-LTR retrotransposons include site-specific elements. However, only two of more than 20 APE-encoding clades, Tx1 and R1, contain site-specific non-LTR elements. Site-specific non-LTR retrotransposons usually target within multi-copy RNA genes, such as rRNA gene (rDNA) clusters, or repetitive genomic sequences, such as telomeric repeats; this behavior may be a symbiotic strategy to reduce the damage to the host genome. Site- and sequence-specificity are variable even among closely related non-LTR elements and appeared to have changed during evolution. In the APE-encoding elements, the primary determinant of the sequence- specific integration is APE itself, which nicks one strand of the target DNA during the initiation of target primed reverse transcription (TPRT). However, other factors, such as interaction between mRNA and the target DNA, and access to the target region in the nuclei also affect the sequence-specificity. In contrast, in the RLE-encoding elements, DNA-binding motifs appear to affect their sequence-specificity, rather than the RLE domain itself. Highly specific integration properties of these site-specific non-LTR elements make them ideal alternative tools for sequence-specific gene delivery, particularly for therapeutic purposes in human diseases.

Origins and evolution of viruses of eukaryotes: The ultimate modularity

Viruses and other selfish genetic elements are dominant entities in the biosphere, with respect to both physical abundance and genetic diversity. Various selfish elements parasitize on all cellular life forms. The relative abundances of different classes of viruses are dramatically different between prokaryotes and eukaryotes. In prokaryotes, the great majority of viruses possess double-stranded (ds) DNA genomes, with a substantial minority of single-stranded (ss) DNA viruses and only limited presence of RNA viruses. In contrast, in eukaryotes, RNA viruses account for the majority of the virome diversity although ssDNA and dsDNA viruses are common as well. Phylogenomic analysis yields tangible clues for the origins of major classes of eukaryotic viruses and in particular their likely roots in prokaryotes. Specifically, the ancestral genome of positive-strand RNA viruses of eukaryotes might have been assembled de novo from genes derived from prokaryotic retroelements and bacteria although a primordial origin of this class of viruses cannot be ruled out. Different groups of double-stranded RNA viruses derive either from dsRNA bacteriophages or from positive-strand RNA viruses. The eukaryotic ssDNA viruses apparently evolved via a fusion of genes from prokaryotic rolling circle-replicating plasmids and positive-strand RNA viruses. Different families of eukaryotic dsDNA viruses appear to have originated from specific groups of bacteriophages on at least two independent occasions. Polintons, the largest known eukaryotic transposons, predicted to also form virus particles, most likely, were the evolutionary intermediates between bacterial tectiviruses and several groups of eukaryotic dsDNA viruses including the proposed order "Megavirales" that unites diverse families of large and giant viruses. Strikingly, evolution of all classes of eukaryotic viruses appears to have involved fusion between structural and replicative gene modules derived from different sources along with additional acquisitions of diverse genes.

A genomewide overexpression screen identifies genes involved in the phosphatidylinositol 3-kinase pathway in the human protozoan parasite Entamoeba histolytica

Entamoeba histolytica is a protozoan parasite that causes amoebic dysentery and liver abscess. E. histolytica relies on motility, phagocytosis, host cell adhesion, and proteolysis of extracellular matrix for virulence. In eukaryotic cells, these processes are mediated in part by phosphatidylinositol 3-kinase (PI3K) signaling. Thus, PI3K may be critical for virulence. We utilized a functional genomics approach to identify genes whose products may operate in the PI3K pathway in E. histolytica. We treated a population of trophozoites that were overexpressing genes from a cDNA library with a near-lethal dose of the PI3K inhibitor wortmannin. This screen was based on the rationale that survivors would be overexpressing gene products that directly or indirectly function in the PI3K pathway. We sequenced the overexpressed genes in survivors and identified a cDNA encoding a Rap GTPase, a protein previously shown to participate in the PI3K pathway. This supports the validity of our approach. Genes encoding a coactosin-like protein, EhCoactosin, and a serine-rich E. histolytica protein (SREHP) were also identified. Cells overexpressing EhCoactosin or SREHP were also less sensitive to a second PI3K inhibitor, LY294002. This corroborates the link between these proteins and PI3K. Finally, a mutant cell line with an increased level of phosphatidylinositol (3,4,5)-triphosphate, the product of PI3K activity, exhibited increased expression of SREHP and EhCoactosin. This further supports the functional connection between these proteins and PI3K in E. histolytica. To our knowledge, this is the first forward-genetics screen adapted to reveal genes participating in a signal transduction pathway in this pathogen.

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.

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.

Targeting novel sites: The N-terminal DNA binding domain of non-LTR retrotransposons is an adaptable module that is implicated in changing site specificities

Restriction-like endonuclease (RLE) bearing non-LTR retrotransposons are site-specific elements that integrate into the genome through target primed reverse transcription (TPRT). RLE-bearing elements have been used as a model system for investigating non-LTR retrotransposon integration. R2 elements target a specific site in the 28S rDNA gene. We previously demonstrated that the two major sub-classes of R2 (R2-A and R2-D) target the R2 insertion site in an opposing manner with regard to the pairing of known DNA binding domains and bound sequences-indicating that the A- and D-clades represent independently derived modes of targeting that site. Elements have been discovered that group phylogenetically with R2 but do not target the canonical R2 site. Here we extend our earlier studies to show that a separate R2-A clade element, which targets a site other than the canonical R2 site, does so by using the N-terminal zinc fingers and Myb motifs. We further extend our targeting studies beyond R2 clade elements by investigating the ability of the N-terminal zinc fingers from the nematode NeSL-1 element to target its integration site. Our data are consistent with the use of an N-terminal DNA binding domain as one of the major targeting determinants used by RLE-bearing non-LTR retrotransposons to secure a protein subunit near the insertion site. This N-terminal DNA binding domain can undergo modifications, allowing the element to target novel sites. The binding orientation of the N-terminal domain relative to the insertion site is quite variable.

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.

Independently derived targeting of 28S rDNA by A- and D-clade R2 retrotransposons: Plasticity of integration mechanism

Restriction-like endonuclease (RLE) bearing non-LTR retrotransposons are site-specific elements that integrate into the genome through a target primed reverse transcription mechanism (TPRT). R2 elements have been used as a model system for investigating non-LTR retrotransposon integration. We previously demonstrated that R2 retrotransposons require two subunits of the element-encoded multifunctional protein to integrate-one subunit bound upstream of the insertion site and one bound downstream. R2 elements have been phylogenetically categorized into four clades: R2-A, B, C and D, that diverged from a common ancestor more than 850 million years ago. All R2 elements target the same sequence within 28S rDNA. The amino-terminal domain of R2Bm, an R2-D clade element, contains a single zinc finger and a Myb motif that are responsible for binding R2 protein downstream of the insertion site. Target site recognition is of interest as it is the first step in the integration reaction and may help elucidate evolutionary history and integration mechanism. The amino-terminal domain of R2-A clade members contains three zinc fingers and a Myb motif. We show here that R2Lp, an R2-A clade member, uses its amino-terminal DNA binding motifs to bind upstream of the insertion site. Because the R2-A and R2-D clade elements recognize 28S rDNA differently, we conclude the A- and D-clades represent independent targeting events to the 28S site. Our results also indicate a certain plasticity of insertional mechanics exists between the two clades.

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.

Target site analysis of RTE1_LA and its AfroSINE partner in the elephant genome

SINEs retrotranspose using their partner LINE's enzymatic machinery. It has recently been proposed that AfroSINEs ending with GGTTT 3' tandem repeats were mobilized by RTE elements ending with CAA 3' tandem repeats in the Afrotherian genome. Using sequences from the elephant genome, we show that AfroSINEs derive from RTE ending with GGTTT-like 3' tandem repeats, a subgroup of RTE1_LA that only reached low copy number, and confirm that they were most likely mobilized by RTE ending with CAA(n) tandem repeats (RTE1_LA-CAA(n)). This partnership is supported by sequence similarity between two regions of the elements, overlap in the timing of their activity, common features of their target site consensus that are not shared by other members of the RTE family, and their high copy number. Detailed analyses of pre-insertion loci reveal that like many other apurinic/apyrimidinic endonuclease encoding elements, RTE1_LA-CAA(n) shows loose target site specificity. In addition, the RTE1_LA-CAA(n) target site consensus shares several structural and primary sequence features with that of LINE1, suggesting that these two elements share close functional similarity in the target primed reverse transcription (TPRT) reaction. Interestingly, although globally similar, the target site consensus of AfroSINE(Anc) and RTE1_LA-CAA(n) differ in several aspects. These differences, not observed among all SINE/LINE pairs so far examined, are most likely due to the fact that AfroSINEs and RTE1_LA-CAA(n) are terminated by a different tandem repeat motif. We propose that these differences reflect constraints imposed by base pairing interactions between the mRNA 3' terminal tandem repeats and the target DNA at the onset of TPRT. So in addition to the endonuclease nicking preference, the mRNA of these elements appears to play an important role in integration site choice through a passive, post-nicking, selective process.

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.

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.

Progress in research on Entamoeba histolytica pathogenesis

Entamoeba histolytica is a protozoan parasite of humans that causes 40,000-100,000 deaths annually. Clinical amoebiasis results from the spread of the normally luminal parasite into the colon wall and beyond; the key development in understanding this complex multistage process has been the publication of the E. histolytica genome, from which has come an explosion in the use of microarrays to examine changes in gene expression that result from changes in growth conditions. The genome has also revealed a unique arrangement of tRNA genes and an extraordinary number of genes for putative virulence factors, many unexpressed under the artificial conditions of growth in culture. The ability to induce apoptosis of mammalian cells and a useful, but as yet little-understood, technique for epigenetic irreversible gene silencing are other exciting developments.

Transcriptional silencing of multiple genes in trophozoites of Entamoeba histolytica

In a previous work we described the transcriptional silencing of the amoebapore A (AP-A) gene (Ehap-a) of Entamoeba histolytica strain HM-1:IMSS. The silencing occurred following transfection with a plasmid containing a 5′ upstream region (473 bp) of Ehap-a that included a truncated segment (140 bp) of a short interspersed nuclear element (SINE1). Silencing remained in effect even after removal of the plasmid (clone G3). Neither short interfering RNA nor methylated DNA were detected, but the chromatin domain of Ehap-a in the gene-silenced trophozoites was modified. Two other similar genes (Ehap-b and one encoding a Saposin-like protein, SAPLIP 1) also became silenced. In the present work we demonstrate the silencing of a second gene of choice, one that encodes the light subunit of the Gal/GalNAc inhibitable lectin (Ehlgl1) and the other, the cysteine proteinase 5 (EhCP-5). This silencing occurred in G3 trophozoites transfected with a plasmid in which the 473 bp 5′ upstream Ehap-a fragment was directly ligated to the second gene. Transcriptional silencing occurred in both the transgene and the chromosomal gene. SINE1 sequences were essential, as was a direct connection between the Ehap-a upstream region and the beginning of the open reading frame of the second gene. Gene silencing did not occur in strain HM-1:IMSS with any of these plasmid constructs. The trophozoites with two silenced genes were virulence-attenuated as were those of clone G3. In addition, trophozoites not expressing Lgl1 and AP-A proteins had a significantly reduced ability to cap the Gal/GalNAc-lectin to the uroid region when incubated with antibodies against the heavy (170 kDa) subunit of the lectin. Lysates of trophozoites lacking cysteine proteinase 5 and AP-A proteins had 30% less cysteine proteinase activity than those of HM-1:IMSS strain or the G3 clone. Silencing of other genes in G3 amoebae could provide a model to study their various functions. In addition, double gene-silenced, virulence-attenuated trophozoites may be an important tool in vaccine development.

The human intestinal parasite Entamoeba histolytica has numerous genes that code for virulence. Silencing the expression of individual genes is useful to determine their roles. In previous work the authors demonstrated the silencing of the gene coding for amoebapore, which is responsible for killing of human cells. They transfected amoebic trophozoites with a plasmid that contained DNA sequences homologous to the promoter region of the amoebapore gene, as well as a portion of a repetitive DNA element (called a short interspersed nuclear element). This construct induced a modification of the chromatin and prevented the expression of the gene. Removal of the plasmid resulted in stable, amoebapore-deficient parasites possessing low virulence. In the present work, Bracha and colleagues show silencing of additional genes following transfection of E. histolytica trophozoites already silenced in amoebapore with a plasmid containing the second gene directly ligated to the upstream region of the amoebapore gene. The DNA sequences that are essential for transferring the silencing from the plasmid to the chromosomal gene copy were identified. Additional virulence genes that the authors irreversibly silenced are those encoding a subunit of a surface lectin that mediates the adherence of the parasite to host cells, and a cysteine proteinase that plays a role in inflammation and invasion of the intestine.

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.

Parasite histories and novel phylogenetic tools: Alternative approaches to inferring parasite evolution from molecular markers

Parasitological research is often contingent on the knowledge of the phylogeny/genealogy of the studied group. Although molecular phylogenetics has proved to be a powerful tool in such investigations, its application in the traditional fashion, based on a tree inference from the primary nucleotide sequences may, in many cases, be insufficient or even improper. These limitations are due to a number of factors, such as a scarcity/ambiguity of phylogenetic information in the sequences, an intricacy of gene relationships at low phylogenetic levels, or a lack of criteria when deciding among several competing coevolutionary scenarios. With respect to the importance of a precise and reliable phylogenetic background in many biological studies, attempts are being made to extend molecular phylogenetics with a variety of new data sources and methodologies. In this review, selected approaches potentially applicable to parasitological research are presented and their advantages as well as drawbacks are discussed. These issues include the usage of idiosyncratic markers (unique features with presumably low probability of homoplasy), such as insertion of mobile elements, gene rearrangements and secondary structure features; the problem of ancestral polymorphism and reticulate relationships at low phylogenetic levels; and the utility of a molecular clock to facilitate discrimination among alternative scenarios in host-parasite coevolution.

An extraordinary retrotransposon family encoding dual endonucleases

Retrotransposons commonly encode a reverse transcriptase (RT), but other functional domains are variable. The acquisition of new domains is the dominant evolutionary force that brings structural variety to retrotransposons. Non-long-terminal-repeat (non-LTR) retrotransposons are classified into two groups by their structure. Early branched non-LTR retrotransposons encode a restriction-like endonuclease (RLE), and recently branched non-LTR retrotransposons encode an apurinic/apyrimidinic endonuclease-like endonuclease (APE). In this study, we report a novel non-LTR retrotransposon family Dualen, identified from the Chlamydomonas reinhardtii genome. Dualen encodes two endonucleases, RLE and APE, with RT, ribonuclease H, and cysteine protease. Phylogenetic analyses of the RT domains revealed that Dualen is positioned at the midpoint between the early-branched and the recently branched groups. In the APE tree, Dualen was branched earlier than the I group and the Jockey group. The ribonuclease H domains among the Dualen family and other non-LTR retrotransposons are monophyletic. Phylogenies of three domains revealed the monophyly of the Dualen family members. The domain structure and the phylogeny of each domain imply that Dualen is a retrotransposon conserving the domain structure just after the acquisition of APE. From these observations, we discuss the evolution of domain structure of non-LTR retrotransposons.

LINEs, SINEs and processed pseudogenes: parasitic strategies for genome modeling

Two major classes of retrotransposons have invaded eukaryotic genomes: the LTR retrotransposons closely resembling the proviral integrated form of infectious retroviruses, and the non-LTR retrotransposons including the widespread, autonomous LINE elements. Here, we review the modeling effects of the latter class of elements, which are the most active in humans, and whose enzymatic machinery is subverted to generate a large series of "secondary" retroelements. These include the processed pseudogenes, naturally present in all eukaryotic genomes possessing non-LTR retroelements, and the very successful SINE elements such as the human Alu sequences which have evolved refined parasitic strategies to efficiently bypass the original "protectionist" cis-preference of LINEs for their own retrotransposition.

Unexpected Diversity and Differential Success of DNA Transposons in Four Species of Entamoeba Protozoans

We report the first comprehensive analysis of transposable element content in the compact genomes (approximately 20 Mb) of four species of Entamoeba unicellular protozoans for which draft sequences are now available. Entamoeba histolytica and Entamoeba dispar, two human parasites, have many retrotransposons, but few DNA transposons. In contrast, the reptile parasite Entamoeba invadens and the free-living Entamoeba moshkovskii contain few long interspersed elements but harbor diverse and recently amplified populations of DNA transposons. Representatives of three DNA transposase superfamilies (hobo/Activator/Tam3, Mutator, and piggyBac) were identified for the first time in a protozoan species in addition to a variety of members of a fourth superfamily (Tc1/mariner), previously reported only from ciliates and Trichomonas vaginalis among protozoans. The diversity of DNA transposons and their differential amplification among closely related species with similar compact genomes are discussed in the context of the biology of Entamoeba protozoans.

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

Autonomous non-long terminal repeat retrotransposons are commonly referred to as long interspersed elements (LINEs). Short non-autonomous elements that borrow the LINE machinery are called SINES. The Entamoeba histolytica genome contains three classes of LINEs and SINEs. Together the EhLINEs/SINEs account for about 6% of the genome. The recognizable functional domains in all three EhLINEs included reverse transcriptase and endonuclease. A novel feature was the presence of two types of members-some with a single long ORF (less frequent) and some with two ORFs (more frequent) in both EhLINE1 and 2. The two ORFs were generated by conserved changes leading to stop codon. Computational analysis of the immediate flanking sequences for each element showed that they inserted in AT-rich sequences, with a preponderance of Ts in the upstream site. The elements were very frequently located close to protein-coding genes and other EhLINEs/SINEs. The possible influence of these elements on expression of neighboring genes needs to be determined.