Publications

2026

Statistics of natural scenes shape contextual modulation in the visual cortex

Journal

Fu, Jiakun; Shrinivasan, Suhas; Baroni, Luca; Ding, Zhuokun; Fahey, Paul G.; Pierzchlewicz, Paweł A.; Karantzas, Nikos; Ponder, Kayla; Froebe, Rachel; Ntanavara, Lydia; Muhammad, Taliah; Willeke, Konstantin F.; Wang, Eric; Ding, Zhiwei; Tran, Dat T.; Papadopoulos, Stelios; Patel, Saumil; Reimer, Jacob; Ecker, Alexander S.; Pitkow, Xaq; Antolík, Ján; Sinz, Fabian H.; Haefner, Ralf M.; Tolias, Andreas S.; Franke, Katrin

Neuron

Abstract

Vision is context dependent, with neuronal responses shaped not only by local features but also by surrounding visual input. While classical studies, using grating stimuli, show that iso-oriented surrounds suppress responses more than orthogonal surrounds, the role of contextual modulation under natural stimulus conditions remains less clear. Using recordings from mouse primary visual cortex (V1), we trained convolutional neural network models to predict neuronal responses to natural images and synthesized surround stimuli that selectively suppressed or facilitated responses to optimal center inputs. In vivo experiments confirmed these predictions. Facilitatory surrounds resembled naturalistic continuations of the optimal center stimulus, consistent with natural image statistics, whereas suppressive surrounds deviated from these predictions. Applying the same approach to macaque V1 revealed similar principles across species. Both models accurately predicted responses to classical grating stimuli. We formalize these results in a normative Bayesian model, showing that neuronal activity for preferred center features reflects posterior beliefs about likely center-surround configurations.

Functional bipartite invariance in mouse primary visual cortex receptive fields

Journal

Ding, Zhiwei; Tran, Dat; Ponder, Kayla; Ding, Zhuokun; Froebe, Rachel; Ntanavara, Lydia; Fahey, Paul G.; Cobos, Erick; Baroni, Luca; Diamantaki, Maria; Wang, Eric Y.; Chang, Andersen; Papadopoulos, Stelios; Fu, Jiakun; Muhammad, Taliah; Papadopoulos, Christos; Cadena, Santiago A.; Evangelou, Alexandros; Willeke, Konstantin; Anselmi, Fabio; Sanborn, Sophia; Antolík, Ján; Froudarakis, Emmanouil; Patel, Saumil; Walker, Edgar Y.; Reimer, Jacob; Sinz, Fabian H.; Ecker, Alexander S.; Franke, Katrin; Pitkow, Xaq; Tolias, Andreas S.

Nature Neuroscience

Abstract

Sensory systems support generalization by representing features that persist under input variation; however, identifying the neuronal basis of these invariances remains difficult due to high-dimensional and nonlinear neural computations. Here we leverage the inception loop paradigm, iterating between large-scale recordings, predictive models and in silico experiments with in vivo verification, to characterize neuronal invariances in mouse primary visual cortex (V1). We synthesize varied exciting inputs (VEIs), dissimilar images that drive target neurons. These VEIs revealed a new bipartite invariance: one subfield encodes a shift-tolerant high-frequency texture and the other encodes a fixed low-frequency pattern. This division aligns with object boundaries defined by spatial frequency differences in highly activating images, suggesting a contribution to segmentation. Analysis of the MICrONS dataset revealed a hierarchy of excitatory neurons in mouse V1 layers 2/3: postsynaptic neurons exhibited greater invariance than their presynaptic inputs, while neurons with lower invariance formed more connections. Together, these results provide insights and scalable methodology for mapping neuronal invariances.

Iso-orientation bias of layer 2/3 connections unifies spontaneous, visually and optogenetically driven V1 dynamics

Journal

Rózsa, Tibor; Cagnol, Rémy; Antolík, Ján

Nature Communications

Abstract

Functionally specific long-range lateral connectivity in layer 2/3 of the primary visual cortex (V1) supports the integration of visual information across visual space and shapes spontaneous, visual and optogenetically driven V1 activity. However, a comprehensive understanding of how these diverse cortical regimes emerge from this underlying cortical circuitry remains elusive. Here we address this gap by showing how the same model assuming moderately iso-orientation biased long-range cortical connectivity architecture explains diverse phenomena, including (i) a range of visually driven phenomena, (ii) modular spontaneous activity, (iii) the propagation of spontaneous cortical waves, and (iv) neural responses to patterned optogenetic stimulation. The model offers testable predictions, including presence of slower and iso-tropic spontaneous wave propagation in layer 4 and non-monotonicity of optogenetically driven cortical response to increasingly larger disk of illumination. We thus offer a holistic framework for studying how cortical circuitry governs information integration across multiple operating regimes.

MEIcoder: Decoding Visual Stimuli from Neural Activity by Leveraging Most Exciting Inputs

Conference

Sobotka, Jan; Baroni, Luca; Antolík, Ján

Advances in Neural Information Processing Systems

Abstract

Decoding visual stimuli from neural population activity is crucial for understanding the brain and for applications in brain-machine interfaces. However, such biological data is often scarce, particularly in primates or humans, where high-throughput recording techniques, such as two-photon imaging, remain challenging or impossible to apply. This, in turn, poses a challenge for deep learning decoding techniques. To overcome this, we introduce MEIcoder, a biologically informed decoding method that leverages neuron-specific most exciting inputs (MEIs), a structural similarity index measure loss, and adversarial training. MEIcoder achieves state-of-the-art performance in reconstructing visual stimuli from single-cell activity in primary visual cortex (V1), especially excelling on small datasets with fewer recorded neurons. Using ablation studies, we demonstrate that MEIs are the main drivers of the performance, and in scaling experiments, we show that MEIcoder can reconstruct high-fidelity natural-looking images from as few as 1,000-2,500 neurons and less than 1,000 training data points. We also propose a unified benchmark with over 160,000 samples to foster future research. Our results demonstrate the feasibility of reliable decoding in early visual system and provide practical insights for neuroscience and neuroengineering applications.

2025

An open-source replication for fast and accessible light propagation modeling in brain tissue

Journal

Berling, David; Střeleček, Jan; Iser, Tomáš; Antolík, Ján

PLOS One

Abstract

Light simulations hold great potential for advancing optical techniques in neuroscience. They facilitate the in-silico refinement of optical stimulator designs and enable simulations of optical recordings from computational brain models, aiding neuroscience in forming a mechanistic understanding of brain circuitry. However, many published light models are inaccessible due to unavailable source code and documentation or are impractical due to excessive computational demands. To address these challenges, we replicate and enhance the efficient and accurate light simulation model by Yona et al. [1], which was previously available only in compiled form accompanied by sparse documentation. In this work, we resolve ambiguities in the original model, correct errors that caused discrepancies between simulations and published results, improve computational efficiency by an order of magnitude, and open-source all the resulting code and detailed documentation. These enhancements enable simulations of cortical volumes exceeding 1 mm3 to run in seconds on standard laptop hardware. Our model software provides an accessible, adaptable, and rapid light simulation tool, which adheres to FAIR principles to ensure future-proof and broad utility for the neuroscience community.

Learning and Aligning Single-Neuron Invariance Manifolds in Visual Cortex

Conference

Bashiri, Mohammad; Baroni, Luca; Antolík, Ján; Sinz, Fabian H.

International Conference on Learning Representations (ICLR) — Oral presentation; equal contribution: Bashiri and Baroni

Abstract

Understanding how sensory neurons exhibit selectivity to certain features and invariance to others is central to uncovering the computational principles underlying robustness and generalization in visual perception. Most existing methods for characterizing selectivity and invariance identify single or finite discrete sets of stimuli. Since these are only isolated measurements from an underlying continuous manifold, characterizing invariance properties accurately and comparing them across neurons with varying receptive field size, position, and orientation, becomes challenging. Consequently, a systematic analysis of invariance types at the population level remains under-explored. Building on recent advances in learning continuous invariance manifolds, we introduce a novel method to accurately identify and align invariance manifolds of visual sensory neurons, overcoming these challenges. Our approach first learns the continuous invariance manifold of stimuli that maximally excite a neuron modeled by a response-predicting deep neural network. It then learns an affine transformation on the pixel coordinates such that the same manifold activates another neuron as strongly as possible, effectively aligning their invariance manifolds spatially. This alignment provides a principled way to quantify and compare neuronal invariances irrespective of receptive field differences. Using simulated neurons, we demonstrate that our method accurately learns and aligns known invariance manifolds, robustly identifying functional clusters. When applied to macaque V1 neurons, it reveals functional clusters of neurons, including simple and complex cells. Overall, our method enables systematic, quantitative exploration of the neural invariance landscape, to gain new insights into the functional properties of visual sensory neurons.

Non-linearity of spatial integration varies across layers of primary visual cortex

Preprint

Cagnol, Rémy; Antolík, Ján; Palmer, Larry A.; Contreras, Diego

bioRxiv

Abstract

The receptive field (RF) of visual cortical neurons is highly dynamic and context-dependent, shaped by both the spatial and temporal properties of stimuli and the complex architecture of cortical circuits. While classical RF mapping through extracellular recordings reveals only the area triggering spiking responses, intracellular recordings reveal a much broader region of subthreshold synaptic input. We investigated how neurons in different cortical layers integrate visual input across space, with a focus on the linearity of spatial summation. Using intracellular recordings, we found that supragranular complex cells integrate input in a highly sublinear manner, in contrast to infragranular complex cells and simple cells, which exhibited near-linear summation. To understand the underlying mechanisms, we employed a large-scale recurrent spiking model of cat primary visual cortex (V1). Modeling results point to the differential patterning of long-range horizontal connections-particularly their targeting of excitatory versus inhibitory neurons-as a potential source of the observed layer-specific integration properties. These findings suggest that RFs emerge from interaction of feedforward, horizontal, and possibly feedback inputs, that are continuous in space, challenging the conventional notions of fixed spatial RF boundaries in early visual processing.

Efficient incorporation of dendrites into a large-scale cortical model reveals their surprising role in sharpening optogenetic responses

Preprint

Berling, David; Rózsa, Tibor; Wybo, Willem; Antolík, Ján

bioRxiv

Abstract

Single-photon optogenetics enables chronic wide-field stimulation of cortex, facilitating large-scale manipulation of neural code to study cortical processing and advance neuroprosthetics. However, access to neural codes organized at fine spatial scales is compromised by the horizontal spread of stimulation-evoked cortical activity. Overcoming this limitation requires a quantitative understanding of the mechanisms contributing to spread, which include light scattering, dendritic activation, and synaptic transmission. We addressed this with morphology-aware simulations of optogenetic stimulation in a functionally-detailed network model of primary visual cortex. We find that synaptic transmission extends activation by 37–50% beyond the illuminated area, while, paradoxically, neuronal morphology sharpens activation, as apical dendrites sample from the superficial cortex, which is less affected by light dispersion. This unexpected sharpening enhances the fidelity of stimulation with spatially distributed patterns. Our study offers guidance for optogenetic interventions targeting topographically organized neural codes and provides a computational testbed to interpret such experiments.

Beyond retinotopy: exploiting native visual representations in cortical neuroprostheses for vision loss remediation

Preprint

Baroni, Luca; Picek, Martin; Patel, Saumil; Tolias, Andreas S.; Antolík, Ján

bioRxiv

Abstract

Cortical prosthetic systems offer a promising path to restoring vision in blindness by stimulating neurons in the visual cortex to evoke visual percepts. To effectively encode visual information, however, it is essential to target neurons based on their functional encoding properties. Current stimulation protocols focus on retinotopic information, ignoring other key encoding properties, and therefore fail to reproduce complex visual percepts. We demonstrate that incorporating orientation selectivity alongside retinotopy to guide stimulation dramatically improves fidelity of the evoked activity to the underlying neural code. We propose a Bottlenecked Rotation-Equivariant CNN (BRCNN) and demonstrate that neural responses can be predicted to a large degree using only retinotopy and orientation preference. Using this model, we design a retinotopy-and-orientation aware stimulation protocol and validate it in a state-of-the-art large-scale optogenetic stimulation simulation framework of primary visual cortex. Our protocol elicits neural activity patterns exhibiting substantially higher correlation with natural vision responses compared to retinotopic-only approaches.

2024

A comprehensive data-driven model of cat primary visual cortex

Journal

Antolík, Ján; Cagnol, Rémy; Rózsa, Tibor; Monier, Cyril; Frégnac, Yves; Davison, Andrew P.

PLOS Computational Biology

Abstract

Knowledge integration based on the relationship between structure and function of the neural substrate is one of the main targets of neuroinformatics and data-driven computational modeling. However, the multiplicity of data sources, the diversity of benchmarks, the mixing of observables of different natures, and the necessity of a long-term, systematic approach make such a task challenging. Here we present a first snapshot of a long-term integrative modeling program designed to address this issue in the domain of the visual system: a comprehensive spiking model of cat primary visual cortex. The presented model satisfies an extensive range of anatomical, statistical and functional constraints under a wide range of visual input statistics. In the presence of physiological levels of tonic stochastic bombardment by spontaneous thalamic activity, the modeled cortical reverberations self-generate a sparse asynchronous ongoing activity that quantitatively matches a range of experimentally measured statistics. When integrating feed-forward drive elicited by a high diversity of visual contexts, the simulated network produces a realistic, quantitatively accurate interplay between visually evoked excitatory and inhibitory conductances; contrast-invariant orientation-tuning width; center surround interactions; and stimulus-dependent changes in the precision of the neural code. This integrative model offers insights into how the studied properties interact, contributing to a better understanding of visual cortical dynamics. It provides a basis for future development towards a comprehensive model of low-level perception.

Optogenetic Stimulation Recruits Cortical Neurons in a Morphology-Dependent Manner

Journal

Berling, David; Baroni, Luca; Chaffiol, Antoine; Gauvain, Gregory; Picaud, Serge; Antolík, Ján

Journal of Neuroscience

Abstract

Single-photon optogenetics enables precise, cell-type–specific modulation of neuronal circuits, making it a crucial tool in neuroscience. Its miniaturization in the form of fully implantable wide-field stimulator arrays enables long-term interrogation of cortical circuits and bears promise for brain–machine interfaces for sensory and motor function restoration. However, achieving selective activation of functional cortical representations poses a challenge, as studies show that targeted optogenetic stimulation results in activity spread beyond one functional domain. While recurrent network mechanisms contribute to activity spread, here we demonstrate with detailed simulations of isolated pyramidal neurons from cats of unknown sex that already neuron morphology causes a complex spread of optogenetic activity at the scale of one cortical column. Since the shape of a neuron impacts its optogenetic response, we find that a single stimulator at the cortical surface recruits a complex spatial distribution of neurons that can be inhomogeneous and vary with stimulation intensity and neuronal morphology across layers. We explore strategies to enhance stimulation precision, finding that optimizing stimulator optics may offer more significant improvements than the preferentially somatic expression of the opsin through genetic targeting. Our results indicate that, with the right optical setup, single-photon optogenetics can precisely activate isolated neurons at the scale of functional cortical domains spanning several hundred micrometers.

2023

Learning invariance manifolds of visual sensory neurons

Conference

Baroni, Luca; Bashiri, Mohammad; Willeke, Konstantin F.; Antolík, Ján; Sinz, Fabian H.

Proceedings of the 1st NeurIPS Workshop on Symmetry and Geometry in Neural Representations

Abstract

Robust object recognition is thought to rely on neural mechanisms that are selective to complex stimulus features while being invariant to others (e.g., spatial location or orientation). To better understand biological vision, it is thus crucial to characterize which features neurons in different visual areas are selective or invariant to. In the past, invariances have commonly been identified by presenting carefully selected hypothesis-driven stimuli which rely on the intuition of the researcher. One example is the discovery of phase invariance in V1 complex cells. However, to identify novel invariances, a data-driven approach is more desirable. Here, we present a method that, combined with a predictive model of neural responses, learns a manifold in the stimulus space along which a target neuron’s response is invariant. Our approach is fully data-driven, allowing the discovery of novel neural invariances, and enables scientists to generate and experiment with novel stimuli along the invariance manifold. We test our method on Gabor-based neuron models as well as on a neural network fitted on macaque V1 responses and show that 1) it successfully identifies neural invariances, and 2) disentangles invariant directions in the stimulus space.

2022

The cortico-thalamic loop attunes competitive lateral interactions across retinotopic and orientation preference maps

Preprint

Guarino, Domenico G.; Davison, Andrew P.; Frégnac, Yves; Antolík, Ján

bioRxiv — Scholar uses year 2022; bioRxiv v2 page is dated 2023-01-24.

Abstract

In the early visual system, corticothalamic feedback projections greatly outnumber thalamocortical feedforward projections. Extensive experimental and modeling work has been devoted to the functional impact of the feedforward pathway, but the role of its denser feedback counterpart remains elusive. Here, we propose a novel unifying framework where thalamic recurrent interactions and corticothalamic feedback act in a closed-loop fashion to attune multiple stimulus representations. At each position of the visual field, the loop puts into competition local representations of the stimulus in thalamus and cortex through direct excitation of narrow topologically-aligned portions of the thalamus, accompanied with peri-geniculate nucleus mediated broad inhibition suppressing the topological surround. We built a detailed conductance-based spiking model incorporating retinal input, lateral geniculate nucleus, peri-geniculate nucleus, primary visual cortex, and all the relevant intra-areal and feedback pathways. For the first time we perform comparative analyses between model configurations with completely or locally inactivated cortico-thalamic feedback, as in the experimental preparations. The model mechanistically explains (i) the existence of intra-thalamic surround suppression, (ii) the sensitivity of thalamic neurons to orientation tuning, (iii) the cortex-dependent center-surround opponency in thalamic cells, (iv) the cortical increase of size and orientation selectivity, (v) the cortically enhanced competition between cross-oriented domains within the hypercolumn, and (vi) the selective suppression of cortical functional connectivity. Our results integrate decades of experimental and theoretical research, supporting the hypothesis that cortico-thalamic loop exerts competitive influence between neighboring regions in the thalamus and cortex, complementing the lateral intra-V1 interactions in center-surround contextual modulation.

2021

Assessment of optogenetically-driven strategies for prosthetic restoration of cortical vision in large-scale neural simulation of V1

Journal

Antolík, Ján; Sabatier, Quentin; Galle, Charlie; Frégnac, Yves; Benosman, Ryad

Scientific Reports

Abstract

The neural encoding of visual features in primary visual cortex (V1) is well understood, with strong correlates to low-level perception, making V1 a strong candidate for vision restoration through neuroprosthetics. However, the functional relevance of neural dynamics evoked through external stimulation directly imposed at the cortical level is poorly understood. Furthermore, protocols for designing cortical stimulation patterns that would induce a naturalistic perception of the encoded stimuli have not yet been established. Here, we demonstrate a proof of concept by solving these issues through a computational model, combining (1) a large-scale spiking neural network model of cat V1 and (2) a virtual prosthetic system transcoding the visual input into tailored light-stimulation patterns which drive in situ the optogenetically modified cortical tissue. Using such virtual experiments, we design a protocol for translating simple Fourier contrasted stimuli (gratings) into activation patterns of the optogenetic matrix stimulator. We then quantify the relationship between spatial configuration of the imposed light pattern and the induced cortical activity. Our simulations in the absence of visual drive (simulated blindness) show that optogenetic stimulation with a spatial resolution as low as 100 μm, and light intensity as weak as 10 ^ 16 photons/s/cm is sufficient to evoke activity patterns in V1 close to those evoked by normal vision.

An anatomically constrained model of V1 simple cells predicts the coexistence of push-pull and broad inhibition

Journal

Taylor, M. Morgan; Contreras, Diego; Destexhe, Alain; Frégnac, Yves; Antolík, Ján

Journal of Neuroscience

Abstract

The spatial organization and dynamic interactions between excitatory and inhibitory synaptic inputs that define the receptive field (RF) of simple cells in the cat primary visual cortex (V1) still raise the following paradoxical issues: (1) stimulation of simple cells in V1 with drifting gratings supports a wiring schema of spatially segregated sets of excitatory and inhibitory inputs activated in an opponent way by stimulus contrast polarity and (2) in contrast, intracellular studies using flashed bars suggest that although ON and OFF excitatory inputs are indeed segregated, inhibitory inputs span the entire RF regardless of input contrast polarity. Here, we propose a biologically detailed computational model of simple cells embedded in a V1-like network that resolves this seeming contradiction. We varied parametrically the RF-correlation-based bias for excitatory and inhibitory synapses and found that a moderate bias of excitatory neurons to synapse onto other neurons with correlated receptive fields and a weaker bias of inhibitory neurons to synapse onto other neurons with anticorrelated receptive fields can explain the conductance input, the postsynaptic membrane potential, and the spike train dynamics under both stimulation paradigms. This computational study shows that the same structural model can reproduce the functional diversity of visual processing observed during different visual contexts.

2018

Arkheia: Data Management and Communication for Open Computational Neuroscience

Journal

Antolík, Ján; Davison, Andrew P.

Frontiers in Neuroinformatics

Abstract

Two trends have been unfolding in computational neuroscience during the last decade. First, a shift of focus to increasingly complex and heterogeneous neural network models, with a concomitant increase in the level of collaboration within the field (whether direct or in the form of building on top of existing tools and results). Second, a general trend in science toward more open communication, both internally, with other potential scientific collaborators, and externally, with the wider public. This multi-faceted development toward more integrative approaches and more intense communication within and outside of the field poses major new challenges for modelers, as currently there is a severe lack of tools to help with automatic communication and sharing of all aspects of a simulation workflow to the rest of the community. To address this important gap in the current computational modeling software infrastructure, here we introduce Arkheia. Arkheia is a web-based open science platform for computational models in systems neuroscience. It provides an automatic, interactive, graphical presentation of simulation results, experimental protocols, and interactive exploration of parameter searches, in a web browser-based application. Arkheia is focused on automatic presentation of these resources with minimal manual input from users. Arkheia is written in a modular fashion with a focus on future development of the platform. The platform is designed in an open manner, with a clearly defined and separated API for database access, so that any project can write its own backend translating its data into the Arkheia database format. Arkheia is not a centralized platform, but allows any user (or group of users) to set up their own repository, either for public access by the general population, or locally for internal use. Overall, Arkheia provides users with an automatic means to communicate information about not only their models but also individual simulation results and the entire experimental context in an approachable graphical manner, thus facilitating the user's ability to collaborate in the field and outreach to a wider audience.

2017

Rapid long-range disynaptic inhibition explains the formation of cortical orientation maps

Journal

Antolík, Ján

Frontiers in Neural Circuits

Abstract

Competitive interactions are believed to underlie many types of cortical processing, ranging from memory formation, attention and development of cortical functional organization (e.g., development of orientation maps in primary visual cortex). In the latter case, the competitive interactions happen along the cortical surface, with local populations of neurons reinforcing each other, while competing with those displaced more distally. This specific configuration of lateral interactions is however in stark contrast with the known properties of the anatomical substrate, i.e., excitatory connections (mediating reinforcement) having longer reach than inhibitory ones (mediating competition). No satisfactory biologically plausible resolution of this conflict between anatomical measures, and assumed cortical function has been proposed. Recently a specific pattern of delays between different types of neurons in cat cortex has been discovered, where direct mono-synaptic excitation has approximately the same delay, as the combined delays of the disynaptic inhibitory interactions between excitatory neurons (i.e., the sum of delays from excitatory to inhibitory and from inhibitory to excitatory neurons). Here we show that this specific pattern of delays represents a biologically plausible explanation for how short-range inhibition can support competitive interactions that underlie the development of orientation maps in primary visual cortex. We demonstrate this statement analytically under simplifying conditions, and subsequently show using network simulations that development of orientation maps is preserved when long-range excitation, direct inhibitory to inhibitory interactions, and moderate inequality in the delays between excitatory and inhibitory pathways is added.

2016

Model constrained by visual hierarchy improves prediction of neural responses to natural scenes

Journal

Antolík, Ján; Hofer, Sonja B.; Bednar, James A.; Mrsic-Flogel, Thomas D.

PLOS Computational Biology

Abstract

Accurate estimation of neuronal receptive fields is essential for understanding sensory processing in the early visual system. Yet a full characterization of receptive fields is still incomplete, especially with regard to natural visual stimuli and in complete populations of cortical neurons. While previous work has incorporated known structural properties of the early visual system, such as lateral connectivity, or imposing simple-cell-like receptive field structure, no study has exploited the fact that nearby V1 neurons share common feed-forward input from thalamus and other upstream cortical neurons. We introduce a new method for estimating receptive fields simultaneously for a population of V1 neurons, using a model-based analysis incorporating knowledge of the feed-forward visual hierarchy. We assume that a population of V1 neurons shares a common pool of thalamic inputs, and consists of two layers of simple and complex-like V1 neurons. When fit to recordings of a local population of mouse layer 2/3 V1 neurons, our model offers an accurate description of their response to natural images and significant improvement of prediction power over the current state-of-the-art methods. We show that the responses of a large local population of V1 neurons with locally diverse receptive fields can be described with surprisingly limited number of thalamic inputs, consistent with recent experimental findings. Our structural model not only offers an improved functional characterization of V1 neurons, but also provides a framework for studying the relationship between connectivity and function in visual cortical areas.

2013

Mechanisms for stable, robust, and adaptive development of orientation maps in the primary visual cortex

Journal

Stevens, Jean-Luc R.; Law, Judith S.; Antolík, Ján; Bednar, James A.

Journal of Neuroscience

Abstract

Development of orientation maps in ferret and cat primary visual cortex (V1) has been shown to be stable, in that the earliest measurable maps are similar in form to the eventual adult map, robust, in that similar maps develop in both dark rearing and in a variety of normal visual environments, and yet adaptive, in that the final map pattern reflects the statistics of the specific visual environment. How can these three properties be reconciled? Using mechanistic models of the development of neural connectivity in V1, we show for the first time that realistic stable, robust, and adaptive map development can be achieved by including two low-level mechanisms originally motivated from single-neuron results. Specifically, contrast-gain control in the retinal ganglion cells and the lateral geniculate nucleus reduces variation in the presynaptic drive due to differences in input patterns, while homeostatic plasticity of V1 neuron excitability reduces the postsynaptic variability in firing rates. Together these two mechanisms, thought to be applicable across sensory systems in general, lead to biological maps that develop stably and robustly, yet adapt to the visual environment. The modeling results suggest that topographic map stability is a natural outcome of low-level processes of adaptation and normalization. The resulting model is more realistic, simpler, and far more robust, and is thus a good starting point for future studies of cortical map development.

The emergence of functional microcircuits in visual cortex

Journal

Ko, Ho; Cossell, Lee; Baragli, Chiara; Antolík, Ján; Clopath, Claudia; Hofer, Sonja B.; Mrsic-Flogel, Thomas D.

Nature

Abstract

Sensory processing occurs in neocortical microcircuits in which synaptic connectivity is highly structured and excitatory neurons form subnetworks that process related sensory information. However, the developmental mechanisms underlying the formation of functionally organized connectivity in cortical microcircuits remain unknown. Here we directly relate patterns of excitatory synaptic connectivity to visual response properties of neighbouring layer 2/3 pyramidal neurons in mouse visual cortex at different postnatal ages, using two-photon calcium imaging in vivo and multiple whole-cell recordings in vitro. Although neural responses were already highly selective for visual stimuli at eye opening, neurons responding to similar visual features were not yet preferentially connected, indicating that the emergence of feature selectivity does not depend on the precise arrangement of local synaptic connections. After eye opening, local connectivity reorganized extensively: more connections formed selectively between neurons with similar visual responses and connections were eliminated between visually unresponsive neurons, but the overall connectivity rate did not change. We propose a sequential model of cortical microcircuit development based on activity-dependent mechanisms of plasticity whereby neurons first acquire feature preference by selecting feedforward inputs before the onset of sensory experience--a process that may be facilitated by early electrical coupling between neuronal subsets--and then patterned input drives the formation of functional subnetworks through a redistribution of recurrent synaptic connections.

Integrated workflows for spiking neuronal network simulations

Journal

Antolík, Ján; Davison, Andrew P.

Frontiers in Neuroinformatics

Abstract

The increasing availability of computational resources is enabling more detailed, realistic modeling in computational neuroscience, resulting in a shift toward more heterogeneous models of neuronal circuits, and employment of complex experimental protocols. This poses a challenge for existing tool chains, as the set of tools involved in a typical modeler's workflow is expanding concomitantly, with growing complexity in the metadata flowing between them. For many parts of the workflow, a range of tools is available; however, numerous areas lack dedicated tools, while integration of existing tools is limited. This forces modelers to either handle the workflow manually, leading to errors, or to write substantial amounts of code to automate parts of the workflow, in both cases reducing their productivity. To address these issues, we have developed Mozaik: a workflow system for spiking neuronal network simulations written in Python. Mozaik integrates model, experiment and stimulation specification, simulation execution, data storage, data analysis and visualization into a single automated workflow, ensuring that all relevant metadata are available to all workflow components. It is based on several existing tools, including PyNN, Neo, and Matplotlib. It offers a declarative way to specify models and recording configurations using hierarchically organized configuration files. Mozaik automatically records all data together with all relevant metadata about the experimental context, allowing automation of the analysis and visualization stages. Mozaik has a modular architecture, and the existing modules are designed to be extensible with minimal programming effort. Mozaik increases the productivity of running virtual experiments on highly structured neuronal networks by automating the entire experimental cycle, while increasing the reliability of modeling studies by relieving the user from manual handling of the flow of metadata between the individual workflow stages.

2011

Development of maps of simple and complex cells in the primary visual cortex

Journal

Antolík, Ján; Bednar, James A.

Frontiers in Computational Neuroscience

Abstract

Hubel and Wiesel (1962) classified primary visual cortex (V1) neurons as either simple, with responses modulated by the spatial phase of a sine grating, or complex, i.e., largely phase invariant. Much progress has been made in understanding how simple-cells develop, and there are now detailed computational models establishing how they can form topographic maps ordered by orientation preference. There are also models of how complex cells can develop using outputs from simple cells with different phase preferences, but no model of how a topographic orientation map of complex cells could be formed based on the actual connectivity patterns found in V1. Addressing this question is important, because the majority of existing developmental models of simple-cell maps group neurons selective to similar spatial phases together, which is contrary to experimental evidence, and makes it difficult to construct complex cells. Overcoming this limitation is not trivial, because mechanisms responsible for map development drive receptive fields (RF) of nearby neurons to be highly correlated, while co-oriented RFs of opposite phases are anti-correlated. In this work, we model V1 as two topographically organized sheets representing cortical layer 4 and 2/3. Only layer 4 receives direct thalamic input. Both sheets are connected with narrow feed-forward and feedback connectivity. Only layer 2/3 contains strong long-range lateral connectivity, in line with current anatomical findings. Initially all weights in the model are random, and each is modified via a Hebbian learning rule. The model develops smooth, matching, orientation preference maps in both sheets. Layer 4 units become simple cells, with phase preference arranged randomly, while those in layer 2/3 are primarily complex cells. To our knowledge this model is the first explaining how simple cells can develop with random phase preference, and how maps of complex cells can develop, using only realistic patterns of connectivity.

2005

Evolutionary tree genetic programming

Conference

Antolík, Ján; Hsu, William H.

Proceedings of the 2005 Genetic and Evolutionary Computation Conference (GECCO '05)

Abstract

We introduce a clustering-based method of subpopulation management in genetic programming (GP) called Evolutionary Tree Genetic Programming (ETGP). The biological motivation behind this work is the observation that the natural evolution follows a tree-like phylogenetic pattern. Our goal is to simulate similar behavior in artificial evolutionary systems such as GP. To test our model we use three common GP benchmarks: the Ant Algorithm, 11-Multiplexer, and Parity problems.The performance of the ETGP system is empirically compared to those of the GP system. Code size and variance are consistently reduced by a small but statistically significant percentage, resulting in a slight speedup in the Ant and 11-Multiplexer problems, while the same comparisons on the Parity problem are inconclusive.

Automatic annotation of medical records

Journal

Antolík, Ján

Studies in Health Technology and Informatics

Abstract

One of the research projects running at the medical informatics department of the Institute of Computer Science AS CR explores the problem of medical information representation and development of electronic health record (EHR). With respect to this effort an interesting problem arises: how to transfer knowledge from a medical record written in a free text form into a structured electronic format represented by the EHR. Currently, this task was solved by writing extraction rules (regular expressions) for every element of information that is to be extracted from the medical record. However, such approach is very time consuming and requires supervision of a skilled programmer whenever the target area of medicine is changed. In this article we explore the possibility to mechanize this process by automatically generating the extraction rules from a pre-annotated corpus of medical records. Since we are currently in the phase of data acquisition and preliminary tests we will not present any final results, rather we will sketch the technologies we intend to use and describe the tools that were developed so far as a part of this project.