Corroboration versus "Strongest Evidence"

Background knowledge comprises accepted (well-corroborated) theories and results. Such theories are taken to be true for the purpose of interpreting evidence when assessing the corroboration of a hypothesis currently in question. Accordingly, background knowledge does not properly include rejected theories, false assumptions, or null models. In particular, regarding a model of random character distribution as "background knowledge" would rule out corroboration of phylogenetic hypotheses, since it would make character data irrelevant to inferring phylogeny. The presence of homoplasy is not grounds for treating characters as if they were randomly distributed, since characters can show strong phylogenetic structure even when they show homoplasy. This means that clique (compatibility) analysis is unjustified, since that method depends crucially on the assumption that characters showing any homoplasy at all are unrelated to phylogeny. Although likelihood does not measure corroboration, corroboration is closely connected to likelihood: for given evidence and background, the most likely trees are also best corroborated. Most parsimonious trees are best corroborated; the apparent clash between parsimony and likelihood is an artifact of the use of unrealistic models in most "maximum likelihood" methods.     

Failed refutations: further comments on parsimony and likelihood methods and their relationship to Popper's degree of corroboration

Kluge's (2001, Syst. Biol. 50:322-330) continued arguments that phylogenetic methods based on the statistical principle of likelihood are incompatible with the philosophy of science described by Karl Popper are based on false premises related to Kluge's misrepresentations of Popper's philosophy. Contrary to Kluge's conjectures, likelihood methods are not inherently verificationist; they do not treat every instance of a hypothesis as confirmation of that hypothesis. The historical nature of phylogeny does not preclude phylogenetic hypotheses from being evaluated using the probability of evidence. The low absolute probabilities of hypotheses are irrelevant to the correct interpretation of Popper's concept termed degree of corroboration, which is defined entirely in terms of relative probabilities. Popper did not advocate minimizing background knowledge; in any case, the background knowledge of both parsimony and likelihood methods consists of the general assumption of descent with modification and additional assumptions that are deterministic, concerning which tree is considered most highly corroborated. Although parsimony methods do not assume (in the sense of entailing) that homoplasy is rare, they do assume (in the sense of requiring to obtain a correct phylogenetic inference) certain things about patterns of homoplasy. Both parsimony and likelihood methods assume (in the sense of implying by the manner in which they operate) various things about evolutionary processes, although violation of those assumptions does not always cause the methods to yield incorrect phylogenetic inferences. Test severity is increased by sampling additional relevant characters rather than by character reanalysis, although either interpretation is compatible with the use of phylogenetic likelihood methods. Neither parsimony nor likelihood methods assess test severity (critical evidence) when used to identify a most highly corroborated tree(s) based on a single method or model and a single body of data; however, both classes of methods can be used to perform severe tests. The assumption of descent with modification is insufficient background knowledge to justify cladistic parsimony as a method for assessing degree of corroboration. Invoking equivalency between parsimony methods and likelihood models that assume no common mechanism emphasizes the necessity of additional assumptions, at least some of which are probabilistic in nature. Incongruent characters do not qualify as falsifiers of phylogenetic hypotheses except under extremely unrealistic evolutionary models; therefore, justifications of parsimony methods as falsificationist based on the idea that they minimize the ad hoc dismissal of falsifiers are questionable. Probabilistic concepts such as degree of corroboration and likelihood provide a more appropriate framework for understanding how phylogenetics conforms with Popper's philosophy of science. Likelihood ratio tests do not assume what is at issue but instead are methods for testing hypotheses according to an accepted standard of statistical significance and for incorporating considerations about test severity. These tests are fundamentally similar to Popper's degree of corroboration in being based on the relationship between the probability of the evidence e in the presence versus absence of the hypothesis h, i.e., between p(e|hb) and p(e|b), where b is the background knowledge. Both parsimony and likelihood methods are inductive in that their inferences (particular trees) contain more information than (and therefore do not follow necessarily from) the observations upon which they are based; however, both are deductive in that their conclusions (tree lengths and likelihoods) follow necessarily from their premises (particular trees, observed character state distributions, and evolutionary models). For these and other reasons, phylogenetic likelihood methods are highly compatible with Karl Popper's philosophy of science and offer several advantages over parsimony methods in this context.     

Analysis of Acorus calamus chloroplast genome and its phylogenetic implications

Determining the phylogenetic relationships among the major lines of angiosperms is a long-standing problem, yet the uncertainty as to the phylogenetic affinity of these lines persists. While a number of studies have suggested that the ANITA (Amborella-Nymphaeales-Illiciales-Trimeniales-Aristolochiales) grade is basal within angiosperms, studies of complete chloroplast genome sequences also suggested an alternative tree, wherein the line leading to the grasses branches first among the angiosperms. To improve taxon sampling in the existing chloroplast genome data, we sequenced the chloroplast genome of the monocot Acorus calamus. We generated a concatenated alignment (89,436 positions for 15 taxa), encompassing almost all sequences usable for phylogeny reconstruction within spermatophytes. The data still contain support for both the ANITA-basal and grasses-basal hypotheses. Using simulations we can show that were the ANITA-basal hypothesis true, parsimony (and distance-based methods with many models) would be expected to fail to recover it. The self-evident explanation for this failure appears to be a long-branch attraction (LBA) between the clade of grasses and the out-group. However, this LBA cannot explain the discrepancies observed between tree topology recovered using the maximum likelihood (ML) method and the topologies recovered using the parsimony and distance-based methods when grasses are deleted. Furthermore, the fact that neither maximum parsimony nor distance methods consistently recover the ML tree, when according to the simulations they would be expected to, when the out-group (Pinus) is deleted, suggests that either the generating tree is not correct or the best symmetric model is misspecified (or both). We demonstrate that the tree recovered under ML is extremely sensitive to model specification and that the best symmetric model is misspecified. Hence, we remain agnostic regarding phylogenetic relationships among basal angiosperm lineages.     

Molecular phylogeny of the marmots (Rodentia: Sciuridae): tests of evolutionary and biogeographic hypotheses

There are 14 species of marmots distributed across the Holarctic, and despite extensive systematic study, their phylogenetic relationships remain largely unresolved. In particular, comprehensive studies have been lacking. A well-supported phylogeny is needed to place the numerous ecological and behavioral studies on marmots in an evolutionary context. To address this situation, we obtained complete cytochrome (cyt) b sequences for 13 of the species and a partial sequence for the 14th. We applied a statistical approach to both phylogeny estimation and hypothesis testing, using parsimony and maximum likelihood-based methods. We conducted statistical tests on a suite of previously proposed hypotheses of phylogenetic relationships and biogeographic histories. The cyt b data strongly support the monophyly of Marmota and a western montane clade in the Nearctic. Although some other scenarios cannot be rejected, the results are consistent with an initial diversification in North America, followed by an invasion and subsequent rapid diversification in the Palearctic. These analyses reject the two major competing hypotheses of M. broweri's phylogenetic relationships--namely, that it is the sister species to M. camtschatica of eastern Siberia, and that it is related closely to M. caligata of the Nearctic. The Alaskan distribution of M. broweri is best explained as a reinvasion from the Palearctic, but a Nearctic origin can not be rejected. Several other conventionally recognized species groups can also be rejected. Social evolution has been homoplastic, with large colonial systems evolving in two groups convergently. The cyt b data do not provide unambiguous resolution of several basal nodes in the Palearctic radiation, leaving some aspects of pelage and karyotypic evolution equivocal.     

Phylogenetic position of cetaceans relative to artiodactyls: reanalysis of mitochondrial and nuclear sequences

By a maximum likelihood analysis of mitochondrial DNA sequences, we examine Graur and Higgins' hypothesis of the Ruminantia/Cetacea clade with Suiformes as an outgroup. Graur and Higgins analyzed these sequences by the neighbor-joining and parsimony methods, as well as by the maximum likelihood method under the assumption that the substitution rate is the same for all sites. The Ruminantia/Suiformes clade assumed by the traditional taxonomy was rejected strongly by this analysis and the Ruminantia/Cetacea clade was supported. Adoption of a more realistic model distinguishing among rates at different codon positions in the maximum likelihood analysis of the same data, however, grossly reduces the significance level on the Graur-Higgins hypothesis. Thus, although the Ruminantia/Suiformes grouping is indeed least likely from Graur and Higgins' data set of mitochondrial DNA, this traditional tree cannot be rejected with statistical significance under the new analysis, and more data are needed to settle the issue. In the same way, we examine Irwin and Arnason's suggestion of the Hippopotamus/Cetacea clade by using cytochrome b and hemoglobins alpha and beta, and it turn out that their suggestion is also fragile. This analysis demonstrates the importance of selecting an appropriate model among the alternatives in the maximum likelihood analysis and of using many different genes from many relevant species in order to make reliable phylogenetic inferences.      

Phylogenomics and the number of characters required for obtaining an accurate phylogeny of eukaryote model species

MOTIVATION: Through the most extensive phylogenomic analysis carried out to date, complete genomes of 11 eukaryotic species have been examined in order to find the homologous of more than 25,000 amino acid sequences. These sequences correspond to the exons of more than 3000 genes and were used as presence/absence characters to test one of the most controversial hypotheses concerning animal evolution, namely the Ecdysozoa hypothesis. Distance, maximum parsimony and Bayesian methods of phylogenetic reconstruction were used to test the hypothesis. RESULTS: The reliability of the ecdysozoa, grouping arthropods and nematodes in a single clade was unequivocally rejected in all the consensus trees. The Coelomata clade, grouping arthropods and chordates, was supported by the highest statistical confidence in all the reconstructions. The study of the dependence of the genomes' tree accuracy on the number of exons used, demonstrated that an unexpectedly larger number of characters are necessary to obtain robust phylogenies. Previous studies supporting ecdysozoa, could not guarantee an accurate phylogeny because the number of characters used was clearly below the minimum required.     

The evolution of generalization? Parasitoid flies and the perils of inferring host range evolution from phylogenies

It is widely assumed that high resource specificity predisposes lineages toward greater likelihood of extinction and lower likelihood of diversification than more generalized lineages. This suggests that host range evolution in parasitic organisms should proceed from generalist to specialist, and specialist lineages should be found at the 'tips' of phylogenies. To test these hypotheses, parsimony and maximum likelihood methods were used to reconstruct the evolution of host range on a phylogeny of parasitoid flies in the family Tachinidae. In contrast to predictions, most reconstructions indicated that generalists were repeatedly derived from specialist lineages and tended to occupy terminal branches of the phylogeny. These results are critically examined with respect to hypotheses concerning the evolution of specialization, the inherent difficulties in inferring host ranges, our knowledge of tachinid-host associations, and the methodological problems associated with ancestral character state reconstruction. Both parsimony and likelihood reconstructions are shown to provide misleading results and it is argued that independent evidence, in addition to phylogenetic trees, is needed to inform models of the evolution of host range and the evolutionary consequences of specialization.     

Multiple data sets,congruence, and hypothesis testing for the phylogeny of basal groups of the lizard genus Sceloporus (Squamata,Phrynosomatidae)

Several data partitions, including nuclear and mitochondrial gene sequences, chromosomes, isoenzymes, and morphological characters, were used to propose a new phylogeny and to test previously published hypotheses about the phylogenetic positions of basal clades of the lizard genus Sceloporus and the relationship of Sceloporus to the former genus "Sator". In accord with earlier studies, our results grouped "Sator" as internal to Sceloporus, and both support a hypothesis of transgulfian vicariance for the origin of the former genus "Sator" on islands in the Sea of Cortez. Robustness of support for internal nodes in our best tree was established though widely used indices (bootstrap proportions, decay values) but also through congruence among independent data partitions. Several deep nodes in the tree recovered by several methods, including equally weighted and differentially weighted parsimony and maximum likelihood models, are only weakly supported by the traditional indices. This methodological concordance is taken as evidence for insensitivity of the deep structure of the topology to alternative assumptions.     

Advantages and Disadvantages of Molecular Phylogenetics: A Case Study of Ascaridoid Nematodes

The advantages of nucleotide sequence data for studying phylogeny have been shown to include number of potential characters available for comparison, rate independence between molecular and morphological evolution, and utility of molecular data for modeling patterns of nucleotide substitution. Potential pitfalls have also been revealed and include difficulties of inferring positional homology, incongruence between organismal and gene genealogies, and low likelihood of recovering the correct phylogeny given certain patterns in the timing of speciation events. Statistical methods for comparing phylogenetic hypotheses have been used to assess the reliability of alternative trees for ascaridoid nematodes. Based on partial ribosomal RNA sequences, tree topologies inconsistent with monophyly of the Ascaridinae were significantly worse by maximum likelihood inference. The topology of the maximum parsimony tree based on full-length sequences of 18S rRNA and 300 nucleotides of Cytochrome oxidase II for 13 ascaridoid species was generally consistent with traditional taxonomic expectations at lower ranks, but inconsistent with most proposed arrangements at higher taxonomic levels.     

The Phylogeny of the Families Lecanoraceae and Bacidiaceae (Lichenized Ascomycota) Inferred from Nuclear SSU rDNA Sequences

Abstract: The phylogeny of the families Lecanoraceae and Bacidiaceae (Lecanorales, Ascomycota) was investigated using 29 nuclear small subunit ribosomal DNA sequences, 9 of which were newly determined. The data set contained 368 variable characters, 234 of which were parsimony-informative. Phylogenetic estimations were performed with maximum parsimony and maximum likelihood optimality criteria. In the most parsimonious and most likely reconstructions, the Bacidiaceae sensu Hafellner 1988 forms a monophyletic group and the Lecanoraceae sensu Hafellner a paraphyletic group. The genera Tephromela and Scoliciosporum appear to belong outside these families. However, the hypothesis that the Lecanoraceae sensu Hafellner is monophyletic cannot be rejected, as indicated by a Kishino-Hasegawa test. Three hypotheses were rejected by Kishino-Hasegawa tests, viz. (1) that the Lecanoraceae and Bacidiaceae together form a monophyletic group; (2) that both the Lecanoraceae (incl. Scoliciosporum ) and Bacidiaceae (incl. Tephromela ) are monophyletic; and (3) that the ascus apex anatomy reflects phylogeny. The suborder Lecanorineae is paraphyletic unless the Stereocaulaceae and Cladoniaceae are included. One or both of the Bacidia and Lecanora types of ascus have probably evolved at least twice.     

Inferring hominoid and early hominid phylogeny using craniodental characters: the role of fossil taxa

Recent discoveries of new fossil hominid species have been accompanied by several phylogenetic hypotheses. All of these hypotheses are based on a consideration of hominid craniodental morphology. However, Collard and Wood (2000) suggested that cladograms derived from craniodental data are inconsistent with the prevailing hypothesis of ape phylogeny based on molecular data. The implication of their study is that craniodental characters are unreliable indicators of phylogeny in hominoids and fossil hominids but, notably, their analysis did not include extinct species. We report here on a cladistic analysis designed to test whether the inclusion of fossil taxa affects the ability of morphological characters to recover the molecular ape phylogeny. In the process of doing so, the study tests both Collard and Wood's (2000) hypothesis of character reliability, and the several recently proposed hypotheses of early hominid phylogeny. One hundred and ninety-eight craniodental characters were examined, including 109 traits that traditionally have been of interest in prior studies of hominoid and early hominid phylogeny, and 89 craniometric traits that represent size-corrected linear dimensions measured between standard cranial landmarks. The characters were partitioned into two data sets. One set contained all of the characters, and the other omitted the craniometric characters. Six parsimony analyses were performed; each data set was analyzed three times, once using an ingroup that consisted only of extant hominoids, a second time using an ingroup of extant hominoids and extinct early hominids, and a third time excluding Kenyanthropus platyops. Results suggest that the inclusion of fossil taxa can play a significant role in phylogenetic analysis. Analyses that examined only extant taxa produced most parsimonious cladograms that were inconsistent with the ape molecular tree. In contrast, analyses that included fossil hominids were consistent with that tree. This consistency refutes the basis for the hypothesis that craniodental characters are unreliable for reconstructing phylogenetic relationships. Regarding early hominids, the relationships of Sahelanthropus tchadensis and Ardipithecus ramidus were relatively unstable. However, there is tentative support for the hypotheses that S. tchadensis is the sister taxon of all other hominids. There is support for the hypothesis that A. anamensis is the sister taxon of all hominids except S. tchadensis and Ar. ramidus. There is no compelling support for the hypothesis that Kenyanthropus platyops shares especially close affinities with Homo rudolfensis. Rather, K. platyops is nested within the Homo + Paranthropus + Australopithecus africanus clade. If K. platyops is a valid species, these relationships suggest that Homo and Paranthropus are likely to have diverged from other hominids much earlier than previously supposed. There is no support for the hypothesis that A. garhi is either the sister taxon or direct ancestor of the genus Homo. Phylogenetic relationships indicate that Australopithecus is paraphyletic. Thus, A. anamensis and A. garhi should be allocated to new genera.     

Molecular phylogeny of the Drosophila tripunctata and closely related species groups (Diptera: Drosophilidae)

We suggest a new phylogenetic hypothesis for the tripunctata radiation based on sequences of mitochondrial genes. Phylogenetic trees were reconstructed by parsimony, maximum likelihood and Bayesian methods. We performed tests for hypotheses of monophyly for taxonomic groups and other specific hypotheses. Results reject the monophyly for the tripunctata group whereas monophyly is not rejected for the tripunctata radiation and other specific groups within the radiation. Although most of the basal nodes were unresolved we were able to identify four clusters within the tripunctata radiation. These results suggest the collection of additional data before a proper taxonomic revision could be proposed.     

Molecular phylogeny of the hominoids: inferences from multiple independent DNA sequence data sets

Consensus on the evolutionary relationships of humans, chimpanzees, and gorillas has not been reached, despite the existence of a number of DNA sequence data sets relating to the phylogeny, partly because not all gene trees from these data sets agree. However, given the well-known phenomenon of gene tree-species tree mismatch, agreement among gene trees is not expected. A majority of gene trees from available DNA sequence data support one hypothesis, but is this evidence sufficient for statistical confidence in the majority hypothesis? All available DNA sequence data sets showing phylogenetic resolution among the hominoids are grouped according to genetic linkage of their corresponding genes to form independent data sets. Of the 14 independent data sets defined in this way, 11 support a human- chimpanzee clade, 2 support a chimpanzee-gorilla clade, and one supports a human-gorilla clade. The hypothesis of a trichotomous speciation event leading to Homo; Pan, and Gorilla can be firmly rejected on the basis of this data set distribution. The multiple-locus test (Wu 1991), which evaluates hypotheses using gene tree-species tree mismatch probabilities in a likelihood ratio test, favors the phylogeny with a Homo-Pan clade and rejects the other alternatives with a P value of 0.002. When the probabilities are modified to reflect effective population size differences among different types of genetic loci, the observed data set distribution is even more likely under the Homo-Pan clade hypothesis. Maximum-likelihood estimates for the time between successive hominoid divergences are in the range of 300,000-2,800,000 years, based on a reasonable range of estimates for long-term hominoid effective population size and for generation time. The implication of the multiple-locus test is that existing DNA sequence data sets provide overwhelming and sufficient support for a human-chimpanzee clade: no additional DNA data sets need to be generated for the purpose of estimating hominoid phylogeny. Because DNA hybridization evidence (Caccone and Powell 1989) also supports a Homo-Pan clade, the problem of hominoid phylogeny can be confidently considered solved.      

The MinMax Squeeze: guaranteeing a minimal tree for population data

We report that for population data, where sequences are very similar to one another, it is often possible to use a two-pronged (MinMax Squeeze) approach to prove that a tree is the shortest possible under the parsimony criterion. Such population data can be in a range where parsimony is a maximum likelihood estimator. This is in sharp contrast to the case with species data, where sequences are much further apart and the problem of guaranteeing an optimal phylogenetic tree is known to be computationally prohibitive for realistic numbers of species, irrespective of whether likelihood or parsimony is the optimality criterion. The Squeeze uses both an upper bound (the length of the shortest tree known) and a lower bound derived from partitions of the columns (the length of the shortest tree possible). If the two bounds meet, the shortest known tree is thus proven to be a shortest possible tree. The implementation is first tested on simulated data sets and then applied to 53 complete human mitochondrial genomes. The shortest possible trees for those data have several significant improvements from the published tree. Namely, a pair of Australian lineages comes deeper in the tree (in agreement with archaeological data), and the non-African part of the tree shows greater agreement with the geographical distribution of lineages.     

Basal actinopterygian relationships: a mitogenomic perspective on the phylogeny of the "ancient fish"

The basal actinopterygians comprise four major lineages (polypteriforms, acipenseriforms, lepisosteids, and Amia) and have been collectively called "ancient fish." We investigated the phylogeny of this group of fishes in relation to teleosts using mitochondrial genomic (mitogenomic) data, and compared this to the various alternative phylogenetic hypotheses that have been proposed previously. In addition to the previously determined complete mitochondrial DNA (mtDNA) sequences from 14 teleosts and two outgroups, we used newly determined mitogenomic sequences of 12 purposefully chosen species representing all the ancient fish lineages plus related teleosts. This data set comprised concatenated nucleotide sequences from 12 protein-coding genes (excluding the ND6 gene and third codon positions) and 22 transfer RNA (tRNA) genes (stem regions only) and these data were subjected to maximum parsimony, maximum likelihood, and Bayesian analyses. The resultant trees from the three methods were well resolved and largely congruent, with most internal branches being supported by high statistical values. Mitogenomic data strongly supported not only the monophyly of the teleosts (osteoglossomorphs and above), but also a sister-group relationship between the teleosts and a clade comprising the acipenseriforms, lepisosteids, and Amia, with the polypteriforms occupying the most basal position in the actinopterygian phylogeny. Although the tree topology differed from any of the previously proposed hypotheses based on morphology, it exhibited congruence with a recently proposed novel hypothesis based on nuclear markers.     

The Role of Functional Analysis in Phylogenetic Inference: Examples from the History of the Xiphosura

SYNOPSIS. Conventional cladistic analyses of phylogeny can beinterpreted as operating at the level of phylogenetic trees.They assume that all "evolutionary steps" (transitions fromone character state to the next, along a morphocline) are independentand equal, and, on that basis, select the cladogram which isconsistent with the most parsimonious trees. Evaluation of theassumptions of independence and equality requires considerationof hypotheses at the levelof scenarios. In some cases, argumentsbased on functional analysis can suggest revised interpretationsof either homology or polarity. If properly formulated, thesearguments can alter the evaluation of parsimony for trees tothe extent that even the choice of cladogram is affected. Thestructure of scenario level arguments is identical to that ofarguments operating at tree level. Examples of phylogeneticinference in the context of xiphosurans (horseshoe crabs), usingboth comparative morphological and functional analysis, illustratethis approach. In different cases, orthodox interpretationsof relationship are either challenged or corroborated. Althoughthe introduction of functional analysis into the process ofphylogenetic inference may appear to compromise the usefulnessof the reconstructed phylogeny for testing hypotheses concerningthe role of natural selection in evolution, it actually increasesthe strength of such tests.     

Understanding mammalian evolution using Bayesian phylogenetic inference

1. Phylogenetic trees are critical in addressing evolutionary hypotheses; however, the reconstruction of a phylogeny is no easy task. This process has recently been made less arduous by using a Bayesian statistical approach. This method offers the advantage that one can determine the probability of some hypothesis (i.e. a phylogeny), conditional on the observed data (i.e. nucleotide sequences). 2. By reconstructing phylogenies using Bayes’ theorem in combination with Markov chain Monte Carlo, phylogeneticists are able to test hypotheses more quickly compared with using standard methods such as neighbour-joining, maximum likelihood or parsimony. Critics of the Bayesian approach suggest that it is not a panacea, and argue that the prior probability is too subjective and the resulting posterior probability is too liberal compared with maximum likelihood. 3. These issues are currently debated in the arena of mammalian evolution. Recently, proponents and opponents of the Bayesian approach have constructed the mammalian phylogeny using different methods under different conditions and with a variety of parameters. These analyses showed the robustness (or lack of) of the Bayesian approach. In the end, consensus suggests that Bayesian methods are robust and give essentially the same answer as maximum likelihood methods but in less time. 4. Approaches based on fossils and molecules typically agree on ordinal-level relationships among mammals but not on higher-level relationships, as Bayesian analyses recognize the African radiation, Afrotheria, and the two Laurasian radiations, Laurasiatheria and Euarchontoglires, whereas fossils did not predict Afrotheria.     

On the best evolutionary rate for phylogenetic analysis

The effect of the evolutionary rate of a gene on the accuracy of phylogeny reconstruction was examined by computer stimulation. The evolutionary rate is measured by the tree length, that is, the expected total number of nucleotide substitutions per site on the phylogeny. DNA sequence data were simulated using both fixed trees with specified branch lengths and random trees with branch lengths generated from a model of cladogenesis. The parsimony and likelihood methods were used for phylogeny reconstruction, and the proportion of correctly recovered branch partitions by each method was estimated. Phylogenetic methods including parsimony appear quite tolerant of multiple substitutions at the same site. The optimum levels of sequence divergence were even higher than upper limits previously suggested for saturation of substitutions, indicating that the problem of saturation may have been exaggerated. Instead, the lack of information at low levels of divergence should be seriously considered in evaluation of a gene's phylogenetic utility, especially when the gene sequence is short. The performance of parsimony, relative to that of likelihood, does not necessarily decrease with the increase of the evolutionary rate.     

Mitochondrial DNA,morphology, and the phylogenetic relationships of Antarctic icefishes (Notothenioidei: Channichthyidae)

The Channichthyidae is a lineage of 16 species in the Notothenioidei, a clade of fishes that dominate Antarctic near-shore marine ecosystems with respect to both diversity and biomass. Among four published studies investigating channichthyid phylogeny, no two have produced the same tree topology, and no published study has investigated the degree of phylogenetic incongruence between existing molecular and morphological datasets. In this investigation we present an analysis of channichthyid phylogeny using complete gene sequences from two mitochondrial genes (ND2 and 16S) sampled from all recognized species in the clade. In addition, we have scored all 58 unique morphological characters used in three previous analyses of channichthyid phylogenetic relationships. Data partitions were analyzed separately to assess the amount of phylogenetic resolution provided by each dataset, and phylogenetic incongruence among data partitions was investigated using incongruence length difference (ILD) tests. We utilized a parsimony-based version of the Shimodaira-Hasegawa test to determine if alternative tree topologies are significantly different from trees resulting from maximum parsimony analysis of the combined partition dataset. Our results demonstrate that the greatest phylogenetic resolution is achieved when all molecular and morphological data partitions are combined into a single maximum parsimony analysis. Also, marginal to insignificant incongruence was detected among data partitions using the ILD. Maximum parsimony analysis of all data partitions combined results in a single tree, and is a unique hypothesis of phylogenetic relationships in the Channichthyidae. In particular, this hypothesis resolves the phylogenetic relationships of at least two species (Channichthys rhinoceratus and Chaenocephalus aceratus), for which there was no consensus among the previous phylogenetic hypotheses. The combined data partition dataset provides substantial statistical power to discriminate among alternative hypotheses of channichthyid relationships. These findings suggest the optimal strategy for investigating the phylogenetic relationships of channichthyids is one that uses all available phylogenetic data in analyses of combined data partitions.     

AFLPs resolve phylogeny and reveal mitochondrial introgression within a species flock of African electric fish (Mormyroidea: Teleostei)

Estimating species phylogeny from a single gene tree can be especially problematic for studies of species flocks in which diversification has been rapid. Here we compare a phylogenetic hypothesis derived from cytochrome b (cyt b) sequences with another based on amplified fragment length polymorphisms (AFLP) for 60 specimens of a monophyletic riverine species flock of mormyrid electric fishes collected in Gabon, west-central Africa. We analyze the aligned cyt b sequences by Wagner parsimony and AFLP data generated from 10 primer combinations using neighbor-joining from a Nei-Li distance matrix, Wagner parsimony, and Dollo parsimony. The different analysis methods yield AFLP tree topologies with few conflicting nodes. Recovered basal relationships in the group are similar between cyt b and AFLP analyses, but differ substantially at many of the more derived nodes. More of the clades recovered with the AFLP characters are consistent with the morphological characters used to designate operational taxonomic units in this group. These results support our hypothesis that the mitochondrial gene tree differs from the overall species phylogeny due at least in part to mitochondrial introgession among lineages. Mapping the two forms of electric organ found in this group onto the AFLP tree suggests that posteriorly innervated electrocytes with nonpenetrating stalks have independently evolved from anteriorly innervated, penetrating-stalk electrocytes at least three times.