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.     

Phylogenetic relationships of Drosophila melanogaster species group deduced from spacer regions of histone gene H2A-H2B

Nucleotide sequences of the spacer region of the histone gene H2A-H2B from 36 species of Drosophila melanogaster species group were determined. The phylogenetic trees were reconstructed with maximum parsimony, maximum likelihood, and Bayesian methods by using Drosophila pseudoobscura as the out group. Our results show that the melanogaster species group clustered in three main lineages: (1). montium subgroup; (2). ananassae subgroup; and (3). the seven oriental subgroups, among which the montium subgroup diverged first. In the third main lineage, suzukii and takahashii subgroups formed a clade, while eugracilis, melanogaster, elegans, ficusphila, and rhopaloa subgroups formed another clade. The bootstrap values at subgroup levels are high. The phylogenetic relationships of these species subgroups derived from our data are very different from those based on some other DNA data and morphology data.     

以ND4L和ND4基因为标记探讨黑腹果蝇种组的系统发育关系

多年来的形态学、染色体组学以及DNA序列几个方面的研究均没有很好地阐明黑腹果蝇种组内的系统发育关系。本实验测定了33个样品的ND4和31个样品的ND4L基因序列,以D．obscuroides为外群,用最大简约法和Bayesian法分别构建进化树。结果表明两种方法构建的拓扑结构一致,而且大部分支系的支持率较高。整个黑腹果蝇种组分成三大谱系：1）montium种亚组;2）ananssae种亚组;3）Oriental种亚组（melanogaster、ficsphila、eugracilis、elegans、suzukii、takahashii）。montium是最早分化的种亚组。在第三谱系中,melanogaster分化得最早;然后依次是ficsphila,eugracilis,elegans;suzukii与takahashii为姐妹种亚组,最后分化。     

Functional constraints of the Cu,Zn superoxide dismutase in species of the Drosophila melanogaster subgroup and phylogenetic analysis

The phylogenetic relationships among the Drosophila melanogaster subgroup species were analyzed using approximately 1550-nucleotide-long sequences of the Cu,Zn SOD gene. Phylogenetic analysis was performed using separately the whole region and the intron sequences of the gene. The resulting phylogenetic trees reveal virtually the same topology, separating the species into distinct clusters. The inferred topology generally agrees with previously proposed classifications based on morphological and molecular data. The amino acid sequences of the Cu,Zn SOD of the D. melanogaster subgroup species reveal a high-conservation pattern. Only 3.9% of the total amino acid sites are variable, and none affects the major structural elements. Comparison of the Drosophila Cu,Zn SOD amino acid sequences with the Cu,Zn SOD of Bos taurus and Xenopus laevis (whose three-dimensional structure has been elucidated) reveals conservation of all the protein's functionally important amino acids and no substitutions that dramatically change the charge or the polarity of the amino acids.     

Molecular evolution of the 5'-flanking regions of the duplicated Amy genes in Drosophila melanogaster species subgroup

The nucleotide sequences of the 5'-flanking regions of the duplicated Amy genes in eight sibling species belonging to the melanogaster species subgroup are analyzed. In Drosophila melanogaster, a region of about 450 bp immediately upstream of the translation initiation site of the two paralogous genes (the proximal and distal genes) has sequence similarities. However, we could not detect any significant sequence similarity in the region more upstream than -450. This result indicates that the coding regions of the ancestral Amy gene were duplicated together with 450 bp of the 5'-flanking region as one unit. Multiple alignment of these 450-bp sequences in the proximal and distal genes of all eight species revealed a mosaic pattern of highly conserved and divergent regions. The conserved regions included almost all the putative regulatory elements identified in previous analyses of the sequences. A phylogenetic analysis of the aligned sequences shows that these 450-bp sequences are clustered into the proximal and the distal groups. As a whole, the divergence between groups in this region is very large in contrast to that in the coding regions. Based on the divergence between groups, the 450-bp region is divided into two subregions. We found that the ratios of the divergence between groups to that within groups differ in the two subregions. From these observations, we discuss a possibility of positive selection acting on the subregion immediately upstream of the Amy coding region to cause divergence of regulatory elements of the paralogous genes.      

The 5S ribosomal genes in the Drosophila melanogaster species subgroup. Nucleotide sequence of a 5S unit from Drosophila simulans and Drosophila teissieri

The 5S genes of the eight species of the D. melanogaster subgroup have been mapped. The spacers, in contrast with coding regions, differ markedly between most species. One 5S gene unit has been sequenced for both D. simulans and D. teissieri. The mature 5S RNA region in these two species is identical to the corresponding region of D. melanogaster. Only 5 nucleotide variations occur between the D. melanogaster and D. simulans 5S gene spacers. The spacer in D. teissieri is very different. Only two segments, located one at each side of the coding region, are clearly homologous to corresponding sequences of D. melanogaster and D. simulans.     

中国特有果蝇 Drosophila curviceps 种亚组在Drosophila immigrans种组中的种系发生地位

果蝇immigrans种组中的curviceps种亚组是1992年新建立的中国特有果蝇类群。该种亚组中的物种主要分布在中国大陆和台湾。目前除了形态学水平的研究外,还没有其他证据支持建立该种亚组的合理性及其起源和种系发生地位。为了在DNA分子水平上探讨果蝇curviceps种亚组在果蝇immigrans种组中的种系发生地位,从而为今后更深入地研究中国特有果蝇,甚至为果蝇亚属的进化遗传学提供理论依据,测定了immigrans种组5个种亚组（nasuta、immigrans、hypocausta、quadrilineata、curviceps)中12个代表物种的rDNA的ITS1和部分Adh基因的序列。其中ITS1序列的长度为513－587bp,共有191个信息位点;Adh基因片段的长度在714－747bp之间,共99个信息位点。考虑到单个分子提供的信息较少,将两个分子的序列综合起来,组成一个较长的复合序列。分别根据ITS1,Adh和两个分子的复合序列排比（Alignment)结果,和最大简约法和邻接法构建分子系统树,其中根据复合序列构建的系统树与形态学研究结果最为一致。分子树显示curviceps种亚组的特种确定单独形成一个分枝,为种亚组级的分类阶元,支持了形态学将其建立为一个新种亚组。根据Kimura距离,估算了复合分子的替换速率约为每百万年1．48％,进而计算出5个种亚组的分 歧年代。结合各物种的地理分布,推测了immigrans种组的进化历史：curviceps种亚组与quadrilineata种亚组的亲缘关系最近,主要分布在中国南部的温带地区。它们之间的分歧时间大约为3．4百万年,是最年轻的两个种亚组。主要分布在苏门答腊及附近的热带地区的hypocausta种亚组的物种是最早分化出来的,与其他种亚组的分歧时间约为9．2百万年。该结果与形态学和生物地理学研究相吻合。值得一提是的,目前归属仍存在争议的物种D.neohypocausta,在分子系统树中与hypocausta种亚组的物种相距较远,而与immiagrasn种亚组的关系较近,但分枝置信度较低（＜50％）。由于还缺乏其他方面的证据,因此D.neohypocausta的归属有待今后的研究来作定论。     

Evolution of the Transposable Element Mariner in Drosophila Species

The distribution of the transposable element mariner was examined in the genus Drosophila. Among the eight species comprising the melanogaster species subgroup, the element is present in D. mauritiana, D. simulans, D. sechellia, D. yakuba and D. teissieri, but it is absent in D. melanogaster, D. erecta and D. orena. Multiple copies of mariner were sequenced from each species in which the element occurs. The inferred phylogeny of the elements and the pattern of divergence were examined in order to evaluate whether horizontal transfer among species or stochastic loss could better account for the discontinuous distribution of the element among the species. The data suggest that the element was present in the ancestral species before the melanogaster subgroup diverged and was lost in the lineage leading to D. melanogaster and the lineage leading to D. erecta and D. orena. This inference is consistent with the finding that mariner also occurs in members of several other species subgroups within the overall melanogaster species group. Within the melanogaster species subgroup, the average divergence of mariner copies between species was lower than the coding region of the alcohol dehydrogenase (Adh) gene. However, the divergence of mariner elements within species was as great as that observed for Adh. We conclude that the relative sequence homogeneity of mariner elements within species is more likely a result of rapid amplification of a few ancestral elements than of concerted evolution. The mariner element may also have had unequal mutation rates in different lineages.     

Evolutionary conservation of the chromosomal configuration and regulation of amylase genes among eight species of the Drosophila melanogaster species subgroup

Nuclear DNA was extracted from each of the eight species comprising the Drosophila melanogaster species subgroup. Southern hybridization of this DNA by using a molecular probe specific for the alpha-amylase coding region showed that the duplicated structure of the amylase locus, first found in D. melanogaster, is conserved among all species of the melanogaster subgroup. Evidence is also presented for the concerted evolution of the duplicated genes within each species. In addition, it is shown that the glucose repression of amylase gene expression, which has been extensively studied in D. melanogaster, is not confined to this species but occurs in all eight members of the species subgroup. Thus, both the duplicated gene structure and the glucose repression of Drosophila amylase gene activity are stable over extended periods of evolutionary time.      

Divergence and heterogeneity of the histone gene repeating units in the Drosophila melanogaster species subgroup

The repeating units of the histone gene cluster containing the H1, H2A, H2B and H4 genes were amplified by PCR from the Drosophila melanogaster species subgroup, i.e., D. yakuba, D. erecta, D. sechellia, D. mauritiana, D. teissieri and D. orena. The PCR products were cloned and their nucleotide sequences of about 4.6-4.8kbp were determined to elucidate the mechanism of molecular evolution of the histone gene family. The heterogeneity among the histone gene repeating units was 0.6% and 0.7% for D. yakuba and D. sechellia, respectively, indicating the same level of heterogeneity as in the H3 gene region of D. melanogaster. Divergence of the genes among species even in the most closely related ones was much greater than the heterogeneity among family members, indicating a concerted mode of evolution for the histone gene repeating units. Among the species in the D. melanogaster species subgroup, the histone gene regions as well as 3rd codon position of the coding region showed nearly the same GC contents. These results suggested that the previous conclusion on analysis of the H3 gene regions, the gene family evolution in a concerted fashion, holds true for the whole histone gene repeating unit.     

Relationships in the Drosophila obscura species group, inferred from mitochondrial cytochrome oxidase II sequences

We compare the sequences for the mitochondrial cytochrome oxidase II gene of 13 species of the Drosophila obscura group. The survey includes six members of the D. affinis subgroup, four of the D. pseudoobscura subgroup, and three of the D. obscura subgroup. In all species, the gene is 688 nucleotides in length, encoding a protein of 229 amino acids plus the first position T of the stop codon. The sequences show the typical high-transition bias for closely related species, but that bias is essentially eliminated for species pairs of > 5% sequence divergence. The phylogenetic relationships in the species group are inferred using both neighbor-joining and maximum parsimony. The two procedures give comparable results, showing that the D. affinis and D. pseudoobscura subgroups are monophyletic groupings that appear to have closer affinities to one another than either has to the D. obscura subgroup. We use transversion distances to estimate times of divergence, on the basis of three different estimates of the time of separation of the D. obscura species group from the D. melanogaster group. If that event occurred 35 Mya, then we can estimate the origin of the nearctic forms at approximately 22 Mya and the separation of the D. affinis and D. pseudoobscura subgroups at approximately 17 Mya.      

Evolution of gypsy endogenous retrovirus in the Drosophila obscura species group

The Ty3/gypsy family of retroelements is closely related to retroviruses, and some of their members have an open reading frame resembling the retroviral gene env. Sequences homologous to the gypsy element from Drosophila melanogaster are widely distributed among Drosophila species. In this work, we report a phylogenetic study based mainly on the analysis of the 5' region of the env gene from several species of the obscura group, and also from sequences already reported of D. melanogaster, Drosophila virilis, and Drosophila hydei. Our results indicate that the gypsy elements from species of the obscura group constitute a monophyletic group which has strongly diverged from the prototypic D. melanogaster gypsy element. Phylogenetic relationships between gypsy sequences from the obscura group are consistent with those of their hosts, indicating vertical transmission. However, D. hydei and D. virilis gypsy sequences are closely related to those of the affinis subgroup, which could be indicative of horizontal transmission.     

A phylogeny of Drosophilidae using the Amyrel gene: questioning the Drosophila melanogaster species group boundaries

In this study, the phylogenetic relationships of 164 species of the family Drosophilidae are discussed, using the Amyrel gene, a member of the α -amylase multigene family. This study focuses on numerous species groups in the subgenera Sophophora and Drosophila of the genus Drosophila but also includes other closely related genera. Nucleotide data were analysed by several methods: maximum parsimony, neighbour joining, maximum likelihood and Bayesian inference. Heterogeneity of base composition (mainly low GC contents in the species groups willistoni and saltans ) has been addressed. In all analyses, the genus Drosophila appeared paraphyletic. The subgenus Sophophora clearly appeared to be a monophyletic group, showing well-resolved clades, with the Neotropical groups arising in a basal position. Here, it is proposed to raise the species subgroups ananassae and montium to the rank of species group, and to restrict the melanogaster species group to the melanogaster subgroup plus the 'Oriental' subgroups, among which the suzukii subgroup is polyphyletic. Some related genera such as Zaprionus , Liodrosophila , Scaptomyza and Hirtodrosophila are clustered with, or inside the subgenus Drosophila , which is therefore paraphyletic and should be reviewed.     

Phylogenetic incongruence in the Drosophila melanogaster species group

Drosophila melanogaster and its close relatives are used extensively in comparative biology. Despite the importance of phylogenetic information for such studies, relationships between some melanogaster species group members are unclear due to conflicting phylogenetic signals at different loci. In this study, we use twelve nuclear loci (eleven coding and one non-coding) to assess the degree of phylogenetic incongruence in this model system. We focus on two nodes: (1) the node joining the Drosophila erecta-Drosophila orena, Drosophila melanogaster-Drosophila simulans, and Drosophila yakuba-Drosophila teissieri lineages, and (2) the node joining the lineages leading to the melanogaster, takahashii, and eugracilis subgroups. We find limited evidence for incongruence at the first node; our data, as well as those of several previous studies, strongly support monophyly of a clade consisting of D. erecta-D. orena and D. yakuba-D. teissieri. By contrast, using likelihood based tests of congruence, we find robust evidence for topological incongruence at the second node. Different loci support different relationships among the melanogaster, takahashii, and eugracilis subgroups, and the observed incongruence is not easily attributable to homoplasy, non-equilibrium base composition, or positive selection on a subset of loci. We argue that lineage sorting in the common ancestor of these three subgroups is the most plausible explanation for our observations. Such lineage sorting may lead to biased estimation of tree topology and evolutionary rates, and may confound inferences of positive selection.     

Changes in the recombinational environment affect divergence in the yellow gene of Drosophila

The complete coding region of the yellow (y) gene was sequenced in different Drosophila species. In the species of the melanogaster subgroup (D. melanogaster, D. simulans, D. mauritiana, D. yakuba, and D. erecta), this gene is located at the tip of the X chromosome in a region with a strong reduction in recombination rate. In contrast, in D. ananassae (included in the ananassae subgroup of the melanogaster group) and in the obscura group species (D. subobscura, D. madeirensis, D. guanche, and D. pseudoobscura), the y gene is located in regions with normal recombination rates. As predicted by the hitchhiking and background selection models, this change in the recombinational environment affected synonymous divergence in the y-gene-coding region. Estimates of the number of synonymous substitutions per site were much lower between the obscura group species and D. ananassae than between the species of the obscura group and the melanogaster subgroup. In fact, a highly significant increase in the rate of synonymous substitution was detected in all lineages leading to the species of the melanogaster subgroup relative to the D. ananassae lineage. This increase can be explained by a higher fixation rate of mutations from preferred to unpreferred codons (slightly deleterious mutations). The lower codon bias detected in all species of the melanogaster subgroup relative to D. ananassae (or to the obscura group species) would be consistent with this proposal. Therefore, at least in Drosophila, changes in the recombination rate in different lineages might cause deviations of the molecular-clock hypothesis and contribute to the overdispersion of the rate of synonymous substitution. In contrast, the change in the recombinational environment of the y gene has no detectable effect on the rate of amino acid replacement in the Yellow protein.     

Molecular evolution of the histone 3 multigene family in the Drosophila melanogaster species subgroup

Molecular evolution of the histone multigene family was studied by cloning and sequencing regions of the histone 3 gene in the Drosophila melanogaster species subgroup. Analysis of the nucleotide substitution pattern showed that in the coding region synonymous changes occurred more frequently to A or T in contrast to the GC-rich base composition, while in the 3' region the nucleotide substitutions were most likely in equilibrium. These results suggested that the base composition at the third codon position of the H3 gene, i.e., codon usage, has been changing to A or T in the Drosophila melanogaster species subgroup.     

Phylogeny of a paradigm lineage: the Drosophila melanogaster species group (Diptera: Drosophilidae)

Although Drosophila melanogaster is a paradigm eukaryote for biology, relationships of this species and the other 174 species in the melanogaster species group are poorly explored and ambiguous. Gene regions of Cytochrome oxidase II (mt:CoII ), Alcohol dehydrogenase ( Adh ) and hunchback ( hb ) were sequenced and analysed phylogenetically to test prior hypotheses of relationships for the group based on chromosomes, morphology, and 28S rRNA gene sequences. A simultaneous cladistic analysis of the three newly sequenced gene regions produced a single well-resolved phylogeny for 49 exemplar species representing eight subgroups. Monophyly of each of the ananassae , melanogaster , montium , and takahashii subgroups is supported; the suzukii subgroup is polyphyletic. This phylogeny is consistent with variation in significant morphological structures, such as the male sex comb on the fore tarsus. The broad range of morphological variation among these species is interpreted and the applicability to evolution and developmental investigations is discussed. This phylogeny facilitates comparative investigations, such as gene family evolution, transposable element transmission, and evolution of morphological structures. © 2002 The Linnean Society of London, Biological Journal of the Linnean Society , 2002, 76 , 21–37.     

The structure of a single unit of ribosomal RNA gene (rDNA) including intergenic subrepeats in the Australian bulldog ant Myrmecia croslandi (Hymenoptera: Formicidae)

A complete single unit of a ribosomal RNA gene (rDNA) of M. croslandi was sequenced. The ends of the 18S, 5.8S and 28S rRNA genes were determined by using the sequences of D. melanogaster rDNAs as references. Each of the tandemly repeated rDNA units consists of coding and non-coding regions whose arrangement is the same as that of D. melanogaster rDNA. The intergenic spacer (IGS) contains, as in other species, a region with subrepeats, of which the sequences are different from those previously reported in other insect species. The length of IGSs was estimated to be 7-12 kb by genomic Southern hybridization, showing that an rDNA repeating unit of M. croslandi is 14-19 kb-long. The sequences of the coding regions are highly conserved, whereas IGS and ITS (internal transcribed spacer) sequences are not. We obtained clones with insertions of various sizes of R2 elements, the target sequence of which was found in the 28S rRNA coding region. A short segment in the IGS that follows the 3' end of the 28S rRNA gene was predicted to form a secondary structure with long stems.