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.     

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.     

Basal relationships in the Drosophila melanogaster species group

The Drosophila melanogaster species group is a popular model for evolutionary studies due to its morphological and ecological diversity and its inclusion of the model species D. melanogaster. However, phylogenetic relationships among major lineages within this species group remain controversial. In this report, the phylogeny of 10 species representing each of the well-supported monophyletic clades in the melanogaster group was studied using the sequences of 14 loci that together comprise 9493 nucleotide positions. Combined Bayesian analysis using gene-specific substitution models produced a 100% credible set of two trees. In the strict consensus of these trees, the ananassae subgroup branches first in the melanogaster species group, followed by the montium subgroup. The remaining lineages form a monophyletic clade in which D. ficusphila and D. elegans branch first, followed by D. biarmipes, D. eugracilis, and the melanogaster subgroup. This strongly supported phylogeny resolves most basal relationships in the melanogaster species group, and provides a framework that can be extended in the future to encompass more species.     

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.     

Phylogenetic relationships and climatic adaptations in the Drosophila takahashii and montium species subgroups

We analyze phylogenetic relationships among temperate, subtropical highland, and subtropical lowland species of the Drosophila takahashii and montium species subgroups based on sequence data of COI and Gpdh genes and discuss the evolution of temperate species in these subgroups with reference to their climatic adaptations. In the takahashii subgroup, D. lutescens (the temperate species) branched off first in the tree based on the combined data set, but D. prostipennis (the subtropical highland species) branched off first in the trees based on single genes. Thus, phylogenetic relationships in this subgroup are still ambiguous. In the montium subgroup, the cool-temperate species are phylogenetically close to the warm-temperate species, and these cool- and warm-temperate species form a cluster with the subtropical highland species. This suggests that perhaps the cool-temperate species derived from the warm-temperate species and the warm-temperate species derived from the subtropical highland species. In comparison with the subtropical lowland species, the subtropical highland species may be better able to colonize temperate areas since, as in the temperate species, they have an ability to develop their ovaries at moderately low temperature. However, the subtropical highland species, as well as the subtropical lowland species, were much less cold tolerant than the temperate species. Therefore, considerable genetic reformation would be required for both the subtropical highland and the subtropical lowland species to adapt to temperate climates.     

The phylogeny of the subgroups within the melanogaster species group: likelihood tests on COI and COII sequences and a Bayesian estimate of phylogeny

The relationships among the majority of the subgroups in the Drosophila melanogaster species group remain unresolved. We present a 2223basepair dataset for mitochondrial cytochrome oxidase I and cytochrome oxidase II for 43 species (including new data from 11 species), sampled to include the major subgroups. After a brief review of competing hypotheses for the ananassae, montium, suzukii, and takahashii subgroups, we combine the two genes based on a new use of the SH test and present KH and SH likelihood comparisons (Kishino and Hasegawa, 1989. J. Mol. Evol. 29, 170-179; Shimodaira and Hasegawa, 1999) to test the monophyly and placement of these subgroups within the larger species group. Although we find insignificant differences between the two suggested placements for the ananassae subgroup, the ananassae is sister to the rest of the subgroups in the melanogaster species group in every investigation. For the takahashii subgroup, although we cannot reject monophyly, the species are so closely related to the suzukii subgroup for these data that the two subgroups often form one clade. Finally, we present a Bayesian estimate of the phylogeny for both genes combined, utilizing a recently published method that allows for different models of evolution for different sites.     

线粒体ND4-ND4L基因在黑腹果蝇种组中的进化特征

本实验对黑腹果蝇种组(melanogaster species group)中8个种亚组33个样品两个线粒体基因ND4和ND4L进行了测序,并分析了ND4基因的序列差异和碱基替换特点,发现近缘物种中存在很明显的转换倾向,而在远缘物种中由于重复替换导致转换数处于饱和状态,我们的实验数据证实了线粒体基因较核基因有较快的进化速度。最后根据D．melanogaster与D．yakuba的遗传距离推算了8个种亚组的分化时间,ananassae种亚组最先分化,然后依次是montium,melanogaster,ficsphila,eugracilis,elegans,suzukii和takahashii最后分化。     

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

黑腹果蝇种组的核型与进化研究

观察了国内黑腹果蝇种组34种果蝇的有丝分裂中期核型,其中首次描述了一些新核型。系统地分析了黑腹果蝇种组8个种亚组之间的核型进化关系及种间亲缘关系。结果是:elegans种亚组的核型为A型;eugracilis、melanogaster和ficusphila种亚组的核型为C型;takahashii和suzukii种亚组的核型为C型和D型;montium种亚组的核型为B、C、C’、D、D’、和E型;ananassae种亚组的核型为F、G和H型。从核型分化的角度可以将黑腹果蝇种组分为5个谱系:elegans,eugracilis-melanogaster-ficusphila,takkahashii-suzukii,montium,ananassae。这与2004年Yang等的观点基本一致,正好从核型进化的角度验证了Yang通过DNA序列分析所得到的结果。差别只在于elegans种亚组,作者把它单独列为一支,认为是祖先种亚组。通过选取同一种果蝇的几个不同地域单雌系的核型分析,结果表明:同一种果蝇的核型存在地域差异。这种差异可能是由于不同生境造成,也可能是本身进化程度的差异,或是两种因素相互作用的结果。     

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.     

Ancestral polymorphism and recent invasion of transposable elements in Drosophila species

ABSTRACT: BACKGROUND: During the evolutionary history of transposable elements, some processes, such as ancestral polymorphisms and horizontal transfer of sequences between species, can produce incongruences in phylogenies. We investigated the evolutionary history of the transposable elements Bari and 412 in the sequenced genomes of the Drosophila melanogaster group and in the sibling species D. melanogaster and D. simulans using traditional phylogenetic and network approaches. RESULTS: The maximum likelihood (ML) phylogenetic analyses revealed incongruences and unresolved relationships for both the Bari and 412 elements. The DNA transposon Bari within the D. ananassae genome is more closely related to the element of the melanogaster complex than to the sequence in D. erecta, which is inconsistent with the species phylogeny. Divergence analysis and the comparison of the rate of synonymous substitutions per synonymous site of the Bari and host gene sequences explain the incongruence as an ancestral polymorphism inherited stochastically by the derived species. Unresolved relationships were observed in the ML phylogeny of both elements involving D. melanogaster, D. simulans and D. sechellia. A network approach was used to attempt to resolve these relationships. The resulting tree suggests recent transfers of both elements between D. melanogaster and D. simulans. The divergence values of the elements between these species support this conclusion. CONCLUSIONS: We showed that an ancestral polymorphism and recent invasion of genomes due to introgression or horizontal transfer between species occurred during the evolutionary history of the Bari and 412 elements in the melanogaster group. These invasions likely occurred in Africa during the Pleistocene, before the worldwide expansion of D. melanogaster and D. simulans.     

Spontaneous mutations affecting glycerol-3-phosphate dehydrogenase enzyme activity in Drosophila melanogaster.

Significant genetic variance in glycerol-3-phosphate dehydrogenase (GPDH) activity was observed between chromosome lines of Drosophila melanogaster that had each accumulated spontaneous mutations for approximately 300 generations. No restriction map variation was found in a 26-kb region surrounding the entire Gpdh gene. The restriction analysis used is capable of detecting insertions/deletions larger than 0.05 kb. The survey would also detect chromosomal recombinations that include the entire Gpdh coding region. Therefore, if the spontaneous mutations that affected the enzyme activity are located inside the Gpdh gene region, then they are base pair substitutions or structural changes that are smaller than the limit in resolution described above.     

Nucleotide Variation at the Gpdh Locus in the Genus Drosophila

The Gpdh locus was sequenced in a broad range of Drosophila species. In contrast to the extreme evolutionary constraint seen at the amino acid level, the synonymous sites evolve at rates comparable to those of other genes. Gpdh nucleotide sequences were used to infer a phylogenetic tree, and the relationships among the species of the obscura group were examined in detail. A survey of nucleotide polymorphism within D. pseudoobscura revealed no amino acid variation in this species. Applying a modified McDonald-Kreitman test, the amino acid divergence between species in the obscura group does not appear to be excessive, implying that drift is adequate to explain the patterns of amino acid change at this locus. In addition, the level of polymorphism at the Gpdh locus in D. pseudoobscura is comparable to that found at other loci, as determined by a Hudson-Kreitman-Aguade test. Thus, the pattern of nucleotide variation within and between species at the Gpdh locus is consistent with a neutral model.     

Evolution of the chromosomal location of rDNA genes in two Drosophila species subgroups: ananassae and melanogaster

The evolution of the chromosomal location of ribosomal RNA gene clusters and the organization of heterochromatin in the Drosophila melanogaster group were investigated using fluorescence in situ hybridization and DAPI staining to mitotic chromosomes. The investigation of 18 species (11 of which were being examined for the first time) belonging to the melanogaster and ananassae subgroups suggests that the ancestral configuration consists of one nucleolus organizer (NOR) on each sex chromosome. This pattern, which is conserved throughout the melanogaster subgroup, except in D. simulans and D. sechellia, was observed only in the ercepeae complex within the ananassae subgroup. Both sex-linked NORs must have been lost in the lineage leading to D. varians and in the ananassae and bipectinata complexes, whereas new sites, characterized by intra-species variation in hybridization signal size, appeared on the fourth chromosome related to heterochromatic rearrangements. Nucleolar material is thought to be required for sex chromosome pairing and disjunction in a variety of organisms including Drosophila. Thus, either remnant sequences, possibly intergenic spacer repeats, are still present in the sex chromosomes which have lost their NORs (as observed in D. simulans and D. sechellia), or an alternative mechanism has evolved.     

PHYLOGENETIC RELATIONSHIPS WITHIN THE COLONIAL VOLVOCALES (CHLOROPHYTA) INFERRED FROM rbcL GENE SEQUENCE DATA

The chloroplast-encoded large subunit of the ribulose-1, 5-bisphosphate carboxylase / oxygenase (rbcL) gene was sequenced from 20 species of the colonial Volvocales (the Volvacaceae, Goniaceae, and Tetrabaenaceae) in order to elucidate phylogenetic relationships within the colonial Volvocales. Eleven hundred twenty-eight base pairs in the coding regions of the (rbcL) gene were analyzed by the neighbor-joining (NJ) method using three kinds of distance estimations, as well as by the maximum parsimony (MP) method. A large group comprising all the anisogamous and oogamous volvocacean species was resolved in the MP tree as well as in the NJ trees based on overall and synonymous substitutions. In all the trees constructed, Basichlamys and Tetrabaena (Tetrabaenaceae) constituted a very robust phylogenetic group. Although not supported by high bootstrap values, the MP tree and the NJ tree based on nonsynonymous substitutions indicated that the Tetrabaenaceae is the sister group to the large group comprising the Volvocaceae and the Goniaceae. In addition, the present analysis strongly suggested that Pandorina and Astrephomene are monophyletic genera whereas Eudorina is nonmonophyletic. These results are essentially consistent with the results of the recent cladistic analyses of morphological data. However, the monophyly of the Volvocaceae previously supported by four morphological synapomorphies is found only in the NJ tree based on nonsynonymous substitutions (with very low bootstrap values). The genus Volvox was clearly resolved as a polyphyletic group with V. rousseletii Pocock separated from other species of Volvox in the rbcL gene comparisons, although this genus represents a monophyletic group in the previous morphological analyses. Furthermore, none of the rbcL gene trees supported the monophyly of the Goniaceae; Astrephomene was placed in various phylogenetic positions .     

Evolution of the AMP-forming acetyl-CoA synthetase gene in the Drosophilidae family

Analysis of the AMP-forming ACS gene was performed in 12 species of the Drosophilidae family. Systematically four introns, aligned at the same positions, were detected, but none of them showed a position similar to those known for species outside the Drosophilidae family. The average length of introns varied from 63 to 75 bp but in two species Drosophila takahashii and D. kikkawai the length of the second intron was 343 and 210 bp, respectively. In coding regions, about 80% of the third codon positions were substituted while first and second positions showed, respectively, 14% and 6% substitutions. Interestingly, the divergence observed at the protein level between species was very low. The phylogenetic tree based on the DNA sequences of the exons was mainly in agreement with taxonomic classification and previous molecular phylogenies except for D. ananassae, which appeared more closely related to D. subobscura and D. funebris than to the species of the melanogaster group.     

Resolution of phylogenetic relationships among recently evolved species as a function of amount of DNA sequence: an empirical study based on woodpeckers (Aves: Picidae)

Synonymous substitutions in the 13 mitochondrial encoded protein genes form a large pool of characters that should approach the ideal for phylogenetic analysis of being independently and identically distributed. Pooling sequences from multiple mitochondrial protein-coding genes should result in statistically more powerful estimates of relationships among species that diverged sufficiently recently that most nucleotide substitutions are synonymous. Cytochrome oxidase I (COI) was sequenced for woodpecker species for which cytochrome b (cyt b) sequences were available. A pairing-design test based on the normal distribution indicated that cyt b evolves more rapidly than COI when all nucleotides are compared but their rates are equal for synonymous substitutions. Nearly all of the phylogenetically informative substitutions among woodpeckers are synonymous. Statistical support for relationships, as measured by bootstrap proportions, increased as the number of nucleotides increased from 1047 (cyt b) to 1512 (COI) to 2559 nucleotides (aggregate data set). Pseudo-bootstrap replicates showed the same trend and increasing the amount of sequence beyond the actual length of 2559 nucleotides to 5120 (2x) resulted in stronger bootstrap support, even though the amount of phylogenetic information was the same. However, the amount of sequence required to resolve an internode depends on the length of the internode and its depth in the phylogeny.     

Molecular phylogeny of glass sponges (Porifera, Hexactinellida): increased taxon sampling and inclusion of the mitochondrial protein-coding gene, cytochrome oxidase subunit I

Marine sponges of the class Hexactinellida (glass sponges) are among the most understudied groups of Porifera, and molecular approaches to investigating their evolution have only recently emerged. Although these first results appeared reliable as they largely corroborated morphology-based hypotheses, they were almost exclusively based on ribosomal RNA genes (rDNA) and should, therefore, be further tested with independent types of genetic data, such as protein-coding genes. To this end, we established the mitochondrial-encoded cytochrome oxidase subunit I gene (COI) as an additional marker, and conducted phylogenetic analyses on DNA- and amino-acid level, as well as a supermatrix analysis based on combined COI DNA and rDNA alignments. Furthermore, we increased taxon sampling compared to previous studies by adding seven additional species. The COI-based phylogenies were largely congruent with the rDNA-based phylogeny but suffered from poor bootstrap support for many nodes. However, addition of the COI sequences to the rDNA data set increased resolution of the overall molecular phylogeny. Thus, although obtaining COI sequences from glass sponges turned out to be quite challenging, this gene appears to be a valuable supplement to rDNA data for molecular evolutionary studies of this group. Some implications of our extended phylogeny for the evolution and systematics of Hexactinellida are discussed.