Identification of key biomarkers and immune infiltration in the synovial tissue of osteoarthritis by bioinformatics analysis.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:31988808", "kind": "science-record", "revision": 7, "content_hash": "sha256:8bf8d962be168f01592f9dd81772880865f4fa44a3ecb2928eb64f5e8460194e", "created_at": "2026-11-01T15:40:55Z", "updated_at": "2026-11-01T22:01:13Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Identification of key biomarkers and immune infiltration in the synovial tissue of osteoarthritis by bioinformatics analysis.", "doi": "10.7717/peerj.8390", "pmid": "31988808", "authors": ["Weisong Cai", "Haohuan Li", "Yubiao Zhang", "Guangtao Han"], "journal": "PeerJ", "year": 2020, "abstract": "BACKGROUND: Osteoarthritis (OA) is the most common chronic degenerative joint disease and is mainly characterized by cartilage degeneration, subcartilage bone hyperplasia, osteophyte formation and joint space stenosis. Recent studies showed that synovitis might also be an important pathological change of OA. However, the molecular mechanisms of synovitis in OA are still not well understood.\nOBJECTIVE: This study was designed to identify key biomarkers and immune infiltration in the synovial tissue of osteoarthritis by bioinformatics analysis.\nMATERIALS AND METHODS: The gene expression profiles of GSE12021, GSE55235 and GSE55457 were downloaded from the GEO database. The differentially expressed genes (DEGs) were identified by the LIMMA package in Bioconductor, and functional enrichment analyses were performed. A protein-protein interaction network (PPI) was constructed, and module analysis was performed using STRING and Cytoscape. The CIBERSORT algorithm was used to analyze the immune infiltration of synovial tissue between OA and normal controls.\nRESULTS: A total of 106 differentially expressed genes, including 68 downregulated genes and 38 upregulated genes, were detected. The PPI network was assessed, and the most significant module containing 14 hub genes was identified. Gene Ontology analysis revealed that the hub genes were significantly enriched in immune cell chemotaxis and cytokine activity. KEGG pathway analysis showed that the hub genes were significantly enriched in the rheumatoid arthritis signaling pathway, IL-17 signaling pathway and cytokine-cytokine receptor interaction signaling pathway. The immune infiltration profiles varied significantly between osteoarthritis and normal controls. Compared with normal tissue, OA synovial tissue contained a higher proportion of memory B cells, naive CD4+ T cells, regulatory T cells, resting dendritic cells and resting mast cells, while naive CD4+ T cells, activated NK cells, activated mast cells and eosinophils contributed to a relatively lower portion (P > 0.05). Finally, the expression levels of 11 hub genes were confirmed by RT-PCR.\nCONCLUSION: The hub genes and the difference in immune infiltration in synovial tissue between osteoarthritis and normal controls might provide new insight for understanding OA development.", "claims": [{"id": "claim-method", "statement": "The corresponding P value of the gene symbols after t test were used, and adjusted P < 0.05 and —logFC—>2 were used as the selection criteria.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:31988808", "locator": "sections[4] \"Identification of DEGs\" paragraph 2", "quote": "The corresponding P value of the gene symbols after t test were used, and adjusted P < 0.05 and —logFC—>2 were used as the selection criteria."}], "state": "complete"}, {"id": "claim-data-access", "statement": "The data is available at NCBI GEO: GSE12021, GSE55235 and GSE55457. The code is available as Supplementary Files.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:31988808", "locator": "sections[19] \"Additional Information and Declarations\" paragraph 5", "quote": "The data is available at NCBI GEO: GSE12021, GSE55235 and GSE55457. The code is available as Supplementary Files."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "complete", "objective": null, "steps": [{"id": "read_published_estimate", "order": 1, "name": "Read the published estimate from the source", "kind": "other", "state": "complete", "blocked_on": null, "description": "The claim, method and evidence span published with the source article.", "duration_hours": null, "parameters": {}, "inputs": ["source:pmid:31988808"], "outputs": ["the published estimate and its evidence"], "instruments": [], "reagents": []}, {"id": "obtain_source_data", "order": 2, "name": "Obtain the source data from the named repository", "kind": "data-preparation", "state": "complete", "blocked_on": null, "description": "The GEO series accession resolves and its sample accessions are listed. This lane did not download the dataset or read its licence.", "duration_hours": null, "parameters": {}, "inputs": ["data:geo:GSE12021"], "outputs": ["the source data with its repository provenance"], "instruments": [], "reagents": []}, {"id": "run_analysis", "order": 3, "name": "Run the published analysis", "kind": "analysis", "state": "complete", "blocked_on": null, "description": "Apply the paper's statistical method to the resolved series and compare the result with the published target.", "duration_hours": null, "parameters": {}, "inputs": ["data:geo:GSE12021"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": []}, "data": {"state": "partial", "assets": [{"id": "source:pmid:31988808", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.7717/peerj.8390", "sha256": "sha256:630dabe1af6d60c2bac945b0836a295903e9ed05d5ac1d6aab2ea6e953538590", "bytes": 28279, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.7717/peerj.8390", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC6970550", "source": "europe_pmc_fulltext", "coverage": "whole_body_without_references", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE12021", "role": "raw", "format": "geo-series", "media_type": "application/vnd.ncbi.geo-series+json", "uri": "https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE12021", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE12021", "version": null, "note": "The data is available at NCBI GEO: GSE12021, GSE55235 and GSE55457. The code is available as Supplementary Files.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The data is available at NCBI GEO: GSE12021, GSE55235 and GSE55457. The code is available as Supplementary Files.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM302859", "GSM302864", "GSM302866", "GSM302870", "GSM302872", "GSM302874", "GSM302876", "GSM302878", "GSM302880", "GSM302881", "GSM302882", "GSM302895", "GSM302896", "GSM302897", "GSM302900", "GSM302905", "GSM302924", "GSM302926", "GSM302927", "GSM302929", "GSM302930", "GSM302932", "GSM302933", "GSM302934", "GSM302943", "GSM302944", "GSM303326", "GSM303339", "GSM303341", "GSM303355", "GSM303356", "GSM303357", "GSM303358", "GSM303359", "GSM303360", "GSM303361", "GSM303362", "GSM303363", "GSM303364", "GSM303365", "GSM303366", "GSM303367", "GSM303368", "GSM303369", "GSM303370", "GSM303377", "GSM303520", "GSM303521", "GSM303522", "GSM303523", "GSM303524", "GSM303525", "GSM303530", "GSM303531", "GSM303532", "GSM303533", "GSM303534"], "sample_count": 57, "group_sizes": null, "group_sizes_basis": null, "platform": "96;97", "taxon": "Homo sapiens", "title": "Identification of inter-individual and gene-specific variances in mRNA expression profiles in the RA SM", "resolution_checked_at": "2026-11-01T15:39:54Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}, {"id": "data:geo:GSE55235", "role": "raw", "format": "geo-series", "media_type": "application/vnd.ncbi.geo-series+json", "uri": "https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE55235", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE55235", "version": null, "note": "The data is available at NCBI GEO: GSE12021, GSE55235 and GSE55457. The code is available as Supplementary Files.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The data is available at NCBI GEO: GSE12021, GSE55235 and GSE55457. The code is available as Supplementary Files.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM1332201", "GSM1332202", "GSM1332203", "GSM1332204", "GSM1332205", "GSM1332206", "GSM1332207", "GSM1332208", "GSM1332209", "GSM1332210", "GSM1332211", "GSM1332212", "GSM1332213", "GSM1332214", "GSM1332215", "GSM1332216", "GSM1332217", "GSM1332218", "GSM1332219", "GSM1332220", "GSM1332221", "GSM1332222", "GSM1332223", "GSM1332224", "GSM1332225", "GSM1332226", "GSM1332227", "GSM1332228", "GSM1332229", "GSM1332230"], "sample_count": 30, "group_sizes": null, "group_sizes_basis": null, "platform": "96", "taxon": "Homo sapiens", "title": "Identification of rheumatoid arthritis and osteoarthritis patients by transcriptome-based rule set generation", "resolution_checked_at": "2026-11-01T15:39:54Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}, {"id": "data:geo:GSE55457", "role": "raw", "format": "geo-series", "media_type": "application/vnd.ncbi.geo-series+json", "uri": "https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE55457", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE55457", "version": null, "note": "The data is available at NCBI GEO: GSE12021, GSE55235 and GSE55457. The code is available as Supplementary Files.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The data is available at NCBI GEO: GSE12021, GSE55235 and GSE55457. The code is available as Supplementary Files.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM1337304", "GSM1337305", "GSM1337306", "GSM1337307", "GSM1337308", "GSM1337309", "GSM1337310", "GSM1337311", "GSM1337312", "GSM1337313", "GSM1337314", "GSM1337315", "GSM1337316", "GSM1337317", "GSM1337318", "GSM1337319", "GSM1337320", "GSM1337321", "GSM1337322", "GSM1337323", "GSM1337324", "GSM1337325", "GSM1337326", "GSM1337327", "GSM1337328", "GSM1337329", "GSM1337330", "GSM1337331", "GSM1337332", "GSM1337333", "GSM1337334", "GSM1337335", "GSM1337336"], "sample_count": 33, "group_sizes": null, "group_sizes_basis": null, "platform": "96", "taxon": "Homo sapiens", "title": "Identification of rheumatoid arthritis and osteoarthritis patients by transcriptome-based rule set generation [Jena]", "resolution_checked_at": "2026-11-01T15:39:54Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}, {"id": "asset:table:pmid:31988808:1", "role": "table", "format": "jats-table", "media_type": "application/vnd.jats.table+json", "uri": null, "sha256": "sha256:0d4e80806540177477f6b8112b9b4ad8108e940c64c03e7640d4e4f16e893550", "bytes": 265, "access": "open", "state": "verified", "reader": "read_table", "reader_status": "available", "identifier": "section[11]/table-wrap[1]#table-1", "version": null, "note": "GO analyses results of DEGs (top 10 according to adjusted P value).“Count” means how many DEGs are involved.", "structure": {"pmcid": "PMC6970550", "section_index": 11, "section_title": "GO function and KEGG pathway enrichment analysis of the DEGs", "row_count": 0, "column_count": 0}, "table": {"label": "Table 1", "caption": "GO analyses results of DEGs (top 10 according to adjusted P value).“Count” means how many DEGs are involved.", "rows": [], "row_count": 0, "column_count": 0, "numeric_cells": 0, "header_rows": 0, "locator": "section[11]/table-wrap[1]#table-1"}}, {"id": "asset:table:pmid:31988808:2", "role": "table", "format": "jats-table", "media_type": "application/vnd.jats.table+json", "uri": null, "sha256": "sha256:bc2a2ba74282f6682ec637b947b983096de30777ebb5c8b1741cca9dfe441ec7", "bytes": 270, "access": "open", "state": "verified", "reader": "read_table", "reader_status": "available", "identifier": "section[11]/table-wrap[2]#table-2", "version": null, "note": "Pathway analyses results of DEGs (top 10 according to adjusted P value).“Count” means how many DEGs are involved.", "structure": {"pmcid": "PMC6970550", "section_index": 11, "section_title": "GO function and KEGG pathway enrichment analysis of the DEGs", "row_count": 0, "column_count": 0}, "table": {"label": "Table 2", "caption": "Pathway analyses results of DEGs (top 10 according to adjusted P value).“Count” means how many DEGs are involved.", "rows": [], "row_count": 0, "column_count": 0, "numeric_cells": 0, "header_rows": 0, "locator": "section[11]/table-wrap[2]#table-2"}}, {"id": "asset:table:pmid:31988808:3", "role": "table", "format": "jats-table", "media_type": "application/vnd.jats.table+json", "uri": null, "sha256": "sha256:3ee03a4e34f5523da129a44b6b53d37fd5dd5fb5fee3a9b2ebba37f0e1bc4d27", "bytes": 275, "access": "open", "state": "verified", "reader": "read_table", "reader_status": "available", "identifier": "section[16]/table-wrap[1]#table-3", "version": null, "note": "GO analyses results of hub genes (top 10 according to adjusted P value).“Count” means how many hub genes are involved,", "structure": {"pmcid": "PMC6970550", "section_index": 16, "section_title": "Discussion", "row_count": 0, "column_count": 0}, "table": {"label": "Table 3", "caption": "GO analyses results of hub genes (top 10 according to adjusted P value).“Count” means how many hub genes are involved,", "rows": [], "row_count": 0, "column_count": 0, "numeric_cells": 0, "header_rows": 0, "locator": "section[16]/table-wrap[1]#table-3"}}, {"id": "asset:table:pmid:31988808:4", "role": "table", "format": "jats-table", "media_type": "application/vnd.jats.table+json", "uri": null, "sha256": "sha256:ecf34ada6a8fde18bef65783d987cb3fe7c05b2cbc18877933fcc730c1e405ed", "bytes": 280, "access": "open", "state": "verified", "reader": "read_table", "reader_status": "available", "identifier": "section[16]/table-wrap[2]#table-4", "version": null, "note": "Pathway analyses results of hub genes (top 10 according to adjusted P value).“Count” means how many hub genes are involved.", "structure": {"pmcid": "PMC6970550", "section_index": 16, "section_title": "Discussion", "row_count": 0, "column_count": 0}, "table": {"label": "Table 4", "caption": "Pathway analyses results of hub genes (top 10 according to adjusted P value).“Count” means how many hub genes are involved.", "rows": [], "row_count": 0, "column_count": 0, "numeric_cells": 0, "header_rows": 0, "locator": "section[16]/table-wrap[2]#table-4"}}], "gaps": ["asset:data:geo:GSE12021 not downloaded", "asset:data:geo:GSE12021 licence not read", "asset:data:geo:GSE55235 not downloaded", "asset:data:geo:GSE55235 licence not read", "asset:data:geo:GSE55457 not downloaded", "asset:data:geo:GSE55457 licence not read"]}, "analysis": {"state": "partial", "methods": [{"id": "method-1", "family": "differential_expression", "specific": "limma_voom", "label": "Differential expression analysis", "dsl_operations": ["normalise_counts", "differential_expression_fit", "fdr_bh"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["count_matrix_with_condition"], "evidence": {"score": 0.7, "band": "medium", "quote": "LIMMA"}}, {"id": "method-2", "family": "network_analysis", "specific": "ppi_network", "label": "Network / pathway analysis", "dsl_operations": ["build_graph", "centrality_scores", "pathway_enrich"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["edge_list"], "evidence": {"score": 0.7, "band": "medium", "quote": "protein-protein interaction"}}, {"id": "method-3", "family": "network_analysis", "specific": "pathway_enrichment", "label": "Network / pathway analysis", "dsl_operations": ["pathway_enrich"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["edge_list"], "evidence": {"score": 0.8, "band": "high", "quote": "pathway enrichment"}}, {"id": "method-4", "family": null, "specific": null, "label": "t test", "dsl_operations": [], "missing_operations": [], "expressible": null, "parameters": {}, "data_shapes": [], "evidence": {"score": 1.0, "band": "verbatim", "quote": "The corresponding P value of the gene symbols after t test were used, and adjusted P < 0.05 and —logFC—>2 were used as the selection criteria."}}], "targets": [], "gaps": ["method.reader_named_unmapped"]}, "capability": {"state": "complete", "requirements": [{"capability_id": "capability:computational_analysis", "label": "Computational analysis platform", "assay_families": ["meta_analysis", "machine_learning", "survival_analysis", "statistical_regression", "pca", "molecular_docking", "network_analysis", "differential_expression", "gwas"], "availability_state": "owned", "tier_unlocked": "R1", "candidate": "Automated Assays Laravel application and simulated adapters"}], "gaps": []}, "readiness": {"state": "complete", "tier": "R1", "blocking_axis": "none", "blocked_on": "nothing:replicable_now", "axes": {"method": {"tier": "R1", "expressible": true, "required_operations": ["normalise_counts", "differential_expression_fit", "fdr_bh"], "missing_operations": [], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: limma_voom."}, "capability": {"tier": "R1", "capability_id": "capability:computational_analysis", "availability_state": "owned", "basis": "The computational analysis is owned. The physical bench the paper used is not, so a physical validation cannot be run here."}, "data": {"tier": "R1", "reachable": true, "accessions": ["GSE12021", "GSE55235", "GSE55457"], "resolved_accessions": ["GSE12021", "GSE55235", "GSE55457"], "statements": ["The data is available at NCBI GEO: GSE12021, GSE55235 and GSE55457. The code is available as Supplementary Files."], "basis": "The paper names GSE12021, GSE55235, GSE55457 in Gene Expression Omnibus (GEO) and the accession resolves in the repository. No dataset was downloaded and no licence was read."}}, "can_run": {"dsl_operations": ["normalise_counts", "differential_expression_fit", "fdr_bh"], "assets": ["source:pmid:31988808", "data:geo:GSE12021", "asset:table:pmid:31988808:1", "asset:table:pmid:31988808:2", "asset:table:pmid:31988808:3", "asset:table:pmid:31988808:4"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE12021"]}, "gaps": []}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-10-31T18:17:53Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 3, "usage": {"prompt_tokens": 20342, "completion_tokens": 1423, "cached_tokens": 793}}, "human_critique": null, "stages": [{"id": "extraction", "state": "complete", "reader": "MiniMax-M3.1-Flash", "provider": "minimax"}], "gaps": ["review.independent_reader", "review.machine_critique", "review.human_critique"]}, "execution": {"state": "absent", "status": "not_run", "mode": null, "runs": [], "gaps": ["execution.not_run"]}, "provenance": {"entities": [{"id": "pmid:31988808", "type": "paper", "label": "pmid:31988808"}, {"id": "source:pmid:31988808", "type": "manuscript", "label": "PMID 31988808 full text"}, {"id": "data:geo:GSE12021", "type": "dataset", "label": "GEO GSE12021"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:31988808", "type": "paper", "label": "pmid:31988808"}, {"id": "source:pmid:31988808", "type": "manuscript", "label": "PMID 31988808 full text"}, {"id": "data:geo:GSE12021", "type": "dataset", "label": "GEO GSE12021"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-31988808", "type": "extraction", "used": ["source:pmid:31988808"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE12021", "type": "accession_resolution", "used": ["data:geo:GSE12021"], "agent": "agent:corpus-assay-engine"}], "agents": [{"id": "agent:corpus-assay-engine", "type": "software", "label": "Corpus assay engine"}, {"id": "agent:minimax-reader", "type": "software", "label": "MiniMax literature reader"}], "source_run_id": "reading-sci-2026-11-01", "dsl_catalogue_fingerprint": "343fb30756ba6ec3", "capability_table_fingerprint": "2cb3e42f047d2ecf", "snapshot_id": "ba34b49490d9bfaf"}, "lineage": {"parent_id": "pmid:31988808", "parent_hash": "sha256:31d612a1d065ad575ffaf5e7bf7714440b93f962e49ace435e1527b055150a6c", "fork_of": null, "generation": 6}, "gaps": ["asset:data:geo:GSE12021 licence not read", "asset:data:geo:GSE12021 not downloaded", "asset:data:geo:GSE55235 licence not read", "asset:data:geo:GSE55235 not downloaded", "asset:data:geo:GSE55457 licence not read", "asset:data:geo:GSE55457 not downloaded", "execution.not_run", "method.reader_named_unmapped", "review.human_critique", "review.independent_reader", "review.machine_critique"], "assays": [{"role": "analysis", "definition_id": "differential_expression", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": ["normalise_counts", "differential_expression_fit", "fdr_bh"], "support_state": "supported", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 0.7, "evidence": ["LIMMA"], "order_basis": "analysis_only", "id": "assay:1"}]}