Virtual tissue engineering and optic pathways: plotting the course of the axons in the retinal nerve fiber layer

As part of a larger project on virtual tissue engineering of the optic pathways, we describe the conditions that guide axons extending from the retina to the optic nerve head and formulate algorithms that meet such conditions. To find the entrance site on the optic nerve head of each axon, we challenge the fibers to comply with current models of axonal pathfinding.

First, we build a retinal map using a single type of retinal ganglion cell (RGC) using density functions from the literature. Dendritic arbors are equated to receptive fields. Shape and size of retinal surface and optic nerve head (ONH) are defined. A computer model relates each soma to the corresponding entry point of its axon into the optic disc. Weights are given to the heuristics that guide the preference entry order in the nerve.

Retinal ganglion cells from the area centralis saturate the temporal section of the disc. Retinal ganglion cells temporal to the area centralis curve their paths surrounding the fovea; some of these cells enter the disc centrally rather than peripherally. Nasal regions of the disc receive mixed axons from the far periphery of the temporal hemiretina, together with axons from the nasal half. The model plots the course of the axon using Bezier curves and compares them with clinical data, for a coincidence level of 86% or higher.

Our model is able to simulate basic data of the early optic pathways including certain singularities and to mimic mechanisms operating during development, such as timing and fasciculation.