Bayesian Ancestral Reconstruction for Bat Echolocation

Total Page:16

File Type:pdf, Size:1020Kb

Bayesian Ancestral Reconstruction for Bat Echolocation Bayesian Ancestral Reconstruction for Bat Echolocation by Joseph Patrick Meagher Thesis Submitted to the University of Warwick for the degree of Doctor of Philosophy Statistics June 2020 Contents List of Tables iv List of Figures v Acknowledgments vii Declarations viii Abstract ix Abbreviations x Chapter 1 Introduction 1 Chapter 2 Literature Review 7 2.1 Some Background on Bats . .7 2.2 The Phylogenetic Comparative Method . 12 2.3 Statistical Models for Data . 16 2.3.1 Gaussian Processes . 16 2.3.2 The Phylogenetic Gaussian Process Framework . 20 2.3.3 Latent Variable Models and Factor Analysis . 24 2.3.4 Functional Data Analysis . 26 2.3.5 Acoustic Signal Processing . 30 2.4 Bayesian Inference . 33 2.4.1 MCMC Methods for Parameter Inference . 35 2.4.2 MCMC Methods for Gaussian Processes . 37 2.4.3 Estimating Model Evidence . 39 Chapter 3 A Phylogenetic Latent Variable Model for Function-valued Traits 41 i 3.1 Introduction . 41 3.2 Methods . 45 3.2.1 The Phylogeny: A Graphical Model for Shared Ancestry . 45 3.2.2 A Phylogenetic Latent Variable Model for Function-valued Traits . 47 3.2.3 Efficient Computation of the Model Likelihood . 48 3.2.4 Prior Specification . 49 3.2.5 Posterior Inference and Model Selection . 53 3.2.6 Ancestral Reconstruction . 56 3.3 Results for a Synthetic Example . 58 3.4 Discussion . 62 Chapter 4 A Generalised Phylogenetic Latent Variable Model 68 4.1 Introduction . 68 4.2 Methods . 70 4.2.1 Data Augmentation . 70 4.2.2 A Generalised Phylogenetic Latent Variable Model . 71 4.2.3 Approximate Posterior Inference . 74 4.2.4 Ancestral Reconstruction . 79 4.3 Results for a Synthetic Example . 80 4.3.1 Model Fitting . 81 4.3.2 Ancestral Reconstruction . 84 4.3.3 Parameter Inference . 86 4.4 Discussion . 87 Chapter 5 Ancestral Reconstruction of the Bat Echolocation Call 93 5.1 Introduction . 93 5.2 A Harmonic Model for Bat Echolocation . 96 5.2.1 Prior Specification . 98 5.2.2 Maximum-a-Posteriori Inference . 102 5.2.3 Fitting the Harmonic Model . 105 5.2.4 A Brief Discussion of the Harmonic Model . 106 5.3 Echolocation Call Reconstruction . 109 5.3.1 Echolocation Call Data . 109 5.3.2 Bat Phylogeny . 110 5.3.3 Echolocation Call Features . 111 5.3.4 Ancestral Reconstruction . 116 5.4 Discussion . 123 ii Chapter 6 Final Remarks 126 Appendix A Tree Traversal Algorithms 132 A.1 Pruned Likelihood Calculation . 132 A.2 Pruned Conditional Distribution . 136 Appendix B Derivations for Variational Inference 139 B.1 Co-ordinate Ascent Variational Inference Updates . 139 B.2 The Evidence Lower Bound . 146 B.3 Predictive Distribution . 151 Appendix C Alternative generalised PLVMs 153 iii List of Tables 3.1 Fixed Parameter and Hyper-Parameter Values for Synthetic Data . 58 3.2 Bayes Factors for Model Comparison . 60 4.1 Fixed Hyper-Parameter Values for Synthetic Phylogenetic Gaussian Processes . 81 5.1 Mexican Bat Echolocation Call Dataset . 110 5.2 MAP estimates and intervals of 90% posterior density for phyloge- netic hyper-parameters of the V-PLVM. 121 iv List of Figures 1.1 The Ancestral Bat Call . .6 2.1 Grouping Organisms over Phylogenies: The definition of monophyletic, paraphyletic, and polyphyletic groups . .9 2.2 The Diversity of Bat Echolocation Calls: Selected call spectrograms 10 2.3 Gaussian Process Regression: Illustration of prior and posterior sam- ples from a Gaussian process . 18 2.4 Comparison of Mat´ernkernels: Illustration of Mat´erncovariance functions and samples from the process for different values of the smoothing parameter ν ......................... 20 2.5 Positions on a Phylogeny . 21 2.6 Registration of Functional Data: A toy example . 29 2.7 Interpolation between Signal Characterisations: A comparison of in- stantaneous frequency and spectrogram representations . 34 3.1 A Taxon-level Phylogeny . 45 3.2 A Phylogeny for Repeated Measurements . 46 3.3 The Phylogenetic Latent Variable Model: A graphical representation 53 3.4 Loadings for a Synthetic Example . 59 3.5 Phylogeny and Trait Observations for a Synthetic Example . 59 3.6 Ancestral Reconstruction of a synthetic Function-Valued Trait . 61 3.7 Sampled Posterior Loading . 63 3.8 Sampled Posterior Phylogenetic Hyper-parameters . 64 3.9 Sampled Posterior Parameters . 65 4.1 The Generalised Phylogenetic Latent Variable Model: A graphical representation . 75 4.2 A Synthetic Collection of Traits on a Phylogeny . 82 4.3 The log Evidence Lower Bound for Generalised PLVMs . 83 v 4.4 Root Ancestral Distribution: The V-PLVM . 85 4.5 Inferred Loading: The V-PLVM . 88 4.6 Inferred Phylogenetic Hyper-parameters: The V-PLVM . 89 4.7 Root Ancestral Distribution: The R-PLVM . 92 5.1 A Harmonic Model for Bat Echolocation: A graphical representation 101 5.2 Fitted Harmonic Models: A selection of bat echolocation calls . 107 5.3 A Phylogeny for Sampled Mexican Bats . 112 5.4 Corrected Fundamental Frequency Curves: A selection of bat echolo- cation calls . 114 5.5 Time Registration of Fundamental Frequency Curves: Pteronotus parnellii ................................. 115 5.6 The log Evidence Lower Bound: Models for the evolution of bat echolocation calls . 119 5.7 Echolocation in Bats Most Recent Common Ancestor . 120 5.8 The Evolution of Bat Echolocation . 122 5.9 Inferred Loadings: A model for the evolution of bat echolocation . 124 A.1 A Toy Phylogeny . 132 C.1 Root Ancestral Distribution: The R-PLVM . 154 C.2 Root Ancestral Distribution: The P-PLVM . 155 C.3 Root Ancestral Distribution: The I-PLVM . 156 vi Acknowledgments Thank you to my supervisors on this project. To Mark Girolami, for offering me the opportunity to undertake this programme of doctoral research; to Kate Jones, for providing such a motivating problem, and to Theo Damoulas, whose mentorship made this work possible. Thank you to both the EPSRC for their generous funding of this project and the Department of Statistics at the University of Warwick for hosting me throughout. Thank you to my friends, old and new, who have been with me on this journey. In particular, I want to thank Arthur, five years is a long time to spend under the same roof, but you made it easy. Thank you to my family; to my grandparents, for whom my admiration only increases with the passage of time; to Ailshe and Liam, may you learn from the mistakes of your brother, whatever you may judge them to be; to my mother, whose constant self-development is as much motivation as it is inspiration; and to my father, for demonstrating the value of consistent enthusiasm. I can never truly thank you for all I've been given in this life, but only do my best to make the most of it. And to Suzy, on to the next chapter. vii Declarations This thesis is submitted to the University of Warwick in support of my application for the degree of Doctor of Philosophy. I declare that it has been composed by myself and that the work contained herein is my own except where explicitly stated otherwise. This work has been completed wholly while in candidature for a research degree at the University of Warwick and has not been submitted for any other degree or professional qualification. Aspects of this research have been published in: • JP Meagher, T Damoulas, KE Jones, and M Girolami. Discussion of \The statistical analysis of acoustic phonetic data: exploring differences between spoken Romance languages", by Davide Pigoli, Pantelis Z Hadjipantelis, John S Coleman, and John AD Aston. Journal of the Royal Statistical Society: Series C (Applied Statistics), 67(5):1103-1145, 2018. • JP Meagher, T Damoulas, KE Jones, and M Girolami. Phylogenetic Gaussian processes for bat echolocation. Statistical Data Science, pages 111-122, 2018. doi:10.1142/q0159. viii Abstract Ancestral reconstruction can be understood as an interpolation between mea- sured characteristics of existing populations to those of their common ancestors. Doing so provides an insight into the characteristics of organisms that lived mil- lions of years ago. Such reconstructions are inherently uncertain, making this an ideal application area for Bayesian statistics. As such, Gaussian processes serve as a basis for many probabilistic models for trait evolution, which assume that mea- sured characteristics, or some transformation of those characteristics, are jointly Gaussian distributed. While these models do provide a theoretical basis for un- certainty quantification in ancestral reconstruction, practical approaches to their implementation have proven challenging. In this thesis, novel Bayesian methods for ancestral reconstruction are developed and applied to bat echolocation calls. This work proposes the first fully Bayesian approach to inference within the Phylogenetic Gaussian Process Regression framework for Function-Valued Traits, producing an ancestral reconstruction for which any uncertainty in this model may be quantified. The framework is then generalised to collections of discrete and continuous traits, and an efficient approximate Bayesian inference scheme proposed, representing the first application of Variational inference techniques to the problem of ancestral re- construction. This efficient approach is then applied to the reconstruction of bat echolocation calls, providing new insights into the developmental pathways of this remarkable characteristic. It is the complexity of bat echolocation that motivates the proposed approach to evolutionary inference, however, the resulting statistical methods are broadly applicable within the field of Evolutionary Biology. ix Abbreviations AM Adaptive Metropolis ARD Automatic Relevance Determination ASIS Ancillarity-Sufficiency Interweaving Strategy
Recommended publications
  • Uncorrelated Genetic Drift of Gene Frequencies and Linkage
    Genet. Res., Camb. (1974), 24, pp. 281-294 281 Printed in Great Britain Uncorrelated genetic drift of gene frequencies and linkage disequilibrium in some models of linked overdominant polymorphisms BY JOSEPH FELSENSTEIN Department of Genetics, University of Washington, Seattle, Washington 98195 {Received 15 April 1974) SUMMARY For large population sizes, gene frequencies p and q at two linked over- dominant loci and the linkage disequilibrium parameter D will remain close to their equilibrium values. We can treat selection and recombination as approximately linear forces on^J, q and D, and we can treat genetic drift as a multivariate normal perturbation with constant variance-covariance matrix. For the additive-multiplicative family of two-locus models, p, q and D are shown to be (approximately) uncorrelated. Expressions for their variances are obtained. When selection coefficients are small the variances of p and q are those previously given by Robertson for a single locus. For small recombination fractions the variance of D is that obtained for neutral loci by Ohta & Kimura. For larger recombination fractions the result differs from theirs, so that for unlinked loci r2~ 2/(3N) instead of l/(2N). For the Lewontin-Kojima and Bodmer symmetric viability models, and for a model symmetric at only one of the loci, a more exact argument is possible. In the asymptotic conditional distribution in these cases, various of p, q and D are uncorrelated, depending on the type of symmetiy in the model. 1. INTRODUCTION Much of the work on linked genes in the last few years has centred on the deter- ministic theory of natural selection of linked polymorphisms.
    [Show full text]
  • Comparative Methods for Reconstructing Ancient Genome Organization Yoann Anselmetti, Nina Luhmann, Sèverine Bérard, Eric Tannier, Cedric Chauve
    Comparative Methods for Reconstructing Ancient Genome Organization Yoann Anselmetti, Nina Luhmann, Sèverine Bérard, Eric Tannier, Cedric Chauve To cite this version: Yoann Anselmetti, Nina Luhmann, Sèverine Bérard, Eric Tannier, Cedric Chauve. Comparative Methods for Reconstructing Ancient Genome Organization. Comparative Genomics: Methods and Protocols, pp.343 - 362, 2017, 10.1007/978-1-4939-7463-4_13. hal-03192460 HAL Id: hal-03192460 https://hal.archives-ouvertes.fr/hal-03192460 Submitted on 8 Apr 2021 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Chapter 13 Comparative Methods for Reconstructing Ancient Genome Organization Yoann Anselmetti, Nina Luhmann, Se`verine Be´rard, Eric Tannier, and Cedric Chauve Abstract Comparative genomics considers the detection of similarities and differences between extant genomes, and, based on more or less formalized hypotheses regarding the involved evolutionary processes, inferring ancestral states explaining the similarities and an evolutionary history explaining the differences. In this chapter, we focus on the reconstruction of the organization of ancient genomes into chromosomes. We review different methodological approaches and software, applied to a wide range of datasets from different kingdoms of life and at different evolutionary depths. We discuss relations with genome assembly, and potential approaches to validate computational predictions on ancient genomes that are almost always only accessible through these predictions.
    [Show full text]
  • Reconstruction of Ancestral RNA Sequences Under Multiple Structural Constraints Olivier Tremblay-Savard1,2, Vladimir Reinharz1 and Jérôme Waldispühl1*
    The Author(s) BMC Genomics 2016, 17(Suppl 10):862 DOI 10.1186/s12864-016-3105-4 RESEARCH Open Access Reconstruction of ancestral RNA sequences under multiple structural constraints Olivier Tremblay-Savard1,2, Vladimir Reinharz1 and Jérôme Waldispühl1* From 14th Annual Research in Computational Molecular Biology (RECOMB) Comparative Genomics Satellite Workshop Montreal, Canada. 11-14 October 2016 Abstract Background: Secondary structures form the scaffold of multiple sequence alignment of non-coding RNA (ncRNA) families. An accurate reconstruction of ancestral ncRNAs must use this structural signal. However, the inference of ancestors of a single ncRNA family with a single consensus structure may bias the results towards sequences with high affinity to this structure, which are far from the true ancestors. Methods: In this paper, we introduce achARNement, a maximum parsimony approach that, given two alignments of homologous ncRNA families with consensus secondary structures and a phylogenetic tree, simultaneously calculates ancestral RNA sequences for these two families. Results: We test our methodology on simulated data sets, and show that achARNement outperforms classical maximum parsimony approaches in terms of accuracy, but also reduces by several orders of magnitude the number of candidate sequences. To conclude this study, we apply our algorithms on the Glm clan and the FinP-traJ clan from the Rfam database. Conclusions: Our results show that our methods reconstruct small sets of high-quality candidate ancestors with better agreement to the two target structures than with classical approaches. Our program is freely available at: http://csb.cs.mcgill.ca/acharnement. Keywords: RNA, Secondary structure, Ancestor reconstruction, Evolution, Phylogeny, Algorithm Background For a long time, most of the attention has been With the development of sequencing technologies given to the reconstruction of ancient protein and DNA emerged the need to elucidate the relationship between sequences, while RNA molecules remained relatively sequences from various organisms.
    [Show full text]
  • Were Ancestral Proteins Less Specific?
    bioRxiv preprint doi: https://doi.org/10.1101/2020.05.27.120261; this version posted May 30, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-ND 4.0 International license. Were ancestral proteins less specific? 1 1,2,3 1,2* Lucas C. Wheeler and Michael J. Harms 2 1. Institute of Molecular Biology, University of Oregon, Eugene OR 97403 3 2. Department of Chemistry and Biochemistry, University of Oregon, Eugene OR 97403 4 3. Department of Ecology and Evolutionary Biology, University of Colorado, Boulder 5 CO 80309 6 7 1 bioRxiv preprint doi: https://doi.org/10.1101/2020.05.27.120261; this version posted May 30, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-ND 4.0 International license. Abstract 8 Some have hypothesized that ancestral proteins were, on average, less specific than their 9 descendants. If true, this would provide a universal axis along which to organize protein 10 evolution and suggests that reconstructed ancestral proteins may be uniquely powerful tools 11 for protein engineering. Ancestral sequence reconstruction studies are one line of evidence 12 used to support this hypothesis. Previously, we performed such a study, investigating the 13 evolution of peptide binding specificity for the paralogs S100A5 and S100A6.
    [Show full text]
  • Ancestral Reconstruction
    TOPIC PAGE Ancestral Reconstruction Jeffrey B. Joy1,2, Richard H. Liang1, Rosemary M. McCloskey1, T. Nguyen1, Art F. Y. Poon1,2* 1 BC Centre for Excellence in HIV/AIDS, Vancouver, British Columbia, Canada, 2 University of British Columbia, Department of Medicine, Vancouver, British Columbia, Canada * [email protected] Introduction Ancestral reconstruction is the extrapolation back in time from measured characteristics of a11111 individuals (or populations) to their common ancestors. It is an important application of phylogenetics, the reconstruction and study of the evolutionary relationships among individu- als, populations, or species to their ancestors. In the context of biology, ancestral reconstruction can be used to recover different kinds of ancestral character states, including the genetic sequence (ancestral sequence reconstruction), the amino acid sequence of a protein, the com- position of a genome (e.g., gene order), a measurable characteristic of an organism (pheno- OPEN ACCESS type), and the geographic range of an ancestral population or species (ancestral range reconstruction). Nonbiological applications include the reconstruction of the vocabulary or Citation: Joy JB, Liang RH, McCloskey RM, Nguyen phonemes of ancient languages [1] and cultural characteristics of ancient societies such as oral T, Poon AFY (2016) Ancestral Reconstruction. PLoS Comput Biol 12(7): e1004763. doi:10.1371/journal. traditions [2] or marriage practices [3]. pcbi.1004763 Ancestral reconstruction relies on a sufficiently realistic model of evolution to accurately recover ancestral states. No matter how well the model approximates the actual evolutionary Published: July 12, 2016 history, however, one's ability to accurately reconstruct an ancestor deteriorates with increasing Copyright: © 2016 Joy et al.
    [Show full text]
  • Workshop on Molecular Evolution (Extended Special Topics Session July 27-August 8,2003 August 8-August 15,2003) Course Director
    Workshop on Molecular Evolution July 27-August 8,2003 (Extended Special Topics Session August 8-August 15,2003) Course Director: Michael P. Cummings, University of Maryland and Marine Biological Laboratory Molecular evolution has become the nexus of many areas of biological research. It both brings together and enriches such areas as biochemistry, molecular biology, microbiology, population genetics, systematics, developmental biology, genomics, bioinformatics, in vitro evolution, and molecular ecology. The Workshop provides an important contribution to these fields in that it promotes interdisciplinary research and interaction, and thus provides a glue that sticks together disparate fields. Due to the wide range of fields addressed by the study of molecular evolution, it is difficult to offer a comprehensive course in a university setting. It is rare for a single institution to maintain expertise in all necessary areas. In contrast, the Workshop is uniquely able to provide necessary breadth and depth by utilizing a large number of faculty with appropriate expertise. Furthermore, the flexible nature of the Workshop allows for rapid adaptation to changes in the dynamic field of molecular evolution. For example, the 2003 Workshop included recently emergent research areas of molecular evolution of development and genomics. The interest in the Workshop remains very strong and is increasing. The number of applications for the 2003 course was 143, continuing the trend of increased applications since 2000. In 2003 there were 60 students participating in the Workshop, which was taught by 19 faculty and 4 teaching assistants. The students came from all over the world (1 7 countries), and represented several career stages: graduate students (57%), postdoctoral researchers (1 3%), faculty/principal investigators (27%), and other (3%).
    [Show full text]
  • Reconstructing Ancestral Gene Content by Coevolution
    Downloaded from genome.cshlp.org on October 3, 2021 - Published by Cold Spring Harbor Laboratory Press Methods Reconstructing ancestral gene content by coevolution Tamir Tuller,1,2,3,4,5 Hadas Birin,1,4 Uri Gophna,2 Martin Kupiec,2 and Eytan Ruppin1,3 1School of Computer Sciences, Tel Aviv University, Ramat Aviv 69978, Israel; 2Department of Molecular Microbiology and Biotechnology, Tel Aviv University, Ramat Aviv 69978, Israel; 3School of Medicine, Tel Aviv University, Ramat Aviv 69978, Israel Inferring the gene content of ancestral genomes is a fundamental challenge in molecular evolution. Due to the statistical nature of this problem, ancestral genomes inferred by the maximum likelihood (ML) or the maximum-parsimony (MP) methods are prone to considerable error rates. In general, these errors are difficult to abolish by using longer genomic sequences or by analyzing more taxa. This study describes a new approach for improving ancestral genome reconstruction, the ancestral coevolver (ACE), which utilizes coevolutionary information to improve the accuracy of such reconstructions over previous approaches. The principal idea is to reduce the potentially large solution space by choosing a single optimal (or near optimal) solution that is in accord with the coevolutionary relationships between protein families. Simulation experiments, both on artificial and real biological data, show that ACE yields a marked decrease in error rate compared with ML or MP. Applied to a large data set (95 organisms, 4873 protein families, and 10,000 coevolutionary relationships), some of the ancestral genomes reconstructed by ACE were remarkably different in their gene content from those reconstructed by ML or MP alone (more than 10% in some nodes).
    [Show full text]
  • Molecular Phylogenetics Next Two Lectures
    Next two lectures: • What is a phylogenetic tree? • How are trees inferred using molecular data? Molecular phylogenetics • How do you assess confidence in trees and clades on trees? Michael Worobey • Why do trees for different data sets sometimes conflict? • What can you do with trees beyond simply inferring relatedness? Molecular phylogenetics fundamentals Tree terminology All of life is related by common ancestry. Recovering this pattern, the "Tree of Life", is one of the primary goals of evolutionary biology. Even at the A tree is a mathematical structure which is used to model the actual evolutionary population level, the phylogenetic tree is indispensable as a tool for estimating history of a group of sequences or organisms. This actual pattern of historical parameters of interest. Likewise at the among species level, it is relationships is the phylogeny or evolutionary tree which we try and indispensable for examining patterns of diversification over time. First, you estimate. A tree consists of nodes connected by branches (also called need to be familiar with some tree terminology. edges). Terminal nodes (also called leaves, OTUs [Operational Taxonomic Units], external nodes or terminal taxa) represent sequences or organisms for which we have data; they may be either extant or extinct. Phylogenetics interlude Phylogenetics interlude • It’s all about ancestor and All the phylogenetics offspring populations, • Mutations happen when you need to know, in lineages branching genetic material is copied 5 minutes… • The ancestor could
    [Show full text]
  • Climatic Niches in Phylogenetic Comparative Studies: a Review Of
    bioRxiv preprint doi: https://doi.org/10.1101/018796; this version posted April 30, 2015. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. 1 Climatic niches in phylogenetic comparative studies: a review of 2 challenges and approaches 1 1 3 Lara Budic and Carsten F. Dormann 4 Running title: Climatic niche in phylogenetic comparative studies 1 5 Faculty of Environment and Natural Resources, University of Freiburg, Tennenbacher Str. 4, 6 79106 Freiburg 7 Corresponding author: 8 Lara Budic 9 e-mail: [email protected] 10 Key words: Climatic niche, phylogenetic comparative studies, ancestral niche reconstruction, 11 evolutionary models 1 bioRxiv preprint doi: https://doi.org/10.1101/018796; this version posted April 30, 2015. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. 12 Abstract Studying the evolution of climatic niches through time in a phylogenetic 13 comparative framework combines species distribution modeling with phylogenies. 14 Phylogenetic comparative studies aid the understanding of the evolution of species' 15 environmental preferences by revealing the underlying evolutionary processes and causes, 16 detecting the differences among groups of species or relative to evolutionary pattern of other 17 phenotypic traits, but also act as a yardstick to gauge the adaptational potential under 18 climate change.
    [Show full text]
  • American Society of Naturalists Honorary Lifetime Membership Awards
    The University of Chicago American Society of Naturalists Honorary Lifetime Membership Awards. Source: The American Naturalist, Vol. 183, No. 6 (June 2014), pp. ii-v Published by: The University of Chicago Press for The American Society of Naturalists Stable URL: http://www.jstor.org/stable/10.1086/676468 . Accessed: 28/09/2014 17:59 Your use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at . http://www.jstor.org/page/info/about/policies/terms.jsp . JSTOR is a not-for-profit service that helps scholars, researchers, and students discover, use, and build upon a wide range of content in a trusted digital archive. We use information technology and tools to increase productivity and facilitate new forms of scholarship. For more information about JSTOR, please contact [email protected]. The University of Chicago Press, The American Society of Naturalists, The University of Chicago are collaborating with JSTOR to digitize, preserve and extend access to The American Naturalist. http://www.jstor.org This content downloaded from 142.103.160.110 on Sun, 28 Sep 2014 17:59:55 PM All use subject to JSTOR Terms and Conditions American Society of Naturalists Honorary Lifetime Membership Awards Jane Lubchenco The American Society of Naturalists is pleased to award step of effective involvement in the political processes Jane Lubchenco an Honorary Lifetime Membership. Jane that direct and fund our endeavors. She has cofounded received her BA in biology from Colorado College, an three organizations: the Leopold Leadership Program, the MS in zoology from the University of Washington, and Communications Partnership for Science and the Sea, her doctorate from Harvard University, studying with and Climate Central.
    [Show full text]
  • Predicting the Ancestral Character Changes in a Tree Is Typically Easier Than Predicting the Root State
    Version dated: September 5, 2013 Predicting ancestral states in a tree Predicting the ancestral character changes in a tree is typically easier than predicting the root state Olivier Gascuel1 and Mike Steel2 1Institut de Biologie Computationnelle, LIRMM, CNRS & Universit´ede Montpellier, France; 2Allan Wilson Centre for Molecular Ecology and Evolution, University of Canterbury, Christchurch, New Zealand Corresponding author: Mike Steel, University of Canterbury, Christchurch, New Zealand; E-mail: [email protected] Abstract.|Predicting the ancestral sequences of a group of homologous sequences related by a phylogenetic tree has been the subject of many studies, and numerous methods have been proposed to this purpose. Theoretical results are available that show that when the arXiv:1309.0926v1 [q-bio.PE] 4 Sep 2013 mutation rate become too large, reconstructing the ancestral state at the tree root is no longer feasible. Here, we also study the reconstruction of the ancestral changes that occurred along the tree edges. We show that, depending on the tree and branch length distribution, reconstructing these changes (i.e. reconstructing the ancestral state of all internal nodes in the tree) may be easier or harder than reconstructing the ancestral root state. However, results from information theory indicate that for the standard Yule tree, 1 the task of reconstructing internal node states remains feasible, even for very high substitution rates. Moreover, computer simulations demonstrate that for more complex trees and scenarios, this result still holds. For a large variety of counting, parsimony-based and likelihood-based methods, the predictive accuracy of a randomly selected internal node in the tree is indeed much higher than the accuracy of the same method when applied to the tree root.
    [Show full text]
  • Reconstructing Large Regions of an Ancestral Mammalian Genome in Silico
    Letter Reconstructing large regions of an ancestral mammalian genome in silico Mathieu Blanchette,1,4,5 Eric D. Green,2 Webb Miller,3 and David Haussler1,5 1Howard Hughes Medical Institute, University of California, Santa Cruz, California 95064, USA; 2National Human Genome Research Institute, National Institutes of Health, Bethesda, Maryland 20892, USA; 3Department of Biology, Pennsylvania State University, University Park, Pennsylvania 16802, USA It is believed that most modern mammalian lineages arose from a series of rapid speciation events near the Cretaceous-Tertiary boundary. It is shown that such a phylogeny makes the common ancestral genome sequence an ideal target for reconstruction. Simulations suggest that with methods currently available, we can expect to get 98% of the bases correct in reconstructing megabase-scale euchromatic regions of an eutherian ancestral genome from the genomes of ∼20 optimally chosen modern mammals. Using actual genomic sequences from 19 extant mammals, we reconstruct 1.1 Mb of ancient genome sequence around the CFTR locus. Detailed examination suggests the reconstruction is accurate and that it allows us to identify features in modern species, such as remnants of ancient transposon insertions, that were not identified by direct analysis. Tracing the predicted evolutionary history of the bases in the reconstructed region, estimates are made of the amount of DNA turnover due to insertion, deletion, and substitution in the different placental mammalian lineages since the common eutherian ancestor, showing considerable variation between lineages. In coming years, such reconstructions may help in identifying and understanding the genetic features common to eutherian mammals and may shed light on the evolution of human or primate-specific traits.
    [Show full text]