Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log out
  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log out
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Cover ArticleResearch Articles, Development/Plasticity/Repair

Mesoscale Architecture Shapes Initiation and Richness of Spontaneous Network Activity

Samora Okujeni, Steffen Kandler and Ulrich Egert
Journal of Neuroscience 5 April 2017, 37 (14) 3972-3987; DOI: https://doi.org/10.1523/JNEUROSCI.2552-16.2017
Samora Okujeni
1Bernstein Center Freiburg and
2Biomicrotechnology, IMTEK—Department of Microsystems Engineering, University of Freiburg, 79110 Freiburg, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Samora Okujeni
Steffen Kandler
1Bernstein Center Freiburg and
2Biomicrotechnology, IMTEK—Department of Microsystems Engineering, University of Freiburg, 79110 Freiburg, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ulrich Egert
1Bernstein Center Freiburg and
2Biomicrotechnology, IMTEK—Department of Microsystems Engineering, University of Freiburg, 79110 Freiburg, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ulrich Egert
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Spontaneous activity in the absence of external input, including propagating waves of activity, is a robust feature of neuronal networks in vivo and in vitro. The neurophysiological and anatomical requirements for initiation and persistence of such activity, however, are poorly understood, as is their role in the function of neuronal networks. Computational network studies indicate that clustered connectivity may foster the generation, maintenance, and richness of spontaneous activity. Since this mesoscale architecture cannot be systematically modified in intact tissue, testing these predictions is impracticable in vivo. Here, we investigate how the mesoscale structure shapes spontaneous activity in generic networks of rat cortical neurons in vitro. In these networks, neurons spontaneously arrange into local clusters with high neurite density and form fasciculating long-range axons. We modified this structure by modulation of protein kinase C, an enzyme regulating neurite growth and cell migration. Inhibition of protein kinase C reduced neuronal aggregation and fasciculation of axons, i.e., promoted uniform architecture. Conversely, activation of protein kinase C promoted aggregation of neurons into clusters, local connectivity, and bundling of long-range axons. Supporting predictions from theory, clustered networks were more spontaneously active and generated diverse activity patterns. Neurons within clusters received stronger synaptic inputs and displayed increased membrane potential fluctuations. Intensified clustering promoted the initiation of synchronous bursting events but entailed incomplete network recruitment. Moderately clustered networks appear optimal for initiation and propagation of diverse patterns of activity. Our findings support a crucial role of the mesoscale architectures in the regulation of spontaneous activity dynamics.

SIGNIFICANCE STATEMENT Computational studies predict richer and persisting spatiotemporal patterns of spontaneous activity in neuronal networks with neuron clustering. To test this, we created networks of varying architecture in vitro. Supporting these predictions, the generation and spatiotemporal patterns of propagation were most variable in networks with intermediate clustering and lowest in uniform networks. Grid-like clustering, on the other hand, facilitated spontaneous activity but led to degenerating patterns of propagation. Neurons outside clusters had weaker synaptic input than neurons within clusters, in which increased membrane potential fluctuations facilitated the initiation of synchronized spike activity. Our results thus show that the intermediate level organization of neuronal networks strongly influences the dynamics of their activity.

  • clustering
  • network structure
  • neuronal networks
  • protein kinase C
  • spatiotemporal pattern
  • spontaneous activity

Introduction

A remarkable feature of neocortical circuitry is the generation of rich spontaneous activity dynamics (Sanchez-Vives and McCormick, 2000; Logothetis et al., 2009; Yanagawa and Mogi, 2009; Sato et al., 2012), which are believed to play an important role in cortical processing and development (Shatz, 1996; Wu et al., 2008; Ringach, 2009; Kilb et al., 2011; Altwegg-Boussac et al., 2014).

Computational network models predict that modular network architectures with highly intrinsically connected subnetworks (i.e., clustered networks) are optimal in generating and sustaining network activity (Kaiser and Hilgetag, 2010; Klinshov et al., 2014) and promote firing-rate variability and state transitions (Litwin-Kumar and Doiron, 2012).

In a simplified view, cortical neurons are indeed arranged in local clusters that have high intrinsic connectivity and connect to other clusters by correlated long-range patchy projections, i.e., neurons in a cluster share common long-range target regions (Ruthazer and Stryker, 1996; Voges et al., 2010). Functionally, this has been linked to the association of neuronal ensembles with similar function, for example the representation of orientation selectivity in the visual cortex (Bosking et al., 1997; Chavane et al., 2011). However, implications of this mesoscale network architecture for spontaneous neuronal dynamics remain speculative.

We tested the theoretical predictions in readily accessible neuronal networks in vitro by modifying the spatial distributions of neurons and neurites. In an abstract way, these networks display a mesoscale architecture with some resemblance to that of the cortex. Cultured neurons spontaneously aggregate to form clusters and grow fasciculated neurites connecting them (Kriegstein and Dichter, 1983; Segev et al., 2003; Robert et al., 2012). Their functional clustering is revealed by the spatial fragmentation of the networks when excitatory synaptic transmission is diminished (Soriano et al., 2008). Functional network reconstruction likewise indicates a mixture of locally clustered and long-range connectivity (Stetter et al., 2012). These networks establish rich spontaneous dynamics (Marom and Shahaf, 2002; van Pelt et al., 2004; Wagenaar et al., 2006a) consisting of propagating synchronized bursting events (SBEs), which are similar to patterns of spontaneous activity characteristic for the developing neocortex (Katz and Shatz, 1996; Wu et al., 2008; Golshani et al., 2009).

We modified the structure of these networks by pharmacological modulation of protein kinase C (PKC), an enzyme regulating neurite growth (Dent and Meiri, 1998), branching (Quinlan and Halpain, 1996; Audesirk et al., 1997; Schrenk et al., 2002; Gundlfinger et al., 2003; Metzger, 2010), fasciculation (Itoh et al., 1989), and cell migration (Kumada and Komuro, 2004; Larsson, 2006). By exploiting the formative effects of PKC modulation, we made networks with more clustered respectively more uniform (defined as toward randomness) arrangement of neurites and somata. Chronic PKC activation (PKC+) promoted neuronal aggregation, neurite entanglement in clusters and axonal fasciculation, and reduced overall neurite densities, suggesting more local and less long-range connectivity. Inhibition of PKC (PKC−) diminished soma clustering and axonal fasciculation and increased neurite densities, suggesting more uniform connectivity.

Differences in the mesoscale architecture crucially affected the generation and spatiotemporal structure of spontaneous activity. In PKC− networks, spontaneous SBEs were elicited at significantly lower rates even though these networks were highly excitable, as indicated by high responsiveness to electrical stimulation. We suggest that reduced clustering and fasciculation diminished recurrent and convergent connectivity patterns that promote background activity integration and SBE initiation. Neurons in clustered networks and particularly those within clusters indeed received stronger synaptic inputs, increasing membrane potential (Vm) fluctuations. Widely distributed burst initiation zones (BIZs) in clustered networks established a much richer repertoire of propagating waves compared with uniform networks where SBEs were elicited only from a few hotspots. However, clustering also resulted in incomplete network recruitment during SBEs, leading to restricted activation patterns and global desynchronization.

Our results support theoretical predictions that clustered network topologies promote generation, maintenance, and richness of spontaneous activity. A balance between local and long-range connectivity furthermore seems necessary for the generation of spontaneous activity dynamics that recruit large parts of the network at the same time.

Materials and Methods

Cell culture techniques.

Primary cortical cell cultures were prepared on different microelectrode arrays (MEAs; Multi Channel Systems; electrode grid layout/pitch distance: 8 × 8/200 μm; 6 × 10/500 μm; 16 × 16/200 μm; 32 × 32/300 μm) and standard coverslips (12 mm diameter, Carl Roth). All substrates were coated with polyethylene-imine (150 μl of 0.2% aqueous solution; Sigma-Aldrich) for cell adhesion. Cell cultures were prepared following Shahaf and Marom (2001). Cortical tissue was prepared from brains of neonatal Wistar rat pups of either sex, minced with a scalpel, and transferred into PBS (Invitrogen). Tissue pieces were incubated with trypsin (isozyme mixture, 0.05%, 15 min at 37°C; Invitrogen) and proteolysis was subsequently stopped with horse serum (20%; Invitrogen). DNase (type IV, 50 μg/ml; Sigma-Aldrich) was added to eliminate cell trapping in DNA strings if needed. Cells were dissociated by trituration with a serological pipette, centrifuged (5 min, 617 × g), and resuspended in growth medium [Minimal Essential Medium supplemented with 5% heat-inactivated horse serum, 0.5–1 mm l-glutamine, 20 mm glucose, and 20 μg/ml gentamycin (all from Invitrogen); 1 ml/pup]. Cells were counted with an automated cell counter (CASY, Schärfe Systems) and seeded at 300,000 cells per culture, resulting in a density of 1500 cells/mm2 at 1 d in vitro (DIV). Networks developed in 1 ml of growth medium in a humidified incubator (5% CO2, 37°C). Animal handling and tissue preparation were done in accordance with the guidelines for animal research at the University of Freiburg.

PKC modulation and disinhibition.

PKC inhibitor Gödecke6976 (Gö6976; 1 μm; Sigma-Aldrich) and PKC agonist phorbol-12-myristate-13-acetate (PMA; 1 μm; Sigma-Aldrich) were dissolved in dimethylsulfoxide (DMSO; Sigma-Aldrich) and added to the culture medium directly after cell preparation. The maximal concentration of DMSO in the growth medium was 0.1%. GABAergic transmission was probed by acute application of the noncompetitive GABAA receptor antagonist picrotoxin (PTX; 10 μm; Sigma-Aldrich) during electrophysiological recordings. Drug washout was performed by a complete medium exchange with fresh medium after ∼4 weeks. Networks growing without PKC modulation (PKCN) were treated the same way to evaluate general washout-induced effects. Cultures were subsequently recorded repeatedly within the next days.

Morphological analyses.

Cell adherence after seeding was documented by phase contrast microscopy at 1 DIV to estimate initial cell densities. Cell positions were determined in phase contrast micrographs by automated detection (2D convolution with a Mexican hat-shaped kernel with an inner radius corresponding to the diameter of neuronal cell bodies and subsequent thresholding of the resulting image). The spatial distribution of neuronal cell bodies was reexamined after 21 DIV based on immunocytochemical staining of neuronal nuclei (NeuN; rabbit-anti-NeuN, 1:500; Abcam) and staining of all cellular nuclei (DAPI; Sigma-Aldrich). Cellular nuclei were detected automatically (similar approach as described above) and neurons were identified as the subset of cells with NeuN immunoreactivity.

To assess minimal possible distances between neuronal somata, we calculated the distance between the centers of pairs of neurons by Delaunay triangulation. Spatial clustering of cell bodies was evaluated by a modified Clark–Evans aggregation index (Clark and Evans, 1954). The clustering index (CI) was calculated as the ratio between average observed and expected (i.e., for random point patterns) nearest-neighbor distance. CI increases from fully clustered (CI = 0 if the distance between nearest neighbors is equal to the minimal distance for all neurons) to random (CI = 1) to grid-like distributions of neurons (CI > 1, depending on the cell density). Note that CI is insensitive to a potentially grid-like arrangement of clusters (i.e., the next level of structural organization), which we do not analyze in this study. Under realistic circumstances, the distances between neurons must be a least the diameter of cell bodies, which leads to longer average nearest-neighbor distances and an underestimation of clustering. To account for this, we simulated random distributions with equal density using the average cell body diameter (10 μm) as minimal possible distance by successively adding neurons to randomly drawn positions and rejecting neurons that would overlap with others (Galli-Resta et al., 1999). We further extended the analysis to assess cluster sizes by calculating the nth nearest-neighbor distances in networks and simulated random distributions.

Neurite morphology was examined by immunocytochemical staining of microtubule-associated protein 2 (MAP2; chicken-anti-MAP2; 1:500; Abcam) and phosphorylated neurofilament 200 kDa (rabbit-anti-neurofilament; 1:10; Abcam) to visualize dendritic and axonal compartments, respectively. To detect axons and dendrites, images for the respective channels were high-pass filtered. In a second step, a peak detection was performed for each row of the resulting image (peaks corresponded to neurites; detection resolution of 2 μm intervals between intersections). To avoid a bias due to strongly orientated neurites, the procedure was repeated for different orientations (in steps of 30°; implemented by image rotation). Neurite counts and intervals of all scans were pooled. Axodendritic intersections were marked at the center-of-mass of patches of colocalized axon-positive and dendrite-positive pixels (at any angle), allowing deviation of 1 pixel width (0.645 μm). To determine the degree of neurite fasciculation, we took the mean interneurite intervals divided by the SD as fasciculation index (FI). This provided a measure similar to the CI for the cell body distribution: FI increases from fully fasciculated (FI = 0; minimal distance between neurites) to random (FI = 1) to grid-like neurite arrangements (FI > 1).

Dendritic morphology was examined in sparse networks by Sholl analysis (Sholl, 1953). All morphological analyses were made with Matlab. Significance is specified by p values determined with the Student's t test (pstt) for independent samples.

Patch-clamp recording and analysis.

Patch pipettes [6.3 ± 1.4 MΩ (mean and SEM)] were filled with a intracellular solution containing the following: potassium d-gluconate (125 mm; Sigma-Aldrich), KCl (20 mm; Sigma-Aldrich), EGTA (5 mm; Carl Roth), Na2-ATP (2 mm; Carl Roth), HEPES (10 mm; Carl Roth), MgCl2 (2 mm; Sigma-Aldrich), CaCl2 (0.5 mm; Sigma-Aldrich), and biocytin (10 mg/ml, Invitrogen), adjusted with KOH to pH 7.4, and with sucrose to 320 mOsm. Patch-clamp recordings in whole-cell configuration were conducted at 37°C (PH01 perfusion heating, Multi Channel Systems; TC02 temperature controller, Multi Channel Systems) and perfusion with carbogenated (95% O2 and 5% CO2; Air Liquide) culture medium without horse serum and without Gö6976 and PMA. Data were sampled at 25 kHz (Micro1401 amplifier and Spike2 software; Cambridge Electronics Design). Up to four neurons were recorded per network for ∼30 min each. Input resistances were determined with hyperpolarizing current pulses (−50 pA). Datasets of ≥20 min were analyzed with Matlab. The SD of Vm was averaged for time bins of 1 s. The first step to detect EPSPs was to bootstrap an EPSP template from the data. Putative events were detected with a threshold set manually for Vm slope. Snippets surrounding detection time points tD were classified by hierarchical clustering and the EPSP template was subsequently derived as the average of a manually chosen cluster resembling the voltage time course of typical isolated EPSPs. Based on this template, EPSPs were detected by correlation of Vm with the onset phase of the template (−5 to 10 ms around tD) and a correlation threshold of 0.5.

Extracellular recording and analyses.

Multiunit spike activity was recorded from microelectrode arrays (MEAs; MEA1060-BC, USB-MEA256-System, and MEA 30-1024-System amplifiers; Multi Channel Systems, 25 kHz sampling frequency, 12 bit) under culture conditions (37°C, 5% CO2) and acquired with MCRack software (Multi Channel Systems; versions 3.3–4.0). Recordings of individual networks lasted ≥1 h. Action potentials (APs) were detected with a threshold set to −5 SDs of the high-pass-filtered baseline signal (Butterworth second-order high-pass filter, 200 Hz cutoff; detection dead time, 2 ms). Stimulation was controlled using a programmable stimulus generator (STG2004, Multi Channel Systems). A single stimulus consisted of a monophasic negative pulse with 400 μs width and 0.6 V amplitude. Stimulation electrode sequences were controlled by MEABench (version 1.0.16; Wagenaar and Potter, 2004). After a baseline recording period of ≥30 min of spontaneous activity, all array electrodes were stimulated consecutively at interstimulus intervals of either 6, 4, or 2 s. The sequence was repeated 20 times in succession.

Raw data from MEA recordings were imported into Matlab using MEA-Tools (Egert et al., 2002) and the Find toolbox (Meier et al., 2008). Spontaneous SBEs were detected as follows: series of spikes with consecutive interspike intervals shorter than a threshold value (100 ms) were detected as bursts. SBEs were defined as periods in which a predefined fraction of electrodes showed simultaneous bursts (10% of active sites, minimally 3, maximally 20 to keep criteria comparable between small and large MEAs). To account for buildup and fading phases of SBEs, spikes within a time windows of 25 ms before and following this SBE core were included into the SBE. Network activity was characterized by the following parameters: the rate of SBEs in the recording period, burst strength as the average number of APs per SBE and active site, average firing rate (AFR) as the gross average firing rate across all active sites, and SBE recruitment as the average fraction of active sites participating in SBEs. For experiments with PTX, we defined as control period the last 1 h section before application of PTX and excluded the first 10 min after application from the analysis to avoid transients. In stimulation experiments, we excluded the first 10 ms after stimuli to blank stimulation artifacts.

Propagation pattern analysis.

To analyze spatiotemporal propagation patterns, SBEs were classified by first-spike rank order, i.e., electrodes were ranked by the relative timing of the first spike recorded at a given electrode during SBEs. Similarity between SBEs was calculated as the Spearman correlation of the resulting electrode sequences. Groups of similar patterns were identified by correlation-matrix-based hierarchical clustering (complete linkage) using one minus the correlation value as the distance measure (Liu et al., 2012). Distances <1 indicate correlated patterns. Distances >1 indicate anticorrelated patterns. To assess the richness of patterns, we determined the number of classes yielded at a given distance threshold between 0.02 and 2. As the contribution of individual classes, we calculated the cumulative fraction of SBEs accounted for by the most frequent pattern classes for series of 100 SBEs.

SBE initiation zones.

BIZs were identified in recordings from large MEAs (16 × 16 and 32 × 32 electrode grid layouts) that spanned almost the entire network area. On these arrays, the median x and y coordinates of the first 10 active sites that fired during SBEs were used to localize the spatial SBE onset position. Note that this approach produced a slight bias of BIZs toward the array center. The spatiotemporal structure of SBE initiation process was further analyzed as firing rate relative to SBE onset time (2 ms bins) and distance to respective BIZs (300 μm bins, i.e., electrode pitch) averaged across 500 consecutive SBEs.

Results

PKC promotes cell migration and clustering of neuronal cell bodies

Neuronal networks were prepared from neonatal rat cortex with ∼600–800 neurons/mm2 surviving at 20 DIV. This variability is unlikely to be relevant for the overall level of network activity (Biffi et al., 2013). Following intense neurite outgrowth, these networks developed spontaneous activity with SBEs as early as 5 DIV.

Under PKCN conditions (i.e., with normal PKC activity), neuronal somata migrated and aggregated, establishing a soma density landscape with high-density clusters within 3 weeks (Fig. 1A). To find the minimal distance between neurons, we calculated the Euclidean distance distribution for the edges derived from a Delaunay triangulation of neuron positions. In all pharmacological conditions, we found a minimal distance of ∼10 μm between cell body centers, which corresponds to their soma diameter (Fig. 1B). The degree of clustering after development was quantified following the Clark–Evans aggregation index (Clark and Evans, 1954), i.e., the ratio of the observed mean nearest-neighbor distance in the network to that expected for a Poisson point process with the same spatial density, modified to take into account the minimal possible distance of 10 μm between soma centers (Fig. 1C). A CI of <1 denotes clustering and a CI >1 denotes dispersion of neuronal somata relative to the expectation. PKCN networks were strongly clustered at 20 DIV compared with almost random spatial distributions in networks shortly after seeding at 1 DIV (CI mean ± SEM: 20 DIV, 0.78 ± 0.03; 1 DIV, 1.0 ± 0.02; n = 6, pstt = 1.3*10−12). We further calculated the cluster index for the nth nearest neighbors to assess the average number of neurons in clusters. At 20 DIV, CI for the nth nearest neighbors (CIn) was significantly reduced approximately up to the 45th nearest neighbor (pstt < 0.01, pairwise testing), indicating typical cluster sizes in this range (Fig. 1D).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Generic network structure and activity dynamics in networks of cultured cortical neurons. A, Staining of neuronal nuclei (NeuN; red) and all cellular nuclei (DAPI; blue) in networks at 20 DIV that had developed under inhibited (PKC−), normal (PKCN), and enhanced PKC (PKC+) activity. Neurons were detected (white circles) and analyzed for their spatial distributions. Scale bars, 100 μm. B, Distributions of minimal distances between soma centers of pairs of neighboring neurons. Distances were ≥10 μm, indicating that the soma diameter was ∼10 μm (dashed line). C, Clark–Evans CI for nearest-neighbor distances indicates significantly weaker clustering in PKC− and stronger clustering in PKC+ networks compared with PKCN networks. D, CIs <1 up to the nth nearest neighbor indicate typical cluster sizes in the range of ≤60 neurons. Boxplots show median, 25th, and 75th percentiles and minimal and maximal values (excluding outliers). ***pstt = 0.001.

To modify the degree of clustering and to investigate its relevance for the spatiotemporal structure of activity in these networks, we increased respectively inhibited PKC activity during network development. Increasing PKC activity by PMA (1 μM; PKC+) during network development promoted neuronal aggregation compared with PKCN networks (CI, 0.67 ± 0.02; n = 5, pstt = 1.1*10−4) and resulted in networks with strikingly regularly spaced, well delineated clusters at 20 DIV (Fig. 1A). Although clustering was increased in PKC+ networks (Fig. 1C), typical cluster sizes were only slightly larger than in PKCN networks, with CIn significantly reduced for the first 50 nearest neighbors (pstt < 0.01 for pairwise testing; Fig. 1D). In contrast, chronic PKC inhibition by Gö6976 diminished cell migration and led to more homogeneously distributed neuronal somata (Fig. 1A) with a significantly lower degree of clustering (CI, 0.90 ± 0.04; n = 5; pstt = 3.6*10−4; Fig. 1C) compared with PKCN networks and PKC+ networks. Remaining cell aggregation resulted in areas with slightly increased density but with weak contrast to the background density of neurons. Again, these denser areas contained ∼50 neurons (Fig. 1D).

PKC promotes fasciculation and local density of neurites

Substantial neurite outgrowth started within hours after seeding. We quantified the density and arrangement of neurites in mature networks at 20 DIV based on immunohistochemical staining for MAP2 in dendrites and neurofilament in axons (Fig. 2A). Neurite densities were determined by counting neurite intersections with straight lines drawn across the network (Fig. 2B,C).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Manipulating PKC activity changes the mesoscopic network architecture. A, Staining of axons (anti-neurofilament; green) and dendrites/cell bodies (anti-MAP2; red) and NeuN (blue). The dashed square indicates the panel positions in Figure 1A and white circles indicate the detected neuronal somata. Under PKCN conditions, axons and dendrites fasciculated and cell bodies clustered. PKC inhibition resulted in more homogeneous axon and dendrite coverage, and soma distributions. Enhanced PKC activity resulted in well delineated clusters with strongly entangled dendrites and axons that interconnected clusters in bundles and looped within clusters. B, C, To quantify fasciculation and density, dendrites (B) and axons (C) were analyzed in high-pass filtered images of red and green channels, respectively, by detecting intensity peaks in image rows. Neurite detections for one row are shown in red. The FI was calculated based on the distribution of intervals between peaks. Note that neurite densities decrease from PKC− over PKCN to PKC+. D, Intersections (small white dots) of axons and dendrites provide a lower bound estimate for putative synaptic sites. The spatial distribution of intersections was concordant to the distributed or clustered arrangement of cell bodies (large circles) and neurites. Scale bars, 100 μm.

PKCN networks had high dendrite densities (125.3 ± 9.3 dendrites/mm, mean ± SEM, n = 6; Fig. 3A) with relatively homogeneous coverage (FI, 0.84 ± 0.01; Fig. 3B). The density of axons (61 ± 8 axons/mm; Fig. 3C) was lower than that of dendrites, and axons also formed bundles between clusters (FI, 0.73 ± 0.09, mean ± SEM; Fig. 3D). Note that this does not mean that axons connected exactly two adjacent clusters. Some axons spanned distances of several millimeters and could pass through several clusters.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Distribution and density of neurites in different PKC activity conditions. A, Dendrite densities in networks decreased significantly with increasing PKC activity. B, Neurite fasciculation increased with PKC activity. C, D, Changes of PKC activity levels had the same and even stronger effects on axons. E, Dendrites were analyzed by Sholl analysis in sparse cultures stained with antibodies against MAP2 at 24 DIV. PKC inhibition increased the average number of dendritic branches at any distance from the soma, thereby producing larger dendrites. Boxplots show median, 25th, and 75th percentiles and minimal and maximal values (excluding outliers). *pstt = 0.05; **pstt = 0.01; ***pstt = 0.001.

Following the clustering of somata, increasing PKC activity had a profound impact on the arrangement of neurites. The overall density of dendrites decreased in PKC+ networks (99 ± 9 dendrites/mm, pstt = 9.8*10−4, n = 5 networks). Furthermore, dendrites were densely aggregated in neuron clusters, leaving regions with lower dendrite density in between clusters (Fig. 2B). Fasciculation of dendrites was significantly increased in PKC+ networks (FI, 0.78 ± 0.04; pstt = 2.7*10−4). The density of axons likewise was significantly reduced compared with PKCN networks (49 ± 5 axons/mm, pstt = 3.3*10−2, n = 5). As in PKCN networks, axons fasciculated and formed bundles between soma clusters (FI, 0.69 ± 0.05). Within clusters, axons strongly ramified and formed local loops.

In contrast to chronic PKC stimulation, developmental PKC inhibition significantly increased dendrite and axon densities (141 ± 3 dendrites/mm, pstt = 1.2*10−2; 131 ± 7 axons/mm, pstt = 6.5*10−7, n = 4) compared with PKCN networks.

To validate the effect at the single-neuron level, we further analyzed the extent of dendritic fields by Sholl analysis in sparse cultures at 24 DIV (Fig. 3E). Chronic PKC inhibition increased dendrite branching and extent, resulting in significantly increased total dendrite length (PKCN, 1449 ± 647 μm; N = 58 neurons; PKC−, 2132 ± 816 μm; N = 78 neurons; pstt = 5.6*10−7). In addition to increasing neurite densities, PKC inhibition significantly diminished dendritic (FI, 0.92 ± 0.01, pstt = 2.5*10−3) and axonal (FI, 0.92 ± 0.002, pstt = 4.2*10−3) fasciculation, leading to a much more homogeneous neurite coverage than in PKCN networks.

Synaptic connectivity

As a lower estimate of the number of potential synapse sites in different network types, we detected intersections between axons and dendrites (Fig. 2D). PKCN networks yielded ∼40,000 intersections/mm2, corresponding to ∼100 sites per neuron. Note that due to the limited spatial resolution of the analysis, in particular in clusters and neurite bundles, this number considerably underestimated the real density of synapses but enabled estimates over large areas of the network. The more homogeneous arrangement of neurites in PKC− networks entailed a corresponding distribution of axodendritic intersections. In PKC− networks, we determined the highest number of axodendritic intersections of ∼70,000 sites/mm2, much higher than in PKCN networks and consistent with the high neurite density. In contrast, the overall density in PKC+ networks of ∼30,000 sites/mm2 was slightly lower than in PKCN networks. As a consequence of the local tangles of neurites, axodendritic intersections were spatially clustered in PKC+ networks.

Changes in the distribution of axodendritic contact sites likely influence the connectivity of the network and thus the synaptic input to individual neurons. Such changes could be critical for the initiation and propagation of spontaneous network activity. We thus analyzed synaptic dynamics in whole-cell patch-clamp recordings during spontaneous network activity between 18 and 30 DIV for neurons grown in different PKC activity conditions (number of neurons and networks: PKCN n = 45 neurons, 16 networks; PKC− n = 45 neurons, 21 networks; PKC+ n = 20 neurons, 10 networks).

Neurons in PKC− and PKCN networks had comparable input resistances (PKCN, 233 ± 14 MΩ, mean ± SEM; PKC−, 225 ± 16 MΩ). In PKC+ networks, however, input resistances were significantly lower (PKC+, 136 ± 19 MΩ; pstt = 6.5*10−4 vs PKCN). Resting potentials varied across neurons and were significantly more negative in PKC− networks (PKC+, −50.4 ± 1.7 mV, mean ± SEM; pstt = 0.08 vs PKCN; PKCN, −53.8 ± 1.0 mV; PKC−, −57.6 ± 1.2 mV, pstt = 0.02 vs PKCN). Average thresholds of AP initiation did not differ significantly between network types (PKCN, −42.5 ± 0.7 mV; PKC−, −41.9 ± 1.4 mV; PKC+, −39.3 ± 2.0 mV). To measure EPSP activity under comparable conditions, the resting potential of all neurons was set to ∼−65 mV by constant current injection.

In all networks, Vm fluctuated around the resting Vm between SBEs and underwent strong depolarization for several hundred milliseconds, during which APs were generated (Fig. 4A). These up states corresponded to SBEs measured extracellularly (Fig. 5A–C). Vm fluctuation between SBEs increased from PKC− to PKCN to PKC+ networks (SD of Vm: PKC−, 0.95 ± 0.12 mV; PKCN, 1.97 ± 0.13 mV; pstt = 8.2*10−8; PKC+, 2.25 ± 0.43 mV).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Intracellular recordings from neurons in different mesoscopic network architectures. A, Vm traces from neurons recorded in PKC−, PKCN, and PKC+ networks at 20 DIV display EPSP activity (black arrows) and strong depolarizations corresponding to SBEs. The AP threshold ∼−40 mV is indicated by the dashed line. Insets, Average shape of nonoverlapping EPSPs in between SBEs. B, C, Average amplitude distributions across neurons reveal a larger fraction of smaller EPSPs in PKC− networks than in PKCN networks. Strongly clustered PKC+ networks had a bimodal amplitude distribution with smaller and large EPSPs. D, Neurons were classified according to their local neuronal neighborhood. EPSP amplitude distributions for neurons in sparse regions did not differ significantly between conditions. E, F, Clustered neurons in PKCN and PKC+ networks had larger EPSP amplitudes than neurons in sparse regions. Boxplots show median, 25th, and 75th percentiles and minimal and maximal values (excluding outliers). *pstt = 0.05; ***pstt = 0.001.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Paired intracellular and MEA recordings. A, B, Vm traces from two neurons near an MEA electrode. C, EPSPs are indicated by vertical ticks beneath the traces and the dashed horizontal line indicates the spiking threshold at ∼−40 mV. Baseline EPSP activity is interrupted by Vm up states, during which these neurons generated spikes. D, E, Vm up states coincide with network-wide SBEs recorded with the MEA. F, Propagation patterns were defined from the order of the first spikes at each electrode after detection of SBE onset, shown here for these two consecutive SBEs. Green line in E indicates MEA electrode marked in C.

To further determine how the embedding of individual neurons into the network affected Vm variability, we classified neurons according to their membership in clusters. This revealed that Vm of neurons within clusters fluctuated more than that of neurons in sparse regions in PKCN networks (SD of Vm: sparse, 1.60 ± 0.14 mV; n = 14; cluster, 2.44 ± 0.27 mV; n = 16; pstt = 5.9*10−3) and PKC+ networks (sparse, 1.19 ± 0.18 mV; n = 6; cluster, 2.87 ± 0.60 mV; n = 13). PKC− networks had no clear clusters to assess such dependence.

As a potential source of these differences, we assessed the effective strength of synaptic inputs in different network types by analyzing identifiable EPSPs during interburst periods. In PKCN networks, EPSPs occurred at 8.1 ± 0.6 Hz (mean ± SEM). EPSPs amplitudes were approximately lognormally distributed (2.3 ± 0.2 mV; Fig. 4B,C). Occasional fast subthreshold excursions of Vm >10 mV possibly corresponded to multisynaptic input from the same neuron. In PKC− networks, EPSPs rates in between SBEs were significantly lower (5.0 ± 0.5 Hz, pstt = 7.8*10−5) with much lower EPSP amplitudes (1.5 ± 0.1 mV, pstt = 6.1*10−4). In contrast, in the strongly clustered PKC+ networks, EPSPs rates were 9.1 ± 0.6 Hz (pstt = 3.7*10−6 vs PKC−) with an average amplitude of 2.7 ± 0.5 mV (pstt = 2.2*10−3 vs PKC− networks).

Since synaptic properties could depend on the local embedding of neurons into the network, we again classified neurons according to their position within or outside a cluster. In neurons outside of clusters, EPSP amplitude distributions and the average EPSP amplitudes were similar in all network types (PKCN, 1.8 ± 0.2 mV; PKC−, 1.5 ± 0.1 mV; PKC+, 1.5 ± 0.1 mV; Fig. 4D,F). Neurons located within clusters, however, had higher average EPSP amplitudes (PKCN, 2.8 ± 0.4 mV; pstt = 0.02; PKC+, 3.3 ± 0.6 mV; Fig. 4E,F), whereas EPSP rates in clustered and nonclustered neurons did not differ significantly.

Clustering promotes spontaneous bursting activity

Vm dynamics in individual neurons reflect the ongoing AP dynamics in recurrent networks subsampled by their particular afferent population. To assess the spontaneous network-wide dynamics, we recorded extracellular spike activity over extended time periods with MEAs. SBEs corresponded to up states measured intracellularly (Fig. 5A–C) and represented the dominating spontaneous activity pattern, clearly standing out from the uncorrelated background activity in all network types. At onset, SBEs typically propagated across the network with a traveling wave front (Fig. 5F). In all networks, typically >80% of spikes were part of SBEs (Fig. 6A; recording at 25–35 DIV; PKCN, 95 ± 1%; n = 22; PKC−, 85 ± 3%, n = 18; PKC+, 97 ± 1%, n = 7). Consistent with the higher rate of Vm up states observed intracellularly, the frequency of SBEs increased with the degree of clustering in networks (Fig. 6B).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Spontaneous SBE activity at 30 DIV in networks recorded with MEAs. A, In all network architectures, spontaneous network activity consisted of SBEs and low levels of spiking in between. Right, Zoom into SBEs marked by the arrow on the left. The spatiotemporal structure of SBE activity differed remarkably in networks with different structures. Spontaneous SBE rates, gross activity levels, and spatiotemporal variability increased with the degree of clustering in networks from PKC− to PKC+ networks. The strength of neuronal recruitment during SBEs and the bursting strength decreased along this axis. B, Spontaneous SBE rates. C, Bursting strength, i.e., the mean number of APs per site and SBE. D, AFRs, i.e., the average AP frequency across all site with activity in the 1 h recording session. E, SBE recruitment, i.e., the average fraction of active sites recruited during SBEs. F, The irregularity of SBE initiation given as the coefficient of variance (CV) of inter-SBE intervals. G, Variability in the intensity of neuronal firing during SBEs as average CV of bursts strengths (#APs/SBE) determined for individual sites. Boxplots show median, 25th, and 75th percentiles and minimal and maximal values (excluding outliers). *pstt = 0.05; **pstt = 0.01; ***pstt = 0.001.

SBE rates increased from PKC− networks (2.0 ± 0.3 SBE/min; pstt = 1.4*10−8 vs PKCN) and PKCN networks (13.0 ± 1.4 SBE/min) to PKC+ networks (29.5 ± 5.0 SBE/min, pstt = 1.8*10−4 vs PKCN).

Furthermore, the temporal variability of SBE initiation significantly increased with the degree of clustering in networks (Fig. 6F). In some networks, SBE rates fluctuated strongly, with periods of strongly increased bursting comparable to superbursts described by Wagenaar et al. (2006a). These superbursts occurred in almost all PKC+ networks, occasionally in the moderately clustered PKCN networks, and rarely in the PKC− networks.

SBEs typically did not involve the entire network. The fraction of the network recruited in individual SBEs decreased with the degree of clustering (Fig. 6E). In PKC− networks, a significantly higher fraction of all active sites participated in individual SBEs, compared with the clustered PKCN and PKC+ networks (PKCN, 67 ± 4%; PKC−, 86 ± 3%; pstt = 1.8*10−3; PKC+, 54 ± 5%; pstt = 1.8*10−3 vs PKC−).

Gross AFRs (Fig. 6D) across recording sessions of ≥1 h were significantly lower in PKC− networks (0.35 ± 0.05 Hz, pstt = 6.6*10−4) compared with PKCN networks (1.50 ± 0.28 Hz). Despite their much higher SBE rate, AFRs in PKC+ networks (1.49 ± 0.15 Hz) were not significantly higher than in PKCN networks. This was due to an inverse correlation between SBE rate and burst strength, i.e., the average number of spikes per SBE and site. In consequence, high SBE rates entailed weaker bursts (PKC−, 10.9 ± 1.9 spikes/site; PKCN, 7.8 ± 1.3 spikes/site; PKC+, 3.5 ± 0.7 spikes/site; pstt = 0.026 vs PKC−; Fig. 6B). The bursting strength was more variable in clustered networks where individual sites displayed varying numbers of spikes during SBEs (Fig. 6G).

Network metastructure influences the richness of activity patterns

The higher SBE variability in terms of burst size, network recruitment, SBE intervals, and the occurrence of superbursts in PKCN and PKC+ networks suggests that clustering promoted the richness of activity, which in this context refers to the spectrum of temporal and spatial patterns of spike activity and their pathways of propagation within a network. Based on MEA recordings, we identified SBE propagation patterns as the rank order of burst onset at individual electrodes (delay of first spike after SBE onset). Such order-based representation of SBE activity, compared with temporal representations that preserve exact delays, was demonstrated to be more robust for reconstructing SBE origins and less dependent on detecting the true first spike of a pattern (Shahaf et al., 2008).

SBEs typically spread with a traveling wave front, which was irregular in PKCN networks, but smooth and regular in the more homogeneous PKC− networks as well as in the strongly clustered PKC+ networks. Propagation pattern similarity was determined as the Spearman correlation of first-spike rank-order sequences across SBEs (Fig. 7A). In PKCN networks, correlation-matrix-based hierarchical clustering (Liu et al., 2012) did not reveal conspicuous pattern classes. In contrast, PKC− as well as PKC+ networks showed clear clusters of highly similar propagation patterns (Fig. 7A–C) with a high intraclass correlation and low or negative (distances >1) interclass correlation.

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Diversity in SBE propagation patterns. A, Spearman correlation matrices for first-spike rank-order sequences in SBEs (N = 100 SBEs, 1 network per condition). B, Corresponding dendrograms for complete linkage clustering. PKC− networks typically displayed a few distinct propagation pattern classes, while PKCN networks were marked by a continuum of gradually differing patterns. Strongly clustered PKC+ networks had many distinct pattern classes. C, Average propagation patterns for classes for a threshold distance of 0.2 (corresponding to a Spearman correlation coefficient of 0.8; B, red line). D, The number of SBE classes depended on the distance threshold applied to group SBEs. PKC− and PKC+ networks formed significantly fewer classes for low and intermediate distance thresholds. E, Fraction of all SBEs of the N-most frequently occurring SBE classes. In PKC− networks a significantly smaller number of pattern classes dominated the activity. Bars above the panel indicates the range of pairs with pstt < 0.01 tested versus PKCN.

To compare the diversity of these propagation patterns across networks, we determined the number of classes with patterns correlated above a given threshold (Fig. 7D). Intermediate distance thresholds between 0.02 and 0.8 (correlation coefficient, 0.2–0.98), yielded significantly fewer classes in PKC− networks than in PKCN networks (p < 0.01), indicating that SBE propagation patterns clustered in few groups with high intragroup similarity in PKC− networks. While this might suggest that diversity increases with clustering, the more strongly clustered network structure found in PKC+ networks did not further promote the diversity of propagation patterns. Instead pattern diversity was reduced and more similar to PKC− patterns with respect to their intraclass similarity, although this was statistically significant only for low classification thresholds between 0.02 and 0.16 (correlation coefficient, 0.84–0.98; average number of classes resulting from classification of 100 SBEs with a threshold of 0.2: PKCN, 59.3 ± 2.7; N = 33; PKC−, 17.5 ± 2.8; pstt = 3.6*10−13; N = 17; PKC+, 46.4 ± 8.0; N = 8). In PKC− networks, SBE activity was dominated by a smaller number of pattern classes than in PKCN and PKC+ networks (Fig. 7E). On average, the five most frequent classes accounted for 89.1 ± 2.8% (pstt = 5.5*10−11 vs PKCN) of the SBEs in PKC− networks but accounted for only 44.5 ± 3.5% in PKCN and 53.8 ± 8.6% in PKC+ networks.

Clustering promotes spatial variability of SBE initiation

SBEs typically initiated locally in BIZs, here defined by the electrodes at which spikes assigned to a burst were detected first (Fig. 8A). With large MEAs that covered networks almost fully, we identified the position of such BIZs. In all networks, BIZs clustered spatially, indicating hot spots of SBE initiation, and had a propensity to be located at the network boundary (Fig. 8B). Consistent with the low diversity of SBE propagation patterns, BIZs were more clustered in PKC− networks than in PKCN and PKC+ networks. BIZs alternated in eliciting SBE (Fig. 8B, right), with clusters of BIZs (hotspots) dominating SBE initiation in particular in PKC− networks. In all network types, spiking activity increased in the vicinity (∼1.5 mm) of the BIZ before SBE onset (Fig. 8C). Interestingly, in PKC+ networks, activity could persist between SBEs in particular in the BIZ.

Figure 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 8.

BIZs and network hotspots. A, SBEs were typically elicited in local areas of the network and spread with a traveling wave front. Note the smoother wave fronts in PKC− and PKC+ networks compared with those in PKCN networks. In these recordings, the MEA covered most of the network. PKC− and PKCN networks were recorded with 1000 electrodes PKC+ networks were recorded with 256 electrodes. Color indicates rank order of the first spikes recorded for an SBE at a given electrode. Left, One sample pattern. Right, Samples of typical patterns placed according to their BIZ, regardless of their frequency of occurrence (downsampled to 8 × 8 tiles with 8 × 8 and 4 × 4 electrodes). B, SBE initiation in the PKC− networks was typically dominated by few BIZs at the network boundary. Right, The trajectory depicts the temporal sequence of BIZs eliciting SBEs (300 SBEs) by connecting BIZ of successive SBEs. In PKCN and PKC+ networks, SBEs were elicited from many distributed and alternating BIZs. Scale bars: A, B, 100 μm. C, Spatiotemporal firing-rate histogram relative to SBE onset (0 ms) and BIZ (0 mm). In all network types, activity accumulated within a distance of 1–2 mm to the BIZ before SBE onset (averages of 500 consecutive SBEs at 30 DIV).

Influence of inhibition on spontaneous dynamics

The influence of inhibition in a network depends on the connectivity of inhibitory neurons and thus may differ in networks with different architecture. To test the influence of inhibition on spontaneous activity dynamics, we disinhibited networks with PTX applied at a concentration blocking ionotropic GABAergic transmission (10 μm; Krishek et al., 1996). PTX application increased spontaneous activity (AFRs) in all network types with the strongest effect in PKC− networks (Student's t test for paired samples: PKCN, +119 ± 33%; N = 18; p = 4.9*10−3; PKC−, +496 ± 108%; N = 18; p = 1.1*10−5; PKC+, +88 ± 45%; N = 8; p = 1.1*10−2; Fig. 9A). More detailed analysis showed that this was due to significantly stronger bursting (spikes/site) during SBEs (PKCN, +980 ± 378%; p = 5.1*10−4; PKC−, +598 ± 96%; p = 4.6*10−7; PKC+, +724 ± 101%; p = 0.069; Fig. 9B). In a result counterintuitive to the expectation that increased excitability following disinhibition promotes SBE initiation, SBE rates significantly decreased in PKCN and PKC+ networks after adding PTX, and remained at the low baseline levels in PKC− networks (PKCN, −62 ± 4%, p = 2.7*10−5; PKC−, −1 ± 12%; p = 0.16; PKC+, −72 ± 6%; p = 1.1*10−2; Fig. 9C).

Figure 9.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 9.

Effects of disinhibition and electrical stimulation. A–C, Blocking GABA receptors (10 μm PTX) in PKCN networks produced stronger bursts at the expense of lowered SBE rates. AFRs significantly increased despite of this trade-off. PTX application likewise increased SBE strength in PKC− networks did not affect SBE rates. In consequence, AFRs significantly increased. D–F, To assess the maximal activity level supported by a network, we electrically stimulated at alternating sites. PKC− networks were highly responsive, producing SBE rates close to the rate of stimulation (dashed line). Higher SBE rates (spontaneous or evoked) were at the cost of weaker bursts with fewer APs per site and SBE but still resulted in significantly increased AFRs. PKCN networks were less responsive to stimulation than PKC− networks and showed no significant increase of AFRs. Boxplots show median, 25th, and 75th percentiles and minimal and maximal values (excluding outliers). **p = 0.01; ***p = 0.001.

These observations suggest that burst strength-dependent network depression and recovery could limit SBE rates. Lower activity, more negative resting Vm,, and the pronounced increase of SBE strength upon disinhibition in PKC− networks could result from increased or more effective inhibition compared with PKCN networks. This then would also predict that their excitability should be lower. To investigate whether PKC− networks were limited by reduced excitability, we stimulated PKC− and PKCN networks electrically.

PKC− networks support high activity levels

To test whether PKC− networks can support higher SBE rates, we stimulated a separate set of networks electrically with interstimulus intervals of either 6, 4, or 2 s at alternating stimulation sites, corresponding to the range of spontaneous SBE activation intervals and the alternating BIZs in PKCN networks. In PKC− networks, electrical stimulation significantly increased AFRs (6 s interval +583 ± 245%, p = 0.030; 4 s interval +757 ± 332%, p = 0.023; 2 s interval +859 ± 344%, p = 0.011; N = 3; mean ± SEM; Student's t test for paired samples; Fig. 9D) and SBE rates (6 s interval, +762 ± 327%; p = 6.1*10−4; 4 s interval, +1117 ± 466%; p = 6.1*10−4; 2 s interval, +1978 ± 755%; p = 9.6*10−4; spontaneous and evoked SBEs; Fig. 9F) beyond the level of spontaneous activity in PKCN networks. Moreover, PKC− networks showed almost complete entrainment of the SBE activity to the stimulation with very high response probability (evoked responses within 100 ms for 91.3, 87.6, and 81.4% of the stimuli given in series of 1200 stimuli at 6, 4, and 2 s intervals, respectively; N = 4) and highly reproducible propagation patterns (similar to spontaneous propagation patterns). Stimulation only slightly increased gross average activity levels in PKCN networks (6 s interval, +31 ± 6%; 4 s interval, +44 ± 17%; 2 s interval, +73 ± 25%; differences were not statistically significant; N = 3; Fig. 9D) and SBE rate (6 s interval, +76 ± 39%; 4 s interval, +95 ± 57%; 2 s interval, +194 ± 112%; differences were not statistically significant; Fig. 9F). PKCN networks were less responsive to stimulation (response probabilities: 6 s interval, 48.8%; 4 s interval, 48.4%; 2 s interval, 45.6%) and responses were much more variable. Additional SBEs started spontaneously between stimuli. Eliciting additional SBEs by stimulation reduced SBE strength in both network types but differences between control and stimulation sessions were not statistically significant for PKCN (PKCN: 6 s interval, −22 ± 12%; 4 s interval, −21 ± 12%; 2 s interval, −32 ± 13%; PKC−: 6 s interval, −21 ± 3%; p = 0.009; 4 s interval, −31 ± 5%; p = 0.018; 2 s interval, −55 ± 6%; p = 0.003; Fig. 9E). In summary, SBE rates in PKC− networks were neither limited by their capacity to sustain higher rates nor by GABAergic inhibition. Electrically driven activity showed a dynamical trade-off between SBE rate and SBE strength.

Removal of PKC modulators after chronic exposure

To identify direct effects of the PKC modulators on the network activity, we replaced the medium with fresh drug-free medium after 4 weeks in all three PKC activity conditions (Fig. 10). Washout-induced changes of burst activity were similar in all groups. Relative differences in spontaneous activity between PKC activity conditions were preserved for ≥3 DIV after drug removal. This supports the idea that the findings above were related to changes in network structure rather than to direct drug effects.

Figure 10.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 10.

Effects of drug washout following chronic treatment. A–F, Complete medium exchange with fresh drug-free medium after 30 d of chronic treatment with PKC modulators produced no systematic changes of the activity dynamics. Changes in PKCN networks indicate an effect of the washout itself, likely due to the strong perturbation or the replacement with fresh medium. Changes in PKC-modulated networks were in the same range or smaller than those in PKCN networks. See Figure 6 for further details.

Discussion

Computational network models suggest that clustered connectivity is beneficial for generation, maintenance (Kaiser and Hilgetag, 2010; Stetter et al., 2012; Klinshov et al., 2014), and variability of spontaneous activity dynamics in neuronal networks (Litwin-Kumar and Doiron, 2012). We addressed this issue in cortical networks in vitro in which we manipulated the degree of clustering and neurite bundling.

Manipulating mesoscale network architecture by PKC modulation

Networks in vitro establish community structures (Girvan and Newman, 2002; Sporns, 2013) exhibiting local neuron clusters with strongly entangled neurites, neurite bundles interconnecting clusters, and long-range axons (Kriegstein and Dichter, 1983; Shefi et al., 2002; Bettencourt et al., 2007; Stetter et al., 2012). We manipulated this architecture by modulating PKC activity. The control of cell migration and neurite outgrowth by PKC has been discussed previously (Itoh et al., 1989; Dent and Meiri, 1998; Kumada and Komuro, 2004; Larsson, 2006; Metzger, 2010). Diminishing PKC activity during network development produced remarkably uniform networks with more randomly distributed somata and less fasciculated neurites than in PKCN networks. The opposite effects resulted from stimulation of PKC, which promoted spatial aggregation of somata and local neurite sprouting within regularly spaced and well delineated clusters.

As with cerebellar Purkinje cells (Metzger, 2010), PKC inhibition increased neurite fields and neurite densities. This presumably increases the probability of a synaptic connection between any pair of neurons, predicting increased availability of long-range connections. Conversely, lower neurite and, in particular, lower axon densities in PKC+ networks should reduce long-range connectivity. Pronounced tangles of neurites within clusters, in turn, likely promote highly recurrent local connectivity. In addition, axon bundles formed mainly between neighboring clusters.

We therefore hypothesize that diminished PKC activity shifted the network toward uniform connectivity, whereas increasing PKC activity promoted recurrent connectivity within and convergent connectivity between local clusters at the expense of long-range connectivity.

Mesoscale network architecture and synaptic dynamics

Anatomical analyses in vivo showed that statistically approximately a quarter of the contact sites between axon and dendrites are realized as functional connections (Stepanyants et al., 2002). Neurite architecture thus is a strong determinant of the realized synaptic connectivity. In our networks, synapses were indeed more uniformly distributed in PKC− networks. Based on these observations and reports that the strength of synapses inversely scales with their number in vitro (Wilson et al., 2007) and in vivo (Turrigiano, 2008), smaller average EPSP amplitudes in PKC− networks might be explained by the higher neurite density resulting in more synaptic connections. Although PKC has been implicated in synaptic plasticity (Saito and Shirai, 2002; Boehm et al., 2006), it is apparently not essential for activity-dependent synaptic plasticity due to other functionally redundant kinases (Wang and Kelly, 1996; Herring and Nicoll, 2016). If so, neurons should retain activity-dependent synaptic plasticity even with reduced PKC activity.

Large EPSP amplitudes in neurons associated with clusters are consistent with findings by Cohen et al. (2008) that in isolated clusters with <10 neurons, the stimulation of single neurons is sufficient to activate the entire ensemble while this is not possible in larger networks. EPSP amplitudes measured in neurons outside of clusters were indeed comparable across conditions, suggesting that their effective synaptic weights were not determined by the respective overall changes in network structure but reflected their local network embedding.

The larger EPSPs in clusters could also be explained by a higher propensity for multiple synaptic contacts between neuron pairs and the highly correlated presynaptic firing attributed to clusters (Helias et al., 2014). The precise synchrony among such connections, however, makes it impossible to distinguish between unitary EPSPs from single EPSPs and compound EPSPs from sister release sites.

Clustering promotes spontaneous SBE generation

Vm fluctuations strongly determine neuronal firing rates (Kuhn et al., 2004), which in a recurrent network translate back again into neuronal input statistics. Giugliano et al. (2004) hypothesized that SBE initiation is a threshold-gated process following activity integration on the network level analogous to the AP generation in neurons. Network architectures that increase background firing-rate fluctuations should thus promote SBE initiation. Along with larger EPSPs and higher EPSP frequencies, average Vm and Vm fluctuations increased with the degree of clustering in networks. Neurons consequently were more likely to reach the spiking threshold, which can be sufficient to activate an entire local cluster (Cohen et al., 2008). Fasciculation of axons that extend from clusters conveys highly correlated output to common target clusters, promoting focused long-range excitation. Neuronal clusters could thus act as local amplifiers or relay stations that promote generation of spontaneous activity and reliable transmission of activity within a network. Consistent with this, SBE rates significantly increased with the degree of clustering in networks. Furthermore, SBEs were elicited at widely distributed sites in clustered networks, suggesting that many local subnetworks exist that trigger SBEs. PKC− networks typically had few hotspots of SBE initiation. We propose that diminished neurite fasciculation fostered divergent rather than convergent connectivity patterns in these networks, reducing correlated input. Interestingly, hotspots were mostly located at the network boundary, which was most prominent in more uniform PKC− networks. In agreement with Gritsun et al. (2012), the boundary may enforce inward connections, thus promoting recurrent and convergent connectivity. Clusters may introduce additional BIZs by pronounced locally recurrent connectivity.

How changes in the network architecture affect overall inhibition levels and activity dynamics is difficult to predict. In PKC− networks excitation–inhibition balance could be more homogeneous than in clustered networks, where this could vary locally. This would introduce hot (and cold) spots that could further contribute to SBE initiation. Although recent studies find that homeostatic synaptic plasticity or spike-frequency adaptation could compensate for local excitation–inhibition imbalance (Barral and D Reyes, 2016; Landau et al., 2016), it is not clear whether this would fully flatten the excitation–inhibition landscape.

Global versus local synchrony

SBEs dominated in all networks, regardless of their particular mesoscale structure. However, SBE rate modulations and the occurrence of super bursts (Wagenaar et al., 2006b) increased with the degree of clustering and neurite fasciculation, which is consistent with network simulations showing that clustered topologies promote firing-rate variability (Litwin-Kumar and Doiron, 2012). In PKC− networks, activity reliably recruited most neurons once SBEs were initiated. In the strongly clustered PKC+ networks, however, SBEs often remained confined to local network areas, as was observed in heavily clustered networks forming without adhesive growth substrates (Teller et al., 2014). We suggest that strong local dynamics in clusters and weak long-distance coupling entails global desynchronization analogous to effects described for weakly coupled oscillators (Strogatz and Mirollo, 1991; Park et al., 2006). Global synchrony thus should depend on the balance between local and global connectivity. Moderately clustered PKCN networks were optimal in terms of allowing high SBE rates conjointly with strong network recruitment in SBEs and seem to possess sufficient recurrent connectivity to sustain activity as well as sufficient long-range connectivity to recruit distant parts of the network.

SBE initiation: a buildup or release process?

The rate of SBE succession may depend on recovery from synaptic depression and a buildup process that involves the amplification of activity in recurrent networks beyond a critical threshold (Tabak et al., 2010). To gain insight into the underlying mechanism, we tested networks by disinhibition and electrical stimulation. As described earlier (Weihberger et al., 2013), disinhibition only slightly increased overall spontaneous activity levels in PKCN networks, since longer bursts were counterbalanced by lower SBE rates. This trade-off suggests that PKCN networks operate within a regime of activity that exploits most resources. In agreement with this, overall activity levels could not be increased much by external stimulation. In PKC− networks, amplified bursting through disinhibition was not balanced by lower SBE rates, suggesting rate limitation by the initiation process rather than by synaptic resources. Consistent with this assumption, electrical stimulation significantly increased activity levels and SBE rates.

Clustering promotes richness of activity patterns

In accordance with Soriano et al. (2008), we found discrete network areas that had a high propensity to trigger SBEs. These BIZs were frequent and widely distributed in clustered networks and less abundant in homogenous PKC− networks, which is consistent with the idea that convergent and recurrent connectivity patterns promote burst initiation. In all networks, SBEs spread with approximately circular wave fronts. Interestingly, PKC+ networks displayed more regular wave fronts than PKCN networks, rather similar to those observed in PKC− networks. One could argue that this is because regular spacing of clusters in PKC+ networks likewise introduces uniformity at the mesoscale level. Since the local direction of propagation depended mainly on the BIZ location, distributed BIZs in clustered networks fostered a rich repertoire of different spatiotemporal propagation patterns.

Conclusion

We showed that modulation of PKC during network development can change the balance between locally clustered and long-range connectivity. Consistent with network simulations, our data indicate that clustering promotes local activity generation but leads to reduced global synchrony. Intriguingly, the richness of structure and activity patterns in PKCN networks are suggestive of small-world and scale-free topologies with their computational advantages. We propose that the coevolution of spontaneous activity levels, richness of activity patterns, and community network structure reflects a fundamental principle of neuronal self-organization.

Footnotes

  • This work was supported by BrainLinks-BrainTools, Cluster of Excellence funded by the German Research Foundation (Grant EXC 1086), and Bernstein Focus: Neurotechnology Freiburg–Tuebingen (FKZ 01GQ0420). We thank Sarah Jarvis, Ehsan Safavieh, and Oliver Weihberger for helpful discussions and gratefully acknowledge technical assistance from Ute Riede, Patrick Pauli, Alexander Giffey, Hanna Kuhn, Nila Mönig, and Patrick Ringwald.

  • Correspondence should be addressed to Samora Okujeni at the above address. E-mail: okujeni{at}bcf.uni-freiburg.de.

References

  1. ↵
    1. Altwegg-Boussac T,
    2. Chavez M,
    3. Mahon S,
    4. Charpier S
    (2014) Excitability and responsiveness of rat barrel cortex neurons in the presence and absence of spontaneous synaptic activity in vivo. J Physiol 592:3577–3595. doi:10.1113/jphysiol.2013.270561 pmid:24732430
    OpenUrlCrossRefPubMed
  2. ↵
    1. Audesirk G,
    2. Cabell L,
    3. Kern M
    (1997) Modulation of neurite branching by protein phosphorylation in cultured rat hippocampal neurons. Brain Res Dev Brain Res 102:247–260. doi:10.1016/S0165-3806(97)00100-4 pmid:9352107
    OpenUrlCrossRefPubMed
  3. ↵
    1. Barral J,
    2. D Reyes A
    (2016) Synaptic scaling rule preserves excitatory-inhibitory balance and salient neuronal network dynamics. Nat Neurosci 19:1690–1696. doi:10.1038/nn.4415 pmid:27749827
    OpenUrlCrossRefPubMed
  4. ↵
    1. Bettencourt LM,
    2. Stephens GJ,
    3. Ham MI,
    4. Gross GW
    (2007) Functional structure of cortical neuronal networks grown in vitro. Phys Rev E Stat Nonlin Soft Matter Phys 75:021915. doi:10.1103/PhysRevE.75.021915 pmid:17358375
    OpenUrlCrossRefPubMed
  5. ↵
    1. Biffi E,
    2. Regalia G,
    3. Menegon A,
    4. Ferrigno G,
    5. Pedrocchi A
    (2013) The influence of neuronal density and maturation on network activity of hippocampal cell cultures: a methodological study. PLoS One 8:e83899. doi:10.1371/journal.pone.0083899 pmid:24386305
    OpenUrlCrossRefPubMed
  6. ↵
    1. Boehm J,
    2. Kang MG,
    3. Johnson RC,
    4. Esteban J,
    5. Huganir RL,
    6. Malinow R
    (2006) Synaptic incorporation of AMPA receptors during LTP is controlled by a PKC phosphorylation site on GluR1. Neuron 51:213–225. doi:10.1016/j.neuron.2006.06.013 pmid:16846856
    OpenUrlCrossRefPubMed
  7. ↵
    1. Bosking WH,
    2. Zhang Y,
    3. Schofield B,
    4. Fitzpatrick D
    (1997) Orientation selectivity and the arrangement of horizontal connections in tree shrew striate cortex. J Neurosci 17:2112–2127. pmid:9045738
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Chavane F,
    2. Sharon D,
    3. Jancke D,
    4. Marre O,
    5. Frégnac Y,
    6. Grinvald A
    (2011) Lateral spread of orientation selectivity in V1 is controlled by intracortical cooperativity. Front Syst Neurosci 5:4. doi:10.3389/fnsys.2011.00004 pmid:21629708
    OpenUrlCrossRefPubMed
  9. ↵
    1. Clark PJ,
    2. Evans FC
    (1954) Distance to nearest neighbor as a measure of spatial relationships in populations. Ecol Soc Am Stable 35:445–453.
    OpenUrl
  10. ↵
    1. Cohen E,
    2. Ivenshitz M,
    3. Amor-Baroukh V,
    4. Greenberger V,
    5. Segal M
    (2008) Determinants of spontaneous activity in networks of cultured hippocampus. Brain Res 1235:21–30. doi:10.1016/j.brainres.2008.06.022 pmid:18602907
    OpenUrlCrossRefPubMed
  11. ↵
    1. Dent EW,
    2. Meiri KF
    (1998) Distribution of phosphorylated GAP-43 (neuromodulin) in growth cones directly reflects growth cone behavior. J Neurobiol 35:287–299. doi:10.1002/(SICI)1097-4695(19980605)35:3%3C287::AID-NEU6%3E3.0.CO%3B2-V pmid:9622012
    OpenUrlCrossRefPubMed
  12. ↵
    1. Egert U,
    2. Knott T,
    3. Schwarz C,
    4. Nawrot M,
    5. Brandt A,
    6. Rotter S,
    7. Diesmann M
    (2002) MEA-Tools: an open source toolbox for the analysis of multi-electrode data with MATLAB. J Neurosci Methods 117:33–42. doi:10.1016/S0165-0270(02)00045-6 pmid:12084562
    OpenUrlCrossRefPubMed
  13. ↵
    1. Galli-Resta L,
    2. Novelli E,
    3. Kryger Z,
    4. Jacobs GH,
    5. Reese BE
    (1999) Modelling the mosaic organization of rod and cone photoreceptors with a minimal-spacing rule. Eur J Neurosci 11:1461–1469. doi:10.1046/j.1460-9568.1999.00555.x pmid:10103140
    OpenUrlCrossRefPubMed
  14. ↵
    1. Girvan M,
    2. Newman ME
    (2002) Community structure in social and biological networks. Proc Natl Acad Sci U S A 99:7821–7826. doi:10.1073/pnas.122653799 pmid:12060727
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Giugliano M,
    2. Darbon P,
    3. Arsiero M,
    4. Lüscher H-R,
    5. Streit J
    (2004) Single-neuron discharge properties and network activity in dissociated cultures of neocortex. J Neurophysiol 92:977–996. pmid:15044515
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Golshani P,
    2. Gonçalves JT,
    3. Khoshkhoo S,
    4. Mostany R,
    5. Smirnakis S,
    6. Portera-Cailliau C
    (2009) Internally mediated developmental desynchronization of neocortical network activity. J Neurosci 29:10890–10899. doi:10.1523/JNEUROSCI.2012-09.2009 pmid:19726647
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Gritsun TA,
    2. le Feber J,
    3. Rutten WL
    (2012) Growth dynamics explain the development of spatiotemporal burst activity of young cultured neuronal networks in detail. PLoS One 7:e43352. doi:10.1371/journal.pone.0043352 pmid:23028450
    OpenUrlCrossRefPubMed
  18. ↵
    1. Gundlfinger A,
    2. Kapfhammer JP,
    3. Kruse F,
    4. Leitges M,
    5. Metzger F
    (2003) Different regulation of Purkinje cell dendritic development in cerebellar slice cultures by protein kinase Calpha and -beta. J Neurobiol 57:95–109. doi:10.1002/neu.10259 pmid:12973831
    OpenUrlCrossRefPubMed
  19. ↵
    1. Helias M,
    2. Tetzlaff T,
    3. Diesmann M
    (2014) The correlation structure of local neuronal networks intrinsically results from recurrent dynamics. PLoS Comput Biol 10:e1003428. doi:10.1371/journal.pcbi.1003428 pmid:24453955
    OpenUrlCrossRefPubMed
  20. ↵
    1. Herring BE,
    2. Nicoll RA
    (2016) Long-term potentiation: from CaMKII to AMPA receptor trafficking. Annu Rev Physiol 78:351–365. doi:10.1146/annurev-physiol-021014-071753 pmid:26863325
    OpenUrlCrossRefPubMed
  21. ↵
    1. Itoh K,
    2. Asou H,
    3. Ikarashi Y,
    4. Maruyama Y
    (1989) Morphological changes and neural cell adhesion molecule expression in mouse cerebrum primary cultures following long-term exposure to phorbol ester. Neurosci Res 6:350–357. doi:10.1016/0168-0102(89)90027-8 pmid:2725992
    OpenUrlCrossRefPubMed
  22. ↵
    1. Kaiser M,
    2. Hilgetag CC
    (2010) Optimal hierarchical modular topologies for producing limited sustained activation of neural networks. Front Neuroinform 4:8. doi:10.3389/fninf.2010.00008 pmid:20514144
    OpenUrlCrossRefPubMed
  23. ↵
    1. Katz LC,
    2. Shatz CJ
    (1996) Synaptic activity and the construction of cortical circuits. Science 274:1133–1138. doi:10.1126/science.274.5290.1133 pmid:8895456
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Kilb W,
    2. Kirischuk S,
    3. Luhmann HJ
    (2011) Electrical activity patterns and the functional maturation of the neocortex. Eur J Neurosci 34:1677–1686. doi:10.1111/j.1460-9568.2011.07878.x pmid:22103424
    OpenUrlCrossRefPubMed
  25. ↵
    1. Klinshov VV,
    2. Teramae JN,
    3. Nekorkin VI,
    4. Fukai T
    (2014) Dense neuron clustering explains connectivity statistics in cortical microcircuits. PLoS One 9:e94292. doi:10.1371/journal.pone.0094292 pmid:24732632
    OpenUrlCrossRefPubMed
  26. ↵
    1. Kriegstein AR,
    2. Dichter MA
    (1983) Morphological classification of rat cortical neurons in cell culture. J Neurosci 3:1634–1647. pmid:6875660
    OpenUrlAbstract
  27. ↵
    1. Krishek BJ,
    2. Moss SJ,
    3. Smart TG
    (1996) A functional comparison of the antagonists bicuculline and picrotoxin at recombinant GABAA receptors. Neuropharmacology 35:1289–1298. doi:10.1016/s0028-3908(96)00089-5 pmid:9014144
    OpenUrlCrossRefPubMed
  28. ↵
    1. Kuhn A,
    2. Aertsen A,
    3. Rotter S
    (2004) Neuronal integration of synaptic input in the fluctuation-driven regime. J Neurosci 24:2345–2356. doi:10.1523/JNEUROSCI.3349-03.2004 pmid:15014109
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Kumada T,
    2. Komuro H
    (2004) Completion of neuronal migration regulated by loss of Ca(2+) transients. Proc Natl Acad Sci U S A 101:8479–8484. doi:10.1073/pnas.0401000101 pmid:15150416
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Landau ID,
    2. Egger R,
    3. Dercksen VJ,
    4. Oberlaender M,
    5. Sompolinsky H
    (2016) The impact of structural heterogeneity on excitation-inhibition balance in cortical networks. Neuron 92:1106–1121. doi:10.1016/j.neuron.2016.10.027 pmid:27866797
    OpenUrlCrossRefPubMed
  31. ↵
    1. Larsson C
    (2006) Protein kinase C and the regulation of the actin cytoskeleton. Cell Signal 18:276–284. doi:10.1016/j.cellsig.2005.07.010 pmid:16109477
    OpenUrlCrossRefPubMed
  32. ↵
    1. Litwin-Kumar A,
    2. Doiron B
    (2012) Slow dynamics and high variability in balanced cortical networks with clustered connections. Nat Neurosci 15:1498–1505. doi:10.1038/nn.3220 pmid:23001062
    OpenUrlCrossRefPubMed
  33. ↵
    1. Liu X,
    2. Zhu XH,
    3. Qiu P,
    4. Chen W
    (2012) A correlation-matrix-based hierarchical clustering method for functional connectivity analysis. J Neurosci Methods 211:94–102. doi:10.1016/j.jneumeth.2012.08.016 pmid:22939920
    OpenUrlCrossRefPubMed
  34. ↵
    1. Logothetis NK,
    2. Murayama Y,
    3. Augath M,
    4. Steffen T,
    5. Werner J,
    6. Oeltermann A
    (2009) How not to study spontaneous activity. Neuroimage 45:1080–1089. doi:10.1016/j.neuroimage.2009.01.010 pmid:19344685
    OpenUrlCrossRefPubMed
  35. ↵
    1. Marom S,
    2. Shahaf G
    (2002) Development, learning and memory in large random networks of cortical neurons: lessons beyond anatomy. Q Rev Biophys 35:63–87. pmid:11997981
    OpenUrlCrossRefPubMed
  36. ↵
    1. Meier R,
    2. Egert U,
    3. Aertsen A,
    4. Nawrot MP
    (2008) FIND—a unified framework for neural data analysis. Neural Netw 21:1085–1093. doi:10.1016/j.neunet.2008.06.019 pmid:18692360
    OpenUrlCrossRefPubMed
  37. ↵
    1. Metzger F
    (2010) Molecular and cellular control of dendrite maturation during brain development. Curr Mol Pharmacol 3:1–11. doi:10.2174/1874467211003010001 pmid:20030626
    OpenUrlCrossRefPubMed
  38. ↵
    1. Park K,
    2. Lai YC,
    3. Gupte S,
    4. Kim JW
    (2006) Synchronization in complex networks with a modular structure. Chaos 16:015105. doi:10.1063/1.2154881 pmid:16599771
    OpenUrlCrossRefPubMed
  39. ↵
    1. Quinlan EM,
    2. Halpain S
    (1996) Emergence of activity-dependent, bidirectional control of microtubule-associated protein MAP2 phosphorylation during postnatal development. J Neurosci 16:7627–7637. pmid:8922419
    OpenUrlAbstract/FREE Full Text
  40. ↵
    1. Ringach DL
    (2009) Spontaneous and driven cortical activity: implications for computation. Curr Opin Neurobiol 19:439–444. doi:10.1016/j.conb.2009.07.005 pmid:19647992
    OpenUrlCrossRefPubMed
  41. ↵
    1. Robert F,
    2. Cloix JF,
    3. Hevor T
    (2012) Ultrastructural characterization of rat neurons in primary culture. Neuroscience 200:248–260. doi:10.1016/j.neuroscience.2011.10.002 pmid:22079571
    OpenUrlCrossRefPubMed
  42. ↵
    1. Ruthazer ES,
    2. Stryker MP
    (1996) The role of activity in the development of long-range horizontal connections in area 17 of the ferret. J Neurosci 16:7253–7269. pmid:8929433
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Saito N,
    2. Shirai Y
    (2002) Protein kinase C gamma (PKC gamma): function of neuron specific isotype. J Biochem 132:683–687. doi:10.1093/oxfordjournals.jbchem.a003274 pmid:12417016
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Sanchez-Vives MV,
    2. McCormick DA
    (2000) Cellular and network mechanisms of rhythmic recurrent activity in neocortex. Nat Neurosci 3:1027–1034. doi:10.1038/79848 pmid:11017176
    OpenUrlCrossRefPubMed
  45. ↵
    1. Sato TK,
    2. Nauhaus I,
    3. Carandini M
    (2012) Traveling waves in visual cortex. Neuron 75:218–229. doi:10.1016/j.neuron.2012.06.029 pmid:22841308
    OpenUrlCrossRefPubMed
  46. ↵
    1. Schrenk K,
    2. Kapfhammer JP,
    3. Metzger F
    (2002) Altered dendritic development of cerebellar Purkinje cells in slice cultures from protein kinase Cgamma-deficient mice. Neuroscience 110:675–689. doi:10.1016/S0306-4522(01)00559-0 pmid:11934475
    OpenUrlCrossRefPubMed
  47. ↵
    1. Segev R,
    2. Benveniste M,
    3. Shapira Y,
    4. Ben-Jacob E
    (2003) Formation of electrically active clusterized neural networks. Phys Rev Lett 90:168101. doi:10.1103/PhysRevLett.90.168101 pmid:12732015
    OpenUrlCrossRefPubMed
  48. ↵
    1. Shahaf G,
    2. Marom S
    (2001) Learning in networks of cortical neurons. J Neurosci 21:8782–8788. pmid:11698590
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Shahaf G,
    2. Eytan D,
    3. Gal A,
    4. Kermany E,
    5. Lyakhov V,
    6. Zrenner C,
    7. Marom S
    (2008) Order-based representation in random networks of cortical neurons. PLoS Comput Biol 4:e1000228. doi:10.1371/journal.pcbi.1000228 pmid:19023409
    OpenUrlCrossRefPubMed
  50. ↵
    1. Shatz CJ
    (1996) Emergence of order in visual system development. Proc Natl Acad Sci U S A 93:602–608. doi:10.1073/pnas.93.2.602 pmid:8570602
    OpenUrlAbstract/FREE Full Text
  51. ↵
    1. Shefi O,
    2. Golding I,
    3. Segev R,
    4. Ben-Jacob E,
    5. Ayali A
    (2002) Morphological characterization of in vitro neuronal networks. Phys Rev E Stat Nonlin Soft Matter Phys 66:021905. doi:10.1103/PhysRevE.66.021905 pmid:12241212
    OpenUrlCrossRefPubMed
  52. ↵
    1. Sholl DA
    (1953) Dendritic organization in the neurons of the visual and motor cortices of the cat. J Anat 87:387–406. pmid:1244622
    OpenUrlPubMed
  53. ↵
    1. Soriano J,
    2. Rodríguez Martínez M,
    3. Tlusty T,
    4. Moses E
    (2008) Development of input connections in neural cultures. Proc Natl Acad Sci U S A 105:13758–13763. doi:10.1073/pnas.0707492105 pmid:18772389
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Sporns O
    (2013) Structure and function of complex brain networks. Dialogues Clin Neurosci 15:247–262. pmid:24174898
    OpenUrlPubMed
  55. ↵
    1. Stepanyants A,
    2. Hof PR,
    3. Chklovskii DB
    (2002) Geometry and structural plasticity of synaptic connectivity. Neuron 34:275–288. doi:10.1016/S0896-6273(02)00652-9 pmid:11970869
    OpenUrlCrossRefPubMed
  56. ↵
    1. Stetter O,
    2. Battaglia D,
    3. Soriano J,
    4. Geisel T
    (2012) Model-free reconstruction of excitatory neuronal connectivity from calcium imaging signals. PLoS Comput Biol 8:e1002653. doi:10.1371/journal.pcbi.1002653 pmid:22927808
    OpenUrlCrossRefPubMed
  57. ↵
    1. Strogatz SH,
    2. Mirollo RE
    (1991) Stability of incoherence in a population of coupled oscillators. J Stat Phys 63:613–635. doi:10.1007/BF01029202
    OpenUrlCrossRef
  58. ↵
    1. Tabak J,
    2. Mascagni M,
    3. Bertram R
    (2010) Mechanism for the universal pattern of activity in developing neuronal networks. J Neurophysiol 103:2208–2221. doi:10.1152/jn.00857.2009 pmid:20164396
    OpenUrlAbstract/FREE Full Text
  59. ↵
    1. Teller S,
    2. Granell C,
    3. De Domenico M,
    4. Soriano J,
    5. Gómez S,
    6. Arenas A
    (2014) Emergence of assortative mixing between clusters of cultured neurons. PLoS Comput Biol 10:e1003796. doi:10.1371/journal.pcbi.1003796 pmid:25188377
    OpenUrlCrossRefPubMed
  60. ↵
    1. Turrigiano GG
    (2008) The self-tuning neuron: synaptic scaling of excitatory synapses. Cell 135:422–435. doi:10.1016/j.cell.2008.10.008 pmid:18984155
    OpenUrlCrossRefPubMed
  61. ↵
    1. van Pelt J,
    2. Wolters PS,
    3. Corner MA,
    4. Rutten WL,
    5. Ramakers GJ
    (2004) Long-term characterization of firing dynamics of spontaneous bursts in cultured neural networks. IEEE Trans Biomed Eng 51:2051–2062. doi:10.1109/TBME.2004.827936 pmid:15536907
    OpenUrlCrossRefPubMed
  62. ↵
    1. Voges N,
    2. Guijarro C,
    3. Aertsen A,
    4. Rotter S
    (2010) Models of cortical networks with long-range patchy projections. J Comput Neurosci 28:137–154. doi:10.1007/s10827-009-0193-z pmid:19866352
    OpenUrlCrossRefPubMed
  63. ↵
    1. Wagenaar DA,
    2. Potter SM
    (2004) A versatile all-channel stimulator for electrode arrays, with real-time control. J Neural Eng 1:39–45. doi:10.1088/1741-2560/1/1/006 pmid:15876621
    OpenUrlCrossRefPubMed
  64. ↵
    1. Wagenaar DA,
    2. Pine J,
    3. Potter SM
    (2006a) An extremely rich repertoire of bursting patterns during the development of cortical cultures. BMC Neurosci 7:11. doi:10.1186/1471-2202-7-11 pmid:16464257
    OpenUrlCrossRefPubMed
  65. ↵
    1. Wagenaar DA,
    2. Nadasdy Z,
    3. Potter SM
    (2006b) Persistent dynamic attractors in activity patterns of cultured neuronal networks. Phys Rev E Stat Nonlin Soft Matter Phys 73:051907. doi:10.1103/PhysRevE.73.051907 pmid:16802967
    OpenUrlCrossRefPubMed
  66. ↵
    1. Wang JH,
    2. Kelly PT
    (1996) The balance between postsynaptic Ca(2+)-dependent protein kinase and phosphatase activities controlling synaptic strength. Learn Mem 3:170–181. doi:10.1101/lm.3.2-3.170 pmid:10456087
    OpenUrlAbstract/FREE Full Text
  67. ↵
    1. Weihberger O,
    2. Okujeni S,
    3. Mikkonen JE,
    4. Egert U
    (2013) Quantitative examination of stimulus-response relations in cortical networks in vitro. J Neurophysiol 109:1764–1774. doi:10.1152/jn.00481.2012 pmid:23274313
    OpenUrlAbstract/FREE Full Text
  68. ↵
    1. Wilson NR,
    2. Ty MT,
    3. Ingber DE,
    4. Sur M,
    5. Liu G
    (2007) Synaptic reorganization in scaled networks of controlled size. J Neurosci 27:13581–13589. doi:10.1523/JNEUROSCI.3863-07.2007 pmid:18077670
    OpenUrlAbstract/FREE Full Text
  69. ↵
    1. Wu JY,
    2. Huang X,
    3. Zhang C
    (2008) Propagating waves of activity in the neocortex: what they are, what they do. Neuroscientist 14:487–502. doi:10.1177/1073858408317066 pmid:18997124
    OpenUrlAbstract/FREE Full Text
  70. ↵
    1. Yanagawa T,
    2. Mogi K
    (2009) Analysis of ongoing dynamics in neural networks. Neurosci Res 64:177–184. doi:10.1016/j.neures.2009.02.011 pmid:19428698
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 37 (14)
Journal of Neuroscience
Vol. 37, Issue 14
5 Apr 2017
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Mesoscale Architecture Shapes Initiation and Richness of Spontaneous Network Activity
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Mesoscale Architecture Shapes Initiation and Richness of Spontaneous Network Activity
Samora Okujeni, Steffen Kandler, Ulrich Egert
Journal of Neuroscience 5 April 2017, 37 (14) 3972-3987; DOI: 10.1523/JNEUROSCI.2552-16.2017

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Mesoscale Architecture Shapes Initiation and Richness of Spontaneous Network Activity
Samora Okujeni, Steffen Kandler, Ulrich Egert
Journal of Neuroscience 5 April 2017, 37 (14) 3972-3987; DOI: 10.1523/JNEUROSCI.2552-16.2017
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • clustering
  • network structure
  • neuronal networks
  • protein kinase C
  • spatiotemporal pattern
  • spontaneous activity

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • Prelimbic Cortex Activity during a Distress Tolerance Task Predicts Cocaine-Seeking Behavior in Male, But Not Female Rats
  • 2-AG-Mediated Control of GABAergic Signaling Is Impaired in a Model of Epilepsy
  • Identification of the Acid-Sensitive Site Critical for Chloral Hydrate (CH) Activation of the Proton-Activated Chloride Channel
Show more Research Articles

Development/Plasticity/Repair

  • Spatiotemporal Developmental Gradient of Thalamic Morphology, Microstructure, and Connectivity fromthe Third Trimester to Early Infancy
  • Loss of motor cortical inputs to the red nucleus after central nervous system disorders in non-human primates
  • Astrocytes transplanted during early postnatal development integrate, mature, and survive long-term in mouse cortex
Show more Development/Plasticity/Repair
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.