Human foveal cone photoreceptor topography and its dependence on eye length
Abstract
We provide the first measures of foveal cone density as a function of axial length in living eyes and discuss the physical and visual implications of our findings. We used a new generation Adaptive Optics Scanning Laser Ophthalmoscope to image cones at and near the fovea in 28 eyes of 16 subjects. Cone density and other metrics were computed in units of visual angle and linear retinal units. The foveal cone mosaic in longer eyes is expanded at the fovea, but not in proportion to eye length. Despite retinal stretching (decrease in cones/mm2), myopes generally have a higher angular sampling density (increase in cones/deg2) in and around the fovea compared to emmetropes, offering the potential for better visual acuity. Reports of deficits in best-corrected foveal vision in myopes compared to emmetropes cannot be explained by increased spacing between photoreceptors caused by retinal stretching during myopic progression.
https://doi.org/10.7554/eLife.47148.001eLife digest
The human eye has many different parts that enable sight. The retina is the light sensitive tissue at the back of the eye, and it contains the fovea, the region that provides the clearest vision. Light must be focused on the retina to create images, a feat achieved by the transparent parts at the front of the eye called the cornea and lens. These parts of the eye are called the optics.
Between birth and adulthood, significant changes take place in the eye. Most noticeably, the distance between the optics and the fovea grows by about seven millimeters. Cells called cone photoreceptors, which provide light sensitivity, migrate and pack into the fovea. Finally, the eye’s optics adjust to maintain a sharp focus. At the same time, the brain is learning how to process inputs from the eyes to generate mental images that realistically correspond to the physical world around it. The development of the eye is fascinating in its complexity, but for more and more people, the process does not go as expected.
Specifically, a growing proportion of the population has eyes that are too long. This means that, for light reflected by far away objects, the eye’s optics form an image in front of the retina instead of on it. As a result, images of distant objects cannot be seen clearly, a condition known as myopia or nearsightedness. Researchers have also discovered that nearsighted people see less clearly than those who do not use glasses, even when given a sharp image to examine at close range. It has been hypothesized that these deficits result from stretching of the retina as the eye becomes bigger.
Until recently, testing this hypothesis by looking at cone photoreceptors directly in the eye was impossible. This is because the optics of all eyes have small imperfections that distort the light passing through them, including any light used to take high resolution microscopic images of the fovea. This hurdle can be overcome using adaptive optics, which means adding a deformable mirror to the instrument being used to image the eye that can adjust to correct the distortion.
Wang et al. use a new generation Adaptive Optics Scanning Laser Ophthalmoscope to check the density of cones at the fovea in relation to the size of the eye. They show that although the center of the fovea has fewer cones when the eye is bigger, this effect is more than offset because the longer eye increases magnification. So, if a near-sighted person wearing contacts and someone who does not need glasses stood side-by-side admiring the full moon, the near-sighted person would most likely have more cones sampling the image and should therefore have a higher resolution view.
These findings rule out reductions in the density of cone photoreceptors as the cause or effect of visual deficits associated with near-sightedness, adding to the understanding of this common condition.
https://doi.org/10.7554/eLife.47148.002Introduction
There has been a rapid increase in prevalence of myopia, of all magnitudes, in the period between 1971–1972 and 1999–2004 (Vitale, 2009). Across sub-populations grouped by race, ethnicity and gender, several studies report axial length of the eye to be the primary variable related to myopia (González Blanco et al., 2008; He et al., 2015; Iyamu et al., 2011). Increased axial length is associated with retinal stretching and thinning of posterior segment layers and the choroid (Fujiwara et al., 2009; Harb et al., 2015) and is associated with sight-threatening, often irreversible pathologies of the retina (Morgan et al., 2012; Verkicharla et al., 2015). Even without any detectable pathology, the structural changes associated with eye growth ought to have functional consequences for vision.
What do we know about functional deficits in myopia?
One might expect that eye growth would stretch the photoreceptor layer and would increase the spacing between cones, causing a longer eye to more coarsely sample an image relative to a shorter eye. However the situation is not that simple; the axial elongation associated with eye growth is accompanied by magnification of the retinal image (Strang et al., 1998). If the enlargement of the retinal image exactly matched the stretching of the cone mosaic, then eyes of different lengths would sample the visual field similarly. In fact, in large scale studies, myopes generally attain reasonably good visual acuity with optical correction (He et al., 2004; Jong et al., 2018).
However, more careful inspection reveals that myopes generally (6 out of 9 studies) have poorer angular resolution and have uniformly (3 out of 3 studies) poorer retinal resolution. Table 1 summarizes published results from psychophysical foveal tasks.
Most notably, Atchison et al. (2006) and Coletta and Watson (2006) show clear deficits in retinal resolution (cyc/mm) with increasing myopia using interferometric methods which bypass the optics of the eye and Rossi et al. (2007) show significant deficits in angular resolution (cyc/deg) in low myopes, even after using adaptive optics to correct for optical blur. All studies that find myopic visual deficits implicate retinal stretching as a possible cause, but what is actually happening structurally at the foveal center during myopic progression is not known. Therefore, the aim of the current study is to more carefully investigate how the length of the eye affects cone density at and near the foveal center.
Models for how photoreceptors change with eye growth
Two types of cone densities will be discussed in this study. Linear density quantifies how many cones are within a fixed area, in square mm, and serves as a way to evaluate physical retinal stretching caused by eye growth. Angular density quantifies how many cones are within one degree visual angle, (the visual angle on the retina is measured from the secondary nodal point of the eye). Angular density serves as a way to evaluate the visual implications of eye growth as it governs the sampling resolution of the eye.
Figure 1 illustrates three models, along the lines of Strang et al. (1998), of how photoreceptor structure might be affected by myopic eye growth. In the first model, called the global expansion model, the retina is proportionally stretched with increasing axial length - cones are more spaced out in longer eyes - and linear density decreases with eye length. Assuming that the secondary nodal point remains at a fixed position relative to the anterior segment, the number of cones within a fixed angular area will remain constant. Therefore, angular cone density will be constant with eye length. In the second model, called the equatorial stretching model, the posterior retina simply moves axially further from the anterior segment of the eye so that the linear density does not change with eye length. Since the retina is moving further from the secondary nodal point, more cones will fall within a fixed angular area and the angular cone density will increase with eye length. The final model, called the over-development model, describes a structural photoreceptor change that mimics the changes that occur during development (Springer and Hendrickson, 2004) whereby the photoreceptors continue to migrate towards the fovea as the eye grows. In this scenario, longer eyes will show both increased linear cone density and an even steeper increase in angular cone density. The model is motivated by observations of increased linear cone density in the foveas of marmosets that underwent lens-induced eye growth (Troilo, 1998).
Previous studies of cone spacing with axial length
The most definitive studies of cone spacing as a function of axial length are done through direct imaging of the retina – wherein sharp images of the cones are enabled through the use of adaptive optics, a set of technologies that actively compensate the blur caused by aberrations of the eye (Liang et al., 1997). Combined with confocal scanning laser ophthalmoscopy (Webb et al., 1987), adaptive optics offers the highest contrast en face images of the foveal photoreceptor mosaic ever recorded in vivo (Dubra et al., 2011; Roorda et al., 2002).
Despite continued advances in image quality, previous studies investigating cone packing and eye length have not made their measurements at the foveal center, the most important region for spatial vision but the most difficult to image owing to the small size of photoreceptors. There are a number of studies on cone packing and eye length (Chui et al., 2008; Elsner et al., 2017; Kitaguchi et al., 2007; Li et al., 2010; Obata and Yanagi, 2014; Park et al., 2013) and here we summarize the published results that are most relevant to our study. Chui et al. (2008) investigated angular and linear cone density at 1 mm and three degrees eccentricity. They found a significant decrease (p<0.05) in linear cone density as a function of eye length at 1 mm (which, by angular distance, is closer to the fovea in a longer eye than in a shorter eye) in all directions except in the nasal retina. They found that the angular cone density at three degrees (which, by linear distance, is closer to the fovea, in a shorter eye than in a longer eye) increased with eye length, but the trends were not significant. Li et al. (2010) made similar measures, but closer to the fovea (from 0.10 mm to 0.30 mm eccentricity). They found that linear cone density decreased with eye length, but the trends were not significant at the smallest eccentricities (0.1 and 0.2 mm). When the data were plotted in angular units and angular distance from the fovea, they found that angular cone density trended toward an increase with eye length but none of the trends were significant. A more recent study measured peak cone densities in the fovea as well as axial length for 22 eyes of 22 subjects (Wilk et al., 2017) but they did not plot peak cone density as a function of axial length, as it was not the aim of their study. We plotted the data they provided in their paper and found that the linear cone density at the foveal center dropped significantly with increases in axial length, similar to what was found by Li et al. (2010) and Chui et al. (2008), but the angular cone density had no dependency on eye length. Summary plots from previous literature are shown in Figure 2AB.
Wilk et al. (2017)’s data were consistent with a global expansion model and Li et al. (2010) and Chui et al. (2008)’s data only leaned toward a model that falls between the global expansion and equatorial stretching models. If the trends found by Li et al. (2010) and Chui et al. (2008) near the fovea were to extend to the foveal center, then myopes would have higher foveal photoreceptor sampling resolution with a consequent potential for better performance on visual tasks compared to emmetropes. As such, the simplest explanation for visual deficits in myopes – increased separation between cones caused by retinal stretching – would have to be ruled out.
With the improvements in resolution of adaptive optics ophthalmoscopes, imaging the smallest cones at the foveal center is now possible in many eyes, enabling a definitive analysis of the cone density at the fovea as a function of eye length.
Results
The experiments were approved by the University of California, Berkeley Committee for the Protection of Human Subjects. All subjects provided informed consent prior to any experimental procedures. Subjects self-reported their eye health so that only healthy individuals with no ocular conditions were included in the study. All eyes were dilated and cyclopleged with 1% Tropicamide and 2.5% Phenylephrine before imaging. We report data from 28 eyes of 16 subjects with a wide range of refractive error and axial length. Age, sex and ethnicity are listed on Table 2.
Biometry data
All the biometric measures used to convert angular dimensions to linear retinal dimensions are listed on Table 2. The strong correlation between refractive error and eye length (p<0.0001) indicates that the myopia in this cohort was predominantly as a result of axial length.
Imaging data
Images of the foveal region, the preferred retinal locus for fixation (PRL) and the fixation stability were recorded with an adaptive optics scanning laser ophthalmoscope (see Materials and methods). The image of one subject (10003L) is shown in Figure 3A. All the cones were resolved with our imaging system. The scatter plot indicates the scatter plot of fixation over the course of a 10 s video. Figure 3B shows the same image with all cones labeled and a color-coded overlay indicating the density. 16,184 labeled cones are shown on the figure. The point of maximum density is indicated by the black cross and the PRL is indicated by the yellow ellipse (best fit to the scatter plot locations in Figure 3A). This eye has a peak linear density of 204,020 cones/mm², and a peak angular density of 15,851 cones/deg². Cone density plots in linear and angular units for all eyes are shown on Figure 3—figure supplements 1 and 2. Original images and a list of the cone locations for each can be downloaded from the Dryad Digital Repository https://datadryad.org/review?doi=doi:10.5061/dryad.nh0fp1b.
Figure 4 shows the linear cone density as a function of linear eccentricity, where the average linear cone density was computed in 25-micron wide annuli centered around the point of peak density. Figure 4—figure supplement 1 are plots of the average cone densities in linear and angular units as a function of eccentricity in microns and arcminutes. Figure 4—figure supplement 2 are plots of the average linear and angular cone densities as a function of horizontal and vertical direction.
-
Figure 4—source data 1
- https://doi.org/10.7554/eLife.47148.015
In order to show the trends of density with axial length Figure 5A and B plot linear and angular cone density as a function of axial length where the colors indicate different eccentricity - red to purple indicate distance from the from fovea towards more parafoveal locations. Figure 5A reveals that peak linear density decreases significantly with axial length and the trend persists and remains significant from the fovea out to 100 microns eccentricity. Axial length accounts for 39% of the variance in the changes in linear cone density. Figure 5B shows the opposite trends when plotted in angular units. Peak angular density increases significantly with axial length and the trend persists and remains significant out to 40 arcminutes eccentricity. Axial length accounts for 32% of the variance in the changes in angular cone density. The plots clearly indicate that although stretching does occur (Figure 5A) it is not a simple global expansion and longer eyes have higher sampling density. The trends hold at and around fovea with statistical significance.
-
Figure 5—source data 1
- https://doi.org/10.7554/eLife.47148.017
A more relevant measure of the impact of eye length on vision is how the angular cone density changes at the PRL, which is often displaced from the location of peak cone density (Li et al., 2010; Putnam et al., 2005; Wilk et al., 2017). If, for example, longer eyes had more displaced PRLs then that could diminish, or even reverse, the trend of increased angular density with eye length reported in Figure 5B. We found that the average displacement between PRL and maximum cone density was 6.1 arcminutes and 30.4 microns. There was no significant linear relationship found between PRL displacement in either angular or linear units vs. axial length. Therefore, the PRL was not more displaced in myopes than in emmetropes from the point of peak cone density. Plots of the cone density at the PRL with axial length show the same trend at the PRL as at the point of maximum cone density (Figure 6A and B).
-
Figure 6—source data 1
- https://doi.org/10.7554/eLife.47148.019
Finally, we explored whether fixational eye movements might have a dependency on axial length. Fixation stability around the PRL had an average standard deviation of 3.94 arcminutes and 19.84 microns. The average area of the best fitting ellipse containing ~68% of the points in the scatterplot (defined as the bivariate contour ellipse area, or BCEA) was 50.7 square arcminutes and 1303 square microns. The plot of BCEA in square microns vs. axial length vs. showed a trend that approached significance (p=0.0596) (Figure 7A), but when we plotted BCEA in square arcminutes vs. axial length, the trend was no longer apparent (p=0.364) (Figure 7B). In other words, if there is any increase in fixational eye movements in microns, it is just a symptom of having a longer eye.
-
Figure 7—source data 1
- https://doi.org/10.7554/eLife.47148.021
Discussion
In this paper, we measure the cone density at and near the foveal center and investigate how it changes as a function of axial length. This is the first comprehensive study of cones in living eyes at the foveal center, the area solely responsible for a human’s fine spatial vision. Our results show that although some expansion does occur (linear cone density decreases with axial length) the angular sampling resolution actually increases, on average, with axial length. Prior to this study, the relationships between cone density and axial length were only made outside of the fovea, the closest being 0.1 mm, or 0.3 degrees (Li et al., 2010). Although an eccentricity of 0.3 degrees might seem close, it is noted that the cone density drops precipitously just outside of the location of peak density (Curcio et al., 1990) as does human vision (Poletti et al., 2013) (Rossi and Roorda, 2010b). There are other factors that govern peak cone density, however; eye length accounts for anywhere between 31% and 39% of the variance in cone density at the foveal center.
Our finding that the slopes of cone density vs. axial length are in opposite directions when plotted in linear (negative slope) and angular (positive slope) units, supports an eye growth model that lies between the global expansion model and an equatorial stretching model. Previous studies from our lab (Li et al., 2010) and also from Chui et al. (2008) leaned in the same direction. None of the cone density studies provide insight into the reasons why the photoreceptor density would behave this way with eye growth, but the results do align with other observations reported in the literature. Specifically, Atchison et al. (2004) used magnetic resonance imaging and found that eyeball dimensions in axial myopes are variable but are generally larger in all directions with a weak tendency to be preferentially greater in the axial direction. Their reported eye growth patterns lie between that illustrated for the global expansion and equatorial stretching models in Figure 1.
Our results differ from Wilk et al. (2017) whose data support a global expansion model (i.e. there is no detectable change in angular cone density with axial length; Figure 2B). But it is important to point out that their study did not set out to address the same question and the number of subjects with long axial lengths was disproportionately low.
Our results also differ from Troilo (1998) who studied retinal cell topography in a marmoset animal myopia model. Higher cone packing densities were observed in the experimentally enlarged eyes compared to normal eyes in the fovea. Their result followed the overdevelopment model, which is the reason why we included it as one of the possible outcomes of our study. In fact, the overdevelopment model is an extension of Springer’s model of development (Springer and Hendrickson, 2004), which offers a biomechanical explanation for how cone packing increases at the foveal center in a developing eye. While our data do not support the overdevelopment model, it does not preclude the existence of biomechanical factors working in opposition to simple global expansion.
The fact that angular cone density (visual sampling resolution) increases with eye length (myopia), at the peak density and at the PRL, means that poorer performance by myopes on resolution tasks cannot be explained by a decrease in photoreceptor sampling. The deficit must arise at a post-receptoral level.
Low-level causes for myopic visual deficits might arise from differences in the connectivity between cones and ganglion cells. Atchison et al. (2006) suggested that abnormal eye growth may be associated with a loss of ganglion cells. Alternately, if ganglion cells pool signals from multiple cones, then they will impose the retinal sampling limit and reduce certain aspects of visual performance (acuity, for example). Recent electron microscopy studies of a human fovea have revealed extensive convergence and divergence connections between photoreceptors and ganglion cells, albeit in an eye from an individual who was born prematurely (Dacey, 2018). These discoveries challenge our current understanding of neural connectivity in the foveal center and force us to consider the possibility of interindividual differences in foveal cone wiring. More experiments are necessary to explore these ideas.
To explain why low myopes did not perform as well on an acuity task as emmetropes, even after correction or bypassing of high order aberrations, Rossi et al. (2007) and Coletta and Watson (2006) both raised the possibility that myopes might have become desensitized to high frequency information (low level myopic amblyopia) as a result of having less exposure to a high contrast visual environment. In this case, it might be possible to train myopes to take advantage of their higher sampling resolution, but one myope in a follow up study by Rossi and Roorda (2010a) never reached the acuity levels of emmetropes in the same study.
Comparisons with previous studies
Peak cone densities
Curcio et al. (1990) measured spatial density of cones and rods in eight explanted whole-mounted human retinas. They found a large range of peak foveal cone densities with an average of 199,000 cones/mm2. Their density measurements were made over a window size of 29 × 45 microns. Our range of peak cone densities measured over a 10 arcmin (~50 micron) sampling window was 118,491 to 204,020 with an average of 163,242 cones/mm2. When we averaged the peak cone density over a circular aperture of 7.5 arcminutes which was more similar to the sampling window that Curcio et al. (1990) used to compute density, we measured peak linear cone densities ranging from 123,611 to 208,606 with an average of 166,854 cones/mm2. See Materials and methods for more discussion on the implications of changing the size of the cone sampling window. Zhang et al. (2015) reported an average peak density of 168,162 cones/mm2 in 40 eyes although they used a much smaller 5 × 5 micron sampling window to measure the peak. Wilk et al. (2017) reported an average peak density of 145,900 cones/mm2 in 22 eyes using a 37 × 37 micron sampling window and Li et al. (2010) reported an average peak density of 150,412 cones/mm2 in four eyes over a sampling window encompassing 150 cones (approximately 37 micron diameter at the foveal center). All reports of cone densities from adaptive optics studies in living eyes are lower than reports from histology. Two possible reasons for this are (i) the excised tissue in Curcio et al. (1990) underwent more shrinkage than estimated or (ii) the adaptive optics reports are subject to selection bias, where individuals with the highest angular cone densities might have been excluded because the image were less well resolved rendering the cones images too difficult to label with confidence. In our study, we attempted to image 73 eyes from 46 subjects and only succeeded in resolving cones across a sufficiently large region at and around the fovea in 28 of them. The reason the images from 45 eyes were not analyzed was due to poor or inconsistent image quality arising from a number of factors: Images from four eyes (three subjects) were not analyzed because their refractive errors were too high (all above –8D) and we ran into the limits of the deformable mirror’s dynamic range. Images from 18 eyes (13 subjects) that were taken early on in the study were not analyzed because the optics of AOSLO were not tuned well enough to resolve foveal cones. Images from four eyes (two subjects) were not analyzed because of uncorrectable image degradation caused by keratoconus and corneal scarring. Images from two eyes (one subject) were not analyzed because of excessive aberrations caused by an orthokeratology refractive correction. The cause of poor or inconsistent image quality among the remaining 17 eyes were varied, including ocular surface dryness, excessive eye motion and small pupils. The average refractive error among these remaining 17 eyes was about the same as the successful eyes.
Anisotropic density distribution
Like Curcio et al. (1990) and Zhang et al. (2015) we found steeper drops in cone density in the superior and inferior directions compared to the nasal and temporal directions. Plots of density along the two cardinal directions are shown on Figure 4—figure supplements 2.
PRL displacements
The distance of the PRL from the foveal center for our study (mean 30.4 microns; range 5.7–83.7; n = 28) roughly agrees with those of Wilk et al. (2017) (mean 63 microns; range 20–263; n = 22), Li et al. (2010) (mean 34 microns; range 3–92; n = 18) and Putnam et al. (2005) (mean 17; range 11–23; n = 5). The differences in cone density between the peak and the PRL were small and the trends (Figures 5 and 6) persisted at both locations.
Spatial vision estimates
The cone array imposes the first retinal sampling limit to human spatial vision (MacLeod et al., 1992; Williams, 1985) and the photoreceptor row-to-row spacing (assuming an hexagonal packing structure) imposes the maximum frequencies that can be relayed to later stages without aliasing. We can compute the sampling limit and estimate the cone center-to-center spacing using the following formulas:
For the densities reported here, the potential spatial frequency resolution limits range from 59.1 to 74.01 cyc/deg (average: 64.9 cyc/deg) at the peak density and 50.4 to 71.9 cyc/deg (average: 62.5 cyc/deg) at the PRL. These correspond to potential acuities ranging from 20/11.9 to 20/8.1 (based on the primary spatial frequency of the three bars of a Snellen E). The cone frequency cut-offs are higher than almost all the interferometric acuity limits reported by Coletta and Watson (2006), even for the emmetropic subjects. The acuities are, however, in the range of those measured from emmetropic subjects after adaptive optics correction by Rossi et al. (2007). The cone center-to-center spacing ranges from 0.59 to 0.47 arcminutes at the peak density and from 0.69 to 0.48 arcminutes at the PRL. A direct comparison of foveal structure and function for each of our subjects was not the scope of this study but will be the topic of future investigation.
Measuring structure and function of cone photoreceptors at the foveal center – the most important region of the human retina – has been one of the more challenging endeavors in vision science. Fortunately, the latest generation of adaptive optics ophthalmoscopes are making it easier and are facilitating new discoveries within this retinal region. The pattern of how cone density changes with eye growth lands somewhere between the global expansion and equatorial stretching models. The cone mosaic in longer eyes is expanded at the fovea, but not in proportion to eye length. Despite retinal stretching, myopes generally have a higher angular sampling density in and around the fovea compared to emmetropes. Reports of reduced best-corrected central visual acuity in myopes compared to emmetropes cannot be explained by decreased photoreceptor density caused by retinal stretching during myopic progression.
Materials and methods
Foveal imaging
Request a detailed protocolWe used our latest generation adaptive optics scanning laser ophthalmoscope (AOSLO) for foveal imaging. The system used a mirror-based, out-of-plane optical design (Dubra et al., 2011), and employed a deformable mirror with a continuous membrane surface and shaped with 97 actuators (DM97, ALPAO, Montbonnot-Saint-Martin, France). The system scans multiple wavelengths simultaneously. Each wavelength was drawn from the same broadband supercontinuum source (SuperK EXTREME, NKT Photonics, Birkerod, Denmark) using a custom-built fiber coupler. Wave aberrations were measured with a custom-built Shack Hartmann wavefront sensor using the 940 nm channel. Images were recorded using the 680 nm channel. 512 × 512 pixel videos were recorded over a 0.9 × 0.9 degree square field for an average sampling resolution of 9.48 pixels per arcminute. Eye alignment and head stabilization was achieved by using either a bite bar or a chin rest with temple pads. At least one 10 s video was recorded at the fovea and at eight more locations where the subjects were instructed to fixate on the corners and sides of the raster, to image an entire foveal region spanning about 1.8 × 1.8 degrees. In order to ensure the best possible focus of the foveal cones, multiple videos were taken over a range of 0.05 D defocus steps to find the sharpest foveal cones. Focus steps were generated by adding a focus shape onto the deformable mirror. Online stabilization and registration algorithms were used to facilitate rapid feedback on the image quality.
Locating the Preferred Retinal Locus of Fixation (PRL)
Request a detailed protocolSteady fixation was achieved at the fovea center by having the subjects fixate on a dark, circular, blinking dot with a diameter of 3.16 arcminutes (30 pixels) in the center of the AOSLO scanning raster. The fixation target was generated by modulating the same 680 nm scanning beam used for imaging and, as such, the target’s location was encoded directly into each frame of the video (Poonja et al., 2005). A scatter plot of the positions of the blinking dot relative to the retina was generated and was fit with a bivariate ellipse using free online MATLAB scripts downloaded from http://www.visiondummy.com/wp-content/uploads/2014/04/error_ellipse.m. The bivariate contour ellipse area (BCEA), which is the area of the best-fitting ellipse encompassing 68% of the points in the scatterplot (Castet and Crossland, 2012) was used to quantify the fixation stability (Figure 7) and the exact location of the PRL within the imaged cone mosaic (Table 2, Figure 3—figure supplements 1 and 2).
Image processing and analysis
Request a detailed protocolHigh quality images were generated from the recorded videos offline using custom software (MATLAB, The MathWorks, Inc, Natick, MA) to measure and correct for distortions caused by eye movements (Stevenson and Roorda, 2005). An open source version of the code is posted at https://github.com/lowvisionresearch/ReVAS (Agaoglu et al., 2018; copy archived at https://github.com/elifesciences-publications/ReVAS). Poor-quality frames were manually excluded and registered frames were averaged into a single high signal-to-noise image. The processed images were stitched together (Photoshop; Adobe Systems, Inc, Mountain View, CA) to create an approximately 1.8-degree montage of the foveal cone mosaic.
We used custom software to identify and label individual cones in the AO retinal images. The program allows the user to select a region of interest and manually add and delete cone labels. A combination of both manual and automated methods (Li and Roorda, 2007) were used to identify cone locations as the current version of the program does not adequately recognize cones in the foveal center where they are dim and smaller (Li et al., 2010). All the cone coordinates were selected and reviewed by two of the authors. In some cases cones were too dim to be seen but there was only a gap in the mosaic (Bruce et al., 2015). If a space that might have been occupied by a cone was dim or dark, we would assume it was a cone and mark its location. We rationalize this for two reasons: First, if there is a gap in the mosaic, then it is likely that a cell is occupying that space, otherwise the adjacent cells would migrate to fill it in (Scoles et al., 2014). Second, in our experience and of others (Pallikaris et al., 2003), cones that appear dark in one visit, can often appear bright in the next. In other cases (uncommon) the contrast was low in some regions or there were interference artifacts in the images (Meadway and Sincich, 2018; Putnam et al., 2010), making the cone locations slightly ambiguous. In these instances, we made manual cone selections based on the assumption that the cones were all similar in size and close-packed into a nearly hexagonal array (Curcio et al., 1990).
Continuous density maps were generated by computing cone density within a circular cone sampling window of 10 arcminutes in diameter around every pixel location across the image. All cones whose centers fell on or within the boundary of the sampling window were counted. The size of the cone sampling windows was chosen for several reasons. First, it encompasses the typical range of eye motion during a typical fixation task. Second, the area was large enough to generate smooth maps, but small enough to resolve local changes. Smaller cone sampling windows showed slightly higher densities (see Discussion) but were also associated with more variability in the identification of the location of peak cone density. Changes in density with eccentricity were generated by computing the density in five-arcminute annuli surrounding the point of peak cone density. For linear density measures we used annuli with 25-micron widths.
Retinal Magnification Factor calculation
Request a detailed protocolThe exact angular dimensions of the AOSLO images were computed by imaging a calibrated model eye in the AOSLO system, but the conversion to linear dimensions on the retinal image requires additional measurements, since the dimensions of each eye governs the actual size of the image on its retina. The conversion from visual angle to retinal distance requires a measurement of the axial length of the eye and an estimation of the location of the secondary nodal point. We did paraxial ray tracing, described by Li et al. (2010), in a four-surface schematic eye model to estimate the location of the secondary nodal point. The corneal first surface radius of curvature, the anterior chamber depth and the axial length were measured for each subject with an IOL Master (Zeiss Meditec, Dublin, CA). The radius of the curvature of the back surface of the cornea was computed as 88.31% of the front surface as per the Gullstrand eye model (Bennett et al., 1994). The indices of refraction of the media and the radii of curvature of the front and back lens surface were taken from the Gullstrand schematic eye (Vojniković and Tamajo, 2013). Once determined, retinal image size is related to visual angle by the equation:
Where I is retinal image size, x is axial length, AN’ is the distance from the corneal apex to the eye’s second nodal point, and θ is the visual angle. As can be seen in Table 2, myopic eyes, which generally have longer focal lengths, have proportionally larger retinal images.
Other methods to compute the retinal magnification factor (RMF) might have been used. We did a paraxial ray trace to find the secondary nodal point using lens parameters from a 6-surface Gullstrand eye model, and found that the RMFs were essentially identical. When we computed the RMF using an equation proposed by Bennett et al. (1994) we found that the RMF was overestimated for hyperopic eyes and underestimated for myopic eyes. This discrepancy arises from the fact that Bennett’s equation is based on the assumption that the focal point of the eye’s optics coincides with the retina, which it does not (Li et al., 2010). Despite this known error, when we redid the analysis using RMFs computed using the Bennett et al. (1994) equation, we still found a significant trend of reducing linear density with axial length, but only at the position of peak cone density. Of course, all the angular density calculations are unaffected by differences in the RMF.
Statistical analysis
Request a detailed protocolGiven the trends of increased angular density as a function of axial length that Li et al. (2010) observed at the location closest to the fovea (slope = 531 cones/deg2 for each mm increase in axial length; standard deviation of the regression errors = 1377 cones/deg2), we estimated that data from approximately 32 eyes, evenly distributed across a range of axial lengths would be sufficient to show if there was a true effect at the fovea. The targeted number was computed using methods outlined by Dupont and Plummer (1998) implemented using free online software (Power and Sample Size Program Version 3.0, January 2009, downloaded from http://biostat.mc.vanderbilt.edu/wiki/Main/PowerSampleSize) with type one error probability of 0.05 and a power of 0.95.
All data collected in this study were analyzed using simple linear regression models in Excel. p-values for all linear regressions are reported and linear trendlines with p-values less than 0.05 are plotted as solid lines and p-values greater than 0.05 as dashed lines.
In this study, we included both the right and left eyes of some of the subjects. When we reran the statistics with a single eye from each subject (n = 16), the results were essentially the same and the same conclusions could be drawn.
Data availability
The following data are available for download at the Dryad Digital Repository: https://dx.doi.org/10.5061/dryad.nh0fp1b. 1. All original images and cone locations. 2. A table of scaling parameters (pixels per degree and retinal magnification factors) for each image. 3. A MATLAB script that can be used to plot cone locations on the original image.
-
Dryad Digital RepositoryData from: Human Foveal Cone Photoreceptor Topography and its Dependence on Eye Length.https://doi.org/10.5061/dryad.nh0fp1b
References
-
Eye shape in Emmetropia and myopiaInvestigative Opthalmology & Visual Science 45:3380–3386.https://doi.org/10.1167/iovs.04-0292
-
Neural and optical limits to visual performance in myopiaVision Research 46:3707–3722.https://doi.org/10.1016/j.visres.2006.05.005
-
Improvements on Littmann's method of determining the size of retinal features by fundus photographyGraefe's Archive for Clinical and Experimental Ophthalmology 232:361–367.https://doi.org/10.1007/BF00175988
-
Normal perceptual sensitivity arising from weakly reflective cone photoreceptorsInvestigative Opthalmology & Visual Science 56:4431–4438.https://doi.org/10.1167/iovs.15-16547
-
Quantifying eye stability during a fixation task: a review of definitions and methodsSeeing and Perceiving 25:449–469.https://doi.org/10.1163/187847611X620955
-
Retinal stretching limits peripheral visual acuity in myopiaVision Research 45:593–605.https://doi.org/10.1016/j.visres.2004.09.016
-
Individual Variations in human cone photoreceptor packing density: variations with refractive errorInvestigative Opthalmology & Visual Science 49:4679–4687.https://doi.org/10.1167/iovs.08-2135
-
Visual performance in high myopiaCurrent Eye Research 9:217–224.https://doi.org/10.3109/02713689009044516
-
Human photoreceptor topographyThe Journal of Comparative Neurology 292:497–523.https://doi.org/10.1002/cne.902920402
-
Discovering visual pathway origins in the center of the human foveola with connectomicsInvestigative Ophthalmology and Visual Science 59:EAbstract#14.
-
Noninvasive imaging of the human rod photoreceptor mosaic using a confocal adaptive optics scanning ophthalmoscopeBiomedical Optics Express 2:1864–1876.https://doi.org/10.1364/BOE.2.001864
-
Power and sample size calculations for studies involving linear regressionControlled Clinical Trials 19:589–601.https://doi.org/10.1016/S0197-2456(98)00037-3
-
Spatial contrast sensitivity of myopic subjectsVision Research 16:437–438.https://doi.org/10.1016/0042-6989(76)90214-5
-
Enhanced depth imaging optical coherence tomography of the choroid in highly myopic eyesAmerican Journal of Ophthalmology 148:445–450.https://doi.org/10.1016/j.ajo.2009.04.029
-
Axial length, corneal radius, and age of myopia onsetOptometry and Vision Science 85:89–96.https://doi.org/10.1097/OPX.0b013e3181622602
-
Choroidal thickness profiles in myopic eyes of young adults in the correction of myopia evaluation trial cohortAmerican Journal of Ophthalmology 160:62–71.https://doi.org/10.1016/j.ajo.2015.04.018
-
Refractive error and visual impairment in urban children in southern chinaInvestigative Opthalmology & Visual Science 45:793–799.https://doi.org/10.1167/iovs.03-1051
-
Altered visual sensitivity in axial high myopia: a local postreceptoral phenomenon?Investigative Opthalmology & Visual Science 47:3695–3702.https://doi.org/10.1167/iovs.05-1569
-
Reduced vision in highly myopic eyes without ocular pathology: the ZOC-BHVI high myopia studyClinical and Experimental Optometry 101:77–83.https://doi.org/10.1111/cxo.12563
-
In vivo measurements of cone photoreceptor spacing in myopic eyes from images obtained by an adaptive optics fundus CameraJapanese Journal of Ophthalmology 51:456–461.https://doi.org/10.1007/s10384-007-0477-7
-
Intersubject variability of foveal cone photoreceptor density in relation to eye lengthInvestigative Opthalmology & Visual Science 51:6858–6867.https://doi.org/10.1167/iovs.10-5499
-
Automated identification of cone photoreceptors in adaptive optics retinal imagesJournal of the Optical Society of America A 24:1358–1363.https://doi.org/10.1364/josaa.24.001358
-
Supernormal vision and high-resolution retinal imaging through adaptive opticsJournal of the Optical Society of America A 14:2884–2892.https://doi.org/10.1364/JOSAA.14.002884
-
Myopia and contrast sensitivity functionCurrent Eye Research 22:81–84.https://doi.org/10.1076/ceyr.22.2.81.5530
-
A visual nonlinearity fed by single conesVision Research 32:347–363.https://doi.org/10.1016/0042-6989(92)90144-8
-
Light propagation and capture in cone photoreceptorsBiomedical Optics Express 9:5543–5565.https://doi.org/10.1364/BOE.9.005543
-
The reflectance of single cones in the living human eyeInvestigative Opthalmology & Visual Science 44:4580–4592.https://doi.org/10.1167/iovs.03-0094
-
Dynamic visual stimulus presentation in an adaptive optics scanning laser ophthalmoscopeJournal of Refractive Surgery 21:S575–S580.https://doi.org/10.1364/fio.2004.fthv3
-
The locus of fixation and the foveal cone mosaicJournal of Vision 5:3–639.https://doi.org/10.1167/5.7.3
-
Adaptive optics scanning laser ophthalmoscopyOptics Express 10:405–412.https://doi.org/10.1364/OE.10.000405
-
The relationship between visual resolution and cone spacing in the human foveaNature Neuroscience 13:156–157.https://doi.org/10.1038/nn.2465
-
In vivo imaging of human cone photoreceptor inner segmentsInvestigative Opthalmology & Visual Science 55:4244–4251.https://doi.org/10.1167/iovs.14-14542
-
Correcting for miniature eye movements in high resolution scanning laser ophthalmoscopyOphthalmic Technologies XV 5688:145–151.
-
The effect of myopia on contrast thresholdsInvestigative Opthalmology & Visual Science 48:2371–2734.https://doi.org/10.1167/iovs.05-1377
-
High myopia does not affect contrast sensitivityCurrent Eye Research 5:635–640.https://doi.org/10.3109/02713688609015130
-
Vision Science and Its Applications: Technical Digest206–209, Changes in retinal morphology following experimentally induced myopia, Vision Science and Its Applications: Technical Digest, 1, Optical Society of America.
-
Current and predicted demographics of high myopia and an update of its associated pathological changesOphthalmic and Physiological Optics 35:465–475.https://doi.org/10.1111/opo.12238
-
Increased prevalence of myopia in the united states between 1971-1972 and 1999-2004Archives of Ophthalmology 127:1632–1639.https://doi.org/10.1001/archophthalmol.2009.303
-
Gullstrand’s optical schematic system of the eye-modified by Vojnikovic & TamajoCollegium Antropologicum 37:41–45.
-
Confocal scanning laser ophthalmoscopeApplied Optics 26:1492–1499.https://doi.org/10.1364/AO.26.001492
-
Aliasing in human foveal visionVision Research 25:195–205.https://doi.org/10.1016/0042-6989(85)90113-0
-
Variability in human cone topography assessed by adaptive optics scanning laser ophthalmoscopyAmerican Journal of Ophthalmology 160:290–300.https://doi.org/10.1016/j.ajo.2015.04.034
Article and author information
Author details
Funding
National Eye Institute (R01EY023591)
- Nicolas Bensaid
- Pavan Tiruveedhula
- Jianqiang Ma
- Austin Roorda
National Eye Institute (K08EY025010)
- Sowmya Ravikumar
National Eye Institute (T35EY007139)
- Yiyi Wang
National Eye Institute (P30EY003176)
- Austin Roorda
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Ethics
Human subjects: The work on human subjects was reviewed and approved by the University of California, Berkeley Institutional Review Board and is described in protocol # 2010-08-2006. All subjects provided informed consent prior to any experimental procedures.
Copyright
© 2019, Wang et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,984
- views
-
- 288
- downloads
-
- 79
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Neuroscience
Hypothalamic kisspeptin (Kiss1) neurons are vital for pubertal development and reproduction. Arcuate nucleus Kiss1 (Kiss1ARH) neurons are responsible for the pulsatile release of gonadotropin-releasing hormone (GnRH). In females, the behavior of Kiss1ARH neurons, expressing Kiss1, neurokinin B (NKB), and dynorphin (Dyn), varies throughout the ovarian cycle. Studies indicate that 17β-estradiol (E2) reduces peptide expression but increases Slc17a6 (Vglut2) mRNA and glutamate neurotransmission in these neurons, suggesting a shift from peptidergic to glutamatergic signaling. To investigate this shift, we combined transcriptomics, electrophysiology, and mathematical modeling. Our results demonstrate that E2 treatment upregulates the mRNA expression of voltage-activated calcium channels, elevating the whole-cell calcium current that contributes to high-frequency burst firing. Additionally, E2 treatment decreased the mRNA levels of canonical transient receptor potential (TPRC) 5 and G protein-coupled K+ (GIRK) channels. When Trpc5 channels in Kiss1ARH neurons were deleted using CRISPR/SaCas9, the slow excitatory postsynaptic potential was eliminated. Our data enabled us to formulate a biophysically realistic mathematical model of Kiss1ARH neurons, suggesting that E2 modifies ionic conductances in these neurons, enabling the transition from high-frequency synchronous firing through NKB-driven activation of TRPC5 channels to a short bursting mode facilitating glutamate release. In a low E2 milieu, synchronous firing of Kiss1ARH neurons drives pulsatile release of GnRH, while the transition to burst firing with high, preovulatory levels of E2 would facilitate the GnRH surge through its glutamatergic synaptic connection to preoptic Kiss1 neurons.
-
- Neuroscience
Specialized chemosensory signals elicit innate social behaviors in individuals of several vertebrate species, a process that is mediated via the accessory olfactory system (AOS). The AOS comprising the peripheral sensory vomeronasal organ has evolved elaborate molecular and cellular mechanisms to detect chemo signals. To gain insight into the cell types, developmental gene expression patterns, and functional differences amongst neurons, we performed single-cell transcriptomics of the mouse vomeronasal sensory epithelium. Our analysis reveals diverse cell types with gene expression patterns specific to each, which we made available as a searchable web resource accessed from https://www.scvnoexplorer.com. Pseudo-time developmental analysis indicates that neurons originating from common progenitors diverge in their gene expression during maturation with transient and persistent transcription factor expression at critical branch points. Comparative analysis across two of the major neuronal subtypes that express divergent GPCR families and the G-protein subunits Gnai2 or Gnao1, reveals significantly higher expression of endoplasmic reticulum (ER) associated genes within Gnao1 neurons. In addition, differences in ER content and prevalence of cubic membrane ER ultrastructure revealed by electron microscopy, indicate fundamental differences in ER function.