The salivary gland chromosomes of Drosophila virilis: a cytological map, pattern of transcription and aspects of chromosome evolution

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.

Heterochromatic genes in Drosophila: a comparative analysis of two genes

Centromeric heterochromatin comprises approximately 30% of the Drosophila melanogaster genome, forming a transcriptionally repressive environment that silences euchromatic genes juxtaposed nearby. Surprisingly, there are genes naturally resident in heterochromatin, which appear to require this environment for optimal activity. Here we report an evolutionary analysis of two genes, Dbp80 and RpL15, which are adjacent in proximal 3L heterochromatin of D. melanogaster. DmDbp80 is typical of previously described heterochromatic genes: large, with repetitive sequences in its many introns. In contrast, DmRpL15 is uncharacteristically small. The orthologs of these genes were examined in D. pseudoobscura and D. virilis. In situ hybridization and whole-genome assembly analysis show that these genes are adjacent, but not centromeric in the genome of D. pseudoobscura, while they are located on different chromosomal elements in D. virilis. Dbp80 gene organization differs dramatically among these species, while RpL15 structure is conserved. A bioinformatic analysis in five additional Drosophila species demonstrates active repositioning of these genes both within and between chromosomal elements. This study shows that Dbp80 and RpL15 can function in contrasting chromatin contexts on an evolutionary timescale. The complex history of these genes also provides unique insight into the dynamic nature of genome evolution.

A Kunitz type protease inhibitor related protein is synthesized in Drosophila prepupal salivary glands and released into the moulting fluid during pupation

From the Drosophila virilis late puff region 31C, we microcloned two neighbouring genes, Kil-1 and Kil-2, that encode putative Kunitz serine protease inhibitor like proteins. The Kil-1 gene is expressed exclusively in prepupal salivary glands. Using a size mutant of the KIL-1 protein and MALDI-TOF analysis, we demonstrate that during pupation this protein is released from the prepupal salivary glands into the pupation fluid covering the surface of the pupa. 3-D-structure predictions are consistent with the known crystal structure of the human Kunitz type protease inhibitor 2KNT. This is the first experimental proof for the extracorporal presence of a distinct Drosophila prepupal salivary gland protein. Possible functions of KIL-1 in the context of the control of proteolytic activities in the pupation fluid are discussed.

Polytene Chromosomes: 70 Years of Genetic Research

Polytene chromosomes were described in 1881 and since 1934 they have served as an outstanding model for a variety of genetic experiments. Using the polytene chromosomes, numerous biological phenomena were discovered. First the polytene chromosomes served as a model of the interphase chromosomes in general. In polytene chromosomes, condensed (bands), decondensed (interbands), genetically active (puffs), and silent (pericentric and intercalary heterochromatin as well as regions subject to position effect variegation) regions were found and their features were described in detail. Analysis of the general organization of replication and transcription at the cytological level has become possible using polytene chromosomes. In studies of sequential puff formation it was found for the first time that the steroid hormone (ecdysone) exerts its action through gene activation, and that the process of gene activation upon ecdysone proceeds as a cascade. Namely on the polytene chromosomes a new phenomenon of cellular stress response (heat shock) was discovered. Subsequently chromatin boundaries (insulators) were discovered to flank the heat shock puffs. Major progress in solving the problems of dosage compensation and position effect variegation phenomena was mainly related to studies on polytene chromosomes. This review summarizes the current status of studies of polytene chromosomes and of various phenomena described using this successful model.

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).

Population, evolutionary and genomic consequences of interference selection

Weakly selected mutations are most likely to be physically clustered across genomes and, when sufficiently linked, they alter each others' fixation probability, a process we call interference selection (IS). Here we study population genetics and evolutionary consequences of IS on the selected mutations themselves and on adjacent selectively neutral variation. We show that IS reduces levels of polymorphism and increases low-frequency variants and linkage disequilibrium, in both selected and adjacent neutral mutations. IS can account for several well-documented patterns of variation and composition in genomic regions with low rates of crossing over in Drosophila. IS cannot be described simply as a reduction in the efficacy of selection and effective population size in standard models of selection and drift. Rather, IS can be better understood with models that incorporate a constant "traffic" of competing alleles. Our simulations also allow us to make genome-wide predictions that are specific to IS. We show that IS will be more severe at sites in the center of a region containing weakly selected mutations than at sites located close to the edge of the region. Drosophila melanogaster genomic data strongly support this prediction, with genes without introns showing significantly reduced codon bias in the center of coding regions. As expected, if introns relieve IS, genes with centrally located introns do not show reduced codon bias in the center of the coding region. We also show that reasonably small differences in the length of intermediate "neutral" sequences embedded in a region under selection increase the effectiveness of selection on the adjacent selected sequences. Hence, the presence and length of sequences such as introns or intergenic regions can be a trait subject to selection in recombining genomes. In support of this prediction, intron presence is positively correlated with a gene's codon bias in D. melanogaster. Finally, the study of temporal dynamics of IS after a change of recombination rate shows that nonequilibrium codon usage may be the norm rather than the exception.

Conservation of gene order, structure and sequence between three closely linked genes in Drosophila melanogaster and Drosophila virilis

In Drosophila melanogaster, the apparently unrelated genes anon-66Da, RpL14, and anon-66Db (from telomere to centromere) are located on a 5547 bp genomic fragment on chromosome arm 3L at cytological position 66D8. The three genes are tightly linked, and flanked by two relatively large genes with unknown function. We have taken a comparative genomic approach to investigate the evolutionary history of the three genes. To this end we isolated a Drosophila virilis 7.3 kb genomic fragment which is homologous to a 5.5 kb genomic region of D. melanogaster. Both fragments map to Muller's element D, namely to section 66D in D. melanogaster and to section 32E in D. virilis, and harbor the genes anon-66Da, RpL14, and anon-66Db. We demonstrate that the three genes exhibit a high conservation of gene topography in general and in detail. While most introns and intergenic regions reveal sequence divergences, there are, however, a number of interspersed conserved sequence motifs. In particular, two introns of the RpL14 gene contain a short, highly conserved 60 nt long sequence located at corresponding positions. This sequence represents a novel Drosophila small nucleolar RNA, which is homologous to human U49. Whereas DNA flanking the three genes shows no significant interspecies homologies, the 3'-flanking region in D. virilis contains sequences from the transposable element Penelope. The Penelope family of transposable elements has been shown to promote chromosomal rearrangements in the D. virilis species group. The presence of Penelope sequences in the D. virilis 7.3 kb genomic fragment may be indicative for a transposon-induced event of transposition which did not yet scramble the order of the three genes but led to the breakdown of sequence identity of the flanking DNA.

Evolution of 5S rRNA gene families in Drosophila

In Drosophila virilis, the three clusters of 5S rRNA genes on chromosome 5 comprise two different gene families (B and C), which differ profoundly in the organization of their spacer sequences. While C-type genes, which are found in two of the clusters, exhibit a true repetitive character, the B-type genes of the third cluster are each embedded in completely different genomic environments. Southern blots of genomic DNA of different D. virilis subspecies, D. hydei and D. melanogaster probed with 5S rRNA gene spacer and coding sequences demonstrate the specificity of C-type sequences for the D. virilis species group. The comparative analysis of flanking sequences of 5S rRNA genes of D. virilis, members of the D. melanogaster species subgroup and of the blowfly Calliphora erythrocephala reveals the existence of conserved sequence motifs both in the 5' upstream and 3' downstream flanking regions. Their possible roles in the control of expression and processing of the 5S rRNA precursor molecule are discussed.

Identification of X chromosomal restriction fragment length polymorphism markers and their use in a gene localization study in Drosophila virilis and D. littoralis

We have identified six restriction fragment length polymorphism (RFLP) markers based on unique gene sequences on the X chromosome of Drosophila virilis and D. littoralis. The markers were localized by in situ hybridization on larval polytene chromosomes, and the conjugation of the X chromosomes of the two species was studied in salivary glands of interspecific hybrid female larvae. The gene arrangement of D. virilis and D. littoralis appeared to be very different at the proximal end of the X chromosome preventing recombination between RFLP markers located in this area. Simple quantitative trait loci (QTL) analysis showed that five of our marker genes (including nonA and Dmca1A, previously found to affect male courtship song in D. melanogaster) are linked with a gene(s) having a major effect on species differences in the male courtship song between D. virilis and D. littoralis. This shows that the song gene(s) may be located inside a large X-chromosomal inversion in D. littoralis (as previously suggested), but that it may also be located on an area between this inversion and the centromere, close to nonA and Dmca1A. Localization of this gene or gene complex will be continued with the aid of our newly identified RFLP markers by making interspecific crosses between D. virilis group species with more similar X chromosomes.

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.

A selective sweep associated with a recent gene transposition in Drosophila miranda

In Drosophila miranda, a chromosome fusion between the Y chromosome and the autosome corresponding to Muller's element C has created a new sex chromosome system. The chromosome attached to the ancestral Y chromosome is transmitted paternally and hence is not exposed to crossing over. This chromosome, conventionally called the neo-Y, and the homologous neo-X chromosome display many properties of evolving sex chromosomes. We report here the transposition of the exuperantia1 (exu1) locus from a neo-sex chromosome to the ancestral X chromosome of D. miranda. Exu1 is known to have several critical developmental functions, including a male-specific role in spermatogenesis. The ancestral location of exu1 is conserved in the sibling species of D. miranda, as well as in a more distantly related species. The transposition of exu1 can be interpreted as an adaptive fixation, driven by a selective advantage conferred by its effect on dosage compensation. This explanation is supported by the pattern of within-species sequence variation at exu1 and the nearby exu2 locus. The implications of this phenomenon for genome evolution are discussed.

Toward a physical map of Drosophila buzzatii. Use of randomly amplified polymorphic dna polymorphisms and sequence-tagged site landmarks

We present a physical map based on RAPD polymorphic fragments and sequence-tagged sites (STSs) for the repleta group species Drosophila buzzatii. One hundred forty-four RAPD markers have been used as probes for in situ hybridization to the polytene chromosomes, and positive results allowing the precise localization of 108 RAPDs were obtained. Of these, 73 behave as effectively unique markers for physical map construction, and in 9 additional cases the probes gave two hybridization signals, each on a different chromosome. Most markers (68%) are located on chromosomes 2 and 4, which partially agree with previous estimates on the distribution of genetic variation over chromosomes. One RAPD maps close to the proximal breakpoint of inversion 2z(3) but is not included within the inverted fragment. However, it was possible to conclude from this RAPD that the distal breakpoint of 2z(3) had previously been wrongly assigned. A total of 39 cytologically mapped RAPDs were converted to STSs and yielded an aggregate sequence of 28,431 bp. Thirty-six RAPDs (25%) did not produce any detectable hybridization signal, and we obtained the DNA sequence from three of them. Further prospects toward obtaining a more developed genetic map than the one currently available for D. buzzatii are discussed.

X chromosome DNA variation in Drosophila virilis

Comparisons of polymorphism patterns between distantly related species are essential in order to determine their generality. However, most work on the genus Drosophila has been done only with species of the subgenus Sophophora. In the present work, we have sequenced one intron and surrounding coding sequences of 6 X-linked genes (chorion protein s36, elav, fused, runt, suppressor of sable and zeste) from 21 strains of wild-type Drosophila virilis (subgenus Drosophila). From these data, we have estimated the average level of DNA polymorphism, inferred the effective population size and population structure of this species, and compared the results with those obtained for other Drosophila species. There is no reduction in variation at two loci close to the centromeric heterochromatin, in contrast to Drosophila melanogaster.

Natural selection on synonymous sites is correlated with gene length and recombination in Drosophila

Evolutionary analysis of codon bias in Drosophila indicates that synonymous mutations are not neutral, but rather are subject to weak selection at the translation level. Here we show that the effectiveness of natural selection on synonymous sites is strongly correlated with the rate of recombination, in accord with the nearly neutral hypothesis. This correlation, however, is apparent only in genes encoding short proteins. Long coding regions have both a lower codon bias and higher synonymous substitution rates, suggesting that they are affected less efficiently by selection. Therefore, both the length of the coding region and the recombination rate modulate codon bias. In addition, the data indicate that selection coefficients for synonymous mutations must vary by a minimum of one or two orders of magnitude. Two hypotheses are proposed to explain the relationship among the coding region length, the codon bias, and the synonymous divergence and polymorphism levels across the range of recombination rates in Drosophila. The first hypothesis is that selection coefficients on synonymous mutations are inversely related to the total length of the coding region. The second hypothesis proposes that interference among synonymous mutations reduces the efficacy of selection on these mutations. We investigated this second hypothesis by carrying out forward simulations of weakly selected mutations in model populations. These simulations show that even with realistic recombination rates, this interference, which we call the "small-scale" Hill-Robertson effect, can have a moderately strong influence on codon bias.

Sequence, molecular organization and products of the Drosophila virilis homologs of the D. melanogaster nested genes lethal(2) tumorous imaginal discs [1(2)tid] and lethal(2) neighbour of tid [1(2)not]

In this study, we describe the isolation of the Drosophila virilis (Dvir) 6201-bp genomic fragment homologous to a 7047-bp genomic region of D. melanogaster (Dmel) that harbors the nested genes lethal(2) tumorous imaginal discs (l(2)tid), lethal(2) neighbour of tid (l(2)not) and lethal(2) relative of tid (l(2)rot). The isolated fragment, which maps at the cytogenetic position 50A5 on chromosome 5, carries the Dvir homologs of the Dmel genes l(2)tid and l(2)not. In both cases, the interspecific comparison of the determined sequences reveals a high homology regarding the protein coding regions and a high degree of evolutionary divergence concerning the intronic parts of the genes. In the two distantly related species, the particular gene within gene arrangement of the two genes is conserved, namely, Dvir tid is located in the intron of Dvir not, on the non-coding DNA strand. Interestingly, the Dvir homolog of the Dmel l(2)rot gene residing in the l(2)not intron on its coding strand, opposite l(2)tid, is not present in the 6201-bp genomic fragment. The protein predicted from the Dvir tid sequence, Dvir Tid58, exhibits 76.5% identity with the putative Tid56 protein of Dmel. The putative Dvir Not58 protein shows 71% identity with its Dmel homolog Not56. The developmental transcript and protein patterns, as well as the characteristics of the protein products encoded by the genes Dvir tid and Dvir not are similar to those identified for their Dmel homologs.

Chromosomal homology and molecular organization of Muller's elements D and E in the Drosophila repleta species group

Thirty-three DNA clones containing protein-coding genes have been used for in situ hybridization to the polytene chromosomes of two Drosophila repleta group species, D. repleta and D. buzzatii. Twenty-six clones gave positive results allowing the precise localization of 26 genes and the tentative identification of another nine. The results were fully consistent with the currently accepted chromosomal homologies and in no case was evidence for reciprocal translocations or pericentric inversions found. Most of the genes mapped to chromosomes 2 and 4 that are homologous, respectively, to chromosome arms 3R and 3L of D. melanogaster (Muller's elements E and D). The comparison of the molecular organization of-these two elements between D. melanogaster and D. repleta (two species that belong to different subgenera and diverged some 62 million years ago) showed an extensive reorganization via paracentric inversions. Using a maximum likelihood procedure, we estimated that 130 paracentric inversions have become fixed in element E after the divergence of the two lineages. Therefore, the evolution rate for element E is approximately one inversion per million years. This value is comparable to previous estimates of the rate of evolution of chromosome X and yields an estimate of 4.5 inversions per million years for the whole Drosophila genome.