IRF8 and its related molecules as potential diagnostic biomarkers or therapeutic candidates and immune cell infiltration characteristics in steroid-induced osteonecrosis of the femoral head.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:36627660", "kind": "science-record", "revision": 4, "content_hash": "sha256:5e444a2fce676728f7f3b9eec0a9c231bab477e315abaf7dfdba446af3858b28", "created_at": "2026-11-01T19:01:11Z", "updated_at": "2026-11-01T22:51:01Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "IRF8 and its related molecules as potential diagnostic biomarkers or therapeutic candidates and immune cell infiltration characteristics in steroid-induced osteonecrosis of the femoral head.", "doi": "10.1186/s13018-022-03381-1", "pmid": "36627660", "authors": ["Xue-Zhen Liang", "Xiao-Chen Liu", "Song Li", "Ming-Tao Wen", "Yan-Rong Chen", "Di Luo", "Bo Xu", "Nian-Hu Li", "Gang Li"], "journal": "Journal of orthopaedic surgery and research", "year": 2023, "abstract": "PURPOSE: Steroid-induced osteonecrosis of the femoral head (SONFH) was a refractory orthopedic hip joint disease in the young and middle-aged people, but the pathogenesis of SONFH remained unclear. We aimed to identify the potential genes and screen potential therapeutic compounds for SONFH.\nMETHODS: The microarray was obtained for blood tissue from the GEO database, and then it identifies differentially expressed genes (DEGs). The DEGs were analyzed to obtain the differences in immune cell infiltration. The gene functional enrichment analysis of SONFH was analyzed. The PPI of DEGs was identified through the STRING database, and the cluster modules and hub genes were ascertained using MCODE and CytoHubba, and the ROC curve of hub genes was analyzed, and the tissues distribution of hub genes was understood by the HPA, Bgee and BioGPS databases. The hub genes and target miRNAs and corresponding upstream lncRNAs were predicted by TargetScan, miRDB and ENCORI database. Subsequently, we used CMap, DGIdb and L1000FWD databases to identify several potential therapeutic molecular compounds for SONFH. Finally, the AutoDockTools Vina, PyMOL and Discovery Studio were employed for molecular docking analyses between compounds and hub genes.\nRESULTS: The microarray dataset GSE123568 was obtained related to SONFH. There were 372 DEGs including 197 upregulated genes and 175 downregulated genes by adjusted P value < 0.01 and |log2FC|> 1. Several significant GSEA enrichment analysis and biological processes and KEGG pathway associated with SONFH were identified, which were significantly related to cytoskeleton organization, nucleobase-containing compound catabolic process, NOD-like receptor signaling pathway, MAPK signaling pathway, FoxO signaling pathway, neutrophil-mediated immunity, neutrophil degranulation and neutrophil activation involved in immune response. Activated T cells CD4 memory, B cells naïve, B cells memory, T cells CD8 and T cells gamma delta might be involved in the occurrence and development of SONFH. Three cluster modules were identified in the PPI network, and eleven hub genes including FPR2, LILRB2, MNDA, CCR1, IRF8, TYROBP, TLR1, HCK, TLR8, TLR2 and CCR2 were identified by Cytohubba, which were differed in bone marrow, adipose tissue and blood, and which had good diagnostic performance in SONFH. We identified IRF8 and 10 target miRNAs that was utilized including Targetsan, miRDB and ENCORI databases and 8 corresponding upstream lncRNAs that was revealed by ENCORI database. IRF8 was detected with consistent expression by qRT-PCR. Based on the CMap, DGIdb and L1000FWD databases, the 11 small molecular compounds that were most strongly therapeutic correlated with SONFH were estradiol, genistein, domperidone, lovastatin, myricetin, fenbufen, rosiglitazone, sirolimus, phenformin, vorinostat and vinblastine. All of 11 small molecules had good binding affinity with the IRF8 in molecular docking.\nCONCLUSION: The occurrence of SONFH was associated with a \"multi-target\" and \"multi-pathway\" pattern, especially related to immunity, and IRF8 and its noncoding RNA were closely related to the development of SONFH. The CMap, DGIdb and L1000FWD databases could be effectively used in a systematic manner to predict potential drugs for the prevention and treatment of SONFH. However, additional clinical and experimental research is warranted.", "claims": [{"id": "claim-method", "statement": "We used “limma” R package to analysis the DEGs.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:36627660", "locator": "sections[3] \"Microarray data quality assessment and identification of DEGs\" paragraph 4", "quote": "We used “limma” R package to analysis the DEGs."}], "state": "complete"}, {"id": "claim-data-access", "statement": "Publicly available datasets were analyzed in this study. These data can be found at the following URL: GSE123568 dataset (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE123568).", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:36627660", "locator": "sections[33] \"Availability of data and materials\" paragraph 2", "quote": "Publicly available datasets were analyzed in this study. These data can be found at the following URL: GSE123568 dataset (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE123568)."}], "state": "complete"}, {"id": "claim-target-result", "statement": "There were 372 DEGs including 197 upregulated genes and 175 downregulated genes by adjusted P value < 0.01 and |log2FC|> 1.", "kind": "quantitative", "target": {"value": 0.01, "unit": "p-value", "tolerance": null, "source": "There were 372 DEGs including 197 upregulated genes and 175 downregulated genes by adjusted P value < 0.01 and |log2FC|> 1."}, "evidence": [{"asset": "source:pmid:36627660", "locator": "sections[0] \"Abstract\" paragraph 3", "quote": "There were 372 DEGs including 197 upregulated genes and 175 downregulated genes by adjusted P value < 0.01 and |log2FC|> 1."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "complete", "objective": "There were 372 DEGs including 197 upregulated genes and 175 downregulated genes by adjusted P value < 0.01 and |log2FC|> 1.", "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:36627660"], "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:GSE123568"], "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:GSE123568"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": []}, "data": {"state": "partial", "assets": [{"id": "source:pmid:36627660", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1186/s13018-022-03381-1", "sha256": "sha256:d825273c54afbcd8bd54375b46eeb7725187da5643b6a457cfef8964fe40e33c", "bytes": 74610, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1186/s13018-022-03381-1", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC9832881", "source": "europe_pmc_fulltext", "coverage": "whole_body", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE123568", "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=GSE123568", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE123568", "version": null, "note": "Publicly available datasets were analyzed in this study. These data can be found at the following URL: GSE123568 dataset (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE123568).", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "Publicly available datasets were analyzed in this study. These data can be found at the following URL: GSE123568 dataset (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE123568).", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM3507251", "GSM3507252", "GSM3507253", "GSM3507254", "GSM3507255", "GSM3507256", "GSM3507257", "GSM3507258", "GSM3507259", "GSM3507260", "GSM3507261", "GSM3507262", "GSM3507263", "GSM3507264", "GSM3507265", "GSM3507266", "GSM3507267", "GSM3507268", "GSM3507269", "GSM3507270", "GSM3507271", "GSM3507272", "GSM3507273", "GSM3507274", "GSM3507275", "GSM3507276", "GSM3507277", "GSM3507278", "GSM3507279", "GSM3507280", "GSM3507281", "GSM3507282", "GSM3507283", "GSM3507284", "GSM3507285", "GSM3507286", "GSM3507287", "GSM3507288", "GSM3507289", "GSM3507290"], "sample_count": 40, "group_sizes": null, "group_sizes_basis": null, "platform": "15207", "taxon": "Homo sapiens", "title": "Identification of Potential Biomarkers for Improving the Precision of Early Detection of Steroid-induced Osteonecrosis of the Femoral Head", "resolution_checked_at": "2026-11-01T15:39:54Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}], "gaps": ["asset:data:geo:GSE123568 not downloaded", "asset:data:geo:GSE123568 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": 1.0, "band": "verbatim", "quote": "We used “limma” R package to analysis the DEGs."}}, {"id": "method-2", "family": "principal_component_analysis", "specific": "pca", "label": "Principal component analysis", "dsl_operations": ["pca_fit"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["wide_numeric_matrix"], "evidence": {"score": 0.8, "band": "high", "quote": "principal components analysis"}}, {"id": "method-3", "family": "clustering", "specific": "hierarchical", "label": "Clustering", "dsl_operations": ["hierarchical_cluster"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["feature_matrix"], "evidence": {"score": 0.8, "band": "high", "quote": "Hierarchical clustering"}}, {"id": "method-4", "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": "PPI"}}, {"id": "method-5", "family": "diagnostic_evaluation", "specific": "roc", "label": "Diagnostic test evaluation", "dsl_operations": ["confusion_counts", "roc_curve", "sensitivity_specificity", "solve_confusion_matrix", "classification_metrics", "published_precision_envelope"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["reference_standard_and_test"], "evidence": {"score": 0.7, "band": "medium", "quote": "ROC curve"}}], "targets": [{"id": "target-published-estimate", "claim_id": "claim-target-result", "operation": "build_graph", "expected_value": 0.01, "unit": "p-value", "tolerance": null, "state": "complete", "note": "The threshold the paper applied to its own result, as the paper states it. A run that reproduces the effect is judged against this threshold."}], "gaps": ["method.reader_differs_from_extractor"]}, "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": ["GSE123568"], "resolved_accessions": ["GSE123568"], "statements": ["Publicly available datasets were analyzed in this study. These data can be found at the following URL: GSE123568 dataset (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE123568)."], "basis": "The paper names GSE123568 in National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database 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:36627660", "data:geo:GSE123568"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE123568"]}, "gaps": []}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-11-01T18:22:46Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 20207, "completion_tokens": 1115, "cached_tokens": 264}}, "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:36627660", "type": "paper", "label": "pmid:36627660"}, {"id": "source:pmid:36627660", "type": "manuscript", "label": "PMID 36627660 full text"}, {"id": "data:geo:GSE123568", "type": "dataset", "label": "GEO GSE123568"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:36627660", "type": "paper", "label": "pmid:36627660"}, {"id": "source:pmid:36627660", "type": "manuscript", "label": "PMID 36627660 full text"}, {"id": "data:geo:GSE123568", "type": "dataset", "label": "GEO GSE123568"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-36627660", "type": "extraction", "used": ["source:pmid:36627660"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE123568", "type": "accession_resolution", "used": ["data:geo:GSE123568"], "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:36627660", "parent_hash": "sha256:b1f83e9c72cdc3a0cc5e7be65019011e1bb2b34e97836d67981f1bdd88e9f2b4", "fork_of": null, "generation": 3}, "gaps": ["asset:data:geo:GSE123568 licence not read", "asset:data:geo:GSE123568 not downloaded", "execution.not_run", "method.reader_differs_from_extractor", "review.human_critique", "review.independent_reader", "review.machine_critique"], "assays": [{"role": "data_generation", "definition_id": "microarray", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.99, "evidence": ["herapeutic compounds for SONFH. METHODS: The microarray was obtained for blood tissue from the GEO d", "etween compounds and hub genes. RESULTS: The microarray dataset GSE123568 was obtained related to SO", "ld for significant statistical significance. Microarray data processing and identification of DEGs T", "ta processing and identification of DEGs The microarray dataset GSE123568 about identification of po"], "order_basis": "paper_text_order:452", "id": "assay:1"}, {"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": 1.0, "evidence": ["We used “limma” R package to analysis the DEGs."], "order_basis": "consumes_earlier_output", "id": "assay:2"}]}