A likelihood look at the supermatrix-supertree controversy

Supermatrix and supertree methods are two strategies advocated for phylogenetic analysis of sequence data from multiple gene loci, especially when some species are missing at some loci. The supermatrix method concatenates sequences from multiple genes into a data supermatrix for phylogenetic analysis, and ignores differences in evolutionary dynamics among the genes. The supertree method analyzes each gene separately and assembles the subtrees estimated from individual genes into a supertree for all species. Most algorithms suggested for supertree construction lack statistical justifications and ignore uncertainties in the subtrees. Instead of supermatrix or supertree, we advocate the use of likelihood function to combine data from multiple genes while accommodating their differences in the evolutionary process. This combines the strengths of the supermatrix and supertree methods while avoiding their drawbacks. We conduct computer simulation to evaluate the performance of the supermatrix, supertree, and maximum likelihood methods applied to two phylogenetic problems: molecular-clock dating of species divergences and reconstruction of species phylogenies. The results confirm the theoretical superiority of the likelihood method. Supertree or separate analyses of data of multiple genes may be useful in revealing the characteristics of the evolutionary process of multiple gene loci, and the information may be used to formulate realistic models for combined analysis of all genes by likelihood.     

Development of the cranium and paired fins in Betta splendens (Teleostei: Percomorpha): Intraspecific variation and interspecific comparisons

The development of the chondrocranium and the relative timing of ossification of the osteocranium is described in the teleost fish Betta splendens from a large series of cleared and differentially stained specimens. General trends in ossification patterns are examined from developmental, phylogenetic, and functional contexts. As in many other vertebrates, dermal bones form before cartilage bones. Ossification sequence conforms to functional need in a very general way, but there are many inconsistencies in the details of order. For example, some bones that are directly involved in feeding ossify no earlier than bones more indirectly involved. Comparisons of ossification sequence within specific cranial regions are made among Betta splendens, Oryzias latipes (Atherinomorpha), and Barbus barbus (Ostariophysi) within a phylogenetic framework. Many evolutionary changes in relative sequence of ossification are evident within regions among these taxa, yet many other sequences are conserved. The logistic difficulty of comparing entire cranial ossification sequences (vs. regional sequences) makes evident the need for new methods for identifying and quantifying sequence changes. Intraspecific variation in order of ossification is described for the first time in teleost fishes. To the extent that ossification sequence varies intraspecifically, conclusions drawn from previous interspecific comparisons are compromised. Understanding the importance of changes in ossification order within and among taxa will require experimental, functional, and evolutionary work. © 1996 Wiley-Liss, Inc.     

Intraphylum diversity and complex evolution of cyanobacterial aminoacyl-tRNA synthetases

A comparative genomic analysis of 35 cyanobacterial strains has revealed that the gene complement of aminoacyl-tRNA synthetases (AARSs) and routes for aminoacyl-tRNA synthesis may differ among the species of this phylum. Several genes encoding AARS paralogues were identified in some genomes. In-depth phylogenetic analysis was done for each of these proteins to gain insight into their evolutionary history. GluRS, HisRS, ArgRS, ThrRS, CysRS, and Glu-Q-RS showed evidence of a complex evolutionary course as indicated by a number of inconsistencies with our reference tree for cyanobacterial phylogeny. In addition to sequence data, support for evolutionary hypotheses involving horizontal gene transfer or gene duplication events was obtained from other observations including biased sequence conservation, the presence of indels (insertions or deletions), or vestigial traces of ancestral redundant genes. We present evidences for a novel protein domain with two putative transmembrane helices recruited independently by distinct AARS in particular cyanobacteria.     

Predicting loss of evolutionary history: Where are we?

The Earth's evolutionary history is threatened by species loss in the current sixth mass extinction event in Earth's history. Such extinction events not only eliminate species but also their unique evolutionary histories. Here we review the expected loss of Earth's evolutionary history quantified by phylogenetic diversity (PD) and evolutionary distinctiveness (ED) at risk. Due to the general paucity of data, global evolutionary history losses have been predicted for only a few groups, such as mammals, birds, amphibians, plants, corals and fishes. Among these groups, there is now empirical support that extinction threats are clustered on the phylogeny; however this is not always a sufficient condition to cause higher loss of phylogenetic diversity in comparison to a scenario of random extinctions. Extinctions of the most evolutionarily distinct species and the shape of phylogenetic trees are additional factors that can elevate losses of evolutionary history. Consequently, impacts of species extinctions differ among groups and regions, and even if global losses are low within large groups, losses can be high among subgroups or within some regions. Further, we show that PD and ED are poorly protected by current conservation practices. While evolutionary history can be indirectly protected by current conservation schemes, optimizing its preservation requires integrating phylogenetic indices with those that capture rarity and extinction risk. Measures based on PD and ED could bring solutions to conservation issues, however they are still rarely used in practice, probably because the reasons to protect evolutionary history are not clear for practitioners or due to a lack of data. However, important advances have been made in the availability of phylogenetic trees and methods for their construction, as well as assessments of extinction risk. Some challenges remain, and looking forward, research should prioritize the assessment of expected PD and ED loss for more taxonomic groups and test the assumption that preserving ED and PD also protects rare species and ecosystem services. Such research will be useful to inform and guide the conservation of Earth's biodiversity and the services it provides.     

Inferring phylogenetic relationships avoiding forbidden rooted triplets

To construct a phylogenetic tree or phylogenetic network for describing the evolutionary history of a set of species is a well-studied problem in computational biology. One previously proposed method to infer a phylogenetic tree/network for a large set of species is by merging a collection of known smaller phylogenetic trees on overlapping sets of species so that no (or as little as possible) branching information is lost. However, little work has been done so far on inferring a phylogenetic tree/network from a specified set of trees when in addition, certain evolutionary relationships among the species are known to be highly unlikely. In this paper, we consider the problem of constructing a phylogenetic tree/network which is consistent with all of the rooted triplets in a given set C and none of the rooted triplets in another given set F. Although NP-hard in the general case, we provide some efficient exact and approximation algorithms for a number of biologically meaningful variants of the problem.     

昆虫卵黄蛋白分子进化的研究进展

卵黄原蛋白（Vg）、卵黄多肽（YP）和小卵黄蛋白（minor YP）是昆虫三类主要的卵黄蛋白,它们之间的同源性一直是研究的重点。本文根据已经解析的Vg,YP和minor YP的氨基酸序列,采用序列比对和系统树分析的方法,并结合国内外对三者同源性研究的基础,对其进化关系进行了分析。结果表明,Vg,YP和minor YP是三类具有不同进化祖先的卵黄蛋白,它们的氨基酸序列相似性较低。Vg在系统进化过程中最为保守,与人类的血清载脂蛋白B(ApoB)具有较高的同源性;YP与脊椎动物的肝脂酶和胰脂酶具有较高的同源性;而minor YP与脊椎动物胃脂肪酶和舌脂肪酶具有较高的同源性。同时,对三者的分子特性做了简单的介绍。     

The colugo (Cynocephalus variegatus,Dermoptera): the primates' gliding sister?

Although a general agreement on the major groups of eutherian orders and their phylogenetic affiliations is emerging, the evolutionary affiliations among the members constituting these groups are still subject to debate. A prominent example is the recently published molecular evidence that challenges the long assumed monophyly of primates, displaying the colugo or flying lemur (Cynocephalus, Dermoptera) as a sister to anthropoid primates (Arnason et al. 2002 ) and positioning them after the prosimian primates (tarsiers and strepsirhines) split off. The phylogenetic analysis of the complete mitochondrial (mt) genome sequence of Cynocephalus variegatus presented in this study first appears to corroborate interpretations of primates as a paraphyletic group. However, more detailed analyses disclosed that mt nucleotide composition and consequently amino acid (AA) composition varied considerably among the species analyzed. This led us to assume that the flying lemur may be incorrectly grouped with anthropoids on the basis of similar mt nucleotide and AA compositions, rather than reflecting the true evolutionary relationship. To reanalyze the flying lemur's evolutionary association with other eutherian orders from a completely different molecular perspective, a molecular cladistic approach was applied. To this end, we determined the presence/absence pattern of transposable elements that provide a nearly homoplasy-free and copious source of molecular evolutionary markers, with well-defined character polarity. We could identify transposable elements, both on a multilocus and single-locus level, being present in all extant primate infraorders but absent in the flying lemur, thus clearly supporting the monophyly of primates by retropositional evidence.     

General properties and phylogenetic utilities of nuclear ribosomal DNA and mitochondrial DNA commonly used in molecular systematics

To choose one or more appropriate molecular markers or gene regions for resolving a particular systematic question among the organisms at a certain categorical level is still a very difficult process. The primary goal of this review, therefore, is to provide a theoretical information in choosing one or more molecular markers or gene regions by illustrating general properties and phylogenetic utilities of nuclear ribosomal DNA (rDNA) and mitochondrial DNA (mtDNA) that have been most commonly used for phylogenetic researches. The highly conserved molecular markers and/or gene regions are useful for investigating phylogenetic relationships at higher categorical levels (deep branches of evolutionary history). On the other hand, the hypervariable molecular markers and/or gene regions are useful for elucidating phylogenetic relationships at lower categorical levels (recently diverged branches). In summary, different selective forces have led to the evolution of various molecular markers or gene regions with varying degrees of sequence conservation. Thus, appropriate molecular markers or gene regions should be chosen with even greater caution to deduce true phylogenetic relationships over a broad taxonomic spectrum.     

介形类系统分类与起源演化的形态、分子和化石证据

介形类(Ostracoda)因其丰富的化石记录和广布的海陆现生代表类群,而被认为是进化生物学中研究生物多样性产生机制和演变历程的颇具潜力的重要模式生物。介形类在甲壳亚门中的谱系发生位置、起源及其内部各类群间的系统关系还存在诸多争议。基于其体制构造的形态学特征,介形类被归入甲壳亚门下的颚足纲(Maxillopoda),但来自18S rDNA序列数据分析却显示Maxillopoda不是单系群。基于化石记录和壳体形态特征,高肌虫(Bradoriida)长期以来被认为是介形类的一个祖先类群,但保存有软躯体的早寒武世化石的研究表明,Bradoriida不是介形类甚至可能也不属于甲壳类。不同的研究者所强调的壳体和肢体形态特征各异,导致介形类最大的现生类群速足目(Podocopida)的四个超科之间的关系也存在诸多推测。壳体和肢体特征在系统演化意义上的不兼容,需要分子生物学等证据的介入。分子、形态和化石证据的积累及各种信息整合是系统演化研究的必然趋势。     

Molecular function prediction for a family exhibiting evolutionary tendencies toward substrate specificity swapping: Recurrence of tyrosine aminotransferase activity in the Iα subfamily

The subfamily Iα aminotransferases are typically categorized as having narrow specificity toward carboxylic amino acids (AATases), or broad specificity that includes aromatic amino acid substrates (TATases). Because of their general role in central metabolism and, more specifically, their association with liver‐related diseases in humans, this subfamily is biologically interesting. The substrate specificities for only a few members of this subfamily have been reported, and the reliable prediction of substrate specificity from protein sequence has remained elusive. In this study, a diverse set of aminotransferases was chosen for characterization based on a scoring system that measures the sequence divergence of the active site. The enzymes that were experimentally characterized include both narrow‐specificity AATases and broad‐specificity TATases, as well as AATases with broader‐specificity and TATases with narrower‐specificity than the previously known family members. Molecular function and phylogenetic analyses underscored the complexity of this family's evolution as the TATase function does not follow a single evolutionary thread, but rather appears independently multiple times during the evolution of the subfamily. The additional functional characterizations described in this article, alongside a detailed sequence and phylogenetic analysis, provide some novel clues to understanding the evolutionary mechanisms at work in this family. Proteins 2013. © 2013 Wiley Periodicals, Inc.     

Anchor-based whole genome phylogeny (ABWGP): a tool for inferring evolutionary relationship among closely related microorganisms [corrected

Phenotypic behavior of a group of organisms can be studied using a range of molecular evolutionary tools that help to determine evolutionary relationships. Traditionally a gene or a set of gene sequences was used for generating phylogenetic trees. Incomplete evolutionary information in few selected genes causes problems in phylogenetic tree construction. Whole genomes are used as remedy. Now, the task is to identify the suitable parameters to extract the hidden information from whole genome sequences that truly represent evolutionary information. In this study we explored a random anchor (a stretch of 100 nucleotides) based approach (ABWGP) for finding distance between any two genomes, and used the distance estimates to compute evolutionary trees. A number of strains and species of Mycobacteria were used for this study. Anchor-derived parameters, such as cumulative normalized score, anchor order and indels were computed in a pair-wise manner, and the scores were used to compute distance/phylogenetic trees. The strength of branching was determined by bootstrap analysis. The terminal branches are clearly discernable using the distance estimates described here. In general, different measures gave similar trees except the trees based on indels. Overall the tree topology reflected the known biology of the organisms. This was also true for different strains of Escherichia coli. A new whole genome-based approach has been described here for studying evolutionary relationships among bacterial strains and species.     

Phylogenetic information and experimental design in molecular systematics.

Despite the widespread perception that evolutionary inference from molecular sequences is a statistical problem, there has been very little attention paid to questions of experimental design. Previous consideration of this topic has led to little more than an empirical folklore regarding the choice of suitable genes for analysis, and to dispute over the best choice of taxa for inclusion in data sets. I introduce what I believe are new methods that permit the quantification of phylogenetic information in a sequence alignment. The methods use likelihood calculations based on Markov-process models of nucleotide substitution allied with phylogenetic trees, and allow a general approach to optimal experimental design. Two examples are given, illustrating realistic problems in experimental design in molecular phylogenetics and suggesting more general conclusions about the choice of genomic regions, sequence lengths and taxa for evolutionary studies.     

Anchor-Based Whole Genome Phylogeny (ABWGP): A Tool for Inferring Evolutionary Relationship among Closely Related Microorganims

Phenotypic behavior of a group of organisms can be studied using a range of molecular evolutionary tools that help to determine evolutionary relationships. Traditionally a gene or a set of gene sequences was used for generating phylogenetic trees. Incomplete evolutionary information in few selected genes causes problems in phylogenetic tree construction. Whole genomes are used as remedy. Now, the task is to identify the suitable parameters to extract the hidden information from whole genome sequences that truly represent evolutionary information. In this study we explored a random anchor (a stretch of 100 nucleotides) based approach (ABWGP) for finding distance between any two genomes, and used the distance estimates to compute evolutionary trees. A number of strains and species of Mycobacteria were used for this study. Anchor-derived parameters, such as cumulative normalized score, anchor order and indels were computed in a pair-wise manner, and the scores were used to compute distance/phylogenetic trees. The strength of branching was determined by bootstrap analysis. The terminal branches are clearly discernable using the distance estimates described here. In general, different measures gave similar trees except the trees based on indels. Overall the tree topology reflected the known biology of the organisms. This was also true for different strains of Escherichia coli. A new whole genome-based approach has been described here for studying evolutionary relationships among bacterial strains and species.     

Utility of nuclear allele networks for the analysis of closely related species in the genus Carabus, subgenus Ohomopterus

Nuclear DNA sequence data for diploid organisms are potentially a rich source of phylogenetic information for disentangling the evolutionary relationships of closely related organisms, but present special phylogenetic problems owing to difficulties arising from heterozygosity and recombination. We analyzed allelic relationships for two nuclear gene regions (phosphoenolpyruvate carboxykinase and elongation factor-1a), along with a mitochondrial gene region (NADH dehydrogenase subunit 5), for an assemblage of closely related species of carabid beetles (Carabus subgenus Ohomopterus). We used a network approach to examine whether the nuclear gene sequences provide substantial phylogenetic information on species relationships and evolutionary history. The mitochondrial gene genealogy strongly contradicted the morphological species boundary as a result of introgression of heterospecific mitochondria. Two nuclear gene regions showed high allelic diversity within species, and this diversity was partially attributable to recombination between various alleles and high variability in the intron region. Shared nuclear alleles among species were rare and were considered to represent shared ancestral polymorphism. Despite the presence of recombination, nuclear allelic networks recovered species monophyly more often and presented genetic differentiation patterns (low to high) among species more clearly. Overall, nuclear gene networks provide clear evidence for separate biological species and information on the phylogenetic relationships among closely related carabid beetles.     

Ribosomal ITS sequences and plant phylogenetic inference

One of the most popular sequences for phylogenetic inference at the generic and infrageneric levels in plants is the internal transcribed spacer (ITS) region of the 18S-5.8S-26S nuclear ribosomal cistron. The prominence of this source of nuclear DNA sequence data is underscored by a survey of phylogenetic publications involving comparisons at the genus level or below, which reveals that of 244 papers published over the last five years, 66% included ITS sequence data. Perhaps even more striking is the fact that 34% of all published phylogenetic hypothesis have been based exclusively on ITS sequences. Notwithstanding the many important contributions of ITS sequence data to phylogenetic understanding and knowledge of genome relationships, a number of molecular genetic processes impact ITS sequences in ways that may mislead phylogenetic inference. These molecular genetic processes are reviewed here, drawing attention to both underlying mechanism and phylogenetic implications. Among the most prevalent complications for phylogenetic inference is the existence in many plant genomes of extensive sequence variation, arising from ancient or recent array duplication events, genomic harboring of pseudogenes in various states of decay, and/or incomplete intra- or inter-array homogenization. These phenomena separately and collectively create a network of paralogous sequence relationships potentially confounding accurate phylogenetic reconstruction. Homoplasy is shown to be higher in ITS than in other DNA sequence data sets, most likely because of orthology/paralogy conflation, compensatory base changes, problems in alignment due to indel accumulation, sequencing errors, or some combination of these phenomena. Despite the near-universal usage of ITS sequence data in plant phylogenetic studies, its complex and unpredictable evolutionary behavior reduce its utility for phylogenetic analysis. It is suggested that more robust insights are likely to emerge from the use of single-copy or low-copy nuclear genes.     

Evidence for phylogenetic congruence among sulfur-oxidizing chemoautotrophic bacterial endosymbionts and their bivalve hosts

Sulfur-oxidizing chemoautotrophic (thioautotrophic) bacteria are now known to occur as endosymbionts in phylogenetically diverse bivalve hosts found in a wide variety of marine environments. The evolutionary origins of these symbioses, however, have remained obscure. Comparative 16S rRNA sequence analysis was used to investigate whether thioautotrophic endosymbionts are monophyletic or polyphyletic in origin and to assess whether phylogenetic relationships inferred among these symbionts reflect those inferred among their hosts. 16S rRNA gene sequences determined for endosymbionts from nine newly examined bivalve species from three families (Vesicomyidae, Lucinidae, and Solemyidae) were compared with previously published 16S rRNA sequences of thioautotrophic symbionts and free-living bacteria. Distance and parsimony methods were used to infer phylogenetic relationships among these bacteria. All newly examined symbionts fall within the gamma subdivision of the Proteobacteria, in clusters containing previously examined symbiotic thioautotrophs. The closest free-living relatives of these symbionts are bacteria of the genus Thiomicrospira. Symbionts of the bivalve superfamily Lucinacea and the family Vesicomyidae each form distinct monophyletic lineages which are strongly supported by bootstrap analysis, demonstrating that host phylogenies inferred from morphological and fossil evidence are congruent with phylogenies inferred for their respective symbionts by molecular sequence analysis. The observed congruence between host and symbiont phylogenies indicates shared evolutionary history of hosts and symbiont lineages and suggests an ancient origin for these symbioses. Correspondence to: D.L. Distel     

MORPHOMETRICS AND CLADISTICS: MEASURING PHYLOGENY IN THE SEA URCHIN ECHINOCARDIUM

A phylogenetic approach to the study of evolutionary patterns is based on taxic homologies (synapomorphies). In contrast, the recognition of evolutionary processes (namely heterochronies) involves analysis of the entire morphology. Recent developments in geometric morphometry permit analysis of morphological similarities grounded in operational homologies. Such morphometric techniques are explored (1) at the level of evolutionary processes, and (2) as a complement in exploration of phylogenetic relationships. To examplify this, we perform a two-part study of the ontogeny and phylogeny of the spatangoid sea urchin Echinocardium. First, a phylogenetic analysis of ten Recent species in the genus is performed on 18 informative characters of the test. Second, morphological divergences among the species are analyzed using procrustean (superimposition) methods based on 49 homologous points. An additive distance tree is built from a matrix of morphometric distances among adult specimens. This tree is fully congruent with the phyletic results. Ontogenetic processes are explored by inserting ontogenetic series into the analysis. A distance tree including the juvenile stages shows that the general evolutionary trend of the genus is peramorphic, but species-to-species comparisons attest that no general clinal trend exists. Our analysis emphasizes the importance of morphometric approaches in evolutionary studies (1) for the understanding of heterochronies; (2) to trace the morphological implications of phylogenetic patterns; and (3) to estimate the impact of homoplasies.     

Phylogenetic relationships among Staphylococcus species and refinement of cluster groups based on multilocus data

ABSTRACT: BACKGROUND: Estimates of relationships among Staphylococcus species have been hampered by poor and inconsistent resolution of phylogenies based largely on single gene analyses incorporating only a limited taxon sample. As such, the evolutionary relationships and hierarchical classification schemes among species have not been confidently established. Here, we address these points through analyses of DNA sequence data from multiple loci (16S rRNA gene, dnaJ, rpoB, and tuf gene fragments) using multiple Bayesian and maximum likelihood phylogenetic approaches that incorporate nearly all recognized Staphylococcus taxa. RESULTS: We estimated the phylogeny of fifty-seven Staphylococcus taxa using partitioned-model Bayesian and maximum likelihood analysis, as well as Bayesian gene-tree species-tree methods. Regardless of methodology, we found broad agreement among methods that the current cluster groups require revision, although there was some disagreement among methods in resolution of higher order relationships. Based on our phylogenetic estimates, we propose a refined classification for Staphylococcus with species being classified into 15 cluster groups (based on molecular data) that adhere to six species groups (based on phenotypic properties) CONCLUSIONS: Our findings are in general agreement with gene tree-based reports of the staphylococcal phylogeny, although we identify multiple previously unreported relationships among species. Our results support the general importance of such multilocus assessments as a standard in microbial studies to more robustly infer relationships among recognized and newly discovered.     

Phylogenetic meta-analysis

Meta-analysis is a powerful statistical technique that combines the results of independent studies to identify general trends. When the species under examination are not independent however, it is also necessary to incorporate phylogenetic information into the analysis. Unfortunately, current meta-analytic approaches cannot account for lack of independence resulting from shared evolutionary history, so a general solution to this problem is lacking. In this article, I derive a model for phylogenetic meta-analysis, so that data across studies may be summarized with evolutionary history explicitly incorporated. The approach takes advantage of common aspects of linear statistical models used by both meta-analysis and the phylogenetic comparative method, thereby allowing them to be analytically combined. In this manner, the correlation structure generated by phylogenetic history can be incorporated directly into the meta-analytic procedure. I illustrate the approach by examining the prevalence of body size clines in mammals. The approach is general, and can also be used to incorporate correlation structure among studies generated by other factors, such as spatial or temporal proximity, or environmental similarity. Therefore, this procedure provides a general statistical template for meta-analytic techniques that can account for attributes that generate nonindependence among studies. Implications of the phylogenetic meta-analysis are discussed.     

Analysis of plasmid genes by phylogenetic profiling and visualization of homology relationships using Blast2Network

Background  

Phylogenetic methods are well-established bioinformatic tools for sequence analysis, allowing to describe the non-independencies of sequences because of their common ancestor. However, the evolutionary profiles of bacterial genes are often complicated by hidden paralogy and extensive and/or (multiple) horizontal gene transfer (HGT) events which make bifurcating trees often inappropriate. In this context, plasmid sequences are paradigms of network-like relationships characterizing the evolution of prokaryotes. Actually, they can be transferred among different organisms allowing the dissemination of novel functions, thus playing a pivotal role in prokaryotic evolution. However, the study of their evolutionary dynamics is complicated by the absence of universally shared genes, a prerequisite for phylogenetic analyses.