{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,7,30]],"date-time":"2026-07-30T15:42:56Z","timestamp":1785426176350,"version":"3.56.0"},"reference-count":38,"publisher":"MDPI AG","issue":"6","license":[{"start":{"date-parts":[[2026,6,15]],"date-time":"2026-06-15T00:00:00Z","timestamp":1781481600000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100005090","name":"Beijing Nova Program","doi-asserted-by":"crossref","award":["H020821500190"],"award-info":[{"award-number":["H020821500190"]}],"id":[{"id":"10.13039\/501100005090","id-type":"DOI","asserted-by":"crossref"}]},{"DOI":"10.13039\/501100005090","name":"Beijing Nova Program","doi-asserted-by":"crossref","award":["Z20230484422"],"award-info":[{"award-number":["Z20230484422"]}],"id":[{"id":"10.13039\/501100005090","id-type":"DOI","asserted-by":"crossref"}]},{"award":["H020821500190"],"award-info":[{"award-number":["H020821500190"]}],"id":[{"id":"https:\/\/ror.org\/033pgsq69","id-type":"ROR","asserted-by":"publisher"}]},{"award":["Z20230484422"],"award-info":[{"award-number":["Z20230484422"]}],"id":[{"id":"https:\/\/ror.org\/033pgsq69","id-type":"ROR","asserted-by":"publisher"}]},{"name":"Guizhou Provincial Clinical Special Project","award":["Qian Ke He Cheng Guo-LC [2026] General 299"],"award-info":[{"award-number":["Qian Ke He Cheng Guo-LC [2026] General 299"]}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Genes"],"abstract":"<jats:p>Background: Osteoarthritis (OA) is a heterogeneous joint disease characterized by cartilage degeneration. The interplay between extracellular matrix (ECM) remodeling, endoplasmic reticulum (ER) stress, and inflammatory signaling in OA pathogenesis remains incompletely understood. This study aimed to identify robust diagnostic biomarkers and explore the mechanistic convergence of key genes in OA cartilage through an integrated transcriptomic framework. Methods: Three independent cartilage transcriptomic datasets (GSE285234, GSE287861, GSE289464) were integrated after ComBat batch correction. Differentially expressed genes (DEGs) were identified using limma, followed by ORA and GSEA for functional enrichment. LASSO logistic regression identified hub genes for a diagnostic model and nomogram, validated by leave-one-out cross-validation (LOOCV). Consensus clustering stratified OA samples into molecular subtypes. Single-cell RNA-sequencing (scRNA-seq) data (GSE169454, GSE220243) were used to validate cell-type-specific expression. Virtual gene knockout (scTenifoldKnk) and pathway analysis inferred downstream functional consequences. Results: Fifty-eight DEGs (predominantly downregulated) were enriched in ECM and ER protein processing pathways. Six hub genes (EIF2S1, GANAB, STT3A, XBP1, MGP, PMP22) showed robust selection stability. The diagnostic model achieved a LOOCV AUC of 0.769, a well-calibrated nomogram, and superior net benefit. Unsupervised clustering revealed two OA subtypes with divergent unfolded protein response (UPR) and TGF-\u03b2 pathway activities. scRNA-seq confirmed hub gene expression in chondrocytes and other joint microenvironment cells. Notably, virtual knockout of five hub genes convergently perturbed IL-17, NF-\u03baB, and chemokine signaling pathways. Conclusions: This study identified and validated a six-gene signature reflecting ECM-ER-inflammatory crosstalk in OA cartilage. The convergent perturbation of inflammatory pathways by functionally distinct hub genes reveals a mechanistic core that may serve as a diagnostic panel and a platform for targeted therapeutic investigation in OA.<\/jats:p>","DOI":"10.3390\/genes17060696","type":"journal-article","created":{"date-parts":[[2026,6,17]],"date-time":"2026-06-17T05:56:31Z","timestamp":1781675791000},"page":"696","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":0,"title":["Integrative Transcriptomic Analysis and Single-Cell Validation Identify a Six-Hub-Gene Signature Converging on Inflammatory Signaling in Osteoarthritis"],"prefix":"10.3390","volume":"17","author":[{"ORCID":"https:\/\/orcid.org\/0009-0009-7957-6849","authenticated-orcid":false,"given":"Xueya","family":"Lv","sequence":"first","affiliation":[{"name":"Department of Orthopedics Medicine, Beijing Jishuitan Hospital, Capital Medical University, Beijing 100035, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Yang","family":"Yu","sequence":"additional","affiliation":[{"name":"Department of Orthopedics, Army No. 82 Group Military Hospital, Baoding 071000, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Jiawen","family":"Fan","sequence":"additional","affiliation":[{"name":"Department of Orthopedics, Army No. 82 Group Military Hospital, Baoding 071000, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Lianjiang","family":"Guo","sequence":"additional","affiliation":[{"name":"Department of Orthopedics, Army No. 82 Group Military Hospital, Baoding 071000, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Xiang","family":"Zhu","sequence":"additional","affiliation":[{"name":"Department of Oncology, Army No. 82 Group Military Hospital, Baoding 071000, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Xingye","family":"Li","sequence":"additional","affiliation":[{"name":"Department of Spine Surgery, Beijing Jishuitan Hospital, Capital Medical University, Beijing 100035, China"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"1968","published-online":{"date-parts":[[2026,6,15]]},"reference":[{"key":"ref_1","unstructured":"GBD 2021 Osteoarthritis Collaborators (2023). Global, regional, and national burden of osteoarthritis, 1990\u20132020 and projections to 2050: A systematic analysis for the Global Burden of Disease Study 2021. Lancet Rheumatol., 5, e508\u2013e522."},{"key":"ref_2","doi-asserted-by":"crossref","first-page":"272","DOI":"10.1080\/03008207.2025.2502591","article-title":"Mechanical insights into fat pads: A comparative study of infrapatellar and suprapatellar fat pads in osteoarthritis","volume":"66","author":"Pettenuzzo","year":"2025","journal-title":"Connect. Tissue Res."},{"key":"ref_3","doi-asserted-by":"crossref","first-page":"2013","DOI":"10.1016\/j.joca.2016.07.017","article-title":"Disease-modifying treatments for osteoarthritis (DMOADs) of the knee and hip: Lessons learned from failures and opportunities for the future","volume":"24","author":"Karsdal","year":"2016","journal-title":"Osteoarthr. Cartil."},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"423","DOI":"10.1136\/annrheumdis-2017-212603","article-title":"Stratification of knee osteoarthritis: Two major patient subgroups identified by genome-wide expression analysis of articular cartilage","volume":"77","author":"Soul","year":"2018","journal-title":"Ann. Rheum. Dis."},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"118","DOI":"10.1093\/biostatistics\/kxj037","article-title":"Adjusting batch effects in microarray expression data using empirical Bayes methods","volume":"8","author":"Johnson","year":"2007","journal-title":"Biostatistics"},{"key":"ref_6","doi-asserted-by":"crossref","unstructured":"Hughes, A., Oxford, A.E., Tawara, K., Jorcyk, C.L., and Oxford, J.T. (2017). Endoplasmic reticulum stress and unfolded protein response in cartilage pathophysiology; contributing factors to apoptosis and osteoarthritis. Int. J. Mol. Sci., 18.","DOI":"10.3390\/ijms18030665"},{"key":"ref_7","doi-asserted-by":"crossref","first-page":"421","DOI":"10.1038\/s41580-020-0250-z","article-title":"Mechanisms, regulation and functions of the unfolded protein response","volume":"21","author":"Hetz","year":"2020","journal-title":"Nat. Rev. Mol. Cell Biol."},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"33","DOI":"10.1038\/nrrheum.2010.196","article-title":"Role of proinflammatory cytokines in the pathophysiology of osteoarthritis","volume":"7","author":"Kapoor","year":"2011","journal-title":"Nat. Rev. Rheumatol."},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"35","DOI":"10.1038\/nri.2017.76","article-title":"Single-cell RNA sequencing to explore immune cell heterogeneity","volume":"18","author":"Papalexi","year":"2018","journal-title":"Nat. Rev. Immunol."},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"100","DOI":"10.1136\/annrheumdis-2017-212863","article-title":"Single-cell RNA-seq analysis reveals the progression of human osteoarthritis","volume":"78","author":"Ji","year":"2019","journal-title":"Ann. Rheum. Dis."},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"23","DOI":"10.1038\/s42003-025-07460-x","article-title":"KLF9-GRK5-HDAC6 axis aggravates osteoarthritis pathogenesis by promoting chondrocyte extracellular matrix degradation and apoptosis","volume":"8","author":"Zhou","year":"2025","journal-title":"Commun. Biol."},{"key":"ref_12","doi-asserted-by":"crossref","unstructured":"Makarczyk, M.J., Zhang, Y., Aguglia, A., Bartholomew, O., Hines, S., Li, K., Sinkar, S., Liu, S., Duvall, C., and Lin, H. (2025). Aging-associated increase of GATA4 levels in articular cartilage is linked to impaired regenerative capacity of chondrocytes and osteoarthritis. bioRxiv.","DOI":"10.7554\/eLife.106224.1"},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"2120","DOI":"10.1007\/s10528-025-11131-1","article-title":"Analysis and ceRNA network construction of differentially expressed lncRNAs and mRNAs in human osteoarthritis cartilage","volume":"64","author":"Zhou","year":"2026","journal-title":"Biochem. Genet."},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"882","DOI":"10.1093\/bioinformatics\/bts034","article-title":"The sva package for removing batch effects and other unwanted variation in high-throughput experiments","volume":"28","author":"Leek","year":"2012","journal-title":"Bioinformatics"},{"key":"ref_15","doi-asserted-by":"crossref","first-page":"e47","DOI":"10.1093\/nar\/gkv007","article-title":"limma powers differential expression analyses for RNA-sequencing and microarray studies","volume":"43","author":"Ritchie","year":"2015","journal-title":"Nucleic Acids Res."},{"key":"ref_16","doi-asserted-by":"crossref","unstructured":"Lin, Y., Golovnina, K., Chen, Z.X., Lee, H.N., Negron, Y.L., Sultana, H., Oliver, B., and Harbison, S.T. (2016). Comparison of normalization and differential expression analyses using RNA-Seq data from 726 individual Drosophila melanogaster. BMC Genom., 17.","DOI":"10.1186\/s12864-015-2353-z"},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"284","DOI":"10.1089\/omi.2011.0118","article-title":"clusterProfiler: An R package for comparing biological themes among gene clusters","volume":"16","author":"Yu","year":"2012","journal-title":"OMICS"},{"key":"ref_18","doi-asserted-by":"crossref","first-page":"15545","DOI":"10.1073\/pnas.0506580102","article-title":"Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles","volume":"102","author":"Subramanian","year":"2005","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"ref_19","doi-asserted-by":"crossref","first-page":"1","DOI":"10.18637\/jss.v033.i01","article-title":"Regularization paths for generalized linear models via coordinate descent","volume":"33","author":"Friedman","year":"2010","journal-title":"J. Stat. Softw."},{"key":"ref_20","doi-asserted-by":"crossref","first-page":"565","DOI":"10.1177\/0272989X06295361","article-title":"Decision curve analysis: A novel method for evaluating prediction models","volume":"26","author":"Vickers","year":"2006","journal-title":"Med. Decis. Mak."},{"key":"ref_21","doi-asserted-by":"crossref","first-page":"1572","DOI":"10.1093\/bioinformatics\/btq170","article-title":"ConsensusClusterPlus: A class discovery tool with confidence assessments and item tracking","volume":"26","author":"Wilkerson","year":"2010","journal-title":"Bioinformatics"},{"key":"ref_22","doi-asserted-by":"crossref","unstructured":"H\u00e4nzelmann, S., Castelo, R., and Guinney, J. (2013). GSVA: Gene set variation analysis for microarray and RNA-seq data. BMC Bioinform., 14.","DOI":"10.1186\/1471-2105-14-7"},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"6105","DOI":"10.2147\/JIR.S385648","article-title":"Single-Cell RNA Sequencing Reveals Transcriptional Changes in the Cartilage of Subchondral Insufficiency Fracture of the Knee","volume":"15","author":"Tang","year":"2022","journal-title":"J. Inflamm. Res."},{"key":"ref_24","doi-asserted-by":"crossref","first-page":"403","DOI":"10.1136\/ard-2022-223227","article-title":"Senescent cell population with ZEB1 transcription factor as its main regulator promotes osteoarthritis in cartilage and meniscus","volume":"82","author":"Swahn","year":"2023","journal-title":"Ann. Rheum. Dis."},{"key":"ref_25","doi-asserted-by":"crossref","first-page":"3573","DOI":"10.1016\/j.cell.2021.04.048","article-title":"Integrated analysis of multimodal single-cell data","volume":"184","author":"Hao","year":"2021","journal-title":"Cell"},{"key":"ref_26","doi-asserted-by":"crossref","first-page":"163","DOI":"10.1038\/s41590-018-0276-y","article-title":"Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage","volume":"20","author":"Aran","year":"2019","journal-title":"Nat. Immunol."},{"key":"ref_27","first-page":"1006","article-title":"scTenifoldKnk: An efficient virtual knockout tool for gene function predictions via single-cell gene regulatory network perturbation","volume":"13","author":"Osorio","year":"2022","journal-title":"Nat. Commun."},{"key":"ref_28","doi-asserted-by":"crossref","first-page":"224","DOI":"10.1186\/ar2592","article-title":"Cartilage homeostasis in health and rheumatic diseases","volume":"11","author":"Goldring","year":"2009","journal-title":"Arthritis Res. Ther."},{"key":"ref_29","doi-asserted-by":"crossref","first-page":"271","DOI":"10.1038\/16729","article-title":"Protein translation and folding are coupled by an endoplasmic-reticulum-resident kinase","volume":"397","author":"Harding","year":"1999","journal-title":"Nature"},{"key":"ref_30","doi-asserted-by":"crossref","first-page":"53","DOI":"10.1016\/j.molcel.2007.06.011","article-title":"XBP1 controls diverse cell type- and condition-specific transcriptional regulatory networks","volume":"27","author":"Zhou","year":"2007","journal-title":"Mol. Cell"},{"key":"ref_31","doi-asserted-by":"crossref","first-page":"57","DOI":"10.1016\/j.ceb.2016.03.021","article-title":"N-linked glycosylation and homeostasis of the endoplasmic reticulum","volume":"41","author":"Cherepanova","year":"2016","journal-title":"Curr. Opin. Cell Biol."},{"key":"ref_32","doi-asserted-by":"crossref","first-page":"142","DOI":"10.1038\/5102","article-title":"Mutations in the gene encoding the human matrix Gla protein cause Keutel syndrome","volume":"21","author":"Munroe","year":"1999","journal-title":"Nat. Genet."},{"key":"ref_33","doi-asserted-by":"crossref","first-page":"673","DOI":"10.1007\/s12035-012-8370-x","article-title":"The PMP22 gene and its related diseases","volume":"47","author":"Li","year":"2013","journal-title":"Mol. Neurobiol."},{"key":"ref_34","doi-asserted-by":"crossref","unstructured":"Zhang, W., Robertson, W.B., Zhao, J., Chen, W., and Xu, J. (2019). Emerging Trend in the Pharmacotherapy of Osteoarthritis. Front. Endocrinol., 10.","DOI":"10.3389\/fendo.2019.00431"},{"key":"ref_35","doi-asserted-by":"crossref","first-page":"9","DOI":"10.1038\/s41413-019-0047-x","article-title":"Pathological mechanisms and therapeutic outlooks for arthrofibrosis","volume":"7","author":"Usher","year":"2019","journal-title":"Bone Res."},{"key":"ref_36","doi-asserted-by":"crossref","first-page":"128","DOI":"10.1097\/EDE.0b013e3181c30fb2","article-title":"Assessing the performance of prediction models: A framework for traditional and novel measures","volume":"21","author":"Steyerberg","year":"2010","journal-title":"Epidemiology"},{"key":"ref_37","doi-asserted-by":"crossref","first-page":"10868","DOI":"10.1038\/s41598-020-67730-y","article-title":"Synovial cell cross-talk with cartilage plays a major role in the pathogenesis of osteoarthritis","volume":"10","author":"Chou","year":"2020","journal-title":"Sci. Rep."},{"key":"ref_38","doi-asserted-by":"crossref","first-page":"748","DOI":"10.1136\/annrheumdis-2016-209757","article-title":"Systemic inhibition of IL-6\/Stat3 signalling protects against experimental osteoarthritis","volume":"76","author":"Latourte","year":"2017","journal-title":"Ann. Rheum. Dis."}],"container-title":["Genes"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/2073-4425\/17\/6\/696\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,6,19]],"date-time":"2026-06-19T04:37:13Z","timestamp":1781843833000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/2073-4425\/17\/6\/696"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,6,15]]},"references-count":38,"journal-issue":{"issue":"6","published-online":{"date-parts":[[2026,6]]}},"alternative-id":["genes17060696"],"URL":"https:\/\/doi.org\/10.3390\/genes17060696","relation":{},"ISSN":["2073-4425"],"issn-type":[{"value":"2073-4425","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,6,15]]}}}