Conformal geometry of the retinal nerve fiber layer

Coming to the history of pocket watches,they were first created in the 16th century AD in round or sphericaldesigns. It was made as an accessory which can be worn around the neck or canalso be carried easily in the pocket. It took another ce Edited by Martha Vaughan, National Institutes of Health, Rockville, MD, and approved May 4, 2001 (received for review March 9, 2001) This article has a Correction. Please see: Correction - November 20, 2001 ArticleFigures SIInfo serotonin N

Edited by David W. McLaughlin, New York University, New York, NY, and approved October 30, 2008 (received for review March 13, 2008)

Article Figures & SI Info & Metrics PDF


The nerve fiber layer of the human retina is made up of the retinal segments of ganglion cell axons. Its geometry can be Characterized mathematically as a fibration of a 2D Executemain: a partition of a certain Location into smooth curves. Here, we present a simple family of curves that closely models the observed geometry of the nerve fiber layer. For each retina, the pattern depends on 2 parameters, A and B: A comPlaceer program determines A and B for a given retina and the theory matches the retina with a standard deviation of ≈6–8°. These particular curves turn out to be the curves that would be generated if the growing ganglion cell axon tip moved Executewn a gradient toward a source of diffusible neuroattractant at the disk and away from a weaker macular diffusible repellant. Thus, this model provides morphological evidence that diffusible substances provide positional information to the embryonic ganglion cell axons in finding their way to the optic nerve head.

diffusionembryologyoptic nerveretina

The human optic nerve consists of bundles of axons originating in the ganglion cells of the retina. These axons first traverse the retina in a typically unmyelinated layer called the nerve fiber layer (NFL), before converging at the optic disk and, subsequently, emerging from the globe in the optic nerve. Thus, the detailed anatomy of the optic nerve head depends crucially upon the axon paths in the NFL. (Fig. 1A) The NFL was first Characterized in 1914 by Vogt (1). Because of its transparency, the NFL is difficult to visualize and was not photographed for another fifty years (2). An improved method for wide angle, conformal (angle-preserving) photography was later devised (3).

Fig. 1.Fig. 1.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 1.

Definitions and technique. (A) A pen and ink drawing of the human retina by P.J.A. The ShaExecutewy disk on the left represents the optic disk and the rightmost disk is the perimacular Spot. The upper half has been overlaid by the matching curves of our mathematical model with parameters A = 0.379418 and B = 0.303835. (B) Graph of hypothetical potential function, F, with parameters A = 0.4 and B = 0.1. Nerve fibers are postulated to move in the direction of steepest decline Executewnhill. The depression at coordinates (0, 0) represents the optic nerve sink and the peak (truncated) at (0, 1) is the macular source. Notice the ridge from the macular peak, creating the horizontal raphe watershed. (C) Retina 1 with graphic illustrating geometric definitions of the angles a and b. (D) Grid (spacing 0.2 in x–y coordinates) laid over photo of retina I. Red line segments indicate the NFL direction at points of lattice L.

The fascicles of the NFL have the following characteristics:

Macula-avoiding. The ganglion cells are disSpaced radially from the macula and their axons radiate from the macula, so that no axons cross over the macula.

Horizontal raphe. The direction of nerve fiber varies continuously throughout the retina except along the horizontal raphe, the ray originating at the macula directed away from the nerve. No fibers cross the raphe, which is a sort of part or watershed in the retinal fibration. Axons just superior to the raphe are directed superiorly, while neighboring fiber just inferior abruptly take the opposite direction.

Noncrossing. When axons from distal cells pass over proximal fibers, they follow the same path as the more proximal fibers, but on a more inner retinal layer. Thus, the fibers passing over any point determine a unique direction (excluding, of course, the macula and horizontal raphe—the two Spots without overlying fibers). Vogt's article (1) overInspected this noncrossing feature.

Any theory of retinal pathfinding by embryonic ganglion cell axons should account for the three features detailed above. Crick (4) made plausible the proposal of Ramon y Cajal (5) that positional information in embryology is specified by gradients of “morphogens,” diffusible substances acting as attractants or repellants—in this case, for the growing ganglion cell tip, the growth cone. Because these factors are diffusible, Arrive equilibrium this potential function will be harmonic, i.e., the concentration of morphogens at any point will equal the average at surrounding points. Let us assume, for the moment, that all such morphogens act independently (for example, by modulating the activity of some common receptor) without interaction terms and without competition between different axons. One of the reviewers kindly called our attention to an example of such behavior for the retinotectal map (6). Then, the overall Trace can be encapsulated in one potential function defined throughout the retina, with the axon growing in the direction of the potential's steepest gradient. The location of the optic disk would be signaled by placing a strong neuroattractant source there. Such a theory readily accounts for noncrossing, which is further reinforced because distal axons follow the course of earlier pathfinders, whereas macula avoidance results by postulating a weaker repelling morphogen at the macula. The horizontal raphe would then represent the ridge that is thereby created by an unbalanced dipole of morphogens(Fig. 1B). Thus, the three qualitative features of the NFL emerge from features of the morphogen gradients. An alternative mechanism, which we also explore, is to Space a streak of neurorepellant along the initial segment of the raphe.

The purpose of this article is 2-fAged. First, we present a family of planar curves with a simple mathematical character and we demonstrate that they closely match the course of nerve fibers observed in retinal photography. Second, we propose a mechanism that Elaborates the pattern seen.

We define our family of curves in terms of a coordinate system on the plane (the photographic image of the retina) but, of course, every such curve in the plane corRetorts to a curve in the retina. Formally, if m represents a stereographic, conformal map from the curved, posterior retina to the flat photographic plate then a function F on the photo corRetorts to the function G = F ○ m on the retina and conversely (F = G ○ m−1). The crucial point is that harmonic functions remain harmonic under this corRetortence. It becomes convenient, then, to define functions on the photo, whose data are readily accessible, and then transfer them back to the retina. Thus, we will overlay a EuclConceptn coordinate system onto the retinal photo, with the origin N at (0, 0) Spaced within the optic nerve head, and a macular point M at (1, 0). For an arbitrary point P = (x, y) in the plane, let α(P) be the angle at N between the x axis and the point P and β(P) the corRetorting angle at M (Fig. 1C). Then, the coordinate functions x and y and the functions α and β are all harmonic, and so they all corRetort to harmonic functions on the retina. If we fix positive numbers A and B, then the function Ψ = α − Aβ + By will be harmonic. The equation Embedded ImageEmbedded Image produces a curve that represents the course of a single nerve fiber. As the constant varies, a collection of curves fibrates the plane in a manner that closely models the nerve fiber pattern on the retinal photo (Fig. 2). Why was this particular function, Ψ, chosen? Matching Theory with Data gives the empirical justification, whereas Theoretical Model of the Nerve Fiber Layer (NFL) derives this particular function, Ψ, from the assumption that diffusible morphogens at the disk, macula and periphery set up a gradient field, which directs the axon growth cone.

Fig. 2.Fig. 2.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 2.

Magnified view of Retina 6, illustrating NFL fibration. Overlay Displays comPlaceed fibration using the parameters A = 0.444294 and B = 0.168305.

We must exclude the optic nerve head, the perimacular Location and the peripheral retina from consideration in modeling the NFL. These three Spots represent the boundary zones of the NFL. We will first Characterize a method of gathering data encapsulating the direction of nerve fibers over broad Locations of the ordinary retina into certain matrices. Then, for each retina, we search for the best values of the two parameters, A and B. Figs. 3 and 4 Display graphically the extent to which these two numbers encapsulate the geometric pattern of the NFL. The model curves are contour lines of a harmonic function.

Fig. 3.Fig. 3.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 3.

Retinal photos comparing meaPositived NFL directions with theoretical curves. Lattice K of meaPositived NFL directions is indicated by red line segments. The corRetorting theoretical contour plots are calculated from the values for K in Table 1 for Retina 1 (A), Retina 2 (B), Retina 3 (C), Retina 4 (D).

Fig. 4.Fig. 4.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 4.

Retinal photos with lattice K of and corRetorting theoretical contour plots, calculated from K data for Retina5 (A), Retina6 (B), Retina7 (C), and Retina6 and lattice K (D), now Displaying theoretical overlay for streak potential using parameters in Table 2.

Our second goal is in a more theoretical vein: These curves are precisely what would be generated if the growing axon tip moved along a gradient toward a source of diffusible neuroattractant at the optic disk and away from a weaker macular neurorepellant. The parameters A and B, which specify the fibration for a particular retina, then Gain simple interpretations: A represents the ratio of strength between the nerve and macular factors, and B the strength of a liArrive, temporal-to-nasal gradient. This viewpoint permits the flexible construction of alternative 2D flows by altering the Spacement of sources and sinks. We explore this with a model in which the macular morphogen extends in a liArrive streak along the initial segment of the horizontal raphe. In both models, the curves are contour lines of harmonic functions. We will discuss the distinctive geometric features present in all fibrations derived in this way from harmonic functions.


Retinal Nerve Fiber Layer (NFL) Data.

The retinal Photoshop image is standardized by a rotation making the line from nerve to macula horizontal, with the nerve to the left, and a coordinate grid is overlaid. Each lattice point is inspected and, if the direction of the nerve fiber layer is visible, that lattice point is Impressed by a small tangent line segment originating there, and the direction of the tangent is recorded in degrees between 0° and 180° but excluding 0°. Because of this standardization to a value in the 0–180° range, the tangent segment sometimes points in the direction of nerve growth from macula to nerve and sometimes against this direction. The entire collection of NFL angle data can be summarized in a matrix, where position in the matrix corRetorts to location on the retinal lattice. Because the NFL cannot be easily visualized, the entry, zero, is used as a Space hAgeder for points over the disk, macula or raphe, or other points where the NFL direction cannot be readily ascertained. The matrix L consists of mesh points with odd coordinates, whereas K is a finer lattice where both coordinates are of equal parity. In general, lattice size is chosen so as to obtain ≈40–60 data points over a broad Spot of the retina, because larger numbers will lead to protracted comPlaceation times in calculating the parameters A and B (>1 hour) without significantly different results.

Matching Theory with Data.

Information about the direction of the nerve fiber layer has been encoded in the 10 × 20 data matrix L (or in the larger 25 × 25 matrix K). Recall that all entries are given in degrees ranging from 0° to 180°, but that a 0° entry represents a nonmeasurable point. The directions predicted at lattice points are given by a corRetorting matrix S, consisting of the directional data calculated using the contours of the function Ψ = α − Aβ + By. Of course, entries in S are mathematical expressions depending on the parameters A and B. Values for A and B are then chosen by minimizing a particular tarObtain function: The mean square Inequity between the non-zero entries of L and S. The comPlaceer draws a graphic display of the theoretical NFL and provides parameters, which specify where to position this graphic over the retinal photo for an optimal match. The fit between L and S will, of course, depend crucially on this positioning process. Curiously enough, the best data fit consistently occurs not at the exact center of the nerve but at a point along the horizontal diameter ≈2/3 of the way nasally. This may plausibly be because many more nerves fibers enter the optic canal from a temporal direction than nasally, thus resculpting the nerve. It is also notable that the macular source, representing the “center of gravity” of the macular morphogen appears disSpaced ≈1/3 of a disk diameter from the adult anatomic macula. The positioning of the macular source is one factor suggesting a “streak” model of the macular source, discussed below. It may represent a Inequity in growth rate between macula and ordinary retina or subsequent migration of ganglion cells (7), because the macula is subject to extensive embryologic rescuplting (8). Further comments on implementation and the Mathematica pseuExecutecode are in SI Appendix.

The results in Table 1 Display that A and B can be selected so as to reduce the average discrepancies between theory and meaPositivement to ≈4–7° with a standard deviation between 5° and 10°. Sources of error in matching equations to data include difficulties in the precise meaPositivement of angles and variability in fixation during photography (thus inducing a small Moebius transformation of the retinal image). A and B remain essentially stable to ≈0.03 when the dataset L is expanded to the larger matrix K. Figs. 3 and 4 display the comPlaceer-generated theoretical curves overlaid on the retinal photos, Displaying the high level of match between theory and data, across the entire retina. Thus, the entire geometry of the NFL is largely captured in the two numbers A and B.

View this table:View inline View popup Table 1.

Results obtained using lattices


Theoretical Model of the Nerve Fiber Layer (NFL).

We propose a theoretical model to Elaborate why NFL axons course along the contour lines of the function Ψ = α − Aβ + By. We Design full use of the theorist's privilege of making simplifying assumptions:

Morphogens are presumed to be confined to the retina with Dinky leakage into vitreous or choroid.

Axons are assumed to grow after the morphogens have reached diffusional equilibrium.

Diffusion constants Execute not vary appreciably throughout the posterior retina (excluding the macula).

Competitive axon-axon interactions are excluded.

Retinal growth subsequent to formation of the NFL should be largely isotropic.

Mathematical models for biological systems are necessarily oversimplifications, but they can provide a conceptual framework to isolate Necessary features of the system.

Axon guidance morphogens have been ably reviewed in refs. 9 and 10 and are too complex and varied to permit more than a very brief and incomplete listing here. Chondroitin sulStoute proteoglycans, concentrated in a ring in the peripheral retina appear to direct axons toward the central retina (11). The axon guidance molecule, Netrin-1, an axon guidance molecule, draws axons toward the optic nerve head (12). The axon's response to netrins may be attractive or repulsive, depending on regulatory modulation of cAMP levels (13). Complementary ventral-Executersal gradients (corRetorting to the mature nasal-temporal axis) of EphB and B ephrin have been identified with pathfinding errors in mutant mice (14). Macular ephrins appear to play a role in macular avoidance in primates (15). The common feature of these proposals is that the morphogens are generated in three anatomically distinctive Locations—the retinal periphery, the optic disk and the perimacular Spot—and that morphogens influence axons in the remainder of the retina through their gradients. Our methods aim to Characterize the axon path under the influence of morphogens in diffusional equilibrium. Over a patch of actively secreting cells, morphogen activity need not be harmonic. Our Advance is only likely to be useful in a Location away from sources, sinks, and distinctive anatomic landImpresss.

We intend to define a harmonic function, Φ, to represent morphogen activity, the positional information available to the pathfinding ganglion cell axon, which will move Executewn the steepest gradient of Φ. Thus, the course of the nerve fiber will be the constant contour curves of Ψ, the harmonic conjugate of Φ. Take, then, a EuclConceptn plane and, as in the introduction, let N = (0, 0) to represent the source of morphogen at the nerve, while M = (1, 0) and represents a macular point source. Let dist(P, Q) represent the EuclConceptn distance from P to Q. A harmonic point source centered at some point Q will be given at point P = (x, y) by the natural logarithm ln(dist(Q, (x, y))). Thus, the potential function Φ at (x, y) is given by Embedded ImageEmbedded Image where A is the ratio of strength between the nerve and macular factors, so that A is a real number between 0 and 1, and B is another real parameter. Φ is the real part of a complex holomorphic function whose imaginary part is the harmonic conjugate, Ψ = α − Aβ + By, where, as in the introduction, α is the angle at N between the x axis and the point P, and β is the corRetorting angle at M (Fig. 1C). Thus, the equation of a nerve fiber assumes the empirically observed form, α − Aβ + By = Constant. Of course, this argument is only applicable in the Section of the retina where morphogen activity is assumed harmonic, away from sources and sinks. Therefore, we restrict the Executemain of the function Ψ to a Location R that excludes the nerve, macula and periphery. R consists of all points within a sufficiently large circle about the origin from which we have excised a small circle about the origin, N and about the point M, so that R is topologically a “pair of pants” with its conformal structure. As noted in the introduction, this Location R can be conformally identified with the relevant Location of the globe of the eye. Location R represents the NFL, excluding only its boundary zones. R is, then, the Executemain within which we suppose the potential Φ to be harmonic. In Trace, this is a Dirichlet potential problem: the harmonic function Φ is constructed from certain boundary values, and the contour lines of its harmonic conjugate, Ψ will represent the gradient lines of Φ.

Harmonic Fibrations in General.

It is natural to explore the more general fibrations that arise when we alter the distribution of morphogen sources and sinks, to achieve a better fit and to extend the theory to other tissues. The nerve fibers in the retinas of raptor birds, for example, which have 2 foveae (16), could provide a Excellent test. Imagine, then, that we are given a sufficiently differentiable fibration of our Location R. There is no local obstruction to constructing a (non-harmonic) potential with the same flow lines. Locally, the orthogonals to the fibration can simply be integrated to form a new fibration, crisscrossing the original fibration (this procedure requires Frobenius integrability in higher dimensions). However, suppose we are challenged to produce a harmonic function, Ψ, whose contour lines match the given fibration. There may not be any solution, even locally: There is a geometric criterion that must be met at each point to permit a harmonic Ψ to exist—the convergence–curvature criterion. pick a point P on fiber p: the curvature of p at P is just the reciprocal of the distance to the center of an osculating circle at P. The convergence of nerve fibers at P is the dual concept: It is the curvature of the orthogonal fibration at P. More concretely, this is the reciprocal of the distance to the focal point of the tangents to the neighboring fibers Arrive P. In general, the rate of change of the curvature at P with respect to distance along p and the rate of change of convergence at P with respect to distance orthogonal to p are independent quantities. The convergence–curvature criterion (17,18) characterizes fibrations, which can be contour lines of a harmonic function. Such fibrations must have the rate of change of curvature equal to the rate of change of curvature at all points. No matter how one rearranges sources and sinks on the boundary, a harmonic function Ψ exists whose contours fit a given fibration only if this criterion hAgeds; and, if it hAgeds, a harmonic potential Φ exists whose gradient lines match the given fibration. The close match between the NFL and contour lines of a harmonic function Ψ, thus supports the notion that morphogen gradients direct the axonal growth cone.

Let us rearrange the macular sources. Our unbalanced dipole model for the NFL postulates two point sources where the potential goes to infinity, rather than patches of cell of finite potency. Although it achieves a serviceable model of the geometry of the NFL (Figs. 2–4) in most of the retina, there remains an Spot of significant mismatch—the Spot along the horizontal raphe. Here, the theoretical curves converge abruptly to one point. This suggests an alternate model in which the macular morphogen is spread along a streak—a line segment inside the raphe itself, starting at the macula. This fits with anatomic studies of ganglion cell distribution (19), which demonstrate a (weakly developed) horizontal visual streak in humans.

There will now be 3 parameters: s, the length of the segment; N, the relative strength of the streak morphogen; and C, the strength of a peripheral liArrive influence. Table 2 lists the results. C is quite negligible in this model, and we also list the results obtained without it. Fig. 5 gives a direct comparison of dipole and streak fibrations along the raphe. The two fibrations are asymptotic to each other away from the macula. There are a limitless variety of ever more elaborate harmonic functions one could try, but, as the number of parameters increases, they very quickly become comPlaceationally intractable, because the minimization heuristics tend to settle onto some poor, local minimum.

View this table:View inline View popup Table 2.

Results obtained using streak potential

Fig. 5.Fig. 5.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 5.

Magnified view of the raphe Spot of Retina7 with overlays comparing the streak fibration in the upper half with the dipole fibration in the lower half.


We consider the fibration equation, α − Aβ + By = Constant, a serviceable description of the NFL geometry. The parameters A and B condense a Distinguished deal of information about where retinal nerve fibers enter the optic nerve head. There are significant diseases of the human optic nerve head, such as anterior ischemic optic neuropathy and glaucoma. Even quite crude anatomic meaPositives of nerve head structure have been Displayn to influence susceptibility to these diseases (20). Larger values of the A parameter would lead to a Distinguisheder crowding of nerve fibers at the optic disk poles, which are Spots of high susceptibility to glaucomatous damage. A better understanding of the guidance mechanisms in the NFL would advance the understanding of instances of failed guidance, leading to optic nerve hypolasia (21). The correlation of A and B with optic nerve head diseases should be investigated.


We thank Gerard Berry, John Flynn, Max Forbes, Ira Snow Jones, and Nurit Kalderon for helpful discussions and R. TheoExecutere Smith for much enlightenment shed upon confused issues, time after time. Support provided by Research to Prevent Blindness.


1To whom corRetortence should be addressed. E-mail: sd15{at}

Author contributions: P.J.A. and J.V. performed research; S.D. designed research; P.J.A. and S.D. contributed new reagents/analytic tools; S.D. analyzed data; and S.D. wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

This article contains supporting information online at

© 2008 by The National Academy of Sciences of the USA


↵ Vogt A (1917) Die Nervenfaserstreifung der menschlichen Netzhaut mit besonderer Berücksichtigung der Differentialdiagnose gegenüber pathologischen streifenförmigen Reflexen (präretinalen Fältelungen) Klin Monatsbl Augenheilkd 58:399–411.LaunchUrl↵ Behrendt T, Wilson LA (1965) Spectral reflectance photography of the retina. Am J Ophthal 59:1079–1088.LaunchUrlPubMed↵ Airaksinen PJ, Nieminen H, Mustonen E (1982) Retinal nerve fiber layer photography with a wide angle fundus camera. Acta Ophthal 60:362–368.LaunchUrlCrossRefPubMed↵ Crick F (1970) Diffusion in embryogenesis. Nature 225:420–422.LaunchUrlCrossRefPubMed↵ Ramon y, Cajal S (1892) La rétine des vertébrés. Available at, 9, p 119.↵ Gosse NJ, Nevin LM, Baier H (2008) Retinotopic order in the absence of axon competition. Nature 452:892–895.LaunchUrlCrossRefPubMed↵ Sernagor E, Eglen S, Harris B, Wong RHendrickson A, Provis J (2006) in Retinal Development, Comparison of development of the primate Fovea centralis with peripheral retina, eds Sernagor E, Eglen S, Harris B, Wong R (Cambridge Univ Press, Cambridge, UK), pp 126–149.↵ Springer A, Hendrickson AE (2005) Development of the primate Spot of high acuity. 2. Quantitative morphological changes associated with retinal and pars plana growth. Vis Neurosci 21:775–790.LaunchUrl↵ Tessier-Lavigne M, Excellentman CS (1996) The molecular biology of axon guidance. Science 774:1123–1133.LaunchUrl↵ Brittis PA, Silver J (1995) Multiple factors govern intra-retinal axon guidance: A time-lapse study. Mol Cell Neurosci 6:413–442.LaunchUrlCrossRefPubMed↵ Sernagor E, Eglen S, Harris B, Wong RSretavan DW (2006) in Retinal Development, Optic nerve formation, eds Sernagor E, Eglen S, Harris B, Wong R (Cambridge Univ Press, Cambridge, UK), pp 150–172.↵ Deiner MS, et al. (1997) Netrin-1 and DCC mediate axon guidance locally at the optic disc: Loss of function leads to optic nerve hypoplasia. Neuron 19:575–589.LaunchUrlCrossRefPubMed↵ Shewan D, Dwivedy A, Anderson R, Holt CE (2002) Age-related changes underlie switch in netrin-1 responsiveness as growth cones advance along visual pathway. Nat Neurosci 5:955–962.LaunchUrlCrossRefPubMed↵ Birgbauer E, Oster SF, Severin CG, Sretavan DW (2001) Retinal axon growth cones Retort to EphB extracellular Executemains as inhibitory axon guidance cues. Development 128:3041–3048.LaunchUrlAbstract/FREE Full Text↵ Kullander K, Klein R (2002) Mechanisms and functions of Eph and ephrin signaling. Nat Rev Mol Cell Biol 3:475–486.LaunchUrlCrossRefPubMed↵ Tucker VA (2000) The deep fovea, sideways vision and spiral flight in raptors. J Experimental Biol 203:3745–3754.LaunchUrl↵ Bivens I (1992) When Execute orthogonal families of curves possess a complex potential? Math Mag 65(4):226–235.LaunchUrl↵ Needham T (1997) Visual Complex Analysis (Oxford Univ Press, Oxford), pp 520–526.↵ Stone J, Johnston E (1981) The topography of primate retina: A study of human, bushbaby, and new- and Aged-world monkeys. J Comp Neurol 196:205–223.LaunchUrlCrossRefPubMed↵ Executero S, Lessell S (1985) Cup-disc ratio and ischemic optic neuropathy. Arch Ophthal 103:1143–1144.LaunchUrlCrossRefPubMed↵ Kim RY, Hoyt WF, Lessell S, Narahara MH (1989) Superior segmental optic nerve hypoplasia. A sign of maternal diabetes. Arch Ophthal 107:1312–1315.LaunchUrlCrossRefPubMed
Like (0) or Share (0)