P1 clones from Drosophila melanogaster as markers to study the chromosomal evolution of Muller's A element in two species of the obscura group of Drosophila

Dynamic turnover of centromeres drives karyotype evolution in Drosophila

Centromeres are the basic unit for chromosome inheritance, but their evolutionary dynamics is poorly understood. We generate high-quality reference genomes for multiple Drosophila obscura group species to reconstruct karyotype evolution. All chromosomes in this lineage were ancestrally telocentric and the creation of metacentric chromosomes in some species was driven by de novo seeding of new centromeres at ancestrally gene-rich regions, independently of chromosomal rearrangements. The emergence of centromeres resulted in a drastic size increase due to repeat accumulation, and dozens of genes previously located in euchromatin are now embedded in pericentromeric heterochromatin. Metacentric chromosomes secondarily became telocentric in the pseudoobscura subgroup through centromere repositioning and a pericentric inversion. The former (peri)centric sequences left behind shrunk dramatically in size after their inactivation, yet contain remnants of their evolutionary past, including increased repeat-content and heterochromatic environment. Centromere movements are accompanied by rapid turnover of the major satellite DNA detected in (peri)centromeric regions.

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.

The High-Quality Genome Sequence of the Oceanic Island Endemic Species Drosophila guanche Reveals Signals of Adaptive Evolution in Genes Related to Flight and Genome Stability

Drosophila guanche is a member of the obscura group that originated in the Canary Islands archipelago upon its colonization by D. subobscura. It evolved into a new species in the laurisilva, a laurel forest present in wet regions that in the islands have only minor long-term weather fluctuations. Oceanic island endemic species such as D. guanche can become model species to investigate not only the relative role of drift and adaptation in speciation processes but also how population size affects nucleotide variation. Moreover, the previous identification of two satellite DNAs in D. guanche makes this species attractive for studying how centromeric DNA evolves. As a prerequisite for its establishment as a model species suitable to address all these questions, we generated a high-quality D. guanche genome sequence composed of 42 cytologically mapped scaffolds, which are assembled into six super-scaffolds (one per chromosome). The comparative analysis of the D. guanche proteome with that of twelve other Drosophila species identified 151 genes that were subject to adaptive evolution in the D. guanche lineage, with a subset of them being involved in flight and genome stability. For example, the Centromere Identifier (CID) protein, directly interacting with centromeric satellite DNA, shows signals of adaptation in this species. Both genomic analyses and FISH of the two satellites would support an ongoing replacement of centromeric satellite DNA in D. guanche.

Partial-arm translocations in evolution of malaria mosquitoes revealed by high-coverage physical mapping of the Anopheles atroparvus genome

Background

Malaria mosquitoes have had a remarkable stability in the number of chromosomes in their karyotype (2n = 6) during 100 million years of evolution. Moreover, autosomal arms were assumed to maintain their integrity even if their associations with each other changed via whole-arm translocations. Here we use high-coverage comparative physical genome mapping of three Anopheles species to test the extent of evolutionary conservation of chromosomal arms in malaria mosquitoes.

Results

In this study, we developed a physical genome map for Anopheles atroparvus, one of the dominant malaria vectors in Europe. Using fluorescence in situ hybridization (FISH) of DNA probes with the ovarian nurse cell polytene chromosomes and synteny comparison, we anchored 56 genomic scaffolds to the An. atroparvus chromosomes. The obtained physical map represents 89.6% of the An. atroparvus genome. This genome has the second highest mapping coverage among Anophelinae assemblies after An. albimanus, which has 98.2% of the genome assigned to its chromosomes. A comparison of the An. atroparvus, An. albimanus, and An. gambiae genomes identified partial-arm translocations between the autosomal arms that break down the integrity of chromosome elements in evolution affecting the structure of the genetic material in the pericentromeric regions. Unlike An. atroparvus and An. albimanus, all chromosome elements of An. gambiae are fully syntenic with chromosome elements of the putative ancestral Anopheles karyotype. We also detected nonrandom distribution of large conserved synteny blocks and confirmed a higher rate of inversion fixation in the X chromosome compared with autosomes.

Conclusions

Our study demonstrates the power of physical mapping for understanding the genome evolution in malaria mosquitoes. The results indicate that syntenic relationships among chromosome elements of Anopheles species have not been fully preserved because of multiple partial-arm translocations.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4663-4) contains supplementary material, which is available to authorized users.

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.

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.

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.

Polytene chromosomal maps of 11 Drosophila species: the order of genomic scaffolds inferred from genetic and physical maps

The sequencing of the 12 genomes of members of the genus Drosophila was taken as an opportunity to reevaluate the genetic and physical maps for 11 of the species, in part to aid in the mapping of assembled scaffolds. Here, we present an overview of the importance of cytogenetic maps to Drosophila biology and to the concepts of chromosomal evolution. Physical and genetic markers were used to anchor the genome assembly scaffolds to the polytene chromosomal maps for each species. In addition, a computational approach was used to anchor smaller scaffolds on the basis of the analysis of syntenic blocks. We present the chromosomal map data from each of the 11 sequenced non-Drosophila melanogaster species as a series of sections. Each section reviews the history of the polytene chromosome maps for each species, presents the new polytene chromosome maps, and anchors the genomic scaffolds to the cytological maps using genetic and physical markers. The mapping data agree with Muller's idea that the majority of Drosophila genes are syntenic. Despite the conservation of genes within homologous chromosome arms across species, the karyotypes of these species have changed through the fusion of chromosomal arms followed by subsequent rearrangement events.

High incidence of interchromosomal transpositions in the evolutionary history of a subset of or genes in Drosophila

In insects, the odorant receptor (Or) multigene family is an intermediate-sized family with genes present in all chromosomes, indicating that duplication followed by interchromosomal transposition played an important role in the early stages of the family evolution. Here, we have explored the occurrence of interchromosomal transpositions in more recent stages through the comparative analysis of a subset of Or genes in Drosophila, where the gene content of chromosomal arms is highly conserved. The studied subset consisted of 11 Or genes located on the left arm of chromosome 3 (Muller's D element) in D. melanogaster. Our study focused on the number and chromosomal arm location of these members of the family across the 12 Drosophila species with complete genome sequences. In contrast to previous results from in situ hybridization comparative mapping that were mainly based on single-copy genes, our study, based on members of a multigene family of moderate size, revealed repeated interchromosomal transposition events and a complex history of some of the studied genes.

Divergence between the Drosophila pseudoobscura and D. persimilis genome sequences in relation to chromosomal inversions

As whole-genome sequence assemblies accumulate, a challenge is to determine how these can be used to address fundamental evolutionary questions, such as inferring the process of speciation. Here, we use the sequence assemblies of Drosophila pseudoobscura and D. persimilis to test hypotheses regarding divergence with gene flow. We observe low differentiation between the two genome sequences in pericentromeric and peritelomeric regions. We interpret this result as primarily a remnant of the correlation between levels of variation and local recombination rate observed within populations. However, we also observe lower differentiation far from the fixed chromosomal inversions distinguishing these species and greater differentiation within and near these inversions. This finding is consistent with models suggesting that chromosomal inversions facilitate species divergence despite interspecies gene flow. We also document heterogeneity among the inverted regions in their degree of differentiation, suggesting temporal differences in the origin of each inverted region consistent with the inversions arising during a process of divergence with gene flow. While this study provides insights into the speciation process using two single-genome sequences, it was informed by lower throughput but more rigorous examinations of polymorphism and divergence. This reliance highlights the need for complementary genomic and population genetic approaches for tackling fundamental evolutionary questions such as speciation.

Selection, recombination and demographic history in Drosophila miranda

Selection, recombination, and the demographic history of a species can all have profound effects on genomewide patterns of variability. To assess the impact of these forces in the genome of Drosophila miranda, we examine polymorphism and divergence patterns at 62 loci scattered across the genome. In accordance with recent findings in D. melanogaster, we find that noncoding DNA generally evolves more slowly than synonymous sites, that the distribution of polymorphism frequencies in noncoding DNA is significantly skewed toward rare variants relative to synonymous sites, and that long introns evolve significantly slower than short introns or synonymous sites. These observations suggest that most noncoding DNA is functionally constrained and evolving under purifying selection. However, in contrast to findings in the D. melanogaster species group, we find little evidence of adaptive evolution acting on either coding or noncoding sequences in D. miranda. Levels of linkage disequilibrium (LD) in D. miranda are comparable to those observed in D. melanogaster, but vary considerably among chromosomes. These patterns suggest a significantly lower rate of recombination on autosomes, possibly due to the presence of polymorphic autosomal inversions and/or differences in chromosome sizes. All chromosomes show significant departures from the standard neutral model, including too much heterogeneity in synonymous site polymorphism relative to divergence among loci and a general excess of rare synonymous polymorphisms. These departures from neutral equilibrium expectations are discussed in the context of nonequilibrium models of demography and selection.

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.

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.

Comparative genome sequencing of Drosophila pseudoobscura: chromosomal, gene, and cis-element evolution

We have sequenced the genome of a second Drosophila species, Drosophila pseudoobscura, and compared this to the genome sequence of Drosophila melanogaster, a primary model organism. Throughout evolution the vast majority of Drosophila genes have remained on the same chromosome arm, but within each arm gene order has been extensively reshuffled, leading to a minimum of 921 syntenic blocks shared between the species. A repetitive sequence is found in the D. pseudoobscura genome at many junctions between adjacent syntenic blocks. Analysis of this novel repetitive element family suggests that recombination between offset elements may have given rise to many paracentric inversions, thereby contributing to the shuffling of gene order in the D. pseudoobscura lineage. Based on sequence similarity and synteny, 10,516 putative orthologs have been identified as a core gene set conserved over 25-55 million years (Myr) since the pseudoobscura/melanogaster divergence. Genes expressed in the testes had higher amino acid sequence divergence than the genome-wide average, consistent with the rapid evolution of sex-specific proteins. Cis-regulatory sequences are more conserved than random and nearby sequences between the species--but the difference is slight, suggesting that the evolution of cis-regulatory elements is flexible. Overall, a pattern of repeat-mediated chromosomal rearrangement, and high coadaptation of both male genes and cis-regulatory sequences emerges as important themes of genome divergence between these species of Drosophila.

A comparative genomic analysis of two distant diptera, the fruit fly, Drosophila melanogaster, and the malaria mosquito, Anopheles gambiae

Genome evolution entails changes in the DNA sequence of genes and intergenic regions, changes in gene numbers, and also changes in gene order along the chromosomes. Genes are reshuffled by chromosomal rearrangements such as deletions/insertions, inversions, translocations, and transpositions. Here we report a comparative study of genome organization in the main African malaria vector, Anopheles gambiae, relative to the recently determined sequence of the Drosophila melanogaster genome. The ancestral lines of these two dipteran insects are thought to have separated approximately 250 Myr, a long period that makes this genome comparison especially interesting. Sequence comparisons have identified 113 pairs of putative orthologs of the two species. Chromosomal mapping of orthologous genes reveals that each polytene chromosome arm has a homolog in the other species. Between 41% and 73% of the known orthologous genes remain linked in the respective homologous chromosomal arms, with the remainder translocated to various nonhomologous arms. Within homologous arms, gene order is extensively reshuffled, but a limited degree of conserved local synteny (microsynteny) can be recognized.

Changes in the recombinational environment affect divergence in the yellow gene of Drosophila

The complete coding region of the yellow (y) gene was sequenced in different Drosophila species. In the species of the melanogaster subgroup (D. melanogaster, D. simulans, D. mauritiana, D. yakuba, and D. erecta), this gene is located at the tip of the X chromosome in a region with a strong reduction in recombination rate. In contrast, in D. ananassae (included in the ananassae subgroup of the melanogaster group) and in the obscura group species (D. subobscura, D. madeirensis, D. guanche, and D. pseudoobscura), the y gene is located in regions with normal recombination rates. As predicted by the hitchhiking and background selection models, this change in the recombinational environment affected synonymous divergence in the y-gene-coding region. Estimates of the number of synonymous substitutions per site were much lower between the obscura group species and D. ananassae than between the species of the obscura group and the melanogaster subgroup. In fact, a highly significant increase in the rate of synonymous substitution was detected in all lineages leading to the species of the melanogaster subgroup relative to the D. ananassae lineage. This increase can be explained by a higher fixation rate of mutations from preferred to unpreferred codons (slightly deleterious mutations). The lower codon bias detected in all species of the melanogaster subgroup relative to D. ananassae (or to the obscura group species) would be consistent with this proposal. Therefore, at least in Drosophila, changes in the recombination rate in different lineages might cause deviations of the molecular-clock hypothesis and contribute to the overdispersion of the rate of synonymous substitution. In contrast, the change in the recombinational environment of the y gene has no detectable effect on the rate of amino acid replacement in the Yellow protein.

Molecular organization of the Drosophila melanogaster Adh chromosomal region in D. repleta and D. buzzatii, two distantly related species of the Drosophila subgenus

The molecular organization of a 1.944-Mb chromosomal region of Drosophila melanogaster around the Adh locus has been analyzed in two repleta group species: D. repleta and D. buzzatii. The extensive genetic and molecular information about this region in D. melanogaster makes it a prime choice for comparative studies of genomic organization among distantly related species. A set of 26 P1 phages from D. melanogaster were successfully hybridized using fluorescence in-situ hybridization (FISH) to the salivary gland chromosomes of both repleta group species. The results show that the Adh region is distributed in D. repleta and D. buzatii over six distant sites of chromosome 3, homologous to chromosomal arm 2L of D. melanogaster (Muller's element B). This observation implies a density of 2.57 fixed breakpoints per Mb in the Adh region and suggests a considerable reorganization of this chromosomal element via the fixation of paracentric inversions. Nevertheless, breakpoint density in the Adh region is three times lower than that estimated for D. repleta chromosome 2, homologous to D. melanogaster 3R (Muller's element E). Differences in the rate of evolution among chromosomal elements are seemingly persistent in the Drosophila genus over long phylogenetic distances.

How malleable is the eukaryotic genome? Extreme rate of chromosomal rearrangement in the genus Drosophila

During the evolution of the genus Drosophila, the molecular organization of the major chromosomal elements has been repeatedly rearranged via the fixation of paracentric inversions. Little detailed information is available, however, on the extent and effect of these changes at the molecular level. In principle, a full description of the rate and pattern of change could reveal the limits, if any, to which the eukaryotic genome can accommodate reorganizations. We have constructed a high-density physical map of the largest chromosomal element in Drosophila repleta (chromosome 2) and compared the order and distances between the markers with those on the homologous chromosomal element (3R) in Drosophila melanogaster. The two species belong to different subgenera (Drosophila and Sophophora, respectively), which diverged 40-62 million years (Myr) ago and represent, thus, the farthest lineages within the Drosophila genus. The comparison reveals extensive reshuffling of gene order from centromere to telomere. Using a maximum likelihood method, we estimate that 114 +/- 14 paracentric inversions have been fixed in this chromosomal element since the divergence of the two species, that is, 0.9-1.4 inversions fixed per Myr. Comparison with available rates of chromosomal evolution, taking into account genome size, indicates that the Drosophila genome shows the highest rate found so far in any eukaryote. Twenty-one small segments (23-599 kb) comprising at least two independent (nonoverlapping) markers appear to be conserved between D. melanogaster and D. repleta. These results are consistent with the random breakage model and do not provide significant evidence of functional constraint of any kind. They support the notion that the Drosophila genome is extraordinarily malleable and has a modular organization. The high rate of chromosomal change also suggests a very limited transferability of the positional information from the Drosophila genome to other insects.

Chromosomal homologies between Drosophila melanogaster and D. funebris determined by in-situ hybridization

Seventeen biotin-labeled DNA sequences were hybridized to polytene chromosomes of Drosophila melanogaster and D. funebris in order to establish chromosomal homologies between these species. Ten probes correspond to cloned DNA sequences from D. melanogaster (RpII 215, MHC, H3-H4, Tor, hsp 68, hsp 28/23, hsp 83, PP1alpha, RpII 140, and ey), four are clones isolated from a D. subobscura genomic library (Xdh, lambdaDsubS3, lambdaDsubG3, and lambdaDsubG4), two are clones from D. funebris (F2 and Adh) and one from D. virilis (ci). The probes were chosen in order to cover all the autosomes, since X-chromosome homologies have already been studied by linkage analysis of morphological mutants. Most probes gave a unique hybridization signal; consequently, our results allow unambiguous inferences about chromosomal homologies. The results show extensive gene rearrangement within all chromosomal elements, probably due to paracentric inversions, but are consistent with Muller's proposal that chromosomal elements have conserved their genetic content during the evolution of Drosophila.

Discordant rates of chromosome evolution in the Drosophila virilis species group

In Drosophila, the availability of polytene chromosome maps and of sets of probes covering most regions of the chromosomes allows a direct comparison of the organization of the genome in different species. In this work, we report the localization, in Drosophila virilis, D. montana, and D. novamexicana, of > 100 bacteriophage PI clones containing approximately 65 kilobase inserts of genomic DNA from D. virilis. Each clone hybridizes with a single euchromatic site in either chromosome 1 or chromosome 3 in D. virilis. From these data, it is possible to estimate the minimum number of inversions required to transform the map positions of the probes in one species into the map positions of the same probes in a related species. The data indicate that, in the D. virilis species group, the X chromosome has up to four times the number of inversions as are observed in chromosome 3. The first photographic polytene chromosome maps for D. montana and D. novamexicana are also presented.

Differentiation of Muller's chromosomal elements D and E in the obscura group of Drosophila

Twenty-two markers located on Muller's elements D or E have been mapped by in situ hybridization in six species of the obscura group of Drosophila and in D. melanogaster. The obscura species can be grouped into a Palearctic cluster (D. subobscura, D. madeirensis and D. guanche) and a Nearctic one (D. pseudoobscura, D. persimilis and D. miranda). Eleven of the probes contain known genes: E74, Acp70A, Est5, hsp28/23, hsp83, emc, hsp70, Xdh, Acph-1, Cec and rp49. The remaining probes are recombinant phages isolated from a D. subobscura genomic library. All these markers hybridize to the putative homologous chromosome or chromosomal arm of elements D and E. Thus, these elements have conserved their genic content during species divergence. Chromosomal homologies proposed previously for each element among the species of the same cluster have been compared with the present results. The distribution of markers within each element has changed considerably as inferred from pairwise comparisons of obscura species included in the two different clusters. Only chromosomal segments defined by closely linked markers have been conserved: one such segment has been detected in element D and three in element E between D subobscura and D. pseudoobscura.