Novel Genetic Variants Associated with Lumbar Spondylosis in Koreans : A Genome-Wide Association Study

Objective : The aim of this study was to identify the susceptibility genes responsible for lumbar spondylosis (LS) in Korean patients. Methods : Data from 1427 subjects were made available for radiographic grading and genome wide association studies (GWAS) analysis. Lateral lumbar spine radiographs were obtained and the various degrees of degenerative change were semi-quantitatively scored. A pilot GWAS was performed using the AffymetrixGenome-Wide Human single-nucleotide polymorphisms (SNPs), 500K array. A total of 352228 SNPs were analyzed and the association between the SNPs and case-control status was analyzed by stepwise logistic regression analyses. Results : The top 100 SNPs with a cutoff p-value of less than $3.7{\times}10^{-4}$ were selected for joint space narrowing, while a cutoff p-value of $6.0{\times}10^{-4}$ was applied to osteophytes and the Kellgren-Lawrence (K-L) osteoarthritis grade. The SNPs with the strongest effect on disc space narrowing, osteophytes, and K-L grade were serine incorporator 1 (rs155467, odds ratio [OR]=17.58, $p=1.6{\times}10^{-4}$), stromal interaction molecule 2 (STIM1, rs210781, OR=5.53, $p=5{\times}10^{-4}$), and transient receptor potential cation channel, subfamily C (rs11224760, OR=3.99, $p=4.8{\times}10^{-4}$), respectively. Leucine-rich repeat-containing G protein-coupled receptor 4 was significantly associated with both disc space narrowing and osteophytes (rs1979400, OR=2.01, $p=1.1{\times}10^{-4}$ for disc space narrowing, OR=1.79, $p=3{\times}10^{-4}$ for osteophytes), while zinc finger and BTB domain containing 7C was significantly and negatively associated with both osteophytes and a K-L grade >2 (rs12457004,OR=0.25, $p=5.8{\times}10^{-4}$ and OR=0.27, $p=5.3{\times}10^{-4}$, respectively). Conclusion : We identified SNPs that potentially contribute to the pathogenesis of LS. This is the first report of a GWAS in an Asian population.