Failure of the ILD to determine data combinability for slow loris phylogeny

Tests for incongruence as an indicator of among-data partition conflict have played an important role in conditional data combination. When such tests reveal significant incongruence, this has been interpreted as a rationale for not combining data into a single phylogenetic analysis. In this study of lorisiform phylogeny, we use the incongruence length difference (ILD) test to assess conflict among three independent data sets. A large morphological data set and two unlinked molecular data sets--the mitochondrial cytochrome b gene and the nuclear interphotoreceptor retinoid binding protein (exon 1)--are analyzed with various optimality criteria and weighting mechanisms to determine the phylogenetic relationships among slow lorises (Primates, Loridae). When analyzed separately, the morphological data show impressive statistical support for a monophyletic Loridae. Both molecular data sets resolve the Loridae as paraphyletic, though with different branching orders depending on the optimality criterion or character weighting used. When the three data partitions are analyzed in various combinations, an inverse relationship between congruence and phylogenetic accuracy is observed. Nearly all combined analyses that recover monophyly indicate strong data partition incongruence (P = 0.00005 in the most extreme case), whereas all analyses that recover paraphyly indicate lack of significant incongruence. Numerous lines of evidence verify that monophyly is the accurate phylogenetic result. Therefore, this study contributes to a growing body of information affirming that measures of incongruence should not be used as indicators of data set combinability.     

Phylogeny of frogs of the Physalaemus pustulosus species group,with an examination of data incongruence

Characters derived from advertisement calls, morphology, allozymes, and the sequences of the small subunit of the mitochondrial ribosomal gene (12S) and the cytochrome oxidase I (COI) mitochondrial gene were used to estimate the phylogeny of frogs of the Physalaemus pustulosus group (Leptodactylidae). The combinability of these data partitions was assessed in several ways: measures of phylogenetic signal, character support for trees, congruence of tree topologies, compatibility of data partitions with suboptimal trees, and homogeneity of data partitions. Combined parsimony analysis of all data equally weighted yielded the same tree as the 12S partition analyzed under parsimony and maximum likelihood. The COI, allozyme, and morphology partitions were generally congruent and compatible with the tree derived from combined data. The call data were significantly different from all other partitions, whether considered in terms of tree topology alone, partition homogeneity, or compatibility of data with trees derived from other partitions. The lack of effect of the call data on the topology of the combined tree is probably due to the small number of call characters. The general incongruence of the call data with other data partitions is consistent with the idea that the advertisement calls of this group of frogs are under strong sexual selection.     

Can three incongruence tests predict when data should be combined?

Advocates of conditional combination have argued that testing for incongruence between data partitions is an important step in data exploration. Unless the partitions have had distinct histories, as in horizontal gene transfer, incongruence means that one or more data support the wrong phylogeny. This study examines the relationship between incongruence and phylogenetic accuracy using three tests of incongruence. These tests were applied to pairs of mitochondrial DNA data partitions from two well-corroborated vertebrate phylogenies. Of the three tests, the most useful was the incongruence length difference test (ILD, also called the partition homogeneity test). This test distinguished between cases in which combining the data generally improved phylogenetic accuracy (P > 0.01) and cases in which accuracy of the combined data suffered relative to the individual partitions (P < 0.001). In contrast, in several cases, the Templeton and Rodrigo tests detected highly significant incongruence (P < 0.001) even though combining the incongruent partitions actually increased phylogenetic accuracy. All three tests identified cases in which improving the reconstruction model would improve the phylogenetic accuracy of the individual partitions.      

Phylogeny of snakes (Serpentes): Combining morphological and molecular data in likelihood,Bayesian and parsimony analyses

Abstract

The phylogeny of living and fossil snakes is assessed using likelihood and parsimony approaches and a dataset combining 263 morphological characters with mitochondrial (2693 bp) and nuclear (1092 bp) gene sequences. The ‘no common mechanism’ (NCMr) and ‘Markovian’ (Mkv) models were employed for the morphological partition in likelihood analyses; likelihood scores in the NCMr model were more closely correlated with parsimony tree lengths. Both models accorded relatively less weight to the molecular data than did parsimony, with the effect being milder in the NCMr model. Partitioned branch and likelihood support values indicate that the mtDNA and nuclear gene partitions agree more closely with each other than with morphology. Despite differences between data partitions in phylogenetic signal, analytic models, and relative weighting, the parsimony and likelihood analyses all retrieved the following widely accepted groups: scolecophidians, alethinophidians, cylindrophiines, macrostomatans (sensu lato) and caenophidians. Anilius alone emerged as the most basal alethinophidian; the combined analyses resulted in a novel and stable position of uropeltines and cylindrophiines as the second‐most basal clade of alethinophidians. The limbed marine pachyophiids, along with Dinilysia and Wonambi, were always basal to all living snakes. Other results stable in all combined analyses include: Xenopeltis and Loxocemus were sister taxa (fide morphology) but clustered with pythonines (fide molecules), and Ungaliophis clustered with a boine‐erycine clade (fide molecules). Tropidophis remains enigmatic; it emerges as a basal alethinophidian in the parsimony analyses (fide molecules) but a derived form in the likelihood analyses (fide morphology), largely due to the different relative weighting accorded to data partitions.     

Attributes and congruence of three molecular data sets: inferring phylogenies among Septoria-related species from woody perennial plants

To improve our understanding of phylogenetic relationships within the anamorphic genus Septoria, three molecular data sets representing 2,417 bp of nuclear and mitochondrial genes were evaluated. Separate gene analyses and combined analyses were performed using first, the maximum parsimony criterion and second, a Bayesian framework. The homogeneity of data partitions was evaluated via a combination of homogeneity partition tests and tree topology incongruence tests before conducting combined analyses. A last incongruence re-evaluation using partitioned Bremer support was performed on the combined tree, which corroborated the previous estimates. After each separate data set attributes were examined, simple explanations were advocated as the causes of the significant incongruences detected. The analysis of multiple gene partitions showed unprecedented phylogenetic resolution within the genus Septoria that supported the results from previously published single gene phylogenies. Specifically, we have delimited distinct but closely related species representing monophyletic groups that frequently correlated with their respective host families. Conversely, the occurrence of well-supported groups including closely related but distinct molecular taxa sampled on unrelated host-plants allowed us to reject, in these particular cases, the co-evolutionary concept expected between a parasite and its host and to discuss alternative evolutionary models recently proposed for these pathogens.     

Phylogenetic relationships of horned lizards (Phrynosoma) based on nuclear and mitochondrial data: evidence for a misleading mitochondrial gene tree

It has proven remarkably difficult to obtain a well-resolved and strongly supported phylogeny for horned lizards (Phrynosoma) because of incongruence between morphological and mitochondrial DNA sequence data. We infer the phylogenetic relationships among all 17 extant Phrynosoma species using >5.1 kb of mtDNA (12S rRNA, 16S rRNA, ND1, ND2, ND4, Cyt b, and associated tRNA genes), and >2.2kb from three nuclear genes (RAG-1, BDNF, and GAPD) for most taxa. We conduct separate and combined phylogenetic analyses of these data using maximum parsimony, maximum likelihood, and Bayesian methods. The phylogenetic relationships inferred from the mtDNA data are congruent with previous mtDNA analyses based on fewer characters and provide strong support for most branches. However, we detected strong incongruence between the mtDNA and nuclear data using comparisons of branch support and Shimodaira-Hasegawa tests, with the (P. platyrhinos+P. goodei) clade identified as the primary source of this conflict. Our analysis of a P. mcalliixP. goodei hybrid suggests that this incongruence is caused by reticulation via introgressive hybridization. Our preferred phylogeny based on an analysis of the combined data (excluding the introgressed mtDNA data) provides a new framework for interpreting character evolution and biogeography within Phrynosoma. In the context of this improved phylogeny we propose a phylogenetic taxonomy highlighting four clades: (1) Tapaja, containing the viviparous short-horned lizards P. ditmarsi, P. hernandesi, P. douglasii, and P. orbiculare; (2) Anota, containing species with prominent cranial horns (P. solare, P. mcallii, and the P. coronatum group); (3) Doliosaurus, containing three species lacking antipredator blood-squirting (P. modestum, P. platyrhinos, and P. goodei); and (4) Brevicauda, containing two viviparous species with extremely short tails that lack blood-squirting (P. braconnieri and P. taurus).     

New clades of euthyneuran gastropods (Mollusca) from 28S rRNA sequences

Recent morphological and molecular results on phylogeny of euthyneuran gastropods, which include opisthobranchs and pulmonates, have greatly diminished previous supposed resolution of their phylogenetic relationships. In addition to recent morphological results, sequences of the D1 and D2 domains of the 28S rRNA are here analyzed by parsimony for 31 euthyneuran species. The molecular and previous morphological data sets were not congruent according to an ILD test, and morphological and molecular data could not be analyzed simultaneously. Consequently Bremer's Combinable Component Consensus was used to obtain a new tree, with the following supported molecular results: monophyly of a new clade of opisthobranchs including actively swimming Euthyneura, i.e., pelagic Gymnosomata and Thecosomata plus benthic Anaspidea; first molecular confirmation of monophylies of Hygrophila, including Chilina, Acteonoidea, and Sacoglossa, which include both shell-bearing species and slugs; and new confirmation of the monophyly of Stylommatophora. Morphological characters which support the new clades obtained here are discussed.     

Assessing phylogenetic resolution among mitochondrial, nuclear, and morphological datasets in Nothonotus darters (Teleostei: Percidae)

External morphological characters are the basis of our understanding of diversity and species relationships in many darter clades. The past decade has seen the publication of many studies utilizing mtDNA sequence data to investigate darter phylogenetics, but only recently have nuclear genes been used to investigate darter relationships. Despite a long tradition of use in darter systematics few studies have examined the phylogenetic utility of external morphological characters in estimating relationships among species in darter clades. We present DNA sequence data from the mitochondrial cytochrome b (cytb) gene, the nuclear encoded S7 intron 1, and discretely coded external morphological characters for all 20 species in the darter clade Nothonotus. Bayesian phylogenetic analyses result in phylogenies that are in broad agreement with previous studies. The cytb gene tree is well resolved, while the nuclear S7 gene tree lacks phylogenetic resolution, node support, and is characterized by a lack of reciprocal monophyly for many of the Nothonotus species. The phylogenies resulting from analysis of the morphological dataset lack resolution, but nodes present are found in the cytb and S7 gene trees. The highest resolution and node support is found in the Bayesian combined data phylogeny. Based on our results we propose continued exploration of the phylogenetic utility of external morphological characters in other darter clades. Given the extensive lack of reciprocal monophyly of species observed in the S7 gene tree we predict that nuclear gene sequences may have limited utility in intraspecific phylogeographic studies of Nothonotus darters.     

Genome-scale phylogenetics: inferring the plant tree of life from 18,896 gene trees

Phylogenetic analyses using genome-scale data sets must confront incongruence among gene trees, which in plants is exacerbated by frequent gene duplications and losses. Gene tree parsimony (GTP) is a phylogenetic optimization criterion in which a species tree that minimizes the number of gene duplications induced among a set of gene trees is selected. The run time performance of previous implementations has limited its use on large-scale data sets. We used new software that incorporates recent algorithmic advances to examine the performance of GTP on a plant data set consisting of 18,896 gene trees containing 510,922 protein sequences from 136 plant taxa (giving a combined alignment length of >2.9 million characters). The relationships inferred from the GTP analysis were largely consistent with previous large-scale studies of backbone plant phylogeny and resolved some controversial nodes. The placement of taxa that were present in few gene trees generally varied the most among GTP bootstrap replicates. Excluding these taxa either before or after the GTP analysis revealed high levels of phylogenetic support across plants. The analyses supported magnoliids sister to a eudicot + monocot clade and did not support the eurosid I and II clades. This study presents a nuclear genomic perspective on the broad-scale phylogenic relationships among plants, and it demonstrates that nuclear genes with a history of duplication and loss can be phylogenetically informative for resolving the plant tree of life.     

Multiple sources of character information and the phylogeny of Hawaiian drosophilids

Relationships among representatives of the five major Hawaiian Drosophila species groups were examined using data from eight different gene regions. A simultaneous analysis of these data resulted in a single most-parsimonious tree that (1) places the adiastola picture-winged subgroup as sister taxon to the other picture-winged subgroups, (2) unites the modified-tarsus species group with flies from the Antopocerus species group, and (3) places the white-tip scutellum species group as the most basal taxon. Because of the different gene sources used in this study, numerous process partitions can be erected within this data set. We examined the incongruence among these various partitions and the ramifications of these data for the taxonomic consensus, prior agreement, and simultaneous analysis approaches to phylogenetic reconstruction. Separate analyses and taxonomic consensus appear to be inadequate methods for dealing with the partitions in this study. Although detection of incongruence is possible and helps elucidate particular areas of disagreement among data sets, separation of partitions on the basis of incongruence is problematic for many reasons. First, analyzing all genes separately and then either presenting them all as possible hypotheses or taking their consensus provides virtually no information concerning the relationships among these flies. Second, despite some evidence of incongruence, there are no clear delineations among the various gene partitions that separate only heterogeneous data. Third, to the extent that problematic genes can be identified, these genes have nearly the same information content, within a combined analysis framework, as the remaining nonproblematic genes. Our data suggest that significant incongruence among data partitions may be isolated to specific relationships and the "false" signal creating this incongruence is most likely to be overcome by a simultaneous analysis. We present a new method, partitioned Bremer support, for examining the contribution of a particular data partition to the topological support of the simultaneous analysis tree.     

Phylogeny of the cricket subfamily Eneopterinae (Orthoptera, Grylloidea, Eneopteridae) based on four molecular loci and morphology

The phylogenetic relationships of 39 species of Eneopterinae crickets are reconstructed using four molecular markers (16S rRNA, 12S rRNA, cytochrome b, 18S rRNA) and a large morphological data set. Phylogenetic analysis via direct optimisation of DNA sequence data using parsimony as optimality criterion is done for six combinations of weighting parameter sets in a sensitivity analysis. The results are discussed in a twofold purpose: first, in term of significance of the molecular markers for phylogeny reconstruction in Ensifera, as our study represents the first molecular phylogeny performed for this insect suborder at this level of diversity; second, in term of corroboration of a previous phylogeny of Eneopterinae, built on morphological data alone. The four molecular markers all convey phylogenetic signal, although variously distributed on the tree. The monophyly of the subfamily, that of three over five tribes, and of 10 over 13 genera, are recovered. Finally, previous hypotheses on the evolution of acoustic devices and signals in the Eneopterinae clade are briefly tested, and supported, by our new data set.     

Phylogeny of the Drosophila saltans species group based on combined analysis of nuclear and mitochondrial DNA sequences

Nucleotide sequences from two nuclear loci, alcohol dehydrogenase and internal transcribed spacer-1 of the nuclear ribosomal DNA repeats, and two mitochondrial genes, cytochrome oxidase I and cytochrome oxidase II, were determined from nine species in the Drosophila saltans species group. The partition homogeneity test and partitioned Bremer support were used to measure incongruence between phylogenetic hypotheses generated from individual partitions. Individual loci were generally congruent with each other and consistent with the previously proposed morphological hypothesis, although they differed in level of resolution. Since extreme conflict between partitions did not exist, the data were combined and analyzed simultaneously. The total evidence method gave a more resolved and highly supported phylogeny, as indicated by bootstrap proportions and decay indices, than did any of the individual analyses. The cordata and elliptica subgroups, considered to have diverged early in the history of the D. saltans group, were sister taxa to the remainder of the saltans group. The sturtevanti subgroup, represented by D. milleri and D. sturtevanti, occupies an intermediate position in this phylogeny. The saltans and parasaltans subgroups are sister clades and occupy the most recently derived portion of the phylogeny. As with previous morphological studies, phylogenetic relationships within the saltans subgroup were not satisfactorily resolved by the molecular data.      

Phylogenetic relationships in the subgenus Mus (genus Mus, family Muridae, subfamily Murinae): examining gene trees and species trees

Over eight kilobases (kb) of sequence from eight genes including two mitochondrial loci, Cyt b and 12S, and six nuclear loci, B2m , Zp3 , Tcp1, Sry, Smcx and Smcy , were used to investigate phylogenetic relationships among 11 taxa representing eight species within the rodent genus Mus . Particular attention was given to discerning relationships among species within the subgenus Mus including members of a Palearctic clade ( M. musculus , M. spicilegus , M. macedonicus and M. spretus ) and members of an Asian clade ( M. caroli , M. cookii and M. cervicolor ), as previous studies using different datasets have produced different topologies for taxa within these two groups. While parsimony analyses of the combined eight-gene dataset yielded a single, fully resolved tree, support values were lower for nodes resolving relationships within the Palearctic and Asian clades than they were elsewhere in the tree. In addition, a maximum likelihood analysis of the same eight-gene dataset yielded different topologies for both the Palearctic and the Asian clades. Both observations are indicative of clade instability. The nature of this instability was explored through a comparison with our previous study in which we included the two mitochondrial loci and only four of the six nuclear genes, and through an analysis of partitioned data, specifically mitochondrial vs. nuclear genes. This study underscores the importance of considering among-site rate variation in phylogeny reconstruction. © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 653–662.     

Evolutionary relationships of the cup-fungus genus Peziza and Pezizaceae inferred from multiple nuclear genes: RPB2, beta-tubulin, and LSU rDNA

To provide a robust phylogeny of Pezizaceae, partial sequences from two nuclear protein-coding genes, RPB2 (encoding the second largest subunit of RNA polymerase II) and beta-tubulin, were obtained from 69 and 72 specimens, respectively, to analyze with nuclear ribosomal large subunit RNA gene sequences (LSU). The three-gene data set includes 32 species of Peziza, and 27 species from nine additional epigeous and six hypogeous (truffle) pezizaceous genera. Analyses of the combined LSU, RPB2, and beta-tubulin data set using parsimony, maximum likelihood, and Bayesian approaches identify 14 fine-scale lineages within Pezizaceae. Species of Peziza occur in eight of the lineages, spread among other genera of the family, confirming the non-monophyly of the genus. Although parsimony analyses of the three-gene data set produced a nearly completely resolved strict consensus tree, with increased confidence, relationships between the lineages are still resolved with mostly weak bootstrap support. Bayesian analyses of the three-gene data, however, show support for several more inclusive clades, mostly congruent with Bayesian analyses of RPB2. No strongly supported incongruence was found among phylogenies derived from the separate LSU, RPB2, and beta-tubulin data sets. The RPB2 region appeared to be the most informative single gene region based on resolution and clade support, and accounts for the greatest number of potentially parsimony informative characters within the combined data set, followed by the LSU and the beta-tubulin region. The results indicate that third codon positions in beta-tubulin are saturated, especially for sites that provide information about the deeper relationships. Nevertheless, almost all phylogenetic signal in beta-tubulin is due to third positions changes, with almost no signal in first and second codons, and contribute phylogenetic information at the "fine-scale" level within the Pezizaceae. The Pezizaceae is supported as monophyletic in analyses of the three-gene data set, but its sister-group relationships is not resolved with support. The results advocate the use of RPB2 as a marker for ascomycete phylogenetics at the inter-generic level, whereas the beta-tubulin gene appears less useful.     

The phylogeny of the social Anelosimus spiders (Araneae: Theridiidae) inferred from six molecular loci and morphology

We use fragments of three nuclear genes (Histone 3, 18SrDNA, and 28SrDNA) and three mitochondrial genes (16SrDNA, ND1, and COI) totalling approximately 4.5kb, in addition to morphological data, to estimate the phylogenetic relationships among Anelosimus spiders, well known for their sociality. The analysis includes 67 individuals representing 23 of the 53 currently recognized Anelosimus species and all species groups previously recognized by morphological evidence. We analyse the data using Bayesian, maximum likelihood, and parsimony methods, considering the genes individually as well as combined (mitochondrial, nuclear, and both combined) in addition to a 'total evidence' analysis including morphology. Most of the data partitions are congruent in agreeing on several fundamental aspects of the phylogeny, and the combined molecular data yield a tree broadly similar to an existing morphological hypothesis. We argue that such congruence among data partitions is an important indicator of support that may go undetected by standard robustness estimators. Our results strongly support Anelosimus monophyly, and the monophyly of the recently revised American 'eximius lineage', although slightly altered by excluding A. pacificus. There was consistent support for the scattering of American Anelosimus species in three clades suggesting intercontinental dispersal. Several recently described species are reconstructed as monophyletic, supporting taxonomic decisions based on morphology and behaviour in this taxonomically difficult group. Corroborating previous results from morphology, the molecular data suggest that social species are scattered across the genus and thus that sociality has evolved multiple times, a significant finding for exploring the causes and consequences of social evolution in this group of organisms.     

Data congruence, paedomorphosis and salamanders

Background

The retention of ancestral juvenile characters by adult stages of descendants is called paedomorphosis. However, this process can mislead phylogenetic analyses based on morphological data, even in combination with molecular data, because the assessment if a character is primary absent or secondary lost is difficult. Thus, the detection of incongruence between morphological and molecular data is necessary to investigate the reliability of simultaneous analyses. Different methods have been proposed to detect data congruence or incongruence. Five of them (PABA, PBS, NDI, LILD, DRI) are used herein to assess incongruence between morphological and molecular data in a case study addressing salamander phylogeny, which comprises several supposedly paedomorphic taxa. Therefore, previously published data sets were compiled herein. Furthermore, two strategies ameliorating effects of paedomorphosis on phylogenetic studies were tested herein using a statistical rigor. Additionally, efficiency of the different methods to assess incongruence was analyzed using this empirical data set. Finally, a test statistic is presented for all these methods except DRI.

Results

The addition of morphological data to molecular data results in both different positions of three of the four paedomorphic taxa and strong incongruence, but treating the morphological data using different strategies ameliorating the negative impact of paedomorphosis revokes these changes and minimizes the conflict. Of these strategies the strategy to just exclude paedomorphic character traits seem to be most beneficial. Of the three molecular partitions analyzed herein the RAG1 partition seems to be the most suitable to resolve deep salamander phylogeny. The rRNA and mtDNA partition are either too conserved or too variable, respectively. Of the different methods to detect incongruence, the NDI and PABA approaches are more conservative in the indication of incongruence than LILD and PBS.

Conclusion

Paedomorphosis induces strong conflicts and can mislead the phylogenetic analyses even in combined analyses. However, different strategies are efficiently minimizing these problems. Though the exploration of different methods to detect incongruence is preferable NDI and PABA are more conservative than the others and NDI is computational less extensive than PABA.     

Phylogeny of bumble bees in the New World subgenus Fervidobombus (Hymenoptera: Apidae): congruence of molecular and morphological data

We present new DNA sequence data (12S, 16S, and opsin gene fragments) and morphological characters of the male genitalia for a phylogenetic analysis of the bumble bee subgenus Fervidobombus. There is no significant incongruence between the three molecular data sets, and little incongruence between the DNA and morphology. Simultaneous analysis of all the data partitions resulted in a tree that was entirely congruent with the All-DNA tree. Optimization of the geographic locations of the taxa onto the tree topology using dispersal/vicariance analysis suggests a complex picture of spread and diversification of Fervidobombus from the Old World into the southern New World. There is a phylogenetic component to their spread into tropical rain forest, as the two primary rain forest species (Bombus transversalis and Bombus pullatus) comprise a monophyletic clade, along with a third species, Bombus atratus, which is widely distributed in South America, including lowland subtropical habitats.     

Influence of data conflict and molecular phylogeny of major clades in Cimicomorphan true bugs (Insecta: Hemiptera: Heteroptera)

Cimicomorpha, which consists of 16 families representing more than 19,400 species, is the largest infraorder in Heteroptera, Insecta. We present the first molecular phylogenetic investigation of family relationships of Cimicomorpha, including 46 taxa from 12 of 16 Cimicomorphan families. Three genes, with a total of 3277 bp of sequence data (nuclear 18S rDNA: 2022 bp, 28S rDNA: 755 bp, and mitochondrial 16S rDNA: 498 bp) were analyzed. Data partitions were analyzed separately and in combination, by employing ML (maximum likelihood), MP (maximum parsimony), and Bayesian methods. As saturation was detected in substitutions of 16S rDNA, influence of data conflict in combined analyses was further explored by three methods: the incongruence length difference (ILD) test, the partitioned Bremer support (PBS), and the partition addition bootstrap alteration approach (PABA). PBS and PABA approaches suggested that 16S rDNA was not very suitable for addressing relationships at this level in Cimicomorpha. Our results also supported the nabid-cimicoid lineage for Cimicoidea proposed by Schuh and Stys [Schuh, R.T., Stys, P., 1991. Phylogenetic analysis of Cimicomorphan family relationships (Heteroptera). J. NY Entomol. Soc. 99 (3), 298-350]. Data incongruence and the utility of the three genes were briefly discussed.     

Phylogeny of the Oriental Drosophila melanogaster species group: a multilocus reconstruction

The melanogaster species group of Drosophila (subgenus Sophophora) has long been a favored model for evolutionary studies because of its morphological and ecological diversity and wide geographic distribution. However, phylogenetic relationships among species and subgroups within this lineage are not well understood. We reconstructed the phylogeny of 17 species representing 7 "oriental" species subgroups, which are especially closely related to D. melanogaster. We used DNA sequences of four nuclear and two mitochondrial loci in an attempt to obtain the best possible estimate of species phylogeny and to assess the extent and sources of remaining uncertainties. Comparison of trees derived from single-gene data sets allowed us to identify several strongly supported clades, which were also consistently seen in combined analyses. The relationships among these clades are less certain. The combined data set contains data partitions that are incongruent with each other. Trees reconstructed from the combined data set and from internally homogenous data sets consisting of three or four genes each differ at several deep nodes. The total data set tree is fully resolved and strongly supported at most nodes. Statistical tests indicated that this tree is compatible with all individual and combined data sets. Therefore, we accepted this tree as the most likely model of historical relationships. We compared the new molecular phylogeny to earlier estimates based on morphology and chromosome structure and discuss its taxonomic and evolutionary implications.