Panoutsopoulou Kalliope, Thiagarajah Shankar, Zengini Eleni, Day-Williams Aaron G, Ramos Yolande Fm, Meessen Jennifer Mta, Huetink Kasper, Nelissen Rob Ghh, Southam Lorraine, Rayner N William, Doherty Michael, Meulenbelt Ingrid, Zeggini Eleftheria, Wilkinson J Mark
Department of Human Genetics, Wellcome Trust Sanger Institute, Hinxton, UK.
Department of Oncology and Metabolism, University of Sheffield, Sheffield, UK.
Ann Rheum Dis. 2017 Jul;76(7):1199-1206. doi: 10.1136/annrheumdis-2016-210373. Epub 2016 Dec 14.
Osteoarthritis (OA) has a strong genetic component but the success of previous genome-wide association studies (GWAS) has been restricted due to insufficient sample sizes and phenotype heterogeneity. Our aim was to examine the effect of clinically relevant endophenotyping according to site of maximal joint space narrowing (maxJSN) and bone remodelling response on GWAS signal detection in hip OA.
A stratified GWAS meta-analysis was conducted in 2118 radiographically defined hip OA cases and 6500 population-based controls. Signals were followed up by analysing differential expression of proximal genes for bone remodelling endophenotypes in 33 pairs of macroscopically intact and OA-affected cartilage.
We report suggestive evidence (p<5×10) of association at 6 variants with OA endophenotypes that would have been missed by using presence of hip OA as the disease end point. For example, in the analysis of hip OA cases with superior maxJSN versus cases with non-superior maxJSN we detected association with a variant in the gene (rs754106, p=1.49×10, OR (95% CIs) 0.70 (0.61 to 0.80)). In the comparison of hypertrophic with non-hypertrophic OA the most significant variant was located between and (rs6766414, p=3.13×10, OR (95% CIs) 1.45 (1.24 to 1.69)). Both of these associations were fully attenuated in non-stratified analyses of all hip OA cases versus population controls (p>0.05). was significantly upregulated in OA-affected versus intact cartilage, particularly in the analysis of hypertrophic and normotrophic compared with atrophic bone remodelling pattern (p=4.2×10).
Our findings demonstrate that stratification of OA cases into more homogeneous endophenotypes can identify genes of potential functional importance otherwise obscured by disease heterogeneity.
骨关节炎(OA)具有很强的遗传成分,但由于样本量不足和表型异质性,以往全基因组关联研究(GWAS)的成效有限。我们的目的是根据最大关节间隙狭窄(maxJSN)部位和骨重塑反应进行临床相关内表型分析,以研究其对髋骨关节炎GWAS信号检测的影响。
对2118例经影像学定义的髋骨关节炎病例和6500例基于人群的对照进行分层GWAS荟萃分析。通过分析33对宏观完整和受OA影响的软骨中骨重塑内表型近端基因的差异表达,对信号进行追踪。
我们报告了6个变体与OA内表型关联的提示性证据(p<5×10),而以髋骨关节炎的存在作为疾病终点则会遗漏这些证据。例如,在分析上颌骨maxJSN较高的髋骨关节炎病例与上颌骨maxJSN不高的病例时,我们检测到与基因中的一个变体相关(rs754106,p=1.49×10,OR(95%CI)0.70(0.61至0.80))。在肥厚性OA与非肥厚性OA的比较中,最显著的变体位于和之间(rs6766414,p=3.13×10,OR(95%CI)1.45(1.24至1.69))。在所有髋骨关节炎病例与人群对照的非分层分析中,这两种关联均完全减弱(p>0.05)。与完整软骨相比,在受OA影响的软骨中显著上调,特别是在肥厚性和正常营养性与萎缩性骨重塑模式的分析中(p=4.2×10)。
我们的研究结果表明,将OA病例分层为更同质的内表型可以识别出潜在功能重要性的基因,否则这些基因会被疾病异质性所掩盖。