Ψ-co-mAFiA: concurrent detection of pseudouridine and m6A in single RNA molecules.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:41002271", "kind": "science-record", "revision": 2, "content_hash": "sha256:60898ea234645d9ac7f75e04440f352efd2d9e03d142d4c6e7c0ed5f4ae14fc7", "created_at": "2026-11-01T20:31:04Z", "updated_at": "2026-11-01T21:35:48Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Ψ-co-mAFiA: concurrent detection of pseudouridine and m6A in single RNA molecules.", "doi": "10.1093/bioinformatics/btaf536", "pmid": "41002271", "authors": ["Adrian Chan", "Isabel S Naarmann-de Vries", "Christoph Dieterich"], "journal": "Bioinformatics (Oxford, England)", "year": 2025, "abstract": "SUMMARY: The development of third-generation sequencing technologies enables the detection of RNA modifications at single-molecule resolution. Specifically for direct RNA sequencing on the ONT platform, we have previously developed an m6A detection algorithm called mAFiA. Here, we present the updated method, now covering all 18 DRACH m6A contexts as well as the identification of pseudouridine sites (Ψ). Our modification level predictions compare favorably with orthogonal methods and respond to knockdown or knock out of writer proteins. The simultaneous detection of multiple modifications on a single RNA molecule opens up the possibility to study cross-modification interactions.\nAVAILABILITY AND IMPLEMENTATION: Ψ-co-mAFiA is available at https://github.com/dieterich-lab/psi-co-mAFiA and licensed under GPLv3.0. An archived version of the software is available on Zenodo at https://doi.org/10.5281/zenodo.16797676.", "claims": [{"id": "claim-method", "statement": "The overall levels of the two modifications in a molecule are positively correlated (Pearson correlation ρ=0.986, P-value = 0.002).", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:41002271", "locator": "sections[3] \"2.1 Algorithm\" paragraph 3", "quote": "The overall levels of the two modifications in a molecule are positively correlated (Pearson correlation ρ=0.986, P-value = 0.002)."}], "state": "complete"}, {"id": "claim-data-access", "statement": "Source code and documentation are available under GPLv3.0 at https://github.com/dieterich-lab/psi-co-mAFiA. An archived version of the software, which was used to produce all results, is available under https://doi.org/10.5281/zenodo.16797676.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:41002271", "locator": "sections[12] \"Data availability\" paragraph 2", "quote": "Source code and documentation are available under GPLv3.0 at https://github.com/dieterich-lab/psi-co-mAFiA. An archived version of the software, which was used to produce all results, is available under https://doi.org/10.5281/zenodo.16797676."}], "state": "complete"}, {"id": "claim-target-result", "statement": "The overall levels of m6A and Ψ within individual RNA molecules are positively correlated, although highly saturated modified reads show an anti-correlation.", "kind": "quantitative", "target": {"value": 0.002, "unit": "p-value", "tolerance": null, "source": "The overall levels of the two modifications in a molecule are positively correlated (Pearson correlation ρ=0.986, P-value = 0.002)."}, "evidence": [{"asset": "source:pmid:41002271", "locator": "sections[3] \"2.1 Algorithm\" paragraph 3", "quote": "The overall levels of the two modifications in a molecule are positively correlated (Pearson correlation ρ=0.986, P-value = 0.002)."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "partial", "objective": "The overall levels of m6A and Ψ within individual RNA molecules are positively correlated, although highly saturated modified reads show an anti-correlation.", "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:41002271"], "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:41002271", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1093/bioinformatics/btaf536", "sha256": "sha256:4440e1b228d5428536459552beb5e1ac09871aa920b46b6a38c6b1f6dd043eec", "bytes": 15152, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1093/bioinformatics/btaf536", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC12552086", "source": "europe_pmc_fulltext", "coverage": "whole_body", "text_hash_scope": "rendered_section_body"}}], "gaps": []}, "analysis": {"state": "complete", "methods": [{"id": "method-1", "family": "correlation_analysis", "specific": "pearson", "label": "Correlation analysis", "dsl_operations": ["pearson_r"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["paired_numeric_columns"], "evidence": {"score": 1.0, "band": "verbatim", "quote": "The overall levels of the two modifications in a molecule are positively correlated (Pearson correlation ρ=0.986, P-value = 0.002)."}}], "targets": [{"id": "target-published-estimate", "claim_id": "claim-target-result", "operation": "pearson_r", "expected_value": 0.002, "unit": "p-value", "tolerance": null, "state": "complete", "note": "The threshold the paper applied to its own result, as the paper states it. A run that reproduces the effect is judged against this threshold."}], "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": ["pearson_r"], "missing_operations": [], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: pearson."}, "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": ["Source code and documentation are available under GPLv3.0 at https://github.com/dieterich-lab/psi-co-mAFiA. An archived version of the software, which was used to produce all results, is available under https://doi.org/10.5281/zenodo.16797676."], "basis": "The paper names no accession to its source data, so the data has to come from the authors."}}, "can_run": {"dsl_operations": ["pearson_r"], "assets": ["source:pmid:41002271"], "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-11-01T20:02:01Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 5466, "completion_tokens": 1363, "cached_tokens": 256}}, "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:41002271", "type": "paper", "label": "pmid:41002271"}, {"id": "source:pmid:41002271", "type": "manuscript", "label": "PMID 41002271 full text"}, {"id": "data:geo:None", "type": "dataset", "label": "GEO None"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:41002271", "type": "paper", "label": "pmid:41002271"}, {"id": "source:pmid:41002271", "type": "manuscript", "label": "PMID 41002271 full text"}, {"id": "data:geo:None", "type": "dataset", "label": "GEO None"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-41002271", "type": "extraction", "used": ["source:pmid:41002271"], "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:41002271", "parent_hash": "sha256:18f21c6156b06dcc2813ebc7d3a0d66883b0440c0b575a724db5422763d65c1d", "fork_of": null, "generation": 1}, "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": "rna_seq", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.93, "evidence": ["molecule resolution. Specifically for direct RNA sequencing on the ONT platform, we have previously deve"], "order_basis": "paper_text_order:205", "id": "assay:1"}, {"role": "data_generation", "definition_id": "transfection", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.85, "evidence": ["rably with orthogonal methods and respond to knockdown or knock out of writer proteins. The simulta"], "order_basis": "paper_text_order:537", "id": "assay:2"}, {"role": "analysis", "definition_id": "correlation_analysis", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": ["pearson_r"], "support_state": "supported", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 1.0, "evidence": ["The overall levels of the two modifications in a molecule are positively correlated (Pearson correlation ρ=0.986, P-value = 0.002)."], "order_basis": "consumes_earlier_output", "id": "assay:3"}]}