Challenges when applying stochastic models to reconstruct the demographic history of populations.

Size: px
Start display at page:

Download "Challenges when applying stochastic models to reconstruct the demographic history of populations."

Transcription

1 Challenges when applying stochastic models to reconstruct the demographic history of populations. Willy Rodríguez Institut de Mathématiques de Toulouse October 11, 2017

2 Outline 1 Introduction 2 Inverse Instantaneous Coalescence Rate (IICR) 3 Applications of the IICR 4 Conclusions

3 Introduction Inverse Instantaneous Coalescence Rate (IICR) Applications of the IICR Conclusions Presentation Outline 1 Introduction 2 Inverse Instantaneous Coalescence Rate (IICR) 3 Applications of the IICR 4 Conclusions Willy Rodríguez The IICR as a summary statistic

4 Goals of Population Genetics Build models that help us understand the main evolutionary events that gave rise to the observed patterns of genetic variation. Evolutionary forces

5 Reconstructing demographic history from DNA Reconstructing demographic history from DNA

6 Wright-Fisher model Wright-Fisher model ( ) Non-overlapping generations. Constant population size (2N genes, 1 indiv. = one gene). Random mating (panmixia).

7 Wright-Fisher model and key concepts Key concepts Coalescence events: Individuals I 3 and I 5 coalesce (have a common ancestor) 4 generations before the present. Probability that two individuals coalesce in the previous 1 generation: 2N. Probability that they do not coalesce: (1 1 2N ). G 2 : Number of generations to reach the ancestor of 2 genes. P(G 2 > l) = (1 1 2N )l. Let 2N be large, G 2 2N D T 2,with T 2 Exp(1).

8 The Coalescent and some approximations

9 The Coalescent and some approximations Kingman, J., The coalescent. A continuous-time Markov chain allowing to describe family relationships between individuals backward in time.

10 The Coalescent and some approximations Kingman, J., The coalescent. A continuous-time Markov chain allowing to describe family relationships between individuals backward in time. Hudson, R. R., The coalescent with recombination. A process allowing to trace back lineages, incorporating recombination events.

11 The Coalescent and some approximations Kingman, J., The coalescent. A continuous-time Markov chain allowing to describe family relationships between individuals backward in time. Hudson, R. R., The coalescent with recombination. A process allowing to trace back lineages, incorporating recombination events. Herbots, H. M. J. D The structured coalescent. An extension of Kingman s coalescent to subdivided populations.

12 The Coalescent and some approximations Kingman, J., The coalescent. A continuous-time Markov chain allowing to describe family relationships between individuals backward in time. Hudson, R. R., The coalescent with recombination. A process allowing to trace back lineages, incorporating recombination events. Herbots, H. M. J. D The structured coalescent. An extension of Kingman s coalescent to subdivided populations. McVean & Cardin, The Sequentially Markovian Coalescent. An approximation to the coalescent with recombination using a Markov chain along the genome.

13 Variable population size. Relation with T 2 Griffiths & Tavaré Distribution of T 2, function of population size change. Population evolving with deterministically varying size. λ(t) = N(t), with 2N: present population size (genes). 2N t 1 Λ(t) = 0 λ(u) du. Distribution of the coalescence times of two genes (T 2 ). P(T 2 > t) = 1 F T2 (t) = e Λ(t) f T2 (t) = ( F T2 (t) ) = 1 λ(t) e Λ(t). Objective: Reconstruct the function λ.

14 Methods for estimating past population size changes Methods for estimating past population size changes MSVAR (1999). Skyline plot (2000). Bayesian skyline plot (2005). dadi (2009) PSMC (2011). DiCal (2013). VarEff (2014). MSMC (2014). stairway plot (2015) PopSizeABC (2016). SMC++ (2017)

15 Methods for estimating past population size changes Methods for estimating past population size changes MSVAR (1999). Skyline plot (2000). Bayesian skyline plot (2005). dadi (2009) PSMC (2011). DiCal (2013). VarEff (2014). MSMC (2014). stairway plot (2015) PopSizeABC (2016). SMC++ (2017)

16 Methods for estimating past population size changes Methods for estimating past population size changes PSMC (2011).

17 The PSMC. Markov Chain along the genome Li and Durbin, 2011.Pairwise Sequentially Markovian Coalescent. The PSMC method Input data: one diploid genome (one diploid individual). Based on the Sequentially Markovian Coalescent (SMC. McVean & Cardin, 2005). Uses a relation between recombinations, mutations and λ k in a model of variable population size. Describes a Hidden Markov Model along the genome, allowing to estimate values of population size in the past. Developed to be applied on long DNA sequences (e.g.: one chromosome.)

18 PSMC. Markov Chain along the genome

19 PSMC. Markov Chain along the genome

20 PSMC. Markov Chain along the genome

21 PSMC. Markov Chain along the genome

22 PSMC. Markov Chain along the genome

23 PSMC. Hidden Markov Chain along the genome Discrete state space, discrete-time HMM. Hidden states: coalescence times (T 2 ) at position a (T 2 [t k, t k+1 ]). Observed states: homozygous or heterozygous at position a. Parameters: scaled mutation rate (θ), scaled recombination rate (ρ), demographic history (λ k ).

24 One example (the psmc) The psmc method (Heng Li & Richard Durbin, 2011) can be affected by population structure.

25 One example (the psmc). A simple structure model n islands model

26 One example (the psmc). A simple structure model What is going on here?

27 Introduction Inverse Instantaneous Coalescence Rate (IICR) Applications of the IICR Conclusions Presentation Outline 1 Introduction 2 Inverse Instantaneous Coalescence Rate (IICR) 3 Applications of the IICR 4 Conclusions Willy Rodríguez The IICR as a summary statistic

28 Coalescence times and population size changes. The panmixia hypothesis (random mating). What if we do not have random mating?

29 Instantaneous Coalescence Rate Griffiths & Tavaré (1994). Variable population size (N e ) N e (t) = N 0 λ(t) P(T 2 > t) = e t 0 1 λ(u) du Inverse Instantaneous Coalescence Rate (IICR). Mazet et al log(p(t 2 > t)) = 1 λ(t) λ(t) = P(T 2>t) f T2 (t) = 1 F T 2 (t) f T2 (t). λ can be evaluated at any t using only the T 2 distribution. Valid for any model. 1 T 2 can be interpreted as a lifetime. λ(t) : instantaneous coalescence rate (hazard function of failure rate). λ(t) may be disconnected with size changes.

30 Based on T 2, it is not possible to distinguish structure from a panmictic model with population size change.

31 Based on T 2, it is not possible to distinguish structure from a panmictic model with population size change. We can estimate the IICR using any method based on the panmixia hypothesis (ex: psmc).

32 Based on T 2, it is not possible to distinguish structure from a panmictic model with population size change. We can estimate the IICR using any method based on the panmixia hypothesis (ex: psmc). We can predict the IICR for any model by doing simulations Take a vector of independent values of T 2. Then: ˆλ(t) = 1 ˆF T2 (t). ˆf T2 (t)

33 Based on T 2, it is not possible to distinguish structure from a panmictic model with population size change. We can estimate the IICR using any method based on the panmixia hypothesis (ex: psmc). We can predict the IICR for any model by doing simulations Take a vector of independent values of T 2. Then: ˆλ(t) = 1 ˆF T2 (t). ˆf T2 (t)

34 Explicit expression for the IICR under the n islands model.

35 Introduction Inverse Instantaneous Coalescence Rate (IICR) Applications of the IICR Conclusions Presentation Outline 1 Introduction 2 Inverse Instantaneous Coalescence Rate (IICR) 3 Applications of the IICR 4 Conclusions Willy Rodríguez The IICR as a summary statistic

36 Trace the IICR for non panmictic models L. Chikhi et al. The IICR (inverse instantaneous coalescence rate) as a summary of genomic diversity: insights into demographic inference and model choice (to appear HDY)

37 Trace the IICR for non panmictic models L. Chikhi et al. The IICR (inverse instantaneous coalescence rate) as a summary of genomic diversity: insights into demographic inference and model choice (to appear HDY)

38 Trace the IICR for non panmictic models L. Chikhi et al. The IICR (inverse instantaneous coalescence rate) as a summary of genomic diversity: insights into demographic inference and model choice (to appear HDY)

39 Considering alternative scenarios

40 Considering alternative scenarios Courve fitting (

41 Considering alternative scenarios Courve fitting (

42 Considering alternative scenarios Courve fitting (

43 Considering alternative scenarios PSMC in human data. Alternative scenario: n islands model with changes in gene flow and recent grow

44 Considering alternative scenarios PSMC in human data. Alternative scenario: n islands model with changes in gene flow and recent grow Summary statistic in an ABC framework

45 Introduction Inverse Instantaneous Coalescence Rate (IICR) Applications of the IICR Conclusions Presentation Outline 1 Introduction 2 Inverse Instantaneous Coalescence Rate (IICR) 3 Applications of the IICR 4 Conclusions Willy Rodríguez The IICR as a summary statistic

46 Conclusions Identifiability problem when using T 2 (T k ) based methods. Reconstructed skyline plots are not meaningless but should be interpreted carefully (as an IICR) in a general case. For some models it is possible to have an explicit expression of the IICR. The IICR can be used as a summary statistic to explain the data under any model, provided we can simulate coalescence times. The IICR can be used as a summary statistic in an ABC framework.

47 Collaborators

48 Acknowledgement

49 Challenges when applying stochastic models to reconstruct the demographic history of populations. Willy Rodríguez Institut de Mathématiques de Toulouse October 11, 2017

50 Challenges when applying stochastic models to reconstruct the demographic history of populations. Willy Rodríguez Institut de Mathématiques de Toulouse October 11, 2017

Gene Genealogies Coalescence Theory. Annabelle Haudry Glasgow, July 2009

Gene Genealogies Coalescence Theory. Annabelle Haudry Glasgow, July 2009 Gene Genealogies Coalescence Theory Annabelle Haudry Glasgow, July 2009 What could tell a gene genealogy? How much diversity in the population? Has the demographic size of the population changed? How?

More information

Stochastic Demography, Coalescents, and Effective Population Size

Stochastic Demography, Coalescents, and Effective Population Size Demography Stochastic Demography, Coalescents, and Effective Population Size Steve Krone University of Idaho Department of Mathematics & IBEST Demographic effects bottlenecks, expansion, fluctuating population

More information

How robust are the predictions of the W-F Model?

How robust are the predictions of the W-F Model? How robust are the predictions of the W-F Model? As simplistic as the Wright-Fisher model may be, it accurately describes the behavior of many other models incorporating additional complexity. Many population

More information

Coalescence times for three genes provide sufficient information to distinguish population structure from population size changes

Coalescence times for three genes provide sufficient information to distinguish population structure from population size changes Coalescence times for three genes provide sufficient information to distinguish population structure from population size changes Simona Grusea, Willy Rodríguez, Didier Pinchon, Lounès Chikhi, Simon Boitard,

More information

Mathematical models in population genetics II

Mathematical models in population genetics II Mathematical models in population genetics II Anand Bhaskar Evolutionary Biology and Theory of Computing Bootcamp January 1, 014 Quick recap Large discrete-time randomly mating Wright-Fisher population

More information

Population Genetics I. Bio

Population Genetics I. Bio Population Genetics I. Bio5488-2018 Don Conrad dconrad@genetics.wustl.edu Why study population genetics? Functional Inference Demographic inference: History of mankind is written in our DNA. We can learn

More information

6 Introduction to Population Genetics

6 Introduction to Population Genetics Grundlagen der Bioinformatik, SoSe 14, D. Huson, May 18, 2014 67 6 Introduction to Population Genetics This chapter is based on: J. Hein, M.H. Schierup and C. Wuif, Gene genealogies, variation and evolution,

More information

Demography April 10, 2015

Demography April 10, 2015 Demography April 0, 205 Effective Population Size The Wright-Fisher model makes a number of strong assumptions which are clearly violated in many populations. For example, it is unlikely that any population

More information

6 Introduction to Population Genetics

6 Introduction to Population Genetics 70 Grundlagen der Bioinformatik, SoSe 11, D. Huson, May 19, 2011 6 Introduction to Population Genetics This chapter is based on: J. Hein, M.H. Schierup and C. Wuif, Gene genealogies, variation and evolution,

More information

Hidden Markov models in population genetics and evolutionary biology

Hidden Markov models in population genetics and evolutionary biology Hidden Markov models in population genetics and evolutionary biology Gerton Lunter Wellcome Trust Centre for Human Genetics Oxford, UK April 29, 2013 Topics for today Markov chains Hidden Markov models

More information

Coalescent based demographic inference. Daniel Wegmann University of Fribourg

Coalescent based demographic inference. Daniel Wegmann University of Fribourg Coalescent based demographic inference Daniel Wegmann University of Fribourg Introduction The current genetic diversity is the outcome of past evolutionary processes. Hence, we can use genetic diversity

More information

A sequentially Markov conditional sampling distribution for structured populations with migration and recombination

A sequentially Markov conditional sampling distribution for structured populations with migration and recombination A sequentially Markov conditional sampling distribution for structured populations with migration and recombination Matthias Steinrücken a, Joshua S. Paul b, Yun S. Song a,b, a Department of Statistics,

More information

arxiv: v1 [q-bio.pe] 9 Nov 2015

arxiv: v1 [q-bio.pe] 9 Nov 2015 Title: On the importance of being structured: instantaneous coalescence rates and a re-evaluation of human evolution Authors: Mazet Olivier 1,, Rodríguez Willy 1,, Grusea Simona 1, Boitard Simon 2,3, and

More information

Selection and Population Genetics

Selection and Population Genetics Selection and Population Genetics Evolution by natural selection can occur when three conditions are satisfied: Variation within populations - individuals have different traits (phenotypes). height and

More information

CS 68: BIOINFORMATICS. Prof. Sara Mathieson Swarthmore College Spring 2018

CS 68: BIOINFORMATICS. Prof. Sara Mathieson Swarthmore College Spring 2018 CS 68: BIOINFORMTICS Prof. Sara Mathieson Swarthmore College Spring 2018 Outline: pr 4 Lab 5 Examples HMM example in population genetics Recap Viterbi lgorithm Forward-Backward lgorithm Posterior Decoding

More information

Computational Systems Biology: Biology X

Computational Systems Biology: Biology X Bud Mishra Room 1002, 715 Broadway, Courant Institute, NYU, New York, USA Human Population Genomics Outline 1 2 Damn the Human Genomes. Small initial populations; genes too distant; pestered with transposons;

More information

A comparison of two popular statistical methods for estimating the time to most recent common ancestor (TMRCA) from a sample of DNA sequences

A comparison of two popular statistical methods for estimating the time to most recent common ancestor (TMRCA) from a sample of DNA sequences Indian Academy of Sciences A comparison of two popular statistical methods for estimating the time to most recent common ancestor (TMRCA) from a sample of DNA sequences ANALABHA BASU and PARTHA P. MAJUMDER*

More information

Endowed with an Extra Sense : Mathematics and Evolution

Endowed with an Extra Sense : Mathematics and Evolution Endowed with an Extra Sense : Mathematics and Evolution Todd Parsons Laboratoire de Probabilités et Modèles Aléatoires - Université Pierre et Marie Curie Center for Interdisciplinary Research in Biology

More information

The Wright-Fisher Model and Genetic Drift

The Wright-Fisher Model and Genetic Drift The Wright-Fisher Model and Genetic Drift January 22, 2015 1 1 Hardy-Weinberg Equilibrium Our goal is to understand the dynamics of allele and genotype frequencies in an infinite, randomlymating population

More information

Bustamante et al., Supplementary Nature Manuscript # 1 out of 9 Information #

Bustamante et al., Supplementary Nature Manuscript # 1 out of 9 Information # Bustamante et al., Supplementary Nature Manuscript # 1 out of 9 Details of PRF Methodology In the Poisson Random Field PRF) model, it is assumed that non-synonymous mutations at a given gene are either

More information

Introduction to population genetics & evolution

Introduction to population genetics & evolution Introduction to population genetics & evolution Course Organization Exam dates: Feb 19 March 1st Has everybody registered? Did you get the email with the exam schedule Summer seminar: Hot topics in Bioinformatics

More information

Diffusion Models in Population Genetics

Diffusion Models in Population Genetics Diffusion Models in Population Genetics Laura Kubatko kubatko.2@osu.edu MBI Workshop on Spatially-varying stochastic differential equations, with application to the biological sciences July 10, 2015 Laura

More information

Genetic Variation in Finite Populations

Genetic Variation in Finite Populations Genetic Variation in Finite Populations The amount of genetic variation found in a population is influenced by two opposing forces: mutation and genetic drift. 1 Mutation tends to increase variation. 2

More information

EVOLUTIONARY DISTANCES

EVOLUTIONARY DISTANCES EVOLUTIONARY DISTANCES FROM STRINGS TO TREES Luca Bortolussi 1 1 Dipartimento di Matematica ed Informatica Università degli studi di Trieste luca@dmi.units.it Trieste, 14 th November 2007 OUTLINE 1 STRINGS:

More information

Frequency Spectra and Inference in Population Genetics

Frequency Spectra and Inference in Population Genetics Frequency Spectra and Inference in Population Genetics Although coalescent models have come to play a central role in population genetics, there are some situations where genealogies may not lead to efficient

More information

ON COMPOUND POISSON POPULATION MODELS

ON COMPOUND POISSON POPULATION MODELS ON COMPOUND POISSON POPULATION MODELS Martin Möhle, University of Tübingen (joint work with Thierry Huillet, Université de Cergy-Pontoise) Workshop on Probability, Population Genetics and Evolution Centre

More information

Scalable Algorithms for Population Genomic Inference

Scalable Algorithms for Population Genomic Inference Scalable Algorithms for Population Genomic Inference Sara Sheehan Electrical Engineering and Computer Sciences University of California at Berkeley Technical Report No. UCB/EECS-2015-94 http://www.eecs.berkeley.edu/pubs/techrpts/2015/eecs-2015-94.html

More information

(Write your name on every page. One point will be deducted for every page without your name!)

(Write your name on every page. One point will be deducted for every page without your name!) POPULATION GENETICS AND MICROEVOLUTIONARY THEORY FINAL EXAMINATION (Write your name on every page. One point will be deducted for every page without your name!) 1. Briefly define (5 points each): a) Average

More information

The Combinatorial Interpretation of Formulas in Coalescent Theory

The Combinatorial Interpretation of Formulas in Coalescent Theory The Combinatorial Interpretation of Formulas in Coalescent Theory John L. Spouge National Center for Biotechnology Information NLM, NIH, DHHS spouge@ncbi.nlm.nih.gov Bldg. A, Rm. N 0 NCBI, NLM, NIH Bethesda

More information

Lecture 18 : Ewens sampling formula

Lecture 18 : Ewens sampling formula Lecture 8 : Ewens sampling formula MATH85K - Spring 00 Lecturer: Sebastien Roch References: [Dur08, Chapter.3]. Previous class In the previous lecture, we introduced Kingman s coalescent as a limit of

More information

Demographic Inference with Coalescent Hidden Markov Model

Demographic Inference with Coalescent Hidden Markov Model Demographic Inference with Coalescent Hidden Markov Model Jade Y. Cheng Thomas Mailund Bioinformatics Research Centre Aarhus University Denmark The Thirteenth Asia Pacific Bioinformatics Conference HsinChu,

More information

Expected coalescence times and segregating sites in a model of glacial cycles

Expected coalescence times and segregating sites in a model of glacial cycles F.F. Jesus et al. 466 Expected coalescence times and segregating sites in a model of glacial cycles F.F. Jesus 1, J.F. Wilkins 2, V.N. Solferini 1 and J. Wakeley 3 1 Departamento de Genética e Evolução,

More information

The mathematical challenge. Evolution in a spatial continuum. The mathematical challenge. Other recruits... The mathematical challenge

The mathematical challenge. Evolution in a spatial continuum. The mathematical challenge. Other recruits... The mathematical challenge The mathematical challenge What is the relative importance of mutation, selection, random drift and population subdivision for standing genetic variation? Evolution in a spatial continuum Al lison Etheridge

More information

Reading for Lecture 13 Release v10

Reading for Lecture 13 Release v10 Reading for Lecture 13 Release v10 Christopher Lee November 15, 2011 Contents 1 Evolutionary Trees i 1.1 Evolution as a Markov Process...................................... ii 1.2 Rooted vs. Unrooted Trees........................................

More information

Population Structure

Population Structure Ch 4: Population Subdivision Population Structure v most natural populations exist across a landscape (or seascape) that is more or less divided into areas of suitable habitat v to the extent that populations

More information

THERE has been interest in analyzing population genomic

THERE has been interest in analyzing population genomic INVESTIGATION The Lengths of Admixture Tracts Mason Liang 1 and Rasmus Nielsen University of California, Berkeley, California 94108 ABSTRACT The distribution of admixture tract lengths has received considerable

More information

Population Genetics: a tutorial

Population Genetics: a tutorial : a tutorial Institute for Science and Technology Austria ThRaSh 2014 provides the basic mathematical foundation of evolutionary theory allows a better understanding of experiments allows the development

More information

O 3 O 4 O 5. q 3. q 4. Transition

O 3 O 4 O 5. q 3. q 4. Transition Hidden Markov Models Hidden Markov models (HMM) were developed in the early part of the 1970 s and at that time mostly applied in the area of computerized speech recognition. They are first described in

More information

The problem Lineage model Examples. The lineage model

The problem Lineage model Examples. The lineage model The lineage model A Bayesian approach to inferring community structure and evolutionary history from whole-genome metagenomic data Jack O Brien Bowdoin College with Daniel Falush and Xavier Didelot Cambridge,

More information

ESTIMATION of recombination fractions using ped- ber of pairwise differences (Hudson 1987; Wakeley

ESTIMATION of recombination fractions using ped- ber of pairwise differences (Hudson 1987; Wakeley Copyright 2001 by the Genetics Society of America Estimating Recombination Rates From Population Genetic Data Paul Fearnhead and Peter Donnelly Department of Statistics, University of Oxford, Oxford, OX1

More information

The nested Kingman coalescent: speed of coming down from infinity. by Jason Schweinsberg (University of California at San Diego)

The nested Kingman coalescent: speed of coming down from infinity. by Jason Schweinsberg (University of California at San Diego) The nested Kingman coalescent: speed of coming down from infinity by Jason Schweinsberg (University of California at San Diego) Joint work with: Airam Blancas Benítez (Goethe Universität Frankfurt) Tim

More information

Evolution in a spatial continuum

Evolution in a spatial continuum Evolution in a spatial continuum Drift, draft and structure Alison Etheridge University of Oxford Joint work with Nick Barton (Edinburgh) and Tom Kurtz (Wisconsin) New York, Sept. 2007 p.1 Kingman s Coalescent

More information

Phylogenomics of closely related species and individuals

Phylogenomics of closely related species and individuals Phylogenomics of closely related species and individuals Matthew Rasmussen Siepel lab, Cornell University In collaboration with Manolis Kellis, MIT CSAIL February, 2013 Short time scales 1kyr-1myrs Long

More information

Neutral Theory of Molecular Evolution

Neutral Theory of Molecular Evolution Neutral Theory of Molecular Evolution Kimura Nature (968) 7:64-66 King and Jukes Science (969) 64:788-798 (Non-Darwinian Evolution) Neutral Theory of Molecular Evolution Describes the source of variation

More information

arxiv: v1 [q-bio.pe] 14 Jan 2016

arxiv: v1 [q-bio.pe] 14 Jan 2016 Efficient Maximum-Likelihood Inference For The Isolation-With-Initial-Migration Model With Potentially Asymmetric Gene Flow Rui J. Costa and Hilde Wilkinson-Herbots arxiv:1601.03684v1 [q-bio.pe] 14 Jan

More information

Modelling Genetic Variations with Fragmentation-Coagulation Processes

Modelling Genetic Variations with Fragmentation-Coagulation Processes Modelling Genetic Variations with Fragmentation-Coagulation Processes Yee Whye Teh, Charles Blundell, Lloyd Elliott Gatsby Computational Neuroscience Unit, UCL Genetic Variations in Populations Inferring

More information

Full Likelihood Inference from the Site Frequency Spectrum based on the Optimal Tree Resolution

Full Likelihood Inference from the Site Frequency Spectrum based on the Optimal Tree Resolution Full Likelihood Inference from the Site Frequency Spectrum based on the Optimal Tree Resolution Raazesh Sainudiin 1, Amandine Véber 2 March 3, 2018 1 Department of Mathematics, Uppsala University, Uppsala,

More information

Using phylogenetics to estimate species divergence times... Basics and basic issues for Bayesian inference of divergence times (plus some digression)

Using phylogenetics to estimate species divergence times... Basics and basic issues for Bayesian inference of divergence times (plus some digression) Using phylogenetics to estimate species divergence times... More accurately... Basics and basic issues for Bayesian inference of divergence times (plus some digression) "A comparison of the structures

More information

Estimating effective population size from samples of sequences: inefficiency of pairwise and segregating sites as compared to phylogenetic estimates

Estimating effective population size from samples of sequences: inefficiency of pairwise and segregating sites as compared to phylogenetic estimates Estimating effective population size from samples of sequences: inefficiency of pairwise and segregating sites as compared to phylogenetic estimates JOSEPH FELSENSTEIN Department of Genetics SK-50, University

More information

How should we organize the diversity of animal life?

How should we organize the diversity of animal life? How should we organize the diversity of animal life? The difference between Taxonomy Linneaus, and Cladistics Darwin What are phylogenies? How do we read them? How do we estimate them? Classification (Taxonomy)

More information

An Introduction to Bioinformatics Algorithms Hidden Markov Models

An Introduction to Bioinformatics Algorithms  Hidden Markov Models Hidden Markov Models Hidden Markov Models Outline CG-islands The Fair Bet Casino Hidden Markov Model Decoding Algorithm Forward-Backward Algorithm Profile HMMs HMM Parameter Estimation Viterbi training

More information

Lecture 14 Chapter 11 Biology 5865 Conservation Biology. Problems of Small Populations Population Viability Analysis

Lecture 14 Chapter 11 Biology 5865 Conservation Biology. Problems of Small Populations Population Viability Analysis Lecture 14 Chapter 11 Biology 5865 Conservation Biology Problems of Small Populations Population Viability Analysis Minimum Viable Population (MVP) Schaffer (1981) MVP- A minimum viable population for

More information

URN MODELS: the Ewens Sampling Lemma

URN MODELS: the Ewens Sampling Lemma Department of Computer Science Brown University, Providence sorin@cs.brown.edu October 3, 2014 1 2 3 4 Mutation Mutation: typical values for parameters Equilibrium Probability of fixation 5 6 Ewens Sampling

More information

An introduction to Approximate Bayesian Computation methods

An introduction to Approximate Bayesian Computation methods An introduction to Approximate Bayesian Computation methods M.E. Castellanos maria.castellanos@urjc.es (from several works with S. Cabras, E. Ruli and O. Ratmann) Valencia, January 28, 2015 Valencia Bayesian

More information

Introduction to Natural Selection. Ryan Hernandez Tim O Connor

Introduction to Natural Selection. Ryan Hernandez Tim O Connor Introduction to Natural Selection Ryan Hernandez Tim O Connor 1 Goals Learn about the population genetics of natural selection How to write a simple simulation with natural selection 2 Basic Biology genome

More information

Processes of Evolution

Processes of Evolution 15 Processes of Evolution Forces of Evolution Concept 15.4 Selection Can Be Stabilizing, Directional, or Disruptive Natural selection can act on quantitative traits in three ways: Stabilizing selection

More information

Genetic Drift in Human Evolution

Genetic Drift in Human Evolution Genetic Drift in Human Evolution (Part 2 of 2) 1 Ecology and Evolutionary Biology Center for Computational Molecular Biology Brown University Outline Introduction to genetic drift Modeling genetic drift

More information

Theoretical Population Biology

Theoretical Population Biology Theoretical Population Biology 87 (013) 6 74 Contents lists available at SciVerse ScienceDirect Theoretical Population Biology journal homepage: www.elsevier.com/locate/tpb Genotype imputation in a coalescent

More information

I of a gene sampled from a randomly mating popdation,

I of a gene sampled from a randomly mating popdation, Copyright 0 1987 by the Genetics Society of America Average Number of Nucleotide Differences in a From a Single Subpopulation: A Test for Population Subdivision Curtis Strobeck Department of Zoology, University

More information

Introduction to Advanced Population Genetics

Introduction to Advanced Population Genetics Introduction to Advanced Population Genetics Learning Objectives Describe the basic model of human evolutionary history Describe the key evolutionary forces How demography can influence the site frequency

More information

Estimating Evolutionary Trees. Phylogenetic Methods

Estimating Evolutionary Trees. Phylogenetic Methods Estimating Evolutionary Trees v if the data are consistent with infinite sites then all methods should yield the same tree v it gets more complicated when there is homoplasy, i.e., parallel or convergent

More information

Populations in statistical genetics

Populations in statistical genetics Populations in statistical genetics What are they, and how can we infer them from whole genome data? Daniel Lawson Heilbronn Institute, University of Bristol www.paintmychromosomes.com Work with: January

More information

How Many Subpopulations is Too Many? Exponential Lower Bounds for Inferring Population Histories

How Many Subpopulations is Too Many? Exponential Lower Bounds for Inferring Population Histories How Many Subpopulations is Too Many? Exponential Lower Bounds for Inferring Population Histories Younhun Kim Frederic Koehler Ankur Moitra Elchanan Mossel Govind Ramnarayan arxiv:8.0377v [q-bio.pe] 7 Nov

More information

Phylogenies & Classifying species (AKA Cladistics & Taxonomy) What are phylogenies & cladograms? How do we read them? How do we estimate them?

Phylogenies & Classifying species (AKA Cladistics & Taxonomy) What are phylogenies & cladograms? How do we read them? How do we estimate them? Phylogenies & Classifying species (AKA Cladistics & Taxonomy) What are phylogenies & cladograms? How do we read them? How do we estimate them? Carolus Linneaus:Systema Naturae (1735) Swedish botanist &

More information

Mathematical Population Genetics II

Mathematical Population Genetics II Mathematical Population Genetics II Lecture Notes Joachim Hermisson March 20, 2015 University of Vienna Mathematics Department Oskar-Morgenstern-Platz 1 1090 Vienna, Austria Copyright (c) 2013/14/15 Joachim

More information

CSci 8980: Advanced Topics in Graphical Models Analysis of Genetic Variation

CSci 8980: Advanced Topics in Graphical Models Analysis of Genetic Variation CSci 8980: Advanced Topics in Graphical Models Analysis of Genetic Variation Instructor: Arindam Banerjee November 26, 2007 Genetic Polymorphism Single nucleotide polymorphism (SNP) Genetic Polymorphism

More information

Mathematical Population Genetics II

Mathematical Population Genetics II Mathematical Population Genetics II Lecture Notes Joachim Hermisson June 9, 2018 University of Vienna Mathematics Department Oskar-Morgenstern-Platz 1 1090 Vienna, Austria Copyright (c) 2013/14/15/18 Joachim

More information

Application of a time-dependent coalescence process for inferring the history of population size changes from DNA sequence data

Application of a time-dependent coalescence process for inferring the history of population size changes from DNA sequence data Proc. Natl. Acad. Sci. USA Vol. 95, pp. 5456 546, May 998 Statistics Application of a time-dependent coalescence process for inferring the history of population size changes from DNA sequence data ANDRZEJ

More information

The phenotype of this worm is wild type. When both genes are mutant: The phenotype of this worm is double mutant Dpy and Unc phenotype.

The phenotype of this worm is wild type. When both genes are mutant: The phenotype of this worm is double mutant Dpy and Unc phenotype. Series 2: Cross Diagrams - Complementation There are two alleles for each trait in a diploid organism In C. elegans gene symbols are ALWAYS italicized. To represent two different genes on the same chromosome:

More information

Modern Evolutionary Classification. Section 18-2 pgs

Modern Evolutionary Classification. Section 18-2 pgs Modern Evolutionary Classification Section 18-2 pgs 451-455 Modern Evolutionary Classification In a sense, organisms determine who belongs to their species by choosing with whom they will mate. Taxonomic

More information

ms a program for generating samples under neutral models

ms a program for generating samples under neutral models ms a program for generating samples under neutral models Richard R. Hudson December 11, 2004 This document describes how to use ms, a program to generate samples under a variety of neutral models. The

More information

Effective population size and patterns of molecular evolution and variation

Effective population size and patterns of molecular evolution and variation FunDamental concepts in genetics Effective population size and patterns of molecular evolution and variation Brian Charlesworth Abstract The effective size of a population,, determines the rate of change

More information

Formalizing the gene centered view of evolution

Formalizing the gene centered view of evolution Chapter 1 Formalizing the gene centered view of evolution Yaneer Bar-Yam and Hiroki Sayama New England Complex Systems Institute 24 Mt. Auburn St., Cambridge, MA 02138, USA yaneer@necsi.org / sayama@necsi.org

More information

Alignment Algorithms. Alignment Algorithms

Alignment Algorithms. Alignment Algorithms Midterm Results Big improvement over scores from the previous two years. Since this class grade is based on the previous years curve, that means this class will get higher grades than the previous years.

More information

A Phylogenetic Network Construction due to Constrained Recombination

A Phylogenetic Network Construction due to Constrained Recombination A Phylogenetic Network Construction due to Constrained Recombination Mohd. Abdul Hai Zahid Research Scholar Research Supervisors: Dr. R.C. Joshi Dr. Ankush Mittal Department of Electronics and Computer

More information

Evolutionary Computation

Evolutionary Computation Evolutionary Computation - Computational procedures patterned after biological evolution. - Search procedure that probabilistically applies search operators to set of points in the search space. - Lamarck

More information

Learning ancestral genetic processes using nonparametric Bayesian models

Learning ancestral genetic processes using nonparametric Bayesian models Learning ancestral genetic processes using nonparametric Bayesian models Kyung-Ah Sohn October 31, 2011 Committee Members: Eric P. Xing, Chair Zoubin Ghahramani Russell Schwartz Kathryn Roeder Matthew

More information

Stationary Distribution of the Linkage Disequilibrium Coefficient r 2

Stationary Distribution of the Linkage Disequilibrium Coefficient r 2 Stationary Distribution of the Linkage Disequilibrium Coefficient r 2 Wei Zhang, Jing Liu, Rachel Fewster and Jesse Goodman Department of Statistics, The University of Auckland December 1, 2015 Overview

More information

Runaway. demogenetic model for sexual selection. Louise Chevalier. Jacques Labonne

Runaway. demogenetic model for sexual selection. Louise Chevalier. Jacques Labonne Runaway demogenetic model for sexual selection Louise Chevalier Master 2 thesis UPMC, Specialization Oceanography and Marine Environments Jacques Labonne UMR Ecobiop INRA - National Institute for Agronomic

More information

Introduction to Artificial Intelligence (AI)

Introduction to Artificial Intelligence (AI) Introduction to Artificial Intelligence (AI) Computer Science cpsc502, Lecture 10 Oct, 13, 2011 CPSC 502, Lecture 10 Slide 1 Today Oct 13 Inference in HMMs More on Robot Localization CPSC 502, Lecture

More information

HIDDEN MARKOV MODELS

HIDDEN MARKOV MODELS HIDDEN MARKOV MODELS Outline CG-islands The Fair Bet Casino Hidden Markov Model Decoding Algorithm Forward-Backward Algorithm Profile HMMs HMM Parameter Estimation Viterbi training Baum-Welch algorithm

More information

Supplementary Material to Full Likelihood Inference from the Site Frequency Spectrum based on the Optimal Tree Resolution

Supplementary Material to Full Likelihood Inference from the Site Frequency Spectrum based on the Optimal Tree Resolution Supplementary Material to Full Likelihood Inference from the Site Frequency Spectrum based on the Optimal Tree Resolution Raazesh Sainudiin, Amandine Véber 1 Pseudo-code First, the global process MakeHistory

More information

Genetics. Metapopulations. Dept. of Forest & Wildlife Ecology, UW Madison

Genetics. Metapopulations. Dept. of Forest & Wildlife Ecology, UW Madison Genetics & Metapopulations Dr Stacie J Robinson Dr. Stacie J. Robinson Dept. of Forest & Wildlife Ecology, UW Madison Robinson ~ UW SJR OUTLINE Metapopulation impacts on evolutionary processes Metapopulation

More information

Evolution of cooperation in finite populations. Sabin Lessard Université de Montréal

Evolution of cooperation in finite populations. Sabin Lessard Université de Montréal Evolution of cooperation in finite populations Sabin Lessard Université de Montréal Title, contents and acknowledgements 2011 AMS Short Course on Evolutionary Game Dynamics Contents 1. Examples of cooperation

More information

Surfing genes. On the fate of neutral mutations in a spreading population

Surfing genes. On the fate of neutral mutations in a spreading population Surfing genes On the fate of neutral mutations in a spreading population Oskar Hallatschek David Nelson Harvard University ohallats@physics.harvard.edu Genetic impact of range expansions Population expansions

More information

THE analysis of population subdivision has been a ble to transform an estimate of F ST into an estimate of

THE analysis of population subdivision has been a ble to transform an estimate of F ST into an estimate of Copyright 2001 by the Genetics Society of America Distinguishing Migration From Isolation: A Markov Chain Monte Carlo Approach Rasmus Nielsen and John Wakeley Department of Organismic and Evolutionary

More information

Evolutionary Models. Evolutionary Models

Evolutionary Models. Evolutionary Models Edit Operators In standard pairwise alignment, what are the allowed edit operators that transform one sequence into the other? Describe how each of these edit operations are represented on a sequence alignment

More information

CLOSED-FORM ASYMPTOTIC SAMPLING DISTRIBUTIONS UNDER THE COALESCENT WITH RECOMBINATION FOR AN ARBITRARY NUMBER OF LOCI

CLOSED-FORM ASYMPTOTIC SAMPLING DISTRIBUTIONS UNDER THE COALESCENT WITH RECOMBINATION FOR AN ARBITRARY NUMBER OF LOCI Adv. Appl. Prob. 44, 391 407 (01) Printed in Northern Ireland Applied Probability Trust 01 CLOSED-FORM ASYMPTOTIC SAMPLING DISTRIBUTIONS UNDER THE COALESCENT WITH RECOMBINATION FOR AN ARBITRARY NUMBER

More information

Metapopulation models for historical inference

Metapopulation models for historical inference Molecular Ecology (2004) 13, 865 875 doi: 10.1111/j.1365-294X.2004.02086.x Metapopulation models for historical inference Blackwell Publishing, Ltd. JOHN WAKELEY Department of Organismic and Evolutionary

More information

Solutions to Even-Numbered Exercises to accompany An Introduction to Population Genetics: Theory and Applications Rasmus Nielsen Montgomery Slatkin

Solutions to Even-Numbered Exercises to accompany An Introduction to Population Genetics: Theory and Applications Rasmus Nielsen Montgomery Slatkin Solutions to Even-Numbered Exercises to accompany An Introduction to Population Genetics: Theory and Applications Rasmus Nielsen Montgomery Slatkin CHAPTER 1 1.2 The expected homozygosity, given allele

More information

On the inadmissibility of Watterson s estimate

On the inadmissibility of Watterson s estimate On the inadmissibility of Watterson s estimate A. Futschik F. Gach Institute of Statistics and Decision Support Systems University of Vienna ISDS-Kolloquium, 2007 Outline 1 Motivation 2 Inadmissibility

More information

CS1820 Notes. hgupta1, kjline, smechery. April 3-April 5. output: plausible Ancestral Recombination Graph (ARG)

CS1820 Notes. hgupta1, kjline, smechery. April 3-April 5. output: plausible Ancestral Recombination Graph (ARG) CS1820 Notes hgupta1, kjline, smechery April 3-April 5 April 3 Notes 1 Minichiello-Durbin Algorithm input: set of sequences output: plausible Ancestral Recombination Graph (ARG) note: the optimal ARG is

More information

Introduction to Wright-Fisher Simulations. Ryan Hernandez

Introduction to Wright-Fisher Simulations. Ryan Hernandez Introduction to Wright-Fisher Simulations Ryan Hernandez 1 Goals Simulate the standard neutral model, demographic effects, and natural selection Start with single sites, and build in multiple sites 2 Hardy-Weinberg

More information

The Impact of Sampling Schemes on the Site Frequency Spectrum in Nonequilibrium Subdivided Populations

The Impact of Sampling Schemes on the Site Frequency Spectrum in Nonequilibrium Subdivided Populations Copyright Ó 2009 by the Genetics Society of America DOI: 10.1534/genetics.108.094904 The Impact of Sampling Schemes on the Site Frequency Spectrum in Nonequilibrium Subdivided Populations Thomas Städler,*,1

More information

genome a specific characteristic that varies from one individual to another gene the passing of traits from one generation to the next

genome a specific characteristic that varies from one individual to another gene the passing of traits from one generation to the next genetics the study of heredity heredity sequence of DNA that codes for a protein and thus determines a trait genome a specific characteristic that varies from one individual to another gene trait the passing

More information

An introduction to mathematical modeling of the genealogical process of genes

An introduction to mathematical modeling of the genealogical process of genes An introduction to mathematical modeling of the genealogical process of genes Rikard Hellman Kandidatuppsats i matematisk statistik Bachelor Thesis in Mathematical Statistics Kandidatuppsats 2009:3 Matematisk

More information

The Structure of Genealogies in the Presence of Purifying Selection: a "Fitness-Class Coalescent"

The Structure of Genealogies in the Presence of Purifying Selection: a Fitness-Class Coalescent The Structure of Genealogies in the Presence of Purifying Selection: a "Fitness-Class Coalescent" The Harvard community has made this article openly available. Please share how this access benefits you.

More information

Wright-Fisher Models, Approximations, and Minimum Increments of Evolution

Wright-Fisher Models, Approximations, and Minimum Increments of Evolution Wright-Fisher Models, Approximations, and Minimum Increments of Evolution William H. Press The University of Texas at Austin January 10, 2011 1 Introduction Wright-Fisher models [1] are idealized models

More information

Probability and Time: Hidden Markov Models (HMMs)

Probability and Time: Hidden Markov Models (HMMs) Probability and Time: Hidden Markov Models (HMMs) Computer Science cpsc322, Lecture 32 (Textbook Chpt 6.5.2) Nov, 25, 2013 CPSC 322, Lecture 32 Slide 1 Lecture Overview Recap Markov Models Markov Chain

More information

Modelling Linkage Disequilibrium, And Identifying Recombination Hotspots Using SNP Data

Modelling Linkage Disequilibrium, And Identifying Recombination Hotspots Using SNP Data Modelling Linkage Disequilibrium, And Identifying Recombination Hotspots Using SNP Data Na Li and Matthew Stephens July 25, 2003 Department of Biostatistics, University of Washington, Seattle, WA 98195

More information