Bioinformatic Analysis Reveals Central Role for Tumor-Infiltrating Immune Cells in Uveal Melanoma Progression.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:34195299", "kind": "science-record", "revision": 7, "content_hash": "sha256:253b94a9454596c87d79656330fc2c206ca21b5b8171309b5af25526b38850e2", "created_at": "2026-11-01T15:41:41Z", "updated_at": "2026-11-01T22:04:05Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Bioinformatic Analysis Reveals Central Role for Tumor-Infiltrating Immune Cells in Uveal Melanoma Progression.", "doi": "10.1155/2021/9920234", "pmid": "34195299", "authors": ["Mieszko Lachota", "Anton Lennikov", "Karl-Johan Malmberg", "Radoslaw Zagozdzon"], "journal": "Journal of immunology research", "year": 2021, "abstract": "Tumor-infiltrating immune cells are capable of effective cancer surveillance, and their abundance is linked to better prognosis in numerous tumor types. However, in uveal melanoma (UM), extensive immune infiltrate is associated with poor survival. This study aims to decipher the role of different tumor-infiltrating cell subsets in UM in order to identify potential targets for future immunotherapeutic treatment. We have chosen the TCGA-UVM cohort as a training dataset and GSE22138 as a testing dataset by mining publicly available databases. The abundance of 22 immune cell types was estimated using CIBERSORTx. Then, to determine the significance of tumor-infiltrating cell subsets in UM, we built a multicell type prognostic signature, which was validated in the testing cohort. The created signature was built upon the negative prognostic role of CD8+ T cells and M0 macrophages and the positive role of neutrophils. Based on the created signature score, we divided the patients into low- and high-risk groups. Kaplan-Meier, Cox, and ROC analyses demonstrated superior performance of our risk score compared to either clinical or pathologic characteristics of both cohorts. Further, we found the molecular pathways associated with cancer immunoevasion and metastasis to be enriched in the high-risk group, explaining both the lack of adequate immune surveillance despite increased infiltration of CD8+ T cells as well as the higher metastatic potential. Genes associated with tryptophan metabolism (IDO1 and KYNU) and metalloproteinases were among the most differentially expressed between the high- and low-risk groups. Our correlation analyses interpreted in context of published in vitro data strongly suggest the central role of CD8+ T cells in shifting the UM tumor microenvironment towards suppressive and metastasis-promoting. Therefore, we propose further investigations of IDO1 and metalloproteinases as novel targets for immunotherapy in lymphocyte-rich metastatic UM patients.", "claims": [{"id": "claim-method", "statement": "Kaplan–Meier analysis was used to examine the prognostic differences between the high- and low-risk groups. The p value was assessed by the two-sided log-rank test.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:34195299", "locator": "sections[8] \"2.6. Statistical Analysis\" paragraph 2", "quote": "Kaplan–Meier analysis was used to examine the prognostic differences between the high- and low-risk groups. The p value was assessed by the two-sided log-rank test."}], "state": "complete"}, {"id": "claim-data-access", "statement": "The gene expression data used in this study can be found in the The Cancer Genome Atlas repository with cohort name GDC TCGA-UVM/TCGA-PANCAN and in the GEO database with cohort name GSE22138.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:34195299", "locator": "sections[20] \"Data Availability\" paragraph 2", "quote": "The gene expression data used in this study can be found in the The Cancer Genome Atlas repository with cohort name GDC TCGA-UVM/TCGA-PANCAN and in the GEO database with cohort name GSE22138."}], "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:34195299"], "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:GSE22138"], "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:GSE22138"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": []}, "data": {"state": "partial", "assets": [{"id": "source:pmid:34195299", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1155/2021/9920234", "sha256": "sha256:0aad2fa524ae75f3052a9513624a7eb18994a00c3e181a6174fc7c20b454eec7", "bytes": 71535, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1155/2021/9920234", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC8214507", "source": "europe_pmc_fulltext", "coverage": "whole_body", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE22138", "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=GSE22138", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE22138", "version": null, "note": "The gene expression data used in this study can be found in the The Cancer Genome Atlas repository with cohort name GDC TCGA-UVM/TCGA-PANCAN and in the GEO database with cohort name GSE22138.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The gene expression data used in this study can be found in the The Cancer Genome Atlas repository with cohort name GDC TCGA-UVM/TCGA-PANCAN and in the GEO database with cohort name GSE22138.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM550623", "GSM550624", "GSM550625", "GSM550626", "GSM550627", "GSM550628", "GSM550629", "GSM550630", "GSM550631", "GSM550632", "GSM550633", "GSM550634", "GSM550635", "GSM550636", "GSM550637", "GSM550638", "GSM550639", "GSM550640", "GSM550641", "GSM550642", "GSM550643", "GSM550644", "GSM550645", "GSM550646", "GSM550647", "GSM550648", "GSM550649", "GSM550650", "GSM550651", "GSM550652", "GSM550653", "GSM550654", "GSM550655", "GSM550656", "GSM550657", "GSM550658", "GSM550659", "GSM550660", "GSM550661", "GSM550662", "GSM550663", "GSM550664", "GSM550665", "GSM550666", "GSM550667", "GSM550668", "GSM550669", "GSM550670", "GSM550671", "GSM550672", "GSM550673", "GSM550674", "GSM550675", "GSM550676", "GSM550677", "GSM550678", "GSM550679", "GSM550680", "GSM550681", "GSM550682", "GSM550683", "GSM550684", "GSM550685"], "sample_count": 63, "group_sizes": null, "group_sizes_basis": null, "platform": "570", "taxon": "Homo sapiens", "title": "Expression Data from Uveal Melanoma primary tumors.", "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:34195299:1", "role": "table", "format": "jats-table", "media_type": "application/vnd.jats.table+json", "uri": null, "sha256": "sha256:1895cd0c2add473e20f69e0b9fc3c097742b539c601036b022502b2124829ddf", "bytes": 1368, "access": "open", "state": "verified", "reader": "read_table", "reader_status": "available", "identifier": "section[11]/table-wrap[1]#tab1", "version": null, "note": "Clinical and pathological characteristics of uveal melanoma patients at the time of diagnosis. Missing data is marked as NA (not available).", "structure": {"pmcid": "PMC8214507", "section_index": 11, "section_title": "3.1. Cohort Characteristics", "row_count": 36, "column_count": 3}, "table": {"label": "Table 1", "caption": "Clinical and pathological characteristics of uveal melanoma patients at the time of diagnosis. Missing data is marked as NA (not available).", "rows": [["Characteristics", "Training cohort (TCGA-UVM, n = 71)", "Testing cohort (GSE22138, n = 60)"], ["Age at diagnosis (years)", "", ""], ["Mean (SD)", "61 (14)", "61 (13)"], ["Min-max", "22-86", "29-85"], ["Gender", "", ""], ["Female", "31 (43.66%)", "22 (36.67%)"], ["Male", "40 (56.34%)", "38 (63.33%)"], ["Clinical stage", "", ""], ["I", "0 (0.00%)", "NA"], ["II", "31 (43.66%)", "NA"], ["III", "36 (50.71%)", "NA"], ["IV", "4 (5.63%)", "NA"], ["T classification", "", ""], ["T1", "0 (0.00%)", "NA"], ["T2", "3 (4.22%)", "NA"], ["T3", "32 (45.07%)", "NA"], ["T4", "34 (47.89%)", "NA"], ["Tx", "2 (2.82%)", "NA"], ["N classification", "", ""], ["N0", "67 (94.37%)", "NA"], ["N1", "0 (0.00%)", "NA"], ["Nx", "4 (5.63%)", "NA"], ["M classification", "", ""], ["M0", "64 (90.14%)", "NA"], ["M1", "3 (4.22%)", "NA"], ["Mx", "4 (5.63%)", "NA"], ["Extrascleral extension", "", ""], ["No", "60 (84.51%)", "47 (78.33%)"], ["Yes", "7 (9.86%)", "4 (6.67%)"], ["Unknown", "4 (5.63%)", "9 (15.00%)"], ["Tumor basal diameter (mm)", "", ""], ["Mean (SD)", "17 (4)", "15 (4)"], ["Range", "8-25", "9-23"], ["Tumor thickness (mm)", "", ""], ["Mean (SD)", "11 (3)", "12 (2)"], ["Range", "4-16", "6-17"]], "row_count": 36, "column_count": 3, "numeric_cells": 39, "header_rows": 1, "locator": "section[11]/table-wrap[1]#tab1"}}, {"id": "asset:table:pmid:34195299:2", "role": "table", "format": "jats-table", "media_type": "application/vnd.jats.table+json", "uri": null, "sha256": "sha256:18d0ffdaf5b7d0e5943602905d0b895a9654e66fa4669a5becdcc1ff15c46a15", "bytes": 1265, "access": "open", "state": "verified", "reader": "read_table", "reader_status": "available", "identifier": "section[16]/table-wrap[1]#tab2", "version": null, "note": "Selected differentially expressed genes (DEGs) between high- and low-risk groups in the training cohort. The genes were selected based on the biological function.", "structure": {"pmcid": "PMC8214507", "section_index": 16, "section_title": "3.6. The Immune Architecture of the Tumors in High- and Low-Risk Group", "row_count": 14, "column_count": 4}, "table": {"label": "Table 2", "caption": "Selected differentially expressed genes (DEGs) between high- and low-risk groups in the training cohort. The genes were selected based on the biological function.", "rows": [["Gene symbol", "Gene name", "Fold change (log2)", "Adjusted p value"], ["KYNU", "Kynureninase", "-1.755", "0.00786"], ["IDO1", "Indoleamine 2,3-dioxygenase 1", "2.326", "0.01879"], ["MMP9", "Matrix metalloproteinase 9", "2.745", "<0.0001"], ["MMP1", "Matrix metalloproteinase 1", "1.761", "0.029"], ["MMP12", "Matrix metalloproteinase 12", "4.562", "0.00065"], ["MMP25", "Matrix metalloproteinase 25", "1.690", "0.00054"], ["ADAMTS2", "ADAM metallopeptidase with thrombospondin type 1 motif 2", "1.696", "0.00251"], ["ADAMDEC1", "ADAM like decysin 1", "2.863", "0.005"], ["ADAM11", "ADAM metallopeptidase domain 11", "0.975", "0.0476"], ["ADAMTS8", "ADAM metallopeptidase with thrombospondin type 1 motif 8", "0.931", "0.0207"], ["ADAMTS4", "ADAM metallopeptidase with thrombospondin type 1 motif 4", "0.928", "0.0047"], ["ADAMTS9", "ADAM metallopeptidase with thrombospondin type 1 motif 9", "0.962", "0.0064"], ["ADAMTS15", "ADAM metallopeptidase with thrombospondin type 1 motif 15", "0.969", "0.0056"]], "row_count": 14, "column_count": 4, "numeric_cells": 42, "header_rows": 1, "locator": "section[16]/table-wrap[1]#tab2"}}], "gaps": ["asset:data:geo:GSE22138 not downloaded", "asset:data:geo:GSE22138 licence not read"]}, "analysis": {"state": "partial", "methods": [{"id": "method-1", "family": "survival_analysis", "specific": "log_rank_test", "label": "Survival analysis", "dsl_operations": ["logrank_test"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["time_to_event_with_censoring"], "evidence": {"score": 1.0, "band": "verbatim", "quote": "Kaplan–Meier analysis was used to examine the prognostic differences between the high- and low-risk groups. The p value was assessed by the two-sided log-rank test."}}, {"id": "method-2", "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": 0.8, "band": "high", "quote": "Wilcoxon"}}, {"id": "method-3", "family": "survival_analysis", "specific": "kaplan_meier", "label": "Survival analysis", "dsl_operations": ["survival_km"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["time_to_event_with_censoring"], "evidence": {"score": 0.7, "band": "medium", "quote": "Kaplan-Meier"}}, {"id": "method-4", "family": "survival_analysis", "specific": "cox_proportional_hazards", "label": "Survival analysis", "dsl_operations": ["cox_ph"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["time_to_event_with_censoring"], "evidence": {"score": 0.8, "band": "high", "quote": "Cox regression"}}, {"id": "method-5", "family": "differential_expression", "specific": "deseq2", "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.8, "band": "high", "quote": "DESeq2"}}], "targets": [], "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": ["logrank_test"], "missing_operations": [], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: log_rank_test."}, "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": ["GSE22138"], "resolved_accessions": ["GSE22138"], "statements": ["The gene expression data used in this study can be found in the The Cancer Genome Atlas repository with cohort name GDC TCGA-UVM/TCGA-PANCAN and in the GEO database with cohort name GSE22138."], "basis": "The paper names GSE22138 in GEO database and the accession resolves in the repository. No dataset was downloaded and no licence was read."}}, "can_run": {"dsl_operations": ["logrank_test"], "assets": ["source:pmid:34195299", "data:geo:GSE22138", "asset:table:pmid:34195299:1", "asset:table:pmid:34195299:2"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE22138"]}, "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:31:39Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 3, "usage": {"prompt_tokens": 57025, "completion_tokens": 2090, "cached_tokens": 784}}, "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:34195299", "type": "paper", "label": "pmid:34195299"}, {"id": "source:pmid:34195299", "type": "manuscript", "label": "PMID 34195299 full text"}, {"id": "data:geo:GSE22138", "type": "dataset", "label": "GEO GSE22138"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:34195299", "type": "paper", "label": "pmid:34195299"}, {"id": "source:pmid:34195299", "type": "manuscript", "label": "PMID 34195299 full text"}, {"id": "data:geo:GSE22138", "type": "dataset", "label": "GEO GSE22138"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-34195299", "type": "extraction", "used": ["source:pmid:34195299"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE22138", "type": "accession_resolution", "used": ["data:geo:GSE22138"], "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:34195299", "parent_hash": "sha256:7028795b9f61df8e04f5ff505aa85240694d3e66359b615fa8794827fa10df44", "fork_of": null, "generation": 6}, "gaps": ["asset:data:geo:GSE22138 licence not read", "asset:data:geo:GSE22138 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": "rna_seq", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.93, "evidence": ["ods 2.4. Gene Expression Data Processing The RNA-seq expression data was HTSeq-counts type. After"], "order_basis": "paper_text_order:2127", "id": "assay:1"}, {"role": "analysis", "definition_id": "survival_analysis", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": ["logrank_test"], "support_state": "supported", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 1.0, "evidence": ["Kaplan–Meier analysis was used to examine the prognostic differences between the high- and low-risk groups. The p value was assessed by the two-sided log-rank test."], "order_basis": "consumes_earlier_output", "id": "assay:2"}]}