Co-Expression Network Analysis Identified Gene Signatures in Osteosarcoma as a Predictive Tool for Lung Metastasis and Survival.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:31333788", "kind": "science-record", "revision": 8, "content_hash": "sha256:55704ad9678affdcd518500493ed89194d60af99c389e9288983014656d6a96b", "created_at": "2026-11-01T15:40:47Z", "updated_at": "2026-11-01T22:00:43Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Co-Expression Network Analysis Identified Gene Signatures in Osteosarcoma as a Predictive Tool for Lung Metastasis and Survival.", "doi": "10.7150/jca.32092", "pmid": "31333788", "authors": ["Honghua Zhang", "Linwei Guo", "Zheng Zhang", "Yunlong Sun", "Honglei Kang", "Chao Song", "Huiyong Liu", "Zhuowei Lei", "Jia Wang", "Baoguo Mi", "Qian Xu", "Hanfeng Guan", "Feng Li"], "journal": "Journal of Cancer", "year": 2019, "abstract": "Osteosarcoma (OS) is the most common primary bone tumor, whose poor prognosis is mainly due to lung metastasis. The aim of this study is to build a practical and valid diagnostic test that can predict the risk of OS metastasis and progression. We performed weighted gene co-expression network analysis (WGCNA) on GSE21257 from the Gene Expression Omnibus (GEO) database, which contains microarray data of biopsies from OS patients. In these modules, the highest association was found between the blue module and metastasis stage (r = -0.52) by Pearson's correlation analysis. Based on Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression, we derived eight clinically significant genes and constructed an eight-gene signature for metastasis status. It showed great efficacy to distinguish metastasis from non-metastasis (AUC = 0.886) and the results were validated in The Cancer Genome Atlas (TCGA) database. Functional enrichment analysis of hub genes showed that their biological processes focused on immune-related pathways, suggesting the important roles of immune cells, immune pathways and the tumor microenvironment in metastasis development. In conclusion, we discovered an efficient gene signature with great efficacy to distinguish metastasis status, which may help improve early diagnosis and treatment, enhancing the clinical outcomes of OS patients. Besides we created an effective protocol to seek for several hub genes in high-throughput data by combining WGCNA and LASSO Cox regression.", "claims": [{"id": "claim-method", "statement": "In the current study, WGCNA was performed on DEG-matched GSE21257, and modules were identified with a dynamic tree-cutting algorithm with a minimum module size of 10 genes, a scale-free topology threshold of 0.9 and merged with a MEDissThres parameter of 0.25.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:31333788", "locator": "sections[5] \"Constructing Dynamic Weighted Gene Co-Expression Network\" paragraph 3", "quote": "In the current study, WGCNA was performed on DEG-matched GSE21257, and modules were identified with a dynamic tree-cutting algorithm with a minimum module size of 10 genes, a scale-free topology threshold of 0.9 and merged with a MEDissThres parameter of 0.25."}], "state": "complete"}, {"id": "claim-data-access", "statement": "Preprocessed gene expression profiles of GSE42352, GSE21257 and GSE36001 were downloaded from the GEO database. GSE42352 is a microarray dataset containing 103 OS cells and 15 mesenchyme stem cells 19.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:31333788", "locator": "sections[3] \"Data Sources and Data Preprocessing\" paragraph 2", "quote": "Preprocessed gene expression profiles of GSE42352, GSE21257 and GSE36001 were downloaded from the GEO database. GSE42352 is a microarray dataset containing 103 OS cells and 15 mesenchyme stem cells 19."}], "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:31333788"], "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:GSE21257"], "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:GSE21257"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": []}, "data": {"state": "partial", "assets": [{"id": "source:pmid:31333788", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.7150/jca.32092", "sha256": "sha256:c3ad4b9c8340512a8141e5e4f9302be6e4d1979c3530527d2d5b4f8dfb4b2308", "bytes": 23095, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.7150/jca.32092", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC6636290", "source": "europe_pmc_fulltext", "coverage": "whole_body_without_references", "text_hash_scope": "rendered_section_body"}}, {"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": "Preprocessed gene expression profiles of GSE42352, GSE21257 and GSE36001 were downloaded from the GEO database. GSE42352 is a microarray dataset containing 103 OS cells and 15 mesenchyme stem cells 19.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "Preprocessed gene expression profiles of GSE42352, GSE21257 and GSE36001 were downloaded from the GEO database. GSE42352 is a microarray dataset containing 103 OS cells and 15 mesenchyme stem cells 19.", "access_statement_kind": "repository_route_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:GSE36001", "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=GSE36001", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE36001", "version": null, "note": "Preprocessed gene expression profiles of GSE42352, GSE21257 and GSE36001 were downloaded from the GEO database. GSE42352 is a microarray dataset containing 103 OS cells and 15 mesenchyme stem cells 19.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "Preprocessed gene expression profiles of GSE42352, GSE21257 and GSE36001 were downloaded from the GEO database. GSE42352 is a microarray dataset containing 103 OS cells and 15 mesenchyme stem cells 19.", "access_statement_kind": "repository_route_statement", "sample_accessions": ["GSM879156", "GSM879157", "GSM879158", "GSM879159", "GSM879160", "GSM879161", "GSM879162", "GSM879163", "GSM879164", "GSM879165", "GSM879166", "GSM879167", "GSM879168", "GSM879169", "GSM879170", "GSM879171", "GSM879172", "GSM879173", "GSM879174", "GSM879175", "GSM879176", "GSM879177", "GSM879178", "GSM879179", "GSM879180"], "sample_count": 25, "group_sizes": null, "group_sizes_basis": null, "platform": "6102", "taxon": "Homo sapiens", "title": "Integrative analysis reveals relationships of genetic and epigenetic alterations in osteosarcoma [gene expression]", "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:GSE42352", "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=GSE42352", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE42352", "version": null, "note": "Preprocessed gene expression profiles of GSE42352, GSE21257 and GSE36001 were downloaded from the GEO database. GSE42352 is a microarray dataset containing 103 OS cells and 15 mesenchyme stem cells 19.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "Preprocessed gene expression profiles of GSE42352, GSE21257 and GSE36001 were downloaded from the GEO database. GSE42352 is a microarray dataset containing 103 OS cells and 15 mesenchyme stem cells 19.", "access_statement_kind": "repository_route_statement", "sample_accessions": ["GSM1038236", "GSM1038237", "GSM1038238", "GSM1038239", "GSM1038240", "GSM1038241", "GSM1038242", "GSM1038243", "GSM1038244", "GSM1038245", "GSM1038246", "GSM1038247", "GSM1038248", "GSM1038249", "GSM1038250", "GSM1038251", "GSM1038252", "GSM1038253", "GSM1038254", "GSM717846", "GSM717847", "GSM717848", "GSM717849", "GSM717850", "GSM717851", "GSM717852", "GSM717853", "GSM717854", "GSM717855", "GSM717856", "GSM717857", "GSM825623", "GSM825624", "GSM825625", "GSM825626", "GSM825627", "GSM825628", "GSM825629", "GSM825630", "GSM825631", "GSM825632", "GSM825633", "GSM825634", "GSM825635", "GSM825636", "GSM825637", "GSM825638", "GSM825639", "GSM825640", "GSM825641", "GSM825642", "GSM825643", "GSM825644", "GSM825645", "GSM825646", "GSM825647", "GSM825648", "GSM825649", "GSM825650", "GSM825651", "GSM825652", "GSM825653", "GSM825654", "GSM825655", "GSM825656", "GSM825657", "GSM825658", "GSM825659", "GSM825660", "GSM825661", "GSM825662", "GSM825663", "GSM825664", "GSM825665", "GSM825666", "GSM825667", "GSM825668", "GSM825669", "GSM825670", "GSM825671", "GSM825672", "GSM825673", "GSM825674", "GSM825675", "GSM825676", "GSM825677", "GSM825678", "GSM825679", "GSM825680", "GSM825681", "GSM825682", "GSM825683", "GSM825684", "GSM825685", "GSM825686", "GSM825687", "GSM825688", "GSM825689", "GSM825690", "GSM825691", "GSM825692", "GSM825693", "GSM825694", "GSM825695", "GSM825696", "GSM825697", "GSM825698", "GSM825699", "GSM825700", "GSM825701", "GSM825702", "GSM825703", "GSM825704", "GSM825705", "GSM825706", "GSM825707", "GSM825708", "GSM825709"], "sample_count": 118, "group_sizes": null, "group_sizes_basis": null, "platform": "10295", "taxon": "Homo sapiens", "title": "Genome-wide gene expression profiling of mesenchymal stem cells, osteosarcoma cells, and osteosarcoma cell lines.", "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:GSE21257 not downloaded", "asset:data:geo:GSE21257 licence not read", "asset:data:geo:GSE36001 not downloaded", "asset:data:geo:GSE36001 licence not read", "asset:data:geo:GSE42352 not downloaded", "asset:data:geo:GSE42352 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": 0.8, "band": "high", "quote": "Kaplan-Meier"}}, {"id": "method-2", "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-3", "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": 0.8, "band": "high", "quote": "log-rank"}}, {"id": "method-4", "family": "correlation_analysis", "specific": "pearson", "label": "Correlation analysis", "dsl_operations": ["pearson_r"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["paired_numeric_columns"], "evidence": {"score": 0.7, "band": "medium", "quote": "Pearson's correlation"}}, {"id": "method-5", "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"}}, {"id": "method-6", "family": "contingency_analysis", "specific": "fisher_exact", "label": "Categorical (contingency table) analysis", "dsl_operations": ["contingency_table", "chi_square_test"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["categorical_cross_tabulation"], "evidence": {"score": 0.8, "band": "high", "quote": "Fisher's exact"}}, {"id": "method-7", "family": null, "specific": null, "label": "Weighted gene co-expression network analysis (WGCNA)", "dsl_operations": [], "missing_operations": [], "expressible": null, "parameters": {}, "data_shapes": [], "evidence": {"score": 1.0, "band": "verbatim", "quote": "In the current study, WGCNA was performed on DEG-matched GSE21257, and modules were identified with a dynamic tree-cutting algorithm with a minimum module size of 10 genes, a scale-free topology threshold of 0.9 and merged with a MEDissThres parameter of 0.25."}}], "targets": [], "gaps": ["method.reader_named_unmapped"]}, "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": ["GSE42352", "GSE21257", "GSE36001"], "resolved_accessions": ["GSE42352", "GSE21257", "GSE36001"], "statements": ["Preprocessed gene expression profiles of GSE42352, GSE21257 and GSE36001 were downloaded from the GEO database. GSE42352 is a microarray dataset containing 103 OS cells and 15 mesenchyme stem cells 19."], "basis": "The paper names GSE42352, GSE21257, GSE36001 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": ["survival_km"], "assets": ["source:pmid:31333788", "data:geo:GSE21257"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE21257"]}, "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:15:37Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": null}, "calls": 3, "usage": {"prompt_tokens": 16847, "completion_tokens": 2382, "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:31333788", "type": "paper", "label": "pmid:31333788"}, {"id": "source:pmid:31333788", "type": "manuscript", "label": "PMID 31333788 full text"}, {"id": "data:geo:GSE21257", "type": "dataset", "label": "GEO GSE21257"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:31333788", "type": "paper", "label": "pmid:31333788"}, {"id": "source:pmid:31333788", "type": "manuscript", "label": "PMID 31333788 full text"}, {"id": "data:geo:GSE21257", "type": "dataset", "label": "GEO GSE21257"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-31333788", "type": "extraction", "used": ["source:pmid:31333788"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE21257", "type": "accession_resolution", "used": ["data:geo:GSE21257"], "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:31333788", "parent_hash": "sha256:d36b91215eb5a58e1c1edd7705c670381837af2f35a6a93a49e435f7bfdc867c", "fork_of": null, "generation": 7}, "gaps": ["asset:data:geo:GSE21257 licence not read", "asset:data:geo:GSE21257 not downloaded", "asset:data:geo:GSE36001 licence not read", "asset:data:geo:GSE36001 not downloaded", "asset:data:geo:GSE42352 licence not read", "asset:data:geo:GSE42352 not downloaded", "execution.not_run", "method.reader_named_unmapped", "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": ["ssion Omnibus (GEO) database, which contains microarray data of biopsies from OS patients. In these"], "order_basis": "paper_text_order:470", "id": "assay:1"}, {"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": 0.8, "evidence": ["Kaplan-Meier"], "order_basis": "consumes_earlier_output", "id": "assay:2"}]}