Quantification of the success of phylogenetic inference in simulations

For phylogenetic simulation studies, the accuracy of topological reconstruction obtained from different data matrices or different methods of phylogenetic inference generally needs to be quantified. Two components of performance within this context are: (1) how the inferred tree topology matches or conflicts with the correct tree topology, and (2) the branch support assigned to both correctly and incorrectly resolved clades. We present a method (averaged overall success of resolution) that incorporates both of these components. Branch support is incorporated in the averaged overall success of resolution by linearly scaling the observed support relative to that conferred by uncontradicted synapomorphies. We believe that this method represents an improvement relative to the commonly used approaches of quantifying the percentage of clades that are correctly resolved in the inferred trees or presenting the Robinson–Foulds distance between the inferred trees and the correct tree. In contrast to Bremer support, the averaged overall success of resolution may be applied equally well to distance, likelihood and parsimony analyses. © The Willi Hennig Society 2006.     

Molecular phylogeny of plethodonine salamanders and hylid frogs: statistical analysis of protein comparisons.

The bootstrapping method of determining confidence in the topology of phylogenetic trees has been applied to electrophoretic protein data for two groups of amphibians: salamanders of two North American genera (Aneides and Plethodon) of the tribe Plethodontini and Holarctic hylid frogs. Some current methods of phylogenetic reconstruction for electrophoretic protein data have been evaluated by comparing the trees obtained from molecular data sets with available morphological data. Molecular data on the phylogenetic relationships of Aneides and Plethodon, data obtained from electrophoretic and immunological studies, indicate that Aneides probably was derived from western Plethodon subsequent to the separation of eastern and western Plethodon. Thus Plethodon very likely is a paraphyletic genus. The extremely low rate of morphological evolution in Plethodon compared with that in Aneides causes difficulty in indicating their evolutionary relationships taxonomically because there are no synapomorphic morphological characters that define either eastern or western Plethodon, whereas there are several for the genus Aneides. Thus molecular data alone probably indicate the evolutionary relationships of the species in these genera. Highton and Larson's (1979) arrangement of species of Plethodon into eight species groups is supported. The topologies of the unweighted pair-group method using arithmetic means (UPGMA) and distance Wagner trees were compared with independent morphological and molecular data on the relationships of the 28 plethodonine species. It was found that UPGMA trees indicate relationships that are more in agreement with other information than are those provided by distance Wagner trees. The use of the bootstrap technique indicates that the topologies of UPGMA trees are better supported statistically than are the topologies of distance Wagner trees. Moreover, different addition criteria produce a variety of distance Wagner trees with different topologies, each with several groupings that are not supported statistically. It is concluded that considerable caution should be used in interpreting the topology of distance Wagner trees. Very similar results were obtained with a second data set on 30 taxa of Holarctic hylid frogs. Trees obtained by the neighbor-joining method are more in agreement with UPGMA phenograms and other data, so this method of phylogenetic reconstruction may be useful to systematists not willing to assume constant rates of evolution.(ABSTRACT TRUNCATED AT 400 WORDS)     

Complete mitochondrial DNA sequence of the fat dormouse, Glis glis: further evidence of rodent paraphyly

The complete mitochondrial genome of the fat dormouse, Glis glis, has been sequenced (16,602 bp). A total of 23 complete mitochondrial mammalian genomes have been taken into account for phylogenetic reconstruction. Phylogenetic analyses were performed with parsimony, distance (stationary Markov model), and maximum-likelihood methods. In all cases, data strongly support the paraphyly of rodents, with dormouse and guinea pig in a different clade from rat and mouse, reaching bootstrap values of 95%. Rodent monophyly and the existence of Glires (Rodentia and Lagomorpha) are weakly supported, with maximum bootstrap values of 11% and 8.6%, respectively. This result agrees with the analyses of isochore patterns in the nuclear genome and the B2 and B2-like retroposons, which show a close relationship between dormice and guinea pigs rather than between dormice and rats and mice.      

The serine repeat antigen (SERA) gene family phylogeny in Plasmodium: the impact of GC content and reconciliation of gene and species trees

Plasmodium falciparum is the parasite responsible for the most acute form of malaria in humans. Recently, the serine repeat antigen (SERA) in P. falciparum has attracted attention as a potential vaccine and drug target, and it has been shown to be a member of a large gene family. To clarify the relationships among the numerous P. falciparum SERAs and to identify orthologs to SERA5 and SERA6 in Plasmodium species affecting rodents, gene trees were inferred from nucleotide and amino acid sequence data for 33 putative SERA homologs in seven different species. (A distance method for nucleotide sequences that is specifically designed to accommodate differing GC content yielded results that were largely compatible with the amino acid tree. Standard-distance and maximum-likelihood methods for nucleotide sequences, on the other hand, yielded gene trees that differed in important respects.) To infer the pattern of duplication, speciation, and gene loss events in the SERA gene family history, the resulting gene trees were then "reconciled" with two competing Plasmodium species tree topologies that have been identified by previous phylogenetic studies. Parsimony of reconciliation was used as a criterion for selecting a gene tree/species tree pair and provided (1) support for one of the two species trees and for the core topology of the amino acid-derived gene tree, (2) a basis for critiquing fine detail in a poorly resolved region of the gene tree, (3) a set of predicted "missing genes" in some species, (4) clarification of the relationship among the P. falciparum SERA, and (5) some information about SERA5 and SERA6 orthologs in the rodent malaria parasites. Parsimony of reconciliation and a second criterion--implied mutational pattern at two key active sites in the SERA proteins-were also seen to be useful supplements to standard "bootstrap" analysis for inferred topologies.     

The role of phylogeny in desert rodent community assembly

Recent advances in molecular genetics and phylogenetic reconstruction have the potential to transform ecology by providing new insights into the historical evolution of ecological communities. This study by Stevens and collaborators complements decades of previous research on desert rodents, by combining data from a field study and a phylogenetic tree for Mojave Desert rodents to address patterns and processes of community assembly. The number of coexisting rodent species is positively correlated, and the average phylogenetic distance among these species is negatively correlated with perennial plant species richness. As rodent species diversity increases along a gradient of increasing environmental heterogeneity, communities are composed of increasingly related species: there is a consistent pattern of phylogenetic structure from over-dispersed through random to clumped. I discuss this pattern in the light of complementary results of previous studies. This paper is noteworthy for calling attention to still unanswered questions about how the historical events of speciation, colonization, extinction, and trait evolution and their relationship to past climates and vegetation have given rise to current patterns of community organization.     

Reconstructing the phylogeny of "Buarremon" brush-finches and near relatives (Aves, Emberizidae) from individual gene trees

Gene trees are often assumed to be equivalent to species trees, but processes such as incomplete lineage sorting can generate incongruence among gene topologies and analyzing multilocus data in concatenated matrices can be prone to systematic errors. Accordingly, a variety of new methods have been developed to estimate species trees using multilocus data sets. Here, we apply some of these methods to reconstruct the phylogeny of Buarremon and near relatives, a group in which phylogenetic analyses of mitochondrial DNA sequences produced results that were inconsistent with relationships implied by a taxonomy based on variation in external phenotype. Gene genealogies obtained for seven loci (one mitochondrial, six nuclear) were varied, with some supporting and some rejecting the monophyly of Buarremon. Overall, our species-tree analyses tended to support a monophyletic Buarremon, but due to lack of congruence between methodologies, resolution of the phylogeny of this group remains uncertain. More generally, our study indicates that the number of individuals sampled can have an important effect on phylogenetic reconstruction, that the use of seven markers does not guarantee obtaining a strongly-supported species tree, and that methods for species-tree reconstruction can produce different results using the same data; these are important considerations for researchers using these new phylogenetic approaches in other systems.     

MOLECULAR SYSTEMATICS OF THE FLORIDEOPHYCEAE (RHODOPHYTA) USING NUCLEAR LARGE AND SMALL SUBUNIT rDNA SEQUENCE DATA

Sequence data are presented for approximately 85% of the nuclear large subunit (LSU) rDNA gene for one member of the Bangiophyceae and 47 members of the Florideophyceae, the latter representing all but one of the currently recognized florideophyte orders. Distance, parsimony, and maximum likelihood analyses of these data were used to generate phylogenetic trees, and bootstrap resampling was implemented to infer robustness for distance and parsimony results. LSU phylogenies were congruent with published nuclear small subunit (SSU) rDNA results in that four higher level florideophyte lineages were resolved: lineage 1, containing the order Hildenbrandiales; lineage 2, recovered only under distance analysis, composed of the orders Acrochaetiales, Balliales, Batrachospermales, Corallinales, Nemaliales, Palmariales, and Rhodogorgonales; lineage 3, containing the Ahnfeltiales; and lineage 4, composed of the orders Bonnemaisoniales, Ceramiales, Gelidiales, Gigartinales, Gracilariales, Halymeniales, Plocamiales, and Rhodymeniales. Analyses were also performed on a combined LSU–SSU data set and an SSU-only data set to account for differences in taxon sampling relative to published studies using this latter gene. Combined LSU–SSU analyses resulted in phylogenetic trees of similar topology and support to those obtained from LSU-only analyses. Phylogenetic trees produced from SSU-only analyses differed somewhat in particulars of branching within lineages 2 and 4 but overall were congruent with the LSU-only and combined LSU–SSU results. We close with a discussion of the phylogenetic potential that the LSU has displayed thus far for resolving relationships within the Florideophyceae.     

长爪沙鼠线粒体DNA控制区全序列测定及分析

目的对长爪沙鼠线粒体DNA控制区全序列进行测定,并对其进行鉴定及进化分析。方法根据长爪沙鼠已知基因序列设计引物,采用PCR产物测序法,对所得的片段进行测序鉴定。结合已公布啮齿类动物D-loop区序列,分析其碱基组成、遗传距离、并基于最小进化法和UPGMA法构建系统进化树。结果获得长爪沙鼠D-loop区序列,其与家鼠、小家鼠和仓鼠平均同源性为58%;碱基组成分析显示,长爪沙鼠与啮齿类动物有相似的碱基组成和碱基偏离,其A-skew和G-skew分别为0.0047和-0.28。进化分析结果显示,长爪沙鼠与家鼠（0.35）、黑家鼠（0.38）和仓鼠（0.39）具有较近的遗传距离,其分化顺序为跳鼠、蔗鼠、长爪沙鼠、仓鼠、家鼠和小家鼠。结论本研究获得长爪沙鼠D-loop区全序列,确定了长爪沙鼠与仓鼠、家鼠、小家鼠及其它啮齿动物的进化关系,为长爪沙鼠进化研究、线粒体的结构和功能研究奠定基础。     

Molecular systematic studies of eubacteria, using sigma70-type sigma factors of group 1 and group 2.

Sigma factors of the sigma70 family were used as a phylogenetic tool to compare evolutionary relationships among eubacteria. Several new sigma factor genes were cloned and sequenced to increase the variety of available sequences. Forty-two group 1 sigma factor sequences of various species were analyzed with the help of a distance matrix method to establish a phylogenetic tree. The tree derived by using sigma factors yielded subdivisions, including low-G+C and high-G+C gram-positive bacteria, cyanobacteria, and the alpha, beta, gamma, and delta subdivisions of proteobacteria, consistent with major bacterial groups found in trees derived from analyses with other molecules. However, some groupings (e.g., the chlamydiae, mycoplasmas, and green sulfur bacteria) are found in different positions than for trees obtained by using other molecular markers. A direct comparison to the most extensively used molecule in systematic studies, small-subunit rRNA, was made by deriving trees from essentially the same species set and using similar phylogenetic methods. Differences and similarities based on the two markers are discussed. Additionally, 31 group 2 sigma factors were analyzed in combination with the group 1 proteins in order to detect functional groupings of these alternative sigma factors. The data suggest that promoters recognized by the major vegetative sigma factors of eubacteria will contain sequence motifs and spacing very similar to those for the sigma70 sigma factors of Escherichia coli.     

Topology improves phylogenetic motif functional site predictions

Prediction of protein functional sites from sequence-derived data remains an open bioinformatics problem. We have developed a phylogenetic motif (PM) functional site prediction approach that identifies functional sites from alignment fragments that parallel the evolutionary patterns of the family. In our approach, PMs are identified by comparing tree topologies of each alignment fragment to that of the complete phylogeny. Herein, we bypass the phylogenetic reconstruction step and identify PMs directly from distance matrix comparisons. In order to optimize the new algorithm, we consider three different distance matrices and 13 different matrix similarity scores. We assess the performance of the various approaches on a structurally nonredundant data set that includes three types of functional site definitions. Without exception, the predictive power of the original approach outperforms the distance matrix variants. While the distance matrix methods fail to improve upon the original approach, our results are important because they clearly demonstrate that the improved predictive power is based on the topological comparisons. Meaning that phylogenetic trees are a straightforward, yet powerful way to improve functional site prediction accuracy. While complementary studies have shown that topology improves predictions of protein-protein interactions, this report represents the first demonstration that trees improve functional site predictions as well.     

Phylogenetic mixtures on a single tree can mimic a tree of another topology

Phylogenetic mixtures model the inhomogeneous molecular evolution commonly observed in data. The performance of phylogenetic reconstruction methods where the underlying data are generated by a mixture model has stimulated considerable recent debate. Much of the controversy stems from simulations of mixture model data on a given tree topology for which reconstruction algorithms output a tree of a different topology; these findings were held up to show the shortcomings of particular tree reconstruction methods. In so doing, the underlying assumption was that mixture model data on one topology can be distinguished from data evolved on an unmixed tree of another topology given enough data and the "correct" method. Here we show that this assumption can be false. For biologists, our results imply that, for example, the combined data from two genes whose phylogenetic trees differ only in terms of branch lengths can perfectly fit a tree of a different topology.     

Phylogenomic analysis resolves the formerly intractable adaptive diversification of the endemic clade of east Asian Cyprinidae (Cypriniformes)

Despite their great diversity and biological importance, evolutionary relationships among the endemic clade of East Asian Cyprinidae remain ambiguous. Understanding the phylogenetic history of this group involves many challenges. For instance, ecomorphological convergence may confound morphology-based phylogenetic inferences, and previous molecular phylogenetic studies based on single genes have often yielded contradictory and poorly supported trees. We assembled a comprehensive data matrix of 100 nuclear gene segments (～ 71132 base pairs) for representative species of the endemic East Asian cyprinid fauna and recovered a robust phylogeny from this genome-wide signal supported by multiple analytical methods, including maximum parsimony, maximum likelihood and Bayesian inference. Relaxed molecular clock analyses indicated species radiations of this clade concentrated at approximately 1.9-7.6 MYA. We provide evidence that the bursts of diversification in this fauna are directly linked to major paleoenvironmental events associated with monsoon evolution occurring from late Miocene to Pliocene. Ancestral state reconstruction reveals convergent morphological characters are hypothesized to be independent products of similar selective pressures in ecosystems. Our study is the first comprehensive phylogenetic study of the enigmatic East-Asian cyprinids. The explicit molecular phylogeny provides a valuable framework for future research in genome evolution, adaptation and speciation of cyprinids.     

用线粒体tRNA基因探讨现存两栖动物三个目间系统发生关系

现存两栖类3个目的系统发生关系仍然没有统一意见,最广泛被接受的假说是单系起源,并且无尾类和有尾类为姐妹群关系而排斥蚓螈类(蛙类假说)。然而,这一假说一直存在争议。我们在测定了泽蛙线粒体基因组全序列的基础上,与已知其他的6种两栖类进行详细的比较分析,同时选择了11种高等脊椎动物的线粒体全基因序列,以硬骨鱼类作外群,用22个tRNA基因合并数据进行系统发生重建分析,结果表明MP、ML树都强力地支持现生两栖类动物为单系群,并且有尾目和蚓螈目为姐妹群关系。这个结果与蛙类假说是相矛盾的,与Bolt(1991)在形态学基础上提出的有尾类和蚓螈类为姐妹群关系的假说相一致,并得到建立在线粒体和核rRNA基因数据基础上的许多分子研究的支持。另外还探讨了本结果与前人的研究不一致的原因,以及利用线粒体全基因序列进行系统发生分析可能存在的偏差。     

DNA barcoding of Murinae (Rodentia: Muridae) and Arvicolinae (Rodentia: Cricetidae) distributed in China

Identification of rodents is very difficult mainly due to high similarities in morphology and controversial taxonomy. In this study, mitochondrial cytochrome oxidase subunit I (COI) was used as DNA barcode to identify the Murinae and Arvicolinae species distributed in China and to facilitate the systematics studies of Rodentia. In total, 242 sequences (31 species, 11 genera) from Murinae and 130 sequences (23 species, 6 genera) from Arvicolinae were investigated, of which 90 individuals were novel. Genetic distance, threshold method, tree‐based method, online BLAST and BLOG were employed to analyse the data sets. There was no obvious barcode gap. The average K2P distance within species and genera was 2.10% and 12.61% in Murinae, and 2.86% and 11.80% in Arvicolinae, respectively. The optimal threshold was 5.62% for Murinae and 3.34% for Arvicolinae. All phylogenetic trees exhibited similar topology and could distinguish 90.32% of surveyed species in Murinae and 82.60% in Arvicolinae with high support values. BLAST analyses yielded similar results with identification success rates of 92.15% and 93.85% for Murinae and Arvicolinae, respectively. BLOG successfully authenticated 100% of detected species except Leopoldamys edwardsi based on the latest taxonomic revision. Our results support the species status of recently recognized Micromys erythrotis, Eothenomys tarquinius and E. hintoni and confirm the important roles of comprehensive taxonomy and accurate morphological identification in DNA barcoding studies. We believe that, when proper analytic methods are applied or combined, DNA barcoding could serve as an accurate and effective species identification approach for Murinae and Arvicolinae based on a proper taxonomic framework.     

Phylogenetic signal, evolutionary process, and rate

A recent advance in the phylogenetic comparative analysis of continuous traits has been explicit, model-based measurement of "phylogenetic signal" in data sets composed of observations collected from species related by a phylogenetic tree. Phylogenetic signal is a measure of the statistical dependence among species' trait values due to their phylogenetic relationships. Although phylogenetic signal is a measure of pattern (statistical dependence), there has nonetheless been a widespread propensity in the literature to attribute this pattern to aspects of the evolutionary process or rate. This may be due, in part, to the perception that high evolutionary rate necessarily results in low phylogenetic signal; and, conversely, that low evolutionary rate or stabilizing selection results in high phylogenetic signal (due to the resulting high resemblance between related species). In this study, we use individual-based numerical simulations on stochastic phylogenetic trees to clarify the relationship between phylogenetic signal, rate, and evolutionary process. Under the simplest model for quantitative trait evolution, homogeneous rate genetic drift, there is no relation between evolutionary rate and phylogenetic signal. For other circumstances, such as functional constraint, fluctuating selection, niche conservatism, and evolutionary heterogeneity, the relationship between process, rate, and phylogenetic signal is complex. For these reasons, we recommend against interpretations of evolutionary process or rate based on estimates of phylogenetic signal.     

Concordance analysis in mitogenomic phylogenetics

Here I advocate the utility of Bayesian concordance analysis as a mechanism for exploring the magnitude and source of phylogenetic signal in concatenated mitogenomic phylogenetic studies. While typically applied to the study of independently evolving gene trees, Bayesian concordance analysis can also be applied to linked, but individually analyzed, gene regions using a prior probability that reflects the expectation of similar phylogenetic reconstructions. For true branches in the mitogenomic tree, concordance factors should represent the number of gene regions that contain phylogenetic signal for a particular clade. As a demonstration of the application of Bayesian concordance analysis to empirical data, I analyzed two different salamander (Hynobiidae and Plethodontidae) mitogenomic data sets using a gene-based partitioning strategy. The results revealed many strongly supported clades in the concatenated trees that have high concordance factors, permitting the inference that these are robustly resolved through phylogenetic signal distributed across the mitogenome. In contrast, a number of strongly supported clades in the concatenated tree received low concordance factors, indicating that their reconstruction is either driven primarily by phylogenetic signal in a small number of gene regions, or that they are inconsistent reconstructions influenced by properties of the data that can produce inaccurate trees (e.g., compositional bias, selection, etc.). Exploration of the Bayesian joint posterior distribution of trees highlighted partitions that contribute phylogenetic information to similar clade reconstructions. This approach was particularly insightful in the hynobiid data, where different combinations of genes were identified that support alternative tree reconstructions. Concatenated analysis of these different subsets of genes highlighted through Bayesian concordance analysis produced strongly supported and contrasting trees, demonstrating the potential for inconsistency in concatenated mitogenomic phylogenetics. The overall results presented here suggest that Bayesian concordance analysis can serve as an effective exploration of the influence of different gene regions in mitogenomic (and other organellar genomic) phylogenetic studies.     

Nonhomogeneous model of sequence evolution indicates independent origins of primary endosymbionts within the enterobacteriales (gamma-Proteobacteria)

Standard methods of phylogenetic reconstruction are based on models that assume homogeneity of nucleotide composition among taxa. However, this assumption is often violated in biological data sets. In this study, we examine possible effects of nucleotide heterogeneity among lineages on the phylogenetic reconstruction of a bacterial group that spans a wide range of genomic nucleotide contents: obligately endosymbiotic bacteria and free-living or commensal species in the gamma-Proteobacteria. We focus on AT-rich primary endosymbionts to better understand the origins of obligately intracellular lifestyles. Previous phylogenetic analyses of this bacterial group point to the importance of accounting for base compositional variation in estimating relationships, particularly between endosymbiotic and free-living taxa. Here, we develop an approach to compare susceptibility of various phylogenetic reconstruction methods to the effects of nucleotide heterogeneity. First, we identify candidate trees of gamma-Proteobacteria groEL and 16S rRNA using approaches that assume homogeneous and stationary base composition, including Bayesian, maximum likelihood, parsimony, and distance methods. We then create permutations of the resulting candidate trees by varying the placement of the AT-rich endosymbiont Buchnera. These permutations are evaluated under the nonhomogeneous and nonstationary maximum likelihood model of Galtier and Gouy, which allows equilibrium base content to vary among examined lineages. Our results show that commonly used phylogenetic methods produce incongruent trees of the Enterobacteriales, and that the placement of Buchnera is especially unstable. However, under a nonhomogeneous model, various groEL and 16S rRNA phylogenies that separate Buchnera from other AT-rich endosymbionts (Blochmannia and Wigglesworthia) have consistently and significantly higher likelihood scores. Blochmannia and Wigglesworthia appear to have evolved from secondary endosymbionts, and represent an origin of primary endosymbiosis that is independent from Buchnera. This application of a nonhomogeneous model offers a computationally feasible way to test specific phylogenetic hypotheses for taxa with heterogeneous and nonstationary base composition.     

Morphological and molecular phylogenetics in the genus Leptopilina (Hymenoptera: Cynipoidea: Eucoilidae)

We have conducted cladistic analyses of the genus Leptopilina , a group of Drosophila parasitoids studied intensively by (behavioural) ecologists. Twenty-three morphological characters were scored in ten Leptopilina and two outgroup species. At the same time, DNA sequences for the second ribosomal internal transcribed spacer (ITS2) were gathered for eight Leptopilina and one outgroup species. Both data sets yielded phylogenetic trees which were largely compatible. A 'total evidence' analysis resulted in a single tree that provides a relatively robust phylogenetic reconstruction of the group, which may serve as a basis for historically interpreting the distribution of ecological and behavioural traits.     

Phylogenetic relationships of Steinernema Travassos, 1927 (Nematoda: Cephalobina: Steinernematidae) based on nuclear, mitochondrial and morphological data

Entomopathogenic nematodes of the genus Steinernema are lethal parasites of insects that are used as biological control agents of several lepidopteran, dipteran and coleopteran pests. Phylogenetic relationships among 25 Steinernema species were estimated using nucleotide sequences from three genes and 22 morphological characters. Parsimony analysis of 28S (LSU) sequences yielded a well-resolved phylogenetic hypothesis with reliable bootstrap support for 13 clades. Parsimony analysis of mitochondrial DNA sequences (12S rDNA and cox 1 genes) yielded phylogenetic trees with a lower consistency index than for LSU sequences, and with fewer reliably supported clades. Combined phylogenetic analysis of the 3-gene dataset by parsimony and Bayesian methods yielded well-resolved and highly similar trees. Bayesian posterior probabilities were high for most clades; bootstrap (parsimony) support was reliable for approximately half of the internal nodes. Parsimony analysis of the morphological dataset yielded a poorly resolved tree, whereas total evidence analysis (molecular plus morphological data) yielded a phylogenetic hypothesis consistent with, but less resolved than trees inferred from combined molecular data. Parsimony mapping of morphological characters on the 3-gene trees showed that most structural features of steinernematids are highly homoplastic. The distribution of nematode foraging strategies on these trees predicts that S. hermaphroditum, S. diaprepesi and S. longicaudum (US isolate) have cruise forager behaviours.     

phangorn: phylogenetic analysis in R

SUMMARY: phangorn is a package for phylogenetic reconstruction and analysis in the R language. Previously it was only possible to estimate phylogenetic trees with distance methods in R. phangorn, now offers the possibility of reconstructing phylogenies with distance based methods, maximum parsimony or maximum likelihood (ML) and performing Hadamard conjugation. Extending the general ML framework, this package provides the possibility of estimating mixture and partition models. Furthermore, phangorn offers several functions for comparing trees, phylogenetic models or splits, simulating character data and performing congruence analyses. AVAILABILITY: phangorn can be obtained through the CRAN homepage http://cran.r-project.org/web/packages/phangorn/index.html. phangorn is licensed under GPL 2.