Cortical Network Models of Firing Rates in the Resting and Active States Predict BOLD Responses

Measurements of blood oxygenation level dependent (BOLD) signals have produced some surprising observations. One is that their amplitude is proportional to the entire activity in a region of interest and not just the fluctuations in this activity. Another is that during sleep and anesthesia the average BOLD correlations between regions of interest decline as the activity declines. Mechanistic explanations of these phenomena are described here using a cortical network model consisting of modules with excitatory and inhibitory neurons, taken as regions of cortical interest, each receiving excitatory inputs from outside the network, taken as subcortical driving inputs in addition to extrinsic (intermodular) connections, such as provided by associational fibers. The model shows that the standard deviation of the firing rate is proportional to the mean frequency of the firing when the extrinsic connections are decreased, so that the mean BOLD signal is proportional to both as is observed experimentally. The model also shows that if these extrinsic connections are decreased or the frequency of firing reaching the network from the subcortical driving inputs is decreased, or both decline, there is a decrease in the mean firing rate in the modules accompanied by decreases in the mean BOLD correlations between the modules, consistent with the observed changes during NREM sleep and under anesthesia. Finally, the model explains why a transient increase in the BOLD signal in a cortical area, due to a transient subcortical input, gives rises to responses throughout the cortex as observed, with these responses mediated by the extrinsic (intermodular) connections.     

An estimation of the absolute number of axons indicates that human cortical areas are sparsely connected

The tracts between cortical areas are conceived as playing a central role in cortical information processing, but their actual numbers have never been determined in humans. Here, we estimate the absolute number of axons linking cortical areas from a whole-cortex diffusion MRI (dMRI) connectome, calibrated using the histologically measured callosal fiber density. Median connectivity is estimated as approximately 6,200 axons between cortical areas within hemisphere and approximately 1,300 axons interhemispherically, with axons connecting functionally related areas surprisingly sparse. For example, we estimate that <5% of the axons in the trunk of the arcuate and superior longitudinal fasciculi connect Wernicke’s and Broca’s areas. These results suggest that detailed information is transmitted between cortical areas either via linkage of the dense local connections or via rare, extraordinarily privileged long-range connections.

Using data from Human Connectome Project to estimate the absolute number of axons linking cortical areas yields surprisingly sparse connectivity; reconciling large-scale functional synchronization with sparse anatomical connectivity presents a challenge for our present understanding of human brain organization.     

A modular latching chain

Many cognitive tasks involve transitions between distinct mental processes, which may range from discrete states to complex strategies. The ability of cortical networks to combine discrete jumps with continuous glides along ever changing trajectories, dubbed latching dynamics, may be essential for the emergence of the unique cognitive capacities of modern humans. Novel trajectories have to be followed in the multidimensional space of cortical activity for novel behaviours to be produced; yet, not everything changes: several lines of evidence point at recurring patterns in the sequence of activation of cortical areas in a variety of behaviours. To extend a mathematical model of latching dynamics beyond the simple unstructured auto-associative Potts network previously analysed, we introduce delayed structured connectivity and hetero-associative connection weights, and we explore their effects on the dynamics. A modular model in the small-world regime is considered, with modules arranged on a ring. The synaptic weights include a standard auto-associative component, stabilizing distinct patterns of activity, and a hetero-associative component, favoring transitions from one pattern, expressed in one module, to the next, in the next module. We then study, through simulations, how structural parameters, like those regulating rewiring probability, noise and feedback connections, determine sequential association dynamics.     

THE EFFECT OF INDUCED HYPOTHERMIA UPON OXYGEN CONSUMPTION IN THE RAT BRAIN

The effect of hypothermia upon cerebral metabolic rate for oxygen (CMR_O2) was studied in artificially ventilated rats, anaesthetized with nitrous oxide. Cerebral blood flow was measured with a modification of the Kety and Schmidt technique using ¹³³xenon. CMR_O2, was found to decrease linearily with temperature in the temperature range 37°C-22°C. At normal temperatures CMR_O2, fell by about 5 per cent per degree C. At a body temperature of 22°C both cerebral blood flow and CMR_O2, were reduced to about 25 per cent of normal.     

Role of omega-3 fatty acids in brain development and function: potential implications for the pathogenesis and prevention of psychopathology

The principle omega-3 fatty acid in brain, docosahexaenoic acid (DHA), accumulates in the brain during perinatal cortical expansion and maturation. Animal studies have demonstrated that reductions in perinatal brain DHA accrual are associated with deficits in neuronal arborization, multiple indices of synaptic pathology including deficits in serotonin and mesocorticolimbic dopamine neurotransmission, neurocognitive deficits, and elevated behavioral indices of anxiety, aggression, and depression. In primates and humans, preterm delivery is associated with deficits in fetal cortical DHA accrual, and children/adolescents born preterm exhibit deficits in cortical gray matter maturation, neurocognitive deficits particularly in the realm of attention, and increased risk for attention-deficit/hyperactivity disorder (ADHD) and schizophrenia. Individuals diagnosed with ADHD or schizophrenia exhibit deficits in cortical gray matter maturation, and medications found to be efficacious in the treatment of these disorders increase cortical and striatal dopamine neurotransmission. These associations in conjunction with intervention trials showing enhanced cortical visual acuity and cognitive outcomes in preterm and term infants fed DHA, suggest that perinatal deficits in brain DHA accrual may represent a preventable neurodevelopmental risk factor for the subsequent emergence of psychopathology.     

Specificity of cortical synaptic connectivity: Emphasis on perspectives gained from quantitative electron microscopy

This report traces the historical development of concepts regarding the specificity of synaptic connectivity in the cerebral cortex as viewed primarily from the perspective of electron microscopy. The occurrence of stereotypical patterns of connection (e.g., contrasting synaptic patterns on the surfaces of spiny vs. non-spiny neurons, the general consistency with which axonal pathways impinge on and originate within specific cortical areas and layers, triadic synaptic relationships) implies that cortical connectivity is highly structured. The high degree of order characterizing many aspects of cortical organization is mirrored by an equally ordered arrangement of synaptic connections between specific types of neurons. This observation is based on quantitative electron microscopic studies of synapses between identified neurons and from the results of correlative anatomical/electrophysiological investigations. The recognition of recurring synaptic patterns and responses between specific neurons has generated increased support for the notion of specificity of synaptic connections at the expense of randomness, but the role of specificity in cortical function is an unresolved question. At the core of cortical processing lie myriad possibilities for computation provided by the wealth of synaptic connections involving each cortical neuron. Specificity, by limiting possibilities for connection, can impose an order on synaptic interactions even as processes of dynamic selection or synaptic remodeling ensure the constant formation and dissolution of cortical circuits. These operations make maximal use of the richness of cortical synaptic connections to produce a highly flexible system, irrespective of the degree of randomness or specificity that obtains for cortical wiring at any particular time.     

A computational study of the role of spike broadening in synaptic facilitation of Hermissenda

Pavlovian conditioning in Hermissenda produces a decrease in voltage-dependent (I_K,A and I_Ca) and Ca²⁺-dependent (I_K,Ca) currents, and an increase in the action potential (AP) duration in type B-photoreceptors. In addition, synaptic connections between B and A photoreceptors and B photoreceptor and type I interneurons are facilitated. The increase in AP duration, produced by decreasing one or more K⁺ currents, may account for synaptic facilitation. The present study examined this issue by using a mathematical model of the B-photoreceptor and the neurosimulator SNNAP. In the model, decreasing g_K,A by 70% increased the duration of the AP in the terminal by 41% and Ca²⁺ influx by 30%. However, if the decrease in g_K,A was combined with a decrease in g_Ca, similar to what has been reported experimentally, the Ca²⁺ influx decreased by 54%. Therefore, the concomitant change in I_Ca counter-acted the broadening-induced increase in Ca²⁺ influx in the synaptic terminal. This result suggests that a spike-duration independent process must contribute to the synaptic facilitation observed following Pavlovian conditioning.     

Heat shock partially dissociates the overlapping modules of the yeast protein-protein interaction network: a systems level model of adaptation

Network analysis became a powerful tool giving new insights to the understanding of cellular behavior. Heat shock, the archetype of stress responses, is a well-characterized and simple model of cellular dynamics. S. cerevisiae is an appropriate model organism, since both its protein-protein interaction network (interactome) and stress response at the gene expression level have been well characterized. However, the analysis of the reorganization of the yeast interactome during stress has not been investigated yet. We calculated the changes of the interaction-weights of the yeast interactome from the changes of mRNA expression levels upon heat shock. The major finding of our study is that heat shock induced a significant decrease in both the overlaps and connections of yeast interactome modules. In agreement with this the weighted diameter of the yeast interactome had a 4.9-fold increase in heat shock. Several key proteins of the heat shock response became centers of heat shock-induced local communities, as well as bridges providing a residual connection of modules after heat shock. The observed changes resemble to a 'stratus-cumulus' type transition of the interactome structure, since the unstressed yeast interactome had a globally connected organization, similar to that of stratus clouds, whereas the heat shocked interactome had a multifocal organization, similar to that of cumulus clouds. Our results showed that heat shock induces a partial disintegration of the global organization of the yeast interactome. This change may be rather general occurring in many types of stresses. Moreover, other complex systems, such as single proteins, social networks and ecosystems may also decrease their inter-modular links, thus develop more compact modules, and display a partial disintegration of their global structure in the initial phase of crisis. Thus, our work may provide a model of a general, system-level adaptation mechanism to environmental changes.     

Spike Timing-Dependent Synaptic Plasticity in Visual Cortex: A Modeling Study

Recent studies show that synaptic modification depends critically on the relative spike timing of pre- and postsynaptic neurons. Here we explore the functional implications of spike timing-dependent synaptic plasticity in the visual cortex using a model circuit with modifiable intracortical excitatory connections. First we simulated the experiments using two-point stimuli, in which two visual stimuli in a topographically represented feature space were repeatedly presented in quick succession, and found that tuning of the cortical neurons was modified in a manner similar to that observed experimentally. We then explored the dependence of results on the model parameter and identified the intracortical parameters that were critical for the magnitude of the shifts and obtained a simple relationship between the amount of shift and (S = (_EXTC^rec_exc)/_INHC^rec_inh). Finally we investigated the effects of moving stimuli in a topographically represented visual space and found that they can effectively induce spike timing-dependent modification of the intracortical connections. It suggests the importance of moving stimuli in dynamic modification of the cortical maps through spike timing-dependent synaptic plasticity.     

Spike suppression in a local cortical circuit induced by transcranial magnetic stimulation

Transcranial magnetic stimulation (TMS) noninvasively interferes with human cortical function, and is widely used as an effective technique for probing causal links between neural activity and cognitive function. However, the physiological mechanisms underlying TMS-induced effects on neural activity remain unclear. We examined the mechanism by which TMS disrupts neural activity in a local circuit in early visual cortex using a computational model consisting of conductance-based spiking neurons with excitatory and inhibitory synaptic connections. We found that single-pulse TMS suppressed spiking activity in a local circuit model, disrupting the population response. Spike suppression was observed when TMS was applied to the local circuit within a limited time window after the local circuit received sensory afferent input, as observed in experiments investigating suppression of visual perception with TMS targeting early visual cortex. Quantitative analyses revealed that the magnitude of suppression was significantly larger for synaptically-connected neurons than for isolated individual neurons, suggesting that intracortical inhibitory synaptic coupling also plays an important role in TMS-induced suppression. A conventional local circuit model of early visual cortex explained only the early period of visual suppression observed in experiments. However, models either involving strong recurrent excitatory synaptic connections or sustained excitatory input were able to reproduce the late period of visual suppression. These results suggest that TMS targeting early visual cortex disrupts functionally distinct neural signals, possibly corresponding to feedforward and recurrent information processing, by imposing inhibitory effects through intracortical inhibitory synaptic connections.     

Characterization of synaptic dysfunction in an in vitro corticostriatal model system of Huntington’s disease

Huntington’s disease (HD) is an autosomal-dominant inherited neurodegenerative disease resulting from expanded amino acid (CAG) repeat in the gene that encodes protein huntingtin (Htt). HD remains incurable for now. A lot of evidence implicates aberrant synaptic connection between cortical and striatal neurons, a key component of HD pathophysilogy, which also leads to cognitive decline and motor disorders. In the present work synaptic activity between cortical and striatal neurons was studied on the corticostriatal co-culture model system of HD. Culture was prepared from HD mouse model YAC128. It was shown that first impairment appears on day 14 in vitro. Interestingly, these alterations occur in cortical neurons. Their activity in YAC128 cultures was higher than in cultures of wild-type neurons. At the same time, there were no differences in morphology of spines in striatal neurons. However, using novel optogenetic approach, we demonstrated that synaptic connections are already dysfunctional in YAC128 cultures. On day 19 in vitro the activity of cortical neurons in YAC128 cultures was reduced, which led to alterations on the post-synaptic side. Dendric spines of medium spiny neurons transformed and disappeared, which is possibly the main reason of neurodegenerative mechanisms during the HD development.     

Extremely Dilute Modular Neuronal Networks: Neocortical Memory Retrieval Dynamics

A model of columnar networks of neocortical association areas is studied. The neuronal network is composed of many Hebbian autoassociators, or modules, each of which interacts with a relatively small number of the others, randomly chosen. Any module encodes and stores a number of elementary percepts, or features. Memory items, or patterns, are peculiar combinations of features sparsely distributed over the multi-modular network. Any feature stored in any module can be involved in several of the stored patterns; feature-sharing is in fact source of local ambiguities and, consequently, a potential cause of erroneous memory retrieval spreading through the model network in pattern completion tasks.The memory retrieval dynamics of the large modular autoassociator is investigated by combining mathematical analysis and numerical simulations. An oscillatory retrieval process is proposed that is very efficient in overcoming feature-sharing drawbacks; it requires a mechanism that modulates the robustness of local attractors to noise, and neuronal activity sparseness such that quiescent and active modules are about equally noisy to any post-synaptic module.Moreover, it is shown that statistical correlation between 'kinds' of features across the set of memory patterns can be exploited to obtain a more efficient achievement of memory retrieval capabilities.It is also shown that some spots of the network cannot be reached by retrieval activity spread if they are not directly cued by the stimulus. The locations of these activity isles depend on the pattern to retrieve, while their extension only depends (in large networks) on statistics of inter-modular connections and stored patterns. The existence of activity isles determines an upper-bound to retrieval quality that does not depend on the specific retrieval dynamics adopted, nor on whether feature-sharing is permitted. The oscillatory retrieval process nearly saturates this bound.     

Ultrastructure of the synapses of the initial axonal segment of the pyramidal neuron

Ultrastructure of the proximal part of the axon in the neurons, identified according to a number of morphological signs as pyramidal, has been studied in the layer III of the cat cerebral hemisphere sensomotor cortex. In sections, tangential to the cortical surface, in the initial axonal segment, a submembranous osmophilic layer and fasciculi of microtubules are revealed. On the initial segment spines are found, they contain cysterns resembling by their structure the spine system of the dendritic spines. Axonal terminals revealed along the axonal distribution are in contact both with the axonal trunk and with the spines. Regarding the initial segment, they are presynaptic, contain oval synaptic vesicles and form symmetric axo-axonal synapses only. In transversal sections axonal terminals are detected, arranging on the surface of the initial segment mostly as single ones, in longitudinal sections they are seen as clusters. Analysing the author's data and those from the literature, a conclusion is made that in intact animals the synaptic contacts at the initial segment of the axon are the only form of axo-axonal synapses in the neocortex.     

Generalized Orientation Responses in the Defective Stage of Schizophrenia

In order to understand the pathogenesis of defective conditions in schizophrenia, and consequently of schizophrenia in general, it is necessary to study various aspects of higher nervous activity with the employment of various physiological techniques. In this connection, the study of orientation reactions, which are a necessary condition for completion of normal cortical reflex connections, offers good prospects (1, 2).     

Autophosphorylation of alphaCaMKII is required for ocular dominance plasticity

Experience is a powerful sculptor of developing neural connections. In the primary visual cortex (V1), cortical connections are particularly susceptible to the effects of sensory manipulation during a postnatal critical period. At the molecular level, this activity-dependent plasticity requires the transformation of synaptic depolarization into changes in synaptic weight. The molecule alpha calcium-calmodulin kinase type II (alphaCaMKII) is known to play a central role in this transformation. Importantly, alphaCaMKII function is modulated by autophosphorylation, which promotes Ca(2+)-independent kinase activity. Here we show that mice possessing a mutant form of alphaCaMKII that is unable to autophosphorylate show impairments in ocular dominance plasticity. These results confirm the importance of alphaCaMKII in visual cortical plasticity and suggest that synaptic changes induced by monocular deprivation are stored specifically in glutamatergic synapses made onto excitatory neurons.     

Modeling of Substance P and 5-HT Induced Synaptic Plasticity in the Lamprey Spinal CPG: Consequences for Network Pattern Generation

Consequences of synaptic plasticity in the lamprey spinal CPG are analyzed by means of simulations. This is motivated by the effects substance P (a tachykinin) and serotonin (5-hydroxytryptamin; 5-HT) have on synaptic transmission in the locomotor network. Activity-dependent synaptic depression and potentiation have recently been shown experimentally using paired intracellular recordings. Although normally activity-dependent plasticity presumably does not contribute to the patterning of network activity, this changes in the presence of the neuromodulators substance P and 5-HT, which evoke significant plasticity. Substance P can induce a faster and larger depression of inhibitory connections but potentiation of excitatory inputs, whereas 5-HT induces facilitation of both inhibitory and excitatory inputs. Changes in the amplitude of the first postsynaptic potential are also seen. These changes could thus be a potential mechanism underlying the modulatory role these substances have on the rhythmic network activity.The aim of the present study has been to implement the activity dependent synaptic depression and facilitation induced by substance P and 5-HT into two alternative models of the lamprey spinal locomotor network, one relying on reciprocal inhibition for bursting and one in which each hemicord is capable of oscillations. The consequences of the plasticity of inhibitory and excitatory connections are then explored on the network level.In the intact spinal cord, tachykinins and 5-HT, which can be endogenously released, increase and decrease the frequency of the alternating left-right burst pattern, respectively. The frequency decreasing effect of 5-HT has previously been explained based on its conductance decreasing effect on K _Ca underlying the postspike afterhyperpolarization (AHP). The present simulations show that short-term synaptic plasticity may have strong effects on frequency regulation in the lamprey spinal CPG. In the network model relying on reciprocal inhibition, the observed effects substance P and 5-HT have on network behavior (i.e., a frequency increase and decrease respectively) can to a substantial part be explained by their effects on the total extent and time dynamics of synaptic depression and facilitation. The cellular effects of these substances will in the 5-HT case further contribute to its network effect.     

Restoration of Glutamatergic Transmission by Dopamine D4 Receptors in Stressed Animals

The prefrontal cortex (PFC), a key brain region for cognitive and emotional processes, is highly regulated by dopaminergic inputs. The dopamine D₄ receptor, which is enriched in PFC, has been implicated in mental disorders, such as attention deficit-hyperactivity disorder and schizophrenia. Recently we have found homeostatic regulation of AMPA receptor-mediated synaptic transmission in PFC pyramidal neurons by the D₄ receptor, providing a potential mechanism for D₄ in stabilizing cortical excitability. Because stress is tightly linked to adaptive and maladaptive changes associated with mental health and disorders, we examined the synaptic actions of D₄ in stressed rats. We found that neural excitability was elevated by acute stress and dampened by repeated stress. D₄ activation produced a potent reduction of excitatory transmission in acutely stressed animals and a marked increase of excitatory transmission in repeatedly stressed animals. These effects of D₄ targeted GluA2-lacking AMPA receptors and relied on the bi-directional regulation of calcium/calmodulin kinase II activity. The restoration of PFC glutamatergic transmission in stress conditions may enable D₄ receptors to serve as a synaptic stabilizer in normal and pathological conditions.     

Spatial Embedding and Wiring Cost Constrain the Functional Layout of the Cortical Network of Rodents and Primates

Mammals show a wide range of brain sizes, reflecting adaptation to diverse habitats. Comparing interareal cortical networks across brains of different sizes and mammalian orders provides robust information on evolutionarily preserved features and species-specific processing modalities. However, these networks are spatially embedded, directed, and weighted, making comparisons challenging. Using tract tracing data from macaque and mouse, we show the existence of a general organizational principle based on an exponential distance rule (EDR) and cortical geometry, enabling network comparisons within the same model framework. These comparisons reveal the existence of network invariants between mouse and macaque, exemplified in graph motif profiles and connection similarity indices, but also significant differences, such as fractionally smaller and much weaker long-distance connections in the macaque than in mouse. The latter lends credence to the prediction that long-distance cortico-cortical connections could be very weak in the much-expanded human cortex, implying an increased susceptibility to disconnection syndromes such as Alzheimer disease and schizophrenia. Finally, our data from tracer experiments involving only gray matter connections in the primary visual areas of both species show that an EDR holds at local scales as well (within 1.5 mm), supporting the hypothesis that it is a universally valid property across all scales and, possibly, across the mammalian class.     

Spectral power and cortical interactions in the beta2 rhythm in normal subjects and schizophrenic patients

Spectral and coherence characteristics of the beta 2 rhythm (20-40 Hz) were compared in patients with acute schizophrenia (the first episode), patients with chronic schizophrenia (disease duration more than two years), and healthy subjects (control group) during cognitive task performance. Examination showed the "excessive" spectral power of this rhythm in the prefrontal cortical areas in patients with acute schizophrenia and its "insufficient" spectral power in all derivations of the right hemisphere in patients with the chronic form as compared to the controls. The similarity of the beta 2 spectral power distribution in different cortical areas measured by covariance method decreased with increase in psychopathologic manifestations. Coherent analysis revealed a substantially lower number of cortico-cortical functional connections in schizophrenic patients than in healthy subjects, absence of interhemispheric connections in patients with acute schizophrenia, and a slight increase in their number if patients with the chronic disease. However, as distinct from the controls revealing many inter- and intrahemispheric connections in all cortical areas, in schizophrenics the interhemispheric connections were observed only in the posterior cortical regions. The results point to a disorder of interhemispheric interaction in schizophrenia.     

Role of GABAA-Mediated Inhibition and Functional Assortment of Synapses onto Individual Layer 4 Neurons in Regulating Plasticity Expression in Visual Cortex

Layer 4 (L4) of primary visual cortex (V1) is the main recipient of thalamocortical fibers from the dorsal lateral geniculate nucleus (LGN_d). Thus, it is considered the main entry point of visual information into the neocortex and the first anatomical opportunity for intracortical visual processing before information leaves L4 and reaches supra- and infragranular cortical layers. The strength of monosynaptic connections from individual L4 excitatory cells onto adjacent L4 cells (unitary connections) is highly malleable, demonstrating that the initial stage of intracortical synaptic transmission of thalamocortical information can be altered by previous activity. However, the inhibitory network within L4 of V1 may act as an internal gate for induction of excitatory synaptic plasticity, thus providing either high fidelity throughput to supragranular layers or transmittal of a modified signal subject to recent activity-dependent plasticity. To evaluate this possibility, we compared the induction of synaptic plasticity using classical extracellular stimulation protocols that recruit a combination of excitatory and inhibitory synapses with stimulation of a single excitatory neuron onto a L4 cell. In order to induce plasticity, we paired pre- and postsynaptic activity (with the onset of postsynaptic spiking leading the presynaptic activation by 10ms) using extracellular stimulation (ECS) in acute slices of primary visual cortex and comparing the outcomes with our previously published results in which an identical protocol was used to induce synaptic plasticity between individual pre- and postsynaptic L4 excitatory neurons. Our results indicate that pairing of ECS with spiking in a L4 neuron fails to induce plasticity in L4-L4 connections if synaptic inhibition is intact. However, application of a similar pairing protocol under GABA_ARs inhibition by bath application of 2μM bicuculline does induce robust synaptic plasticity, long term potentiation (LTP) or long term depression (LTD), similar to our results with pairing of pre- and postsynaptic activation between individual excitatory L4 neurons in which inhibitory connections are not activated. These results are consistent with the well-established observation that inhibition limits the capacity for induction of plasticity at excitatory synapses and that pre- and postsynaptic activation at a fixed time interval can result in a variable range of plasticity outcomes. However, in the current study by virtue of having two sets of experimental data, we have provided a new insight into these processes. By randomly mixing the assorting of individual L4 neurons according to the frequency distribution of the experimentally determined plasticity outcome distribution based on the calculated convergence of multiple individual L4 neurons onto a single postsynaptic L4 neuron, we were able to compare then actual ECS plasticity outcomes to those predicted by randomly mixing individual pairs of neurons. Interestingly, the observed plasticity profiles with ECS cannot account for the random assortment of plasticity behaviors of synaptic connections between individual cell pairs. These results suggest that connections impinging onto a single postsynaptic cell may be grouped according to plasticity states.