The sizes of different neuronal populations within the CNS are precisely controlled, but whether neuronal number is coordinated between cell types is unknown. We examined the covariance structure of 12 different retinal cell types across 30 genetically distinct lines of mice, finding minimal covariation when comparing synaptically connected or developmentally related cell types. Variation mapped to one or more genomic loci for each cell type, but rarely were these shared, indicating minimal genetic coregulation of final number. Multiple genes, therefore, participate in the specification of the size of every population of retinal neuron, yet genetic variants work largely independent of one another during development to modulate those numbers, yielding substantial variability in the convergence ratios between pre- and postsynaptic populations. Density-dependent cellular interactions in the outer plexiform layer overcome this variability to ensure the formation of neuronal circuits that maintain constant retinal coverage and complete afferent sampling.

The present study interrogated a quantitative trait locus (QTL) on Chr 4 associated with the population sizes of two types of bipolar cell in the mouse retina. This locus was identified by quantifying the number of rod bipolar cells and Type 2 cone bipolar cells across a panel of recombinant inbred (RI) strains of mice derived from two inbred laboratory strains, C57BL/6J (B6/J) and A/J, and mapping a proportion of that variation in cell number, for each cell type, to this shared locus. There, we identified the candidate gene X Kell blood group precursor related family member 8 homolog (Xkr8). While Xkr8 has no documented role in the retina, we localize robust expression in the mature retina via in situ hybridization, confirm its developmental presence via immunolabeling, and show that it is differentially regulated during the postnatal period between the B6/J and A/J strains using qPCR. Microarray analysis, derived from whole eye mRNA from the entire RI strain set, demonstrates significant negative correlation of Xkr8 expression with the number of each of these two types of bipolar cells, and the variation in Xkr8 expression across the strains maps a cis-eQTL, implicating a regulatory variant discriminating the parental genomes. Xkr8 plasmid electroporation during development yielded a reduction in the number of bipolar cells in the retina, while sequence analysis of Xkr8 in the two parental strain genomes identified a structural variant in the 3' UTR that may disrupt mRNA stability, and two SNPs in the promoter that create transcription factor binding sites. We propose that Xkr8, via its participation in mediating cell death, plays a role in the specification of bipolar cell number in the retina.

The population sizes of different retinal cell types vary between different strains of mice, and that variation can be mapped to genomic loci in order to identify its polygenic origin. In some cases, controlling genes act independently, whereas in other instances, they exhibit epistasis. Here, we identify an epistatic interaction revealed through the mapping of quantitative trait loci from a panel of recombinant inbred strains of mice. The population of retinal horizontal cells exhibits a twofold variation in number, mapping to quantitative trait loci on chromosomes 3 and 13, where these loci are shown to interact epistatically. We identify a prospective genetic interaction underlying this, mediated by the bHLH transcription factor Neurog2, at the chromosome 3 locus, functioning to repress the LIM homeodomain transcription factor Isl1, at the chromosome 13 locus. Using single and double conditional knockout mice, we confirm the countervailing actions of each gene, and validate in vitro a crucial role for two single nucleotide polymorphisms in the 5'UTR of Isl1, one of which yields a novel E-box, mediating the repressive action of Neurog2.

Programmed cell death contributes to the histogenesis of the nervous system, and is believed to be modulated through the sustaining effects of afferents and targets during the period of synaptogenesis. Cone bipolar cells undergo programmed cell death during development, and we confirm that the numbers of three different types are increased when the pro-apoptotic Bax gene is knocked out. When their cone afferents are selectively eliminated, or when the population of retinal ganglion cells is increased, however, cone bipolar cell number remains unchanged. Programmed cell death of the cone bipolar cell populations, therefore, may be modulated cell-intrinsically rather than via interactions with these synaptic partners.

The retina consists of ordered arrays of individual types of neurons for processing vision. Here, we show that such order is necessary for intrinsically photosensitive retinal ganglion cells (ipRGCs) to function as irradiance detectors. We found that during development, ipRGCs undergo proximity-dependent Bax-mediated apoptosis. Bax mutant mice exhibit disrupted ipRGC spacing and dendritic stratification with an increase in abnormally localized synapses. ipRGCs are the sole conduit for light input to circadian photoentrainment, and either their melanopsin-based photosensitivity or ability to relay rod/cone input is sufficient for circadian photoentrainment. Remarkably, the disrupted ipRGC spacing does not affect melanopsin-based circadian photoentrainment but severely impairs rod/cone-driven photoentrainment. We demonstrate reduced rod/cone-driven cFos activation and electrophysiological responses in ipRGCs, suggesting that impaired synaptic input to ipRGCs underlies the photoentrainment deficits. Thus, for irradiance detection, developmental apoptosis is necessary for the spacing and connectivity of ipRGCs that underlie their functioning within a neural network.

Sequence variants modulating gene function or expression affect various heritable traits, including the number of neurons within a population. The present study employed a forward-genetic approach to identify candidate causal genes and their sequence variants controlling the number of one type of retinal neuron, the AII amacrine cell. Data from twenty-six recombinant inbred (RI) strains of mice derived from the parental C57BL/6J (B6/J) and A/J laboratory strains were used to identify genomic loci regulating cell number. Large variation in cell number is present across the RI strains, from a low of ∼57,000 cells to a high of ∼87,000 cells. Quantitative trait locus (QTL) analysis revealed three prospective controlling genomic loci, on Chromosomes (Chrs) 9, 11, and 19, each contributing additive effects that together approach the range of variation observed. Composite interval mapping validated two of these loci, and chromosome substitution strains, in which the A/J genome for Chr 9 or 19 was introgressed on a B6/J genetic background, showed increased numbers of AII amacrine cells as predicted by those two QTL effects. Analysis of the respective genomic loci identified candidate controlling genes defined by their retinal expression, their established biological functions, and by the presence of sequence variants expected to modulate gene function or expression. Two candidate genes, Dtx4 on Chr 19, being a regulator of Notch signaling, and Dixdc1 on Chr 9, a modulator of the WNT-β-catenin signaling pathway, were explored in further detail. Postnatal overexpression of Dtx4 was found to reduce the frequency of amacrine cells, while Dixdc1 knockout retinas contained an excess of AII amacrine cells. Sequence variants in each gene were identified, being the likely sources of variation in gene expression, ultimately contributing to the final number of AII amacrine cells.

The present study has taken advantage of publicly available cell type specific mRNA expression databases in order to identify potential genes participating in the development of retinal AII amacrine cells. We profile two such genes, Delta/Notch-like EGF repeat containing (Dner) and nuclear factor I/A (Nfia), that are each heavily expressed in AII amacrine cells in the mature mouse retina, and which conjointly identify this retinal cell population in its entirety when using antibodies to DNER and NFIA. DNER is present on the plasma membrane, while NFIA is confined to the nucleus, consistent with known functions of each of these two proteins. DNER also identifies some other subsets of retinal ganglion and amacrine cell types, along with horizontal cells, while NFIA identifies a subset of bipolar cells as well as Muller glia and astrocytes. During early postnatal development, NFIA labels astrocytes on the day of birth, AII amacrine cells at postnatal (P) day 5, and Muller glia by P10, when horizontal cells also transiently exhibit NFIA immunofluorescence. DNER, by contrast, is present in ganglion and amacrine cells on P1, also labeling the horizontal cells by P10. Developing AII amacrine cells exhibit accumulating DNER labeling at the dendritic stalk, labeling that becomes progressively conspicuous by P10, as it is in maturity. This developmental time course is consistent with a prospective role for each gene in the differentiation of AII amacrine cells.

The two populations of cholinergic amacrine cells in the inner nuclear layer (INL) and the ganglion cell layer (GCL) differ in their spatial organization in the mouse retina, but the basis for this difference is not understood. The present investigation examined this issue in six strains of mice that differ in their number of cholinergic cells, addressing how the regularity, packing, and spacing of these cells varies as a function of strain, layer, and density. The number of cholinergic cells was lower in the GCL than in the INL in all six strains. The nearest neighbor and Voronoi domain regularity indexes as well as the packing factor were each consistently lower for the GCL. While these regularity indexes and the packing factor were largely stable across variation in density, the effective radius was inversely related to density for both the GCL and INL, being smaller and more variable in the GCL. Consequently, despite the lower densities in the GCL, neighboring cells were more likely to be positioned closer to one another than in the higher-density INL, thereby reducing regularity and packing. This difference in the spatial organization of cholinergic cells may be due to the cells in the GCL having been passively displaced by fascicles of optic axons and an expanding retinal vasculature during development. In support of this interpretation, we show such displacement of cholinergic somata relative to their dendritic stalks and a decline in packing efficiency and regularity during postnatal development that is more severe for the GCL.

The precise specification of cellular fate is thought to ensure the production of the correct number of neurons within a population. Programmed cell death may be an additional mechanism controlling cell number, believed to refine the proper ratio of pre- to post-synaptic neurons for a given species. Here, we consider the size of three different neuronal populations in the rod pathway of the mouse retina: rod photoreceptors, rod bipolar cells, and AII amacrine cells. Across a collection of 28 different strains of mice, large variation in the numbers of all three cell types is present. The variation in their numbers is not correlated, so that the ratio of rods to rod bipolar cells, as well as rod bipolar cells to AII amacrine cells, varies as well. Establishing connectivity between such variable pre- and post-synaptic populations relies upon plasticity that modulates process outgrowth and morphological differentiation, which we explore experimentally for both rod bipolar and AII amacrine cells in a mouse retina with elevated numbers of each cell type. While both rod bipolar dendritic and axonal arbors, along with AII lobular arbors, modulate their areal size in relation to local homotypic cell densities, the dendritic appendages of the AII amacrine cells do not. Rather, these processes exhibit a different form of plasticity, regulating the branching density of their overlapping arbors. Each form of plasticity should ensure uniformity in retinal coverage in the presence of the independent specification of afferent and target cell number.

Multiple factors regulate the differentiation of neuronal morphology during development, including interactions with afferents, targets, and homotypic neighbors, as well as cell-intrinsic transcriptional regulation. Retinal bipolar cells provide an exemplary model system for studying the control of these processes, as there are 15 transcriptionally and morphologically distinct types, each extending their dendritic and axonal arbors in respective strata within the synaptic layers of the retina. Here we have examined the role of the transcription factor Sox5 in the control of the morphological differentiation of one type of cone bipolar cell (CBC), the Type 7 cell. We confirm selective expression of SOX5 in this single bipolar cell type, emerging at the close of the first post-natal week, prior to morphological differentiation. Conditional knockout mice were generated by crossing a bipolar cell-specific cre-expressing line with mice carrying floxed Sox5 alleles, as well as the Gustducin-gfp reporter which labels Type 7 CBCs. Loss of SOX5 was confirmed in the bipolar cell stratum, in GFP+ Type 7 cells. Such SOX5-deficient Type 7 cells differentiate axonal and dendritic arbors that are each reduced in areal extent. The axonal arbors exhibit sprouting in the inner plexiform layer (IPL), thereby extending their overall radial extent, while the dendritic arbors connect with fewer cone pedicles in the outer plexiform layer, showing an increase in the average number of dendritic contacts at each pedicle. SOX5-deficient Type 7 CBCs should therefore exhibit smaller receptive fields derived from fewer if now hyper-innervated pedicles, transmitting their signals across a broader depth through the IPL.

Sox2 is a transcriptional regulator that is highly expressed in retinal astrocytes, yet its function in these cells has not previously been examined. To understand its role, we conditionally deleted Sox2 from the population of astrocytes and examined the consequences on retinal development. We found that Sox2 deletion does not alter the migration of astrocytes, but it impairs their maturation, evidenced by the delayed upregulation of glial fibrillary acidic protein (GFAP) across the retina. The centro-peripheral gradient of angiogenesis is also delayed in Sox2-CKO retinas. In the mature retina, we observed lasting abnormalities in the astrocytic population evidenced by the sporadic loss of GFAP immunoreactivity in the peripheral retina as well as by the aberrant extension of processes into the inner retina. Blood vessels in the adult retina are also under-developed and show a decrease in the frequency of branch points and in total vessel length. The developmental relationship between maturing astrocytes and angiogenesis suggests a causal relationship between the astrocytic loss of Sox2 and the vascular architecture in maturity. We suggest that the delay in astrocytic maturation and vascular invasion may render the retina hypoxic, thereby causing the abnormalities we observe in adulthood. These studies uncover a novel role for Sox2 in the development of retinal astrocytes and indicate that its removal can lead to lasting changes to retinal homeostasis.

Many types of retinal neuron modulate the distribution of their processes to ensure a uniform coverage of the retinal surface. Dendritic field area, for instance, is inversely related to the variation in cellular density for many cell types, observed either across retinal eccentricity or between different strains of mice that differ in cell number. Dopaminergic amacrine (DA) cells, by contrast, have dendritic arbors that bear no spatial relationship to the presence of their immediate homotypic neighbors, yet it remains to be determined whether their coverage upon the retina, as a population, is conserved across variation in their total number. The present study assessed the overall density of the dopaminergic plexus in the inner plexiform layer in the presence of large variation in the total number of DA cells, as well as their retinal dopamine content, to determine whether either of these features is conserved. We first compared these traits between two strains of mice (C57BL/6J and A/J) that exhibit a two-fold difference in DA cell number. We subsequently examined these same traits in littermate mice for which the pro-apoptotic Bax gene was either intact or knocked out, yielding a five-fold difference in DA cell number. In both comparisons, we found greater plexus density and DA content in the strain or condition with the greater number of DA cells. The population of DA cells, therefore, does not appear to self-regulate its process coverage to achieve a constant density as the DA mosaic is established during development, nor its functional dopamine content in maturity.

Genetic variants modulate the numbers of various retinal cell types in mice. For instance, there is minimal variation in the number of rod bipolar cells (RBCs) in two inbred strains of mice (A/J and C57BL/6J), yet their F1 offspring contain significantly more RBCs than either parental strain. To investigate the genetic source of this variation, we mapped the variation in the number of RBCs across 24 genetically distinct recombinant inbred (RI) strains (the AXB/BXA strain-set), seeking to identify quantitative trait loci (QTL). We then sought to identify candidate genes and potential casual variants at those genomic loci. Variation in RBC number mapped to three genomic loci, each modulating cell number in excess of one-third of the range observed across the RI strains. At each of these loci, we identified candidate genes containing variants that might alter gene function or expression. The latter genes were also analyzed using a transcriptome database, revealing a subset for which expression correlated with variation in RBC number. Using an electroporation strategy, we demonstrate that early postnatal expression of one of them, Ggct (gamma-glutamyl cyclotransferase), modulates bipolar cell number. We identify candidate regulatory variants for this gene, finding a large structural variant (SV) in the putative promoter that reduces expression using a luciferase assay. This SV reducing Ggct expression in vitro is likely the causal variant within the gene associated with the variation in Ggct expression in vivo, implicating it as a quantitative trait variant (QTV) participating in the control of RBC number.

Retinal neurons are often arranged as non-random distributions called "mosaics," as their somata minimize proximity to neighboring cells of the same type. The horizontal cells serve as an example of such a mosaic, but little is known about the developmental mechanisms that underlie their patterning. To identify genes involved in this process, we have used three different spatial statistics to assess the patterning of the horizontal cell mosaic across a panel of genetically distinct recombinant inbred strains. To avoid the confounding effect of cell density, which varies twofold across these different strains, we computed the "real/random regularity ratio," expressing the regularity of a mosaic relative to a randomly distributed simulation of similarly sized cells. To test whether this latter statistic better reflects the variation in biological processes that contribute to horizontal cell spacing, we subsequently compared the genomic linkage for each of these two traits, the regularity index, and the real/random regularity ratio, each computed from the distribution of nearest neighbor (NN) distances and from the Voronoi domain (VD) areas. Finally, we compared each of these analyses with another index of patterning, the packing factor. Variation in the regularity indexes, as well as their real/random regularity ratios, and the packing factor, mapped quantitative trait loci to the distal ends of Chromosomes 1 and 14. For the NN and VD analyses, we found that the degree of linkage was greater when using the real/random regularity ratio rather than the respective regularity index. Using informatic resources, we narrowed the list of prospective genes positioned at these two intervals to a small collection of six genes that warrant further investigation to determine their potential role in shaping the patterning of the horizontal cell mosaic.

Individual types of retinal neurons are distributed to minimize proximity to neighboring cells. Many of these same cell types extend dendrites to provide coverage of the retinal surface. These two cardinal features of retinal mosaics are disrupted, for certain cell types, in mice deficient for the Down syndrome cell adhesion molecule, Dscam, exhibiting an aberrant clustering of somata and fasciculation of dendrites. The Dscam mutant mouse retina also exhibits excess numbers of these same cell types. The present study compared these two features in Dscam mutant retinas with the Bax knockout retina, in which excess numbers of two of these cell types, the melanopsin-positive retinal ganglion cells (MRGCs) and the dopaminergic amacrine cells (DACs), are also present. Whole retinas were immunolabeled for both populations, and every labeled soma was plotted. For the MRGCs, we found a gene dosage effect for Dscam, with the Dscam+/- retinas showing smaller increases in cell number, clustering, and fasciculation. Curiously, Bax-/- retinas, showing numbers of MRGCs intermediate to those found in the Dscam-/- and Dscam+/- retinas, also had clustering and fasciculation phenotypes that were intermediate to retinas with those genotypes. DACs, by comparison, showed changes in both the Dscam-/- and the Bax-/- retinas that did not correlate with their increases in DAC number. The fasciculation phenotype in the Dscam-/- retina was particularly prominent despite only modest clustering. These results demonstrate that the somal clustering and fasciculation observed in the Dscam mutant retina are not unique to Dscam deficiency and are manifested distinctively by different retinal cell types.

To determine the role of homotypic interactions between neighboring dopaminergic amacrine (DA) cells upon dendritic morphogenesis, the morphology of single cells was examined relative to the positioning of all neighboring homotypic cells. For each labeled cell, the dendritic field was reconstructed, its Voronoi domain was calculated, and the two were related. The dendritic fields of DA cells were observed to be large, sparse, and highly irregular. Dendrites readily overlapped those of neighboring cells, showing no evidence for dendritic tiling or inter-digitation consistent with homotypic repulsion or avoidance. Furthermore, a direct comparison of dendritic field area with the Voronoi domain area of the same cell showed no evidence for dendritic growth being constrained or biased by the local distribution of homotypic neighbors in wild-type retinas. A comparison of the processes of adjacent filled cells confirmed their immediate proximity to one another within the inner plexiform layer, indicating that they do not engage in mutual avoidance by coursing at different depths. Together, these results suggest that the morphogenesis of DA cells is independent of homotypic interactions. However, in the absence of the pro-apoptotic Bax gene, which yields a fourfold increase in DA cell number, a small but significant reduction in dendritic field size was obtained, although not so great as would be predicted by the increase in density. The present results are considered in light of recent studies on the role of cell adhesion molecules expressed by developing DA cells.

Diabetic retinopathy is the most common microvascular complication of diabetes and is a major cause of blindness, but an understanding of the pathogenesis of the disease has been hampered by a lack of accurate animal models. Here, we explore the dynamics of retinal cellular changes in the Nile rat (Arvicanthis niloticus), a carbohydrate-sensitive model for type 2 diabetes. The early retinal changes in diabetic Nile rats included increased acellular capillaries and loss of pericytes that correlated linearly with the duration of diabetes. These vascular changes occurred in the presence of microglial infiltration but in the absence of retinal ganglion cell loss. After a prolonged duration of diabetes, the Nile rat also exhibits a spectrum of retinal lesions commonly seen in the human condition including vascular leakage, capillary non-perfusion, and neovascularization. Our longitudinal study documents a range and progression of retinal lesions in the diabetic Nile rat remarkably similar to those observed in human diabetic retinopathy, and suggests that this model will be valuable in identifying new therapeutic strategies.

Amacrine cells of the retina are conspicuously variable in their morphologies, their population demographics, and their ensuing functions. Vesicular glutamate transporter 3 (VGluT3) amacrine cells are a recently characterized type of amacrine cell exhibiting local dendritic autonomy. The present analysis has examined three features of this VGluT3 population, including their density, local distribution, and dendritic spread, to discern the extent to which these are interrelated, using male and female mice. We first demonstrate that Bax-mediated cell death transforms the mosaic of VGluT3 cells from a random distribution into a regular mosaic. We subsequently examine the relationship between cell density and mosaic regularity across recombinant inbred strains of mice, finding that, although both traits vary across the strains, they exhibit minimal covariation. Other genetic determinants must therefore contribute independently to final cell number and to mosaic order. Using a conditional KO approach, we further demonstrate that Bax acts via the bipolar cell population, rather than cell-intrinsically, to control VGluT3 cell number. Finally, we consider the relationship between the dendritic arbors of single VGluT3 cells and the distribution of their homotypic neighbors. Dendritic field area was found to be independent of Voronoi domain area, while dendritic coverage of single cells was not conserved, simply increasing with the size of the dendritic field. Bax-KO retinas exhibited a threefold increase in dendritic coverage. Each cell, however, contributed less dendrites at each depth within the plexus, intermingling their processes with those of neighboring cells to approximate a constant volumetric density, yielding a uniformity in process coverage across the population.SIGNIFICANCE STATEMENT Different types of retinal neuron spread their processes across the surface of the retina to achieve a degree of dendritic coverage that is characteristic of each type. Many of these types achieve a constant coverage by varying their dendritic field area inversely with the local density of like-type neighbors. Here we report a population of retinal amacrine cells that do not develop dendritic arbors in relation to the spatial positioning of such homotypic neighbors; rather, this cell type modulates the extent of its dendritic branching when faced with a variable number of overlapping dendritic fields to approximate a uniformity in dendritic density across the retina.

The nuclear factor one (NFI) transcription factor genes Nfia, Nfib, and Nfix are all enriched in late-stage retinal progenitor cells, and their loss has been shown to retain these progenitors at the expense of later-generated retinal cell types. Whether they play any role in the specification of those later-generated fates is unknown, but the expression of one of these, Nfia, in a specific amacrine cell type may intimate such a role. Here, Nfia conditional knockout (Nfia-CKO) mice (both sexes) were assessed, finding a massive and largely selective absence of AII amacrine cells. There was, however, a partial reduction in type 2 cone bipolar cells (CBCs), being richly interconnected to AII cells. Counts of dying cells showed a significant increase in Nfia-CKO retinas at postnatal day (P)7, after AII cell numbers were already reduced but in advance of the loss of type 2 CBCs detected by P10. Those results suggest a role for Nfia in the specification of the AII amacrine cell fate and a dependency of the type 2 CBCs on them. Delaying the conditional loss of Nfia to the first postnatal week did not alter AII cell number nor differentiation, further suggesting that its role in AII cells is solely associated with their production. The physiological consequences of their loss were assessed using the ERG, finding the oscillatory potentials to be profoundly diminished. A slight reduction in the b-wave was also detected, attributed to an altered distribution of the terminals of rod bipolar cells, implicating a role of the AII amacrine cells in constraining their stratification.SIGNIFICANCE STATEMENT The transcription factor NFIA is shown to play a critical role in the specification of a single type of retinal amacrine cell, the AII cell. Using an Nfia-conditional knockout mouse to eliminate this population of retinal neurons, we demonstrate two selective bipolar cell dependencies on the AII cells; the terminals of rod bipolar cells become mis-stratified in the inner plexiform layer, and one type of cone bipolar cell undergoes enhanced cell death. The physiological consequence of this loss of the AII cells was also assessed, finding the cells to be a major contributor to the oscillatory potentials in the electroretinogram.

Retinal OFF bipolar cells show distinct connectivity patterns with photoreceptors in the wild-type mouse retina. Some types are cone-specific while others penetrate further through the outer plexiform layer (OPL) to contact rods in addition to cones. To explore dendritic stratification of OFF bipolar cells in the absence of rods, we made use of the 'cone-full' Nrl-/- mouse retina in which all photoreceptor precursor cells commit to a cone fate including those which would have become rods in wild-type retinas. The dendritic distribution of OFF bipolar cell types was investigated by confocal and electron microscopic imaging of immunolabeled tissue sections. The cells' dendrites formed basal contacts with cone terminals and expressed the corresponding glutamate receptor subunits at those sites, indicating putative synapses. All of the four analyzed cell populations showed distinctive patterns of vertical dendritic invasion through the OPL. This disparate behavior of dendritic extension in an environment containing only cone terminals demonstrates type-dependent specificity for dendritic outgrowth in OFF bipolar cells: rod terminals are not required for inducing dendritic extension into distal areas of the OPL.