Identification of Key Genes and Immune Infiltrate in Nonalcoholic Steatohepatitis: A Bioinformatic Analysis.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:34552988", "kind": "science-record", "revision": 8, "content_hash": "sha256:765f38675577012d6227c7da1674dbe5d562ad317fa2f1aea84db9887598e228", "created_at": "2026-11-01T15:41:53Z", "updated_at": "2026-11-01T22:47:34Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Identification of Key Genes and Immune Infiltrate in Nonalcoholic Steatohepatitis: A Bioinformatic Analysis.", "doi": "10.1155/2021/7561645", "pmid": "34552988", "authors": ["Zhen-Yu Jiang", "Yi Zhou", "Lu Zhou", "Shao-Wei Li", "Bang-Mao Wang"], "journal": "BioMed research international", "year": 2021, "abstract": "BACKGROUND: Nonalcoholic steatohepatitis (NASH) can progress to cirrhosis and hepatic carcinoma and is closely associated with changes in the neurological environment. The discovery of new biomarkers would aid in the treatment of NASH.\nMETHODS: Data GSE89632 were downloaded from the Gene Expression Omnibus (GEO) database, and R package \"limma\" was used to identify differentially expressed genes (DEGs) for NASH vs. normal tissues. The STRING database was used to construct a protein-protein interaction (PPI) network, and the Cytoscape software program (Version 3.80) was used to visualize the PPI network and identify key genes. The immune infiltration of NASH was determined using the R package \"CIBERSORT\".\nRESULTS: We screened 41 DEGs. GO and KEGG enrichment analyses of the DEGs revealed the enrichment of pathways related to NAFLD steatosis and inflammation. A PPI network analysis was also performed on the DEGs, and seven genes (MYC, CXCL8, FOS, SOCS1, SOCS3, IL6, and PTGS2) were identified as hub genes. An immune infiltration assessment revealed that macrophages M2, memory resting CD4+ T cells, and γΔ T cells play important roles in the immune microenvironment of NASH, which may be mediated by the seven identified hub genes.", "claims": [{"id": "claim-method", "statement": "We downloaded the gene expression matrix of GSE89632 and the annotation file of GPL14951. First, we replace the ID of the probe with the gene symbol, and if they have the same probe, the max one will be retained. According to the aim of our study, we retained expression profiling data for only 24 HCs and 19 NASH patients. As the data had already been processed in a “quantile normalized” fashion, we proceeded directly to the next analysis. The “limma” package was used to identify differentially expressed genes (DEGs) between NASH patients and HCs; genes were considered DEGs if they had |logFC | >2 and P value <0.05.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:34552988", "locator": "sections[4] \"2.2. Data Process and DEGs of NASH Identified\" paragraph 2", "quote": "We downloaded the gene expression matrix of GSE89632 and the annotation file of GPL14951. First, we replace the ID of the probe with the gene symbol, and if they have the same probe, the max one will be retained. According to the aim of our study, we retained expression profiling data for only 24 HCs and 19 NASH patients. As the data had already been processed in a “quantile normalized” fashion, we proceeded directly to the next analysis. The “limma” package was used to identify differentially expressed genes (DEGs) between NASH patients and HCs; genes were considered DEGs if they had |logFC | >2 and P value <0.05."}], "state": "complete"}, {"id": "claim-data-access", "statement": "The data used to support the findings of this study are included within the article.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:34552988", "locator": "sections[17] \"Data Availability\" paragraph 2", "quote": "The data used to support the findings of this study are included within the article."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "partial", "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:34552988"], "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": "blocked", "blocked_on": "data:accession_absent", "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": [], "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": [], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": ["step:obtain_source_data blocked on data:accession_absent"]}, "data": {"state": "partial", "assets": [{"id": "source:pmid:34552988", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1155/2021/7561645", "sha256": "sha256:6abf8015faaa223e3c9f35b65156fd379dfd665ee3e32d191ca1266361ba9798", "bytes": 41609, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1155/2021/7561645", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC8452393", "source": "europe_pmc_fulltext", "coverage": "whole_body", "text_hash_scope": "rendered_section_body"}}, {"id": "asset:table:pmid:34552988:1", "role": "table", "format": "jats-table", "media_type": "application/vnd.jats.table+json", "uri": null, "sha256": "sha256:653cd4c92360a042ef0ae71b65641364f561d5d952f1de212e077346cd1c8771", "bytes": 421, "access": "open", "state": "verified", "reader": "read_table", "reader_status": "available", "identifier": "section[11]/table-wrap[1]#tab1", "version": null, "note": "Top 10 DEGs in the network ranked by the MCC method.", "structure": {"pmcid": "PMC8452393", "section_index": 11, "section_title": "3.3. Identification of the DEG PPI Network and Hub Genes", "row_count": 11, "column_count": 3}, "table": {"label": "Table 1", "caption": "Top 10 DEGs in the network ranked by the MCC method.", "rows": [["Rank", "Name", "Score"], ["1", "IL6", "1101"], ["2", "FOS", "1082"], ["3", "MYC", "926"], ["4", "CXCL8", "774"], ["5", "PTGS2", "769"], ["6", "SOCS3", "744"], ["7", "SOCS1", "726"], ["8", "FOSB", "294"], ["9", "FOSL1", "264"], ["9", "JUNB", "264"]], "row_count": 11, "column_count": 3, "numeric_cells": 20, "header_rows": 1, "locator": "section[11]/table-wrap[1]#tab1"}}], "gaps": []}, "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": "We downloaded the gene expression matrix of GSE89632 and the annotation file of GPL14951. First, we replace the ID of the probe with the gene symbol, and if they have the same probe, the max one will be retained. According to the aim of our study, we retained expression profiling data for only 24 HCs and 19 NASH patients. As the data had already been processed in a “quantile normalized” fashion, we proceeded directly to the next analysis. The “limma” package was used to identify differentially expressed genes (DEGs) between NASH patients and HCs; genes were considered DEGs if they had |logFC | >2 and P value <0.05."}}, {"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": "GO enrichment"}}], "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": "blocked", "tier": "R4", "blocking_axis": "data", "blocked_on": "data:accession_absent", "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": "R4", "reachable": false, "accessions": [], "resolved_accessions": [], "statements": ["The data used to support the findings of this study are included within the article."], "basis": "The paper names no accession to its source data, so the data has to come from the authors."}}, "can_run": {"dsl_operations": ["normalise_counts", "differential_expression_fit", "fdr_bh"], "assets": ["source:pmid:34552988", "asset:table:pmid:34552988:1"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:None"]}, "gaps": ["data:accession_absent"]}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-10-31T19:49:45Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 3, "usage": {"prompt_tokens": 33237, "completion_tokens": 1313, "cached_tokens": 21760}}, "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:34552988", "type": "paper", "label": "pmid:34552988"}, {"id": "source:pmid:34552988", "type": "manuscript", "label": "PMID 34552988 full text"}, {"id": "data:geo:None", "type": "dataset", "label": "GEO None"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:34552988", "type": "paper", "label": "pmid:34552988"}, {"id": "source:pmid:34552988", "type": "manuscript", "label": "PMID 34552988 full text"}, {"id": "data:geo:None", "type": "dataset", "label": "GEO None"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-34552988", "type": "extraction", "used": ["source:pmid:34552988"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-None", "type": "accession_resolution", "used": ["data:geo:None"], "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:34552988", "parent_hash": "sha256:7a84abf10a952e14fba15150f27357d8ae605f65dbf6053bb662d301ade36891", "fork_of": null, "generation": 7}, "gaps": ["execution.not_run", "review.human_critique", "review.independent_reader", "review.machine_critique", "step:obtain_source_data blocked on data:accession_absent"], "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": 1.0, "evidence": ["We downloaded the gene expression matrix of GSE89632 and the annotation file of GPL14951. First, we replace the ID of the probe with the gene symbol, and if they have the same probe, the max one will be retained. According to the aim of our study, we retained expression profiling data for only 24 HCs and 19 NASH patients. As the data had already been processed in a “quantile normalized” fashion, we proceeded directly to the next analysis. The “limma” package was used to identify differentially expressed genes (DEGs) between NASH patients and HCs; genes were considered DEGs if they had |logFC | >2 and P value <0.05."], "order_basis": "analysis_only", "id": "assay:1"}]}