Additional support for Afrotheria and Paenungulata,the performance of mitochondrial versus nuclear genes,and the impact of data partitions with heterogeneous base composition

We concatenated sequences for four mitochondrial genes (12S rRNA, tRNA valine, 16S rRNA, cytochrome b) and four nuclear genes [aquaporin, alpha 2B adrenergic receptor (A2AB), interphotoreceptor retinoid-binding protein (IRBP), von Willebrand factor (vWF)] into a multigene data set representing 11 eutherian orders (Artiodactyla, Hyracoidea, Insectivora, Lagomorpha, Macroscelidea, Perissodactyla, Primates, Proboscidea, Rodentia, Sirenia, Tubulidentata). Within this data set, we recognized nine mitochondrial partitions (both stems and loops, for each of 12S rRNA, tRNA valine, and 16S rRNA; and first, second, and third codon positions of cytochrome b) and 12 nuclear partitions (first, second, and third codon positions, respectively, of each of the four nuclear genes). Four of the 21 partitions (third positions of cytochrome b, A2AB, IRBP, and vWF) showed significant heterogeneity in base composition across taxa. Phylogenetic analyses (parsimony, minimum evolution, maximum likelihood) based on sequences for all 21 partitions provide 99-100% bootstrap support for Afrotheria and Paenungulata. With the elimination of the four partitions exhibiting heterogeneity in base composition, there is also high bootstrap support (89-100%) for cow + horse. Statistical tests reject Altungulata, Anagalida, and Ungulata. Data set heterogeneity between mitochondrial and nuclear genes is most evident when all partitions are included in the phylogenetic analyses. Mitochondrial-gene trees associate cow with horse, whereas nuclear-gene trees associate cow with hedgehog and these two with horse. However, after eliminating third positions of A2AB, IRBP, and vWF, nuclear data agree with mitochondrial data in supporting cow + horse. Nuclear genes provide stronger support for both Afrotheria and Paenungulata. Removal of third positions of cytochrome b results in improved performance for the mitochondrial genes in recovering these clades.     

Highly Congruent Molecular Support for a Diverse Superordinal Clade of Endemic African Mammals

A solution to higher level mammalian phylogeny is going to depend on the congruent establishment of superordinal groupings followed by a linking together of these clades. We present congruent and convincing evidence from four disparate nuclear protein coding genes and from a tandem alignment of the 12S–16S mitochondrial region, for a superordinal clade of endemic African mammals that includes elephant shrews, aardvarks, golden mole, elephants, sirenians, and hyraxes. Because of strong support for golden mole as part of this clade, the Insectivora are rendered paraphyletic or polyphyletic, with constrained monophyly of the insectivores judged significantly worse in the vast majority of tests. Branching arrangement within this clade remains highly uncertain; however, a tandem alignment of the protein coding genes suggests elephant shrew is the earliest African lineage. None of the individual data sets or combinations of data sets support the widely held view of a mirorder Tethytheria (Sirenia/Proboscidea), although only a tandem alignment of protein coding and mitochondrial loci significantly rejects this association. The majority of the data sets and analyses provide strong support for Caviomorpha as part of a monophyletic Rodentia.     

A retroposon analysis of Afrotherian phylogeny

Recent comprehensive studies of DNA sequences support the monophyly of Afrotheria, comprising elephants, sirenians (dugongs and manatees), hyraxes, tenrecs, golden moles, aardvarks, and elephant shrews, as well as that of Paenungulata, comprising elephants, sirenians, and hyraxes. However, phylogenetic relationships among paenungulates, as well as among nonpaenungulates, have remained ambiguous. Here we applied an extensive retroposon analysis to these problems to support the monophyly of aardvarks, tenrecs, and golden moles, with elephant shrews as their sister group. Regarding phylogenetic relationships in Paenungulata, we could characterize only one informative locus, although we could isolate many insertions specific to each of three lineages, namely, Proboscidea, Sirenia, and Hyracoidea. These data prompted us to reexamine phylogenetic relationships among Paenungulata using 19 nuclear gene sequences resulting in three different analyses, namely, short interspersed element (SINE) insertions, nuclear sequence analyses, and morphological cladistics, supporting different respective phylogenies. We concluded that these three lineages diverged very rapidly in a very short evolutionary period, with the consequence that ancestral polymorphism present in the last common ancestor of Paenungulata results in such incongruence. Our results suggest the rapid fixation of many large-scale morphological synapomorphies for Tethytheria; implications of this in relation to the morphological evolution in Paenungulata are discussed.     

Paenungulata: a comparison of the hemoglobin sequences from elephant, hyrax, and manatee

Inspection of the amino acid differences among hemoglobin sequences of a wide range of mammalian species suggested that at alpha 19, alpha 110, alpha 111, beta 23, beta 44, and beta 56, synapomorphies group manatee (Trichechus inungius, Sirenia), Indian and African elephant (Elephas maximus and Loxodonta africana, Proboscidea), and rock hyrax (Procavia habessinica, Hyracoidea) into a monophyletic clade. Results obtained by parsimony analysis provide evidence for this grouping--and thus support for the genealogical validity of Simpson's superorder Paenungulata, which contains as the extant orders Proboscidea, Sirenia, and Hyracoidea. All of the 39 most, or nearly most, parsimonious of 10,395 trees constructed from a tandemly combined alpha- and beta- hemoglobin sequence for 103 vertebrate species (of which 79 were mammals from 16 extant orders), depicted Paenungulata as one of the most anciently separated branches of Eutheria. It was found on examining thousands of alternative trees that to not group Proboscidea, Hyracoidea, and Sirenia in a monophyletic clade required at least four additional substitutions.      

Chromosome painting among Proboscidea, Hyracoidea and Sirenia: support for Paenungulata (Afrotheria, Mammalia) but not Tethytheria

Despite marked improvements in the interpretation of systematic relationships within Eutheria, particular nodes, including Paenungulata (Hyracoidea, Sirenia and Proboscidea), remain ambiguous. The combination of a rapid radiation, a deep divergence and an extensive morphological diversification has resulted in a limited phylogenetic signal confounding resolution within this clade both at the morphological and nucleotide levels. Cross-species chromosome painting was used to delineate regions of homology between Loxodonta africana (2n=56), Procavia capensis (2n=54), Trichechus manatus latirostris (2n=48) and an outgroup taxon, the aardvark (Orycteropus afer, 2n=20). Changes specific to each lineage were identified and although the presence of a minimum of 11 synapomorphies confirmed the monophyly of Paenungulata, no change characterizing intrapaenungulate relationships was evident. The reconstruction of an ancestral paenungulate karyotype and the estimation of rates of chromosomal evolution indicate a reduced rate of genomic repatterning following the paenungulate radiation. In comparison to data available for other mammalian taxa, the paenungulate rate of chromosomal evolution is slow to moderate. As a consequence, the absence of a chromosomal character uniting two paenungulates (at the level of resolution characterized in this study) may be due to a reduced rate of chromosomal change relative to the length of time separating successive divergence events.     

The chromosomes of afrotheria and their bearing on Mammalian genome evolution

Afrotheria is the clade of placental mammals that, together with Xenarthra, Euarchontoglires and Laurasiatheria, represents 1 of the 4 main recognized supraordinal eutherian clades. It reunites 6 orders of African origin: Proboscidea, Sirenia, Hyracoidea, Macroscelidea, Afrosoricida and Tubulidentata. The apparently unlikely relationship among such disparate morphological taxa and their possible basal position at the base of the eutherian phylogenetic tree led to a great deal of attention and research on the group. The use of biomolecular data was pivotal in Afrotheria studies, as they were the basis for the recognition of this clade. Although morphological evidence is still scarce, a plethora of molecular data firmly attests to the phylogenetic relationship among these mammals of African origin. Modern cytogenetic techniques also gave a significant contribution to the study of Afrotheria, revealing chromosome signatures for the group as a whole, as well as for some of its internal relationships. The associations of human chromosomes HSA1/19 and 5/21 were found to be chromosome signatures for the group and provided further support for Afrotheria. Additional chromosome synapomorphies were also identified linking elephants and manatees in Tethytheria (the associations HSA2/3, 3/13, 8/22, 18/19 and the lack of HSA4/8) and elephant shrews with the aardvark (HSA2/8, 3/20 and 10/17). Herein, we review the current knowledge on Afrotheria chromosomes and genome evolution. The already available data on the group suggests that further work on this apparently bizarre assemblage of mammals will provide important data to a better understanding on mammalian genome evolution.     

Phylogenetic assessment of molecular and morphological data for eutherian mammals

The interordinal relationships of eutherian (placental) mammals were evaluated by a phylogenetic analysis of four published data sets (three sequences and one morphological). The nature and degree of support and conflict for particular groups were assessed by separate bootstrap and homogeneity tests, which were followed by combined analyses of the sequence and morphological data. Between orders, strong support (i.e., > or = 95% bootstrap scores) was found for a paraphyletic Artiodactyla (relative to Cetacea) and a monophyletic Cetartiodactyla (Artiodactyla and Cetacea) and Paenungulata (Hyracoidea, Proboscidea, and Sirenia). In turn, some reasonable to strong evidence (> or = 85%) was obtained for Hyracoidea with Sirenia, Dermoptera with Scandentia, Glires (Lagomorpha with Rodentia), and Afrotheria (Amblysomus, Macroscelidea, Paenungulata, and Tubulidentata). Otherwise, no other interordinal clades were supported at these reasonable to strong levels. This overall lack of resolution for eutherian interordinal clusters agrees with other studies that suggest further progress will continue to be slow and difficult. Further resolution will require the integration of more recently published data, the continued sampling of taxa and characters, and the use of more powerful methods of data analysis.     

The platypus is in its place: nuclear genes and indels confirm the sister group relation of monotremes and Therians

Morphological data supports monotremes as the sister group of Theria (extant marsupials + eutherians), but phylogenetic analyses of 12 mitochondrial protein-coding genes have strongly supported the grouping of monotremes with marsupials: the Marsupionta hypothesis. Various nuclear genes tend to support Theria, but a comprehensive study of long concatenated sequences and broad taxon sampling is lacking. We therefore determined sequences from six nuclear genes and obtained additional sequences from the databases to create two large and independent nuclear data sets. One (data set I) emphasized taxon sampling and comprised five genes, with a concatenated length of 2,793 bp, from 21 species (two monotremes, six marsupials, nine placentals, and four outgroups). The other (data set II) emphasized gene sampling and comprised eight genes and three proteins, with a concatenated length of 10,773 bp or 3,669 amino acids, from five taxa (a monotreme, a marsupial, a rodent, human, and chicken). Both data sets were analyzed by parsimony, minimum evolution, maximum likelihood, and Bayesian methods using various models and data partitions. Data set I gave bootstrap support values for Theria between 55% and 100%, while support for Marsupionta was at most 12.3%. Taking base compositional bias into account generally increased the support for Theria. Data set II exclusively supported Theria, with the highest possible values and significantly rejected Marsupionta. Independent phylogenetic evidence in support of Theria was obtained from two single amino acid deletions and one insertion, while no supporting insertions and deletions were found for Marsupionta. On the basis of our data sets, the time of divergence between Monotremata and Theria was estimated at 231-217 MYA and between Marsupialia and Eutheria at 193-186 MYA. The morphological evidence for a basal position of Monotremata, well separated from Theria, is thus fully supported by the available molecular data from nuclear genes.     

Base compositional bias and phylogenetic analyses: a test of the "flying DNA" hypothesis

Phylogenetic methods can produce biased estimates of phylogeny when base composition varies along different lineages. Pettigrew (1994, Curr. Biol. 4:277-280) has suggested that base composition bias is responsible for the apparent support for the monophyly of bats (Chiroptera: megabats and microbats) from several different nuclear and mitochondrial genes. Pettigrew's "flying DNA" hypothesis makes several predictions: (1) that metabolic constraints associated with flying result in elevated levels of adenine and thymine throughout the genome of both megabats and microbats, (2) that the resulting base compositional bias in bats is sufficient to mislead phylogenetic methods and account for the support for bat monophyly from several nuclear and mitochondrial genes, and (3) that phylogenetic analysis using pairwise distances corrected for compositional bias should eliminate the support for bat monophyly. We tested these predictions by analyzing DNA sequences from two nuclear and three mitochondrial genes. The predicted base compositional bias does not appear to exist in some of the genes, and in other genes the differences in AT content are very small. Analyses under a wide diversity of criteria and models of evolution, including analyses that take base composition into account (using log-determinant distances), all strongly support bat monophyly. Moreover, simulation analyses indicate that even extreme bias toward AT-base composition in bats would be insufficient to explain the observed levels of support for bat monophyly. These analyses provide no support for the "flying DNA" hypothesis, whereas the monophyly of bats appears to be well supported by the DNA sequence data.     

Molecular Evidence for the Major Clades of Placental Mammals

Higher-level relationships among placental mammals, as well as the historical biogeography of this group against the backdrop of continental fragmentation and reassembly, remain poorly understood. Here, we analyze two independent molecular data sets that represent all placental orders. The first data set includes six genes (A2AB, IRBP, vWF, 12S rRNA, tRNA valine, 16S rRNA; total = 5.71 kb) for 26 placental taxa and two marsupials; the second data set includes 2.95 kb of exon 11 of the BRCA1 gene for 51 placental taxa and four marsupials. We also analyzed a concatenation of these data sets (8.66 kb) for 26 placentals and one marsupial. Unrooted and rooted analyses were performed with parsimony, distance methods, maximum likelihood, and a Bayesian approach. Unrooted analyses provide convincing support for a fundamental separation of placental orders into groups with southern and northern hemispheric origins according to the current fossil record. On rooted trees, one or both of these groups are monophyletic depending on the position of the root. Maximum likelihood and Bayesian analyses with the BRCA1 and combined 8.66 kb data sets provide strong support for the monophyly of the northern hemisphere group (Boreoeutheria). Boreoeutheria is divided into Laurasiatheria (Carnivora + Cetartiodactyla + Chiroptera + Eulipotyphla + Perissodactyla + Pholidota) and Euarchonta (Dermoptera + Primates + Scandentia) + Glires (Lagomorpha + Rodentia). The southern hemisphere group is either monophyletic or paraphyletic, depending on the method of analysis used. Within this group, Afrotheria (Proboscidea + Sirenia + Hyracoidea + Tubulidentata + Macroscelidea + Afrosoricida) is monophyletic. A unique nine base-pair deletion in exon 11 of the BRCA1 gene also supports Afrotheria monophyly. Given molecular dates that suggest that the southern hemisphere group and Boreoeutheria diverged in the Early Cretaceous, a single trans-hemispheric dispersal event may have been of fundamental importance in the early history of crown-group Eutheria. Parallel adaptive radiations have subsequently occurred in the four major groups: Laurasiatheria, Euarchonta + Glires, Afrotheria, and Xenarthra.     

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.     

Phylogenetic utility of different types of molecular data used to infer evolutionary relationships among stalk-eyed flies (Diopsidae)

A phylogenetic hypothesis of relationships among 33 species of stalk-eyed flies was generated from a molecular data set comprising three mitochondrial and three nuclear gene regions. A combined analysis of all the data equally weighted produced a single most-parsimonious cladogram with relatively strong support at the majority of nodes. The phylogenetic utility of different classes of molecular data was also examined. In particular, using a number of different measures of utility in both a combined and separate analysis framework, we focused on the distinction between mitochondrial and nuclear genes and between faster-evolving characters and slower-evolving characters. For the first comparison, by nearly any measure of utility, the nuclear genes are substantially more informative for resolving diopsid relationships than are the mitochondrial genes. The nuclear genes exhibit less homoplasy, are less incongruent with one another and with the combined data, and contribute more support to the combined analysis topology than do the mitochondrial genes. Results from the second comparison, however, provide little evidence of a clear difference in utility. Despite indications of rapid divergence and saturation, faster-evolving characters in both the nuclear and mitochondrial data sets still provide substantial phylogenetic signal. In general, inclusion of the more rapidly evolving data consistently improves the congruence among partitions.     

Base Compositional Bias and Phylogenetic Analyses: A Test of the “Flying DNA” Hypothesis

Phylogenetic methods can produce biased estimates of phylogeny when base composition varies along different lineages. Pettigrew (1994,Curr. Biol.4:277–280) has suggested that base composition bias is responsible for the apparent support for the monophyly of bats (Chiroptera: megabats and microbats) from several different nuclear and mitochondrial genes. Pettigrew's “flying DNA” hypothesis makes several predictions: (1) that metabolic constraints associated with flying result in elevated levels of adenine and thymine throughout the genome of both megabats and microbats, (2) that the resulting base compositional bias in bats is sufficient to mislead phylogenetic methods and account for the support for bat monophyly from several nuclear and mitochondrial genes, and (3) that phylogenetic analysis using pairwise distances corrected for compositional bias should eliminate the support for bat monophyly. We tested these predictions by analyzing DNA sequences from two nuclear and three mitochondrial genes. The predicted base compositional bias does not appear to exist in some of the genes, and in other genes the differences in AT content are very small. Analyses under a wide diversity of criteria and models of evolution, including analyses that take base composition into account (using log-determinant distances), all strongly support bat monophyly. Moreover, simulation analyses indicate that even extreme bias toward AT-base composition in bats would be insufficient to explain the observed levels of support for bat monophyly. These analyses provide no support for the “flying DNA” hypothesis, whereas the monophyly of bats appears to be well supported by the DNA sequence data.     

Data partitions and complex models in Bayesian analysis: the phylogeny of Gymnophthalmid lizards

Phylogenetic studies incorporating multiple loci, and multiple genomes, are becoming increasingly common. Coincident with this trend in genetic sampling, model-based likelihood techniques including Bayesian phylogenetic methods continue to gain popularity. Few studies, however, have examined model fit and sensitivity to such potentially heterogeneous data partitions within combined data analyses using empirical data. Here we investigate the relative model fit and sensitivity of Bayesian phylogenetic methods when alternative site-specific partitions of among-site rate variation (with and without autocorrelated rates) are considered. Our primary goal in choosing a best-fit model was to employ the simplest model that was a good fit to the data while optimizing topology and/or Bayesian posterior probabilities. Thus, we were not interested in complex models that did not practically affect our interpretation of the topology under study. We applied these alternative models to a four-gene data set including one protein-coding nuclear gene (c-mos), one protein-coding mitochondrial gene (ND4), and two mitochondrial rRNA genes (12S and 16S) for the diverse yet poorly known lizard family Gymnophthalmidae. Our results suggest that the best-fit model partitioned among-site rate variation separately among the c-mos, ND4, and 12S + 16S gene regions. We found this model yielded identical topologies to those from analyses based on the GTR+I+G model, but significantly changed posterior probability estimates of clade support. This partitioned model also produced more precise (less variable) estimates of posterior probabilities across generations of long Bayesian runs, compared to runs employing a GTR+I+G model estimated for the combined data. We use this three-way gamma partitioning in Bayesian analyses to reconstruct a robust phylogenetic hypothesis for the relationships of genera within the lizard family Gymnophthalmidae. We then reevaluate the higher-level taxonomic arrangement of the Gymnophthalmidae. Based on our findings, we discuss the utility of nontraditional parameters for modeling among-site rate variation and the implications and future directions for complex model building and testing.     

Molecular phylogenetics of the allodapine bee genus Braunsapis: A-T bias and heterogeneous substitution parameters

Extreme AT bias in Hymenopteran mitochondrial genes have created difficulties for molecular phylogenetic analyses, especially for older divergences where multiple substitutions can erode signal. Heterogeneity in the evolutionary rates of different codon positions and different genes also appears to have been a major problem in resolving ancient divergences in allodapine bees. Here we examine the phylogeny of relatively recent divergences in the allodapine bee genus Braunsapis. We examined heterogeneity in nucleotide substitution parameters for one nuclear gene and codon positions in two mitochondrial genes, exploring various phylogenetic analyses for recovering relationships among species from Africa, Madagascar, southern Asia, and Australia. We explored maximum parsimony, maximum likelihood, Log determinant and Bayesian analyses. Broad topological features of best fit trees tended to be similar for equivalent data sets (e.g., total, or with 3rd mt positions excluded), regardless of the analytic method used (e.g., maximum likelihood or Bayesian). Analyses that used the total data set without modelling partitions separately gave unlikely results, indicating that the Malagasy species was most closely related to Australian species. However, analyses that excluded 3rd mitochondrial positions, or modelled partitions separately, suggested that the Malagasy species falls within the African clade. The unlikely topologies apparently result from long branch attraction, and this problem is ameliorated where modelling allows more realistic estimates of base composition and evolutionary rates for 3rd mitochondrial positions. However, we found that even when codon positions are modelled separately, estimated evolutionary rates for 3rd mitochondrial positions are likely to underestimate true rates. Long branch attraction and multiple substitutions are likely to be much more difficult to circumvent in analyses that explore older, generic-level, divergences in allodapine bees where overwriting is expected to be much more extreme. Our results indicate an African origin for Braunsapis, followed by a single, very early, dispersal event into Asia and then by a later dispersal event into Australia. The Malagasy species is derived from within the African clade.     

Dynamically heterogenous partitions and phylogenetic inference: an evaluation of analytical strategies with cytochrome b and ND6 gene sequences in cranes

ki ctes over whether molecular sequence data should be partitioned for phylogenetic analysis often confound two types of heterogeneity among partitions. We distinguish historical heterogeneity (i.e., different partitions have different evolutionary relationships) from dynamic heterogeneity (i.e., different partitions show different patterns of sequence evolution) and explore the impact of the latter on phylogenetic accuracy and precision with a two-gene, mitochondrial data set for cranes. The well-established phylogeny of cranes allows us to contrast tree-based estimates of relevant parameter values with estimates based on pairwise comparisons and to ascertain the effects of incorporating different amounts of process information into phylogenetic estimates. We show that codon positions in the cytochrome b and NADH dehydrogenase subunit 6 genes are dynamically heterogenous under both Poisson and invariable-sites + gamma-rates versions of the F84 model and that heterogeneity includes variation in base composition and transition bias as well as substitution rate. Estimates of transition-bias and relative-rate parameters from pairwise sequence comparisons were comparable to those obtained as tree-based maximum likelihood estimates. Neither rate-category nor mixed-model partitioning strategies resulted in a loss of phylogenetic precision relative to unpartitioned analyses. We suggest that weighted-average distances provide a computationally feasible alternative to direct maximum likelihood estimates of phylogeny for mixed-model analyses of large, dynamically heterogenous data sets.     

A partitioned likelihood analysis of swallowtail butterfly phylogeny (Lepidoptera:Papilionidae)

Although it is widely agreed that data from multiple sources are necessary to confidently resolve phylogenetic relationships, procedures for accommodating and incorporating heterogeneity in such data remain underdeveloped. We explored the use of partitioned, model-based analyses of heterogeneous molecular data in the context of a phylogenetic study of swallowtail butterflies (Lepidoptera: Papilionidae). Despite substantial basic and applied study, phylogenetic relationships among the major lineages of this prominent group remain contentious. We sequenced 3.3 kb of mitochondrial and nuclear DNA (2.3 kb of cytochrome oxidase I and II and 1.0 kb of elongation factor-1 alpha, respectively) from 22 swallowtails, including representatives of Baroniinae, Parnassiinae, and Papilioninae, and from several moth and butterfly outgroups. Using parsimony, we encountered considerable difficulty in resolving the deepest splits among these taxa. We therefore chose two outgroups with undisputed relationships to each other and to Papilionidae and undertook detailed likelihood analyses of alternative topologies. Following from previous studies that have demonstrated substantial heterogeneity in the evolutionary dynamics among process partitions of these genes, we estimated evolutionary parameters separately for gene-based and codon-based partitions. These values were then used as the basis for examining the likelihoods of possible resolutions and rootings under several partitioned and unpartitioned likelihood models. Partitioned models gave markedly better fits to the data than did unpartitioned models and supported different topologies. However, the most likely topology varied from model to model. The most likely ingroup topology under the best-fitting, six-partition GTR + gamma model favors a paraphyletic Parnassiinae. However, when examining the likelihoods of alternative rootings of this tree relative to rootings of the classical hypothesis, two rootings of the latter emerge as most likely. Of these two, the most likely rooting is within the Papilioninae, although a rooting between Baronia and the remaining Papilionidae is only nonsignificantly less likely.     

Phylogenetic relationships of artiodactyls and cetaceans as deduced from the comparison of cytochrome b and 12S rRNA mitochondrial sequences

A data set of complete mitochondrial cytochrome b and 12S rDNA sequences is presented here for 17 representatives of Artiodactyla and Cetacea, together with potential outgroups (two Perissodactyla, two Carnivora, two Tethytheria, four Rodentia, and two Marsupialia). We include seven sequences not previously published from Hippopotamidae (Ancodonta) and Camelidae (Tylopoda), yielding a total of nearly 2.1 kb for both genes combined. Distance and parsimony analyses of each gene indicate that 11 clades are well supported, including the artiodactyl taxa Pecora, Ruminantia (with low 12S rRNA support), Tylopoda, Suina, and Ancodonta, as well as Cetacea, Perissodactyla, Carnivora, Tethytheria, Muridae, and Caviomorpha. Neither the cytochrome b nor the 12S rDNA genes resolve the relationships between these major clades. The combined analysis of the two genes suggests a monophyletic Cetacea +Artiodactyla clade (defined as "Cetartiodactyla"), whereas Perissodactyla, Carnivora, and Tethytheria fall outside this clade. Perissodactyla could represent the sister taxon of Cetartiodactyla, as deduced from resampling studies among outgroup lineages. Cetartiodactyla includes five major lineages: Ruminantia, Tylopoda, Suina, Ancodonta, and Cetacea, among which the phylogenetic relationships are not resolved. Thus, Suiformes do not appear to be monophyletic, justifying their split into the Suina and Ancodonta infraorders. An association between Cetacea and Hippopotamidae is supported by the cytochrome b gene but not by the 12S rRNA gene. Calculation of divergence dates suggests that the Cetartiodactyla could have diverged from other Ferungulata about 60 MYA.      

Searching for the closest living relative(s) of tetrapods through evolutionary analyses of mitochondrial and nuclear data

The phylogenetic relationships of the African lungfish (Protopterus dolloi) and the coelacanth (Latimeria chalumnae) with respect to tetrapods were analyzed using complete mitochondrial genome DNA sequences. A lungfish + coelancanth clade was favored by maximum parsimony (although this result is dependent on which transition:transversion weights are applied), and a lungfish + tetrapod clade was supported by neighbor-joining and maximum-likelihood analyses. These two hypotheses received the strongest statistical and bootstrap support to the exclusion of the third alternative, the coelacanth + tetrapod sister group relationship. All mitochondrial protein coding genes combined favor a lungfish + tetrapod grouping. We can confidently reject the hypothesis that the coelacanth is the closest living relative of tetrapods. When the complete mitochondrial sequence data were combined with nuclear 28S rRNA gene data, a lungfish + coelacanth clade was supported by maximum parsimony and maximum likelihood, but a lungfish + tetrapod clade was favored by neighbor-joining. The seeming conflicting results based on different data sets and phylogenetic methods were typically not statistically strongly supported based on Kishino-Hasegawa and Templeton tests, although they were often supported by strong bootstrap values. Differences in rate of evolution of the different mitochondrial genes (slowly evolving genes such as the cytochrome oxidase and tRNA genes favored a lungfish + coelacanth clade, whereas genes of relatively faster substitution rate, such as several NADH dehydrogenase genes, supported a lungfish + tetrapod grouping), as well as the rapid radiation of the lineages back in the Devonian, rather than base compositional biases among taxa seem to be directly responsible for the remaining uncertainty in accepting one of the two alternate hypotheses.     

Shaking the Diptera tree of life: performance analysis of nuclear and mitochondrial sequence data partitions

The value of mitochondrial versus nuclear gene sequence data in phylogenetic analysis has received much attention without yielding definitive conclusions. Theoretical arguments and empirical data suggest a lower phylogenetic utility than equivalent nuclear gene sequences, but there are many examples of important progress made using mitochondrial sequences. We undertook a systematic performance analysis of mitochondrial and nuclear sequence partitions taken from a representative sample of dipteran species. When analysed alone, mitochondrial genes generally performed less well than nuclear genes; however, these genes resolved some branches for which nuclear genes failed. Moreover, the combined use of mitochondrial and nuclear sequences produced superior results without artifacts for nodes where mitochondrial and nuclear gene data generated conflicting topologies. These findings strongly advocate the inclusion of mitochondrial sequences, even in deep phylogeny reconstruction. Comparison of tree support between our and previous analyses identified robustly supported high‐confidence clades in the Diptera, but also revealed problematic groupings in need of further analysis.