Statistical structure of lateral connections in the primary visual cortex
© Hunt et al; licensee BioMed Central Ltd. 2011
Received: 12 April 2010
Accepted: 16 September 2010
Published: 26 January 2011
The statistical structure of the visual world offers many useful clues for understanding how biological visual systems may understand natural scenes. One particularly important early process in visual object recognition is that of grouping together edges which belong to the same contour. The layout of edges in natural scenes have strong statistical structure. One such statistical property is that edges tend to lie on a common circle, and this 'co-circularity' can predict human performance at contour grouping. We therefore tested the hypothesis that long-range excitatory lateral connections in the primary visual cortex, which are believed to be involved in contour grouping, display a similar co-circular structure.
By analyzing data from tree shrews, where information on both lateral connectivity and the overall structure of the orientation map was available, we found a surprising diversity in the relevant statistical structure of the connections. In particular, the extent to which co-circularity was displayed varied significantly.
Overall, these data suggest the intriguing possibility that V1 may contain both co-circular and anti-cocircular connections.
Grouping edges which belong to the same object or contour is a vital part of object recognition. Biological vision systems excel at this task, yet it is still extremely challenging for artificial vision systems . In humans, contour detection develops well after birth , suggesting that learning from visual experience plays an important role.
The layout of edges in natural scenes has strong statistical structure. One well-known property is co-linearity. Co-linearity is important in contour grouping and texture discrimination tasks [3–5]. A more recently discovered generalization of co-linearity displayed by natural scenes is co-circularity: edges tend to be tangent to a common circle more often than would be expected by chance . Furthermore, the degree of co-circularity in a contour can be used to predict human contour detection performance . This raises the question of what biological substrate underlies this effect of co-circularity on performance.
An obvious candidate is the pattern of excitatory lateral connections in the primary visual cortex (V1). These connections are primarily found in layer 2/3 and are longer range than the inhibitory connections found predominately in other layers . In humans these connections develop after birth, and their development coincides with improvement in contour detection [2, 9, 10]. Furthermore they have been shown to play a role in co-linear facilitation [11–14].
It has been previously demonstrated that long-range lateral connections connect preferentially to patches of similar orientation and primarily along the axis of orientation in V1 [15–18] (although others have claimed this preference is weak ). This is consistent with their role in co-linear facilitation. Although previous work has demonstrated that the variance in some of these data is not adequately explained by co-linearity alone , whether there is a more general tendency towards co-circularity in these connections has not been investigated.
Here we reanalyzed previously published functional anatomical data regarding lateral connections in tree-shrew , using a noise-resistant measure of co-circularity that we recently introduced to study the structure of orientation maps . Surprisingly, we found a large variation in the statistics of lateral connections between animals. One unifying explanation for both this and previous results is that lateral connections in V1 are connected in a variety of ways, both co-circularly and anti-cocircularly. We suggest reasons why this may be a desirable arrangement for analyzing the structure of visual scenes.
We reanalyzed data from four tree shrews previously published by Bosking et al. . The experimental methods used are described in detail in that paper and we provide only a brief overview here. Tree shrew orientation preference maps were obtained using optical imaging. Additionally, 540 nm light was used to map surface blood vessels used for alignment. Biocytin was then injected into a specific site in V1 and the animal was sacrificed 16 hours later. Slices of V1 were imaged to locate the biocytin bouton and the surface blood vessels. The blood vessel information was then used to align the orientation preference maps with the bouton images giving overlaid information on the underlying connectivity from the injection site on the animal. The original experiment used a total of ten cases, however, we were only able to recover the data for four cases.
For the original analysis  electrophysiological recordings taken during the injection were used to determine orientation. We determined the orientation of the injection site using the orientation value of the underlying pixels that were measured using optical imaging. This differed from the electrophysiological recordings by an average of 15° (range 0-23°) from the optical orientation value (averaged across the animals).
Characterizing the lateral connection statistics
The orientation preferences and are defined in real space where the stimulus is presented, relative to zero degrees along a horizontal x-axis. They are always taken on the y-positive side of the horizontal and thus range from 0° to 180°. The relative position of the cortical points ϕ is defined relative to the horizontal axis of the retinotopic map and is calculated modulo 180°. Except where noted, calculations of orientations and angular positions were done using modulo 180° arithmetic.
The wavelength λ of each map was calculated as the mean of the Fourier spectrum averaged over all directions as previously described . As there is variation in column spacing between individuals , and the imaging of each map may be at different magnifications, the use of λ as the unit of length aims to ensure that the quantification of reduced symmetry is not affected by overall changes in map scaling. The wavelength gives a rough measure of the size of iso-orientation patches in the map.
We calculated a probability distribution for the connection statistics for each animal p(α, β, r). However, if ϕ is not uniformly distributed this can introduce spurious correlations. We corrected for this by calculating a four-dimensional probability p(, β, r, ϕ) and then marginalized along ϕ. This was done by discretizing both distance r and orientations , β and ϕ and binning the connections. Given the full distribution it was then possible to calculate the relevant marginal distributions such as p(|β, r) and p(β|, r).
Ddiff(r) has the intuitively appealing property that it is 45° when no co-circularity is present (that is, when the connections from the injection site have no systematic dependence with orientation and angle) and decreases if the connections are connected co-circularly. To calculate Ddiff(r) we binned r in intervals of 1.5λ. The first bin, referred to hereafter as r = 0, contained pairs where 0 ≤ r <0.75λ, the r = 1.5λ bin contained the pairs where 0.75λ ≤ r <2.25λ and so on. This bin spacing balances the requirements of having a reasonable number of pairs in each bin to reduce noise, while allowing fine enough resolution to investigate trends with r.
In addition, the distribution of co-circularity from each injection/bouton pair can provide useful information since various distributions of ddiff, the co-circularity of individual bouton sites (see eq. 2), would result in the same mean Ddiff(r). In addition to calculating the average co-circularity Ddiff(r), the values of ddiff = | + − 2ϕ| mod 180° for each pair at distance r were therefore binned to create a histogram H (r, d). The r bins were the same as used for calculating Ddiff(r) and the ddiff bins were spaced 15° apart. Since we were not interested in the absolute number of points at different distances, we calculated and plotted the conditional probability P(d|r). This allows a basic characterization of the distribution of co-circularity values.
Co-circularity with offset
for = 0 to 80° in 10° increments. Since co-circularity is invariant under ϕ → ϕ + 90° searching over this range covers all possible values for . We also calculated (r) for each value of r, which was defined as the value of which minimized Ddiff(r, ).
Testing for co-circularity
In the original analysis of the tree shrew data it was demonstrated that lateral connections are denser along the axis of injection site orientation and connect preferentially to sites of similar orientation. This evidence demonstrates that lateral connections are likely to prefer co-linearity. In order to test whether the co-circularity we found in the horizontal connections was simply due to this co-linearity, we recalculated Ddiff(r) (eq. 4), but excluded all points lying near the axis of the injection site. This was done by excluding all injection/bouton pairs with β <30° or β >150°.
Controls and statistics
It was important to have a reliable control to assess the significance of any co-circularity measured in the maps. We are interested in the significance of any deviation from 45° that occurs in Ddiff(r). One comparison was simply examining the values of Ddiff(r, ) at ≠ 0 since these are expected by our hypothesis to be less co-circular then with the true origin ( = 0). This means that we would expect to be near 0 when connections are the most co-circular.
We created 99 control cases (N = 99) for each animal by adding Gaussian noise (σ = 100 pixels) to the bouton positions and calculated Ddiff(r) j for each of these controls (j denotes the index of the control). The added noise is high enough to overwhelm any co-circular preferences in the connections since it is much larger than the wavelength of the maps, which means that it is larger than the size of iso-orientation patches in the maps.
where H (.) is the Heaviside function (H (x) = 0 for x <0, H (x) = 1 for x ≥ 0). We considered P <0.05 to be significant.
We used an existing image dataset  for all natural scene analysis. These images are uncompressed, deblurred and have linearized intensity values and thus should have no artifacts due to compression or other steps in the image acquisition. This image set was also used when co-circularity in natural scenes was originally demonstrated . The images are 12-bit greyscale with a resolution of 1536 × 1024 pixels and an angular resolution of approximately 1 min of arc per pixel. The orientation and orientation strength of each pixel were calculated using steerable filters provided by matlabPyrTools . After extracting orientation and orientation strength we calculated Ddiff(r) for each image analogously to how it was calculated for lateral connections. The exception was that when averaging ddiff we weighted by the orientation strength of the pixels in the pair and we excluded pixels with very low orientation strength (essentially un-orientated pixels). This was done differently because we did not have information on orientation selectivity for most the animal data. We then calculated a natural scene value by re-averaging the values of Ddiff(r) for each scene to create an average across all scenes. We did this double averaging to ensure that the calculation was not overly preferenced towards particular types of scenes which have more orientation density (such as tree scenes).
Alignment and topography
Figure 2 shows the aligned orientation preference maps with overlaid bouton and injection sites for all animals. In each case the vertical meridian was selected using previously published results of tree shrew topography  and the map rotated to make it vertical. All pixels representing receptive field positions in the ipsilateral visual field (that is, those to the right of the vertical meridian) were removed and these sections were not used in subsequent analysis due to the very different magnification factor in this region of the cortex. For analysis, the remaining region of V1 was treated as isotropic with a constant magnification factor. We have previously demonstrated that this approximation does not significantly affect our measures of co-circularity .
Number of included boutons for each animal.
Co-circularity with offset is a generalization of co-circularity. Edges are co-circular when tangent to a common circle and co-circularity with offset is when edges occur at an orientation to the tangent of a common circle (so = 0 is co-circularity). Co-circularity with non-zero is not found in natural scenes and thus we would predict that if co-circularity in the lateral connections is significant, then it should be present more strongly than co-circularity with any non-zero offset . We measured co-circularity with offset at a range of non-zero offsets to test this hypothesis (Figure 4).
Case ts9509 is strongly anti-cocircular at longer ranges (Figure 4) and this tendency remains even if co-linear connections are removed (Figure 5). Case ts9514 is strongly co-circular, although when co-linearity is excluded this tendency becomes less strong (and it is anti-cocircular at large r). The two other cases show inconsistent results at different wavelengths, although case ts9531 shows strong co-circularity for the longest connections.
Examining the results in the four cases at different values of confirms that there is significant variation in the statistics of wiring between the four animals. Case 9509 has values of (r) which are near 45°, while case 9514 shows (r) near 0° (the other two cases show inconclusive results). These results demonstrate that the variation between co-circularity and anti-cocircularity between the cases is significant and are an additional verification that the wiring on these connections is not simply co-circular as our initial hypothesis predicted. As we discuss later, these findings indicate that the results cannot be explained by minor variations in topographical alignment between the animals.
Additionally, we directly examined the probability distribution of edges in both natural scenes and lateral connections (data not shown). We found connections are isotropic for short-range connections (as also seen in the Ddiff results). Long-range connections have strong biases for both orientation and direction. However, these statistics are probably due primarily to the previously studied elongation along the axis of the injection site orientation , rather than any specific arrangement of edges, as they are also present in controls (where noise is added to the bouton positions).
We examined lateral connections in V1 and compared the statistics of their connections with the edge arrangement of natural scenes. These connections develop after eye-opening [2, 9] and therefore may be driven by natural scene input. As found previously, at short distances the lateral connectivity is isotropic (and probably mostly inhibitory ). However, long-range lateral connections are clearly anisotropic. We hypothesised that these connections might display similar connectivity statistics to edges in natural scenes. However, we found that the connectivity statistics showed a range of behaviour for the different individuals and both co-circularity and anti-cocircularity were present. This study is, as far as we know, the first to directly examine lateral connectivity in V1 for evidence of co-circular connections.
Some limitations of our results are as follows. Firstly, due to archiving problems, the original data from the study of  was only available in fragmentary form and we were only able to reconstruct sufficient information for our analysis from four tree shrews. Secondly, the experimental data has a significant amount of noise due to error in alignment of optical imaging with histology. Because co-circularity is a second order property these errors are compounded when calculating co-circularity statistics. An additional possible source of noise is based on more theoretical considerations and inherent to any experimental study of these connections. Suppose the brain has some desired function for specifying lateral connectivity (for instance some spread of co-circularity), it is unlikely this function is realised without a significant level of noise in the connectivity. Because of this, even without any experimental uncertainty, it may be difficult to recover a good approximation to the original generating function. Assuming the noise in connectivity is independent between coloumns, this can be overcome by increasing the number of injection sites, however, in this study only four sites could be used. The restricted captive environment in which the tree shrews were raised was probably not as rich as a truly natural environment, however, at the level of statistical characterisation considered here these differences are not likely to be consequential. Because of these limitations, it is difficult to make definitive conclusions from our results. However, this data is still able to provide some important insights into functional connectivity in V1.
The variations between the four animals are not explainable by minor variations in topography. We examined the value of co-circularity for each of the animals for various values of , which is the value of co-circularity that would be present if the assigned topography was rotated by . In order for, for instance, case ts9509 to be in agreement with case ts9514 they would have to contain a combined error of 45° in the assigned orientation of the topography. Other work has demonstrated that the topography of tree shrews has little inter-individual variability and consistent local structure . Additionally, our previous work has shown that the co-circularity statistics we used are robust to other deviations in topography such as anisotropy.
Previous work  has shown that co-linearity alone does not fully explain the anisotropy in the original paper of Bosking et al.  and some broader form of connectivity is required. Here, we postulated that co-circularity similar to that found in natural scenes would explain this. However, we found that there is a large variation in co-circularity statistics between different injections sites (and animals). In particular case ts9509 shows strong anti-cocircularity, while case ts9514 shows strong co-circularity.
Several other studies have examined horizontal connections in V1 in a variety of animals. Early work in tree shrew  found long-range (around 2.5 mm) lateral connectivity and some indications of orientation selectivity in the connections, as did similar work in cats . There is conflicting information on whether connections in V1 are strongly iso-linear  (tree shrew) or mostly isotropic  (cat). However, there is significant evidence from psychophysics [11, 12] that lateral connections play a role in some forms of co-linear facilitation. It has also been demonstrated that connections are refined during development  and there are correlations in co-linear neuron spiking . In summary, there is strong evidence for co-linear lateral connections, but also indications that other types of connections are present.
Previous work has shown that an important functional role of lateral connections is inhibitory, such as iso-orientation suppression . However, the majority of connections we have considered here, particularly at longer distances (r >500 μm), are believed to be excitatory [8, 17]. Excitatory connections are known to be connected co-linearly and to be important for co-linear facilitation [11–13]. There is some evidence that V1 cells in some species may also be involved in border-ownership computation , however, since we are not aware of any characterisation of the type of connectivity necessary for these computations, we did not test the proposition here.
We suggest that there may be a role for both anti-cocircular and co-circular lateral connections. Edges in natural scenes are preferentially co-circular (with a small but consistent bias). For tasks such as co-linear facilitation , excitatory co-linear lateral connections are necessary. This sort of facilitation is thought to be involved in contour integration . However, from another point of view, it is precisely when input conflicts with what is expected, (in this case co-circularity), that it is providing novel and important information (indeed, it has been shown that high entropy areas of visual scenes are usually at the center of gaze ). From this point of view edges that are anti-cocircular are the most salient, and excitatory connections to facilitate their detection are a reasonable postulate. It is possible that the brain uses a combination of these strategies to both detect the unexpected while facilitating description of the more expected case (co-circularity). Our findings are consistent with such a bifurcated connection strategy.
An interesting experimental followup to this work would be to examine the lateral connectivity of animals reared under unusual rearing conditions (as done for orientation maps in ). Since these connections develop after eye opening it is likely that these rearing conditions would have an effect on lateral connectivity (indeed previous work has indicated that strabismus during development has an effect on connectivity ). This could provide further insight into the developmental mechanisms and statistics underlying lateral connections and the degree to which these connections are influenced by natural visual scenes during development.
We thank Len White and David Fitzpatrick for assistance and helpful discussions. This research was supported by grants from Australian Research Council (DP0878939) and the National Health and Medical Research Council (Project grant 525459).
- Chellapilla K, Larson K, Simard P, Czerwinski M: Building segmentation based human-friendly human interactive proofs (HIPs). Proceedings of the Second International Workshop on Human Interactive Proofs. 2005, 2: 1-26. full_text.View ArticleGoogle Scholar
- Kovács I, Kozma P, Fehér A, Benedek G: Late maturation of visual spatial integration in humans. Proc Natl Acad Sci USA. 1999, 96 (21): 12204-12209. 10.1073/pnas.96.21.12204.PubMed CentralView ArticlePubMedGoogle Scholar
- Schmidt KE, Goebel R, Löwel S, Singer W: The perceptual grouping criterion of colinearity is reflected by anisotropies of connections in the primary visual cortex. Eur J Neurosci. 1997, 9 (5): 1083-1089. 10.1111/j.1460-9568.1997.tb01459.x.View ArticlePubMedGoogle Scholar
- Roelfsema PR: Cortical algorithms for perceptual grouping. Annu Rev Neurosci. 2006, 29: 203-227. 10.1146/annurev.neuro.29.051605.112939.View ArticlePubMedGoogle Scholar
- Schwarzkopf D, Kourtzi Z: Experience shapes the utility of natural statistics for perceptual contour integration. Curr Biol. 2008, 18: 1162-1167. 10.1016/j.cub.2008.06.072.View ArticlePubMedGoogle Scholar
- Sigman M, Cecchi GA, Gilbert CD, Magnasco MO: On a common circle: natural scenes and Gestalt rules. Proc Natl Acad Sci USA. 2001, 98 (4): 1935-1940. 10.1073/pnas.031571498.PubMed CentralView ArticlePubMedGoogle Scholar
- Geisler WS, Perry JS, Super BJ, Gallogly DP: Edge co-occurrence in natural images predicts contour grouping performance. Vis Res. 2001, 41 (6): 711-724. 10.1016/S0042-6989(00)00277-7.View ArticlePubMedGoogle Scholar
- Kisvárday ZF, Tóth E, Rausch M, Eysel UT: Orientation-specific relationship between populations of excitatory and inhibitory lateral connections in the visual cortex of the cat. Cereb Cortex. 1997, 7 (7): 605-618. 10.1093/cercor/7.7.605.View ArticlePubMedGoogle Scholar
- Burkhalter A, Bernardo K, Charles V: Development of local circuits in human visual cortex. J Neurosci. 1993, 13 (5): 1916-PubMedGoogle Scholar
- Hou C, Pettet MW, Sampath V, Candy TR, Norcia AM: Development of the spatial organization and dynamics of lateral interactions in the human visual system. J Neurosci. 2003, 23 (25): 8630-8640.PubMedGoogle Scholar
- Polat U, Sagi D: Lateral interactions between spatial channels: suppression and facilitation revealed by lateral masking experiments. Vis Res. 1993, 33 (7): 993-999. 10.1016/0042-6989(93)90081-7.View ArticlePubMedGoogle Scholar
- Cass JR, Spehar B: Dynamics of collinear contrast facilitation are consistent with long-range horizontal striate transmission. Vis Res. 2005, 45 (21): 2728-2739. 10.1016/j.visres.2005.03.010.View ArticlePubMedGoogle Scholar
- Cass J, Alais D: The mechanisms of collinear integration. J Vis. 2006, 6 (9): 915-922. 10.1167/6.9.5.View ArticlePubMedGoogle Scholar
- Field D, Hayes A: Contour integration and the lateral connections of V1 neurons. The Visual Neurosciences. Edited by: Chalupa LM, Werner JS. 2004, Cambridge, MA: MIT Press, 1-9.Google Scholar
- Mitchison G, Crick F: Long axons within the striate cortex: their distribution, orientation, and patterns of connection. Proc Natl Acad Sci USA. 1982, 79 (11): 3661-3665. 10.1073/pnas.79.11.3661.PubMed CentralView ArticlePubMedGoogle Scholar
- Malach R, Amir Y, Harel M, Grinvald A: Relationship between intrinsic connections and functional architecture revealed by optical imaging and in vivo targeted biocytin injections in primate striate cortex. Proc Nat Acad Sci USA. 1993, 90 (22): 10469-10473. 10.1073/pnas.90.22.10469.PubMed CentralView ArticlePubMedGoogle Scholar
- Bosking WH, Zhang Y, Schofield B, Fitzpatrick D: Orientation selectivity and the arrangement of horizontal connections in tree shrew striate cortex. J Neurosci. 1997, 17 (6): 2112-2127.PubMedGoogle Scholar
- Ts'o DY, Gilbert CD, Wiesel TN: Relationships between horizontal interactions and functional architecture in cat striate cortex as revealed by cross-correlation analysis. J Neurosci. 1986, 6 (4): 1160-1170.PubMedGoogle Scholar
- Ben-Shahar O, Zucker S: Geometrical computations explain projection patterns of long-range horizontal connections in visual cortex. Neural Comput. 2004, 16 (3): 445-476. 10.1162/089976604772744866.View ArticlePubMedGoogle Scholar
- Hunt JJ, Giacomantonio CE, Tang H, Mortimer D, Jaffer S, Vorobyov V, Ericksson G, Sengpiel F, Goodhill GJ: Natural scene statistics and the structure of orientation maps in the visual cortex. Neuroimage. 2009, 47: 157-172. 10.1016/j.neuroimage.2009.03.052.PubMed CentralView ArticlePubMedGoogle Scholar
- Bosking WH, Kretz R, Pucak ML, Fitzpatrick D: Functional specificity of callosal connections in tree shrew striate cortex. J Neurosci. 2000, 20 (6): 2346-2359.PubMedGoogle Scholar
- Carreira-Perpin MA, Goodhill GJ: Influence of lateral connections on the structure of cortical maps. J Neurophysiol. 2004, 92 (5): 2947-2959. 10.1152/jn.00281.2004.View ArticleGoogle Scholar
- Kaschube M, Wolf F, Geisel T, Lowel S: Genetic influence of quantitative features of neocortical architecture. J Neurosci. 2002, 22: 7206-7217.PubMedGoogle Scholar
- van Hateren JH, van der Schaaf A: Independent component filters of natural images compared with simple cells in primary visual cortex. Proc Biol Sci. 1998, 265 (1394): 359-366. 10.1098/rspb.1998.0303.PubMed CentralView ArticlePubMedGoogle Scholar
- Simoncelli E, Freeman W: The steerable pyramid: a flexible architecture for multi-scale derivative computation. 2nd Annual Intl. Conf. on Image Processing. 1995Google Scholar
- Bosking WH, Crowley JC, Fitzpatrick D: Spatial coding of position and orientation in primary visual cortex. Nat Neurosci. 2002, 5 (9): 874-882. 10.1038/nn908.View ArticlePubMedGoogle Scholar
- Rockland KS, Lund JS: Widespread periodic intrinsic connections in the tree shrew visual cortex. Science. 1982, 215 (4539): 1532-1534. 10.1126/science.7063863.View ArticlePubMedGoogle Scholar
- Gilbert CD, Wiesel TN: Columnar specificity of intrinsic horizontal and corticocortical connections in cat visual cortex. J Neurosci. 1989, 9 (7): 2432-2442.PubMedGoogle Scholar
- Li Z: A saliency map in primary visual cortex. Trends Cogn Sci. 2002, 6: 9-16. 10.1016/S1364-6613(00)01817-9.View ArticlePubMedGoogle Scholar
- Tamura H, Sato H, Katsuyama N, Hata Y, Tsumoto T: Less segregated processing of visual information in V2 than in V1 of the monkey visual cortex. Eur J Neurosci. 1996, 8 (2): 300-309. 10.1111/j.1460-9568.1996.tb01214.x.View ArticlePubMedGoogle Scholar
- Reinagel P, Zador AM: Natural scene statistics at the centre of gaze. Network. 1999, 10 (4): 341-350. 10.1088/0954-898X/10/4/304.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.