Genomic Control of Retinal Cell Number: Challenges, Protocol, and Results

Quantitative trait loci on chromosomes 9 and 19 modulate AII amacrine cell number in the mouse retina

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.

Interrelationships between Cellular Density, Mosaic Patterning, and Dendritic Coverage of VGluT3 Amacrine Cells

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.

From random to regular: Variation in the patterning of retinal mosaics

The various types of retinal neurons are each positioned at their respective depths within the retina where they are believed to be assembled as orderly mosaics, in which like-type neurons minimize proximity to one another. Two common statistical analyses for assessing the spatial properties of retinal mosaics include the nearest neighbor analysis, from which an index of their "regularity" is commonly calculated, and the density recovery profile derived from autocorrelation analysis, revealing the presence of an exclusion zone indicative of anti-clustering. While each of the spatial statistics derived from these analyses, the regularity index and the effective radius, can be useful in characterizing such properties of orderly retinal mosaics, they are rarely sufficient for conveying the natural variation in the self-spacing behavior of different types of retinal neurons and the extent to which that behavior generates uniform intercellular spacing across the mosaic. We consider the strengths and limitations of these and other spatial statistical analyses for assessing the patterning in retinal mosaics, highlighting a number of misconceptions and their frequent misuse. Rather than being diagnostic criteria for determining simply whether a population is "regular," they should be treated as descriptive statistics that convey variation in the factors that influence neuronal positioning. We subsequently apply multiple spatial statistics to the analysis of eight different mosaics in the mouse retina, demonstrating conspicuous variability in the degree of patterning present, from essentially random to notably regular. This variability in patterning has both a developmental as well as a functional significance, reflecting the rules governing the positioning of different types of neurons as the architecture of the retina is assembled, and the distinct mechanisms by which they regulate dendritic growth to generate their characteristic coverage and connectivity.

Dopaminergic amacrine cell number, plexus density, and dopamine content in the mouse retina: Strain differences and effects of Bax gene disruption

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 Control of Rod Bipolar Cell Number in the Mouse Retina

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.

DNER and NFIA are expressed by developing and mature AII amacrine cells in the mouse retina

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.

Random spatial patterning of cone bipolar cell mosaics in the mouse retina

Retinal bipolar cells spread their dendritic arbors to tile the retinal surface, extending them to the tips of the dendritic fields of their homotypic neighbors, minimizing dendritic overlap. Such uniform nonredundant dendritic coverage of these populations would suggest a degree of spatial order in the properties of their somal distributions, yet few studies have examined the patterning in retinal bipolar cell mosaics. The present study examined the organization of two types of cone bipolar cells in the mouse retina, the Type 2 cells and the Type 4 cells, and compared their spatial statistical properties with those of the horizontal cells and the cholinergic amacrine cells, as well as to random simulations of cells matched in density and constrained by soma size. The Delauney tessellation of each field was computed, from which nearest neighbor distances and Voronoi domain areas were extracted, permitting a calculation of their respective regularity indexes (RIs). The spatial autocorrelation of the field was also computed, from which the effective radius and packing factor (PF) were determined. Both cone bipolar cell types were found to be less regular and less efficiently packed than either the horizontal cells or cholinergic amacrine cells. Furthermore, while the latter two cell types had RIs and PFs in excess of those for their matched random simulations, the two types of cone bipolar cells had spatial statistical properties comparable to random distributions. An analysis of single labeled cone bipolar cells revealed dendritic arbors frequently skewed to one side of the soma, as would be expected from a randomly distributed population of cells with dendrites that tile. Taken together, these results suggest that, unlike the horizontal cells or cholinergic amacrine cells which minimize proximity to one another, cone bipolar cell types are constrained only by their physical size.

Genomic control of neuronal demographics in the retina

The mature retinal architecture is composed of various types of neuron, each population differing in size and constrained to particular layers, wherein the cells achieve a characteristic patterning in their local organization. These demographic features of retinal nerve cell populations are each complex traits controlled by multiple genes affecting different processes during development, and their genetic determinants can be dissected by correlating variation in these traits with their genomic architecture across recombinant-inbred mouse strains. Using such a resource, we consider how the variation in the numbers of twelve different types of retinal neuron are independent of one another, including those sharing transcriptional regulation as well as those that are synaptically-connected, each mapping to distinct genomic loci. Using the populations of two retinal interneurons, the horizontal cells and the cholinergic amacrine cells, we present in further detail examples where the variation in neuronal number, as well as the variation in mosaic patterning or in laminar positioning, each maps to discrete genomic loci where allelic variants modulating these features must be present. At those loci, we identify candidate genes which, when rendered non-functional, alter those very demographic properties, and in turn, we identify candidate coding or regulatory variants that alter protein structure or gene expression, respectively, being prospective contributors to the variation in phenotype. This forward-genetic approach provides an alternative means for dissecting the molecular genetic control of neuronal population dynamics, with each genomic locus serving as a causal anchor from which we may ultimately understand the developmental principles responsible for the control of those traits.