Single-cell atlas of the human brain vasculature across development, adulthood and disease.

A broad range of brain pathologies critically relies on the vasculature, and cerebrovascular disease is a leading cause of death worldwide. However, the cellular and molecular architecture of the human brain vasculature remains incompletely understood1. Here we performed single-cell RNA sequencing analysis of 606,380 freshly isolated endothelial cells, perivascular cells and other tissue-derived cells from 117 samples, from 68 human fetuses and adult patients to construct a molecular atlas of the developing fetal, adult control and diseased human brain vasculature. We identify extensive molecular heterogeneity of the vasculature of healthy fetal and adult human brains and across five vascular-dependent central nervous system (CNS) pathologies, including brain tumours and brain vascular malformations. We identify alteration of arteriovenous differentiation and reactivated fetal as well as conserved dysregulated genes and pathways in the diseased vasculature. Pathological endothelial cells display a loss of CNS-specific properties and reveal an upregulation of MHC class II molecules, indicating atypical features of CNS endothelial cells. Cell-cell interaction analyses predict substantial endothelial-to-perivascular cell ligand-receptor cross-talk, including immune-related and angiogenic pathways, thereby revealing a central role for the endothelium within brain neurovascular unit signalling networks. Our single-cell brain atlas provides insights into the molecular architecture and heterogeneity of the developing, adult/control and diseased human brain vasculature and serves as a powerful reference for future studies.

Highly local activation of inhibition at the seizure wavefront in vivo.

The propagation of a seizure wavefront in the cortex divides an intensely firing seizure core from a low-firing seizure penumbra. Seizure propagation is currently thought to generate strong activation of inhibition in the seizure penumbra that leads to its decreased neuronal firing. However, the direct measurement of neuronal excitability during seizures has been difficult to perform in vivo. We used simultaneous optogenetics and calcium imaging (all-optical interrogation) to characterize real-time neuronal excitability in an acute mouse model of seizure propagation. We find that single-neuron excitability is decreased in close proximity to the seizure wavefront but becomes increased distal to the seizure wavefront. This suggests that inhibitory neurons of the seizure wavefront create a proximal circumference of hypoexcitability but do not influence neuronal excitability in the penumbra.

Control of working memory by phase-amplitude coupling of human hippocampal neurons.

Retaining information in working memory is a demanding process that relies on cognitive control to protect memoranda-specific persistent activity from interference1,2. However, how cognitive control regulates working memory storage is unclear. Here we show that interactions of frontal control and hippocampal persistent activity are coordinated by theta-gamma phase-amplitude coupling (TG-PAC). We recorded single neurons in the human medial temporal and frontal lobe while patients maintained multiple items in their working memory. In the hippocampus, TG-PAC was indicative of working memory load and quality. We identified cells that selectively spiked during nonlinear interactions of theta phase and gamma amplitude. The spike timing of these PAC neurons was coordinated with frontal theta activity when cognitive control demand was high. By introducing noise correlations with persistently active neurons in the hippocampus, PAC neurons shaped the geometry of the population code. This led to higher-fidelity representations of working memory content that were associated with improved behaviour. Our results support a multicomponent architecture of working memory1,2, with frontal control managing maintenance of working memory content in storage-related areas3-5. Within this framework, hippocampal TG-PAC integrates cognitive control and working memory storage across brain areas, thereby suggesting a potential mechanism for top-down control over sensory-driven processes.

Evaluating the Current State of Epilepsy Care in the Province of Ontario.

There are numerous challenges pertaining to epilepsy care across Ontario, including Epilepsy Monitoring Unit (EMU) bed pressures, surgical access and community supports. We sampled the current clinical, community and operational state of Ontario epilepsy centres and community epilepsy agencies post COVID-19 pandemic. A 44-item survey was distributed to all 11 district and regional adult and paediatric Ontario epilepsy centres. Qualitative responses were collected from community epilepsy agencies. Results revealed ongoing gaps in epilepsy care across Ontario, with EMU bed pressures and labour shortages being limiting factors. A clinical network advising the Ontario Ministry of Health will improve access to epilepsy care.

In-silico testing of new pharmacology for restoring inhibition and human cortical function in depression.

Reduced inhibition by somatostatin-expressing interneurons is associated with depression. Administration of positive allosteric modulators of α5 subunit-containing GABAA receptor (α5-PAM) that selectively target this lost inhibition exhibit antidepressant and pro-cognitive effects in rodent models of chronic stress. However, the functional effects of α5-PAM on the human brain in vivo are unknown, and currently cannot be assessed experimentally. We modeled the effects of α5-PAM on tonic inhibition as measured in human neurons, and tested in silico α5-PAM effects on detailed models of human cortical microcircuits in health and depression. We found that α5-PAM effectively recovered impaired cortical processing as quantified by stimulus detection metrics, and also recovered the power spectral density profile of the microcircuit EEG signals. We performed an α5-PAM dose-response and identified simulated EEG biomarker candidates. Our results serve to de-risk and facilitate α5-PAM translation and provide biomarkers in non-invasive brain signals for monitoring target engagement and drug efficacy.

Diversity-induced trivialization and resilience of neural dynamics.

Heterogeneity is omnipresent across all living systems. Diversity enriches the dynamical repertoire of these systems but remains challenging to reconcile with their manifest robustness and dynamical persistence over time, a fundamental feature called resilience. To better understand the mechanism underlying resilience in neural circuits, we considered a nonlinear network model, extracting the relationship between excitability heterogeneity and resilience. To measure resilience, we quantified the number of stationary states of this network, and how they are affected by various control parameters. We analyzed both analytically and numerically gradient and non-gradient systems modeled as non-linear sparse neural networks evolving over long time scales. Our analysis shows that neuronal heterogeneity quenches the number of stationary states while decreasing the susceptibility to bifurcations: a phenomenon known as trivialization. Heterogeneity was found to implement a homeostatic control mechanism enhancing network resilience to changes in network size and connection probability by quenching the system's dynamic volatility.

Hardware-Efficient 1D CNN for Patient-Specific Early Seizure Detection.

Closed-loop brain-implantable neuromodulation devices are a new treatment option for patients with refractory epilepsy. Seizure detection algorithms implemented on such devices are subject to strict power and area constraints. Deep learning methods, though very powerful, tend to have high computational complexity and thus are typically impractical for resource-constrained neuromodulation devices. In this paper, we propose a compact and hardware-efficient one-dimensional convolutional neural network (1D CNN) structure for patient-specific early seizure detection. Feature extraction techniques and a novel initialization method based on the forward-chaining training and testing scheme are used to improve model performance. Our compact model achieves similar accuracy to that of support vector machines, the state-of-the-art method for seizure detection, while consuming over 20x less power.

Emergence Patterns from General Anesthesia after Epilepsy Surgery: An Observational Pilot Study.

Objective Emergence from anesthesia starts from the limbic structures and then spreads outwards to brainstem, reticular activating systems, and then to the cortex. Epilepsy surgery often involves resection of limbic structures and hence may disrupt the pattern of emergence. The aim of this study was to explore the pattern of emergence from anesthesia following epilepsy surgery and to determine associated variables affecting the emergence pattern. Setting and Design Tertiary care center, prospective observational study. Materials and Methods We conducted a prospective observation pilot study on adult patients undergoing anterior temporal lobectomy and amygdalohippocampectomy for epilepsy. Anesthesia management was standardized in all patients, and they were allowed to wake up with "no touch" technique. Primary outcome of the study was the pattern of emergence (normal emergence, agitated emergence, or slow emergence) from anesthesia. Secondary outcomes were to explore the differences in preoperative neuropsychological profile and limbic structure volumes between the different patterns of emergence. Quantitative variables were analyzed using Student's t -test. Qualitative variables were analyzed using chi-square test. Results Twenty-nine patients completed the study: 9 patients (31%) had agitated emergence, and 20 patients had normal emergence. Among the agitated emergence, 2 patients had Riker scale of 7 indicating violent emergence. Patient demographics, anesthetic used, neuropsychological profile, and limbic structure volumes were similar between normal emergence and agitated emergence groups. However, two patients who had severe agitation (Riker scale of 7) had the lowest intelligence quotient. Conclusion Our pilot study showed that emergence agitation is not uncommon in patients undergoing epilepsy surgery. However, due to smaller sample size, the role of preoperative neuropsychologic profile and hippocampal volumes in predicting the pattern of emergence is inconclusive.

Abstract representations emerge in human hippocampal neurons during inference behavior.

Humans have the remarkable cognitive capacity to rapidly adapt to changing environments. Central to this capacity is the ability to form high-level, abstract representations that take advantage of regularities in the world to support generalization 1 . However, little is known about how these representations are encoded in populations of neurons, how they emerge through learning, and how they relate to behavior 2,3 . Here we characterized the representational geometry of populations of neurons (single-units) recorded in the hippocampus, amygdala, medial frontal cortex, and ventral temporal cortex of neurosurgical patients who are performing an inferential reasoning task. We find that only the neural representations formed in the hippocampus simultaneously encode multiple task variables in an abstract, or disentangled, format. This representational geometry is uniquely observed after patients learn to perform inference, and consisted of disentangled directly observable and discovered latent task variables. Interestingly, learning to perform inference by trial and error or through verbal instructions led to the formation of hippocampal representations with similar geometric properties. The observed relation between representational format and inference behavior suggests that abstract/disentangled representational geometries are important for complex cognition.

Distinctive biophysical features of human cell-types: insights from studies of neurosurgically resected brain tissue.

Electrophysiological characterization of live human tissue from epilepsy patients has been performed for many decades. Although initially these studies sought to understand the biophysical and synaptic changes associated with human epilepsy, recently, it has become the mainstay for exploring the distinctive biophysical and synaptic features of human cell-types. Both epochs of these human cellular electrophysiological explorations have faced criticism. Early studies revealed that cortical pyramidal neurons obtained from individuals with epilepsy appeared to function "normally" in comparison to neurons from non-epilepsy controls or neurons from other species and thus there was little to gain from the study of human neurons from epilepsy patients. On the other hand, contemporary studies are often questioned for the "normalcy" of the recorded neurons since they are derived from epilepsy patients. In this review, we discuss our current understanding of the distinct biophysical features of human cortical neurons and glia obtained from tissue removed from patients with epilepsy and tumors. We then explore the concept of within cell-type diversity and its loss (i.e., "neural homogenization"). We introduce neural homogenization to help reconcile the epileptogenicity of seemingly "normal" human cortical cells and circuits. We propose that there should be continued efforts to study cortical tissue from epilepsy patients in the quest to understand what makes human cell-types "human".

Hyperpolarization-Activated Cation Channels Shape the Spiking Frequency Preference of Human Cortical Layer 5 Pyramidal Neurons.

Discerning the contribution of specific ionic currents to complex neuronal dynamics is a difficult, but important, task. This challenge is exacerbated in the human setting, although the widely characterized uniqueness of the human brain compared with preclinical models necessitates the direct study of human neurons. Neuronal spiking frequency preference is of particular interest given its role in rhythm generation and signal transmission in cortical circuits. Here, we combine the frequency-dependent gain (FDG), a measure of spiking frequency preference, and novel in silico analyses to dissect the contributions of individual ionic currents to the suprathreshold features of human layer 5 (L5) neurons captured by the FDG. We confirm that a contemporary model of such a neuron, primarily constrained to capture subthreshold activity driven by the hyperpolarization-activated cyclic nucleotide gated (h-) current, replicates key features of the in vitro FDG both with and without h-current activity. With the model confirmed as a viable approximation of the biophysical features of interest, we applied new analysis techniques to quantify the activity of each modeled ionic current in the moments before spiking, revealing unique dynamics of the h-current. These findings motivated patch-clamp recordings in analogous rodent neurons to characterize their FDG, which confirmed that a biophysically detailed model of these neurons captures key interspecies differences in the FDG. These differences are correlated with distinct contributions of the h-current to neuronal activity. Together, this interdisciplinary and multispecies study provides new insights directly relating the dynamics of the h-current to suprathreshold spiking frequency preference in human L5 neurons.

Investigation of neurophysiologic and functional connectivity changes following glioma resection using magnetoencephalography.

Background: In patients with glioma, clinical manifestations of neural network disruption include behavioral changes, cognitive decline, and seizures. However, the extent of network recovery following surgery remains unclear. The aim of this study was to characterize the neurophysiologic and functional connectivity changes following glioma surgery using magnetoencephalography (MEG). Methods: Ten patients with newly diagnosed intra-axial brain tumors undergoing surgical resection were enrolled in the study and completed at least two MEG recordings (pre-operative and immediate post-operative). An additional post-operative recording 6-8 weeks following surgery was obtained for six patients. Resting-state MEG recordings from 28 healthy controls were used for network-based comparisons. MEG data processing involved artifact suppression, high-pass filtering, and source localization. Functional connectivity between parcellated brain regions was estimated using coherence values from 116 virtual channels. Statistical analysis involved standard parametric tests. Results: Distinct alterations in spectral power following tumor resection were observed, with at least three frequency bands affected across all study subjects. Tumor location-related changes were observed in specific frequency bands unique to each patient. Recovery of regional functional connectivity occurred following glioma resection, as determined by local coherence normalization. Changes in inter-regional functional connectivity were mapped across the brain, with comparable changes in low to mid gamma-associated functional connectivity noted in four patients. Conclusion: Our findings provide a framework for future studies to examine other network changes in glioma patients. We demonstrate an intrinsic capacity for neural network regeneration in the post-operative setting. Further work should be aimed at correlating neurophysiologic changes with individual patients' clinical outcomes.

Intrinsic neural diversity quenches the dynamic volatility of neural networks.

Heterogeneity is the norm in biology. The brain is no different: Neuronal cell types are myriad, reflected through their cellular morphology, type, excitability, connectivity motifs, and ion channel distributions. While this biophysical diversity enriches neural systems' dynamical repertoire, it remains challenging to reconcile with the robustness and persistence of brain function over time (resilience). To better understand the relationship between excitability heterogeneity (variability in excitability within a population of neurons) and resilience, we analyzed both analytically and numerically a nonlinear sparse neural network with balanced excitatory and inhibitory connections evolving over long time scales. Homogeneous networks demonstrated increases in excitability, and strong firing rate correlations-signs of instability-in response to a slowly varying modulatory fluctuation. Excitability heterogeneity tuned network stability in a context-dependent way by restraining responses to modulatory challenges and limiting firing rate correlations, while enriching dynamics during states of low modulatory drive. Excitability heterogeneity was found to implement a homeostatic control mechanism enhancing network resilience to changes in population size, connection probability, strength and variability of synaptic weights, by quenching the volatility (i.e., its susceptibility to critical transitions) of its dynamics. Together, these results highlight the fundamental role played by cell-to-cell heterogeneity in the robustness of brain function in the face of change.

Flexible 3D printed microwires and 3D microelectrodes for heart-on-a-chip engineering.

We developed a heart-on-a-chip platform that integrates highly flexible, vertical, 3D micropillar electrodes for electrophysiological recording and elastic microwires for the tissue's contractile force assessment. The high aspect ratio microelectrodes were 3D-printed into the device using a conductive polymer, poly(3,4-ethylenedioxythiophene):poly(styrene sulfonate) (PEDOT:PSS). A pair of flexible, quantum dots/thermoplastic elastomer nanocomposite microwires were 3D printed to anchor the tissue and enable continuous contractile force assessment. The 3D microelectrodes and flexible microwires enabled unobstructed human iPSC-based cardiac tissue formation and contraction, suspended above the device surface, under both spontaneous beating and upon pacing with a separate set of integrated carbon electrodes. Recording of extracellular field potentials using the PEDOT:PSS micropillars was demonstrated with and without epinephrine as a model drug, non-invasively, along within situmonitoring of tissue contractile properties and calcium transients. Uniquely, the platform provides integrated profiling of electrical and contractile tissue properties, which is critical for proper evaluation of complex, mechanically and electrically active tissues, such as the heart muscle under both physiological and pathological conditions.

Control of working memory maintenance by theta-gamma phase amplitude coupling of human hippocampal neurons.

Retaining information in working memory (WM) is a demanding process that relies on cognitive control to protect memoranda-specific persistent activity from interference. How cognitive control regulates WM storage, however, remains unknown. We hypothesized that interactions of frontal control and hippocampal persistent activity are coordinated by theta-gamma phase amplitude coupling (TG-PAC). We recorded single neurons in the human medial temporal and frontal lobe while patients maintained multiple items in WM. In the hippocampus, TG-PAC was indicative of WM load and quality. We identified cells that selectively spiked during nonlinear interactions of theta phase and gamma amplitude. These PAC neurons were more strongly coordinated with frontal theta activity when cognitive control demand was high, and they introduced information-enhancing and behaviorally relevant noise correlations with persistently active neurons in the hippocampus. We show that TG-PAC integrates cognitive control and WM storage to improve the fidelity of WM representations and facilitate behavior.

In-silico EEG biomarkers of reduced inhibition in human cortical microcircuits in depression.

Reduced cortical inhibition by somatostatin-expressing (SST) interneurons has been strongly associated with treatment-resistant depression. However, due to technical limitations it is impossible to establish experimentally in humans whether the effects of reduced SST interneuron inhibition on microcircuit activity have signatures detectable in clinically-relevant brain signals such as electroencephalography (EEG). To overcome these limitations, we simulated resting-state activity and EEG using detailed models of human cortical microcircuits with normal (healthy) or reduced SST interneuron inhibition (depression), and found that depression microcircuits exhibited increased theta, alpha and low beta power (4-16 Hz). The changes in depression involved a combination of an aperiodic broadband and periodic theta components. We then demonstrated the specificity of the EEG signatures of reduced SST interneuron inhibition by showing they were distinct from those corresponding to reduced parvalbumin-expressing (PV) interneuron inhibition. Our study thus links SST interneuron inhibition level to distinct features in EEG simulated from detailed human microcircuits, which can serve to better identify mechanistic subtypes of depression using EEG, and non-invasively monitor modulation of cortical inhibition.

Beamforming Seizures from the Temporal Lobe Using Magnetoencephalography.

BACKGROUND: Surgical treatment of drug-resistant temporal lobe epilepsy (TLE) depends on proper identification of the seizure onset zone (SOZ) and differentiation of mesial, temporolimbic seizure onsets from temporal neocortical seizure onsets. Noninvasive source imaging using electroencephalography (EEG) and magnetoencephalography (MEG) can provide accurate information on interictal spike localization; however, EEG and MEG have low sensitivity for epileptiform activity restricted to deep temporolimbic structures. Moreover, in mesial temporal lobe epilepsy (MTLE), interictal spikes frequently arise in neocortical foci distant from the SOZ, rendering interictal spike localization potentially misleading for presurgical planning. METHODS: In this study, we used two different beamformer techniques applied to the MEG signal of ictal events acquired during EEG-MEG recordings in six patients with TLE (three neocortical, three MTLE) in whom the ictal source localization results could be compared to ground truth SOZ localizations determined from intracranial EEG and/or clinical, neuroimaging, and postsurgical outcome evidence. RESULTS: Beamformer analysis proved to be highly accurate in all cases and was able to identify focal SOZs in mesial, temporolimbic structures. In three patients, interictal spikes were absent, too complex for dipole modeling, or localized to anterolateral temporal neocortex distant to a mesial temporal SOZ, and thus unhelpful in presurgical investigation. CONCLUSIONS: MEG beamformer source reconstruction is suitable for analysis of ictal events in TLE and can complement or supersede the traditional analysis of interictal spikes. The method outlined is applicable to any type of epileptiform event, expanding the information value of MEG and broadening its utility for presurgical recording in epilepsy.

Age-dependent increased sag amplitude in human pyramidal neurons dampens baseline cortical activity.

Aging involves various neurobiological changes, although their effect on brain function in humans remains poorly understood. The growing availability of human neuronal and circuit data provides opportunities for uncovering age-dependent changes of brain networks and for constraining models to predict consequences on brain activity. Here we found increased sag voltage amplitude in human middle temporal gyrus layer 5 pyramidal neurons from older subjects and captured this effect in biophysical models of younger and older pyramidal neurons. We used these models to simulate detailed layer 5 microcircuits and found lower baseline firing in older pyramidal neuron microcircuits, with minimal effect on response. We then validated the predicted reduced baseline firing using extracellular multielectrode recordings from human brain slices of different ages. Our results thus report changes in human pyramidal neuron input integration properties and provide fundamental insights into the neuronal mechanisms of altered cortical excitability and resting-state activity in human aging.

An in vitro whole-cell electrophysiology dataset of human cortical neurons.

BACKGROUND: Whole-cell patch-clamp electrophysiology is an essential technique for understanding how single neurons translate their diverse inputs into a functional output. The relative inaccessibility of live human cortical neurons for experimental manipulation has made it difficult to determine the unique features of how human cortical neurons differ from their counterparts in other species. FINDINGS: We present a curated repository of whole-cell patch-clamp recordings from surgically resected human cortical tissue, encompassing 118 neurons from 35 individuals (age range, 21-59 years; 17 male, 18 female). Recorded human cortical neurons derive from layers 2 and 3 (L2&3), deep layer 3 (L3c), or layer 5 (L5) and are annotated with a rich set of subject and experimental metadata. For comparison, we also provide a limited set of comparable recordings from 21-day-old mice (11 cells from 5 mice). All electrophysiological recordings are provided in the Neurodata Without Borders (NWB) format and are available for further analysis via the Distributed Archives for Neurophysiology Data Integration online repository. The associated data conversion code is made publicly available and can help others in converting electrophysiology datasets to the open NWB standard for general reuse. CONCLUSION: These data can be used for novel analyses of biophysical characteristics of human cortical neurons, including in cross-species or cross-lab comparisons or in building computational models of individual human neurons.