The folding fingerprint of visual cortex reveals the timing of human V1 and V2

Dynamics of retinotopic spatial attention revealed by multifocal MEG

Visual focal attention is both fast and spatially localized, making it challenging to investigate using human neuroimaging paradigms. Here, we used a new multivariate multifocal mapping method with magnetoencephalography (MEG) to study how focal attention in visual space changes stimulus-evoked responses across the visual field. The observer's task was to detect a color change in the target location, or at the central fixation. Simultaneously, 24 regions in visual space were stimulated in parallel using an orthogonal, multifocal mapping stimulus sequence. First, we used univariate analysis to estimate stimulus-evoked responses in each channel. Then we applied multivariate pattern analysis to look for attentional effects on the responses. We found that attention to a target location causes two spatially and temporally separate effects. Initially, attentional modulation is brief, observed at around 60-130 ms post stimulus, and modulates responses not only at the target location but also in adjacent regions. A later modulation was observed from around 200 ms, which was specific to the location of the attentional target. The results support the idea that focal attention employs several processing stages and suggest that early attentional modulation is less spatially specific than late.

A visual encoding model links magnetoencephalography signals to neural synchrony in human cortex

Synchronization of neuronal responses over large distances is hypothesized to be important for many cortical functions. However, no straightforward methods exist to estimate synchrony non-invasively in the living human brain. MEG and EEG measure the whole brain, but the sensors pool over large, overlapping cortical regions, obscuring the underlying neural synchrony. Here, we developed a model from stimulus to cortex to MEG sensors to disentangle neural synchrony from spatial pooling of the instrument. We find that synchrony across cortex has a surprisingly large and systematic effect on predicted MEG spatial topography. We then conducted visual MEG experiments and separated responses into stimulus-locked and broadband components. The stimulus-locked topography was similar to model predictions assuming synchronous neural sources, whereas the broadband topography was similar to model predictions assuming asynchronous sources. We infer that visual stimulation elicits two distinct types of neural responses, one highly synchronous and one largely asynchronous across cortex.

A population receptive field model of the magnetoencephalography response

Computational models which predict the neurophysiological response from experimental stimuli have played an important role in human neuroimaging. One type of computational model, the population receptive field (pRF), has been used to describe cortical responses at the millimeter scale using functional magnetic resonance imaging (fMRI) and electrocorticography (ECoG). However, pRF models are not widely used for non-invasive electromagnetic field measurements (EEG/MEG), because individual sensors pool responses originating from several centimeter of cortex, containing neural populations with widely varying spatial tuning. Here, we introduce a forward-modeling approach in which pRFs estimated from fMRI data are used to predict MEG sensor responses. Subjects viewed contrast-reversing bar stimuli sweeping across the visual field in separate fMRI and MEG sessions. Individual subject's pRFs were modeled on the cortical surface at the millimeter scale using the fMRI data. We then predicted cortical time series and projected these predictions to MEG sensors using a biophysical MEG forward model, accounting for the pooling across cortex. We compared the predicted MEG responses to observed visually evoked steady-state responses measured in the MEG session. We found that pRF parameters estimated by fMRI could explain a substantial fraction of the variance in steady-state MEG sensor responses (up to 60% in individual sensors). Control analyses in which we artificially perturbed either pRF size or pRF position reduced MEG prediction accuracy, indicating that MEG data are sensitive to pRF properties derived from fMRI. Our model provides a quantitative approach to link fMRI and MEG measurements, thereby enabling advances in our understanding of spatiotemporal dynamics in human visual field maps.

Modulation of the Earliest Component of the Human VEP by Spatial Attention: An Investigation of Task Demands

Spatial attention modulations of initial afferent activity in area V1, indexed by the first component “C1” of the human visual evoked potential, are rarely found. It has thus been suggested that early modulation is induced only by special task conditions, but what these conditions are remains unknown. Recent failed replications—findings of no C1 modulation using a certain task that had previously produced robust modulations—present a strong basis for examining this question. We ran 3 experiments, the first to more exactly replicate the stimulus and behavioral conditions of the original task, and the second and third to manipulate 2 key factors that differed in the failed replication studies: the provision of informative performance feedback, and the degree to which the probed stimulus features matched those facilitating target perception. Although there was an overall significant C1 modulation of 11%, individually, only experiments 1 and 2 showed reliable effects, underlining that the modulations do occur but not consistently. Better feedback induced greater P1, but not C1, modulations. Target-probe feature matching had an inconsistent influence on modulation patterns, with behavioral performance differences and signal-overlap analyses suggesting interference from extrastriate modulations as a potential cause.

Spatial attention modulates earliest visual processing: An electrical neuroimaging study

Several studies showed that shifting of visuospatial attention modulates sensory processing at multiple levels of the visual pathways and beyond, including the occipital striate cortices level. However, inconsistent findings have been reported thus leaving these issues still disputed. 21 participants took part to the present study (the EEG signals of 4 of them were discarded due to artifacts). We used ERPs and their neural sources to investigate whether shifting spatial attention in space across the horizontal meridian of the visual field affected striate cortices activation at the earliest latency. Time-series of scalp topographical maps indicated that, unlike ERPs to attentional-neutral central cues, ERPs to attention-directing local cues showed earliest polarity inversions as a function of stimulated field and processing latency range considered, at occipital-parietal electrodes. In between 60-75 ms, attentional shifting cues elicited a positivity for both visual fields, whereas at a later latency (75–90 ms) they elicited a positivity and a negativity for the upper and lower visual hemifields, respectively. Computed neural sources included striate, besides extrastriate, cortices for both visual hemifields and latency ranges. Conjointly, behavioral responses to targets were faster when they were preceded by local than by neutral cues, and when presented in the upper than the lower hemifield. Our findings support the hypothesis that attention shifts may affect early sensory processing in visual cortices.

Sparse EEG/MEG source estimation via a group lasso

Non-invasive recordings of human brain activity through electroencephalography (EEG) or magnetoencelphalography (MEG) are of value for both basic science and clinical applications in sensory, cognitive, and affective neuroscience. Here we introduce a new approach to estimating the intra-cranial sources of EEG/MEG activity measured from extra-cranial sensors. The approach is based on the group lasso, a sparse-prior inverse that has been adapted to take advantage of functionally-defined regions of interest for the definition of physiologically meaningful groups within a functionally-based common space. Detailed simulations using realistic source-geometries and data from a human Visual Evoked Potential experiment demonstrate that the group-lasso method has improved performance over traditional ℓ2 minimum-norm methods. In addition, we show that pooling source estimates across subjects over functionally defined regions of interest results in improvements in the accuracy of source estimates for both the group-lasso and minimum-norm approaches.

Early retinotopic responses to violations of emotion-location associations may depend on conscious awareness

Reports of modulations of early visual processing suggest that retinotopic visual cortex may actively predict upcoming stimuli. We tested this idea by showing healthy human participants images of human faces at fixation, with different emotional expressions predicting stimuli in either the upper or the lower visual field. On infrequent test trials, emotional faces were followed by combined stimulation of upper and lower visual fields, thus violating previously established associations. Results showed no effects of such violations at the level of the retinotopic C1 of the visual evoked potential over the full sample. However, when separating participants who became aware of these associations from those who did not, we observed significant group differences during extrastriate processing of emotional faces, with inverse solution results indicating stronger activity in unaware subjects throughout the ventral visual stream. Moreover, within-group comparisons showed that the same peripheral stimuli elicited differential activity patterns during the C1 interval, depending on which stimulus elements were predictable. This effect was selectively observed in manipulation-aware subjects. Our results provide preliminary evidence for the notion that early visual processing stages implement predictions of upcoming events. They also point to conscious awareness as a moderator of predictive coding.

From evoked potentials to cortical currents: Resolving V1 and V2 components using retinotopy constrained source estimation without fMRI

Despite evoked potentials' (EP) ubiquity in research and clinical medicine, insights are limited to gross brain dynamics as it remains challenging to map surface potentials to their sources in specific cortical regions. Multiple sources cancellation due to cortical folding and cross-talk obscures close sources, e.g. between visual areas V1 and V2. Recently retinotopic functional magnetic resonance imaging (fMRI) responses were used to constrain source locations to assist separating close sources and to determine cortical current generators. However, an fMRI is largely infeasible for routine EP investigation. We developed a novel method that replaces the fMRI derived retinotopic layout (RL) by an approach where the retinotopy and current estimates are generated from EEG or MEG signals and a standard clinical T1-weighted anatomical MRI. Using the EEG-RL, sources were localized to within 2 mm of the fMRI-RL constrained localized sources. The EEG-RL also produced V1 and V2 current waveforms that closely matched the fMRI-RL's (n = 2) r(1,198)  = 0.99, P < 0.0001. Applying the method to subjects without fMRI (n = 4) demonstrates it generates waveforms that agree closely with the literature. Our advance allows investigators with their current EEG or MEG systems to create a library of brain models tuned to individual subjects' cortical folding in retinotopic maps, and should be applicable to auditory and somatosensory maps. The novel method developed expands EP's ability to study specific brain areas, revitalizing this well-worn technique. Hum Brain Mapp 37:1696-1709, 2016. © 2016 Wiley Periodicals, Inc.

Visual field asymmetries in visual evoked responses

Behavioral responses to visual stimuli exhibit visual field asymmetries, but cortical folding and the close proximity of visual cortical areas make electrophysiological comparisons between different stimulus locations problematic. Retinotopy-constrained source estimation (RCSE) uses distributed dipole models simultaneously constrained by multiple stimulus locations to provide separation between individual visual areas that is not possible with conventional source estimation methods. Magnetoencephalography and RCSE were used to estimate time courses of activity in V1, V2, V3, and V3A. Responses to left and right hemifield stimuli were not significantly different. Peak latencies for peripheral stimuli were significantly shorter than those for perifoveal stimuli in V1, V2, and V3A, likely related to the greater proportion of magnocellular input to V1 in the periphery. Consistent with previous results, sensor magnitudes for lower field stimuli were about twice as large as for upper field, which is only partially explained by the proximity to sensors for lower field cortical sources in V1, V2, and V3. V3A exhibited both latency and amplitude differences for upper and lower field responses. There were no differences for V3, consistent with previous suggestions that dorsal and ventral V3 are two halves of a single visual area, rather than distinct areas V3 and VP.

How to use fMRI functional localizers to improve EEG/MEG source estimation

EEG and MEG have excellent temporal resolution, but the estimation of the neural sources that generate the signals recorded by the sensors is a difficult, ill-posed problem. The high spatial resolution of functional MRI makes it an ideal tool to improve the localization of the EEG/MEG sources using data fusion. However, the combination of the two techniques remains challenging, as the neural generators of the EEG/MEG and BOLD signals might in some cases be very different. Here we describe a data fusion approach that was developed by our team over the last decade in which fMRI is used to provide source constraints that are based on functional areas defined individually for each subject. This mini-review describes the different steps that are necessary to perform source estimation using this approach. It also provides a list of pitfalls that should be avoided when doing fMRI-informed EEG/MEG source imaging. Finally, it describes the advantages of using a ROI-based approach for group-level analysis and for the study of sensory systems.

Towards obtaining spatiotemporally precise responses to continuous sensory stimuli in humans: a general linear modeling approach to EEG

Noninvasive investigation of human sensory processing with high temporal resolution typically involves repeatedly presenting discrete stimuli and extracting an average event-related response from scalp recorded neuroelectric or neuromagnetic signals. While this approach is and has been extremely useful, it suffers from two drawbacks: a lack of naturalness in terms of the stimulus and a lack of precision in terms of the cortical response generators. Here we show that a linear modeling approach that exploits functional specialization in sensory systems can be used to rapidly obtain spatiotemporally precise responses to complex sensory stimuli using electroencephalography (EEG). We demonstrate the method by example through the controlled modulation of the contrast and coherent motion of visual stimuli. Regressing the data against these modulation signals produces spatially focal, highly temporally resolved response measures that are suggestive of specific activation of visual areas V1 and V6, respectively, based on their onset latency, their topographic distribution and the estimated location of their sources. We discuss our approach by comparing it with fMRI/MRI informed source analysis methods and, in doing so, we provide novel information on the timing of coherent motion processing in human V6. Generalizing such an approach has the potential to facilitate the rapid, inexpensive spatiotemporal localization of higher perceptual functions in behaving humans.

Modulation of early cortical processing during divided attention to non-contiguous locations

We often face the challenge of simultaneously attending to multiple non-contiguous regions of space. There is ongoing debate as to how spatial attention is divided under these situations. Whereas, for several years, the predominant view was that humans could divide the attentional spotlight, several recent studies argue in favor of a unitary spotlight that rhythmically samples relevant locations. Here, this issue was addressed by the use of high-density electrophysiology in concert with the multifocal m-sequence technique to examine visual evoked responses to multiple simultaneous streams of stimulation. Concurrently, we assayed the topographic distribution of alpha-band oscillatory mechanisms, a measure of attentional suppression. Participants performed a difficult detection task that required simultaneous attention to two stimuli in contiguous (undivided) or non-contiguous parts of space. In the undivided condition, the classic pattern of attentional modulation was observed, with increased amplitude of the early visual evoked response and increased alpha amplitude ipsilateral to the attended hemifield. For the divided condition, early visual responses to attended stimuli were also enhanced, and the observed multifocal topographic distribution of alpha suppression was in line with the divided attention hypothesis. These results support the existence of divided attentional spotlights, providing evidence that the corresponding modulation occurs during initial sensory processing time-frames in hierarchically early visual regions, and that suppressive mechanisms of visual attention selectively target distracter locations during divided spatial attention.

Optimization of retinotopy constrained source estimation constrained by prior

Studying how the timing and amplitude of visual evoked responses (VERs) vary between visual areas is important for understanding visual processing but is complicated by difficulties in reliably estimating VERs in individual visual areas using noninvasive brain measurements. Retinotopy constrained source estimation (RCSE) addresses this challenge by using multiple, retinotopically mapped stimulus locations to simultaneously constrain estimates of VERs in visual areas V1, V2, and V3, taking advantage of the spatial precision of fMRI retinotopy and the temporal resolution of magnetoencephalography (MEG) or electroencephalography (EEG). Nonlinear optimization of dipole locations, guided by a group-constrained RCSE solution as a prior, improved the robustness of RCSE. This approach facilitated the analysis of differences in timing and amplitude of VERs between V1, V2, and V3, elicited by stimuli with varying luminance contrast in a sample of eight adult humans. The V1 peak response was 37% larger than that of V2 and 74% larger than that of V3, and also ~10-20 ms earlier. Normalized contrast response functions were nearly identical for the three areas. Results without dipole optimization, or with other nonlinear methods not constrained by prior estimates were similar but suffered from greater between-subject variability. The increased reliability of estimates offered by this approach may be particularly valuable when using a smaller number of stimulus locations, enabling a greater variety of stimulus and task manipulations.

The cruciform model of striate generation of the early VEP, re-illustrated, not revoked: a reply to Ales et al. (2013)

Here we summarize the points raised in our dialog with Ales and colleagues on the cortical generators of the early visual evoked potential (VEP), and offer observations on the results of additional simulations that were run in response to our original comment. For small stimuli placed at locations in the upper and lower visual field for which the human VEP has been well characterized, simulated scalp projections of each of the visual areas V1, V2 and V3 invert in polarity. However, the empirically measured, earliest VEP component, "C1," matches the simulated V1 generators in terms of polarity and topography, but not the simulated V2 and V3 generators. We thus conclude that, 1) consistent with the title of Ales et al. (2010a), polarity inversion on its own is not a sufficient criterion for inferring neuroelectric sources in primary visual cortex; but 2) inconsistent with additional claims made in Ales et al. (2010a), the simulated topographies provide additional evidence for - not against - the tenet that the C1 component is generated in V1.

On determining the intracranial sources of visual evoked potentials from scalp topography: a reply to Kelly et al. (this issue)

The cruciform model posits that if a Visual Evoked Potential component originates in cortical area V1, then stimuli placed in the upper versus lower visual field will generate responses with opposite polarity at the scalp. In our original paper (Ales et al., 2010b) we showed that the cruciform model provides an insufficient criterion for identifying V1 sources. This conclusion was reached on the basis of simulations that used realistic 3D models of early visual areas to simulate scalp topographies expected for stimuli of different sizes and shapes placed in different field locations. The simulations indicated that stimuli placed in the upper and lower visual field produce polarity inverting scalp topographies for activation of areas V2 and V3, but not for area V1. As a consequence of the non-uniqueness of the polarity inversion criterion, we suggested that past studies using the cruciform model had not adequately excluded contributions from sources outside V1. In their comment on our paper, Kelly et al. (this issue) raise several concerns with this suggestion. They claim that our initial results did not use the proper stimulus locations to constitute a valid test of the cruciform model. Kelly et al., also contend that the cortical source of the initial visually evoked component (C1) can be identified based on latency and polarity criteria derived from intracranial recordings in non-human primates. In our reply we show that simulations using the suggested critical stimulus locations are consistent with our original findings and thus do not change our conclusions regarding the use of the polarity inversion criterion. We further show that the anatomical assumptions underlying the putatively optimal locations are not consistent with available V1 anatomical data. We then address the non-human primate data, describing how differences in stimuli across studies and species confound an effective utilization of the non-human primate data for interpreting human evoked potential responses. We also show that, considered more broadly, the non-human primate literature shows that multiple visual areas onset simultaneously with V1. We suggest several directions for future research that will further clarify how to make the best use of scalp data for inferring cortical sources.

Generation of the VESPA response to rapid contrast fluctuations is dominated by striate cortex: evidence from retinotopic mapping

The VESPA (visual-evoked spread spectrum analysis) method derives an impulse response function of the visual system from scalp electroencephalographic (EEG) data using the controlled modulation of some feature of a visual stimulus. Recent research using VESPA responses to modulations of stimulus contrast has provided new insights into both early visual attention mechanisms and the specificity of visual-processing deficits in schizophrenia. To allow a fuller interpretation of these and future findings, it is necessary to further characterize the VESPA in terms of its underlying cortical generators. To that end, we here examine spatio-temporal variations in the components of the VESPA as a function of stimulus location. We found that the first two VESPA components (C1/P1) each have a posterior dorsal midline focus and reverse in polarity across the horizontal meridian, consistent with retinotopic projections to calcarine cortex (V1) for the stimulus locations tested. Furthermore, the focal scalp topography of the VESPA was strikingly constant across the entire C1-P1 timeframe (50-120 ms) for each stimulus location, with negligible global scalp activity visible at the zero-crossing dividing the two. This indicates a common focal source underpinning both components, which was further supported by a significant correlation between C1 and P1 amplitudes across subjects (r=0.54; p<0.05). These results, along with factors implicit in the method of derivation of the contrast-VESPA, lead us to conclude that these responses are dominated by activity from striate cortex. We discuss the implications of this finding for previous and future research using the VESPA.

Improved method for retinotopy constrained source estimation of visual-evoked responses

Retinotopy constrained source estimation (RCSE) is a method for noninvasively measuring the time courses of activation in early visual areas using magnetoencephalography (MEG) or electroencephalography (EEG). Unlike conventional equivalent current dipole or distributed source models, the use of multiple, retinotopically mapped stimulus locations to simultaneously constrain the solutions allows for the estimation of independent waveforms for visual areas V1, V2, and V3, despite their close proximity to each other. We describe modifications that improve the reliability and efficiency of this method. First, we find that increasing the number and size of visual stimuli results in source estimates that are less susceptible to noise. Second, to create a more accurate forward solution, we have explicitly modeled the cortical point spread of individual visual stimuli. Dipoles are represented as extended patches on the cortical surface, which take into account the estimated receptive field size at each location in V1, V2, and V3 as well as the contributions from contralateral, ipsilateral, dorsal, and ventral portions of the visual areas. Third, we implemented a map fitting procedure to deform a template to match individual subject retinotopic maps derived from functional magnetic resonance imaging (fMRI). This improves the efficiency of the overall method by allowing automated dipole selection, and it makes the results less sensitive to physiological noise in fMRI retinotopy data. Finally, the iteratively reweighted least squares (IRLS) method was used to reduce the contribution from stimulus locations with high residual error for robust estimation of visual evoked responses.

Increasing the accuracy of electromagnetic inverses using functional area source correlation constraints

Estimating cortical current distributions from electroencephalographic (EEG) or magnetoencephalographic data is a difficult inverse problem whose solution can be improved by the addition of priors on the associated neural responses. In the context of visual activation studies, we propose a new approach that uses a functional area constrained estimator (FACE) to increase the accuracy of the reconstructions. It derives the source correlation matrix from a segmentation of the cortex into areas defined by retinotopic maps of the visual field or by functional localizers obtained independently by fMRI. These areas are computed once for each individual subject and the associated estimators can therefore be reused for any new study on the same participant. The resulting FACE reconstructions emphasize the activity of sources within these areas or enforce their intercorrelations. We used realistic Monte-Carlo simulations to demonstrate that this approach improved our estimates of a diverse set of source configurations. Reconstructions obtained from a real EEG dataset demonstrate that our priors improve the localization of the cortical areas involved in horizontal disparity processing.

Retinotopic mapping of the primary visual cortex - a challenge for MEG imaging of the human cortex

Magnetoencephalography (MEG) can be used to reconstruct neuronal activity with high spatial and temporal resolution. However, this reconstruction problem is ill-posed, and requires the use of prior constraints in order to produce a unique solution. At present there are a multitude of inversion algorithms, each employing different assumptions, but one major problem when comparing the accuracy of these different approaches is that often the true underlying electrical state of the brain is unknown. In this study, we explore one paradigm, retinotopic mapping in the primary visual cortex (V1), for which the ground truth is known to a reasonable degree of accuracy, enabling the comparison of MEG source reconstructions with the true electrical state of the brain. Specifically, we attempted to localize, using a beanforming method, the induced responses in the visual cortex generated by a high contrast, retinotopically varying stimulus. Although well described in primate studies, it has been an open question whether the induced gamma power in humans due to high contrast gratings derives from V1 rather than the prestriate cortex (V2). We show that the beanformer source estimate in the gamma and theta bands does vary in a manner consistent with the known retinotopy of V1. However, these peak locations, although retinotopically organized, did not accurately localize to the cortical surface. We considered possible causes for this discrepancy and suggest that improved MEG/magnetic resonance imaging co-registration and the use of more accurate source models that take into account the spatial extent and shape of the active cortex may, in future, improve the accuracy of the source reconstructions.

Effects of attentional load on early visual processing depend on stimulus timing

A growing number of studies suggest that early visual processing is not only affected by low-level perceptual attributes but also by higher order cognitive factors such as attention or emotion. Using high-density electroencephalography, we recently demonstrated that attentional load of a task at fixation reduces the response of primary visual cortex to irrelevant peripheral stimuli, as indexed by the C1 component. In the latter study, peripheral stimuli were always presented during intervals without task-relevant stimuli. Here, we use a similar paradigm but present central task stimuli and irrelevant peripheral stimuli simultaneously while keeping all other stimulus characteristics constant. Results show that rather than to suppress responses to peripheral stimulation, high attentional load elicits higher C1 amplitudes under these conditions. These findings suggest that stimulus timing can profoundly alter the effects of attentional load on the earliest stages of processing in human visual cortex.

Top-down effects on early visual processing in humans: a predictive coding framework

An increasing number of human electroencephalography (EEG) studies examining the earliest component of the visual evoked potential, the so-called C1, have cast doubts on the previously prevalent notion that this component is impermeable to top-down effects. This article reviews the original studies that (i) described the C1, (ii) linked it to primary visual cortex (V1) activity, and (iii) suggested that its electrophysiological characteristics are exclusively determined by low-level stimulus attributes, particularly the spatial position of the stimulus within the visual field. We then describe conflicting evidence from animal studies and human neuroimaging experiments and provide an overview of recent EEG and magnetoencephalography (MEG) work showing that initial V1 activity in humans may be strongly modulated by higher-level cognitive factors. Finally, we formulate a theoretical framework for understanding top-down effects on early visual processing in terms of predictive coding.

Phase delays within visual cortex shape the response to steady-state visual stimulation

Although the spatial organization of visual areas can be revealed by functional Magnetic Resonance Imaging (fMRI), the synoptic, non-invasive access to the temporal characteristics of the information flow amongst distributed visual processes remains a technical and methodological challenge. Using frequency-encoded steady-state visual stimulation together with a combination of time-resolved functional magnetic source imaging from magnetoencephalography (MEG) and anatomical magnetic resonance imaging (MRI), this study evidences maps of visuotopic sustained oscillatory neural responses distributed across the visual cortex. Our results further reveal relative phase delays across responding striate and extra-striate visual areas, which thereby shape the chronometry of neural processes amongst these regions. The methodology developed in this study points at further developments in time-resolved analyses of distributed visual processes in the millisecond range, and to new ways of exploring the dynamics of functional processes within the human visual cortex non-invasively.

V1 is not uniquely identified by polarity reversals of responses to upper and lower visual field stimuli

The cruciform hypothesis states that if a visual evoked potential component originates in V1, then stimuli placed in the upper versus lower visual fields will generate responses with opposite polarity at the scalp. This diagnostic has been used by many studies as a definitive marker of V1 sources. To provide an empirical test of the validity of the cruciform hypothesis, we generated forward models of cortical areas V1, V2 and V3 that were based on realistic estimates of the 3-D shape of these areas and the shape and conductivity of the brain, skull and scalp. Functional MRI was used to identify the location of early visual areas and anatomical MRI data was used to construct detailed cortical surface reconstructions and to generate boundary element method forward models of the electrical conductivity of each participant's head. These two data sets for each subject were used to generate simulated scalp activity from the dorsal and ventral subdivisions of each visual area that correspond to the lower and upper visual field representations, respectively. The predicted topographies show that sources in V1 do not fully conform to the cruciform sign-reversal. Moreover, contrary to the model, retinotopic visual areas V2 and V3 show polarity reversals for upper and lower field stimuli. The presence of a response polarity inversion for upper versus lower field stimuli is therefore an insufficient criterion for identifying responses as having originated in V1.