Pattern formation of the cone mosaic in the zebrafish retina: a cell rearrangement model

Neurodegeneration, Neuroprotection and Regeneration in the Zebrafish Retina

Neurodegenerative retinal diseases, such as glaucoma and diabetic retinopathy, involve a gradual loss of neurons in the retina as the disease progresses. Central nervous system neurons are not able to regenerate in mammals, therefore, an often sought after course of treatment for neuronal loss follows a neuroprotective or regenerative strategy. Neuroprotection is the process of preserving the structure and function of the neurons that have survived a harmful insult; while regenerative approaches aim to replace or rewire the neurons and synaptic connections that were lost, or induce regrowth of damaged axons or dendrites. In order to test the neuroprotective effectiveness or the regenerative capacity of a particular agent, a robust experimental model of retinal neuronal damage is essential. Zebrafish are being used more often in this type of study because their eye structure and development is well-conserved between zebrafish and mammals. Zebrafish are robust genetic tools and are relatively inexpensive to maintain. The large array of functional and behavioral tests available in zebrafish makes them an attractive model for neuroprotection studies. Some common insults used to model retinal disease and study neuroprotection in zebrafish include intense light, chemical toxicity and mechanical damage. This review covers the existing retinal neuroprotection and regeneration literature in the zebrafish and highlights their potential for future studies.

Dynamical pattern selection of growing cellular mosaic in fish retina

A Markovian lattice model for photoreceptor cells is introduced to describe the growth of mosaic patterns on fish retina. The radial stripe pattern observed in wild-type zebrafish is shown to be selected naturally during retina growth, against the geometrically equivalent circular stripe pattern. The mechanism of such dynamical pattern selection is clarified on the basis of both numerical simulations and theoretical analyses, which find that the successive emergence of local defects plays a critical role in the realization of the wild-type pattern.

Investigation of Adaptive Optics Imaging Biomarkers for Detecting Pathological Changes of the Cone Mosaic in Patients with Type 1 Diabetes Mellitus

PURPOSE

To investigate a set of adaptive optics (AO) imaging biomarkers for the assessment of changes of the cone mosaic spatial arrangement in patients with type 1 diabetes mellitus (DM1).

METHODS

16 patients with ≥20/20 visual acuity and a diagnosis of DM1 in the past 8 years to 37 years and 20 age-matched healthy volunteers were recruited in this study. Cone density, cone spacing and Voronoi diagrams were calculated on 160x160 μm images of the cone mosaic acquired with an AO flood illumination retinal camera at 1.5 degrees eccentricity from the fovea along all retinal meridians. From the cone spacing measures and Voronoi diagrams, the linear dispersion index (LDi) and the heterogeneity packing index (HPi) were computed respectively. Logistic regression analysis was conducted to discriminate DM1 patients without diabetic retinopathy from controls using the cone metrics as predictors.

RESULTS

Of the 16 DM1 patients, eight had no signs of diabetic retinopathy (noDR) and eight had mild nonproliferative diabetic retinopathy (NPDR) on fundoscopy. On average, cone density, LDi and HPi values were significantly different (P<0.05) between noDR or NPDR eyes and controls, with these differences increasing with duration of diabetes. However, each cone metric alone was not sufficiently sensitive to discriminate entirely between membership of noDR cases and controls. The complementary use of all the three cone metrics in the logistic regression model gained 100% accuracy to identify noDR cases with respect to controls.

CONCLUSION

The present set of AO imaging biomarkers identified reliably abnormalities in the spatial arrangement of the parafoveal cones in DM1 patients, even when no signs of diabetic retinopathy were seen on fundoscopy.

A Cellular Potts Model of single cell migration in presence of durotaxis

Cell migration is a fundamental biological phenomenon during which cells sense their surroundings and respond to different types of signals. In presence of durotaxis, cells preferentially crawl from soft to stiff substrates by reorganizing their cytoskeleton from an isotropic to an anisotropic distribution of actin filaments. In the present paper, we propose a Cellular Potts Model to simulate single cell migration over flat substrates with variable stiffness. We have tested five configurations: (i) a substrate including a soft and a stiff region, (ii) a soft substrate including two parallel stiff stripes, (iii) a substrate made of successive stripes with increasing stiffness to create a gradient and (iv) a stiff substrate with four embedded soft squares. For each simulation, we have evaluated the morphology of the cell, the distance covered, the spreading area and the migration speed. We have then compared the numerical results to specific experimental observations showing a consistent agreement.

Functional architecture of the retina: development and disease

Structure and function are highly correlated in the vertebrate retina, a sensory tissue that is organized into cell layers with microcircuits working in parallel and together to encode visual information. All vertebrate retinas share a fundamental plan, comprising five major neuronal cell classes with cell body distributions and connectivity arranged in stereotypic patterns. Conserved features in retinal design have enabled detailed analysis and comparisons of structure, connectivity and function across species. Each species, however, can adopt structural and/or functional retinal specializations, implementing variations to the basic design in order to satisfy unique requirements in visual function. Recent advances in molecular tools, imaging and electrophysiological approaches have greatly facilitated identification of the cellular and molecular mechanisms that establish the fundamental organization of the retina and the specializations of its microcircuits during development. Here, we review advances in our understanding of how these mechanisms act to shape structure and function at the single cell level, to coordinate the assembly of cell populations, and to define their specific circuitry. We also highlight how structure is rearranged and function is disrupted in disease, and discuss current approaches to re-establish the intricate functional architecture of the retina.

Regeneration of cone photoreceptors when cell ablation is primarily restricted to a particular cone subtype

We sought to characterize the regenerated cells, if any, when photoreceptor ablation was mostly limited to a particular cone subtype. This allowed us to uniquely assess whether the remaining cells influence specification of regenerating photoreceptors. The ability to replace lost photoreceptors via stem cell therapy holds promise for treating many retinal degenerative diseases. Zebrafish are potent for modelling this because they have robust regenerative capacity emanating from endogenous stem cells, and abundant cone photoreceptors including multiple spectral subtypes similar to human fovea. We ablated the homolog of the human S-cones, the ultraviolet-sensitive (UV) cones, and tested the hypothesis that the photoreceptors regenerating in their place take on identities matching those expected from normal cone mosaic development. We created transgenic fish wherein UV cones can be ablated by addition of a prodrug. Thus photoreceptors developed normally and only the UV cones expressed nitroreductase; the latter converts the prodrug metronidazole to a cell-autonomous neurotoxin. A significant increase in proliferation of progenitor cell populations (p<0.01) was observed when cell ablation was primarily limited to UV cones. In control fish, we found that BrdU primarily incorporated into rod photoreceptors, as expected. However the majority of regenerating photoreceptors became cones when retinal cell ablation was predominantly restricted to UV cones: a 2-fold increase in the relative abundance of cones (p = 0.008) was mirrored by a 35% decrease in rods. By primarily ablating only a single photoreceptor type, we show that the subsequent regeneration is biased towards restoring the cognate photoreceptor type. We discuss the hypothesis that, after cone death, the microenvironment formed by the remaining retinal cells may be influential in determining the identity of regenerating photoreceptors, though other interpretations are plausible. Our novel animal model provides control of ablation that will assist in identifying mechanisms required to replace cone photoreceptors clinically to restore daytime vision.

Crb apical polarity proteins maintain zebrafish retinal cone mosaics via intercellular binding of their extracellular domains

Cone photoreceptors are assembled by unknown mechanisms into geometrically regular mosaics in many vertebrate species. The formation and maintenance of photoreceptor mosaics are speculated to require differential cell-cell adhesion. However, the molecular basis for this theory has yet to be identified. The retina and many other tissues express Crumbs (Crb) polarity proteins. The functions of the extracellular domains of Crb proteins remain to be understood. Here we report cell-type-specific expression of the crb2a and crb2b genes at the cell membranes of photoreceptor inner segments and Müller cell apical processes in the zebrafish retina. We demonstrate that the extracellular domains of Crb2a and Crb2b mediate a cell-cell adhesion function, which plays an essential role in maintaining the integrity of photoreceptor layer and cone mosaics. Because Crb proteins are expressed in many types of epithelia, the Crb-based cell-cell adhesion may underlie cellular patterning in other epithelium-derived tissues as well.

Multi-scale modeling of tissues using CompuCell3D

The study of how cells interact to produce tissue development, homeostasis, or diseases was, until recently, almost purely experimental. Now, multi-cell computer simulation methods, ranging from relatively simple cellular automata to complex immersed-boundary and finite-element mechanistic models, allow in silico study of multi-cell phenomena at the tissue scale based on biologically observed cell behaviors and interactions such as movement, adhesion, growth, death, mitosis, secretion of chemicals, chemotaxis, etc. This tutorial introduces the lattice-based Glazier-Graner-Hogeweg (GGH) Monte Carlo multi-cell modeling and the open-source GGH-based CompuCell3D simulation environment that allows rapid and intuitive modeling and simulation of cellular and multi-cellular behaviors in the context of tissue formation and subsequent dynamics. We also present a walkthrough of four biological models and their associated simulations that demonstrate the capabilities of the GGH and CompuCell3D.

Ontogeny of cone photoreceptor mosaics in zebrafish

Cone photoreceptors in fish are typically arranged into a precise, reiterated pattern known as a "cone mosaic." Cone mosaic patterns can vary in different fish species and in response to changes in habitat, yet their function and the mechanisms of their development remain speculative. Zebrafish (Danio rerio) have four cone subtypes arranged into precise rows in the adult retina. Here we describe larval zebrafish cone patterns and investigate a previously unrecognized transition between larval and adult cone mosaic patterns. Cone positions were determined in transgenic zebrafish expressing green fluorescent protein (GFP) in their UV-sensitive cones, by the use of multiplex in situ hybridization labelling of various cone opsins. We developed a "mosaic metric" statistical tool to measure local cone order. We found that ratios of the various cone subtypes in larval and adult zebrafish were statistically different. The cone photoreceptors in larvae form a regular heterotypic mosaic array; i.e., the position of any one cone spectral subtype relative to the other cone subtypes is statistically different from random. However, the cone spectral subtypes in larval zebrafish are not arranged in continuous rows as in the adult. We used cell birth dating to show that the larval cone mosaic pattern remains as a distinct region within the adult retina and does not reorganize into the adult row pattern. In addition, the abundance of cone subtypes relative to other subtypes is different in this larval remnant compared with that of larvae or canonical adult zebrafish retina. These observations provide baseline data for understanding the development of cone mosaics via comparative analysis of larval and adult cone development in a model species.

Avian cone photoreceptors tile the retina as five independent, self-organizing mosaics

The avian retina possesses one of the most sophisticated cone photoreceptor systems among vertebrates. Birds have five types of cones including four single cones, which support tetrachromatic color vision and a double cone, which is thought to mediate achromatic motion perception. Despite this richness, very little is known about the spatial organization of avian cones and its adaptive significance. Here we show that the five cone types of the chicken independently tile the retina as highly ordered mosaics with a characteristic spacing between cones of the same type. Measures of topological order indicate that double cones are more highly ordered than single cones, possibly reflecting their posited role in motion detection. Although cones show spacing interactions that are cell type-specific, all cone types use the same density-dependent yardstick to measure intercone distance. We propose a simple developmental model that can account for these observations. We also show that a single parameter, the global regularity index, defines the regularity of all five cone mosaics. Lastly, we demonstrate similar cone distributions in three additional avian species, suggesting that these patterning principles are universal among birds. Since regular photoreceptor spacing is critical for uniform sampling of visual space, the cone mosaics of the avian retina represent an elegant example of the emergence of adaptive global patterning secondary to simple local interactions between individual photoreceptors. Our results indicate that the evolutionary pressures that gave rise to the avian retina's various adaptations for enhanced color discrimination also acted to fine-tune its spatial sampling of color and luminance.

Multicell simulations of development and disease using the CompuCell3D simulation environment

Mathematical modeling and computer simulation have become crucial to biological fields from genomics to ecology. However, multicell, tissue-level simulations of development and disease have lagged behind other areas because they are mathematically more complex and lack easy-to-use software tools that allow building and running in silico experiments without requiring in-depth knowledge of programming. This tutorial introduces Glazier-Graner-Hogeweg (GGH) multicell simulations and CompuCell3D, a simulation framework that allows users to build, test, and run GGH simulations.

Development and adult morphology of the eye lens in the zebrafish

The zebrafish has become an important vertebrate model organism to study the development of the visual system. Mutagenesis projects have resulted in the identification of hundreds of eye mutants. Analysis of the phenotypes of these mutants relies on in depth knowledge of the embryogenesis in wild-type animals. While the morphological events leading to the formation of the retina and its connections to the central nervous system have been described in great detail, the characterization of the development of the eye lens is still incomplete. In the present study, we provide a morphological description of embryonic and larval lens development as well as adult lens morphology in the zebrafish. Our analyses show that, in contrast to other vertebrate species, the zebrafish lens delaminates from the surface ectoderm as a solid cluster of cells. Detachment of the prospective lens from the surface ectoderm is facilitated by apoptosis. Primary fibre cell elongation occurs in a circular fashion resulting in an embryonic lens nucleus with concentric shells of fibres. After formation of a monolayer of lens epithelial cells, differentiation and elongation of secondary lens fibres result in a final lens morphology similar to that of other vertebrate species. As in other vertebrates, secondary fibre cell differentiation includes the programmed degradation of nuclei, the interconnection of adjacent fibres via protrusions at the fibre cells' edges and the establishment of gap junctions between lens fibre cells. The very close spacing of the nuclei of the differentiating secondary fibres in a narrow zone close to the equatorial epithelium, however, suggests that secondary fibre cell differentiation deviates from that described for mammalian or avian lenses. In summary, while there are similarities in the development and final morphology of the zebrafish lens with mammalian and avian lenses, there are also significant differences, suggesting caution when extrapolating findings on the zebrafish to, for example, human lens development or function.

Have we achieved a unified model of photoreceptor cell fate specification in vertebrates?

How does a retinal progenitor choose to differentiate as a rod or a cone and, if it becomes a cone, which one of their different subtypes? The mechanisms of photoreceptor cell fate specification and differentiation have been extensively investigated in a variety of animal model systems, including human and non-human primates, rodents (mice and rats), chickens, frogs (Xenopus) and fish. It appears timely to discuss whether it is possible to synthesize the resulting information into a unified model applicable to all vertebrates. In this review we focus on several widely used experimental animal model systems to highlight differences in photoreceptor properties among species, the diversity of developmental strategies and solutions that vertebrates use to create retinas with photoreceptors that are adapted to the visual needs of their species, and the limitations of the methods currently available for the investigation of photoreceptor cell fate specification. Based on these considerations, we conclude that we are not yet ready to construct a unified model of photoreceptor cell fate specification in the developing vertebrate retina.

Remodeling of the cone photoreceptor mosaic during metamorphosis of flounder (Pseudopleuronectes americanus)

The retinal cone mosaic of the winter flounder, Pseudopleuronectes americanus, is extensively remodeled during metamorphosis when its visual system shifts from monochromatic to trichromatic. Here we describe the reorganization and re-specification of existing cone subtypes in which larval cones alter their spatial arrangement, morphology, and opsin expression to determine whether mechanisms controlling cell birth, mosaic position, and opsin selection are coordinated or independent. We labeled dividing cells with tritiated ((3)H) thymidine prior to mosaic remodeling to determine whether existing cone photoreceptors change phenotype. We also used in situ hybridization to identify mosaic type and opsin expression in transitional retinas to understand the sequence of transformation. Our data indicate that in the winter flounder retina the choice of new opsin species and the cellular rearrangement of the mosaic proceed independently. The production of the precise cone mosaic arrangement is not due to a stereotyped series of sequential cellular inductions, but rather might be the product of a set of distinct, flexible processes that rely on plasticity in cell phenotype.

Possibility of tissue separation caused by cell adhesion

Adhesive molecules are suggested to play an important role when a single tissue is separated into two in developmental processes, illustrated by tissue-specific cadherins in the neural tube formation of amphibians. In this paper, we study the possibility for tissue separation to be carried out only by differential cell adhesion and random cell movement without any other morphogenetic mechanisms. We consider a two-dimensional regular triangular lattice filled with cells of three types (black, white, and gray). In the initial state, a cluster of black cells and a cluster of white cells are in contact and are surrounded by gray cells. Nearest-neighbor cells exchange their location at random, but the movement occurs faster if it increases the total adhesion. We considered separation to be successful if, in the final state, black cells and white cells kept their clusters but two clusters lost their direct contact with each other as gray cells are inserted between them. The maximum total adhesion (MTA) rule conjectures that the spatial pattern achieving maximum total adhesion might be that obtained in the final state. In the computer simulation, the runs for successful separation satisfied the condition predicted by the MTA rule. However, the condition for successful separation was more restricted than that predicted by the MTA rule. For some combinations of adhesions, it took an extremely long time to accomplish tissue separation. Finally, we discuss the role of homophilic adhesion molecules (such as cadherins) in the tissue separation processes, and show that the new expression of homophilic adhesion molecules cannot perform tissue separation without the change in other morphogenetic processes.

Difference in the retinal cone mosaic pattern between zebrafish and medaka: cell-rearrangement model

In fish retina, four kinds of photoreceptor cells (or cones) are two-dimensionally arranged in a very regular manner, forming cone mosaics. Mosaic pattern differs between species--two typical patterns are "row mosaic" and "square mosaic", exemplified by the cone mosaics in zebrafish and in medaka, respectively. In this paper, we study a cell-rearrangement model. Cells with pre-fixed fate exchange their locations between nearest neighbors and form regular mosaic patterns spontaneously, if the adhesive force between nearest neighbors and between next-nearest neighbors depend on their cell types in an appropriate manner. The same model can produce both row and square mosaic patterns. However, if the cell-cell interaction is restricted to nearest neighbors only, the square mosaic (medaka pattern) cannot be generated, showing the importance of interaction between next-nearest neighbors. In determining whether row mosaic (zebrafish pattern) or square mosaic (medaka pattern) is to be formed, two shape factors are very important, which control the way adhesions in different geometric relations are combined. We also developed theoretical analysis of the parameter ranges for the row mosaic and the square mosaic to have higher total adhesion than alternative spatial patterns.

Influence of cell fate mechanisms upon retinal mosaic formation: a modelling study

Many types of retinal neurone are arranged in a spatially regular manner so that the visual scene is uniformly sampled. Several mechanisms are thought to be involved in the development of regular cellular positioning. One early-acting mechanism is the lateral inhibition of neighbouring cells from acquiring the same fate, mediated by Delta-Notch signalling. We have used computer modelling to test whether lateral inhibition might transform an initial population of undifferentiated cells into more regular populations of two types of differentiated cells. Initial undifferentiated cells were positioned randomly, subject only to a minimal distance constraint. Each undifferentiated cell then acquired either primary or secondary fate using one of several lateral inhibition mechanisms. Mosaic regularity was assessed using the regularity index and the packing factor. We found that for irregular undifferentiated mosaics, the arrangement of resulting primary (but not secondary) fate cells was more regular than in the initial undifferentiated population. However, for regular undifferentiated mosaics, no further increases in the regularity of the primary fate mosaics were observed. We have used this model to test the specific hypothesis that on- and off-centre retinal ganglion cells emerge from an initial, irregular undifferentiated population of ganglion cells. Lateral inhibition can subdivide an initially irregular population into two types of cell that are mildly regular. However, lateral inhibition alone is insufficient to produce mosaics of the same regularity as observed experimentally. Likewise, and in contrast to earlier reports, cell death alone is insufficient to match the regularity of experimental mosaics. We conclude that lateral inhibition can transform irregular distributions into regular mosaics, upon which subsequent processes (such as lateral cell movement or cell death) can further refine mosaic regularity.