Network Models Predict That Pyramidal Neuron Hyperexcitability and Synapse Loss in the dlPFC Lead to Age-Related Spatial Working Memory Impairment in Rhesus Monkeys

Behavioral studies have shown spatial working memory impairment with aging in several animal species, including humans. Persistent activity of layer 3 pyramidal dorsolateral prefrontal cortex (dlPFC) neurons during delay periods of working memory tasks is important for encoding memory of the stimulus. In vitro studies have shown that these neurons undergo significant age-related structural and functional changes, but the extent to which these changes affect neural mechanisms underlying spatial working memory is not understood fully. Here, we confirm previous studies showing impairment on the Delayed Recognition Span Task in the spatial condition (DRSTsp), and increased in vitro action potential firing rates (hyperexcitability), across the adult life span of the rhesus monkey. We use a bump attractor model to predict how empirically observed changes in the aging dlPFC affect performance on the Delayed Response Task (DRT), and introduce a model of memory retention in the DRSTsp. Persistent activity—and, in turn, cognitive performance—in both models was affected much more by hyperexcitability of pyramidal neurons than by a loss of synapses. Our DRT simulations predict that additional changes to the network, such as increased firing of inhibitory interneurons, are needed to account for lower firing rates during the DRT with aging reported in vivo. Synaptic facilitation was an essential feature of the DRSTsp model, but it did not compensate fully for the effects of the other age-related changes on DRT performance. Modeling pyramidal neuron hyperexcitability and synapse loss simultaneously led to a partial recovery of function in both tasks, with the simulated level of DRSTsp impairment similar to that observed in aging monkeys. This modeling work integrates empirical data across multiple scales, from synapse counts to cognitive testing, to further our understanding of aging in non-human primates.

Introduction

Spatial working memory in rhesus monkeys, as in humans, is mediated by the action potential firing activity of neurons in the dorsolateral prefrontal cortex (dlPFC) (reviews: Funahashi, 2017; Constantinidis and Qi, 2018; Miller et al., 2018). And, as in humans (Albert, 1993; Salthouse et al., 2003; Fisk and Sharp, 2004; Rhodes, 2004; Sorel and Pennequin, 2008), the effects of aging on working memory is heterogeneous—while a significant proportion of rhesus monkeys become increasingly impaired on tasks of spatial working memory during normal aging (“unsuccessful agers”), a significant number of “successful agers” show no signs of impairment (Lacreuse et al., 2005; Moore et al., 2006, 2017). Furthermore, for impaired subjects, impairment typically begins fairly early in the aging process, during early middle age (Moore et al., 2005, 2006, 2017). This distribution of impaired and spared monkeys across the lifespan enables assessment of brain changes associated not simply with aging, but with cognitive performance per se in this model of normal human aging. This approach has been taken in studies of normal aging by our group and others (reviews: Peters, 2007; Luebke et al., 2010, 2015; Wang et al., 2011; Peters and Kemper, 2012), in which monkeys were assessed on working memory tasks such as the Delayed Response Task (DRT) and the Delayed Recognition Span Task in the spatial condition (DRSTsp), and their brains subsequently examined for a wide variety of parameters. Thus, declines in spatial working memory have been associated with many sub-lethal changes to the structure and function of neurons, glial cells, and white matter pathways in the dlPFC of the rhesus monkey (Peters, 2007, 2009; Peters et al., 2008; Bowley et al., 2010; Luebke et al., 2010; Shobin et al., 2017). Reductions in synapses and increased dystrophy of white matter pathways begin in early middle age; for example, Peters et al. (2008) showed a continuous decrease in the number of excitatory and inhibitory synapses, detectable even in middle age in the monkey dlPFC.

In addition to well-documented structural changes, functional alterations to the electrophysiological properties of supragranular neurons in the aging monkey dlPFC have been reported in several studies (Chang et al., 2005; Wang et al., 2011; Coskren et al., 2015). Layer 3 dlPFC pyramidal neurons exhibit significantly increased evoked action potential (AP) firing rates in vitro for aged compared to young monkeys (Chang et al., 2005; Coskren et al., 2015). In contrast, another study (Wang et al., 2011) reported that in vivo firing rates of dlPFC DELAY neurons decreased with aging during the DRT, for monkeys performing the task well. Since persistent firing patterns of dlPFC pyramidal neurons encode precisely tuned spatial and temporal information during a working memory task, exploring how these findings might complement each other may reveal new predictions about spatial working memory decline.

While the growing body of literature over the past 20 years has documented numerous changes to neurons, glial cells and white matter pathways in the aging brain that are associated with cognitive changes, which changes are the key determinants of age-related cognitive declines has not yet been firmly established (Hof and Morrison, 2004; Luebke et al., 2010; Morrison and Baxter, 2012; Peters and Kemper, 2012; Konar et al., 2016; Motley et al., 2018; Cleeland et al., 2019). This lack of insight is partially due to the difficulty in first, selectively targeting individual variables (such as firing rate or synapse number) in humans or experimental subjects without, second, also altering upstream or downstream effectors. A powerful way to gain insight into which age-related changes are most consequential for dlPFC network behavior during normal aging is through computational models that are constrained by empirical data. This is the approach we use here through a systematic study of the “bump attractor” neural network model, and its extension to memory retention in a new task similar to the DRSTsp.

The bump attractor network has been used to model spatial working memory in the DRT (Compte et al., 2000; Wang et al., 2011, 2013; Wimmer et al., 2014; Wu et al., 2016 for review). The bump attractor exhibits persistent firing activity after an initial stimulus disappears, encoding a memory of the stimulus location. The bump attractor was used to show that an age-related loss of synaptic strength could account for reduced firing rates of dlPFC DELAY neurons (Wang et al., 2011), but the effects of increased excitability of pyramidal neurons seen in vitro on the function of bump attractor models has not been examined. We have previously used computational modeling to predict which intrinsic electrophysiological and morphological properties of individual pyramidal neurons contribute to the action potential firing rate increases seen in aging (Coskren et al., 2015; Rumbell et al., 2016). Here, we extend these previous single neuron studies using network-level modeling to predict the effect of empirically observed physiological and structural changes in the aging rhesus dlPFC on cognitive behavior. We obtained whole-cell patch clamp recordings of dlPFC pyramidal neurons from behaviorally characterized monkeys across the adult lifespan. Empirical results indicated that previously reported physiological changes seen in aging are already present in middle age, and are correlated with cognitive impairment. This modeling work aims to connect age-related changes in non-human primates across multiple empirical scales, from synapse counts and physiology of single neurons up to network output and cognitive performance. We predict that the increased firing rates and reduced synapse density observed with aging may partially compensate one another, but still are sufficient to induce substantial cognitive deficits in the DRT and DRSTsp.

Materials and Methods

Experimental Subjects

The rhesus monkeys (Macaca mulatta) studied here were a part of a larger study of normal aging of the brain. Data were obtained from a total of 9 young (8.9 ± 0.5 years old), 18 middle-aged (18.2 ± 0.5 years old), and 12 aged (24.7 ± 0.7 years old) monkeys (Table 1). Prior to the electrophysiological experiments, all but 3 young monkeys completed cognitive testing of the DRSTsp (Moore et al., 2017). Monkeys were obtained from the Yerkes National Primate Research Center at Emory University (Atlanta, GA, USA), and housed in the Laboratory Animal Science Center (LASC) at the Boston University School of Medicine (BUSM) under strict accordance with the guidelines established by the NIH Guide for the Care and Use of Laboratory Animals and the US Public Health Service Policy on Humane Car and Use of Laboratory Animals. Both Yerkes National Primate Research Center and BUSM LASC are fully accredited by the Association for Assessment and Accreditation of Laboratory Animal Care.

TABLE 1

Table 1. Experimental subjects.

Behavioral Assessment

A total of 36 subjects were tested on the DRSTsp, 6 young, 18 middle-aged, and 12 aged monkeys. The DRSTsp is a short-term memory task in which the subject distinguishes a novel cue (spatial location) from an increasing set of recently presented, familiar cues (spatial locations) (Moore et al., 2017 and Figure 1A). All stimuli were identical brown disks, each covering one of 18 small wells arranged on a testing board in a 3 × 6 matrix. On the first trial one of the wells was baited with a reward then covered with a disk. Once a screen was raised, the monkey moved the disk to receive the reward. The screen was lowered before the second trial, when that same well (with no reward) was covered with a disk and a second well was baited with a reward and covered with an identical disk (10 s delay). Subsequent trials include the previous cue plus one new cue, until the monkey is unable to identify novel cues (spatial locations). Monkeys were tested for 10 consecutive days. The cognitive score recorded is the recognition span: the mean number of unique cues the monkey identified correctly before making a mistake. Some of the subjects were assessed on the DRSTsp more than once over a period of several years as part of a longitudinal study of aging. The DRSTsp score used in the present study for longitudinally assessed subjects was the one most proximal to the date the subject was sacrificed.

FIGURE 1

Figure 1. Cognitive performance and physiological variables vary with aging. (A) DRSTsp setup. Brown circles represent disks covering the wells arranged on a board in a 3 × 6 matrix. (B) DRSTsp performance vs. age. (C) Input resistance vs. age; regression includes subject as random effect. Data shown as mean ± SEM for each subject. (D) FR vs. input resistance for + 130 pA current step; these variables were highly correlated at all levels of current injection (not shown). (E) Left: Trains of APs evoked by 2-s depolarizing current steps for representative young, middle-aged, and aged neurons. Scale bar: 50 mV and 500 ms. Current steps were applied from resting membrane potentials, shown as arrows. Right: FR vs. injected current. (F) DRSTsp performance vs. AP firing rate in response to a +130 pA current injection. In all graphs, young, middle-aged, and aged subjects are shown in black, blue, and red, respectively. Data were obtained from a total of 9 young, 18 middle-aged, and 12 aged monkeys (Table 1). All but 3 young monkeys completed cognitive testing of the Delayed Recognition Span Task in the spatial condition. Whole-cell patch clamp recordings from 324 pyramidal neurons from 38 subjects were used.

Preparation of Cortical Slices and Whole-Cell Patch Clamp Experiments

Monkeys were sacrificed as described in our previous publications (Amatrudo et al., 2012; Medalla and Luebke, 2015; Medalla et al., 2017). Briefly, following perfusion with Krebs-Henseleit solution, a block of tissue (~1 cm³) was removed from the left lateral prefrontal cortex, and sectioned in ice-cold oxygenated Ringer's solution (concentrations, in mM: 26 NaHCO₃, 124 NaCl, 2 KCl, 3 KH₂PO₄, 10 glucose, 1.3 MgCl₂, pH 7.4; Sigma-Aldrich) into 300 μm-thick coronal slices with a vibrating microtome. Slices were immediately transferred into room temperature, oxygenated Ringer's solution, for a minimum of 1 h equilibration prior to recordings.

Individual slices were placed into submersion-type recording chambers (Warner Instruments) positioned on Nikon E600 infrared-differential interference contrast (IR-DIC) microscopes. During recordings, slices were continuously superfused with oxygenated, room-temperature Ringer's solution at a flow rate of 2–2.5 ml/min.

Standard tight-seal, whole-cell patch clamp recordings were performed on layer 3 pyramidal neurons as described previously (Amatrudo et al., 2012; Luebke et al., 2015; Medalla and Luebke, 2015; Medalla et al., 2017). Electrodes were fabricated using a Flaming/Brown micropipette puller (MODEL P87, Sutter Instruments), and filled with potassium methane sulfonate (KMS)-based internal solution (concentration, mM: 122 KCH₃SO₃, 2 MgCl₂, 5 EGTA, 10 NaHEPES, 1% biocytin, pH 7.4; Sigma-Aldrich). Experiments were performed using EPC-9 or EPC-10 amplifiers (HEKA Elektronik) controlled with PatchMaster acquisition software (HEKA Elektronik). Signals were low pass-filtered at 10 kHz, and access resistance was monitored and maintained throughout each experiment. Current clamp protocols were implemented as previously described (Amatrudo et al., 2012; Luebke et al., 2015; Medalla and Luebke, 2015; Medalla et al., 2017), to assess the input resistance (in MΩ) and repetitive firing rate (FR, measured in Hz) of the neurons.

Whole-cell patch clamp recordings from 324 layer 3 dlPFC pyramidal neurons from 38 subjects were used in this study. In total, 40 neurons for the young, 188 for the middle-aged and 96 for the aged subjects.

The Bump Attractor Network Model

The “bump attractor” network model was originally designed to simulate a spatial working memory task for awake behaving monkeys known as the DRT or the Oculomotor Delay Task (Figure 2A, Goldman-Rakic, 1995; Compte et al., 2000; Wang et al., 2011; Wimmer et al., 2014). After the monkey fixates on the center of a computer screen (for 1 s in our simulations), a light appears during the cue period in one of eight directions around a circle (occurring at 0° and lasting 0.5 s in our simulations). After the cue period the light turns off, and the monkey fixates on the center of the screen during a delay period (2 s in our simulations). During the response period (0.7 s in our simulations), the monkey shifts its gaze (makes a saccade) to indicate its memory of the location of the stimulus. The selective firing of pyramidal neurons in the dlPFC to different spatial locations (their preferred directions) is thought to be the mechanism by which stimulus location is encoded. Thus, monkeys performing this task successfully have neurons that exhibit persistent neural activity during the delay period tuned to the stimulus location (Goldman-Rakic, 1995; Compte et al., 2000; Wimmer et al., 2014). This is achieved in the bump attractor model with excitatory neurons tuned to each spatial location (so that we can label each neuron by that location—its “preferred direction”), with the synaptic strength between pairs of excitatory neurons determined by the difference between the neurons' preferred directions.

FIGURE 2

Figure 2. The bump attractor model of the DRT. (A) The DRT includes fixation, cue, delay, and response periods of 1, 0.5, 2, and 0.7 s, respectively. The ×, black dots, and green dot, respectively, represent the fixation point, the 8 possible stimulus locations, and actual stimulus location (0°) in the simulations. Arrow indicates the saccade direction during the response period. (B) The recurrent excitatory-inhibitory network. (C) The f-I curve for young model excitatory neurons. (D) Ring connectivity among excitatory neurons. M_ji is the connectivity matrix element between excitatory neurons j and i, depending on the difference between their preferred directions (θ_j and θ_i, respectively). Each neuron is strongly connected to its nearest neighbors (red line), with the connection strength decaying with distance (yellow and green lines). The neuron with preferred direction 0° receives the strongest stimulus current, which follows a Gaussian distribution centered at 0°, generating a bump of activity during the delay period (persistent activity tuned to the stimulus location) when excitation and inhibition in the network are well-balanced. (E) Sample DRT model output (side and top views) when the 0° stimulus location is encoded correctly: excitatory neuron FR vs. simulation time with neurons labeled by their preferred direction.

The bump attractor network used here assumed a firing-rate model for individual neuronal dynamics, modified from Wimmer et al. (2014). The code was written in Matlab and is available on ModelDB (McDougal et al., 2017) at Accession #256610. The local cortical network was composed of N_E = 640 excitatory (pyramidal) neurons (80%) and N_I = 160 inhibitory neurons (interneurons, 20%) (Abeles, 1991; Braitenberg and Schütz, 1991) (Figure 2B). Each excitatory and inhibitory neuron received three general types of synaptic inputs: from excitatory and inhibitory connections with other neurons in the same network, and excitatory connections from other unspecified areas external to the network. The modeled excitatory neurons also received an external input current representative of the “cue.” The synaptic excitation of each neuron was modeled with two distinct equations, separately representing synaptic inputs mediated by AMPA and NMDA receptors. The voltage dependence of the NMDA currents was omitted from the present study, since Compte et al. (2000) showed that the slow kinetics of NMDA-mediated synaptic transmission was the most important feature for generating stable persistent activity.

The equations for FR of each excitatory and inhibitory neuron (r_E and r_I) were given by

$\begin{array}{l}{\tau }_E\frac{d{r}_E^j}{dt}=-r_E^j+f_E\left(I_E^j\right)+{\sigma }_E{\xi }_E^j\left(t\right)\\ {\tau }_I\frac{d{r}_I^k}{dt}=-r_I^k+f_I\left(I_I^k\right)+{\sigma }_I{\xi }_I^k\left(t\right),& (1)\end{array}$

for j = 1, …, N_E and k = 1, …, N_I. The parameters τ_E and τ_I represented the membrane time constants for excitatory and inhibitory neurons, respectively, (in ms), with f(I) the FR activation function (or f-I curve, in Hz), total synaptic input currents I_E and I_I (in pA), and Gaussian white noise inputs ${\sigma }_E{\xi }_E^j(t)$ and ${\sigma }_I{\xi }_I^k(t)$ (in Hz), where σ_E and σ_I were the standard deviations. In this study, τ_E = 20 ms, τ_I = 10 ms, σ_E = 1, and σ_I = 3. The activation function for the firing rate of each neuron took the piecewise form (Brunel, 2003; Wimmer et al., 2014; Roxin and Compte, 2016):

$\begin{array}{ll}f\left(I\right)=\{\begin{array}{cc}0& I\le 0\\ {\nu }_c{\left(\frac{I}{I_c}\right)}^2& 0<I<I_c\\ 2{\nu }_c\sqrt{\frac{I}{I_c}-\frac{3}{4}}& I_c\le I\end{array}& (2)\end{array}$

The parameter I_c represented a current near rheobase (measured in pA), with corresponding firing rate ν_c (in Hz) (see Figure 2C). Thus, the parameter pairs {ν_ce, I_ce} and {ν_ci, I_ci} uniquely defined the f-I curves for excitatory and inhibitory neurons. Here, I_ci = 20 pA and ν_ci = 50 Hz, values approximately in the regime for Chandelier cells in the monkey PFC (Zaitsev et al., 2009; Povysheva et al., 2013). The parameters ν_ce and I_ce were chosen to fit our empirical data as described below. See also the Supplementary Information for a general, less excitable form of both f-I curves.

The total synaptic input currents (in pA) for each neuron were

$\begin{array}{l}{I}_E^j\left(t\right)=I_{Ea}^j\left(t\right)+I_{En}^j\left(t\right)+I_{Eg}^j\left(t\right)+I_{0E}+I_{stim}^j\\ I_I^k\left(t\right)=I_{Ia}^k\left(t\right)+I_{In}^k\left(t\right)+I_{Ig}^k\left(t\right)+I_{0I},\end{array}$

with time-dependent synaptic currents given by

$\begin{array}{r}{\tau }_a\frac{d{I}_{Ea}^j}{dt}=-I_{Ea}^j+{\displaystyle \sum _{i=1}^{N_E}}{G}_{EEa}{M}_{ji}{r}_E^i\\ {\tau }_n\frac{d{I}_{En}^j}{dt}=-I_{En}^j+{\displaystyle \sum _{i=1}^{N_E}}{G}_{EEn}{M}_{ji}{r}_E^i\\ {\tau }_g\frac{d{I}_{Eg}^j}{dt}=-I_{Eg}^j-G_{IE}{\overline{r}}_I\\ {\tau }_a\frac{d{I}_{Ia}^k}{dt}=-I_{Ia}^k+G_{EIa}{\overline{r}}_E\\ {\tau }_n\frac{d{I}_{In}^k}{dt}=-I_{In}^k+G_{EIn}{\overline{r}}_E\\ {\tau }_g\frac{d{I}_{Ig}^k}{dt}=-I_{Ig}^k-G_{II}{\overline{r}}_I.& (3)\end{array}$

The subindices a and n represented the AMPA and NMDA synaptic receptors for the glutamatergic synapses, and g the GABA_A receptor. Parameters τ_a, τ_n, and τ_g were the synaptic decay time constants for AMPA, NMDA, and GABA_A (Hestrin et al., 1990; Spruston et al., 1995; Salin and Prince, 1996; Xiang et al., 1998). We set τ_a = 2 ms, τ_n = 100 ms, τ_g = 10 ms. The quantity $\overline{r}$ represented the mean FR value of each neuron type, and the G_pq parameters with {p, q} = {E, I} were the constant synaptic weights. The connectivity matrix M was a circular Gaussian function describing the translation-invariant connections among the excitatory neurons. Labeling each neuron by its preferred direction, the matrix element M_ji for two neurons j and i depended on the difference between their preferred directions. That is, M_ji∝ |θ_j − θ_i|, with

$\begin{array}{l}{\theta }_j=\left(\frac{2j}{N_E}-1\right)\pi \end{array}$

being the preferred stimulus direction of the excitatory neuron j, where a complete circumference was divided into N_E angles (see Figure 2D). The quantities I_0E and I_0I represented constant excitatory input currents from other brain areas to excitatory and inhibitory neurons, respectively, and $I_{stim}^j$ the input current to each excitatory neuron due to stimuli during cue periods of the tasks described below. Here, I_0E = 80 pA and I_0I = 15 pA. The parameters G_pq were varied using the Latin Hypercube Sampling design, as described below.

The input stimulus currents for each excitatory neuron representing a stimulus at 0° was modeled as in Wimmer et al. (2014) as

$\begin{array}{l}{I}_{stim}^j=\frac{I_{st}{e}^{c·cos{\theta }_j}}{{\displaystyle {\sum }_{j=1}^{N_E}}{e}^{c·cos{\theta }_j}},\end{array}$

which provided the bump-shaped input current to all excitatory neurons (Figure 2D). Parameter I_st defined the strength of the stimulus and c the concentration of excitatory neuron to excitatory neuron connectivity. In our simulations, I_st = 40,000 pA, and c = 1.5 for the DRT model and c = 20 for the DRSTsp model. A typical simulation of the DRT is shown in Figure 2E.

Modeling Effects of Aging on Individual Neurons

Increased AP Firing Rate (Hyperexcitability) of Pyramidal Neurons With Aging

We fit parameters ν_ce and I_ce of the excitatory neurons f-I curve to in vitro FR data from the subjects described above. We computed the mean FR vs. different input currents of all neurons from each monkey, then computed a grand mean of FR vs. input current for monkeys in each age group (Figure 3B). Fits of the parameters ν_ce and I_ce to these data were estimated manually using the FindFit function from Mathematica (Wolfram Research, Champaign, IL). The value I_ce = 98 pA fit all three age groups well, and also agreed with rheobase values for regular spiking pyramidal neurons with low input resistance (Zaitsev et al., 2012). We used ν_ce = 5, 7, and 9 Hz to fit the f-I curve to data from the young, middle-aged, and aged groups, respectively (Figure 3B). Thus, to fit in vitro FR data it was sufficient to increase the single parameter ν_ce. Assuming the magnitude of changes to f-I curves with aging would be similar at physiological vs. room temperature, we used these parameter values directly in our “young,” “middle-aged,” and “aged” network models. We also repeated some DRT simulations after fitting a generalized, less excitable firing rate activation function to the young, middle-aged, and aged monkey data (Supplementary Information and Supplementary Figure 1a).

FIGURE 3

Figure 3. DRT model output varied across the parameter space. (A) From left to right: Example model output showing (1) under-excitation; (2) TPA-S; (3) partial over-excitation; and (4) over-excitation. (B) The model f-I curve fit (solid lines) to empirical AP firing rates of pyramidal neurons of young, middle-aged, and aged subjects, averaged for each age group (black, blue, and red, respectively, shown as mean ± SEM at each injection level). (C) DRT model output for 4,200 points of the parameter space LHS, shown in 3D projections across the subspace of the excitatory (G_EEa and G_EEn) and the inhibitory (G_IE) synaptic weights of pyramidal neurons as their excitability increased (ν_ce = 3, 5, 7, and 9 Hz). The numbers in the 3D projection for ν_ce = 5 Hz indicate each type of network performance shown in (A), represented as blue, green, yellow, and red dots, respectively, in each graph. (D) Number of points in each LHS maintaining TPA-S as ν_ce increased. Dark gray bar indicates the young model cohort. (E) Projections across the (G_EEa, G_EEn) subspace showed that points maintaining TPA-S shifted to lower (G_EEa, G_EEn) values as ν_ce increased. (F) Mean synaptic weights (G_EEa, G_EEn, G_IE, G_EIa, G_EIn, and G_II) for all points maintaining TPA-S as ν_ce increased. G_EEa, G_EEn, and G_IE values shown as purple, pink, and red open circles, respectively; G_EIa, G_EIn, and G_II values shown as dark green, light green, and open blue triangles, respectively. S.E.M. bars lie beneath the symbols.

Loss of Excitatory and Inhibitory Synapses With Aging

Peters et al. (2008) reported a loss of the number of asymmetric and symmetric synapses (excitatory and inhibitory, respectively) in Layers 2/3 of area 46 of rhesus monkey as adult monkeys age: 10% in middle-aged monkeys, and 30% in aged. For the “young monkey model cohort” DRT studies below we represented these respective losses as a 10% decrease in the synaptic weights G_EEa, G_EEn, and G_IE in middle-aged models, and 30% decrease in these parameters in aged models. For the “young monkey model cohort” DRSTsp studies, we decreased these parameters in a semi-continuous way for the simulated “aged” networks.

Exploring the DRT Model Parameter Space

We conducted two general types of parameter exploration studies for the simulated DRT. The first were sweeps across parameter space using a space-filling Latin Hypercube Sampling (LHS) design, as in (Johnson et al., 1990; Rumbell et al., 2016). The LHS design identified 4,200 points (networks) across the parameter space of synaptic weights (G_EEa, G_EEn, G_IE, G_EIa, G_EIn, and G_II), representing semi-random combinations of the synaptic weights homogeneously distributed across the 6-dimensional space (Rumbell et al., 2016) (see Table 2). Bounds of the parameter space were chosen so that the mean firing rate of excitatory neurons in all “young” model networks performing the DRT successfully had mean firing rates of excitatory neurons near 35 Hz, a realistic value for the PFC (see Figure 5, black curves). We examined the results of network simulations throughout the LHS as the excitatory firing rate parameter ν_ce varied from 2 to 10 Hz (Figure 3). All other parameters remained fixed for all simulations. For each of the 4,200 parameter combinations in the LHS, the DRT was simulated 7 times to account for variations due to noise in the firing rate equations (Equation 1), and results were averaged across all repetitions. We identified five main categories of network behavior: maintaining tuned persistent activity (TPA) to any angle orientation until the end of the delay period; maintaining TPA tuned to the original stimulus location at 0° (TPA-S) until the end of the delay period; under-excited networks (either unable to maintain TPA until the end of the delay, or unable to generate a bump of activity at all); over-excited networks (all neurons firing at a similar rate by the beginning of the cue period and for the remainder of the task); and partially over-excited networks (initially maintaining TPA, but becoming over-excited sometime during or after the cue period).

TABLE 2

Table 2. Parameter ranges for each synaptic weight (in pA·s) used to generate the Latin Hypercube Sample in the DRT model.

For the second type of parameter exploration, we first defined a young model cohort as points in the LHS associated with ν_ce = 5 Hz which maintained TPA-S until the end of the delay period in at least one of the 7 simulation repetitions (shown as the ν_ce = 5 Hz bar in Figure 3D). We then perturbed parameters of the young model cohort to simulate the effects of morphological and physiological changes that have been observed with aging, or that we propose as novel possibilities (Figure 4). Eight kinds of perturbations were made to the young model cohort: (1) The increased excitability condition: increasing ν_ce from 5 Hz to 7 and 9 Hz, keeping all else constant, to simulate the effect of increasing excitability of individual neurons as monkeys age beyond young adult. (2) The synapse loss condition: reducing the values of G_EEa, G_EEn, and G_IE by either 10 or 30%, simulating the effect of fewer excitatory and inhibitory synapses onto pyramidal neurons for middle-aged and aged monkeys, respectively. (3) Combined increased excitability and synapse loss condition: simultaneously perturbing the excitability and synapse parameters from (1) and (2). Below this perturbation is also called the “observed data conditions.” (4) The observed data conditions (combined increased excitability and synapse loss conditions) as in (3), plus increasing the excitability of the inhibitory neurons in the same proportion as in the excitatory neurons: increasing ν_ci from 50 Hz to 70 and 90 Hz for middle-aged and aged simulations, respectively. Each of these four perturbations was then repeated after adding short-term synaptic facilitation (as described below) to the networks of the young model cohort. As with the LHS, each young cohort network simulation was repeated 7 times with different randomized seeds and results were averaged.

FIGURE 4

Figure 4. Perturbing the DRT young model cohort. Young cohort defined as networks from the LHS maintaining TPA-S when ν_ce = 5 Hz (679 networks, dark gray bar in Figure 3D, and bar for young simulations in all panels in this figure). For middle-aged and aged simulations: (A) Applied the synapse condition to young cohort, holding the excitability parameters constant. (B) Applied the excitability condition, holding synapse parameters constant. (C) Applied both observed data conditions simultaneously (excitability and synapse conditions). (D–F) Analogous perturbations as in (A–C), after adding synaptic facilitation to each network.

During in vivo recordings, Wang et al. (2011) found a similar reduction of the firing activity of the DELAY neurons during the cue and delay periods of a DRT with aging, in monkeys performing the task well, and they used computational modeling of the DRT to show that an age-related loss of synaptic strength could account for that reduction. Here we studied if this was still true when the excitability of the excitatory neurons increased at the same time (Figure 5A). To do this, we examined the mean firing rate of excitatory neurons for the networks in the LHS maintaining TPA-S as ν_ce increased from young to aged models (for ν_ce = 5, 7, and 9 Hz). In these simulations we reduced the stimulus scaling factor, I_st, by 10 and 30% for middle aged and aged networks, respectively. This made firing rates during the cue period more consistent with data from Wang et al. (2011), without significantly affecting firing rates during the delay period (driven primarily by G_EEa and G_EEn).

FIGURE 5

Figure 5. Firing rates and total synaptic current of excitatory neurons in DRT model networks maintaining TPA-S. (A) Left, mean FR of excitatory neurons vs. simulated time for all sampled networks maintaining TPA-S when ν_ce = 5, 7, and 9 Hz (networks shown in Figure 3E). Shown are neurons with preferred directions matching the 0° stimulus; gray curve shows the mean FR for neurons with a 180° preferred direction, which fired in none of these networks. Right, mean total synaptic input current to excitatory neurons with the 0° preferred direction vs. simulation time for the same networks. (B) Mean FR of excitatory neurons for all networks maintaining TPA-S for the model cohort perturbations when the observed data conditions were applied (Figures 4C,F). Left: without facilitation. Right: with facilitation. (C) Top row: DRT model cohort with applied both observed data conditions and increased excitability of inhibitory neurons. Bottom row: Mean FR of excitatory neurons for all networks maintaining TPA-S. Left: without facilitation. Right: with facilitation. In all panels, young, middle-aged, and aged groups are shown in black, blue, and red, respectively; shaded regions indicate 95% confidence intervals about each mean.

Simulating Memory Retention in a DRSTsp—Like Task

Most of the computational models developed to understand the neural mechanisms underlying working memory are based on delayed memory tasks with just one stimulus to remember during the delay period. A small but growing number of computational studies simulate behavioral tasks that involve the maintenance of several stimuli in memory. Developing models able to maintain the memory of multiple stimuli is challenging, due to interference produced by the simultaneous activation of the many neurons encoding the different stimuli. One way to avoid this problem is to have each stimulus encoded by sparse bumps of activity that does not overlap with firing activity encoding other stimuli (Amit et al., 2003). The levels of inhibition (Edin et al., 2009) and excitation (Wei et al., 2012) in the network are two factors that control the working memory capacity for maintaining multiple stimuli. The balance between them influences the number of stimuli that can be remembered and how bumps fail (fading out vs. merging). Continuous attractor models require fine-tuning of the parameters to control the balance of excitation and inhibition. Such fine tuning can be relaxed by introducing short-term synaptic plasticity. Mongillo et al. (2008) proposed that calcium-mediated synaptic facilitation could be the neural mechanism underlying working memory. Rolls et al. (2013) showed that synaptic facilitation can make the system more robust to model parameters. Both found that synaptic facilitation increases the multi-item working memory capacity in a discrete attractor network. Mi et al. (2017) also found that working memory capacity increased with the time constant of synaptic depression. Some bump attractor models of the DRT have included synaptic facilitation (Itskov et al., 2011; Hansel and Mato, 2013), but in the past it has often been omitted.

We extended the simulation setup of the DRT to create a model of memory retention in a simplified task that is similar to the DRSTsp. Our original bump attractor model had great difficulty maintaining the memory of multiple stimuli simultaneously, due to interference between the stimulus responses (Amit et al., 2003; Barak and Tsodyks, 2014) and because during the first cue and delay periods, all neurons not involved in the bump activity were highly inhibited. Thus, when introducing subsequent stimuli, the stimulus input current (equal for all stimuli) was insufficient to induce subsequent bumps of firing activity. This difficulty was addressed by adding short-term synaptic facilitation in the excitatory-to-excitatory synaptic connections. As in Itskov et al. (2011), we assumed that the overall behavior of the synapses was facilitating, and modified the first two equations in Equation (3) as

$\begin{array}{l}{\tau }_a\frac{d{I}_{Ea}^j}{dt}=-I_{Ea}^j+{\displaystyle \sum _{i=1}^{N_E}}{G}_{EEa}{M}_{ji}{u}^i{r}_E^i\\ {\tau }_n\frac{d{I}_{En}^j}{dt}=-I_{En}^j+{\displaystyle \sum _{i=1}^{N_E}}{G}_{EEn}{M}_{ji}{u}^i{r}_E^i,\end{array}$

with the facilitation variable u given by

$\begin{array}{l}{\tau }_f\frac{du}{dt}=-\left(u-U\right)+{\tau }_{f}U{r}_E\left(1-u\right).\end{array}$

Here u is a utilization variable defining the fraction of total neuronal resources that each spike is able to access, representing increased calcium levels at presynaptic terminals (Mongillo et al., 2008). Without firing, u decays to its baseline value U with time constant τ_f. In all our simulations U = 0.001 and τ_f = 1, 500 ms.

Behavioral data showed a maximum DRSTsp score of 3.61 for one of the middle-aged subjects (Figure 1B). The cognitive score recorded is the recognition span: the mean number of stimuli the monkey identified correctly before making a mistake. Thus, we designed our computational task so that perfect performance gave a simulated DRSTsp score (hereafter called DRST_sim) of four, corresponding to successful encoding of four stimuli equally spaced around the ring from the DRT (Figure 6A). Our DRSTsp model simulated a “pre-cue” period of 1 s, followed by alternating cue and delay periods (lasting 0.5 and 2 s, respectively). The pre-cue and delay periods represented times in the real DRSTsp when the screen was down while appropriate wells were baited. In each subsequent cue period, a new stimulus was presented in addition to any stimuli that were presented before. In the first cue period a stimulus was presented at 0°; in the second, stimuli were presented at 0° and 90°. In the third cue period stimuli were presented at 0°, 90°, and −90°, and in the fourth stimuli were presented at 0°, 90°, −90°, and 180°. Stimuli were identical each time at each presentation, with I_st = 40,000 pA and c = 20. After the first cue period, if a network simulation encoded all stimuli during cue period n and maintained all n−1 stimuli throughout the previous delay period (by maintaining a bump for each stimulus with maximum FR of at least 5 Hz), we assumed that the novel stimulus was identified correctly. In this way, our task models the ability of our network to retain the presented stimuli, but does not model the process of choosing which stimulus among all presented is the novel one. Our task also assumes that each trial period lasts the same amount of time, whereas monkeys performing the DRSTsp terminate each successful trial by moving the disk corresponding to their choice before the maximum trial time expires. Sample networks with DRST_sim = 4, 3, 2, and 1 are shown in Figure 6B.

FIGURE 6

Figure 6. Memory retention model of the DRSTsp shows impairment under neuronal hyperexcitability. The ability to maintain the memory of a successively increasing number of stimuli strongly depended on the excitability of individual pyramidal neurons, but changes in the synaptic weighs partially compensated for the increased neuronal excitability. (A) Simulated task included a pre-cue period, with four cue periods separated by three delay periods. Green circles indicate the stimuli presented in each cue period. (B) Examples of DRSTsp model output, with FR shown in color (log scale) for neurons labeled by their preferred direction vs. simulated time. From left to right: DRST_sim = 4, optimal performance, all four stimuli were encoded correctly and recalled; DRST_sim = 3, where the first three stimuli were encoded correctly, but persistent activity representing the third was lost during the subsequent delay period; and DRST_sim = 2 and DRST_sim = 1 where only the first two or one stimuli were encoded correctly and recalled. Cue periods shown between white lines; arrowheads indicate when persistent activity was lost. (C) DRSTsp model output for 4,200 points of the parameter space LHS, shown in 3D projections across the (G_EEa, G_EEn, G_IE) subspace as ν_ce increased. Points with DRST_sim = 4, 3, 2, and 1 are shown as green, red, yellow, and blue dots, respectively. (D) Histogram of DRST_sim values in the LHS as ν_ce increased (ν_ce = 5, 7, and 9 Hz shown as black, blue, and red bars, respectively). (E) The first three simulations from (B), repeated for ν_ce = 7 (top row) and 9 Hz (bottom row). (F) Mean synaptic weights for all LHS points with DRST_sim = 4 for ν_ce = 5, 7, and 9 Hz. Top: G_EEa, G_EEn, and G_IE shown as purple, pink, and red open circles. Bottom: G_EIa, G_EIn, and G_II shown as dark green, light green, and open blue triangles. S.E.M. bars lie beneath the symbols.

With this new DRST_sim established, we again used LHS to explore how task performance varied across the parameter space of the six synaptic weights (G_EEa, G_EEn, G_IE, G_EIa, G_EIn, and G_II; see Table 3) and for ν_ce = 5, 7, and 9 Hz. As for the DRT, each LHS had 4,200 points (Figure 6C); we chose the bounds of the parameter space so that all 4,200 networks encoded at least the first stimulus for all three ν_ce values. Simulations were performed only once, since synaptic facilitation made the network performance much less sensitive to random fluctuations.

TABLE 3

Table 3. Parameter ranges for each synaptic weight (in pA·s) used to generate the Latin Hypercube Sample in the DRSTsp model.

We also defined a young monkey model cohort for the DRSTsp model to be a group of 40 randomly chosen points in the LHS with ν_ce = 5 Hz (Figure 6D, young case: black bars), half with DRST_sim = 4 and the other half with DRST_sim = 3. We then “aged” the networks in this cohort by applying increased excitability and synaptic loss conditions, and increasing the excitability of the inhibitory neurons, as follows. We set bounds for ν_ce just beyond our young and aged values (4.5 and 9.5 Hz) and generated 40 randomly distributed values between them. We then defined two piecewise linear functions that interpolated the synapse loss conditions from the DRT (0, 10, and 30% loss when ν_ce = 5, 7, 9 Hz, respectively, slope 5% loss/Hz when ν_ce < 7 Hz and 10% loss/Hz when ν_ce ≥ 7), and computed a corresponding level of synapse loss for each of the 40 ν_ce values. We computed corresponding values for ν_ci by multiplying each ν_ce value by 10, consistent with the relationship between these two parameters elsewhere in our study. We computed an associated age for each ν_ce value similarly, interpolating the mean age for the three empirical groups (8.9, 18.2, and 24.7 years for young, middle aged, and aged groups) when ν_ce = 5, 7, 9 Hz, respectively. We induced the increased excitability and synapse loss conditions of these 40 values separately and together, plus increased excitability of inhibitory neurons, and then computed the corresponding DRST_sim for each transformed member of the young cohort (6,400 simulations in all: 160 for each point in the cohort). We then calculated the mean DRST_sim for each simulated age group (young group, 7.5 < years < 13; middle-aged group, 13 ≤ years < 21; and aged group, 21 < years < 26). To summarize the change in performance across the simulated age span (Figure 7), we calculated the relative change in the mean DRST_sim from the mean ages of the young simulations to the aged simulations (9.9 and 24.2 years, respectively) under each perturbation condition.

FIGURE 7

Figure 7. Perturbing the DRSTsp young model cohort. Each dot indicates the mean DRST_sim of the 40 young cohort points. (A) Synapse condition only; slight improvement (+4%) from the young to the aged simulations. (B) Excitability condition only; strong impairment (−45%) from young to aged. (C) Applied both observed conditions simultaneously, resulting in just a −14% impairment. (D) Applied both observed data conditions and increased excitability of inhibitory neurons; −17% impairment.

Statistics

All analyses were performed in RStudio 1.1.463 (RStudio Team, 2015; Mangiafico, 2016). We treated age as a continuous variable in most analyses, using standard linear regression as well as generalized linear mixed-effects models (GLMMs) to explore relationships between output variables vs. dependent variables. Because electrophysiological data were collected from several neurons per subject, the subject was treated as a random effect blocking factor in relevant analyses (Grafen and Hails, 2002; Darian-Smith et al., 2013). For the GLMMs, marginal R² captured the variance explained by fixed factors; conditional R² captured the variance explained by the whole model including the random effects (Nakagawa and Schielzeth, 2013). In analyses with DRSTsp score (one data value per subject) vs. physiological variables, physiological variables were first averaged within each subject.

Age was treated as a categorical variable (young, middle-aged, and aged) when analyzing relationships firing rate vs. age (Figure 1E). Specifically, we used two-way repeated measures ANOVA (including subject as a random effect) to determine whether there were differences in FR vs. age group, with injected current as the repeated measure. Before performing the ANOVA, quantile-quantile plots were used to identify extreme outliers from normality for each age group and level of injected current, and then Shapiro-Wilk's test of normality was applied. Post-hoc analyses were conducted using Tukey-adjusted comparison of estimated marginal means. For all analyses, the significance level was α = 0.05.

Results

Rhesus Monkeys Show DRSTsp Impairment and Increased AP Firing Rates of dlPFC Pyramidal Neurons With Aging

We examined relationships between performance on the DRSTsp, age of the monkeys, and physiological variables assessed (Figure 1). Among the 39 subjects analyzed for this study, 36 of them completed behavioral testing on the DRSTsp (6 young, 18 middle-aged, and 12 aged). There was a significant decrease in performance on the DRSTsp with age [linear regression, DRST = −0.032^*Age + 2.97, R² = 0.18, F_{(1, 34)} = 7.43, p = 0.0101; Figure 1B] replicating past findings on impaired working memory with age in monkeys (Herndon et al., 1997; Moore et al., 2006, 2017; Wang et al., 2011). In vitro slices were prepared from the dlPFC of the same monkeys that were tested behaviorally for whole cell patch clamp recordings. Data from 324 neurons met our criteria for inclusion. Our statistical analysis used random effects modeling to control for the variation of physiological variables within monkeys (Materials and Methods). Input resistance (R_n) increased with aging (R_n = 4.02^*Age + 89.1 MΩ, marginal R² = 0.08, conditional R² = 0.20, p < 0.0001; Figure 1C), while resting membrane potential did not (p = 0.11). There was a significant positive relationship between input resistance and AP firing rate at all levels of current injection (p < 0.00016 and conditional R² ≥ 0.21 for each level, Figure 1D).

Next, we compared depolarizing current step-evoked AP firing rate in neurons from young, middle-aged, and aged subjects (Figure 1E). Before using the two-way repeated measures ANOVA, outliers from normality were removed for AP firing rate vs. current injection level for each age group. There was a significant effect of age group on firing rate [F_{(2, 22)} = 8.63, p = 0.0017] but no significant interaction between age group and injected current [F_{(2, 799)} = 2.76, p = 0.064]. Firing rates were significantly higher for the aged vs. young monkeys [Tukey-adjusted comparison of estimated marginal means, t₍₂₂₎ = 4.14, p = 0.0012]; firing rates for middle-aged monkeys did not differ from the other two groups [young vs. middle-aged, t₍₂₂₎ = 2.23, p = 0.072; middle-aged vs. aged, t₍₂₂₎ = 1.99, p = 0.139; Figure 1E]. Hereafter, we refer to this increased firing rate of current step-evoked APs as hyperexcitability of pyramidal neurons.

Given significant relationships between aging and both DRSTsp performance and the physiological variables, we next examined whether both aging and physiological variables might together affect DRSTsp performance. Since R_n and firing rate were significantly related at each injection level, we only included a single physiological variable in the analysis: firing rate in response to the 130 pA current injection. There was a significant negative linear relationship between DRSTsp span and firing rate [DRSTsp = −0.055^*FR + 2.95 Hz, R² = 0.19, F_{(1, 33)} = 7.71, p = 0.0090, Figure 1F]. To examine the effect of age and firing rate on DRSTsp span, we fitted a regression model that included age, firing rate, and their interaction ($R^2=0.26,F_{(3,31)}=3.68,p=0.0225$). The negative effect of the firing rate on DRSTsp was significant (p = 0.0441) and it did not depend on the age of the animal (age by firing rate interaction effect p = 0.1340). The effect of age alone was not significant either (p = 0.0904). These analyses lend support to the hypothesis that firing rate affects DRSTsp performance more than age does, though this support is somewhat tenuous for the following reasons. First, the delay between DRSTsp assessment and physiological recording of neurons from these valuable subjects can last several months to a year or more, and may mask a relationship between these variables. Second, some of these subjects were employed in longitudinal studies for other projects and tested on the DRSTsp a few times; we cannot rule out a practice effect in these subjects. Thus, we turned next to computational modeling to create a controlled environment for exploring how changes in AP firing rates (independent of a subject's age) might affect performance on spatial working memory tasks.

Model 1: The Delayed Response Task (DRT)

An Optimal Level of Pyramidal Neuronal Excitability Maximized Successful Performance in the DRT Model

Here, we used computational modeling to explore how these empirical results of physiological changes and a past study on reduced synapse counts with aging might together impact cognitive function in middle-aged and aged monkeys. We chose a continuous bump attractor network model which has been used by several previous research groups (Compte et al., 2000; Wang et al., 2011, 2013; Wimmer et al., 2014). The network model simulates an empirically measured behavior (the oculomotor DRT administered to macaque monkeys) through a simple but elegant model network that can be manipulated in specific ways (Figure 2; Materials and Methods). During the cue period of the DRT simulations, all excitatory neurons received an external input current, tuned to 0° direction, proportional to their orientation preference (Figure 2D). In the subsequent delay period, the main input to model neurons came through recurrent activity within the network. If persistent neural activity is maintained at the end of the delay period—in the form of a bump-shaped FR profile for a subset of the excitatory neurons—then the orientation angle corresponding to the center of the bump of activity at the end of the delay is assumed to represent the stimulus location encoded by the network (Figure 2E).

Varying model parameters created a distinct network configuration, and each simulation represented a model DRT trial. Successful DRT trials maintained TPA (and in particular TPA-S) until the end of the delay period (Figures 2E, 3A2). Failed DRT trials can occur in several ways, including under-excitation, a loss of TPA during the cue or delay periods (Figure 3A1); generating TPA early in the cue or delay periods but ending with full network over-excitation (Figure 3A3); and an over-excited network throughout the simulation (Figure 3A4).

An important feature of the bump attractor model we chose to use (Wimmer et al., 2014) is the ability to control the f-I curve of individual neurons by two simple parameters. Adjusting the parameter ν_ce was sufficient to fit the f-I curves of excitatory model neurons to the mean in vitro AP firing rates of pyramidal neurons reported above (Figure 3B and Equation 2). Young, middle aged, and aged f-I curves used ν_ce values of 5, 7, and 9 Hz respectively. These fits were accurate for low input current, but became more excitable than the empirical f-I curves for high inputs. However, our results did not depend on the specific shape of the f-I curves: similar results were found when using less excitable analogs of the firing rate function (Supplementary Figure 1). Synaptic weight parameters represented the number of excitatory and inhibitory synapses onto excitatory (G_EEa, G_EEn, and G_IE) and inhibitory (G_EIa, G_EIn, and G_II) neurons. Varying these synaptic weights allowed us to model published data showing an age-related reduction in the excitatory and inhibitory synapses with aging in the neuropil of layer 2/3 of rhesus monkey dlPFC (Peters et al., 2008). In a series of simulations, we kept ν_ce fixed at distinct values and then used the Latin Hypercube Sampling (Materials and Methods) to choose points throughout the parameter space generated by the six synaptic weights. Each point was simulated seven times, as if administering multiple DRT trials to each virtual monkey. We examined whether points across the parameter space maintained TPA throughout the delay period, or had a different outcome (Figure 3C).

Increasing the excitability of model excitatory neurons strongly affected the ability of the DRT model networks to maintain TPA. The 3D plots in Figure 3C summarize DRT model results for 4,200 points across the parameter space of synaptic weights as ν_ce increased. Blue circles represent “under-excited” networks, where no TPA was detected at the end of the delay period (e.g., Figure 3A1). Green circles represent networks that maintained TPA for the entire delay period (e.g., Figure 3A2)—corresponding to “virtual monkeys” that encoded some stimulus location. Remaining circles represent partially over-excited networks (yellow; e.g., Figure 3A3), which started out with TPA that led to firing in all neurons, or were already over-excited networks during the fixation period before the cue ever appeared (red; e.g., Figure 3A4). Below we focus on the green points, which maintained TPA throughout the delay period.

The greatest likelihood of maintaining TPA, and in particular TPA-S, across the synaptic weight parameter space occurred for intermediate values of ν_ce–neither too low nor too high (Figure 3D). The optimal value among those we tested was ν_ce = 5 Hz. Values of ν_ce below this optimum led to many networks across parameter space that were under-excited; above the optimum many networks were over-excited (Figure 3C). In either case, the networks were much less likely to encode the stimulus successfully until the end of the delay period, and this reduced likelihood was visible even in our “middle-aged” networks (ν_ce = 7 Hz). Thus, if the excitability of single pyramidal neurons increased above the optimal value without complementary changes in other parameters, the network's ability to maintain TPA, in particular to the stimulus location—and successful working memory performance—was impaired. These findings also held with the less excitable form of the f-I curves fit to the data (Supplementary Figures 1b,d).

Synaptic Changes Partially Compensated for Increased Pyramidal Neuron Excitability in Successful DRT Model Networks

Points across parameter space that did maintain TPA-S as ν_ce increased had lower values of the excitatory synapse weights for excitatory neurons (G_EEa and G_EEn) and higher values of the inhibitory synapse weights (G_IE) for these neurons (see Figure 3F; standard error bars lie beneath the symbols). Figure 3E shows the mean values of each synaptic weight for all sampled points maintaining TPA-S as ν_ce increased. There was no clear relationship between ν_ce and the synaptic weight parameters onto interneurons (G_EIa, G_EIn, and G_II). However, excitatory parameters G_EEa and G_EEn decreased noticeably as ν_ce increased before leveling off after the middle-aged ν_ce value, while G_IE increased in a nearly linear fashion. To maintain TPA in the DRT model (particularly when tuned to the stimulus location), the increased excitability of excitatory neurons was compensated by altering synaptic weights onto those same neurons: with less synaptic excitation and more synaptic inhibition. Yet despite the compensatory synaptic weight changes, fewer networks were able to maintain TPA-S as ν_ce increased (Figure 3D). See also Supplementary Figure 1c for the generalized f-I curve.

These LHS simulations revealed that network function can be maintained somewhat as parameters vary freely to compensate pyramidal neuron hyperexcitability with aging, but the predicted compensation is not consistent with the empirically observed loss of both excitatory and inhibitory synapses on pyramidal neurons with aging. To explore how concomitant changes to firing rates and synapses might affect network function, we first considered all networks maintaining TPA-S when ν_ce = 5 Hz as corresponding to a simulated cohort of cognitively healthy young monkeys (the 679 model networks with ν_ce = 5 Hz in Figure 3D, called the “young model cohort”). We then examined the performance of the model networks in this cohort on the DRT after perturbing parameters in proportions that matched empirically observed changes with aging (Figure 4). First is the hyperexcitability for pyramidal neurons in vitro described above and in our past work (Chang et al., 2005; Coskren et al., 2015). Second, Peters et al. (2008) reported a ~10% loss of excitatory and inhibitory synapses in the rhesus dlPFC neuropil by middle-age, and a ~30% loss of both synapse types in aged subjects. This finding is consistent with other studies of synapse and spine loss with aging in the rhesus dlPFC (Hof et al., 2002; Duan et al., 2003; Young et al., 2014).

Accordingly, we induced “middle-aged” and “aged” conditions for pyramidal neuron excitability in the young cohort by increasing ν_ce to 7 and 9 Hz, respectively. To model synapse loss, we reduced G_EEa, G_EEn, and G_IE in the young model cohort by 10 and 30%, respectively to create either a middle-aged or aged synapse condition. We then induced the two conditions separately and together (forming the “observed data conditions”), and determined whether each transformed model still maintained TPA-S during DRT simulations (Figures 4A–C). A similar synapse loss condition for excitatory synapses was modeled by Wang et al. (2011) to explain an age-related decrease in firing rates during in vivo recordings of the DELAY neurons for monkeys successfully performing the DRT. Our study builds on Wang et al. (2011) by adding inhibitory synapse loss and f-I curve changes with aging, and examining how these alterations affect task performance for a range of parameter values.

DRT performance of the young model cohort was much more sensitive to the excitability aging condition than to the synapse condition (Figures 4A,B). Even under the strongest (aged) synapse condition, over half the networks still maintained TPA-S (Figure 4A). In contrast, under the excitability condition most of the networks had already lost their capacity for TPA-S by middle age (Figure 4B). In the parameter space exploration above, pyramidal neuron hyperexcitability was partially compensated by decreasing excitatory synapse weights and increasing inhibitory ones. Thus, while Peters et al. (2008) found that inhibitory synapse counts actually decrease with aging, we might still expect a partial recovery of function when the young cohort underwent the excitability and synapse aging conditions simultaneously (“observed data conditions,” Figure 4C). However, this was not the case, as almost none of the middle aged and aged models maintained TPA-S after these perturbations. This does not seem to match reality, since Wang et al. (2011) noted no difficulty in identifying middle-aged and aged monkeys who perform the DRT successfully, and impairment on the more difficult DRSTsp is noticeable but not severe by middle age (Figure 1B).

Since external inputs in vivo are mediated by synapses whose numbers likely also decrease with aging as in Peters et al. (2008), we tested whether reducing the input stimulus current affected the results for the conditions applied in Figures 4A,C (Supplementary Figures 2a,b). We reduced the stimulus scaling factor, I_st, by 10 and 30% for middle aged and aged networks, respectively, in addition to applying the synapse condition (Supplementary Figure 2a) and both observed data conditions (Supplementary Figure 2b). The stimulus reduction did cause a few more of the young cohort models to fail in the middle age and aged conditions, but the overall trend of the results and the number of impaired networks was similar. In particular, DRT performance of the young model cohort was still much more sensitive to the excitability aging condition than to the synapse condition with stimulus reduction (compare Figure 4B and Supplementary Figure 2a).

Recent studies have shown that including short-term synaptic facilitation in discrete and continuous attractor models can reduce the fine-tuning needed to maintain performance (Mongillo et al., 2008; Itskov et al., 2011; Rolls et al., 2013). We added short-term synaptic facilitation to our middle-aged and aged DRT model networks (same parameter values in both groups), as a baseline exploration of how this mechanism might affect DRT performance with aging. The results were qualitatively similar to those without facilitation, except the model cohort networks were affected less by the aging perturbations (Figures 4D–F and Supplementary Figures 2c,d including stimulus reduction). Under the synapse condition, almost all networks maintained TPA-S both for middle-aged and aged simulations, with even a slight increase from middle-aged to aged when the input stimulus was constant. Networks were still more sensitive to the excitability aging condition than to the synapse condition (Figure 4D and Supplementary Figure 2c vs. Figure 4E), but under the excitability condition nearly all the middle-aged and about half the aged networks maintained TPA-S. Applying both observed aging conditions together with facilitation led to more networks maintaining TPA-S (Figure 4F and Supplementary Figure 2d), suggesting that synaptic facilitation can boost DRT performance.

Examining Firing Rates During the Simulated DRT

Wang et al. (2011) examined in vivo firing rates of area 46 (dlPFC) pyramidal neurons in rhesus monkeys across the adult age span while they performed the DRT. For monkeys performing the task well, they found the firing rate of DELAY neurons during the cue and delay periods decreased with aging, both for neurons whose preferred and anti-preferred directions matched the cue presented. In our sampling of parameter space (Figure 3) we found that DRT model networks maintaining TPA-S—representing monkeys that perform the DRT successfully—had lower firing rates during both the cue and delay periods as ν_ce increased (Figure 5A, left). This was particularly true for neurons whose preferred direction matched the stimulus, with a substantial decrease from young and middle-aged to aged model networks. This seemed counterintuitive at first, since increasing ν_ce makes individual excitatory neurons more excitable in response to a given input current (Figure 3B). Indeed, hyperexcitability of individual neurons often did lead to over-excitation throughout the network (yellow and red points in Figure 3C). However, the networks that maintained TPA-S as ν_ce increased compensated for the increased excitability with a reduction in G_EEa and G_EEn and increase in G_IE (Figure 3F), leading to lower firing rates. Cue period FR in our model was driven much more by the stimulus current than by synaptic weights, and decreased less across the age groups than the delay period FR did. As such, reducing the stimulus current in the sampled middle aged and aged networks as in Supplementary Figure 2 was necessary to match the Wang et al. (2011) cue period results (Materials and Methods). As a whole, these parameter space explorations show a way to unify our in vitro data and in vivo experiments, provided that the total input current to pyramidal neurons during the DRT, I_E, is lower in middle-aged and aged monkeys than in young monkeys (Figure 5A, right).

We also examined firing rates in the model cohort simulations perturbed by the observed data conditions with and without facilitation (Figures 4C,F). The main distinction between the parameter space sampling simulations and the cohort simulations was how synaptic inhibition onto the excitatory neurons, G_IE, varied as the excitability parameter (ν_ce) increased: G_IE increased in the parameter space explorations but decreased in the cohort simulations due to observed inhibitory synapse loss with aging. Indeed, the weakened inhibitory connections in the middle-aged and aged cohorts led to a FR increase during the cue and delay periods of those networks compared to the young cohort (Figure 5B, left), opposing the Wang et al. (2011) results. The increase was even greater when synaptic facilitation was included (Figure 5B, right).

Our young cohort simulations did not match data from the in vitro physiology, synapse counts, and in vivo physiology findings simultaneously, suggesting that one or more age-related changes to the system may be missing from our current model. We simulate one such possibility here, to demonstrate how the current model can help us test hypotheses and incorporate new data as they become available. Specifically, our parameter space explorations showed that an overall increase in inhibition onto the excitatory neurons can compensate for increased firing rates of excitatory neurons. Thus, we examined how network function would change if the observed data conditions were accompanied by a comparable age-related increase of inhibitory neuron firing rates—an empirical question whose answer is not yet known. In models without synaptic facilitation, more middle-aged and aged networks maintained TPA-S when inhibitory neuron firing rates increased with aging (Figure 4C vs. Figure 5C, top left). Firing rates of excitatory neurons in those simulations also decreased with aging as observed in vivo, consistent with the Wang et al. data (Wang et al., 2011; Figure 5C, bottom left). However, results differed for models that included synaptic facilitation: fewer middle-aged and aged networks maintained TPA-S when inhibitory neuron firing increased (Figure 4F vs. Figure 5C, top right), while (as before) firing rates of excitatory neurons during the DRT increased with aging (Figure 5C, bottom right). These results suggest that our DRT model might need further refinement to reflect reality better—for example, by also incorporating synaptic depression (Mi et al., 2017) or changes in synaptic plasticity with aging. Regardless, our model allows us to quantify how much the known age-related changes at the single cell level affect cognitive performance, and to predict the extent to which other relevant processes might contribute to working memory impairment with aging.

Model 2: Working Memory Retention in a DRSTsp-Like Task

Pyramidal Neuron Hyperexcitability Led to Reduced Memory Capacity in Simulated DRSTsp Networks

These DRT modeling results add to the wide range of literature using the bump attractor model to study spatial working memory, but the DRT was not administered to the monkeys from this study. Spatial working memory capacity in these monkeys was assessed with the more complex DRSTsp. To see what predictions our model might have for these data, we extended our DRT model to simulate memory retention in a task with three important similarities to the DRSTsp: (1) the network encoded several stimuli simultaneously, with stimuli that are (2) spatially distinct; and (3) presented to the network successively over time. The changes made to extend the DRT model simulations to a task like the DRSTsp are shown in Figure 6A, and a simulated cognitive score (DRST_sim) was assigned to the outcome of each simulation. Examples of model networks with a DRST_sim of 4 (perfect performance) down to 1 are shown in Figure 6B. The second network achieved a DRST_sim of 3 because it encoded stimuli 1, 2, and 3 correctly but lost persistent activity corresponding to stimulus 3 in the subsequent delay period (red arrowhead). Likewise, the third network (DRST_sim = 2) lost persistent activity corresponding to stimulus 2 before stimulus 3 was applied; the final network (DRST_sim = 1) encoded stimulus 1 properly but not stimulus 2.

We examined how DRST_sim varied across the parameter space of synaptic weights for the young, middle-aged and aged network models (ν_ce = 5, 7, and 9 Hz; Figure 6C). Networks achieving a DRST_sim of 1 (worst), 2, 3, and 4 (best) are shown, respectively, as blue, yellow, red, and green. As in the DRT, the value of ν_ce strongly affected the ability of the DRSTsp model networks to maintain persistent activity tuned to the location of the different stimuli as the number of delay periods and unique stimuli increased. There were many fewer green points in the parameter regime explored, and more blue points, as ν_ce increased (As ν_ce increased, networks with the best DRST_sim performance had lower values of G_EEa, G_EEn, and G_IE than those shown here.). Figure 6D summarizes these counts within our regime, showing a clear leftward shift in the distribution of DRST_sim across the parameter space (poorer performance) as ν_ce increased.

Figure 6E shows how the model networks of Figure 6B were affected by increasing ν_ce from 5 (in Figure 6B) to 7 and 9 Hz. In general, as ν_ce increased, the initial 0° stimulus was encoded by a stronger bump (higher firing rate) for the entire DRSTsp simulation, and subsequent stimuli encoded by comparatively weaker bumps. This led to an under-excitation of the network during the second and third delay periods, losing the ability to maintain later stimuli. An example can be seen in the second column of Figures 6B,E. For the young monkey case (Figure 6B, ν_ce = 5 Hz), the maximum firing rate of the first bump was around 75 Hz; this and the width of the first bump decreased slightly as the simulation continued. Also, as noted above, persistent activity for stimulus 3 was lost during the third delay period. For the middle-aged case (Figure 6E top row, ν_ce = 7 Hz), the maximum firing rate of the first bump maintained a maximum firing rate near 150 Hz and a similar width throughout the simulation, but the second bump was lost early in the subsequent delay period. This trend was even stronger for the aged case (Figure 6E bottom row, ν_ce = 9 Hz), where the first bump remained wide with a maximum firing rate of 250 Hz throughout the simulation but none of the subsequent stimuli were encoded.

To see how synaptic parameters compensated to maintain DRSTsp performance as ν_ce increased, we looked at the mean values of the synaptic weights for all networks with perfect task performance (DRST_sim = 4). Values of both the excitatory and inhibitory synaptic weights for excitatory neurons were lower as ν_ce increased (Figure 6F top; standard error bars lie beneath the symbols). Excitatory synaptic weights for inhibitory neurons decreased and inhibitory synaptic weights increased as ν_ce increased (Figure 6F bottom; standard error bars lie beneath the symbols). Overall, for the DRSTsp model, both panels showed that synaptic excitation as well as inhibition to excitatory neurons decreased to partially compensate for the increased excitability of individual excitatory neurons.

“Aging” the DRSTsp Model Networks Led to Working Memory Impairment as Observed Empirically

As with the DRT model, we defined a simulated “young monkey cohort” for the DRSTsp model and “aged” the model networks in this cohort by applying increased excitability and synaptic loss conditions (Figure 7), choosing 40 levels of parameter perturbations throughout the simulated age span. We induced the two observed data conditions separately and together—and to match Figure 5C top right, also increased the excitability of inhibitory neurons (Since reducing the stimulus current had a minimal effect on performance during the synapse condition perturbations of the DRT, we kept the stimulus current constant across all DRSTsp simulations.). We then computed the corresponding DRST_sim for each transformed member of the young cohort. To summarize the change in performance across the simulated age span, we calculated the relative change in the mean DRST_sim from the young to the aged models under each perturbation condition. This facilitated comparison between performance on the simulated task (four stimuli spaced evenly around a ring), and the real DRSTsp (18 possible stimuli, organized spatially into three rows of six wells each).

All four cases showed a strong linear relationship between DRST_sim and parameter perturbations. As in the DRT model above, the DRSTsp performance of the young model cohort was much more sensitive to the excitability condition (Figure 7B) than to the synapse condition (Figure 7A). While the synapse condition showed a 4% improvement with aging for the DRSTsp model (R² = 0.98, p < 0.0001), the excitability-only condition showed a 45% decrease in DRST_sim with aging (R² = 0.98, p < 0.0001). Applying the two observed conditions together led to a 14% impairment with aging (Figure 7C; R² = 0.82, p < 0.0001); applying the observed data conditions plus increasing the excitability of inhibitory neurons led to a 17% impairment (Figure 7D; R² = 0.81, p < 0.0001). Thus, both cases of applying the observed data conditions were comparable to the 20% impairment exhibited by the empirical data (Figure 1B). These results are also consistent with Figure 6F: the empirically-observed decrease in synaptic excitation and inhibition to excitatory neurons partially compensated the increased excitability of the excitatory neurons, so that DRST_sim impairment with aging was less severe. The model also suggests that hyperexcitability of inhibitory neurons would impact performance on the DRT more than it would the DRSTsp.

Discussion

This study used computational modeling to explore how single pyramidal neuron hyperexcitability and excitatory and inhibitory synapse loss observed with aging in layer 3 of the rhesus monkey dlPFC might contribute to spatial working memory impairment. First, we presented new empirical data from rhesus monkeys across the adult age span, consistent with past studies, showing that working memory impairment and increased FRs of neurons occur in some middle-aged monkeys. This replicated previous studies showing a decline in cognitive performance from young to aged monkeys (Herndon et al., 1997; Moore et al., 2003, 2006; Konar et al., 2016; Motley et al., 2018), and a concomitant increase in R_n and FR of dlPFC pyramidal neurons (Chang et al., 2005; Coskren et al., 2015). These data motivated the modeling, and were used to constrain a bump attractor network model of the DRT (Wimmer et al., 2014) and our network model of working memory retention in the DRSTsp. These models predict that the observed age-related firing rate increase and synapse loss in rhesus dlPFC (Peters et al., 2008) are sufficient to induce significant impairment of spatial working memory function.

There is evidence that persistent neural activity during delay periods of a working memory task is important for encoding the memory of a stimulus (reviewed in Constantinidis et al., 2018). When age-related alterations to the excitability properties of neurons comprising the working memory network occur, the network's capacity for maintaining persistent activity—keeping activity tuned to the original stimulus—is impaired (Wang et al., 2011). To provide insight into the relationship between neuronal firing rate and synaptic changes seen in aging and working memory decline, we examined the parameter space of a bump attractor model of the DRT, then extended the model to simulate how aging might affect memory retention of multiple stimuli in a more complex task—the DRSTsp. The region of parameter space that simulated the DRT successfully was largest when the pyramidal neuron excitability parameter (ν_ce) was tuned to fit the young FR data. The network model had its own built-in compensatory mechanisms: when pyramidal neuron excitability increased, the region with successful DRT models shifted to lower excitatory and higher inhibitory synapse weights. This tradeoff allowed networks to control the overall excitatory:inhibitory balance needed to maintain TPA. These compensations were consistent with empirical data showing an age-related decrease in excitatory synapses, but not consistent with the reported decrease in inhibitory synapses (Peters et al., 2008). However, this compensatory mechanism led to lower firing rates during the delay period, as well as during the cue period when the stimulus current was also reduced, unifying our in vitro observations with the in vivo data (Wang et al., 2011). We predict that the network compensates higher FR in individual pyramidal neurons with aging by lower total synaptic input current to each neuron, and that the stimulus current could also be decreasing with aging.

To test how hyperexcitability of pyramidal neurons and loss of excitatory and inhibitory synapses with aging affected network output, we ran a second set of simulations: “aging” model neurons by decreasing both excitatory and inhibitory synaptic strengths and increasing the FR of the excitatory neurons. DRT model performance (represented by the number of networks performing the task successfully) was substantially impaired in middle-aged and aged simulations, more than reported in a sample of rhesus monkeys (Wang et al., 2011), and the FR increase affected DRT performance much more than the loss of synapses did. Adding synaptic facilitation restored successful function to many middle aged and aged DRT models. However, only one of the aging cohort scenarios we tested captured the hyperexcitable pyramidal neuron f-I curves observed in vitro, dlPFC synapse loss observed with electron microscopy, and the in vivo FR decrease during the DRT cue and delay periods from Wang et al. (2011): the one without synaptic facilitation but added hyperexcitability of inhibitory neurons to restore excitatory:inhibitory balance. One next step is to examine the excitability of dlPFC interneurons with aging, but our simulations raise several other questions about the DRT. Is there an age-related reduction in overall input to dlPFC during the cue period? How central a role does synaptic facilitation play during the task? How do age-related changes in synaptic plasticity, discussed below, affect DRT performance? Might neuromodulation (Arnsten et al., 2012; Davis et al., 2017; Wang et al., 2019) during working memory task performance reduce the pyramidal neuron hyperexcitability observed in vitro? Computational models provide an essential means for testing and refining hypotheses as future data become available.

Finally, we extended the DRT bump attractor model to create a model of working memory retention in a task like the DRSTsp. Other network models of working memory have simulated several stimuli at once (Edin et al., 2009; Wei et al., 2012; Mi et al., 2017). Yet, while the DRSTsp has been used for many years to evaluate spatial working memory capacity in rhesus monkeys (Herndon et al., 1997; Luebke et al., 2004; Chang et al., 2005; Moore et al., 2005, 2017; Peters et al., 2008; Luebke and Amatrudo, 2012)—including in the present study—to our knowledge this is the first computational model of this task. Increasing pyramidal neuron FR strengthened the response to the initial spatial cue and reduced the ability of the DRSTsp model networks to encode memories of subsequent spatial cues, so that middle-aged and aged networks had lower DRST_sim than young simulations. As with the DRT model, “aging” the young DRSTsp model cohort led to severe impairment when increasing pyramidal neuron excitability alone, but not synapse loss alone. The maintenance of the DRSTsp model performance under synapse loss alone is consistent with Peters et al. (2008), which found no correlation between the numerical density of synapses and DRSTsp impairment. Simultaneously varying synapse loss and pyramidal neuron hyperexcitability within ranges observed with aging, both with and without hyperexcitability of interneurons, the amount of DRST_sim impairment was similar to that seen in DRSTsp scores across the adult rhesus monkey life span. Thus, our model predicts that pyramidal neuron hyperexcitability and synapse loss may be sufficient to explain empirically observed levels of spatial working memory impairment. Future experiments will examine whether the amount of synapse loss reported in Peters et al. (2008), is also seen in the dlPFC of a subset of monkeys from this study, and whether there are correlations between those data, AP firing rates, and spontaneous postsynaptic currents.

There is an open debate about the role of synaptic mechanisms for working memory (Constantinidis et al., 2018; Lundqvist et al., 2018). While delay activity is recognized as important to working memory, a recent perspective posited that sparse spiking activity and synaptic plasticity between spike times—rather than asynchronous persistent activity—might be the mechanism actually responsible for memory maintenance (Lundqvist et al., 2018). This view is supported by several modeling studies (Sandberg et al., 2003; Mongillo et al., 2008; Lundqvist et al., 2011). Our models lend further support to an essential role for both persistent activity and synaptic facilitation for working memory tasks involving multiple cues. Synaptic plasticity (both facilitation and depression) help overcome limitations of continuous attractor models, such as filtering out distractors (Compte et al., 2000) and remembering multiple items (Itskov et al., 2011; Wei et al., 2012; Hansel and Mato, 2013; Rolls et al., 2013; Mi et al., 2017) as in our DRSTsp model. Bump attractor models of visuo-spatial working memory have shown that short-term synaptic facilitation in recurrent connections of excitatory neurons slows down the drift of the bump often seen in continuous attractor models (Itskov et al., 2011). Hansel and Mato (2013) showed that short-term synaptic facilitation might lead to the direction selectivity of synaptic weights and network bistability in a DRT model. In vitro experiments have shown synaptic facilitation in PFC pyramidal neurons (Hempel et al., 2000; Wang et al., 2006). Several other studies have shown age-related changes in synaptic plasticity (not modeled here), some correlating these changes with learning and memory deficits (Norris et al., 1998a,b; Bach et al., 1999; Rosenzweig and Barnes, 2003; Thibault et al., 2007, 2013; Gant et al., 2011, 2014, 2015).

This network-level study represents an extension of our past modeling studies at the single-neuron level on parameters affecting firing properties of individual pyramidal neurons with aging (Coskren et al., 2015; Luebke et al., 2015; Rumbell et al., 2016), in different cortical regions (Amatrudo et al., 2012), and in neurodegeneration (Goodliffe et al., 2018). The elegance of these network models lies in their ability to represent empirically observed changes in FR and synaptic weights with aging as perturbations to a small number of network parameters. Based on our findings with two f-I curves, replacing our firing rate model neurons with a spiking neuron model (Stein, 1967; Hodgkin and Huxley, 1990; Izhikevich, 2003; Teeter et al., 2018) in our DRT and DRSTsp networks should yield qualitatively similar results. Our future network modeling will incorporate age-related changes in passive and active channel parameters predicted in our and other past studies (Wang et al., 2011; Coskren et al., 2015; Rumbell et al., 2016), and other empirically observed changes in aging dlPFC including increased slow after hyperpolarization current (Luebke and Amatrudo, 2012); white matter changes (review: Luebke et al., 2010; Kubicki et al., 2019); and dysregulation of calcium homeostasis (Foster, 2007; Toescu and Vreugdenhil, 2010; Oliveira and Bading, 2011). Age-related changes to synaptic strength and neurotransmitters reduced the probability of short-term memory recall from a discrete attractor network model (Rolls and Deco, 2015); similar perturbations in a spiking-network version of our DRT and DRSTsp models would be extremely insightful.

The main limitation of our current DRSTsp model—assuming that the network correctly identifies the novel stimulus among several retained—is also a limitation of the classical bump attractor as a whole. Experiments have shown that decision-making and attentional tasks involve connections between dlPFC and other cortical regions including posterior parietal cortex, anterior cingulate cortex, and secondary motor areas (Medalla and Barbas, 2010; Oemisch et al., 2015; Constantinidis and Qi, 2018; Marvel et al., 2019; Tsunada et al., 2019). These connections may represent different aspects of a task like DRSTsp, such as the monkey's internal model of itself and the task; its level of attention; signals of reward and error; and motor traces that might correspond to rehearsal strategies between active stimulus stages of the DRSTsp. Recent theoretical studies have shown that model neurons with mixed selectivity, responding to sensory stimuli as well as internal states of the system, can perform complex cognitive tasks (Maass et al., 2002; Salinas, 2004; Rigotti et al., 2010; Chaisangmongkon et al., 2017). Such units could be well-suited for extending the DRSTsp model introduced here to simulate choice of the novel spatial cue and related internal signals, plus the full 3 × 6 array of wells presented to the monkeys during the DRSTsp. There is evidence that aging affects information flow between cortical regions, which can contribute to working memory impairment (Engle et al., 2016; Lee et al., 2016; Proskovec et al., 2016; King et al., 2018; Koen et al., 2019). Any such age effects must be examined in concert with the changes to dlPFC synapses and neuronal excitability, to predict the extent to which these concomitant changes compound vs. partially compensate each other.

In summary, our computational models of two spatial memory tasks predicted that empirically observed changes in pyramidal neuron excitability and synapse loss lead to spatial working memory impairment with aging. Computational studies such as this provide a means to evaluate connections between physiological and behavioral—and in vitro vs. in vivo—experiments, extending the impact provided by animal studies.

Data Availability Statement

Our simulation code is publicly available at: http://modeldb.yale.edu, Accession #256610.

Ethics Statement

The animal study was reviewed and approved by Boston University Institutional Animal Care and Use Committee (IACUC).

Author Contributions

SI developed and implemented the models, analyzed the results, carried out the statistical analysis of the data, and wrote the manuscript. JL performed the experiments, wrote the manuscript, and designed the study. WC performed the experiments. DD designed the statistical analysis of the data. CW carried out the statistical analysis of the data, wrote the paper, and designed the study.

Funding

This work was supported by NIH/NIA grants R01 AG059028 and R56 AG049870.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Drs. Tara Moore and Doug Rosene for the behavioral data and brain tissue from which recordings were made. We thank Drs. Maria Medalla, Patrick Hof, and Patrick Coskren for productive discussions. This manuscript has been released as a Pre-Print at https://www.biorxiv.org/content/10.1101/745901v2 (Ibañez et al., 2019).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncom.2019.00089/full#supplementary-material

References

Abeles, M. (1991). Corticonics: Neural Circuits of the Cerebral Cortex. New York, NY: Cambridge University Press. doi: 10.1017/CBO9780511574566

CrossRef Full Text | Google Scholar

Albert, M. (1993). Neuropsychological and neurophysiological changes in healthy adult humans across the age range. Neurobiol. Aging 14, 623–625. doi: 10.1016/0197-4580(93)90049-H

PubMed Abstract | CrossRef Full Text | Google Scholar

Amatrudo, J. M., Weaver, C. M., Crimins, J. L., Hof, P. R., Rosene, D. L., and Luebke, J. I. (2012). Influence of highly distinctive structural properties on the excitability of pyramidal neurons in monkey visual and prefrontal cortices. J. Neurosci. 32, 13644–13660. doi: 10.1523/JNEUROSCI.2581-12.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Amit, D. J., Bernacchia, A., and Yakovlev, V. (2003). Multiple-object working memory - a model for behavioral performance. Cereb. Cortex 13, 435–443. doi: 10.1093/cercor/13.5.435

PubMed Abstract | CrossRef Full Text | Google Scholar

Arnsten, A. F. T., Wang, M. J., and Paspalas, C. D. (2012). Neuromodulation of thought: flexibilities and vulnerabilities in prefrontal cortical network synapses. Neuron 76, 223–239. doi: 10.1016/j.neuron.2012.08.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Bach, M. E., Barad, M., Son, H., Zhuo, M., Lu, Y. F., Shih, R., et al. (1999). Age-related defects in spatial memory are correlated with defects in the late phase of hippocampal long-term potentiation in vitro and are attenuated by drugs that enhance the cAMP signaling pathway. Proc. Natl. Acad. Sci. U.S.A. 96, 5280–5285. doi: 10.1073/pnas.96.9.5280

PubMed Abstract | CrossRef Full Text | Google Scholar

Barak, O., and Tsodyks, M. (2014). Working models of working memory. Curr. Opin. Neurobiol. 25, 20–24. doi: 10.1016/j.conb.2013.10.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowley, M. P., Cabral, H., Rosene, D. L., and Peters, A. (2010). Age changes in myelinated nerve fibers of the cingulate bundle and corpus callosum in the rhesus monkey. J. Comp. Neurol. 518, 3046–3064. doi: 10.1002/cne.22379

PubMed Abstract | CrossRef Full Text | Google Scholar

Braitenberg, V., and Schütz, A. (1991). Anatomy of the Cortex. Berlin: Springer-Verlag. doi: 10.1007/978-3-662-02728-8

CrossRef Full Text | Google Scholar

Brunel, N. (2003). Dynamics and plasticity of stimulus-selective persistent activity in cortical network models. Cereb. Cortex 13, 1151–1161. doi: 10.1093/cercor/bhg096

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaisangmongkon, W., Swaminathan, S. K., Freedman, D. J., and Wang, X.-J. (2017). Computing by robust transience: how the fronto-parietal network performs sequential, category-based decisions. Neuron 93, 1504–1517. doi: 10.1016/j.neuron.2017.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, Y.-M., Rosene, D. L., Killiany, R. J., Mangiamele, L. A., and Luebke, J. I. (2005). Increased action potential firing rates of layer 2/3 pyramidal cells in the prefrontal cortex are significantly related to cognitive performance in aged monkeys. Cereb. Cortex 15, 409–418. doi: 10.1093/cercor/bhh144

PubMed Abstract | CrossRef Full Text | Google Scholar

Cleeland, C., Pipingas, A., Scholey, A., and White, D. (2019). Neurochemical changes in the aging brain: a systematic review. Neurosci. Biobehav. Rev. 98, 306–319. doi: 10.1016/j.neubiorev.2019.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Compte, A., Brunel, N., Goldman-Rakic, P. S., and Wang, X. J. (2000). Synaptic mechanisms and network dynamics underlying spatial working memory in a cortical network model. Cereb. Cortex 10, 910–923. doi: 10.1093/cercor/10.9.910

PubMed Abstract | CrossRef Full Text | Google Scholar

Constantinidis, C., Funahashi, S., Lee, D., Murray, J. D., Qi, X.-L., Wang, M., et al. (2018). Persistent spiking activity underlies working memory. J. Neurosci. 38, 7020–7028. doi: 10.1523/JNEUROSCI.2486-17.2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Constantinidis, C., and Qi, X.-L. (2018). Representation of spatial and feature information in the monkey dorsal and ventral prefrontal cortex. Front. Integr. Neurosci. 12:31. doi: 10.3389/fnint.2018.00031

PubMed Abstract | CrossRef Full Text | Google Scholar

Coskren, P. J., Luebke, J. I., Kabaso, D., Wearne, S. L., Yadav, A., Rumbell, T., et al. (2015). Functional consequences of age-related morphologic changes to pyramidal neurons of the rhesus monkey prefrontal cortex. J. Comput. Neurosci. 38, 263–283. doi: 10.1007/s10827-014-0541-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Darian-Smith, C., Lilak, A., and Alarcón, C. (2013). Corticospinal sprouting occurs selectively following dorsal rhizotomy in the macaque monkey. J. Comp. Neurol. 521, 2359–2372. doi: 10.1002/cne.23289

PubMed Abstract | CrossRef Full Text | Google Scholar

Davis, S. W., Luber, B., Murphy, D. L. K., Lisanby, S. H., and Cabeza, R. (2017). Frequency-specific neuromodulation of local and distant connectivity in aging and episodic memory function. Hum. Brain Mapp. 38, 5987–6004. doi: 10.1002/hbm.23803

PubMed Abstract | CrossRef Full Text | Google Scholar

Duan, H., Wearne, S. L., Rocher, A. B., Macedo, A., Morrison, J. H., and Hof, P. R. (2003). Age-related dendritic and spine changes in corticocortically projecting neurons in macaque monkeys. Cereb. Cortex 13, 950–961. doi: 10.1093/cercor/13.9.950

PubMed Abstract | CrossRef Full Text | Google Scholar

Edin, F., Klingberg, T., Johansson, P., McNab, F., Tegnér, J., and Compte, A. (2009). Mechanism for top-down control of working memory capacity. Proc. Natl. Acad. Sci. U.S.A. 106, 6802–6807. doi: 10.1073/pnas.0901894106

PubMed Abstract | CrossRef Full Text | Google Scholar

Engle, J. R., Machado, C. J., Permenter, M. R., Vogt, J. A., Maurer, A. P., Bulleri, A. M., et al. (2016). Network patterns associated with navigation behaviors are altered in aged nonhuman primates. J. Neurosci. 36, 12217–12227. doi: 10.1523/JNEUROSCI.4116-15.2016

PubMed Abstract | CrossRef Full Text | Google Scholar

Fisk, J. E., and Sharp, C. A. (2004). Age-related impairment in executive functioning: updating, inhibition, shifting, and access. J. Clin. Exp. Neuropsychol. 26, 874–890. doi: 10.1080/13803390490510680

PubMed Abstract | CrossRef Full Text | Google Scholar

Foster, T. C. (2007). Calcium homeostasis and modulation of synaptic plasticity in the aged brain. Aging Cell 6, 319–325. doi: 10.1111/j.1474-9726.2007.00283.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Funahashi, S. (2017). Working memory in the prefrontal cortex. Brain Sci. 7:49. doi: 10.3390/brainsci7050049

CrossRef Full Text | Google Scholar

Gant, J. C., Blalock, E. M., Chen, K.-C., Kadish, I., Porter, N. M., Norris, C. M., et al. (2014). FK506-binding protein 1b/12.6: a key to aging-related hippocampal Ca2+ dysregulation? Eur. J. Pharmacol. 739, 74–82. doi: 10.1016/j.ejphar.2013.10.070

PubMed Abstract | CrossRef Full Text | Google Scholar

Gant, J. C., Chen, K.-C., Kadish, I., Blalock, E. M., Thibault, O., Porter, N. M., et al. (2015). Reversal of aging-related neuronal Ca2+ dysregulation and cognitive impairment by delivery of a transgene encoding FK506-binding protein 12.6/1b to the hippocampus. J. Neurosci. 35, 10878–10887. doi: 10.1523/JNEUROSCI.1248-15.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Gant, J. C., Chen, K.-C., Norris, C. M., Kadish, I., Thibault, O., Blalock, E. M., et al. (2011). Disrupting function of FK506-binding protein 1b/12.6 induces the Ca²⁺-dysregulation aging phenotype in hippocampal neurons. J. Neurosci. 31, 1693–1703. doi: 10.1523/JNEUROSCI.4805-10.2011

CrossRef Full Text | Google Scholar

Goldman-Rakic, P. S. (1995). Cellular basis of working memory. Neuron 14, 477–485. doi: 10.1016/0896-6273(95)90304-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Goodliffe, J. W., Song, H., Rubakovic, A., Chang, W., Medalla, M., Weaver, C. M., et al. (2018). Differential changes to D1 and D2 medium spiny neurons in the 12-month-old Q175+/- mouse model of Huntington's disease. PLoS ONE 13:e0200626. doi: 10.1371/journal.pone.0200626

PubMed Abstract | CrossRef Full Text | Google Scholar

Grafen, A., and Hails, R. (2002). Modern Statistics for the Life Sciences. Oxford; New York, NY: Oxford University Press.

Google Scholar

Hansel, D., and Mato, G. (2013). Short-term plasticity explains irregular persistent activity in working memory tasks. J. Neurosci. 33, 133–149. doi: 10.1523/JNEUROSCI.3455-12.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Hempel, C. M., Hartman, K. H., Wang, X. J., Turrigiano, G. G., and Nelson, S. B. (2000). Multiple forms of short-term plasticity at excitatory synapses in rat medial prefrontal cortex. J. Neurophysiol. 83, 3031–3041. doi: 10.1152/jn.2000.83.5.3031

PubMed Abstract | CrossRef Full Text | Google Scholar

Herndon, J. G., Moss, M. B., Rosene, D. L., and Killiany, R. J. (1997). Patterns of cognitive decline in aged rhesus monkeys. Behav. Brain Res. 87, 25–34. doi: 10.1016/S0166-4328(96)02256-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Hestrin, S., Sah, P., and Nicoll, R. A. (1990). Mechanisms generating the time course of dual component excitatory synaptic currents recorded in hippocampal slices. Neuron 5, 247–253. doi: 10.1016/0896-6273(90)90162-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Hodgkin, A. L., and Huxley, A. F. (1990). A quantitative description of membrane current and its application to conduction and excitation in nerve. Bull. Math. Biol. 52, 25–71. doi: 10.1016/S0092-8240(05)80004-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Hof, P. R., Duan, H., Page, T. L., Einstein, M., Wicinski, B., He, Y., et al. (2002). Age-related changes in GluR2 and NMDAR1 glutamate receptor subunit protein immunoreactivity in corticocortically projecting neurons in macaque and patas monkeys. Brain Res. 928, 175–186. doi: 10.1016/S0006-8993(01)03345-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Hof, P. R., and Morrison, J. H. (2004). The aging brain: morphomolecular senescence of cortical circuits. Trends Neurosci. 27, 607–613. doi: 10.1016/j.tins.2004.07.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Ibañez, S., Luebke, J. I., Chang, W., Draguljić, D., and Weaver, C. M. (2019). Network models predict that pyramidal neuron hyperexcitability and synapse loss in the dlPFC lead to age-related spatial working memory impairment in rhesus monkeys. BioRxiv 745901. doi: 10.1101/745901

CrossRef Full Text | Google Scholar

Itskov, V., Hansel, D., and Tsodyks, M. (2011). Short-term facilitation may stabilize parametric working memory trace. Front. Comput. Neurosci. 5:40. doi: 10.3389/fncom.2011.00040

PubMed Abstract | CrossRef Full Text | Google Scholar

Izhikevich, E. M. (2003). Simple model of spiking neurons. IEEE Trans. Neural Netw. 14, 1569–1572. doi: 10.1109/TNN.2003.820440

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, M. E., Moore, L. M., and Ylvisaker, D. (1990). Minimax and maximin distance designs. J. Stat. Plan Inference. 26, 131–148. doi: 10.1016/0378-3758(90)90122-B

CrossRef Full Text | Google Scholar

King, B. R., van Ruitenbeek, P., Leunissen, I., Cuypers, K., Heise, K.-F., Santos Monteiro, T., et al. (2018). Age-related declines in motor performance are associated with decreased segregation of large-scale resting state brain networks. Cereb. Cortex 28, 4390–4402. doi: 10.1093/cercor/bhx297

PubMed Abstract | CrossRef Full Text | Google Scholar

Koen, J. D., Hauck, N., and Rugg, M. D. (2019). The relationship between age, neural differentiation, and memory performance. J. Neurosci. 39, 149–162. doi: 10.1523/JNEUROSCI.1498-18.2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Konar, A., Singh, P., and Thakur, M. K. (2016). Age-associated cognitive decline: insights into molecular switches and recovery avenues. Aging Dis. 7, 121–129. doi: 10.14336/AD.2015.1004

PubMed Abstract | CrossRef Full Text | Google Scholar

Kubicki, M., Baxi, M., Pasternak, O., Tang, Y., Karmacharya, S., Chunga, N., et al. (2019). Lifespan trajectories of white matter changes in rhesus monkeys. Cereb. Cortex 29, 1584–1593. doi: 10.1093/cercor/bhy056

PubMed Abstract | CrossRef Full Text | Google Scholar

Lacreuse, A., Kim, C. B., Rosene, D. L., Killiany, R. J., Moss, M. B., Moore, T. L., et al. (2005). Sex, age, and training modulate spatial memory in the rhesus monkey (Macaca mulatta). Behav. Neurosci. 119, 118–126. doi: 10.1037/0735-7044.119.1.118

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, A., Tan, M., and Qiu, A. (2016). Distinct aging effects on functional networks in good and poor cognitive performers. Front. Aging Neurosci. 8:215. doi: 10.3389/fnagi.2016.00215

PubMed Abstract | CrossRef Full Text | Google Scholar

Luebke, J. I., and Amatrudo, J. M. (2012). Age-related increase of sI(AHP) in prefrontal pyramidal cells of monkeys: relationship to cognition. Neurobiol. Aging 33, 1085–1095. doi: 10.1016/j.neurobiolaging.2010.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Luebke, J. I., Barbas, H., and Peters, A. (2010). Effects of normal aging on prefrontal area 46 in the rhesus monkey. Brain Res. Rev. 62:212. doi: 10.1016/j.brainresrev.2009.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Luebke, J. I., Chang, Y.-M., Moore, T. L., and Rosene, D. L. (2004). Normal aging results in decreased synaptic excitation and increased synaptic inhibition of layer 2/3 pyramidal cells in the monkey prefrontal cortex. Neuroscience 125, 277–288. doi: 10.1016/j.neuroscience.2004.01.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Luebke, J. I., Medalla, M., Amatrudo, J. M., Weaver, C. M., Crimins, J. L., Hunt, B., et al. (2015). Age-related changes to layer 3 pyramidal cells in the rhesus monkey visual cortex. Cereb. Cortex 25, 1454–1468. doi: 10.1093/cercor/bht336

PubMed Abstract | CrossRef Full Text | Google Scholar

Lundqvist, M., Herman, P., and Lansner, A. (2011). Theta and gamma power increases and alpha/beta power decreases with memory load in an attractor network model. J. Cogn. Neurosci. 23, 3008–3020. doi: 10.1162/jocn_a_00029

PubMed Abstract | CrossRef Full Text | Google Scholar

Lundqvist, M., Herman, P., and Miller, E. K. (2018). Working memory: delay activity, yes! persistent activity? Maybe not. J. Neurosci. 38, 7013–7019. doi: 10.1523/JNEUROSCI.2485-17.2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Maass, W., Natschläger, T., and Markram, H. (2002). Real-time computing without stable states: a new framework for neural computation based on perturbations. Neural Comput. 14, 2531–2560. doi: 10.1162/089976602760407955

PubMed Abstract | CrossRef Full Text | Google Scholar

Mangiafico, S. (2016). Summary and Analysis of Extension Program Evaluation in R, version 1.15.0. Available online at: http://rcompanion.org/handbook/ (accessed January 5, 2019).

Google Scholar

Marvel, C. L., Morgan, O. P., and Kronemer, S. I. (2019). How the motor system integrates with working memory. Neurosci. Biobehav. Rev. 102, 184–194. doi: 10.1016/j.neubiorev.2019.04.017

PubMed Abstract | CrossRef Full Text | Google Scholar

McDougal, R. A., Morse, T. M., Carnevale, T., Marenco, L., Wang, R., Migliore, M., et al. (2017). Twenty years of ModelDB and beyond: building essential modeling tools for the future of neuroscience. J. Comput. Neurosci. 42:1. doi: 10.1007/s10827-016-0623-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Medalla, M., and Barbas, H. (2010). Anterior cingulate synapses in prefrontal areas 10 and 46 suggest differential influence in cognitive control. J. Neurosci. 30, 16068–16081. doi: 10.1523/JNEUROSCI.1773-10.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Medalla, M., Gilman, J. P., Wang, J.-Y., and Luebke, J. I. (2017). Strength and diversity of inhibitory signaling differentiates primate anterior cingulate from lateral prefrontal cortex. J. Neurosci. 37, 4717–4734. doi: 10.1523/JNEUROSCI.3757-16.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

Medalla, M., and Luebke, J. I. (2015). Diversity of glutamatergic synaptic strength in lateral prefrontal versus primary visual cortices in the rhesus monkey. J. Neurosci. 35, 112–127. doi: 10.1523/JNEUROSCI.3426-14.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Mi, Y., Katkov, M., and Tsodyks, M. (2017). Synaptic correlates of working memory capacity. Neuron 93, 323–330. doi: 10.1016/j.neuron.2016.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, E. K., Lundqvist, M., and Bastos, A. M. (2018). Working memory 2.0. Neuron 100, 463–475. doi: 10.1016/j.neuron.2018.09.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Mongillo, G., Barak, O., and Tsodyks, M. (2008). Synaptic theory of working memory. Science 319, 1543–1546. doi: 10.1126/science.1150769

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, T. L., Bowley, B., Shultz, P., Calderazzo, S., Shobin, E., Killiany, R. J., et al. (2017). Chronic curcumin treatment improves spatial working memory but not recognition memory in middle-aged rhesus monkeys. Geroscience 39, 571–584. doi: 10.1007/s11357-017-9998-2

CrossRef Full Text | Google Scholar

Moore, T. L., Killiany, R. J., Herndon, J. G., Rosene, D. L., and Moss, M. B. (2003). Impairment in abstraction and set shifting in aged rhesus monkeys. Neurobiol. Aging 24, 125–134. doi: 10.1016/S0197-4580(02)00054-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, T. L., Killiany, R. J., Herndon, J. G., Rosene, D. L., and Moss, M. B. (2006). Executive system dysfunction occurs as early as middle-age in the rhesus monkey. Neurobiol. Aging 27, 1484–1493. doi: 10.1016/j.neurobiolaging.2005.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, T. L., Schettler, S. P., Killiany, R. J., Herndon, J. G., Luebke, J. I., Moss, M. B., et al. (2005). Cognitive impairment in aged rhesus monkeys associated with monoamine receptors in the prefrontal cortex. Behav. Brain Res. 160, 208–221. doi: 10.1016/j.bbr.2004.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Morrison, J. H., and Baxter, M. G. (2012). The aging cortical synapse: hallmarks and implications for cognitive decline. Nat. Rev. Neurosci. 13, 240–250. doi: 10.1038/nrn3200

CrossRef Full Text | Google Scholar

Motley, S. E., Grossman, Y. S., Janssen, W. G. M., Baxter, M. G., Rapp, P. R., Dumitriu, D., et al. (2018). Selective loss of thin spines in area 7a of the primate intraparietal sulcus predicts age-related working memory impairment. J. Neurosci. 38, 10467–10478. doi: 10.1523/JNEUROSCI.1234-18.2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakagawa, S., and Schielzeth, H. (2013). A general and simple method for obtaining R² from generalized linear mixed-effects models. Methods Ecol. Evol. 4, 133–142. doi: 10.1111/j.2041-210x.2012.00261.x

CrossRef Full Text | Google Scholar

Norris, C. M., Halpain, S., and Foster, T. C. (1998a). Reversal of age-related alterations in synaptic plasticity by blockade of L-type Ca2+ channels. J. Neurosci. 18, 3171–3179. doi: 10.1523/JNEUROSCI.18-09-03171.1998

PubMed Abstract | CrossRef Full Text | Google Scholar

Norris, C. M., Halpain, S., and Foster, T. C. (1998b). Alterations in the balance of protein kinase/phosphatase activities parallel reduced synaptic strength during aging. J. Neurophysiol. 80, 1567–1570. doi: 10.1152/jn.1998.80.3.1567

PubMed Abstract | CrossRef Full Text | Google Scholar

Oemisch, M., Westendorff, S., Everling, S., and Womelsdorf, T. (2015). Interareal spike-train correlations of anterior cingulate and dorsal prefrontal cortex during attention shifts. J. Neurosci. 35, 13076–13089. doi: 10.1523/JNEUROSCI.1262-15.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Oliveira, A. M. M., and Bading, H. (2011). Calcium signaling in cognition and aging-dependent cognitive decline. Biofactors 37, 168–174. doi: 10.1002/biof.148

PubMed Abstract | CrossRef Full Text | Google Scholar

Peters, A. (2007). Golgi, cajal, and the fine structure of the nervous system. Brain Res. Rev. 55, 256–263. doi: 10.1016/j.brainresrev.2006.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Peters, A. (2009). The effects of normal aging on myelinated nerve fibers in monkey central nervous system. Front. Neuroanat. 3:11. doi: 10.3389/neuro.05.011.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Peters, A., and Kemper, T. (2012). A review of the structural alterations in the cerebral hemispheres of the aging rhesus monkey. Neurobiol. Aging 33, 2357–2372. doi: 10.1016/j.neurobiolaging.2011.11.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Peters, A., Sethares, C., and Luebke, J. I. (2008). Synapses are lost during normal aging in the primate prefrontal cortex. Neuroscience 152, 970–981. doi: 10.1016/j.neuroscience.2007.07.014

CrossRef Full Text | Google Scholar

Povysheva, N. V., Zaitsev, A. V., Gonzalez-Burgos, G., and Lewis, D. A. (2013). Electrophysiological heterogeneity of fast-spiking interneurons: chandelier versus basket cells. PLoS ONE 8:e70553. doi: 10.1371/journal.pone.0070553

PubMed Abstract | CrossRef Full Text | Google Scholar

Proskovec, A. L., Heinrichs-Graham, E., and Wilson, T. W. (2016). Aging modulates the oscillatory dynamics underlying successful working memory encoding and maintenance. Hum. Brain Mapp. 37, 2348–2361. doi: 10.1002/hbm.23178

PubMed Abstract | CrossRef Full Text | Google Scholar

Rhodes, M. G. (2004). Age-related differences in performance on the Wisconsin card sorting test: a meta-analytic review. Psychol. Aging 19, 482–494. doi: 10.1037/0882-7974.19.3.482

PubMed Abstract | CrossRef Full Text | Google Scholar

Rigotti, M., Ben Dayan Rubin, D., Wang, X.-J., and Fusi, S. (2010). Internal representation of task rules by recurrent dynamics: the importance of the diversity of neural responses. Front. Comput. Neurosci. 4:24. doi: 10.3389/fncom.2010.00024

PubMed Abstract | CrossRef Full Text | Google Scholar

Rolls, E. T., and Deco, G. (2015). Stochastic cortical neurodynamics underlying the memory and cognitive changes in aging. Neurobiol. Learn Mem. 118, 150–161. doi: 10.1016/j.nlm.2014.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Rolls, E. T., Dempere-Marco, L., and Deco, G. (2013). Holding multiple items in short term memory: a neural mechanism. PLoS ONE 8:e61078. doi: 10.1371/journal.pone.0061078

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosenzweig, E. S., and Barnes, C. A. (2003). Impact of aging on hippocampal function: plasticity, network dynamics, and cognition. Prog. Neurobiol. 69, 143–179. doi: 10.1016/S0301-0082(02)00126-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Roxin, A., and Compte, A. (2016). Oscillations in the bistable regime of neuronal networks. Phys. Rev. E 94:012410. doi: 10.1103/PhysRevE.94.012410

PubMed Abstract | CrossRef Full Text | Google Scholar

RStudio Team (2015). RStudio: Integrated Development for R. RStudio, Inc., Boston, MA. Available online at: http://www.rstudio.com/

Google Scholar

Rumbell, T. H., Draguljić, D., Yadav, A., Hof, P. R., Luebke, J. I., and Weaver, C. M. (2016). Automated evolutionary optimization of ion channel conductances and kinetics in models of young and aged rhesus monkey pyramidal neurons. J. Comput. Neurosci. 41, 65–90. doi: 10.1007/s10827-016-0605-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Salin, P. A., and Prince, D. A. (1996). Spontaneous GABAA receptor-mediated inhibitory currents in adult rat somatosensory cortex. J. Neurophysiol. 75, 1573–1588. doi: 10.1152/jn.1996.75.4.1573

PubMed Abstract | CrossRef Full Text | Google Scholar

Salinas, E. (2004). Fast remapping of sensory stimuli onto motor actions on the basis of contextual modulation. J. Neurosci. 24, 1113–1118. doi: 10.1523/JNEUROSCI.4569-03.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Salthouse, T. A., Atkinson, T. M., and Berish, D. E. (2003). Executive functioning as a potential mediator of age-related cognitive decline in normal adults. J. Exp. Psychol. Gen. 132, 566–594. doi: 10.1037/0096-3445.132.4.566

PubMed Abstract | CrossRef Full Text | Google Scholar

Sandberg, A., Tegnér, J., and Lansner, A. (2003). A working memory model based on fast Hebbian learning. Network 14, 789–802. doi: 10.1088/0954-898X_14_4_309

PubMed Abstract | CrossRef Full Text | Google Scholar

Shobin, E., Bowley, M. P., Estrada, L. I., Heyworth, N. C., Orczykowski, M. E., Eldridge, S. A., et al. (2017). Microglia activation and phagocytosis: relationship with aging and cognitive impairment in the rhesus monkey. Geroscience 39, 199–220. doi: 10.1007/s11357-017-9965-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Sorel, O., and Pennequin, V. (2008). Aging of the planning process: the role of executive functioning. Brain Cogn. 66, 196–201. doi: 10.1016/j.bandc.2007.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Spruston, N., Jonas, P., and Sakmann, B. (1995). Dendritic glutamate receptor channels in rat hippocampal CA3 and CA1 pyramidal neurons. J. Physiol. 482, 325–352. doi: 10.1113/jphysiol.1995.sp020521

PubMed Abstract | CrossRef Full Text | Google Scholar

Stein, R. B. (1967). Some models of neuronal variability. Biophys. J. 7, 37–68. doi: 10.1016/S0006-3495(67)86574-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Teeter, C., Iyer, R., Menon, V., Gouwens, N., Feng, D., Berg, J., et al. (2018). Generalized leaky integrate-and-fire models classify multiple neuron types. Nat. Commun. 9:709. doi: 10.1038/s41467-017-02717-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Thibault, O., Anderson, K. L., DeMoll, C., Brewer, L. D., Landfield, P. W., and Porter, N. M. (2013). Hippocampal calcium dysregulation at the nexus of diabetes and brain aging. Eur. J. Pharmacol. 719, 34–43. doi: 10.1016/j.ejphar.2013.07.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Thibault, O., Gant, J. C., and Landfield, P. W. (2007). Expansion of the calcium hypothesis of brain aging and Alzheimer's disease: minding the store. Aging Cell 6, 307–317. doi: 10.1111/j.1474-9726.2007.00295.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Toescu, E. C., and Vreugdenhil, M. (2010). Calcium and normal brain ageing. Cell Calcium 47, 158–164. doi: 10.1016/j.ceca.2009.11.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsunada, J., Cohen, Y., and Gold, J. I. (2019). Post-decision processing in primate prefrontal cortex influences subsequent choices on an auditory decision-making task. Elife 8:e46770. doi: 10.7554/eLife.46770

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M., Datta, D., Enrwright, J., Galvin, V., Yang, S. T., Paspalas, C., et al. (2019). A novel dopamine D1 receptor agonist excites delay-dependent working memory-related neuronal firing in primate dorsolateral prefrontal cortex. Neuropharmacology 150, 46–58. doi: 10.1016/j.neuropharm.2019.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M., Gamo, N. J., Yang, Y., Jin, L. E., Wang, X.-J., Laubach, M., et al. (2011). Neuronal basis of age-related working memory decline. Nature 476, 210–213. doi: 10.1038/nature10243

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M., Yang, Y., Wang, C.-J., Gamo, N. J., Jin, L. E., Mazer, J. A., et al. (2013). NMDA receptors subserve persistent neuronal firing during working memory in dorsolateral prefrontal cortex. Neuron 77, 736–749. doi: 10.1016/j.neuron.2012.12.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Markram, H., Goodman, P. H., Berger, T. K., Ma, J., and Goldman-Rakic, P. S. (2006). Heterogeneity in the pyramidal network of the medial prefrontal cortex. Nat. Neurosci. 9, 534–542. doi: 10.1038/nn1670

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, Z., Wang, X.-J., and Wang, D.-H. (2012). From distributed resources to limited slots in multiple-item working memory: a spiking network model with normalization. J. Neurosci. 32, 11228–11240. doi: 10.1523/JNEUROSCI.0735-12.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Wimmer, K., Nykamp, D. Q., Constantinidis, C., and Compte, A. (2014). Bump attractor dynamics in prefrontal cortex explains behavioral precision in spatial working memory. Nat. Neurosci. 17, 431–439. doi: 10.1038/nn.3645

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, S., Wong, K. Y. M., Fung, C. C. A., Mi, Y., and Zhang, W. (2016). Continuous attractor neural networks: candidate of a canonical model for neural information representation [version 1; peer review: 2 approved]. F1000Research 5:156. doi: 10.12688/f1000research.7387.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiang, Z., Huguenard, J. R., and Prince, D. A. (1998). GABA_A receptor-mediated currents in interneurons and pyramidal cells of rat visual cortex. J. Physiol. 506, 715–730. doi: 10.1111/j.1469-7793.1998.715bv.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Young, M. E., Ohm, D. T., Dimitriu, D., Rapp, P. R., and Morrison, J. H. (2014). Differential effects of aging on dendritic spines in visual cortex and prefrontal cortex of the rhesus monkey. Neuroscience 274, 33–43. doi: 10.1016/j.neuroscience.2014.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Zaitsev, A. V., Povysheva, N. V., Gonzalez-Burgos, G., and Lewis, D. A. (2012). Electrophysiological classes of layer 2/3 pyramidal cells in monkey prefrontal cortex. J. Neurophysiol. 108, 595–609. doi: 10.1152/jn.00859.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Zaitsev, A. V., Povysheva, N. V., Gonzalez-Burgos, G., Rotaru, D., Fish, K. N., Krimer, L. S., et al. (2009). Interneuron diversity in layers 2-3 of monkey prefrontal cortex. Cereb. Cortex 19, 1597–1615. doi: 10.1093/cercor/bhn198

PubMed Abstract | CrossRef Full Text | Google Scholar