A Bayesian hierarchical model to estimate DNA methylation conservation in colorectal tumors.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:34487148", "kind": "science-record", "revision": 3, "content_hash": "sha256:05129ce188e8884e7d94964dcaa7d697ce90d03f1588dbb2181bfa117627a3bf", "created_at": "2026-11-01T17:36:03Z", "updated_at": "2026-11-01T21:26:09Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "A Bayesian hierarchical model to estimate DNA methylation conservation in colorectal tumors.", "doi": "10.1093/bioinformatics/btab637", "pmid": "34487148", "authors": ["Kevin A Murgas", "Yanlin Ma", "Lidea K Shahidi", "Sayan Mukherjee", "Andrew S Allen", "Darryl Shibata", "Marc D Ryser"], "journal": "Bioinformatics (Oxford, England)", "year": 2021, "abstract": "MOTIVATION: Conservation is broadly used to identify biologically important (epi)genomic regions. In the case of tumor growth, preferential conservation of DNA methylation can be used to identify areas of particular functional importance to the tumor. However, reliable assessment of methylation conservation based on multiple tissue samples per patient requires the decomposition of methylation variation at multiple levels.\nRESULTS: We developed a Bayesian hierarchical model that allows for variance decomposition of methylation on three levels: between-patient normal tissue variation, between-patient tumor-effect variation and within-patient tumor variation. We then defined a model-based conservation score to identify loci of reduced within-tumor methylation variation relative to between-patient variation. We fit the model to multi-sample methylation array data from 21 colorectal cancer (CRC) patients using a Monte Carlo Markov Chain algorithm (Stan). Sets of genes implicated in CRC tumorigenesis exhibited preferential conservation, demonstrating the model's ability to identify functionally relevant genes based on methylation conservation. A pathway analysis of preferentially conserved genes implicated several CRC relevant pathways and pathways related to neoantigen presentation and immune evasion. Our findings suggest that preferential methylation conservation may be used to identify novel gene targets that are not consistently mutated in CRC. The flexible structure makes the model amenable to the analysis of more complex multi-sample data structures.\nAVAILABILITY AND IMPLEMENTATION: The data underlying this article are available in the NCBI GEO Database, under accession code GSE166212. The R analysis code is available at https://github.com/kevin-murgas/DNAmethylation-hierarchicalmodel.\nSUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.", "claims": [{"id": "claim-method", "statement": "We developed a Bayesian hierarchical model that allows for variance decomposition of methylation on three levels: between-patient normal tissue variation, between-patient tumor-effect variation and within-patient tumor variation.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:34487148", "locator": "sections[0] \"Abstract\" paragraph 2", "quote": "We developed a Bayesian hierarchical model that allows for variance decomposition of methylation on three levels: between-patient normal tissue variation, between-patient tumor-effect variation and within-patient tumor variation."}], "state": "complete"}, {"id": "claim-data-access", "statement": "The data underlying this article are available in the NCBI GEO Database, under accession code GSE166212. The R analysis code is available at https://github.com/kevin-murgas/DNAmethylation-hierarchicalmodel.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:34487148", "locator": "sections[0] \"Abstract\" paragraph 3", "quote": "The data underlying this article are available in the NCBI GEO Database, under accession code GSE166212. The R analysis code is available at https://github.com/kevin-murgas/DNAmethylation-hierarchicalmodel."}], "state": "complete"}, {"id": "claim-target-result", "statement": "Sets of genes implicated in CRC tumorigenesis exhibited preferential conservation, demonstrating the model’s ability to identify functionally relevant genes based on methylation conservation.", "kind": "qualitative", "target": null, "evidence": [{"asset": "source:pmid:34487148", "locator": "sections[0] \"Abstract\" paragraph 2", "quote": "Sets of genes implicated in CRC tumorigenesis exhibited preferential conservation, demonstrating the model’s ability to identify functionally relevant genes based on methylation conservation."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "complete", "objective": "Sets of genes implicated in CRC tumorigenesis exhibited preferential conservation, demonstrating the model’s ability to identify functionally relevant genes based on methylation conservation.", "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:34487148"], "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:GSE166212"], "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:GSE166212"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": []}, "data": {"state": "partial", "assets": [{"id": "source:pmid:34487148", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1093/bioinformatics/btab637", "sha256": "sha256:7b29f2d6f0c30061f1e70df87a97d3090a6af10cd37ca647d5e501e4bc339583", "bytes": 1834, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1093/bioinformatics/btab637", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC10060692", "source": "ncbi_pmc_fulltext", "coverage": "whole_body", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE166212", "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=GSE166212", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE166212", "version": null, "note": "The data underlying this article are available in the NCBI GEO Database, under accession code GSE166212. The R analysis code is available at https://github.com/kevin-murgas/DNAmethylation-hierarchicalmodel.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The data underlying this article are available in the NCBI GEO Database, under accession code GSE166212. The R analysis code is available at https://github.com/kevin-murgas/DNAmethylation-hierarchicalmodel.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM5065983", "GSM5065984", "GSM5065985", "GSM5065986", "GSM5065987", "GSM5065988", "GSM5065989", "GSM5065990", "GSM5065991", "GSM5065992", "GSM5065993", "GSM5065994", "GSM5065995", "GSM5065996", "GSM5065997", "GSM5065998", "GSM5065999", "GSM5066000", "GSM5066001", "GSM5066002", "GSM5066003", "GSM5066004", "GSM5066005", "GSM5066006", "GSM5066007", "GSM5066008", "GSM5066009", "GSM5066010", "GSM5066011", "GSM5066012", "GSM5066013", "GSM5066014", "GSM5066015", "GSM5066016", "GSM5066017", "GSM5066018", "GSM5066019", "GSM5066020", "GSM5066021", "GSM5066022", "GSM5066023", "GSM5066024", "GSM5066025", "GSM5066026", "GSM5066027", "GSM5066028", "GSM5066029", "GSM5066030"], "sample_count": 48, "group_sizes": null, "group_sizes_basis": null, "platform": "21145", "taxon": "Homo sapiens", "title": "Colorectal Cancer DNA Methylation", "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:GSE166212 not downloaded", "asset:data:geo:GSE166212 licence not read"]}, "analysis": {"state": "partial", "methods": [{"id": "method-1", "family": null, "specific": null, "label": "Bayesian hierarchical model", "dsl_operations": [], "missing_operations": [], "expressible": null, "parameters": {}, "data_shapes": [], "evidence": {"score": 1.0, "band": "verbatim", "quote": "We developed a Bayesian hierarchical model that allows for variance decomposition of methylation on three levels: between-patient normal tissue variation, between-patient tumor-effect variation and within-patient tumor variation."}}], "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": "blocked", "tier": "R4", "blocking_axis": "method", "blocked_on": "method:not_determinable", "axes": {"method": {"tier": "R4", "expressible": false, "required_operations": [], "missing_operations": [], "basis": "The reader quoted the method the paper names (Bayesian hierarchical model), but the corpus's method extractor could not map it onto a DSL operation, so the method is published but not expressible."}, "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": ["GSE166212"], "resolved_accessions": ["GSE166212"], "statements": ["The data underlying this article are available in the NCBI GEO Database, under accession code GSE166212. The R analysis code is available at https://github.com/kevin-murgas/DNAmethylation-hierarchicalmodel."], "basis": "The paper names GSE166212 in NCBI GEO Database and the accession resolves in the repository. No dataset was downloaded and no licence was read."}}, "can_run": {"dsl_operations": [], "assets": ["source:pmid:34487148", "data:geo:GSE166212"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE166212"]}, "gaps": ["method:not_determinable"]}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-11-01T17:02:44Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 1938, "completion_tokens": 925, "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:34487148", "type": "paper", "label": "pmid:34487148"}, {"id": "source:pmid:34487148", "type": "manuscript", "label": "PMID 34487148 full text"}, {"id": "data:geo:GSE166212", "type": "dataset", "label": "GEO GSE166212"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:34487148", "type": "paper", "label": "pmid:34487148"}, {"id": "source:pmid:34487148", "type": "manuscript", "label": "PMID 34487148 full text"}, {"id": "data:geo:GSE166212", "type": "dataset", "label": "GEO GSE166212"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-34487148", "type": "extraction", "used": ["source:pmid:34487148"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE166212", "type": "accession_resolution", "used": ["data:geo:GSE166212"], "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:34487148", "parent_hash": "sha256:810c546632748896733cb52b0358625a654a11e4ebb349de552d0fa80fd43663", "fork_of": null, "generation": 2}, "gaps": ["asset:data:geo:GSE166212 licence not read", "asset:data:geo:GSE166212 not downloaded", "execution.not_run", "method.reader_named_unmapped", "review.human_critique", "review.independent_reader", "review.machine_critique"], "assays": [{"role": "data_generation", "definition_id": "methylation_array", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.93, "evidence": ["variation. We fit the model to multi-sample methylation array data from 21 colorectal cancer (CRC) patient"], "order_basis": "paper_text_order:898", "id": "assay:1"}]}