Identification and characterization of two consistent osteoarthritis subtypes by transcriptome and clinical data integration

OBJECTIVE

To identify OA subtypes based on cartilage transcriptomic data in cartilage tissue and characterize their underlying pathophysiological processes and/or clinically relevant characteristics.

METHODS

This study includes n = 66 primary OA patients (41 knees and 25 hips), who underwent a joint replacement surgery, from which macroscopically unaffected (preserved, n = 56) and lesioned (n = 45) OA articular cartilage were collected [Research Arthritis and Articular Cartilage (RAAK) study]. Unsupervised hierarchical clustering analysis on preserved cartilage transcriptome followed by clinical data integration was performed. Protein-protein interaction (PPI) followed by pathway enrichment analysis were done for genes significant differentially expressed between subgroups with interactions in the PPI network.

RESULTS

Analysis of preserved samples (n = 56) resulted in two OA subtypes with n = 41 (cluster A) and n = 15 (cluster B) patients. The transcriptomic profile of cluster B cartilage, relative to cluster A (DE-AB genes) showed among others a pronounced upregulation of multiple genes involved in chemokine pathways. Nevertheless, upon investigating the OA pathophysiology in cluster B patients as reflected by differentially expressed genes between preserved and lesioned OA cartilage (DE-OA-B genes), the chemokine genes were significantly downregulated with OA pathophysiology. Upon integrating radiographic OA data, we showed that the OA phenotype among cluster B patients, relative to cluster A, may be characterized by higher joint space narrowing (JSN) scores and low osteophyte (OP) scores.

CONCLUSION

Based on whole-transcriptome profiling, we identified two robust OA subtypes characterized by unique OA, pathophysiological processes in cartilage as well as a clinical phenotype. We advocate that further characterization, confirmation and clinical data integration is a prerequisite to allow for development of treatments towards personalized care with concurrently more effective treatment response.

Radiographic endophenotyping in hip osteoarthritis improves the precision of genetic association analysis

OBJECTIVE

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.

METHODS

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.

RESULTS

We report suggestive evidence (p<5×10-6) 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 LRCH1 gene (rs754106, p=1.49×10-7, 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 STT3B and GADL1 (rs6766414, p=3.13×10-6, 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). STT3B 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-4).

CONCLUSIONS

Our findings demonstrate that stratification of OA cases into more homogeneous endophenotypes can identify genes of potential functional importance otherwise obscured by disease heterogeneity.

Validity of joint space width measurements in hand osteoarthritis

OBJECTIVE

To investigate the validity of joint space width (JSW) measurements in millimeters (mm) in hand osteoarthritis (OA) patients by comparison to controls, grading of joint space narrowing (JSN), and clinical features.

METHODS

Hand radiographs of 235 hand OA patients (mean age 65 years, 83% women) and 471 controls were used. JSW was measured with semi-automated image analysis software in the distal, proximal interphalangeal and metacarpal joints (DIPJs, PIPJs and MCPJs). JSN (grade 0-3) was assessed using the osteoarthritis research society international (OARSI) atlas. Associations between the two methods and clinical determinants (presence of pain, nodes and/or erosions, decreased mobility) were assessed using Generalized Estimating Equations with adjustments for age, sex, body mass index (BMI) and mean width of proximal phalanx.

RESULTS

JSW was measured in 5631 joints with a mean JSW of 0.98 mm (standard deviation (SD) 0.21), being the smallest for DIPJs (0.70 (SD 0.25)) and largest for MCPJs (1.40 (SD 0.25)). The JSN=0 group had a mean JSW of 1.28 mm (SD 0.34), the JSN=3 group 0.17 mm (SD 0.23). Controls had larger JSW than hand OA patients (P-value<0.001). In hand OA, females had smaller JSW than men (β -0.08, (95% confidence interval (95% CI) -0.15 to -0.01)) and lower JSW was associated with the presence of pain, nodes, erosions and decreased mobility (adjusted β -0.21 (95% CI -0.27, -0.16), -0.37 (-0.40, -0.34), -0.61 (-0.68, -0.54) and -0.46 (-0.68, -0.24) respectively). These associations were similar for JSN in grades.

CONCLUSION

In hand OA the quantitative JSW measurement is a valid method to measure joint space and shows a good relation with clinical features.