| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
From the Division of Viral and Rickettsial Diseases, National Center for Infectious Diseases, Centers for Disease Control and Prevention, U.S. Department of Health and Human Services, Atlanta, Georgia
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Currently, array technology is most useful in establishing broad patterns of gene expression and in screening for differential gene expression. Validation of expression differences is accomplished with an alternate method such as Northern blot hybridization or RNase protection assay. However, these assays are time-consuming, labor-intensive, and require large amounts of RNA (>5 µg total RNA). Conventional reverse transcription-polymerase chain reaction (RT-PCR) can be done with smaller amounts of RNA (2040 ng), but quantification is difficult and relies on endpoint analysis of the PCR product.8, 9, 10 Real-time (kinetic) PCR evaluates product accumulation during the log-linear phase of the reaction and is currently the most accurate and reproducible approach to gene quantification.9, 10 In this study, we explored the applicability of kinetic RT-PCR as a rapid procedure for the validation of a number of differentially expressed genes identified by HDFA. Because of our interest in the interaction of human papillomavirus (HPV) on cellular gene expression, we used the HDFA expression profiles of two subclones differing in the integration status of HPV (integrated or mixed episomal/integrated) as a model system to test our validation approach. We found that a two-step RT-PCR using SYBR Green I dye detection with product verification by melting curve analysis is rapid, quantitative, and applicable to samples with limited amount of RNA. The method was robust enough to validate relative changes in the expression of a number of genes with varying abundance of transcripts.
| Materials and Methods |
|---|
|
|
|---|
-irradiated (5000 rads) Swiss Mouse 3T3 fibroblast
feeder layers in F-medium (3:1 F12 and Dulbeccos modified Eagles
medium) with 5% fetal bovine serum (FBS).11
CaSki,
a human cervical cancer cell line, was obtained from American Type
Culture Collection (Manassas, VA). CaSki monolayers were grown in
RPMI-1640 medium with 10% FBS and 2.5 mmol/L L-glutamine. Cells were
incubated at 37°C in 5% CO2 and harvested at
60 to 70% confluence. Cultures were washed with phosphate-buffered
saline, followed by 0.02% EDTA to remove the feeder cells. All monolayers were lysed with guanidinium thiocyanate for RNA extraction.12 The total RNA from each sample was divided in half: one half for HDFA after poly(A)+ RNA isolation by using the Oligotex mRNA kit (Qiagen, Santa Clarita, CA) and the other half for HDFA validation by LightCycler (Roche Molecular Biochemicals, Indianapolis, IN). RNA quality and quantity were evaluated by UV spectrophotometry and denaturing formaldehyde agarose gel electrophoresis.13
Gene Expression Profiling by HDFA
Probe synthesis and hybridization conditions optimized for
chemiluminescent detection with HDFA were used as previously
described.14
In brief, cDNA probes were synthesized in a
20 µl RT reaction with 1 µg of
poly(A)+ RNA,
oligo(dT)1218, random hexamers,
digoxigenin-dUTP (Roche Molecular Biochemicals), and SuperScript II
reverse transcriptase enzyme (Life Technologies, Gaithersburg, MD). One
half of the labeled cDNA was used to hybridize the Atlas Human Cancer
cDNA Expression Array (Clontech, Palo Alto, CA). After an overnight
hybridization at 42°C, membranes were washed and hybridization
signals were detected with anti-digoxigenin/alkaline phosphatase
conjugate and CDP-Star substrate. Membranes were exposed to LumiFilm
(Roche Molecular Biochemicals) for 12 minutes after incubating with the
substrate for 1 hour.
The films were scanned and the images were analyzed using BioNumerics (Applied Maths, Kortrijk, Belgium).15 Briefly, images were acquired and converted to Tagged Image File (.tif) format using a flatbed scanner. These array image files were then analyzed in BioNumerics software that subtracted background and normalized intensity on the basis of the lowest negative control as 0 and the highest positive control as 100. The data were copied into Microsoft Excel for further analysis. The lower limit of reliable detection was defined by calculating a threshold value equal to the average intensity of 3 negative controls plus 5 times the SD. Intensities above this threshold were considered positive signals.
Validation of Relative Gene Expression by Kinetic RT-PCR
cDNA Synthesis
Fifty micrograms of total RNA from each sample were treated with
DNase I (0.4 units/µg RNA) according to instructions of the
MessageClean kit (GenHunter Corp., Nashville, TN). One microgram of
DNase-I-treated total RNA was used for cDNA synthesis (20 µl), using
conditions described previously except that random hexamers and
digoxigenin-dUTP were omitted and all dNTPs were maintained at 0.5
mmol/L.14
Primers
Gene-specific primers corresponding to the PCR targets on the
Atlas Human Cancer cDNA Expression Array were obtained from Clontech.
Preliminary experiments were done with each primer pair and CaSki cDNA
to determine the annealing temperature that yielded the greatest amount
of specific product with melting temperature (Tm)
separable from primer-dimer Tm. The acquisition
temperature was set 1 to 2°C below the Tm of
the specific PCR product.16
The experimentally determined
annealing and fluorescent signal acquisition temperatures for each gene
tested in this series of experiments are given in Table 1
.
|
PCR Assay Conditions
DNA Master SYBR Green I mix (containing Taq DNA
polymerase, dNTP, MgCl2, and SYBR Green I dye;
Roche Molecular Biochemicals) was incubated with TaqStart Antibody for
5 minutes at room temperature before the addition of primers and cDNA
template. Each reaction (20 µl) contained 2 µl of the respective
cDNA dilution, primers at 0.4 µmol/L, and MgCl2
at 4 µmol/L. The amplification program consisted of 1 cycle of
95°C with 60-second hold ("hot start") followed by 50 cycles of
95°C with 0-second hold, specified annealing temperature with
5-second hold, 72°C with 18-second hold, and specified acquisition
temperature with 2-second hold (Table 1)
. Amplification was followed by
melting curve analysis using the program run for one cycle at 95°C
with 0-second hold, 65°C with 10-second hold, and 95°C with
0-second hold at the step acquisition mode. A negative control without
cDNA template was run with every assay to assess the overall
specificity.
LightCycler Data Analysis
Unless otherwise mentioned, each assay included duplicate
reactions for each dilution and was repeated once. A relative value for
the initial target concentration in each reaction was determined on the
basis of the kinetic approach9, 10
using the LightCycler
software, version 3. The mean concentration of
glyceraldehyde-3-phosphate dehydrogenase (G3PDH) was used to control
for input RNA because it is considered a stable housekeeping gene and
was detected at the same level in both subclones by HDFA and
LightCycler. The mean G3PDH concentration was determined once for each
cDNA sample and used to normalize all other genes tested from the same
cDNA sample. The relative change in gene expression was recorded as the
ratio of normalized target concentrations in the 1:200 cDNA dilution.
| Results |
|---|
|
|
|---|
|
|
Differences were also seen in the magnitude of differential expression detected by the two methods. For example, among the 14 differentially expressed genes confirmed by LightCycler, 10 genes showed expression differences greater than that determined by HDFA. Notable among these are the expression differences between subclones for fibronectin (eightfold by HDFA and 20-fold by LightCycler), vimentin (twofold by HDFA and 23-fold by LightCycler), and leukocyte interferon-inducible peptide (twofold by HDFA and sevenfold by LightCycler).
| Discussion |
|---|
|
|
|---|
Hot-start PCR mediated by TaqStart antibody and fluorescent signal acquisition at temperatures just below the Tm of specific product allowed SYBR Green I dye-based kinetic-PCR to be sensitive and specific. In this study, the dye-based detection format resulted in an average CV of 12 to 18% for genes with high and low hybridization intensities. This average CV is well below the reported average of 25 to 35% for gene quantification with LightCycler. Most genes could be detected in dilutions of cDNA as low as 1:200 to 1:20,000 from a standard 20-µl RT reaction with 1 µg total RNA. For genes with low expression levels, the amount of cDNA may be increased 10-fold more than suggested in this protocol. Because the assay requires small volumes of cDNA, nearly 100 to 1000 genes can be validated with 1 µg of total RNA. By contrast, validation with Northern blot or RNase protection assay requires at least 5 µg of total RNA per assay, approximately 5000 times more RNA than for the LightCycler assay reported here.
A random selection of genes with varying expression levels detected by HDFA was evaluated by kinetic RT-PCR. Overall, 17 of 24 genes (71%) were confirmed by the LightCycler assay. Both the hybridization intensity and the relative level of gene expression influenced the likelihood that HDFA differences were validated. Genes with weak hybridization signals (<30) and less than fourfold difference in expression were least likely to be validated by the LightCycler. On the other hand, several genes (8 of 10 genes) with high hybridization intensity (>30) and only two- to fourfold difference in expression were validated by the LightCycler. The largest group of differentially expressed genes in many studies using DNA arrays are those with two- to fourfold differences in expression.1, 6, 7 Our data using kinetic RT-PCR suggests that these genes cannot be eliminated as false nor be accepted as true without validation by a secondary procedure.
For genes confirmed as differentially expressed by LightCycler RT-PCR, the level of gene expression differences could be quite different by the two methods. A notable example of such a discrepancy was vimentin gene expression (23-fold differences by LightCycler as opposed to twofold difference by HDFA). The vimentin PCR product includes nucleotides 11641604, a region spanning coding through the 3' UTR. This sequence has 80 to 90% homology over a stretch of 100 nucleotides with other members of the intermediate filament gene family, such as desmin, plasticin, and internexin. This observation suggests that true expression differences for specific members of gene families may be masked by cross-hybridization in microarrays.
Our findings support the use of microarrays as screening tools and emphasize the need for validation of microarray results. The strength of LightCycler assay as a secondary validation procedure lies in its potential to quantify relative change in expression of large number of genes with limited RNA rapidly and precisely.
| Acknowledgments |
|---|
| Footnotes |
|---|
N. T. was supported by the U.S. Department of Energy and CDC Interagency Agreement administered by the Research Participation Program of the Oak Ridge Institute for Science and Education.
Accepted for publication October 27, 2000.
| References |
|---|
|
|
|---|
, ß, or
using oligonucleotide arrays. Proc Natl Acad Sci USA 1998, 95:15623-15628This article has been cited by other articles:
![]() |
A. P.A. van Montfoort, J. P.M. Geraedts, J. C.M. Dumoulin, A. P.M. Stassen, J. L.H. Evers, and T. A.Y. Ayoubi Differential gene expression in cumulus cells as a prognostic indicator of embryo viability: a microarray analysis Mol. Hum. Reprod., March 1, 2008; 14(3): 157 - 168. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Ding, L. Xu, B. D. Jovanovic, I. B. Helenowski, D. L. Kelly, W. J. Catalona, X. J. Yang, M. Pins, and R. C. Bergan The Methodology Used to Measure Differential Gene Expression Affects the Outcome J. Biomol. Tech., December 1, 2007; 18(5): 321 - 330. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Jakob, P. Satorhelyi, C. Lange, V. F. Wendisch, B. Silakowski, S. Scherer, and K. Neuhaus Gene Expression Analysis of Corynebacterium glutamicum Subjected to Long-Term Lactic Acid Adaptation J. Bacteriol., August 1, 2007; 189(15): 5582 - 5590. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Ezaki, R.-J. Guo, H. Li, A. B. Reynolds, and J. P. Lynch The homeodomain transcription factors Cdx1 and Cdx2 induce E-cadherin adhesion activity by reducing beta- and p120-catenin tyrosine phosphorylation Am J Physiol Gastrointest Liver Physiol, July 1, 2007; 293(1): G54 - G65. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. P. O'Grady, H. Mulcahy, J. O'Callaghan, C. Adams, and F. O'Gara Pseudomonas aeruginosa Infection of Airway Epithelial Cells Modulates Expression of Kruppel-Like Factors 2 and 6 via RsmA-Mediated Regulation of Type III Exoenzymes S and Y. Infect. Immun., October 1, 2006; 74(10): 5893 - 5902. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. J. Rampil, D. H. Moller, and A. H. Bell Isoflurane modulates genomic expression in rat amygdala. Anesth. Analg., May 1, 2006; 102(5): 1431 - 1438. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. Hessner, B. Xiang, S. Jia, R. Geoffrey, S. Holmes, L. Meyer, S. Muheisen, and X. Wang Three-color cDNA microarrays with prehybridization quality control yield gene expression data comparable to that of commercial platforms Physiol Genomics, March 13, 2006; 25(1): 166 - 178. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Boyer, W. L. Allen, E. G. McLean, P. M. Wilson, A. McCulla, S. Moore, D. B. Longley, C. Caldas, and P. G. Johnston Pharmacogenomic identification of novel determinants of response to chemotherapy in colon cancer. Cancer Res., March 1, 2006; 66(5): 2765 - 2777. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Cunha-Neto, V. J. Dzau, P. D. Allen, D. Stamatiou, L. Benvenutti, M. L. Higuchi, N. S. Koyama, J. S. Silva, J. Kalil, and C.-C. Liew Cardiac Gene Expression Profiling Provides Evidence for Cytokinopathy as a Molecular Mechanism in Chagas' Disease Cardiomyopathy Am. J. Pathol., August 1, 2005; 167(2): 305 - 313. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. A White and L. A Salamonsen A guide to issues in microarray analysis: application to endometrial biology Reproduction, July 1, 2005; 130(1): 1 - 13. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. W. Tothill, A. Kowalczyk, D. Rischin, A. Bousioutas, I. Haviv, R. K. van Laar, P. M. Waring, J. Zalcberg, R. Ward, A. V. Biankin, et al. An Expression-Based Site of Origin Diagnostic Method Designed for Clinical Application to Cancer of Unknown Origin Cancer Res., May 15, 2005; 65(10): 4031 - 4040. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. T. Lin, M. B. Connelly, C. Amolo, S. Otani, and D. S. Yaver Global Transcriptional Response of Bacillus subtilis to Treatment with Subinhibitory Concentrations of Antibiotics That Inhibit Protein Synthesis Antimicrob. Agents Chemother., May 1, 2005; 49(5): 1915 - 1926. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. E. C. Lutsiak, R. T. Semnani, R. De Pascalis, S. V. S. Kashmiri, J. Schlom, and H. Sabzevari Inhibition of CD4+25+ T regulatory cell function implicated in enhanced immune response by low-dose cyclophosphamide Blood, April 1, 2005; 105(7): 2862 - 2868. [Abstract] [Full Text] [PDF] |
||||
![]() |
S.-T. Chiu, F.-J. Hsieh, S.-W. Chen, C.-L. Chen, H.-F. Shu, and H. Li Clinicopathologic Correlation of Up-regulated Genes Identified Using cDNA Microarray and Real-time Reverse Transcription-PCR in Human Colorectal Cancer Cancer Epidemiol. Biomarkers Prev., February 1, 2005; 14(2): 437 - 443. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Seifarth, O. Frank, U. Zeilfelder, B. Spiess, A. D. Greenwood, R. Hehlmann, and C. Leib-Mosch Comprehensive Analysis of Human Endogenous Retrovirus Transcriptional Activity in Human Tissues with a Retrovirus-Specific Microarray J. Virol., January 1, 2005; 79(1): 341 - 352. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. D. Sgarlato, C. L. Eastman, and H. H. Sussman Panel of Genes Transcriptionally Up-regulated in Squamous Cell Carcinoma of the Cervix Identified by Representational Difference Analysis, Confirmed by Macroarray, and Validated by Real-Time Quantitative Reverse Transcription-PCR Clin. Chem., January 1, 2005; 51(1): 27 - 34. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Pusztai and K. R. Hess Clinical trial design for microarray predictive marker discovery and assessment Ann. Onc., December 1, 2004; 15(12): 1731 - 1737. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. P. Ponnampalam, G. C. Weston, A. C. Trajstman, B. Susil, and P. A.W. Rogers Molecular classification of human endometrial cycle stages by transcriptional profiling Mol. Hum. Reprod., December 1, 2004; 10(12): 879 - 893. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Katahira, T. Takayama, K. Miyanishi, T. Hayashi, T. Ikeda, Y. Takahashi, R. Takimoto, T. Matsunaga, J. Kato, and Y. Niitsu Plasma Glutathione S-Transferase P1-1 as a Prognostic Factor in Patients with Advanced Non-Hodgkin's Lymphoma (Stages III and IV) Clin. Cancer Res., December 1, 2004; 10(23): 7934 - 7940. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. R. Miles, C. E. Farin, K. F. Rodriguez, J. E. Alexander, and P. W. Farin Angiogenesis and Morphometry of Bovine Placentas in Late Gestation from Embryos Produced In Vivo or In Vitro Biol Reprod, December 1, 2004; 71(6): 1919 - 1926. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Zhao, A. Chow, J. Powers, G. Fajardo, and D. Bernstein Microarray analysis of gene expression after transverse aortic constriction in mice Physiol Genomics, September 16, 2004; 19(1): 93 - 105. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. O-charoenrat, V. Rusch, S. G. Talbot, I. Sarkaria, A. Viale, N. Socci, I. Ngai, P. Rao, and B. Singh Casein Kinase II Alpha Subunit and C1-Inhibitor Are Independent Predictors of Outcome in Patients with Squamous Cell Carcinoma of the Lung Clin. Cancer Res., September 1, 2004; 10(17): 5792 - 5803. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Wu, P. de Kievit, L. Vahlkamp, D. Pijnenburg, M. Smit, M. Dankers, D. Melchers, M. Stax, P. J. Boender, C. Ingham, et al. Quantitative assessment of a novel flow-through porous microarray for the rapid analysis of gene expression profiles Nucleic Acids Res., August 27, 2004; 32(15): e123 - e123. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Biedermann, N. Dandachi, M. Trattner, G. Vogl, H. Doppelmayr, E. More, A. Staudach, O. Dietze, and C. Hauser-Kronberger Comparison of Real-Time PCR Signal-Amplified In Situ Hybridization and Conventional PCR for Detection and Quantification of Human Papillomavirus in Archival Cervical Cancer Tissue J. Clin. Microbiol., August 1, 2004; 42(8): 3758 - 3765. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. S. Keller, T. Ezaki, R.-J. Guo, and J. P. Lynch Cdx1 or Cdx2 expression activates E-cadherin-mediated cell-cell adhesion and compaction in human COLO 205 cells Am J Physiol Gastrointest Liver Physiol, July 1, 2004; 287(1): G104 - G114. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Tsuchiya, J. M. Dhahbi, X. Cui, P. L. Mote, A. Bartke, and S. R. Spindler Additive regulation of hepatic gene expression by dwarfism and caloric restriction Physiol Genomics, May 19, 2004; 17(3): 307 - 315. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Slobedman, J. L. Stern, A. L. Cunningham, A. Abendroth, D. A. Abate, and E. S. Mocarski Impact of Human Cytomegalovirus Latent Infection on Myeloid Progenitor Cell Gene Expression J. Virol., April 15, 2004; 78(8): 4054 - 4062. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Dhahbi, H.-J. Kim, P. L. Mote, R. J. Beaver, and S. R. Spindler Temporal linkage between the phenotypic and genomic responses to caloric restriction PNAS, April 13, 2004; 101(15): 5524 - 5529. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Therianos, M. Zhu, E. Pyun, and P. D. Coleman Single-Channel Quantitative Multiplex Reverse Transcriptase-Polymerase Chain Reaction for Large Numbers of Gene Products Differentiates Nondemented from Neuropathological Alzheimer's Disease Am. J. Pathol., March 1, 2004; 164(3): 795 - 806. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Konig, D. Baldessari, N. Pollet, C. Niehrs, and R. Eils Reliability of gene expression ratios for cDNA microarrays in multiconditional experiments with a reference design Nucleic Acids Res., February 13, 2004; 32(3): e29 - e29. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. D. Sloan, B. Nicholson, V. Urquidi, and S. Goodison Detection of Differentially Expressed Genes in an Isogenic Breast Metastasis Model using RNA Arbitrarily Primed-Polymerase Chain Reaction Coupled with Array Hybridization (RAP-Array) Am. J. Pathol., January 1, 2004; 164(1): 315 - 323. [Abstract] [Full Text] [PDF] |
||||
![]() |
S D Jenson, R S Robetorye, S D Bohling, J A Schumacher, J W Morgan, M S Lim, and K S J Elenitoba-Johnson Validation of cDNA microarray gene expression data obtained from linearly amplified RNA Mol. Pathol., December 1, 2003; 56(6): 307 - 312. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Liu, A. R. Bauskin, J. Zaunders, D. A. Brown, S. Pankurst, P. J. Russell, and S. N. Breit Macrophage Inhibitory Cytokine 1 Reduces Cell Adhesion and Induces Apoptosis in Prostate Cancer Cells Cancer Res., August 15, 2003; 63(16): 5034 - 5040. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. G. Bottone Jr., J. M. Martinez, J. B. Collins, C. A. Afshari, and T. E. Eling Gene Modulation by the Cyclooxygenase Inhibitor, Sulindac Sulfide, in Human Colorectal Carcinoma Cells: POSSIBLE LINK TO APOPTOSIS J. Biol. Chem., July 3, 2003; 278(28): 25790 - 25801. [Abstract] [Full Text] [PDF] |
||||
![]() |
D.-S. Bae, R. J. Handa, R. S. H. Yang, and J. A. Campain Gene Expression Patterns as Potential Molecular Biomarkers for Malignant Transformation in Human Keratinocytes Treated with MNNG, Arsenic, or a Metal Mixture Toxicol. Sci., July 1, 2003; 74(1): 32 - 42. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Michalak and M. A. F. Noor Genome-Wide Patterns of Expression in Drosophila Pure Species and Hybrid Males Mol. Biol. Evol., July 1, 2003; 20(7): 1070 - 1076. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Naciff, G. J. Overmann, S. M. Torontali, G. J. Carr, J. P. Tiesman, B. D. Richardson, and G. P. Daston Gene Expression Profile Induced by 17{alpha}-Ethynyl Estradiol in the Prepubertal Female Reproductive System of the Rat Toxicol. Sci., April 1, 2003; 72(2): 314 - 330. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Pedotti, J. J. DeVoss, S. Youssef, D. Mitchell, J. Wedemeyer, R. Madanat, H. Garren, P. Fontoura, M. Tsai, S. J. Galli, et al. Multiple elements of the allergic arm of the immune response modulate autoimmune demyelination PNAS, February 18, 2003; 100(4): 1867 - 1872. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Yang, H. Haddad, C. Tomas, K. Alsaker, and E. T. Papoutsakis A segmental nearest neighbor normalization and gene identification method gives superior results for DNA-array analysis PNAS, February 4, 20 |