Bayesian phylogenetics using an RNA substitution model applied to early mammalian evolution

We study the phylogeny of the placental mammals using molecular data from all mitochondrial tRNAs and rRNAs of 54 species. We use probabilistic substitution models specific to evolution in base paired regions of RNA. A number of these models have been implemented in a new phylogenetic inference software package for carrying out maximum likelihood and Bayesian phylogenetic inferences. We describe our Bayesian phylogenetic method which uses a Markov chain Monte Carlo algorithm to provide samples from the posterior distribution of tree topologies. Our results show support for four primary mammalian clades, in agreement with recent studies of much larger data sets mainly comprising nuclear DNA. We discuss some issues arising when using Bayesian techniques on RNA sequence data.     

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.     

Pleistocene diversification of living squirrel monkeys (Saimiri spp.) inferred from complete mitochondrial genome sequences

In order to enhance our understanding of the evolutionary history of squirrel monkeys (Saimiri spp.), we newly sequenced and analyzed data from seven complete mitochondrial genomes representing six squirrel monkey taxa. While previous studies have lent insights into the taxonomy and phylogeny of the genus, phylogenetic relationships and divergence date estimates among major squirrel monkey clades remain unclear. Using maximum likelihood and Bayesian procedures, we inferred a highly resolved phylogenetic tree with strong support for a sister relationship between Saimiri boliviensis and all other Saimiri, for monophyly of Saimiri oerstedii and Saimiri sciureus sciureus, and for Saimiri sciureus macrodon as the sister lineage to the S. oerstedii/S. s. sciureus clade. We inferred that crown lineages for extant squirrel monkeys diverged around 1.5 million years ago (MYA) in the Pleistocene Epoch, with other major clades diverging between 0.9 and 1.1 MYA. Our results suggest a relatively recent timeline of squirrel monkey evolution and challenge previous conceptions about the diversification of the genus and its expansion into Central America.     

Phylogenetics in Caenorhabditis elegans: an analysis of divergence and outcrossing

This study establishes a phylogenetic framework for the natural geographic isolates of the widely studied nematode species Caenorhabditis elegans. Virtually complete mitochondrial genomes are sequenced from 27 C. elegans natural isolates to characterize mitochondrial divergence patterns and to investigate the evolutionary history of the C. elegans hermaphrodite lineages. Phylogenetic analysis of mitochondrial sequences reveals the presence of two major C. elegans hermaphrodite clades (designated clade I and clade II). Fifty-six nuclear loci, widely distributed across the five autosomes and the X chromosome, are also analyzed in a subset of the C. elegans isolates to evaluate nuclear divergence patterns and the extent of mating between different strains. A comparison of the phylogenetic tree derived from mitochondrial data with the phylogenetic tree derived from nuclear data reveals only one inconsistency in the distribution of isolates into clades I and II, suggesting that mating between divergent C. elegans strains is an infrequent event in the wild.     

Complete mitochondrial DNA genome sequences show that modern birds are not descended from transitional shorebirds

To test the hypothesis put forward by Feduccia of the origin of modern birds from transitional birds, we sequenced the first two complete mitochondrial genomes of shorebirds (ruddy turnstone and blackish oystercatcher) and compared their sequences with those of already published avian genomes. When corrected for rate heterogeneity across sites and non-homogeneous nucleotide compositions among lineages in maximum likelihood (ML), the optimal tree places palaeognath birds as sister to the neognaths including shorebirds. This optimal topology is a re-rooting of recently published ordinal-level avian trees derived from mitochondrial sequences. Using a penalized likelihood (PL) rate-smoothing process in conjunction with dates estimated from fossils, we show that the basal splits in the bird tree are much older than the Cretaceous-Tertiary (K-T) boundary, reinforcing previous molecular studies that rejected the derivation of modern birds from transitional shorebirds. Our mean estimate for the origin of modern birds at about 123 million years ago (Myr ago) is quite close to recent estimates using both nuclear and mitochondrial genes, and supports theories of continental break-up as a driving force in avian diversification. Not only did many modern orders of birds originate well before the K-T boundary, but the radiation of major clades occurred over an extended period of at least 40 Myr ago, thus also falsifying Feduccia's rapid radiation scenario following a K-T bottleneck.     

A bayesian analysis of metazoan mitochondrial genome arrangements

Genome arrangements are a potentially powerful source of information to infer evolutionary relationships among distantly related taxa. Mitochondrial genome arrangements may be especially informative about metazoan evolutionary relationships because (1) nearly all animals have the same set of definitively homologous mitochondrial genes, (2) mitochondrial genome rearrangement events are rare relative to changes in sequences, and (3) the number of possible mitochondrial genome arrangements is huge, making convergent evolution of genome arrangements appear highly unlikely. In previous studies, phylogenetic evidence in genome arrangement data is nearly always used in a qualitative fashion-the support in favor of clades with similar or identical genome arrangements is considered to be quite strong, but is not quantified. The purpose of this article is to quantify the uncertainty among the relationships of metazoan phyla on the basis of mitochondrial genome arrangements while incorporating prior knowledge of the monophyly of various groups from other sources. The work we present here differs from our previous work in the statistics literature in that (1) we incorporate prior information on classifications of metazoans at the phylum level, (2) we describe several advances in our computational approach, and (3) we analyze a much larger data set (87 taxa) that consists of each unique, complete mitochondrial genome arrangement with a full complement of 37 genes that were present in the NCBI (National Center for Biotechnology Information) database at a recent date. In addition, we analyze a subset of 28 of these 87 taxa for which the non-tRNA mitochondrial genomes are unique where the assumption of our inversion-only model of rearrangement is more plausible. We present summaries of Bayesian posterior distributions of tree topology on the basis of these two data sets.     

A comprehensive analysis of mammalian mitochondrial genome base composition and improved phylogenetic methods

Phylogenetic analysis of mammalian species using mitochondrial protein genes has proved to be problematic in many previous studies. The high mutation rate of mitochondrial DNA and unusual base composition of several species has prompted us to conduct a detailed study of the composition of 69 mammalian mitochondrial genomes. Most major changes in base composition between lineages can be attributed to shifts between the proportions of C and T on the L-strand. These changes are significant at all codon positions and are shown to affect amino acid composition. Correlated changes in the base composition of the RNA loops and stems are also observed. Following up from previous studies, we investigate changes in the base composition of all 12 H-strand proteins and find that variability in proportions of C and T is correlated with location on the genome. Variation in base composition across genes and species is known to adversely affect the performance of phylogenetic inference methods. We have, therefore, developed a customized three-state general time-reversible DNA substitution model, implemented in the PHASE phylogenetic inference package, which lumps C and T into a composite pyrimidine state. We compare the phylogenetic tree obtained using the new three-state model with that obtained using a standard four-state model. Results using the three-state model are more congruent with recent studies using large sets of nuclear genes and help resolve some of the apparent conflicts between studies using nuclear and mitochondrial proteins.     

Long-branch attraction phenomenon and the impact of among-site rate variation on rodent phylogeny

The phylogenetic relationships among major lineages of rodents is one of the issues most debated by both paleontologists and molecular biologists. In the present study, we have analyzed all complete mammalian mitochondrial genomes available in the databases, including five rodent species (rat, mouse, dormouse, squirrel and guinea-pig). Phylogenetic analyses were performed on H-strand amino acid sequences by means of maximum-likelihood and on H-strand protein-coding and ribosomal genes by means of distance methods. Also, log-likelihood ratio tests were applied to different tree topologies under the assumption of rodent monophyly, paraphyly or polyphyly. The analyses significantly rejected rodent monophyly and showed the existence of two differentiated clades, one containing non-murids (dormouse, squirrel and guinea-pig) and the other containing murids (rat and mouse). Long-branch attraction between murids and the outgroups could not be responsible for the existence of two different rodent clades, as no significant differences in evolutionary rate have been observed, except in the case of the squirrel, which shows a lower rate. The impact of among-site rate variation models on the phylogeny of rodents has been evaluated using the gamma distribution model. Results have shown that relationships among rodents remained unchanged, and the general topology of the tree was not affected, even though some branches were not properly resolved, most likely due to a lack of fit between estimated and real rate heterogeneity parameters.     

Mitochondrial data are not suitable for resolving placental mammal phylogeny

Mitochondrial data have traditionally been used in reconstructing a variety of species phylogenies. The low rates of recombination and thorough characterization of mitochondrial data across vertebrate species make it a particularly attractive phylogenetic marker. The relatively low number of fully sequenced mammal genomes and the lack of extensive sampling within Superorders have posed a serious problem for reaching agreement on the placement mammal species. The use of mitochondrial data sequences from large numbers of mammals could serve to circumvent the taxon-sampling deficit. Here we assess the suitability of mitochondrial data as a phylogenetic marker in mammal phylogenetics. MtDNA datasets of mammal origin have been filtered as follows: (i) we have sampled sparsely across the phylogenetic tree, (ii) we have constrained our sampling to genes with high taxon coverage, (iii) we have categorised rates across sites in a phylogeny independent manner and have removed fast evolving sites, and (iv), we have sampled from very shallow divergence times to reduce phylogenetic conflict. However, topologies obtained using these filters are not consistent with previous studies and are discordant across different genes. Individual mitochondrial genes, and indeed all mitochondrial genes analysed as a supermatrix, resulted in poor resolution of the species phylogeny. Overall, our study highlights the limitations of mitochondrial data, not only for resolving deep divergences and but also for shallow divergences in the mammal phylogeny.     

Mitochondrial versus nuclear gene sequences in deep-level mammalian phylogeny reconstruction

Both mitochondrial and nuclear gene sequences have been employed in efforts to reconstruct deep-level phylogenetic relationships. A fundamental question in molecular systematics concerns the efficacy of different types of sequences in recovering clades at different taxonomic levels. We compared the performance of four mitochondrial data sets (cytochrome b, cytochrome oxidase II, NADH dehydrogenase subunit I, 12S rRNA-tRNA-16S rRNA) and eight nuclear data sets (exonic regions of alpha-2B adrenergic receptor, aquaporin, ss-casein, gamma-fibrinogen, interphotoreceptor retinoid binding protein, kappa-casein, protamine, von Willebrand Factor) in recovering deep-level mammalian clades. We employed parsimony and minimum-evolution with a variety of distance corrections for superimposed substitutions. In 32 different pairwise comparisons between these mitochondrial and nuclear data sets, we used the maximum set of overlapping taxa. In each case, the variable-length bootstrap was used to resample at the size of the smaller data set. The nuclear exons consistently performed better than mitochondrial protein and rRNA-tRNA coding genes on a per-residue basis in recovering benchmark clades. We also concatenated nuclear genes for overlapping taxa and made comparisons with concatenated mitochondrial protein-coding genes from complete mitochondrial genomes. The variable-length bootstrap was used to score the recovery of benchmark clades as a function of the number of resampled base pairs. In every case, the nuclear concatenations were more efficient than the mitochondrial concatenations in recovering benchmark clades. Among genes included in our study, the nuclear genes were much less affected by superimposed substitutions. Nuclear genes having appropriate rates of substitution should receive strong consideration in efforts to reconstruct deep-level phylogenetic relationships.     

Phylogeny and evolution of Cervidae based on complete mitochondrial genomes

Mitochondrial DNA sequences can be used to estimate phylogenetic relationships among animal taxa and for molecular phylogenetic evolution analysis. With the development of sequencing technology, more and more mitochondrial sequences have been made available in public databases, including whole mitochondrial DNA sequences. These data have been used for phylogenetic analysis of animal species, and for studies of evolutionary processes. We made phylogenetic analyses of 19 species of Cervidae, with Bos taurus as the outgroup. We used neighbor joining, maximum likelihood, maximum parsimony, and Bayesian inference methods on whole mitochondrial genome sequences. The consensus phylogenetic trees supported monophyly of the family Cervidae; it was divided into two subfamilies, Plesiometacarpalia and Telemetacarpalia, and four tribes, Cervinae, Muntiacinae, Hydropotinae, and Odocoileinae. The divergence times in these families were estimated by phylogenetic analysis using the Bayesian method with a relaxed molecular clock method; the results were consistent with those of previous studies. We concluded that the evolutionary structure of the family Cervidae can be reconstructed by phylogenetic analysis based on whole mitochondrial genomes; this method could be used broadly in phylogenetic evolutionary analysis of animal taxa.     

Mitochondrial genomes of two luminous beetles, Rhagophthalmus lufengensis and R. ohbai (Arthropoda, Insecta, Coleoptera)

We determined the mitochondrial DNA (mtDNA) sequences of two luminous beetles (Arthropoda, Insecta, Coleoptera), Rhagophthalmus lufengensis from Yunnan, China and Rhagophthalmus ohbai from Yaeyama Island, Japan. We identified all the 37 mtDNA genes of R. lufengensis (15,982 bp) and the 34 genes of R. ohbai (15,704 bp). R. lufengensis and R. ohbai genomes have higher A + T contents than other coleopteran genomes although the gene arrangements are similar. Interestingly, in a study of the evolutionary relationship among R. lufengensis, R. ohbai and the firefly Pyrocoelia rufa, the phylogenetic tree inferred from lrRNA genes from mitochondrial genomes indicates a biogeographic relationship among the bioluminescent insects in East Asia and the phylogenetic tree inferred from luciferase-related genes from nuclear genomes shows an appropriate relationship among coleopterans, reflecting the evolutionary origin of bioluminescence. Thus, the mtDNAs of luminescent beetles can provide an insight into their evolutionary origin and biogeography.     

Arthropod Phylogenetics in Light of Three Novel Millipede (Myriapoda: Diplopoda) Mitochondrial Genomes with Comments on the Appropriateness of Mitochondrial Genome Sequence Data for Inferring Deep Level Relationships

Background

Arthropods are the most diverse group of eukaryotic organisms, but their phylogenetic relationships are poorly understood. Herein, we describe three mitochondrial genomes representing orders of millipedes for which complete genomes had not been characterized. Newly sequenced genomes are combined with existing data to characterize the protein coding regions of myriapods and to attempt to reconstruct the evolutionary relationships within the Myriapoda and Arthropoda.

Results

The newly sequenced genomes are similar to previously characterized millipede sequences in terms of synteny and length. Unique translocations occurred within the newly sequenced taxa, including one half of the Appalachioria falcifera genome, which is inverted with respect to other millipede genomes. Across myriapods, amino acid conservation levels are highly dependent on the gene region. Additionally, individual loci varied in the level of amino acid conservation. Overall, most gene regions showed low levels of conservation at many sites. Attempts to reconstruct the evolutionary relationships suffered from questionable relationships and low support values. Analyses of phylogenetic informativeness show the lack of signal deep in the trees (i.e., genes evolve too quickly). As a result, the myriapod tree resembles previously published results but lacks convincing support, and, within the arthropod tree, well established groups were recovered as polyphyletic.

Conclusions

The novel genome sequences described herein provide useful genomic information concerning millipede groups that had not been investigated. Taken together with existing sequences, the variety of compositions and evolution of myriapod mitochondrial genomes are shown to be more complex than previously thought. Unfortunately, the use of mitochondrial protein-coding regions in deep arthropod phylogenetics appears problematic, a result consistent with previously published studies. Lack of phylogenetic signal renders the resulting tree topologies as suspect. As such, these data are likely inappropriate for investigating such ancient relationships.     

A mitochondrial genome phylogeny of Diptera: whole genome sequence data accurately resolve relationships over broad timescales with high precision

Abstract Mitochondrial genomes provide a promising new tool for understanding deep‐level insect phylogenetics, but have yet to be evaluated for their ability to resolve intraordinal relationships. We tested the utility of mitochondrial genome data for the resolution of relationships within Diptera, the insect order for which the most data are available. We sequenced an additional three genomes, from a syrphid, nemestrinid and tabanid, representing three additional dipteran clades, ‘aschiza’, non‐heteroneuran muscomorpha and ‘basal brachyceran’, respectively. We assessed the influence of optimality criteria, gene inclusion/exclusion, data recoding and partitioning strategies on topology and nodal support within Diptera. Our consensus phylogeny of Diptera was largely consistent with previous phylogenetic hypotheses of the order, except that we did not recover a monophyletic Muscomorpha (Nesmestrinidae grouped with Tabanidae) or Acalyptratae (Drosophilidae grouped with Calliphoridae). The results were very robust to optimality criteria, as parsimony, likelihood and Bayesian approaches yielded very similar topologies, although nodal support varied. The addition of ribosomal and transfer RNA genes to the protein coding genes traditionally used in mitochondrial genome phylogenies improved the resolution and support, contrary to previous suggestions that these genes would evolve too quickly or prove too difficult to align to provide phylogenetic signal at deep nodes. Strategies to recode data, aimed at reducing homoplasy, resulted in a decrease in tree resolution and branch support. Bayesian analyses were highly sensitive to partitioning strategy: biologically realistic partitions into codon groups produced the best results. The implications of this study for dipteran systematics and the effective approaches to using mitochondrial genome data are discussed. Mitochondrial genomes resolve intraordinal relationships within Diptera accurately over very wide time ranges (1–200 million years ago) and genetic distances, suggesting that this may be an excellent data source for deep‐level studies within other, less studied, insect orders.     

Complete mitochondrial DNA sequences of the green turtle and blue-tailed mole skink: statistical evidence for archosaurian affinity of turtles.

Turtles have highly specialized morphological characteristics, and their phylogenetic position has been under intensive debate. Previous molecular studies have not established a consistent and statistically well supported conclusion on this issue. In order to address this, complete mitochondrial DNA sequences were determined for the green turtle and the blue-tailed mole skink. These genomes possess an organization of genes which is typical of most other vertebrates, such as placental mammals, a frog, and bony fishes, but distinct from organizations of alligators and snakes. Molecular evolutionary rates of mitochondrial protein sequences appear to vary considerably among major reptilian lineages, with relatively rapid rates for snake and crocodilian lineages but slow rates for turtle and lizard lineages. In spite of this rate heterogeneity, phylogenetic analyses using amino acid sequences of 12 mitochondrial proteins reliably established the Archosauria (birds and crocodilians) and Lepidosauria (lizards and snakes) clades postulated from previous morphological studies. The phylogenetic analyses further suggested that turtles are a sister group of the archosaurs, and this untraditional relationship was provided with strong statistical evidence by both the bootstrap and the Kishino-Hasegawa tests. This is the first statistically significant molecular phylogeny on the placement of turtles relative to the archosaurs and lepidosaurs. It is therefore likely that turtles originated from a Permian-Triassic archosauromorph ancestor with two pairs of temporal fenestrae behind the skull orbit that were subsequently lost. The traditional classification of turtles in the Anapsida may thus need to be reconsidered.     

Whole genome-based phylogenetic analysis of free-living microorganisms.

A phylogenetic 'tree of life' has been constructed based on the observed presence and absence of families of protein-encoding genes observed in 11 complete genomes of free-living microorganisms. Past attempts to reconstruct the evolutionary relation-ships of microorganisms have been limited to sets of genes rather than complete genomes. Despite apparent rampant lateral gene transfer among microorganisms, these results indicate a single robust underlying evolutionary history for these organisms. Broadly, the tree produced is very similar to the small subunit rRNA tree although several additional phylogenetic relationships appear to be resolved, including the relationship of Archaeoglobus to the methanogens studied. This result is in contrast to notions that a robust phylogenetic reconstruction of microorganisms is impossible due to their genomes being composed of an incomprehensible amalgam of genes with complicated histories and suggests that this style of genome-wide phylogenetic analysis could become an important method for studying the ancient diversification of life on Earth. Analyses using informational and operational subsets of the genes showed that this 'tree of life' is not dependent on the phylogenetically more consistent informational genes.     

Combined nuclear and mitochondrial DNA sequences resolve generic relationships within the Cracidae (Galliformes,Aves)

The Cracidae is one of the most endangered and distinctive bird families in the Neotropics, yet the higher relationships among taxa remain uncertain. The molecular phylogeny of its 11 genera was inferred using 10,678 analyzable sites (5,412 from seven different mitochondrial segments and 5,266 sites from four nuclear genes). We performed combinability tests to check conflicts in phylogenetic signals of separate genes and genomes. Phylogenetic analysis showed that the unrooted tree of ((curassows, horned guan) (guans, chachalacas)) was favored by most data partitions and that different data partitions provided support for different parts of the tree. In particular, the concatenated mitochondrial DNA (mtDNA) genes resolved shallower nodes, whereas the combined nuclear sequences resolved the basal connections among the major clades of curassows, horned guan, chachalacas, and guans. Therefore, we decided that for the Cracidae all data should be combined for phylogenetic analysis. Maximum parsimony (MP), maximum likelihood (ML), and Bayesian analyses of this large data set produced similar trees. The MP tree indicated that guans are the sister group to (horned guan, (curassows, chachalacas)), whereas the ML and Bayesian analysis recovered a tree where the horned guan is a sister clade to curassows, and these two clades had the chachalacas as a sister group. Parametric bootstrapping showed that alternative trees previously proposed for the cracid genera are significantly less likely than our estimate of their relationships. A likelihood ratio test of the hypothesis of a molecular clock for cracid mtDNA sequences using the optimal ML topology did not reject rate constancy of substitutions through time. We estimated cracids to have originated between 64 and 90 million years ago (MYA), with a mean estimate of 76 MYA. Diversification of the genera occurred approximately 41-3 MYA, corresponding with periods of global climate change and other Earth history events that likely promoted divergences of higher level taxa.     

Congruent mammalian trees from mitochondrial and nuclear genes using Bayesian methods

Analyses of mitochondrial and nuclear gene sequences have often produced different mammalian tree topologies, undermining confidence in the merit of molecular approaches with respect to "traditional" morphological classification. The recent sequencing of the complete mitochondrial genomes of two additional rodents (Spalax judaei and Jaculus jaculus) and one lagomorph (Ochotona princeps) has prompted us to reinvestigate the issue. Using Bayesian phylogenetics, we found phylogenetic relationships between mammalian species highly congruent with previous results based on nuclear genes. Our results show the existence of four primary lineages of placental mammals: Xenarthra, Afrotheria, Laurasiatheria, and Euarchontoglires. Relationships between and within these lineages strongly suggest that the gene trees may also be congruent with the underlying species phylogeny.     

The complete mitochondrial genome of the mantid shrimp Harpiosquilla harpax, and a phylogenetic investigation of the Decapoda using mitochondrial sequences

The complete mitochondrial DNA sequence was determined for the mantid shrimp Harpiosquilla harpax. These data demonstrate that the H. harpax mitochondrial genome is a 15,714 bp circular molecule and encodes the typical 37 metazoan mitochondrial genes (13 protein-coding, 22 tRNA, and two rRNA genes). The gene arrangement of H. harpax is consistent with that of the putative arthropod ancestral gene order as depicted by Limulus polyphemus. H. harpax was employed as an outgroup taxon for a phylogenetic investigation of the Decapoda using sequences from complete mitochondrial genomes. Whilst our results are largely in agreement with current taxonomic treatments, the relationships indicated among the reptantian decapods are novel. Our results provide strong statistical support for a sister-group relationship between the Achelata and the Astacida. These findings not only refute previous phylogenetic hypotheses, but also have serious implications for the interpretation of morphological and developmental evolution in the Decapoda. In addition we also investigated the effects of outgroup selection on the resolution of ingroup relationships. We found outgroup choice to significantly influence tree topology thus reinforcing the importance of appropriate outgroup selection in phylogenetic studies.     

Sampling properties of DNA sequence data in phylogenetic analysis

We inferred phylogenetic trees from individual genes and random samples of nucleotides from the mitochondrial genomes of 10 vertebrates and compared the results to those obtained by analyzing the whole genomes. Individual genes are poor samples in that they infrequently lead to the whole-genome tree. A large number of nucleotide sites is needed to exactly determine the whole-genome tree. A relatively small number of sites, however, often results in a tree close to the whole-genome tree. We found that blocks of contiguous sites were less likely to lead to the whole-genome tree than samples composed of sites drawn individually from throughout the genome. Samples of contiguous sites are not representative of the entire genome, a condition that violates a basic assumption of the bootstrap method as it is applied in phylogenetic studies.