Using theoretical models to analyse neural development.

The development of the nervous system is an extremely complex and dynamic process. Through the continuous interplay of genetic information and changing intra- and extracellular environments, the nervous system constructs itself from precursor cells that divide and form neurons, which migrate, differentiate and establish synaptic connections. Our understanding of neural development can be greatly assisted by mathematical and computational modelling, because it allows us to bridge the gap between system-level dynamics and the lower level cellular and molecular processes. This Review shows the potential of theoretical models to examine many aspects of neural development.

Development of dendritic tonic GABAergic inhibition regulates excitability and plasticity in CA1 pyramidal neurons.

Synaptic plasticity rules change during development: while hippocampal synapses can be potentiated by a single action potential pairing protocol in young neurons, mature neurons require burst firing to induce synaptic potentiation. An essential component for spike timing-dependent plasticity is the backpropagating action potential (BAP). BAP along the dendrites can be modulated by morphology and ion channel composition, both of which change during late postnatal development. However, it is unclear whether these dendritic changes can explain the developmental changes in synaptic plasticity induction rules. Here, we show that tonic GABAergic inhibition regulates dendritic action potential backpropagation in adolescent, but not preadolescent, CA1 pyramidal neurons. These developmental changes in tonic inhibition also altered the induction threshold for spike timing-dependent plasticity in adolescent neurons. This GABAergic regulatory effect on backpropagation is restricted to distal regions of apical dendrites (>200 µm) and mediated by α5-containing GABA(A) receptors. Direct dendritic recordings demonstrate α5-mediated tonic GABA(A) currents in adolescent neurons which can modulate BAPs. These developmental modulations in dendritic excitability could not be explained by concurrent changes in dendritic morphology. To explain our data, model simulations propose a distally increasing or localized distal expression of dendritic α5 tonic inhibition in mature neurons. Overall, our results demonstrate that dendritic integration and plasticity in more mature dendrites are significantly altered by tonic α5 inhibition in a dendritic region-specific and developmentally regulated manner.

A simple rule for dendritic spine and axonal bouton formation can account for cortical reorganization after focal retinal lesions.

Lasting alterations in sensory input trigger massive structural and functional adaptations in cortical networks. The principles governing these experience-dependent changes are, however, poorly understood. Here, we examine whether a simple rule based on the neurons' need for homeostasis in electrical activity may serve as driving force for cortical reorganization. According to this rule, a neuron creates new spines and boutons when its level of electrical activity is below a homeostatic set-point and decreases the number of spines and boutons when its activity exceeds this set-point. In addition, neurons need a minimum level of activity to form spines and boutons. Spine and bouton formation depends solely on the neuron's own activity level, and synapses are formed by merging spines and boutons independently of activity. Using a novel computational model, we show that this simple growth rule produces neuron and network changes as observed in the visual cortex after focal retinal lesions. In the model, as in the cortex, the turnover of dendritic spines was increased strongest in the center of the lesion projection zone, while axonal boutons displayed a marked overshoot followed by pruning. Moreover, the decrease in external input was compensated for by the formation of new horizontal connections, which caused a retinotopic remapping. Homeostatic regulation may provide a unifying framework for understanding cortical reorganization, including network repair in degenerative diseases or following focal stroke.

Spine calcium transients induced by synaptically-evoked action potentials can predict synapse location and establish synaptic democracy.

CA1 pyramidal neurons receive hundreds of synaptic inputs at different distances from the soma. Distance-dependent synaptic scaling enables distal and proximal synapses to influence the somatic membrane equally, a phenomenon called "synaptic democracy". How this is established is unclear. The backpropagating action potential (BAP) is hypothesised to provide distance-dependent information to synapses, allowing synaptic strengths to scale accordingly. Experimental measurements show that a BAP evoked by current injection at the soma causes calcium currents in the apical shaft whose amplitudes decay with distance from the soma. However, in vivo action potentials are not induced by somatic current injection but by synaptic inputs along the dendrites, which creates a different excitable state of the dendrites. Due to technical limitations, it is not possible to study experimentally whether distance information can also be provided by synaptically-evoked BAPs. Therefore we adapted a realistic morphological and electrophysiological model to measure BAP-induced voltage and calcium signals in spines after Schaffer collateral synapse stimulation. We show that peak calcium concentration is highly correlated with soma-synapse distance under a number of physiologically-realistic suprathreshold stimulation regimes and for a range of dendritic morphologies. Peak calcium levels also predicted the attenuation of the EPSP across the dendritic tree. Furthermore, we show that peak calcium can be used to set up a synaptic democracy in a homeostatic manner, whereby synapses regulate their synaptic strength on the basis of the difference between peak calcium and a uniform target value. We conclude that information derived from synaptically-generated BAPs can indicate synapse location and can subsequently be utilised to implement a synaptic democracy.

External drive to inhibitory cells induces alternating episodes of high- and low-amplitude oscillations.

Electrical oscillations in neuronal network activity are ubiquitous in the brain and have been associated with cognition and behavior. Intriguingly, the amplitude of ongoing oscillations, such as measured in EEG recordings, fluctuates irregularly, with episodes of high amplitude alternating with episodes of low amplitude. Despite the widespread occurrence of amplitude fluctuations in many frequency bands and brain regions, the mechanisms by which they are generated are poorly understood. Here, we show that irregular transitions between sub-second episodes of high- and low-amplitude oscillations in the alpha/beta frequency band occur in a generic neuronal network model consisting of interconnected inhibitory and excitatory cells that are externally driven by sustained cholinergic input and trains of action potentials that activate excitatory synapses. In the model, we identify the action potential drive onto inhibitory cells, which represents input from other brain areas and is shown to desynchronize network activity, to be crucial for the emergence of amplitude fluctuations. We show that the duration distributions of high-amplitude episodes in the model match those observed in rat prefrontal cortex for oscillations induced by the cholinergic agonist carbachol. Furthermore, the mean duration of high-amplitude episodes varies in a bell-shaped manner with carbachol concentration, just as in mouse hippocampus. Our results suggest that amplitude fluctuations are a general property of oscillatory neuronal networks that can arise through background input from areas external to the network.

Altered temporal correlations in parietal alpha and prefrontal theta oscillations in early-stage Alzheimer disease.

Encoding and retention of information in memory are associated with a sustained increase in the amplitude of neuronal oscillations for up to several seconds. We reasoned that coordination of oscillatory activity over time might be important for memory and, therefore, that the amplitude modulation of oscillations may be abnormal in Alzheimer disease (AD). To test this hypothesis, we measured magnetoencephalography (MEG) during eyes-closed rest in 19 patients diagnosed with early-stage AD and 16 age-matched control subjects and characterized the autocorrelation structure of ongoing oscillations using detrended fluctuation analysis and an analysis of the life- and waiting-time statistics of oscillation bursts. We found that Alzheimer's patients had a strongly reduced incidence of alpha-band oscillation bursts with long life- or waiting-times (< 1 s) over temporo-parietal regions and markedly weaker autocorrelations on long time scales (1-25 seconds). Interestingly, the life- and waiting-times of theta oscillations over medial prefrontal regions were greatly increased. Whereas both temporo-parietal alpha and medial prefrontal theta oscillations are associated with retrieval and retention of information, metabolic and structural deficits in early-stage AD are observed primarily in temporo-parietal areas, suggesting that the enhanced oscillations in medial prefrontal cortex reflect a compensatory mechanism. Together, our results suggest that amplitude modulation of neuronal oscillations is important for cognition and that indices of amplitude dynamics of oscillations may prove useful as neuroimaging biomarkers of early-stage AD.

Impact of dendritic size and dendritic topology on burst firing in pyramidal cells.

Neurons display a wide range of intrinsic firing patterns. A particularly relevant pattern for neuronal signaling and synaptic plasticity is burst firing, the generation of clusters of action potentials with short interspike intervals. Besides ion-channel composition, dendritic morphology appears to be an important factor modulating firing pattern. However, the underlying mechanisms are poorly understood, and the impact of morphology on burst firing remains insufficiently known. Dendritic morphology is not fixed but can undergo significant changes in many pathological conditions. Using computational models of neocortical pyramidal cells, we here show that not only the total length of the apical dendrite but also the topological structure of its branching pattern markedly influences inter- and intraburst spike intervals and even determines whether or not a cell exhibits burst firing. We found that there is only a range of dendritic sizes that supports burst firing, and that this range is modulated by dendritic topology. Either reducing or enlarging the dendritic tree, or merely modifying its topological structure without changing total dendritic length, can transform a cell's firing pattern from bursting to tonic firing. Interestingly, the results are largely independent of whether the cells are stimulated by current injection at the soma or by synapses distributed over the dendritic tree. By means of a novel measure called mean electrotonic path length, we show that the influence of dendritic morphology on burst firing is attributable to the effect both dendritic size and dendritic topology have, not on somatic input conductance, but on the average spatial extent of the dendritic tree and the spatiotemporal dynamics of the dendritic membrane potential. Our results suggest that alterations in size or topology of pyramidal cell morphology, such as observed in Alzheimer's disease, mental retardation, epilepsy, and chronic stress, could change neuronal burst firing and thus ultimately affect information processing and cognition.

Envisioning the reward.

The primary visual cortex (area V1) is for vision. At least, that is what most researchers believe. However, in a recent issue of Science, Shuler and Bear demonstrate a correlate of reward timing in area V1. This surprising result indicates that brain circuits for reward processing are more extensive than expected and that area V1 has more functionality than previously thought.

A simple rule for axon outgrowth and synaptic competition generates realistic connection lengths and filling fractions.

Neural connectivity at the cellular and mesoscopic level appears very specific and is presumed to arise from highly specific developmental mechanisms. However, there are general shared features of connectivity in systems as different as the networks formed by individual neurons in Caenorhabditis elegans or in rat visual cortex and the mesoscopic circuitry of cortical areas in the mouse, macaque, and human brain. In all these systems, connection length distributions have very similar shapes, with an initial large peak and a long flat tail representing the admixture of long-distance connections to mostly short-distance connections. Furthermore, not all potentially possible synapses are formed, and only a fraction of axons (called filling fraction) establish synapses with spatially neighboring neurons. We explored what aspects of these connectivity patterns can be explained simply by random axonal outgrowth. We found that random axonal growth away from the soma can already reproduce the known distance distribution of connections. We also observed that experimentally observed filling fractions can be generated by competition for available space at the target neurons--a model markedly different from previous explanations. These findings may serve as a baseline model for the development of connectivity that can be further refined by more specific mechanisms.

Activity-dependent structural plasticity.

Plasticity in the brain reaches far beyond a mere changing of synaptic strengths. Recent time-lapse imaging in the living brain reveals ongoing structural plasticity by forming or breaking of synapses, motile spines, and re-routing of axonal branches in the developing and adult brain. Some forms of structural plasticity do not follow Hebbian- or anti-Hebbian paradigms of plasticity but rather appear to contribute to the homeostasis of network activity. Four decades of lesion studies have brought up a wealth of data on the mutual interdependence of neuronal activity, neurotransmitter release and neuronal morphogenesis and network formation. Here, we review these former studies on structural plasticity in the context of recent experimental studies. We compare spontaneous and experience-dependent structural plasticity with lesion-induced (reactive) structural plasticity that occurs during development and in the adult brain. Understanding the principles of neural network reorganization on a structural level is relevant for a deeper understanding of long-term memory formation as well as for the treatment of neurological diseases such as stroke.

Low-frequency stimulation induces stable transitions in stereotypical activity in cortical networks.

Reverberating spontaneous synchronized brain activity is believed to play an important role in neural information processing. Whether and how external stimuli can influence this spontaneous activity is poorly understood. Because periodic synchronized network activity is also prominent in in vitro neuronal cultures, we used cortical cultures grown on multielectrode arrays to examine how spontaneous activity is affected by external stimuli. Spontaneous network activity before and after low-frequency electrical stimulation was quantified in several ways. Our results show that the initially stable pattern of stereotypical spontaneous activity was transformed into another activity pattern that remained stable for at least 1 h. The transformations consisted of changes in single site and culture-wide network activity as well as in the spatiotemporal dynamics of network bursting. We show for the first time that low-frequency electrical stimulation can induce long-lasting alterations in spontaneous activity of cortical neuronal networks. We discuss whether the observed transformations in network activity could represent a switch in attractor state.

Genetic contributions to long-range temporal correlations in ongoing oscillations.

The amplitude fluctuations of ongoing oscillations in the electroencephalographic (EEG) signal of the human brain show autocorrelations that decay slowly and remain significant at time scales up to tens of seconds. We call these long-range temporal correlations (LRTC). Abnormal LRTC have been observed in several brain pathologies, but it has remained unknown whether genetic factors influence the temporal correlation structure of ongoing oscillations. We recorded the ongoing EEG during eyes-closed rest in 390 monozygotic and dizygotic twins and investigated the temporal structure of ongoing oscillations in the alpha- and beta-frequency bands using detrended fluctuation analysis (DFA). The strength of LRTC was more highly correlated in monozygotic than in dizygotic twins. Statistical analysis attributed up to approximately 60% of the variance in DFA to genetic factors, indicating a high heritability for the temporal structure of amplitude fluctuations in EEG oscillations. Importantly, the DFA and EEG power were uncorrelated. LRTC in ongoing oscillations are robust, heritable, and independent of power, suggesting that LRTC and oscillation power are governed by distinct biophysical mechanisms and serve different functions in the brain. We propose that the DFA method is an important complement to classical spectral analysis in fundamental and clinical research on ongoing oscillations.

Inverse relationship between adult hippocampal cell proliferation and synaptic rewiring in the dentate gyrus.

Adult neurogenesis is a key feature of the hippocampal dentate gyrus (DG). Neurogenesis is accompanied by synaptogenesis as new cells become integrated into the circuitry of the hippocampus. However, little is known to what extent the embedding of new neurons rewires the pre-existing network. Here we investigate synaptic rewiring in the DG of gerbils (Meriones unguiculatus) under different rates of adult cell proliferation caused by different rearing conditions as well as juvenile methamphetamine treatment. Surprisingly, we found that an increased cell proliferation reduced the amount of synaptic rewiring. To help explain this unexpected finding, we developed a novel model of dentate network formation incorporating neurogenesis and activity-dependent synapse formation and remodelling. In the model, we show that homeostasis of neuronal activity can account for the inverse relationship between cell proliferation and synaptic rewiring.

Avalanche dynamics of human brain oscillations: relation to critical branching processes and temporal correlations.

Human brain oscillations fluctuate erratically in amplitude during rest and exhibit power-law decay of temporal correlations. It has been suggested that this dynamics reflects self-organized activity near a critical state. In this framework, oscillation bursts may be interpreted as neuronal avalanches propagating in a network with a critical branching ratio. However, a direct comparison of the temporal structure of ongoing oscillations with that of activity propagation in a model network with critical connectivity has never been made. Here, we simulate branching processes and characterize the activity propagation in terms of avalanche life-time distributions and temporal correlations. An equivalent analysis is introduced for characterizing ongoing oscillations in the alpha-frequency band recorded with magnetoencephalography (MEG) during rest. We found that models with a branching ratio near the critical value of one exhibited power-law scaling in life-time distributions with similar scaling exponents as observed in the MEG data. The models reproduced qualitatively the power-law decay of temporal correlations in the human data; however, the correlations in the model appeared on time scales only up to the longest avalanche, whereas human data indicate persistence of correlations on time scales corresponding to several burst events. Our results support the idea that neuronal networks generating ongoing alpha oscillations during rest operate near a critical state, but also suggest that factors not included in the simple classical branching process are needed to account for the complex temporal structure of ongoing oscillations during rest on time scales longer than the duration of individual oscillation bursts.

Dynamic Hebbian Cross-Correlation Learning Resolves the Spike Timing Dependent Plasticity Conundrum.

Spike Timing-Dependent Plasticity has been found to assume many different forms. The classic STDP curve, with one potentiating and one depressing window, is only one of many possible curves that describe synaptic learning using the STDP mechanism. It has been shown experimentally that STDP curves may contain multiple LTP and LTD windows of variable width, and even inverted windows. The underlying STDP mechanism that is capable of producing such an extensive, and apparently incompatible, range of learning curves is still under investigation. In this paper, it is shown that STDP originates from a combination of two dynamic Hebbian cross-correlations of local activity at the synapse. The correlation of the presynaptic activity with the local postsynaptic activity is a robust and reliable indicator of the discrepancy between the presynaptic neuron and the postsynaptic neuron's activity. The second correlation is between the local postsynaptic activity with dendritic activity which is a good indicator of matching local synaptic and dendritic activity. We show that this simple time-independent learning rule can give rise to many forms of the STDP learning curve. The rule regulates synaptic strength without the need for spike matching or other supervisory learning mechanisms. Local differences in dendritic activity at the synapse greatly affect the cross-correlation difference which determines the relative contributions of different neural activity sources. Dendritic activity due to nearby synapses, action potentials, both forward and back-propagating, as well as inhibitory synapses will dynamically modify the local activity at the synapse, and the resulting STDP learning rule. The dynamic Hebbian learning rule ensures furthermore, that the resulting synaptic strength is dynamically stable, and that interactions between synapses do not result in local instabilities. The rule clearly demonstrates that synapses function as independent localized computational entities, each contributing to the global activity, not in a simply linear fashion, but in a manner that is appropriate to achieve local and global stability of the neuron and the entire dendritic structure.

Mathematical modelling and numerical simulation of the morphological development of neurons.

BACKGROUND: The morphological development of neurons is a very complex process involving both genetic and environmental components. Mathematical modelling and numerical simulation are valuable tools in helping us unravel particular aspects of how individual neurons grow their characteristic morphologies and eventually form appropriate networks with each other. METHODS: A variety of mathematical models that consider (1) neurite initiation (2) neurite elongation (3) axon pathfinding, and (4) neurite branching and dendritic shape formation are reviewed. The different mathematical techniques employed are also described. RESULTS: Some comparison of modelling results with experimental data is made. A critique of different modelling techniques is given, leading to a proposal for a unified modelling environment for models of neuronal development. CONCLUSION: A unified mathematical and numerical simulation framework should lead to an expansion of work on models of neuronal development, as has occurred with compartmental models of neuronal electrical activity.

Competition in neurite outgrowth and the development of nerve connections.

During the development of the nervous system, neurons form their characteristic morphologies and become assembled into synaptically connected networks. In both neuronal morphogenesis and the development of nerve connections, competition plays an important role. Although the notion of competition is commonly used in neurobiology, there is little understanding of the nature of the competitive process and the underlying molecular and cellular mechanisms. In this chapter, we review a model of competition between outgrowing neurites, as well as various models of competition that have been proposed for the refinement of connections that takes place in the development of the neuromuscular and visual systems. We describe in detail a model that links competition in the development of nerve connections with the underlying actions and biochemistry of neurotrophic factors.

Role of synaptic inhibition in spatiotemporal patterning of cortical activity.

Developmental upregulation of the GABAA receptor alpha1 subunit causes a faster decay of GABAergic inhibitory postsynaptic currents (IPSCs) in the visual cortex around the time of eye opening. In alpha1 deficient mice, a juvenile type of GABAA receptors is retained during maturation. As a result the decay time of the IPSCs is longer in alpha1-/- mice than in WT mice during the whole life span of the mice. Hence they form a valuable mouse model for studies on cellular aspects of neuronal network functioning. Using voltage sensitive dye imaging methods, we monitored the spatiotemporal excitation patterning in visual cortex slices upon local stimulation of the network. We found that in the alpha1-/- mice, the ability of the network to fire synchronously at gamma-frequencies (20-50 Hz) is diminished. This finding indicates that early onset of GABA synapse maturation is required for the normal neuronal network function in the maturating visual cortex.

H-Channels Affect Frequency, Power and Amplitude Fluctuations of Neuronal Network Oscillations.

Oscillations in network activity are ubiquitous in the brain and are involved in diverse cognitive functions. Oscillation characteristics, such as power, frequency, and temporal structure, depend on both network connectivity and intrinsic cellular properties, such as ion channel composition. An important class of channels, with key roles in regulating cell excitability, are h-channels. The h-current (Ih) is a slow, hyperpolarization-activated, depolarizing current that contributes to neuronal resonance and membrane potential. The impact of Ih on network oscillations, however, remains poorly understood. To elucidate the network effects of Ih, we used a computational model of a generic oscillatory neuronal network consisting of inhibitory and excitatory cells that were externally driven by excitatory action potentials and sustained depolarizing currents. We found that Ih increased the oscillation frequency and, in combination with external action potentials, representing input from areas outside the network, strongly decreased the synchrony of firing. As a consequence, the oscillation power and the duration of episodes during which the network exhibited high-amplitude oscillations were greatly reduced in the presence of Ih. Our results suggest that modulation of Ih or impaired expression of h-channels, as observed in epilepsy, could, by affecting oscillation dynamics, markedly alter network-level activity and potentially influence oscillation-dependent cognitive processes such as learning, memory and attention.

Homeostatic structural plasticity can account for topology changes following deafferentation and focal stroke.

After brain lesions caused by tumors or stroke, or after lasting loss of input (deafferentation), inter- and intra-regional brain networks respond with complex changes in topology. Not only areas directly affected by the lesion but also regions remote from the lesion may alter their connectivity-a phenomenon known as diaschisis. Changes in network topology after brain lesions can lead to cognitive decline and increasing functional disability. However, the principles governing changes in network topology are poorly understood. Here, we investigated whether homeostatic structural plasticity can account for changes in network topology after deafferentation and brain lesions. Homeostatic structural plasticity postulates that neurons aim to maintain a desired level of electrical activity by deleting synapses when neuronal activity is too high and by providing new synaptic contacts when activity is too low. Using our Model of Structural Plasticity, we explored how local changes in connectivity induced by a focal loss of input affected global network topology. In accordance with experimental and clinical data, we found that after partial deafferentation, the network as a whole became more random, although it maintained its small-world topology, while deafferentated neurons increased their betweenness centrality as they rewired and returned to the homeostatic range of activity. Furthermore, deafferentated neurons increased their global but decreased their local efficiency and got longer tailed degree distributions, indicating the emergence of hub neurons. Together, our results suggest that homeostatic structural plasticity may be an important driving force for lesion-induced network reorganization and that the increase in betweenness centrality of deafferentated areas may hold as a biomarker for brain repair.