CSD drives wave-like calcium activity in meningeal afferents

(A) Mice received a trigeminal ganglion injection of an AAV to express GCaMP6s in trigeminal meningeal afferents. After 8-10 weeks, following the implantation of a headpost and a cranial window, mice were habituated to head restraint and subjected to two-photon calcium imaging while head-fixed on a running wheel to study the effect of pinprick-triggered CSD on the activity of meningeal afferents. (B) Example of a CSD-associated meningeal calcium wave that spreads across the field of view, with local segments of long afferent fibers becoming sequentially activated as the wave progresses (arrowheads). M: medial, L: lateral, A: anterior, P: posterior. (C) Summary of speed and direction of CSD-associated meningeal calcium waves, typically from anterior (‘Ant.’) (closer to where CSD was triggered anterior to the cranial window) to more posterior locations (‘Post.’). Speed estimates were obtained using the analysis method described in Figure S1A. On average, the wave progressed at 3.8±0.2 mm/min. (D) Map of 18 ROIs belonging to a single meningeal afferent fiber visible in B. (E) Activity heatmap of the afferent ROIs indicated in D illustrates progressive activation in response to CSD. (F) In contrast, the same afferent ROIs become activated simultaneously during a locomotion bout. (G) The pace of the CSD-associated afferent calcium wave was much slower than the spread of activity along the same afferent fibers during locomotion-evoked activity pre-CSD (**** P < 0.0001, paired, two-tailed t-test). (H) Example heatmaps of afferent activity observed during CSD showing different timecourse and magnitudes when compared to the activity observed during a locomotion bout. (I) Comparisons across all FOVs indicate a higher proportion of afferents exhibiting acute activation during the CSD vs. during locomotion (**** P < 0.0001, iterated bootstrap). (J) A higher proportion of afferents also displayed increased magnitudes of activation (* P < 0.05, paired, two-tailed t-test). See also Movie S1 and supplementary Fig. 1B.

CSD-related persistent changes in the ongoing activity of meningeal afferents

(A) Example heatmap of standardized ongoing activity (fraction of time afferents exhibited calcium events when the mouse is not locomoting) for all afferent fibers from a single FOV during baseline and up to 120 min following CSD. Data shows concatenated 1 min bins of activity. Afferents were either activated, suppressed, or unaffected by CSD. Note the delayed activation and immediate suppression in two small subsets of fibers. (B) Mean activity time course of the activated and suppressed afferents from the same population depicted in A. (C) Pie chart depicting the breakdown of the afferent subpopulations based on their change in ongoing activity following CSD. Most afferents were not affected (orange), while two smaller populations either exhibited prolonged activation (maroon) or suppression (blue) of ongoing activity following CSD (n=8 FOVs). (D) Afferents exhibiting prolonged activation had a longer onset latency than those exhibiting suppression (**** P < 0.0001, Mann Whitney t-test. Error bars: SEM). (E) The duration of increases in ongoing activity and suppressions in activity were similar (P = 0.97, two-tailed t-test. Error bars: SEM).

Locomotion and related meningeal deformations changes post CSD

(A) In head-fixed mice, wheel running drives meningeal scaling, shearing, and positive Z-shifts (i.e., meningeal movement toward the skull). (B, C) When compared to the baseline period, CSD was associated with an overall decrease in wheel running (** P < 0.01, paired t-test, n = 9) and locomotion bout rate (* P < 0.05, Wilcoxon, signed-rank test). (D, E) CSD, however, did not affect bout duration (P = 0.50, paired t-test) or bout peak velocity (P = 0.18, paired t-test). (F, G, H). CSD also did not affect subsequent locomotion bout-evoked meningeal scaling, shearing, or Z-shift (P = 0.56; P = 0.55, P = 0.18, paired t-tests, respectively, n = 9 for scale and shear, n = 7 for Z-shift). Bars depict the mean.

CSD leads to sensitization of meningeal afferents to local deformation signals

(A) Example of meningeal afferent sensitization following CSD. Locomotion and its related Z-shift (bottom traces) are comparable before (left) and after (right) CSD, but afferent fibers exhibit greater responses after CSD (heatmaps, top panels). (B) Example GLM fit of afferent activity before CSD. A raw calcium activity trace recorded pre-CSD (gray traces, z-scored; σ: 1 standard deviation) is plotted along with the model fit (purple). The deviance explained (‘dev exp’) is a metric of GLM fit quality and is above the threshold (0.05) for classifying an afferent’s activity as reasonably well-fit by the GLM. (C) GLM β coefficients used as a metric of the coupling between the Z-shift and the activity of the example afferent shown in B across different delays. A maximal coefficient at zero delay indicates the alignment of activity with Z-shifts. Note the greater afferent activation per unit Z-shift after CSD relative to baseline, indicative of an augmented or sensitized response. (D) Pie chart indicating the numbers and distribution of all afferents that were well-fit by deformation and/or locomotion signals either before and/or after CSD. Afferents categorized as sensitized if they (i) had significant GLM fits both pre- and post-CSD, as well as higher coefficients for a given deformation and/or locomotion predictor post-CSD (purple) or (ii) were well-fit only post-CSD (magenta). Two small subsets of afferents categorized as desensitized had worse GLM fits post-CSD (mustard) or were no longer well-fit post-CSD (orange). The incidence of afferent sensitization exceeded that of desensitization (P < 0.001, χ2 test). (E, F) Comparisons of pre- and post-CSD GLM coefficients for each of the deformation and locomotion predictors. Data are shown for sensitized afferents with well-fit models pre- and post-CSD (corresponding to the purple population in D) and for afferents with well-fit models only post-CSD (i.e., silent pre-CSD, corresponding to the magenta population in D). Mouse velocity coefficients were close to 0 in all cases and are not shown. In the two sensitized afferent populations, only coefficients related to deformation predictors increased post-CSD (** p< 0.01, *** p< 0.001, **** p< 0.0001, Wilcoxon sign rank test with correction for multiple analyses). (G, H) The response bias of sensitized afferents to meningeal deformation was further observed when comparing these GLMs to restricted GLMs that included only the group of deformation predictors or the group of locomotion predictors. The deviance explained by the deformation response component (estimated as the difference between the full GLM and the GLM lacking deformation variables) was significantly greater than for the locomotion response component in sensitized afferents that were well-fit pre- and post-CSD and for those that were well-fit only post-CSD (*** P < 0.001 and **** P < 0.0001, Wilcoxon test for G and H respectively). Bars depict mean; error bars indicate SEM. (I) Among the sensitized afferents with enhanced sensitivity to deformation variables, we observed a similar sensitization to scale, shear, and Z-shift variables. Bars depict mean; error bars indicate SEM. (J) There was no difference in the incidence of sensitized afferents among afferents that showed prolonged activation, prolonged suppression, or no change in ongoing activity post-CSD (P = 0.9, χ2 test; cf. Figure 2).

Further analyses of hemodynamics and calcium activity related to the CSD wave

(A) Example cropped subregion of a field-of-view (FOV) demonstrating a typical pial vascular response to CSD, including an initial constriction (middle) followed by dilation (right). (B) Analysis used for detecting the CSD-associated meningeal calcium wave. Top: Time derivative of the FOV-averaged fluorescence signal over the first 90 s of imaging after the cortical pinprick. An initial period was defined between the derivative’s maximum (orange circle) and minimum (orange X) and was used to define the baseline fluorescence value (Fpre). Bottom: Normalized FOV fluorescence signal. The maximum within the initial period (purple circle) is taken as the peak of the CSD-associated wave, and a threshold value (dashed line) is defined as 10% of that peak value. The threshold crossings before (blue circle) and after (blue x) the peak define the wave’s onset and offset, respectively.