Sequence relationship of retrotransposable elements R1 and R2 within and between divergent insect species

Hybridogenesis and a potential case of R2 non-LTR retrotransposon horizontal transmission in Bacillus stick insects (Insecta Phasmida)

Horizontal transfer (HT) is an event in which the genetic material is transferred from one species to another, even if distantly related, and it has been demonstrated as a possible essential part of the lifecycle of transposable elements (TEs). However, previous studies on the non-LTR R2 retrotransposon, a metazoan-wide distributed element, indicated its vertical transmission since the Radiata-Bilateria split. Here we present the first possible instances of R2 HT in stick insects of the genus Bacillus (Phasmida). Six R2 elements were characterized in the strictly bisexual subspecies B. grandii grandii, B. grandii benazzii and B. grandii maretimi and in the obligatory parthenogenetic taxon B. atticus. These elements were compared with those previously retrieved in the facultative parthenogenetic species B. rossius. Phylogenetic inconsistencies between element and host taxa, and age versus divergence analyses agree and support at least two HT events. These HT events can be explained by taking into consideration the complex Bacillus reproductive biology, which includes also hybridogenesis, gynogenesis and androgenesis. Through these non-canonical reproductive modes, R2 elements may have been transferred between Bacillus genomes. Our data suggest, therefore, a possible role of hybridization for TEs survival and the consequent reshaping of involved genomes.

Endonuclease domain of non-LTR retrotransposons: loss-of-function mutants and modeling of the R2Bm endonuclease

Non-LTR retrotransposons are an important class of mobile elements that insert into host DNA by target-primed reverse transcription (TPRT). Non-LTR retrotransposons must bind to their mRNA, recognize and cleave their target DNA, and perform TPRT at the site of DNA cleavage. As DNA binding and cleavage are such central parts of the integration reaction, a better understanding of the endonuclease encoded by non-LTR retrotransposons is needed. This paper explores the R2 endonuclease domain from Bombyx mori using in vitro studies and in silico modeling. Mutations in conserved sequences located across the putative PD-(D/E)XK endonuclease domain reduced DNA cleavage, DNA binding and TPRT. A mutation at the beginning of the first α-helix of the modeled endonuclease obliterated DNA cleavage and greatly reduced DNA binding. It also reduced TPRT when tested on pre-cleaved DNA substrates. The catalytic K was located to a non-canonical position within the second α-helix. A mutation located after the fourth β-strand reduced DNA binding and cleavage. The motifs that showed impaired activity form an extensive basic region. The R2 biochemical and structural data are compared and contrasted with that of two other well characterized PD-(D/E)XK endonucleases, restriction endonucleases and archaeal Holliday junction resolvases.

Integration, Regulation, and Long-Term Stability of R2 Retrotransposons

R2 elements are sequence specific non-LTR retrotransposons that exclusively insert in the 28S rRNA genes of animals. R2s encode an endonuclease that cleaves the insertion site and a reverse transcriptase that uses the cleaved DNA to prime reverse transcription of the R2 transcript, a process termed target primed reverse transcription. Additional unusual properties of the reverse transcriptase as well as DNA and RNA binding domains of the R2 encoded protein have been characterized. R2 expression is through co-transcription with the 28S gene and self-cleavage by a ribozyme encoded at the R2 5' end. Studies in laboratory stocks and natural populations of Drosophila suggest that R2 expression is tied to the distribution of R2-inserted units within the rDNA locus. Most individuals have no R2 expression because only a small fraction of their rRNA genes need to be active, and a contiguous region of the locus free of R2 insertions can be selected for activation. However, if the R2-free region is not large enough to produce sufficient rRNA, flanking units - including those inserted with R2 - must be activated. Finally, R2 copies rapidly turnover within the rDNA locus, yet R2 has been vertically maintained in animal lineages for hundreds of millions of years. The key to this stability is R2's ability to remain dormant in rDNA units outside the transcribed regions for generations until the stochastic nature of the crossovers that drive the concerted evolution of the rDNA locus inevitably reshuffle the inserted and uninserted units, resulting in transcription of the R2-inserted units.

Non-LTR R2 element evolutionary patterns: phylogenetic incongruences, rapid radiation and the maintenance of multiple lineages

Retrotransposons of the R2 superclade specifically insert within the 28S ribosomal gene. They have been isolated from a variety of metazoan genomes and were found vertically inherited even if their phylogeny does not always agree with that of the host species. This was explained with the diversification/extinction of paralogous lineages, being proved the absence of horizontal transfer. We here analyze the widest available collection of R2 sequences, either newly isolated from recently sequenced genomes or drawn from public databases, in a phylogenetic framework. Results are congruent with previous analyses, but new important issues emerge. First, the N-terminal end of the R2-B clade protein, so far unknown, presents a new zinc fingers configuration. Second, the phylogenetic pattern is consistent with an ancient, rapid radiation of R2 lineages: being the estimated time of R2 origin (850–600 Million years ago) placed just before the metazoan Cambrian explosion, the wide element diversity and the incongruence with the host phylogeny could be attributable to the sudden expansion of available niches represented by host’s 28S ribosomal genes. Finally, we detect instances of coexisting multiple R2 lineages showing a non-random phylogenetic pattern, strongly similar to that of the “library” model known for tandem repeats: a collection of R2s were present in the ancestral genome and then differentially activated/repressed in the derived species. Models for activation/repression as well as mechanisms for sequence maintenance are also discussed within this framework.

A population genetic model for the maintenance of R2 retrotransposons in rRNA gene loci

R2 retrotransposable elements exclusively insert into the tandemly repeated rRNA genes, the rDNA loci, of their animal hosts. R2 elements form stable long-term associations with their host, in which all individuals in a population contain many potentially active copies, but only a fraction of these individuals show active R2 retrotransposition. Previous studies have found that R2 RNA transcripts are processed from a 28S co-transcript and that the likelihood of R2-inserted units being transcribed is dependent upon their distribution within the rDNA locus. Here we analyze the rDNA locus and R2 elements from nearly 100 R2-active and R2-inactive individuals from natural populations of Drosophila simulans. Along with previous findings concerning the structure and expression of the rDNA loci, these data were incorporated into computer simulations to model the crossover events that give rise to the concerted evolution of the rRNA genes. The simulations that best reproduce the population data assume that only about 40 rDNA units out of the over 200 total units are actively transcribed and that these transcribed units are clustered in a single region of the locus. In the model, the host establishes this transcription domain at each generation in the region with the fewest R2 insertions. Only if the host cannot avoid R2 insertions within this 40-unit domain are R2 elements active in that generation. The simulations also require that most crossover events in the locus occur in the transcription domain in order to explain the empirical observation that R2 elements are seldom duplicated by crossover events. Thus the key to the long-term stability of R2 elements is the stochastic nature of the crossover events within the rDNA locus, and the inevitable expansions and contractions that introduce and remove R2-inserted units from the transcriptionally active domain.

Selfish transposable elements survive in eukaryotic genomes despite the elaborate mechanisms developed by the hosts to limit their activity. One accessible system that simplifies the complex interactions between element and host involves the R2 elements, which exclusively insert in the tandemly arranged rRNA genes. R2 exhibits remarkable stability in animal lineages even though each insertion inactivates one rRNA gene. Here we determine the size of the rDNA locus and R2 number in natural isolates of Drosophila simulans. Combined with previous data concerning the expression and regulation of R2, we develop a detailed population genetic model for rRNA gene and R2 evolution that duplicates all properties of the rRNA loci in natural populations. Critical components of the model are that only a contiguous 40 unit array of rRNA gene units are needed for transcription, that R2 elements are active only when present in this transcription domain, and that most of the crossovers in the rDNA loci occur in this domain. These results suggest that the key to the long-term survival of R2 is the redistribution of rDNA units in the locus brought about by the crossovers that maintain sequence identity in all rDNA units.

B chromosome in the beetle Coprophanaeus cyanescens (Scarabaeidae): emphasis in the organization of repetitive DNA sequences

BACKGROUND

To contribute to the knowledge of coleopteran cytogenetics, especially with respect to the genomic content of B chromosomes, we analyzed the composition and organization of repetitive DNA sequences in the Coprophanaeus cyanescens karyotype. We used conventional staining and the application of fluorescence in situ hybridization (FISH) mapping using as probes C0t-1 DNA fraction, the 18S and 5S rRNA genes, and the LOA-like non-LTR transposable element (TE).

RESULTS

The conventional analysis detected 3 individuals (among 50 analyzed) carrying one small metacentric and mitotically unstable B chromosome. The FISH analysis revealed a pericentromeric block of C0t-1 DNA in the B chromosome but no 18S or 5S rDNA clusters in this extra element. Using the LOA-like TE probe, the FISH analysis revealed large pericentromeric blocks in eight autosomal bivalents and in the B chromosome, and a pericentromeric block extending to the short arm in one autosomal pair. No positive hybridization signal was observed for the LOA-like element in the sex chromosomes.

CONCLUSIONS

The results indicate that the origin of the B chromosome is associated with the autosomal elements, as demonstrated by the hybridization with C0t-1 DNA and the LOA-like TE. The present study is the first report on the cytogenetic mapping of a TE in coleopteran chromosomes. These TEs could have been involved in the origin and evolution of the B chromosome in C. cyanescens.

Gypsy, RTE and Mariner transposable elements populate Eyprepocnemis plorans genome

We analyze here the presence and abundance of three types of transposable elements (TEs), i.e. Gypsy, RTE and Mariner, in the genome of the grasshopper Eyprepocnemis plorans. PCR experiments allowed amplification, cloning and sequencing of these elements (EploGypI, EploRTE5, EploMar20) from the E. plorans genome. Fluorescent in situ hybridization (FISH) showed that all three elements are restricted to euchromatic regions, thus being absent from the pericentromeric region of all A chromosomes, which contain a satellite DNA (satDNA) and ribosomal DNA (rDNA), and being very scarce in B chromosomes mostly made up of these two types of repetitive DNA. FISH suggested that EploGypI is the most abundant and EploMar20 is the least abundant, with EploRTE5 showing intermediate abundance. An estimation of copy number, by means of quantitative PCR, showed that EploGypI is, by far, the most abundant element, followed by EploRTE5 and EploMar20, in consistency with FISH results. RNA isolation and PCR experiments on complementary DNA (cDNA) showed the presence of transcripts for the three TE elements. The implications of the preferential location of these TE elements into euchromatin, the significance of TE abundance in the giant genome of this species, and a possible relationship between TEs and B chromosome mutability, are discussed.

Ribosomal DNA organization before and after magnification in Drosophila melanogaster

In all eukaryotes, the ribosomal RNA genes are stably inherited redundant elements. In Drosophila melanogaster, the presence of a Ybb(-) chromosome in males, or the maternal presence of the Ribosomal exchange (Rex) element, induces magnification: a heritable increase of rDNA copy number. To date, several alternative classes of mechanisms have been proposed for magnification: in situ replication or extra-chromosomal replication, either of which might act on short or extended strings of rDNA units, or unequal sister chromatid exchange. To eliminate some of these hypotheses, none of which has been clearly proven, we examined molecular-variant composition and compared genetic maps of the rDNA in the bb(2) mutant and in some magnified bb(+) alleles. The genetic markers used are molecular-length variants of IGS sequences and of R1 and R2 mobile elements present in many 28S sequences. Direct comparison of PCR products does not reveal any particularly intensified electrophoretic bands in magnified alleles compared to the nonmagnified bb(2) allele. Hence, the increase of rDNA copy number is diluted among multiple variants. We can therefore reject mechanisms of magnification based on multiple rounds of replication of short strings. Moreover, we find no changes of marker order when pre- and postmagnification maps are compared. Thus, we can further restrict the possible mechanisms to two: replication in situ of an extended string of rDNA units or unequal exchange between sister chromatids.

Fosmid library end sequencing reveals a rarely known genome structure of marine shrimp Penaeus monodon

Background

The black tiger shrimp (Penaeus monodon) is one of the most important aquaculture species in the world, representing the crustacean lineage which possesses the greatest species diversity among marine invertebrates. Yet, we barely know anything about their genomic structure. To understand the organization and evolution of the P. monodon genome, a fosmid library consisting of 288,000 colonies and was constructed, equivalent to 5.3-fold coverage of the 2.17 Gb genome. Approximately 11.1 Mb of fosmid end sequences (FESs) from 20,926 non-redundant reads representing 0.45% of the P. monodon genome were obtained for repetitive and protein-coding sequence analyses.

Results

We found that microsatellite sequences were highly abundant in the P. monodon genome, comprising 8.3% of the total length. The density and the average length of microsatellites were evidently higher in comparison to those of other taxa. AT-rich microsatellite motifs, especially poly (AT) and poly (AAT), were the most abundant. High abundance of microsatellite sequences were also found in the transcribed regions. Furthermore, via self-BlastN analysis we identified 103 novel repetitive element families which were categorized into four groups, i.e., 33 WSSV-like repeats, 14 retrotransposons, 5 gene-like repeats, and 51 unannotated repeats. Overall, various types of repeats comprise 51.18% of the P. monodon genome in length. Approximately 7.4% of the FESs contained protein-coding sequences, and the Inhibitor of Apoptosis Protein (IAP) gene and the Innexin 3 gene homologues appear to be present in high abundance in the P. monodon genome.

Conclusions

The redundancy of various repeat types in the P. monodon genome illustrates its highly repetitive nature. In particular, long and dense microsatellite sequences as well as abundant WSSV-like sequences highlight the uniqueness of genome organization of penaeid shrimp from those of other taxa. These results provide substantial improvement to our current knowledge not only for shrimp but also for marine crustaceans of large genome size.

Molecular and cytological characterization of repetitive DNA sequences from the centromeric heterochromatin of Sciara coprophila

Sciara coprophila (Diptera, Nematocera) constitutes a classic model to analyze unusual chromosome behavior such as the somatic elimination of paternal X chromosomes, the elimination of the whole paternal, plus non-disjunction of the maternal X chromosome at male meiosis. The molecular organization of the heterochromatin in S. coprophila is mostly unknown except for the ribosomal DNA located in the X chromosome pericentromeric heterochromatin. The characterization of the centromeric regions, thus, is an essential and required step for the establishment of S. coprophila as a model system to study fundamental mechanisms of chromosome segregation. To accomplish such a study, heterochromatic sections of the X chromosome centromeric region from salivary glands polytene chromosomes were microdissected and microcloned. Here, we report the identification and characterization of two tandem repeated DNA sequences from the pericentromeric region of the X chromosome, a pericentromeric RTE element and an AT-rich centromeric satellite. These sequences will be important tools for the cloning of S. coprophila centromeric heterochromatin using libraries of large genomic clones.

R2 dynamics in Triops cancriformis (Bosc, 1801) (Crustacea, Branchiopoda, Notostraca): turnover rate and 28S concerted evolution

The R2 retrotransposon is here characterized in bisexual populations of the European crustacean Triops cancriformis. The isolated element matches well with the general aspects of the R2 family and it is highly differentiated from that of the congeneric North American Triops longicaudatus. The analysis of 5' truncations indicates that R2 dynamics in T. cancriformis populations show a high turnover rate as observed in Drosophila simulans. For the first time in the literature, though, individuals harboring truncation variants, but lacking the complete element, are found. Present results suggest that transposition-mediated deletion mechanisms, possibly involving genomic turnover processes acting on rDNAs, can dramatically decrease the copy number or even delete R2 from the ribosomal locus. The presence of R2 does not seem to impact on the nucleotide variation of inserted 28S rDNA with respect to the uninserted genes. On the other hand, a low level of polymorphism characterizes rDNA units because new 28S variants continuously spread across the ribosomal array. Again, the interplay between transposition-mediated deletion and molecular drive may explain this pattern.

Origin of nascent lineages and the mechanisms used to prime second-strand DNA synthesis in the R1 and R2 retrotransposons of Drosophila

Comparative analysis of 12 Drosophila genomes reveals insights into the evolution and mechanism of integration of R1 and R2 retrotransposons.

Background

Most arthropods contain R1 and R2 retrotransposons that specifically insert into the 28S rRNA genes. Here, the sequencing reads from 12 Drosophila genomes have been used to address two questions concerning these elements. First, to what extent is the evolution of these elements subject to the concerted evolution process that is responsible for sequence homogeneity among the different copies of rRNA genes? Second, how precise are the target DNA cleavages and priming of DNA synthesis used by these elements?

Results

Most copies of R1 and R2 in each species were found to exhibit less than 0.2% sequence divergence. However, in many species evidence was obtained for the formation of distinct sublineages of elements, particularly in the case of R1. Analysis of the hundreds of R1 and R2 junctions with the 28S gene revealed that cleavage of the first DNA strand was precise both in location and the priming of reverse transcription. Cleavage of the second DNA strand was less precise within a species, differed between species, and gave rise to variable priming mechanisms for second strand synthesis.

Conclusions

These findings suggest that the high sequence identity amongst R1 and R2 copies is because all copies are relatively new. However, each active element generates its own independent lineage that can eventually populate the locus. Independent lineages occur more often with R1, possibly because these elements contain their own promoter. Finally, both R1 and R2 use imprecise, rapidly evolving mechanisms to cleave the second strand and prime second strand synthesis.

The R2 mobile element of Rhynchosciara americana: molecular, cytological and dynamic aspects

Ribosomal RNA genes are encoded by large units clustered (18S, 5S, and 28S) in the nucleolar organizer region in several organisms. Sometimes additional insertions are present in the coding region for the 28S rDNA. These insertions are specific non-long terminal repeat retrotransposons that have very restricted integration targets within the genome. The retrotransposon present in the genome of Rhynchosciara americana, RaR2, was isolated by the screening of a genomic library. Sequence analysis showed the presence of conserved regions, such as a reverse transcriptase domain and a zinc finger motif in the amino terminal region. The insertion site was highly conserved in R. americana and a phylogenetic analysis showed that this element belongs to the R2 clade. The chromosomal localization confirmed that the RaR2 mobile element was inserted into a specific site in the rDNA gene. The expression level of RaR2 in salivary glands during larval development was determined by quantitative RT-PCR, and the increase of relative expression in the 3P of the fourth instar larval could be related to intense gene activity characteristic of this stage. 5'-Truncated elements were identified in different DNA samples. Additionally, in three other Rhynchosciara species, the R2 element was present as a full-length element.

The pattern of R2 retrotransposon activity in natural populations of Drosophila simulans reflects the dynamic nature of the rDNA locus

The pattern and frequency of insertions that enable transposable elements to remain active in a population are poorly understood. The retrotransposable element R2 exclusively inserts into the 28S rRNA genes where it establishes long-term, stable relationships with its animal hosts. Previous studies with laboratory stocks of Drosophila simulans have suggested that control over R2 retrotransposition resides within the rDNA loci. In this report, we sampled 180 rDNA loci of animals collected from two natural populations of D. simulans. The two populations were found to have similar patterns of R2 activity. About half of the rDNA loci supported no or very low levels of R2 transcripts with no evidence of R2 retrotransposition. The remaining half of the rDNA loci had levels of R2 transcripts that varied in a continuous manner over almost a 100-fold range and did support new retrotransposition events. Structural analysis of the rDNA loci in 18 lines that spanned the range of R2 transcript levels in these populations revealed that R2 number and rDNA locus size varied 2-fold; however, R2 activity was not readily correlated with either of these parameters. Instead R2 activity was best correlated with the distribution of elements within the rDNA locus. Loci with no activity had larger contiguous blocks of rDNA units free of R2-insertions. These data suggest a model in which frequent recombination within the rDNA locus continually redistributes R2-inserted units resulting in changing levels of R2 activity within individual loci and persistent R2 activity within the population.

Ribosomal RNA gene insertions in the R2 site of Rhynchosciara (Diptera: Sciaridae)

Ribosomal RNA genes of most insects are interrupted by R1/R2 retrotransposons. The occurrence of R2 retrotransposons in sciarid genomes was studied by PCR and Southern blot hybridization in three Rhynchosciara species and in Trichosia pubescens. Amplification products with the expected size for non-truncated R2 elements were only obtained in Rhynchosciara americana. The rDNA in this species is located in the proximal end of the X mitotic chromosome but in the salivary gland is associated with all four polytene chromosomes. Approximately 50% of the salivary gland rDNA of most R. americana larval groups analysed had an insertion in the R2 site, while no evidence for the presence of R1 elements was found. In-situ hybridization results showed that rDNA repeat units containing R2 take part in the structure of the extrachromosomal rDNA. Also, rDNA resistance to Bal 31 digestion could be interpreted as evidence for nonlinear rDNA as part of the rDNA in the salivary gland. Insertions in the rDNA of three other sciarid species were not detected by Southern blot and in-situ hybridization, suggesting that rDNA retrotransposons are significantly under-represented in their genomes in comparison with R. americana. R2 elements apparently restricted to R. americana correlate with an increased amount of repetitive DNA in its genome in contrast to other Rhynchosciara species. The results obtained in this work together with previous results suggest that evolutionary changes in the genus Rhynchosciara occurred by differential genomic occupation not only of satellite DNA but possibly also of rDNA retrotransposons.

Identifying the predator complex of Homalodisca vitripennis (Hemiptera: Cicadellidae): a comparative study of the efficacy of an ELISA and PCR gut content assay

A growing number of ecologists are using molecular gut content assays to qualitatively measure predation. The two most popular gut content assays are immunoassays employing pest-specific monoclonal antibodies (mAb) and polymerase chain reaction (PCR) assays employing pest-specific DNA. Here, we present results from the first study to simultaneously use both methods to identify predators of the glassy winged sharpshooter (GWSS), Homalodisca vitripennis (Germar) (Hemiptera: Cicadellidae). A total of 1,229 arthropod predators, representing 30 taxa, were collected from urban landscapes in central California and assayed first by means of enzyme-linked immunosorbent assay (ELISA) using a GWSS egg-specific mAb and then by PCR using a GWSS-specific DNA marker that amplifies a 197-base pair fragment of its cytochrome oxidase gene (subunit I). The gut content analyses revealed that GWSS remains were present in 15.5% of the predators examined, with 18% of the spiders and 11% of the insect predators testing positive. Common spider predators included members of the Salticidae, Clubionidae, Anyphaenidae, Miturgidae, and Corinnidae families. Common insect predators included lacewings (Neuroptera: Chrysopidae), praying mantis (Mantodea: Mantidae), ants (Hymenoptera: Formicidae), assassin bugs (Hemiptera: Reduviidae), and damsel bugs (Hemiptera: Nabidae). Comparison of the two assays indicated that they were not equally effective at detecting GWSS remains in predator guts. The advantages of combining the attributes of both types of assays to more precisely assess field predation and the pros and cons of each assay for mass-screening predators are discussed.

CR1 clade of non-LTR retrotransposons from Maculinea butterflies (Lepidoptera: Lycaenidae): evidence for recent horizontal transmission

Background

Non-long terminal repeat (non-LTR) retrotransposons are mobile genetic elements that propagate themselves by reverse transcription of an RNA intermediate. Non-LTR retrotransposons are known to evolve mainly via vertical transmission and random loss. Horizontal transmission is believed to be a very rare event in non-LTR retrotransposons. Our knowledge of distribution and diversity of insect non-LTR retrotransposons is limited to a few species – mainly model organisms such as dipteran genera Drosophila, Anopheles, and Aedes. However, diversity of non-LTR retroelements in arthropods seems to be much richer. The present study extends the analysis of non-LTR retroelements to CR1 clade from four butterfly species of genus Maculinea (Lepidoptera: Lycaenidae).

The lycaenid genus Maculinea, the object of interest for evolutionary biologists and also a model group for European biodiversity studies, possesses a unique, specialized myrmecophilous lifestyle at larval stage. Their caterpillars, after three weeks of phytophagous life on specific food plants drop to the ground where they are adopted to the ant nest by Myrmica foraging workers.

Results

We found that the genome of Maculinea butterflies contains multiple CR1 lineages of non-LTR retrotransposons, including those from MacCR1A, MacCR1B and T1Q families. A comparative analysis of RT nucleotide sequences demonstrated an extremely high similarity among elements both in interspecific and intraspecific comparisons. CR1A-like elements were found only in family Lycaenidae. In contrast, MacCR1B lineage clones were extremely similar to CR1B non-LTR retrotransposons from Bombycidae moths: silkworm Bombyx mori and Oberthueria caeca.

Conclusion

The degree of coding sequence similarity of the studied elements, their discontinuous distribution, and results of divergence-versus-age analysis make it highly unlikely that these sequences diverged at the same time as their host taxa. The only reasonable alternative explanation is horizontal transfer. In addition, phylogenetic markers for population analysis of Maculinea could be developed based on the described non-LTR retrotransposons.

HeT-A and TART, two Drosophila retrotransposons with a bona fide role in chromosome structure for more than 60 million years

Drosophila telomeres have been maintained by retrotransposition for at least 60 MY, which predates the separation of extant species of this genus. Studies of D. melanogaster, D. yakuba, and D. virilis show that, in Drosophila, telomeres are composed of two non-LTR retrotransposons, HeT-A and TART. Far from being static, HeT-A and TART evolve faster than Drosophila euchromatic genes. In spite of their high rate of sequence change, HeT-A and TART maintain their basic structures and unusual individual features. The maintenance of their separate identities suggests that HeT-A and TART cooperate either in the process of retrotransposition onto the chromosome end, or in the formation of telomere chromatin by transposed DNA copies. The telomeric retrotransposons and the Drosophila genome constitute an example of a robust symbiotic relationship between mobile elements and the genome.

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.

An ancestral MADS-box gene duplication occurred before the divergence of plants and animals

Changes in genes encoding transcriptional regulators can alter development and are important components of the molecular mechanisms of morphological evolution. MADS-box genes encode transcriptional regulators of diverse and important biological functions. In plants, MADS-box genes regulate flower, fruit, leaf, and root development. Recent sequencing efforts in Arabidopsis have allowed a nearly complete sampling of the MADS-box gene family from a single plant, something that was lacking in previous phylogenetic studies. To test the long-suspected parallel between the evolution of the MADS-box gene family and the evolution of plant form, a polarized gene phylogeny is necessary. Here we suggest that a gene duplication ancestral to the divergence of plants and animals gave rise to two main lineages of MADS-box genes: TypeI and TypeII. We locate the root of the eukaryotic MADS-box gene family between these two lineages. A novel monophyletic group of plant MADS domains (AGL34 like) seems to be more closely related to previously identified animal SRF-like MADS domains to form TypeI lineage. Most other plant sequences form a clear monophyletic group with animal MEF2-like domains to form TypeII lineage. Only plant TypeII members have a K domain that is downstream of the MADS domain in most plant members previously identified. This suggests that the K domain evolved after the duplication that gave rise to the two lineages. Finally, a group of intermediate plant sequences could be the result of recombination events. These analyses may guide the search for MADS-box sequences in basal eukaryotes and the phylogenetic placement of new genes from other plant species.

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.

Relationships between transposable elements based upon the integrase-transposase domains: is there a common ancestor?

The integrase domain of RNA-mediated elements (class I) and the transposase domain of DNA-mediated transposable elements (class II) were compared. A number of elements contain the DDE signature, which plays an important role in their integration. The possible relationships between mariner-Tc1 and IS elements, retrotransposons, and retroviruses were analyzed from an alignment of this region. The mariner-Tc1 superfamily, and LTR retrotransposons and retroviruses were found to be monophyletic groups. However, the IS elements of bacteria were found in several groups. These results were used to propose an evolutionary history that suggests a common ancestor for some integrases and transposases.

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.

A new non-LTR retrotransposon provides evidence for multiple distinct site-specific elements in Crithidia fasciculata miniexon arrays

We have identified a new member of the family of trypanosome site-specific retrotransposons, using a degenerate oligonucleotide PCR strategy. The 9595 bp element, termed Crithidia retrotransposable element 2 (CRE2), was cloned and found to be inserted in the tandemly arrayed miniexon genes of Crithidia fasciculata. The element is flanked by 29 bp target site duplications but lacks the 3' poly dA tract characteristic of most other non-long terminal repeat retrotransposons. The amino terminal region of the single 2518-codon open reading frame contains a putative metal-binding motif and a proline-rich region similar to gag-like domains of other retrotransposons. The carboxy terminal region of this open reading frame shares sequence homology with the reverse transcriptase and putative endonuclease regions of three previously described trypanosomatid site-specific retrotransposons. All four of these retrotransposons are specifically inserted between nucleotides 11 and 12 of the highly conserved 39mer sequence of the miniexon gene. Most copies of CRE2 and the previously characterized CRE1 are located on different sized chromosomes. Additional CRE-related sequences were identified by screening Crithidia libraries. These results suggest that a particular sequence in the C. fasciculata miniexon repeat is the target for multiple distinct site-specific retrotransposon insertions.

RNA template requirements for target DNA-primed reverse transcription by the R2 retrotransposable element

R2 is a non-long terminal repeat-retrotransposable element that inserts specifically in the 28S rRNA gene of most insects. The single protein encoded by R2 has been shown to contain both site-specific endonuclease and reverse transcriptase activities. Integration of the element involves cleavage of one strand of the 28S target DNA and the utilization of the exposed 3' hydroxyl group to prime the reverse transcription of the R2 RNA transcript. We have characterized the RNA requirement of this target DNA-primed reverse transcription reaction and found that the 250 nucleotides corresponding to the 3' untranslated region of the R2 transcript were necessary and sufficient for the reaction. To investigate the sequence requirements at the site of reverse transcription initiation, a series of RNA templates that contained substitutions and deletions at the extreme 3' end of the RNA were tested. The R2 templates used most efficiently had 3' ends which corresponded to the precise boundary of the R2 element with the 28S gene found in vivo. Transcripts containing short polyadenylated tails (8 nucleotides) were not utilized efficiently. R2 RNAs that were truncated at their 3' ends by 3 to 6 nucleotides were used less efficiently as templates and then only after the R2 reverse transcriptase had added extra, apparently nontemplated, nucleotides to the target DNA. The ability of the reverse transcriptase to add additional nucleotides to the target DNA before engaging the RNA template might be a mechanism for the generation of poly(A) or simple repeat sequences found at the 3' end of most non-long terminal repeat-retrotransposable elements.

Vertical transmission of the retrotransposable elements R1 and R2 during the evolution of the Drosophila melanogaster species subgroup

R1 and R2 are non-long-terminal repeat retrotransposable elements that insert into specific sequences of insect 28S ribosomal RNA genes. These elements have been extensively described in Drosophila melanogaster. To determine whether these elements have been horizontally or vertically transmitted, we characterized R1 and R2 elements from the seven other members of the melanogaster species subgroup by genomic blotting and nucleotide sequencing. Each species was found to have homogeneous families of R1 and R2 elements with the exception of erecta and orena, which have no R2 elements. The DNA sequences of multiple R1 and R2 copies from each species indicated nucleotide divergence within each species averaged only 0.48% for R1 and 0.35% for R2, well below the level of divergence among the species. Most copies of R1 and R2 (40 of 47) sequenced from the seven species were potentially functional, as indicated by the absence of premature termination codons or translational frameshifts that would destroy the open reading frame of the element. The sequence relationships of both the R1 and R2 elements from the various members of the melanogaster subgroup closely followed that of the species phylogeny, suggesting that R1 and R2 have been stably maintained by vertical transmission since the origin of this species subgroup 17-20 million years ago. The remarkable stability of R1 and R2, compared to what has been suggested for transposable elements that insert at multiple locations in these same species, may be due to their unique specificity for sites in the rRNA gene locus. Under low copy number conditions, when it is essential for any mobile element to transpose, the insertion specificities of R1 and R2 ensure uniform developmentally regulated target sites that can be occupied with little or no detrimental effect on the host.