Genetic polymorphisms inside and outside the MHC improve prediction of AS radiographic severity in addition to clinical variables

How to translate genetic findings into clinical applications in spondyloarthritis?

Spondyloarthritis (SpA) is characterized by a strong genetic predisposition evidenced by the identification of up to 50 susceptibility loci, in addition to HLA-B27, the major genetic factor associated with the disease. These loci have not only deepened our understanding of disease pathogenesis but also offer the potential to improve disease management. Diagnostic delay is a major issue in SpA. HLA-B27 testing is widely used as diagnostic biomarker in SpA but its predictive value is limited. Several attempts have been made to develop more sophisticated polygenic risk score (PRS). However, these scores currently offer very little improvement as compared to HLA-B27 and are still difficult to implement in clinical routine. Genetics might also help to predict disease outcome including treatment response. Several genetic variants have been reported to be associated with radiographic damage or with poor response to TNF blockers, unfortunately with lack of coherence across studies. Large-scale studies should be conducted to obtain more robust findings. Genetic and genomic evidence in complex diseases can be further used to support the identification of new drug targets and to repurpose existing drugs. Although not fully driven by genetics, development of IL-17 blockers has been facilitated by the discovery of the association between IL23R variants and SpA. Development of recent approaches combining GWAS findings with functional genomics will help to prioritize new drug targets in the future. Although very promising, translational genetics in SpA remains challenging and will require a multidisciplinary approach that integrates genetics, genomics, immunology, and clinical research.

Clinical and genetic factors associated with radiographic damage in patients with ankylosing spondylitis

OBJECTIVES

To identify clinical and genetic factors associated with severe radiographic damage in patients with ankylosing spondylitis (AS).

METHODS

We newly generated genome-wide single nucleotide polymorphism data (833K) for 444 patients with AS. The severity of radiographic damage was assessed using the modified Stoke Ankylosing Spondylitis Spinal Score (mSASSS). To identify clinical and genetic factors associated with severe radiographic damage, multiple linear regression analyses were performed. Human AS-osteoprogenitor and control-osteoprogenitor cells were used for functional validation.

RESULTS

The significant clinical factors of final mSASSS were baseline mSASSS (β=0.796, p=3.22×10-75), peripheral joint arthritis (β=-0.246, p=6.85×10-6), uveitis (β=0.157, p=1.95×10-3), and smoking (β=0.130, p=2.72×10-2) after adjusting for sex, age and disease duration. After adjusting significant clinical factors, the Ryanodine receptor 3 (RYR3) gene was associated with severe radiographic damage (p=1.00×10-6). For pathway analysis, the PI3K-Akt signalling pathway was associated with severe radiographic damage in AS (p=2.21×10-4, false discovery rate=0.040). Treatment with rhodamine B, a ligand of RYR3, dose-dependently induced matrix mineralisation of AS osteoprogenitors. However, the rhodamine B-induced accelerated matrix mineralisation was not definitive in control osteoprogenitors. Knockdown of RYR3 inhibited matrix mineralisation in SaOS2 cell lines.

CONCLUSIONS

This study identified clinical and genetic factors that contributed to better understanding of the pathogenesis and biology associated with radiographic damage in AS.

Bone Involvement in Patients with Spondyloarthropathies

Spondyloarthropathies (SpA) are common systemic inflammatory rheumatic diseases, in which, as in other rheumatic diseases, levels of markers of bone resorption are elevated, leading to bone loss and elevated risk of vertebral fractures. However, the diseases are also associated with new bone formation in the spine, the so-called syndesmophytes. We tried to unravel the pathogenesis of formation and growth of syndesmophytes and evaluated new diagnostic and treatment options. After a successful meeting of the Working Group on Rheumatic Diseases at the ECTS 2020, we (WL and CR) were excited about the quality of the speakers (CM, JH, AG, and GL) and their complimentary lectures. Given the relative lack of reviews on spondyloarthropathies and bone, we decided to work together on a comprehensive review that might be interesting for basic scientists and clinically relevant for clinicians. Radiographic progression in axSpA is linked to several risk factors, like male sex, smoking, HLA-B-27, increased levels of CRP, presence of syndesmophytes, and marked inflammation on MRI. The potential role of mechanical stress in the context of physically demanding jobs has been also suggested to promote structural damages. Different treatment options from NSAIDs to biologic agents like TNF inhibitors (TNFi) or IL-17inhibitors (IL-17i) result in a reduction of inflammation and symptoms. However, all these different treatment options failed to show clear and reproducible results on inhibition on syndesmophyte formation. The majority of data are available on TNFi, and some studies suggested an effect in subgroups of patients with ankylosing spondylitis. Less information is available on NSAIDs and IL-17i. Since IL-17i have been introduced quite recently, more studies are expected. IL-17 inhibitors (Il-17i) potently reduce signs and symptoms, but serum level of IL-17 is not elevated, therefore, IL-17 probably has mainly a local effect. The failure of anti-IL-23 in axSpA suggests that IL-17A production could be independent from IL-23. It may be upregulated by TNFα, resulting in lower expression of DKK1 and RANKL and an increase in osteogenesis. In active AS markers of bone resorption are increased, while bone formation markers can be increased or decreased. Bone Turnover markers and additional markers related to Wnt such as DKK1, sclerostin, and RANKL are valuable for elucidating bone metabolism on a group level and they are not (yet) able to predict individual patient outcomes. The gold standard for detection of structural lesions in clinical practice is the use of conventional radiographics. However, the resolution is low compared to the change over time and the interval for detecting changes are 2 years or more. Modern techniques offer substantial advantages such as the early detection of bone marrow edema with MRI, the fivefold increased detection rate of new or growing syndesmophytes with low-dose CT, and the decrease in 18F-fluoride uptake during treatment with TNFα-inhibitors (TNFi) in a pilot study in 12 AS patients. Detection of bone involvement by new techniques, such as low-dose CT, MRI and 18-Fluoride PET-scans, and bone turnover markers, in combination with focusing on high-risk groups such as patients with early disease, elevated CRP, syndesmophytes at baseline, male patients and patients with HLA-B27 + are promising options for the near future. However, for optimal prevention of formation of syndesmophytes we need more detailed insight in the pathogenesis of bone formation in axSpA and probably more targeted therapies.

Prediction of Ankylosing Spondylitis in the HUNT Study by a Genetic Risk Score Combining 110 Single-nucleotide Polymorphisms of Genome-wide Significance

OBJECTIVE

The genetic component of ankylosing spondylitis (AS) development is ∼90%. Of the known heritability, ∼20% is explained by HLA-B27, and 113 identified AS-associated single-nucleotide polymorphisms (SNP) account for ∼7.4%. The objectives were to construct a weighted genetic risk score (wGRS) using currently known genome-wide susceptibility SNP, and to evaluate its predictive ability for AS in the Norwegian population-based Nord-Trøndelag Health Study (HUNT).

METHODS

AS cases (n = 164) and controls (n = 49,032) were from the second (1995-1997) and third (2006-2008) waves of the HUNT study, to which the entire adult population of the northern region of Trøndelag was invited. A wGRS based on 110 SNP weighted by published OR for AS was constructed, representing each person's carriage of all risk variants. Logistic regression models including the wGRS alone or in combination with HLA-B27 carrier state and other adjustment variables (sex, age, smoking, body mass index, and hypertension) were developed. Discrimination among models was compared using area under the curve (AUC).

RESULTS

The wGRS was associated with AS (OR 1.7, 95% CI 1.4-2.1), but showed low discrimination (AUC 0.62, 95% CI 0.58-0.67). HLA-B27 was significantly associated with AS (OR 50, 95% CI 32-81), showing high discrimination (AUC 0.88, 95% CI 0.85-0.90). Combining the wGRS and HLA-B27 improved prediction (AUC 0.90, 95% CI 0.87-0.92; p < 0.001 vs wGRS alone, p < 0.01 vs HLA-B27 alone). Further inclusion of adjustment variables gave a small improvement (AUC 0.91, 95% CI 0.89-0.94; p = 0.03).

CONCLUSION

Prediction in a population-based setting based on all currently known AS susceptibility SNP was better than HLA-B27 carrier state alone, although the improvement was small and of uncertain clinical value.

Neural EGFL like 1 as a potential pro-chondrogenic, anti-inflammatory dual-functional disease-modifying osteoarthritis drug

Arthritis, an inflammatory condition that causes pain and cartilage destruction in joints, affects over 54.4 million people in the US alone. Here, for the first time, we demonstrated the emerging role of neural EGFL like 1 (NELL-1) in arthritis pathogenesis by showing that Nell-1-haploinsufficient (Nell-1^+/6R) mice had accelerated and aggravated osteoarthritis (OA) progression with elevated inflammatory markers in both spontaneous primary OA and chemical-induced secondary OA models. In the chemical-induced OA model, intra-articular injection of interleukin (IL)1β induced more severe inflammation and cartilage degradation in the knee joints of Nell-1^+/6R mice than in wildtype animals. Mechanistically, in addition to its pro-chondrogenic potency, NELL-1 also effectively suppressed the expression of inflammatory cytokines and their downstream cartilage catabolic enzymes by upregulating runt-related transcription factor (RUNX)1 in mouse and human articular cartilage chondrocytes. Notably, NELL-1 significantly reduced IL1β-stimulated inflammation and damage to articular cartilage in vivo. In particular, NELL-1 administration markedly reduced the symptoms of antalgic gait observed in IL1β-challenged Nell-1^+/6R mice. Therefore, NELL-1 is a promising pro-chondrogenic, anti-inflammatory dual-functional disease-modifying osteoarthritis drug (DMOAD) candidate for preventing and suppressing arthritis-related cartilage damage.

A Single Nucleotide Polymorphism in the Il17ra Promoter Is Associated with Functional Severity of Ankylosing Spondylitis

The aim of this study was to identify new genetic variants associated with the severity of ankylosing spondylitis (AS). We sequenced the exome of eight patients diagnosed with AS, selected on the basis of the severity of their clinical parameters. We identified 27 variants in exons and regulatory regions. The contribution of candidate variants found to AS severity was validated by genotyping two Spanish cohorts consisting of 180 cases/300 controls and 419 cases/656 controls. Relationships of SNPs and clinical variables with the Bath Ankylosing Spondylitis Disease Activity and Functional Indices BASDAI and BASFI were analyzed. BASFI was standardized by adjusting for the duration of the disease since the appearance of the first symptoms. Refining the analysis of SNPs in the two cohorts, we found that the rs4819554 minor allele G in the promoter of the IL17RA gene was associated with AS (p<0.005). This variant was also associated with the BASFI score. Classifying AS patients by the severity of their functional status with respect to BASFI/disease duration of the 60th, 65th, 70th and 75th percentiles, we found the association increased from p60 to p75 (cohort 1: p<0.05 to p<0.01; cohort 2: p<0.01 to p<0.005). Our findings indicate a genetic role for the IL17/ILRA axis in the development of severe forms of AS.

Improved prediction of knee osteoarthritis progression by genetic polymorphisms: the Arthrotest Study

OBJECTIVE

The aim of this study was to develop a genetic prognostic tool to predict radiographic progression towards severe disease in primary knee OA (KOA) patients.

METHODS

This investigation was a cross-sectional, retrospective, multicentric association study in 595 Spanish KOA patients. Caucasian patients aged ≥40 years at the time of diagnosis of primary KOA of Kellgren-Lawrence grade 2 or 3 were included. Patients who progressed to Kellgren-Lawrence score 4 or who were referred for total knee replacement within 8 years after diagnosis were classified as progressors to severe disease. Clinical variables of the initial stages of the disease (gender, BMI, age at diagnosis, OA in the contralateral knee, and OA in other joints) were registered as potential predictors. Single nucleotide polymorphisms and clinical variables with an association of P < 0.05 were included in the multivariate analysis using forward logistic regression.

RESULTS

A total of 23 single nucleotide polymorphisms and the time of primary KOA diagnosis were significantly associated with KOA severe progression in the exploratory cohort (n = 220; P < 0.05). The predictive accuracy of the clinical variables was limited: area under the curve (AUC) = 0.66. When genetic variables were added to the clinical model (full model), the prediction of KOA progression was significantly improved (AUC = 0.82). Combining only genetic variables (rs2073508, rs10845493, rs2206593, rs10519263, rs874692, rs7342880, rs780094 and rs12009), a predictive model with good accuracy was also obtained (AUC = 0.78). The predictive ability for KOA progression of the full model was confirmed on the replication cohort (two-sample Z-test; n = 62; P = 0.190).

CONCLUSION

An accurate prognostic tool to predict primary KOA progression has been developed based on genetic and clinical information from OA patients.

Bone morphogenetic protein 6 polymorphisms are associated with radiographic progression in ankylosing spondylitis

BACKGROUND AND OBJECT

Nearly 25 genetic loci associated with susceptibility to ankylosing spondylitis (AS) have been identified by several large studies. However, there have been limited studies to identify the genes associated with radiographic severity of the disease. Thus we investigated which genes involved in bone formation pathways might be associated with radiographic severity in AS.

METHODS

A total of 417 Korean AS patients were classified into two groups based on the radiographic severity as defined by the modified Stoke' Ankylosing Spondylitis Spinal Score (mSASSS) system. Severe AS was defined by the presence of syndesmophytes and/or fusion in the lumbar or cervical spine (n = 195). Mild AS was defined by the absence of any syndesmophyte or fusion (n = 170). A total of 251 single nucleotide polymorphisms (SNPs) within 52 genes related to bone formation were selected and genotyped. Odds ratios (OR) and 95% confidence interval (95% CI) were analysed by multivariate logistic regression controlling for age at onset of symptoms, sex, disease duration, and smoking status as covariates.

RESULTS

We identified new loci of bone morphogenetic protein 6 (BMP6) associated with radiographic severity in patients with AS that passed false discovery rate threshold. Two SNPs in BMP6 were significantly associated with radiologic severity [rs270378 (OR 1.97, p = 6.74 × 10(-4)) and rs1235192 [OR 1.92, p = 1.17 × 10(-3)]) adjusted by covariates.

CONCLUSION

This is the first study to demonstrate that BMP6 is associated with radiographic severity in AS, supporting the role wingless-type like/BMP pathway on radiographic progression in AS.

Association study of genes related to bone formation and resorption and the extent of radiographic change in ankylosing spondylitis

OBJECTIVE

To identify genetic associations with severity of radiographic damage in ankylosing spondylitis (AS).

METHOD

We studied 1537 AS cases of European descent; all fulfilled the modified New York Criteria. Radiographic severity was assessed from digitised lateral radiographs of the cervical and lumbar spine using the modified Stoke Ankylosing Spondylitis Spinal Score (mSASSS). A two-phase genotyping design was used. In phase 1, 498 single nucleotide polymorphisms (SNPs) were genotyped in 688 cases; these were selected to capture >90% of the common haplotypic variation in the exons, exon-intron boundaries, and 5 kb flanking DNA in the 5' and 3' UTR of 74 genes involved in anabolic or catabolic bone pathways. In phase 2, 15 SNPs exhibiting p<0.05 were genotyped in a further cohort of 830 AS cases; results were analysed both separately and in combination with the discovery phase data. Association was tested by contingency tables after separating the samples into 'mild' and 'severe' groups, defined as the bottom and top 40% by mSASSS, adjusted for gender and disease duration.

RESULTS

Experiment-wise association was observed with the SNP rs8092336 (combined OR 0.32, p=1.2×10(-5)), which lies within RANK (receptor activator of NFκB), a gene involved in osteoclastogenesis, and in the interaction between T cells and dendritic cells. Association was also found with the SNP rs1236913 in PTGS1 (prostaglandin-endoperoxide synthase 1, cyclooxygenase 1), giving an OR of 0.53 (p=2.6×10(-3)). There was no observed association between radiographic severity and HLA-B*27.

CONCLUSIONS

These findings support roles for bone resorption and prostaglandins pathways in the osteoproliferative changes in AS.

Advances in the genetics of spondyloarthritis and clinical implications

Spondyloarthritis (SpA) is a group of inflammatory rheumatic diseases whose main clinical feature is inflammation of the axial spine. Articular, periarticular, and extra-articular manifestations can also occur, depending on the type of spondyloarthritis. The most common clinical subsets of SpA are ankylosing spondylitis (AS) and psoriatic arthritis (PsA). SpA is a major health challenge given the propensity to affect young adults and the potential requirement for lifelong treatment. Although the precise etiology of SpA is unknown, there is mounting evidence that these diseases are a result of complex interplay of genetic, environmental, and immunological factors. In this review on SpA, we will discuss genetic variants with genome-wide significance, highlight potential clinical application of genetic variants, and discuss challenges in further elucidating the genetic basis of SpA.

The genetics of ankylosing spondylitis and axial spondyloarthritis

Ankylosing spondylitis (AS) and spondyloarthritis are strongly genetically determined. The long-standing association with HLA-B27 is well described, although the mechanism by which that association induces AS remains uncertain. Recent developments include the description of HLA-B27 tag single nucleotide polymorphisms in European and Asian populations. An increasing number of non-MHC genetic associations have been reported, which provided amongst other things the first evidence of the involvement of the IL-23 pathway in AS. The association with ERAP1 is now known to be restricted to HLA-B27 positive disease. Preliminary studies on the genetics of axial spondyloarthritis demonstrate a lower HLA-B27 carriage rate compared with AS. Studies with larger samples and including non-European ethnic groups are likely to further advance the understanding of the genetics of AS and spondyloarthritis.

Both baseline clinical factors and genetic polymorphisms influence the development of severe functional status in ankylosing spondylitis

Functional severity in ankylosing spondylitis (AS) patients is variable and difficult to predict early. The aim of our study was to assess whether a combination of baseline clinical factors and genetic markers may predict the development of severe functional status in AS. We performed a cross-sectional association study on AS patients included in the Spanish National Registry of Spondyloarthropathies—REGISPONSER. Bath Ankylosing Spondylitis Functional Index (BASFI) was standardized by adjusting for disease duration since the first symptoms (BASFI/t). We considered as severe functional status the values of BASFI/t in the top of the 60th (p60), 65th (p65), 70th (p70), and 75th (p75) percentile. We selected 384 single nucleotide polymorphisms (SNPs) distributed in 190 genes to be analyzed. The study cohort included 456 patients with mean age 50.8(±10.5) years and with mean disease duration since first symptoms 24.7 (±10.1) years. Older age at disease onset and neck pain at baseline showed statistical significant association with severe BASFI/t. Polymorphisms associated in the allele frequencies test with severe BASFI/t in all classifications were: rs2542151 (p60 [P = .04], p65 [P = .04], p70 [P = .001] and p75 [P = .001]) and rs2254441 (p60 [P = .004], p65 [P = .02], p70 [P = .01] and p75 [P<.001]).. Genotype association, after adjustment for covariates, found an association in three of the four patients' classifications for rs2542151 and in two of the classifications for rs2254441.Forward logistic regression did not identify any model with a good predictive power for severe functional development.

In our study we identified clinical factors and 24 polymorphisms associated with development of severe functional status in AS patients. Validation of these results in other cohorts is required.