Low-dimensional criticality embedded in high-dimensional awake brain dynamics.

Whether cortical neurons operate in a strongly or weakly correlated dynamical regime determines fundamental information processing capabilities and has fueled decades of debate. We offer a resolution of this debate; we show that two important dynamical regimes, typically considered incompatible, can coexist in the same local cortical circuit by separating them into two different subspaces. In awake mouse motor cortex, we find a low-dimensional subspace with large fluctuations consistent with criticality-a dynamical regime with moderate correlations and multi-scale information capacity and transmission. Orthogonal to this critical subspace, we find a high-dimensional subspace containing a desynchronized dynamical regime, which may optimize input discrimination. The critical subspace is apparent only at long timescales, which explains discrepancies among some previous studies. Using a computational model, we show that the emergence of a low-dimensional critical subspace at large timescales agrees with established theory of critical dynamics. Our results suggest that the cortex leverages its high dimensionality to multiplex dynamical regimes across different subspaces.

Sensory input to cortex encoded on low-dimensional periphery-correlated subspaces.

As information about the world is conveyed from the sensory periphery to central neural circuits, it mixes with complex ongoing cortical activity. How do neural populations keep track of sensory signals, separating them from noisy ongoing activity? Here, we show that sensory signals are encoded more reliably in certain low-dimensional subspaces. These coding subspaces are defined by correlations between neural activity in the primary sensory cortex and upstream sensory brain regions; the most correlated dimensions were best for decoding. We analytically show that these correlation-based coding subspaces improve, reaching optimal limits (without an ideal observer), as noise correlations between cortex and upstream regions are reduced. We show that this principle generalizes across diverse sensory stimuli in the olfactory system and the visual system of awake mice. Our results demonstrate an algorithm the cortex may use to multiplex different functions, processing sensory input in low-dimensional subspaces separate from other ongoing functions.

Odor modality is transmitted to cortical brain regions from the olfactory bulb.

Odor perception is the impetus for important animal behaviors with two predominate modes of processing: odors pass through the front of the nose (orthonasal) while inhaling and sniffing, or through the rear (retronasal) during exhalation and while eating. Despite the importance of olfaction for an animal's well-being and that ortho and retro naturally occur, it is unknown how the modality (ortho vs. retro) is even transmitted to cortical brain regions, which could significantly affect how odors are processed and perceived. Using multielectrode array recordings in tracheotomized anesthetized rats, which decouples ortho-retro modality from breathing, we show that mitral cells in rat olfactory bulb can reliably and directly transmit orthonasal versus retronasal modality with ethyl butyrate, a common food odor. Drug manipulations affecting synaptic inhibition via GABAA lead to worse decoding of ortho versus retro, independent of whether overall inhibition increases or decreases, suggesting that the olfactory bulb circuit may naturally favor encoding this important aspect of odors. Detailed data analysis paired with a firing rate model that captures population trends in spiking statistics shows how this circuit can encode odor modality. We have not only demonstrated that ortho/retro information is encoded to downstream brain regions but also used modeling to demonstrate a plausible mechanism for this encoding; due to synaptic adaptation, it is the slower time course of the retronasal stimulation that causes retronasal responses to be stronger and less sensitive to inhibitory drug manipulations than orthonasal responses.NEW & NOTEWORTHY Whether ortho (sniffing odors) versus retro (exhalation and eating) is encoded from the olfactory bulb to other brain areas is not completely known. Using multielectrode array recordings in anesthetized rats, we show that the olfactory bulb transmits this information downstream via spikes. Altering inhibition degrades ortho/retro information on average. We use theory and computation to explain our results, which should have implications on cortical processing considering that only food odors occur retronasally.

Low dimensional criticality embedded in high dimensional awake brain dynamics.

Whether cortical neurons operate in a strongly or weakly correlated dynamical regime determines fundamental information processing capabilities and has fueled decades of debate. Here we offer a resolution of this debate; we show that two important dynamical regimes, typically considered incompatible, can coexist in the same local cortical circuit by separating them into two different subspaces. In awake mouse motor cortex, we find a low-dimensional subspace with large fluctuations consistent with criticality - a dynamical regime with moderate correlations and multi-scale information capacity and transmission. Orthogonal to this critical subspace, we find a high-dimensional subspace containing a desynchronized dynamical regime, which may optimize input discrimination. The critical subspace is apparent only at long timescales, which explains discrepancies among some previous studies. Using a computational model, we show that the emergence of a low-dimensional critical subspace at large timescale agrees with established theory of critical dynamics. Our results suggest that cortex leverages its high dimensionality to multiplex dynamical regimes across different subspaces.

Collapse of complexity of brain and body activity due to excessive inhibition and MeCP2 disruption.

Complex body movements require complex dynamics and coordination among neurons in motor cortex. Conversely, a long-standing theoretical notion supposes that if many neurons in motor cortex become excessively synchronized, they may lack the necessary complexity for healthy motor coding. However, direct experimental support for this idea is rare and underlying mechanisms are unclear. Here we recorded three-dimensional body movements and spiking activity of many single neurons in motor cortex of rats with enhanced synaptic inhibition and a transgenic rat model of Rett syndrome (RTT). For both cases, we found a collapse of complexity in the motor system. Reduced complexity was apparent in lower-dimensional, stereotyped brain-body interactions, neural synchrony, and simpler behavior. Our results demonstrate how imbalanced inhibition can cause excessive synchrony among movement-related neurons and, consequently, a stereotyped motor code. Excessive inhibition and synchrony may underlie abnormal motor function in RTT.

Differences in olfactory bulb mitral cell spiking with ortho- and retronasal stimulation revealed by data-driven models.

The majority of olfaction studies focus on orthonasal stimulation where odors enter via the front nasal cavity, while retronasal olfaction, where odors enter the rear of the nasal cavity during feeding, is understudied. The coding of retronasal odors via coordinated spiking of neurons in the olfactory bulb (OB) is largely unknown despite evidence that higher level processing is different than orthonasal. To this end, we use multi-electrode array in vivo recordings of rat OB mitral cells (MC) in response to a food odor with both modes of stimulation, and find significant differences in evoked firing rates and spike count covariances (i.e., noise correlations). Differences in spiking activity often have implications for sensory coding, thus we develop a single-compartment biophysical OB model that is able to reproduce key properties of important OB cell types. Prior experiments in olfactory receptor neurons (ORN) showed retro stimulation yields slower and spatially smaller ORN inputs than with ortho, yet whether this is consequential for OB activity remains unknown. Indeed with these specifications for ORN inputs, our OB model captures the salient trends in our OB data. We also analyze how first and second order ORN input statistics dynamically transfer to MC spiking statistics with a phenomenological linear-nonlinear filter model, and find that retro inputs result in larger linear filters than ortho inputs. Finally, our models show that the temporal profile of ORN is crucial for capturing our data and is thus a distinguishing feature between ortho and retro stimulation, even at the OB. Using data-driven modeling, we detail how ORN inputs result in differences in OB dynamics and MC spiking statistics. These differences may ultimately shape how ortho and retro odors are coded.

Odor-evoked increases in olfactory bulb mitral cell spiking variability.

The spiking variability of neural networks has important implications for how information is encoded to higher brain regions. It has been well documented by numerous labs in many cortical and motor regions that spiking variability decreases with stimulus onset, yet whether this principle holds in the OB has not been tested. In stark contrast to this common view, we demonstrate that the onset of sensory input can cause an increase in the variability of neural activity in the mammalian OB. We show this in both anesthetized and awake rodents. Furthermore, we use computational models to describe the mechanisms of this phenomenon. Our findings establish sensory evoked increases in spiking variability as a viable alternative coding strategy.

Probing spatial inhomogeneity of cholinergic changes in cortical state in rat.

Acetylcholine (ACh) plays an essential role in cortical information processing. Cholinergic changes in cortical state can fundamentally change how the neurons encode sensory input and motor output. Traditionally, ACh distribution in cortex and associated changes in cortical state have been assumed to be spatially diffuse. However, recent studies demonstrate a more spatially inhomogeneous structure of cholinergic projections to cortex. Moreover, many experimental manipulations of ACh have been done at a single spatial location, which inevitably results in spatially non-uniform ACh distribution. Such non-uniform application of ACh across the spatial extent of a cortical microcircuit could have important impacts on how the firing of groups of neurons is coordinated, but this remains largely unknown. Here we describe a method for applying ACh at different spatial locations within a single cortical circuit and measuring the resulting differences in population neural activity. We use two microdialysis probes implanted at opposite ends of a microelectrode array in barrel cortex of anesthetized rats. As a demonstration of the method, we applied ACh or neostigmine in different spatial locations via the microdialysis probes while we concomitantly recorded neural activity at 32 locations with the microelectrode array. First, we show that cholinergic changes in cortical state can vary dramatically depending on where the ACh was applied. Second, we show that cholinergic changes in cortical state can vary dramatically depending on where the state-change is measured. These results suggests that previous work with single-site recordings or single-site ACh application should be interpreted with some caution, since the results could change for different spatial locations.

Strong neuron-to-body coupling implies weak neuron-to-neuron coupling in motor cortex.

Cortical neurons can be strongly or weakly coupled to the network in which they are embedded, firing in sync with the majority or firing independently. Both these scenarios have potential computational advantages in motor cortex. Commands to the body might be more robustly conveyed by a strongly coupled population, whereas a motor code with greater information capacity could be implemented by neurons that fire more independently. Which of these scenarios prevails? Here we measure neuron-to-body coupling and neuron-to-population coupling for neurons in motor cortex of freely moving rats. We find that neurons with high and low population coupling coexist, and that population coupling was tunable by manipulating inhibitory signaling. Importantly, neurons with different population coupling tend to serve different functional roles. Those with strong population coupling are not involved with body movement. In contrast, neurons with high neuron-to-body coupling are weakly coupled to other neurons in the cortical population.

A theoretical framework for analyzing coupled neuronal networks: Application to the olfactory system.

Determining how synaptic coupling within and between regions is modulated during sensory processing is an important topic in neuroscience. Electrophysiological recordings provide detailed information about neural spiking but have traditionally been confined to a particular region or layer of cortex. Here we develop new theoretical methods to study interactions between and within two brain regions, based on experimental measurements of spiking activity simultaneously recorded from the two regions. By systematically comparing experimentally-obtained spiking statistics to (efficiently computed) model spike rate statistics, we identify regions in model parameter space that are consistent with the experimental data. We apply our new technique to dual micro-electrode array in vivo recordings from two distinct regions: olfactory bulb (OB) and anterior piriform cortex (PC). Our analysis predicts that: i) inhibition within the afferent region (OB) has to be weaker than the inhibition within PC, ii) excitation from PC to OB is generally stronger than excitation from OB to PC, iii) excitation from PC to OB and inhibition within PC have to both be relatively strong compared to presynaptic inputs from OB. These predictions are validated in a spiking neural network model of the OB-PC pathway that satisfies the many constraints from our experimental data. We find when the derived relationships are violated, the spiking statistics no longer satisfy the constraints from the data. In principle this modeling framework can be adapted to other systems and be used to investigate relationships between other neural attributes besides network connection strengths. Thus, this work can serve as a guide to further investigations into the relationships of various neural attributes within and across different regions during sensory processing.

State-dependent intrinsic predictability of cortical network dynamics.

The information encoded in cortical circuit dynamics is fleeting, changing from moment to moment as new input arrives and ongoing intracortical interactions progress. A combination of deterministic and stochastic biophysical mechanisms governs how cortical dynamics at one moment evolve from cortical dynamics in recently preceding moments. Such temporal continuity of cortical dynamics is fundamental to many aspects of cortex function but is not well understood. Here we study temporal continuity by attempting to predict cortical population dynamics (multisite local field potential) based on its own recent history in somatosensory cortex of anesthetized rats and in a computational network-level model. We found that the intrinsic predictability of cortical dynamics was dependent on multiple factors including cortical state, synaptic inhibition, and how far into the future the prediction extends. By pharmacologically tuning synaptic inhibition, we obtained a continuum of cortical states with asynchronous population activity at one extreme and stronger, spatially extended synchrony at the other extreme. Intermediate between these extremes we observed evidence for a special regime of population dynamics called criticality. Predictability of the near future (10-100 ms) increased as the cortical state was tuned from asynchronous to synchronous. Predictability of the more distant future (>1 s) was generally poor, but, surprisingly, was higher for asynchronous states compared to synchronous states. These experimental results were confirmed in a computational network model of spiking excitatory and inhibitory neurons. Our findings demonstrate that determinism and predictability of network dynamics depend on cortical state and the time-scale of the dynamics.

Maximizing Sensory Dynamic Range by Tuning the Cortical State to Criticality.

Modulation of interactions among neurons can manifest as dramatic changes in the state of population dynamics in cerebral cortex. How such transitions in cortical state impact the information processing performed by cortical circuits is not clear. Here we performed experiments and computational modeling to determine how somatosensory dynamic range depends on cortical state. We used microelectrode arrays to record ongoing and whisker stimulus-evoked population spiking activity in somatosensory cortex of urethane anesthetized rats. We observed a continuum of different cortical states; at one extreme population activity exhibited small scale variability and was weakly correlated, the other extreme had large scale fluctuations and strong correlations. In experiments, shifts along the continuum often occurred naturally, without direct manipulation. In addition, in both the experiment and the model we directly tuned the cortical state by manipulating inhibitory synaptic interactions. Our principal finding was that somatosensory dynamic range was maximized in a specific cortical state, called criticality, near the tipping point midway between the ends of the continuum. The optimal cortical state was uniquely characterized by scale-free ongoing population dynamics and moderate correlations, in line with theoretical predictions about criticality. However, to reproduce our experimental findings, we found that existing theory required modifications which account for activity-dependent depression. In conclusion, our experiments indicate that in vivo sensory dynamic range is maximized near criticality and our model revealed an unanticipated role for activity-dependent depression in this basic principle of cortical function.

Retronasal odor concentration coding in glomeruli of the rat olfactory bulb.

The mammalian olfactory system processes odorants presented orthonasally (inhalation through the nose) and also retronasally (exhalation), enabling identification of both external as well as internal objects during food consumption. There are distinct differences between ortho- and retronasal air flow patterns, psychophysics, multimodal integration, and glomerular responses. Recent work indicates that rats can also detect odors retronasally, that rats can associate retronasal odors with tastes, and that their olfactory bulbs (OBs) can respond to retronasal odorants but differently than to orthonasal odors. To further characterize retronasal OB input activity patterns, experiments here focus on determining the effects of odor concentration on glomerular activity by monitoring calcium activity in the dorsal OB of rats using a dextran-conjugated calcium-sensitive dye in vivo. Results showed reliable concentration-response curves that differed between odorants, and recruitment of additional glomeruli, as odor concentration increased. We found evidence of different concentration-response functions between glomeruli, that in turn depended on odor. Further, the relation between dynamics and concentration differed remarkably among retronasal odorants. These dynamics are suggested to reduce the odor map ambiguity based on response amplitude. Elucidating the coding of retronasal odor intensity is fundamental to the understanding of feeding behavior and the neural basis of flavor. These data further establish and refine the rodent model of flavor neuroscience.

Millisecond, micron precision multi-whisker detector.

The neural mechanisms of somatosensory information processing in the rodent vibrissae system are a topic of intense debate and research. Certain hypotheses emphasize the importance of stick-slip whisker motion, high-frequency resonant vibrations, and/or the ability to decode complex textures. Other hypotheses focus on the importance of integrating information from multiple whiskers. Tests of the former require measurements of whisker motion that achieve high spatiotemporal accuracy without altering the mechanical properties of whiskers. Tests of the latter require the ability to monitor the motion of multiple whiskers simultaneously. Here we present a device that achieves both these requirements for two-dimensional whisker motion in the plane perpendicular to the whiskers. Moreover, the system we present is significantly less expensive (<$2.5 k) and simpler to build than alternative devices which achieve similar detection capabilities. Our system is based on two laser diodes and two linear cameras. It attains millisecond temporal precision and micron spatial resolution. We developed automated algorithms for processing the data collected by our device and benchmarked their performance against manual detection by human visual inspection. By this measure, our detection was successful with less than 10 µm deviation between the automated and manual detection, on average. Here, we demonstrate its utility in anesthetized rats by measuring the motion of multiple whiskers in response to an air puff.

Direct behavioral evidence for retronasal olfaction in rats.

The neuroscience of flavor perception is becoming increasingly important to understand abnormal feeding behaviors and associated chronic diseases such as obesity. Yet, flavor research has mainly depended on human subjects due to the lack of an animal model. A crucial step towards establishing an animal model of flavor research is to determine whether the animal uses the retronasal mode of olfaction, an essential element of flavor perception. We designed a go- no go behavioral task to test the rat's ability to detect and discriminate retronasal odorants. In this paradigm, tasteless aqueous solutions of odorants were licked by water-restricted head-fixed rats from a lick spout. Orthonasal contamination was avoided by employing a combination of a vacuum around the lick-spout and blowing clean air toward the nose. Flow models support the effectiveness of both approaches. The licked odorants were successfully discriminated by rats. Moreover, the tasteless odorant amyl acetate was reliably discriminated against pure distilled water in a concentration-dependent manner. The results from this retronasal odor discrimination task suggest that rats are capable of smelling retronasally. This direct behavioral evidence establishes the rat as a useful animal model for flavor research.

Retronasal odor representations in the dorsal olfactory bulb of rats.

Animals perceive their olfactory environment not only from odors originating in the external world (orthonasal route) but also from odors released in the oral cavity while eating food (retronasal route). Retronasal olfaction is crucial for the perception of food flavor in humans. However, little is known about the retronasal stimulus coding in the brain. The most basic questions are if and how route affects the odor representations at the level of the olfactory bulb (OB), where odor quality codes originate. We used optical calcium imaging of presynaptic dorsal OB responses to odorants in anesthetized rats to ask whether the rat OB could be activated retronasally, and how these responses compare to orthonasal responses under similar conditions. We further investigated the effects of specific odorant properties on orthonasal versus retronasal response patterns. We found that at a physiologically relevant flow rate, retronasal odorants can effectively reach the olfactory receptor neurons, eliciting glomerular response patterns that grossly overlap with those of orthonasal responses, but differ from the orthonasal patterns in the response amplitude and temporal dynamics. Interestingly, such differences correlated well with specific odorant properties. Less volatile odorants yielded relatively smaller responses retronasally, but volatility did not affect relative temporal profiles. More polar odorants responded with relatively longer onset latency and time to peak retronasally, but polarity did not affect relative response magnitudes. These data provide insight into the early stages of retronasal stimulus coding and establish relationships between orthonasal and retronasal odor representations in the rat OB.

Taste quality and intensity of 100 stimuli as reported by rats: the taste-location association task.

The interpretation of neural activity related to sensory stimulation requires an understanding of the subject's perception of the stimulation. Previous methods used to evaluate the perception of chemosensory stimuli by rodents have distinct limitations. We developed a novel behavioral paradigm, the taste-location association task, to complement these methods. First we tested if rats are able to learn associations between five basic taste stimuli and their spatial locations. This spatial task was based on four prototypical tastants and water. All four rats trained to perform the task reached levels of performance well above chance. Control trials demonstrated that the rats used only taste cues. Further, the learned stimulus set was resistant to interference, allowing for generalization experiments performed subsequently. We tested the rats' gustatory generalizations of 100 tastants to the five trained stimuli, both regarding their taste qualities as well as intensity ratings. The taste profiles generated by these experiments contribute to the understanding of how perception of the specific taste stimuli relate to the perception of the five basic taste qualities in intact behaving rats. In this large taste space we found that intensity plays a major role. Furthermore, umami stimuli were not reported as being similar to other basic tastants. Our new paradigm enables neurophysiological studies of taste-based learning and memory in awake, freely moving animals.

Evidence that the sweetness of odors depends on experience in rats.

Humans describe their perception of certain odorants in terms of taste qualities (e.g., sweet). It has also been found that in humans, novel odorants can rapidly and irreversibly acquire a taste, even after just a single pairing with a taste. It remains unclear whether flavor objects in general, and odor-taste generalizations in particular, are experience-dependent. Interactions might result from a failure by humans to sufficiently analyze the olfactory and gustatory components of compound flavorants. Here, we tested odor-taste generalizations in rats with or without paired exposure to an odorant and a tastant. We evaluated the generalization of conditioned odor aversion to tastants by rats. Our findings suggest that rats behave toward putatively tasteless retronasal odorants as if they were sweet only after prior paired experience of the odorant with a sweet tastant. These data support the hypothesis that taste-like qualities of odors are learned and are not innate. Furthermore, the present results suggest that acquisition of a taste quality by an odor need not depend on higher cognitive abilities. This study helps to establish the rat as a model for the study of behavioral neuroscience of flavor.

Voltage-gated channel properties of epithelial cells in porcine vomeronasal organ.

Bipolar vomeronasal sensory neurons (VSNs) in the vomeronasal organ (VNO) are believed to detect pheromones in most mammals. The vomeronasal sensory epithelium (VSE) is composed of VSNs and supporting cells. There are morphological differences in VNOs between species. Many electrophysiological experiments have been performed on rodent VSEs but few on other mammals. We therefore investigated voltage-gated channel properties of cells in the porcine VSE using slice whole-cell voltage-clamp techniques. In immunohistochemical study of the porcine VSE, most PGP9.5-immunoreactive cells were found between the middle and basal region, and negative cells were distributed in the apical to middle region. Depolarizing pulses to epithelial cells from -90mV produced transient inward Na+ channel currents and sustained outward K+ channel currents with various amplitudes. The distribution of cells having high and low Na+ current densities was mostly consistent with the histological distribution of VSNs and supporting cells, respectively. The half-inactivation voltage of voltage-gated Na+ channels in supporting cells was 26mV more negative than that in VSNs. Voltage-gated K+ channel currents in both cell types were suppressed by tetraethylammonium to the same extent. VSNs possessed TTX-sensitive voltage-gated Na+ channels and Ni2+ -sensitive T-type Ca2+ channels. These results suggest that the histological distribution of porcine vomeronasal epithelial cells is more similar to the dog and goat than to rodents, and that the electrophysiological characteristics of porcine vomeronasal epithelial cells are similar to those of rodents. It is also suggested that porcine VSNs detecting pheromones generate action potentials through these channels.

T-type Ca2+ channels contribute to IBMX/forskolin- and K(+)-induced Ca(2+) transients in porcine olfactory receptor neurons.

T-type Ca(2+) channels are low-voltage-activated Ca(2+) channels that control Ca(2+) entry in excitable cells during small depolarization above resting potentials. Using Ca(2+) imaging with a laser scanning confocal microscope we investigated the involvement of T-type Ca(2+) channels in IBMX/forskolin- and sparingly elevated extracellular K(+)-induced Ca(2+) transients in freshly isolated porcine olfactory receptor neurons (ORNs). In the presence of mibefradil (10microM) or Ni(2+) (100microM), the selective T-type Ca(2+) channel inhibitors, IBMX/forskolin-induced Ca(2+) transients in the soma were either strongly (>60%) inhibited or abolished completely. However, the Ca(2+) transients in the knob were only partially (<60%) inhibited. Ca(2+) transients induced by 30mM K(+) were also partially ( approximately 60%) inhibited at both the knob and soma. Furthermore, ORNs responded to as little as a 2.5mM increase in the extracellular K(+) concentration (7.5mM K(+)), and such responses were completely inhibited by mibefradil or Ni(2+). These results reveal functional expression of T-type Ca(2+) channels in porcine ORNs, and suggest a role for these channels in the spread Ca(2+) transients from the knob to the soma during activation of the cAMP cascade following odorant binding to G-protein-coupled receptors on the cilia/knob of ORNs.