Extrahippocampal Seizure and Memory Circuits Overlap

Abstract Seizures cause retrograde amnesia. We have previously demonstrated that seizures erode recently formed memories through shared ensembles and mechanisms in the CA1 region of the hippocampus. Here, we tested whether seizure circuits overlap spatial memory circuits outside of the CA. Spatial memory is consolidated by the hippocampal-cortical coupling that are connected via multiple pathways. We tested whether a seizure invades structures involved in memory consolidation by using the activity reporter TRAP2 mice. T-maze alternation learning activated neurons in the dentate gyrus (DG), mediodorsal thalamus (MD), retrosplenial cortex (RSC), and medial prefrontal cortex (mPFC). This spatial memory relies on the plasticity of the AMPA receptor GluA1 subunit. GluA1 knock-out (KO)/TRAP2 mice did not learn to alternate, and structures interposed between the hippocampus and the cortex were not active. A seizure prevented the recall of alternation memory and activated memory-labeled structures. There was a widespread overlap between learning-activated ensembles and seizure-activated neurons, which likely contributes to retrograde amnesia.


Introduction
Seizures cause retrograde amnesia, and the mechanisms are poorly understood. Seizures cause transient declarative memory deficits that occur even in patients who do not have underlying hippocampal pathology, in contrast to fixed epilepsy-memory deficits (Hoppe, 2007). Electroconvulsive therapy-induced seizures also cause an amnestic effect on recently acquired event memories (Duncan, 1949;Chorover and Schiller, 1965;Squire et al., 2001). We seek to understand the mechanisms of seizure-induced retrograde amnesia.
The parietal, RSC, anterior cingulate cortex (ACC), and mPFC display ripple oscillations (100-150 Hz) concurrent with hippocampal ripples, which are strengthened during sleep following learning (Khodagholy et al., 2017). Pathologic interictal epileptiform discharges (IEDs) in the hippocampus are correlated with decreased ripple occurrence in the mPFC (Gelinas et al., 2016). On the other hand, hippocampal IEDs evoke mPFC spindle oscillations (9-16 Hz) via induction of cortical DOWN state, which is characterized by d waves associated with neuronal hyperpolarization in the deep cortical layers, possibly leading to consolidation deficits. These studies suggested that epileptiform discharges interrupt physiological mechanisms coupling the hippocampus to the cortex, which may impair cognition in patients with temporal lobe epilepsy.
Mapping seizure and declarative memory-activated neurons is an alternative approach to understanding seizure-induced amnesia. We mapped neurons engaged in storing memory and its retrieval using activity reporter mice (Josselyn and Tonegawa, 2020). Seizure networks can be mapped with a cellular resolution using activity reporter mice as well (Dabrowska et al., 2019;Adotevi et al., 2020;Brodovskaya et al., 2021). Seizures like memory often travel from the hippocampus to the cortex (Adotevi et al., 2020). We tested whether seizures activate the memory circuit to propagate.
We used TRAP2 mice that express a fluorescent protein under the immediate early gene promoter c-Fos (Tasaka et al., 2018) to tag spatial memory-activated neurons beyond the classical CA1 region permanently and in the extrahippocampal circuit. The T-maze spatial learning task is dependent on the plasticity of the GluA1 subunit of AMPA receptors. We used GluA1 knock-out (KO) mice crossed with TRAP2 mice to confirm that labeled activated neurons were specific to spatial memory and not just locomotion or spatial navigation. A seizure shortly after a mouse has learned to alternate for a reward on the T-maze degraded performance the following day (Naik et al., 2021). We tested whether seizure-activated neuronal ensembles overlapped spatial memory-activated neurons by performing TRAP2 and ARC dual-labeling. We propose that seizures cause retrograde amnesia by engaging the circuits of memory consolidation.

Animal model
All studies were performed following the approved protocols of the Animal Care and Use Committee. For tagging activated neuronal ensembles following learning, we used TRAP2 mice generated by crossing mice expressing tdTomato from the Rosa locus [(B6.Cg-Gt (ROSA)26Sor tm9(CAG-tdTomato)Hze /J, ME, #007909; The Jackson Laboratory] to mice expressing Cre-ER under cFos promoter (Fos 2A-iCreER , #030323; Tasaka et al., 2018). AMPA receptor KO mice (GluA1 global KO) were obtained from Seeburg (Max-Plank Institute for Medical Research, Heidelberg, Germany) and have been characterized in detail previously (Reisel et al., 2002;Adotevi et al., 2020). To permanently label the cFos-expressing neurons following learning in GluA1 KO mice, we crossed GluA1 global KO mice with TRAP2 mice to generate GluA1 KO/TRAP2 mice. Mice of both sexes were maintained on 12 h light (6 A.M. to 6 P.M.)/12 h dark (6 P.M. to 6 A.M.) cycle in a homecage (12 Â 6 Â 5 in) with corn cob bedding and had ad libitum access to food and water. For genotyping, KAPA Biosystems kit was used.

Delayed discrete trial rewarded alternation task in Tmaze
We trained TRAP2 mice (six to eight weeks, n = 40) individually housed to alternate on a T-maze for a reward of sweetened condensed milk. We food-deprived mice to no less than 80% of their beginning body weight a week before to motivate them for the reward. After a week of food restriction and familiarization with sweetened milk, we trained mice in a delayed discrete trial rewarded alternation task for 5 d (Naik et al., 2021). Days 1 and 2 served as a training session wherein no delay was introduced between sample and choice runs; however, on day 3 and afterward, a delay of 20 s was introduced between sample and choice trials. Percent correct choices were recorded for each mouse and also sample and choice latencies.
To confirm the behavioral phenotype of hyperactivity in the GluA1 KO mice, we tested the GluA1 KO/TRAP2 and littermate controls (TRAP2) in an Open Field Activity Monitor (Nodulus Instruments). The mice were individually tested for 20 min, and activity maps were used for phenotyping. This behavioral phenotyping was used in addition to the genotypic characterization of the GluA1 KO mice. Data are presented as heatmaps of the activity in the Open Field Activity Monitor.

Labeling training activated neuronal ensembles
To permanently label the activated population of neurons following training and learning in the T-maze task, separate cohorts of TRAP2 mice were injected with fastacting 4-hydroxy tamoxifen (4-OHT) 60 min after the training on days 1-3. We imaged learning-activated neuronal ensembles in CA1, dentate gyrus (DG), and RSC regions across multiple training sessions. Separate cohorts of mice were injected with 4-OHT on days 1-3. To avoid transfer of handling-induced neuronal activation, the mice were returned to home cages only an hour after 4-OHT injection. Mice were transcardially perfused with 4% paraformaldehyde (PFA) in 0.1 M PBS a week after 4-OHT injections and stored overnight in 4% PFA. The brains were cut coronally (40 mm), stained with DAPI or NeuN, and imaged at 10Â/0.45 NA magnification on a Nikon Eclipse Ti2 fully motorized confocal microscope. Imaris version 9.8.0 (Bitplane Scientific) was used to create composite images. Unbiased manual counting of tdT 1ve neurons was done individually by two people blinded to the treatment in every slice. tdT 1ve neurons with well-defined cell bodies and processes in the granule cell layer of DG were counted only. Scrolling over the whole Z-stack was performed to avoid missing any tagged neurons. Fiji ImageJ was used for automated counting of tdT 1ve neurons in RSC, mPFC, MD, and paraventricular thalamus (PVT). Reference atlases from Allen as well as Paxinos were used to define the boundaries of the target brain regions (DG, mPFC, RSC, PVT, and MD) in coronal slices.
TRAPing seizure-activated neuronal populations TRAP2 mice (six to eight weeks old) were injected with PTZ (40 mg/kg, i.p.) and observed for a generalized tonicclonic seizure. Mice that experienced a seizure were injected with 4-OHT 60 min after the seizure. Mice that received saline injections served as a control. All mice injected with 4-OHT were perfused after a week and processed for imaging as mentioned above. Tagged neuronal populations following a seizure were counted in the brain regions like DG, RSC, mPFC, MD, PVT in every slice, and presented as mean 6 SEM.
Dual-labeling of neuronal ensembles activated following rewarded alternation learning and seizure First, using 4-OHT injections on day 2 of the T-maze, the learning-activated neuronal ensembles were permanently tagged with tdTomato. The same mouse received a PTZ injection a week later to cause a seizure. Mice were killed 60 min after the seizure. Arc immunohistochemistry was performed using the antibody (1:2000, Rb polyclonal, 156003, Synaptic Systems). Dual-labeling results were confirmed and validated using an anti-cFos antibody (1:1000, Rb polyclonal, ab190289, Abcam) as well in TRAP2 mice. Those mice that were injected with 4-OHT in a home-cage and later given a seizure served as controls for the dual-labeling experiment. To test whether the observed overlap exceeded the probability of random overlap, we conducted the following calculation (Roy et al., 2017): overlap probability = ((day 2 learning-activated tdT 1ve neurons in a region/Total DAPI 1ve neurons in that region) Â (seizure-activated tdT 1ve neurons in the same region/Total DAPI 1ve neurons in that region)) Â 100%.

Statistical analysis
Data are expressed as means 6 SEM. Graphs were prepared in Adobe Photoshop CC, images were created with BioRender.com, and statistical analysis was done using GraphPad Prism 8.0 software (GraphPad Software Inc.). Two-tailed Student's t tests were used to compare saline and seizure groups. Multiple t tests were used for comparisons among two treatment groups with multiple brain regions under analysis. One-way ANOVA was used for multiple group comparisons. Black asterisks in figures indicate significance by Student's t test. Data were tested for normality before the use of t tests and multiple comparison tests.

Results
We trained activity reporter TRAP2 mice to alternate on a T-maze for a reward to map neuronal ensembles recruited while learning a spatial memory task. On the second day of training, mice acquired the rewarded alteration memory because they demonstrated improved performance on the following and subsequent days ( Fig. 1A; Naik et al., 2021). Therefore, memory encoding occurred on day 2, and memory retrieval occurred on day 3. To tag and identify learning-activated neuronal ensembles, we injected separate cohorts of TRAP2 mice with 4-OHT on training days 1-3, 60 min after training (Fig. 1B). We selected 60 min because c-Fos mRNA expression peaks during that time (Dabrowska et al., 2019). Previous studies find CA1 pyramidal neuron activation during this spatial memory task (Naik et al., 2021). Here, we tested the activation of neurons in the structures outside of the CA1 and extrahippocampal systems that receive projections from the CA1 (Fig. 1C).

Dentate recruitment
We first analyzed spatial memory-activated neurons in the DG that plays a role in pattern separation and forwards the information to the CA1 region of the hippocampus. The acquisition and retrieval of episodic memories engage dentate granule cells (DGCs; Liu et al., 2012;Tonegawa et al., 2015;Josselyn and Tonegawa, 2020).
After T-maze training, we counted tagged DGC neurons manually and by an automated counting method in every slice and found a two-fold increase on day 2 in the number of tagged DGCs than on the previous day ( Fig. 1F,G, day 1: 390.0 6 51.23; day 2: 859.0 6 114.90; day 3: 411.3 6 78.22, one-way ANOVA, n = 6 mice, p = 0.0063). The density of tagged DGCs was higher in anterior dorsal than posterior ventral regions (Fig. 1H). Neurons were clustered near the hinge joining the suprapyramidal and infra-pyramidal blades. There were more tagged DGCs in the supra-pyramidal than in the infrapyramidal blade. The dorsal slices had tagged DGCs in both supra-pyramidal and infra-pyramidal blades, whereas the ventral slices had tagged DGCs only in the supra-blade (Fig. 1F,G). Consistent with the standard model of memory consolidation, day 3 group had fewer labeled DGCs compared with the previous day ( Fig. 1G, p = 0.019). The CA1 pyramidal neuronal recruitment follows a similar trend across training sessions (Naik et al., 2021).

Learning-activated RSC ensembles
The RSC is positioned at the interface between hippocampal formation and cortical sensory regions and receives direct input from CA1 and subiculum (Sugar et al., 2011). This connectivity makes it a candidate for encoding and retrieving spatial memory (Todd et al., 2019;Vann et al., 2009). We found a threefold increase in tagged neuron population on day 2 group compared with day 1 of training ( Fig. 1D,E, day 1, 1069.0 6 146.60, day 2, 3466.0 6 356.20, day 3, 2808.0 6 414.0, one-way ANOVA, n = 6, p , 0.05, F = 11.07). Labeled neurons in RSC were mainly seen in layers 2/3 and 5 of granular and dysgranular regions and crossed the rostrocaudal axis (Fig. 1D). The dorsal part of the RSC had a higher density of tagged cells than the ventral. On day 3, the RSC was as active as observed on day 2 (Fig. 1E).

The PFC and MD
We then assessed whether rewarded alternation learning activated neurons in the dorsal mPFC, which consists of prelimbic (PL), infra-limbic (IL), and ACC and receives direct projections from the hippocampus (Gelinas et al., 2016). Compared with day 1, the fraction of tdTomato1ve neurons in mPFC was significantly more the following day ( Fig. 2A,C, day 1, 508.33 6 30.41, day 2, 3127.20 6 351.23, day 3, 1977.33 6 246.28, one-way ANOVA, n = 6, p , 0.0001, F = 30.19). There were learning-tagged mPFC neuronal ensembles Figure 1. Rewarded alternation task learning triggers the formation of hippocampal and cortical ensembles. A, TRAP2 mice learn to alternate for reward on day 2 and show successful alternation with high accuracy on day 3. B, A schematic of the experimental design for tagging learning activated neuronal ensembles following spatial memory task. C, A hypothetical circuit diagram for spatial memory engram of brain structures with anatomic connectivity and established role in spatial memory. D, Activity-dependent neuronal activation following T-maze training in the RSC on days 1-3 with tdTomato1ve neurons in layers 2/3 and 5 of both granular and dysgranular regions. E, A significantly larger fraction of RSC neurons was tagged on day 2 and many remained active on day 3. F, Learning activated DGCs across training sessions from day 1 to 3. A large fraction of DGCs was tagged on day 2, whereas by day 3, the number dropped. G, Mean number of DGCs labeled on day 2 compared with day 1. H, Location maps show DGCs activation on T-maze training day 2 across the anterior-posterior axis (blue dots represent every tagged DGC). Data presented as mean 6 SEM, ***p , 0.001, ANOVA and post hoc Tukey's test.
Research Article: New Research on day 3 with similar recruitment levels as observed on day 2 ( Fig. 2A,C).
The MD thalamus has extensive reciprocal connections with the mPFC (Mitchell and Chakraborty, 2013), amygdala (Kim and Cho, 2017), and EC (Yang et al., 2019), all of which receive projections from the CA1. Optogenetic inhibition of MD terminals that project to mPFC impairs T-maze task performance (Mitchell and Chakraborty, 2013; Saalmann, 2014; Figure 2. Activation of cortical modules in a working memory paradigm. A, Representative coronal sections showing learningtagged neuronal ensembles in the mPFC, which consists of PL, ILA, and ACC cortices across training days 1-3. There were more tagged neurons on days 2 and 3 compared with day 1. B, Coronal sections through the thalamus with tdTomato (red) and DAPI (blue) staining. On days 2 and 3, learning-triggered neuronal ensembles were evident in MD and LD, which were absent on the first day of training. As a control structure, the PVT was quiescent with no changes across the training sessions. C, Number of tagged neuronal ensembles in mPFC showing a 7-and 5-fold increase on days 2 and 3 compared with day 1. D, MD had many learningtagged neurons on days 2 and 3. E, Labeled neuronal population in PVT remained unchanged across days 1-3. Data are mean 6 SEM, n = 6 per group, **p , 0.01, ****p , 0.0001, ANOVA post hoc Tukey's test. Bolkan et al., 2017). We found very few labeled neurons in MD on day 1 and a larger labeled population on the second day (Fig. 2B,D). Interestingly, during the retrieval on day 3, the MD was as active as on day 2 (Fig. 2B,D). We also studied the PVT as a control region with no role in spatial learning/memory (Fig. 2B,E). We did not find any differences in tdTomato1ve neuronal population in PVT across training days 1-3 of rewarded alternation task (Fig. 2B,E).
Poor recruitment of hippocampal/cortical neuronal ensembles in GluA1 KO mice We then confirmed that neuronal activation during the T-maze task was specific to the spatial working memory neuronal activation and not just because of the locomotion or spatial navigation during the task. The T-maze spatial working memory task is dependent on the plasticity of the GluA1 subunit of AMPA receptors (Sanderson et al., 2008). Mice that lack GluA1 subunit have impaired spatial working memory that is GluA1-dependent, short-term memory. In contrast, these mice show normal spatial reference memory acquisition such as hidden-platform water-maze task that requires GluA1-independent, long-term memory (Reisel et al., 2002;Sanderson et al., 2008).
We tested global GluA1 KO mice for discrete delayed trial rewarded alternation task in T-maze and found that these mice did not learn to alternate. The accuracy rates of successful alternation for these mice were at the chance levels for all training days (see the learning curve; Fig. 3A). The littermate controls showed a typical learning curve for the task with alternation accuracy rates of above 75% by day 3 and stable expression on the following days of probe testing. When tested in the open field apparatus, the GluA1 KO mice were hyperactive independent of sex with no signs of habituation (M4-M7; Fig. 3B). In contrast, control mice habituated shortly after exploring novel environment (M1-M3; Fig. 3B).
We next determined neuronal ensembles recruited following training in GluA1 KO mice. To visualize activated neurons in GluA1 KO mice, we crossed global GluA1 KO with TRAP2 to create GluA1 KO/TRAP2 mice. GluA1 KO/ TRAP mice were injected with 4-OHT 60 min after training on day 2 to tag the activated population following learning in T-maze (Fig. 3F). Visual examination revealed far fewer labeled neurons in these mice on day 2 than the TRAP2 controls (Fig. 3D,E). The hippocampus of GluA1 KO/TRAP2 mice showed significantly low recruitment of neuronal ensembles, with hardly two or three tagged pyramidal neurons in the entire CA1 compared with littermate TRAP2 controls ( Fig. 3C-E). Also, far fewer DGCs were labeled in these mice on the second training day than in controls ( Fig. 3C-E). There was a reduced engagement of RSC neuronal ensembles in GluA1 KO/TRAP2 mice (Fig. 3C-E, n = 6, p , 0.001). The mPFC of these mice was also devoid of any tagged neurons (Fig. 3G). There was poor recruitment of MD compared with controls ( Fig. 3C-E). However, in contrast to this scarce global engagement, the tagged neuronal population in the PVT of GluA1 KO mice was similar to that in controls ( Fig. 3C-E). Thus, failure to learn the alternation rule in GluA1 KO mice was associated with poor recruitment of memory ensembles in these crucial hippocampal and cortical sites otherwise engaged in TRAP2 littermate controls.

Seizure network
Previous studies show that a single generalized tonicclonic seizure immediately after day 2 of training-induced transient retrograde amnesia and impaired T-maze performance accuracy on day 3 (Naik et al., 2021). We tested whether a seizure also engaged memory circuits outside of the CA1 by TRAPing neurons (Fig. 4). Two investigators blinded to the treatment performed unbiased quantification of tagged neurons following either a seizure (PTZ, 40 mg/kg) or saline in TRAP2 mice across brain structures activated during spatial memory formation on day 2 (Fig. 1C). Visual examination showed a larger population of neurons labeled in mice that experienced a seizure than in saline controls (Fig. 4). There were significantly more tagged DGCs in the DG of mice that experienced a seizure (Fig. 4B,E, saline, 790.0 6 119.3, seizure, 3091.0 6 768.70, t test, n = 6 mice each, p = 0.018) than saline controls (Fig. 4A,E). Both infra-blade and supra-blade of the DG of the seizure group mice had many tdT 1ve granule cells across anterior-posterior axis.
We quantified seizure-activated tagged neuronal populations in RSC and mPFC. A single seizure tagged more neurons in the RSC than saline controls (Fig. 4B,D, saline, 1120.0 6 106.20, seizure, 6560.0 6 450.00, n = 6, t test, p , 0.01). Tagged neurons were present in layers 2/3 and 5 in both dorsal and ventral parts of the RSC in seizuretreated mice (Fig. 4B). There were tagged neurons across the anterior-posterior axis that involved both granular and dysgranular parts of the RSC. In addition to RSC, we found significantly more tagged neurons in the mPFC following a single seizure (Fig. 4F). Besides, a qualitative examination showed a large fraction of seizure-tagged neurons in brain regions like the piriform cortex (PIR), MD, EC, amygdala [basolateral amygdala (BLA), basomedial amygdala (BMA)], and others ( Fig. 4B) compared with saline controls (Fig. 3A).

Seizure stimulates spatial memory-activated neurons
We next tested whether seizure-activated neuronal ensembles overlap spatial learning-activated neurons by performing dual-labeling. We first permanently tagged the learning-activated neurons in TRAP2 mice on day 2 of training, and a week later, the same mice received a seizure. The mice were killed 1 h after the seizure to detect Arc-labeled neurons (see the experimental design in Fig.  5A). The control mice were injected with 4-OHT in a home-cage to tag baseline neuronal activation and later received a seizure.
Only 2% of learning-tagged DGCs expressed seizure-activated ARC (Figs. 5B,C, 6A). The control mice showed a similar overlap (Fig. 5B), which could be a result of the basal activity in the DG. We also calculated joint probability of overlap for double labeling of tdTomato and ARC (overlap probability = ((day 2 learning-activated tdT 1ve neurons in a region/Total DAPI 1ve neurons in that region) Â (seizure-activated tdT 1ve neurons in in the same region/Total DAPI 1ve neurons in that region)) Â 100%). The probability of random overlap for DG was 0.0011% [((859/501 240) Â (3091/501 240)) Â 100%; Keller et al., 2018]. We next studied the extra-hippocampal structures, the RSC, and mPFC.
We found a 33% overlap between learning-activated neuronal ensembles and seizure-tagged neurons in the RSC (Fig. 5B,D). We observed co-labeled neurons in layers 2/3 and 5 of the dorsal RSC (Fig. 6B). This overlap in RSC in mice following day 2 training and a seizure a week later was significantly higher than observed in control mice (Fig. 5B). The probability overlap for RSC was 0.23% [((3466/98,148) Â (6560/98,148)) Â 100%; Keller et al., 2018;Murakami et al., 2018]. We also found a 35% overlap among memory and seizure tagged neurons in the mPFC, which was significantly larger than in the control mice (Fig. 5B,E; S1C). In contrast, the probability Figure 3. GluA1 KO/TRAP2 mice fail to learn rewarded alternation and show poor recruitment of memory ensembles. A, Learning curve for mice in the T-maze task shows GluA1 KO/TRAP2 mice fail to learn to alternate for reward even after prolonged training, whereas littermate TRAP2 controls show a normal learning curve. B, Heatmaps of activity-monitoring in open field test (test session of 15 min) in control TRAP2 mice (mouse 1-3) and GluA1 KO mice (mouse 4-7) confirm hyperactivity phenotype in GluA1 KO/TRAP2 mice. C, Quantification shows a significantly fewer number of learning-tagged neurons in GluA1 KO mice in RSC, DG, CA1, mPFC, and MD, but normal engagement in PVT compared with control TRAP2 mice. D, E, Confocal images show representative coronal slices (40 mm) at the level of the hippocampus with tdTomato (red) and NeuN (green) labeling from T-maze day 2 TRAP2 group mice (D) and T-maze day 2 GluA1 KO/TRAP2 group (E). Magnified images show more spatial memory-activated neurons on day 2 in TRAP2 controls than in GluA1 KO/TRAP2 mice in the RSC, hippocampus (CA1, DG), and MD thalamic nucleus. F, A schematic of the experimental design. G, mPFC is not engaged following learning on day 2 of T-maze in mice lacking GluA1. Data are presented as mean 6 SEM, n = 6 mice for each group, **p , 0.01, t test for each pair TRAP2 controls versus GluA1 KO/TRAP2. overlap was 0.070% [((3127/121,000) Â (3255/121,000)) Â 100%; Keller et al., 2018;Murakami et al., 2018]. Interestingly, in the control brain region (PVT) with no apparent spatial memory role, we did not find any overlap among the seizure and learning-activated ensembles (Fig. 5B,F).

Discussion
We find that seizures and spatial learning activate the overlapping set of neurons widely distributed in the brain beyond the hippocampal subfields, causing retrograde amnesia. We mapped for the first-time neurons activated by spatial memory using activity reporter mice and showed their representation beyond the classical hippocampal CA1 in the brain regions like RSC, MD, and mPFC on the cellular level. The T-maze spatial learning task is dependent on the plasticity of the GluA1 subunit of AMPA receptors (Sanderson et al., 2008), so we used GluA1 KO mice crossed to TRAP2 mice to test whether widespread neuronal activation was specific to the spatial memory task and not because of the locomotion or spatial navigation. GluA1 KO/TRAP2 mice failed to learn the spatial memory task despite extended training sessions and lacked tagged neurons in the RSC, MD, and mPFC. A single seizure prevented the recall of alternation memory (Naik et al., 2021). Previous work demonstrated that seizures invade the circuits of memory-activated neurons in the CA1 region of the hippocampus, altering their physiology through saturated LTP (Naik et al., 2021). In this study, seizure-activated and spatial learning-activated ensembles also overlap in the regions outside of the CA1 and cortical structures, causing the failure of hippocampal-cortical interaction and the perturbation of systems consolidation.
We hypothesize that anatomic connectivity, together with neuronal excitability, inhibitory surround, and neuronal plasticity, shape the seizure circuit. The pattern of neuronal activation during seizures might be a stochastic process guided by the increased plasticity of certain cells. Neuronal plasticity level could increase or decrease the probability of cells being recruited by a seizure. Our previous work demonstrated that global removal of the GluA1   (Adotevi et al., 2020). GluA1 KO reduced mortality, severity, and duration of status epilepticus. This indicates that multiple plasticity mechanisms are critical for the occurrence and maintenance of seizures. Seizure-activated neurons are more excitable and have potentiated synapses (Naik et al., 2021). Synaptic connectivity and neuronal excitability optimize the efficacy of memory retrieval, contributing to memory formation and consolidation (Pignatelli et al., 2019). In addition, frontal lobe seizures follow its anatomic projections through the basal ganglia, specifically through the indirect pathway, the neurons of which are more excitable than those of the direct pathway .
The inhibitory tone might also play a role in shaping the seizure circuit. During kindling, neuronal firing and Figure 6. Cortical-neuronal ensembles activated following learning and seizure overlap. Representative images show learningtagged neurons (tdTomato1ve, red) on day 2 of T-maze in (A) DG, (B) mPFC, and (C) RSC. ARC1ve neurons represent neuronal population activated following a single tonic-clonic seizure showed low levels of co-localization in DGCs, whereas mPFC and RSC showed higher levels of overlap between learning-tagged and seizure-tagged neuronal ensembles.
continued design for dual-labeling memory-activated and seizure-activated neurons in the same mouse. Mice that were injected with 4-OHT in homecage and later received a PTZ seizure served as control for dual-labeling experiment. B, Co-localization of memory and seizure activated neurons in DG, RSC, mPFC, and PVT regions with high overlap in mPFC and RSC and no overlap in DG and PVT. C, Memory-activated neuronal ensembles permanently tagged with tdTomato (red) and seizure-activated ARC-immunolabeled neurons (green) in DGC showed low levels of overlap among the ensembles. D, Active neuronal ensembles in the RSC on day 2 (tdTomato) and seizure-tagged neurons (ARC) showed overlap. E, Representative images of mPFC neurons labeled on day 2 of training (tdTomato) and seizure-tagged neurons (ARC) showed a large fraction of tdT 1ve neurons co-labeled for ARC. F, In PVT, day 2 training-activated (tdTomato) and seizure-activated (ARC) neuronal ensembles did not overlap. Data are mean 6 SEM, n = 6 mice for day 2 with seizures, n = 3 for home-cage with seizure, **p , 0.01, t test for each pair.
synchrony are more correlated along the lamellar hippocampal axis that contains excitatory connections than along the septotemporal axis that has inhibitory connections (Ren et al., 2021). Propagation of seizures through the hippocampus is gated by strong inhibition of the DG (Krook-Magnuson et al., 2015). Similarly, in the neocortex, the inhibitory surround shapes the seizure spread (Schevon et al., 2012). Cortical seizures propagate along the neural pathways but with large temporal variability that depends on local inhibitory tone (Wenzel et al., 2017). Also, frontal lobe seizure spread through the basal ganglia, where 95% of the striatal medium spiny neurons are GABAergic 23 . Some seizures required several hundred milliseconds to arrive at the striatum or substantia nigra reticulata, whereas others took seconds to appear in those structures, suggesting stochastic spread based on latency . Internal inhibition within these structures might prevent immediate ictal onset within them.
The standard consolidation model proposes that the hippocampus initially stores memory traces that later fade from the hippocampus and get transferred to the neocortical modules, becoming hippocampus-independent (Frankland and Bontempi, 2005;Tse et al., 2007;Kesner and Rolls, 2015;Squire et al., 2015;Barry and Maguire, 2019). Selective cell activity that begins in the hippocampus entrains other brain regions in the cortex, with seizures disrupting this pattern. Previous studies reported the CA1 engram for spatial memory (Naik et al., 2021). Other studies have also reported that the hippocampal engram in the CA1 maps experiences (Choi et al., 2018;Tanaka, 2021;Tanaka and McHugh, 2018). An important hippocampal subfield, the DG separates patterns and plays a role in memory storage (Liu et al., 2012;Squire et al., 2015;Tonegawa et al., 2015;Josselyn and Tonegawa, 2020). Like CA1, DG activation was transient on day 3, suggesting the memory trace has become hippocampus independent. Hippocampal and neocortical coupling during consolidation is primarily derived from electrophysiological recordings of sharp-wave ripples, spindle oscillations, and u -g coupling present during encoding or sleep (Buzsáki, 2015;Joo and Frank, 2018). However, direct visualization of neuronal ensembles following learning and consolidation remained challenging until the activity-dependent IEG reporter mice were generated (Guenthner et al., 2013;Josselyn and Tonegawa, 2020). We report a map of spatial memory-activated neurons that spans beyond the hippocampal CA1 in the mPFC, MD, and RSC.
The mPFC connects with the dorsal hippocampus and is an important hub in the consolidation process. mPFC role during recent acquisitions and retrieval of remote memories is supported by both human and nonhuman studies (Gelinas et al., 2016;Tanaka, 2021). We find mPFC was already active during early learning (on day 2 of training), contrary to the standard consolidation theory, which proposes delayed recruitment of cortical ensembles. However, recent studies have also found that fear engrams appear early in cortical areas (Josselyn and Tonegawa, 2020). Surprisingly, unlike CA1 and DG, where we observed a drop in tagged neuronal ensembles on day 3, mPFC was as active on day 3 as on the second day of training. The mPFC neuronal ensembles may facilitate successful performance and accuracy in mice on day 3.
We report similar activation dynamics in the RSC, which was significantly engaged on day 2 of training. The RSC, like PFC, was re-engaged on day 3 supporting that cortical modules are active during recent retrieval. Previous studies have shown that in addition to the presence of head direction cells in RSC, it acts as an interface between hippocampal formation and cortical sensory regions, supporting its crucial role for encoding and retrieving spatial memory (Sugar et al., 2011;Buzsáki, 2015). We also found activation of MD nucleus during early learning on day 2 and retrieval on day 3 of the T-maze task. This co-activation is not surprising since MD has reciprocal connections with the mPFC and exhibits minimal connectivity with either motor or sensory pathways (Mitchell and Chakraborty, 2013). Our study supports the role of MD in the neuronal network activated by spatial memory. In contrast to MD, the tagged neuronal population in the PVT remained unchanged with learning, supporting that the tagged neuronal ensembles specifically in structures, i.e., DG, RSC, mPFC, MD, and CA1, are a part of the neuronal network complex for the memory trace. Wide representation of spatial learning-activated neurons reported here is similar to the fear memory studies that demonstrate that aversive learning is also widely distributed in the brain (Joo and Frank, 2018).
Because the T-maze spatial working memory task depends on the GluA1 subunit plasticity of AMPA receptors (Sanderson et al., 2008), we used mice lacking the GluA1 subunit, who have deficits in hippocampus-dependent spatial working memory tasks but intact spatial reference memory (Reisel et al., 2002;Sanderson et al., 2008). We bred GluA1 KO with TRAP2 mice to validate that specific neuronal ensemble represents the memory trace. These mice failed to learn rewarded alternation task and showed a hyperactive-phenotype as reported previously (Reisel et al., 2002;Sanderson et al., 2008). They did not have learning-tagged neuronal ensembles in any learning-associated brain structures, CA1, DG, mPFC, RSC, and MD. However, the PVT region was as active as in the control TRAP2 mice based on the presence of tagged neurons. Although the absence of GluA1 may increase neuronal activation threshold, previous studies determined that global removal of the GluA1 subunit did not affect PTZ seizure threshold or susceptibility (Adotevi et al., 2020). Loss of plasticity mechanisms disrupted the learning and consolidation patterns.
We used the dual-labeling technique and looked for the overlap between memory and seizure-activated neurons. There was a greater overlap in cortical areas, where RSC showed a 33% and mPFC 35% overlap among the memory and seizure-tagged neurons, exceeding the probability levels for double labeling (0.23% and 0.070%, respectively; Keller et al., 2018;Murakami et al., 2018). An overlap of this level in RSC and mPFC suggests that a seizure would erode the learning ensemble information postacquisition, resulting in memory impairment during recent retrieval. However, the hippocampal subfield CA1 showed 11% overlap, whereas only 2% of tagged neurons overlapped in the DG (with 0.0018% and 0.0011% probability overlaps, respectively; Keller et al., 2018). The absence of any overlap in the PVT nucleus, otherwise active during both experiences, also supports its limited role in spatial memory neuronal ensemble.
The DG overlap was less than the overlap in other brain regions. This likely is because of the sparse coding function of the DG (Piatti et al., 2013). Previous studies demonstrated that activity levels in the DG are lower than in the upstream cortical areas (Piatti et al., 2013). DGCs receive inputs from a smaller population in the EC, and its representation is thought to be expanded onto the larger population of DGCs (Piatti et al., 2013). The basal activation of DGCs is estimated to be 1-2% of the total population, similar to our findings (Piatti et al., 2013). Sparse coding in the DG could also be one of the reasons for the lack of overlap between task learning-activated and seizure-activated cells. Neural activity in the DG is sparse both in terms of the proportion of active neurons and their mean firing rates (Piatti et al., 2013). Studies that mapped neuronal activation during novel environment exploration, however, found that the DG exhibited significantly more activated neurons when compared with the familiar (Mazurkiewicz et al., 2022). Therefore, different types of memories can activate different circuits. It is possible that seizures also interfere with encoding a novel environmental context.
In summary, we demonstrate that seizure and memory occupy the same space in the brain. Increased plasticity, more active pathways, or a combination of both, increased the probability of cells being recruited by a seizure. Seizure and memory consolidation pathways overlapped in the neuronal network of extrahippocampal structures, which may explain seizure-induced retrograde amnesia.