THE GENETIC BASIS OF EVOLUTION OF THE MALE COURTSHIP SOUNDS IN THE DROSOPHILA VIRILIS GROUP

When courting, males of the Drosophila virilis group vibrate their wings and emit species-specific courtship sounds consisting of trains of polycyclic sound pulses. To analyze the genetic basis of evolutionary changes in the sounds we made an F₁ diallel set of reciprocal crosses between the members of the virilis phylad of the group (two stocks of D. virilis and one of D. americana americana, D. a. texana, D. novamexicana, and D. lummei). We also crossed the D. virilis stocks with the members of the montana phylad of the same group (D. kanekoi, D. littoralis, D. borealis, D. flavomontana, D. lacicola, and D. montana) and made a backcross (D. virilis x D. littoralis) x D. virilis using a D. virilis marker stock (b; sv t tb gp; cd; pe). The sounds of the hybrids were analyzed using the following parameters: the length of a pulse train (PTL), the number of pulses in a train (PN), the interpulse interval (IPI), the length of a pulse (PL), the number of cycles in a pulse (CN), and the length of a cycle (CL). In the virilis phylad, the differences between species appeared to be determined mainly by autosomal genes in each sound trait. The heritabilities (narrow-/broad-sense) obtained from the diallel tables were the following: PTL 0.662/0.817, PN 0.651/0.841, IPI 0.193/0.546, PL 0.408/0.552, CN 0.425/0.719, and CL 0.361/0.764. The direction of dominance is for longer PTL, higher PN and CN, and shorter IPI and CL. PL shows ambidirectional dominance. In the sounds of the virilis phylad species, PTL and PL seem to be phenotypically the most important parameters, since their components (PN and IPI for PTL, CN and CL for PL) are negatively correlated. In crosses between D. virilis and D. littoralis or D. flavomontana reciprocal hybrids differed from each other in PTL, IPI, PL, and CN indicating X-chromosomal or cytoplasmic inheritance. In the backcrosses between D. virilis and D. littoralis the role of the X chromosome was ascertained to be decisive. We conclude that an X-chromosomal major change allowing variation in IPI has occurred during the separation of the two D. virilis group phylads, the long IPI allowing variation also in PL (and CN). The evolution of the sounds in the virilis phylad has probably gone towards longer and denser pulse trains, while in the montana phylad the sounds have evolved in different directions.     

One- and two-dimensional polyacrylamide gel analysis of the heat shock proteins of the virilis group of Drosophila

The heat shock proteins of the virilis group of Drosophila are analyzed by one- and two-dimensional polyacrylamide gel analysis. This group consists of the two closely related but distinct virilis and montana phylads. The analysis reveals that some of the heat shock proteins are highly conserved among the two phylads while others are not. The 83-, 72-, and 69-kdalton proteins comigrate in all species examined. There is, however, a noticable trend toward greater molecular weight variability in the smaller heat shock proteins. In general, the heat shock protein patterns within each phylad follow the proposed phylogenetic relationships with some exceptions. D. ezoana and D. littoralis, both members of the montana phylad, exhibit heat shock protein patterns more similar to those of the virilis phylad. The data also demonstrate that the montana phylad has almost two times the heat shock allele members that the virilis phylad has. It is also shown that F₁ and F₂ hybrid flies of crosses between Drosophila species having different patterns of heat shock proteins show Mendelian segregation of alleles. After several generations of inbred growth, however, the pattern of heat shock protein synthesis in reciprocal hybrids each resembles that of the paternal parent. The implications of these findings are discussed.This research was supported in part by Damon Runyon-Walter Winchell Grant DRG-233F to R.M.S. and NIH Grant GM 27611 to R.V.S. R.V.S. is the recipient of an NIH Research Career Development Award.     

Morphological analysis of male mating organ in the Drosophila virilis species group: a multivariate approach

The shape of the male mating organ differs among 11 closely related species of the Drosophila virilis species group. Multivariate analyses of variation of a suite of 35 morphological traits (indices) describing the phallus shape were carried out in order to characterize interspecies variability of the traits. An overwhelming majority of the traits displayed species‐specific variability. The main result of the investigation was the revelation of the differences involved in the traits studied in the evolution of the D. virilis species group. The structure of species‐specific variability of some traits was discovered to be in accordance with the generally accepted taxonomy of the species group, while that of other traits required isolation of D. virilisper se from lummei phylad (former virilis phylad), and confirmed separation of montana phylad into three subphylads: montana proper, littoralis and kanekoi. Several subsets of traits having separate variability were determined in different parts of the male mating organ which correspond to spots of evolutionarily significant variability.     

Variation of wing shape in the Drosophila virilis species group (Diptera: Drosophilidae)

Wing shape variation was analysed with geometric morphometric methods in 17 laboratory strains, representing 11 closely related species (including two subspecies) of the Drosophila virilis group: D. virilis, D. lummei, D. novamexicana, D. americana americana, D. americana texana, D. montana, D. lacicola, D. flavomontana, D. borealis, D. littoralis, D. ezoana and D. kanekoi. Overall shape estimated using Procrustes coordinates of 14 landmarks was highly variable among strains and very similar in females and males. The landmarks in the distal part of the wing showed higher variation across strains than those in the proximal part. Procrustes distances between species were not consistent with phylogenetic distances previously suggested for the virilis group. Moreover, Procrustes distances between strains within species and within two major phylads (virilis and montana) were comparable with those between species and between phylads, respectively. The most different from other members of the group was the endemic D. kanekoi species, currently viewed as separate subphylad within the montana phylad. Allometric effects were found to be partly responsible for shape differences between the strains. Three most significant shape transformations were considered using the relative warp analysis and the strains were ordinated in accordance with transformation values. The pattern of relative warp scores could be easily interpreted only for the third warp explaining about 13% of shape variation. It separated the largest species, D. montana, D. ezoana and D. kanekoi, from other ones and was mainly associated with shape changes in the proximal region of the wing. The results of the present work suggest that wing shape in the virilis species group is not related to the speciation process. The observed proximal‐distal contrasts and allometric effects are in agreement with data of other studies, in which wing shape variation was analysed within Drosophila species.     

A RAPD fingerprinting of sibling species of the <Emphasis Type="Italic">Drosophila virilis</Emphasis> group

A comparative analysis of the sibling species of Drosophila virilis was performed by RAPD-PCR technique using a set of random primers. The degree of relatedness was studied by cluster analysis (UPGMA) and multidimensional scaling. The resulting pattern of species relationships contradicts the classical taxonomy. The main result of the cluster analysis is that D. virilis does not cluster with the remaining three species of its phylad, while according to multidimensional scaling, D. virilis is equidistant from all the species of its group, from both the species of its phylad and the species of the montana phylad. The montana phylad is extremely heterogeneous; moreover, the species D. littoralis, D. ezoana, and D. kanekoi appear to be closer to the virilis phylad than to the other species of the montana phylad, wherein these species are traditionally included. The phylogenetic relationships between the studied species discovered using RAPD fingerprinting comply with the results obtained using protein markers and quantitative traits.     

The DAIBAM MITE element is involved in the origin of one fixed and two polymorphic Drosophila virilis phylad inversions

Chromosomal inversions can originate from breakage and repair by non-homologous end-joining. Nevertheless, they can also originate from ectopic recombination between transposable elements located on the same chromosome inserted in opposite orientations. Here, we show that a MITE element (DAIBAM), previously involved in the origin of one Drosophila americana polymorphic inversion, is also involved in the origin of one fixed inversion between D. virilis and D. americana and another D. americana polymorphic inversion. Therefore, DAIBAM is responsible for at least 20% of the chromosomal rearrangements that are observed within and between species of the virilis phylad (D. virilis, D. lummei, D. novamexicana and D. americana), having thus played a significant role in the chromosomal evolution of this group of closely related species.     

Sequence evolution of theGpdh gene in theDrosophila virilis species group

The nucleotide sequence of theGpdh gene from six taxa,D. virilis, D. lummei, D. novamexicana, D. a. americana, D. a. texana andD. ezoana, belonging to thevirilis species group was determined to examine details of evolutionary change in the structure of theGpdh gene. TheGpdh gene is comprised of one 5 non-translated region, eight exons, seven introns and three 3 non-translated regions. Exon/intron organization was identical in all the species examined, but different from that of mammals. Interspecific nucleotide divergence in the entireGpdh gene followed the common pattern: it was low in the exon, high in the intron and intermediate in the non-translated regions. The degree of nucleotide divergence differed within these regions, suggesting that selection exerts constraints differentially on nucleotide change of theGpdh gene. A phylogenetic tree of thevirilis phylad constructed from nucleotide variation of total sequence was consistent with those obtained from other data.Nucleotide sequences for theGpdh gene ofD. lummei, D. novamexicana, D. a. americana, D. a. texana andD. ezoana have been submitted to GenBank with accession numbers D50087, D50088, D50089, D50090 and D50091.     

A multilocus microsatellite phylogeny of the Drosophila virilis group

We used a set of 48 polymorphic microsatellites derived from Drosophila virilis to infer phylogenetic relationships in the D. virilis clade. Consistent with previous studies, D. virilis and D. lummei were the most basal species of the group. Within the D. montana phylad, the phylogenetic relationship could not be resolved. Special attention was given to the differentiation between D. americana texana, D. americana americana and D. novamexicana. Significant differences between these three groups were detected by F(ST) analyses. Similarly, a model-based clustering method for multilocus genotype data also provided strong support for the presence of three differentiated groups. This genome-wide differentiation between D. americana texana and D. americana americana contrasts with previous analyses based on DNA sequence data.     

A RAPD fingerprinting of sibling species of the Drosophila virilis group

A comparative analysis of the sibling species of Drosophila virilis was performed by RAPD-PCR technique using a set of random primers. The degree of relatedness was studied by cluster analysis (UPGMA) and multi-dimensional scaling. The resulting pattern of species relationships contradicts the classical taxonomy. The main result of the cluster analysis is that D. virilis does not cluster with the remaining three species of its phylad, while according to multidimensional scaling, D. virilis is equidistant from all the species of its group, from both the species of its phylad and the species of the montana phylad. The montana phylad is extremely heterogeneous; moreover, the species D. littoralis, D. ezoana, and D. kanekoi appear to be closer to the virilis phylad than to the other species of the montana phylad, wherein these species are traditionally included. The phylogenetic relationships between the studied species discovered using RAPD fingerprinting comply with the results obtained using protein markers and quantitative traits.     

Content of histone H1 and histone phosphorylation in relation to the higher order structures of chromatin in Drosophila

The amount of histone H1 relative to core histones has been determined in three Drosophila species (D. melanogaster, D. texana and D. virilis) in chromatin from several tissues differing in chromatin structure and genetic activity. Low levels of H1 were found in relatively undifferentiated, early embryos as well as in a line of cultured cells. In late embryos the content of H1 was highest in D. virilis which possesses larger amounts of and a partially more compacted constitutive heterochromatin than the two other species. Polytene chromatin from larval salivary glands showed increased levels of H1 compared with diploid chromatin and the degree of phosphorylation of this histone was relatively low. The degree of phosphorylation of H2A was found to be drastically reduced in polytene as compared with diploid embryonic chromatin, which parallels the extensive underreplication of constitutive heterochromatin. Also, in diploid chromatin a qualitative correlation was observed between the relative amounts of heterochromatin and the levels of H2A phosphorylation. These findings suggest a connection between H2A phosphorylation and heavy compaction of interphase chromatin.     

Sites of the 5s Ribosomal Genes in Drosophila. I. the Multiple Clusters in the VIRILIS Group

Drosophila melanogaster ¹²⁵I-5S RNA was annealed to salivary gland preparations of 6 species in the virilis group of Drosophila. Two patterns of annealing were found. D. virilis, D. montana and D. borealis showed three 5S gene clusters on chromosome 5; S_d–f and W_c–j were strongly labeled, but X_a–e was weakly labeled. D. montana and D. borealis have a greater percentage of their total 5S cistrons at S _d–f than does D. virilis. D. americana americana, D. americana texana and D. novamexicana showed 2 sites labeled; no label was seen at S_d–f while W_c–j was weakly labeled and X_a–e was strongly labeled, the reciprocal of the previous pattern in the W-X region. Hybrids between D. a. americana and D. virilis showed no difference in chromosome banding at the sites of the 5S clusters despite their pattern differences. D. a texana x D. virilis, on the other hand, did show a difference in staining the X_a–e region. These patterns fall squarely into the biosystematic groupings deduced by many previous workers.     

Phylogeographic patterns in Drosophila montana

The Drosophila virilis species group offers valuable opportunities for studying the roles of chromosomal re-arrangements and mating signals in speciation. The 13 species are divided into two subgroups, the montana and virilis 'phylads'. There is greater differentiation among species within the montana phylad in both karyotype and acoustic signals than exists among members of the virilis phylad. Drosophila montana is a divergent species which is included in the montana phylad. Here, we analyse the phylogeography of D. montana to provide a framework for understanding divergence of acoustic signals among populations. We analysed mitochondrial sequences corresponding to the cytochrome oxidase I and cytochrome oxidase II genes, as well as 16 microsatellite loci, from 108 lines of D. montana covering most of the species' range. The species shows a clear genetic differentiation between North American and Scandinavian populations. Microsatellite allele frequencies and mitochondrial DNA haplotypes gave significant FST values between populations from Canada, USA and Finland. A Bayesian analysis of population structure based on the microsatellite frequencies showed four genetically distinct groups, corresponding to these three populations plus a small sample from Japan. A network based on mitochondrial haplotypes showed two Finnish clades of very different shape and variability, and another clade with all sequences from North America and Japan. All D. montana populations showed evidence of demographic expansion but the patterns inferred by coalescent analysis differed between populations. The divergence times between Scandinavian and North American clades were estimated to range from 450,000 to 900,000 years with populations in Canada and the USA possibly representing descendants of different refugial populations. Long-term separation of D. montana populations could have provided the opportunity for differentiation observed in male signal traits, especially carrier frequency of the song, but relaxation of sexual selection during population expansion may have been necessary.     

The effect of distamycin A on heterochromatin condensation of Drosophila chromosomes

Larval neuroblasts of four species of Drosophila (melanogaster, hydei, virilis and funebris) were treated with distamycin A, a DNA ligand which induces distinct undercondensation in AT-rich heterochromatin. For each species the patterns of undercondensation were correlated with distribution of quinacrine-bright regions and of satellite DNAs. An overlapping of distamycin A-sensitive and quinacrine-bright heterochromatic regions was demonstrated for D. melanogaster and D. virilis, but not for D. hydei and D. funebris. Distamycin A undercondensation is thus a further criterion for resolving heterochromatin into different parts and enables identification of the steps of the condensation process within heterochromatic regions.     

Nontranscribed spacers in Drosophila ribosomal DNA

Ribosomal DNA nontranscribed spacers in Drosophila virilis DNA have been examined in some detail by restriction site analysis of cloned segments of rDNA, nucleic acid hybridizations involving unfractionated rDNA, and base composition estimates. The overall G+C content of the spacer is 27–28%; this compares with 39% for rDNA as a whole, 40% for main band DNA, and 26% for the D. virilis satellites. Much of the spacer is comprised of 0.25 kb repeats revealed by digestion with Msp I, Fnu DII or Rsd I, which terminate very near the beginning of the template for the ribosomal RNA precursor. The spacers are heterogeneous in length among rDNA repeats, and this is largely accounted for by variation among rDNA units in the number of 0.25 kb elements per spacer. Despite its high A+T content and the repetitive nature of much of the spacer, and the proximity of rDNA and heterochromatin in Drosophila, pyrimidine tract analysis gave no indication of relatedness between the spacer and satellite DNA sequences. Species of Drosophila closely related to D. virilis have rDNA spacers that are homologous with those in D. virilis to the extent that hybridization of a cloned spacer segment of D. virilis rDNA to various DNA is comparable with hybridization to homologous DNA, and distributions of restriction enzyme cleavage sites are very similar (but not identical) among spacers of the various species. There is spacer length heterogeneity in the rDNA of all species, and each species has a unique major rDNA spacer length. Judging from Southern blot hybridization, D. hydei rDNA spacers have 20–30% sequence homology with D. virilis rDNA spacers, and a repetitive component is similarly sensitive to Msp I and Fnu DII digestion, D. melanogaster rDNA spacers have little or no homology with counterparts in D. virilis rDNA, despite a similar content of 0.25 kb repetitive elements. In contrast, sequences in rDNA that encode 18S and 28S ribosomal RNA have been highly conserved during the divergence of Drosophila species; this is inferred from interspecific hybridizations involving ribosomal RNA and a comparison of distributions of restriction enzyme cleavage sites in rDNA.Dedicated to Professor Wolfgang Beermann on the occasion of his sixtieth birthday     

Detection and location of three simple sequence DNAs in polytene chromosomes from virilis group species of Drosophila

In vitro synthesized RNAs complementary to the three satellite DNAs of Drosophila virilis have been used in a series of in situ hybridization experiments with polytene chromosomes from virilis group species. Gall and Atherton (1974) demonstrated that each of the satellites of D. virilis is comprised of many repeats of a distinct, seven base pair long, simple sequence. With few exceptions, copies of each of these simple sequences are detected in the chromocenters of all virilis group species. This is true even in species which do not possess satellite DNAs at buoyant densities corresponding to those of the satellite DNAs of D. virilis. Small quantities of the three simple sequences are also detected in euchromatic arms of several different species. The same euchromatic location may contain detectable copies of one, two, or all three simple sequence DNAs. The amounts of simple sequences at each location in the euchromatin may vary between species, between different stocks of the same species, and even between individuals of the same stock. The simple sequences located in the euchromatin appear to undergo DNA replication during formation of polytene chromosomes unlike those in heterochromatin. The locations of the euchromatic sequences are not the results of single chromosomal inversion events involving heterochromatic and euchromatic breakpoints.     

Histone content in relation to amount of heterochromatin and developmental stage in three species of Drosophila

Relative amounts of various histone fractions in Drosophila chromatin were estimated densitometrically on electrophoretic gel separations. Several consistent and highly significant differences were obtained between larval and adult chromatin. The arginine-rich histones showed the most conspicuous changes: higher amounts of H4 in larvae, higher H3 in adults. The level of modification of these histones was clearly higher in larval than in adult chromatin. The modification of the two slower subfractions of H4 involved, in all probability, phosphorylation as well as acetylation. In all types of Drosophila chromatin studied 50% or more of the H2a molecules were phosphorylated—a remarkably high proportion. The species differences observed in relative amounts of histone were consistent in both stages of development. D. melanogaster differed from D. hydei and D. virilis in all histones except H2b, while the latter two species were generally similar. The interspecific variation in histone pattern was generally not correlated to differences in content of heterochromatin. The level of modification of H3 was, however, presumably an exception, as it was significantly lower for both larvae and adults in D. virilis than in the other two species. These differ from D. virilis in containing appreciably lower proportions of heterochromatic chromosome segments.     

The DAIBAM MITE element is involved in the origin of one fixed and two polymorphic Drosophila virilis phylad inversions

Chromosomal inversions can originate from breakage and repair by non-homologous end-joining. Nevertheless, they can also originate from ectopic recombination between transposable elements located on the same chromosome inserted in opposite orientations. Here, we show that a MITE element (DAIBAM), previously involved in the origin of one Drosophila americana polymorphic inversion, is also involved in the origin of one fixed inversion between D. virilis and D. americana and another D. americana polymorphic inversion. Therefore, DAIBAM is responsible for at least 20% of the chromosomal rearrangements that are observed within and between species of the virilis phylad (D. virilis, D. lummei, D. novamexicana and D. americana), having thus played a significant role in the chromosomal evolution of this group of closely related species.     

Localization of <Emphasis Type="Italic">Dras1</Emphasis> gene in polytene chromosomes of sibling species of <Emphasis Type="Italic">Drosophila virilis</Emphasis> group

The Dras1 gene was mapped by in situ hybridization to polytene chromosomes of several sibling species of the Drosophila virilis group and their hybrids. A 1037-bp fragment of Dras1 gene from the D. virilis genome was used as the probe. The gene sequence was localized in the region of a 25 A-B disk in chromosome 2 (in accordance with the D. virilis polytene chromosome map (Gubenko and Evgen’ev, 1984).     

Isolation and characterization of microsatellites in Drosophila virilis and their cross species amplification in members of the D. virilis group

We report the isolation and cross species amplification of 42 Drosophila virilis microsatellite loci. Nine loci were isolated from mapped P1 bacteriophage clones and 33 were obtained from genomic DNA or GenBank searches. Cross species amplification was tested for all members of the D. virilis group. The amplification success was high (varying from 45% to 100%) and most of the loci were polymorphic. This set of loci can be applied for several genetic studies such as mapping behavioural quantitative trait loci (QTL) and for studying population structure in a phylogeographical framework in D. virilis group species.