Optimizing retinal ganglion cell nuclear staining for automated cell counting.

The retinal ganglion cells (RGCs) serve as the critical pathway for transmitting visual information from the retina to the brain, yet they can be dramatically impacted by diseases such as glaucoma. When investigating disease processes affecting RGCs in mouse models, accurately quantifying affected cells becomes essential. However, the use of pan RGC markers like RBPMS or THY1 presents challenges in accurate total cell counting. While Brn3a serves as a reliable RGC nuclear marker for automated counting, it fails to encompass all RGC subtypes in mice. To address this limitation and enable precise automated counting, our research endeavors to develop a method for labeling nuclei in all RGC subtypes. Investigating RGC subtypes labeled with the nuclear marker POU6F2 revealed that numerous RGCs unlabeled by Brn3a were, in fact, labeled with POU6F2. We hypothesize that using antibodies against both Brn3a and POU6F2 would label virtually all RGC nuclei in the mouse retina. Our experiments confirmed that staining retinas with both markers resulted in the labeling of all RGCs. Additionally, when using the cell body marker RBPMS known to label all mouse RGCs, all RBPMS-labeled cells also exhibited Brn3a or POU6F2 labeling. This combination of Brn3a and POU6F2 antibodies provides a pan-RGC nuclear stain, facilitating accurate automated counting by labeling cell nuclei in the retina.

Commonalities of optic nerve injury and glaucoma-induced neurodegeneration: Insights from transcriptome-wide studies.

Glaucoma is a collection of diseases that lead to an irreversible vision loss due to damage of retinal ganglion cells (RGCs). Although the underlying events leading to RGC death are not fully understood, recent research efforts are beginning to define the genetic changes that play a critical role in the initiation and progression of glaucomatous injury and RGC death. Several genetic and experimental animal models have been developed to mimic glaucomatous neurodegeneration. These models differ in many respects but all result in the loss of RGCs. Assessing transcriptional changes across different models could provide a more complete perspective on the molecular drivers of RGC degeneration. For the past several decades, changes in the retinal transcriptome during neurodegeneration process were defined using microarray methods, RNA sequencing and now single cell RNA sequencing. It is understood that these methods have strengths and weaknesses due to technical differences and variations in the analytical tools used. In this review, we focus on the use of transcriptome-wide expression profiling of the changes occurring as RGCs are lost across different glaucoma models. Commonalities of optic nerve crush and glaucoma-induced neurodegeneration are identified and discussed.

Genomic locus modulating corneal thickness in the mouse identifies POU6F2 as a potential risk of developing glaucoma.

Central corneal thickness (CCT) is one of the most heritable ocular traits and it is also a phenotypic risk factor for primary open angle glaucoma (POAG). The present study uses the BXD Recombinant Inbred (RI) strains to identify novel quantitative trait loci (QTLs) modulating CCT in the mouse with the potential of identifying a molecular link between CCT and risk of developing POAG. The BXD RI strain set was used to define mammalian genomic loci modulating CCT, with a total of 818 corneas measured from 61 BXD RI strains (between 60-100 days of age). The mice were anesthetized and the eyes were positioned in front of the lens of the Phoenix Micron IV Image-Guided OCT system or the Bioptigen OCT system. CCT data for each strain was averaged and used to QTLs modulating this phenotype using the bioinformatics tools on GeneNetwork (www.genenetwork.org). The candidate genes and genomic loci identified in the mouse were then directly compared with the summary data from a human POAG genome wide association study (NEIGHBORHOOD) to determine if any genomic elements modulating mouse CCT are also risk factors for POAG.This analysis revealed one significant QTL on Chr 13 and a suggestive QTL on Chr 7. The significant locus on Chr 13 (13 to 19 Mb) was examined further to define candidate genes modulating this eye phenotype. For the Chr 13 QTL in the mouse, only one gene in the region (Pou6f2) contained nonsynonymous SNPs. Of these five nonsynonymous SNPs in Pou6f2, two resulted in changes in the amino acid proline which could result in altered secondary structure affecting protein function. The 7 Mb region under the mouse Chr 13 peak distributes over 2 chromosomes in the human: Chr 1 and Chr 7. These genomic loci were examined in the NEIGHBORHOOD database to determine if they are potential risk factors for human glaucoma identified using meta-data from human GWAS. The top 50 hits all resided within one gene (POU6F2), with the highest significance level of p = 10-6 for SNP rs76319873. POU6F2 is found in retinal ganglion cells and in corneal limbal stem cells. To test the effect of POU6F2 on CCT we examined the corneas of a Pou6f2-null mice and the corneas were thinner than those of wild-type littermates. In addition, these POU6F2 RGCs die early in the DBA/2J model of glaucoma than most RGCs. Using a mouse genetic reference panel, we identified a transcription factor, Pou6f2, that modulates CCT in the mouse. POU6F2 is also found in a subset of retinal ganglion cells and these RGCs are sensitive to injury.

Using BXD mouse strains in vision research: A systems genetics approach.

We illustrate the growing power of the BXD family of mice (recombinant inbred strains from a cross of C57BL/6J and DBA/2J mice) and companion bioinformatic tools to study complex genome-phenome relations related to glaucoma. Over the past 16 years, our group has integrated powerful murine resources and web-accessible tools to identify networks modulating visual system traits-from photoreceptors to the visual cortex. Recent studies focused on retinal ganglion cells and glaucoma risk factors, including intraocular pressure (IOP), central corneal thickness (CCT), and susceptibility of cellular stress. The BXD family was exploited to define key gene variants and then establish linkage to glaucoma in human cohorts. The power of this experimental approach to precision medicine is highlighted by recent studies that defined cadherin 11 (Cdh11) and a calcium channel (Cacna2d1) as genes modulating IOP, Pou6f2 as a genetic link between CCT and retinal ganglion cell (RGC) death, and Aldh7a1 as a gene that modulates the susceptibility of RGCs to death after elevated IOP. The role of three of these gene variants in glaucoma is discussed, along with the pathways activated in the disease process.

RNA sequencing profiling of the retina in C57BL/6J and DBA/2J mice: Enhancing the retinal microarray data sets from GeneNetwork.

Purpose: The goal of the present study is to provide an independent assessment of the retinal transcriptome signatures of C57BL/6J (B6) and DBA/2J (D2) mice, and to enhance existing microarray data sets for accurately defining the allelic differences in the BXD recombinant inbred strains. Methods: Retinas from B6 and D2 mice (three of each) were used for the RNA sequencing (RNA-seq) analysis. Transcriptome features were examined for both strains. Differentially expressed genes between the two strains were identified, and bioinformatic analysis was performed to analyze the transcriptome differences between the B6 and D2 strains, including Gene Ontology (GO) analysis, Phenotype and Reactome enrichment, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The RNA-seq data were then directly compared with one of the microarray data sets (Department of Defense [DoD] Retina Normal Affy MoGene 2.0 ST RMA Gene Level Microarray Database) hosted on GeneNetwork. Results: RNA-seq provided an in-depth analysis of the transcriptome of the B6 and D2 retinas with a total of more than 30,000,000 reads per sample. More than 70% of the reads were uniquely mapped, resulting in a total of 18,100 gene counts for all six samples. A total of 1,665 genes were differentially expressed, with 858 of these more highly expressed in the B6 retinas and 807 more highly expressed in the D2 retinas. Several molecular pathways were differentially active between the two strains, including the retinoic acid metabolic process, endoplasmic reticulum lumen, extracellular matrix (ECM) organization, and the PI3K-Akt signaling pathway. The most enriched KEGG pathways were the pentose and glucuronate interconversions pathway, the cytochrome P450 pathway, the protein digestion and absorption pathway, and the ECM-receptor interaction pathway. Each of these pathways had a more than fourfold enrichment. The DoD Normal Retina Microarray Database provided expression profiling for 26,191 annotated transcripts for B6 mouse, D2 mouse, and 53 BXD strains. A total of 13,793 genes in this microarray data set were comparable to the RNA-seq data set. For the B6 and D2 retinas, the RNA-seq data and the microarray data were highly correlated with each other (Pearson's r=0.780 for the B6 mice and 0.784 for D2 mice). These results suggest that the microarray data set can reliably detect differentially expressed genes between the B6 and D2 retinas, with an overall accuracy of 91.1%. Examples of true positive and false positive genes are provided. Conclusions: Retinal transcriptome features of B6 and D2 mouse strains provide a useful reference for a better understanding of the mouse retina. Generally, the microarray database presented on GeneNetwork shows good agreement with the RNA-seq data, but we note that any allelic difference between B6 and D2 mice should be verified with the latter.

Optic nerve regeneration in the mouse is a complex trait modulated by genetic background.

Purpose: The present study is designed to identify the influences of genetic background on optic nerve regeneration using the two parental strains (C57BL/6J and DBA/2J) and seven BXD recombinant inbred mouse strains. Methods: To study regeneration in the optic nerve, Pten was knocked down in the retinal ganglion cells using adenoassociated virus (AAV) delivery of shRNA, and a mild inflammatory response was induced with an intravitreal injection of zymosan with CPT-cAMP. The axons of the retinal ganglion cells were damaged by optic nerve crush (ONC). Following a 12-day survival period, regenerating axons were labeled by cholera toxin B, and 2 days later, the regenerating axons within the optic nerve were examined. The number of axons at 0.5 mm and 1 mm from the crush site were counted. In addition, we measured the distance that five axons had grown down the nerve and the longest distance a single axon reached. Results: The analysis revealed a considerable amount of differential axonal regeneration across the seven BXD strains and the parental strains. There was a statistically significant difference (p=0.014 Mann-Whitney U test) in the regenerative capacity in the number of axons reaching 0.5 mm from a low of 236.1±24.4 axons in the BXD102 mice to a high of 759.8±79.2 axons in the BXD29 mice. There were also statistically significant differences (p=0.014 Mann-Whitney U test) in the distance axons traveled. Looking at a minimum of five axons, the shortest distance was 787.2±46.5 µm in the BXD102 mice, and the maximum distance was 2025.5±223.3 µm in the BXD29 mice. Conclusions: Differences in genetic background can have a profound effect on axonal regeneration causing a threefold increase in the number of regenerating axons at 0.5 mm from the crush site and a 2.5-fold increase in the distance traveled by at least five axons in the damaged optic nerve.

The genetic dissection of Myo7a gene expression in the retinas of BXD mice.

Purpose: Usher syndrome (US) is characterized by a loss of vision due to retinitis pigmentosa (RP) and deafness. US has three clinical subtypes, but even within each subtype, the severity varies. Myosin VIIA, coded by Myo7a, has been identified as one of the causal genes of US. This study aims to identify pathways and other genes through which Myo7a interacts to affect the presentation of US symptoms. Methods: In this study, we used the retinal tissue of BXD recombinant inbred (RI) mice to examine the expression of Myo7a and perform genetic mapping. Expression quantitative trait locus (eQTL), single nucleotide polymorphism (SNP), and gene correlation analysis were performed using GeneNetwork. Gene set enrichment analysis was performed using WebGestalt, and gene network construction was performed using the Gene Cohesion Analysis Tool. Results: We found Myo7a to be cis-regulated, with varied levels of expression across BXD strains. Here, we propose a genetic network with 40 genes whose expression is highly correlated with Myo7a. Among these genes, six have been linked to retinal diseases, three to deafness, and five share a transcription factor with Myo7a. Gene ontology and pathway analysis revealed a strong connection among ion channel activity, Myo7a, and US. Conclusions: Although Myo7a is a causal gene of US type I, this gene works with many other genes and pathways to affect the severity of US. Many of the genes found in the genetic network, pathways, and gene ontology categories of Myo7a are related to either deafness or blindness. Further investigation is needed to examine the specific relationships between these genes, which may assist in the treatment of US.

Genomic loci modulating retinal ganglion cell death following elevated IOP in the mouse.

The present study was designed to identify genomic loci modulating the susceptibility of retinal ganglion cells (RGC) to elevated intraocular pressure (IOP) in the BXD recombinant inbred mouse strain set. IOP was elevated by injecting magnetic microspheres into the anterior chamber and blocking the trabecular meshwork using a handheld magnet to impede drainage. The IOP was then measured over the next 21 days. Only animals with IOP greater than 25 mmHg for two consecutive days or an IOP above 30 mmHg on a single day after microsphere-injection were used in this study. On day 21, mice were sacrificed and the optic nerve was processed for histology. Axons were counted for both the injected and the control eye in 49 BXD strains, totaling 181 normal counts and 191 counts associated with elevated IOP. The axon loss for each strain was calculated and the data were entered into genenetwork.org. The average number of normal axons in the optic nerve across all strains was 54,788 ±â€¯16% (SD), which dropped to 49,545 ±â€¯20% in animals with artificially elevated IOP. Interval mapping demonstrated a relatively similar genome-wide map for both conditions with a suggestive Quantitative Trait Locus (QTL) on proximal Chromosome 3. When the relative axon loss was used to generate a genome-wide interval map, we identified one significant QTL (p < 0.05) on Chromosome 18 between 53.6 and 57 Mb. Within this region, the best candidate gene for modulating axon loss was Aldh7a1. Immunohistochemistry demonstrated ALDH7A1 expression in mouse RGCs. ALDH7A1 variants were not significantly associated with glaucoma in the NEIGHBORHOOD GWAS dataset, but this enzyme was identified as part of the butanoate pathway previously associated with glaucoma risk. Our results suggest that genomic background influences susceptibility to RGC degeneration and death in an inducible glaucoma model.

Differential Exon Expression in a Large Family of Retinal Genes Is Regulated by a Single Trans Locus.

Transcription and RNA processing can generate many variant mRNAs (isoforms) from a given genomic locus. The more we learn about RNA processing the more we realize how complex it can be. Examining the expression profiles of individual exons, we observed that specific exons were differentially expressed across a large number of genes in mice. We found that each isoform or exon is independently expressed compared to other exons from the same gene and regulated separately in trans. Each trans locus was identified by mapping using linkage analysis in a large mouse recombinant inbred strain set. We present evidence for a limited number of these master regulatory loci in the retina. One major locus controls about half the expression of the individual exons and resides on Chromosome 4, between 133 and 136 Mb.

Transcriptome networks in the mouse retina: An exon level BXD RI database.

PURPOSE: Differences in gene expression provide diverse retina phenotypes and may also contribute to susceptibility to injury and disease. The present study defines the transcriptome of the retina in the BXD RI strain set, using the Affymetrix Mouse Gene 2.0 ST array to investigate all exons of traditional protein coding genes, non-coding RNAs, and microRNAs. These data are presented in a highly interactive database on the GeneNetwork website. METHODS: In the Normal Retina Database, the mRNA levels of the transcriptome from retinas was quantified using the Affymetrix Mouse Gene 2.0 ST array. This database consists of data from male and female mice. The data set includes a total of 52 BXD RI strains, the parental strains (C57BL/6J and DBA/2J), and a reciprocal cross. RESULTS: In combination with GeneNetwork, the Department of Defense (DoD) Congressionally Directed Medical Research Programs (CDMRP) Normal Retina Database provides a large resource for mapping, graphing, analyzing, and testing complex genetic networks. Protein-coding and non-coding RNAs can be used to map quantitative trait loci (QTLs) that contribute to expression differences among the BXD strains and to establish links between classical ocular phenotypes associated with differences in the genomic sequence. Using this resource, we extracted transcriptome signatures for retinal cells and defined genetic networks associated with the maintenance of the normal retina. Furthermore, we examined differentially expressed exons within a single gene. CONCLUSIONS: The high level of variation in mRNA levels found among the BXD RI strains makes it possible to identify expression networks that underline differences in retina structure and function. Ultimately, we will use this database to define changes that occur following blast injury to the retina.

ImagePAD, a novel counting application for the Apple iPad, used to quantify axons in the mouse optic nerve.

The present article introduces a new and easy to use counting application for the Apple iPad. The application "ImagePAD" takes advantage of the advanced user interface features offered by the Apple iOS platform, simplifying the rather tedious task of quantifying features in anatomical studies. For example, the image under analysis can be easily panned and zoomed using iOS-supported multi-touch gestures without losing the spatial context of the counting task, which is extremely important for ensuring count accuracy. This application allows one to quantify up to 5 different types of objects in a single field and output the data in a tab-delimited format for subsequent analysis. We describe two examples of the use of the application: quantifying axons in the optic nerve of the C57BL/6J mouse and determining the percentage of cells labeled with NeuN or ChAT in the retinal ganglion cell layer. For the optic nerve, contiguous images at 60× magnification were taken and transferred onto an Apple iPad. Axons were counted by tapping on the touch-sensitive screen using ImagePAD. Nine optic nerves were sampled and the number of axons in the nerves ranged from 38,872 axons to 50,196 axons with an average of 44,846 axons per nerve (SD = 3980 axons).

POU6F2, a risk factor for glaucoma, myopia and dyslexia, labels specific populations of retinal ganglion cells.

Pou6f2 is a genetic connection between central corneal thickness (CCT) in the mouse and a risk factor for developing primary open-angle glaucoma. POU6F2 is also a risk factor for several conditions in humans, including glaucoma, myopia, and dyslexia. Recent findings demonstrate that POU6F2-positive retinal ganglion cells (RGCs) comprise a number of RGC subtypes in the mouse, some of which also co-stain for Cdh6 and Hoxd10. These POU6F2-positive RGCs appear to be novel of ON-OFF directionally selective ganglion cells (ooDSGCs) that do not co-stain with CART or SATB2 (typical ooDSGCs markers). These POU6F2-positive cells are sensitive to damage caused by elevated intraocular pressure. In the DBA/2J mouse glaucoma model, heavily-labeled POU6F2 RGCs decrease by 73% at 8 months of age compared to only 22% loss of total RGCs (labeled with RBPMS). Additionally, Pou6f2-/- mice suffer a significant loss of acuity and spatial contrast sensitivity along with an 11.4% loss of total RGCs. In the rhesus macaque retina, POU6F2 labels the large parasol ganglion cells that form the magnocellular (M) pathway. The association of POU6F2 with the M-pathway may reveal in part its role in human glaucoma, myopia, and dyslexia.

A crystallin gene network in the mouse retina.

The present study was designed to examine the regulation of crystallin genes and protein in the mouse retina using the BXD recombinant inbred (RI) strains. Illumina Sentrix BeadChip Arrays (MouseWG-6v2) were used to analyze mRNA levels in 75 BXD RI strains along with the parental strains (C57Bl/6J and DBA/2J), and the reciprocal crosses in the Hamilton Eye Institute (HEI) Retina Dataset (www.genenetwork.org). Protein levels were investigated using immunoblots to quantify levels of proteins and indirect immunohistochemistry to define the distribution of protein. Algorithms in the Genomatix program were used to identify transcription factor binding sites common to the regulatory sequences in the 5' regions of co-regulated set of crystallin and other genes as compared to a set of control genes. As subset of genes, including many encoding lens crystallins is part of a tightly co-regulated network that is active in the retina. Expression of this crystallin network appears to be binary in nature, being expressed either at relatively low levels or being highly upregulated. Relative to a control set of genes, the 5' regulatory sequences of the crystallin network genes show an increased frequency of a set of common transcription factor-binding sites, the most common being those of the Maf family. Chromatin immunoprecipitation of human lens epithelial cells (HLEC) and rat retinal ganglion cells (RGC) confirmed the functionality of these sites, showing that MafA binds the predicted sites of CRYGA and CRYGD in HLE and CRYAB, CRYGA, CRYBA1, and CRYBB3 in RGC cells. In the retina there is a highly correlated group of genes containing many members of the α- ß- and γ-crystallin families. These genes can be dramatically upregulated in the retina. One transcription factor that appears to be involved in this coordinated expression is the MAF family transcription of factors associated with both lens and extralenticular expression of crystallin genes.

A practical approach to optic nerve crush in the mouse.

Our goal is to provide an instructional resource to help others wishing to use the optic nerve crush (ONC) as an experimental procedure. The process is described beginning with the anesthesia, followed by positioning of the mouse, the surgery itself, and post-surgical care. We analyzed the effect of ONC on retinal blood flow, using optical coherence tomography doppler. This procedure produces a consistent loss of cells in the ganglion cell layer, using whole mounts of retina stained with TO-PRO-3. An instructional video is presented that demonstrates a simple surgical approach to effectively crush the optic nerve of the mouse.

Dose-dependent treatment of optic nerve crush by exogenous systemic mutant erythropoietin.

The goal of the present study was to determine the minimum concentration of systemic erythropoietin-R76E required for neuroprotection in the retina. Erythropoietin (EPO) exhibits neuroprotective effects in both in vitro and in vivo models of neuronal cell death although its classical function is the regulation of red blood cell production. It can cross the blood brain barrier and therefore can be delivered systemically to affect the retina. However, long-term treatment with exogenous erythropoietin causes polycythemia. To decrease this potentially lethal effect, we generated and tested a modified form that contains a single arginine to glutamate mutation at the 76th position (EPO-R76E). In previous studies, this mutant protected retinal neurons in mouse models of retinal degeneration and glaucoma with similar efficacy as wild-type EPO. However, EPO-R76E has attenuated erythropoietic activity, therefore, neuroprotection can be achieved without causing a significant rise in hematocrit. BALB/cByJ mice received a single intramuscular injection of recombinant adeno-associated virus carrying enhanced green fluorescent protein, Epo, or Epo-R76E. To result in continuous production of four different doses of EPO-R76E, two doses of two different serotypes (2/5 and 2/8) were used. Mice were subjected to optic nerve crush and analysis was performed thirty days later. EPO-R76E showed dose-dependent protection of the retinal ganglion cell bodies, but was unable to prevent axonal degeneration. Furthermore, EPO-R76E induced a dose-dependent rise in the hematocrit that was still attenuated as compared to wild-type EPO.

A mouse model of ocular blast injury that induces closed globe anterior and posterior pole damage.

We developed and characterized a mouse model of primary ocular blast injury. The device consists of: a pressurized air tank attached to a regulated paintball gun with a machined barrel; a chamber that protects the mouse from direct injury and recoil, while exposing the eye; and a secure platform that enables fine, controlled movement of the chamber in relation to the barrel. Expected pressures were calculated and the optimal pressure transducer, based on the predicted pressures, was positioned to measure output pressures at the location where the mouse eye would be placed. Mice were exposed to one of three blast pressures (23.6, 26.4, or 30.4 psi). Gross pathology, intraocular pressure, optical coherence tomography, and visual acuity were assessed 0, 3, 7, 14, and 28 days after exposure. Contralateral eyes and non-blast exposed mice were used as controls. We detected increased damage with increased pressures and a shift in the damage profile over time. Gross pathology included corneal edema, corneal abrasions, and optic nerve avulsion. Retinal damage was detected by optical coherence tomography and a deficit in visual acuity was detected by optokinetics. Our findings are comparable to those identified in Veterans of the recent wars with closed eye injuries as a result of blast exposure. In summary, this is a relatively simple system that creates injuries with features similar to those seen in patients with ocular blast trauma. This is an important new model for testing the short-term and long-term spectrum of closed globe blast injuries and potential therapeutic interventions.

EDL-291, a novel isoquinoline, presents antiglioblastoma effects in vitro and in vivo.

To investigate the effectiveness of EDL-291, a 6,7-dimethoxy-1-[4-(4-methoxypyridin-3-yl)benzyl]-1,2,3,4-tetrahydroisoquinoline dihydrochloride compound, in inhibiting the survival of glioblastoma in vitro and in vivo. Dose-response curves were generated to determine the EC50 in rat and human glioblastoma cell lines by treatment with different dilutions of EDL-291. To evaluate the architecture of the glioblastoma cells after treatment with EDL-291, the rat and human glioblastoma cells were stained with Mito Tracker Green FM. To determine whether autophagy was induced in EDL-291-treated glioblastoma cells, both rat and human glioblastoma cell lines were stained with acridine orange and light chain-3 immunoblots were performed. The efficacy of EDL-291 was monitored in vivo using a rat glioblastoma model. Rat glioblastoma cells were transplanted into an intracranial rat model, followed by infusions of saline, a low dose of EDL-291 (20 mg/kg for the first half hour, followed by 40 mg/kg EDL-291 in saline for 4 h), or a high dose of EDL-291 (60 mg/kg for the first half hour, followed by 90 mg/kg EDL-291 for 4 h). EDL-291 inhibits glioblastoma in vitro by destroying the mitochondria as shown with Mito Tracker Green FM. Acridine orange staining and light chain-3 immunoblots suggest that autophagy is induced when glioblastoma cells are treated with EDL-291. In vivo, a low dosage of EDL-291 is sufficient and effective in reducing glioblastoma tumor size. EDL-291 selectively induces cell death in rat and human glioblastoma cell lines by the induction of autophagy. EDL-291 exhibits antiglioblastoma effects both in vitro and in vivo.

Genetic networks in the mouse retina: growth associated protein 43 and phosphatase tensin homolog network.

PURPOSE: The present study examines the structure and covariance of endogenous variation in gene expression across the recently expanded family of C57BL/6J (B) X DBA/2J (D) Recombinant Inbred (BXD RI) strains of mice. This work is accompanied by a highly interactive database that can be used to generate and test specific hypotheses. For example, we define the genetic network regulating growth associated protein 43 (Gap43) and phosphatase tensin homolog (Pten). METHODS: The Hamilton Eye Institute (HEI) Retina Database within GeneNetwork features the data analysis of 346 Illumina Sentrix BeadChip Arrays (mouse whole genome-6 version 2). Eighty strains of mice are presented, including 75 BXD RI strains, the parental strains (C57BL/6J and DBA/2J), the reciprocal crosses, and the BALB/cByJ mice. Independent biologic samples for at least two animals from each gender were obtained with a narrow age range (48 to 118 days). Total RNA was prepared followed by the production of biotinylated cRNAs, which were pipetted into the Mouse WG-6V2 arrays. The data was globally normalized with rank invariant and stabilization (2z+8). RESULTS: The HEI Retina Database is located on the GeneNetwork website. The database was used to extract unique transcriptome signatures for specific cell types in the retina (retinal pigment epithelial, amacrine, and retinal ganglion cells). Two genes associated with axonal outgrowth (Gap43 and Pten) were used to display the power of this new retina database. Bioinformatic tools located within GeneNetwork in conjunction with the HEI Retina Database were used to identify the unique signature Quantitative Trait Loci (QTLs) for Gap43 and Pten on chromosomes 1, 2, 12, 15, 16, and 19. Gap43 and Pten possess networks that are similar to ganglion cell networks that may be associated with axonal growth in the mouse retina. This network involves high correlations of transcription factors (SRY sex determining region Y-box 2 [Sox2], paired box gene 6 [Pax6], and neurogenic differentiation 1 [Neurod1]), and genes involved in DNA binding (proliferating cell nuclear antigen [Pcna] and zinc finger, BED-type containing 4 [Zbed4]), as well as an inhibitor of DNA binding (inhibitor of DNA binding 2, dominant negative helix-loop-helix protein [Id2]). Furthermore, we identified the potential upstream modifiers on chromosome 2 (teashirt zinc finger homeobox 2 [Tshz2], RNA export 1 homolog [Rae1] and basic helix-loop-helix domain contatining, class B4 [Bhlhb4]) on chromosome 15 (RAB, member of RAS oncogene family-like 2a [Rabl2a], phosphomannomutase 1 [Pmm1], copine VIII [Cpne8], and fibulin 1 [Fbln1]). CONCLUSIONS: The endogenous variation in mRNA levels among BXD RI strains can be used to explore and test expression networks underlying variation in retina structure, function, and disease susceptibility. The Gap43 and Pten network highlights the covariance of gene expression and forms a molecular network associated with axonal outgrowth in the adult retina.

Complex interactions of Tyrp1 in the eye.

PURPOSE: To use a systems genetics approach to construct and analyze co-expression networks that are causally linked to mutations in a key pigementation gene, tyrosinase-related protein 1 (Tyrp1), that is associated both with oculocutaneous albinism type 3 (OCA3) in humans and with glaucoma in mice. METHODS: Gene expression patterns were measured in whole eyes of a large family of BXD recombinant inbred (RI) mice derived from parental lines that encode for wildtype (C57BL/6J) and mutant (DBA/2J) Tyrp1. Protein levels of Tyrp1 were measured in whole eyes and isolated irides. Bioinformatics analyses were performed on the expression data along with our archived sequence data. Separate data sets were generated which were comprised of strains that harbor either wildtype or mutant Tyrp1 and each was mined individually to identify gene networks that covary significantly with each isoform of Tyrp1. Ontology trees and network graphs were generated to probe essential function and statistical significance of covariation. Genes with strong covariance in wildtype mice were assembled into genome-wide heatmaps for cohorts carrying either wildtype or mutant Tyrp1. RESULTS: Single nucleotide polymorphism (SNP) analysis verified the presence of the Tyrp1b mutation in the Tyrp1 gene. Message levels were greater in BXD strains with the mutant Tyrp1. Interval mapping of these BXD mice revealed a strong expression quantitative trait locus (eQTL) on Chr 4 at the location of the gene itself. Composite mapping revealed a suggestive eQTL on Chr 9 at the location of myosin-Va (Myo5a), mutations in which are known as dilute. Enriched biologic processes associated with wildtype Tyrp1 included pigmentation, melanin biosynthetic process, and mesenchymal cell development, while associations with the mutant gene included categories of neural crest cell development, protein metabolic processes and glycoprotein metabolic processes. Genome-wide heatmaps revealed strong candidate cis-eQTLs on Chr 4 at Tyrp1 and on Chr 9 at Myo5a in all mice. In the wildtype data set, Tyrp1 was an upstream regulator of six pigmentation and two mesenchyme genes. In addition, five genes, including Tyrp1, were at least partially regulated by Myo5a. Analyses of the strains harboring the mutant gene revealed significant loss of correlation to traditional genes and gain of correlation to genes with little or no functional relationship. CONCLUSIONS: These findings indicate that the Tyrp1(b) mutation modifies the pathways and gene networks in which Tyrp1 functions. Our results also indicate direct and indirect regulatory control of Tyrp1 and other pigmentation and mesenchymal genes by Myo5a. Lastly, we find that the mutations reduce the ability of Tyrp1 to regulate expression of other genes that participate in pigmentation metabolism.

The effects of a CD81 null mutation on retinal pigment epithelium in mice.

The present study examines the effects of Cd81-null mutation on the development of the retinal pigment epithelium (RPE), specifically cell size and number of cells with multiple nuclei. The outlines of RPE in retinal flat mounts were stained with rhodamine-labeled phalloidin and RPE nuclei with Hoechst stain. The RPE layer was sampled to define the number of cells, the size of the RPE cells and the number of nuclei within the cells. The Cd81-null mutation caused an increase in the number of cells within the RPE layer. The cells were smaller than those in the wild type mice. Furthermore there was an increase in the number of mono-nucleated cells. In the posterior portion of the eye there was a significant increase in the number of multi-nucleated cells. The data indicate that CD81 plays a significant role in the final stages of RPE development, controlling cell number and overall developmental pattern.