Yoshinao Koike, Masahiko Takahata, Masahiro Nakajima, Nao Otomo, Hiroyuki Suetsugu, Xiaoxi Liu, Tsutomu Endo, Shiro Imagama, Kazuyoshi Kobayashi, Takashi Kaito, Satoshi Kato, Yoshiharu Kawaguchi, Masahiro Kanayama, Hiroaki Sakai, Takashi Tsuji, Takeshi Miyamoto, Hiroyuki Inose, Toshitaka Yoshii, Masafumi Kashii, Hiroaki Nakashima, Kei Ando, Yuki Taniguchi, Kazuhiro Takeuchi, Shuji Ito, Kohei Tomizuka, Keiko Hikino, Yusuke Iwasaki, Yoichiro Kamatani, Shingo Maeda, Hideaki Nakajima, Kanji Mori, Atsushi Seichi, Shunsuke Fujibayashi, Tsukasa Kanchiku, Kei Watanabe, Toshihiro Tanaka, Kazunobu Kida, Sho Kobayashi, Masahito Takahashi, Kei Yamada, Hiroshi Takuwa, Hsing-Fang Lu, Shumpei Niida, Kouichi Ozaki, Yukihide Momozawa, Genetic Study Group of Investigation Committee on Ossification of the Spinal Ligaments, Masashi Yamazaki, Atsushi Okawa, Morio Matsumoto, Norimasa Iwasaki, Chikashi Terao, and Shiro Ikegawa
Ossification of the posterior longitudinal ligament of the spine (OPLL) is an intractable disease leading to severe neurological deficits. Its etiology and pathogenesis are primarily unknown. The relationship between OPLL and comorbidities, especially type 2 diabetes (T2D) and high body mass index (BMI), has been the focus of attention; however, no trait has been proven to have a causal relationship. We conducted a meta-analysis of genome-wide association studies (GWASs) using 22,016 Japanese individuals and identified 14 significant loci, 8 of which were previously unreported. We then conducted a gene-based association analysis and a transcriptome-wide Mendelian randomization approach and identified three candidate genes for each. Partitioning heritability enrichment analyses observed significant enrichment of the polygenic signals in the active enhancers of the connective/bone cell group, especially H3K27ac in chondrogenic differentiation cells, as well as the immune/hematopoietic cell group. Single-cell RNA sequencing of Achilles tendon cells from a mouse Achilles tendon ossification model confirmed the expression of genes in GWAS and post-GWAS analyses in mesenchymal and immune cells. Genetic correlations with 96 complex traits showed positive correlations with T2D and BMI and a negative correlation with cerebral aneurysm. Mendelian randomization analysis demonstrated a significant causal effect of increased BMI and high bone mineral density on OPLL. We evaluated the clinical images in detail and classified OPLL into cervical, thoracic, and the other types. GWAS subanalyses identified subtype-specific signals. A polygenic risk score for BMI demonstrated that the effect of BMI was particularly strong in thoracic OPLL. Our study provides genetic insight into the etiology and pathogenesis of OPLL and is expected to serve as a basis for future treatment development.