SummaryObjectivesJoint degeneration in osteoarthritis (OA) is characterised by damage and loss of articular cartilage. The pattern of loss is consistent with damage occurring only where the mechanical loading is high. We have investigated using RNA-sequencing (RNA-seq) and systems analyses the changes that occur in damaged OA cartilage by comparing it with intact cartilage from the same joint.MethodsCartilage was obtained from eight OA patients undergoing total knee replacement. RNA was extracted from cartilage on the damaged distal medial condyle (DMC) and the intact posterior lateral condyle (PLC). RNA-seq was performed to identify differentially expressed genes (DEGs) and systems analyses applied to identify dysregulated pathways.ResultsIn the damaged OA cartilage, there was decreased expression of chondrogenic genes SOX9, SOX6, COL11A2, COL9A1/2/3, ACAN and HAPLN1; increases in non-chondrogenic genes COL1A1, COMP and FN1; an altered pattern of secreted proteinase expression; but no expression of major inflammatory cytokines. Systems analyses by PhenomeExpress revealed significant sub-networks of DEGs including mitotic cell cycle, Wnt signalling, apoptosis and matrix organisation that were influenced by a core of altered transcription factors (TFs), FOSL1, AHR, E2F1 and FOXM1.ConclusionsGene expression changes in damaged cartilage suggested a signature non-chondrogenic response of altered matrix protein and secreted proteinase expression. There was evidence of a damage response in this late OA cartilage, which surprisingly showed features detected experimentally in the early response of cartilage to mechanical overload. PhenomeExpress analysis identified a hub of DEGs linked by a core of four differentially regulated TFs.

IntroductionOsteoarthritis (OA) is a heterogeneous and complex disease. We have used a network biology approach based on genome-wide analysis of gene expression in OA knee cartilage to seek evidence for pathogenic mechanisms that may distinguish different patient subgroups.MethodsResults from RNA-Sequencing (RNA-Seq) were collected from intact knee cartilage at total knee replacement from 44 patients with OA, from 16 additional patients with OA and 10 control patients with non-OA. Results were analysed to identify patient subsets and compare major active pathways.ResultsThe RNA-Seq results showed 2692 differentially expressed genes between OA and non-OA. Analysis by unsupervised clustering identified two distinct OA groups: Group A with 24 patients (55%) and Group B with 18 patients (41%). A 10 gene subgroup classifier was validated by RT-qPCR in 16 further patients with OA. Pathway analysis showed increased protein expression in both groups. PhenomeExpress analysis revealed group differences in complement activation, innate immune responses and altered Wnt and TGFβ signalling, but no activation of inflammatory cytokine expression. Both groups showed suppressed circadian regulators and whereas matrix changes in Group A were chondrogenic, in Group B they were non-chondrogenic with changes in mechanoreceptors, calcium signalling, ion channels and in cytoskeletal organisers. The gene expression changes predicted 478 potential biomarkers for detection in synovial fluid to distinguish patients from the two groups.ConclusionsTwo subgroups of knee OA were identified by network analysis of RNA-Seq data with evidence for the presence of two major pathogenic pathways. This has potential importance as a new basis for the stratification of patients with OA for drug trials and for the development of new targeted treatments.