The 3' ends of tRNA-derived SINEs originated from the 3' ends of LINEs: a new example from the bovine genome

Comparison of ruminal microbiota, IL-1β gene variation, and tick incidence between Holstein × Gyr and Holstein heifers in grazing system

Introduction

The variation in bacterial communities among breeds has been previously reported and may be one of the reasons why Holstein × Gyr dairy heifers have better development in grazing systems in tropical conditions. This study aimed to explore the ruminal microbiota composition, the IL-1β gene variation, tick incidence, and blood parameters of Holstein × Gyr (½ Holstein × ½ Gyr) and Holstein heifers grazing intensely managed Guinea grass (Panicum maximum Jacq. cv. Mombaça).

Methods

Sixteen heifers were divided into two groups consisting of 8 Holstein × Gyr and 8 Holstein heifers. The experimental period was comprised of 3 periods of 21 days. Ruminal samples were taken via the stomach tube technique. The sequencing of the V4 hypervariable region of the 16S rRNA gene was performed using the Illumina MiSeq platform. Counting and collection of ticks were conducted each 21 days. Blood and skeletal muscle tissue biopsies were performed at the end of the experiment.

Results

Firmicutes were the most abundant phyla present in both breed rumen samples and Bacteroidota showed differences in relative abundance between breed groups, with greater values for Holstein heifers (p < 0.05 with FDR correction). The 10 most abundant unique OTUs identified in each breed included several OTUs of the genus Prevotella. Holstein heifers had a greater tick count and weight (9.8 ticks/animal and 1.6 g/animal, respectively) than Holstein × Gyr (2.56 ticks/animal and 0.4 g/animal, respectively). We found nucleotide substitutions in the IL-1β gene that might be related to adaptation and resistance phenotypes to tick infestation in Holstein × Gyr heifers. Blood concentrations of urea, albumin, insulin-like growth factor 1, triiodothyronine, and thyroxine were greater in Holstein × Gyr than in Holstein heifers.

Conclusion

Adaptations in Holstein × Gyr heifers such as ruminal microbiota, tick resistance, nucleotide substitutions in IL-1β gene, and hormone concentration suggest a better energy metabolism and thermoregulation resulting in better performance in tropical grazing systems.

How to start a LINE: 5' switching rejuvenates LINE retrotransposons in tobacco and related Nicotiana species

By contrast to their conserved mammalian counterparts, plant long interspersed nuclear elements (LINEs) are highly variable, splitting into many low-copy families. Curiously, LINE families from the retrotransposable element (RTE) clade retain a stronger sequence conservation and hence reach higher copy numbers. The cause of this RTE-typical property is not yet understood, but would help clarify why some transposable elements are removed quickly, whereas others persist in plant genomes. Here, we bring forward a detailed study of RTE LINE structure, diversity and evolution in plants. For this, we argue that the nightshade family is the ideal taxon to follow the evolutionary trajectories of RTE LINEs, given their high abundance, recent activity and partnership to non-autonomous elements. Using bioinformatic, cytogenetic and molecular approaches, we detect 4029 full-length RTE LINEs across the Solanaceae. We finely characterize and manually curate a core group of 458 full-length LINEs in allotetraploid tobacco, show an integration event after polyploidization and trace hybridization by RTE LINE composition of parental genomes. Finally, we reveal the role of the untranslated regions (UTRs) as causes for the unique RTE LINE amplification and evolution pattern in plants. On the one hand, we detected a highly conserved motif at the 3' UTR, suggesting strong selective constraints acting on the RTE terminus. On the other hand, we observed successive rounds of 5' UTR cycling, constantly rejuvenating the promoter sequences. This interplay between exchangeable promoters and conserved LINE bodies and 3' UTR likely allows RTE LINEs to persist and thrive in plant genomes.

The non-LTR retrotransposons of Entamoeba histolytica: genomic organization and biology

Genome sequence analysis of Entamoeba species revealed various classes of transposable elements. While E. histolytica and E. dispar are rich in non-long terminal repeat (LTR) retrotransposons, E. invadens contains predominantly DNA transposons. Non-LTR retrotransposons of E. histolytica constitute three families of long interspersed nuclear elements (LINEs), and their short, nonautonomous partners, SINEs. They occupy ~ 11% of the genome. The EhLINE1/EhSINE1 family is the most abundant and best studied. EhLINE1 is 4.8 kb, with two ORFs that encode functions needed for retrotransposition. ORF1 codes for the nucleic acid-binding protein, and ORF2 has domains for reverse transcriptase (RT) and endonuclease (EN). Most copies of EhLINEs lack complete ORFs. ORF1p is expressed constitutively, but ORF2p is not detected. Retrotransposition could be demonstrated upon ectopic over expression of ORF2p, showing that retrotransposition machinery is functional. The newly retrotransposed sequences showed a high degree of recombination. In transcriptomic analysis, RNA-Seq reads were mapped to individual EhLINE1 copies. Although full-length copies were transcribed, no full-length 4.8 kb transcripts were seen. Rather, sense transcripts mapped to ORF1, RT and EN domains. Intriguingly, there was strong antisense transcription almost exclusively from the RT domain. These unique features of EhLINE1 could serve to attenuate retrotransposition in E. histolytica.

Population analysis of retrotransposons in giraffe genomes supports RTE decline and widespread LINE1 activity in Giraffidae

BACKGROUND

The majority of structural variation in genomes is caused by insertions of transposable elements (TEs). In mammalian genomes, the main TE fraction is made up of autonomous and non-autonomous non-LTR retrotransposons commonly known as LINEs and SINEs (Long and Short Interspersed Nuclear Elements). Here we present one of the first population-level analysis of TE insertions in a non-model organism, the giraffe. Giraffes are ruminant artiodactyls, one of the few mammalian groups with genomes that are colonized by putatively active LINEs of two different clades of non-LTR retrotransposons, namely the LINE1 and RTE/BovB LINEs as well as their associated SINEs. We analyzed TE insertions of both types, and their associated SINEs in three giraffe genome assemblies, as well as across a population level sampling of 48 individuals covering all extant giraffe species.

RESULTS

The comparative genome screen identified 139,525 recent LINE1 and RTE insertions in the sampled giraffe population. The analysis revealed a drastically reduced RTE activity in giraffes, whereas LINE1 is still actively propagating in the genomes of extant (sub)-species. In concert with the extremely low activity of the giraffe RTE, we also found that RTE-dependent SINEs, namely Bov-tA and Bov-A2, have been virtually immobile in the last 2 million years. Despite the high current activity of the giraffe LINE1, we did not find evidence for the presence of currently active LINE1-dependent SINEs. TE insertion heterozygosity rates differ among the different (sub)-species, likely due to divergent population histories.

CONCLUSIONS

The horizontally transferred RTE/BovB and its derived SINEs appear to be close to inactivation and subsequent extinction in the genomes of extant giraffe species. This is the first time that the decline of a TE family has been meticulously analyzed from a population genetics perspective. Our study shows how detailed information about past and present TE activity can be obtained by analyzing large-scale population-level genomic data sets.

Complete CSN1S2 Characterization, Novel Allele Identification and Association With Milk Fatty Acid Composition in River Buffalo

The αs2-casein is one of the phosphoproteins secreted in all ruminants' milk, and it is the most hydrophilic of all caseins. However, this important gene (CSN1S2) has not been characterized in detail in buffaloes with only two alleles detected (reported as alleles A and B), and no association studies with milk traits have been carried out unlike what has been achieved for other species of ruminants. In this study, we sequenced the whole gene of two Mediterranean river buffalo homozygotes for the presence/absence of the nucleotide C (g.7539G>C) realized at the donor splice site of exon 7 and, therefore, responsible for the skipping of the same exon at mRNA level (allele B). A high genetic variability was found all over the two sequenced CSN1S2 alleles. In particular, 74 polymorphic sites were found in introns, six in the promoter, and three SNPs in the coding region (g.11072C>T, g.12803A>T, and g.14067A>G) with two of them responsible for amino acid replacements. Considering this genetic diversity, those found in the database and the SNP at the donor splice site of exon 7, it is possible to deduce at least eight different alleles (CSN1S2 A, B, B1, B2, C, D, E, and F) responsible for seven different possible translations of the buffalo αs2-casein. Haplotype data analysis suggests an evolutionary pathway of buffalo CSN1S2 gene consistent with our proposal that the published allele CSN1S2 A is the ancestral αs2-CN form, and the B2 probably arises from interallelic recombination (single crossing) between the alleles D and B (or B1). The allele CSN1S2 C is of new identification, while CSN1S2 B, B1, and B2 are deleted alleles because all are characterized by the mutation g.7539G>C. Two SNPs (g.7539G>C and g.14067A>G) were genotyped in 747 Italian buffaloes, and major alleles had a relative frequency of 0.83 and 0.51, respectively. An association study between these SNPs and milk traits including fatty acid composition was carried out. The SNP g.14067A>G showed a significant association (P < 0.05) on the content of palmitic acid in buffalo milk, thus suggesting its use in marker-assisted selection programs aiming for the improvement of buffalo milk fatty acid composition.

Existence of Bov-B LINE Retrotransposons in Snake Lineages Reveals Recent Multiple Horizontal Gene Transfers with Copy Number Variation

Transposable elements (TEs) are dynamic elements present in all eukaryotic genomes. They can “jump” and amplify within the genome and promote segmental genome rearrangements on both autosomes and sex chromosomes by disruption of gene structures. The Bovine-B long interspersed nuclear element (Bov-B LINE) is among the most abundant TE-retrotransposon families in vertebrates due to horizontal transfer (HT) among vertebrate lineages. Recent studies have shown multiple HTs or the presence of diverse Bov-B LINE groups in the snake lineage. It is hypothesized that Bov-B LINEs are highly dynamic and that the diversity reflects multiple HTs in snake lineages. Partial sequences of Bov-B LINE from 23 snake species were characterized. Phylogenetic analysis resolved at least two Bov-B LINE groups that might correspond to henophidian and caenophidian snakes; however, the tree topology differed from that based on functional nuclear and mitochondrial gene sequences. Several Bov-B LINEs of snakes showed greater than 80% similarity to sequences obtained from insects, whereas the two Bov-B LINE groups as well as sequences from the same snake species classified in different Bov-B LINE groups showed sequence similarities of less than 80%. Calculation of estimated divergence time and pairwise divergence between all individual Bov-B LINE copies suggest invasion times ranging from 79.19 to 98.8 million years ago in snakes. Accumulation of elements in a lineage-specific fashion ranged from 9 × 10⁻⁶% to 5.63 × 10⁻²% per genome. The genomic proportion of Bov-B LINEs varied among snake species but was not directly associated with genome size or invasion time. No differentiation in Bov-B LINE copy number between males and females was observed in any of the snake species examined. Incongruence in tree topology between Bov-B LINEs and other snake phylogenies may reflect past HT events. Sequence divergence of Bov-B LINEs between copies suggests that recent multiple HTs occurred within the same evolutionary timeframe in the snake lineage. The proportion of Bov-B LINEs varies among species, reflecting species specificity in TE invasion. The rapid speciation of snakes, coinciding with Bov-B LINE invasion in snake genomes, leads us to better understand the effect of Bov-B LINEs on snake genome evolution.

Divergence of 3' ends as a driver of short interspersed nuclear element (SINE) evolution in the Salicaceae

Short interspersed nuclear elements (SINEs) are small, non-autonomous and heterogeneous retrotransposons that are widespread in plants. To explore the amplification dynamics and evolutionary history of SINE populations in representative deciduous tree species, we analyzed the genomes of the six following Salicaceae species: Populus deltoides, Populus euphratica, Populus tremula, Populus tremuloides, Populus trichocarpa, and Salix purpurea. We identified 11 Salicaceae SINE families (SaliS-I to SaliS-XI), comprising 27 077 full-length copies. Most of these families harbor segmental similarities, providing evidence for SINE emergence by reshuffling or heterodimerization. We observed two SINE groups, differing in phylogenetic distribution pattern, similarity and 3' end structure. These groups probably emerged during the 'salicoid duplication' (~65 million years ago) in the Salix-Populus progenitor and during the separation of the genus Salix (45-65 million years ago), respectively. In contrast to conserved 5' start motifs across species and SINE families, the 3' ends are highly variable in sequence and length. This extraordinary 3'-end variability results from mutations in the poly(A) tail, which were fixed by subsequent amplificational bursts. We show that the dissemination of newly evolved 3' ends is accomplished by a displacement of older motifs, leading to various 3'-end subpopulations within the SaliS families.

LINEs Contribute to the Origins of Middle Bodies of SINEs besides 3' Tails

Short interspersed elements (SINEs), which are nonautonomous transposable elements, require the transposition machinery of long interspersed elements (LINEs) to mobilize. SINEs are composed of two or more independently originating parts. The 5′ region is called the “head” and is derived mainly from small RNAs, and the 3′ region (“tail”) originates from the 3′ region of LINEs and is responsible for being recognized by counterpart LINE proteins. The origin of the middle “body” of SINEs is enigmatic, although significant sequence similarities among SINEs from very diverse species have been observed. Here, a systematic analysis of the similarities among SINEs and LINEs deposited on Repbase, a comprehensive database of eukaryotic repeat sequences was performed. Three primary findings are described: 1) The 5′ regions of only two clades of LINEs, RTE and Vingi, were revealed to have contributed to the middle parts of SINEs; 2) The linkage of the 5′ and 3′ parts of LINEs can be lost due to occasional tail exchange of SINEs; and 3) The previously proposed Ceph-domain was revealed to be a fusion of a CORE-domain and a 5′ part of RTE clade of LINE. Based on these findings, a hypothesis that the 5′ parts of bipartite nonautonomous LINEs, which possess only the 5′ and 3′ regions of the original LINEs, can contribute to the undefined middle part of SINEs is proposed.

Evolutionary modes of emergence of short interspersed nuclear element (SINE) families in grasses

Short interspersed nuclear elements (SINEs) are non-autonomous transposable elements which are propagated by retrotransposition and constitute an inherent part of the genome of most eukaryotic species. Knowledge of heterogeneous and highly abundant SINEs is crucial for de novo (or improvement of) annotation of whole genome sequences. We scanned Poaceae genome sequences of six important cereals (Oryza sativa, Triticum aestivum, Hordeum vulgare, Panicum virgatum, Sorghum bicolor, Zea mays) and Brachypodium distachyon to examine the diversity and evolution of SINE populations. We comparatively analyzed the structural features, distribution, evolutionary relation and abundance of 32 SINE families and subfamilies within grasses, comprising 11 052 individual copies. The investigation of activity profiles within the Poaceae provides insights into their species-specific diversification and amplification. We found that Poaceae SINEs (PoaS) fall into two length categories: simple SINEs of up to 180 bp and dimeric SINEs larger than 240 bp. Detailed analysis at the nucleotide level revealed that multimerization of related and unrelated SINE copies is an important evolutionary mechanism of SINE formation. We conclude that PoaS families diversify by massive reshuffling between SINE families, likely caused by insertion of truncated copies, and provide a model for this evolutionary scenario. Twenty-eight of 32 PoaS families and subfamilies show significant conservation, in particular either in the 5' or 3' regions, across Poaceae species and share large sequence stretches with one or more other PoaS families.

Conserved 3' UTR stem-loop structure in L1 and Alu transposons in human genome: possible role in retrotransposition

BACKGROUND

In the process of retrotransposition LINEs use their own machinery for copying and inserting themselves into new genomic locations, while SINEs are parasitic and require the machinery of LINEs. The exact mechanism of how a LINE-encoded reverse transcriptase (RT) recognizes its own and SINE RNA remains unclear. However it was shown for the stringent-type LINEs that recognition of a stem-loop at the 3'UTR by RT is essential for retrotransposition. For the relaxed-type LINEs it is believed that the poly-A tail is a common recognition element between LINE and SINE RNA. However polyadenylation is a property of any messenger RNA, and how the LINE RT recognizes transposon and non-transposon RNAs remains an open question. It is likely that RNA secondary structures play an important role in RNA recognition by LINE encoded proteins.

RESULTS

Here we selected a set of L1 and Alu elements from the human genome and investigated their sequences for the presence of position-specific stem-loop structures. We found highly conserved stem-loop positions at the 3'UTR. Comparative structural analyses of a human L1 3'UTR stem-loop showed a similarity to 3'UTR stem-loops of the stringent-type LINEs, which were experimentally shown to be recognized by LINE RT. The consensus stem-loop structure consists of 5-7 bp loop, 8-10 bp stem with a bulge at a distance of 4-6 bp from the loop. The results show that a stem loop with a bulge exists at the 3'-end of Alu. We also found conserved stem-loop positions at 5'UTR and at the end of ORF2 and discuss their possible role.

CONCLUSIONS

Here we presented an evidence for the presence of a highly conserved 3'UTR stem-loop structure in L1 and Alu retrotransposons in the human genome. Both stem-loops show structural similarity to the stem-loops of the stringent-type LINEs experimentally confirmed as essential for retrotransposition. Here we hypothesize that both L1 and Alu RNA are recognized by L1 RT via the 3'-end RNA stem-loop structure. Other conserved stem-loop positions in L1 suggest their possible functions in protein-RNA interactions but to date no experimental evidence has been reported.

A genome survey sequencing of the Java mouse deer (Tragulus javanicus) adds new aspects to the evolution of lineage specific retrotransposons in Ruminantia (Cetartiodactyla)

Ruminantia, the ruminating, hoofed mammals (cow, deer, giraffe and allies) are an unranked artiodactylan clade. Around 50-60 million years ago the BovB retrotransposon entered the ancestral ruminantian genome through horizontal gene transfer. A survey genome screen using 454-pyrosequencing of the Java mouse deer (Tragulus javanicus) and the lesser kudu (Tragelaphus imberbis) was done to investigate and to compare the landscape of transposable elements within Ruminantia. The family Tragulidae (mouse deer) is the only representative of Tragulina and phylogenetically important, because it represents the earliest divergence in Ruminantia. The data analyses show that, relative to other ruminantian species, the lesser kudu genome has seen an expansion of BovB Long INterspersed Elements (LINEs) and BovB related Short INterspersed Elements (SINEs) like BOVA2. In comparison the genome of Java mouse deer has fewer BovB elements than other ruminants, especially Bovinae, and has in addition a novel CHR-3 SINE most likely propagated by LINE-1. By contrast the other ruminants have low amounts of CHR SINEs but high numbers of actively propagating BovB-derived and BovB-propagated SINEs. The survey sequencing data suggest that the transposable element landscape in mouse deer (Tragulina) is unique among Ruminantia, suggesting a lineage specific evolutionary trajectory that does not involve BovB mediated retrotransposition. This shows that the genomic landscape of mobile genetic elements can rapidly change in any lineage.

Mechanism by which a LINE protein recognizes its 3' tail RNA

LINEs mobilize their own copies via retrotransposition. LINEs can be divided into two types. One is a stringent type, which constitutes a majority of LINEs. The other is a relaxed type. To elucidate the molecular mechanism of retrotransposition, we used here two different zebrafish LINEs belonging to the stringent type. By using retrotransposition assays, we demonstrated that proteins (ORF2) encoded by an individual LINE recognize the cognate 3′ tail sequence of the LINE RNA strictly. By conducting in vitro binding assays with a variety of ORF2 proteins, we demonstrated that the region between the endonuclease and reverse transcriptase domains in ORF2 is the site at which the proteins bind the stem-loop structure of the 3′ tail RNA, showing that the strict recognition of the stem-loop structure by the cognate ORF2 protein is an important step in retrotransposition. This recognition can be bipartite, involving the general recognition of the stem by cTBR (conserved tail-binding region) of ORF2 and the specific recognition of the loop by vTBR (variable tail-binding region). This is the first report that clearly characterized the RNA-binding region in ORF2, providing the generality for the recognition mechanism of the RNA tail by the ORF2 protein encoded by LINEs.

RNA-Mediated Gene Duplication and Retroposons: Retrogenes, LINEs, SINEs, and Sequence Specificity

A substantial number of “retrogenes” that are derived from the mRNA of various intron-containing genes have been reported. A class of mammalian retroposons, long interspersed element-1 (LINE1, L1), has been shown to be involved in the reverse transcription of retrogenes (or processed pseudogenes) and non-autonomous short interspersed elements (SINEs). The 3′-end sequences of various SINEs originated from a corresponding LINE. As the 3′-untranslated regions of several LINEs are essential for retroposition, these LINEs presumably require “stringent” recognition of the 3′-end sequence of the RNA template. However, the 3′-ends of mammalian L1s do not exhibit any similarity to SINEs, except for the presence of 3′-poly(A) repeats. Since the 3′-poly(A) repeats of L1 and Alu SINE are critical for their retroposition, L1 probably recognizes the poly(A) repeats, thereby mobilizing not only Alu SINE but also cytosolic mRNA. Many flowering plants only harbor L1-clade LINEs and a significant number of SINEs with poly(A) repeats, but no homology to the LINEs. Moreover, processed pseudogenes have also been found in flowering plants. I propose that the ancestral L1-clade LINE in the common ancestor of green plants may have recognized a specific RNA template, with stringent recognition then becoming relaxed during the course of plant evolution.

Expansion of CORE-SINEs in the genome of the Tasmanian devil

Background

The genome of the carnivorous marsupial, the Tasmanian devil (Sarcophilus harrisii, Order: Dasyuromorphia), was sequenced in the hopes of finding a cure for or gaining a better understanding of the contagious devil facial tumor disease that is threatening the species’ survival. To better understand the Tasmanian devil genome, we screened it for transposable elements and investigated the dynamics of short interspersed element (SINE) retroposons.

Results

The temporal history of Tasmanian devil SINEs, elucidated using a transposition in transposition analysis, indicates that WSINE1, a CORE-SINE present in around 200,000 copies, is the most recently active element. Moreover, we discovered a new subtype of WSINE1 (WSINE1b) that comprises at least 90% of all Tasmanian devil WSINE1s. The frequencies of WSINE1 subtypes differ in the genomes of two of the other Australian marsupial orders. A co-segregation analysis indicated that at least 66 subfamilies of WSINE1 evolved during the evolution of Dasyuromorphia. Using a substitution rate derived from WSINE1 insertions, the ages of the subfamilies were estimated and correlated with a newly established phylogeny of Dasyuromorphia. Phylogenetic analyses and divergence time estimates of mitochondrial genome data indicate a rapid radiation of the Tasmanian devil and the closest relative the quolls (Dasyurus) around 14 million years ago.

Conclusions

The radiation and abundance of CORE-SINEs in marsupial genomes indicates that they may be a major player in the evolution of marsupials. It is evident that the early phases of evolution of the carnivorous marsupial order Dasyuromorphia was characterized by a burst of SINE activity. A correlation between a speciation event and a major burst of retroposon activity is for the first time shown in a marsupial genome.

SINEs

Short interspersed elements (SINEs) are mobile genetic elements that invade the genomes of many eukaryotes. Since their discovery about 30 years ago, many gaps in our understanding of the biology and function of SINEs have been filled. This review summarizes the past and recent advances in the studies of SINEs. The structure and origin of SINEs as well as the processes involved in their amplification, transcription, RNA processing, reverse transcription, and integration of a SINE copy into the genome are considered. Then we focus on the significance of SINEs for the host genomes. While these genomic parasites can be deleterious to the cell, the long-term being in the genome has made SINEs a valuable source of genetic variation providing regulatory elements for gene expression, alternative splice sites, polyadenylation signals, and even functional RNA genes.

Polypteridae (Actinopterygii: Cladistia) and DANA-SINEs insertions

SINE sequences are interspersed throughout virtually all eukaryotic genomes and greatly outnumber the other repetitive elements. These sequences are of increasing interest for phylogenetic studies because of their diagnostic power for establishing common ancestry among taxa, once properly characterized. We identified and characterized a peculiar family of composite tRNA-derived short interspersed SINEs, DANA-SINEs, associated with mutational activities in Danio rerio, in a group of species belonging to one of the most basal bony fish families, the Polypteridae, in order to investigate their own inner specific phylogenetic relationships. DANA sequences were identified, sequenced and then localized, by means of fluorescent in situ hybridization (FISH), in six Polypteridae species (Polypterus delhezi, P. ornatipinnis, P. palmas, P. buettikoferi P. senegalus and Erpetoichthys calabaricus) After cloning, the sequences obtained were aligned for phylogenetic analysis, comparing them with three Dipnoan lungfish species (Protopterus annectens, P. aethiopicus, Lepidosiren paradoxa), and Lethenteron reissneri (Petromyzontidae)was used as outgroup. The obtained overlapping MP, ML and NJ tree clustered together the species belonging to the two taxonomically different Osteichthyans groups: the Polypteridae, by one side, and the Protopteridae by the other, with the monotypic genus Erpetoichthys more distantly related to the Polypterus genus comprising three distinct groups: P. palmas and P. buettikoferi, P. delhezi and P. ornatipinnis and P. senegalus. In situ hybridization with DANA probes marked along the whole chromosome arms in the metaphases of all the Polypteridae species examined.

5S rRNA-derived and tRNA-derived SINEs in fruit bats

Most short retroposons (SINEs) descend from cellular tRNA of 7SL RNA. Here, four new SINEs were found in megabats (Megachiroptera) but neither in microbats nor in other mammals. Two of them, MEG-RS and MEG-RL, descend from another cellular RNA, 5S rRNA; one (MEG-T2) is a tRNA-derived SINE; and MEG-TR is a hybrid tRNA/5S rRNA SINE. Insertion locus analysis suggests that these SINEs were active in the recent fruit bat evolution. Analysis of MEG-RS and MEG-RL in comparison with other few 5S rRNA-derived SINEs demonstrates that the internal RNA polymerase III promoter is their most invariant region, while the secondary structure is more variable. The mechanisms underlying the modular structure of these and other SINEs as well as their variation are discussed. The scenario of evolution of MEG SINEs is proposed.

Retroposed SNOfall--a mammalian-wide comparison of platypus snoRNAs

Diversification of mammalian species began more than 160 million years ago when the egg-laying monotremes diverged from live bearing mammals. The duck-billed platypus (Ornithorhynchus anatinus) and echidnas are the only potential contemporary witnesses of this period and, thereby, provide a unique insight into mammalian genome evolution. It has become clear that small RNAs are major regulatory agents in eukaryotic cells, and the significant role of non-protein-coding (npc) RNAs in transcription, processing, and translation is now well accepted. Here we show that the platypus genome contains more than 200 small nucleolar (sno) RNAs among hundreds of other diverse npcRNAs. Their comparison among key mammalian groups and other vertebrates enabled us to reconstruct a complete temporal pathway of acquisition and loss of these snoRNAs. In platypus we found cis- and trans-duplication distribution patterns for snoRNAs, which have not been described in any other vertebrates but are known to occur in nematodes. An exciting novelty in platypus is a snoRNA-derived retroposon (termed snoRTE) that facilitates a very effective dispersal of an H/ACA snoRNA via RTE-mediated retroposition. From more than 40,000 detected full-length and truncated genomic copies of this snoRTE, at least 21 are processed into mature snoRNAs. High-copy retroposition via multiple host gene-promoted transcription units is a novel pathway for combining housekeeping function and SINE-like dispersal and reveals a new dimension in the evolution of novel snoRNA function.

Bov-B-mobilized SINEs in vertebrate genomes

Two new short retroposon families (SINEs) have been found in the genome of springhare Pedetes capensis (Rodentia). One of them, Ped-1, originated from 5S rRNA, while the other one, Ped-2, originated from tRNA-derived SINE ID. In contrast to most currently active mammalian SINEs mobilized by L1 long retrotransposon (LINE), Ped-1 and Ped-2 are mobilized by Bov-B, a LINE family of the widely distributed RTE clade. The 3' part of these SINEs originates from two sequences in the 5' and 3' regions of Bov-B. Such bipartite structure of the LINE-derived part has been revealed in all Bov-B-mobilized SINEs known to date (AfroSINE, Bov-tA, Mar-1, and Ped-1/2), which distinguishes them from other SINEs with only a 3' LINE-derived part. Structural analysis and the distribution of Bov-B LINEs and partner SINEs supports the horizontal transfer of Bov-B, while the SINEs emerged independently in lineages with this LINE.

Poxviruses as possible vectors for horizontal transfer of retroposons from reptiles to mammals

Poxviruses (Poxviridae) are a family of double-stranded DNA viruses with no RNA stage. Members of the genus Orthopoxvirus (OPV) are highly invasive and virulent. It was recently shown that the taterapox virus (TATV) from a West African rodent is the sister of camelpox virus and therefore belongs to the clade closest to the variola virus (VARV), the etiological agent of smallpox. Although these OPVs are among the most dreaded pathogens on Earth, our current knowledge of their genomes, their origins, and their possible hosts is still very limited. Here, we report the horizontal transfer of a retroposon (known only from reptilian genomes) to the TATV genome. After isolating and analyzing different subfamilies of short interspersed elements (SINEs) from lizards and snakes, we identified a highly poisonous snake (Echis ocellatus) from West Africa as the closest species from which the SINE sequence discovered in the TATV genome (TATV-SINE) was transferred to the virus. We discovered direct repeats derived from the virus flanking the TATV-SINE, and the absence of any snake-derived DNA flanking the SINE. These data provide strong evidence that the TATV-SINE was actually transferred within the snake to the viral genome by retrotransposition and not by any horizontal transfer at the DNA level. We propose that the snake is another host for TATV, suggesting that VARV-related epidemiologically relevant viruses may have derived from our cold-blooded ancestors and that poxviruses are possible vectors for horizontal transfer of retroposons from reptiles to mammals.

Genomic expansion of the Bov-A2 retroposon relating to phylogeny and breed management

Bov-A2 is a retroposon that is widely distributed among the genomes of ruminants (e.g., cow, deer, giraffe, pronghorn, musk deer, and chevrotain). This retroposon is composed of two monomers, called Bov-A units, which are joined by a linker sequence. The structure and origin of Bov-A2 has been well characterized but a genome-level exploration of this retroposon has not been implemented. In this study we performed an extensive search for Bov-A2 using all available genome sequence data on Bos taurus. We found unique Bov-A2-derived sequences that were longer than Bov-A2 due to amplification of three to six Bov-A units arranged in tandem. Detailed analysis of these elongated Bov-A2-derived sequences revealed that they originated through unequal crossing-over of Bov-A2. We found a large number of these elongated Bov-A2-derived sequences in cattle genomes, indicating that unequal crossing-over of Bov-A2 occurred very frequently. We found that this type of elongation is not observed in wild bovine and is therefore specific to the domesticated cattle genome. Furthermore, at specific loci, the number of Bov-A units was also polymorphic between alleles, implying that the elongation of Bov-A units might have occurred very recently. For these reasons, we speculate that genomic instability in bovine genomes can lead to extensive unequal crossing-over of Bov-A2 and levels of polymorphism might be generated in part by repeated outbreeding.

DNA cleavage and Trp53 differentially affect SINE transcription

Among the cellular responses observed following treatment with DNA-damaging agents is the activation of Short Interspersed Elements (SINEs; retrotransposable genetic elements that comprise over 10% of the human genome). By placing a human SINE (the Alu element) into murine cells, we have previously shown that DNA-damaging agents such as etoposide can induce both upregulation of SINE transcript levels and SINE retrotransposition. A similarly cytotoxic (but not genotoxic) exposure to vincristine was not associated with SINE activation. Here we demonstrate that multiple other genotoxic exposures are associated with upregulation of SINE transcript levels. By comparing the effects of similarly cytotoxic doses of the topoisomerase II inhibitors etoposide and merbarone, we confirm that DNA strand breakage is specifically associated with SINE induction. By evaluating transcription rate and RNA stability, we demonstrate that SINE induction by genotoxic exposure is associated with transcriptional induction and not with transcript stabilization. Finally we demonstrate that SINE induction by genotoxic stress is mediated by a Trp53-independent pathway, and in fact that Trp53 plays an inhibitory role in attenuating the transcriptional induction of SINE elements following exposure to a genotoxic agent. Together these data support a model in which initial DNA damage can trigger genomic instability due to SINE activation, a response which may be amplified in cancer cells lacking functional TP53.

Repetitive sequence environment distinguishes housekeeping genes

Housekeeping genes are expressed across a wide variety of tissues. Since repetitive sequences have been reported to influence the expression of individual genes, we employed a novel approach to determine whether housekeeping genes can be distinguished from tissue-specific genes by their repetitive sequence context. We show that Alu elements are more highly concentrated around housekeeping genes while various longer (>400-bp) repetitive sequences ("repeats"), including Long Interspersed Nuclear Element-1 (LINE-1) elements, are excluded from these regions. We further show that isochore membership does not distinguish housekeeping genes from tissue-specific genes and that repetitive sequence environment distinguishes housekeeping genes from tissue-specific genes in every isochore. The distinct repetitive sequence environment, in combination with other previously published sequence properties of housekeeping genes, was used to develop a method of predicting housekeeping genes on the basis of DNA sequence alone. Using expression across tissue types as a measure of success, we demonstrate that repetitive sequence environment is by far the most important sequence feature identified to date for distinguishing housekeeping genes.

Sauria SINEs: Novel short interspersed retroposable elements that are widespread in reptile genomes

SINEs are short interspersed retrotransposable elements that invade new genomic sites. Their retrotransposition depends on reverse transcriptase and endonuclease activities encoded by partner LINEs (long interspersed elements). Recent genomic research has demonstrated that retroposons account for at least 40% of the human genome. Hitherto, more than 30 families of SINEs have been characterized in mammalian genomes, comprising approximately 4600 extant species; the distribution and extent of SINEs in reptilian genomes, however, are poorly documented. With more than 7400 species of lizards and snakes, Squamata constitutes the largest and most diverse group of living reptiles. We have discovered and characterized a novel SINE family, Sauria SINEs, whose members are widely distributed among genomes of lizards, snakes, and tuataras. Sauria SINEs comprise a 5' tRNA-related region, a tRNA-unrelated region, and a 3' tail region (containing short tandem repeats) derived from LINEs. We distinguished eight Sauria SINE subfamilies in genomes of four major squamate lineages and investigated their evolutionary relationships. Our data illustrate the overall efficacy of Sauria SINEs as novel retrotransposable markers for elucidation of squamate evolutionary history. We show that all Sauria SINEs share an identical 3' sequence with Bov-B LINEs and propose that they utilize the enzymatic machinery of Bov-B LINEs for their own retrotransposition. This finding, along with the ubiquity of Bov-B LINEs previously demonstrated in squamate genomes, suggests that these LINEs have been an active partner of Sauria SINEs since this SINE family was generated more than 200 million years ago.

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.

SINEs and LINEs: symbionts of eukaryotic genomes with a common tail

Many SINEs and LINEs have been characterized to date, and examples of the SINE and LINE pair that have the same 3' end sequence have also increased. We report the phylogenetic relationships of nearly all known LINEs from which SINEs are derived, including a new example of a SINE/LINE pair identified in the salmon genome. We also use several biological examples to discuss the impact and significance of SINEs and LINEs in the evolution of vertebrate genomes.

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.

From the margins of the genome: mobile elements shape primate evolution

As is the case with mammals in general, primate genomes are inundated with repetitive sequence. Although much of this repetitive content consists of "molecular fossils" inherited from early mammalian ancestors, a significant portion of this material comprises active mobile element lineages. Despite indications that these elements played a major role in shaping the architecture of the genome, there remain many unanswered questions surrounding the nature of the host-element relationship. Here we review advances in our understanding of the host-mobile element dynamic and its overall impact on primate evolution.

Diversity and clustered distribution of retrotransposable elements in the compact genome of the pufferfish Tetraodon nigroviridis

We report the characterization and chromosomal distribution of retroelements in the compact genome of the pufferfish Tetraodon nigroviridis. We have reconstructed partial/complete retroelement sequences, established their phylogenetic relationship to other known eukaryotic retrotransposons, and performed double-color FISH analyses to gain new insights into their patterns of chromosomal distribution. We could identify 43 different reverse transcriptase retrotransposons belonging to the three major known subclasses (14 non-LTR retrotransposons from seven clades, 25 LTR retrotransposons representing the five major known groups, and four Penelope-like elements), and well as two SINEs (non-autonomous retroelements). Such a diversity of retrotransposable elements, which seems to be relatively common in fish but not in mammals, is astonishing in such a compact genome. The total number of retroelements was approximately 3000, roughly representing only 2.6% of the genome of T. nigroviridis. This is much less than in other vertebrate genomes, reflecting the compact nature of the genome of this pufferfish. Major differences in copy number were observed between different clades, indicating differential success in invading and persisting in the genome. Some retroelements displayed evidence of recent activity. Finally, FISH analysis showed that retrotransposable elements preferentially accumulate in specific heterochromatic regions of the genome of T. nigroviridis, revealing a degree of genomic compartmentalization not observed in the human genome.

The changing tails of a novel short interspersed element in Aedes aegypti: genomic evidence for slippage retrotransposition and the relationship between 3' tandem repeats and the poly(dA) tail

A novel family of tRNA-related SINEs named gecko was discovered in the yellow fever mosquito, Aedes aegypti. Approximately 7200 copies of gecko were distributed in the A. aegypti genome with a significant bias toward A + T-rich regions. The 3' end of gecko is similar in sequence and identical in secondary structure to the 3' end of MosquI, a non-LTR retrotransposon in A. aegypti. Nine conserved substitutions and a deletion separate gecko into two groups. Group I includes all gecko that end with poly(dA) and a copy that ends with AGAT repeats. Group II comprises gecko elements that end with CCAA or CAAT repeats. Members within each group cannot be differentiated when the 3' repeats are excluded in phylogenetic and sequence analyses, suggesting that the alterations of 3' tails are recent. Imperfect poly(dA) tail was recorded in group I and partial replication of the 3' tandem repeats was frequently observed in group II. Genomic evidence underscores the importance of slippage retrotransposition in the alteration and expansion of the tandem repeat during the evolution of gecko sequences, although we do not rule out postinsertion mechanisms that were previously invoked to explain the evolution of Alu-associated microsatellites. We propose that the 3' tandem repeats and the poly(dA) tail may be generated by similar mechanisms during retrotransposition of both SINEs and non-LTR retrotransposons and thus the distinction between poly(dA) retrotransposons such as L1 and non-poly(dA) retrotransposons such as I factor may not be informative.

BoS: A Large and Diverse Family of Short Interspersed Elements (SINEs) in Brassica oleracea

Short interspersed elements (SINEs) are nonautonomous non-LTR retrotransposons that populate eukaryotic genomes. Numerous SINE families have been identified in animals, whereas only a few have been described in plants. Here we describe a new family of SINEs, named BoS, that is widespread in Brassicaceae and present at approximately 2000 copies in Brassica oleracea. In addition to sharing a modular structure and target site preference with previously described SINEs, BoS elements have several unusual features. First, the head regions of BoS RNAs can adopt a distinct hairpin-like secondary structure. Second, with 15 distinct subfamilies, BoS represents one of the most diverse SINE families described to date. Third, several of the subfamilies have a mosaic structure that has arisen through the exchange of sequences between existing subfamilies, possibly during retrotransposition. Analysis of BoS subfamilies indicate that they were active during various time periods through the evolution of Brassicaceae and that active elements may still reside in some Brassica species. As such, BoS elements may be a valuable tool as phylogenetic makers for resolving outstanding issues in the evolution of species in the Brassicaceae family.

Isolation and Characterization of Active LINE and SINEs from the Eel

Long interspersed elements (LINEs) and short interspersed elements (SINEs) are retrotransposons. These elements can mobilize by the "copy-and-paste" mechanism, in which their own RNA is reverse-transcribed into complementary DNA (cDNA). LINEs and SINEs not only are components of eukaryotic genomes but also drivers of genomic evolution. Thus, studies of the amplification mechanism of LINEs and SINEs are important for understanding eukaryotic genome evolution. Here we report the characterization of one LINE family (UnaL2) and two SINE families (UnaSINE1 and UnaSINE2) from the eel (Anguilla japonica) genome. UnaL2 is approximately 3.6 kilobases (kb) and encodes only one open reading frame (ORF). UnaL2 belongs to the stringent type--thought to be a major group of LINEs--and can mobilize in HeLa cells. We also show that UnaL2 and the two UnaSINEs have similar 3' tails, and that both UnaSINE1 and UnaSINE2 can be mobilized by UnaL2 in HeLa cells. These elements are thus useful for delineating the amplification mechanism of stringent type LINEs as well as that of SINEs.

The Evolution of SINEs and LINEs in the genus Chironomus (Diptera)

Genomic DNA amplification from 51 species of the family Chironomidae shows that most contain relatives of NLRCth1 LINE and CTRT1 SINE retrotransposons first found in Chironomus thummi. More than 300 cloned PCR products were sequenced. The amplified region of the reverse transcriptase gene in the LINEs is intact and highly conserved, suggesting active elements. The SINEs are less conserved, consistent with minimal/no selection after transposition. A mitochondrial gene phylogeny resolves the Chironomus genus into six lineages (Guryev et al. 2001). LINE and SINE phylogenies resolve five of these lineages, indicating their monophyletic origin and vertical inheritance. However, both the LINE and the SINE tree topologies differ from the species phylogeny, resolving the elements into "clusters I-IV" and "cluster V" families. The data suggest a descent of all LINE and SINE subfamilies from two major families. Based on the species phylogeny, a few LINEs and a larger number of SINEs are cladisitically misplaced. Most misbranch with LINEs or SINEs from species with the same families of elements. From sequence comparisons, cladistically misplaced LINEs and several misplaced SINEs arose by convergent base substitutions. More diverged SINEs result from early transposition and some are derived from multiple source SINEs in the same species. SINEs from two species (C. dorsalis, C. pallidivittatus), expected to belong to the clusters I-IV family, branch instead with cluster V family SINEs; apparently both families predate separation of cluster V from clusters I-IV species. Correlation of the distribution of active SINEs and LINEs, as well as similar 3' sequence motifs in CTRT1 and NLRCth1, suggests coevolving retrotransposon pairs in which CTRT1 transposition depends on enzymes active during NLRCth1 LINE mobility.

Solution structure of an RNA stem-loop derived from the 3' conserved region of eel LINE UnaL2

The eel long interspersed element (LINE) UnaL2 and its partner short interspersed element (SINE) share a conserved 3' tail containing a stem-loop that is critical for their retrotransposition. Presumably, the first step of retrotransposition is the recognition of their 3' tails by UnaL2-encoded reverse transcriptase. The solution structure of a 17-nucleotide RNA derived from the 3' tail of UnaL2 was determined by NMR. The GGAUA loop forms a specific structure in which the uridine is exposed to solvent with the third and fifth adenosines stacked. A sharp turn in the phosphodiester backbone occurs between the second guanosine and third adenosine. When the uridine is mutated (but not deleted), all mutants form the loop structure, indicating that the loop structure requires an exposed fourth residue. The retrotransposition assay in HeLa cells revealed that retrotransposition requires the second guanosine, although any nucleoside functions at the fourth position, suggesting that UnaL2 reverse transcriptase specifically recognizes the 5' side of the GGANA loop.

A new pair of CR1-like LINE and tRNA-derived SINE elements in Podarcis sicula genome

We have identified and characterized a new pair of LINE and SINE elements, called Lucy-1 CR1-like LINE and P.s.1/SINE, respectively, in Podarcis sicula genome. The 3'-tail region in the 3' untranslated region (UTR) of Lucy-1 element is almost identical to the of P.s.1/SINE element. This identity suggests that the P.s.1/SINE element, during evolution, has gained the 3'-end sequence of the Lucy-1 element and has exclusively recruited the enzymatic machinery of its partner CR1 LINE for retroposition. Moreover, the complex molecular organization around Lucy-1 insertion site is discussed and we found that Lucy-1 insertion is associated with the calcium binding transporter gene. Our results confirm that the retrotransposons can be an additional source of genomic diversification and the evolution of the retrotransposable elements can be a vector force shaping genomes by reassorting DNA domains thus forming a new DNA arrangement.

Mosquito transposable elements

The completion of the genome assembly for the African malaria mosquito, Anopheles gambiae, and continuing genomic efforts for the yellow fever mosquito, Aedes aegypti, have allowed the use of bioinformatics tools to identify and characterize a diverse array of transposable elements (TEs) in these and other mosquito genomes. An overview of the types and number of both RNA-mediated and DNA-mediated TEs that are found in mosquito genomes is presented. A number of novel and interesting TEs from these species are discussed in more detail. These findings have significant implications for our understanding of mosquito genome evolution and for future modifications of natural mosquito populations through the use of TE-mediated genetic transformation.

Sequencing and expression of the second allele of the interleukin-1beta1 gene in rainbow trout (Oncorhynchus mykiss): identification of a novel SINE in the third intron

A lambda clone containing a rainbow trout IL-1beta1 gene was isolated by a PCR screening strategy from a genomic library cloned in lambda GEM-11, and an EcoRI fragment from this clone was fully sequenced, and contained 1680 bp 5'-flanking sequence, the whole IL-1beta1 gene open reading frame, and the 3'-flanking region with two potential poly A signals and poly A sites. This clone encoded a protein that shared 99.8% identity to the previously published trout IL-1beta1 cDNA sequence, with only three base substitutions. The main difference was that this clone had an additional complete HpaI SINE insertion in the 3rd intron making intron III 211 bp larger (834 bp via 623 bp). Thus this sequence was designated as allele B (Big intron III) of IL-1beta1 and the previously reported sequence as allele S (Short intron III). Three lines of evidence (allele specific PCR, cloning and sequencing, and direct sequencing of PCR products) revealed that allele B was constitutively expressed and could respond to stimulation with lipopolysaccharide or trout recombinant IL-1beta. Searching of the GenBank database with the HpaI SINE sequence resulted in three additional HpaI loci being identified in rainbow trout. Another SINE retroposition was also identified in the same intron of both alleles of IL-1beta1 by comparison with the trout IL-1beta2 gene. This novel SINE sequence, sharing high homology with the HpaI SINE at the 3'-end region, is present in EST databases of several species including human, mouse and fish. The consensus of this novel SINE shares 57 to 61% identities to tRNA-Leu from different species. Another older retroposition event in the same intron of IL-1beta1 has also been hypothesised, recognised as six adenines, that may function as a RNA polIII terminator. A model for the IL-1beta1 allele formation is proposed. Following the earliest retroposition into one of the two IL-1beta genes that resulted from a genome duplication in salmonids, the proper environment for successive PV SINE retroposition was created. A recent retroposition of the HpaI SINE in IL-1beta1 resulted in the formation of the two alleles of IL-1beta1. Examination of the SINEs insertion and their host gene microenvironments revealed that the SINE retroposition does not appear random, both in the site selection and the direction of insertion. The mechanism governing this outcome is discussed.

Gene discovery in the Entamoeba invadens genome

Entamoeba invadens, a parasite of reptiles, is a model for the study of encystation by the human enteric pathogen Entamoeba histolytica, because E. invadens form cysts in axenic culture. With approximately 0.5-fold sequence coverage of the genome, we were able to get insights into E. invadens gene and genome features. Overall, the E. invadens genome displays many of the features that are emerging from ongoing genome sequencing efforts in E. histolytica. At the nucleotide level the E. invadens genome has on average 60% sequence identity with that of E. histolytica. The presence of introns in E. invadens was predicted with similar consensus (GTTTGT em leader A/TAG) sequences to those identified in E. histolytica and Entamoeba dispar. Sequences highly repeated in the genome of E. histolytica (rRNAs, tRNAs, CXXC-rich proteins, and Leu-rich repeat proteins) were found to be highly repeated in the E. invadens genome. Numerous proteins homologous to those implicated in amoebic virulence, (Gal/GalNAc lectins, amoebapores, and cysteine proteinases) and drug resistance (p-glycoproteins) were identified. Homologs of proteins involved in cell cycle, vesicular trafficking and signal transduction were identified, which may be involved in en/excystation and cell growth of E. invadens. Finally, multiple copies of a number of E. invadens genes coding for predicted enzymes involved in core metabolism and the targets of anti-amoebic drugs were identified.

A patchwork interspersed sequence is present in a high copy number in the sheep genome

We have isolated a new interspersed sequence present in a high copy number in the ovine genome. This patchwork sequence, named 3.79 AS1, is part of a larger element encompassing similarities to constant region of reverse transcriptase and to art2 shared with the Bovine Dimer Driven Family (BDDF). The 3.79 AS1 sequence includes homologies to amplification promoting sequences (APS), to a potential origin of bidirectional DNA replication (OBR), to the Alu core sequence motif GGAGGC required for RNA polymerase III promoter function and to the ATGGCTGCCAT sequence that has been shown to be able to induce amplification-dependent transformation in murine cells. Fluorescent in situ hybridization experiments using probes derived from both ends of the 3.79 AS1 sequence showed a widespread signal over all sheep chromosomes, except the Y chromosome. We propose that the structural features of the 3.79 AS1 patchwork sequence, that is likely to be a subfamily of Bov B LINE that invaded the Artiodactyl genome prior to the separation of the Bovidae species, facilitated its massive amplification and dispersion in the ovine genome.

LINEs mobilize SINEs in the eel through a shared 3' sequence

We characterized members of the LINE (UnaL2) and SINE (UnaSINE1) families from the eel genome and found that these LINE/SINE partners share similar 3' tails. A retrotransposition assay in HeLa cells demonstrated that the 3' conserved tail of UnaL2 is necessary for its retrotransposition. This 3' tail is recognized in trans by the UnaL2 reverse transcriptase at a surprisingly high rate, and that of UnaSINE1 can also be recognized, thus providing experimental evidence that a SINE can be mobilized by the retrotransposition machinery of a partner LINE. We also demonstrated that short repeats at the 3' end of UnaL2 are required for retrotransposition suggesting that UnaL2 retrotransposes in a manner reminiscent of the reverse transcriptase activity of telomerases.

Mammalian retroelements

The eukaryotic genome has undergone a series of epidemics of amplification of mobile elements that have resulted in most eukaryotic genomes containing much more of this 'junk' DNA than actual coding DNA. The majority of these elements utilize an RNA intermediate and are termed retroelements. Most of these retroelements appear to amplify in evolutionary waves that insert in the genome and then gradually diverge. In humans, almost half of the genome is recognizably derived from retroelements, with the two elements that are currently actively amplifying, L1 and Alu, making up about 25% of the genome and contributing extensively to disease. The mechanisms of this amplification process are beginning to be understood, although there are still more questions than answers. Insertion of new retroelements may directly damage the genome, and the presence of multiple copies of these elements throughout the genome has longer-term influences on recombination events in the genome and more subtle influences on gene expression.

Mosaic structure and retropositional dynamics during evolution of subfamilies of short interspersed elements in African cichlids

The African cichlid (AFC) family of short interspersed elements (SINEs) is found in the genomes of cichlid fish. The alignment of the sequences of 70 members of this family, isolated from such fish in Africa, revealed the presence of correlated changes in specific nucleotides (diagnostic nucleotides) that allowed us to categorize the various members into six subfamilies, which were designated Af1 through Af6. Dividing the SINE consensus sequence into a 5'-head and 3'-tail region, these subfamilies were defined by various combinations of four types of head region (A-D) and three types of tail region [X, Y, and (YX)], with each region of each type including unique diagnostic nucleotides. The observed structures of the subfamilies Af1 through Af6 were AX, AY, CY, A(YX), BY, and DX, respectively. The formation of such structures might have involved the shuffling of head or tail regions among preexisting and existing (or both) subfamilies of the AFC family (and, probably, even another SINE family or a pseudogene for a tRNA in the case of the Af6 subfamily) by recombination at the so-called core region during the course of evolution. By plotting the timing of the retroposition of individual members of each subfamily on a phylogenetic tree of AFCs, we found that the Af3 and Af6 subfamilies became active only recently in the evolutionary history of these fish. The integrity of the 3'-tails of SINEs, which are, apparently, recognized by reverse transcriptase, has been reported to be indispensable for retention of retropositional activity. Therefore, we postulate that recombination might have been involved in the apparent recent activation of the retroposition of the Af3 and Af6 subfamilies via introduction of active tails (types Y and X, respectively) into potential ancestral sequences that might have had inactive tails. If this hypothesis is correct, shuffling of tail regions among subfamilies by recombination at the core region might have played a role in the recycling of dead copies of AFC SINEs.

Alu repeats and human genomic diversity

During the past 65 million years, Alu elements have propagated to more than one million copies in primate genomes, which has resulted in the generation of a series of Alu subfamilies of different ages. Alu elements affect the genome in several ways, causing insertion mutations, recombination between elements, gene conversion and alterations in gene expression. Alu-insertion polymorphisms are a boon for the study of human population genetics and primate comparative genomics because they are neutral genetic markers of identical descent with known ancestral states.

Pstl repeat: a family of short interspersed nucleotide element (SINE)-like sequences in the genomes of cattle, goat, and buffalo

The PstI family of elements are short, highly repetitive DNA sequences interspersed throughout the genome of the Bovidae. We have cloned and sequenced some members of the PstI family from cattle, goat, and buffalo. These elements are approximately 500 bp, have a copy number of 2 x 10(5) - 4 x 10(5), and comprise about 4% of the haploid genome. Studies of nucleotide sequence homology indicate that the buffalo and goat PstI repeats (type II) are similar types of short interspersed nucleotide element (SINE) sequences, but the cattle PstI repeat (type I) is considerably more divergent. Additionally, the goat PstI sequence showed significant sequence homology with bovine serine tRNA, and is therefore likely derived from serine tRNA. Interestingly, Southern hybridization suggests that both types of SINEs (I and II) are present in all the species of Bovidae. Dendrogram analysis indicates that cattle PstI SINE is similar to bovine Alu-like SINEs. Goat and buffalo SINEs formed a separate cluster, suggesting that these two types of SINEs evolved separately in the genome of the Bovidae.

Evolutionary dynamics and evolutionary history in the RTE clade of non-LTR retrotransposons

This study examined the evolutionary dynamics of Bov-B LINEs in vertebrates and the evolution of the RTE clade of non-LTR retrotransposons. The first full-length reptilian Bov-B LINE element is described; it is 3.2 kb in length, with a structural organization typical of the RTE clade of non-LTR retrotransposons. The long-term evolution of Bov-B LINEs was studied in 10 species of Squamata by analysis of a PCR-amplified 1.8-kb fragment encoding part of apurinic/apyrimidinic endonuclease, the intervening domain, and the palm/fingers subdomain of reverse transcriptase. A very high level of conservation in Squamata Bov-B long interspersed nuclear elements has been found, reaching 86% identity in the nearly 600 amino acids of ORF2. The same level of conservation exists between the ancestral snake lineage and Ruminantia. Such a high level is exceptional when compared with the level of conservation observed in nuclear and mitochondrial proteins and in other transposable elements. The RTE clade has been found to be much more widely distributed than previously thought, and novel representatives have been discovered in plants, brown algae, annelids, crustaceans, mollusks, echinoderms, and teleost fishes. Evolutionary relationships in the RTE clade were deduced at the amino acid level from three separate regions of ORF2. By using different independent methods, including the divergence-versus-age analysis, several examples of horizontal transfer in the RTE clade were recognized, with important implications for the existence of HT in non-LTR retrotransposons.

Gene organization of bovine BCNT that contains a portion corresponding to an endonuclease domain derived from an RTE-1 (Bov-B LINE), non-LTR retrotransposable element: duplication of an intramolecular repeat unit downstream of the truncated RTE-1

BCNT (a protein named after Bucentaur or craniofacial development protein 1) has a unique structure in Ruminantia. Bovine BCNT contains a region of the endonuclease domain derived from a truncated RTE-1 (previously called Bov-B LINE), a non-LTR retrotransposable repetitive element, and two repeat units (intramolecular repeat, IR) each with 40 amino acids in the C-terminal region. In contrast the human and mouse BCNT proteins contain one repeat unit and lack the RTE-1-derived portion. The 3' UTR of bovine bcnt cDNA also contains an approximately 300-bp portion homologous to the 3'-part of RTE-1. We examined the bovine bcnt genomic DNA sequence to understand how the bovine bcnt gene has been organized. The sequence of 3' UTR homologous portion was found to more closely resemble the Art2 element than the bovine RTE-1. By PCR screening a bovine/hamster hybrid somatic cell panel, the bovine bcnt gene was mapped to chromosome 18, syntenic human chromosome 16q on which human BCNT is located. The bcnt genomic DNA sequence corresponding to the cDNA downstream of a RTE-1 derived portion reveals that each IR unit is flanked by both 5'-side and 3'-side introns and that 3'-UTR consists of one exon. The alignment of the above sequence with a bovine RTE-1 did not show any significant homology downstream of the endonuclease domain. On the other hand, the alignment of the intron sequences with each other revealed that the six sequential homologous segments ranging in size from 40 to 453 bp existed over a 1 kb long sequence between both the 5'- and 3'-side introns flanking each bovine IR unit. In addition, both the 174-bp of 5'-side intron and 80-bp of 3'-side intron neighboring each 120-bp IR exon are significantly homologous among the two bovine IRs, human IR and mouse IR. These results suggest that a truncated bovine RTE-1 was inserted into the intron upstream of an IR unit of an ancestor bcnt gene and that a duplication of a relatively long region that includes IR occurred in the bovine genome.

Maque, a family of extremely short interspersed repetitive elements: characterization, possible mechanism of transposition, and evolutionary implications

Database analysis revealed a novel family of very short interspersed repetitive elements named Maque in the African malaria mosquito, Anopheles gambiae. Past mobility of Maque was demonstrated by evidence of its insertion that resulted in a target duplication. Approximately 220 copies of Maque were present in the A. gambiae genome. Although only approximately 60 bp long, Maque has the appearance of a distinct transposition unit. Eleven of the 12 Maque elements found in the database were flanked by 9-14 bp direct repeats, indicating that their transposition was relatively recent. Sequence comparison and phylogenetic analyses suggest that there are at least two subgroups within the Maque family, suggesting that they may have been originated from more than one source. Five of the 12 Maque elements had at least one other repetitive element nearby. Three of the Maque elements were found near genes. However, Maque was not found in the coding regions of genes or in any of the expressed sequence tags (ESTs), which is consistent with its significantly biased distribution toward A + T rich regions. Several characteristics of Maque indicate that it is likely a non-autonomous retro-element. The evolutionary origin of Maque and the differences between Maque and other known retro-elements including short interspersed repetitive elements (SINEs) are discussed. A hypothesis is proposed in which short sequences containing just the reverse transcriptase recognition signal (RTRS) could potentially contribute to exon shuffling and the genesis of some primordial SINEs.

Insect globin gene polymorphisms: intronic minisatellites and a retroposon interrupting exon 1 of homologous globin genes in Chironomus (Diptera)

Exon 1 of globin gene ct-13RT in clone lambdagb2-1 from Chironomus thummi contains a 444nt SINE (CTRT1). Based on in situ hybridization to polytene salivary gland chromosomes, C. thummi (ct), C. piger (cp) and C. tentans (ctn) contain copies of CTRT1 at multiple chromosomal loci. Genomic PCR amplifications reveal interrupted (ct-13RT) and uninterrupted (ct-13) alleles of the globin gene in the German population of C. thummi maintained in our laboratory, and only uninterrupted alleles or their homologs in different populations of C. thummi, C. piger and C. tentans. PCR amplification did generate different length fragments from cp-13 gene homologs in natural and laboratory C. piger populations that were due to variation in the length of minisatellite expansions of the central introns of the genes rather than a CTRT1-like SINE. Within minisatellite arrays, aligned homologs were more similar than paralogs in a single population, indicating that a tandem cluster of these repeats predated separation of the C. piger populations. The ct-13 genes of several C. thummi populations lack the minisatellites, suggesting their origin in C. piger only after the thummi/piger split. CTRT1 transposition into a ct-13 allele is even more recent, occurring after separation of German and other European C. thummi populations. The nearly intact ct-13RT and comparison with its intact ct-13 allele support a very recent transposition of the CTRT1 SINE into one of at least two already diverged ct-13 globin gene alleles. PCR analysis of DNA from individual adults in C. thummi shows a 1:2:1 distribution of ct-13/ct-13:ct-13/ct-13RT:ct-13RT/ct-13RT genotypes, consistent with a neutral spread of the ct-13RT allele since transposition, and indicating that the hemoglobin encoded by ct-13 is not necessary for survival, at least in a laboratory population of C. thummi.