{
 "dandiset": "001692",
 "page": "dandisets/001692/",
 "generated": "2026-10-19T12:04:40+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
 },
 "attempts": [
  {
   "id": "001692-bouabid-2026-01",
   "dandiset": "001692",
   "dandiset_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"
   },
   "tier": "agent_verified",
   "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, first section; Figure 1e",
     "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": "dandiset",
     "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, first section; Figure 1f",
     "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": "dandiset",
     "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",
     "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": "dandiset",
     "where": null,
     "proposed": null
    },
    {
     "claim": "fig1f-positive-posterior",
     "statement": "Dominant positive correlations lie in the posterior striatum.",
     "source": "Results, first section; Figure 1f",
     "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": "dandiset",
     "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, first section",
     "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": "dandiset",
     "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, first section",
     "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": "dandiset",
     "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, second section",
     "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": "dandiset",
     "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, second section; Figure 2d",
     "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": "dandiset",
     "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, second section",
     "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": "dandiset",
     "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, second section",
     "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": "dandiset",
     "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",
     "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": "dandiset",
     "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, third section; Figure 3b",
     "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": "dandiset",
     "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",
     "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": "dandiset",
     "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",
     "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": "dandiset",
     "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, fourth section; Figure 4a",
     "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": "dandiset",
     "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, fourth section; Figure 4c",
     "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": "dandiset",
     "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, fourth section; Figure 4g-h",
     "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": "dandiset",
     "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",
     "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": "dandiset",
     "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",
     "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": "dandiset",
     "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, fifth section; Figure 5d",
     "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": "dandiset",
     "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, fifth section",
     "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": "dandiset",
     "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, fifth section; Figure 5e-g",
     "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": "dandiset",
     "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, sixth section; Figure 6a",
     "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": "dandiset",
     "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",
     "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": "dandiset",
     "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, sixth section; Figure 6e-h; Methods, generalized linear model",
     "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": "dandiset",
     "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, seventh section; Extended Data Figure 10b",
     "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": "dandiset",
     "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, seventh section; Figure 7e",
     "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": "dandiset",
     "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, seventh section",
     "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": "dandiset",
     "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",
     "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": "dandiset",
     "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",
     "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",
     "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",
     "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",
     "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.",
   "dandisets": [
    "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"
    }
   ]
  }
 ]
}