A Glycerophospholipid Metabolism-Based Prognostic Model Guides Osteosarcoma Therapy.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:41172875", "kind": "science-record", "revision": 4, "content_hash": "sha256:6aea5a7f046028bc24fed99516b41432265e5c96be31c7ed120c509343ca3b66", "created_at": "2026-11-01T20:31:13Z", "updated_at": "2026-11-01T22:57:23Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "A Glycerophospholipid Metabolism-Based Prognostic Model Guides Osteosarcoma Therapy.", "doi": "10.1016/j.tranon.2025.102572", "pmid": "41172875", "authors": ["Yong Wen", "Donglian Wang", "Hongshen Wang", "Jian Zhao", "Li Yao", "Jin Zhang", "Bolai Chen"], "journal": "Translational oncology", "year": 2026, "abstract": "BACKGROUND: Glycerophospholipid metabolism regulates membrane synthesis, signaling, and cell homeostasis, but its role in osteosarcoma (OS) progression, immunity, and therapy remains unclear.\nMETHODS: We integrated bulk transcriptome data from TARGET (n = 85) and GEO cohorts (n = 124) with single-cell RNA-seq data (GSE162454, six OS samples). Glycerophospholipid metabolism activity was quantified using five algorithms to generate a GAS score. Prognostic genes were identified via univariate Cox regression, followed by Lasso-Cox modeling. Downstream analyses included cell communication, immune infiltration, mutation profiling, drug sensitivity, and immunotherapy response.\nRESULTS: Twelve cell types were identified; malignant OS cells showed the lowest glycerophospholipid metabolism. Nine genes correlated with GAS (r=±0.2, p < 0.01), six were used to construct the model, which outperformed 42 published signatures (C-index>0.7). High GAS patients exhibited elevated immune gene expression and copy number alterations (p < 0.05). GAS remained an independent prognostic factor (HR=3.58, 95 %CI:2.23-5.74, p < 0.001). Drug prediction highlighted lovastatin, simvastatin, and tamatinib for high-risk patients, and higher GAS associated with poor immunotherapy response (p < 0.01).\nCONCLUSION: GAS-based stratification not only robustly predicts OS prognosis but also reveals interactions between glycerophospholipid metabolism and tumor immunity, guiding personalized therapeutic strategies and highlighting novel drug candidates for high-risk patients.", "claims": [{"id": "claim-method", "statement": "Survival analyses were conducted using Kaplan-Meier and Cox regression models.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:41172875", "locator": "sections[9] \"Statistical analysis\" paragraph 2", "quote": "Survival analyses were conducted using Kaplan-Meier and Cox regression models."}], "state": "complete"}, {"id": "claim-data-access", "statement": "All datasets utilized or generated during the study are fully available within this article.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:41172875", "locator": "sections[25] \"Data availability\" paragraph 2", "quote": "All datasets utilized or generated during the study are fully available within this article."}], "state": "complete"}, {"id": "claim-target-result", "statement": "GAS remained an independent prognostic factor (HR=3.58, 95 %CI:2.23–5.74, p < 0.001).", "kind": "quantitative", "target": {"value": 95.0, "unit": "per cent", "tolerance": null, "source": "GAS remained an independent prognostic factor (HR=3.58, 95 %CI:2.23–5.74, p < 0.001)."}, "evidence": [{"asset": "source:pmid:41172875", "locator": "sections[0] \"Abstract\" paragraph 10", "quote": "GAS remained an independent prognostic factor (HR=3.58, 95 %CI:2.23–5.74, p < 0.001)."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "complete", "objective": "GAS remained an independent prognostic factor (HR=3.58, 95 %CI:2.23–5.74, p < 0.001).", "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:41172875"], "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:GSE16091"], "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:GSE16091"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": []}, "data": {"state": "partial", "assets": [{"id": "source:pmid:41172875", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1016/j.tranon.2025.102572", "sha256": "sha256:bb80b3b6d0aacab157fc8435d65592976f9264f6dc83fc25de9e7f2093fbd7ee", "bytes": 30990, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1016/j.tranon.2025.102572", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC12613023", "source": "europe_pmc_fulltext", "coverage": "whole_body_without_references", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE16091", "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=GSE16091", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE16091", "version": null, "note": "All datasets utilized or generated during the study are fully available within this article.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "All datasets utilized or generated during the study are fully available within this article.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM402718", "GSM402719", "GSM402720", "GSM402721", "GSM402722", "GSM402723", "GSM402724", "GSM402725", "GSM402726", "GSM402727", "GSM402728", "GSM402729", "GSM402730", "GSM402731", "GSM402732", "GSM402733", "GSM402734", "GSM402735", "GSM402736", "GSM402737", "GSM402738", "GSM402739", "GSM402740", "GSM402741", "GSM402742", "GSM402743", "GSM402744", "GSM402745", "GSM402746", "GSM402747", "GSM402748", "GSM402749", "GSM402750", "GSM402751"], "sample_count": 34, "group_sizes": null, "group_sizes_basis": null, "platform": "96", "taxon": "Homo sapiens", "title": "Gene expression profiles of human osteosarcoma, set2", "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:GSE162454", "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=GSE162454", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE162454", "version": null, "note": "All datasets utilized or generated during the study are fully available within this article.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "All datasets utilized or generated during the study are fully available within this article.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM4952363", "GSM4952364", "GSM4952365", "GSM5155198", "GSM5155199", "GSM5155200"], "sample_count": 6, "group_sizes": null, "group_sizes_basis": null, "platform": "24676", "taxon": "Homo sapiens", "title": "Single-cell RNA sequencing of human osteosarcoma", "resolution_checked_at": "2026-11-01T20:16:08Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}, {"id": "data:geo:GSE21257", "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=GSE21257", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE21257", "version": null, "note": "All datasets utilized or generated during the study are fully available within this article.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "All datasets utilized or generated during the study are fully available within this article.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM530667", "GSM530899", "GSM531283", "GSM531284", "GSM531285", "GSM531286", "GSM531287", "GSM531288", "GSM531289", "GSM531290", "GSM531291", "GSM531292", "GSM531293", "GSM531294", "GSM531295", "GSM531296", "GSM531297", "GSM531298", "GSM531299", "GSM531300", "GSM531301", "GSM531302", "GSM531303", "GSM531304", "GSM531305", "GSM531306", "GSM531307", "GSM531308", "GSM531309", "GSM531310", "GSM531311", "GSM531312", "GSM531313", "GSM531314", "GSM531319", "GSM531320", "GSM531321", "GSM531322", "GSM531323", "GSM531324", "GSM531325", "GSM531326", "GSM531327", "GSM531328", "GSM531329", "GSM531330", "GSM531331", "GSM531332", "GSM531333", "GSM531334", "GSM531335", "GSM531351", "GSM531352"], "sample_count": 53, "group_sizes": null, "group_sizes_basis": null, "platform": "10295", "taxon": "Homo sapiens", "title": "Genome-wide gene expression profiling on pre-chemotherapy biopsies of osteosarcoma patients who developed metastases within 5yrs (n=34) and patients who did not develop metastases within 5yrs (n=19)", "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:GSE39055", "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=GSE39055", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE39055", "version": null, "note": "All datasets utilized or generated during the study are fully available within this article.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "All datasets utilized or generated during the study are fully available within this article.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM954790", "GSM954791", "GSM954792", "GSM954793", "GSM954794", "GSM954795", "GSM954796", "GSM954797", "GSM954798", "GSM954799", "GSM954800", "GSM954801", "GSM954802", "GSM954803", "GSM954804", "GSM954805", "GSM954806", "GSM954807", "GSM954808", "GSM954809", "GSM954810", "GSM954811", "GSM954812", "GSM954813", "GSM954814", "GSM954815", "GSM954816", "GSM954817", "GSM954818", "GSM954819", "GSM954820", "GSM954821", "GSM954822", "GSM954823", "GSM954824", "GSM954825", "GSM954826"], "sample_count": 37, "group_sizes": null, "group_sizes_basis": null, "platform": "14951", "taxon": "Homo sapiens", "title": "mRNA profiling and clinical outcomes in human osteosarcoma (biopsies)", "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:GSE16091 not downloaded", "asset:data:geo:GSE16091 licence not read", "asset:data:geo:GSE162454 not downloaded", "asset:data:geo:GSE162454 licence not read", "asset:data:geo:GSE21257 not downloaded", "asset:data:geo:GSE21257 licence not read", "asset:data:geo:GSE39055 not downloaded", "asset:data:geo:GSE39055 licence not read"]}, "analysis": {"state": "partial", "methods": [{"id": "method-1", "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": 1.0, "band": "verbatim", "quote": "Survival analyses were conducted using Kaplan-Meier and Cox regression models."}}, {"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": "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.7, "band": "medium", "quote": "Cox regression"}}, {"id": "method-4", "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": "PCA"}}, {"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.8, "band": "high", "quote": "ROC analysis"}}, {"id": "method-6", "family": "contingency_analysis", "specific": "chi_square", "label": "Categorical (contingency table) analysis", "dsl_operations": ["chi_square_test"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["categorical_cross_tabulation"], "evidence": {"score": 0.8, "band": "high", "quote": "chi-square"}}], "targets": [{"id": "target-published-estimate", "claim_id": "claim-target-result", "operation": "chi_square_test", "expected_value": 95.0, "unit": "per cent", "tolerance": null, "state": "complete", "note": "The published estimate as the reader quoted it from the paper. No tolerance is published because the paper states none."}, {"id": "target-published-estimate-hazard-ratio", "claim_id": "claim-target-result", "operation": "chi_square_test", "expected_value": 3.58, "unit": "hazard ratio", "tolerance": null, "state": "complete", "note": "The published estimate as the reader quoted it from the paper. No tolerance is published because the paper states none."}, {"id": "target-published-estimate-p-value", "claim_id": "claim-target-result", "operation": "chi_square_test", "expected_value": 0.001, "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": ["survival_km"], "missing_operations": [], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: kaplan_meier."}, "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": ["GSE16091", "GSE21257", "GSE39055", "GSE162454"], "resolved_accessions": ["GSE16091", "GSE21257", "GSE39055", "GSE162454"], "statements": ["All datasets utilized or generated during the study are fully available within this article."], "basis": "The paper names GSE16091, GSE21257, GSE39055, GSE162454 in GEO and the accession resolves in the repository. No dataset was downloaded and no licence was read."}}, "can_run": {"dsl_operations": ["survival_km"], "assets": ["source:pmid:41172875", "data:geo:GSE16091"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE16091"]}, "gaps": []}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-11-01T20:04:50Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 8677, "completion_tokens": 1432, "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:41172875", "type": "paper", "label": "pmid:41172875"}, {"id": "source:pmid:41172875", "type": "manuscript", "label": "PMID 41172875 full text"}, {"id": "data:geo:GSE16091", "type": "dataset", "label": "GEO GSE16091"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:41172875", "type": "paper", "label": "pmid:41172875"}, {"id": "source:pmid:41172875", "type": "manuscript", "label": "PMID 41172875 full text"}, {"id": "data:geo:GSE16091", "type": "dataset", "label": "GEO GSE16091"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-41172875", "type": "extraction", "used": ["source:pmid:41172875"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE16091", "type": "accession_resolution", "used": ["data:geo:GSE16091"], "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:41172875", "parent_hash": "sha256:287b6017fee3758c3a8b0f852f02b2b470d3e6dc5ced6699f608763c3ad78c55", "fork_of": null, "generation": 3}, "gaps": ["asset:data:geo:GSE16091 licence not read", "asset:data:geo:GSE16091 not downloaded", "asset:data:geo:GSE162454 licence not read", "asset:data:geo:GSE162454 not downloaded", "asset:data:geo:GSE21257 licence not read", "asset:data:geo:GSE21257 not downloaded", "asset:data:geo:GSE39055 licence not read", "asset:data:geo:GSE39055 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.99, "evidence": [") and GEO cohorts (n = 124) with single-cell RNA-seq data (GSE162454, six OS samples). Glyceropho", "remains unclear. METHODS: We integrated bulk transcriptome data from TARGET (n = 85) and GEO cohorts (n"], "order_basis": "paper_text_order:260", "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.93, "evidence": ["RGET (n = 85) and GEO cohorts (n = 124) with single-cell RNA-seq data (GSE162454, six OS samples). Glyceropho"], "order_basis": "paper_text_order:331", "id": "assay:2"}, {"role": "analysis", "definition_id": "survival_analysis", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": ["survival_km"], "support_state": "supported", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 1.0, "evidence": ["Survival analyses were conducted using Kaplan-Meier and Cox regression models."], "order_basis": "consumes_earlier_output", "id": "assay:3"}]}