Engrailed 1 promotes immune evasion and chemoresistance in glioma through single cell and CeRNA network analyses.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:41513764", "kind": "science-record", "revision": 3, "content_hash": "sha256:9b4cad40662f1ec615d544a4e91076ba22a91b5cf2ef8f563b60a25f9a7d5522", "created_at": "2026-11-01T20:31:29Z", "updated_at": "2026-11-01T22:15:17Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Engrailed 1 promotes immune evasion and chemoresistance in glioma through single cell and CeRNA network analyses.", "doi": "10.1038/s41598-026-35553-y", "pmid": "41513764", "authors": ["Zhen Jia", "Yanchao Wang", "Jieqin Yao", "Liuyue Zhang", "Xingbo Liang", "Tao Wen", "Wenle Li", "Zhenzhe Li"], "journal": "Scientific reports", "year": 2026, "abstract": "Immune evasion is a major barrier to effective glioma treatment, yet the molecular drivers contributing to this process remain insufficiently understood. Engrailed-1 (EN1), a developmental transcription factor, has recently emerged as a potential oncogenic regulator, but its role in glioma immune modulation is unclear. This study aimed to comprehensively investigate the function of EN1 in glioma, with a particular focus on its involvement in immune evasion, by integrating bulk and single-cell transcriptomic analyses together with ceRNA regulatory network construction and experimental validation. EN1 expression, diagnostic value, and prognostic significance were assessed using TCGA, GTEx, and CGGA datasets. Immune infiltration and microenvironmental features were evaluated through CIBERSORT, ESTIMATE, and WGCNA. Single-cell RNA-seq data (GSE182109) were used to characterize the cell-type distribution and developmental trajectory of EN1. A NEAT1/miR-9-5p/miR-128-3p/EN1 ceRNA axis was constructed using multi-database predictions. Functional assays including Western blot, CCK-8, and Transwell experiments were performed to validate EN1's effects in glioblastoma cells. EN1 was markedly overexpressed in glioma and associated with poorer survival. Elevated EN1 expression correlated with increased infiltration of immunosuppressive cells, reduced tumor purity, and higher immune checkpoint expression. Single-cell analysis revealed progressive EN1 upregulation along malignant cell pseudotime. Drug sensitivity prediction suggested that the EN1-high group may have reduced sensitivity to temozolomide and additional chemotherapeutic agents. The identified NEAT1/miR-9-5p/miR-128-3p/EN1 ceRNA loop suggested a regulatory mechanism contributing to EN1 activation. EN1 knockdown significantly suppressed glioblastoma cell proliferation and invasion in vitro. EN1 is associated with glioma aggressiveness, immune microenvironmental features, and predicted chemotherapeutic response, and a putative NEAT1/miR-9-5p/miR-128-3p/EN1 axis may contribute to EN1 dysregulation. These findings identify EN1 as a promising biomarker and potential therapeutic target for improving glioma treatment.", "claims": [{"id": "claim-method", "statement": "To determine the diagnostic performance of EN1, we performed receiver operating characteristic (ROC) curve analysis, calculating the area under the ROC curve (AUC). Tumor data from TCGA and normal tissue data from GTEx were used.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:41513764", "locator": "sections[5] \"Diagnostic value analysis\" paragraph 2", "quote": "To determine the diagnostic performance of EN1, we performed receiver operating characteristic (ROC) curve analysis, calculating the area under the ROC curve (AUC). Tumor data from TCGA and normal tissue data from GTEx were used."}], "state": "complete"}, {"id": "claim-data-access", "statement": "Data availabilityThe bulk RNA-seq and corresponding clinical data analyzed in this study were obtained from The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/), including the TCGA-GBM and TCGA-LGG cohorts. Normal brain tissue expression data were retrieved from the Genotype-Tissue Expression (GTEx) project (https://gtexportal.org). An independent validation cohort was obtained from the Chinese Glioma Genome Atlas (CGGA, http://www.cgga.org.cn/), specifically the CGGA_693 dataset. Single-cell RNA sequencing data were downloaded from the Gene Expression Omnibus (GEO) database under accession number GSE182109. All datasets used in this study are publicly available and fully de-identified.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:41513764", "locator": "sections[33] \"notes\" paragraph 1", "quote": "Data availabilityThe bulk RNA-seq and corresponding clinical data analyzed in this study were obtained from The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/), including the TCGA-GBM and TCGA-LGG cohorts. Normal brain tissue expression data were retrieved from the Genotype-Tissue Expression (GTEx) project (https://gtexportal.org). An independent validation cohort was obtained from the Chinese Glioma Genome Atlas (CGGA, http://www.cgga.org.cn/), specifically the CGGA_693 dataset. Single-cell RNA sequencing data were downloaded from the Gene Expression Omnibus (GEO) database under accession number GSE182109. All datasets used in this study are publicly available and fully de-identified."}], "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:41513764"], "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:GSE182109"], "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:GSE182109"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": []}, "data": {"state": "partial", "assets": [{"id": "source:pmid:41513764", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1038/s41598-026-35553-y", "sha256": "sha256:e9a50ad7b23af4a5a6dad6b0730c565752406e4ae37d16341505df4a4be48298", "bytes": 33675, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1038/s41598-026-35553-y", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC12873234", "source": "europe_pmc_fulltext", "coverage": "whole_body_without_references", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE182109", "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=GSE182109", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE182109", "version": null, "note": "Data availabilityThe bulk RNA-seq and corresponding clinical data analyzed in this study were obtained from The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/), including the TCGA-GBM and TCGA-LGG cohorts. Normal brain tissue expression data were retrieved from the Genotype-Tissue Expression (GTEx) project (https://gtexportal.org). An independent validation cohort was obtained from the Chinese Glioma Genome Atlas (CGGA, http://www.cgga.org.cn/), specifically the CGGA_693 dataset. Single-cell RNA sequencing data were downloaded from the Gene Expression Omnibus (GEO) database under accession number GSE182109. All datasets used in this study are publicly available and fully de-identified.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "Data availabilityThe bulk RNA-seq and corresponding clinical data analyzed in this study were obtained from The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/), including the TCGA-GBM and TCGA-LGG cohorts. Normal brain tissue expression data were retrieved from the Genotype-Tissue Expression (GTEx) project (https://gtexportal.org). An independent validation cohort was obtained from the Chinese Glioma Genome Atlas (CGGA, http://www.cgga.org.cn/), specifically the CGGA_693 dataset. Single-cell RNA sequencing data were downloaded from the Gene Expression Omnibus (GEO) database under accession number GSE182109. All datasets used in this study are publicly available and fully de-identified.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM5518596", "GSM5518597", "GSM5518598", "GSM5518599", "GSM5518600", "GSM5518601", "GSM5518602", "GSM5518603", "GSM5518604", "GSM5518605", "GSM5518606", "GSM5518607", "GSM5518608", "GSM5518609", "GSM5518610", "GSM5518611", "GSM5518612", "GSM5518613", "GSM5518614", "GSM5518615", "GSM5518616", "GSM5518617", "GSM5518618", "GSM5518619", "GSM5518620", "GSM5518621", "GSM5518622", "GSM5518623", "GSM5518624", "GSM5518625", "GSM5518626", "GSM5518627", "GSM5518628", "GSM5518629", "GSM5518630", "GSM5518631", "GSM5518632", "GSM5518633", "GSM5518634", "GSM5518635", "GSM5518636", "GSM5518637", "GSM5518638", "GSM5518639"], "sample_count": 44, "group_sizes": null, "group_sizes_basis": null, "platform": "20301;28038", "taxon": "Homo sapiens", "title": "A glioma immune cell atlas identifies prognostic macrophage signatures and a novel immunotherapy target, S100A4", "resolution_checked_at": "2026-11-01T19:09:46Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}], "gaps": ["asset:data:geo:GSE182109 not downloaded", "asset:data:geo:GSE182109 licence not read"]}, "analysis": {"state": "partial", "methods": [{"id": "method-1", "family": "group_comparison", "specific": "anova", "label": "Group comparison", "dsl_operations": ["group_summary", "welch_t", "bootstrap_mean_difference", "group_comparison_anova"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["two_or_more_groups_continuous"], "evidence": {"score": 0.8, "band": "high", "quote": "ANOVA"}}, {"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 component analysis"}}, {"id": "method-3", "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 curve"}}, {"id": "method-4", "family": null, "specific": null, "label": "receiver operating characteristic (ROC) curve analysis", "dsl_operations": [], "missing_operations": [], "expressible": null, "parameters": {}, "data_shapes": [], "evidence": {"score": 1.0, "band": "verbatim", "quote": "To determine the diagnostic performance of EN1, we performed receiver operating characteristic (ROC) curve analysis, calculating the area under the ROC curve (AUC). Tumor data from TCGA and normal tissue data from GTEx were used."}}], "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": ["group_summary", "welch_t", "bootstrap_mean_difference", "group_comparison_anova"], "missing_operations": [], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: anova."}, "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": ["GSE182109"], "resolved_accessions": ["GSE182109"], "statements": ["Data availabilityThe bulk RNA-seq and corresponding clinical data analyzed in this study were obtained from The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/), including the TCGA-GBM and TCGA-LGG cohorts. Normal brain tissue expression data were retrieved from the Genotype-Tissue Expression (GTEx) project (https://gtexportal.org). An independent validation cohort was obtained from the Chinese Glioma Genome Atlas (CGGA, http://www.cgga.org.cn/), specifically the CGGA_693 dataset. Single-cell RNA sequencing data were downloaded from the Gene Expression Omnibus (GEO) database under accession number GSE182109. All datasets used in this study are publicly available and fully de-identified."], "basis": "The paper names GSE182109 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": ["group_summary", "welch_t", "bootstrap_mean_difference", "group_comparison_anova"], "assets": ["source:pmid:41513764", "data:geo:GSE182109"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE182109"]}, "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:09:50Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 9949, "completion_tokens": 1805, "cached_tokens": 128}}, "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:41513764", "type": "paper", "label": "pmid:41513764"}, {"id": "source:pmid:41513764", "type": "manuscript", "label": "PMID 41513764 full text"}, {"id": "data:geo:GSE182109", "type": "dataset", "label": "GEO GSE182109"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:41513764", "type": "paper", "label": "pmid:41513764"}, {"id": "source:pmid:41513764", "type": "manuscript", "label": "PMID 41513764 full text"}, {"id": "data:geo:GSE182109", "type": "dataset", "label": "GEO GSE182109"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-41513764", "type": "extraction", "used": ["source:pmid:41513764"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE182109", "type": "accession_resolution", "used": ["data:geo:GSE182109"], "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:41513764", "parent_hash": "sha256:0ff53fbd0e6c8cfaffd2bd361df4529ed210924f3368e5f26aa662cf0f9f3c5f", "fork_of": null, "generation": 2}, "gaps": ["asset:data:geo:GSE182109 licence not read", "asset:data:geo:GSE182109 not downloaded", "execution.not_run", "method.reader_named_unmapped", "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": ["CIBERSORT, ESTIMATE, and WGCNA. Single-cell RNA-seq data (GSE182109) were used to characterize t", "seq data processing and analysis Single-cell RNA-seq data (GSE182109) were processed using Seurat", "evasion, by integrating bulk and single-cell transcriptomic analyses together with ceRNA regulatory netw"], "order_basis": "paper_text_order:567", "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": ["ated through CIBERSORT, ESTIMATE, and WGCNA. Single-cell RNA-seq data (GSE182109) were used to characterize t", "thods scRNA-seq data processing and analysis Single-cell RNA-seq data (GSE182109) were processed using Seurat", "Materials and methods scRNA-seq data processing and analysis Single-cell RNA"], "order_basis": "paper_text_order:892", "id": "assay:2"}, {"role": "data_generation", "definition_id": "western_blot", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.93, "evidence": ["ase predictions. Functional assays including Western blot, CCK-8, and Transwell experiments were perfo"], "order_basis": "paper_text_order:1140", "id": "assay:3"}, {"role": "data_generation", "definition_id": "apoptosis_assay", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.93, "evidence": ["s. Functional assays including Western blot, CCK-8, and Transwell experiments were performed to"], "order_basis": "paper_text_order:1154", "id": "assay:4"}, {"role": "data_generation", "definition_id": "transfection", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.99, "evidence": ["echanism contributing to EN1 activation. EN1 knockdown significantly suppressed glioblastoma cell p", "using the Monocle3 package. Cell culture and transfection Normal human astrocyte cell line (HA1800) an", "S and 1% penicillin-streptomycin. Cells were transfected with EN1 shRNA or scrambled control shRNA us", "treptomycin. Cells were transfected with EN1 shRNA or scrambled control shRNA using Lipofectami"], "order_basis": "paper_text_order:1847", "id": "assay:5"}, {"role": "data_generation", "definition_id": "cell_culture", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.77, "evidence": ["ll Biology (Shanghai, China). All cells were cultured in standard DMEM supplemented with 10% FBS and"], "order_basis": "paper_text_order:3183", "id": "assay:6"}, {"role": "analysis", "definition_id": "group_comparison", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": ["group_summary", "welch_t", "bootstrap_mean_difference", "group_comparison_anova"], "support_state": "supported", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 0.8, "evidence": ["ANOVA"], "order_basis": "consumes_earlier_output", "id": "assay:7"}]}