Search

Abstract The hippocampus and entorhinal cortex exhibit rich oscillatory patterns critical for cognitive functions. In the hippocampal region CA1, specific gamma-frequency oscillations, timed at different phases of the ongoing theta rhythm, are hypothesized to facilitate the integration of information from varied sources and contribute to distinct cognitive processes. Here, we show that gamma elements -a multidimensional characterization of transient gamma oscillatory episodes- occur at any frequency or phase relative to the ongoing theta rhythm across all CA1 layers in male mice. Despite their low power and stochastic-like nature, individual gamma elements still carry behavior-related information and computational modeling suggests that they reflect neuronal firing. Our findings challenge the idea of rigid gamma sub-bands, showing that behavior shapes ensembles of irregular gamma elements that evolve with learning and depend on hippocampal layers. Widespread gamma diversity, beyond randomness, may thus reflect complexity, likely functional but invisible to classic average-based analyses.

In carefully designed experimental paradigms, cognitive scientists interpret the mean event-related potentials (ERP) in terms of cognitive operations. However, the huge signal variability from one trial to the next, questions the representability of such mean events. We explored here whether this variability is an unwanted noise, or an informative part of the neural response. We took advantage of the rapid changes in the visual system during human infancy and analyzed the variability of visual responses to central and lateralized faces in 2-to 6-month-old infants compared to adults using high-density electroencephalography (EEG). We observed that neural trajectories of individual trials always remain very far from ERP components, only moderately bending their direction with a substantial temporal jitter across trials. However, single trial trajectories displayed characteristic patterns of acceleration and deceleration when approaching ERP components, as if they were under the active influence of steering forces causing transient attraction and stabilization. These dynamic events could only partly be accounted for by induced microstate transitions or phase reset phenomena. Importantly, these structured modulations of response variability, both between and within trials, had a rich sequential organization, which in infants, was modulated by the task difficulty and age. Our approaches to characterize Event Related Variability (ERV) expand on classic ERP analyses and provide the first evidence for the functional role of ongoing neural variability in human infants.

The variability in synaptic connectivity observed at the cerebellar granule cell - Purkinje cell connection in mice accounts for motor behavior traits at the individual level, suggesting that cerebellar networks encode internal models underlying individual-specific motor adaptation.

AbstractNeural computation is associated with the emergence, reconfiguration, and dissolution of cell assemblies in the context of varying oscillatory states. Here, we describe the complex spatiotemporal dynamics of cell assemblies through temporal network formalism. We use a sliding window approach to extract sequences of networks of information sharing among single units in hippocampus and entorhinal cortex during anesthesia and study how global and node-wise functional connectivity properties evolve through time and as a function of changing global brain state (theta vs. slow-wave oscillations). First, we find that information sharing networks display, at any time, a core-periphery structure in which an integrated core of more tightly functionally interconnected units links to more loosely connected network leaves. However the units participating to the core or to the periphery substantially change across time windows, with units entering and leaving the core in a smooth way. Second, we find that discrete network states can be defined on top of this continuously ongoing liquid core-periphery reorganization. Switching between network states results in a more abrupt modification of the units belonging to the core and is only loosely linked to transitions between global oscillatory states. Third, we characterize different styles of temporal connectivity that cells can exhibit within each state of the sharing network. While inhibitory cells tend to be central, we show that, otherwise, anatomical localization only poorly influences the patterns of temporal connectivity of the different cells. Furthermore, cells can change temporal connectivity style when the network changes state. Altogether, these findings reveal that the sharing of information mediated by the intrinsic dynamics of hippocampal and entorhinal cortex cell assemblies have a rich spatiotemporal structure, which could not have been identified by more conventional time- or state-averaged analyses of functional connectivity.

Electrophysiology recordings are frequently affected by artifacts (e.g., subject motion or eye movements), which reduces the number of available trials and affects the statistical power. When artifacts are unavoidable and data are scarce, signal reconstruction algorithms that allow for the retention of sufficient trials become crucial. Here, we present one such algorithm that makes use of large spatiotemporal correlations in neural signals and solves the low-rank matrix completion problem, to fix artifactual entries. The method uses a gradient descent algorithm in lower dimensions to learn the missing entries and provide faithful reconstruction of signals. We carried out numerical simulations to benchmark the method and estimate optimal hyperparameters for actual EEG data. The fidelity of reconstruction was assessed by detecting event-related potentials (ERP) from a highly artifacted EEG time series from human infants. The proposed method significantly improved the standardized error of the mean in ERP group analysis and a between-trial variability analysis compared to a state-of-the-art interpolation technique. This improvement increased the statistical power and revealed significant effects that would have been deemed insignificant without reconstruction. The method can be applied to any time-continuous neural signal where artifacts are sparse and spread out across epochs and channels, increasing data retention and statistical power.

The hippocampal formation is one of the brain systems in which the functional roles of coordinated oscillations in information representation and communication are better studied. Within this circuit, neuronal oscillations are conceived as a mechanism to precisely coordinate upstream and downstream neuronal ensembles, underlying dynamic exchange of information. Within a global reference framework provided by theta (θ) oscillations, different gamma-frequency (γ) carriers would temporally segregate information originating from different sources, thereby allowing networks to disambiguate convergent inputs. Two γ sub-bands were thus defined according to their frequency (slow γ, 30–80 Hz; medium γ, 60–120 Hz) and differential power distribution across CA1 dendritic layers. According to this prevalent model, layer-specific γ oscillations in CA1 would reliably identify the temporal dynamics of afferent inputs and may therefore aid in identifying specific memory processes (encoding for medium γ vs. retrieval for slow γ). However, this influential view, derived from time-averages of either specific γ sub-bands or different projection methods, might not capture the complexity of CA1 θ-γ interactions. Recent studies investigating γ oscillations at the θ cycle timescale have revealed a more dynamic and diverse landscape of θ-γ motifs, with many θ cycles containing multiple γ bouts of various frequencies. To properly capture the hippocampal oscillatory complexity, we have argued in this review that we should consider the entirety of the data and its multidimensional complexity. This will call for a revision of the actual model and will require the use of new tools allowing the description of individual γ bouts in their full complexity.

A critical challenge in current research on Alzheimer’s disease (AD) is to clarify the relationship between network dysfunction and the emergence of subtle memory deficits in itspreclinical stage. The AppNL-F/MAPT double knock-in (dKI) model with humanized β-amyloid peptide (Aβ) and tau was used to investigate both memory and network dysfunctions at an early stage. Young male dKI mice (2 to 6 months) were tested in three tasks taxing different aspects of recognition memory affected in preclinical AD. An early deficit first appeared in the object-place association task at the age of 4 months, when increased levels of β-CTF and Aβ were detected in both the hippocampus and the medial temporal cortex, and tau pathology was found only in the medial temporal cortex. Object-place task-dependent c-Fos activation was then analyzed in 22 subregions across the medial prefrontal cortex, claustrum, retrosplenial cortex, and medial temporal lobe. Increased c-Fos activation was detected in the entorhinal cortex and the claustrum of dKI mice. During recall, network efficiency was reduced across cingulate regions with a major disruption of information flow through the retrosplenial cortex. Our findings suggest that early perirhinal-entorhinal pathology is associated with abnormal activity which may spread to downstream regions such as the claustrum, the medial prefrontal cortex and ultimately the key retrosplenial hub which relays information from frontal to temporal lobes. The similarity between our findings and those reported in preclinical stages of AD suggests that the AppNL-F/MAPT dKI model has a high potential for providing key insights into preclinical AD.

Flexibility is a key feature of psychological health, allowing the individual to dynamically adapt to changing environmental demands, which is impaired in many psychiatric disorders like obsessive–compulsive disorder (OCD). Adequately responding to varying demands requires the brain to switch between different patterns of neural activity, which are represented by different brain network configurations (functional connectivity patterns). Here, we operationalize neural flexibility as the dissimilarity between consecutive connectivity matrices of brain regions (jump length). In total, 132 fMRI scans were obtained from 17 patients that were scanned four to five times during inpatient psychotherapy, and from 17 controls that were scanned at comparable time intervals. Significant negative correlations were found between the jump lengths and the symptom severity scores of OCD, depression, anxiety, and stress, suggesting that high symptom severity corresponds to inflexible brain functioning. Further analyses revealed that impaired reconfiguration (pattern stability) of the brain seems to be more related to general psychiatric impairment rather than to specific symptoms, e.g., of OCD or depression. Importantly, the group × time interaction of a repeated measures ANOVA was significant, as well as the post-hoc paired t-tests of the patients (first vs. last scan). The results suggest that psychotherapy is able to significantly increase the neural flexibility of patients. We conclude that psychiatric symptoms like anxiety, stress, depression, and OCD are associated with an impaired adaptivity of the brain. In general, our results add to the growing evidence that dynamic functional connectivity captures meaningful properties of brain functioning.

Functional Connectivity (FC) during resting-state or task conditions is not static but inherently dynamic. Yet, there is no consensus on whether fluctuations in FC may resemble isolated transitions between discrete FC states rather than continuous changes. This quarrel hampers advancing the study of dynamic FC. This is unfortunate as the structure of fluctuations in FC can certainly provide more information about developmental changes, aging, and progression of pathologies. We merge the two perspectives and consider dynamic FC as an ongoing network reconfiguration, including a stochastic exploration of the space of possible steady FC states. The statistical properties of this random walk deviate both from a purely “order-driven” dynamics, in which the mean FC is preserved, and from a purely “randomness-driven” scenario, in which fluctuations of FC remain uncorrelated over time. Instead, dynamic FC has a complex structure endowed with long-range sequential correlations that give rise to transient slowing and acceleration epochs in the continuous flow of reconfiguration. Our analysis for fMRI data in healthy elderly revealed that dynamic FC tends to slow down and becomes less complex as well as more random with increasing age. These effects appear to be strongly associated with age-related changes in behavioural and cognitive performance.

Dynamic Functional Connectivity (dFC) in the resting state (rs) is considered as a correlate of cognitive processing. Describing dFC as a flow across morphing connectivity configurations, our notion of dFC speed quantifies the rate at which FC networks evolve in time. Here we probe the hypothesis that variations of rs dFC speed and cognitive performance are selectively interrelated within specific functional subnetworks.In particular, we focus on Sleep Deprivation (SD) as a reversible model of cognitive dysfunction. We found that whole-brain level (global) dFC speed significantly slows down after 24h of SD. However, the reduction in global dFC speed does not correlate with variations of cognitive performance in individual tasks, which are subtle and highly heterogeneous. On the contrary, we found strong correlations between performance variations in individual tasks –including Rapid Visual Processing (RVP, assessing sustained visual attention)– and dFC speed quantified at the level of functional sub-networks of interest. Providing a compromise between classic static FC (no time) and global dFC (no space), modular dFC speed analyses allow quantifying a different speed of dFC reconfiguration independently for sub-networks overseeing different tasks. Importantly, we found that RVP performance robustly correlates with the modular dFC speed of a characteristic frontoparietal module.

At the macroscale, the brain operates as a network of interconnected neuronal populations, which display coordinated rhythmic dynamics that support interareal communication. Understanding how stimulation of different brain areas impacts such activity is important for gaining basic insights into brain function and for further developing therapeutic neurmodulation. However, the complexity of brain structure and dynamics hinders predictions regarding the downstream effects of focal stimulation. More specifically, little is known about how the collective oscillatory regime of brain network activity-in concert with network structure-affects the outcomes of perturbations. Here, we combine human connectome data and biophysical modeling to begin filling these gaps. By tuning parameters that control collective system dynamics, we identify distinct states of simulated brain activity and investigate how the distributed effects of stimulation manifest at different dynamical working points. When baseline oscillations are weak, the stimulated area exhibits enhanced power and frequency, and due to network interactions, activity in this excited frequency band propagates to nearby regions. Notably, beyond these linear effects, we further find that focal stimulation causes more distributed modifications to interareal coherence in a band containing regions' baseline oscillation frequencies. Importantly, depending on the dynamical state of the system, these broadband effects can be better predicted by functional rather than structural connectivity, emphasizing a complex interplay between anatomical organization, dynamics, and response to perturbation. In contrast, when the network operates in a regime of strong regional oscillations, stimulation causes only slight shifts in power and frequency, and structural connectivity becomes most predictive of stimulation-induced changes in network activity patterns. In sum, this work builds upon and extends previous computational studies investigating the impacts of stimulation, and underscores the fact that both the stimulation site, and, crucially, the regime of brain network dynamics, can influence the network-wide responses to local perturbations.

The capability of cortical regions to flexibly sustain an "ignited" state of activity has been discussed in relation to conscious perception or hierarchical information processing. Here, we investigate how the intrinsic propensity of different regions to get ignited is determined by the specific topological organisation of the structural connectome. More specifically, we simulated the resting-state dynamics of mean-field whole-brain models and assessed how dynamic multistability and ignition differ between a reference model embedding a realistic human connectome, and alternative models based on a variety of randomised connectome ensembles. We found that the strength of global excitation needed to first trigger ignition in a subset of regions is substantially smaller for the model embedding the empirical human connectome. Furthermore, when increasing the strength of excitation, the propagation of ignition outside of this initial core-which is able to self-sustain its high activity-is way more gradual than for any of the randomised connectomes, allowing for graded control of the number of ignited regions. We explain both these assets in terms of the exceptional weighted core-shell organisation of the empirical connectome, speculating that this topology of human structural connectivity may be attuned to support enhanced ignition dynamics.

Functional Connectivity, describing the interaction between brain regions beyond their anatomical interconnection, is highly dynamic even when no task is performed (“resting state”) and it remains a methodological challenge to properly describe its changes in time without strong assumptions. We have developed a framework to describe the dynamics of Functional Connectivity (dFC) estimated from brain activity time-series as a as a smooth reconfiguration process, combining “liquid” and “coordinated” aspects. Our framework considers dFC as a complex random walk in the space of possible functional networks. Unlike other previous approaches, our method does not require the explicit extraction of discrete connectivity states but tracks changes in a continuous time fashion. • We introduced several dFC random walk metrics. First, dFC speed analyses extract the distribution of the time-resolved rate of reconfiguration of FC along time. These distributions have a clear peak (typical dFC speed, that can already serve as a biomarker) and fat tails (denoting deviations from Gaussianity that can be detected by suitable scaling analyses of FC network streams). • Second, meta-connectivity (MC) analyses identify groups of functional links whose fluctuations co-vary in time and that define veritable dFC modules organized along specific dFC meta-hub controllers (differing from conventional FC modules and hubs). The decomposition of whole-brain dFC by MC allows performing dFC speed analyses separately for each of the detected dFC modules. • We present here blocks and pipelines for dFC random walk analyses that are made easily available through a dedicated MATLABⓇ toolbox (dFCwalk), openly downloadable. Although we applied such analyses mostly to fMRI resting state data, in principle our methods can be extended to any type of neural activity (from Local Field Potentials to EEG, MEG, fNIRS, etc.) or even non-neural time-series.

We investigated interactions within chimera states in a phase oscillator network with two coupled subpopulations. To quantify interactions within and between these subpopulations, we estimated the corresponding (delayed) mutual information that—in general—quantifies the capacity or the maximum rate at which information can be transferred to recover a sender's information at the receiver with a vanishingly low error probability. After verifying their equivalence with estimates based on the continuous phase data, we determined the mutual information using the time points at which the individual phases passed through their respective Poincaré sections. This stroboscopic view on the dynamics may resemble, e.g., neural spike times, that are common observables in the study of neuronal information transfer. This discretization also increased processing speed significantly, rendering it particularly suitable for a fine-grained analysis of the effects of experimental and model parameters. In our model, the delayed mutual information within each subpopulation peaked at zero delay, whereas between the subpopulations it was always maximal at non-zero delay, irrespective of parameter choices. We observed that the delayed mutual information of the desynchronized subpopulation preceded the synchronized subpopulation. Put differently, the oscillators of the desynchronized subpopulation were “driving” the ones in the synchronized subpopulation. These findings were also observed when estimating mutual information of the full phase trajectories. We can thus conclude that the delayed mutual information of discrete time points allows for inferring a functional directed flow of information between subpopulations of coupled phase oscillators.

Flexible information routing underlies the function of many biological and artificial networks. Here, the authors present a theoretical framework that shows how information can be flexibly routed across networks using collective reference dynamics and how local changes may induce remote rerouting.

Peer review is the cornerstone of scholarly publishing and it is essential that peer reviewers are appointed on the basis of their expertise alone. However, it is difficult to check for any bias in the peer-review process because the identity of peer reviewers generally remains confidential. Here, using public information about the identities of 9000 editors and 43000 reviewers from the Frontiers series of journals, we show that women are underrepresented in the peer-review process, that editors of both genders operate with substantial same-gender preference (homophily), and that the mechanisms of this homophily are gender-dependent. We also show that homophily will persist even if numerical parity between genders is reached, highlighting the need for increased efforts to combat subtler forms of gender bias in scholarly publishing.

Tuning curves are the functions that relate the responses of sensory neurons to various values within one continuous stimulus dimension (such as the orientation of a bar in the visual domain or the frequency of a tone in the auditory domain). They are commonly determined by fitting a model e.g. a Gaussian or other bell-shaped curves to the measured responses to a small subset of discrete stimuli in the relevant dimension. However, as neuronal responses are irregular and experimental measurements noisy, it is often difficult to determine reliably the appropriate model from the data. We illustrate this general problem by fitting diverse models to representative recordings from area MT in rhesus monkey visual cortex during multiple attentional tasks involving complex composite stimuli. We find that all models can be well-fitted, that the best model generally varies between neurons and that statistical comparisons between neuronal responses across different experimental conditions are affected quantitatively and qualitatively by specific model choices. As a robust alternative to an often arbitrary model selection, we introduce a model-free approach, in which features of interest are extracted directly from the measured response data without the need of fitting any model. In our attentional datasets, we demonstrate that data-driven methods provide descriptions of tuning curve features such as preferred stimulus direction or attentional gain modulations which are in agreement with fit-based approaches when a good fit exists. Furthermore, these methods naturally extend to the frequent cases of uncertain model selection. We show that model-free approaches can identify attentional modulation patterns, such as general alterations of the irregular shape of tuning curves, which cannot be captured by fitting stereotyped conventional models. Finally, by comparing datasets across different conditions, we demonstrate effects of attention that are cell- and even stimulus-specific. Based on these proofs-of-concept, we conclude that our data-driven methods can reliably extract relevant tuning information from neuronal recordings, including cells whose seemingly haphazard response curves defy conventional fitting approaches.

To understand the proximate and ultimate causes that shape acoustic communication in animals, objective characterizations of the vocal repertoire of a given species are critical, as they provide the foundation for comparative analyses among individuals, populations and taxa. Progress in this field has been hampered by a lack of standard in methodology, however. One problem is that researchers may settle on different variables to characterize the calls, which may impact on the classification of calls. More important, there is no agreement how to best characterize the overall structure of the repertoire in terms of the amount of gradation within and between call types. Here, we address these challenges by examining 912 calls recorded from wild chacma baboons (Papio ursinus). We extracted 118 acoustic variables from spectrograms, from which we constructed different sets of acoustic features, containing 9, 38, and 118 variables; as well 19 factors derived from principal component analysis. We compared and validated the resulting classifications of k-means and hierarchical clustering. Datasets with a higher number of acoustic features lead to better clustering results than datasets with only a few features. The use of factors in the cluster analysis resulted in an extremely poor resolution of emerging call types. Another important finding is that none of the applied clustering methods gave strong support to a specific cluster solution. Instead, the cluster analysis revealed that within distinct call types, subtypes may exist. Because hard clustering methods are not well suited to capture such gradation within call types, we applied a fuzzy clustering algorithm. We found that this algorithm provides a detailed and quantitative description of the gradation within and between chacma baboon call types. In conclusion, we suggest that fuzzy clustering should be used in future studies to analyze the graded structure of vocal repertoires. Moreover, the use of factor analyses to reduce the number of acoustic variables should be discouraged.

Neuronal dynamics are fundamentally constrained by the underlying structural network architecture, yet much of the details of this synaptic connectivity are still unknown even in neuronal cultures in vitro. Here we extend a previous approach based on information theory, the Generalized Transfer Entropy, to the reconstruction of connectivity of simulated neuronal networks of both excitatory and inhibitory neurons. We show that, due to the model-free nature of the developed measure, both kinds of connections can be reliably inferred if the average firing rate between synchronous burst events exceeds a small minimum frequency. Furthermore, we suggest, based on systematic simulations, that even lower spontaneous inter-burst rates could be raised to meet the requirements of our reconstruction algorithm by applying a weak spatially homogeneous stimulation to the entire network. By combining multiple recordings of the same in silico network before and after pharmacologically blocking inhibitory synaptic transmission, we show then how it becomes possible to infer with high confidence the excitatory or inhibitory nature of each individual neuron.

A systematic assessment of global neural network connectivity through direct electrophysiological assays has remained technically infeasible, even in simpler systems like dissociated neuronal cultures. We introduce an improved algorithmic approach based on Transfer Entropy to reconstruct structural connectivity from network activity monitored through calcium imaging. We focus in this study on the inference of excitatory synaptic links. Based on information theory, our method requires no prior assumptions on the statistics of neuronal firing and neuronal connections. The performance of our algorithm is benchmarked on surrogate time series of calcium fluorescence generated by the simulated dynamics of a network with known ground-truth topology. We find that the functional network topology revealed by Transfer Entropy depends qualitatively on the time-dependent dynamic state of the network (bursting or non-bursting). Thus by conditioning with respect to the global mean activity, we improve the performance of our method. This allows us to focus the analysis to specific dynamical regimes of the network in which the inferred functional connectivity is shaped by monosynaptic excitatory connections, rather than by collective synchrony. Our method can discriminate between actual causal influences between neurons and spurious non-causal correlations due to light scattering artifacts, which inherently affect the quality of fluorescence imaging. Compared to other reconstruction strategies such as cross-correlation or Granger Causality methods, our method based on improved Transfer Entropy is remarkably more accurate. In particular, it provides a good estimation of the excitatory network clustering coefficient, allowing for discrimination between weakly and strongly clustered topologies. Finally, we demonstrate the applicability of our method to analyses of real recordings of in vitro disinhibited cortical cultures where we suggest that excitatory connections are characterized by an elevated level of clustering compared to a random graph (although not extreme) and can be markedly non-local.

Anatomic connections between brain areas affect information flow between neuronal circuits and the synchronization of neuronal activity. However, such structural connectivity does not coincide with effective connectivity (or, more precisely, causal connectivity), related to the elusive question "Which areas cause the present activity of which others?". Effective connectivity is directed and depends flexibly on contexts and tasks. Here we show that dynamic effective connectivity can emerge from transitions in the collective organization of coherent neural activity. Integrating simulation and semi-analytic approaches, we study mesoscale network motifs of interacting cortical areas, modeled as large random networks of spiking neurons or as simple rate units. Through a causal analysis of time-series of model neural activity, we show that different dynamical states generated by a same structural connectivity motif correspond to distinct effective connectivity motifs. Such effective motifs can display a dominant directionality, due to spontaneous symmetry breaking and effective entrainment between local brain rhythms, although all connections in the considered structural motifs are reciprocal. We show then that transitions between effective connectivity configurations (like, for instance, reversal in the direction of inter-areal interactions) can be triggered reliably by brief perturbation inputs, properly timed with respect to an ongoing local oscillation, without the need for plastic synaptic changes. Finally, we analyze how the information encoded in spiking patterns of a local neuronal population is propagated across a fixed structural connectivity motif, demonstrating that changes in the active effective connectivity regulate both the efficiency and the directionality of information transfer. Previous studies stressed the role played by coherent oscillations in establishing efficient communication between distant areas. Going beyond these early proposals, we advance here that dynamic interactions between brain rhythms provide as well the basis for the self-organized control of this "communication-through-coherence", making thus possible a fast "on-demand" reconfiguration of global information routing modalities.

Visually induced neuronal activity in V1 displays a marked gamma-band component which is modulated by stimulus properties. It has been argued that synchronized oscillations contribute to these gamma-band activity. However, analysis of Local Field Potentials (LFPs) across different experiments reveals considerable diversity in the degree of oscillatory behavior of this induced activity. Contrast-dependent power enhancements can indeed occur over a broad band in the gamma frequency range and spectral peaks may not arise at all. Furthermore, even when oscillations are observed, they undergo temporal decorrelation over very few cycles. This is not easily accounted for in previous network modeling of gamma oscillations. We argue here that interactions between cortical layers can be responsible for this fast decorrelation. We study a model of a V1 hypercolumn, embedding a simplified description of the multi-layered structure of the cortex. When the stimulus contrast is low, the induced activity is only weakly synchronous and the network resonates transiently without developing collective oscillations. When the contrast is high, on the other hand, the induced activity undergoes synchronous oscillations with an irregular spatiotemporal structure expressing a synchronous chaotic state. As a consequence the population activity undergoes fast temporal decorrelation, with concomitant rapid damping of the oscillations in LFPs autocorrelograms and peak broadening in LFPs power spectra. We show that the strength of the inter-layer coupling crucially affects this spatiotemporal structure. We predict that layer VI inactivation should induce global changes in the spectral properties of induced LFPs, reflecting their slower temporal decorrelation in the absence of inter-layer feedback. Finally, we argue that the mechanism underlying the emergence of synchronous chaos in our model is in fact very general. It stems from the fact that gamma oscillations induced by local delayed inhibition tend to develop chaos when coupled by sufficiently strong excitation.