{"id":"0d400dd7-76aa-4b06-b2c8-0306286fc822","arxiv_id":"2505.06583","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":1.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A pedagogical review of persistent homology with an illustrative Betti number computation (β0=1, β1=3, β2=0) on a supercoiled DNA structure.","lead":"This paper is a tutorial that explains persistent homology and applies it to the 3D structure of a supercoiled DNA molecule. It is aimed at readers new to algebraic topology who want a plain-language entry into topological data analysis and its biological uses.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central application's input data cannot be a DNA molecule: the paper downloads a '3-1m-supercoiled DNA' structure from AlphaFold and extracts Cα atoms, but AlphaFold contains only protein structures, so the reported Betti numbers are likely computed on a protein mislabeled as DNA.","rationale":"The reader's weakest_assumption correctly identifies the protein/DNA conflation as the central threat to the application's claim. My independent reading of Section 5.1 confirms this: the authors say they take Cα coordinates from AlphaFold for a 'DNA' molecule, which is structurally impossible. This is not merely a labeling slip; if the input point cloud is not DNA, then the computed Betti numbers say nothing about the 3-1m-supercoiled DNA structure that the abstract and title advertise. I also note a second, independent inconsistency: the authors describe H1 features as short-lived/minor in the barcode discussion, then interpret β1=3 as 'three significant loops.' This undermines the strength of the interpretation even if the data-source issue were fixed. Both issues are correctable, and the pedagogical introduction to PH itself is largely standard and useful, so a conditional verdict remains appropriate: accept only after the authors either supply a valid DNA data source (with accession) and reconcile the significance claim, or reposition the application as a protein example. Since the reader's verdict already imposes this condition, I recommend no change to the verdict outcome, though I would emphasize the data-source problem as the primary blocker.","tokens_in":14073,"tokens_out":3922,"duration_ms":39021,"concrete_test":"Request the exact AlphaFold accession or PDB ID used in Section 5.1 and verify whether the entry exists and corresponds to a DNA molecule. Since AlphaFold contains only protein predictions, if the entry is a protein or does not exist, the reported Betti numbers cannot be attributed to 3-1m-supercoiled DNA. As a complementary check, re-run the described pipeline on a genuine supercoiled DNA structure (e.g., a PDB entry for supercoiled DNA) and see whether the Betti numbers match Table 3.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 5.1 states that the 3D structure of the '3-1m-supercoiled DNA molecules sequence' is downloaded from the AlphaFold database and that the xyz coordinates of carbon-alpha (Cα) atoms are extracted. This conflates two incompatible objects: AlphaFold is a protein structure database, and Cα atoms are protein backbone atoms; DNA has no amino acids and no Cα atoms. The only way the reported 111-point point cloud can exist is if the authors actually downloaded a protein structure and mislabeled it as DNA. If so, the central claim—that β0=1, β1=3, β2=0 describe the topology of 3-1m-supercoiled DNA, with β1=3 revealing three significant loops—is unsupported: the computation may be internally correct for some protein, but it does not describe the claimed biological object. Additionally, the paper's own barcode interpretation (Section 5.1, Figure 15) says H1 features have 'relatively short lifespans, suggesting minor or less significant loops,' which directly contradicts the later claim that β1=3 reveals 'three significant loops or cycles.' Thus both the data-source validity and the significance interpretation are load-bearing; either one, if wrong, invalidates the application's conclusion.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a pedagogical introduction to persistent homology, moving from topological spaces, homeomorphisms, homotopy, and metric spaces through simplicial complexes, homology, and persistence diagrams/barcodes, and then applies the pipeline to a structure labeled '3-1m-supercoiled DNA'. The application extracts xyz coordinates of carbon-alpha atoms from AlphaFold, computes a Vietoris–Rips filtration with GUDHI/Ripser, and reports Betti numbers β0=1, β1=3, β2=0, interpreting β1=3 as three significant loops or cycles.","tokens_in":14435,"tokens_out":3932,"duration_ms":39135,"significance":"The pedagogical sections (Sections 2–4) contain mostly standard, correctly stated definitions and will likely be useful to readers without an algebraic-topology background. The paper also provides concrete implementation pointers and a worked computational example, which increase its value as an introductory tutorial. However, the application section contains a load-bearing inconsistency about the biological object being analyzed, and the interpretation of the H1 features is internally contradictory. Because the stated application is the paper's main demonstration, these issues need to be resolved before the manuscript can be accepted.","major_comments":[{"comment":"The application's input object is described inconsistently. Section 5 opens by stating that the method is applied to a protein, but §5.1 identifies the same structure as '3-1m-supercoiled DNA molecules' and says its 3D structure was downloaded from the AlphaFold database using carbon-alpha (Cα) atom coordinates. AlphaFold contains protein structures, and Cα atoms are protein backbone atoms; DNA has no amino acids and no Cα atoms. As written, the reported Betti numbers in Table 3 cannot be claimed to describe supercoiled DNA. The authors must either correct the biological object and data source or explicitly re-frame the example as a protein-structure analysis.","section":"§5, Step 1 and §5.1"},{"comment":"The interpretation of the H1 features is internally contradictory. The text accompanying Figure 15 says that H1 features have 'relatively short lifespans, suggesting minor or less significant loops,' while the discussion after Table 3 says that β1=3 reveals 'the presence of three significant loops or cycles.' Since the significance of loops is the central biological conclusion of the application, these two statements must be reconciled with a concrete persistence criterion.","section":"§5.1, Figure 15 and Table 3"},{"comment":"The computation is not reproducible from the manuscript as written. No filtration parameters, persistence pairs, or code output are given—only the final Betti numbers. The claim that the three H1 features are 'significant' therefore rests on an unspecified threshold. Please include the persistence diagram coordinates or a stated persistence threshold that justifies labeling β1=3 as three significant cycles.","section":"§5.1"}],"minor_comments":[{"comment":"The terms 'protein' and 'DNA' are used interchangeably; the authors should choose a consistent terminology for the analyzed object.","section":"§5"},{"comment":"The heading 'Vietoris-Risp filtration' should read 'Vietoris–Rips filtration'.","section":"§5.1"},{"comment":"The caption contains the typo 'residus'; it should be 'residues'.","section":"Figure 12"},{"comment":"The boundary maps are denoted δ_k earlier, but Definition 11 uses d_k; the notation should be unified.","section":"Definition 11"}],"recommendation":"major_revision","confidential_remarks":"The pedagogical core of the manuscript is sound, but the application section's data-source inconsistency is substantial. If the authors cannot verify the biological identity of the structure they analyzed, the application should be removed or replaced with a verifiable example. The internal contradiction between the persistence-barcode interpretation and the Betti-number discussion also needs explicit resolution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Sarah — quick take on arXiv:2505.06583. It's a teaching-oriented introduction to persistent homology, not a research contribution. The math in Sections 2–4 is standard and largely correct: definitions of topological spaces, simplicial complexes, chain complexes, homology and Betti numbers are set out cleanly, with worked examples and barcode/persistence-diagram explanations. I checked the boundary-map example and the chain-complex definitions; no errors there. As a primer for someone with no algebraic topology background, it does its job. The exposition is the strongest part.\n\nThe application section is where things fall apart. The paper says it's analyzing '3-1m-supercoiled DNA' but the method extracts Cα carbon-alpha coordinates from AlphaFold. AlphaFold contains protein structures only, and Cα atoms are protein backbone atoms; DNA has no amino acids. So the 111-point point cloud has to be a protein molecule, mislabeled as DNA. That is a factual error, and it makes the reported Betti numbers vacuous for the claimed object. The internal contradiction the stress-test flagged is real too: the barcode discussion in Figure 15 says H1 features have 'relatively short lifespans, suggesting minor or less significant loops,' but Table 3 and the text call β1=3 'three significant loops or cycles.' You cannot have both. Also, the GitHub link is just 'this GitHub repository' with no URL, so the computation is not reproducible as written.\n\nNone of this kills the pedagogical value. The math sections would be a fine basis for a tutorial article after the application is either fixed (use a real protein, or use a DNA structure from PDB with C1'/P atoms and correct the interpretation) or removed/rewritten as a generic protein example. The contradiction in significance should be resolved; it's not just wording, it's the difference between noise and structure. I'd like to see the code link and the AlphaFold identifier.\n\nMy call: this deserves a serious referee, on the understanding that the referee's main job will be to require the application section to be made consistent and reproducible. It is not a desk reject. But it is not a paper to cite for any biological claim. In the reading group it could be useful if we're ever teaching TDA to non-specialists.","headline":"A clear, correct pedagogical survey of persistent homology whose biological application is compromised by a protein/DNA mix-up and a contradictory interpretation of H1 significance; fixable, and worth fixing.","tokens_in":14838,"tokens_out":2007,"would_cite":false,"duration_ms":19362,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["55N31","55-01","92-08"],"pacs":[],"model":"deepseek-v4-flash","headline":"Persistent homology, applied to a 3-1 supercoiled DNA structure, reduces its shape to Betti numbers β0=1, β1=3, β2=0 — one connected component, three significant loops, and no voids.","keywords":["persistent homology","topological data analysis","Betti numbers","supercoiled DNA","Vietoris–Rips filtration","persistence barcode","pedagogical introduction","protein structure"],"falsifier":"Fetch the structure used in the example and inspect its database entry: if the sequence is a protein rather than the named 3-1m-supercoiled DNA, the point cloud is not a faithful DNA representation. Recomputing the Vietoris–Rips persistence diagram on a correctly labelled DNA double-helix model would settle whether β1=3 is a genuine property of supercoiling or an artefact of the coordinate source.","tokens_in":13861,"feed_emoji":"🧬","tokens_out":7572,"duration_ms":68342,"temperature":0.7,"pith_summary":"This paper sets out to make persistent homology accessible to readers with no prior algebraic topology, and to demonstrate it on a real biological object. It builds the pipeline from first principles — topological spaces, simplicial complexes, chain groups, homology, Betti numbers — and then applies it to a 3-1 supercoiled DNA structure. The worked example reports Betti numbers β0=1, β1=3, β2=0, which the authors read as one connected component, three significant loops or cycles, and no enclosed voids. The broader claim is that this multi-scale topological summary captures structure that ordinary machine learning on coordinates would miss, and that the method is learnable by nonspecialists.","feed_headline":"Topology finds 3 loops in a supercoiled DNA structure","feed_subtitle":"One connected component, three persistent cycles, zero voids: a step-by-step tutorial turns coordinates into interpretable Betti numbers.","key_machinery":"The carrying object is the Vietoris–Rips filtration built from pairwise Euclidean distances among the coordinate points. As the scale radius grows, simplices are added and topological features are born and die; homology in each dimension counts cycles modulo boundaries. The output is a persistence diagram and a barcode, whose long bars designate the meaningful features, and the final numeric summary is the set of Betti numbers. The structure theorem for persistence modules justifies encoding persistent homology as a multiset of intervals, and the stability theorem supports treating long-lived features as reliable.","core_discovery":"The paper's central claim is that a 3-1m-supercoiled DNA molecule, represented as a point cloud of alpha-carbon coordinates and filtered with a Vietoris–Rips construction, carries the persistent homology signature β0=1, β1=3, β2=0. The intended reading is that the molecule is a single connected piece with three long-lived one-dimensional cycles — twists or knot-like loops that may be biologically meaningful — and no cavities. Around this example the paper also claims that the full persistent-homology pipeline, from pairwise distances to barcodes to Betti numbers, is usable by people with no prior topological training.","pith_inferences":["My inference: if the coordinate set is valid but the molecule is actually a protein rather than the named DNA, then the reported 'DNA' signature is in fact a protein-folding signature; re-labelling the object preserves the computation but changes the biological interpretation.","My inference: the three one-dimensional cycles likely correspond to superhelical turns or plectonemic crossings; comparing persistence diagrams of relaxed versus supercoiled forms of the same sequence would test whether β1 tracks supercoiling density.","My inference: the paper's own tutorial pipeline suggests a natural extension it does not run — vectorising the persistence diagram into a persistence landscape or persistence image and feeding it to a classifier to distinguish topoisomers automatically."],"forward_implications":["Persistent homology reduces a complex three-dimensional biomolecular geometry to a small set of interpretable numbers that can serve as features in clustering or classification pipelines.","Because long-lived topological features are stable under small perturbations of the input coordinates, the three reported loops are more trustworthy than short-lived fluctuations.","The same step-by-step recipe — point cloud, filtration, barcode, Betti numbers — transfers directly to other biological structures such as proteins, RNA, or molecular dynamics ensembles.","For nonspecialists, the tutorial format lowers the barrier to applying topological data analysis without first taking a course in algebraic topology.","The computed signature (one component, three loops, zero voids) gives a concrete topological fingerprint of the supercoiled state that could be compared across molecular conformations."],"supporting_citations":[{"why":"supplies the 3D coordinate data that become the point cloud in the DNA example.","marker":"[29]"},{"why":"sets out the filtration constructions, including Vietoris–Rips complexes, used to build the sequence of simplicial complexes.","marker":"[23]"},{"why":"provides the barcode and persistence-diagram framework used for visualising feature birth and death.","marker":"[8]"},{"why":"underwrites the structure theorem that justifies summarising persistent homology as a multiset of intervals.","marker":"[6]"},{"why":"supports the claim that small perturbations change persistence diagrams only slightly, making long-lived features reliable.","marker":"[7]"},{"why":"gives the protein-folding application context the paper extends to a supercoiled structure.","marker":"[18]"}],"fun_headline_variants":["Persistent homology reveals 3 loops in supercoiled DNA","Three cycles, no voids: topology of supercoiled DNA","Betti numbers for DNA: one component, three loops, zero voids","From coordinates to Betti numbers: a DNA tutorial","Supercoiled DNA: β1=3, a persistent homology result"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole application rests on treating the downloaded 3D coordinates as a valid protein-style model of the supercoiled DNA molecule; if those coordinates are mislabelled or the alpha-carbon convention does not apply to DNA, the Betti numbers describe a different object.","fun_headline_variants_meta":{"raw":{"variants":["Persistent homology reveals 3 loops in supercoiled DNA","Three cycles, no voids: topology of supercoiled DNA","Betti numbers for DNA: one component, three loops, zero voids","From coordinates to Betti numbers: a DNA tutorial","Supercoiled DNA: β1=3, a persistent homology result"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000913,"raw_usage":{"total_tokens":3856,"prompt_tokens":814,"completion_tokens":3042,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":430,"completion_tokens_details":{"reasoning_tokens":2954}},"tokens_in":430,"tokens_out":3042,"duration_ms":22146,"temperature":1.0,"reasoning_tokens":2954,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:37:27.120192+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fetch the structure used in the example and inspect its database entry: if the sequence is a protein rather than the named 3-1m-supercoiled DNA, the point cloud is not a faithful DNA representation. Recomputing the Vietoris–Rips persistence diagram on a correctly labelled DNA double-helix model would settle whether β1=3 is a genuine property of supercoiling or an artefact of the coordinate source.","supporting_citations":[{"cited_title":"Springer, 2016","cited_arxiv_id":null,"evidence_quote":"underwrites the structure theorem that justifies summarising persistent homology as a multiset of intervals."},{"cited_title":"Stability of persistence diagrams.Discrete & Computational Geometry, 37(1):103–120, 2007","cited_arxiv_id":null,"evidence_quote":"supports the claim that small perturbations change persistence diagrams only slightly, making long-lived features reliable."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the protein-folding application context the paper extends to a supercoiled structure."}],"review_version":1}