Analysis of Drosophila species genome size and satellite DNA content reveals significant differences among strains as well as between species

The size of eukaryotic genomes can vary by several orders of magnitude, yet genome size does not correlate with the number of genes nor with the size or complexity of the organism. Although "whole"-genome sequences, such as those now available for 12 Drosophila species, provide information about euchromatic DNA content, they cannot give an accurate estimate of genome sizes that include heterochromatin or repetitive DNA content. Moreover, genome sequences typically represent only one strain or isolate of a single species that does not reflect intraspecies variation. To more accurately estimate whole-genome DNA content and compare these estimates to newly assembled genomes, we used flow cytometry to measure the 2C genome values, relative to Drosophila melanogaster. We estimated genome sizes for the 12 sequenced Drosophila species as well as 91 different strains of 38 species of Drosophilidae. Significant differences in intra- and interspecific 2C genome values exist within the Drosophilidae. Furthermore, by measuring polyploid 16C ovarian follicle cell underreplication we estimated the amount of satellite DNA in each of these species. We found a strong correlation between genome size and amount of satellite underreplication. Addition and loss of heterochromatin satellite repeat elements appear to have made major contributions to the large differences in genome size observed in the Drosophilidae.     

Evolution of protein-coding genes in Drosophila

Several contributing factors have been implicated in evolutionary rate heterogeneity among proteins, but their evolutionary mechanisms remain poorly characterized. The recently sequenced 12 Drosophila genomes provide a unique opportunity to shed light on these unresolved issues. Here, we focus on the role of natural selection in shaping evolutionary rates. We use the Drosophila genomic data to distinguish between factors that increase the strength of purifying selection on proteins and factors that affect the amount of positive selection experienced by proteins. We confirm the importance of translational selection in shaping protein evolution in Drosophila and show that factors such as tissue bias in expression, gene essentiality, intron number, and recombination rate also contribute to evolutionary rate variation among proteins.     

The 19 genomes of Drosophila: a BAC library resource for genus-wide and genome-scale comparative evolutionary research

The genus Drosophila has been the subject of intense comparative phylogenomics characterization to provide insights into genome evolution under diverse biological and ecological contexts and to functionally annotate the Drosophila melanogaster genome, a model system for animal and insect genetics. Recent sequencing of 11 additional Drosophila species from various divergence points of the genus is a first step in this direction. However, to fully reap the benefits of this resource, the Drosophila community is faced with two critical needs: i.e., the expansion of genomic resources from a much broader range of phylogenetic diversity and the development of additional resources to aid in finishing the existing draft genomes. To address these needs, we report the first synthesis of a comprehensive set of bacterial artificial chromosome (BAC) resources for 19 Drosophila species from all three subgenera. Ten libraries were derived from the exact source used to generate 10 of the 12 draft genomes, while the rest were generated from a strategically selected set of species on the basis of salient ecological and life history features and their phylogenetic positions. The majority of the new species have at least one sequenced reference genome for immediate comparative benefit. This 19-BAC library set was rigorously characterized and shown to have large insert sizes (125-168 kb), low nonrecombinant clone content (0.3-5.3%), and deep coverage (9.1-42.9×). Further, we demonstrated the utility of this BAC resource for generating physical maps of targeted loci, refining draft sequence assemblies and identifying potential genomic rearrangements across the phylogeny.     

Conservation of the MHC-like region throughout evolution

Identification of conserved regions between the genomes of distant species is a crucial step in the reconstruction of the genomic organization of their last common ancestor. Here we confirm for the first time with robust evidence, the existence of a region of conserved synteny between the human genome and the Drosophila genome. This evolutionarily conserved synteny involves the human MHC and paralogous regions, and we identified 19 conserved genes between these two species in a Drosophila genomic region of less than 2 Mb. The statistical analysis of the distribution of these 19 genes between the Drosophila and human genomes shows that it cannot be explained by chance. Our study constitutes a first step towards the reconstruction of the genome of Urbilateria (the ancestor of all bilaterian) and allows for a better understanding of the evolutionary history of our genome as well as other metazoan genomes.     

Chromatin domain boundary element search tool for Drosophila

Chromatin domain boundary elements prevent inappropriate interaction between distant or closely spaced regulatory elements and restrict enhancers and silencers to correct target promoters. In spite of having such a general role and expected frequent occurrence genome wide, there is no DNA sequence analysis based tool to identify boundary elements. Here, we report chromatin domain Boundary Element Search Tool (cdBEST), to identify boundary elements. cdBEST uses known recognition sequences of boundary interacting proteins and looks for 'motif clusters'. Using cdBEST, we identified boundary sequences across 12 Drosophila species. Of the 4576 boundary sequences identified in Drosophila melanogaster genome, >170 sequences are repetitive in nature and have sequence homology to transposable elements. Analysis of such sequences across 12 Drosophila genomes showed that the occurrence of repetitive sequences in the context of boundaries is a common feature of drosophilids. We use a variety of genome organization criteria and also experimental test on a subset of the cdBEST boundaries in an enhancer-blocking assay and show that 80% of them indeed function as boundaries in vivo. These observations highlight the role of cdBEST in better understanding of chromatin domain boundaries in Drosophila and setting the stage for comparative analysis of boundaries across closely related species.     

Phylogenetic taxonomy in Drosophila: Problems and prospects

The genus Drosophila is one of the best-studied model systems in modern biology, with twelve fully sequenced genomes available. In spite of the large number of genetic and genomic resources, little is known concerning the phylogenetic relationships, ecology, and evolutionary history of all but a few species. Recent molecular systematic studies have shown that this genus is comprised of at least three independent lineages and that several other genera are actually embedded within Drosophila. This genus accounts for over 2000 described, and many more undescribed, species. While some Drosophila researchers are advocating dividing this genus into three or more separate genera, others favor maintaining Drosophila as a single large genus. With the recent sequencing of the genomes of multiple Drosophila species and their expanding use in comparative biology, it is critical that the Drosophila research community understands the taxonomic framework underlying the naming and relationships of these species. The subdivision of this genus has significant biological implications, ranging from the accurate annotation of single genes to understanding how ecological adaptations have occurred over the history of the group.     

Chromosomal rearrangement inferred from comparisons of 12 Drosophila genomes

The availability of 12 complete genomes of various species of genus Drosophila provides a unique opportunity to analyze genome-scale chromosomal rearrangements among a group of closely related species. This article reports on the comparison of gene order between these 12 species and on the fixed rearrangement events that disrupt gene order. Three major themes are addressed: the conservation of syntenic blocks across species, the disruption of syntenic blocks (via chromosomal inversion events) and its relationship to the phylogenetic distribution of these species, and the rate of rearrangement events over evolutionary time. Comparison of syntenic blocks across this large genomic data set confirms that genetic elements are largely (95%) localized to the same Muller element across genus Drosophila species and paracentric inversions serve as the dominant mechanism for shuffling the order of genes along a chromosome. Gene-order scrambling between species is in accordance with the estimated evolutionary distances between them and we find it to approximate a linear process over time (linear to exponential with alternate divergence time estimates). We find the distribution of synteny segment sizes to be biased by a large number of small segments with comparatively fewer large segments. Our results provide estimated chromosomal evolution rates across this set of species on the basis of whole-genome synteny analysis, which are found to be higher than those previously reported. Identification of conserved syntenic blocks across these genomes suggests a large number of conserved blocks with varying levels of embryonic expression correlation in Drosophila melanogaster. On the other hand, an analysis of the disruption of syntenic blocks between species allowed the identification of fixed inversion breakpoints and estimates of breakpoint reuse and lineage-specific breakpoint event segregation.     

New Drosophila P-like elements and reclassification of Drosophila P-elements subfamilies

Genomic searches for P-like transposable elements were performed (1) in silico in the 12 available Drosophila genomes and (2) by PCR using degenerate primers in 21 Neotropical Drosophila species. In silico searches revealed P-like sequences only in Drosophila persimilis and Drosophila willistoni. Sixteen new P-like elements were obtained by PCR. These sequences were added to sequences of previously described P-like elements, and a phylogenetic analysis was performed. The subfamilies of P-elements described in the literature (Canonical, M, O, T, and K) were included in the reconstructed tree, and all were monophyletic. However, we suggest that some subfamilies can be enlarged, other subdivided, and some new subfamilies may be proposed, totalizing eleven subfamilies, most of which contain new P-like sequences. Our analyses support the monophyly of P-like elements in Drosophilidae. We suggest that, once these elements need host-specific factors to be mobilizable, the horizontal transfer (HT) of P-like elements may be inhibited among more distant taxa. Nevertheless, HT among Drosophilidae species appears to be a common phenomenon.     

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.     

Maintenance of the copy number of retrotransposon MDG3 in the Drosophila melanogaster genome]

The genomes of laboratory stocks and natural population of Drosophila melanogaster contain 8-12 copies of retrotransposon MDG3 detected by in situ hybridization. Construction of genotypes with decreased MDG3 copy number using X-chromosome and chromosome 3 free of MDG3 copies results in appearance of hybrid genomes carrying up to 7-10 copies, instead of 2-4 copies expected. New MDG3 copies are detected in different genome regions, including the 42B hot spot of their location. The chromosomes, where new clusters of MDG3 were observed, carry conserved "parental pattern" of MDG1 arrangement. The data obtained suggest the existence of genomic mechanism for maintenance of retrotransposon copy number on a definite level.     

Parametric alignment of Drosophila genomes

The classic algorithms of Needleman-Wunsch and Smith-Waterman find a maximum a posteriori probability alignment for a pair hidden Markov model (PHMM). To process large genomes that have undergone complex genome rearrangements, almost all existing whole genome alignment methods apply fast heuristics to divide genomes into small pieces that are suitable for Needleman-Wunsch alignment. In these alignment methods, it is standard practice to fix the parameters and to produce a single alignment for subsequent analysis by biologists. As the number of alignment programs applied on a whole genome scale continues to increase, so does the disagreement in their results. The alignments produced by different programs vary greatly, especially in non-coding regions of eukaryotic genomes where the biologically correct alignment is hard to find. Parametric alignment is one possible remedy. This methodology resolves the issue of robustness to changes in parameters by finding all optimal alignments for all possible parameters in a PHMM. Our main result is the construction of a whole genome parametric alignment of Drosophila melanogaster and Drosophila pseudoobscura. This alignment draws on existing heuristics for dividing whole genomes into small pieces for alignment, and it relies on advances we have made in computing convex polytopes that allow us to parametrically align non-coding regions using biologically realistic models. We demonstrate the utility of our parametric alignment for biological inference by showing that cis-regulatory elements are more conserved between Drosophila melanogaster and Drosophila pseudoobscura than previously thought. We also show how whole genome parametric alignment can be used to quantitatively assess the dependence of branch length estimates on alignment parameters.     

Telomere maintenance in Drosophila: Rapid transposon evolution at chromosome ends

The maintenance of terminal sequences is an important role of the telomere, since it prevents the loss of internal regions that encode essential genes. In most eukaryotes, this is accomplished by the telomerase. However, telomere length can also be maintained by other mechanisms, such as homologous recombination and transposition of telomeric retrotransposons to the chromosome ends. A remarkable situation is the case of Drosophila, where telomerase was lost, and thus telomeres managed to be maintained by occasional retrotransposition of telomeric elements to the receding ends. In the recent analysis of 12 Drosophila genomes, ¬¬the multiplicity of autonomous and non-autonomous telomere-specific retrotransposons has revealed extensive and rapid evolution of telomeric DNA. The phylogenetic relationship among these telomeric retrotransposons is congruent with the species phylogeny, suggesting that they have been vertically transmitted from a common ancestor. In this review, we also suggest that the formation of a non-canonical DNA structure at Drosophila telomeres could be the way to protect the ends.     

Sequences homologous to the suffix from the Drosophila melanogaster genome in the DNA of other organisms]

Sequences homologous to the suffix repeat involved in 3'maturation of many Drosophila melanogaster genes have been searched in various organisms by molecular hybridization technique. Several species from evolutionary close and distant to Drosophila melanogaster taxons have been found to contain homology with the suffix in their genomes.     

Characterization of a species-specific satellite DNA of Drosophila buzzatii

The pBuM189 satellite DNA family was found to be species specific for Drosophila buzzatii . It consists of slightly AT-rich tandemly arranged repeats with a high copy number in the genome and shows a very high level of intraspecific sequence similarity. pBuM189 repeats cannot be detected in the genomes of closely related species such as Drosophila serido , Drosophila borborema and Drosophila koepferae . The data support the marginal systematic position of D. buzzatii within the buzzatii cluster of the Drosophila repleta group.     

Wolbachia genomes: insights into an intracellular lifestyle

The genome sequence of the Wolbachia endosymbiont that infects the nematode Brugia malayi has recently been determined together with three partial Wolbachia genomes from different Drosophila species. These data along with the previously published Wolbachia genome from Drosophila melanogaster provide new insights into how this endosymbiont has managed to become so successful.     

Evolutionary perspective on innate immune recognition

Analysis of human and Drosophila genomes demonstrates an ancient origin of innate immunity and the diversity of the mechanisms of innate immune recognition.     

The evolutionary dynamics of the Helena retrotransposon revealed by sequenced Drosophila genomes

Background  

Several studies have shown that genomes contain a mixture of transposable elements, some of which are still active and others ancient relics that have degenerated. This is true for the non-LTR retrotransposon Helena, of which only degenerate sequences have been shown to be present in some species (Drosophila melanogaster), whereas putatively active sequences are present in others (D. simulans). Combining experimental and population analyses with the sequence analysis of the 12 Drosophila genomes, we have investigated the evolution of Helena, and propose a possible scenario for the evolution of this element.     

Genes encoding vitamin-K epoxide reductase are present in Drosophila and trypanosomatid protists

Vitamin-K epoxide reductase is encoded by the VKORC1 gene in mammals and other vertebrates, which also have a paralog, VKORC1L1. Single homologs are present in basal deuterostome and insect genomes, including Drosophila, and three trypanosomatid protists. VKOR is therefore an ancient gene/protein that can be studied in the Drosophila model system.     

The latest buzz in comparative genomics

A second species of fruit fly has just been added to the growing list of organisms with complete and annotated genome sequences. The publication of the Drosophila pseudoobscura sequence provides a snapshot of how genomes have changed over tens of millions of years and sets the stage for the analysis of more fly genomes.     

Genome comparisons highlight similarity and diversity within the eukaryotic kingdoms

In 2000, the number of completely sequenced eukaryotic genomes increased to four. The addition of Drosophila and Arabidopsis into this cohort permits additional insights into the processes that have shaped evolution. Analysis and comparisons of both completed genomes and partially sequenced genomes have already shed light on mechanisms such as gene duplication and gene loss that have long been hypothesized to be major forces in speciation. Indeed, duplicate gene pairs in Saccharomyces, Arabidopsis, Caenorhabditis and Drosophila are high: 30%, 60%, 48% and 40%, respectively. Evidence of horizontal gene-transfer, thought to be a major evolutionary force in bacteria, has been found in Arabidopsis. The release of the 'first draft' of the human genome sequence in 2000 heralds a new stage of biological study. Understanding the as-yet-unannotated human genome will be largely based on conclusions, techniques and tools developed during the analysis and comparison of the genome of these four model organisms.