N6-methyladenosine-related non-coding RNAs are potential prognostic and immunotherapeutic responsiveness biomarkers for bladder cancer.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:34950253", "kind": "science-record", "revision": 5, "content_hash": "sha256:39ce8918782f344fddc81a4af54ef35dd2371008b4ddadc775f8ebc86cdcd4c2", "created_at": "2026-11-01T17:36:39Z", "updated_at": "2026-11-01T22:48:07Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "N6-methyladenosine-related non-coding RNAs are potential prognostic and immunotherapeutic responsiveness biomarkers for bladder cancer.", "doi": "10.1007/s13167-021-00259-w", "pmid": "34950253", "authors": ["Miaolong Lu", "Hailun Zhan", "Bolong Liu", "Dongyang Li", "Wenbiao Li", "Xuelian Chen", "Xiangfu Zhou"], "journal": "The EPMA journal", "year": 2021, "abstract": "BACKGROUND: Bladder cancer (BC) is a commonly occurring malignant tumor of the urinary system, demonstrating high global morbidity and mortality rates. BC currently lacks widely accepted biomarkers and its predictive, preventive, and personalized medicine (PPPM) is still unsatisfactory. N6-methyladenosine (m6A) modification and non-coding RNAs (ncRNAs) have been shown to be effective prognostic and immunotherapeutic responsiveness biomarkers and contribute to PPPM for various tumors. However, their role in BC remains unclear.\nMETHODS: m6A-related ncRNAs (lncRNAs and miRNAs) were identified through a comprehensive analysis of TCGA, starBase, and m6A2Target databases. Using TCGA dataset (training set), univariate and least absolute shrinkage and selection operator (LASSO) regression analyses were performed to develop an m6A-related ncRNA-based prognostic risk model. Kaplan-Meier analysis of overall survival (OS) and receiver operating characteristic (ROC) curves were used to verify the prognostic evaluation power of the risk model in the GSE154261 dataset (testing set) from Gene Expression Omnibus (GEO). A nomogram containing independent prognostic factors was developed. Differences in BC clinical characteristics, m6A regulators, m6A-related ncRNAs, gene expression patterns, and differentially expressed genes (DEGs)-associated molecular networks between the high- and low-risk groups in TCGA dataset were also analyzed. Additionally, the potential applicability of the risk model in the prediction of immunotherapeutic responsiveness was evaluated based on the \"IMvigor210CoreBiologies\" data set.\nRESULTS: We identified 183 m6A-related ncRNAs, of which 14 were related to OS. LASSO regression analysis was further used to develop a prognostic risk model that included 10 m6A-related ncRNAs (BAALC-AS1, MIR324, MIR191, MIR25, AC023509.1, AL021707.1, AC026362.1, GATA2-AS1, AC012065.2, and HCP5). The risk model showed an excellent prognostic evaluation performance in both TCGA and GSE154261 datasets, with ROC curve areas under the curve (AUC) of 0.62 and 0.83, respectively. A nomogram containing 3 independent prognostic factors (risk score, age, and clinical stage) was developed and was found to demonstrate high prognostic prediction accuracy (AUC = 0.83). Moreover, the risk model could also predict BC progression. A higher risk score indicated a higher pathological grade and clinical stage. We identified 1058 DEGs between the high- and low-risk groups in TCGA dataset; these DEGs were involved in 3 molecular network systems, i.e., cellular immune response, cell adhesion, and cellular biological metabolism. Furthermore, the expression levels of 8 m6A regulators and 12 m6A-related ncRNAs were significantly different between the two groups. Finally, this risk model could be used to predict immunotherapeutic responses.\nCONCLUSION: Our study is the first to explore the potential application value of m6A-related ncRNAs in BC. The m6A-related ncRNA-based risk model demonstrated excellent performance in predicting prognosis and immunotherapeutic responsiveness. Based on this model, in addition to identifying high-risk patients early to provide them with focused attention and targeted prevention, we can also select beneficiaries of immunotherapy to deliver personalized medical services. Furthermore, the m6A-related ncRNAs could elucidate the molecular mechanisms of BC and lead to a new direction for the improvement of PPPM for BC.\nSUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s13167-021-00259-w.", "claims": [{"id": "claim-method", "statement": "Using TCGA dataset (training set), univariate and least absolute shrinkage and selection operator (LASSO) regression analyses were performed to develop an m6A-related ncRNA–based prognostic risk model. Kaplan-Meier analysis of overall survival (OS) and receiver operating characteristic (ROC) curves were used to verify the prognostic evaluation power of the risk model in the GSE154261 dataset (testing set) from Gene Expression Omnibus (GEO).", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:34950253", "locator": "sections[0] \"Abstract\" paragraph 2", "quote": "Using TCGA dataset (training set), univariate and least absolute shrinkage and selection operator (LASSO) regression analyses were performed to develop an m6A-related ncRNA–based prognostic risk model. Kaplan-Meier analysis of overall survival (OS) and receiver operating characteristic (ROC) curves were used to verify the prognostic evaluation power of the risk model in the GSE154261 dataset (testing set) from Gene Expression Omnibus (GEO)."}], "state": "complete"}, {"id": "claim-data-access", "statement": "Data and download URLs involved in this study had been described in detail in the materials and methods section.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:34950253", "locator": "sections[25] \"Availability of data and material\" paragraph 2", "quote": "Data and download URLs involved in this study had been described in detail in the materials and methods section."}], "state": "complete"}, {"id": "claim-target-result", "statement": "The risk model showed an excellent prognostic evaluation performance in both TCGA and GSE154261 datasets, with ROC curve areas under the curve (AUC) of 0.62 and 0.83, respectively.", "kind": "qualitative", "target": null, "evidence": [{"asset": "source:pmid:34950253", "locator": "sections[0] \"Abstract\" paragraph 3", "quote": "The risk model showed an excellent prognostic evaluation performance in both TCGA and GSE154261 datasets, with ROC curve areas under the curve (AUC) of 0.62 and 0.83, respectively."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "complete", "objective": "The risk model showed an excellent prognostic evaluation performance in both TCGA and GSE154261 datasets, with ROC curve areas under the curve (AUC) of 0.62 and 0.83, respectively.", "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:34950253"], "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:GSE154261"], "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:GSE154261"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": []}, "data": {"state": "partial", "assets": [{"id": "source:pmid:34950253", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1007/s13167-021-00259-w", "sha256": "sha256:c7132b104aa47568e8c9a23b283c4ec1cdc4145f6409a3a1efe7ed0bc0972cd5", "bytes": 57828, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1007/s13167-021-00259-w", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC8648947", "source": "europe_pmc_fulltext", "coverage": "whole_body", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE154261", "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=GSE154261", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE154261", "version": null, "note": "Data and download URLs involved in this study had been described in detail in the materials and methods section.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "Data and download URLs involved in this study had been described in detail in the materials and methods section.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM4667908", "GSM4667909", "GSM4667910", "GSM4667911", "GSM4667912", "GSM4667913", "GSM4667914", "GSM4667915", "GSM4667916", "GSM4667917", "GSM4667918", "GSM4667919", "GSM4667920", "GSM4667921", "GSM4667922", "GSM4667923", "GSM4667924", "GSM4667925", "GSM4667926", "GSM4667927", "GSM4667928", "GSM4667929", "GSM4667930", "GSM4667931", "GSM4667932", "GSM4667933", "GSM4667934", "GSM4667935", "GSM4667936", "GSM4667937", "GSM4667938", "GSM4667939", "GSM4667940", "GSM4667941", "GSM4667942", "GSM4667943", "GSM4667944", "GSM4667945", "GSM4667946", "GSM4667947", "GSM4667948", "GSM4667949", "GSM4667950", "GSM4667951", "GSM4667952", "GSM4667953", "GSM4667954", "GSM4667955", "GSM4667956", "GSM4667957", "GSM4667958", "GSM4667959", "GSM4667960", "GSM4667961", "GSM4667962", "GSM4667963", "GSM4667964", "GSM4667965", "GSM4667966", "GSM4667967", "GSM4667968", "GSM4667969", "GSM4667970", "GSM4667971", "GSM4667972", "GSM4667973", "GSM4667974", "GSM4667975", "GSM4667976", "GSM4667977", "GSM4667978", "GSM4667979", "GSM4667980", "GSM4667981", "GSM4667982", "GSM4667983", "GSM4667984", "GSM4667985", "GSM4667986", "GSM4667987", "GSM4667988", "GSM4667989", "GSM4667990", "GSM4667991", "GSM4667992", "GSM4667993", "GSM4667994", "GSM4667995", "GSM4667996", "GSM4667997", "GSM4667998", "GSM4667999", "GSM4668000", "GSM4668001", "GSM4668002", "GSM4668003", "GSM4668004", "GSM4668005", "GSM4668006"], "sample_count": 99, "group_sizes": null, "group_sizes_basis": null, "platform": "20301", "taxon": "Homo sapiens", "title": "Identification of Differential Tumor Subtypes of T1 Bladder Cancer", "resolution_checked_at": "2026-11-01T17:30:42Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}], "gaps": ["asset:data:geo:GSE154261 not downloaded", "asset:data:geo:GSE154261 licence not read"]}, "analysis": {"state": "complete", "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": "Using TCGA dataset (training set), univariate and least absolute shrinkage and selection operator (LASSO) regression analyses were performed to develop an m6A-related ncRNA–based prognostic risk model. Kaplan-Meier analysis of overall survival (OS) and receiver operating characteristic (ROC) curves were used to verify the prognostic evaluation power of the risk model in the GSE154261 dataset (testing set) from Gene Expression Omnibus (GEO)."}}, {"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.8, "band": "high", "quote": "Cox regression"}}, {"id": "method-3", "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"}}, {"id": "method-4", "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.7, "band": "medium", "quote": "ROC curve"}}], "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": "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": ["GSE154261"], "resolved_accessions": ["GSE154261"], "statements": ["Data and download URLs involved in this study had been described in detail in the materials and methods section."], "basis": "The paper names GSE154261 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": ["survival_km"], "assets": ["source:pmid:34950253", "data:geo:GSE154261"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE154261"]}, "gaps": []}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-11-01T17:13:01Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 16817, "completion_tokens": 2452, "cached_tokens": 264}}, "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:34950253", "type": "paper", "label": "pmid:34950253"}, {"id": "source:pmid:34950253", "type": "manuscript", "label": "PMID 34950253 full text"}, {"id": "data:geo:GSE154261", "type": "dataset", "label": "GEO GSE154261"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:34950253", "type": "paper", "label": "pmid:34950253"}, {"id": "source:pmid:34950253", "type": "manuscript", "label": "PMID 34950253 full text"}, {"id": "data:geo:GSE154261", "type": "dataset", "label": "GEO GSE154261"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-34950253", "type": "extraction", "used": ["source:pmid:34950253"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE154261", "type": "accession_resolution", "used": ["data:geo:GSE154261"], "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:34950253", "parent_hash": "sha256:a9d018163437d55b65bb682f1619a2b3c855ad38dd162a12bdcc067ad8d0bfa1", "fork_of": null, "generation": 4}, "gaps": ["asset:data:geo:GSE154261 licence not read", "asset:data:geo:GSE154261 not downloaded", "execution.not_run", "review.human_critique", "review.independent_reader", "review.machine_critique"], "assays": [{"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": ["Using TCGA dataset (training set), univariate and least absolute shrinkage and selection operator (LASSO) regression analyses were performed to develop an m6A-related ncRNA–based prognostic risk model. Kaplan-Meier analysis of overall survival (OS) and receiver operating characteristic (ROC) curves were used to verify the prognostic evaluation power of the risk model in the GSE154261 dataset (testing set) from Gene Expression Omnibus (GEO)."], "order_basis": "analysis_only", "id": "assay:1"}]}