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.     

Nucleotide variation at the no-on-transient A gene in Drosophila littoralis

The no-on-transient A (nonA) gene encodes a putative RNA-binding protein, and mutations in this gene are known to affect vision, male courtship song and viability in Drosophila melanogaster. Here we have sequenced the coding region of the nonA gene of Drosophila littoralis and compared it with those of Drosophila virilis and D. melanogaster. All portions of nonA appeared to be conserved between D. littoralis and D. virilis, while the 5' region of the gene of these two species showed high divergence from that of a more distantly-related species, D. melanogaster. The same was true for the glycine repeat regions. No significant deviation from neutrality was observed in the analysis of intraspecific nucleotide variation in 5' or 3' region of the nonA gene in D. littoralis population. Also, comparison of D. littoralis sequences with homologous sequence of D. virilis suggests that the gene is evolving neutrally in D. virilis group. Divergence of the 5' regions between D. virilis group species and D. melanogaster could be a result of positive selection, but this finding is obscured by the long divergence time of the species groups.     

Inheritance of species differences in female receptivity and song requirement between Drosophila virilis and D. montana

Females of two Drosophila virilis group species, D. virilis and D. montana, have different requirements for the courting males. In the present study we have examined species differences in female receptivity and male courtship song requirement using females' acceptance signal instead of copulation for measuring female readiness to mate. Behavior of D. virilis and D. montana females and F1 and backcross hybrid females was observed in a single-pair courtships with D. virilis and D. montana males and normal and wingless (mute) F1 hybrid males. D. virilis females were very receptive and they commonly accepted the courtship of males unable to produce courtship song. D. montana females, on the contrary, had a low receptivity and these females accepted the courting male only after hearing his song. Interspecific F1 and backcross (BCm) females resembled D. virilis more than D. montana in their receptivity. These females, however, resembled D. montana in their song requirement. These findings suggest that female song requirement and female receptivity are determined by different genetic factors.     

Chromosomal rearrangement In(2)TY and linkage maps of the second chromosome of Drosophila virilis

A chromosomal rearrangement In(2)TY, found in a natural population of Drosophila virilis, was examined cytologically, and its genetical characterization was carried out by the use of some marker isozymes and visible marker loci. From the interspecific hybridization between D. virilis and D. americana or D. texana, we deduced that the origin of In(2)TY was from a member of a closely related species of virilis.     

Unusual arrangement of the hsp68 locus in the virilis species group of Drosophila implicates evolutionary loss of an hsp68 gene.

Unlike all other Drosophila species studied to date, species in the virilis group of Drosophila have 2 complete copies of hsp68 arranged in inverted head-to-head orientation. Evidence for this conclusion includes Southern blots for D. virilis, D. lummei, and D. montana, PCR analysis of the former 2 species, in situ hybridization in D. virilis x D. lummei hybrids, and the complete nucleotide sequence of the locus in D. lummei. This organization resembles the primitive state of hsp70 in Diptera. Moreover, the Hsp68 peptide sequence for D. virilis and D. lummei is intermediate between that of Hsp70 and Hsp68 from other Drosophila spp. Therefore, we suggest that the hsp68 locus may have arisen via duplication of the hsp70 locus (or vice versa) early in the history of the genus Drosophila, with 1 hsp68 copy subsequently lost in most other Drosophila species groups.     

Localization of the Dras1 sequence to polytene chromosomes in sibling species of the Drosophila virilis group

Drasl gene was mapped by in situ hybridization to polytene chromosomes of several sibling species of the Drosophila virilis group and hybrids between them. A 1037 bp fragment of the Drasl gene of the D. virilis genome was used as a probe. The gene sequence is localized to the region of the disk 25 A-B on the chromosome 2 of the polytene chromosome map of D. virilis.     

Molecular phylogeny of the Drosophila virilis section (Diptera: Drosophilidae) based on mitochondrial and nuclear sequences

Regardless of the well-documented virilis species group, most groups of the Drosophila virilis section have not been completely studied at molecular level since it was suggested. Therefore, phylogenetic relationships among and within species groups of the virilis section are generally unknown. In present paper, the complete mitochondrial ND2 gene and fragment of COI gene in combination with a nuclear gene, Adh coding region, were used to derive the most extensive molecular phylogeny to date for the Drosophila virilis section. A total of 111 individuals covering 61 species were sampled in this study. Novel phylogenetic findings included (1) support for the paraphyly of the melanica and robusta species group and at least two subgroups of the robusta species group, the lacertosa and okadai subgroups, were distinguished as paraphyletic taxa. In addition, (2) present results revealed the sister relationship between D. moriwakii and the robusta subgroup, conflicting with current taxonomy regarding D. moriwakii, which was shifted from the robusta species group to the melanica group. (3) In contrast to the robusta and melanica species groups, monophyly of the polychaeta species group, the angor group and the virilis group was confirmed, respectively. However, the monophyletic quadrisetata species group was resolved with uncertainty. (4) Our analyses of combined data set suggested close relationship between the quadrisetata species group and the unpublished clefta group, and the okadai subgroup is sister to the clade comprising of the quadrisetata and clefta species groups. Within the virilis section, D. fluvialis and three tropical species groups, the polychaeta group, the angor group and the repleta group, are found to branch off earlier than other ingroup taxa. This suggests that the virilis section might have originated in the Old World tropics. Besides, the derived status of the close affinities of the quadrisetata group, the clefta group, and the melanica and robusta groups is probably the result of their adaptation to forests between subtropical and cool-temperate climate. Based on the consideration of the phylogenetic placement of the species of the virilis section, we suggest that at least five independent migrations occurred from the Old World to the New World.     

A 9.6 kb intervening sequence in D. virilis rDNA, and sequence homology in rDNA interruptions of diverse species of Drosophila and other diptera.

A large proportion of the 28S ribosomal RNA genes in Drosophila virilis are interrupted by a DNA sequence 9.6 kilobase pairs long. As regards both its presence and its position in the 28S gene (about two thirds of the way in), the D. virilis rDNA intervening sequence is similar to that found in D. melanogaster rDNA, but lengths differ markedly between the two species. Degrees of nucleotide sequence homology have been detected bewteen rDNA interruptions of the two species. This homology extends to putative rDNA intervening sequences in diverse higher diptera (other Drosophila species, the house fly and the flesh fly), but hybridization of cloned D. melanogaster and D. virilis rDNA interruption segments to DNA of several lower diptera has been negative. As is the case with melanogaster rDNA interruptions, segments of the virilis rDNA intervening sequence hybridize with non-rDNA components of the virilis genome, and interspecific homology may involve these non-rDNA sequences as well as rDNA interruptions. There is, however, evidence from buoyant density fractionation of DNA that the distributions of interruption-related sequences are distinct in D. melanogaster and D. virilis genomes. Moreover, thermal denaturation studies have indicated differing extents of homology between hybridizable sequences in D. virilis DNA and different segments of the D. melanogaster rDNA intervening sequence. We infer from our studies that rDNA intervening sequences are prevalent among higher diptera; that in the course of the evolution of these organisms, elements of the intervening sequences have been moderately to highly conserved; and that this conservation extends in at least two distantly related species of Drosophila to similar sequences found elsewhere in the genomes.     

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 P1 clones containing ~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.     

Drosophila virilis has long and highly polymorphic microsatellites

Comparative genomics is a powerful approach to inference of the dynamics of genome evolution. Most information about the evolution of microsatellites in the genus Drosophila has been obtained from Drosophila melanogaster. For comparison, we collected microsatellite data for the distantly related species Drosophila virilis. Screening about 0.5 Mb of nonredundant genomic sequence from GenBank, we identified 239 dinucleotide microsatellites. On average, D. virilis dinucleotides were significantly longer than D. melanogaster microsatellites (7.69 repeats vs. 6.75 repeats). Similarly, direct cloning of microsatellites resulted in a higher mean repeat number in D. virilis than in D. melanogaster (12.7 repeats vs. 12.2 repeats). Characterization of 11 microsatellite loci mapping to division 40-49 on the fourth chromosome of D. virilis indicated that D. virilis microsatellites are more variable than those of D. melanogaster.     

A genealogical view of chromosomal evolution and species delimitation in the Drosophila virilis species subgroup

Chromosomal arrangement was a historically important character used for defining taxonomic boundaries. The Drosophila virilis species group exhibits a series of chromosomal rearrangements, and the resulting differences among karyotypes were primary characters originally used to define taxa within the group. However, some chromosomally divergent forms have not been sufficiently resolved in phylogenetic reconstructions of DNA sequences from several nuclear genes. Sequences of mitochondrial regions have the potential for finer-scale resolution of closely related taxa; therefore, sequences of two mitochondrial genes were used to examine phylogenetic relationships within the chromosomally variable virilis subgroup. Sequences were obtained from multiple strains of the Palearctic species, D. virilis and D. lummei, and the Nearctic species, D. novamexicana and two chromosomal forms of D. americana. Analyses support the recent emergence of the different chromosomal forms in North America. However, none of these chromosomally divergent forms exhibit reciprocal monophyly of their mtDNA sequences, which is the requirement for attaining genealogical species status.     

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.     

Invasion of Drosophila virilis by the Penelope transposable element

The Penelope family of transposable elements (TEs) is broadly distributed in most species of the virilis species group of Drosophila. This element plays a pivotal role in hybrid dysgenesis in Drosophila virilis, in which at least four additional TE families are also activated. Here we present evidence that the Penelope family of elements has recently invaded D. virilis. This evidence includes: (1) a patchy geographical distribution, (2) genomic locations mainly restricted to euchromatic chromosome arms in various geographical strains, and (3) a high level of nucleotide similarity among members of the family. Two samples from a Tashkent (Middle Asia) population of D. virilis provide further support for the invasion hypothesis. The 1968 Tashkent strain is free of Penelope sequences, but all individuals collected from a 1997 population carry at least five Penelope copies. Furthermore, a second TE, Ulysses, has amplified and spread in this population. These results provide evidence for the Penelope invasion of a D. virilis natural population and the mobilization of unrelated resident transposons following the invasion.     

DNA sequence divergence in the Drosophila virilis group

DNA sequence divergence was analyzed in some sibling species of the Drosophila virilis group. Clones comprising about 0.1% of the genome DNA were selected at random from a D. virilis library for a comparative study on DNA from D. lummei, D. novamexicana, D. borealis, and D. lacicola. Blot hybridization experiments indicated that about 70% of DNA from D. lummei and D. novamexicana and less than 50% of DNA from D. borealis and D. lacicola share sequences that are homologous to DNA in D. virilis. This finding is in excellent agreement with the genealogical tree based on cytological studies (Throckmorton 1982). - Four plasmids with inserts which are present in one or a few copies per genome were hybridized in situ to polytene chromosomes. These experiments demonstrate that (1) homologous "unique" DNA sequences are localized exclusively in homologous bands and (2) homologous bands that appear to be identical in different species may contain different DNA sequences.     

The relationships among the species of the Drosophila virilis group inferred from the gene Ras1 sequences

Comparative analysis of a group of closely related Drosophila species (D. virilis, D. lummei, D. novamexicana, D. americana texana, D. flavomontana, D. montana, D. borealis, D. lacicola, D. littoralis, D. kanekoi, and D. ezoana) was conducted based on an incomplete sequence of gene Ras1. The pattern of the relationships among the species corresponded to that expected from analysis of morphological and cytogenetic characters. Statistical data favoring neutrality of the substitutions examined in the Ras1 gene are presented. This character of the gene Ras1 evolution confers more reliability to reconstruction of phylogenetic relationships among closely related species. The resultant tree for main phylads of the group is as follows: (D. virilis, D. lummei, D. montana, D. ezoana).     

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.     

Autonomy of association frequencies of telomeric regions in salivary gland chromosomes during interspecific hybridization]

The study of association frequency and combination of salivary gland chromosomes participating in telomeric associations in interspecific hybrids between related species of "virilis" group of Drosophila (D. virilis, D. texana and D. imeretensis Sokolov) enables to conclude that each homologous chromosome behaves as an individual unit in the process. Besides, it is shown that the frequency of association of an autonomous property of telomere, i. e. it is not changed in interspecific hybrids in comparison with parental species.