ORIGINAL RESEARCH article
Analysis of Age-Dependent Alterations in Excitability Properties of CA1 Pyramidal Neurons in an APPPS1 Model of Alzheimer’s Disease
- 1Institute of Biophysics, National Research Council, Palermo, Italy
- 2Université Côte d’Azur, CNRS, IPMC, Valbonne, France
- 3Biomedical Center, Ludwig Maximilian University of Munich, Munich, Germany
Age-dependent accumulation of amyloid-β, provoking increasing brain amyloidopathy, triggers abnormal patterns of neuron activity and circuit synchronization in Alzheimer’s disease (AD) as observed in human AD patients and AD mouse models. Recent studies on AD mouse models, mimicking this age-dependent amyloidopathy, identified alterations in CA1 neuron excitability. However, these models generally also overexpress mutated amyloid precursor protein (APP) and presenilin 1 (PS1) and there is a lack of a clear correlation of neuronal excitability alterations with progressive amyloidopathy. The active development of computational models of AD points out the need of collecting such experimental data to build a reliable disease model exhibiting AD-like disease progression. We therefore used the feature extraction tool of the Human Brain Project (HBP) Brain Simulation Platform to systematically analyze the excitability profile of CA1 pyramidal neuron in the APPPS1 mouse model. We identified specific features of neuron excitability that best correlate either with over-expression of mutated APP and PS1 or increasing Aβ amyloidopathy. Notably, we report strong alterations in membrane time constant and action potential width and weak alterations in firing behavior. Also, using a CA1 pyramidal neuron model, we evidence amyloidopathy-dependent alterations in Ih. Finally, cluster analysis of these recordings showed that we could reliably assign a trace to its correct group, opening the door to a more refined, less variable analysis of AD-affected neurons. This inter-disciplinary analysis, bringing together experimentalists and modelers, helps to further unravel the neuronal mechanisms most affected by AD and to build a biologically plausible computational model of the AD brain.
Increasing accumulation of the beta-amyloid peptides (Aβ) in the brain, a pathological phenomenon termed progressive amyloidopathy, is thought to be an important factor underlying the decline of cognitive functions in Alzheimer’s disease (AD) (Selkoe and Hardy, 2016). In human AD patients, this Aβ accumulation is particularly precocious in the hippocampus, a key brain region involved in memory processing, and worsens with aging. For this reason, Aβ-dependent hippocampus dysfunction has been intensively investigated in pre-clinical studies to identify the basis of memory loss in AD. Progressive age-dependent accumulation of Aβ with the formation of amyloid plaques, as observed in AD patients, is well reproduced in a number of transgenic mouse lines (Morrissette et al., 2009; Ashe and Zahs, 2010). To generate these models, the amyloid precursor protein (APP) is generally overexpressed as a transgene alone or together with another transgene expressing presenilin 1 (PS1), a member of the γ-secretase complex responsible for the cleavage of APP and the production of Aβ. To increase the generation of progressive amyloidopathy in these models, both APP and PS1 harbor mutations that are found in familial Alzheimer’s disease. These pre-clinical models exhibit age-dependent hippocampus dysfunction that correlate with hippocampus-dependent memory loss (Gong et al., 2004; Trinchese et al., 2004; Ma et al., 2012; Kartalou et al., 2020). These models are thus pertinent to investigate the impact of AD-like progressive amyloidopathy on hippocampus function.
In these models, most neurophysiological studies focused on the analysis of synapse function in the hippocampus. These studies identified alterations in several forms of excitatory synapse plasticity, such as long-term potentiation and depression (reviewed in Marchetti and Marie, 2011). More recently, other aspects of hippocampus activity are being investigated in the context of AD, in particular disturbances in intrinsic excitability and network synchrony at the basis of brain oscillatory activities, two other facets of hippocampus function crucial for adequate memory processing (Hanslmayr et al., 2019). Changes in network activity have been observed in human AD patients, who also suffer from a higher incidence of seizures (Amatniek et al., 2006; Scarmeas et al., 2009). Correlating with these clinical observations, work in AD mouse models evidences alterations in hippocampus network synchrony (Goutagny et al., 2013; Cayzac et al., 2015). At a single neuron level, neuron excitability has been characterized by measuring either passive membrane properties, single action potential (AP) kinetics or intrinsic excitability properties of hippocampal CA1 pyramidal neurons in a few of these APP or APP/PS1 over-expressing models when they were at a particular stage of progressive amyloidopathy (Brown et al., 2011; Kaczorowski et al., 2011; Verret et al., 2012; Wykes et al., 2012; Kerrigan et al., 2014; Šišková et al., 2014; Tamagnini et al., 2015). Only three of these studies analyzed specific aspects of neuronal excitability profiles at two ages in the same model to correlate progressive Aβ accumulation to alterations in some of these parameters (Brown et al., 2011; Kaczorowski et al., 2011; Šišková et al., 2014). When these studies are compared, no clear consensus emerges regarding alterations in the different parameters of neuron excitability. This could be at least partially due to the fact that each study investigated a different mouse model of amyloidopathy, moreover at different ages, representing different levels of amyloidopathy progression.
A thorough age-dependent investigation of these neuronal excitability parameters in the same model for correlative analysis of phenotypes with progressing Aβ amyloidopathy has yet to be reported. This type of analysis is particularly relevant for several reasons. First, when reporting phenotypes of these mice, authors usually point to Aβ accumulation as the cause of these phenotypes. Yet, there is increasing evidence that both APP and PS1 proteins per se, regulate neuron function even under physiological conditions (Saura et al., 2004; Weyer et al., 2011; Pousinha et al., 2017; Barthet et al., 2018; Lee et al., 2020). Overexpression of these proteins in the mouse models could therefore in itself perturb the different parameters analyzed, but these phenotypes should be observed independently of Aβ accumulation. It is thus of interest to more clearly correlate alterations in neuronal excitability to progression of Aβ amyloidopathy that occurs with aging in these models. This type of correlative analysis cannot definitely assign functional alterations to Aβ accumulation, but it should help us pinpoint phenotypes that are most likely to be due to this neuropathological hallmark. With the active development of computational modeling of AD (Golriz Khatami et al., 2019), it will be necessary to obtain such experimental data to build a reliable disease model exhibiting AD-like disease progression.
To this end, we performed a patch-clamp analysis of the passive membrane properties, single action potential kinetics (AP) and intrinsic excitability profile of CA1 pyramidal neurons of the hippocampus in an APPPS1 mouse model at three different ages representing amyloidopathy free, weak amyloidopathy and strong amyloidopathy stages of the disease and compared them to control wild-type (WT) littermates. This APPPS1 mouse model (Radde et al., 2006) is particularly useful to study the progression of amyloidopathy as these mice exhibit alterations in Aβ CSF levels, Aβ brain load and Tau CSF levels in a temporal sequence and magnitude of Aβ and Tau changes observed in the CSF of patients with sporadic and dominantly inherited AD (Maia et al., 2013). We used the feature extraction tool available on the Brain Simulation Platform of the Human Brain Project to analyze these data and report the results obtained for 14 features linked to either passive membrane properties, single action potential kinetics (AP) or intrinsic excitability profile. We also took advantage of the latest CA1 pyramidal neuron model developed by Migliore et al. (2018) to predict additional alterations in underlying currents that were not directly accessible by feature extraction of the recorded traces. Finally, using these extracted features, we performed a cluster analysis to predict if we could assign traces of the different conditions to their respective groups.
Materials and Methods
APPPS1 mice carrying transgenes expressing human APP bearing the KM670/671NL “Swedish” mutation and human PS1 bearing the L166P mutation under control of the Thy1 promotor were used in this study (Radde et al., 2006)1. The animals were generated and maintained on a C57BL/6J genetic background (Charles River, France). APP/PS1 and WT male littermates were used at three ages: 1 month (3–4 weeks), 3–4 months, and 9–10 months of age. All experiments were done according to policies on the care and use of laboratory animals of the European Committees Council Directive (86/609/EEC). The protocols were approved by the French Research Ministry following evaluation by a specialized ethics committee (APAFIS#6855-20 16091615385487 v5). All efforts were made to minimize animal suffering and reduce the number of animals used. The animals were housed six per cage under controlled laboratory conditions with a 12-h dark light cycle, a temperature of 22 ± 2°C. Animals had free access to standard rodent diet and tap water.
Biochemical Analysis of Aβ Load
Four APPPS1 mice of 1 month of age, 6 APPPS1 mice of 3–4 months of age, 5 APPPS1 mice of 9–10 months of age and 4 WT mice of 9–10 months of age were sacrificed and their hippocampi were dissected out and snap-frozen in liquid nitrogen to process for biochemical analysis of Aβ load. Briefly, DEA (0.2% Diethylamine in NaCl 50 mM, pH 10) and RIPA lysates [Tris-HCl (pH 7.5) 20 mM, NaCl 150 mM, Na2EDTA 1 mM, NP-40 1%, sodium deoxycholate 1%, sodium pyrophosphate 2.5 mM] were prepared from the frozen tissues. The later was centrifuged at 14,000 g (60 min at 4°C) and the remaining pellet was homogenized in 70% formic acid (FA fraction). The FA fraction was neutralized with 20 × 1 M Tris-HCl buffer at pH 9.5 and used for Aβ analysis. For Aβ detection by immunoblotting, proteins were separated on Tris-Tricine (10–20%, Thermo Fisher Scientific, Germany) gels, transferred to nitrocellulose membranes (0.1 μm, GE Healthcare, United States) which were boiled for 5 min in PBS and subsequently incubated with the blocking solution containing 0.2% I-Block (Thermo Fisher Scientific, Germany) and 0.1% Tween 20 (Merck, United States) in PBS for 1 h, followed by overnight incubation with 2 mg/ml 2D8 antibody in the blocking solution. The rat monoclonal 2D8 antibody against Aβ was described before (Shirotani et al., 2007). Antibody detection was performed using the corresponding anti-rat-IgG-HRP conjugated secondary antibody (Santa Cruz) and chemiluminescence detection reagent ECL (Thermo Fisher Scientific, Germany). Aβ contained in FA fractions was quantified by a sandwich immunoassay using the Meso Scale Aβ Triplex 6E10 plates and Discovery SECTOR Imager 2400 as described previously (Page et al., 2008). Samples were measured in triplicates.
Mice were anesthetized using Ketamine 150 mg/kg/Xylazine 10 mg/kg solution before transcardiac perfusion with the oxygenated (95% O2/5% CO2) ice-cold sucrose cutting solution containing (in mM): KCl 2.5, NaH2PO4 1.25, MgSO4 10, CaCl2 0.5, glucose 11, sucrose 234, NaHCO3 26. The brain was immediately removed and the dissected hippocampus was mounted in an agar support and placed in the previous oxygenated ice-cold sucrose cutting solution. Horizontal slices of the dorsal hippocampus with 250 μm were obtained using standard procedure (Pousinha et al., 2019), placed in an aCSF chamber containing (in mM): NaCl 119, KCl 2.5, NaH2PO4 1.25, MgSO4 1.3, CaCl2 2.5, NaHCO3 26, glucose 11 at 37°C for 1 h and then kept at room temperature. All chemicals were from Sigma-Aldrich (Lyon, France).
For recordings, slices were transferred to a recording chamber containing continuous circulating (2 mL/min), oxygenated and warm (32°C) aCSF. Slices were visualized on an upright microscope with IR-DIC illumination and epi-fluorescence (Scientifica Ltd, Uckfield, United Kingdom). Current-clamp experiments were obtained using a Multiclamp 700B (Molecular Devices, Sunnyvale, CA, United States). Signals were collected and stored using a Digidata 1440 A converter and pCLAMP 10.2 software (Molecular Devices, Sunnyvale, CA, United States). Whole cell current clamp recordings were made using 4–6 MΩ fire-polished glass electrodes and a tight seal (>1 GΩ) on the cell body of the selected neuron was obtained. Internal solution contained (in mM): K-D-gluconate 135, NaCl 5, MgCl2 2, HEPES 10, EGTA 0.5, MgATP 2, NaGTP 0.4. The resting membrane potential (Vm) was first measured in the absence of any spontaneous firing, and only cells more negative than –50 mV were considered. The number of abnormally depolarized neurons were not different in the various conditions tested (about 5% in each condition) and most likely reflected a technical issue such as poor seal. The firing frequency was obtained by clamping neurons at Vh = –65 mV and then injecting pulses of increased intensity in steps of 50 pA (from –200 to 400 pA, 400 ms duration).
Electrophysiological features were extracted from individual experimental traces using the Feature extraction tool available on the Brain Simulation Platform of the Human Brain Project2, which exploits the open source Electrophysiological Feature Extraction Library (eFEL3). Based on our judgment and to allow comparison with previous works, we decided to select 14 features. A total of 1200 traces were analyzed, using 14 different current injections. They were organized in 92 sets of recordings, with 14 neurons from WT-1m animals, 12 neurons from WT-4m, 19 neurons from WT-10m, 11 neurons from AD-1m, 15 neurons from AD-4m, and 21 neurons from AD-10m. Recordings were obtained from 10 control mice (3 WT-1m, 3 WT-4m, and 4 WT-10m) and 9 APPPS1 mice (3 AD-1m, 3 AD-4m, and 3 AD-10m). To calculate spike times, action potential threshold was set at –10 mV.
All simulations were carried out using the NEURON simulation environment (v7.7.2, Hines and Carnevale, 1997). For all simulations we used a reconstructed morphology of a mouse CA1 pyramidal neuron (from neuromorpho.org, cell id NMO_60522, Peng et al., 2016), including only passive properties and an Ih taken from a previous model of rat CA1 neurons (Migliore et al., 2018). Temperature was set at 34°C. On apical dendrites, the Ih channel was distributed following a sigmoidal increase as a function of distance from soma (Migliore et al., 2018). Model and simulation files will be uploaded to the ModelDB database4 (accession no. 266848).
The Matlab Classification Learner app was used to perform the classification of the traces. A 10-fold cross validation was performed, training the models on 90% of the data and validating on the remaining part. Splitting during the validation process was done randomly, and after 10 repeats all samples have been left out once. The classifiers that gave the best accuracy in most of the cases were Ensembles of Bagged or Boosted Trees, which combines several weak learners to produce better predictive performance than using a single weak learner. Bagging (Bootstrap Aggregation) is used in order to reduce the variance of a decision tree. Several subsets of data from training sample are created by choosing randomly with replacement. Each collection of subsets was used to train their decision trees. Boosting is another ensemble technique to create a collection of predictors. With this method, learners are formed sequentially with early learners fitting simple models to the data and then analyzing the results for errors. When an input was misclassified, its weight was increased so that a next iteration can more likely classify it correctly. With both methods, a strong learner is formed at the end of the classification process.
Results in all figures are presented as box plots containing median, quartiles, minimum/maximum values and outliers, if present. Statistical tests were performed using built-in functions in Sigmaplot v13.0 or Graphprism packages. In the great majority of cases, data were normally distributed, and a student t-test was used to compare results for Control and AD groups independently for each age and current injection. For samples that did not follow a normal distribution, a Mann–Whitney Rank-sum test was used. No multiple pairwise comparisons were carried out for statistical analyses. Each group of age was tested independently from the others. For ELISA data analysis, the data followed a normal distribution and were analyzed with a two-way ANOVA followed by Tukey’s multiple comparisons. A difference with a p ≤ 0.05 was considered statistically significant. All statistics are presented in Supplementary Tables 1–4.
To identify alterations in neuron excitability that correlate to different stages of Aβ amyloidopathy in the APPPS1 model, we recorded CA1 pyramidal neurons of mice at 1, 3–4, and 9–10 months of age. At these ages, the mice exhibit vastly different Aβ amyloidopathy profiles (Table 1; Radde et al., 2006; Serneels et al., 2009; Gengler et al., 2010; Maia et al., 2013; Sebastian Monasor et al., 2020). At 1 month of age, expression of human APP and PS1 transgenes in neurons is already high (starting at 2 weeks of age), but Aβ40, Aβ42 levels and the Aβ40/Aβ42 ratio are still very low. Amyloid plaques are completely absent and cognitive functions are still normal. This age group (APPPS1 and WT littermates) can be considered Aβ amyloidopathy free and is labeled WT-1m and AD-1m here after. At 3–4 months of age, Aβ accumulation is rising in the brain (both Aβ40 and Aβ42) and the ratio Aβ40/Aβ42 is fivefold that of control mice (Table 1). Amyloid plaques start to appear in the hippocampus, but these mice do not yet display hippocampus-dependent cognitive deficits. This age group (APPPS1 and WT littermates) thus exhibit markers mimicking weak Aβ amyloidopathy and is labeled WT-4m and AD-4m here after. At 9–10 months of age, strong Aβ amyloidopathy is reached in this model (Table 1) with very high brain levels of Aβ40 and Aβ42 and a stagnating Aβ40/Aβ42 ratio value. Strong amyloid plaque load is observed in the hippocampus and these mice are now exhibiting impaired cognitive functions when tested in hippocampus-dependent memory tasks (Radde et al., 2006; Serneels et al., 2009; Kartalou et al., 2020). This age group (APPPS1 and WT littermates) thus exhibit markers mimicking strong Aβ amyloidopathy, reaching levels observed in AD human patients, correlating to cognitive deficits, and is labeled WT-10m and AD-10m here after. To confirm that the chosen age groups displayed increasing Aβ amyloidopathy in the hippocampus, we measured Aβ load by extracting the FA fraction (RIPA insoluble) of these hippocampi. We used this fraction to detect Aβ by immunoblotting with the 2D8 antibody (Figure 1A) and quantified the levels of three main forms of Aβ (Aβ38, Aβ40, and Aβ42) by ELISA (Figure 1B). As expected, the AD-1m did not exhibit any accumulation of Aβ. In AD-4m mice, Aβ load started to rise, especially the Aβ42 form. In AD-10m, there was strong accumulation of Aβ, both of the Aβ40 and Aβ42 forms, which is in accordance with the original observation in this transgenic mouse line (Radde et al., 2006). No Aβ was detected in WT-10m mice by immunoblotting (Figure 1A) nor by ELISA (data not shown).
Figure 1. Analysis of Aβ load. (A) Immunoblot analysis of Aβ levels (FA fraction) in brains of APPPS1 mice (AD-1m, AD-4m, AD-10m) and WT control mice (WT-10m) (three mice of each group). While the increase in hippocampal Aβ is evident from the FA fractions obtained from AD-10m mice compared to AD-4m mice, no signal can be detected in AD-1m mice and WT-10m mice. (B) ELISA quantification of Aβ38/Aβ40/Aβ42 in FA extracts of AD-1m (four mice), AD-4m (six mice), and AD-10m (five mice). Statistical analysis was performed by Two-way ANOVA. Aβ40: ****p ≤ 0.0001; Aβ42: + +p ≤ 0.01; + + + +p ≤ 0.0001. See full statistical analysis in Supplementary Table 1.
For these three age groups, we recorded the excitability profile of hippocampal CA1 pyramidal neurons in response to increasing somatic current injections by whole-cell patch clamp electrophysiology. We performed an analysis of the electrophysiological traces (see section “Materials and Methods”), and in the following subsections we present the results for selected features relative to passive membrane properties, the first action potential in a train, and the firing behavior. Results for features not discussed in the main text are reported in Supplementary Tables 2–4.
Membrane Time Constant Is Differentially Affected in Aβ Amyloidopathy-Free Mice and in Mice With Strong Aβ Amyloidopathy
We first analyzed the passive membrane properties of CA1 pyramidal neurons in each condition by analyzing the input resistance (Rin), Resting membrane potential (RMP, measured in the experiments), and the membrane time constant. These measures allow us to identify if neurons have modified their general response to current changes, as would occur in a simple passive electrical circuit. Alterations in these passive properties can influence the integration of incoming electric signals. Values of these measures are reported in Figure 2 for all groups. We calculated Rin in two different ways, both based on Ohm’s law: (1) using the native resting potential of a neuron and the holding current needed to bring it to –65 mV; (2) using the steady-state hyperpolarization reached during small negative somatic current injections. In both cases, we did not observe significant differences between genotypes in each age group (pRin = 0.133, 0.207, and 0.369 and pRMP = 0.291, 0.516, and 0.06 after t-test, for 1, 3–4, and 9–10 months respectively). In Figure 2A we plot the values obtained using the first method. RMP was also stable within the different age groups (Figure 2B). The membrane time constant was calculated for different negative current injections (–200, –150, –100, and –50 pA) by fitting the first part of the membrane potential, from the start of the current injection up to the peak of the sag, as schematically shown in Figure 2C for two typical traces from WT-1m (Figure 2C, black line) and AD-1m (Figure 2C, orange line) in response to a –50 pA current injection. By calculating the time constant at different current steps, we had the possibility to discover a feature that has been so far neglected in most experiments, namely a modulation of the membrane time constant with the hyperpolarization level. Typical traces obtained for a –50 pA injection are plotted in Figure 2D, and suggest that AD-1m mice exhibited a significantly slower membrane time constant compared to WT-1m littermates (Figure 2D, left). This alteration normalized in AD-4m mice, and became significantly faster in AD-10m mice, with respect to their aged-matched WT littermates (Figure 2D, right). Quantified results are reported in Figure 2E. Together, these data suggest that this membrane property is differently perturbed in young AD mice (without observable Aβ amyloidopathy) and in old mice (with strong Aβ amyloidopathy), with faster membrane time constant correlating with increased Aβ load.
Figure 2. Differences in passive properties. (A) Input resistance of individual neurons for the different conditions. Average values are shown as closed symbols. (B) Native resting potential of individual neurons for the different conditions. Average values are shown as closed symbols. (C) Typical recordings for a –50 pA current injection from a WT-1m (black trace) or AD-1m (orange) mice. The gray line represents the exponential time course best fitting the initial part of the experimental trace and used to calculate τ. The inset shows the equation used for the fit. (D) Typical traces (left, WT-1m and AD-1m; right WT-10m and AD-10m) used for the calculation of the membrane time constant for a –50 pA input current. Gray lines represent the exponential decay used to calculate τ values; experimental traces have been rescaled to have the same minimum value. (E) Box plot of membrane time constant as a function of the current injection. The red line and the * symbol indicate for which current the values for WT and AD groups were significantly different (p ≤ 0.05). See full statistical analysis for all parameters in Supplementary Table 2.
Ih Current Is Affected in APPPS1 Mice With Strong Aβ Amyloidopathy
In Figure 2E one can also observe that, for all genotypes and ages, the time constant becomes approximately 20% faster with stronger hyperpolarization. These data provide a clear indication of the involvement of a shunting, hyperpolarization activated, current in these neurons. Furthermore, under all conditions, we also observed a prominent sag in all hyperpolarizing current injections. A sag is defined as the difference between the minimal voltage reached during a hyperpolarizing current step and the steady state voltage at the end of the stimulus (Figure 3B), and it is a typical characteristic of cells possessing the Ih current. For this reason, we decided to investigate in more detail the sag amplitude, as it can provide an indirect measure of alterations in Ih peak value or kinetics. A box plot of sag values as a function of hyperpolarizing current injection for each age group is shown in Figure 3A. Typical examples of traces exhibiting a sag are plotted in Figures 3B–D. The AD-1m and AD-4m mice did not show any difference in the sag amplitude (Figure 3A, left and middle panels, p > 0.05 in all cases). AD-10m mice (Figure 3A, right panel) exhibited a significant difference for –200 pA (5.6 ± 0.3 in WT-10m vs. 4.5 ± 0.3 in AD-10m, p = 0.034) and at –50 pA (1.7 ± 0.1 in WT-10m vs. 1.4 ± 0.1 in AD-10m, p = 0.013). This difference could be caused by modifications in Ih current.
Figure 3. Differences in Ih current properties. (A) Box plot of sag amplitude as a function of current injection; See full statistical analysis for sag amplitude in Supplementary Table 2. (B–D) Comparison between model (colored) and experimental (gray) traces under Control and AD conditions for 1, 3–4, and 9–10 months, panels (B–D), respectively; (left and middle) best fit of experimental traces for each condition; (right) kinetic curve for each genotype. Dotted lines represent steady state and solid lines represent time constant of activation of Ih current. *p ≤ 0.05.
Obtaining direct measures of Ih kinetic properties would have been technically challenging in the context of the experiments carried out in this work. For this reason, in order to figure out which Ih kinetic property might be affected in AD-10m neurons, we decided instead to implement a computational model using a three-dimensional morphological reconstruction of a mouse CA1 pyramidal neuron, including only passive properties and Ih. We started from an Ih model file developed for rats (Migliore et al., 2018, modelDB acc.n. 244688 and HBP Live Papers5) and previously validated against a number of experimental findings (Migliore and Migliore, 2012). The kinetic parameters were optimized for mice under control conditions. As reference experimental traces, we used typical recordings at –150, –100, and –50 pA. The standard Run Fitter tool available in NEURON was used to simultaneously optimize the membrane resistance (Rm), the leak reversal potential (e_pas), the Ih peak conductance and the kinetic parameters modeling the activation curve (the shape factor, kl and Vl1/2,) and its time constant (Vt1/2, a0t, zetat, and gmt). The optimization was carried out to find the best fit of a typical set of experimental traces for each condition. As can be seen in the left and middle panels of Figure 3B (for 1 m animals), Figure 3C (4 m animals), and Figure 3D (10 m animals), the optimizations converged into a good solution, with an error lower than 0.3 mV. The final parameters are reported in Table 2, and the Ih kinetics for each case is plotted in the right panels of Figures 3B–D. The role of Ih in modulating subthreshold membrane potential is demonstrated in Supplementary Figure 1, where we show the time constant for models in which we simulated the application of an Ih blocker.
To try to explain the alterations in the sag amplitude by alterations in Ih kinetics, we evaluated differences between WT-10m and AD-10m cells, in the physiological voltage range above –90 mV (Figure 3D). The model (see Table 2) for the AD-10m traces, suggested an approximately 50% increase in the peak Ih channel conductance and an activation kinetics (kl) 50% smoother than in WT-10m. These differences will generate opposite effects in the overall Ih current, and explain why for intermediate current injections (–100 and –150 pA) there was no difference in the average sag amplitude. Furthermore, for small hyperpolarizing currents (i.e., approximately above –80 mV), the activation time constant was much faster under AD conditions. Taken together, the model suggests that in AD-10m mice the amyloidopathy correlates with an increase of the Ih peak conductance, and a faster channel activation. Of note, the model suggested no alterations in Ih parameters for AD-1m mice, correlating with normal sag amplitudes at this age. However, in AD-4m mice the model suggested an increase in the peak Ih channel conductance and a smoother Ih activation kinetics (Table 2), although with a sag amplitude similar to control (Figure 3A). This latter effect could be explained with transitory compensatory mechanisms that are lost with aging and amyloidopathy progression.
Action Potential Width Is Differently Affected in the Different Age Groups of APPPS1 Mice
Modifications of single AP kinetics could strongly influence information transfer in the brain. We thus next focused our attention on the shape of the first AP generated by positive current injections (50–400 pA). Out of the different features analyzed for this AP (see Supplementary Table 3 for full analysis data), AP width of APPPS1 mice showed significant differences at all ages, and we noticed a reversal of this alteration correlating with progression of Aβ amyloidopathy (Figures 4A–C). Indeed, while AD-1m mice exhibited significantly wider AP width, this phenotype was opposite in AD-4m and AD-10m mice, where AP width was significantly narrower than the respective WT control littermates.
Figure 4. Differences in the first Action Potential of train. (A) Box plot of first AP width in the different conditions as a function of current injection. The red line and symbol in each plot indicate the currents for which the average values were significantly different. *p ≤ 0.05. (B) Comparison of typical first APs under different conditions for a 300 pA input current. (C) Phase plots for the first AP for a 300 pA current injection in the different conditions, calculated from all traces at 300 pA. See full statistical analysis for single AP analysis in Supplementary Table 3.
Additional analysis of AP shape was performed to investigate if modifications in Na+ channels function could be evidenced and responsible for this modification in AP width. For this purpose, we constructed phase-plane plots of dV/dt vs. V for the first AP in a train generated during 300 pA current steps, and analyzed the corresponding average dV/dt values (Figures 4B,C). No difference was observed in the different groups (p = 0.548, 0.710, 0.955 after t-test for 1, 3–4, and 9–10 months respectively), suggesting that sodium channel function is not altered at any of these ages in this mouse model (Supplementary Table 3). Our electrophysiological recordings did not allow us to analyze potential alterations in the function of potassium channels, which could also modify AP width.
A Few Parameters of Firing Behavior Are Perturbed in APPPS1 Mice
We next analyzed firing behavior of these neurons, i.e., how they respond to increasing depolarizing current with trains of APs. Alterations in the normal firing behavior of neurons will modify information processing in the brain. We found that many electrophysiological features of this firing behavior, were not perturbed in the AD model over the entire range of current and ages tested (see Supplementary Table 4 for results on additional features). We do not discuss these features here further.
We describe here features that were evaluated in other AD models in previous studies (Brown et al., 2011; Kerrigan et al., 2014; Tamagnini et al., 2015), that are the spike count (Figures 5A,B) and the instantaneous frequency (at 300 pA injection step, Figures 6A,B). We did not evidence any major alterations in spike count for most comparisons (Figure 5B). We however observed that spike count was lower at high current steps in AD-1m mice, suggesting lower excitability at this age. This phenotype normalized in AD-4m mice and reversed AD-10m with spike count becoming higher at high current steps in these mice, suggesting higher excitability. As these long high current steps do not reflect natural patterns of activation in vivo, we also analyzed instantaneous frequency at 300 pA for the first 11 spikes for each group. This instantaneous frequency analysis is a feature that is more relevant to in vivo-like situations. The instantaneous frequency was significantly lower than WT littermates in AD-1m mice from pulse 7 to 10, normalized in AD-4m mice and then was significantly lower again in AD-10m mice, but this time for most pulses analyzed. Together these data suggest that AD-1m and AD-10m neurons exhibit different firing behaviors from their WT littermates.
Figure 5. Differences in spike count. (A) Overlap of typical traces for the different conditions. (B) Box plot of number of spikes as a function of somatic current injection. The red symbols in the graphs indicate currents for which the average values were significantly different. *p ≤ 0.05. See full statistical analysis for these features in Supplementary Table 4.
Figure 6. Differences in firing frequency. (A) Enlargement of the traces to illustrate first 11 spikes in typical experimental traces at 300 pA for each age, under Control and AD conditions. (B) Box plot of instantaneous firing frequency calculated from the first 10 inter-spike intervals, for the different conditions. The red line and symbols in the graphs indicate currents for which the average values were significantly different. *p ≤ 0.05. See full statistical analysis for these features in Supplementary Table 4.
Finally, we hypothesized that the variability of the experimental traces under the different conditions may be caused by the differential impact of ongoing amyloidopathy on individual neurons. This could result in a significant number of neurons in APPPS1 mice still exhibiting control properties, while others exhibit AD-linked alterations. From this point of view, it would be useful to have a way to distinguish if a trace recorded in an APPPS1 mice belongs to a cell already affected by the ongoing pathology or not. For this reason, we performed a classification task using the 11 current-dependent features measured in this work. In Table 3, we report the best classification accuracy for each case. We first investigated the possibility to distinguish traces recorded from animal under control or AD; the Ensemble of Bagged Trees gave 75.1% classification accuracy for traces at 1 m, 73.3% at 4 m and 75.1% at 10 m. Next, we considered the possibility to identify traces affected by AD, after lumping together traces from all ages. We found that the classification accuracy was lower (69.5%), but still above chance level. This was somewhat expected since AD-1m were supposed to be very similar to control, and can thus reduce the likelihood to be distinguished from control. We also tested if it was possible to predict the age of an animal from a single trace. After training and validating on all the traces from control, using an Ensemble Boosted Tree classifier, we obtained an accuracy of 82.4%, which raised to 84.7% in predicting the age using all traces from the AD groups.
Taken together, these results suggest that it is possible to successfully classify a trace as belonging to an animal affected by AD, with an accuracy well above chance. Interestingly, the classification can also reliably predict the age of the animal from a single trace, suggesting that age itself is also an important determinant of the electrophysiological features exhibited by mouse hippocampal pyramidal neurons, both under control and AD conditions.
In this study, we took advantage of the new feature extraction tool of the HBP platform to systematically analyze 14 features of whole-cell patch clamp recordings of CA1 pyramidal neurons at three ages of APPPS1 mice, representing three stages of Aβ amyloidopathy, and compared these features to their respective WT control. One could also extract age-dependent alterations in these parameters using the WT recordings of our analysis. Our study was, however, not specifically designed for this type of evaluation. Indeed, recordings were interleaved between WT and APPPS1 mice within the same age group, usually on the same day of recording, allowing for accurate comparison of genotype effect per age group. To analyze age-dependent alterations in WT mice, we would have designed the study differently by interleaving recordings of different WT age groups within the same days. For this reason, we preferred not to analyze or discuss potential age-dependent alterations in neuron excitability measured in WT mice. Aging-related alterations in neuron excitability of the hippocampus have been investigated previously as reviewed in Oh et al. (2016). We report here strong alterations in membrane time constant and AP width, and weak alterations in Ih current and spiking behavior that correlate with different levels of amyloidopathy in this APPPS1 mouse model. Of note, the membrane time constant, the AP width and spike count exhibited an interesting phenotypic pattern in APPPS1 mice in that they reversed with increasing amyloidopathy. Membrane time constant and AP width were increased in amyloidopathy-free (1 months old) APPS1 mice, but then decreased in mice with weak and strong amyloidopathy. In the opposite way, spike count was decreased in amyloidopathy-free mice, but then increased in APPPS1 mice with strong amyloidopathy. Potential explanations for these findings will be discussed in detail below. We took advantage of these extracted features to perform a cluster analysis and we reliably (69–85%) identified traces to their respective group.
The time constant of the membrane is traditionally associated with the passive properties of a neuron, but its value can also be modulated by any active current that is significantly open around resting potential. Modulation of the membrane time constant, as a function of membrane voltage, can significantly alter the synaptic temporal integration window and the balance between excitation and inhibition, as a function of the synaptic activation frequency. This could at least be partly due to modulation of the hyperpolarization-activated cation current (Ih), that contributes to the sag amplitude, and which is expressed at high density in CA1 pyramidal neurons playing a major role in modulating subthreshold signals (Magee, 1998). A progressive change in the Ih properties, in response to an increase in Aβ load, could explain the observed switch of the membrane time constant as mice get older. From this point of view, the Ih can play a predominant role. In our case, the model suggested that, keeping constant the mice age, the passive properties (Rm and e_pas) did not change much under the control and AD conditions. However, the model supports the notion that amyloidopathy in 4-month animals seems to be a crucial step in the modification of neuronal activity, with some compensatory mechanism in effect. Indeed, while in control neurons, independently from the age and in AD-1m mice, the Ih peak conductance remains constant, its value increases by 50% (increasing the current) with amyloidopathy in AD-4m and AD-10m (Table 2) but, at the same time, the activation is modified in such a way to reduce the overall current.
The main difference observed for AP waveform analysis is the alterations in AP width. The narrowing of this width (∼10–15%) with progression of amyloidopathy is one of the unique features that has been consistently observed in CA1 pyramidal neurons in other AD mouse models at ages representing weak and strong amyloidopathy, such as in the PDAPP mice (Kerrigan et al., 2014; Tamagnini et al., 2015), in the PSAPP mice (Brown et al., 2011) and in the CRND8 mice (Wykes et al., 2012). This narrower AP width could stem from faster sodium channel inactivation, alterations in fast gating potassium channels with roles in AP repolarization, or loss of a contribution from voltage-gated calcium channels to AP waveform. Wykes et al. (2012) proposed the increased expression in the potassium channel Kv3.1b as a potential mechanism in the CRND8 mouse model. There is also accumulating evidence that in vitro and in vivo exposure of hippocampal neurons to Aβ modifies voltage-gated sodium, potassium and calcium channels (Pannaccione et al., 2007; Durán-González et al., 2013; Mayordomo-Cava et al., 2015; Gavello et al., 2018; Ciccone et al., 2019; Wang et al., 2019). We can thus reasonably suggest that this consistent phenotype observed in the different AD mouse models is likely to be due to Aβ accumulation and relevant to the human AD pathology.
Surprisingly, when we analyzed the APPPS1 mice at a very young age when no detectable increase in Aβ levels in the brain were reported (Radde et al., 2006; Maia et al., 2013), we observed wider AP width, a complete opposite phenotype to the one observed with weak and strong amyloidopathy. This very early phenotype was not previously observed in other studies, but recording of other mouse models at this very young amyloidopathy-free age has not yet been reported. This phenotype could be due to the early effect of mutated APP and PS1 transgene expressions per se. As mentioned above, voltage-gated sodium, potassium, and calcium channels govern AP width. A first analysis of voltage-gated sodium channel kinetics with the phase-plots did not show alterations in these channels at any ages in these mice. Our recordings did not allow us to evaluate potential alterations in voltage-gated potassium or calcium channels. However, there is evidence that overexpression of APP per se increases sodium channel currents and cell surface expression (Liu et al., 2015; Li et al., 2016). Also, the voltage-gated sodium channel beta2-subunit, an auxiliary subunit of the voltage-gated sodium channel involved both in channel gating and cell surface expression of α subunits of the sodium channel, was shown to be a substrate for PS1-dependent gamma-secretase-mediated cleavage (Kim et al., 2005). Furthermore, overexpression of a mutated PS1 (albeit the deltaE9 mutation) decreased potassium channel current and cell surface expression (Plant et al., 2002). It is therefore reasonable to speculate that overexpression of mutated APP and PS1 levels could significantly alter AP kinetics by acting on sodium and potassium channels. Independent of Aβ accumulation, these types of alterations could contribute to the widening of the AP width observed in very young APPPS1 mice. With age, the brain of these mice might be slowly adapting to cellular alterations provoked by this overexpression within the first months. With progressive amyloidopathy, AP width then reverses phenotype, becoming abnormally narrower as observed in the other AD mouse models with ongoing amyloidopathy. Together these data suggest that AP width alterations are closely linked to APP processing and the underlying molecular mechanisms should be further investigated to identify this relationship.
Spike count was decreased in young amyloidopathy-free mice, suggesting hypoexcitability, and increased in old APPPS1 mice with strong amyloidopathy, suggesting hyperexcitability. These phenotypes were however only observed at high current steps (350 and 400 pA) and not at lower current steps (50–300 pA), which are likely to represent more physiological conditions. Regarding these findings in strong amyloidopathy conditions, our data confirm other studies using a similar protocol with the lack of phenotype at weak current injections (100–300 pA) as observed in the PDAPP model (Kerrigan et al., 2014; Tamagnini et al., 2015). They however do not support a study in PSAPP mice where hypoexcitability was observed (Brown et al., 2011). Yet, all three previous studies identified higher instantaneous frequency with weak and strong amyloidopathy in these other two models (Brown et al., 2011; Kerrigan et al., 2014; Tamagnini et al., 2015). Here, we report lower instantaneous frequency in 9–10 months old APPPS1 mice. At present, we cannot explain this discrepancy. Nevertheless, altogether these data suggest that CA1 pyramidal neurons exposed to strong amyloidopathy are quite resistant to disruption in firing behavior.
Surprisingly, in our model, both spike count and instantaneous frequency measures transiently normalized in 3–4 months old mice exhibiting weak amyloidopathy. We could speculate that in this APPPS1 model early overexpression of mutated APP and PS1 perturb homeostatic processes that drive neuron excitability and that the brain transiently copes with these alterations restoring the required homeostatic balance. With aging and increasing amyloidopathy, the network cannot cope with the increasing cellular modifications slowly degrading homeostatic processes with increasing perturbations of firing behavior. We currently have not proof of this possibility in this model, but the general concept of loss of firing homeostasis in the context of AD has been well described in a recent review (Styr and Slutsky, 2018).
Previous studies reporting the activity profile of CA1 pyramidal neurons in vivo in models similar to our APPPS1 model, i.e., overexpressing mutated forms of APP and of PS1, support our data correlating high levels of Aβ-amyloidopathy, including Aβ plaques, with increased excitability. Indeed, female and male mice (10–14 months) heterozygous for both transgenes, K670N/M671L-mutated amyloid precursor protein (APPswe), and the M146V mutated presenilin 1 (PS1) under control of a neuron-specific promoter (Willuweit et al., 2009), which display strong Aβ deposition, exhibited increased action potential frequency in CA1 pyramidal neurons in vivo (Šišková et al., 2014). These authors analyzed these mice at 2–4 months, before Aβ plaque appearance, and did not observe alterations in CA1 firing behavior in vivo at this age. However, we do not know the levels of Aβ accumulation in 2–4 months old mice of this AD model, and these could represent early Aβ accumulation with rise in soluble forms of Aβ, similar to our 3–4 months APPPS1 mice where we also did not evidence any alterations in firing behavior. Furthermore, another similar mouse model, the APP23xPS45 double transgenic mouse model overexpressing the human APP with the Swedish (670/671) mutation and the human G384A-mutated presenilin 1 (PS1) (Busche et al., 2008) was analyzed for alterations in spontaneous calcium transients in vivo as a measure of neuron activity in the CA1 region of the hippocampus (Busche et al., 2012). In this study, they report that Aβ accumulation, before plaques and after plaque deposition, correlates with increased hyperactivity of CA1 pyramidal neurons. This study did not analyze AD mice before an increase of soluble Aβ. They also report that an acute application of soluble synthetic AβS26C dimer (100 nM) itself is sufficient to increase CA1 pyramidal neuron hyperactivity in WT mice. Together these data confirm our finding that a strong increase in Aβ load correlates with hyperactivity of CA1 pyramidal neurons.
To summarize our findings, we noted strong age-dependent alterations in membrane time constant and AP width and weak age-dependent alterations in Ih current and firing behavior correlating with levels of Aβ amyloidopathy. It is important to note that we have no evidence that Aβ accumulation is mediating these phenotypes. Unfortunately, using these types of APP and PS1 overexpressing AD models, it will be difficult to sort out what phenotype is directly due to Aβ accumulation, and what phenotype is due to overexpression of mutated APP or PS1 as detailed above, or to accumulation of other cleaved APP peptides. For example, the intracellular C-terminal domain (AICD) of APP is also overproduced in these APP-based AD models (Pousinha et al., 2017) and there is evidence that it can also modulate per se CA1 pyramidal neuron excitability (Pousinha et al., 2019). A throughout evaluation of a new type of AD mouse models exhibiting progressive amyloidopathy, the APP knock-in models (Saito et al., 2014), could be useful to address some of these issues. Indeed, in this model of Aβ amyloidopathy, which progresses at a much slower rate than in APP over-expressing models, the endogenous mouse APP gene was replaced by a humanized form of the APP gene containing familial AD mutations. The human APP gene is thus driven by the endogenous mouse APP promoter avoiding APP overexpression.
One of the hallmarks of AD is the formation of more or less spatially sparse and distributed plaques that progressively hinder neuronal functions. This implies that, at any given time, there can be a significant population of still “normal” neurons in an APPPS1 animal. One limitation of our work is that, in old APPS1 mice, that harbor Aβ plaques, we could not directly relate the location of the patched neurons to the vicinity of these plaques. It is possible that, under this condition, we have patched neurons that were differently affected due to their location with respect to these plaques. Yet, the clustering analysis point out to an interesting possibility. This approach may help in preferentially selecting for analysis only those neurons that can be classified as affected by amyloidopathy, reducing the overall variability in the measured differences between control and AD animals.
We identified specific features of neuron excitability that correlated best with either with over-expression of mutated APP and PS1 (in AD-1m mice) or increasing Aβ amyloidopathy (in AD-4m and AD-10m mice). By improving our knowledge of the ionic mechanisms that are most affected by AD, and by making experimentally testable predictions on the role of Ih, the results can be usefully exploited by both modelers and experimentalists. Modelers can use them for computational modeling studies on the hippocampal network under AD conditions, and experimentalists can use our clustering analysis to improve the quality of their analyses on the effects of amyloidopathy.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: data and models are available at the live paper section of ebrains (https://humanbrainproject.github.io/hbp-bsp-live-papers/2021/vitale_et_al_2021/vitale_et_al_2021.html). Model is also available on Model DB accession no. 266848 (https://senselab.med.yale.edu/ModelDB/showmodel.cshtml?model=266848#tabs-1).
The animal study was reviewed and approved by the French Research Ministry.
AS-P and HM maintained the mouse colony, performed the genotyping, and performed the electrophysiology experiments. PV and RM analyzed the data. CL performed the cluster analysis. MW performed the biochemical analysis of Aβ load. HM and MM designed the study and wrote manuscript with help from all other authors. All authors contributed to the article and approved the submitted version.
Association France Alzheimer to AS-P and HM, the EU Joint Program–Neurodegenerative Disease Research (JPND) Project CIRCPROT (jointly funded by ANR and MIUR and EU Horizon 2020 Grant Agreement 643417) to HM and MM, the EU Horizon 2020 Framework Program for Research and Innovation (Specific Grant 945539, Human Brain Project SGA3) to MM, the Flag ERA JTC 2019 (MILEDI Project) to HM and MM; the UCAJEDI Investments in the Future project managed by the National Research Agency (ANR-15-IDEX-01; ComputaBrain project) to HM, Fenix computing and storage resources under the Specific Grant Agreement No. 800858 (Human Brain Project ICEI), and a grant from the Swiss National Supercomputing Centre (CSCS) under project ID ich011 to MM.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Heike Hampel and Brigitte Nuscher, Biomedical Center, Ludwig Maximilian University of Munich, Germany, for technical help. We gratefully acknowledge the computing time granted by CINECA (Italy) on the MARCONI HPC system.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2021.668948/full#supplementary-material
- ^ https://www.alzforum.org/research-models/appps1
- ^ https://www.humanbrainproject.eu/en/hbp-platforms
- ^ https://github.com/BlueBrain/eFEL
- ^ https://senselab.med.yale.edu/modeldb/
- ^ https://humanbrainproject.github.io/hbp-bsp-live-papers/2018/migliore_et_al_2018/migliore_et_al_2018.html
Amatniek, J. C., Hauser, W. A., Del Castillo-Castaneda, C., Jacobs, D. M., Marder, K., Bell, K., et al. (2006). Incidence and predictors of seizures in patients with Alzheimer’s disease. Epilepsia 47, 867–872. doi: 10.1111/j.1528-1167.2006.00554.x
Barthet, G., Jordà-Siquier, T., Rumi-Masante, J., Bernadou, F., Müller, U., and Mulle, C. (2018). Presenilin-mediated cleavage of APP regulates synaptotagmin-7 and presynaptic plasticity. Nat. Commun. 9:4780. doi: 10.1038/s41467-018-06813-x
Brown, J. T., Chin, J., Leiser, S. C., Pangalos, M. N., and Randall, A. D. (2011). Altered intrinsic neuronal excitability and reduced Na+ currents in a mouse model of Alzheimer’s disease. Neurobiol. Aging 32, 2109–2114. doi: 10.1016/j.neurobiolaging.2011.05.025
Busche, M. A., Chen, X., Henning, H. A., Reichwald, J., Staufenbiel, M., Sakmann, B., et al. (2012). Critical role of soluble amyloid-β for early hippocampal hyperactivity in a mouse model of Alzheimer’s disease. Proc. Natl. Acad. Sci. U S A. 109, 8740–8745. doi: 10.1073/pnas.1206171109
Busche, M. A., Eichhoff, G., Adelsberger, H., Abramowski, D., Wiederhold, K. H., Haass, C., et al. (2008). Clusters of hyperactive neurons near amyloid plaques in a mouse model of Alzheimer’s disease. Science 321, 1686–1689. doi: 10.1126/science.1162844
Cayzac, S., Mons, N., Ginguay, A., Allinquant, B., Jeantet, Y., and Cho, Y. H. (2015). Altered hippocampal information coding and network synchrony in APP-PS1 mice. Neurobiol. Aging 36, 3200–3213. doi: 10.1016/j.neurobiolaging.2015.08.023
Ciccone, R., Franco, C., Piccialli, I., Boscia, F., Casamassa, A., de Rosa, V., et al. (2019). Amyloid β-Induced Upregulation of Nav1.6 Underlies Neuronal Hyperactivity in Tg2576 Alzheimer’s Disease Mouse Model. Sci. Rep. 9:13592. doi: 10.1038/s41598-019-50018-1
Durán-González, J., Michi, E. D., Elorza, B., Perez-Córdova, M. G., Pacheco-Otalora, L. F., Touhami, A., et al. (2013). Amyloid β peptides modify the expression of antioxidant repair enzymes and a potassium channel in the septohippocampal system. Neurobiol. Aging 34, 2071–2076. doi: 10.1016/j.neurobiolaging.2013.02.005
Gavello, D., Calorio, C., Franchino, C., Cesano, F., Carabelli, V., Carbone, E., et al. (2018). Early Alterations of Hippocampal Neuronal Firing Induced by Abeta42. Cereb. Cortex 28, 433–446. doi: 10.1093/cercor/bhw377
Gengler, S., Hamilton, A., and Holscher, C. (2010). Synaptic plasticity in the hippocampus of a APP/PS1 mouse model of Alzheimer’s disease is impaired in old but not young mice. PLoS One 5:e9764. doi: 10.1371/journal.pone.0009764
Golriz Khatami, S., Robinson, C., Birkenbihl, C., Domingo-Fernández, D., Hoyt, C. T., and Hofmann-Apitius, M. (2019). Challenges of Integrative Disease Modeling in Alzheimer’s Disease. Front. Mol. Biosci. 6:158. doi: 10.3389/fmolb.2019.00158
Gong, B., Vitolo, O. V., Trinchese, F., Liu, S., Shelanski, M., and Arancio, O. (2004). Persistent improvement in synaptic and cognitive functions in an Alzheimer mouse model after rolipram treatment. J. Clin. Invest. 114, 1624–1634. doi: 10.1172/JCI22831
Goutagny, R., Gu, N., Cavanagh, C., Jackson, J., Chabot, J. G., Quirion, R., et al. (2013). Alterations in hippocampal network oscillations and theta-gamma coupling arise before Aβ overproduction in a mouse model of Alzheimer’s disease. Eur. J. Neurosci. 37, 1896–1902. doi: 10.1111/ejn.12233
Kaczorowski, C. C., Sametsky, E., Shah, S., Vassar, R., and Disterhoft, J. F. (2011). Mechanisms underlying basal and learning-related intrinsic excitability in a mouse model of Alzheimer’s disease. Neurobiol. Aging 32, 1452–1465. doi: 10.1016/j.neurobiolaging.2009.09.003
Kartalou, G. I., Salgueiro-Pereira, A. R., Endres, T., Lesnikova, A., Casarotto, P., Pousinha, P., et al. (2020). Anti-Inflammatory Treatment with FTY720 Starting after Onset of Symptoms Reverses Synaptic Deficits in an AD Mouse Model. Int. J. Mol. Sci. 21:8957. doi: 10.3390/ijms21238957
Kerrigan, T. L., Brown, J. T., and Randall, A. D. (2014). Characterization of altered intrinsic excitability in hippocampal CA1 pyramidal cells of the Aβ-overproducing PDAPP mouse. Neuropharmacology 79, 515–524. doi: 10.1016/j.neuropharm.2013.09.004
Kim, D. Y., Ingano, L. A. M., Carey, B. W., Pettingell, W. H., and Kovacs, D. M. (2005). Presenilin/gamma-secretase-mediated cleavage of the voltage-gated sodium channel beta2-subunit regulates cell adhesion and migration. J. Biol. Chem. 280, 23251–23261. doi: 10.1074/jbc.M412938200
Lee, S. H., Kang, J., Ho, A., Watanabe, H., Bolshakov, V. Y., and Shen, J. (2020). APP Family Regulates Neuronal Excitability and Synaptic Plasticity but Not Neuronal Survival. Neuron 108, 676.e–690.e. doi: 10.1016/j.neuron.2020.08.011
Li, S., Wang, X., Ma, Q.-H., Yang, W.-L., Zhang, X.-G., Dawe, G. S., et al. (2016). Amyloid precursor protein modulates Nav1.6 sodium channel currents through a Go-coupled JNK pathway. Sci. Rep. 6:39320. doi: 10.1038/srep39320
Liu, C., Tan, F. C. K., Xiao, Z.-C., and Dawe, G. S. (2015). Amyloid precursor protein enhances Nav1.6 sodium channel cell surface expression. J. Biol. Chem. 290, 12048–12057. doi: 10.1074/jbc.M114.617092
Ma, T., Du, X., Pick, J. E., Sui, G., Brownlee, M., and Klann, E. (2012). Glucagon-like peptide-1 cleavage product GLP-1(9-36) amide rescues synaptic plasticity and memory deficits in Alzheimer’s disease model mice. J. Neurosci. 32, 13701–13708. doi: 10.1523/JNEUROSCI.2107-12.2012
Maia, L. F., Kaeser, S. A., Reichwald, J., Hruscha, M., Martus, P., Staufenbiel, M., et al. (2013). Changes in amyloid-β and Tau in the cerebrospinal fluid of transgenic mice overexpressing amyloid precursor protein. Sci. Transl. Med. 5:194re2. doi: 10.1126/scitranslmed.3006446
Mayordomo-Cava, J., Yajeya, J., Navarro-López, J. D., and Jiménez-Díaz, L. (2015). Amyloid-β(25-35) Modulates the Expression of GirK and KCNQ Channel Genes in the Hippocampus. PLoS One 10:e0134385. doi: 10.1371/journal.pone.0134385
Magee, J. C. (1998). Dendritic hyperpolarization-activated currents modify the integrative properties of hippocampal CA1 pyramidal neurons. J. Neurosci. 18, 7613–7624. doi: 10.1523/JNEUROSCI.18-19-07613.1998
Migliore, M., and Migliore, R. (2012). Know your current Ih: interaction with a shunting current explains the puzzling effects of its pharmacological or pathological modulations. PLoS One 7:e36867. doi: 10.1371/journal.pone.0036867
Migliore, R., Lupascu, C. A., Bologna, L. L., Romani, A., Courcol, J.-D., Antonel, S., et al. (2018). The physiological variability of channel density in hippocampal CA1 pyramidal cells and interneurons explored using a unified data-driven modeling workflow. PLoS Comput. Biol. 14:e1006423. doi: 10.1371/journal.pcbi.1006423
Morrissette, D. A., Parachikova, A., Green, K. N., and LaFerla, F. M. (2009). Relevance of transgenic mouse models to human Alzheimer disease. J. Biol. Chem. 284, 6033–6037. doi: 10.1074/jbc.R800030200
Oh, M. M., Simkin, D., and Disterhoft, J. F. (2016). Intrinsic Hippocampal Excitability Changes of Opposite Signs and Different Origins in CA1 and CA3 Pyramidal Neurons Underlie Aging-Related Cognitive Deficits. Front. Syst. Neurosci. 10:00052. doi: 10.3389/fnsys.2016.00052
Page, R. M., Baumann, K., Tomioka, M., Pérez-Revuelta, B. I., Fukumori, A., Jacobsen, H., et al. (2008). Generation of Abeta38 and Abeta42 is independently and differentially affected by familial Alzheimer disease-associated presenilin mutations and gamma-secretase modulation. J. Biol. Chem. 283, 677–683. doi: 10.1074/jbc.M708754200
Pannaccione, A., Boscia, F., Scorziello, A., Adornetto, A., Castaldo, P., Sirabella, R., et al. (2007). Up-regulation and increased activity of KV3.4 channels and their accessory subunit MinK-related peptide 2 induced by amyloid peptide are involved in apoptotic neuronal death. Mol. Pharmacol. 72, 665–673. doi: 10.1124/mol.107.034868
Peng, Y., Lu, Z., Li, G., Piechowicz, M., Anderson, M., Uddin, Y., et al. (2016). The autism-associated MET receptor tyrosine kinase engages early neuronal growth mechanism and controls glutamatergic circuits development in the forebrain. Mol. Psychiatry 21, 925–935. doi: 10.1038/mp.2015.182
Plant, L. D., Boyle, J. P., Thomas, N. M., Hipkins, N. J., Benedikz, E., Hooper, N. M., et al. (2002). Presenilin-1 mutations alter K+ currents in the human neuroblastoma cell line, SH-SY5Y. Neuroreport 13, 1553–1556. doi: 10.1097/00001756-200208270-00013
Pousinha, P. A., Mouska, X., Bianchi, D., Temido-Ferreira, M., Rajão-Saraiva, J., Gomes, R., et al. (2019). The Amyloid Precursor Protein C-Terminal Domain Alters CA1 Neuron Firing, Modifying Hippocampus Oscillations and Impairing Spatial Memory Encoding. Cell. Rep. 29, 317.e–331.e. doi: 10.1016/j.celrep.2019.08.103
Pousinha, P. A., Mouska, X., Raymond, E. F., Gwizdek, C., Dhib, G., Poupon, G., et al. (2017). Physiological and pathophysiological control of synaptic GluN2B-NMDA receptors by the C-terminal domain of amyloid precursor protein. Elife 6:25659. doi: 10.7554/eLife.25659
Radde, R., Bolmont, T., Kaeser, S. A., Coomaraswamy, J., Lindau, D., Stoltze, L., et al. (2006). Abeta42-driven cerebral amyloidosis in transgenic mice reveals early and robust pathology. EMBO Rep. 7, 940–946. doi: 10.1038/sj.embor.7400784
Saura, C. A., Choi, S. Y., Beglopoulos, V., Malkani, S., Zhang, D., Shankaranarayana Rao, B. S., et al. (2004). Loss of presenilin function causes impairments of memory and synaptic plasticity followed by age-dependent neurodegeneration. Neuron. 42, 23–36. doi: 10.1016/s0896-6273(04)00182-5
Scarmeas, N., Honig, L. S., Choi, H., Cantero, J., Brandt, J., Blacker, D., et al. (2009). Seizures in Alzheimer disease: who, when, and how common? Arch. Neurol. 66, 992–997. doi: 10.1001/archneurol.2009.130
Sebastian Monasor, L., Müller, S. A., Colombo, A. V., Tanrioever, G., König, J., Roth, S., et al. (2020). Fibrillar Aβ triggers microglial proteome alterations and dysfunction in Alzheimer mouse models. Elife 9:e54083. doi: 10.7554/eLife.54083
Serneels, L., Van Biervliet, J., Craessaerts, K., Dejaegere, T., Horré, K., Van Houtvin, T., et al. (2009). gamma-Secretase heterogeneity in the Aph1 subunit: relevance for Alzheimer’s disease. Science 324, 639–642. doi: 10.1126/science.1171176
Shirotani, K., Tomioka, M., Kremmer, E., Haass, C., and Steiner, H. (2007). Pathological activity of familial Alzheimer’s disease-associated mutant presenilin can be executed by six different gamma-secretase complexes. Neurobiol. Dis. 27, 102–107. doi: 10.1016/j.nbd.2007.04.011
Šišková, Z., Justus, D., Kaneko, H., Friedrichs, D., Henneberg, N., Beutel, T., et al. (2014). Dendritic structural degeneration is functionally linked to cellular hyperexcitability in a mouse model of Alzheimer’s disease. Neuron 84, 1023–1033. doi: 10.1016/j.neuron.2014.10.024
Tamagnini, F., Novelia, J., Kerrigan, T. L., Brown, J. T., Tsaneva-Atanasova, K., and Randall, A. D. (2015). Altered intrinsic excitability of hippocampal CA1 pyramidal neurons in aged PDAPP mice. Front. Cell. Neurosci. 9:372. doi: 10.3389/fncel.2015.00372
Trinchese, F., Liu, S., Battaglia, F., Walter, S., Mathews, P. M., and Arancio, O. (2004). Progressive age−related development of Alzheimer−like pathology in APP/PS1 mice. Ann Neurol. 55, 801–814. doi: 10.1002/ana.20101
Verret, L., Mann, E. O., Hang, G. B., Barth, A. M. I., Cobos, I., Ho, K., et al. (2012). Inhibitory interneuron deficit links altered network activity and cognitive dysfunction in Alzheimer model. Cell 149, 708–721. doi: 10.1016/j.cell.2012.02.046
Wang, Y.-X., Xia, Z.-H., Jiang, X., Li, L.-X., An, D., Wang, H.-G., et al. (2019). Genistein Inhibits Aβ25-35-Induced Neuronal Death with Changes in the Electrophysiological Properties of Voltage-Gated Sodium and Potassium Channels. Cell. Mol. Neurobiol. 39, 809–822. doi: 10.1007/s10571-019-00680-w
Weyer, S. W., Klevanski, M., Delekate, A., Voikar, V., Aydin, D., Hick, M., et al. (2011). APP and APLP2 are essential at PNS and CNS synapses for transmission, spatial learning and LTP. EMBO J. 30, 2266–2280. doi: 10.1038/emboj.2011.119
Willuweit, A., Velden, J., Godemann, R., Manook, A., Jetzek, F., Tintrup, H., et al. (2009). Early-onset and robust amyloid pathology in a new homozygous mouse model of Alzheimer’s disease. PLoS One. 4:e7931. doi: 10.1371/journal.pone.0007931
Keywords: hippocampus, electrophysiological features, clustering analysis, computational modeling, amyloidopathy
Citation: Vitale P, Salgueiro-Pereira AR, Lupascu CA, Willem M, Migliore R, Migliore M and Marie H (2021) Analysis of Age-Dependent Alterations in Excitability Properties of CA1 Pyramidal Neurons in an APPPS1 Model of Alzheimer’s Disease. Front. Aging Neurosci. 13:668948. doi: 10.3389/fnagi.2021.668948
Received: 17 February 2021; Accepted: 14 May 2021;
Published: 11 June 2021.
Edited by:Stephen D. Ginsberg, Nathan Kline Institute for Psychiatric Research, United States
Reviewed by:Mark J. Wall, University of Warwick, United Kingdom
Ksenia Kastanenka, Massachusetts General Hospital and Harvard Medical School, United States
Krishnan Padmanabhan, University of Rochester, United States
Copyright © 2021 Vitale, Salgueiro-Pereira, Lupascu, Willem, Migliore, Migliore and Marie. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†These authors share first authorship
‡These authors share senior authorship