VEGF-C prophylaxis favors lymphatic drainage and modulates neuroinflammation in a stroke model

Meningeal lymphatic vessels (MLVs) promote tissue clearance and immune surveillance in the central nervous system (CNS). Vascular endothelial growth factor-C (VEGF-C) regulates MLV development and maintenance and has therapeutic potential for treating neurological disorders. Herein, we investigated the effects of VEGF-C overexpression on brain fluid drainage and ischemic stroke outcomes in mice. Intracerebrospinal administration of an adeno-associated virus expressing mouse full-length VEGF-C (AAV-mVEGF-C) increased CSF drainage to the deep cervical lymph nodes (dCLNs) by enhancing lymphatic growth and upregulated neuroprotective signaling pathways identified by single nuclei RNA sequencing of brain cells. In a mouse model of ischemic stroke, AAV-mVEGF-C pretreatment reduced stroke injury and ameliorated motor performances in the subacute stage, associated with mitigated microglia-mediated inflammation and increased BDNF signaling in brain cells. Neuroprotective effects of VEGF-C were lost upon cauterization of the dCLN afferent lymphatics and not mimicked by acute post-stroke VEGF-C injection. We conclude that VEGF-C prophylaxis promotes multiple vascular, immune, and neural responses that culminate in a protection against neurological damage in acute ischemic stroke.

Transient Impairment in Microglial Function Causes Sex-Specific Deficits in Synaptic and Hippocampal Function in Mice Exposed to Early Adversity

Abnormal development and function of the hippocampus are two of the most consistent findings in humans and rodents exposed to early life adversity, with males often being more affected than females. Using the limited bedding (LB) paradigm as a rodent model of early life adversity, we found that male adolescent mice that had been exposed to LB exhibit significant deficits in contextual fear conditioning and synaptic connectivity in the hippocampus, which are not observed in females. This is linked to altered developmental refinement of connectivity, with LB severely impairing microglial-mediated synaptic pruning in the hippocampus of male and female pups on postnatal day 17 (P17), but not in adolescent P33 mice when levels of synaptic engulfment by microglia are substantially lower. Since the hippocampus undergoes intense synaptic pruning during the second and third weeks of life, we investigated whether microglia are required for the synaptic and behavioral aberrations observed in adolescent LB mice. Indeed, transient ablation of microglia from P13-21, in normally developing mice caused sex-specific behavioral and synaptic abnormalities similar to those observed in adolescent LB mice. Furthermore, chemogenetic activation of microglia during the same period reversed the microglial-mediated phagocytic deficits at P17 and restored normal contextual fear conditioning and synaptic connectivity in adolescent LB male mice. Our data support an additional contribution of astrocytes in the sex-specific effects of LB, with increased expression of the membrane receptor MEGF10 and enhanced synaptic engulfment in hippocampal astrocytes of 17-day-old LB females, but not in LB male littermates. This finding suggests a potential compensatory mechanism that may explain the relative resilience of LB females. Collectively, these studies highlight a novel role for glial cells in mediating sex-specific hippocampal deficits in a mouse model of early-life adversity.

Cross-hemicord spinal fiber reorganization associates with cortical sensory and motor network expansion in the rat model of hemicontusion cervical spinal cord injury

Magnetic resonance imaging plays an important role in characterizing microstructural changes and reorganization after traumatic injuries to the nervous system. In this study, we tested the feasibility of ex-vivo spinal cord diffusion tensor imaging (DTI) in combination with in vivo brain functional MRI to characterize spinal reorganization and its supraspinal association after a hemicontusion cervical spinal cord injury (SCI). DTI parameters (fractional anisotropy [FA], mean diffusion [MD]) and fiber orientation changes related to reorganization in the contused cervical spinal cord were compared to sham specimens. Altered fiber density and fiber directions occurred across the ipsilateral and contralateral hemicords but with only ipsilateral FA and MD changes. The hemicontusion SCI resulted in ipsilateral fiber breaks, voids and vivid fiber reorientations along the injury epicenter. Fiber directional changes below the injury level were primarily inter-hemispheric, indicating prominent below-level cross-hemispheric reorganization. In vivo resting state functional connectivity of the brain from the respective rats before obtaining the spinal cord samples indicated spatial expansion and increased connectivity strength across both the sensory and motor networks after SCI. The consistency of the neuroplastic changes along the neuraxis (both brain and spinal cord) at the single-subject level, indicates that distinctive reorganizational relationships exist between the spinal cord and the brain post-SCI.

Fluorescent liposomal nanocarriers for targeted drug delivery in ischemic stroke therapy

Ischemic stroke causes acute CNS injury and long-term disability, with limited treatment options such as surgical clot removal or clot-busting drugs. Neuroprotective therapies are needed to protect vulnerable brain regions. The purinergic receptor P2X4 is activated during stroke and exacerbates post-stroke damage. The chemical compound 5-(3-Bromophenyl)-1,3-dihydro-2H-Benzofuro[3,2-e]-1,4-diazepin-2-one (5BDBD) inhibits P2X4 and has shown neuroprotective effects in rodents. However, it is difficult to formulate for systemic delivery to the CNS. The current manuscript reports for the first time, the synthesis and characterization of 5BDBD PEGylated liposomal formulations and evaluates their feasibility to treat stroke in a preclinical mice model. A PEGylated liposomal formulation of 5BDBD was synthesized and characterized, with encapsulation efficacy of >80%, and release over 48 hours. In vitro and in vivo experiments with Nile red encapsulation showed cytocompatibility and CNS infiltration of nanocarriers. Administered 4 or 28 hours after stroke onset, the nanoformulation provided significant neuroprotection, reducing infarct volume by ∼50% compared to controls. It outperformed orally-administered 5BDBD with a lower dose and shorter treatment duration, suggesting precise delivery by nanoformulation improves outcomes. The fluorescent nanoformulations may serve as a platform for delivering and tracking therapeutic agents for stroke treatment.

Spatiotemporal features of neurovascular (un)coupling with stimulus-induced activity and hypercapnia challenge in cerebral cortex and olfactory bulb

Carbon dioxide (CO2) is traditionally considered as metabolic waste, yet its regulation is critical for brain function. It is well accepted that hypercapnia initiates vasodilation, but its effect on neuronal activity is less clear. Distinguishing how stimulus- and CO2-induced vasodilatory responses are (dis)associated with neuronal activity has profound clinical and experimental relevance. We used an optical method in mice to simultaneously image fluorescent calcium (Ca2+) transients from neurons and reflectometric hemodynamic signals during brief sensory stimuli (i.e., hindpaw, odor) and CO2 exposure (i.e., 5%). Stimuli-induced neuronal and hemodynamic responses swiftly increased within locally activated regions exhibiting robust neurovascular coupling. However, hypercapnia produced slower global vasodilation which was temporally uncoupled to neuronal deactivation. With trends consistent across cerebral cortex and olfactory bulb as well as data from GCaMP6f/jRGECO1a mice (i.e., green/red Ca2+ fluorescence), these results unequivocally reveal that stimuli and CO2 generate comparable vasodilatory responses but contrasting neuronal responses. In summary, observations of stimuli-induced regional neurovascular coupling and CO2-induced global neurovascular uncoupling call for careful appraisal when using CO2 in gas mixtures to affect vascular tone and/or neuronal excitability, because CO2 is both a potent vasomodulator and a neuromodulator.

A multi-laboratory preclinical trial in rodents to assess treatment candidates for acute ischemic stroke

Human diseases may be modeled in animals to allow preclinical assessment of putative new clinical interventions. Recent, highly publicized failures of large clinical trials called into question the rigor, design, and value of preclinical assessment. We established the Stroke Preclinical Assessment Network (SPAN) to design and implement a randomized, controlled, blinded, multi-laboratory trial for the rigorous assessment of candidate stroke treatments combined with intravascular thrombectomy. Efficacy and futility boundaries in a multi-arm multi-stage statistical design aimed to exclude from further study highly effective or futile interventions after each of four sequential stages. Six independent research laboratories performed a standard focal cerebral ischemic insult in five animal models that included equal numbers of males and females: young mice, young rats, aging mice, mice with diet-induced obesity, and spontaneously hypertensive rats. The laboratories adhered to a common protocol and efficiently enrolled 2615 animals with full data completion and comprehensive animal tracking. SPAN successfully implemented treatment masking, randomization, prerandomization inclusion and exclusion criteria, and blinded assessment of outcomes. The SPAN design and infrastructure provide an effective approach that could be used in similar preclinical, multi-laboratory studies in other disease areas and should help improve reproducibility in translational science.

VEGF-C promotes brain-derived fluid drainage, confers neuroprotection, and improves stroke outcomes

Meningeal lymphatic vessels promote tissue clearance and immune surveillance in the central nervous system (CNS). Vascular endothelium growth factor-C (VEGF-C) is essential for meningeal lymphatic development and maintenance and has therapeutic potential for treating neurological disorders, including ischemic stroke. We have investigated the effects of VEGF-C overexpression on brain fluid drainage, single cell transcriptome in the brain, and stroke outcomes in adult mice. Intra-cerebrospinal fluid administration of an adeno-associated virus expressing VEGF-C (AAV-VEGF-C) increases the CNS lymphatic network. Post-contrast T1 mapping of the head and neck showed that deep cervical lymph node size and drainage of CNS-derived fluids were increased. Single nuclei RNA sequencing revealed a neuro-supportive role of VEGF-C via upregulation of calcium and brain-derived neurotrophic factor (BDNF) signaling pathways in brain cells. In a mouse model of ischemic stroke, AAV-VEGF-C pretreatment reduced stroke injury and ameliorated motor performances in the subacute stage. AAV-VEGF-C thus promotes CNS-derived fluid and solute drainage, confers neuroprotection, and reduces ischemic stroke damage.

Short abstract

Intrathecal delivery of VEGF-C increases the lymphatic drainage of brain-derived fluids confers neuroprotection, and improves neurological outcomes after ischemic stroke.

A consensus protocol for functional connectivity analysis in the rat brain

Task-free functional connectivity in animal models provides an experimental framework to examine connectivity phenomena under controlled conditions and allows for comparisons with data modalities collected under invasive or terminal procedures. Currently, animal acquisitions are performed with varying protocols and analyses that hamper result comparison and integration. Here we introduce StandardRat, a consensus rat functional magnetic resonance imaging acquisition protocol tested across 20 centers. To develop this protocol with optimized acquisition and processing parameters, we initially aggregated 65 functional imaging datasets acquired from rats across 46 centers. We developed a reproducible pipeline for analyzing rat data acquired with diverse protocols and determined experimental and processing parameters associated with the robust detection of functional connectivity across centers. We show that the standardized protocol enhances biologically plausible functional connectivity patterns relative to previous acquisitions. The protocol and processing pipeline described here is openly shared with the neuroimaging community to promote interoperability and cooperation toward tackling the most important challenges in neuroscience.

Concerted roles of LRRTM1 and SynCAM 1 in organizing prefrontal cortex synapses and cognitive functions

Multiple trans-synaptic complexes organize synapse development, yet their roles in the mature brain and cooperation remain unclear. We analyzed the postsynaptic adhesion protein LRRTM1 in the prefrontal cortex (PFC), a region relevant to cognition and disorders. LRRTM1 knockout (KO) mice had fewer synapses, and we asked whether other synapse organizers counteract further loss. This determined that the immunoglobulin family member SynCAM 1 controls synapse number in PFC and was upregulated upon LRRTM1 loss. Combined LRRTM1 and SynCAM 1 deletion substantially lowered dendritic spine number in PFC, but not hippocampus, more than the sum of single KO impairments. Their cooperation extended presynaptically, and puncta of Neurexins, LRRTM1 partners, were less abundant in double KO (DKO) PFC. Electrophysiology and fMRI demonstrated aberrant neuronal activity in DKO mice. Further, DKO mice were impaired in social interactions and cognitive tasks. Our results reveal concerted roles of LRRTM1 and SynCAM 1 across synaptic, network, and behavioral domains.

Decreased but diverse activity of cortical and thalamic neurons in consciousness-impairing rodent absence seizures

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.

Lateralized Supraspinal Functional Connectivity Correlate with Pain and Motor Dysfunction in Rat Hemicontusion Cervical Spinal Cord Injury

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.

Thalamic activations in rat brain by fMRI during tactile (forepaw, whisker) and non-tactile (visual, olfactory) sensory stimulations

The thalamus is a crucial subcortical hub that impacts cortical activity. Tracing experiments in animals and post-mortem humans suggest rich morphological specificity of the thalamus. Very few studies reported rodent thalamic activations by functional MRI (fMRI) as compared to cortical activations for different sensory stimuli. Here, we show different portions of the rat thalamus in response to tactile (forepaw, whisker) and non-tactile (visual, olfactory) sensory stimuli with high field fMRI (11.7T) using a custom-build quadrature surface coil to capture high sensitivity signals from superficial and deep brain regions simultaneously. Results demonstrate reproducible thalamic activations during both tactile and non-tactile stimuli. Forepaw and whisker stimuli activated broader regions within the thalamus: ventral posterior lateral (VPL), ventral posterior medial (VPM), lateral posterior mediorostral (LPMR) and posterior medial (POm) thalamic nuclei. Visual stimuli activated dorsal lateral geniculate nucleus (DLG) of the thalamus but also parts of the superior/inferior colliculus, whereas olfactory stimuli activated specifically the mediodorsal nucleus of the thalamus (MDT). BOLD activations in LGN and MDT were much stronger than in VPL, VPM, LPMR and POm. These fMRI-based thalamic activations suggest that forepaw and whisker (i.e., tactile) stimuli engage VPL, VPM, LPMR and POm whereas visual and olfactory (i.e., non-tactile) stimuli, respectively, recruit DLG and MDT exclusively.

The Stroke Preclinical Assessment Network: Rationale, Design, Feasibility, and Stage 1 Results

Cerebral ischemia and reperfusion initiate cellular events in brain that lead to neurological disability. Investigating these cellular events provides ample targets for developing new treatments. Despite considerable work, no such therapy has translated into successful stroke treatment. Among other issues-such as incomplete mechanistic knowledge and faulty clinical trial design-a key contributor to prior translational failures may be insufficient scientific rigor during preclinical assessment: nonblinded outcome assessment; missing randomization; inappropriate sample sizes; and preclinical assessments in young male animals that ignore relevant biological variables, such as age, sex, and relevant comorbid diseases. Promising results are rarely replicated in multiple laboratories. We sought to address some of these issues with rigorous assessment of candidate treatments across 6 independent research laboratories. The Stroke Preclinical Assessment Network (SPAN) implements state-of-the-art experimental design to test the hypothesis that rigorous preclinical assessment can successfully reduce or eliminate common sources of bias in choosing treatments for evaluation in clinical studies. SPAN is a randomized, placebo-controlled, blinded, multilaboratory trial using a multi-arm multi-stage protocol to select one or more putative stroke treatments with an implied high likelihood of success in human clinical stroke trials. The first stage of SPAN implemented procedural standardization and experimental rigor. All participating research laboratories performed middle cerebral artery occlusion surgery adhering to a common protocol and rapidly enrolled 913 mice in the first of 4 planned stages with excellent protocol adherence, remarkable data completion and low rates of subject loss. SPAN stage 1 successfully implemented treatment masking, randomization, prerandomization inclusion/exclusion criteria, and blinded assessment to exclude bias. Our data suggest that a large, multilaboratory, preclinical assessment effort to reduce known sources of bias is feasible and practical. Subsequent SPAN stages will evaluate candidate treatments for potential success in future stroke clinical trials using aged animals and animals with comorbid conditions.

Supraspinal sensorimotor and pain-related reorganization after a hemicontusion rat cervical spinal cord injury

Since the presence of pain impedes motor recovery in individuals with spinal cord injury (SCI), it is necessary to understand their supraspinal substrates in translational animal models. Using functional magnetic resonance imaging (fMRI) in a rat model of hemicontusion cervical SCI, supraspinal changes were mapped and correlated with sensorimotor behavioral outcomes. Female adult rats underwent sham or SCI using a 2.5 mm impactor and 150 kDyne force. SCI permanently impaired motor activity in only the ipsilesional forelimb along with thermal hyperalgesia at 5 and 6 wks. Spinal MRI at 8 wks after SCI showed ipsilateral T1 and T2 lesions with no discernable lesions across shams. fMRI mapping during electrical forepaw stimulation indicated SCI-induced sensorimotor reorganization with an expansion of the contralesional forelimb representation. Resting state fMRI based functional connectivity density (FCD), a marker of regional neuronal hubs increased or decreased across brain regions involved in nociception. FCD increases after SCI were in the primary and secondary somatosensory cortices (S1 and S2), anterior cingulate cortex (ACC), insula and the prefrontal cortex (PFC) and decreases were across the hippocampus, thalamus, hypothalamus and amygdala in SCI. Resting state functional connectivity (RSFC) assessments from the FCD altered regions of interest indicated cortico-cortical RSFC increases and cortico-insular, cortico-thalamic and cortico-hypothalamic RSFC decreases after SCI. Hippocampus, amygdala and thalamus showed decreased RSFC with most cortical regions and between themselves except the hippocampus-amygdala network, which showed increased RSFC after SCI. While select nociceptive region's intrinsic activity associated strongly with evoked pain behaviors after SCI (eg., PFC, ACC, hippocampus, thalamus, hypothalamus, M1 and S1BF) other nociceptive regions had weaker associations (eg., amygdala, insula, auditory cortex, S1FL, S1HL, S2 and M2), but differed significantly in their intrinsic activities between sham and SCI. The weaker associated nociceptive regions may possibly encode both the evoked and affective components of pain.

Small loci of astroglial glutamine synthetase deficiency in the postnatal brain cause epileptic seizures and impaired functional connectivity

OBJECTIVE

The astroglial enzyme glutamine synthetase (GS) is deficient in small loci in the brain in adult patients with different types of focal epilepsy; however, the role of this deficiency in the pathogenesis of epilepsy has been difficult to assess due to a lack of sufficiently sensitive and specific animal models. The aim of this study was to develop an in vivo approach for precise and specific deletions of the GS gene in the postnatal brain.

METHODS

We stereotaxically injected various adeno-associated virus (AAV)-Cre recombinase constructs into the hippocampal formation and neocortex in 22-70-week-old GS^flox/flox mice to knock out the GS gene in a specific and focal manner. The mice were subjected to seizure threshold determination, continuous video-electroencephalographic recordings, advanced in vivo neuroimaging, and immunocytochemistry for GS.

RESULTS

The construct AAV8-glial fibrillary acidic protein-green fluorescent protein-Cre eliminated GS in >99% of astrocytes in the injection center with a gradual return to full GS expression toward the periphery. Such focal GS deletion reduced seizure threshold, caused spontaneous recurrent seizures, and diminished functional connectivity.

SIGNIFICANCE

These results suggest that small loci of GS deficiency in the postnatal brain are sufficient to cause epilepsy and impaired functional connectivity. Additionally, given the high specificity and precise spatial resolution of our GS knockdown approach, we anticipate that this model will be extremely useful for rigorous in vivo and ex vivo studies of astroglial GS function at the brain-region and single-cell levels.

Association Between Magnetic Resonance Imaging-Based Spinal Morphometry and Sensorimotor Behavior in a Hemicontusion Model of Incomplete Cervical Spinal Cord Injury in Rats

Aim: Structural connectivity in the reorganizing spinal cord after injury dictates functional connectivity and hence the neurological outcome. As magnetic resonance imaging (MRI)-based structural parameters are mostly accessible across spinal cord injury (SCI) patients, we studied MRI-based spinal morphological changes and their relationship to neurological outcome in the rat model of cervical SCI. Introduction: Functional connectivity assessments on patients with SCI rely heavily on MRI-based approaches to investigate the complete neural axis (both spinal cord and brain). Hence, underlying MRI-based structural and morphometric changes in the reorganizing spinal cord and their relationship to neurological outcomes is crucial for meaningful interpretation of functional connectivity changes across the neural axis. Methods: Young adult rats, aged 1.5 months, underwent a precise mechanical impact hemicontusion incomplete cervical SCI at the C4/C5 level, after which sensorimotor behavioral assessments were tracked during the reorganization period of 1-6 weeks, followed by MRI of the cervical spinal cord at 8 weeks after SCI. Results: A significant ipsilesional forelimb motor debilitation was observed from 1 to 6 weeks after injury. Heat sensitivity testing (Hargreaves) showed ipsilesional forelimb hypersensitivity at 5 and 6 weeks after SCI. MRI of the cervical spine showed ipsilateral T1- and T2-weighted lesions across all SCI rats compared with no significant lesions in sham rats. Morphometric assessments of the lesional and nonlesional changes showed the diverse nature of their interindividual variability in the SCI receiving rats. While the various T1 and T2 MRI lesional volumes associated weakly or moderately with neurological outcome, the nonlesional spinal morphometric changes associated much more strongly. The results have important implications for interpreting functional MRI-based functional connectivity after SCI by providing vital underlying structural changes and their relative neurological impact. Impact statement Functional connectivity assessments on patients with SCI relies heavily upon MRI based approaches. Hence, underlying MRI based structural and morphometric changes in the reorganizing spinal cord and its relationship to neurological outcomes is vital for meaningful interpretation of functional connectivity changes across the complete neural axis (both spinal cord and the brain).

APOE genotype-dependent pharmacogenetic responses to rapamycin for preventing Alzheimer's disease

The ε4 allele of Apolipoprotein (APOE4) is the strongest genetic risk factor for Alzheimer's disease (AD), the most common form of dementia. Cognitively normal APOE4 carriers have developed amyloid beta (Aβ) plaques and cerebrovascular, metabolic and structural deficits decades before showing the cognitive impairment. Interventions that can inhibit Aβ retention and restore the brain functions to normal would be critical to prevent AD for the asymptomatic APOE4 carriers. A major goal of the study was to identify the potential usefulness of rapamycin (Rapa), a pharmacological intervention for extending longevity, for preventing AD in the mice that express human APOE4 gene and overexpress Aβ (the E4FAD mice). Another goal of the study was to identify the potential pharmacogenetic differences in response to rapamycin between the E4FAD and E3FAD mice, the mice with human APOE ε3 allele. We used multi-modal MRI to measure in vivo cerebral blood flow (CBF), neurotransmitter levels, white matter integrity, water content, cerebrovascular reactivity (CVR) and somatosensory response; used behavioral assessments to determine cognitive function; used biochemistry assays to determine Aβ retention and blood-brain barrier (BBB) functions; and used metabolomics to identify brain metabolic changes. We found that in the E4FAD mice, rapamycin normalized bodyweight, restored CBF (especially in female), BBB activity for Aβ transport, neurotransmitter levels, neuronal integrity and free fatty acid level, and reduced Aβ retention, which were not observe in the E3FAD-Rapa mice. In contrast, E3FAD-Rapa mice had lower CVR responses, lower anxiety and reduced glycolysis in the brain, which were not seen in the E4FAD-Rapa mice. Further, rapamycin appeared to normalize lipid-associated metabolism in the E4FAD mice, while slowed overall glucose-associated metabolism in the E3FAD mice. Finally, rapamycin enhanced overall water content, water diffusion in white matter, and spatial memory in both E3FAD and E4FAD mice, but did not impact the somatosensory responses under hindpaw stimulation. Our findings indicated that rapamycin was able to restore brain functions and reduce AD risk for young, asymptomatic E4FAD mice, and there were pharmacogenetic differences between the E3FAD and E4FAD mice. As the multi-modal MRI methods used in the study are readily to be used in humans and rapamycin is FDA-approved, our results may pave a way for future clinical testing of the pharmacogenetic responses in humans with different APOE alleles, and potentially using rapamycin to prevent AD for asymptomatic APOE4 carriers.

Kaempferol Treatment after Traumatic Brain Injury during Early Development Mitigates Brain Parenchymal Microstructure and Neural Functional Connectivity Deterioration at Adolescence

Targeting mitochondrial ion homeostasis using Kaempferol, a mitochondrial Ca2+ uniporter channel activator, improves energy metabolism and behavior soon after a traumatic brain injury (TBI) in developing rats. Because of broad TBI pathophysiology and brain mitochondrial heterogeneity, Kaempferol-mediated early-stage behavioral and brain metabolic benefits may accrue from diverse sources within the brain. We hypothesized that Kaempferol influences TBI outcome by differentially impacting the neural, vascular, and synaptic/axonal compartments. After TBI at early development (P31), functional magnetic resonance imaging and diffusion tensor imaging (DTI) were applied to determine imaging outcomes at adolescence (2 months post-injury). Vehicle and Kaempferol treatments were made at 1, 24, and 48 h post-TBI, and their effects were assessed at adolescence. A significant increase in neural connectivity was observed after Kaempferol treatment as assessed by the spatial extent and strength of the somatosensory cortical and hippocampal resting-state functional connectivity (RSFC) networks. However, no significant RSFC changes were observed in the thalamus. DTI measures of fractional anisotropy (FA) and apparent diffusion coefficient, representing synaptic/axonal and microstructural integrity, showed significant improvements after Kaempferol treatment, with highest changes in the frontal and parietal cortices and hippocampus. Kaempferol treatment also increased corpus callosal FA, indicating measurable improvement in the interhemispheric structural connectivity. TBI prognosis was significantly altered at adolescence by early Kaempferol treatment, with improved neural connectivity, neurovascular coupling, and parenchymal microstructure in select brain regions. However, Kaempferol failed to improve vasomotive function across the whole brain, as measured by cerebrovascular reactivity. The differential effects of Kaempferol treatment on various brain functional compartments support diverse cellular-level mitochondrial functional outcomes in vivo.

Orthonasal versus retronasal glomerular activity in rat olfactory bulb by fMRI

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.

Alterations of Parenchymal Microstructure, Neuronal Connectivity, and Cerebrovascular Resistance at Adolescence after Mild-to-Moderate Traumatic Brain Injury in Early Development

Traumatic brain injury (TBI) is a leading cause of morbidity in children. To investigate outcome of early developmental TBI during adolescence, a rat model of fluid percussion injury was developed, where previous work reported deficits in sensorimotor behavior and cortical blood flow at adolescence.1 Based on the nonlocalized outcome, we hypothesized that multiple neurophysiological components of brain function, namely neuronal connectivity, synapse/axonal microstructural integrity, and neurovascular function, are altered and magnetic resonance imaging (MRI) methods could be used to determine regional alterations. Adolescent outcomes of developmental TBI were studied 2 months after injury, using functional MRI (fMRI) and diffusion tensor imaging (DTI). fMRI-based resting-state functional connectivity (RSFC), representing neural connectivity, was significantly altered between sham and TBI. RSFC strength decreased in the cortex, hippocampus, and thalamus, accompanied by decrease in spatial extent of their corresponding RSFC networks and interhemispheric asymmetry. Cerebrovascular reactivity to arterial CO2 changes diminished after TBI across both hemispheres, with a more pronounced decrease in the ipsilateral hippocampus, thalamus, and motor cortex. DTI measures of fractional anisotropy and apparent diffusion coefficient, reporting on axonal and microstructural integrity of the brain, indicated similar interhemispheric asymmetry, with highest change in the ipsilateral hippocampus and regions adjoining the ipsilateral thalamus, hypothalamus, and amygdala. TBI-induced corpus callosal microstructural alterations indicated measurable changes in interhemispheric structural connectivity. Hippocampus, thalamus, and select cortical regions were most consistently affected in multiple imaging markers. The multi-modal MRI results demonstrate cortical and subcortical alterations in neural connectivity, cerebrovascular resistance, and parenchymal microstructure in the adolescent brain, indicating the highly diffuse and persistent nature of the lateral fluid percussion TBI early in development.

Spontaneous activity forms a foundation for odor-evoked activation maps in the rat olfactory bulb

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.

Metabolic demands of neural-hemodynamic associated and disassociated areas in brain

Interpretation of regional blood oxygenation level-dependent (BOLD) responses in functional magnetic resonance imaging (fMRI) is contingent on whether local field potential (LFP) and multi-unit activity (MUA) is either dissociated or associated. To examine whether neural-hemodynamic associated and dissociated areas have different metabolic demands, we recorded sensory-evoked responses of BOLD signal, blood flow (CBF), and blood volume (CBV), which with calibrated fMRI provided oxidative metabolism (CMRO2) from rat's ventral posterolateral thalamic nucleus (VPL) and somatosensory forelimb cortex (S1FL) and compared these neuroimaging signals to neurophysiological recordings. MUA faithfully recorded evoked latency differences between VPL and S1FL because evoked MUA in these regions were similar in magnitude. Since evoked LFP was significantly attenuated in VPL, we extracted the time courses of the weaker thalamic LFP to compare with the stronger cortical LFP using wavelet transform. BOLD and CBV responses were greater in S1FL than in VPL, similar to LFP regional differences. CBF and CMRO2 responses were both comparably larger in S1FL and VPL. Despite different levels of CBF-CMRO2 and LFP-MUA couplings in VPL and S1FL, the CMRO2 was well matched with MUA in both regions. These results suggest that neural-hemodynamic associated and dissociated areas in VPL and S1FL can have similar metabolic demands.

New horizons in neurometabolic and neurovascular coupling from calibrated fMRI

Neurovascular coupling relates changes in neuronal activity to constriction/dilation of microvessels. However neurometabolic coupling, which is less well known, relates alterations in neuronal activity with metabolic demands. The link between the blood oxygenation level dependent (BOLD) signal and neural activity opened doors for functional MRI (fMRI) to be a powerful neuroimaging tool in the neurosciences. But due to the complex makeup of BOLD contrast, researchers began to investigate the relationship between BOLD signal and blood flow and/or volume changes during functional brain activation, which together provided the tools to measure oxygen consumption on the basis of the biophysical model of BOLD. This field is called calibrated fMRI, thereby allowed probing of both neurometabolic and neurovascular couplings for a variety of health conditions in animals and humans. Calibrated fMRI may provide brain disorder biomarkers that could be used for monitoring effective therapies.

Afatinib plus Cetuximab Delays Resistance Compared to Single-Agent Erlotinib or Afatinib in Mouse Models of TKI-Naïve EGFR L858R-Induced Lung Adenocarcinoma

PURPOSE

The EGFR tyrosine kinase inhibitors (TKIs), erlotinib and afatinib, have transformed the treatment of advanced EGFR-mutant lung adenocarcinoma. However, almost all patients who respond develop acquired resistance on average approximately 1 year after starting therapy. Resistance is commonly due to a secondary mutation in EGFR (EGFR(T790M)). We previously found that the combination of the EGFR TKI afatinib and the EGFR antibody cetuximab could overcome EGFR(T790M)-mediated resistance in preclinical models. This combination has shown a 29% response rate in a clinical trial in patients with acquired resistance to first-generation TKIs. An outstanding question is whether this regimen is beneficial when used as first-line therapy.

EXPERIMENTAL DESIGN

Using mouse models of EGFR-mutant lung cancer, we tested whether the combination of afatinib plus cetuximab delivered upfront to mice with TKI-naïve EGFR(L858R)-induced lung adenocarcinomas delayed tumor relapse and drug-resistance compared with single-agent TKIs.

RESULTS

Afatinib plus cetuximab markedly delayed the time to relapse and incidence of drug-resistant tumors, which occurred in only 63.6% of the mice, in contrast to erlotinib or afatinib treatment where 100% of mice developed resistance. Mechanisms of tumor escape observed in afatinib plus cetuximab resistant tumors include the EGFR(T790M) mutation and Kras mutations. Experiments in cell lines and xenografts confirmed that the afatinib plus cetuximab combination does not suppress the emergence of EGFR(T790M).

CONCLUSIONS

These results highlight the potential of afatinib plus cetuximab as an effective treatment strategy for patients with TKI-naïve EGFR-mutant lung cancer and indicate that clinical trial development in this area is warranted.

Comparison of glomerular activity patterns by fMRI and wide-field calcium imaging: Implications for principles underlying odor mapping

Functional imaging signals arise from distinct metabolic and hemodynamic events at the neuropil, but how these processes are influenced by pre- and post-synaptic activities need to be understood for quantitative interpretation of stimulus-evoked mapping data. The olfactory bulb (OB) glomeruli, spherical neuropil regions with well-defined neuronal circuitry, can provide insights into this issue. Optical calcium-sensitive fluorescent dye imaging (OICa(2+)) reflects dynamics of pre-synaptic input to glomeruli, whereas high-resolution functional magnetic resonance imaging (fMRI) using deoxyhemoglobin contrast reveals neuropil function within the glomerular layer where both pre- and post-synaptic activities contribute. We imaged odor-specific activity patterns of the dorsal OB in the same anesthetized rats with fMRI and OICa(2+) and then co-registered the respective maps to compare patterns in the same space. Maps by each modality were very reproducible as trial-to-trial patterns for a given odor, overlapping by ~80%. Maps evoked by ethyl butyrate and methyl valerate for a given modality overlapped by ~80%, suggesting activation of similar dorsal glomerular networks by these odors. Comparison of maps generated by both methods for a given odor showed ~70% overlap, indicating similar odor-specific maps by each method. These results suggest that odor-specific glomerular patterns by high-resolution fMRI primarily tracks pre-synaptic input to the OB. Thus combining OICa(2+) and fMRI lays the framework for studies of OB processing over a range of spatiotemporal scales, where OICa(2+) can feature the fast dynamics of dorsal glomerular clusters and fMRI can map the entire glomerular sheet in the OB.

Quantitative β mapping for calibrated fMRI

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.

Role of mitochondrial calcium uptake homeostasis in resting state fMRI brain networks

Mitochondrial Ca(2+) uptake influences both brain energy metabolism and neural signaling. Given that brain mitochondrial organelles are distributed in relation to vascular density, which varies considerably across brain regions, we hypothesized different physiological impacts of mitochondrial Ca(2+) uptake across brain regions. We tested the hypothesis by monitoring brain "intrinsic activity" derived from the resting state functional MRI (fMRI) blood oxygen level dependent (BOLD) fluctuations in different functional networks spanning the somatosensory cortex, caudate putamen, hippocampus and thalamus, in normal and perturbed mitochondrial Ca(2+) uptake states. In anesthetized rats at 11.7 T, mitochondrial Ca(2+) uptake was inhibited or enhanced respectively by treatments with Ru360 or kaempferol. Surprisingly, mitochondrial Ca(2+) uptake inhibition by Ru360 and enhancement by kaempferol led to similar dose-dependent decreases in brain-wide intrinsic activities in both the frequency domain (spectral amplitude) and temporal domain (resting state functional connectivity; RSFC). The fact that there were similar dose-dependent decreases in the frequency and temporal domains of the resting state fMRI-BOLD fluctuations during mitochondrial Ca(2+) uptake inhibition or enhancement indicated that mitochondrial Ca(2+) uptake and its homeostasis may strongly influence the brain's functional organization at rest. Interestingly, the resting state fMRI-derived intrinsic activities in the caudate putamen and thalamic regions saturated much faster with increasing dosage of either drug treatment than the drug-induced trends observed in cortical and hippocampal regions. Regional differences in how the spectral amplitude and RSFC changed with treatment indicate distinct mitochondrion-mediated spontaneous neuronal activity coupling within the various RSFC networks determined by resting state fMRI.

Brain region and activity-dependent properties of M for calibrated fMRI

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.

Distribution of temperature changes and neurovascular coupling in rat brain following 3,4-methylenedioxymethamphetamine (MDMA, "ecstasy") exposure

(+/-)3,4-methylenedioxymethamphetamine (MDMA, "ecstasy") is an abused psychostimulant that produces strong monoaminergic stimulation and whole-body hyperthermia. MDMA-induced thermogenesis involves activation of uncoupling proteins (UCPs), primarily a type specific to skeletal muscle (UCP-3) and absent from the brain, although other UCP types are expressed in the brain (e.g. thalamus) and might contribute to thermogenesis. Since neuroimaging of brain temperature could provide insights into MDMA action, we measured spatial distributions of systemically administered MDMA-induced temperature changes and dynamics in rat cortex and subcortex using a novel magnetic resonance method, Biosensor Imaging of Redundant Deviation in Shifts (BIRDS), with an exogenous temperature-sensitive probe (thulium ion and macrocyclic chelate 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetramethyl-1,4,7,10-tetraacetate (DOTMA(4-))). The MDMA-induced temperature rise was greater in the cortex than in the subcortex (1.6 ± 0.4 °C versus 1.3 ± 0.4 °C) and occurred more rapidly (2.0 ± 0.2 °C/h versus 1.5 ± 0.2 °C/h). MDMA-induced temperature changes and dynamics in the cortex and body were correlated, although the body temperature exceeded the cortex temperature before and after MDMA. Temperature, neuronal activity, and blood flow (CBF) were measured simultaneously in the cortex and subcortex (i.e. thalamus) to investigate possible differences of MDMA-induced warming across brain regions. MDMA-induced warming correlated with increases in neuronal activity and blood flow in the cortex, suggesting that the normal neurovascular response to increased neural activity was maintained. In contrast to the cortex, a biphasic relationship was seen in the subcortex (i.e. thalamus), with a decline in CBF as temperature and neural activity rose, transitioning to a rise in CBF for temperature above 37 °C, suggesting that MDMA affected CBF and neurovascular coupling differently in subcortical regions. Considering that MDMA effects on CBF and heat dissipation (as well as potential heat generation) may vary regionally, neuroprotection may require different cooling strategies.

Rhythmic 3-4Hz discharge is insufficient to produce cortical BOLD fMRI decreases in generalized seizures

Absence seizures are transient episodes of impaired consciousness accompanied by 3-4 Hz spike-wave discharge on electroencephalography (EEG). Human functional magnetic resonance imaging (fMRI) studies have demonstrated widespread cortical decreases in the blood oxygen-level dependent (BOLD) signal that may play an important role in the pathophysiology of these seizures. Animal models could provide an opportunity to investigate the fundamental mechanisms of these changes, however they have so far failed to consistently replicate the cortical fMRI decreases observed in human patients. This may be due to important differences between human seizures and animal models, including a lack of cortical development in rodents or differences in the frequencies of rodent (7-8 Hz) and human (3-4 Hz) spike-wave discharges. To examine the possible contributions of these differences, we developed a ferret model that exhibits 3-4 Hz spike-wave seizures in the presence of a sulcated cortex. Measurements of BOLD fMRI and simultaneous EEG demonstrated cortical fMRI increases during and following spike-wave seizures in ferrets. However unlike human patients, significant fMRI decreases were not observed. The lack of fMRI decreases was consistent across seizures of different durations, discharge frequencies, and anesthetic regimes, and using fMRI analysis models similar to human patients. In contrast, generalized tonic-clonic seizures under the same conditions elicited sustained postictal fMRI decreases, verifying that the lack of fMRI decreases with spike-wave was not due to technical factors. These findings demonstrate that 3-4 Hz spike-wave discharge in a sulcated animal model does not necessarily produce fMRI decreases, leaving the mechanism for this phenomenon open for further investigation.

DYNAmic Multi-coIl TEchnique (DYNAMITE) shimming of the rat brain at 11.7 T

The in vivo rat model is a workhorse in neuroscience research, preclinical studies and drug development. A repertoire of MR tools has been developed for its investigation; however, high levels of B0 magnetic field homogeneity are required for meaningful results. The homogenization of magnetic fields in the rat brain, i.e. shimming, is a difficult task because of a multitude of complex, susceptibility-induced field distortions. Conventional shimming with spherical harmonic (SH) functions is capable of compensating for shallow field distortions in limited areas, e.g. in the cortex, but performs poorly in difficult-to-shim subcortical structures or for the entire brain. Based on the recently introduced multi-coil approach for magnetic field modeling, the DYNAmic Multi-coIl TEchnique (DYNAMITE) is introduced for magnetic field shimming of the in vivo rat brain and its benefits for gradient-echo echo-planar imaging (EPI) are demonstrated. An integrated multi-coil/radiofrequency (MC/RF) system comprising 48 individual localized DC coils for B0 shimming and a surface transceive RF coil has been developed that allows MR investigations of the anesthetized rat brain in vivo. DYNAMITE shimming with this MC/RF set-up is shown to reduce the B0 standard deviation to a third of that achieved with current shim technology employing static first- through third-order SH shapes. The EPI signal over the rat brain increased by 31%, and a 24% gain in usable EPI voxels could be realized. DYNAMITE shimming is expected to critically benefit a wide range of preclinical and neuroscientific MR research. Improved magnetic field homogeneity, together with the achievable large brain coverage of this method, will be crucial when signal pathways, cortical circuitry or the brain's default network are studied. Together with the efficiency gains of MC-based shimming compared with SH approaches demonstrated recently, DYNAMITE shimming has the potential to replace conventional SH shim systems in small-bore animal scanners.

Decreased resting functional connectivity after traumatic brain injury in the rat

Traumatic brain injury (TBI) contributes to about 10% of acquired epilepsy. Even though the mechanisms of post-traumatic epileptogenesis are poorly known, a disruption of neuronal networks predisposing to altered neuronal synchrony remains a viable candidate mechanism. We tested a hypothesis that resting state BOLD-fMRI functional connectivity can reveal network abnormalities in brain regions that are connected to the lesioned cortex, and that these changes associate with functional impairment, particularly epileptogenesis. TBI was induced using lateral fluid-percussion injury in seven adult male Sprague-Dawley rats followed by functional imaging at 9.4T 4 months later. As controls we used six sham-operated animals that underwent all surgical operations but were not injured. Electroencephalogram (EEG)-functional magnetic resonance imaging (fMRI) was performed to measure resting functional connectivity. A week after functional imaging, rats were implanted with bipolar skull electrodes. After recovery, rats underwent pentyleneterazol (PTZ) seizure-susceptibility test under EEG. For image analysis, four pairs of regions of interests were analyzed in each hemisphere: ipsilateral and contralateral frontal and parietal cortex, hippocampus, and thalamus. High-pass and low-pass filters were applied to functional imaging data. Group statistics comparing injured and sham-operated rats and correlations over time between each region were calculated. In the end, rats were perfused for histology. None of the rats had epileptiform discharges during functional imaging. PTZ-test, however revealed increased seizure susceptibility in injured rats as compared to controls. Group statistics revealed decreased connectivity between the ipsilateral and contralateral parietal cortex and between the parietal cortex and hippocampus on the side of injury as compared to sham-operated animals. Injured animals also had abnormal negative connectivity between the ipsilateral and contralateral parietal cortex and other regions. Our data provide the first evidence on abnormal functional connectivity after experimental TBI assessed with resting state BOLD-fMRI.

Mapping phosphorylation rate of fluoro-deoxy-glucose in rat brain by (19)F chemical shift imaging

(19)F magnetic resonance spectroscopy (MRS) studies of 2-fluoro-2-deoxy-d-glucose (FDG) and 2-fluoro-2-deoxy-d-glucose-6-phosphate (FDG-6P) can be used for directly assessing total glucose metabolism in vivo. To date, (19)F MRS measurements of FDG phosphorylation in the brain have either been achieved ex vivo from extracted tissue or in vivo by unusually long acquisition times. Electrophysiological and functional magnetic resonance imaging (fMRI) measurements indicate that FDG doses up to 500 mg/kg can be tolerated with minimal side effects on cerebral physiology and evoked fMRI-BOLD responses to forepaw stimulation. In halothane-anesthetized rats, we report localized in vivo detection and separation of FDG and FDG-6P MRS signals with (19)F 2D chemical shift imaging (CSI) at 11.7 T. A metabolic model based on reversible transport between plasma and brain tissue, which included a non-saturable plasma to tissue component, was used to calculate spatial distribution of FDG and FDG-6P concentrations in rat brain. In addition, spatial distribution of rate constants and metabolic fluxes of FDG to FDG-6P conversion were estimated. Mapping the rate of FDG to FDG-6P conversion by (19)F CSI provides an MR methodology that could impact other in vivo applications such as characterization of tumor pathophysiology.

Increased resting functional connectivity in spike-wave epilepsy in WAG/Rij rats

PURPOSE

Functional magnetic resonance imaging (fMRI)-based resting functional connectivity is well suited for measuring slow correlated activity throughout brain networks. Epilepsy involves chronic changes in normal brain networks, and recent work demonstrated enhanced resting fMRI connectivity between the hemispheres in childhood absence epilepsy. An animal model of this phenomenon would be valuable for investigating fundamental mechanisms and testing therapeutic interventions.

METHODS

We used fMRI-based resting functional connectivity for studying brain networks involved in absence epilepsy. Wistar Albino Glaxo rats from Rijswijk (WAG/Rij) exhibit spontaneous episodes of staring and unresponsiveness accompanied by spike-wave discharges (SWDs) resembling human absence seizures in behavior and electroencephalography (EEG). Simultaneous EEG-fMRI data in epileptic WAG/Rij rats in comparison to nonepileptic Wistar controls were acquired at 9.4 T. Regions showing cortical fMRI increases during SWDs were used to define reference regions for connectivity analysis to investigate whether chronic seizure activity is associated with changes in network resting functional connectivity.

KEY FINDINGS

We observed high degrees of cortical-cortical correlations in all WAG/Rij rats at rest (when no SWDs were present), but not in nonepileptic controls. Strongest connectivity was seen between regions most intensely involved in seizures, mainly in the bilateral somatosensory and adjacent cortices. Group statistics revealed that resting interhemispheric cortical-cortical correlations were significantly higher in WAG/Rij rats compared to nonepileptic controls.

SIGNIFICANCE

These findings suggest that activity-dependent plasticity may lead to long-term changes in epileptic networks even at rest. The results show a marked difference between the epileptic and nonepileptic animals in cortical-cortical connectivity, indicating that this may be a useful interictal biomarker associated with the epileptic state.

Functional MRI and neural responses in a rat model of Alzheimer's 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 functional state impacts spontaneous neocortical activity and resting state FMRI

Mitochondrial Ca²⁺ 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²⁺ uptake capacities on intrinsic brain activity using clinically relevant markers. From a multi-compartment estimate of probable baseline Ca²⁺ ranges in the brain, we hypothesized that reduced or enhanced mitochondrial Ca²⁺ 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²⁺ 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²⁺ uptake capacity was perturbed by pharmacologically inhibiting or enhancing the mitochondrial Ca²⁺ 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²⁺ uptake capacities lead to diminished resting state modes of brain functional connectivity.

Lactate preserves neuronal metabolism and function following antecedent recurrent hypoglycemia

Hypoglycemia occurs frequently during intensive insulin therapy in patients with both type 1 and type 2 diabetes and remains the single most important obstacle in achieving tight glycemic control. Using a rodent model of hypoglycemia, we demonstrated that exposure to antecedent recurrent hypoglycemia leads to adaptations of brain metabolism so that modest increments in circulating lactate allow the brain to function normally under acute hypoglycemic conditions. We characterized 3 major factors underlying this effect. First, we measured enhanced transport of lactate both into as well as out of the brain that resulted in only a small increase of its contribution to total brain oxidative capacity, suggesting that it was not the major fuel. Second, we observed a doubling of the glucose contribution to brain metabolism under hypoglycemic conditions that restored metabolic activity to levels otherwise only observed at euglycemia. Third, we determined that elevated lactate is critical for maintaining glucose metabolism under hypoglycemia, which preserves neuronal function. These unexpected findings suggest that while lactate uptake was enhanced, it is insufficient to support metabolism as an alternate substrate to replace glucose. Lactate is, however, able to modulate metabolic and neuronal activity, serving as a "metabolic regulator" instead.

Role of ongoing, intrinsic activity of neuronal populations for quantitative neuroimaging of functional magnetic resonance imaging-based networks

A primary objective in neuroscience is to determine how neuronal populations process information within networks. In humans and animal models, functional magnetic resonance imaging (fMRI) is gaining increasing popularity for network mapping. Although neuroimaging with fMRI-conducted with or without tasks-is actively discovering new brain networks, current fMRI data analysis schemes disregard the importance of the total neuronal activity in a region. In task fMRI experiments, the baseline is differenced away to disclose areas of small evoked changes in the blood oxygenation level-dependent (BOLD) signal. In resting-state fMRI experiments, the spotlight is on regions revealed by correlations of tiny fluctuations in the baseline (or spontaneous) BOLD signal. Interpretation of fMRI-based networks is obscured further, because the BOLD signal indirectly reflects neuronal activity, and difference/correlation maps are thresholded. Since the small changes of BOLD signal typically observed in cognitive fMRI experiments represent a minimal fraction of the total energy/activity in a given area, the relevance of fMRI-based networks is uncertain, because the majority of neuronal energy/activity is ignored. Thus, another alternative for quantitative neuroimaging of fMRI-based networks is a perspective in which the activity of a neuronal population is accounted for by the demanded oxidative energy (CMR(O2)). In this article, we argue that network mapping can be improved by including neuronal energy/activity of both the information about baseline and small differences/fluctuations of BOLD signal. Thus, total energy/activity information can be obtained through use of calibrated fMRI to quantify differences of ΔCMR(O2) and through resting-state positron emission tomography/magnetic resonance spectroscopy measurements for average CMR(O2).

Pitfalls in Fractal Time Series Analysis: fMRI BOLD as an Exemplary Case

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.

Analysis of time and space invariance of BOLD responses in the rat visual system

Neuroimaging studies of functional magnetic resonance imaging (fMRI) and electrophysiology provide the linkage between neural activity and the blood oxygenation level-dependent (BOLD) response. Here, BOLD responses to light flashes were imaged at 11.7T and compared with neural recordings from superior colliculus (SC) and primary visual cortex (V1) in rat brain--regions with different basal blood flow and energy demand. Our goal was to assess neurovascular coupling in V1 and SC as reflected by temporal/spatial variances of impulse response functions (IRFs) and assess, if any, implications for general linear modeling (GLM) of BOLD responses. Light flashes induced high magnitude neural/BOLD responses reproducibly from both regions. However, neural/BOLD responses from SC and V1 were markedly different. SC signals followed the boxcar shape of the stimulation paradigm at all flash rates, whereas V1 signals were characterized by onset/offset transients that exhibited different flash rate dependencies. We find that IRF(SC) is generally time-invariant across wider flash rate range compared with IRF(V1), whereas IRF(SC) and IRF(V1) are both space invariant. These results illustrate the importance of measured neural signals for interpretation of fMRI by showing that GLM of BOLD responses may lead to misinterpretation of neural activity in some cases.

Fractal analysis of spontaneous fluctuations of the BOLD signal in rat brain

Analysis of task-evoked fMRI data ignores low frequency fluctuations (LFF) of the resting-state the BOLD signal, yet LFF of the spontaneous BOLD signal is crucial for analysis of resting-state connectivity maps. We characterized the LFF of resting-state BOLD signal at 11.7T in α-chloralose and domitor anesthetized rat brain and modeled the spontaneous signal as a scale-free (i.e., fractal) distribution of amplitude power (|A|²) across a frequency range (f) compatible with an |A(f)|² ∝ 1/f(β) model where β is the scaling exponent (or spectral index). We compared β values from somatosensory forelimb area (S1FL), cingulate cortex (CG), and caudate putamen (CPu). With α-chloralose, S1FL and CG β values dropped from ~0.7 at in vivo to ~0.1 at post mortem (p<0.0002), whereas CPu β values dropped from ~0.3 at in vivo to ~0.1 at post mortem (p<0.002). With domitor, cortical (S1FL, CG) β values were slightly higher than with α-chloralose, while subcortical (CPu) β values were similar with α-chloralose. Although cortical and subcortical β values with both anesthetics were significantly different in vivo (p<0.002), at post mortem β values in these regions were not significantly different and approached zero (i.e., range of -0.1 to 0.2). Since a water phantom devoid of susceptibility gradients had a β value of zero (i.e., random), we conclude that deoxyhemoglobin present in voxels post-sacrifice still impacts tissue water diffusion. These results suggest that in the anesthetized rat brain the LFF of BOLD signal at 11.7T follow a general 1/f(β) model of fractality where β is a variable responding to physiology. We describe typical experimental pitfalls which may elude detection of fractality in the resting-state BOLD signal.

Neurovascular and Neurometabolic Couplings in Dynamic Calibrated fMRI: Transient Oxidative Neuroenergetics for Block-Design and Event-Related Paradigms

Functional magnetic resonance imaging (fMRI) with blood-oxygenation level dependent (BOLD) contrast is an important tool for mapping brain activity. Interest in quantitative fMRI has renewed awareness in importance of oxidative neuroenergetics, as reflected by cerebral metabolic rate of oxygen consumption(CMRO2), for supporting brain function. Relationships between BOLD signal and the underlying neurophysiological parameters have been elucidated to allow determination of dynamic changes inCMRO2 by "calibrated fMRI," which require multi-modal measurements of BOLD signal along with cerebral blood flow (CBF) and volume (CBV). But how doCMRO2 changes, steady-state or transient, derived from calibrated fMRI compare with neural activity recordings of local field potential (LFP) and/or multi-unit activity (MUA)? Here we discuss recent findings primarily from animal studies which allow high magnetic fields studies for superior BOLD sensitivity as well as multi-modal CBV and CBF measurements in conjunction with LFP and MUA recordings from activated sites. A key observation is that while relationships between neural activity and sensory stimulus features range from linear to non-linear, associations between hyperemic components (BOLD, CBF, CBV) and neural activity (LFP, MUA) are almost always linear. More importantly, the results demonstrate good agreement between the changes inCMRO2 and independent measures of LFP or MUA. The tight neurovascular and neurometabolic couplings, observed from steady-state conditions to events separated by <200 ms, suggest rapid oxygen equilibration between blood and tissue pools and thus calibrated fMRI at high magnetic fields can provide high spatiotemporal mapping ofCMRO2 changes.

Cerebral oxygen demand for short-lived and steady-state events

Because of the importance of oxidative energetics for cerebral function, extraction of oxygen consumption (CMR(O2)) from blood oxygenation level-dependent (BOLD) signal using multi-modal measurements of blood flow (CBF) and volume (CBV) has become an accepted functional magnetic resonance imaging (fMRI) technique. This approach, termed calibrated fMRI, is based on a biophysical model which describes tissue oxygen extraction at steady-state. A problem encountered for calculating dynamic CMR(O2) relates to concerns whether the conventional BOLD model can be applied transiently. In particular, it is unclear whether calculation of CMR(O2) differs between short and long stimuli. Linearity was experimentally demonstrated between BOLD-related components and neural activity, thereby making it possible to use calibrated fMRI in a dynamic manner. We used multi-modal fMRI and electrophysiology, in alpha-chloralose anesthetized rats during forepaw stimulation to show that respective transfer functions (of BOLD, CBV, CBF) generated by deconvolution with neural activity are time invariant, for events in the millisecond to minute range. These results allowed extraction of a significant component of the BOLD signal that can be ascribed to CMR(O2) transients. We discuss the importance of minimizing residual signal, represented by the difference between modeled and raw signals, in convolution analysis of multi-modal signals.

Tactile and non-tactile sensory paradigms for fMRI and neurophysiologic studies in rodents

Functional magnetic resonance imaging (fMRI) has become a popular functional imaging tool for human studies. Future diagnostic use of fMRI depends, however, on a suitable neurophysiologic interpretation of the blood oxygenation level dependent (BOLD) signal change. This particular goal is best achieved in animal models primarily due to the invasive nature of other methods used and/or pharmacological agents applied to probe different nuances of neuronal (and glial) activity coupled to the BOLD signal change. In the last decade, we have directed our efforts towards the development of stimulation protocols for a variety of modalities in rodents with fMRI. Cortical perception of the natural world relies on the formation of multi-dimensional representation of stimuli impinging on the different sensory systems, leading to the hypothesis that a sensory stimulus may have very different neurophysiologic outcome(s) when paired with a near simultaneous event in another modality. Before approaching this level of complexity, reliable measures must be obtained of the relatively small changes in the BOLD signal and other neurophysiologic markers (electrical activity, blood flow) induced by different peripheral stimuli. Here we describe different tactile (i.e., forepaw, whisker) and non-tactile (i.e., olfactory, visual) sensory paradigms applied to the anesthetized rat. The main focus is on development and validation of methods for reproducible stimulation of each sensory modality applied independently or in conjunction with one another, both inside and outside the magnet. We discuss similarities and/or differences across the sensory systems as well as advantages they may have for studying essential neuroscientific questions. We envisage that the different sensory paradigms described here may be applied directly to studies of multi-sensory interactions in anesthetized rats, en route to a rudimentary understanding of the awake functioning brain where various sensory cues presumably interrelate.

Frequency-dependent tactile responses in rat brain measured by functional MRI

We measured frequency-dependent functional MRI (fMRI) activations (at 11.7 T) in the somatosensory cortex with whisker and forepaw stimuli in the same alpha-chloralose anesthetized rats. Whisker and forepaw stimuli were attained by computer-controlled pulses of air puffs and electrical currents, respectively. Air puffs deflected (+/-2 mm) the chosen whisker(s) in the right snout in the rostral to caudal direction, and electrical currents (2 mA amplitude, 0.3 ms duration) stimulated the left forepaw with subcutaneous copper electrodes placed between the second and fourth digits. In the same subject, unimodal stimulation of whisker and forepaw gave rise to significant blood oxygen level-dependent (BOLD) signal increases in corresponding contralateral somatosensory areas of whisker barrel field (S1BF) and forelimb (S1FL), respectively, with no significant spatial overlap between these regions. The BOLD responses in S1(BF) and S1(FL) regions were found to be differentially variable with frequency of each stimulus type. In the S1BF, a linear increase in the BOLD response was observed with whisker stimulation frequency of up to approximately 12 Hz, beyond which the response seemed to saturate (and/or slightly attenuate) up to the maximum frequency studied (i.e. 30 Hz). In the S1FL, the magnitude of the BOLD response was largest at forepaw stimulation frequency between 1.5 and 3 Hz, beyond which the response diminished with little or no activity at frequencies higher than 20 Hz. The volume of tissue activated by each stimulus type followed a similar pattern to that of the stimulation frequency dependence. These results of bimodal whisker and forepaw stimuli in the same subject may provide a framework to study interactions of different tactile modules, with both fMRI and neurophysiology (i.e. inside and outside the magnet).

Negative BOLD with large increases in neuronal activity

Blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) is widely used in neuroscience to study brain activity. However, BOLD fMRI does not measure neuronal activity directly but depends on cerebral blood flow (CBF), cerebral blood volume (CBV), and cerebral metabolic rate of oxygen (CMRO(2)) consumption. Using fMRI, CBV, CBF, neuronal recordings, and CMRO(2) modeling, we investigated how the signals are related during seizures in rats. We found that increases in hemodynamic, neuronal, and metabolic activity were associated with positive BOLD signals in the cortex, but with negative BOLD signals in hippocampus. Our data show that negative BOLD signals do not necessarily imply decreased neuronal activity or CBF, but can result from increased neuronal activity, depending on the interplay between hemodynamics and metabolism. Caution should be used in interpreting fMRI signals because the relationship between neuronal activity and BOLD signals may depend on brain region and state and can be different during normal and pathological conditions.

NMDA and non-NMDA receptors stimulation causes differential oxidative stress in rat cortical slices

Glutamate receptor activated neuronal cell death is attributed to a massive influx of Ca(2+) and subsequent formation of reactive oxygen species (ROS) but the relative contribution of NMDA and non-NMDA sub-types of glutamate receptors in excitotoxicity is not known. In the present study, we have examined the role of NMDA and non-NMDA receptors in glutamate-induced neuronal injury in cortical slices from young (20+/-2 day) and adult (80+/-5 day) rats. Treatment of slices with glutamate receptor agonists NMDA, AMPA and KA elicited the formation of reactive oxygen species (ROS) and neuronal cell death. In young slices, NMDA receptor stimulation caused a higher ROS formation and neurotoxicity, but KA was more effective in producing ROS and cell death in adult slices. AMPA exhibited an intermediate effect on ROS formation and toxicity in both the age groups. A significant protection in glutamate mediated ROS formation and neurotoxicity was observed in presence of NMDA or/and non-NMDA receptors antagonists APV and NBQX, respectively. This further confirms the involvement of both NMDA and non-NMDA receptors in glutamate mediated neurotoxicity. In adult slices, we did not find positive correlation between ligand induced neurotoxicity and mitochondrial depolarization. Though, NMDA and KA stimulation produced differential effect on ROS formation and neurotoxicity in young and adult slices, the mitochondrial depolarization was higher and comparable on NMDA stimulation in both the age groups as compared to KA, suggesting that the mitochondrial depolarization may not be a good indicator for neurotoxicity. Our results demonstrate that both NMDA and non-NMDA sub-types of glutamate receptors are involved in glutamate mediated neurotoxicity but their relative contribution is highly dependent on the age of the animal.

Xanthine oxidase, nitric oxide synthase and phospholipase A2 produce reactive oxygen species via mitochondria

The formation of reactive oxygen species (ROS) has been suggested to be associated with excitotoxicity but the involvement of cytoplasmic enzymes in ROS formation is not clearly known. In the present study, we examined the role of xanthine oxidase (XO), nitric oxide synthase (NOS) and phospholipase A(2) (PLA(2)) in glutamate-induced oxidative stress in rat cortical slices. Glutamate-induced ROS formation and mitochondrial depolarization were measured in rat cortical slices in presence of allopurinol, L-NAME and 4-bromophenacylbromide, the specific inhibitors of XO, NOS and PLA(2), respectively. Upon stimulation of slices with glutamate, a significant increase in ROS formation and mitochondrial depolarization was observed. However, pretreatment of slices with allopurinol, L-NAME and 4-bromophenacylbromide inhibited the glutamate-induced ROS formation and mitochondrial depolarization. The glutamate-induced ROS formation was dependent on the concentration of these inhibitors and also on the duration of the treatment. Allopurinol was found to be less effective as compared to L-NAME and 4-bromophenacylbromide. The combined treatment of slices with these enzyme inhibitors showed further inhibition in ROS formation and mitochondrial depolarization. The inhibition in ROS formation as well as mitochondrial depolarization by allopurinol, L-NAME and 4-bromophenacylbromide clearly suggests that the activation of XO, NOS and PLA(2) by calcium during glutamate receptor stimulation may release some chemicals which depolarize mitochondria resulting in ROS formation.

Glutamate-induced differential mitochondrial response in young and adult rats

Excitatory amino acid glutamate is involved in neurotransmission in the nervous system but it becomes a potent neurotoxin under variety of conditions. However, the molecular mechanism of excitotoxicity is not known completely. We have studied the influence of glutamate on intracellular calcium and mitochondrial functions in cortical slices from young and adult rats. The slices from both the age groups exhibited comparable intracellular calcium changes upon glutamate stimulation. Glutamate treatment caused a decrease in adenosine 5'-diphosphate/adenosine 5'-triphosphate (ADP/ATP) and an increase in nicotinamide adenine dinucleotide/nicotinamide adenine dinucleotide reduced form (NAD/NADH) ratio in both the age groups but the magnitude and the nature of temporal change was different. Glutamate-induced decrease in ATP/ADP and increase in NAD/NADH ratio was significantly higher in slices from the adult as compared to the young rats. The slices from young rats elicited slightly higher mitochondrial depolarization than adult rats. However, the formation of reactive oxygen species (ROS) and lactate dehydrogenase (LDH) release were significantly higher in adult rats as compared to young rats. Glutamate-induced mitochondrial depolarization, ROS formation and LDH release were highly dependent on the presence of Ca(2+) in the extracellular medium. The treatment of slices with mitochondrial inhibitors rotenone and oligomycin inhibited ROS formation and LDH release substantially. Our results suggest that the glutamate-induced increase in intracellular calcium is not the only factor responsible for neuronal cell death but the mitochondrial functions could be crucial in excitotoxicity.