Retrotransposable elements R1 and R2 interrupt the rRNA genes of most insects

Divergent non-LTR retrotransposon lineages from the genomes of scorpions (Arachnida: Scorpiones)

We screened across the taxonomic diversity of order Scorpiones (22 species belonging to 21 genera and 10 families) for the presence of seven different clades of non-LTR retrotransposons in their genomes using PCR with newly designed clade-specific consensus-degenerate hybrid oligonucleotide primers. Scorpion genomes were found to contain four known non-LTR retrotransposon clades: R1, I, Jockey, and CR1. In total, 35 fragments of reverse transcriptase genes of new elements from 22 scorpion species were obtained and analyzed for three clades, Jockey, I, and CR1. Phylogenies of different clades of elements were built using amino acid sequences inferred from 33 non-LTR retrotransposon clones. Distinct evolutionary lineages, with several major groups of the non-LTR retroelements were identified, showing significant variation. Four lineages were revealed in Jockey clade. The phylogeny of I clade showed strong support for the monophyletic origin of such group of elements in scorpions. Three separate lineages can be distinguished in the phylogenetic tree of CR1 clade. The large fraction of the isolated elements appeared to be defective.

Transposable element orientation bias in the Drosophila melanogaster genome

Nonrandom distributions of transposable elements can be generated by a variety of genomic features. Using the full D. melanogaster genome as a model, we characterize the orientations of different classes of transposable elements in relation to the directionality of genes. DNA-mediated transposable elements are more likely to be in the same orientation as neighboring genes when they occur in the nontranscribed region's that flank genes. However, RNA-mediated transposable elements located in an intron are more often oriented in the direction opposite to that of the host gene. These orientation biases are strongest for genes with highly biased codon usage, probably reflecting the ability of such loci to respond to weak positive or negative selection. The leading hypothesis for selection against transposable elements in the coding orientation proposes that transcription termination poly(A) signal motifs within retroelements interfere with normal gene transcription. However, after accounting for differences in base composition between the strands, we find no evidence for global selection against spurious transcription termination signals in introns. We therefore conclude that premature termination of host gene transcription due to the presence of poly(A) signal motifs in retroelements might only partially explain strand-specific detrimental effects in the D. melanogaster genome.

Competition between R1 and R2 transposable elements in the 28S rRNA genes of insects

R1 and R2 are non-LTR retrotransposons that insert in the 28S rRNA genes of arthropods. R1 elements insert into a site that is 74 bp downstream of the R2 insertion site, thus the presence of an R2 in the same 28S gene may inhibit the expression of R1. Consistent with such a suggestion, the R1 elements of Drosophila melanogaster have a strong bias against inserting into 28S genes already containing an R2 element. R2 elements, on the other hand, are only 2-3 fold inhibited from inserting into a 28S gene already containing an R1. D. melanogaster R1 elements are unusual in that they generate a 23-bp deletion of the target site upstream of the insertion. Using in vitro assays developed to study R2 integration, we show that the presence of R1 sequences 51 bp downstream of the R2 insertion site changes the nucleosomal structure that can be formed by the R2 target site. The R2 endonuclease is inhibited from cleaving these altered nucleosomes. We suggest that R1 elements have been selected to make this large deletion of the 28S gene to block the insertion of an upstream R2 element. These findings are consistent with the model that R1 and R2 are in competition for the limited number of insertion sites available within their host's genome.

APE-type non-LTR retrotransposons: determinants involved in target site recognition

Non-long terminal repeat (Non-LTR) retrotransposons represent a diverse and widely distributed group of transposable elements and an almost ubiquitous component of eukaryotic genomes that has a major impact on evolution. Their copy number can range from a few to several million and they often make up a significant fraction of the genomes. The members of the dominating subtype of non-LTR retrotransposons code for an endonuclease with homology to apurinic/apyrimidinic endonucleases (APE), and are thus termed APE-type non-LTR retrotransposons. In the last decade both the number of identified non-LTR retrotransposons and our knowledge of biology and evolution of APE-type non-LTR retrotransposons has increased tremendously.

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.

Ticks have R2 retrotransposons but not the consensus transposon target site of other arthropods

Some copies of the large subunit rRNA genes (LSU rDNA) of most arthropods studied to date are inactivated by R-element retrotransposons at a specific target region that is highly conserved in sequence across all kingdoms of organisms. Here we report finding R2 elements in low copy numbers in the LSU rDNA of hard and soft ticks. Although the elements were inserted at the same LSU rDNA location as in insects, there were substitutions in the consensus R2 endonuclease cleavage site in the ticks and some other parasitiform mites. The substituted region comprises a critical contact point with small subunit rRNA, but in vitro structure probing analysis revealed novel, presumably stabilizing base-pairing.

Large-scale analysis of adeno-associated virus vector integration sites in normal human cells

The integration sites of viral vectors used in human gene therapy can have important consequences for safety and efficacy. However, an extensive evaluation of adeno-associated virus (AAV) vector integration sites has not been completed, despite the ongoing use of AAV vectors in clinical trials. Here we have used a shuttle vector system to isolate and analyze 977 unique AAV vector-chromosome integration junctions from normal human fibroblasts and describe their genomic distribution. We found a significant preference for integrating within CpG islands and the first 1 kb of genes, but only a slight overall preference for transcribed sequences. Integration sites were clustered throughout the genome, including a major preference for integration in ribosomal DNA repeats, and 13 other hotspots that contained three or more proviruses within a 500-kb window. Both junctions were localized from 323 proviruses, allowing us to characterize the chromosomal deletions, insertions, and translocations associated with vector integration. These studies establish a profile of insertional mutagenesis for AAV vectors and provide unique insight into the chromosomal distribution of DNA strand breaks that may facilitate integration.

Long-term inheritance of the 28S rDNA-specific retrotransposon R2

R2 is a non-long-terminal-repeat (LTR) retrotransposon that inserts specifically into 28S rDNA. R2 has been identified in many species of arthropods and three species of chordates. R2 may be even more widely distributed in animals, and its origin may be traceable to early animal evolution. In this study, we identified R2 elements in medaka fish, White Cloud Mountain minnow, Reeves' turtle, hagfish, sea lilies, and some arthropod species, using degenerate polymerase chain reaction methods. We also identified two R2 elements from the public genomic sequence database of the bloodfluke Schistosoma mansoni. One of the two bloodfluke R2 elements has two zinc-finger motifs at the N-terminus; this differs from other known R2 elements, which have one or three zinc-finger motifs. Phylogenetic analysis revealed that the whole phylogeny of R2 can be divided into 11 parts (subclades), in which the local R2 phylogeny and the corresponding host phylogeny are consistent. Divergence-versus-age analysis revealed that there is no reliable evidence for the horizontal transfer of R2 but supports the proposition that R2 has been vertically transferred since before the divergence of the deuterostomes and protostomes. The seeming inconsistency between the R2 phylogeny and the phylogeny of their hosts is due to the existence of paralogous lineages. The number of N-terminal zinc-finger motifs is consistent with the deep phylogeny of R2 and indicates that the common ancestor of R2 had three zinc-finger motifs at the N-terminus. This study revealed the long-term vertical inheritance and the ancient origin of sequence specificity of R2, both of which seem applicable to some other non-LTR retrotransposons.

A site-specific insertion sequence in flax genotrophs induced by environment

A single-copy 5.7 kilobase (kb) DNA fragment, termed Linum Insertion Sequence 1 (LIS-1), has been identified and characterized. This is one of the DNA changes associated with the environmentally induced heritable changes resulting in stable lines termed genotrophs in flax (Linum usitatissimum). The insertion sequence and its insertion site have been cloned from genomic libraries and sequenced. PCR products across the insertion and surrounding regions have also been cloned and sequenced. The 5.7 kb DNA fragment is inserted into a 3.7 kb EcoRI fragment in the plastic line (Pl) with the generation of a 3 base pair duplication at the insertion site, as well as additional sequence changes. The identical insertion was also found in other genotrophs and flax varieties. The intact element was not present in Pl but appeared to be generated by a reproducible series of complex rearrangements or insertion events. LIS-1 is the result of a targeted, highly specific, complex insertion event that occurs during the formation of some of the genotrophs, and occurs naturally in many flax and linseed varieties.

Functional roles of 3'-terminal structures of template RNA during in vivo retrotransposition of non-LTR retrotransposon, R1Bm

R1Bm is a non-LTR retrotransposon found specifically within 28S rRNA genes of the silkworm. Different from other non-LTR retrotransposons encoding two open reading frames (ORFs), R1Bm structurally lacks a poly (A) tract at its 3' end. To study how R1Bm initiates reverse transcription from the poly (A)-less template RNA, we established an in vivo retrotransposition system using recombinant baculovirus, and characterized retrotransposition activities of R1Bm. Target-primed reverse transcription (TPRT) of R1Bm occurred from the cleavage site generated by endonuclease (EN). The 147 bp of 3'-untranslated region (3'UTR) was essential for efficient retrotransposition of R1Bm. Even using the complete R1Bm element, however, reverse transcription started from various sites of the template RNA mostly with 5'-UG-3' or 5'-UGU-3' at their 3' ends, which are presumably base-paired with 3' end of the EN-digested 28S rDNA target sequence, 5'-AGTAGATAGGGACA-3'. When the downstream sequence of 28S rDNA target was added to the 3' end of R1 unit, reverse transcription started exactly from the 3' end of 3'UTR and retrotransposition efficiency increased. These results indicate that 3'-terminal structure of template RNA including read-through region interacts with its target rDNA sequences of R1Bm, which plays important roles in initial process of TPRT in vivo.

Characterization of active R2 retrotransposition in the rDNA locus of Drosophila simulans

The rRNA gene (rDNA) loci of all arthropod lineages contain non-LTR retrotransposable elements that have evolved to specifically insert into the 28S rRNA genes. Extensive in vitro experiments have been conducted to investigate the mechanism of R2 retrotransposition but little is known of the insertion frequency or cellular factors that might regulate R2 activity. In this article, isofemale lines obtained from a population of Drosophila simulans were surveyed for recent R2 insertions. Within most lines, all individuals showed the same collection of R2 insertions, providing no evidence for recent R2 activity. However, in a few of the isofemale lines, virtually all individuals differed in their R2 insertion profiles. The descendants of individual pairs of flies from these "active lines" rapidly accumulated new insertions. The frequent insertion of new R2 elements was associated with the elimination of old R2 elements from the rDNA locus. The existence of lines in which R2 retrotransposes frequently and lines in which the elements appear dormant suggests that cellular mechanisms that can regulate the activity of R2 exist. Retrotransposition activity was correlated with the number of full-length R2 elements but not with the size of the rDNA locus or the number of uninserted units.

Evolution of the Transposable Element Pokey in the Ribosomal DNA of Species in the Subgenus Daphnia (Crustacea: Cladocera)

Pokey is a member of the piggyBac (previously called the TTAA-specific) family of transposons and inserts into a conserved region of the large subunit ribosomal RNA gene. This location is a "hot spot" for insertional activity, as it is known to contain other arthropod transposable elements. However, Pokey is unique in that it is the first DNA transposon yet known to insert into this region. All other insertions are class I non-LTR retrotransposons. This study surveyed variation in Pokey elements through phylogenetic analysis of the 3' ends of Pokey elements from ribosomal DNA (rDNA) in species from the nominate subgenus of the genus Daphnia (Crustacea: Cladocera). The results suggest that Pokey has been stably, vertically inherited within rDNA over long periods of evolutionary time. No evidence was found to support horizontal transfer, which commonly occurs in other DNA transposons, such as P and mariner. Furthermore, Pokey has diverged into sublineages that have persisted across speciation events in some groups. In addition, a new highly divergent paralogous Pokey element was discovered in the rDNA of one species.

Crystal structure of the endonuclease domain encoded by the telomere-specific long interspersed nuclear element, TRAS1

The telomere-specific long interspersed nuclear element, TRAS1, encodes an endonuclease domain, TRAS1-EN, which specifically cleaves the telomeric repeat targets (TTAGG)n of insects and (TTAGGG)n of vertebrates. To elucidate the sequence-specific recognition properties of TRAS1-EN, we determined the crystal structure at 2.4-A resolution. TRAS1-EN has a four-layered alpha/beta sandwich structure; its topology is similar to apurinic/apyrimidinic endonucleases, but the beta-hairpin (beta10-beta11) at the edge of the DNA-binding surface makes an extra loop that distinguishes TRAS1-EN from cellular apurinic/apyrimidinic endonucleases. A protein-DNA complex model suggests that the beta10-beta11 hairpin fits into the minor groove, enabling interaction with the telomeric repeats. Mutational studies of TRAS1-EN also indicated that the Asp-130 and beta10-beta11 hairpin structure are involved in specific recognition of telomeric repeats.

R1 and R2 retrotransposition and deletion in the rDNA loci on the X and Y chromosomes of Drosophila melanogaster

The non-LTR retrotransposons R1 and R2 insert into the 28S rRNA genes of arthropods. Comparisons among Drosophila lineages have shown that these elements are vertically inherited, while studies within species have indicated a rapid turnover of individual copies (elimination of old copies and the insertion of new copies). To better understand the turnover of R1 and R2, 200 retrotranspositions and nearly 100 eliminations have been scored in the Harwich mutation-accumulation lines of Drosophila melanogaster. Because the rDNA arrays in D. melanogaster are present on the X and Y chromosomes and no exchanges were detected in these lines, it was possible to show that R1 retrotranspositions occur predominantly in the male germ line, while R2 retrotranspositions were more evenly divided between the germ lines of both sexes. The rate of elimination of elements from the Y rDNA array was twice that of the X rDNA array with both chromosomal loci containing regions where the rate of elimination was on average eight times higher. Most R1 and R2 eliminations appear to occur by large intrachromosomal events (i.e., loop-out events) that involve multiple rDNA units. These findings are interpreted in light of the known abundance of R1 and R2 elements in the X and Y rDNA loci of D. melanogaster.

The structure of a single unit of ribosomal RNA gene (rDNA) including intergenic subrepeats in the Australian bulldog ant Myrmecia croslandi (Hymenoptera: Formicidae)

A complete single unit of a ribosomal RNA gene (rDNA) of M. croslandi was sequenced. The ends of the 18S, 5.8S and 28S rRNA genes were determined by using the sequences of D. melanogaster rDNAs as references. Each of the tandemly repeated rDNA units consists of coding and non-coding regions whose arrangement is the same as that of D. melanogaster rDNA. The intergenic spacer (IGS) contains, as in other species, a region with subrepeats, of which the sequences are different from those previously reported in other insect species. The length of IGSs was estimated to be 7-12 kb by genomic Southern hybridization, showing that an rDNA repeating unit of M. croslandi is 14-19 kb-long. The sequences of the coding regions are highly conserved, whereas IGS and ITS (internal transcribed spacer) sequences are not. We obtained clones with insertions of various sizes of R2 elements, the target sequence of which was found in the 28S rRNA coding region. A short segment in the IGS that follows the 3' end of the 28S rRNA gene was predicted to form a secondary structure with long stems.

Integration of the 5' end of the retrotransposon, R2Bm, can be complemented by homologous recombination

R2Bm is a non-long-terminal-repeat (non-LTR) retrotransposon that was identified at a specific target site in the 28S rRNA genes of the silkworm, Bombyx mori. Although in vitro analysis has revealed that the 3' end of R2Bm is integrated into the target site by means of target-primed reverse transcription (TPRT), the mechanism of the 5' end integration is not well understood. We established a novel in vivo system to assay the insertion mechanism of R2Bm using a cultured cell line, C65, and a baculovirus, AcNPV, as host and vector, respectively. The 3' end of R2Bm integrated at the target site in the rRNA genes of C65 cells when an AcNPV containing both the full-length 3' UTR and the entire open reading frame (ORF) of R2Bm was introduced while the 5' end integration was incorrect. The 5' end of R2Bm was integrated, however, when the 28S gene sequence upstream of the R2Bm target site was added to the R2Bm sequence. Thus, in our assay, homologous sequences were likely essential for the successful integration of the entire R2Bm into the host cell genome. We also demonstrated that the failure to integrate caused by a frame-shifted ORF was rescued by co-infection with a helper virus that contained only the R2Bm ORF. This indicates that R2 retrotransposition can be complemented in trans. These findings suggest that the host's mechanism for DNA repair may be necessary for the integration of the 5' end of R2Bm and that R2Bm protein may only have the ability to integrate the 3' end of the element by TPRT.

Effects of recombination rate and gene density on transposable element distributions in Arabidopsis thaliana

Transposable elements (TEs) comprise a major component of eukaryotic genomes, and exhibit striking deviations from random distribution across the genomes studied, including humans, flies, nematodes, and plants. Although considerable progress has been made in documenting these patterns, the causes are subject to debate. Here, we use the genome sequence of Arabidopsis thaliana to test for the importance of competing models of natural selection against TE insertions. We show that, despite TE accumulation near the centromeres, recombination does not generally correlate with TE abundance, suggesting that selection against ectopic recombination does not influence TE distribution in A. thaliana. In contrast, a consistent negative correlation between gene density and TE abundance, and a strong under-representation of TE insertions in introns suggest that selection against TE disruption of gene expression is playing a more important role in A. thaliana. High rates of self-fertilization may reduce the importance of recombination rate in genome structuring in inbreeding organisms such as A. thaliana and Caenorhabditis elegans.

R2 retrotransposition on assembled nucleosomes depends on the translational position of the target site

R2 retrotransposons insert into the 28S rRNA genes of insects. Integration occurs by specific cleavage of the target site and utilization of the released DNA end to prime reverse transcription of the RNA transcript. Specificity of the protein to the target site is dependent upon nucleotide sequence recognition extending from 35 bp upstream to 15 bp downstream of the cleavage site. In this report, we show that sequence recognition and cleavage by the R2 protein can occur while the target site is assembled into nucleosomes. Reconstitution of DNA fragments containing the 28S gene sequence into a set of nucleosomes with different translational frames revealed that the R2 site adopted the same rotational orientation with respect to the histone octamer. Binding and cleavage by the R2 protein were most efficient when the upstream binding site for the R2 protein was near a nucleosome end. Interaction of the R2 protein with the nucleosome disrupted the histone:DNA contacts in the 50 bp region directly bound by R2, but did not modify the remainder of the nucleosome structure.

Ribosomal DNA variation in the grasshopper, Dichroplus elongatus

We report an RFLP analysis of ribosomal DNA variation in natural populations of the grasshopper, Dichroplus elongatus, previously analyzed for mitochondrial DNA variation. DNA samples were digested with five restriction enzymes, BamHI, EcoRI, HindIII, PstI, and XbaI. BamHI was the only enzyme that showed no variation. The remaining enzymes showed fragment size variation at both intra- and interpopulation levels. Stepwise regression analysis revealed that the average number of length variants per individual is significantly associated with altitude. Moreover, the same analysis indicated that the frequency of some restriction variants exhibits a significant regression on both geographic and climatic variables. The intra- and interpopulation variability of rDNA was analysed by Lynch's and Hedrick's similarity indices using presence or absence of a fragment and band intensities, respectively. The corresponding neighbour-joining (N-J) trees based on Lynch's and Hedrick's genetic distances resulted in similar topologies. However, these trees were not in agreement with the N-J dendrogram obtained from mtDNA data previously reported by Clemente et al. (2000). The disagreement between mtDNA and rDNA phenograms along with the observed correlation between rDNA variability and geographical and climatic variables suggest some form of selection, besides genetic drift and migration, is involved in the pattern of rDNA variation.

Rates of R1 and R2 retrotransposition and elimination from the rDNA locus of Drosophila melanogaster

R1 and R2 elements are non-LTR retrotransposons that insert specifically into the 28S rRNA genes of arthropods. The process of concerted evolution of the rDNA locus should give rise to rapid turnover of these mobile elements compared to elements that insert at sites throughout a genome. To estimate the rate of R1 and R2 turnover we have examined the insertion of new elements and elimination of old elements in the Harwich mutation accumulation lines of Drosophila melanogaster, a set of inbred lines maintained for >350 generations. Nearly 300 new insertion and elimination events were observed in the 19 Harwich lines. The retrotransposition rate for R1 was 18 times higher than the retrotransposition rate for R2. Both rates were within the range previously found for retrotransposons that insert outside the rDNA loci in D. melanogaster. The elimination rates of R1 and R2 from the rDNA locus were similar to each other but over two orders of magnitude higher than that found for other retrotransposons. The high rates of R1 and R2 elimination from the rDNA locus confirm that these elements must maintain relatively high rates of retrotransposition to ensure their continued presence in this locus.

Dynalign: an algorithm for finding the secondary structure common to two RNA sequences

With the rapid increase in the size of the genome sequence database, computational analysis of RNA will become increasingly important in revealing structure-function relationships and potential drug targets. RNA secondary structure prediction for a single sequence is 73 % accurate on average for a large database of known secondary structures. This level of accuracy provides a good starting point for determining a secondary structure either by comparative sequence analysis or by the interpretation of experimental studies. Dynalign is a new computer algorithm that improves the accuracy of structure prediction by combining free energy minimization and comparative sequence analysis to find a low free energy structure common to two sequences without requiring any sequence identity. It uses a dynamic programming construct suggested by Sankoff. Dynalign, however, restricts the maximum distance, M, allowed between aligned nucleotides in the two sequences. This makes the calculation tractable because the complexity is simplified to O(M(3)N(3)), where N is the length of the shorter sequence. The accuracy of Dynalign was tested with sets of 13 tRNAs, seven 5 S rRNAs, and two R2 3' UTR sequences. On average, Dynalign predicted 86.1 % of known base-pairs in the tRNAs, as compared to 59.7 % for free energy minimization alone. For the 5 S rRNAs, the average accuracy improves from 47.8 % to 86.4 %. The secondary structure of the R2 3' UTR from Drosophila takahashii is poorly predicted by standard free energy minimization. With Dynalign, however, the structure predicted in tandem with the sequence from Drosophila melanogaster nearly matches the structure determined by comparative sequence analysis.

The transposable elements of the Drosophila melanogaster euchromatin: a genomics perspective

BACKGROUND

Transposable elements are found in the genomes of nearly all eukaryotes. The recent completion of the Release 3 euchromatic genomic sequence of Drosophila melanogaster by the Berkeley Drosophila Genome Project has provided precise sequence for the repetitive elements in the Drosophila euchromatin. We have used this genomic sequence to describe the euchromatic transposable elements in the sequenced strain of this species.

RESULTS

We identified 85 known and eight novel families of transposable element varying in copy number from one to 146. A total of 1,572 full and partial transposable elements were identified, comprising 3.86% of the sequence. More than two-thirds of the transposable elements are partial. The density of transposable elements increases an average of 4.7 times in the centromere-proximal regions of each of the major chromosome arms. We found that transposable elements are preferentially found outside genes; only 436 of 1,572 transposable elements are contained within the 61.4 Mb of sequence that is annotated as being transcribed. A large proportion of transposable elements is found nested within other elements of the same or different classes. Lastly, an analysis of structural variation from different families reveals distinct patterns of deletion for elements belonging to different classes.

CONCLUSIONS

This analysis represents an initial characterization of the transposable elements in the Release 3 euchromatic genomic sequence of D. melanogaster for which comparison to the transposable elements of other organisms can begin to be made. These data have been made available on the Berkeley Drosophila Genome Project website for future analyses.

Dynamics of R1 and R2 elements in the rDNA locus of Drosophila simulans

The mobile elements R1 and R2 insert specifically into the rRNA gene locus (rDNA locus) of arthropods, a locus known to undergo concerted evolution, the recombinational processes that preserve the sequence homogeneity of all repeats. To monitor how rapidly individual R1 and R2 insertions are turned over in the rDNA locus by these processes, we have taken advantage of the many 5' truncation variants that are generated during the target-primed reverse transcription mechanism used by these non-LTR retrotransposons for their integration. A simple PCR assay was designed to reveal the pattern of the 5' variants present in the rDNA loci of individual X chromosomes in a population of Drosophila simulans. Each rDNA locus in this population was found to have a large, unique collection of 5' variants. Each variant was present at low copy number, usually one copy per chromosome, and was seldom distributed to other chromosomes in the population. The failure of these variants to spread to other units in the same rDNA locus suggests a strong recombinational bias against R1 and R2 that results in the individual copies of these elements being rapidly lost from the rDNA locus. This bias suggests a significantly higher frequency of R1 and R2 retrotransposition than we have previously suggested.

Characterization of a retrotransposon-like element from Entamoeba histolytica

The protozoan parasite Entamoeba histolytica is the causative agent of amoebiasis. The genome organization of this organism is not well understood. We had earlier reported the presence of a multicopy sequence, HMc, in E. histolytica. Subsequent analysis showed that HMc is a member of a retrotransposon family that we have named the E. histolytica retrotransposon-like element (EhRLE). Four other members of this family have been characterized. The EhRLE family is distributed across all chromosomes of the parasite. There are 140 copies, which show minor sequence variation with respect to one another (2--4% from the consensus sequence). From a sequence analysis of five members of the EhRLE family, the complete EhRLE unit is estimated to be 4086 bp in length. It has a 27-mer inverted repeat at its ends. A pairwise comparison with sequences in the database showed a highly significant match of a part of EhRLE with reverse transcriptases (RT), especially those encoded by non-long terminal repeat retrotransposons. There are stop codons in all the five EhRLEs, but a continuous open reading frame of 464 amino acids could be reconstructed by comparing the sequences of several EhRLEs. The reconstructed sequence showed a much better identity with RT as compared with any of the original EhRLE sequences. The non-pathogenic species, Entamoeba dispar, also contains this element, with 85% sequence identity with EhRLE. The data suggest that EhRLE may be a retrotransposon, but many of its members are probably nonfunctional due to the accumulation of mutations.

Structural and phylogenetic analysis of TRAS, telomeric repeat-specific non-LTR retrotransposon families in Lepidopteran insects

TRAS1 is a non-LTR retrotransposon inserted specifically into the telomeric repeat (TTAGG)(n) in the silkworm, Bombyx mori. To characterize the evolutionary origin of TRAS-like elements, we identified seven TRAS families (TRAS3, TRAS4, TRAS5, TRAS6, TRASY, TRASZ, and TRASW) from B. mori and four elements from two Lepidoptera, Dictyoploca japonica (TRASDJ) and Samia cynthia ricini (TRASSC3, TRASSC4, and TRASSC9). More than 2,000 copies of various Bombyx TRAS elements accumulated within (TTAGG)(n) sequences as unusual but orderly tandem repeats. The 5' and 3' regions were highly conserved within each class of Bombyx TRAS elements without truncation. This suggests that distinct classes of TRAS have been maintained independently by retrotransposition into (TTAGG)(n). The phylogenetic tree of site-specific retroelements showed that nine TRAS families in Lepidoptera constitute a single phylogenetic group that is closely related to the R1 family that inserts specifically into arthropod 28S rDNA. The higher amino acid sequence identity from endonuclease (EN) to reverse transcriptase (RT) domains between TRAS groups (about 37%-70%) than among TRAS elements and R1Bm (about 25%-30%), may reflect the presence of some DNA structure responsible for their target specificity. Sequence comparison from EN to RT domains among non-LTR elements revealed several regions conserved only within TRAS elements. We found a highly conserved region that resembles the Myb-like DNA-binding structure, between the EN and RT domains. These regions may be involved in site-specific integration of TRAS elements into the (TTAGG)(n) telomeric repeats.

Multiple lineages of R1 retrotransposable elements can coexist in the rDNA loci of Drosophila

R1 non-long terminal repeat retrotransposable elements insert specifically into the 28S rRNA genes of arthropods. One aspect of R1 evolution that has been difficult to explain is the presence of divergent lineages of R1 in the rDNA loci of the same species. Multiple lineages should compete for a limited number of insertion sites, in addition to being subject to the concerted evolution processes homogenizing the rRNA genes. The presence of multiple lineages suggests either the ability of the elements to overcome these factors and diverge within rDNA loci, or the introduction of new lineages by horizontal transmission. To address this issue, we attempted to characterize the complete set of R1 elements in the rDNA locus from five Drosophila species groups (melanogaster, obscura, testacea, quinaria, and repleta). Two major R1 lineages, A and B, that diverged about 100 MYA were found to exist in Drosophila. Elements of the A lineage were found in all 35 Drosophila species tested, while elements of the B lineage were found in only 11 species from three species groups. Phylogenetic analysis of the R1 elements, supported by comparison of their rates of nucleotide sequence substitution, revealed that both the A and the B lineages have been maintained by vertical descent. The B lineage was less stable and has undergone numerous, independent elimination events, while the A lineage has diverged into three sublineages, which were, in turn, differentially stable. We conclude that while the differential retention of multiple lineages greatly complicates its phylogenetic history, the available R1 data continue to be consistent with the strict vertical descent of these elements.

Identification of Waldo-A and Waldo-B, two closely related non-LTR retrotransposons in Drosophila

We have identified two novel, closely related subfamilies of non-long-terminal-repeat (non-LTR) retrotransposons in Drosophila melanogaster, the Waldo-A and Waldo-B subfamilies, that are in the same lineage as site-specific LTR retrotransposons of the R1 clade. Both contain potentially active copies with two large open reading frames, having coding capacities for a nucleoprotein as well as endonuclease and reverse transcriptase activities. Many copies are truncated at the 5' end, and most are surrounded by target site duplications of variable lengths. Elements of both subfamilies have a nonrandom distribution in the genome, often being inserted within or very close to (CA)(n) arrays. At the DNA level, the longest elements of Waldo-A and Waldo-B are 69% identical on their entire length, except for the 5' untranslated regions, which have a mosaic organization, suggesting that one arose from the other following new promoter acquisition. This event occurred before the speciation of the D. melanogaster subgroup of species, since both Waldo-A and Waldo-B coexist in other species of this subgroup.

Sequence-specific recognition and cleavage of telomeric repeat (TTAGG)(n) by endonuclease of non-long terminal repeat retrotransposon TRAS1

The telomere of the silkworm Bombyx mori consists of (TTAGG/CCTAA)(n) repeats and harbors a large number of telomeric repeat-specific non-long terminal repeat retrotransposons, such as TRAS1 and SART1. To understand how these retrotransposons recognize and integrate into the telomeric repeat in a sequence-specific manner, we expressed the apurinic-apryrimidinic endonuclease-like endonuclease domain of TRAS1 (TRAS1 EN), which is supposed to digest the target DNA, and characterized its enzymatic properties. Purified TRAS1 EN could generate specific nicks on both strands of the telomeric repeat sequence between T and A of the (TTAGG)(n) strand (bottom strand) and between C and T of the (CCTAA)(n) strand (top strand). These sites are consistent with insertion sites expected from the genomic structure of boundary regions of TRAS1. Time course studies of nicking activities on both strands revealed that the cleavages on the bottom strand preceded those on the top strand, supporting the target-primed reverse transcription model. TRAS1 EN could cleave the telomeric repeats specifically even if it was flanked by longer tracts of nontelomeric sequence, indicating that the target site specificity of the TRAS1 element was mainly determined by its EN domain. Based on mutation analyses, TRAS1 EN recognizes less than 10 bp around the initial cleavage site (upstream 7 bp and downstream 3 bp), and the GTTAG sequence especially is essential for the cleavage reaction on the bottom strand (5'. TTAGGTT downward arrow AGG. 3'). TRAS1 EN, the first identified endonuclease digesting telomeric repeats, may be used as a genetic tool to shorten the telomere in insects and some other organisms.

Retrotransposition of the I factor, a non-long terminal repeat retrotransposon of Drosophila, generates tandem repeats at the 3' end

Non-long terminal repeat (LTR) retrotransposons or LINEs transpose by reverse transcription of an RNA intermediate and are thought to use the 3' hydroxyl of a chromosomal cleavage to initiate synthesis of the first strand of the cDNA. Many of them terminate in a poly(dA) sequence at the 3' end of the coding strand although some, like the I factor of Drosophila melanogaster, have 3' ends formed by repeats of the trinucleotide TAA. We report results showing that I factor transcripts end a few nucleotides downstream of the TAA repeats and that these extra nucleotides are not integrated into chromosomal DNA during retrotransposition. We also show that the TAA repeats are not required for transposition and that I elements containing mutations affecting the TAA sequences generate transposed copies ending with tandem repeats of various types. Our results suggest that during integration the 3' end of the I factor RNA template can pair with nucleotides at the target site and that tandem duplications are generated by the reverse transcriptase of the I factor in a manner that is reminiscent of the activity of the reverse transcriptases of telomerases. Reverse transcriptases of other non-LTR retrotransposons may function in a similar way.

Target specificity of the endonuclease from the Xenopus laevis non-long terminal repeat retrotransposon, Tx1L

Elements of the Tx1L family are non-long terminal repeat retrotransposons (NLRs) that are dispersed in the genome of Xenopus laevis. Essentially all genomic copies of Tx1L are found inserted at a specific site within another family of transposable elements (Tx1D). This suggests that Tx1L is a site-specific retrotransposon. Like many (but not all) other NLRs, the Xenopus element encodes an apparent endonuclease that is related in sequence to the apurinic-apyrimidinic endonucleases that participate in DNA repair. This enzyme is thought to introduce the single-strand break in target DNA that initiates transposition by the target-primed reverse transcription (TPRT) mechanism. To explore the issue of target specificity more fully, we expressed the polypeptide encoded by the endonuclease domain of open reading frame 2 from Tx1L (Tx1L EN) and characterized its cleavage capabilities. This endonuclease makes a specific nick in the bottom strand precisely at one end of the presumed Tx1L target duplication. Because this activity leaves a 5'-phosphate and 3'-hydroxyl at the nick, it has the location and chemistry required to initiate new insertion events by TPRT. Tx1L EN does not make a specific cut at a preferred target site for Tx1D elements, ruling out the alternative possibility that the composite Tx1L-Tx1D element moves as a unit under the control of functions encoded by Tx1L. Further characterization revealed that the endonuclease remains active for many hours at room temperature and that it is capable of enzymatic turnover. Scanning substitution mutagenesis located the recognition site for Tx1L EN within 10 bp surrounding the primary nick site. Implications of these features for natural transposition events are discussed.

Retrotransposable elements R1 and R2 in the rDNA units of Drosophila mercatorum: abnormal abdomen revisited

R1 and R2 retrotransposable elements are stable components of the 28S rRNA genes of arthropods. While each retrotransposition event leads to incremental losses of rDNA unit expression, little is known about the selective consequences of these elements on the host genome. Previous reports suggested that in the abnormal abdomen (aa) phenotype of Drosophila mercatorum, high levels of rDNA insertions (R1) in conjunction with the under-replication locus (ur), enable the utilization of different ecological conditions via a population level shift to younger age. We have sequenced the R1 and R2 elements of D. mercatorum and show that the levels of R1- and R2-inserted rDNA units were inaccurately scored in the original studies of aa, leading to several misinterpretations. In particular, contrary to earlier reports, aa flies differentially underreplicate R1- and R2-inserted rDNA units, like other species of Drosophila. However, aa flies do not undergo the lower level of underreplication of their functional rDNA units (general underreplication) that is seen in wild-type strains. The lack of general underreplication is expected to confer a selective advantage and, thus, can be interpreted as an adaptation to overcome high levels of R1 and R2 insertions. These results allow us to reconcile some of the apparently contradictory effects of aa and the bobbed phenotype found in other species of Drosophila.

Conserved features at the 5 end of Drosophila R2 retrotransposable elements: implications for transcription and translation

R2 non-LTR retrotransposable elements insert site-specifically into the 28S ribosomal genes of insects. The sequence of the 5' end of full-length R2 elements from thirteen species of Drosophila were compared. Sequences within the 5' untranslated region (5' UTR) revealed little to suggest the presence of a promoter. Protein translation initiates within the 5' UTR and requires the bypassing of a highly conserved termination codon preceding the single R2 open reading frame. This bypassing probably involves a conserved RNA secondary structure which brings a potential initiation codon into close proximity to this termination codon. The most highly conserved sequence within the 5' UTR has properties similar to internal ribosomal entry sites. Based on these findings, we propose that R2elements are co-transcribed with the 28S gene and are translated as part of a large ribosomal subunit.

Site-specific insertion of a SINE-like element, Cp1, into centromeric tandem repeats from Chironomus pallidivittatus

A SINE-like dispersed element, Cp1, from the dipteran Chironomus pallidivittatus was found to show site-specific insertion into two different centromeric tandem repeats. The insertions result in identical target site duplications of nine base-pairs. In contrast, extracentromeric Cp1 elements, which are polymorphic and degenerate, are previously known to be surrounded by different target site duplications. The intracentromeric Cp1 is uniform in structure and contains a single pol III unit, upstream of which 87 bp arms of a palindrome surround a 103 bp unique sequence. The numbers of Cp1 elements per centromere were determined in microdissected material and were found to be in the range of five to ten units per centromere. The well-defined insertion properties, correlated to chromosomal localization, suggest that Cp1 is likely to be a component of importance for the centromere. Similarities of Cp1 and its parts to functionally identified centromeres in Saccharomyces cerevisiae and Schizosaccharomyces pombe are discussed.

Molecular phylogenetics at the Felsenstein zone: approaching the Strepsiptera problem using 5.8S and 28S rDNA sequences

Recent efforts to reconstruct the phylogenetic position of the insect order Strepsiptera have elicited a major controversy in molecular phylogenetics. We sequenced the 5.8S rDNA and major parts of the 28S rDNA 5' region of the strepsipteran species Stylops melittae. Their evolutionary dynamics were analyzed together with previously published insect rDNA sequences to identify tree estimation bias risks and to explore additional sources of phylogenetic information. Several major secondary structure changes were found as being autapomorphic for the Diptera, the Strepsiptera, or the Archaeognatha. Besides elevated substitution rates a significant AT bias was present in dipteran and strepsipteran 28S rDNA which, however, was restricted to stem sites in the Diptera while also affecting single-stranded sites in the Strepsiptera. When dipteran taxa were excluded from tree estimation all methods consistently supported the placement of Strepsiptera to within the Holometabola. When dipteran taxa were included maximum likelihood continued to favor a sister-group relationship of Strepsiptera with Mecoptera while remaining methods strongly supported a sister-group relationship with Diptera. Parametric bootstrap analysis revealed maximum likelihood as a consistent estimator if rate heterogeneity across sites was taken into account. Though the position of Strepsiptera within Holometabola remains elusive, we conclude that the evolution of dipteran and strepsipteran rDNA involved similar yet independent changes of substitution parameters.

RNA-induced changes in the activity of the endonuclease encoded by the R2 retrotransposable element

R2 is a non-long terminal repeat retrotransposable element that inserts itself site specifically in the 28S rRNA genes of arthropods. The 120-kDa protein encoded by R2 has been shown to cleave one strand of the 28S gene at the target site and to use the 3' hydroxyl group generated from this nick to prime reverse transcription of its own RNA. This reaction has been termed target-primed reverse transcription (TPRT). Cleavage of the second DNA strand can occur in the presence or absence of reverse transcription but requires RNA. In this study, more sensitive in vitro assays have enabled further characterization of these reactions. R2 protein is capable of only a single round of TPRT because, once bound to the target DNA, it does not dissociate at physiological ionic strengths. Analysis of the role of RNA in the DNA cleavage reaction has revealed that the binding of RNA induces the R2 protein to form a multimeric complex. While larger complexes may form, the active component appears to be a dimer based on sedimentation studies and the change in stoichiometry of the cleavage reaction from a 1:1 ratio of protein subunit to target DNA in the absence of RNA to a 2:1 ratio of subunit to DNA target in the presence of RNA. Nonspecific RNA can also induce formation of this RNA-protein (RNP) complex, but the association of the protein with R2 RNA is stronger as revealed by its stability in 0.4 M NaCl. Finally, formation of the RNP complex gives rise to a 150-fold increase in the ability of the R2 endonuclease to find the target site. The specificity of this RNP complex is sufficiently great that it can find the 28S gene target site and conduct the TPRT reaction with total genomic DNA.

Transposable elements and genome organization: a comprehensive survey of retrotransposons revealed by the complete Saccharomyces cerevisiae genome sequence

We conducted a genome-wide survey of Saccharomyces cerevisiae retrotransposons and identified a total of 331 insertions, including 217 Ty1, 34 Ty2, 41 Ty3, 32 Ty4, and 7 Ty5 elements. Eighty-five percent of insertions were solo long terminal repeats (LTRs) or LTR fragments. Overall, retrotransposon sequences constitute >377 kb or 3.1% of the genome. Independent evolution of retrotransposon sequences was evidenced by the identification of a single-base pair insertion/deletion that distinguishes the highly similar Ty1 and Ty2 LTRs and the identification of a distinct Ty1 subfamily (Ty1'). Whereas Ty1, Ty2, and Ty5 LTRs displayed a broad range of sequence diversity (typically ranging from 70%-99% identity), Ty3 and Ty4 LTRs were highly similar within each element family (most sharing >96% nucleotide identity). Therefore, Ty3 and Ty4 may be more recent additions to the S. cerevisiae genome and perhaps entered through horizontal transfer or past polyploidization events. Distribution of Ty elements is distinctly nonrandom: 90% of Ty1, 82% of Ty2, 95% of Ty3, and 88% of Ty4 insertions were found within 750 bases of tRNA genes or other genes transcribed by RNA polymerase III. tRNA genes are the principle determinant of retrotransposon distribution, and there is, on average, 1.2 insertions per tRNA gene. Evidence for recombination was found near many Ty elements, particularly those not associated with tRNA gene targets. For these insertions, 5'- and 3'-flanking sequences were often duplicated and rearranged among multiple chromosomes, indicating that recombination between retrotransposons can influence genome organization. S. cerevisiae offers the first opportunity to view organizational and evolutionary trends among retrotransposons at the genome level, and we hope our compiled data will serve as a starting point for further investigation and for comparison to other, more complex genomes.

Retrotransposon R1Bm endonuclease cleaves the target sequence

The R1Bm element, found in the silkworm Bombyx mori, is a member of a group of widely distributed retrotransposons that lack long terminal repeats. Some of these elements are highly sequence-specific and others, like the human L1 sequence, are less so. The majority of R1Bm elements are associated with ribosomal DNA (rDNA). R1Bm inserts into 28S rDNA at a specific sequence; after insertion it is flanked by a specific 14-bp target site duplication of the 28S rDNA. The basis for this sequence specificity is unknown. We show that R1Bm encodes an enzyme related to the endonuclease found in the human L1 retrotransposon and also to the apurinic/apyrimidinic endonucleases. We expressed and purified the enzyme from bacteria and showed that it cleaves in vitro precisely at the positions in rDNA corresponding to the boundaries of the 14-bp target site duplication. We conclude that the function of the retrotransposon endonucleases is to define and cleave target site DNA.

Transposable element-host interactions: regulation of insertion and excision

Transposable elements propagate by inserting into new locations in the genomes of the hosts they inhabit. Their transposition might thus negatively affect the fitness of the host, suggesting the requirement for a tight control in the regulation of transposable element mobilization. The nature of this control depends on the structure of the transposable element. DNA elements encode a transposase that is necessary, and in most cases sufficient, for mobilization. In general, regulation of these elements depends on intrinsic factors with little direct input from the host. Retrotransposons require an RNA intermediate for transposition, and their frequency of mobilization is controlled at multiple steps by the host genome by regulating both their expression levels and their insertional specificity. As a result, a symbiotic relationship has developed between transposable elements and their host. Examples are now emerging showing that transposons can contribute significantly to the well being of the organisms they populate.

A new family of site-specific retrotransposons, SART1, is inserted into telomeric repeats of the silkworm, Bombyx mori

The telomeres of the silkworm, Bombyx mori, consist of pentanucleotide repeats (TTAGG)n . We previously characterized the non-LTR element TRAS1, which terminates with oligo (A) in a head to tail orientation at the exact position (between A and C) of the (CCTAA) n repeats. Here we characterized another family of telomere-specific non-LTR retrotransposon named SART1. The SART1 family was inserted at another site of the (TTAGG) n in a reverse orientation from that of TRAS1. The complete unit of SART1, 6.7 kb in length with a poly (A) stretch, contains two open reading frames encoding putative gag and pol products, overlapping by 54 bp in the -1 reading frame. Most of the 600 SART1 copies in the silkworm haploid genome are completely conserved in structure without 5'truncation. All SART1 sequences analyzed were inserted at the same position (between T and A) within the (TTAGG) n repeats. Fluorescence in situ hybridization showed that many of the SART1 copies were localized in the chromosomal ends. A phylogenetic tree showed that the SART1, TRAS1 and two other site-specific elements, R1 and RT, which insert into 28S ribosomal RNA genes in insects, belong to the same group. Based on the orientation for the chromosomal insertion and structural similarities, these elements could be further classified into two subgroups, R1/TRAS1 and RT/SART1, suggesting that the target specificity of the two telomere-associated elements was changed independently.

An unusual form of transcriptional silencing in yeast ribosomal DNA

Generalized transcriptional repression of large chromosomal regions in Saccharomyces cerevisiae occurs at the silent mating loci and at telomeres and is mediated by the silent information regulator (SIR) genes. We have identified a novel form of transcriptional silencing in S. cerevisiae in the ribosomal DNA (rDNA) tandem array. Ty1 retrotransposons marked with a weakened URA3 gene (Ty1-mURA3) efficiently integrated into rDNA. The mURA3 marker in rDNA was transcriptionally silenced in a SIR2-dependent manner. MET15 and LEU2 were also partially silenced, indicating that rDNA silencing may be quite general. Deletion of SIR4 enhanced mURA3 and MET15 silencing, but deletion of SIR1 or SIR3 did not affect silencing, indicating that the mechanism of silencing differs from that at telomeres and silent mating loci. Deletion of SIR2 resulted in increased psoralen cross-linking of the rDNA in vivo, suggesting that a specific chromatin structure in rDNA down-regulates polymerase II promoters.

Target site selection in transposition

Transposable elements are discrete mobile DNA segments that can insert into non-homologous target sites. Diverse patterns of target site selectivity are observed: Some elements display considerable target site selectivity and others display little obvious selectivity, although none appears to be truly "random." A variety of mechanisms for target site selection are used: Some elements use direct interactions between the recombinase and target DNA whereas other elements depend upon interactions with accessory proteins that communicate both with the target DNA and the recombinase. The study of target site selectivity is useful in probing recombination mechanisms, in studying genome structure and function, and also in providing tools for genome manipulation.

Secondary structure model of the RNA recognized by the reverse transcriptase from the R2 retrotransposable element

RNA transcripts corresponding to the 250-nt 3' untranslated region of the R2 non-LTR retrotransposable element are recognized by the R2 reverse transcriptase and are sufficient to serve as templates in the target DNA-primed reverse transcription (TPRT) reaction. The R2 protein encoded by the Bombyx mori R2 can recognize this region from both the B. mori and Drosophila melanogaster R2 elements even though these regions show little nucleotide sequence identity. A model for the RNA secondary structure of the 3' untranslated region of the D. melanogaster R2 retrotransposon was developed by sequence comparison of 10 species aided by free energy minimization. Chemical modification experiments are consistent with this prediction. A secondary structure model for the 3' untranslated region of R2 RNA from the R2 element from B. mori was obtained by a combination of chemical modification data and free energy minimization. These two secondary structure models, found independently, share several common sites. This study shows the utility of combining free energy minimization, sequence comparison, and chemical modification to model an RNA secondary structure.

Downstream 28S gene sequences on the RNA template affect the choice of primer and the accuracy of initiation by the R2 reverse transcriptase

R2 non-long terminal repeat retrotransposable elements insert at a unique site in the 28S rRNA genes of insects. The protein encoded by the single open reading frame of R2 is capable of conducting the initial steps of its integration in vitro. The protein nicks the noncoding strand of the 28S target DNA (the strand which serves as a template for RNA synthesis) and uses the 3' hydroxyl group exposed by this nick to prime reverse transcription of the R2 RNA template. This target-primed reverse transcription (TPRT) reaction requires that the RNA template contains the 250-nucleotide 3' untranslated region of the R2 element. If this RNA template ends at the precise 3' end of the R2 element, then extra nucleotides, which we refer to as nontemplated nucleotides, are added to the target before cDNA synthesis. The presence of downstream 28S gene sequences on the RNA template reduces the total efficiency but eliminates these nontemplated additions, resulting in nearly 90% of all TPRT products reproducing the 3' junctions seen in vivo. Templates with 5 to 10 nucleotides of the 28S sequence are used most efficiently in this in vitro TPRT reaction. The requirement for downstream 28S rRNA sequences probably explains why the R2 elements of most insects differ from the majority of non-long terminal repeat retrotransposons in that they do not contain an A-rich repeat at their 3' junction with the target DNA. The presence of downstream sequences on these in vitro R2 templates also revealed that the R2 reverse transcriptase can prime cDNA synthesis by using the 3' end of another RNA molecule. This RNA-primed cDNA synthesis is not based on sequence complementarity between the RNA primer and the R2 template. The ability to use the 3' end of a noncomplementary RNA molecule has also been seen with the reverse transcriptase of the mitochondrial Mauriceville plasmid of Neurospora crassa.

Evolution of multiple families of non-LTR retrotransposons in phlebotomine sandflies

In this paper we report on the diversity and distribution of a set of non-LTR retrotransposon (RTP) reverse transcriptase (RT) sequences isolated from phlebotomine sandflies, and their potential for investigating the evolutionary histories of members of this subfamily of flies (Diptera:Psychodidae, Phlebotominae). The phlebotomine RT sequence families derived from one species were as different from each other as they were from RT sequences derived from other species. When each was used to probe Southern blots of sandfly genomic DNA they hybridized only to the species of source and, usually, to others of the same subgenus, but not to DNA from other subgenera-a hybridization pattern consistent with vertical evolution. There was considerable intraspecific variation in hybridization pattern, suggesting the RTs were part of non-LTR RTPs that are (or were recently) subject to flux in genomic position and copy number. Most of the RT families detected in phlebotomines are monophyletic with respect to previously described RTs, and all are monophyletic with RTs of the F/Jockey (Drosophila melanogaster) type of RTP. Orthologous sequences were isolated from the closely related species Phlebotomus perniciosus and P. tobbi (subgenus Larroussius), and different populations of P. perniciosus. The level of sequence divergence among these orthologous RTs, the subgeneric distribution of each RT family, and the intraspecific variation in hybridization pattern of many of them, indicate this class of sequence will provide genetic markers at the sub-generic level.

Evolutionary conservation and molecular characteristics of repetitive sequences of Drosophila koepferae

Thirteen middle repetitive DNA clones obtained from the genome of Drosophila koepferae have been tested for their evolutionary conservation in the other seven species of the buzzatii and martensis clusters (repleta group). All but two of these clones exhibit qualitatively similar patterns of hybridization in the eight species. The average interspecific hybridization signal is 85 per cent of that found intraspecifically, ranging from 73 to 93 per cent. Partial sequencing of six of these clones has shown sequences related to the retrotransposon Gypsy, first characterized in D. melanogaster, as well as to the Anopheles gambiae LINE elements T1Ag and Q. A fragment of a hitherto unknown, short inverted repeat transposable element has also been found. The evolutionary conservation of repetitive D. koepferae sequences seems to be related to the high proportion of simple DNA and inactive mobile elements in the genome of this species.

Integration of the yeast retrotransposon Ty1 is targeted to regions upstream of genes transcribed by RNA polymerase III

Retroviruses and their relatives, the LTR-containing retrotransposons, integrate newly replicated cDNA copies of their genomes into the genomes of their hosts using element-encoded integrases. Although target site selection is not well understood for this general class of elements, it is becoming clear that some elements target their integration events to very specific regions of their host genomes. Evidence is accumulating that the yeast retrotransposon Ty1 behaves in this manner. Ty1 is found frequently adjacent to tRNA genes in the yeast genome and experimental evidence implicates these regions as preferred integration sites. To determine the basis for Ty1 targeting, we developed an in vivo integration assay using a Ty1 donor plasmid and a second target plasmid that could be used to measure the relative frequency of Ty1 integration into sequences cloned from various regions of the yeast genome. Targets containing genes transcribed by RNA polymerase III (Pol III) were up to several hundredfold more active as integration targets than "cold" sequences lacking such genes. High-frequency targeting was dependent on Pol III transcription, and integration was "region specific," occurring exclusively upstream of the transcription start sites of these genes. Thus, Ty1 has evolved a powerful targeting mechanism, requiring Pol III transcription to integrate its DNA at very specific locations within the yeast genome.

R4, a non-LTR retrotransposon specific to the large subunit rRNA genes of nematodes

A 4.7 kb sequence-specific insertion in the 26S ribosomal RNA gene of Ascaris lumbricoides, named R4, is shown to be a non-long terminal repeat (non-LTR) retrotransposable element. The R4 element inserts at a site in the large subunit rRNA gene which is midway between two other sequence-specific non-LTR retrotransposable elements, R1 and R2, found in most insect species. Based on the structure of its open reading frame and the sequence of its reverse transcriptase domain, R4 elements do not appear to be a family of R1 or R2 elements that have changed their insertion site. R4 is most similar in structure and in sequence to the element Dong, which is not specialized for insertion into rRNA units. Thus R4 represents a separate non-LTR retrotransposable element that has become specialized for insertion in the rRNA genes of its host. Using oligonucleotide primers directed to a conserved region of the reverse transcriptase encoding domain, insertions in the R4 site were also amplified from Parascaris equorum and Haemonchus contortus. Why several non-LTR retrotransposable elements have become specialized for insertion into a short (87 bp) region of the large subunit rRNA gene is discussed.

Isolation and characterization of a retroelement from B chromosome (PSR) in the parasitic wasp Nasonia vitripennis

Molecular characterization of the paternal-sex-ratio (PSR) chromosome in Nasonia vitripennis (Hymenoptera: Pteromalidae) has led to the isolation of a dispersed repetitive element. The element is a LTR-containing retrotransposon which has been named NATE (NAsonia Transposable Element). NATE has direct terminal repeats and has an internal amino acid sequence similar to reverse transcriptases of other retroelements. Phylogenetic analysis indicates NATE is a member of the Gypsy/Ty3 group of retrotransposons, and represents the first isolated from Hymenoptera. Five closely related copies of NATE were isolated from the PSR chromosome, but cross-hybridizing elements were not detected on the autosomes of N. vitripennis. Strongly cross-hybridizing elements were, however, detected in two other Nasonia species. This observed distribution of NATE is interesting, because the supernumerary PSR chromosome may be derived from the genome of a sibling species of N. vitripennis.