Hamilton Eye Institute Mouse Eye M430v2 (April06) RMA Data Set modify this page

Accession number: GN107

    Summary:

SUPERCEDED EYE DATA SET. The HEIMED April 2006 data set provides estimates of mRNA expression in whole eyes of 71 lines of young adult mice generated using 132 Affymetrix M430 2.0 arrays. Data were generated at UTHSC with support from a grant from Dr. Barrett Haik, Director of the Hamilton Eye Institute (HEI). We used pooled RNA samples, usually two independent pools; one male, one female, for each straion. This data set was processed using the RMA protocol. To simplify comparison among different transforms, RMA values of each array have been adjusted to an average expression of 8 units and a standard deviation of 2 units.

    About the cases used to generate this set of data:

We used a set of 55 BXD recombinant inbred strains, 14 conventional inbred strains including C57BL/6J (B6) and DBA/2J (D2), and their reciprocal F1s. BXD strains were generated by crossing C57BL/6J with DBA/2J. The BXDs are particularly useful for systems genetics because both parental strains have been sequenced (8x coverage of B6 and 1.5x coverage of D2). Physical maps in WebQTL incorporate approximately 2 million B6 vs D2 SNPs from Celera Genomics and from the Perlegen-NIEHS sequencing effort. BXD1 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.

Please note that BXD24/TyJ (JAX stock number 000031) used in this study is also known as BXD24b/TyJ and has nearly complete retinal degeneration. BXD24a/TyJ, a 1988 F80 stock that has now been rederived, does not have retinal degeneration (stock number 005243) and is an ideal coisogenic control.

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.

Mouse Diversity Panel (MDP). In addition to the BXD strains, we have profiled a MDP consisting 14 inbred strains and a pair of reciprocal F1 hybrids; B6D2F1 and D2B6F1. These strains were selected for several reasons:

  • genetic and phenotypic diversity, including use by the Phenome Project
  • their use in making genetic reference populations including recombinant inbred strains, cosomic strains, congenic and recombinant congenic strains
  • their use by the Complex Trait Consortium to make the Collaborative Cross (Nairobi/Wellcome, Oak Ridge/DOE, and Perth/UWA)
  • genome sequence data from three sources (NHGRI, Celera, and Perlegen-NIEHS)
  • availability from The Jackson Laboratory

All eight parents of the Collaborative Cross (129S1/SvImJ, A/J, C57BL/6J, CAST/EiJ, NOD/LtJ, NZO/HlLtJ, PWK/PhJ, and WSB/EiJ) have been included in the MDP. Twelve MDP strains have been sequenced, or are currently being resequenced by Perlegen for the NIEHS.

  1. 129S1/SvImJ
        Collaborative Cross strain sequenced by NIEHS; background for many knockouts (R1 ES cell line); Phenome Project A list. This strain (JAX No 002448, aka 129S1/Sv-++Kitl/+) also carries hypopigmentation mutations (white bellied chinchilla) of the tyrosinase gene on Chr 7 and a mutant allele of the steel (Kitl) gene.
  2. A/J
        Collaborative Cross strain sequenced by Perlegen/NIEHS; parent of the AXB/BXA panel. A tyrosinase-negative albino (c) mutant
  3. BALB/cByJ
         Sequenced by NIEHS; maternal parent of the CXB panel; Phenome Project A list. A tyrosinase-negative albino (c) mutant
  4. C3H/HeJ
        Sequenced by Perlegen/NIEHS; paternal parent of the BXH panel; Phenome Project A list. Important to note for this Eye dataset, C3H/HeJ is a Pdeb6 mutant with near total photoreceptor loss at maturity.
  5. C57BL/6J
        Sequenced by NHGRI; parental strain of AXB/BXA, BXD, and BXH; Phenome Project A list
  6. CAST/Ei
        Collaborative Cross strain sequenced by NIEHS; Phenome Project A list
  7. DBA/2J
        Sequenced by Perlegen/NIEHS and Celera; paternal parent of the BXD panel; Phenome Project A list
  8. KK/HlJ
        Sequenced by Perlegen/NIEHS
  9. LG/J
        Paternal parent of the LGXSM panel
  10. NOD/LtJ
        Collaborative Cross strain sequenced by NIEHS; Phenome Project B list; diabetic
  11. NZO/HlLtJ
        Collaborative Cross strain
  12. PWD/PhJ
        Sequenced by Perlegen/NIEHS; parental strain for a consomic set by Forjet and colleagues
  13. PWK/PhJ
        Collaborative Cross strain; Phenome Project D list
  14. WSB/EiJ
        Collaborative Cross strain sequenced by NIEHS; Phenome Project C list
  15. B6D2F1 and D2B6F1, aka F1 in some graphs and tables
    F1 hybrids generated by crossing C57BL/6J with DBA/2J. These reciprocal F1 can be used to detect some imprinted genes.

    About the tissue used to generate this set of data:

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.

Each array was hybridized with a pool of RNA from 4 to 8 eyes from 2 to 4 animals. RNA was extracted at UTHSC by Zhiping Jia. If tissue was saved for RNA extraction at a later time, eyes were placed directly in RNAlater (Ambion, Inc.) and treated per the manufacturer’s directions. If eyes were used for immediate RNA extraction then we proceeded immediately to the next steps.

Dissecting and preparing eyes for RNA extraction

  1. Place eyes for RNA extraction in RNA STAT-60 (Tel-Test Inc.) and process per manufacturer’s instructions (in brief form below).
  2. Store RNA in 75% ethanol at –80 deg. C until use.

Total RNA was extracted with RNA STAT-60 (Tel-Test Inc.) 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. centrifuged 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 VA Medical Center, Memphis, Rheumatology Disease Research Core Center lead by Drs. John Stuart and Weikuan Gu. All processing steps were performed by Dr. Yan Jiao. In brief, samples were purified using a standard sodium acetate in alcohol method (recommended by Affymetrix). The RNA quality was checked using a 1% agarose gel. The 18S and 28S bands had to be clear and the 28S band had to be more prominent. RNA concentation was measured using a spectrophotometer. The 260/280 ratios had to be greater than 1.7, and the majority were 1.8 or higher. We used a total of 8 micrograms of RNA as starting amount for cDNA synthesis using a standard Eberwine T7 polymerase method (Superscript II RT, Invitrogen Inc., Affy Part No 900431, GeneChip Expression 3' Amplification One-Cyle cDNA Synthesis Kit). The Affymetrix IVT labeling kit (Affy 900449) was used to generate labeled cRNA. At this point the cRNA was evaluated again using both the 260/280 ratio (values of 2.0 or above were acceptable) and 1% agarose gel inspection of the product (a size range from 200 to 7000 bp is considered suitable for use). We used 45 micrograms of labeled cRNA for fragmentation. Those samples that passed both QC steps (<10% usually fail) were then sheared using a fragmentation buffer included in the Affymetrix GeneChip Sample Cleanup Module (Part No.900371). After fragmentation, samples were either stored at -80 deg. centrigrade until use (one third roughly) or were immediately used for hybridization.

Dealing with Ocular Pigmentation: Variable ocular pigmentation is a potential confound in a study of the whole eye transcriptome. Even the most careful RNA preparations taken from brown and beige colored mice tend to have faint residual pigmentation that affects hybridization signal. The key determinant of this interesting effect is the Tyrp1 (brown) locus on Chr 4 at about 80 Mb. Loci on Chr 4 that map at this location should be considered with skepticism and reviewed carefully. To address this problem Yan Jiao purified total RNA a second time using the Qiagen RNeasy MinElute Cleanup Kit (Cat No. 74204). This was done for 8 colored samples (R2534E2, R2578E2, R1441E2, R2537E2, R2536E2, R2589E2, R2539E2) in the first batch of arrays (the November 05 data set) of which 5 were finally included in this data set (cases in which the strain ID is labeled with asterisks in the table). This same protocol was used for all samples in the second batch added in April 2006.

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. The first and second batches of array data, collectively represents a reasonably well balanced sample of males and females belonging to 62 strains, but without within-strain-by-sex replication. Six strains are represented only by male sample pools (BXD15, 28, 29, 55, 98, and DBA/2J. Four strains are represented only by a female pool sample (BXD1, 27, 73 and 86). Please use the probe sets 1427262_at (Xist, high in females) and probe set 1426438_at (Ddx3y, high in males) as quantitative surrogates for the sex balance in this data set.

Batch Structure: This data set consists of a two batches: the original batch that makes up the November 2005 data set and a new batch of 63 arrays (R0857E through R2649E, and R2682E through R2742E, non-consecutive identifiers) run in January 2006 by Dr. Yan Jiao. The arrays in the two batches are from two different lots. All arrays in the second batch were from Lot 4016879 (expiration date 12.28.06). We started working with a total of 140 arrays that passed initial crude quality control based on RNA quality and initial Affymetrix report file information such as 3'/5' ratio, scale factor, and percent present calls. A total of 130 arrays were finally approved for inclusion in this April 2006 data set. The complex normalization procedure is described below.

The table below summarizes information on strain, age, sex, original CEL filename, Affymetrix quality control values, and source of mice. Columns labeled "PDNN 2Z outlier" and "RMA 2Z outlier" list the fraction of probe sets with values that deviated more than 2 z units from the mean. Scale factor, background average, present, absent, marginal and 3'/5' ratios for actin and Gapdh were collated from the Affymetrix Report (RPT) files.

IN PROGRESS: PLEASE NOTE THAT THIS TABLE IS NOW BEING UPDATED TO INCLUDE BATCH 2 OF EARLY 2006.

ID

tube ID

group_type

 Strain

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

1

R2533E1

GDP

129S1/SvImJ

60

M

R2533E.CEL

0.025

0.028

2.11

94

57.90%

40.50%

1.60%

1.37

0.78

UTM RW

2

R2595E1

GDP

129S1/SvImJ

59

F

R2595E.CEL

0.033

0.036

1.79

115

61.00%

37.50%

1.50%

1.46

0.77

UTM RW

3

R0754E2

GDP

A/J

60

M

R0754E.CEL

0.027

0.03

2.72

86

59.80%

38.70%

1.50%

1.36

0.76

JAX

4

R2546E1

GDP

A/J

66

F

R2545E.CEL

0.024

0.029

1.99

96

58.60%

39.70%

1.70%

1.47

0.78

UTM RW

5

R2601E1

GDP BXD

B6D2F1

73

F

R2601E.CEL

0.007

0.008

2.55

92

58.90%

39.60%

1.50%

1.44

0.78

UTM RW

6

R2602E1

GDP BXD

B6D2F1

73

M

R2602E.CEL

0.003

0.008

2.60

84

59.70%

38.80%

1.50%

1.37

0.78

UTM RW

7

R1672E1

GDP

BALB/cByJ

83

M

R1672E.CEL

0.043

0.039

2.22

111

59.90%

38.60%

1.50%

1.26

0.80

JAX

8

R1676E1

GDP

BALB/cByJ

83

F

R1676E.CEL

0.083

0.085

2.69

98

58.90%

39.60%

1.50%

1.46

0.74

JAX

9

R2581E1

BXD

BXD11

65

F

R2581E.CEL

0.009

0.021

1.94

89

62.10%

36.40%

1.60%

1.55

0.81

UTM RW

10

R2543E1

BXD

BXD12

63

M

R2543E.CEL

0.018

0.017

1.61

118

58.60%

39.90%

1.60%

1.43

0.77

UTM RW

11

R2586E1

BXD

BXD13

60

F

R2586E.CEL

0.259

0.258

2.01

74

56.40%

42.00%

1.60%

2.85

3.81

Glenn

12

R2557E1

BXD

BXD14

60

F

R2557E.CEL

0.012

0.027

1.83

99

62.50%

36.10%

1.40%

1.31

0.78

Glenn

13

R2567E1

BXD

BXD16

60

M

R2567E.CEL

0.048

0.058

2.24

82

56.70%

41.60%

1.70%

1.37

0.75

Glenn

14

R2559E1

BXD

BXD18

59

M

R2559E.CEL

0.01

0.012

1.65

104

60.80%

37.70%

1.50%

1.27

0.78

Glenn

15

R2560E1

BXD

BXD19

60

F

R2560E.CEL

0.009

0.012

1.79

98

60.90%

37.50%

1.60%

1.35

0.80

Glenn

16

R2597E1

BXD

BXD2

61

M

R2597E.CEL

0.005

0.012

2.37

94

60.30%

38.30%

1.50%

1.34

0.77

Glenn

17

R2584E1

BXD

BXD20

59

F

R2584E.CEL

0.011

0.017

2.07

84

59.30%

39.10%

1.60%

1.40

0.76

Glenn

18

R2541E2

BXD

BXD21

61

M

R2541E2.CEL

0.049

0.084

2.63

125

56.00%

42.40%

1.50%

1.29

0.78

UTM RW

19

R2553E1

BXD

BXD22

58

F

R2553E.CEL

0.004

0.01

1.95

111

59.90%

38.50%

1.50%

1.28

0.76

Glenn

20

R2558E1

BXD

BXD23

60

F

R2558E-2.CEL

0.018

0.027

1.91

115

59.90%

38.80%

1.40%

1.20

0.82

Glenn

21

R2589E2

BXD

BXD24

59

M

R2589E2.CEL

0.132

0.176

2.61

112

57.50%

40.90%

1.60%

1.24

0.80

Glenn

22

R2573E1

BXD

BXD25

67

F

R2573E-2.CEL

0.055

0.063

3.15

72

57.90%

40.70%

1.40%

1.77

0.97

UAB

23

R2562E1

BXD

BXD29

60

M

R2562E.CEL

0.007

0.01

1.65

116

59.90%

38.40%

1.70%

1.37

0.79

Glenn

24

R2598E1

BXD

BXD31

61

M

R2598E.CEL

0.006

0.013

1.99

106

60.90%

37.60%

1.50%

1.27

0.78

UTM RW

25

R2563E1

BXD

BXD32

63

F

R2563E.CEL

0.023

0.025

1.55

102

61.90%

36.70%

1.40%

1.50

0.80

UTM RW

26

R2542E1

BXD

BXD33

67

F

R2542E.CEL

0.058

0.062

2.13

97

56.50%

41.80%

1.60%

1.91

0.93

UTM RW

27

R2585E1

BXD

BXD34

60

M

R2585E.CEL

0.024

0.032

2.64

75

58.30%

40.00%

1.70%

1.25

0.77

Glenn

28

R2532E1

BXD

BXD38

62

M

R2532E.CEL

0.002

0.006

2.04

94

59.80%

38.70%

1.50%

1.37

0.80

UTM RW

29

R2574E1

BXD

BXD39

70

F

R2574E.CEL

0.003

0.008

1.98

91

61.20%

37.30%

1.50%

1.39

0.78

UTM RW

30

R2590E1

BXD

BXD40

60

M

R2590E.CEL

0.007

0.012

2.71

77

59.10%

39.30%

1.50%

1.40

0.77

Glenn

31

R2596E1

BXD

BXD42

59

M

R2596E.CEL

0.016

0.03

2.63

108

59.00%

39.60%

1.50%

1.24

0.80

Glenn

32

R2605E1

BXD

BXD43

79

M

R2607E.CEL

0.006

0.01

1.82

131

60.50%

38.20%

1.30%

1.32

0.80

UTM RW

33

R2594E1

BXD

BXD44

63

F

R2594E.CEL

0.014

0.024

1.77

117

59.80%

38.80%

1.40%

1.35

0.85

UTM RW

34

R2592E1

BXD

BXD45

62

M

R2592E.CEL

0.005

0.011

1.85

106

60.10%

38.60%

1.30%

1.43

0.85

UTM RW

35

R2606E1

BXD

BXD48

78

M

R2606E.CEL

0.007

0.015

2.56

106

58.90%

39.70%

1.40%

1.35

0.83

UTM RW

36

R2591E1

BXD

BXD5

60

F

R2591E.CEL

0.052

0.014

1.70

136

58.50%

40.00%

1.50%

1.33

0.78

Glenn

37

R2603E1

BXD

BXD51

66

F

R2603E.CEL

0.007

0.02

2.49

115

57.70%

40.80%

1.50%

1.24

0.79

UTM RW

38

R2570E1

BXD

BXD6

65

F

R2570E.CEL

0.013

0.017

1.99

87

58.50%

40.00%

1.50%

1.46

0.76

UTM RW

39

R2534E2

BXD

BXD61

70

F

R2534E2.CEL

0.03

0.058

2.47

118

57.90%

40.60%

1.50%

1.42

0.79

UTM RW

40

R2611E1

BXD

BXD64

68

M

R2611E.CEL

0.067

0.068

2.29

92

58.00%

40.50%

1.50%

1.57

1.06

UTM RW

41

R2583E1

BXD

BXD65

60

M

R2583E.CEL

0.027

0.03

2.49

70

56.90%

41.50%

1.60%

1.67

1.01

UTM RW

42

R2536E2

BXD

BXD66

64

F

R2536E2.CEL

0.067

0.139

2.74

109

56.10%

42.30%

1.70%

1.28

0.79

UTM RW

43

R2551E1

BXD

BXD68

67

F

R2551E.CEL

0.294

0.291

2.49

92

54.30%

44.10%

1.60%

2.91

1.55

UTM RW

44

R2593E1

BXD

BXD69

59

F

R2593E.CEL

0.027

0.038

1.67

128

59.20%

39.50%

1.30%

1.47

0.92

UTM RW

45

R2537E2

BXD

BXD70

59

M

R2537E2.CEL

0.049

0.092

2.93

99

58.00%

40.50%

1.60%

1.29

0.75

UTM RW

46

R2565E1

BXD

BXD75

61

F

R2565E.CEL

0.118

0.124

1.79

102

58.00%

40.50%

1.50%

2.31

3.47

UTM RW

47

R2538E1

BXD

BXD8

77

F

R2538E.CEL

0.033

0.056

1.91

102

61.20%

37.30%

1.50%

1.52

0.79

UTM RW

48

R2579E1

BXD

BXD80

65

F

R2579E.CEL

0.013

0.026

2.42

72

59.20%

39.40%

1.50%

1.73

0.82

UTM RW

49

R2540E1

BXD

BXD87

63

M

R2540E.CEL

0.014

0.034

2.33

93

61.10%

37.40%

1.40%

1.22

0.81

UTM RW

50

R2545E1

BXD

BXD89

67

M

R2546E.CEL

0.266

0.257

1.67

105

56.20%

42.30%

1.50%

3.60

9.84

UTM RW

51

R2569E1

BXD

BXD9

67

M

R2569E.CEL

0.256

0.239

1.75

87

55.10%

43.40%

1.50%

2.82

3.14

UTM RW

52

R2578E2

BXD

BXD90

61

F

R2578E2.CEL

0.041

0.062

2.79

92

58.60%

39.80%

1.60%

1.52

0.77

UTM RW

53

R2554E1

BXD

BXD96

67

M

R2554E.CEL

0.005

0.008

2.18

93

60.20%

38.30%

1.50%

1.46

0.77

UTM RW

54

R2577E1

BXD

BXD97

55

M

R2577E.CEL

0.065

0.069

2.07

77

59.50%

39.10%

1.40%

1.87

1.29

UTM RW

55

R1700E1

GDP

C3H/HeJ

83

F

R1700E.CEL

0.152

0.168

2.98

69

60.80%

37.90%

1.40%

1.48

0.78

UTM RW

56

R1704E1

GDP

C3H/HeJ

83

M

R1704E.CEL

0.154

0.165

2.58

88

60.10%

38.60%

1.30%

1.38

0.84

UTM RW

57

R0872E2

GDP BXD

C57BL/6J

66

M

R0872E.CEL

0.014

0.023

3.13

89

58.90%

39.60%

1.50%

1.30

0.79

UTM RW

58

R2607E1

GDP BXD

C57BL/6J

67

F

R2605E.CEL

0.008

0.018

2.43

115

58.60%

40.00%

1.40%

1.31

0.76

UTM RW

59

R2564E1

GDP

CAST/Ei

64

F

R2564E.CEL

0.124

0.105

1.94

89

58.50%

39.90%

1.60%

1.60

0.77

JAX

60

R2580E1

GDP

CAST/Ei

64

M

R2580E.CEL

0.123

0.109

2.09

95

58.20%

40.10%

1.70%

1.40

0.76

JAX

61

R2600E1

GDP BXD

D2B6F1

72

F

R2600E.CEL

0.008

0.02

2.47

95

58.10%

40.20%

1.70%

1.41

0.78

UTM RW

62

R2604E1

GDP BXD

D2B6F1

69

M

R2604E.CEL

0.005

0.014

2.66

90

59.40%

39.20%

1.50%

1.28

0.79

UTM RW

63

R2572E1

GDP BXD

DBA/2J

65

M

R2572E.CEL

0.091

0.106

2.41

79

55.50%

42.90%

1.60%

1.37

0.79

UTM RW

64

R2636E1

GDP

KK/HIJ

64

F

R2636E.CEL

0.044

0.043

2.61

93

58.90%

39.50%

1.50%

1.39

0.76

UTM RW

65

R2637E1

GDP

KK/HIJ

64

M

R2637E.CEL

0.056

0.036

2.19

103

59.40%

39.00%

1.50%

1.30

0.79

UTM RW

66

R0999E1

GDP

LG/J

57

F

R0999E.CEL

0.021

0.023

2.45

82

59.40%

39.10%

1.50%

1.38

0.79

UTM RW

67

R1004E1

GDP

LG/J

65

M

R1004E.CEL

0.025

0.028

2.44

92

58.70%

39.80%

1.50%

1.38

0.79

UTM RW

68

R1688E1

GDP

NOD/LtJ

66

F

R1688E.CEL

0.028

0.033

2.66

98

58.60%

39.90%

1.50%

1.26

0.80

JAX

69

R2566E1

GDP

NOD/LtJ

76

M

R2566E-2.CEL

0.036

0.04

3.03

69

59.80%

38.80%

1.50%

1.38

0.75

UTM RW

70

R2535E1

GDP

NZO/H1LtJ

62

F

R2535E.CEL

0.037

0.062

1.89

86

60.40%

38.20%

1.40%

1.41

0.85

JAX

71

R2550E1

GDP

NZO/HILtJ

96

M

R2550E.CEL

0.025

0.029

1.79

87

60.70%

37.80%

1.50%

1.52

0.82

JAX

72

R2634E1

GDP

PWD/PhJ

62

F

R2635E.CEL

0.126

0.114

3.29

90

55.90%

42.50%

1.60%

1.57

0.81

JAX

73

R2635E1

GDP

PWD/PhJ

62

M

R2634E.CEL

0.15

0.137

3.72

80

54.20%

44.10%

1.70%

1.53

0.85

JAX

74

R2544E1

GDP

PWK/PhJ

63

F

R2544E.CEL

0.174

0.175

2.20

108

54.90%

43.50%

1.70%

1.36

0.82

JAX

75

R2549E1

GDP

PWK/PhJ

83

M

R2549E.CEL

0.103

0.087

2.28

84

57.30%

41.20%

1.50%

1.57

0.83

JAX

76

R2368E1

GDP

WSB/EI

67

F

R2368E.CEL

0.041

0.047

2.57

86

59.50%

39.10%

1.40%

1.29

0.74

UTM RW

77

R2704E

BXD

BXD1

59

F

R2704E.CEL

0.029

0.03

2.066

139.61

56.60%

41.90%

1.50%

1.31

0.81

GU

78

R2612E

BXD

BXD11

70

M

R2612E.CEL

0.101

0.112

1.83

142.03

58.20%

40.50%

1.40%

1.78

0.81

GU

79

R2742E

BXD

BXD12

71

F

R2742E.CEL

0.073

0.077

2.127

134.14

57.00%

41.60%

1.40%

1.64

0.78

GU

80

R1086E

BXD

BXD23

55

M

R1086E.CEL

0.043

0.034

2.233

125.05

58.60%

39.90%

1.50%

1.43

0.77

GU

81

R2716E

BXD

BXD15

60

M

R2716E.CEL

0.035

0.037

2.015

150.83

56.40%

42.10%

1.60%

1.42

0.81

GU

82

R2711E

BXD

BXD16

61

F

R2711E.CEL

0.032

0.021

1.953

118.53

59.00%

39.60%

1.50%

1.45

0.8

GU

83

R2720E

BXD

BXD18

59

F

R2720E.CEL

0.014

0.019

2.32

99.93

59.50%

39.00%

1.50%

1.33

0.77

GU

84

R2713E

BXD

BXD19

60

M

R2713E.CEL

0.055

0.021

1.67

120.82

60.20%

38.30%

1.50%

1.45

0.8

GU

85

R1231E

BXD

BXD2

64

F

R1231E.CEL

0.044

0.037

2.197

138.73

57.30%

41.30%

1.40%

1.41

0.77

GU

86

R2731E

BXD

BXD20

60

M

R2731E.CEL

0.017

0.019

1.825

147

59.00%

39.50%

1.50%

1.4

0.8

GU

87

R2702E

BXD

BXD21

59

F

R2702E.CEL

0.009

0.008

1.811

128.65

59.40%

39.10%

1.40%

1.26

0.8

GU

88

R2700E

BXD

BXD22

59

M

R2700E.CEL

0.01

0.015

1.858

102.96

61.50%

37.10%

1.30%

1.48

0.79

GU

89

R1128E

BXD

BXD14

65

M

R1128E.CEL

0.037

0.038

2.366

118.39

57.30%

41.30%

1.40%

1.45

0.81

GU

90

R2719E

BXD

BXD24

123

F

R2719E.CEL

0.112

0.111

1.47

140.38

61.50%

37.20%

1.30%

1.38

0.79

GU

91

R2683E

BXD

BXD25

58

M

R2683E.CEL

0.068

0.068

1.777

115.64

58.30%

40.30%

1.40%

2.01

0.79

GU

92

R2703E

BXD

BXD27

60

F

R2703E.CEL

0.008

0.012

1.263

134.78

62.60%

36.10%

1.40%

1.44

0.78

GU

93

R2721E

BXD

BXD28

60

M

R2721E.CEL

0.04

0.048

2.065

157.39

56.10%

42.40%

1.50%

1.31

0.81

GU

94

R1258E

BXD

BXD31

57

F

R1258E.CEL

0.037

0.036

2.063

117.09

59.00%

39.50%

1.50%

1.54

0.78

GU

95

R1216E

BXD

BXD32

76

M

R1216E.CEL

0.05

0.049

2.23

111.99

58.80%

39.80%

1.40%

1.35

0.79

GU

96

R857E

BXD

BXD33

77

M

R857E.CEL

0.078

0.108

1.737

113.98

61.90%

36.70%

1.30%

1.6

0.77

GU

97

R859E

BXD

BXD90

72

M

R859E.CEL

0.028

0.02

1.847

152.22

57.90%

40.70%

1.40%

1.36

0.77

GU

98

R1207E

BXD

BXD66

83

M

R1207E.CEL

0.017

0.012

1.681

136.86

60.40%

38.10%

1.50%

1.45

0.77

GU

99

R2710E

BXD

BXD38

55

F

R2710E.CEL

0.033

0.031

2.112

122.1

58.80%

39.80%

1.40%

1.37

0.78

GU

100

R2695E

BXD

BXD39

59

M

R2695E.CEL

0.018

0.016

1.638

122.7

60.80%

37.80%

1.50%

1.42

0.8

GU

101

R2699E

BXD

BXD40

59

F

R2699E.CEL

0.014

0.015

1.827

105.23

61.70%

36.90%

1.40%

1.42

0.81

GU

102

R2696E

BXD

BXD42

58

F

R2696E.CEL

0.01

0.017

1.622

118.95

62.00%

36.60%

1.50%

1.53

0.79

GU

103

R943E-2

BXD

BXD64

56

F

R943E-2.CEL

0.024

0.021

1.591

141.34

60.10%

38.40%

1.50%

1.32

0.76

GU

104

R967E

BXD

BXD48

64

F

R967E.CEL

0.101

0.052

1.948

130.95

57.30%

41.20%

1.50%

1.63

0.81

GU

105

R2714E

BXD

BXD5

58

M

R2714E.CEL

0.047

0.014

1.404

144.35

60.60%

37.90%

1.50%

1.43

0.79

GU

106

R1042E

BXD

BXD51

62

M

R1042E.CEL

0.028

0.027

2.352

104.12

58.70%

39.90%

1.40%

1.53

0.82

GU

107

R2690E

BXD

BXD55

65

M

R2690E.CEL

0.081

0.067

1.887

164.01

56.10%

42.30%

1.60%

1.43

0.8

GU

108

R2694E

BXD

BXD6

58

M

R2694E.CEL

0.012

0.018

1.983

97.23

61.60%

37.10%

1.30%

1.39

0.82

GU

109

R975E

BXD

BXD70

64

F

R975E.CEL

0.028

0.024

1.841

137.97

58.00%

40.50%

1.40%

1.36

0.79

GU

110

R2684E

BXD

BXD61

62

M

R2684E.CEL

0.031

0.032

2.01

131.03

57.00%

41.50%

1.50%

1.34

0.78

GU

111

R994E

BXD

BXD43

60

F

R994E.CEL

0.013

0.014

1.966

113.12

60.80%

37.80%

1.40%

1.66

0.8

GU

112

R2610E

BXD

BXD44

68

M

R2610E.CEL

0.013

0.009

1.814

142.91

59.00%

39.50%

1.50%

1.35

0.8

GU

113

R2689E

BXD

BXD65

63

F

R2689E.CEL

0.008

0.008

1.721

142.44

59.90%

38.60%

1.50%

1.38

0.76

GU

114

R2727E

BXD

BXD69

65

M

R2727E.CEL

0.01

0.008

1.578

143.86

60.30%

38.30%

1.40%

1.34

0.77

GU

115

R2726E

BXD

BXD68

64

M

R2726E.CEL

0.125

0.025

1.811

153.09

58.70%

39.80%

1.50%

1.39

0.78

GU

116

R2732E

BXD

BXD45

63

F

R2732E.CEL

0.039

0.036

2.154

122.45

56.50%

42.10%

1.40%

1.8

0.83

GU

117

R2709E

BXD

BXD8

61

M

R2709E.CEL

0.012

0.011

1.99

99.79

60.90%

37.60%

1.50%

1.42

0.76

GU

118

R2686E

BXD

BXD80

61

M

R2686E.CEL

0.046

0.05

2.342

119.63

56.00%

42.60%

1.50%

1.38

0.79

GU

119

R2692E

BXD

BXD85

63

F

R2692E.CEL

0.006

0.007

1.423

160.87

60.20%

38.30%

1.40%

1.46

0.79

GU

120

R2715E

BXD

BXD85

91

M

R2715E.CEL

0.007

0.008

1.488

142.6

61.20%

37.30%

1.40%

1.5

0.78

GU

121

R1405E

BXD

BXD86

58

F

R1405E.CEL

0.053

0.052

2.351

119.34

56.40%

42.20%

1.40%

1.64

0.81

GU

122

R2724E

BXD

BXD87

63

F

R2724E.CEL

0.013

0.019

1.906

113.71

60.70%

37.90%

1.40%

1.45

0.79

GU

123

R1451E

BXD

BXD34

61

F

R1451E.CEL

0.01

0.009

1.843

140.05

59.00%

39.50%

1.50%

1.42

0.81

GU

124

R1433E

BXD

BXD89

63

F

R1433E.CEL

0.029

0.026

2.241

115.86

57.70%

40.80%

1.50%

1.41

0.78

GU

125

R2733E

BXD

BXD96

67

F

R2733E.CEL

0.024

0.054

1.7

113.99

62.10%

36.60%

1.30%

1.4

0.78

GU

126

R2649E

BXD

BXD97

74

F

R2649E.CEL

0.029

0.032

2.343

119.04

57.50%

41.20%

1.40%

1.53

0.8

GU

127

R2688E

BXD

BXD98

67

M

R2688E.CEL

0.032

0.03

1.772

145.24

58.50%

40.00%

1.50%

1.48

0.81

GU

128

R877E

BXD

BXD13

76

M

R877E.CEL

0.026

0.067

1.558

125.63

61.20%

37.50%

1.20%

1.42

0.81

GU

129

R1397E-re

BXD

BXD75

58

M

R1397E-re.CEL

0.032

0.01

1.449

189.71

59.60%

39.00%

1.40%

1.39

0.82

GU

130

R2779E

BXD

BXD73

64

F

R2779E.CEL

0.012

0.038

1.746

121.11

59.60%

39.00%

1.40%

1.5

0.8

GU

131

R2708E

BXD

BXD9

60

F

R2708E.CEL

0.024

0.045

1.966

126.46

57.70%

40.70%

1.50%

1.4

0.84

GU

132

R2547E1

GDP

WSB/Ei

67

M

R2547E.CEL

0.041

0.039

2.14

90

58.20%

40.10%

1.60%

1.32

0.77

UTM RW

    About downloading this data set:

This data set is not yet available as a bulk download. Please contact Robert W. Williams to request special data access.

    About the array platfrom:

Affymetrix Mouse Genome 430 2.0 arrays: The 430 2.0 array consists of 992936 25-nucleotide probes that estimate the expression of approximately 39,000 transcripts (many are near duplicates). The array sequences were selected late in 2002 using Unigene Build 107. The array nominally contains the same probe sequences as the old M430A and 430B array pair. However, we have found that roughy 75000 probes differ between those on A and B arrays and those on the new 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. The CEL files were processed using the RMA protocol. We processed the two batches together in RMA.
  • 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 arrays 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 7: Finally, when appropriate, 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.

After RMA processing all arrays were rank-order normalized. This second round of quantile normalization removes much residual non-linearity across arrays and forces every array to have the same distribution of values as the mean of all arrays. Comparative array data quality was then evaluated in DataDesk. Outlier arrays were flagged by visual inspection in DataDesk, usually by means of an analysis of scatter plots and more quantitatively by generating a correlation matrix of all arrays. Those arrays with mean correlation <0.96 versus all other arrays indicates trouble or a biological outlier). In some cases, outliers were expected, such as samples from strains with retinal degeneration (C3H/HeJ and BXD24) and samples from wild subspecies such as CAST/Ei, PWD/Ph, and PWK/Ph. However, when arrays were anomolous both within strain and across strains, they were often simply discarded. We tended to keep arrays that "conformed" to the expectation. The assumption in these cases is that anomolous data are much more likely due to experimental problem and errors than to informative biological variation. Approximately 8 arrays total were discarded in batches 1 and 2 combined.

After this process, the acceptable set of arrays was renormalized using all step as above, starting with the original RMA procedure, etc.

We then categorized arrays into XXX major "technical groups" depending on expression patterns as noted in scatterplots. This process of defining technical groups was done in DataDesk by manually "typing" arrays. These technical groups are apparently due to subtle within-batch effect that we do not yet understand and that cannot be corrected by quantile normalization. These XXX major technical groups are not obviously related to strain, sex, age, or any other known biological effect or variable. They are also not obviously related to any of the Affymetrix QC data types (3'/5' ratios, gain, etc.). Once the technical groups were defined, we forced the means of each probe set in the XX technical groups to the same value. This simple process partially removes a technical error of unknown origin in large expression array data sets.

We reviewed the final data set using a new method developed by RW Williams, Jeremy Peirce, and Hongqiang Li. For the full set of 140 arrays that passed standard QC protocols described above, we computed the strain means for the BXD strains, B6, D2, and F1s. Using this set of strain means we then computed LRS scores for all 45101 probe sets and counted the number of transcripts that generated QTLs with LRS values greater than 50. This value (e.g. 1000) represented the QTL harvest for the full data set. We then dropped a single array from the data set (n = 139 arrays), recomputed strain means, and recomputed the number of transcripts with LRS scores great than 50. This value is expected to typically reduce the number of QTLs that reach the criterion level (e.g., 950 QTLs > 50). This process was repeated for every array to obtain an array-specific difference value--the effect of removing that array on the total QTL count. For example, the loss of a single array might cause a decrease in 50 QTLs (1000-950). Values ranged from -90 (good0 to +38 (bad). This procedure is similar in some ways to a jackknife protocol, although we are not using this procedure to esimate an error term, but rather as a final method to polish a data set. By applying this procedure we discovered that a set of XX (7?) arrays could be excluded while simultaneously improving the total number of QTLs with values above 50.

During this final process we discovered that nearly XX arrays in the second batch had been mislabeled at some point in processing. We computed the correct strain membership of each array using a large number of Mendelian probe sets (more than 50) and comparing their match to standard SNP and microsatellite markers and the original array data set of November 2005. This allowed us to rescue a large number of arrays that were of very high quality.

    Data source acknowledgment:

Support for acquisition of microarray data sets was generously provided by Dr. Barrrett Haik, Chair of the Department of Ophthalmology, and director of the Hamilton Eye Institute. Support for the continued development of GeneNetwork was provided by a NIDA/NIMH/NIAAA Human Brain Project grant. All arrays were processed at the VA Medical Center, Memphis by Weikuan Gu.

    Information about this text file:

This text file originally generated by RWW, May 26, 2006. Updated by RWW, May 27, 2006.