Phylogenetic analysis of the repleta species group of the genus Drosophila using multiple sources of characters

The species in the repleta group of the genus Drosophila have been placed into five subgroups-the mulleri, hydei, mercatorum, repleta, and fasciola subgroups. Each subgroup has been further subdivided into complexes and clusters. Extensive morphological and cytological analyses of the members of this species group have formed the foundation for the proposed relationships among the members of the repleta species group. Fifty-four taxa, including 46 taxa belonging to the repleta species group, were sequenced for fragments of four genes-16S ribosomal DNA (16S), cytochrome oxidase II (COII), and nitrogen dehydrogenase 1 (ND1) of the mitochondrial genome and a region of the hunchback (hb) nuclear gene. We also generated a partial data set of elongation factor 1-alpha (Ef1alpha) sequences for a subset of taxa. Our analysis used both DNA characters and chromosomal inversion data. The phylogenetic hypothesis we obtained supports many of the traditionally accepted clades within the mulleri subgroup, but the monophyly of taxonomic groups outside of this subgroup appears not to be supported. Phylogenetic analysis revealed one well-supported, highly resolved clade that consists of closely related members of the mulleri and buzzatii complexes. The remaining taxa, a wide assortment of taxonomic groups, ranging from members of other species groups to members of several subgroups and members of three species complexes from the mulleri subgroup are found in poorly supported arrangements at the base of the tree.     

The Evolutionary History of DROSOPHILA BUZZATII. III. Cytogenetic Relationships between Two Sibling Species of the Buzzatii Cluster

Drosophila buzzatii has been found sympatric in Argentina with a closely-related sibling species, D. serido. The biogeographical, reproductive and chromosomal data allow us to combine these species into an evolutionary unit, the buzzatii cluster. Salivary gland chromosomes also have been used to determine their phylogenetic relationships with other closely related species, showing that the buzzatii cluster species share two inversions—2d² and 2s⁶—with the species of the martensis cluster. Both clusters arose from South American populations of the ancestor of the mulleri complex, and we propose to include D. buzzatii and D. serido in the mulleri complex of the repleta group.     

DNA Sequence Comparison among Closely Related Drosophila Species in the MULLERI Complex

DNA hybridization was used to establish DNA sequence relationships among seven Drosophila species. Single-copy DNA was isolated from four species within the Drosophila mulleri complex, D. mojavensis, D. arizonensis, D. ritae and D. starmeri. These single-copy DNAs were used as tracers to be hybridized with each other and one additional member of the mulleri complex, D. aldrichi, a member of a closely related complex, D. hydei, and a distantly related species, D. melanogaster. Two methods have been used to determine the relatedness between these species: (1) the extent of duplex formed as measured by binding to hydroxyapatite and (2) the thermal stability of the duplexed DNA. Moderately repetitive DNA was purified from these species and used similarly to determine the divergence of this family of sequences. The rate of nucleotide substitution was estimated to be 0.2 +/-, 0.1% base pair change per million years for both single-copy and middle-repetitive DNAs. The size of the D. arizonensis genome, a representative of the mulleri complex, was calculated to be 2.2 X 10(8) base pairs from its kinetic complexity similar to that of D. hydei. The relative amounts (18%) and average reiteration frequency (100 copies) of the middle-repetitive DNA are similar for all Drosophila species studied. Finally, the data are presented in a phylogenetic tree.     

Multiple events of horizontal transfer of the Minos transposable element between Drosophila species

In this study the Minos element was analyzed in 26 species of the repleta group and seven species of the saltans group of the genus Drosophila. The PCR and Southern blot analysis showed a wide occurrence of the Minos transposable element among species of the repleta and the saltans groups and also a low number of insertions in both genomes. Three different analyses, nucleotide divergence, historical associations, and comparisons between substitution rates (d(N) and d(S)) of Minos and Adh host gene sequences, suggest the occurrence of horizontal transfer between repleta and saltans species. These data reinforce and extend the Arca and Savakis [Genetica 108 (2000) 263] results and suggest five events of horizontal transfer to explain the present Minos distribution: between D. saltans and the ancestor of the mulleri and the mojavensis clusters; between D. hydei and the ancestor of the mulleri and the mojavensis clusters; between D. mojavensis and D. aldrichi; between D. buzzatii and D. serido; and between D. spenceri and D. emarginata. An alternative explanation would be that repeated events of horizontal transfer involving D. hydei, which is a cosmopolitan species that diverged from the others repleta species as long as 14Mya, could have spread Minos within the repleta group and to D. saltans. The data presented in this article support a model in which distribution of Minos transposon among Drosophila species is determined by horizontal transmission balanced by vertical inactivation and extinction.     

Drosophila auraria复合种Period基因Thr-Gly区段的初步研究

测定了金色果蝇复合种（Drosophila auraria species complex）5个姐妹种（sibling species）和D.rufa的period(per)基因的Thr-Gly区段序列。该区段序列分析表明：DNA序列的碱其组成拥有果蝇其他基因的共同特点;颠换数多于转换数,两两种种间的颠换率与转换率的总比值为2～5,密码子第3位的颠换与转换的比值为2.5～5;同义替换/异义替换（Ks/Ka）值远大于10,且有的物种间根本不存在非同义突变,低的Ka值说明该复合种的per基因Thr-Gly区段在进化过程中可能承受着较强的选择压力。运用所得的核苷酸序列构建Drosophila auraria复合种的系统发育树,为澄清该类群的系统演化关系提供了新的线索。     

Phylogenetic relationships among holarctic populations of Chironomus entis and Chironomus plumosus in view of possible horizontal transfer of mitochondrial genes

In eight Holarctic populations of two typical chironomid sibling species of the plumosus group, Chrionomus entis and Chironomus plumosus, nucleotides sequences of mitochondrial (cytb) and nuclear (gb2b) gene regions were examined. The phylogenetic trees reflecting the evolutionary histories of the nuclear and mitochondrial markers exhibited significant differences. On the tree based on the nuclear gene sequences the populations clustered according to their species affiliation, whereas on the tree based on the mitochondrial gene sequences the populations were grouped according to their geographic position. This discrepancy is probably explained by mitochondrial gene flow between sympatric species with incomplete reproductive isolation (sibling species). Based on our results together with the earlier data on nuclear and mitochondrial gene sequences of some other species from the phylogenetic group plumosus, a scheme of phylogenetic relationships within this group is proposed. This scheme is in many ways different from the traditional view on the evolutionary relationships among species of the plumosus group.     

The causes of phylogenetic conflict in a classic Drosophila species group

Bifurcating phylogenies are frequently used to describe the evolutionary history of groups of related species. However, simple bifurcating models may poorly represent the evolutionary history of species that have been exchanging genes. Here, we show that the history of three well-known closely related species, Drosophila pseudoobscura, D. persimilis and D. p. bogotana, is not well represented by a bifurcating phylogenetic tree. The phylogenetic relationships among these species vary widely between different genomic regions. Much of this phylogenetic variation can be explained by the potential of different genomic regions to introgress between species, as measured in laboratory studies. We argue that the utility of multiple markers in species-level phylogenetic studies can be greatly enhanced by knowledge of genomic location and, in the case of hybridizing species, by knowledge of the functional or linkage relationships among the markers and regions of the genome that reduce hybrid fitness.     

Evolution of the nad3-rps12 Gene Cluster in Angiosperm Mitochondria: Comparison of Edited and Unedited Sequences

We have analyzed the nad3-rps12 locus for eight angiosperms in order to compare the utility of mitochondrial DNA and edited mRNA sequences in phylogenetic reconstruction. The two coding regions, containing from 25 to 35 editing sites in the various plants, have been concatenated in order to increase the significance of the analysis. Differing from the corresponding chloroplast sequences, unedited mitochondrial DNA sequences seem to evolve under a quasi-neutral substitution process which undifferentiates the nucleotide substitution rates for the three codon positions. By using complete gene sequences (all codon positions) we found that genomic sequences provide a classical angiosperm phylogenetic tree with a clear-cut grouping of monocotyledons and dicotyledons with Magnoliidae at the basal branch of the tree. Conversely, owing to their low nucleotide substitution rates, edited mRNA sequences were found not to be suitable for studying phylogenetic relationships among angiosperms. Received: 24 January 1996 / Accepted: 5 June 1996     

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.     

Molecular phylogeny of the subgenus Drosophila (Diptera, Drosophilidae) with an emphasis on Neotropical species and groups: a nuclear versus mitochondrial gene approach

The genus Drosophila has played an essential role in many biological studies during the last 100 years but much controversy and many incompletely addressed issues still remain to be elucidated regarding the phylogeny of this genus. Because information on the Neotropical species contained in the subgenus Drosophila is particularly incomplete, with this taxonomic group being underrepresented in many studies, we designed a study to answer some evolutionary questions related to these species. We subjected at least 41 Drosophilidae taxa to a phylogenetic analysis using a 516-base pair (bp) fragment of the alpha-methyldopa (Amd) nuclear gene and a 672 bp fragment of the mitochondrial cytochrome oxidase subunit II (COII) gene both individually and in combination. We found that the subgenus Drosophila is paraphyletic and subdivided into two main clusters: the first containing species traditionally placed in the virilis-repleta radiation and the second assembling species of the immigrans-Hirtodrosophila radiation. Inside the first of these clusters we could detect the monophyly of both the flavopilosa (the sister-clade of the annulimana group) and the mesophragmatica (closely related to the repleta group) species groups. Concerning the immigrans-Hirtodrosophila lineage, Zaprionus, Liodrosophila, Samoaia, and Hirtodrosophila were the early offshoots, followed by the immigrans, quinaria, testacea, and funebris species groups. The tripunctata radiation appears to be a derived clade, composed of a paraphyletic tripunctata group, intimately interposed with members of the cardini, guarani, and guaramunu species groups. Overall, the COII gene yielded a poor phylogenetic performance when compared to the Amd gene, the evolutionary hypothesis of which agreed with the total evidence tree. This phenomenon can be explained by the fast saturation of transitional substitutions in COII, due to strong biases in both base composition and substitution patterns, as also by its great among-site rate variation heterogeneity.     

Toward a comprehensive phylogeny for mammalian and avian herpesviruses

With the aim of deriving a definitive phylogenetic tree for as many mammalian and avian herpesvirus species as possible, alignments were made of amino acid sequences from eight conserved and ubiquitously present genes of herpesviruses, with 48 virus species each represented by at least one gene. Phylogenetic trees for both single-gene and concatenated alignments were evaluated thoroughly by maximum-likelihood methods, with each of the three herpesvirus subfamilies (the Alpha-, Beta-, and Gammaherpesvirinae) examined independently. Composite trees were constructed starting with the top-scoring tree based on the broadest set of genes and supplemented by addition of virus species from trees based on narrower gene sets, to give finally a 46-species tree; branching order for three regions within the tree remained unresolved. Sublineages of the Alpha- and Betaherpesvirinae showed extensive cospeciation with host lineages by criteria of congruence in branching patterns and consistency in extent of divergence. The Gammaherpesvirinae presented a more complex picture, with both higher and lower substitution rates in different sublineages. The final tree obtained represents the most detailed view to date of phylogenetic relationships in any family of large-genome viruses.     

Inferring phylogenies from RAD sequence data

Reduced-representation genome sequencing represents a new source of data for systematics, and its potential utility in interspecific phylogeny reconstruction has not yet been explored. One approach that seems especially promising is the use of inexpensive short-read technologies (e.g., Illumina, SOLiD) to sequence restriction-site associated DNA (RAD)--the regions of the genome that flank the recognition sites of restriction enzymes. In this study, we simulated the collection of RAD sequences from sequenced genomes of different taxa (Drosophila, mammals, and yeasts) and developed a proof-of-concept workflow to test whether informative data could be extracted and used to accurately reconstruct "known" phylogenies of species within each group. The workflow consists of three basic steps: first, sequences are clustered by similarity to estimate orthology; second, clusters are filtered by taxonomic coverage; and third, they are aligned and concatenated for "total evidence" phylogenetic analysis. We evaluated the performance of clustering and filtering parameters by comparing the resulting topologies with well-supported reference trees and we were able to identify conditions under which the reference tree was inferred with high support. For Drosophila, whole genome alignments allowed us to directly evaluate which parameters most consistently recovered orthologous sequences. For the parameter ranges explored, we recovered the best results at the low ends of sequence similarity and taxonomic representation of loci; these generated the largest supermatrices with the highest proportion of missing data. Applications of the method to mammals and yeasts were less successful, which we suggest may be due partly to their much deeper evolutionary divergence times compared to Drosophila (crown ages of approximately 100 and 300 versus 60 Mya, respectively). RAD sequences thus appear to hold promise for reconstructing phylogenetic relationships in younger clades in which sufficient numbers of orthologous restriction sites are retained across species.     

Interspecific divergence, intrachromosomal recombination, and phylogenetic utility of Y-chromosomal genes in Drosophila

Reconstruction of phylogenetic relationships among recently diverged species is complicated by three general problems: segregation of polymorphisms that pre-date species divergence, gene flow during and after speciation, and intra-locus recombination. In light of these difficulties, the Y chromosome offers several important advantages over other genomic regions as a source of phylogenetic information. These advantages include the absence of recombination, rapid coalescence, and reduced opportunity for interspecific introgression due to hybrid male sterility. In this report, we test the phylogenetic utility of Y-chromosomal sequences in two groups of closely related and partially inter-fertile Drosophila species. In the D. bipectinata species complex, Y-chromosomal loci unambiguously recover the phylogeny most consistent with previous multi-locus analysis and with reproductive relationships, and show no evidence of either post-speciation gene flow or persisting ancestral polymorphisms. In the D. simulans species complex, the situation is complicated by the duplication of at least one Y-linked gene region, followed by intrachromosomal recombination between the duplicate genes that scrambles their genealogy. We suggest that Y-chromosomal sequences are a useful tool for resolving phylogenetic relationships among recently diverged species, especially in male-heterogametic organisms that conform to Haldane's rule. However, duplication of Y-linked genes may not be uncommon, and special care should be taken to distinguish between orthologous and paralogous sequences.     

Phylogenetic Relationships among Holarctic Populations of Chironomus entis and Chironomus plumosus in View of Possible Horisontal Transfer of Mitochondrial Genes

In eight Holarctic populations of two typical chironomid sibling species of the plumosus group, Chironomus entisandChironomus plumosus, nucleotide sequences of mitochondrial (cytb) and nuclear (gb2b) gene regions were examined. The phylogenetic trees reflecting the evolutionary histories of the nuclear and mitochondrial markers exhibited significant differences. On the tree based on the nuclear gene sequences the populations clustered according to their species affiliation, whereas on the tree based on the mitochondrial gene sequences the populations were grouped according to their geographic position. This discrepancy is probably explained by mitochondrial gene flow between sympatric species with incomplete reproductive isolation (sibling species). Based on our results together with the earlier data on nuclear and mitochondrial gene sequences of some other species from the phylogenetic group plumosus, a scheme of phylogenetic relationships within this group is proposed. This scheme is in many ways different from the traditional view on the evolutionary relationships among species of the plumosus group.     

Origins and divergence times of mammalian class II MHC gene clusters

The class I and II major histocompatibility complex (MHC) genes are apparently subject to evolution by a birth-and-death process. The rate of gene turnover is much slower in the latter genes than in the former. In placental mammals, the class II region can be subdivided into different orthologous subregions or gene clusters (DR, DQ, DO, and DN), but the origins and evolutionary relationships of these gene clusters are not well established. Here we report the results of our study of the times of origin and evolutionary relationships of these gene clusters in mammals. Our analysis suggests that both class II alpha-chain and beta-chain gene clusters are shared by placental mammals and marsupials, but the gene clusters from nonmammalian species are paralogous to mammalian gene clusters. We estimated the times of divergence between gene clusters in placental mammals using the linearized tree and distance regression methods. Our results indicate that most gene clusters originated 170-200 million years (MY) ago, but that DO beta-chain genes diverged from the other beta-chain gene clusters approximately 210-260 MY ago. The phylogenetic trees for the alpha- and beta-chain genes were not congruent, suggesting that the evolutionary history of the class II gene clusters is more complex than previously thought.     

Assessing the performance of single-copy genes for recovering robust phylogenies

Phylogenies involving nonmodel species are based on a few genes, mostly chosen following historical or practical criteria. Because gene trees are sometimes incongruent with species trees, the resulting phylogenies may not accurately reflect the evolutionary relationships among species. The increase in availability of genome sequences now provides large numbers of genes that could be used for building phylogenies. However, for practical reasons only a few genes can be sequenced for a wide range of species. Here we asked whether we can identify a few genes, among the single-copy genes common to most fungal genomes, that are sufficient for recovering accurate and well-supported phylogenies. Fungi represent a model group for phylogenomics because many complete fungal genomes are available. An automated procedure was developed to extract single-copy orthologous genes from complete fungal genomes using a Markov Clustering Algorithm (Tribe-MCL). Using 21 complete, publicly available fungal genomes with reliable protein predictions, 246 single-copy orthologous gene clusters were identified. We inferred the maximum likelihood trees using the individual orthologous sequences and constructed a reference tree from concatenated protein alignments. The topologies of the individual gene trees were compared to that of the reference tree using three different methods. The performance of individual genes in recovering the reference tree was highly variable. Gene size and the number of variable sites were highly correlated and significantly affected the performance of the genes, but the average substitution rate did not. Two genes recovered exactly the same topology as the reference tree, and when concatenated provided high bootstrap values. The genes typically used for fungal phylogenies did not perform well, which suggests that current fungal phylogenies based on these genes may not accurately reflect the evolutionary relationships among species. Analyses on subsets of species showed that the phylogenetic performance did not seem to depend strongly on the sample. We expect that the best-performing genes identified here will be very useful for phylogenetic studies of fungi, at least at a large taxonomic scale. Furthermore, we compare the method developed here for finding genes for building robust phylogenies with previous ones and we advocate that our method could be applied to other groups of organisms when more complete genomes are available.