Download datasets and supplementary data files |
---|
Summary
This is a subtractive dataset. The Normal retina dataset was subtracted from the ONC data set probe by probe to create a data set of the changes occurring following ONC. This data set can be used to define gene changes following ONC. It is not compatible with most of the bioinformatic tools available on GeneNetwork.
HEI Retina Illumina V6.2 (April 2010) RankInv was normalized and scaled by William E. Orr and uploaded by Arthur Centeno and Xiaodong Zhou on April 7, 2010. This data set consists of either 69 BXD strains (Normal data set) or 75 BXD strains (Full data set), C57BL/6J, DBA/2J, both reciprocal F1s, and BALB/cByJ. A total of either 74 strains (Normal data set) or 80 strains (Full data set) were quantified.
COMMENT on FULL versus NORMAL data sets: For many general uses there is no significant difference between FULL and NORMAL data sets. However, the FULL data set includes strains with high endogenous Gfap mRNA expression, indicative of reactive gliosis. For that reason, and to compare to OPTIC NERVE CRUSH (ONC), we removed data from six strains to make the NORMAL data set.
The NORMAL data set exludes data from BXD24, BXD32, BXD49, BXD70, BXD83, and BXD89. BXD24 has known retinal degeneration and is now known officially as BXD24/TyJ-Cep290/J, JAX Stock number 000031. BXD32 has mild retinal degeneration. The NORMAL data set does include BXD24a, now also known as BXD24/TyJ (JAX Stock number 005243).
The data are now open and available for analysis.
Please cite: Freeman NE, Templeton JP, Orr WE, Lu L, Williams RW, Geisert EE (2011) Genetic networks in the mouse retina: Growth Associated Protein 43 and Phosphate Tensin Homology network. Molecular Vision 17:1355-1372. Full Text PDF or HTML
This is rank invariant expression data that has been normalized using what we call a 2z+8 scale, but without special correction for batch effects. The data for each strains were computed as the mean of four samples per strain. Expression values on a log2 scale range from 6.25 to 18.08 (11.83 units), a nominal range of approximately 3600-fold. After taking the log2 of the original non-logged expression estimates, we convert data within an array to a z score. We then multiply the z score by 2. Finally, we add 8 units to ensure that no values are negative. The result is a scale with a mean of 8 units and a standard deviation of 2 units. A two-fold difference in expression is equivalent roughly to 1 unit on this scale.
The lowest level of expression is 6.254 for ILMN_2747167 (Arhgap11a) from HEI Retina Illumina V6.2 (April 2010) RankInv **. Lowest single data about 5.842.
The highest level of expression is 18.077 for ILMN_2516699 (Ubb). Highest single value is about 18.934.
Other Related Publications
- Geisert EE, Lu L, Freeman-Anderson NE, Templeton JP, Nassr M, Wang X, Gu W, Jiao Y, Williams RW.:Gene expression in the mouse eye: an online resource for genetics using 103 strains of mice. Molecular Vision 2009 Aug 31;15:1730-63, (Link)
- Geisert EE, Jr., Williams RW: The Mouse Eye Transcriptome: Cellular Signatures, Molecular Networks, and Candidate Genes for Human Disease. In Eye, Retina, and Visual System of the Mouse. Edited by Chalupa LM, Williams RW. Cambridge: The MIT Press; 2008:659-674
- Peirce JL, Lu L, Gu J, Silver LM, Williams RW: A new set of BXD recombinant inbred lines from advanced intercross populations in mice. BMC Genet 2004, 5:7. (Link)
- Templeton JP, Nassr M, Vazquez-Chona F, Freeman-Anderson NE, Orr WE, Williams RW, Geisert EE: Differential response of C57BL/6J mouse and DBA/2J mouse to optic nerve crush. BMC Neurosci. 2009, July 30;10:90.(Link)
Other Data Sets Users of these mouse retina data may also find the following complementary resources useful:
- NEIBank collection of ESTs and SAGE data.
- RetNet: the Retinal Information Network--tables of genes and loci causing inherited retinal diseases
- Mouse Retina SAGE Library from the Cepko laboratory. This site provides extensive developmental data from as early as embryonic day E12.5.
- Digital reference of ophthalmology from Columbia provides high quality photographs of human ocular diseases, case studies, and short explanations. This reference does not have a molecular focus.
- Mouse Retinal Developmental Gene Expression data sets from the Friedlander laboratory. This site provides extensive developmental data using the Affymetrix U74 v 2 array (predecessor of the M430).
- Data sets on differential gene expression in anatomical compartments of the human eye from Pat Brown's lab. View expression signatures for different ocular tissues using the geneXplorer 2.0.
Experiment design
Expression profiling by array
We used pooled RNA samples of retinas, usually two independent pools--two male, two female pool--for most lines of mice.
All normalization was performed by William E. Orr in the HEI Vision Core Facility
- Computed the log base 2 of each raw signal value
- Calculated the mean and standard Deviation of each Mouse WG-6 v2.0 array
- Normalized each array using the formula, 2 (z-score of log2 [intensity]) The result is to produce arrays that have a mean of 8, a variance of 4, and a standard deviation of 2. The advantage is that a two-fold difference in expression level corresponds approximately to a 1 unit difference.
- computed the mean of the values for the set of microarrays for each strain. Technical replicates were averaged before computing the mean for independent biological samples.
About cases
Almost all animals are young adults between 60 and 90 days of age (Table 1, minimum age is 48 and maximum age is 118 days). We measured expression in conventional inbred strains, BXD recombinant inbred (RI) strains, and reciprocal F1s between C57BL/6J and DBA/2J.
BXD strains:
- The first 32 of these strains are from the Taylor series of BXD strains generated at the Jackson Laboratory by Benjamin A. Taylor. BXD1 through BXD32 were started in the late 1970s, whereas BXD33 through 42 were started in the 1990s.
- In 2004, BXD24/TyJ developed a spontaneous mutation, rd16 which resulted in retinal degeneration and was renamed BXD24b/TyJ (BXD24 in this database). The strain, BXD24a, was cryo-recovered in 2004 from 1988 embryo stocks (F80) and does not exhibit retinal degeneration. In 2009, BXD24b was renamed BXD24/TyJ-Cep290rd16/J by JAX Labs to reflect the discovery of the genetic basis of the mutation. At the same time BXD24a was then referred to just as BXD24/TyJ by Jax Labs, but still called BXD24a in this dataset.
- The other 36 BXD strains (BXD43 and higher) were bred by Lu Lu, Jeremy Peirce, Lee M. Silver, and Robert W. Williams starting in 1997 using B6D2 generation 10 advanced intercross progeny. This modified breeding protocol doubles the number of recombinations per BXD strain and improves mapping resolution (Peirce et al. 2004). All of the Taylor series of BXD strains and many of the new BXD strains are available from the Jackson Laboratory. All of the new BXD strains (BXD43 and higher) are also available directly from Lu Lu and colleagues at the University of Tennessee Health Science Center in Memphis, TN, USA. BXD24/TyJ is now known as BXD24b/TyJ and has nearly complete retinal degeneration. BXD24a/TyJ, a 1988 F80 inbred stock that has been rederived from cryogenic storage, does not have retinal degeneration (stock number 005243) and is an ideal coisogenic control, but is not included in the HEI data set.
What Makes the G2 HEI Retina Database different from the HEI Retina Database Examination of Gfap expression across all of the strains in the HEI Retinal Dataset, reveals that some strains express very high levels of Gfap relative to others. For example, BXD24 expresses Gfap at a 9-fold higher level, than BXD22. It has been established that BXD24 acquired a mutation in Cep290 that results in early onset photoreceptor degeneration (Chang et al., 2006). This degeneration results in reactive gliosis throughout the retina. In addition to BXD24, other BXD strains expressed very high levels of Gfap including: BXD32, BXD49, BXD70, BXD83 and BXD89. For the G2 dataset all of these strains with potential reactive gliosis were removed from the dataset.
About tissue
Tissue preparation protocol. Animal were killed by rapid cervical dislocation. Retinas were removed immediately and placed in RNAlater at room temperature. Two retinas from one mouse were stored in a single tube.
Each array was hybridized with a pool of cRNA from 2 retinas (1 mouse). Natalie Freeman-Anderson extracted RNA at UTHSC.
Dissecting and preparing eyes for RNA extraction
Retinas for RNA extraction were placed in RNA STAT-60 (Tel-Test Inc.) and processed per manufacturer’s instructions (in brief form below). Total RNA was extracted with RNA STAT-60 (Tel-Test Inc.) according to the manufacturer's instructions. Briefly we:
- Homogenize tissue samples in the RNA STAT-60 (1 ml/50 to 100 mg tissue via syringe)
- Allow the homogenate to stand for 5-10 min at room temperature
- Add 0.2 ml of chloroform per 1 ml RNA STAT-60
- Mix the sample vigorously for 15 sec and let the sample incubate at room temperature for 5-10 min
- Centrifuge at 12,000 g for 1 hr at 4°C
- Transfer the aqueous phase to a clean centrifuge tube
- Add 0.5 ml of isopropanol per 1 ml RNA STAT-60
- Vortex and incubate the sample at -20°C for 1 hr or overnight
- Centrifuge at 12,000 g for 1 hr
- Remove the supernatant and wash the RNA pellet with 75% ethanol
- Remove ethanol, let air dry (5-10 min)
- Dissolve the pellet in 50 μl of nuclease free water.
About platform
Illumina MouseWG-6 v2.0 arrays: The Illumina Sentrix Mouse-6 BeadChip uses 50-nucleotide probes to interrogate approximately 46,000 sequences from the mouse transcriptome. For each array, the RNA was pooled from two retinas.
About data processing
Values of all 45,281 probe sets in this data set range from a low of 6.25 (Rho GTPase activating protein 11A, Arhgap11a, probe ID ILMN_2747167) to a high of 18.08 (Ubiquitin B, Ubb, probe ID ILMN_2516699). This corresponds to 11.83 units or a 1 to 3641 dynamic range of expression (2^11.83). We normalized raw signal values using Beadstudio’s rank invariant normalization algorithm. BXD62 was the strain used as the control group
Sample Processing: Drs. Natalie E. Freeman-Anderson and Justin P. Templeton extracted the retinas from the mice and Dr. Natalie Freeman-Anderson processed all samples in the HEI Vision Core Facility. The tissue was homogenized and extracted according to the RNA-Stat-60 protocol as described by the manufacturer (Tel-Test, Friendswood, TX) listed above. The quality and purity of RNA was assessed using an Agilent Bioanalyzer 2100 system. The RNA from each sample was processed with the Illumina TotalPrep RNA Amplification Kit (Ambion, Austin, TX) to produce labeled cRNA. The cRNA for each sample was then hybridized to an Illumina Sentrix® Mouse-6-V2 BeadChip (Illumina, San Diego, CA)
Quality control analysis of the raw image data was performed using the Illumina BeadStudio software. MIAME standards were used for all microarray data. Rank invariant normalization with BeadStudio software was used to calculate the data. Once this data was collected, the data was globally normalized across all samples using the formula 2 (z-score of log2 [intensity]) + 8.
Replication, sex, and sample balance: Our goal was to obtain data for independent biological sample pools from both sexes for most lines of mice. The four batches of arrays included in this final data set, collectively represent a reasonably well-balanced sample of males and females, in general without within-strain-by-sex replication.
Table 1: HEI Retina case IDs, including sample tube ID, strain, age, sex, and source of mice
|
Contributors
Eldon E. Geisert, Lu Lu, Natalie E. Freeman-Anderson, Justin P. Templeton, Robert W. Williams
Acknowledgment
The HEI Retinal Database is supported by National Eye Institute Grants:
- R01EY017841 (Dr. Eldon E. Geisert, PI)
- P030EY13080 (NEI Vision Core Grant), and
- A Unrestricted Grand from Research to Prevent Blindness (Dr. Barrett Haik, PI)
Notes
This data set is available as a bulk download in several formats. The data are available as either strain means or the individual arrays. Due to the involved normalization procedures required to correct for batch effects we strongly recommend not using the raw CEL files without special statistical procedures.