Inferring the Probability of the Derived Vs. the Ancestral Allelic State at a Polymorphic Site

Total Page:16

File Type:pdf, Size:1020Kb

Inferring the Probability of the Derived Vs. the Ancestral Allelic State at a Polymorphic Site | INVESTIGATION Inferring the Probability of the Derived vs. the Ancestral Allelic State at a Polymorphic Site Peter D. Keightley1 and Benjamin C. Jackson Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, EH9 3FL, United Kingdom ABSTRACT It is known that the allele ancestral to the variation at a polymorphic site cannot be assigned with certainty, and that the most frequently used method to assign the ancestral state—maximum parsimony—is prone to misinference. Estimates of counts of sites that have a certain number of copies of the derived allele in a sample (the unfolded site frequency spectrum, uSFS) made by parsimony are therefore also biased. We previously developed a maximum likelihood method to estimate the uSFS for a focal species using information from two outgroups while assuming simple models of nucleotide substitution. Here, we extend this approach to allow multiple outgroups (implemented for three outgroups), potentially any phylogenetic tree topology, and more complex models of nucleotide substitution. We find, however, that two outgroups and the Kimura two-parameter model are adequate for uSFS inference in most cases. We show that using parsimony to infer the ancestral state at a specific site seriously breaks down in two situations. The first is where the outgroups provide no information about the ancestral state of variation in the focal species. In this case, nucleotide variation will be underestimated if such sites are excluded. The second is where the minor allele in the focal species agrees with the allelic state of the outgroups. In this situation, parsimony tends to overestimate the probability of the major allele being derived, because it fails to account for the fact that sites with a high frequency of the derived allele tend to be rare. We present a method that corrects this deficiency and is capable of providing nearly unbiased estimates of ancestral state probabilities on a site-by-site basis and the uSFS. KEYWORDS nucleotide polymorphism; ancestral allele; derived allele; unfolded site frequency spectrum; parsimony; misinference ANY population genetic and quantitative genetic anal- derived alleles across sites, the unfolded site frequency spec- Mysis methods require the assignment of ancestral vs. trum (uSFS) (e.g., Boyko et al. 2008; Schneider et al. 2011; derived states at polymorphic nucleotide sites. For example, Tataru et al. 2017). Fay and Wu (2000) and Zeng et al. (2006) proposed statis- The minor allele at a site or counts of numbers of minor tics, H and E, that compare the numbers of high, intermedi- alleles at a group of sites (the folded site frequency spectrum) ate, and low frequency derived variants, which can then be can be observed directly fromsequence polymorphism data.In used to distinguish between different modes of natural selec- contrast, the derived vs. the ancestral allele at a site cannot be tion and demographic change. A number of methods have known with certainty, because at least one outgroup is re- also been developed to infer selection and demographic quired for inference, and there is the possibility of more than change based on the complete distribution of counts of one mutation separating the focal species from the outgroup. This also implies that the uSFS cannot be known precisely. Copyright © 2018 Keightley and Jackson For the purpose of ancestral state inference, rule-based max- doi: https://doi.org/10.1534/genetics.118.301120 Manuscript received March 23, 2018; accepted for publication May 14, 2018; imum parsimony is the most frequently applied method in published Early Online May 16, 2018. molecular evolutionary genetics (e.g., Voight et al. 2006; Available freely online through the author-supported open access option. Dreszer et al. 2007; Keinan et al. 2007; Sabeti et al. 2007; This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which 1000 Genomes Project Consortium 2010, 2015; Lohse and permits unrestricted use, distribution, and reproduction in any medium, provided the Barton 2011; Langley et al. 2012; Schmidt et al. 2017). It has original work is properly cited. Supplemental material available at Figshare: https://doi.org/10.25386/genetics. been recognized, however, that parsimony potentially pro- 6275915. duces misleading results (Felsenstein 1981; Collins et al. 1Corresponding author: Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Charlotte Auerbach Rd., Edinburgh EH9 3FL, United 1994; Eyre-Walker 1998). Of particular relevance here is that Kingdom. E-mail: peter.keightley@ed.ac.uk sites that have a low frequency of the derived allele are Genetics, Vol. 209, 897–906 July 2018 897 usually more common than sites that have a high frequency infer ancestral state probabilities for each polymorphic site of the derived allele. This implies that misinference tends in the data. We evaluate the performance of the new ap- to upwardly bias counts of high frequency derived alleles proach by simulations, apply it to data from the Drosophila (Baudry and Depaulis 2003; Hernandez et al. 2007). Population Genomics Project (DPGP) as a test case, and rein- There is a related problem concerning the assignment of fer the ancestral state probabilities for a population of the ancestral states, which does not seem to have been addressed. 1000 Genomes Project in humans, which were previously If ancestral states are assigned on site-by-site basis, potentially inferred by a parsimony-related approach. useful information is ignored. For example, consider the case of a single outgroup species that is uninformative about the ancestral allele of the variation in a focal species at a site. It is Materials and Methods more likely, however, that the ancestral allele at the site is the Following Keightley et al. (2016), uSFS inference is carried low frequency allele, if sites with a high frequency of the out in two-steps. Evolutionary rate parameters are estimated derived allele are uncommon in the data set as a whole (as is from all sites in the data (including polymorphic and mono- usually the case). morphic sites) in step 1. In step 2 the uSFS is computed, Matsumoto et al. (2015) pointed out that ancestral states conditional on the evolutionary rate parameter estimates. are not observable, that a single best ancestral reconstruction Information from steps 1 and 2 is then combined in a third is not advisable, and that assuming one can bias molecular step to infer the ancestral state probability for each polymor- evolutionary inference. This was developed by Jackson et al. phic site. (2017), who assigned the ancestral state probability at a site fi as the inferred probability of the node for the common an- Representation of the data and some de nitions cestor of the focal species and the closest outgroup, obtained Suppose we have sampled m orthologous gene copies at a set using PAML (Yang 2007), while ignoring polymorphism of sites from a population of a focal species. The uSFS we data. However, this does not optimally weight information require to estimate therefore contains m 2 1 elements, ex- coming from the focal site itself and from the data as a whole. cluding the elements where the ancestral or derived allele is Inference of ancestral states on a site-by-site basis has been fixed. We assume that we have randomly sampled a single problematic, but there has been progress in inferring the uSFS. gene copy at each site in one or more outgroup species. We Hernandez et al. (2007) developed a context-dependent sub- assume that the tree topology relating the species is known stitution model using a single outgroup to infer the ancestral and does not vary among sites (Figure 1). In the analysis we state at a polymorphic site in a focal species, and then imple- assume that the nucleotide variation within the focal species fi mented a step to correct for ancestral misidenti cation. Their coalesces within the branch labeled b1. The consequences of simulations suggested, however, that the approach only par- polymorphism in the outgroup species and violation of the tially corrects for ancestral misidentification, depending on assumptions of an invariant tree topology and coalescence the divergence between the focal species and the outgroup. within branch b1 are investigated in simulations. The ob- Schneider et al. (2011) developed a probabilistic method served nucleotide configuration for a site is the count of each to infer the uSFS on a site-by-site basis, but did not use infor- of the four nucleotides in the focal species (labeled X, Y for a mation from the frequencies of polymorphisms across all sites, biallelic site), along with the state for each outgroup (A, C, G, so results from this method are biased. Keightley et al. (2016) or T). Let the number of outgroups = n (in Figure 1, n = 3), developed a maximum likelihood (ML) method that addresses and denote the outgroups o1, o2...on. Assuming an unrooted the deficiency in Schneider et al. (2011) and simulations sug- tree (as in Figure 1), the number of branches in the tree is gested that it is capable of correctly inferring the uSFS. It uses a therefore b =2n 2 1. two-stage process in which the evolutionary rates are estimated Models of nucleotide substitution by ML and then, assuming the rates, estimates the uSFS ele- ments by ML, while correctly weighting information from in- The JC model, K2 model, and a model allowing six symmet- formative and uninformative sites. However, the method is rical rates (R6; Figure 2) are considered. All substitution limited to two outgroups, assumes simple substitution models models require the estimation of evolutionary rates (i.e., [for one outgroup, the Kimura two-parameter (K2) model; for mean number of nucleotide changes per site) for each – two outgroups, the Jukes Cantor (JC) model], and is not read- branch, K1...Kb.
Recommended publications
  • Genetics Society News
    July 2008 . ISSUE 59 GENETICS SOCIETY NEWS www.genetics.org.uk IN THIS ISSUE Genetics Society News is edited by Steve Russell. Items for future issues should be sent to Steve Russell, preferably by email to • Genetics Society Epigenetics Meeting s.russell@gen.cam.ac.uk, or hard copy to Department of Genetics, • Genetics Society Sponsored Meetings University of Cambridge, Downing Street, Cambridge CB2 3EH. The Newsletter is published twice a year, with copy dates of 1st June and • Travel, Fieldwork and Studentship Reports 26th November. • John Evans: an Appreciation Cocoons of the parasitoid wasp Cotesia vestalis on cabbage leaf in Taiwan. From the • Twelve Galton Lectures fieldwork report by Jetske G. de Boer on page 36. • My Favourite Paper A WORD FROM THE EDITOR A word from the editor ow soon until the $1000 based on the results of tests we genome is actually with barely understand! Here in the Hus and individual UK there is currently a sequencing is widespread? The moratorium, adhered to by publication of increasing most insurers, on the use of numbers of individual human genetic testing information for genome sequences suggests assessing life insurance that we should start to consider applications. It is important some of the implications that this remains in place and associated with the availability its effectiveness is reviewed of personal genetic well before the current information. In this issue we moratorium expires in 2011. present two articles reflecting The Human Genetics on his issue: a report from a Commission Genetics Society sponsored (http://www.hgc.gov.uk) meeting recently held in monitor issues relating to Cambridge organised by The genetic discrimination in the Triple Helix, an international UK and are a point of contact undergraduate organisation, as for those with any concerns in the Millennium Technology Prize.
    [Show full text]
  • Estimation of the Spontaneous Mutation Rate in Heliconius Melpomene
    Estimation of the Spontaneous Mutation Rate in Heliconius melpomene The Harvard community has made this article openly available. Please share how this access benefits you. Your story matters Citation Keightley, Peter D., Ana Pinharanda, Rob W. Ness, Fraser Simpson, Kanchon K. Dasmahapatra, James Mallet, John W. Davey, and Chris D. Jiggins. 2014. “Estimation of the Spontaneous Mutation Rate in Heliconius melpomene.” Molecular Biology and Evolution 32 (1): 239-243. doi:10.1093/molbev/msu302. http://dx.doi.org/10.1093/ molbev/msu302. Published Version doi:10.1093/molbev/msu302 Citable link http://nrs.harvard.edu/urn-3:HUL.InstRepos:13890662 Terms of Use This article was downloaded from Harvard University’s DASH repository, and is made available under the terms and conditions applicable to Other Posted Material, as set forth at http:// nrs.harvard.edu/urn-3:HUL.InstRepos:dash.current.terms-of- use#LAA Estimation of the Spontaneous Mutation Rate in Heliconius melpomene Peter D. Keightley,*,1 Ana Pinharanda,2 Rob W. Ness,1 Fraser Simpson,3 Kanchon K. Dasmahapatra,3,4 James Mallet,3,5 John W. Davey,2 and Chris D. Jiggins2 1Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom 2Department of Zoology, University of Cambridge, Cambridge, United Kingdom 3Department of Genetics, Evolution and Environment, University College London, London, United Kingdom 4Department of Biology, University of York, York, United Kingdom 5Department of Organismic and Evolutionary Biology, Harvard University *Corresponding author: E-mail: peter.keightley@ed.ac.uk. Associate editor:JohnParsch Abstract We estimated the spontaneous mutation rate in Heliconius melpomene by genome sequencing of a pair of parents and 30oftheiroffspring,basedontheratioofnumberofdenovoheterozygotes to the number of callable site-individuals.
    [Show full text]
  • Letter Estimation of the Spontaneous Mutation Rate In
    Estimation of the Spontaneous Mutation Rate in Heliconius melpomene Peter D. Keightley,*,1 Ana Pinharanda,2 Rob W. Ness,1 Fraser Simpson,3 Kanchon K. Dasmahapatra,3,4 James Mallet,3,5 John W. Davey,2 and Chris D. Jiggins2 1Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom 2Department of Zoology, University of Cambridge, Cambridge, United Kingdom 3Department of Genetics, Evolution and Environment, University College London, London, United Kingdom 4Department of Biology, University of York, York, United Kingdom 5Department of Organismic and Evolutionary Biology, Harvard University *Corresponding author: E-mail: peter.keightley@ed.ac.uk. Associate editor:JohnParsch Abstract We estimated the spontaneous mutation rate in Heliconius melpomene by genome sequencing of a pair of parents and 30oftheiroffspring,basedontheratioofnumberofdenovoheterozygotes to the number of callable site-individuals. We detected nine new mutations, each one affecting a single site in a single offspring. This yields an estimated mutation Downloaded from rate of 2.9 Â 10À9 (95% confidence interval, 1.3 Â 10À9–5.5 Â 10À9), which is similar to recent estimates in Drosophila melanogaster, the only other insect species in which the mutation rate has been directly estimated. We infer that recent effective population size of H. melpomene is about 2 million, a substantially lower value than its census size, suggesting a role for natural selection reducing diversity. We estimate that H. melpomene diverged from its Mullerian€ comimic H. erato about 6 Ma, a somewhat later date than estimatesbasedonalocalmolecularclock. http://mbe.oxfordjournals.org/ Key words: mutation, Heliconius, genome sequencing. Understanding the process of spontaneous mutation is cen- as synonymous sites (Drake et al.
    [Show full text]
  • Directory 2016/17 the Royal Society of Edinburgh
    cover_cover2013 19/04/2016 16:52 Page 1 The Royal Society of Edinburgh T h e R o Directory 2016/17 y a l S o c i e t y o f E d i n b u r g h D i r e c t o r y 2 0 1 6 / 1 7 Printed in Great Britain by Henry Ling Limited, Dorchester, DT1 1HD ISSN 1476-4334 THE ROYAL SOCIETY OF EDINBURGH DIRECTORY 2016/2017 PUBLISHED BY THE RSE SCOTLAND FOUNDATION ISSN 1476-4334 The Royal Society of Edinburgh 22-26 George Street Edinburgh EH2 2PQ Telephone : 0131 240 5000 Fax : 0131 240 5024 email: rse@royalsoced.org.uk web: www.royalsoced.org.uk Scottish Charity No. SC 000470 Printed in Great Britain by Henry Ling Limited CONTENTS THE ORIGINS AND DEVELOPMENT OF THE ROYAL SOCIETY OF EDINBURGH .....................................................3 COUNCIL OF THE SOCIETY ..............................................................5 EXECUTIVE COMMITTEE ..................................................................6 THE RSE SCOTLAND FOUNDATION ..................................................7 THE RSE SCOTLAND SCIO ................................................................8 RSE STAFF ........................................................................................9 LAWS OF THE SOCIETY (revised October 2014) ..............................13 STANDING COMMITTEES OF COUNCIL ..........................................27 SECTIONAL COMMITTEES AND THE ELECTORAL PROCESS ............37 DEATHS REPORTED 26 March 2014 - 06 April 2016 .....................................................43 FELLOWS ELECTED March 2015 ...................................................................................45
    [Show full text]
  • Letter Estimation of the Spontaneous Mutation Rate in Heliconius
    Estimation of the Spontaneous Mutation Rate in Heliconius melpomene Peter D. Keightley,*,1 Ana Pinharanda,2 Rob W. Ness,1 Fraser Simpson,3 Kanchon K. Dasmahapatra,3,4 James Mallet,3,5 John W. Davey,2 and Chris D. Jiggins2 1Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom 2Department of Zoology, University of Cambridge, Cambridge, United Kingdom 3Department of Genetics, Evolution and Environment, University College London, London, United Kingdom 4Department of Biology, University of York, York, United Kingdom 5Department of Organismic and Evolutionary Biology, Harvard University *Corresponding author: E-mail: peter.keightley@ed.ac.uk. Associate editor:JohnParsch Abstract We estimated the spontaneous mutation rate in Heliconius melpomene by genome sequencing of a pair of parents and 30oftheiroffspring,basedontheratioofnumberofdenovoheterozygotes to the number of callable site-individuals. Downloaded from We detected nine new mutations, each one affecting a single site in a single offspring. This yields an estimated mutation rate of 2.9 Â 10À9 (95% confidence interval, 1.3 Â 10À9–5.5 Â 10À9), which is similar to recent estimates in Drosophila melanogaster, the only other insect species in which the mutation rate has been directly estimated. We infer that recent effective population size of H. melpomene is about 2 million, a substantially lower value than its census size, suggesting a role for natural selection reducing diversity. We estimate that H. melpomene diverged from its Mullerian€ comimic H. erato http://mbe.oxfordjournals.org/ about 6 Ma, a somewhat later date than estimatesbasedonalocalmolecularclock. Key words: mutation, Heliconius, genome sequencing. Understanding the process of spontaneous mutation is cen- as synonymous sites (Drake et al.
    [Show full text]
  • T:\Year Book\2003\Review\Pagema
    The Royal Society of Edinburgh Review 2003 (Session 2001-2002) The Royal Society of Edinburgh Review 2003 THE ROYAL SOCIETY OF EDINBURGH REVIEW OF THE SESSION 2001-2002 The Royal Society of Edinburgh 22-26 George Street Edinburgh, EH2 2PQ Telephone : 0131 240 5000 Fax : 0131 240 5024 email : rse@royalsoced.org.uk Scottish Charity No SC000470 Printed in Great Britain by J W Arrowsmith Ltd, Bristol, BS3 2NT Cover illustration by Aird McKinstrie. Design by Jennifer Cameron THE ROYAL SOCIETY OF EDINBURGH REVIEW OF THE SESSION 2001-2002 PUBLISHED BY THE RSE SCOTLAND FOUNDATION ISSN 1476-4342 CONTENTS Proceedings of the Ordinary Meetings .................................................. 3 Proceedings of the Annual Statutory Meeting ...................................... 5 Trustees report to 31 March 2002 ...................................................... 17 Auditor’s Report and Accounts ........................................................... 29 Schedule of Investments .................................................................... 51 Activities Prize Lectures ............................................................................... 55 Lectures ....................................................................................... 71 Conferences, Symposia, Workshops and Exhibitions .................... 127 Publications ................................................................................ 149 Fundraising ................................................................................. 151 The Scottish Science Advisory
    [Show full text]
  • Inferring Parameters of the Distribution of Fitness Effects of New Mutations When Beneficial Mutations Are Strongly Advantageous
    G3: Genes|Genomes|Genetics Early Online, published on May 5, 2020 as doi:10.1534/g3.120.401052 1 Inferring parameters of the distribution of 2 fitness effects of new mutations when 3 beneficial mutations are strongly 4 advantageous and rare 5 Tom R. Booker1,2 6 7 1. Department of Forest and Conservation Sciences, University of British Columbia, Vancouver, 8 Canada. 9 2. Biodiversity Research Centre, University of British Columbia, Vancouver, Canada 10 11 Correspondence: booker@zoology.ubc.ca 12 13 14 15 16 17 18 19 1 © The Author(s) 2020. Published by the Genetics Society of America. 20 Abstract 21 Characterising the distribution of fitness effects (DFE) for new mutations is central in evolutionary 22 genetics. Analysis of molecular data under the McDonald-Kreitman test has suggested that 23 adaptive substitutions make a substantial contribution to between-species divergence. Methods 24 have been proposed to estimate the parameters of the distribution of fitness effects for positively 25 selected mutations from the unfolded site frequency spectrum (uSFS). Such methods perform well 26 when beneficial mutations are mildly selected and frequent. However, when beneficial mutations 27 are strongly selected and rare, they may make little contribution to standing variation and will 28 thus be difficult to detect from the uSFS. In this study, I analyse uSFS data from simulated 29 populations subject to advantageous mutations with effects on fitness ranging from mildly to 30 strongly beneficial. As expected, frequent, mildly beneficial mutations contribute substantially to 31 standing genetic variation and parameters accurately recovered from the uSFS. However, when 32 advantageous mutations are strongly selected and rare, there are very few segregating in 33 populations at any one time.
    [Show full text]
  • 1 the Nature of Genetic Variation for Complex Traits Revealed by GWAS and Regional Heritability Mapping Analyses Armando Caballe
    Genetics: Early Online, published on October 19, 2015 as 10.1534/genetics.115.177220 The nature of genetic variation for complex traits revealed by GWAS and Regional Heritability Mapping analyses Armando Caballero*, Albert Tenesa§,† and Peter D. Keightley‡ * Departamento de Bioquímica, Genética e Inmunología, Facultad de Biología, Universidad de Vigo, Vigo, Spain § The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, United Kingdom † MRC HGU at the MRC IGMM, University of Edinburgh, Western General Hospital, Crewe Road South, Edinburgh, United Kingdom ‡ Institute of Evolutionary Biology, University of Edinburgh, Charlotte Auerbach Rd, Edinburgh, EH9 3FL, United Kingdom E-mail: armando@uvigo.es E-mail: Albert.Tenesa@ed.ac.uk E-mail: Peter.Keightley@ed.ac.uk 1 Copyright 2015. Running Head: Genetic basis of missing heritability Type of manuscript: Contributed paper Key words: quantitative trait variation; complex traits; missing heritability; fitness; additive genetic variance Word count: 9442 ADDRESS FOR CORRESPONDENCE: Armando Caballero Dpto. de Bioquímica, Genética e Inmunología, Facultad de Biología, Universidad de Vigo, 36310 Vigo, SPAIN E-mail: armando@uvigo.es Tel: +34 986 812 568. Fax: +34 986 812 556 2 Abstract We use computer simulations to investigate the amount of genetic variation for complex traits that can be revealed by single SNP Genome Wide Association Studies (GWAS) or Regional Heritability Mapping (RHM) analyses based on full genome sequence data or SNP chips. We model a large population subject to mutation, recombination, selection, and drift, assuming a pleiotropic model of mutations sampled from a bivariate distribution of effects of mutations on a quantitative trait and fitness.
    [Show full text]
  • Causes of Natural Variation in Fitness: Evidence from Studies Of
    Correction EVOLUTION Correction for “Causes of natural variation in fitness: Evidence from studies of Drosophila populations,” by Brian Charlesworth, which appeared in issue 6, February 10, 2015, of Proc Natl Acad Sci USA (112:1662–1669; first published January 8, 2015; 10.1073/ pnas.1423275112). The authors note that on page 1666, right column, first full paragraph, line 6, “an estimate of t of 0.044” should instead appear as “an estimate of t over both classes of mutation of 0.0073.” The authors also note that on page 1666, right column, fourth full paragraph, line 3, “although even this is likely to be a sub- stantial underestimate” should instead appear as “although the overall value is likely to be somewhat less than 0.01.” www.pnas.org/cgi/doi/10.1073/pnas.1502053112 CORRECTION www.pnas.org PNAS | March 3, 2015 | vol. 112 | no. 9 | E1049 Downloaded by guest on September 30, 2021 Causes of natural variation in fitness: Evidence from INAUGURAL ARTICLE studies of Drosophila populations Brian Charlesworth1 Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3JT, United Kingdom This contribution is part of the special series of Inaugural Articles by members of the National Academy of Sciences elected in 2013. Contributed by Brian Charlesworth, December 12, 2014 (sent for review September 18, 2014; reviewed by James Fry and David Houle) DNA sequencing has revealed high levels of variability within most quantitative genetics, respectively. The first approach sheds light species. Statistical methods based on population genetics theory on the general nature of the fitness effects of the DNA sequence have been applied to the resulting data and suggest that most variants found in natural populations, but says little about how mutations affecting functionally important sequences are deleteri- these fitness effects are caused.
    [Show full text]
  • Estimation of the Spontaneous Mutation Rate in Heliconius Melpomene' Molecular Biology and Evolution, Vol
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Edinburgh Research Explorer Edinburgh Research Explorer Estimation of the Spontaneous Mutation Rate in Heliconius melpomene Citation for published version: Keightley, PD, Pinharanda, A, Ness, RW, Simpson, F, Dasmahapatra, KK, Mallet, J, Davey, JW & Jiggins, CD 2014, 'Estimation of the Spontaneous Mutation Rate in Heliconius melpomene' Molecular Biology and Evolution, vol. 32, no. 1, pp. 239-243. DOI: 10.1093/molbev/msu302 Digital Object Identifier (DOI): 10.1093/molbev/msu302 Link: Link to publication record in Edinburgh Research Explorer Document Version: Publisher's PDF, also known as Version of record Published In: Molecular Biology and Evolution Publisher Rights Statement: The Author 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons. org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited General rights Copyright for the publications made accessible via the Edinburgh Research Explorer is retained by the author(s) and / or other copyright owners and it is a condition of accessing these publications that users recognise and abide by the legal requirements associated with these rights. Take down policy The University of Edinburgh has made every reasonable effort to ensure that Edinburgh Research Explorer content complies with UK legislation. If you believe that the public display of this file breaches copyright please contact openaccess@ed.ac.uk providing details, and we will remove access to the work immediately and investigate your claim.
    [Show full text]
  • Selective Interference Among Deleterious Mutations Favours Sex And
    Selective interference among deleterious mutations favours sex and recombination in finite populations regardless of the nature of epistasis Peter D. Keightley1 and Sarah P. Otto2 1Institute of Evolutionary Biology, University of Edinburgh, West Mains Road, Edinburgh EH9 3JT, UK 2Department of Zoology, University of British Columbia, 6270 University Blvd., Vancouver BC V6T 1Z4, Canada Corresponding author Prof. Peter Keightley Institute of Evolutionary Biology University of Edinburgh West Mains Road Edinburgh EH9 3JT UK Email: keightley.nature2006@spambob.net 1 Sex and recombination are widespread, but explaining these phenomena has been one of the most difficult problems in evolutionary biology. Recombination is advantageous when different individuals in a population carry different advantageous alleles1,2. By bringing together advantageous alleles onto the same chromosome, recombination speeds up the process of adaptation1,3,4,5 and opposes the fixation of harmful mutations via Muller's ratchet4,5. Nevertheless, adaptive substitutions favour sex and recombination only if the rate of adaptive mutation is high1,6, and Muller's ratchet operates only in small or asexual populations7. Here, by tracking the fate of modifier alleles that alter the frequency of sex and recombination, we demonstrate that background selection against deleterious mutant alleles provides a stochastic advantage to sex and recombination that increases with population size. The advantage arises because, with low levels of recombination, selection at other loci severely reduces the effective population size and genetic variance in fitness at a focal locus8 (the Hill- Robertson effect), making a population less able to respond to selection and to rid itself of deleterious mutations. Sex and recombination reveal the hidden genetic variance in fitness by combining chromosomes of intermediate fitness to create chromosomes that are relatively free of (or loaded with) deleterious mutations.
    [Show full text]
  • Estimation of the Spontaneous Mutation Rate Per Nucleotide Site in a Drosophila Melanogaster Full-Sib Family
    Genetics: Early Online, published on November 8, 2013 as 10.1534/genetics.113.158758 Estimation of the spontaneous mutation rate per nucleotide site in a Drosophila melanogaster full-sib family Peter D. Keightley, Rob W. Ness, Daniel L. Halligan and Penelope R. Haddrill Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, West Mains Road, Edinburgh EH9 3JT, UK Corresponding author: Peter Keightley, Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, West Mains Road, Edinburgh EH9 3JT, UK 1 Copyright 2013. Summary We employed deep genome sequencing of two parents and 12 of their offspring to estimate the mutation rate per site per generation in a full-sib family of D. melanogaster recently sampled from a natural population. Sites that were homozygous for the same allele in the parents and heterozygous in one or more offspring were categorized as candidate mutations and subjected to detailed analysis. In 1.23 x 109 callable sites from 12 individuals, we confirmed six single nucleotide mutations. We estimated the false negative rate in the experiment by generating synthetic mutations using the empirical distributions of numbers of non-reference bases at heterozygous sites in the offspring. The proportion of synthetic mutations at callable sites that we failed to detect was less than 1%, implying that the false negative rate was extremely low. Our estimate of the point mutation rate is 2.8 x 10-9 (95% confidence interval = 1.0 x 10-9 – 6.1 x 10-9) per site per generation, which is at the low end of the range of previous estimates, and suggests an effective population size for the species of ~1.4 x 106.
    [Show full text]