Memory-relevant nap sleep physiology in healthy and pathological aging

STUDY OBJECTIVES

Aging is associated with detrimental changes in sleep physiology, a process accelerated in Alzheimer's disease. Fine-tuned temporal interactions of non-rapid eye movement slow oscillations and spindles were shown to be particularly important for memory consolidation, and to deteriorate in healthy older adults. Whether this oscillatory interaction further decline in early stages of Alzheimer's disease such as mild cognitive impairment has not been investigated to date, but may have important therapeutic implications.

METHODS

Here, we assessed differences in sleep architecture and memory-relevant slow oscillation, sleep spindles and their functional coupling during a 90-min nap between healthy young and older adults, and in older patients with mild cognitive impairment. Furthermore, associations of nap-sleep characteristics with sleep-dependent memory performance change were evaluated.

RESULTS

We found significant differences between young and older healthy adults, and between young adults and patients with mild cognitive impairment, but not between healthy older adults and patients for several sleep metrics, including slow oscillation-spindle coupling. Moreover, sleep-dependent retention of verbal memories was significantly higher in young healthy adults versus older adults with and without mild cognitive impairment, but no difference between the two older groups was observed. Associations with sleep metrics were only found for pre-nap memory performances.

CONCLUSIONS

In conclusion, our results indicate changes in nap sleep physiology and sleep-related memory consolidation in older adults with and without mild cognitive impairment. Thus, interventions targeted at improving sleep physiology may help to reduce memory decline in both groups, but our study does not indicate additional benefits for patients with mild cognitive impairment.

CLINICAL TRAIL REGISTRATION

Effects of Brain Stimulation During Daytime Nap on Memory Consolidation in Younger, Healthy Subjects: https://clinicaltrials.gov/ct2/show/NCT01840865; NCT01840865. Effects of Brain Stimulation During a Daytime Nap on Memory Consolidation in Older Adults; https://clinicaltrials.gov/ct2/show/study/NCT01840839?term=01840839&draw=2&rank=1; NCT01840839. Effects of Brain Stimulation During a Daytime Nap on Memory Consolidation in Patients With Mild Cognitive Impairment; https://clinicaltrials.gov/ct2/show/NCT01782365?term=01782365&draw=2&rank=1; NCT01782365.

Resting-State Neural Firing Rate Is Linked to Cardiac-Cycle Duration in the Human Cingulate and Parahippocampal Cortices

Stimulation and functional imaging studies have revealed the existence of a large network of cortical regions involved in the regulation of heart rate. However, very little is known about the link between cortical neural firing and cardiac-cycle duration (CCD). Here, we analyze single-unit and multiunit data obtained in humans at rest, and show that firing rate covaries with CCD in 16.7% of the sample (25 of 150). The link between firing rate and CCD was most prevalent in the anterior medial temporal lobe (entorhinal and perirhinal cortices, anterior hippocampus, and amygdala), where 36% (18 of 50) of the units show the effect, and to a lesser extent in the mid-to-anterior cingulate cortex (11.1%, 5 of 45). The variance in firing rate explained by CCD ranged from 0.5 to 11%. Several lines of analysis indicate that neural firing influences CCD, rather than the other way around, and that neural firing affects CCD through vagally mediated mechanisms in most cases. These results show that part of the spontaneous fluctuations in firing rate can be attributed to the cortical control of the cardiac cycle. The fine tuning of the regulation of CCD represents a novel physiological factor accounting for spontaneous variance in firing rate. It remains to be determined whether the "noise" introduced in firing rate by the regulation of CCD is detrimental or beneficial to the cognitive information processing carried out in the parahippocampal and cingulate regions.SIGNIFICANCE STATEMENT Fluctuations in heart rate are known to be under the control of cortical structures, but spontaneous fluctuations in cortical firing rate, or "noise," have seldom been related to heart rate. Here, we analyze unit activity in humans at rest and show that spontaneous fluctuations in neural firing in the medial temporal lobe, as well as in the mid-to-anterior cingulate cortex, influence heart rate. This phenomenon was particularly pronounced in the entorhinal and perirhinal cortices, where it could be observed in one of three neurons. Our results show that part of spontaneous firing rate variability in regions best known for their cognitive role in spatial navigation and memory corresponds to precise physiological regulations.

Weak electric fields promote resonance in neuronal spiking activity: Analytical results from two-compartment cell and network models

Transcranial brain stimulation and evidence of ephaptic coupling have sparked strong interests in understanding the effects of weak electric fields on the dynamics of neuronal populations. While their influence on the subthreshold membrane voltage can be biophysically well explained using spatially extended neuron models, mechanistic analyses of neuronal spiking and network activity have remained a methodological challenge. More generally, this challenge applies to phenomena for which single-compartment (point) neuron models are oversimplified. Here we employ a pyramidal neuron model that comprises two compartments, allowing to distinguish basal-somatic from apical dendritic inputs and accounting for an extracellular field in a biophysically minimalistic way. Using an analytical approach we fit its parameters to reproduce the response properties of a canonical, spatial model neuron and dissect the stochastic spiking dynamics of single cells and large networks. We show that oscillatory weak fields effectively mimic anti-correlated inputs at the soma and dendrite and strongly modulate neuronal spiking activity in a rather narrow frequency band. This effect carries over to coupled populations of pyramidal cells and inhibitory interneurons, boosting network-induced resonance in the beta and gamma frequency bands. Our work contributes a useful theoretical framework for mechanistic analyses of population dynamics going beyond point neuron models, and provides insights on modulation effects of extracellular fields due to the morphology of pyramidal cells.

The elongated spatial structure of pyramidal neurons, which possess large apical dendrites, plays an important role for the integration of synaptic inputs and mediates sensitivity to weak extracellular electric fields. Modeling studies at the population level greatly contribute to our mechanistic understanding but face a methodological challenge because morphologically detailed neuron models are too complex for use in noisy, in-vivo like conditions and large networks in particular. Here we present an analytical approach based on a two-compartment spiking neuron model that can distinguish synaptic inputs at the apical dendrite from those at the somatic region and accounts for an extracellular field in a biophysically minimalistic way. We devised efficient methods to approximate the responses of a spatially more detailed pyramidal neuron model, and to study the spiking dynamics of single neurons and sparsely coupled large networks in the presence of fluctuating inputs. Using these methods we focused on how responses are affected by oscillatory weak fields. Our results suggest that ephaptic coupling may play a mechanistic role for oscillations of population activity and indicate the potential to entrain networks by weak electric stimulation.

Low-dimensional spike rate models derived from networks of adaptive integrate-and-fire neurons: Comparison and implementation

The spiking activity of single neurons can be well described by a nonlinear integrate-and-fire model that includes somatic adaptation. When exposed to fluctuating inputs sparsely coupled populations of these model neurons exhibit stochastic collective dynamics that can be effectively characterized using the Fokker-Planck equation. This approach, however, leads to a model with an infinite-dimensional state space and non-standard boundary conditions. Here we derive from that description four simple models for the spike rate dynamics in terms of low-dimensional ordinary differential equations using two different reduction techniques: one uses the spectral decomposition of the Fokker-Planck operator, the other is based on a cascade of two linear filters and a nonlinearity, which are determined from the Fokker-Planck equation and semi-analytically approximated. We evaluate the reduced models for a wide range of biologically plausible input statistics and find that both approximation approaches lead to spike rate models that accurately reproduce the spiking behavior of the underlying adaptive integrate-and-fire population. Particularly the cascade-based models are overall most accurate and robust, especially in the sensitive region of rapidly changing input. For the mean-driven regime, when input fluctuations are not too strong and fast, however, the best performing model is based on the spectral decomposition. The low-dimensional models also well reproduce stable oscillatory spike rate dynamics that are generated either by recurrent synaptic excitation and neuronal adaptation or through delayed inhibitory synaptic feedback. The computational demands of the reduced models are very low but the implementation complexity differs between the different model variants. Therefore we have made available implementations that allow to numerically integrate the low-dimensional spike rate models as well as the Fokker-Planck partial differential equation in efficient ways for arbitrary model parametrizations as open source software. The derived spike rate descriptions retain a direct link to the properties of single neurons, allow for convenient mathematical analyses of network states, and are well suited for application in neural mass/mean-field based brain network models.

Characterizing the dynamics of biophysically modeled, large neuronal networks usually involves extensive numerical simulations. As an alternative to this expensive procedure we propose efficient models that describe the network activity in terms of a few ordinary differential equations. These systems are simple to solve and allow for convenient investigations of asynchronous, oscillatory or chaotic network states because linear stability analyses and powerful related methods are readily applicable. We build upon two research lines on which substantial efforts have been exerted in the last two decades: (i) the development of single neuron models of reduced complexity that can accurately reproduce a large repertoire of observed neuronal behavior, and (ii) different approaches to approximate the Fokker-Planck equation that represents the collective dynamics of large neuronal networks. We combine these advances and extend recent approximation methods of the latter kind to obtain spike rate models that surprisingly well reproduce the macroscopic dynamics of the underlying neuronal network. At the same time the microscopic properties are retained through the single neuron model parameters. To enable a fast adoption we have released an efficient Python implementation as open source software under a free license.

The Influence of Mexican Hat Recurrent Connectivity on Noise Correlations and Stimulus Encoding

Noise correlations are a common feature of neural responses and have been observed in many cortical areas across different species. These correlations can influence information processing by enhancing or diminishing the quality of the neural code, but the origin of these correlations is still a matter of controversy. In this computational study we explore the hypothesis that noise correlations are the result of local recurrent excitatory and inhibitory connections. We simulated two-dimensional networks of adaptive spiking neurons with local connection patterns following Gaussian kernels. Noise correlations decay with distance between neurons but are only observed if the range of excitatory connections is smaller than the range of inhibitory connections (“Mexican hat” connectivity) and if the connection strengths are sufficiently strong. These correlations arise from a moving blob-like structure of evoked activity, which is absent if inhibitory interactions have a smaller range (“inverse Mexican hat” connectivity). Spatially structured external inputs fixate these blobs to certain locations and thus effectively reduce noise correlations. We further investigated the influence of these network configurations on stimulus encoding. On the one hand, the observed correlations diminish information about a stimulus encoded by a network. On the other hand, correlated activity allows for more precise encoding of stimulus information if the decoder has only access to a limited amount of neurons.

Extending Integrate-and-Fire Model Neurons to Account for the Effects of Weak Electric Fields and Input Filtering Mediated by the Dendrite

Transcranial brain stimulation and evidence of ephaptic coupling have recently sparked strong interests in understanding the effects of weak electric fields on the dynamics of brain networks and of coupled populations of neurons. The collective dynamics of large neuronal populations can be efficiently studied using single-compartment (point) model neurons of the integrate-and-fire (IF) type as their elements. These models, however, lack the dendritic morphology required to biophysically describe the effect of an extracellular electric field on the neuronal membrane voltage. Here, we extend the IF point neuron models to accurately reflect morphology dependent electric field effects extracted from a canonical spatial “ball-and-stick” (BS) neuron model. Even in the absence of an extracellular field, neuronal morphology by itself strongly affects the cellular response properties. We, therefore, derive additional components for leaky and nonlinear IF neuron models to reproduce the subthreshold voltage and spiking dynamics of the BS model exposed to both fluctuating somatic and dendritic inputs and an extracellular electric field. We show that an oscillatory electric field causes spike rate resonance, or equivalently, pronounced spike to field coherence. Its resonance frequency depends on the location of the synaptic background inputs. For somatic inputs the resonance appears in the beta and gamma frequency range, whereas for distal dendritic inputs it is shifted to even higher frequencies. Irrespective of an external electric field, the presence of a dendritic cable attenuates the subthreshold response at the soma to slowly-varying somatic inputs while implementing a low-pass filter for distal dendritic inputs. Our point neuron model extension is straightforward to implement and is computationally much more efficient compared to the original BS model. It is well suited for studying the dynamics of large populations of neurons with heterogeneous dendritic morphology with (and without) the influence of weak external electric fields.

How extracellular electric fields—as generated endogenously or through transcranial brain stimulation—affect the dynamics of neuronal populations is of great interest but not well understood. To study neuronal activity at the network level single-compartment neuron models have been proven very successful, because of their computational efficiency and analytical tractability. Unfortunately, these models lack the dendritic morphology to biophysically account for the effects of electric fields, and for changes in synaptic integration due to morphology alone. Here, we consider a canonical, spatially extended model neuron and characterize its responses to fluctuating synaptic input as well as an oscillatory, weak electric field. In order to accurately reproduce these responses we analytically derive an extension for the popular integrate-and-fire point neuron models. We show that the dendritic cable acts as a filter for the synaptic input current, which depends on the input location, and that an electric field modulates the neuronal spike rate strongest at a certain (preferred) field frequency. These phenomena can be successfully reproduced using integrate-and-fire models, extended by a small number of components that are straightforward to implement. The extended point models are thus well suited for studying populations of coupled neurons with different morphology, exposed to extracellular electric fields.

Controlling statistical moments of stochastic dynamical networks

We consider a general class of stochastic networks and ask which network nodes need to be controlled, and how, to stabilize and switch between desired metastable (target) states in terms of the first and second statistical moments of the system. We first show that it is sufficient to directly interfere with a subset of nodes which can be identified using information about the graph of the network only. Then we develop a suitable method for feedback control which acts on that subset of nodes and preserves the covariance structure of the desired target state. Finally, we demonstrate our theoretical results using a stochastic Hopfield network and a global brain model. Our results are applicable to a variety of (model) networks and further our understanding of the relationship between network structure and collective dynamics for the benefit of effective control.

How adaptation currents change threshold, gain, and variability of neuronal spiking

Many types of neurons exhibit spike rate adaptation, mediated by intrinsic slow K(+) currents, which effectively inhibit neuronal responses. How these adaptation currents change the relationship between in vivo like fluctuating synaptic input, spike rate output, and the spike train statistics, however, is not well understood. In this computational study we show that an adaptation current that primarily depends on the subthreshold membrane voltage changes the neuronal input-output relationship (I-O curve) subtractively, thereby increasing the response threshold, and decreases its slope (response gain) for low spike rates. A spike-dependent adaptation current alters the I-O curve divisively, thus reducing the response gain. Both types of an adaptation current naturally increase the mean interspike interval (ISI), but they can affect ISI variability in opposite ways. A subthreshold current always causes an increase of variability while a spike-triggered current decreases high variability caused by fluctuation-dominated inputs and increases low variability when the average input is large. The effects on I-O curves match those caused by synaptic inhibition in networks with asynchronous irregular activity, for which we find subtractive and divisive changes caused by external and recurrent inhibition, respectively. Synaptic inhibition, however, always increases the ISI variability. We analytically derive expressions for the I-O curve and ISI variability, which demonstrate the robustness of our results. Furthermore, we show how the biophysical parameters of slow K(+) conductances contribute to the two different types of an adaptation current and find that Ca(2+)-activated K(+) currents are effectively captured by a simple spike-dependent description, while muscarine-sensitive or Na(+)-activated K(+) currents show a dominant subthreshold component.

Adaptation controls synchrony and cluster states of coupled threshold-model neurons

We analyze zero-lag and cluster synchrony of delay-coupled nonsmooth dynamical systems by extending the master stability approach, and apply this to networks of adaptive threshold-model neurons. For a homogeneous population of excitatory and inhibitory neurons we find (i) that subthreshold adaptation stabilizes or destabilizes synchrony depending on whether the recurrent synaptic excitatory or inhibitory couplings dominate, and (ii) that synchrony is always unstable for networks with balanced recurrent synaptic inputs. If couplings are not too strong, synchronization properties are similar for very different coupling topologies, i.e., random connections or spatial networks with localized connectivity. We generalize our approach for two subpopulations of neurons with nonidentical local dynamics, including bursting, for which activity-based adaptation controls the stability of cluster states, independent of a specific coupling topology.

How adaptation shapes spike rate oscillations in recurrent neuronal networks

Neural mass signals from in-vivo recordings often show oscillations with frequencies ranging from <1 to 100 Hz. Fast rhythmic activity in the beta and gamma range can be generated by network-based mechanisms such as recurrent synaptic excitation-inhibition loops. Slower oscillations might instead depend on neuronal adaptation currents whose timescales range from tens of milliseconds to seconds. Here we investigate how the dynamics of such adaptation currents contribute to spike rate oscillations and resonance properties in recurrent networks of excitatory and inhibitory neurons. Based on a network of sparsely coupled spiking model neurons with two types of adaptation current and conductance-based synapses with heterogeneous strengths and delays we use a mean-field approach to analyze oscillatory network activity. For constant external input, we find that spike-triggered adaptation currents provide a mechanism to generate slow oscillations over a wide range of adaptation timescales as long as recurrent synaptic excitation is sufficiently strong. Faster rhythms occur when recurrent inhibition is slower than excitation and oscillation frequency increases with the strength of inhibition. Adaptation facilitates such network-based oscillations for fast synaptic inhibition and leads to decreased frequencies. For oscillatory external input, adaptation currents amplify a narrow band of frequencies and cause phase advances for low frequencies in addition to phase delays at higher frequencies. Our results therefore identify the different key roles of neuronal adaptation dynamics for rhythmogenesis and selective signal propagation in recurrent networks.

Can the human lumbar posterior columns be stimulated by transcutaneous spinal cord stimulation? A modeling study

Stimulation of different spinal cord segments in humans is a widely developed clinical practice for modification of pain, altered sensation, and movement. The human lumbar cord has become a target for modification of motor control by epidural and, more recently, by transcutaneous spinal cord stimulation. Posterior columns of the lumbar spinal cord represent a vertical system of axons and when activated can add other inputs to the motor control of the spinal cord than stimulated posterior roots. We used a detailed three-dimensional volume conductor model of the torso and the McIntyre-Richard-Grill axon model to calculate the thresholds of axons within the posterior columns in response to transcutaneous lumbar spinal cord stimulation. Superficially located large-diameter posterior column fibers with multiple collaterals have a threshold of 45.4 V, three times higher than posterior root fibers (14.1 V). With the stimulation strength needed to activate posterior column axons, posterior root fibers of large and small diameters as well as anterior root fibers are coactivated. The reported results inform on these threshold differences, when stimulation is applied to the posterior structures of the lumbar cord at intensities above the threshold of large-diameter posterior root fibers.