Because the human brain consumes a disproportionate fraction of the resting body's energy, positron emission tomography (PET) measurements of absolute glucose metabolism (CMRglc) can serve as disease biomarkers. Global mean normalization (GMN) of PET data reveals disease-based differences from healthy individuals as fractional changes across regions relative to a global mean. To assess the impact of GMN applied to metabolic data, we compared CMRglc with and without GMN in healthy awake volunteers with eyes closed (i.e., control) against specific physiological/clinical states, including healthy/awake with eyes open, healthy/awake but congenitally blind, healthy/sedated with anesthetics, and patients with disorders of consciousness. Without GMN, global CMRglc alterations compared to control were detected in all conditions except in congenitally blind where regional CMRglc variations were detected in the visual cortex. However, GMN introduced regional and bidirectional CMRglc changes at smaller fractions of the quantitative delocalized changes. While global information was lost with GMN, the quantitative approach (i.e., a validated method for quantitative baseline metabolic activity without GMN) not only preserved global CMRglc alterations induced by opening eyes, sedation, and varying consciousness but also detected regional CMRglc variations in the congenitally blind. These results caution the use of GMN upon PET-measured CMRglc data in health and disease.

Based on the hypothesis that brain plaques and tangles can affect cortical function in Alzheimer's disease (AD), we investigated functional responses in an AD rat model (called the Samaritan Alzheimer's rat achieved by ventricular infusion of amyloid peptide) and age-matched healthy control. High-field functional magnetic resonance imaging (fMRI) and extracellular neural activity measurements were applied to characterize sensory-evoked responses. Electrical stimulation of the forepaw led to BOLD and neural responses in the contralateral somatosensory cortex and thalamus. In AD brain we noted much smaller BOLD activation patterns in the somatosensory cortex (i.e., about 50% less activated voxels compared to normal brain). While magnitudes of BOLD and neural responses in the cerebral cortex were markedly attenuated in AD rats compared to normal rats (by about 50%), the dynamic coupling between the BOLD and neural responses in the cerebral cortex, as assessed by transfer function analysis, remained unaltered between the groups. However thalamic BOLD and neural responses were unaltered in AD brain compared to controls. Thus cortical responses in the AD model were indeed diminished compared to controls, but the thalamic responses in the AD and control rats were quite similar. Therefore these results suggest that Alzheimer's disease may affect cortical function more than subcortical function, which may have implications for interpreting altered human brain functional responses in fMRI studies of Alzheimer's disease.

Mitochondrial Ca(2+) uptake, central to neural metabolism and function, is diminished in aging whereas enhanced after acute/sub-acute traumatic brain injury. To develop relevant translational models for these neuropathologies, we determined the impact of perturbed mitochondrial Ca(2+) uptake capacities on intrinsic brain activity using clinically relevant markers. From a multi-compartment estimate of probable baseline Ca(2+) ranges in the brain, we hypothesized that reduced or enhanced mitochondrial Ca(2+) uptake capacity would decrease or increase spontaneous neuronal activity respectively. As resting state fMRI-BOLD fluctuations and stimulus-evoked BOLD responses have similar physiological origins [1] and stimulus-evoked neuronal and hemodynamic responses are modulated by mitochondrial Ca(2+) uptake capacity [2], [3] respectively, we tested our hypothesis by measuring hemodynamic fluctuations and spontaneous neuronal activities during normal and altered mitochondrial functional states. Mitochondrial Ca(2+) uptake capacity was perturbed by pharmacologically inhibiting or enhancing the mitochondrial Ca(2+) uniporter (mCU) activity. Neuronal electrical activity and cerebral blood flow (CBF) fluctuations were measured simultaneously and integrated with fMRI-BOLD fluctuations at 11.7T. mCU inhibition reduced spontaneous neuronal activity and the resting state functional connectivity (RSFC), whereas mCU enhancement increased spontaneous neuronal activity but reduced RSFC. We conclude that increased or decreased mitochondrial Ca(2+) uptake capacities lead to diminished resting state modes of brain functional connectivity.

The postnatal subventricular zone (SVZ) contains proliferating neural progenitor cells in close proximity to blood vessels. Insults and drug treatments acutely stimulate cell proliferation in the SVZ, which was assessed by labeling cells entering S phase. Although G1-to-S progression is metabolically demanding on a minute-to-hour time scale, it remains unknown whether increased SVZ cell proliferation is accompanied by a local hemodynamic response. This neurovascular coupling provides energy substrates to active neuronal assemblies. Transcardial dye perfusion revealed the presence of capillaries throughout the SVZ that constrict upon applications of the thromboxane A(2) receptor agonist U-46119 in acute brain slice preparations. We then monitored in vivo blood flow using laser Doppler flowmetry via a microprobe located either in the SVZ or a mature network. U-46119 injections into the lateral ventricle decreased blood flow in the SVZ and the striatum, which are near the ventricle. A 1-hour ventricular injection of epidermal and basic fibroblast growth factor (EGF and bFGF) significantly increased the percentage of Sox2 transcription factor-positive cells in S phase 1.5 hours post-injection. This increase was accompanied by a sustained rise in blood flow in the SVZ but not in the striatum. Direct growth factor injections into the cortex did not alter local blood flow, ruling out direct effects on capillaries. These findings suggest that an acute increase in the number of G1-to-S cycling SVZ cells is accompanied by neurometabolic-vascular coupling, which may provide energy and nutrient for cell cycle progression.

The metabolic and hemodynamic dependencies of the blood oxygenation level-dependent (BOLD) signal form the basis for calibrated fMRI, where the focus is on oxidative energy demanded by neural activity. An important part of calibrated fMRI is the power-law relationship between the BOLD signal and the deoxyhemoglobin concentration, which in turn is related to the ratio between oxidative demand (CMRO2) and blood flow (CBF). The power-law dependence between BOLD signal and deoxyhemoglobin concentration is signified by a scaling exponent β. Until recently most studies assumed a β value of 1.5, which is based on numerical simulations of the extravascular BOLD component. Since the basal value of CMRO2 and CBF can vary from subject-to-subject and/or region-to-region, a method to independently measure β in vivo should improve the accuracy of calibrated fMRI results. We describe a new method for β mapping through characterizing R2' - the most sensitive relaxation component of BOLD signal (i.e., the reversible magnetic susceptibility component that is predominantly of extravascular origin at high magnetic field) - as a function of intravascular magnetic susceptibility induced by an FDA-approved superparamagnetic contrast agent. In α-chloralose anesthetized rat brain, at 9.4 T, we measured β values of ~0.8 uniformly across large neocortical swathes, with lower magnitude and more heterogeneity in subcortical areas. Comparison of β maps in rats anesthetized with medetomidine and α-chloralose revealed that β is independent of neural activity levels at these resting states. We anticipate that this method for β mapping can help facilitate calibrated fMRI for clinical studies.

Odorants can reach olfactory receptor neurons (ORNs) by two routes: orthonasally, when volatiles enter the nasal cavity during inhalation/sniffing, and retronasally, when food volatiles released in the mouth pass into the nasal cavity during exhalation/eating. Previous work in humans has shown that both delivery routes of the same odorant can evoke distinct perceptions and patterns of neural responses in the brain. Each delivery route is known to influence specific responses across the dorsal region of the glomerular sheet in the olfactory bulb (OB), but spatial distributions across the entire glomerular sheet throughout the whole OB remain largely unexplored. We used functional MRI (fMRI) to measure and compare activations across the entire glomerular sheet in rat OB resulting from both orthonasal and retronasal stimulations of the same odors. We observed reproducible fMRI activation maps of the whole OB during both orthonasal and retronasal stimuli. However, retronasal stimuli required double the orthonasal odor concentration for similar response amplitudes. Regardless, both the magnitude and spatial extent of activity were larger during orthonasal versus retronasal stimuli for the same odor. Orthonasal and retronasal response patterns show overlap as well as some route-specific dominance. Orthonasal maps were dominant in dorsal-medial regions, whereas retronasal maps were dominant in caudal and lateral regions. These different whole OB encodings likely underlie differences in odor perception between these biologically important routes for odorants among mammals. These results establish the relationships between orthonasal and retronasal odor representations in the rat OB.

Physiological processes-such as, the brain's resting-state electrical activity or hemodynamic fluctuations-exhibit scale-free temporal structuring. However, impacts common in biological systems such as, noise, multiple signal generators, or filtering by transport function, result in multimodal scaling that cannot be reliably assessed by standard analytical tools that assume unimodal scaling. Here, we present two methods to identify breakpoints or crossovers in multimodal multifractal scaling functions. These methods incorporate the robust iterative fitting approach of the focus-based multifractal formalism (FMF). The first approach (moment-wise scaling range adaptivity) allows for a breakpoint-based adaptive treatment that analyzes segregated scale-invariant ranges. The second method (scaling function decomposition method, SFD) is a crossover-based design aimed at decomposing signal constituents from multimodal scaling functions resulting from signal addition or co-sampling, such as, contamination by uncorrelated fractals. We demonstrated that these methods could handle multimodal, mono- or multifractal, and exact or empirical signals alike. Their precision was numerically characterized on ideal signals, and a robust performance was demonstrated on exemplary empirical signals capturing resting-state brain dynamics by near infrared spectroscopy (NIRS), electroencephalography (EEG), and blood oxygen level-dependent functional magnetic resonance imaging (fMRI-BOLD). The NIRS and fMRI-BOLD low-frequency fluctuations were dominated by a multifractal component over an underlying biologically relevant random noise, thus forming a bimodal signal. The crossover between the EEG signal components was found at the boundary between the δ and θ bands, suggesting an independent generator for the multifractal δ rhythm. The robust implementation of the SFD method should be regarded as essential in the seamless processing of large volumes of bimodal fMRI-BOLD imaging data for the topology of multifractal metrics free of the masking effect of the underlying random noise.

Afferent nociceptive activity in the reorganizing spinal cord after SCI influences supraspinal regions to establish pain. Clinical evidence of poor motor functional recovery in SCI patients with pain, led us to hypothesize that sensory-motor integration transforms into sensory-motor interference to manifest pain. This was tested by investigating supraspinal changes in a rat model of hemicontusion cervical SCI. Animals displayed ipsilateral forelimb motor dysfunction and pain, which persisted at 6 weeks after SCI. Using resting state fMRI at 8 weeks after SCI, RSFC across 14 ROIs involved in nociception, indicated lateral differences with a relatively weaker right-right connectivity (deafferented-contralateral) compared to left-left (unaffected-ipsilateral). However, the sensory (S1) and motor (M1/M2) networks showed greater RSFC using right hemisphere ROI seeds when compared to left. Voxel seeds from the somatosensory forelimb (S1FL) and M1/M2 representations reproduced the SCI-induced sensory and motor RSFC enhancements observed using the ROI seeds. Larger local connectivity occurred in the right sensory and motor networks amidst a decreasing overall local connectivity. This maladaptive reorganization of the right (deafferented) hemisphere localized the sensory component of pain emerging from the ipsilateral forepaw. A significant expansion of the sensory and motor network s overlap occurred globally after SCI when compared to sham, supporting the hypothesis that sensory and motor interference manifests pain. Voxel-seed based analysis revealed greater sensory and motor network overlap in the left hemisphere when compared to the right. This left predominance of the overlap suggested relatively larger pain processing in the unaffected hemisphere, when compared to the deafferented side.

Large-scale functional networks have been characterized in both rodent and human brains, typically by analyzing fMRI-BOLD signals. However, the relationship between fMRI-BOLD and underlying neural activity is complex and incompletely understood, which poses challenges to interpreting network organization obtained using this technique. Additionally, most work has assumed a disjoint functional network organization (i.e., brain regions belong to one and only one network). Here, we employ wide-field Ca2+ imaging simultaneously with fMRI-BOLD in mice expressing GCaMP6f in excitatory neurons. We determine cortical networks discovered by each modality using a mixed-membership algorithm to test the hypothesis that functional networks exhibit overlapping organization. We find that there is considerable network overlap (both modalities) in addition to disjoint organization. Our results show that multiple BOLD networks are detected via Ca2+ signals, and networks determined by low-frequency Ca2+ signals are only modestly more similar to BOLD networks. In addition, the principal gradient of functional connectivity is nearly identical for BOLD and Ca2+ signals. Despite similarities, important differences are also detected across modalities, such as in measures of functional connectivity strength and diversity. In conclusion, Ca2+ imaging uncovers overlapping functional cortical organization in the mouse that reflects several, but not all, properties observed with fMRI-BOLD signals.

Calibrated fMRI extracts changes in oxidative energy demanded by neural activity based on hemodynamic and metabolic dependencies of the blood oxygenation level-dependent (BOLD) response. This procedure requires the parameter M, which is determined from the dynamic range of the BOLD signal between deoxyhemoglobin (paramagnetic) and oxyhemoglobin (diamagnetic). Since it is unclear if the range of M-values in human calibrated fMRI is due to regional/state differences, we conducted a 9.4T study to measure M-values across brain regions in deep (α-chloralose) and light (medetomidine) anesthetized rats, as verified by electrophysiology. Because BOLD signal is captured differentially by gradient-echo (R2*) and spin-echo (R2) relaxation rates, we measured M-values by the product of the fMRI echo time and R2' (i.e., the reversible magnetic susceptibility component), which is given by the absolute difference between R2* and R2. While R2' mapping was shown to be dependent on the k-space sampling method used, at nominal spatial resolutions achieved at high magnetic field of 9.4T the M-values were quite homogenous across cortical gray matter. However cortical M-values varied in relation to neural activity between brain states. The findings from this study could improve precision of future calibrated fMRI studies by focusing on the global uniformity of M-values in gray matter across different resting activity levels.

This article will be positioned on our previous work demonstrating the importance of adhering to a carefully selected set of criteria when choosing the suitable method from those available ensuring its adequate performance when applied to real temporal signals, such as fMRI BOLD, to evaluate one important facet of their behavior, fractality. Earlier, we have reviewed on a range of monofractal tools and evaluated their performance. Given the advance in the fractal field, in this article we will discuss the most widely used implementations of multifractal analyses, too. Our recommended flowchart for the fractal characterization of spontaneous, low frequency fluctuations in fMRI BOLD will be used as the framework for this article to make certain that it will provide a hands-on experience for the reader in handling the perplexed issues of fractal analysis. The reason why this particular signal modality and its fractal analysis has been chosen was due to its high impact on today's neuroscience given it had powerfully emerged as a new way of interpreting the complex functioning of the brain (see "intrinsic activity"). The reader will first be presented with the basic concepts of mono and multifractal time series analyses, followed by some of the most relevant implementations, characterization by numerical approaches. The notion of the dichotomy of fractional Gaussian noise and fractional Brownian motion signal classes and their impact on fractal time series analyses will be thoroughly discussed as the central theme of our application strategy. Sources of pitfalls and way how to avoid them will be identified followed by a demonstration on fractal studies of fMRI BOLD taken from the literature and that of our own in an attempt to consolidate the best practice in fractal analysis of empirical fMRI BOLD signals mapped throughout the brain as an exemplary case of potentially wide interest.

Fluctuations in spontaneous activity have been observed by many neuroimaging techniques, but because these resting-state changes are not evoked by stimuli, it is difficult to determine how they relate to task-evoked activations. We conducted multi-modal neuroimaging scans of the rat olfactory bulb, both with and without odor, to examine interaction between spontaneous and evoked activities. Independent component analysis of spontaneous fluctuations revealed resting-state networks, and odor-evoked changes revealed activation maps. We constructed simulated activation maps using resting-state networks that were highly correlated to evoked activation maps. Simulated activation maps derived by intrinsic optical signal (IOS), which covers the dorsal portion of the glomerular sheet, significantly differentiated one odor's evoked activation map from the other two. To test the hypothesis that spontaneous activity of the entire glomerular sheet is relevant for representing odor-evoked activations, we used functional magnetic resonance imaging (fMRI) to map the entire glomerular sheet. In contrast to the IOS results, the fMRI-derived simulated activation maps significantly differentiated all three odors' evoked activation maps. Importantly, no evoked activation maps could be significantly differentiated using simulated activation maps produced using phase-randomized resting-state networks. Given that some highly organized resting-state networks did not correlate with any odors' evoked activation maps, we posit that these resting-state networks may characterize evoked activation maps associated with odors not studied. These results emphasize that fluctuations in spontaneous activity form a foundation for active processing, signifying the relevance of resting-state mapping to functional neuroimaging.

Childhood maltreatment is associated with a wide range of psychopathologies including anxiety that emerge in childhood and in many cases persist in adulthood. Increased amygdala activation in response to threat and abnormal amygdala connectivity with frontolimbic brain regions, such as the hippocampus and the prefrontal cortex, are some of the most consistent findings seen in individuals exposed to childhood maltreatment. The underlying mechanisms responsible for these changes are difficult to study in humans but can be elucidated using animal models of early-life stress. Such studies are especially powerful in the mouse where precise control of the genetic background and the stress paradigm can be coupled with resting-state fMRI (rsfMRI) to map abnormal connectivity in circuits that regulate anxiety. To address this issue we first compared the effects of two models of early-life stress, limited bedding (LB) and unpredictable postnatal stress (UPS), on anxiety-like behavior in juvenile and adult mice. We found that UPS, but not LB, causes a robust increase in anxiety in juvenile and adult male mice. Next, we used rsfMRI to compare frontolimbic connectivity in control and UPS adult male mice. We found increased amygdala-prefrontal cortex and amygdala-hippocampus connectivity in UPS. The strength of the amygdala-hippocampal and amygdala-prefrontal cortex connectivity was highly correlated with anxiety-like behavior in the open-field test and elevated plus maze. These findings are the first to link hyperconnectivity in frontolimbic circuits and increased anxiety in a mouse model of early-life stress, allowing for more mechanistic understanding of parallel findings in humans.

Achieving a comprehensive understanding of brain function requires multiple imaging modalities with complementary strengths. We present an approach for concurrent widefield optical and functional magnetic resonance imaging. By merging these modalities, we can simultaneously acquire whole-brain blood-oxygen-level-dependent (BOLD) and whole-cortex calcium-sensitive fluorescent measures of brain activity. In a transgenic murine model, we show that calcium predicts the BOLD signal, using a model that optimizes a gamma-variant transfer function. We find consistent predictions across the cortex, which are best at low frequency (0.009-0.08 Hz). Furthermore, we show that the relationship between modality connectivity strengths varies by region. Our approach links cell-type-specific optical measurements of activity to the most widely used method for assessing human brain function.

Large-scale functional networks have been characterized in both rodent and human brains, typically by analyzing fMRI-BOLD signals. However, the relationship between fMRI-BOLD and underlying neural activity is complex and incompletely understood, which poses challenges to interpreting network organization obtained using this technique. Additionally, most work has assumed a disjoint functional network organization (i.e., brain regions belong to one and only one network). Here, we employed wide-field Ca2+ imaging simultaneously with fMRI-BOLD in mice expressing GCaMP6f in excitatory neurons. We determined cortical networks discovered by each modality using a mixed-membership algorithm to test the hypothesis that functional networks are overlapping rather than disjoint. Our results show that multiple BOLD networks are detected via Ca2+ signals; there is considerable network overlap (both modalities); networks determined by low-frequency Ca2+ signals are only modestly more similar to BOLD networks; and, despite similarities, important differences are detected across modalities (e.g., brain region "network diversity"). In conclusion, Ca2+ imaging uncovered overlapping functional cortical organization in the mouse that reflected several, but not all, properties observed with fMRI-BOLD signals.

Schizophrenia is a debilitating neuropsychiatric disorder typically diagnosed from late adolescence to adulthood. Subthreshold behavioral symptoms (e.g., cognitive deficits and substance abuse) often precede the clinical diagnosis of schizophrenia. However, these prodromal symptoms have not been consistently associated with structural and functional brain biomarkers, limiting the chance of early diagnosis of schizophrenia.

Impaired consciousness in temporal lobe seizures has a major negative impact on quality of life. The prevailing view holds that this disorder impairs consciousness by seizure spread to the bilateral temporal lobes. We propose instead that seizures invade subcortical regions and depress arousal, causing impairment through decreases rather than through increases in activity. Using functional magnetic resonance imaging in a rodent model, we found increased activity in regions known to depress cortical function, including lateral septum and anterior hypothalamus. Importantly, we found suppression of intralaminar thalamic and brainstem arousal systems and suppression of the cortex. At a cellular level, we found reduced firing of identified cholinergic neurons in the brainstem pedunculopontine tegmental nucleus and basal forebrain. Finally, we used enzyme-based amperometry to demonstrate reduced cholinergic neurotransmission in both cortex and thalamus. Decreased subcortical arousal is a critical mechanism for loss of consciousness in focal temporal lobe seizures.

Fluctuations in resting-state cerebral hemodynamics show scale-free behavior over two distinct scaling ranges. Changes in such bimodal (multi) fractal pattern give insight to altered cerebrovascular or neural function. Our main goal was to assess the distribution of local scale-free properties characterizing cerebral hemodynamics and to disentangle the influence of aging on these multifractal parameters. To this end, we obtained extended resting-state records (N = 214) of oxyhemoglobin (HbO), deoxyhemoglobin (HbR) and total hemoglobin (HbT) concentration time series with continuous-wave near-infrared spectroscopy technology from the brain cortex. 52 healthy volunteers were enrolled in this study: 24 young (30.6 ± 8.2 years), and 28 elderly (60.5 ± 12.0 years) subjects. Using screening tests on power-law, multifractal noise, and shuffled data sets we evaluated the presence of true multifractal hemodynamics reflecting long-range correlation (LRC). Subsequently, scaling-range adaptive bimodal signal summation conversion (SSC) was performed based on standard deviation (σ) of signal windows across a range of temporal scales (s). Building on moments of different order (q) of the measure, σ(s), multifractal SSC yielded generalized Hurst exponent function, H(q), and singularity spectrum, D(h) separately for a fast and slow component (the latter dominating the highest temporal scales). Parameters were calculated reflecting the estimated measure at s = N (focus), degree of LRC [Hurst exponent, H(2) and maximal Hölder exponent, hmax] and measuring strength of multifractality [full-width-half-maximum of D(h) and ΔH15 = H(-15)-H(15)]. Correlation-based signal improvement (CBSI) enhanced our signal in terms of interpreting changes due to neural activity or local/systemic hemodynamic influences. We characterized the HbO-HbR relationship with the aid of fractal scale-wise correlation coefficient, rσ(s) and SSC-based multifractal covariance analysis. In the majority of subjects, cerebral hemodynamic fluctuations proved bimodal multifractal. In case of slow component of raw HbT, hmax, and Ĥ(2) were lower in the young group explained by a significantly increased rσ(s) among elderly at high temporal scales. Regarding the fast component of CBSI-pretreated HbT and that of HbO-HbR covariance, hmax, and focus were decreased in the elderly group. These observations suggest an attenuation of neurovascular coupling reflected by a decreased autocorrelation of the neuronal component concomitant with an accompanying increased autocorrelation of the non-neuronal component in the elderly group.

Functional magnetic resonance imaging (fMRI) techniques using the blood-oxygen level-dependent (BOLD) signal have shown great potential as clinical biomarkers of disease. Thus, using these techniques in preclinical rodent models is an urgent need. Calibrated fMRI is a promising technique that can provide high-resolution mapping of cerebral oxygen metabolism (CMRO2). However, calibrated fMRI is difficult to use in rodent models for several reasons: rodents are anesthetized, stimulation-induced changes are small, and gas challenges induce noisy CMRO2 predictions. We used, in mice, a relaxometry-based calibrated fMRI method which uses cerebral blood flow (CBF) and the BOLD-sensitive magnetic relaxation component, R2', the same parameter derived in the deoxyhemoglobin-dilution model of calibrated fMRI. This method does not use any gas challenges, which we tested on mice in both awake and anesthetized states. As anesthesia induces a whole-brain change, our protocol allowed us to overcome the former limitations of rodent studies using calibrated fMRI. We revealed 1.5-2 times higher CMRO2, dependent upon brain region, in the awake state versus the anesthetized state. Our results agree with alternative measurements of whole-brain CMRO2 in the same mice and previous human anesthesia studies. The use of calibrated fMRI in rodents has much potential for preclinical fMRI.

Absence seizures are brief episodes of impaired consciousness, behavioral arrest, and unresponsiveness, with yet-unknown neuronal mechanisms. Here we report that an awake female rat model recapitulates the behavioral, electroencephalographic, and cortical functional magnetic resonance imaging characteristics of human absence seizures. Neuronally, seizures feature overall decreased but rhythmic firing of neurons in cortex and thalamus. Individual cortical and thalamic neurons express one of four distinct patterns of seizure-associated activity, one of which causes a transient initial peak in overall firing at seizure onset, and another which drives sustained decreases in overall firing. 40-60 s before seizure onset there begins a decline in low frequency electroencephalographic activity, neuronal firing, and behavior, but an increase in higher frequency electroencephalography and rhythmicity of neuronal firing. Our findings demonstrate that prolonged brain state changes precede consciousness-impairing seizures, and that during seizures distinct functional groups of cortical and thalamic neurons produce an overall transient firing increase followed by a sustained firing decrease, and increased rhythmicity.