# Elevated pyramidal cell firing orchestrates arteriolar vasoconstriction through COX-2-derived prostaglandin E2 signaling

Benjamin Le Gac, Marine Tournissac, Esther Belzic, Sandrine Picaud, Isabelle Dusart, Hédi Soula, Dongdong Li, Serge Charpak, Bruno Cauli

► To cite this version:

Benjamin Le Gac, Marine Tournissac, Esther Belzic, Sandrine Picaud, Isabelle Dusart, et al.. Elevated pyramidal cell firing orchestrates arteriolar vasoconstriction through COX-2-derived prostaglandin E2 signaling. eLife, 2025, 13, 10.7554/eLife.102424. [hal-05349112](https://hal.science/hal-05349112v1)

HAL Id: hal-05349112

https://hal.science/hal-05349112v1

Submitted on 5 Nov 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Elevated pyramidal cell firing orchestrates arteriolar vasoconstriction through COX-2-derived prostaglandin E2 signaling

Benjamin Le Gac¹,², Marine Tournissac³, Esther Belzic¹,², Sandrine Picaud¹,², Isabelle Dusart¹,², Hédi Soula⁴, Dongdong Li¹,², Serge Charpak³, Bruno Cauli¹,²*

¹Sorbonne Université, CNRS, Inserm, Neuro-SU, Paris, France; ²Sorbonne Université, CNRS, Inserm, Institut de Biologie Paris-Seine, Paris, France; ³Sorbonne Université, CNRS, Inserm, Institut de la Vision, F-75012, Paris, France; ⁴Sorbonne Université, INSERM, Nutrition and Obesities: Systemic Approaches, NutriOmics, Research Unit, Paris, France

## eLife Assessment

This study presents **important** findings on the role of pyramidal cells driving vasoconstriction in brain arteries through a COX-2/PGE2 pathway, with additional contributions from NPY (interneurons) and 20-HETE (astrocytes). Optogenetic stimulation of cortical pyramidal neurons induces vasoconstriction, potentially leading to oxygen and nutrient undersupply in regions with sustained activation - a mechanism potentially relevant under pathological conditions. The authors provide **convincing** evidence from brain slice experiments and some in vivo data from anesthetized animals, carefully discussing the strengths and limitations of both approaches.

## Abstract
Neurovascular coupling, linking neuronal activity to cerebral blood flow, is essential for brain function and underpins functional brain imaging. Whereas mechanisms involved in vasodilation are well-documented, those controlling vasoconstriction remain overlooked. This study unravels the mechanisms by which pyramidal cells elicit arteriole vasoconstriction. Using patch-clamp recording, vascular and Ca²⁺ imaging in mouse cortical slices, we show that strong optogenetic activation of layer II/III pyramidal cells induces vasoconstriction, correlating with firing frequency and somatic Ca²⁺ increase. Ex vivo and in vivo pharmacological investigations indicate that this vasoconstriction predominantly recruits prostaglandin E2 through the cyclooxygenase-2 pathway, and activation of EP1 and EP3 receptors. We also present evidence that specific interneurons releasing neuropeptide Y, and astrocytes, through 20-hydroxyeicosatetraenoic acid, contribute to this process. By revealing the mechanisms by which pyramidal cells lead to vasoconstriction, our findings shed light on the complex regulation of neurovascular coupling.

## Introduction

The brain critically depends on the uninterrupted blood supply provided by a dense vasculature (Schmid et al., 2019). Cerebral blood flow (CBF) is locally and temporally controlled by neuronal activity, by an essential process called neurovascular coupling (NVC), and is impaired in early stages of numerous neurological disorders (Iadecola, 2017). NVC also serves as the physiological basis for

<footnote>*For correspondence: bruno.cauli@upmc.fr</footnote>

<footnote>**Competing interest:** The authors declare that no competing interests exist.</footnote>

<footnote>**Funding:** See page 21</footnote>

<footnote>**Preprint posted**
13 August 2024</footnote>

<footnote>**Sent for Review**
10 October 2024</footnote>

<footnote>**Reviewed preprint posted**
13 December 2024</footnote>

<footnote>**Reviewed preprint revised**
06 March 2025</footnote>

<footnote>**Version of Record published**
23 April 2025</footnote>

<footnote>**Reviewing Editor:** Hana Uhlirova, Faculty of Mechanical Engineering, Brno University of Technology and Institute of Physical Engineering, Czech Republic</footnote>

<footnote>© Copyright Le Gac et al. This article is distributed under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use and redistribution provided that the original author and source are credited.</footnote>
functional brain imaging widely used to map neuronal activity. Neuronal activity increases CBF within seconds (Iadecola, 2017). In the cerebral cortex, the hyperemic response linked to neural activity is supported by dynamically controlled vasodilation that is spatially and temporally constrained by vaso-constriction in a second phase (Devor et al., 2007). Conversely, vasoconstriction and decreased CBF usually correlate with reduced neuronal activity (Devor et al., 2007; Shmuel et al., 2002).

Mounting evidence indicates that the positive correlation between neuronal activity and CBF is not always maintained under physiological conditions: (i) robust sensory-evoked vasodilation can occur in the absence of substantial neuronal response (O'Herron et al., 2016), (ii) conversely, pronounced neuronal activity is not systematically associated with increased hemodynamics (Ma et al., 2016), (iii) CBF is decreased in several cortical areas despite local increase in neuronal activity (Devor et al., 2008), and (iv) optogenetic simulation of inhibitory GABAergic interneurons results in vasodilation (Uhlirova et al., 2016). Furthermore, in pathological conditions with intense neuronal activity such as epileptic seizures, a sustained hypoperfusion induced by vasoconstriction is observed (Farrell et al., 2016; Tran et al., 2020).

NVC is achieved by the synthesis and release of vasoactive messengers within the neurovascular unit (Iadecola, 2017), which act on the contractility of mural cells (smooth muscle cells and pericytes) to control vessel caliber and CBF along the vascular tree (Rungta et al., 2018). Pial and penetrating arterioles, which have a higher density of contractile mural cells and control their diameter faster than capillaries (Hartmann et al., 2021; Hill et al., 2015; Rungta et al., 2021; Rungta et al., 2018), play a key role in regulating CBF.

Different experimental approaches, each with their advantages and limitations, have allowed the identification of several mediators of NVC (Grutzendler and Nedergaard, 2019; Iadecola and Nedergaard, 2007). Ex vivo brain slices provide a well-controlled environment, ideal for pharmacological investigations to dissect the underlying mechanisms. However, they lack connectivity and blood flow which provides both vascular tone and natural oxygenation and therefore require in vivo validation. Conversely, pharmacological studies are more challenging with in vivo preparations. Awake animals allow physiologically relevant context with largely undisturbed network and neuromodulatory activity. However, this preparation is subject to brain state changes which may affect network activity, metabolism, and vascular physiology (Grutzendler and Nedergaard, 2019), potentially complexifying the analysis of specific mechanisms of NVC. Although chronic anesthetized animals have reduced network and neuromodulation activity, the NVC response is only slowed (Rungta et al., 2021), providing a valuable model for validating ex vivo observations.

Messengers of vasodilation released by excitatory neurons, GABAergic interneurons, astrocytes, or endothelial cells, include nitric oxide, K+, arachidonic acid derivatives such as prostaglandin E2 (PGE2; Iadecola, 2017), or more recently glutamate (Zhang et al., 2024). Despite its physio pathological importance, vasoconstriction is less understood with fewer cell types and vasoactive messengers that have been identified. It is now generally accepted that GABAergic interneurons are key players in vasoconstriction by releasing neuropeptide Y (NPY; Cauli et al., 2004; Uhlirova et al., 2016). Under certain conditions, astrocytes can also induce vasoconstriction via 20-hydroxyeicosatetraenoic acid (20-HETE) (Mulligan and MacVicar, 2004) or high K+ concentration (Girouard et al., 2010). However, the involvement of pyramidal cells in vasoconstriction has been overlooked.

PGE2 has emerged as a bimodal messenger of NVC, similar to K+ (Girouard et al., 2010) and glutamate (Zhang et al., 2024), that can induce either vasodilation (Gordon et al., 2008; Lacroix et al., 2015; Lecrux et al., 2011; Mishra et al., 2016) or vasoconstriction (Dabertrand et al., 2013; Rosehart et al., 2021) depending on its concentration and/or site of action along the vascular tree. Under physiological conditions, PGE2 is produced during NVC by either astrocytes (Mishra et al., 2016) or pyramidal cells (Lacroix et al., 2015) via the rate-limiting synthesizing enzymes cyclooxygenase-1 (COX-1) or -2 (COX-2), respectively. Since COX-2-expressing pyramidal cells can release glutamate and PGE2, both of which induce vasoconstriction at high concentrations (Dabertrand et al., 2013; Rosehart et al., 2021; Zhang et al., 2024), pyramidal cells may be responsible for vasoconstriction when their spiking activity is high.

To test this hypothesis, we used ex vivo and in vivo approaches in combination with optogenetics to precisely control pyramidal cell firing in the mouse barrel cortex while monitoring the resulting arborolar response. We found that pyramidal cells induce vasoconstriction at high stimulation frequency and about half of them express all the transcripts required for a cell autonomous synthesis of the
vasoconstrictor messengers PGE2 and prostaglandin F2α (PGF2α). Pharmacological investigations
revealed that this neurogenic vasoconstriction depends on COX-2-derived PGE2 via the direct acti-
vation of vascular EP1 and EP3 receptors. It also involves the recruitment of intermediary NPY inter-
neurons acting on the Y1 receptor, and, to a lesser extent astrocytes, via 20-HETE and COX-1-derived
PGE2. Thus, our study reveals the mechanisms by which high-frequency pyramidal cell firing leads to
vasoconstriction.

Results

Pyramidal cells induce vasoconstriction at high firing frequency

To determine if pyramidal cells action potential (AP) firing can induce vasoconstriction in a frequency-
dependent manner, we used optogenetics to induce AP firing while monitoring the resulting vascular
response in cortical slices. We used Emx1-cre;Ai32 transgenic mice expressing the H134R variant of
channelrhodopsin-2 (ChR2) in the cortical glutamatergic neurons (Gorski et al., 2002), conferring
robust pyramidal cell photoexcitability (Madisen et al., 2012). Wide-field photostimulation of cortical
slices was achieved in layers I to III (Figure 1—figure supplement 1A) using 10 s trains of 5ms light
pulses (see Materials and methods) delivered at five different frequencies (1, 2, 5, 10 and 20 Hz,
Figure 1A).

First, we ensured the efficiency of the photostimulation paradigm by recording layer II-III pyra-
midal cells in whole-cell current clamp mode (Figure 1A). We observed that optogenetic stimulation
resulted in the firing of an initial AP that was followed by a train of spikes whose amplitude and
frequency transiently decreased before reaching a steady state (Figure 1A, upper traces). Consistent
with the kinetic properties of the H134R ChR2 variant (Lin et al., 2009) and the intrinsic firing prop-
erties of pyramidal cells (Karagiannis et al., 2009), the steady-state firing frequency matched the
photostimulation frequency up to 5 Hz but was lower at higher frequencies (Figure 1A, steady-state
spike success rate: 100 ± 0% at 1, 2, and 5 Hz, 70 ± 11% at 10 Hz and 55 ± 12% at 20 Hz). These
observations demonstrate efficient pyramidal cell activation over a wide range of photostimulation
frequencies.

To test the hypothesis that neuronal activity induces vasoconstriction, we analyzed the optoge-
netically induced response of penetrating arterioles. Layer I arterioles were imaged for 30 min in
cortical slices (Figure 1—figure supplement 2; Table 1) without preconstriction to facilitate obser-
vation of vasoconstriction (Cauli et al., 2004). Examination of the evoked vascular response over
30 min (Figure 1—figure supplement 2) showed that increasing the frequency of photostimulation
shifted the overall vascular response from a barely discernible delayed response between 1 Hz and
5 Hz to a sustained vasoconstriction at 10 Hz and above which began less than 2 min after photo-
stimulation (10 Hz: 1.4±0.4 min; 20 Hz: 1.6±0.5 min). Most vessels (n=8 of 10 arterioles) showed a
strong and rapid vasoconstriction at 20 Hz. On average, this response peaked at 6.8±2.4 min, much
earlier than at lower frequencies, which typically required more than 10 min to reach a maximum
(Figure 1—figure supplement 2C, 1 Hz: 15.6±4.0 min; 2 Hz: 13.2±2.3 min; 5 Hz: 16.0±3.6 min; 10 Hz:
15.7±2.4 min). Because the vascular response shifted to reliable vasoconstriction, with onset and peak
in less than 2 and 10 min, respectively, similar to previous observations in cortical slices (Cauli et al.,
2004), when the frequency of photostimulation was increased to 20 Hz, we defined the first 10 min
of recording as the vasoconstriction time frame for subsequent comparisons and analyses. While
photostimulation at 1–5 Hz failed to elicit fast reliable vascular responses (Figure 1C and D), 10 Hz
photostimulation predominantly induced vasoconstriction (n=4 of 5 arterioles, Figure 1C and D, area
under the curve (AUC)=-1.7 ± 1.1 x 10³ %.s, n=5). This response was even more pronounced at 20 Hz,
as all arterioles showed vasoconstriction of high magnitude (Figure 1B-D; AUC = -3.7 ± 0.7 x 10³ %.s,
F(4, 30)=6.135, p=9.89 x 10⁻⁴, one-way ANOVA, n=10 arterioles). This difference was particularly striking
when comparing the magnitude at 20 Hz (Figure 1D) with those at 1 Hz (t(12) = -3.48, p=0.0407,
t-test), 2 Hz (t(18) = -4.09, p=0.0250, t-test) and 5 Hz (t(14) = -3.7, p=0.0346, t-test). Intense optogenetic
stimulation of pyramidal cells has been shown to elicit cortical spreading depression (Chung et al.,
2019; Pham et al., 2024), which induces vasoconstriction (Zhang et al., 2024) and fast cell swelling
(Zhou et al., 2010). We ruled out this possibility by showing that the rate of change in light trans-
mittance associated with cell swelling remained below that of cortical spreading depression (Zhou
et al., 2010; Table 1). On the other hand, ChR2-independent vascular changes induced by high light
<figure><img src="image_3.png" /><figcaption>**Figure 1.** The occurrence and strength of vasoconstriction depends on the photostimulation frequency of pyramidal cells. (A) Representative examples of the voltage responses of a layer II-III pyramidal cell (upper traces light grey to black traces) induced by photostimulations (470 nm, 10 s train, 5ms pulses) delivered at 1, 2, 5, 10, and 20 Hz (cyan lower traces) and mean spike success rate (middle trace, n=4 cells from 3 mice). The SEMs envelope</figcaption></figure>

Figure 1 continued on next page
Figure 1 continued

the mean traces. The red dashed lines represent a spike success rate of 100%. (B) Representative example showing IR-DGC pictures of a layer I penetrating arteriole (1) before a 20 Hz photostimulation, (2) at the maximal diameter decrease, and (3) after 10 min of recording. Pial surface is upward. Yellow calipers represent the measured diameters. White dashed lines indicate the initial position of the vessel wall. Scale bar: 25 µm. (C) Kinetics of arteriolar diameter changes induced by photostimulation (vertical cyan bars) at 1 Hz (n=4 arterioles from 3 mice), 2 Hz (n=10 arterioles from 8 mice), 5 Hz (n=6 arterioles from 6 mice), 10 Hz (n=5 arterioles from 5 mice), and 20 Hz (n=10 arteriole from 9 mice). The SEMs envelope the mean traces. The blue trace represents the kinetics of the diameter changes of the arteriole shown in (B). (D) Effects of the different photostimulation frequencies on AUC of vascular responses during 10 min of recording. Data are presented as the individual values and mean ± SEM. * statistically different from 20 Hz stimulation with p<0.05.

The online version of this article includes the following source data and figure supplement(s) for figure 1:

**Source data 1.** Detection of spikes per light pulse interval from different cells used to determine spike success rate in **Figure 1A**.

**Source data 2.** Diameter measurements (µm) of individual arterioles used to determine diameter changes in **Figure 1C** and **Figure 1—figure supplement 2**.

**Figure supplement 1.** Vasoconstriction induced by widefield photostimulation is specific of ChR2 expression in pyramidal cells.

**Figure supplement 1—source data 1.** Diameter measurements (µm) of individual arterioles used to determine diameter changes in **Figure 1—figure supplement 1C**.

**Figure supplement 2.** Vasoconstrictions occurred during the 30 first minutes after pyramidal cells photoactivation.

intensity have been reported (Rungta et al., 2017). We verified that 20 Hz photostimulation did not induce a vascular response in wild-type mice that do not express ChR2 (Figure 1—figure supplement 1). Taken together, our observations indicate that photostimulation of pyramidal cells produces a frequency-dependent vasoconstriction.

<table><caption>**Table 1.** Morphological and physiological properties, and neurovascular responses of diving arterioles used in the analysis of the frequency-dependence of the polarity of neurovascular response evoked by pyramidal cells.</caption><thead><tr><th>Frequency</th><th>1 Hz</th><th>2 Hz</th><th>5 Hz</th><th>10 Hz</th><th>20 Hz</th></tr></thead><tbody><tr><td>Number of arterioles</td><td>n=4</td><td>n=10</td><td>n=6</td><td>n=5</td><td>n=10</td></tr><tr><td></td><td>1.2±0.2</td><td>1.6±0.2</td><td>1.6±0.2</td><td>1.3±0.2</td><td>1.0±0.1</td></tr><tr><td></td><td colspan="5">F (4, 30)=2.161<br>p=0.098</td></tr><tr><td>Resting stability (%)</td><td>n.s.</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>3.6±0.8</td><td>3.8±0.3</td><td>3.2±0.5</td><td>4.1±0.8</td><td>4.0±0.2</td></tr><tr><td></td><td colspan="5">F (4, 30)=0.656<br>p=0.627</td></tr><tr><td>Wall thickness (µm)</td><td>n.s.</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>0.5±0.2</td><td>0.0±0.5</td><td>0.2±1</td><td>-1.7±1.1</td><td>-3.7±0.7</td></tr><tr><td></td><td colspan="5">F (4, 30)=6.135<br>p=0.00099</td></tr><tr><td>Area under the curve after photostimulation (AUC; x10³ %s)</td><td>***</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>0.19±0.09</td><td>0.64±0.14</td><td>0.22±0.04</td><td>0.43±0.19</td><td>0.8±0.11</td></tr><tr><td>Maximal dΔT/dt (%s<sup>-1</sup>)</td><td colspan="5">All &lt;2 %s<sup>-1</sup></td></tr></tbody></table>

Data are mean ± SEM, one-way ANOVA F test and corresponding exact p-value. n.s., not statistically different and ***: p<0.001.

The online version of this article includes the following source data for table 1:

**Source data 1.** Properties of individual arterioles used for **Table 1**.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Photostimulation of pyramidal cells elicits a time-locked firing and a frequency-dependent calcium increase. (A) Voltage response (top trace) and kinetics of relative fluorescence changes (red bottom trace) induced by photostimulation at 20 Hz. Insets, IR-DGC (top), Rhod2 fluorescence (bottom) pictures of an imaged layer II/III pyramidal cell. The somatic region of interest is outlined in white. Pial surface is upward. Scale bar: 20 µm. (B) Mean relative variations of Ca²⁺ fluorescence in response to photostimulation at 2 Hz (grey, n=9 cells from 5 mice) and 20 Hz (black, n=9 cells from 5 mice). Dashed line represents the baseline. The vertical cyan bar indicates the duration of photostimulation. SEMs envelope the mean traces. Inset, Maximum increase in relative fluorescence changes induced immediately after photostimulation, indicated by the black arrow. The data are shown as the individual values and mean ± SEM. * statistically different with p<0.05.</figcaption></figure>

The online version of this article includes the following source data and figure supplement(s) for figure 2:

**Source data 1.** Somatic fluorescence measurements (A.U.) used to determine fluorescence changes in **Figure 2**.

**Figure supplement 1.** Photostimulation of pyramidal cells does not evoke recurrent spiking network activity.

**Figure supplement 1—source data 1.** Detection of spikes per second used to determine the mean firing frequency in **Figure 2—figure supplement 1C and D**.

## Optogenetic stimulation induces a frequency-dependent, gradual increase in somatic calcium that precedes the vascular response

These observations raise the questions of how pyramidal neurons can induce vasoconstriction at higher AP-firing rates. It is generally accepted that the synthesis and/or release of vasodilatory substances requires an increase in intracellular Ca²⁺ in the releasing cells (Attwell et al., 2010; Cauli and Hamel, 2010), but little is known about the release of vasoconstricting substances. We therefore determined whether an increase in somatic Ca²⁺ concentration in cortical neurons was also dependent on photostimulation frequency. We combined optogenetic stimulation with whole-cell current clamp recording and intracellular Ca²⁺ imaging using Rhod-2 delivered by patch pipette (Figure 2A). Excitation of this red Ca²⁺ indicator at 585 nm did not induce any voltage response in the recorded pyramidal cells (Figure 2A), as expected from the action spectrum of ChR2 (Lin et al., 2009). In contrast, photostimulation at 470 nm elicited a train of spikes accompanied by a somatic Ca²⁺ increase that decayed for tens of seconds after photostimulation without triggering any significant recurrent spiking activity (Figure 2—figure supplement 1). The Ca²⁺ response evoked by 20 Hz photostimulation was more than twice of that evoked by 2 Hz photostimulation (Figure 2B; 2 Hz: ΔF/F₀=34.5 ± 3.7 %, n=9 cells, vs. 20 Hz:ΔF/F₀=79.5 ± 17.7 %, n=9 cells; t (16)=2.485, p=0.024397), while the average number of evoked spikes was about five times higher (Figure 1A, Figure 2—figure supplement 1). These results demonstrate a frequency-dependent increase in intracellular Ca²⁺ induced by photostimulation, that precedes vasoconstriction. We therefore aimed to understand the molecular mechanisms linking neuronal activity to vasoconstriction.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Optogenetically-induced vasoconstriction requires AP firing and partially glutamatergic transmission. Effect of TTX (1 µM, brown, n=6 arterioles from 5 mice) and cocktail antagonists of AMPA/kainate (DNQX, 10 µM), NMDA (D-AP5, 50 µM), mGluR1 (LY367385, 100 µM) and mGluR5 (MPEP, 50 µM) receptors (gray, n=10 arterioles from 6 mice) on (A) kinetics and (B) magnitude of arteriolar vasoconstriction induced by 20 Hz photostimulation (cyan bar). The SEMs envelope the mean traces. Dashed lines represent the initial diameter. The shaded traces correspond to the kinetics of arteriolar vasoconstriction in control condition (Figure 1C – 20 Hz). Data are presented as the individual values and mean ± SEM. * and *** statistically different from control condition (Figure 1D – 20 Hz) with p<0.05 and p<0.001, respectively.</figcaption></figure>

The online version of this article includes the following source data and figure supplement(s) for figure 3:

**Source data 1.** Diameter measurements (µm) of individual arterioles used to determine diameter changes in **Figure 3A**.

**Figure supplement 1.** Basal network activity and tonic glutamate do not influence resting vascular tone.

**Figure supplement 1—source data 1.** Diameter measurements (µm) of individual arterioles used to determine diameter changes in **Figure 3—figure supplement 1**.

## Vasoconstriction induced by pyramidal cells requires AP firing and is partially dependent on glutamatergic transmission

In pyramidal cells, APs induce both somatic Ca²⁺ elevation (Smetters et al., 1999) and glutamate release. To determine whether spiking activity is required for vasoconstriction induced by 20 Hz photostimulation, we blocked APs with the voltage-activated sodium channel blocker tetrodotoxin (TTX, 1 µM, n=6 arterioles). This treatment completely abolished the vasoconstriction evoked by 20 Hz photostimulation (Figure 3; AUC = 0.4 ± 0.4 x 10³ %·s, t(14) = 5.57, p=8.6656 x 10⁻⁶). These data indicate that APs are mandatory for neurogenic vascular response and may involve glutamate release.

Indeed, high levels of glutamate released from pyramidal cells may activate NPY-expressing interneurons or astrocytes through activation of ionotropic or group I metabotropic glutamate receptors (Girouard et al., 2010; Mulligan and MacVicar, 2004; Uhlirova et al., 2016). It may also directly activate NMDA receptors on arteriolar smooth muscle cells, resulting in a large intracellular Ca²⁺ increase and subsequent vasoconstriction (Zhang et al., 2024). To test the hypothesis that glutamate from pyramidal cells, either directly or indirectly, results in vasoconstriction, we blocked glutamatergic transmission by antagonizing AMPA/kainate, NMDA and group I metabotropic receptors expressed by cortical neurons (Tasic et al., 2016; Zeisel et al., 2015) and juvenile astrocytes (Sun et al., 2013). Glutamate receptor antagonists reduced the magnitude of vasoconstriction (–1.6±0.4 x 10³ %·s, t(18) = 3.28, p=0.0160) by approximately half (Figure 3). Taken together, our data suggest that photostimulation of pyramidal cells elicits a frequency-dependent vasoconstriction that requires AP firing and partially involves glutamatergic transmission. We therefore sought to elucidate the glutamate-independent vasoactive pathway underlying this neurogenic vascular response.

## Pyramidal cells express the mRNAs for a cell autonomous PGE2 and PGF2α synthesis

Several arachidonic acid metabolites produced after intracellular Ca²⁺ elevation, including PGF2α, but also PGE2, exert dose-dependent vasoconstrictive effects (Dabertrand et al., 2013; Rosehart et al., 2021; Zonta et al., 2003). These prostaglandins could therefore be progressively released as the frequency of photostimulation and somatic Ca²⁺ increase and thereby promote vasoconstriction. Layer II-III pyramidal cells have been shown to produce PGE2 (Lacroix et al., 2015). To determine whether
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Layer II-III pyramidal cells express PGE2 and PGF2α synthesizing enzymes. (A) Voltage responses of a layer II-III pyramidal cell induced by injection of current (bottom traces). In response to a just-above-threshold current pulse, the neuron fired long-lasting action potentials with little frequency adaptation (middle black trace). Near saturation, it exhibits the pronounced spike amplitude accommodation and marked frequency adaptation characteristic of regular spiking cells (upper grey trace). (B) Agarose gel analysis of the scRT-PCR products of the pyramidal cell shown in (A) revealing expression of vGluT1, COX-2, mPGES2, cPGES, PM-PGFS and CBR1. Φ x174 digested by *Hae*III (Φ in bps) was used as molecular weight marker (C) Histogram summarizing the single-cell detection rate of PGE2 and PGf2α synthesizing enzymes in layer II-III pyramidal cells (n=16 cells from 6 mice). PGES (green zone) corresponds to mPGES1, mPGES2 and/or cPGES and PGFS (blue zone) to PM-PGFS, CBR1 and/or AKR1B3. (D) Co-expression of PGE2 and PGf2α synthesizing enzymes in pyramidal cells. The box size is proportional to the detection rate. Note the absence of co-expression between COX-1 (purple) and COX-2 (red). Co-expression of a PGES (left, green) and a PGFS (right, blue) with COX-1 (up) and COX-2 (bottom).</figcaption></figure>

The online version of this article includes the following source data and figure supplement(s) for figure 4:

**Source data 1.** Original files of the full raw unedited gels shown in **Figure 4B**.

**Source data 2.** Uncropped gels shown in **Figure 4B** with relevant lanes labeled.

**Figure supplement 1.** Sensitivity of the RT-mPCR protocol.

**Figure supplement 1—source data 1.** Original file of the full raw unedited gel shown in **Figure 4—figure supplement 1**.

**Figure supplement 1—source data 2.** Uncropped gels shown in **Figure 4—figure supplement 1** with relevant lanes labeled.

the synthesizing enzymes of PGE2 and PGF2α are present in pyramidal cells, we performed single-cell RT-PCR after patch-clamp recording (Devienne et al., 2018). Sixteen layer II-III pyramidal cells were visually identified based on the triangular shape of their soma and a prominent apical dendrite. Their glutamatergic phenotype was confirmed both by their stereotypical regular spiking firing pattern (Figure 4A) and also by the expression of the vesicular glutamate transporter, vGlut1, and neither of the two GABA synthesizing enzymes, thus excluding possible contamination by GABAergic interneurons (Figure 4B; Karagiannis et al., 2009). The rate-limiting enzymes of prostaglandin synthesis, cyclooxygenase-1 (COX-1) and -2 (COX-2), were detected in 25% (n=4 of 16 cells) and 31% (n=5 of 16 cells) of pyramidal cells, respectively, (Figure 4B–D) but were never co-expressed (Figure 4D). Although the differential expression of COX-1 or COX-2 allowed the definition of three non-overlapping molecular subpopulations of pyramidal cells, they did not show distinctive electrophysiological features
(Table 2). The cytosolic enzyme responsible for synthesizing PGE2 (cPGES) was observed in most pyramidal cells (Figure 4B and C; 88%, n=14 of 16 cells). In addition, the microsomal PGES, mPGES1 and mPGES2 were detected in 6% (Figure 4C, n=1 of 16 cells) and 38% (Figure 4B and C; n=6 of 16 cells) of pyramidal cells, respectively, and were always co-expressed with cPGES. The PGF2α terminal-synthesizing enzyme AKR1B3 was observed in the majority of neurons (Figure 4B and C; n=11 of 16 cells, 69%). Occasionally, it was co-detected with the prostamide/prostaglandin F synthase (PM-PGFS, Figure 4C, n=5 of 16 cells, 31%) and the PGE2 converting enzyme carbonyl reductase 1 (CBR1, Figure 4C; n=3 of 16 cells, 19%). CBR1 was consistently detected alongside at least one PGES. Pyramidal cells positive for COX-1 also expressed PGES (Figure 4D), with most of them also co-expressing PGFS (n=3 out of 4). All neurons positive for COX-2 co-expressed both PGES and PGFS (Figure 4D). These molecular observations suggest that subpopulations accounting for about half of layer II-III pyramidal cells express all the transcripts necessary for the synthesis of PGE2 and PGF2α.

## Prostaglandins underpin vasoconstriction ex vivo and in vivo

To investigate whether prostaglandins could mediate neurogenic vasoconstriction, we inhibited their synthesis. In cortical slices, the non-selective COX inhibitor indomethacin (5 µM, Table 3, Figure 5A and B) completely abolished the vascular response (n=10 arterioles, AUC = 0.0 ± 0.2 x 10³ %.s, t₁₈ = 5.86, P=3.4385 x 10⁻⁵). To verify that high-frequency stimulation of pyramidal cells also induces vasoconstriction in vivo, 10 Hz photostimulation was reproduced in anesthetized Emx1-cre;Ai32 mice. Pial arterioles diameter measured by two-photon line-scan imaging (Figure 5C and D) revealed that pyramidal cells induce vasodilation (1st phase) followed by sustained vasoconstriction (2nd phase, Figure 5D and E). The constriction phase was inhibited by indomethacin (injected intravenously (i.v.)), indicating the involvement of prostaglandins in the vasoconstriction (AUC<sub>Ctrl</sub> = 291.5 ± 92.4 %.s vs. AUC<sub>Indo.</sub>=332.4 ± 184.4 %.s, U<sub>(5,4)</sub> = 0, p=0.0159, Figure 5D–F), which confirms our ex vivo observations. To determine whether they originated from COX-1 or COX-2 activity, we utilized selective inhibitors in cortical slices. The vasoconstriction magnitude was reduced by the COX-1 inhibitor SC-560 (100 nM, Table 3, n=10 arterioles, -1.4±0.7 x 10³ %.s, t<sub>(18)</sub> = 3.54, p=9.3396 x 10⁻³, Figure 5A and B). The COX-2 inhibitor NS-398 (10 µM, Table 3, n=7 arterioles) completely abolished pyramidal cell-induced vasoconstriction in a more potent manner (Figure 5A and B; AUC = 0.1 ± 0.3 x 10³ %.s, t<sub>(15)</sub> = 5.45, p=1.0853 x 10⁻⁵), mimicking the ex vivo effect of indomethacin. These observations suggest that prostaglandins, derived mainly from COX-2 activity, and to a lesser extent from COX-1 activity, mediate pyramidal cell-induced vasoconstriction.

## PGE2 mediates vasoconstriction by acting primarily on EP1 receptor

To determine the nature of the prostaglandins and their receptors, we selectively antagonized the vasoconstrictor receptors of PGE2, EP1 or EP3, or the FP receptor of PGF2α. The magnitude of vasoconstriction was reduced by the selective EP1 receptor antagonist ONO-8130 (10 nM, Table 3, n=9 arterioles, Figure 5G and H, 0.3±0.3 x 10³ %.s, t<sub>(17)</sub> = 6.01, p=2.8451 x 10⁻⁶), and to a lesser extent, by the EP3 receptor antagonist L-798,106 (1 µM, Table 3, n=9 arterioles, Figure 5G and H, -0.9±0.4 x 10³ %.s, t<sub>(17)</sub> = 4.30, p=8.0261 x 10⁻⁴). Impairing FP receptor signaling with AL-8810 (10 µM, Table 3, n=9 arterioles, Figure 5G and H) tended to reduce the evoked vasoconstriction, however, it did not reach statistical significance (AUC = -1.9 ± 0.4 x 10³ %.s, t<sub>(17)</sub> = 2.82, p=0.0533). Additionally, the preconstricted state induced by this weak partial FP agonist (Sharif and Klimko, 2019; Figure 5—figure supplement 1F) resulted in a diameter reduction of approximately 4% (diameter before application 21.8±2.5 µm vs. during 20.9±2.6 µm, n=9 arterioles, t<sub>(8)</sub> = 2.374, p=0.0457, paired t-test), which underestimated the optogenetic vascular response. Taken together, these results indicate that pyramidal cell photoactivation induces vasoconstriction through the release of PGE2 originating mainly from COX-2. This effect primarily acts on the EP1 receptor and, to a lesser extent, on the EP3 receptor.

To test a direct effect of PGE2 through vascular EP1/EP3 activation, we determined whether exogenous agonists of PGE2 receptors could mimic the vasoconstriction induced by pyramidal cell photostimulation. Similar to increasing photostimulation frequencies, exogenous application of PGE2 induced vasoconstriction in a dose-dependent manner which persisted for several minutes after removal (Figure 5—figure supplement 2A and D). Likewise, 10 µM sulprostone, an EP1/EP3 agonist with an EC₅₀ comparable to that of PGE2 (Boie et al., 1997), mimicked the vasoconstriction induced by 1–10 µM PGE2 (Figure 5—figure supplement 2). Application of 10 µM PGE2 in the presence of
<table><caption>**Table 2.** Electrophysiological properties of pyramidal cells recorded during single-cell RT-PCR experiments.</caption><thead><tr><th></th><th>COXs-negative<br>(n=7)</th><th>COX-1 positive<br>(n=4)</th><th>COX-2 positive<br>(n=5)</th></tr></thead><tbody><tr><td colspan="4"><strong>Passive properties</strong></td></tr><tr><td>Resting potential (mV)</td><td>-82.0±2.2</td><td>-84.7±4.2</td><td>-82.0±6.1</td></tr><tr><td>Input resistance (MΩ)</td><td>329±53.7</td><td>360.8±63.7</td><td>314.2±79.7</td></tr><tr><td>Time constant (ms)</td><td>50.7±6.4</td><td>47.3±7.2</td><td>47.74±10.1</td></tr><tr><td>Membrane capacitance (pF)</td><td>161.4±14.7</td><td>133.2±7.8</td><td>159.4±23.0</td></tr><tr><td>Sag index (%)</td><td>11.3±3.8</td><td>6.7±1.1</td><td>6.9±1.8</td></tr><tr><td colspan="4"><strong>Just above threshold properties</strong></td></tr><tr><td>Rheobase (pA)</td><td>52.7±8.7</td><td>51.7±15.5</td><td>62.3±18.2</td></tr><tr><td>First spike latency (ms)</td><td>295.2±50.8</td><td>271.6±62.9</td><td>178.7±55.6</td></tr><tr><td>Adaptation (Hz/s)</td><td>-3.1±0.9</td><td>-2.6±0.3</td><td>-3.3±1.3</td></tr><tr><td>Minimal frequency (Hz)</td><td>5.5±0.7</td><td>4.6±0.4</td><td>6.3±1.2</td></tr><tr><td colspan="4"><strong>Firing properties</strong></td></tr><tr><td>Accommodation (mV)</td><td>16.4±4.6</td><td>24.8±4.9</td><td>11.9±3.7</td></tr><tr><td>Amplitude of early adaptation (Hz)</td><td>62.1±13.3</td><td>89.8±8.5</td><td>62.3±16.4</td></tr><tr><td>Time constant of early adaptation (ms)</td><td>29.7±4.4</td><td>28.3±2.2</td><td>43.1±18.7</td></tr><tr><td>Late adaptation (Hz/s)</td><td>-11.6±2.3</td><td>-9.6±3.9</td><td>-10.7±1.7</td></tr><tr><td>Maximal frequency (Hz)</td><td>23.4±1.7</td><td>27.6±3.3</td><td>22.0±2.9</td></tr><tr><td colspan="4"><strong>Action potentials properties</strong></td></tr><tr><td>1<sup>st</sup> spike amplitude (mV)</td><td>94.8±1.7</td><td>93.0±5.0</td><td>88.6±1.9</td></tr><tr><td>1<sup>st</sup> spike duration (ms)</td><td>1.8±0.1</td><td>1.9±0.1</td><td>1.8±0.1</td></tr><tr><td>2<sup>nd</sup> spike amplitude (mV)</td><td>91.6±1.9</td><td>91.9±4.5</td><td>85.5±2.5</td></tr><tr><td>2<sup>nd</sup> spike duration (ms)</td><td>1.9±0.1</td><td>1.9±0.1</td><td>1.9±0.1</td></tr><tr><td>Amplitude Reduction (%)</td><td>3.4±0.4</td><td>1.1±0.7</td><td>3.6±1.4</td></tr><tr><td>Duration Increase (%)</td><td>5.3±0.9</td><td>3.8±1.5</td><td>7.5±1.6</td></tr><tr><td colspan="4"><strong>AHP and ADP properties</strong></td></tr><tr><td>1<sup>st</sup> spike fast AHP (mV)</td><td>-8.7±0.8</td><td>-8.7±0.6</td><td>-7.9±1.1</td></tr><tr><td>1<sup>st</sup> spike ADP (mV)</td><td>0.2±0.1</td><td>0.2±0.2</td><td>0±0</td></tr><tr><td>1<sup>st</sup> spike medium AHP (mV)</td><td>-13.8±1.4</td><td>-15.5±0.6</td><td>-13.7±0.8</td></tr><tr><td>1<sup>st</sup> spike fast AHP latency (ms)</td><td>8.2±1.0</td><td>7.7±1.4</td><td>9.8±1.5</td></tr><tr><td>1<sup>st</sup> spike ADP latency (ms)</td><td>4.6±2.2</td><td>2.0±2.0</td><td>0±0</td></tr><tr><td>1<sup>st</sup> spike, medium AHP latency (ms)</td><td>49.2±3.3</td><td>48.2±5.4</td><td>54.8±6.1</td></tr><tr><td>2<sup>nd</sup> spike fast AHP (mV)</td><td>-9.8±1.1</td><td>-8.5±0.5</td><td>-9.6±1.0</td></tr><tr><td>2<sup>nd</sup> spike ADP (mV)</td><td>0±0</td><td>0.1±0.1</td><td>0±0</td></tr><tr><td>2<sup>nd</sup> spike medium AHP (mV)</td><td>-16.1±1.2</td><td>-17.1±0.6</td><td>-15.8±0.4</td></tr><tr><td>2<sup>nd</sup> spike, fast AHP latency (ms)</td><td>8.6±0.9</td><td>7.1±0.7</td><td>11.7±1.3</td></tr><tr><td colspan="4">F (2.13)=4.063<br>p=0.0426</td></tr><tr><td colspan="4">*</td></tr><tr><td colspan="4">No significant difference in multiple comparisons.</td></tr></tbody></table>

Table 2 continued on next page
Table 2 continued

<table><thead><tr><th></th><th>COXs-negative<br>(n=7)</th><th>COX-1 positive<br>(n=4)</th><th>COX-2 positive<br>(n=5)</th></tr></thead><tbody><tr><td>2<sup>nd</sup> spike ADP latency (ms)</td><td>1.2±1.2</td><td>1.7±1.7</td><td>0±0</td></tr><tr><td>2<sup>nd</sup> spike, medium AHP latency (ms)</td><td>57.1±5.8</td><td>51.9±4.6</td><td>58.3±6.7</td></tr></tbody></table>

Data are presented as mean ± SEM. Statistical analyses were performed using a one-way ANOVA (F test) or a Kruskal-Wallis test, depending on the result of the Shapiro-Wilk normality test. If a significant result was found, the corresponding statistics are reported, and post-hoc multiple comparisons are performed.

The online version of this article includes the following source data for table 2:

**Source data 1.** Electrophysiological and molecular properties of pyramidal cells used for **Table 2**.

TTX did not impair the evoked vasoconstriction (Figure 5—figure supplement 2). These observations suggest that PGE2 and its EP1 and EP3 receptors mediate a sustained neurogenic vasoconstriction and that once PGE2 is released, its constrictive effect is independent of AP firing.

Astrocytes through 20-HETE and NPY interneurons are indirect intermediates of pyramidal cell-induced vasoconstriction

In addition to smooth muscle cells, PGE2 released by pyramidal cells can also activate astrocytes and neurons (Clasadonte et al., 2011; Di Cesare et al., 2006), which also express its receptors (Tasic et al., 2016; Zeisel et al., 2015). To assess whether astrocytes could mediate the PGE2-dependent vasoconstriction, we first targeted the large conductance Ca²⁺-activated (BK) channels and the 20-HETE pathways, both of which mediate astrocyte-derived vasoconstriction dependent on glutamatergic transmission (Girouard et al., 2010; Mulligan and MacVicar, 2004). Blockade of BK channels with paxilline (1 µM, Table 3, n=10 arterioles) did not impair the vascular response (Figure 6; AUC = -3.7 ± 0.3 x 10³ %s, t₁₈ = 0.03, p=1). Selective inhibition of the 20-HETE synthesizing enzyme, CYP450 ω-hydroxylase, with HET-0016 (100 nM, n=10 arterioles, Table 3) reduced the magnitude of the evoked vasoconstriction (Figure 6; AUC = -1.6 ± 0.7 x 10³ %s, t₁₈ = 3.32, p=0.0160). These data suggest that astrocytes partially mediate the vasoconstriction induced by pyramidal cells via 20-HETE

<table><caption>**Table 3.** IC₅₀ and concentrations of inhibitors, blocker and antagonists used in tissue.</caption><thead><tr><th rowspan="2">Inhibitor/antagonist</th><th colspan="3">Concentrations used for inhibition/antagonism</th></tr><tr><th>In vitro IC<sub>50</sub></th><th>Preparation</th><th>Concentration</th></tr></thead><tbody><tr><td>Indomethacin</td><td>COX-1: 22 nM; <i>Lora et al., 1998</i><br>COX-2: 87 nM</td><td>Mouse brain slices; <i>Lacroix et al., 2015</i></td><td>5 µM</td></tr><tr><td>SC560</td><td>COX-1: 9 nM; <i>Smith et al., 1998</i><br>COX-2: 6.3 µM</td><td>Mouse brain slices; <i>Lacroix et al., 2015</i></td><td>100 nM</td></tr><tr><td>NS-398</td><td>COX-1: 50 µM; <i>Lora et al., 1998</i><br>COX-2: 0.6 µM</td><td>Mouse brain slices; <i>Lacroix et al., 2015</i></td><td>10 µM</td></tr><tr><td>ONO-8130</td><td>EP1 receptors: 9.3 nM; <i>Säfholm et al., 2013a</i></td><td>isolated guinea pig trachea; <i>Säfholm et al., 2013b</i></td><td>10 nM</td></tr><tr><td>L798,106</td><td>EP3 receptors: 0.3 nM (Ki); <i>Juteau et al., 2001</i></td><td>Isolated mouse mesenteric arteries; <i>Chia et al., 2011</i></td><td>1 µM</td></tr><tr><td>AL8810</td><td>FP receptors: 426 nM (Ki); <i>Griffin et al., 1999</i></td><td>Isolated porcine retinal arterioles; <i>Oversø Hansen et al., 2015</i></td><td>10 µM</td></tr><tr><td>Paxilline</td><td>BK channels: 97 nM; <i>Tammaro et al., 2004</i></td><td>Mouse brain slices; <i>Girouard et al., 2010</i></td><td>1 µM</td></tr><tr><td>HET-0016</td><td>CYP4A isoforms: 35 nM; <i>Miyata et al., 2001</i></td><td>Mouse brain slices; <i>Blanco et al., 2008</i></td><td>100 nM</td></tr><tr><td>BIBP3226</td><td>Y1 receptors: 26 nM; <i>Rudolf et al., 1994</i></td><td>Mouse brain slices; <i>Sun et al., 2003</i></td><td>1 µM</td></tr></tbody></table>
<figure><img src="image_7.png" /><figcaption>**Figure 5.** GE2 mostly derived from COX-2 activity and its EP1 and EP3 receptors mediates vasoconstriction induced by optogenetically activated pyramidal cells. (A, B) Ex vivo effects of the COX1/2 inhibitor indomethacin (magenta, n=10 arterioles from 9 mice), the COX-1 inhibitor SC-560 (purple, n=10 arterioles from 7 mice), and the COX-2 inhibitor NS-398 (red, n=7 arterioles from 6 mice) on kinetics (A) and AUC (B) of arteriolar vasoconstriction induced by 20 Hz photostimulation (vertical cyan bar). In vivo experiments are highlighted by a black frame. (C) Optogenetic stimulation was achieved in vivo with an optic fiber through a chronic cranial window over the barrel cortex. (D) Left, diameter of pial arterioles labeled with fluorescein dextran (i.v.) was measured with line-scan crossing the vessel (white line). Right, Representative examples of vascular response upon photostimulation (10 Hz, 10 s) under control (top) and indomethacin condition (bottom). (E) Diameter changes upon photostimulation under control (black; n=5 arterioles, 4 mice) or indomethacin (magenta; n=4 arterioles, 4 mice) conditions. (F) Area under the curve of the diameter change in control (black) or indomethacin (magenta) conditions calculated between 20 and 40 s (unpaired, two-tailed Mann Whitney test, * p<0.05). (G, H) Effects of the EP1, EP3 and FP antagonists, ONO-8130 (10 nM, dark green, n=9 arterioles from 7 mice), L798,106 (1 µM, light green, n=9 arterioles from 5 mice) and AL8810 (10 µM, dark blue, n=9 arterioles from 7 mice), respectively, on kinetics (G) and AUC (H) of arteriolar vasoconstriction induced by 20 Hz photostimulation. The data are shown as the individual values and mean ± SEM. Dashed line represents the baseline. The SEMs envelope the mean traces. The shaded traces in A and G correspond to the control condition (from **Figure 1C** – 20 Hz). *, ** and *** statistically different from 20 Hz control condition with p<0.05, 0.01 and 0.001, respectively.</figcaption></figure>

The online version of this article includes the following source data and figure supplement(s) for figure 5:

**Source data 1.** Diameter measurements (µm) of individual arterioles used to determine diameter changes in **Figure 5A, B, G and H**.

**Source data 2.** Diameter changes of individual arterioles shown in **Figure 5E** under control condition.

**Source data 3.** Diameter changes of individual arterioles shown in **Figure 5E** after indomethacin treatment.

**Figure supplement 1.** Tonic PGE2 does not affect basal vascular tone.

**Figure supplement 1—source data 1.** Diameter measurements (µm) of individual arterioles used to determine diameter changes in **Figure 5—figure supplement 1**.

**Figure supplement 2.** PGE2 dose-dependently induces vasoconstriction.

**Figure supplement 2—source data 1.** Diameter measurements (µm) of individual arterioles used to determine diameter changes in **Figure 5—figure supplement 2**.
<figure><img src="image_8.png" /><figcaption>**Figure 6.** NPY Y1 receptors activation and 20-HETE synthesis mediates the vasoconstriction induced by pyramidal neurons. Effects of paxilline (1 µM, orange, n=10 arterioles from 6 mice), HET-0016 (100 nM, blue-grey, n=10 arterioles from 7 mice) and BIBP3226 (1 µM, yellow, n=10 arterioles from 6 mice) on (A) kinetics and (B) AUC of arteriolar vasoconstriction induced by 20 Hz photostimulation (vertical blue bar). Dashed line represents the baseline. The SEMs envelope the mean traces. The shaded traces correspond to the control condition (Figure 1C – 20 Hz). The data are shown as the individual values and mean ± SEM. * and *** statistically different from 20 Hz control condition with p<0.05 and 0.001.</figcaption></figure>

The online version of this article includes the following source data and figure supplement(s) for figure 6:

**Source data 1.** Diameter measurements (µm) of individual arterioles used to determine diameter changes in *Figure 6*.

**Figure supplement 1.** Vasoconstrictive pathways do not influence resting vascular tone.

**Figure supplement 1—source data 1.** Diameter measurements (µm) of individual arterioles used to determine diameter changes in **Figure 6—figure supplement 1**.

but not via K⁺ release. We next determined whether NPY, a potent vasoconstrictor (Cauli et al., 2004), was involved in neurogenic vasoconstriction. Antagonism of the NPY Y1 receptors by BIBP3226 (1 µM, Table 3, n=10 arterioles) abolished neurogenic vasoconstriction (Figure 6; AUC = -0.3 ± 0.2 x 10³ %·s, t₁₈ = 5.28, p=1.9512 x 10⁻⁵). These results suggest that neurogenic vasoconstriction induced by pyramidal cell photostimulation involves NPY release and the activation of Y1 receptors (Cauli et al., 2004; Karagiannis et al., 2009; Uhlirova et al., 2016) and astrocytes via 20-HETE in a glutamatergic-dependent and -independent manner.

Discussion

This study establishes that pyramidal cell activity leads to arteriolar vasoconstriction, and that the magnitude of the vasoconstriction depends on AP firing frequency and correlates with a graded increase in pyramidal cell somatic Ca²⁺. This vascular response partially involves glutamatergic transmission through direct and indirect mechanisms on arteriolar smooth muscle cells. Ex vivo and in vivo observations revealed that PGE2, predominantly produced by layer II-III COX-2 pyramidal cells, and its EP1 and EP3 receptors play a crucial role in neurogenic vasoconstriction. Pharmacological evidence indicates that some interneurons, via NPY release and activation of Y1 receptors, and to a lesser extent, astrocytes through 20-HETE and possibly COX-1 derived PGE2 play an intermediary role in this process (Figure 7).

We found that increasing the frequency of photostimulation in an ex vivo preparation caused nearby arteriole to go from a barely discernible response to robust vasoconstriction. In contrast, in vivo observations in anesthetized animals, with slower NVC compared to awake animals (Rungta et al., 2021; Uhlirova et al., 2016), have shown that the optogenetic stimulation of pyramidal cells results in a biphasic response: a fast hyperemic/vasodilatory response (Kahn et al., 2013; Lacroix et al., 2015; Scott and Murphy, 2012), which can be followed by a pronounced vasoconstriction (Figure 5; Uhlirova et al., 2016). The slow kinetics of the vascular response observed ex vivo is comparable with previous observations in slices (Cauli et al., 2004; Rancillac et al., 2006), and is likely due to the lower recording temperature compared to in vivo, which slows the synthesis of vasoactive mediators (Rancillac et al., 2006) and downstream reactions. The difficulty in observing vasodilation in cortical slices may be due to relaxed arterioles which favor vasoconstriction (Blanco et al., 2008). The evidence that neurogenic vasoconstriction is frequency-dependent (Figure 1) and that pharmacologically induced vasoconstriction persists in preconstructed (Girouard et al., 2010) or pressurized
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Possible pathways of vasoconstriction induced by pyramidal neurons. 20 Hz photostimulation induces activation of pyramidal neurons expressing channelrhodopsin-2 (ChR2<sub>H134R</sub>) and increases intracellular calcium (Ca<sup>2+</sup>). Arachidonic acid (AA) is released from membrane phospholipids (MPL) by phospholipases (PL) activated by intracellular Ca<sup>2+</sup> and is metabolized by type-1 and type-2 cyclooxygenases (COX-1 and COX-2) and prostaglandin E2 synthases (PGES) to produce prostaglandin E2 (PGE2). Three non-exclusive pathways can be proposed for arteriolar vasoconstriction in layer I: (1) PGE2 released into the extracellular space may act directly on arteriolar EP1 and EP3 receptors to induce smooth muscle cell constriction. (2) Glutamate released from pyramidal cells may activate neuropeptide Y (NPY) interneurons and NPY is released to act on vascular and neuronal Y1 receptors to constrict smooth muscle cells and promote glutamate release, respectively. Glutamate can also activate astrocytes to induce constriction through the 20-HETE and the COX-1/PGE2 pathways. (3) PGE2 may act on pre- and postsynaptic EP2 receptors to facilitate glutamate release and NPY interneuron activation, respectively.</figcaption></figure>

arterioles (*Dabertrand et al., 2013*), suggests that the neurogenic vasoconstriction primarily depends on a high pyramidal cell firing rate rather than on vascular tone.

Most previous observations did not report a decreased in CBF induced by optogenetic stimulation of pyramidal cells in vivo (*Lacroix et al., 2015; Scott and Murphy, 2012*). This may be attributed to differences in the photostimulation paradigm and/or the specific subtype of pyramidal cells that were stimulated. In our study, we used 10 s of photostimulation both in vivo and ex vivo. Earlier studies have employed shorter photostimulation times, lasting no more than 1 s (*Lacroix et al., 2015; Scott and Murphy, 2012; Uhlirova et al., 2016*), which may have resulted in an insufficient number of elicited APs to induce robust vasoconstriction. Furthermore, we observed that neurogenic vasoconstriction is highly dependent on COX-2, which is primarily expressed in layer II-III pyramidal cells (*Lacroix et al., 2015; Tasic et al., 2016; Zeisel et al., 2015*). In our study, photostimulation of almost all pyramidal cells in Emx1-Cre;Ai32 mice (*Gorski et al., 2002; Madisen et al., 2012*) likely resulted in the release of more COX-2 metabolites. Thy1-ChR2 mice used in previous studies (*Scott and Murphy, 2012*;
Uhlirova et al., 2016), on the other hand, express primarily ChR2 in layer V pyramidal cells (Kahn et al., 2013) which more rarely express COX-2.

Our ex vivo and in vivo observations revealed that PGE2, primarily derived from COX-2, plays a critical role in neurogenic vasoconstriction by activating EP1 and EP3 receptors expressed by vascular smooth muscle cells (Zhang et al., 2024). Previous studies have shown that COX-2 pyramidal cells, when activated in vivo by sensory stimulation or ex vivo, induce a NMDA-dependent increase in CBF and vasodilation through PGE2 and EP2/EP4 receptors (Lacroix et al., 2015; Lecrux et al., 2011; Niwa et al., 2000). Differences in the levels and/or sites of action of released PGE2 may explain the absence of secondary vasoconstriction. In Emx1-Cre;Ai32 mice (Gorski et al., 2002; Madisen et al., 2012), optogenetic stimulation may have activated a greater number of COX-2 pyramidal cells and resulted in a higher local release of PGE2 compared to sensory stimulation. Furthermore, since PGE2 is barely catabolized in the cerebral cortex (Alix et al., 2008), most of its removal occurs across the blood-brain barrier by specific transporters. The lack of blood perfusion in brain slices may impair this clearance mechanism, leading to PGE2 accumulation. It is noteworthy that the PGE2-induced vasoconstriction persisted after its removal (Figure 5—figure supplement 2). A high level of PGE2 may have facilitated the activation of the EP1 receptor, which has a lower affinity than the EP2/EP4 receptors (Boie et al., 1997). Additionally, it may have promoted the rapid desensitization of the dilatory EP4 receptor (Desai et al., 2000) thereby favoring vasoconstriction. Furthermore, PGE2 can induce either EP1-dependent arteriolar dilation or constriction depending on whether it is locally applied to capillaries or arterioles. Constriction prevails when both segments are exposed (Rosehart et al., 2021). Our photostimulation focused on superficial penetrating arterioles, which lack a capillary network in their close vicinity (Kasischke et al., 2011). This may have facilitated the direct EP1-mediated arteriolar constriction (Dabertrand et al., 2013; Rosehart et al., 2021). Overall, these observations suggest that COX-2 pyramidal cells can sequentially promote both vasodilation and vasoconstriction through the release of PGE2, depending on the context.

Consistent with previous reports in rodents (Lacroix et al., 2015; Tasic et al., 2016; Yamagata et al., 1993; Zeisel et al., 2015), the transcripts of the rate-limiting enzymes COX-1 and COX-2, were detected in subpopulations of mouse layer II-III pyramidal cells, respectively. COX-1/2 expression was found to be systematically associated with at least one PGE2 synthesizing enzyme. The major isoforms were cPGES and mPGES2, with the latter being less prevalent (Lacroix et al., 2015; Tasic et al., 2016; Zeisel et al., 2015). The low detection rate of mPGES1, an isoform co-induced with COX-2 by various stimuli (Takemiya et al., 2007; Yamagata et al., 2001), reflects its low constitutive basal expression level. The presence of PM-PGFS, CBR1 and AKR1B3 in layer II-III pyramidal cells is consistent with single-cell RNAseq data (Tasic et al., 2016; Zeisel et al., 2015). The expression of a PGFS was systematically observed in COX-2 positive pyramidal cells and in a majority of COX-1 positive neurons, similar to PGES. These observations collectively indicate that subpopulations of layer II-III pyramidal cells express the mRNAs required for PGE2 and PGF2α synthesis derived from COX-1 or COX-2 activity. Our pharmacological observations did not reveal a contribution of PGF2α in neurogenic vasoconstriction, despite the potential ability of pyramidal cells to produce it. This is likely because PGF2α is only detectable in pyramidal neurons under conditions where COX-2 is over-expressed (Takei et al., 2012).

Pyramidal cells may have an indirect effect on vascular activity through the activation of intermediate cell types, in addition to the direct vascular effects of PGE2 and glutamate (Zhang et al., 2024). Consistent with previous observations, we found that glutamate transmission from pyramidal cells is involved to some extent (Uhlirova et al., 2016). Additionally, we found that the NPY Y1 receptor plays a key role in neurogenic vasoconstriction. It is likely that glutamatergic transmission contributed to NPY release, considering that NPY GABAergic interneurons express a wide range of ionotropic and metabotropic glutamate receptors (Tasic et al., 2016; Zeisel et al., 2015). Consistently, the Y1 receptor has been shown to be involved in vasoconstriction induced by sensory and optogenetic stimulation of GABAergic interneurons (Uhlirova et al., 2016). Activation of group I metabotropic receptors in perivascular astrocytes has been shown to promote vasoconstriction via BK channel-dependent K⁺ release (Girouard et al., 2010) or 20-HETE (Mulligan and MacVicar, 2004). However, the neurogenic vasoconstriction was not affected by the blockade of BK channels, which rules out this astrocytic pathway. In contrast, the inhibition of ω-hydroxylase partially reduced neurogenic vasoconstriction, suggesting the involvement of 20-HETE. Additionally, astrocytes may also have contributed
to vasoconstriction through the release of PGE2 derived from COX-1 (Attwell et al., 2016), as indicated by its mild impairment under SC-560.

The observation that both EP1 and Y1 antagonists abolished the vasoconstriction suggests that PGE2 and NPY may act in series and/or in a more complex manner involving their neuronal receptors. One possibility is that PGE2 activates NPY interneurons via the EP1 receptor. However, NPY interneurons barely express its transcript (Tasic et al., 2016; Zeisel et al., 2015) and PGE2 constricts arterioles independently of AP firing, suggesting a direct vascular effect of PGE2. Nevertheless, PGE2 may have facilitated NPY release via pre- and postsynaptic EP2-signaling which have been shown to facilitate glutamate release (Sang et al., 2005) and to induce neuronal firing (Clasadonte et al., 2011), respectively. On the other hand, in addition to smooth muscle cells, Y1 receptors are also enriched in pyramidal neurons (Smith et al., 2019), including COX-2 positive ones (Tasic et al., 2016), and this receptor has been shown to increase extracellular glutamate in the hippocampus (Meurs et al., 2012). By promoting glutamate and possibly PGE2 release, neuronal activation of the Y1 receptor by NPY may also have favored direct (i.e. PGE2) and indirect (i.e. 20-HETE) vasoconstrictive pathways. The combined activation of vascular and neuronal Y1 receptors may explain the complete blockage of optogenetically induced vasoconstriction by its antagonist BIBP3226. In ex vivo relaxed arterioles, where vasoconstriction is favored (Blanco et al., 2008), Gq or Gs signaling of EP1 or Y1 receptors, respectively, appears sufficient to induce vasoconstriction. In vivo, where blood flow both induces myogenic tone and allows PGE2 clearance, NPY and PGE2 could also synergistically promote vasoconstriction by decreasing and increasing cAMP and Ca2+ levels, respectively, in smooth muscle cells. PGE2 and NPY may also exert temporally distinct vasoconstrictor effects. Indeed, exogenous application of NPY induces a rapid and transient vasoconstriction that returns to baseline levels after removal (Cauli et al., 2004), whereas PGE2-induced vasoconstriction is slower and more persistent (Figure 5—figure supplement 1). The more transient effect of NPY likely reflects the presence of multiple NPY-degrading enzymes (Wagner et al., 2015) and/or the desensitization of the Y1 receptor (Tsurumaki et al., 2003; Gicquiaux et al., 2002; Tsurumaki et al., 2002) which is not the case for PGE2 (Alix et al., 2008) and its vasoconstrictor receptors.

In awake mice, synchronous pyramidal cell activity occurs in the absence of any stimulus during the so-called resting state, but it is observed at a much lower frequency than that which triggers vasoconstriction and is associated with increased blood volume (Ma et al., 2016). Therefore, neurogenic vasoconstriction described here is unlikely to occur under these conditions. Brief sensory stimulation increases pyramidal cell activity and largely causes vasodilation in both awake and anesthetized animals (Rungta et al., 2021). This hyperemic response can be followed by delayed vasoconstriction (Devor et al., 2007) and involves NPY/Y1 receptor signaling (Uhlirova et al., 2016), similar to the mechanisms reported here. It remains unclear whether PGE2 signaling is also involved in this secondary response. During prolonged sensory stimulation the evoked hyperemic response appears to be more restricted to the activated area at the end of the stimulation than at the beginning (Berwick et al., 2008). It is possible that neurogenic vasoconstriction contributes to the later spatial confinement of the vascular response. The time-locked photostimulation of virtually all pyramidal cells leading to vasoconstriction would have resulted in hypersynchrony, a phenomenon that can be observed during sleep/wake transitions (Asadi-Pooya and Sperling, 2019). A decrease in hemodynamics has been reported during the transition from rapid eye movement sleep to wakefulness (Gheres et al., 2023; Tsai et al., 2021), possibly involving neurogenic vasoconstriction. Hypersynchrony is also observed in pathological conditions such as epileptic seizures (Jiruska et al., 2013) and in early stages of Alzheimer's disease (Bezzina et al., 2015; Palop et al., 2007). Although vasoconstriction observed in epilepsy (Farrell et al., 2016) exhibits similarities to the neurogenic vasoconstriction described herein, there are notable differences between the two. Like neurogenic vasoconstriction, seizure-induced hypoperfusion is dependent on COX-2 (Farrell et al., 2016; Tran et al., 2020) and, to some extent on PGE2 (Farrell et al., 2016), likely through EP1 and/or EP3 receptors. However, epileptic seizures induce the overexpression of both COX-2 and mPGES1 (Takemiya et al., 2007; Yamagata et al., 1993) as well as the ectopic expression of NPY (Baraban, 2004). Similar transcriptional upregulations have also been reported in Alzheimer's disease (Bezzina et al., 2015; Chaudhry et al., 2008; Palop et al., 2007; Pasinetti and Aisen, 1998) Additionally, PGF2α synthesis by COX-2 pyramidal cells is also observed during seizures (Takei et al., 2012). Taken together, these observations suggest that the mechanisms governing
neurogenic vasoconstriction are exacerbated in pathological hypersynchrony and may represent potential therapeutic targets.

This neurogenic vasoconstriction, observed during strong pyramidal cell activity, may seem counterintuitive as it would lead to an undersupply of energy substrates despite a high energy demand. However, vasoconstriction has been reported contralateral to the main activated area, despite bilateral increases in neuronal activity and blood glucose uptake (Devor et al., 2008), suggesting that neurogenic vasoconstriction plays a physiological role. Through glutamate uptake by astrocytes, neuronal activity stimulates blood glucose uptake and lactate release (Pellerin and Magistretti, 1994; Voutsinos-Porche et al., 2003). In addition to its role as an oxidative energy substrate for cortical neurons, lactate is also a signaling molecule that enhances their spiking activity (Karagiannis et al., 2021). Therefore, uncontrolled lactate supply and metabolism could potentially lead to deleterious hyperactivity (Cauli et al., 2023; Sada et al., 2015). Thus, the purpose of neurogenic vasoconstriction may be to restrict energy delivery to prevent an overexcitation of the cortical network.

Here, using multidisciplinary approaches, we describe a new mechanism of vasoconstriction that depends on a high firing rate of pyramidal cells. This neurogenic vasoconstriction primarily involves the release of COX-2-derived PGE2 and activation of EP1 and EP3 receptors. It is mediated by direct effects on vascular smooth muscle cells but also by indirect mechanisms involving NPY release from GABAergic interneurons and astrocytes by 20-HETE synthesis. In contrast to previously described mechanisms of neurogenic vasoconstriction, that have been mostly associated with GABAergic interneurons and neuronal inhibition (Cauli et al., 2004; Devor et al., 2007; Krawchuk et al., 2020; Lee et al., 2021; Uhlirova et al., 2016), our data suggest the involvement of glutamatergic excitatory neurons and increased neuronal activity. This finding will help to update the interpretation of the functional brain imaging signals used to map network activity in health and disease (ladecola, 2017; Zhang and Raichle, 2010). This excitatory form of neurogenic vasoconstriction may also help to understand the etiopathogenesis of epilepsy (Farrell et al., 2016; Tran et al., 2020) and Alzheimer's disease (Palop and Mucke, 2010) in which increased cortical network activity and hypoperfusion often overlap.

# Materials and methods

## Animals

Homozygous Emx1-Cre mice Jackson Laboratory, stock #005628, B6.129S2-Emx1<sup>tm1(cre)Kij</sup>/J (Gorski et al., 2002) were crossed with homozygous Ai32 mice [Jackson Laboratory, stock #012569, B6;129S-Gt(ROSA)26Sor<sup>tm32(CAG-COP4*H134R/EFP)Hze</sup>/J (Madisen et al., 2012)] to obtain heterozygous Emx1<sup>cre/WT</sup>;Ai32<sup>ChR2/WT</sup> mice for optogenetic stimulations. C57BL/6RJ mice were used for PGE2 and sulprostone exogenous applications, control optogenetic experiments and single-cell RT-PCR. 16-21 postnatal day-old females and males were used for all ex vivo experiments. Female Emx1<sup>cre/WT</sup>;Ai32<sup>ChR2/WT</sup> mice, 3- to 5-month-old, were used for in vivo experiments.

All experimental procedures using animals were carried out in strict accordance with French regulations (Code Rural R214/87 to R214/130) and conformed to the ethical guidelines of the European Communities Council Directive of September 22, 2010 (2010/63/UE). Mice were fed ad libitum and housed in a 12 hr light/dark cycle. In vivo experiments were done in accordance with the Institut national de la santé et de la recherche médicale (Inserm) animal care and approved by the ethical committee Charles Darwin (Comité national de réflexion éthique sur l'expérimentation animale – n°5) (protocol number #27135 2020091012114621).

## Ex vivo slice preparation

Mice were deeply anesthetized by isoflurane (IsoVet, Piramal Healthcare UK or IsoFlo, Axience) evaporation in an induction box then euthanized by decapitation. The brain was quickly removed and placed in cold (~4 °C), oxygenated artificial cerebrospinal fluid (aCSF) containing (in mM): 125 NaCl, 2.5 KCl, 1.25 NaH2PO4, 2 CaCl2, 1 MgCl2, 26 NaHCO3, 10 glucose, 15 sucrose and 1 kynurenic acid (Sigma-Aldrich). 300-µm-thick coronal slices containing the barrel cortex were cut with a vibratome (VT1000s; Leica) and were allowed to recover at room temperature for at least 45 min with oxygenated aCSF (95% O2/5% CO2; Devienne et al., 2018). The slices were then transferred to a submerged recording
chamber and perfused continuously at room temperature (20–25°C) at a rate of 2 ml/min with oxygen-
ated aCSF lacking kynurenic acid.

Whole-cell recordings

Patch pipettes (5.5±0.2 MΩ) pulled from borosilicate glass were filled with 8 µl of RNase free internal solution containing (in mM): 144 K-gluconate, 3 MgCl2, 0.5 EGTA, 10 HEPES, pH 7.2 (285/295 mOsm). For electrophysiological recordings combined with calcium imaging, EGTA was replaced by 200 µM Rhod-2 (20777, Cayman chemicals). Whole-cell recordings were performed using a patch-clamp amplifier (Axopatch 200B, MDS). Data were filtered at 5–10 kHz and digitized at 50 kHz using an acquisition board (Digidata 1440, MDS) attached to a personal computer running pCLAMP 10.2 software package (MDS). Electrophysiological properties were determined in current-clamp mode (Karagiannis et al., 2009). Membrane potential values were corrected for theoretical liquid junction potential (−15.6 mV). Resting membrane potential of neurons was measured immediately after passing in whole-cell configuration. Only neurons with a resting membrane potential more hyperpolarized than −60 mV were analyzed further.

Optogenetic stimulation

Optogenetic stimulation was achieved through the objective using a 470 nm light emitting device (LED, CoolLED, Precise Excite) attached to the epifluorescence port of a BX51WI microscope (Olympus) and a set of multiband filters consisting of an excitation filter (HC 392/474/554/635, Semrock), a dichroic mirror (BS 409/493/573/652, Semrock), and an emission filter (HC 432/515/595/730, Semrock). Photo-stimulation consisted of a 10 s train of 5ms light pulses at an intensity of 38 mW/mm² and delivered at five different frequencies (1, 2, 5, 10, and 20 Hz).

Infrared imaging

Blood vessels and cells were observed in slices under infrared illumination with Doddt gradient contrast optics (IR-DGC, Luigs and Neumann) using a double-port upright microscope (BX51WI, Olympus) and a collimated light emitting device (LED; 780 nm; ThorLabs) as the transmitted light source, a 40 X (LUMPlanF/IR, 40 X/0.80 W, Olympus) or a 60 X (LUMPlan FL/IR 60 X/0.90 W, Olympus) objective and a digital camera (OrcaFlash 4.0, Hamamatsu) attached to the front port of the microscope. Penetrating arterioles in layer I were selected by IR-DGC videomicroscopy based on their well-defined luminal diameter (10–40 µm), their length remaining in the focal plane for at least 50 µm (Lacroix et al., 2015), and the thickness of their wall (4.1±0.1 µm, n=176 blood vessels). A resting period of at least 30 min (Zonta et al., 2003) was observed after slice transfer. After light-induced responses, arteriolar contractility was tested by the application of aCSF containing the thromboxane A2 agonist, U46619 (100 nM; Cauli et al., 2004) or K⁺ enriched solution (composition in mM: 77.5 NaCl, 50 KCl, 1.25 NaH2PO4, 2 CaCl2, 1 MgCl2, 26 NaHCO3, 10 glucose, 15 sucrose). Vessels that did not constrict with these applications were discarded. Only one arteriole was monitored per slice receiving a single optogenetic or pharmacological stimulation. IR-DGC images were acquired at 0.1 Hz for pharmacological applications and at 1 Hz for optogenetic experiments using Imaging Workbench 6.1 software(Indec Biosystems). The focal plane was continuously maintained on-line using IR-DGC images of cells as anatomical landmarks (Lacroix et al., 2015).

Calcium imaging

Visually and electrophysiologically identified layer II-III pyramidal cells were filled with the calcium-
sensitive dye Rhod-2 (200 µM, Cayman chemical, 20777) using patch pipettes. Optical recordings
were made at least 15 min after passing in whole-cell configuration to allow for somatic diffusion
of the dye. Rhod-2 was excited with a 585 nm LED (Cool LED, Precise Excite) at an intensity of
0.56 mW/mm² and the filter set used for optogenetic stimulation using the Imaging Workbench 6.1
software(Indec Byosystems). IR-DGC and fluorescence images were acquired by alternating epifluo-
rescence and transmitted light sources. IR-DGC and fluorescence were respectively sampled at 5 Hz
and 1 Hz during baseline and optogenetic stimulation, respectively, and at 1 Hz and 0.2 Hz after
photostimulation. During photostimulation, bleed-through occurred in the Rhod-2 channel due to the
fluorescence of the EYFP-ChR2 transgene (Madisen et al., 2012). Therefore, the Ca²⁺ response could
not be reliably analyzed during this period. To compensate for potential x-y drifts, all images were
registered off-line using the 'StackReg' plug-in (Thévenaz et al., 1998) of the ImageJ 1.53 software. To define somatic regions of interest (ROIs), the soma was manually delineated from IR-DGC images. Fluorescence intensity changes (ΔF/F₀) were expressed as the ratio (F-F₀)/F₀ where F is the mean fluorescence intensity in the ROI at a given time point, and F₀ is the mean fluorescence intensity in the same ROI during the 30 s control baseline.

## Drugs

All pharmacological compounds were bath applied after a 5-min baseline, and vascular dynamics were recorded during bath application. The following drugs were dissolved in water: D-(-)-2-amino-5-phosphonopentanoic acid (D-AP5, 50 µM, Hello Bio, HB0225), 6,7-dinitroquinoxaline-2,3-dione (DNQX, 10 µM, Hello Bio, HB0262), LY367385 (100 µM, Hello Bio, HB0398) and BIBP3226 (1 µM, Tocris, 2707). Tetrodotoxin (TTX, 1 µM, L8503, Latoxan) was dissolved in 90% acetic acid. PGE2 (HB3460, Hello Bio), sulprostone (10 µM, Cayman chemical, 14765), 2-methyl-6-(phenylethynyl)pyridine (MPEP, 50 µM, Hello Bio, HB0426) and 9,11-dideoxy-9α,11α-methanoepoxy prostaglandin F2α (U-46619, 100 nM, Enzo, BML-PG023) were dissolved in ethanol. Indomethacin (5 µM, Sigma-Aldrich, I7378), SC-560 (100 nM, Sigma-Aldrich, S2064-5MG), NS-398 (10 µM, Enzo, BML-EI261), ONO-8130 (10 nM, Tocris, 5406), L-798,106 (1 µM, Cayman chemical, 11129), AL8810 (10 µM, Cayman chemical, 16735), paxilline (10 µM, Tocris, 2006) and HET0016 (100 µM, Merck, SML2416-5MG) in DMSO. Acetic acid, ethanol and DMSO doses were always used below 0.1%. Synthesis inhibitors and BIBP3226 were applied at least 30 min before optogenetic stimulation. TTX was applied at least 15 min before, while glutamate receptor antagonists and paxilline were applied at least 10 and 5 min before, respectively.

## Vascular reactivity analysis

To compensate for potential x-y drifts, all images were realigned off-line using the 'StackReg' plug-in Thévenaz et al., 1998 of the ImageJ 1.53 software. Luminal diameter was measured in layer I on registered images using custom analysis software developed in MATLAB (MathWorks) (Lacroix et al., 2015). To avoid potential drawbacks due to vessel instability, only arterioles with a stable diameter were analyzed further. Arterioles were considered stable if the relative standard deviation of their diameter during the baseline period was less than 5% (Lacroix et al., 2015). Comparison of the mean arteriolar diameter during the 5 min baseline and the 5 min final pharmacological treatments revealed that all drugs, except AL8810, had no effect on resting diameter (Figure 3—figure supplement 1, Figure 5—figure supplement 1, Figure 6—figure supplement 1).

Diameter changes (ΔD/D₀) were expressed as (D_t - D₀)/D₀ where D_t is the diameter at the time t and D₀ is the mean diameter during the baseline period. To eliminate sharp artifacts due to transient loss of focus, diameter change traces were smoothed using a sliding three-point median filter. The overall vascular response over time was captured by the area under the curve of diameter changes after photostimulation. To determine the onset of vasoconstriction, a Z-score was calculated from the diameter change traces using the formula: Z = (x - μ)/ σ, where both the mean μ and the standard deviation σ were calculated from the values before photostimulation. Onset of vasoconstriction was defined as the time after the start of photostimulation at which the Z-score exceeded or fell below -a value of -1.96 (95% criteria) for 10 s. If a vessel showed no vasoconstriction, the onset was arbitrarily set at 1800s. Graphs were generated using R software version 4.3.0 (R Development Core Team, 2023) and Matplotlib package (Caswell, 2023).

## Intrinsic optical signals analysis

Variations in IR light transmittance (ΔT) (Zhou et al., 2010) were determined using ImageJ 1.53 software according to: ΔT = (T_t - T_0)/ T_0 where T_t is the light transmittance at a time t and T_0 is the average light transmittance during the baseline period of a squared region of interest of 100 µm x 100 µm manually delineated in layer I. The rate of ΔT change was determined as the first derivative of ΔT (dΔT/dt, where ΔT is the change in light transmittance and t is time). Slices that showed a maximum rate of increase of dΔT/dt greater than 2% /s, indicating the occurrence of spreading depression (Zhou et al., 2010), were excluded.
## Surgery

Chronic cranial windows were implanted one week after the head bar surgery as previously described (Tournissac et al., 2022). We used a 100 µm thick glass coverslip over the barrel cortex (~4 mm²). Before two-photon experiments, a recovery period of 7–10 days minimum was respected.

## In vivo two-photon imaging and photostimulation

For two-photon excitation, we used a femtosecond laser (Mai Tai eHP; SpectraPhysics) with a dispersion compensation module (Deepsee; SpectraPhysics) emitting 70-fs pulses at 80 MHz. The laser power was attenuated by an acousto-optical modulator (AA Optoelectronic, MT110-B50-A1.5-IR-Hk). Scanning was performed with Galvanometric scanner (GS) mirrors (8315KM60B; Cambridge Technology). Fluorescein was excited at 920 nm and the emitted light was collected with a LUMFLN60XW (Olympus, 1.1 NA) water immersion objective. Collected photons were sorted using a dichroic mirror centered at 570 nm, a FF01-525/50 nm filter (Semrock) and a GaAsP (Hamamatsu) photomultipliers tube. Customized LabView software was used to control the system. Line-scans were drawn across pial vessels to measure the change in arterioles diameter, which are not compromised by the fluorescence from the ChR2-EYFP transgene in the parenchyma (Madisen et al., 2012), are less affected by potential movement in the x-y plan than in penetrating arterioles, and whose dilation dynamics are similar in the somatosensory cortex (Rungta et al., 2021).

Mice were anesthetized with a mixture of ketamine and medetomidine (100 and 0.5 mg/kg, respectively, intraperitoneally (i.p.)) during imaging sessions. Body temperature was maintained at 36.5°C with a retro-controlled heating pad. Fluorescein dextran (70 kDa) was injected i.v. through a retro-orbital injection to label brain vessels. Mice received continuous air through a nose cone supplemented with oxygen to reach a final concentration of 30% O₂.

Photostimulation was delivered with a 473 nm laser (Coblot MLD, Sweden) through an optical fiber placed above the glass coverslip and directed at the pial artery of interest. Each photostimulation consisted of a 10-second train of 5 ms light pulses at an intensity of 1 mW delivered at 10 Hz, with a 5-minute interstimulus interval to allow full recovery to baseline. Indomethacin (10 mg/kg, #15425529, Thermo Fisher Scientific) was administered i.v. through a retro-orbital injection.

## Imaging analysis

Pial arteriole diameter change was determined with line-scan acquisitions and a home-made Matlab script as previously described (Rungta et al., 2018). Trials from the same vessel were averaged (with a 0.1 s interpolation) for analysis. Area under the curve and statistics were performed using GraphPad Prim (version 6).

## Cytoplasm harvesting and single-cell RT-PCR

At the end of the whole-cell recording, which lasted less than 15 min, the cytoplasmic content was collected in the recording pipette by applying a gentle negative pressure. The pipette's content was expelled into a test tube and RT was performed in a final volume of 10 µl as described previously (Devienne et al., 2018). The scRT-PCR protocol was designed to probe simultaneously the expression of prostaglandins synthesizing enzymes and neuronal markers (Lacroix et al., 2015). Prostaglandins synthesizing enzymes included COX-1 and COX-2, the terminal PGE2 synthases (PGES): mPGES1, mPGES2 and cPGES, the terminal PGF2αsynthases (PGFS): PM-PGFS (Prxl2b) and AKR1B3 and the carbonyl reductase CBR1. Neuronal markers included the vesicular glutamate transporter, vGluT1, and the two isoforms of glutamic acid decarboxylase, GAD65 and GAD67. Two-step amplification was performed essentially as described (Devienne et al., 2018). First, cDNAs present in the 10 µl reverse transcription reaction were simultaneously amplified with all external primer pairs listed in Appendix 1—key resources table. Taq polymerase (2.5 U; QIAGEN) and external primers mix (20 pmol each) were added to the manufacturer's buffer (final volume, 100 µl), and 20 cycles (95 °C, 30 s; 60 °C, 30 s; and 72 °C, 35 s) of PCR were performed. Second rounds of PCR were performed using 1 µl of the first PCR product as a template. In this second round, each cDNA was amplified individually using its specific nested primer pair (Appendix 1—key resources table) by performing 35 PCR cycles (as described above). 10 µl of each individual PCR product were run on a 2% agarose gel stained with ethidium bromide using ΦX174 digested by HaeIII as a molecular weight marker. The efficiency of the protocol was validated using 500 pg of total forebrain RNAs (Figure 4—figure supplement 1).
## Statistical analyses

Statistical analyses were performed using GraphPad Prism version 7.00 for Windows (GraphPad Software, La Jolla California USA, https://www.graphpad.com/) and R software version 4.3.0 (R Development Core Team, 2023). Normality of distribution was assessed using the Shapiro-Wilk tests. Equality of variance was assessed using Brown-Forsythe tests for comparisons between groups and using F-tests for comparisons with a control group. Parametric tests were only used if these criteria were met. Statistical significance of morphological and physiological properties of penetrating arterioles was determined using one-way ANOVA for comparison between groups. Statistical significance of calcium was determined using two-tailed unpaired t-tests and Statistical significance of vascular responses were appreciated using Tukey posthoc tests for the different frequencies conditions and using Dunnett's posthoc tests for the different pharmacological conditions compared to the 20 Hz condition without pharmacological compound. False discovery rate correction was used for multiple comparisons. Statistical significance of vascular diameter for drug applications was determined using two-tailed paired t-tests. Statistical significance on all figures uses the following convention: \*p<0.05, \*\*p<0.01 and \*\*\*p<0.001.

## Acknowledgements

The authors thank Dr Rebecca Piskorowski for constructive criticism of the manuscript. We acknowledge the invaluable support of the animal facilities of IBPS (RongIBPS) and Institut de la vision for their expert care and maintenance of the animals used in this study. Financial support was provided by grants from the Agence Nationale pour la Recherche (ANR-17-CE37-0010-03, BC; CE37_2020_TF-fUS-CADASIL, SC; ANR-20-CE14-0025, DL; ANR-23-CE14-0038-01, BC), the Fondation Alzheimer France (M21JRCN009, SC) and the i-Bio initiative of Sorbonne University (BC). BLG and EB were supported by fellowships from Fondation pour la Recherche sur Alzheimer and MT by a fellowship from the Fondation pour la Recherche Médicale (SPF201909009103)

## Additional information

### Funding

<table><thead><tr><td>Funder</td><td>Grant reference number</td><td>Author</td></tr></thead><tbody><tr><td>Fondation pour la Recherche sur Alzheimer</td><td></td><td>Benjamin Le Gac<br>Esther Belzic</td></tr><tr><td>Agence Nationale de la Recherche</td><td>ANR-17-CE37-0010-03</td><td>Bruno Cauli</td></tr><tr><td>Agence Nationale de la Recherche</td><td>CE37_2020_TF-fUS-CADASIL</td><td>Serge Charpak</td></tr><tr><td>Agence Nationale de la Recherche</td><td>ANR-20-CE14-0025</td><td>Dongdong Li</td></tr><tr><td>Agence Nationale de la Recherche</td><td>ANR-23-CE14-0038-01</td><td>Bruno Cauli</td></tr><tr><td>FONDATION ALZHEIMER</td><td>M21JRCN009</td><td>Serge Charpak</td></tr><tr><td>i-Bio Initiative</td><td></td><td>Bruno Cauli</td></tr><tr><td>Fondation pour la Recherche Médicale</td><td>SPF201909009103</td><td>Marine Tournissac</td></tr></tbody></table>

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

### Author contributions

Benjamin Le Gac, Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing; Marine Tournissac, Conceptualization, Data curation, Formal analysis, Investigation, Visualization, Methodology,
Writing – review and editing; Esther Belzic, Data curation, Formal analysis, Funding acquisition, Investi-
gation, Writing – review and editing; Sandrine Picaud, Formal analysis, Investigation; Isabelle Dusart,
Writing – review and editing; Hédi Soula, Data curation; Dongdong Li, Conceptualization, Data cura-
tion, Funding acquisition, Investigation, Methodology, Writing – review and editing; Serge Charpak,
Conceptualization, Resources, Data curation, Formal analysis, Supervision, Funding acquisition, Investi-
gation, Methodology, Writing – review and editing; Bruno Cauli, Conceptualization, Resources, Data
curation, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization,
Methodology, Writing – original draft, Project administration, Writing – review and editing

Author ORCIDs

Benjamin Le Gac https://orcid.org/0000-0002-4703-6217

Marine Tournissac https://orcid.org/0000-0002-1990-5075

Isabelle Dusart https://orcid.org/0000-0002-3211-4323

Hédi Soula https://orcid.org/0000-0001-5306-9712

Dongdong Li https://orcid.org/0000-0002-6731-4771

Serge Charpak https://orcid.org/0000-0002-5516-1245

Bruno Cauli https://orcid.org/0000-0003-1471-4621

Ethics

All experimental procedures using animals were carried out in strict accordance with French regu-
lations (Code Rural R214/87 to R214/130) and conformed to the ethical guidelines of the European
Communities Council Directive of September 22, 2010 (2010/63/UE). Mice were fed ad libitum and
housed in a 12-hour light/dark cycle. In vivo experiments were done in accordance with the Institut
national de la santé et de la recherche médicale (Inserm) animal care and approved by the ethical
committee Charles Darwin (Comité national de réflexion éthique sur l'expérimentation animale - n°5)
protocoll number #27135 2020091012114621).

Peer review material

Reviewer #1 (Public review): https://doi.org/10.7554/eLife.102424.3.sa1

Reviewer #2 (Public review): https://doi.org/10.7554/eLife.102424.3.sa2

Author response https://doi.org/10.7554/eLife.102424.3.sa3

Additional files

Supplementary files

MDAR checklist

Source code 1. Matlab script for blood vessel analysis. The method has been described in *Lacroix et al., 2015*.

Data availability

All data generated or analyzed during this study are included in the manuscript and supporting files.
Source data files have been provided for Figures 1 to 6. Source code for blood vessel analysis is
provided in Source code 1.

References

Alix E, Schmitt C, Strazielle N, Ghersi-Egea JF. 2008. Prostaglandin E2 metabolism in rat brain: role of the blood-brain interfaces. *Cerebrospinal Fluid Research* **5**:5. DOI: https://doi.org/10.1186/1743-8454-5-5, PMID: 18318891

Asadi-Pooya AA, Sperling MR. 2019. Normal awake, drowsy, and sleep EEG patterns that might be overinterpreted as abnormal. *Journal of Clinical Neurophysiology* **36**:250–256. DOI: https://doi.org/10.1097/WNP.0000000000000585, PMID: 31274687

Attwell D, Buchan AM, Charpak S, Lauritzen M, Macvicar BA, Newman EA. 2010. Glial and neuronal control of brain blood flow. *Nature* **468**:232–243. DOI: https://doi.org/10.1038/nature09613, PMID: 21068832

Attwell D, Mishra A, Hall CN, O'Farrell FM, Dalkara T. 2016. What is a pericyte? *Journal of Cerebral Blood Flow and Metabolism* **36**:451–455. DOI: https://doi.org/10.1177/0271678X15610340, PMID: 26661200

Baraban SC. 2004. Neuropeptide Y and epilepsy: recent progress, prospects and controversies. *Neuropeptides* **38**:261–265. DOI: https://doi.org/10.1016/j.npep.2004.04.006, PMID: 15337378

Berwick J, Johnston D, Jones M, Martindale J, Martin C, Kennerley AJ, Redgrave P, Mayhew JEW. 2008. Fine detail of neurovascular coupling revealed by spatiotemporal analysis of the hemodynamic response to single
whisker stimulation in rat barrel cortex. *Journal of Neurophysiology* **99**:787-798. DOI: https://doi.org/10.1152/jn.00658.2007, PMID: 18046008

Bezzina C, Verret L, Juan C, Remaud J, Halley H, Rampon C, Dahan L. 2015. Early onset of hypersynchronous network activity and expression of a marker of chronic seizures in the Tg2576 mouse model of Alzheimer's disease. *PLOS ONE* **10**:e0119910. DOI: https://doi.org/10.1371/journal.pone.0119910, PMID: 25768013

Blanco VM, Stern JE, Filosa JA. 2008. Tone-dependent vascular responses to astrocyte-derived signals. *American Journal of Physiology-Heart and Circulatory Physiology* **294**:H2855-H2863. DOI: https://doi.org/10.1152/ajpheart.91451.2007

Boie Y, Stocco R, Sawyer N, Slipetz DM, Ungrin MD, Neuschäfer-Rube F, Püschel GP, Metters KM, Abramovitz M. 1997. Molecular cloning and characterization of the four rat prostaglandin E2 prostanoid receptor subtypes. *European Journal of Pharmacology* **340**:227-241. DOI: https://doi.org/10.1016/s0014-2999(97)01383-6, PMID: 9537820

Cabezas C, Irinopoulou T, Cauli B, Poncer JC. 2013. Molecular and functional characterization of GAD67-expressing, newborn granule cells in mouse dentate gyrus. *Frontiers in Neural Circuits* **7**:60. DOI: https://doi.org/10.3389/fncir.2013.00060, PMID: 23565079

Caswell TA. 2023. Matplotlib/matplotlib: REL: v3.7.2. v3.7.2. Zenodo. DOI: https://doi.org/10.5281/zenodo.8118151

Cauli B, Tong XK, Rancillac A, Serluca N, Lambolez B, Rossier J, Hamel E. 2004. Cortical GABA interneurons in neurovascular coupling: relays for subcortical vasoactive pathways. *The Journal of Neuroscience* **24**:8940-8949. DOI: https://doi.org/10.1523/JNEUROSCI.3065-04.2004, PMID: 15483113

Cauli B, Hamel E. 2010. Revisiting the role of neurons in neurovascular coupling. *Frontiers in Neuroenergetics* **2**:9. DOI: https://doi.org/10.3389/fnene.2010.00009, PMID: 20616884

Cauli B, Dusart I, Li D. 2023. Lactate as a determinant of neuronal excitability, neuroenergetics and beyond. *Neurobiology of Disease* **184**:106207. DOI: https://doi.org/10.1016/j.nbd.2023.106207, PMID: 37331530

Chaudhry UA, Zhuang H, Crain BJ, Doré S. 2008. Elevated microsomal prostaglandin-E synthase-1 in Alzheimer's disease. *Alzheimer's & Dementia* **4**:6-13. DOI: https://doi.org/10.1016/j.jalz.2007.10.015, PMID: 18631945

Chia E, Kagota S, Wijekoon EP, McGuire JJ. 2011. Protection of protease-activated receptor 2 mediated vasodilatation against angiotensin II-induced vascular dysfunction in mice. *BMC Pharmacology* **11**:10. DOI: https://doi.org/10.1186/1471-2210-11-10, PMID: 21955547

Chung DY, Sadeghian H, Qin T, Lule S, Lee H, Karakaya F, Goins S, Oka F, Yaseen MA, Houben T, Tolner EA, van den Maagdenberg AMJM, Whalen MJ, Sakadžic S, Ayata C. 2019. Determinants of optogenetic cortical spreading depolarizations. *Cerebral Cortex* **29**:1150-1161. DOI: https://doi.org/10.1093/cercor/bhy021, PMID: 29425263

Clasadonte J, Poulain P, Hanchate NK, Corfas G, Ojeda SR, Prevot V. 2011. Prostaglandin E₂ release from astrocytes triggers gonadotropin-releasing hormone (GnRH) neuron firing via EP2 receptor activation. *PNAS* **108**:16104-16109. DOI: https://doi.org/10.1073/pnas.1107533108

Dabertrand F, Hannah RM, Pearson JM, Hill-Eubanks DC, Brayden JE, Nelson MT. 2013. Prostaglandin E2, a postulated astrocyte-derived neurovascular coupling agent, constricts rather than dilates parenchymal arterioles. *Journal of Cerebral Blood Flow and Metabolism* **33**:479-482. DOI: https://doi.org/10.1038/jcbfm.2013.9, PMID: 23385200

Desai S, April H, Nwaneshiu du C, Ashby B. 2000. Comparison of agonist-induced internalization of the human EP2 and EP4 prostaglandin receptors: role of the carboxyl terminus in EP4 receptor sequestration. *Molecular Pharmacology* **58**:1279-1286. DOI: https://doi.org/10.1124/mol.58.6.1279, PMID: 11093764

Devienne G, Le Gac B, Piquet J, Cauli B. 2018. Single cell multiplex reverse transcription polymerase chain reaction after patch-clamp. *Journal of Visualized Experiments* **1**:57627. DOI: https://doi.org/10.3791/57627, PMID: 29985318

Devor A, Tian P, Nishimura N, Teng IC, Hillman EMC, Narayanan SN, Ulbert I, Boas DA, Kleinfeld D, Dale AM. 2007. Suppressed neuronal activity and concurrent arteriolar vasoconstriction may explain negative blood oxygenation level-dependent signal. *The Journal of Neuroscience* **27**:4452-4459. DOI: https://doi.org/10.1523/JNEUROSCI.0134-07.2007, PMID: 17442830

Devor A, Hillman EMC, Tian P, Waeber C, Teng IC, Ruvinskaya L, Shalinsky MH, Zhu H, Haslinger RH, Narayanan SN, Ulbert I, Dunn AK, Lo EH, Rosen BR, Dale AM, Kleinfeld D, Boas DA. 2008. Stimulus-induced changes in blood flow and 2-deoxyglucose uptake dissociate in ipsilateral somatosensory cortex. *The Journal of Neuroscience* **28**:14347-14357. DOI: https://doi.org/10.1523/JNEUROSCI.4307-08.2008, PMID: 19118167

Di Cesare A, Del Piccolo P, Zacchetti D, Grohovaz F. 2006. EP2 receptor stimulation promotes calcium responses in astrocytes via activation of the adenylyl cyclase pathway. *Cellular and Molecular Life Sciences* **63**:2546-2553. DOI: https://doi.org/10.1007/s00018-006-6262-9, PMID: 17058034

Farrell JS, Gaxiola-Valdez I, Wolff MD, David LS, Dika HI, Geeraert BL, Wang R, Singh S, Spanswick SC, Dunn JF, Antle MC, Federico P, Campbell Teskey G. 2016. Postictal behavioural impairments are due to a severe prolonged hypoperfusion/ hypoxia event that is COX-2 dependent. *eLife* **1**:e2001. DOI: https://doi.org/10.7554/eLife.19352.001

Férézou I, Cauli B, Hill EL, Rossier J, Hamel E, Lambolez B. 2002. 5-HT3 receptors mediate serotonergic fast synaptic excitation of neocortical vasoactive intestinal peptide/cholecystokinin interneurons. *The Journal of Neuroscience* **22**:7389-7397. DOI: https://doi.org/10.1523/JNEUROSCI.22-17-07389.2002, PMID: 12196560

Gheres KW, Ünsal HS, Han X, Zhang Q, Turner KL, Zhang N, Drew PJ. 2023. Arousal state transitions occlude sensory-evoked neurovascular coupling in neonatal mice. *Communications Biology* **6**:738. DOI: https://doi.org/10.1038/s42003-023-05121-5, PMID: 37460780
Gicquiaux H, Lecat S, Gaire M, Dieterlen A, Mély Y, Takeda K, Bucher B, Galzi JL. 2002. Rapid internalization and recycling of the human neuropeptide Y (1) receptor. *The Journal of Biological Chemistry* **277**:6645–6655. DOI: https://doi.org/10.1074/jbc.M107224200, PMID: 11741903

Girouard H, Bonev AD, Hannah RM, Meredith A, Aldrich RW, Nelson MT. 2010. Astrocytic endfoot Ca2+ and BK channels determine both arteriolar dilation and constriction. *PNAS* **107**:3811–3816. DOI: https://doi.org/10.1073/pnas.0914722107, PMID: 20133576

Gordon GRJ, Choi HB, Rungta RL, Ellis-Davies GCR, MacVicar BA. 2008. Brain metabolism dictates the polarity of astrocyte control over arterioles. *Nature* **456**:745–749. DOI: https://doi.org/10.1038/nature07525, PMID: 18971930

Gorski JA, Talley T, Qiu M, Puelles L, Rubenstein JLR, Jones KR. 2002. Cortical excitatory neurons and glia, but not GABAergic neurons, are produced in the Emx1-expressing lineage. *The Journal of Neuroscience* **22**:6309–6314. DOI: https://doi.org/10.1523/JNEUROSCI.22-15-06309.2002, PMID: 12151506

Griffin BW, Klimko P, Crider JY, Sharif NA. 1999. AL-8810: a novel prostaglandin F2 alpha analog with selective antagonist effects at the prostaglandin F2 alpha (FP) receptor. *The Journal of Pharmacology and Experimental Therapeutics* **290**:1278–1284 PMID: 10454504.

Grutzendler J, Nedergaard M. 2019. Cellular control of brain capillary blood flow: in vivo imaging veritas. *Trends in Neurosciences* **42**:528–536. DOI: https://doi.org/10.1016/j.tins.2019.05.009, PMID: 31255380

Hartmann DA, Berthiaume A-A, Grant RI, Harrill SA, Koski T, Tieu T, McDowell KP, Faino AV, Kelly AL, Shih AY. 2021. Brain capillary pericytes exert a substantial but slow influence on blood flow. *Nature Neuroscience* **24**:633–645. DOI: https://doi.org/10.1038/s41593-020-00793-2, PMID: 33603231

Hill RA, Tong L, Yuan P, Murikinati S, Gupta S, Grutzendler J. 2015. Regional blood flow in the normal and ischemic brain is controlled by arteriolar smooth muscle cell contractility and not by capillary pericytes. *Neuron* **87**:95–110. DOI: https://doi.org/10.1016/j.neuron.2015.06.001, PMID: 26119027

Iadecola C, Nedergaard M. 2007. Glial regulation of the cerebral microvasculature. *Nature Neuroscience* **10**:1369–1376. DOI: https://doi.org/10.1038/nn2003, PMID: 17965657

Iadecola C. 2017. The neurovascular unit coming of age: a journey through neurovascular coupling in health and disease. *Neuron* **96**:17–42. DOI: https://doi.org/10.1016/j.neuron.2017.07.030, PMID: 28957666

Jiruska P, de Curtis M, Jefferys JGR, Schevon CA, Schiff SJ, Schindler K. 2013. Synchronization and desynchronization in epilepsy: controversies and hypotheses. *The Journal of Physiology* **591**:787–797. DOI: https://doi.org/10.1113/jphysiol.2012.239590

Juteau H, Gareau Y, Labelle M, Sturino CF, Sawyer N, Tremblay N, Lamontagne S, Carrière MC, Denis D, Metters KM. 2001. Structure–activity relationship of cinnamic acylsulfonamide analogues on the human EP3 prostanoid receptor. *Bioorganic & Medicinal Chemistry* **9**:1977–1984. DOI: https://doi.org/10.1016/S0968-0896(01)00110-9

Kahn I, Knoblich U, Desai M, Bernstein J, Graybiel AM, Boyden ES, Buckner RL, Moore CI. 2013. Optogenetic drive of neocortical pyramidal neurons generates fMRI signals that are correlated with spiking activity. *Brain Research* **1511**:33–45. DOI: https://doi.org/10.1016/j.brainres.2013.03.011, PMID: 23523914

Karagiannis A, Gallopin T, Dávid C, Battaglia D, Geoffroy H, Rossier J, Hillman EMC, Staiger JF, Cauli B. 2009. Classification of NPY-expressing neocortical interneurons. *The Journal of Neuroscience* **29**:3642–3659. DOI: https://doi.org/10.1523/JNEUROSCI.0058-09.2009, PMID: 19295167

Karagiannis A, Gallopin T, Lacroix A, Plaisier F, Piquet J, Geoffroy H, Hepp R, Naudé J, Le Gac B, Egger R, Lambolez B, Li D, Rossier J, Staiger JF, Imamura H, Seino S, Roeper J, Cauli B. 2021. Lactate is an energy substrate for rodent cortical neurons and enhances their firing activity. *eLife* **1**:e1424. DOI: https://doi.org/10.7554/eLife.71424, PMID: 34766906

Kasischke KA, Lambert EM, Panepento B, Sun A, Gelbard HA, Burgess RW, Foster TH, Nedergaard M. 2011. Two-photon NADH imaging exposes boundaries of oxygen diffusion in cortical vascular supply regions. *Journal of Cerebral Blood Flow and Metabolism* **31**:68–81. DOI: https://doi.org/10.1038/jcbfm.2010.158, PMID: 20859293

Krawchuk MB, Ruff CF, Yang X, Ross SE, Vazquez AL. 2020. Optogenetic assessment of VIP, PV, SOM and NOS inhibitory neuron activity and cerebral blood flow regulation in mouse somato-sensory cortex. *Journal of Cerebral Blood Flow and Metabolism* **40**:1427–1440. DOI: https://doi.org/10.1177/0271678X19870105, PMID: 31418628

Lacroix A, Toussay X, Anenberg E, Lecrux C, Ferreirós N, Karagiannis A, Plaisier F, Chausson P, Jarlier F, Burgess SA, Hillman E, Tegeder I, Murphy TH, Hamel E, Cauli B, Ferreiro N, Burgess SA, Hillman E, Tegeder I, Murphy TH, et al. 2015. COX-2-derived prostaglandin e2 produced by pyramidal neurons contributes to neurovascular coupling in the rodent cerebral cortex. *The Journal of Neuroscience* **35**:11791–11810. DOI: https://doi.org/10.1523/JNEUROSCI.0651-15.2015, PMID: 26311764

Lecrux C, Toussay X, Kocharyan A, Fernandes P, Neupane S, Lévesque M, Plaisier F, Shmuel A, Cauli B, Hamel E. 2011. Pyramidal neurons are "neurogenic hubs" in the neurovascular coupling response to whisker stimulation. *The Journal of Neuroscience* **31**:9836–9847. DOI: https://doi.org/10.1523/JNEUROSCI.4943-10.2011, PMID: 21734275

Lee JMJ, Stile CL, Bice AR, Rosenthal ZP, Yan P, Snyder AZ, Lee JMJ, Bauer AQ. 2021. Opposed hemodynamic responses following increased excitation and parvalbumin-based inhibition. *Journal of Cerebral Blood Flow & Metabolism* **41**:841–856. DOI: https://doi.org/10.1177/0271678X20930831

Lin JY, Lin MZ, Steinbach P, Tsien RY. 2009. Characterization of engineered channelrhodopsin variants with improved properties and kinetics. *Biophysical Journal* **96**:1803–1814. DOI: https://doi.org/10.1016/j.bpj.2008.11.034, PMID: 19254539
Lora M, Denault JB, Leduc R, de Brum-Fernandes AJ. 1998. Systematic pharmacological approach to the characterization of NSAIDs. *Prostaglandins, Leukotrienes, and Essential Fatty Acids* **59**:55-62. DOI: https://doi.org/10.1016/s0952-3278(98)90052-7, PMID: 9758208

Ma Y, Shaik MA, Kozberg MG, Kim SH, Portes JP, Timerman D, Hillman EMC. 2016. Resting-state hemodynamics are spatiotemporally coupled to synchronized and symmetric neural activity in excitatory neurons. *PNAS* **113**:E8463-E8471. DOI: https://doi.org/10.1073/pnas.1525369113, PMID: 27974609

Madisen L, Mao T, Koch H, Zhuo J, Berenyi A, Fujisawa S, Hsu Y-WA, Gu X, Zanella S, Kidney J, Gu H, Mao Y, Hooks BM, Boyden ES, Buzsáki G, Ramirez JM, Jones AR, Svoboda K, Han X, Turner EE, et al. 2012. A toolbox of Cre-dependent optogenetic transgenic mice for light-induced activation and silencing. *Nature Neuroscience* **15**:793-802. DOI: https://doi.org/10.1038/nn.3078, PMID: 22446880

Meurs A, Portelli J, Clinckers R, Balasubramaniam A, Michotte Y, Smolders I. 2012. Neuropeptide Y increases in vivo hippocampal extracellular glutamate levels through Y1 receptor activation. *Neuroscience Letters* **510**:143-147. DOI: https://doi.org/10.1016/j.neulet.2012.01.023, PMID: 22285263

Mishra A, Reynolds JP, Chen Y, Gourine AV, Rusakov DA, Attwell D. 2016. Astrocytes mediate neurovascular signaling to capillary pericytes but not to arterioles. *Nature Neuroscience* **19**:1619-1627. DOI: https://doi.org/10.1038/nn.4428, PMID: 27775719

Miyata N, Taniguchi K, Seki T, Ishimoto T, Sato-Watanabe M, Yasuda Y, Doi M, Kametani S, Tomishima Y, Ueki T, Sato M, Kameo K. 2001. HET0016, a potent and selective inhibitor of 20-HETE synthesizing enzyme. *British Journal of Pharmacology* **133**:325-329. DOI: https://doi.org/10.1038/sj.bjp.0704101, PMID: 11375247

Mulligan SJ, MacVicar BA. 2004. Calcium transients in astrocyte endfeet cause cerebrovascular constrictions. *Nature* **431**:195-199. DOI: https://doi.org/10.1038/nature02827, PMID: 15356633

Niwa K, Araki E, Morham SG, Ross ME, Iadecola C. 2000. Cyclooxygenase-2 contributes to functional hyperemia in whisker-barrel cortex. *The Journal of Neuroscience* **20**:763-770. DOI: https://doi.org/10.1523/JNEUROSCI.20-02-00763.2000, PMID: 10632605

O'Herron P, Chhatbar PY, Levy M, Shen Z, Schramm AE, Lu Z, Kara P. 2016. Neural correlates of single-vessel haemodynamic responses in vivo. *Nature* **534**:378-382. DOI: https://doi.org/10.1038/nature17965, PMID: 27281215

Oversø Hansen P, Kringelholt S, Simonsen U, Bek T. 2015. Hypoxia-induced relaxation of porcine retinal arterioles in vitro depends on inducible NO synthase and EP4 receptor stimulation in the perivascular retina. *Acta Ophthalmologica* **93**:457-463. DOI: https://doi.org/10.1111/aos.12669, PMID: 25619924

Palop JJ, Chin J, Roberson ED, Wang J, Thwin MT, Bien-Ly N, Yoo J, Ho KO, Yu GQ, Kreitzer A, Finkbeiner S, Noebels JL, Mucke L. 2007. Aberrant excitatory neuronal activity and compensatory remodeling of inhibitory hippocampal circuits in mouse models of Alzheimer's disease. *Neuron* **55**:697-711. DOI: https://doi.org/10.1016/j.neuron.2007.07.025, PMID: 17785178

Palop JJ, Mucke L. 2010. Amyloid-beta-induced neuronal dysfunction in Alzheimer's disease: from synapses toward neural networks. *Nature Neuroscience* **13**:812-818. DOI: https://doi.org/10.1038/nn.2583, PMID: 20581818

Pasinetti GM, Aisen PS. 1998. Cyclooxygenase-2 expression is increased in frontal cortex of Alzheimer's disease brain. *Neuroscience* **87**:319-324. DOI: https://doi.org/10.1016/s0306-4522(98)00218-8, PMID: 9740394

Pellerin L, Magistretti PJ. 1994. Glutamate uptake into astrocytes stimulates aerobic glycolysis: a mechanism coupling neuronal activity to glucose utilization. *PNAS* **91**:10625-10629. DOI: https://doi.org/10.1073/pnas.91.22.10625, PMID: 7938003

Perrenoud Q, Geoffroy H, Gauthier B, Rancillac A, Alfonsi F, Kessaris N, Rossier J, Vitalis T, Gallopin T. 2012. Characterization of Type I and Type II nNOS-expressing interneurons in the barrel cortex of mouse. *Frontiers in Neural Circuits* **6**:36. DOI: https://doi.org/10.3389/fncir.2012.00036, PMID: 22754499

Pham C, Komaki Y, Deàs-Just A, Le Gac B, Mouffle C, Franco C, Chaperon A, Vialou V, Tsurugizawa T, Cauli B, Li D. 2024. Astrocyte aquaporin mediates a tonic water efflux maintaining brain homeostasis. *eLife* **13**:RP95873. DOI: https://doi.org/10.7554/eLife.95873, PMID: 39508543

Rancillac A, Rossier J, Guille M, Tong XK, Geoffroy H, Amatore C, Arbault S, Hamel E, Cauli B. 2006. Glutamatergic control of microvascular tone by distinct gaba neurons in the cerebellum. *The Journal of Neuroscience* **26**:6997-7006. DOI: https://doi.org/10.1523/JNEUROSCI.5515-05.2006, PMID: 16807329

R Development Core Team. 2023. R: A language and environment for statistical computing. Vienna, Austria. R Foundation for Statistical Computing. https://www.r-project.org

Rosehart AC, Longden TA, Weir N, Fontaine JT, Joutel A, Dabertrand F. 2021. Prostaglandin $E_2$ dilates intracerebral arterioles when applied to capillaries: implications for small vessel diseases. *Frontiers in Aging Neuroscience* **13**:695965. DOI: https://doi.org/10.3389/fnagi.2021.695965, PMID: 34483880

Rudolf K, Eberlein W, Engel W, Wieland HA, Willim KD, Entzeroth M, Wienen W, Beck-Sickinger AG, Doods HN. 1994. The first highly potent and selective non-peptide neuropeptide Y Y1 receptor antagonist: BIBP3226. *European Journal of Pharmacology* **271**:R11-R13. DOI: https://doi.org/10.1016/0014-2999(94)90822-2

Rungta RL, Osmanski BF, Boido D, Tanter M, Charpak S. 2017. Light controls cerebral blood flow in naive animals. *Nature Communications* **8**:14191. DOI: https://doi.org/10.1038/ncomms14191, PMID: 28139643

Rungta RL, Chaigneau E, Osmanski BF, Charpak S. 2018. Vascular compartmentalization of functional hyperemia from the synapse to the pia. *Neuron* **99**:362-375.. DOI: https://doi.org/10.1016/j.neuron.2018.06.012, PMID: 29937277

Rungta RL, Zuend M, Aydin A, Weber B, Charpak S, Boido D. 2021. Vascular arbors in layer II / III somatosensory cortex. *Commun Biol* **94**:855. DOI: https://doi.org/10.1038/s42003-021-02382-w
Sada N, Lee S, Katsu T, Otsuki T, Inoue T. 2015. Epilepsy treatment targeting LDH enzymes with a stiripentol analog to treat epilepsy. *Science* **347**:1362–1367. DOI: https://doi.org/10.1126/science.aaa1299, PMID: 25792327

Säfholm J, Dahlén SE, Adner M. 2013a. Antagonising EP1 and EP2 receptors reveal that the TP receptor mediates a component of antigen-induced contraction of the guinea pig trachea. *European Journal of Pharmacology* **718**:277–282. DOI: https://doi.org/10.1016/j.ejphar.2013.08.021, PMID: 24012781

Säfholm J, Dahlén SE, Delin I, Maxey K, Stark K. Cardell LO, Adner M. 2013b. PGE2 maintains the tone of the guinea pig trachea through a balance between activation of contractile EP1 receptors and relaxant EP2 receptors. *British Journal of Pharmacology* **168**:794–806. DOI: https://doi.org/10.1111/j.1476-5381.2012.02189.x, PMID: 22934927

Sang N, Zhang J, Marcheselli V, Bazan NG, Chen C. 2005. Postsynaptically synthesized prostaglandin E2 (PGE2) modulates hippocampal synaptic transmission via a presynaptic PGE2 EP2 receptor. *The Journal of Neuroscience* **25**:9858–9870. DOI: https://doi.org/10.1523/JNEUROSCI.2392-05.2005, PMID: 16251433

Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, Preibisch S, Rueden C, Saalfeld S, Schmid B, Tinevez JY, White DJ, Hartenstein V, Eliceiri K, Tomancak P, Cardona A. 2012. Fiji: an open-source platform for biological-image analysis. *Nature Methods* **9**:676–682. DOI: https://doi.org/10.1038/nmeth.2019, PMID: 22743772

Schmid F, Barrett M, Jenny P, Weber B. 2019. Vascular density and distribution in neocortex. *NeuroImage* **197**:792–805. DOI: https://doi.org/10.1016/j.neuroimage.2017.06.046, PMID: 28669910

Scott NA, Murphy TH. 2012. Hemodynamic responses evoked by neuronal stimulation via channelrhodopsin-2 can be independent of intracortical glutamatergic synaptic transmission. *PLOS ONE* **7**:e29859. DOI: https://doi.org/10.1371/journal.pone.0029859, PMID: 22253807

Sharif NA, Klimko PG. 2019. Prostaglandin FP receptor antagonists: discovery, pharmacological characterization and therapeutic utility. *British Journal of Pharmacology* **176**:1059–1078. DOI: https://doi.org/10.1111/bph.14335, PMID: 29679483

Shmuel A, Yacoub E, Pfeuffer J, Van de Moortele PF, Adriany G, Hu X, Ugurbil K. 2002. Sustained negative BOLD, blood flow and oxygen consumption response and its coupling to the positive response in the human brain. *Neuron* **36**:1195–1210. DOI: https://doi.org/10.1016/s0896-6273(02)01061-9, PMID: 12495632

Smetters D, Majewska A, Yuste R. 1999. Detecting action potentials in neuronal populations with calcium imaging. *Methods* **18**:215–221. DOI: https://doi.org/10.1006/meth.1999.0774, PMID: 10356353

Smith CJ, Zhang Y, Koboldt CM, Muhammad J, Zweifel BS, Shaffer A, Talley JJ, Masferrer JL, Seibert K, Isakson PC. 1998. Pharmacological analysis of cyclooxygenase-1 in inflammation. *PNAS* **95**:13313–13318. DOI: https://doi.org/10.1073/pnas.95.22.13313, PMID: 9789085

Smith SJ, Sümbül U, Graybuck LT, Collman F, Seshamani S, Gala R, Gliko O, Elabbady L, Miller JA, Bakken TE, Rossier J, Yao Z, Lein E, Zeng H, Tasic B, Hawrylycz M. 2019. Single-cell transcriptomic evidence for dense intracortical neuropeptide networks. *eLife* **8**:e7889. DOI: https://doi.org/10.7554/eLife.47889, PMID: 31710287

Sun QQ, Baraban SC, Prince DA, Huguenard JR. 2003. Target-specific neuropeptide Y-ergic synaptic inhibition and its network consequences within the mammalian thalamus. *The Journal of Neuroscience* **23**:9639–9649. DOI: https://doi.org/10.1523/JNEUROSCI.23-29-09639.2003, PMID: 14573544

Sun W, McConnell E, Pare JF, Xu Q, Chen M, Peng W, Lovatt D, Han X, Smith Y, Nedergaard M. 2013. Glutamate-dependent neuroglial calcium signaling differs between young and adult brain. *Science* **339**:197–200. DOI: https://doi.org/10.1126/science.1226740, PMID: 23307741

Takei S, Hasegawa-Ishii S, Uekawa A, Chiba Y, Umegaki H, Hosokawa M, Woodward DF, Watanabe K, Shimada A. 2012. Immunohistochemical demonstration of increased prostaglandin $F_2\alpha$ levels in the rat hippocampus following kainic acid-induced seizures. *Neuroscience* **218**:295–304. DOI: https://doi.org/10.1016/j.neuroscience.2012.05.013, PMID: 22609937

Takemiya T, Matsumura K, Yamagata K. 2007. Roles of prostaglandin synthesis in excitotoxic brain diseases. *Neurochemistry International* **51**:112–120. DOI: https://doi.org/10.1016/j.neuint.2007.05.009, PMID: 17629358

Tammaro P, Smith AL, Hutchings SR, Smirnov SV. 2004. Pharmacological evidence for a key role of voltage-gated K+ channels in the function of rat aortic smooth muscle cells. *British Journal of Pharmacology* **143**:303–317. DOI: https://doi.org/10.1038/sj.bjp.0705957, PMID: 15326038

Tasic B, Menon V, Nguyen TN, Kim TK, Jarsky T, Yao Z, Levi B, Gray LT, Sorensen SA, Dolbeare T, Bertagnolli D, Goldy J, Shapovalova N, Parry S, Lee C, Smith K, Bernard A, Madisen L, Sunkin SM, Hawrylycz M, et al. 2016. Adult mouse cortical cell taxonomy revealed by single cell transcriptomics. *Nature Neuroscience* **19**:335–346. DOI: https://doi.org/10.1038/nn.4216, PMID: 26727548

Thévenaz P, Ruttimann UE, Unser M. 1998. A pyramid approach to subpixel registration based on intensity. *IEEE Transactions on Image Processing* **7**:27–41. DOI: https://doi.org/10.1109/83.650848, PMID: 18267377

Tournissac M, Boido D, Omnès M, Houssen YG, Ciobanu L, Charpak S. 2022. Cranial window for longitudinal and multimodal imaging of the whole mouse cortex. *Neurophotonics* **9**:031921. DOI: https://doi.org/10.1117/1.NPh.9.3.031921, PMID: 36159711

Tran CHT, George AG, Teskey GC, Gordon GR. 2020. Seizures elevate glyovascular unit Ca2+ and cause sustained vasoconstriction. *JCI Insight* **5**:e136469. DOI: https://doi.org/10.1172/jci.insight.136469, PMID: 33004688

Tsai CJ, Nagata T, Liu CY, Suganuma T, Kanda T, Miyazaki T, Liu K, Saitoh T, Nagase H, Lazarus M, Vogt KE, Yanagisawa M, Hayashi Y. 2021. Cerebral capillary blood flow upsurge during REM sleep is mediated by A2a receptors. *Cell Reports* **36**:109558. DOI: https://doi.org/10.1016/j.celrep.2021.109558, PMID: 34407410
Tsurumaki T, Yamaguchi T, Higuchi H. 2002. Marked neuropeptide Y-induced contractions via NPY-Y1 receptor and its desensitization in rat veins. *Vascular Pharmacology* **39**:325-333. DOI: https://doi.org/10.1016/s1537-1891(03)00044-2, PMID: 14567071

Tsurumaki T, Muraoka O, Yamaguchi T, Higuchi H. 2003. Neuropeptide Y-induced contraction and its desensitization through the neuropeptide Y receptor subtype in several rat veins. *Journal of Cardiovascular Pharmacology* **41 Suppl 1**:S23-S27 PMID: 12688392.

Uhlirova H, Kılıç K, Tian P, Thunemann M, Desjardins M, Saisan PA, Sakadžić S, Ness TV, Mateo C, Cheng Q, Weldy KL, Razoux F, Vandenberghe M, Cremonesi JA, Ferri CG, Nizar K, Sridhar VB, Steed TC, Abashin M, Fainman Y, et al. 2016. Cell type specificity of neurovascular coupling in cerebral cortex. eLife **5**:e4315. DOI: https://doi.org/10.7554/eLife.14315, PMID: 27244241

Voutsinos-Porche B, Bonvento G, Tanaka K, Steiner P, Welker E, Chatton JY, Magistretti PJ, Pellerin L. 2003. Glial glutamate transporters mediate a functional metabolic crosstalk between neurons and astrocytes in the mouse developing cortex. *Neuron* **37**:275-286. DOI: https://doi.org/10.1016/s0896-6273(02)01170-4, PMID: 12546822

Wagner L, Wolf R, Zeitschel U, Rossner S, Petersén Å, Leavitt BR, Kästner F, Rothermundt M, Gärtner U-T, Gündel D, Schlenzig D, Frerker N, Schade J, Manhart S, Rahfeld J-U, Demuth H-U, von Hörsten S. 2015. Proteolytic degradation of neuropeptide Y (NPY) from head to toe: Identification of novel NPY-cleaving peptidases and potential drug interactions in CNS and Periphery. *Journal of Neurochemistry* **135**:1019-1037. DOI: https://doi.org/10.1111/jnc.13378, PMID: 26442809

Yamagata K, Andreasson KI, Kaufmann WE, Barnes CA, Worley PF. 1993. Expression of a mitogen-inducible cyclooxygenase in brain neurons: regulation by synaptic activity and glucocorticoids. *Neuron* **11**:371-386. DOI: https://doi.org/10.1016/0896-6273(93)90192-t, PMID: 8352945

Yamagata K, Matsumura K, Inoue W, Shiraki T, Suzuki K, Yasuda S, Sugiura H, Cao C, Watanabe Y, Kobayashi S. 2001. Coexpression of microsomal-type prostaglandin E synthase with cyclooxygenase-2 in brain endothelial cells of rats during endotoxin-induced fever. *The Journal of Neuroscience* **21**:2669-2677. DOI: https://doi.org/10.1523/JNEUROSCI.21-08-02669.2001, PMID: 11306620

Zeisel A, Muñoz-Manchado AB, Codeluppi S, Lönnerberg P, La Manno G, Juréus A, Marques S, Munguba H, He L, Betsholtz C, Rolny C, Castelo-Branco G, Hjerling-Leffler J, Linnarsson S. 2015. Brain structure: cell types in the mouse cortex and hippocampus revealed by single-cell RNA-seq. *Science* **347**:1138-1142. DOI: https://doi.org/10.1126/science.aaa1934, PMID: 25700174

Zhang DY, Raichle ME. 2010. Disease and the brain's dark energy. *Nature Reviews Neurology* **6**:15-28. DOI: https://doi.org/10.1038/nrneurol.2009.198

Zhang D, Ruan J, Peng S, Li J, Hu X, Zhang Y, Zhang T, Ge Y, Zhu Z, Xiao X, Zhu Y, Li X, Li T, Zhou L, Gao Q, Zheng G, Zhao B, Li X, Zhu Y, Wu J, et al. 2024. Synaptic-like transmission between neural axons and arteriolar smooth muscle cells drives cerebral neurovascular coupling. *Nature Neuroscience* **27**:232-248. DOI: https://doi.org/10.1038/s41593-023-01515-0, PMID: 38168932

Zhou N, Gordon GRJ, Feighan D, MacVicar BA. 2010. Transient swelling, acidification, and mitochondrial depolarization occurs in neurons but not astrocytes during spreading depression. *Cerebral Cortex* **20**:2614-2624. DOI: https://doi.org/10.1093/cercor/bhq018, PMID: 20176688

Zonta M, Angulo MC, Gobbo S, Rosengarten B, Hossmann KA, Pozzan T, Carmignoto G. 2003. Neuron-to-astrocyte signaling is central to the dynamic control of brain microcirculation. *Nature Neuroscience* **6**:43-50. DOI: https://doi.org/10.1038/nn980, PMID: 12469126
# Appendix 1

## Appendix 1—key resources table

<table><thead><tr><th>Reagent type (species) or resource</th><th>Designation</th><th>Source or reference</th><th>Identifiers</th><th>Additional information</th></tr></thead><tbody><tr><td>Strain, strain background (Mus musculus, male and female)</td><td>C57BL/6RJ, Wild type</td><td>Janvier Labs</td><td>C57BL/6RJ</td><td></td></tr><tr><td>Strain, strain background (Mus musculus, male and female)</td><td>B6.129P2- Emx1<sup>tm1(cre)Kj</sup>/J, Emx1<sup>Cre/Cre</sup></td><td>PMID:12151506; Gorski et al., 2002</td><td>RRID:IMSR_JAX:005628</td><td></td></tr><tr><td>Strain, strain background (Mus musculus, male and female)</td><td>B6.129P2- Gt(ROSA)26Sor<sup>tm32(CAG-COP4*H134R/EYFP)Hze</sup>, Gt(ROSA)26Sor ChR2(H134R)-EYFP/ ChR2(H134R)-EYFP</td><td>PMID:22446880; Madisen et al., 2012</td><td>RRID:IMSR_JAX:024109</td><td></td></tr><tr><td>Sequence-based reagent</td><td>Slc17a7 external sense PCR primer (vGluT1)</td><td>PMID:23565079; Cabezas et al., 2013</td><td></td><td>GGCTCCTTTTTCTGGGGCTAC</td></tr><tr><td>Sequence-based reagent</td><td>Slc17a7 external antisense PCR primer (vGluT1)</td><td>PMID:23565079; Cabezas et al., 2013</td><td></td><td>CCAGCCGACTCCGTTCTAAG</td></tr><tr><td>Sequence-based reagent</td><td>Slc17a7 internal sense PCR primer (vGluT1)</td><td>PMID:23565079; Cabezas et al., 2013</td><td></td><td>ATTCGCAGCCAACAGGGTCT</td></tr><tr><td>Sequence-based reagent</td><td>Slc17a7 internal antisense PCR primer (vGluT1)</td><td>PMID:23565079; Cabezas et al., 2013</td><td></td><td>TGGCAAGCAGGGTATGTGAC</td></tr><tr><td>Sequence-based reagent</td><td>Gad2 external sense PCR primer (GAD 65)</td><td>PMID:19295167; Karagiannis et al., 2009</td><td></td><td>CCAAAAGTTCACGGGCGG</td></tr><tr><td>Sequence-based reagent</td><td>Gad2 external antisense PCR primer (GAD 65)</td><td>PMID:19295167; Karagiannis et al., 2009</td><td></td><td>TCCTCCAGATTTTGCGGTTG</td></tr><tr><td>Sequence-based reagent</td><td>Gad2 internal sense PCR primer (GAD 65)</td><td>PMID:22754499; Perrenoud et al., 2012</td><td></td><td>CACCTGCGACCAAAAACCCT</td></tr><tr><td>Sequence-based reagent</td><td>Gad2 internal antisense PCR primer (GAD 65)</td><td>PMID:22754499; Perrenoud et al., 2012</td><td></td><td>GATTTTGCGGTTGGTCTGCC</td></tr><tr><td>Sequence-based reagent</td><td>Gad1 external sense PCR primer (GAD 67)</td><td>PMID:12196560; Férézou et al., 2002</td><td></td><td>TACGGGGTTCGCACAGGTC</td></tr><tr><td>Sequence-based reagent</td><td>Gad1 external antisense PCR primer (GAD 67)</td><td>PMID:12196560; Cabezas et al., 2013</td><td></td><td>CCCAGGCAGCATCCACAT</td></tr><tr><td>Sequence-based reagent</td><td>Gad1 internal sense PCR primer (GAD 67)</td><td>PMID:23565079; Cabezas et al., 2013</td><td></td><td>CCCAGAAGTGAAGACAAAAGGC</td></tr><tr><td>Sequence-based reagent</td><td>Gad1 internal antisense PCR primer (GAD 67)</td><td>PMID:23565079; Cabezas et al., 2013</td><td></td><td>AATGCTCCGTAAACAGTCGTGC</td></tr><tr><td>Sequence-based reagent</td><td>Ptgs1 external sense PCR primer (COX-1)</td><td>This paper</td><td></td><td>ATCCCTGTTGTTACTATCCGTGC</td></tr><tr><td>Sequence-based reagent</td><td>Ptgs1 external antisense PCR primer (COX-1)</td><td>This paper</td><td></td><td>TGTGGGGCAGTCTTTGGGTA</td></tr><tr><td>Sequence-based reagent</td><td>Ptgs1 internal sense PCR primer (COX-1)</td><td>This paper</td><td></td><td>AGGGTGTCTGTGTCCGCTTT</td></tr></tbody></table>

Appendix 1 Continued on next page
Appendix 1 Continued

<table><thead><tr><th>Reagent type (species) or resource</th><th>Designation</th><th>Source or reference</th><th>Identifiers</th><th>Additional information</th></tr></thead><tbody><tr><td>Sequence-based reagent</td><td>Ptgs1 internal antisense PCR primer (COX-1)</td><td>This paper</td><td></td><td>GGCTGGGGATAAGGTTGGAC</td></tr><tr><td>Sequence-based reagent</td><td>Ptgs2 external sense PCR primer (COX-2)</td><td>PMID:21734275; Lecrux et al., 2011</td><td></td><td>CTGAAGCCCACCCCAAACAC</td></tr><tr><td>Sequence-based reagent</td><td>Ptgs2 external antisense PCR primer (COX-2)</td><td>PMID:29985318; Devienne et al., 2018</td><td></td><td>CCTTATTTCCCTTCACACCCAT</td></tr><tr><td>Sequence-based reagent</td><td>Ptgs2 internal sense PCR primer (COX-2)</td><td>PMID:29985318; Devienne et al., 2018</td><td></td><td>AACAACATCCCCTTCCTGCG</td></tr><tr><td>Sequence-based reagent</td><td>Ptgs2 internal antisense PCR primer (COX-2)</td><td>PMID:29985318; Devienne et al., 2018</td><td></td><td>TGGGAGTTGGGCAGTCATCT</td></tr><tr><td>Sequence-based reagent</td><td>Ptges external sense PCR primer (mPGES1)</td><td>This paper</td><td></td><td>GCCTGGTGATGGAGAGCG</td></tr><tr><td>Sequence-based reagent</td><td>Ptges external antisense PCR primer (mPGES1)</td><td>This paper</td><td></td><td>GGAGCGAAGGCGTGGGTT</td></tr><tr><td>Sequence-based reagent</td><td>Ptges internal sense PCR primer (mPGES1)</td><td>This paper</td><td></td><td>AGATGAGGCTGCGGAAGAAG</td></tr><tr><td>Sequence-based reagent</td><td>Ptges internal antisense PCR primer (mPGES1)</td><td>This paper</td><td></td><td>CACGAAGCCGAGGAAGAGGA</td></tr><tr><td>Sequence-based reagent</td><td>Ptges2 external sense PCR primer (mPGES1)</td><td>This paper</td><td></td><td>CGACTTCCACTCCCTGCC</td></tr><tr><td>Sequence-based reagent</td><td>Ptges2 external antisense PCR primer (mPGES2)</td><td>This paper</td><td></td><td>CATCTCCTCCGTCCTGGCTT</td></tr><tr><td>Sequence-based reagent</td><td>Ptges2 internal sense PCR primer (mPGES2)</td><td>This paper</td><td></td><td>GAGGTGAATCCCGTGAGAAGG</td></tr><tr><td>Sequence-based reagent</td><td>Ptges2 internal antisense PCR primer (mPGES2)</td><td>This paper</td><td></td><td>TTCCTTCCCGCCATACATCT</td></tr><tr><td>Sequence-based reagent</td><td>Ptges3 external sense PCR primer (cPGES)</td><td>This paper</td><td></td><td>TCCAAGCATAAAAGAACAGACAGA</td></tr><tr><td>Sequence-based reagent</td><td>Ptges3 external antisense PCR primer (cPGES)</td><td>This paper</td><td></td><td>TGGCATCTTTTCATCATCACTGTC</td></tr><tr><td>Sequence-based reagent</td><td>Ptges3 internal sense PCR primer (cPGES)</td><td>This paper</td><td></td><td>TAACAAAGGAAAGGGCAAAGC</td></tr><tr><td>Sequence-based reagent</td><td>Ptges3 internal antisense PCR primer (cPGES)</td><td>This paper</td><td></td><td>CATCATCTGCTCCATCTACTTCTG</td></tr><tr><td>Sequence-based reagent</td><td>Prxl2b external sense PCR primer (PM-PGFS)</td><td>This paper</td><td></td><td>AGGAGTTTCTGGATGGTGGTTAC</td></tr><tr><td>Sequence-based reagent</td><td>Prxl2b external antisense PCR primer (PM-PGFS)</td><td>This paper</td><td></td><td>CACCTCCCACACACCTCTTCAT</td></tr><tr><td>Sequence-based reagent</td><td>Prxl2b internal sense PCR primer (PM-PGFS)</td><td>This paper</td><td></td><td>ACCTGTTCGTGATGTAGCCTCC</td></tr><tr><td>Sequence-based reagent</td><td>Prxl2b internal antisense PCR primer (PM-PGFS)</td><td>This paper</td><td></td><td>CTGGGGTGGCTTGCTGGA</td></tr><tr><td>Sequence-based reagent</td><td>Akr1b1 external sense PCR primer (Akr1b3)</td><td>This paper</td><td></td><td>CAGAATGAGAAGGAGGTGGGA</td></tr><tr><td>Sequence-based reagent</td><td>Akr1b1 external antisense PCR primer (Akr1b3)</td><td>This paper</td><td></td><td>TTGAAGTTGGAGACACCGATTG</td></tr></tbody></table>

Appendix 1 Continued on next page
Appendix 1 Continued

<table><thead><tr><th>Reagent type (species) or resource</th><th>Designation</th><th>Source or reference</th><th>Identifiers</th><th>Additional information</th></tr></thead><tbody><tr><td>Sequence-based reagent</td><td>Akr1b1 internal sense PCR primer (Akr1b3)</td><td>This paper</td><td></td><td>CAAGGAGCAGGTGGTGAAGC</td></tr><tr><td>Sequence-based reagent</td><td>Akr1b1 internal antisense PCR primer (Akr1b3)</td><td>This paper</td><td></td><td>CATAGCCGTCCAAGTGTCCA</td></tr><tr><td>Sequence-based reagent</td><td>Cbr1 external sense PCR primer (CBR1)</td><td>This paper</td><td></td><td>AACCCGCAGAGCATTCGC</td></tr><tr><td>Sequence-based reagent</td><td>Cbr1 external antisense PCR primer (CBR1)</td><td>This paper</td><td></td><td>GCCAACCTTCTTCCGCAT</td></tr><tr><td>Sequence-based reagent</td><td>Cbr1 internal sense PCR primer (CBR1)</td><td>This paper</td><td></td><td>CAATGACGACACCCCCTTCC</td></tr><tr><td>Sequence-based reagent</td><td>Cbr1 internal antisense PCR primer (CBR1)</td><td>This paper</td><td></td><td>CTCCTCTGTGATGGTCTCGCTT</td></tr><tr><td>Chemical compound, drug</td><td>Rhod-2</td><td>Cayman chemical</td><td>20777</td><td></td></tr><tr><td>Chemical compound, drug</td><td>9,11-dideoxy-9α,11α-methanoepoxy prostaglandin F2α</td><td>Enzo</td><td>BML-PG023</td><td></td></tr><tr><td>Chemical compound, drug</td><td>Tetrodotoxin</td><td>Latoxan</td><td>L8503</td><td></td></tr><tr><td>Chemical compound, drug</td><td>D-(-)-2-amino-5-phosphonopentanoic acid</td><td>Hello Bio</td><td>HB0225</td><td></td></tr><tr><td>Chemical compound, drug</td><td>6,7-dinitroquinoxaline-2,3-dione</td><td>Hello Bio</td><td>HB0262</td><td></td></tr><tr><td>Chemical compound, drug</td><td>LY367385</td><td>Hello Bio</td><td>HB0398</td><td></td></tr><tr><td>Chemical compound, drug</td><td>2-methyl-6-(phenylethynyl)pyridine</td><td>Hello Bio</td><td>HB0426</td><td></td></tr><tr><td>Chemical compound, drug</td><td>Indomethacin</td><td>Sigma-Aldrich</td><td>I7378</td><td></td></tr><tr><td>Chemical compound, drug</td><td>SC-560</td><td>Sigma-Aldrich</td><td>S2064</td><td></td></tr><tr><td>Chemical compound, drug</td><td>NS-398</td><td>Enzo</td><td>BML-EI261</td><td></td></tr><tr><td>Chemical compound, drug</td><td>ONO-8130</td><td>Tocris</td><td>5406</td><td></td></tr><tr><td>Chemical compound, drug</td><td>L-798,106</td><td>Cayman chemical</td><td>11129</td><td></td></tr><tr><td>Chemical compound, drug</td><td>AL8810</td><td>Cayman chemical</td><td>16735</td><td></td></tr><tr><td>Chemical compound, drug</td><td>PGE2</td><td>Hello Bio</td><td>HB3460</td><td></td></tr><tr><td>Chemical compound, drug</td><td>Sulprostone</td><td>Cayman chemical</td><td>14765</td><td></td></tr><tr><td>Chemical compound, drug</td><td>BIBP3226</td><td>Tocris</td><td>2707</td><td></td></tr><tr><td>Chemical compound, drug</td><td>HET0016</td><td>Merck</td><td>SML2416</td><td></td></tr><tr><td>Chemical compound, drug</td><td>paxilline</td><td>Tocris</td><td>2006</td><td></td></tr></tbody></table>

Appendix 1 Continued on next page
Appendix 1 Continued

<table><thead><tr><th>Reagent type (species) or resource</th><th>Designation</th><th>Source or reference</th><th>Identifiers</th><th>Additional information</th></tr></thead><tbody><tr><td>Chemical compound, drug</td><td>Dithiothreitol</td><td>VWR</td><td>443852 A</td><td></td></tr><tr><td>Chemical compound, drug</td><td>Primer "random"</td><td>Roche</td><td>11034731001</td><td></td></tr><tr><td>Chemical compound, drug</td><td>dNTPs</td><td>GE Healthcare Life Sciences</td><td>28-4065-52</td><td></td></tr><tr><td>Chemical compound, drug</td><td>Mineral Oil</td><td>Sigma-Aldrich</td><td>M5904</td><td></td></tr><tr><td>Chemical compound, drug</td><td>RNasin Ribonuclease Inhibitors</td><td>Promega</td><td>N2511</td><td></td></tr><tr><td>Chemical compound, drug</td><td>SuperScript II Reverse Transcriptase</td><td>Invitrogen</td><td>18064014</td><td></td></tr><tr><td>Chemical compound, drug</td><td>Taq DNA Polymerase</td><td>Qiagen</td><td>201205</td><td></td></tr><tr><td>Software, algorithm</td><td>Pclamp v 10.2</td><td>Molecular Devices</td><td>RRID:SCR_011323</td><td></td></tr><tr><td>Software, algorithm</td><td>Matlab v 2018b</td><td>MathWorks</td><td>RRID:SCR_001622</td><td></td></tr><tr><td>Software, algorithm</td><td>GraphPad Prism v 7</td><td>GraphPad</td><td>RRID:SCR_002798</td><td></td></tr><tr><td>Software, algorithm</td><td>ImagingWorkbench v 6.1</td><td>INDEC Systems</td><td></td><td></td></tr><tr><td>Software, algorithm</td><td>FIJI</td><td>PMID:22743772; Schindelin et al., 2012</td><td>RRID:SCR_002285</td><td></td></tr><tr><td>Software, algorithm</td><td>R v 4.3.0</td><td>R Core Team</td><td>RRID:SCR_001905</td><td></td></tr><tr><td>Other</td><td>Vibratome</td><td>Leica</td><td>VT1000S<br>RRID:SCR_016495</td><td></td></tr><tr><td>Other</td><td>Upright microscope</td><td>Olympus</td><td>BX51WI</td><td></td></tr><tr><td>Other</td><td>Dual port module</td><td>Olympus</td><td>WI-DPMC</td><td></td></tr><tr><td>Other</td><td>60 x Objective</td><td>Olympus</td><td>LUMPlan FI /IR<br>60 x/0.90 W</td><td></td></tr><tr><td>Other</td><td>40 x Objective</td><td>Olympus</td><td>LUMPlan FI /IR<br>40 x/0.80 W</td><td></td></tr><tr><td>Other</td><td>sCMOS camera</td><td>Hamamatsu</td><td>ORCA-Flash4.0</td><td></td></tr><tr><td>Other</td><td>Axopatch 200B</td><td>Molecular Devices</td><td>RRID:SCR_018866</td><td></td></tr><tr><td>Other</td><td>Digidata 1440 A</td><td>Molecular Devices</td><td>RRID:SCR_021038</td><td></td></tr><tr><td>Other</td><td>S900 stimulator</td><td>Dagan corporation</td><td></td><td></td></tr><tr><td>Other</td><td>pE-2</td><td>CoolLED</td><td></td><td></td></tr><tr><td>Other</td><td>Excitation filter</td><td>Semrock</td><td>HC 392/474/554/635</td><td></td></tr><tr><td>Other</td><td>Dichroic mirror</td><td>Semrock</td><td>BS 409/493/573/652</td><td></td></tr><tr><td>Other</td><td>Emission filter</td><td>Semrock</td><td>HC 432/515/595/730</td><td></td></tr><tr><td>Other</td><td>780 nm Collimated LED</td><td>Thorlabs</td><td>M780L3-C1</td><td></td></tr><tr><td>Other</td><td>Dodt Gradient Contrast</td><td>Luigs and Neumann</td><td>200-100 200 0155</td><td></td></tr><tr><td>Other</td><td>Beam splitter</td><td>Semrock</td><td>725 DCSPXR</td><td></td></tr><tr><td>Other</td><td>Analogic CCD camera</td><td>Sony</td><td>XC ST-70 CE</td><td></td></tr></tbody></table>