Datgan, a reusable software system for facile interrogation and visualization of complex transcription profiling data
© Howell et al; licensee BioMed Central Ltd. 2011
Received: 9 June 2011
Accepted: 24 August 2011
Published: 24 August 2011
We introduce Glaucoma Discovery Platform (GDP), an online environment for facile visualization and interrogation of complex transcription profiling datasets for glaucoma. We also report the availability of Datgan, the suite of scripts that was developed to construct GDP. This reusable software system complements existing repositories such as NCBI GEO or EBI ArrayExpress as it allows the construction of searchable databases to maximize understanding of user-selected transcription profiling datasets.
Datgan scripts were used to construct both the underlying data tables and the web interface that form GDP. GDP is populated using data from a mouse model of glaucoma. The data was generated using the DBA/2J strain, a widely used mouse model of glaucoma. The DBA/2J-Gpnmb + strain provided a genetically matched control strain that does not develop glaucoma. We separately assessed both the retina and the optic nerve head, important tissues in glaucoma. We used hierarchical clustering to identify early molecular stages of glaucoma that could not be identified using morphological assessment of disease. GDP has two components. First, an interactive search and retrieve component provides the ability to assess gene(s) of interest in all identified stages of disease in both the retina and optic nerve head. The output is returned in graphical and tabular format with statistically significant differences highlighted for easy visual analysis. Second, a bulk download component allows lists of differentially expressed genes to be retrieved as a series of files compatible with Excel. To facilitate access to additional information available for genes of interest, GDP is linked to selected external resources including Mouse Genome Informatics and Online Medelian Inheritance in Man (OMIM).
Datgan-constructed databases allow user-friendly access to datasets that involve temporally ordered stages of disease or developmental stages. Datgan and GDP are available from http://glaucomadb.jax.org/glaucoma.
Transcription profiling is a powerful tool for understanding biological process and the role they play in the pathogenesis of disease. For complex diseases, it is necessary to assess many different samples, resulting in very large amounts of data that are cumbersome to analyze and understand. Specific analyses often require significant computing power, time and analytical expertise. These needs hinder detailed interrogation of deposited datasets by many members of the scientific community. Transcription profiling datasets are deposited in central databases such as Gene Expression Omnibus (GEO) at the National Center for Biotechnology Information (NCBI)  and ArrayExpress at the European Bioinformatics Institute (EBI) . GEO and ArrayExpress are designed to efficiently store large data sets and to provide mechanisms that allow the scientific community to query, locate, review and download experiments of interest. Although useful, we have found that these platforms are not ideal for optimal interrogation of large and complex datasets. We became aware of the need for a more optimized and facile environment for analyzing complex datasets when studying glaucoma, a complex and asynchronous neurodegenerative disease. Our datasets consists of multiple, temporally ordered stages of disease and existing datasets do not allow simultaneous searching and retrieving of differentially expressed genes in multiple stages of disease compared to a no-disease control.
To meet this need, we developed Glaucoma Discovery Platform (GDP), an online environment to visualize and interrogate large and complex expression datasets in a user-friendly manner. GDP allows simultaneous querying of multiple genes, viewing results across multiple datasets and assessing expression differences for multiple probe sets for each gene. To our knowledge, no other resource is available that provides this combination of user-friendly functionality provided by GDP. We deployed GDP to derive maximum benefit from an extensive transcription profiling study of glaucoma. The study focused on identifying early molecular stages that occur prior to significant damage (described in more detail below and ). More than 100 different samples from at least 50 mice were profiled and 70 pairwise group comparisons were made. GDP utilizes a web-based user interface to access the underlying gene expression profiling data, which is organized on a MySQL database server. GDP has greatly facilitated our understanding of these glaucoma data. Its user-friendly interface provides easy and instant data interrogation without a need for specialized knowledge or training. It thereby allows general access to our glaucoma datasets and is freely available for the benefit of the wider scientific community.
We also report the availability of Datgan (Welsh, meaning 'to express'), the suite of scripts that was used to construct both the underlying data tables and the web interface that form GDP. Datgan was written in the python language. It is packaged as a reusable software system for constructing and populating discovery platforms to visualize user-determined sets of transcription profiling data. Datgan is organized to allow biologists with some experience running scripts (or a systems administrator) to establish their own personalized discovery platforms and can be readily adapted to incorporate transcription profiling data generated using RNA-seq.
Construction and content
Populating GDP with profiling data for glaucoma
Glaucoma is a complex, neurodegenerative disorder affecting 70 million people worldwide and is associated with the death of retinal ganglion cells (RGCs) and the associated degeneration of the optic nerve . DBA/2J is a widely used mouse model of glaucoma that shows hallmarks of human glaucoma including age-related IOP elevation, optic nerve excavation and regional patterns of RGC loss [5–10]. DBA/2J mice develop glaucoma as a result of a disease of the iris that leads to an elevation in IOP. The disease of the iris is caused by mutations in two genes, Gpnmb R150X and Tyrp1 b [5, 6]. DBA/2J-Gpnmb + mice have a functioning Gpnmb gene and serve as a genetically matched control strain that does not develop glaucoma . An important insult occurs to RGC axons at the optic nerve head in DBA/2J glaucoma . However, other compartments of the RGC also are likely to undergo early changes in glaucoma such as the RGC soma  and synapses . The mechanisms involved in these early changes are not well understood.
Building GDP using Datgan scripts
Visualizing and interrogating profiling data with GDP
Home page/new search (Figure 3a)
The home page provides an overview of the database and permanent headers for convenient links to the bulk download page and the detailed users quick guide. The bulk downloads page allows all differentially expressed genes and associated information for individual comparisons to be downloaded in excel-friendly format. The quick search feature allows a single gene to be searched in all datasets. Information on downloading Datgan is also provided.
The results summary page (Figure 3b)
Results are returned in both graphical and tabular format. Any gene names that were not recognized from the search page are listed and links to MGI are given (for clarification of official gene symbols). A 'summary of results' table is shown indicating the number of genes found for each dataset/tissue selected on the search page. Below the summary of results are separate tables and associated graphs for each dataset/tissue searched. Each table contains the gene names searched, the official gene symbol (or description), the representative probe set, and the fold change with associated q value for each stage of glaucoma, compared to the chosen reference group. For ease of identification, significantly differentially expressed values are shown in white on a red background. Genes are ranked in the table based on the chosen option in the detailed search tool (q value or gene symbol). The graph provides a visual display of the fold changes across all stages of disease for the genes of interest, up to a maximum of ten. For searches containing more than 10 genes, the first 10 genes in the table are displayed in the graph. Check boxes are available to the left of each gene to allow the user to select up to 10 genes of choice to view in the graph. Results from the table can be exported as comma separated values (csv, Excel-friendly format). Links to useful publicly available databases are also provided (see External resources below).
Probe set details page (Figure 3c)
Each gene on the Affymetrix 430 v2 array is interrogated by a probe set of eleven 25 mer probes and a gene can have multiple probe sets . Each probe set corresponds to a particular region of a gene and finding that multiple probes sets for the same gene behave similarly can add confidence to a result. Alternatively, for those probe sets corresponding to alternatively spliced exons/transcripts, multiple probe sets can give insight into the behavior of splice variants . Unfortunately, some probe sets were designed to early versions of gene sequences (prior to accurate genome sequence). This results in some probes in a probe set not being identical to the gene sequence. These probe sets will not accurately reflect the level of transcript for these genes. For each probe set, a link is provided to the Ensembl website where mapping information for each probe set is provided. The values for the representative probe set for each gene is displayed in the results summary page. It is possible to view details for all probe sets for a given gene, to view probe sets for selected genes or for all genes in the table. The probe set detail page provides a table/graph for each gene.
Expression values page (Figure 3d)
Values provided in the summary pages are relative fold changes or log2 fold changes compared to a selected reference. Each stage of disease contains biological replicates, and the fold change for a given probe set (for a given gene) is calculated as the average of the normalized expression value of all replicates for a defined stage of disease, relative to the reference group. The normalized expression value reflects the relative abundance of the gene in the tissue interrogated (either retina or optic nerve head) with respect to all other genes. Those transcripts with the highest normalized expression value are most abundant. Conversely, those with the lowest normalized expression value are less abundant. An expression level of approximately 4 or less may be considered to represent a gene that is likely not expressed in the assessed tissue (i.e. is close to background levels).
The relative abundance of transcripts corresponding to particular probe sets is important information. For instance, a small fold change in some lowly expressed gene may have greater biological importance than a small fold change in some abundant gene. Knowing the variability in expression levels for different genes in individual eyes within groups also is important. Within any stage of disease (determined either morphologically or molecularly), transcripts within individual biological replicates may behave differently. Genes with low variability within groups may be better targets for intervention strategies than variable genes. The expression values corresponding to each probe set for individual samples (replicates) can be accessed from the probe details page. The expression levels are displayed as a histogram. The average (± 1 standard deviation) for the reference group is indicated.
Summary of external resources accessed directly from GDP
Home page address
Mouse Genome Informatics
Gene catalogues, links to additional resources including phenotype data, genome browsers and SNP databases
Gene catalogues, links to additional resources including genome browsers and SNP databases
Online Medelian Inheritance in Man
The Gene Ontology@MGI
Controlled vocabulary for functional attributes
Kyoto Encyclopedia of Genes and Genomes
Ensembl Genome Browser
Determine alignment of probe sets onto genome sequence
In this section, we describe the workflow for extracting the data relevant to members of the TNF superfamily. This serves both to reinforce key features of the database described above and as an example that could be followed to interrogate any gene (or group of genes) of interest.
Assessing members of the tumor necrosis factor (TNF) superfamily
TNF (formerly TNFα) has been suggested to be important in retinal ganglion cell loss during glaucoma [24, 25]. Therefore, we assessed all members of the TNF superfamily in both the molecular ONH dataset (Figure 2B, most sensitive dataset of early disease changes in the optic nerve head), and the molecular retina dataset (Figure 2C, most sensitive dataset of early disease changes in the retina).
Step 1: Performing detailed search (Figure 4)
First, we determined the best search option. In this case, a wild card search using 'TNF*' will retrieve all TNF superfamily members as no a priori knowledge about which genes are present in the TNF superfamily is required. Second, we selected the datasets and tissues of interest - molecular ONH dataset, tissue 'ONH' and molecular retina dataset, tissue 'retina'. Finally, we selected the reference group; in this case, D2-Gpnmb + control eyes. In this example, we did not limit our search based on fold change or q value, and chose the default options for the layout of the results page (display expression differences as fold change rather than log2 of fold change, and ordered genes based on significance of gene expression differences (q value) rather than gene symbol).
Step 2: Visualizing the results (Figure 5)
All 'Tnf'-related genes in the molecular optic nerve head dataset
Representative Probe Set ID
Num Probe Sets
FC for Stage 1
QV for Stage 1
FC for Stage 2
QV for Stage 2
FC for Stage 3
QV for Stage 3
FC for Stage 4
QV for Stage 4
FC for Stage 5
QV for Stage 5
Step 3: Interpreting the results (Figure 6)
Tnf shows only a modest increase in expression in stage 4 (1.2 fold, q value = 0.0413) in the molecular ONH dataset. The five most significant genes are Tnf, alpha-induced protein 8-like 2 (Tnfaip8l2), Tnf superfamily receptor 1a (Tnfrsf1a, formerly Tnfr1) and 1b (Tnfrsf1b, formerly Tnfr2), Tnf alpha-induced protein 2 (Tnfaip2) and Tnf superfamily member 9 (Tnfrsf9) (Figure 5C, Table 2). Tnfrsf9 has the highest fold change of all TNF-related genes in the molecular ONH dataset, with a 2.9 fold expression difference compared to D2-Gpnmb + controls in stage 3. Of direct relevance to DBA/2J glaucoma, Tnfrsf9 is involved in the proliferation of monocytes that are precursors of microglia and macrophages [26–28]. This information was retrieved from either OMIM or Entrez Gene. Microglia/macrophages have been shown to increase in the optic nerve head and retina early in glaucoma [3, 29].
The differing expression levels for the three Tnfrsf9 probe sets may result from detection of alternative splice forms by the different probe sets or from errors in the original probe set designs. The apparent differences in expression levels between Tnfrsf9 probe sets are likely due to probe sequence errors. Detailed data for all probes within each probe set is available in Ensembl and can be accessed directly from GDP. Only for probe set 1460469_at do all eleven probes match identically within the Tnfrsf9 gene (Figure 6c). Ensembl does not report an alignment for any probes in the 1421481_at probe set. For the remaining probe set, 1428034_a_at, three of the eleven probes have mismatches, and another three probes include intronic sequences that would not be included in the messenger RNA. This explains a lack of expression for this probe set as an average of the expression level of all probes is taken. Given this detailed probe information, it is clear that the expression level of Tnfrsf9 is only truly represented by the 1460469_at probe set.
In summary, using GDP, we have readily performed a detailed search of all TNF superfamily members and readily identified all genes that change significantly compared to controls in two out of three different datasets. In particular, the expression of Tnfrsf9 increases in the optic nerve head at very early stages of glaucoma and may be related to an accumulation of microglia. These results are of great interest to the glaucoma community and are immediately and freely available. They highlight major advantages of accessing these gene expression datasets through GDP, and the utility of Datgan for providing similar environments for many other datasets.
Here, we describe GDP, constructed using the reusable software system Datgan, for visualizing and interrogating complex gene expression profiling data. Datgan is freely available for download (via http://www.simonjohnlab.org) and can be used to construct platforms for viewing and interrogating any set of microarray or RNA-seq based transcription profiling datasets. Datgan will be particularly useful for complex diseases that have variable onset and progression. As exemplified by GDP, platforms developed using Datgan allow general access and improved understanding of complex datasets.
Use of GDP is proving very important for understanding our glaucoma datasets. For publication in most journals, all profiling datasets are required to be deposited in either GEO (NCBI) or ArrayExpress (EBI). All data and associated files for our study were deposited in NCBI GEO (Accession number GSE26299). Data deposition in these archives is necessary to ensure access to all the raw data for interested parties. However, these archives were not designed to allow detailed interrogation of multiple genes or biological pathways across multiple datasets. GDP includes this functionality and massively facilitates interrogation of our glaucoma datasets by us, and the scientific community, ultimately allowing the most benefit to be derived from the data. Other gene expression profiling studies relevant to glaucoma, which have been or will be carried out, also could be incorporated into GDP. Many investigators, irrespective of computational expertise, can perform detailed 'on the fly' searches quickly and easily.
Datgan is specifically designed to allow additional datasets to be easily incorporated after the initial online environment has been constructed. For glaucoma, other gene expression profiling studies have been carried out that could be included into GDP. These include studies that have profiled additional animal models of glaucoma [30–32]. In addition, comparing cell- and tissue- specific transcription profiling studies (e.g ) in non-diseased settings would enable hypotheses to be made about which cells are changing early in glaucoma. Finally, useful predictions could be made by comparing studies that have profiled other relevant diseases such as other neurodegenerative disorders or diseases affecting retinal ganglion cells.
Datgan is a powerful software package for developing user-friendly platforms to visualize transcription profiling data. Implemented as GDP, it is already proving an essential tool for interrogating the molecular pathogenesis of glaucoma.
Availability and Requirements
The database can be accessed via http://glaucomadb.jax.org/glaucoma. Web browsers: Tested on the major web browsers including the latest versions of Firefox, Safari and Chrome. It has also been tested on smart phones and tablet computers. Operating system(s): Tested on SUSE Linux Enterprise Server 11 & OpenSUSE 11.0, and should work on other versions of Linux and Mac OSX. Tools required to establish Datgan-derived databases include: Ruby 1.8.7, Rails 2.3.5, Python 2.6, Percona-Server-5.5.13 (should work with most MySQL), apache2, Mongrel Web Server 1.1.5. Source code details can be obtained from the Datgan tab at the top of the home page.
Acknowledgements and Funding
We acknowledge the Scientific Services at The Jackson Laboratory, particularly Tim Stearns of Computational Sciences who performed the computational comparisons of the datasets uploaded to GDP. This work was funded in part by American Health Assistance Foundation (GRH), Glaucoma Research Foundation (GRH), EY018606 (RTL), Research to Prevent Blindness Career Development Award (RTL), EY011721 (SWMJ), The Barbara and Joseph Cohen Foundation, and The Partridge Foundation. SWMJ is an Investigator of the Howard Hughes Medical Institute.
- Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Edgar R: NCBI GEO: mining tens of millions of expression profiles--database and tools update. Nucleic Acids Res. 2007, 35 (Database): D760-765. 10.1093/nar/gkl887.PubMedView ArticleGoogle Scholar
- Parkinson H, Kapushesky M, Shojatalab M, Abeygunawardena N, Coulson R, Farne A, Holloway E, Kolesnykov N, Lilja P, Lukk M, et al: ArrayExpress--a public database of microarray experiments and gene expression profiles. Nucleic Acids Res. 2007, 35 (Database): D747-750. 10.1093/nar/gkl995.PubMedView ArticleGoogle Scholar
- Howell GR, Macalinao DG, Sousa GS, Walden M, Soto I, Kneeland SL, Barbay JM, King BL, Marchant JK, Hibbs M, et al: Molecular clustering identifies complement and endothelin induction as early events in a mouse model of glaucoma. J.Clinical Investigation. 121 (4): 1429-1444.
- Quigley HA, Broman AT: The number of people with glaucoma worldwide in 2010 and 2020. Br J Ophthalmol. 2006, 90 (3): 262-267. 10.1136/bjo.2005.081224.PubMedPubMed CentralView ArticleGoogle Scholar
- Anderson MG, Smith RS, Hawes NL, Zabaleta A, Chang B, Wiggs JL, John SW: Mutations in genes encoding melanosomal proteins cause pigmentary glaucoma in DBA/2J mice. Nat Genet. 2002, 30 (1): 81-85. 10.1038/ng794.PubMedView ArticleGoogle Scholar
- Chang B, Smith RS, Hawes NL, Anderson MG, Zabaleta A, Savinova O, Roderick TH, Heckenlively JR, Davisson MT, John SW: Interacting loci cause severe iris atrophy and glaucoma in DBA/2J mice. Nat Genet. 1999, 21 (4): 405-409. 10.1038/7741.PubMedView ArticleGoogle Scholar
- Howell GR, Libby RT, Jakobs TC, Smith RS, Phalan FC, Barter JW, Barbay JM, Marchant JK, Mahesh N, Porciatti V, et al: Axons of retinal ganglion cells are insulted in the optic nerve early in DBA/2J glaucoma. J Cell Biol. 2007, 179 (7): 1523-1537. 10.1083/jcb.200706181.PubMedPubMed CentralView ArticleGoogle Scholar
- Jakobs TC, Libby RT, Ben Y, John SW, Masland RH: Retinal ganglion cell degeneration is topological but not cell type specific in DBA/2J mice. J Cell Biol. 2005, 171 (2): 313-325. 10.1083/jcb.200506099.PubMedPubMed CentralView ArticleGoogle Scholar
- Libby RT, Anderson MG, Pang IH, Robinson ZH, Savinova OV, Cosma IM, Snow A, Wilson LA, Smith RS, Clark AF, et al: Inherited glaucoma in DBA/2J mice: pertinent disease features for studying the neurodegeneration. Vis Neurosci. 2005, 22 (5): 637-648.PubMedView ArticleGoogle Scholar
- Schlamp CL, Li Y, Dietz JA, Janssen KT, Nickells RW: Progressive ganglion cell loss and optic nerve degeneration in DBA/2J mice is variable and asymmetric. BMC Neurosci. 2006, 7: 66-10.1186/1471-2202-7-66.PubMedPubMed CentralView ArticleGoogle Scholar
- Howell GR, Libby RT, Marchant JK, Wilson LA, Cosma IM, Smith RS, Anderson MG, John SW: Absence of glaucoma in DBA/2J mice homozygous for wild-type versions of Gpnmb and Tyrp1. BMC Genet. 2007, 8 (1): 45-PubMedPubMed CentralView ArticleGoogle Scholar
- Soto I, Oglesby E, Buckingham BP, Son JL, Roberson ED, Steele MR, Inman DM, Vetter ML, Horner PJ, Marsh-Armstrong N: Retinal ganglion cells downregulate gene expression and lose their axons within the optic nerve head in a mouse glaucoma model. J Neurosci. 2008, 28 (2): 548-561. 10.1523/JNEUROSCI.3714-07.2008.PubMedView ArticleGoogle Scholar
- Stevens B, Allen NJ, Vazquez LE, Howell GR, Christopherson KS, Nouri N, Micheva KD, Mehalow AK, Huberman AD, Stafford B, et al: The classical complement cascade mediates CNS synapse elimination. Cell. 2007, 131 (6): 1164-1178. 10.1016/j.cell.2007.10.036.PubMedView ArticleGoogle Scholar
- Cui X, Churchill GA: Statistical tests for differential expression in cDNA microarray experiments. Genome Biol. 2003, 4 (4): 210-10.1186/gb-2003-4-4-210.PubMedPubMed CentralView ArticleGoogle Scholar
- Libby RT, Li Y, Savinova OV, Barter J, Smith RS, Nickells RW, John SW: Susceptibility to neurodegeneration in a glaucoma is modified by Bax gene dosage. PLoS Genet. 2005, 1 (1): 17-26. 10.1371/journal.pgen.0010017.PubMedView ArticleGoogle Scholar
- Bult CJ, Eppig JT, Kadin JA, Richardson JE, Blake JA: The Mouse Genome Database (MGD): mouse biology and model systems. Nucleic Acids Res. 2008, 36 (Database): D724-728.PubMedView ArticleGoogle Scholar
- Blake JA, Bult CJ, Eppig JT, Kadin JA, Richardson JE: The Mouse Genome Database genotypes::phenotypes. Nucleic Acids Res. 2009, 37 (Database): D712-719. 10.1093/nar/gkn886.PubMedView ArticleGoogle Scholar
- Sayers EW, Barrett T, Benson DA, Bryant SH, Canese K, Chetvernin V, Church DM, DiCuccio M, Edgar R, Federhen S, et al: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2009, 37 (Database): D5-15. 10.1093/nar/gkn741.PubMedView ArticleGoogle Scholar
- Aoki KF, Kanehisa M: Using the KEGG database resource. Curr Protoc Bioinformatics. 2005, Chapter 1: Unit 1 12Google Scholar
- Churchill GA: Using ANOVA to analyze microarray data. Biotechniques. 2004, 37 (2): 173-175. 177PubMedGoogle Scholar
- Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31 (4): e15-10.1093/nar/gng015.PubMedPubMed CentralView ArticleGoogle Scholar
- Stalteri MA, Harrison AP: Interpretation of multiple probe sets mapping to the same gene in Affymetrix GeneChips. BMC Bioinformatics. 2007, 8: 13-10.1186/1471-2105-8-13.PubMedPubMed CentralView ArticleGoogle Scholar
- Shaw DR: Searching the Mouse Genome Informatics (MGI) resources for information on mouse biology from genotype to phenotype. Curr Protoc Bioinformatics. 2009, Chapter 1: Unit1 7Google Scholar
- Nakazawa T, Nakazawa C, Matsubara A, Noda K, Hisatomi T, She H, Michaud N, Hafezi-Moghadam A, Miller JW, Benowitz LI: Tumor necrosis factor-alpha mediates oligodendrocyte death and delayed retinal ganglion cell loss in a mouse model of glaucoma. J Neurosci. 2006, 26 (49): 12633-12641. 10.1523/JNEUROSCI.2801-06.2006.PubMedView ArticleGoogle Scholar
- Tezel G, Li LY, Patil RV, Wax MB: TNF-alpha and TNF-alpha receptor-1 in the retina of normal and glaucomatous eyes. Invest Ophthalmol Vis Sci. 2001, 42 (8): 1787-1794.PubMedGoogle Scholar
- Jiang D, Chen Y, Schwarz H: CD137 induces proliferation of murine hematopoietic progenitor cells and differentiation to macrophages. J Immunol. 2008, 181 (6): 3923-3932.PubMedView ArticleGoogle Scholar
- Langstein J, Michel J, Fritsche J, Kreutz M, Andreesen R, Schwarz H: CD137 (ILA/4-1BB), a member of the TNF receptor family, induces monocyte activation via bidirectional signaling. J Immunol. 1998, 160 (5): 2488-2494.PubMedGoogle Scholar
- Langstein J, Michel J, Schwarz H: CD137 induces proliferation and endomitosis in monocytes. Blood. 1999, 94 (9): 3161-3168.PubMedGoogle Scholar
- Bosco A, Steele MR, Vetter ML: Early microglia activation in a mouse model of chronic glaucoma. J Comp Neurol. 519 (4): 599-620.
- Steele MR, Inman DM, Calkins DJ, Horner PJ, Vetter ML: Microarray analysis of retinal gene expression in the DBA/2J model of glaucoma. Invest Ophthalmol Vis Sci. 2006, 47 (3): 977-985. 10.1167/iovs.05-0865.PubMedView ArticleGoogle Scholar
- Johnson EC, Jia L, Cepurna WO, Doser TA, Morrison JC: Global changes in optic nerve head gene expression after exposure to elevated intraocular pressure in a rat glaucoma model. Invest Ophthalmol Vis Sci. 2007, 48 (7): 3161-3177. 10.1167/iovs.06-1282.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang Z, Quigley HA, Pease ME, Yang Y, Qian J, Valenta D, Zack DJ: Changes in Gene Expression in Experimental Glaucoma and Optic Nerve Transection: The Equilibrium between Protective and Detrimental Mechanisms. Invest Ophthalmol Vis Sci. 2007, 48 (12): 5539-5548. 10.1167/iovs.07-0542.PubMedView ArticleGoogle Scholar
- Cahoy JD, Emery B, Kaushal A, Foo LC, Zamanian JL, Christopherson KS, Xing Y, Lubischer JL, Krieg PA, Krupenko SA, et al: A transcriptome database for astrocytes, neurons, and oligodendrocytes: a new resource for understanding brain development and function. J Neurosci. 2008, 28 (1): 264-278. 10.1523/JNEUROSCI.4178-07.2008.PubMedView ArticleGoogle Scholar
- Anderson MG, Libby RT, Gould DB, Smith RS, John SW: High-dose radiation with bone marrow transfer prevents neurodegeneration in an inherited glaucoma. Proc Natl Acad Sci USA. 2005, 102 (12): 4566-4571. 10.1073/pnas.0407357102.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.