Molecular Phylogeny of the <Emphasis Type="Italic">Drosophila melanogaster</Emphasis> Species Subgroup

Although molecular and phenotypic evolution have been studied extensively in Drosophila melanogaster and its close relatives, phylogenetic relationships within the D. melanogaster species subgroup remain unresolved. In particular, recent molecular studies have not converged on the branching orders of the D. yakuba–D. teissieri and D. erecta–D. orena species pairs relative to the D. melanogaster–D. simulans–D. mauritiana–D. sechellia species complex. Here, we reconstruct the phylogeny of the melanogaster species subgroup using DNA sequence data from four nuclear genes. We have employed vectorette PCR to obtain sequence data for orthologous regions of the Alcohol dehydrogenase (Adh), Alcohol dehydrogenase related (Adhr), Glucose dehydrogenase (Gld), and rosy (ry) genes (totaling 7164 bp) from six melanogaster subgroup species (D. melanogaster, D. simulans, D. teissieri, D. yakuba, D. erecta, and D. orena) and three species from subgroups outside the melanogaster species subgroup [D. eugracilis (eugracilis subgroup), D. mimetica (suzukii subgroup), and D. lutescens (takahashii subgroup)]. Relationships within the D. simulans complex are not addressed. Phylogenetic analyses employing maximum parsimony, neighbor-joining, and maximum likelihood methods strongly support a D. yakuba–D. teissieri and D. erecta–D. orena clade within the melanogaster species subgroup. D. eugracilis is grouped closer to the melanogaster subgroup than a D. mimetica–D. lutescens clade. This tree topology is supported by reconstructions employing simple (single parameter) and more complex (nonreversible) substitution models. Present address (Ryan M. David): University of Texas Health Science Center at San Antonio, 7703 Floyd Curl Drive, San Antonio, TX 78284, USA     

Discrepancy in divergence of the mitochondrial and nuclear genomes ofDrosophila teissieri andDrosophila yakuba

Summary Restriction sites were compared in the mitochondrial DNA (mtDNA) molecules from representatives of two closely related species of fruit flies: nine strains ofDrosophila teissieri and eight strains ofDrosophila yakuba. Nucleotide diversities amongD. teissieri strains and amongD. yakuba strains were 0.07% and 0.03%, respectively, and the nucleotide distance between the species was 0.22%. Also determined was the nucleotide sequence of a 2305-nucleotide pari (ntp) segment of the mtDNA molecule ofD. teissieri that contains the noncoding adenine+thymine (A+T)-rich region (1091 ntp) as well as the genes for the mitochondrial small-subunit rRNA, tRNA^f-met, tRNA^gln, and tRNA^ile, and portions of the ND2 and tRNA^val genes. This sequence differs from the corresponding segment of theD. yakuba mtDNA by base substitutions at 0.1% and 0.8% of the positions in the coding and noncoding regions, respectively. The higher divergence due to base substitutions in the A+T-rich region is accompanied by a greater number of insertions/deletions than in the coding regions. From alignment of theD. teissieri A+T-rich sequence with those ofD. yakuba andDrosophila virilis, it appears that the 40% of this sequence that lies adjacent to the tRNA^ile gene has been highly conserved. Divergence between the entireD. teissieri andD. yakuba mtDNA molecules, estimated from the sequences, was 0.3%; this value is close to the value (0.22%) obtained from the restriction analysis, but 10 times lower than the value estimated from published DNA hybridization results. From consideration of the relationships of mitochondrial nucleotide distance and allozyme genetic distance found among seven species of theDrosophila melanogaster subgroup, the mitochondrial nucleotide distance observed forD. teissieri andD. yakuba is anomalously low in relation to the nuclear genetic distance.     

Adjacent chromosomal regions can evolve at very different rates: Evolution of theDrosophila 68C glue gene cluster

Summary The 68C puff is a highly transcribed region of theDrosophila melanogaster salivary gland polytene chromosomes. Three different classes of messenger RNA originate in a 5000-bp region in the puff; each class is translated to one of the salivary gland glue proteins sgs-3, sgs-7, or sgs-8. These messenger RNA classes are coordinately controlled, with each RNA appearing in the third larval instar and disappearing at the time of puparium formation. Their disappearance is initiated by the action of the steroid hormone ecdysterone. In the work reported here, we studied evolution of this hormone-regulated gene cluster in themelanogaster species subgroup ofDrosophila. Genome blot hybridization experiments showed that five other species of this subgroup have DNA sequences that hybridize toD. melanogaster 68C sequences, and that these sequences are divided into a highly conserved region, which does not contain the glue genes, and an extraordinarily diverged region, which does. Molecular cloning of this DNA fromD. simulans, D. erecta, D. yakuba, andD. teissieri confirmed the division of the region into a slowly and a rapidly evolving protion, and also showed that the rapidly evolving region of each species codes for third instar larval salivary gland RNAs homologous to theD. melanogaster glue mRNAs. The highly conserved region is at least 13,000 bp long, and is not known to code for any RNAs.     

Evidence for interspecific transfer of the transposable element mariner betweenDrosophila andZaprionus

Summary The transposable element mariner occurs widely in themelanogaster species group ofDrosophila. However, in drosophilids outside of themelanogaster species group, sequences showing strong DNA hybridization with mariner are found only in the genusZaprionus. the mariner sequence obtained fromZaprionus tuberculatus is 97% identical with that fromDrosophila mauritiana, a member of themelanogaster species subgroup, whereas a mariner sequence isolated fromDrosophila tsacasi is only 92% identical with that fromD. mauritiana. BecauseD. tsacasi is much more closely related toD. mauritiana than isZaprionus, the presence of mariner inZaprionus may result from horizontal transfer. In order to confirm lack of a close phylogenetic relationship between the genusZaprionus and themelanogaster species group, we compared the alcohol dehydrogenase (Adh) sequences among these species. The results show that the coding region of Adh is only 82% identical betweenZ. tuberculatus andD. mauritiana, as compared with 90% identical betweenD. tsacasi andD. mauritiana. Furthermore, the mariner gene phylogeny obtained by maximum likelihood and maximum parsimony analyses is discordant with the species phylogeny estimated by using the Adh genes. The only inconsistency in the mariner gene phylogeny is in the placement of theZaprionus mariner sequence, which clusters with mariner fromDrosophila teissieri andDrosophila yakuba in themelanogaster species subgroup. These results strongly suggest horizontal transfer.     

Evolution of the transposable element mariner in the Drosophila melanogaster species group

The population biology and molecular evolution of the transposable element mariner has been studied in the eight species of the melanogaster subgroup of the Drosophila subgenus Sophophora. The element occurs in D. simulans, D. mauritiana, D. sechellia, D. teissieri, and D. yakuba, but is not found in D. melanogaster, D. erecta, or D. orena. Sequence comparisons suggest that the mariner element was present in the ancestor of the species subgroup and was lost in some of the lineages. Most species contain both active and inactive mariner elements. A deletion of most of the 3 end characterizes many elements in D. teissieri, but in other species the inactive elements differ from active ones only by simple nucleotide substitutions or small additions/deletions. Active mariner elements from all species are quite similar in nucleotide sequence, although there are some-species-specific differences. Many, but not all, of the inactive elements are also quite closely related. The genome of D. mauritiana contains 20–30 copies of mariner, that of D. simulans 0–10, and that of D. sechellia only two copies (at fixed positions in the genome). The mariner situation in D. sechellia may reflect a reduced effective population size owing to the restricted geographical range of this species and its ecological specialization to the fruit of Morinda citrifolia.     

Drosophila mitochondrial DNA: Conserved sequences in the A+T-rich region and supporting evidence for a secondary structure model of the small ribosomal RNA

Summary The sequence of a segment of theDrosophila virilis mitochondrial DNA (mtDNA) molecule that contains the A+T-rich region, the small rRNA gene, the tRNA^f-met, tRNA^gln, and tRNA^ile genes, and portions of the ND2 and tRNA^val genes is presented and compared with the corresponding segment of theD. yakuba mtDNA molecule. The A+T-rich regions ofD. virilis andD. yakuba contain two correspondingly located sequences of 49 and 276/274 nucleotides that appear to have been conserved during evolution. In each species the replication origin of the mtDNA molecule is calculated to lie within a region that overlaps the larger conserved sequence, and within this overlap is found a potential hairpin structure. Substitutions between the larger conserved sequences of the A+T-rich regions, the small mt-rRNA genes, and the ND2 genes are biased in favor of transversions, 71–97% of which are AT changes. There is a 13.8 times higher frequency of nucleotide differences between the 5 halves than between the 3 halves of theD. virilis andD. yakuba small mt-rRNA genes. Considerations of the effects of observed substitutions and deletion/insertions on possible nucleotide pairing within the small mt-rRNA genes ofD. virilis andD. yakuba strongly support the secondary structure model for theDrosophila small mt-rRNA that we previously proposed.     

Structural evolution of the Drosophila 5S ribosomal genes

We compare the 5S gene structure from nine Drosophila species. New sequence data (5S genes of D. melanogaster, D. mauritiana, D. sechellia, D. yakuba, D. erecta, D. orena, and D. takahashii) and already-published data (5S genes of D. melanogaster, D. simulans, and D. teissieri) are used in these comparisons. We show that four regions within the Drosophila 5S genes display distinct rates of evolution: the coding region (120 bp), the 5-flanking region (54–55 bp), the 3-flanking region (21–22 bp), and the internal spacer (149–206 bp). Intra- and interspecific heterogeneity is due mainly to insertions and deletions of 6–17-bp oligomers. These small rearrangements could be generated by fork slippages during replication and could produce rapid sequence divergence in a limited number of steps. Correspondence to: M. Wegnez     

Evolution of the threonine-glycine repeat region of the period gene in the melanogaster species subgroup of drosophila

Summary The Threonine-Glycine (Thr-Gly) region of the period gene (per) in Drosophila was compared in the eight species of the D. melanogaster subgroup. This region can be divided into a diverged variable-length segment which is flanked by more conserved sequences. The number of amino acids encoded in the variable-length region ranges from 40 in D. teissieri to 69 in D. mauritiana. This is similar to the range found within natural populations of D. melanogaster. It was possible to derive a Thr-Gly allele of one species from that of another by invoking hypothetical Thr-Gly intermediates. A phylogeny based on the more conserved flanking sequences was produced. The results highlighted some of the problems which are encountered when highly polymorphic genes are used to infer phylogenies of closely related species.     

Fruiting studies in species ofCapronia (Herpotrichiellaceae)

A recently described protocol for thein vitro production of ascomata was employed to determine the sexual incompatibility systems of five species ofCapronia. The formation of mature ascomata in isolates derived from single ascospores demonstrated thatC. epimyces, C. mansonii, andC. munkii n. sp. are homothallic. In contrast, fertile ascomata were observed only in mass-ascospore isolates and pairwise crosses between specific single-ascospore isolates inC. dactylotricha n. sp. andC. moravica. TheExophiala anamorphs ofC. dactylotricha andC. munkii are described and aPhialophora-like synanamorph is reported for the former species. Germinating ascospores ofC. munkii formed conidiogenous cells directly, while the ascospores of the remaining species germinated to produce germ tubes and hyphae. The application of the terms microcyclic conidiation to secondary conidium production and sclerotial body and stroma to the multicellular structures produced by species ofCapronia andExophiala are discussed.     

Development time and pupation behavior in theDrosophila melanogaster subgroup (Diptera: Drosophilidae)

This study is an in-depth analysis of intersexual, intraspecific, and interspecific variability in larvopupal developmental time, pupation site preference, and larval and pupal survival of a number of isofemale lines of the speciesDrosophila mauritiana, D. melanogaster, D. sechellia, D. simulans, D. teissieri, andD. yakuba. There was no significant sex differences in pupation height, but females eclosed significantly earlier than males in all species. In addition, the suggestion of a strong negative correlation between larval developmental time and pupation height could not be confirmed in this study. The hypothesis that differences in pupation height provide a basis for niche partitioning between closely related species with overlapping distributions was tested by three planned orthogonal contrast analyses of variance. First, the two speciesD. teissieri andD. yakuba, with largely overlapping distribution, were significantly different in pupation height. Second, the two allopatric, nonoverlapping island speciesD. mauritiana andD. sechellia did not significantly differ in pupation height. However, the absence of a significant difference in the final contrast between the two cosmopolitan speciesD. melanogaster andD. simulans, which are often found together, makes us cautious to accept the hypothesis.     

An approach to study the evolution of theDrosophila 5S ribosomal genes using P-element transformation

Summary The P-element-mediated gene transfer system was used to introduceDrosophila teissieri 5S genes into theDrosophila melanogaster genome. Eight transformedD. melanogaster strains that carryD. teissieri 5S mini-clusters consisting of 9–21 adjacent 5S units were characterized. No genetic exchanges betweenD. melanogaster andD. teissieri 5S clusters were detected over a 2-year survey of the eight strains. The occurrence of small rearrangements within theD. melanogaster 5S cluster was demonstrated in one of the transformed strains.     

The LSP1-α gene is not dosage compensated in the Drosophila melanogaster species subgroup

The X-linked subunit of larval serum protein 1 (LSP1-) is shown to lack dosage compensation in six members of the melanogaster species subgroup, viz., Drosophila melanogaster, D. simulans, D. mauritiana, D. erecta, D. yakuba, and D. teissieri, by quantitative filter hybridization and by electrophoretic and autoradiographic analyses of fat body proteins. These results support the hypothesis that there is little selection pressure on the LSP1- gene to acquire dosage compensation.H.W.B. acknowledges a Commonwealth Scholarship. This work was supported by an SRC grant to D.B.R.     

The β-tubulin gene family evolution in theDrosophila montium subgroup of themelanogaster species group

The 1-, 2-, and 3-tubulin genes have been mapped by in situ hybridization on the polytene chromosomes of 11 selected species (15 strains) belonging to theDrosophila montium subgroup. Although the hybridization pattern among the strains of the same species does not differ, this pattern is significantly different among the species. The -tubulin genes in themontium subgroup seem to be organized in a cluster, or in a semi-cluster, or are completely dispersed. The clustered arrangement is found in the North-Oriental sibling speciesD. auraria, D. triauraria, andD. quadraria. The semi-clustered arrangement, wherein the 1 and 2 genes are located at the same locus while 3 is at a different one, appears in the South-Oriental speciesD. bicomuta, D. serrata, andD. birchii, as well as in the Afrotropical speciesD. diplacantha andD. seguyi. The complete separation of the genes is observed in the Indian speciesD. kikkawai andD. jambulina and in the Afrotropical speciesD. vulcana. Based on the above results, a possible mode of evolution of the -tubulin genes in the montium subgroup is attempted. In addition, phylogenetic relationships among themontium species are discussed. Correspondence to: Z.G. Scouras     

A comparative ultrastructural study of ‘glue’ production and secretion of the salivary glands in different species of theDrosophila melanogaster group

Summary Salivary gland cells of members of theDrosophila melanogaster group (from four different subgroups) were examined electron microscopically and histochemically during the late larval period of development. The secretory product, which is supposed to be utilized as glue at the time of puparium formation, appears, by analogy to Palade and Jamieson's results, to be synthesized partially in the rough endoplasmic reticulum (RER) and partially in the Golgi complex. The latter is also the usual site of the packaging of the product into secretory granules, except in the case of one of the secretory granule components ofD. lucipennis. The phylogenetic relationships among the subgroups, implied by the morphological appearance of the secretory granules, fit well with the existing phylogenetic relationships within the group. The secretory granules of each species have their own morphological features; granules of species of the same subgroup share some of these features. Secretion occurs from the cells via exocytosis during which the morphology of the secretory granules changes. Light microscope examination of PAS (Periodic Acid-Schiff reaction) stained glands shows a strong positive reaction in most species, with the exception of the species of thesuzukii subgroup which show a weak, or a negative reaction (D. rajasekari). Electron histochemical localization of polysaccharides in the secretory granules was possible inD. melanogaster and the species of theananassae subgroup.     

Globin evolution in the genusXenopus: Comparative analysis of cDNAs coding for adult globin polypeptides ofXenopus borealis andXenopus tropicalis

Summary Globin mRNAs ofXenopus borealis andXenopus tropicalis have been cloned and sequenced. The nucleotide and derived amino acid sequences were compared with each other and with already available data fromXenopus laevis. This analysis rendered clear evidence that the common ancestor ofX. laevis andX. borealis, but not ofX. tropicalis, had lost one amino acid of the -globins prior to a genome duplication event that preceded the segregation of the former two species. Replacement-site substitutions were used to calculate a rough time scale of genome duplication and species segregation. The results suggest an ancient separation between theX. laevis and theX. tropicalis groups occurring approximately 110–120 million years ago. Analysis of the amino acid chains demonstrated various alterations. However, some functional domains, like heme-binding sites and₁₂ contact sites, were subject to a high degree of conservation, indicating the existence of functional constraints on them also in the genusXenopus.     

The evolutionary genetics of thehobo transposable element in theDrosophila melanogaster complex

Hobo elements are a family of transposable elements found inDrosophila melanogaster and its three sibling species:D. simulans, D. mauritiana andD. sechellia. Studies inD. melanogaster have shown thathobo may be mobilized, and that the genetic effects of such mobilizations included the general features of hybrid dysgenesis: mutations, chromosomal rearrangements and gonadal dysgenis in F1 individuals. At the evolutionary level somehobo-hybridizing sequences have also been found in the other members of themelanogaster subgroup and in many members of the relatedmontium subgroup. Surveys of older collected strains ofD. melanogaster suggest that completehobo elements were absent prior to 50 years ago and that they have recently been introduced into this species by horizontal transfer. In this paper we review our findings and those of others, in order to precisely describe the geographical distribution and the evolutionary history ofhobo in theD. melanogaster complex. Studies of the DNA sequences reveal a different level of divergence between the groupD. melanogaster, D. simulans andD. mauritiana and the fourth speciesD. sechellia. The hypothesis of multiple transfers in the recent past into theD. melanogaster complex from a common outside source is discussed.     

Drosophila virilis histone gene clusters lacking H1 coding segments

Summary Approximately 30–40% ofDrosophila virilis DNA complementary to clonedDrosophila histone genes is reduced to 3.4-kilobase-pair (kbp) segments by Bgl I or Bgl II digestion. The core histone genes of a 3.4-kbp Bgl II segment cloned in the plasmid pDv3/3.4 have the same order as theD. melanogaster core histone genes in the plasmid cDm500: . Nonetheless, pDv3/3.4 and cDm500 have different histone gene configurations: In pDv3/3.4, the region between the H2B and H3 genes contains 0.35 kbp and cannot encode histone H1; in cDm500, the region contains 2.0 kbp and encodes histone H1. The lack of an H1 gene between the H2B and H3 genes in 30–40% ofD. virilis histone gene clusters suggests that changes in histone gene arrays have occurred during the evolution ofDrosophila. The ancestors of modernDrosophila may have possessed multiple varieties of histone gene clusters, which were subsequently lost differentially in thevirilis andmelanogaster lineages. Alternatively, they may have possessed a single variety, which was rearranged during evolution. The H1 genes ofD. virilis andD. melanogaster did not cross-hybridize in vitro under conditions that maintain stable duplexes between DNAs that are 75% homologous. Consequently,D. virilis H1 genes could not be visualized by hybridization to an H1-specific probe and thus remain unidentified. Our observations suggest that the coding segments in the H1 genes ofD. virilis andD. melanogaster are >25% divergent. Our estimate of sequence divergence in the H1 genes ofD. virilis andD. melanogaster seems high until one considers that the coding sequences of cloned H1 genes from the closely related speciesD. melanogaster andD. simulans are 5% divergent.     

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.     

DNA base composition within the genusScenedesmus (Chlorophyta)

The DNA base compositions of 60 strains ofScenedesmus were determined and found to be a valuable indicator for the differentiation within this genus (except for the very closely related species of the subsect.Desmodesmus). The range for allScenedesmus species was 49.9–69.3 mol% GC for nuclear DNA and 36.8–39.9 mol% GC for chloroplast DNA. The separation of the genusTetradesmus cannot be verified by GC values, becauseS. obliquus andS. (Tetradesmus)wisconsinensis have a similar GC content.S. (Chlorella)ultrasquamata, S. costato-granulatus (sect.Costato-granulati) andS. lunatus are separated clearly from all other species of the genus because of their high GC content.     

Seed oil fatty acid patterns of theAconitum-Delphinium-Helleborus complex (Ranunculaceae)

Many members ofRanunculaceae contain unusual fatty acids in their seed oils. This leads to rather typical genus-specific fatty acid patterns or fingerprints in these seed oils. The members of theDelphinioideae and/orHelleboroideae, however, do not contain highly unusual fatty acids. Nevertheless, their seed oil fatty acid fingerprints are also fairly typical and genus-specific, and the patterns found are rather consistent throughout several species of one genus. It was found that species ofAconitum do not contain fatty acids with 20 carbon atoms.Delphinium, Consolida, Helleborus, Nigella and others do contain C₂₀ fatty acids. In allHelleborus species, for example, there was a consistent C₂₀ fatty acid pattern of 20:020:120:2>20:3. Species ofNigella andGaridella contain high levels,Helleborus low levels, of 20:2n-6 in their seed oils.Delphinium andAconitum both contain low levels of 18:3n-3, whereasHelleborus spp. consistently show high levels of this fatty acid. The genus-specific fatty acid patterns found are discussed, and a correlation with the subfamily and tribe affiliation of the genera investigated here is attempted.