An Efficient and Robust Statistical Modeling Approach to Discover Differentially Expressed Genes Using Genomic Expression Profiles

Table 2.

Top 25 Genes More Highly Expressed in ALL Than in AML

Gene description Probe Difference S.E. Z-score P-value Fold diff t-test t rank ut-test ut rank W. r.s. W rank8
C-myb U22376 −3183.79 429.30 −7.42 <0.0001 5.39 0.1827 7 <0.0001 1 0.0371 2
p48 X74262 −1115.97 152.35 −7.33 <0.0001 5.24 0.2268 8 <0.0001 2 0.0911 6
Proteasome iota chain X59417 −3331.23 481.79 −6.91 <0.0001 3.94 0.3777 12 <0.0001 3 0.3756 23
Myosin light chain (alkali) M31211 −408.99 59.50 −6.87 <0.0001 3.57 0.1584 5 <0.0001 4 0.2854 19
Macmarcks HG1612– −2512.37 372.76 −6.74 <0.0001 2.87 0.2555 9 <0.0001 5 0.0432 3
 HT1612
Transcription factor 3  (E2A) M65214 −471.84 70.04 −6.74 <0.0001 1.56 0.1759 6 <0.0001 6 0.8348 >59
Inducible protein L47738 −1055.39 159.50 −6.62 <0.0001 6.64 0.7782 23 0.0012 7 0.0584 5
MB-1 (CD79b) U05259 −3399.02 523.02 −6.50 <0.0001 12.26 1.0000 >28 0.0016 8 1.0000 >59
Crystallin zeta (quinone  reductase) L13278 −118.93 18.33 −6.49 <0.0001 27.36 0.4366 16 0.0018 9 0.5629 37
Transcriptional activator  hSNF2b D26156 −766.23 118.24 −6.48 <0.0001 2.35 0.3940 13 0.0019 10 0.4304 18
Acyl-Coenzyme A  dehydrogenase, C-4 to  C-12 straight chain M91432 −669.68 104.25 −6.42 <0.0001 4.29 1.0000 >28 0.0020 11 0.1626 11
Oncoprotein 18 M31303 −2013.15 314.83 −6.39 <0.0001 2.12 0.0309 2 0.0029 13 0.1874 13
Thymopoietin beta U09087 −125.08 19.56 −6.39 <0.0001 3.43 0.8446 26 0.0023 12 0.3277 22
Cyclin D3 M92287 −3025.10 484.75 −6.24 <0.0001 3.86 1.0000 >28 0.0035 14 0.0911 7
Serine kinase SRPK2 U88666 −105.40 16.98 −6.21 <0.0001 2.08 0.4101 15 0.0044 17 0.7333 44
Transcription factor 3  (E2A) M31523 −1044.46 169.32 −6.17 <0.0001 4.38 1.0000 >28 0.0043 15 0.0371 1
Adenosine triphosphatase,  calcium Z69881 −1809.52 293.38 −6.17 <0.0001 7.26 1.0000 >28 0.0043 16 0.6423 39
IEF SSP 9502 L07758 −278.58 45.28 −6.15 <0.0001 2.06 0.5292 19 0.0052 18 0.2483 18
Minichromosome  maintenance deficient 3 D38073 −598.54 97.78 −6.12 <0.0001 2.85 0.8160 24 0.0055 20 0.4921 29
Cytoplasmic dynein  light chain 1 U32944 −1349.78 221.68 −6.09 <0.0001 5.16 1.0000 >28 0.0054 19 1.0000 >59
Aldehyde reductase 1 X15414 −818.57 135.98 −6.02 <0.0001 2.26 0.0312 20 0.0090 24 0.4305 27
Spectrin, alpha, non-erythrocytic 1  (alpha-fodrin) J05243 −732.24 122.22 −5.99 <0.0001 6.26 1.0000 >28 0.0078 21 0.1626 12
Rabaptin-5 Y08612 −220.83 36.90 −5.98 <0.0001 2.20 1.0000 >28 0.0079 23 1.0000 >59
Topoisomerase (DNA)  II beta Z15115 −2927.58 490.37 −5.97 <0.0001 3.20 1.0000 >28 0.0079 22 0.0678 5
HKR-T1 S50223 −287.95 48.38 −5.95 <0.0001 5.68 0.5208 18 0.0097 25 1.0000 >59
  • Dataset from Golub et al. (1999).

  • P-value computed from Z-score using a modified Bonferroni's correction.

  • P-value obtained from t-test with equivalent variances using a modified Bonferroni's correction.

  • Relative ranking by significance values obtained from t-test with equivalent variances. >28 indicates that the gene was not ranked because the P-value was 1.0.

  • P-value obtained from t-test with unequal variances using a modified Bonferroni's correction.

  • Relative ranking by significance values obtained from t-test with unequal variances.

  • P-value obtained from Wilcoxon rank sum using a modified Bonferroni's correction.

  • Relative ranking by Wilcoxon significance values. >59 indicates that the gene was not ranked because the P-value was 1.0.

This Article

  1. Genome Res. 11: 1227-1236

Preprint Server