Abstract
Fragile X syndrome (FXS) is the most common source of intellectual disability and autism. Extensive studies have been performed on the network and behavioral correlates of the syndrome, but our knowledge about intrinsic conductance changes is still limited. In this study, we show a differential effect of FMRP knockout in different subsections of hippocampus using whole-cell patch clamp in mouse hippocampal slices. We observed no significant change in spike numbers in the CA1 region of hippocampus, but a significant increase in CA3, in juvenile mice. However, in adult mice we see a reduction in spike number in the CA1 with no significant difference in CA3. In addition, we see increased variability in spike numbers in CA1 cells following a variety of steady and modulated current step protocols. This effect emerges in adult mice (8 weeks) but not juvenile mice (4 weeks). This increased spiking variability was correlated with reduced spike number and with elevated AHP. The increased AHP arose from elevated SK currents (small conductance calcium-activated potassium channels), but other currents involved in medium AHP, such as Ih and M, were not significantly different. We obtained a partial rescue of the cellular variability phenotype when we blocked SK current using the specific blocker apamin. Our observations provide a single-cell correlate of the network observations of response variability and loss of synchronization, and suggest that the elevation of SK currents in FXS may provide a partial mechanistic explanation for this difference.
Significance Statement
Fragile X syndrome leads to a range of intellectual disability effects and autism. We have found a differential effect of FMRP KO in different subsections of hippocampus where it caused an increased spiking in CA3 in juveniles and reduced spiking in CA1, in adults. We have also found that even individual neurons with this mutation exhibit increased variability in their activity patterns. Importantly, this effect emerges after 6 weeks of age in mice. We showed that a specific ion channel protein, the SK channel, was partially responsible, and blockage of these channels led to a partial restoration of cellular activity. This is interesting as it provides a possible molecular link between activity variability in single cells and reported irregularity in network activity.
Introduction
Fragile X syndrome (FXS) is an autism spectrum disorder arising from increased repeats of CGG trinucleotide in the FMR1 gene. This leads to silencing of the gene and hence to the absence of the FMRP protein (O’Donnell and Warren, 2002). FMRP is a transcription factor that affects multiple downstream genes including ion channels (Brown et al., 2001; Darnell et al., 2011). Studies have shown that the absence of FMRP affects the functioning of ion channels either by regulating the number of channels (Meredith et al., 2007; Strumbos et al., 2010; Gross et al., 2011; Lee et al., 2011a; Routh et al., 2013; Ferron et al., 2014; Zhu et al., 2018) or by directly interacting with the channels (Brown et al., 2010; Deng et al., 2013). Multiple studies have shown that potassium channels have a significant effect on spike precision (Fricker and Miles, 2000; Gittelman, 2006; Cudmore et al., 2010), and many of these potassium channels are transcription hits for FMRP. This leads to the hypothesis that FXS may alter the functioning of one or multiple potassium channels, leading to effects on spike precision.
In recent years a number of in vivo hippocampal recording studies have shown that there is poor correlation of spiking activity between cells, and abnormal theta phase–gamma phase coupling in FXS mice (Radwan et al., 2016; Arbab et al., 2018a,b; Talbot et al., 2018). In medial prefrontal cortex, variability in calcium (Ca2+) responses has also been observed, leading to impaired spike timing-dependent plasticity (STDP) (Meredith et al., 2007).These studies have led to the “discoordination hypothesis” for FXS (Talbot et al., 2018). This hypothesis states that neurons in FXS are uncorrelated and have aberrant network discharges. In apparent contradiction to this hypothesis, neurons showed hyperconnectivity and synchronization in cortical networks of FXS model mice (Testa-Silva et al., 2012; Gonçalves et al., 2013).
Synchronicity is an emergent property of a network and is a function of both network connectivity and intrinsic properties. Specifically, potassium conductance has been shown to have significant effects on spike precision and network synchrony (Fricker and Miles, 2000; Pfeuty et al., 2003; Deister et al., 2009; Cudmore et al., 2010; Gastrein et al., 2011; Hou et al., 2012). Modeling studies have also shown that conductance that mediates spike frequency adaptation helps to synchronize network firing (Crook et al., 1998). Ih currents, IM currents, and SK currents have been shown to mediate adaptation for phase precession in hippocampal neurons, (Orbán et al., 2006; Kwag et al., 2014) spiking reliability in insect neurons (Gabbiani and Dewell, 2018) and also to maintain regular firing in robust nucleus of acropallium in zebra finches (Hou et al., 2012). Thus FXS mutations, through their effects on potassium channels, may lead to several of the observed network dysfunctions in FXS models.
In the present study, we tested alterations in intrinsic properties of the cell that might underlie the above network observations of uncorrelated activity. First, we observed that cells from FXS knock-out (KO) hippocampal slices were more unreliable, both at within-cell and between-cell levels, than controls. Consistent with the developmental manifestations of FXS, this effect was absent in juvenile mice and emerged only after ∼6 weeks. Decreased spiking in KO mice was correlated with increased variability. At the mechanistic level, we found that the spiking changes were due to elevated medium afterhyperpolarization (mAHP), which in turn was due to elevated SK currents. In contrast to these observations of decreased spiking in the CA1 region of adult animals, we observed increased spiking in the CA3 region of juvenile animals (Deng et al., 2019), but not in the CA1 region. Consistent with earlier observations (Deng et al., 2019), we did not see any significant change in the spatial distribution of SK channels in either CA3 or CA1. Finally, we tested the outcome of the blockers of SK and found that that they indeed led to a partial rescue of the observed phenotype of increased spiking variability and reduced timing precision in adult animals.
Materials and Methods
Animals
C57BL/6 strain of Fmr1tmCgr male mice were used for the experiments. All experimental procedures were approved by the National Centre for Biological Sciences ethics committee [Project ID: NCBS-IAE-2017/04(N)]. The animals were housed in the institute animal house where they were maintained on a 12 h light/dark cycle. The animals used were from an older animal group in the range of 6–8 weeks of age; the younger group was 3–4 weeks of age.
Slice preparation
Mice were anesthetized with halothane. Their head was decapitated after they were killed by cervical dislocation. Hippocampal slices were made in the ice-cold aCSF of the following composition: 115 mm NaCl, 25 mm glucose, 25.5 mm NaHCO3, 1.05 mm NaH2PO4, 3.3 mm KCl, 2 mm CaCl2, and 1 mm MgCl2. 400-μm-thick slices were made using a VT1200S vibratome and then incubated at room temperature for 1 h in the aCSF, which was constantly bubbled with 95% O2 and 5% CO2. Subsequently, the slices were transferred to the recording chamber where they were maintained at an elevated temperature of 30–34°C for the recordings.
Electrophysiology
CA1 neurons were identified under an upright differential interference contrast microscope (BX1WI microscope, Olympus) using a 40× objective (water immersion lens, 0.9 numerical aperture, LUMPLFLN, 40×). 2–4 MΩ pipettes were pulled from thick-walled borosilicate glass capillaries on a P-1000 Flaming Micropipette Puller (Sutter Instrument). The pipettes were filled with internal solution of the following composition for whole-cell current-clamp recordings: 120 mm potassium gluconate, 20 mm KCl, 0.2 mm EGTA, 4 mm NaCl,10 mm HEPES buffer, 10 mm phosphocreatine, 4 mm Mg-ATP, and 0.3 mm Na-GTP, at (pH 7.4 and 295 mOsm). For voltage-clamp recordings, the same composition of internal solution was used with the one change: 120 mm potassium gluconate was substituted with 120 mm potassium methylsulphate. Cells were recorded if they had a resting potential of ≤60 mV. We also required that they exhibit stable firing with little or no depolarization block for lower current inputs. Series resistance and input resistance were continuously monitored during the protocols, and the cell was discarded if these parameters changed by >25%.
Protocol for measuring spike variability and analysis
All spike variability and precision experiments were performed in current-clamp mode. A step input current stimulus of 150 pA DC for 900 ms was used for the majority of the recordings. In some cases, as indicated in the text, frozen noise and sinusoidal input currents were also used, riding on a baseline current step of 150 pA, and again for a duration of 900 ms. The SDs of noise used in noise protocols were 10, 25, 50, and 100 pA with a time cutoff of τ = 3 ms. For sinusoidal currents, SDs of 50 and 100 pA were used at 5 Hz. All the protocols were repeated for 25 trials. Within-cell spike variability (CVw) was computed as coefficient of variation (CV) in spike numbers across 25 trials. Spike numbers computed for all 25 trials were used to find the SD in spike numbers, which was then divided by the mean number of spikes to find CVw, as follows: (1)
where n is the cell index, σn is the SD across trials within cell n, and μn is the mean number of spikes across all the trials within cell n.
The array of floating point numbers (i.e., vector for CVw) was computed for all recorded KO and wild-type (WT) cells, respectively, and was used for comparing the responses (Fig. 1E). The distributions in CVw were compared using the Wilcoxon test.
The analysis of between-cell spiking variability (CVb) was derived from an article that reports experiments on network synchrony (Talbot et al., 2018). In contrast to our CVw analysis, where variability in spiking was assessed within each cell, CVb measures the variability across trials between cells. For computing CVb, spike numbers on corresponding trial numbers were compared between different cells. Specifically, we constructed a vector for each cell whose entries were the number of spikes for successive trials. These vectors were compared across different cells within the same genotype and were used to find CVb. The following formula was used for the computation: (2)
where i is the trial index, σi is the SD over trial i between cells, and μi is the mean number of spikes on trial i between all cells. The CVb vectors were computed for WT and KO populations, respectively (Fig. 1F), and compared using Wilcoxon statistics.
Measurement of spike timing precision
Measurement of spike timing precision: CVs
Spike times were calculated for all the spikes in all trials. The interspike interval (ISI) was computed for the first two spikes and for spikes that bracketed a time of 500 ms after the start of the step input protocol (Fig. 1I). For each cell m, we obtained the coefficient of variation, CVs0_m and CVs500_m for the first ISI and the ISI at 500 ms, respectively, as follows: (3)
where m is the cell index, σm is the SD of ISI across trials within cell m, and μm is the mean ISI across all the trials for cell m. We calculated CVs500_m in a similar manner.
By repeating this calculation over all cells, we obtained the following four vectors: CVs0_WT, CVs500_WT, CVs0_KO, and CVs500_KO. The distributions in these four vectors were compared using a two-way ANOVA (Fig. 1J).
Jitter slope
The jitter slope analysis was taken from the study by Bacci and Huguenard (2006). Jitter was defined as σti, which is the SD of spike timing for spike number i, compared over all trials for a given cell. In this manner, the vector σt was obtained by considering all the spike numbers. Note that not all trials had the same number of spikes, so the final few spike numbers may have had fewer samples in σt.
The jitter slope was obtained as a regression fit for the vector σt, with the x-axis defined by the corresponding vector for mean timings μt. This slope was calculated for spikes falling in the last 400 ms of the current pulse since this is where the SK channel is activated. As a pre hoc criterion, cells were required to have a minimum of two spike time points to be selected for this calculation. Note that this criterion excluded several cells that had strong suppression of firing toward the end of the trial, presumably due to SK channel activation. Thus, this measure was compromised by being insensitive to precisely those cells with the largest effects in the KO.
Measurement of frequency versus current curves
The measurement of frequency versus current (f–I) curves was performed in current clamp. For plotting f–I curves, each cell was given input currents from −100 to 400 pA at its resting membrane potential (RMP). The mean ± SEM of the number of spikes elicited for each input current was plotted on the y-axis, and the corresponding input current on the x-axis.
AHP measurement
AHPs were measured by holding the cell at −70 mV in current-clamp mode, using some holding current. In each trial, the holding current was set to 0 for 200 ms, and then step input for 100 ms was delivered. The current during this step input was adjusted so that the cell produced five spikes in 100 ms. This stimulus was repeated for 15 trials, and the voltage traces of all the trials were averaged. AHPs were analyzed as the hyperpolarization after the spike train. mAHP was computed as the average hyperpolarization 50 ms after the spike train, and slow AHP (sAHP) was computed 200 ms after the spike train (see Fig. 5A,B).
Immunofluorescence experiments and analysis
The 3- to 4-week-old or 6- to 8-week-old animals were anesthetized, and brains were dissected out. Brains were fixed in 4% PFA for 48 h at 4°C. Postfixation, brains were washed with PBS, and 30-μm-thick slices were cut on a VT12000S vibratome. Slices were kept in cryoprotectant solution and stored at −20°C for further use. Slices were resuspended in PBS solution for the experiment. 50 mM ammonium chloride (NH4Cl) was used to quench autofluorescence from the samples. After multiple steps of washing and permeabilization in 0.3% PBS with Triton X-100, cells were incubated with 5% normal goat serum (NGS) for blocking for 1–2 h. Followed by washes, slices were incubated in primary antibody for the SK2 channel, overnight at 4°C. Primary antibody was purchased from Alomone Labs (catalog #APC–028) and were used in 1:300 dilution with 2.5% NGS. Further, slices were incubated in secondary antibody at 1:500 dilution for 1–2 h at room temperature. After several washes and DAPI staining (Sigma-Aldrich) for 15 min in PBS, slices were mounted in Prolong Gold Mounting Media and imaged with an Olympus FV3000 confocal microscope.
Measurement of M and Ih currents
For measuring the Ih sag produced in voltage change, cells were held at −70 mV and given a step hyperpolarization of 250 pA for 500 ms in current-clamp mode. Voltage deflection was calculated by finding the percentage change between maximum and steady-state voltages during the hyperpolarizing pulse [(Vmax − Vss)/Vmax] * 100. Each cell was presented with 15 trials, and Ih sag was averaged across all the trials to get a value for a cell.
TTX (0.5 μm) was added in the bath while measuring M currents in voltage clamp. For measuring M currents, cells were held at −20 mV, and M currents were deactivated in steps from −20 to −80 mV with a difference of 10 mV in a 900 ms time duration. The difference between the steady-state current after the start of the protocol (from the last 250 ms) and the baseline of the cell was quantified as the current elicited due to the protocol. Specific M current blocker XE991 (30 μm) was perfused to get the postblocker trace. The difference between pre- and post-XE991 perfusion was quantified as the M current.
Measurement of SK currents
TTX (0.5 μm) was added into the bath while measuring SK currents in voltage clamp. For measuring SK currents, cells were held at −55 mV while increasing voltage steps were given from −55 to +35 mV, in increments of 10 mV, for 800 ms duration. Each voltage step was preceded with a small hyperpolarization from −55 to −65 mV for 100 ms. Cell capacitance compensation was performed during the recordings. The SK currents were estimated from the amplitude of the outward current 25 ms after the voltage was returned to the holding potential as the initial part of the trace will be contaminated with capacitance currents (see Fig. 6A,B). Apamin perfusion (100 nm) was performed to identify these currents. The same protocol was used post-apamin perfusion. The difference in current between pre-apamin and post-apamin perfusion was quantified to find the apamin-sensitive SK current.
Data acquisition and analysis
Recordings were performed using an Axopatch 200 B Amplifier and digitized using a Digidata 1400A digitizer. Current-clamp recordings were filtered at 10 kHz and sampled at 20 kHz. Voltage-clamp recordings were conducted at a gain of 5, filtered at 3 kHz, and sampled at 20 kHz. All analyses were performed using custom-written code in MATLAB (R2013a) and GraphPad Prism 6. Statistical tests were performed after checking the normality of the data. We plotted the data to check for normality. Even if the data looked normal, a Wilcoxon rank sum test was used for most analysis for more stringent results. For multiple comparisons, ANOVA was used with Sidak’s post hoc test or t test with Bonferroni correction as post hoc tests. All data are plotted as the mean ± SEM.
Drugs and chemicals used
All toxins were purchased from Sigma-Aldrich except for TTX, which was from Hello Bio. TTX (stock solution of 0.5 mm) and apamin (stock solution of 1 mm) were made with Milli-Q water. XE991 (stock solution of 10 mm) was made in DMSO. All stocks were stored at −20°C and were diluted to the working concentration immediately before the experiments.
Experiment design and statistical analysis
We have used hippocampal slice preparation from male mice in this study, where CA1 and CA3 regions were used for electrical recordings. Different current inputs were injected into the cell soma via recording electrode, and elicited spikes were recorded in current-clamp mode. Spiking variability was assessed by finding the CV for within a cell and between cells. Spike precision was assessed by comparing CV ISI occurring in the initial part of the trail with the ISI occurring around 500ms post commencement of the trial. SK currents were measured in voltage-clamp mode and were found to be partially responsible for increased variability. In addition, we found a differential effect of FMRP KO in different subsections of hippocampus. In young animals, we observed KOs having hyperexcitability for CA3 cells with no difference in CA1 cells. In older animals, we saw that the hyperexcitability seen in CA3 goes away while there is an emergence of hypoexcitability in CA1 cells. All analysis and plotting was performed in MATLAB and GraphPad 6. The normality of the data was checked by plotting it. Even if the distribution followed a normal distribution profile, a Wilcoxon test was used. For multiple comparisons, two-way or three-way ANOVA with Sidak’s post hoc test or t test using a Bonferroni correction were used as the post hoc test.
Results
FMR1 KO cells exhibit increased variability in spike numbers and impaired spike precision
We performed whole-cell patch clamp recordings in brain slices from CA1 cells of dorsal hippocampus in male WT and KO littermate animals (6–8 weeks old). Somatic current injection was used to inject various waveforms of current (Fig. 1A). These included a square current step, frozen noise riding on a current step, and a sinusoidal wave form riding on a current step (Fig. 2). In each case, we measured the variability in spiking over the 900 ms of current injection, in the following two respects: first, the variability within a neuron, over repeated trials (CVw); and, second, the variability between neurons, comparing matching trials between cells (CVb; see Materials and Methods).
We computed the co-efficient of variation (CVw) of spike numbers for each of the current stimulus protocols (see Materials and Methods). We observed that the CV of spike numbers across trials within cells (CVw) was significantly higher for KO cells compared with WT cells, but only for the step current protocol (CVw for WT, 0.038 ± 0.005, n = 30 cells; CVw for KO, 0.057 ± 0.006, n = 29 cells; p = 0.0048; Wilcoxon rank sum test; Fig. 1B–D). For other protocols, where frozen current noise or sinusoidal waveforms were riding on top of a current pulse, spikes were very reliable over different trials for both WTs and KOs (Fig. 2A–E). We then estimated the variability between cells (CVb). To do this, we compared spike numbers between different WT and KO cells by matching corresponding trials, such that we compared spiking between the first trials of WT between all the WT cells, second trials of WT between all cells, and so on forth, within the same genotype. Similar comparisons were made for KO cells as well (see Materials and Methods). From this procedure, we found that between-cell variability was significantly elevated in KO cells compared with WT cells (CVb for WT cells, 0.29 ± 0.003, n = 30 cells; CVb for KO cells, 0.41 ± 0.002, n = 29 cells; p < 0.001; Wilcoxon rank sum test; Fig. 1F).We plotted the CVw of spike numbers within the cell for WT and KO populations and observed that KOs have a much wider spread compared with WT (SD of the CVw population for WT cells = 0.026, n = 30 cells; SD of the CVw population for KO cells = 0.03, n = 29 cells; Fig. 1G,H). This implied that KO cells have a more heterogeneous population than WT cells. In KO mice, some cells were reliable (i.e., they produced a similar number of spikes across multiple trials but some did not). Based on this analysis (Fig. 1G,H), we found that 51.7% of KO cells had CVw values that did not overlap with those of WT cells.
Spike timing precision was computed by finding CVs in ISIs (Materials and Methods) and by jitter slope analysis (Materials and Methods). We compared CVs0 with CVs500 to assess the contribution of SK channels, which are shut at t = 0, but open by t = 500. We chose not to use the last ISI for this comparison because the timing of the last spike was very variable in KO cells. We then compared CVs at t = 0 and t = 500, between WT and KO cells using two-way ANOVA. We found that there was no significant difference between WT and KO cells for CVs0. However, at t = 500, we found that CVs500_KO was significantly elevated compared with CVs500_WT (CVs0_WT = 0.11 ± 0.014; CVs0_KO = 0.09 ± 0.007; CVs500_WT = 0.13 ± 0.01; CVs500_KO = 0.23 ± 0.03; F(1,113) = 5.78; p = 0.018 for WT versus KO differences; two-way ANOVA, t test with Bonferroni corrections; Fig. 1J). The jitter slope analysis did not report a significant difference in spike imprecision between WT and KO. We used the last 400 ms of the trial for the analysis as this part showed peak SK current activity. The analysis was weakened because of our criterion that a cell should have minimum of two spikes for regression fit. Hence, several KO cells were rejected, although these were probably the ones that probably had the strongest effect of elevated SK channel activity (Fig. 1K,L).
We observed increased within-cell variability (i.e., CVw), between-cell variability (i.e., CVb) of spike numbers and increased CVs500 of spike timings, which are all indicative of abnormal coding in KO neurons. Thus, these initial recordings showed that there was increased heterogeneity in spiking in CA1 pyramidal neurons of KO mice.
KO mice exhibit increased CVb for frozen noise stimuli
Neuronal coding takes place against a noisy background of continuous synaptic input (Faisal et al., 2008). Frozen noise stimuli have therefore been used in several studies to examine single-neuron computations and are considered more physiological (Stevens and Zador, 1998; Carandini, 2004; Fujisawa et al., 2004; Han and Boyden, 2007). We injected noise waveforms in the soma using a design similar to that of Mainen and Sejnowski (1995). The stimulus waveforms were mainly either of random noise or sinuosoidal input currents (for details, see Materials and Methods). The RMP was controlled by maintaining constant current injection into the cell while recording. All cells were held at a constant potential of −65 mV for all of the protocols. Cells were accepted only if the fluctuations in their potential while recordings were smaller than ±2.5 mV. For noisy inputs, spikes across multiple trials aligned very precisely for both WT cells (Fig. 2A,C) and KO cells (Fig. 2B,D). On computation of CV of spike numbers across trials within a cell (i.e., CVw), there were no statistically significant differences for any of the protocols (Fig. 2E). However, between-cell variability differences (CVb) were still significantly elevated for KO cells, for all stimulus protocols (Fig. 2F). We performed a three-way ANOVA for all of the protocols on the number of spikes produced in respective protocols for WT and KO cells. The three factors were as follows: across trials, across cells, and WT versus KO cells. ANOVA results were significant across cells (p < 0.001) and for WT versus KO (p < 0.001), but not for across trials (p = 1) for all the protocols. Thus, between-cell variability (CVb) reveals robust phenotypic differences, which are also observed between WT and KO cells for the input noise stimulus. The same is not true for across-trial within-cell (CVw) differences.
Spike imprecision diverges between KO and WT only in adult mice
Fragile X is known to have an age-dependent effect where some disease phenotypes get worse with age (Larson, 2005; Cornish et al., 2008; Choi et al., 2010). We asked whether similar effects of age might manifest in the phenotype of the firing imprecision of KO cells. We therefore repeated our experimental protocols in young (3- to 4-week-old) male WT and KO littermate mice. Regular spike firing was observed in both KO and WT cells for step input current injection (Fig. 3A–C). As before, we assessed the variability of spike numbers across multiple trials within the cell (CVw).There was no significant difference between WT and KO cells in these younger mice (CVw for WT, 0.041 ± 0.006, n = 26 cells; CVw for KO, 0.03 ± 0.002, n = 19 cells; p = 0.81, Wilcoxon rank sum test; Fig. 3A–D). Spike precision was assessed as above by finding the CVs500. By this measure too, no significant difference in the precision of spike timings was found for KO compared with WT (CVs500_WT younger, 0.133 ± 0.02, n = 26 cells; CVs500_KO younger, 0.103 ± 0.01, n = 19 cells; p = 0.82, Wilcoxon rank sum test; Fig. 3E). We used a two-way ANOVA interaction model for assessing age-dependent differences in spiking variability across different trials in a cell (CVw; F(1,100) = 4.97, p = 0.028 for WT and KO differences; F(1,100) = 8.05, p = 0.005 for interaction factor; two-way ANOVA with interaction model, t test with Bonferroni corrections post hoc; Fig. 3D). We also estimated differences in spike timing precision across different trials in a cell (CVs500; F(1,99) = 7.47, p = 0.007 for age factor; F(1,99) = 9.46, p = 0.003 for interaction factor; two-way ANOVA with interaction model, t test with Bonferroni corrections; Fig. 3E). Spike timing precision differences were also assessed by jitter slope analysis (F(1,92) = 6.83, p = 0.011 for interaction factor; F(1,92) = 2.2, p = 0.14 for age factor; two-way ANOVA with interaction model, t test with Bonferroni corrections; Fig. 3F). ANOVA indicated a significant interaction between the two factors age and genotype of the animal. Using post hoc tests, we concluded that for KO cells there is a worsening of phenotype for both variability in spiking number (CVw) across multiple trials within a cell and in spike timing precision (CVs500 of ISI) across multiple trials within a cell (Fig. 3D–F; Table 1, detailed statistics). However, the same was not seen in WT cells. A two-way ANOVA was performed on the number of spikes produced in the 900 ms stimulus period used in the step depolarization protocol (Fig. 1A). In this ANOVA, the two factors were young versus old animals and WT versus KO. There was no significant change between the ages of the animals for both WT and KO (F(1,10) = 0.31; p = 0.6; two-way ANOVA, t test with Bonferroni correction), but there was a significant reduction in spike number for KO animals compared with WT animals (F(1,101) = 4.7; p = 0.032; two-way ANOVA, t test with Bonferroni correction; Fig. 3G). However, the f–I curves were not significantly different between WT and KO animals in either in the older (8 weeks) or younger (4 weeks) age groups, when the slope of the f–I curves (gain) was used as a readout (mean gain for f–I curve for older WT animals, 0.16 ± 0.007 spikes · s−1 · pA−1, n = 27 cells; mean gain for f–I curve for older KO animals, 0.16 ± 0.008 spikes · s−1 · pA−1, n = 25 cells; mean gain for f–I curve for younger WT animals, 0.133 ± 0.01 spikes · s−1 · pA−1, n = 26 cells; mean gain for f–I curve younger KO animals, 0.136 ± 0.013 spikes · s−1 · pA−1, n = 19 cells; F(1,94) = 0.03, p = 0.86 for WT and KO differences; two-way ANOVA, t test with Bonferroni correction; Fig. 3I,J). When comparisons were made of f–I curves of WT and KOs cells, there was reduced excitability seen for some input currents in older animals (mean number of spikes for an input current of 125 pA in WT animals, 22.07 ± 1.74; n = 27 cells; mean number of spikes for input current of 125 pA in KO animals, 15.92 ± 2.1; n = 25 cells; F(1,1070) = 31.03, p < 0.001; two-way ANOVA; Fig. 3L), but not in younger animals (mean number of spikes for input current of 125 pA in WT animals, 21.8 ± 2.9; n = 26 cells; mean number of spikes for input current of 125 pA in KO animals, 22.8 ± 2.9; n = 19 cells, F(1,923) = 0.92, p = 0.33; two-way ANOVA, t test with Bonferroni correction as post hoc; Fig. 3K). Other properties like spike threshold were not significantly different between WT and KO cells (spike threshold for the first spike in WT, −53.6 ± 0.56 mV, n = 30 cells; spike threshold for first spike KO, −53.8 ± 0.62 mV, n = 29 cells; spike threshold for the last spike WT, −41.9 ± 0.66 mV, n = 30 cells; spike threshold for the last spike KO, −43.5 ± 0.66 mV, n = 29 cells; F(1,115) = 1.98, p = 0.16, two-way ANOVA, t test with Bonferroni correction as post hoc test. The spread in the number of spikes produced by KO neurons was greater than observed in WT neurons, further substantiating the point of increased variability in KO neurons (SD of spike numbers for WT = 7.2; n = 750 spikes from 30 cells 25 trials each; SD of spike numbers for KO = 8.1; n = 725 spikes from 29 cells 25 trials each; Fig. 3H). We performed a three-way ANOVA in which the factors were as follows: across trials, across cells, and WT versus KO for older age group data. The ANOVA results were highly significant for the variability in spike numbers between cells (CVb; F(29,1420) = 32.04, p < 0.001) and differences in spike numbers between WT and KO for older animals (F(1,1420) = 209.9, p < 0.001), but across-trial comparisons did not exhibit a significant difference (F(24,1420) = 0.15, p = 1).Thus, there is an increase in spike variability with lowered spiking precision for adult KO mice compared with juvenile KO mice. This divergence that is observed for KO animals with age does not happen for WT animals.
Differential effect of FMRP KO in different brain regions
Mouse models of FXS are known to exhibit increased excitability (Contractor et al., 2015; Deng and Klyachko, 2016; Luque et al., 2017). In contrast to previous studies, we observed that excitability was reduced in CA1 neurons in adult mice, but not in juveniles. Deng et al. (2019) have observed increased CA3 excitability in juveniles, which was in contrast to what we observed in CA1, namely, that there is no significant effect on excitability of CA1 pyramidal cells at this time point. To address this contrasting result, we repeated the protocol from Deng et al. (2019) in our FXS model and compared outcomes between CA3 and CA1 (Fig. 4). In agreement with their study, we observed significantly increased spiking in CA3 cells of FXS mice for some of the input holding currents (number of spikes in 60 s for −48 mV holding potential for WT, 131.9 ± 17.42, n = 10 cells; number of spikes in 60 s for −48 mV holding potential for KO, 186.3 ± 27.27, n = 10 cells; F(1,72) = 5.35, p = 0.023 for difference between WT and KO; two-way ANOVA, Sidak’s post hoc test; Fig. 4A,B). We performed similar experiments in CA1 pyramidal cells, using the same protocol as in Figure 4, A and B. We did not find any difference in excitability between WT and KO (number of spikes in 60 s for −48 mV holding potential for WT, 355.58 ± 40.93, n = 12 cells; number of spikes in 60 s for −48 mV holding potential for KO, 358.67 ± 22.56, n = 12 cells; F(1,88) = 1.06, p = 0.31 for difference between WT and KO; two-way ANOVA, Sidak’s post hoc test; Fig. 4C,D). On seeing this differential effect of FMRP KO in different regions of hippocampus, we were interested to see whether this might be a function of age. Therefore, we recorded similar protocols from CA3 and CA1 of older animals. We found that the hyperexcitability seen in younger CA3 cells goes away, leading to no significant difference in the spiking activity in these cells for WT versus KO (number of spikes in 60 s for −48 mV holding potential for WT, 313.6 ± 83.8, n = 10 cells; number of spikes in 60 s for −48 mV holding potential for KO, 181.1 ± 54.9, n = 10 cells; F(1,75) = 1.16, p = 0.28 for difference between WT and KO; two-way ANOVA, Sidak’s post hoc test; Fig. 4E,F). In agreement with our previous results, as shown in Figure 3, we observed reduced excitability in CA1 KO cells compared with WT cells (number of spikes in 60 s for −50 mV holding potential for WT, 135.6 ± 35.2, n = 8 cells; number of spikes in 60 s for −50 mV holding potential for KO, 54 ± 20.4, n = 8 cells; F(1,59) = 7.42, p = 0.009 for difference between WT and KO; two-way ANOVA, Sidak’s post hoc test; Fig. 4G,H). Using these experiments, we were able to replicate previous published results of Deng et al. (2019). In addition, we found that there were contrasting and age-dependent effects of FMRP KO on the spiking properties of CA1 and CA3 pyramidal neurons.
mAHP currents are elevated in older KO animals compared with WT but no significant changes in Ih and M currents
We next investigated the mechanisms for increased variability in KO animals, on the basis of the above observations of increased variability and fewer spikes (Fig. 3D,G,L) in the later part of the current pulse. In particular, we hypothesized that mAHP and sAHP currents were likely candidates to mediate such slower-onset effects. We asked whether these currents might be altered in FXS cells.
To test this hypothesis, we compared AHPs for WT versus KO in older and younger animals. Cells were held at −70 mV to record AHPs. We quantified AHPs by examining the hyperpolarization elicited with a spike train of five spikes in 100 ms. Hyperpolarization was quantified at mAHP peak (50 ms after the spike train) and at sAHP peak (200 ms after the spike train). We found that mAHP was elevated only in older KOs but not in younger KOs when compared with their WT littermates (mAHP for older WTs, −0.18 ± 0.14, n = 12 cells; mAHP for older KOs, −0.71 ± 0.14; n = 13 cells; mAHP for younger WTs, −0.6 ± 0.19, n = 14 cells; mAHP for younger KOs, −1.04 ± 0.4, n = 18 cells; F(1,46) = 4.8, p = 0.033 for WT and KO differences; two-way ANOVA, t test with Bonferroni correction; Fig. 5C, Table 1, detailed statistics) . We found similar results for sAHP measurements as well. We found that sAHP was significantly elevated only for older KO animals compared with their WT littermates but not significantly different for younger KO animals compared with their younger WT littermates (sAHP for older WTs, 0.05 ± 0.09, n = 12 cells; sAHP for older KOs, 0.4 ± 0.11; n = 13 cells; sAHP for younger WTs, 0.07 ± 0.12, n = 14 cells; sAHP for younger KOs, 0.5 ± 0.17, n = 18 cells; F(1,54) = 8.84, p = 0.004 for WT and KO differences; two-way ANOVA, t test with Bonferroni correction; Fig. 5D). As there was no significant difference in spike variability and spike time precision in younger KO animals, we hypothesized that the observed spiking variability in older KO animals might be due to elevated mAHPs, which has a strong change with age in KO animals. Further, the stimulus period used in the previous experiments (Fig. 1) was of 900 ms, which is the peak time of activation for mAHP mediating ion channels.
The mAHP in CA1 cells is mediated by Ih, M, and SK currents. A previous study has shown that Ih voltage sag is not significantly altered in soma of FXS cells (Brager et al., 2012). M currents are known to affect the RMP and spike adaptation index (Madison et al., 1987; Guan et al., 2011; Nigro et al., 2014; Hönigsperger et al., 2015). Neither RMP nor spike adaptation index were significantly changed in FXS cells, indicating that M currents might not be responsible for the observed phenotype. To further test this idea, Ih voltage sag and M currents were recorded from soma of CA1 cells and compared between WT and KO cells for adult animals. As there was no difference in mAHP at juvenile stage for KO cells, we did not measure their Ih voltage sag and M currents. Ih channel activity was recorded using a hyperpolarization pulse. The voltage sag produced in response to the pulse, was quantified according to a method adapted from Brager et al. (2012; see Materials and Methods; Fig. 5E). As observed previously, there was no significant difference in Ih sag between WT and KO CA1 cells, at somatic level (percentage of sag for WT, 17.7 ± 1.64, n = 10 cells; percentage of sag for KO, 18.87 ± 2.51, n = 10 cells; p = 0.91, Wilcoxon rank sum test; Fig. 5F).
To measure M currents, we used the protocol adapted from Shah et al., 2002 (Materials and Methods). According to the protocol, the cells were held at −20 mV and M currents were deactivated in steps from −20 to −80 mV. M current-specific blocker XE991 was used to obtain the postblocker trace. Postblocker trace was subtracted from preblocker trace to obtain XE991-sensitive M currents (Fig. 5G). We found no significant differences in M currents between WT and KO CA1 cells. (mean M current for step −80 mV for WT, −40.9 ± 6.33, n = 19 cells; mean M current for step −80 mV for KO, 45.94 ± 6.87, n = 19 cells; F(1,258) = 2.54, p = 0.11 between WT and KO; two-way ANOVA, t test with Bonferroni correction; Fig. 5H). Overall, these observations suggest that neither Ih voltage sag nor M currents contribute to the observed elevation in mAHP in KO animals.
SK currents are elevated in older KO animals compared with WT
SK currents are also known to contribute substantially to mAHP (Sah and Clements, 1999; Stocker et al., 1999). It has also been shown that Ca2+ influx from VGCCs is elevated in FXS cells (Deng et al., 2011, 2013), indicating elevated Ca2+ levels intracellularly and hence increased SK channel activation. Thus, we hypothesized that SK currents are elevated in FXS cells, leading to the observed elevated mAHP in these cells. To test this, we used a procedure from the study by Carlen (1997). Outward currents were evoked using step depolarization of the cell from −55 to +5 mV and above (Fig. 6A,B), in voltage-clamp mode. Apamin was bath applied to obtain post-apamin perfusion trace. The difference in current before and after apamin perfusion was quantified as SK currents. These currents were significantly higher in KO animals compared with WT animals for the adult age group (SK current for older WTs at +15 mV step, 54.2 ± 9, n = 13 cells; SK current for older KOs at +15 mV step, 131.73 ± 25.5, n = 13 cells; F(1,240) = 33.3, p < 0.001 between WTs and KOs, two-way ANOVA, t test with Bonferroni correction; Fig. 6C). However, this phenotype was seen only in adult WT versus adult KO animals. When similar experiments were performed in juveniles of both genotypes, there was no significant difference in SK currents (SK current for younger WTs at +15 mV step, 61.93 ± 11.9, n = 11 cells; SK current for younger KOs at +15 mV step, 51.2 ± 8.9, n = 9 cells; F(1,189) = 0.51, p = 0.5 between WTs and KOs; two-way ANOVA, t test with Bonferroni correction; Fig. 6D). To delineate the developmental change in SK channels levels between WT and KO animals, comparisons were made within the same genotype between different age groups. There was no significant change in the levels of SK currents with age for WT animals (SK current for older WTs at +15 mV step, 54.2 ± 9, n = 13 cells; SK current for younger WTs at +15 mV step, 61.93 ± 11.9, n = 11 cells; F(1,229) = 0.01, p = 0.9 between WTs and KOs; two-way ANOVA, t test with Bonferroni correction; Fig. 6E). On the contrary there was a significant increment in levels of SK currents for KO animals as they progressed from juvenile to adult stage (SK current for older KOs at +15 mV step, 131.73 ± 25.5, n = 13 cells; SK current for younger KOs at +15 mV step, 51.2 ± 8.9, n = 9 cells; F(1,209) = 21.02, p < 0.001 between WTs and KOs; two-way ANOVA, t test with Bonferroni correction; Fig. 6F). Together, these observations show that there is an increase in the levels of SK currents with age in KO animals, and we note that this correlates with the worsening of spike reliability between adult and juvenile KO animals, as described in Figure 3.
SK channel expression is unchanged in CA1 and CA3 regions for KO compared with WT animals
We have observed that SK current levels are elevated in KO CA1 cells compared with WT cells (Fig. 6). This elevation might arise due to an increase in the levels of SK channels or due to a change in SK channel kinetics. The latter might arise due to a change in the kinetics of the SK channel itself or due to changes in the activity of other channels that have an effect on SK channels. The putative candidates might be Ca2+ channels known to be affected in FXS (Meredith et al., 2007; Ferron et al., 2014; Brown et al., 2001; Contractor et al., 2015). Previous published results have not shown SK channels to be targets of FMRP for transcription regulation, decreasing the likelihood that the levels will be affected in the CA1 region (Darnell et al., 2011; Deng et al., 2019). To specifically test this point, immunofluorescence experiments were used to assess the expression of SK channels. We performed in situ immunolabeling to measure the intensity of the SK2 isoform in CA3 and CA1 regions of the hippocampus using a specific antibody (Materials and Methods; Fig. 7A). We analyzed levels of SK2 protein in different regions and subsections of the hippocampus for both WT and KO animals. Similar experiments were performed on both old (adult) and young (juvenile) slices. In young animals, we have seen an increased number of spikes in CA3 with no significant change in CA1 (Fig. 4). However, the levels of SK2 protein were not significantly altered in KO slices for either CA3 (N = 4 animals for WTs; N = 4 animals for KOs; F(1,24) = 0.48, p = 0.5 for WT vs KO for CA3; two-way ANOVA, Sidak’s post hoc test; Fig. 7B) or CA1 (N = 4 animals for WTs; N = 4 animals for KOs; F(1,24) = 0.21, p = 0.6 for WT vs KO cells for CA1; two-way ANOVA, Sidak’s post hoc test; Fig. 7C). For old animals, we have seen a reduction in spike number in the CA1 region (Fig. 3G,L), but no changes in SK2 levels could be found in either CA1 or CA3 (CA1: N = 4 animals for WTs; N = 4 animals for KOs; F(1,24) = 0.25, p = 0.61 for WT vs KO for CA1; two-way ANOVA, Sidak’s post hoc test; Fig. 7E; CA3: N = 4 animals for WTs; N = 4 animals for KOs; F(1,24) = 0.3; p = 0.6 for WT vs KO for CA3; two-way ANOVA, Sidak’s post hoc test; Fig. 7D). Overall, these labeling studies suggest that the KO does not affect SK expression in either CA1 or CA3. The latter observation is in agreement with previous results (Deng et al., 2019).
Apamin partially rescues the phenotype of increased variability in KO animals
Having shown that elevated and variable SK currents play a role in cellular firing variability, we then asked whether blockage of the SK current might rescue the phenotype of increased variability. We observed that WT and KO populations responded differently on blocking SK currents using apamin. In WT cells, the mean of within-cell variability (CVw) between trials increased (CVw for WT cells pre-apamin perfusion, 0.03 ± 0.002, n = 13 cells; CVw for WT cells post-apamin perfusion, 0.04 ± 0.005, n = 13 cells; Fig. 8A,B,E). However, for KO cells, the mean of within-cell variability (CVw) between trials decreased (CVw for KO cells pre-apamin perfusion, 0.06 ± 0.008, n = 13 cells; CVw for KO cells post-apamin perfusion, 0.05 ± 0.009, n = 13 cells; Fig. 8C,D,E). Apamin treatment had no significant effect on within-cell variability (CVw) for either WT or KO. However, it leads to a partial rescue such that post-apamin perfusion KO CVw values were not significantly different from pre-apamin perfusion WT CVw values. A complete rescue by apamin would have resulted in a significant reduction of the CVw of KO cells (which did not happen) as well as a convergence of variability between KO and WT cells (which did) (F(1,49) = 0, p = 0.96 for the effect of apamin perfusion; F(1,48) = 8.74, p = 0.005 for WT and KO differences; two-way ANOVA, t test with Bonferroni correction; Fig. 8E). We performed similar analysis for between-cell (CVb) variability. Apamin perfusion led to a partial rescue in the CVb parameter as well, such that there was a significant reduction in the CVb of KO cells on apamin perfusion, but no convergence of variability between WT and KO cells (CVb for WT cells pre-apamin perfusion, 0.413 ± 0.002; CVb for WT cells post-apamin perfusion, 0.414 ± 0.003; CVb for KO cells pre-apamin perfusion, 0.47 ± 0.003; CVb for KO cells post-apamin perfusion, 0.43 ± 0.002; F(1,96) = 157.5, p < 0.001 for WT and KO differences; F(1,96) = 53.1; p < 0.001 for effect of apamin on CVb; two-way ANOVA, t test with Bonferroni correction; Fig. 8F). CVs500 also showed a partial rescue with apamin perfusion (CVs500_WT pre-apamin perfusion, 0.15 ± 0.02; CVs500_WT post-apamin perfusion, 0.16 ± 0.02; CVs500_KO pre-apamin perfusion, 0.35 ± 0.06; CVs500_KO post-apamin perfusion, 0.21 ± 0.04; F(1,49) = 9.63; p = 0.003 for WT and KO differences; F(1,49) = 2.87, p = 0.09 for effect of apamin; two-way ANOVA, t test with Bonferroni correction; Fig. 8G). However, the jitter slope analysis was not able to recapitulate the rescue effects (jitter slope for WT cells pre-apamin perfusion, 0.04 ± 0.007; jitter slope for WT cells post-apamin perfusion, 0.04 ± 0.006; jitter slope for KO cells pre-apamin perfusion, 0.09 ± 0.03; jitter slope for KO cells post-apamin perfusion, 0.04 ± 0.007; F(1,41) = 2.1; p = 0.15 for effect of apamin; F(1,41) = 3.6; p = 0.06 for WT and KO differences; two-way ANOVA, t test with Bonferroni correction (Fig. 8H). Overall, we find that the blockage of SK currents by apamin partially rescues the variability phenotypes, which is consistent with our interpretation that the elevation in SK currents has a role in variability.
Blocking calcium has no rescue affect
From the above observations, blocking SK currents led to only a partial rescue in the phenotype of increased imprecision and unreliability for spikes in KO cells. We therefore asked whether Ca2+, as an upstream modulator of SK, could also contribute to a partial rescue. To do this, we used BAPTA to chelate intracellular Ca2+. We added BAPTA (10 mm) into the intracellular recording pipette and assessed spike precision 10 min after patching. We found that chelating calcium using BAPTA led to a further increase in spike variability instead of rescue. Spike reliability within cell was measured (CVw) in these experiments, and it was found that spikes became significantly more variable post-BAPTA administration than pre-BAPTA for both WT and FXS cells (CVw for WT cells without BAPTA administration, 0.038 ± 0.005, n = 30 cells; CVw for WT cells with BAPTA administration, 0.12 ± 0.02, n = 12 cells; p < 0.001; CVw for KO cells without BAPTA administration, 0.057 ± 0.006, n = 29 cells; CVw for KO cells with BAPTA administration, 0.16 ± 0.04; F(1,81) = 40.5, p < 0.001 for the effect of BAPTA; two-way ANOVA, t test with Bonferroni correction; Fig. 8I). Since Ca2+ has a large number of ion channels as other targets, we concluded that it was not a good target for the rescue of the variability phenotype.
Discussion
We have used a mouse model of autism, the fragile X knock-out mouse, to investigate cellular correlates of autism in hippocampal CA1 pyramidal neurons. We find that the reliability of spiking is impaired in FXS adult animals. We observed an increase in spiking variability both within a cell over multiple trials and between cells in matching trials, in FXS animals. The variability appeared between 6 and 8 weeks of age, consistent with current understanding of autism as a neurodevelopmental disorder. We observed that the increased spiking variability was accompanied by reduced spiking and elevated mAHP currents. We further dissected out the channel contributions to the elevated mAHP and found that SK currents were elevated in adult FXS neurons. However, levels of SK channels were not different between WT and KO mice. To test the hypothesis that elevated SK currents are responsible for unreliable and imprecise spiking in FXS neurons, we used the specific SK blocker apamin and obtained a partial rescue of the phenotype of increased spike variability. In addition to this, we discovered a differential effect of FMRP KO in different hippocampal regions. In agreement with previous published results, we found that there is an increase in excitability in CA3 region of hippocampus in juvenile animals. At the same age, CA1 cells do not show any difference in excitability between WT and KO animals. As the animal becomes older, we observe that the increased excitability seen in CA3 goes away, and for CA1 cells in this age group there is an observation of reduced excitability. This is the first study that has shown that FMRP KO can lead to very different outcomes in different regions of the same brain structure, further highlighting the complexity of the disorder (Figs. 9, 10).
Lowered spike reliability and spike imprecision at cellular intrinsic level in FXS KO hippocampal cells
Our study demonstrates that there is a cell-intrinsic contribution to spiking variability in a mouse model of autism and implicates increased currents from SK channels as part of the mechanism. This finding takes significance with the failure of mGluR antagonists and blockers to significantly alleviate FXS symptoms in human patients. These findings point to the likelihood that the mGluR theory (Erickson et al., 2017) may be incomplete. Specifically, our finding of within-cell variability and between-cell variability in FXS cells shows that cells in the mutant animals may form a heterogeneous population of intrinsically variable neurons, the combined effects of which may lead to network-level incoherence.
This variability in neuronal responses over multiple trials may have implications for multiple behaviors. For example, in learning, the stimulus typically needs to be repeated over multiple trials. If the response of a cell at its intrinsic level as well as between cells is variable, the precision of input spiking as well as output activity would be reduced. This would adversely impact learning rules such as STDP and associativity (Markram et al., 1997; Bi and Poo, 1998). Further, elevated SK currents will have a direct impact on the recorded spontaneous EPSPs at the soma (Ngo-Anh et al., 2005).
SK channel currents elevation may underlie previously observed age-dependent changes in KO cellular physiology
FXS has many clinical and experimental attributes of a neurodevelopmental syndrome (Nimchinsky et al., 2001; Cruz-Martín et al., 2010; Meredith et al., 2012). An age-dependent decline in learning has been shown in the Drosophila fragile X model (Choi et al., 2010). Theta burst stimulation (TBS) LTP has been showed to produced exaggerated LTP in anterior piriform cortex for older KO mice (12–18 months of age) but not in younger mice (<6 months of age; Larson, 2005). Similar results were seen in prefrontal cortex where TBS produced a significant impairment in LTP for 12-month-old mice but not in 2-month-old animals (Martin et al., 2016). Even in humans, it has been shown that the males who have a premutation (have CGG repeats <200) predominantly demonstrated fragile X-associated tremor/ataxia syndrome among patients who were >50 years old (Cornish et al., 2008). At the molecular level, an age-dependent change has also been shown in protein levels in FXS from postnatal day 17 (P17) mice to adults (Tang et al., 2015).
We observed a similar age-dependent change in the spike variability phenotype, which was absent in young KO animals, but manifested in older (>6 weeks old) mice (Fig. 3D). Age-dependent changes in SK channels have been previously reported such that the levels of SK channels are elevated in older animals compared with younger ones (Blank et al., 2003). This complements our finding that the AHP difference between WT and KO animals arose in older animals (Fig. 5) and that the SK channel underlay the difference in the older animals. There is also evidence of age-related changes in Ca2+ levels and Ca2+ uptake mechanisms, such that basal Ca2+ levels are elevated in aging brains (Das and Ghosh, 1996; Kirischuk et al., 1996; Foster and Norris, 1997). Such altered Ca2+ levels would also have a direct impact on the functionality of SK channels. We were not able to directly investigate this possible upstream effect due to the very diverse outcomes of modulating calcium signaling (Fig. 8I).
SK channels may be an FMRP target with multiple physiologic and behavioral outcomes
FMRP targets multiple downstream pathways. It not only is a translational regulator but its N terminus is known to interact with multiple other proteins affecting their functioning (Brown et al., 2010; Deng et al., 2013). Thus, mutations in FMRP affect multiple aspects of cellular functioning. Previous RNA sequencing results have not shown that SK channels are elevated or affected in FXS (Brown et al., 2001; Darnell et al., 2011). However, the functioning of SK currents might be subject to other indirect effects, such as changes in calcium channel distribution, or in their functioning (Meredith et al., 2007; Ferron et al., 2014). Alterations in SK channels have substantial implications for multiple other phenotypes at the cellular and behavioral level. Previous studies have found that FXS KOs have impaired LTP in CA1 hippocampus (Lauterborn et al., 2007; Lee et al., 2011b; Tian et al., 2017). We suggest that this impairment might be partially due to elevated SK currents, which are known to shunt EPSPs and also to increase NMDA Mg2+ block (Ngo-Anh et al., 2005). Hence, our observation of elevated SK currents is consistent with the reported increase in LTP threshold in FXS (Lee et al., 2011b). These effects on LTP may have knock-on effects on memory phenotypes in FXS models. For example, a mild impairment in learning due to a change in platform location in the Morris water maze test has been found in FXS models (D'Hooge et al., 1997; Dobkin et al., 2000). SK currents have been shown to be reduced after learning, and this reduction leads to decreased variability in spiking events (Sourdet et al., 2003).
Comparison with previously published studies
Multiple previous studies have found hyperexcitability in hippocampus in FXS model mice (Deng and Klyachko, 2016; Luque et al., 2017). We do not find a strong significant difference between KO and WT f–I curves (Fig. 3I–L). However, the previously published studies were conducted in the FVB strain of mice and the present study was performed in BL/6 mice. Strain-dependent differences in FXS have previously been found (Paradee et al., 1999; Dobkin et al., 2000; Lee et al., 2011b; Routh et al., 2013). However, Brager et al. (2012) found a lowering of input resistance for KO CA1 cells indicative of reduced excitability. Telias et al. (2015) have found a phenotype similar to ours in human hippocampal culture cells. In their study, they found that FXS cells were incapable of producing multiple spikes on sustained depolarization at multiple holding potentials, similar to what we observe in this study. It is interesting to note that most excitability analyses and f–I curves use 500 ms current steps, whereas in our study a difference arises after 600 ms and leads to a reduction in excitability in the KO (Figs. 1B, 3G). This suggests that a nuanced interpretation of such excitability studies is desirable, as the responses may have complex temporal features not reducible to a single measure of excitability. Another important point of difference is the age of the animals. Previous published studies have been performed on relatively young animals (<P30) compared with the present study, where adult animals (6–8 weeks of age) have also been used for the experiments. Our own results (Fig. 4B,C) show that there is a substantial change in spiking variability as well as AHPs as the animal matures. Other studies have also reported an age-dependent effect of the expression of plasticity phenotypes in FXS mice (Larson, 2005; Martin et al., 2016). Fragile X syndrome is known to be a developmental disease where symptoms and many other biochemical and other pathways are known to undergo changes with the age of the animal. Hence, it may be important to factor in age dependence and possibly use a wider range of ages to re-examine some of the conclusions that have been made from FXS models.
Our study is also the first to examine the effect of variability and spiking asynchrony, which again we find to be an age-dependent effect. We postulate that variability in the firing of individual cells, which previous studies have not explicitly examined, may be a further complicating factor in analyzing baseline changes arising due to the FMRP mutation.
FMRP KO has different physiologic outcomes in different regions of the hippocampus brain regions
In the present study, we make the important and unexpected finding that FMRP KO has differential cell-physiologic effects on different subsections of the same brain region. In an attempt to replicate previous results, we found that there is hyperexcitability in pyramidal neurons of the CA3 region in juvenile mice, (as previously reported by Deng et al., 2019) but there was no significant difference in spiking in the CA1 region at the same age. Further, the phenotype reverses for CA1 cells in adult FMRP KO animals such that pyramidal neurons have reduced excitability, whereas for CA3 regions the hyperexcitability goes away. We suggest that these differential effects arise from the balance of effects among the direct effects of FMRP on SK channels, the indirect effects of FMRP on SK via calcium, and the effects of age, again potentially via calcium.
Direct effects of FMRP on SK
Deng et al. (2019) have shown that the absence of FMRP leads to reduced activity of SK channels. They have shown that a physical interaction between FMRP and SK leads to the activation of SK channels. This suggests that FMRP is a direct activator of SK channels. In the case of mutant animals where FMRP is absent, there is a reduction in SK currents leading to hyperexcitability, which has been observed both in the present and previous study (Fig. 10). We see this effect in CA3 cells for KO animals in the younger age group.
Indirect effects on SK via calcium
SK channels are dependent on intracellular calcium for their activation. Calcium channels are known targets of FMRP for transcription regulators (Brown et al., 2001; Darnell et al., 2011). Calcium influx via L-type Ca2+ channels are known to be elevated in neural progenitor cells (Danesi et al., 2018). Cortical culture neurons have been shown to have significantly elevated basal calcium levels in FXS neurons (Castagnola et al., 2018), and there is increased calcium influx in CA3 cells (Deng et al., 2011, 2013). In reference to human patients, multiple point mutations in Ca2+ genes are associated with autism (Contractor et al., 2015; Pinggera et al., 2015; Nguyen et al., 2018). There are multiple other studies in other neuron types that show increased intracellular calcium in autism models (Laumonnier et al., 2006; Krey and Dolmetsch, 2007; Tessier and Broadie, 2011; Pinggera et al., 2015; Limpitikul et al., 2016). These studies are in agreement with our observations in the present study where we see increased SK currents in CA1 cells, which is expected if the intracellular calcium levels are elevated in FXS neurons (Fig. 10). This increased calcium influx can explain the reduced excitability phenotypes seen in older CA3 and CA1 cells.
Effects of age on calcium
Calcium binding proteins and intracellular calcium dynamics are known to undergo development-dependent alterations in FXS neurons (Tessier and Broadie, 2011). Calcium levels and Ca2+ uptake mechanisms are also known to undergo changes with aging in animals (Das and Ghosh, 1996; Kirischuk et al., 1996; Foster and Norris, 1997). Though there are limited studies available in reference to a change in calcium channel levels with age in autism models, we hypothesize that a change in calcium dynamics is consistent with the fact that autism is a neurodevelopment disorder (Fig. 10). We see this effect in both CA3 and CA1 regions for KO animals, where in CA3 the increased excitability seen in younger animals levels off as the animal becomes adult. On the other hand, because of the increase in calcium, CA1 cells in which excitability was balanced in the young age group goes in the direction of hypoexcitability in older age groups.
Based on these findings, we propose a schematic model to explain these differences (Fig. 10). In the young age group for KO CA3 cells, the absence of FMRP leads to hyperexcitability, as explained above. As the animal becomes older, this effect goes away due to increased calcium influx with age, leading to SK activation and balancing out hyperexcitability. On the other hand, for young KO CA1 cells, both WT and KO start with maximum increased excitability where it is not significantly different. As the animal becomes older, increased calcium influx moves KO cells toward hypoexcitability.
Defects in single-cell physiology may lead to network and population asynchrony
Multiple studies have shown that variability or lack of coordination in spiking exists in FXS at the network level (Arbab et al., 2018b; Talbot et al., 2018). In these tetrode-based studies, variability has been shown at the level of cellular firing patterns in the in vivo network; hence, it is difficult to separate cell-intrinsic properties from synaptic and network phenomena. In our study, we use the CVb parameter to gain an understanding of the correlation between cells at a population level. Using this analysis, we have tried to estimate the variability that exists between individual neuron responses when actually neurons were recorded one at a time. In contrast to tetrode recordings or field response recordings, where the recorded response is a result of multiple neurons and hence is a population response, CVb remains an intrinsic response of neurons depending on its active channels. CVb is an indirect readout of network correlations and should not be inferred as giving actual network correlation information. Clear mechanisms are lacking for this reported asynchrony. However, prior researchers have postulated that synaptic mechanisms may be responsible for the observed asynchrony (Talbot et al., 2018). Our study provides an alternative hypothesis for the mechanism, namely, that the effect is cell intrinsic and is in part due to altered kinetics of SK channels. Other studies provide evidence in favor of a role for SK channels in mediating precision and synchronous firing. For example, Deister et al. (2009) showed that SK channels were responsible for mediating precision in autonomous firing in globus pallidus neurons. Schultheiss et al. (2010) showed that dendritic SK channels also control synchronization. Another study by Kleiman-Weiner et al. (2009) has shown that SK channels play a role in mediating thalamocortical oscillations. A study by Combe et al. (2018) shows that SK channels contribute to phase locking for the spiking of CA1 cells of the hippocampus to slow gamma frequency inputs from Schaffer collaterals. Thus, the alterations in SK channel kinetics in FXS mutants as established by our study might have substantial impacts on phase coding and synchronization of hippocampal cell populations.
A theory of FXS pathophysiology
An alternate theory to mGluR theory in explanation of FXS pathophysiology is the discoordination theory. The idea behind this theory is that many of the cognitive impairments seen in FXS are due to altered coordination in the spiking of cells in networks, which might be due either to an imbalance in excitation-inhibition in the networks (El Idrissi et al., 2005; Heulens et al., 2012) or to an intrinsic incoherent tendency of the neurons (Brock et al., 2002; Fenton, 2016; Radwan et al., 2016). Though this theory is preliminary, it may offer a promising alternative point of view from the traditional mGluR theory. Some studies have shown that intrinsic changes in cells during FXS affect their excitatory/inhibitory balance (Zhang et al., 2014; Lee et al., 2017) and that modulation of these channels leads to the rescue of the FXS phenotype. This present study has probed the dysfunction of a specific ion channel, SK, which was shown to affect the spiking variability at the single-cell level. It is an addition to the increasing list of ion channels that are known to be affected in FXS and may have an important effect on FXS etiology.
Acknowledgments
Acknowledgments: We thank Suranjana Pal, Bishal Basak, Vidya Ramesh, Sarfaraz Nawaz for help with Immunofluorescence experiments. We also thank the National Centre for Biological Sciences Animal Facility and Central Imaging and Flow Cytrometry (CIFF) facility for assisting in animal handling and imaging, respectively. In addition, we thank Giselle Fernandes, Shreya Das Sharma, Aviral Goel, Sathyaa Subramaniyam, Aanchal Bhatia, and Sahil Moza for useful suggestions.
Footnotes
The authors declare no competing financial interests.
D.D. was funded by National Centre for Biological Sciences/Tata Institute of Fundamental Research Intramural funds during the project. National Mouse Resource and Department of Biotechnology provided funds for mouse maintenance and experimental equipment.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.