Cytochrome b and Bayesian inference of whale phylogeny

In the mid 1990s cytochrome b and other mitochondrial DNA data reinvigorated cetacean phylogenetics by proposing many novel and provocative hypotheses of cetacean relationships. These results sparked a revision and reanalysis of morphological datasets, and the collection of new nuclear DNA data from numerous loci. Some of the most controversial mitochondrial hypotheses have now become benchmark clades, corroborated with nuclear DNA and morphological data; others have been resolved in favor of more traditional views. That major conflicts in cetacean phylogeny are disappearing is encouraging. However, most recent papers aim specifically to resolve higher-level conflicts by adding characters, at the cost of densely sampling taxa to resolve lower-level relationships. No molecular study to date has included more than 33 cetaceans. More detailed molecular phylogenies will provide better tools for evolutionary studies. Until more genes are available for a high number of taxa, can we rely on readily available single gene mitochondrial data? Here, we estimate the phylogeny of 66 cetacean taxa and 24 outgroups based on Cytb sequences. We judge the reliability of our phylogeny based on the recovery of several deep-level benchmark clades. A Bayesian phylogenetic analysis recovered all benchmark clades and for the first time supported Odontoceti monophyly based exclusively on analysis of a single mitochondrial gene. The results recover the monophyly of all but one family level taxa within Cetacea, and most recently proposed super- and subfamilies. In contrast, parsimony never recovered all benchmark clades and was sensitive to a priori weighting decisions. These results provide the most detailed phylogeny of Cetacea to date and highlight the utility of both Bayesian methodology in general, and of Cytb in cetacean phylogenetics. They furthermore suggest that dense taxon sampling, like dense character sampling, can overcome problems in phylogenetic reconstruction.     

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.     

c-mos variation in songbirds: molecular evolution, phylogenetic implications, and comparisons with mitochondrial differentiation

Nucleotide sequences from the c-mos proto-oncogene have previously been used to reconstruct the phylogenetic relationships between distantly related vertebrate taxa. To explore c-mos variation at shallower levels of avian divergence, we compared c-mos sequences from representative passerine taxa that span a range of evolutionary differentiation, from basal passerine lineages to closely allied genera. Phylogenetic reconstructions based on these c-mos sequences recovered topologies congruent with previous DNA-DNA hybridization-based reconstructions, with many nodes receiving high support, as indicated by bootstrap and reliability values. One exception was the relationship of Acanthisitta to the remaining passerines, where the c-mos-based searches indicated a three-way polytomy involving the Acanthisitta lineage and the suboscine and oscine passerine clades. We also compared levels of c-mos and mitochondrial differentiation across eight oscine passerine taxa and found that c-mos nucleotide substitutions accumulate at a rate similar to that of transversion substitutions in mitochondrial protein-coding genes. These comparisons suggest that nuclear-encoded loci such as c-mos provide a temporal window of phylogenetic resolution that overlaps the temporal range where mitochondrial protein-coding sequences have their greatest utility and that c-mos substitutions and mtDNA transversions can serve as complementary, informative, and independent phylogenetic markers for the study of avian relationships.     

The phylogeny of acorn weevils (genus Curculio) from mitochondrial and nuclear DNA sequences: the problem of incomplete data

We considered the contribution of two mitochondrial and two nuclear data sets for the phylogenetic reconstruction of 22 species of seed beetles in the genus Curculio (Coleoptera: Cuculionidae). A phylogenetic tree from representatives found on various hosts was inferred from a combined data set of mitochondrial DNA cytochrome oxidase subunit I, mitochondrial cytochrome b, nuclear elongation factor 1alpha, and nuclear phosphoglycerate mutase, used for the first time as a molecular marker. Separate parsimony analyses of each data set showed that individual gene trees were mainly congruent and often complementary in the support of clades but the analysis was complicated by failure of PCR amplification of nuclear genes for many taxa and hence missing data entries. When the four gene partitions were combined in a simultaneous analysis despite the missing data, this increased the resolution and taxonomic coverage compared to the individual source trees. Alternative approaches of combining the information via supertree methodology produced a comparatively less resolved tree, and hence seem inferior to combining data matrices even in cases where numerous taxa are missing. The molecular data suggest a classification of the European species into two species groups that are in accordance with morphological characteristics but the data do no support any of the previously recognised American species groups.     

Phylogenetic relationships and divergence times among mustelids (Mammalia: Carnivora) based on nucleotide sequences of the nuclear interphotoreceptor retinoid binding protein and mitochondrial cytochrome b genes

Phylogenetic relationships among 20 species-group taxa of Mustelidae, representing Mustelinae (Mustela, Martes, Gulo), Lutrinae (Enhydra), and Melinae (Meles), were examined using nucleotide sequences of the nuclear interphotoreceptor retinoid binding protein (IRBP) and mitochondrial cytochrome b genes. Neighbor-joining and maximum-parsimony phylogenetic analyses on these genes separately and combined were conducted. While IRBP performed better than cytochrome b in recovering more-inclusive clades, cytochrome b demonstrated more resolving power in recovering less-inclusive clades. Strong support was found for a close affinity of Enhydra with Mustela to the exclusion of Martes and Gulo (causing Mustelinae to be paraphyletic); the most-basal position of Mustela vison within Mustela, followed by Mustela erminea; an association of Mustela lutreola, Mustela itatsi, Mustela sibirica, and the subgenus Putorius (including Mustela putorius and Mustela eversmanii), to the exclusion of Mustela nivalis and Mustela altaica; and a basal position of Mustela itatsi to a clade containing Mustela sibirica and Putorius. Whereas cytochrome b strongly supported Mustela lutreola as the sister species to Putorius, IRBP strongly supported its basal placement to the Mustela itatsi-Mustela sibirica-Putorius clade. The low level of sequence divergence in cytochrome b between Mustela lutreola and Putorius is therefore a result of interspecific mitochondrial introgression between these taxa, rather than a recent origin of Mustela lutreola in a close relationship to Putorius. Time estimates inferred from IRBP and cytochrome b for mustelid divergence events are mostly in agreement with the fossil record.     

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.     

Use of mitogenomic information in teleostean molecular phylogenetics: a tree-based exploration under the maximum-parsimony optimality criterion

We explored the phylogenetic utility and limits of the individual and concatenated mitochondrial genes for reconstructing the higher-level relationships of teleosts, using the complete (or nearly complete) mitochondrial DNA sequences of eight teleosts (including three newly determined sequences), whose relative phylogenetic positions were noncontroversial. Maximum-parsimony analyses of the nucleotide and amino acid sequences of 13 protein-coding genes from the above eight teleosts, plus two outgroups (bichir and shark), indicated that all of the individual protein-coding genes, with the exception of ND5, failed to recover the expected phylogeny, although unambiguously aligned sequences from 22 concatenated transfer RNA (tRNA) genes (stem regions only) recovered the expected phylogeny successfully with moderate statistical support. The phylogenetic performance of the 13 protein-coding genes in recovering the expected phylogeny was roughly classified into five groups, viz. very good (ND5, ND4, COIII, COI), good (COII, cyt b), medium (ND3, ND2), poor (ND1, ATPase 6), and very poor (ND4L, ND6, ATPase 8). Although the universality of this observation was unclear, analysis of successive concatenation of the 13 protein-coding genes in the same ranking order revealed that the combined data sets comprising nucleotide sequences from the several top-ranked protein-coding genes (no 3rd codon positions) plus the 22 concatenated tRNA genes (stem regions only) best recovered the expected phylogeny, with all internal branches being supported by bootstrap values >90%. We conclude that judicious choice of mitochondrial genes and appropriate data weighting, in conjunction with purposeful taxonomic sampling, are prerequisites for resolving higher-level relationships in teleosts under the maximum-parsimony optimality criterion.     

Phylogeny and comparative substitution rates of frogs inferred from sequences of three nuclear genes

Phylogenetic relationships among major clades of anuran amphibians were studied using partial sequences of three nuclear protein coding genes, Rag-1, Rag-2, and rhodopsin in 26 frog species from 18 families. The concatenated nuclear data set comprised 2,616 nucleotides and was complemented by sequences of the mitochondrial 12S and 16S rRNA genes for analyses of evolutionary rates. Separate and combined analyses of the nuclear markers supported the monophyly of modern frogs (Neobatrachia), whereas they did not provide support for the monophyly of archaic frog lineages (Archaeobatrachia), contrary to previous studies based on mitochondrial data. The Neobatrachia contain two well supported clades that correspond to the subfamilies Ranoidea (Hyperoliidae, Mantellidae, Microhylidae, Ranidae, and Rhacophoridae) and Hyloidea (Bufonidae, Hylidae, Leptodactylidae, and Pseudidae). Two other families (Heleophrynidae and Sooglossidae) occupied basal positions and probably represent ancient relicts within the Neobatrachia, which had been less clearly indicated by previous mitochondrial analyses. Branch lengths of archaeobatrachians were consistently shorter in all separate analyses, and nonparametric rate smoothing indicated accelerated substitution rates in neobatrachians. However, relative rate tests confirmed this tendency only for mitochondrial genes. In contrast, nuclear gene sequences from our study and from an additional GenBank survey showed no clear phylogenetic trends in terms of differences in rates of molecular evolution. Maximum likelihood trees based on Rag-1 and using only one neobatrachian and one archaeobatrachian sequence, respectively, even had longer archaeobatrachian branches averaged over all pairwise comparisons. More data are necessary to understand the significance of a possibly general assignation of short branches to basal and species-poor taxa by tree-reconstruction algorithms.     

Phylogeny of Trichoptera (caddisflies): characterization of signal and noise within multiple datasets

Trichoptera are holometabolous insects with aquatic larvae that, together with the Lepidoptera, make up the Amphiesmenoptera. Despite extensive previous morphological work, little phylogenetic agreement has been reached about the relationship among the three suborders--Annulipalpia, Spicipalpia, and Integripalpia--or about the monophyly of Spicipalpia. In an effort to resolve this conflict, we sequenced fragments of the large and small subunit nuclear ribosomal RNAs (1078 nt; D1, D3, V4-5), the nuclear elongation factor 1 alpha gene (EF-1 alpha; 1098 nt), and a fragment of mitochondrial cytochrome oxidase I (COI; 411 nt). Seventy adult and larval morphological characters were reanalyzed and added to molecular data in a combined analysis. We evaluated signal and homoplasy in each of the molecular datasets and attempted to rank the particular datasets according to how appropriate they were for inferring relationships among suborders. This evaluation included testing for conflict among datasets, comparing tree lengths among alternative hypotheses, measuring the left-skew of tree-length distributions from maximally divergent sets of taxa, evaluating the recovery of expected clades, visualizing whether or not substitutions were accumulating with time, and estimating nucleotide compositional bias. Although all these measures cast doubt on the reliability of the deep-level signal coming from the nucleotides of the COI and EF-1 alpha genes, these data could still be included in combined analyses without overturning the results from the most conservative marker, the rRNA. The different datasets were found to be evolving under extremely different rates. A site-specific likelihood method for dealing with combined data with nonoverlapping parameters was proposed, and a similar weighting scheme under parsimony was evaluated. Among our phylogenetic conclusions, we found Annulipalpia to be the most basal of the three suborders, with Spicipalpia and Integripalpia forming a clade. Monophyly of Annulipalpia and Integripalpia was confirmed, but the relationships among spicipalpians remain equivocal.     

Neotomine-peromyscine rodent systematics based on combined analyses of nuclear and mitochondrial DNA sequences

Recently, sequences from two nuclear genes (exon 6 of the dentin matrix protein 1 gene and intron 7 of the beta-fibrinogen gene) and one mitochondrial gene (cytochrome b gene) were used independently in an attempt to resolve phylogenetic relationships within the neotomine-peromyscine complex. Although these studies provided testable hypotheses regarding this group of rodents, the affinities of certain tribes and genera remain uncertain. To elucidate these relationships, the three data partitions were tested for heterogeneity and then concatenated according to conditional data combination and total evidence approaches. Support was found for five clades, four of which correspond to well recognized tribes (the Neotomini, Peromyscini=Reithrodontomyini, Baiomyini, and Tylomyini). Recommendations are made regarding the recognition of Ochrotomys as a tribe of its own, the Ochrotomyini, paralleling other recent findings. The Peromyscini, Baiomyini, and Ochrotomyini are unresolved in relation to each other, but as a whole are sister to the Neotomini. The Tylomyini is basal to all clades. It appears that combined data from the nuclear and mitochondrial genes (analyzing all three partitions simultaneously) resulted in the best phylogenetic hypothesis regarding the complex.     

Multiple nuclear insertions of mitochondrial cytochrome b sequences in callitrichine primates

We report the presence of four nuclear paralogs of a 380-bp segment of cytochrome b in callitrichine primates (marmosets and tamarins). The mitochondrial cytochrome b sequence and each nuclear paralog were obtained from several species, allowing multiple comparisons of rates and patterns of substitution both between mitochondrial and nuclear sequences and among nuclear sequences. The mitochondrial DNA had high overall rates of molecular evolution and a strong bias toward substitutions at third codon positions. Rates of molecular evolution among the nuclear sequences were low and constant, and there were small differences in substitution patterns among the nuclear clades which were probably attributable to the small number of sites involved. A novel method of phylogenetic reconstruction based on the large difference in rates of evolution at different codon positions among mitochondrial and nuclear clades was used to determine whether different nuclear paralogs represent independent transposition events or duplications following a single insertion. This method is generally applicable in cases where differences in pattern of molecular evolution are known, and it showed that at least three of the four nuclear clades represent independent transposition events. The insertion events giving rise to two of the nuclear clades predate the divergence of the callitrichines, whereas those leading to the other two nuclear clades may have occurred in the common ancestor of marmosets.     

Lizards from the end of the world: phylogenetic relationships of the Liolaemus lineomaculatus section (Squamata: Iguania: Liolaemini)

The Liolaemus lineomaculatus section is a geographically widely distributed group of lizards from the Patagonian region of southern South America, and includes 18 described species representing the most southerly distributed Liolaemus taxa (the genus includes 228 species and extends from Tierra del Fuego north to south-central Peru). Despite high species diversity, the phylogenetic relationships of this section are unknown. In the present work we sampled all described species in the L. lineomaculatus section as well as currently undescribed candidate species to reconstruct the first complete phylogenetic hypothesis for the clade. Our data set included four anonymous nuclear loci, three nuclear protein-coding loci, and two mitochondrial genes. We compared results obtained with three different phylogenetic methods for the concatenated data set (Maximum Parsimony, Maximum Likelihood and Bayesian Inference) with a coalescent-based species tree approach (BEST), and recovered congruent, strongly-supported topological arrangements across all methods. We identified four main clades within the L. lineomaculatus section: the lineomaculatus, magellanicus, somuncurae, and kingii+archeforus groups, for which we estimated divergence times. We discuss the taxonomic implications of these results and how the future integration of phylogeographic, niche modeling and morphological approaches will allow testing biogeographical hypotheses in this clade.     

STBase: One Million Species Trees for Comparative Biology

Comprehensively sampled phylogenetic trees provide the most compelling foundations for strong inferences in comparative evolutionary biology. Mismatches are common, however, between the taxa for which comparative data are available and the taxa sampled by published phylogenetic analyses. Moreover, many published phylogenies are gene trees, which cannot always be adapted immediately for species level comparisons because of discordance, gene duplication, and other confounding biological processes. A new database, STBase, lets comparative biologists quickly retrieve species level phylogenetic hypotheses in response to a query list of species names. The database consists of 1 million single- and multi-locus data sets, each with a confidence set of 1000 putative species trees, computed from GenBank sequence data for 413,000 eukaryotic taxa. Two bodies of theoretical work are leveraged to aid in the assembly of multi-locus concatenated data sets for species tree construction. First, multiply labeled gene trees are pruned to conflict-free singly-labeled species-level trees that can be combined between loci. Second, impacts of missing data in multi-locus data sets are ameliorated by assembling only decisive data sets. Data sets overlapping with the user’s query are ranked using a scheme that depends on user-provided weights for tree quality and for taxonomic overlap of the tree with the query. Retrieval times are independent of the size of the database, typically a few seconds. Tree quality is assessed by a real-time evaluation of bootstrap support on just the overlapping subtree. Associated sequence alignments, tree files and metadata can be downloaded for subsequent analysis. STBase provides a tool for comparative biologists interested in exploiting the most relevant sequence data available for the taxa of interest. It may also serve as a prototype for future species tree oriented databases and as a resource for assembly of larger species phylogenies from precomputed trees.     

Toward the resolution of an explosive radiation--a multilocus phylogeny of oceanic dolphins (Delphinidae)

Oceanic dolphins (Delphinidae) are the product of a rapid radiation that yielded ～36 extant species of small to medium-sized cetaceans that first emerged in the Late Miocene. Although they are a charismatic group of organisms that have become poster children for marine conservation, many phylogenetic relationships within Delphinidae remain elusive due to the slow molecular evolution of the group and the difficulty of resolving short branches from successive cladogenic events. Here I combine existing and newly generated sequences from four mitochondrial (mt) genes and 20 nuclear (nu) genes to reconstruct a well-supported phylogenetic hypothesis for Delphinidae. This study compares maximum-likelihood and Bayesian inference methods of several data sets including mtDNA, combined nuDNA, gene trees of individual nuDNA loci, and concatenated mtDNA+nuDNA. In addition, I contrast these standard phylogenetic analyses with the species tree reconstruction method of Bayesian concordance analysis (BCA). Despite finding discordance between mtDNA and individual nuDNA loci, the concatenated matrix recovers a completely resolved and robustly supported phylogeny that is also broadly congruent with BCA trees. This study strongly supports groupings such as Delphininae, Lissodelphininae, Globicephalinae, Sotalia+Delphininae, Steno+Orcaella+Globicephalinae, and Leucopleurus acutus, Lagenorhynchus albirostris, and Orcinus orca as basal delphinid taxa.     

Resolving an ancient, rapid radiation in Saxifragales

Despite the prior use of approximately 9000 bp, deep-level relationships within the angiosperm clade, Saxifragales remain enigmatic, due to an ancient, rapid radiation (89.5 to 110 Ma based on the fossil record). To resolve these deep relationships, we constructed several new data sets: (1) 16 genes representing the three genomic compartments within plant cells (2 nuclear, 10 plastid, 4 mitochondrial; aligned, analyzed length = 21,460 bp) for 28 taxa; (2) the entire plastid inverted repeat (IR; 26,625 bp) for 17 taxa; (3) "total evidence" (50,845 bp) for both 17 and 28 taxa (the latter missing the IR). Bayesian and ML methods yielded identical topologies across partitions with most clades receiving high posterior probability (pp = 1.0) and bootstrap (95% to 100%) values, suggesting that with sufficient data, rapid radiations can be resolved. In contrast, parsimony analyses of different partitions yielded conflicting topologies, particularly with respect to the placement of Paeoniaceae, a clade characterized by a long branch. In agreement with published simulations, the addition of characters increased bootstrap support for the putatively erroneous placement of Paeoniaceae. Although having far fewer parsimony-informative sites, slowly evolving plastid genes provided higher resolution and support for deep-level relationships than rapidly evolving plastid genes, yielding a topology close to the Bayesian and ML total evidence tree. The plastid IR region may be an ideal source of slowly evolving genes for resolution of deep-level angiosperm divergences that date to 90 My or more. Rapidly evolving genes provided support for tip relationships not recovered with slowly evolving genes, indicating some complementarity. Age estimates using penalized likelihood with and without age constraints for the 28-taxon, total evidence data set are comparable to fossil dates, whereas estimates based on the 17-taxon data are much older than implied by the fossil record. Hence, sufficient taxon density, and not simply numerous base pairs, is important in reliably estimating ages. Age estimates indicate that the early diversification of Saxifragales occurred rapidly, over a time span as short as 6 million years. Between 25,000 and 50,000 bp were needed to resolve this radiation with high support values. Extrapolating from Saxifragales, a similar number of base pairs may be needed to resolve the many other deep-level radiations of comparable age in angiosperms.     

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.     

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.     

Inter-generic relationships of the crows, jays, magpies and allied groups (Aves: Corvidae) based on nucleotide sequence data

Phylogenetic relationships were studied based on DNA sequences obtained from all recognized genera of the family Corvidae sensu stricto . The aligned data set consists 2589 bp obtained from one mitochondrial and two nuclear genes. Maximum parsimony, maximum-likelihood, and Bayesian inference analyses were used to estimate phylogenetic relationships. The analyses were done for each gene separately, as well as for all genes combined. An analysis of a taxonomically expanded data set of cytochrome b sequences was performed in order to infer the phylogenetic positions of six genera for which nuclear genes could not be obtained. Monophyly of the Corvidae is supported by all analyses, as well as by the occurrence of a deletion of 16 bp in the β-fibrinogen intron in all ingroup taxa. Temnurus and Pyrrhocorax are placed as the sister group to all other corvids, while Cissa and Urocissa appear as the next clade inside them. Further up in the tree, two larger and well-supported clades of genera were recovered by the analyses. One has an entirely New World distribution (the New World jays), while the other includes mostly Eurasian (and one African) taxa. Outside these two major clades are Cyanopica and Perisoreus whose phylogenetic positions could not be determined by the present data. A biogeographic analysis of our data suggests that the Corvidae underwent an initial radiation in Southeast Asia. This is consistent with the observation that almost all basal clades in the phylogenetic tree consist of species adapted to tropical and subtropical forest habitats.     

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

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