The Esterase and PHD Domains in CR1-Like Non-LTR Retrotransposons
Genetic Information Research Institute, Mountain View, Californiav.*w|0, 百拇医药
Abstractv.*w|0, 百拇医药
Most active non-LTR (long terminal repeat) retrotransposons carry two open reading frames (ORFs) encoding ORF1p and ORF2p proteins. The ORF2p proteins are relatively well studied and are known to contain endonuclease/reverse transcriptase domains. At the same time, the biological function of ORF1p proteins remains poorly understood, except in that they nonspecifically bind single-stranded mRNA/DNA molecules. CR1-like elements form the most widely distributed clade/superfamily of non-LTR retrotransposons. We found that ORF1p proteins encoded by diverse CR1-like elements contain conserved esterase domain (ES) or plant homeodomain (PHD). This indicates that CR1-like ORF1p proteins are either lipolytic enzymes or are involved in protein-protein interactions related to chromatin remodeling. Sequence conservation of ES suggests that interaction with cellular membranes is an important phase in life circles of CR1-like elements. Presumably such interaction helps in penetrating host cells. As a consequence, the presence of multiple young CR1 families characterized by ~ 10% intrafamily and 40% interfamily identities may be explained by a relatively frequent horizontal transfer of these CR1-like elements. Unexpectedly, ES links together non-LTR retrotransposons and single-stranded RNA viruses like influenza C and coronaviruses, which are known to depend on their own ES.
Key Words: non-LTR retrotransposon • CR1 clade • ORF1p • esterase • PHD homeodomain9^6, 百拇医药
Introduction9^6, 百拇医药
Genomes of all known eukaryotes are populated by transposable elements (TEs) capable of intragenomic multiplication or transposition . For example, recognizable fossils of TEs constitute approximately 45% and 12% of the Homo sapiens and Arabidospis thaliana genomes, respectively. Eukaryotic TEs can be divided into the following four classes: endogenous retroviruses and long terminal repeat (LTR) retrotransposons , non-LTR retrotransposons, including so-called LINEs, SINEs, and processed pseudogenes cut-and-paste DNA transposons , and rolling-circle DNA transposons . Duplication of a retrotransposon depends on reverse transcription and endonucleolytic cleavage, both of which are catalyzed by a reverse transcriptase (RT) and endonuclease domains of a polyprotein encoded by itself or by other retrotransposons. Primed by an endonucleolytic nick at the host DNA, an mRNA molecule, expressed during transcription of the retrotransposon DNA, is reverse transcribed and inserted in the genome. At present, known non-LTR retrotransposons belong to ~ 10 superfamilies or clades identified on the basis of phylogenetic studies of their protein sequences . Non-LTR retrotransposons form a clade if they share a common ancestor that is not shared by any other non-LTR retrotransposons outside the clade. Most non-LTR retrotransposons carry two long open reading frames, ORF1 and ORF2, which encode ORF1p and ORF2p proteins, respectively. ORF2p includes the RT, apurinic/apyrimidinic endonuclease (APE) or restriction-enzyme-like endonuclease domains. In some retroelements, ORF2p also includes a ribonuclease H domain. Whereas both the structure and the function of ORF2p are relatively well understood, properties of ORF1p remain obscure, in part because of the lack of significant similarity between ORF1p and proteins with known enzymatic functions. To date, the only structural elements discovered in different ORF1p proteins are the nonspecific zinc finger, leucine zipper, and coiled coil motifs . Experimental data also suggest that ORF1p proteins from the L1 and I clades bind single-stranded RNA–DNA. Overall, the role of ORF1 proteins in non-LTR retrotransposons is uncertain, although there are indications linking ORF1p to retroviral nucleocapsid proteins involved in packaging retroviral RNA and in other important steps of a retroviral "life cycle" .
CR1 is one of the most abundant and widely distributed clades of non-LTR retrotransposons . Most CR1 elements are severely truncated at their 5' ends. Therefore, it was found only recently that they are non-LTR retrotransposons populating genomes of birds, amphibians, and fishes ; lizards and turtles ; mammals ; and invertebrates .@l, 百拇医药
In this paper we report new full-length CR1-like elements from zebrafish, medaka, and fruit fly. We show that ORF1-encoded proteins in various CR1-like non-LTR retrotransposons include conserved plant homeodomain (PHD) and esterase domains. Given the conservation of the PHD and esterase domains in highly divergent CR1-like retrotransposons from different species, including those split several hundred million years ago, we assume that the PHD and esterase activities of the ORF1-encoded proteins were necessary for survival of these retrotransposons. Interestingly, as for the CR1-like non-LTR retrotransposons, the life cycle of enveloped negative-stranded and positive-stranded RNA viruses in birds and mammals depends on their own esterase.
Materials and Methods#v;/f, 百拇医药
All non-LTR retroelements reported here were found by using various methods of computational analysis. Starting with a compilation of known transposable elements collected in Repbase Update at , we identified their copies in DNA sequences deposited in GenBank. The identification began with comparing DNA and protein sequences of known TEs against the sequenced portion of the Danio rerio and Drosophila melanogaster genomes by using CENSOR and TBLASTN , respectively.#v;/f, 百拇医药
Using the majority rule applied to the corresponding set of multiple aligned copies of retrotransposons, we built their consensus sequences. Copies of TEs not produced directly by transpositions, such as those created by chromosomal duplications or redundant sequencing, were discarded based on the similarities between their flanking regions.#v;/f, 百拇医药
Distantly related proteins were identified using PsiBLAST . Multiple alignments of protein sequences were created by CLUSTAL-W . Alignments of DNA sequences were performed using the VMALN2 and PALN2 programs developed at the Genetic Information research Institute, Mountain View, California (GIRI). Phylogenetic analysis was conducted using MEGA 2.1 . Protein domains described in this article were identified using the Family Pairwise Search (FPS) algorithm and the SUPERFAMILY protein assignments server ; . Scoring of the protein sequences by FPS was performed against Pfam, a collection of protein family alignments reconstructed using hidden Markov models . Assignment by SUPERFAMILY has been performed using SCOPE, a library of protein superfamilies .
Sequences of retrotransposons reported here were deposited in the Repbase Update in the sections designated for fruit fly, zebrafish, human, vertebrates and invertebrates.ap, 百拇医药
Resultsap, 百拇医药
CR1 Elements from Zebrafishap, 百拇医药
Screening reverse transcriptase-like sequences in the publicly available DNA sequences, covering 1% of the D. rerio genome revealed multiple copies of non-LTR elements that belong to the CR1 clade. Cluster analysis of these sequences indicates that CR1-like elements in the D. rerio genome belong to over 10 young and diverse families characterized by 5% intrafamily and 35% interfamily nucleotide divergence (unpublished data. We assembled three consensus sequences that belong to families named CR1-1_DR, CR1-2_DR, and CR1-3_DR .ap, 百拇医药
fig.ommittedap, 百拇医药
FIG. 1. Schematic structure of complete CR1-like retrotransposons from fishes and insects. CR1-1_DR, CR1-2_DR and CR1-3_DR—are consensus sequences of retrotransposons that belong to the three families of retrotransposons identified in the Danio rerio genome. Maui and Rex1 are the consensus sequences of two retrotransposons from CR1-like families present in the Fugu rubripes genome. CR1_OL is a slightly damaged element identified in the Oryzias latipes genome. Horizontally shaded boxes mark ORF1s and ORF2s. ORF2s encode proteins composed of the apurinic/apyrimidinic endonuclease (APE) and reverse transcriptase (RT) domains. Proteins encoded by ORF1s are composed of putative zinc finger/leucine zipper (ZL) motifs, the plant homeodomain (PHD) and the esterase (ES) domains. Black squares, diamonds and hexagons indicate different unclassified domains. The 3' termini of all retrotransposons, excluding CR1-2_DR and CR1_DM, are shown starting from the polyadenylation signal, followed by terminal microsatellite repeats composed of different 4–7-bp units repeated 2–8 times. The average number of the repetitions is shown as a subscript index
The 4985-bp CR1-1_DR consensus sequence was built from 10 copies and it is ~ 98% identical with them. One partially truncated CR1-1_DR copy has been reported recently as the CR1DR1 element . The CR1-1_DR consensus sequence harbors two ORFs encoding 300-aa CR1-1DR1p and 1000-aa CR1-1DR2p proteins .cn, 百拇医药
The second family of zebrafish CR1-like elements is represented by a 4238-bp CR1-2_DR consensus sequence assembled from another 10 copies that are also ~ 98% identical to the consensus. Originally, a 600-bp fragment, 98% identical to a portion of the CR1-2_DR consensus sequence (positions 3422–4062), was reported as a LINE element . Recently, a 2900-bp CR1-2_DR copy was deposited in Repbase Update as the CR1DR2 element , which is 98% identical to the coding region of the CR1-2_DR consensus sequence (positions 1111–4008). Surprisingly, the CR1-2_DR consensus sequence includes a long 1110-bp 5' UTR region corresponding to ORF1 in various CR1-like elements. The only ORF in CR1-2_DR encodes a 965-aa CR1-2_DRp protein composed of the APE and RT domains.
Finally, the 5047-bp CR1-3_DR consensus sequence was built from seven different copies; they are ~ 95% identical with the consensus. CR1-3_DR carries ORF1 (positions 254-1801 ) and ORF2 (positions 1805–4717) encoding, respectively, a 516-aa CR1-3_DR1p protein and a 971-aa CR1-3_DR2p protein. As expected, CR1-3_DR2p is composed of the APE and RT domains.k, 百拇医药
CR1 Retrotransposon from Fruit Flyk, 百拇医药
A TBLASTN-based screening of CR1-like reverse transcriptases encoded by the Drosophila melanogaster genome revealed a rather abundant family of CR1-like non-LTR retrotransposons, hereafter named CR1_DM. The 4470-bp CR1_DM consensus sequence was constructed from 20 copies that were 10% divergent from one another. It contains ORF1 and ORF2 respectively encoding a 355-aa CR1_DM1p protein and a 964-aa CR1_DM2 protein . Approximately 100 copies of CR1_DM are present in a sequenced portion of the D. melanogaster genome that covers mainly euchromatin regions representing 70% of the genome. Multiple subfamilies of CR1_DM are present in the genome (unpublished data).
CR1 Elements from Medaka and Blood Flukel8)', 百拇医药
We also characterized a full-length 4985-bp copy of a CR1-like element in the Oryzias latipes (medaka fish) genome, called CR1-1_OL . Its 5' and 3' boundaries (GenBank, positions 330–5314) are labeled by 00-bp direct repeats composed of an 18-bp minisatellite unit. The CR1-1_OL element has been inserted into the genome relatively recently. Its ORF1 encodes a 271-aa CR1-1_OL1p protein (positions 322-1134), and ORF2 (positions 1352–4221) is corrupted by only two false frame shifts and one false stop codon.l8)', 百拇医药
Using genome survey sequences (GSS) from GenBank, we built a 3032-bp consensus sequence of the Schistosoma mansoni SR1 retrotransposon that is 700-bp longer than its sequence reported previously . The consensus sequence encodes a 950-aa SR1p protein (positions 36–2885) composed of the APE and RT domains. The extended region encodes APE. Available sequence data do not permit obtaining of any further 5'-extension of the SR1 consensus sequence, and we cannot prove or disprove the existence of ORF1p encoded by this element.
Diversity of the 3' Tails&y*ixw\, 百拇医药
Studies of DNA sequences flanking CR1-like elements presented in this article have revealed characteristics similar to those of known CR1-like elements reported previously (Haas et al. 1997; Kajikawa, Ohshima, and Okada 1997; Poulter, Butler, and Ormandy 1999). These elements do not generate target site duplications, and their 3' tails are composed of microsatellites . It appears that different families of CR1-elements, even those that populate the same genome, are characterized by different microsatellites that are specific for the each family. For example, the 3' termini of CR1-1_DR elements are composed of (ATTGA)n which follows GCTTGA and the polyadenylation signal. The 3' termini of CR1-2_DRs contain (AAATGT)n and they do not have any polyadenylation signal. In contrast, the 3' termini of CR1-3_DR elements are composed of the polyadenylation signal followed by (CTTGC)n.&y*ixw\, 百拇医药
It has not yet been proved, however, whether 3' microsatellite tails of CR1-like retrotransposons are their real termini or genomic microsatellites that served as targets during insertions of the retrotransposons. To resolve this question, we identified several CR1-like elements inserted into copies of other known TEs (unpublished data) that do not contain the microsatellites at positions targeted by the insertions. This observation suggests that the 3' microsatellites have been inserted into the genome together with CR1-like elements, and they can be considered to be distinctive hallmarks or signatures of different families. Presumably, these signatures depend on slightly different family-specific enzymatic activities encoded by the CR1-like elements. It is likely that generation of microsatellites at the 3' ends of CR1-like elements is a result of nontemplated additions by CR1-like reverse transcriptases, as shown experimentally for the I and R2 non-LTR retrotransposons .
Phylogenetic Analysishzq7+/, http://www.100md.com
shows a phylogeny of ORF2p proteins encoded by CR1-like non-LTR retrotransposons and several other elements that belong to non-CR1 clades. The phylogenetic analysis strongly suggests that the CR1 clade is composed of three major subclades.hzq7+/, http://www.100md.com
fig.ommittedhzq7+/, http://www.100md.com
FIG. 2. Phylogeny of the CR1-like non-LTR retrotransposons based on their endonuclease and reverse transcriptase domains. The phylogenetic tree also includes several retrotransposons from the Jockey, LOA, I, and L1 clades. Numbers next to each node indicate bootstrap values calculated as percentages of similar topologies out of 1,000 replicas for the neighbor-joining method. The names of non-LTR retrotransposons families and their host species are shown adjacent to the tree nodes. A scale of distances between the protein sequences is indicated. Solid triangles denote retrotransposons whose ORF1s code for the esterase. GenBank proteins identification numbers are as follows: Jockey (134083), Juan-C (1079026), Doc (8823), Lian (7511795), I (903726), CR1_BF (17529698), CR1 (2331059), CR1_PS (6576738), Q (11359829), T1 (159644), L1 (2072977). Sequences of the remaining retrotransposons have been deposited in the following sections of Repbase Update: humrep.ref (L2 and L3), dmrep.ref (CR1_DM, BAGGINS1, IVK), fugrep.ref (Maui, REX1), zebrep.ref (CR1-1_DR, CR1-2_DR, CR1-3_DR) and invrep.ref (SR1)
The first CR1 subclade, called CR1-I, includes CR1 and CR1_PS from chicken and turtle, L3 from mammals, and SR1 and CR1_BF from blood fluke and lancelet, respectively. Given the tree topology, distances, and bootstrap values, it is highly likely that, pending additional sequence data, this subclade will be split into at least three minor subclades. If so, SR1 and CR1_BF represent the potential minor subclades.7k*](t0, 百拇医药
The second major CR1 subclade, called CR1-II, includes CR1-elements identified in insects only: Q and T1 from the African malaria mosquito, and CR1_DM from the fruit fly genome. Actually, T1 element was the first element from the CR1 clade identified as a non-LTR retrotransposon . CR1 replaced T1 as the name of the clade after classification introduced by .7k*](t0, 百拇医药
Finally, the third major subclade, called CR1_III, includes L2 from mammals, Maui from pufferfish, and the CR1-1_DR, CR1-2_DR, and CR1-3_DR families from zebrafish. Interestingly, L2 and CR1-2_DR form a distinctive group separated from the other members of the third subclade by the 100% bootstrap value, and neither L2 nor CR1-2_DR encodes ORF1p-like proteins. In addition, the REX1 element may also represent a major subclade. It was suggested recently that L2-like and REX1-like elements form two novel clades of non-LTR retrotransposons called L2 and REX1, respectively . We do not think that the introduction of these two clades is strongly supported by available data. Moreover, we report here (see below) that the esterase domain is conserved in ORF1p encoded by different elements that belong to the CR1-I and CR1-III (L2) major subclades .
The PHD Domainz'$(e(}, 百拇医药
Computational analysis of the OFR1 proteins encoded by the zebrafish CR1-1_DR, CR1-2_DR, and CR1-3_DR elements failed to identify any zinc finger/leucine zipper motifs (ZL) similar to those present in the CR1, CR1_PS, and Maui retrotransposons from the chicken, turtle, and pufferfish genomes, respectively . However, ORF1p in CR1_OL from medaka fish harbors one motif distantly similar to ZL . Because stop codons and frame shifts that distort ORF2 encoded by the only available CR1_OL copy are present, it is likely that the originally intact ZL has also been damaged by mutations.z'$(e(}, 百拇医药
fig.ommittedz'$(e(}, 百拇医药
FIG. 3. Zinc finger motifs in proteins encoded by ORF1s of different CR1-like non-LTR retrotransposons. C denotes cysteine; L, leucine; H, histidine; X, any residue; the subscript index indicates the number of the amino acid residues marked by it. A, Putative zinc finger/leucine zipper domains in CR1 , CR1_PS (), Maui and CR1_OL from chicken, turtle, pufferfish, and medaka, respectively. B, ORF1 proteins encoded by the fruit fly CR1_DM and the malaria mosquito Q1 and T non-LTR retrotransposons that harbor the PHD domain. Conserved cysteine and histidine residues matching the PHD consensus sequence are highlighted. Numbers at the beginning and the end of the amino acid sequences indicate positions of the corresponding amino acid residues in the protein sequences deposited in GenBank and Repbase Update
Surprisingly, N-terminal portions of the ORF1 proteins encoded by the fly CR1_DM, and mosquito Q and T1 elements include motifs that fit the consensus sequence of a unique zinc finger domain called the PHD (plant homeodomain) or LAP (leukemia-associated protein) domain . The PHD domain has been identified in proteins primarily associated with chromatin and involved in chromatin-mediated transcription control . As indicated by a PsiBLAST search, all three ORF1p proteins from CR1_DM, T1, and Q are similar to one another over a span of ~ 300 aa, although the overall sequence identity is only ~ 20%. Additionally, presence of the PHD domain in the ORF1 proteins encoded by CR1_DM, T1, and Q was confirmed by other computational methods. Given an E-value < 0.01, PHD was the only domain identified in these proteins by the Family Pairwise Search algorithm and SUPERFAMILY . Therefore, the presence of PHD in the highly divergent N-terminal portions of the ORF1p proteins encoded by insect CR1-like elements should be considered to be a strong indication that this domain is important for their life cycle.
The exact function of the PHD domain is not yet known, but it is thought to be involved in protein–protein interactions and to be of importance for the assembly or activity of multicomponent complexes involved in transcriptional activation or repression . Multiple lines of evidence suggest that PHD domain proteins can be targeted to DNA only indirectly via protein–protein interactions . Therefore, it is unlikely that the zinc fingers encoded by ORF1s in CR1-like elements are involved directly in DNA or RNA binding, as proposed earlier for the putative zinc finger/leucine zipper domains in CR1_PS .es+, http://www.100md.com
The Esterase Domaines+, http://www.100md.com
On the basis of a BLASTP search, we identified only three GenBank proteins similar to CR1-3_DR1p (E < 0.01). They are the ORF1 proteins from the chicken CR1 , turtle CR1_PS , and pufferfish Maui retrotransposons. Because only a central portion of CR1-3_DR1p (positions 168–327) is similar to the ORF1p proteins (22% to 32% identity), we used it as a separate query for a PsiBLAST search (E < 0.005). After several iterations, the central portion converged with ~ 150 eukaryotic and prokaryotic proteins from the esterase/acetylhydrolase superfamily . The same classification of CR1-3_DRp1 (E < 10-10) was also supported by the SUPERFAMILY genome assignments server . shows a multiple alignment of ORF1p proteins from CR1 retrotransposons and several prokaryotic and eukaryotic esterases. Two esterases included in the alignment were comprehensively studied experimentally: PAF-AH, a brain acetylhydrolase from cow , and RGAE, a rhamnogalacturonan acetylesterase from fungi . The most conserved structural hallmark of esterases is a catalytic triad composed of properly arranged serine, histidine, and aspartic acid residues . Different order and spacing of amino acid residues from the catalytic triad define several families of esterases ). Presumably, the esterase domain (ES) encoded by the CR1-like ORF1 proteins belongs to a specific family called GDSL , SGNH , or the rhamnogalacturonan acetylesterase family . This family is characterized by GDS, GXND, and DXXH conserved motifs . It has been shown experimentally that serine from the first motif and aspartic acid plus histidine from the third motif belong to the catalytic triad. Strikingly, all three motifs and the catalytic triad are perfectly conserved in the highly divergent ORF1p encoded by CR1-elements from the chicken, turtle, medaka, pufferfish, and zebrafish genomes . The alignment also includes ES found in the ancient L3 retrotransposon fossilized in the human genome (see next section). Additionally, we found ES conserved in putative ORF1p proteins encoded by CR1-like elements fossilized in the crocodile, frog, and salmon genomes (unpublished data).
fig.ommittedq?, 百拇医药
FIG. 4. Multiple sequence alignment of the putative conserved esterase domains encoded by ORF1s in CR1-like non-LTR retrotransposons and other esterases. Solid arrowheads mark the catalytic serine-asparate-histidine triad. Ambiguous amino acids are denoted by Xs. GenBank protein identification numbers are as follows: CR1_PS (6576737), CR1 (2331058), Maui (4378024), NeuA (13876786, CMP-N-acetylneuramic acid synthetase from Streptococcus agalactiae), TesA (267107, Acyl-CoA thioesterase I from Escherichia coli), RGAE (7766904), rhamnogalacturonan acetylesterase from Aspergillus aculeatus), PAF-AH (2624421, platelet-activating factor acetylhydrolase from Bos taurus). Amino acid sequences of ORF1p proteins encoded by CR1-1_DR, CR1-2_DR, CR1-3_DR, CR1-1_TN, CR1_OL and L3 are deposited in Repbase Updateq?, 百拇医药
L3, the Most Ancient Transposable Element Ever Reconstructed In Silicoq?, 百拇医药
Using TBLASTN, we found that the 169-aa CR1_PS esterase domain matches eight proteins (10-10 " E " 10-2) encoded by human GenBank sequences . Apparently none of these proteins are functional because of multiple stop codons. DNA sequences encoding these proteins were extracted from the corresponding GenBank sequences. Based on a multiple alignment of the extracted fragments, a 700-bp consensus sequence was assembled, which was ~ 73% identical to the fragments. Remarkably, a 120-aa protein encoded by the consensus sequence was much more similar to the CR1_PS esterase domain than any of the eight ORFs (E = 10-33, 52% identity). Moreover, the consensus sequence ORF was not interrupted by stop codons. Applying the BLASTN search with the consensus sequence as a query, we identified 13 fragments in the human genome similar to the consensus, including the original eight fragments. However, screening of assembled human chromosomes by CENSOR revealed ~ 300 copies similar to the consensus sequence. Each copy has been expanded up to 7 kb in both directions. After pairwise alignment of the expanded sequences with each other, we eliminated chromosomal duplications more than 80% identical to each other. Subsequently, the final set of 220 sequences was screened for known TEs using CENSOR and Repbase Update. As we expected, there was a striking correlation between the extracted esterase-like fragments and L3. The CR1-like esterase domain was followed by remnants of L3 in 53 sequences, and both the esterase and L3 were in the same orientation. This is consistent with the fact that the L3 element is an ancient CR1-like non-LTR retrotransposon whose reverse transcriptase is closest to CR1_PS2p. However, L3 is so old that only its relatively short 3' portion was recovered previously as the ~ 1.8-kb L3 consensus sequence . Using the set of 220 masked sequences, we iteratively built three consensus sequences that represent missing parts of the ancient L3 . They encode domains closest to those present in CR1_PS, including the very beginning of ORF1p, esterase, endonuclease, and a middle portion of ORF2p. These domains are also present in rodents and other mammals (unpublished data).
fig.ommittedh%-'*, http://www.100md.com
FIG. 5. Flowchart of the identification of ORF1p encoded by the ancient L3 retrotransposon fossilized in the human genome. Arrows indicate information flow directions. Rectangles illustrate different computational processes indicated by corresponding program names. Parallelograms indicate specific sets of data. Cans symbolize databases. GenBank accession numbers of sequences containing DNA regions which encode protein sequences TBLASTN-similar to CR1_PS1p are indicated together with corresponding E values. The "smiley face" marks the esterase domain found in different CR1-like elementsh%-'*, http://www.100md.com
The L3 consensus sequence reported here is composed of four separate segments ~ 65% identical to L3 copies. Given this diversity and the age (over 200 Myr), we can recover only the most conserved parts of L3s. Our data indicate that the esterase domain was functional in CR1-like retrotransposons that had been multiplied in ancestors of all mammals. It should be pointed out that many L3 copies are interrupted by L2 elements inserted randomly at different positions (unpublished data). Therefore, L3 elements are older than L2 elements. As the oldest transposable element identified in the human genome, L3 can be an extremely useful reference sequence in evolutionary studies.
Discussiontz19, http://www.100md.com
Despite the abundance and variety of non-LTR retrotransposons, many aspects of their life cycles and evolution are not known. Overall, they are viewed quite mechanistically as genomic parasites vertically transmitted during evolution of eukaryotic genomes ). Our functional understanding of non-LTR retrotransposons is definitely lagging behind what we know about endogenous retroviruses/LTR retrotransposons (1997). In particular, we know little about the function of proteins encoded by ORF1 in different families of non-LTR retrotransposons. Our discovery of two specific and conserved domains, PHD (plant homeodomain) and ES (esterase, a lipolytic acetylhydrolase) may represent a breakthrough in this respect. Based on a classification by SCOP , the esterase/acetylhydrolase superfamily is composed of four families: (I) esterase, (II) esterase domain of hemagglutinin glycoprotein HEF1, (III) acetylhydrolase, and (IV) rhamnogalacturonan acetylesterase. Presumably, families III and IV can be joined together as the so-called GDSL or SGNH families . GDSL includes secreted and outer membrane–bound esterases, acetylhydrolases, and arylesterases . Usually, these enzymes remove acetyl or fatty acids from complex polysaccharides, viral glycoproteins, and cellular proteins interacting with membranes and involved in cell signaling or the regulation of the immune system.
Esterase is important for the life cycles of enveloped negative-stranded and positive-stranded RNA viruses infecting birds and mammals. For example, esterase domains are included in membrane glycoproteins, so called hemagglutinin-esterase or the HEF1 proteins, which are present on the surfaces of influenza C and coronaviruses and toroviruses ). During the infection process, hemagglutinins interact with sialic acid molecules bound to the cell receptors. This interaction is followed by entrance into the cell of virus particles that cannot be efficient unless ester bonds formed between the hemagglutinin glycoproteins and the receptor sialic acids are cut by esterase .6, 百拇医药
It is known that esterases perform enzymatic depalmitoylation of viral glycoproteins and various cellular proteins. As a result, fatty acids (usually palmitate) covalently attached to cysteines near C-termini of palmitoylated proteins are cleaved off. It is thought that palmitoylation can affect a protein's affinity for membranes, subcellular localization, and interactions with membrane proteins. Rhamnogalacturonan acetylesterase from fungi (RGAE, ) catalyzes degradation of polysaccharides that constitute a cell wall in the plant host .
The esterase domain is conserved in the highly divergent ORF1p proteins of CR1-like elements from the chicken, turtle, fish, and human genomes (and putatively from the frog and crocodile genomes). This underscores its functional importance for life cycles of CR1-like elements. Surprisingly, its function is not linked directly to any known stages of a non-LTR retrotransposon life cycle. It is really difficult to understand why the esterase was preserved by non-LTR retrotransposons whose evolution is thought to follow, usually, a "vertical transmission" model . However, a regular horizontal transfer/transmission of CR1-like elements would favor esterases involved in penetration of cell membranes. Interestingly, the chicken and zebrafish genomes harbor multiple CR1-like families of approximately the same age. Six of them have been identified in the chicken genome . Three families of CR1-like elements from the zebrafish are reported in this article. All three have been retrotransposed relatively recently because of a low ~ 5% to 10% nucleotide divergence between elements that belong to the same family. However, there is an enormous ~ 40% divergence between elements that belong to any of different families residing in the same genome. It is conceivable that these families have invaded the host independently, and most of their diversity was acquired in some other hosts.
The PHD domain is another specific domain that we identified in the ORF1p proteins encoded by the CR1_DM, T1, and Q non-LTR retrotransposons from fruit fly and African malaria mosquito, respectively . These elements form the only well-defined CR1 sub-clade whose members do not code for esterase . As for esterase, the PHD domain is conserved in highly divergent proteins, and its function is not related to DNA/RNA binding. The PHD domain is thought to be involved in protein–protein interactions related to chromatin remodeling . Therefore, it is possible that the PHD domain in CR1-like retrotransposons is necessary for both efficient retrotransposition and minimization of potentially harmful insertions of retrotransposons into the host genome by providing dynamic regulatory feedback between chromatin structure, expression of reverse transcriptase/endonuclease by retrotransposons, and their target-specificity. Interestingly, T1 and Q elements are most abundant in paracentromeric heterochromatin . Similar abundance of different TEs in paracentromeric heterochromatin has been observed in other species . It is possible that most of the TEs inserted accidentally into paracentromeric heterochromatin were fixed, whereas most of their relatives inserted originally into euchromatin have been lost. It is also possible, however, that insertion of some TEs can be channeled to heterochromatin regions by PHD-like regulatory elements, which may suppress transcription of retrotransposons at stages when most of the euchromatin is open . It is striking that some gypsy-like LTR retrotransposons have acquired "chromodomain" which, like to PHD, is also involved in chromatin remodeling .
Interestingly, the PHD domain was acquired by Kaposi's sarcoma–associated herpesvirus . The N-terminal PHD domain in MIR proteins encoded by the herpesvirus is directly involved in recruiting cellular proteins that regulate endocytosis of host immune recognition proteins . As for the herpesvirus, CR1-like elements might have recruited the PHD domain to evade the host defense. Such evasion may be potentially important if these elements regularly trespass host cells.y33, 百拇医药
One may design other interesting models employing function of PHD and ES in ORF1p proteins. However, our main goal was to identify new domains in the ORF1p proteins and to underscore the complexity of the life cycle of non-LTR retrotransposons concealed by the popular "vertical transmission" model.y33, 百拇医药
Acknowledgementsy33, 百拇医药
We are grateful to Jolanta Walichiewicz and Michael Jurka for help with illustrations, and for editing the manuscript. We thank reviewers of the manuscript for useful comments. This work was supported by National Institutes of Health grant 2 P41 LM06252-04A1.
Literature Citedw.v[j/j, 百拇医药
Aasland, R., T. J. Gibson, and A. F. Stewart. 1995. The PHD finger: implications for chromatin-mediated transcriptional regulation. Trends Biochem. Sci 20:56-59.w.v[j/j, 百拇医药
Aasland, R., and A. F. Stewart. 1995. The chromo shadow domain, a second chromo domain in heterochromatin-binding protein 1, HP1. Nucleic Acids Res 23:3168-3174.w.v[j/j, 百拇医药
Altschul, S. F., T. L. Madden, A. A. Schaffer, J. Zhang, Z. Zhang, W. Miller, and D. J. Lipman. 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25:3389-3402.w.v[j/j, 百拇医药
Arpigny, J. L., and K. E. Jaeger. 1999. Bacterial lipolytic enzymes: classification and properties. Biochem. J 343:177-183.w.v[j/j, 百拇医药
Bailey, T. L., and W. N. Grundy. 1999. Classifying proteins by family using the product of correlated p-values, pp. 10–14. in P. Istrail, P. Pevzner, and M. Waterman, eds. Proceedings of the Third International Conference on Computational Molecular Biology (RECOMB99). ACM, New York.w.v[j/j, 百拇医药
Bateman, A., E. Birney, L. Cerruti, R. Durbin, L. Etwiller, S. R. Eddy, S. Griffiths-Jones, K. L. Howe, M. Marshall, and E. L. Sonnhammer. 2002. The Pfam protein families database. Nucleic Acids Res 30:276-280.r?y+?/, http://www.100md.com
Berg, D. E., and M. H. Howe. 1987. Mobile DNA. American Society for Microbiology Press, Washington, DC.r?y+?/, http://www.100md.com
Besansky, N. J. 1990. A retrotransposable element from the mosquito Anopheles gambiae. Mol. Cell. Biol 10:863-871.r?y+?/, http://www.100md.com
Besansky, N. J., J. A. Bedell, and O. Mukabayire. 1994. Q: a new retrotransposon from the mosquito Anopheles gambiae. Insect Mol. Biol 3:49-56.r?y+?/, http://www.100md.com
Burch, J. B., D. L. Davis, and N. B. Haas. 1993. Chicken repeat 1 elements contain a pol-like open reading frame and belong to the non-long terminal repeat class of retrotransposons. Proc. Natl. Acad. Sci. USA 90:8199-8203.r?y+?/, http://www.100md.com
Capili, A. D., D. C. Schultz, I. F. Rauscher, and K. L. Borden. 2001. Solution structure of the PHD domain from the KAP-1 corepressor: structural determinants for PHD, RING and LIM zinc-binding domains. EMBO J 20:165-177.
Capy, P., C. Bazin, D. Higuet, and T. Langin. 1998. Dynamics and evolution of transposable elements. Chapman & Hall, New York.|/0sm, 百拇医药
Chaboissier, M. C., D. Finnegan, and A. Bucheton. 2000. Retrotransposition of the I factor, a non-long terminal repeat retrotransposon of Drosophila, generates tandem repeats at the 3' end. Nucleic Acids Res 28:2467-2472.|/0sm, 百拇医药
Coffin, J. M., S. H. Hughes, and H. E. Varmus. 1997. Retroviruses. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, New York.|/0sm, 百拇医药
Coscoy, L., D. J. Sanchez, and D. Ganem. 2001. A novel class of herpesvirus-encoded membrane-bound E3 ubiquitin ligases regulates endocytosis of proteins involved in immune recognition. J. Cell Biol 155:1265-1273.|/0sm, 百拇医药
Craig, N. L. 1995. Unity in transposition reactions. Science 270:253-254.|/0sm, 百拇医药
Dalrymple, B. P., D. H. Cybinski, I. Layton, C. S. McSweeney, G. P. Xue, Y. J. Swadling, and J. B. Lowry. 1997. Three Neocallimastix patriciarum esterases associated with the degradation of complex polysaccharides are members of a new family of hydrolases. Microbiology 143:2605-2614.
Dawson, A., E. Hartswood, T. Paterson, and D. J. Finnegan. 1997. A LINE-like transposable element in Drosophila, the I factor, encodes a protein with properties similar to those of retroviral nucleocapsids. EMBO J 16:4448-4455.*5f2$, 百拇医药
Drablos, F., and S. B. Petersen. 1997. Identification of conserved residues in family of esterase and lipase sequences. Methods Enzymol 284:28-61.*5f2$, 百拇医药
Drew, A. C., and P. J. Brindley. 1997. A retrotransposon of the non-long terminal repeat class from the human blood fluke Schistosoma mansoni. Similarities to the chicken-repeat-1-like elements of vertebrates. Mol. Biol. Evol 14:602-610.*5f2$, 百拇医药
Dunphy, J. T., and M. E. Linder. 1998. Signalling functions of protein palmitoylation. Biochim. Biophys. Acta 1436:245-261.*5f2$, 百拇医药
Eickbush, D. G., D. D. Luan, and T. H. Eickbush. 2000. Integration of Bombyx mori R2 sequences into the 28S ribosomal RNA genes of Drosophila melanogaster. Mol. Cell. Biol 20:213-223.*5f2$, 百拇医药
Gibbons, R. J., S. Bachoo, D. J. Picketts, et al 1997. Mutations in transcriptional regulator ATRX establish the functional significance of a PHD-like domain. Nat. Genet 17:146-148.
Gough, J., K. Karplus, R. Hughey, and C. Chothia. 2001. Assignment of homology to genome sequences using a library of hidden Markov models that represent all proteins of known structure. J. Mol. Biol 313:903-919.#ls@, 百拇医药
Haas, N. B., J. M. Grabowski, J. North, J. V. Moran, H. H. Kazazian, and J. B. Burch. 2001. Subfamilies of CR1 non-LTR retrotransposons have different 5' UTR sequences but are otherwise conserved. Gene 265:175-183.#ls@, 百拇医药
Haas, N. B., J. M. Grabowski, A. B. Sivitz, and J. B. Burch. 1997. Chicken repeat 1 (CR1) elements, which define an ancient family of vertebrate non-LTR retrotransposons, contain two closely spaced open reading frames. Gene 197:305-309.#ls@, 百拇医药
Herrler, G., R. Rott, H. D. Klenk, H. P. Muller, A. K. Shukla, and R. Schauer. 1985. The receptor-destroying enzyme of influenza C virus is neuraminate-O-acetylesterase. EMBO J 4:1503-1506.#ls@, 百拇医药
Ho, Y. S., L. Swenson, U. Derewenda, et al 1997. Brain acetylhydrolase that inactivates platelet-activating factor is a G-protein-like trimer. Nature 385:89-93.
Hohjoh, H., and M. F. Singer. 1997. Sequence-specific single-strand RNA binding protein encoded by the human LINE-1 retrotransposon. EMBO J 16:6034-6043.(x!q, 百拇医药
Holmes, S. E., M. F. Singer, and G. D. Swergold. 1992. Studies on p40, the leucine zipper motif-containing protein encoded by the first open reading frame of an active human LINE-1 transposable element. J. Biol. Chem 267:19765-19768.(x!q, 百拇医药
Jekosch, K. 2002. CR1-like repeat from Danio rerio. Repbase Reports 2:7-8(x!q, 百拇医药
Jurka, J. 2000. Repbase Update: a database and an electronic journal of repetitive elements. Trends Genet 16:418-420.(x!q, 百拇医药
Jurka, J., and V. V. Kapitonov. 1999a. L3, humrep, Repbase Update .(x!q, 百拇医药
Jurka, J., and 1999b. Sectorial mutagenesis by transposable elements. Genetica 107:239-248.(x!q, 百拇医药
Jurka, J., P. Klonowski, V. Dagman, and P. Pelton. 1996. CENSOR—a program for identification and elimination of repetitive elements from DNA sequences. Comput. Chem 20:119-121.(x!q, 百拇医药
Kajikawa, M., K. Ohshima, and N. Okada. 1997. Determination of the entire sequence of turtle CR1: the first open reading frame of the turtle CR1 element encodes a protein with a novel zinc finger motif. Mol. Biol. Evol 14:1206-1217.1, 百拇医药
Kapitonov, V. V., and J. Jurka. 1999. Molecular paleontology of transposable elements from Arabidopsis thaliana. Genetica 107:27-37.1, 百拇医药
Kapitonov, V. V., and 2001. Rolling-circle transposons in eukaryotes. Proc. Natl. Acad. Sci. USA 98:8714-8719.1, 百拇医药
Kehle, J., D. Beuchle, S. Treuheit, B. Christen, J. A. Kennison, M. Bienz, and J. Muller. 1998. dMi-2, a hunchback-interacting protein that functions in polycomb repression. Science 282:1897-1900.1, 百拇医药
Koipally, J., A. Renold, J. Kim, and K. Georgopoulos. 1999. Repression by Ikaros and Aiolos is mediated through histone deacetylase complexes. EMBO J 18:3090-3100.1, 百拇医药
Kolosha, V. O., and S. L. Martin. 1997. In vitro properties of the first ORF protein from mouse LINE-1 support its role in ribonucleoprotein particle formation during retrotransposition. Proc. Natl. Acad. Sci. USA 94:10155-10160.
Kumar, S., K. Tamura, I. B. Jakobsen, and M. Nei. 2001. MEGA2: molecular evolutionary genetics analysis software. Bioinformatics 17:1244-1245.%t)}@@, 百拇医药
Lander, E. S., L. M. Linton, B. Birren, et al 2001. Initial sequencing and analysis of the human genome. Nature 409:860-921.%t)}@@, 百拇医药
Lovsin, N., F. Gubensek, and D. Kordi. 2001. Evolutionary dynamics in a novel L2 clade of non-LTR retrotransposons in Deuterostomia. Mol. Biol. Evol 18:2213-2224.%t)}@@, 百拇医药
Lyngso, C., G. Bouteiller, C. K. Damgaard, D. Ryom, S. Sanchez-Munoz, P. L. Norby, B. J. Bonven, and P. Jorgensen. 2000. Interaction between the transcription factor SPBP and the positive cofactor RNF4. An interplay between protein binding zinc fingers. J. Biol. Chem 275:26144-26149.%t)}@@, 百拇医药
Malik, H. S., W. D. Burke, and T. H. Eickbush. 1999. The age and evolution of non-LTR retrotransposable elements. Mol. Biol. Evol 16:793-805.%t)}@@, 百拇医药
Malik, H. S., and T. H. Eickbush. 1999. Modular evolution of the integrase domain in the Ty3/Gypsy class of LTR retrotransposons. J. Virol 73:5186-5190.
Martin, S. L., and F. D. Bushman. 2001. Nucleic acid chaperone activity of the ORF1 protein from the mouse LINE-1 retrotransposon. Mol. Cell. Biol 21:467-475.3j\e/i#, 百拇医药
Molgaard, A., S. Kauppinen, and S. Larsen. 2000. Rhamnogalacturonan acetylesterase elucidates the structureand function of a new family of hydrolases. Structure Fold Des 8:373-383.3j\e/i#, 百拇医药
Mukabayire, O., and N. J. Besansky. 1996. Distribution of T1, Q, Pegasus and mariner transposable elements on the polytene chromosomes of PEST, a standard strain of Anopheles gambiae. Chromosoma 104:585-595.3j\e/i#, 百拇医药
Murzin, A. G., S. E. Brenner, T. Hubbard, and C. Chothia. 1995. SCOP: a structural classification of proteins database for the investigation of sequences and structures. J. Mol. Biol 247:536-540.3j\e/i#, 百拇医药
Okada, N., M. Hamada, I. Ogiwara, and K. Ohshima. 1997. SINEs and LINEs share common 3' sequences: a review. Gene 205:229-243.3j\e/i#, 百拇医药
Poulter, R., M. Butler, and J. Ormandy. 1999. A LINE element from the pufferfish (fugu) Fugu rubripes which shows similarity to the CR1 family of non-LTR retrotransposons. Gene 227:169-179.
Rosenthal, P. B., X. Zhang, F. Formanowski, W. Fitz, C. H. Wong, H. Meier-Ewert, J. J. Skehel, and D. C. Wiley. 1998. Structure of the haemagglutinin-esterase-fusion glycoprotein of influenza C virus. Nature 396:92-96.s/-c?, http://www.100md.com
Saha, V., T. Chaplin, A. Gregorini, P. Ayton, and B. D. Young. 1995. The leukemia-associated-protein (LAP) domain, a cysteine-rich motif, is present in a wide range of proteins, including MLL, AF10, and MLLT6 proteins. Proc. Natl. Acad. Sci. USA 92:9737-9741.s/-c?, http://www.100md.com
Smit, A. F. 1996. The origin of interspersed repeats in the human genome. Curr. Opin. Genet. Dev 6:743-748.s/-c?, http://www.100md.com
Smit, A. F. 2000. L3, humrep. Repbase Update .s/-c?, http://www.100md.com
Smit, A. F. 2001. REX1_FURC. Repbase Update (fugrep.ref) .s/-c?, http://www.100md.com
Thompson, J. D., D. G. Higgins, and T. J. Gibson. 1994. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22:4673-4680.s/-c?, http://www.100md.com
Vandergon, T. L., and M. Reitman. 1994. Evolution of chicken repeat 1 (CR1) elements: evidence for ancient subfamilies and multiple progenitors. Mol. Biol. Evol 11:886-898.s/-c?, http://www.100md.com
Volff, J. N., C. Korting, and M. Schartl. 2000. Multiple lineages of the non-LTR retrotransposon Rex1 with varying success in invading fish genomes. Mol. Biol. Evol 17:1673-1684.s/-c?, http://www.100md.com
Weiner, A. M. 2000. Do all SINEs lead to LINEs?. Nat. Genet 24:332-333.s/-c?, http://www.100md.com
Wurzer, W. J., K. Obojes, and R. Vlasak. 2002. The sialate-4-O-acetylesterases of coronaviruses related to mouse hepatitis virus: a proposal to reorganize group 2 Coronaviridae. J. Gen. Virol 83:395-402.s/-c?, http://www.100md.com
Yochum, G. S., and D. E. Ayer. 2001. Pf1, a novel PHD zinc finger protein that links the TLE corepressor to the mSin3A-histone deacetylase complex. Mol. Cell. Biol 21:4110-4118.s/-c?, http://www.100md.com
Accepted for publication September 2, 2002.(Vladimir V. Kapitonov and Jerzy Jurka)
Abstractv.*w|0, 百拇医药
Most active non-LTR (long terminal repeat) retrotransposons carry two open reading frames (ORFs) encoding ORF1p and ORF2p proteins. The ORF2p proteins are relatively well studied and are known to contain endonuclease/reverse transcriptase domains. At the same time, the biological function of ORF1p proteins remains poorly understood, except in that they nonspecifically bind single-stranded mRNA/DNA molecules. CR1-like elements form the most widely distributed clade/superfamily of non-LTR retrotransposons. We found that ORF1p proteins encoded by diverse CR1-like elements contain conserved esterase domain (ES) or plant homeodomain (PHD). This indicates that CR1-like ORF1p proteins are either lipolytic enzymes or are involved in protein-protein interactions related to chromatin remodeling. Sequence conservation of ES suggests that interaction with cellular membranes is an important phase in life circles of CR1-like elements. Presumably such interaction helps in penetrating host cells. As a consequence, the presence of multiple young CR1 families characterized by ~ 10% intrafamily and 40% interfamily identities may be explained by a relatively frequent horizontal transfer of these CR1-like elements. Unexpectedly, ES links together non-LTR retrotransposons and single-stranded RNA viruses like influenza C and coronaviruses, which are known to depend on their own ES.
Key Words: non-LTR retrotransposon • CR1 clade • ORF1p • esterase • PHD homeodomain9^6, 百拇医药
Introduction9^6, 百拇医药
Genomes of all known eukaryotes are populated by transposable elements (TEs) capable of intragenomic multiplication or transposition . For example, recognizable fossils of TEs constitute approximately 45% and 12% of the Homo sapiens and Arabidospis thaliana genomes, respectively. Eukaryotic TEs can be divided into the following four classes: endogenous retroviruses and long terminal repeat (LTR) retrotransposons , non-LTR retrotransposons, including so-called LINEs, SINEs, and processed pseudogenes cut-and-paste DNA transposons , and rolling-circle DNA transposons . Duplication of a retrotransposon depends on reverse transcription and endonucleolytic cleavage, both of which are catalyzed by a reverse transcriptase (RT) and endonuclease domains of a polyprotein encoded by itself or by other retrotransposons. Primed by an endonucleolytic nick at the host DNA, an mRNA molecule, expressed during transcription of the retrotransposon DNA, is reverse transcribed and inserted in the genome. At present, known non-LTR retrotransposons belong to ~ 10 superfamilies or clades identified on the basis of phylogenetic studies of their protein sequences . Non-LTR retrotransposons form a clade if they share a common ancestor that is not shared by any other non-LTR retrotransposons outside the clade. Most non-LTR retrotransposons carry two long open reading frames, ORF1 and ORF2, which encode ORF1p and ORF2p proteins, respectively. ORF2p includes the RT, apurinic/apyrimidinic endonuclease (APE) or restriction-enzyme-like endonuclease domains. In some retroelements, ORF2p also includes a ribonuclease H domain. Whereas both the structure and the function of ORF2p are relatively well understood, properties of ORF1p remain obscure, in part because of the lack of significant similarity between ORF1p and proteins with known enzymatic functions. To date, the only structural elements discovered in different ORF1p proteins are the nonspecific zinc finger, leucine zipper, and coiled coil motifs . Experimental data also suggest that ORF1p proteins from the L1 and I clades bind single-stranded RNA–DNA. Overall, the role of ORF1 proteins in non-LTR retrotransposons is uncertain, although there are indications linking ORF1p to retroviral nucleocapsid proteins involved in packaging retroviral RNA and in other important steps of a retroviral "life cycle" .
CR1 is one of the most abundant and widely distributed clades of non-LTR retrotransposons . Most CR1 elements are severely truncated at their 5' ends. Therefore, it was found only recently that they are non-LTR retrotransposons populating genomes of birds, amphibians, and fishes ; lizards and turtles ; mammals ; and invertebrates .@l, 百拇医药
In this paper we report new full-length CR1-like elements from zebrafish, medaka, and fruit fly. We show that ORF1-encoded proteins in various CR1-like non-LTR retrotransposons include conserved plant homeodomain (PHD) and esterase domains. Given the conservation of the PHD and esterase domains in highly divergent CR1-like retrotransposons from different species, including those split several hundred million years ago, we assume that the PHD and esterase activities of the ORF1-encoded proteins were necessary for survival of these retrotransposons. Interestingly, as for the CR1-like non-LTR retrotransposons, the life cycle of enveloped negative-stranded and positive-stranded RNA viruses in birds and mammals depends on their own esterase.
Materials and Methods#v;/f, 百拇医药
All non-LTR retroelements reported here were found by using various methods of computational analysis. Starting with a compilation of known transposable elements collected in Repbase Update at , we identified their copies in DNA sequences deposited in GenBank. The identification began with comparing DNA and protein sequences of known TEs against the sequenced portion of the Danio rerio and Drosophila melanogaster genomes by using CENSOR and TBLASTN , respectively.#v;/f, 百拇医药
Using the majority rule applied to the corresponding set of multiple aligned copies of retrotransposons, we built their consensus sequences. Copies of TEs not produced directly by transpositions, such as those created by chromosomal duplications or redundant sequencing, were discarded based on the similarities between their flanking regions.#v;/f, 百拇医药
Distantly related proteins were identified using PsiBLAST . Multiple alignments of protein sequences were created by CLUSTAL-W . Alignments of DNA sequences were performed using the VMALN2 and PALN2 programs developed at the Genetic Information research Institute, Mountain View, California (GIRI). Phylogenetic analysis was conducted using MEGA 2.1 . Protein domains described in this article were identified using the Family Pairwise Search (FPS) algorithm and the SUPERFAMILY protein assignments server ; . Scoring of the protein sequences by FPS was performed against Pfam, a collection of protein family alignments reconstructed using hidden Markov models . Assignment by SUPERFAMILY has been performed using SCOPE, a library of protein superfamilies .
Sequences of retrotransposons reported here were deposited in the Repbase Update in the sections designated for fruit fly, zebrafish, human, vertebrates and invertebrates.ap, 百拇医药
Resultsap, 百拇医药
CR1 Elements from Zebrafishap, 百拇医药
Screening reverse transcriptase-like sequences in the publicly available DNA sequences, covering 1% of the D. rerio genome revealed multiple copies of non-LTR elements that belong to the CR1 clade. Cluster analysis of these sequences indicates that CR1-like elements in the D. rerio genome belong to over 10 young and diverse families characterized by 5% intrafamily and 35% interfamily nucleotide divergence (unpublished data. We assembled three consensus sequences that belong to families named CR1-1_DR, CR1-2_DR, and CR1-3_DR .ap, 百拇医药
fig.ommittedap, 百拇医药
FIG. 1. Schematic structure of complete CR1-like retrotransposons from fishes and insects. CR1-1_DR, CR1-2_DR and CR1-3_DR—are consensus sequences of retrotransposons that belong to the three families of retrotransposons identified in the Danio rerio genome. Maui and Rex1 are the consensus sequences of two retrotransposons from CR1-like families present in the Fugu rubripes genome. CR1_OL is a slightly damaged element identified in the Oryzias latipes genome. Horizontally shaded boxes mark ORF1s and ORF2s. ORF2s encode proteins composed of the apurinic/apyrimidinic endonuclease (APE) and reverse transcriptase (RT) domains. Proteins encoded by ORF1s are composed of putative zinc finger/leucine zipper (ZL) motifs, the plant homeodomain (PHD) and the esterase (ES) domains. Black squares, diamonds and hexagons indicate different unclassified domains. The 3' termini of all retrotransposons, excluding CR1-2_DR and CR1_DM, are shown starting from the polyadenylation signal, followed by terminal microsatellite repeats composed of different 4–7-bp units repeated 2–8 times. The average number of the repetitions is shown as a subscript index
The 4985-bp CR1-1_DR consensus sequence was built from 10 copies and it is ~ 98% identical with them. One partially truncated CR1-1_DR copy has been reported recently as the CR1DR1 element . The CR1-1_DR consensus sequence harbors two ORFs encoding 300-aa CR1-1DR1p and 1000-aa CR1-1DR2p proteins .cn, 百拇医药
The second family of zebrafish CR1-like elements is represented by a 4238-bp CR1-2_DR consensus sequence assembled from another 10 copies that are also ~ 98% identical to the consensus. Originally, a 600-bp fragment, 98% identical to a portion of the CR1-2_DR consensus sequence (positions 3422–4062), was reported as a LINE element . Recently, a 2900-bp CR1-2_DR copy was deposited in Repbase Update as the CR1DR2 element , which is 98% identical to the coding region of the CR1-2_DR consensus sequence (positions 1111–4008). Surprisingly, the CR1-2_DR consensus sequence includes a long 1110-bp 5' UTR region corresponding to ORF1 in various CR1-like elements. The only ORF in CR1-2_DR encodes a 965-aa CR1-2_DRp protein composed of the APE and RT domains.
Finally, the 5047-bp CR1-3_DR consensus sequence was built from seven different copies; they are ~ 95% identical with the consensus. CR1-3_DR carries ORF1 (positions 254-1801 ) and ORF2 (positions 1805–4717) encoding, respectively, a 516-aa CR1-3_DR1p protein and a 971-aa CR1-3_DR2p protein. As expected, CR1-3_DR2p is composed of the APE and RT domains.k, 百拇医药
CR1 Retrotransposon from Fruit Flyk, 百拇医药
A TBLASTN-based screening of CR1-like reverse transcriptases encoded by the Drosophila melanogaster genome revealed a rather abundant family of CR1-like non-LTR retrotransposons, hereafter named CR1_DM. The 4470-bp CR1_DM consensus sequence was constructed from 20 copies that were 10% divergent from one another. It contains ORF1 and ORF2 respectively encoding a 355-aa CR1_DM1p protein and a 964-aa CR1_DM2 protein . Approximately 100 copies of CR1_DM are present in a sequenced portion of the D. melanogaster genome that covers mainly euchromatin regions representing 70% of the genome. Multiple subfamilies of CR1_DM are present in the genome (unpublished data).
CR1 Elements from Medaka and Blood Flukel8)', 百拇医药
We also characterized a full-length 4985-bp copy of a CR1-like element in the Oryzias latipes (medaka fish) genome, called CR1-1_OL . Its 5' and 3' boundaries (GenBank, positions 330–5314) are labeled by 00-bp direct repeats composed of an 18-bp minisatellite unit. The CR1-1_OL element has been inserted into the genome relatively recently. Its ORF1 encodes a 271-aa CR1-1_OL1p protein (positions 322-1134), and ORF2 (positions 1352–4221) is corrupted by only two false frame shifts and one false stop codon.l8)', 百拇医药
Using genome survey sequences (GSS) from GenBank, we built a 3032-bp consensus sequence of the Schistosoma mansoni SR1 retrotransposon that is 700-bp longer than its sequence reported previously . The consensus sequence encodes a 950-aa SR1p protein (positions 36–2885) composed of the APE and RT domains. The extended region encodes APE. Available sequence data do not permit obtaining of any further 5'-extension of the SR1 consensus sequence, and we cannot prove or disprove the existence of ORF1p encoded by this element.
Diversity of the 3' Tails&y*ixw\, 百拇医药
Studies of DNA sequences flanking CR1-like elements presented in this article have revealed characteristics similar to those of known CR1-like elements reported previously (Haas et al. 1997; Kajikawa, Ohshima, and Okada 1997; Poulter, Butler, and Ormandy 1999). These elements do not generate target site duplications, and their 3' tails are composed of microsatellites . It appears that different families of CR1-elements, even those that populate the same genome, are characterized by different microsatellites that are specific for the each family. For example, the 3' termini of CR1-1_DR elements are composed of (ATTGA)n which follows GCTTGA and the polyadenylation signal. The 3' termini of CR1-2_DRs contain (AAATGT)n and they do not have any polyadenylation signal. In contrast, the 3' termini of CR1-3_DR elements are composed of the polyadenylation signal followed by (CTTGC)n.&y*ixw\, 百拇医药
It has not yet been proved, however, whether 3' microsatellite tails of CR1-like retrotransposons are their real termini or genomic microsatellites that served as targets during insertions of the retrotransposons. To resolve this question, we identified several CR1-like elements inserted into copies of other known TEs (unpublished data) that do not contain the microsatellites at positions targeted by the insertions. This observation suggests that the 3' microsatellites have been inserted into the genome together with CR1-like elements, and they can be considered to be distinctive hallmarks or signatures of different families. Presumably, these signatures depend on slightly different family-specific enzymatic activities encoded by the CR1-like elements. It is likely that generation of microsatellites at the 3' ends of CR1-like elements is a result of nontemplated additions by CR1-like reverse transcriptases, as shown experimentally for the I and R2 non-LTR retrotransposons .
Phylogenetic Analysishzq7+/, http://www.100md.com
shows a phylogeny of ORF2p proteins encoded by CR1-like non-LTR retrotransposons and several other elements that belong to non-CR1 clades. The phylogenetic analysis strongly suggests that the CR1 clade is composed of three major subclades.hzq7+/, http://www.100md.com
fig.ommittedhzq7+/, http://www.100md.com
FIG. 2. Phylogeny of the CR1-like non-LTR retrotransposons based on their endonuclease and reverse transcriptase domains. The phylogenetic tree also includes several retrotransposons from the Jockey, LOA, I, and L1 clades. Numbers next to each node indicate bootstrap values calculated as percentages of similar topologies out of 1,000 replicas for the neighbor-joining method. The names of non-LTR retrotransposons families and their host species are shown adjacent to the tree nodes. A scale of distances between the protein sequences is indicated. Solid triangles denote retrotransposons whose ORF1s code for the esterase. GenBank proteins identification numbers are as follows: Jockey (134083), Juan-C (1079026), Doc (8823), Lian (7511795), I (903726), CR1_BF (17529698), CR1 (2331059), CR1_PS (6576738), Q (11359829), T1 (159644), L1 (2072977). Sequences of the remaining retrotransposons have been deposited in the following sections of Repbase Update: humrep.ref (L2 and L3), dmrep.ref (CR1_DM, BAGGINS1, IVK), fugrep.ref (Maui, REX1), zebrep.ref (CR1-1_DR, CR1-2_DR, CR1-3_DR) and invrep.ref (SR1)
The first CR1 subclade, called CR1-I, includes CR1 and CR1_PS from chicken and turtle, L3 from mammals, and SR1 and CR1_BF from blood fluke and lancelet, respectively. Given the tree topology, distances, and bootstrap values, it is highly likely that, pending additional sequence data, this subclade will be split into at least three minor subclades. If so, SR1 and CR1_BF represent the potential minor subclades.7k*](t0, 百拇医药
The second major CR1 subclade, called CR1-II, includes CR1-elements identified in insects only: Q and T1 from the African malaria mosquito, and CR1_DM from the fruit fly genome. Actually, T1 element was the first element from the CR1 clade identified as a non-LTR retrotransposon . CR1 replaced T1 as the name of the clade after classification introduced by .7k*](t0, 百拇医药
Finally, the third major subclade, called CR1_III, includes L2 from mammals, Maui from pufferfish, and the CR1-1_DR, CR1-2_DR, and CR1-3_DR families from zebrafish. Interestingly, L2 and CR1-2_DR form a distinctive group separated from the other members of the third subclade by the 100% bootstrap value, and neither L2 nor CR1-2_DR encodes ORF1p-like proteins. In addition, the REX1 element may also represent a major subclade. It was suggested recently that L2-like and REX1-like elements form two novel clades of non-LTR retrotransposons called L2 and REX1, respectively . We do not think that the introduction of these two clades is strongly supported by available data. Moreover, we report here (see below) that the esterase domain is conserved in ORF1p encoded by different elements that belong to the CR1-I and CR1-III (L2) major subclades .
The PHD Domainz'$(e(}, 百拇医药
Computational analysis of the OFR1 proteins encoded by the zebrafish CR1-1_DR, CR1-2_DR, and CR1-3_DR elements failed to identify any zinc finger/leucine zipper motifs (ZL) similar to those present in the CR1, CR1_PS, and Maui retrotransposons from the chicken, turtle, and pufferfish genomes, respectively . However, ORF1p in CR1_OL from medaka fish harbors one motif distantly similar to ZL . Because stop codons and frame shifts that distort ORF2 encoded by the only available CR1_OL copy are present, it is likely that the originally intact ZL has also been damaged by mutations.z'$(e(}, 百拇医药
fig.ommittedz'$(e(}, 百拇医药
FIG. 3. Zinc finger motifs in proteins encoded by ORF1s of different CR1-like non-LTR retrotransposons. C denotes cysteine; L, leucine; H, histidine; X, any residue; the subscript index indicates the number of the amino acid residues marked by it. A, Putative zinc finger/leucine zipper domains in CR1 , CR1_PS (), Maui and CR1_OL from chicken, turtle, pufferfish, and medaka, respectively. B, ORF1 proteins encoded by the fruit fly CR1_DM and the malaria mosquito Q1 and T non-LTR retrotransposons that harbor the PHD domain. Conserved cysteine and histidine residues matching the PHD consensus sequence are highlighted. Numbers at the beginning and the end of the amino acid sequences indicate positions of the corresponding amino acid residues in the protein sequences deposited in GenBank and Repbase Update
Surprisingly, N-terminal portions of the ORF1 proteins encoded by the fly CR1_DM, and mosquito Q and T1 elements include motifs that fit the consensus sequence of a unique zinc finger domain called the PHD (plant homeodomain) or LAP (leukemia-associated protein) domain . The PHD domain has been identified in proteins primarily associated with chromatin and involved in chromatin-mediated transcription control . As indicated by a PsiBLAST search, all three ORF1p proteins from CR1_DM, T1, and Q are similar to one another over a span of ~ 300 aa, although the overall sequence identity is only ~ 20%. Additionally, presence of the PHD domain in the ORF1 proteins encoded by CR1_DM, T1, and Q was confirmed by other computational methods. Given an E-value < 0.01, PHD was the only domain identified in these proteins by the Family Pairwise Search algorithm and SUPERFAMILY . Therefore, the presence of PHD in the highly divergent N-terminal portions of the ORF1p proteins encoded by insect CR1-like elements should be considered to be a strong indication that this domain is important for their life cycle.
The exact function of the PHD domain is not yet known, but it is thought to be involved in protein–protein interactions and to be of importance for the assembly or activity of multicomponent complexes involved in transcriptional activation or repression . Multiple lines of evidence suggest that PHD domain proteins can be targeted to DNA only indirectly via protein–protein interactions . Therefore, it is unlikely that the zinc fingers encoded by ORF1s in CR1-like elements are involved directly in DNA or RNA binding, as proposed earlier for the putative zinc finger/leucine zipper domains in CR1_PS .es+, http://www.100md.com
The Esterase Domaines+, http://www.100md.com
On the basis of a BLASTP search, we identified only three GenBank proteins similar to CR1-3_DR1p (E < 0.01). They are the ORF1 proteins from the chicken CR1 , turtle CR1_PS , and pufferfish Maui retrotransposons. Because only a central portion of CR1-3_DR1p (positions 168–327) is similar to the ORF1p proteins (22% to 32% identity), we used it as a separate query for a PsiBLAST search (E < 0.005). After several iterations, the central portion converged with ~ 150 eukaryotic and prokaryotic proteins from the esterase/acetylhydrolase superfamily . The same classification of CR1-3_DRp1 (E < 10-10) was also supported by the SUPERFAMILY genome assignments server . shows a multiple alignment of ORF1p proteins from CR1 retrotransposons and several prokaryotic and eukaryotic esterases. Two esterases included in the alignment were comprehensively studied experimentally: PAF-AH, a brain acetylhydrolase from cow , and RGAE, a rhamnogalacturonan acetylesterase from fungi . The most conserved structural hallmark of esterases is a catalytic triad composed of properly arranged serine, histidine, and aspartic acid residues . Different order and spacing of amino acid residues from the catalytic triad define several families of esterases ). Presumably, the esterase domain (ES) encoded by the CR1-like ORF1 proteins belongs to a specific family called GDSL , SGNH , or the rhamnogalacturonan acetylesterase family . This family is characterized by GDS, GXND, and DXXH conserved motifs . It has been shown experimentally that serine from the first motif and aspartic acid plus histidine from the third motif belong to the catalytic triad. Strikingly, all three motifs and the catalytic triad are perfectly conserved in the highly divergent ORF1p encoded by CR1-elements from the chicken, turtle, medaka, pufferfish, and zebrafish genomes . The alignment also includes ES found in the ancient L3 retrotransposon fossilized in the human genome (see next section). Additionally, we found ES conserved in putative ORF1p proteins encoded by CR1-like elements fossilized in the crocodile, frog, and salmon genomes (unpublished data).
fig.ommittedq?, 百拇医药
FIG. 4. Multiple sequence alignment of the putative conserved esterase domains encoded by ORF1s in CR1-like non-LTR retrotransposons and other esterases. Solid arrowheads mark the catalytic serine-asparate-histidine triad. Ambiguous amino acids are denoted by Xs. GenBank protein identification numbers are as follows: CR1_PS (6576737), CR1 (2331058), Maui (4378024), NeuA (13876786, CMP-N-acetylneuramic acid synthetase from Streptococcus agalactiae), TesA (267107, Acyl-CoA thioesterase I from Escherichia coli), RGAE (7766904), rhamnogalacturonan acetylesterase from Aspergillus aculeatus), PAF-AH (2624421, platelet-activating factor acetylhydrolase from Bos taurus). Amino acid sequences of ORF1p proteins encoded by CR1-1_DR, CR1-2_DR, CR1-3_DR, CR1-1_TN, CR1_OL and L3 are deposited in Repbase Updateq?, 百拇医药
L3, the Most Ancient Transposable Element Ever Reconstructed In Silicoq?, 百拇医药
Using TBLASTN, we found that the 169-aa CR1_PS esterase domain matches eight proteins (10-10 " E " 10-2) encoded by human GenBank sequences . Apparently none of these proteins are functional because of multiple stop codons. DNA sequences encoding these proteins were extracted from the corresponding GenBank sequences. Based on a multiple alignment of the extracted fragments, a 700-bp consensus sequence was assembled, which was ~ 73% identical to the fragments. Remarkably, a 120-aa protein encoded by the consensus sequence was much more similar to the CR1_PS esterase domain than any of the eight ORFs (E = 10-33, 52% identity). Moreover, the consensus sequence ORF was not interrupted by stop codons. Applying the BLASTN search with the consensus sequence as a query, we identified 13 fragments in the human genome similar to the consensus, including the original eight fragments. However, screening of assembled human chromosomes by CENSOR revealed ~ 300 copies similar to the consensus sequence. Each copy has been expanded up to 7 kb in both directions. After pairwise alignment of the expanded sequences with each other, we eliminated chromosomal duplications more than 80% identical to each other. Subsequently, the final set of 220 sequences was screened for known TEs using CENSOR and Repbase Update. As we expected, there was a striking correlation between the extracted esterase-like fragments and L3. The CR1-like esterase domain was followed by remnants of L3 in 53 sequences, and both the esterase and L3 were in the same orientation. This is consistent with the fact that the L3 element is an ancient CR1-like non-LTR retrotransposon whose reverse transcriptase is closest to CR1_PS2p. However, L3 is so old that only its relatively short 3' portion was recovered previously as the ~ 1.8-kb L3 consensus sequence . Using the set of 220 masked sequences, we iteratively built three consensus sequences that represent missing parts of the ancient L3 . They encode domains closest to those present in CR1_PS, including the very beginning of ORF1p, esterase, endonuclease, and a middle portion of ORF2p. These domains are also present in rodents and other mammals (unpublished data).
fig.ommittedh%-'*, http://www.100md.com
FIG. 5. Flowchart of the identification of ORF1p encoded by the ancient L3 retrotransposon fossilized in the human genome. Arrows indicate information flow directions. Rectangles illustrate different computational processes indicated by corresponding program names. Parallelograms indicate specific sets of data. Cans symbolize databases. GenBank accession numbers of sequences containing DNA regions which encode protein sequences TBLASTN-similar to CR1_PS1p are indicated together with corresponding E values. The "smiley face" marks the esterase domain found in different CR1-like elementsh%-'*, http://www.100md.com
The L3 consensus sequence reported here is composed of four separate segments ~ 65% identical to L3 copies. Given this diversity and the age (over 200 Myr), we can recover only the most conserved parts of L3s. Our data indicate that the esterase domain was functional in CR1-like retrotransposons that had been multiplied in ancestors of all mammals. It should be pointed out that many L3 copies are interrupted by L2 elements inserted randomly at different positions (unpublished data). Therefore, L3 elements are older than L2 elements. As the oldest transposable element identified in the human genome, L3 can be an extremely useful reference sequence in evolutionary studies.
Discussiontz19, http://www.100md.com
Despite the abundance and variety of non-LTR retrotransposons, many aspects of their life cycles and evolution are not known. Overall, they are viewed quite mechanistically as genomic parasites vertically transmitted during evolution of eukaryotic genomes ). Our functional understanding of non-LTR retrotransposons is definitely lagging behind what we know about endogenous retroviruses/LTR retrotransposons (1997). In particular, we know little about the function of proteins encoded by ORF1 in different families of non-LTR retrotransposons. Our discovery of two specific and conserved domains, PHD (plant homeodomain) and ES (esterase, a lipolytic acetylhydrolase) may represent a breakthrough in this respect. Based on a classification by SCOP , the esterase/acetylhydrolase superfamily is composed of four families: (I) esterase, (II) esterase domain of hemagglutinin glycoprotein HEF1, (III) acetylhydrolase, and (IV) rhamnogalacturonan acetylesterase. Presumably, families III and IV can be joined together as the so-called GDSL or SGNH families . GDSL includes secreted and outer membrane–bound esterases, acetylhydrolases, and arylesterases . Usually, these enzymes remove acetyl or fatty acids from complex polysaccharides, viral glycoproteins, and cellular proteins interacting with membranes and involved in cell signaling or the regulation of the immune system.
Esterase is important for the life cycles of enveloped negative-stranded and positive-stranded RNA viruses infecting birds and mammals. For example, esterase domains are included in membrane glycoproteins, so called hemagglutinin-esterase or the HEF1 proteins, which are present on the surfaces of influenza C and coronaviruses and toroviruses ). During the infection process, hemagglutinins interact with sialic acid molecules bound to the cell receptors. This interaction is followed by entrance into the cell of virus particles that cannot be efficient unless ester bonds formed between the hemagglutinin glycoproteins and the receptor sialic acids are cut by esterase .6, 百拇医药
It is known that esterases perform enzymatic depalmitoylation of viral glycoproteins and various cellular proteins. As a result, fatty acids (usually palmitate) covalently attached to cysteines near C-termini of palmitoylated proteins are cleaved off. It is thought that palmitoylation can affect a protein's affinity for membranes, subcellular localization, and interactions with membrane proteins. Rhamnogalacturonan acetylesterase from fungi (RGAE, ) catalyzes degradation of polysaccharides that constitute a cell wall in the plant host .
The esterase domain is conserved in the highly divergent ORF1p proteins of CR1-like elements from the chicken, turtle, fish, and human genomes (and putatively from the frog and crocodile genomes). This underscores its functional importance for life cycles of CR1-like elements. Surprisingly, its function is not linked directly to any known stages of a non-LTR retrotransposon life cycle. It is really difficult to understand why the esterase was preserved by non-LTR retrotransposons whose evolution is thought to follow, usually, a "vertical transmission" model . However, a regular horizontal transfer/transmission of CR1-like elements would favor esterases involved in penetration of cell membranes. Interestingly, the chicken and zebrafish genomes harbor multiple CR1-like families of approximately the same age. Six of them have been identified in the chicken genome . Three families of CR1-like elements from the zebrafish are reported in this article. All three have been retrotransposed relatively recently because of a low ~ 5% to 10% nucleotide divergence between elements that belong to the same family. However, there is an enormous ~ 40% divergence between elements that belong to any of different families residing in the same genome. It is conceivable that these families have invaded the host independently, and most of their diversity was acquired in some other hosts.
The PHD domain is another specific domain that we identified in the ORF1p proteins encoded by the CR1_DM, T1, and Q non-LTR retrotransposons from fruit fly and African malaria mosquito, respectively . These elements form the only well-defined CR1 sub-clade whose members do not code for esterase . As for esterase, the PHD domain is conserved in highly divergent proteins, and its function is not related to DNA/RNA binding. The PHD domain is thought to be involved in protein–protein interactions related to chromatin remodeling . Therefore, it is possible that the PHD domain in CR1-like retrotransposons is necessary for both efficient retrotransposition and minimization of potentially harmful insertions of retrotransposons into the host genome by providing dynamic regulatory feedback between chromatin structure, expression of reverse transcriptase/endonuclease by retrotransposons, and their target-specificity. Interestingly, T1 and Q elements are most abundant in paracentromeric heterochromatin . Similar abundance of different TEs in paracentromeric heterochromatin has been observed in other species . It is possible that most of the TEs inserted accidentally into paracentromeric heterochromatin were fixed, whereas most of their relatives inserted originally into euchromatin have been lost. It is also possible, however, that insertion of some TEs can be channeled to heterochromatin regions by PHD-like regulatory elements, which may suppress transcription of retrotransposons at stages when most of the euchromatin is open . It is striking that some gypsy-like LTR retrotransposons have acquired "chromodomain" which, like to PHD, is also involved in chromatin remodeling .
Interestingly, the PHD domain was acquired by Kaposi's sarcoma–associated herpesvirus . The N-terminal PHD domain in MIR proteins encoded by the herpesvirus is directly involved in recruiting cellular proteins that regulate endocytosis of host immune recognition proteins . As for the herpesvirus, CR1-like elements might have recruited the PHD domain to evade the host defense. Such evasion may be potentially important if these elements regularly trespass host cells.y33, 百拇医药
One may design other interesting models employing function of PHD and ES in ORF1p proteins. However, our main goal was to identify new domains in the ORF1p proteins and to underscore the complexity of the life cycle of non-LTR retrotransposons concealed by the popular "vertical transmission" model.y33, 百拇医药
Acknowledgementsy33, 百拇医药
We are grateful to Jolanta Walichiewicz and Michael Jurka for help with illustrations, and for editing the manuscript. We thank reviewers of the manuscript for useful comments. This work was supported by National Institutes of Health grant 2 P41 LM06252-04A1.
Literature Citedw.v[j/j, 百拇医药
Aasland, R., T. J. Gibson, and A. F. Stewart. 1995. The PHD finger: implications for chromatin-mediated transcriptional regulation. Trends Biochem. Sci 20:56-59.w.v[j/j, 百拇医药
Aasland, R., and A. F. Stewart. 1995. The chromo shadow domain, a second chromo domain in heterochromatin-binding protein 1, HP1. Nucleic Acids Res 23:3168-3174.w.v[j/j, 百拇医药
Altschul, S. F., T. L. Madden, A. A. Schaffer, J. Zhang, Z. Zhang, W. Miller, and D. J. Lipman. 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25:3389-3402.w.v[j/j, 百拇医药
Arpigny, J. L., and K. E. Jaeger. 1999. Bacterial lipolytic enzymes: classification and properties. Biochem. J 343:177-183.w.v[j/j, 百拇医药
Bailey, T. L., and W. N. Grundy. 1999. Classifying proteins by family using the product of correlated p-values, pp. 10–14. in P. Istrail, P. Pevzner, and M. Waterman, eds. Proceedings of the Third International Conference on Computational Molecular Biology (RECOMB99). ACM, New York.w.v[j/j, 百拇医药
Bateman, A., E. Birney, L. Cerruti, R. Durbin, L. Etwiller, S. R. Eddy, S. Griffiths-Jones, K. L. Howe, M. Marshall, and E. L. Sonnhammer. 2002. The Pfam protein families database. Nucleic Acids Res 30:276-280.r?y+?/, http://www.100md.com
Berg, D. E., and M. H. Howe. 1987. Mobile DNA. American Society for Microbiology Press, Washington, DC.r?y+?/, http://www.100md.com
Besansky, N. J. 1990. A retrotransposable element from the mosquito Anopheles gambiae. Mol. Cell. Biol 10:863-871.r?y+?/, http://www.100md.com
Besansky, N. J., J. A. Bedell, and O. Mukabayire. 1994. Q: a new retrotransposon from the mosquito Anopheles gambiae. Insect Mol. Biol 3:49-56.r?y+?/, http://www.100md.com
Burch, J. B., D. L. Davis, and N. B. Haas. 1993. Chicken repeat 1 elements contain a pol-like open reading frame and belong to the non-long terminal repeat class of retrotransposons. Proc. Natl. Acad. Sci. USA 90:8199-8203.r?y+?/, http://www.100md.com
Capili, A. D., D. C. Schultz, I. F. Rauscher, and K. L. Borden. 2001. Solution structure of the PHD domain from the KAP-1 corepressor: structural determinants for PHD, RING and LIM zinc-binding domains. EMBO J 20:165-177.
Capy, P., C. Bazin, D. Higuet, and T. Langin. 1998. Dynamics and evolution of transposable elements. Chapman & Hall, New York.|/0sm, 百拇医药
Chaboissier, M. C., D. Finnegan, and A. Bucheton. 2000. Retrotransposition of the I factor, a non-long terminal repeat retrotransposon of Drosophila, generates tandem repeats at the 3' end. Nucleic Acids Res 28:2467-2472.|/0sm, 百拇医药
Coffin, J. M., S. H. Hughes, and H. E. Varmus. 1997. Retroviruses. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, New York.|/0sm, 百拇医药
Coscoy, L., D. J. Sanchez, and D. Ganem. 2001. A novel class of herpesvirus-encoded membrane-bound E3 ubiquitin ligases regulates endocytosis of proteins involved in immune recognition. J. Cell Biol 155:1265-1273.|/0sm, 百拇医药
Craig, N. L. 1995. Unity in transposition reactions. Science 270:253-254.|/0sm, 百拇医药
Dalrymple, B. P., D. H. Cybinski, I. Layton, C. S. McSweeney, G. P. Xue, Y. J. Swadling, and J. B. Lowry. 1997. Three Neocallimastix patriciarum esterases associated with the degradation of complex polysaccharides are members of a new family of hydrolases. Microbiology 143:2605-2614.
Dawson, A., E. Hartswood, T. Paterson, and D. J. Finnegan. 1997. A LINE-like transposable element in Drosophila, the I factor, encodes a protein with properties similar to those of retroviral nucleocapsids. EMBO J 16:4448-4455.*5f2$, 百拇医药
Drablos, F., and S. B. Petersen. 1997. Identification of conserved residues in family of esterase and lipase sequences. Methods Enzymol 284:28-61.*5f2$, 百拇医药
Drew, A. C., and P. J. Brindley. 1997. A retrotransposon of the non-long terminal repeat class from the human blood fluke Schistosoma mansoni. Similarities to the chicken-repeat-1-like elements of vertebrates. Mol. Biol. Evol 14:602-610.*5f2$, 百拇医药
Dunphy, J. T., and M. E. Linder. 1998. Signalling functions of protein palmitoylation. Biochim. Biophys. Acta 1436:245-261.*5f2$, 百拇医药
Eickbush, D. G., D. D. Luan, and T. H. Eickbush. 2000. Integration of Bombyx mori R2 sequences into the 28S ribosomal RNA genes of Drosophila melanogaster. Mol. Cell. Biol 20:213-223.*5f2$, 百拇医药
Gibbons, R. J., S. Bachoo, D. J. Picketts, et al 1997. Mutations in transcriptional regulator ATRX establish the functional significance of a PHD-like domain. Nat. Genet 17:146-148.
Gough, J., K. Karplus, R. Hughey, and C. Chothia. 2001. Assignment of homology to genome sequences using a library of hidden Markov models that represent all proteins of known structure. J. Mol. Biol 313:903-919.#ls@, 百拇医药
Haas, N. B., J. M. Grabowski, J. North, J. V. Moran, H. H. Kazazian, and J. B. Burch. 2001. Subfamilies of CR1 non-LTR retrotransposons have different 5' UTR sequences but are otherwise conserved. Gene 265:175-183.#ls@, 百拇医药
Haas, N. B., J. M. Grabowski, A. B. Sivitz, and J. B. Burch. 1997. Chicken repeat 1 (CR1) elements, which define an ancient family of vertebrate non-LTR retrotransposons, contain two closely spaced open reading frames. Gene 197:305-309.#ls@, 百拇医药
Herrler, G., R. Rott, H. D. Klenk, H. P. Muller, A. K. Shukla, and R. Schauer. 1985. The receptor-destroying enzyme of influenza C virus is neuraminate-O-acetylesterase. EMBO J 4:1503-1506.#ls@, 百拇医药
Ho, Y. S., L. Swenson, U. Derewenda, et al 1997. Brain acetylhydrolase that inactivates platelet-activating factor is a G-protein-like trimer. Nature 385:89-93.
Hohjoh, H., and M. F. Singer. 1997. Sequence-specific single-strand RNA binding protein encoded by the human LINE-1 retrotransposon. EMBO J 16:6034-6043.(x!q, 百拇医药
Holmes, S. E., M. F. Singer, and G. D. Swergold. 1992. Studies on p40, the leucine zipper motif-containing protein encoded by the first open reading frame of an active human LINE-1 transposable element. J. Biol. Chem 267:19765-19768.(x!q, 百拇医药
Jekosch, K. 2002. CR1-like repeat from Danio rerio. Repbase Reports 2:7-8(x!q, 百拇医药
Jurka, J. 2000. Repbase Update: a database and an electronic journal of repetitive elements. Trends Genet 16:418-420.(x!q, 百拇医药
Jurka, J., and V. V. Kapitonov. 1999a. L3, humrep, Repbase Update .(x!q, 百拇医药
Jurka, J., and 1999b. Sectorial mutagenesis by transposable elements. Genetica 107:239-248.(x!q, 百拇医药
Jurka, J., P. Klonowski, V. Dagman, and P. Pelton. 1996. CENSOR—a program for identification and elimination of repetitive elements from DNA sequences. Comput. Chem 20:119-121.(x!q, 百拇医药
Kajikawa, M., K. Ohshima, and N. Okada. 1997. Determination of the entire sequence of turtle CR1: the first open reading frame of the turtle CR1 element encodes a protein with a novel zinc finger motif. Mol. Biol. Evol 14:1206-1217.1, 百拇医药
Kapitonov, V. V., and J. Jurka. 1999. Molecular paleontology of transposable elements from Arabidopsis thaliana. Genetica 107:27-37.1, 百拇医药
Kapitonov, V. V., and 2001. Rolling-circle transposons in eukaryotes. Proc. Natl. Acad. Sci. USA 98:8714-8719.1, 百拇医药
Kehle, J., D. Beuchle, S. Treuheit, B. Christen, J. A. Kennison, M. Bienz, and J. Muller. 1998. dMi-2, a hunchback-interacting protein that functions in polycomb repression. Science 282:1897-1900.1, 百拇医药
Koipally, J., A. Renold, J. Kim, and K. Georgopoulos. 1999. Repression by Ikaros and Aiolos is mediated through histone deacetylase complexes. EMBO J 18:3090-3100.1, 百拇医药
Kolosha, V. O., and S. L. Martin. 1997. In vitro properties of the first ORF protein from mouse LINE-1 support its role in ribonucleoprotein particle formation during retrotransposition. Proc. Natl. Acad. Sci. USA 94:10155-10160.
Kumar, S., K. Tamura, I. B. Jakobsen, and M. Nei. 2001. MEGA2: molecular evolutionary genetics analysis software. Bioinformatics 17:1244-1245.%t)}@@, 百拇医药
Lander, E. S., L. M. Linton, B. Birren, et al 2001. Initial sequencing and analysis of the human genome. Nature 409:860-921.%t)}@@, 百拇医药
Lovsin, N., F. Gubensek, and D. Kordi. 2001. Evolutionary dynamics in a novel L2 clade of non-LTR retrotransposons in Deuterostomia. Mol. Biol. Evol 18:2213-2224.%t)}@@, 百拇医药
Lyngso, C., G. Bouteiller, C. K. Damgaard, D. Ryom, S. Sanchez-Munoz, P. L. Norby, B. J. Bonven, and P. Jorgensen. 2000. Interaction between the transcription factor SPBP and the positive cofactor RNF4. An interplay between protein binding zinc fingers. J. Biol. Chem 275:26144-26149.%t)}@@, 百拇医药
Malik, H. S., W. D. Burke, and T. H. Eickbush. 1999. The age and evolution of non-LTR retrotransposable elements. Mol. Biol. Evol 16:793-805.%t)}@@, 百拇医药
Malik, H. S., and T. H. Eickbush. 1999. Modular evolution of the integrase domain in the Ty3/Gypsy class of LTR retrotransposons. J. Virol 73:5186-5190.
Martin, S. L., and F. D. Bushman. 2001. Nucleic acid chaperone activity of the ORF1 protein from the mouse LINE-1 retrotransposon. Mol. Cell. Biol 21:467-475.3j\e/i#, 百拇医药
Molgaard, A., S. Kauppinen, and S. Larsen. 2000. Rhamnogalacturonan acetylesterase elucidates the structureand function of a new family of hydrolases. Structure Fold Des 8:373-383.3j\e/i#, 百拇医药
Mukabayire, O., and N. J. Besansky. 1996. Distribution of T1, Q, Pegasus and mariner transposable elements on the polytene chromosomes of PEST, a standard strain of Anopheles gambiae. Chromosoma 104:585-595.3j\e/i#, 百拇医药
Murzin, A. G., S. E. Brenner, T. Hubbard, and C. Chothia. 1995. SCOP: a structural classification of proteins database for the investigation of sequences and structures. J. Mol. Biol 247:536-540.3j\e/i#, 百拇医药
Okada, N., M. Hamada, I. Ogiwara, and K. Ohshima. 1997. SINEs and LINEs share common 3' sequences: a review. Gene 205:229-243.3j\e/i#, 百拇医药
Poulter, R., M. Butler, and J. Ormandy. 1999. A LINE element from the pufferfish (fugu) Fugu rubripes which shows similarity to the CR1 family of non-LTR retrotransposons. Gene 227:169-179.
Rosenthal, P. B., X. Zhang, F. Formanowski, W. Fitz, C. H. Wong, H. Meier-Ewert, J. J. Skehel, and D. C. Wiley. 1998. Structure of the haemagglutinin-esterase-fusion glycoprotein of influenza C virus. Nature 396:92-96.s/-c?, http://www.100md.com
Saha, V., T. Chaplin, A. Gregorini, P. Ayton, and B. D. Young. 1995. The leukemia-associated-protein (LAP) domain, a cysteine-rich motif, is present in a wide range of proteins, including MLL, AF10, and MLLT6 proteins. Proc. Natl. Acad. Sci. USA 92:9737-9741.s/-c?, http://www.100md.com
Smit, A. F. 1996. The origin of interspersed repeats in the human genome. Curr. Opin. Genet. Dev 6:743-748.s/-c?, http://www.100md.com
Smit, A. F. 2000. L3, humrep. Repbase Update .s/-c?, http://www.100md.com
Smit, A. F. 2001. REX1_FURC. Repbase Update (fugrep.ref) .s/-c?, http://www.100md.com
Thompson, J. D., D. G. Higgins, and T. J. Gibson. 1994. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22:4673-4680.s/-c?, http://www.100md.com
Vandergon, T. L., and M. Reitman. 1994. Evolution of chicken repeat 1 (CR1) elements: evidence for ancient subfamilies and multiple progenitors. Mol. Biol. Evol 11:886-898.s/-c?, http://www.100md.com
Volff, J. N., C. Korting, and M. Schartl. 2000. Multiple lineages of the non-LTR retrotransposon Rex1 with varying success in invading fish genomes. Mol. Biol. Evol 17:1673-1684.s/-c?, http://www.100md.com
Weiner, A. M. 2000. Do all SINEs lead to LINEs?. Nat. Genet 24:332-333.s/-c?, http://www.100md.com
Wurzer, W. J., K. Obojes, and R. Vlasak. 2002. The sialate-4-O-acetylesterases of coronaviruses related to mouse hepatitis virus: a proposal to reorganize group 2 Coronaviridae. J. Gen. Virol 83:395-402.s/-c?, http://www.100md.com
Yochum, G. S., and D. E. Ayer. 2001. Pf1, a novel PHD zinc finger protein that links the TLE corepressor to the mSin3A-histone deacetylase complex. Mol. Cell. Biol 21:4110-4118.s/-c?, http://www.100md.com
Accepted for publication September 2, 2002.(Vladimir V. Kapitonov and Jerzy Jurka)