Abstract
Shadows cast by retinal blood vessels are represented in striate cortex of the squirrel monkey. Their pattern was exploited to generate a true retinotopic map of V1. For calibration, retinal landmarks were projected onto a tangent screen to measure their visual field location. Next, the retina was warped onto striate cortex, distorting it as necessary to match each retinal vessel to its cortical representation. Maps from four hemispheres of two normal adult squirrel monkeys were created and used to derive expressions for cortical magnification factor (M). A mean map was produced by averaging the individual maps. To address the controversial issue of whether the ratio of retinal ganglion cell (RGC) density to M is constant at all eccentricities, we stained a retinal whole mount from one of the two monkeys for Nissl substance. A ganglion cell density map was compiled by sampling the concentration of cells at 171 retinal points. Allowance was made for displaced amacrine cells and for the centripetal displacement of RGCs from central photoreceptors. After these corrections the V1 surface area and RGC density were compared at each eccentricity. The cortical representation of the macula was found to be amplified, even beyond the magnification expected from its high density of RGCs. For example, the central 4° of visual field were allotted 27% of the surface area of V1 but were supplied by only 12% of RGCs. We conclude that, in monkey striate cortex, more tissue is allocated per ganglion cell for the analysis of information emanating from the macula as compared with the peripheral retina.
- blind spot
- monocular crescent
- retina
- blood vessel
- flat-mount
- ocular dominance column
- cytochrome oxidase
- magnification factor
- anisotropy
- displaced amacrine cell
- GABA
- retinal ganglion cell
- cone
- macula
- Henle fiber layer
Introduction
Where explorers go, they depict new territory with maps. The first map of any area in the brain was prepared by an ophthalmologist, Tatsuji Inouye (1909). He correlated visual field scotomas with the location of cranial gunshot wounds in survivors of the Russo-Japanese war to compile a map of the visual field representation in the primary visual cortex. Similar studies were conducted in veterans of later wars (Holmes and Lister, 1916; Spalding, 1952; Teuber et al., 1960). More recently, maps have been made by imaging focal occipital lesions in patients with visual field defects (McAuley and Russell, 1979; Kattah et al., 1981; Spector et al., 1981;Horton and Hoyt, 1991). Maps also can be prepared by using positron emission tomography (Fox et al., 1987) or functional magnetic resonance imaging (Sereno et al., 1995; DeYoe et al., 1996; Engel et al., 1997) in normal subjects. The main problem with the maps resulting from these approaches is their limited spatial resolution.
More precise maps can be prepared in animals by recording from cells with microelectrodes. Talbot and Marshall (1941) mapped the occipital surface (operculum) in the macaque and quantified “cortical discrimination” in terms of degrees of visual angle per millimeter of cortex. Daniel and Whitteridge (1961) plotted the reciprocal, “cortical magnification factor,” after making electrode penetrations through different regions of striate cortex. They found that cortical magnification fell by several log units from the foveal representation to the monocular crescent representation, reflecting the greater allocation of cortex to central vision. From their measurement of cortical magnification, they predicted the shape of V1 and the layout of its retinotopic map. Subsequent electrophysiological studies have elaborated on their conclusions (Cowey, 1964; Allman and Kaas, 1971; Dow et al., 1981; Van Essen et al., 1984; Rosa et al., 1997).
The resolution of microelectrode recordings is limited by eye movements, receptive field scatter, inaccuracy in plotting receptive fields, and errors in measuring coordinates in the visual field. It is a daunting task to reconstruct multiple electrode penetrations, especially through folded cortex, so that each recording site can be matched accurately with its corresponding visual field locus. As pointed out by Van Essen et al. (1984), the deoxyglucose method is inherently more precise than microelectrode recordings for the construction of a retinotopic map. After [14C]deoxyglucose administration the retina can be stimulated with a ring and ray stimulus to produce a metabolic imprint in the cortex (Tootell et al., 1988). However, deoxyglucose maps cannot be prepared easily from buried, folded cortex. In monkeys much of striate cortex lies hidden in the calcarine fissure, beyond the reach of the deoxyglucose method. This tissue is also inaccessible to optical imaging (Blasdel and Campbell, 2001).
Recently, we have shown that retinal blood vessels are represented in striate cortex of the squirrel monkey (Adams and Horton, 2002). They can be detected by staining flat-mounts of the cortex for cytochrome oxidase (CO) after enucleation of one eye. In this paper we used this phenomenon to generate a retinotopic map by matching each segment of the retinal vascular tree to its cortical image. The resulting map is extremely accurate, because it is based on alignment between a multitude of cortical and retinal landmarks. It is immune to many sources of error because it is literally a retinotopic, rather than a visuotopic, map.
With the aid of this retinotopic map, we have addressed a number of questions long debated by cartographers. (1) Is the mathematical expression for cortical magnification factor a complex log function? (2) Is cortical magnification factor equal along the vertical and horizontal meridia? (3) Is local cortical magnification factor equal along isopolar and isoeccentricity contours at any locus in the visual field representation? (4) Is the ratio of retinal ganglion cell (RGC) density to cortical magnification factor a constant at all eccentricities?
Materials and Methods
This paper describes a mean retinotopic map generated from four hemispheres of two normal squirrel monkeys (Saimiri sciureus) named Monkey Q and Monkey P. A description of our methods for visualizing the cortical representation of angioscotomas was presented previously (Adams and Horton, 2002). The steps required to create a retinotopic map of striate cortex from angioscotoma representations are described in Results. All experimental procedures were approved by the University of California, San Francisco, Committee on Animal Research.
Retinal ganglion cell density. The left eye of Monkey P was enucleated under general anesthesia with 2% isoflurane. After fixation by immersion in 2% paraformaldehyde, a whole mount was prepared by dissecting the retina from the underlying pigment epithelium (Stone, 1981). Six radial relieving cuts permitted the retina to be mounted on a 3 × 2“ subbed slide, photoreceptor side down. It was stained for Nissl substance via an aqueous chrome alum-gallocyanine method (Kiernan, 1990), coverslipped under glycerol, and sealed with nail varnish. A ring and ray pattern was superimposed onto the whole mount, using 14 calibration points determined in vivo by projection of retinal landmarks (e.g., vessel crossings, bifurcations) onto a tangent screen. Counting fields were selected in the middle of each compartment (see Fig. 13D,E) and located at high power, using retinal blood vessels as guides.
A 40× Plan-Apochromat oil immersion lens (Zeiss, Thornwood, NY) and a Spot RT digital camera (Diagnostic Instruments, Sterling Heights, MI) were used to display an image of the RGC layer onto a computer monitor at a net magnification of 2000×. Parafoveally, where the RGC layer was thickest, it was necessary to focus the microscope 60 μm into the retina to visualize the deepest cells. The outline of each RGC nucleus located within a circular field was traced manually onto a sheet of transparent film taped to the monitor screen. Nuclei were drawn only if their center fell within the boundaries of the window printed onto the transparency film. After all cells were identified in a given plane, we focused more superficially, adding nuclear profiles as they came into view. After reaching the inner limiting membrane, we focused back through the tissue to make sure that every nuclear profile was outlined. Then the total number of outlined profiles on the transparency sheet was counted. Sampling windows were 2000 μm2at 0.5, 1.5, 3, 4.5, and 6°; 8000 μm2at 12°; and 20,000 μm2 at 20, 28, 41, and 61°. Larger sampling windows were used at higher eccentricity to compensate for the lower density of ganglion cells. The raw counts were converted to cells/mm2.
Counts were converted from cells/mm2 to cells/degrees squared. This was done by comparing the distance (in μm) between the fovea and landmarks in the whole mount to their distance (in degrees) measured in vivo by projection of landmarks onto a tangent screen. This approach eliminated the issue of postmortem retinal shrinkage. Local retinal magnification factor averaged 161 μm/degree within the central 24°. Beyond this eccentricity, accurate retinal photography and projection of landmarks were impossible. We therefore scaled a wide-angle schematic eye (Drasdo and Fowler, 1974) to our whole mount to derive values for retinal magnification factor of 159 μm/degree at 28°, 149 μm/degree at 41°, and 116 μm/degree at 61°. It is worth mentioning that shrinkage of striate cortex was not measured but was probably negligible because animals were perfused with only 1% paraformaldehyde.
There is a danger of underestimating RGC concentration centrally, because cells are packed so densely in the macula of whole mounts, especially after dehydration. To avoid this pitfall, Wässle et al. (1990) cut out the central retina from whole mounts and prepared serial cross sections. This approach, however, requires a mathematical approximation to account for nuclear fragments (Abercrombie, 1946).Curcio et al. (1987) achieved accurate cell counts in whole mounts by examining unstained, hydrated tissue with Nomarski differential interference contrast (DIC) microscopy. In hydrated tissue the ganglion cell layer remains 50–60 μm thick in the macula, rather than collapsing to 20 μm (Curcio and Allen, 1990). If one avoids shrinkage of the retina, it is possible to obtain accurate ganglion cell counts in the macula of whole mounts. We used an aqueous-based Nissl stain, which allowed us to maintain our retinal whole mount in a hydrated state and to use conventional light-field optics to distinguish cell types.
Displaced amacrine cell density. So that displaced amacrine cells from RGCs could be differentiated, cross sections of squirrel monkey retina from another eye were stained for GABA immunoreactivity (Hendrickson et al., 1985; Hendry et al., 1987; Wässle et al., 1987, 1990). The tissue was postfixed in 2% glutaraldehyde for 24 hr, embedded in plastic, and sectioned at 3 μm. Sections were air-dried on slides, and the plastic was dissolved in NaOH (Lane and Europa, 1965). Then the retina was incubated with ICN Biomedical Clone 5A9 monoclonal anti-GABA antibody (Irvine, CA) for 2 d at a dilution of 1:30. A Vectastain Elite ABC kit (Vector Laboratories, Burlingame, CA) was used to visualize anti-GABA immunoreactivity with a diaminobenzidine peroxidase reaction product.
Monocular crescent. Peripherally, we determined the eccentricity of the monocular crescent border by placing an anesthetized squirrel monkey in a stereotaxic apparatus and using a light source mounted on an arc perimeter. As the light entered the extreme nasal field, we noted when the corneal light reflex appeared at the nasal limbus. This occurred at an eccentricity of 72° from the midline. The animal's interpupillary distance was 18 mm, and the corneal apex was situated 4 mm posterior to the frontal plane defined by the bridge of the nose. In the human and macaque the monocular crescent defined by the ora serrata of the peripheral temporal retina is smaller than the monocular crescent defined by the nasal bridge. Because the squirrel monkey's eyes are set more closely and the nasal bridge is lower, we assumed that the eccentricity of the monocular crescent border corresponded closely to that of the temporal retinal perimeter. The validity of this assumption was tested by determining the fraction of striate cortex occupied by the monocular crescent representation. In squirrel monkey it measured only 3.2% (n = 4), whereas in macaque (Horton and Hocking, 1996b) and human (our unpublished data) it measures 5%. It was significantly smaller in squirrel monkey (p < 0.0001, Student's two-tailed t test), reflecting the fact that it begins at ∼72°, whereas in macaque and human it begins at ∼60°.
Results
A retinotopic map from angioscotomas
Figure 1 illustrates the process of mapping the representation of retinal blood vessels in striate cortex for Monkey P. Both retinas were photographed and calibrated by projection of landmarks onto a tangent screen. After removal of the left eye, striate cortex was flat-mounted and stained for CO to visualize the pattern of ocular dominance columns. Emanating from each optic disc representation were 7–10 serpentine lines, representing major blood vessels in the retina. As explained previously (Adams and Horton, 2002), deprivation caused by shadows from retinal blood vessels induces the rearrangement of geniculocortical afferents after birth. The pattern formed by these remodeled geniculocortical afferents can be detected in adults by monocular enucleation, followed by CO staining.
Figure 2 highlights the steps involved in generating a retinotopic map by using cortical angioscotoma representations. A ring and ray pattern representing the visual hemifield, calibrated in the living animal, was superimposed on the photograph of the nasal retina. Correspondence points, marked with yellow squares, were placed at points at which the location could be determined unambiguously in the retina and the cortex. Ideal correspondence points were provided by vessel bifurcations, crossings, or inflection points. Features such as the borders of V1, the optic disc, and the monocular crescent provided additional retinotopic data.
After correspondence points were identified, they were warped onto the corresponding cortical elements by using a computer program, Elastic Reality 3.0 (Avid Technology, Tewksbury, MA), which performs a linear interpolation between points. The resulting distortion of the ring and ray pattern produced a retinotopic map of each cortex. Each map was based on ∼60 correspondence points, plus intercalated segments.
No angioscotomas are represented at eccentricities <4°, because retinal blood vessels become too narrow to cast a dense shadow (Adams and Horton, 2002). Cowey (1964) mapped the central 4°, exposed on the smooth convexity of the squirrel monkey occipital cortex, by plotting the receptive field location of 48 electrode recording sites. His detailed retinotopic map of the central 4° was inserted into our map by scaling to the V1/V2 border, fovea, and location of the medial fold of the operculum.
A retinotopic map was prepared from both hemispheres of two monkeys, Monkey Q and Monkey P (Fig. 3). For each map we digitized 172 points, corresponding to each ring–ray intersection. The average map was prepared by reflecting the right cortices onto the left cortices. The individual maps were not scaled or distorted in any way. All four maps were aligned at the fovea (explaining the absence of error bars at the foveal point). They then were rotated to align the point at which the horizontal meridian intersects the monocular crescent (explaining the absence of a vertical error bar at the point marked by the asterisk). We plotted the meanx and y coordinate of each point and connected them with spline curves.
We chose two animals with similar cortical shapes and areas to minimize the size of the error bars. There exists substantial variation from animal to animal in the surface area of striate cortex (Van Essen et al., 1984) and in the shape of striate cortex. Both factors cause variation in the retinotopic map. To construct a maximally precise map, we deliberately chose animals with similar cortices.
The mean map had a surface area of 798 mm2, close to the average area of squirrel monkey V1 reported previously (Horton and Hocking, 1996a). The monocular crescent occupied 23.5 mm2, or 2.9% of the total striate area. This was less than in the macaque, in which the monocular crescent averages 5.0% of striate cortex (Horton and Hocking, 1996a). In binocular cortex the upper field quadrant measured 417 mm2 (54%), and the lower field quadrant measured 357 mm2 (46%). We found no evidence for a greater representation of the lower field, as suggested for macaques (Van Essen et al., 1984; Tootell et al., 1988). The total perimeter of striate cortex was 104.3 mm. The upper vertical meridian measured 48.5 mm, the lower vertical meridian measured 44.5 mm, and the horizontal meridian measured 30.5 mm (each measured from the fovea to the border of the monocular crescent representation). The vertical meridia were 1.52 times longer than the horizontal meridia. This exceeds the ratio of 1.25 found in the macaque (Daniel and Whitteridge, 1961; Tootell et al., 1988), reflecting the more circular shape of striate cortex in the squirrel monkey.
Cortical magnification factor
The most striking feature of the map was the exaggerated representation of the central visual field. This was quantified by measuring M, the cortical magnification factor (Daniel and Whitteridge, 1961). At a particular eccentricity M can be measured in two different ways. (1) Areal magnification (Ma) is the cortical surface area devoted to 1 square degree of visual field. Its units are mm2/degrees squared. (2) Linear magnification is the cortical distance corresponding to 1° in the visual field. Its units are mm/degree. Linear magnification commonly is measured either along a ray (Mp) or a ring (Me).
To calculate areal magnification, we measured (in mm2) the cortical surface area of each compartment in each of the four retinotopic maps. Next we calculated (in degrees squared) the area (i.e., solid angle) of each corresponding compartment in the visual field. Dividing the two values gave the areal magnification in mm2/degrees squared.
To calculate linear magnification along polar rays (Mp), we measured the distance along each ray between two eccentricity rings. This distance was divided by the angle subtended by the corresponding ray segment in the visual field (in degrees). This angle was easy to calculate because each ray is analogous to a longitude line on a globe. For example, the visual angle between the 32° and 16° ring was equal to 16°.
To calculate linear magnification along eccentricity rings (Me), we measured the distance along each ring between two polar rays (analogous to a latitude line on a globe). This distance was divided by the visual angle subtended by the corresponding ring segment in the visual field (in degrees).
To show how magnification varies over the retinotopic map, we plotMa,Mp, andMe against eccentricity in Figure4. Mawas averaged over all 18 compartments in one-half of a spherical segment and plotted with SE bars at an eccentricity midway between two rings (Fig. 4A). BecauseMa declines in an approximately logarithmic way, this slightly overestimatesMa but has the advantage of making no assumptions about its mathematical function (see Tootell et al., 1988).Mp was averaged over all 19 ray segments between two neighboring rings and plotted midway between them (Fig. 4B). Me was averaged over all 18 segments along a given ring and plotted exactly at its eccentricity (Fig. 4C).
The first map of primary visual cortex, compiled by Inouye (1909), is shown in Figure 5A. Polar angle is represented by horizontal lines, and eccentricity is represented by vertical lines. Although Inouye may not have understood its mathematical implications fully (see translation by Glickstein and Fahle, 2000), the map is recognizable as the complex logarithmic transformation of the visual field (Schwartz, 1977, 1980). It is physically impossible, because the fovea is represented as a line. Consequently, cortical magnification approaches infinity as eccentricity drops to 0°. Electrophysiological recordings by Talbot and Marshall (1941) later showed convincingly that the isopolar rays converge to a point at the representation of the fovea. All subsequent maps have agreed on this principle.
Using retinotopic data from different sources, various authors have attempted to model mathematically how M varies as a function of eccentricity (Daniel and Whitteridge, 1961; Fischer, 1973; Hubel and Freeman, 1977; Schwartz, 1977; Van Essen et al., 1984; LeVay et al., 1985; Schein and de Monasterio, 1987; Johnston, 1989). The general formula for magnification factor (Van Essen et al., 1984) is: where E is eccentricity and a, b, and c are constants. The constant a is a scaling factor. The constant b was introduced by Schwartz (1977). It prevents M from climbing to infinity as Eapproaches 0°, thus correcting the error in Inouye's map.
If the visual field map in striate cortex conformed to the complex logarithmic mapping function of Schwartz (1984), the constantc would equal −1 for linear magnification and −2 for areal magnification. We fit curves to our data to derive the constants in the magnification factor expression (Fig. 4). Either the value ofc was restricted to obey the complex logarithmic function, or it was allowed to vary according to the value that yielded the best fit (Van Essen et al., 1984).
The value of c for Mp was close to −1, indicating that Mpdecreased logarithmically with increasing eccentricity. However, the value of c for Me was less than would be expected for a purely conformal map. In a conformal map the upper and lower vertical meridia should be parallel in the periphery (Fig. 5B). In the squirrel monkey (Fig. 3) the upper and lower vertical meridia draw back together in the periphery after reaching a maximum separation at 8°. This shortens the representation of the isoeccentricity rings, reducing c to less than −1 for Me and less than −2 for Ma. Macaque V1 has a more oval shape than squirrel monkey V1. Therefore, the vertical meridia are more parallel in the midperiphery. Schwartz proposed that the complex log model fits macaque V1 out to 30°. Beyond that eccentricity the vertical meridia converge, so the map deviates from the shape of the theoretical map shown in Figure 5B (Schwartz, 1977, 1980,1983).
In fact, even in macaque the vertical meridia are never truly parallel. As in squirrel monkey, the cortex reaches a maximum width at 8° and then begins to taper (Horton and Hocking, 1996b, 1998). Therefore, even the macaque does not conform to the complex log model. Although one can modify further the equation for cortical magnification factor in an attempt to approximate the actual shape of striate cortex (seeSchwartz, 1983), it is doubtful that the visual field representation is governed by any strict mathematical principle in any species.
The cortical map is described as isotropic if values forMe andMp are equal at any given locus. The curves in Figure 4, B and C, depicting mean values for Me andMp are not identical. To compare their values, we calculated the anisotropy index,Mp/Me(Van Essen et al., 1984), for each compartment in the cortical map (Fig. 6). To make this calculation, one must derive values for Mp andMe at the center of each compartment. For Mp we averaged the value ofMp for two adjacent ray segments. ForMe we dividedMa by meanMp for each compartment.
Figure 6 shows that the anisotropy index varies systematically across the visual field map. As in the macaque (Van Essen et al., 1984), it tends to be greater along the vertical meridian (mean value, 1.56) than along the horizontal meridian (mean value, 0.88). This occurs because the path length of each polar ray increases systematically the closer its location to the vertical meridian. The anisotropy index also becomes larger at more peripheral eccentricities (e.g., >16°). This occurs because the convergence of polar rays near the monocular crescent representation shortens the isoeccentricity lines, reducing the relative value of Me.
To illustrate how the anisotropy index affects the representation of the visual field, we projected a field of identical squares onto striate cortex (Fig. 6), using the same warping technique used to create our mean map. Squares near the foveal representation are magnified enormously, as predicted from theMa ratio of ∼10,000:1 from fovea/periphery (Fig. 4A). Along the horizontal meridian the squares appear elongated in the cortex, with their long axes generally oriented perpendicular to the horizontal meridian. This distortion (approximately speaking, from squares to rectangular trapezoids) is produced by the decrement in global cortical magnification from center to periphery. Along the vertical meridia the squares also are elongated, but with their long axes parallel to the vertical meridia. Stretching occurs in this direction despite declining global cortical magnification from center to periphery because of the high anisotropy index along the vertical meridian. In squirrel monkeyBlasdel and Campbell (2001) have used optical imaging with square wave gratings to show this anisotropy along the vertical meridian.
As mentioned above, higher anisotropy indices result from the greater length of isopolar rays near the vertical meridia. To examine this effect more closely, we graphed Mpversus polar angle between each adjacent pair of isoeccentricity rings (Fig. 7A). Each curve is “hammock” shaped; that is, Mpdeclines toward the horizontal meridian at each eccentricity.Me increases slightly toward the horizontal meridian (data not shown), thereby maintainingMa nearly constant at each polar angle (Fig. 7B) for a given eccentricity. In simplest terms the map compartments (Fig. 3) are more rectangular near the vertical meridian than near the horizontal meridian, but their net area is quite constant. We did not observe an increase inMa near the horizontal meridian in our mean map, as reported by Van Essen et al. (1984) for the macaque. However, there appears to be real variability in this property among individual squirrel monkeys, judging from the variation we have observed in the shape of the blind spot representation (oval versus circular) in normal animals.
Visual field domain of a cortical patch
In many species, especially primates, striate cortex contains a regular array of cytochrome oxidase patches (Horton and Hubel, 1981;Carroll and Wong-Riley, 1984; Horton, 1984). Farias et al. (1997) have emphasized that in macaques their spatial density is quite constant throughout binocular cortex. The same is true in the squirrel monkey. Figure 8 shows a section through layer 3 from the right striate cortex of Monkey Q, an animal used to prepare our retinotopic map. There were 2534 patches in 792 mm2 of binocular cortex. The patches seemed to be evenly distributed. To confirm this impression, we measured their density at different eccentricities. It ranged between 2.77 and 3.95 patches/mm2. By calculating the area of each spherical zone in the visual field, we determined the number of patches/degree squared. These data are graphed in Figure9. Because the density of patches in the cortex is quite constant, the graph bears a close resemblance to the graph of Ma versus eccentricity (Fig.4A).
The reciprocal of patches/degree squared yields a function called the “blob (patch) domain” (Tootell et al., 1988). This corresponds to the amount of visual field served by a CO patch at a given eccentricity. In the squirrel monkey it ranges from 0.008 degrees squared in the central 1° to 29.3 degrees squared in the periphery (Fig. 9). The most vivid impression of how patch domain varies with eccentricity can be obtained by back-projecting the CO patches (Fig. 8) onto the visual field (Fig. 10). Analogous back-projections have been performed for the ocular dominance columns (Hubel and Freeman, 1977; LeVay et al., 1985). The back-projection was accomplished by warping the cortical image of the patches onto the polar coordinates of the visual hemifield, essentially reversing the process used to create cortical maps (Fig. 2). In the peripheral field the patches are sparse because they each cover a large territory, whereas near the central field they are packed so densely that the image must be magnified to see them. Their distribution mirrors the gradient in Ma. Because the patches vary in their shape and pattern, it is difficult to infer any underlying principle from their local appearance in the visual field projection.
Retinal ganglion cell density
Ganglion cells are the output elements of the retina. Peripherally, the ganglion cell layer contains only a scattering of cells, but centrally it forms a tightly packed sheet, several cells deep. It is controversial if the ratio of RGC density to V1 surface area is constant at all eccentricities (for review, see Pointer, 1986;Tolhurst and Ling, 1988). Previous analyses have been hampered by difficulties in measuring ganglion cell density, compensating for Henle fiber layer displacement, translating retinal points onto the visual field, and mapping accurately striate cortex. To make matters worse, ganglion cell densities derived from one animal have been compared with cortical maps prepared in another, introducing additional variation.
We decided to reexamine the relationship between ganglion cell density and Ma, because our experimental approach has minimized three obstacles mentioned above. First, we made ganglion cell counts in one of the monkeys used to generate our retinotopic map. Second, retinal landmarks were projected onto a tangent screen, allowing us to determine experimentally the visual field coordinates of retinal compartments (some extrapolation was required beyond 24°). Third, the accuracy of our cortical map allowed for valid measurements of Ma.
Figure 11A shows a Nissl-stained whole mount of the left retina from one of the two monkeys used to generate the mean map. In the temporal retina the border is visible between the neurosensory retina and the pars plana. In the living squirrel monkey it corresponds to an eccentricity of ∼72°. A circle of this radius defines the extent of the visual hemifield seen with both eyes (the binocular visual field). We therefore placed the 72° isoeccentricity ring at this border. Nasal retina located outside this ring serves the monocular crescent of the visual field. It was not included in our analysis.
In Figure 11B a montage of fundus pictures (Fig.1A) was superimposed on the whole mount by matching the location of the fovea, optic disc, and major blood vessels. There was a good correspondence between blood vessel landmarks in the whole mount and the photograph. Because these landmarks were calibrated by projection onto a tangent screen, we could use the retinal photograph to determine the location of rings within the central 24°. In this zone retinal magnification factor is nearly constant (Lapuerta and Schein, 1995). In this squirrel monkey we obtained a mean value of 161 μm/degree. Beyond 24° retinal magnification declines because of the optical characteristics of the eye (Hughes, 1976; Holden and Fitzke, 1988; Bennett et al., 1994). With our fundus camera we could not photograph the peripheral retina or project its blood vessels onto a tangent screen. This made it impossible to locate the 32° and 50° isoeccentricity rings experimentally. Experimental measurement of peripheral retinal magnification factor is an unsolved technical challenge. In the primate there has been only one attempt, making laser lesions at peripheral eccentricities and then examining the retina histologically after enucleation (Frisén and Schöldström, 1977). As an alternative, investigators have relied on ray tracing in schematic eyes to project the peripheral visual field onto the retina. We used a widely accepted model of the optical function of the eye, the wide-angle schematic eye of Drasdo and Fowler (1974), to position the 32° and 50° isoeccentricity rings. This was done by scaling Drasdo and Fowler's projection of the visual field [see Drasdo and Fowler (1974), their Fig. 6] to the retina illustrated in Figure 11B.
The ring–ray pattern was superimposed onto the retina, allowing us to establish the relationship between coordinates in the visual field and physical location on the whole mount (Fig. 11C). In the far periphery, where relieving cuts were located, it was necessary to bend the rays. To accomplish this, we split the ray pattern down the middle of each relieving cut. Then the seams were warped, contracting the peripheral rays onto the retina (Fig. 11D).
Figure 12 shows representative fields of cells at 0.5, 6, 20, and 61°. Ganglion cells are recognizable by their large round nuclei, prominent nucleoli, and relative size (DeBruyn et al., 1980; Leventhal et al., 1981; Stone and Johnston, 1981; Perry et al., 1984). They can be differentiated easily from spindle-shaped endothelial cells lining blood vessels. Astrocytes and microglia are rarely present in the ganglion cell layer (Curcio and Allen, 1990). At any given eccentricity they are smaller than ganglion cells and have an oval, dark nucleus without a prominent nucleolus. Ganglion cell density was sampled at 162 locations, between each pair of isoeccentricity rings at every 20° of polar angle (see yellow dots in Fig. 11D,E). The size of the counting window was adjusted with eccentricity so that greater areas of retina were sampled peripherally where cell density is less. In total, 13,616 cells were counted throughout the retina. The peak ganglion cell density was between 3 and 6°. To sample this interval more closely, we made nine additional cell density measurements at 4.5°.
Displaced amacrine cells
Displaced amacrine cells are indistinguishable from ganglion cells in Nissl and therefore contaminated our counts. It was necessary to adjust the raw ganglion cell counts to correct for them. This was done by GABA immunostaining the displaced amacrine cells in a different retina and subtracting their proportion from the ganglion cell counts obtained at each eccentricity. A belt of retina from the fovea to the nasal pars plana was cut from a whole mount. The 2 × 10 mm piece of retina was folded in half, embedded in plastic, sectioned at 3 μm, and processed for GABA.
Figure 13 shows a single immunostained section passing through the fovea and the optic disc. The inner plexiform layer was heavily labeled. Contrary to a previous study in the squirrel monkey showing homogeneous labeling (Brecha, 1983), there were three distinct bands of heavy immunoreactivity, as described in the macaque (Hendrickson et al., 1985; Nishimura et al., 1985). Darkly labeled cells were present in both the inner nuclear layer and the ganglion cell layer (Koontz et al., 1989; Grünert and Wässle, 1990; Koontz and Hendrickson, 1990; Wässle et al., 1990; Crooks and Kolb, 1992). GABA immunohistochemistry is useful for identifying displaced amacrine cells, because 90% stain positively (Wässle et al., 1990; Koontz et al., 1993), whereas ganglion cells are nonreactive (Grünert and Wässle, 1990;Wässle et al., 1990). We counted displaced amacrine cells and ganglion cells at nine eccentricities, centered between isoeccentricity rings, using a 40× oil immersion objective with DIC. The mean proportion of GABA-positive cells at each eccentricity was determined by examining 11 sections, spaced 150–200 μm apart. For the 12° measurement only six sections were examined because the optic disc was present in five sections. We measured displaced amacrine cell proportion, rather than absolute density, eliminating the need for correction factors that are required when cells are counted in cross sections (Abercrombie, 1946).
Four examples of counting windows are shown in Figure13A–D. The amacrine cells constituted a variable percentage of the ganglion cell layer at each eccentricity: 0.5° (0%), 1.5° (4.7%), 3° (3.5%), 4.5° (5.5%), 6° (7.6%), 12° (21.1%), 20° (23.0%), 28° (26.9%), 41° (22.9%), 61° (23.1%). The raw ganglion cell density was multiplied by the proportion of displaced amacrine cells to obtain the corrected ganglion cell density between each pair of isoeccentricity rings (Fig.14).
Centrifugal displacement of ganglion cells in the macula
The primate fovea is highly specialized to maximize visual acuity. One of its principal adaptations is the centrifugal displacement of inner retinal elements, allowing photoreceptors unimpeded access to rays of light. The axons of photoreceptors fan out radially from the fovea to contact parafoveal bipolar cells, forming the Henle fiber layer of the macula. Bipolar cell axon terminals also are displaced centrifugally, contacting ganglion cells at an even greater eccentricity. The net result is that a photoreceptor located at 0.5°, for example, supplies a ganglion cell situated at 2.5°. Calculation of ganglion cell density in the retina cannot be related meaningfully to visual eccentricity without correcting for this displacement (Boycott et al., 1987; Perry and Cowey, 1988). To perform this correction, we measured the displacement of ganglion cells from their respective photoreceptors in the right retina from the animal whose RGCs we counted, one of the two monkeys used to make the mean retinotopic map (Adams and Horton, 2002).
Under general anesthesia and neuromuscular blockade, the right eye was photographed and landmarks were projected onto a tangent screen for eccentricity calibration. Then the right eye was enucleated, the anterior segment removed, and the fundus rephotographed in vitro (Fig. 15A). A ring pattern was superimposed onto the retina, calibrated with vascular landmarks identified before enucleation. Serial plastic 1 μm sections were cut along the vertical meridian through the fovea and stained with toluidine blue. Figure 15B shows a single section, cut along the vertical white line in the retinal photograph. Three major vessels are indicated by colored arrows. They could be identified clearly in the retinal photograph. To translate eccentricities on the tissue section from millimeters to degrees (thereby eliminating the problem of postmortem tissue shrinkage), we matched these blood vessels to their location in the calibrated retinal photograph.
We then measured the displacement of ganglion cells from their photoreceptors located at 1, 2, 4, and 8°. Applying the approach used by Schein (1988), we followed single fibers from photoreceptor nuclei through the Henle fiber layer to their pedicle in the outer plexiform layer (Fig. 15C). It was seldom possible to follow a single fiber the entire distance, but adjacent fibers followed a similar “grain,” allowing us to trace their trajectory accurately. Our measurements were made along only one dimension, the inferior vertical meridian. Perry and Cowey (1988) have shown that the Henle fibers are slightly longer in nasal retina than temporal retina, with those along the vertical meridian being intermediate in length. Thus our decision to sample along the vertical meridian was designed to provide a measurement of the mean displacement from photoreceptor to bipolar cell.
A small additional correction must be allowed for the radial displacement of bipolar cell axons onto ganglion cells (Polyak, 1941;Perry and Cowey, 1988; Schein, 1988). These axons could not be traced in our tissue cross sections. To estimate the displacement, we projected the long axis of the oval-shaped bipolar cells directly to the middle of the ganglion cell layer (Fig. 15C). This resulted in an additional 10% displacement. Perry and Cowey (1988)have measured this displacement in Golgi preparations of macaque retinal whole mounts and found that it accounts for 10% of the total photoreceptor-to-ganglion cell displacement, consistent with our estimate. Figure 16 shows the relationship between the eccentricity of photoreceptors and their corresponding ganglion cells in the squirrel monkey.
Relation between Ma and ganglion cell density as a function of eccentricity
We have calculated the density of ganglion cells in the retina (Fig. 14) and the displacement of ganglion cells from their photoreceptors (Fig. 16). By combining these data sets, we found that it is possible to correct for the displacement of ganglion cells to determine their functional density versus eccentricity. Via the principles elaborated by Schein (1988), the ganglion cells, in effect, must be shifted centrally into smaller annuli. The result is a measure of retinal output (effective ganglion cell density) at each eccentricity. This calculation was done in stages, as described below:
First, the centrifugal displacement from photoreceptors to their ganglion cells was applied to the isoeccentricity rings as follows: 1° ring → 3.10°, 2° ring → 4.05°, 4° ring → 5.45°, 8° ring → 9.25° (Fig. 16). This step translated the isoeccentricity rings from their optical position in the photoreceptor layer to their effective position in the ganglion cell layer. The rings beyond 8° were not expanded because they fell outside the region of centrifugal displacement of ganglion cells.
Second, the number of cells between the new, displaced locations of the isoeccentricity rings in the ganglion cell layer was calculated. This number was the product of the retinal surface area of each annulus (mm2) and the mean density of ganglion cells within that annulus (cells/mm2). Surface area was measured directly from the retina in Figure11D, and mean ganglion cell density was taken from the graph in Figure 14.
Third, the number of ganglion cells obtained in step 2 was divided by the area (mm2) of the corresponding annulus between undisplaced isoeccentricity rings, yielding the effective density of ganglion cells as a function of visual field eccentricity.
The result of these three steps is shown in Table1 and graphed in Figure17. The transformation produced a huge increase in the effective density of ganglion cells centrally. For example, only 488 ganglion cells were located physically within the central 1°, but 25,605 ganglion cells received their input from photoreceptors lying within this zone. The actual density of 5743 cells/mm2 (488 cells/0.085 mm2) increased to an effective density of 301,259 cells/mm2 (25,605/0.085 mm2). The 1–2° annulus contained 3665 ganglion cells, but its photoreceptors supplied 35,258 cells, a density increase from 14,316 to 137,727 cells/mm2. More peripheral annuli were affected less, because the displacement of photoreceptors from ganglion cells was smaller. The effective density of ganglion cells in the 4–8°and 8–16° annuli actually declined, because correction for the displacement transferred more cells out of the annuli than in.
There were just over a million ganglion cells in the squirrel monkey retina, close to the number in macaques and humans (Perry et al., 1984;Wässle et al., 1989; Curcio and Allen, 1990). Table 1 shows the percentage of the total ganglion cell population serving each annulus in the visual field. For example, ∼12% (2.5 + 3.4 + 6.4%) of the ganglion cells was devoted to representing photoreceptors located in the central 4° of the visual field. Approximately 30% of the ganglion cells served the central 8°.
Is the fraction of ganglion cells situated between each pair of isoeccentricity rings equal to the fraction of surface area corresponding to that annulus on the cortical map? If so, each RGC possesses a constant areal domain in striate cortex, regardless of its eccentricity. We measured the surface area (mm2) of each annulus in the mean map (Fig. 3). These data were compared with the effective ganglion cell number in each annulus (Table 1). Fig.18 shows the cumulative percentage of ganglion cells and cortical surface area as a function of increasing eccentricity. The percentage of ganglion cells within each annulus does not equal the percentage of surface area occupied by its corresponding annulus in the cortical map. For example, the cortical representation of the macula accounts for a greater percentage of surface area than one would predict from its percentage of ganglion cells. Specifically, the central 8° of the visual field occupies approximately one-half (48%) of striate cortex, yet the central 8° of retina contains only approximately one-third (30%) of the ganglion cells. The central 4° occupies 27% of the cortex but accounts for only 12% of the ganglion cells.
Discussion
Over the past century neuroscientists have toiled to compile an accurate retinotopic map of the primary visual cortex. This effort has been motivated by a conviction that a more exact map will lead to a better understanding of how the brain converts the external world into an internal representation for visual perception. In this paper we analyzed a mean retinotopic map derived from the representation of the vascular tree of the retina in the primary visual cortex of the squirrel monkey (Adams and Horton, 2002).
From this map several new insights have emerged. We derived mathematical expressions for cortical magnification factor along isopolar and isoeccentricity coordinates of the visual field map. Analysis of these functions showed that a modified complex logarithmic function (Schwartz, 1984) described accurately the cortical magnification along isopolar rays but overestimated the cortical magnification along isoeccentric rings. A map of local anisotropy (Mp/Me) in the cortical representation of the visual field was compiled, which showed a distinctive pattern of distortion related to systematic differences in Mp andMe along the vertical versus horizontal meridia. The CO patches in the cortex were back-projected onto the visual field to determine patch domain and density at each eccentricity. Finally, we sampled ganglion cell density at various eccentricities, using the retina from an animal whose cortices contributed to our generation of the mean map. Appropriate corrections were made for displaced amacrine cells and for centrifugal displacement of ganglion cells from their photoreceptors. Comparison of cortical magnification with ganglion cell density showed that the cortical territory corresponding to each ganglion cell is not constant at all eccentricities. Approximately twice the cortical tissue is devoted to the processing of visual information emanating from ganglion cells situated in the macula as compared with ganglion cells located in the periphery.
Accuracy of the retinotopic map
The precision of our map was obtained by matching a large number (ranging between 55 and 89) of corresponding retinal and cortical points in each hemisphere. The alignment was strengthened further by matching the curved line segments (usually blood vessels) connecting retinal points to their cortical representation. The mean map (Fig. 3) was based on a total retinal/cortical match of 295 points and 310 line segments. Because the map was essentially a metabolic imprint of the retina, there was no relative error in location between points. In contrast, maps derived from microelectrode sampling have the problems of residual eye movements, receptive field scatter, inaccuracy in plotting receptive field centers, and difficulty in reconstructing electrode tracks. These factors combine to create error in the location of each retinotopic point, and this error varies unpredictably from point to point.
Our retinotopic map contained several potential sources of error peculiar to our approach. The process of rendering the folded occipital lobe into a flat-mount creates distortion. This distortion is minimal, because the preparation of a flat-mount entails unfolding, with little stretching or compressing, but some local distortion must occur. We have no good way to measure this distortion. Ideally, one might place a flexible grid pattern over the cortical surface of the flat-mount, fold it back up to resemble its original shape, and see how the grid pattern is distorted. Short of performing such an experiment, we can say only that too much distortion simply causes the tissue to tear, creating a natural relieving cut. With practice, one can flatten V1 with no tears, indicating that the tissue is not being deformed excessively.
The calibration of each retina was achieved by projecting the location of ∼15 separate retinal landmarks onto a tangent screen, each with an error of ± 0.1°. This accuracy was achieved by using a fundus camera, a method that is more precise than a reversing ophthalmoscope. A correction could be made to convert linear measurements on the tangent screen to eccentricity in degrees. Unfortunately, the squirrel monkey eye is small and could not be photographed readily beyond 24°. For placement of the 32 and 50° isoeccentricity rings, we were forced to rely on Drasdo's schematic eye, scaled from the human to the squirrel monkey. Because this schematic eye may not be perfectly accurate and the optics of the squirrel monkey and human eye likely differ, this extrapolation probably introduced a small error. However, the peripheral retina is compressed into a small region of cortex, so a minor error in the location of these isoeccentricity rings could not have resulted in a displacement of more than a few hundred micrometers in the cortical map. Another potential error stemmed from our reliance on optical measurements to conclude that the monocular crescent begins at 72°. A behavioral measurement in an alert, trained monkey might have been more accurate.
The biggest source of error in our map resulted from the fact that blood vessels in the central 4° of the retina could not be discerned in the cortex. To fill this area, we inserted Cowey's map (1964) of the exposed cortical surface in the squirrel monkey. Fortuitously, this region corresponds to the central 4° of the visual hemifield, so it enabled us to localize the 1, 2, and 4° rings. Cowey's map was based on dozens of points sampled in a total of 47 monkeys. It is much easier to map accurately exposed cortex than buried cortex. For these reasons we believe his map is highly reliable. It was drawn, however, on a rear surface view of the intact occipital lobe. The shape of our flat-mount was slightly different, requiring a judgment about the best fit in transferring the map. This introduced uncertainty, especially in the location of the horizontal meridian. In our map it was placed in a position that did not divide the central 4° into equal upper and lower quadrants. Instead, the upper field representation occupied 122 mm2 (58%), and the lower field representation occupied 89 mm2 (42%). This discrepancy was caused by difficulty amalgamating our map with Cowey's and accounts for the higher anisotropy index of the lower field representation within the central 4° in Figure 6 as well as the higher density of patches shown within this sector in Figure 10.
Previous investigations have concluded that the lower field representation occupies more cortex than the upper field representation (Van Essen et al., 1984; Tootell et al., 1988). In humans the visual acuity, contrast sensitivity, temporal resolution, and reaction time are slightly better in the lower visual field (Low, 1943; Payne, 1965;Millodot and Lamont, 1974; Skrandies, 1985a,b). It is tempting to infer that these differences are attributable to greater cortical representation of the lower visual field. However, beyond 8° we found no difference in the areas occupied by the upper and lower visual fields. We found no evidence, either, for an enlarged representation of the lower visual field in the central 4°. Admittedly, our data could be biased because we may have mislocated the horizontal meridian in the central field slightly, as stipulated in the preceding paragraph. On the other hand, Van Essen and colleagues sampled large regions quite sparsely in their electrophysiological map [seeVan Essen et al. (1984), their Fig. 5], and the deoxyglucose images of Tootell and colleagues lack a wedge of lower operculum containing upper field representation [see Tootell et al. (1988), their Fig. 3]. Given these uncertainties, we believe that it is premature to conclude that any asymmetry exists in upper versus lower field representations.
Magnification factor
Calculation of cortical magnification across the visual field revealed how Ma,Me, andMp vary with eccentricity. We compared their functions with those predicted from the “complex log” map (for review, see Schwartz, 1984). The latter was found to fit well the average isopolar magnification factor (Fig. 4B). However, the cortex differed from the shape predicted by the complex log model (Fig. 5), because the upper and lower vertical meridia converged peripherally. As a result, the complex log model overestimated Ma andMe peripherally (Fig.4A,C).
Cortical anisotropy occurs whenMp andMe are unequal locally. A perfectly conformal map has no local anisotropy. The anisotropy map (Fig. 6) showed that Mp tends to be greater than Me along the vertical meridian and in the far periphery, owing to the shape of the cortex. Interestingly, the local variations ofMp andMe combine so thatMa remains approximately constant as a function of polar angle at each eccentricity (Fig. 7B).
Retinal ganglion cells
The central visual field occupies a huge portion of striate cortex. This phenomenon arises, to some extent, from the high density of ganglion cells in the macula. Our data show, however, that even more cortex is allocated to the macula than one would expect from its concentration of ganglion cells (Fig. 18). Thus there is not a constant relationship between ganglion cell density andMa from fovea to periphery. Instead, relatively more cortex is available for analysis of signals originating from central ganglion cells. This explains, in part, why macular cortex is labeled more faintly after intraocular [3H]proline injection than peripheral cortex [see Adams and Horton (2002), their Fig. 3].
Previous investigators (Malpeli and Baker, 1975; Myerson et al., 1977;Perry and Cowey, 1985; Silveira et al., 1989, 1993; Azzopardi and Cowey, 1993, 1996; Popovic and Sjöstrand, 2001) have arrived at this conclusion. It has been questioned by others (Albus, 1975; Schein and de Monasterio, 1987; Wässle et al., 1989, 1990), who have argued for the principle that ganglion cells supply the same cortical territory at all eccentricities. If our cortical map is accurate, the only important source of potential error lies in our measurement of ganglion cell density. An underestimate centrally, or an overestimate peripherally, might have led us astray. As mentioned in Materials and Methods, maintaining our retinal whole mount in a hydrated state made it possible to count accurately the small, densely packed ganglion cells in the macula. Thus we do not believe that we undercounted ganglion cells centrally. To avoid overestimating their density peripherally, we used GABA immunohistochemistry to identify displaced amacrine cells. We found that 22–27% of neurons in the peripheral retina were GABA-positive, in agreement with a previous study that found 20% of cells GABA-positive at 50° (Koontz et al., 1993). In contrast, Wässle et al. (1990) found that displaced amacrine cells outnumbered ganglion cells by a factor of 5 in the peripheral temporal retina. This surprisingly high ratio resulted in a low figure for peripheral ganglion cell density, leading the authors to conclude that there is no selective amplification of the macula in the cortex.
Correction for the centrifugal displacement of ganglion cells from photoreceptors introduces another potential source of error in calculating the allocation of cortical territory to ganglion cells. However, even a large error would not undermine our basic conclusion that central vision is over-represented. By 8° eccentricity the ganglion cells are nearly back in register with their photoreceptors. The central 8° contain 28% of the ganglion cells but correspond to 48% of the cortex (Table 1). Thus, ignoring the relative displacement of ganglion cells within the central 8°, one can see that this subpopulation of ganglion cells, taken as a whole, is over-represented in the cortex.
Approximately 10% of ganglion cells project to other targets than the lateral geniculate nucleus (Schiller and Malpeli, 1977; de Monasterio, 1978a,b; Perry and Cowey, 1984; Perry et al., 1984). If these nongeniculate targets place less emphasis on central vision, a greater percentage of peripheral ganglion cells may project to them. This source of potential error is small, however, because the overwhelming majority of ganglion cells projects to the lateral geniculate body. We recalculated our ganglion cells densities, assuming a “worst case scenario” (namely, that all nongeniculate projections come from peripheral retina), and found only a modest shift in the ratio between ganglion cell density and Ma, leaving our basic conclusion unchanged.
Strict scaling of cortical magnification factor to the density of RGCs is an appealing principle. However, like many principles of biology, it is true only to a first approximation. We conclude that the macula occupies a disproportionate amount of the surface area of V1. This phenomenon implies that the cortex handles processing of the macula in a special manner, by according it extra tissue for sensory analysis. This is another of the many specialized features of the macula that has evolved to optimize acuity.
Footnotes
This work was supported by the National Eye Institute (Grant RO1 EY10217 and Core Grant EY02162), That Man May See, The Bunter Fund, and the Lew W. Wasserman Merit Award from Research to Prevent Blindness. The California Regional Primate Research Center is supported by National Institutes of Health Base Grant RR00169. We thank the California Regional Primate Research Center (especially Jenny Short and David Robb). We also thank Davina Hocking and Irmgard Wood for help with tissue processing and Robin Troyer for assistance with animal care and surgery. Philip Sabes, Lawrence Sincich, and Lars Frisén provided valuable comments on this manuscript.
Correspondence should be addressed to Dr. Daniel L. Adams, Beckman Vision Center, University of California, San Francisco, 10 Kirkham Street, San Francisco, CA 94143-0730. E-mail: dadams{at}itsa.ucsf.edu.