Multiple microarray analyses identify key genes associated with the development of Non-Small Cell Lung Cancer from Chronic Obstructive Pulmonary Disease.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:33442399", "kind": "science-record", "revision": 7, "content_hash": "sha256:eec0fc774cfd20c104dd3917d17c4d7ed27e6284fade1bfc9ea92adfeee47911", "created_at": "2026-11-01T15:41:22Z", "updated_at": "2026-11-01T22:02:56Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Multiple microarray analyses identify key genes associated with the development of Non-Small Cell Lung Cancer from Chronic Obstructive Pulmonary Disease.", "doi": "10.7150/jca.51264", "pmid": "33442399", "authors": ["Lemeng Zhang", "Jianhua Chen", "Hua Yang", "Changqie Pan", "Haitao Li", "Yongzhong Luo", "Tianli Cheng"], "journal": "Journal of Cancer", "year": 2021, "abstract": "Introduction: Chronic obstructive pulmonary disease (COPD) is an independent risk factor of non-small cell lung cancer (NSCLC). This study aimed to analyze the key genes and potential molecular mechanisms that are involved in the development from COPD to NSCLC. Methods: Expression profiles of COPD and NSCLC in GSE106899, GSE12472, and GSE12428 were downloaded from the Gene Expression Omnibus (GEO) database, followed by identification of the differentially expressed genes (DEGs) between COPD and NSCLC. Based on the identified DEGs, functional pathway enrichment and lung carcinogenesis-related networks analyses were performed and further visualized with Cytoscape software. Then, principal component analysis (PCA), cluster analysis, and support vector machines (SVM) verified the ability of the top modular genes to distinguish COPD from NSCLC. Additionally, the corrections between these key genes and clinical staging of NSCLC were studied using the UALCAN and HPA websites. Finally, a prognostic risk model was constructed based on multivariate Cox regression analysis. Kaplan-Meier survival curves of the top modular genes on the training and verification sets were generated. Results: A total of 2350, 1914, and 1850 DEGs were obtained from GSE106899, GSE12472, and GSE12428 datasets, respectively. Following analysis of protein-protein interaction networks, the identified modular gene signatures containing H2AFX, MCM2, MCM3, MCM7, POLD1, and RPA1 were identified as markers for discrimination between COPD and NSCLC. The modular gene signatures were mainly enriched in the processes of DNA replication, cell cycle, mismatch repair, and others. Besides, the expression levels of these genes were significantly higher in NSCLC than in COPD, which was further verified by the immunohistochemistry. In addition, the high expression levels of H2AFX, MCM2, MCM7, and POLD1 correlate with poor prognosis of lung adenocarcinoma (LUAD). The Cox regression prognostic risk model showed the similar results and the predictive ability of this model is independent of other clinical variables. Conclusions: This study revealed several key modules that closely relate to NSCLC with underlying disease COPD, which provide a deeper understanding of the potential mechanisms underlying the malignant development from COPD to NSCLC. This study provides valuable prognostic factors in high-risk lung cancer patients with COPD.", "claims": [{"id": "claim-method", "statement": "Multivariate Cox regression analysis was performed on the modular gene signatures, and the prognosis risk model was constructed according to the following formula:", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:33442399", "locator": "sections[14] \"Construction of the prognostic risk model based on the Cox regression analysis\" paragraph 3", "quote": "Multivariate Cox regression analysis was performed on the modular gene signatures, and the prognosis risk model was constructed according to the following formula:"}], "state": "complete"}, {"id": "claim-data-access", "statement": "The raw data, supporting the conclusions of this manuscript, will be made available by contacting the corresponding author.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:33442399", "locator": "sections[33] \"Data Availability Statement\" paragraph 2", "quote": "The raw data, supporting the conclusions of this manuscript, will be made available by contacting the corresponding author."}], "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:33442399"], "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:33442399", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.7150/jca.51264", "sha256": "sha256:f1e7545acc9d933c466d6e8d20e4d001291e8bd9edab35cb92af18dd55125197", "bytes": 53991, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.7150/jca.51264", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC7797649", "source": "europe_pmc_fulltext", "coverage": "whole_body", "text_hash_scope": "rendered_section_body"}}, {"id": "asset:table:pmid:33442399:1", "role": "table", "format": "jats-table", "media_type": "application/vnd.jats.table+json", "uri": null, "sha256": "sha256:7464369201892d67e535a8f0becdd315796758ac33cbda8836ad4d198efeea8b", "bytes": 941, "access": "open", "state": "verified", "reader": "read_table", "reader_status": "available", "identifier": "section[26]/table-wrap[1]#T1", "version": null, "note": "Univariate and multivariate cox regression analysis of clinical phenotypes on LUAD samples", "structure": {"pmcid": "PMC7797649", "section_index": 26, "section_title": "Univariate and multivariate Cox regression analysis of clinical variables and risk group", "row_count": 12, "column_count": 7}, "table": {"label": "Table 1", "caption": "Univariate and multivariate cox regression analysis of clinical phenotypes on LUAD samples", "rows": [["Variables", "Univariate analysis", "Multivariate analysis"], ["HR", "95 % CI", "P value", "HR", "95 % CI", "P value"], ["Gender", "", "", "", "", "", ""], ["Male/Female", "2.068", "1.264-3.382", "0.003807702", "1.852", "1.124-3.054", "0.01565781"], ["Age", "", "", "", "", "", ""], [">60/≤60", "1.273", "0.752-2.158", "0.368987051", "1.293", "0.754-2.216", "0.35075578"], ["TNM stage", "", "", "", "", "", ""], ["Stage III-IV/Stage I-II", "3.019", "1.760-5.179", "0.00006", "2.414", "1.374-4.241", "0.00216557"], ["Smoke years", "", "", "", "", "", ""], [">15/≤15", "1.194", "0.569-2.503", "0.639280745", "0.990", "0.464-2.112", "0.98036512"], ["Risk group", "", "", "", "", "", ""], ["High/Low", "2.412", "1.320-4.407", "0.003807702", "2.217", "1.198-4.102", "0.01124920"]], "row_count": 12, "column_count": 7, "numeric_cells": 34, "header_rows": 2, "locator": "section[26]/table-wrap[1]#T1"}}], "gaps": []}, "analysis": {"state": "complete", "methods": [{"id": "method-1", "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": 1.0, "band": "verbatim", "quote": "Multivariate Cox regression analysis was performed on the modular gene signatures, and the prognosis risk model was constructed according to the following formula:"}}, {"id": "method-2", "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-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": "logRank"}}, {"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.7, "band": "medium", "quote": "principal component analysis"}}, {"id": "method-5", "family": "machine_learning", "specific": "support_vector_machine", "label": "Machine learning", "dsl_operations": ["train_classifier", "cross_validate", "roc_auc"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["labelled_feature_matrix"], "evidence": {"score": 0.7, "band": "medium", "quote": "support vector machines"}}, {"id": "method-6", "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-7", "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.7, "band": "medium", "quote": "pathway 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": ["cox_ph"], "missing_operations": [], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: cox_proportional_hazards."}, "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 raw data, supporting the conclusions of this manuscript, will be made available by contacting the corresponding author."], "basis": "The paper names no accession to its source data, so the data has to come from the authors."}}, "can_run": {"dsl_operations": ["cox_ph"], "assets": ["source:pmid:33442399", "asset:table:pmid:33442399: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:21:59Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 3, "usage": {"prompt_tokens": 42587, "completion_tokens": 2705, "cached_tokens": 792}}, "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:33442399", "type": "paper", "label": "pmid:33442399"}, {"id": "source:pmid:33442399", "type": "manuscript", "label": "PMID 33442399 full text"}, {"id": "data:geo:None", "type": "dataset", "label": "GEO None"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:33442399", "type": "paper", "label": "pmid:33442399"}, {"id": "source:pmid:33442399", "type": "manuscript", "label": "PMID 33442399 full text"}, {"id": "data:geo:None", "type": "dataset", "label": "GEO None"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-33442399", "type": "extraction", "used": ["source:pmid:33442399"], "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:33442399", "parent_hash": "sha256:f464f9a9f44449573eb3173be6d4932bcba581cd75a7bc94dcb6aaa660e9e608", "fork_of": null, "generation": 6}, "gaps": ["execution.not_run", "review.human_critique", "review.independent_reader", "review.machine_critique", "step:obtain_source_data blocked on data:accession_absent"], "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.99, "evidence": ["Multiple microarray analyses identify key genes associated with", "Materials and Methods Retrieval of microarray datasets and the workflow of this study “COP", "abase (http://www.ncbi.nlm.nih.gov/geo/) 17. Microarray data that met the following criteria were in"], "order_basis": "paper_text_order:0", "id": "assay:1"}, {"role": "data_generation", "definition_id": "ihc", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.93, "evidence": ["n in COPD, which was further verified by the immunohistochemistry. In addition, the high expression levels of"], "order_basis": "paper_text_order:1897", "id": "assay:2"}, {"role": "analysis", "definition_id": "survival_analysis", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": ["cox_ph"], "support_state": "supported", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 1.0, "evidence": ["Multivariate Cox regression analysis was performed on the modular gene signatures, and the prognosis risk model was constructed according to the following formula:"], "order_basis": "consumes_earlier_output", "id": "assay:3"}]}