Predictive Mathematical Models for the Spread and Treatment of Hyperoxia-induced Photoreceptor Degeneration in Retinitis Pigmentosa

of hyperoxia-induced photoreceptor degeneration in retinitis pigmentosa

T he group of retinal diseases known collectively as retinitis pigmentosa (RP) is the world's most common inherited retinal degeneration. 1Affecting approximately 1 in 4000 individuals worldwide, 2,3 RP causes a progressive loss of visual function, culminating in complete blindness by approximately the age of 60. 3 While a number of hypotheses have been proposed as to the disease etiology, at present no treatments are clinically available to halt or reverse disease progression. 4][10][11][12] These patches then expand and coalesce, leading to a number of distinctive spatio-temporal patterns of retinal degeneration, correlating with visual field loss patterns.Grover et al. 13 classified visual field loss into three classes (Patterns 1, 2, and 3), which may be subdivided into six subclasses: Patterns 1A, 1B, 2A, 2B, 2C, and 3 (see Fig. 2).Pattern 1A involves a constriction of the peripheral visual field, while 1B also includes a parafoveal/perifoveal ring scotoma (blind spot).Pattern 2 begins with either a loss of nasal or superior nasal vision (2A), out from which an arcuate (bow shaped) scotoma winds inferiorly through the midperiphery to the temporal side, or with a superior temporal (2C) or complete superior (2B) restriction, out from which an arcuate scotoma winds inferiorly through the midperiphery to the nasal side.In each case, only the central and inferior visual field is preserved.Lastly, in Pattern 3, a complete or partial midperipheral ring scotoma develops initially.The scotoma then expands into either the superior or inferior periphery, leaving a central island of vision, together with a U-or n-shaped peripheral visual field, the arms of the 'U' or 'n' subsequently retracting.For all patterns, central vision is generally the best preserved, being lost in the final stage of the disease.][16][17][18][19][20][21][22][23] While the initial loss of photoreceptors may be attributed to genetic mutations, 2,3 the cause of the expansion of degenerate patches remains a matter of speculation.Five main hypotheses have been proposed to explain this phenomenon: the oxygen toxicity, [5][6][7] rod trophic factor, [24][25][26][27][28][29][30] toxic substance, 31 microglia, 32 and metabolic dysregulation (mTOR) 33,34 hypotheses.The oxygen toxicity hypothesis suggests that the mutationinduced loss of photoreceptors causes an increase in the retinal oxygen tension, such that conditions become toxic (hyperoxic) for the remaining photoreceptors.This results in a positive feedback loop, as subsequent photoreceptor loss further increases oxygen levels, which are maintained due to the inability of the choroid (the retina's chief oxygen supply) to autoregulate. 5,35,36yperoxia-induced photoreceptor degeneration may be combated using a combination of antioxidant and trophic factor treatments. Antioxidants neutralize the excess reactive oxygen species generated under hyperoxia, 37 while trophic factors are thought to increase photoreceptor resistance to apoptosis under hyperoxic conditions. 38,391][42][43][44] Both treatments have been shown to reduce retinal degeneration in RP. 35,38,39,43,[45][46][47] Several mathematical models have been developed to describe the rod trophic factor, [48][49][50][51][52][53][54] mTOR, 54 and toxic substance [55][56][57] hypotheses (see Roberts et al. 58 for a detailed review of the current state-of-the-art in retinal modeling, including models of RP).In an earlier study we developed a one-dimensional (1D) model (fovea to ora serrata) to describe retinal degeneration under the oxygen toxicity hypothesis. 59o the best of our knowledge, this was the first model to consider this hypothesis, or to account for the spatial distribution of rods and cones.In this paper, we extend this model to two dimensions (2D), incorporating the azimuthal dimension, so that we can consider scenarios where radial symmetry is broken.Comparing our modeling predictions with the patterns of retinal degeneration that can be inferred from Grover et al., 13 we judge the strengths and weaknesses of the oxygen toxicity hypothesis by our models' ability to recapitulate the characteristic patterns associated with RP.We also use our models to predict the effects of treatment with antioxidants and trophic factors under the oxygen toxicity hypothesis.

METHODS
In what follows we provide a succinct description of our three mathematical models.The first model considers only hyperoxia-induced photoreceptor degeneration, while the second and third include mutation-induced rod and cone degeneration  respectively.For further details on the assumptions behind the models see Roberts et al. 59 The equations were solved using the finite element method, using the Portable, Extensible Toolkit for Scientific Computation libraries (available in the public domain, www.mcs.anl.gov/petsc/)together with piecewise linear basis functions (see Roberts 60 for further details).
Our models focus on the outer retina, which we define as the region between Bruch's membrane (BM)-where the retina abuts the inner capillary layer of the choroid known as the choriocapillaris (CC)-and the photoreceptor inner segments (ISs) inclusive.This region also contains the retinal pigment epithelium and the photoreceptor outer segments (OSs).RP typically occurs as a rod-cone dystrophy, in which rod photoreceptors are affected earlier and more severely than cone photoreceptors; however, cone-rod varieties also exist. 2,3ur models account for the heterogeneous distribution of photoreceptors across the retina, distinguishing between rod and cone photoreceptors when considering the effects of mutation-induced rod or cone loss in the rod-cone and conerod forms of RP.
Adopting a spherical polar coordinate system, we model the eye as a spherical cap (see Fig. 3), orientating the coordinate axes such that the center of the fovea lies at 08 eccentricity (h ¼ 0 [rad], noting that all angles are given in radians).We model the retina in 2D: polar, h È [0, H] (rad), and azimuthal, / È [0, 2p) (rad), neglecting the radial dimension, r (m), because the retinal thickness (80-320 [lm], see Webvision, available in the pubic domain, http://webvision.med.utah.edu/) is two orders of magnitude smaller than its radius of curvature (» 1.2 [cm] 61 ).The ora serrata is assumed to lie at h ¼ H » 1.33 (rad), based on data supplied by Curcio et al. 62 (noting that, unlike the retina depicted in Fig. 3a, the retinas measured by Curcio et al. 62 do not extend beyond the equator at h ¼ p/2 [rad] in the sector of the eye from which these measurements are taken).We assume the portion of the retina we are modeling comprises an effectively depth-averaged section between BM and the outer tips of the photoreceptor ISs.
Oxygen diffuses freely through the retina, is consumed by retinal tissue and exchanged with the CC, while photoreceptors maintain their healthy density or regenerate lost biomass under normoxia (healthy oxygen levels) and degenerate under hyperoxic conditions.We define partial differential equations for the retinal oxygen concentration, c(h,/,t) (mol m -3 ), and the photoreceptor density, p(h,/,t) (photoreceptors m -2 ), over time, t (s), as follows where ]c/]t and ]p/]t are the rate of change of the oxygen concentration and photoreceptor density respectively over time.The variable p(h,/,t) can also be thought of as representing OS biomass density (multiplied by the appropriate scaling factor to maintain dimensional consistency), this interpretation being more appropriate in those cases where regrowth occurs since new photoreceptors cannot be regrown in the adult retina (see Table 1 for a summary of the model variables).We note that D (m 2 s -1 ) is the (constant 63 ) diffusivity of oxygen and R (m) is the radial position of the retina.Oxygen consumption is assumed to have Michaelis-Menten kinetics (an increasing, saturating function of the retinal oxygen concentration), where Q (mol s -1 [tissue unit] -1 ) is the maximum rate of oxygen consumption under light-adapted conditions (when consumption is lowest and the risk of hyperoxic damage is highest), 64 assuming that rods and cones have the same oxygen demand, 65 and c (mol m -3 ) is the oxygen concentration at which oxygen uptake is half maximal.The parameter a (m -1 ) is the ratio of unit surface area to unit volume (ensuring dimensional consistency), c ch (mol m -3 ) is the oxygen concentration in the CC, b (m s -1 ) is the effective permeability of the CC vessels and BM to oxygen, h (m -1 ) is the capillary surface area per unit volume of tissue and d (s -1 ) is the rate of hyperoxia-induced photoreceptor degeneration.We remark that the rates of oxygen supply and uptake are large enough to maintain a heterogeneous oxygen profile despite the smoothing effect of diffusion.It has been necessary to reduce the values of both Q and b by a factor of 10 from those used in Roberts et al. 59 to render the simulations computationally feasible.This results in a faster rate of propagation of hyperoxic degeneration, but does not alter the spatial pattern that develops, which is our focus here (see Roberts 60 for details).OS biomass regrowth is assumed to occur logistically with intrinsic growth rate l (s -1 ) and carrying capacity pðhÞ (photoreceptors m -2 ), where This functional form was chosen because of its qualitative similarity to the healthy human photoreceptor distribution, where the first two terms capture the cone profile and the last term captures the rod profile.The values for the parameters B 1 (photoreceptors m -2 ), B 2 (photoreceptors m -2 ), B 3 (photoreceptors m -2 rad -1 ), b 1 (rad -1 ), b 2 (rad -1 ), and b 3 (rad -1 ) were obtained by using the Matlab curve fitting toolbox (Math-Works, Natick, MA, USA) to fit Equation 3 to the mean of eight experimentally measured healthy adult photoreceptor distributions across the temporal horizontal meridian, using data provided by Curcio et al. 62 (see Table 2).We note that, for simplicity, we assume that the photoreceptor distribution depends only upon h and hence is axisymmetric about the zaxis (see Fig. 3).The functions k1 ðcÞ (dimensionless) and k2 ðcÞ (dimensionless) that appear in Equation 2 are defined as wherein c crit (mol m -3 ) represents the 'hyperoxic threshold', above which photoreceptors degenerate and below which they remain healthy or regenerate biomass.The parameter S c (rad m 3 mol -1 ) determines the sharpness of the transition between normoxia and hyperoxia.We use hyperbolic (tanh) functions here, rather than the Heaviside step functions used in our earlier work, 59 in order to make the governing equations smooth, improving the convergence of the numerical scheme.This remains a biologically realistic choice for these terms.We close the system by imposing the following boundary and initial conditions: The boundary conditions are zero-flux, such that there is no net gain or loss of oxygen to or from the retina along its boundary at the ora serrata.The initial oxygen distribution, c init (h,/) (mol m -3 ), is the steady-state oxygen concentration (the distribution to which oxygen settles over time) corresponding to the initial photoreceptor profile (that is, the steady-state solution to Equations 1 and 5 with p(h,/) ¼ p init (h,/) [photoreceptors m -2 ]).The function 0 F(h,/) 1 (dimensionless) is used to remove a patch of photoreceptors and takes one of the following three forms: which correspond to a healthy retina, the removal of an annulus (ring) of photoreceptors and the removal of a disc of photoreceptors respectively (see Supplementary Fig. S1 for graphs depicting each of these initial conditions).Degenerate discs represent the patchy loss observed in the early stages of RP, while annuli represent a later disease stage.The parameters S (dimensionless) and Ŝ (rad -1 ) determine the sharpness of the edge of an annulus or disc respectively, while h 1 (rad) and h 2 (rad) are the eccentricities of the inner and outer boundaries of the annulus.The disc is centered at ðh; /Þ ¼ ðh c ; /c Þ (rad) (where we write the azimuthal position with a 'hat' to distinguish it from the rate of mutation-induced cone degeneration below) with radius governed by w (rad).The sin 2 (h) term is required in order to preserve arc length in the azimuthal direction.We choose hyperbolic (tanh) functions, rather than step functions (as used in our earlier work) 59 in order to smooth the transition between the healthy and degenerate retina, so that a coarser finite element mesh can be used to resolve the initial wavefront of degeneration.See Table 2 for parameter values.

Mutation-induced Rod Degeneration
Modifying our model to account for mutation-induced rod degeneration (corresponding to the rod-cone form of RP), we denote rod density by p r (h,/,t) (photoreceptors m -2 ) and cone density by p c (h,/,t) (photoreceptors m -2 ), to obtain the following system: where / r (s -1 ) is the rate of mutation-induced rod degeneration and rods are assumed to be incapable of recovering biomass.We close Equations 8 to 10 by imposing zero-flux oxygen boundary conditions at h ¼ H (rad) (given by Equation 5) and the following initial conditions: where F(h,/) is defined in Equation 7and the healthy rod and cone distributions, pr ðhÞ(photoreceptors m -2 ) and pc ðhÞ (photoreceptors m -2 ), are given by As before, the initial oxygen concentration, c init (h,/), is the steady-state oxygen concentration corresponding to the initial photoreceptor profile (that is, the steady-state solution to Equations 5 and 8 with p r ðh;

Mutation-induced Cone Degeneration
Similarly, we may account for mutation-induced cone degeneration (corresponding to the cone-rod form of RP), obtaining Equation 8 together with ]p r ]t ¼ lp r 1 À p r pr ðhÞ k1 ðcÞ À dp r k2 ðcÞ; where / c (s -1 ) is the rate of mutation-induced cone degeneration.The value chosen for / c is twice that used for / r in the present study (see Table 2), and twice that used for / r and / c in Roberts et al., 59 in order to reduce the computation time, while remaining within the range of biologically realistic values for this parameter, the spatial pattern of degeneration Where two sets of units are stated, the first set is consistent with the models, while the second, in brackets, is that specified in the references.
remaining the same as for the lower value.We close the system by imposing the boundary and initial conditions given by Equations 5 and 11 respectively.

Treatment With Antioxidants or Trophic Factors
For each of the cases above, we may apply treatment in the form of antioxidants or trophic factors.These treatments are assumed to increase the hyperoxic threshold, such that c crit becomes a function of time, t.We set c crit ðtÞ ¼ c crit1 þ c crit2 Hðt À t crit Þ, where c crit1 (mol m -3 ), c crit2 (mol m -3 ) and t crit (s) are positive constants and H is a Heaviside step function, such that c crit ðtÞ ¼ c crit1 for t < t crit and c crit ðtÞ ¼ c crit1 þ c crit2 for t ‡ t crit .We use a Heaviside function rather than a hyperbolic (tanh) function here because smoothing is not required in this case for numerical stability.

Degenerate Annulus
Our model predicts that the photoreceptor distribution may evolve in one of four ways following the removal of an annulus (ring) of photoreceptors, in the absence of mutation-induced rod or cone degeneration (see Supplementary Fig. S2).The degenerate annulus may (1) remain stationary (Supplementary Fig. S2a), (2) spread centrally (Supplementary Fig. S2b), (3) spread peripherally (Supplementary Fig. S2c), or (4) expand both centrally and peripherally (Supplementary Fig. S2d).The evolution of the annulus is determined by the initial positions of the inner and outer boundaries of the annulus.Boundaries remain stationary where they border regions with a high local photoreceptor density and propagate outward where they border regions with a low local photoreceptor density.We note that these results are consistent with earlier results from a simpler 1D model (see Roberts et al. 59 ).

Degenerate Disc
Simulations predict that, in the absence of mutation-induced rod or cone degeneration, the evolution of the photoreceptor distribution following the removal of a disc of photoreceptors is determined mainly by the disc's eccentricity, h c (where h c ¼ 0 [rad] at the foveal center and h c ¼ H [rad] at the ora serrata), and, to a lesser extent, by its radius, governed by w (larger values of w corresponding to a larger radius, see Fig. 4a).(See Roberts et al. 59 for an asymptotic analysis showing that the qualitative behavior of the model will be maintained within realistic parameter ranges due to the separation of scales between various parameter groupings.)As with the case of a degenerate annulus, a degenerate disc expands where its boundary borders regions in which the local photoreceptor density is low, and remains stationary otherwise.Photoreceptors recover for 0.3 3 H h c 0.4 3 H (rad) (Fig. 4d), while photoreceptor loss within the disc, which is only partial initially, is consolidated (that is photoreceptor loss becomes complete in this region) for h c ¼ 0.7 3 H (rad) (Fig. 4e) for all disc radii examined.The disc expands into the central retina for h c ¼ 0 (rad) (Fig. 4b) and around the retinal periphery for h c > 0.7 3 H (rad) (Fig. 4f) for all disc radii examined.Degeneration is consolidated within the disc for w ¼ 0.025 3 H (rad) and w ¼ 0.0125 3 H (rad) when h c ¼ 0.1 3 H (rad); however, when w ¼ 0.05 3 H (rad), degeneration spreads into the central retina (Fig. 4c).Lastly, the disc either completely recovers, consolidates completely or recovers in some areas while consolidating in others for h c ¼ 0.2 3 H (rad), hc ¼ 0.5 3 H (rad) and hc ¼ 0.6 3 H (rad), depending upon its radius.

Mutation-induced Rod Degeneration
Our model simulations suggest that mutation-induced rod degeneration will cause hyperoxic photoreceptor degeneration to initiate first at the ora serrata and later in the parafoveal/ perifoveal region (see Fig. 5).Waves of hyperoxic degeneration propagate centrally from the ora serrata and peripherally from the macula, meeting in the midperiphery, such that only the foveal cones remain at steady-state.  2 for parameter values.
Simulations further predict that mutation-induced rod degeneration can cause expansion of degenerate regions such as annuli and discs that would otherwise have remained fixed in size.Comparison of the results presented in Figures 4e and 6 illustrates this behavior: a degenerate disc which remains fixed in size in the absence of mutation-induced rod degeneration (see Fig. 4e) expands in size when mutation-induced rod degeneration is active (see Fig. 6).

Mutation-induced Cone Degeneration
Simulations suggest that mutation-induced cone degeneration will stimulate hyperoxia-driven degeneration of photoreceptors in the central retina (see Fig. 7).Degeneration initiates in the fovea and spreads peripherally, coming to a halt as it enters the midperiphery.

Treatment With Antioxidants or Trophic Factors
Our models predict that treatment with antioxidants and/or trophic factors may either delay, halt, or partially reverse hyperoxia-induced photoreceptor degeneration depending upon the strength (c crit2 ) and timing (t crit ) of treatment, the extent of the degenerate region(s) when treatment is applied and the presence or absence of mutation-induced rod or cone degeneration (see Fig. 8).In the absence of mutation-induced rod or cone degeneration, treatment may halt degeneration (Figs.8a, 8c) or enable partial recovery (Fig. 8b).By contrast, in the presence of mutation-induced cone degeneration treatment may at best halt (Fig. 8e) degeneration and in the presence of mutation-induced rod degeneration delay it (Fig. 8d).

DISCUSSION
Mathematical modeling enables us to isolate biological mechanisms in a manner that would be difficult or impossible experimentally.In this paper we have used spatially resolved 2D models to predict those patterns of retinal degeneration that could be produced if oxygen toxicity were the sole cause for the spread of photoreceptor degeneration in RP.Comparing our results against the taxonomy of visual field loss patterns developed by Grover et al., 13 we find that our models recapitulate several patterns.Mutation-induced rod degeneration is predicted to produce a pattern in which degeneration initiates at the periphery and in the parafoveal/perifoveal regions and spreads into the midperiphery, similar to Pattern 1 and especially Pattern 1B.Likewise, if a patch of photoreceptor loss protrudes far enough into the peripheral retina, degeneration sweeps around the periphery of the retina, similar to the mid-late stage of Pattern 3. Patterns that involve a preferential loss of midperipheral photoreceptors, such as Pattern 2 and the initial stage of Pattern 3, are theoretically inaccessible to our models and hence to the oxygen toxicity mechanism.This is because the high density of photoreceptors and hence high rate of oxygen consumption in this region causes the propagation of hyperoxia-induced photoreceptor degeneration to halt (see Roberts et al. 59 for a detailed analysis).For the same reason, the central cone island is preserved, in agreement with   We note that the step-like appearance of the curves in some panels is the result of multiple regions falling below the p average threshold simultaneously.This effect is more pronounced for radiallysymmetric solutions (a, d, e) and is exaggerated in (e) due to the magnification of the y-axis.Treatment is applied at t crit ¼ 0.63 (yr) in (a), t crit ¼ 1.27 (yr) in (b), t crit ¼ 3.17 (yr) in (c), t crit ¼ 9.51 (yr) in (d), and t crit ¼ 15.84 (yr) in (e).See Table 2 for the remaining parameter values.
Grover et al., 13 though its eventual loss cannot be explained unless mutation-induced cone loss is active.Similarly, the hyperoxic photoreceptor degeneration triggered by mutationinduced cone loss stalls at it enters the midperiphery because this region is rod-dominated.Thus our theoretical models demonstrate, in a novel manner, the strengths and weaknesses of the oxygen toxicity hypothesis: while this mechanism is sufficient to generate some of the patterns of degeneration seen in vivo, other mechanisms are necessary to explain the remaining patterns.
Our models illustrate how treatment with antioxidants and/ or trophic factors could halt (or partially reverse) the spread of photoreceptor degeneration in those cases where a patch of photoreceptors has been lost and mutation-induced rod and cone degeneration are inactive.In these cases the treatments are effective because they raise the hyperoxic threshold and, hence, decrease the critical photoreceptor density at which degeneration ceases to propagate.As a result, degenerate patches may remain stable, without expanding, in a larger proportion of the retina.While partial recovery may be feasible where photoreceptor IS are still intact, regeneration is likely to be limited, as many photoreceptors will have degenerated past the point of recovery.In the case of mutation-induced rod degeneration, treatment can at best delay degeneration, as the steady attrition of rods will eventually bring the local photoreceptor density below the critical level required for the propagation of hyperoxic degeneration.Treatment may, however, halt hyperoxic degeneration associated with mutation-induced cone degeneration.This is because cone density is low outside of the fovea, the remainder of the retina being rod-dominated.These results suggest that, where possible, antioxidant/trophic factor treatment should be targeted to regions of the retina in which the photoreceptor density is low (e.g., using ocular gene therapy 4,43 ) as these regions are most susceptible to hyperoxic degeneration.
We plan to extend this work in a number of directions in future studies.We will incorporate heterogeneity in photoreceptor density in the azimuthal dimension and consider patterns of degeneration in other mammals such as rats, mice, rabbits, cats, pigs, and monkeys.This will allow us to explore further the strengths and weaknesses of the oxygen toxicity hypothesis and may help us to explain the distinct spatiotemporal patterns of photoreceptor degeneration associated with different species.We will also consider a 2D model, spanning the depth of the retina in 1D (as in Roberts et al. 74 ) and the region between the foveal center and the ora serrata in the other (neglecting the azimuthal dimension), to account for the oxygen distribution across the retinal depth and the effects of photoreceptor OS shrinkage upon the propagation of retinal degeneration.Lastly, we will use the present modeling framework to consider other disease mechanisms in RP, in particular the trophic factor, toxic substance, microglia, and mTOR hypotheses, both in isolation and in combination.
More experimental and clinical data are needed to validate and inform future modeling studies.The most useful data would be those derived from more detailed longitudinal studies in animals and humans, recording the spatio-temporal development of degenerate patches (using optical coherence tomography) and photoreceptor topography (using adaptive optics scanning laser ophthalmoscopy) across the retina.This would enable us to conduct a more detailed quantitative validation of our models.Measurements of the oxygen distribution across the retinal depth at a range of eccentricities and azimuthal angles at different stages of the disease would also be valuable in validating and extending our models, as would detailed quantitative data describing the effects of antioxidant and trophic factor treatment over time.
Our theoretical study has clarified the explanatory scope of the oxygen toxicity hypothesis as a disease mechanism underlying RP.Hyperoxia may well play a role in the spread of retinal degeneration, but other mechanisms are required to complete the picture.Our models further demonstrate the potential of antioxidant and trophic factor therapies in reducing hyperoxic disease progression, while highlighting the complementary insights that can be attained by combining a mathematical modeling approach with experimental and clinical studies.

FIGURE 1 .
FIGURE 1. Schematic showing the different regions of the retina (in the left eye, as viewed from the front).

FIGURE 2 .
FIGURE 2. Patterns of visual field loss as described by Grover et al. 13 Vision loss can be classified into three cases and six subcases.Large gray arrows indicate transitions between stages of vision loss and small red arrows indicate the direction of scotoma propagation.See text for details.

FIGURE 3 .
FIGURE 3. Diagrams to show the model geometry and 2D visualization.(a) Diagram of the human eye, viewed in the transverse plane.Model equations are solved on a 2D domain that extends from the center of the fovea (h ¼ 0 [rad]) to the ora serrata (h ¼ H [rad]). Diagram reproduced, with modifications, from http://www.nei.nih.gov/health/coloboma/coloboma,courtesy: National Eye Institute, National Institutes of Health (NEI/ NIH).(b) Equations are solved on a spherical cap (the healthy photoreceptor distribution is shown here), ignoring the optic nerve head and projected onto the xy-plane for visualization.This has the effect of preserving area near the center of the retina and compressing the area in the polar direction toward the periphery.

FIGURE 4 .
FIGURE 4. Sensitivity analysis demonstrating the degenerate retinal area resulting from an initial disc of photoreceptor loss.(a) Graph to show the effects of degenerate disc eccentricity, h c , and radius, governed by w, upon the resultant degenerate area at t ¼ 3 3 10 8 (s) (» 9.5 [yr]).The system has reached steady-state by this time for h c 0.7 3 H (rad) and is close to steady-state for h c > 0.7 3 H (rad).The inset shows a magnification of the bottom-left corner of the graph, demonstrating the degenerate area resulting from discs for which h c /H 0.2.An area is considered degenerate when photoreceptors have dropped beneath 10% of their average density, where p average ¼ 1.11 3 10 5 (photoreceptors mm -2 ).Degeneration is most severe for h c < 0.2 3 H (rad) and h c > 0.7 3 H (rad), while degeneration recovers for 0.3 3 H h c 0.4 3 H (rad). Degeneration is also more severe for larger values of w. (b-f) Panels show snapshots of the simulations at different stages in the disease progression as the initial radius and location of the degenerate patch vary.The plots show the projection of the solution from the retinal surface onto the (x,y) -plane (cm) (see Fig. 3b).Note that plots are magnified in (b, c) to clarify the dynamics.(b) h c ¼ 0 (rad), (c) h c ¼ 0.1 3 H (rad), (d) h c ¼ 0.3 3 H (rad), (e) h c ¼ 0.7 3 H (rad), and (f) h c ¼ 0.8 3 H (rad). (b-e) w ¼ 0.05 3 H (rad) and (f) w ¼ 0.0125 3 H (rad). See Table 2 for the remaining parameter values.

FIGURE 5 .
FIGURE 5. Simulation results illustrating the pattern of retinal degeneration associated with mutation-induced rod degeneration.Panels show the evolving rod (a) and cone (b) distributions over time.Note that we plot the natural logarithm of the (dimensionless) scaled cone density, pc ¼ p c =p average , where p average ¼ 1.11 3 10 5 (photoreceptors mm -2 ), in (b), with the lower bound of the color range fixed to the minimum value of logðp c ðh; /; tÞÞ at t ¼ 0 (yr) to improve visual clarity.Hyperoxia-induced photoreceptor loss initiates first at the ora serrata and later in the parafoveal/ perifoveal region, spreading into the midperiphery.Only the central cone island remains at steady-state.See Table2for parameter values.

FIGURE 6 .
FIGURE 6. Simulation results showing how mutation-induced rod degeneration can stimulate expansion of a pre-existing degenerate disc.Disc eccentricity h c ¼ 0.7 3 H (rad) and radius parameter w ¼ 0.05 3 H (rad) are the same as in Figure 4e, where the disc does not expand.Panels show the evolving rod (a) and cone (b) distributions over time.(See Fig. 5 legend for details of the plotting procedure.)SeeTable 2 for the remaining parameter values.

FIGURE 7 .
FIGURE 7. Simulation results illustrating the pattern of retinal degeneration associated with mutation-induced cone degeneration.Panels show the evolving rod (a) and cone (b) distributions over time.(See Fig. 5 legend for details of the plotting procedure.)Hyperoxiainduced photoreceptor loss initiates in the fovea, spreading outward and stalling as it enters the midperiphery.SeeTable 2 for parameter values.

FIGURE 8 .
FIGURE 8. Graphs to show the effect of treatment with antioxidants and/or trophic factors upon the degenerate retinal area.An area is considered degenerate when photoreceptors have dropped beneath 10% of their average density, where p average ¼ 1.11 3 10 5 (photoreceptors mm -2 ).(a) Degenerate annulus with (h 1 ,h 2 ) ¼ (0.08, 0.9) 3H (rad), (b) degenerate disc with (h c ,w) ¼ (0.1, 0.05) 3 H (rad), (c) degenerate disc with (h c ,w) ¼ (0.8, 0.0125) 3 H (rad), (d) mutation-induced rod degeneration, and (e) mutation-induced cone degeneration.Treatment is predicted to delay the spread of degeneration in (d), halt its spread in (a, c, e), and partially reverse degeneration in (b).We note that the step-like appearance of the curves in some panels is the result of multiple regions falling below the p average threshold simultaneously.This effect is more pronounced for radiallysymmetric solutions (a, d, e) and is exaggerated in (e) due to the magnification of the y-axis.Treatment is applied at t crit ¼ 0.63 (yr) in (a), t crit ¼ 1.27 (yr) in (b), t crit ¼ 3.17 (yr) in (c), t crit ¼ 9.51 (yr) in (d), and t crit ¼ 15.84 (yr) in (e).See Table2for the remaining parameter values.

TABLE 1 .
Summary of the Physical Variables Used in the Mathematical Models (seeEquations 1-14)

TABLE 2 .
Parameter Values for the Three Mathematical Models (see Equations 1-7; Equations 8-12; and Equations 13,14), Given to a Maximum of Three Significant Figures