To explore the independent effects of various clinical and genetic risk factors contributing to pain susceptibility in PD patients, we applied a multivariate stepwise regression model that simultaneously considered the effects of onset age, disease duration, sex, depression, motor symptom severity, exonic nonsynonymous SNPs of candidate genes and haplotypes of the COMT gene. Here, COMT is linked to major depressive disorder.