Insight into the fundamental trade-offs of diffusion MRI from polarization-sensitive optical coherence tomography in ex vivo human brain

In the first study comparing high angular resolution diffusion MRI (dMRI) in the human brain to axonal orientation measurements from polarization-sensitive optical coherence tomography (PSOCT), we compare the accuracy of orientation estimates from various dMRI sampling schemes and reconstruction methods. We find that, if the reconstruction approach is chosen carefully, single-shell dMRI data can yield the same accuracy as multi-shell data, and only moderately lower accuracy than a full Cartesian-grid sampling scheme. Our results suggest that current dMRI reconstruction approaches do not benefit substantially from ultra-high b-values or from very large numbers of diffusion-encoding directions. We also show that accuracy remains stable across dMRI voxel sizes of 1 mm or smaller but degrades at 2 mm, particularly in areas of complex white-matter architecture. We also show that, as the spatial resolution is reduced, axonal configurations in a dMRI voxel can no longer be modeled as a small set of distinct axon populations, violating an assumption that is sometimes made by dMRI reconstruction techniques. Our findings have implications for in vivo studies and illustrate the value of PSOCT as a source of ground-truth measurements of white-matter organization that does not suffer from the distortions typical of histological techniques.

Two-photon microscopic imaging of capillary red blood cell flux in mouse brain reveals vulnerability of cerebral white matter to hypoperfusion

Despite the importance of understanding the regulation of microvascular blood flow in white matter, no data on subcortical capillary blood flow parameters are available, largely due to the lack of appropriate imaging methods. To address this knowledge gap, we employed two-photon microscopy using a far-red fluorophore Alexa680 and photon-counting detection to measure capillary red blood cell (RBC) flux in both cerebral gray and white matter, in isoflurane-anesthetized mice. We have found that in control animals, baseline capillary RBC flux in the white matter was significantly higher than in the adjacent cerebral gray matter. In response to mild hypercapnia, RBC flux in the white matter exhibited significantly smaller fractional increase than in the gray matter. Finally, during global cerebral hypoperfusion, RBC flux in the white matter was reduced significantly in comparison to the controls, while RBC flux in the gray matter was preserved. Our results suggest that blood flow in the white matter may be less efficiently regulated when challenged by physiological perturbations as compared to the gray matter. Importantly, the blood flow in the white matter may be more susceptible to hypoperfusion than in the gray matter, potentially exacerbating the white matter deterioration in brain conditions involving global cerebral hypoperfusion.

Using the General Linear Model to Improve Performance in fNIRS Single Trial Analysis and Classification: A Perspective

Within a decade, single trial analysis of functional Near Infrared Spectroscopy (fNIRS) signals has gained significant momentum, and fNIRS joined the set of modalities frequently used for active and passive Brain Computer Interfaces (BCI). A great variety of methods for feature extraction and classification have been explored using state-of-the-art Machine Learning methods. In contrast, signal preprocessing and cleaning pipelines for fNIRS often follow simple recipes and so far rarely incorporate the available state-of-the-art in adjacent fields. In neuroscience, where fMRI and fNIRS are established neuroimaging tools, evoked hemodynamic brain activity is typically estimated across multiple trials using a General Linear Model (GLM). With the help of the GLM, subject, channel, and task specific evoked hemodynamic responses are estimated, and the evoked brain activity is more robustly separated from systemic physiological interference using independent measures of nuisance regressors, such as short-separation fNIRS measurements. When correctly applied in single trial analysis, e.g., in BCI, this approach can significantly enhance contrast to noise ratio of the brain signal, improve feature separability and ultimately lead to better classification accuracy. In this manuscript, we provide a brief introduction into the GLM and show how to incorporate it into a typical BCI preprocessing pipeline and cross-validation. Using a resting state fNIRS data set augmented with synthetic hemodynamic responses that provide ground truth brain activity, we compare the quality of commonly used fNIRS features for BCI that are extracted from (1) conventionally preprocessed signals, and (2) signals preprocessed with the GLM and physiological nuisance regressors. We show that the GLM-based approach can provide better single trial estimates of brain activity as well as a new feature type, i.e., the weight of the individual and channel-specific hemodynamic response function (HRF) regressor. The improved estimates yield features with higher separability, that significantly enhance accuracy in a binary classification task when compared to conventional preprocessing—on average +7.4% across subjects and feature types. We propose to adapt this well-established approach from neuroscience to the domain of single-trial analysis and preprocessing wherever the classification of evoked brain activity is of concern, for instance in BCI.

Erratum: Prolonged monitoring of cerebral blood flow and autoregulation with diffuse correlation spectroscopy in neurocritical care patients

[This corrects the article DOI: 10.1117/1.NPh.5.4.045005.].

Awake chronic mouse model of targeted pial vessel occlusion via photothrombosis

Animal models of stroke are used extensively to study the mechanisms involved in the acute and chronic phases of recovery following stroke. A translatable animal model that closely mimics the mechanisms of a human stroke is essential in understanding recovery processes as well as developing therapies that improve functional outcomes. We describe a photothrombosis stroke model that is capable of targeting a single distal pial branch of the middle cerebral artery with minimal damage to the surrounding parenchyma in awake head-fixed mice. Mice are implanted with chronic cranial windows above one hemisphere of the brain that allow optical access to study recovery mechanisms for over a month following occlusion. Additionally, we study the effect of laser spot size used for occlusion and demonstrate that a spot size with small axial and lateral resolution has the advantage of minimizing unwanted photodamage while still monitoring macroscopic changes to cerebral blood flow during photothrombosis. We show that temporally guiding illumination using real-time feedback of blood flow dynamics also minimized unwanted photodamage to the vascular network. Finally, through quantifiable behavior deficits and chronic imaging we show that this model can be used to study recovery mechanisms or the effects of therapeutics longitudinally.

Improved physiological noise regression in fNIRS: A multimodal extension of the General Linear Model using temporally embedded Canonical Correlation Analysis

For the robust estimation of evoked brain activity from functional Near-Infrared Spectroscopy (fNIRS) signals, it is crucial to reduce nuisance signals from systemic physiology and motion. The current best practice incorporates short-separation (SS) fNIRS measurements as regressors in a General Linear Model (GLM). However, several challenging signal characteristics such as non-instantaneous and non-constant coupling are not yet addressed by this approach and additional auxiliary signals are not optimally exploited. We have recently introduced a new methodological framework for the unsupervised multivariate analysis of fNIRS signals using Blind Source Separation (BSS) methods. Building onto the framework, in this manuscript we show how to incorporate the advantages of regularized temporally embedded Canonical Correlation Analysis (tCCA) into the supervised GLM. This approach allows flexible integration of any number of auxiliary modalities and signals. We provide guidance for the selection of optimal parameters and auxiliary signals for the proposed GLM extension. Its performance in the recovery of evoked HRFs is then evaluated using both simulated ground truth data and real experimental data and compared with the GLM with short-separation regression. Our results show that the GLM with tCCA significantly improves upon the current best practice, yielding significantly better results across all applied metrics: Correlation (HbO max. +45%), Root Mean Squared Error (HbO max. −55%), F-Score (HbO up to 3.25-fold) and p-value as well as power spectral density of the noise floor. The proposed method can be incorporated into the GLM in an easily applicable way that flexibly combines any available auxiliary signals into optimal nuisance regressors. This work has potential significance both for conventional neuroscientific fNIRS experiments as well as for emerging applications of fNIRS in everyday environments, medicine and BCI, where high Contrast to Noise Ratio is of importance for single trial analysis.

Portable System for Time-Domain Diffuse Correlation Spectroscopy

We introduce a portable system for clinical studies based on time-domain diffuse correlation spectroscopy (DCS). After evaluating different lasers and detectors, the final system is based on a pulsed laser with about 550 ps pulsewidth, a coherence length of 38 mm, and two types of single-photon avalanche diodes (SPAD). The higher efficiency of the red-enhanced SPAD maximizes detection of the collected light, increasing the signal-to-noise ratio, while the better timing response of the CMOS SPAD optimizes the selection of late photons and increases spatial resolution. We discuss component selection and performance, and we present a full characterization of the system, measurement stability, a phantom-based validation study, and preliminary in vivo results collected from the forearms and the foreheads of four healthy subjects. With this system, we are able to resolve blood flow changes 1 cm below the skin surface with improved depth sensitivity and spatial resolution with respect to continuous wave DCS.

Development of a beam propagation method to simulate the point spread function degradation in scattering media

Scattering is one of the main issues that limit the imaging depth in deep tissue optical imaging. To characterize the role of scattering, we have developed a forward model based on the beam propagation method and established the link between the macroscopic optical properties of the media and the statistical parameters of the phase masks applied to the wavefront. Using this model, we have analyzed the degradation of the point-spread function of the illumination beam in the transition regime from ballistic to diffusive light transport. Our method provides a wave-optic simulation toolkit to analyze the effects of scattering on image quality degradation in scanning microscopy. Our open-source implementation is available at https://github.com/BUNPC/Beam-Propagation-Method.

Contribution of speckle noise in near-infrared spectroscopy measurements

Near-infrared spectroscopy (NIRS) is widely used in biomedical optics with applications ranging from basic science, such as in functional neuroimaging, to clinical, as in pulse oximetry. Despite the relatively low absorption of tissue in the near-infrared, there is still a significant amount of optical attenuation produced by the highly scattering nature of tissue. Because of this, designers of NIRS systems have to balance source optical power and source–detector separation to maximize the signal-to-noise ratio (SNR). However, theoretical estimations of SNR neglect the effects of speckle. Speckle manifests as fluctuations of the optical power received at the detector. These fluctuations are caused by interference of the multiple random paths taken by photons in tissue. We present a model for the NIRS SNR that includes the effects of speckle. We performed experimental validations with a NIRS system to show that it agrees with our model. Additionally, we performed computer simulations based on the model to estimate the contribution of speckle noise for different collection areas and source–detector separations. We show that at short source–detector separation, speckle contributes most of the noise when using long coherence length sources. Considering this additional noise is especially important for hybrid applications that use NIRS and speckle contrast simultaneously, such as in diffuse correlation spectroscopy.

Efficient non-degenerate two-photon excitation for fluorescence microscopy

Non-degenerate two-photon excitation (ND-TPE) has been explored in two-photon excitation microscopy. However, a systematic study of the efficiency of ND-TPE to guide the selection of fluorophore excitation wavelengths is missing. We measured the relative non-degenerate two-photon absorption cross-section (ND-TPACS) of several commonly used fluorophores (two fluorescent proteins and three small-molecule dyes) and generated 2-dimensional ND-TPACS spectra. We observed that the shape of a ND-TPACS spectrum follows that of the corresponding degenerate two-photon absorption cross-section (D-TPACS) spectrum, but is higher in magnitude. We found that the observed enhancements are higher than theoretical predictions.

More homogeneous capillary flow and oxygenation in deeper cortical layers correlate with increased oxygen extraction

Our understanding of how capillary blood flow and oxygen distribute across cortical layers to meet the local metabolic demand is incomplete. We addressed this question by using two-photon imaging of resting-state microvascular oxygen partial pressure (PO₂) and flow in the whisker barrel cortex in awake mice. Our measurements in layers I-V show that the capillary red-blood-cell flux and oxygenation heterogeneity, and the intracapillary resistance to oxygen delivery, all decrease with depth, reaching a minimum around layer IV, while the depth-dependent oxygen extraction fraction is increased in layer IV, where oxygen demand is presumably the highest. Our findings suggest that more homogeneous distribution of the physiological observables relevant to oxygen transport to tissue is an important part of the microvascular network adaptation to local brain metabolism. These results will inform the biophysical models of layer-specific cerebral oxygen delivery and consumption and improve our understanding of the diseases that affect cerebral microcirculation.

Spatio-temporal dynamics of cerebral capillary segments with stalling red blood cells

Optical coherence tomography (OCT) allows label-free imaging of red blood cell (RBC) flux within capillaries with high spatio-temporal resolution. In this study, we utilized time-series OCT-angiography to demonstrate interruptions in capillary RBC flux in mouse brain in vivo. We noticed ∼7.5% of ∼200 capillaries had at least one stall in awake mice with chronic windows during a 9-min recording. At any instant, ∼0.45% of capillaries were stalled. Average stall duration was ∼15 s but could last over 1 min. Stalls were more frequent and longer lasting in acute window preparations. Further, isoflurane anesthesia in chronic preparations caused an increase in the number of stalls. In repeated imaging, the same segments had a tendency to stall again over a period of one month. In awake animals, functional stimulation decreased the observance of stalling events. Stalling segments were located distally, away from the first couple of arteriolar-side capillary branches and their average RBC and plasma velocities were lower than nonstalling capillaries within the same region. This first systematic analysis of capillary RBC stalls in the brain, enabled by rapid and continuous volumetric imaging of capillaries with OCT-angiography, will lead to future investigations of the potential role of stalling events in cerebral pathologies.

Non-invasive methods for measuring vascular changes in neurovascular headaches

Vascular changes during spontaneous headache attacks have been studied over the last 30 years. The interest in cerebral vessels in headache research was initially due to the hypothesis of cerebral vessels as the pain source. Here, we review the knowledge gained by measuring the cerebral vasculature during spontaneous primary headache attacks with the use of single photon emission tomography (SPECT), positron emission tomography (PET), magnetic resonance imaging (MRA) and transcranial Doppler (TCD). Furthermore, the use of near-infrared spectroscopy in headache research is reviewed. Existing TCD studies of migraine and other headache disorders do not provide solid evidence for cerebral blood flow velocity changes during spontaneous attacks of migraine headache. SPECT studies have clearly shown cortical vascular changes following migraine aura and the differences between migraine with aura compared to migraine without aura. PET studies have shown focal activation in brain structures related to headache, but whether the changes are specific to different primary headaches have yet to be demonstrated. MR angiography has shown precise changes in large cerebral vessels during spontaneous migraine without aura attacks. Future development in more precise imaging methods may further elucidate the pathophysiological mechanisms in primary headaches.

Normalized field autocorrelation function-based optical coherence tomography three-dimensional angiography

Optical coherence tomography angiography (OCTA) has been widely used for en face visualization of the microvasculature, but is challenged for real three-dimensional (3-D) topologic imaging due to the "tail" artifacts that appear below large vessels. Further, OCTA is generally incapable of differentiating descending arterioles from ascending venules. We introduce a normalized field autocorrelation function-based OCTA (g1-OCTA), which minimizes the tail artifacts and is capable of distinguishing penetrating arterioles from venules in the 3-D image. g1   (  τ  )   is calculated from repeated optical coherence tomography (OCT) acquisitions for each spatial location. The decay amplitude of g1   (  τ  )   is retrieved to represent the dynamics for each voxel. To account for the small g1   (  τ  )   decay in capillaries where red blood cells are flowing slowly and discontinuously, Intralipid is injected to enhance the OCT signal. We demonstrate that the proposed technique realizes 3-D OCTA with negligible tail projections and the penetrating arteries are readily identified. In addition, compared to regular OCTA, the proposed g1-OCTA largely increased the depth-of-field. This technique provides a more accurate rendering of the vascular 3-D anatomy and has the potential for more quantitative characterization of vascular networks.

Choosing a laser for laser speckle contrast imaging

The use of laser speckle contrast imaging (LSCI) has expanded rapidly for characterizing the motion of scattering particles. Speckle contrast is related to the dynamics of the scattering particles via a temporal autocorrelation function, but the quality of various elements of the imaging system can adversely affect the quality of the signal recorded by LSCI. While it is known that the laser coherence affects the speckle contrast, it is generally neglected in in vivo LSCI studies and was not thoroughly addressed in a practical matter. In this work, we address the question of how the spectral width of the light source affects the speckle contrast both experimentally and through numerical simulations. We show that commonly used semiconductor laser diodes have a larger than desired spectral width that results in a significantly reduced speckle contrast compared with ideal narrow band lasers. This results in a reduced signal-to-noise ratio for estimating changes in the motion of scattering particles. We suggest using a volume holographic grating stabilized laser diode or other diodes that have a spectrum of emitted light narrower than ≈1 nm to improve the speckle contrast.

Dependence of the MR signal on the magnetic susceptibility of blood studied with models based on real microvascular networks

PURPOSE

The primary goal of this study was to estimate the value of β , the exponent in the power law relating changes of the transverse relaxation rate and intra-extravascular local magnetic susceptibility differences as Δ R 2 ∗ ∝ ( Δ χ ) β . The secondary objective was to evaluate any differences that might exist in the value of β obtained using a deoxyhemoglobin-weighted Δ χ distribution versus a constant Δ χ distribution assumed in earlier computations. The third objective was to estimate the value of β that is relevant for methods based on susceptibility contrast agents with a concentration of Δ χ higher than that used for BOLD fMRI calculations.

METHODS

Our recently developed model of real microvascular anatomical networks is used to extend the original simplified Monte-Carlo simulations to compute β from the first principles.

RESULTS

Our results show that β = 1 for most BOLD fMRI measurements of real vascular networks, as opposed to earlier predictions of β = 1 .5 using uniform Δ χ distributions. For perfusion or fMRI methods based on contrast agents, which generate larger values for Δ χ , β = 1 for B 0 ≤ 9.4 T, whereas at 14 T β can drop below 1 and the variation across subjects is large, indicating that a lower concentration of contrast agent with a lower value of Δ χ is desired for experiments at high B₀ .

CONCLUSION

These results improve our understanding of the relationship between R₂ ^* and the underlying microvascular properties. The findings will help to infer the cerebral metabolic rate of oxygen and cerebral blood volume from BOLD and perfusion MRI, respectively.

Awake Mouse Imaging: From Two-Photon Microscopy to Blood Oxygen Level-Dependent Functional Magnetic Resonance Imaging

BACKGROUND

Functional magnetic resonance imaging (fMRI) in awake behaving mice is well positioned to bridge the detailed cellular-level view of brain activity, which has become available owing to recent advances in microscopic optical imaging and genetics, to the macroscopic scale of human noninvasive observables. However, though microscopic (e.g., two-photon imaging) studies in behaving mice have become a reality in many laboratories, awake mouse fMRI remains a challenge. Owing to variability in behavior among animals, performing all types of measurements within the same subject is highly desirable and can lead to higher scientific rigor.

METHODS

We demonstrated blood oxygenation level-dependent fMRI in awake mice implanted with long-term cranial windows that allowed optical access for microscopic imaging modalities and optogenetic stimulation. We started with two-photon imaging of single-vessel diameter changes (n = 1). Next, we implemented intrinsic optical imaging of blood oxygenation and flow combined with laser speckle imaging of blood flow obtaining a mesoscopic picture of the hemodynamic response (n = 16). Then we obtained corresponding blood oxygenation level-dependent fMRI data (n = 5). All measurements could be performed in the same mice in response to identical sensory and optogenetic stimuli.

RESULTS

The cranial window did not deteriorate the quality of fMRI and allowed alternation between imaging modalities in each subject.

CONCLUSIONS

This report provides a proof of feasibility for multiscale imaging approaches in awake mice. In the future, this protocol could be extended to include complex cognitive behaviors translatable to humans, such as sensory discrimination or attention.

Two-photon phosphorescence lifetime microscopy of retinal capillary plexus oxygenation in mice

Impaired oxygen delivery and/or consumption in the retinal tissue underlies the pathophysiology of many retinal diseases. However, the essential tools for measuring oxygen concentration in retinal capillaries and studying oxygen transport to retinal tissue are still lacking. We show that two-photon phosphorescence lifetime microscopy can be used to map absolute partial pressures of oxygen (pO2) in the retinal capillary plexus. Measurements were performed at various retinal depths in anesthetized mice under systemic normoxic and hyperoxic conditions. We used a newly developed two-photon phosphorescent oxygen probe, based on a two-photon absorbing platinum tetraphthalimidoporphyrin, and commercially available optics without correction for optical aberrations of the eye. The transverse and axial distances within the tissue volume were calibrated using a model of the eye's optical system. We believe this is the first demonstration of in vivo depth-resolved imaging of pO2 in retinal capillaries. Application of this method has the potential to advance our understanding of oxygen delivery on the microvascular scale and help elucidate mechanisms underlying various retinal diseases.

High-density microfibers as a potential optical interface to reach deep brain regions

OBJECTIVE

Optical techniques for recording and manipulating neural activity have traditionally been constrained to superficial brain regions due to light scattering. New techniques are needed to extend optical access to large 3D volumes in deep brain areas, while retaining local connectivity.

APPROACH

We have developed a method to implant bundles of hundreds or thousands of optical microfibers, each with a diameter of 8 μm. During insertion, each fiber moves independently, following a path of least resistance. The fibers achieve near total internal reflection, enabling optically interfacing with the tissue near each fiber aperture.

MAIN RESULTS

At a depth of 3 mm, histology shows fibers consistently splay over 1 mm in diameter throughout the target region. Immunohistochemical staining after chronic implants reveals neurons in close proximity to the fiber tips. Models of photon fluence indicate that fibers can be used as a stimulation light source to precisely activate distinct patterns of neurons by illuminating a subset of fibers in the bundle. By recording fluorescent beads diffusing in water, we demonstrate the recording capability of the fibers.

SIGNIFICANCE

Our histology, modeling and fluorescent bead recordings suggest that the optical microfibers may provide a minimally invasive, stable, bidirectional interface for recording or stimulating genetic probes in deep brain regions-a hyper-localized form of fiber photometry.

Cardiac pulsatility mapping and vessel type identification using laser speckle contrast imaging

Systemic flow variations caused by the cardiac cycle can play a role or be an important marker in both normal and pathological conditions. The shape, magnitude and propagation speed of the flow pulse reflect mechanical properties of the vasculature and are known to vary significantly with vascular diseases. Most conventional techniques are not capable of imaging cardiac activity in the microcirculation due to spatial and/or temporal resolution limitations and instead make inferences about propagation speed by making measurements at two points along an artery. Here, we apply laser speckle contrast imaging to images with high spatial resolution in the high frequency harmonics of cardiac activity in the cerebral cortex of a mouse. We reveal vessel dependent variation in the cardiac pulse activity and use this information to automatically identify arteries and veins.

Morphine Attenuates fNIRS Signal Associated With Painful Stimuli in the Medial Frontopolar Cortex (medial BA 10)

Functional near infrared spectroscopy (fNIRS) is a non-invasive optical imaging method that provides continuous measure of cortical brain functions. One application has been its use in the evaluation of pain. Previous studies have delineated a deoxygenation process associated with pain in the medial anterior prefrontal region, more specifically, the medial Brodmann Area 10 (BA 10). Such response to painful stimuli has been consistently observed in awake, sedated and anesthetized patients. In this study, we administered oral morphine (15 mg) or placebo to 14 healthy male volunteers with no history of pain or opioid abuse in a crossover double blind design, and performed fNIRS scans prior to and after the administration to assess the effect of morphine on the medial BA 10 pain signal. Morphine is the gold standard for inhibiting nociceptive processing, most well described for brain effects on sensory and emotional regions including the insula, the somatosensory cortex (the primary somatosensory cortex, S1, and the secondary somatosensory cortex, S2), and the anterior cingulate cortex (ACC). Our results showed an attenuation effect of morphine on the fNIRS-measured pain signal in the medial BA 10, as well as in the contralateral S1 (although observed in a smaller number of subjects). Notably, the extent of signal attenuation corresponded with the temporal profile of the reported plasma concentration for the drug. No clear attenuation by morphine on the medial BA 10 response to innocuous stimuli was observed. These results provide further evidence for the role of medial BA 10 in the processing of pain, and also suggest that fNIRS may be used as an objective measure of drug-brain profiles independent of subjective reports.

Prolonged monitoring of cerebral blood flow and autoregulation with diffuse correlation spectroscopy in neurocritical care patients

Monitoring of cerebral blood flow (CBF) and autoregulation are essential components of neurocritical care, but continuous noninvasive methods for CBF monitoring are lacking. Diffuse correlation spectroscopy (DCS) is a noninvasive diffuse optical modality that measures a CBF index ( CBF i ) in the cortex microvasculature by monitoring the rapid fluctuations of near-infrared light diffusing through moving red blood cells. We tested the feasibility of monitoring CBF i with DCS in at-risk patients in the Neurosciences Intensive Care Unit. DCS data were acquired continuously for up to 20 h in six patients with aneurysmal subarachnoid hemorrhage, as permitted by clinical care. Mean arterial blood pressure was recorded synchronously, allowing us to derive autoregulation curves and to compute an autoregulation index. The autoregulation curves suggest disrupted cerebral autoregulation in most patients, with the severity of disruption and the limits of preserved autoregulation varying between subjects. Our findings suggest the potential of the DCS modality for noninvasive, long-term monitoring of cerebral perfusion, and autoregulation.

Extracting individual neural activity recorded through splayed optical microfibers

Previously introduced bundles of hundreds or thousands of microfibers have the potential to extend optical access to deep brain regions, sampling fluorescence activity throughout a three-dimensional volume. Each fiber has a small diameter ( 8    μ m ) and follows a path of least resistance, splaying during insertion. By superimposing the fiber sensitivity profile for each fiber, we model the interface properties for a simulated neural population. Our modeling results suggest that for small ( < 200 ) bundles of fibers, each fiber will collect fluorescence from a small number of nonoverlapping neurons near the fiber apertures. As the number of fibers increases, the bundle delivers more uniform excitation power to the region, moving to a regime where fibers collect fluorescence from more neurons and there is greater overlap between neighboring fibers. Under these conditions, it becomes feasible to apply source separation to extract individual neural contributions. In addition, we demonstrate a source separation technique particularly suited to the interface. Our modeling helps establish performance expectations for this interface and provides a framework for estimating neural contributions under a range of conditions.

In vivo imaging and analysis of cerebrovascular hemodynamic responses and tissue oxygenation in the mouse brain

Cerebrovascular dysfunction has an important role in the pathogenesis of multiple brain disorders. Measurement of hemodynamic responses in vivo can be challenging, particularly as techniques are often not described in sufficient detail and vary between laboratories. We present a set of standardized in vivo protocols that describe high-resolution two-photon microscopy and intrinsic optical signal (IOS) imaging to evaluate capillary and arteriolar responses to a stimulus, regional hemodynamic responses, and oxygen delivery to the brain. The protocol also describes how to measure intrinsic NADH fluorescence to understand how blood O₂ supply meets the metabolic demands of activated brain tissue, and to perform resting-state absolute oxygen partial pressure (pO₂) measurements of brain tissue. These methods can detect cerebrovascular changes at far higher resolution than MRI techniques, although the optical nature of these techniques limits their achievable imaging depths. Each individual procedure requires 1-2 h to complete, with two to three procedures typically performed per animal at a time. These protocols are broadly applicable in studies of cerebrovascular function in healthy and diseased brain in any of the existing mouse models of neurological and vascular disorders. All these procedures can be accomplished by a competent graduate student or experienced technician, except the two-photon measurement of absolute pO₂ level, which is better suited to a more experienced, postdoctoral-level researcher.

Phasor analysis of NADH FLIM identifies pharmacological disruptions to mitochondrial metabolic processes in the rodent cerebral cortex

Investigating cerebral metabolism in vivo at a microscopic level is essential for understanding brain function and its pathological alterations. The intricate signaling and metabolic dynamics between neurons, glia, and microvasculature requires much more detailed understanding to better comprehend the mechanisms governing brain function and its disease-related changes. We recently demonstrated that pharmacologically-induced alterations to different steps of cerebral metabolism can be distinguished utilizing 2-photon fluorescence lifetime imaging of endogenous reduced nicotinamide adenine dinucleotide (NADH) fluorescence in vivo. Here, we evaluate the ability of the phasor analysis method to identify these pharmacological metabolic alterations and compare the method's performance with more conventional nonlinear curve-fitting analysis. Visualization of phasor data, both at the fundamental laser repetition frequency and its second harmonic, enables resolution of pharmacologically-induced alterations to mitochondrial metabolic processes from baseline cerebral metabolism. Compared to our previous classification models based on nonlinear curve-fitting, phasor-based models required fewer parameters and yielded comparable or improved classification accuracy. Fluorescence lifetime imaging of NADH and phasor analysis shows utility for detecting metabolic alterations and will lead to a deeper understanding of cerebral energetics and its pathological changes.