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

Chromosome-scale assemblies reveal the structural evolution of African cichlid genomes

BACKGROUND

African cichlid fishes are well known for their rapid radiations and are a model system for studying evolutionary processes. Here we compare multiple, high-quality, chromosome-scale genome assemblies to elucidate the genetic mechanisms underlying cichlid diversification and study how genome structure evolves in rapidly radiating lineages.

RESULTS

We re-anchored our recent assembly of the Nile tilapia (Oreochromis niloticus) genome using a new high-density genetic map. We also developed a new de novo genome assembly of the Lake Malawi cichlid, Metriaclima zebra, using high-coverage Pacific Biosciences sequencing, and anchored contigs to linkage groups (LGs) using 4 different genetic maps. These new anchored assemblies allow the first chromosome-scale comparisons of African cichlid genomes. Large intra-chromosomal structural differences (∼2-28 megabase pairs) among species are common, while inter-chromosomal differences are rare (<10 megabase pairs total). Placement of the centromeres within the chromosome-scale assemblies identifies large structural differences that explain many of the karyotype differences among species. Structural differences are also associated with unique patterns of recombination on sex chromosomes. Structural differences on LG9, LG11, and LG20 are associated with reduced recombination, indicative of inversions between the rock- and sand-dwelling clades of Lake Malawi cichlids. M. zebra has a larger number of recent transposable element insertions compared with O. niloticus, suggesting that several transposable element families have a higher rate of insertion in the haplochromine cichlid lineage.

CONCLUSION

This study identifies novel structural variation among East African cichlid genomes and provides a new set of genomic resources to support research on the mechanisms driving cichlid adaptation and speciation.

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.

Evolution of Two Short Interspersed Elements in Callorhinchus milii (Chondrichthyes, Holocephali) and Related Elements in Sharks and the Coelacanth

Short interspersed elements (SINEs) are non-autonomous retrotransposons. Although they usually show fast evolutionary rates, in some instances highly conserved domains (HCDs) have been observed in elements with otherwise divergent sequences and from distantly related species. Here, we document the life history of two HCD-SINE families in the elephant shark Callorhinchus milii, one specific to the holocephalan lineage (CmiSINEs) and another one (SacSINE1-CM) with homologous elements in sharks and the coelacanth (SacSINE1s, LmeSINE1s). The analyses of their relationships indicated that these elements share the same 3′-tail, which would have allowed both elements to rise to high copy number by exploiting the C. milii L2-2_CM long interspersed element (LINE) enzymes. Molecular clock analysis on SINE activity in C. milii genome evidenced two replication bursts occurring right after two major events in the holocephalan evolution: the end-Permian mass extinction and the radiation of modern Holocephali. Accordingly, the same analysis on the coelacanth homologous elements, LmeSINE1, identified a replication wave close to the split age of the two extant Latimeria species. The genomic distribution of the studied SINEs pointed out contrasting results: some elements were preferentially sorted out from gene regions, but accumulated in flanking regions, while others appear more conserved within genes. Moreover, data from the C. milii transcriptome suggest that these SINEs could be involved in miRNA biogenesis and may be targets for miRNA-based regulation.

Rare horizontal transmission does not hide long-term inheritance of SINE highly conserved domains in the metazoan evolution

Transposable elements (TEs) are self-replicating, mobile DNA sequences which constitute a significant fraction of eukaryotic genomes. They are generally considered selfish DNA, as their replication and random insertion may have deleterious effects on genome functionalities, although some beneficial effects and evolutionary potential have been recognized. Short interspersed elements (SINEs) are non-autonomous TEs with a modular structure: a small RNA-related head, a body, and a long interspersed element-related tail. Despite their high turnover rate and de novo emergence, the body may retain highly conserved domains (HCDs) shared among divergent SINE families: in metazoans, at least nine HCD-SINEs have been recognized. Data mining on public molecular databases allowed the retrieval of 16 new HCD-SINE families from cnidarian, molluscs, arthropods, and vertebrates. Tracking the ancestry of HCDs on the metazoan phylogeny revealed that some of them date back to the Radiata–Bilateria split. Moreover, phylogenetic and age versus divergence analyses of the most ancient HCDs suggested that long-term vertical inheritance is the rule, with few horizontal transfer events. We suggest that the evolutionary conservation of HCDs may be linked to their potential to serve as recombination hotspots. This indirectly affects host genomes by maintaining active and diverse SINE lineages, whose insertions may impact (either positively or negatively) on the evolution of the genome.

RUDI, a short interspersed element of the V-SINE superfamily widespread in molluscan genomes

Short interspersed elements (SINEs) are non-autonomous retrotransposons that are widespread in eukaryotic genomes. They exhibit a chimeric sequence structure consisting of a small RNA-related head, an anonymous body and an AT-rich tail. Although their turnover and de novo emergence is rapid, some SINE elements found in distantly related species retain similarity in certain core segments (or highly conserved domains, HCD). We have characterized a new SINE element named RUDI in the bivalve molluscs Ruditapes decussatus and R. philippinarum and found this element to be widely distributed in the genomes of a number of mollusc species. An unexpected structural feature of RUDI is the HCD domain type V, which was first found in non-amniote vertebrate SINEs and in the SINE from one cnidarian species. In addition to the V domain, the overall sequence conservation pattern of RUDI elements resembles that found in ancient AmnSINE (~310 Myr old) and Au SINE (~320 Myr old) families, suggesting that RUDI might be among the most ancient SINE families. Sequence conservation suggests a monophyletic origin of RUDI. Nucleotide variability and phylogenetic analyses suggest long-term vertical inheritance combined with at least one horizontal transfer event as the most parsimonious explanation for the observed taxonomic distribution.

Lake Tanganyika--a 'melting pot' of ancient and young cichlid lineages (Teleostei: Cichlidae)?

A long history of research focused on the East Africa cichlid radiations (EAR) revealed discrepancies between mtDNA and nuclear phylogenies, suggesting that interspecific hybridisation may have been significant during the radiation of these fishes. The approximately 250 cichlid species of Lake Tanganyika have their roots in a monophyletic African cichlid assemblage, but controversies remain about the precise phylogenetic origin and placement of different lineages and consequently about L. Tanganyika colonization scenarios. 3312 AFLP loci and the mitochondrial ND2 gene were genotyped for 91 species representing almost all major lacustrine and riverine haplotilapiine east African cichlid lineages with a focus on L. Tanganyika endemics. Explicitly testing for the possibility of ancient hybridisation events, a comprehensive phylogenetic network hypothesis is proposed for the origin and diversification of L. Tanganyika cichlids. Inference of discordant phylogenetic signal strongly suggests that the genomes of two endemic L. Tanganyika tribes, Eretmodini and Tropheini, are composed of an ancient mixture of riverine and lacustrine lineages. For the first time a strong monophyly signal of all non-haplochromine mouthbrooding species endemic to L. Tanganyika (“ancient mouthbrooders”) was detected. Further, in the genomes of early diverging L. Tanganyika endemics Trematocarini, Bathybatini, Hemibatini and Boulengerochromis genetic components of other lineages belonging to the East African Radiation appear to be present. In combination with recent palaeo-geological results showing that tectonic activity in the L. Tanganyika region resulted in highly dynamic and heterogeneous landscape evolution over the Neogene and Pleistocene, the novel phylogenetic data render a single lacustrine basin as the geographical cradle of the endemic L. Tanganyika cichlid lineages unlikely. Instead a scenario of a pre-rift origin of several independent L. Tanganyika precursor lineages which diversified in ancient rivers and precursor lakes and then amalgamated in the extant L. Tanganyika basin is put forward as an alternative: the 'melting pot Tanganyika' hypothesis.

Inverse PCR-based method for isolating novel SINEs from genome

Short interspersed elements (SINEs) are moderately repetitive DNA sequences in eukaryotic genomes. Although eukaryotic genomes contain numerous SINEs copy, it is very difficult and laborious to isolate and identify them by the reported methods. In this study, the inverse PCR was successfully applied to isolate SINEs from Opsariichthys bidens genome in Eastern Asian Cyprinid. A group of SINEs derived from tRNA(Ala) molecular had been identified, which were named Opsar according to Opsariichthys. SINEs characteristics were exhibited in Opsar, which contained a tRNA(Ala)-derived region at the 5' end, a tRNA-unrelated region, and AT-rich region at the 3' end. The tRNA-derived region of Opsar shared 76 % sequence similarity with tRNA(Ala) gene. This result indicated that Opsar could derive from the inactive or pseudogene of tRNA(Ala). The reliability of method was tested by obtaining C-SINE, Ct-SINE, and M-SINEs from Ctenopharyngodon idellus, Megalobrama amblycephala, and Cyprinus carpio genomes. This method is simpler than the previously reported, which successfully omitted many steps, such as preparation of probes, construction of genomic libraries, and hybridization.

Conserved domains and SINE diversity during animal evolution

Eukaryotic genomes harbour a number of mobile genetic elements (MGEs); moving from one genomic location to another, they are known to impact on the host genome. Short interspersed elements (SINEs) are well-represented, non-autonomous retroelements and they are likely the most diversified MGEs. In some instances, sequence domains conserved across unrelated SINEs have been identified; remarkably, one of these, called Nin, has been conserved since the Radiata-Bilateria splitting. Here we report on two new domains: Inv, derived from Nin, identified in insects and in deuterostomes, and Pln, restricted to polyneopteran insects. The identification of Inv and Pln sequences allowed us to retrieve new SINEs, two in insects and one in a hemichordate. The diverse structural combination of the different domains in different SINE families, during metazoan evolution, offers a clearer view of SINE diversity and their frequent de novo emergence through module exchange, possibly underlying the high evolutionary success of SINEs.

Characterization of a novel SINE superfamily from invertebrates: "Ceph-SINEs" from the genomes of squids and cuttlefish

Five tRNA-derived short interspersed repetitive elements (SINEs), named SepiaSINE, Sepioth-SINE1, Sepioth-SINE2A, Sepioth-SINE2B and OegopSINE, were isolated from the genomes of three decabrachian species [Sepia officinalis (order Sepiida), Sepiotheuthis lessoniana (suborder Myopsida), and Mastigoteuthis cordiformes (suborder Oegopsida)], by random sequencing and genome screening. In addition, two tRNA-derived SINEs, named IdioSINE1 and IdioSINE2, were further detected from EST (expressed sequence tag) data of Idiosepius paradoxus (order Idiosepiida), using a GenBank FASTA search with a conserved sequence of the SepiaSINE as the query. All the isolated SINEs had a common and unique highly conserved 149-bp sequence in their central structures (Sepioth-SINE2B and IdioSINEs, however, had a continuous 73-bp deletion in the conserved region.), and are therefore grouped as the fourth SINE superfamily "Ceph-SINEs", following the CORE-SINE, V-SINE, and DeuSINE superfamilies. Our analysis suggested that the central conserved region called the "Ceph-domain" might have originated before the diversification of cephalopods (505 myr ago). A sequence alignment of Sepioth-SINE1, Sepioth-SINE2A, and Sepioth-SINE2B demonstrated that Sepioth-SINE2A has a chimeric structure shared with two other SINEs. The above relationship suggests possible template switching in the central conserved domain during reverse transcription for the birth of Sepioth-SINE2A, providing the possibility that the presence of the conserved domain contributed to yield a variety of SINEs during evolution. Furthermore, the distributions of the isolated SINEs showed that order Sepiida, suborders Oegopsida and Myopsida, and order Idiosepiida have their own independent SINE(s), and suggest that order Sepiida can be largely separated into two groups, with clarification of the phylogenetic relatedness between subfamily Sepioteuthinae and the other loliginid squids.

New SINE families from rice, OsSN, with poly(A) at the 3' ends

A database search of the sequences flanking a member of rice retrotransposon RIRE7 revealed that a 298-bp sequence in the region downstream of the member is a repetitive sequence interspersed in the genome of Oryza sativa cv. Nipponbare. Most of the repetitive sequences were flanked by a direct repeat of a target-site sequence, about 14 bp in length. The consensus sequence, 293 bp in length, had no regions encoding any proteins but had sequence motifs of an internal promoter of RNA polymerase III. These indicate that the sequence is a retroposon SINE, designated OsSN1 (Oryza sativa SINE1). OsSN1 is a new rice SINE, because it has no homology with any of the three p-SINE families previously identified from rice, and because it has a stretch of A at the 3' end, unlike p-SINE and any other Gramineae SINEs which have a stretch of T at the 3' end. The Nipponbare genome was found to have many members related to OsSN1, forming two additional new SINE families (designated OsSN2 and OsSN3). OsSN2 and OsSN3 are highly homologous to the 3' and 5' regions of OsSN1, respectively. This suggests that OsSN1 has a mosaic structure, which is generated by sequence exchange (or shuffling) between ancestral OsSN2 and OsSN3. Despite the absence of homology in the 3' regions between OsSN1 (or OsSN2) and OsSN3, a sequence, 5'-TTCTC-3', is commonly present in the region preceding the A stretch at the 3' end. This sequence together with the A stretch and a stem-loop structure found in the region near the A stretch are assumed to be important for retroposition. OsSN members were present in strains of Oryza species, as were p-SINE members. Some of the members showed insertion polymorphism at the respective loci among the rice strains. p-SINE had such polymorphic members, which are useful for classification and phylogenetic analysis of various strains of Oryza species. The polymorphic members of OsSN were more frequently found than those of p-SINE, and therefore, such members are likely to be useful for extensive taxonomic and phylogenetic studies on various rice strains.

Comparative analysis reveals signatures of differentiation amid genomic polymorphism in Lake Malawi cichlids

Low coverage survey sequencing shows that although Lake Malawi cichlids are phenotypically and behaviorally diverse, they appear genetically like a subdivided population.

Background

Cichlid fish from East Africa are remarkable for phenotypic and behavioral diversity on a backdrop of genomic similarity. In 2006, the Joint Genome Institute completed low coverage survey sequencing of the genomes of five phenotypically and ecologically diverse Lake Malawi species. We report a computational and comparative analysis of these data that provides insight into the mechanisms that make closely related species different from one another.

Results

We produced assemblies for the five species ranging in aggregate length from 68 to 79 megabase pairs, identified putative orthologs for more than 12,000 human genes, and predicted more than 32,000 cross-species single nucleotide polymorphisms (SNPs). Nucleotide diversity was lower than that found among laboratory strains of the zebrafish. We collected around 36,000 genotypes to validate a subset of SNPs within and among populations and across multiple individuals of about 75 Lake Malawi species. Notably, there were no fixed differences observed between focal species nor between major lineages. Roughly 3% to 5% of loci surveyed are statistical outliers for genetic differentiation (F_ST) within species, between species, and between major lineages. Outliers for F_ST are candidate genes that may have experienced a history of natural selection in the Malawi lineage.

Conclusion

We present a novel genome sequencing strategy, which is useful when evolutionary diversity is the question of interest. Lake Malawi cichlids are phenotypically and behaviorally diverse, but they appear genetically like a subdivided population. The unique structure of Lake Malawl cichlid genomes should facilitate conceptually new experiments, employing SNPs to identity genotype-phenotype association, using the entire species flock as a mapping panel.

Parasite histories and novel phylogenetic tools: Alternative approaches to inferring parasite evolution from molecular markers

Parasitological research is often contingent on the knowledge of the phylogeny/genealogy of the studied group. Although molecular phylogenetics has proved to be a powerful tool in such investigations, its application in the traditional fashion, based on a tree inference from the primary nucleotide sequences may, in many cases, be insufficient or even improper. These limitations are due to a number of factors, such as a scarcity/ambiguity of phylogenetic information in the sequences, an intricacy of gene relationships at low phylogenetic levels, or a lack of criteria when deciding among several competing coevolutionary scenarios. With respect to the importance of a precise and reliable phylogenetic background in many biological studies, attempts are being made to extend molecular phylogenetics with a variety of new data sources and methodologies. In this review, selected approaches potentially applicable to parasitological research are presented and their advantages as well as drawbacks are discussed. These issues include the usage of idiosyncratic markers (unique features with presumably low probability of homoplasy), such as insertion of mobile elements, gene rearrangements and secondary structure features; the problem of ancestral polymorphism and reticulate relationships at low phylogenetic levels; and the utility of a molecular clock to facilitate discrimination among alternative scenarios in host-parasite coevolution.

The rise and fall of the CR1 subfamily in the lineage leading to penguins

The evolution of penguins has been investigated extensively, although inconclusively, by morphologists, biogeographers and molecular phylogeneticists. We investigated this issue using retroposon analysis of insertions of CR1, which is a member of the LINE (long interspersed element) family, in the genomes of penguins and penguin relatives. The retroposon method is a powerful tool for identifying monophyletic groups. Because retroposons often show different relative frequencies of retroposition during evolution, it is first necessary to identify a certain subgroup that was specifically active during the period when the species in question diverged. Hence, we systematically analyzed many CR1 members isolated from penguin and penguin-related genomes. These CR1s are divided into at least three distinct subgroups that share diagnostic nucleotide insertions and/or deletions, namely, penguin CR1 Sph I, Sph II type A and Sph II type B. The analysis of the inserted retroposons by PCR revealed that different CR1 subfamilies or types had amplified at different rates among different periods during penguin evolution. Namely, the penguin CR1 Sph I subfamily had higher rates of retroposition in a common ancestor of all orders examined in this study or at least in a common ancestor of all extant penguins, and the subfamily Sph II type A also had the same tendency. Therefore, these CR1 members can be used to elucidate the phylogenetic relationships of Sphenisciformes (penguins) among different avian orders. In contrast, the penguin CR1 Sph II type B subfamily had higher rates of retroposition just before and after the emergence of the extant genera in Spheniscidae, suggesting that they are useful for elucidating the intra-relationships among extant penguins. This is the first report for the characterization among the members of CR1 family in avian genomes excluding those of chickens. Hence, this work will be a cornerstone for elucidating the phylogenetic relationships in penguin evolution using the retroposon method.

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.

Parallelism of amino acid changes at the RH1 affecting spectral sensitivity among deep-water cichlids from Lakes Tanganyika and Malawi

Many examples of the appearance of similar traits in different lineages are known during the evolution of organisms. However, the underlying genetic mechanisms have been elucidated in very few cases. Here, we provide a clear example of evolutionary parallelism, involving changes in the same genetic pathway, providing functional adaptation of RH1 pigments to deep-water habitats during the adaptive radiation of East African cichlid fishes. We determined the RH1 sequences from 233 individual cichlids. The reconstruction of cichlid RH1 pigments with 11-cis-retinal from 28 sequences showed that the absorption spectra of the pigments of nine species were shifted toward blue, tuned by two particular amino acid replacements. These blue-shifted RH1 pigments might have evolved as adaptations to the deep-water photic environment. Phylogenetic evidence indicates that one of the replacements, A292S, has evolved several times independently, inducing similar functional change. The parallel evolution of the same mutation at the same amino acid position suggests that the number of genetic changes underlying the appearance of similar traits in cichlid diversification may be fewer than previously expected.

The species flocks of East African cichlid fishes: recent advances in molecular phylogenetics and population genetics

With more than 3,000 species, the fish family Cichlidae is one of the most species-rich families of vertebrates. Cichlids occur in southern and central America, Africa, Madagascar, and India. The hotspot of their biodiversity is East Africa, where they form adaptive radiations composed of hundreds of endemic species in several lakes of various sizes and ages. The unparalleled species richness of East African cichlids has been something of a conundrum for evolutionary biologists and ecologists, since it has been in doubt whether these hundreds of species arose by allopatric speciation or whether it is necessary to invoke somewhat less traditional models of speciation, such as micro-allopatric, peripatric, or even sympatric speciation or evolution through sexual selection mediated by female choice. Ernst Mayr's analyses of these evolutionary uniquely diverse species assemblages have contributed to a more direct approach to this problem and have led to a deeper understanding of the patterns and processes that caused the formation of these huge groups of species. We review here recent molecular data on population differentiation and phylogenetics, which have helped to unravel, to some extent, the patterns and processes that led to the formation and ecological maintenance of cichlid species flocks. It is becoming apparent that sexually selected traits do play an important role in speciation in micro-allopatric or even sympatric settings. Species richness seems to be roughly correlated with the surface area, but not the age, of the lakes. We observe that the oldest lineages of a species flock of cichlids are often less species-rich and live in the open water or deepwater habitats. While the species flocks of the Lake Malawai and the Lake Victoria areas were shown to be monophyletic, the cichlid assemblage of Lake Tanganyika seems to consist of several independent species flocks. Cichlids emerge as an evolutionary model system in which many fundamental questions in evolution and ecology can be tested successfully, yet for other fish species flocks the relative importance of alternative mechanisms of speciation is likely to differ from that in cichlid fish.

First application of the SINE (short interspersed repetitive element) method to infer phylogenetic relationships in reptiles: an example from the turtle superfamily Testudinoidea

Although turtles (order Testudines) constitute one of the major reptile groups, their phylogenetic relationships remain largely unresolved. Hence, we attempted to elucidate their phylogeny using the SINE (short interspersed repetitive element) method, in which the sharing of a SINE at orthologous loci is indicative of synapomorphy. First, a detailed characterization of the tortoise polIII/SINE was conducted using 10 species from eight families of hidden-necked turtles (suborder Cryptodira). Our analysis of 382 SINE sequences newly isolated in the present study revealed two subgroups, namely Cry I and Cry II, which were distinguishable according to diagnostic nucleotides in the 3' region. Furthermore, four (IA-ID) and five (IIA-IIE) different SINE types were identified within Cry I and Cry II subgroups, respectively, based on features of insertions/deletions located in the middle of the SINE sequences. The relative frequency of occurrence of the subgroups and the types of SINEs in this family were highly variable among different lineages of turtles, suggesting active differential retroposition in each lineage. Further application of the SINE method using the most retrotranspositionally active types, namely IB and IC, challenged the established phylogenetic relationships of Bataguridae and its related families. The data for 11 orthologous loci demonstrated a close relationship between Bataguridae and Testudinidae, as well as the presence of the three clades within Bataguridae. Although the SINE method has been used to infer the phylogenies of a number of vertebrate groups, it has never been applied to reptiles. The present study represents the first application of this method to a phylogenetic analysis of this class of vertebrates, and it provides detailed information on the SINE subgroups and types. This information may be applied to the phylogenetic resolution of relevant turtle lineages.

On the roles of repetitive DNA elements in the context of a unified genomic-epigenetic system

Repetitive DNA sequences comprise a substantial portion of most eukaryotic and some prokaryotic chromosomes. Despite nearly forty years of research, the functions of various sequence families as a whole and their monomer units remain largely unknown. The inability to map specific functional roles onto many repetitive DNA elements (REs), coupled with the taxon-specificity of sequence families, have led many to speculate that these genomic components are "selfish" replicators generating genomic "junk." The purpose of this paper is to critically examine the selfishness, evolutionary effects, and functionality of REs. First, a brief overview of the range of ideas pertaining to RE function is presented. Second, the argument is presented that the selfish DNA "hypothesis" is actually a narrative scheme, that it serves to protect neo-Darwinian assumptions from criticism, and that this story is untestable and therefore not a hypothesis. Third, attempts to synthesize the selfish DNA concept with complex systems models of the genome and RE functionality are critiqued. Fourth, the supposed connection between RE-induced mutations and macroevolutionary events are stated to be at variance with empirical evidence and theoretical considerations. Hypotheses that base phylogenetic transitions in repetitive sequence changes thus remain speculative. Fifth and finally, the case is made for viewing REs as integrally functional components of chromosomes, genomes, and cells. It is argued throughout that a new conceptual framework is needed for understanding the roles of repetitive DNA in genomic/epigenetic systems, and that neo-Darwinian "narratives" have been the primary obstacle to elucidating the effects of these enigmatic components of chromosomes.