ORIGINAL RESEARCH article
Volume 15 - 2021 | https://doi.org/10.3389/fncir.2021.718270
The Role of Hub Neurons in Modulating Cortical Dynamics
- 1The Edmond & Lily Safra Center for Brain Sciences, The Hebrew University of Jerusalem, Jerusalem, Israel
- 2Department of Neurobiology, The Hebrew University of Jerusalem, Jerusalem, Israel
- 3Blue Brain Project, École Polytechnique Fédérale de Lausanne (EPFL), Campus Biotech, Geneva, Switzerland
Many neurodegenerative diseases are associated with the death of specific neuron types in particular brain regions. What makes the death of specific neuron types particularly harmful for the integrity and dynamics of the respective network is not well understood. To start addressing this question we used the most up-to-date biologically realistic dense neocortical microcircuit (NMC) of the rodent, which has reconstructed a volume of 0.3 mm3 and containing 31,000 neurons, ∼37 million synapses, and 55 morphological cell types arranged in six cortical layers. Using modern network science tools, we identified hub neurons in the NMC, that are connected synaptically to a large number of their neighbors and systematically examined the impact of abolishing these cells. In general, the structural integrity of the network is robust to cells’ attack; yet, attacking hub neurons strongly impacted the small-world topology of the network, whereas similar attacks on random neurons have a negligible effect. Such hub-specific attacks are also impactful on the network dynamics, both when the network is at its spontaneous synchronous state and when it was presented with synchronized thalamo-cortical visual-like input. We found that attacking layer 5 hub neurons is most harmful to the structural and functional integrity of the NMC. The significance of our results for understanding the role of specific neuron types and cortical layers for disease manifestation is discussed.
Research at the macro- and peso-scale brain anatomy has demonstrated a clear connection between structure and function. Indeed, the global network structure of the brain was shown to be altered in diseases such as schizophrenia (Rubinov and Bullmore, 2013), bipolar disorder (Syan et al., 2018) and others (Stam, 2014). Yet, pathology takes place at the micro-scale, at the cellular and synaptic level architecture of neuronal microcircuits. How the connectomics at this level shape the dynamics and functionality of biological circuits is indeed a key question in neuroscience (Abbott et al., 2020; Amsalem et al., 2020; Turner et al., 2020). Of particular interest is the impact of structural disruption of the connectome, whether due to natural aging or due to disease. These two types of disruption are rather different. Whereas a non-selective general reduction in the number of cells was found in the aging brain, recent studies showed selective cell vulnerability associated with certain pathologies. For example, a significant decrease in the number of specific cell types in cortical areas were found in Alzheimer’s disease (Stranahan and Mattson, 2010; Fu et al., 2018; Murray et al., 2018), multiple sclerosis (Schirmer et al., 2019) and Parkinson’s disease (Hammond et al., 2007).
Experimental investigation of the role of specific cell populations in the neocortex has advanced significantly in recent years. Optogenetic methods (Deisseroth, 2015) together with genetic dissection of specific neurons (Luo et al., 2018) enables precise recording and manipulation (silencing and activating) of specific neuronal populations. Manipulating the neural activity of specific cell types during in vivo experiments is presently used to affect animal behavior (Guo et al., 2015; Carrillo-Reid et al., 2019; Robinson et al., 2020), but the effects of such cell-type-specific manipulations on circuit dynamics are rarely characterized at the network scale (Cardin et al., 2009; Pouille et al., 2009; Adesnik and Scanziani, 2010; Xue et al., 2014; Bitzenhofer et al., 2017). Consequently, we currently lack understanding of the role of particular cell populations, e.g., in the case of specific diseases, in shaping neural network dynamics and eventually network functionality.
To address this gap, we hereby utilized network theoretical approaches to explore the correlates between microcircuitry structure and function, as theoretical approaches provided new insights into the relationship between structure and function of neural networks (Newman, 2003; Avena-Koenigsberger et al., 2017; Nolte et al., 2020). Toward this end we simulated the most up-to-date, biologically realistic, dense digital reconstruction of a neocortical microcircuit, NMC (Markram et al., 2015). This 0.3 mm3 cortical circuit contains some 31,000 neurons, ∼37 million excitatory and inhibitory synapses, and 55 morphological cell types (m-types). This model circuit enables an unprecedented opportunity to directly investigate the impact of network structure on system dynamics by introducing cell-specific and layer-specific attack/damage while measuring the collective neural activity under different physiological conditions.
This cortical microcircuit displays structured architecture with multiple emerging network features at multiple levels, from the local single-cell level, through clusters and motifs of two and three neurons, to the global network-wide architectures (Gal et al., 2017). In particular, at the single-cell level, the circuit contained highly connected neurons (“hub neurons”) that either received or made significantly more connections than average. Moreover, these hub neurons belong to a surprisingly small number of cell-type subclasses and are densely interconnected among themselves forming a cell-type-specific core of hubs (“rich club”). Among all 55 m-types that constitute the mouse somatosensory cortex, the hubs belonged mainly to pyramidal cells from the deeper layers 4–6. Hub neurons with a high number of outgoing connections (“out hubs”) tend to arise from the intermediate layers (4–5), and those with a high number of incoming connections (“in hubs”) were positioned more deeply (in layers 5–6). Surprisingly, the existence of such cell-type-specific wiring specificities was found to be essentially unavoidable, emerging mostly from the asymmetrical structure of individual cortical neurons (Gal et al., 2019). However, determining the functional role of such ubiquitous wiring specificities in shaping the cortical circuit dynamics remains elusive (Setareh et al., 2017; Luccioli et al., 2018). Indeed, the question emerges: is the death (“attack”) of these hub cells, or of specific cell types, in realistic cortical microcircuits more harmful to the network dynamics as compared to that of other cell types?
Using network science tools, we first analyzed how the circuit structural connectivity was impacted after removing these cells. At the global level, the structural connectivity demonstrated small-world topology with an average synaptic distance of 2.5 synapses separating any two neurons. In the present study, we measured the increase in the average synaptic distance (mean shortest path in a network; see section “Materials and Methods”) before and after hub attacks; this measure is related to the efficiency of information flow across the network (Watts and Strogatz, 1998). Another respective measure of the network topology is the clustering coefficient, which measures by the tendency of nodes to cluster together (see section “Materials and Methods”). When a network simultaneously displays a short-averaged path length and a high clustering level it is termed as a “small-world network” (Watts and Strogatz, 1998). Reduction in the “small-worldness” of the networks might imply a reduction in efficiency of information exchange and capacity for associative memory (Bullmore and Sporns, 2009).
We next used several measures to evaluate network dynamics such as mean firing rate, coefficient of variation (CV), SPIKE-synchronization (Kreuz et al., 2015), etc. Each measure provides us with different aspects of the circuit activity and, hence, helps us understand how targeted attacks on hub neurons are more disruptive to the network functionality than attacking the same number of neurons randomly. Combining these results with structural network measures following cells’ attack sheds new light on the robustness of the NMC to a variety of attacks and, at the same time, on the functional sensitivity of the network to some of these attacks. These findings provided important insights into the impact of the death of specific cell types (e.g., due to certain diseases) on the dynamics and functionality of local cortical microcircuits.
Hub Neuron Attacks Impact the NMC Small-World Topology
To explore the structural and functional impact of attacking highly-connected hub neurons we started by ranking all neurons according to their total degree (total number of pre- and post-synaptic cells connected to a given neuron). We then removed (attacked) different quantities of these neurons, starting from the highest degree hub cells to the lowest degree.
To quantify the structural effect, we first measured the overall connectedness of the network, as captured by the size of the network’s largest connected component (giant component size; see section “Materials and Methods”). We highlight two extreme outcomes on network architecture following hubs attack (Figure 1A). On one extreme, the removal of hub neurons may completely break the network into multiple smaller unconnected components (Figure 1A top; “breakable”). On the other extreme, the removal of hub neurons will not break the network and the remaining neurons will remain in one connected giant component (Figure 1A bottom; “unbreakable”). A finer structural feature utilized here, which relates to the efficiency of network communication and computation, is the small-world topology of the network. This measure relies on two opposing requirements: a short path length between any pair of nodes/neurons (Figure 1B, right; see section “Materials and Methods”) and clustered interconnectivity, c, within groups of nodes (Figure 1B, left; see section “Materials and Methods”). Thus, the “small-worldness” of a neural network reflects the degree in which it balances the needs for global integration and local segregation of neural information (Sporns, 2013a).
Figure 1. Structural disruption of cortical microcircuits following attacks on hub neurons. (A) Schematic illustration of two network architectures having different sensitivity to a similar targeted hubs attack. (Top) The case of breakable network whereby an attack on 3 hub neurons disintegrates the network into 3 separate sub-networks. (Bottom) The case of an unbreakable network which, after an attack on 3 hub neurons, remains fully connected. (B) Schematic illustration of the two features used to define small-world networks. (Right) The shortest path length connecting two nodes (cells “a” and “b”, overplayed on the modeled neocortical microcircuit). (Left) The local clustering coefficient of a node i, ci, which is the density of connections among the neighbors of this node. In the example shown, node i, has 8 neighbors; among them only 6 connections out of all 8⋅(8–1) = 56 possible connections (ci = 6/56 ≅ 0.11). (C) The size of the giant component in the NMC as a function of the number of hubs attacked (Attack strength), black line, compared to the corresponding random attacks (orange line), demonstrating that the NMC network is “unbreakable.” (D) Mean path-length and (E) mean clustering coefficient following hub attacks (black) versus random attack (orange). These two features are particularly sensitive to hub attack. Light orange depicts 95% confidence interval in all figures.
We found that the giant component of the NMC is not broken by hub-attack as its size only negligibly decreased when a large number of hub neurons are attacked (Figure 1C, black line; see section “Materials and Methods”), showing that the global integrity of the circuit is robust to such attacks. Moreover, the effect was similar to respective random cell attacks (Figure 1C, orange line) in which a similar number of cells were selected randomly (see section “Materials and Methods”). This implies that the NMC circuit is of the unbreakable type. In contrast, the small-world topology of the circuit is more sensitive to attacks on hub neurons. For example, attacking 5,000 hub neurons increased the path length from 2.48 to 2.69 (8% increase, Figure 1D) while reducing c from 0.029 to 0.025 (12% reduction, Figure 1E). To test whether this disrupted small-world topology is expected by chance and resulting merely due to the number of eliminated nodes, we performed control random attacks with matching number of nodes (Figures 1C–E). We found that the disruption of both path-length and clustering due to hub attacks were significantly stronger than that expected from the random attacks (p < 0.001 for both, two-tailed Wilcoxon rank sum test; n1 = 10, n2 = 100; see section “Materials and Methods”). Additionally, the observed disruptions were found significant (P < 0.001, two-tailed Wilcoxon rank sum test; n1 = 10, n2 = 100; see section “Materials and Methods”) compared to that of randomly attacked networks with similar numbers of eliminated edges (Supplementary Figure 1).
Hub Neurons Are Key for Network Synchrony
The structural analysis has uncovered the disruptive effect of hub attacks on the small-world properties of the neocortical microcircuit. However, the functional implications of such structural changes are not trivial. To elucidate the functional impact of hub neurons on network dynamics, we simulated NMC networks following different cells’ attacks and examined various functional features of network activity.
It has been shown (Markram et al., 2015) that at an extracellular calcium concentration of 1.4 mM the NMC network generates spontaneously synchronous bursts at ∼1 Hz (Figure 2A and see section “Materials and Methods”). At this synchronous state, the cells’ mean firing rate is 3.7 Hz, their CV is 2.16 and SPIKE-synchronization measure is 0.25 (see section “Materials and Methods”). We simulated the circuit after removing 2,977 hub neurons (Figure 2B) or 2,977 random neurons (Figure 2D). Removal of hub neurons markedly reduced the number of bursts in the network (Figure 2B), whereas the bursting properties of the network were unaltered for the respective random attack (Figure 2D). Increasing further the number of attacked hub neurons to 7,993 completely abolished the bursting activity of the network, shifting it to the asynchronous state (Figure 2C) whereas it did not change the burstiness in the case of random cells removal (Figure 2E). Figure 2F summarizes the change in burst number due to different strengths of attacks. In addition to the reduction in burst activity due to removal of hub cells from the circuit, the coefficient of variation (std divided by mean ISI), CV, firing rate and spike-synchronization all showed a stronger reduction compared to removal of random neurons (Figures 2G–I).
Figure 2. Effect of hub attacks on the dynamics of the NMC network in the synchronous regime. (A) Raster plot (top) and time histogram (bottom) of the NMC spiking activity during spontaneous bursting state (see section “Materials and Methods”). In this state, all cortical layers tend to burst synchronously at about 1 Hz, with lower layers starting to fire earlier. (B,C) Same as (A) after attacking 2,977 and 7,993 hub neurons, respectively. (D,E) As in (B,C) but for respective random attacks. (F,G) Impact of hubs attack (black) versus random attack (orange) on network activity as a function of attack strength. (F) Impact on the number of bursts/sec. (G) Impact on average firing rate. (H) The impact on the coefficient of variation. (I) On global SPIKE synchronization measure. For all measures, hub attack is significantly more impactful (light orange depicts 95% confidence interval).
Because hub neurons are mostly excitatory (Gal et al., 2017), hub attacks primarily remove excitatory neurons from the network. Indeed, in all analyzed hub attacks, ranging up to 15,000 neurons, the percentage of excitatory neurons in the attacked neurons was above 97%, but when attacking random nodes, we converge to the full circuit distribution of E/I neurons (85% excitatory). To address this discrepancy, we performed random attacks that matched both the number and the E/I identity of nodes (Supplementary Figure 2 and see section “Materials and Methods”). Indeed, this attack was more disruptive than the completely random attacks, but still less than the hub attacks. We also show that the number of edges attacked is not the main factor for this effect (Supplementary Figure 2). This analysis demonstrates that, in the synchronous state, the impact of a neuron on the generations of collective synchrony in the NMC is more affected by their embedding in the network (“hubiness”) rather than by their physiological effect (the network E/I distribution).
Finally, we repeated the above analysis also for the asynchronous state, which is induced by setting the calcium concentration in the NMC simulations to 1.25 mM (Markram et al., 2015). At this state, the difference between attacking hubs versus random neurons is still significant, but not prominent as in the asynchronous case (Supplementary Figure 3; see section “Materials and Methods”). Indeed, the asynchronous state is, in general, more robust to cell-attacks.
Functional Implication of Layer-Specific Hub Attacks
We showed that hub neurons are more effective in driving network synchrony. These hubs potentially belong to multiple cell-types at the different layers. To further detail the impact layer-specific hubs, we measured the connectivity among excitatory and inhibitory cells within and between layers (Figure 3A). To compactly examine the connectivity among all 55 cell types in the NMC, we employed a force-directed graph drawing algorithm, whose 55 nodes depict the cell types whereas edge strengths correspond to the pairwise connection probability. In this presentation, tightly connected nodes will tend to appear closer. Inspecting the original network, the existence of large cell-type groups and clusters can be seen within each layer (Figure 3B large nodes). After attacking 15,000 hub neurons, several changes were prominent (Figure 3C). The network layout had spread more widely, indicating that the strength of the connections between cell types is reduced. Additionally, the large nodes of L5 almost disappeared, hinting to their possible impact on disruption of network functionality (see below).
Figure 3. Sensitivity of network dynamic to hub attack is mostly attributed to attacks on L5 PCs. (A) Percentages of synapses for excitatory (brown-red arrows) and inhibitory (blue-green arrows) connections in the NMC microcircuit (layer 1 omitted). Arrow width and corresponding numbers indicate the percentage of total synapses formed by this pathway (omitted for pathways with <1% of synapses). The total percentage of plotted synapses is 98%; the remaining 2% originate in layer 1. Rectangle sizes are proportional to the sizes of the corresponding number of excitatory or inhibitory populations. (B) The connectivity among layers and cell-types is visualized using the force-directed graph algorithm (see section “Materials and Methods”). The network is composed of 55 morphological cell-types (circular nodes, colors match to that in A); edges strength corresponds to pairwise connection probabilities. Strongly connected cell-types are displayed closer in space. (C) Same as (B), but for the network following an attack on 15,000 hub neurons. Note the disappearance of the large nodes in Layer 5. (D–G) Different quantification of network activity for layer-specific hub attack, global hub attack and random nodes attack. (D) Number of bursts/sec, (E) average network firing rate, (F) coefficient of variation, (G) global SPIKE synchronization measure (see section “Materials and Methods”). Note that attack of L5 pyramidal cells is the most disruptive layer attack.
To test the importance of the different layers for network activity we performed layer-specific hub attacks by simulating the circuit while removing hub neurons in specific layers. We found that, in general, attacking L5 hub neurons is the most disruptive attack as it caused the largest change in all the functional measures used (Figures 3D,F,G, but see the impact of L4 in Figure 3E). When revisiting the high-level circuit connectivity (Figure 3A) one observes that L5 excitatory cells are the most interconnected population in the circuit per cell (the overall percentage of incoming and outgoing connections); this is probably the reason for the high functional influence of L5 attack. L6 attacks also resulted in a large change of the functional measurements except for the number of bursts; similarly, L6 excitatory cells are also relatively highly interconnected (Figure 3A). We also note that when attacking L2/3 there is a slight increase in the network firing rate, this might be due to the relatively large part of inhibitory neurons in this layer.
We summarize this section by noting that attacking hub neurons is most disruptive to network dynamics when hub neurons are attacked at all layers (Figures 3D–G, dashed black line). Layer 5 is the most-sensitive layer to such an attack. As neurons in different cortical layers belong to different genetic types (Gouwens et al., 2019; Yuste et al., 2020), this result shows that, although the same number of neurons might be degraded due to different pathologies that target specific genetic cell types (e.g., L5 thick-tufted pyramidal cells), they will have a very different impact on the overall dynamics of the cortical network and, thus, on the manifestation of specific diseases.
Functional Implication of Hub Cells on Thalamic Input Processing
The above sections have demonstrated that hub attacks are significantly more effective in disrupting circuit-wide synchronization in the spontaneous synchronized case. In this section, we set to test whether this observation is general enough, and also valid for the case where the synchronized activity is generated by realistic sensory input. Toward this end we innervated the NMC circuit by 574 thalamic fibers (the thalamo-cortical, TC, input). These TC axons project mostly to neurons in lower layers 3 and 5 (Figure 4A), where some neurons might receive up to 750 thalamic synapses. Each reconstructed axon is making synapses on dendrites that are adjacent to its path (Figure 4B; see section “Materials and Methods”), functionally impacting a vertically confined space (Amsalem et al., 2020).
Figure 4. Network synchrony due to thalamic input is disrupted following attack on hub neurons. (A) Schematic illustration of the NMC circuit, each neuron is color-coded by the number of thalamic synapses it receives. (B) The spatial distribution of the thalamic input is illustrated by showing the cortical postsynaptic neurons receiving inputs from two exemplar TC axons (green and blue axons and respective colored cortical neurons). (C) Mean thalamic activity (bottom) simulating their response moving bars at 1 Hz (top, see section “Materials and Methods”). Bin size = 45 ms. (D) Raster plots of the circuit responding to the thalamic input. This circuit fires asynchronously at its spontaneous state ([Ca2+]0 = 1.25 mM; see Supplementary Figure 3A and section “Materials and Methods”). Red, inhibitory neurons; black, excitatory neurons. (E,F) The response of the circuit to the thalamic input after attack of 2,500 random nodes (E) and 2,500 hub neurons (F). (G) Circuit’s SPIKE-synchronization profile (Kreuz et al., 2015) in response to thalamic input is sharpened in the intact and randomly-attacked networks as compared to that of the TC input itself (blue line) and it decreased dramatically after hub attacks (black line), even below that of the TC input. (H–K) Different quantifications of network activity under thalamic input for hub and random attacks (as in Figures 2F,G).
For these simulations we set the calcium concentration value to 1.25 mM; this results in the network being in a spontaneous asynchronous state (as in Supplementary Figure 3A; see also Figure 15 and Supplementary Figure 12 in Markram et al., 2015). We then simulated a grating drifting at 1 Hz by generating the firing rate of the thalamic axons from an inhomogeneous Poisson process with a time-varying rate that followed a sinusoidal function (Figure 4C; see section “Materials and Methods”). The circuit responded by following the oscillatory input firing with highly time-locked synchronized bursts of spikes at 1 Hz (Figure 4D). We repeated the simulation following the attack on 2,500 random neurons (Figure 4E) or 2,500 hub neurons (Figure 4F). In the random attack the circuit continued to follow the oscillatory input and the response remained synchronized to the input. However, following attack on hub neurons, the circuit response was much less synchronized.
We next quantified the circuit activity using SPIKE-synchronization time profile in response to the thalamic input for different cases (Figure 4G). The blue line in this figure shows the spike-synchronization measure of the TC axons whereas the gray line depicts the spike-synchronization of the cortical neurons. We found that the circuit strongly sharpens the synchronicity of the thalamic input (Figure 4G, compare gray to blue line), and that attacking random nodes (in this case 2,500) only slightly reduced this sharpening (orange line). In contrast, attacking (2,500) hub neurons reduced the synchronization in response to the TC input dramatically (Figure 4G, black line); this case is even less synchronized than the thalamic input itself (Figure 4G, compare black line to blue line). We further conducted a complete set of simulations while attacking random or hubs neurons, and quantified different functional features (Figures 4H–K). Hub attacks were much more destructive compared to the random attacks and caused a larger change of all measures. Interestingly removing hub neurons reduced the SPIKE-synchronization profile to a value which is lower than that of the input, showing the strong dependence of the circuit ability to follow and sharpen synchronized input on hub neurons. These results highlight the importance of hub neurons in processing sensory input, clearly demonstrating that the integrity of the cortical hub neurons is critical for the fast and reliable response of the cortical circuit to sensory information.
We introduced in this work a network-based approach to investigate the relation between a cortical circuit structure to its function by removing cells according to different criteria in a highly detailed simulation. Additionally, we proposed a network-based approach for identifying potentially interesting neurons. Our analysis shows that the importance of a neuron in maintaining synchronous activity is more affected by their embedding in the network (e.g., the neurons’ in/out degree, its “hubness”) rather than strictly by their physiology.
To examine the importance of different cells and layers, we simulated the network activity after removing hub cells globally or from specific layers and compared the results to control models where random cells were removed. We discovered that hub neuron attacks have the largest change of structural network measures, leading to loss of the small-world properties of the simulated NMC (Figures 1C–E). Accordingly, attacks on this population resulted in the largest decrease in the network synchrony, firing rate and number of bursts (Figure 2). The attack changed the network response from spontaneously synchronous to asynchronous, resulting in no bursts and reduced CV, as noted by another recent research (Nolte et al., 2019).
Among attacks targeting specific layers, mimicking a more biologically plausible scenario, we found that attacking L5 hub neurons resulted in the largest effect on all functional measures (Figures 3D–G) and therefore is the most distributive attack. We believe that this phenomenon is rooted in the high interconnectivity of L5 excitatory neurons (Figure 3A). Interestingly the specific genetic profile of L5 excitatory neurons is widely used to optogenetically target and record from this subset of neurons (de Vries et al., 2020), and open the possibility of examining our predictions by specifically silencing this subpopulation while recording from the neocortex.
Hub neurons play an important role not only in maintaining spontaneous network oscillatory activity but also in the processing of sensory input from other brain areas. When thalamic (sensory) drifting sinusoidal input impinged on our modeled cortical circuit, the circuit not only followed the thalamic synchrony, but resulted in activity that was more synchronized. We then found that attacking hub cells caused a significant reduction of the synchrony of the cortical column with respect to the oscillation of the TC input (Figures 4G–K), eventually resulting in activity that is less synchronized than the input. Nevertheless, random attacks seem to have a minor effect on thalamic input processing by the network, demonstrating yet again the robustness of the cortical microcircuit to random cell death. These results highlight the functional role of hub neurons in fast processing of sensory information in the cortical microcircuit.
In a previous study we have demonstrated that both in-hubs and out-hubs neurons belong to a small subset of cortical neuron types (Gal et al., 2017; Figures 3A,B). It was recently shown that, although any disease-associated genes are expressed in multiple cell types, the pathologic variants affect primarily specific cell types (Guan et al., 2021), and see also related recent work of disease map and specific cell type in the retina by Siegert et al. (2012). Therefore, it is highly likely that dysfunction of specific genes might abolish specific hub cells. We further emphasize that hub neurons (both in-hubs and out-hubs) require more metabolic energy as they receive/form a larger number of synapses and this might entail earlier death (e.g., in neurodegenerative diseases) compared to non-hub neurons (see e.g., Pathak et al., 2013).
Our findings can also be seen as a demonstration of how network science theory is implemented in realistic networks. As the modeled cortical microcircuit was shown to display small-world topology with a high clustering coefficient and short mean path, we now can systematically characterize how the targeting of fundamental components of the network, the hubs, indeed leads to loss of this topology in the neural microcircuit, and results in major functional disruptions. Robustness of complex networks to random attacks and vulnerability to removal of highly connected nodes was previously shown theoretically (Watts and Strogatz, 1998; Boccaletti et al., 2006). Here, we provided concrete evidence that hubs fulfill their theoretical key role in affecting the dynamics of the network, in a highly detailed biological model circuit. The loss of small-world topology among brain regions was shown to be related to neurodegenerative diseases using fMRI data (Sporns, 2013b). Our theoretical experiments suggest that a decrease in network clustering and an increase in the mean short path length can cause functional failures also at the microscale level of resolution.
To conclude, while neurons are usually characterized according to genetic markers, morphology and physiology (Berg et al., 2020; Yuste et al., 2020), we showed how a specific structural measure such as the number of synapses that defines hub cells, has a direct effect on the network functionality. As hub cells are specific subtypes of neurons, it is possible to use in vivo cell-specific knockout experiments to explore the behavioral implications of the neural network functional disruption suggested in our work.
Materials and Methods
NMC Connectedness Measures
The general connectedness of the NMC was characterized by identifying connected components of the network. A strongly connected component is a group of nodes in which any node is reachable from any other node through a directed path (a series of nodes and directed edges). Intuitively, a strongly connected component reflects a group of recurrently interlinked neurons that could give rise to an anatomical module with functional specialization.
The first property of the small-world analysis is based on the length of the shortest path lij between pairs of nodes in the network. A path length between two nodes in the network is expressed as the number of connections along that path. To generalize this property for the entire network, the characteristic path length (l) of a network was used, which is the mean shortest path length averaged over all pairs of neurons.
Clearly, this measure is well-defined in connected networks where any node is reachable from any other node. The NMC network initially contained a single giant component of 31,329 ± 5 neurons that were mutually reachable (99.95% ± 0.01% of all neurons; see above).
The second property of the small-world analysis is captured by the tendency of nodes to cluster together. The local clustering of individual nodes measures the level at which the neighbors of a node are interconnected among themselves. Let the binary (unweighted) adjacency matrix of a directed network be denoted by A; then the local clustering coefficient ci of a node i is defined as,
where depicts the total degree (in-degree + out-degree) of node i. Essentially, in directed networks, this definition reflects the ratio of the number of triangles among a node and its neighbors to the number of all possible triangles that could have been formed (Fagiolo, 2007). The value of ci ranges from 0 (none of the neighbors are connected to each other) to 1 (all neighbors are mutually connected). The network-wide clustering coefficient (c) is computed by averaging over all local clustering coefficients.
Hubs Versus the Reference Random Attacks
To selectively attack the highly connected hub neurons we started by ranking all neurons according to their total degree (total number of pre- and post-synaptic cells). Then, we performed several attacks at different strengths (the number of removed hubs). In each attack, given required number of cells to attack (s, attack strength), we removed the top degreed cells.
To test the significance of the results observed in hub attacked networks we performed several types of random attacks for comparison. In the first, most naïve random attack, for each hub attack at strength s we performed 10 matching random attacks in which randomly selected s neurons were attacked. In the second, for each hub attack we counted the number of excitatory sE and inhibitory sI neurons that were attacked (s = sE + sI); we then performed 10 matching random attacks in which sE excitatory neurons and sI inhibitory neurons were selected randomly. We note that a few hub neurons might also be selected, by chance, in these random selection controls.
Structural Analysis Statistical Tests
To compare the structural disruption of hub attacks to that of random node attacks (Figure 1) we compared the structural metric values (mean path length and clustering coefficient) of the strongest attacks. Specifically, we took the ten strongest hub attacks (n1 = 10), and the ten matching random attacks for each strength (n2 = 100). For both metrics (path length or clustering coefficient) a two-tailed Wilcoxon rank sum test indicated that the disruption was greater for hub attacks than for matching control (n1 = 10, n2 = 100, P < 0.001).
For random edge comparisons (Supplementary Figure 1) we took the 100 strongest attacks (n2 = 100) and compared to the 10 closest hub attacks (n1 = 10). In agreement with the previous control, also here, a two-tailed Wilcoxon rank sum test indicated that the disruption was greater for hub attacks than for matching control (n1 = 10, n2 = 100, P < 0.001), for both metrics.
Dense Model of Neocortical Microcircuit (NMC)
Simulations were performed on a previously published model of the NMC of a two-year-old rat of a neocortical microcircuit. Full details on the constructing of the circuits and its simulation methods were described in Markram et al. (2015). The microcircuit (Figure 1B) consisted of 31,346 biophysical Hodgkin-Huxley 3D reconstructed NEURON models with around 7.8 million synaptic connections forming around 36.4 million synapses. Synaptic connectivity between 55 distinct morphological types of neurons (m-types) was predicted algorithmically and constrained by experimental data (Reimann et al., 2015). The densities of ion-channels on morphologically-detailed neuron models were optimized to reproduce the behavior of different electrical neuron types (e-types) and synaptic dynamics recorded in vitro (Van Geit et al., 2016). Simulations were run on HPE SGI 8600 supercomputer (BlueBrain 5) using NEURON (Carnevale and Hines, 2006) and CoreNEURON (Kumbhar et al., 2019).
Simulation of Baseline Spontaneous Activity
To account for the missing long-range connections and missing neuromodulators, neurons were depolarized with a noisy somatic current injection of 100% of first spike threshold (Markram et al., 2015, Figure 15). In addition, synapses spontaneous release probability was modified by setting the extracellular calcium concentration [Ca2+]o. Two conditions were tested, [Ca2+]o of 1.25 and 1.4 mM each positioning the circuit in different activity regimes (Markram et al., 2015, Figure 15). Synaptic conductances and kinetics are as in Markram et al. (2015). Each attack was simulated twice with different randomization of the noisy step currents and timing of the spontaneous synaptic release, each time for 10 s.
Simulation Oscillatory Thalamocortical Input to the NMC
The oscillatory thalamocortical input to the NMC (Figure 4) was generated by following the same principle as in Amsalem et al. (2016) – the spike times of the axons from inhomogeneous Poisson process with time-varying rate as,
where t is time in seconds, f the frequency of the oscillatory input was set to 1 Hz, B5.83 is a factor, so that the mean firing rate of the oscillatory part would be equal to 7 Hz, and b is the baseline spontaneous firing rate which was set to 1 Hz. For Supplementary Figure 1 each simulation was of 10 s, and conducted 5 times with different randomly generated thalamic input.
We detected bursts by extracting the multivariate SPIKE-synchronization profile (Mulansky and Kreuz, 2016) for each simulation, smoothing the result using a running mean filter of ∼200 ms and then counting the number of events larger than half the maximal synchronization (but at least larger than 0.15).
Force-Directed Graph Layout
Visualization of complex networks in an informative and meaningful way is a challenging task. How to position a large number of nodes, densely interconnected with non-obvious organization, in a two-dimensional layout that can expose inherent symmetries and structures such as hubs and clusters?
To provide a layout in which the distance between nodes (cell types) is more or less proportional to their edge weight (connection probability), we employed a Force-directed graph drawing algorithm. The algorithm is based on a physical model that assigns different forces among the nodes. On one hand, to promote attraction between connected nodes spring-like attractive forces, which depend on the distance and edge weight, are simulated. On the other hand, to avoid overlapping of nodes, repulsive forces (such as Coulomb’s law between electrically charged particles) are simulated to separate all pairs of nodes. By iteratively determining all the forces and moving the nodes accordingly, the system gets closer to an equilibrium where all forces add up to zero, and the position of the nodes stays stable.
Here, we used the implementation from D3.js library1.
Figures were created using Matplotlib (Hunter, 2007). For analysis we used Python and Numpy (Harris et al., 2020).
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.
EG, OA, and IS conceived the study and wrote the manuscript. EG and OA carried out the simulations and the analysis. AS, ML, FS, and HM participated in discussions and helped writing the manuscript. HM and FS developed the in silico microcircuit and provided the respective simulation framework and data. All authors contributed to the article and approved the submitted version.
This work was made possible through the Patrick and Lina Drahi Foundation (PLFA), the Blue Brain Project, a research center of the École polytechnique fédérale de Lausanne (EPFL), from the Swiss government’s ETH Board of the Swiss Federal Institutes of Technology, the Gatsby Charitable Foundation, and the NIH Grant Agreement U01MH114812.
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.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
We thank Noam Kahlon for the helpful discussions related to this project and for conducting the initial hub-attack experiments.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncir.2021.718270/full#supplementary-material
Supplementary Figure 1 | Hub attack disruption is stronger than random attacks with matching number of edges. Same as Figure 1, only the attack strength (x-axis) is measured by the number of removed edges.
Supplementary Figure 2 | Matched E/I ratio to hub attack and matched with number of edges (A–D) hub attack (black), random attack (orange) and random attack with E/I ratio that is matched to the hub attacks (green) effects on network activity as a function of attack strength. (A) Average network firing rate, (B) the number of bursts/sec, (C) coefficient of variation, and (D) global SPIKE-synchronization measure (see section “Materials and Methods” for details about the different attacks and measures). For all measures, hub attack is much more impactful. Panels (E–H) same as in Figure 2, only the attack strength (x-axis) is measured by the number of removed edges.
Supplementary Figure 3 | Effect of hub attacks on the dynamics of the NMC network in the asynchronous regime. (A) Raster plot (top) and PSTH (bottom) of the NMC during spontaneous asynchronous state (see section “Materials and Methods”). (B,C) Same as (A) after attacking 3,134 and 8,149 hub neurons, respectively. (D,E) As in (B,C) but for respective random attacks. (F–G) Impact of hubs (black) versus random (orange) attacks on network activity as a function of attack strength. (F) On the number of bursts/sec, (G) on average network firing rate, (H) on coefficient of variation, (I) on global SPIKE-synchronization measure (see section “Materials and Methods”). For all measures, hub attack is much more impactful.
Abbott, L. F., Bock, D. D., Callaway, E. M., Denk, W., Dulac, C., Fairhall, A. L., et al. (2020). The mind of a mouse. Cell 182, 1372–1376. doi: 10.1016/j.cell.2020.08.010
Adesnik, H., and Scanziani, M. (2010). Lateral competition for cortical space by layer-specific horizontal circuits. Nature 464, 1155–1160. doi: 10.1038/nature08935
Amsalem, O., King, J., Reimann, M., Ramaswamy, S., Muller, E., Markram, H., et al. (2020). Dense computer replica of cortical microcircuits unravels cellular underpinnings of auditory surprise response. bioRxiv [Preprint] doi: 10.1101/2020.05.31.126466
Amsalem, O., Van Geit, W., Muller, E., Markram, H., and Segev, I. (2016). From neuron biophysics to orientation selectivity in electrically coupled networks of neocortical L2/3 large basket cells. Cereb. Cortex 26, 3655–3668. doi: 10.1093/cercor/bhw166
Avena-Koenigsberger, A., Misic, B., and Sporns, O. (2017). Communication dynamics in complex brain networks. Nat. Rev. Neurosci. 19, 17–33. doi: 10.1038/nrn.2017.149
Berg, J., Sorensen, S. A., Ting, J. T., Miller, J. A., Chartrand, T., Buchin, A., et al. (2020). Human cortical expansion involves diversification and specialization of supragranular intratelencephalic-projecting neurons. bioRxiv [Preprint] doi: 10.1101/2020.03.31.018820
Bitzenhofer, S. H., Ahlbeck, J., Wolff, A., Wiegert, J. S., Gee, C. E., Oertner, T. G., et al. (2017). Layer-specific optogenetic activation of pyramidal neurons causes beta–gamma entrainment of neonatal networks. Nat. Commun. 8:14563. doi: 10.1038/ncomms14563
Boccaletti, S., Latora, V., Moreno, Y., Chavez, M., and Hwang, D. U. (2006). Complex networks: structure and dynamics. Phys. Rep. 424, 175–308. doi: 10.1016/j.physrep.2005.10.009
Bullmore, E., and Sporns, O. (2009). Complex brain networks: graph theoretical analysis of structural and functional systems. Nat. Rev. Neurosci. 10, 186–198. doi: 10.1038/nrn2575
Cardin, J. A., Carlén, M., Meletis, K., Knoblich, U., Zhang, F., Deisseroth, K., et al. (2009). Driving fast-spiking cells induces gamma rhythm and controls sensory responses. Nature 459, 663–667. doi: 10.1038/nature08002
Carnevale, N. T., and Hines, M. L. (2006). The NEURON Book. Cambridge: Cambridge University Press, doi: 10.1017/CBO9780511541612
Carrillo-Reid, L., Han, S., Yang, W., Akrouh, A., and Yuste, R. (2019). Controlling visually guided behavior by holographic recalling of cortical ensembles. Cell 178, 447–457.e5. doi: 10.1016/j.cell.2019.05.045
de Vries, S. E. J., Lecoq, J. A., Buice, M. A., Groblewski, P. A., Ocker, G. K., Oliver, M., et al. (2020). A large-scale standardized physiological survey reveals functional organization of the mouse visual cortex. Nat. Neurosci. 23, 138–151. doi: 10.1038/s41593-019-0550-9
Deisseroth, K. (2015). Optogenetics: 10 years of microbial opsins in neuroscience. Nat. Neurosci. 18, 1213–1225. doi: 10.1038/nn.4091
Fagiolo, G. (2007). Clustering in complex directed networks. Phys. Rev. E - Stat. Nonlinear, Soft Matter Phys. 76:026107. doi: 10.1103/PhysRevE.76.026107
Fu, H., Hardy, J., and Duff, K. E. (2018). Selective vulnerability in neurodegenerative diseases. Nat. Neurosci. 21, 1350–1358. doi: 10.1038/s41593-018-0221-2
Gal, E., London, M., Globerson, A., Ramaswamy, S., Reimann, M. W., Muller, E., et al. (2017). Rich cell-type-specific network topology in neocortical microcircuitry. Nat. Neurosci. 20, 1004–1013.
Gal, E., Perin, R., Markram, H., London, M., and Segev, I. (2019). Neuron geometry underlies a universal local architecture in neuronal networks. bioRxiv [Preprint] doi: 10.1101/656058
Gouwens, N. W., Sorensen, S. A., Berg, J., Lee, C., Jarsky, T., Ting, J., et al. (2019). Classification of electrophysiological and morphological neuron types in the mouse visual cortex. Nat. Neurosci. 22, 1182–1195. doi: 10.1038/s41593-019-0417-0
Guan, J., Lin, Y., Wang, Y., Gao, J., and Ji, G. (2021). An analytical method for the identification of cell type-specific disease gene modules. J. Transl. Med. 19:20. doi: 10.1186/s12967-020-02690-5
Guo, J. Z., Graves, A. R., Guo, W. W., Zheng, J., Lee, A., Rodríguez-Gonzá Lez, J., et al. (2015). Cortex commands the performance of skilled movement. Elife 4:10774. doi: 10.7554/eLife.10774
Hammond, C., Bergman, H., and Brown, P. (2007). Pathological synchronization in Parkinson’s disease: networks, models and treatments. Trends Neurosci. 30, 357–364. doi: 10.1016/j.tins.2007.05.004
Harris, C. R., Millman, K. J., van der Walt, S. J., Gommers, R., Virtanen, P., Cournapeau, D., et al. (2020). Array programming with NumPy. Nature 585, 357–362. doi: 10.1038/s41586-020-2649-2
Hunter, J. D. (2007). Matplotlib: A 2D graphics environment. Comput. Sci. Eng. 9, 90–95. doi: 10.1109/MCSE.2007.55
Kreuz, T., Bozanic, N., and Mulansky, M. (2015). SPIKE-Synchronization: a parameter-free and time-resolved coincidence detector with an intuitive multivariate extension. BMC Neurosci. 16:170. doi: 10.1186/1471-2202-16-S1-P170
Kumbhar, P., Hines, M., Fouriaux, J., Ovcharenko, A., King, J., Delalondre, F., et al. (2019). CoreNEURON: An Optimized Compute Engine for the NEURON Simulator. Front. Neuroinform. 13:63. doi: 10.3389/fninf.2019.00063
Luccioli, S., Angulo-Garcia, D., Cossart, R., Malvache, A., Módol, L., Sousa, V. H., et al. (2018). Modeling driver cells in developing neuronal networks. PLoS Comput. Biol. 14:e1006551. doi: 10.1371/journal.pcbi.1006551
Luo, L., Callaway, E. M., and Svoboda, K. (2018). Genetic dissection of neural circuits: a decade of progress. Neuron. 98, 256–281. doi: 10.1016/j.neuron.2018.03.040
Markram, H., Muller, E., Ramaswamy, S., Reimann, M. W., Abdellah, M., Sanchez, C. A., et al. (2015). Reconstruction and simulation of neocortical microcircuitry. Cell 163, 456–492. doi: 10.1016/j.cell.2015.09.029
Mulansky, M., and Kreuz, T. (2016). PySpike—A Python library for analyzing spike train synchrony. SoftwareX 5, 183–189.
Murray, H. C., Swanson, M. E. V., Dieriks, B. V., Turner, C., Faull, R. L. M., and Curtis, M. A. (2018). Neurochemical characterization of PSA-NCAM+ cells in the human brain and phenotypic quantification in alzheimer’s disease entorhinal cortex. Neuroscience 372, 289–303. doi: 10.1016/j.neuroscience.2017.12.019
Newman, M. E. J. (2003). The structure and function of complex networks. SIAM Rev. 45, 167–256.
Nolte, M., Gal, E., Markram, H., and Reimann, M. W. (2020). Impact of higher order network structure on emergent cortical activity. Netw. Neurosci. 4, 292–314. doi: 10.1162/netn_a_00124
Nolte, M., Reimann, M. W., King, J. G., Markram, H., and Muller, E. B. (2019). Cortical reliability amid noise and chaos. Nat. Commun. 10:3792. doi: 10.1038/s41467-019-11633-8
Pathak, D., Berthet, A., and Nakamura, K. (2013). Energy failure: does it contribute to neurodegeneration? Ann. Neurol. 74:506. doi: 10.1002/ANA.24014
Pouille, F., Marin-Burgin, A., Adesnik, H., Atallah, B. V., and Scanziani, M. (2009). Input normalization by global feedforward inhibition expands cortical dynamic range. Nat. Neurosci. 12, 1577–1585. doi: 10.1038/nn.2441
Reimann, M. W., King, J. G., Muller, E. B., Ramaswamy, S., and Markram, H. (2015). An algorithm to predict the connectome of neural microcircuits. Front. Comput. Neurosci. 9:120. doi: 10.3389/fncom.2015.00120
Robinson, N. T. M., Descamps, L. A. L., Russell, L. E., Buchholz, M. O., Bicknell, B. A., Antonov, G. K., et al. (2020). Targeted activation of hippocampal place cells drives memory-guided spatial behavior. Cell 183, 1586.e–1599.e. doi: 10.1016/j.cell.2020.09.061
Rubinov, M., and Bullmore, E. (2013). Schizophrenia and abnormal brain network hubs. Dialogues Clin. Neurosci. 15, 339–349. doi: 10.31887/dcns.2013.15.3/mrubinov
Schirmer, L., Velmeshev, D., Holmqvist, S., Kaufmann, M., Werneburg, S., Jung, D., et al. (2019). Neuronal vulnerability and multilineage diversity in multiple sclerosis. Nature 573, 75–82. doi: 10.1038/s41586-019-1404-z
Setareh, H., Deger, M., Petersen, C. C. H., and Gerstner, W. (2017). Cortical Dynamics in Presence of Assemblies of Densely Connected Weight-Hub Neurons. Front. Comput. Neurosci. 11:52. doi: 10.3389/fncom.2017.00052
Siegert, S., Cabuy, E., Scherf, B. G., Kohler, H., Panda, S., Le, Y.-Z., et al. (2012). Transcriptional code and disease map for adult retinal cell types. Nat. Neurosci. 2012 153, 487–495. doi: 10.1038/nn.3032
Sporns, O. (2013a). Network attributes for segregation and integration in the human brain. Curr. Opin. Neurobiol. 23, 162–171. doi: 10.1016/j.conb.2012.11.015
Sporns, O. (2013b). Structure and function of complex brain networks. Dialogues Clin. Neurosci. 15, 247–262. doi: 10.1137/S003614450342480
Stam, C. J. (2014). Modern network science of neurological disorders. Nat. Rev. Neurosci. 15, 683–695.
Stranahan, A. M., and Mattson, M. P. (2010). Selective vulnerability of neurons in layer II of the entorhinal cortex during aging and Alzheimer’s disease. Neural Plast. 2010:108190. doi: 10.1155/2010/108190
Syan, S. K., Minuzzi, L., Smith, M., Costescu, D., Allega, O. R., Hall, G. B. C., et al. (2018). Brain structure and function in women with Comorbid bipolar and premenstrual dysphoric disorder. Front. Psychiatry 8:301. doi: 10.3389/fpsyt.2017.00301
Turner, N. L., Macrina, T., Alexander Bae, J., Yang, R., Wilson, A. M., Schneider-Mizell, C., et al. (2020). Multiscale and multimodal reconstruction of cortical structure and function. bioRxiv [Preprint] doi: 10.1101/2020.10.14.338681
Van Geit, W., Gevaert, M., Chindemi, G., Rössert, C., Courcol, J.-D., Muller, E. B., et al. (2016). BluePyOpt: leveraging open source software and cloud infrastructure to optimise model parameters in neuroscience. Front. Neuroinform. 10:17. doi: 10.3389/fninf.2016.00017
Watts, D. J., and Strogatz, S. H. (1998). Collective dynamics of “small-world” networks. Nature 393, 440–442. doi: 10.1038/30918
Xue, M., Atallah, B. V., and Scanziani, M. (2014). Equalizing excitation-inhibition ratios across visual cortical neurons. Nature 511, 596–600. doi: 10.1038/nature13321
Yuste, R., Hawrylycz, M., Aalling, N., Aguilar-Valles, A., Arendt, D., Arnedillo, R. A., et al. (2020). A community-based transcriptomics classification and nomenclature of neocortical cell types. Nat. Neurosci. 23, 1456–1468. doi: 10.1038/s41593-020-0685-8
Keywords: network science, cortical microcircuitry, network synchronicity, scale free network, hub neurons, small world topology, cell ablation strategy
Citation: Gal E, Amsalem O, Schindel A, London M, Schürmann F, Markram H and Segev I (2021) The Role of Hub Neurons in Modulating Cortical Dynamics. Front. Neural Circuits 15:718270. doi: 10.3389/fncir.2021.718270
Received: 31 May 2021; Accepted: 24 August 2021;
Published: 24 September 2021.
Edited by:Yoshiyuki Kubota, National Institute for Physiological Sciences (NIPS), Japan
Reviewed by:Yasuhiro Tsubo, Ritsumeikan University, Japan
Diana Martinez, Cooper Medical School of Rowan University, United States
Copyright © 2021 Gal, Amsalem, Schindel, London, Schürmann, Markram and Segev. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Eyal Gal, firstname.lastname@example.org; Oren Amsalem, email@example.com; Idan Segev, firstname.lastname@example.org
†These authors have contributed equally to this work and share first authorship