Pituitary tumor-transforming gene 1 regulates the patterning of retinal mosaics

Pituitary Tumor Transforming Gene 1 Orchestrates Gene Regulatory Variation in Mouse Ventral Midbrain During Aging

Dopaminergic neurons in the midbrain are of particular interest due to their role in diseases such as Parkinson’s disease and schizophrenia. Genetic variation between individuals can affect the integrity and function of dopaminergic neurons but the DNA variants and molecular cascades modulating dopaminergic neurons and other cells types of ventral midbrain remain poorly defined. Three genetically diverse inbred mouse strains – C57BL/6J, A/J, and DBA/2J – differ significantly in their genomes (∼7 million variants), motor and cognitive behavior, and susceptibility to neurotoxins. To further dissect the underlying molecular networks responsible for these variable phenotypes, we generated RNA-seq and ChIP-seq data from ventral midbrains of the 3 mouse strains. We defined 1000–1200 transcripts that are differentially expressed among them. These widespread differences may be due to altered activity or expression of upstream transcription factors. Interestingly, transcription factors were significantly underrepresented among the differentially expressed genes, and only one transcription factor, Pttg1, showed significant differences between all three strains. The changes in Pttg1 expression were accompanied by consistent alterations in histone H3 lysine 4 trimethylation at Pttg1 transcription start site. The ventral midbrain transcriptome of 3-month-old C57BL/6J congenic Pttg1^–/– mutants was only modestly altered, but shifted toward that of A/J and DBA/2J in 9-month-old mice. Principle component analysis (PCA) identified the genes underlying the transcriptome shift and deconvolution of these bulk RNA-seq changes using midbrain single cell RNA-seq data suggested that the changes were occurring in several different cell types, including neurons, oligodendrocytes, and astrocytes. Taken together, our results show that Pttg1 contributes to gene regulatory variation between mouse strains and influences mouse midbrain transcriptome during aging.

Spatiotemporal gene expression patterns reveal molecular relatedness between retinal laminae

In several areas of the central nervous system, neurons are regionally organized into groups or layers that carry out specific activities. In this form of patterning, neurons of distinct types localize their cell bodies to just one or a few of the layers within a structure. However, little is known about whether diverse neuron types within a lamina share molecular features that coordinate their organization. To begin to identify such candidates, we used the laminated murine retina to screen 92 lacZ reporter lines available through the Knockout Mouse Project. Thirty-two of these displayed reporter expression in restricted subsets of inner retina neurons. We then identified the spatiotemporal expression patterns of these genes at key developmental stages. This uncovered several that were heavily enriched in development but reduced in adulthood, including the transcriptional regulator Hmga1. An additional set of genes displayed maturation associated laminar enrichment. Among these, we identified Bbox1 as a novel gene that specifically labels all neurons in the ganglion cell layer but is largely excluded from otherwise molecularly similar neurons in the inner retina. Finally, we established Dbn1 as a new marker enriched in amacrines and Fmnl3 as a marker for subsets of αRGCs. Together, these data provide a spatiotemporal map for laminae-specific molecules and suggest that diverse neuron types within a lamina share coordinating molecular features that may inform their fate or function.

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.

Spatial structure arising from chase-escape interactions with crowding

Movement of individuals, mediated by localised interactions, plays a key role in numerous processes including cell biology and ecology. In this work, we investigate an individual-based model accounting for various intraspecies and interspecies interactions in a community consisting of two distinct species. In this framework we consider one species to be chasers and the other species to be escapees, and we focus on chase-escape dynamics where the chasers are biased to move towards the escapees, and the escapees are biased to move away from the chasers. This framework allows us to explore how individual-level directional interactions scale up to influence spatial structure at the macroscale. To focus exclusively on the role of motility and directional bias in determining spatial structure, we consider conservative communities where the number of individuals in each species remains constant. To provide additional information about the individual-based model, we also present a mathematically tractable deterministic approximation based on describing the evolution of the spatial moments. We explore how different features of interactions including interaction strength, spatial extent of interaction, and relative density of species influence the formation of the macroscale spatial patterns.

Mouse models and strain-dependency of Chédiak-Higashi syndrome-associated neurologic dysfunction

Chédiak-Higashi syndrome (CHS) is a lethal disorder caused by mutations in the LYST gene that involves progressive neurologic dysfunction. Lyst-mutant mice exhibit neurologic phenotypes that are sensitive to genetic background. On the DBA/2J-, but not on the C57BL/6J-background, Lyst-mutant mice exhibit overt tremor phenotypes associated with loss of cerebellar Purkinje cells. Here, we tested whether assays for ataxia could measure this observed strain-dependency, and if so, establish parameters for empowering phenotype- and candidate-driven approaches to identify genetic modifier(s). A composite phenotypic scoring system distinguished phenotypes in Lyst-mutants and uncovered a previously unrecognized background difference between wild-type C57BL/6J and DBA/2J mice. Accelerating rotarod performance also distinguished phenotypes in Lyst-mutants, but at more advanced ages. These results establish that genetic background, Lyst genotype, and age significantly influence the severity of CHS-associated neurologic deficits. Purkinje cell quantifications likewise distinguished phenotypes of Lyst-mutant mice, as well as background differences between wild-type C57BL/6J and DBA/2J mice. To aid identification of potential genetic modifier genes causing these effects, we searched public datasets for cerebellar-expressed genes that are differentially expressed and/or contain potentially detrimental genetic variants. From these approaches, Nos1, Prdx2, Cbln3, Gnb1, Pttg1 were confirmed to be differentially expressed and leading candidates.

Conditional Deletion of AP-2α and AP-2β in the Developing Murine Retina Leads to Altered Amacrine Cell Mosaics and Disrupted Visual Function

Purpose

The combined action of the activating protein-2 (AP-2) transcription factors, AP-2α and AP-2β, is important in early retinal development, specifically in the formation of horizontal cells. However, in previous studies, it was not possible to analyze postnatal development and function of additional retinal subtypes.

Methods

We used a double conditional deletion of AP-2α and AP-2β from the retina to further examine the combinatory role of these genes in retinal cell patterning and function in postnatal adult mice as measured by Voronoi domain area and nearest-neighbor distance spatial analyses and ERGs, respectively.

Results

Conditional deletion of both AP-2α and AP-2β from the retina resulted in a variety of abnormalities, including the absence of horizontal cells, defects in the photoreceptor ribbons in which synapses failed to form, along with evidence of aberrant amacrine cell arrangement. Although no significant changes in amacrine cell population numbers were observed in the double mutants, significant irregularities in the mosaic patterning of amacrine cells was observed as demonstrated by both Voronoi domain areas and nearest-neighbor distances analyses. These changes were further accompanied by an alteration in the retinal response to light as recorded by ERGs. In particular, in the double-mutant mice lacking AP-2α and AP-2β, the b-wave amplitude, representative of interneuron signal processing, was significantly reduced compared with control littermates.

Conclusions

Together these findings demonstrate the requirement for both AP-2α and AP-2β in proper amacrine mosaic patterning and a normal functional light response in the retina.

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.

Gene expression in retinal ischemic post-conditioning

PURPOSE

The pathophysiology of retinal ischemia involves mechanisms including inflammation and apoptosis. Ischemic post-conditioning (Post-C), a brief non-lethal ischemia, induces a long-term ischemic tolerance, but the mechanisms of ischemic post-conditioning in the retina have only been described on a limited basis. Accordingly, we conducted this study to determine the molecular events in retinal ischemic post-conditioning and to identify targets for therapeutic strategies for retinal ischemia.

METHODS

To determine global molecular events in ischemic post-conditioning, a comprehensive study of the transcriptome of whole retina was performed. We utilized RNA sequencing (RNA-Seq), a recently developed, deep sequencing technique enabling quantitative gene expression, with low background noise, dynamic detection range, and discovery of novel genes. Rat retina was subjected to ischemia in vivo by elevation of intraocular pressure above systolic blood pressure. At 24 h after ischemia, Post-C or sham Post-C was performed by another, briefer period of ischemia, and 24 h later, retinas were collected and RNA processed.

RESULTS

There were 71 significantly affected pathways in post-conditioned/ischemic vs. normals and 43 in sham post conditioned/ischemic vs. normals. Of these, 28 were unique to Post-C and ischemia. Seven biological pathways relevant to ischemic injury, in Post-C as opposed to sham Post-C, were examined in detail. Apoptosis, p53, cell cycle, JAK-STAT, HIF-1, MAPK and PI3K-Akt pathways significantly differed in the number as well as degree of fold change in genes between conditions.

CONCLUSION

Post-C is a complex molecular signaling process with a multitude of altered molecular pathways. We identified potential gene candidates in Post-C. Studying the impact of altering expression of these factors may yield insight into new methods for treating or preventing damage from retinal ischemic disorders.

Reproducibility and replicability of rodent phenotyping in preclinical studies

The scientific community is increasingly concerned with the proportion of published “discoveries” that are not replicated in subsequent studies. The field of rodent behavioral phenotyping was one of the first to raise this concern, and to relate it to other methodological issues: the complex interaction between genotype and environment; the definitions of behavioral constructs; and the use of laboratory mice and rats as model species for investigating human health and disease mechanisms. In January 2015, researchers from various disciplines gathered at Tel Aviv University to discuss these issues. The general consensus was that the issue is prevalent and of concern, and should be addressed at the statistical, methodological and policy levels, but is not so severe as to call into question the validity and the usefulness of model organisms as a whole. Well-organized community efforts, coupled with improved data and metadata sharing, have a key role in identifying specific problems and promoting effective solutions. Replicability is closely related to validity, may affect generalizability and translation of findings, and has important ethical implications.

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.

Loss of CPEB3 Upregulates MEGF10 to Impair Mosaic Development of ON Starburst Amacrine Cells

Cytoplasmic polyadenylation element binding protein 3 (CPEB3) regulates target RNA translation in neurons. Here, we examined CPEB3 distribution and function in the mouse retina. CPEB3 is expressed in retinal neurons, including those located in the inner nuclear layer (INL) and ganglion cell layer (GCL) but not in cone and rod photoreceptors in the outer nuclear layer (ONL). A previous study found CPEB3 expressed in cholinergic starburst amacrine cells (SACs). We first examined these cells and observed aberrant SAC mosaicism in CPEB3-knockout (KO) retinas. Retinal neurons showed orderly spatial arrangements. Many individual subtypes are organized non-randomly in patterns called mosaics. Despite CPEB3 being expressed in both populations of SACs, OFF SACs in the INL and ON SACs in the GCL, aberrant mosaic regularity was observed in only ON SACs of CPEB3-KO retinas. Molecular characterization revealed that translation of multiple epidermal growth factor 10 (Megf10) RNA is suppressed by CPEB3 during the first week of postnatal development, when MEGF10 is primarily expressed in SACs and mediates homotypic repulsive interactions to define intercellular spacing of SACs. Thus, elevated MEGF10 expression in the absence of the translational repressor CPEB3 may account for the defective spatial organization of ON SACs. Our findings uncover for the first time that translational control plays a role in shaping retinal mosaic arrangement.

Wide-field diffuse amacrine cells in the monkey retina contain immunoreactive Cocaine- and Amphetamine-Regulated Transcript (CART)

The goals of this study were to localize the neuropeptide Cocaine- and Amphetamine-Regulated Transcript (CART) in primate retinas and to describe the morphology, neurotransmitter content and synaptic connections of the neurons that contain it. Using in situ hybridization, light and electron microscopic immunolabeling, CART was localized to GABAergic amacrine cells in baboon retinas. The CART-positive cells had thin, varicose dendrites that gradually descended through the inner plexiform layer and ramified extensively in the innermost stratum. They resembled two types of wide-field diffuse amacrine cells that had been described previously in macaque retinas using the Golgi method and also A17, serotonin-accumulating and waterfall cells of other mammals. The CART-positive cells received synapses from rod bipolar cell axons and made synapses onto the axons in a reciprocal configuration. The CART-positive cells also received synapses from other amacrine cells. Some of these were located on their primary dendrites, and the presynaptic cells there included dopaminergic amacrine cells. Although some CART-positive somas were localized in the ganglion cell layer, they did not contain the ganglion cell marker RNA binding protein with multiple splicing (RBPMS). Based on these results and electrophysiological studies in other mammals, the CART-positive amacrine cells would be expected to play a major role in the primary rod pathway of primates, providing feedback inhibition to rod bipolar cells.

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.

Spatial moment dynamics for collective cell movement incorporating a neighbour-dependent directional bias

The ability of cells to undergo collective movement plays a fundamental role in tissue repair, development and cancer. Interactions occurring at the level of individual cells may lead to the development of spatial structure which will affect the dynamics of migrating cells at a population level. Models that try to predict population-level behaviour often take a mean-field approach, which assumes that individuals interact with one another in proportion to their average density and ignores the presence of any small-scale spatial structure. In this work, we develop a lattice-free individual-based model (IBM) that uses random walk theory to model the stochastic interactions occurring at the scale of individual migrating cells. We incorporate a mechanism for local directional bias such that an individual's direction of movement is dependent on the degree of cell crowding in its neighbourhood. As an alternative to the mean-field approach, we also employ spatial moment theory to develop a population-level model which accounts for spatial structure and predicts how these individual-level interactions propagate to the scale of the whole population. The IBM is used to derive an equation for dynamics of the second spatial moment (the average density of pairs of cells) which incorporates the neighbour-dependent directional bias, and we solve this numerically for a spatially homogeneous case.

Spatial structure arising from neighbour-dependent bias in collective cell movement

Mathematical models of collective cell movement often neglect the effects of spatial structure, such as clustering, on the population dynamics. Typically, they assume that individuals interact with one another in proportion to their average density (the mean-field assumption) which means that cell-cell interactions occurring over short spatial ranges are not accounted for. However, in vitro cell culture studies have shown that spatial correlations can play an important role in determining collective behaviour. Here, we take a combined experimental and modelling approach to explore how individual-level interactions give rise to spatial structure in a moving cell population. Using imaging data from in vitro experiments, we quantify the extent of spatial structure in a population of 3T3 fibroblast cells. To understand how this spatial structure arises, we develop a lattice-free individual-based model (IBM) and simulate cell movement in two spatial dimensions. Our model allows an individual's direction of movement to be affected by interactions with other cells in its neighbourhood, providing insights into how directional bias generates spatial structure. We consider how this behaviour scales up to the population level by using the IBM to derive a continuum description in terms of the dynamics of spatial moments. In particular, we account for spatial correlations between cells by considering dynamics of the second spatial moment (the average density of pairs of cells). Our numerical results suggest that the moment dynamics description can provide a good approximation to averaged simulation results from the underlying IBM. Using our in vitro data, we estimate parameters for the model and show that it can generate similar spatial structure to that observed in a 3T3 fibroblast cell population.

Identification of candidate genes that underlie the QTL on chromosome 1 that mediates genetic differences in stress-ethanol interactions

Alcoholism, stress, and anxiety are strongly interacting heritable, polygenetic traits. In a previous study, we identified a quantitative trait locus (QTL) on murine chromosome (Chr) 1 between 23.0 and 31.5 Mb that modulates genetic differences in the effects of ethanol on anxiety-related phenotypes. The goal of the present study was to extend the analysis of this locus with a focus on identifying candidate genes using newly available data and tools. Anxiety-like behavior was evaluated with an elevated zero maze following saline or ethanol injections (1.8 g/kg) in C57BL/6J, DBA2J, and 72 BXD strains. We detected significant effects of strain and treatment and their interaction on anxiety-related behaviors, although surprisingly, sex was not a significant factor. The Chr1 QTL is specific to the ethanol-treated cohort. Candidate genes in this locus were evaluated using now standard bioinformatic criteria. Collagen 19a1 (Col19a1) and family sequence 135a (Fam135a) met most criteria but have lower expression levels and lacked biological verification and, therefore, were considered less likely candidates. In contrast, two other genes, the prenylated protein tyrosine phosphate family member Ptp4a1 (protein tyrosine phosphate 4a1) and the zinc finger protein Phf3 (plant homeoDomain finger protein 3) met each of our bioinformatic criteria and are thus strong candidates. These findings are also of translational relevance because both Ptp4a1 and Phf3 have been nominated as candidates genes for alcohol dependence in a human genome-wide association study. Our findings support the hypothesis that variants in one or both of these genes modulate heritable differences in the effects of ethanol on anxiety-related behaviors.

The patterning of retinal horizontal cells: normalizing the regularity index enhances the detection of genomic linkage

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.