Genome-wide pathway-based association study implicates complement system in the development of Kashin-Beck disease in Han Chinese

Kashin-Beck disease (KBD) is a chronic osteochondropathy. The pathogenesis of KBD remains unknown. To identify relevant biological pathways for KBD, we conducted a genome-wide pathway-based association study (GWPAS) following by replication analysis, totally using 2743 Chinese Han adults. A modified gene set enrichment algorithm was used to detect association between KBD and 963 biological pathways. Cartilage gene expression analysis and serum complement measurement were performed to evaluate the functional relevance of identified pathway with KBD. We found that the Complement and Coagulation Cascades (CACC) pathway was significantly associated with KBD (P value=3.09×10^-5, false-discovery rate=0.042). Within the CACC pathway, the most significant association was observed at rs1656966 (P value=1.97×10^-4) of KNG1 gene. Further replication study observed that rs1656966 (P value=0.037) was significantly associated with KBD in an independent validation sample of 1026 subjects. Gene expression analysis observed that CFD (ratio=3.39±2.68), A2M (ratio=3.67±5.63), C5 (ratio=2.65±2.52) and CD46 (ratio=2.29±137) genes of the CACC pathway were up-regulated in KBD articular cartilage compared to healthy articular cartilage. The serum level of complement C5 in KBD patients were significantly higher than that in healthy controls (P value=0.038). Our study is the first to suggest that complement system-related CACC pathway contributed to the development of KBD.