Molecular evolution of Bov-B LINEs in vertebrates

Genome Evolution and the Future of Phylogenomics of Non-Avian Reptiles

Non-avian reptiles comprise a large proportion of amniote vertebrate diversity, with squamate reptiles-lizards and snakes-recently overtaking birds as the most species-rich tetrapod radiation. Despite displaying an extraordinary diversity of phenotypic and genomic traits, genomic resources in non-avian reptiles have accumulated more slowly than they have in mammals and birds, the remaining amniotes. Here we review the remarkable natural history of non-avian reptiles, with a focus on the physical traits, genomic characteristics, and sequence compositional patterns that comprise key axes of variation across amniotes. We argue that the high evolutionary diversity of non-avian reptiles can fuel a new generation of whole-genome phylogenomic analyses. A survey of phylogenetic investigations in non-avian reptiles shows that sequence capture-based approaches are the most commonly used, with studies of markers known as ultraconserved elements (UCEs) especially well represented. However, many other types of markers exist and are increasingly being mined from genome assemblies in silico, including some with greater information potential than UCEs for certain investigations. We discuss the importance of high-quality genomic resources and methods for bioinformatically extracting a range of marker sets from genome assemblies. Finally, we encourage herpetologists working in genomics, genetics, evolutionary biology, and other fields to work collectively towards building genomic resources for non-avian reptiles, especially squamates, that rival those already in place for mammals and birds. Overall, the development of this cross-amniote phylogenomic tree of life will contribute to illuminate interesting dimensions of biodiversity across non-avian reptiles and broader amniotes.

Evolutionary impact of transposable elements on genomic diversity and lineage-specific innovation in vertebrates

Since their discovery, a growing body of evidence has emerged demonstrating that transposable elements are important drivers of species diversity. These mobile elements exhibit a great variety in structure, size and mechanisms of transposition, making them important putative actors in organism evolution. The vertebrates represent a highly diverse and successful lineage that has adapted to a wide range of different environments. These animals also possess a rich repertoire of transposable elements, with highly diverse content between lineages and even between species. Here, we review how transposable elements are driving genomic diversity and lineage-specific innovation within vertebrates. We discuss the large differences in TE content between different vertebrate groups and then go on to look at how they affect organisms at a variety of levels: from the structure of chromosomes to their involvement in the regulation of gene expression, as well as in the formation and evolution of non-coding RNAs and protein-coding genes. In the process of doing this, we highlight how transposable elements have been involved in the evolution of some of the key innovations observed within the vertebrate lineage, driving the group's diversity and success.

[STRUCTURE AND DISTRIBUTION OF THE RETROTRANSPOSON BOV-B LINE]

The classification of mobile elements was discussed. Special attention was devoted to the retroelement of the LINE group: retrotransposon Bov-B LINE. The history of its origin and distribution in the nature was considered. The results of the phenomenon of horizontal transition of the retrotransposon Bov-B LINE between evolutionally distant classes were discussed.

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.

The genome sequence of taurine cattle: a window to ruminant biology and evolution

To understand the biology and evolution of ruminants, the cattle genome was sequenced to about sevenfold coverage. The cattle genome contains a minimum of 22,000 genes, with a core set of 14,345 orthologs shared among seven mammalian species of which 1217 are absent or undetected in noneutherian (marsupial or monotreme) genomes. Cattle-specific evolutionary breakpoint regions in chromosomes have a higher density of segmental duplications, enrichment of repetitive elements, and species-specific variations in genes associated with lactation and immune responsiveness. Genes involved in metabolism are generally highly conserved, although five metabolic genes are deleted or extensively diverged from their human orthologs. The cattle genome sequence thus provides a resource for understanding mammalian evolution and accelerating livestock genetic improvement for milk and meat production.

Molecular evolution of toxin genes in Elapidae snakes

The venom of the sea krait, Laticauda semifasciata, consists primarily of two toxic proteins, phospholipase A(2) (PLA(2)) and a three-finger-structure toxin. We have cloned both toxic protein genes, including the upstream region. PLA(2) genes contain three types of inserted sequences: an AG-rich region, a chicken repeat 1-like long interspersed nucleotide element sequence and an intron II 3' side repeat sequence. The molecular divergence of L. semifasciata PLA(2) genes was defined on the basis of the inserted sequences and their sequence homology. The length of intron I in the three-finger-structure toxin genes differs from species to species. The alignment analysis of intron I of the three-finger-structure toxin genes revealed that the intron I sequence of the ancestral gene comprised ten genetic regions. A hypothetical evolutionary process for the three-finger-structure toxin genes has also been developed.

Transcribed Tc1-like transposons in salmonid fish

BACKGROUND

Mobile genetic elements comprise a substantial fraction of vertebrate genomes. These genes are considered to be deleterious, and in vertebrates they are usually inactive. High throughput sequencing of salmonid fish cDNA libraries has revealed a large number of transposons, which remain transcribed despite inactivation of translation. This article reports on the structure and potential role of these genes.

RESULTS

A search of EST showed the ratio of transcribed transposons in salmonid fish (i.e., 0.5% of all unique cDNA sequences) to be 2.4-32 times greater than in other vertebrate species, and 68% of these genes belonged to the Tc1-family of DNA transposons. A phylogenetic analysis of reading frames indicate repeated transposition of distantly related genes into the fish genome over protracted intervals of evolutionary time. Several copies of two new DNA transposons were cloned. These copies showed relatively little divergence (11.4% and 1.9%). The latter gene was transcribed at a high level in rainbow trout tissues, and was present in genomes of many phylogenetically remote fish species. A comparison of synonymous and non-synonymous divergence revealed remnants of divergent evolution in the younger gene, while the older gene evolved in a neutral mode. From a 1.2 MB fragment of genomic DNA, the salmonid genome contains approximately 10(5) Tc1-like sequences, the major fraction of which is not transcribed. Our microarray studies showed that transcription of rainbow trout transposons is activated by external stimuli, such as toxicity, stress and bacterial antigens. The expression profiles of Tc1-like transposons gave a strong correlation (r2 = 0.63-0.88) with a group of genes implicated in defense response, signal transduction and regulation of transcription.

CONCLUSION

Salmonid genomes contain a large quantity of transcribed mobile genetic elements. Divergent or neutral evolution within genomes and lateral transmission can account for the diversity and sustained persistence of Tc1-like transposons in lower vertebrates. A small part of transposons remain transcribed and their transcription is enhanced by responses to acute conditions.

The origin and evolution of mosquito APE retroposons

The detection of horizontal transfer is important to understanding the origin and spread of transposable elements and in assessing their impact on genetic diversity. The occurrence of the phenomenon is not in doubt for two of the three major groups of elements, but is disputed for retroposons, largely on the grounds of data paucity and overreliance on divergence estimates between host species. We present here the most wide-ranging retroposon data set assembled to date for a species group, the mosquitoes. The results provide no evidence for horizontal transfer events and show conclusively that four previously reported events, involving Juan-A, Juan-C, T1, and Q, did not occur. We propose that the origin of all known mosquito retroposons can be attributed to vertical inheritance and that retroposons have therefore been a persistent source of genetic diversity in mosquito genomes since the emergence of the taxon. Furthermore, the data confirm that the unprecedented levels of retroposon diversity previously reported in Anopheles gambiae extends to at least seven other species representing five genera and all three mosquito subfamilies. Most notably, this included the L1 elements, which are not known in other insects. A number of novel well-defined monophyletic groups were also identified, particularly, JM2 and JM3 within the Jockey clade, which included sequences from seven and five mosquito species, respectively. As JM3 does not contain an Anopheles element, this represents a good example of stochastic loss and the best out of at least four found in this study. This exceptionally diverse data set when compared with the wealth of data available for the many unrelated species with which mosquitoes have intimate contact through blood feeding ought to be fertile ground for the discovery of horizontal transfer events. The absence of positive results therefore supports the view that retroposon horizontal transfer does not occur or is far more exceptional than for other types of transposable elements.

A processed pseudogene contributes to apparent mule deer prion gene heterogeneity

Pathogenesis and transmission of the prion disorders (transmissible spongiform encephalopathies, TSEs) are mediated by a modified isoform of the prion protein (PrP). Prion protein gene (PRNP) alleles associated with relative susceptibility to TSE have been identified in sheep, humans and possibly elk. Comparable data have not been derived for mule deer, a species susceptible to the TSE chronic wasting disease (CWD). Initial analysis of the open reading frame (ORF) in exon 3 of the mule deer PRNP gene revealed polymorphisms in all 145 samples analyzed, with 10 potential polymorphic sites. Because 144/145 (99.3%) of the samples were heterozygous for a coding change (N/S) at codon 138 (bp 412) and a non-coding polymorphism at bp 418, and individual deer with three or four different alleles were identified a possible gene duplication was indicated. Analysis of BAC clones containing mule deer PRNP genes revealed a full length functional gene and a processed pseudogene. The pseudogene was characteristic of previously described retroelements, in that it lacks introns and is flanked by repeat sequences. Three alleles of the functional gene were identified, with coding changes only at codons 20 (D/G) and 225 (S/F). Determination of PRNP functional gene alleles from 47 CWD-positive mule deer showed the predominant allele encoded 20D225S (frequency 0.85). When alleles were grouped by coding changes in the functional gene, four of the six possible peptide combinations were identified in infected deer. Three pseudogene alleles with coding changes in exon 3 were identified in the mule deer samples examined. Because the TSEs appear to be "protein only" disorders, the presence of an untranslated pseudogene is not expected to affect disease resistance. Therefore, selection of a genotyping method specific for the functional gene is critical for large-scale studies to identify the role of the PRNP gene in susceptibility to CWD in mule deer.

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.

Evolutionary dynamics in a novel L2 clade of non-LTR retrotransposons in Deuterostomia

The evolution of the novel L2 clade of non-long terminal repeat (LTR) retrotransposons and their evolutionary dynamics in Deuterostomia has been examined. The short-term evolution of long interspersed nuclear element 2s (LINE2s) has been studied in 18 reptilian species by analysis of a PCR amplified 0.7-kb fragment encoding the palm/fingers subdomain of reverse transcriptase (RT). Most of the reptilian LINE2s examined are inactive since they contain multiple stop codons, indels, or frameshift mutations that disrupt the RT. Analysis of reptilian LINE2s has shown a high degree of sequence divergence and an unexpectedly large number of deletions. The evolutionary dynamics of LINE2s in reptiles has been found to be complex. LINE2s are shown to form a novel clade of non-LTR retrotransposons that is well separated from the CR1 clade. This novel L2 clade is more widely distributed than previously thought, and new representatives have been discovered in echinoderms, insects, teleost fishes, Xenopus, Squamata, and marsupials. There is an apparent absence of LINE2s from different vertebrate classes, such as cartilaginous fishes, Archosauria (birds and crocodiles), and turtles. Whereas the LINE2s are present in echinoderms and teleost fishes in a conserved form, in most tetrapods only highly degenerated pseudogenes can be found. The predominance of inactive LINE2s in Tetrapoda indicates that, in the host genomes, only inactive copies are still present. The present data indicate that the vertical inactivation of LINE2s might have begun at the time of Tetrapoda origin, 400 MYA. The evolutionary dynamics of the L2 clade in Deuterostomia can be described as a gradual vertical inactivation in Tetrapoda, stochastic loss in Archosauria and turtles, and strict vertical transmission in echinoderms and teleost fishes.

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.