Mapping Temporally Varying Quantitative Trait Loci in Time-to-Failure Experiments

Existing methods for mapping quantitative trait loci (QTL) in time-to-failure experiments assume that the QTL effect is constant over the course of the study. This assumption may be violated when the gene(s) underlying the QTL are up- or downregulated on a biologically meaningful timescale. In such situations, models that assume a constant effect can fail to detect QTL in a whole-genome scan. To investigate this possibility, we utilize an extension of the Cox model (EC model) within an interval-mapping framework. In its simplest form, this model assumes that the QTL effect changes at some time point t0 and follows a linear function before and after this change point. The approximate time point at which this change occurs is estimated. Using simulated and real data, we compare the mapping performance of the EC model to the Cox proportional hazards (CPH) model, which explicitly assumes a constant effect. The results show that the EC model detects time-dependent QTL, which the CPH model fails to detect. At the same time, the EC model recovers all of the QTL the CPH model detects. We conclude that potentially important QTL may be missed if their time-dependent effects are not accounted for.

[1]  D. Schoenfeld The asymptotic properties of nonparametric tests for comparing survival distributions , 1981 .

[2]  Karl W Broman,et al.  Mapping quantitative trait loci in the case of a spike in the phenotype distribution. , 2003, Genetics.

[3]  R. Wu,et al.  Functional mapping — how to map and study the genetic architecture of dynamic complex traits , 2006, Nature Reviews Genetics.

[4]  P. Grambsch,et al.  Proportional hazards tests and diagnostics based on weighted residuals , 1994 .

[5]  Jun Yan Survival Analysis: Techniques for Censored and Truncated Data , 2004 .

[6]  D. Lin,et al.  Semiparametric Methods for Mapping Quantitative Trait Loci with Censored Data , 2005, Biometrics.

[7]  E. Lehmann,et al.  Nonparametrics: Statistical Methods Based on Ranks , 1976 .

[8]  G. Casella,et al.  Mapping epigenetic quantitative trait loci (QTL) altering a developmental trajectory. , 2002, Genome.

[9]  M. Schemper Cox Analysis of Survival Data with Non‐Proportional Hazard Functions , 1992 .

[10]  M. Tableman,et al.  Survival Analysis Using S : Analysis of Time-to-Event Data , 2003 .

[11]  G. Mcclearn,et al.  Genetic architecture of fast- and slow-twitch skeletal muscle weight in 200-day-old mice of the C57BL/6J and DBA/2J lineage. , 2003, Physiological genomics.

[12]  R. Doerge,et al.  Empirical threshold values for quantitative trait mapping. , 1994, Genetics.

[13]  H. Toutenburg,et al.  Lehmann, E. L., Nonparametrics: Statistical Methods Based on Ranks, San Francisco. Holden‐Day, Inc., 1975. 480 S., $ 22.95 . , 1977 .

[14]  D. Schoenfeld,et al.  Sample-size formula for the proportional-hazards regression model. , 1983, Biometrics.

[15]  V. Ducrocq,et al.  Interval mapping methods for detecting QTL affecting survival and time-to-event phenotypes. , 2005, Genetical research.

[16]  H. Snoeck,et al.  Genetically determined variation in the number of phenotypically defined hematopoietic progenitor and stem cells and in their response to early-acting cytokines. , 2002, Blood.

[17]  Michal Abrahamowicz,et al.  Time-Dependent Hazard Ratio: Modeling and Hypothesis Testing with Application in Lupus Nephritis , 1996 .

[18]  G. Vogler,et al.  QTL influencing baseline hematocrit in the C57BL/6J and DBA/2J lineage: age-related effects , 2006, Mammalian Genome.

[19]  G. Casella,et al.  Functional Mapping of Quantitative Trait Loci Underlying Growth Trajectories Using a Transform‐Both‐Sides Logistic Model , 2004, Biometrics.

[20]  W. Li,et al.  Mapping of quantitative trait loci based on growth models , 2002, Theoretical and Applied Genetics.

[21]  Fei Zou,et al.  Mapping Quantitative Trait Loci With Censored Observations , 2004, Genetics.

[22]  R. Mosher,et al.  Multigenic control of Listeria monocytogenes susceptibility in mice , 2001, Nature Genetics.

[23]  David A. Schoenfeld,et al.  Partial residuals for the proportional hazards regression model , 1982 .

[24]  G. Casella,et al.  Functional mapping of quantitative trait loci underlying the character process: a theoretical framework. , 2002, Genetics.

[25]  J. Cheverud,et al.  Quantitative trait loci for murine growth. , 1996, Genetics.

[26]  David P. Harrington,et al.  Modified Kolmogorov-Smirnov Test Procedures with Application to Arbitrarily Right-Censored Data , 1980 .