Comprehensive Bioinformatics Analysis Reveals Hub Genes and Inflammation State of Rheumatoid Arthritis.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:32802866", "kind": "science-record", "revision": 8, "content_hash": "sha256:08652bbdc37d4e2b8e41a4ba48a42e4b73f7d47742a956ae0c5d0fd464484e61", "created_at": "2026-11-01T15:41:08Z", "updated_at": "2026-11-01T22:45:13Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Comprehensive Bioinformatics Analysis Reveals Hub Genes and Inflammation State of Rheumatoid Arthritis.", "doi": "10.1155/2020/6943103", "pmid": "32802866", "authors": ["Conglin Ren", "Mingshuang Li", "Weibin Du", "Jianlan Lü", "Yang Zheng", "Haipeng Xu", "Renfu Quan"], "journal": "BioMed research international", "year": 2020, "abstract": "Rheumatoid arthritis (RA) is an autoimmune disease characterized by erosive arthritis, which has not been thoroughly cured yet, and standardized treatment is helpful for alleviating clinical symptoms. Here, various bioinformatics analysis tools were comprehensively utilized, aiming to identify critical biomarkers and possible pathogenesis of RA. Three gene expression datasets profiled by microarray were obtained from GEO database. Dataset GSE55235 and GSE55457 were merged for subsequent analyses. We identified differentially expressed genes (DEGs) in RStudio with limma package, performing functional enrichment analysis based on GSEA software and clusterProfiler package. Next, protein-protein interaction (PPI) network was set up through STRING database and Cytoscape. Moreover, CIBERSORT website was used to assess the inflammatory state of RA. Finally, we validated the candidate hub genes with dataset GSE77298. As a result, we identified 106 DEGs (72 upregulated and 34 downregulated genes). Through GO, KEGG, and GSEA analysis, we found that DEGs were mainly involved in immune response and inflammatory signaling pathway. With the help of Cytoscape software and MCODE plug-in, the most prominent subnetwork was screened out, containing 14 genes and 45 edges. For ROC curve analysis, eight genes with AUC >0.80 were considered as hub genes of RA. In conclusion, compared with healthy controls, the DEGs and their closely related biological functions were analyzed, and we held that chemokines and immune cells infiltration promote the progression of rheumatoid arthritis. Targeting the eight biomarkers we identified may be useful for the diagnosis and treatment of rheumatoid arthritis.", "claims": [{"id": "claim-method", "statement": "Differentially expressed genes (DEGs) were screened out by limma package [16] based on the comparison of expression values between HC samples and RA samples.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:32802866", "locator": "sections[4] \"2.2. Identification of DEGs\" paragraph 2", "quote": "Differentially expressed genes (DEGs) were screened out by limma package [16] based on the comparison of expression values between HC samples and RA samples."}], "state": "complete"}, {"id": "claim-data-access", "statement": "The data that support the findings of this study are openly available.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:32802866", "locator": "sections[19] \"Data Availability\" paragraph 2", "quote": "The data that support the findings of this study are openly available."}], "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:32802866"], "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:GSE55235"], "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:GSE55235"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": []}, "data": {"state": "partial", "assets": [{"id": "source:pmid:32802866", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1155/2020/6943103", "sha256": "sha256:4abcb4aebb758c074b0b02b13b23d54d7496398005e309e79aaa31e8763b7c46", "bytes": 34267, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1155/2020/6943103", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC7424395", "source": "europe_pmc_fulltext", "coverage": "whole_body", "text_hash_scope": "rendered_section_body"}}, {"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 that support the findings of this study are openly available.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The data that support the findings of this study are openly available.", "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 that support the findings of this study are openly available.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The data that support the findings of this study are openly available.", "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": "data:geo:GSE77298", "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=GSE77298", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE77298", "version": null, "note": "The data that support the findings of this study are openly available.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The data that support the findings of this study are openly available.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM2048265", "GSM2048266", "GSM2048267", "GSM2048268", "GSM2048269", "GSM2048270", "GSM2048271", "GSM2048272", "GSM2048273", "GSM2048274", "GSM2048275", "GSM2048276", "GSM2048277", "GSM2048278", "GSM2048279", "GSM2048280", "GSM2048281", "GSM2048282", "GSM2048283", "GSM2048284", "GSM2048285", "GSM2048286", "GSM2048287"], "sample_count": 23, "group_sizes": null, "group_sizes_basis": null, "platform": "570", "taxon": "Homo sapiens", "title": "Synovial biopsies of rheumatoid arthritis and healthy controls", "resolution_checked_at": "2026-11-01T15:39:55Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}, {"id": "asset:table:pmid:32802866:1", "role": "table", "format": "jats-table", "media_type": "application/vnd.jats.table+json", "uri": null, "sha256": "sha256:e44ad0186f8a74951e40cbffc5ff1a29e3bbe0add71fa18cbec1881d1e9bdea4", "bytes": 1880, "access": "open", "state": "verified", "reader": "read_table", "reader_status": "available", "identifier": "section[11]/table-wrap[1]#tab1", "version": null, "note": "GO analysis results of DEGs (top 8 terms of each category were listed).", "structure": {"pmcid": "PMC7424395", "section_index": 11, "section_title": "3.2. GO and KEGG Enrichment Analysis", "row_count": 25, "column_count": 5}, "table": {"label": "Table 1", "caption": "GO analysis results of DEGs (top 8 terms of each category were listed).", "rows": [["Ontology", "ID", "Description", "Adj. p value", "Count"], ["BP", "GO:0050900", "Leukocyte migration", "2.33E-08", "18"], ["BP", "GO:0050851", "Antigen receptor-mediated signaling pathway", "2.33E-08", "15"], ["BP", "GO:0002449", "Lymphocyte mediated immunity", "6.91E-08", "15"], ["BP", "GO:0051249", "Regulation of lymphocyte activation", "6.91E-08", "17"], ["BP", "GO:0050853", "B cell receptor signaling pathway", "2.54E-07", "10"], ["BP", "GO:0002429", "Immune response-activating cell surface receptor signaling pathway", "3.05E-07", "16"], ["BP", "GO:0030098", "Lymphocyte differentiation", "4.36E-07", "14"], ["BP", "GO:0006959", "Humoral immune response", "4.36E-07", "14"], ["CC", "GO:0009897", "External side of plasma membrane", "1.28E-08", "17"], ["CC", "GO:0042571", "Immunoglobulin complex, circulating", "7.02E-06", "7"], ["CC", "GO:0019814", "Immunoglobulin complex", "8.86E-06", "9"], ["CC", "GO:0042101", "T cell receptor complex", "0.008265", "5"], ["CC", "GO:0001772", "Immunological synapse", "0.012398", "3"], ["CC", "GO:0072562", "Blood microparticle", "0.012678", "5"], ["CC", "GO:0042613", "MHC class II protein complex", "0.030215", "2"], ["CC", "GO:0098802", "Plasma membrane receptor complex", "0.033455", "6"], ["MF", "GO:0034987", "Immunoglobulin receptor binding", "8.86E-06", "7"], ["MF", "GO:0008009", "Chemokine activity", "1.34E-05", "6"], ["MF", "GO:0045236", "CXCR chemokine receptor binding", "1.51E-05", "4"], ["MF", "GO:0042379", "Chemokine receptor binding", "5.51E-05", "6"], ["MF", "GO:0003823", "Antigen binding", "7.28E-05", "8"], ["MF", "GO:0005125", "Cytokine activity", "7.80E-05", "9"], ["MF", "GO:0048018", "Receptor ligand activity", "0.000201", "12"], ["MF", "GO:0001664", "G protein-coupled receptor binding", "0.002229", "8"]], "row_count": 25, "column_count": 5, "numeric_cells": 72, "header_rows": 1, "locator": "section[11]/table-wrap[1]#tab1"}}, {"id": "asset:table:pmid:32802866:2", "role": "table", "format": "jats-table", "media_type": "application/vnd.jats.table+json", "uri": null, "sha256": "sha256:42a873cd58917c922cbbceb0a1b5af3e8b8ae3268ab4406f9dd49a0a8366122d", "bytes": 842, "access": "open", "state": "verified", "reader": "read_table", "reader_status": "available", "identifier": "section[11]/table-wrap[2]#tab2", "version": null, "note": "Top 10 pathways of KEGG analysis.", "structure": {"pmcid": "PMC7424395", "section_index": 11, "section_title": "3.2. GO and KEGG Enrichment Analysis", "row_count": 11, "column_count": 4}, "table": {"label": "Table 2", "caption": "Top 10 pathways of KEGG analysis.", "rows": [["ID", "Description", "Adj. p value", "Count"], ["hsa04060", "Cytokine-cytokine receptor interaction", "1.44E-05", "13"], ["hsa04061", "Viral protein interaction with cytokine and cytokine receptor", "0.000459", "7"], ["hsa04062", "Chemokine signaling pathway", "0.002564", "8"], ["hsa04659", "Th17 cell differentiation", "0.003911", "6"], ["hsa03320", "PPAR signaling pathway", "0.005839", "5"], ["hsa04658", "Th1 and Th2 cell differentiation", "0.011245", "5"], ["hsa04657", "IL-17 signaling pathway", "0.011245", "5"], ["hsa04640", "Hematopoietic cell lineage", "0.01248", "5"], ["hsa04620", "Toll-like receptor signaling pathway", "0.013885", "5"], ["hsa05020", "Prion diseases", "0.026379", "3"]], "row_count": 11, "column_count": 4, "numeric_cells": 32, "header_rows": 1, "locator": "section[11]/table-wrap[2]#tab2"}}], "gaps": ["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", "asset:data:geo:GSE77298 not downloaded", "asset:data:geo:GSE77298 licence not read"]}, "analysis": {"state": "complete", "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": "Differentially expressed genes (DEGs) were screened out by limma package [16] based on the comparison of expression values between HC samples and RA samples."}}, {"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": "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": [], "gaps": []}, "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": ["GSE55235", "GSE55457", "GSE77298"], "resolved_accessions": ["GSE55235", "GSE55457", "GSE77298"], "statements": ["The data that support the findings of this study are openly available."], "basis": "The paper names GSE55235, GSE55457, GSE77298 in 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:32802866", "data:geo:GSE55235", "asset:table:pmid:32802866:1", "asset:table:pmid:32802866:2"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE55235"]}, "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:22:00Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 3, "usage": {"prompt_tokens": 28375, "completion_tokens": 1498, "cached_tokens": 9736}}, "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:32802866", "type": "paper", "label": "pmid:32802866"}, {"id": "source:pmid:32802866", "type": "manuscript", "label": "PMID 32802866 full text"}, {"id": "data:geo:GSE55235", "type": "dataset", "label": "GEO GSE55235"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:32802866", "type": "paper", "label": "pmid:32802866"}, {"id": "source:pmid:32802866", "type": "manuscript", "label": "PMID 32802866 full text"}, {"id": "data:geo:GSE55235", "type": "dataset", "label": "GEO GSE55235"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-32802866", "type": "extraction", "used": ["source:pmid:32802866"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE55235", "type": "accession_resolution", "used": ["data:geo:GSE55235"], "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:32802866", "parent_hash": "sha256:ed7d9dfe132ee83a93e509ad16afc61e9d6b37e4054cce2cb47ab44ea7996c01", "fork_of": null, "generation": 7}, "gaps": ["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", "asset:data:geo:GSE77298 licence not read", "asset:data:geo:GSE77298 not downloaded", "execution.not_run", "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.93, "evidence": [". Three gene expression datasets profiled by microarray were obtained from GEO database. Dataset GSE"], "order_basis": "paper_text_order:450", "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": ["Differentially expressed genes (DEGs) were screened out by limma package [16] based on the comparison of expression values between HC samples and RA samples."], "order_basis": "consumes_earlier_output", "id": "assay:2"}]}