Methodological development of molecular endotype discovery from synovial fluid of individuals with knee osteoarthritis: the STEpUP OA Consortium
Authors: Deng Y. , Perry T.A. , Hulley P. , Maciewicz R.A. , Mitchelmore J. , Perry D. , Larsson S. , Brachat S. , Struglics A. , Appleton C.T. , Kluzek S. , Arden N. K. , Felson D. , Marsden B. , Tom B.D.M. , Bondi L. , Kapoor M. , Batchelor V. , Mackay-Alderson J. , Kumar V. , Lohmander L. S. , Welting T. J. , Walsh D. A. , Valdes A.M. ,   the STEpUP OA Consortium, Vincent T. L. , Watt F. E. , Jostins-Dean L. 
Published: August 2023

Objectives

To develop and validate a pipeline for quality controlled (QC) protein data for largescale analysis of synovial fluid (SF), using SomaLogic technology.

Design

Knee SF and associated clinical data were from partner cohorts. SF samples were centrifuged, supernatants stored at −80 °C, then analysed by SomaScan Discovery Plex V4.1 (>7000 SOMAmers/proteins).

Setting

An international consortium of 9 academic and 8 commercial partners (STEpUP OA).

Participants

1746 SF samples from 1650 individuals comprising OA, joint injury, healthy controls and inflammatory arthritis controls, divided into discovery (n=1045) and replication (n=701) datasets.

Primary and secondary outcome measures

An optimised approach to standardisation was developed iteratively, monitoring reliability and precision (comparing coefficient of variation [%CV] of ‘pooled’ SF samples between plates and correlation with prior immunoassay for 9 analytes). Pre-defined technical confounders were adjusted for (by Limma) and batch correction was by ComBat. Poorly performing SOMAmers and samples were filtered. Variance in the data was determined by principal component (PC) analysis. Data were visualised by Uniform Manifold Approximation and Projection (UMAP).

Results

Optimal SF standardisation aligned with that used for plasma, but without median normalisation. There was good reliability (<20 %CV for >80% of SOMAmers in pooled samples) and overall good correlation with immunoassay. PC1 accounted for 48% of variance and strongly correlated with individual SOMAmer signal intensities (median correlation coefficient 0.70). These could be adjusted using an ‘intracellular protein score’. PC2 (7% variance) was attributable to processing batch and was batch-corrected by ComBat. Lesser effects were attributed to other technical confounders. Data visualisation by UMAP revealed clustering of injury and OA cases in overlapping but distinguishable areas of high-dimensional proteomic space.

Conclusions

We define a standardised approach for SF analysis using the SOMAscan platform and identify likely ‘intracellular’ protein as being a major driver of variance in the data.

Strengths and limitations

This is the largest number of individual synovial fluid samples analysed by a high content proteomic platform (SomaLogic technology)

SomaScan offers reliable, precise relative SF data following standardisation for over 6000 proteins

Significant variance in the data was driven by a protein signal which is likely intracellular in origin: it is not yet clear whether this is due to technical considerations, normal cell turnover or relevant pathological processes

Adjusting for confounding factors might conceal the true structure of the data and reduce the ability to detect ‘molecular endotypes’ within disease groups