* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Download Normalization of gene expression measurements in tumor tissues
Epitranscriptome wikipedia , lookup
Gene therapy wikipedia , lookup
Transposable element wikipedia , lookup
Metagenomics wikipedia , lookup
Primary transcript wikipedia , lookup
Gene nomenclature wikipedia , lookup
Vectors in gene therapy wikipedia , lookup
Pathogenomics wikipedia , lookup
Essential gene wikipedia , lookup
Gene desert wikipedia , lookup
Epigenetics in learning and memory wikipedia , lookup
Public health genomics wikipedia , lookup
Short interspersed nuclear elements (SINEs) wikipedia , lookup
Cancer epigenetics wikipedia , lookup
History of genetic engineering wikipedia , lookup
Epigenetics of neurodegenerative diseases wikipedia , lookup
Gene therapy of the human retina wikipedia , lookup
Quantitative trait locus wikipedia , lookup
RNA interference wikipedia , lookup
X-inactivation wikipedia , lookup
Oncogenomics wikipedia , lookup
Non-coding RNA wikipedia , lookup
RNA silencing wikipedia , lookup
Minimal genome wikipedia , lookup
Biology and consumer behaviour wikipedia , lookup
Genome evolution wikipedia , lookup
Epigenetics of diabetes Type 2 wikipedia , lookup
Microevolution wikipedia , lookup
Polycomb Group Proteins and Cancer wikipedia , lookup
Genome (book) wikipedia , lookup
Genomic imprinting wikipedia , lookup
Site-specific recombinase technology wikipedia , lookup
Ridge (biology) wikipedia , lookup
Long non-coding RNA wikipedia , lookup
Designer baby wikipedia , lookup
Therapeutic gene modulation wikipedia , lookup
Nutriepigenomics wikipedia , lookup
Artificial gene synthesis wikipedia , lookup
Gene expression programming wikipedia , lookup
Epigenetics of human development wikipedia , lookup
Mir-92 microRNA precursor family wikipedia , lookup
Laboratory Investigation (2005) 85, 154–159 & 2005 USCAP, Inc All rights reserved 0023-6837/05 $30.00 www.laboratoryinvestigation.org Normalization of gene expression measurements in tumor tissues: comparison of 13 endogenous control genes Jacques B de Kok1, Rian W Roelofs1, Belinda A Giesendorf1, Jeroen L Pennings2, Erwin T Waas3, Ton Feuth4, Dorine W Swinkels1 and Paul N Span5 1 Department of Clinical Chemistry, University Medical Centre Nijmegen, Nijmegen, The Netherlands; National Institute for Public Health and Environment, Bilthoven, The Netherlands; 3Department of Surgery, University Medical Centre Nijmegen, Nijmegen, The Netherlands; 4Department of Epidemiology and Biostatistics, University Medical Centre Nijmegen, Nijmegen, The Netherlands and 5Department of Chemical Endocrinology, University Medical Centre Nijmegen, Nijmegen, The Netherlands 2 For interpretation of quantitative gene expression measurements in clinical tumor samples, a normalizer is necessary to correct expression data for differences in cellular input, RNA quality, and RT efficiency between samples. In many studies, a single housekeeping gene is used for normalization. However, no unequivocal single reference gene (with proven invariable expression between cells) has been identified yet. As the best alternative, the mean expression of multiple housekeeping genes can be used for normalization. In this study, no attempt was made to determine the gold-standard gene for normalization, but to identify the best single housekeeping gene that could accurately replace the measurement of multiple genes. Expression patterns of 13 frequently used housekeeping genes were determined in 80 normal and tumor samples from colorectal, breast, prostate, skin, and bladder tissues with real-time quantitative RT-PCR. These genes included, large ribosomal protein, b-actin, cyclophilin A, glyceraldehyde-3-phosphate dehydrogenase, phosphoglycerokinase 1, b-2microglobin, b-glucuronidase, hypoxanthine ribosyltransferase (HPRT), TATA-box-binding protein, transferrin receptor, porphobilinogen deaminase, ATP synthase 6, and 18S ribosomal RNA. Principal component analysis was used to analyze these expression patterns, independent of the level of expression. Our approach identified HPRT as the single best reference gene that could be used as an accurate and economic alternative for the measurement of multiple housekeeping genes. We recommend this gene for future studies to standardize gene expression measurements in cancer research and tumor diagnostics until a definite gold standard has been determined. Laboratory Investigation (2005) 85, 154–159, advance online publication, 15 November 2004; doi:10.1038/labinvest.3700208 Keywords: normalization; housekeeping gene; reference gene; control gene; real-time quantitative RT-PCR Quantitative RT-PCR is a frequently applied method to identify genes that correlate with tumor diagnosis or disease prognosis.1,2 In this assay, it is mandatory to relate gene expression to the amount of tissue analyzed, that is, to normalize for optimal comparison of expression levels between tissue samples. Inaccurate normalization results in inadequate quantification and spurious conclusions. Usually, cellular maintenance genes, the so-called housekeeping genes, are selected to normalize for Correspondence: Dr JB de Kok, PhD, Clinical Chemistry/564 AKC, University Medical Centre Nijmegen, PO Box 9101, 6500 HB Nijmegen, The Netherlands. E-mail: [email protected] Received 18 August 2004; revised and accepted 1 October 2004; published online 15 November 2004 the variability between clinical samples. These genes regulate basic and ubiquitous cellular functions and code, for example, for components of the cytoskeleton (b-actin), major histocompatibility complex (b-2-microglobulin), glycolytic pathway (glyceraldehyde-3-phosphate dehydrogenase (GAPDH), phosphoglycerokinase 1), metabolic salvage of nucleotides (hypoxanthine ribosyltransferase), protein folding (cyclophilin), or synthesis of ribosome subunits (rRNA). In many experiments, the expression of these genes is assumed invariable between cells of different samples and used as normalizer without proper validation. Other studies have shown, however, that their individual expression may vary as a result of neoplastic growth,3,4 hypoxia,5 or experimental treatment,6–8 and may seriously influence correct interpretation of results. Normalization of gene expression JB de Kok et al The ‘gold-standard’ normalization gene is still the subject of debate. Currently, the best option is to measure the expression of multiple housekeeping genes and normalize using their mean expression. In cancer research, only few studies attempted to investigate the variation in expression of housekeeping genes between tissue samples. Mostly, only two or three candidate genes were compared.3,4,9–11 In the present study, we attempted to identify a single housekeeping gene that could replace the measurement of multiple genes. Therefore, expression patterns of 13 ordinarily used housekeeping genes were investigated in 80 epithelial tissue samples (normal and tumor tissues). These samples represented colorectal, breast, prostate, skin, and bladder tissues, with tumors ranging from noninvasive to metastatic carcinomas. Quantitative RT-PCR was used for the generation of expression data. Principal component analysis (PCA), linear regression and difference-plot analysis were used to interpret data. Materials and methods increasing TNM stage:14 three T2N0, three T3Nx, and three T4NxMx-tumors. For the skin, three normal tissues, three normal nevi, three dysplastic nevi, three primary melanomas, and four melanoma metastases were selected. For the bladder, four normal urothelia were selected together with 12 urothelial cell carcinomas:14,15 two pTa grade 1, two pTa grade 2, two pT1 grade 2, two pT2 grade 3, two pT3 grade 3, and two pT4 grade 3 with distant metastates. Total RNA was isolated from tissues and reverse transcribed to cDNA as described previously.16 155 Housekeeping Genes A total of 13 housekeeping genes, used in many studies, were selected for gene expression analysis (Table 1). All genes are constitutively expressed in various tissues.9,17,18 To our knowledge, all genes have independent functions in cellular maintenance, and regulation of their expression is assumed not to be related directly. Only GAPDH and PGK share an identical biochemical process in the cell, namely glycolysis. Tissues A total of 16 frozen tissues were selected for each organ (colon, breast, prostate, skin, bladder) from the tissue banks of the Departments of Pathology, Urology, Surgery, and Chemical Endocrinology (UMC Nijmegen, The Netherlands). For the colon, four normal tissues were selected together with 12 tumor tissues with increasing Dukes stage:12 three Dukes A, three Dukes B, three Dukes C, and three Dukes D tumors. For the breast, two normal tissues were selected together with two fibroadenomas and 12 tumors with increasing pathological stage:13 two Tis (carcinoma in situ), one stage 1, one stage 2A, two stage 2B, two stage 3A, two stage 3B, and two stage 4 tumors with metastases. For prostate, four normal tissues were selected together with three benign prostate hyperplasias, and nine tumor tissues with Gene Expression Measurements Measurement of the expression of all genes, except for PBGD and ATP6, was performed using the human endogenous control plate (Applied Biosystems). The human endogenous control plate is a 96-well PCR plate divided into 12 columns of eight identical wells. A total of 11 columns contain TaqMan primers and probes for the detection of 11 different housekeeping genes. The 12th column contains a fixed copy number of an artificial positive control sequence (APC) and corresponding primers and probe. If PCR inhibitors from the sample are present in the PCR reaction, the signal generated by the APC sequence will diminish. Samples containing inhibitors are excluded from further Table 1 Selected housekeeping genes for gene expression analysis Code Gene name Abbreviation Cellular function G1 G2 G3 G4 G5 G6 G7 G8 G9 G10 G11 G12 G13 Large ribosomal protein b-actin Cyclophilin A Glyceraldehyde-3-phosphate dehydrogenase Phosphoglycerokinase 1 b-2-microglobulin b-glucuronidase Hypoxanthine ribosyltransferase TATA-box-binding protein Transferrin receptor Porphobilinogen deaminase ATP synthase 6 18S ribosomal RNA LRP BACT CYC GAPDH PGK B2M BGUS HPRT TBP TfR PBGD ATP6 rRNA Transcription Cytoskeleton Serine-threonine phosphatase inhibitor Glycolysis enzyme Glycolysis enzyme Major histocompatibility complex Exoglycosidase in lysosomes Metabolic salvage of purines Transcription by RNA polymerases Cellular iron uptake Heme synthesis Oxydative phosphorylation Ribosome subunit Laboratory Investigation (2005) 85, 154–159 Normalization of gene expression JB de Kok et al 156 analysis. In the absence of inhibitors, the APC sequence should present a constant cycle threshold (Ct) value. A measure of 10 ml of cDNA of each tissue sample (n ¼ 80) was added to 315 ml H2O and 325 ml Universal Master Mix (Applied Biosystems). This mixture was distributed over a row of 12 wells, with 50 ml mixture in each well. Consequently, 11 housekeeping genes and the APC sequence were amplified for each individual tissue sample with equal input of cDNA for each gene measurement. In total, 10 96-well PCR plates were used for all five tissue-types (5 16 samples). Increase of fluorescence was measured at every PCR cycle in each well by the ABI Prism 7700 sequence detection system (Applied Biosystems). Measurement of the expression of ATP6 and PBGD genes was described previously,9,19 also using an equal input of cDNA. The expression of PBGD and ATP6 could not be measured in prostate tissues due to insufficient sample. After PCR, the number of PCR cycles to reach the fluorescence threshold in each sample was defined as the cycle threshold (Ct). Ct values are proportional to the negative logarithm of the initial amount of input cDNA.20 Ct values of 13 housekeeping genes in one tissue sample were directly related, since the input of cDNA was equal for each PCR reaction. Statistical Analysis Analysis for patterns of gene expression within the set of 13 housekeeping genes was performed with PCA using the SAS package (SAS Institute Inc.). PCA is a linear dimensionality reduction technique, which identifies (orthogonal) directions of maximum variance in the original data. The data are projected in a lower dimension, formed of (a subset of) components with the highest variance.21,22 Principal components are therefore linear combinations of the original variables, orthogonal, and ordered with respect to their variance. The first principal component has the largest variance and represents better than any other linear combination the general differences/similarities between the housekeeping genes. In the analysis of the expression of different housekeeping genes, the housekeeping genes were used as the variables and the Ct values of the tissue samples as the observations. This will group housekeeping genes with similar variation in Ct (ie expression) patterns between all tissue samples. All PCA results were analyzed in two dimensions. To assess the reliability of the PCA analysis on the whole data set we used a split-sample technique to replicate the results. For this, the set of 80 tissues was divided into two subsets, comparable with respect to type of organ and degree of malignancy. Laboratory Investigation (2005) 85, 154–159 Results Quality Control In each endogenous control plate, Ct values for the APC sequence were constant with a within-run standard deviation of Cto0.3, demonstrating that PCR inhibitors were not present and that Ct values of all housekeeping genes could be compared. Only one out of 80 samples (a metastatic melanoma tissue sample) showed a significantly higher Ct-value (ie lower expression) from the others. In addition, one fibroadenoma of the breast contained a very low amount of RNA. Both samples were omitted from further analysis. Principle of Expression-Profile Analysis In our study, the cDNA generated from each sample was split into 13 equal parts for the measurement of 13 different housekeeping genes. The standardized cDNA input for the measurement of each gene linked expression patterns (represented by Ct values) in a sample directly. This provided the basis for our analyses. For example, when the expression patterns of two housekeeping genes were closely related, Ct values increased (or decreased) simultaneously between a series of tissue samples (depending on differences of cDNA amounts between the samples). Subsequently, housekeeping genes could be identified, which show similar expression patterns in all tissue samples. Closely related genes, which fluctuate concurrently, form a cluster with PCA. Principal Component Analysis The relationship between the expression patterns of 13 housekeeping genes was investigated mathematically with PCA of the Ct values of all 78 tissues. In the two-dimensional analysis, the first and second principal components explained 76 and 12% of the variability (88% in total). Analysis with more components was considered superfluous. Three different patterns of expression were observed (Figure 1). One pattern of expression was represented by ATP6 (G12), another by rRNA (G13), and the third pattern encompassed a cluster of 11 housekeeping genes (G1–G11). After repeating the PCA analysis on two halves of data sets (splitsample method), we observed the same clustering of genes in both data sets (data not shown). Linear Regression and Difference Plot The possibility that the expression of a single gene may represent the mean expression of the whole cluster of 11 genes was illustrated by plotting the expression (Ct value) of each individual housekeeping gene vs the mean expression (mean Ct) of Normalization of gene expression JB de Kok et al 157 the remaining housekeeping genes for each sample. After linear regression, coefficients of correlation (r2) were calculated, representing the relationship of the expression of each individual housekeeping gene with the mean expression of the other 10 genes. Coefficients of correlation were presented for each tissue-type separately and for all tissue-types together (Table 2). The power of the data for each tissue-type separately is lower (ie has a higher deviation) than the coefficient of correlation for all tissues together, since only 15 or 16 tissue samples were used for analysis per tissue-type. Clearly, hypoxanthine ribosyltransferase (HPRT) shows a very high correlation coefficient for each tissue-type separately (r2Z0.90), and the best score for all tissues together (r2 ¼ 0.93). The absolute relationship between the expression of HPRT and the mean expression of the other genes can be quantified using a difference plot (Figure 2). Here, the difference between the Ct of HPRT and the mean Ct is plotted against the mean Ct. The mean difference (Ct ¼ 3.0) represents the constant difference in level of expression between HPRT and the mean expression of the other genes: HPRT has an overall lower expression than the mean expression (3.0 extra PCR cycles are needed to reach 1.2 1 fluorescence threshold). Important, however, is the deviation around the mean expression. This represents the accuracy of the HPRT expression measurement compared to the mean expression. For HPRT, two times the standard deviation from the mean difference is 1.3 Ct, showing that measurement of HPRT expression maximally deviates 71.3 Ct ( ¼ PCR cycles) from the expression measurement of the other genes together (for 95% of the tissue samples). The mean difference in expression with corresponding accuracy (2 s.d.) is presented for all genes in Table 3. Discussion A major difficulty in studying differential gene expression is how to normalize for heterogeneity between tissue samples. Ideally, the gene transcript number is corrected for the number of cells analyzed. Unfortunately, counting of cells is not practical for solid epithelial tissues. Alternatively, total RNA can be quantified and used for normalization.5,7,9,10,23,24 However, cellular RNA content may increase with increasing tumor aneuploidy.25,26 Moreover, the majority of total RNA consists of ribosomal and mitochondrial RNA, hence it is not surprising that the expression of rRNA and ATP6 (a mitochondrial-coded gene) correlates with total G12 0.6 0.4 0.2 0 −0.2 0.2 0.4 0.6 0.8 G9 G3 G4 G11 G10G8 G6 G2 G7 1 G1 −0.4 Difference (Ct - mean Ct) Component 2 0.8 G5 G13 −0.6 −0.8 Component 1 Figure 1 PCA of the expression data (Ct values) of 13 housekeeping genes in 78 tissue samples. The distance between genes (G1–G13) is expressed in relative units and inversely related to similarities in expression pattern. Hypoxanthine Ribosyltransferase 5 +2SD 4 3 Mean difference 2 −2SD 1 0 20 22 24 26 28 30 32 34 Mean Ct Figure 2 Difference plot, representing the accuracy of the expression of HPRT compared to the mean expression of 10 other housekeeping genes in all tissue samples. Table 2 Coefficients of correlation between the expression of each individual gene and the mean expression of the remaining 10 genes Gene code Abbreviation G1 G2 G3 G4 G5 G6 G7 G8 G9 G10 G11 LRP BACT CYC GAPDH PGK B2M BGUS HPRT TBP TfR PBGD Colon Breast Prostate Skin Bladder All tissues 0.96 0.92 0.81 0.74 0.91 0.62 0.74 0.92 0.93 0.89 0.87 0.55 0.86 0.66 0.82 0.65 0.53 0.83 0.90 0.72 0.88 0.93 0.44 0.86 0.90 0.80 0.95 0.62 0.96 0.92 0.83 0.97 NA 0.70 0.94 0.95 0.79 0.96 0.70 0.97 0.95 0.92 0.73 0.93 0.62 0.84 0.92 0.90 0.85 0.71 0.84 0.95 0.75 0.81 0.61 0.74 0.82 0.85 0.76 0.88 0.76 0.91 0.93 0.75 0.81 0.81 Data are represented for tissue-types separately and for all tissues together. NA ¼ not analysed. Laboratory Investigation (2005) 85, 154–159 Normalization of gene expression JB de Kok et al 158 Table 3 Independent relationship between the expression of each single housekeeping gene compared to the mean expression of the other 10 genes Gene code Abbreviation G1 G2 G3 G4 G5 G6 G7 G8 G9 G10 G11 LRP BACT CYC GAPDH PGK B2M BGUS HPRT TBP TfR PBGD Mean difference (Ct) Accuracy (2 s.d.) 2.6 5.5 2.3 4.2 0.3 2.2 1.7 3.0 5.4 2.0 0.6 2.8 1.9 2.0 2.3 1.9 2.5 1.4 1.3 2.3 2.0 1.7 The mean difference (Ct–mean Ct) represents the difference in height of expression, whereas the accuracy represents 2 the standard deviation. RNA amount in some studies.3,9 This does not prove that these genes are invariably expressed during tumorigenesis. In addition, total RNA quantification does not correct for RNA quality differences (eg due to storage), nor for differences in reverse transcriptase efficiencies between samples.27 It is now generally accepted that gene expression levels should be normalized to an invariably expressed internal control gene that reflects differences in cellular input, RNA quality, and RT efficiency. However, how can such a gene be found when no gold standard is available to refer its expression to? To circumvent this circular problem we used a different approach. We hypothesized that the mean expression of a large set of housekeeping genes with independent cellular functions would accurately reflect optimal normalization. Therefore, we chose 13 housekeeping genes with independent functions in cellular maintenance. This independence is important, because selection of genes that share identical biochemical pathways could bias analysis. A large set of epithelial tissue samples was selected (n ¼ 80), representing five different organs, each spanning a range from normal tissues to metastatic carcinomas. Expression of each housekeeping gene was measured in the tissue samples with real-time quantitative PCR. No attempt was made to determine the gold-standard gene, but to select a single gene that could replace multiple gene measurements. PCA showed that the expression pattern of ATP6 and rRNA significantly differed from a large cluster of 11 genes and from each other. The difference in expression patterns between the three gene clusters may be associated with the different RNA polymerases that regulate their transcription: RNA polymerase I transcribes rRNA, mitochondrial RNA polymerase transcribes ATP6, and RNA polymerase II transcribes the other 11 genes. Since (1) PCA was not influenced by differences in expression levels between genes but recognized patterns only Laboratory Investigation (2005) 85, 154–159 and (2) regulation of expression of the 13 housekeeping genes was (assumed to be) unrelated, we concluded that ATP6 and rRNA were poor normalization genes and both were excluded from further analysis. Within the large cluster of the 11 remaining genes, the expression pattern of HPRT most accurately reflected the mean expression pattern of the other 10 genes. In the difference plot, an accuracy of 71.3 PCR cycles (Ct) implies that normalization may differ approximately two-fold using the expression of HPRT instead of the mean expression of all genes. This demonstrates that the expression of HPRT could be used as a simple and economic alternative for the measurement of 11 genes simultaneously. Although the combination of more than one gene may improve accuracy, for most research applications this will not be necessary. In exceptional situations, when only little RNA of a clinical sample is available (eg after microdissection) and the gene of interest has a high expression, HPRT expression may be measured (close to) negative and accurate normalization is not possible. In this situation, a high-expression housekeeping gene can be selected from Table 3. For instance, b-actin has a mean expression that is 360-fold higher than HPRT (8.5 PCR cycles difference to reach fluorescence threshold) and shows a good coefficient of correlation for colorectal and skin tissues (Table 2). Alternatively, when the gene of interest has intermediate or low expression, or in all situations when sufficient cells are present, a low-copy housekeeping gene controls best for RNA isolation efficiency, RNA quality, and RT-efficiency. HPRT then remains a good choice, since its expression is relatively low. In a recent study by Vandesompele et al,28 a similar method was presented using the geometric mean of the expression of 10 housekeeping genes. The design of their method was not suitable to identify a single normalization gene, but rather combinations of at least two genes that best represented the geometric mean. Even though most of their samples were derived from cell lines (neuroblastoma) and cell cultures (fibroblasts), hence very different in nature from our epithelial tissue samples, HPRT was frequently present in the best combination of genes. In addition, b-2-microglobulin was one of the worst scoring genes in their samples. This also corroborated our data. Our study clearly demonstrates that the expression of the HPRT gene accurately reflects the mean expression of multiple commonly used housekeeping genes. For differential expression studies in cancer research, HPRT will be the most economic and accurate choice as single normalization gene. Acknowledgements We thank Dr van Muijen and Dr Verhaegh of the Departments of Pathology and Experimental Normalization of gene expression JB de Kok et al Urology, respectively (UMC Nijmegen, The Netherlands), for the donation of clinical samples from their collection. 14 References 15 1 de Kok JB, Verhaegh GW, Roelofs RW, et al. DD3(PCA3), a very sensitive and specific marker to detect prostate tumors. Cancer Res 2002;62:2695–2698. 2 Span PN, Manders P, Heuvel JJ, et al. Expression of the transcription factor Ets-1 is an independent prognostic marker for relapse-free survival in breast cancer. Oncogene 2002;21:8506–8509. 3 Goidin D, Mamessier A, Staquet MJ, et al. Ribosomal 18S RNA prevails over glyceraldehyde-3-phosphate dehydrogenase and beta-actin genes as internal standard for quantitative comparison of mRNA levels in invasive and noninvasive human melanoma cell subpopulations. Anal Biochem 2001;295:17–21. 4 Kunth K, Hofler H, Atkinson MJ. Quantification of messenger RNA expression in tumors: which standard should be used for best RNA normalization? Verh Dtsch Ges Pathol 1994;78:226–230. 5 Zhong H, Simons JW. Direct comparison of GAPDH, beta-actin, cyclophilin, and 28S rRNA as internal standards for quantifying RNA levels under hypoxia. Biochem Biophys Res Commun 1999;259:523–526. 6 Thellin O, Zorzi W, Lakaye B, et al. Housekeeping genes as internal standards: use and limits. J Biotechnol 1999;75:291–295. 7 Schmittgen TD, Zakrajsek BA. Effect of experimental treatment on housekeeping gene expression: validation by real-time, quantitative RT-PCR. J Biochem Biophys Methods 2000;46:69–81. 8 Wu YY, Rees JL. Variation in epidermal housekeeping gene expression in different pathological states. Acta Derm Venereol 2000;80:2–3. 9 Gerard CJ, Andrejka LM, Macina RA. Mitochondrial ATP synthase 6 as an endogenous control in the quantitative RT-PCR analysis of clinical cancer samples. Mol Diagn 2000;5:39–46. 10 Finnegan MC, Goepel JR, Hancock BW, et al. Investigation of the expression of housekeeping genes in nonHodgkin’s lymphoma. Leuk Lymphoma 1993;10: 387–393. 11 Degan M, Mazzocco FT, Di Francia R, et al. Normalizing complementary DNA by quantitative reverse transcriptase-polymerase chain reaction of beta 2microglobulin: molecular monitoring of minimal residual disease in acute promyelocytic leukemia. Diagn Mol Pathol 2000;9:98–109. 12 Jass JR, Sobin LH (eds). Histological typing of intestinal tumours. In: International Histological Classification of Tumours, 2nd edn. Springer-Verlag: New York, 1989. 13 Sobin LH Wittekind CH (eds). International Union Against Cancer (UICC) TNM Classification of 16 17 18 19 20 21 22 23 24 25 26 27 28 Malignant Tumours, 5th edn. Wiley-Liss: New York, 1997. Spiessl B, Beahrs OH, Hermanek P, et al. TNM Atlas: Illustrated Guide to the TNM/pTNM Classification of Malignant Tumors, 3rd edn, 2nd revision. SpringerVerlag: Berlin, 1992, pp 274–279. Mostofi FK, Sobin LH, Torloni H. Histologic Typing of Urinary Bladders, Vol. 10. World Health Organisation: Geneva, 1973. de Kok JB, Schalken JA, Aalders TW, et al. Quantitative measurement of telomerase reverse transcriptase (hTERT) mRNA in urothelial cell carcinomas. Int J Cancer 2000;87:217–220. User Manual for TaqMan Human Endogenous Control Plate. Perkin Elmer Corporation: Eden Prairie, 1999. Chretien S, Dubart A, Beaupain D, et al. Alternative transcription and splicing of the human porphobilinogen deaminase gene result either in tissue-specific or in housekeeping expression. Proc Natl Acad Sci USA 1988;85:6–10. Mensink E, van de Locht A, Schattenberg A, et al. Quantitation of minimal residual disease in Philadelphia chromosome positive chronic myeloid leukaemia patients using real-time quantitative RT-PCR. Br J Haematol 1998;102:768–774. Gibson UE, Heid CA, Williams PM. A novel method for real time quantitative RT-PCR. Genome Res 1996;6: 995–1001. Raychaudhuri S, Stuart JM, Altman RB. Principal components analysis to summarize microarray experiments: application to sporulation time series. Pac Symp Biocomput 2000; 455–466. Jolliffe IT. Principal Component Analysis. Springer: New York, 2001. Savonet V, Maenhaut C, Miot F, et al. Pitfalls in the use of several ‘housekeeping’ genes as standards for quantitation of mRNA: the example of thyroid cells. Anal Biochem 1997;247:165–167. Tricarico C, Pinzani P, Bianchi S, et al. Quantitative real-time reverse transcription polymerase chain reaction: normalization to rRNA or single housekeeping genes is inappropriate for human tissue biopsies. Anal Biochem 2002;309:293–300. el Naggar AK, Batsakis JG, Teague K, et al. Acridine orange flow cytometric analysis of renal cell carcinoma. Clinicopathologic implications of RNA content. Am J Pathol 1990;137:275–280. Enker WE, Kimmel M, Cibas ES, et al. DNA/RNA content and proliferative fractions of colorectal carcinomas: a five-year prospective study relating flow cytometry to survival. J Natl Cancer Inst 1991;83: 701–707. Freeman WM, Walker SJ, Vrana KE. Quantitative RTPCR: pitfalls and potential. Biotechniques 1999;26: 112–115. Vandesompele J, De Preter K, Pattyn F, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 2002;3:1–12. 159 Laboratory Investigation (2005) 85, 154–159