{'fraction_non_alphanumeric': 0.040544913279834324, 'fraction_numerical': 0.03148459746311157, 'mean_word_length': 4.835536253776435, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 5, 'https://': 1, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The CRTAC1 Protein in Plasma Is Associated With Osteoarthritis and Predicts Progression to Joint Replacement: A Large‐Scale Proteomics Scan in Iceland

Objective Biomarkers for diagnosis and progression of osteoarthritis (OA) are lacking. This study was undertaken to identify circulating biomarkers for OA that could predict disease occurrence and/or progression to joint replacement. Methods Using the SomaScan platform, we measured 4,792 proteins in plasma from 37,278 individuals, of whom 12,178 individuals had OA and 2,524 had undergone joint replacement. We performed a case–control study for identification of potential protein biomarkers for hip, knee, and/or hand OA, and a prospective study for identification of biomarkers for joint replacement. Results Among the large panel of plasma proteins assessed, cartilage acidic protein 1 (CRTAC1) was the most strongly associated with both OA diagnosis (odds ratio 1.46 [95% confidence interval 1.41–1.52] for knee OA, odds ratio 1.36 [95% confidence interval 1.29–1.43] for hip OA, and odds ratio 1.33 [95% confidence interval 1.26–1.40] for hand OA) and progression to joint replacement (hazard ratio 1.40 [95% confidence interval 1.30–1.51] for knee replacement and hazard ratio 1.31 [95% confidence interval 1.19–1.45] for hip replacement). Patients with OA who were in the highest quintile of risk of joint replacement, based on known risk factors (i.e., age, sex, and body mass index) and plasma CRTAC1 level, were 16 times more likely to undergo knee replacement within 5 years of plasma sample collection than those in the lowest quintile, and 6.5 times more likely to undergo hip replacement. CRTAC1 was not associated with other types of inflammatory arthritis. A specific protein profile was identified in those patients who had undergone joint replacement prior to plasma sample collection. Conclusion Through a hypothesis‐free approach, we identified CRTAC1 in plasma as a novel promising candidate biomarker for OA that is both associated with occurrence of OA and predictive of progression to joint replacement. This biomarker might also be useful in the selection of suitable patients for clinical trial enrollment.


INTRODUCTION
Osteoarthritis (OA) is a major global health burden affecting >300 million people worldwide (1). Prevalence of OA is increasing since the prevalence of its main risk factors, obesity and older age, continue to rise (2)(3)(4). The disease can have a great impact on quality of life due to pain and loss of joint function, mostly affecting the knees, hips, and hands. OA is a heterogeneous disease with variable radiographic and clinical features (5), and there are presently no measures available for early diagnosis before destructive Supported by deCODE genetics, Inc./Amgen Inc. A portion of the study was conducted using the UK Biobank Resource under application number 23359.
changes are observable on radiographs. The current therapeutic approach consists of pain medications, lifestyle changes with weight reduction and exercise, and joint arthroplasty for severe disease (3), since no disease-modifying drugs are available. Thus, there is an unmet need for treatment that could effectively slow or halt OA progression.
A biomarker that is associated with OA disease occurrence and/or progression would help to identify cases earlier or monitor the disease course. Lack of such biomarkers has hindered development of effective therapy for this common disease. There are several studies on candidate biomarkers for OA (6)(7)(8), with somewhat inconclusive findings, possibly reflecting heterogeneity in the disease definition, study design, or size. However, meta-analyses of the 2 most extensively studied biomarkers for OA, serum cartilage oligomeric matrix protein (COMP) and urinary C-terminal telopeptide of type II collagen (CTX-II), show correlation with both disease occurrence and progression (9)(10)(11)(12)(13). These biomarkers are, however, not currently used in clinical settings.
In this study, we conducted a large-scale screen of biomarkers for OA, examining 4,792 human proteins in plasma from 39,155 individuals. We aimed to find a biomarker that is specific for OA development and/or disease progression, preferably a single biomarker since this is most feasible in the clinical setting, but also a more comprehensive tool that could yield better prediction. We included plasma proteins and known risk factors for OA (i.e., age, sex, and body mass index [BMI]) in the association and prediction models. To improve the models, we also included polygenic risk scores for OA, accounting for the role of genetics, as >80 independent sequence variants have been associated with the risk of OA to date (14,15).

SUBJECTS AND METHODS
Study population. We generated proteomics data from plasma collected from 39,155 Icelanders between August 2000 and January 2019, using 4,963 slow off-rate modified aptamers (SOMAmers) (SomaLogic; https://somal ogic.com/) that measured 4,792 individual human plasma proteins (Supplementary Figure 1 We designed 3 main studies from these 39,155 samples: 1) case-control association studies of OA, including 6,136 subjects with primary OA of the hips, knees, or hands (excluding those with prior joint replacement) and 13,789-16,001 controls without OA, 2) a prospective study of progression to joint replacement, including 10,701 subjects with OA, and 3) a case-control association study of prior joint replacement, including 1,452 cases and 35,826 controls. Subjects with missing information regarding age, sex, or BMI were excluded from all analyses, as were those who had rheumatoid arthritis (RA) or were age <40 years at the time of plasma sampling from analyses 1 and 2 ( Figure 1 and Supplementary Methods).
In a secondary analysis of hand OA severity, we used 2 data sets. One included patients with hand OA for whom the severity of disease was evaluated by radiography and clinical findings. The second data set included individuals who participated in the population-based deCODE Health study (16) and had available high-quality hand photographs for analysis of hand OA (17,18). Subjects who participated in the deCODE Health study also answered a simple question regarding pain. Details of these analyses are available in Supplementary Methods.
The study was approved by the National Bioethics Committee of Iceland (VSN_14-015v8, VSN_14-148, and VSN_15-214), and was conducted in accordance with requirements issued by the Data Protection Authority of Iceland. All participants were of Icelandic descent.

Comparison of SOMAmer measurements with Olink assay.
To test the specificity of the SomaScan measurements, we measured proteins using a proximity extension assay method (Olink Bioscience) (19). Proteins were measured using a cardiometabolic panel of 92 unique proteins, including cartilage acidic protein 1 (CRTAC1) and COMP. Of 200 samples, all of which had been assayed by SomaScan, 199 passed quality control. The median correlation between proteins was 0.757, ranging from −0.128 to 0.949, with an interquartile range of 0.541 to 0.829.

Polygenic risk scores.
A polygenic risk score integrates a large fraction of the genetics contributing to a disease. We generated the polygenic risk scores based on effect estimates from OA genome-wide association study (GWAS) data in the UK Biobank (www.ukbio bank.ac.uk) (14) and calculated the risk scores for genotyped Icelanders, using 600,000 common variants and LDpred software, essentially as previously described (20). Details are provided in the Supplementary Methods (http://onlin e libr ary.wiley.com/doi/10.1002/art.41793/ abstract).

Genome-wide association of plasma proteins (protein quantitative trait locus [QTL] study).
To identify a set of independent DNA sequence variants that are associated with each protein, we performed a genome-wide association analysis (significance threshold P < 5 × 10 −8 ) with a subsequent recursive Statistical analysis. Plasma protein levels were adjusted for the age of the individual at the time of plasma collection, sex, collection site, and the storage age of the sample. After adjustment, the plasma protein levels were rank-transformed onto the standard normal distribution with a mean of 0 and SD of 1. All modeling was completed with standardized protein levels, whereas raw unadjusted levels of CRTAC1 were used in visual representations (e.g., correlation with hand OA severity in Figure 2 and Kaplan-Meier curves in Figure 3).
Differences in plasma protein levels after joint replacement as well as their association with hip, knee, and hand OA were estimated with logistic regression, and area under the curve (AUC) was estimated using the bootstrap method. The risk of hip and knee replacement from the time of sample collection was estimated with Cox proportional hazards regression and visualized with Kaplan-Meier curves. Association of plasma protein levels with hand OA severity was estimated by linear regression. We used Bonferroni correction to adjust for multiple testing, yielding the following P values for each model: 1) P ≤ 3.3 × 10 −6 in analyses of association with OA, accounting for 3 phenotypes and 4,983 aptamers (which capture 4,792 proteins); 2) P ≤ 1.0 × 10 −5 in analyses of association with joint replacement, accounting for 1 phenotype and 4,983 aptamers; and 3) P ≤ 5.0 × 10 −6 in analyses of risk of progression to joint replacement, accounting for 2 phenotypes and 4,983 aptamers. A detailed description of the statistical methods is provided in the Supplementary Methods (http://onlin elibr ary.wiley.com/doi/10.1002/art.41793/ abstract).
Data availability. Data supporting the findings of this study are presented herein and in the supplementary data files, and also can be provided upon request from the corresponding author.

RESULTS
Association of plasma proteins with OA diagnosis. We tested whether any of the measured plasma proteins were associated with OA, adjusting the analysis for known disease risk factors (i.e., older age, female sex, and higher BMI). This adjustment is critical to exclude confounders, as a vast number of the 4,792 measured proteins are associated with these traits (e.g., 96% of the proteins tested were significantly associated with BMI in our data set).
We identified 45 proteins associated with knee OA, 7 proteins associated with hip OA, and 44 proteins associated with hand OA, that met the threshold for significance when accounting for multiple testing (P ≤ 3.3 × 10 −6 ; accounting for 3 phenotypes and 4,983 protein aptamers) and after excluding proteins that were also associated with inflammatory arthritis (Supplementary Table 1, available on the Arthritis & Rheumatology website at http://onlin elibr ary.wiley. Figure 2. Correlation between cartilage acidic protein 1 (CRTAC1) levels in plasma and hand osteoarthritis severity scores. Correlation between unstandardized levels of CRTAC1 in plasma and hand osteoarthritis severity scores was estimated based on digital photographs (n = 5,445). Photographs were taken at the same time plasma samples were collected from participants in the deCODE Health study (18) and scored for the presence and severity of osteoarthritis, resulting in an aggregate score. Patients with rheumatoid arthritis were excluded. Data are presented as box plots, where the boxes represent the 25th to 75th percentiles, the lines within the boxes represent the median, and the lines outside the boxes represent the 10th to 90th percentiles. Circles represent outliers. All patients were age >40 years. RFU = relative fluorescence units.   (Table 1 and Supplementary Tables 1 and 2, available on the Arthritis & Rheumatology website at http://onlin elibr ary.wiley.com/doi/10.1002/art.41793/ abstract). The plasma CRTAC1 level was more strongly associated with OA by several orders of magnitude compared to any of the other proteins measured. COMP, an investigative biomarker for OA, was the second most significantly associated protein (OR 1.31 for knee OA [P = 8.5 × 10 −43 ], 1.26 for hip OA [P = 4.5 × 10 −20 ], and 1.19 for hand OA [P = 2.3 × 10 −10 ]). COMP was somewhat correlated with CRTAC1 in our data (r 2 = 0.362 [95% confidence interval 0.355-0.370]). In models conditioning on the CRTAC1 levels, the association of COMP with OA was no longer significant (P = 0.013, P = 0.029, and P = 0.97 for knee, hip, and hand OA, respectively), whereas the association of CRTAC1 with OA remained significant in models conditioning on the COMP levels (P = 5.3 × 10 −47 , P = 5.3 × 10 −18 , and P = 3.0 × 10 −18 for knee, hip, and hand OA, respectively).
The majority of the associated proteins correlated nominally with all 3 forms of OA. For 5 of these proteins, CRTAC1, COMP, cartilage intermediate-layer protein, retinoblastoma-like 2, and cytokine-like 1, the correlations with all 3 forms of OA were significant (Supplementary Table 1).
Approximately half of our samples were derived from a cancer study (The Icelandic Cancer Project, see Figure 1 and Supplementary Methods), which may have affected the results of our association analysis. However, when we tested the association of the proteins with any type of cancer, CRTAC1 levels ranked as 1,337. Furthermore, adjusting the association analysis for cancer did not significantly change the results for CRTAC1 (OR 1.46 for knee OA, 1.36 for hip OA, and 1.32 for hand OA) compared to not adjusting for cancer (OR 1.46 for knee OA, 1.33 for hip OA, and 1.26 for hand OA). Likewise, although other comorbidities were more common in OA cases compared to controls (Table 2), these conditions were associated with lower levels of CRTAC1, whereas OA was associated with increased levels of CRTAC1 (Supplementary Table 3, http://onlin elibr ary.wiley.com/doi/10.1002/art.41793/ abstract). This highlights that CRTAC1 is a biomarker for OA and this association is not influenced by underlying comorbidities.
Correlation of plasma CRTAC1 levels with joint pain and OA disease severity. We further analyzed whether plasma levels of CRTAC1, the most significantly associated protein in OA, correlated with joint pain (a characteristic of OA) or with OA disease severity in the deCODE Health study, where the OA-related information was gathered at the time of plasma collection. Information on joint pain was gathered through a simple yes/no question on pain (answered "yes" by 649 of 2,708 patients), and hand OA severity was scored based on high-quality digital photographs (17,18)  CRTAC1 levels were associated with both joint pain (OR 1.14, P = 0.0046) and the hand OA severity score (OR 1.08, P = 8.7 × 10 −11 ), with a direction of effects consistent with that of the association with OA overall. Figure 2 shows an increase in raw, unstandardized CRTAC1 levels (not adjusted for age or sex) with the hand OA severity score. With each increase in score, CRTAC1 levels were increased by 5.7% (P = 2.2 × 10 −48 ).
We also assessed whether CRTAC1 levels correlated with thumb and finger OA severity scores, defined by radiography and clinical assessment in OA patients (n = 2,347) (Supplementary Table 5, http://onlin elibr ary.wiley.com/doi/10.1002/art.41793/ abstract). We observed that CRTAC1 levels were associated with severity scores in the finger joints (OR 1.12, P = 0.0030), specifically the thumb carpometacarpal joint (OR 1.09, P = 0.022) of these patients (Table 1). We noted that CRTAC1 was the only protein associated with both OA joint pain and severity of hand OA (Supplementary Table 1).
We evaluated whether CRTAC1 levels increased with the number of OA-affected joint types. In these analyses, we observed that the CRTAC1 plasma levels increased by 0. Plasma CRTAC1 levels as a predictor of joint replacement in patients with OA. Joint replacement is a treatment for severe OA symptoms. We used Cox proportional hazards regression to investigate whether any of the 4,792 proteins measured in plasma could predict the time until joint replacement among patients with OA (n = 10,701), limiting the study group to those whose plasma samples were collected before the joint replacement (described below). Since in many cases the specific diagnosis of knee OA or hip OA was not recorded in the electronic health record registries used in this study until the time of joint replacement surgery, we included any patients diagnosed as having OA in this analysis and not only those who were coded as M16.0, M16.1, M17.0, and M17.1 as was the case in the knee and hip association analysis (Figure 1). CRTAC1 was 1 of only 2 proteins that could significantly predict hip replacement, whereas 8 proteins could predict knee replacement (Supplementary Table 6 . Therefore, of all the proteins tested, CRTAC1 was most significantly associated with OA diagnosis and was the best predictor of progression to joint replacement. We assessed the likelihood of joint replacement over time using plasma CRTAC1 levels and classic risk factors for OA (i.e., age, sex, and BMI), stratified into quintiles of risk scores, and compared the risk of joint replacement for each quintile of the predicted risk using the Kaplan-Meier estimator. Patients with OA who were in the highest quintile of risk were 16 times more likely to undergo knee replacement within 5 years of plasma sample collection ‡ than those in the lowest quintile, and 6.5 times more likely to undergo hip replacement (Figure 3). Using unadjusted CRTAC1 levels in the absence of classic risk factors, those in the highest quintile of risk were 3 times more likely to undergo knee replacement and 2.5 times more likely to undergo hip replacement within 5 years of plasma sample collection than those who were in the lowest quintile, which demonstrated that independent of age, sex, or BMI, the CRTAC1 levels alone substantially increased the likelihood of joint replacement (Supplementary Figure 3  * Except where indicated otherwise, values are the number (%). All subjects in both studies were age >40 years at the time of plasma collection, had information on body mass index (BMI) within 5 years of plasma collection, had not undergone joint replacement, were not diagnosed as having rheumatoid arthritis (RA), and had protein levels measured using SOMAScan. NA = not applicable; IQR = interquartile range. † In the case-control association study, patients were required to have a primary diagnosis of hand, knee, or hip osteoarthritis (OA) for inclusion. ‡ In the progression to joint replacement study, a diagnosis of any OA among the individuals undergoing surgery was sufficient for inclusion. § Missing genotype is only an exclusion criterion in the full models. ‡ Correction added after online publication 18 October 2021: the text was changed from "within 5 years of diagnosis" to "within 5 years of plasma sample collection."

| 2031
Generalizability of CRTAC1 as a biomarker. To address the stability of CRTAC1 levels as a biomarker for OA, we assessed its association with age, sex, and BMI (adjusted for time from sample collection, and source of sample collection). CRTAC1 levels did not correlate with BMI (P = 0.71), but they did increase slightly with age (8.24 relative fluorescence units per year, P < 1 × 10 −300 ) (Figure 4). In an association analysis stratified by sex, age, and BMI, we showed that the effect of CRTAC1 on OA risk did not differentiate significantly between any of the groups tested (Supplementary Table 7, http://onlin elibr ary.wiley.com/ doi/10.1002/art.41793/ abstract), demonstrating that CRTAC1 was an independent risk factor for OA.

Potential role of CRTAC1 in the pathogenesis of OA.
To determine whether CRTAC1 levels in plasma contributed to the pathogenesis of OA, we performed a GWAS of DNA sequence variants associated with CRTAC1 levels (protein QTL study) and identified 8 CRTAC1 protein QTL variants, one of which was at the CRTAC1 locus. We then tested those protein QTL variants for association with OA in a combined study population from Iceland and the UK Biobank (14), with a total of 17,151 patients with hip OA, 23,877 with knee OA, and 9,773 with hand OA, compared to >560,000 controls. None of these sequence variants were associated with OA (Supplementary Table 8, http://onlin e libr ary.wiley.com/doi/10.1002/art.41793/ abstract). Furthermore, in large meta-analyses of OA, no association with sequence variants at the CRTAC1 locus has been described (14,15). This indicates that the observed increase in CRTAC1 in plasma was unlikely to cause OA. Moreover, adjusting the CRTAC1 association for the OA polygenic risk scores, which integrated a large fraction of the contribution of genetics to the disease, had little impact on CRTAC1 association with OA (Supplementary Figure 5, http://onlin elibr ary. wiley.com/doi/10.1002/art.41793/ abstract). This suggests that CRTAC1 captured a different component of the disease.
Models for OA association and prediction of progression to joint replacement. Although a single, easily interpreted biomarker is currently most feasible to use in clinical settings, a model incorporating other factors besides CRTAC1 alone may provide more accurate association/prediction of OA development and progression to joint replacement. We therefore aimed to build the best available association model for OA diagnosis and for prediction of progression to joint replacement by including all independently associated proteins and the polygenic risk scores for OA in the model, in addition to age, sex, and BMI. Exclusion/ inclusion criteria were the same as before, with the exception that those without genotypes were also excluded. To avoid overfitting of the models, we split the data set in two and used three-fourths of the data to build the model and one-fourth to test it. We performed backward stepwise selection of associated variables (i.e., proteins, polygenic risk score, age, sex, and BMI) in a single logistic regression model and likewise, to determine a set of proteins that predict OA progression to joint replacement, we used Cox proportional hazards regression. We demonstrated both a good fit between the training and test sets and adequate power of our study to detect association with OA (Supplementary Figure 6, http://onlin elibr ary.wiley.com/doi/10.1002/art.41793/ abstract).
To evaluate the discriminatory power of our models, we applied a receiver operating characteristic curve and calculated the AUC. In the full model, the AUC was 0.762 for hip OA, 0.770 for knee OA, and 0.800 for hand OA, and the cumulative AUC estimate for joint replacement was 0.661 for hip replacement and 0.734 for knee replacement (Supplementary Table 9 The polygenic risk score and the classic risk factors, age and sex, contributed to all of the models. BMI, however, did not contribute to the hip replacement model. As expected, given the composition of the sample set, classic risk factors contributed the most in all of the models. However, adding the proteins and polygenic risk scores significantly improved the association of the models with OA (P = 1.5 × 10 −101 for hand OA, P = 3.1 × 10 −75 for hip OA, and P = 3.8 × 10 −203 for knee OA) and prediction of joint replacement (P = 4.8 × 10 −111 for knee replacement and P = 9.3 × 10 −71 for hip replacement) (Supplementary Table 9).
These OA risk profiles were strongly associated with joint pain, and the hand risk score was strongly associated with hand OA severity (Supplementary Table 10 Table 11, http://onlin elibr ary. wiley.com/doi/10.1002/art.41793/ abstract) and were, therefore, not specific for OA. Furthermore, use of the full models required measuring up to 35 proteins for both OA association and prediction of joint replacement, and genotypes for the polygenic risk score, which would be challenging in clinical settings.
As described above regarding CRTAC1, we also analyzed whether the levels of the plasma proteins selected in the association and prediction models could cause OA by testing protein QTL sequence variants (n = 267) (Supplementary Note, Supplementary Tables 12 and 13, and Supplementary Figure 5, http://onlin elibr ary. wiley.com/doi/10.1002/art.41793/ abstract). As with CRTAC1, our results indicated that plasma levels of these proteins are unlikely to cause OA, but their altered levels are consequences of the disease.
Induction of a specific protein profile by joint replacement. We observed a striking difference in levels of OAassociated proteins between patients who had undergone a joint replacement before plasma collection (n = 1,452) and those who had not (n = 4,229). In our data, a prior joint replacement was associated with 137 plasma proteins (Supplementary Table 14, http://onlin elibr ary.wiley.com/doi/10.1002/art.41793/ abstract), a minority of which were associated with OA in patients who had not undergone joint replacement (n = 16 proteins). This indicated that the joint replacement itself was associated with a plasma protein profile that was a consequence of the joint implant rather than the disease. Therefore, we excluded patients who had undergone joint replacement prior to plasma collection from all the OA association studies described above. The protein with the strongest association with prior joint replacement, CUB domain-containing protein 1 (CDCP1) (P = 1.1 × 10 −198 , OR 2.67), was not associated with OA without joint replacement. Like most proteins associated with joint replacement but not with OA, CDCP1 levels changed very little during the 20 years preceding joint replacement, followed by an increase in levels at the time of arthroplasty (β = 0.866, P = 5.7 × 10 −16 ) and a continued steady increase thereafter (β = 0.072 per year, P = 8.7 × 10 −43 ), as demonstrated in data from different individuals at different time points after joint replacement. In contrast, the proteins that were also associated with OA demonstrated an increase in levels prior to joint replacement that either continued (e.g., ASB9) or decreased (e.g., CRTAC1) after surgery (Supplementary Figure 11

DISCUSSION
In this study of a large panel of plasma proteins (4,792 measured), we determined that plasma CRTAC1 levels were the most strongly associated with OA diagnosis, and CRTAC1 was the best predictor of OA progression to joint replacement. Patients with OA who were in the highest quintile of risk based on known risk factors (i.e., age, sex, and BMI) and plasma CRTAC1 levels were 16 times more likely to undergo knee replacement within 5 years, and 6.5 times more likely to undergo hip replacement than those in the lowest quintile. Furthermore, CRTAC1 levels in plasma were also associated with joint pain and hand OA severity. Importantly, CRTAC1 levels are not associated with other inflammatory joint diseases such as RA, gout, or PsA. We determined that the protein panels of up to 27 proteins, together with the OA polygenic risk scores (in addition to age, sex, and BMI) in the models, performed better than CRTAC1 alone (in combination with age, sex, and BMI). However, these expanded associations and prediction models were not specific to OA. This highlights CRTAC1 as a promising biomarker for OA, since we found that, among all of the disease subphenotypes studied, this protein was not only strongly associated with OA, but also specific for OA with no association with other inflammatory joint diseases. We also | 2033 demonstrated that CRTAC1 is a risk factor for OA independent of age, sex, and BMI. Based on CRTAC1 QTL genetic analysis, we demonstrated that CRTAC1 levels in plasma are unlikely to cause OA.
The function of CRTAC1 is not fully understood. It is a glycosylated, calcium-binding extracellular matrix protein that is highly conserved from bacteria to humans (21), that was initially identified as a marker of chondrocyte development from mesenchymal stem cells (22). This protein is found in greater concentrations in knee synovial fluid from OA patients compared to controls (23), and also found at higher levels in OA cartilage (24). Furthermore, CRTAC1 is also involved in apoptosis and pyroptosis of human lens epithelial cells (25,26).
In addition to these findings, we observed that joint replacement was associated with a protein profile in plasma that lacks association with OA and persists long after joint replacement surgery. This profile may therefore reflect the body's reaction to the joint implant. The plasma levels of these proteins could be candidates for prediction of prosthesis survival time or early prosthesis failure.
To our knowledge, this study comprises the largest biomarker study of OA to date. We simultaneously investigated 4,792 plasma proteins in 39,155 individuals of whom 12,178 had OA. This study thus has greatly increased power over previous studies to identify circulating biomarkers for OA. We also demonstrated that plasma levels of COMP, one of the most investigated biomarkers for OA to date (9)(10)(11)(12), were associated with hand, hip, and knee OA and predicted knee replacement. Thus, our study supports previous findings indicating that COMP is a likely biomarker for OA. However, we showed that COMP was a substantially weaker predictor than CRTAC1, and because it is also associated with weight, it is a less suitable biomarker. We did not have access to measurements of the other much-studied biomarker for OA, urinary CTX-II, which is a degradation product of type II collagen found in urine. We therefore could not test whether urinary CTX-II correlated with CRTAC1 or whether it was independently associated with OA.
Limitations of this study include the lack of information on age at the onset of OA symptoms or age at first confirmed diagnosis, which reflected the registry nature of the data. Thus, we were unable to determine whether the CRTAC1 levels were already altered at a very early stage of the disease. Furthermore, in the association analysis, we included all patients diagnosed as having OA by March 2020, irrespective of when their plasma sample was collected. This group is therefore heterogeneous, i.e., at the time the plasma sample was collected, some patients had advanced disease whereas others were symptom-free. Moreover, we did not have access to other population studies for replication of our findings. Follow-up studies in other populations of various ancestries, in addition to studies addressing technical issues, are extremely important for validation of CRTAC1 as a general biomarker for OA.
In this large hypothesis-free study on OA and joint replacement, we identified a novel biomarker, CRTAC1, with a strong, and specific, effect on both disease risk and disease severity. It had a stronger effect than biomarkers that were previously established as having the best predictive abilities. If validated in further studies, CRTAC1 will be a promising biomarker both for clinical use and for therapeutic development, where there is currently a large unmet need.