Coordinated Conditional Simulation with SLINK and SUP of Many Markers Linked or Associated to a Trait in Large Pedigrees

Simulation of genotypes in pedigrees is an important tool to evaluate the power of a linkage or an association study and to assess the empirical significance of results. SLINK is a widely-used package for pedigree simulations, but its implementation has not previously been described in a published paper. SLINK was initially derived from the LINKAGE programs. Over the 20 years since its release, SLINK has been modified to incorporate faster algorithms, notably from the linkage analysis package FASTLINK, also derived from LINKAGE. While SLINK can simulate genotypes on pedigrees of high complexity, one limitation of SLINK, as with most methods based on peeling algorithms to evaluate pedigree likelihoods, is the small number of linked markers that can be generated. The software package SUP includes an elegant wrapper for SLINK that circumvents the limitation on number of markers by using descent markers generated by SLINK to simulate a much larger number of markers on the same chromosome, linked and possibly associated with a trait locus. We have released new coordinated versions of SLINK (3.0; available from http://watson.hgen.pitt.edu) and SUP (v090804; available from http://mlemire.freeshell.org/software or http://watson.hgen.pitt.edu) that integrate the two software packages. Thereby, we have removed some of the previous limitations on the joint functionality of the programs, such as the number of founders in a pedigree. We review the history of SLINK and describe how SLINK and SUP are now coordinated to permit the simulation of large numbers of markers linked and possibly associated with a trait in large pedigrees.

[1]  Oliver A. Ryder,et al.  Pedigree analysis by computer simulation , 1986 .

[2]  M. King,et al.  Possible linkage of the estrogen receptor gene to breast cancer in a family with late-onset disease. , 1991, American journal of human genetics.

[3]  J. Ott Computer-simulation methods in human linkage analysis. , 1989, Proceedings of the National Academy of Sciences of the United States of America.

[4]  A. Schäffer,et al.  Linkage of autosomal-dominant common variable immunodeficiency to chromosome 4q , 2006, European Journal of Human Genetics.

[5]  J. Ott,et al.  Strategies for multilocus linkage analysis in humans. , 1984, Proceedings of the National Academy of Sciences of the United States of America.

[6]  T. Hudson,et al.  Transmission-ratio distortion and allele sharing in affected sib pairs: a new linkage statistic with reduced bias, with application to chromosome 6q25.3. , 2004, American journal of human genetics.

[7]  W. Fung,et al.  Detection of parent‐of‐origin effects using general pedigree data , 2010, Genetic epidemiology.

[8]  Daniel F. Gudbjartsson,et al.  Allegro, a new computer program for multipoint linkage analysis , 2000, Nature genetics.

[9]  S E Hodge,et al.  The power to detect linkage in complex disease by means of simple LOD-score analyses. , 1998, American journal of human genetics.

[10]  A A Schäffer,et al.  Faster sequential genetic linkage computations. , 1993, American journal of human genetics.

[11]  Mike Schmidt,et al.  Statistical Applications in Genetics and Molecular Biology Extension of the SIMLA Package for Generating Pedigrees with Complex Inheritance Patterns : Environmental Covariates , Gene-Gene and Gene-Environment Interaction , 2011 .

[12]  S. Leal,et al.  SimPed: A Simulation Program to Generate Haplotype and Genotype Data for Pedigree Structures , 2005, Human Heredity.

[13]  G. Abecasis,et al.  Merlin—rapid analysis of dense genetic maps using sparse gene flow trees , 2002, Nature Genetics.

[14]  T. Burns,et al.  Case-Deletion Diagnostics for Maximum Likelihood Multipoint Quantitative Trait Locus Linkage Analysis , 2009, Human Heredity.

[15]  Dan Geiger,et al.  Exact genetic linkage computations for general pedigrees , 2002, ISMB.

[16]  Mathieu Lemire,et al.  SUP: an extension to SLINK to allow a larger number of marker loci to be simulated in pedigrees conditional on trait values , 2006, BMC Genetics.

[17]  R. Elston,et al.  A general model for the genetic analysis of pedigree data. , 1971, Human heredity.

[18]  M. Boehnke,et al.  Estimating the power of a proposed linkage study for a complex genetic trait. , 1989, American journal of human genetics.

[19]  M. Leboyer,et al.  Performance of linkage analysis under misclassification error when the genetic model is unknown , 1989, Genetic epidemiology.

[20]  M Speer,et al.  Chromosome‐based method for rapid computer simulation in human genetic linkage analysis , 1993, Genetic epidemiology.

[21]  J. Ott Analysis of Human Genetic Linkage , 1985 .

[22]  Derek Gordon,et al.  When a Case Is Not a Case: Effects of Phenotype Misclassification on Power and Sample Size Requirements for the Transmission Disequilibrium Test with Affected Child Trios , 2009, Human Heredity.

[23]  W. Knowler,et al.  Design and Analysis of Genetic Association Studies to Finely Map a Locus Identified by Linkage Analysis: Assessment of the Extent to Which an Association Can Account for the Linkage , 2007, Annals of human genetics.