Marker-marker interaction analyses in case-control sample were performed by mixed model backward logistic regression (Zuo et al 2007). The phenotypes served as the dependent variable and sex, age and genotypes were covariates.