Assessing the Accuracy of Ancestral Protein Reconstruction Methods

The phylogenetic inference of ancestral protein sequences is a powerful technique for the study of molecular evolution, but any conclusions drawn from such studies are only as good as the accuracy of the reconstruction method. Every inference method leads to errors in the ancestral protein sequence, resulting in potentially misleading estimates of the ancestral protein's properties. To assess the accuracy of ancestral protein reconstruction methods, we performed computational population evolution simulations featuring near-neutral evolution under purifying selection, speciation, and divergence using an off-lattice protein model where fitness depends on the ability to be stable in a specified target structure. We were thus able to compare the thermodynamic properties of the true ancestral sequences with the properties of “ancestral sequences” inferred by maximum parsimony, maximum likelihood, and Bayesian methods. Surprisingly, we found that methods such as maximum parsimony and maximum likelihood that reconstruct a “best guess” amino acid at each position overestimate thermostability, while a Bayesian method that sometimes chooses less-probable residues from the posterior probability distribution does not. Maximum likelihood and maximum parsimony apparently tend to eliminate variants at a position that are slightly detrimental to structural stability simply because such detrimental variants are less frequent. Other properties of ancestral proteins might be similarly overestimated. This suggests that ancestral reconstruction studies require greater care to come to credible conclusions regarding functional evolution. Inferred functional patterns that mimic reconstruction bias should be reevaluated

Mitochondrial DNA clocks and the phylogeny of Danaus butterflies

Molecular clocks based on sequence change in mitochondrial (mt) DNA have been useful for placing molecular phytogenies in their historical context, thereby enhancing evolutionary insight. Nonetheless, despite their importance to phylogeographers, the methodology is controversial. Here we report on two mitochondrial clocks for the butterfly genus Danaus based on sequences from the cytochrome c oxidase subunit I (COI) and small subunit 12S rRNA (12S) genes. Both clocks are, within the context of Danaus, reliable time-keepers, mutually consistent and, respectively, in agreement with a crustacean COI clock and a molluscan 12S clock. Though we have no fossils with which directly to calibrate sequence divergence rates for Danaus, the 12S molluscan and COI crustacean clocks chosen for comparison were calibrated to radiometrically dated geomorphological events. Our results indicate that the Danaus COI clock evolves approximately four times faster than the 12S clock. Differences between rates of sequence change in terminal sister-taxa are small and likelihood ratio tests do not reject a hypothesis that evolution has been clock-like. The species Danaus chrysippus is paraphyletic and, therefore, invalid. Danaus probably split from its sister-genus Tirumala around 4.9 ± 0.3 million years ago in the early Pliocene

Quantum mechanics of lattice gas automata. I. One particle plane waves and potentials

Classical lattice gas automata effectively simulate physical processes such as diffusion and fluid flow (in certain parameter regimes) despite their simplicity at the microscale. Motivated by current interest in quantum computation we recently defined quantum lattice gas automata; in this paper we initiate a project to analyze which physical processes these models can effectively simulate. Studying the single particle sector of a one dimensional quantum lattice gas we find discrete analogues of plane waves and wave packets, and then investigate their behaviour in the presence of inhomogeneous potentials.Comment: 19 pages, plain TeX, 14 PostScript figures included with epsf.tex (ignore the under/overfull \vbox error messages), two additional large figures available upon reques

Lyapunov spectral analysis of a nonequilibrium Ising-like transition

By simulating a nonequilibrium coupled map lattice that undergoes an Ising-like phase transition, we show that the Lyapunov spectrum and related dynamical quantities such as the dimension correlation length~$\xi_\delta$ are insensitive to the onset of long-range ferromagnetic order. As a function of lattice coupling constant~$g$ and for certain lattice maps, the Lyapunov dimension density and other dynamical order parameters go through a minimum. The occurrence of this minimum as a function of~$g$ depends on the number of nearest neighbors of a lattice point but not on the lattice symmetry, on the lattice dimensionality or on the position of the Ising-like transition. In one-space dimension, the spatial correlation length associated with magnitude fluctuations and the length~$\xi_\delta$ are approximately equal, with both varying linearly with the radius of the lattice coupling.Comment: 29 pages of text plus 15 figures, uses REVTeX macros. Submitted to Phys. Rev. E

From quantum cellular automata to quantum lattice gases

A natural architecture for nanoscale quantum computation is that of a quantum cellular automaton. Motivated by this observation, in this paper we begin an investigation of exactly unitary cellular automata. After proving that there can be no nontrivial, homogeneous, local, unitary, scalar cellular automaton in one dimension, we weaken the homogeneity condition and show that there are nontrivial, exactly unitary, partitioning cellular automata. We find a one parameter family of evolution rules which are best interpreted as those for a one particle quantum automaton. This model is naturally reformulated as a two component cellular automaton which we demonstrate to limit to the Dirac equation. We describe two generalizations of this automaton, the second of which, to multiple interacting particles, is the correct definition of a quantum lattice gas.Comment: 22 pages, plain TeX, 9 PostScript figures included with epsf.tex (ignore the under/overfull \vbox error messages); minor typographical corrections and journal reference adde

Differential effects of glucagon-like peptide-1 receptor agonists on heart rate

Abstract While glucagon-like peptide-1 receptor agonists (GLP-1 RAs) are known to increase heart rate (HR), it is insufficiently recognized that the extent varies greatly between the various agonists and is affected by the assessment methods employed. Here we review published data from 24-h time-averaged HR monitoring in healthy individuals and subjects with type 2 diabetes mellitus (T2DM) treated with either short-acting GLP-1 RAs, lixisenatide or exenatide, or long-acting GLP-1 RAs, exenatide LAR, liraglutide, albiglutide, or dulaglutide (N\ua0=\ua01112; active-treatment arms). HR effects observed in two independent head-to-head trials of lixisenatide and liraglutide (N\ua0=\ua0202; active-treatment arms) are also reviewed. Short-acting GLP-1 RAs, exenatide and lixisenatide, are associated with a transient (1\u201312\ua0h) mean placebo- and baseline-adjusted 24-h HR increase of 1\u20133\ua0beats per minute (bpm). Conversely, long-acting GLP-1 RAs are associated with more pronounced increases in mean 24-h HR; the highest seen with liraglutide and albiglutide at 6\u201310\ua0bpm compared with dulaglutide and exenatide LAR at 3\u20134\ua0bpm. For both liraglutide and dulaglutide, HR increases were recorded during both the day and at night. In two head-to-head comparisons, a small, transient mean increase in HR from baseline was observed with lixisenatide; liraglutide induced a substantially greater increase that remained significantly elevated over 24\ua0h. The underlying mechanism for increased HR remains to be elucidated; however, it could be related to a direct effect at the sinus node and/or stimulation of the sympathetic nervous system, with this effect related to the duration of action of the respective GLP-1 RAs. In conclusion, this review indicates that the effects on HR differ within the class of GLP-1 RAs: short-acting GLP-1 RAs are associated with a modest and transient HR increase before returning to baseline levels, while some long-acting GLP-1 RAs are associated with a more pronounced and sustained increase during the day and night. Findings from recently completed trials indicate that a GLP-1 RA-induced increase in HR, regardless of magnitude, does not present an increased cardiovascular risk for subjects with T2DM, although a pronounced increase in HR may be associated with adverse clinical outcomes in those with advanced heart failure

A Comparison of Phylogenetic Network Methods Using Computer Simulation

Background: We present a series of simulation studies that explore the relative performance of several phylogenetic network approaches (statistical parsimony, split decomposition, union of maximum parsimony trees, neighbor-net, simulated history recombination upper bound, median-joining, reduced median joining and minimum spanning network) compared to standard tree approaches, (neighbor-joining and maximum parsimony) in the presence and absence of recombination. Principal Findings: In the absence of recombination, all methods recovered the correct topology and branch lengths nearly all of the time when the substitution rate was low, except for minimum spanning networks, which did considerably worse. At a higher substitution rate, maximum parsimony and union of maximum parsimony trees were the most accurate. With recombination, the ability to infer the correct topology was halved for all methods and no method could accurately estimate branch lengths. Conclusions: Our results highlight the need for more accurate phylogenetic network methods and the importance of detecting and accounting for recombination in phylogenetic studies. Furthermore, we provide useful information for choosing a network algorithm and a framework in which to evaluate improvements to existing methods and nove

From Gene Trees to Organismal Phylogeny in Prokaryotes:The Case of the γ-Proteobacteria

The rapid increase in published genomic sequences for bacteria presents the first opportunity to reconstruct evolutionary events on the scale of entire genomes. However, extensive lateral gene transfer (LGT) may thwart this goal by preventing the establishment of organismal relationships based on individual gene phylogenies. The group for which cases of LGT are most frequently documented and for which the greatest density of complete genome sequences is available is the γ-Proteobacteria, an ecologically diverse and ancient group including free-living species as well as pathogens and intracellular symbionts of plants and animals. We propose an approach to multigene phylogeny using complete genomes and apply it to the case of the γ-Proteobacteria. We first applied stringent criteria to identify a set of likely gene orthologs and then tested the compatibilities of the resulting protein alignments with several phylogenetic hypotheses. Our results demonstrate phylogenetic concordance among virtually all (203 of 205) of the selected gene families, with each of the exceptions consistent with a single LGT event. The concatenated sequences of the concordant families yield a fully resolved phylogeny. This topology also received strong support in analyses aimed at excluding effects of heterogeneity in nucleotide base composition across lineages. Our analysis indicates that single-copy orthologous genes are resistant to horizontal transfer, even in ancient bacterial groups subject to high rates of LGT. This gene set can be identified and used to yield robust hypotheses for organismal phylogenies, thus establishing a foundation for reconstructing the evolutionary transitions, such as gene transfer, that underlie diversity in genome content and organization