Abstract
Dravet syndrome (DS) is an epileptic encephalopathy that still lacks biomarkers for epileptogenesis and its treatment. Dysfunction of NaV1.1 sodium channels, which are chiefly expressed in inhibitory interneurons, explains the epileptic phenotype. Understanding the network effects of these cellular deficits may help predict epileptogenesis. Here, we studied θ-γ coupling as a potential marker for altered inhibitory functioning and epileptogenesis in a DS mouse model. We found that cortical θ-γ coupling was reduced in both male and female juvenile DS mice and persisted only if spontaneous seizures occurred. θ-γ Coupling was partly restored by cannabidiol (CBD). Locally disrupting NaV1.1 expression in the hippocampus or cortex yielded early attenuation of θ-γ coupling, which in the hippocampus associated with fast ripples, and which was replicated in a computational model when voltage-gated sodium currents were impaired in basket cells (BCs). Our results indicate attenuated θ-γ coupling as a promising early indicator of inhibitory dysfunction and seizure risk in DS.
Significance Statement
Biomarkers for Dravet syndrome (DS), a severe type of childhood epilepsy, are currently not available. Despite advances in our understanding of the cellular pathophysiology of DS, the effects on brain network dynamics largely remain to be uncovered. In the present work we show that decreased θ-γ coupling precedes and associates with seizure activity in a DS mouse model, which could be replicated by brain region-specific ablation of NaV1.1, a sodium channel affected in the majority of DS patients. Computational modeling of sodium channel dysfunction in inhibitory interneurons yielded a similar decrease in θ-γ coupling. Together, these results suggest that impaired θ-γ coupling is a promising early indicator of inhibitory dysfunction and seizure risk in DS.
Introduction
Dravet syndrome (DS) is an epileptic encephalopathy characterized by early-onset seizures, followed by severe cognitive and behavioral deficits (Wolff et al., 2006; Dravet, 2011). The majority of DS cases result from a heterozygous loss-of-function mutation in the SCN1A gene, which encodes the pore-forming α1 subunit of voltage-gated sodium channel type 1 (NaV1.1). In mice with a heterozygous knock-out of Scn1a, core clinical features of DS including spontaneous seizures (Yu et al., 2006), cognitive deficits and autism-related behavior (Han et al., 2012) are replicated. Voltage-dependent sodium currents are reduced in GABAergic inhibitory interneurons, while unaffected in hippocampal pyramidal (EX) neurons of DS mice (Yu et al., 2006). In the hippocampus and cortex, NaV1.1 is mostly expressed in interneurons (Ogiwara et al., 2007) and loss of NaV1.1 function in inhibitory populations is sufficient to reproduce the epileptic phenotype of DS mice (Cheah et al., 2012). These studies thus implicate loss of inhibitory functioning in DS pathophysiology.
How loss of NaV1.1 function affects network dynamics, however, remains unclear, while such knowledge is crucial to understanding epileptogenesis and cognitive deficits. Although excitability of cortical inhibitory interneurons was decreased in brain slices of DS mice, in vivo cortical local field potential (LFP) and spontaneous interneuron firing appeared unaffected (De Stasi et al., 2016). Fast-spiking parvalbumin-positive interneurons, a class of inhibitory interneurons that is affected in DS mice (Tai et al., 2014; De Stasi et al., 2016), importantly contribute to LFP γ oscillations (Freund and Katona, 2007; Buzsáki and Wang, 2012). Parvalbuminergic neurons show a particularly high propensity to fire phase-locked to θ oscillations in the hippocampus (Csicsvari et al., 1999) and cortex (Sirota et al., 2008), suggesting a key role of these interneurons in θ-γ cross-frequency coupling (Buzsáki and Wang, 2012). Coupling of γ oscillations to θ rhythm is indeed decreased when synaptic inhibition onto parvalbumin-positive neurons is impaired (Wulff et al., 2009). Recently, local loss of hippocampal NaV1.1 function was shown to impair inhibitory firing and coupling of γ amplitude to θ phase in rats (Sakkaki et al., 2020). How this relates to seizure development, hippocampal and cortical dynamics in DS mice, however, remains to be explored.
Here, we studied spontaneous cortical network dynamics during development in DS mice for early changes in θ or γ oscillations and their cross-frequency phase-amplitude coupling. By inducing local heterozygous or homozygous knock-out of Scn1a in the hippocampus or cortex, we addressed the effect of partial or total ablation of NaV1.1, respectively, within these structures on network dynamics in relation to seizure development. We found that θ-γ coupling was decreased in DS mice compared with wild-type (WT), which, also upon local ablation of NaV1.1, preceded seizure activity. Decreased θ-γ coupling persisted only in animals with spontaneous seizures, and associated with hippocampal fast ripples. In a computational hippocampal network model, loss of voltage-gated sodium currents in inhibitory neurons reproduced the main experimental findings, including a pronounced impairment of θ-γ coupling. These results suggest that θ-γ coupling may serve as an early indicator of inhibitory dysfunction and seizure risk in DS.
Materials and Methods
Animals
Mice with a deletion of exon 8 of the Scn1a gene were created as described previously (Jansen et al., 2020). These global Scn1a-/+ mice (referred to as “DS mice”) were backcrossed to C57BL/6J mice for at least five generations. Both male and female DS mice and WT littermates were used for experiments. In the conditional Scn1a mouse line, exon 8 of the Scn1a gene is flanked by LoxP sites (Jansen et al., 2020), which allows for region-specific ablation after exposure to Cre recombinase. Heterozygous and homozygous floxed (Scn1afl/+ or Scn1afl/fl, respectively) mice and WT littermates were used for the viral injection experiments. Animals had free access to food and water and were kept under standard housing conditions (temperature of 22 ± 1.5°C, 12/12 h light/dark cycle). All experiments were approved by local and national ethical committees following recommendations of the European Communities Council Directive (2010/63/EU) and conducted in accordance with ARRIVE guidelines.
Surgery and viral infection
For chronic recordings in DS mice and WT littermates, LFP electrodes (75-μm platinum/iridium; PT6718, Advent Research Materials) were implanted at P20–P21 under isoflurane anesthesia (induction 4%; maintenance 1.5%). Electrodes were implanted in primary visual cortex [V1; −3.5 mm anteroposterior (AP), ±2.4 mm mediolateral (ML), and −0.5 mm dorsoventral (DV) relative to bregma] and primary motor cortex (M1; +1.5 mm AP, ±1.8 mm ML, −0.5 mm DV). A subset of WT and DS mice was implanted with an additional LFP electrode in the right dorsal hippocampus (−2.0 mm AP, +2.2 mm ML, −1.4 mm DV).
To induce local knock-out of Scn1a, an AAV vector expressing mCherry-Cre under the EF1α-promotor (viral prep # 55632, a gift from Karl Deisseroth, Addgene) was injected in postnatal day (P)42–P56 Scn1afl/+ and Scn1afl/fl mice. Injections (500 nl per injection, 50 nl/min) were made bilaterally in dorsal (−2.0 mm AP, ±2.2 mm ML, −1.4 mm DV) and ventral (−3.0 mm AP, ±3.1 mm ML, −2.6 mm DV) hippocampus, or in V1, followed by implantation of LFP electrodes. For cortically injected animals, LFP electrodes were implanted in V1 and M1. The use of heterozygous and homozygous floxed mice was predicted to result in a partial and total ablation of NaV1.1 following local expression of Cre, respectively. WT littermates that underwent the same procedure served as controls. For cortical injections, in a subset of Scn1afl/+ mice, a different AAV vector was used (expressing eGFP-Cre; viral prep #105545, a gift from James M. Wilson, Addgene) to allow comparison with previously published data (Jansen et al., 2020).
Pharmacology
Cannabidiol (CBD; Tocris) was dissolved in DMSO, which was diluted in a vehicle solution (5% DMSO: 5% cremophor: 90% saline) immediately before injection. P35–P36 DS mice received an intraperitoneal injection of CBD (100 mg/kg body weight) or vehicle on two consecutive days between 8 and 10 A.M., in a randomized fashion and by an experimenter blinded to the treatment.
Immunohistochemistry
After experiments, animals were euthanized by CO2 and transcardially perfused with PBS and 4% paraformaldehyde. Following removal and postfixation of brain tissue, coronal sections (20 μm) were prepared. Sections were heated in 10 mmol/l sodium citrate buffer with 0.05% Tween for 10 min at 80°C. Subsequently, a blocking solution containing 10% normal goat serum was applied, and sections were incubated in rabbit anti-NaV1.1 (1:200; Alomone Labs), followed by incubation in goat anti-rabbit Cy2 (1:200; Jackson ImmunoResearch). Mounting was performed in glycerol/PBS (1:1) containing 12.5 mg/ml sodium azide and 1 μl/ml Hoechst-33258. Sections were examined using an epifluorescence or confocal microscope with appropriate filter sets.
Data acquisition and analysis
DS mice and WT littermates were videotaped from P18 to P19 in a PhenoTyper home cage (Noldus) to allow detection of early spontaneous seizures, until the day of surgery at P20–P21. Immediately after surgery, both these mice and virally injected mice were connected to a seven-channel commutator in a Faraday cage for continuous recordings of LFP and video. Signals were preamplified (3×), filtered (0.05–500 Hz), and amplified (200× and 400× for hippocampal and cortical recordings, respectively) using custom-build hardware, and digitized at 5 kHz (Power 1401 and Spike2 software, CED). All LFP recordings were post hoc inspected for epileptiform activity. The Racine scale (Racine, 1972) was used for scoring of behavior from video recordings during discharges that lasted >5 s. Seizures that progressed to stage 3–5 were included in analyses, as these seizures could also be reliably detected when only video recordings were available (i.e., for P18–P21 mice). For detection of interictal spikes (IISs) following local knock-out of Scn1a, local LFP power in the 30- to 100-Hz range was calculated every 24 h for 1 h of non-REM sleep, a vigilance state during which IIS are prevalent in DS mice (Kalume et al., 2015). Candidate events exceeding a threshold of three SDs of signal were visually inspected, and animals were considered showing IIS if more than five IIS were observed during the analyzed 1-h window on that day. To detect fast ripples, hippocampal LFP recordings were bandpass filtered (infinite impulse response Butterworth filter of second order) between 200 and 600 Hz and candidate events were detected by thresholding (>10 SDs), which effectively distinguished fast ripples from IIS. Only events that occurred during behavioral immobility, assessed from the reference signal and locomotor activity recorded by an infrared motion detection sensor in the Faraday cage, were visually inspected, as muscle artifacts were often found in the 200- to 600-Hz frequency range during active wakefulness.
To assess cortical network dynamics during REM sleep and active wakefulness, LFP recorded during the morning (6–12 A.M.) was analyzed, starting at P23, allowing 2–3 d of habituation to the recording setup. Vigilance state was determined per 5-s epoch using V1 LFP, the reference signal and locomotor activity recorded by the infrared motion detection sensor. Epochs during seizures and/or postictal periods (1 h after the end of a seizure) were excluded from further analyses. REM sleep was defined by a θ (5–10 Hz) to δ (1–4 Hz) power ratio of >2.5, in the absence of locomotor activity or motion artifacts. Active wakefulness was defined by a θ-δ power ratio of >2 during epochs containing locomotor activity and/or high variance in the reference signal. Epochs were artifact-rejected, low-pass-filtered (Chebyshev IIR eighth-order filter) and down-sampled to 1 kHz. A total of 10 min of randomly selected epochs was used for further analyses. Power spectra were computed by applying a Hamming window over each epoch, followed by a fast Fourier transform and averaging of the resulting power spectra. Phase-amplitude coupling was assessed using a measure referred to as the modulation index (MI), calculated as described previously (Tort et al., 2010). In short, LFP epochs were bandpass-filtered for phase (3–14 Hz, 2-Hz bandwidths, steps of 1 Hz) and amplitude (40–200 or 40–300 Hz, 4-Hz bandwidths, steps of 2 Hz) frequencies. The Hilbert transform was applied to extract phase and amplitude information. Each phase frequency was binned in 18 bins of 20°, and mean amplitude per amplitude frequency was calculated inside each bin. A distribution was subsequently obtained by normalization of the mean amplitude in each phase bin by the sum over all bins. To compare the resulting phase-amplitude distribution with a uniform distribution, the Kullback–Leibler distance was calculated. The MI was obtained by normalizing this metric to the logarithm of the number of phase bins. For group analyses, θ-γ coupling was expressed in a single value by averaging over all θ and γ frequencies, or over specific γ bands, as indicated. All analyses were performed using custom-written MATLAB (versions used: 2018b and 2019b, MathWorks) scripts.
Computer modeling
The basic equations for the membrane potential of individual neurons and various ion channels were adopted from (Kopell et al., 2010). The model network consists of five excitatory (EX) neurons, five oriens-lacunosum moleculare (OLM) interneurons, and five fast-spiking basket cell (BC) interneurons, with all-to-all connections. We noticed that increasing the network size did not change the conclusions from the model as long as the synaptic connections were properly scaled (data not shown). The equations for individual cells were modified and extended to incorporate the dynamics of various ion species in the intracellular and extracellular spaces of the neurons and neurotransmitter homeostasis using a formalism previously developed (Cressman et al., 2009; Hübel et al., 2017; Ullah, 2019). Full details about the model are given in Extended Data Fig. 10-1. Here, we describe the main equations.
The membrane potential, Vm, of each EX, BC, and OLM neuron in the network is controlled by various Na+ (INa), K+ (IK), and Cl– (ICl) currents, Na+ and K+ currents because of the Na+/K+-ATPase (Ipump), and random inputs from neurons that were not a part of the network, modeled as a Gaussian distribution (Istoch), that is,
where Cm is the membrane capacitance and superscripts EX, BC, and OLM correspond to EX, BC, and OLM neurons, respectively. Various currents used for the different neuron types are given in the section Extended Data Fig. 10-1, membrane potential.
To model the dynamics of various ionic species, we considered the extracellular space (ECS) as a separate compartment surrounding each cell. Each neuron exchanged ions with its ECS compartment through active and passive currents, and the Na+/K+-ATPase. The ECS compartment can also exchange K+ with the glial compartment, and perfusion solution (or vasculature in intact brain; Cressman et al., 2009; Ullah and Schiff, 2010; Wei et al., 2014; Ullah et al., 2015). Including diffusion of ions between the ECS compartments of neighboring neurons did not affect the results and was therefore ignored.
The change in extracellular K+ concentration ([K+]o) is a function of IK, Ipump, uptake by glia surrounding the neuron (Iglia), and the diffusion between the neuron and bath perfusate (Idiff).
Where β is the ratio of intracellular space (ICS) to ECS. 
is the conversion factor from current units to flux units. F, 
, and ωi are the Faraday's constant, surface area of the neuron, and intracellular volume of the neuron, respectively. The factor 2 in front of Ipump is because of the fact that the Na+/K+ pump extrudes two K+ in exchange for three Na+ ions.
The rate of change of intracellular Na+ concentration ([Na+]i) is controlled by INa and Ipump (Cressman et al., 2009), that is,
The fluxes used in Equations 2, 3 are listed in Extended Data Fig. 10-1, ion concentrations. The concentrations of intracellular K+ ([K+]i), intracellular Cl– ([Cl–]i), extracellular Na+ ([Na+]o), and extracellular Cl– ([Cl–]o) are given by conservation laws and are also listed in the same section.
In addition to ion concentration dynamics, our model also incorporates the details of neurotransmitter homeostasis. Glutamate and GABA in the synaptic cleft and ECS are regulated by (1) spiking-related release from EX (for glutamate), BC and OLM (for GABA), (2) uptake by presynaptic and postsynaptic neurons and glia through glutamate and GABA transporters, (3) diffusion between the cleft and ECS, and (4) the recycling into releasable pools. The model also includes the binding and unbinding of glutamate to NMDA and AMPA receptors, and GABA to GABAA receptors, and the currents because of these receptors. Details of the equations are given in Extended Data Fig. 10-1, neurotransmission.
The LFP is modeled by considering a non-spiking EX neuron in the network as described previously (Kopell et al., 2010). This neuron receives the same synaptic inputs as the active EX neurons in the network, but it does not generate action potentials because of the absence of the main currents responsible for action potential generation. The model LFP thus reflects subthreshold voltage changes in the EX-cell population. Similar to the in vivo data analysis, 10 min of simulated LFP data were used to obtain power spectra and phase-amplitude coupling metrics, as described earlier.
Experimental design and statistical analysis
Data from 89 mice were included, including 15 WT and 26 DS mice for video-EEG recordings, nine DS mice for pharmacological experiments, and 13 WT, 14 Scn1afl/+, and 12 Scn1afl/fl mice for AAV vector injections. Data are presented as mean ± SEM, unless indicated otherwise. Statistical testing was conducted in MATLAB or GraphPad Prism (GraphPad Software). Based on whether data were paired and their distribution normality, which was assessed by the Anderson–Darling test, data were compared using the Welch's t test, paired t test, Mann–Whitney test, or Wilcoxon test for single comparisons. For multiple comparisons, a one or two-way ANOVA with Tukey's test, repeated-measures ANOVA with Dunnett's test or Kruskal–Wallis with Dunn's test was used. Sample sizes were not based on statistical methods, but instead informed by those frequently used in the field; p < 0.05 was considered to indicate significance.
Results
Modulation of cortical γ amplitude by θ phase is attenuated in juvenile DS mice
Cortical LFP and video were recorded in freely behaving Scn1a-/+ (DS) and WT mice starting from P20–P21. Modulation of cortical γ oscillations by θ rhythm is low during vigilance states lacking prominent θ oscillations, i.e., during quiet wakefulness and non-REM sleep (Scheffzük et al., 2011). We therefore limited our analyses to epochs during vigilance states with prominent cortical θ oscillations, i.e., REM sleep (examples in Fig. 1A) and active wakefulness. Average V1 LFP power spectral density for θ (5–10 Hz), low γ (40–90 Hz), and high γ (90–160 Hz) frequencies during REM sleep at P23 was similar between genotypes (Fig. 1B,C), although a decrease in θ peak frequency was observed in DS mice (Fig. 1D). Comodulogram analyses of V1 LFP in WT mice showed that θ phase modulated the amplitude of low and high γ (Fig. 1E). Although this pattern was also present in DS mice, the modulation of γ power by θ phase was significantly impaired for both frequency ranges (Fig. 1F). The θ phase at which γ amplitude peaked was similar between genotypes, with high γ peaking significantly later in the θ cycle than low γ for both WT and DS mice (Fig. 1G,H). Average REM sleep duration was similar between WT and DS mice (48.9 ± 4.2 and 46.6 ± 3.0 s, respectively; t(26) = 0.45, p = 0.66, t test), thereby excluding decreased REM sleep duration as a potential cause for impaired θ-γ coupling strength (Bandarabadi et al., 2019). Similar results were obtained from epochs during active wakefulness (Fig. 2), although in both WT and DS mice, modulation of high γ oscillations was largely absent, and modulation of both low and high γ frequencies was lower when compared with REM sleep (Fig. 2F–I). These data indicate that cortical θ-γ coupling is attenuated in juvenile DS mice, which cannot be explained by changes in θ or γ power.
DS mice display decreased cortical θ-γ cross-frequency coupling at P23. A, Example traces of raw, θ-filtered (5–10 Hz) and γ-filtered (40–160 Hz) LFP in the V1 during REM sleep in a WT and DS mouse. B, Average V1 LFP power spectra during REM sleep in WT (black) and DS (red) mice. C, No significant differences in power were present between genotypes within θ, low (40–90 Hz) or high γ (90–160 Hz) frequency ranges. Peak θ frequency was significantly reduced in DS mice (D; detail of mean power in θ range, black and red dots indicate average peak θ frequency; t(27) = 2.5, *p = 0.020, Welch's t test). E, Average phase-amplitude comodulograms of V1 during REM sleep in WT (n = 9) and DS (n = 19) mice at P23. F, θ-γ Coupling was significantly reduced in DS mice, for both low and high γ frequencies (t(8.9) = 3.8 and t(9.2) = 2.9, respectively, *p = 0.019 and **p = 0.004, Welch's t test). G, Normalized γ amplitude per 20° phase bin of the θ cycle for V1 during REM sleep. No significant differences were present between WT and DS mice. H, High γ amplitude peaked significantly later in the θ cycle than low γ in both genotypes (*p < 0.005, Wilcoxon test).
Deficient cortical θ-γ coupling during active wakefulness in DS mice. A, Average V1 LFP power spectra during wakefulness in WT (black) and DS (red) mice. No significant differences were present between genotypes in power within θ (5–10 Hz), low γ (40–90 Hz), or high γ (90–160 Hz) frequency ranges (B). C, Peak θ frequency was significantly reduced in DS mice (t(22.7) = 3.5, *p = 0.002, Welch's t test). D, Average phase-amplitude comodulograms of V1 LFP during active wakefulness in WT (n = 9) and DS (n = 19) mice at P23. E, θ-Low γ coupling was significantly reduced in DS mice (t(7.9) = 2.7, *p = 0.030, Welch's t test). For both WT (F, G) and DS (H, I) mice, θ modulation of γ amplitude was largely limited to low γ frequencies during active wakefulness, while modulation of both low and high γ frequencies was lower when compared with REM sleep.
Impaired θ-γ coupling in the hippocampus of juvenile DS mice
Previous studies indicate that θ oscillations recorded in the neocortex are volume-conducted currents originating from the hippocampus (Gerbrandt et al., 1978; Sirota et al., 2008). The hippocampus, a brain region that importantly contributes to the epileptic phenotype of DS mice (Liauard et al., 2013; Stein et al., 2019; Jansen et al., 2020), also displays θ-γ phase-amplitude coupling (Bragin et al., 1995; Buzsáki et al., 2003). To establish whether hippocampal θ-γ coupling is disturbed in DS mice, we assessed hippocampal LFP in WT (n = 6) and DS (n = 7) mice during REM sleep at P23. Similar to results obtained from V1, θ and γ power were not significantly different between genotypes (Fig. 3A,B), whereas θ peak frequency was reduced in DS mice (Fig. 3A, inset). Variability in γ power and θ-γ coupling appeared increased when compared with cortical LFP, which may be explained by slight variations in electrode position within the hippocampus since patterns of θ-γ coupling are different among hippocampal sublayers (Schomburg et al., 2014). Hippocampal θ-γ coupling was, however, robustly decreased in DS mice (Fig. 3C,D).
Deficient θ-γ coupling in the hippocampus of DS mice. A, Power spectral density of LFP recorded in the dorsal hippocampus during REM sleep in WT (black) and DS (red) mice. Peak θ frequency was significantly reduced in DS mice (inset; t(6.9) = 3.2, *p = 0.011, Welch's t test). No significant differences were present between genotypes in power within θ (5–10 Hz) and γ (40–300 Hz) frequency ranges (B). C, Average phase-amplitude comodulograms of the same LFP signals. D, θ-γ Coupling was significantly reduced in DS mice (t(8.1) = 3.64, *p = 0.006, Welch's t test).
Cortical θ-γ coupling normalizes in DS mice without spontaneous seizures
The profound cortical θ-γ coupling (for both low and high γ) observed during REM sleep led us to limit subsequent analyses to this vigilance state. First, we assessed whether developmental changes in cortical θ-γ coupling occurred over the three-week recording period by analyzing V1 LFP during REM sleep at 2, 7, 14, and 21 d after surgery (i.e., P23, P28, P35, and P42; Fig. 4). θ-γ Coupling increased over this period in both WT and DS mice (Fig. 4A,B). Similar to the observations at P23, DS mice showed decreased θ-γ coupling when compared with WT mice at P42 [MI values of 5.2 ± 0.7 × 10−4 and 10.4 ± 1.5 × 10−4, respectively; calculated for the total γ range (40–160 Hz); t(11) = 3.2, p = 0.009, Welch's t test].
Persistently reduced cortical θ-γ coupling in DS mice that show spontaneous seizures. A, Average phase-amplitude comodulograms of V1 showing progression of cortical θ-γ coupling during postnatal weeks 4–7. B, Average MI, calculated for θ (5–10 Hz) and total γ (40–160 Hz), was markedly increased in both WT and DS mice at P42, when compared with P23 (t(8) = 7.0 and t(15) = 6.4, respectively, *p < 0.0001, paired t test). C, Total frequency of spontaneous seizures in DS mice (n = 19). A subset of DS mice (n = 3) died during the recording period, between P24–P26. D, Progression of cortical θ-γ cross-frequency coupling for WT (black) and DS mice with (gray) and without (red) spontaneous seizures over the recording period. Note that at P23 the MI was significantly decreased compared with WT for both groups of DS mice, while at later ages this decrease was maintained only in DS mice with seizures for both low (40–90 Hz) and high (90–160 Hz) γ (F(2,22) = 9.5, *p < 0.05, **p < 0.005, repeated-measures ANOVA with Tukey's test).
The recording period included postnatal week 4, during which spontaneous seizures are most prevalent in DS mice (Kalume et al., 2013). Accordingly, here, spontaneous seizures in DS mice were most often detected between P21–P28 (Fig. 4C). Three DS mice had a fatal seizure (during P24–P26), whereas a subset of DS mice (seven out of 19) remained seizure-free over the three-week recording period. As early seizures, induced by hyperthermia or a chemoconvulsant, were shown to aggravate the epileptic phenotype and cognitive deficits in a knock-in DS mouse model (Salgueiro-Pereira et al., 2019), we assessed whether θ-γ coupling progressed differently in DS mice that did or did not express seizures during the recording period. Although seizure-free DS mice did show a decrease in cortical θ-γ coupling for both low and high γ at baseline (P23), no significant differences were present at P28, P35, and P42 in these mice compared with WT littermates (Fig. 4D). In contrast, θ-γ coupling was persistently decreased, for both low and high γ, in DS mice that showed at least one seizure during the three-week recording period (seizure count per mouse: 5.1 ± 1.5, range 1–19), when compared with WT mice (Fig. 4D).
Decreased progression of cortical θ-γ coupling in DS mice is improved by CBD
To test whether antiseizure drug treatment could improve attenuated cortical θ-γ coupling, a separate group of DS mice (n = 9) was used to assess the effects of CBD. CBD decreased seizure frequency in patients with DS (Devinsky et al., 2017), while in DS mice it decreased seizure frequency and severity and increased inhibitory neurotransmission in the hippocampus (Kaplan et al., 2017). Here, in DS mice, the progression of cortical θ-γ coupling ceased after P35, contrasting the steady progression of θ-γ coupling observed in WT mice (Fig. 5A,B). DS mice received a single intraperitoneal injection of CBD (100 mg/kg) or vehicle at P35 using a randomized crossover design at an interval of 24 h. θ-γ Coupling was assessed in V1 LFP during REM sleep up to 5 h after intraperitoneal injection, during which plasma and brain CBD levels are high in mice (Deiana et al., 2012). When compared with vehicle, acute treatment of DS mice with CBD resulted in a modest but consistent increase in θ-γ coupling (Fig. 5C,D).
Acute treatment with CBD improves progression of θ-γ coupling in DS mice. WT mice showed steady progression of average θ-γ cross-frequency coupling over time [A; calculated for total γ (40–160 Hz); F(3,8) = 15.5], while for DS mice progression stagnated after P35 (B; F(3,15) = 27.5, *p < 0.05, **p < 0.005, repeated-measures ANOVA with Tukey's test). Note that analyses presented in A, B are performed on the same groups as included in Figure 2, excluding DS mice that died during the recording (n = 3). C, Average phase-amplitude comodulograms of V1 following treatment with vehicle (VEH) or CBD (100 mg/kg) in DS mice between P35–P42. D, Average MI, calculated for θ (5–10 Hz) and total γ (40–160 Hz), was significantly increased following treatment with CBD (t(8) = 5.1, *p = 0.0009, ratio paired t test).
Attenuation of cortical θ-γ coupling precedes seizures following hippocampal-specific knock-out of Scn1a
The early attenuation of cortical θ-γ coupling that was observed in DS mice at P23 occurred regardless of the epileptic phenotype. Thus, the deficit in cross-frequency coupling may be directly related to the inhibitory dysfunction caused by impaired sodium currents in the DS model (Yu et al., 2006). Alternatively, attenuated coupling may result from a secondary effect of inhibitory dysfunction on cortical connectivity during critical period plasticity (Takesian and Hensch, 2013; Cardin, 2018). We aimed to address this issue by introducing the Scn1a defect at adult age, using local brain injections of an AAV vector expressing mCherry-Cre in P42–P56 mice in which exon 8 of the Scn1a gene is floxed by LoxP sites (in one or both copies of the gene, i.e., Scn1afl/+ or Scn1afl/fl, respectively). This approach additionally allowed assessment of the effects of brain region-specific NaV1.1 dysfunction on θ-γ coupling.
Hippocampal ablation of NaV1.1 results in an increased susceptibility to seizures (Stein et al., 2019; Jansen et al., 2020). Here, we studied the effect of exposure to Cre recombinase in the dorsal and ventral hippocampus (Fig. 6A) of Scn1afl/fl mice, resulting in complete NaV1.1 ablation, and Scn1afl/+ mice, to study also partial NaV1.1 ablation (Fig. 6B) which may be more relevant in light of the haploinsufficiency observed in patients and Scn1a-/+ DS mice. Within 28 d after injection, all Scn1afl/fl mice and the majority of Scn1afl/+ mice developed IISs (n = 6/6 and n = 6/7, respectively; Fig. 6D) and spontaneous seizures (n = 6/6 and n = 5/7, respectively; Fig. 6C,E), whereas neither were observed in WT controls injected with AAV-mCherry-Cre (n = 7). Statistical comparison indicated a significant increase in seizure count for Scn1afl/fl and Scn1afl/+ mice when compared with controls (18.8 ± 6.8 and 9.4 ± 3.8 seizures, respectively; H(2) = 11.7, p < 0.05, Kruskal–Wallis with Dunn's test), but seizure count did not differ significantly between the two genotypes (p = 0.68). Average onset of IIS was 11.5 d (range 9–14), while the first seizure was observed 13.1 d (range 10–18) after injection. Cortical θ-low γ coupling was impaired already on day 7 after injection for Scn1afl/fl mice (Fig. 6F,G). Also, θ-high γ coupling was impaired at this time in Scn1afl/fl and Scn1afl/+ mice (Fig. 6F,H). Modulation of both γ frequency ranges did not recover over subsequent weeks (Fig. 6G,H). Power changes were limited to an increase in γ power in Scn1afl/fl mice on day 21 after injection (1.5 ± 0.1 × 10−3 and 2.2 ± 0.3 × 10−3 mV2/Hz, on days 2 and 21, respectively; t(6) = 3.1, p = 0.042, paired t test). These data indicate that attenuation of cortical θ-γ coupling precedes IIS and spontaneous seizures following hippocampal deletion of NaV1.1, independent of changes in power spectral density.
Chronic seizure activity following hippocampal ablation of NaV1.1 is preceded by decreases in cortical θ-γ coupling. A, Example of dorsal and ventral hippocampal areas (dHC and vHC, respectively) infected by AAV-mCherry-Cre (red; Hoechst in blue). Scale bar: 500 µm. The targeted area included the CA1 and CA2 region (boxed area expanded in bottom inset; dashed line indicates border CA1/2; scale bar: 200 µm). B, Detail of hippocampal CA1/2 region showing reduced NaV1.1 staining in cells infected by AAV-mCherry-Cre (arrowheads indicate double-labeled cells) in Scn1afl/+ (HET) and Scn1afl/fl (HOM) when compared with WT. Scale bar: 50 µm. C, Spontaneous seizure in a HOM mouse on day 19 after hippocampal AAV-mCherry-Cre. D, Proportion of HET (n = 7) and HOM (n = 6) mice with IISs, after hippocampal AAV-mCherry-Cre on day 0. No IIS were observed in WT (n = 7). E, Average frequency of spontaneous seizures following AAV-mCherry-Cre. F, Average phase-amplitude comodulograms of V1 at different time points after injection. G, H, MI, normalized to values obtained on day 2, was significantly decreased for both low γ (40–90 Hz; G) and high γ (90–160 Hz; H) frequencies in HET and HOM mice, when compared with WT mice (F(2,17) = 14.5 and F(2,17) = 24.9, respectively, *p < 0.05, **p < 0.005; two-way ANOVA with Tukey's test).
Attenuation of θ-high γ coupling occurs independently of seizures following cortical-specific knock-out of Scn1a
In DS mice, both hippocampal and cortical regions may contribute to seizures (Jansen et al., 2020), and inhibitory dysfunction has been demonstrated in both structures (Yu et al., 2006; Tai et al., 2014). We therefore studied the effect of local cortical knock-out of Scn1a by injecting AAV-mCherry-Cre in V1 in Scn1afl/+ and Scn1afl/fl mice, and WT controls (Fig. 7A,B). Interestingly, all Scn1afl/fl mice developed IIS during the recording period, whereas this was the case for only a subset of Scn1afl/+ mice (n = 6/6 and n = 2/7, respectively; Fig. 7D). Spontaneous seizures were observed in a subset of Scn1afl/fl mice (n = 2/6), but never in Scn1afl/+ mice (Fig. 7E). Previously, we reported a higher incidence of seizures in Scn1afl/fl mice (n = 6/6) following cortical injections with an AAV vector expressing GFP-Cre under the CMV promotor (Jansen et al., 2020). As the discrepancy with our current data suggests a greater potency of the latter AAV vector, we tested whether this vector (AAV-GFP-Cre) would result in seizures following cortical injections in 6 Scn1afl/+ mice. Again, none of the Scn1afl/+ mice developed seizures, indicating that homozygous, but not heterozygous, knock-out of Scn1a in V1 predisposes animals to develop spontaneous seizures. Early attenuation of θ-γ coupling, evident 7 d after injection, was restricted to the high γ frequency range and persisted in Scn1afl/fl mice but recovered in Scn1afl/+ mice (Fig. 7F–H). Power changes were limited to an increase in γ power in Scn1afl/+ mice on day 21 after injection (1.8 ± 0.2 × 10−3 and 2.4 ± 0.2 × 10−3 mV2/Hz, on days 2 and 21, respectively; t(6) = 12.1, p < 0.0001, paired t test), while θ peak frequency was decreased in both Scn1afl/+ and Scn1afl/fl on day 21 (difference of −0.2 ± 0.1 and −0.4 ± 0.1 Hz, respectively; t(6) = 2.6, p = 0.038, and t(5) = 3.9, p = 0.012, paired t test). These observations suggest that cortical NaV1.1 is required for modulation of high γ by θ.
Cortical ablation of NaV1.1 specifically decreases θ-high γ coupling. A, Example of the V1 area infected by AAV-mCherry-Cre (red; Hoechst in blue). Scale bar: 500 µm. B, Detail of V1 showing reduced NaV1.1 staining in cells infected by AAV-mCherry-Cre (arrowheads indicate double-labeled cells) in Scn1afl/+ (HET) and Scn1afl/fl (HOM) when compared with WT. Scale bar: 50 µm. C, Spontaneous seizure recorded in an Scn1afl/fl mouse on day 18 after cortical AAV-mCherry-Cre. D, Proportion of HET (n = 7) and HOM (n = 6) mice with IISs, following cortical AAV-mCherry-Cre on day 0. No IIS were observed in WT mice (n = 6). E, Average frequency of spontaneous seizures following injection of AAV-mCherry-Cre. F, Average phase-amplitude comodulograms of V1 at different time points after injection. G, H, MI, normalized to values obtained on day 2, was significantly decreased for high γ (H), but not low γ, in HET and HOM mice when compared with WT mice (F(2,16) = 18.5 and F(2,16) = 1.3, respectively, *p < 0.05, **p < 0.005; two-way ANOVA with Tukey's test).
Hippocampal fast ripples and attenuated θ-γ coupling precede epilepsy onset following hippocampal-specific knock-out of Scn1a
In addition to cortical LFP, mice that were injected in the hippocampus had local LFP recorded from AAV injection sites. Hippocampal θ-γ coupling, θ and γ power, and θ peak frequency decreased over time following injection in both Scn1afl/+ and Scn1afl/fl mice (Fig. 8) and comodulograms showed modulation of different γ frequency ranges, including those above 150 Hz (Fig. 9A). Interestingly, the γ frequency range that showed the strongest modulation by θ phase decreased over time following local Scn1a knock-out (Figs. 8A, 9A). As such, the decrease in θ-γ coupling (Fig. 8B) was caused by a strong reduction in θ-modulated high γ, while low γ was not significantly affected (Fig. 8C,D).
Deficient hippocampal θ-γ coupling and power following local ablation of NaV1.1. A, Average phase-amplitude comodulograms of LFP obtained from dorsal hippocampus during days 2 (left) and 14 (right) after hippocampal injection of mCherry-Cre. B, θ-γ Coupling was reduced in Scn1afl/+ and Scn1afl/fl mice on day 14 (total θ-γ MI normalized to day 2; F(2,17) = 13.9, *p < 0.01, one-way ANOVA with Dunnett's test). C, MI normalized per 2-Hz γ frequency bin revealed that MI was specifically reduced for γ frequencies >150 Hz, while showing large variation for the lower γ range. D, Statistical comparison revealed no difference for θ-modulated low γ (F(2,17) = 1.3, p = 0.31), but a significant reduction for high γ (F(2,17) = 11.4, p = 0.001) for both Scn1afl/+ and Scn1afl/fl mice (*p = 0.004 and *p = 0.001, ANOVA with Dunnett's test). E, Hippocampal LFP power spectral density during REM sleep in WT (n = 7), Scn1afl/+ (HET, n = 7), and Scn1afl/fl (HOM, n = 6) mice during days 2 (black), 14 (red), and 21 (gray) after hippocampal injection of mCherry-Cre. Power was significantly reduced only in HET and HOM mice. Asterisks of corresponding color indicate significant differences of power in θ (5–10 Hz), low γ (40–150 Hz), or high γ (150–300 Hz) on day 2 (*p < 0.05, paired t tests).
θ-γ Coupling is affected more in areas that show fast ripples following hippocampal NaV1.1 ablation. A, Example phase-amplitude comodulograms of LFP in the dorsal and ventral hippocampus (dHC and vHC, respectively) following hippocampal injection of AAV-mCherry-Cre on day 0 in an Scn1afl/fl mouse. B, Examples of fast ripples observed in the ventral hippocampus of an Scn1afl/fl (left) and Scn1afl/+ (right) mouse. Fast ripples (detailed in insets, corresponding to the shaded area) were observed unilaterally (vHC ipsi) and had a high power in the 250- to 500-Hz range, as shown by the spectrogram of the ipsilateral vHC signal at the bottom. Their onset preceded or coincided with IISs on the contralateral side (vHC contra). C, For hippocampal LFP that contained fast ripples (ipsi), θ-γ coupling, assessed by the MI (averaged for the total γ frequency range, 40–300 Hz), was significantly decreased in the 12 h preceding the first seizure (preictal) when compared with day 2 after injection (baseline; *p = 0.031, Wilcoxon test). In parallel, the γ frequency showing peak MI was significantly reduced (D; p = 0.031, Wilcoxon test). MI changes of contralateral LFP that did not contain fast ripples (contra) did not reach statistical significance.
Hippocampal fast ripples (250–500 Hz) appear to occur specifically in the seizure onset zone (Jefferys et al., 2012), independent of cell loss associated with hippocampal sclerosis (Jiruska et al., 2010). To study the association between fast ripples and local θ-γ coupling following hippocampal-specific knock-out of Scn1a, we inspected hippocampal LFP recordings for fast ripples over 12 h preceding the first spontaneous seizure. In 7 out of 13 animals (n = 3 Scn1afl/+ and n = 4 Scn1afl/fl mice), we found fast ripples that were superimposed on interictal discharges (Fig. 9B) in either the ventral (n = 4) or dorsal (n = 3) hippocampus. Following exclusion of one animal with bilateral fast ripples, comparison of hippocampal θ-γ coupling ipsilateral and contralateral to the fast ripples revealed a profound reduction in average MI for the ipsilateral electrode (Fig. 9C). The difference between baseline and preictal MI was also significantly greater for the ipsilateral compared with the contralateral electrode (−7.8 ± 1.5 × 10−4 and −2.1 ± 0.6 × 10−4, respectively; p = 0.002, Mann–Whitney test). Similarly, the γ frequency that showed the strongest modulation by θ phase was decreased in the hours preceding seizure activity, specifically for the ipsilateral electrode (Fig. 9D). The drop in γ frequency between baseline and the preictal period was not significantly different between sides (−79 ± 19 Hz and −36 ± 18 Hz for ipsilateral and contralateral, respectively; p = 0.14, Mann–Whitney test). These data confirm θ-γ coupling as a marker for epileptogenesis in these mice, as its attenuation precedes the first spontaneous seizure in an area corresponding to the seizure onset zone.
Impairment of voltage-gated sodium currents in inhibitory neurons attenuates θ-γ coupling and θ frequency in a computational model
Our in vivo experiments indicate that θ-γ coupling decreases before the first spontaneous seizure, suggesting that epilepsy-related remodeling is not responsible for this deficit in network functioning. To test whether decreased θ-γ coupling may result from loss of voltage-gated sodium currents in inhibitory neurons, we developed a hippocampal network model in which we varied the amount of voltage-gated sodium channel (NaV) function by decreasing their sodium conductance in modified Hodgkin–Huxley neurons representing fast-spiking BC neurons (Wang and Buzsáki, 1996) and OLM neurons (Tort et al., 2007). Synaptic connections between BC, OLM, and EX neurons were all-to-all (Fig. 10A) and five neurons of each type were included in the model. LFP generated by the model was dominated by θ oscillations both at baseline (BC: 10 mS/cm2, OLM: 35 mS/cm2) and after decreasing sodium conductance in inhibitory neurons (BC: 6 mS/cm2, OLM: 28 mS/cm2; Fig. 10B). θ Peak frequency however was shifted to the left in case of impaired NaV function when compared with baseline, while power in the γ range was comparable (9.49 × 10−5 and 9.51 × 10−5, respectively; Fig. 10C). Modulation of γ by θ was present in both conditions, but reduced in case of impaired NaV function, which was most evident for γ oscillations >100 Hz (Fig. 10D). Impairment of inhibitory NaV function in this computational model thus faithfully replicates both the reduced θ frequency and reduced θ-γ coupling observed in DS mice.
Modeling impaired NaV functioning in inhibitory neurons replicates the attenuated θ-γ coupling and decreased θ frequency observed after ablation of NaV1.1. A, Schematic of the computational network consisting of five EX, OLM, and fast-spiking BC neurons each with all-to-all connectivity. Full details about the model are provided in Extended Data Figure 10-1. B, Example traces of raw, θ-filtered (θ; 5–10 Hz), low γ-filtered (low γ; 40–100 Hz), and high γ (100–200 Hz) filtered LFP from the model with normal (baseline) and impaired NaV function in both OLM neurons (35 and 28 mS/cm2, respectively) and BC neurons (10 and 6 mS/cm2, respectively). C, Average LFP power spectra from the same simulations used in B with normal (black) or impaired (red) inhibitory NaV function, showing intact γ power and a reduced θ peak frequency in the latter (C, inset). D, Phase-amplitude comodulograms of these simulations show impaired θ-modulated γ. E, Stepwise impairment of NaV function specifically in BC neurons (blue, left plot) or OLM neurons (green, right plot) resulted in a leftward shift in peak θ frequency for OLM neurons (black: 35 mS/cm2, bright green: 28 mS/cm2, steps of 1.75 mS/cm2), but not for BC neurons (black: 10 mS/cm2, bright blue: 6 mS/cm2, steps of 1 mS/cm2). F, Average θ-γ MI calculated for total γ (left), low γ (center), and high γ (right) with decreasing levels of NaV function in BC neurons (blue), OLM neurons (green) and both BC and OLM neurons (black).
Extended Data Figure 10-1
Currents and fluxes regulating membrane potential, ion concentrations, and neurotransmitter homeostasis for the computational neuronal network model. Download Figure 10-1, PDF file.
To delineate the specific roles of the two populations of inhibitory interneurons in this model, we next systematically decreased NaV function in either BC or OLM neurons. We found that the leftward shift in θ power results from NaV dysfunction in the OLM population (Fig. 10E). Decreasing NaV function in BC or OLM neurons had opposite effects on θ-γ coupling: decreasing NaV function in BC neurons alone impaired θ-γ MI, while decreasing NaV function in OLM neurons increased MI (Fig. 10F). Simultaneous reduction of NaV function in both BC and OLM neurons resulted in a progressive decrease in θ-γ MI over the high γ (100–200 Hz) range, while a slight increase followed by a decrease occurred over the low γ (40–100 Hz) range. Thus, θ-modulated high γ is more strongly affected than low γ in this computational model, which reproduces the in vivo results obtained following local Scn1a knock-out.
Discussion
We here report in juvenile DS mice impaired θ-γ phase-amplitude coupling, that persisted in mice that developed spontaneous seizures but recovered to WT levels in mice that remained seizure-free. Impaired θ-γ coupling was improved following acute CBD treatment. Hippocampal or cortical-specific knock-out of Scn1a impaired cortical θ-γ coupling, which preceded spontaneous seizures. Attenuated modulation of high γ by θ in the hippocampus, which colocalized with fast ripple activity preceding the start of chronic epileptic activity, was reproduced in a computational model in which NaV dysfunction was selectively introduced in inhibitory neurons. Together, our findings indicate that θ-γ coupling may serve as an early indicator of inhibitory dysfunction and epileptogenesis in DS.
Dysfunction of inhibitory interneurons has been demonstrated in DS mouse models (Tai et al., 2014; Rubinstein et al., 2015; De Stasi et al., 2016; Goff and Goldberg, 2019), and may depend on age window or brain region (Favero et al., 2018; Almog et al., 2019). Reports on how inhibitory dysfunction affects network dynamics in DS mice appear more ambiguous. For instance, decreased δ power during non-REM sleep was reported in DS mice (Kalume et al., 2015), but changes in sleep-related δ power were not present in another DS model (Papale et al., 2013). Also, interneurons were hypoexcitable, but LFP power was unchanged in DS mice during the preepileptic period (De Stasi et al., 2016). Clinically, patients with DS often have a normal background EEG at one to two years of age, while seizures already manifest at this time (Bureau and Dalla Bernardina, 2011; Wirrell et al., 2017). We did not detect differences in θ or γ power during REM sleep and active wakefulness in DS versus WT mice. This contrasts with a previous study reporting decreased θ power during REM sleep in DS mice (Kalume et al., 2015), which may be explained by differences in mouse strain background and the cortical area studied. Also, we used a strict cutoff in θ/δ ratio to delineate epochs of REM sleep and active wakefulness, which may have rendered our analyses less sensitive to subtle differences in θ power. Nevertheless, our results indicate that θ-γ coupling is a more sensitive measure of network dysfunction in DS than power spectral density.
Knock-out of Scn1a appears to specifically affect inhibitory populations, with reduced sodium currents lowering their excitability (Yu et al., 2006; Ogiwara et al., 2007; Tai et al., 2014; Favero et al., 2018), although increased sodium currents have been reported in hippocampal excitatory neurons (Mistry et al., 2014). Such hyperexcitability of hippocampal excitatory neurons may be attributed to seizure-induced remodeling in DS mice (Dutton et al., 2017; Salgueiro-Pereira et al., 2019). Importantly, chronic epilepsy may impair θ-γ coupling, as demonstrated for the hippocampal CA1 region in a rat model of temporal lobe epilepsy (Lopez-Pigozzi et al., 2016). Similarly, impairments in θ-γ coupling were recently found in the hippocampus of rats following local knock-down of Scn1a, whereas no epileptiform activity was observed (Sakkaki et al., 2020). Although such activity may not be detected since recordings were performed during wakefulness, when IIS are rare (Kalume et al., 2015), NaV1.1 dysfunction was induced only in CA1 of the dorsal hippocampus (Sakkaki et al., 2020), which may prove insufficient to cause epileptiform activity. In our study, the majority of the hippocampus was targeted, such that subregional effects of NaV1.1 dysfunction cannot be discerned. DS mice and mice with local hippocampal or cortical knock-out of Scn1a showed early impairment of θ-γ coupling before, or in the absence of, chronic epileptic activity. This impairment was transient in seizure-free DS mice, but persisted in mice with spontaneous seizure activity. Thus, our data indicate that θ-γ coupling is impaired early in all DS mice, but persistence of impaired θ-γ coupling at later stages is dependent on seizure-induced network remodeling.
Over postnatal weeks 4–7, WT and DS mice showed an increase in cortical θ-γ coupling. In mice, this is a developmental period of high cortical plasticity that starts with maturation of fast-spiking parvalbumin-positive interneurons (Takesian and Hensch, 2013; Cardin, 2018). In seizure-free DS mice, we found that θ-γ coupling is impaired at P23, but normalized to WT levels at P28–P42. This evolution roughly corresponds to developmental data from cortical parvalbumin-positive interneurons that show decreased firing only at P18–P21, but not at P35–P56, suggesting that their dysfunction contributes to the onset, but not chronification, of epilepsy in DS (Favero et al., 2018). Here, local cortical heterozygous knock-out of Scn1a resulted in transient attenuation of θ-γ coupling, which may reflect such temporary inhibitory dysfunction. CBD, which increases inhibitory neurotransmission in DS mice (Kaplan et al., 2017), increased θ-γ coupling, supporting a link between inhibitory functioning and θ-γ coupling. Evidence that θ rhythms modulate a greater proportion of inhibitory interneurons when compared with excitatory neurons (Csicsvari et al., 1999; Sirota et al., 2008) suggests that inhibitory dysfunction particularly impacts θ-γ coupling. This is further attested by impaired modulation of hippocampal γ amplitude by θ phase in mice that lacked inhibition onto parvalbumin-positive interneurons, while θ power and peak frequency was reduced (Wulff et al., 2009). Similarly, we found a reduced θ peak frequency in DS mice, and following local NaV1.1 ablation. Decreased NaV1.1 expression was also observed in a mouse model of Alzheimer's disease, as well as in patients, providing an explanation for dysfunction of parvalbumin-positive interneurons (Verret et al., 2012). Since impaired θ-γ coupling has been demonstrated in Alzheimer's models (Ittner et al., 2014; Zhang et al., 2016), similar mechanisms may underlie such rhythmopathies in Alzheimer's disease and DS models.
In our computational model, decreases in θ peak frequency could be attributed to NaV dysfunction in OLM neurons. In animal models of temporal lobe epilepsy, hippocampal θ frequency is decreased (Dugladze et al., 2007; Kilias et al., 2018) and degeneration of OLM neurons, resulting in loss of dendritic inhibition, has been observed (Cossart et al., 2001). Thus, decreased θ frequency appears to be a general feature of the epileptic hippocampus, and dysfunction and/or loss of OLM neurons may underlie this phenomenon. Although (parvalbuminergic) BC neurons appear to maintain their function in the epileptic hippocampus (Cossart et al., 2001; Wittner et al., 2005), dysfunction of these neurons has been related to impaired θ-γ coupling (Lopez-Pigozzi et al., 2016), which we confirmed here by modeling BC NaV dysfunction. Together, our findings suggest that decreased θ frequency and θ-γ coupling reflect the impaired inhibitory functioning observed in DS mice. Network dynamics underlying cross-frequency coupling, however, are complex, and different factors may affect θ-γ coupling in other disease models.
We distinguished low and high γ bands, supported by evidence that mechanisms underlying these oscillations are different (Buzsáki and Wang, 2012). High γ (above 80–90 Hz) peaks at a different phase of the θ cycle than low γ (Belluscio et al., 2012), as confirmed by our data, and strongly correlates with neuronal spiking (Ray and Maunsell, 2011). Local NaV1.1 dysfunction impaired local θ-modulated high, but not low, γ oscillations, which was computationally replicated. This suggests that NaV1.1 dysfunction impairs modulation of local neuronal firing by θ, as previously shown for CA1 neurons following local NaV1.1 knock-down (Sakkaki et al., 2020). The critical role for fast-spiking BC in regulating spike timing by somatic inhibition (Cardin, 2018) may explain these early deficits in θ-high γ coupling.
Increased hippocampal θ-γ coupling is associated with memory recall in rodents (Tort et al., 2009) and humans (Axmacher et al., 2010; Heusser et al., 2016). Memory impairment is frequently observed in patients with temporal lobe epilepsy (Helmstaedter and Elger, 2009; Vlooswijk et al., 2010). Although we did not relate changes in θ-γ coupling to learning, previous studies indicate that context-dependent and spatial memory are impaired in DS mouse models (Han et al., 2012; Salgueiro-Pereira et al., 2019). Cognitive deficits were only present in adult DS mice that had experienced seizures at juvenile age (Salgueiro-Pereira et al., 2019), which may be related to persistent deficits in θ-γ coupling that we only observed in adult DS mice with juvenile-onset seizures. The hippocampus appears to critically contribute to spatial memory deficits and the epileptic phenotype in DS mice (Liauard et al., 2013; Stein et al., 2019). Our data suggest that impaired θ-γ coupling precedes seizures and associates unilaterally with fast ripple activity in the hippocampus, a marker for epileptogenesis (Jefferys et al., 2012). Further studies are required to translate these findings to other epilepsy models. Reductions in θ frequency in such models have been suggested to underlie cognitive deficits in epilepsy (Shuman et al., 2017; Kilias et al., 2018), but in our experiments occurred regardless of epileptic activity. Similarly, local knock-down of NaV1.1 in the medial septum, which provides input to the hippocampus, decreased θ frequency and disrupted spatial memory in rats (Bender et al., 2016) without inducing seizure activity (Bender et al., 2013), indicating that regions other than the hippocampus may contribute to impaired spatial memory and θ rhythmogenesis in DS. To recapitulate rhythmopathies in DS, the precision of computational models may therefore be improved by including inputs such as provided by the medial septum.
Footnotes
This work was supported by the Dutch National Epilepsy Foundation Grant 2017-10, Medical NeuroDelta, and the National Institutes of Health Grant R01AG053988.
The authors declare no competing financial interests.
- Correspondence should be addressed to Nico A. Jansen at n.a.jansen{at}lumc.nl
