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.     

Distribution and evolution of mobile elements in the virilis species group of Drosophila

The distributions of Penelope and Ulysses, two transposable elements that can induce hybrid dysgenesis, were studied in several species groups of Drosophila. No significant hybridization to Penelope and Ulysses probes was detected by Southern blot analyses of species outside the virilis group. In contrast, both element families have had a long residence in all species of the virilis species group, as indicated by their strong presence in the heterochromatic chromocenter. Except for D. kanekoi, D. lummei, and some strains of D. virilis, species of the group carry full-sized, and at least potentially functional, copies of both element families. Consistent with the occurrence of recent transposition, Penelope and Ulysses elements are located at different chromosomal sites in different geographical strains of the same species. A total of 79 Penelope and 47 Ulysses euchromatic insertion sites were localized to chromosomal subsections in species of the virilis group. Highly significant deviations from independence of the distributions of Penelope and Ulysses and previously established inversion breakpoints were documented, suggesting that these transposable elements may have played an important role in genomic reorganization and evolution of the virilis species group, which is especially rich in karyotypic variation. Received: 13 April 1999; in revised form: 20 July 1999 / Accepted: 27 July 1999     

Microsatellite-based species identification method for Drosophila virilis group species

Species of the D. virilis group are widely used in evolutionary research, but the individuals of different species are difficult to distinguish from each other morphologically. We constructed a fast and easy microsatellite-based identification method for the species of the group occurring sympatrically in northern Europe. The neighbor joining tree based on 14 microsatellite loci also gave a good resolution of the species divergence pattern in the whole group.     

Aggregation pheromones in five taxa of the Drosophila virilis species group

ABSTRACT. Aggregation pheromones have been demonstrated in the closely related taxa: Drosophila americana americana Spencer, D. a. texana Patterson, D. novamexicana Patterson, and D. lummei Hackman. These pheromones function much as has been reported previously for D. virilis Sturtevant. The compounds are produced by sexually mature males, but both sexes respond in a wind-tunnel olfactometer. In all species except D. lummei , a 21-carbon alkene is an important pheromone component. In D. virilis the hydrocarbon is (Z)-10-heneicosene (Z10–21), but in D. a. americana, D. a. texana and D. novamexicana it is (Z)-9-heneicosene (Z9-21). All these taxa respond best to the heneicosene which they produce. D. lummei possesses no heneicosenes but, curiously, responds well to both Z9-21 and Z10-21. All species possess five male-specific esters which were previously discovered in D. virilis : methyl tiglate, ethyl tiglate, isopropyl tiglate, methyl hexanoate and ethyl hexanoate. Ethyl tiglate is the most abundant in each case. Responses to the esters vary among the taxa, ranging from highly significant in D. lummei , particularly to ethyl tiglate, to not demonstrable in D. a. americana. Variability in ester response has also been demonstrated between two strains of D. virilis. In all cases the crude male-derived pheromone is synergistic with an extract of fermented willow bark, on which oviposition is said to occur.     

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.     

Resolving the phylogenetic relationships and evolutionary history of the Drosophila virilis group using multilocus data

The Drosophila virilis group is one of the major lineages of Drosophila previously recognised and it has been used as a model for different types of studies. It comprises 13 species whose phylogenetic relationships are not well resolved. In the present study, six nuclear genes (Adh, fused, Gpdh, NonA, CG9631 and CG7219) and the mitochondrial ribosomal RNA genes (12S-16S) have been used to estimate the evolutionary tree of the group using different methods of phylogenetic reconstruction. Different competing evolutionary hypotheses have also been compared using the Approximately Unbiased test to further evaluate the robustness of the inferred trees. Results are, in general, consistent with previous studies in recovering the four major lineages of the group (D. virilis phylad, Drosophila montana subphylad, Drosophila kanekoi subphylad and Drosophila littoralis subphylad), although D. kanekoi, D. littoralis and Drosophila ezoana are here inferred to be more closely related to the D. virilis phylad than to the D. montana subphylad. The age of the crown group, estimated with a Bayesian method that assumes a relaxed molecular clock, is placed in the late Miocene (～ 10 Mya). The oldest lineages also appeared during this period (～ 7.5 to ～ 8.9 Mya), while the ages of the basal nodes of the montana subphylad and the virilis phylad are located in the early Pliocene (～ 4.9 and ～ 4.1 Mya). Major cladogenesis events correlate to geological and palaeoclimatic occurrences that most likely affected the freshwater and deciduous forests where these species are found. The inferred biogeographical history of the group, based on the statistical dispersal-vicariance analysis, indicates that the last common ancestor of the group had a Holarctic distribution from which the North American and the Eurasian lineages evolved as a result of a vicariant event.     

Mitochondrial DNA polymorphism in the natural populations of the Drosophila virilis species group

Restriction fragment length polymorphism (RFLP) analysis has been used to evaluate mitochondrial DNA (mtDNA) variation in 12 sibling species forming the Drosophila virilis species group. The variation thresholds corresponding to the interspecific and interstrain levels have been determined. The results indicate that interspecific hybridization has significantly contributed to the evolutionary history of the virilis species group.     

The Evolutionary History of the Transposable Element Penelope in the Drosophila virilis Group of Species

We have used phylogenetic techniques to study the evolutionary history of the Penelope transposable element in the Drosophila virilis species group. Two divergent types of Penelope have been detected, one previously described, clade I, and a new one which we have termed clade III. The phylogeny of some copies of the Penelope clade I element was partially consistent with the species phylogeny of the D. montana subphylad, suggesting cospeciation and allowing the estimation of the evolutionary rate of Penelope. Divergence times of elements found in different species are younger than the age of the species, suggesting horizontal transfer events. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Dmitri Petrov]     

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.     

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.     

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

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.     

Association of telomeric regions of salivary gland chromosomes in Drosophila species of the virilis group]

The frequency of participating salivary gland chromosomes in telomeric associations and specific combinations of associated chromosomes were studied in the following species belonging to "virilis" group of Drosophila: D. virilis, D. texana, D. littoralis Sokolov and D. imeretensis Sokolov. In all species studied predominantly telomeres of two chromosomes form an association. In females of the species mentioned the X-chromosome takes part in association twice as frequent as in males. The total length of the chromosome and the presence of subterminal inversions may sometimes influence the frequency of participating the chromosome in associations. However, the data obtained enable to postulate that internal homology of telomeric regions plays the main role in the process. The equal frequency of participating in telomeric associations for definite chromosomes in species studied may resemble the phylogenetic relation of the species.     

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.