# Impaired astrocytic Ca2+ signaling in awake-behaving Alzheimer’s disease transgenic mice

1. Knut Sindre Åbjørsbråten 1
2. Gry HE Syverstad Skaaraas 2
3. Céline Cunen 3 4
4. Daniel M Bjørnstad 1
5. Kristin M Gullestad Binder 1
6. Laura Bojarskaite 1 5
7. Vidar Jensen 1
8. Lars NG Nilsson 6
9. Shreyas B Rao 2
10. Wannan Tang 1 7
11. Gudmund Horn Hermansen 3
12. Erlend A Nagelhus 1
13. Ole Petter Ottersen 8
14. Reidun Torp 2
15. Rune Enger rune.enger@medisin.uio.no 1

1. GliaLab at the Letten Centre, Division of Anatomy, Department of Molecular Medicine, Institute of Basic Medical Sciences, University of Oslo
2. Division of Anatomy, Department of Molecular Medicine, Institute of Basic Medical Sciences, University of Oslo
3. Statistics and Data Science group, Department of Mathematics, Faculty of Mathematics and Natural Sciences, University of Oslo
4. Norwegian Computing Center
5. Department of Neurology, Oslo University Hospital
6. Department of Pharmacology, University of Oslo and Oslo University Hospital
7. Department of Clinical and Molecular Medicine, Norwegian University of Science and Technology
8. Office of the President, Karolinska Institutet

2022-07-14

- Research Article

- Neuroscience

- Alzheimer's disease
- calcium imaging
- norepinephrine
- astrocytes
- Mouse

- publisher-id 75055
- doi 10.7554/eLife.75055
- elocation-id e75055

## Abstract

Increased astrocytic Ca2+ signaling has been shown in Alzheimer’s disease mouse models, but to date no reports have characterized behaviorally induced astrocytic Ca2+ signaling in such mice. Here, we employ an event-based algorithm to assess astrocytic Ca2+ signals in the neocortex of awake-behaving tg-ArcSwe mice and non-transgenic wildtype littermates while monitoring pupil responses and behavior. We demonstrate an attenuated astrocytic Ca2+ response to locomotion and an uncoupling of pupil responses and astrocytic Ca2+ signaling in 15-month-old plaque-bearing mice. Using the genetically encoded fluorescent norepinephrine sensor GRABNE, we demonstrate a reduced norepinephrine signaling during spontaneous running and startle responses in the transgenic mice, providing a possible mechanistic underpinning of the observed reduced astrocytic Ca2+ responses. Our data points to a dysfunction in the norepinephrine–astrocyte Ca2+ activity axis, which may account for some of the cognitive deficits observed in Alzheimer’s disease.

## Introduction

Since astrocytic Ca2+ signals were first discovered in the early 1990s, they have been the object of numerous studies exploring their roles in brain physiology and pathophysiology. Importantly, such signals have been shown to occur in response to a wide array of neurotransmitters and trigger the release of substances that affect neuronal signaling and the vasculature. A growing body of evidence suggests that astrocytic Ca2+ signals play important roles in higher brain functions such as memory formation and cortical processing, mediated in part through the neuromodulatory systems 3 Adamsky et al.2018 34 Kol et al.2020 50 Poskanzer and Yuste 2016 49 Poskanzer and Yuste 2011 74 Ye et al.2020 45 Paukert et al.2014.

Astrocytic Ca2+ signaling in an Alzheimer’s disease (AD) mouse model was first described by 35 Kuchibhotla et al.2009, who found pathological Ca2+ waves originating at amyloid plaques, and a general increase in astrocytic Ca2+ signaling. Later, it was demonstrated that plaque-associated astrocytic hyperactivity was mediated through activation of metabotropic purine receptors 16 Delekate et al.2014 52 Reichenbach et al.2018. These studies were performed under anesthesia, which severely attenuates physiological Ca2+ signals 65 Thrane et al.2012.

The field of astrocytic Ca2+ signaling is undergoing a revolution as developments in optical imaging and genetically encoded fluorescent sensors now allow us to monitor these signals in awake-behaving mice, without the confounding effects of anesthesia 61 Srinivasan et al.2015 9 Bojarskaite et al.2020. Such studies have revealed exceedingly rich and complex astrocytic Ca2+ signaling ranging from large activations of nearly all astrocytes in a field of view (FOV) under locomotion and startle responses due to noradrenergic activity 17 Ding et al.2013 45 Paukert et al.2014, to small, localized signals occurring spontaneously or as a response to local neuronal activity 6 Bindocci et al.2017 62 Stobart et al.2018 61 Srinivasan et al.2015. New analytical tools now also enable us to accurately quantify and describe these signals 70 Wang et al.2019 7 Bjørnstad et al.2021.

The brain noradrenergic system is crucial for mediating responses to external environmental stimuli, optimizing central nervous system performance and thus for arousal and cognition. Interestingly, norepinephrine is one of the main drivers of astrocytic Ca2+ signaling activating astrocytic α1-adrenergic receptors 17 Ding et al.2013 65 Thrane et al.2012 45 Paukert et al.2014. Notably, noradrenergic signaling in relation to locomotion and startle responses causes global increases of astrocytic Ca2+ 17 Ding et al.2013 61 Srinivasan et al.2015. The main downstream effects of this astrocytic Ca2+ activity are not yet well understood, but thought to support neurons metabolically 44 O’Donnell et al.2012, or through dynamic changes of the size and composition of the extracellular fluids 69 Wang et al.2012. Nonetheless, astrocytes are believed to be key actuators of the noradrenergic system 44 O’Donnell et al.2012 68 Wahis and Holt 2021. Importantly, NA signaling is assumed to play a pivotal role in AD pathology 71 Weinshenker 2018 46 Peterson and Li 2018 28 Holland et al.2021. Profound loss of noradrenergic neurons is a hallmark of AD that occurs early in the development of the disorder, likely accounting for some of the symptoms of AD, including loss of cognitive function 41 Matchett et al.2021. However, the mechanisms and signaling pathways connecting perturbed NA and cognition are only rudimentarily understood.

As astrocytes are tightly associated with NA signaling and AD is characterized by perturbed NA system, we set out to investigate whether the norepinephrine–astrocyte Ca2+ signaling axis was perturbed in ~15-month-old unanesthetized awake-behaving tg-ArcSwe mice compared to wildtype (WT) littermates. The tg-ArcSwe mice carry two mutations in the amyloid precursor protein gene, the Arctic (E693G) and Swedish (KM670/6701NL) mutations, and exhibit amyloid-β (Aβ) deposits, a hallmark of AD 40 Lord et al.2006 73 Yang et al.2011 36 Lillehaug et al.2014 47 Philipson et al.2009. The mice have a robust phenotype with emergence of Aβ plaques at 6–7 months of age both in the parenchyma and in the walls of blood vessels, which are biochemically similar to human Aβ plaques 47 Philipson et al.2009. The prominent cerebral amyloid angiopathy (CAA) in close proximity to astrocytic endfeet makes the model suitable for studying astrocytic signaling 73 Yang et al.2011.

We demonstrate an attenuated behaviorally induced Ca2+ signaling in cortical astrocytes during locomotion. As noradrenergic signaling is known to trigger astrocytic Ca2+ signaling in startle responses and locomotion, and pupil responses are regarded a faithful, although indirect, readout of noradrenergic signaling in the brain during physiological conditions 53 Reimer et al.2016 75 Zuend et al.2020 15 Costa and Rudebeck 2016, we compared the pupil responses to the astrocytic Ca2+ signals during running and startle and found a positive correlation in WT mice. No such correlation was present in the AD mice. Using the fluorescent extracellular norepinephrine sensor GRABNE 19 Feng et al.2019, we found attenuated norepinephrine signaling in relation to locomotion and startle responses in the AD mice, likely explaining the observed reduced astrocytic Ca2+ signaling. Such perturbed behaviorally induced norepinephrine-linked astrocytic Ca2+ signaling may account for some of the cognitive deficiencies observed in AD patients.

## Results

### Widespread reactive astrogliosis in tg-ArcSwe mice

AD transgenic mouse models with high levels of Aβ deposition exhibit reactive astrogliosis 54 Rodríguez-Arellano et al.2016. As astrogliosis might affect astrocytic Ca2+ activity 58 Shigetomi et al.2019 55 Sano et al.2019 and has not been quantified in neocortex of tg-ArcSwe mice yet, we measured glial fibrillary acidic protein (GFAP) expression by immunofluorescence and mRNA levels, and performed morphometric analyses of astrocytes in tg-ArcSwe mice and WT littermates (Figure 1).

Figure 1.

#### Widespread reactive astrogliosis in tg-ArcSwe mice.

(A) Representative micrographs of 15-month-old wildtype (WT) and tg-ArcSwe mice somatosensory cortex labeled with anti-GFAP antibodies (green) and anti-NeuN antibodies (blue). Scale bar: 40 μm. (B) Gfap mRNA expression was considerably higher in tg-ArcSwe mice (p<0.0001, n = 8 animals in both groups, Mann–Whitney U-test). (C) The area fraction of GFAP labeling was significantly higher in tg-ArcSwe mice compared to controls (n = 6 tg-ArcSwe and 4 WT littermates). (D) Astrocytes were isolated (inset in A) and analyzed with Simple Neurite Tracer (SNT) plugin in Fiji ImageJ. Tg-ArcSwe mice displayed an increase in the total number of branches (p=0.02), the length and volume of processes (p=0.0095 and 0.095, respectively, two-way ANOVA) of the GFAP-labeled astrocytes (n = 8 astrocytes from each animal; six tg-ArcSwe and four WT littermates), as well as an increased number of branching points and intersections at 16–32 μm distance from the nucleus (p<0.0001, two-way ANOVA).

We found a strong increase in levels of Gfap mRNA (81.58 ± 19.56 for WT; 512.13 ± 66.95 for tg-ArcSwe, p<0.0001, n = 8 animals in both groups) in tg-ArcSwe compared to WT littermates (Figure 1B). This was supported by a significantly higher GFAP labeling fraction in tg-ArcSwe animals compared to WT littermates (Figure 1C, 24.03% ± 4.80% in tg-ArcSwe, n = 6 mice vs. 4.77% ± 1.58% in WT, n = 4 mice, p=0.02). Morphometric analyses displayed a significantly higher number of labeled astrocytic processes (38.63 ± 3.61 in tg-ArcSwe, n = 6 mice; 23.84 ± 3.72 in WT, n = 4 mice; p=0.02), as well as total length (837.9 µm ± 55.9 µm for tg-ArcSwe; 506.0 µm ± 62.5 µm for WT; p=0.0095) and volume of processes (990.9 µm3 ± 93.8 µm3 for tg-ArcSwe; 401.5 µm3 ± 82.9 µm3 for WT; p=0.0095) in tg-ArcSwe mice compared to WT littermates (Figure 1D). In addition, astrocytes in the tg-ArcSwe mice exhibited significantly more branching points at the distance 16–32 µm from the nucleus (p<0.0001).

### Two-photon imaging of awake-behaving tg-ArcSwe mice

To characterize astrocytic Ca2+ signaling in awake tg-ArcSwe mice and nontransgenic littermates, we employed two-photon microscopy of cortical layer 1–3 astrocytes in the somatosensory cortex expressing GCaMP6f. Recombinant adeno-associated virus (rAAV) was used to deliver the genetic construct, and the GFAP promoter was used to target astrocytes (Figure 2A). Aβ plaques were visualized in vivo by methoxy-X04 delivered by intraperitoneal injection (Figure 2A). Methoxy-X04 enters the brain and specifically stains parenchymal Aβ plaques and cerebrovascular deposits 33 Klunk et al.2002, and is widely used for in vivo imaging of transgenic mice with Aβ 16 Delekate et al.2014 35 Kuchibhotla et al.2009 43 Meyer-Luehmann et al.2008. Both tg-ArcSwe mice and littermates were injected with methoxy-X04 to rule out any confounding effects of methoxy-X04 on astrocytic Ca2+ activity. Imaging was performed at ~30 Hz frame rate to capture fast populations of astrocytic Ca2+ transients with simultaneous surveillance video recording of mouse behavior, movement of the treadmill, as well as pupil diameter (Figure 2B) to monitor the level of arousal 53 Reimer et al.2016. The mice were allowed to spontaneously move on a custom-built disc-shaped treadmill, and all mice exhibited both periods of quiet wakefulness (absence of locomotion) and running (Figure 2—figure supplement 1). Astrocytic Ca2+ signals were analyzed using a newly developed event-based Ca2+ signal analysis toolkit, outlining so-called regions of activity (ROAs), combined with manually drawn regions of interest (ROIs) defining astrocytic subcompartments 7 Bjørnstad et al.2021 9 Bojarskaite et al.2020. This method detects Ca2+ signals in a pixel-by-pixel fashion by (1) estimating the noise level per pixel over time, (2) binarizing the data into signal/non-signal by noise level-based thresholds per pixel, and (3) connecting adjoining active pixels in space and time resulting in x-y-t ROAs. For the remainder of this article, we primarily report the ROA density, which is the fraction (in %) of the compartment analyzed with signal at any given time. See ‘Materials and methods’ for more details.

Figure 2.

#### Experimental setup and astrocytic Ca2+ signaling in quiet wakefulness.

(A) A craniotomy was performed and virus encoding GCaMP6f was injected into the somatosensory cortex of tg-ArcSwe mice and wildtype (WT) littermates. After 3 weeks of recovery, mice were habituated to head fixation on a disc-shaped treadmill, allowing the mice to move freely at will. Methoxy-X04 was injected 24 hr prior to imaging to visualize amyloid-β (Aβ) plaques. During imaging, both locomotor activity and pupil responses were recorded. At 150 s in 300 s recordings, the mice were subjected to an air puff to induce a startle response. (B) A representative imaging trial. Simultaneous recording of astrocytic Ca2+ signaling, pupil diameter, and mouse locomotion during quiet wakefulness, spontaneous locomotion, and air puff-induced startle. (C) The mice were ~15 months of age during experiments, at a time when they exhibited dense-core Aβ plaques. Left image: confocal micrograph of an Aβ plaque (methoxy-X04, blue) and astrocytes (anti-GFAP, green). Scale bar: 40 μm. Right image: electron micrograph showing a dense Aβ plaque (purple overlay), autophagic vacuoles (yellow overlay), and relatively normal neuropil morphology (green overlay). Scale bar: 2 μm. (D) Astrocytic Ca2+ signals during quiet wakefulness (absence of locomotion) in the form of regions of activity (ROAs) displayed in an x-y-t 3D rendering where red regions denote signal. Box-and-whisker plots representing overall Ca2+ signals in quiet wakefulness in tg-ArcSwe mice and littermates. (E) 3D visualization of the imaging plane relative to Aβ plaques, with lines representing shortest distance from plaque to region of activity (ROI). We found a low correlation between distance to nearest plaque and the overall level of astrocytic Ca2+ signaling in quiet wakefulness. (F) Example of pathological astrocytic Ca2+ wave emanating from an Aβ plaque. Methoxy-X04 visualized by 2PL (upper left) and time series from the same field of view (FOV) of astrocytic Ca2+ imaging with a pathological signaling event outlined in red. Scale bar: 50 μm.

A .csv file containing Ca2+ response data (region of activity [ROA] density and frequency), locomotion data, and covariates used for statistical modeling for quiet wakefulness.

Figure 2—figure supplement 1.

#### Locomotor behavior in wildtype (WT) and tg-ArcSwe mice.

▲ indicates air-puff. (A) Behavioral classification for each individual trial, grouped by mouse. See ‘Materials and methods’ for classification and filtering algorithm used. (B) Running wheel speed for each individual trial, grouped by mouse. 20 degrees of rotation of the running wheel per second was the threshold used to classify running. (C) Percentage of total time a mouse spent above a given speed across all trials. (D) Percentage of total time a mouse was running (>20°/s). Raw running data not filtered by the run classification algorithm (p=0.35). (E–G) Post-classification run episode characteristics.

Figure 2—figure supplement 2.

#### Robustness and sensitivity analyses for the uncoupling between pupil dilation and the region of activity (ROA) density rise rate.

All analyses are with respect to the interaction effect between pupil dilation and genotype. (A) Leave-one-trial-out analysis. The full model is repeatedly estimated removing one trial at the time from the dataset (as a type of jackknife resampling). The estimated interaction coefficients are shown as a histogram. Note that the histogram indicates less variation in the estimated interaction when comparing the variance or empirical quantiles of the histogram with the corresponding confidence interval obtained from the linear mixed effect model (see solid and dashed lines in plot B). (B) Leave-one-mouse-out analysis. In each iteration, one mouse, M10-M24, is removed from the estimation of the model. The estimated interaction effects with corresponding confidence intervals are then compared to the estimate and confidence interval from the full model (solid and dashed horizontal lines). Some mice influence the results more than others, but there appears to be no systematic bias associated with a single mouse. (C) The p-value for the interaction effect for different window lengths. The window length is measured in number of frames. It is clear that the significance of the interaction effect does not depend strongly on a particular choice of window length. In the analysis, we used a window of length 50 frames.

Figure 2—figure supplement 3.

#### GCaMP6f expression level assessed by GFP labeling.

(A) Median GFP-labeling fluorescence indicative of GCaMP6f expression was similar in the two groups.

Dots represent median values per mouse. (B) Mean area fraction of the field of view (FOV) with GFP labeling was comparable between the two groups.

Figure 2—figure supplement 4.

#### Regions of activity (ROAs) during quiet wakefulness exhibit the same characteristics in both genotypes.

During quiet wakefulness episodes, active regions detected by the ROA algorithm. Compartment association is determined by the center of the ROA. (A–C) Mean duration of active regions during quiet wakefulness in different cellular compartments. Right: GLME statistical analysis estimates. (D–F) Mean maximum size of a detected region. Right: GLME statistical analysis estimates.

The potential differences in astrocytic Ca2+ signaling between the two genotypes during spontaneous runs and startle responses, as well as the relationship between Ca2+ signaling and pupil dilation, were investigated through statistical models. These results are presented in the following sections, and the methods are described in detail in ‘Statistical analyses’. The Ca2+ signal responses were modeled by linear mixed effects regression models. The coefficients of primary interest were the effect of genotype, and the interaction between the effect of pupil dilation and genotype (i.e., the parameter indicating potential differences in the effect of pupil dilation on Ca2+ signaling between the two genotypes). In addition, the models included some technical and biological covariates, like the μm per pixel (images acquired with two different levels of magnification), the depth of the measurements, and the maximal speed of mouse locomotion in the time window studied. Further, the models include random intercepts for each mouse. The sensitivity of our result to these modeling choices was assessed by various robustness checks, primarily by iteratively removing individual mice and trials to see to which degree this influenced p-values and estimates (see ‘Statistical analyses’ and Figure 2—figure supplement 2).

### Astrocytes close to Aβ plaques express GCaMP6f

The tg-ArcSwe mice were imaged at ~15 months of age. At this age, they present with Aβ plaques throughout the cortical mantle, and score poorly on behavioral tasks (14 Codita et al.2010; 36 Lillehaug et al.2014; 40 Lord et al.2006; Figure 2C). Aβ plaques were characterized by loss of cells and severely perturbed tissue morphology, including autophagic vacuoles (Figure 2C). Even so, relatively normal cellular morphology was present at short distances away from Aβ plaques, and astrocytes faithfully expressed the GCaMP6f Ca2+ sensor 3 weeks after viral transduction (Figure 2A and C, Figure 2—figure supplement 3). Of note, GCaMP6f expression as judged by average GFP labeling and area fraction of GFP labeling in confocal micrographs were similar in the two groups (Figure 2—figure supplement 3).

### Astrocytic Ca2+ signals in quiet wakefulness are preserved in tg-ArcSwe mice

In quiet wakefulness (defined as absence of locomotion), the gross level of astrocytic Ca2+ signaling was similar in mutant mice and their littermates as measured by ROA frequency, ROA density (the active fraction of a compartment), as well as event size and duration in the full FOV and across the different astrocytic subcompartments (Figure 2D, Figure 2—figure supplement 4). We were neither able to detect a clear correlation in overall astrocytic Ca2+ signaling measured by ROA density and the distance from nearest Aβ plaque (in 3D) (slope = –0.00020, Figure 2E). However, we found some examples of long-lasting pathological Ca2+ waves as reported previously in anesthetized mice 16 Delekate et al.2014 35 Kuchibhotla et al.2009. Such Ca2+ waves were found in 10–15% of recordings from tg-ArcSwe mice (Figure 2F), and the number of such clear pathological events was very low compared to the overall astrocytic Ca2+ signaling.

### Uncoupling between pupil dilation and astrocytic Ca2+ responses during spontaneous running in tg-ArcSwe mice

Locomotor behavior is known to be strongly correlated with astrocytic Ca2+ signaling 45 Paukert et al.2014 9 Bojarskaite et al.2020 61 Srinivasan et al.2015, putatively through the activation of the noradrenergic and cholinergic neuromodulatory systems in conjunction with local network activity 32 Kjaerby et al.2017. To investigate whether the physiological astrocytic Ca2+ responses were preserved in the tg-ArcSwe mice, they were allowed to move freely on a custom built disc-shaped treadmill 9 Bojarskaite et al.2020. All mice exhibited both running and behavioral quiescence, and the level of running between the two genotypes was comparable (Figure 2—figure supplement 1). Running was accompanied by an increase in pupil size and a brisk increase in astrocytic Ca2+ signaling typically involving most of the astrocytes in the FOV in both genotypes (Figure 3A,B). Overall, astrocytic Ca2+ responses to spontaneous locomotion were somewhat reduced in tg-ArcSwe mice compared to WT littermates: When astrocytic Ca2+ signals were analyzed using a linear mixed effects regression model (Table 1), a lower ROA density rise rate was found in tg-ArcSwe mice when assessing the full FOV (0.31 in WT vs. 0.20 in tg-ArcSwe, p=0.032), and astrocytic processes (0.32 in WT vs. 0.20 in tg-ArcSwe, p=0.032), whereas the Ca2+ responses were not significantly different in astrocytic somata and endfeet (0.46 in WT vs. 0.39 in tg-ArcSwe, p=0.23, and 0.35 in WT vs. 0.30 in tg-ArcSwe, p=0.11, respectively) (Figure 3C). Max ROA density values in WT vs. tg-ArcSwe were significantly different when assessing the full FOV (0.63 vs. 0.42, p=0.033), near significantly different when assessing astrocytic processes and endfeet (0.72 vs. 0.53, p=0.053 for processes and 0.69 vs. 0.51, p=0.068 for endfeet), and not significantly different for astrocytic somata (0.81 vs. 0.71, p=0.25) (Figure 3C). For mean ROA density values, see Figure 3—figure supplement 1.

Figure 3.

#### Uncoupling of pupil dilation and astrocytic Ca2+ responses during spontaneous running.

(A) Left: Representative x-y-t rendering of regions of activity (ROAs) during quiet wakefulness and spontaneous locomotion in wildtype (WT) littermates and tg-ArcSwe mice, where red regions denote signal. A spontaneous run occurred at about 60 s. The slight increase in ROA density at 25 s is due to small locomotor activity that was not large enough to be defined as a run. Right: mean time course of (top-to-bottom) locomotion expressed as wheel speed, relative increase in pupil diameter, and astrocytic Ca2+ signals expressed as ROA density. The locomotion trace shows mean response across all trials with 95% confidence interval, pupil trace, and ROA density traces show median across all trials ± median absolute deviation. (B) Same as (A), but for tg-ArcSwe mice. A spontaneous run occurred at about 40 s. (C) Hierarchical plots showing median value per genotype in upper level, median value per mouse in middle level, and median level per trial in lower level, in the two genotypes. Upper row: plots showing the median levels of ROA density rise rate in (left-to-right) the full FOV, astrocytic somata, astrocytic processes, and astrocytic endfeet. Lower row: Same as upper row, but showing median levels of max ROA density. (D) Scatterplots of ROA density rise rate (left) and max ROA density (right) vs. relative increase in pupil size upon spontaneous running. See Table 2 for statistical analyses and p-values.

A .csv file containing Ca2+ response data (region of activity [ROA] density and frequency), locomotion data, and covariates used for statistical modeling for spontaneous running.

Figure 3—figure supplement 1.

#### Mean region of activity (ROA) activation in active wakefulness.

(A) Hierarchical plot showing medians of mean region of activity (ROA) density per trial, mouse and genotype in the field of view (FOV), and astrocytic subcellular compartments during spontaneous runs.

(B) Same as (A) but during startle responses. (C) Scatterplots of mean ROA density vs. relative increase in pupil size in spontaneous running and startle responses.

Table 1.

#### Estimated statistical model for region of activity (ROA) density rise rate in the field of view (FOV) during spontaneous running.

μm per pixel has a small/no effect, but we include it for consistency. Deeper measurements tend to have lower ROA density rise rate. Given that the mice are running, a higher speed does not appear to be associated with a higher ROA density rise rate (n = 109, in 10 mice).

| Estimate | p-Value
--- | --- | ---
Intercept (TG) | 0.1969 |
Genotype (WT) | 0.1085 | 0.0318
μm per pixel | –0.0055 | >0.10
Depth | –0.0432 | 0.0002
Maximum speed | 0 | >0.10

Pupil responses are known to be a faithful indirect indicator of activity in the locus coeruleus in mice, even though also the cholinergic neuromodulatory system plays a role for sustained pupil dilation 53 Reimer et al.2016. As norepinephrine is known to be a potent trigger of astrocytic Ca2+ signaling 5 Bekar et al.2008 61 Srinivasan et al.2015 45 Paukert et al.2014, one would expect to find a correlation between pupil dilations and astrocytic Ca2+ signals. This was indeed the case for WT mice (Table 2): When comparing ROA density rise rate and pupil dilation, we found a positive slope in line with previous reports (75 Zuend et al.2020; Figure 3D). In transgenic mice, this correlation was lost, or even reversed (Figure 3D, Table 2: slope of 0.80 in WT vs. –0.35 in tg-ArcSwe, p=0.007), demonstrating an uncoupling between pupil responses and astrocytic Ca2+ signaling in tg-ArcSwe mice. Similarly, when assessing max ROA density, we found a positive slope in WT, which was lost in tg-ArcSwe (Figure 3D, 2.06 in WT vs. –0.58 in tg-ArcSwe, p=0.00039). Similar slopes were found when comparing mean ROA density vs. pupil dilation (Figure 3—figure supplement 1C).

Table 2.

#### Estimated statistical model for region of activity (ROA) density rise rate in the field of view (FOV) and pupil dilation during spontaneous running.

When including pupil dynamics in the model (cf. Table 1 without pupil dynamics), the effect of genotype was no longer significant. However, the relationship between pupil dilation and ROA density rise rate was different in the two genotypes, with the wildtype (WT) mice having a significantly higher slope. μm per pixel and depth of the recording had no significant effect on the model, but we included them for consistency. The slope for the WT group is found by adding the slope in the TG group with the interaction term (e.g., –0.3537 + 1.1562 = 0.80) (n = 44, in 10 mice).

| Estimate | p-Value
--- | --- | ---
Intercept (TG) | 0.2297 |
Genotype (WT) | –0.1245 | >0.10
Pupil dilation (TG) | –0.3537 | >0.10
Pupil dilation × genotype (WT) | 1.1562 | 0.007
μm per pixel (binary) | 0.0306 | >0.10
Depth | –0.0409 | 0.0996
Maximum speed | –0.0011 | 0.0351

### Uncoupling between pupil dilation and astrocytic Ca2+ responses during startle in tg-ArcSwe mice

Another main trigger for astrocytic Ca2+ signals are startle responses that also are mediated through an activation of the noradrenergic system 17 Ding et al.2013 61 Srinivasan et al.2015. Even though typically triggering running, the startle response could also trigger freezing behavior and is thought to activate different subcortical networks than spontaneous locomotor behavior 12 Caggiano et al.2018 20 Ferreira-Pinto et al.2018 21 Grillner and El Manira 2020. To investigate whether the startle-evoked astrocytic Ca2+ responses were preserved in the tg-ArcSwe mice, they were subjected to 10 air puffs delivered at 10 Hz directed to the vibrissa, nasal, and facial region contralaterally to the recording side once per trial at 150 s in a 300 s two-photon imaging recording. Trials in which the mouse was spontaneously running at or immediately before the air puff were excluded from the analyses. We found no signs of habituation to the stimulus in terms of behavioral response (Figure 4—figure supplement 1). Interestingly, tg-ArcSwe mice were more prone to start running during startle responses than WT littermates (Figure 2—figure supplement 1), consistent with previous reports of enhanced startle response in other mouse models of AD 42 McCool et al.2003. The level of pupil dilation was however similar in the two genotypes (0.17 vs. 0.12 relative increase in pupil size in WT vs. tg-ArcSwe, respectively, p=0.36, 86 trials) (Figure 4A and B). Regarding astrocyte Ca2+ response to startle, we found no differences between the genotypes. When modeled with a mixed effects linear regression model (Table 3), for ROA density rise rate, Ca2+ responses in the full FOV and in all astrocytic subcompartments were statistically similar between the genotypes (FOV: 0.34 in WT vs. 0.23 in tg-ArcSwe, p=0.09; somata: 0.25 in WT vs. 0.23 in tg-ArcSwe, p=0.8; processes: 0.42 in WT vs. 0.31, p=0.09; endfeet: 0.36 in WT vs. 0.25, p=0.10) (Figure 4C). Similarly, max ROA density was similar for the full FOV (0.36 in WT vs. 0.35 in tg-ArcSwe, p=0.40), astrocytic somata (0.37 in WT vs. 0.35 in tg-ArcSwe, p=0.83), astrocytic processes (0.41 in WT vs. 0.34 in tg-ArcSwe, p=0.39), and astrocytic endfeet (0.41 in WT vs. 0.35 in tg-ArcSwe, p=0.48) (Figure 4C). For mean ROA density values, see Figure 3—figure supplement 1.

Figure 4.

#### Uncoupling of pupil dilation and astrocytic Ca2+ responses during startle.

(A) In wildtype (WT) mice, in an imaging trial of 300 s duration, mice were subjected to an air puff to the face/vibrissa contralateral to the imaging window at 150 s. This caused an increase in locomotor activity, pupil dilation, and a pronounced increase in astrocytic Ca2+ signaling. Region of activity (ROAs) are presented as an x-y-t 3D rendering where red regions denote signal. Locomotion trace shows mean response across trials with 95% confidence interval, pupil trace, and ROA density trace show median activation across trials ± median absolute deviation. (B) Same as (A), but for tg-ArcSwe mice. (C) Hierarchical plots showing median value per genotype in upper level, median value per mouse in middle level, and median level per trial in lower level, in the two genotypes. Upper row: plots showing the median level of the ROA density rise rate in (left-to-right) the full field of view (FOV), astrocytic somata, astrocytic processes, and astrocytic endfeet. Lower row: same as upper row but for max ROA density. (D) Scatterplots of ROA density rise rate (left) and max ROA density (right) vs. relative increase in pupil size. See Tables 3 and 4 for statistical analyses and p-values.

A .csv file containing Ca2+ response data (region of activity [ROA] density and frequency), locomotion data, and covariates used for statistical modeling for startle responses.

Figure 4—figure supplement 1.

#### Behavioral response to air puff across trials.

(A) Neither tg-ArcSwe nor wildtype (WT) mice displayed any noticeable habituation to air puffs across all trials and experimental days.

(B) Within a single experimental day, there was no habituation to air puffs across the different trials.

Figure 4—figure supplement 2.

#### The dynamic range of pupil dilation and constrictions in the two genotypes.

Dots represent median values per mouse.

Figure 4—figure supplement 3.

#### Method for assessing region of activity (ROA) density and pupil dynamics.

(A) Computation of the region of activity (ROA) density rise rate, and max ROA density, in particular startle period.

The vertical dashed lines indicate the time window that had the maximal increase in ROA density, which in this example took place over 49 frames. (B) Computation of the pupil dilation in a particular startle period, here we find (0.268–0.202)/0.202 = 0.33.

Table 3.

#### Estimated statistical model for region of activity (ROA) density rise rate in the field of view (FOV) during startle.

Depth from the cortical surface and maximal speed during the startle response had a significant effect in the model, whereas the μm per pixel (level of magnification of the image recording) had no significant effect (n = 117 trials, in 11 mice).

| Estimate | p-Value
--- | --- | ---
Intercept (TG) | 0.2295 |
Genotype (WT) | 0.111 | 0.0933
μm per pixel (binary) | –0.0034 | >0.10
Depth | –0.0696 | 0.0004
Maximum speed | 0.0009 | 0.022

Table 4.

#### Estimated statistical model for region of activity (ROA) density rise rate in the field of view (FOV) and pupil dilation during startle.

When including the pupil dilation in the statistical model, the two genotypes are clearly different. Moreover, depth from the cortical surface had a significant effect in the model, whereas optical zoom and maximum speed had a nonsignificant effect (n = 95 trials, in 11 mice). The slope for the wildtype (WT) group is found by adding the slope in the TG group with the interaction term (e.g., –0.0129 + 0.9235 = 0.91).

| Estimate | p-Value
--- | --- | ---
Intercept (TG) | 0.2199 |
Genotype (WT) | –0.0279 | >0.1
Pupil dilation (TG) | –0.0129 | >0.1
Pupil dilation × genotype (WT) | 0.9235 | 0.00043
Optical zoom | 0.0036 | >0.1
Depth | –0.0679 | 0.00063
Maximum speed | 0.0006 | 0.0587

However, the relationship between astrocytic Ca2+ responses and pupillary responses were different in the two genotypes (Figure 4D). WT mice displayed a clear positive slope, while tg-ArcSwe exhibited a slope close to zero (0.91 vs. -0.013 p = 0.00043 for ROA density rise rate, and 1.27 vs. 0.083 p = 0.0043 and max ROA density, respectively), suggesting a potential uncoupling between pupillary and astrocytic Ca2+ responses during startle in tg-ArcSwe mice, similar to what was observed during spontaneous running (Figure 3D). Similar slopes were found when comparing mean ROA density vs. pupil dilation (Figure 3—figure supplement 1).

### Attenuated noradrenergic signaling in tg-ArcSwe mice

To investigate the cause of our observed attenuated astrocytic Ca2+ responses during locomotion and startle responses, and uncoupling of pupil responses and astrocytic Ca2+ responses, we took advantage of the newly developed genetically encoded fluorescent extracellular norepinephrine sensor, the GRABNE sensor 19 Feng et al.2019. The sensor was delivered by rAAV transduction of the GRABNE (version 1m) construct and targeted to membranes of neurons using the hSyn1 promoter. The mice were subjected to the same imaging protocols as for Ca2+ imaging, and we recorded norepinephrine fluorescence, pupil dynamics, and locomotor activity. Differences between the genotypes were modeled by linear mixed effects regression models as outlined for Ca2+ imaging (results presented in Tables 5 and 6). The coefficients modeled were the effect of genotype, and the interaction between the effect of pupil dilation and genotype (i.e., the parameter indicating potential differences in the effect of pupil dilation on norepinephrine signaling between the two genotypes). During both spontaneous running and startle responses, we observed increases in norepinephrine fluorescence in WT mice (Figure 5). This effect was less prominent in tg-ArcSwe mice, and similar to Ca2+ imaging, we found a decoupling between pupil responses and noradrenergic signaling. In WT mice, for spontaneous runs, the relative increase in pupil diameter explained a large proportion of change in norepinephrine signaling (Table 5). In tg-ArcSwe mice, pupil responses and norepinephrine signaling were less correlated. The results for startle responses were similar (Table 6). These findings suggest that the reduced Ca2+ signaling and uncoupling between pupil responses and astrocytic Ca2+ signaling in tg-ArcSwe mice could be due to attenuated noradrenergic release in cortex and lower correlation between pupil responses and noradrenergic signaling.

Table 5.

#### Estimated statistical model for change in GRABNE fluorescence and pupil dilation during spontaneous runs.

The coupling between pupil responses and noradrenergic signaling is weaker in tg-ArcSwe mice compared to wildtype (WT) littermates (p=0.043, n = 100 running episodes, 43 trials, in 7 mice). The slope for the WT group is found by adding the slope in the TG group with the interaction term (e.g., 0.0030 + 0.043 = 0.046).

| Estimate | p-Value
--- | --- | ---
Intercept (TG) | 0.003 |
Genotype (WT) | –0.0012 | >0.1
Pupil dilation (TG) | 0.00022 | >0.1
Pupil dilation × genotype (WT) | 0.043 | 0.042
Increase in speed | –0.00000401 | >0.1

Table 6.

#### Estimated statistical model for change in GRABNE fluorescence and pupil dilation during startle.

The coupling between pupil responses and noradrenergic signaling during startle responses is weaker in tg-ArcSwe mice compared to wildtype (WT) littermates (p=0.037, n = 59 trials, in 7 mice). The slope for the WT group is found by adding the slope in the TG group with the interaction term (e.g., 0.0078 + 0.037 = 0.045).

| Estimate | p-Value
--- | --- | ---
Intercept (TG) | 0.0078 |
Genotype (WT) | –0.0063 | >0.1
Pupil dilation (TG) | –0.013 | >0.1
Pupil dilation × genotype (WT) | 0.037 | 0.036
Increase in speed | –0.0000029 | >0.1

Figure 5.

#### Attenuated norepinephrine responses in tg-ArcSwe mice.

Mice expressing the fluorescent extracellular norepinephrine sensor GRABNE were subjected to two-photon imaging trials similar to the astrocytic Ca2+ imaging experiments. (A) Mean traces ± standard error of locomotor responses, pupil responses and GRABNE fluorescence during spontaneous runs (left column) and during startle (right column). (B) Hierarchical plots showing median ΔF/F0 of GRABNE fluorescence per (from bottom to top) trial, mouse and genotype, in spontaneous runs and startle responses.

A .mat file containing GRABNE fluorescence traces during spontaneous running and startle as well as locomotion data and covariates used for statistical modeling.

## Discussion

Astrocytic Ca2+ signaling is emerging as a key component of signal processing in the brain and figures prominently in brain state transitions and memory formation 50 Poskanzer and Yuste 2016 49 Poskanzer and Yuste 2011 3 Adamsky et al.2018 9 Bojarskaite et al.2020 66 Vaidyanathan et al.2021. Aberrant astrocytic Ca2+ signaling could hence be implicated in the perturbed cognition seen in dementia. Indeed, previous studies have shown increased astrocytic Ca2+ signaling and spreading pathological Ca2+ waves in AD mouse models 35 Kuchibhotla et al.2009 16 Delekate et al.2014 52 Reichenbach et al.2018. Based on these studies in anesthetized animals and other studies 25 Haughey and Mattson 2003 37 Lim et al.2013 1 Abramov et al.2003 2 Abramov et al.2004 67 Verkhratsky 2019, the current concept is that aberrant signaling to some degree is spatially coupled to amyloid deposits – the pathological hallmark of AD.

Methodological advances now allow astrocytic Ca2+ signals to be studied in awake animals. Benefitting from this opportunity, we show that tg-ArcSwe mice sustain a pattern of behaviorally induced astrocytic Ca2+ signaling similar to that found in littermate controls. However, the signals are weaker than in controls and do not display the correlation with pupil responses typically seen in WT animals. The behaviorally induced Ca2+ signals bear no clear spatial correlation to Aβ plaques. We conclude that elimination of anesthesia unveils a new dimension of Ca2+ signaling, superimposed on the locally induced signals described in previous studies. The uniform attenuation of the behaviorally induced signals in tg-ArcSwe mice and their uncoupling from pupil responses point to a potential perturbation of neuromodulatory control.

Astrocytic Ca2+ signaling in awake-behaving mice is dominated by norepinephrine-induced astrocytic Ca2+ signals across the cortical mantle in relation to locomotor or startle responses 61 Srinivasan et al.2015 17 Ding et al.2013. The downstream effects of these synchronized, global Ca2+ signals coupled to arousal are not fully understood, but they may play a role in altering the levels of extracellular K+ 69 Wang et al.2012, release of gliotransmitters such as glutamate or ATP 32 Kjaerby et al.2017 26 Haydon and Nedergaard 2015, or metabolic supply 75 Zuend et al.2020. Independent of what the exact downstream mechanisms are, there are reasons to believe that astrocytes serve as actuators of the noradrenergic neuromodulatory system, presumably exerting noradrenergic effects on neural network processing and ultimately affecting cognitive function 74 Ye et al.2020 28 Holland et al.2021 50 Poskanzer and Yuste 2016.

The perturbed astrocytic Ca2+ response during locomotion in AD mice seems to be at least partially explained by a blunted norepinephrine increase upon startle and locomotion as evaluated by GRABNE fluorescence. The primary norepinephrine nucleus in the brain, locus coeruleus, is known to be affected early in AD, both in humans and in animal models 71 Weinshenker 2018 29 Jacobs et al.2021 10 Braak and Del Tredici 2011, and perturbed noradrenergic function is likely an important factor both in disease progression in AD patients and accounting for the cognitive decline of AD patients 71 Weinshenker 2018 46 Peterson and Li 2018 28 Holland et al.2021. Restoring noradrenergic signaling in AD mice with perturbed noradrenergic signaling pharmacologically or by pharmacogenetics has been demonstrated to be beneficial 28 Holland et al.2021. In our study, we find a largely retained level of pupil response in the AD mice during startle responses and a similar dynamic range of pupil responses (Figure 4—figure supplement 2), which under physiological circumstances would suggest a functional noradrenergic system, as a strong correlation between locus coeruleus activity and pupil responses has been established 53 Reimer et al.2016 15 Costa and Rudebeck 2016. The underlying connectivity that enables the pupils to be faithful readouts of the noradrenergic system is to the best of our knowledge not fully established, even though spinal projections from locus coeruleus have been demonstrated 23 Hancock and Fougerousse 1976 15 Costa and Rudebeck 2016 38 Liu et al.2017. Our finding of a decoupling between norepinephrine-mediated astrocytic Ca2+ signaling and arousal-induced pupillary dilation could be due to perturbed connectivity between locus coeruleus and the relevant nuclei or projections controlling activity in the superior cervical ganglion and consequently the sympathetic projections to the pupils.

Given the clear differences in norepinephrine signaling between the two genotypes, one would have expected even more pronounced differences in astrocytic Ca2+ signaling. This discrepancy could possibly be explained by either the lack of sensitivity of the GRABNE sensor in the current experimental paradigm to pick up small increases in cortical norepinephrine or potentially by a compensatory upregulation of adrenergic receptors or other receptors in reactive astrocytes in the tg-ArcSwe mice. In line with this, Reichenbach et al. demonstrated an age/disease stage-related upregulation of P2R1-receptors in an AD mouse model that correlated with the level of reactive astrogliosis 52 Reichenbach et al.2018. Moreover, potentially the expression of α1-adrenergic receptors in astrocytes is changed in AD, similar to reports of how α2 and β2 receptors are upregulated in microvessels and neurons in AD patients, respectively 30 Kalaria et al.1989 31 Kalaria and Harik 1989. Future studies probing the expression levels and patterns of adrenergic receptors and downstream signaling pathways in astrocytes in AD models are warranted.

Although still sparsely investigated, both attenuated and increased Ca2+ signaling have been demonstrated in various models of reactive astrogliosis, but so far without a clear understanding of the mechanisms involved 58 Shigetomi et al.2019. Widespread reactive astrogliosis is nonetheless highly likely to perturb the physiological signaling of astrocytes as prominent changes in the expression of key receptors and molecules of intracellular pathways are known to occur 22 Habib et al.2020 18 Escartin et al.2021. In this study, we employed 15-month-old tg-ArcSwe mice and WT littermates. These mice exhibit a relatively high plaque burden and represent a relatively advanced stage of AD. It is possible that the attenuated astrocytic Ca2+ responses we observe is a feature of advanced disease, and not earlier stages of the disease development. Similarly, the sparsity of pathological ‘hyperactive’ astrocytic Ca2+ signaling events in our data compared to previous studies 16 Delekate et al.2014 35 Kuchibhotla et al.2009 52 Reichenbach et al.2018 could stem from the other studies being performed at ages that represents earlier stages in the disease progression. Further studies across AD models and stages of disease progression mice are warranted to conclude.

Previous studies have shown mixed results regarding an effect of plaque–astrocyte distance, with reports of both hyperactivity close to plaques 16 Delekate et al.2014 and no such correlation except for plaques serving as initiation sites for pathological intercellular Ca2+ waves 35 Kuchibhotla et al.2009. In our awake-behaving mice, we found no clear correlation between overall Ca2+ activity and 3D distance reconstruction of plaque positions relative to the imaging plane (Figure 2E). The apparent discrepancies in the literature and this study may be due to different AD mouse models, different age groups investigated, different Ca2+ indicators employed, or lastly, the effects of removing anesthesia allowing for a much richer repertoire of behaviorally induced astrocytic Ca2+ signaling to emerge, effectively masking a potential weak correlation.

Astrocytes have a highly complex and specialized morphology, and this morphology is known to change in reactive astrogliosis 18 Escartin et al.2021. The majority of astrocytic processes are much smaller than what can be clearly delineated by non-super-resolution optical microscopy, but to what extent these small processes are altered in reactive astrogliosis is unknown, although the astrocytic territories are known to be preserved 72 Wilhelmsson et al.2006. Moreover, GFAP, which was used as the promoter for the GCaMP6f expression, was highly upregulated in tg-ArcSwe mice. This difference in GFAP expression did not translate into obvious differences in sensor expression levels as judged by GFP labeling of brain slices (Figure 2—figure supplement 3). We cannot rule out that gliosis-induced morphological changes in combination with potential subtle differences in GCaMP6f expression and potentially an elevated baseline intracellular Ca2+ concentration 35 Kuchibhotla et al.2009 may influence our results. However, the ability for the whole FOV to be activated during spontaneous running and startle responses and the lack of correlation (if anything a negative correlation) between distance to the nearest Aβ plaque and gross level of Ca2+ signaling suggest that our findings are not due to degree or pattern of sensor expression.

This study underscores the importance of studying astrocytic Ca2+ signals in unanesthetized mice and to carefully consider animal behavior when interpreting astrocytic Ca2+ dynamics. By lifting the confounding effects of anesthesia, we found that the astrocytic hyperactivity previously reported in AD mouse models was only a part of the total picture. At first glance, the physiological Ca2+ responses were remarkably well preserved and not characterized by a general increase in astrocytic Ca2+ signaling. However, behavior like quiet wakefulness and locomotion are not static entities, and the degree of activation of all relevant parameters needs to be taken into account with statistical modeling to be able to conclude if there are relevant differences between the genotypes. We were able to demonstrate attenuated Ca2+ dynamics and norepinephrine signaling in the tg-ArcSwe mice, and given the growing spectrum of roles ascribed to the norepinephrine–astrocyte Ca2+ signaling in higher brain functions, the present findings may highlight one potential cause for the cognitive decline of AD patients.

## Materials and methods

### Animals

Tg-ArcSwe mice carry a human APP cDNA with the Arctic (p. E693G) and Swedish (p. KM670/671NL) mutations where the human AβPP gene is inherited only from male mice to ensure a more uniform onset of Aβ deposition 36 Lillehaug et al.2014. The transgenic animals develop parenchymal Aβ plaques from 6 months of age and CAA from 8 months of age 73 Yang et al.2011 40 Lord et al.2006. Six tg-ArcSwe mice and five WT littermates (both males and females) were used in this study. Genotyping was done as previously described with primers annealing to the Thy1 promoter and the human APP transgene 40 Lord et al.2006. The animals were housed under standard conditions with 12 hr dark–light cycles and unrestricted access to food and water. All animal procedures were in accordance with the National Institutes of Health Guide for the care and use of laboratory animals and approved by the Norwegian Food Safety Authority (project number: FOTS #11983).

### Viral transduction and delivery of fluorophores

Serotype 2/1 rAAV from plasmid construct pAAV-GFAP-GCaMP6f 13 Chen et al.2013 was generated (rAAV titers about 1.0–6.0 × 1012 viral genomes/ml) and used for visualizing astrocytic Ca2+ signaling. GCaMP6f was amplified by PCR from pGP-CMV-GCaMP6f with 5′ Bam_HI and 3′ _Hin_dIII, and sub-cloned into the recombinant rAAV vector pAAV-6P-SEWB (57 Shevtsova et al.2005) for generating pAAV-_SYN-GCaMP6f. The human GFAP promoter 27 Hirrlinger et al.2009 was inserted with Mlu_I and _Bam_HI into pAAV-_SYN-GCaMP6f construct for obtaining pAAV-GFAP-GCaMP6f. Serotype 2/1 rAAVs from pAAV-GFAP-GCaMP6f were produced 63 Tang et al.2015 and purified by AVB Sepharose affinity chromatography 60 Smith et al.2009, following titration with real-time PCR (rAAV titers about 1.0–6.0 × 1012 viral genomes/ml, TaqMan Assay, Applied Biosystems). For norepinephrine imaging experiments, mice were injected with ssAAV-9/2-hSyn1-GRAB-NE1m-WPRE-hGHp(A) (University of Zürich Viral Vector Facility). To visualize Aβ plaques, 7 mg/g methoxy-X04 59 Skaaraas et al.2021 dissolved in 0.1 M phosphate buffer saline (PBS) was injected intraperitoneally 24 hr prior to imaging. Also, WT littermates were injected with methoxy-X04 to rule out any confounding effects in one of the groups. We found that one injection provided enough signal to outline Aβ plaques for up to 3 days.

### Surgical preparation

Mice were anesthetized with isoflurane (3% for initiation, then 1–1.5% for maintenance) in room air enriched with 20% pure oxygen, and given buprenorphine 0.1 mg/kg s.c. preemptively for analgesia. Bupivacaine was administered subcutaneously over the skull and left for 10 min before a boat-shaped skin flap was removed. After removing the skin, a 2.5-mm-diameter craniotomy was drilled over the somatosensory cortex with center coordinates 3.5 mm lateral and –1.5 mm posterior to Bregma. Virus was injected (70 nl at 35 nl/min at 200 µm depth from the brain’s surface) at three evenly spaced locations positioned to stay clear of large blood vessels, and a glass plug consisting of two coverslips glued together was placed in the craniotomy, slightly pressing the dura to prevent dural overgrowth 9 Bojarskaite et al.2020. The surrounding area of the skull was sealed with cyanoacrylate glue and a layer of dental cement. Postoperatively the mice were given meloxicam 2 mg/kg for 2 days. Only animals with normal postoperative recovery were included in the study. Mice were left to recover for a minimum of 2 weeks before habituation to head-fixation and imaging.

### Two-photon microscopy

After recovery, the animals were imaged in layers 1–3 of the barrel cortex using a two-photon microscope (Ultima IV from Bruker/Prairie Technologies, Nikon 16 × 0.8 NA water-immersion objective model CFI75 LWD 16XW, Spectra-Physics InSight DeepSee laser, Peltier cooled photomultiplier tubes model 7422PA-40 by Hamamatsu Photonics K.K.). An excitation wavelength of 990 nm was used for GCaMP6f imaging, and separate recordings were made using 790 nm excitation light for imaging methoxy-X04, and 890 nm for Texas-Red labeled dextran for visualization of the vasculature. GRABNE fluorescence was recorded at 920 nm excitation wavelength. Image time series were recorded at 30 Hz. 3D volume recordings of the morphology and plaque locations were performed in a subset of experiments. All images were recorded at 512 × 512 pixels with a resolution of approximately either 0.42 or 0.67 μm per pixel (FOV of 215 × 215 μm, or 343 × 343 μm, respectively). For imaging, the mice were head-fixed to a custom-built stage that allowed free locomotion on a wheel attached to a rotary encoder. During experiments, the mice were illuminated by an IR LED diode and monitored with an IR-sensitive camera. A second camera was used to capture the pupil dynamics (see below). Air puffs to elicit startle responses were delivered by a Picospritzer III (Parker). Instrument synchronization and data acquisition were performed with a custom-made LabVIEW 2015 (National Instruments) virtual instrument.

### Behavioral analysis

Mouse locomotion was recorded with a rotary encoder connected to the running wheel. Locomotion signal was captured using a National Instruments data acquisition card using a counter task in the NI Max software, activated through a custom LabVIEW (2015) VI. Data were processed with custom MATLAB scripts to classify running and quiet wakefulness. Criteria for run and quiet wakefulness episodes were validated by manual observation and defined as follows: running was defined as continuous segment of at least 4 s forward wheel motion at over 30°/s, no movement faster than 20°/s the last 10 s prior to start of running. Spontaneous running was not defined within 30 s following air puff. Quiet wakefulness was defined as continuous segment of at least 10 s duration with less than 2°/s locomotion, as well as no locomotor activity faster than 2°/s for 15 s before segment start. Quiet wakefulness episodes were not defined within 30 s following air puff. Quiet wakefulness periods were reviewed manually through the IR-surveillance video recordings to ensure animals were awake when sitting still.

### Pupillometry

Pupil size was recorded with a Basler Dart USB camera (daA1600-60um) with a 25 mm fixed focal length lens and 2× fixed focal length lens extender (Edmund optics, items #59–871 and #54–356). As two-photon imaging must be conducted in the dark where the pupil would be fully dilated, we illuminated the eye using 470 nm blue light fiber, with a shielding to avoid light contamination, to slightly constrict the pupil. Pupil size was delineated for time periods of isolated runs and startle using a custom-built tool developed in MATLAB 2020a. The ratio of pupil to eye-size ratio was calculated and used in the analysis.

### Image processing and analysis

Recordings were exclusively processed using MATLAB 2018a to 2020b. Imaging data were corrected for movement artifacts using the NoRMCOrre movement correction software 48 Pnevmatikakis and Giovannucci 2017. We used our recently published imaging analysis toolbox ‘Begonia’ to remove motion artifacts, mark ROIs, and perform event-based Ca2+ signal detection (ROA method). Astrocytic Ca2+ signals were detected using the ROA method, and the ROA density (i.e., fraction [in %] of the compartment analyzed being active) and ROA frequency in the whole FOV or in anatomical subcompartments were calculated. Methods and algorithms related to this processing pipeline are described in detail in 7 Bjørnstad et al.2021. In brief, movement-corrected image time-series data of astrocytic Ca2+ activity were preprocessed by variance stabilization (square root transform) and smoothing in the temporal and spatial domain to achieve a signal-to-noise ratio of at least 9. All image time series were for this reason smoothed with 14 frames in the temporal domain (0.45 s) and two pixels in the x–y domain. Four times the median standard deviation of the variance stabilized data plus the baseline value for a pixel was then used to calculate a threshold for signal detection for each pixel. Once the time series had been binarized into signal and non-signal, adjoining x–y–t voxels (in time and space) were connected to ROAs. ROAs with a spatial extent of less than three pixels per frame and duration of less than three frames in time were not included in the analyses. For GRABNE~ fluorescence, data were first downsampled from 30 Hz to 5 Hz and corrected for vascular artifacts (dips in fluorescence) due to increases in blood flow by software kindly shared by Prof. Kira Poskanzer at University of California San Francisco with colleagues and Prof. Guoqiang Yu with colleagues at Virginia Tech.

### 3D Aβ plaque mapping

Before any recording, z-stacks of images at 5 μm intervals taken at 512 × 512 pixels were recorded while illuminating the tissue with 790 nm laser light. The stacks started approximately at or below dura mater and extended 100–200 μm down. All other recordings were undertaken inside this mapped volume to ensure plaques outside the imaging plane were accounted for. For each imaged time series inside the volume, we recorded an additional single image at 790 nm to ensure the precise locations of plaques were known. Plaque locations were detected by binarizing the stack at a manually set threshold at which the morphology of plaques was visible. The binarized single images were used to manually align the 2D time-series data with the 3D plaque volume. The resulting 3D binary image had small points removed using the MATLAB function bwareaopen. Distances from plaque to nearest ROI (RP) were calculated using the 3D distance formula

R P=√(x r o i−x p l q)2+(y r o i−y p l q)2+(z r o i−z p l q)2

### Tissue processing

All animals were sacrificed after final imaging procedures at the age of 18 months. The animals were anesthetized with Isofluran Baxter (IsoFlo, Abbot Laboratories) and intraperitoneally injected with ZRF mixture (zolazepam 3.8 mg/ml, tiletamine 3.8 mg/ml, xylazine 0.45 mg/ml, and fentanyl 2.6 µg/ml) before transcardial perfusion with 4°C 2% dextran in 0.1 M phosphate buffer (PB) for approximately 30 s, immediately followed by 4% formaldehyde (FA) in PB for 10 min at the speed of 6 ml/min. Following perfusion, the brains were extracted and post-fixed by immersion in the fixative at 4°C overnight protected from light. The tissue was stored in 0.1% FA in PB at 4°C protected from light until further processing. Cryoprotective steps in graded sucrose solution (10, 20, and 30% sucrose in PB) were performed before the brains were cut on a freeze microtome (Thermo Scientific Microm KS 34) in 40 µm free-floating coronary sections and stored in 0.1% FA in PB at 4°C protected from light until usage. In addition, eight tg-ArcSwe animals and eight WT littermates were sacrificed at 12 months of age and used for qPCR analysis. These animals were anesthetized as described above and decapitated. The brains were extracted, and the left hemisphere was dissected into the frontal cortex, hippocampus, cerebellum, and the rest of the brain. This tissue was frozen and stored in –80°C pending analysis.

### RNA isolation and real-time PCR

48 hr prior to RNA extraction, the samples were suspended in RNAlater-ICE (Ambion; Cat# AM7030). To isolate total RNA from the frontal cortex tissue samples, the RNeasy Mini Kit (QIAGEN, Hilden, Germany), including the on-column DNase digestion, was used. The RNA concentration and integrity were determined using a NanoDrop 2000c spectrophotometer (Thermo Fisher Scientific) and ethidium bromide visualization after agarose gel electrophoresis, respectively. Following the manufacturer’s protocol, 1 µg of total RNA was reverse-transcribed into cDNA with Oligo (dT)15 using the GoScript Reverse Transcription System (Promega, Madison, USA, Cat# A5001). All the cDNA samples were diluted in Tris-EDTA buffer (pH 8.0) to a final concentration of 2.5 ng/μl. Real-time PCR was carried out in a total volume of 20 µl, containing 2× AB Power SYBR Green PCR Master Mix (Thermo Fisher Scientific) with gene-specific primers (at a final concentration of 200  nM) and 2 µl cDNA samples. Amplification was performed on the StepOnePlus system (Applied Biosystems) with the following conditions: 95°C for 10 min, followed by 40 cycles of 95°C for 15 s and 60°C for 1 min, followed by melting curve analysis to check for unspecific products. Each sample was run in duplicates. Using the NormFinder software 4 Andersen et al.2004, HPRT1 was determined as an internal control for normalization of the gene expression. The primers were designed online using Primer BLAST and setting the amplicon size to a maximum of 200  bp. The primers designed span exon–exon junctions, and standards prepared as previously described 51 Rao et al.2021. Details of the Gfap forward and reverse primer are presented in Table 7.

Table 7.

#### Primer used for mRNA quantification.

Gene | Protein name | Accession | Forward (5′) | Reverse (3′)
--- | --- | --- | --- | ---
Gfap | Glial fibrillary acidic protein | NM_010277.3 | GCACTCAATACGAGGCAGTG | GCTCTAGGGACTCGTTCGTG

### Immunohistochemistry

One section from each animal was chosen for quantification of astrogliosis and washed in PBS 0.01 M for 10 min, followed by two times in 0.1% TritonX100 in PBS (PBST) for 5 min. The PBST was removed and blocking (10% normal donkey serum [NDS], 1% bovine serum albumin [BSA], 0.5% Triton X100 in PBS) performed for 1 hr at room temperature. This was directly followed by incubation overnight at room temperature with primary antibodies (GFAP; host: mouse; diluted 1:1000; Sigma-Aldrich; Cat# MAB360. GFP; host: chicken; diluted 1:2000; Abcam; Cat# ab13970) diluted in antibody solution (ABS; 3% NDS, 1% BSA, 0.1% Triton X100 in PBS). The following day the sections were rinsed in 0.1% PBST two times for 1 min, followed by three times for 10 min. Secondary antibodies (CY5 donkey anti-mouse; Jackson ImmunoResearch Labs; Cat# 715-175-151; CY3 donkey anti-chicken; Jackson ImmunoResearch Labs; Cat# 703-165-155) were spun in a centrifuge for 10 min at 13,000 rpm, diluted 1:500 in ABS, and the sections incubated for 1 hr at room temperature. After the second incubation, the sections were washed in PBS for 10 min three times. Propidium iodide (diluted 1:5000 in 0.01 M PBS; Sigma-Aldrich; Cat# 04511 [Cellstain double staining kit]) for nuclear staining was added for 10 min, before rinsing the sections twice for 5 min in PBS. All sections were transferred to distilled water and mounted with ProLong Gold antifade reagent (Thermo Fisher Scientific; Cat# P36934). They were stored in –20°C protected from light until confocal imaging. For electron microscopy, EM standard procedure was followed for embedding and preparing of the tissue 73 Yang et al.2011. For enhancing the contrast, uranyl acetate (Fluorochem) in double-distilled water and lead citrate was used. The sections were examined in a transmission electron microscope (TECNAI 12).

### Confocal imaging

GFP-positive astrocytes (indicating GCaMP6f expression) were used to locate the appropriate cortical area used for in vivo imaging. Only sections where we could locate positive GFP staining within the cortex corresponding to the image area were chosen for confocal imaging and further analysis. All single-plane and z-stack images were acquired using a Zeiss LSM 710 confocal microscope. To provide overview and verify that the correct area was identified, one tile scan of 3072 × 3072 pixels was achieved with a ×40 objective (1.20; water korr M27) using three channels (CY2, CY3, and CY5) for WTs and four channels (Dapi, CY2, CY3, and CY5) for tg-ArcSwe. Next, one z-stack of 2048 × 2048 pixels (×40 objective; 1.20; water korr M27) from within the GFP-positive area was acquired from cortical layer 2. All single-plane and z-stack images were obtained with identical settings. In addition, a z-stack of 1024 × 1024 pixels was obtained (×40 objective; 1.20; water korr M27) outside the located GFP-positive stained area. No postprocessing on the analyzed images was performed.

### 3D reconstruction analysis of astrocytes

This procedure was adapted from the detail outlined in 64 Tavares et al.2017 that utilizes the free software plugin Simple Neurite Tracer (SNT) of Fiji ImageJ 39 Longair et al.2011. Only astrocytes with a single nucleus where at least ⅔ of the circumference was covered by GFAP staining were selected for 3D reconstruction. Astrocytes with processes touching the borders of FOV were omitted. Eight randomly selected astrocytes from each z-stack (two from each image quadrant) obtained within the GFP-positive imaging area were used for analysis. To quantify and visualize the morphological complexity of the astrocytes, we analyzed the number of processes, total length of processes in µm, process thickness (µm3), and number of intersections (provided from the Sholl analysis).

### GFAP and GFP area fraction analyses

Z-projections based on the CY5 channel to visualize GFAP/GFP-positive astrocytes of all stack images were rendered using the Fiji ImageJ software 56 Schindelin et al.2012. The images were blinded to the analyst, converted to 8-bit images and the scale removed. The threshold was manually adjusted so that only what was considered to be GFAP-positive staining was red, and a percentage value of positive staining in the image was obtained.

### Statistical analyses

Astrocytic Ca2+ signals were studied by means of the ROA density – the fraction of the compartment with activity per time unit. Here, the area may be the entire FOV or we may limit ourselves to the area identified as belonging to cellular subcompartments; the astrocytic processes, somata, or endfeet. In each trial, we had time series of ROA density lasting around 300 s (approximately 9000 frames). Within each trial, the startle period was defined as starting from the air puff at 150 s and lasting 600 frames (~20 s). In addition, one or more time periods within the trial could be identified as spontaneous runs or as periods of quiet wakefulness. A first important question concerns the choice of summary statistics adequately describing the Ca2+ response in such periods of interest (runs and startle), which are dynamic behavioral states that entail both acceleration, steady locomotion and deceleration. We have studied the mean and max ROA density and the ROA density rise rate, which is defined as the maximal increase in ROA density over a maximum of 50 frames. Initial explorations indicated that the main results are fairly robust to the choice of window length, and 50 frames appeared to be a sensible choice compared to the kinetics of astrocytic Ca2+ signals. The ROA density rise rate is meant to capture some of the dynamics in astrocytic Ca2+ signaling and can be understood as the maximum acceleration inside the time period of interest. The rise rate has a high correlation with the maximum, and if the length of the window is increased sufficiently these two statistics tend to become almost identical (since most traces are close to zero at some point in the trial). See Figure 4—figure supplement 3A, which displays both the max ROA density and the rise rate in an example.

Pupil size measurements had a coarser time resolution than for the astrocytic Ca2+ signals; see example in Figure 4—figure supplement 3B, showing the pupil sizes around a startle response. The pattern in the figure – a sharp increase in pupil size after air puff, before a gradual decline – was quite typical. Therefore, we chose to only consider the measurements in a time window of 6.67 s on each side of the air puff (or start of running). We defined the pupil dilation as the relative increase in the ratio of pupil diameter to eye diameter inside this window. In other words, we calculate the average ratio before the air puff (or start of running) and after the air puff, compute the difference, and divide by the average before the air puff.

For GRABNE analyses, time series were segmented similar to the Ca2+ imaging trials. Due to the slower response, change in GRABNE fluorescence was calculated slightly different than the increases in Ca2+ signaling as the difference between the median GRABNE ΔF/F0 in the baseline period (10 s immediately prior to locomotion/startle) and during locomotion/starle. Pupil responses and locomotor responses were defined within the same time intervals for these trials.

#### Interpreting hierarchical plots

We have chosen to present some of our data in the form of hierarchical plots (see Figure 3B and C, Figure 4B and C, Figure 5B, and figure supplements). These plots allow the reader to assess the degree of separation between the genotypes and at the same time get an impression of the variation at different levels of the analysis – in our case, the variation between different mice of the same genotype and the variation between repeated measurements on the same mouse. At the lowest level – the trial level – we have points representing the observations themselves. For the spontaneous runs, there are sometimes more than one run per trial and in that case the points are the median ROA density rise rate in these runs. The maximal number of spontaneous runs in a single trial was 5 (average: 1.6 runs per trial). At the middle level, we have the median ROA density rise rate for each mouse, and at the top level the median ROA density rise rate for each genotype. The lines between the levels indicate which observations belong to each mouse and to each genotype respectively. We have made similar plots for the max ROA density also. The hierarchical plots are a useful tool for exploratory analysis. They are also meant to promote transparency in scientific reporting and highlight the importance of intra-group variation. Still, it is important to realize that the impression conveyed by the plots might not be identical to the results from statistical modeling. In the plots, we do not include the influence of various technical and biological covariates that one typically would include in a statistical model. Some of the variation between observations belonging to the same mouse and between different mice of the same genotype may be explained by such covariates, as we will see in the next section.

#### Modeling

Statistical analyses were conducted in R (version 4.1.1). The ROA density rise rate was modeled by linear mixed effect regression models, which were fitted using the glmmTMB package 11 Brooks et al.2017. We conducted two sets of analyses: (1) to investigate potential differences in ROA density rise rate between the two genotypes and (2) to investigate the relationship between ROA density rise rate and the pupil dilation, including potential differences between the two genotypes with respect to this relationship. For (1), the coefficient of primary interest is the one belonging to the genotype variable, while for (2) we are interested in the effect of pupil dilation on the ROA density rise rate, as well as the interaction between this pupil effect and the genotype. In both of these sets of analyses, we adjusted for the following fixed effect covariates: the level of magnification (μm per pixel) (two levels), the depth of the measurements (in μm), and the maximal speed in the relevant time window. We included random intercepts for each mouse (five WT and six tg-ArcSwe). We analyzed the ROA density max and mean values with similar models as the ones described here for the rise rate. For (1), the number of observations ranged between 77 and 117 trials with startle data (depending on the subcompartment) and between 72 and 109 episodes of spontaneous running, while for (2) the number of observations ranged between 60 and 86 for the startle data and between 35 and 44 for the spontaneous runs. There were less observations for the (2) analyses because some episodes/trials had missing or incomplete pupil measurements. Similar modeling was performed for GRABNE analyses.

The sensitivity of our result to these modeling choices was assessed by various robustness checks (see next section). The adequacy of model assumptions was investigated by residual plots 24 Hartig 2022. In the cases where the residual plots indicated deviations from the assumption of constant residual variance, we extended the model by allowing the residual variance to vary as a function of genotype. The reported p-values are based on the t-distribution, with degrees of freedom as provided from the glmmTMB package. No corrections for multiple comparisons were applied.

When analyzing the relationship between distance from the nearest plaque and astrocytic Ca2+ signaling (Figure 2E), we considered the mean ROA density in each ROI (n = 10,988) in the quiet wakefulness episodes. If present, any effect of plaque proximity on Ca2+ signaling should be discernible among ROIs observed in the same episode. The dashed lines in Figure 2E show the effect of distance on Ca2+ signaling within each episode, and they form an uncertain picture: in some episodes, there is a weak positive relationship, with seemingly higher mean ROA density further away from plaques, while in many episodes there is a negative relationship, with somewhat higher mean ROA density close to plaques. The overall line is found by fitting a linear mixed effect model with mean ROA density as the response, with a fixed effect of distance and with each episode having its own random intercept and slope for the distance effect.

For the results presented in Figure 1, unless otherwise stated, the data are presented as mean ± standard error of the mean (SEM). A p-value equal to or less than 0.05 was considered statistically significant. Mann–Whitney U-test was used to analyze the number of processes, total length of processes in µm, process thickness in µm3, and area fraction of positive GFAP staining. Two-way ANOVA followed by Sidak post-hoc comparison was used to analyze the number of intersections. Statistical analysis was performed in GraphPad Prism version 8.0.1 for Windows (GraphPad Software). For qPCR analysis, mean copy number per ng of total RNA was compared between genotypes by Mann–Whitney U-test in SPSS Statistics 26 (SPSS).

#### Robustness checks

In order to check the robustness of the various statistical analyses, the stability of estimates and p-values was examined with respect to the length of window, for the rise rate response variable, and also sensitivity to individual mice and trials.

The type of sensitivity analysis performed is illustrated (see Figure 2—figure supplement 2) for the analysis of the uncoupling between pupil dilation and astrocytic Ca2+ responses (see Figures 3 and 4 for details about the data). Again, for the purpose of illustration, we focus on the interaction between pupil dilation and genotype for the ROA density rise rate response variable; see Figure 2—figure supplement 2 for details. Similar analyses were performed for the main statistical models; none of these sensitivity analyses provided clear or strong evidence for any change in the main conclusions.

### Data and source code availability

A complete dataset of raw and processed data including two-photon microscopy traces, behavioral monitoring, and pupil tracking data is provided at https://doi.org/10.11582/2021.00100. Data analyses were performed with the Begonia toolkit 7 Bjørnstad et al.2021, which is available at https://github.com/GliaLab/Begonia (copy archived at swh:1:rev:30447e22ca87fe7a2308b0029240bc4491a74518, 8 Bjørnstad 2022).

## References

1. 1. AY Abramov
   2. L Canevari
   3. MR Duchen
   2003 Changes in intracellular calcium and glutathione in astrocytes as the primary mechanism of amyloid neurotoxicity 23 The Journal of Neuroscience 5088 5095
2. 1. AY Abramov
   2. L Canevari
   3. MR Duchen
   2004 Calcium signals induced by amyloid beta peptide and their consequences in neurons and astrocytes in culture 1742 Biochimica et Biophysica Acta 81 87
3. 1. A Adamsky
   2. A Kol
   3. T Kreisel
   4. A Doron
   5. N Ozeri-Engelhard
   6. T Melcer
   7. R Refaeli
   8. H Horn
   9. L Regev
   10. M Groysman
   11. M London
   12. I Goshen
   2018 Astrocytic activation generates de novo neuronal potentiation and memory enhancement 174 Cell 59 71
4. 1. CL Andersen
   2. JL Jensen
   3. TF Ørntoft
   2004 Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets 64 Cancer Research 5245 5250
5. 1. LK Bekar
   2. W He
   3. M Nedergaard
   2008 Locus coeruleus alpha-adrenergic-mediated activation of cortical astrocytes in vivo 18 Cerebral Cortex 2789 2795
6. 1. E Bindocci
   2. I Savtchouk
   3. N Liaudet
   4. D Becker
   5. G Carriero
   6. A Volterra
   2017 Three-dimensional Ca2+ imaging advances understanding of astrocyte biology 356 Science
7. 1. DM Bjørnstad
   2. KS Åbjørsbråten
   3. E Hennestad
   4. C Cunen
   5. GH Hermansen
   6. L Bojarskaite
   7. KH Pettersen
   8. K Vervaeke
   9. R Enger
   2021 Begonia—a two-photon imaging analysis pipeline for astrocytic ca2+ signals 15 Frontiers in Cellular Neuroscience
8. 1. DM Bjørnstad
   2022 Begonia https://archive.softwareheritage.org/swh:1:dir:2e8c4af17a6dc52b248460ed52bd9bc28f389f09;origin=https://github.com/GliaLab/Begonia;visit=swh:1:snp:853be548f2e884d045707f4eef209bc70a60b2ce;anchor=swh:1:rev:30447e22ca87fe7a2308b0029240bc4491a74518
9. 1. L Bojarskaite
   2. DM Bjørnstad
   3. KH Pettersen
   4. C Cunen
   5. GH Hermansen
   6. KS Åbjørsbråten
   7. AR Chambers
   8. R Sprengel
   9. K Vervaeke
   10. W Tang
   11. R Enger
   12. EA Nagelhus
   2020 Astrocytic Ca2+ signaling is reduced during sleep and is involved in the regulation of slow wave sleep 11 Nature Communications 1 16
10. 1. H Braak
    2. K Del Tredici
    2011 The pathological process underlying Alzheimer’s disease in individuals under thirty 121 Acta Neuropathologica 171 181
11. 1. M Brooks
    2. K Kristensen
    3. K Benthem
    4. A Magnusson
    5. C Berg
    6. A Nielsen
    7. H Skaug
    8. M Mächler
    9. B Bolker
    2017 Glmmtmb balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling 9 The R Journal
12. 1. V Caggiano
    2. R Leiras
    3. H Goñi-Erro
    4. D Masini
    5. C Bellardita
    6. J Bouvier
    7. V Caldeira
    8. G Fisone
    9. O Kiehn
    2018 Midbrain circuits that set locomotor speed and gait selection 553 Nature 455 460
13. 1. T-W Chen
    2. TJ Wardill
    3. Y Sun
    4. SR Pulver
    5. SL Renninger
    6. A Baohan
    7. ER Schreiter
    8. RA Kerr
    9. MB Orger
    10. V Jayaraman
    11. LL Looger
    12. K Svoboda
    13. DS Kim
    2013 Ultrasensitive fluorescent proteins for imaging neuronal activity 499 Nature 295 300
14. 1. A Codita
    2. A Gumucio
    3. L Lannfelt
    4. P Gellerfors
    5. B Winblad
    6. AH Mohammed
    7. LNG Nilsson
    2010 Impaired behavior of female tg-ArcSwe APP mice in the IntelliCage: A longitudinal study 215 Behavioural Brain Research 83 94
15. 1. VD Costa
    2. PH Rudebeck
    2016 More than meets the eye: The relationship between pupil size and locus coeruleus activity 89 Neuron 8 10
16. 1. A Delekate
    2. M Füchtemeier
    3. T Schumacher
    4. C Ulbrich
    5. M Foddis
    6. GC Petzold
    2014 Metabotropic P2Y1 receptor signalling mediates astrocytic hyperactivity in vivo in an Alzheimer’s disease mouse model 5 Nature Communications
17. 1. F Ding
    2. J O’Donnell
    3. AS Thrane
    4. D Zeppenfeld
    5. H Kang
    6. L Xie
    7. F Wang
    8. M Nedergaard
    2013 α1-Adrenergic receptors mediate coordinated Ca2+ signaling of cortical astrocytes in awake, behaving mice 54 Cell Calcium 387 394
18. 1. C Escartin
    2. E Galea
    3. A Lakatos
    4. JP O’Callaghan
    5. GC Petzold
    6. A Serrano-Pozo
    7. C Steinhäuser
    8. A Volterra
    9. G Carmignoto
    10. A Agarwal
    11. NJ Allen
    12. A Araque
    13. L Barbeito
    14. A Barzilai
    15. DE Bergles
    16. G Bonvento
    17. AM Butt
    18. W-T Chen
    19. M Cohen-Salmon
    20. C Cunningham
    21. B Deneen
    22. B De Strooper
    23. B Díaz-Castro
    24. C Farina
    25. M Freeman
    26. V Gallo
    27. JE Goldman
    28. SA Goldman
    29. M Götz
    30. A Gutiérrez
    31. PG Haydon
    32. DH Heiland
    33. EM Hol
    34. MG Holt
    35. M Iino
    36. KV Kastanenka
    37. H Kettenmann
    38. BS Khakh
    39. S Koizumi
    40. CJ Lee
    41. SA Liddelow
    42. BA MacVicar
    43. P Magistretti
    44. A Messing
    45. A Mishra
    46. AV Molofsky
    47. KK Murai
    48. CM Norris
    49. S Okada
    50. SHR Oliet
    51. JF Oliveira
    52. A Panatier
    53. V Parpura
    54. M Pekna
    55. M Pekny
    56. L Pellerin
    57. G Perea
    58. BG Pérez-Nievas
    59. FW Pfrieger
    60. KE Poskanzer
    61. FJ Quintana
    62. RM Ransohoff
    63. M Riquelme-Perez
    64. S Robel
    65. CR Rose
    66. JD Rothstein
    67. N Rouach
    68. DH Rowitch
    69. A Semyanov
    70. S Sirko
    71. H Sontheimer
    72. RA Swanson
    73. J Vitorica
    74. I-B Wanner
    75. LB Wood
    76. J Wu
    77. B Zheng
    78. ER Zimmer
    79. R Zorec
    80. MV Sofroniew
    81. A Verkhratsky
    2021 Reactive astrocyte nomenclature, definitions, and future directions 24 Nature Neuroscience 312 325
19. 1. J Feng
    2. C Zhang
    3. JE Lischinsky
    4. M Jing
    5. J Zhou
    6. H Wang
    7. Y Zhang
    8. A Dong
    9. Z Wu
    10. H Wu
    11. W Chen
    12. P Zhang
    13. J Zou
    14. SA Hires
    15. JJ Zhu
    16. G Cui
    17. D Lin
    18. J Du
    19. Y Li
    2019 A genetically encoded fluorescent sensor for rapid and specific in vivo detection of norepinephrineA Genetically Encoded Fluorescent Sensor for Rapid and Specific In Vivo Detection of Norepinephrine 102 Neuron 745 761
20. 1. MJ Ferreira-Pinto
    2. L Ruder
    3. P Capelli
    4. S Arber
    2018 Connecting circuits for supraspinal control of locomotion 100 Neuron 361 374
21. 1. S Grillner
    2. A El Manira
    2020 Current principles of motor control, with special reference to vertebrate locomotion 100 Physiological Reviews 271 320
22. 1. N Habib
    2. C McCabe
    3. S Medina
    4. M Varshavsky
    5. D Kitsberg
    6. R Dvir-Szternfeld
    7. G Green
    8. D Dionne
    9. L Nguyen
    10. JL Marshall
    11. F Chen
    12. F Zhang
    13. T Kaplan
    14. A Regev
    15. M Schwartz
    2020 Disease-associated astrocytes in Alzheimer’s disease and aging 23 Nature Neuroscience 701 706
23. 1. MB Hancock
    2. CL Fougerousse
    1976 Spinal projections from the nucleus locus coeruleus and nucleus subcoeruleus in the cat and monkey as demonstrated by the retrograde transport of horseradish peroxidase 1 Brain Research Bulletin 229 234
24. 1. F Hartig
    2022 Residual diagnostics for hierarchical models http://florianhartig.github.io/DHARMa/
25. 1. NJ Haughey
    2. MP Mattson
    2003 Alzheimer’s amyloid β-peptide enhances atp/gap junction-mediated calcium-wave propagation in astrocytes 3 NeuroMolecular Medicine 173 180
26. 1. PG Haydon
    2. M Nedergaard
    2015 How do astrocytes participate in neural plasticity?7 Cold Spring Harbor Perspectives in Biology
27. 1. J Hirrlinger
    2. A Scheller
    3. PG Hirrlinger
    4. B Kellert
    5. W Tang
    6. MC Wehr
    7. S Goebbels
    8. A Reichenbach
    9. R Sprengel
    10. MJ Rossner
    11. F Kirchhoff
    12. DP Martin
    2009 Split-cre complementation indicates coincident activity of different genes in vivo 4 PLOS ONE
28. 1. N Holland
    2. TW Robbins
    3. JB Rowe
    2021 The role of noradrenaline in cognition and cognitive disorders 144 Brain 2243 2256
29. 1. HIL Jacobs
    2. JA Becker
    3. K Kwong
    4. N Engels-Domínguez
    5. PC Prokopiou
    6. KV Papp
    7. M Properzi
    8. OL Hampton
    9. F Uquillas
    10. JS Sanchez
    11. DM Rentz
    12. G El Fakhri
    13. MD Normandin
    14. JC Price
    15. DA Bennett
    16. RA Sperling
    17. KA Johnson
    2021 In vivo and neuropathology data support locus coeruleus integrity as indicator of Alzheimer’s disease pathology and cognitive decline 13 Science Translational Medicine
30. 1. RN Kalaria
    2. AC Andorn
    3. M Tabaton
    4. PJ Whitehouse
    5. SI Harik
    6. JR Unnerstall
    1989 Adrenergic receptors in aging and Alzheimer’s disease: increased beta 2-receptors in prefrontal cortex and hippocampus 53 Journal of Neurochemistry 1772 1781
31. 1. RN Kalaria
    2. SI Harik
    1989 Increased alpha 2- and beta 2-adrenergic receptors in cerebral microvessels in Alzheimer disease 106 Neuroscience Letters 233 238
32. 1. C Kjaerby
    2. R Rasmussen
    3. M Andersen
    4. M Nedergaard
    2017 Does global astrocytic calcium signaling participate in awake brain state transitions and neuronal circuit function?42 Neurochemical Research 1810 1822
33. 1. WE Klunk
    2. BJ Bacskai
    3. CA Mathis
    4. ST Kajdasz
    5. ME McLellan
    6. MP Frosch
    7. ML Debnath
    8. DP Holt
    9. Y Wang
    10. BT Hyman
    2002 Imaging Abeta plaques in living transgenic mice with multiphoton microscopy and methoxy-X04, A systemically administered Congo red derivative 61 Journal of Neuropathology and Experimental Neurology 797 805
34. 1. A Kol
    2. A Adamsky
    3. M Groysman
    4. T Kreisel
    5. M London
    6. I Goshen
    2020 Astrocytes contribute to remote memory formation by modulating hippocampal-cortical communication during learning 23 Nature Neuroscience 1229 1239
35. 1. KV Kuchibhotla
    2. CR Lattarulo
    3. BT Hyman
    4. BJ Bacskai
    2009 Synchronous hyperactivity and intercellular calcium waves in astrocytes in Alzheimer mice 323 Science 1211 1215
36. 1. S Lillehaug
    2. GH Syverstad
    3. LNG Nilsson
    4. JG Bjaalie
    5. TB Leergaard
    6. R Torp
    2014 Brainwide distribution and variance of amyloid-beta deposits in tg-ArcSwe mice 35 Neurobiology of Aging 556 564
37. 1. D Lim
    2. A Iyer
    3. V Ronco
    4. AA Grolla
    5. PL Canonico
    6. E Aronica
    7. AA Genazzani
    2013 Amyloid beta deregulates astroglial mGluR5-mediated calcium signaling via calcineurin and Nf-kB 61 Glia 1134 1145
38. 1. Y Liu
    2. C Rodenkirch
    3. N Moskowitz
    4. B Schriver
    5. Q Wang
    2017 Dynamic lateralization of pupil dilation evoked by locus coeruleus activation results from sympathetic, not parasympathetic, contributions 20 Cell Reports 3099 3112
39. 1. MH Longair
    2. DA Baker
    3. JD Armstrong
    2011 Simple Neurite Tracer: open source software for reconstruction, visualization and analysis of neuronal processes 27 Bioinformatics 2453 2454
40. 1. A Lord
    2. H Kalimo
    3. C Eckman
    4. X-Q Zhang
    5. L Lannfelt
    6. LNG Nilsson
    2006 The Arctic Alzheimer mutation facilitates early intraneuronal Abeta aggregation and senile plaque formation in transgenic mice 27 Neurobiology of Aging 67 77
41. 1. BJ Matchett
    2. LT Grinberg
    3. P Theofilas
    4. ME Murray
    2021 The mechanistic link between selective vulnerability of the locus coeruleus and neurodegeneration in Alzheimer’s disease 141 Acta Neuropathologica 631 650
42. 1. MF McCool
    2. GB Varty
    3. RA Del Vecchio
    4. TM Kazdoba
    5. EM Parker
    6. JC Hunter
    7. LA Hyde
    2003 Increased auditory startle response and reduced prepulse inhibition of startle in transgenic mice expressing a double mutant form of amyloid precursor protein 994 Brain Research 99 106
43. 1. M Meyer-Luehmann
    2. TL Spires-Jones
    3. C Prada
    4. M Garcia-Alloza
    5. A Calignon
    6. A Rozkalne
    7. J Koenigsknecht-Talboo
    8. DM Holtzman
    9. BJ Bacskai
    10. BT Hyman
    2008 Rapid appearance and local toxicity of amyloid-beta plaques in a mouse model of Alzheimer’s disease 451 Nature 720 724
44. 1. J O’Donnell
    2. D Zeppenfeld
    3. E McConnell
    4. S Pena
    5. M Nedergaard
    2012 Norepinephrine: A neuromodulator that boosts the function of multiple cell types to optimize CNS performance 37 Neurochemical Research 2496 2512
45. 1. M Paukert
    2. A Agarwal
    3. J Cha
    4. VA Doze
    5. JU Kang
    6. DE Bergles
    2014 Norepinephrine controls astroglial responsiveness to local circuit activity 82 Neuron 1263 1270
46. 1. AC Peterson
    2. C-SR Li
    2018 Noradrenergic dysfunction in alzheimer’s and parkinson’s diseases-an overview of imaging studies 10 Frontiers in Aging Neuroscience
47. 1. O Philipson
    2. P Hammarström
    3. KPR Nilsson
    4. E Portelius
    5. T Olofsson
    6. M Ingelsson
    7. BT Hyman
    8. K Blennow
    9. L Lannfelt
    10. H Kalimo
    11. LNG Nilsson
    2009 A highly insoluble state of Abeta similar to that of Alzheimer’s disease brain is found in Arctic APP transgenic mice 30 Neurobiology of Aging 1393 1405
48. 1. EA Pnevmatikakis
    2. A Giovannucci
    2017 NoRMCorre: An online algorithm for piecewise rigid motion correction of calcium imaging data 291 Journal of Neuroscience Methods 83 94
49. 1. KE Poskanzer
    2. R Yuste
    2011 Astrocytic regulation of cortical UP states 108 PNAS 18453 18458
50. 1. KE Poskanzer
    2. R Yuste
    2016 Astrocytes regulate cortical state switching in vivo 113 PNAS
51. 1. SB Rao
    2. N Skauli
    3. N Jovanovic
    4. S Katoozi
    5. A Frigeri
    6. SC Froehner
    7. ME Adams
    8. OP Ottersen
    9. M Amiry-Moghaddam
    2021 Orchestrating aquaporin-4 and connexin-43 expression in brain: Differential roles of α1- and β1-syntrophin 8 Biochim Biophys Acta Biomembr
52. 1. N Reichenbach
    2. A Delekate
    3. B Breithausen
    4. K Keppler
    5. S Poll
    6. T Schulte
    7. J Peter
    8. M Plescher
    9. JN Hansen
    10. N Blank
    11. A Keller
    12. M Fuhrmann
    13. C Henneberger
    14. A Halle
    15. GC Petzold
    2018 P2Y1 receptor blockade normalizes network dysfunction and cognition in an Alzheimer’s disease model 215 The Journal of Experimental Medicine 1649 1663
53. 1. J Reimer
    2. MJ McGinley
    3. Y Liu
    4. C Rodenkirch
    5. Q Wang
    6. DA McCormick
    7. AS Tolias
    2016 Pupil fluctuations track rapid changes in adrenergic and cholinergic activity in cortex 7 Nature Communications
54. 1. JJ Rodríguez-Arellano
    2. V Parpura
    3. R Zorec
    4. A Verkhratsky
    2016 Astrocytes in physiological aging and Alzheimer’s disease 323 Neuroscience 170 182
55. 1. F Sano
    2. E Shigetomi
    3. Y Shinozaki
    4. H Tsuzukiyama
    5. K Saito
    6. K Mikoshiba
    7. K Sugita
    8. M Aihara
    9. S Koizumi
    2019 Reactive Astrocyte-Driven Epileptogenesis Is Induced by Microglia Initially Activated Following Status Epilepticus bioRxiv
56. 1. J Schindelin
    2. I Arganda-Carreras
    3. E Frise
    4. V Kaynig
    5. M Longair
    6. T Pietzsch
    7. S Preibisch
    8. C Rueden
    9. S Saalfeld
    10. B Schmid
    11. J-Y Tinevez
    12. DJ White
    13. V Hartenstein
    14. K Eliceiri
    15. P Tomancak
    16. A Cardona
    2012 Fiji: an open-source platform for biological-image analysis 9 Nature Methods 676 682
57. 1. Z Shevtsova
    2. JMI Malik
    3. U Michel
    4. M Bähr
    5. S Kügler
    2005 Promoters and serotypes: targeting of adeno-associated virus vectors for gene transfer in the rat central nervous system in vitro and in vivo 90 Experimental Physiology 53 59
58. 1. E Shigetomi
    2. K Saito
    3. F Sano
    4. S Koizumi
    2019 Aberrant calcium signals in reactive astrocytes: A key process in neurological disorders 20 International Journal of Molecular Sciences
59. 1. GHS Skaaraas
    2. C Melbye
    3. MA Puchades
    4. DSY Leung
    5. D Jacobsen
    6. SB Rao
    7. OP Ottersen
    8. TB Leergaard
    9. R Torp
    2021 Cerebral amyloid angiopathy in a mouse model of alzheimer’s disease associates with upregulated angiopoietin and downregulated hypoxia-inducible factor 83 Journal of Alzheimer’s Disease 1651 1663
60. 1. RH Smith
    2. JR Levy
    3. RM Kotin
    2009 A simplified baculovirus-AAV expression vector system coupled with one-step affinity purification yields high-titer rAAV stocks from insect cells 17 Molecular Therapy 1888 1896
61. 1. R Srinivasan
    2. BS Huang
    3. S Venugopal
    4. AD Johnston
    5. H Chai
    6. H Zeng
    7. P Golshani
    8. BS Khakh
    2015 Ca(2+) signaling in astrocytes from Ip3r2(-/-) mice in brain slices and during startle responses in vivo 18 Nature Neuroscience 708 717
62. 1. JL Stobart
    2. KD Ferrari
    3. MJP Barrett
    4. C Glück
    5. MJ Stobart
    6. M Zuend
    7. B Weber
    2018 Cortical circuit activity evokes rapid astrocyte calcium signals on a similar timescale to neurons 98 Neuron 726 735
63. 1. W Tang
    2. K Szokol
    3. V Jensen
    4. R Enger
    5. CA Trivedi
    6. Ø Hvalby
    7. PJ Helm
    8. LL Looger
    9. R Sprengel
    10. EA Nagelhus
    2015 Stimulation-evoked Ca2+ signals in astrocytic processes at hippocampal CA3-CA1 synapses of adult mice are modulated by glutamate and ATP 35 The Journal of Neuroscience 3016 3021
64. 1. G Tavares
    2. M Martins
    3. JS Correia
    4. VM Sardinha
    5. S Guerra-Gomes
    6. SP Neves
    7. F Marques
    8. N Sousa
    9. JF Oliveira
    2017 Employing an open-source tool to assess astrocyte tridimensional structure 222 Brain Structure & Function 1989 1999
65. 1. AS Thrane
    2. V Rangroo Thrane
    3. D Zeppenfeld
    4. N Lou
    5. Q Xu
    6. EA Nagelhus
    7. M Nedergaard
    2012 General anesthesia selectively disrupts astrocyte calcium signaling in the awake mouse cortex 109 PNAS 18974 18979
66. 1. TV Vaidyanathan
    2. M Collard
    3. S Yokoyama
    4. ME Reitman
    5. KE Poskanzer
    2021 Cortical astrocytes independently regulate sleep depth and duration via separate GPCR pathways 10 eLife
67. 1. A Verkhratsky
    2019 Astroglial calcium signaling in aging and alzheimer’s disease 11 Cold Spring Harbor Perspectives in Biology
68. 1. J Wahis
    2. MG Holt
    2021 Astrocytes, noradrenaline, α1-adrenoreceptors, and neuromodulation: Evidence and unanswered questions 15 Frontiers in Cellular Neuroscience
69. 1. F Wang
    2. NA Smith
    3. Q Xu
    4. T Fujita
    5. A Baba
    6. T Matsuda
    7. T Takano
    8. L Bekar
    9. M Nedergaard
    2012 Astrocytes modulate neural network activity by Ca 5 Science Signaling
70. 1. Y Wang
    2. NV DelRosso
    3. TV Vaidyanathan
    4. MK Cahill
    5. ME Reitman
    6. S Pittolo
    7. X Mi
    8. G Yu
    9. KE Poskanzer
    2019 Accurate quantification of astrocyte and neurotransmitter fluorescence dynamics for single-cell and population-level physiology 22 Nature Neuroscience 1936 1944
71. 1. D Weinshenker
    2018 Long road to ruin: Noradrenergic dysfunction in neurodegenerative disease 41 Trends in Neurosciences 211 223
72. 1. U Wilhelmsson
    2. EA Bushong
    3. DL Price
    4. BL Smarr
    5. V Phung
    6. M Terada
    7. MH Ellisman
    8. M Pekny
    2006 Redefining the concept of reactive astrocytes as cells that remain within their unique domains upon reaction to injury 103 PNAS 17513 17518
73. 1. J Yang
    2. LK Lunde
    3. P Nuntagij
    4. T Oguchi
    5. LMA Camassa
    6. LNG Nilsson
    7. L Lannfelt
    8. Y Xu
    9. M Amiry-Moghaddam
    10. OP Ottersen
    11. R Torp
    2011 Loss of astrocyte polarization in the tg-ArcSwe mouse model of Alzheimer’s disease 27 Journal of Alzheimer’s Disease 711 722
74. 1. L Ye
    2. M Orynbayev
    3. X Zhu
    4. EY Lim
    5. RR Dereddi
    6. A Agarwal
    7. DE Bergles
    8. MA Bhat
    9. M Paukert
    2020 Ethanol abolishes vigilance-dependent astroglia network activation in mice by inhibiting norepinephrine release 11 Nature Communications 1 20
75. 1. M Zuend
    2. AS Saab
    3. MT Wyss
    4. KD Ferrari
    5. L Hösli
    6. ZJ Looser
    7. JL Stobart
    8. J Duran
    9. JJ Guinovart
    10. LF Barros
    11. B Weber
    2020 Arousal-induced cortical activity triggers lactate release from astrocytes 2 Nature Metabolism 179 191