Evolutionary dynamics of transposable elements in the short-tailed opossum Monodelphis domestica

Evolution of Endogenous Retroviruses in the Subfamily of Caprinae

The interest in endogenous retroviruses (ERVs) has been fueled by their impact on the evolution of the host genome. In this study, we used multiple pipelines to conduct a de novo exploration and annotation of ERVs in 13 species of the Caprinae subfamily. Through analyses of sequence identity, structural organization, and phylogeny, we defined 28 ERV groups within Caprinae, including 19 gamma retrovirus groups and 9 beta retrovirus groups. Notably, we identified four recent and potentially active groups prevalent in the Caprinae genomes. Additionally, our investigation revealed that most long noncoding genes (lncRNA) and protein-coding genes (PC) contain ERV-derived sequences. Specifically, we observed that ERV-derived sequences were present in approximately 75% of protein-coding genes and 81% of lncRNA genes in sheep. Similarly, in goats, ERV-derived sequences were found in approximately 74% of protein-coding genes and 75% of lncRNA genes. Our findings lead to the conclusion that the majority of ERVs in the Caprinae genomes can be categorized as fossils, representing remnants of past retroviral infections that have become permanently integrated into the genomes. Nevertheless, the identification of the Cap_ERV_20, Cap_ERV_21, Cap_ERV_24, and Cap_ERV_25 groups indicates the presence of relatively recent and potentially active ERVs in these genomes. These particular groups may contribute to the ongoing evolution of the Caprinae genome. The identification of putatively active ERVs in the Caprinae genomes raises the possibility of harnessing them for future genetic marker development.

Differential Conservation and Loss of CR1 Retrotransposons in Squamates Reveals Lineage-Specific Genome Dynamics across Reptiles

Transposable elements (TEs) are repetitive DNA sequences which create mutations and generate genetic diversity across the tree of life. In amniotic vertebrates, TEs have been mainly studied in mammals and birds, whose genomes generally display low TE diversity. Squamates (Order Squamata; ~11,000 extant species of lizards and snakes) show as much variation in TE abundance and activity as they do in species and phenotypes. Despite this high TE activity, squamate genomes are remarkably uniform in size. We hypothesize that novel, lineage-specific dynamics have evolved over the course of squamate evolution to constrain genome size across the order. Thus, squamates may represent a prime model for investigations into TE diversity and evolution. To understand the interplay between TEs and host genomes, we analyzed the evolutionary history of the CR1 retrotransposon, a TE family found in most tetrapod genomes. We compared 113 squamate genomes to the genomes of turtles, crocodilians, and birds, and used ancestral state reconstruction to identify shifts in the rate of CR1 copy number evolution across reptiles. We analyzed the repeat landscapes of CR1 in squamate genomes and determined that shifts in the rate of CR1 copy number evolution are associated with lineage-specific variation in CR1 activity. We then used phylogenetic reconstruction of CR1 subfamilies across amniotes to reveal both recent and ancient CR1 subclades across the squamate tree of life. The patterns of CR1 evolution in squamates contrast other amniotes, suggesting key differences in how TEs interact with different host genomes and at different points across evolutionary history.

Expanding the evo-devo toolkit: generation of 3D mammary tissue from diverse mammals

The varying pathways of mammary gland development across species and evolutionary history are underexplored, largely due to a lack of model systems. Recent progress in organoid technology holds the promise of enabling in-depth studies of the developmental adaptations that have occurred throughout the evolution of different species, fostering beneficial phenotypes. The practical application of this technology for mammary glands has been mostly confined to rodents and humans. In the current study, we have successfully created next-generation 3D mammary gland organoids from eight eutherian mammals and the first branched organoid of a marsupial mammary gland. Using mammary organoids, we identified a role for ROCK protein in regulating branching morphogenesis, a role that manifests differently in organoids from different mammals. This finding demonstrates the utility of the 3D organoid model for understanding the evolution and adaptations of signaling pathways. These achievements highlight the potential for organoid models to expand our understanding of mammary gland biology and evolution, and their potential utility in studies of lactation or breast cancer.

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.

Paleozoic Protein Fossils Illuminate the Evolution of Vertebrate Genomes and Transposable Elements

Genomes hold a treasure trove of protein fossils: fragments of formerly protein-coding DNA, which mainly come from transposable elements (TEs) or host genes. These fossils reveal ancient evolution of TEs and genomes, and many fossils have been exapted to perform diverse functions important for the host's fitness. However, old and highly-degraded fossils are hard to identify, standard methods (e.g. BLAST) are not optimized for this task, and few Paleozoic protein fossils have been found. Here, a recently optimized method is used to find protein fossils in vertebrate genomes. It finds Paleozoic fossils predating the amphibian/amniote divergence from most major TE categories, including virus-related Polinton and Gypsy elements. It finds 10 fossils in the human genome (8 from TEs and 2 from host genes) that predate the last common ancestor of all jawed vertebrates, probably from the Ordovician period. It also finds types of transposon and retrotransposon not found in human before. These fossils have extreme sequence conservation, indicating exaptation: some have evidence of gene-regulatory function, and they tend to lienearest to developmental genes. Some ancient fossils suggest "genome tectonics", where two fragments of one TE have drifted apart by up to megabases, possibly explaining gene deserts and large introns. This paints a picture of great TE diversity in our aquatic ancestors, with patchy TE inheritance by later vertebrates, producing new genes and regulatory elements on the way. Host-gene fossils too have contributed anciently-conserved DNA segments. This paves the way to further studies of ancient protein fossils.

Constitutive Heterochromatin in Eukaryotic Genomes: A Mine of Transposable Elements

Transposable elements (TEs) are abundant components of constitutive heterochromatin of the most diverse evolutionarily distant organisms. TEs enrichment in constitutive heterochromatin was originally described in the model organism Drosophila melanogaster, but it is now considered as a general feature of this peculiar portion of the genomes. The phenomenon of TE enrichment in constitutive heterochromatin has been proposed to be the consequence of a progressive accumulation of transposable elements caused by both reduced recombination and lack of functional genes in constitutive heterochromatin. However, this view does not take into account classical genetics studies and most recent evidence derived by genomic analyses of heterochromatin in Drosophila and other species. In particular, the lack of functional genes does not seem to be any more a general feature of heterochromatin. Sequencing and annotation of Drosophila melanogaster constitutive heterochromatin have shown that this peculiar genomic compartment contains hundreds of transcriptionally active genes, generally larger in size than that of euchromatic ones. Together, these genes occupy a significant fraction of the genomic territory of heterochromatin. Moreover, transposable elements have been suggested to drive the formation of heterochromatin by recruiting HP1 and repressive chromatin marks. In addition, there are several pieces of evidence that transposable elements accumulation in the heterochromatin might be important for centromere and telomere structure. Thus, there may be more complexity to the relationship between transposable elements and constitutive heterochromatin, in that different forces could drive the dynamic of this phenomenon. Among those forces, preferential transposition may be an important factor. In this article, we present an overview of experimental findings showing cases of transposon enrichment into the heterochromatin and their positive evolutionary interactions with an impact to host genomes.

A reference genome for the critically endangered woylie, Bettongia penicillata ogilbyi

Biodiversity is declining globally, and Australia has one of the worst extinction records for mammals. The development of sequencing technologies means that genomic approaches are now available as important tools for wildlife conservation and management. Despite this, genome sequences are available for only 5% of threatened Australian species. Here we report the first reference genome for the woylie (Bettongia penicillata ogilbyi), a critically endangered marsupial from Western Australia, and the first genome within the Potoroidae family. The woylie reference genome was generated using Pacific Biosciences HiFi long-reads, resulting in a 3.39 Gbp assembly with a scaffold N50 of 6.49 Mbp and 86.5% complete mammalian BUSCOs. Assembly of a global transcriptome from pouch skin, tongue, heart and blood RNA-seq reads was used to guide annotation with Fgenesh++, resulting in the annotation of 24,655 genes. The woylie reference genome is a valuable resource for conservation, management and investigations into disease-induced decline of this critically endangered marsupial.

Genome-wide analyses of the relict gull (Larus relictus): insights and evolutionary implications

BACKGROUND

The relict gull (Larus relictus), was classified as vulnerable on the IUCN Red List and is a first-class national protected bird in China. Genomic resources for L. relictus are lacking, which limits the study of its evolution and its conservation.

RESULTS

In this study, based on the Illumina and PacBio sequencing platforms, we successfully assembled the genome of L. relictus, one of the few known reference genomes in genus Larus. The size of the final assembled genome was 1.21 Gb, with a contig N50 of 8.11 Mb. A total of 18,454 genes were predicted from the assembly results, with 16,967 (91.94%) of these genes annotated. The genome contained 92.52 Mb of repeat sequence, accounting for 7.63% of the assembly. A phylogenetic tree was constructed using 4902 single-copy orthologous genes, which showed L. relictus had closest relative of L. smithsonianus, with divergence time of 14.7 Mya estimated between of them. PSMC analyses indicated that L. relictus had been undergoing a long-term population decline during 0.01-0.1 Mya with a small effective population size fom 8800 to 2200 individuals.

CONCLUSIONS

This genome will be a valuable genomic resource for a range of genomic and conservation studies of L. relictus and will help to establish a foundation for further studies investigating whether the breeding population is a complex population. As the species is threatened by habitat loss and fragmentation, actions to protect L. relictus are suggested to alleviate the fragmentation of breeding populations.

Genetic Diversity and Population Structures in Chinese Miniature Pigs Revealed by SINE Retrotransposon Insertion Polymorphisms, a New Type of Genetic Markers

Simple Summary

Our previous studies suggested that the short interspersed nuclear element (SINE) retrotransposon insertion polymorphisms (RIPs), as a new type of molecular marker developed very recently, are ideal molecular markers and have the potential to be used for population genetic analysis and molecular breeding in pigs and possibly it can be extended to other livestock animals as well. However, no report is available for the application of SINE RIPs in population genetic analysis in livestock, including pigs. Here, we evaluated 30 SINE RIPs in several indigenous Chinese miniature pig breeds, including three subpopulations of Bama pigs (BM-cov, BM-clo, and BM-inb). BM-cov is a subpopulation conserved in the national conservation farm, and BM-clo is a closed population maintained over 30 years with only 2 boars and 14 sows imported from its original area, while BM-inb herd is an 18 generation continuous inbreeding line based on the BM-clo population. To our knowledge, it is the first time to report the genetic diversity, breed differentiation, and population structures for these populations by using SINE RIPs, and which suggests the feasibility of SINE RIPs in pig genetic analysis.

Abstract

RIPs have been developed as effective genetic markers and popularly applied for genetic analysis in plants, but few reports are available for domestic animals. Here, we established 30 new molecular markers based on the SINE RIPs, and applied them for population genetic analysis in seven Chinese miniature pigs. The data revealed that the closed herd (BM-clo), inbreeding herd (BM-inb) of Bama miniature pigs were distinctly different from the BM-cov herds in the conservation farm, and other miniature pigs (Wuzhishan, Congjiang Xiang, Tibetan, and Mingguang small ear). These later five miniature pig breeds can further be classified into two clades based on a phylogenetic tree: one included BM-cov and Wuzhishan, the other included Congjiang Xiang, Tibetan, and Mingguang small ear, which was well-supported by structure analysis. The polymorphic information contents estimated by using SINE RIPs are lower than the predictions based on microsatellites. Overall, the genetic distances and breed-relationships between these populations revealed by 30 SINE RIPs generally agree with their evolutions and geographic distributions. We demonstrated the potential of SINE RIPs as new genetic markers for genetic monitoring and population structure analysis in pigs, which can even be extended to other livestock animals.

The whale shark genome reveals how genomic and physiological properties scale with body size

The endangered whale shark (Rhincodon typus) is the largest fish on Earth and a long-lived member of the ancient Elasmobranchii clade. To characterize the relationship between genome features and biological traits, we sequenced and assembled the genome of the whale shark and compared its genomic and physiological features to those of 83 animals and yeast. We examined the scaling relationships between body size, temperature, metabolic rates, and genomic features and found both general correlations across the animal kingdom and features specific to the whale shark genome. Among animals, increased lifespan is positively correlated to body size and metabolic rate. Several genomic traits also significantly correlated with body size, including intron and gene length. Our large-scale comparative genomic analysis uncovered general features of metazoan genome architecture: Guanine and cytosine (GC) content and codon adaptation index are negatively correlated, and neural connectivity genes are longer than average genes in most genomes. Focusing on the whale shark genome, we identified multiple features that significantly correlate with lifespan. Among these were very long gene length, due to introns being highly enriched in repetitive elements such as CR1-like long interspersed nuclear elements, and considerably longer neural genes of several types, including connectivity, activity, and neurodegeneration genes. The whale shark genome also has the second slowest evolutionary rate observed in vertebrates to date. Our comparative genomics approach uncovered multiple genetic features associated with body size, metabolic rate, and lifespan and showed that the whale shark is a promising model for studies of neural architecture and lifespan.

Hagfish genome reveals parallel evolution of 7SL RNA-derived SINEs

Background

Short interspersed elements (SINEs) are ubiquitous components of eukaryotic genomes. SINEs are composite transposable elements that are mobilized by non-long terminal repeat (non-LTR) retrotransposons, also called long interspersed elements (LINEs). The 3′ part of SINEs usually originated from that of counterpart non-LTR retrotransposons. The 5′ part of SINEs mostly originated from small RNA genes. SINE1 is a group of SINEs whose 5′ part originated from 7SL RNA, and is represented by primate Alu and murine B1. Well-defined SINE1 has been found only from Euarchontoglires, a group of mammals, in contrast to the wide distribution of SINE2, which has a tRNA-derived sequence, from animals to plants to protists. Both Alu and B1 are mobilized by L1-type non-LTR retrotransposons, which are the only lineage of autonomous non-LTR retrotransposons active in these mammalian lineages.

Results

Here a new lineage of SINE1 is characterized from the seashore hagfish Eptatretus burgeri genome. This SINE1 family, designated SINE1-1_EBu, is young, and is transposed by RTE-type non-LTR retrotransposon, not L1-type. Comparison with other SINE families from hagfish indicated the birth of SINE1-1_EBu through chimera formation of a 7SL RNA-derived sequence and an older tRNA-derived SINE family. It reveals parallel evolution of SINE1 in two vertebrate lineages with different autonomous non-LTR retrotransposon partners. The comparison between two SINE1 lineages supports that the RNA secondary structure of the Alu domain of 7SL RNA is required for the efficient retrotransposition.

Conclusions

The hagfish SINE1 is the first evident SINE1 family found outside of Euarchontoglires. Independent evolution of SINE1 with similar RNA secondary structure originated in 7SL RNA indicates the functional importance of 7SL RNA-derived sequence in the proliferation of SINEs.

Variation in base composition underlies functional and evolutionary divergence in non-LTR retrotransposons

Background

Non-LTR retrotransposons often exhibit base composition that is markedly different from the nucleotide content of their host’s gene. For instance, the mammalian L1 element is AT-rich with a strong A bias on the positive strand, which results in a reduced transcription. It is plausible that the A-richness of mammalian L1 is a self-regulatory mechanism reflecting a trade-off between transposition efficiency and the deleterious effect of L1 on its host. We examined if the A-richness of L1 is a general feature of non-LTR retrotransposons or if different clades of elements have evolved different nucleotide content. We also investigated if elements belonging to the same clade evolved towards different base composition in different genomes or if elements from different clades evolved towards similar base composition in the same genome.

Results

We found that non-LTR retrotransposons differ in base composition among clades within the same host but also that elements belonging to the same clade differ in base composition among hosts. We showed that nucleotide content remains constant within the same host over extended period of evolutionary time, despite mutational patterns that should drive nucleotide content away from the observed base composition.

Conclusions

Our results suggest that base composition is evolving under selection and may be reflective of the long-term co-evolution between non-LTR retrotransposons and their host. Finally, the coexistence of elements with drastically different base composition suggests that these elements may be using different strategies to persist and multiply in the genome of their host.

Transposable elements as genetic accelerators of evolution: contribution to genome size, gene regulatory network rewiring and morphological innovation

In the current era, as a growing number of genome sequence assemblies have been reported in animals, an in-depth analysis of transposable elements (TEs) is one of the most fundamental and essential studies for evolutionary genomics. Although TEs have, in general, been regarded as non-functional junk/selfish DNA, parasitic elements or harmful mutagens, studies have revealed that TEs have had a substantial and sometimes beneficial impact on host genomes in several ways. First, TEs are themselves diverse and thus provide lineage-specific characteristics to the genomes. Second, because TEs constitute a substantial fraction of animal genomes, they are a major contributing factor to evolutionary changes in genome size and composition. Third, host organisms have co-opted many repetitive sequences as genes, cis-regulatory elements and chromatin domain boundaries, which alter gene regulatory networks and in addition are partly involved in morphological evolution, as has been well documented in mammals. Here, I review the impact of TEs on various aspects of the genome, such as genome size and diversity in animals, as well as the evolution of gene networks and genome architecture in mammals. Given that a number of TE families probably remain to be discovered in many non-model organisms, unknown TEs may have contributed to gene networks in a much wider variety of animals than considered previously.

Development of retrotransposons insertion polymorphic markers and application in the genetic variation evaluation of Chinese Bama miniature pigs

Retrotransposons are genetic elements that can amplify themselves in a genome and are abundant in many eukaryotic organisms. In this study, we established some new short interspersed nuclear elements (SINE) and endogenous retroviruses (ERV) retrotransposons insertion polymorphism (RTIP) markers based on BLAT alignment tool strategy, and followed by PCR evaluation. We investigated the genetic variations among four subpopulations of Chinese Bama miniature pigs (BM), including BM in national conservation farm (BM-cov), BM inbreeding population (BM-inb) and BM closed Herd (BM-clo) in Guangxi University, and BM in the Experimental pig farm of Yangzhou University (BM-yzu). Genetic distance, polymorphism information content (PIC) and heterozygosity (He) of these markers in four of BM subpopulations were measured. Twelve SINE and twenty-eight ERV polymorphic molecular markers were identified in the four subpopulations. The BM-cov pigs represented the highest He and PIC, which indicated that BM-cov pigs maintain relatively highly genetic diversity. BM-inb pigs represented the lowest He and PIC indicating less variation and a high degree of inbreeding. Microsatellites polymorphism in four BM populations also well supported the results of these RTIP markers. In summary, retrotransposons insertion polymorphic markers could be a useful tool for population genetic variation analysis. Current SINE and ERV variation data may also provide a reference guide for the conservation and utilization of the BM miniature pig resource.

Convergence in LINE-1 nucleotide variations can benefit redundantly forming triplexes with lncRNA in mammalian X-chromosome inactivation

Background

Associations between X-inactive transcript (Xist)–long noncoding RNA (lncRNA) and chromatin are critical intermolecular interactions in the X-chromosome inactivation (XCI) process. Despite high-resolution analyses of the Xist RNA-binding sites, specific interaction sequences are yet to be identified. Based on elusive features of the association between Xist RNA and chromatin and the possible existence of multiple low-affinity binding sites in Xist RNA, we defined short motifs (≥5 nucleotides), termed as redundant UC/TC (r-UC/TC) or AG (r-AG) motifs, which may help in the mediation of triplex formation between the lncRNAs and duplex DNA.

Results

The study showed that r-UC motifs are densely dispersed throughout mouse and human Xist/XIST RNAs, whereas r-AG motifs are even more densely dispersed along opossum RNA-on-the-silent X (Rsx) RNA, and also along both full-length and truncated long interspersed nuclear elements (LINE-1s, L1s) of the three species. Predicted secondary structures of the lncRNAs showed that the length range of these sequence motifs available for forming triplexes was even shorter, mainly 5- to 9-nucleotides long. Quartz crystal microbalance (QCM) measurements and Monte Carlo (MC) simulations indicated that minimum-length motifs can reinforce the binding state by increasing the copy number of the motifs in the same RNA or DNA molecule. Further, r-AG motifs in L1s had a similar length-distribution pattern, regardless of the similarities in the length or sequence of L1s across the three species; this also applies to high-frequency mutations in r-AG motifs, which suggests convergence in L1 sequence variations.

Conclusions

Multiple short motifs in both RNA and duplex DNA molecules could be brought together to form triplexes with either Hoogsteen or reverse Hoogsteen hydrogen bonding, by which their associations are cooperatively enhanced. This novel triplex interaction could be involved in associations between lncRNA and chromatin in XCI, particularly at the sites of L1s. Potential binding of Xist/XIST/Rsx RNAs specifically at L1s is most likely preserved through the r-AG motifs conserved in mammalian L1s through convergence in L1 nucleotide variations and by maintaining a particular r-UC/r-AG motif ratio in each of these lncRNAs, irrespective of their poorly conserved sequences.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0173-4) contains supplementary material, which is available to authorized users.

Retrotransposons evolution and impact on lncRNA and protein coding genes in pigs

BACKGROUND

Retrotransposons are the major determinants of genome sizes and they have shaped both genes and genomes in mammalian organisms, but their overall activity, diversity, and evolution dynamics, particularly their impact on protein coding and lncRNA genes in pigs remain largely unknown.

RESULTS

In the present study, we performed de novo detection of retrotransposons in pigs by using multiple pipelines, four distinct families of pig-specific L1 s classified into 51 distinct subfamilies and representing four evolution models and three expansion waves of pig-specific SINEs represented by three distinct families were identified. ERVs were classified into 18 families and found two most "modern" subfamilies in the pig genome. The transposition activity of pig L1 was verified by experiment, the sense and antisense promoter activities of young L1 5'UTRs and ERV LTRs and expression profiles of young retrotransposons in multiple tissues and cell lines were also validated. Furthermore, retrotransposons had an extensive impact on lncRNA and protein coding genes at both the genomic and transcriptomic levels. Most protein coding and lncRNA (> 80%) genes contained retrotransposon insertions, and about half of protein coding genes (44.30%) and one-fourth (24.13%) of lncRNA genes contained the youngest retrotransposon insertions. Nearly half of protein coding genes (43.78%) could generate chimeric transcripts with retrotransposons. Significant distribution bias of retrotransposon composition, location, and orientation in lncRNA and protein coding genes, and their transcripts, were observed.

CONCLUSIONS

In the current study, we characterized the classification and evolution profile of retrotransposons in pigs, experimentally proved the transposition activity of the young pig L1 subfamily, characterized the sense and antisense expression profiles and promoter activities of young retrotransposons, and investigated their impact on lncRNA and protein coding genes by defining the mobilome landscapes at the genomic and transcriptomic levels. These findings help provide a better understanding of retrotransposon evolution in mammal and their impact on the genome and transcriptome.

Identification of transposable elements fused in the exonic region of the olive flounder genome

Transposable elements (TEs) are mobile genetic sequences that comprise a large portion of vertebrate genomes. The olive flounder (Paralichthys olivaceus) is a valuable marine resource in East Asia. The scope of most genomic studies on the olive flounder is limited to its immunology as their focus is the prevention of mass mortality of this species. Thus, for a broader understanding of the species, its genomic information is consistently in demand. Transcripts sequences were acquired from transcriptome analysis using gill tissues of 12 olive flounders. Distribution of TEs inserted in exonic region of the olive flounder genome was analyzed using RepeatMasker ( http://www.repeatmasker.org/ ). We found 1140 TEs in the exonic region of the genome and long interspersed nuclear elements (LINEs) and long terminal repeats (LTRs) insertions occurred with forward orientation preferences. Transposons belonging to the hAt, Gypsy, and LINE 1 (L1) subfamilies were the most abundant DNA transposons, LTRs, and long interspersed elements (LINEs), respectively. Finally, we carried out a gene ontology analysis to determine the function of TE-fused genes. These results provide some genomic information about TEs that is useful for future research on changes in properties and functions of genes by TEs in the olive flounder genome.

Human transposable elements in Repbase: genomic footprints from fish to humans

Repbase is a comprehensive database of eukaryotic transposable elements (TEs) and repeat sequences, containing over 1300 human repeat sequences. Recent analyses of these repeat sequences have accumulated evidences for their contribution to human evolution through becoming functional elements, such as protein-coding regions or binding sites of transcriptional regulators. However, resolving the origins of repeat sequences is a challenge, due to their age, divergence, and degradation. Ancient repeats have been continuously classified as TEs by finding similar TEs from other organisms. Here, the most comprehensive picture of human repeat sequences is presented. The human genome contains traces of 10 clades (L1, CR1, L2, Crack, RTE, RTEX, R4, Vingi, Tx1 and Penelope) of non-long terminal repeat (non-LTR) retrotransposons (long interspersed elements, LINEs), 3 types (SINE1/7SL, SINE2/tRNA, and SINE3/5S) of short interspersed elements (SINEs), 1 composite retrotransposon (SVA) family, 5 classes (ERV1, ERV2, ERV3, Gypsy and DIRS) of LTR retrotransposons, and 12 superfamilies (Crypton, Ginger1, Harbinger, hAT, Helitron, Kolobok, Mariner, Merlin, MuDR, P, piggyBac and Transib) of DNA transposons. These TE footprints demonstrate an evolutionary continuum of the human genome.

A Role for the Mutagenic DNA Self-Catalyzed Depurination Mechanism in the Evolution of 7SL-Derived RNAs

The Alu element, the most prevalent SINE (short interspersed element) in the human genome, is one of the many RNA-encoding genes that evolved from the 7SL RNA gene. During analysis of the evolution of 7SL-derived RNAs, two distinct evolutionary intermediates capable of self-catalyzed DNA depurination (SDP) were identified. These SDP sequences spontaneously create apurinic sites that can result in increased mutagenesis due to their error-prone repair. This DNA self-depurination mechanism has been shown both in vitro and in vivo to lead to substitution and short frameshift mutations at a frequency that far exceeds their occurrence due to random errors in DNA replication. In both evolutionary intermediates, the same self-depurination sequence overlaps motifs necessary for successful transcription and SRP9/14 (signal recognition particle) binding; hence, mutations in this region could disrupt RNA activity. Yet, the 7SL-derived RNAs that arose from the elements capable of SDP show significant diversity in this region, and every new sequence retains the transcription and SRP9/14-binding motifs, even as it has lost the SDP sequence. While some (but not all) of the mutagenesis can be alternatively attributed to CpG decay, the very fact that the self-depurinating sequences are selectively discarded in all cases suggests that this was evolutionarily motivated to prevent further destructive mutagenesis by the SDP mechanism.

Heterochromatin variation and LINE-1 distribution in Artibeus (Chiroptera, Phyllostomidae) from Central Amazon, Brazil

Species in the subgenus Artibeus Leach, 1821 are widely distributed in Brazil. Conserved karyotypes characterize the group with identical diploid number and chromosome morphology. Recent studies suggested that the heterochromatin distribution and accumulation patterns can vary among species. In order to assess whether variation can also occur within species, we have analyzed the chromosomal distribution of constitutive heterochromatin in A. planirostris (Spix, 1823) and A. lituratus (Olfers, 1818) from Central Amazon (North Brazil) and contrasted our findings with those reported for other localities in Brazil. In addition, Ag-NOR staining and FISH with 18S rDNA, telomeric, and LINE-1 probes were performed to assess the potential role that these different repetitive markers had in shaping the current architecture of heterochromatic regions. Both species presented interindividual variation of constitutive heterochromatin. In addition, in A. planirostris the centromeres of most chromosomes are enriched with LINE-1, colocated with pericentromeric heterochromatin blocks. Overall, our data indicate that amplification and differential distribution of the investigated repetitive DNAs might have played a significant role in shaping the chromosome architecture of the subgenus Artibeus.

Evolution and Diversity of Transposable Elements in Vertebrate Genomes

Transposable elements (TEs) are selfish genetic elements that mobilize in genomes via transposition or retrotransposition and often make up large fractions of vertebrate genomes. Here, we review the current understanding of vertebrate TE diversity and evolution in the context of recent advances in genome sequencing and assembly techniques. TEs make up 4-60% of assembled vertebrate genomes, and deeply branching lineages such as ray-finned fishes and amphibians generally exhibit a higher TE diversity than the more recent radiations of birds and mammals. Furthermore, the list of taxa with exceptional TE landscapes is growing. We emphasize that the current bottleneck in genome analyses lies in the proper annotation of TEs and provide examples where superficial analyses led to misleading conclusions about genome evolution. Finally, recent advances in long-read sequencing will soon permit access to TE-rich genomic regions that previously resisted assembly including the gigantic, TE-rich genomes of salamanders and lungfishes.

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.

When Genomics Is Not Enough: Experimental Evidence for a Decrease in LINE-1 Activity During the Evolution of Australian Marsupials

The autonomous transposable element LINE-1 is a highly abundant element that makes up between 15% and 20% of therian mammal genomes. Since their origin before the divergence of marsupials and placental mammals, LINE-1 elements have contributed actively to the genome landscape. A previous in silico screen of the Tasmanian devil genome revealed a lack of functional coding LINE-1 sequences. In this study we present the results of an in vitro analysis from a partial LINE-1 reverse transcriptase coding sequence in five marsupial species. Our experimental screen supports the in silico findings of the genome-wide degradation of LINE-1 sequences in the Tasmanian devil, and identifies a high frequency of degraded LINE-1 sequences in other Australian marsupials. The comparison between the experimentally obtained LINE-1 sequences and reference genome assemblies suggests that conclusions from in silico analyses of retrotransposition activity can be influenced by incomplete genome assemblies from short reads.

[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.

Massive contribution of transposable elements to mammalian regulatory sequences

Barbara McClintock discovered the existence of transposable elements (TEs) in the late 1940s and initially proposed that they contributed to the gene regulatory program of higher organisms. This controversial idea gained acceptance only much later in the 1990s, when the first examples of TE-derived promoter sequences were uncovered. It is now known that half of the human genome is recognizably derived from TEs. It is thus important to understand the scope and nature of their contribution to gene regulation. Here, we provide a timeline of major discoveries in this area and discuss how transposons have revolutionized our understanding of mammalian genomes, with a special emphasis on the massive contribution of TEs to primate evolution. Our analysis of primate-specific functional elements supports a simple model for the rate at which new functional elements arise in unique and TE-derived DNA. Finally, we discuss some of the challenges and unresolved questions in the field, which need to be addressed in order to fully characterize the impact of TEs on gene regulation, evolution and disease processes.

Comparative and evolutionary insights into CD4 gene across mammalian and avian taxa

The present day genetic architecture of a species bears much significance to its closely related species which is due to species-specific differences, shaped by different evolutionary forces across time scale. With the availability of whole genome sequence of several closely related species, it is now possible to infer evolutionary patterns of genes and genomes in specific lineages. To this respect, CD4 gene, primarily responsible for defensive mechanism in human, is conserved across a few taxa, and thus, comparative genomic studies could be useful for better understanding of host-pathogen biology. Comparative and evolutionary analyses were performed in eleven taxa (10 mammalian and avian) with different statistical algorithms. Phylogenetic inferences revealed recent divergence of human and chimpanzee, and pig was found to be diverged from rest of the taxa significantly. Additionally, gene length, microsatellites, and secondary structures were observed across taxa. The genetic architecture of CD4 gene and its evolutionary history in different mammalian taxa provide crucial evidence in support of the fact that this gene might have been evolving at a similar rate to other human immune system genes. Future population-based study and structural modeling would unravel the differential ability to interact with HIV virus and influence immune system in humans.

The devil is in the details: Transposable element analysis of the Tasmanian devil genome

The third marsupial genome was sequenced from the Tasmanian devil (Sarcophilus harrisii), a species that currently is driven to extinction by a rare transmissible cancer. The transposable element (TE) landscape of the Tasmanian devil genome revealed that the main driver of retrotransposition the Long INterspersed Element 1 (LINE1) seem to have become inactivated during the past 12 million years. Strangely, the Short INterspersed Elements (SINE), that normally hijacks the LINE1 retrotransposition system, became inactive prior to LINE1 at around 30 million years ago. The SINE inactivation was in vitro verified in several species. Here I discuss that the apparent LINE1 inactivation might be caused by a genome assembly artifact. The repetitive fraction of any genome is highly complex to assemble and the observed problems are not unique to the Tasmanian devil genome.

Organization and Chromosomal Distribution of Histone Genes and Transposable Rex Elements in the Genome of Astyanax bockmanni (Teleostei, Characiformes)

An important feature of eukaryotic organisms is the number of different repetitive DNA sequences in their genome, a feature not observed in prokaryotes. These sequences are considered to be important components for understanding evolutionary mechanisms and the karyotypic differentiation processes. Thus, we aimed to physically map the histone genes and transposable elements of the Rex family in 6 fish populations of Astyanax bockmanni. FISH results using a histone H1 gene probe showed fluorescent clusters in 2 chromosome pairs in all 6 samples analyzed. In contrast, FISH with a histone H3 probe showed conspicuous blocks in 4 chromosomes in 5 of the 6 populations analyzed. The sixth population revealed 7 chromosomes marked with this probe. Probes for the transposable elements Rex1 and Rex6 showed small sites dispersed on most chromosomes of the 6 populations, and the Rex3 element is located in a big block concentrated in only 1 acrocentric chromosome of 2 populations. As for the other populations, a Rex3 probe showed large blocks in more than 1 chromosome. Fish from Alambari and Campo Novo Stream have Rex3 elements dispersed along most of the chromosomes. Additionally, the conspicuous signals of Rex1, Rex3, and Rex6 were identified in the acrocentric B microchromosome of A. bockmanni found only in individuals of the Alambari River. Thus, we believe that different mechanisms drive the spread of repetitive sequences among the populations analyzed, which appear to be organized differently in the genome of A. bockmanni. The presence of transposable elements in the B chromosome also suggests that these sequences could play a role in the origin and maintenance of the supernumerary element in the genome of this species.

Targeted Capture of Phylogenetically Informative Ves SINE Insertions in Genus Myotis

Identification of retrotransposon insertions in nonmodel taxa can be technically challenging and costly. This has inhibited progress in understanding retrotransposon insertion dynamics outside of a few well-studied species. To address this problem, we have extended a retrotransposon-based capture and sequence method (ME-Scan [mobile element scanning]) to identify insertions belonging to the Ves family of short interspersed elements (SINEs) across seven species of the bat genus Myotis. We identified between 120,000 and 143,000 SINE insertions in six taxa lacking a draft genome by comparing to the M. lucifugus reference genome. On average, each Ves insertion was sequenced to 129.6 × coverage. When mapped back to the M. lucifugus reference genome, all insertions were confidently assigned within a 10-bp window. Polymorphic Ves insertions were identified in each taxon based on their mapped locations. Using cross-species comparisons and the identified insertion positions, a presence–absence matrix was created for approximately 796,000 insertions. Dollo parsimony analysis of more than 85,000 phylogenetically informative insertions recovered strongly supported, monophyletic clades that correspond with the biogeography of each taxa. This phylogeny is similar to previously published mitochondrial phylogenies, with the exception of the placement of M. vivesi. These results support the utility of our variation on ME-Scan to identify polymorphic retrotransposon insertions in taxa without a reference genome and for large-scale retrotransposon-based phylogenetics.

Interspersed repeats in the horse (Equus caballus); spatial correlations highlight conserved chromosomal domains

The interspersed repeat content of mammalian genomes has been best characterized in human, mouse and cow. In this study, we carried out de novo identification of repeated elements in the equine genome and identified previously unknown elements present at low copy number. The equine genome contains typical eutherian mammal repeats, but also has a significant number of hybrid repeats in addition to clade-specific Long Interspersed Nuclear Elements (LINE). Equus caballus clade specific LINE 1 (L1) repeats can be classified into approximately five subfamilies, three of which have undergone significant expansion. There are 1115 full-length copies of these equine L1, but of the 103 presumptive active copies, 93 fall within a single subfamily, indicating a rapid recent expansion of this subfamily. We also analysed both interspersed and simple sequence repeats (SSR) genome-wide, finding that some repeat classes are spatially correlated with each other as well as with G+C content and gene density. Based on these spatial correlations, we have confirmed that recently-described ancestral vs. clade-specific genome territories can be defined by their repeat content. The clade-specific Short Interspersed Nuclear Element correlations were scattered over the genome and appear to have been extensively remodelled. In contrast, territories enriched for ancestral repeats tended to be contiguous domains. To determine if the latter territories were evolutionarily conserved, we compared these results with a similar analysis of the human genome, and observed similar ancestral repeat enriched domains. These results indicate that ancestral, evolutionarily conserved mammalian genome territories can be identified on the basis of repeat content alone. Interspersed repeats of different ages appear to be analogous to geologic strata, allowing identification of ancient vs. newly remodelled regions of mammalian genomes.

Disentangling the relationship of the Australian marsupial orders using retrotransposon and evolutionary network analyses

The ancestors to the Australian marsupials entered Australia around 60 (54-72) Ma from Antarctica, and radiated into the four living orders Peramelemorphia, Dasyuromorphia, Diprotodontia, and Notoryctemorphia. The relationship between the four Australian marsupial orders has been a long-standing question, because different phylogenetic studies have not been able to consistently reconstruct the same topology. Initial in silico analysis of the Tasmanian devil genome and experimental screening in the seven marsupial orders revealed 20 informative transposable element insertions for resolving the inter- and intraordinal relationships of Australian and South American orders. However, the retrotransposon insertions support three conflicting topologies regarding Peramelemorphia, Dasyuromorphia, and Notoryctemorphia, indicating that the split between the three orders may be best understood as a network. This finding is supported by a phylogenetic reanalysis of nuclear gene sequences, using a consensus network approach that allows depicting hidden phylogenetic conflict, otherwise lost when forcing the data into a bifurcating tree. The consensus network analysis agrees with the transposable element analysis in that all possible topologies regarding Peramelemorphia, Dasyuromorphia, and Notoryctemorphia in a rooted four-taxon topology are equally well supported. In addition, retrotransposon insertion data support the South American order Didelphimorphia being the sistergroup to all other living marsupial orders. The four Australian orders originated within 3 Myr at the Cretaceous-Paleogene boundary. The rapid divergences left conflicting phylogenetic information in the genome possibly generated by incomplete lineage sorting or introgressive hybridization, leaving the relationship among Australian marsupial orders unresolvable as a bifurcating process millions of years later.

U6 snRNA Pseudogenes: Markers of Retrotransposition Dynamics in Mammals

Transposable elements comprise more than 45% of the human genome and long interspersed nuclear element 1 (LINE-1 or L1) is the only autonomous mobile element remaining active. Since its identification, it has been proposed that L1 contributes to the mobilization and amplification of other cellular RNAs and more recently, experimental demonstrations of this function has been described for many transcripts such as Alu, a nonautonomous mobile element, cellular mRNAs, or small noncoding RNAs. Detailed examination of the mobilization of various cellular RNAs revealed distinct pathways by which they could be recruited during retrotransposition; template choice or template switching. Here, by analyzing genomic structures and retrotransposition signatures associated with small nuclear RNA (snRNA) sequences, we identified distinct recruiting steps during the L1 retrotransposition cycle for the formation of snRNA-processed pseudogenes. Interestingly, some of the identified recruiting steps take place in the nucleus. Moreover, after comparison to other vertebrate genomes, we established that snRNA amplification by template switching is common to many LINE families from several LINE clades. Finally, we suggest that U6 snRNA copies can serve as markers of L1 retrotransposition dynamics in mammalian genomes.

terMITEs: miniature inverted-repeat transposable elements (MITEs) in the termite genome (Blattodea: Termitoidae)

Transposable elements (TEs) are discrete DNA sequences which are able to replicate and jump into different genomic locations. Miniature inverted-repeats TEs (MITEs) are non-autonomous DNA elements whose origin is still poorly understood. Recently, some MITEs were found to contain core repeats that can be arranged in tandem arrays; in some instances, these arrays have even given rise to satellite DNAs in the (peri)centromeric region of the host chromosomes. I report the discovery and analysis of three new MITEs found in the genome of several termite species (hence the name terMITEs) in two different families. For two of the MITEs (terMITE1-Tc1/mariner superfamily; terMITE2-piggyBac superfamily), evidence of past mobility was retrieved. Moreover, these two MITEs contained core repeats, 16 bp and 114 bp long respectively, exhibiting copy number variation. In terMITE2, the tandem duplication appeared associated with element degeneration, in line with a recently proposed evolutionary model on MITEs and the origin of tandem arrays. Concerning their genomic distribution, terMITE1 and terMITE3 appeared more frequently inserted close to coding regions while terMITE2 was mostly associated with TEs. Although MITEs are commonly distributed in coding regions, terMITE2 distribution is in line with that of other insects' piggyBac-related elements and of other small TEs found in termite genomes. This has been explained through insertional preference rather than through selective processes. Data presented here add to the knowledge on the poorly exploited polyneopteran genomes and will provide an interesting framework in which to study TEs' evolution and host's life history traits.

Evolutionary histories of transposable elements in the genome of the largest living marsupial carnivore, the Tasmanian devil

The largest living carnivorous marsupial, the Tasmanian devil (Sarcophilus harrisii), is the sole survivor of a lineage originating about 12 Ma. We set out to investigate the spectrum of transposable elements found in the Tasmanian devil genome, the first high-coverage genome of an Australian marsupial. Marsupial genomes have been shown to have the highest amount of transposable elements among vertebrates. We analyzed the horizontally transmitted DNA transposons OC1 and hAT-1_MEu in the Tasmanian devil genome. OC1 is present in all carnivorous marsupials, while having a very limited distribution among the remaining Australian marsupial orders. In contrast, hAT-1_MEu is present in all Australian marsupial orders, and has so far only been identified in a few placental mammals. We screened 158 introns for phylogenetically informative retrotransposons in the order Dasyuromorphia, and found that the youngest SINE (Short INterspersed Element), WSINE1, is no longer active in the subfamily Dasyuridae. The lack of detectable WSINE1 activity in this group may be due to a retrotransposon inactivation event approximately 30 Ma. We found that the Tasmanian devil genome contains a relatively low number of continuous full-length LINE-1 (Long INterspersed Element 1, L1) retrotransposons compared with the opossum genome. Furthermore, all L1 elements in the Tasmanian devil appeared to be nonfunctional. Hidden Markov Model approaches suggested that other potential sources of functional reverse transcriptase are absent from the genome. We discuss the issues associated with assembling long, highly similar L1 copies from short read Illumina data and describe how assembly artifacts can potentially lead to erroneous conclusions.

Evolution and gene capture in ancient endogenous retroviruses - insights from the crocodilian genomes

Background

Crocodilians are thought to be hosts to a diverse and divergent complement of endogenous retroviruses (ERVs) but a comprehensive investigation is yet to be performed. The recent sequencing of three crocodilian genomes provides an opportunity for a more detailed and accurate representation of the ERV diversity that is present in these species. Here we investigate the diversity, distribution and evolution of ERVs from the genomes of three key crocodilian species, and outline the key processes driving crocodilian ERV proliferation and evolution.

Results

ERVs and ERV related sequences make up less than 2% of crocodilian genomes. We recovered and described 45 ERV groups within the three crocodilian genomes, many of which are species specific. We have also revealed a new class of ERV, ERV4, which appears to be common to crocodilians and turtles, and currently has no characterised exogenous counterpart. For the first time, we formally describe the characteristics of this ERV class and its classification relative to other recognised ERV and retroviral classes. This class shares some sequence similarity and sequence characteristics with ERV3, although it is phylogenetically distinct from the other ERV classes. We have also identified two instances of gene capture by crocodilian ERVs, one of which, the capture of a host KIT-ligand mRNA has occurred without the loss of an ERV domain.

Conclusions

This study indicates that crocodilian ERVs comprise a wide variety of lineages, many of which appear to reflect ancient infections. In particular, ERV4 appears to have a limited host range, with current data suggesting that it is confined to crocodilians and some lineages of turtles. Also of interest are two ERV groups that demonstrate evidence of host gene capture. This study provides a framework to facilitate further studies into non-mammalian vertebrates and highlights the need for further studies into such species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-014-0071-2) contains supplementary material, which is available to authorized users.

Genomic landscape of human, bat, and ex vivo DNA transposon integrations

The integration and fixation preferences of DNA transposons, one of the major classes of eukaryotic transposable elements, have never been evaluated comprehensively on a genome-wide scale. Here, we present a detailed study of the distribution of DNA transposons in the human and bat genomes. We studied three groups of DNA transposons that integrated at different evolutionary times: 1) ancient (>40 My) and currently inactive human elements, 2) younger (<40 My) bat elements, and 3) ex vivo integrations of piggyBat and Sleeping Beauty elements in HeLa cells. Although the distribution of ex vivo elements reflected integration preferences, the distribution of human and (to a lesser extent) bat elements was also affected by selection. We used regression techniques (linear, negative binomial, and logistic regression models with multiple predictors) applied to 20-kb and 1-Mb windows to investigate how the genomic landscape in the vicinity of DNA transposons contributes to their integration and fixation. Our models indicate that genomic landscape explains 16-79% of variability in DNA transposon genome-wide distribution. Importantly, we not only confirmed previously identified predictors (e.g., DNA conformation and recombination hotspots) but also identified several novel predictors (e.g., signatures of double-strand breaks and telomere hexamer). Ex vivo integrations showed a bias toward actively transcribed regions. Older DNA transposons were located in genomic regions scarce in most conserved elements-likely reflecting purifying selection. Our study highlights how DNA transposons are integral to the evolution of bat and human genomes, and has implications for the development of DNA transposon assays for gene therapy and mutagenesis applications.

Jumping the fine LINE between species: horizontal transfer of transposable elements in animals catalyses genome evolution

Horizontal transfer (HT) is the transmission of genetic material between non-mating species, a phenomenon thought to occur rarely in multicellular eukaryotes. However, many transposable elements (TEs) are not only capable of HT, but have frequently jumped between widely divergent species. Here we review and integrate reported cases of HT in retrotransposons of the BovB family, and DNA transposons, over a broad range of animals spanning all continents. Our conclusions challenge the paradigm that HT in vertebrates is restricted to infective long terminal repeat (LTR) retrotransposons or retroviruses. This raises the possibility that other non-LTR retrotransposons, such as L1 or CR1 elements, believed to be only vertically transmitted, can horizontally transfer between species. Growing evidence indicates that the process of HT is much more general across different TEs and species than previously believed, and that it likely shapes eukaryotic genomes and catalyses genome evolution.

Exaptation of transposable elements into novel cis-regulatory elements: is the evidence always strong?

Transposable elements (TEs) are mobile genetic sequences that can jump around the genome from one location to another, behaving as genomic parasites. TEs have been particularly effective in colonizing mammalian genomes, and such heavy TE load is expected to have conditioned genome evolution. Indeed, studies conducted both at the gene and genome levels have uncovered TE insertions that seem to have been co-opted--or exapted--by providing transcription factor binding sites (TFBSs) that serve as promoters and enhancers, leading to the hypothesis that TE exaptation is a major factor in the evolution of gene regulation. Here, we critically review the evidence for exaptation of TE-derived sequences as TFBSs, promoters, enhancers, and silencers/insulators both at the gene and genome levels. We classify the functional impact attributed to TE insertions into four categories of increasing complexity and argue that so far very few studies have conclusively demonstrated exaptation of TEs as transcriptional regulatory regions. We also contend that many genome-wide studies dealing with TE exaptation in recent lineages of mammals are still inconclusive and that the hypothesis of rapid transcriptional regulatory rewiring mediated by TE mobilization must be taken with caution. Finally, we suggest experimental approaches that may help attributing higher-order functions to candidate exapted TEs.

How do mammalian transposons induce genetic variation? A conceptual framework: the age, structure, allele frequency, and genome context of transposable elements may define their wide-ranging biological impacts

In this essay, we discuss new insights into the wide-ranging impacts of mammalian transposable elements (TE) on gene expression and function. Nearly half of each mammalian genome is comprised of these mobile, repetitive elements. While most TEs are ancient relics, certain classes can move from one chromosomal location to another even now. Indeed, striking recent data show that extensive transposition occurs not only in the germline over evolutionary time, but also in developing somatic tissues and particular human cancers. While occasional germline TE insertions may contribute to genetic variation, many other, similar TEs appear to have little or no impact on neighboring genes. However, the effects of somatic insertions on gene expression and function remain almost completely unknown. We present a conceptual framework to understand how the ages, allele frequencies, molecular structures, and especially the genomic context of mammalian TEs each can influence their various possible functional consequences.

Widespread horizontal transfer of retrotransposons

In higher organisms such as vertebrates, it is generally believed that lateral transfer of genetic information does not readily occur, with the exception of retroviral infection. However, horizontal transfer (HT) of protein coding repetitive elements is the simplest way to explain the patchy distribution of BovB, a long interspersed element (LINE) about 3.2 kb long, that has been found in ruminants, marsupials, squamates, monotremes, and African mammals. BovB sequences are a major component of some of these genomes. Here we show that HT of BovB is significantly more widespread than believed, and we demonstrate the existence of two plausible arthropod vectors, specifically reptile ticks. A phylogenetic tree built from BovB sequences from species in all of these groups does not conform to expected evolutionary relationships of the species, and our analysis indicates that at least nine HT events are required to explain the observed topology. Our results provide compelling evidence for HT of genetic material that has transformed vertebrate genomes.

Mapping the prion protein distribution in marsupials: insights from comparing opossum with mouse CNS

The cellular form of the prion protein (PrP^C) is a sialoglycoprotein widely expressed in the central nervous system (CNS) of mammalian species during neurodevelopment and in adulthood. The location of the protein in the CNS may play a role in the susceptibility of a species to fatal prion diseases, which are also known as the transmissible spongiform encephalopathies (TSEs). To date, little is known about PrP^C distribution in marsupial mammals, for which no naturally occurring prion diseases have been reported. To extend our understanding of varying PrP^C expression profiles in different mammals we carried out a detailed expression analysis of PrP^C distribution along the neurodevelopment of the metatherian South American short-tailed opossum (Monodelphis domestica). We detected lower levels of PrP^C in white matter fiber bundles of opossum CNS compared to mouse CNS. This result is consistent with a possible role for PrP^C in the distinct neurodevelopment and neurocircuitry found in marsupials compared to other mammalian species.

SINEBase: a database and tool for SINE analysis

SINEBase (http://sines.eimb.ru) integrates the revisited body of knowledge about short interspersed elements (SINEs). A set of formal definitions concerning SINEs was introduced. All available sequence data were screened through these definitions and the genetic elements misidentified as SINEs were discarded. As a result, 175 SINE families have been recognized in animals, flowering plants and green algae. These families were classified by the modular structure of their nucleotide sequences and the frequencies of different patterns were evaluated. These data formed the basis for the database of SINEs. The SINEBase website can be used in two ways: first, to explore the database of SINE families, and second, to analyse candidate SINE sequences using specifically developed tools. This article presents an overview of the database and the process of SINE identification and analysis.

Molecular characterization and genetic variability at κ-casein gene (CSN3) in camels

κ-casein is a glycosilated protein belonging to a family of phosphoproteins (αs1, β, αs2, κ) that represent the major protein component in mammalian milk. κ-casein plays an essential role in the casein micelle stabilization, determining the size and the specific function. In the present paper, we report for the first time the characterization of the nucleotide sequence of the whole κ-casein-encoding gene (CSN3) plus 1045 nucleotides at the 5' flanking region in Camelus dromedarius. The promoter region and the complete cDNA were also provided for the first time in Camelus bactrianus. The gene is spread over 9.3kb and consists of 5 exons varying in length from 33bp (exon 3) to 494bp (exon 4), and 4 introns from 1200bp (intron 3) to 2928bp (intron 2). Highly conserved sequences, located in the 5' flanking region, have been found. The regulatory regions of camels seems to be more related to equids than to other compared species. 17 polymorphic sites have been detected, one of these (g.1029T>C) is responsible for the creation of a new putative consensus sequence for the transcription factor HNF-1. In general, these SNPs are the first reported in camels for casein loci. Finally, seven interspersed repeated elements were also identified at intronic level.

LTRsift: a graphical user interface for semi-automatic classification and postprocessing of de novo detected LTR retrotransposons

Background

Long terminal repeat (LTR) retrotransposons are a class of eukaryotic mobile elements characterized by a distinctive sequence similarity-based structure. Hence they are well suited for computational identification. Current software allows for a comprehensive genome-wide de novo detection of such elements. The obvious next step is the classification of newly detected candidates resulting in (super-)families. Such a de novo classification approach based on sequence-based clustering of transposon features has been proposed before, resulting in a preliminary assignment of candidates to families as a basis for subsequent manual refinement. However, such a classification workflow is typically split across a heterogeneous set of glue scripts and generic software (for example, spreadsheets), making it tedious for a human expert to inspect, curate and export the putative families produced by the workflow.

Results

We have developed LTRsift, an interactive graphical software tool for semi-automatic postprocessing of de novo predicted LTR retrotransposon annotations. Its user-friendly interface offers customizable filtering and classification functionality, displaying the putative candidate groups, their members and their internal structure in a hierarchical fashion. To ease manual work, it also supports graphical user interface-driven reassignment, splitting and further annotation of candidates. Export of grouped candidate sets in standard formats is possible. In two case studies, we demonstrate how LTRsift can be employed in the context of a genome-wide LTR retrotransposon survey effort.

Conclusions

LTRsift is a useful and convenient tool for semi-automated classification of newly detected LTR retrotransposons based on their internal features. Its efficient implementation allows for convenient and seamless filtering and classification in an integrated environment. Developed for life scientists, it is helpful in postprocessing and refining the output of software for predicting LTR retrotransposons up to the stage of preparing full-length reference sequence libraries. The LTRsift software is freely available at http://www.zbh.uni-hamburg.de/LTRsift under an open-source license.