Photostimulation setup:

The excitation and stimulation light pass through a FV30-NDM690 dichroic mirror with two notch filters, at 458 nm and 552 nm, to excite TexasRed, EYFP and ChR2 within the mouse. The emitted light passes through the objective, is reflected off the FV30-NDM690 dichroic mirror, and passes through a 650 nm barrier filter before reaching a 570 nm long pass filter (LPF) separating emitted light from EYFP and TexasRed, which respectively pass through 495-540 nm and 575-630 nm barrier filters to be collected via GaAsP detectors.

Computational analysis pipeline:

A. The stacks of 2PFM slices were registered using ANTS rigid registration and aligned to the reference time point. B. Images were upsampled using bicubic interpolation to an isotropic resolution of 0.99×0.99×0.99 μm. C. An ensemble of UNETR deep learning models with dropout generated segmentation masks at each time point, producing probability maps. D. The mean and standard deviation of the probability of each pixel being vasculature were computed and used to create binary vascular segmentation masks. E. The union over the vascular segmentation masks for all time points was computed, and background pixel clusters within vessel masks were removed. F. The vascular segmentation mask was thinned down to centerlines and rendered as a graph, where edges were vessel segments connecting branch points (nodes). This skeleton was overlaid on the vasculature channel from which the neuron channel was subtracted. G. The plane orthogonal to the tangent to the vessel’s travel direction was computed every micrometre along the centerline. H,I. 1D signal intensity profiles at each centerline vertex were computed in the orthogonal plane every 10°. J. The boundary for each profile was placed at the minimum of the signal gradient for that signal intensity profile. K. The raw intensity image with the detected boundary points, where outlier boundary points (in green) were defined as points over 2 standard deviations from the mean and excluded. L. Visualization of the changes in vertex-wise radii on a sample vascular network.

Model performance metrics:

The Dice, precision, recall, mean surface distance, and HD95 distance for the vascular (A) and neuron (B) channels. Each model was evaluated on the same test dataset composed of 9 images (250×507×507 μm each) from 6 mice. A Wilcoxon signed-rank test was used to compare the model’s performance on each performance metric for images from the test dataset. * p < 0.05, ** p < 0.005, and *** p < 0.0005. p-values were not adjusted.

Visual model comparison:

A. Raw images of the vascular channel with the neuron channel subtracted to facilitate vessel visualization. The first and last stacks in each row span from the cortical surface to 250 μm below the surface, while the middle stack spans from 250 μm below the surface to 500 μm below the surface. All images were from the test dataset, which was unseen during model training. B. Ground truth segmentation masks for the vasculature were generated by a rater who utilized ilastik-assisted manual segmentation. C. Ilastik predictions generated via a random forest model. D. Binary segmentation masks generated by an ensemble of 3D UNet Models. E. Binary segmentation masks generated by an ensemble of 3D UNETR models.

Estimation of simulated radii changes:

A. An image in the plane orthogonal to the local tangent to a capillary with the detected boundary (in blue) and with the estimated radius of 2.28 μm. On the right, this image was resized (upsampling, via bicubic interpolation, by 1.10 times) to simulate dilation. B. The plot shows correspondence between the estimated radius following scaling and the simulated level of scaling. C. An image in the plane orthogonal to the local tangent of a capillary with the detected boundary (in blue) and with the estimated radius of 3.65 μm. On the right, Gaussian noise with a sigma of 205.36 SU was added to the image. D. The estimated % change in the vessel’s radius after the addition of varying levels of Gaussian noise, demonstrating the robustness of the radius estimated to noise.

Vascular graph examples:

A. Baseline variability in vessel diameter estimated by the standard deviation of each vessel’s mean radius across baseline time frames. B. Mean change in the vessel radius induced by optogenetic stimulation. C. Mean change in the vertexwise radius, allowing the visualization of heterogeneity of radius changes within each vessel. D. Distance from each vertex to the closest pyramidal neuron. Each row corresponds to the vascular graph of a different mouse.

Vertex-wise radii along vessel lengths of a sample artery, capillary and venule at baseline vs. post-stimulation:

A. MIP of an artery, vein, and capillary segments before (left) and after (right) optogenetic stimulation with 458nm light at 1.1 mW/mm2. The artery and capillary dilated by 1.33±0.86 μm and 0.42±0.39 μm, respectively (for both p<1e-4, Mann-Whitney U test), whereas there was no significant change in the venular caliber upon photostimulation (p=0.22, Mann-Whitney U test). B. Estimates of the vertex-wise radius obtained along each of the three vessels’ centrelines, before and after stimulation. C. Vertex-wise radii changes in response to optogenetic stimulation. D. The vertex-wise distance from the vascular surface to the closest YFP-expressing neuron.

S1FL Vascular Network Morphological Properties

Optogenetic activation-induced changes in vessel-wise microvascular radii:

Capillary responses included both dilatations, shown in A. and constrictions, shown in B., with changes in the magnitude of the capillary response with increased photostimulation power. * p < 0.05, ** p < 0.005, and *** p < 0.0005. p-values were not adjusted. C. Changes to capillary radii are displayed in relation to the closest pyramidal neurons. The proportion of vessels constricting increased with the higher intensity of blue light stimulation, and constrictions tended to occur further away from pyramidal neurons than did dilations. D. Mean cortical depth of responding capillaries showed a tendency for dilators to be closer to the surface and for constrictors to be deeper in the tissue.

Microvascular Network Coordination Following Optogenetic Stimulation:

A. Graph representation of a vascular network of 425 vascular segments from a single image stack. Vessel segments are depicted as nodes of the graph; vascular segments that are joined at junctions are connected by edges. Nodes are coloured by the change in the mean vessel-wise radius following photostimulation with 458 nm light at 4.3 mW/mm2. B. Assortativity of photostimulation-induced changes in mean capillary radius increased with increasing photostimulation power. C. Photostimulation-induced changes in the efficiency of the capillary network. The capillary network efficiency changed by a median -0.16 PΩ-1 (IQR: -0.39 to 0.10 PΩ-1) in response to green light; -0.14 PΩ-1 (IQR: -0.55 to 0.27 PΩ-1) in response to lower intensity bluelight; and 0.22 PΩ-1 (IQR = -0.43;1.47 PΩ-1) in response to higher intensity blue light.. There was a significant increase (p = 0.03) in the capillary network efficiency post 458-nm light at 4.3 mW/mm2, when compared to that following the control green illumination. The measurements came from 72 paired acquisitions of 32 image stacks acquired in 17 mice (9M/8F). * p < 0.05, ** p < 0.005, and *** p < 0.0005. p-values were not adjusted.

Details of responder vessels

Attrition:

Flow chart of animal numbers at each step of the experiment.

Vascular Network Characteristics:

A. Probability density of the extracted mean radii of vascular segments. B. Probability density of the lengths of extracted vessel segments between branch points or terminal ends. C. Probability density of the mean vessel segment depths. D. Probability density of the depths of vessel branch points. Terminal nodes were excluded from this probability density. 12555 vessels, and 6421 vascular junctions from 17 mice (9M/8F) were used to estimate these PDFs.

Large Vessels Radius Changes Following Optogenetic Stimulation:

All vascular responses in large vessels (radius > 5μm) separated into dilators A. and constrictors B. showing an increase in the magnitude of response as stimulation power increases. C. Radius changes to vascular radii in relation to the closest pyramidal neurons (within 80 μm). D. Mean depth of responding large vessels.

Physiological Monitoring Data

Data Augmentations:

Model performance comparisons: