|
|
||||||||
Affymetrix, Inc., Santa Clara, California 95051
| ABSTRACT |
|---|
|
|
|---|
housekeeping genes; DNA chips; human gene expression; cellular maintenance genes
| INTRODUCTION |
|---|
|
|
|---|
| MATERIALS AND METHODS |
|---|
|
|
|---|
Fragmentation, array hybridization, and scanning.
Labeled target was fragmented by incubation at 94°C for 35 min in the presence of 40 mM Tris-acetate, pH 8.1, 100 mM potassium acetate, and 30 mM magnesium acetate. The hybridization solution consisted of 20 µg fragmented cRNA and 0.1 mg/ml sonicated herring sperm DNA, in 1x MES buffer (containing 100 mM MES, 1 M Na+, 20 mM EDTA, and 0.01% Tween 20). The hybridization mixture was heated to 99°C for 5 min followed by incubation at 45°C for 5 min before injection of the sample into the probe array cartridge. All preparations and hybridizations were performed in duplicate and were carried out at 45°C for 1617 h with mixing on a rotisserie at 60 rpm. Following hybridization, the solutions were removed, arrays were rinsed with 1x MES. Subsequent washing and staining of the arrays was carried out using the GeneChip Fluidics station protocol EukGE_WS2. The EukGE_WS2 protocol included two posthybridization washes, staining, and a poststain wash. The first wash consisted of 10 cycles of 2 mixes per cycle with nonstringent wash buffer (6x SSPE, 0.01% Tween 20, and 0.005% antifoam) at 25°C. The second wash consisted of 4 cycles of 15 mixes per cycle with stringent wash buffer (100 mM MES, 0.1 M Na+, and 0.01% Tween 20) at 50°C. The probe arrays were stained for 10 min in streptavidin-phycoerythrin solution (SAPE) [1x MES solution, 0.005% antifoam, 10 µg/ml SAPE (Molecular Probes, Eugene, OR), and 2 µg/µl acetylated BSA (Sigma, St. Louis, MO)] at 25°C. The poststain wash consisted of 10 cycles of 4 mixes per cycle at 25°C. The probe arrays were treated for 10 min in antibody solution [1x MES solution, 0.005% antifoam, 2 µg/µl acetylated BSA, 0.1 µg/µl normal goat IgG (Sigma Chemical), 3 µg/µl antibody (goat), and antistreptavidin, biotinylated (Vector Laboratories, Burlingame, CA)] at 25°C. The final wash consisted of 15 cycles of 4 mixes per cycle at 30°C. Following washing and staining, probe arrays were scanned twice (multiple image scan) at 3-µm resolution using the GeneChip System confocal scanner made for Affymetrix by Hewlett-Packard.
Probe arrays.
The arrays were synthesized using light-directed combinatorial chemistry as described previously (9, 10). The HuGeneFL GeneChip probe arrays used for the current study contain probe sets representing 7,129 genes. The oligonucleotides are 25 bases in length. Probes are complementary and correspond to human genes registered in Unigene, GenBank, and The Institute for Genomic Research Database (TIGR). Each probe set has oligonucleotides that are identical to sequence in the gene and oligonucleotides that contain a homomeric (base transversion) mismatch at the central base position of the oligomer used for measuring cross hybridization. Probes are selected with a bias toward the 3' region of each gene. Probe pairs representing human genes such as glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ß-actin, transferrin receptor, and transcription factor ISGF-3 serve as internal controls for monitoring RNA integrity. In addition, the probe arrays contain oligonucleotides representing sequences of bacterial genes, BioB, BioC, and BioD, and one phage gene, Cre, as quantitative standards. Copy numbers are determined by correlating the known concentrations of the spiked standards with their hybridization intensities as described previously (16). Copies per cell are calculated based on the assumption that the average transcript length is 1 kb and there are 300,000 transcripts per cell.
Analysis.
All samples were prepared and hybridized in duplicate. Only those transcripts detected as present in duplicate hybridizations or absent in duplicate hybridizations are reported. Of the transcripts present in duplicate hybridizations, the hybridization values were within twofold. The values from the duplicate hybridizations were averaged. GeneChip 3.0 software was used to scan and analyze the data. Microsoft Excel and Microsoft Access were also used for data analysis.
| RESULTS |
|---|
|
|
|---|
(GenBank accession no. J04617) is detected at high abundance and the transcript for E2 ubiquitin (U39317) is detected in low abundance. Of the 535 genes, 288 of the transcripts vary in expression level by 5- to 10-fold, 134 transcripts vary by 11- to 19-fold, and 69 vary by greater than 19-fold. (For a list of the 535 genes expressed in all 11 tissues and the 47 transcripts expressed at the same levels, see Tables 1 and 2 of the supplementary material.1)
The majority of the maintenance transcripts detected were present in moderate levels.
The distributions of transcripts detected in all 11 tissues sorted by tissue type and abundance level are shown in Fig. 1. The subset of transcripts expressed in each of the tissues, the maintenance transcripts, sorted by tissue type and abundance level are shown in Fig. 2. Most transcripts detected in any one cell type are detected at low levels,
5 copies per cell. The majority of the maintenance transcripts detected are at moderate levels, 1050 copies per cell. This may suggest that most maintenance transcripts are produced in excess and regulation occurs during translation or protein modification and/or delivery. Alternatively, because we chose to study whole organs, in which the variety of transcripts produced could be substantially complex, we may not be detecting all of the low-abundance messages.
|
|
In a study of only adult tissues, 695 transcripts are expressed in all 7 tissues, with a subset of 241 genes expressed at the same level; 333 of the genes vary in expression level by 5- to 10-fold. Forty genes expressed in all 7 tissues differ in transcript levels by greater than 19-fold, and of these, 8 differ by more than 50-fold, including COX7A muscle isoform (GenBank accession no. M83186) varying by 52-fold (highest in heart, lowest in kidney, pancreas, and testis), lectin (J04456) varying by 58-fold (highest in uterus, lowest in kidney and pancreas), myosin heavy chain (AF001548) varying by 61-fold (highest in uterus, lowest in brain and pancreas), elongation factor-1
(Z21507) varying by 69-fold (highest in pancreas, lowest in lung and kidney), RNA polymerase II elongation protein (Z47087) varying by 70-fold (highest in brain, lowest in pancreas), extracellular mRNA for glutathione peroxidase (D00632) varying by 78-fold (highest in kidney, lowest in brain, pancreas, and testis), 14-9-9 protein
-chain (D78577) varying by 81-fold (highest in brain, lowest in testis), and L-arginine:glycine amidinotransferase (S68805) varying by 133-fold (highest in pancreas, lowest in heart and lung). A set of genes frequently used as controls in standard expression analysis were found to vary in expression level by 7- to 23-fold; these include ß-actin (M10277) varying by 7-fold, with highest expression in brain and uterus and lowest expression in heart, and GAPDH (M33197) varying by 8-fold, with highest expression in brain, heart, and kidney and lowest in pancreas. Another form of ß-actin (X00351) varies by 22-fold, with highest expression in uterus and lowest in pancreas.
-Actin (X13839) varies by 23-fold, and
-actin (M19283) varies by 9-fold.
In fetal tissues, we found 767 transcripts expressed in all four tissues, 397 of which are expressed at the same level, 310 vary in expression level by 5- to 10-fold, 45 vary by 11- to 19-fold, and 15 vary by more than 19-fold. (See Tables 3 and 4 of the supplementary material for a list of the 695 shared adult transcripts and the 768 shared fetal transcripts.)
Tissue-specific transcripts.
In the same experiments, we identified genes expressed uniquely in each of the tissues. For instance, in adult heart there were 3 transcripts not detected in the other 10 tissues, muscle glycogen synthase (GenBank accession no. J04501), MLC-1V/Sb isoform (M24248), and cytokine inducible nuclear protein (X83703). Not surprisingly, we found the greatest number of uniquely expressed genes in fetal tissues. Transcript numbers for brain were lower than anticipated, probably due to the difficulty of obtaining nondegraded brain RNA and the complexity of whole brain tissue. (See Table 5 of the supplementary material for a list of the genes uniquely expressed in a comparison of the tissues.)
| DISCUSSION |
|---|
|
|
|---|
Determining biologically relevant differences in expression levels.
What is a biologically relevant difference in expression level? From a functional perspective, protein activity is the most critical measure of biological significance. We know that biological systems are complex and regulation of gene expression is quite removed from the ultimate destiny of the gene product, protein activity. Transcription, posttranscriptional modification, translation, posttranslational modification, and transport processes are not 100% efficient (1, 4, 11, 12, 20, 25, 29, 33). The cell must be able to tolerate and compensate for processing inefficiencies. The system must be flexible and, in most cases, probably produce an excess of transcript. Our data suggest that this may be the case for the maintenance genes. In a comparison of the abundance levels of all of the transcripts detected in all of the tissues with the abundance levels of the maintenance transcripts alone, we found that the majority of transcripts are expressed in low abundance, less than five copies per cell, whereas the maintenance transcripts are present in moderate levels, 1050 copies per cell. Of course, some genes must be tightly controlled at the transcription step, but for the group of proteins responsible for basic cellular maintenance and survival, tight regulation at the transcription level is probably too risky. Studies of Saccharomyces cerevisiae and S. pombe support this line of reasoning. Twenty percent of genes in S. cerevisiae show noisy oscillations throughout the cell cycle, and in S. pombe it has been demonstrated that transcription is present in the absence of cell cycle progression and cellular concentrations of transcripts vary by two- to fourfold (3, 15). Here, we report genes as expressed at the same level if they are expressed in all 11 tissues at levels within fourfold. For most genes, differences less than fourfold are probably not biologically significant, but there is not sufficient data to conclude that a five- or sixfold difference is more biologically significant than a three- or fourfold difference (5, 14).
Until recently the technical challenge of accurately measuring small differences in gene expression has been practically insurmountable; consequently, there is little evidence to support the importance of small differences. For a subset of genes, it is likely that small differences have biological relevance, such as the genes encoding proteins that function differently when bound to high-affinity vs. low-affinity receptors or gene products triggering cellular cascades (2, 7, 18, 26). What is a biologically significant fold difference at the mRNA level? With so few data, it is difficult to know what a biologically significant difference in expression level is, but with the increase in sensitivity made possible by array technology and the development of other competing methods, we are surely about to find out.
| ACKNOWLEDGMENTS |
|---|
GeneChip is a registered trademark of Affymetrix, Inc.
Address for reprint requests and other correspondence: J. A. Warrington, Affymetrix, Inc., 3380 Central Expressway, Santa Clara, CA 95051 (E-mail: janet_warrington{at}affymetrix.com).
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D. Levens How the c-myc Promoter Works and Why It Sometimes Does Not J Natl Cancer Inst Monographs, July 1, 2008; 2008(39): 41 - 43. [Abstract] [Full Text] [PDF] |
||||
![]() |
F Salway, P. Day, W. Ollier, and T. Peakman Levels of 5' RNA tags in plasma and buffy coat from EDTA blood increase with time Int. J. Epidemiol., April 1, 2008; 37(suppl_1): i11 - i15. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Appolloni, F. Calzolari, G. Corte, R. Perris, and P. Malatesta Six3 Controls the Neural Progenitor Status in the Murine CNS Cereb Cortex, March 1, 2008; 18(3): 553 - 562. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Kidd, B. Nadler, S. Mane, G. Eick, M. Malfertheiner, M. Champaneria, R. Pfragner, and I. Modlin GeneChip, geNorm, and gastrointestinal tumors: novel reference genes for real-time PCR Physiol Genomics, August 20, 2007; 30(3): 363 - 370. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Greenwood, D. J. Jagger, L.-C. Huang, N. Hoya, P. R. Thorne, S. S. Wildman, B. F. King, K. Pak, A. F. Ryan, and G. D. Housley P2X receptor signaling inhibits BDNF-mediated spiral ganglion neuron development in the neonatal rat cochlea Development, April 1, 2007; 134(7): 1407 - 1417. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Ivanga, Y. Labrie, E. Calvo, P. Belleau, C. Martel, V. Luu-The, J. Morissette, F. Labrie, and F. Durocher Temporal analysis of E2 transcriptional induction of PTP and MKP and downregulation of IGF-I pathway key components in the mouse uterus Physiol Genomics, March 14, 2007; 29(1): 13 - 23. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Synnergren, T. L. Giesler, S. Adak, R. Tandon, K. Noaksson, A. Lindahl, P. Nilsson, D. Nelson, B. Olsson, M. C.O. Englund, et al. Differentiating Human Embryonic Stem Cells Express a Unique Housekeeping Gene Signature Stem Cells, February 1, 2007; 25(2): 473 - 480. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. C. Hsu, J. Chang, T. Wang, E. Steingrimsson, M. K. Magnusson, and K. Bergsteinsdottir Statistically designing microarrays and microarray experiments to enhance sensitivity and specificity Brief Bioinform, January 1, 2007; 8(1): 22 - 31. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Kullberg, M. A. Nilsson, U. Arnason, E. H. Harley, and A. Janke Housekeeping Genes for Phylogenetic Analysis of Eutherian Relationships Mol. Biol. Evol., August 1, 2006; 23(8): 1493 - 1503. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Borriello, R. Capparelli, M. Bianco, D. Fenizia, F. Alfano, F. Capuano, D. Ercolini, A. Parisi, S. Roperto, and D. Iannelli Genetic Resistance to Brucella abortus in the Water Buffalo (Bubalus bubalis) Infect. Immun., April 1, 2006; 74(4): 2115 - 2120. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Nishida, M. Yoshioka, and J. St-Amand Regulation of hypothalamic gene expression by glucocorticoid: implications for energy homeostasis Physiol Genomics, March 13, 2006; 25(1): 96 - 104. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Ishii, A. M. Wallace, X. Zhang, J. Gosselink, R. T. Abboud, J. C. English, P. D. Pare, and A. J. Sandford Stability of housekeeping genes in alveolar macrophages from COPD patients Eur. Respir. J., February 1, 2006; 27(2): 300 - 306. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Luu-The, P. Tremblay, and F. Labrie Characterization of Type 12 17{beta}-Hydroxysteroid Dehydrogenase, an Isoform of Type 3 17{beta}-Hydroxysteroid Dehydrogenase Responsible for Estradiol Formation in Women Mol. Endocrinol., February 1, 2006; 20(2): 437 - 443. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Chen, D. A. Rider, and R. Ruan Identification of Valid Housekeeping Genes and Antioxidant Enzyme Gene Expression Change in the Aging Rat Liver J. Gerontol. A Biol. Sci. Med. Sci., January 1, 2006; 61(1): 20 - 27. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Chen, Y. Wen, Z. Zhang, Y. Guo, J. A. Warrington, and M. L. Polan Microarray analysis of differentially expressed genes in vaginal tissues from women with stress urinary incontinence compared with asymptomatic women Hum. Reprod., January 1, 2006; 21(1): 22 - 29. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. W. Liu, S. T. Chen, and H. P. Liu Choice of endogenous control for gene expression in nonsmall cell lung cancer Eur. Respir. J., December 1, 2005; 26(6): 1002 - 1008. [Abstract] [Full Text] [PDF] |
||||
![]() |
C.L. Tower, S.L. Chappell, K. Morgan, N. Kalsheker, P.N. Baker, and L.J. Morgan Transforming growth factor {beta}1 regulates angiotensin II type I receptor gene expression in the extravillous trophoblast cell line SGHPL-4 Mol. Hum. Reprod., December 1, 2005; 11(12): 847 - 852. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Nicot, J.-F. Hausman, L. Hoffmann, and D. Evers Housekeeping gene selection for real-time RT-PCR normalization in potato during biotic and abiotic stress J. Exp. Bot., November 1, 2005; 56(421): 2907 - 2914. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Buzzeo, S. Enkemann, R. Nimmanapalli, M. Alsina, M. G. Lichtenheld, W. S. Dalton, and D. M. Beaupre Characterization of a R115777-Resistant Human Multiple Myeloma Cell Line with Cross-Resistance to PS-341 Clin. Cancer Res., August 15, 2005; 11(16): 6057 - 6064. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. de la Grange, M. Dutertre, N. Martin, and D. Auboeuf FAST DB: a website resource for the study of the expression regulation of human gene products Nucleic Acids Res., July 28, 2005; 33(13): 4276 - 4284. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. V. Jongeneel, M. Delorenzi, C. Iseli, D. Zhou, C. D. Haudenschild, I. Khrebtukova, D. Kuznetsov, B. J. Stevenson, R. L. Strausberg, A. J.G. Simpson, et al. An atlas of human gene expression from massively parallel signature sequencing (MPSS) Genome Res., July 1, 2005; 15(7): 1007 - 1014. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. D. Barber, D. W. Harmer, R. A. Coleman, and B. J. Clark GAPDH as a housekeeping gene: analysis of GAPDH mRNA expression in a panel of 72 human tissues Physiol Genomics, May 11, 2005; 21(3): 389 - 395. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. J. Wild, R. C. Krieg, J. Seidl, R. Stoehr, K. Reher, C. Hofmann, J. Louhelainen, A. Rosenthal, A. Hartmann, C. Pilarsky, et al. RNA expression profiling of normal and tumor cells following photodynamic therapy with 5-aminolevulinic acid-induced protoporphyrin IX in vitro Mol. Cancer Ther., April 1, 2005; 4(4): 516 - 528. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Helm, S. A. Enkemann, D. Coppola, J. S. Barthel, S. T. Kelley, and T. J. Yeatman Dedifferentiation Precedes Invasion in the Progression from Barrett's Metaplasia to Esophageal Adenocarcinoma Clin. Cancer Res., April 1, 2005; 11(7): 2478 - 2485. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Yanai, H. Benjamin, M. Shmoish, V. Chalifa-Caspi, M. Shklar, R. Ophir, A. Bar-Even, S. Horn-Saban, M. Safran, E. Domany, et al. Genome-wide midrange transcription profiles reveal expression level relationships in human tissue specification Bioinformatics, March 1, 2005; 21(5): 650 - 659. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Pelletier, V. Luu-The, S. Li, J. Ouellet, and F. Labrie Cellular Localization of mRNA Expression of Enzymes Involved in the Formation and Inactivation of Hormonal Steroids in the Mouse Prostate J. Histochem. Cytochem., October 1, 2004; 52(10): 1351 - 1356. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. B. Plotkin, H. Robins, and A. J. Levine Tissue-specific codon usage and the expression of human genes PNAS, August 24, 2004; 101(34): 12588 - 12591. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. L. Andersen, J. L. Jensen, and T. F. Orntoft Normalization of Real-Time Quantitative Reverse Transcription-PCR Data: A Model-Based Variance Estimation Approach to Identify Genes Suited for Normalization, Applied to Bladder and Colon Cancer Data Sets Cancer Res., August 1, 2004; 64(15): 5245 - 5250. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Zhang, A. Moseley, A. G. Jegga, A. Gupta, D. P. Witte, M. Sartor, M. Medvedovic, S. S. Williams, C. Ley-Ebert, L. M. Coolen, et al. Neural system-enriched gene expression: relationship to biological pathways and neurological diseases Physiol Genomics, July 8, 2004; 18(2): 167 - 183. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. T. Clark, M. S. Bodnar, M. Fox, R. T. Rodriquez, M. J. Abeyta, M. T. Firpo, and R. A. R. Pera Spontaneous differentiation of germ cells from human embryonic stem cells in vitro Hum. Mol. Genet., April 1, 2004; 13(7): 727 - 739. [Abstract] [Full Text] [PDF] |
||||
![]() |
E.-S. Han, Y. Wu, R. McCarter, J. F. Nelson, A. Richardson, and S. G. Hilsenbeck Reproducibility, Sources of Variability, Pooling, and Sample Size: Important Considerations for the Design of High-Density Oligonucleotide Array Experiments J. Gerontol. A Biol. Sci. Med. Sci., April 1, 2004; 59(4): B306 - B315. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. Abeyta, A. T. Clark, R. T. Rodriguez, M. S. Bodnar, R. A. R. Pera, and M. T. Firpo Unique gene expression signatures of independently-derived human embryonic stem cell lines Hum. Mol. Genet., March 15, 2004; 13(6): 601 - 608. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. R. Klevecz, J. Bolen, G. Forrest, and D. B. Murray From the Cover: A genomewide oscillation in transcription gates DNA replication and cell cycle PNAS, February 3, 2004; 101(5): 1200 - 1205. [Abstract] [Full Text] [PDF] |
||||
| < |