Abstract
Malformations of cortical development (MCD) is associated with a wide range of developmental delay and drug resistant epilepsy in children. By using resting-state functional magnetic resonance imaging (RS-fMRI) and event-related spectral perturbation (ERSP) of cortical electroencephalography (EEG) data, we tried to investigate the neural changes of spatiotemporal functional connectivity (FC) and fast oscillation (FO) dynamics in a rat model of methylazoxymethanol (MAM)-induced MCD. A total of 28 infant rats with prenatal exposure to MAM and those of age matched 28 controls with prenatal saline exposure were used. RS-fMRI were acquired at postnatal day 15 (P15) and 29 (P29), and correlation coefficient analysis of eleven region of interests (ROI) was done to find the differences of functional networks between four groups. Two hour-cortical EEGs were also recorded at P15 and P29 and the ERSP of gamma (30–80 Hz) and ripples (80–200 Hz) were analyzed. The rats with MCD showed significantly delayed development of superior colliculus-brainstem network compared to control rats at P15. In contrast to marked maturation of default mode network (DMN) in controls from P15 to P29, there was no clear development in MCD rats. The MCD rats showed significantly higher cortical gamma and ripples-ERSP at P15 and lower cortical ripples-ERSP at P29 than those of control rats. This study demonstrated delayed development of FC and altered cortical FO dynamics in rats with malformed brain. The results should be further investigated in terms of the epileptogenesis and cognitive dysfunction in patients with MCD.
Introduction
Malformations of cortical development (MCD) is a group of disorders that has disruption at any step of human brain development including cell proliferation, neuronal migration, or post-migrational cortical organization and connectivity (; ). The clinical feature and course of MCD is widely varied from normal to severe cognitive impairment, and epilepsy is one of the most common presenting symptoms. Moreover, near 90% of patients with epilepsy caused by MCD has been reported to be drug-resistant ().
Although there has rapid evolution of molecular biology, genetics, and neuroimaging techniques in clinics, studies with animal models are still needed for the investigation of MCD pathogenesis due to limitation of clinical studies involving patients or human tissues. There are several animal models of MCD () and one such model is the prenatal methylazoxymethanol (MAM)-induced MCD model. Prenatal exposure of MAM disrupts cell migration resulting the malformations mostly in the hippocampal CA1 and CA2 region, to a lesser extent in the striatum, thalamus, hypothalamus, and cerebral cortex, which are similar to MCDs in human patients (; ; ). Spontaneous epileptic seizures have not been reported so far but several studies indicated a lower threshold for seizures or epileptiform activities (; ; ). MAM-treated rats also had behavioral alteration even at infancy which is consistent with cognitive impairment in patients with MCD (; ).
Recently, resting-state functional magnetic resonance imaging (RS-fMRI) is widely used for mapping large-scale brain networks, the intrinsic functional connectivity (FC) of brain without external stimuli in humans and animal models (; ). Although RS-fMRI on small animals is still scarce, they have the potential to identify the pathologic brain circuits responsible for a neurological disorder. We hypothesized that the structural brain malformation caused the disrupted FC which is associated with behavioral alteration and susceptibility to seizures in MCD rats, and tried to find the abnormal FC of MCD rats using RS-fMRI.
Unfortunately, RS-fMRI is indirect and a low time resolution method to investigate local neuronal activities associated with malformed cortex. As electroencephalography (EEG) is another measure of complex structural-functional dynamics which reflects a direct neuronal activity in high temporal resolution, the interest of combining these two modalities, RS-fMRI and EEG, has grown markedly in last two decades (; ). Classically, cortical EEG has been used to find the epileptic networks or seizure onset zone in human and animal models. High frequency oscillations (HFOs) including gamma oscillation has been suggested as a marker for physiologic phenomenon of cognition () or pathologic seizure onset zone recently (). To quantify the changes of these HFO, several measurement indices are used and event-related spectral perturbation (ERSP) was one of them (; ). Previously, our group also showed the elevated fast oscillation-event-related spectral perturbations (FO-ERSPs) of the pathologic brain lesions in patients with hypsarrhythmia ().
To understand the mechanism of epileptogenicity and cognitive impairment of the brain with MCD, prenatally MAM-treated rats were used in this study. It is hypothesized that prenatally MAM-exposed MCD rats show disrupted resting state FC and abnormal FO dynamics during their brain development. To test this hypothesis and quantitatively measure the changes of the neural system, the FC maps of RS-fMRI data and the cortical FO dynamics were compared between rats with MCD and controls at their age of P15 and P29.
Materials and Methods
Animals
All experiments were approved by the by the Institutional Animal Care and Use Committee of the University of Ulsan College of Medicine and conducted in accordance with the Revised Guide for the Care and Use of Laboratory Animals. Pregnant Sprague–Dawley rat dams (n = 10, Orient Bio Inc., Seoul, South Korea) were acquired on gestational day (G) 14 and housed individually in the animal facility during the remainder of their pregnancy under a 12-h light/dark cycle with free access to food and water. They were injected intraperitoneally with either 0.9% physiological saline (control group, n = 5) or two doses of MAM (15 mg/kg/dose, MRIGlobal, Kansas City, MO, United States) in 10 mL/kg saline at 0800 and 1800 (MCD group, n = 5) on G15. Delivery occurred consistently on G22, which was considered postnatal day (P) 0 for the offspring. The rats used in these experiments are listed in Supplementary Table S1.
Magnetic Resonance Imaging (MRI) Acquisition
Each of eight MCD and control rats was used for imaging data at P15 and P29. All MRI studies were conducted using a 7.0 T/160-mm small-animal imaging system (Bruker Pharmascan, Ettlingen, Germany) with 400 mT/m gradient system and a surface coil for reception and 72 mm volume coil for transmission. During scanning, rats were anesthetized with mixture of 1.0% isoflurane and room air delivered by nose cone and their respiratory rate, electrocardiogram, and the rectal temperature were monitored using small-animal physiological monitoring devices. High-resolution anatomical T2 weighted images were acquired with rapid acquisition with relaxation enhancement (RARE) sequence [TR/TE = 4000/33.0 ms, rare factor = 8, slice thickness = 0.8 mm (total 20 slices), matrix size = 256 × 256, and field of view (FOV) = 25 × 25 mm]. Functional MR data were acquired using a single shot gradient echo based echo-planar image (GE-EPI) sequence positioned parallel to the anterior–posterior commissure plane over the entire brain (TR/TE = 1000/16.734 ms, flip angle = 35′, 300 repetitions, FOV = 25 × 25 mm, matrix size = 96 × 96, no inter slice gap with 10 axial slices).
Thickness and Size Measurement of Neocortex, Hippocampus, and Ventricle
To assess the morphological changes of MCD rats at different developmental ages, the neocortical depths, hippocampal, and ventricle size were analyzed with high resolution T2 weighted images using SPM12 and ImageJ1. To measure the same cortical point of each rat, the anatomical images at of all subjects were overlaid and co-registered each other using SPM12. The cortical thickness of motor, somatosensory, and insular cortices in 3 mm posterior from the bregma, the areas of both hippocampi at 5 mm posterior from the bregma, and the lateral ventricles at 4 mm posterior from the bregma were measured (Figure 1A).
FIGURE 1
Resting-State Networks Analysis
Preprocessing
Several steps of preprocessing were performed using MRIcron software2 and AFNI software (Analysis of Functional NeuroImages3). Preprocessing steps included: (i) Brain extraction: Brain extraction was done manually using MRIcron software. (ii) Registration to standard space: The registrations of functional image to anatomical data and anatomical data to a high-resolution rat brain template were carried out using (FMRIB software library) FSL’s flirt (FMRIB’s Linear Image Registration Tool4) using 12° of freedom affine transformation. The rat brain template was selected for a T2 weighted anatomical scan with field of view of 136 mm × 102 mm × 71 mm and spatial resolution of 1.25 mm × 1.25 mm × 1.25 mm mm (Scalable Brain Atlas5). (iii) Data cleanup: Discarding the first three volumes, slice timing correction, detrending, despiking, motion correction, ventricular, and global signal regression. (iv) Band-pass filtering: Functional images were band-pass filtered between 0.01 and 0.1 Hz. (v) Spatial smoothing: We applied the spatial smoothing using Gaussian kernel FWHM of 0.6 mm to identify relatively large-scale networks across the whole brain of a young rat.
Network Analysis
Resting state analysis was performed using spatial ICA (independent component analysis) by the tool MELODIC (Multivariate Exploratory Linear Optimized Decomposition into Independent Components) from FSL. ICA is being successfully applied method in diverse range of neuroscience for automatically separating various “independent” sources (; ). To find the common spatial patterns of each four groups, we used group-ICA (; ) on the entire set of fMRI data concatenated across the subject of each group. The extracting components were set to 40, which seem to be a reasonable number to provide sub-network separation without incurring into mathematical granularity that produces individual structures and components (; ).
To assess FC change of the individual part with structural abnormalities, five neocortical areas [primary/secondary visual cortex/posterior parietal cortices (V1/V2/PPC), somatosensory cortex (SS), motor cortex (Motor), auditory cortex (Au), and anterior/posterior cingulate cortices (aCg/pCg)], five subcortical areas [hippocampus, basal ganglia, thalamus, hypothalamus, and superior colliculus (SC)] and brainstem were designated the region of interests (ROI) and the Pearson correlation coefficient (CC) among the mean signal intensity time courses of total 11 ROI were calculated. We represent the 11 ROI in Supplementary Figure S1. There was a large gap in the brain size and morphology between P15 and P29, we used the different masks for those P15 and P29 rats. A Fisher’s r-to-z transformation was applied to each correlation map to obtain an approximately normal distribution of the FC values. Two CC maps were compared between those with rats with MCD and controls at each timepoint P15 and P29 to show the abnormal FC of MCD rats during early developmental period. To adjust the random effect of clustering from the same mother and repeated measure effect at two timepoints, we performed the linear mixed model analysis.
Cortical Electroencephalography (EEG) Recording and Analysis
For intracranial EEG recording, electrodes were surgically implanted in each of ten MCD and control rats under sedation with ketamine/xylazine (50/7 mg/kg in 10 mL/kg saline IP) at P13 and P27, respectively. Two cortical electrodes were implanted over both somatosensory cortices in each of five MCD and control rats. They are connected to a multiple socket and secured to the skull with dental acrylic. At P15 and P29, 2-h non-sedated EEGs of cortical electrodes were recorded with simultaneous and synchronized video using the Twin EEG system (Grass Technologies Corp.). The sampling rate was 400 Hz with 0.1 Hz high-pass filter. The 80 s of each artifact-free, resting state EEG were collected and reassessed. The extracted EEGs were extended temporally with 1, 30, and 80 Hz high pass filters to detect visible gamma (30–80 Hz) and ripples (80–200 Hz) with sensitivities of 50 and 25 μV/mm, respectively. For quantitative estimation, ERSP of gamma and ripples from each epoch was calculated using the EEGLAB toolbox of MATLAB 2017b. ERSP (log) was used to show dynamic brain changes, with the zero point in each epoch set as the baseline. In each 1-s epoch, ERSPs were analyzed with fast Fourier transform and Hanning window tapering. The ERSP formula for averaged estimates across data trials (n trials) is defined below:
where Fk(f,t) is the spectral estimate of trial k at frequency f and time t.
The cortical gamma and ripples-ERSP extracted from 80 epoch in each rat were compared between those with MCD rats and controls at different ages, P15 and P29 using linear mixed model analysis with adjusting random effects of clustering from same mother and repetitive measured values (80 epoch) from each rat.
Results
Comparisons of Cortical Thickness, Hippocampal, and Ventricle Sizes
Using in vivo MRI, we could demonstrate the widespread cortical thinning of MCD rats of this study in motor, somatosensory and insular cortices when compared to control rats (Figure 1, unpaired t-test, P15-motor: t = 7.001, df = 12.507, p = 0.001; P15-somatosensory: t = 6.032, df = 10.761, p = 0.001; P15-insular: t = 3.348, df = 10.864, p = 0.007; P29-motor: t = 4.703, df = 11.703, p = 0.001; P29-somatosensory: t = 7.605, df = 10.748, p = 0.001; P29-insular: t = 8.576, df = 12.916, p = 0.001) as previously shown (). The two-dimensional measurements of bilateral hippocampi and lateral ventricles also revealed significant reduction in those of MCD rats at P15 and P29 (Unpaired t-test; P15-hippo: t = 6.632, df = 11.854, p = 0.001; P15-ventricle: t = 5.939, df = 8.115, p = 0.001; P29-hippo: t = 11.187, df = 13.060, p = 0.001; P29-ventricle: t = 19.702, df = 9.175, p = 0.001).
Resting-State Networks Identification
Using group-ICA, four functional networks including default mode, motor-somatosensory, basal ganglia-hypothalamus-hippocampus, and brainstem network were found in each group (Supplementary Figure S2). To quantitatively measure the group difference of FC between the cortical and subcortical areas, the correlation coefficient (CC) among 11 ROIs were analyzed and compared to each other. At P15, the whole brain FC was lower in both groups (Figures 2A,B) and the intra-cortical, intra-subcortical positive FC and neocortical-subcortical negative connectivity became stronger in both groups at P29 (Figures 2C,D). The MCD rats showed significantly reduced SC-brainstem connectivity compared to control rats at P15 and this difference disappeared at P29 (Figures 2E,F; P15: t = 4.45, p = 0.047; P29: t = −0.03, p = 0.976). When we compare the FC of P15 and P29 in each group, control rats showed a significant increase in aCg/pCg-Au FC (Figure 2G, t = −4.97, p = 0.038) and V1/V2/PPC-Au FC (Figure 2G, t = −6.13, p = 0.026) and MCD rats showed a significant increase in SC-brainstem FC (Figure 2H, t = −5.4, p = 0.033) from P15 to P29. To show the effect magnitude of cortical malformation and development, the CC of these significantly changed networks were shown in Figure 3.
FIGURE 2
FIGURE 3
Cortical EEG Analysis
With 1 Hz high pass filters, visual analysis of cortical EEG showed relatively increase of fast activities in MCD rats especially at P15. The examples of EEGs on visual analysis are shown in Figure 4.
FIGURE 4
To quantify the time-related shift of FO, the mean cortical gamma and ripples-ERSP of rats with MCD or controls were compared at P15 and P29. The MCD rats at P15 showed significantly higher gamma- and ripple-ERSPs than controls (Figures 4E,F, cortical gamma: 4.41 ± 0.40 vs. 4.36 ± 0.32, p = 0.021, df = 278.906; cortical ripples: 4.43 ± 0.53 vs. 4.38 ± 0.24, p = 0.029, df = 201.113). Whereas cortical ripples-ERSP of MCD rats at P29 was significantly lower than that of control rats (Figure 4F, 4.25 ± 0.61 vs. 4.52 ± 0.39; p < 0.001; df = 119.808).
Discussion
Malformations of cortical development are commonly associated with refractory epilepsy and cognitive impairment in children, yet the brain network alteration and developmental electrophysiological changes of this structurally abnormal brain are not clearly defined. Using a rat model of MCD, we analyzed the resting-state FC and FO dynamics of MCD brains to identify the developmental changes of functional networks and cortical electrical activities caused by the malformed brain in this study.
Both in humans and animal models, RS-fMRI has been used to study FC during rest (; ; ) and this spontaneous fluctuation of BOLD signal was seen as similar pattern within same species (e.g., rats) or different species such as human and rodents (). In this study, intrinsic connectivity networks (ICN) were successfully identified in all four groups (Supplementary Figure S2) which was consistent with the previously established essential networks in the anesthetized (; ) or awake adult rats (), and 2-week-old infant rats (). Furthermore, this study also demonstrated the serial representative FC pattern at two different developmental ages.
To gain insight for the impaired network properties in diseased brain, it is necessary to examine whether there is significant difference in the degree of variability between the control and disease group (; ). At first, we could confirm the widespread structural disruptions in MCD rats compared to controls in a previous study (Figure 1; ). Using these animals with malformed brain, we analyzed the correlation coefficient among the major structures of brain (11 ROIs) to investigate the individual changes of FC and we found abnormal functional network development in rats with MCDs. During the developmental period from P15 to P29, positive correlation of short-range FC (intra-neocortical and intra-subcortical networks) and negative correlation of long-range FC became stronger (Figures 2A–D) as shown in previous human studies (; ). The rats with MCD showed reduced SC-brainstem connectivity compared to control rats at P15 but achieved catch-up development of SC-brainstem connectivity at P29 (Figures 2E,H, 3A).
Superior colliculus has strong descending projection to brainstem seizure circuitry, including direct projections to nucleus reticularis pontis oralis (). The activation of SC exerted broad-spectrum anticonvulsant actions () and lesions of SC markedly attenuated seizure in genetically epilepsy-prone rats (). Consistent with previous studies which showed the important role of SC-brainstem circuit in generation of the spasms (; ; ; ), this delayed development of SC-brainstem network in MCD rats can be associated with the spasms susceptibility in these animals ().
In addition, the development of default mode network (DMN), aCg/pCg-Au and V1/V2/PPC-Au of MCD rats were observed only in control rats at this specific time periods (Figures 2G, 3B,C). The DMN is known as the basal network since many goal-oriented tasks deactivate this network () and a similar network has also been observed in awake and anesthetized rats (; ; ). Many neurological and psychiatric disorders including schizophrenia (), Alzheimer’s disease (), autism () have been linked to DMN dysregulation. Thus, the developmental pattern of the DMN in control rats of this study can be considered as normal maturation of DMN and no significant changes of DMN-FC of the MCD rats can be interpreted as the delayed integration of DMN component which can be associated with the malformation of the consistent regions and the cognitive deficit and poor behavioral performances in these rats (; ; ).
We also monitored direct electrophysiological signals with cortical electrodes to focus the local neuronal activities associated with this malformed cortex because RS-fMRI has the limitations of poor time-resolution and inability to detect neural electrical activities. Among the electrical activities of brain, the FO (over 30 Hz) are alleged as have a critical role in integration of neural networks during cognitive processes (; ) and these FO including gamma and ripples activities are associated with synchronized activation of reciprocally connected excitatory pyramidal neurons and inhibitory interneurons (; ). Moreover, the abnormal regulation of high frequency band activities were found in human patients with epileptogenic lesions (; ), schizophrenia (; ), or autism (), which suggests the FO as a biomarker of neuropsychiatric disease ().
We hypothesized that the abnormal neural network of MCD rats will lead to disrupted FO dynamics and the rats with MCD showed significant elevated gamma and ripples-ERSP at P15 and decreased ripple-ERSP at P29. The ERSP has a strong ability to detect the time-related shift of the specific band frequency (; ; ) and these elevated FO-ERSP of infant rats with MCD shows impaired FO dynamics in these rats. In the brains of normally developed infant rats, the density of glutamatergic synapses is too low to make physiological HFOs and epileptic HFOs can be triggered by an additional drive provided by excitatory GABA (). Thus, the increased FO-ERSP of rat with MCD at P15 in this study reflect the dysregulation of FO during their early development, which can be driven by abnormal excitatory GABA mediated depolarization and can be associated with the cognitive impairment and seizure susceptibility of those rats (; ).
On the other hands, physiological HFOs in rats appear only after the developmental GABA switch from excitation to inhibition and the formations of sufficient glutamatergic synapses after P15 (). Previous immunohistochemical studies with rats with exposure to MAM in utero (; ; ) showed impaired neuronal network formations in these animals and the decreased ripple ERSP of rat with MCD at P29 may be the results.
Unfortunately, the direct correlation between disrupted large-scale FC and dynamic changes of local neuronal circuit could not be assessed due to the lack of the simultaneous recording of RS-fMRI and cortical EEG. However, this is the first study which identifies the change of FC and FO dynamics in young MCD rats at two different developmental periods. The rats with MCD showed delayed functional networks development and impaired cortical oscillatory dynamics during their infancy and childhood. These developmental in vivo imaging and EEG changes can be suggested as biomarkers for malformed brain. Further, these features can be used to get insights for the pathophysiology of cognitive dysfunction and epilepsy associated with MCD.
Statements
Data availability statement
The datasets generated for this study are available on request to the corresponding authors.
Ethics statement
All experiments were approved by the Institutional Animal Care and Use Committee of the University of Ulsan College of Medicine and conducted in accordance with the Revised Guide for the Care and Use of Laboratory Animals.
Author contributions
M-JK, M-SY, W-HS, and T-SK conceived and designed the research and interpreted results of the experiments. M-JK, ML, E-JK, and YJ performed the experiments. M-JK and YJ analyzed the data. M-JK and M-SY prepared the figures and drafted the manuscript. M-SY, W-HS, and T-SK edited and revised the manuscript. All authors approved the final version of the manuscript.
Funding
This work was supported by the Basic Science Research Program through the National Research Foundation of Korea funded by the Ministry of Education (NRF-2018R1D1A1A02085427) and a grant (2020IL0008) from Asan Institute for Life Sciences, Seoul, South Korea.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2020.00711/full#supplementary-material
References
1
AhnaouA.HuysmansH.Van De CasteeleT.DrinkenburgW. (2017). Cortical high gamma network oscillations and connectivity: a translational index for antipsychotics to normalize aberrant neurophysiological activity.Transl. Psychiatry7:1285.
2
BajicD.CraigM. M.BorsookD.BecerraL. (2016). Probing intrinsic resting-state networks in the infant rat brain.Front. Behav. Neurosci.10:192. 10.3389/fnbeh.2016.00192
3
BarabanS. C.SchwartzkroinP. A. (1995). Electrophysiology of CA1 pyramidal neurons in an animal model of neuronal migration disorders: prenatal methylazoxymethanol treatment.Epilepsy Res.22145–156. 10.1016/0920-1211(95)00045-3
4
BaraditsM.KakusziB.BalintS.FullajtarM.ModL.BitterI.et al (2018). Alterations in resting-state gamma activity in patients with schizophrenia: a high-density EEG study.Eur. Arch. Psychiatry Clin. Neurosci.269429–437. 10.1007/s00406-018-0889-z
5
BarkovichA. J.GuerriniR.KuznieckyR. I.JacksonG. D.DobynsW. B. (2012). A developmental and genetic classification for malformations of cortical development: update 2012.Brain1351348–1369. 10.1093/brain/aws019
6
BecerraL.PendseG.ChangP. C.BishopJ.BorsookD. (2011). Robust reproducible resting state networks in the awake rodent brain.PLoS One6:e25701. 10.1371/journal.pone.0025701
7
BeckmannC. F.SmithS. M. (2004). Probabilistic independent component analysis for functional magnetic resonance imaging.IEEE Trans. Med. Imaging23137–152. 10.1109/tmi.2003.822821
8
CalhounV. D.AdaliT.PearlsonG. D.PekarJ. J. (2001). A method for making group inferences from functional MRI data using independent component analysis.Hum. Brain Mapp.14140–151. 10.1002/hbm.1048
9
ChuganiH. T.ShewmonD. A.SankarR.ChenB. C.PhelpsM. E. (1992). Infantile spasms: II. Lenticular nuclei and brain stem activation on positron emission tomography.Ann. Neurol.31212–219. 10.1002/ana.410310212
10
ColciaghiF.FinardiA.FrascaA.BalossoS.NobiliP.CarrieroG.et al (2011). Status epilepticus-induced pathologic plasticity in a rat model of focal cortical dysplasia.Brain1342828–2843. 10.1093/brain/awr045
11
DelormeA.MakeigS. (2004). EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis.J. Neurosci. Methods1349–21. 10.1016/j.jneumeth.2003.10.009
12
DennisE. L.ThompsonP. M. (2013). Mapping connectivity in the developing brain.Int. J. Dev. Neurosci.31525–542.
13
FairD. A.DosenbachN. U.ChurchJ. A.CohenA. L.BrahmbhattS.MiezinF. M.et al (2007). Development of distinct control networks through segregation and integration.Proc. Natl. Acad. Sci. U.S.A.10413507–13512. 10.1073/pnas.0705843104
14
FoxM. D.SnyderA. Z.VincentJ. L.CorbettaM.Van EssenD. C.RaichleM. E. (2005). The human brain is intrinsically organized into dynamic, anticorrelated functional networks.Proc. Natl. Acad. Sci. U.S.A.1029673–9678. 10.1073/pnas.0504136102
15
Gonzalez-BurgosG.ChoR. Y.LewisD. A. (2015). Alterations in cortical network oscillations and parvalbumin neurons in schizophrenia.Biol. Psychiatry771031–1040. 10.1016/j.biopsych.2015.03.010
16
GorgesM.RoselliF.MullerH. P.LudolphA. C.RascheV.KassubekJ. (2017). Functional connectivity mapping in the animal model: principles and applications of resting-state fMRI.Front. Neurol.8:200. 10.3389/fneur.2017.00200
17
GreiciusM. D.SrivastavaG.ReissA. L.MenonV. (2004). Default-mode network activity distinguishes Alzheimer’s disease from healthy aging: evidence from functional MRI.Proc. Natl. Acad. Sci. U.S.A.1014637–4642. 10.1073/pnas.0308627101
18
GuerriniR.DobynsW. B. (2014). Malformations of cortical development: clinical features and genetic causes.Lancet Neurol.13710–726. 10.1016/s1474-4422(14)70040-7
19
HernanA. E.MahoneyJ. M.CurryW.RichardG.LucasM. M.MasseyA.et al (2018). Environmental enrichment normalizes hippocampal timing coding in a malformed hippocampus.PLoS One13:e0191488. 10.1371/journal.pone.0191488
20
HillaryF. G.RomanC. A.VenkatesanU.RajtmajerS. M.BajoR.CastellanosN. D. (2015). Hyperconnectivity is a fundamental response to neurological disruption.Neuropsychology2959–75. 10.1037/neu0000110
21
HutchisonR. M.MirsattariS. M.JonesC. K.GatiJ. S.LeungL. S. (2010). Functional networks in the anesthetized rat brain revealed by independent component analysis of resting-state FMRI.J. Neurophysiol.1033398–3406. 10.1152/jn.00141.2010
22
HutchisonR. M.WomelsdorfT.AllenE. A.BandettiniP. A.CalhounV. D.CorbettaM.et al (2013). Dynamic functional connectivity: promise, issues, and interpretations.Neuroimage80360–378. 10.1016/j.neuroimage.2013.05.079
23
JamesC. J.HesseC. W. (2005). Independent component analysis for biomedical signals.Physiol. Meas.26R15–R39.
24
JiruskaP.Alvarado-RojasC.SchevonC. A.StabaR.StaceyW.WendlingF.et al (2017). Update on the mechanisms and roles of high-frequency oscillations in seizures and epileptic disorders.Epilepsia581330–1339. 10.1111/epi.13830
25
JobstB. C.EngelJ.Jr. (2015). Is it time to replace epileptic spikes with fast ripples?Neurology85114–115. 10.1212/wnl.0000000000001750
26
KagaK.MarshR. R.FukuyamaY. (1982). Auditory brain stem responses in infantile spasms.Int. J. Pediatr. Otorhinolaryngol.457–67. 10.1016/0165-5876(82)90078-7
27
KennedyD. P.RedcayE.CourchesneE. (2006). Failing to deactivate: resting functional abnormalities in autism.Proc. Natl. Acad. Sci. U.S.A.1038275–8280. 10.1073/pnas.0600674103
28
KimE. H.YumM. S.LeeM.KimE. J.ShimW. H.KoT. S. (2017). A new rat model of epileptic spasms based on methylazoxymethanol-induced malformations of cortical development.Front. Neurol.8:271. 10.3389/fneur.2017.00271
29
KimM. J.YumM. S.YehH. R.KoT. S. (2018). Fast oscillation dynamics during hypsarrhythmia as a localization biomarker.J. Neurophysiol.119679–687. 10.1152/jn.00662.2017
30
LadoF. A.MosheS. L. (2002). Role of subcortical structures in the pathogenesis of infantile spasms: what are possible subcortical mediators?Int. Rev. Neurobiol.49115–140. 10.1016/s0074-7742(02)49010-1
31
Le Van QuyenM.KhalilovI.Ben-AriY. (2006). The dark side of high-frequency oscillations in the developing brain.Trends Neurosci.29419–427. 10.1016/j.tins.2006.06.001
32
LewisD. A.CurleyA. A.GlausierJ. R.VolkD. W. (2012). Cortical parvalbumin interneurons and cognitive dysfunction in schizophrenia.Trends Neurosci.3557–67. 10.1016/j.tins.2011.10.004
33
LiangZ.KingJ.ZhangN. (2011). Uncovering intrinsic connectional architecture of functional networks in awake rat brain.J. Neurosci.313776–3783. 10.1523/jneurosci.4557-10.2011
34
LodgeD. J.BehrensM. M.GraceA. A. (2009). A loss of parvalbumin-containing interneurons is associated with diminished oscillatory activity in an animal model of schizophrenia.J. Neurosci.292344–2354. 10.1523/jneurosci.5419-08.2009
35
LuH.ZouQ.GuH.RaichleM. E.SteinE. A.YangY. (2012). Rat brains also have a default mode network.Proc. Natl. Acad. Sci. U.S.A.1093979–3984. 10.1073/pnas.1200506109
36
LucasM. M.Lenck-SantiniP. P.HolmesG. L.ScottR. C. (2011). Impaired cognition in rats with cortical dysplasia: additional impact of early-life seizures.Brain1341684–1693. 10.1093/brain/awr087
37
LuhmannH. J. (2016). Models of cortical malformation–Chemical and physical.J. Neurosci. Methods26062–72. 10.1016/j.jneumeth.2015.03.034
38
MakeigS.DebenerS.OntonJ.DelormeA. (2004). Mining event-related brain dynamics.Trends Cogn. Sci.8204–210.
39
MandinoF.CerriD. H.GarinC. M.StraathofM.Van TilborgG. A. F.ChakravartyM. M.et al (2019). Animal functional magnetic resonance imaging: trends and path toward standardization.Front. Neuroinform.13:78. 10.3389/fninf.2019.00078
40
MerrillM. A.CloughR. W.JobeP. C.BrowningR. A. (2003). Role of the superior colliculus and the intercollicular nucleus in the brainstem seizure circuitry of the genetically epilepsy-prone rat.Epilepsia44305–314. 10.1046/j.1528-1157.2003.31802.x
41
NowakM.ZichC.StaggC. J. (2018). Motor cortical gamma oscillations: what have we learnt and where are we headed?Curr. Behav. Neurosci. Rep.5136–142. 10.1007/s40473-018-0151-z
42
PellockJ. M.HrachovyR.ShinnarS.BaramT. Z.BettisD.DlugosD. J.et al (2010). Infantile spasms: a U.S. consensus report.Epilepsia512175–2189. 10.1111/j.1528-1167.2010.02657.x
43
RaichleM. E.MacleodA. M.SnyderA. Z.PowersW. J.GusnardD. A.ShulmanG. L. (2001). A default mode of brain function.Proc. Natl. Acad. Sci. U.S.A.98676–682.
44
RaichleM. E.MintunM. A. (2006). Brain work and brain imaging.Annu. Rev. Neurosci.29449–476. 10.1146/annurev.neuro.29.051605.112819
45
RedgraveP.MitchellI. J.DeanP. (1987). Descending projections from the superior colliculus in rat: a study using orthograde transport of wheatgerm-agglutinin conjugated horseradish peroxidase.Exp. Brain Res.68147–167.
46
RojasD. C.WilsonL. B. (2014). gamma-band abnormalities as markers of autism spectrum disorders.Biomark. Med.8353–368. 10.2217/bmm.14.15
47
SchwarzA. J.GassN.SartoriusA.RisterucciC.SpeddingM.SchenkerE.et al (2013). Anti-correlated cortical networks of intrinsic connectivity in the rat brain.Brain Connect.3503–511. 10.1089/brain.2013.0168
48
SierakowiakA.MonnotC.AskiS. N.UppmanM.LiT. Q.DambergP.et al (2015). Default mode network, motor network, dorsal and ventral basal ganglia networks in the rat brain: comparison to human networks using resting state-fMRI.PLoS One10:e0120345. 10.1371/journal.pone.0120345
49
SinghS. C. (1977). Ectopic neurones in the hippocampus of the postnatal rat exposed to methylazoxymethanol during foetal development.Acta Neuropathol.40111–116. 10.1007/bf00688698
50
SmithS. M.BeckmannC. F.AnderssonJ.AuerbachE. J.BijsterboschJ.DouaudG.et al (2013). Resting-state fMRI in the human connectome project.Neuroimage80144–168.
51
SmithaK. A.Akhil RajaK.ArunK. M.RajeshP. G.ThomasB.KapilamoorthyT. R.et al (2017). Resting state fMRI: a review on methods in resting state connectivity analysis and resting state networks.Neuroradiol. J.30305–317. 10.1177/1971400917697342
52
SoperC.WickerE.KulickC. V.N’gouemoP.ForcelliP. A. (2016). Optogenetic activation of superior colliculus neurons suppresses seizures originating in diverse brain networks.Neurobiol. Dis.87102–115. 10.1016/j.nbd.2015.12.012
53
SpatzM.LaqueurG. L. (1968). Transplacental chemical induction of microencephaly in two strains of rats.I Proc. Soc. Exp. Biol. Med.129705–710. 10.3181/00379727-129-33404
54
SumiyoshiA.RieraJ. J.OgawaT.KawashimaR. (2011). A mini-cap for simultaneous EEG and fMRI recording in rodents.Neuroimage541951–1965. 10.1016/j.neuroimage.2010.09.056
55
UpadhyayJ.BakerS. J.ChandranP.MillerL.LeeY.MarekG. J.et al (2011). Default-mode-like network activation in awake rodents.PLoS One6:e27839. 10.1371/journal.pone.0027839
56
WhiteR. S.SiegelS. J. (2016). Cellular and circuit models of increased resting-state network gamma activity in schizophrenia.Neuroscience32166–76. 10.1016/j.neuroscience.2015.11.011
57
Whitfield-GabrieliS.ThermenosH. W.MilanovicS.TsuangM. T.FaraoneS. V.MccarleyR. W.et al (2009). Hyperactivity and hyperconnectivity of the default network in schizophrenia and in first-degree relatives of persons with schizophrenia.Proc. Natl. Acad. Sci. U.S.A.1061279–1284. 10.1073/pnas.0809141106
58
WirsichJ.RidleyB.BessonP.JirsaV.BenarC.RanjevaJ. P.et al (2017). Complementary contributions of concurrent EEG and fMRI connectivity for predicting structural connectivity.Neuroimage161251–260. 10.1016/j.neuroimage.2017.08.055
59
ZijlmansM.JiruskaP.ZelmannR.LeijtenF. S.JefferysJ. G.GotmanJ. (2012). High-frequency oscillations as a new biomarker in epilepsy.Ann. Neurol.71169–178. 10.1002/ana.22548
Summary
Keywords
malformations of cortical development, resting state functional magnetic resonance images, functional connectivity, event-related spectral perturbation, gamma, ripples
Citation
Kim M-J, Yum M-S, Jo Y, Lee M, Kim E-J, Shim W-H and Ko T-S (2020) Delayed Functional Networks Development and Altered Fast Oscillation Dynamics in a Rat Model of Cortical Malformation. Front. Neurosci. 14:711. doi: 10.3389/fnins.2020.00711
Received
16 January 2020
Accepted
12 June 2020
Published
18 August 2020
Volume
14 - 2020
Edited by
Xu Lei, Southwest University, China
Reviewed by
Daniele Corbo, University of Brescia, Italy; Joanes Grandjean, Radboud University Medical Center, Netherlands
Updates
Copyright
© 2020 Kim, Yum, Jo, Lee, Kim, Shim and Ko.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mi-Sun Yum, misun.yum@gmail.comWoo-Hyun Shim, ssimu81@gmail.com
†These authors have contributed equally to this work
This article was submitted to Brain Imaging Methods, a section of the journal Frontiers in Neuroscience
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.