Baseline oxygen consumption decreases with cortical depth.

The cerebral cortex is organized in cortical layers that differ in their cellular density, composition, and wiring. Cortical laminar architecture is also readily revealed by staining for cytochrome oxidase-the last enzyme in the respiratory electron transport chain located in the inner mitochondrial membrane. It has been hypothesized that a high-density band of cytochrome oxidase in cortical layer IV reflects higher oxygen consumption under baseline (unstimulated) conditions. Here, we tested the above hypothesis using direct measurements of the partial pressure of O2 (pO2) in cortical tissue by means of 2-photon phosphorescence lifetime microscopy (2PLM). We revisited our previously developed method for extraction of the cerebral metabolic rate of O2 (CMRO2) based on 2-photon pO2 measurements around diving arterioles and applied this method to estimate baseline CMRO2 in awake mice across cortical layers. To our surprise, our results revealed a decrease in baseline CMRO2 from layer I to layer IV. This decrease of CMRO2 with cortical depth was paralleled by an increase in tissue oxygenation. Higher baseline oxygenation and cytochrome density in layer IV may serve as an O2 reserve during surges of neuronal activity or certain metabolically active brain states rather than reflecting baseline energy needs. Our study provides to our knowledge the first quantification of microscopically resolved CMRO2 across cortical layers as a step towards better understanding of brain energy metabolism.

Scalable Thousand Channel Penetrating Microneedle Arrays on Flex for Multimodal and Large Area Coverage BrainMachine Interfaces.

The Utah array powers cutting-edge projects for restoration of neurological function, such as BrainGate, but the underlying electrode technology has itself advanced little in the last three decades. Here, advanced dual-side lithographic microfabrication processes is exploited to demonstrate a 1024-channel penetrating silicon microneedle array (SiMNA) that is scalable in its recording capabilities and cortical coverage and is suitable for clinical translation. The SiMNA is the first penetrating microneedle array with a flexible backing that affords compliancy to brain movements. In addition, the SiMNA is optically transparent permitting simultaneous optical and electrophysiological interrogation of neuronal activity. The SiMNA is used to demonstrate reliable recordings of spontaneous and evoked field potentials and of single unit activity in chronically implanted mice for up to 196 days in response to optogenetic and to whisker air-puff stimuli. Significantly, the 1024-channel SiMNA establishes detailed spatiotemporal mapping of broadband brain activity in rats. This novel scalable and biocompatible SiMNA with its multimodal capability and sensitivity to broadband brain activity will accelerate the progress in fundamental neurophysiological investigations and establishes a new milestone for penetrating and large area coverage microelectrode arrays for brain-machine interfaces.

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.

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.

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.

Correlation Structure in Micro-ECoG Recordings is Described by Spatially Coherent Components.

Electrocorticography (ECoG) is becoming more prevalent due to improvements in fabrication and recording technology as well as its ease of implantation compared to intracortical electrophysiology, larger cortical coverage, and potential advantages for use in long term chronic implantation. Given the flexibility in the design of ECoG grids, which is only increasing, it remains an open question what geometry of the electrodes is optimal for an application. Conductive polymer, PEDOT:PSS, coated microelectrodes have an advantage that they can be made very small without losing low impedance. This makes them suitable for evaluating the required granularity of ECoG recording in humans and experimental animals. We used two-dimensional (2D) micro-ECoG grids to record intra-operatively in humans and during acute implantations in mouse with separation distance between neighboring electrodes (i.e., pitch) of 0.4 mm and 0.2/0.25 mm respectively. To assess the spatial properties of the signals, we used the average correlation between electrodes as a function of the pitch. In agreement with prior studies, we find a strong frequency dependence in the spatial scale of correlation. By applying independent component analysis (ICA), we find that the spatial pattern of correlation is largely due to contributions from multiple spatially extended, time-locked sources present at any given time. Our analysis indicates the presence of spatially structured activity down to the sub-millimeter spatial scale in ECoG despite the effects of volume conduction, justifying the use of dense micro-ECoG grids.

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.

Interstitial solute transport in 3D reconstructed neuropil occurs by diffusion rather than bulk flow.

The brain lacks lymph vessels and must rely on other mechanisms for clearance of waste products, including amyloid [Formula: see text] that may form pathological aggregates if not effectively cleared. It has been proposed that flow of interstitial fluid through the brain's interstitial space provides a mechanism for waste clearance. Here we compute the permeability and simulate pressure-mediated bulk flow through 3D electron microscope (EM) reconstructions of interstitial space. The space was divided into sheets (i.e., space between two parallel membranes) and tunnels (where three or more membranes meet). Simulation results indicate that even for larger extracellular volume fractions than what is reported for sleep and for geometries with a high tunnel volume fraction, the permeability was too low to allow for any substantial bulk flow at physiological hydrostatic pressure gradients. For two different geometries with the same extracellular volume fraction the geometry with the most tunnel volume had [Formula: see text] higher permeability, but the bulk flow was still insignificant. These simulation results suggest that even large molecule solutes would be more easily cleared from the brain interstitium by diffusion than by bulk flow. Thus, diffusion within the interstitial space combined with advection along vessels is likely to substitute for the lymphatic drainage system in other organs.

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.

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 ΔR2∗∝(Δχ)ß . 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 B0≤ 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 B0 . CONCLUSION: These results improve our understanding of the relationship between R2* 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.

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.

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.

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.

Leveraging Genomic Annotations and Pleiotropic Enrichment for Improved Replication Rates in Schizophrenia GWAS.

Most of the genetic architecture of schizophrenia (SCZ) has not yet been identified. Here, we apply a novel statistical algorithm called Covariate-Modulated Mixture Modeling (CM3), which incorporates auxiliary information (heterozygosity, total linkage disequilibrium, genomic annotations, pleiotropy) for each single nucleotide polymorphism (SNP) to enable more accurate estimation of replication probabilities, conditional on the observed test statistic ("z-score") of the SNP. We use a multiple logistic regression on z-scores to combine information from auxiliary information to derive a "relative enrichment score" for each SNP. For each stratum of these relative enrichment scores, we obtain nonparametric estimates of posterior expected test statistics and replication probabilities as a function of discovery z-scores, using a resampling-based approach that repeatedly and randomly partitions meta-analysis sub-studies into training and replication samples. We fit a scale mixture of two Gaussians model to each stratum, obtaining parameter estimates that minimize the sum of squared differences of the scale-mixture model with the stratified nonparametric estimates. We apply this approach to the recent genome-wide association study (GWAS) of SCZ (n = 82,315), obtaining a good fit between the model-based and observed effect sizes and replication probabilities. We observed that SNPs with low enrichment scores replicate with a lower probability than SNPs with high enrichment scores even when both they are genome-wide significant (p < 5x10-8). There were 693 and 219 independent loci with model-based replication rates ≥80% and ≥90%, respectively. Compared to analyses not incorporating relative enrichment scores, CM3 increased out-of-sample yield for SNPs that replicate at a given rate. This demonstrates that replication probabilities can be more accurately estimated using prior enrichment information with CM3.

Quantifying the microvascular origin of BOLD-fMRI from first principles with two-photon microscopy and an oxygen-sensitive nanoprobe.

The blood oxygenation level-dependent (BOLD) contrast is widely used in functional magnetic resonance imaging (fMRI) studies aimed at investigating neuronal activity. However, the BOLD signal reflects changes in blood volume and oxygenation rather than neuronal activity per se. Therefore, understanding the transformation of microscopic vascular behavior into macroscopic BOLD signals is at the foundation of physiologically informed noninvasive neuroimaging. Here, we use oxygen-sensitive two-photon microscopy to measure the BOLD-relevant microvascular physiology occurring within a typical rodent fMRI voxel and predict the BOLD signal from first principles using those measurements. The predictive power of the approach is illustrated by quantifying variations in the BOLD signal induced by the morphological folding of the human cortex. This framework is then used to quantify the contribution of individual vascular compartments and other factors to the BOLD signal for different magnet strengths and pulse sequences.

Non-degenerate 2-photon excitation in scattering medium for fluorescence microscopy.

Non-degenerate 2-photon excitation (ND-2PE) of a fluorophore with two laser beams of different photon energies offers an independent degree of freedom in tuning of the photon flux for each beam. This feature takes advantage of the infrared wavelengths used in degenerate 3-photon excitation (D-3PE) microscopy to achieve increased penetration depths, while preserving a relatively high 2-photon excitation cross section in comparison to that of D-3PE. Here, using spatially and temporally aligned Ti:Sapphire laser and optical parametric oscillator beams operating at near infrared (NIR) and short-wavelength infrared (SWIR) optical frequencies, we employ ND-2PE and provide a practical demonstration that a constant fluorophore emission intensity is achievable deeper into a scattering medium using ND-2PE as compared to the commonly used degenerate 2-photon excitation (D-2PE).

Robust and fragile aspects of cortical blood flow in relation to the underlying angioarchitecture.

We review the organizational principles of the cortical vasculature and the underlying patterns of blood flow under normal conditions and in response to occlusion of single vessels. The cortex is sourced by a two-dimensional network of pial arterioles that feeds a three-dimensional network of subsurface microvessels in close proximity to neurons and glia. Blood flow within the surface and subsurface networks is largely insensitive to occlusion of a single vessel within either network. However, the penetrating arterioles that connect the pial network to the subsurface network are bottlenecks to flow; occlusion of even a single penetrating arteriole results in the death of a 500 µm diameter cylinder of cortical tissue despite the potential for collateral flow through microvessels. This pattern of flow is consistent with that calculated from a full reconstruction of the angioarchitecture. Conceptually, collateral flow is insufficient to compensate for the occlusion of a penetrating arteriole because penetrating venules act as shunts of blood that flows through collaterals. Future directions that stem from the analysis of the angioarchitecture concern cellular-level issues, in particular the regulation of blood flow within the subsurface microvascular network, and system-level issues, in particular the role of penetrating arteriole occlusions in human cognitive impairment.

In vivo stimulus-induced vasodilation occurs without IP3 receptor activation and may precede astrocytic calcium increase.

Calcium-dependent release of vasoactive gliotransmitters is widely assumed to trigger vasodilation associated with rapid increases in neuronal activity. Inconsistent with this hypothesis, intact stimulus-induced vasodilation was observed in inositol 1,4,5-triphosphate (IP3) type-2 receptor (R2) knock-out (KO) mice, in which the primary mechanism of astrocytic calcium increase-the release of calcium from intracellular stores following activation of an IP3-dependent pathway-is lacking. Further, our results in wild-type (WT) mice indicate that in vivo onset of astrocytic calcium increase in response to sensory stimulus could be considerably delayed relative to the simultaneously measured onset of arteriolar dilation. Delayed calcium increases in WT mice were observed in both astrocytic cell bodies and perivascular endfeet. Thus, astrocytes may not play a role in the initiation of blood flow response, at least not via calcium-dependent mechanisms. Moreover, an increase in astrocytic intracellular calcium was not required for normal vasodilation in the IP3R2-KO animals.

Special Section Guest Editorial: Frontiers in Neurophotonics.

The editorial presents the two-part Special Section on Frontiers in Neurophotonics.

Time-Lagged Functional Ultrasound for Multi-Parametric Cerebral Hemodynamic Imaging.

We introduce an ultrasound speckle decorrelation-based time-lagged functional ultrasound technique (tl-fUS) for the quantification of the relative changes in cerebral blood flow speed (rCBF [Formula: see text]), cerebral blood volume (rCBV) and cerebral blood flow (rCBF) during functional stimulations. Numerical simulations, phantom validations, and in vivo mouse brain experiments were performed to test the capability of tl-fUS to parse out and quantify the ratio change of these hemodynamic parameters. The blood volume change was found to be more prominent in arterioles compared to venules and the peak blood flow changes were around 2.5 times the peak blood volume change during brain activation, agreeing with previous observations in the literature. The tl-fUS shows the ability of distinguishing the relative changes of rCBFspeed, rCBV, and rCBF, which can inform specific physiological interpretations of the fUS measurements.