Feedforward and feedback frequency-dependent interactions in a large-scale laminar network of the primate cortex

The neural oscillatory mechanism underlying human brain fingerprint recognition using a portable EEG acquisition device

In recent years, brainprint recognition has emerged as a novel method of personal identity verification. Although studies have demonstrated the feasibility of this technology, some limitations hinder its further development into the society, such as insufficient efficiency (extended wear time for multi-channel EEG cap), complex experimental paradigms (more time in learning and completing experiments), and unclear neurobiological characteristics (lack of intuitive biomarkers and an inability to eliminate the impact of noise on individual differences). Overall, these limitations are due to the incomplete understanding of the underlying neural mechanisms. Therefore, this study aims to investigate the neural mechanisms behind brainwave recognition and simplify the operation process. We recorded prefrontal resting-state EEG data from 40 participants, which is followed up over nine months using a single-channel portable brainwave device. We found that portable devices can effectively and stably capture the characteristics of different subjects in the alpha band (8-13Hz) over long periods, as well as capturing their individual differences (no alpha peak, 1 alpha peak, or 2 alpha peaks). Through correlation analysis, alpha-band activity can reveal the uniqueness of the subjects compared to others within one minute. We further used a descriptive model to dissect the oscillatory and non-oscillatory components in the alpha band, demonstrating the different contributions of fine oscillatory features to individual differences (especially amplitude and bandwidth). Our study validated the feasibility of portable brainwave devices in brainwave recognition and the underlying neural oscillation mechanisms. The fine characteristics of various alpha oscillations will contribute to the accuracy of brainwave recognition, providing new insights for the development of future brainwave recognition technology.

Analysis of the alpha activity envelope in electroencephalography in relation to the ratio of excitatory to inhibitory neural activity

Alpha waves, one of the major components of resting and awake cortical activity in human electroencephalography (EEG), are known to show waxing and waning, but this phenomenon has rarely been analyzed. In the present study, we analyzed this phenomenon from the viewpoint of excitation and inhibition. The alpha wave envelope was subjected to secondary differentiation. This gave the positive (acceleration positive, Ap) and negative (acceleration negative, An) values of acceleration and their ratio (Ap-An ratio) at each sampling point of the envelope signals for 60 seconds. This analysis was performed on 36 participants with Alzheimer's disease (AD), 23 with frontotemporal dementia (FTD) and 29 age-matched healthy participants (NC) whose data were provided as open datasets. The mean values of the Ap-An ratio for 60 seconds at each EEG electrode were compared between the NC and AD/FTD groups. The AD (1.41 ±0.01 (SD)) and FTD (1.40 ±0.02) groups showed a larger Ap-An ratio than the NC group (1.38 ±0.02, p<0.05). A significant correlation between the envelope amplitude of alpha activity and the Ap-An ratio was observed at most electrodes in the NC group (Pearson's correlation coefficient, r = -0.92 ±0.15, mean for all electrodes), whereas the correlation was disrupted in AD (-0.09 ±0.21, p<0.05) and disrupted in the frontal region in the FTD group. The present method analyzed the envelope of alpha waves from a new perspective, that of excitation and inhibition, and it could detect properties of the EEG, Ap-An ratio, that have not been revealed by existing methods. The present study proposed a new method to analyze the alpha activity envelope in electroencephalography, which could be related to excitatory and inhibitory neural activity.

A large-scale neuronal network modelling study: Stimulus size modulates gamma oscillations in the primary visual cortex by long-range connections

Stimulus size modulation of neuronal firing activity is a fundamental property of the primary visual cortex. Numerous biological experiments have shown that stimulus size modulation is affected by multiple factors at different spatiotemporal scales, but the exact pathways and mechanisms remain incompletely understood. In this paper, we establish a large-scale neuronal network model of primary visual cortex with layer 2/3 to study how gamma oscillation properties are modulated by stimulus size and especially how long-range connections affect the modulation as realistic neuronal properties and spatial distributions of synaptic connections are considered. It is shown that long-range horizontal synaptic connections are sufficient to produce dimensional modulation of firing rates and gamma oscillations. In particular, with increasing grating stimulus size, the firing rate increases and then decreases, the peak frequency of gamma oscillations decreases and the spectral power increases. These are consistent with biological experimental observations. Furthermore, we explain in detail how the number and spatial distribution of long-range connections affect the size modulation of gamma oscillations by using the analysis of neuronal firing activity and synaptic current fluctuations. Our results provide a mechanism explanation for size modulation of gamma oscillations in the primary visual cortex and reveal the important and unique role played by long-range connections, which contributes to a deeper understanding of the cognitive function of gamma oscillations in visual cortex.

Sequential deactivation across the thalamus-hippocampus-mPFC pathway during loss of consciousness

How consciousness is lost in states such as sleep or anesthesia remains a mystery. To gain insight into this phenomenon, we conducted concurrent recordings of electrophysiology signals in the anterior cingulate cortex and whole-brain functional magnetic resonance imaging (fMRI) in rats exposed to graded propofol, undergoing the transition from consciousness to unconsciousness. Our results reveal that upon the loss of consciousness (LOC), as indicated by the loss of righting reflex, there is a sharp increase in low-frequency power of the electrophysiological signal. Additionally, simultaneously measured fMRI signals exhibit a cascade of deactivation across a pathway including the hippocampus, thalamus, and medial prefrontal cortex (mPFC) surrounding the moment of LOC, followed by a broader increase in brain activity across the cortex during sustained unconsciousness. Furthermore, sliding window analysis demonstrates a temporary increase in synchrony of fMRI signals across the hippocampus-thalamus-mPFC pathway preceding LOC. These data suggest that LOC might be triggered by sequential activities in the hippocampus, thalamus and mPFC, while wide-spread activity increases in other cortical regions commonly observed during anesthesia-induced unconsciousness might be a consequence, rather than a cause of LOC. Taken together, our study identifies a cascade of neural events unfolding as the brain transitions into unconsciousness, offering critical insight into the systems-level neural mechanisms underpinning LOC.

Thalamic contributions to the state and contents of consciousness

Consciousness can be conceptualized as varying along at least two dimensions: the global state of consciousness and the content of conscious experience. Here, we highlight the cellular and systems-level contributions of the thalamus to conscious state and then argue for thalamic contributions to conscious content, including the integrated, segregated, and continuous nature of our experience. We underscore vital, yet distinct roles for core- and matrix-type thalamic neurons. Through reciprocal interactions with deep-layer cortical neurons, matrix neurons support wakefulness and determine perceptual thresholds, whereas the cortical interactions of core neurons maintain content and enable perceptual constancy. We further propose that conscious integration, segregation, and continuity depend on the convergent nature of corticothalamic projections enabling dimensionality reduction, a thalamic reticular nucleus-mediated divisive normalization-like process, and sustained coherent activity in thalamocortical loops, respectively. Overall, we conclude that the thalamus plays a central topological role in brain structures controlling conscious experience.

Human connectome topology directs cortical traveling waves and shapes frequency gradients

Traveling waves and neural oscillation frequency gradients are pervasive in the human cortex. While the direction of traveling waves has been linked to brain function and dysfunction, the factors that determine this direction remain elusive. We hypothesized that structural connectivity instrength gradients - defined as the gradually varying sum of incoming connection strengths across the cortex - could shape both traveling wave direction and frequency gradients. We confirm the presence of instrength gradients in the human connectome across diverse cohorts and parcellations. Using a cortical network model, we demonstrate how these instrength gradients direct traveling waves and shape frequency gradients. Our model fits resting-state MEG functional connectivity best in a regime where instrength-directed traveling waves and frequency gradients emerge. We further show how structural subnetworks of the human connectome generate opposing wave directions and frequency gradients observed in the alpha and beta bands. Our findings suggest that structural connectivity instrength gradients affect both traveling wave direction and frequency gradients.

Cytoarchitectonic gradients of laminar degeneration in behavioral variant frontotemporal dementia

Behavioral variant frontotemporal dementia (bvFTD) is a clinical syndrome primarily caused by either tau (bvFTD-tau) or TDP-43 (bvFTD-TDP) proteinopathies. We previously found lower cortical layers and dorsolateral regions accumulate greater tau than TDP-43 pathology; however, patterns of laminar neurodegeneration across diverse cytoarchitecture in bvFTD is understudied. We hypothesized that bvFTD-tau and bvFTD-TDP have distinct laminar distributions of pyramidal neurodegeneration along cortical gradients, a topologic order of cytoarchitectonic subregions based on increasing pyramidal density and laminar differentiation. Here, we tested this hypothesis in a frontal cortical gradient consisting of five cytoarchitectonic types (i.e., periallocortex, agranular mesocortex, dysgranular mesocortex, eulaminate-I isocortex, eulaminate-II isocortex) spanning anterior cingulate, paracingulate, orbitofrontal, and mid-frontal gyri in bvFTD-tau (n=27), bvFTD-TDP (n=47), and healthy controls (HC; n=32). We immunostained all tissue for total neurons (NeuN; neuronal-nuclear protein) and pyramidal neurons (SMI32; non-phosphorylated neurofilament) and digitally quantified NeuN-immunoreactivity (ir) and SMI32-ir in supragranular II-III, infragranular V-VI, and all I-VI layers in each cytoarchitectonic type. We used linear mixed-effects models adjusted for demographic and biologic variables to compare SMI32-ir between groups and examine relationships with the cortical gradient, long-range pathways, and clinical symptoms. We found regional and laminar distributions of SMI32-ir expected for HC, validating our measures within the cortical gradient framework. While SMI32-ir loss was not related to the cortical gradient in bvFTD-TDP, SMI32-ir progressively decreased along the cortical gradient of bvFTD-tau and included greater SMI32-ir loss in supragranular eulaminate-II isocortex in bvFTD-tau vs bvFTD-TDP ( p =0.039). In a structural model for long-range laminar connectivity between infragranular mesocortex and supragranular isocortex, we found a larger laminar ratio of mesocortex-to-isocortex SMI32-ir in bvFTD-tau vs bvFTD-TDP ( p =0.019), suggesting select long-projecting pathways may contribute to isocortical-predominant degeneration in bvFTD-tau. In cytoarchitectonic types with the highest NeuN-ir, we found lower SMI32-ir in bvFTD-tau vs bvFTD-TDP ( p =0.047), suggesting pyramidal neurodegeneration may occur earlier in bvFTD-tau. Lastly, we found that reduced SMI32-ir related to behavioral severity and frontal-mediated letter fluency, not temporal-mediated confrontation naming, demonstrating the clinical relevance and specificity of frontal pyramidal neurodegeneration to bvFTD-related symptoms. Our data suggest loss of neurofilament-rich pyramidal neurons is a clinically relevant feature of bvFTD that selectively worsens along a frontal cortical gradient in bvFTD-tau, not bvFTD-TDP. Therefore, tau-mediated degeneration may preferentially involve pyramidal-rich layers that connect more distant cytoarchitectonic types. Moreover, the hierarchical arrangement of cytoarchitecture along cortical gradients may be an important neuroanatomical framework for identifying which types of cells and pathways are differentially involved between proteinopathies.

Action of GABAB receptor on local network oscillation in somatosensory cortex of oral part: focusing on NMDA receptor

The balance of activity between glutamatergic and GABAergic networks is particularly important for oscillatory neural activities in the brain. Here, we investigated the roles of GABAB receptors in network oscillation in the oral somatosensory cortex (OSC), focusing on NMDA receptors. Neural oscillation at the frequency of 8-10 Hz was elicited in rat brain slices after caffeine application. Oscillations comprised a non-NMDA receptor-dependent initial phase and a later NMDA receptor-dependent oscillatory phase, with the oscillator located in the upper layer of the OSC. Baclofen was applied to investigate the actions of GABAB receptors. The later NMDA receptor-dependent oscillatory phase completely disappeared, but the initial phase did not. These results suggest that GABAB receptors mainly act on NMDA receptor, in which metabotropic actions of GABAB receptors may contribute to the attenuation of NMDA receptor activities. A regulatory system for network oscillation involving GABAB receptors may be present in the OSC.

A ubiquitous spectrolaminar motif of local field potential power across the primate cortex

The mammalian cerebral cortex is anatomically organized into a six-layer motif. It is currently unknown whether a corresponding laminar motif of neuronal activity patterns exists across the cortex. Here we report such a motif in the power of local field potentials (LFPs). Using laminar probes, we recorded LFPs from 14 cortical areas across the cortical hierarchy in five macaque monkeys. The laminar locations of recordings were histologically identified by electrolytic lesions. Across all areas, we found a ubiquitous spectrolaminar pattern characterized by an increasing deep-to-superficial layer gradient of high-frequency power peaking in layers 2/3 and an increasing superficial-to-deep gradient of alpha-beta power peaking in layers 5/6. Laminar recordings from additional species showed that the spectrolaminar pattern is highly preserved among primates-macaque, marmoset and human-but more dissimilar in mouse. Our results suggest the existence of a canonical layer-based and frequency-based mechanism for cortical computation.

Phase of neural oscillations as a reference frame for attention-based routing in visual cortex

Selective attention allows the brain to efficiently process the image projected onto the retina, selectively focusing neural processing resources on behaviorally relevant visual information. While previous studies have documented the crucial role of the action potential rate of single neurons in relaying such information, little is known about how the activity of single neurons relative to their neighboring network contributes to the efficient representation of attended stimuli and transmission of this information to downstream areas. Here, we show in the dorsal visual pathway of monkeys (medial superior temporal area) that neurons fire spikes preferentially at a specific phase of the ongoing population beta (∼20 Hz) oscillations of the surrounding local network. This preferred spiking phase shifts towards a later phase when monkeys selectively attend towards (rather than away from) the receptive field of the neuron. This shift of the locking phase is positively correlated with the speed at which animals report a visual change. Furthermore, our computational modeling suggests that neural networks can manipulate the preferred phase of coupling by imposing differential synaptic delays on postsynaptic potentials. This distinction between the locking phase of neurons activated by the spatially attended stimulus vs. that of neurons activated by the unattended stimulus, may enable the neural system to discriminate relevant from irrelevant sensory inputs and consequently filter out distracting stimuli information by aligning the spikes which convey relevant/irrelevant information to distinct phases linked to periods of better/worse perceptual sensitivity for higher cortices. This strategy may be used to reserve the narrow windows of highest perceptual efficacy to the processing of the most behaviorally relevant information, ensuring highly efficient responses to attended sensory events.

Hierarchical rhythmic propagation of corticothalamic interactions for consciousness: A computational study

Clarifying the mechanisms of loss and recovery of consciousness in the brain is a major challenge in neuroscience, and research on the spatiotemporal organization of rhythms at the brain region scale at different levels of consciousness remains scarce. By applying computational neuroscience, an extended corticothalamic network model was developed in this study to simulate the altered states of consciousness induced by different concentration levels of propofol. The cortex area containing oscillation spread from posterior to anterior in four successive time stages, defining four groups of brain regions. A quantitative analysis showed that hierarchical rhythm propagation was mainly due to heterogeneity in the inter-brain region connections. These results indicate that the proposed model is an anatomically data-driven testbed and a simulation platform with millisecond resolution. It facilitates understanding of activity coordination across multiple areas of the conscious brain and the mechanisms of action of anesthetics in terms of brain regions.

Cell type-specific connectome predicts distributed working memory activity in the mouse brain

Recent advances in connectomics and neurophysiology make it possible to probe whole-brain mechanisms of cognition and behavior. We developed a large-scale model of the multiregional mouse brain for a cardinal cognitive function called working memory, the brain's ability to internally hold and process information without sensory input. The model is built on mesoscopic connectome data for interareal cortical connections and endowed with a macroscopic gradient of measured parvalbumin-expressing interneuron density. We found that working memory coding is distributed yet exhibits modularity; the spatial pattern of mnemonic representation is determined by long-range cell type-specific targeting and density of cell classes. Cell type-specific graph measures predict the activity patterns and a core subnetwork for memory maintenance. The model shows numerous attractor states, which are self-sustained internal states (each engaging a distinct subset of areas). This work provides a framework to interpret large-scale recordings of brain activity during cognition, while highlighting the need for cell type-specific connectomics.

Distinct Patterns of Connectivity between Brain Regions Underlie the Intra-Modal and Cross-Modal Value-Driven Modulations of the Visual Cortex

Past reward associations may be signaled from different sensory modalities; however, it remains unclear how different types of reward-associated stimuli modulate sensory perception. In this human fMRI study (female and male participants), a visual target was simultaneously presented with either an intra- (visual) or a cross-modal (auditory) cue that was previously associated with rewards. We hypothesized that, depending on the sensory modality of the cues, distinct neural mechanisms underlie the value-driven modulation of visual processing. Using a multivariate approach, we confirmed that reward-associated cues enhanced the target representation in early visual areas and identified the brain valuation regions. Then, using an effective connectivity analysis, we tested three possible patterns of connectivity that could underlie the modulation of the visual cortex: a direct pathway from the frontal valuation areas to the visual areas, a mediated pathway through the attention-related areas, and a mediated pathway that additionally involved sensory association areas. We found evidence for the third model demonstrating that the reward-related information in both sensory modalities is communicated across the valuation and attention-related brain regions. Additionally, the superior temporal areas were recruited when reward was cued cross-modally. The strongest dissociation between the intra- and cross-modal reward-driven effects was observed at the level of the feedforward and feedback connections of the visual cortex estimated from the winning model. These results suggest that, in the presence of previously rewarded stimuli from different sensory modalities, a combination of domain-general and domain-specific mechanisms are recruited across the brain to adjust the visual perception.SIGNIFICANCE STATEMENT Reward has a profound effect on perception, but it is not known whether shared or disparate mechanisms underlie the reward-driven effects across sensory modalities. In this human fMRI study, we examined the reward-driven modulation of the visual cortex by visual (intra-modal) and auditory (cross-modal) reward-associated cues. Using a model-based approach to identify the most plausible pattern of inter-regional effective connectivity, we found that higher-order areas involved in the valuation and attentional processing were recruited by both types of rewards. However, the pattern of connectivity between these areas and the early visual cortex was distinct between the intra- and cross-modal rewards. This evidence suggests that, to effectively adapt to the environment, reward signals may recruit both domain-general and domain-specific mechanisms.

Effects of optogenetic and visual stimulation on gamma activity in the visual cortex

Studying brain functions and activity during gamma oscillations can be a challenge because it requires careful planning to create the necessary conditions for a controlled experiment. Such an experiment consists of placing the brain into a gamma state and investigating cognitive processing with a careful design. Cortical oscillations in the gamma frequency range (30-80 Hz) play an essential role in a variety of cognitive processes, including visual processing and cognition. The present study aims to investigate the effects of a visual stimulus on the primary visual cortex under gamma oscillations. Specifically, we sought to explore the behavior of gamma oscillations triggered by optogenetic stimulation in the II and IV layers of the visual cortex, both with and without concurrent visual stimulation. Our results show that optogenetic stimulation increases the power of gamma oscillation in both layers of the visual cortex. However, the combined stimuli resulted in a reduction of gamma power in layer II and an increase and reinforcement in gamma power in layer IV. Modelling the results with the Wilson-Cowan model suggests changes in the input of the excitatory population due to the combined stimuli. In addition, our analysis of the data using the Lempel-Ziv complexity method supports our interpretations from the modeling. Thus, our results suggest that optogenetic stimulation enhances low gamma power in both layers of the visual cortex, while simultaneous visual stimulation has differing effects on the two layers, reducing gamma power in layer II and increasing it in layer IV.

Inter-areal transmission of multiple neural signals through frequency-division-multiplexing communication

Inter-areal information transmission in the brain cortex relates to cognitive functions. Researches used to pay attention to activity pattern transmission, signals gating, or routing in neuronal networks. However, the underlying mechanism of simultaneous transmission of multiple neural signals in the same channel across networks remains unclear. In this work, we construct a two-layer feedforward neuronal network (sender-receiver) with each layer's intrinsic rhythms consisting of slow- (low-frequency) and fast- gamma rhythms (high-frequency), investigating how to realize simultaneous transmission of multiple signals in neuronal systems. With the aid of resonance and frequency analysis, it is shown that low- and high-frequency signals can be transmitted simultaneously in such a feedforward network through frequency division multiplexing (FDM) communication. The transmission performance is related to the local resonance, connectivity, as well as background noise. Moreover, low- and high-frequency signals can also be gated or selected with appropriate adjustments of recurrent connection strength and delay, and background noise. Our model might provide a novel insight into the underlying mechanism of complex signals communication between different cortex areas.

Structural basis of envelope and phase intrinsic coupling modes in the cerebral cortex

Intrinsic coupling modes (ICMs) can be observed in ongoing brain activity at multiple spatial and temporal scales. Two families of ICMs can be distinguished: phase and envelope ICMs. The principles that shape these ICMs remain partly elusive, in particular their relation to the underlying brain structure. Here we explored structure-function relationships in the ferret brain between ICMs quantified from ongoing brain activity recorded with chronically implanted micro-ECoG arrays and structural connectivity (SC) obtained from high-resolution diffusion MRI tractography. Large-scale computational models were used to explore the ability to predict both types of ICMs. Importantly, all investigations were conducted with ICM measures that are sensitive or insensitive to volume conduction effects. The results show that both types of ICMs are significantly related to SC, except for phase ICMs when using measures removing zero-lag coupling. The correlation between SC and ICMs increases with increasing frequency which is accompanied by reduced delays. Computational models produced results that were highly dependent on the specific parameter settings. The most consistent predictions were derived from measures solely based on SC. Overall, the results demonstrate that patterns of cortical functional coupling as reflected in both phase and envelope ICMs are both related, albeit to different degrees, to the underlying structural connectivity in the cerebral cortex.

Spontaneous variations in arousal modulate subsequent visual processing and local field potential dynamics in the ferret during quiet wakefulness

Behavioral states affect neuronal responses throughout the cortex and influence visual processing. Quiet wakefulness (QW) is a behavioral state during which subjects are quiescent but awake and connected to the environment. Here, we examined the effects of pre-stimulus arousal variability on post-stimulus neural activity in the primary visual cortex and posterior parietal cortex in awake ferrets, using pupil diameter as an indicator of arousal. We observed that the power of stimuli-induced alpha (8-12 Hz) decreases when the arousal level increases. The peak of alpha power shifts depending on arousal. High arousal increases inter- and intra-areal coherence. Using a simplified model of laminar circuits, we show that this connectivity pattern is compatible with feedback signals targeting infragranular layers in area posterior parietal cortex and supragranular layers in V1. During high arousal, neurons in V1 displayed higher firing rates at their preferred orientations. Broad-spiking cells in V1 are entrained to high-frequency oscillations (>80 Hz), whereas narrow-spiking neurons are phase-locked to low- (12-18 Hz) and high-frequency (>80 Hz) rhythms. These results indicate that the variability and sensitivity of post-stimulus cortical responses and coherence depend on the pre-stimulus behavioral state and account for the neuronal response variability observed during repeated stimulation.

Noninvasive modulation of predictive coding in humans: causal evidence for frequency-specific temporal dynamics

Increasing evidence indicates that the brain predicts sensory input based on past experiences, importantly constraining how we experience the world. Despite a growing interest on this framework, known as predictive coding, most of such approaches to multiple psychological domains continue to be theoretical or primarily provide correlational evidence. We here explored the neural basis of predictive processing using noninvasive brain stimulation and provide causal evidence of frequency-specific modulations in humans. Participants received 20 Hz (associated with top-down/predictions), 50 Hz (associated with bottom-up/prediction errors), or sham transcranial alternating current stimulation on the left dorsolateral prefrontal cortex while performing a social perception task in which facial expression predictions were induced and subsequently confirmed or violated. Left prefrontal 20 Hz stimulation reinforced stereotypical predictions. In contrast, 50 Hz and sham stimulation failed to yield any significant behavioral effects. Moreover, the frequency-specific effect observed was further supported by electroencephalography data, which showed a boost of brain activity at the stimulated frequency band. These observations provide causal evidence for how predictive processing may be enabled in the human brain, setting up a needed framework to understand how it may be disrupted across brain-related conditions and potentially restored through noninvasive methods.

The fractal brain: scale-invariance in structure and dynamics

The past 40 years have witnessed extensive research on fractal structure and scale-free dynamics in the brain. Although considerable progress has been made, a comprehensive picture has yet to emerge, and needs further linking to a mechanistic account of brain function. Here, we review these concepts, connecting observations across different levels of organization, from both a structural and functional perspective. We argue that, paradoxically, the level of cortical circuits is the least understood from a structural point of view and perhaps the best studied from a dynamical one. We further link observations about scale-freeness and fractality with evidence that the environment provides constraints that may explain the usefulness of fractal structure and scale-free dynamics in the brain. Moreover, we discuss evidence that behavior exhibits scale-free properties, likely emerging from similarly organized brain dynamics, enabling an organism to thrive in an environment that shares the same organizational principles. Finally, we review the sparse evidence for and try to speculate on the functional consequences of fractality and scale-freeness for brain computation. These properties may endow the brain with computational capabilities that transcend current models of neural computation and could hold the key to unraveling how the brain constructs percepts and generates behavior.

Frontal-to-visual information flow explains predictive motion tracking

Predictive tracking demonstrates our ability to maintain a line of vision on moving objects even when they temporarily disappear. Models of smooth pursuit eye movements posit that our brain achieves this ability by directly streamlining motor programming from continuously updated sensory motion information. To test this hypothesis, we obtained sensory motion representation from multivariate electroencephalogram activity while human participants covertly tracked a temporarily occluded moving stimulus with their eyes remaining stationary at the fixation point. The sensory motion representation of the occluded target evolves to its maximum strength at the expected timing of reappearance, suggesting a timely modulation of the internal model of the visual target. We further characterize the spatiotemporal dynamics of the task-relevant motion information by computing the phase gradients of slow oscillations. We discovered a predominant posterior-to-anterior phase gradient immediately after stimulus occlusion; however, at the expected timing of reappearance, the axis reverses the gradient, becoming anterior-to-posterior. The behavioral bias of smooth pursuit eye movements, which is a signature of the predictive process of the pursuit, was correlated with the posterior division of the gradient. These results suggest that the sensory motion area modulated by the prediction signal is involved in updating motor programming.

Functional interactions among neurons within single columns of macaque V1

Recent developments in high-density neurophysiological tools now make it possible to record from hundreds of single neurons within local, highly interconnected neural networks. Among the many advantages of such recordings is that they dramatically increase the quantity of identifiable, functional interactions between neurons thereby providing an unprecedented view of local circuits. Using high-density, Neuropixels recordings from single neocortical columns of primary visual cortex in nonhuman primates, we identified 1000s of functionally interacting neuronal pairs using established crosscorrelation approaches. Our results reveal clear and systematic variations in the synchrony and strength of functional interactions within single cortical columns. Despite neurons residing within the same column, both measures of interactions depended heavily on the vertical distance separating neuronal pairs, as well as on the similarity of stimulus tuning. In addition, we leveraged the statistical power afforded by the large numbers of functionally interacting pairs to categorize interactions between neurons based on their crosscorrelation functions. These analyses identified distinct, putative classes of functional interactions within the full population. These classes of functional interactions were corroborated by their unique distributions across defined laminar compartments and were consistent with known properties of V1 cortical circuitry, such as the lead-lag relationship between simple and complex cells. Our results provide a clear proof-of-principle for the use of high-density neurophysiological recordings to assess circuit-level interactions within local neuronal networks.

Connectivity concepts in neuronal network modeling

Sustainable research on computational models of neuronal networks requires published models to be understandable, reproducible, and extendable. Missing details or ambiguities about mathematical concepts and assumptions, algorithmic implementations, or parameterizations hinder progress. Such flaws are unfortunately frequent and one reason is a lack of readily applicable standards and tools for model description. Our work aims to advance complete and concise descriptions of network connectivity but also to guide the implementation of connection routines in simulation software and neuromorphic hardware systems. We first review models made available by the computational neuroscience community in the repositories ModelDB and Open Source Brain, and investigate the corresponding connectivity structures and their descriptions in both manuscript and code. The review comprises the connectivity of networks with diverse levels of neuroanatomical detail and exposes how connectivity is abstracted in existing description languages and simulator interfaces. We find that a substantial proportion of the published descriptions of connectivity is ambiguous. Based on this review, we derive a set of connectivity concepts for deterministically and probabilistically connected networks and also address networks embedded in metric space. Beside these mathematical and textual guidelines, we propose a unified graphical notation for network diagrams to facilitate an intuitive understanding of network properties. Examples of representative network models demonstrate the practical use of the ideas. We hope that the proposed standardizations will contribute to unambiguous descriptions and reproducible implementations of neuronal network connectivity in computational neuroscience.

Compensatory mechanism of attention-deficit/hyperactivity disorder recovery in resting state alpha rhythms

Alpha rhythms in the human electroencephalogram (EEG), oscillating at 8-13 Hz, are located in parieto-occipital cortex and are strongest when awake people close their eyes. It has been suggested that alpha rhythms were related to attention-related functions and mental disorders (e.g., Attention-deficit/hyperactivity disorder (ADHD)). However, many studies have shown inconsistent results on the difference in alpha oscillation between ADHD and control groups. Hence it is essential to verify this difference. In this study, a dataset of EEG recording (128 channel EGI) from 87 healthy controls (HC) and 162 ADHD (141 persisters and 21 remitters) adults in a resting state with their eyes closed was used to address this question and a three-gauss model (summation of baseline and alpha components) was conducted to fit the data. To our surprise, the power of alpha components was not a significant difference among the three groups. Instead, the baseline power of remission and HC group in the alpha band is significantly stronger than that of persister groups. Our results suggest that ADHD recovery may have compensatory mechanisms and many abnormalities in EEG may be due to the influence of behavior rather than the difference in brain signals.

Rate and oscillatory switching dynamics of a multilayer visual microcircuit model

The neocortex is organized around layered microcircuits consisting of a variety of excitatory and inhibitory neuronal types which perform rate- and oscillation-based computations. Using modeling, we show that both superficial and deep layers of the primary mouse visual cortex implement two ultrasensitive and bistable switches built on mutual inhibitory connectivity motives between somatostatin, parvalbumin, and vasoactive intestinal polypeptide cells. The switches toggle pyramidal neurons between high and low firing rate states that are synchronized across layers through translaminar connectivity. Moreover, inhibited and disinhibited states are characterized by low- and high-frequency oscillations, respectively, with layer-specific differences in frequency and power which show asymmetric changes during state transitions. These findings are consistent with a number of experimental observations and embed firing rate together with oscillatory changes within a switch interpretation of the microcircuit.

Visual evoked feedforward-feedback traveling waves organize neural activity across the cortical hierarchy in mice

Sensory processing is distributed among many brain regions that interact via feedforward and feedback signaling. Neuronal oscillations have been shown to mediate intercortical feedforward and feedback interactions. Yet, the macroscopic structure of the multitude of such oscillations remains unclear. Here, we show that simple visual stimuli reliably evoke two traveling waves with spatial wavelengths that cover much of the cerebral hemisphere in awake mice. 30-50 Hz feedforward waves arise in primary visual cortex (V1) and propagate rostrally, while 3-6 Hz feedback waves originate in the association cortex and flow caudally. The phase of the feedback wave modulates the amplitude of the feedforward wave and synchronizes firing between V1 and parietal cortex. Altogether, these results provide direct experimental evidence that visual evoked traveling waves percolate through the cerebral cortex and coordinate neuronal activity across broadly distributed networks mediating visual processing.

The quest for multiscale brain modeling

Addressing the multiscale organization of the brain, which is fundamental to the dynamic repertoire of the organ, remains challenging. In principle, it should be possible to model neurons and synapses in detail and then connect them into large neuronal assemblies to explain the relationship between microscopic phenomena, large-scale brain functions, and behavior. It is more difficult to infer neuronal functions from ensemble measurements such as those currently obtained with brain activity recordings. In this article we consider theories and strategies for combining bottom-up models, generated from principles of neuronal biophysics, with top-down models based on ensemble representations of network activity and on functional principles. These integrative approaches are hoped to provide effective multiscale simulations in virtual brains and neurorobots, and pave the way to future applications in medicine and information technologies.

Thalamic deep brain stimulation paradigm to reduce consciousness: Cortico-striatal dynamics implicated in mechanisms of consciousness

Anesthetic manipulations provide much-needed causal evidence for neural correlates of consciousness, but non-specific drug effects complicate their interpretation. Evidence suggests that thalamic deep brain stimulation (DBS) can either increase or decrease consciousness, depending on the stimulation target and parameters. The putative role of the central lateral thalamus (CL) in consciousness makes it an ideal DBS target to manipulate circuit-level mechanisms in cortico-striato-thalamic (CST) systems, thereby influencing consciousness and related processes. We used multi-microelectrode DBS targeted to CL in macaques while recording from frontal, parietal, and striatal regions. DBS induced episodes of abnormally long, vacant staring with low-frequency oscillations here termed vacant, perturbed consciousness (VPC). DBS modulated VPC likelihood in a frequency-specific manner. VPC events corresponded to decreases in measures of neural complexity (entropy) and integration (Φ*), proposed indices of consciousness, and substantial changes to communication in CST circuits. During VPC, power spectral density and coherence at low frequencies increased across CST circuits, especially in thalamo-parietal and cortico-striatal pathways. Decreased consciousness and neural integration corresponded to shifts in cortico-striatal network configurations that dissociated parietal and subcortical structures. Overall, the features of VPC and implicated networks were similar to those of absence epilepsy. As this same multi-microelectrode DBS method–but at different stimulation frequencies–can also increase consciousness in anesthetized macaques, it can be used to flexibly address questions of consciousness with limited confounds, as well as inform clinical investigations of other consciousness disorders.

Deep brain stimulation (DBS) is an effective treatment for a number of neurological and psychiatric disorders. Recent evidence shows that thalamic DBS can increase consciousness under general anesthesia, and has the potential to be an effective treatment for disorders of consciousness like coma. However, past studies also imply that thalamic DBS can decrease consciousness by producing periods of inactivity with vacant staring. In this study, we use thalamic DBS to induce vacant periods of perturbed consciousness (VPC) while recording from frontal, parietal and subcortical brain areas. Much of modern research on consciousness focuses on cortical contributions, debating if frontal or more posterior areas are more important. We show via machine learning, information theoretic and functional connectivity measures that VPC involves altered neural communication between parietal cortex and subcortical regions, challenging a cortico-centric focus. Further, these mechanisms resemble those associated with absence epilepsy, suggesting a link between this disorder of consciousness and VPC. Finally, occurrence of VPC in our study was stimulation frequency-dependent and influenced by the experimental history of DBS. This provides insight to mechanisms of thalamic DBS, informative for clinical manipulations of consciousness.

Theory of the Multiregional Neocortex: Large-Scale Neural Dynamics and Distributed Cognition

The neocortex is a complex neurobiological system with many interacting regions. How these regions work together to subserve flexible behavior and cognition has become increasingly amenable to rigorous research. Here, I review recent experimental and theoretical work on the modus operandi of a multiregional cortex. These studies revealed several general principles for the neocortical interareal connectivity, low-dimensional macroscopic gradients of biological properties across cortical areas, and a hierarchy of timescales for information processing. Theoretical work suggests testable predictions regarding differential excitation and inhibition along feedforward and feedback pathways in the cortical hierarchy. Furthermore, modeling of distributed working memory and simple decision-making has given rise to a novel mathematical concept, dubbed bifurcation in space, that potentially explains how different cortical areas, with a canonical circuit organization but gradients of biological heterogeneities, are able to subserve their respective (e.g., sensory coding versus executive control) functions in a modularly organized brain.

Normalizing the brain connectome for communication through synchronization

Networks in neuroscience determine how brain function unfolds, and their perturbations lead to psychiatric disorders and brain disease. Brain networks are characterized by their connectomes, which comprise the totality of all connections, and are commonly described by graph theory. This approach is deeply rooted in a particle view of information processing, based on the quantification of informational bits such as firing rates. Oscillations and brain rhythms demand, however, a wave perspective of information processing based on synchronization. We extend traditional graph theory to a dual, particle-wave, perspective, integrate time delays due to finite transmission speeds, and derive a normalization of the connectome. When applied to the database of the Human Connectome Project, it explains the emergence of frequency-specific network cores including the visual and default mode networks. These findings are robust across human subjects (N = 100) and are a fundamental network property within the wave picture. The normalized connectome comprises the particle view in the limit of infinite transmission speeds and opens the applicability of graph theory to a wide range of novel network phenomena, including physiological and pathological brain rhythms. These two perspectives are orthogonal, but not incommensurable, when understood within the novel, here-proposed, generalized framework of structural connectivity.

Bipolar reasoning in feedback pathways

Instead of the conventional 0 and 1 values, bipolar reasoning uses -1, 0, +1 to describe double-sided judgements in which neutral elements are halfway between positive and negative evaluations (e.g., "uncertain" lies between "impossible" and "totally sure"). We discuss the state-of-the-art in bipolar logics and recall two medieval forerunners, i.e., William of Ockham and Nicholas of Autrecourt, who embodied a bipolar mode of thought that is eminently modern. Starting from the trivial observation that "once a wheat sheaf is sealed and tied up, the packed down straws display the same orientation", we work up a new theory of the bipolar nature of networks, suggesting that orthodromic (i.e., feedforward, bottom-up) projections might be functionally coupled with antidromic (i.e., feedback, top-down) projections via the mathematical apparatus of presheaves/globular sets. When an entrained oscillation such as a neuronal spike propagates from A to B, changes in B might lead to changes in A, providing unexpected antidromic effects. Our account points towards the methodological feasibility of novel neural networks in which message feedback is guaranteed by backpropagation mechanisms endowed in the same feedforward circuits. Bottom-up/top-down transmission at various coarse-grained network levels provides fresh insights in far-flung scientific fields such as object persistence, memory reinforcement, visual recognition, Bayesian inferential circuits and multidimensional activity of the brain. Implying that axonal stimulation by external sources might backpropagate and modify neuronal electric oscillations, our theory also suggests testable previsions concerning the optimal location of transcranial magnetic stimulation's coils in patients affected by drug-resistant epilepsy.

Extracting Dynamical Understanding From Neural-Mass Models of Mouse Cortex

New brain atlases with high spatial resolution and whole-brain coverage have rapidly advanced our knowledge of the brain's neural architecture, including the systematic variation of excitatory and inhibitory cell densities across the mammalian cortex. But understanding how the brain's microscale physiology shapes brain dynamics at the macroscale has remained a challenge. While physiologically based mathematical models of brain dynamics are well placed to bridge this explanatory gap, their complexity can form a barrier to providing clear mechanistic interpretation of the dynamics they generate. In this work, we develop a neural-mass model of the mouse cortex and show how bifurcation diagrams, which capture local dynamical responses to inputs and their variation across brain regions, can be used to understand the resulting whole-brain dynamics. We show that strong fits to resting-state functional magnetic resonance imaging (fMRI) data can be found in surprisingly simple dynamical regimes—including where all brain regions are confined to a stable fixed point—in which regions are able to respond strongly to variations in their inputs, consistent with direct structural connections providing a strong constraint on functional connectivity in the anesthetized mouse. We also use bifurcation diagrams to show how perturbations to local excitatory and inhibitory coupling strengths across the cortex, constrained by cell-density data, provide spatially dependent constraints on resulting cortical activity, and support a greater diversity of coincident dynamical regimes. Our work illustrates methods for visualizing and interpreting model performance in terms of underlying dynamical mechanisms, an approach that is crucial for building explanatory and physiologically grounded models of the dynamical principles that underpin large-scale brain activity.

Spatially distributed computation in cortical circuits

The traditional view of neural computation in the cerebral cortex holds that sensory neurons are specialized, i.e., selective for certain dimensions of sensory stimuli. This view was challenged by evidence of contextual interactions between stimulus dimensions in which a neuron's response to one dimension strongly depends on other dimensions. Here, we use methods of mathematical modeling, psychophysics, and electrophysiology to address shortcomings of the traditional view. Using a model of a generic cortical circuit, we begin with the simple demonstration that cortical responses are always distributed among neurons, forming characteristic waveforms, which we call neural waves. When stimulated by patterned stimuli, circuit responses arise by interference of neural waves. Results of this process depend on interaction between stimulus dimensions. Comparison of modeled responses with responses of biological vision makes it clear that the framework of neural wave interference provides a useful alternative to the standard concept of neural computation.

Corticocortical innervation subtypes of layer 5 intratelencephalic cells in the murine secondary motor cortex

Feedback projections from the secondary motor cortex (M2) to the primary motor and sensory cortices are essential for behavior selection and sensory perception. Intratelencephalic (IT) cells in layer 5 (L5) contribute feedback projections to diverse cortical areas. Here we show that L5 IT cells participating in feedback connections to layer 1 (L1) exhibit distinct projection patterns, genetic profiles, and electrophysiological properties relative to other L5 IT cells. An analysis of the MouseLight database found that L5 IT cells preferentially targeting L1 project broadly to more cortical regions, including the perirhinal and auditory cortices, and innervate a larger volume of striatum than the other L5 IT cells. We found experimentally that in upper L5 (L5a), ER81 (ETV1) was found more often in L1-preferring IT cells, and in IT cells projecting to perirhinal/auditory regions than those projecting to primary motor or somatosensory regions. The perirhinal region-projecting L5a IT cells were synaptically connected to each other and displayed lower input resistance than contra-M2 projecting IT cells including L1-preferring and nonpreferring cells. Our findings suggest that M2-L5a IT L1-preferring cells exhibit stronger ER81 expression and broader cortical/striatal projection fields than do cells that do not preferentially target L1.

Mechanisms of distributed working memory in a large-scale network of macaque neocortex

Neural activity underlying working memory is not a local phenomenon but distributed across multiple brain regions. To elucidate the circuit mechanism of such distributed activity, we developed an anatomically constrained computational model of large-scale macaque cortex. We found that mnemonic internal states may emerge from inter-areal reverberation, even in a regime where none of the isolated areas is capable of generating self-sustained activity. The mnemonic activity pattern along the cortical hierarchy indicates a transition in space, separating areas engaged in working memory and those which do not. A host of spatially distinct attractor states is found, potentially subserving various internal processes. The model yields testable predictions, including the idea of counterstream inhibitory bias, the role of prefrontal areas in controlling distributed attractors, and the resilience of distributed activity to lesions or inactivation. This work provides a theoretical framework for identifying large-scale brain mechanisms and computational principles of distributed cognitive processes.

Feedforward and feedback influences through distinct frequency bands between two spiking-neuron networks

Several studies on brain signals suggested that bottom-up and top-down influences are exerted through distinct frequency bands among visual cortical areas. It was recently shown that theta and gamma rhythms subserve feedforward, whereas the feedback influence is dominated by the alpha-beta rhythm in primates. A few theoretical models for reproducing these effects have been proposed so far. Here we show that a simple but biophysically plausible two-network motif composed of spiking-neuron models and chemical synapses can exhibit feedforward and feedback influences through distinct frequency bands. Different from previous studies, this kind of model allows us to study directed influences not only at the population level, by using a proxy for the local field potential, but also at the cellular level, by using the neuronal spiking series.

Gamma rhythms in the visual cortex: functions and mechanisms

Gamma-band activity, peaking around 30–100 Hz in the local field potential's power spectrum, has been found and intensively studied in many brain regions. Although gamma is thought to play a critical role in processing neural information in the brain, its cognitive functions and neural mechanisms remain unclear or debatable. Experimental studies showed that gamma rhythms are stochastic in time and vary with visual stimuli. Recent studies further showed that multiple rhythms coexist in V1 with distinct origins in different species. While all these experimental facts are a challenge for understanding the functions of gamma in the visual cortex, there are many signs of progress in computational studies. This review summarizes and discusses studies on gamma in the visual cortex from multiple perspectives and concludes that gamma rhythms are still a mystery. Combining experimental and computational studies seems the best way forward in the future.

Predisposition to domain-wide maladaptive changes in predictive coding in auditory phantom perception

Tinnitus is hypothesised to be a predictive coding problem. Previous research indicates lower sensitivity to prediction errors (PEs) in tinnitus patients while processing auditory deviants corresponding to tinnitus-specific stimuli. However, based on research with patients with hallucinations and no psychosis we hypothesise tinnitus patients may be more sensitive to PEs produced by auditory stimuli that are not related to tinnitus characteristics. Specifically in patients with minimal to no hearing loss, we hypothesise a more top-down subtype of tinnitus that may be driven by maladaptive changes in an auditory predictive coding network. To test this, we use an auditory oddball paradigm with omission of global and local deviants, a measure that is previously shown to empirically characterise hierarchical prediction errors (PEs). We observe: (1) increased predictions characterised by increased pre-stimulus response and increased alpha connectivity between the parahippocampus, dorsal anterior cingulate cortex and parahippocampus, pregenual anterior cingulate cortex and posterior cingulate cortex; (2) increased PEs characterised by increased P300 amplitude and gamma activity and increased theta connectivity between auditory cortices, parahippocampus and dorsal anterior cingulate cortex in the tinnitus group; (3) increased overall feed-forward connectivity in theta from the auditory cortex and parahippocampus to the dorsal anterior cingulate cortex; (4) correlations of pre-stimulus theta activity to tinnitus loudness and alpha activity to tinnitus distress. These results provide empirical evidence of maladaptive changes in a hierarchical predictive coding network in a subgroup of tinnitus patients with minimal to no hearing loss. The changes in pre-stimulus activity and connectivity to non-tinnitus specific stimuli suggest that tinnitus patients not only produce strong predictions about upcoming stimuli but also may be predisposed to stimulus a-specific PEs in the auditory domain. Correlations with tinnitus-related characteristics may be a biomarker for maladaptive changes in auditory predictive coding.

A mechanism for inter-areal coherence through communication based on connectivity and oscillatory power

Inter-areal coherence between field potentials is a widespread phenomenon in cortex. Coherence has been hypothesized to reflect phase-synchronization between oscillators and flexibly gate communication according to behavioral and cognitive demands. We reveal an alternative mechanism where coherence is not the cause but the consequence of communication and naturally emerges because spiking activity in a sending area causes post-synaptic potentials both in the same and in other areas. Consequently, coherence depends in a lawful manner on power and phase-locking in the sender and connectivity. Changes in oscillatory power explained prominent changes in fronto-parietal and LGN-V1 coherence across behavioral conditions. Optogenetic experiments and excitatory-inhibitory network simulations identified afferent synaptic inputs rather than spiking entrainment as the principal determinant of coherence. These findings suggest that unique spectral profiles of different brain areas automatically give rise to large-scale coherence patterns that follow anatomical connectivity and continuously reconfigure as a function of behavior and cognition.

Multiple gamma rhythms carry distinct spatial frequency information in primary visual cortex

Gamma rhythms in many brain regions, including the primary visual cortex (V1), are thought to play a role in information processing. Here, we report a surprising finding of 3 narrowband gamma rhythms in V1 that processed distinct spatial frequency (SF) signals and had different neural origins. The low gamma (LG; 25 to 40 Hz) rhythm was generated at the V1 superficial layer and preferred a higher SF compared with spike activity, whereas both the medium gamma (MG; 40 to 65 Hz), generated at the cortical level, and the high gamma HG; (65 to 85 Hz), originated precortically, preferred lower SF information. Furthermore, compared with the rates of spike activity, the powers of the 3 gammas had better performance in discriminating the edge and surface of simple objects. These findings suggest that gamma rhythms reflect the neural dynamics of neural circuitries that process different SF information in the visual system, which may be crucial for multiplexing SF information and synchronizing different features of an object.

A dopamine gradient controls access to distributed working memory in the large-scale monkey cortex

Dopamine is required for working memory, but how it modulates the large-scale cortex is unknown. Here, we report that dopamine receptor density per neuron, measured by autoradiography, displays a macroscopic gradient along the macaque cortical hierarchy. This gradient is incorporated in a connectome-based large-scale cortex model endowed with multiple neuron types. The model captures an inverted U-shaped dependence of working memory on dopamine and spatial patterns of persistent activity observed in over 90 experimental studies. Moreover, we show that dopamine is crucial for filtering out irrelevant stimuli by enhancing inhibition from dendrite-targeting interneurons. Our model revealed that an activity-silent memory trace can be realized by facilitation of inter-areal connections and that adjusting cortical dopamine induces a switch from this internal memory state to distributed persistent activity. Our work represents a cross-level understanding from molecules and cell types to recurrent circuit dynamics underlying a core cognitive function distributed across the primate cortex.

Structural Attributes and Principles of the Neocortical Connectome in the Marmoset Monkey

The marmoset monkey has become an important primate model in Neuroscience. Here, we characterize salient statistical properties of interareal connections of the marmoset cerebral cortex, using data from retrograde tracer injections. We found that the connectivity weights are highly heterogeneous, spanning 5 orders of magnitude, and are log-normally distributed. The cortico-cortical network is dense, heterogeneous and has high specificity. The reciprocal connections are the most prominent and the probability of connection between 2 areas decays with their functional dissimilarity. The laminar dependence of connections defines a hierarchical network correlated with microstructural properties of each area. The marmoset connectome reveals parallel streams associated with different sensory systems. Finally, the connectome is spatially embedded with a characteristic length that obeys a power law as a function of brain volume across rodent and primate species. These findings provide a connectomic basis for investigations of multiple interacting areas in a complex large-scale cortical system underlying cognitive processes.

Cortical interaction of bilateral inputs is similar for noxious and innocuous stimuli but leads to different perceptual effects

The cerebral integration of somatosensory inputs from multiple sources is essential to produce adapted behaviors. Previous studies suggest that bilateral somatosensory inputs interact differently depending on stimulus characteristics, including their noxious nature. The aim of this study was to clarify how bilateral inputs evoked by noxious laser stimuli, noxious shocks, and innocuous shocks interact in terms of perception and brain responses. The experiment comprised two conditions (right-hand stimulation and concurrent stimulation of both hands) in which painful laser stimuli, painful shocks and non-painful shocks were delivered. Perception, somatosensory-evoked potentials (P45, N100, P260), laser-evoked potentials (N1, N2 and P2) and event-related spectral perturbations (delta to gamma oscillation power) were compared between conditions and stimulus modalities. The amplitude of negative vertex potentials (N2 or N100) and the power of delta/theta oscillations were increased in the bilateral compared with unilateral condition, regardless of the stimulus type (P < 0.01). However, gamma oscillation power increased for painful and non-painful shocks (P < 0.01), but not for painful laser stimuli (P = 0.08). Despite the similarities in terms of brain activity, bilateral inputs interacted differently for painful stimuli, for which perception remained unchanged, and non-painful stimuli, for which perception increased. This may reflect a ceiling effect for the attentional capture by noxious stimuli and warrants further investigations to examine the regulation of such interactions by bottom-up and top-down processes.

Network Path Convergence Shapes Low-Level Processing in the Visual Cortex

Hierarchical counterstream via feedforward and feedback interactions is a major organizing principle of the cerebral cortex. The counterstream, as a topological feature of the network of cortical areas, is captured by the convergence and divergence of paths through directed links. So defined, the convergence degree (CD) reveals the reciprocal nature of forward and backward connections, and also hierarchically relevant integrative properties of areas through their inward and outward connections. We asked if topology shapes large-scale cortical functioning by studying the role of CD in network resilience and Granger causal coupling in a model of hierarchical network dynamics. Our results indicate that topological synchronizability is highly vulnerable to attacking edges based on CD, while global network efficiency depends mostly on edge betweenness, a measure of the connectedness of a link. Furthermore, similar to anatomical hierarchy determined by the laminar distribution of connections, CD highly correlated with causal coupling in feedforward gamma, and feedback alpha-beta band synchronizations in a well-studied subnetwork, including low-level visual cortical areas. In contrast, causal coupling did not correlate with edge betweenness. Considering the entire network, the CD-based hierarchy correlated well with both the anatomical and functional hierarchy for low-level areas that are far apart in the hierarchy. Conversely, in a large part of the anatomical network where hierarchical distances are small between the areas, the correlations were not significant. These findings suggest that CD-based and functional hierarchies are interrelated in low-level processing in the visual cortex. Our results are consistent with the idea that the interplay of multiple hierarchical features forms the basis of flexible functional cortical interactions.

A Complex-Valued Oscillatory Neural Network for Storage and Retrieval of Multidimensional Aperiodic Signals

Recurrent neural networks with associative memory properties are typically based on fixed-point dynamics, which is fundamentally distinct from the oscillatory dynamics of the brain. There have been proposals for oscillatory associative memories, but here too, in the majority of cases, only binary patterns are stored as oscillatory states in the network. Oscillatory neural network models typically operate at a single/common frequency. At multiple frequencies, even a pair of oscillators with real coupling exhibits rich dynamics of Arnold tongues, not easily harnessed to achieve reliable memory storage and retrieval. Since real brain dynamics comprises of a wide range of spectral components, there is a need for oscillatory neural network models that operate at multiple frequencies. We propose an oscillatory neural network that can model multiple time series simultaneously by performing a Fourier-like decomposition of the signals. We show that these enhanced properties of a network of Hopf oscillators become possible by operating in the complex-variable domain. In this model, the single neural oscillator is modeled as a Hopf oscillator, with adaptive frequency and dynamics described over the complex domain. We propose a novel form of coupling, dubbed "power coupling," between complex Hopf oscillators. With power coupling, expressed naturally only in the complex-variable domain, it is possible to achieve stable (normalized) phase relationships in a network of multifrequency oscillators. Network connections are trained either by Hebb-like learning or by delta rule, adapted to the complex domain. The network is capable of modeling N-channel electroencephalogram time series with high accuracy and shows the potential as an effective model of large-scale brain dynamics.

Cerebellar Purkinje cells can differentially modulate coherence between sensory and motor cortex depending on region and behavior

Activity of sensory and motor cortices is essential for sensorimotor integration. In particular, coherence between these areas may indicate binding of critical functions like perception, motor planning, action, or sleep. Evidence is accumulating that cerebellar output modulates cortical activity and coherence, but how, when, and where it does so is unclear. We studied activity in and coherence between S1 and M1 cortices during whisker stimulation in the absence and presence of optogenetic Purkinje cell stimulation in crus 1 and 2 of awake mice, eliciting strong simple spike rate modulation. Without Purkinje cell stimulation, whisker stimulation triggers fast responses in S1 and M1 involving transient coherence in a broad spectrum. Simultaneous stimulation of Purkinje cells and whiskers affects amplitude and kinetics of sensory responses in S1 and M1 and alters the estimated S1-M1 coherence in theta and gamma bands, allowing bidirectional control dependent on behavioral context. These effects are absent when Purkinje cell activation is delayed by 20 ms. Focal stimulation of Purkinje cells revealed site specificity, with cells in medial crus 2 showing the most prominent and selective impact on estimated coherence, i.e., a strong suppression in the gamma but not the theta band. Granger causality analyses and computational modeling of the involved networks suggest that Purkinje cells control S1-M1 phase consistency predominantly via ventrolateral thalamus and M1. Our results indicate that activity of sensorimotor cortices can be dynamically and functionally modulated by specific cerebellar inputs, highlighting a widespread role of the cerebellum in coordinating sensorimotor behavior.

Computational exploration of dynamic mechanisms of steady state visual evoked potentials at the whole brain level

Periodic visual stimulation can induce stable steady-state visual evoked potentials (SSVEPs) distributed in multiple brain regions and has potential applications in both neural engineering and cognitive neuroscience. However, the underlying dynamic mechanisms of SSVEPs at the whole-brain level are still not completely understood. Here, we addressed this issue by simulating the rich dynamics of SSVEPs with a large-scale brain model designed with constraints of neuroimaging data acquired from the human brain. By eliciting activity of the occipital areas using an external periodic stimulus, our model was capable of replicating both the spatial distributions and response features of SSVEPs that were observed in experiments. In particular, we confirmed that alpha-band (8-12 Hz) stimulation could evoke stronger SSVEP responses; this frequency sensitivity was due to nonlinear entrainment and resonance, and could be modulated by endogenous factors in the brain. Interestingly, the stimulus-evoked brain networks also exhibited significant superiority in topological properties near this frequency-sensitivity range, and stronger SSVEP responses were demonstrated to be supported by more efficient functional connectivity at the neural activity level. These findings not only provide insights into the mechanistic understanding of SSVEPs at the whole-brain level but also indicate a bright future for large-scale brain modeling in characterizing the complicated dynamics and functions of the brain.

Consciousness depends on integration between parietal cortex, striatum, and thalamus

The neural substrates of consciousness remain elusive. Competing theories that attempt to explain consciousness disagree on the contribution of frontal versus posterior cortex and omit subcortical influences. This lack of understanding impedes the ability to monitor consciousness, which can lead to adverse clinical consequences. To test substrates and measures of consciousness, we recorded simultaneously from frontal cortex, parietal cortex, and subcortical structures, the striatum and thalamus, in awake, sleeping, and anesthetized macaques. We manipulated consciousness on a finer scale using thalamic stimulation, rousing macaques from continuously administered anesthesia. Our results show that, unlike measures targeting complexity, a measure additionally capturing neural integration (Φ∗) robustly correlated with changes in consciousness. Machine learning approaches show parietal cortex, striatum, and thalamus contributed more than frontal cortex to decoding differences in consciousness. These findings highlight the importance of integration between parietal and subcortical structures and challenge a key role for frontal cortex in consciousness.

Excitatory-inhibitory tone shapes decision strategies in a hierarchical neural network model of multi-attribute choice

We are constantly faced with decisions between alternatives defined by multiple attributes, necessitating an evaluation and integration of different information sources. Time-varying signals in multiple brain areas are implicated in decision-making; but we lack a rigorous biophysical description of how basic circuit properties, such as excitatory-inhibitory (E/I) tone and cascading nonlinearities, shape attribute processing and choice behavior. Furthermore, how such properties govern choice performance under varying levels of environmental uncertainty is unknown. We investigated two-attribute, two-alternative decision-making in a dynamical, cascading nonlinear neural network with three layers: an input layer encoding choice alternative attribute values; an intermediate layer of modules processing separate attributes; and a final layer producing the decision. Depending on intermediate layer E/I tone, the network displays distinct regimes characterized by linear (I), convex (II) or concave (III) choice indifference curves. In regimes I and II, each option's attribute information is additively integrated. In regime III, time-varying nonlinear operations amplify the separation between offer distributions by selectively attending to the attribute with the larger differences in input values. At low environmental uncertainty, a linear combination most consistently selects higher valued alternatives. However, at high environmental uncertainty, regime III is more likely than a linear operation to select alternatives with higher value. Furthermore, there are conditions where readout from the intermediate layer could be experimentally indistinguishable from the final layer. Finally, these principles are used to examine multi-attribute decisions in systems with reduced inhibitory tone, leading to predictions of different choice patterns and overall performance between those with restrictions on inhibitory tone and neurotypicals.