Differential vulnerability of neuronal subpopulations of the subiculum in a mouse model for mesial temporal lobe epilepsy

Abstract

Selective loss of inhibitory interneurons (INs) that promotes a shift toward an excitatory predominance may have a critical impact on the generation of epileptic activity. While research on mesial temporal lobe epilepsy (MTLE) has mostly focused on hippocampal changes, including IN loss, the subiculum as the major output region of the hippocampal formation has received less attention. The subiculum has been shown to occupy a key position in the epileptic network, but data on cellular alterations are controversial. Using the intrahippocampal kainate (KA) mouse model for MTLE, which recapitulates main features of human MTLE such as unilateral hippocampal sclerosis and granule cell dispersion, we identified cell loss in the subiculum and quantified changes in specific IN subpopulations along its dorso-ventral axis. We performed intrahippocampal recordings, FluoroJade C-staining for degenerating neurons shortly after status epilepticus (SE), fluorescence in situ hybridization for glutamic acid decarboxylase (Gad) 67 mRNA and immunohistochemistry for neuronal nuclei (NeuN), parvalbumin (PV), calretinin (CR) and neuropeptide Y (NPY) at 21 days after KA. We observed remarkable cell loss in the ipsilateral subiculum shortly after SE, reflected in lowered density of NeuN+ cells in the chronic stage when epileptic activity occurred in the subiculum concomitantly with the hippocampus. In addition, we show a position-dependent reduction of Gad67-expressing INs by ∼50% (along the dorso-ventral as well as transverse axis of the subiculum). This particularly affected the PV- and to a lesser extent CR-expressing INs. The density of NPY-positive neurons was increased, but the double-labeling for Gad67 mRNA expression revealed that an upregulation or de novo expression of NPY in non-GABAergic cells with a concomitant reduction of NPY-positive INs underlies this observation. Our data suggest a position- and cell type-specific vulnerability of subicular INs in MTLE, which might contribute to hyperexcitability of the subiculum, reflected in epileptic activity.

1. Introduction

Mesial temporal lobe epilepsy (MTLE) is among the most frequent forms of pharmacoresistant epilepsies. To elucidate the underlying pathophysiological mechanisms, the hippocampus as site of origin of focal seizures has been intensely investigated. This revealed histological alterations including hippocampal sclerosis, granule cell dispersion and mossy fiber sprouting as hallmarks of MTLE in patients and animal models (Sutula et al., 1989; Houser, 1990; Sloviter, 1994; Suzuki et al., 1995; Haas et al., 2002; Thom et al., 2002; Blümcke et al., 2013; Häussler et al., 2016). In contrast, the subiculum received less attention, despite playing a pivotal role in the network as a major output region of the hippocampal formation.

With its large proportion of bursting cells (Taube, 1993; Greene and Totterdell, 1997) and a network that facilitates synchronization of neuronal activity (Behr and Heinemann, 1996; Fiske et al., 2020), the subiculum is an ideal candidate region for the generation of epileptic activity. Moreover, the subiculum becomes subject of considerable deafferentation due to extensive cell loss in the CA1 region in patients (Blümcke et al., 2013) and animal models (Suzuki et al., 1995; Janz et al., 2017a) which might induce network imbalance.

Indeed, using acute brain slices from MTLE patients undergoing resective surgery, it has been proposed that the subiculum is the origin of spontaneous inter-ictal activity (Cohen et al., 2002; Wozny et al., 2005), which was confirmed by in vivo recordings of patients with pharmacoresistant MTLE (Fabó et al., 2008). Interestingly, the highest probability for generation of paroxysmal high frequency oscillations (HFOs) switches from the hippocampus to the subiculum with increasing grades of hippocampal sclerosis (Tóth et al., 2021). On the histological level, pathological changes in the subiculum are less evident than those seen the sclerotic hippocampus: Overall neuronal numbers are only mildly affected in patient tissue (Arellano et al., 2004; Andrioli et al., 2007), but data on parvalbumin (PV)-expressing interneurons (INs) reach from no loss (Arellano et al., 2004) to more than 50% loss (Andrioli et al., 2007). Functionally, epileptic activity in the subiculum has been suggested to arise from impaired inhibition (Cohen et al., 2002).

Studies in animal models for MTLE have shown decreased overall neuronal density (Knopp et al., 2005; Drexel et al., 2012b) and loss of several IN populations [PV, calretinin (CR), somatostatin (SOM), and others; (He et al., 2010; Drexel et al., 2011, 2012a)]. However, the induction of MTLE by the systemic injection of a convulsant drug causes damage to the hippocampus and subiculum of both hemispheres and does not exactly reproduce the large patient cohort with unilateral hippocampal sclerosis (Blümcke et al., 2013).

To investigate changes in the subiculum of both hemispheres under conditions of severe focal, unilateral hippocampal sclerosis, we used the intrahippocampal kainate (ihKA) mouse model which reproduces the human characteristics particularly well (Suzuki et al., 1995; Depaulis and Hamelin, 2015). Our previous studies revealed that despite the focal KA injection, epileptic activity propagates to different subregions of the hippocampus proper, along large parts of the dorso-ventral (DV) axis, and to the contralateral hippocampus and the entorhinal cortex with different intensity (Froriep et al., 2012; Häussler et al., 2012; Janz et al., 2017b; Kilias et al., 2018, 2023; Tulke et al., 2019). It remains, however, unclear whether epileptic activity also affects the subiculum and whether this is associated with histological changes in this animal model.

Here, we show that not only generalized seizures but also subclinical epileptic activity affects the subiculum and is associated with neuronal loss, in particular loss of different IN populations to type- and position-specific degrees, depending on the hemisphere, the DV level and the position along the transverse axis of the subiculum.

2. Materials and methods

2.1. Animals

Experiments were performed in adult (>8 weeks) male negative littermates of transgenic mouse lines bred on C57Bl/6 background [B6(SJL)-Tg(Amigo2-icre/ERT2)1Ehs, C57BL/6-Tg(Thy1-eGFP)-M-Line, C57BL/6-Tg(Rbp4-cre)KL100Gsat or C57BL/6-Tg(Rbp4-cre)KL100Gsat x Ai32(RCL-ChR2)(H134R)/EYFP; CEMT Freiburg]. All lines have been shown to develop the full pattern of histological and electrophysiological characteristics of the ihKA model in a comparable manner in our earlier studies (Paschen et al., 2020; Kilias et al., 2023, Paschen/Kleis personal communication, own observations). Mice were housed under standard conditions [12 h-light/dark, room temperature (RT), food and water ad libitum]. Experiments were carried out following the guidelines of the European Community’s Council Directive of 22 September 2010 (2010/63/EU), approved by the regional council (Regierungspräsidium Freiburg). In total, n = 36 mice were used (n_NaCl = 14, n_KA = 22).

2.2. Surgery and in vivo recordings

Surgery was performed under deep anesthesia [ketamine hydrochloride 100 mg/kg body weight, xylazine 5 mg/kg, atropine 0.1 mg/kg, intraperitoneal (i.p.) injection] with analgesic treatment [buprenorphine hydrochloride 0.1 mg/kg, carprofen 5 mg/kg, subcutaneous (s.c.) injection; lidocaine 2 mg/kg, s.c. to scalp]. Mice were mounted in a stereotaxic frame and injections of 50 nL KA (20 mM in 0.9% NaCl, Tocris, Bristol, UK) or 0.9% NaCl into the right dorsal hippocampus were performed using a Nanoject III (Drummond Scientific Company, Broomall, PA, USA) at coordinates anteroposterior (AP) = −2.0 mm, mediolateral (ML) = −1.4 mm and DV = −1.8 mm from bregma. In KA-injected mice, behavioral status epilepticus (SE) with head nodding, convulsions, rotations and immobility was verified as previously described (Tulke et al., 2019). Postsurgical analgesia was performed with carprofen for 2 days.

One week after ihKA/ihNaCl injection, a subgroup of mice (n_NaCl = 5, n_KA = 6) was implanted with wire electrodes [platinum-iridium, teflon-insulated, ∅125 μm, World Precision Instruments (WPI), Sarasota, FL, USA] into the dentate gyrus (DG; AP = −2.0 mm, ML = ± 1.4 mm, DV = −1.9 mm) and the subiculum (AP = −3.0 mm, ML = ± 1.4 mm, DV = −1.4 mm) of one or both hemispheres (anesthesia see above, analgesia with buprenorphine+carprofen for 2 days). Jewelers’ screws above the frontal cortex served as ground/reference. A connector was permanently fixed to the skull with cyanoacrylate and dental cement. Freely moving mice were recorded (6−8 sessions, ∼2 h each) between 14 and 33 days after ihKA connected to a miniature preamplifier (Multi Channel Systems, Reutlingen, Germany). Signals were amplified (500-fold, band-pass filter 1 Hz–5 kHz) and digitized [sampling rate 10 kHz, Spike2 software; Cambridge Electronic Design (CED), Cambridge, UK]. Mice that died after surgery, with misplaced electrodes or with unsuccessful kainate injection (verified by post-hoc Nissl staining in frontal brain slices) were excluded (n_NaCl = 2, n_KA = 3).

2.3. Perfusion and slice preparation

Mice were transcardially perfused under deep anesthesia with 0.9% saline followed by paraformaldehyde [PFA, 4% in 0.1 M phosphate buffer (PB), pH 7.4]. Brains were post-fixed in PFA overnight (4°C), cryoprotected (30% sucrose) and frozen in isopentane. Horizontal sections (50 μm) were prepared on a cryostat and collected in PB for immunohistochemistry or in 2× saline sodium citrate (SSC) for fluorescence in situ hybridization (FISH).

2.4. FluoroJade C

FluoroJade C (FJC, Merck Millipore, Darmstadt, Germany) staining was used to monitor cell death 2 and 4 days after ihKA injection. Sections were mounted on gelatin-coated slides and dried overnight at RT, re-hydrated (1% NaOH in 100% Ethanol for 5 min; 70% Ethanol for 2 min, distilled water for 2 min), incubated in 0.06% potassium permanganate (10 min), rinsed and transferred in 0.0001% FJC solution (10 min, darkness), followed by washing, xylene immersion, and coverslipping with Hypermount.

2.5. Immunohistochemistry

After preincubation with 0.25% Triton X−100 and 10% normal horse/goat serum in 0.1 M PB (30 min, free-floating), sections were incubated overnight (4°C) with primary antibodies: guinea pig anti-NeuN (1:500; #266004, Synaptic Systems, Göttingen, Germany), rabbit anti-calretinin (CR; 1:500; #7699/4), rabbit anti-parvalbumin (PV; 1:1000; #PV27, all Swant, Burgdorf, Switzerland), and rabbit anti-NPY (1:1000; #ab30914, Abcam, Cambridge, UK). Secondary antibodies, marked with Cy2, Cy5 (both 1:200; Jackson ImmunoResearch, West Grove, PA, USA) or Alexa555 (1:500; Thermo Fisher Scientific, Waltham, MA, USA) were applied (2−3 h, RT). Sections were counterstained with 4’,6-diamidino-2-phenylindole (DAPI, 1:10,000). When combining FISH and immunohistochemistry, Triton X-100 was omitted. Sections were coverslipped with fluorescence mounting medium (Immu-Mount, Thermo Shandon Ltd., Schwerte, Germany or Dako, Hamburg, Germany).

2.6. Fluorescence in situ hybridization

Expression of Gad67 mRNA was investigated by FISH with digoxigenin (DIG)-labeled cRNA probes. Probes were generated by in vitro transcription from appropriate plasmids (Kulik et al., 2003). Sections were hybridized overnight at 45°C, followed by detection with an anti-DIG antibody and tyramide signal amplification as described (Marx et al., 2013; Tulke et al., 2019).

2.7. Microscopy and quantification

Sections were imaged with an epifluorescence microscope (AxioImager 2, Zen software, Zeiss, Göttingen, Germany; 10-fold Plan Apochromat objective, NA = 0.45; AxioCam 506). Cell counting was performed with ImageJ (version 1.52a; NIH) on all available brain sections from DV −1.2 to −3.2 mm relative to bregma. The region of interest (ROI, principal cell layers of the subiculum) was drawn according to a brain atlas (Franklin and Paxinos, 2008). Cells were counted manually using the cell counter plugin. Criteria for selection were: signal intensity clearly distinguishable from background, cell shape and presence of a DAPI-stained nucleus. Blinding to the treatment was impossible due to KA-induced alterations in the hippocampus. To compare control and KA-injected mice along the DV axis, sections were assigned to the following subgroups: dorsal to the injection site level (−1.2 to −1.8 mm ventral to bregma–usually three sections/mouse), approximately at the injection site (−1.9 to −2.1 mm–2 sections/mouse) and ventral (−2.2 to maximally −3.2 mm – 4−6 sections/mouse). Cell density was calculated for the ROI (cells/mm²) and averaged within subgroups per animal. Values for the ipsi- and contralateral subiculum of NaCl-injected mice were grouped since we did not observe any differences. FluoroJade C-positive cells were detected in sections at the injection site as regional maxima in an intensity landscape using the image processing toolbox in Matlab (R2020b, The MathWorks Inc., Natick, MA, USA). A rectangular window spanning from the alveus to the molecular layer of 200 μm width was used as region of interest and divided into 20 equally sized bins to normalize for differences in width of the subiculum between individual sections. Counts for every bin were averaged across mice.

2.8. Statistics

Statistical analysis was performed with GraphPad Prism (version 9.3.1, GraphPad Software, Boston, MA, USA). After confirmation of normality with a Shapiro-Wilk test, a one-way analysis of variance (ANOVA) followed by Tukey’s multiple comparisons test was performed, otherwise a Kruskal–Wallis test was followed by Dunn’s multiple comparisons test.

3. Results

3.1. Epileptic activity occurs in the dentate gyrus and in the subiculum

In control mice (n_NaCl = 3), local field potential (LFP) patterns during immobility consisted of irregular activity and occasional dentate spikes in the DG and sharp-wave ripples (SWR) in the subiculum (Figure 1A). In ihKA mice (n_KA = 3) we recorded recurrent epileptic activity in the DG, as described previously (Häussler et al., 2012; Tulke et al., 2019), but also in the subiculum (Figure 1B). Fast ripples superimposed on large population spikes were evident in the DG (Figure 1B1) and in the subiculum where they resembled SWRs (Figure 1B2). In 2/3 mice generalized convulsive seizures occurred during the recording sessions which involved the ipsilateral DG and subiculum, and, in addition, the contralateral side where large epileptic population spikes were measurable (Figures 1C, C1). In agreement with previous observations (Lisgaras and Scharfman, 2022), it was impossible to determine the exact onset site whereas the post-ictal depression phase started with a delay on the contralateral side (Figure 1C2).

FIGURE 1

3.2. Position- and layer-specific neuronal loss in the subiculum

To determine whether the subiculum is structurally altered in chronic MTLE, we performed immunohistochemistry for NeuN at 21 days after ihKA (Figure 1, n_NaCl = 7 mice, n_KA = 7 mice). We will refer to the part of the subiculum adjoining CA1 as proximal and the area adjacent to the presubiculum as distal and regard the pyramidal cell layer along its radial axis using the alveus and the cell-poor molecular layer as borders, according to the brain atlas (Franklin and Paxinos, 2008). We compared three levels along the DV axis in horizontal sections: dorsal to the injection site (Figure 1D), around the injection site (Figure 1E) and ventral to the injection site (Figure 1F).

Extensive cell loss in areas CA1 and CA3 of the hippocampus at the dorsal and the injection site level, accompanied by salient granule cell dispersion at the level of the injection site confirmed the successful KA injection (Figures 1J, K), as described previously (Suzuki et al., 1995; Heinrich et al., 2006; Marx et al., 2013). Regarding the dorsal level, a reduced density of NeuN⁺ cells was visible in proximal and distal parts of the subiculum (Figure 1J) in comparison to controls (Figure 1G). At the level of the injection site NeuN⁺ cells were also thinned out in the proximal and distal ipsilateral subiculum (Figures 1H, P, K, Q), mainly in the lower half adjacent to the molecular layer (Figure 1Q). Moreover, many of the remaining neurons in the subiculum in both hemispheres displayed weaker NeuN staining than controls (Figures 1P-R). The latter was also visible in the ventral subiculum, yet without any clear sign of cell loss (Figure 1L) compared to control (Figure 1I). In the contralateral subiculum neuronal density was comparable to controls at all levels (Figures 1M-O, R).

To determine whether the reduced density was due to SE-induced neurodegeneration, we performed FJC staining at 2 and 4 days after ihKA injection. At 2 days, degenerating neurons were densely arranged in ipsilateral CA1 (Figures 2A, B) and all over the ipsilateral hippocampus (data not shown), as described previously (Marx et al., 2013). In the subiculum, FJC⁺ cells were abundant throughout the transverse axis of the pyramidal cell layer at the dorsal position, whereas at the injection site and the ventral position, they tended to be denser in the lower half located adjacent to the molecular layer than close to the alveus (Figures 2A-C, n_KA = 2 mice). At 4 days after ihKA, FJC⁺ cells showed a similar distribution, but with overall lower density (Figures 2D-F, mean spatially resolved cell density at the level of the injection site in Figure 2E1, n_KA = 4 mice). On the contralateral side there were only very few FJC⁺ cells both in the hippocampus and the subiculum at all levels for both time points (Figures 2G-I, 2 days not shown).

FIGURE 2

Together, early cell loss and reduced neuronal density in the chronic stage confirm ihKA-induced structural changes in the ipsilateral subiculum, yet much less prominent than in the hippocampus. Next, we focused on selected IN populations.

3.3. Loss of GABAergic interneurons in the subiculum

In controls, GABAergic INs, as displayed by FISH for Gad67 mRNA, were densely arranged throughout the pyramidal cell layer of the proximal and distal subiculum (Figure 2J-L, n_NaCl = 4 mice). Quantitative analysis did not reveal any major differences between the dorsal subiculum and the level of the injection site (Figures 2J-L, S, T) but in the ventral subiculum, the density of Gad67 expressing cells was slightly higher (Figure 2U).

At 21 days after ihKA (n_KA = 4 mice), Gad67 mRNA⁺ INs were strongly diminished throughout the ipsilateral dorsal subiculum (Figure 2M) and at the level of the injection site (Figure 2N), and in some mice they still appeared less dense further ventrally (Figure 2O). In addition, the Gad67 FISH was less intense indicating reduced mRNA levels per cell. In the contralateral subiculum, the density of Gad67⁺ cells was comparable to control at all levels but the staining intensity was also lower (Figures 2P-R). Quantification in the ipsilateral subiculum revealed the loss of ∼50% of Gad67⁺ INs at the dorsal position and at the level of the injection site (Figures 2S, T; dorsal: one-way ANOVA p = 0.002, Tukey’s multiple comparison test NaCl-KAi p = 0.002, NaCl-KAc p = 0.453, KAi-KAc p = 0.008; injection site: one-way ANOVA p = 0.024, Tukey’s multiple comparison test NaCl-KAi p = 0.029, NaCl-KAc p = 0.906, KAi-KAc p = 0.057). Ventrally, mean values did not differ from controls but individuals ranged from a strong loss to normal density (Figure 2U; one-way ANOVA p = 0.305).

3.4. Loss of parvalbumin- and calretinin-positive interneurons in the ipsilateral subiculum

In control mice (n_NaCl = 5), INs expressing the Ca²⁺-binding protein PV were scattered throughout the pyramidal cell layer of the subiculum with their fibers building a honeycomb-like pattern reaching into the molecular layer (Figures 3A-C). Their density was comparable at all DV positions (Figures 3J-L).

FIGURE 3

At 21 days after ihKA (n_KA = 6 mice), we observed a reduction of PV⁺ INs at all DV positions of the ipsilateral subiculum (Figures 3D-F), but most pronounced dorsally (Figure 3D) with PV⁺ cells located close to the molecular layer being slightly more affected than those close to the alveus (Figures 3D-F). Quantitative analysis confirmed a significant reduction of PV⁺ INs in the dorsal ipsilateral subiculum by ∼40% (Figure 3J; one-way ANOVA p = 0.0002, Tukey’s multiple comparison test NaCl-KAi p = 0.0008, NaCl-KAc p = 0.955, KAi-KAc p = 0.0004) and at the injection site by ∼36% (Figure 3K; one-way ANOVA p = 0.026, Tukey’s multiple comparison test NaCl-KAi p = 0.025, NaCl-KAc p = 0.634, KAi-KAc p = 0.113). The reduction of mean density in the ventral subiculum did not reach significance (Figure 3L; one-way ANOVA p = 0.051). On the contralateral side, the distribution and density of PV⁺ cells were comparable to control at all positions (Figures 3G-L).

Next, we analyzed the distribution of cells expressing the Ca²⁺-binding protein CR. These INs were much sparser than PV⁺ INs and scattered throughout the pyramidal cell layer of all DV levels in controls (n_NaCl = 5 mice, Figures 4A-C).

FIGURE 4

At 21 days after ihKA we did not observe any major changes in the dorsal subiculum of the ipsi- or contralateral hemisphere, respectively, (n_KA = 6 mice; Figures 4D, G, J; one-way ANOVA p = 0.138). In contrast, at the level of the injection site the density of CR-expressing cells was reduced throughout the pyramidal cell layer of the ipsilateral (Figures 4E, H) and contralateral subiculum (Figures 4E, H, K, one-way ANOVA p = 0.007, Tukey’s multiple comparison test NaCl-KAi p = 0.010, NaCl-KAc p = 0.012, KAi-KAc p = 0.996). In ventral parts, CR⁺ cells were reduced in some mice in both hemispheres (Figures 4F, I), whereas others showed normal distribution of CR⁺ cells. This is reflected in high variability in the quantification, but no significant differences (Figure 4L, one-way ANOVA p = 0.123).

3.5. Reduction of neuropeptide Y-expressing interneurons but compensatory NPY upregulation

In control mice (n_NaCl = 5), NPY⁺ neurons were sparsely distributed throughout the pyramidal cell layer mainly of the proximal and to a much lesser extent the distal subiculum (Figure 5A) with increasing mean density from dorsal to ventral (Figures 5S1-U1). Their intense arborizations extended into the molecular layer (Figures 5B, C).

FIGURE 5

At 21 days after ihKA (n_KA = 6 mice), NPY staining was strongly increased in most mice reflected by a much higher density of NPY⁺ somata and axonal processes at all DV and proximo-distal levels of the ipsilateral (Figures 5G-I) and contralateral subiculum (Figures 5M-O). Shapes, sizes and the staining intensity of NPY⁺ cells appeared very diverse after KA injection (Figure 5H1). Quantitative analysis revealed a high variability in the density of NPY⁺ cells ranging from unchanged density to a fourfold increase along the entire DV axis of the ipsilateral subiculum, but the means were not significantly different (Figures 5S1-U1; one-way ANOVA dorsal p = 0.068, injection site p = 0.163, ventral p = 0.174). A similar variability with strongly elevated density in some but not all mice was also observable in the contralateral subiculum.

Is there an increase of NPY-expressing INs in the subiculum? In our previous study regarding the hippocampus proper we found a transient upregulation in excitatory cells alongside a decreased number of NPY⁺ INs (Marx et al., 2013). To test this in the subiculum, we combined NPY immunohistochemistry with FISH for Gad67 mRNA in a subset of samples. In controls, the majority of NPY⁺ neurons were also Gad67-positive (n_NaCl = 3 mice; Figures 5D-F, E1), but a small number of NPY⁺/Gad67^– cells (Figures 5E, E1) explains the overall smaller numbers of double-stained cells compared to NPY alone in the quantification of controls (Figures 5S2-U2).

In contrast, at 21 days after ihKA, a subset of the NPY⁺ cells co-expressed Gad67 mRNA, whereas most were negative (NPY⁺/Gad67^–), indicating the upregulation of NPY in non-GABAergic cells of both hemispheres (n_KA = 3 mice, Figures 5J-L, P-R, K1, Q1). In the contralateral subiculum, the density and distribution of NPY⁺/Gad67⁺ double-positive neurons were comparable to controls (Figures 5S2-U2). In the ipsilateral subiculum, quantification of NPY⁺/Gad67⁺ cells revealed a tendency toward lower numbers of double-labeled cells which was significant at the level of the injection site (Figures 5S2-U2; dorsal one-way ANOVA p = 0.441, injection site Kruskal-Wallis test p = 0.011, Dunn’s multiple comparison test NaCl-KAi p = 0.034, NaCl-KAc p = 0.408, KAi-KAc p = 0.890, ventral one-way ANOVA p = 0.184). Together our data show a loss of NPY⁺ INs and NPY upregulation in the subiculum in non-GABAergic cells.

4. Discussion

Our study provides a comprehensive insight into the pathological changes of the subiculum in a focal MTLE model with unilateral hippocampal sclerosis. Degenerating neurons in the ipsilateral subiculum were detected already 2 days after ihKA and neuronal density was persistently lower in chronic epilepsy. Specifically, we found the loss of GABAergic cells, mainly of PV⁺ INs and to a lesser extent of CR⁺ and NPY⁺ INs with location-dependent patterns for the individual populations. In contrast, NPY was upregulated in non-GABAergic cells throughout the subiculum of both hemispheres. Together, our data show that in the ihKA mouse model, the subiculum is involved in epileptic activity and undergoes epilepsy-associated histological changes of similar nature but much less pronounced than in the dorsal ipsilateral hippocampus.

The degeneration of neurons early after ihKA injection is in agreement with the hypothesis that cell death in the subiculum contributes to epileptogenesis in the network beyond the hippocampus. Neuronal degeneration was nearly equally pronounced at the dorsal position and at the level of the injection site. With increasing distance from the injection site, the subiculum was better preserved in agreement with our previous studies investigating the hippocampus (Häussler et al., 2012; Marx et al., 2013; Janz et al., 2017b). Interestingly, the part of the pyramidal cell layer closer to the molecular layer was more affected than that close to the alveus, in particular for GAD- and PV-positive INs, suggesting a position-dependent differential vulnerability. Indeed, expression of molecular markers, intrinsic firing patterns and in- and outward connectivity of principal cells in the subiculum differ along the proximo-distal as well as the transverse axis (Cembrowski et al., 2018a,b), which might render them differentially vulnerable. A similar position-dependent diversity of neuronal properties among IN populations is also conceivable. Otherwise, different interaction patterns with the individual principal cell populations in the individual layers might also lead to selective vulnerability. Finally, external inputs, e.g., those expressing the vesicular glutamate transporter 2 (vGlut2) and Zink transporter 3 (ZnT3) are not equally distributed along the DV and transversal axes of the subiculum (Ishihara and Fukuda, 2016) and might lead to different degrees of overexcitation during status epilepticus and in chronic epilepsy.

In comparison to the hippocampus, the distribution of IN populations in the subiculum is less well described. We therefore chose Gad67 as a general IN marker and selected three markers of specific populations to first map their distribution in controls under our experimental conditions (Figure 6). INs expressing PV were the largest fraction with somata throughout the pyramidal cell layer and processes extending into the molecular layer. We found PV⁺ INs distributed nearly equally in the proximal and distal subiculum, which is in accordance with previous work (Fujise et al., 1995). CR⁺ cells were much sparser and scattered throughout the pyramidal cell layer of the proximal and distal subiculum. Cells expressing NPY, of which most were GABAergic in controls, showed a gradient from dorsal to ventral and were more numerous in the proximal than in the distal subiculum. To complement the analysis of IN populations expressing specific Ca²⁺-binding proteins, we initially also performed immunostaining for the calbindin, which is not only a marker for INs, but also for subicular principal cell populations (Fujise et al., 1995; Ishihara and Fukuda, 2016). However, the staining patterns under epileptic conditions were too diffuse for analysis and therefore discarded. For the present study, we did not map INs expressing SOM, cholecystokinin (CCK) or vasoactive intestinal peptide (VIP), which based on the density of GAD⁺ cells might comprise around 25−30% of INs and will be of relevance in future studies.

FIGURE 6

In ihKA mice, PV⁺ INs were lost in the dorsal subiculum and at the injection site (40 and 36% reduction, respectively), which is similar or eventually slightly less than in systemic rat and mouse MTLE models [30−60% loss, (Knopp et al., 2005; He et al., 2010; Drexel et al., 2011)]. In contrast, at the ventral position, the density of PV⁺ INs varied between normal and reduced density which is in agreement with observations we made earlier in the hippocampus (Marx et al., 2013) and might depend on the intensity of SE and the DV extent of hippocampal sclerosis. This could eventually serve as an explanation why findings on PV⁺ INs in the subiculum of humans are controversial (Arellano et al., 2004; Andrioli et al., 2007), as this might depend on the degree of sclerosis and the localization of the analyzed tissue along the hippocampal length axis. In accordance with the unilateral injection and unchanged numbers of Gad67 mRNA-expressing cells, there was no major loss of PV⁺ INs in the contralateral subiculum, suggesting that generalized seizures that also affect the contralateral subiculum as shown by our intracranial recordings, do not lead to major IN loss. It is, thus, more likely that ipsilateral damage associated with direct cytotoxic effects of KA is detrimental to the ipsilateral PV⁺ cell population. For CR⁺ INs we found a reduction at the injection site whereas at the dorsal site their density was nearly normal, which differed from Gad67⁺ and PV⁺ populations. Generally, controversial data on CR exist for the human subiculum in MTLE ranging from a reduction (Maglóczky et al., 2000) to normal density (Blümcke et al., 1996), and high variability was also evident in systemic rat models (Knopp et al., 2008; Drexel et al., 2011). A solely local effect focused to the injection site is an unlikely explanation, as the reduction also affects ventral parts in some mice. Instead, the heterogeneous nature of CR-expressing neurons might play a role and require more elaborate investigation of subpopulation vulnerability with different double labeling approaches along the DV axis. Interestingly, the reduction also affected the contralateral subiculum despite no or only little detectable cell loss and no reduction of Gad67⁺ cells. This suggests that not necessarily a loss of CR⁺ cells, but also position- or neuron type-dependent downregulation of CR might play a role. To investigate this, fate mapping of those neurons in mouse lines with CR-dependent expression of Cre recombinase will be necessary to determine whether CR⁺ cells are lost or CR expression is downregulated, and to find out, whether this is associated with altered Ca²⁺ buffering capacity or might be compensated by upregulation of other Ca²⁺-binding proteins. Double-labeling for NPY and Gad67 mRNA revealed that these INs were also vulnerable to the ihKA injection, albeit less than PV⁺ INs in agreement with observations in the hippocampus (Marx et al., 2013). As mentioned earlier, we did not analyze populations expressing SOM, VIP or CCK, but the strong loss of Gad67 mRNA-expressing cells in the dorsal subiculum suggests also reductions for these populations.

How might the reduced IN density affect population activity in the subiculum? It has been shown that during SWRs in the LFP subicular pyramidal cells either increase or decrease their firing and that these divergent cell types correspond to bursting and regular firing neuron types, respectively (Böhm et al., 2015). Interestingly, regular firing cells receive stronger synaptic inhibition during SWRs, in particular from PV⁺ INs which themselves are strongly activated during SWRs. Given the reduction in PV⁺ (and other) INs after ihKA injection, a shift toward more cells in bursting mode and hence from physiological SWRs to pathologic epileptic bursts is conceivable. In addition, despite the loss of CA1 input to the subiculum, projections from CA2 persist (own unpublished observations) and inputs from the entorhinal cortex, which strongly target the subiculum (Steward, 1976; Cembrowski et al., 2018a) are also very likely to persist as we have shown it for the DG (Janz et al., 2017b). In coincidence with reduced feed-forward inhibition these excitatory projections might strongly activate the subiculum. Such high excitatory pressure might also cause compensatory mechanisms to restore network balance. Indeed, we observed the upregulation of NPY after ihKA, comparable to epileptic rats (Drexel et al., 2012a). Various in vitro and in vivo studies in animal models (Woldbye et al., 1996; Klapstein and Colmers, 1997), as well as in human tissue (Patrylo et al., 1999; Ledri et al., 2015) have shown NPY acting as endogenous anti-convulsant in the hippocampal network. The upregulation did not only affect the ipsilateral, but also the contralateral subiculum, which presents with no or only little cell loss. This indicates that whereas cell loss might be an early effect due to direct effects of KA or cell-intrinsic processes triggered by hyperexcitation during status epilepticus, compensatory processes might apply to protect neurons in both hemispheres from generalized seizures or to reduce overall excitability. Determining whether subicular NPY plays a role in seizure suppression in the subiculum will require experiments with local knock-out of NPY or its receptors.

In summary, our data point toward an important role of the subiculum in the epileptic network by undergoing a 50% loss of INs and compensatory processes. The fact that the subiculum shows changes similar in nature as the hippocampus, but less intense, might render it into a hub which is disinhibited enough, and on the other hand presents sufficient neuronal preservation to contribute to the generation and/or the propagation of epileptic activity from the hippocampus to other parts of the brain.

References

• 1

  AndrioliA.Alonso-NanclaresL.ArellanoJ. I.DeFelipeJ. (2007). Quantitative analysis of parvalbumin-immunoreactive cells in the human epileptic hippocampus.Neuroscience149131–143. 10.1016/j.neuroscience.2007.07.029

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  ArellanoJ. I.MuñozA.Ballesteros-YáñezI.SolaR. G.DeFelipeJ. (2004). Histopathology and reorganization of chandelier cells in the human epileptic sclerotic hippocampus.Brain12745–64. 10.1093/brain/awh004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  BehrJ.HeinemannU. (1996). Low Mg2+ induced epileptiform activity in the subiculum before and after disconnection from rat hippocampal and entorhinal cortex slices.Neurosci. Lett.20525–28. 10.1016/0304-3940(96)12360-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  BlümckeI.BeckH.NitschR.EickhoffC.SchefflerB.CelioM. R.et al (1996). Preservation of calretinin-immunoreactive neurons in the hippocampus of epilepsy patients with Ammon’s horn sclerosis.J. Neuropathol. Exp. Neurol.55329–341.

  • Google Scholar

• 5

  BlümckeI.ThomM.AronicaE.ArmstrongD. D.BartolomeiF.BernasconiA. (2013). International consensus classification of hippocampal sclerosis in temporal lobe epilepsy: a task force report from the ILAE commission on Diagnostic methods.Epilepsia541315–1329.

  • Google Scholar

• 6

  BöhmC.PengY.MaierN.WintererJ.PouletJ. F. A.GeigerJ. R. P.et al (2015). Functional diversity of subicular principal cells during hippocampal ripples.J. Neurosci.3513608–13618. 10.1523/JNEUROSCI.5034-14.2015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  CembrowskiM. S.PhillipsM. G.DiLisioS. F.ShieldsB. C.WinnubstJ.ChandrashekarJ.et al (2018a). Dissociable structural and functional hippocampal outputs via distinct subiculum cell classes.Cell1731280–1292.e18.

  • Google Scholar

• 8

  CembrowskiM. S.WangL.LemireA. L.CopelandM.DiLisioS. F.ClementsJ.et al (2018b). The subiculum is a patchwork of discrete subregions.elife7:e37701. 10.7554/eLife.37701

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  CohenI.NavarroV.ClemenceauS.BaulacM.MilesR. (2002). On the origin of interictal activity in human temporal lobe epilepsy in vitro.Science2981418–1421.

  • Google Scholar

• 10

  DepaulisA.HamelinS. (2015). Animal models for mesiotemporal lobe epilepsy: the end of a misunderstanding?Rev. Neurol.171217–226. 10.1016/j.neurol.2015.01.558

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  DrexelM.KirchmairE.Wieselthaler-HölzlA.PreidtA. P.SperkG. (2012a). Somatostatin and neuropeptide Y neurons undergo different plasticity in parahippocampal regions in kainic acid-induced epilepsy.J. Neuropathol. Exp. Neurol.71312–329. 10.1097/NEN.0b013e31824d9882

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12

  DrexelM.PreidtA. P.SperkG. (2012b). Sequel of spontaneous seizures after kainic acid-induced status epilepticus and associated neuropathological changes in the subiculum and entorhinal cortex.Neuropharmacology63806–817. 10.1016/j.neuropharm.2012.06.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  DrexelM.PreidtA. P.KirchmairE.SperkG. (2011). Parvalbumin interneurons and calretinin fibers arising from the thalamic nucleus reuniens degenerate in the subiculum after kainic acid-induced seizures.Neuroscience189316–329. 10.1016/j.neuroscience.2011.05.021

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14

  FabóD.MaglóczkyZ.WittnerL.PékA.ErossL.CzirjákS.et al (2008). Properties of in vivo interictal spike generation in the human subiculum.Brain131485–499. 10.1093/brain/awm297

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  FiskeM. P.AnstötzM.WeltyL. J.MaccaferriG. (2020). The intrinsic cell type-specific excitatory connectivity of the developing mouse subiculum is sufficient to generate synchronous epileptiform activity.J. Physiol.5981965–1985. 10.1113/JP279561

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  FranklinK. B. J.PaxinosG. (2008). The Mouse Brain in Stereotaxic Coordinates, 3rd Edn. San Diego: Academic Press.

  • Google Scholar

• 17

  FroriepU. P.KumarA.Cosandier-RiméléD.HäusslerU.KiliasA.HaasC. A.et al (2012). Altered θ coupling between medial entorhinal cortex and dentate gyrus in temporal lobe epilepsy.Epilepsia531937–1947. 10.1111/j.1528-1167.2012.03662.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  FujiseN.HunzikerW.HeizmannC. W.KosakaT. (1995). Distribution of the calcium binding proteins, calbindin D-28K and parvalbumin, in the subicular complex of the adult mouse.Neurosci. Res.2289–107. 10.1016/0168-0102(95)00882-t

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  GreeneJ. R.TotterdellS. (1997). Morphology and distribution of electrophysiologically defined classes of pyramidal and nonpyramidal neurons in rat ventral subiculum in vitro.J. Comp. Neurol.380395–408. 10.1002/(sici)1096-9861(19970414)380:3<395::aid-cne8>3.0.co;2-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  HaasC. A.DudeckO.KirschM.HuszkaC.KannG.PollakS.et al (2002). Role for reelin in the development of granule cell dispersion in temporal lobe epilepsy.J. Neurosci.225797–5802.

  • Google Scholar

• 21

  HäusslerU.BielefeldL.FroriepU. P.WolfartJ.HaasC. A. (2012). Septotemporal position in the hippocampal formation determines epileptic and neurogenic activity in temporal lobe epilepsy.Cereb. Cortex2226–36. 10.1093/cercor/bhr054

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  HäusslerU.RinasK.KiliasA.EgertU.HaasC. A. (2016). Mossy fiber sprouting and pyramidal cell dispersion in the hippocampal CA2 region in a mouse model of temporal lobe epilepsy.Hippocampus26577–588. 10.1002/hipo.22543

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  HeD. F.MaD. L.TangY. C.EngelJ.BraginA.TangF. R. (2010). Morpho-physiologic characteristics of dorsal subicular network in mice after pilocarpine-induced status epilepticus.Brain Pathol.2080–95. 10.1111/j.1750-3639.2009.00243.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  HeinrichC.NittaN.FlubacherA.MüllerM.FahrnerA.KirschM.et al (2006). Reelin deficiency and displacement of mature neurons, but not neurogenesis, underlie the formation of granule cell dispersion in the epileptic hippocampus.J. Neurosci.264701–4713. 10.1523/JNEUROSCI.5516-05.2006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  HouserC. R. (1990). Granule cell dispersion in the dentate gyrus of humans with temporal lobe epilepsy.Brain Res.535195–204.

  • Google Scholar

• 26

  IshiharaY.FukudaT. (2016). Immunohistochemical investigation of the internal structure of the mouse subiculum.Neuroscience337242–266.

  • Google Scholar

• 27

  JanzP.SavanthrapadianS.HäusslerU.KiliasA.NestelS.KretzO.et al (2017a). Synaptic remodeling of entorhinal input contributes to an aberrant hippocampal network in temporal lobe epilepsy.Cereb. Cortex272348–2364. 10.1093/cercor/bhw093

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  JanzP.SchwaderlappN.HeiningK.HäusslerU.KorvinkJ. G.von ElverfeldtD.et al (2017b). Early tissue damage and microstructural reorganization predict disease severity in experimental epilepsy.elife6:e25742. 10.7554/eLife.25742

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  KiliasA.HäusslerU.HeiningK.FroriepU. P.HaasC. A.EgertU. (2018). Theta frequency decreases throughout the hippocampal formation in a focal epilepsy model.Hippocampus28375–391. 10.1002/hipo.22838

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  KiliasA.TulkeS.BarheierN.RutherP.HäusslerU. (2023). Integration of the CA2 region in the hippocampal network during epileptogenesis.Hippocampus33223–240. 10.1002/hipo.23479

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  KlapsteinG. J.ColmersW. F. (1997). Neuropeptide Y suppresses epileptiform activity in rat hippocampus in vitro.J. Neurophysiol.781651–1661. 10.1152/jn.1997.78.3.1651

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  KnoppA.FrahmC.FidzinskiP.WitteO. W.BehrJ. (2008). Loss of GABAergic neurons in the subiculum and its functional implications in temporal lobe epilepsy.Brain1311516–1527. 10.1093/brain/awn095

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  KnoppA.KiviA.WoznyC.HeinemannU.BehrJ. (2005). Cellular and network properties of the subiculum in the pilocarpine model of temporal lobe epilepsy.J. Comp. Neurol.483476–488. 10.1002/cne.20460

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34

  KulikA.VidaI.LujánR.HaasC. A.López-BenditoG.ShigemotoR.et al (2003). Subcellular localization of metabotropic GABA(B) receptor subunits GABA(B1a/b) and GABA(B2) in the rat hippocampus.J. Neurosci.2311026–11035. 10.1523/JNEUROSCI.23-35-11026.2003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  LedriM.SørensenA. T.MadsenM. G.ChristiansenS. H.LedriL. N.CifraA.et al (2015). Differential effect of neuropeptides on excitatory synaptic transmission in human epileptic hippocampus.J. Neurosci.359622–9631.

  • Google Scholar

• 36

  LisgarasC. P.ScharfmanH. E. (2022). Robust chronic convulsive seizures, high frequency oscillations, and human seizure onset patterns in an intrahippocampal kainic acid model in mice.Neurobiol. Dis.166:105637. 10.1016/j.nbd.2022.105637

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  MaglóczkyZ.WittnerL.BorhegyiZ.HalászP.VajdaJ.CzirjákS.et al (2000). Changes in the distribution and connectivity of interneurons in the epileptic human dentate gyrus.Neuroscience967–25. 10.1016/s0306-4522(99)00474-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  MarxM.HaasC. A.HäusslerU. (2013). Differential vulnerability of interneurons in the epileptic hippocampus.Front. Cell Neurosci.7:167. 10.3389/fncel.2013.00167

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  PaschenE.ElguetaC.HeiningK.VieiraD. M.KleisP.OrcinhaC.et al (2020). Hippocampal low-frequency stimulation prevents seizure generation in a mouse model of mesial temporal lobe epilepsy. Elife9:e54518. 10.7554/eLife.54518

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  PatryloP. R.van den PolA. N.SpencerD. D.WilliamsonA. (1999). NPY inhibits glutamatergic excitation in the epileptic human dentate gyrus.J. Neurophysiol.82478–483. 10.1152/jn.1999.82.1.478

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  SloviterR. S. (1994). On the relationship between neuropathology and pathophysiology in the epileptic hippocampus of humans and experimental animals.Hippocampus4250–253.

  • Google Scholar

• 42

  StewardO. (1976). Topographic organization of the projections from the entorhinal area to the hippocampal formation of the rat.J. Comp. Neurol.167285–314. 10.1002/cne.901670303

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  SutulaT.CascinoG.CavazosJ.ParadaI.RamirezL. (1989). Mossy fiber synaptic reorganization in the epileptic human temporal lobe.Ann. Neurol.26321–330.

  • Google Scholar

• 44

  SuzukiF.JunierM. P.GuilhemD.SørensenJ. C.OntenienteB. (1995). Morphogenetic effect of kainate on adult hippocampal neurons associated with a prolonged expression of brain-derived neurotrophic factor.Neuroscience64665–674. 10.1016/0306-4522(94)00463-f

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  TaubeJ. S. (1993). Electrophysiological properties of neurons in the rat subiculum in vitro.Exp. Brain Res.96304–318.

  • Google Scholar

• 46

  ThomM.SisodiyaS. M.BeckettA.MartinianL.LinW.-R.HarknessW.et al (2002). Cytoarchitectural abnormalities in hippocampal sclerosis.J. Neuropathol. Exp. Neurol.61510–519.

  • Google Scholar

• 47

  TóthE.BokodiV.SomogyváriZ.MaglóczkyZ.WittnerL.UlbertI.et al (2021). Laminar distribution of electrically evoked hippocampal short latency ripple activity highlights the importance of the subiculum in vivo in human epilepsy, an intraoperative study.Epilepsy Res.169:106509. 10.1016/j.eplepsyres.2020.106509

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  TulkeS.HaasC. A.HäusslerU. (2019). Expression of brain-derived neurotrophic factor and structural plasticity in the dentate gyrus and CA2 region correlate with epileptiform activity.Epilepsia601234–1247. 10.1111/epi.15540

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  WoldbyeD. P.MadsenT. M.LarsenP. J.MikkelsenJ. D.BolwigT. G. (1996). Neuropeptide Y inhibits hippocampal seizures and wet dog shakes.Brain Res.737162–168.

  • Google Scholar

• 50

  WoznyC.KnoppA.LehmannT.-N.HeinemannU.BehrJ. (2005). The subiculum: a potential site of ictogenesis in human temporal lobe epilepsy.Epilepsia46(Suppl. 5) 17–21.

  • Google Scholar