Pharmacological Activation of Piezo1 Channels Enhances Astrocyte-Neuron Communication via NMDA Receptors in the Murine Neocortex.

The Piezo1 mechanosensitive ion channel is abundant on several elements of the central nervous system including astrocytes. It has been already demonstrated that activation of these channels is able to elicit calcium waves on astrocytes, which contributes to the release of gliotransmitters. Astrocyte- and N-methyl-D-aspartate (NMDA) receptor-dependent slow inward currents (SICs) are hallmarks of astrocyte-neuron communication. These currents are triggered by glutamate released as gliotransmitter, which in turn activates neuronal NMDA receptors responsible for this inward current having slower kinetics than any synaptic events. In this project, we aimed to investigate whether Piezo1 activation and inhibition is able to alter spontaneous SIC activity of murine neocortical pyramidal neurons. When the Piezo1 opener Yoda1 was applied, the SIC frequency and the charge transfer by these events in a minute time was significantly increased. These changes were prevented by treating the preparations with the NMDA receptor inhibitor D-AP5. Furthermore, Yoda1 did not alter the spontaneous EPSC frequency and amplitude when SICs were absent. The Piezo1 inhibitor Dooku1 effectively reverted the actions of Yoda1 and decreased the rise time of SICs when applied alone. In conclusion, activation of Piezo1 channels is able to alter astrocyte-neuron communication. Via enhancement of SIC activity, astrocytic Piezo1 channels have the capacity to determine neuronal excitability.

Cell-type-resolved mosaicism reveals clonal dynamics of the human forebrain.

Debate remains around the anatomical origins of specific brain cell subtypes and lineage relationships within the human forebrain1-7. Thus, direct observation in the mature human brain is critical for a complete understanding of its structural organization and cellular origins. Here we utilize brain mosaic variation within specific cell types as distinct indicators for clonal dynamics, denoted as cell-type-specific mosaic variant barcode analysis. From four hemispheres and two different human neurotypical donors, we identified 287 and 780 mosaic variants, respectively, that were used to deconvolve clonal dynamics. Clonal spread and allele fractions within the brain reveal that local hippocampal excitatory neurons are more lineage-restricted than resident neocortical excitatory neurons or resident basal ganglia GABAergic inhibitory neurons. Furthermore, simultaneous genome transcriptome analysis at both a cell-type-specific and a single-cell level suggests a dorsal neocortical origin for a subgroup of DLX1+ inhibitory neurons that disperse radially from an origin shared with excitatory neurons. Finally, the distribution of mosaic variants across 17 locations within one parietal lobe reveals that restriction of clonal spread in the anterior-posterior axis precedes restriction in the dorsal-ventral axis for both excitatory and inhibitory neurons. Thus, cell-type-resolved somatic mosaicism can uncover lineage relationships governing the development of the human forebrain.

Activity-Dependent Ectopic Spiking in Parvalbumin-Expressing Interneurons of the Neocortex.

Canonically, action potentials of most mammalian neurons initiate at the axon initial segment (AIS) and propagate bidirectionally: orthodromically along the distal axon and retrogradely into the soma and dendrites. Under some circumstances, action potentials may initiate ectopically, at sites distal to the AIS, and propagate antidromically along the axon. These "ectopic action potentials" (EAPs) have been observed in experimental models of seizures and chronic pain, and more rarely in nonpathological forebrain neurons. Here we report that a large majority of parvalbumin-expressing (PV+) interneurons in the upper layers of mouse neocortex, from both orbitofrontal and primary somatosensory areas, fire EAPs after sufficient activation of their somata. Somatostatin-expressing interneurons also fire EAPs, though less robustly. Ectopic firing in PV+ cells occurs in varying temporal patterns and can persist for several seconds. PV+ cells evoke strong synaptic inhibition in pyramidal neurons and interneurons and play critical roles in cortical function. Our results suggest that ectopic spiking of PV+ interneurons is common and may contribute to both normal and pathological network functions of the neocortex.

Morphological correlates of pyramidal cell axonal myelination in mouse and human neocortex.

The axons of neocortical pyramidal neurons are frequently myelinated. Heterogeneity in the topography of axonal myelination in the cerebral cortex has been attributed to a combination of electrophysiological activity, axonal morphology, and neuronal-glial interactions. Previously, we showed that axonal segment length and caliber are critical local determinants of fast-spiking interneuron myelination. However, the factors that determine the myelination of individual axonal segments along neocortical pyramidal neurons remain largely unexplored. Here, we used structured illumination microscopy to examine the extent to which axonal morphology is predictive of the topography of myelination along neocortical pyramidal neurons. We identified critical thresholds for axonal caliber and interbranch distance that are necessary, but not sufficient, for myelination of pyramidal cell axons in mouse primary somatosensory cortex (S1). Specifically, we found that pyramidal neuron axonal segments with a caliber < 0.24 µm or interbranch distance < 18.10 µm are rarely myelinated. Moreover, we further confirmed that these findings in mice are similar for human neocortical pyramidal cell myelination (caliber < 0.25 µm, interbranch distance < 19.00 µm), suggesting that this mechanism is evolutionarily conserved. Taken together, our findings suggest that axonal morphology is a critical correlate of the topography and cell-type specificity of neocortical myelination.

Developmental trajectories of GABAergic cortical interneurons are sequentially modulated by dynamic FoxG1 expression levels.

GABAergic inhibitory interneurons, originating from the embryonic ventral forebrain territories, traverse a convoluted migratory path to reach the neocortex. These interneuron precursors undergo sequential phases of tangential and radial migration before settling into specific laminae during differentiation. Here, we show that the developmental trajectory of FoxG1 expression is dynamically controlled in these interneuron precursors at critical junctures of migration. By utilizing mouse genetic strategies, we elucidate the pivotal role of precise changes in FoxG1 expression levels during interneuron specification and migration. Our findings underscore the gene dosage-dependent function of FoxG1, aligning with clinical observations of FOXG1 haploinsufficiency and duplication in syndromic forms of autism spectrum disorders. In conclusion, our results reveal the finely tuned developmental clock governing cortical interneuron development, driven by temporal dynamics and the dose-dependent actions of FoxG1.

Brain Sample Preparation After Performing in Utero Electroporation to Proliferating and Differentiated Cells in the Developing Mouse Neocortex.

In utero electroporation (IUE) enables labeling and manipulating specific types of cells by introducing DNA plasmids with desired promoters. After the surgery, mouse brains are fixed at any stage and analyzed after staining using specific antibodies. Here, we describe the flow of the IUE experiment from the preparation to microscopic observations.

Conditioning and pseudoconditioning differently change intrinsic excitability of inhibitory interneurons in the neocortex.

Many studies indicate a broad role of various classes of GABAergic interneurons in the processes related to learning. However, little is known about how the learning process affects intrinsic excitability of specific classes of interneurons in the neocortex. To determine this, we employed a simple model of conditional learning in mice where vibrissae stimulation was used as a conditioned stimulus and a tail shock as an unconditioned one. In vitro whole-cell patch-clamp recordings showed an increase in intrinsic excitability of low-threshold spiking somatostatin-expressing interneurons (SST-INs) in layer 4 (L4) of the somatosensory (barrel) cortex after the conditioning paradigm. In contrast, pseudoconditioning reduced intrinsic excitability of SST-LTS, parvalbumin-expressing interneurons (PV-INs), and vasoactive intestinal polypeptide-expressing interneurons (VIP-INs) with accommodating pattern in L4 of the barrel cortex. In general, increased intrinsic excitability was accompanied by narrowing of action potentials (APs), whereas decreased intrinsic excitability coincided with AP broadening. Altogether, these results show that both conditioning and pseudoconditioning lead to plastic changes in intrinsic excitability of GABAergic interneurons in a cell-specific manner. In this way, changes in intrinsic excitability can be perceived as a common mechanism of learning-induced plasticity in the GABAergic system.

Female-specific dysfunction of sensory neocortical circuits in a mouse model of autism mediated by mGluR5 and estrogen receptor α.

Little is known of the brain mechanisms that mediate sex-specific autism symptoms. Here, we demonstrate that deletion of the autism spectrum disorder (ASD)-risk gene, Pten, in neocortical pyramidal neurons (NSEPten knockout [KO]) results in robust cortical circuit hyperexcitability selectively in female mice observed as prolonged spontaneous persistent activity states. Circuit hyperexcitability in females is mediated by metabotropic glutamate receptor 5 (mGluR5) and estrogen receptor α (ERα) signaling to mitogen-activated protein kinases (Erk1/2) and de novo protein synthesis. Pten KO layer 5 neurons have a female-specific increase in mGluR5 and mGluR5-dependent protein synthesis. Furthermore, mGluR5-ERα complexes are generally elevated in female cortices, and genetic reduction of ERα rescues enhanced circuit excitability, protein synthesis, and neuron size selectively in NSEPten KO females. Female NSEPten KO mice display deficits in sensory processing and social behaviors as well as mGluR5-dependent seizures. These results reveal mechanisms by which sex and a high-confidence ASD-risk gene interact to affect brain function and behavior.

New Empirical Bayes Models to Jointly Analyze Multiple RNA-Sequencing Data in a Hypophosphatasia Disease Study.

Hypophosphatasia is a rare inherited metabolic disorder caused by the deficiency of tissue-nonspecific alkaline phosphatase. More severe and early onset cases present symptoms of muscle weakness, diminished motor coordination, and epileptic seizures. These neurological manifestations are poorly characterized. Thus, it is urgent to discover novel differentially expressed genes for investigating the genetic mechanisms underlying the neurological manifestations of hypophosphatasia. RNA-sequencing data offer a high-resolution and highly accurate transcript profile. In this study, we apply an empirical Bayes model to RNA-sequencing data acquired from the spinal cord and neocortex tissues of a mouse model, individually, to more accurately estimate the genetic effects without bias. More importantly, we further develop two integration methods, weighted gene approach and weighted Z method, to incorporate two RNA-sequencing data into a model for enhancing the effects of genetic markers in the diagnostics of hypophosphatasia disease. The simulation and real data analysis have demonstrated the effectiveness of our proposed integration methods, which can maximize genetic signals identified from the spinal cord and neocortex tissues, minimize the prediction error, and largely improve the prediction accuracy in risk prediction.

Functional synergy of a human-specific and an ape-specific metabolic regulator in human neocortex development.

Metabolism has recently emerged as a major target of genes implicated in the evolutionary expansion of human neocortex. One such gene is the human-specific gene ARHGAP11B. During human neocortex development, ARHGAP11B increases the abundance of basal radial glia, key progenitors for neocortex expansion, by stimulating glutaminolysis (glutamine-to-glutamate-to-alpha-ketoglutarate) in mitochondria. Here we show that the ape-specific protein GLUD2 (glutamate dehydrogenase 2), which also operates in mitochondria and converts glutamate-to-αKG, enhances ARHGAP11B's ability to increase basal radial glia abundance. ARHGAP11B + GLUD2 double-transgenic bRG show increased production of aspartate, a metabolite essential for cell proliferation, from glutamate via alpha-ketoglutarate and the TCA cycle. Hence, during human evolution, a human-specific gene exploited the existence of another gene that emerged during ape evolution, to increase, via concerted changes in metabolism, progenitor abundance and neocortex size.

A coupled neural field model for the standard consolidation theory.

The standard consolidation theory states that short-term memories located in the hippocampus enable the consolidation of long-term memories in the neocortex. In other words, the neocortex slowly learns long-term memories with a transient support of the hippocampus that quickly learns unstable memories. However, it is not clear yet what could be the neurobiological mechanisms underlying these differences in learning rates and memory time-scales. Here, we propose a novel modeling approach of the standard consolidation theory, that focuses on its potential neurobiological mechanisms. In addition to synaptic plasticity and spike frequency adaptation, our model incorporates adult neurogenesis in the dentate gyrus as well as the difference in size between the neocortex and the hippocampus, that we associate with distance-dependent synaptic plasticity. We also take into account the interconnected spatial structure of the involved brain areas, by incorporating the above neurobiological mechanisms in a coupled neural field framework, where each area is represented by a separate neural field with intra- and inter-area connections. To our knowledge, this is the first attempt to apply neural fields to this process. Using numerical simulations and mathematical analysis, we explore the short-term and long-term dynamics of the model upon alternance of phases of hippocampal replay and retrieval cue of an external input. This external input is encodable as a memory pattern in the form of a multiple bump attractor pattern in the individual neural fields. In the model, hippocampal memory patterns become encoded first, before neocortical ones, because of the smaller distances between the bumps of the hippocampal memory patterns. As a result, retrieval of the input pattern in the neocortex at short time-scales necessitates the additional input delivered by the memory pattern of the hippocampus. Neocortical memory patterns progressively consolidate at longer times, up to a point where their retrieval does not need the support of the hippocampus anymore. At longer times, perturbation of the hippocampal neural fields by neurogenesis erases the hippocampus pattern, leading to a final state where the memory pattern is exclusively evoked in the neocortex. Therefore, the dynamics of our model successfully reproduces the main features of the standard consolidation theory. This suggests that neurogenesis in the hippocampus and distance-dependent synaptic plasticity coupled to synaptic depression and spike frequency adaptation, are indeed critical neurobiological processes in memory consolidation.

Progressive Circuit Hyperexcitability in Mouse Neocortical Slice Cultures with Increasing Duration of Activity Silencing.

Forebrain neurons deprived of activity become hyperactive when activity is restored. Rebound activity has been linked to spontaneous seizures in vivo following prolonged activity blockade. Here, we measured the time course of rebound activity and the contributing circuit mechanisms using calcium imaging, synaptic staining, and whole-cell patch clamp in organotypic slice cultures of mouse neocortex. Calcium imaging revealed hypersynchronous activity increasing in intensity with longer periods of deprivation. While activity partially recovered 3 d after slices were released from 5 d of deprivation, they were less able to recover after 10 d of deprivation. However, even after the longer period of deprivation, activity patterns eventually returned to baseline levels. The degree of deprivation-induced rebound was age-dependent, with the greatest effects occurring when silencing began in the second week. Pharmacological blockade of NMDA receptors indicated that hypersynchronous rebound activity did not require activation of Hebbian plasticity. In single-neuron recordings, input resistance roughly doubled with a concomitant increase in intrinsic excitability. Synaptic imaging of pre- and postsynaptic proteins revealed dramatic reductions in the number of presumptive synapses with a larger effect on inhibitory than excitatory synapses. Putative excitatory synapses colocalizing PSD-95 and Bassoon declined by 39 and 56% following 5 and 10 d of deprivation, but presumptive inhibitory synapses colocalizing gephyrin and VGAT declined by 55 and 73%, respectively. The results suggest that with prolonged deprivation, a progressive reduction in synapse number is accompanied by a shift in the balance between excitation and inhibition and increased cellular excitability.

Control of neuronal excitation-inhibition balance by BMP-SMAD1 signalling.

Throughout life, neuronal networks in the mammalian neocortex maintain a balance of excitation and inhibition, which is essential for neuronal computation1,2. Deviations from a balanced state have been linked to neurodevelopmental disorders, and severe disruptions result in epilepsy3-5. To maintain balance, neuronal microcircuits composed of excitatory and inhibitory neurons sense alterations in neural activity and adjust neuronal connectivity and function. Here we identify a signalling pathway in the adult mouse neocortex that is activated in response to increased neuronal network activity. Overactivation of excitatory neurons is signalled to the network through an increase in the levels of BMP2, a growth factor that is well known for its role as a morphogen in embryonic development. BMP2 acts on parvalbumin-expressing (PV) interneurons through the transcription factor SMAD1, which controls an array of glutamatergic synapse proteins and components of perineuronal nets. PV-interneuron-specific disruption of BMP2-SMAD1 signalling is accompanied by a loss of glutamatergic innervation in PV cells, underdeveloped perineuronal nets and decreased excitability. Ultimately, this impairment of the functional recruitment of PV interneurons disrupts the cortical excitation-inhibition balance, with mice exhibiting spontaneous epileptic seizures. Our findings suggest that developmental morphogen signalling is repurposed to stabilize cortical networks in the adult mammalian brain.

High Contrast PET Imaging of Subcortical and Allocortical Amyloid-ß in Early Alzheimer's Disease Using [11C]AZD2184.

Background: Deposits of amyloid-ß (Aß) appear early in Alzheimer's disease (AD). Objective: The aim of the present study was to compare the presence of cortical and subcortical Aß in early AD using positron emission tomography (PET). Methods: Eight cognitively unimpaired (CU) subjects, 8 with mild cognitive impairment (MCI) and 8 with mild AD were examined with PET and [11C]AZD2184. A data driven cut-point for Aß positivity was defined by Gaussian mixture model of isocortex binding potential (BPND) values. Results: Sixteen subjects (3 CU, 5 MCI and 8 AD) were Aß-positive. BPND was lower in subcortical and allocortical regions compared to isocortex. Fifteen of the 16 Aß-positive subjects displayed Aß binding in striatum, 14 in thalamus and 10 in allocortical regions. Conclusions: Aß deposits appear to be widespread in early AD. It cannot be excluded that deposits appear simultaneously throughout the whole brain which has implications for improved diagnostics and disease monitoring.

Hippocampal ripples coincide with "up-state" and spindles in retrosplenial cortex.

During NREM sleep, hippocampal sharp-wave ripple (SWR) events are thought to stabilize memory traces for long-term storage in downstream neocortical structures. Within the neocortex, a set of distributed networks organized around retrosplenial cortex (RS-network) interact preferentially with the hippocampus purportedly to consolidate those traces. Transient bouts of slow oscillations and sleep spindles in this RS-network are often observed around SWRs, suggesting that these two activities are related and that their interplay possibly contributes to memory consolidation. To investigate how SWRs interact with the RS-network and spindles, we combined cortical wide-field voltage imaging, Electrocorticography, and hippocampal LFP recordings in anesthetized and sleeping mice. Here, we show that, during SWR, "up-states" and spindles reliably co-occur in a cortical subnetwork centered around the retrosplenial cortex. Furthermore, retrosplenial transient activations and spindles predict slow gamma oscillations in CA1 during SWRs. Together, our results suggest that retrosplenial-hippocampal interaction may be a critical pathway of information exchange between the cortex and hippocampus.

Cannabidiol Exerts Anticonvulsant Effects Alone and in Combination with Δ9-THC through the 5-HT1A Receptor in the Neocortex of Mice.

Cannabinoids have shown potential in drug-resistant epilepsy treatment; however, we lack knowledge on which cannabinoid(s) to use, dosing, and their pharmacological targets. This study investigated (i) the anticonvulsant effect of Cannabidiol (CBD) alone and (ii) in combination with Delta-9 Tetrahydrocannabinol (Δ9-THC), as well as (iii) the serotonin (5-HT)1A receptor's role in CBD's mechanism of action. Seizure activity, induced by 4-aminopyridine, was measured by extracellular field recordings in cortex layer 2/3 of mouse brain slices. The anticonvulsant effect of 10, 30, and 100 µM CBD alone and combined with Δ9-THC was evaluated. To examine CBD's mechanism of action, slices were pre-treated with a 5-HT1A receptor antagonist before CBD's effect was evaluated. An amount of ≥30 µM CBD alone exerted significant anticonvulsant effects while 10 µM CBD did not. However, 10 µM CBD combined with low-dose Δ9-THC (20:3 ratio) displayed significantly greater anticonvulsant effects than either phytocannabinoid alone. Furthermore, blocking 5-HT1A receptors before CBD application significantly abolished CBD's effects. Thus, our results demonstrate the efficacy of low-dose CBD and Δ9-THC combined and that CBD exerts its effects, at least in part, through 5-HT1A receptors. These results could address drug-resistance while providing insight into CBD's mechanism of action, laying the groundwork for further testing of cannabinoids as anticonvulsants.

The growth factor EPIREGULIN promotes basal progenitor cell proliferation in the developing neocortex.

Neocortex expansion during evolution is linked to higher numbers of neurons, which are thought to result from increased proliferative capacity and neurogenic potential of basal progenitor cells during development. Here, we show that EREG, encoding the growth factor EPIREGULIN, is expressed in the human developing neocortex and in gorilla cerebral organoids, but not in the mouse neocortex. Addition of EPIREGULIN to the mouse neocortex increases proliferation of basal progenitor cells, whereas EREG ablation in human cortical organoids reduces proliferation in the subventricular zone. Treatment of cortical organoids with EPIREGULIN promotes a further increase in proliferation of gorilla but not of human basal progenitor cells. EPIREGULIN competes with the epidermal growth factor (EGF) to promote proliferation, and inhibition of the EGF receptor abrogates the EPIREGULIN-mediated increase in basal progenitor cells. Finally, we identify putative cis-regulatory elements that may contribute to the observed inter-species differences in EREG expression. Our findings suggest that species-specific regulation of EPIREGULIN expression may contribute to the increased neocortex size of primates by providing a tunable pro-proliferative signal to basal progenitor cells in the subventricular zone.

Rapid volume pulsations of the extracellular space accompany epileptiform activity in trauma-injured neocortex and depend on the sodium-bicarbonate cotransporter NBCe1.

Post traumatic epilepsy (PTE) is a treatment-resistant consequence of traumatic brain injury (TBI). Recently, it has been revealed that epileptiform activity in acute chemoconvulsant seizure models is accompanied by transient shrinkages of extracellular space (ECS) called rapid volume pulsations (RVPs). Shrinkage of the ECS surrounding neurons and glia may contribute to ictogenic hyperexcitability and hypersynchrony during the chronic phase of TBI. Here, we identify the phenomenon of RVPs occurring spontaneously in rat neocortex at ≥ 3 weeks after injury in the controlled cortical impact (CCI) model for PTE. We further report that blocking the electrogenic action of the astrocytic cotransporter NBCe1 with 4,4'-diisothiocyano-2,2'-stilbenedisulfonic acid (DIDS) eliminates both RVPs and epileptiform activity in ex-vivo CCI neocortical brain slices. We conclude that NBCe1-mediated extracellular volume shrinkage may represent a new target for therapeutic intervention in PTE.

An operating principle of the cerebral cortex, and a cellular mechanism for attentional trial-and-error pattern learning and useful classification extraction.

A feature of the brains of intelligent animals is the ability to learn to respond to an ensemble of active neuronal inputs with a behaviorally appropriate ensemble of active neuronal outputs. Previously, a hypothesis was proposed on how this mechanism is implemented at the cellular level within the neocortical pyramidal neuron: the apical tuft or perisomatic inputs initiate "guess" neuron firings, while the basal dendrites identify input patterns based on excited synaptic clusters, with the cluster excitation strength adjusted based on reward feedback. This simple mechanism allows neurons to learn to classify their inputs in a surprisingly intelligent manner. Here, we revise and extend this hypothesis. We modify synaptic plasticity rules to align with behavioral time scale synaptic plasticity (BTSP) observed in hippocampal area CA1, making the framework more biophysically and behaviorally plausible. The neurons for the guess firings are selected in a voluntary manner via feedback connections to apical tufts in the neocortical layer 1, leading to dendritic Ca2+ spikes with burst firing, which are postulated to be neural correlates of attentional, aware processing. Once learned, the neuronal input classification is executed without voluntary or conscious control, enabling hierarchical incremental learning of classifications that is effective in our inherently classifiable world. In addition to voluntary, we propose that pyramidal neuron burst firing can be involuntary, also initiated via apical tuft inputs, drawing attention toward important cues such as novelty and noxious stimuli. We classify the excitations of neocortical pyramidal neurons into four categories based on their excitation pathway: attentional versus automatic and voluntary/acquired versus involuntary. Additionally, we hypothesize that dendrites within pyramidal neuron minicolumn bundles are coupled via depolarization cross-induction, enabling minicolumn functions such as the creation of powerful hierarchical "hyperneurons" and the internal representation of the external world. We suggest building blocks to extend the microcircuit theory to network-level processing, which, interestingly, yields variants resembling the artificial neural networks currently in use. On a more speculative note, we conjecture that principles of intelligence in universes governed by certain types of physical laws might resemble ours.

Epigenomic tomography for probing spatially defined chromatin state in the brain.

Spatially resolved epigenomic profiling is critical for understanding biology in the mammalian brain. Single-cell spatial epigenomic assays were developed recently for this purpose, but they remain costly and labor intensive for examining brain tissues across substantial dimensions and surveying a collection of brain samples. Here, we demonstrate an approach, epigenomic tomography, that maps spatial epigenomes of mouse brain at the scale of centimeters. We individually profiled neuronal and glial fractions of mouse neocortex slices with 0.5 mm thickness. Tri-methylation of histone 3 at lysine 27 (H3K27me3) or acetylation of histone 3 at lysine 27 (H3K27ac) features across these slices were grouped into clusters based on their spatial variation patterns to form epigenomic brain maps. As a proof of principle, our approach reveals striking dynamics in the frontal cortex due to kainic-acid-induced seizure, linked with transmembrane ion transporters, exocytosis of synaptic vesicles, and secretion of neurotransmitters. Epigenomic tomography provides a powerful and cost-effective tool for characterizing brain disorders based on the spatial epigenome.