What transposable elements tell us about genome organization and evolution: the case of Drosophila

Transposable elements (TEs) have been identified in every organism in which they have been looked for. The sequencing of large genomes, such as the human genome and those of Drosophila, Arabidopsis, Caenorhabditis, has also shown that they are a major constituent of these genomes, accounting for 15% of the genome of Drosophila, 45% of the human genome, and more than 70% in some plants and amphibians. Compared with the 1% of genomic DNA dedicated to protein-coding sequences in the human genome, this has prompted various researchers to suggest that the TEs and the other repetitive sequences that constitute the so-called "noncoding DNA", are where the most stimulating discoveries will be made in the future (Bromham, 2002). We are therefore getting further and further from the original idea that this DNA was simply "junk DNA", that owed its presence in the genome entirely to its capacity for selfish transposition. Our understanding of the structures of TEs, their distribution along the genomes, their sequence and insertion polymorphisms within genomes, and within and between populations and species, their impact on genes and on the regulatory mechanisms of genetic expression, their effects on exon shuffling and other phenomena that reshape the genome, and their impact on genome size has increased dramatically in recent years. This leads to a more general picture of the impact of TEs on genomes, though many copies are still mainly selfish or junk DNA. In this review we focus mainly on discoveries made in Drosophila, but we also use information about other genomes when this helps to elucidate the general processes involved in the organization, plasticity, and evolution of genomes.     

Deriving the genomic tree of life in the presence of horizontal gene transfer: conditioned reconstruction

The horizontal gene transfer (HGT) being inferred within prokaryotic genomes appears to be sufficiently massive that many scientists think it may have effectively obscured much of the history of life recorded in DNA. Here, we demonstrate that the tree of life can be reconstructed even in the presence of extensive HGT, provided the processes of genome evolution are properly modeled. We show that the dynamic deletions and insertions of genes that occur during genome evolution, including those introduced by HGT, may be modeled using techniques similar to those used to model nucleotide substitutions that occur during sequence evolution. In particular, we show that appropriately designed general Markov models are reasonable tools for reconstructing genome evolution. These studies indicate that, provided genomes contain sufficiently many genes and that the Markov assumptions are met, it is possible to reconstruct the tree of life. We also consider the fusion of genomes, a process not encountered in gene sequence evolution, and derive a method for the identification and reconstruction of genome fusion events. Genomic reconstructions of a well-defined classical four-genome problem, the root of the multicellular animals, show that the method, when used in conjunction with paralinear/logdet distances, performs remarkably well and is relatively unaffected by the recently discovered big genome artifact.     

The cobweb of life revealed by genome-scale estimates of horizontal gene transfer

With the availability of increasing amounts of genomic sequences, it is becoming clear that genomes experience horizontal transfer and incorporation of genetic information. However, to what extent such horizontal gene transfer (HGT) affects the core genealogical history of organisms remains controversial. Based on initial analyses of complete genomic sequences, HGT has been suggested to be so widespread that it might be the “essence of phylogeny” and might leave the treelike form of genealogy in doubt. On the other hand, possible biased estimation of HGT extent and the findings of coherent phylogenetic patterns indicate that phylogeny of life is well represented by tree graphs. Here, we reexamine this question by assessing the extent of HGT among core orthologous genes using a novel statistical method based on statistical comparisons of tree topology. We apply the method to 40 microbial genomes in the Clusters of Orthologous Groups database over a curated set of 297 orthologous gene clusters, and we detect significant HGT events in 33 out of 297 clusters over a wide range of functional categories. Estimates of positions of HGT events suggest a low mean genome-specific rate of HGT (2.0%) among the orthologous genes, which is in general agreement with other quantitative of HGT. We propose that HGT events, even when relatively common, still leave the treelike history of phylogenies intact, much like cobwebs hanging from tree branches.     

Organizational heterogeneity of vertebrate genomes

Genomes of higher eukaryotes are mosaics of segments with various structural, functional, and evolutionary properties. The availability of whole-genome sequences allows the investigation of their structure as "texts" using different statistical and computational methods. One such method, referred to as Compositional Spectra (CS) analysis, is based on scoring the occurrences of fixed-length oligonucleotides (k-mers) in the target DNA sequence. CS analysis allows generating species- or region-specific characteristics of the genome, regardless of their length and the presence of coding DNA. In this study, we consider the heterogeneity of vertebrate genomes as a joint effect of regional variation in sequence organization superimposed on the differences in nucleotide composition. We estimated compositional and organizational heterogeneity of genome and chromosome sequences separately and found that both heterogeneity types vary widely among genomes as well as among chromosomes in all investigated taxonomic groups. The high correspondence of heterogeneity scores obtained on three genome fractions, coding, repetitive, and the remaining part of the noncoding DNA (the genome dark matter--GDM) allows the assumption that CS-heterogeneity may have functional relevance to genome regulation. Of special interest for such interpretation is the fact that natural GDM sequences display the highest deviation from the corresponding reshuffled sequences.     

Mitochondrial DNA of Vitis vinifera and the issue of rampant horizontal gene transfer

The mitochondrial genome of grape (Vitis vinifera), the largestorganelle genome sequenced so far, is presented. The genomeis 773,279 nt long and has the highest coding capacity amongknown angiosperm mitochondrial DNAs (mtDNAs). The proportionof promiscuous DNA of plastid origin in the genome is also thelargest ever reported for an angiosperm mtDNA, both in absoluteand relative terms. In all, 42.4% of chloroplast genome of Vitishas been incorporated into its mitochondrial genome. In orderto test if horizontal gene transfer (HGT) has also contributedto the gene content of the grape mtDNA, we built phylogenetictrees with the coding sequences of mitochondrial genes of grapeand their homologs from plant mitochondrial genomes. Many incongruentgene tree topologies were obtained. However, the extent of incongruencebetween these gene trees is not significantly greater than thatobserved among optimal trees for chloroplast genes, the commonancestry of which has never been in doubt. In both cases, weattribute this incongruence to artifacts of tree reconstruction,insufficient numbers of characters, and gene paralogy. Thisfinding leads us to question the recent phylogenetic interpretationof Bergthorsson et al. (2003, 2004) and Richardson and Palmer(2007) that rampant HGT into the mtDNA of Amborella best explainsphylogenetic incongruence between mitochondrial gene trees forangiosperms. The only evidence for HGT into the Vitis mtDNAfound involves fragments of two coding sequences stemming fromtwo closteroviruses that cause the leaf roll disease of thisplant. We also report that analysis of sequences shared by bothchloroplast and mitochondrial genomes provides evidence fora previously unknown gene transfer route from the mitochondrionto the chloroplast.     

The fundamental units,processes and patterns of evolution,and the Tree of Life conundrum

Background

The elucidation of the dominant role of horizontal gene transfer (HGT) in the evolution of prokaryotes led to a severe crisis of the Tree of Life (TOL) concept and intense debates on this subject.

Concept

Prompted by the crisis of the TOL, we attempt to define the primary units and the fundamental patterns and processes of evolution. We posit that replication of the genetic material is the singular fundamental biological process and that replication with an error rate below a certain threshold both enables and necessitates evolution by drift and selection. Starting from this proposition, we outline a general concept of evolution that consists of three major precepts.1. The primary agency of evolution consists of Fundamental Units of Evolution (FUEs), that is, units of genetic material that possess a substantial degree of evolutionary independence. The FUEs include both bona fide selfish elements such as viruses, viroids, transposons, and plasmids, which encode some of the information required for their own replication, and regular genes that possess quasi-independence owing to their distinct selective value that provides for their transfer between ensembles of FUEs (genomes) and preferential replication along with the rest of the recipient genome.2. The history of replication of a genetic element without recombination is isomorphously represented by a directed tree graph (an arborescence, in the graph theory language). Recombination within a FUE is common between very closely related sequences where homologous recombination is feasible but becomes negligible for longer evolutionary distances. In contrast, shuffling of FUEs occurs at all evolutionary distances. Thus, a tree is a natural representation of the evolution of an individual FUE on the macro scale, but not of an ensemble of FUEs such as a genome.3. The history of life is properly represented by the "forest" of evolutionary trees for individual FUEs (Forest of Life, or FOL). Search for trends and patterns in the FOL is a productive direction of study that leads to the delineation of ensembles of FUEs that evolve coherently for a certain time span owing to a shared history of vertical inheritance or horizontal gene transfer; these ensembles are commonly known as genomes, taxa, or clades, depending on the level of analysis. A small set of genes (the universal genetic core of life) might show a (mostly) coherent evolutionary trend that transcends the entire history of cellular life forms. However, it might not be useful to denote this trend "the tree of life", or organismal, or species tree because neither organisms nor species are fundamental units of life.

Conclusion

A logical analysis of the units and processes of biological evolution suggests that the natural fundamental unit of evolution is a FUE, that is, a genetic element with an independent evolutionary history. Evolution of a FUE on the macro scale is naturally represented by a tree. Only the full compendium of trees for individual FUEs (the FOL) is an adequate depiction of the evolution of life. Coherent evolution of FUEs over extended evolutionary intervals is a crucial aspect of the history of life but a "species" or "organismal" tree is not a fundamental concept.

Reviewers

This articles was reviewed by Valerian Dolja, W. Ford Doolittle, Nicholas Galtier, and William Martin

     

RIP: the evolutionary cost of genome defense

Repeat-induced point mutation (RIP) is a homology-based process that mutates repetitive DNA and frequently leads to epigenetic silencing of the mutated sequences through DNA methylation. Consistent with the hypothesis that RIP serves to control selfish DNA, an analysis of the Neurospora crassa genome sequence reveals a complete absence of intact mobile elements. As in most eukaryotes, the centromeric regions of N. crassa are rich in sequences that are related to transposable elements; however, in N crassa these sequences have been heavily mutated. The analysis of the N. crassa genome sequence also reveals that RIP has impacted genome evolution significantly through gene duplication, which is considered to be crucial for the evolution of new functions. Most if not all paralogs in N. crassa duplicated and diverged before the emergence of RIP. Thus, RIP illustrates the extraordinary extent to which genomes will go to defend themselves against mobile genetic elements.     

A functional bacteria-derived restriction modification system in the mitochondrion of a heterotrophic protist

The overarching trend in mitochondrial genome evolution is functional streamlining coupled with gene loss. Therefore, gene acquisition by mitochondria is considered to be exceedingly rare. Selfish elements in the form of self-splicing introns occur in many organellar genomes, but the wider diversity of selfish elements, and how they persist in the DNA of organelles, has not been explored. In the mitochondrial genome of a marine heterotrophic katablepharid protist, we identify a functional type II restriction modification (RM) system originating from a horizontal gene transfer (HGT) event involving bacteria related to flavobacteria. This RM system consists of an HpaII-like endonuclease and a cognate cytosine methyltransferase (CM). We demonstrate that these proteins are functional by heterologous expression in both bacterial and eukaryotic cells. These results suggest that a mitochondrion-encoded RM system can function as a toxin–antitoxin selfish element, and that such elements could be co-opted by eukaryotic genomes to drive biased organellar inheritance.

This study reveals that a functional type II restriction modification system of flavobacterial ancestry has been horizontally transferred into the mitochondrion of a marine protist and is capable of encoding potent function, perhaps allowing it to play a role in inter-organellar warfare or protection against further integration of foreign DNA.     

Poxvirus genome evolution by gene gain and loss

The poxviruses (Poxviridae) are a family of viruses with double-stranded DNA genomes and substantial numbers (often >200) of genes per genome. We studied the patterns of gene gain and loss over the evolutionary history of 17 poxvirus complete genomes. A phylogeny based on gene family presence/absence showed good agreement with families based on concatenated amino acid sequences of conserved single-copy genes. Gene duplications in poxviruses were often lineage specific, and the most extensively duplicated viral gene families were found in only a few of the genomes analyzed. A total of 34 gene families were found to include a member in at least one of the poxvirus genomes analyzed and at least one animal genome; in 16 (47%) of these families, there was evidence of recent horizontal gene transfer (HGT) from host to virus. Gene families with evidence of HGT included several involved in host immune defense mechanisms (the MHC class I, interleukin-10, interleukin-24, interleukin-18, the interferon gamma receptor, and tumor necrosis factor receptor II) and others (glutaredoxin and glutathione peroxidase) involved in resistance of cells to oxidative stress. Thus "capture" of host genes by HGT has been a recurrent feature of poxvirus evolution and has played an important role in adapting the virus to survive host antiviral defense mechanisms.     

WindowMasker: window-based masker for sequenced genomes

MOTIVATION: Matches to repetitive sequences are usually undesirable in the output of DNA database searches. Repetitive sequences need not be matched to a query, if they can be masked in the database. RepeatMasker/Maskeraid (RM), currently the most widely used software for DNA sequence masking, is slow and requires a library of repetitive template sequences, such as a manually curated RepBase library, that may not exist for newly sequenced genomes. RESULTS: We have developed a software tool called WindowMasker (WM) that identifies and masks highly repetitive DNA sequences in a genome, using only the sequence of the genome itself. WM is orders of magnitude faster than RM because WM uses a few linear-time scans of the genome sequence, rather than local alignment methods that compare each library sequence with each piece of the genome. We validate WM by comparing BLAST outputs from large sets of queries applied to two versions of the same genome, one masked by WM, and the other masked by RM. Even for genomes such as the human genome, where a good RepBase library is available, searching the database as masked with WM yields more matches that are apparently non-repetitive and fewer matches to repetitive sequences. We show that these results hold for transcribed regions as well. WM also performs well on genomes for which much of the sequence was in draft form at the time of the analysis. AVAILABILITY: WM is included in the NCBI C++ toolkit. The source code for the entire toolkit is available at ftp://ftp.ncbi.nih.gov/toolbox/ncbi_tools++/CURRENT/. Once the toolkit source is unpacked, the instructions for building WindowMasker application in the UNIX environment can be found in file src/app/winmasker/README.build. SUPPLEMENTARY INFORMATION: Supplementary data are available at ftp://ftp.ncbi.nlm.nih.gov/pub/agarwala/windowmasker/windowmasker_suppl.pdf     

Molecular analysis of the genomic organization of the Hymenoptera Diadromus pulchellus and Eupelmus vuilleti

The genomic organization of two parasitic wasps was analyzed by DNA reassociation. Cot curves revealed a pattern with three types of components. A highly repetitive DNA, accounting for 15 to 25% of the genome, was identified as satellite DNA. The moderately repetitive DNA corresponds to 26 to 42% of the genome in both species, and shows large variations in complexity, repetitive frequency and a number of sub-components between males and females. These variations are seen as resulting from DNA amplification during somatic and sexual differentiation. Dot blot analyses show that such DNA amplifications concern several types of structural and regulatory genes. The presence of repeated mobile elements was studied by the Roninson method to compare the repeated sequence patterns of Diadromus pulchellus and Eupelmus vuilleti with those of Drosophila melanogaster. The occurrence and organization of mobile elements in these Hymenoptera differ from those of the neighboring order of Diptera. The repetitive and unique components define very large genomes (1 to 3 × 10⁹ base pairs). The genomic organization in Parasitica appears to be an extreme drosophilan type. We propose that the germinal genome of these parasitic wasps is primarily composed of satellite DNA blocks and very long stretches of unique sequences, separated by a few repeated and/or variously deleted, interspersed elements of each mobile element family.     

Repetitive DNA variation and pivotal-differential evolution of wild wheats.

Several polyploid species in the genus Triticum contain a U genome derived from the diploid T. umbellulatum. In these species, the U genome is considered to be unmodified from the diploid based on chromosome pairing analysis, and it is referred to as pivotal. The additional genome(s) are considered to be modified, and they are thus referred to as differential genomes. The M genome derived from the diploid T. comosum is found in many U genome polyploids. In this study, we cloned three repetitive DNA sequences found primarily in the U genome and two repetitive DNA sequences found primarily in the M genome. We used these to monitor variation for these sequences in a large set of species containing U and M genomes. Investigation of sympatric and allopatric accessions of polyploid species did not show repetitive DNA similarities among sympatric species. This result does not support the idea that the polyploid species are continually exchanging genetic information through introgression. However, it is also possible that repetitive DNA is not a suitable means of addressing the question of introgression. The U genomes of both diploid and polyploid U genome species were similar regarding hybridization patterns observed with U genome probes. Much more variation was found both among diploid T. comosum accessions and polyploids containing M genomes. The observed variation supports the cytogenetic evidence that the M genome is more variable than the U genome. It also raises the possibility that the differential nature of the M genome may be due to variation within the diploid T. comosum, as well as among polyploid M genome species and accessions.     

A novel repetitive DNA sequence in the genus Oryza

Summary Repetitive DNA sequences in the genus Oryza (rice) represent a large fraction of the nuclear DNA. The isolation and characterization of major repetitive DNA sequences will lead to a better understanding of rice genome organization and evolution. Here we report the characterization of a novel repetitive sequence, CC-1, from the CC genome. This repetitive sequence is present as long tandem arrays with a repeat unit 194 bp in length in the CC-diploid genome but 172 bp in length in the BBCC and CCDD tetraploid genomes. This repetitive sequence is also present, though at lower copy numbers, in the AA and BB genomes, but is absent in the EE and FF genomes. Hybridization experiments revealed considerable differences both in copy numbers and in restriction fragment patterns of CC-1 both between and within rice species. The results support the hypothesis that the CC genome is more closely related to the AA genome than to the BB genome, and most distantly related to the EE and FF genomes.     

Reiterated DNA sequences in Rhizobium and Agrobacterium spp.

Repeated DNA sequences are a general characteristic of eucaryotic genomes. Although several examples of DNA reiteration have been found in procaryotic organisms, only in the case of the archaebacteria Halobacterium halobium and Halobacterium volcanii [C. Sapienza and W. F. Doolittle, Nature (London) 295:384-389, 1982], has DNA reiteration been reported as a common genomic feature. The genomes of two Rhizobium phaseoli strains, one Rhizobium meliloti strain, and one Agrobacterium tumefaciens strain were analyzed for the presence of repetitive DNA. Rhizobium and Agrobacterium spp. are closely related soil bacteria that interact with plants and that belong to the taxonomical family Rhizobiaceae. Rhizobium species establish a nitrogen-fixing symbiosis in the roots of legumes, whereas Agrobacterium species is a pathogen in different plants. The four strains revealed a large number of repeated DNA sequences. The family size was usually small, from 2 to 5 elements, but some presented more than 10 elements. Rhizobium and Agrobacterium spp. contain large plasmids in addition to the chromosomes. Analysis of the two Rhizobium strains indicated that DNA reiteration is not confined to the chromosome or to some plasmids but is a property of the whole genome.     

Isolation of A/D and C genome specific dispersed and clustered repetitive DNA sequences from Avena sativa.

DNA gel-blot and in situ hybridization with genome-specific repeated sequences have proven to be valuable tools in analyzing genome structure and relationships in species with complex allopolyploid genomes such as hexaploid oat (Avena sativa L., 2n = 6x = 42; AACCDD genome). In this report, we describe a systematic approach for isolating genome-, chromosome-, and region-specific repeated and low-copy DNA sequences from oat that can presumably be applied to any complex genome species. Genome-specific DNA sequences were first identified in a random set of A. sativa genomic DNA cosmid clones by gel-blot hybridization using labeled genomic DNA from different Avena species. Because no repetitive sequences were identified that could distinguish between the A and D gneomes, sequences specific to these two genomes are refereed to as A/D genome specific. A/D or C genome specific DNA subfragments were used as screening probes to identify additional genome-specific cosmid clones in the A. sativa genomic library. We identified clustered and dispersed repetitive DNA elements for the A/D and C genomes that could be used as cytogenetic markers for discrimination of the various oat chromosomes. Some analyzed cosmids appeared to be composed entirely of genome-specific elements, whereas others represented regions with genome- and non-specific repeated sequences with interspersed low-copy DNA sequences. Thus, genome-specific hybridization analysis of restriction digests of random and selected A. sativa cosmids also provides insight into the sequence organization of the oat genome.     

Mobile DNAs: the poacher turned gamekeeper

The selfish DNA hypothesis imagines the genome as an ecological community, a collection of interacting DNA sequences with differing evolutionary origins and potentially different interests. We are now finding out more about the ways in which host sequences can enlist the help of formerly parasitic DNAs.