Microscale Physiological Events on the Human Cortical Surface

Despite ongoing advances in our understanding of local single-cellular and network-level activity of neuronal populations in the human brain, extraordinarily little is known about their "intermediate" microscale local circuit dynamics. Here, we utilized ultra-high-density microelectrode arrays and a rare opportunity to perform intracranial recordings across multiple cortical areas in human participants to discover three distinct classes of cortical activity that are not locked to ongoing natural brain rhythmic activity. The first included fast waveforms similar to extracellular single-unit activity. The other two types were discrete events with slower waveform dynamics and were found preferentially in upper cortical layers. These second and third types were also observed in rodents, nonhuman primates, and semi-chronic recordings from humans via laminar and Utah array microelectrodes. The rates of all three events were selectively modulated by auditory and electrical stimuli, pharmacological manipulation, and cold saline application and had small causal co-occurrences. These results suggest that the proper combination of high-resolution microelectrodes and analytic techniques can capture neuronal dynamics that lay between somatic action potentials and aggregate population activity. Understanding intermediate microscale dynamics in relation to single-cell and network dynamics may reveal important details about activity in the full cortical circuit.

Diversity is our Strength

Neurophotonics Editor in Chief Anna Devor reflects on need to cherish and cultivate diversity within the global neuroscience and neurophotonics community by creating inclusive environments to welcome young students in sharing the joy of science.

A Suite of Neurophotonic Tools to Underpin the Contribution of Internal Brain States in fMRI

Recent developments in optical microscopy, applicable for large-scale and longitudinal imaging of cortical activity in behaving animals, open unprecedented opportunities to gain a deeper understanding of neurovascular and neurometabolic coupling during different brain states. Future studies will leverage these tools to deliver foundational knowledge about brain state-dependent regulation of cerebral blood flow and metabolism as well as regulation as a function of brain maturation and aging. This knowledge is of critical importance to interpret hemodynamic signals observed with functional magnetic resonance imaging (fMRI).

All-Optical Electrophysiology in hiPSC-Derived Neurons With Synthetic Voltage Sensors

Voltage imaging and “all-optical electrophysiology” in human induced pluripotent stem cell (hiPSC)-derived neurons have opened unprecedented opportunities for high-throughput phenotyping of activity in neurons possessing unique genetic backgrounds of individual patients. While prior all-optical electrophysiology studies relied on genetically encoded voltage indicators, here, we demonstrate an alternative protocol using a synthetic voltage sensor and genetically encoded optogenetic actuator that generate robust and reproducible results. We demonstrate the functionality of this method by measuring spontaneous and evoked activity in three independent hiPSC-derived neuronal cell lines with distinct genetic backgrounds.

Welcome New Members of the Neurophotonics Editorial Board!

Neurophotonics Editor in Chief Anna Devor welcomes an expanded editorial board and announces three Special Sections for the coming year.

Imaging localized fast optical signals of neural activation with optical coherence tomography in awake mice

We report optical coherence tomography (OCT) imaging of localized fast optical signals (FOSs) arising from whisker stimulation in awake mice. The activated voxels were identified by fitting the OCT intensity signal time course with a response function over a time scale of a few hundred milliseconds after the whisker stimulation. The significantly activated voxels were shown to be localized to the expected brain region for whisker stimulation. The ability to detect functional stimulus-evoked, depth-resolved FOS with intrinsic contrast from the cortex provides a new tool for neural activity studies.

Chronic Cranial Windows for Long Term Multimodal Neurovascular Imaging in Mice

Chronic cranial windows allow for longitudinal brain imaging experiments in awake, behaving mice. Different imaging technologies have their unique advantages and combining multiple imaging modalities offers measurements of a wide spectrum of neuronal, glial, vascular, and metabolic parameters needed for comprehensive investigation of physiological and pathophysiological mechanisms. Here, we detail a suite of surgical techniques for installation of different cranial windows targeted for specific imaging technologies and their combination. Following these techniques and practices will yield higher experimental success and reproducibility of results.

Editorial: Neurophotonics for you

Neurophotonics Editor in Chief Anna Devor outlines a vision for the journal.

Impact of Brain Surface Boundary Conditions on Electrophysiology and Implications for Electrocorticography

Volume conduction of electrical potentials in the brain is highly influenced by the material properties and geometry of the tissue and recording devices implanted into the tissue. These effects are very large in EEG due to the volume conduction through the skull and scalp but are often neglected in intracranial electrophysiology. When considering penetrating electrodes deep in the brain, the assumption of an infinite and homogenous medium can be used when the sources are far enough from the brain surface and the electrodes to minimize the boundary effect. When the electrodes are recording from the brain's surface the effect of the boundary cannot be neglected, and the large surface area and commonly used insulating materials in surface electrode arrays may further increase the effect by altering the nature of the boundary in the immediate vicinity of the electrodes. This gives the experimenter some control over the spatial profiles of the potentials by appropriate design of the electrode arrays. We construct a simple three-layer model to describe the effect of material properties and geometry above the brain surface on the electric potentials and conduct empirical experiments to validate this model. A laminar electrode array is used to measure the effect of insulating and relatively conducting layers above the cortical surface by recording evoked potentials alternating between a dried surface and saline covering layer, respectively. Empirically, we find that an insulating boundary amplifies the potentials relative to conductive saline by about a factor of 4, and that the effect is not constrained to potentials that originate near the surface. The model is applied to predict the influence of array design and implantation procedure on the recording amplitude and spatial selectivity of the surface electrode arrays.

Understanding the genetic determinants of the brain with MOSTest

Regional brain morphology has a complex genetic architecture, consisting of many common polymorphisms with small individual effects. This has proven challenging for genome-wide association studies (GWAS). Due to the distributed nature of genetic signal across brain regions, multivariate analysis of regional measures may enhance discovery of genetic variants. Current multivariate approaches to GWAS are ill-suited for complex, large-scale data of this kind. Here, we introduce the Multivariate Omnibus Statistical Test (MOSTest), with an efficient computational design enabling rapid and reliable inference, and apply it to 171 regional brain morphology measures from 26,502 UK Biobank participants. At the conventional genome-wide significance threshold of α = 5 × 10-8, MOSTest identifies 347 genomic loci associated with regional brain morphology, more than any previous study, improving upon the discovery of established GWAS approaches more than threefold. Our findings implicate more than 5% of all protein-coding genes and provide evidence for gene sets involved in neuron development and differentiation.

Spatially resolved estimation of metabolic oxygen consumption from optical measurements in cortex

Significance: The cerebral metabolic rate of oxygen (CMRO 2 ) is an important indicator of brain function and pathology. Knowledge about its magnitude is also required for proper interpretation of the blood oxygenation level-dependent (BOLD) signal measured with functional MRI. Despite the need for estimatingCMRO 2 , no gold standard exists. Traditionally, the estimation ofCMRO 2 has been pursued with somewhat indirect approaches combining several different types of measurements with mathematical modeling of the underlying physiological processes. The recent ability to measure the level of oxygen (pO 2 ) in cortex with two-photon resolution in in vivo conditions has provided a more direct way for estimatingCMRO 2 , but has so far only been used to estimate the averageCMRO 2 close to cortical penetrating arterioles in rats. Aim: The aim of this study was to propose a method to provide spatial maps ofCMRO 2 based on two-photonpO 2 measurements. Approach: The method has two key steps. First, thepO 2 maps are spatially smoothed to reduce the effects of noise in the measurements. Next, the Laplace operator (a double spatial derivative) in two spatial dimensions is applied on the smoothedpO 2 maps to obtain spatially resolvedCMRO 2 estimates. Result: The smoothing introduces a bias, and a balance must be found where the effects of the noise are sufficiently reduced without introducing too much bias. In this model-based study, we explored this balance using synthetic model-based data, that is, data where the spatial maps ofCMRO 2 were preset and thus known. The correspondingpO 2 maps were found by solving the Poisson equation, which relatesCMRO 2 andpO 2 . MATLAB code for using the method is provided. Conclusion: Through this model-based study, we propose a method for estimatingCMRO 2 with high spatial resolution based on measurements ofpO 2 in cerebral cortex.

Secondary effects on brain physiology caused by focused ultrasound-mediated disruption of the blood-brain barrier

Focused ultrasound (FUS) combined with microbubbles is a non-invasive method for targeted, reversible disruption of the blood-brain barrier (FUS-BBB opening). This approach holds great promise for improving delivery of therapeutics to the brain. In order to achieve this clinically important goal, the approach necessarily breaks a protective barrier, temporarily, which plays a fundamental role in maintaining a homeostatic environment in the brain. Preclinical and clinical research has identified a set of treatment parameters under which this can be performed safely, whereby the BBB is disrupted to the point of being permeable to normally non-penetrant agents without causing significant acute damage to endothelial or neuronal cells. Much of the early work in this field focused on engineering questions around how to achieve optimal delivery of therapeutics via BBB disruption. However, there is increasing interest in addressing biological questions related to whether and how various aspects of neurophysiology might be affected when this fundamental protective barrier is compromised by the specific mechanisms of FUS-BBB opening. Improving our understanding of these secondary effects is becoming vital now that FUS-BBB opening treatments have entered clinical trials. Such information would help to safely expand FUS-BBB opening protocols into a wider range of drug delivery applications and may even lead to new types of treatments. In this paper, we will critically review our current knowledge of the secondary effects caused by FUS-BBB opening on brain physiology, identify areas that remain understudied, and discuss how a better understanding of these processes can be used to safely advance FUS-BBB opening into a wider range of clinical applications.

Comparing the fundamental imaging depth limit of two-photon, three-photon, and non-degenerate two-photon microscopy

We have systematically characterized the degradation of imaging quality with depth in deep brain multi-photon microscopy, utilizing our recently developed numerical model that computes wave propagation in scattering media. The signal-to-background ratio (SBR) and the resolution determined by the width of the point spread function are obtained as functions of depth. We compare the imaging quality of two-photon (2PM), three-photon (3PM), and non-degenerate two-photon microscopy (ND-2PM) for mouse brain imaging. We show that the imaging depth of 2PM and ND-2PM are fundamentally limited by the SBR, while the SBR remains approximately invariant with imaging depth for 3PM. Instead, the imaging depth of 3PM is limited by the degradation of the resolution, if there is sufficient laser power to maintain the signal level at large depth. The roles of the concentration of dye molecules, the numerical aperture of the input light, the anisotropy factor g, noise level, input laser power, and the effect of temporal broadening are also discussed.

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.

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.

Complex Oscillatory Waves Emerging from Cortical Organoids Model Early Human Brain Network Development

Structural and transcriptional changes during early brain maturation follow fixed developmental programs defined by genetics. However, whether this is true for functional network activity remains unknown, primarily due to experimental inaccessibility of the initial stages of the living human brain. Here, we developed human cortical organoids that dynamically change cellular populations during maturation and exhibited consistent increases in electrical activity over the span of several months. The spontaneous network formation displayed periodic and regular oscillatory events that were dependent on glutamatergic and GABAergic signaling. The oscillatory activity transitioned to more spatiotemporally irregular patterns, and synchronous network events resembled features similar to those observed in preterm human electroencephalography. These results show that the development of structured network activity in a human neocortex model may follow stable genetic programming. Our approach provides opportunities for investigating and manipulating the role of network activity in the developing human cortex.

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.

Computational Modeling of Genetic Contributions to Excitability and Neural Coding in Layer V Pyramidal Cells: Applications to Schizophrenia Pathology

Pyramidal cells in layer V of the neocortex are one of the most widely studied neuron types in the mammalian brain. Due to their role as integrators of feedforward and cortical feedback inputs, they are well-positioned to contribute to the symptoms and pathology in mental disorders-such as schizophrenia-that are characterized by a mismatch between the internal perception and external inputs. In this modeling study, we analyze the input/output properties of layer V pyramidal cells and their sensitivity to modeled genetic variants in schizophrenia-associated genes. We show that the excitability of layer V pyramidal cells and the way they integrate inputs in space and time are altered by many types of variants in ion-channel and Ca²⁺ transporter-encoding genes that have been identified as risk genes by recent genome-wide association studies. We also show that the variability in the output patterns of spiking and Ca²⁺ transients in layer V pyramidal cells is altered by these model variants. Importantly, we show that many of the predicted effects are robust to noise and qualitatively similar across different computational models of layer V pyramidal cells. Our modeling framework reveals several aspects of single-neuron excitability that can be linked to known schizophrenia-related phenotypes and existing hypotheses on disease mechanisms. In particular, our models predict that single-cell steady-state firing rate is positively correlated with the coding capacity of the neuron and negatively correlated with the amplitude of a prepulse-mediated adaptation and sensitivity to coincidence of stimuli in the apical dendrite and the perisomatic region of a layer V pyramidal cell. These results help to uncover the voltage-gated ion-channel and Ca²⁺ transporter-associated genetic underpinnings of schizophrenia phenotypes and biomarkers.

Selective Formation of Porous Pt Nanorods for Highly Electrochemically Efficient Neural Electrode Interfaces

The enhanced electrochemical activity of nanostructured materials is readily exploited in energy devices, but their utility in scalable and human-compatible implantable neural interfaces can significantly advance the performance of clinical and research electrodes. We utilize low-temperature selective dealloying to develop scalable and biocompatible one-dimensional platinum nanorod (PtNR) arrays that exhibit superb electrochemical properties at various length scales, stability, and biocompatibility for high performance neurotechnologies. PtNR arrays record brain activity with cellular resolution from the cortical surfaces in birds and nonhuman primates. Significantly, strong modulation of surface recorded single unit activity by auditory stimuli is demonstrated in European Starling birds as well as the modulation of local field potentials in the visual cortex by light stimuli in a nonhuman primate and responses to electrical stimulation in mice. PtNRs record behaviorally and physiologically relevant neuronal dynamics from the surface of the brain with high spatiotemporal resolution, which paves the way for less invasive brain-machine interfaces.

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.

Measurement of the relative non-degenerate two-photon absorption cross-section for fluorescence microscopy

In non-degenerate two-photon microscopy (ND-TPM), the required energy for fluorescence excitation occurs via absorption of two photons of different energies derived from two synchronized pulsed laser beams. ND-TPM is a promising imaging technology offering flexibility in the choice of the photon energy for each beam. However, a formalism to quantify the efficiency of two-photon absorption (TPA) under non-degenerate excitation, relative to the resonant degenerate excitation, is missing. Here, we derive this formalism and experimentally validate our prediction for a common fluorophore, fluorescein. An accurate quantification of non-degenerate TPA is important to optimize the choice of photon energies for each fluorophore.

Alterations in Schizophrenia-Associated Genes Can Lead to Increased Power in Delta Oscillations

Genome-wide association studies have implicated many ion channels in schizophrenia pathophysiology. Although the functions of these channels are relatively well characterized by single-cell studies, the contributions of common variation in these channels to neurophysiological biomarkers and symptoms of schizophrenia remain elusive. Here, using computational modeling, we show that a common biomarker of schizophrenia, namely, an increase in delta-oscillation power, may be a direct consequence of altered expression or kinetics of voltage-gated ion channels or calcium transporters. Our model of a circuit of layer V pyramidal cells highlights multiple types of schizophrenia-related variants that contribute to altered dynamics in the delta-frequency band. Moreover, our model predicts that the same membrane mechanisms that increase the layer V pyramidal cell network gain and response to delta-frequency oscillations may also cause a deficit in a single-cell correlate of the prepulse inhibition, which is a behavioral biomarker highly associated with schizophrenia.

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.