Molecular approach to the phylogenetics of sea spiders (Arthropoda: Pycnogonida) using partial sequences of nuclear ribosomal DNA

The phylogenetic relationships among major evolutionary lineages of the sea spiders (subphylum Pycnogonida) were investigated using partial sequences of nuclear DNA, 18S, and 28S ribosomal genes. Topological differences were obtained with separate analyses of 18S and 28S, and estimates of phylogeny were found to be significantly different between a combined molecular data set (18S and 28S) and a subset of a morphological data matrix analyzed elsewhere. Colossendeidae played a major role in the conflicts; it was closely related to Callipallenidae or Nymphonidae with 18S or 28S, respectively, but related to Ammotheidae according to morphological characters. Austrodecidae was defined as a basal taxon for Pycnogonida by these molecular data. The 18S sequences were surprisingly conserved among pycnogonid taxa, suggesting either an unusual case of slow evolution of the gene, or an unexpected recent divergence of pycnogonid lineages. Notwithstanding difficulties such as non-optimal taxon sampling, this is the first attempt to reconstruct the pycnogonid phylogeny based on DNA. Continued studies of sequences and other characters should increase the reliability of the analyses and our understanding of the phylogenetics of sea spiders.     

DNA在鸟类分子系统发育研究中的应用

鸟类分子系统发育研究中常用的DNA技术有DNA杂交、RFLP和DNA序列分析等。DNA杂交技术曾在鸟类中有过大规模的应用,并由此诞生了一套新的鸟类分类系统。在鸟类的RFLP分析中,用的最多的靶序列是线粒体DNA。DNA序列分析技术被认为是进行分子系统发育研究最有效、最可靠的方法。在DNA序列分析中,线粒体基因应用最广泛,但由于其自身的一些不足,近年来,不少学者把目光投向了核基因,将线粒体基因和核基因结合起来进行系统发育研究。目前在鸟类分子系统发育中,应用较多的核基因是scnDNA,其内含子可以用于中等阶元水平的系统研究,而外显子主要用于高等阶元的系统研究。除了分子标记自身的问题之外,鸟类分子系统发育研究中还存在着方法上的问题,包括分子标记的选择,样本数量以及数据处理等。今后鸟类分子系统发育研究应该更加注重方法的标准化。     

Taxon influence index: assessing taxon-induced incongruities in phylogenetic inference

Understanding the evolutionary history of species is at the core of molecular evolution and is done using several inference methods. The critical issue is to quantify the uncertainty of the inference. The posterior probabilities in Bayesian phylogenetic inference and the bootstrap values in frequentist approaches measure the variability of the estimates due to the sampling of sites from genes and the sampling of genes from genomes. However, they do not measure the uncertainty due to taxon sampling. Taxa that experienced molecular homoplasy, recent selection, a spur of evolution, and so forth may disrupt the inference and cause incongruences in the estimated phylogeny. We define a taxon influence index to assess the influence of each taxon on the phylogeny. We found that although most taxa have a weak influence on the phylogeny, a small fraction of influential taxa strongly alter it even in clades only loosely related to them. We conclude that highly influential taxa should be given special attention and sampling them more thoroughly can lead to more dependable phylogenies.     

Molecular phylogeny of glass sponges (Porifera, Hexactinellida): increased taxon sampling and inclusion of the mitochondrial protein-coding gene, cytochrome oxidase subunit I

Marine sponges of the class Hexactinellida (glass sponges) are among the most understudied groups of Porifera, and molecular approaches to investigating their evolution have only recently emerged. Although these first results appeared reliable as they largely corroborated morphology-based hypotheses, they were almost exclusively based on ribosomal RNA genes (rDNA) and should, therefore, be further tested with independent types of genetic data, such as protein-coding genes. To this end, we established the mitochondrial-encoded cytochrome oxidase subunit I gene (COI) as an additional marker, and conducted phylogenetic analyses on DNA- and amino-acid level, as well as a supermatrix analysis based on combined COI DNA and rDNA alignments. Furthermore, we increased taxon sampling compared to previous studies by adding seven additional species. The COI-based phylogenies were largely congruent with the rDNA-based phylogeny but suffered from poor bootstrap support for many nodes. However, addition of the COI sequences to the rDNA data set increased resolution of the overall molecular phylogeny. Thus, although obtaining COI sequences from glass sponges turned out to be quite challenging, this gene appears to be a valuable supplement to rDNA data for molecular evolutionary studies of this group. Some implications of our extended phylogeny for the evolution and systematics of Hexactinellida are discussed.     

Mitochondrial phylogeny of Anura (Amphibia): a case study of congruent phylogenetic reconstruction using amino acid and nucleotide characters

We explore whether phylogenetic analyses of the same sequence data set at the amino acid and nucleotide level are able to recover congruent topologies, as well as the advantages and limitations of both alternative approaches. As a case study, mitochondrial protein-coding genes were used to discern among competing hypotheses on the phylogenetic relationships of major anuran amphibian lineages. To properly address this phylogenetic question, the complete nucleotide sequences of the mitochondrial genomes of two archaeobatrachian species, Ascaphus truei and Pelobates cultripes, were determined anew. Bayesian and maximum likelihood phylogenetic inferences of the same sequence data set were performed based on both amino acid and nucleotide characters, with the latter analysed either as codons or as a reduced data set of first+second (P12) codon positions. In addition, likelihood-based ratio tests were performed to evaluate the support of alternative topologies. The different data sets arrived at congruent and highly supported topologies, suggesting a similar phylogenetic resolving power of the two character types provided that correctly selected sites and appropriate evolutionary models are used. The reconstructed anuran mitochondrial phylogeny supports the paraphyly of Archaeobatrachia, with Ascaphus as sister group to all the remaining anurans, and Pelobates as sister group of Neobatrachia. However, the employed tree reconstruction methods and likelihood-based ratio tests seemed to be negatively affected by the fast evolving sequences of neobatrachians, suggesting that the phylogeny of Anura here presented is not definitive, and needs further investigation using an extended taxon sampling.     

Lack of resolution in the animal phylogeny: closely spaced cladogeneses or undetected systematic errors?

A recent phylogenomic study reported that the animal phylogeny was unresolved despite the use of 50 genes. This lack of resolution was interpreted as "a positive signature of closely spaced cladogenetic events." Here, we propose that this lack of resolution is rather due to the mutual cancellation of the phylogenetic signal (historical) and the nonphylogenetic signal (due to systematic errors) that results from inadequate taxon sampling and/or model of sequence evolution. Starting with a data set of comparable size, we use 3 different strategies to reduce the nonphylogenetic signal: 1) increasing the number of species; 2) replacing a fast-evolving species by a slowly evolving one; and 3) using a better model of sequence evolution. In all cases, the phylogenetic resolution is markedly improved, in agreement with our hypothesis that the originally reported lack of resolution was artifactual.     

Probing the relationships of the branchiopod crustaceans

The Branchiopoda display extraordinary variation in body form, even within the morphologically diverse crustaceans. To fully understand the origin and evolution of these morphological reconfigurations, a robust phylogeny of the group is essential. To infer the affinities among branchiopods, we employed two approaches to taxon and gene sampling, presented new sequence data from three genes, incorporated previously published sequence data from three additional genes, and utilized comprehensive techniques of phylogeny reconstruction. The results provided support for a number of longstanding hypotheses concerning the relationships among the orders. For example, we obtained support for the Cladoceromorpha and Gymnomera, and favoured a unique arrangement of the cladoceran orders. A few affinities remain to be resolved, particularly at the base of the Phyllopoda and within the Anomopoda. However, the results suggest that increased gene sampling is recommended for future investigations of branchiopod systematics.     

Mitochondrial genomes in the Hymenoptera and their utility as phylogenetic markers

Abstract.  In this study, we assessed the ability of mitochondrial genome sequences to recover a test phylogeny of five hymenopteran taxa from which phylogenetic relationships are well accepted. Our analyses indicated that the test phylogeny was well recovered in all nucleotide Bayesian analyses when all the available holometabolan (i.e. outgroup) taxa were included, but only in Bayesian analyses excluding third codon positions when only the hymenopteran representatives and a single outgroup were included. This result suggests that taxon sampling of the outgroup might be as important as taxon sampling of the ingroup when recovering hymenopteran phylogenetic relationships using whole mitochondrial genomes. Parsimony analyses were more sensitive to both taxon sampling and the analytical model than Bayesian analyses, and analyses using the protein dataset did not recover the test phylogeny. In general, mitochondrial genomes did not resolve the position of the Hymenoptera within the Holometabola with confidence, suggesting that an increased taxon sampling, both within the Holometabola and among outgroups, is necessary.     

Phylogenetic analysis, structural evolution and functional divergence of the 12-oxo-phytodienoate acid reductase gene family in plants

Background

Whenever different data sets arrive at conflicting phylogenetic hypotheses, only testable causal explanations of sources of errors in at least one of the data sets allow us to critically choose among the conflicting hypotheses of relationships. The large (28S) and small (18S) subunit rRNAs are among the most popular markers for studies of deep phylogenies. However, some nodes supported by this data are suspected of being artifacts caused by peculiarities of the evolution of these molecules. Arthropod phylogeny is an especially controversial subject dotted with conflicting hypotheses which are dependent on data set and method of reconstruction. We assume that phylogenetic analyses based on these genes can be improved further i) by enlarging the taxon sample and ii) employing more realistic models of sequence evolution incorporating non-stationary substitution processes and iii) considering covariation and pairing of sites in rRNA-genes.

Results

We analyzed a large set of arthropod sequences, applied new tools for quality control of data prior to tree reconstruction, and increased the biological realism of substitution models. Although the split-decomposition network indicated a high noise content in the data set, our measures were able to both improve the analyses and give causal explanations for some incongruities mentioned from analyses of rRNA sequences. However, misleading effects did not completely disappear.

Conclusion

Analyses of data sets that result in ambiguous phylogenetic hypotheses demand for methods, which do not only filter stochastic noise, but likewise allow to differentiate phylogenetic signal from systematic biases. Such methods can only rely on our findings regarding the evolution of the analyzed data. Analyses on independent data sets then are crucial to test the plausibility of the results. Our approach can easily be extended to genomic data, as well, whereby layers of quality assessment are set up applicable to phylogenetic reconstructions in general.     

Molecular systematics of rhizobia based on maximum likelihood and Bayesian phylogenies inferred from rrs, atpD, recA and nifH sequences, and their use in the classification of Sesbania microsymbionts from Venezuelan wetlands

A well-resolved rhizobial species phylogeny with 51 haplotypes was inferred from a combined atpD + recA data set using Bayesian inference with best-fit, gene-specific substitution models. Relatively dense taxon sampling for the genera Rhizobium and Mesorhizobium was achieved by generating atpD and recA sequences for six type and 24 reference strains not previously available in GenBank. This phylogeny was used to classify nine nodule isolates from Sesbania exasperata, S. punicea and S. sericea plants native to seasonally flooded areas of Venezuela, and compared with a PCR-RFLP analysis of rrs plus rrl genes and large maximum likelihood rrs and nifH phylogenies. We show that rrs phylogenies are particularly sensitive to strain choice due to the high levels of sequence mosaicism found at this locus. All analyses consistently identified the Sesbania isolates as Mesorhizobium plurifarium or Rhizobium huautlense. Host range experiments on ten legume species coupled with plasmid profiling uncovered potential novel biovarieties of both species. This study demonstrates the wide geographic and environmental distribution of M. plurifarium, that R. galegae and R. huautlense are sister lineages, and the synonymy of R. gallicum, R. mongolense and R. yanglingense. Complex and diverse phylogeographic, inheritance and host-association patterns were found for the symbiotic nifH locus. The results and the analytical approaches used herein are discussed in the context of rhizobial taxonomy and molecular systematics.     

Global versus Chinese perspectives on the phylogeny of the N-fixing clade

There has been increasing interest in integrating a regional tree of life with community assembly rules in the ecological research. This raises questions regarding the impacts of taxon sampling strategies at the regional versus global scales on the topology. To address this concern, we constructed two trees for the nitrogen-fixing clade: (i) a genus-level global tree including 1023 genera; and (ii) a regional tree comprising 303 genera, with taxon sampling limited to China. We used the supermatrix approach and performed maximum likelihood analyses on combined matK, rbcL, and trnL-F plastid sequences. We found that the topology of the global and the regional tree of the N-fixing clade were generally congruent. However, whereas relationships among the four orders obtained with the global tree agreed with the accepted topology obtained in focused analyses with more genes, the regional topology obtained different relationships, albeit weakly supported. At a finer scale, the phylogenetic position of the family Myricaceae was found to be sensitive to sampling density. We expect that internal support throughout the phylogeny could be improved with denser taxon sampling. The taxon sampling approach (global vs. regional) did not have a major impact on fine-level branching patterns of the N-fixing clade. Thus, a well-resolved phylogeny with relatively dense taxon sampling strategy at the regional scale appears, in this case, to be a good representation of the overall phylogenetic pattern and could be used in ecological research. Otherwise, the regional tree should be adjusted according to the correspondingly reliable global tree.     

Phylogenomics of <Emphasis Type="Italic">Bartheletia paradoxa</Emphasis> reveals its basal position in Agaricomycotina and that the early evolutionary history of basidiomycetes was rapid and probably not strictly bifurcating

The higher level phylogeny of fungi has been addressed in previous studies, but for those analyses, either taxon sampling or gene sampling was low, or some basal lineages important for the inference of basidiomycete phylogeny were lacking. Here, a phylogenomic analysis based on highly conserved genes and including the enigmatic species Bartheletia paradoxa from Ginkgo biloba is presented. While phylogenetic analyses including also less conserved parts of core eukaryotic genes yielded a basal position for the extremophile genus Wallemia with low support, an exclusion of highly variable parts of these genes suggested Bartheletia paradoxa as the most basal member of the Agaricomycotina, but again with low support. Network analyses suggest a network-like evolution at the base of the Basidiomycota, supported by phylogenies based on single genes and gene clusters with shared topology. When further removing noise by removing poorly resolving genes, strong but not maximum support was obtained for Bartheletia paradoxa being the sister lineage to all other Agaricomycotina. We speculate that the lack of support for the early splits in Agaricomycotina and Basidiomycota can probably be explained by rapid radiation, linked to major evolutionary developments, such as, in the case of Basidiomycota, the advent of basidia in the last common ancestor.     

Phylogenetic assessment of introns and SINEs within the Y chromosome using the cat family felidae as a species tree

The cat family Felidae was used as a species tree to assess the phylogenetic performance of genes, and their embedded SINE elements, within the nonrecombining region of the Y chromosome (NRY). Genomic segments from single-copy X-Y homologs SMCY, UBE1Y, and ZFY (3,604 bp) were amplified in 36 species of cat. These genes are located within the X-degenerate region of the NRY and are thought to be molecular "fossils" that ceased conventional recombination with the X chromosome early within the placental mammal evolution. The pattern and tempo of evolution at these three genes is significant in light of the recent, rapid evolution of the family over approximately 12 Myr and provides exceptional support for each of the eight recognized felid lineages, as well as clear diagnostic substitutions identifying nearly all species. Bootstrap support and Bayesian posterior probabilities are uniformly high for defining each of the eight monophyletic lineages. Further, the preferential use of specific target-site motifs facilitating SINE insertion is empirically supported by sequence analyses of SINEs embedded within the three genes. Target-site insertion is thought to explain the contradiction between intron phylogeny and results of the SMCY SINE phylogeny that unites distantly related species. Overall, our data suggest X-degenerate genes within the NRY are singularly powerful markers and offer a valuable patrilineal perspective in species evolution.     

Origin of land plants revisited in the light of sequence contamination and missing data

Knowing the closest relatives of land plants is key to understanding the complex adaptations to terrestrial life. Unfortunately, multi-gene analyses have yielded highly incongruent results, suggesting, for instance, Charales [1], Zygnematales [2,3], or Coleochaete[4] as the sister-group of land plants. Such controversy may result from the real history of life, in particular closely spaced speciation events, incomplete lineage sorting, gene duplication or horizontal gene transfer. In such cases, the solution resides in improved taxon sampling and sophisticated models of evolution [5]. However, we will show that the quality of data used to infer the phylogeny may also play a major role. In particular, the inclusion of contaminant sequences from other species, and of genes with incomplete taxon sampling explains a large part of the discrepancies observed between various studies [2-4]. The use of a carefully checked and almost complete dataset suggests that land plants are closely related to a group composed of Zygnematales and Coleochaetales.     

Complete cpDNA genome sequence of Smilax china and phylogenetic placement of Liliales--influences of gene partitions and taxon sampling

The complete nucleotide sequence of the chloroplast genome (cpDNA) of Smilax china L. (Smilacaceae) is reported. It is the first complete cp genome sequence in Liliales. Genomic analyses were conducted to examine the rate and pattern of cpDNA genome evolution in Smilax relative to other major lineages of monocots. The cpDNA genomic sequences were combined with those available for Lilium to evaluate the phylogenetic position of Liliales and to investigate the influence of taxon sampling, gene sampling, gene function, natural selection, and substitution rate on phylogenetic inference in monocots. Phylogenetic analyses using sequence data of gene groups partitioned according to gene function, selection force, and total substitution rate demonstrated evident impacts of these factors on phylogenetic inference of monocots and the placement of Liliales, suggesting potential evolutionary convergence or adaptation of some cpDNA genes in monocots. Our study also demonstrated that reduced taxon sampling reduced the bootstrap support for the placement of Liliales in the cpDNA phylogenomic analysis. Analyses of sequences of 77 protein genes with some missing data and sequences of 81 genes (all protein genes plus the rRNA genes) support a sister relationship of Liliales to the commelinids-Asparagales clade, consistent with the APG III system. Analyses of 63 cpDNA protein genes for 32 taxa with few missing data, however, support a sister relationship of Liliales (represented by Smilax and Lilium) to Dioscoreales-Pandanales. Topology tests indicated that these two alignments do not significantly differ given any of these three cpDNA genomic sequence data sets. Furthermore, we found no saturation effect of the data, suggesting that the cpDNA genomic sequence data used in the study are appropriate for monocot phylogenetic study and long-branch attraction is unlikely to be the cause to explain the result of two well-supported, conflict placements of Liliales. Further analyses using sufficient nuclear data remain necessary to evaluate these two phylogenetic hypotheses regarding the position of Liliales and to address the causes of signal conflict among genes and partitions.     

Multiple data sets, high homoplasy, and the phylogeny of softshell turtles (Testudines: Trionychidae)

We present a phylogenetic hypothesis and novel, rank-free classification for all extant species of softshell turtles (Testudines:Trionychidae). Our data set included DNA sequence data from two mitochondrial protein-coding genes and a approximately 1-kb nuclear intron for 23 of 26 recognized species, and 59 previously published morphological characters for a complimentary set of 24 species. The combined data set provided complete taxonomic coverage for this globally distributed clade of turtles, with incomplete data for a few taxa. Although our taxonomic sampling is complete, most of the modern taxa are representatives of old and very divergent lineages. Thus, due to biological realities, our sampling consists of one or a few representatives of several ancient lineages across a relatively deep phylogenetic tree. Our analyses of the combined data set converge on a set of well-supported relationships, which is in accord with many aspects of traditional softshell systematics including the monophyly of the Cyclanorbinae and Trionychinae. However, our results conflict with other aspects of current taxonomy and indicate that most of the currently recognized tribes are not monophyletic. We use this strong estimate of the phylogeny of softshell turtles for two purposes: (1) as the basis for a novel rank-free classification, and (2) to retrospectively examine strategies for analyzing highly homoplasious mtDNA data in deep phylogenetic problems where increased taxon sampling is not an option. Weeded and weighted parsimony, and model-based techniques, generally improved the phylogenetic performance of highly homoplasious mtDNA sequences, but no single strategy completely mitigated the problems of associated with these highly homoplasious data. Many deep nodes in the softshell turtle phylogeny were confidently recovered only after the addition of largely nonhomoplasious data from the nuclear intron.     

Phylogeny, rate variation, and genome size evolution of Pelargonium (Geraniaceae)

The phylogeny of 58 Pelargonium species was estimated using five plastid markers (rbcL, matK, ndhF, rpoC1, trnL-F) and one mitochondrial gene (nad5). The results confirmed the monophyly of three major clades and four subclades within Pelargonium but also indicate the need to revise some sectional classifications. This phylogeny was used to examine karyotype evolution in the genus: plotting chromosome sizes, numbers and 2C-values indicates that genome size is significantly correlated with chromosome size but not number. Accelerated rates of nucleotide substitution have been previously detected in both plastid and mitochondrial genes in Pelargonium, but sparse taxon sampling did not enable identification of the phylogenetic distribution of these elevated rates. Using the multigene phylogeny as a constraint, we investigated lineage- and locus-specific heterogeneity of substitution rates in Pelargonium for an expanded number of taxa and demonstrated that both plastid and mitochondrial genes have had accelerated substitution rates but with markedly disparate patterns. In the plastid, the exons of rpoC1 have significantly accelerated substitution rates compared to its intron and the acceleration was mainly due to nonsynonymous substitutions. In contrast, the mitochondrial gene, nad5, experienced substantial acceleration of synonymous substitution rates in three internal branches of Pelargonium, but this acceleration ceased in all terminal branches. Several lineages also have dN/dS ratios significantly greater than one for rpoC1, indicating that positive selection is acting on this gene, whereas the accelerated synonymous substitutions in the mitochondrial gene are the result of elevated mutation rates.     

Phylogeny and temporal diversification of darters (Percidae: Etheostomatinae)

Discussions aimed at resolution of the Tree of Life are most often focused on the interrelationships of major organismal lineages. In this study, we focus on the resolution of some of the most apical branches in the Tree of Life through exploration of the phylogenetic relationships of darters, a species-rich clade of North American freshwater fishes. With a near-complete taxon sampling of close to 250 species, we aim to investigate strategies for efficient multilocus data sampling and the estimation of divergence times using relaxed-clock methods when a clade lacks a fossil record. Our phylogenetic data set comprises a single mitochondrial DNA (mtDNA) gene and two nuclear genes sampled from 245 of the 248 darter species. This dense sampling allows us to determine if a modest amount of nuclear DNA sequence data can resolve relationships among closely related animal species. Darters lack a fossil record to provide age calibration priors in relaxed-clock analyses. Therefore, we use a near-complete species-sampled phylogeny of the perciform clade Centrarchidae, which has a rich fossil record, to assess two distinct strategies of external calibration in relaxed-clock divergence time estimates of darters: using ages inferred from the fossil record and molecular evolutionary rate estimates. Comparison of Bayesian phylogenies inferred from mtDNA and nuclear genes reveals that heterospecific mtDNA is present in approximately 12.5% of all darter species. We identify three patterns of mtDNA introgression in darters: proximal mtDNA transfer, which involves the transfer of mtDNA among extant and sympatric darter species, indeterminate introgression, which involves the transfer of mtDNA from a lineage that cannot be confidently identified because the introgressed haplotypes are not clearly referable to mtDNA haplotypes in any recognized species, and deep introgression, which is characterized by species diversification within a recipient clade subsequent to the transfer of heterospecific mtDNA. The results of our analyses indicate that DNA sequences sampled from single-copy nuclear genes can provide appreciable phylogenetic resolution for closely related animal species. A well-resolved near-complete species-sampled phylogeny of darters was estimated with Bayesian methods using a concatenated mtDNA and nuclear gene data set with all identified heterospecific mtDNA haplotypes treated as missing data. The relaxed-clock analyses resulted in very similar posterior age estimates across the three sampled genes and methods of calibration and therefore offer a viable strategy for estimating divergence times for clades that lack a fossil record. In addition, an informative rank-free clade-based classification of darters that preserves the rich history of nomenclature in the group and provides formal taxonomic communication of darter clades was constructed using the mtDNA and nuclear gene phylogeny. On the whole, the appeal of mtDNA for phylogeny inference among closely related animal species is diminished by the observations of extensive mtDNA introgression and by finding appreciable phylogenetic signal in a modest sampling of nuclear genes in our phylogenetic analyses of darters.     

Multidimensional vector space representation for convergent evolution and molecular phylogeny

With growing amounts of genome data and constant improvement of models of molecular evolution, phylogenetic reconstruction became more reliable. However, our knowledge of the real process of molecular evolution is still limited. When enough large-sized data sets are analyzed, any subtle biases in statistical models can support incorrect topologies significantly because of the high signal-to-noise ratio. We propose a procedure to locate sequences in a multidimensional vector space (MVS), in which the geometry of the space is uniquely determined in such a way that the vectors of sequence evolution are orthogonal among different branches. In this paper, the MVS approach is developed to detect and remove biases in models of molecular evolution caused by unrecognized convergent evolution among lineages or unexpected patterns of substitutions. Biases in the estimated pairwise distances are identified as deviations (outliers) of sequence spatial vectors from the expected orthogonality. Modifications to the estimated distances are made by minimizing an index to quantify the deviations. In this way, it becomes possible to reconstruct the phylogenetic tree, taking account of possible biases in the model of molecular evolution. The efficacy of the modification procedure was verified by simulating evolution on various topologies with rate heterogeneity and convergent change. The phylogeny of placental mammals in previous analyses of large data sets has varied according to the genes being analyzed. Systematic deviations caused by convergent evolution were detected by our procedure in all representative data sets and were found to strongly affect the tree structure. However, the bias correction yielded a consistent topology among data sets. The existence of strong biases was validated by examining the sites of convergent evolution between the hedgehog and other species in mitochondrial data set. This convergent evolution explains why it has been difficult to determine the phylogenetic placement of the hedgehog in previous studies.