常用系统发育树构建算法和软件鸟瞰

系统发育树又称进化树、生命树等,在达尔文的"进化论"一书中首次出现,之后系统发育树的重构被广大生物学家所接受。该文阐述了构建系统发育树的基本流程,对目前用于构建系统发育树的四类算法(距离法、最大简约法、最大似然法和贝叶斯法)进行了详细地分析和比较,并介绍了一些常用系统发育树构建和分析软件(PHYLIP、MEGA、MrBayes)的特点。     

基于孕激素受体基因的类人猿亚目的分子系统学研究

通过对类人猿亚目中部分种类的孕激素受体基因进行分析,重建类人猿亚目的 系统发育关系.扩增并测定了来源于14个属的类人猿亚目物种的孕激素受体编码区序列,并基于这一序列数据,分别采用邻接法、最大简约法和最大似然法重建了系统发育关系.除了阔鼻下目,3种方法构建的系统发生树的拓扑结构类似且各节点支持率高.重建的人猿超科和猴超科内部亲缘关系支持多数人所认可的分类系统.本研究为黑猩猩和人的姐妹群关系提供了证据,提示黑猩猩比大猩猩或其他猿猴更接近人类.阔鼻下目中蜘蛛猴科、卷尾猴科和僧面猴科三者之间的系统发育关系在本研究中未得到很好辨析.     

基于Cytb和COⅠ基因部分序列的15种麻蝇之间的系统发育关系(双翅目:麻蝇科)(英文)

本研究通过测序Cytb基因和COⅠ基因的部分序列来推定15种麻蝇之间的系统发育关系。在世界麻蝇名录中,本研究的15种麻蝇能够代表麻蝇属Sarcophaga的6个亚属。连接序列(972bp)被用于系统发育分析;分析方法包括了了最大简约法、最大似然法以及贝叶斯法。我们的结果提示了亚麻蝇亚属Parasarcophaga、别麻蝇亚属Boettcherisca以及红麻蝇亚属Liopygia的单系性,同时也表明蛇麻蝇亚属Liosarcophaga和德麻蝇亚属Pandelleisca并不是单源的。不过,目前的研究并不能分辨野德麻蝇S.(Pandelleisca)similis和峨眉叉麻蝇S.(Robineauella)coei的系统发育位置。此外,最大简约分析和似然功能分析在scopariiformis-iwuensis进化枝和polystylata-hui进化枝的关系上产生了不一致的系统发育推断。因此,后续研究不仅需要其他的分子标记,也需要更多的分类取样。     

外群选择对隧蜂科（膜翅目：蜜蜂总科）系统重建的影响

外群用于给树附根和推断祖先性状状态。通常,来自内群的姐妹群中的多个分类单元被共同选择作为外群。为了在经验上验证这一方法, 我们采用了3种外群选择策略： 姐妹群中的单一分类单元, 姐妹群中的多个分类单元和连续姐妹群中的多个分类单元。以隧蜂科（膜翅目： 蜜蜂总科）的系统发育重建为例, 我们评估了这3种策略对树拓扑结构的影响, 包括最大似然树、 最大简约树和贝叶斯树。初步结果表明： 相比其他两种策略, 采用姐妹群中的多个分类单元作为外群更有利于系统发育重建得到现已被广泛认可的隧蜂科系统发育关系; 相比最大似然法和贝叶斯法, 虽然隧蜂科系统发育关系没有被很好地解决, 但最大简约法在不同外群选择策略下得到了较为一致的拓扑结构     

贝叶斯推论及其在百合目分子系统学中的应用

贝叶斯推论作为进化生物学的最新进展,在适用复杂模型、大型数据集、计算速度和结果容易解释等方面明显优于其它算法。本文简要介绍了贝叶斯推论原理及其在分子进化和系统发育研究中的重要性,并使用该方法对百合目主要类群的系统发育关系进行了重建。结果显示,百合目rbcL基因最适合的DNA进化模型为GTR I G,贝叶斯法与距离法和最大简约法构建的系统发育树拓扑结构相似,没有显著差异,但是分辨率和支持率明显比后者高。贝叶斯分析结果显示,百合目内划分的7个科,除Smilacaceae科外,其余各科均为高后验概率(PP=1·0)支持的单系类群;文中作者还对各科间的系统关系进行了探讨。     

基于Cytb和COⅠ基因部分序列的15种麻蝇之间的系统发育关系(双翅目：麻蝇科)

本研究通过测序Cytb基因和COⅠ基因的部分序列来推定15种麻蝇之间的系统发育关系。在世界麻蝇名录中,本研究的15种麻蝇能够代表麻蝇属Sarcophaga的6个亚属。连接序列（972 bp）被用于系统发育分析;分析方法包括了了最大简约法、最大似然法以及贝叶斯法。我们的结果提示了亚麻蝇亚属Parasarcophaga、别麻蝇亚属Boettcherisca以及红麻蝇亚属Liopygia的单系性,同时也表明蛇麻蝇亚属Liosarcophaga和德麻蝇亚属Pandelleisca并不是单源的。不过,目前的研究并不能分辨野德麻蝇S. (Pandelleisca) similis和峨眉叉麻蝇S. (Robineauella) coei的系统发育位置。此外,最大简约分析和似然功能分析在scopariiformis-iwuensis进化枝和polystylata-hui进化枝的关系上产生了不一致的系统发育推断。因此,后续研究不仅需要其他的分子标记,也需要更多的分类取样。     

系统发育研究中直同源预测方法的选择策略

测序技术的不断发展和价格的日益降低使得系统发育组学更深层次的研究成为可能。在系统发育组学分析中,至关重要的步骤是直同源预测,这是因为进行系统发育重建的先决条件是进行比对的基因是直同源的。这里我们简单地回顾了直同源的定义和直同源预测的不同方法,与此同时,我们还给出了一些选择更合适的直同源预测方法的建议。     

基于线粒体ND2基因的中国斑翅蝗科部分种类分子系统学研究（直翅目：蝗总科）

本文的目的是通过对斑翅蝗科部分种类的线粒体ND2基因进行分析,重建斑翅蝗科昆虫的系统发育关系,并探讨分子系统发育关系和传统分类结果的异同。扩增并测定了我国斑翅蝗科10属16种蝗虫的线粒体ND2全基因1 023 bp的序列,对序列的碱基组成、转换颠换、系统发育信号等进行了分析。并基于ND2全基因序列数据,分别采用邻接法（NJ）、最简约法（MP）、最大似然法（ML）和贝叶斯法重建了10属16种蝗虫的系统发育关系。结果表明：斑翅蝗科蝗虫ND2全基因A+T含量平均为74.6%;痂蝗亚科和异痂蝗亚科没能得到区分,建议合并为一个亚科;而斑翅蝗亚科和飞蝗亚科的分类地位还存在争议。     

线粒体D-loop序列变异与东方鲀属鱼类系统发育

东方鲀属的红鳍东方鲀（Takifugu rubripes）是后基因组时代的一种重要模式生物。本研究中,利用东方鲀属11种鱼类（18尾）的D-loop基因序列,对东方鲀属鱼类的系统发育关系进行研究。经序列比对排定后,分析中D-loop序列有841个位点,其中395个位点为可变位点,267个位点为系统发育信息位点。分别采用邻接法（NJ）、最大简约法（MP）、最大似然法（ML）和贝叶斯方法构建了分子系统树。研究结果表明：（1）东方鲀属鱼类为一单系类群;（2）由横纹东方鲀（T. oblongus）和铅点东方鲀（T. alboplumbeus）构成的姊妹群位于这个类群的基部。此外,本属鱼类物种分类现状还需要进一步的澄清。     

基于ITS序列分析仲彬草属植物的亲缘关系

以旱雀麦为外类群,用PAUP 4.0b10软件并采用最大简约法和邻接法对11份仲彬草属物种的ITS区序列进行系统发育分析,两种方法得到的系统发育树基本一致。结果表明:(1)整个ITS序列长度变异范围为596~601 bp;G C含量在所有ITS中的变化范围为61.20%~62.44%;序列间的遗传分化距离为0.003~0.033,平均值为0.015;(2)疏花仲彬草和塔克拉干仲彬草2个物种聚为一支,位于系统发育树的底部,在最大简约法和邻接法分析中分别获得78%和82%的自展支持率,它们之间的亲缘关系较近;(3)形态相似、地理分布一致的物种有聚在一起的倾向,表现出较近的亲缘关系;(4)ITS区序列分析的结果与细胞学、形态学的研究结果基本一致,因此ITS区序列分析能反映仲彬草属种间关系。     

Molecular systematics of the Eastern Fence Lizard (Sceloporus undulatus): a comparison of Parsimony, Likelihood, and Bayesian approaches

Phylogenetic analysis of large datasets using complex nucleotide substitution models under a maximum likelihood framework can be computationally infeasible, especially when attempting to infer confidence values by way of nonparametric bootstrapping. Recent developments in phylogenetics suggest the computational burden can be reduced by using Bayesian methods of phylogenetic inference. However, few empirical phylogenetic studies exist that explore the efficiency of Bayesian analysis of large datasets. To this end, we conducted an extensive phylogenetic analysis of the wide-ranging and geographically variable Eastern Fence Lizard (Sceloporus undulatus). Maximum parsimony, maximum likelihood, and Bayesian phylogenetic analyses were performed on a combined mitochondrial DNA dataset (12S and 16S rRNA, ND1 protein-coding gene, and associated tRNA; 3,688 bp total) for 56 populations of S. undulatus (78 total terminals including other S. undulatus group species and outgroups). Maximum parsimony analysis resulted in numerous equally parsimonious trees (82,646 from equally weighted parsimony and 335 from weighted parsimony). The majority rule consensus tree derived from the Bayesian analysis was topologically identical to the single best phylogeny inferred from the maximum likelihood analysis, but required approximately 80% less computational time. The mtDNA data provide strong support for the monophyly of the S. undulatus group and the paraphyly of "S. undulatus" with respect to S. belli, S. cautus, and S. woodi. Parallel evolution of ecomorphs within "S. undulatus" has masked the actual number of species within this group. This evidence, along with convincing patterns of phylogeographic differentiation suggests "S. undulatus" represents at least four lineages that should be recognized as evolutionary species.     

Applying a multiobjective metaheuristic inspired by honey bees to phylogenetic inference

The development of increasingly popular multiobjective metaheuristics has allowed bioinformaticians to deal with optimization problems in computational biology where multiple objective functions must be taken into account. One of the most relevant research topics that can benefit from these techniques is phylogenetic inference. Throughout the years, different researchers have proposed their own view about the reconstruction of ancestral evolutionary relationships among species. As a result, biologists often report different phylogenetic trees from a same dataset when considering distinct optimality principles. In this work, we detail a multiobjective swarm intelligence approach based on the novel Artificial Bee Colony algorithm for inferring phylogenies. The aim of this paper is to propose a complementary view of phylogenetics according to the maximum parsimony and maximum likelihood criteria, in order to generate a set of phylogenetic trees that represent a compromise between these principles. Experimental results on a variety of nucleotide data sets and statistical studies highlight the relevance of the proposal with regard to other multiobjective algorithms and state-of-the-art biological methods.     

Reconstruction of ancestral sequences by the inferential method,a tool for protein engineering studies

This paper describes the inferential method, an approach for reconstructing protein and nucleotide sequences of ancestral species, starting from known, homologous, contemporary sequences. The method requires knowledge of the topology of the phylogenetic tree, whose nodes are the species to whom the reconstructed sequences belong.The method has been tested by computer simulation of speciation and nucleotide substitutions, starting from a single ancestral sequence, and by subsequent reconstruction of nodal sequences. Results have shown that reconstructions obtained by the inferential method are affected by limited error frequencies, which (1) are proportional to the squares of nucleotide substitution rates and of internodal distances, and (2) are little influenced by non-uniformity of transformation rates of nucleotides.Furthermore, good agreement of the results has been obtained by comparing protein-sequence reconstructions carried out with the inferential method with those obtained using the maximum parsimony method in two different cases: e.g., a reconstruction of simulated sequences and a reconstruction of mammalian ribonuclease sequences.Abbreviations used MP maximum parsimony method - ML maximum likelihood method - IM inferential method - MY millions of years - N-tree natural-like phylogenetic tree - E-tree equibranched phylogenetic tree - EA percentage number of erroneous amino acids in a reconstructed sequence - EC percentage number of erroneous codons in a reconstructed sequence - t _n time interval between a P- and its - F-sequence nucleotides and amino acids are indicated by their I.U.B. codes (N.C.-I.U.B., 1985) Correspondence to: A. Di Donato     

Reanalysis of Murphy et al.'s data gives various mammalian phylogenies and suggests overcredibility of Bayesian trees

Murphy and colleagues reported that the mammalian phylogeny was resolved by Bayesian phylogenetics. However, the DNA sequences they used had many alignment gaps and undetermined nucleotide sites. We therefore reanalyzed their data by minimizing unshared nucleotide sites and retaining as many species as possible (13 species). In constructing phylogenetic trees, we used the Bayesian, maximum likelihood (ML), maximum parsimony (MP), and neighbor-joining (NJ) methods with different substitution models. These trees were constructed by using both protein and DNA sequences. The results showed that the posterior probabilities for Bayesian trees were generally much higher than the bootstrap values for ML, MP, and NJ trees. Two different Bayesian topologies for the same set of species were sometimes supported by high posterior probabilities, implying that two different topologies can be judged to be correct by Bayesian phylogenetics. This suggests that the posterior probability in Bayesian analysis can be excessively high as an indication of statistical confidence and therefore Murphy et al.'s tree, which largely depends on Bayesian posterior probability, may not be correct.     

Inconsistency of phylogenetic estimates from concatenated data under coalescence

Although multiple gene sequences are becoming increasingly available for molecular phylogenetic inference, the analysis of such data has largely relied on inference methods designed for single genes. One of the common approaches to analyzing data from multiple genes is concatenation of the individual gene data to form a single supergene to which traditional phylogenetic inference procedures - e.g., maximum parsimony (MP) or maximum likelihood (ML) - are applied. Recent empirical studies have demonstrated that concatenation of sequences from multiple genes prior to phylogenetic analysis often results in inference of a single, well-supported phylogeny. Theoretical work, however, has shown that the coalescent can produce substantial variation in single-gene histories. Using simulation, we combine these ideas to examine the performance of the concatenation approach under conditions in which the coalescent produces a high level of discord among individual gene trees and show that it leads to statistically inconsistent estimation in this setting. Furthermore, use of the bootstrap to measure support for the inferred phylogeny can result in moderate to strong support for an incorrect tree under these conditions. These results highlight the importance of incorporating variation in gene histories into multilocus phylogenetics.     

New approaches to phylogenetic tree search and their application to large numbers of protein alignments

Phylogenetic tree estimation plays a critical role in a wide variety of molecular studies, including molecular systematics, phylogenetics, and comparative genomics. Finding the optimal tree relating a set of sequences using score-based (optimality criterion) methods, such as maximum likelihood and maximum parsimony, may require all possible trees to be considered, which is not feasible even for modest numbers of sequences. In practice, trees are estimated using heuristics that represent a trade-off between topological accuracy and speed. I present a series of novel algorithms suitable for score-based phylogenetic tree reconstruction that demonstrably improve the accuracy of tree estimates while maintaining high computational speeds. The heuristics function by allowing the efficient exploration of large numbers of trees through novel hill-climbing and resampling strategies. These heuristics, and other computational approximations, are implemented for maximum likelihood estimation of trees in the program Leaphy, and its performance is compared to other popular phylogenetic programs. Trees are estimated from 4059 different protein alignments using a selection of phylogenetic programs and the likelihoods of the tree estimates are compared. Trees estimated using Leaphy are found to have equal to or better likelihoods than trees estimated using other phylogenetic programs in 4004 (98.6%) families and provide a unique best tree that no other program found in 1102 (27.1%) families. The improvement is particularly marked for larger families (80 to 100 sequences), where Leaphy finds a unique best tree in 81.7% of families.     

LisBeth: New cladistics for phylogenetics and biogeography

Within phylogenetics, two methods are known to implement cladistics: parsimony or maximum parsimony (MP) and three-item analysis (3ia). Despite the lack of suitable software, 3ia is occasionally used in systematic, and more regularly, in historical biogeography. Here, we present LisBeth, the first and only phylogenetic/biogeographic program freely available that uses the 3ia approach and offer some insights into its theoretical propositions. LisBeth does not rely on the conventional taxon/character matrix. Instead, characters are represented as rooted trees. LisBeth performs 3ia analyses based on maximum congruence of three-item statements and calculates the intersection tree (which differs from usual consensus). In biogeography, it applies the transparent method to handle widespread taxa and implements paralogy-free subtree analysis to remove redundant distributions. For the sake of interoperability, LisBeth may import/export characters from/to matrix in NEXUS format, allowing comparison with other cladistic programs. LisBeth also imports phylogenetic characters from Xper² knowledge bases.     

No speed dating please! Patterns of social preference in male and female house mice

Glass sponges (Class Hexactinellida) are important components of deep-sea ecosystems and are of interest from geological and materials science perspectives. The reconstruction of their phylogeny with molecular data has only recently begun and shows a better agreement with morphology-based systematics than is typical for other sponge groups, likely because of a greater number of informative morphological characters. However, inconsistencies remain that have far-reaching implications for hypotheses about the evolution of their major skeletal construction types (body plans). Furthermore, less than half of all described extant genera have been sampled for molecular systematics, and several taxa important for understanding skeletal evolution are still missing. Increased taxon sampling for molecular phylogenetics of this group is therefore urgently needed. However, due to their remote habitat and often poorly preserved museum material, sequencing all 126 currently recognized extant genera will be difficult to achieve. Utilizing morphological data to incorporate unsequenced taxa into an integrative systematics framework therefore holds great promise, but it is unclear which methodological approach best suits this task. Here, we increase the taxon sampling of four previously established molecular markers (18S, 28S, and 16S ribosomal DNA, as well as cytochrome oxidase subunit I) by 12 genera, for the first time including representatives of the order Aulocalycoida and the type genus of Dactylocalycidae, taxa that are key to understanding hexactinellid body plan evolution. Phylogenetic analyses suggest that Aulocalycoida is diphyletic and provide further support for the paraphyly of order Hexactinosida; hence these orders are abolished from the Linnean classification. We further assembled morphological character matrices to integrate so far unsequenced genera into phylogenetic analyses in maximum parsimony (MP), maximum likelihood (ML), Bayesian, and morphology-based binning frameworks. We find that of these four approaches, total-evidence analysis using MP gave the most plausible results concerning congruence with existing phylogenetic and taxonomic hypotheses, whereas the other methods, especially ML and binning, performed more poorly. We use our total-evidence phylogeny of all extant glass sponge genera for ancestral state reconstruction of morphological characters in MP and ML frameworks, gaining new insights into the evolution of major hexactinellid body plans and other characters such as different spicule types. Our study demonstrates how a comprehensive, albeit in some parts provisional, phylogeny of a larger taxon can be achieved with an integrative approach utilizing molecular and morphological data, and how this can be used as a basis for understanding phenotypic evolution. The datasets and associated trees presented here are intended as a resource and starting point for future work on glass sponge evolution.     

基于COI基因序列的九种真蝽属昆虫的分子系统学研究(半翅目:蝽科)

通过对真蝽属Pentatoma 9种昆虫线粒体COI基因约798bp的序列进行分子进化分析,并以同蝽科宽铗同蝽Acanthosoma labiduroides为外群,采用最大简约法、最大似然法和邻接法构建了分子系统树,来探讨真蝽属的系统发育关系.研究结果支持褐真蝽群P. semiannulata-group的划分,绿角真蝽Pentatoma viridicornuta应划归到褐真蝽群P. Semiannulata-group;红足真蝽群中的角肩真蝽P. angulata与红足真蝽P. rufipes遗传距离较小,它们是否为1个物种值得关注;真蝽属各群之间的系统发育关系以及是否可分为3个属值得进一步研究.     

Avian higher-level phylogeny: well-supported clades and what we can learn from a phylogenetic analysis of 2954 morphological characters

It has been shown that increased character sampling betters the accuracy of phylogenetic reconstructions in the case of molecular data. A recently published analysis of avian higher-level phylogenetics based on 2954 morphological characters now provides an empirical example to test whether this is also true in the case of morphological characters. Several clades are discussed which are supported by multiple analyses of mutually independent molecular data (sequences of nuclear genes on different chromosomes and mitochondrial genes) as well as morphological apomorphies, but did not result from parsimony analysis of the large morphological data set. Incorrect character scorings in that analysis notwithstanding, it is concluded that in the case of morphological data, increased character sampling does not necessarily better the accuracy of a phylogenetic reconstruction. Because morphological characters usually have a strongly varying complexity, many simple and homoplastic characters may overrule fewer ones of greater phylogenetic significance in large data sets, thus producing a low ratio of phylogenetic signal to 'noise' in the data.