This paper is concerned with a striking visual experience: that of seeing geometric visual hallucinations. Hallucinatory images were classi¢ed by Klu« ver into four groups called form constants comprising (i) gratings, lattices, fretworks, ¢ligrees, honeycombs and chequer-boards, (ii) cobwebs, (iii) tunnels, funnels, alleys, cones and vessels, and (iv) spirals. This paper describes a mathematical investigation of their origin based on the assumption that the patterns of connection between retina and striate cortex (henceforth referred to as V1)öthe retinocortical mapöand of neuronal circuits in V1, both local and lateral, determine their geometry.In the ¢rst part of the paper we show that form constants, when viewed in V1 coordinates, essentially correspond to combinations of plane waves, the wavelengths of which are integral multiples of the width of a human Hubel^Wiesel hypercolumn, ca. 1.33^2 mm. We next introduce a mathematical description of the large-scale dynamics of V1 in terms of the continuum limit of a lattice of interconnected hypercolumns, each of which itself comprises a number of interconnected iso-orientation columns. We then show that the patterns of interconnection in V1 exhibit a very interesting symmetry, i.e. they are invariant under the action of the planar Euclidean group E(2)öthe group of rigid motions in the planeö rotations, re£ections and translations. What is novel is that the lateral connectivity of V1 is such that a new group action is needed to represent its properties: by virtue of its anisotropy it is invariant with respect to certain shifts and twists of the plane. It is this shift^twist invariance that generates new representations of E(2). Assuming that the strength of lateral connections is weak compared with that of local connections, we next calculate the eigenvalues and eigenfunctions of the cortical dynamics, using Rayleigh^Schro« dinger perturbation theory. The result is that in the absence of lateral connections, the eigenfunctions are degenerate, comprising both even and odd combinations of sinusoids in , the cortical label for orientation preference, and plane waves in r, the cortical position coordinate.`Switching-on' the lateral interactions breaks the degeneracy and either even or else odd eigenfunctions are selected. These results can be shown to follow directly from the Euclidean symmetry we have imposed.In the second part of the paper we study the nature of various even and odd combinations of eigenfunctions or planforms, the symmetries of which are such that they remain invariant under the particular action of E(2) we have imposed. These symmetries correspond to certain subgroups of E(2), the so-called axial subgroups. Axial subgroups are important in that the equivariant branching lemma indicates that when a symmetrical dynamical system becomes unstable, new solutions emerge which have symmetries corresponding to the axial subgroups of the underlying symmetry group. This is precisely the case studied in this paper. Thus we study the various planforms that emer...

Efficiently identifying and quantifying disease- or treatment-related changes in the abundance of proteins is an important area of research for the pharmaceutical industry. Here we describe an automated, label-free method for finding differences in complex mixtures using complete LC-MS data sets, rather than subsets of extracted peaks or features. The method selectively finds statistically significant differences in the intensity of both high-abundance and low-abundance ions, accounting for the variability of measured intensities and the fact that true differences will persist in time. The method was used to compare two complex peptide mixtures with known peptide differences. This controlled experiment allowed us to assess the validity of each difference found and so to analyze the method's sensitivity and specificity. The method detects both presence versus absence and a 2-fold change in peptide concentration near the limit of detection of the instrument used, where chromatographic peaks may not be sufficiently well defined to be detected in individual samples. The method is more sensitive and gives fewer false positives than subtractive methods that ignore signal variability. Differential mass spectrometry combined with targeted MS/MS analysis of only identified differences may save both computation time and human effort compared to shotgun proteomics approaches.

To dissect common human diseases such as obesity and diabetes, a systematic approach is needed to study how genes interact with one another, and with genetic and environmental factors, to determine clinical end points or disease phenotypes. Bayesian networks provide a convenient framework for extracting relationships from noisy data and are frequently applied to large-scale data to derive causal relationships among variables of interest. Given the complexity of molecular networks underlying common human disease traits, and the fact that biological networks can change depending on environmental conditions and genetic factors, large datasets, generally involving multiple perturbations (experiments), are required to reconstruct and reliably extract information from these networks. With limited resources, the balance of coverage of multiple perturbations and multiple subjects in a single perturbation needs to be considered in the experimental design. Increasing the number of experiments, or the number of subjects in an experiment, is an expensive and time-consuming way to improve network reconstruction. Integrating multiple types of data from existing subjects might be more efficient. For example, it has recently been demonstrated that combining genotypic and gene expression data in a segregating population leads to improved network reconstruction, which in turn may lead to better predictions of the effects of experimental perturbations on any given gene. Here we simulate data based on networks reconstructed from biological data collected in a segregating mouse population and quantify the improvement in network reconstruction achieved using genotypic and gene expression data, compared with reconstruction using gene expression data alone. We demonstrate that networks reconstructed using the combined genotypic and gene expression data achieve a level of reconstruction accuracy that exceeds networks reconstructed from expression data alone, and that fewer subjects may be required to achieve this superior reconstruction accuracy. We conclude that this integrative genomics approach to reconstructing networks not only leads to more predictive network models, but also may save time and money by decreasing the amount of data that must be generated under any given condition of interest to construct predictive network models.

Many observers see geometric visual hallucinations after taking hallucinogens such as LSD, cannabis, mescaline or psilocybin; on viewing bright flickering lights; on waking up or falling asleep; in "near-death" experiences; and in many other syndromes. Klüver organized the images into four groups called form constants: (I) tunnels and funnels, (II) spirals, (III) lattices, including honeycombs and triangles, and (IV) cobwebs. In most cases, the images are seen in both eyes and move with them. We interpret this to mean that they are generated in the brain. Here, we sum- nectivity between hypercolumns exhibits symmetries, rendering it invariant under the action of the Euclidean group E(2), composed of reflections and translations in the plane, and a (novel) shift-twist action. Using this symmetry, we show that the various patterns of activity that spontaneously emerge when V1's spatially uniform resting state becomes unstable correspond to the form constants when transformed to the visual field using the retino-cortical map. The results are sensitive to the detailed specification of the lateral connectivity and suggest that the cortical mechanisms that generate geometric visual hallucinations are closely related to those used to process edges, contours, surfaces, and textures.

It is not clear how information related to cognitive or psychological processes is carried by or represented in the responses of single neurons. One provocative proposal is that precisely timed spike patterns play a role in carrying such information. This would require that these spike patterns have the potential for carrying information that would not be available from other measures such as spike count or latency. We examined exactly timed (1-ms precision) triplets and quadruplets of spikes in the stimulus-elicited responses of lateral geniculate nucleus (LGN) and primary visual cortex (V1) neurons of the awake fixating rhesus monkey. Large numbers of these precisely timed spike patterns were found. Information theoretical analysis showed that the precisely timed spike patterns carried only information already available from spike count, suggesting that the number of precisely timed spike patterns was related to firing rate. We therefore examined statistical models relating precisely timed spike patterns to response strength. Previous statistical models use observed properties of neuronal responses such as the peristimulus time histogram, interspike interval, and/or spike count distributions to constrain the parameters of the model. We examined a new stochastic model, which unlike previous models included all three of these constraints and unlike previous models predicted the numbers and types of observed precisely timed spike patterns. This shows that the precise temporal structures of stimulus-elicited responses in LGN and V1 can occur by chance. We show that any deviation of the spike count distribution, no matter how small, from a Poisson distribution necessarily changes the number of precisely timed spike patterns expected in neural responses. Overall the results indicate that the fine temporal structure of responses can only be interpreted once all the coarse temporal statistics of neural responses have been taken into account.