Eye M430v2 (November05) PDNN modify this page

Accession number: GN94

    Summary:

CAUTION: DO NOT USE THIS PDNN TRANSFORM of the HEIMED EYE Database. USE RMA INSTEAD. This April 2005 data freeze provides estimates of mRNA expression in adult eye from 50 lines of mice including C57BL/6J, DBA/2J, their F1 hybrids, and 47 BXD recombinant inbred strains. Data were generated at UTHSC. Samples were hybridized in small pools (n = 3) to Affymetrix M430 2.0 arrays. This particular data set was processed using the PDNN method of Zhang. To simplify comparison among transforms, PDNN values of each array were adjusted to an average of 8 units and a standard deviation of 2 units.

    About the cases used to generate this set of data:

We have used a set of BXD recombinant inbred strains generated by crossing C57BL/6J (B6 or B) with DBA/2J (D2 or D). The BXDs are particularly useful for systems genetics because both parental strains have been sequenced (8x coverage of B6 and 1.5x coverage of D). Physical maps in WebQTL incorporate approximately 2 million B vs D SNPs from Celera Genomics and from the Perlegen-NIEHS sequencing effort. BXD2 through BXD32 were bred by Benjamin A. Taylor starting in the late 1970s. BXD33 through 42 were bred by Taylor in the 1990s. These strains are available from The Jackson Laboratory. BXD43 through BXD99 were bred by Lu Lu, Jeremy Peirce, Lee M. Silver, and Robert W. Williams in the late 1990s and early 2000s using advanced intercross progeny (Peirce et al. 2004). Many of the 50 new BXD strains are available from Lu Lu and colleagues.

All stock was obtained originally from The Jackson Laboratory between 1998 and 2003. Most BXD animals were born and housed at the University of Tennessee Health Science Center. Some cases were bred at the University of Memphis (Douglas Matthews) or the University of Alabama (John Mountz and Hui-Chen Hsu).

    About the tissue used to generate this set of data:

The INIA M430 brain Database (April05) consists of 105 Affymetrix 430A and 430B microarray pairs. Each pair was hybridized in sequence (A array first, B array second) with a pool of brain tissue (forebrain minus olfactory bulb, plus the entire midbrain) taken from three adult animals of closely matched age and the same sex. RNA was extracted at UTHSC by Lu Lu, Zhiping Jia, and Hongtao Zhai at UTHSC.

Tissue preparation protocol. Animal were killed by rapid cervical dislocation. Eyes were removed immediately and placed in RNAlater at room temperature. Usually six eyes from animals with a common sex, age, and strain were stored in a single tube. The body was sprayed lightly with 70% ethanol to wet the hair. the following standard approach was used to extract the brain:

  1. Using small surgical scissors make an incision under the skin on the dorsal side of the neck. Cut the skin overlying the skull close to the midsagittal plane towards the nose. Pull and reflect the skin to expose the entire dorsal skull.
  2. Slip the points of the scissors through into the cisterna magna just caudal to the cerebellum and gently enlarge this opening until is it possible to cut through the skull overlying the cerebellum.
  3. Cut rostrally through the skull along the midsagittal line almost all the way to the nasal opening, taking care not to damage the dorsal surface of the brain.
  4. Approximately midway along this incision, make a lateral cut. Repeat along the incision and peel back the resulting strips of skull.
  5. Using small forceps, free the olfactory bulbs rostrally and ventrally, taking care to retain their connection to the rest of the forebrain.
  6. Gently lift the brain from the base the skull starting from the olfactory bulbs, pulling the brain toward a nearly vertical position. Cut the optic and trigeminal nerves. Separate the brain from the spinal cord about 2 mm distal to the medulla.
  7. Spread the hemispheres of the forebrain gently with forceps and then cut from dorsal to ventral using a straight scalpel, separating the hemispheres from each other (but not from the cerebellum). Take care to retain both paraflocculi.
At this point the protocol divides. If tissue is to be saved for RNA extraction at a later time, the whole brain is placed directly in RNAlater (Ambion, Inc.) and treated per the manufacturer’s directions. Step 7 is still very important because RNAlater may not fully penetrate the forebrain if the lobes are not separated. If tissue is to be used for immediate RNA extraction, one lobe of the forebrain is removed for processing and the rest of the brain is stored in RNAlater.

Dissecting and preparing forebrain and midbrain for RNA extraction

  1. Remove the left or right hemisphere of the forebrain and midbrain (referred to here as the forebrain for simplicity), either fresh or preserved in RNAlater by cutting from the caudal border of the inferior colliculus on the dorsal side and extending the cut ventrally to the basis pedunculi and the pons (cut just rostral of the pons) on the ventral side. See steps 7 and 8 here
  2. Place tissue for RNA extraction in RNA STAT-60 (Tel-Test Inc.) and process per manufacturer’s instructions (in brief form below).
  3. Store RNA in 75% ethanol at –80 deg. C until use.

Total RNA was extracted with RNA STAT-60 (Tel-Test) according to the manufacturer’s instructions. Briefly we:

  1. homogenize tissue samples in the RNA STAT-60 (1 ml/50 to 100 mg tissue)
  2. allowed the homogenate to stand for 5 min at room temperature
  3. added 0.2 ml of chloroform per 1 ml RNA STAT-60
  4. shook the sample vigorously for 15 sec and let the sample sit at room temperature for 3 min
  5. centrifuged at 12,000 G for 15 min
  6. transfered the aqueous phase to a fresh tube
  7. added 0.5 ml of isopropanol per 1 ml RNA STAT-60
  8. vortexed and allowed sample to stand at room temperature for 5-10 min
  9. centrifugeed at 12,000 G for 10-15 min
  10. removed the supernatant and washed the RNA pellet with 75% ethanol
  11. stored the pellet in 75% ethanol at -80 deg C until use

Sample Processing. Samples were processed in the INIA Bioanalytical Core at the W. Harry Feinstone Center of Excellence, The University of Memphis, lead by Thomas R. Sutter. All processing steps were performed by Shirlean Goodwin. In brief, samples were quality control checked for RNA purity using 260/280 ratios (samples had to be greater than 1.8, but the majority were 1.9 or higher). RNA integrity was assessed using the Agilent Bioanalyzer 2100. We required an RNA integrity number (RIN) of greater than 8, based on the intensity ratio and amplitude of 18S and 28S rRNA signals. The standard Eberwine T7 polymerase method was used to catalyze the synthesis of cDNA template from polyA-tailed RNA using Superscript II RT (Invitrogen Inc.). The Enzo LIfe Sciences, Inc., BioArray High Yield RNA Transcript Labeling Kit (T7, Part No. 42655) was used to synthesize labeled cRNA. At this point the cRNA was evaluated again using both the 260/280 ratio (values of 2.0 or above are acceptable) and the Bioanalyzer output (a dark cRNA smear on the 2100 output centered roughly between 600 and 2000 nt is required). Those samples that passed both QC steps (10% usually fail) were then sheared using a fragment buffer included in the Affymetrix GeneChip Sample Cleanup Module (Part No.900371). After fragmentation, samples were either stored at -80 deg. centrigrade until use or were immediately injected onto the array.

Replication and Sample Balance: Our goal is to obtain data for independent biological sample pools from at least one of sample from each sex for all BXD strains. We have not yet achieved this goal. Ten of 45 strains are still represented by single sex samples: BXD2 (F), BXD8 (F), BXD15 (F), BXD18 (F), BXD25 (F), BXD29 (F), BXD33 (M), BXD45 (F), BXD77 (M), and BXD90 (M). Eleven strains are represented by three independent samples with the following breakdown by sex: C57BL/6J (1F 2M), DBA/2J (2F 2M), B6D2F1 (2F 2M) + D2B6F1 (1F 1M), BXD6 (2F 1M), BXD13 (2F 1M), BXD14 (1F 2M), BXD28 (2F 1M), BXD34 (1F 2M), BXD36 (1F 2M), BXD38 (1F 2M), BXD42 (1F 2M).

Batch Structure: Before running the first batch of 30 pairs of array (dated Jan04), we ran four test samples (Nov03). The main batch of 30 includes the four test samples (four technical replicates). The Nov03 data was combined with the Jan04 data and was treated as a single batch that consists of one male and one female pool from C57BL/6J, DBA/2J, the B6D2F1 hybrid, 11 female BXD samples, and 11 male BXD samples. The second large batch was run February 2005 (Feb05) and consists of 71 pairs of arrays. Batch effects were corrected at the individual probe level as described below.

The table below summarizes information on tubeID, group_type, strain, age, sex, original CEL filename, and source of mice.

Idtube IDgroup typeStrain age sex original CEL filename PDNN 2Z outlier RMA 2Z outlier scale factor background average present absent marginal AFFX-b-ActinMur(3'/5') AFFX-GapdhMur(3'/5') source
1R2607E1GDP BXDC57BL/6J67FR2605E.CEL0.0050.0092.428115.120.5860.40.0141.310.76UTM RW
2R0872E2GDP BXDC57BL/6J66MR0872E.CEL0.0130.0123.12888.580.5890.3960.0151.30.79UTM RW
3R2572E1GDP BXDDBA/2J65MR2572E.CEL0.0410.0512.40679.070.5550.4290.0161.370.79UTM RW
4R2601E1GDP BXDB6D2F173FR2601E.CEL0.0030.0042.54591.960.5890.3960.0151.440.78UTM RW
5R2602E1GDP BXDB6D2F173MR2602E.CEL0.0010.0042.59984.440.5970.3880.0151.370.78UTM RW
6R2600E1GDP BXDD2B6F172FR2600E.CEL0.0030.0082.4794.750.5810.4020.0171.410.78UTM RW
7R2604E1GDP BXDD2B6F169MR2604E.CEL0.0030.0072.65789.630.5940.3920.0151.280.79UTM RW
8R2597E1BXDBXD261MR2597E.CEL0.0030.0072.37493.560.6030.3830.0151.340.77Glenn
9R2591E1BXDBXD560FR2591E.CEL0.0510.0091.7136.480.5850.40.0151.330.78Glenn
10R2570E1BXDBXD665FR2570E.CEL0.0020.0061.98786.730.5850.40.0151.460.76UTM RW
11R2538E1BXDBXD877FR2538E.CEL0.0370.0281.905101.980.6120.3730.0151.520.79UTM RW
12R2569E1BXDBXD967MR2569E.CEL0.0140.0271.75387.360.5510.4340.0152.823.14UTM RW
13R2581E1BXDBXD1165FR2581E.CEL0.0060.0121.94188.550.6210.3640.0161.550.81UTM RW
14R2543E1BXDBXD1263MR2543E.CEL0.0360.0071.605117.690.5860.3990.0161.430.77UTM RW
15R2586E1BXDBXD1360FR2586E.CEL0.0200.0352.00673.610.5640.420.0162.853.81Glenn
16R2557E1BXDBXD1460FR2557E.CEL0.0140.0171.8398.760.6250.3610.0141.310.78Glenn
17R2567E1BXDBXD1660MR2567E.CEL0.0160.0252.23982.350.5670.4160.0171.370.75Glenn
18R2559E1BXDBXD1859MR2559E.CEL0.0350.0061.654103.680.6080.3770.0151.270.78Glenn
19R2560E1BXDBXD1960FR2560E.CEL0.0260.0071.79298.330.6090.3750.0161.350.8Glenn
20R2584E1BXDBXD2059FR2584E.CEL0.0030.0072.0783.820.5930.3910.0161.40.76Glenn
21R2541E2BXDBXD2161MR2541E2.CEL0.0490.0362.625125.080.560.4240.0151.290.78UTM RW
22R2553E1BXDBXD2258FR2553E.CEL0.0030.0051.952111.30.5990.3850.0151.280.76Glenn
23R2558E1BXDBXD2360FR2558E-2.CEL0.0130.0151.908114.530.5990.3880.0141.20.82Glenn
24R2589E2BXDBXD2459MR2589E2.CEL0.0980.0982.606112.190.5750.4090.0161.240.8Glenn
25R2573E1BXDBXD2567FR2573E-2.CEL0.0090.0183.15371.880.5790.4070.0141.770.97UAB
26R2562E1BXDBXD2860FR2562E.CEL0.0030.0051.649116.350.5990.3840.0171.370.79Glenn
27R2561E1BXDBXD2960FR2561E.CEL0.0190.0291.95293.320.5830.4020.0152.191Glenn
28R2598E1BXDBXD3161MR2598E.CEL0.0030.0061.989106.480.6090.3760.0151.270.78UTM RW
29R2563E1BXDBXD3263FR2563E.CEL0.0080.0111.547101.520.6190.3670.0141.50.8UTM RW
30R2542E1BXDBXD3367FR2542E.CEL0.0100.0162.12897.080.5650.4180.0161.910.93UTM RW
31R2585E1BXDBXD3460MR2585E.CEL0.0070.0142.6475.130.5830.40.0171.250.77Glenn
32R2532E1BXDBXD3862MR2532E.CEL0.0020.0032.03893.650.5980.3870.0151.370.8UTM RW
33R2574E1BXDBXD3970FR2574E.CEL0.0010.0041.98190.640.6120.3730.0151.390.78UTM RW
34R2590E1BXDBXD4060MR2590E.CEL0.0040.0072.70877.30.5910.3930.0151.40.77Glenn
35R2596E1BXDBXD4259MR2596E.CEL0.0130.0172.632108.460.590.3960.0151.240.8Glenn
36R2605E1BXDBXD4379MR2607E.CEL0.0030.0061.817131.220.6050.3820.0131.320.8UTM RW
37R2594E1BXDBXD4463FR2594E.CEL0.0040.0091.766117.330.5980.3880.0141.350.85UTM RW
38R2592E1BXDBXD4562MR2592E.CEL0.0020.0041.85106.160.6010.3860.0131.430.85UTM RW
39R2606E1BXDBXD4878MR2606E.CEL0.0030.0102.556106.160.5890.3970.0141.350.83UTM RW
40R2603E1BXDBXD5166FR2603E.CEL0.0030.0092.488115.160.5770.4080.0151.240.79UTM RW
41R2534E2BXDBXD6170FR2534E2.CEL0.0300.0282.473117.760.5790.4060.0151.420.79UTM RW
42R2611E1BXDBXD6468MR2611E.CEL0.0130.0222.29291.990.580.4050.0151.571.06UTM RW
43R2583E1BXDBXD6560MR2583E.CEL0.0050.0102.49270.430.5690.4150.0161.671.01UTM RW
44R2536E2BXDBXD6664FR2536E2.CEL0.0390.0652.74108.620.5610.4230.0171.280.79UTM RW
45R2551E1BXDBXD6867FR2551E.CEL0.0370.0392.49392.380.5430.4410.0162.911.55UTM RW
46R2593E1BXDBXD6959FR2593E.CEL0.0080.0131.672127.60.5920.3950.0131.470.92UTM RW
47R2537E2BXDBXD7059MR2537E2.CEL0.0460.0442.9398.660.580.4050.0161.290.75UTM RW
48R2565E1BXDBXD7561FR2565E.CEL0.0090.0171.79101.680.580.4050.0152.313.47UTM RW
49R2579E1BXDBXD8065FR2579E.CEL0.0050.0102.41972.130.5920.3940.0151.730.82UTM RW
50R2540E1BXDBXD8763MR2540E.CEL0.0130.0162.33393.150.6110.3740.0141.220.81UTM RW
51R2545E1BXDBXD8967MR2546E.CEL0.0460.0461.667104.760.5620.4230.0153.69.84UTM RW
52R2578E2BXDBXD9061FR2578E2.CEL0.0330.0342.78592.270.5860.3980.0161.520.77UTM RW
53R2554E1BXDBXD9667MR2554E.CEL0.0040.0042.17793.020.6020.3830.0151.460.77UTM RW
54R2577E1BXDBXD9755MR2577E.CEL0.0190.0162.0776.580.5950.3910.0141.871.29UTM RW
55R2595E1GDP129S1/SvImJ59FR2595E.CEL0.0170.0211.792115.390.610.3750.0151.460.77UTM RW
56R2533E1GDP129S1/SvImJ60MR2533E.CEL0.0210.0132.10793.550.5790.4050.0161.370.78UTM RW
57R2546E1GDPA/J66FR2545E.CEL0.0180.0141.98995.590.5860.3970.0171.470.78UTM RW
58R0754E2GDPA/J60MR0754E.CEL0.0140.0162.71885.630.5980.3870.0151.360.76JAX
59R1676E1GDPBALB/cByJ83FR1676E.CEL0.0420.0412.68598.370.5890.3960.0151.460.74JAX
60R1672E1GDPBALB/cByJ83MR1672E.CEL0.0220.0222.216110.520.5990.3860.0151.260.8JAX
61R1700E1GDPC3H/HeJ83FR1700E.CEL0.0900.0922.97868.770.6080.3790.0141.480.78UTM RW
62R1704E1GDPC3H/HeJ83MR1704E.CEL0.0860.0892.58188.290.6010.3860.0131.380.84UTM RW
63R2564E1GDPCAST/Ei64FR2564E.CEL0.0780.0641.93788.890.5850.3990.0161.60.77JAX
64R2580E1GDPCAST/Ei64MR2580E.CEL0.0760.0672.08994.640.5820.4010.0171.40.76JAX
65R2636E1GDPKK/HIJ64FR2636E.CEL0.0230.0262.6193.10.5890.3950.0151.390.76UTM RW
66R2637E1GDPKK/HIJ64MR2637E.CEL0.0390.0202.189102.780.5940.390.0151.30.79UTM RW
67R0999E1GDPLG/J57FR0999E.CEL0.0120.0122.44882.090.5940.3910.0151.380.79UTM RW
68R1004E1GDPLG/J65MR1004E.CEL0.0130.0152.43891.710.5870.3980.0151.380.79UTM RW
69R1688E1GDPNOD/LtJ66FR1688E.CEL0.0170.0192.66497.650.5860.3990.0151.260.8JAX
70R2566E1GDPNOD/LtJ76MR2566E-2.CEL0.0190.0253.03169.440.5980.3880.0151.380.75UTM RW
71R2550E1GDPNZO/HlLtJ96MR2550E.CEL0.0230.0151.79487.160.6070.3780.0151.520.82JAX
72R2535E1GDPNZO/HlLtJ62FR2535E.CEL0.0460.0251.89385.670.6040.3820.0141.410.85JAX
73R2634E1GDPPWD/PhJ62FR2635E.CEL0.0770.0693.29289.80.5590.4250.0161.570.81JAX
74R2635E1GDPPWD/PhJ62MR2634E.CEL0.0880.0813.72280.050.5420.4410.0171.530.85JAX
75R2544E1GDPPWK/PhJ63FR2544E.CEL0.1060.1002.196107.510.5490.4350.0171.360.82JAX
76R2549E1GDPPWK/PhJ83MR2549E.CEL0.0650.0482.27583.80.5730.4120.0151.570.83JAX
77R2368E1GDPWSB/EiJ67FR2368E.CEL0.0250.0282.56785.70.5950.3910.0141.290.74UTM RW
78R2547E1GDPWSB/EiJ67MR2547E.CEL0.0320.0212.13590.040.5820.4010.0161.320.77UTM RW

    About the array platfrom :

Affymetrix Mouse Genome 430v2: The 430A and B array pairs collectively consist of 992936 25-nucleotide probes that estimate the expression of approximately 39,000 transcripts (many are essentially duplicates). The array sequences were selected late in 2002 using Unigene Build 107. The arrays nominally contain the same probe sequence as the 430 2.0 series. However, we have found that roughy 75000 probes differ between those on A and B arrays and those on the 430 2.0.

    About data processing:

Probe (cell) level data from the CEL file: These CEL values produced by GCOS are 75% quantiles from a set of 91 pixel values per cell.
  • Step 1: We added an offset of 1.0 unit to each cell signal to ensure that all values could be logged without generating negative values. We then computed the log base 2 of each cell.
  • Step 2: We performed a quantile normalization of the log base 2 values for the total set of 105 arrays (processed as two batches) using the same initial steps used by the RMA transform.
  • Step 3: We computed the Z scores for each cell value.
  • Step 4: We multiplied all Z scores by 2.
  • Step 5: We added 8 to the value of all Z scores. The consequence of this simple set of transformations is to produce a set of Z scores that have a mean of 8, a variance of 4, and a standard deviation of 2. The advantage of this modified Z score is that a two-fold difference in expression level corresponds approximately to a 1 unit difference.
  • Step 6: We eliminated much of the systematic technical variance introduced by the two batches (n = 34 and n = 71 array pairs) at the probe level. To do this we calculated the ratio of each batch mean to the mean of both batches and used this as a single multiplicative probe-specific batch correction factor. The consequence of this simple correction is that the mean probe signal value for each batch is the same.
  • Step 7a: The 430A and 430B arrays include a set of 100 shared probe sets (a total of 2200 probes) that have identical sequences. These probes and probe sets provide a way to calibrate expression of the 430A and 430B arrays to a common scale. To bring the two arrays into alignment, we regressed Z scores of the common set of probes to obtain a linear regression correction to rescale the 430B arrays to the 430A array. In our case this involved multiplying all 430B Z scores by the slope of the regression and adding or subtracting a small offset. The result of this step is that the mean of the 430A expression is fixed at a value of 8, whereas that of the 430B chip is typically reduced to 7. The average of the merged 430A and 430B array data set is approximately 7.5.
  • Step 7b: We recentered the merged 430A and 430B data sets to a mean of 8 and a standard deviation of 2. This involved reapplying Steps 3 through 5.
  • Step 8: Finally, we computed the arithmetic mean of the values for the set of microarrays for each strain. Technical replicates were averaged before computing the mean for independent biological samples. Note, that we have not (yet) corrected for variance introduced by differences in sex, age, source of animals, or any interaction terms. We have not corrected for background beyond the background correction implemented by Affymetrix in generating the CEL file. We eventually hope to add statistical controls and adjustments for some of these variables.

    Data source acknowledgment:

Support for acquisition of microarray data were generously provided by the NIAAA and its INIA grant program to RWW, Thomas Sutter, and Daniel Goldowitz (U01AA013515, U01AA013499-03S1, U01AA013488, U01AA013503-03S1). Support for the continued development of the GeneNetwork and WebQTL was provided by a NIMH Human Brain Project grant (P20MH062009). All arrays were processed at the University of Memphis by Thomas Sutter and colleagues with support of the INIA Bioanalytical Core.

    Information about this text file:

This text file originally generated by RWW, YHQ, and EJC, Oct 2004. Updated by RWW, Nov 5, 2004; April 7, 2005; RNA/tissue preparation protocol updatedby JLP, Sept 2, 2005; Sept 26, 2005.