Development of feedforward receptive field structure of a simple cell and its contribution to the orientation selectivity: a modeling study

Recent experimental studies of hetero-synaptic interactions in various systems have shown the role of signaling in the plasticity, challenging the conventional understanding of Hebb's rule. It has also been found that activity plays a major role in plasticity, with neurotrophins acting as molecular signals translating activity into structural changes. Furthermore, role of synaptic efficacy in biasing the outcome of competition has also been revealed recently. Motivated by these experimental findings we present a model for the development of simple cell receptive field structure based on the competitive hetero-synaptic interactions for neurotrophins combined with cooperative hetero-synaptic interactions in the spatial domain. We find that with proper balance in competition and cooperation, the inputs from two populations (ON/OFF) of LGN cells segregate starting from the homogeneous state. We obtain segregated ON and OFF regions in simple cell receptive field. Our modeling study supports the experimental findings, suggesting the role of synaptic efficacy and the role of spatial signaling. We find that using this model we obtain simple cell RF, even for positively correlated activity of ON/OFF cells. We also compare different mechanism of finding the response of cortical cell and study their possible role in the sharpening of orientation selectivity. We find that degree of selectivity improvement in individual cells varies from case to case depending upon the structure of RF field and type of sharpening mechanism.     

A model for the coordinated development of columnar systems in primate striate cortex

The existence of patchy regions in primate striate cortex in which orientation selectivity is reduced, and which lie in the centers of ocular dominance stripes is well established (Hubel and Livingstone 1981). Analysis of functional maps obtained with voltage sensitive dyes (Blasdel and Salama 1986) has suggested that regions where the spatial rate of change of orientation preference is high, tend to be aligned either along the centers of ocular dominance stripes, or to intersect stripe borders at right angles. In this paper I present results from a developmental model which show that a tendency for orientation selectivity to develop more slowly in the centers of ocular dominance stripes would lead to the observed relationships between the layout of ocular dominance and the map of orientation gradient. This occurs despite the fact that there is no direct connection between the measures of preferred orientation (from which the gradient map is derived) and orientation selectivity (which is independent of preferred orientation). I also show that in both the monkey and the model, orientation singularities have an irregular distribution, but tend to be concentrated in the centers of the ocular dominance stripes. The average density of singularities is about 3/ ², where is the period of the orientation columns. The results are based on an elaboration of previous models (Swindale 1980, 1982) which show how, given initially disordered starting conditions, lateral interactions that are short-range excitatory and long-range inhibitory can lead to the development of patterns of orientation or ocular dominance that resemble those found in monkey striate cortex. To explain the coordinated development of the two kinds of column, it is proposed that there is an additional tendency in development for the rate of increase in orientation selectivity to be reduced in the centers of emerging ocular dominance stripes. This might come about if a single factor modulates plasticity in each cell, or column of cells. Thus plasticity may be turned off first in regions in the centers of ocular dominance stripes where relatively extreme and therefore stable ocular dominance values are achieved early in development. Consequently it will be hard for cells in these columns to modify other properties such as orientation preference or selectivity.     

Synaptic integration by V1 neurons depends on location within the orientation map

Neurons in the primary visual cortex (V1) are organized into an orientation map consisting of orientation domains arranged radially around "pinwheel centers" at which the representations of all orientations converge. We have combined optical imaging of intrinsic signals with intracellular recordings to estimate the subthreshold inputs and spike outputs of neurons located near pinwheel centers or in orientation domains. We find that neurons near pinwheel centers have subthreshold responses to all stimulus orientations but spike responses to only a narrow range of orientations. Across the map, the selectivity of inputs covaries with the selectivity of orientations in the local cortical network, while the selectivity of spike outputs does not. Thus, the input-output transformation performed by V1 neurons is powerfully influenced by the local structure of the orientation map.     

Spontaneous symmetry breaking in self–organizing neural fields

We extend the theory of self-organizing neural fields in order to analyze the joint emergence of topography and feature selectivity in primary visual cortex through spontaneous symmetry breaking. We first show how a binocular one-dimensional topographic map can undergo a pattern forming instability that breaks the underlying symmetry between left and right eyes. This leads to the spatial segregation of eye specific activity bumps consistent with the emergence of ocular dominance columns. We then show how a 2-dimensional isotropic topographic map can undergo a pattern forming instability that breaks the underlying rotation symmetry. This leads to the formation of elongated activity bumps consistent with the emergence of orientation preference columns. A particularly interesting property of the latter symmetry breaking mechanism is that the linear equations describing the growth of the orientation columns exhibits a rotational shift-twist symmetry, in which there is a coupling between orientation and topography. Such coupling has been found in experimentally generated orientation preference maps     

A reaction-diffusion model to capture disparity selectivity in primary visual cortex

Decades of experimental studies are available on disparity selective cells in visual cortex of macaque and cat. Recently, local disparity map for iso-orientation sites for near-vertical edge preference is reported in area 18 of cat visual cortex. No experiment is yet reported on complete disparity map in V1. Disparity map for layer IV in V1 can provide insight into how disparity selective complex cell receptive field is organized from simple cell subunits. Though substantial amounts of experimental data on disparity selective cells is available, no model on receptive field development of such cells or disparity map development exists in literature. We model disparity selectivity in layer IV of cat V1 using a reaction-diffusion two-eye paradigm. In this model, the wiring between LGN and cortical layer IV is determined by resource an LGN cell has for supporting connections to cortical cells and competition for target space in layer IV. While competing for target space, the same type of LGN cells, irrespective of whether it belongs to left-eye-specific or right-eye-specific LGN layer, cooperate with each other while trying to push off the other type. Our model captures realistic 2D disparity selective simple cell receptive fields, their response properties and disparity map along with orientation and ocular dominance maps. There is lack of correlation between ocular dominance and disparity selectivity at the cell population level. At the map level, disparity selectivity topography is not random but weakly clustered for similar preferred disparities. This is similar to the experimental result reported for macaque. The details of weakly clustered disparity selectivity map in V1 indicate two types of complex cell receptive field organization.     

The influence of cortical feature maps on the encoding of the orientation of a short line

The inhomogeneous distribution of the receptive fields of cortical neurons influences the cortical representation of the orientation of short lines seen in visual images. We construct a model of the response of populations of neurons in the human primary visual cortex by combining realistic response properties of individual neurons and cortical maps of orientation and location preferences. The encoding error, which characterizes the difference between the parameters of a visual stimulus and their cortical representation, is calculated using Fisher information as the square root of the variance of a statistically efficient estimator. The error of encoding orientation varies considerably with the location and orientation of the short line stimulus as modulated by the underlying orientation preference map. The average encoding error depends only weakly on the structure of the orientation preference map and is much smaller than the human error of estimating orientation measured psychophysically. From this comparison we conclude that the actual mechanism of orientation perception does not make efficient use of all the information available in the neuronal responses and that it is the decoding of visual information from neuronal responses that limits psychophysical performance. Action Editor: Terrence Sejnowski     

Receptive field properties of near neighbor orientation selective neurons in the visual cortex: a modeling study

The primary visual cortex is organized into clusters of cells having similar receptive fields (RFs). A purely feedforward model has been shown to produce realistic simple cell receptive fields. The modeled cells capture a wide range of receptive field properties of orientation selective cortical cells. We have analyzed the responses of 78 nearby cell pairs to study which RF properties are clustered. Orientation preference shows strongest clustering. Orientation tuning width (hwhh) and tuning height (spikes/sec) at the preferred orientation are not as tightly clustered. Spatial frequency is also not as tightly clustered and RF phase has the least clustering. Clustering property of orientation preference, orientation tuning height and width depend on the location of cells in the orientation map. No such location dependence is observed for spatial frequency and RF phase. Our results agree well with experimental data.     

Sequence of Genes Replicated in Salmonella typhimurium as Examined by Transduction Techniques

5-Bromouracil (BU) was pulsed into the genome of synchronously growing cells of an F(-) strain of Salmonella typhimurium of LT2. BU-labeled genes were transduced with P22 phage to a series of recipient auxotrophs. When BU was incorporated early in the replication cycle, the transducing markers that had hybrid densities were those that lie between 9 and 12 o'clock on the genetic map. When BU was incorporated during the terminal period of the synchronous cycle, the transducing particles had hybrid densities for genes that lie from 1 to 8 o'clock clockwise. When phages were prepared on cells in which the middle period was BU-labeled, transducing particles with hybrid densities appeared for genes that lie in two separate regions: between 7 and 9 o'clock and between 12 and 2 o'clock. Analysis of the map sequences of the transduced BU genes, the relative frequency of transduction for each marker, and the time sequence of replication led to a hypothesis that the origin of replication is near the isoleucine-valine gene on the chromosome map. As for direction of replication, several models were considered, including the possibility that replication may proceed in both directions in the same chromosome. It was also found that the aroB, cysG, and strA genes are cotransduced and strA and aroC are also cotransduced. The relative order of the four genes was found to be aroB, cysG, strA, aroC, but the orientation in the circular map was not determined.     

V1 neurons: in tune with the neighbors

In this issue of Neuron, Nauhaus et al. use a combination of optical imaging and multiple electrode recording to demonstrate that the orientation tuning of single cells in primary visual cortex is reliably related to the local structure of the orientation preference map in both cats and monkeys.     

The M1 and pre-M1 segments contribute differently to ion selectivity in ASICs and ENaCs

The ability to discriminate between different ionic species, termed ion selectivity, is a key feature of ion channels and forms the basis for their physiological function. Members of the degenerin/epithelial sodium channel (DEG/ENaC) superfamily of trimeric ion channels are typically sodium selective, but to a surprisingly variable degree. While acid-sensing ion channels (ASICs) are weakly sodium selective (sodium:potassium ratio ∼10:1), ENaCs show a remarkably high preference for sodium over potassium (>500:1). This discrepancy may be expected to originate from differences in the pore-lining second transmembrane segment (M2). However, these show a relatively high degree of sequence conservation between ASICs and ENaCs, and previous functional and structural studies could not unequivocally establish that differences in M2 alone can account for the disparate degrees of ion selectivity. By contrast, surprisingly little is known about the contributions of the first transmembrane segment (M1) and the preceding pre-M1 region. In this study, we used conventional and noncanonical amino acid–based mutagenesis in combination with a variety of electrophysiological approaches to show that the pre-M1 and M1 regions of mASIC1a channels are major determinants of ion selectivity. Mutational investigations of the corresponding regions in hENaC show that these regions contribute less to ion selectivity, despite affecting ion conductance. In conclusion, our work suggests that the remarkably different degrees of sodium selectivity in ASICs and ENaCs are achieved through different mechanisms. These results further highlight how M1 and pre-M1 are likely to differentially affect pore structure in these related channels.     

Modeling the role of gap junctions between excitatory neurons in the developing visual cortex

Recent experiments in the developing mammalian visual cortex have revealed that gap junctions couple excitatory cells and potentially influence the formation of chemical synapses. In particular, cells that were coupled by a gap junction during development tend to share an orientation preference and are preferentially coupled by a chemical synapse in the adult cortex, a property that is diminished when gap junctions are blocked. In this work, we construct a simplified model of the developing mouse visual cortex including spike-timing-dependent plasticity of both the feedforward synaptic inputs and recurrent cortical synapses. We use this model to show that synchrony among gap-junction-coupled cells underlies their preference to form strong recurrent synapses and develop similar orientation preference; this effect decreases with an increase in coupling density. Additionally, we demonstrate that gap-junction coupling works, together with the relative timing of synaptic development of the feedforward and recurrent synapses, to determine the resulting cortical map of orientation preference.     

Distribution of Orientation Selectivity in Recurrent Networks of Spiking Neurons with Different Random Topologies

Neurons in the primary visual cortex are more or less selective for the orientation of a light bar used for stimulation. A broad distribution of individual grades of orientation selectivity has in fact been reported in all species. A possible reason for emergence of broad distributions is the recurrent network within which the stimulus is being processed. Here we compute the distribution of orientation selectivity in randomly connected model networks that are equipped with different spatial patterns of connectivity. We show that, for a wide variety of connectivity patterns, a linear theory based on firing rates accurately approximates the outcome of direct numerical simulations of networks of spiking neurons. Distance dependent connectivity in networks with a more biologically realistic structure does not compromise our linear analysis, as long as the linearized dynamics, and hence the uniform asynchronous irregular activity state, remain stable. We conclude that linear mechanisms of stimulus processing are indeed responsible for the emergence of orientation selectivity and its distribution in recurrent networks with functionally heterogeneous synaptic connectivity.     

Coordinated Optimization of Visual Cortical Maps (I) Symmetry-based Analysis

In the primary visual cortex of primates and carnivores, functional architecture can be characterized by maps of various stimulus features such as orientation preference (OP), ocular dominance (OD), and spatial frequency. It is a long-standing question in theoretical neuroscience whether the observed maps should be interpreted as optima of a specific energy functional that summarizes the design principles of cortical functional architecture. A rigorous evaluation of this optimization hypothesis is particularly demanded by recent evidence that the functional architecture of orientation columns precisely follows species invariant quantitative laws. Because it would be desirable to infer the form of such an optimization principle from the biological data, the optimization approach to explain cortical functional architecture raises the following questions: i) What are the genuine ground states of candidate energy functionals and how can they be calculated with precision and rigor? ii) How do differences in candidate optimization principles impact on the predicted map structure and conversely what can be learned about a hypothetical underlying optimization principle from observations on map structure? iii) Is there a way to analyze the coordinated organization of cortical maps predicted by optimization principles in general? To answer these questions we developed a general dynamical systems approach to the combined optimization of visual cortical maps of OP and another scalar feature such as OD or spatial frequency preference. From basic symmetry assumptions we obtain a comprehensive phenomenological classification of possible inter-map coupling energies and examine representative examples. We show that each individual coupling energy leads to a different class of OP solutions with different correlations among the maps such that inferences about the optimization principle from map layout appear viable. We systematically assess whether quantitative laws resembling experimental observations can result from the coordinated optimization of orientation columns with other feature maps.     

Is the development of orientation selectivity instructed by activity?

Is the development of orientation selectivity in visual cortex instructed by the patterns of neural activity of input neurons? We review evidence as to the role of activity, review models of activity-instructed development, and discuss how these models can be tested. The models can explain the normal development of simple cells with binocularly matched orientation preferences, the effects of monocular deprivation and reverse suture on the orientation map, and the development of a full intracortical circuit sufficient to explain mature response properties including the contrast-invariance of orientation tuning. Existing experiments are consistent with the models, in that (a) selective blockade of ON-center ganglion cells, which will degrade or eliminate the information predicted to drive development of orientation selectivity, in fact prevents development of orientation selectivity; and (b) the spontaneous activities of inputs serving the two eyes are correlated in the lateral geniculate nucleus at appropriate developmental times, as was predicted to be required to achieve binocular matching of preferred orientations. However, definitive tests remain to be done to firmly establish the instructive rather than simply permissive role of activity and determine whether the retinotopically and center type-specific patterns of activity predicted by the models actually exist. We conclude by critically examining alternative scenarios for the development of orientation selectivity and maps, including the idea that maps are genetically prespecified.     

Low red cell production may protect against severe anemia during a malaria infection—Insights from modeling

The malaria parasite causes lysis of red blood cells, resulting in anemia, a major cause of mortality and morbidity. Intuitively, one would expect the production of red blood cells to increase in order to compensate for this loss. However, it has been observed that this response is weaker than would be expected. Furthermore, iron supplementation for iron deficient children in malaria endemic regions can paradoxically adversely affect the clinical outcome of malaria infection. A possible explanation may lie in the preference that some malaria parasites show for infecting immature red blood cells (reticulocytes). In the presence of a parasite preference for immature red cells, a rise in red cell production can ‘fuel the fire’ of infection by increasing the availability of the parasite's preferred target cell.We present a mathematical model of red blood cell production and infection in order to explore this hypothesis. We assess the effect of varying the reticulocyte replacement rate and preference of the parasite for reticulocytes on four key outcome measures assessing anemia and parasitemia.For a given level of parasite preference for reticulocytes we uncover an optimal erythropoietic response which minimizes disease severity. Increasing red blood cell production much above this optimum confers no benefit to the patient, and in fact can increase the degree of anemia and parasitemia. These conclusions are consistent with epidemiological studies demonstrating that both iron deficiency and anemia are protective against severe malaria, whilst iron supplementation in malaria endemic regions is with an increased number of malaria related adverse effects. Thus, suppression of red blood cell production, rather than being an unfortunate side effect of inflammation, may be a host protective effect against severe malarial anemia.     

Large-scale modeling of the primary visual cortex: influence of cortical architecture upon neuronal response.

A large-scale computational model of a local patch of input layer 4 [Formula: see text] of the primary visual cortex (V1) of the macaque monkey, together with a coarse-grained reduction of the model, are used to understand potential effects of cortical architecture upon neuronal performance. Both the large-scale point neuron model and its asymptotic reduction are described. The work focuses upon orientation preference and selectivity, and upon the spatial distribution of neuronal responses across the cortical layer. Emphasis is given to the role of cortical architecture (the geometry of synaptic connectivity, of the ordered and disordered structure of input feature maps, and of their interplay) as mechanisms underlying cortical responses within the model. Specifically: (i) Distinct characteristics of model neuronal responses (firing rates and orientation selectivity) as they depend upon the neuron's location within the cortical layer relative to the pinwheel centers of the map of orientation preference; (ii) A time independent (DC) elevation in cortico-cortical conductances within the model, in contrast to a "push-pull" antagonism between excitation and inhibition; (iii) The use of asymptotic analysis to unveil mechanisms which underly these performances of the model; (iv) A discussion of emerging experimental data. The work illustrates that large-scale scientific computation--coupled together with analytical reduction, mathematical analysis, and experimental data, can provide significant understanding and intuition about the possible mechanisms of cortical response. It also illustrates that the idealization which is a necessary part of theoretical modeling can outline in sharp relief the consequences of differing alternative interpretations and mechanisms--with final arbiter being a body of experimental evidence whose measurements address the consequences of these analyses.     

Footprinting studies of DNA-sequence recognition by nogalamycin.

We have studied the DNA sequence binding preference of the antitumour antibiotic nogalamycin by DNase-I footprinting using a variety of DNA fragments. The DNA fragments were obtained by cloning synthetic oligonucleotides into longer DNA fragments and were designed to contain isolated ligand-binding sites surrounded by repetitive sequences such as (A)n.(T)n and (AT)n. Within regions of (A)n.(T)n, clear footprints are observed with low concentrations of nogalamycin (< 5 microM), with apparent binding affinities for tetranucleotide sequences which decrease in the order TGCA > AGCT = ACGT > TCGA. In contrast, within regions of (AT)n, the ligand binds best to AGCT; binding to TCGA and TGCA is no stronger than to alternating AT. Within (ATT)n, the preference is for ACGT > TCGA. Although each of these binding sites contains all four base pairs, there is no apparent consensus sequence, suggesting that the selectivity is affected by local DNA dynamic and structural effects. At higher drug concentrations (> 25 microM), nogalamycin prevents DNAse-I cleavage of (AT)n but shows no interaction with regions of (AC)n.(GT)n. Regions of (A)n.(T)n, which are poorly cut by DNase I, show enhanced rates of cleavage in the presence of low concentrations of nogalamycin, but are protected from cleavage at higher concentrations. We suggest that this arises because drug binding to adjacent regions distorts the DNA to a structure which is more readily cut by the enzyme and which is better able to bind further ligand molecules.     

Application of Kohonen's self-organizing feature map algorithm to cortical maps of orientation and direction preference

Cortical maps of orientation preference in cats, ferrets and monkeys contain numerous half-rotation point singularities. Experimental data have shown that direction preference also has a smooth representation in these maps, with preferences being for the most part orthogonal to the axis of preferred orientation. As a result, the orientation singularities induce an extensive set of linear fractures in the direction map. These fractures run between and connect nearby point orientation singularities. Their existence appears to pose a puzzle for theories that postulate that cortical maps maximize continuity of representation, because the fractures could be avoided if the orientation map contained full-rotation singularities. Here we show that a dimension-reduction model of cortical map formation, which implements principles of continuity and completeness, produces an arrangement of linear direction fractures connecting point orientation singularities which is similar to that observed experimentally. We analyse the behaviour of this model and suggest reasons why the model maps contain half-rotation rather than full-rotation orientation singularities.