{
 "dataset": {
  "archive": "dandi",
  "id": "001692",
  "label": "DANDI:001692"
 },
 "page": "datasets/dandi/001692/",
 "generated": "2026-10-20T19:09:16+00:00",
 "caution": "A failure to reproduce is not a finding that the paper is wrong. It may be a failure of the attempt.",
 "outcome_counts": {
  "reproduced": 18,
  "partially_reproduced": 8,
  "not_reproduced": 3,
  "blocked": 0
 },
 "outcome_counts_with_outside_data": {
  "reproduced": 0,
  "partially_reproduced": 0,
  "not_reproduced": 0,
  "blocked": 0
 },
 "outcome_counts_elsewhere": {
  "reproduced": 3,
  "partially_reproduced": 0,
  "not_reproduced": 0,
  "blocked": 1
 },
 "outcome_counts_links": {
  "reproduced": 0,
  "partially_reproduced": 0,
  "not_reproduced": 0,
  "blocked": 0
 },
 "attempts": [
  {
   "id": "001692-bouabid-2026-01",
   "dataset": {
    "archive": "dandi",
    "id": "001692",
    "version": "0.260303.1724",
    "doi": "10.48324/dandi.001692/0.260303.1724"
   },
   "key": "dandi:001692",
   "label": "DANDI:001692",
   "dataset_version": "0.260303.1724",
   "paper": "10.64898/2026.01.20.700614",
   "paper_title": "An anatomical hotspot for striatal dopamine-acetylcholine interactions during reward and movement",
   "date": "2026-10-19",
   "kind": "agent",
   "accountable": "CatalystNeuro",
   "agent": {
    "harness": "Claude Code (GitHub Actions)",
    "model": "claude-opus-5-5",
    "protocol": "pilot-7"
   },
   "runner": {
    "platform": "GitHub Actions",
    "url": "https://github.com/catalystneuro/reproduction-registry/actions/runs/39187568966"
   },
   "tier": "agent_verified",
   "deprecated": "",
   "deprecated_on": null,
   "entry": "001692-bouabid-2026",
   "entry_title": "Striatal dopamine-acetylcholine hotspot (Bouabid, Vu et al. 2026)",
   "description": "Tests the findings of Bouabid, Vu et al. (2026) on dandiset 001692, which holds multi-fiber photometry of dopamine and acetylcholine across the striatum of 17 mice. It reproduces the spontaneous cross-correlations and their anterior dorsolateral hotspot, paired transients, reward, Pavlovian cue, locomotion and optogenetic stimulation analyses, and the control mice, and tests the ex vivo slice statistics of Figure 8 on the paper's source data on Zenodo.",
   "notebooks": [
    {
     "path": "reproduce_bouabid_2026.ipynb",
     "title": "Reproducing Bouabid, Vu et al. (2026) from dandiset 001692",
     "executed": "notebooks/001692-bouabid-2026/reproduce_bouabid_2026.html"
    }
   ],
   "outcomes": [
    {
     "claim": "fig1e-dominant-polarity",
     "statement": "Across 174 sites in 5 mice, the dominant DA-ACh cross-correlation (the largest magnitude Pearson r within 1 s lags) during spontaneous behavior is negative at 132 sites, 127 of them with DA leading, and positive at 42 sites, 31 of them with ACh leading.",
     "source": "Results, Spatially organized coupling of transient DA-ACh dynamics during behavior, paragraph 1; Figure 1e",
     "quote": "The polarity of the strongest correlation varied significantly across striatal sites, with some exhibiting dominant positive correlations (42/174 sites, 31 with ACh leading) and others dominant negative (132/174 sites, 127 with DA leading)",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "Of 174 sites, 133 have a dominant negative correlation (128 with DA leading) and 41 a dominant positive one (32 with ACh leading); 6 have their extreme at lag 0. All four counts are within 10% of the paper's. Every mouse has more negative than positive sites. Correlation over concatenated sessions gives 133, 128, 41 and 30.",
     "deviations": [
      "One UG27 reward session (240313r) lacks the artifact-removed green series and is skipped, as the code skips empty series; it is not a spontaneous session, so this claim is unaffected."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig1f-adls-hotspot",
     "statement": "The dominant correlations form a significant hotspot, identified by local Moran's I, in the anterior dorsolateral striatum (aDLS), where the strongest correlations are negative.",
     "source": "Results, Spatially organized coupling of transient DA-ACh dynamics during behavior, paragraph 1; Figure 1f",
     "quote": "This analysis identified a significant 'hotspot' in the aDLS where the dominant correlations were highest, and the polarity of the strongest correlations were negative",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "partially_reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "The hotspot procedure finds a significant hotspot of 1.01 mm3 with mean smoothed r -0.368 (hull mean -0.185). Its centroid (AP 0.63, ML 2.02, DV -3.14 mm) is anterior and dorsal of the median site (0.10, 2.04, -3.21 mm) but 0.02 mm medial of it, which fails the lateral part of the criterion fixed beforehand. With the 99th percentile of the null the hotspot has 518 voxels and its centroid is at ML 2.09 mm, lateral of the median. 28 sites from all five mice are inside. The location test against the median site is a coarse stand-in for the paper's maps, in which the hotspot lies in aDLS; the shortfall is small. Finding a hotspot is not evidence on its own, as the procedure finds one in 67 of 100 shuffles of the values across sites; none of these is as large (at most 482 voxels) or has a mean as low.",
     "deviations": [
      "The hotspot procedure is simplified from the paper's code (linear in place of natural-neighbor interpolation, 100 um voxels in place of 50 um, the convex hull of the sites in place of the Allen striatum mask, 200 in place of 10000 shuffles).",
      "The paper does not define the aDLS location in numbers; the criterion compares the centroid with the median site coordinates."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig1h-in-vs-out",
     "statement": "Sites inside the aDLS hotspot have more negative dominant correlations than sites outside (Wilcoxon rank-sum p = 1.67e-16).",
     "source": "Figure 1h",
     "quote": "Sites inside the hotspot show significantly stronger negative correlations",
     "statistics": [
      {
       "test": "Wilcoxon rank-sum",
       "p": 1.67e-16
      }
     ],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "The 28 sites inside have median dominant r -0.41 against -0.14 for the 146 outside, rank-sum p = 1.9e-13 (the paper 1.67e-16). The test is partly circular, as in the paper.",
     "deviations": [
      "The hotspot is the simplified one of fig1f-adls-hotspot."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig1f-positive-posterior",
     "statement": "Dominant positive correlations lie in the posterior striatum.",
     "source": "Results, Spatially organized coupling of transient DA-ACh dynamics during behavior, paragraph 1; Figure 1f",
     "quote": "These maps revealed a clear topography, with dominant negative correlations concentrated in the anterior dorsolateral striatum (aDLS) and dominant positive correlations in the posterior striatum",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "Sites with a positive dominant r have median AP -0.40 mm (41 sites) against 0.24 mm for negative sites (133), one-sided rank-sum p = 4.3e-07. The direction holds in all five mice and is significant in UG27, UG29 and UG31.",
     "deviations": [
      "The paper shows this on maps; the test of AP coordinates was chosen here."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "ed2-hotspot-per-mouse",
     "statement": "The hotspot is consistent across mice; in each mouse, sites inside the group hotspot have more negative dominant correlations than sites outside.",
     "source": "Extended Data Figure 2d; Results, Spatially organized coupling of transient DA-ACh dynamics during behavior, paragraph 1",
     "quote": "The presence and anatomical localization of this hotspot was remarkably consistent across mice (Extended Data Fig. 2).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "All five mice have lower median r inside the group hotspot; p < 0.05 in UG27, UG28, UG30 and UG31, and p = 0.084 (exact p = 0.087) in UG29, which has 2 sites inside.",
     "deviations": [
      "The per-mouse hotspots of Extended Data Figure 2b-c are not computed."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "ed1-controls-no-organization",
     "statement": "Mice expressing non-functional sensors (a mutant ACh sensor and tdTomato, n = 3) show only sparse, weak positive red-green correlations with no anatomical organization; no hotspot is found.",
     "source": "Extended Data Figure 1i-k; Results, Spatially organized coupling of transient DA-ACh dynamics during behavior, paragraph 1",
     "quote": "showing sparse weak positive correlations with no anatomical organisation. ... No significant hotspot was identified for these mice.",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "not_reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "In AD4 to AD6 (ACh4.3mut and tdTomato) the hotspot procedure finds a hotspot with mean smoothed r +0.33 (551 voxels, centroid AP 0.68, ML 1.29, DV -4.14 mm), so by the paper's test the claim of no organization fails. The procedure also finds a hotspot in 73 of 100 shuffles of the control values across sites, but none of those reaches 551 voxels (31 have a mean of at least +0.33), so by the paper's own test the control correlations are organized in space. The correlations are mostly positive (68% of 124 sites, median 0.20), as stated, but not weak; the median absolute dominant r is 0.238 against 0.240 in the sensor mice, and 65% of control sites exceed 0.2 against 54%. One control site has no defined correlation (a constant red trace) and is dropped; 19 site-sessions have constant red traces. Keeping the reward periods gives median r 0.22. Two deviations could add correlated red-green signal and so bear on this outcome; the controls use the high-pass series without artifact removal (none is in their files), and only their 15 reward sessions, where the paper used all ITI data.",
     "deviations": [
      "The 15 uncued reward sessions of the control mice are used, with 6 s after rewards and first licks left out; the paper used the ITI data of all their sessions, which the files' event tables would allow but which would take most of the time limit to stream.",
      "The control files hold no artifact-removed series; the high-pass series is used."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig2a-directionality",
     "statement": "For coincident transients of opposite polarity (ACh trough with DA peak, ACh peak with DA trough) the DA transient leads; for same-polarity pairs (peak with peak, trough with trough) the ACh transient leads.",
     "source": "Figure 2a; Results, Spatially organized coupling of transient DA-ACh dynamics during behavior, paragraph 2",
     "quote": "opposite signed pairs (cases 1 and 2) had timings indicating DA leading, while timings of same sign pairs (cases 3 and 4) indicated ACh leading (Fig. 2a).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "Median directionality indices are -0.27 (ACh trough with DA peak), -0.30 (ACh peak with DA trough), +0.33 (peak with peak) and +0.12 (trough with trough) over 174 sites, each with Wilcoxon p below 1e-24, and each mouse has the same four signs.",
     "deviations": [
      "Pairs are matched by the Methods' rule (each first transient takes its earliest partner, each partner used once), a simplification of the code's two-way matching.",
      "The code keeps only sites flagged as having signal in both channels (`sig`); these flags are not in the files and all 174 included sites are used, as the paper reports 174."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig2-magnitude-correlations",
     "statement": "Magnitudes of paired transients are negatively correlated at many sites for opposite-sign pairs (121 of 174 for ACh troughs with DA peaks, 81 of 174 for ACh peaks with DA troughs) and positively for same-sign pairs (97 of 174 peak-peak, 78 of 174 trough-trough), with no site significant in the other direction.",
     "source": "Results, Spatially organized coupling of transient DA-ACh dynamics during behavior, paragraph 3; Figure 2d",
     "quote": "For opposite sign event pairs, a majority of sites had significant negative magnitude correlations ... 121/174 sites significant for ACh troughs vs DA peaks; 81/174 sites for ACh peaks vs DA troughs ... for same sign pairs, a majority had positive correlations (97/174 significant for peak-peak correlations, 78/174 for trough-trough",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "partially_reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "Significant sites of the stated sign are 114, 69, 97 and 78 against 121, 81, 97 and 78; the second is 3.9 below its tolerance. The counts move by several sites with the random draws (the first run gave 115, 70, 90 and 76). No site is significant in the opposite direction, as stated. All sites have at least 100 pairs. The code's computation (Pearson r over all pairs, p < 0.05) finds 142 to 149 sites per type, so the paper's counts fit the Methods' bootstrap.",
     "deviations": [
      "1000 in place of 10000 repetitions of the 100-pair bootstrap."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig2d-magnitude-in-out",
     "statement": "Magnitude correlations of opposite-sign pairs, but not same-sign pairs, are stronger inside the aDLS hotspot than outside (rank-sum p < 0.05 for opposite-sign pairs, p > 0.05 for same-sign pairs).",
     "source": "Figure 2d; Results, Spatially organized coupling of transient DA-ACh dynamics during behavior, paragraph 3",
     "quote": "Magnitude correlations for opposite sign transient pairs, but not same sign pairs, were significantly stronger for sites inside the aDLS global correlation hotspot region (Fig. 1f) than outside",
     "statistics": [
      {
       "about": "opposite-sign pairs, inside against outside the hotspot",
       "test": "Wilcoxon rank-sum",
       "p": 0.05,
       "p_relation": "<"
      },
      {
       "about": "same-sign pairs, inside against outside the hotspot",
       "test": "Wilcoxon rank-sum",
       "p": 0.05,
       "p_relation": ">"
      }
     ],
     "kind": "finding",
     "inferred": "",
     "outcome": "partially_reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "ACh troughs with DA peaks are more negative inside (median -0.33 against -0.23, p = 0.0001); ACh peaks with DA troughs do not differ (-0.19 against -0.16, p = 0.11). Same-sign pairs do not differ (p = 0.11 and 0.22), as stated.",
     "deviations": [],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig2b-pairing-overlaps-hotspot",
     "statement": "The frequency of paired transients is highest in aDLS, in a region overlapping the correlation hotspot, for all four pair types (volume overlap test, p < 0.01 for each).",
     "source": "Figure 2b-c; Results, Spatially organized coupling of transient DA-ACh dynamics during behavior, paragraph 2",
     "quote": "The frequencies of all paired events were highest in the aDLS within a region significantly overlapping with the hotspot identified by correlation analysis",
     "statistics": [
      {
       "test": "volume overlap test",
       "p": 0.01,
       "p_relation": "<"
      }
     ],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "The occurrence hotspot of each pair type overlaps the correlation hotspot by 780 to 826 voxels against null medians of 88 to 111 (overlap p 0.005, 0.015, 0.005 and below 0.005; the paper 0.005, 0.005, 0.002 and below 0.001). Occurrence is higher inside than outside for each type (p < 0.001). With the chance occurrence subtracted (DA and ACh transients from different sessions of the same site, as in the Methods), every site keeps an excess for every type, the excess is higher inside (p < 0.001), and its hotspots overlap the correlation hotspot with p 0.030, 0.005, 0.015 and below 0.005.",
     "deviations": [
      "200 random volumes in the overlap test, so p cannot fall below 0.005.",
      "The chance occurrence is the mean of 10 shuffles of sessions, in place of 10000."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "ed3e-transient-rates-in-hotspot",
     "statement": "DA and ACh transient rates are higher at sites inside the global correlation hotspot than outside.",
     "source": "Extended Data Figure 3d-e",
     "quote": "Sites inside the hotspot showed significantly higher transient rates",
     "statistics": [
      {
       "test": "Wilcoxon rank-sum",
       "p": 0.05,
       "p_relation": "<"
      }
     ],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "All four transient types are more frequent inside the hotspot (for example DA peaks 0.226 against 0.176 per second of usable data), p of 0.0013 or less. Rates are now per second of data that is not NaN and not within 6 s after a reward; the first run divided by the whole session and gave the same result.",
     "deviations": [],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig3b-reward-correlation",
     "statement": "After unpredicted rewards, trial-by-trial DA peak magnitude correlates with ACh; dominant correlations are negative at 114 of 174 sites with DA leading (mean latency 177.9 +/- 96.7 ms) and positive at 58 of 174 with ACh leading (-41.2 +/- 218.9 ms).",
     "source": "Results, Conserved spatial structure of distinct anti-correlated dynamics to unpredicted rewards and Pavlovian cues following conditioning and extinction, paragraph 1; Figure 3b",
     "quote": "sites with dominant negative correlations (114/174 sites) indicated the DA peak leading ACh at a latency aligning with the ACh dip (177.9 + 96.7ms mean latency), and sites with dominant positive correlations indicated the DA peak lagging ACh (58/174 sites)",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "partially_reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "122 sites are dominant negative (mean latency 172.1 ms, 120 with DA leading) and 52 dominant positive (mean latency -59.8 ms, 37 with ACh leading). Three values are within tolerance; the positive count is 0.2 below its range (52.2 to 63.8). Negative sites outnumber positive ones in every mouse.",
     "deviations": [
      "The latency tolerances are the standard errors implied by the paper's spreads.",
      "One UG27 reward session (240313r) is skipped for lacking the artifact-removed green series."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig3d-reward-in-out",
     "statement": "Reward-evoked DA-ACh correlations are significantly more negative inside the global hotspot than outside.",
     "source": "Figure 3d; Results, Conserved spatial structure of distinct anti-correlated dynamics to unpredicted rewards and Pavlovian cues following conditioning and extinction, paragraph 1",
     "quote": "Correlations were significantly more negative inside than outside the hotspot (Fig. 3d).",
     "statistics": [
      {
       "test": "Wilcoxon rank-sum",
       "p": 0.05,
       "p_relation": "<"
      }
     ],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "Median r -0.52 inside against -0.31 outside, p = 7.3e-09; same sign in all five mice, significant in three.",
     "deviations": [],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig3e-reward-hotspot-overlap",
     "statement": "Reward anti-correlations form an aDLS hotspot that significantly overlaps the hotspot of spontaneous correlations (volume overlap test, p = 0.008).",
     "source": "Figure 3c, 3e; Results, Conserved spatial structure of distinct anti-correlated dynamics to unpredicted rewards and Pavlovian cues following conditioning and extinction, paragraph 1",
     "quote": "The anatomical organization of these correlations indicated a 'hotspot' of dominant negative correlations in the aDLS which significantly overlapped with the hotspot identified from spontaneous events (Fig. 3c,e).",
     "statistics": [
      {
       "test": "volume overlap test",
       "p": 0.008
      }
     ],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "The reward hotspot has mean -0.53 (1022 voxels, centroid AP 0.54, ML 2.07, DV -3.26 mm) and shares 721 voxels with the spontaneous hotspot against a null median of 116, overlap p = 0.01 (the paper 0.008).",
     "deviations": [
      "200 random volumes in the overlap test."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig4-cue-sequence",
     "statement": "After Pavlovian conditioning, the conditioned light cue evokes an ACh peak, then a DA peak, then an ACh dip at a majority of sites (49 of 75).",
     "source": "Results, Conserved spatial structure of distinct anti-correlated dynamics to unpredicted rewards and Pavlovian cues following conditioning and extinction, paragraph 2; Figure 4a",
     "quote": "conditioned cues evoked a temporal sequence of ACh peak\u2192 DA peak\u2192 ACh dip at a majority of sites (49/75 sites).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "not_reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "The sequence is found at 58 of 75 sites under both readings (code percentiles 55 and 52), above the tolerance of 44.1 to 53.9. The qualitative statement holds; a majority of sites shows the sequence overall and in each mouse (12 of 17, 19 of 24, 27 of 34). The criterion counted only a shortfall as partial; a count above the range is not reproduced by its wording.",
     "deviations": [
      "The paper does not define how the sequence is identified; significant ACh peak, DA peak and ACh dip with ordered extremes in the mean trace was chosen here.",
      "Learning sessions are identified from the events (every light cue rewarded); the last 3 per mouse are used, following the code's intent (see findings on the code).",
      "Tone trials cannot be analysed (no tone events); the paper's figure uses light trials."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig4c-cue-correlation",
     "statement": "For conditioned light cues, ACh correlates with the trial-by-trial DA peak negatively at lagging times at 39 of 75 sites (ACh lagging, 202.28 +/- 128.4 ms) and positively at leading times at 26 of 75 (ACh leading, -79.1 +/- 224.0 ms).",
     "source": "Results, Conserved spatial structure of distinct anti-correlated dynamics to unpredicted rewards and Pavlovian cues following conditioning and extinction, paragraph 2; Figure 4c",
     "quote": "Correlations of the ACh \u0394F/F with the peak DA \u0394F/F were negative at lagging timepoints (39/75 sites; ACh lagging DA; 202.28 + 128.4ms) and positive at leads (26/75 sites; ACh leading DA; - 79.1 + 224.0ms)",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "partially_reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "All 75 sites have a significant cue DA peak, so the two readings agree. 41 sites are negative (mean latency 206.0 ms, 37 with ACh lagging) and 34 positive (-83.3 ms, 22 with ACh leading). The positive count is above its range (23.4 to 28.6); the paper's 39 plus 26 make 65, which fits neither reading here. The code's percentiles also give 75 sites, and the null draws from 68868 to 74196 candidate times per mouse, so neither explains the difference. The difference remains unexplained.",
     "deviations": [
      "The denominator is ambiguous (75 sites or those with a significant DA peak); both readings are shown and agree here."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig4h-extinction-correlation",
     "statement": "After cue-specific extinction (light cue rewarded on 20% of trials), cue-evoked DA dips correlate negatively with lagging ACh (DA dip then ACh late peak) at 41 of 74 sites, mean lag 117.9 +/- 176.2 ms.",
     "source": "Results, Conserved spatial structure of distinct anti-correlated dynamics to unpredicted rewards and Pavlovian cues following conditioning and extinction, paragraph 3; Figure 4g-h",
     "quote": "the timing relationships for correlations with DA dip magnitudes closely resembled those for the DA peak sites with dominant negative correlations indicating ACh lagging (DA dip\u2192Ach late peak) (Fig. 4h; 41/74 sites; mean lag 117.9 + 176.2ms).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "Under reading (a), 38 of 61 sites with a significant DA dip are negative with mean lag 118.4 ms (the paper 41 of 74, 117.9 ms). Over all 75 sites, 47 negative with mean lag 114.7 ms. Only 28 of the 38 have ACh lagging, and AD2's negative sites have a mean lag of -79.4 ms. The code's polarity (DA peak) gives 44 sites with mean lag 137.6 ms.",
     "deviations": [
      "The DA trough is used, as the legend and Methods state; the code correlates the DA peak in this step, shown as a variant.",
      "The paper's denominator (74) and the legend's 48 sites in 5 mice do not match each other or the 75 sites in 3 mice of the files; two readings are shown."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig4d-cue-in-out",
     "statement": "Conditioned-cue DA-ACh correlations are more negative at sites inside the global aDLS hotspot than outside.",
     "source": "Figure 4d-e; Results, Conserved spatial structure of distinct anti-correlated dynamics to unpredicted rewards and Pavlovian cues following conditioning and extinction, paragraph 2",
     "quote": "The spatial structure of DA-ACh correlations was also similar, with a significant aDLS hotspot of dominant anti-correlations overlapping with the hotspot for paired spontaneous events (Fig. 4b-e).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "The 9 AD1 to AD3 sites inside the global hotspot have median r -0.31 against 0.02 for the 66 outside, p = 0.024 (exact p = 0.022). The cue hotspot's own overlap with the global hotspot is not significant (p = 0.175), unlike the paper's Figure 4e.",
     "deviations": [
      "Sites of the Pavlovian cohort are placed in the global hotspot by their coordinates on the shared grid."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig4j-extinction-in-out",
     "statement": "After extinction, DA dip to ACh correlations are more negative at sites inside the global aDLS hotspot than outside.",
     "source": "Figure 4i-j; Results, Conserved spatial structure of distinct anti-correlated dynamics to unpredicted rewards and Pavlovian cues following conditioning and extinction, paragraph 3",
     "quote": "The spatial organization of dominant negative correlations was also aligned, significantly overlapping with the aDLS correlation hotspot (Fig. 4g,i,j; Extended Data Fig. 5e,f).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "Among dip sites, median r -0.35 inside (8 sites) against -0.21 outside (53), p = 0.035 (exact p = 0.034); all sites p = 0.018. The extinction hotspot's overlap with the global one is not significant (p = 0.17).",
     "deviations": [],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig5d-initiation-responses",
     "statement": "At locomotion initiations, ACh rises and then falls at 131 of 174 sites; DA dips at 119 of 174 sites, rises at 73, and both rises and dips at 32.",
     "source": "Results, Spatially organized anti-correlations corresponding to phase specific locomotion signals, paragraph 1; Figure 5d",
     "quote": "initiations from rest were associated with rapid ACh increases at many sites, followed by decreases (131/174 sites, Fig. 5a-c). In contrast, DA release at some sites (119/174 sites) rapidly decreased, and at others in the same mouse, rapidly increased (73/174 sites) ... only a minority of sites had both DA increases and decreases (32/174 sites).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "partially_reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "DA dips at 129 sites and rises at 77 are within 10% (the paper 119 and 73); both at 46 is not (32). ACh rise then fall at 107 sites under reading (a) and a rise at 147 under reading (b), neither within 117.9 to 144.1. Code percentiles give 114 dips, 61 rises and 26 both. The code's bout of 2.2 s gives 132, 74 and 44. UG27 contributes 501 initiations against 92 to 149 for the others; its wheel velocity exceeds the threshold 49% of the time.",
     "deviations": [
      "A zero-phase Butterworth filter replaces MATLAB's lowpass.",
      "1000 bootstrap samples in place of 10000.",
      "The Methods' bout of 2.27 s is used; the code's 2.2 s is a variant."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig5-latency-order",
     "statement": "Initiation-related signals occur in sequence; DA dips first (median latency 111.11 ms), then ACh peaks (222.22 ms) and DA peaks (227.78 ms), then ACh dips (611.11 ms).",
     "source": "Results, Spatially organized anti-correlations corresponding to phase specific locomotion signals, paragraph 1",
     "quote": "DA dips were earliest (median latency 111.11, IQR 55.56-708.33 ms), followed by ACh peaks and DA peaks (median latencies 222.22, IQR 222.22-277.78 ms; and 227.78, IQR -388.89-444.44 ms, respectively), then ACh dips (median latency 611.11, IQR 500.00-722.22 ms).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "Medians are 111.1 ms (DA dips), 277.8 ms (ACh peaks), 277.8 ms (DA peaks) and 666.7 ms (ACh dips), in the stated order; the ACh peak and dip medians are exactly one sample later than the paper's, at the edge of the tolerance. The per-mouse order is irregular (DA peak medians from -388.9 to 527.8 ms).",
     "deviations": [
      "The Methods do not say how latency is measured; the time of the extreme of the triggered average in -0.5 to 1 s was chosen here."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig5e-initiation-topography",
     "statement": "Initiation-related DA dips are concentrated in the aDLS and DA peaks in dorsomedial and posterior striatum; ACh dips are localized to the aDLS.",
     "source": "Results, Spatially organized anti-correlations corresponding to phase specific locomotion signals, paragraph 1; Figure 5e-g",
     "quote": "early DA dips were concentrated in the aDLS, followed by DA peaks in the dorsomedial and posterior striatum (Fig. 5e,f; Extended Data Fig. 6a). ... later ACh dips were localized to the aDLS (Fig. 5e,g).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "not_reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "Dominant DA dips are more anterior than DA rises (median AP 0.36 against -0.48 mm, p < 0.0001) but more medial (ML 1.88 against 2.30 mm) and ventral; ACh-dip sites (148 of 174) are not more anterior or lateral than the 26 others (p = 0.33 and 0.25). All 28 hotspot sites have a dominant DA dip and an ACh dip. The paper describes maps; with ACh dips at 85% of sites a contrast between sites with and without dips may not capture where dips are strongest.",
     "deviations": [
      "The paper shows the topography as smoothed maps; site-level coordinate tests were chosen here."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig6a-initiation-correlation",
     "statement": "At sites with initiation-related DA dips (119 sites), trial-by-trial DA dip magnitudes correlate negatively with ACh with DA leading (mean latency 193.80 +/- 177.68 ms).",
     "source": "Results, Spatially organized anti-correlations corresponding to phase specific locomotion signals, paragraph 3; Figure 6a",
     "quote": "dominant negative correlations were present between DA and ACh with timings indicating DA leading (Fig. 6a; mean latency 193.80 +/- 177.68ms).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "Of 129 sites with a DA dip (the paper 119), 98 are dominant negative (88 with ACh lagging), mean latency 188.8 ms under reading (a); 98.2 ms under reading (b). Most dip sites are negative in every mouse. The trough search now ignores NaNs; before the fix a trough within 2 samples of a NaN was missed. 2197 of 30341 site-trials have a NaN in the range searched, so at most that many were affected. Before the fix the result was 99 sites and 182.4 ms.",
     "deviations": [
      "The latency tolerance is the standard error implied by the paper's spread."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig6c-initiation-in-out",
     "statement": "Initiation DA dip to ACh correlations are concentrated in the aDLS and more negative inside the global hotspot than outside.",
     "source": "Figure 6b-d; Results, Spatially organized anti-correlations corresponding to phase specific locomotion signals, paragraph 3",
     "quote": "these correlations were concentrated in the aDLS, the area expressing the largest initiation evoked DA dips, and this region significantly overlapped with the previously identified hotspot (Fig. 6b-d).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "Median r -0.35 inside (28 sites) against -0.21 outside (101), p = 1.7e-05 (6.1e-08 before the fix of the trough search).",
     "deviations": [],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig6g-invigoration-in-out",
     "statement": "At locomotion invigorations, GLM-estimated DA peak magnitudes correlate negatively with ACh at lagging times, and these correlations are concentrated in the aDLS, more negative inside the global hotspot (131 sites with invigoration DA peaks).",
     "source": "Results, Spatially organized anti-correlations corresponding to phase specific locomotion signals, paragraph 3; Figure 6e-h; Methods, generalized linear model",
     "quote": "Again, dominant negative correlations were observed at lags relative to the DA peak and were selectively concentrated in the aDLS, significantly overlapping with the previously identified anatomical hotspot (Fig. 6e-h;).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "The GLM finds invigoration DA peaks at 160 sites (the paper 131); 105 are dominant negative (97 with ACh lagging, mean lag 179.4 ms); median r -0.34 inside (28) against -0.21 outside (132), p = 4.5e-05.",
     "deviations": [
      "Ordinary least squares in numpy in place of MATLAB fitglm (normal family, identical model)."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig7-da-release",
     "statement": "Optogenetic stimulation of midbrain DA neurons evokes significant DA release at 179 of 180 sites.",
     "source": "Results, Optogenetic dopamine stimulation drives spatially selective ACh suppression in the hotspot, paragraph 1; Extended Data Figure 10b",
     "quote": "Significant stimulation evoked release was measured at all sites but one (99%, 179 of 180 fibers; Extended Data Fig. 10b).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "179 of 180 sites show a significant DA peak after ITI stimulation (178 with the code's percentiles); 115 ITI trains in total (20 per mouse, 15 in ADS6).",
     "deviations": [],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig7e-ach-suppression",
     "statement": "DA stimulation evokes significant ACh suppression at a subset of sites (116 of 180, 64%).",
     "source": "Results, Optogenetic dopamine stimulation drives spatially selective ACh suppression in the hotspot, paragraph 1; Figure 7e",
     "quote": "significant ACh suppression following DA stimulation was observed only in a subset of fibers (64%, 116 of 180 sites), and these responses were spatially concentrated in the aDLS (Fig. 7b-e).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "partially_reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "147 of 180 sites (82%) show a significant ACh dip against 116 (64%); 128 (71%) with the code's percentiles. The 33 sites without a dip are mostly in ADS19 (17 of 27) and ADS16 (9 of 37).",
     "deviations": [
      "1000 bootstrap samples in place of 10000."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "ed9-cohort-hotspot",
     "statement": "The stimulation cohort replicates the aDLS anti-correlation hotspot in spontaneous release outside stimulation periods; its dominant correlations are more negative inside the original hotspot, and its own hotspot overlaps the original.",
     "source": "Extended Data Figure 9a-d; Results, Optogenetic dopamine stimulation drives spatially selective ACh suppression in the hotspot, paragraph 1",
     "quote": "We confirmed that an aDLS DA-ACh anti-correlation hotspot was present in this separate cohort of mice which significantly overlapped with that reported in previous cohorts (Extended Data Fig. 9a-d).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "partially_reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "The 9 cohort sites inside the original hotspot are more negative (median -0.42 against -0.30, p = 0.049, exact p = 0.048; OR exclusion p = 0.041). The cohort's own hotspot is negative (mean -0.44, 1536 voxels, centroid AP 0.71, ML 2.00, DV -3.45 mm) and shares 564 voxels with the original against a null median of 65, but overlap p = 0.085 with 200 random volumes. 167 of 180 cohort sites are dominant negative. The AND exclusion leaves out 6.5% of samples and the OR exclusion 32.7%. Of the samples within 3 s after a stimulation, AND keeps 42% and OR none. With OR the cohort hotspot is similar (mean -0.41, centroid AP 0.72, ML 2.08, DV -3.43 mm), with overlap p = 0.085.",
     "deviations": [
      "The code's AND of the three exclusions is followed; OR is a variant."
     ],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig7f-ach-dip-in-out",
     "statement": "Stimulation-evoked ACh dips are larger inside the cohort's correlation hotspot than outside, and their own hotspot overlaps it.",
     "source": "Figure 7d, 7f-g; Results, Optogenetic dopamine stimulation drives spatially selective ACh suppression in the hotspot, paragraph 1",
     "quote": "The spatial concentration of ACh stimulation-evoked dips significantly overlapped with the anti-correlation hotspot for spontaneous release events ... ACh \u0394F/F dip magnitudes were significantly higher, on average, inside than out (Fig. 7f,g).",
     "statistics": [],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "ACh minima are deeper inside the cohort hotspot (median -0.0179 against -0.0149 dF/F, 17 against 163 sites, p = 0.0005). The AND exclusion keeps 42% of the samples within 3 s after a stimulation, so its cohort hotspot is built partly from stimulation responses; with the hotspot of the OR exclusion, which leaves them out, the result holds (19 against 161 sites, median -0.0188 against -0.0149, p = 0.0003). Against the original hotspot p = 0.094; the ACh-minimum hotspot's overlap with the cohort hotspot has p = 0.065, so the claim's second part is not supported at p < 0.05; the criterion treats it as a variant. Inside sites are lower in all six mice, significantly in ADS19.",
     "deviations": [],
     "disagreed": false,
     "data": "dataset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig8b-quinpirole-region",
     "statement": "The D2 agonist quinpirole suppresses evoked ACh release in DLS and NAcC, more strongly in DLS; two-way ANOVA region effect F(1,40) = 105.9, region by stimulus interaction F(4,40) = 3.014, stimulus effect F(4,40) = 17.80.",
     "source": "Figure 8b; Results, Differential D2R-mediated inhibition across striatal regions inside and outside the aDLS hotspot, paragraph 1",
     "quote": "Quinpirole suppressed evoked ACh release in both DLS and NAcC (Fig. 8a,b), but the magnitude of suppression was significantly greater in DLS across all stimulation pulses (Fig. 8b).",
     "statistics": [
      {
       "about": "region by stimulus number interaction",
       "test": "two-way ANOVA",
       "statistic": "F(4,40) = 3.014",
       "p": 0.029,
       "n": "5 slices per region from 5 mice"
      },
      {
       "about": "region effect",
       "test": "two-way ANOVA",
       "statistic": "F(1,40) = 105.9",
       "p": 0.0001,
       "p_relation": "<",
       "n": "5 slices per region from 5 mice"
      },
      {
       "about": "stimulus number effect",
       "test": "two-way ANOVA",
       "statistic": "F(4,40) = 17.80",
       "p": 0.0001,
       "p_relation": "<",
       "n": "5 slices per region from 5 mice"
      }
     ],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "Region F 105.87, interaction F 3.014 and stimulus F 17.80 match the paper; DLS mean 25.1 against NAcC 68.1. Type III sums of squares are used, as Prism computes them; Type II gives the same values here.",
     "deviations": [],
     "disagreed": false,
     "data": "elsewhere",
     "where": "Zenodo doi:10.5281/zenodo.18202873, file Bouabid_Vu (2026)_Source_Data_Fig8.xlsx",
     "proposed": null
    },
    {
     "claim": "fig8e-cocaine-region",
     "statement": "Cocaine reduces evoked ACh release more in DLS than in NAcC; two-way ANOVA stimulus by region interaction F(4,55) = 10.75, stimulus effect F(4,55) = 36.37, region effect F(1,55) = 199.9.",
     "source": "Figure 8e; Results, Differential D2R-mediated inhibition across striatal regions inside and outside the aDLS hotspot, paragraph 2",
     "quote": "cocaine robustly reduced ACh release at the second through fourth pulses in DLS, while producing only modest suppression in NAcC (Fig. 8d,e).",
     "statistics": [
      {
       "about": "stimulus number by region interaction",
       "test": "two-way ANOVA",
       "statistic": "F(4,55) = 10.75",
       "p": 0.0001,
       "p_relation": "<"
      },
      {
       "about": "stimulus number effect",
       "test": "two-way ANOVA",
       "statistic": "F(4,55) = 36.37",
       "p": 0.0001,
       "p_relation": "<"
      },
      {
       "about": "region effect",
       "test": "two-way ANOVA",
       "statistic": "F(1,55) = 199.9",
       "p": 0.0001,
       "p_relation": "<"
      }
     ],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "With Type III sums of squares, interaction F 10.75, stimulus F 36.37 and region F 199.86 match the paper (10.75, 36.37 and 199.9); DLS 32.8 against NAcC 74.1. The design is unbalanced (8 DLS and 5 NAcC slices). The first run used Type II sums of squares, which give stimulus F 47.68, and scored the claim partially reproduced.",
     "deviations": [],
     "disagreed": true,
     "data": "elsewhere",
     "where": "Zenodo doi:10.5281/zenodo.18202873, file Bouabid_Vu (2026)_Source_Data_Fig8.xlsx",
     "proposed": null
    },
    {
     "claim": "fig8fg-d2-block",
     "statement": "The D2 antagonist L-741,626 attenuates the cocaine-induced suppression in both regions; D2R block effect F(1,55) = 231.9 in DLS and F(1,40) = 32.57 in NAcC.",
     "source": "Figure 8f-g; Results, Differential D2R-mediated inhibition across striatal regions inside and outside the aDLS hotspot, paragraph 2",
     "quote": "This effect was significantly attenuated by the D2R antagonist L-741,626 in both regions (Fig. 8f,g)",
     "statistics": [
      {
       "about": "D2R block effect in DLS",
       "test": "two-way ANOVA",
       "statistic": "F(1,55) = 231.9",
       "p": 0.0001,
       "p_relation": "<"
      },
      {
       "about": "D2R block effect in NAcC",
       "test": "two-way ANOVA",
       "statistic": "F(1,40) = 32.57",
       "p": 0.0001,
       "p_relation": "<"
      }
     ],
     "kind": "finding",
     "inferred": "",
     "outcome": "reproduced",
     "blocked_by": null,
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "Block F 231.87 (DLS) and 32.57 (NAcC) match; the antagonist raises ACh in both (32.8 to 78.6, 74.1 to 107.9). With Type III sums of squares the DLS stimulus F is 34.46, as in the legend (Type II 45.04).",
     "deviations": [],
     "disagreed": false,
     "data": "elsewhere",
     "where": "Zenodo doi:10.5281/zenodo.18202873, file Bouabid_Vu (2026)_Source_Data_Fig8.xlsx",
     "proposed": null
    },
    {
     "claim": "fig8c-cocaine-da-similar",
     "statement": "Cocaine elevates evoked DA release to a similar degree in DLS and NAcC (two-way ANOVA, cocaine by region interaction F(8,36) = 0.20, p = 0.98).",
     "source": "Figure 8c; Results, Differential D2R-mediated inhibition across striatal regions inside and outside the aDLS hotspot, paragraph 2",
     "quote": "we enhanced DA availability with the reuptake inhibitor cocaine, which elevates DA release by similar degree in both regions (Fig. 8c).",
     "statistics": [
      {
       "about": "cocaine by region interaction",
       "test": "two-way ANOVA",
       "statistic": "F(8,36) = 0.20",
       "p": 0.98,
       "n": "3 slices from 3 mice"
      },
      {
       "about": "cocaine effect",
       "test": "two-way ANOVA",
       "statistic": "F(8,36) = 5.83",
       "p": 0.0001,
       "p_relation": "<",
       "n": "3 slices from 3 mice"
      },
      {
       "about": "region effect",
       "test": "two-way ANOVA",
       "statistic": "F(1,36) = 0.17",
       "p": 0.67,
       "n": "3 slices from 3 mice"
      }
     ],
     "kind": "finding",
     "inferred": "",
     "outcome": "blocked",
     "blocked_by": "methods_underspecified",
     "notebook": "reproduce_bouabid_2026.ipynb",
     "notes": "The sheet Figure_8c_TP holds a voltammogram (voltage and current). Figure_8c_LP and _RP hold mean DA time courses (24 time points from 0.625 to 3.5 s) for three mice per region and two conditions. The legend's ANOVA has 8 degrees of freedom for cocaine and 36 residual, which fits neither a two-level cocaine factor nor the 24 time points, so the test cannot be repeated. Descriptively, the cocaine peak is 1.48 times the control peak on average in both regions (per mouse 1.31 to 1.77 in DLS, 1.10 to 2.09 in NAcC).",
     "deviations": [],
     "disagreed": false,
     "data": "elsewhere",
     "where": "Zenodo doi:10.5281/zenodo.18202873, file Bouabid_Vu (2026)_Source_Data_Fig8.xlsx",
     "proposed": null
    }
   ],
   "findings": [
    {
     "kind": "metadata_error",
     "statement": "Every photometry series states a sampling rate of 9 Hz (starting_time rate 8.999 to 9.0). Each series has as many samples as the behavior series, which steps at 55.6 ms (18 Hz), and the Methods and code use 18 Hz; at the stated rate the recordings would last twice as long as their timestamps."
    },
    {
     "kind": "metadata_insufficient",
     "statement": "No field in the files states which experiment (spontaneous, unpredicted reward, Pavlovian learning or extinction, control, stimulation) a session belongs to, or which Pavlovian phase. Cohorts can be told only from subject names and indicators, and phases only by counting rewarded light cues (1.0, 0.8 or 0.2). The paper's Zenodo record has a table of this mapping outside the dandiset."
    },
    {
     "kind": "data_missing",
     "statement": "The Pavlovian files have light cue events but no tone cue events, although each session has 30 tone cues (the 68 rewards of a learning session are 30 light, 30 tone and 8 uncued). Tone trials cannot be analysed."
    },
    {
     "kind": "data_missing",
     "statement": "Session sub-UG27_ses-240313r lacks the high-pass, artifact-removed green series that all other UG sessions hold, so it cannot enter the analyses that use it."
    },
    {
     "kind": "metadata_error",
     "statement": "The device description of ExcitationSource470 describes 470 nm and 405 nm (isosbestic) light alternated at 100 Hz, and that of ExcitationSource570 reads \"405 nm\"; the Methods alternate 470 and 570 nm at 36 Hz. The mutant sensor is named ACh4.3mut in the files where the Methods and Key Resources Table name ACh3.0-mut. The session description gives 9 uL rewards at 5 to 30 s, where the Methods give 5 uL at 10 to 30 s, and stores the micro sign mis-encoded."
    },
    {
     "kind": "metadata_insufficient",
     "statement": "The code keeps only sites with signal in both channels (a `sig` flag per site) for the transient analyses; the flags are not in the files. All 174 included sites gave the paper's count, so this did not change the result here."
    },
    {
     "kind": "code_mismatch",
     "statement": "In fig4_pav.m (commit 833fa3d) the loop over a phase's last three sessions loads exp_dirs{d} in place of task_dirs{d}, so it reads the mouse's first three recordings for every phase, and the extinction correlation is called with polarity 1 (the DA peak) although the legend and Methods describe the DA trough. The notebook follows the evident intent and the Methods."
    },
    {
     "kind": "code_mismatch",
     "statement": "The code departs from the Methods in several parameters; triggered-average significance uses the 0.5th and 99.5th percentiles (Methods 2.5th and 97.5th), the minimum locomotion bout is 2.2 s (Methods 2.27 s), paired-transient magnitude correlations are Pearson r over all pairs (Methods a 100-pair bootstrap), and the stimulation cohort's exclusion of reward and stimulation periods combines the three conditions with AND, so a sample is dropped only when all three hold. hotspot_comparison.m places the second hotspot's random volumes with the first grid's dimensions."
    },
    {
     "kind": "paper_inconsistency",
     "statement": "For extinction, the Results give 41 of 74 sites and the Figure 4g legend 48 sites in 5 mice, while the Methods give 3 Pavlovian mice and the files hold 75 sites in 3 mice. For learned cues the Results give 39 of 75 negative and 26 of 75 positive sites, which sum to 65."
    },
    {
     "kind": "paper_inconsistency",
     "statement": "In the Figure 8 source data, the Figure 8b DLS group holds 6 slice labels from 4 mouse labels (slice_B_4 holds only stimulus 5) where the legend states 5 slices from 5 mice."
    },
    {
     "kind": "other",
     "statement": "The ex vivo slice data of Figure 8 are not in the dandiset; their source data are in the paper's Zenodo record (doi:10.5281/zenodo.18202873)."
    }
   ],
   "verifications": [
    {
     "date": "2026-10-19",
     "kind": "agent",
     "by": "Independent agent run by CatalystNeuro (same harness, separate context, did not write the notebook)",
     "reviewed": "The executed notebook of the first run with the conclusion cells and the summary removed, the helper module and the claims file. It also read the paper and the Zenodo Figure 8 source data, and queried the DANDI asset list.",
     "agreed": 32,
     "of": 33,
     "disagreements": [
      {
       "claim": "fig8e-cocaine-region",
       "outcome": "reproduced",
       "reason": "Cell 81 used Type II sums of squares, giving stimulus F = 47.68. Rerun with Type III (GraphPad Prism, which the Methods name, uses Type III; the design is unbalanced at 8 against 5 slices) gives 36.374, 10.754 and 199.857, all within tolerance (B2)."
      }
     ]
    },
    {
     "date": "2026-10-19",
     "kind": "agent",
     "by": "Independent agent run by CatalystNeuro (same harness, separate context, did not write the notebook)",
     "reviewed": "The full executed notebook of the revised run (121 cells), the helper module, the claims file and the attempt record. It also read the paper, checking the Figure 8 legend statistics against it.",
     "agreed": 33,
     "of": 33,
     "disagreements": []
    }
   ]
  }
 ],
 "entries": [
  {
   "name": "001692-bouabid-2026",
   "type": "reproduction",
   "title": "Striatal dopamine-acetylcholine hotspot (Bouabid, Vu et al. 2026)",
   "description": "Tests the findings of Bouabid, Vu et al. (2026) on dandiset 001692, which holds multi-fiber photometry of dopamine and acetylcholine across the striatum of 17 mice. It reproduces the spontaneous cross-correlations and their anterior dorsolateral hotspot, paired transients, reward, Pavlovian cue, locomotion and optogenetic stimulation analyses, and the control mice, and tests the ex vivo slice statistics of Figure 8 on the paper's source data on Zenodo.",
   "datasets": [
    "dandi:001692"
   ],
   "papers": [
    "10.64898/2026.01.20.700614"
   ],
   "derived_from": [],
   "notebooks": [
    {
     "path": "reproduce_bouabid_2026.ipynb",
     "title": "Reproducing Bouabid, Vu et al. (2026) from dandiset 001692",
     "executed": "notebooks/001692-bouabid-2026/reproduce_bouabid_2026.html"
    }
   ]
  }
 ]
}
