Non-equilibrium allele frequency spectra via spectral methods

Sergio Lukić, Jody Hey, Kevin Chen

Research output: Contribution to journalArticlepeer-review

30 Scopus citations

Abstract

A major challenge in the analysis of population genomics data consists of isolating signatures of natural selection from background noise caused by random drift and gene flow. Analyses of massive amounts of data from many related populations require high-performance algorithms to determine the likelihood of different demographic scenarios that could have shaped the observed neutral single nucleotide polymorphism (SNP) allele frequency spectrum. In many areas of applied mathematics, Fourier Transforms and Spectral Methods are firmly established tools to analyze spectra of signals and model their dynamics as solutions of certain Partial Differential Equations (PDEs). When spectral methods are applicable, they have excellent error properties and are the fastest possible in high dimension; see Press et al. (2007). In this paper we present an explicit numerical solution, using spectral methods, to the forward Kolmogorov equations for a Wright-Fisher process with migration of K populations, influx of mutations, and multiple population splitting events.

Original languageEnglish (US)
Pages (from-to)203-219
Number of pages17
JournalTheoretical Population Biology
Volume79
Issue number4
DOIs
StatePublished - Jun 2011

All Science Journal Classification (ASJC) codes

  • Ecology, Evolution, Behavior and Systematics

Keywords

  • Diffusion theory
  • Kolmogorov forward equation
  • Population genetics
  • Spectral methods

Fingerprint

Dive into the research topics of 'Non-equilibrium allele frequency spectra via spectral methods'. Together they form a unique fingerprint.

Cite this