The limited molecular classifications and disease signatures of osteoarthritis (OA) impede the development of prediagnosis and targeted therapeutics for OA patients. To classify and understand the subtypes of OA, we collected three types of tissue including cartilage, subchondral bone, and synovium from multiple clinical centers and constructed an extensive transcriptome atlas of OA patients. By applying unsupervised clustering analysis to the cartilage transcriptome, OA patients were classified into four subtypes with distinct molecular signatures: a glycosaminoglycan metabolic disorder subtype (C1), a collagen metabolic disorder subtype (C2), an activated sensory neuron subtype (C3), and an inflammation subtype (C4). Through ligand-receptor crosstalk analysis of the three knee tissue types, we linked molecular functions with the clinical symptoms of different OA subtypes. For example, the Gene Ontology functional term of vasculature development was enriched in the subchondral bone-cartilage crosstalk of C2 and the cartilage-subchondral bone crosstalk of C4, which might lead to severe osteophytes in C2 patients and apparent joint space narrowing in C4 patients. Based on the marker genes of the four OA subtypes identified in this study, we modeled OA subtypes with two independent published RNA-seq datasets through random forest classification. The findings of this work contradicted traditional OA diagnosis by medical imaging and revealed distinct molecular subtypes in knee OA patients, which may allow for precise diagnosis and treatment of OA.

骨关节炎（OA）的有限分子分类和疾病特征阻碍了OA患者的预诊断和靶向治疗的发展。为了分类和了解OA的亚型，我们从多个临床中心收集了三种类型的组织，包括软骨，软骨下骨和滑膜，并构建了OA患者的广泛转录组图谱。通过对软骨转录组进行无监督聚类分析，将OA患者分为具有不同分子特征的四种亚型：糖胺聚糖代谢异常亚型（C1），胶原代谢异常亚型（C2），激活的感觉神经元亚型（C3）和炎症亚型（C4）。通过对三种膝盖组织类型的配体-受体串扰分析，我们将分子功能与不同OA亚型的临床症状联系起来。例如，脉管系统发育的基因本体论功能术语丰富了C2的软骨下骨软骨串扰和C4的软骨-软骨下骨串扰，这可能导致C2患者的骨赘严重，而C4患者的关节间隙明显缩小。基于本研究中鉴定的四种OA亚型的标记基因，我们通过随机森林分类，使用两个独立出版的RNA-seq数据集对OA亚型进行建模。这项工作的发现与医学影像学对传统OA的诊断相矛盾，并且揭示了膝OA患者的独特分子亚型，这可能有助于OA的精确诊断和治疗。C2的软骨下骨软骨串扰和C4的软骨-软骨下骨串扰丰富了脉管系统发育的基因本体学功能术语，这可能导致C2患者的骨赘严重，C4患者的关节间隙明显缩小。基于本研究中鉴定的四种OA亚型的标记基因，我们通过随机森林分类，使用两个独立出版的RNA-seq数据集对OA亚型进行建模。这项工作的发现与医学影像学对传统OA的诊断相矛盾，并且揭示了膝OA患者的独特分子亚型，这可能有助于OA的精确诊断和治疗。C2的软骨下骨软骨串扰和C4的软骨-软骨下骨串扰丰富了脉管系统发育的基因本体学功能术语，这可能导致C2患者的骨赘严重，C4患者的关节间隙明显缩小。基于本研究中鉴定的四种OA亚型的标记基因，我们通过随机森林分类，使用两个独立出版的RNA-seq数据集对OA亚型进行建模。这项工作的发现与医学影像学对传统OA的诊断相矛盾，并且揭示了膝OA患者的独特分子亚型，这可能有助于OA的精确诊断和治疗。这可能导致C2患者出现严重的骨赘，而C4患者的关节间隙明显缩小。基于本研究中鉴定的四种OA亚型的标记基因，我们通过随机森林分类，使用两个独立出版的RNA-seq数据集对OA亚型进行建模。这项工作的发现与医学影像学对传统OA的诊断相矛盾，并揭示了膝OA患者的独特分子亚型，这可能有助于OA的精确诊断和治疗。这可能导致C2患者出现严重的骨赘，而C4患者的关节间隙明显缩小。基于本研究中鉴定的四种OA亚型的标记基因，我们通过随机森林分类，使用两个独立出版的RNA-seq数据集对OA亚型进行建模。这项工作的发现与医学影像学对传统OA的诊断相矛盾，并且揭示了膝OA患者的独特分子亚型，这可能有助于OA的精确诊断和治疗。