Chromosomal elements evolve at different rates in the Drosophila genome

Computational Sequence Analysis of Inversion Breakpoint Regions in the Cactophilic Drosophila mojavensis Lineage

We investigated rates of chromosomal evolution in Drosophila mojavensis using whole-genome sequence information from D. mojavensis, Drosophila buzzatii, and Drosophila virilis. Drosophila mojavensis is a cactophilic species of the repleta group living under extreme ecological conditions in the deserts of the Southwestern United States and Northwestern México. The genome of D. buzzatii, another member of the repleta group, was recently sequenced and the largest scaffolds anchored to all chromosomes using diverse procedures. Chromosome organization between D. mojavensis and D. buzzatii was compared using MUMmer and GRIMM software. Our results corroborate previous cytological analyses that indicated chromosome 2 differed between these 2 species by 10 inversions, chromosomes X and 5 differed by one inversion each, and chromosome 4 was homosequential. In contrast, we found that chromosome 3 differed by 5 inversions instead of the expected 2 that were previously inferred by cytological analyses. Thirteen of these inversions occurred in the D. mojavensis lineage: 12 are fixed and one of them is a polymorphic inversion previously described in populations from Sonora and Baja California, México. We previously investigated the breakpoints of chromosome 2 inversions fixed in D. mojavensis. Here we characterized the breakpoint regions of the 5 inversions found in chromosome 3 in order to infer the molecular mechanism that generated each inversion and its putative functional consequences. Overall, our results reveal a number of gene alterations at the inversion breakpoints with putative adaptive consequences that point to natural selection as the cause for fast chromosomal evolution in D. mojavensis.

Muller "Elements" in Drosophila: How the Search for the Genetic Basis for Speciation Led to the Birth of Comparative Genomics

The concept of synteny, or conservation of genes on the same chromosome, traces its origins to the early days of Drosophila genetics. This discovery emerged from comparisons of linkage maps from different species of Drosophila with the goal of understanding the process of speciation. H. J. Muller published a landmark article entitled Bearings of the "Drosophila" work on systematics, where he synthesized genetic and physical map data and proposed a model of speciation and chromosomal gene content conservation. These models have withstood the test of time with the advent of molecular genetic analysis from protein to genome level variation. Muller's ideas provide a framework to begin to answer questions about the evolutionary forces that shape the structure of the genome.

Inversion evolutionary rates might limit the experimental identification of inversion breakpoints in non-model species

Chromosomal inversions are structural changes that alter gene order but generally not gene content in the affected region. In Drosophila, extensive cytological studies revealed the widespread character of inversion polymorphism, with evidence for its adaptive character. In Drosophila subobscura, polymorphism affects both its four large autosomal elements and its X (A) chromosome. The characterization of eight of these autosomal inversions breakpoints revealed that most of them originated through the staggered-breaks mechanism. Here, we have performed chromosomal walks to identify the breakpoints of two X-chromosome widely distributed inversions -A₂ and A₁- of D. subobscura. Inversion A₂ is considered a warm-adapted arrangement that exhibits parallel latitudinal clines in the species ancestral distribution area and in both American subcontinents, whereas inversion A₁ is only present in the Palearctic region where it presents an east-west cline. The duplication detected at the A₂ inversion breakpoints is consistent with its origin by the staggered-breaks mechanism. Inversion A₁ breakpoints could not be molecularly identified even though they could be narrowly delimited. This result points to chromosome walking limitations when using as a guide the genome of other species. Limitations stem from the rate of evolution by paracentric inversions, which in Drosophila is highest for the X chromosome.

Genomics of parallel adaptation at two timescales in Drosophila

Two interesting unanswered questions are the extent to which both the broad patterns and genetic details of adaptive divergence are repeatable across species, and the timescales over which parallel adaptation may be observed. Drosophila melanogaster is a key model system for population and evolutionary genomics. Findings from genetics and genomics suggest that recent adaptation to latitudinal environmental variation (on the timescale of hundreds or thousands of years) associated with Out-of-Africa colonization plays an important role in maintaining biological variation in the species. Additionally, studies of interspecific differences between D. melanogaster and its sister species D. simulans have revealed that a substantial proportion of proteins and amino acid residues exhibit adaptive divergence on a roughly few million years long timescale. Here we use population genomic approaches to attack the problem of parallelism between D. melanogaster and a highly diverged conger, D. hydei, on two timescales. D. hydei, a member of the repleta group of Drosophila, is similar to D. melanogaster, in that it too appears to be a recently cosmopolitan species and recent colonizer of high latitude environments. We observed parallelism both for genes exhibiting latitudinal allele frequency differentiation within species and for genes exhibiting recurrent adaptive protein divergence between species. Greater parallelism was observed for long-term adaptive protein evolution and this parallelism includes not only the specific genes/proteins that exhibit adaptive evolution, but extends even to the magnitudes of the selective effects on interspecific protein differences. Thus, despite the roughly 50 million years of time separating D. melanogaster and D. hydei, and despite their considerably divergent biology, they exhibit substantial parallelism, suggesting the existence of a fundamental predictability of adaptive evolution in the genus.

Both local adaptation on short timescales and the long-term accumulation of adaptive differences between species have recently been investigated using comparative genomic and population genomic approaches in several species. However, the repeatability of adaptive evolution at the genetic level is poorly understood. Here we attack this problem by comparing patterns of long and short-term adaptation in Drosophila melanogaster to patterns of adaptation on two timescales in a highly diverged congener, Drosophila hydei. We found, despite the fact that these species diverged from a common ancestor roughly 50 million years ago, the population genomics of latitudinal allele frequency differentiation shows that there is a substantial shared set of genes likely playing a role in the short term adaptive divergence of populations in both species. Analyses of longer-term adaptive protein divergence for the D. hydei-D. mojavensis and D. melanogaster-D. simulans clades reveal a striking level of parallel adaptation. This parallelism includes not only the specific genes/proteins that exhibit adaptive evolution, but extends even to the magnitudes of the selective effects on interspecific protein differences.

Dense gene physical maps of the non-model species Drosophila subobscura

The comparative analysis of genetic and physical maps as well as of whole genome sequences had revealed that in the Drosophila genus, most structural rearrangements occurred within chromosomal elements as a result of paracentric inversions. Genome sequence comparison would seem the best method to estimate rates of chromosomal evolution, but the high-quality reference genomes required for this endeavor are still scanty. Here, we have obtained dense physical maps for Muller elements A, C, and E of Drosophila subobscura, a species with an extensively studied rich and adaptive chromosomal polymorphism. These maps are based on 462 markers: 115, 236, and 111 markers for elements A, C, and E, respectively. The availability of these dense maps will facilitate genome assembly and will thus greatly contribute to obtaining a good reference genome, which is a required step for D. subobscura to attain the model species status. The comparative analysis of these physical maps and those obtained from the D. pseudoobscura and D. melanogaster genomes allowed us to infer the number of fixed inversions and chromosomal evolutionary rates for each pairwise comparison. For all three elements, rates inferred from the more closely related species were higher than those inferred from the more distantly related species, which together with results of relative-rate tests point to an acceleration in the D. subobscura lineage at least for elements A and E.

Chromosome evolution in malaria mosquitoes inferred from physically mapped genome assemblies

Polymorphic inversions in mosquitoes are distributed nonrandomly among chromosomes and are associated with ecological, behavioral, and physiological adaptations related to pathogen transmission. Despite their significance, the patterns and mechanism of genome rearrangements are not well understood. Recent sequencing and physical mapping of the genomes for 16 Anopheles mosquito species provided an opportunity to study chromosome evolution at the highest resolution. New studies revealed that fixed rearrangement accumulated [Formula: see text]3 times faster on the X chromosome than on autosomes. The highest densities of transposable elements (TEs) and satellites of different sizes have also been found on the X chromosome, suggesting a mechanism for the inversion generation. The high rate of X chromosome rearrangements is in sharp contrast with the paucity of polymorphic inversions on the X in the majority of anopheline species. This paper highlights the advances in understanding chromosome evolution in malaria vectors and discusses possible future directions in studying mechanisms and biological roles of genome rearrangements.

Genomic composition and evolution of Aedes aegypti chromosomes revealed by the analysis of physically mapped supercontigs

BACKGROUND

An initial comparative genomic study of the malaria vector Anopheles gambiae and the yellow fever mosquito Aedes aegypti revealed striking differences in the genome assembly size and in the abundance of transposable elements between the two species. However, the chromosome arms homology between An. gambiae and Ae. aegypti, as well as the distribution of genes and repetitive elements in chromosomes of Ae. aegypti, remained largely unexplored because of the lack of a detailed physical genome map for the yellow fever mosquito.

RESULTS

Using a molecular landmark-guided fluorescent in situ hybridization approach, we mapped 624 Mb of the Ae. aegypti genome to mitotic chromosomes. We used this map to analyze the distribution of genes, tandem repeats and transposable elements along the chromosomes and to explore the patterns of chromosome homology and rearrangements between Ae. aegypti and An. gambiae. The study demonstrated that the q arm of the sex-determining chromosome 1 had the lowest gene content and the highest density of minisatellites. A comparative genomic analysis with An. gambiae determined that the previously proposed whole-arm synteny is not fully preserved; a number of pericentric inversions have occurred between the two species. The sex-determining chromosome 1 had a higher rate of genome rearrangements than observed in autosomes 2 and 3 of Ae. aegypti.

CONCLUSIONS

The study developed a physical map of 45% of the Ae. aegypti genome and provided new insights into genomic composition and evolution of Ae. aegypti chromosomes. Our data suggest that minisatellites rather than transposable elements played a major role in rapid evolution of chromosome 1 in the Aedes lineage. The research tools and information generated by this study contribute to a more complete understanding of the genome organization and evolution in mosquitoes.

Differential occurrence of chromosome inversion polymorphisms among Muller's elements in three species of the tripunctata group of Drosophila, including a species with fast chromosomal evolution

Detailed chromosome maps with reliable homologies among chromosomes of different species are the first step to study the evolution of the genetic architecture in any set of species. Here, we present detailed photo maps of the polytene chromosomes of three closely related species of the tripunctata group (subgenus Drosophila): Drosophila mediopunctata, D. roehrae, and D. unipunctata. We identified Muller's elements in each species, using FISH, establishing reliable chromosome homologies among species and D. melanogaster. The simultaneous analysis of chromosome inversions revealed a distribution pattern for the inversion polymorphisms among Muller's elements in the three species. Element E is the most polymorphic, with many inversions in each species. Element C follows; while the least polymorphic elements are B and D. While interesting, it remains to be determined how general this pattern is among species of the tripunctata group. Despite previous studies showing that D. mediopunctata and D. unipunctata are phylogenetically closer to each other than to D. roehrae, D. unipunctata shows rare karyotypic changes. It has two chromosome fusions: an additional heterochromatic chromosome pair and a pericentric inversion in the X chromosome. This especial conformation suggests a fast chromosomal evolution that deserves further study.

Testing chromosomal phylogenies and inversion breakpoint reuse in Drosophila. The martensis cluster revisited

The chromosomal relationships of the four martensis cluster species are among the most complex and intricate within the entire Drosophila repleta group, due to the so-called sharing of inversions. Here, we have revised these relationships using comparative mapping of bacterial artificial chromosome (BAC) clones on the salivary gland chromosomes. A physical map of chromosome 2 of Drosophila uniseta (one of the cluster members) was generated by in situ hybridization of 82 BAC clones from the physical map of the Drosophila buzzatii genome (an outgroup that represents the ancestral arrangement). By comparing the marker positions, we determined the number, order, and orientation of conserved chromosomal segments between chromosome 2 of D. buzzatii and D. uniseta. GRIMM software was used to infer that a minimum of five chromosomal inversions are necessary to transform the chromosome 2 of D. buzzatii into that of D. uniseta. Two of these inversions have been overlooked in previous cytological analyses. The five fixed inversions entail two breakpoint reuses because only nine syntenic segments and eight interruptions were observed. We tested for the presence of the five inversions fixed in D. uniseta in the other three species of the martensis cluster by in situ hybridization of eight breakpoint-bearing BAC clones. The results shed light on the chromosomal phylogeny of the martensis cluster, yet leave a number of questions open.

A physical map for an Asian malaria mosquito, Anopheles stephensi

Physical mapping is a useful approach for studying genome organization and evolution as well as for genome sequence assembly. The availability of polytene chromosomes in malaria mosquitoes provides a unique opportunity to develop high-resolution physical maps. We report a 0.6-Mb-resolution physical map consisting of 422 DNA markers hybridized to 379 chromosomal sites of the Anopheles stephensi polytene chromosomes. This makes An. stephensi second only to Anopheles gambiae in density of a physical map among malaria mosquitoes. Three hundred sixty-three (363) probes hybridized to single chromosomal sites, whereas 59 clones yielded multiple signals. This physical map provided a suitable basis for comparative genomics, which was used for determining inversion breakpoints, duplications, and origin of novel genes across species.

Genome landscape and evolutionary plasticity of chromosomes in malaria mosquitoes

BACKGROUND

Nonrandom distribution of rearrangements is a common feature of eukaryotic chromosomes that is not well understood in terms of genome organization and evolution. In the major African malaria vector Anopheles gambiae, polymorphic inversions are highly nonuniformly distributed among five chromosomal arms and are associated with epidemiologically important adaptations. However, it is not clear whether the genomic content of the chromosomal arms is associated with inversion polymorphism and fixation rates.

METHODOLOGY/PRINCIPAL FINDINGS

To better understand the evolutionary dynamics of chromosomal inversions, we created a physical map for an Asian malaria mosquito, Anopheles stephensi, and compared it with the genome of An. gambiae. We also developed and deployed novel Bayesian statistical models to analyze genome landscapes in individual chromosomal arms An. gambiae. Here, we demonstrate that, despite the paucity of inversion polymorphisms on the X chromosome, this chromosome has the fastest rate of inversion fixation and the highest density of transposable elements, simple DNA repeats, and GC content. The highly polymorphic and rapidly evolving autosomal 2R arm had overrepresentation of genes involved in cellular response to stress supporting the role of natural selection in maintaining adaptive polymorphic inversions. In addition, the 2R arm had the highest density of regions involved in segmental duplications that clustered in the breakpoint-rich zone of the arm. In contrast, the slower evolving 2L, 3R, and 3L, arms were enriched with matrix-attachment regions that potentially contribute to chromosome stability in the cell nucleus.

CONCLUSIONS/SIGNIFICANCE

These results highlight fundamental differences in evolutionary dynamics of the sex chromosome and autosomes and revealed the strong association between characteristics of the genome landscape and rates of chromosomal evolution. We conclude that a unique combination of various classes of genes and repetitive DNA in each arm, rather than a single type of repetitive element, is likely responsible for arm-specific rates of rearrangements.

FISH mapping of microsatellite loci from Drosophila subobscura and its comparison to related species

Microsatellites are highly polymorphic markers that are distributed through all the genome being more abundant in non-coding regions. Whether they are neutral or under selection, these markers if localized can be used as co-dominant molecular markers to explore the dynamics of the evolutionary processes. Their cytological localization can allow identifying genes under selection, inferring recombination from a genomic point of view, or screening for the genomic reorganizations occurring during the evolution of a lineage, among others. In this paper, we report for the first time the localization of microsatellite loci by fluorescent in situ hybridization on Drosophila polytene chromosomes. In Drosophila subobscura, 72 dinucleotide microsatellite loci were localized by fluorescent in situ hybridization yielding unique hybridization signals. In the sex chromosome, microsatellite distribution was not uniform and its density was higher than in autosomes. We identified homologous segments to the sequence flanking the microsatellite loci by browsing the genome sequence of Drosophila pseudoobscura and Drosophila melanogaster. Their localization supports the conservation of Muller's chromosomal elements among Drosophila species and the existence of multiple intrachromosomal rearrangements within each evolutionary lineage. Finally, the lack of microsatellite repeats in the homologous D. melanogaster sequences suggests convergent evolution for high microsatellite density in the distal part of the X chromosome.

Cloning and sequencing of the breakpoint regions of inversion 5g fixed in Drosophila buzzatii

Chromosomal inversions are ubiquitous in Drosophila both as intraspecific polymorphisms and interspecific differences. Many gaps still remain in our understanding of the mechanisms that generate them. Previous work has shown that in Drosophila buzzatii, three polymorphic inversions were generated by ectopic recombination between copies of the transposon Galileo. In this study, we have characterized the breakpoint regions of inversion 5g, fixed in D. buzzatii and absent in Drosophila koepferae and other closely related species. A novel approach comprising four experimental steps was used. First, D. buzzatii BAC clones encompassing the breakpoints were identified and their ends sequenced. Then, breakpoint regions were mapped at high resolution in the Drosophila mojavensis genome sequence. Finally, breakpoint regions were isolated by polymerase chain reaction in D. buzzatii and D. koepferae and sequenced. Our aim was to shed light on the mechanism that generated inversion 5g and specifically to test for an implication of the transposon Galileo. No evidence implicates Galileo or other transposable elements in the origin of inversion 5g that was generated most likely by two independent breaks and non-homologous end-joining repair. Our results show that different inversion-generating mechanisms may coexist within the same lineage and suggest a hypothesis for the evolutionary time and mode of their operation.

BAC library construction and BAC end sequencing of five Drosophila species: the comparative map with the D. melanogaster genome

We constructed and characterized arrayed bacterial artificial chromosome (BAC) libraries of five Drosophila species (D. melanogaster, D. simulans, D. sechellia, D. auraria, and D. ananassae), which are genetically well characterized in the studies of meiosis, evolution, population genetics, and developmental biology. The BAC libraries comprise 8,000 to 12,500 clones for each species, estimated to cover the most of the genomes. We sequenced both ends of most of these BAC clones with a success rate of 91%. Of these, 53,701 clones consisting of non-repetitive BAC end sequences (BESs) were mapped with reference of the public D. melanogaster genome sequences. The BES mapping estimated that the BAC libraries of D. auraria and D. ananassae covered 47% and 57% of the D. melanogaster genome, respectively, and those of D. melanogaster, D. sechellia, and D. simulans covered 94-97%. The low coverage by BESs of D. auraria and D. ananassae may be due to the high sequence divergence with D. melanogaster. From the comparative BES mapping, 111 possible breakpoints of chromosomal rearrangements were identified in these four species. The breakpoints of the major chromosome rearrangement between D. simulans and D. melanogaster on the third chromosome were determined within 20 kb in 84E and 30 kb in 93E/F. Corresponding breakpoints were also identified in D. sechellia. The BAC clones described here will be an important addition to the Drosophila genomic resources.

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.

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.

Independent origins of new sex-linked chromosomes in the melanica and robusta species groups of Drosophila

Background

Recent translocations of autosomal regions to the sex chromosomes represent important systems for identifying the evolutionary forces affecting convergent patterns of sex-chromosome heteromorphism. Additions to the sex chromosomes have been reported in the melanica and robusta species groups, two sister clades of Drosophila. The close relationship between these two species groups and the similarity of their rearranged karyotypes motivates this test of alternative hypotheses; the rearranged sex chromosomes in both groups are derived through a common origin, or the rearrangements are derived through at least two independent origins. Here we examine chromosomal arrangement in representatives of the melanica and the robusta species groups and test these alternative hypotheses using a phylogenetic approach.

Results

Two mitochondrial and two nuclear gene sequences were used to reconstruct phylogenetic relationships of a set of nine ingroup species having fused and unfused sex chromosomes and representing a broad sample of both species groups. Different methods of phylogenetic inference, coupled with concurrent cytogenetic analysis, indicate that the hypothesis of independent origins of rearranged sex chromosomes within each species group is significantly more likely than the alternative hypothesis of a single common origin. An estimate tightly constrained around 8 My was obtained for the age of the rearranged sex chromosomes in the melanica group; however, a more loosely constrained estimate of 10–15 My was obtained for the age of the rearrangement in the robusta group.

Conclusion

Independent acquisition of new chromosomal arms by the sex chromosomes in the melanica and robusta species groups represents a case of striking convergence at the karyotypic level. Our findings indicate that the parallel divergence experienced by newly sex-linked genomic regions in these groups represents an excellent system for studying the tempo of sex chromosome evolution.

A BAC-based integrated linkage map of the silkworm Bombyx mori

An integrated map of the Bombyx mori genome has been constructed using 361.1 Mb of BAC contigs and singletons together with a genetic map containing 1689 independent genes and synteny among Apis, Tribolium, and Bombyx was examined.

Background

In 2004, draft sequences of the model lepidopteran Bombyx mori were reported using whole-genome shotgun sequencing. Because of relatively shallow genome coverage, the silkworm genome remains fragmented, hampering annotation and comparative genome studies. For a more complete genome analysis, we developed extended scaffolds combining physical maps with improved genetic maps.

Results

We mapped 1,755 single nucleotide polymorphism (SNP) markers from bacterial artificial chromosome (BAC) end sequences onto 28 linkage groups using a recombining male backcross population, yielding an average inter-SNP distance of 0.81 cM (about 270 kilobases). We constructed 6,221 contigs by fingerprinting clones from three BAC libraries digested with different restriction enzymes, and assigned a total of 724 single copy genes to them by BLAST (basic local alignment search tool) search of the BAC end sequences and high-density BAC filter hybridization using expressed sequence tags as probes. We assigned 964 additional expressed sequence tags to linkage groups by restriction fragment length polymorphism analysis of a nonrecombining female backcross population. Altogether, 361.1 megabases of BAC contigs and singletons were integrated with a map containing 1,688 independent genes. A test of synteny using Oxford grid analysis with more than 500 silkworm genes revealed six versus 20 silkworm linkage groups containing eight or more orthologs of Apis versus Tribolium, respectively.

Conclusion

The integrated map contains approximately 10% of predicted silkworm genes and has an estimated 76% genome coverage by BACs. This provides a new resource for improved assembly of whole-genome shotgun data, gene annotation and positional cloning, and will serve as a platform for comparative genomics and gene discovery in Lepidoptera and other insects.

Principles of genome evolution in the Drosophila melanogaster species group

That closely related species often differ by chromosomal inversions was discovered by Sturtevant and Plunkett in 1926. Our knowledge of how these inversions originate is still very limited, although a prevailing view is that they are facilitated by ectopic recombination events between inverted repetitive sequences. The availability of genome sequences of related species now allows us to study in detail the mechanisms that generate interspecific inversions. We have analyzed the breakpoint regions of the 29 inversions that differentiate the chromosomes of Drosophila melanogaster and two closely related species, D. simulans and D. yakuba, and reconstructed the molecular events that underlie their origin. Experimental and computational analysis revealed that the breakpoint regions of 59% of the inversions (17/29) are associated with inverted duplications of genes or other nonrepetitive sequences. In only two cases do we find evidence for inverted repetitive sequences in inversion breakpoints. We propose that the presence of inverted duplications associated with inversion breakpoint regions is the result of staggered breaks, either isochromatid or chromatid, and that this, rather than ectopic exchange between inverted repetitive sequences, is the prevalent mechanism for the generation of inversions in the melanogaster species group. Outgroup analysis also revealed evidence for widespread breakpoint recycling. Lastly, we have found that expression domains in D. melanogaster may be disrupted in D. yakuba, bringing into question their potential adaptive significance.

The organization of genes on chromosomes changes over evolutionary time. In some organisms, such as fruit flies and mosquitoes, inversions of chromosome regions are widespread. This has been associated with adaptation to environmental pressures and speciation. However, the mechanisms by which inversions are generated at the molecular level are poorly understood. The prevailing view involves the interactions of sequences that are moderately repeated in the genome. Here, we use molecular and computational methods to study 29 inversions that differentiate the chromosomes of three closely related fruit fly species. We find little support for a causal role of repetitive sequences in the origin of inversions and, instead, detect the presence of inverted duplications of ancestrally unique sequences (generally protein-coding genes) in the breakpoint regions of many inversions. This leads us to propose an alternative model in which the generation of inversions is coupled with the generation of duplications of flanking sequences. Additionally, we find evidence for genomic regions that are prone to breakage, being associated with inversions generated independently during the evolution of the ancestors of existing species.

Chromosomal inversion breakpoints were compared between three closely related Drosophila species. Many are associated with inverted gene duplications, suggesting that the prevalent mechanism for their generation involves staggered breakpoints.

The tiny eukaryote Ostreococcus provides genomic insights into the paradox of plankton speciation

The smallest known eukaryotes, at approximately 1-mum diameter, are Ostreococcus tauri and related species of marine phytoplankton. The genome of Ostreococcus lucimarinus has been completed and compared with that of O. tauri. This comparison reveals surprising differences across orthologous chromosomes in the two species from highly syntenic chromosomes in most cases to chromosomes with almost no similarity. Species divergence in these phytoplankton is occurring through multiple mechanisms acting differently on different chromosomes and likely including acquisition of new genes through horizontal gene transfer. We speculate that this latter process may be involved in altering the cell-surface characteristics of each species. In addition, the genome of O. lucimarinus provides insights into the unique metal metabolism of these organisms, which are predicted to have a large number of selenocysteine-containing proteins. Selenoenzymes are more catalytically active than similar enzymes lacking selenium, and thus the cell may require less of that protein. As reported here, selenoenzymes, novel fusion proteins, and loss of some major protein families including ones associated with chromatin are likely important adaptations for achieving a small cell size.

Evolution of secretin family GPCR members in the metazoa

Background

Comparative approaches using protostome and deuterostome data have greatly contributed to understanding gene function and organismal complexity. The family 2 G-protein coupled receptors (GPCRs) are one of the largest and best studied hormone and neuropeptide receptor families. They are suggested to have arisen from a single ancestral gene via duplication events. Despite the recent identification of receptor members in protostome and early deuterostome genomes, relatively little is known about their function or origin during metazoan divergence. In this study a comprehensive description of family 2 GPCR evolution is given based on in silico and expression analyses of the invertebrate receptor genes.

Results

Family 2 GPCR members were identified in the invertebrate genomes of the nematodes C. elegans and C. briggsae, the arthropods D. melanogaster and A. gambiae (mosquito) and in the tunicate C. intestinalis. This suggests that they are of ancient origin and have evolved through gene/genome duplication events. Sequence comparisons and phylogenetic analyses have demonstrated that the immediate gene environment, with regard to gene content, is conserved between the protostome and deuterostome receptor genomic regions. Also that the protostome genes are more like the deuterostome Corticotrophin Releasing Factor (CRF) and Calcitonin/Calcitonin Gene-Related Peptide (CAL/CGRP) receptors members than the other family 2 GPCR members. The evolution of family 2 GPCRs in deuterostomes is characterised by acquisition of new family members, with SCT (Secretin) receptors only present in tetrapods. Gene structure is characterised by an increase in intron number with organismal complexity with the exception of the vertebrate CAL/CGRP receptors.

Conclusion

The family 2 GPCR members provide a good example of gene duplication events occurring in tandem with increasing organismal complexity during metazoan evolution. The putative ancestral receptors are proposed to be more like the deuterostome CAL/CGRP and CRF receptors and this may be associated with their fundamental role in calcium regulation and the stress response, both of which are essential for survival.

Testing chromosomal phylogenies and inversion breakpoint reuse in Drosophila

A combination of cytogenetic and bioinformatic procedures was used to test the chromosomal phylogeny relating Drosophila buzzatii with D. repleta. Chromosomes X and 2, harboring most of the inversions fixed between these two species, were analyzed. First, chromosomal segments conserved during the divergence of the two species were identified by comparative in situ hybridization to the D. repleta chromosomes of 180 BAC clones from a BAC-based physical map of the D. buzzatii genome. These conserved segments were precisely delimited with the aid of clones containing inversion breakpoints. Then GRIMM software was used to estimate the minimum number of rearrangements necessary to transform one genome into the other and identify all possible rearrangement scenarios. Finally, the most plausible inversion trajectory was tested by hybridizing 12 breakpoint-bearing BAC clones to the chromosomes of seven other species in the repleta group. The results show that chromosomes X and 2 of D. buzzatii and D. repleta differ by 12 paracentric inversions. Nine of them are fixed in chromosome 2 and entail two breakpoint reuses. Our results also show that the cytological relationship between D. repleta and D. mercatorum is closer than that between D. repleta and D. peninsularis, and we propose that the phylogenetic relationships in this lineage of the repleta group be reconsidered. We also estimated the rate of rearrangement between D. repleta and D. buzzatii and conclude that rates within the genus Drosophila vary substantially between lineages, even within a single species group.

Comparative polytene chromosome maps of D. montana and D. virilis

Chromosomal inversion polymorphism was characterized in Finnish Drosophila montana populations. A total of 14 polymorphic inversions were observed in Finnish D. montana of which nine had not been described before. The number of polymorphic inversions in each chromosome was not significantly different from that expected, assuming equal chance of occurrence in the euchromatic genome. There was, however, no correlation between the number of polymorphic inversions and that of fixed inversions in each chromosome. Therefore, a simple neutral model does not explain the evolutionary dynamics of inversions. Furthermore, in contrast to results obtained by others, no significant correlation was found between the two transposable elements (TEs) Penelope and Ulysses and inversion breakpoints in D. montana. This result suggests that these TEs were not involved in the creation of the polymorphic inversions seen in D. montana. A comparative analysis of D. montana and Drosophila virilis polytene chromosomes 4 and 5 was performed with D. virilis bacteriophage P1 clones, thus completing the comparative studies of the two species.

A genomic comparison of faster-sex, faster-X, and faster-male evolution between Drosophila melanogaster and Drosophila pseudoobscura

A genomic comparison of Drosophila melanogaster and Drosophila pseudoobscura provides a unique opportunity to investigate factors involved in sequence divergence. The chromosomal arrangements of these species include an autosomal segment in D. melanogaster which is homologous to part of the X chromosome in D. pseudoobscura. Using orthologues to calculate rates of nonsynonymous (d(N)) substitutions, we found genes on the X chromosome to be significantly more diverged than those on the autosomes, but it is not true for segment 3L-XR which is autosomal in D. melanogaster (3L) and X-linked in D. pseudoobscura (XR). We also found that the median d(N) values for genes having reproductive functions in either the male, the female, or both sexes are higher than those for sequences without reproductive function and even higher for sequences involved in male-specific function. These estimates of divergence for male sex-related sequences are most likely underestimates, as the very rapidly evolving reproductive genes would tend to lose homology sooner and thus not be included in the comparison of orthologues. We also noticed a high proportion of male reproductive genes among the othologous genes with the highest rates of d(N). Reproductive genes with and without an orthologue in D. pseudoobscura were compared among D. melanogaster, D. simulans, and D. yakuba and it was found that there were in fact higher rates of divergence in the group without a D. pseudoobscura orthologue. These results, from widely separated taxa, bolster the thesis that sexual system genes experience accelerated rates of change in comparison to nonsexual genes in evolution and speciation.

Rates and patterns of chromosomal evolution in Drosophila pseudoobscura and D. miranda

Comparisons of gene orders between species permit estimation of the rate of chromosomal evolution since their divergence from a common ancestor. We have compared gene orders on three chromosomes of Drosophila pseudoobscura with its close relative, D. miranda, and the distant outgroup species, D. melanogaster, by using the public genome sequences of D. pseudoobscura and D. melanogaster and approximately 50 in situ hybridizations of gene probes in D. miranda. We find no evidence for extensive transfer of genes among chromosomes in D. miranda. The rates of chromosomal rearrangements between D. miranda and D. pseudoobscura are far higher than those found before in Drosophila and approach those for nematodes, the fastest rates among higher eukaryotes. In addition, we find that the D. pseudoobscura chromosome with the highest level of inversion polymorphism (Muller's element C) does not show an unusually fast rate of evolution with respect to chromosome structure, suggesting that this classic case of inversion polymorphism reflects selection rather than mutational processes. On the basis of our results, we propose possible ancestral arrangements for the D. pseudoobscura C chromosome, which are different from those in the current literature. We also describe a new method for correcting for rearrangements that are not detected with a limited set of markers.

A BAC-based physical map of the Drosophila buzzatii genome

Large-insert genomic libraries facilitate cloning of large genomic regions, allow the construction of clone-based physical maps, and provide useful resources for sequencing entire genomes. Drosophila buzzatii is a representative species of the repleta group in the Drosophila subgenus, which is being widely used as a model in studies of genome evolution, ecological adaptation, and speciation. We constructed a Bacterial Artificial Chromosome (BAC) genomic library of D. buzzatii using the shuttle vector pTARBAC2.1. The library comprises 18,353 clones with an average insert size of 152 kb and an approximately 18x expected representation of the D. buzzatii euchromatic genome. We screened the entire library with six euchromatic gene probes and estimated the actual genome representation to be approximately 23x. In addition, we fingerprinted by restriction digestion and agarose gel electrophoresis a sample of 9555 clones, and assembled them using FingerPrint Contigs (FPC) software and manual editing into 345 contigs (mean of 26 clones per contig) and 670 singletons. Finally, we anchored 181 large contigs (containing 7788 clones) to the D. buzzatii salivary gland polytene chromosomes by in situ hybridization of 427 representative clones. The BAC library and a database with all the information regarding the high coverage BAC-based physical map described in this paper are available to the research community.

Conservation of regulatory sequences and gene expression patterns in the disintegrating Drosophila Hox gene complex

Homeotic (Hox) genes are usually clustered and arranged in the same order as they are expressed along the anteroposterior body axis of metazoans. The mechanistic explanation for this colinearity has been elusive, and it may well be that a single and universal cause does not exist. The Hox-gene complex (HOM-C) has been rearranged differently in several Drosophila species, producing a striking diversity of Hox gene organizations. We investigated the genomic and functional consequences of the two HOM-C splits present in Drosophila buzzatii. Firstly, we sequenced two regions of the D. buzzatii genome, one containing the genes labial and abdominal A, and another one including proboscipedia, and compared their organization with that of D. melanogaster and D. pseudoobscura in order to map precisely the two splits. Then, a plethora of conserved noncoding sequences, which are putative enhancers, were identified around the three Hox genes closer to the splits. The position and order of these enhancers are conserved, with minor exceptions, between the three Drosophila species. Finally, we analyzed the expression patterns of the same three genes in embryos and imaginal discs of four Drosophila species with different Hox-gene organizations. The results show that their expression patterns are conserved despite the HOM-C splits. We conclude that, in Drosophila, Hox-gene clustering is not an absolute requirement for proper function. Rather, the organization of Hox genes is modular, and their clustering seems the result of phylogenetic inertia more than functional necessity.

Duplicative and conservative transpositions of larval serum protein 1 genes in the genus Drosophila

Interspecific comparative molecular analyses of transposed genes and their flanking regions can help to elucidate the time, direction, and mechanism of gene transposition. In the Drosophila melanogaster genome, three Larval serum protein 1 (Lsp1) genes (alpha, beta and gamma) are present and each of them is located on a different chromosome, suggesting multiple transposition events. We have characterized the molecular organization of Lsp1 genes in D. buzzatii, a species of the Drosophila subgenus and in D. pseudoobscura, a species of the Sophophora subgenus. Our results show that only two Lsp1 genes (beta and gamma) exist in these two species. The same chromosomal localization and genomic organization, different from that of D. melanogaster, is found in both species for the Lsp1beta and Lsp1gamma genes. Overall, at least two duplicative and two conservative transpositions are necessary to explain the present chromosomal distribution of Lsp1 genes in the three Drosophila species. Clear evidence for implication of snRNA genes in the transposition of Lsp1beta in Drosophila has been found. We suggest that an ectopic exchange between highly similar snRNA sequences was responsible for the transposition of this gene. We have also identified the putative cis-acting regulatory regions of these genes, which seemingly transposed along with the coding sequences.

Insertion-deletion biases and the evolution of genome size

Numerous theories have been proposed to account for the pronounced differences in the quantity of non-coding DNA among eukaryotic genomes, but the current repertoire remains incomplete because the only explicit mechanisms it provides involve DNA gain. It has been proposed more recently that biases in spontaneous insertions and deletions (indels) can lead to genome shrinkage by mutational mechanisms alone. The present article provides the first detailed critical discussion of this approach, and covers three different ideas related to it: (1) the general notion of DNA loss by deletion bias, (2) the "DNA loss hypothesis" which supposes that variation in genome size can be attributed to differences in DNA loss rate, and (3) the "mutational equilibrium model" which attempts to describe the long-term evolution of genome size. The mutational equilibrium model is found to be problematic, and it is noted that DNA loss by small indels is too slow in real time to determine variation in genome size above a relatively low threshold. Some alternative explanations for the observed patterns are provided, and the critique also identifies some potential problems with the current dataset. These include a failure to cite a more detailed (and somewhat contradictory) mammalian dataset, a questionable use of arithmetic means with highly skewed data, and important discrepancies among the particular DNA sequences so far analyzed. Overall, evolutionary reductions in genome size are considered important, but the specific mechanism relating to small deletion bias is far too weak to be accepted as a primary determinant of genome size variation in general.

Hierarchical scaffolding with Bambus

The output of a genome assembler generally comprises a collection of contiguous DNA sequences (contigs) whose relative placement along the genome is not defined. A procedure called scaffolding is commonly used to order and orient these contigs using paired read information. This ordering of contigs is an essential step when finishing and analyzing the data from a whole-genome shotgun project. Most recent assemblers include a scaffolding module; however, users have little control over the scaffolding algorithm or the information produced. We thus developed a general-purpose scaffolder, called Bambus, which affords users significant flexibility in controlling the scaffolding parameters. Bambus was used recently to scaffold the low-coverage draft dog genome data. Most significantly, Bambus enables the use of linking data other than that inferred from mate-pair information. For example, the sequence of a completed genome can be used to guide the scaffolding of a related organism. We present several applications of Bambus: support for finishing, comparative genomics, analysis of the haplotype structure of genomes, and scaffolding of a mammalian genome at low coverage. Bambus is available as an open-source package from our Web site.

The Drosophila melanogaster genome

Drosophila's importance as a model organism made it an obvious choice to be among the first genomes sequenced, and the Release 1 sequence of the euchromatic portion of the genome was published in March 2000. This accomplishment demonstrated that a whole genome shotgun (WGS) strategy could produce a reliable metazoan genome sequence. Despite the attention to sequencing methods, the nucleotide sequence is just the starting point for genome-wide analyses; at a minimum, the genome sequence must be interpreted using expressed sequence tag (EST) and complementary DNA (cDNA) evidence and computational tools to identify genes and predict the structures of their RNA and protein products. The functions of these products and the manner in which their expression and activities are controlled must then be assessed-a much more challenging task with no clear endpoint that requires a wide variety of experimental and computational methods. We first review the current state of the Drosophila melanogaster genome sequence and its structural annotation and then briefly summarize some promising approaches that are being taken to achieve an initial functional annotation.

Rearrangements and chromosomal evolution

Comparisons of the genome sequences of related species suggests varying patterns of chromosomal rearrangements in different evolutionary lineages. In this review, I focus on the quantitative characterization of rearrangement processes and discuss specific inventories that have been compiled to date. Of particular interest are the statistical distribution of the lengths of inverted or locally transposed chromosome fragments (notably very short ones), inhomogeneities in susceptibility to evolutionary breakpoints in chromosomal regions, the relative importance of genome doubling in the history of multicellular eukaryotes, and of lateral transfer versus gene gain and loss in prokaryotes. These developments provide challenges to computational biologists to refine, revise and scale up mathematical models and algorithms for analyzing genome rearrangements.

Low occurrence of gene transposition events during the evolution of the genus Drosophila

The role played by gene transpositions during the evolution of eukaryotic genomes is still poorly understood and indeed has been analyzed in detail only in nematodes. In Drosophila, a limited number of transpositions have been detected by comparing the chromosomal location of genes between different species. The relative importance of gene transposition versus other types of chromosomal rearrangements, for example, inversions, has not yet been evaluated. Here, we use physical mapping to perform an extensive search for long-distance gene transpositions and assess their impact during the evolution of the Drosophila genome. We compare the relative order of 297 molecular markers that cover 60% of the euchromatic fraction of the genome between two related Drosophila species and conclude that the frequency of gene transpositions is very low, namely one order of magnitude lower than that of nematodes. In addition, gene transpositions seem to be events almost exclusively associated with genes of repetitive nature such as the Histone gene complex (HIS-C).

Inversions and gene order shuffling in Anopheles gambiae and A. funestus

In tropical Africa, Anopheles funestus is one of the three most important malaria vectors. We physically mapped 157 A. funestus complementary DNAs (cDNAs) to the polytene chromosomes of this species. Sequences of the cDNAs were mapped in silico to the A. gambiae genome as part of a comparative genomic study of synteny, gene order, and sequence conservation between A. funestus and A. gambiae. These species are in the same subgenus and diverged about as recently as humans and chimpanzees. Despite nearly perfect preservation of synteny, we found substantial shuffling of gene order along corresponding chromosome arms. Since the divergence of these species, at least 70 chromosomal inversions have been fixed, the highest rate of rearrangement of any eukaryote studied to date. The high incidence of paracentric inversions and limited colinearity suggests that locating genes in one anopheline species based on gene order in another may be limited to closely related taxa.

Assessing the impact of comparative genomic sequence data on the functional annotation of the Drosophila genome

BACKGROUND

It is widely accepted that comparative sequence data can aid the functional annotation of genome sequences; however, the most informative species and features of genome evolution for comparison remain to be determined.

RESULTS

We analyzed conservation in eight genomic regions (apterous, even-skipped, fushi tarazu, twist, and Rhodopsins 1, 2, 3 and 4) from four Drosophila species (D. erecta, D. pseudoobscura, D. willistoni, and D. littoralis) covering more than 500 kb of the D. melanogaster genome. All D. melanogaster genes (and 78-82% of coding exons) identified in divergent species such as D. pseudoobscura show evidence of functional constraint. Addition of a third species can reveal functional constraint in otherwise non-significant pairwise exon comparisons. Microsynteny is largely conserved, with rearrangement breakpoints, novel transposable element insertions, and gene transpositions occurring in similar numbers. Rates of amino-acid substitution are higher in uncharacterized genes relative to genes that have previously been studied. Conserved non-coding sequences (CNCSs) tend to be spatially clustered with conserved spacing between CNCSs, and clusters of CNCSs can be used to predict enhancer sequences.

CONCLUSIONS

Our results provide the basis for choosing species whose genome sequences would be most useful in aiding the functional annotation of coding and cis-regulatory sequences in Drosophila. Furthermore, this work shows how decoding the spatial organization of conserved sequences, such as the clustering of CNCSs, can complement efforts to annotate eukaryotic genomes on the basis of sequence conservation alone.