Characterization of genomic imprinting effects and patterns with parametric accelerated failure time model

Genomic imprinting, a genetic phenomenon of non-equivalent allele expression that depends on parental origins, has been ubiquitously observed in nature. It does not only control the traits of growth and development but also may be responsible for survival traits. Based on the accelerated failure time model, we construct a general parametric model for mapping the imprinted QTL (iQTL). Within the framework of interval mapping, maximum likelihood estimation of iQTL parameters is implemented via EM algorithm. The imprinting patterns of the detected iQTL are statistically tested according to a series of null hypotheses. BIC model selection criterion is employed to choose an optimal baseline hazard function with maximum likelihood and parsimonious parameters. Simulations are used to validate the proposed mapping procedure. A published dataset from a mouse model system was used to illustrate the proposed framework. Results show that among the five commonly used survival distributions, Log-logistic distribution is the optimal baseline hazard function for mapping QTL of hyperoxic acute lung injury (HALI) survival; under the log-logistic distribution, four QTLs were identified, in which only one QTL was inherited in Mendelian fashion, whereas others were imprinted in different imprinting patterns.