Fast and flexible simulation of DNA sequence data.
Simulation of genomic sequences under the coalescent with recombination has conventionally been impractical for regions beyond tens of megabases. This work presents an algorithm, implemented as the program MaCS (Markovian Coalescent Simulator), that can efficiently simulate haplotypes under any arbitrary model of population history. We present several metrics comparing the performance of MaCS with other available simulation programs. Practical usage of MaCS is demonstrated through a comparison of measures of linkage disequilibrium between generated program output and real genotype data from populations considered to be structured.
Chen
Gary K
GK
Department of Preventive Medicine, University of Southern California, Los Angeles, California 90033, USA.
Marjoram
Paul
P
Wall
Jeffrey D
JD
D000465:Algorithms; D003198:Computer Simulation; D004247:DNA / Q000235:genetics*; D055106:Genome-Wide Association Study / Q000706:statistics & numerical data; D023281:Genomics / Q000706:statistics & numerical data; D006801:Humans; D015810:Linkage Disequilibrium; D008390:Markov Chains; D020641:Polymorphism, Single Nucleotide; D011995:Recombination, Genetic; D017422:Sequence Analysis, DNA / Q000706:statistics & numerical data*; D012984:Software
D004247:DNA