Single-Cell Transcriptome Analysis Reveals Different Immune Signatures in HPV- and HPV + Driven Human Head and Neck Squamous Cell Carcinoma.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:36157879", "kind": "science-record", "revision": 4, "content_hash": "sha256:b73e26a888ea77142653870e9c2f2a1f2d7d0fb98a829fa0c77d2545bfd3f01f", "created_at": "2026-11-01T18:15:26Z", "updated_at": "2026-11-01T22:07:39Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Single-Cell Transcriptome Analysis Reveals Different Immune Signatures in HPV- and HPV + Driven Human Head and Neck Squamous Cell Carcinoma.", "doi": "10.1155/2022/2079389", "pmid": "36157879", "authors": ["Simin Li", "Yang Wang", "Rui Sun", "Debora Franceschi", "Hongying Pan", "Chenxuan Wei", "Anthony Chukwunonso Ogbuehi", "Bernd Lethaus", "Vuk Savkovic", "Sebastian Gaus", "Rüdiger Zimmerer", "Dirk Ziebolz", "Gerhard Schmalz", "Xiao Jiang"], "journal": "Journal of immunology research", "year": 2022, "abstract": "Background: Head and neck squamous cell carcinoma (HNSCC) is a significant health problem and related to poor long-term outcomes, indicating more research to be done to deeply understand the underlying pathways.\nObjective: This current study aimed in the assessment of the viral- (especially human papilloma virus [HPV]) and carcinogen-driven head and neck squamous cell carcinoma (HNSCC) microenvironment based on single-cell sequencing analysis.\nMethods: Data were downloaded from GEO database (GSE139324), including 131224 cells from 18 HP- HNSCC patients and 8 HPV+ HNSCC patients. Following data normalization, all highly variable genes in single cells were identified, and batch correction was applied. Differentially expressed genes were identified using Wilcoxon rank sum test. A gene enrichment analysis was performed in each cell cluster using KEGG analysis. Single-cell pseudotime trajectories were constructed with MONOCLE (version 2.6.4). Cell-cell interactions were analyzed with CellChat R package. Additionally, cell-cell communication patterns in key signal pathways were compared in different tissue groups. A hidden Markov model (HMM) was used to predict gene expression states (on or off) throughout pseudotime. Five-year overall survival outcomes were compared in both HPV+ and HPV- subsets.\nResults: 20,978 high-quality individual cells passed quality control. RNA-seq data were used from 522 HNSCC primary tumor samples. 1,137 differentially expressed genes between HPV+ and HPV- HNSCC patients were investigated. 96 differentially expressed genes were associated with overall survival and highly enriched in B cell associated biological process. Cell composition differed between types of samples. MHC-I, MHC-II, and MIF signaling pathways were found to be most relevant. Within these pathways, some cells were either signal receiver or signal sender, depending on sample type, respectively. Six genes were obtained, AREG and TGFBI (upregulation), CD27, CXCR3, MS4A1, and CD19 (downregulation), whose expression and HPV types were highly associated with worse overall survival. AREG and TGFBI were pDC marker genes, CXCR3 and CD27 were significantly expressed in T cell-related cells, while MS4A1 and CD19 were mainly expressed in B naïve cells.\nConclusions: This study revealed dynamic changes in cell percentage and heterogeneity of cell subtypes of HNSCC. AREG, TGFBI, CD27, CXCR3, MS4A1, and CD19 were associated with worse overall survival in HPV-related HNSCC. Especially B-cell related pathways were revealed as particularly relevant in HPV-related HNSCC. These findings are a basis for the development of biomarkers and therapeutic targets in respective patients.", "claims": [{"id": "claim-method", "statement": "Differentially expressed genes were identified using Wilcoxon rank sum test.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:36157879", "locator": "sections[0] \"Abstract\" paragraph 3", "quote": "Differentially expressed genes were identified using Wilcoxon rank sum test."}], "state": "complete"}, {"id": "claim-data-access", "statement": "The data analyzed during the current study are available in the TCGA database with the accession numbers TCGA-HNSCC, and GEO repository (GSE139324 dataset).", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:36157879", "locator": "sections[19] \"Data Availability\" paragraph 2", "quote": "The data analyzed during the current study are available in the TCGA database with the accession numbers TCGA-HNSCC, and GEO repository (GSE139324 dataset)."}], "state": "complete"}, {"id": "claim-target-result", "statement": "Six genes were obtained, AREG and TGFBI (upregulation), CD27, CXCR3, MS4A1, and CD19 (downregulation), whose expression and HPV types were highly associated with worse overall survival.", "kind": "qualitative", "target": null, "evidence": [{"asset": "source:pmid:36157879", "locator": "sections[0] \"Abstract\" paragraph 4", "quote": "Six genes were obtained, AREG and TGFBI (upregulation), CD27, CXCR3, MS4A1, and CD19 (downregulation), whose expression and HPV types were highly associated with worse overall survival."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "partial", "objective": "Six genes were obtained, AREG and TGFBI (upregulation), CD27, CXCR3, MS4A1, and CD19 (downregulation), whose expression and HPV types were highly associated with worse overall survival.", "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:36157879"], "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:GSE139324"], "outputs": ["the source data with its repository provenance"], "instruments": [], "reagents": []}, {"id": "run_analysis", "order": 3, "name": "Run the published analysis", "kind": "analysis", "state": "blocked", "blocked_on": "method:group_comparison_wilcoxon", "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:GSE139324"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": ["step:run_analysis blocked on method:group_comparison_wilcoxon"]}, "data": {"state": "partial", "assets": [{"id": "source:pmid:36157879", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1155/2022/2079389", "sha256": "sha256:acd609f01a64760dbecbc67e909e940cd20e5261eb043a35a49108fade1e997a", "bytes": 67218, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1155/2022/2079389", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC9507777", "source": "europe_pmc_fulltext", "coverage": "whole_body", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE139324", "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=GSE139324", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE139324", "version": null, "note": "The data analyzed during the current study are available in the TCGA database with the accession numbers TCGA-HNSCC, and GEO repository (GSE139324 dataset).", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The data analyzed during the current study are available in the TCGA database with the accession numbers TCGA-HNSCC, and GEO repository (GSE139324 dataset).", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM4138110", "GSM4138111", "GSM4138112", "GSM4138113", "GSM4138114", "GSM4138115", "GSM4138116", "GSM4138117", "GSM4138118", "GSM4138119", "GSM4138120", "GSM4138121", "GSM4138122", "GSM4138123", "GSM4138124", "GSM4138125", "GSM4138126", "GSM4138127", "GSM4138128", "GSM4138129", "GSM4138130", "GSM4138131", "GSM4138132", "GSM4138133", "GSM4138134", "GSM4138135", "GSM4138136", "GSM4138137", "GSM4138138", "GSM4138139", "GSM4138140", "GSM4138141", "GSM4138142", "GSM4138143", "GSM4138144", "GSM4138145", "GSM4138146", "GSM4138147", "GSM4138148", "GSM4138149", "GSM4138150", "GSM4138151", "GSM4138152", "GSM4138153", "GSM4138154", "GSM4138155", "GSM4138156", "GSM4138157", "GSM4138158", "GSM4138159", "GSM4138160", "GSM4138161", "GSM4138162", "GSM4138163", "GSM4138164", "GSM4138165", "GSM4138166", "GSM4138167", "GSM4138168", "GSM4138169", "GSM4138170", "GSM4138171", "GSM4138172"], "sample_count": 63, "group_sizes": null, "group_sizes_basis": null, "platform": "18573", "taxon": "Homo sapiens", "title": "Immune landscape of viral- and carcinogen-drived head and neck cancer", "resolution_checked_at": "2026-11-01T18:46:34Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}], "gaps": ["asset:data:geo:GSE139324 not downloaded", "asset:data:geo:GSE139324 licence not read"]}, "analysis": {"state": "partial", "methods": [{"id": "method-1", "family": "group_comparison", "specific": "wilcoxon", "label": "Group comparison", "dsl_operations": [], "missing_operations": ["group_comparison_wilcoxon"], "expressible": false, "parameters": {}, "data_shapes": ["two_or_more_groups_continuous"], "evidence": {"score": 1.0, "band": "verbatim", "quote": "Differentially expressed genes were identified using Wilcoxon rank sum test."}}, {"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"}}], "targets": [], "gaps": ["method.group_comparison_wilcoxon"]}, "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": false, "required_operations": [], "missing_operations": ["group_comparison_wilcoxon"], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: wilcoxon."}, "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": ["GSE139324"], "resolved_accessions": ["GSE139324"], "statements": ["The data analyzed during the current study are available in the TCGA database with the accession numbers TCGA-HNSCC, and GEO repository (GSE139324 dataset)."], "basis": "The paper names GSE139324 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": [], "assets": ["source:pmid:36157879", "data:geo:GSE139324"], "steps": []}, "must_build": {"dsl_operations": [{"operation": "group_comparison_wilcoxon", "unlocks_papers": null, "priority": null}], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE139324"]}, "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:06:51Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 18573, "completion_tokens": 1690, "cached_tokens": 256}}, "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:36157879", "type": "paper", "label": "pmid:36157879"}, {"id": "source:pmid:36157879", "type": "manuscript", "label": "PMID 36157879 full text"}, {"id": "data:geo:GSE139324", "type": "dataset", "label": "GEO GSE139324"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:36157879", "type": "paper", "label": "pmid:36157879"}, {"id": "source:pmid:36157879", "type": "manuscript", "label": "PMID 36157879 full text"}, {"id": "data:geo:GSE139324", "type": "dataset", "label": "GEO GSE139324"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-36157879", "type": "extraction", "used": ["source:pmid:36157879"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE139324", "type": "accession_resolution", "used": ["data:geo:GSE139324"], "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:36157879", "parent_hash": "sha256:6fab494649446c4e5d5dfe630dcb7a18df4a7b94b8b4986110bd679affe5df78", "fork_of": null, "generation": 3}, "gaps": ["asset:data:geo:GSE139324 licence not read", "asset:data:geo:GSE139324 not downloaded", "execution.not_run", "method.group_comparison_wilcoxon", "review.human_critique", "review.independent_reader", "review.machine_critique", "step:run_analysis blocked on method:group_comparison_wilcoxon"], "assays": [{"role": "data_generation", "definition_id": "rna_seq", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.99, "evidence": ["Single-Cell Transcriptome Analysis Reveals Different Immune Signatures", "ity individual cells passed quality control. RNA-seq data were used from 522 HNSCC primary tumor", "l Survival Compared with HPV- HNSCC Patients RNA-seq data were used from 522 HNSCC primary tumor", "mmune microenvironment of HNSCC, single-cell RNA sequencing analysis was performed on the samples from d"], "order_basis": "paper_text_order:0", "id": "assay:1"}, {"role": "data_generation", "definition_id": "scrna_seq", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "unsupported", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.99, "evidence": ["define the immune microenvironment of HNSCC, single-cell RNA sequencing analysis was performed on the samples from d", "immune cells were identified (Figure 1(a)). Single-cell RNA sequencing results from immune cells in HNSCC. (a) UMAP"], "order_basis": "paper_text_order:2945", "id": "assay:2"}, {"role": "analysis", "definition_id": "group_comparison", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 1.0, "evidence": ["Differentially expressed genes were identified using Wilcoxon rank sum test."], "order_basis": "consumes_earlier_output", "id": "assay:3"}]}