Macroevolutionary relationships of species of Drosophila melanogaster group based on mtDNA sequences

The phylogenetic relationships among the Drosophila melanogaster group species were analyzed using approximately 1700 nucleotide-long sequences of the mitochondrial DNA. Phylogenetic analysis was performed using this region consisting of a part of the cytochrome b (cytb) coding gene, the entire coding sequences of tRNA-Leu, tRNA-Ser and the first subunit of NADH dehydrogenase (NADH1), and a part of the 16S-rRNA gene. The study of these sequences showed that this region of mtDNA is very invariable, as regards with the type of the genes that it contains, as well as the order that they are located on it. The resulting phylogenetic trees reveal a topology that separates the species into three main ancestral lines, leading to the following subgroups: (a) ananassae subgroup, (b) montium subgroup, and (c) melanogaster and Oriental subgroups. The inferred topology complements and generally agrees with previously proposed classifications based on morphological and molecular data.     

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.     

The Heat-Shock Gene hsp83 of Drosophila auraria: Genomic Organization, Nucleotide Sequence, and Long Antiparallel Coupled ORFs (LAC ORFs)

The genomic organization of the hsp83 gene of Drosophila auraria, a far-eastern endemic species belonging to the montium subgroup of the melanogaster species group, is presented here. Based on in situ hybridization on polytene chromosomes, cDNA and genomic clone mapping, nucleotide sequencing, and genomic Southern analysis, hsp83 is shown to be present as a single-copy gene at locus 64B on the 3L chromosome arm in D. auraria. This gene is organized into two exons separated by a 929-bp intron. The first exon represents the mRNA leader sequence and is not translated, while the coding region, having a length of 2,151 bp, is solely included in the second exon. Nucleotide sequence comparisons of D. auraria hsp83 with homologous sequences from other organisms show high conservation of the coding region (88–92% identity) in the genus Drosophila, in addition to the conserved genomic organization of two-exons–one-intron, of comparable size and arrangement. A phylogenetic tree based on the protein sequences of homologous genes from representative organisms is in accord with the accredited phylogenetic position of D. auraria. In the hsp83 gene region, a second case of long antiparallel coupled open reading frames (LAC ORFs) for this species was found. The antiparallel to the hsp83 gene ORF is 1,554 bases long, while the two ORFs overlap has a size of 1,548 bp. The anti-hsp83 ORF does not show significant homology to any known gene sequences. In addition, no similar LAC ORF structures were found in homologous gene regions of other organisms. Received: 18 April 1997 / Accepted: 1 August 1997     

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     

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.     

Identification of One Intron Loss and Phylogenetic Evolution of Dfak Gene in the Drosophila melanogaster Species Group

Intron loss and its evolutionary significance have been noted in Drosophila. The current study provides another example of intron loss within a single-copy Dfak gene in Drosophila. By using polymerase chain reaction (PCR), we amplified about 1.3 kb fragment spanning intron 5–10, located in the position of Tyr kinase (TyK) domain of Dfak gene from Drosophila melanogaster species group, and observed size difference among the amplified DNA fragments from different species. Further sequencing analysis revealed that D. melanogaster and D. simulans deleted an about 60 bp of DNA fragment relative to other 7 Drosophila species, such as D. elegans, D. ficusphila, D. biarmipes, D. takahashii, D. jambulina, D. prostipennis and D. pseudoobscura, and the deleted fragment located precisely in the position of one intron. The data suggested that intron loss might have occurred in the Dfak gene evolutionary process of D. melanogaster and D. simulans of Drosophila melanogaster species group. In addition, the constructed phylogenetic tree based on the Dfak TyK domains clearly revealed the evolutionary relationships between subgroups of Drosophila melanogaster species group, and the intron loss identified from D. melanogaster and D. simulans provides a unique diagnostic tool for taxonomic classification of the melanogaster subgroup from other group of genus Drosophila.     

Evolution of thedec-1 eggshell locus inDrosophila. I. Restriction site mapping and limited sequence comparison in themelanogaster species subgroup

Summary We have analyzed 18 kb of DNA in and upstream of thedefective chorion-1 (dec-1) locus of the eight known species of themelanogaster species subgroup ofDrosophila. The restriction maps ofD. simulans, D. mauritiana, D. sechellia, D. erecta, andD. orena are shown to have basically the restriction map ofD. melanogaster, whereas the maps ofD. teissieri andD. yakuba were more difficult to align. However, the basic amount of DNA and sequence arrangement appear to have been conserved in these species. A small deletion of varying length (65–200 bp) is found in a repeated sequence of the central transcribed region ofD. melanogaster, D. simulans, andD. erecta. Restriction site mapping indicated that thedec-1 gene is highly conserved in themelanogaster species subgroup. However, sequence comparison revealed that the amount of nucleotide and amino acid substitution in the repeated region is much larger than in the 5 translated region. The 5 flanking region showed noticeable restriction site polymorphisms between species. Based on calculations from the restriction maps a dendrogram was derived that supports earlier published phylogenetic relationships within themelanogaster species subgroup except that theerecta-orena pair is placed closer to themelanogaster complex than toD. teissieri andD. yakuba.     

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.     

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.     

Distribution and conservation of the foldback transposable element inDrosophila

Summary Foldback elements are a family of transposable elements described inDrosophila melanogaster. The members of this dispersed repetitive family have terminal inverted repeats that sometimes flank a central region. The inverted repeats of all the family members are homologous.The study of the distribution and conservation of the foldback elements in differentDrosophila species shows that this distribution is different from that of the hybrid dysgenesis systems (PM and IR). Sequences homologous to foldback elements were observed by Southern blots and in situ hybridization in all species of themelanogaster subgroup and in some species of themontium andtakahashii subgroups. The element was probably already present before the radiation of these subgroups. No evidence of horizontal transmission of the foldback element could be observed.     

Photorhabdus temperata subsp. stackebrandtii subsp. nov. (Enterobacteriales: Enterobacteriaceae)

The bacterial symbiont of the entomopathogenic nematode Heterorhabditis bacteriophora strain GPS11 was characterized by 16S rRNA gene sequence and physiological traits. The phylogenetic tree built upon 16S rRNA gene sequences clustered the GPS11 bacterial isolate with Photorhabdus temperata strains which have been previously isolated from Heterorhabditis species. The phylogenetic tree further identified four subgroups in P. temperata, and the relationships among these subgroups were confirmed by gyrase subunit B (gyrB) gene sequence analysis. The subgroup containing the GPS11 bacterial isolate differs from other subgroups in sequences of 16S rRNA and gyrB gene, physiological traits, nematode host species, and geographic origin. Therefore, the subgroup comprising the GPS11 bacterial isolate is proposed here as a new subspecies: Photorhabdus temperata subsp. stackebrandtii subsp. nov. (type strain GPS11). The type strain has been deposited in ATCC and DSMZ collections.     

Functional conservation of a glucose-repressible amylase gene promoter from Drosophila virilis in Drosophila melanogaster

Summary Previous studies have demonstrated that the expression of the -amylase gene is repressed by dietary glucose in Drosophila melanogaster. Here, we show that the -amylase gene of a distantly related species, D. virilis, is also subject to glucose repression. Moreover, the cloned amylase gene of D. virilis is shown to be glucose repressible when it is transiently expressed in D. melanogaster larvae. This cross-species, functional conservation is mediated by a 330-bp promoter region of the D. virilis amylase gene. These results indicate that the promoter elements required for glucose repression are conserved between distantly related Drosophila species. A sequence comparison between the amylase genes of D. virilis and D. melanogaster shows that the promoter sequences diverge to a much greater degree than the coding sequences. The amylase promoters of the two species do, however, share small clusters of sequence similarity, suggesting that these conserved cis-acting elements are sufficient to control the glucose-regulated expression of the amylase gene in the genus Drosophila.Offprint requests to: D.A. Hickey     

Association Between Levels of Coding Sequence Divergence and Gene Misregulation in <Emphasis Type="Italic">Drosophila</Emphasis> Male Hybrids

Previous studies have shown widespread conservation of gene expression levels between species of the Drosophila melanogaster subgroup as well as a positive correlation between coding sequence divergence and expression level divergence between species. Meanwhile, large-scale misregulation of gene expression level has been described in interspecific sterile hybrids between D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. Using data from gene expression analysis involving D. simulans, D. melanogaster, and their hybrids, we observed a significant positive correlation between protein sequence divergence and gene expression differences between hybrids and their parental species. Furthermore, we demonstrate that underexpressed misregulated genes in hybrids are evolving more rapidly at the protein sequence level than nonmisregulated genes or overexpressed misregulated genes, highlighting the possible effects of sexual and natural selection as male-biased genes and nonessential genes are the principal gene categories affected by interspecific hybrid misregulation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Carlo G. Artieri and Wilfried Haerty contributed equally to this publication.     

小麦抗白粉病相关基因GST克隆与表达

以从小麦抗白粉病相关基因差异表达分析中获得的EST-3 (Genbank序列号EX567360)为标签,采用电子克隆的方法对其进行延伸,并对电子克隆结果进行半定量RT-PCR验证,最后对白粉菌不同侵染时间进行了表达分析.经RT-PCR扩增,EST-3表达的带型变化趋势与其在抑制性消减杂交SSH-cDNA的差异显示情况一致,且RT-PCR获得的序列与电子克隆的序列一致性达98％.生物信息学分析表明,该序列是由875 bp核苷酸组成的,具有完整的开放阅读框架,编码蛋白为229个氨基酸,GenBank序列号JK841279,含有一个N端和C端谷胱甘肽硫转移酶结构域,该序列与小麦谷胱甘肽硫转移酶基因(GST)一致性较高,达97％.表达分析结果显示,白粉菌侵染24 h表达受到抑制,48 h开始表达,侵染72 h表达最强,96 h又开始下降,表明GST基因属于白粉菌诱导型相关基因,参与小麦对白粉病的应答反应.     

The sex-lethal gene homologue in Chrysomya rufifacies is highly conserved in sequence and exon-intron organization

A great variety of sex determination mechanisms exists in insect species. In Drosophila melanogaster sex is determined by the ratio between X chromosomes and autosomes, while in the blowfly Chrysomya rufifacies it is maternally determined. A cascade of genes which are involved in sex determination has been identified in D. melanogaster with the Sex-lethal gene (Sxl) as the key gene. We screened genomic libraries of C. rufifacies with a probe of the Sxl gene from D. melanogaster and isolated a genomic region that included most of the homologous gene. DNA- and protein-sequence comparison showed a high percent identity between the Chrysomya and the Drosophila gene. Up to 90% identity of the amino acid sequences was found in the region that contained the RNA-binding domains. The degree of identity is much lower outside of this functionally important region (18% identity). cDNA analysis showed a highly conserved exon-intron structure between the two species, although sex-specific splicing as used in D. melanogaster for the regulation of Sxl activity, could not be detected in C. rufifacies.     

Evolutionary History and Mode of the <Emphasis Type="Italic">amylase</Emphasis> Multigene Family in <Emphasis Type="Italic">Drosophila</Emphasis>

Previous studies indicate that the tandemly repeated members of the amylase (Amy) gene family evolved in a concerted manner in the melanogaster subgroup and in some other species. In this paper, we analyzed all of the 49 active and complete Amy gene sequences in Drosophila, mostly from subgenus Sophophora. Phylogenetic analysis indicated that the two types of diverged Amy genes in the Drosophila montium subgroup and Drosophila ananassae, which are located in distant chromosomal regions from each other, originated independently in different evolutionary lineages of the melanogaster group after the split of the obscura and melanogaster groups. One of the two clusters was lost after duplication in the melanogaster subgroup. Given the time, 24.9 mya, of divergence between the obscura and the melanogaster groups (Russo et al. 1995), the two duplication events were estimated to occur at about 13.96 ± 1.93 and 12.38 ± 1.76 mya in the montium subgroup and D. ananassae, respectively. An accelerated rate of amino acid changes was not observed in either lineage after these gene duplications. However, the G+C contents at the third codon positions (GC3) decreased significantly along one of the two Amy clusters both in the montium subgroup and in D. ananassae right after gene duplication. Furthermore, one of the two types of the Amy genes with a lower GC3 content has lost a specific regulatory element within the montium subgroup species and D. ananassae. While the tandemly repeated members evolved in a concerted manner, the two types of diverged Amy genes in Drosophila experienced frequent gene duplication, gene loss, and divergent evolution following the model of a birth-and-death process.     

Evolution of the LINE-like I element in the Drosophila melanogaster species subgroup

LINE-like retrotransposons, the so-called I elements, control the system of I-R (inducer-reactive) hybrid dysgenesis in Drosophila melanogaster. I elements are present in many Drosophila species. It has been suggested that active, complete I elements, located at different sites on the chromosomes, invaded natural populations of D. melanogaster recently (1920–1970). But old strains lacking active I elements have only defective I elements located in the chromocenter. We have cloned I elements from D. melanogaster and the melanogaster subgroup. In D. melanogaster, the nucleotide sequences of chromocentral I elements differed from those on chromosome arms by as much as 7%. All the I elements of D. mauritiana and D. sechellia are more closely related to the chromosomal I elements of D. melanogaster than to the chromocentral I elements in any species. No sequence difference was observed in the surveyed region between two chromosomal I elements isolated from D. melanogaster and one from D. simulans. These findings strongly support the idea that the defective chromocentral I elements of D. melanogaster originated before the species diverged and the chromosomal I elements were eliminated. The chromosomal I elements reinvaded natural populations of D. melanogaster recently, and were possibly introduced from D. simulans by horizontal transmission.     

Evolutionary pattern of the gtwin retrotransposon in the Drosophila melanogaster subgroup

The gtwin retrotransposon was recently discovered in the Drosophila melanogaster genome and it is evolutionarily closer to gypsy endogenous retrovirus. This study has identified gtwin homologous sequences in the genome of D. simulans, D. sechellia, D. erecta and D. yakuba by performing homology searches against the public genome database of Drosophila species. The phylogenetic analyses of the gtwin env gene sequences of these species have shown some incongruities with the host species phylogeny, suggesting some horizontal transfer events for this retroelement. Moreover, we reported the existence of DNA sequences putatively encoding full-length Env proteins in the genomes of Drosophila species other than D. melanogaster. The results suggest that the gtwin element may be an infectious retrovirus able to invade the genome of new species, supporting the gtwin evolutionary picture shown in this work.     

Mitochondrial gene sequences alone or combined with ITS region sequences provide firm molecular criteria for the classification of Lecanicillium species

The recent revision of Verticillium sect. Prostrata led to the introduction of the genus Lecanicillium, which comprises the majority of the entomopathogenic strains. Sixty-five strains previously classified as Verticillium lecanii or Verticillium sp. from different geographical regions and hosts were examined and their phylogenetic relationships were determined using sequences from three mitochondrial (mt) genes [the small rRNA subunit (rns), the NADH dehydrogenase subunits 1 (nad1) and 3 (nad3)] and the ITS region. In general, single gene phylogenetic trees differentiated and placed the strains examined in well-supported (by BS analysis) groups of L. lecanii, L. longisporum, L. muscarium, and L. nodulosum, although in some cases a few uncertainties still remained. nad1 was the most informative single gene in phylogenetic analyses and was also found to contain group I introns with putative open reading frames (ORFs) encoding for GIY–YIG endonucleases. The combined use of mt gene sequences resolved taxonomic uncertainties arisen from ITS analysis and, alone or in combination with ITS sequences, helped in placing uncharacterised Verticillium lecanii and Verticillium sp. firmly into Lecanicillium species. Combined gene data from all the mt genes and all the mt genes and the ITS region together, were very similar. Furthermore, a relaxed correlation with host specificity—at least for Homoptera—was indicated for the rns and the combined mt gene sequences. Thus, the usefulness of mt gene sequences as a convenient molecular tool in phylogenetic studies of entomopathogenic fungi was demonstrated.     

Phylogenetic position of someChlorella species within the chlorococcales based upon complete small-subunit ribosomal RNA sequences

Summary Complete small-subunit rRNA (16S-like rRNA) coding region sequences were determined for eight species of the Chlorococcales (Chlorophyceae). The genera investigated includePrototheca, Ankistrodesmus, Scenedesmus, and fiveChlorella species. Distance matrix methods were used to infer a phylogenetic tree that describes evolutionary relationships between several plant and green algal groups. The tree exhibits a bifurcation within the Chlorococcales consistent with the division into Oocystaceae and Scenedesmaceae, but three of the fiveChlorella species are more similar to other algae than toChlorella vulgaris. All of the sequences contain primary and secondary structural features that are characteristic of 16S-like rRNAs of chlorophytes and higher plants.Anikstrodesmus stipitatus, however, contains a 394-bp group I intervening sequence in its 16S-like rRNA coding region.