Control of gene expression during T cell activation: alternate regulation of mRNA transcription and mRNA stability
- Chris Cheadle†1Email author,
- Jinshui Fan†2,
- Yoon S Cho-Chung1,
- Thomas Werner3,
- Jill Ray4,
- Lana Do4,
- Myriam Gorospe2 and
- Kevin G Becker5
© Cheadle et al; licensee BioMed Central Ltd. 2005
Received: 18 February 2005
Accepted: 20 May 2005
Published: 20 May 2005
Microarray technology has become highly valuable for identifying complex global changes in gene expression patterns. The effective correlation of observed changes in gene expression with shared transcription regulatory elements remains difficult to demonstrate convincingly. One reason for this difficulty may result from the intricate convergence of both transcriptional and mRNA turnover events which, together, directly influence steady-state mRNA levels.
In order to investigate the relative contribution of gene transcription and changes in mRNA stability regulation to standard analyses of gene expression, we used two distinct microarray methods which individually measure nuclear gene transcription and changes in polyA mRNA gene expression. Gene expression profiles were obtained from both polyA mRNA (whole-cell) and nuclear run-on (newly transcribed) RNA across a time course of one hour following the activation of human Jurkat T cells with PMA plus ionomycin. Comparative analysis revealed that regulation of mRNA stability may account for as much as 50% of all measurements of changes in polyA mRNA in this system, as inferred by the absence of any corresponding regulation of nuclear gene transcription activity for these groups of genes. Genes which displayed dramatic elevations in both mRNA and nuclear run-on RNA were shown to be inhibited by Actinomycin D (ActD) pre-treatment of cells while large numbers of genes regulated only through altered mRNA turnover (both up and down) were ActD-resistant. Consistent patterns across the time course were observed for both transcribed and stability-regulated genes.
We propose that regulation of mRNA stability contributes significantly to the observed changes in gene expression in response to external stimuli, as measured by high throughput systems.
Virtually all microarray studies to-date have measured changes in steady-state mRNA levels by harvesting total cellular RNA and using it to generate probes through a variety of strategies including end-labeling of purified mRNA , incorporating a label into the first strand cDNA made from mRNA , or attaching a T7 RNA polymerase promoter during cDNA synthesis, then labeling of the resulting RNA . More recently, several groups have demonstrated the feasibility of hybridizing metabolically labeled mRNAs directly from nuclear run-on (NRO) reactions to nylon filter microarrays in order to investigate nascent transcripts [1, 4–6]. Schuhmacher et al. , in particular, used a B cell line carrying a conditional, tetracycline-regulated myc gene, and found that myc induction resulted in only a small overlap in regulated mRNAs at 4 hours post-induction when comparing polyA mRNA and NRO RNA on microarrays. This early work provided evidence that transcriptional activation of genes does not necessarily lead to a corresponding increase of their steady-state mRNA levels.
More recently, our laboratory has examined the relationship between newly transcribed (NRO) RNA and polyA mRNA in a stress model using human non-small lung carcinoma H1299 cells. In response to a variety of stresses (ultraviolet light, heat shock, or prostaglandin), we found that approximately half of the observed changes in mRNA levels of stress-regulated genes were accompanied by a corresponding increase or decrease in gene transcription as measured by NRO. The remaining half of stress-altered changes in gene expression was largely attributable to changes in mRNA turnover, thus suggesting that, on a global level, changes in mRNA turnover profoundly influence gene expression patterns .
Several questions, however, remained to be answered from these earlier studies. Since in both the myc induction and stress experiments, as mentioned above, measurements of both newly transcribed and polyA mRNAs were made at a single time point, there existed a reasonable possibility of temporal disjunctions between the timing of mRNA new gene synthesis and the rates of accumulation of mRNA in the cell. The second question remaining unanswered is whether or not highly significant levels of mRNA stability regulation (> 50% of all measured gene expression in the stress example) is common to different biological model systems. In order to begin to address these questions we investigated changes at the levels of transcription and total cellular mRNA abundance simultaneously across a time course of activation using Jurkat T cells.
T-cell activation is one of the most widely studied models of cellular response to exogenous stimulation. The initial events include rapid signaling via protein-protein interactions, phosphorylation/dephosphorylation of target signaling molecules, and release of Ca2+ from intracellular stores. Subsequent activation of signal transduction cascades culminates in the implementation of gene expression patterns characteristic of the immune response. Initial microarray studies using T cells have focused on gene expression changes occurring several hours after activation [7–13], even though earlier work using more traditional methods had defined the commitment period for T-cell activation, including alteration in gene expression patterns, as occurring between 1–2 h after exposure to the activating agent . In order to investigate these earlier gene expression changes we chose to examine a time course of activation spanning the first hour after stimulation. While a recent study by Garcia-Martinez et al.  in yeast using a similar approach demonstrated large shifts in mRNA stability following a glucose-to-galactose shift, the work presented here is the first systematic accounting of the changes in both gene transcription and mRNA stability in response to a major cellular activation event over a defined time period in higher eukaryotes.
The final, relatively minor groups of regulated genes included genes which were regulated at both the transcriptional and the polyA mRNA levels in either the same (8.4%) or opposite directions (9%). The relatively low concordance between transcriptional production of mRNA and its measured appearance in polyA mRNA levels was somewhat surprising, although clear examples of coordinated step-wise production were noted for some key genes, as for example, the early response genes EGR1 and ETR101 (previously shown to be induced at 30 minutes by phorbol ester treatment of a human promyelocytic leukemia cell line  and, the apoptosis-related genes DAP (death-associated protein, mediator of interferon-gamma-induced apoptosis) and CASP3, as well as the immune response signal transducer and activator of transcription (STAT) 6 (Fig. 3A).
Dramatic activation of immune response, immediate-early response genes, and apoptosis-related genes was observed in the nuclear run-on RNA as early as 5 minutes following activation (Fig. 3A). This group of genes included immediate-early response genes commonly up-regulated during cellular activation (ETR101, Myb, Myc, and genes of the JUN and EGR families), genes specifically associated with an early response in immune cells (IL6, IL8, STAT4, STAT6, PPP3C, NFKBIA, IRF5, and CD69), as well as genes involved in regulating apoptosis (BCL2A1, CASP3, CASP9, CASP10, and DAP). Many of these same genes were eventually up-regulated in polyA mRNA later in the time course and their increase in expression after one hour was independently validated by single end-point PCR validation using GS320 technology (Fig. 3B).
A particularly interesting example of the dichotomy between transcription and changes in polyA mRNA levels was seen in the production of the mRNAs encoding NFKB1 (NF-kappa B), a key mediator of the transcriptional control of genes involved in the immune response and acute phase reactions, and its inhibitor, NFKBIA (NF-kappa B inhibitor A). Both NFKB1 and NFKBIA have previously been shown by microarray analysis to be significantly induced in polyA mRNA between 3–4 hrs following phorbol or lectin activation of either Jurkat or human peripheral blood lymphocytes [7, 8, 11, 17]. As demonstrated here (Fig. 3C), the production of NFKB1 mRNA clearly increases between 30 minutes and one hour at the transcriptional level without a detectable corresponding increase in polyA mRNA during that time (subsequent PCR analysis did show some increase at the steady-state level between 0 and 60 minutes for the NFKB1 gene but this increase failed to meet the significance thresholds set for the microarray analysis). NFKBIA, on the other hand, is rapidly induced transcriptionally to a maximum level by 30 minutes, returning essentially to baseline within one hour. Meanwhile, NFKBIA steady-state levels can be seen to gradually rise across the first hour of the time course. Analysis of the dynamics of gene expression for NFKB1 and its inhibitor as deduced from conventional microarrays might suggest that by one hour NFKB1 production had not yet begun (contradicted by the NRO data here), and also that the mRNA for the inhibitor of NKFKB1 is steadily increasing (when, in fact, it is clear from NRO data that virtually all increases in the production of NFKBIA have concluded by one hour). These data provide a clear example of how information from nuclear run-on microarrays can enhance studies of gene activation and feedback mechanisms.
The absence of significant regulation at the transcriptional level during a time in which gene expression is strongly perturbed at the polyA mRNA level during T cell activation reveals that many genes are being regulated through changes in stability. The lack of detectable transcriptional regulation of large numbers of steady-state mRNA gene changes is particularly striking insofar as the NRO measurements are likely to be even more sensitive to changes in gene expression than polyA mRNA measurements since they are a direct measure of newly synthesized mRNA. This finding is consistent with the work of Raghavan et. al. [12, 20] who also noted the presence of numerous transcripts exhibiting stimulus-dependent changes in mRNA decay in human T lymphocytes treated for 3 hours with anti-CD3 and/or anti-CD28 antibodies. These experiments were carried out by arresting transcription with Actinomycin D (in the presence or absence of activation), and mRNA turnover rates were globally measured by applying polyA mRNA to microarrays. Similarly, in this current work, patterns of Actinomycin D resistant changes in gene expression among groups of genes significantly regulated only at the polyA mRNA level supports the observation that large numbers of genes are regulated primarily by stability during T cell activation.
The smaller, but still substantial, group of genes which displayed significant transcriptional regulation without any corresponding changes in polyA mRNA, may either be a reflection of the enhanced sensitivity of transcriptional detection or, perhaps, result from a persistent lag between changes in transcriptional output and their reflection in steady-state mRNA levels. Such a lag might arise from differences in scale of the absolute size of mRNA pools between newly transcribed and polyA mRNA. In this scenario, changes in the amounts of mRNA that are readily detectable by nuclear run-on may be too small to have an immediate detectable impact on the steady-state RNA pools, possibly for long periods of time. Another possibility is that the nascent mRNAs of some of these genes are so rapidly degraded that they never significantly impact on polyA mRNA levels at all.
As previously noted (Fig. 6), up-regulation of gene expression was a dominant trend for transcriptionally regulated genes while down-regulation was dramatically favored for genes regulated at the polyA mRNA level. From the standpoint of an overall cellular economy of gene expression, effective control at the level of transcription can be achieved primarily by turning on new gene transcription, while at the whole-cell level, rapid and effective regulation can be achieved by massive shifts in the stability of existing mRNA pools. In fact, the single most dramatic regulatory event experimentally observed was the rapid clearance of large pools of steady-state mRNA, presumably as result of a sudden and demanding change in cellular conditions.
Questions remain as to the extent that this type of regulation may occur under differing biological scenarios. Global regulation of mRNA stability has thus far been systematically studied under conditions of stress and cellular activation. It remains to be determined whether or not changes in mRNA stability are responsible for altering gene expression programs in response to other biological conditions. The existence of these regulatory paradigms may require a reevaluation of the common model of the control of gene expression which essentially invokes the turning on and off of gene transcription in order to explain changes in polyA mRNA levels.
Although stability regulation might reasonably be considered as a cellular measure to bridge the gap between the very rapid events of signal transduction and the longer term induction of cellular programs involving the coordinated transcription of batteries of new gene synthesis, the current data actually suggests otherwise. One of the goals in this work was to examine changes in gene expression due to altered transcription and those due to altered turnover across a time course allowing sufficient time for the resolution of time disparities between new mRNA synthesis and their appearance as measurable cellular pools. Sufficient time for effective transcription and even translation was clearly demonstrated by the example of NFKB1 mRNA induction followed by the transcription of its downstream targets. Patterns of stability regulated genes were, in this model system, non-random and surprisingly persistent throughout the entire time course suggesting that stability regulation was a major component in the control of gene expression for significant periods of time. The presence of highly active and global regulation of polyA mRNA levels by stability-altering mechanisms suggests a new significant level of regulation in the control of gene expression which spans wide phylogenetic distances from yeast  to humans, emphasizes a possible parsimonious role for new gene synthesis in response to changing cellular environments, and may help to explain some of the difficulties encountered in attempts to comprehensively correlate clustered changes in polyA mRNA with common promoter regulatory elements.
Cell culture and treatment
The human Jurkat E6-1 T cell line was cultured in RPMI-1640 medium (GIBCO-BRL, Gaithersburg, MD) supplemented with 10% FCS (Hyclone, Logan, UT) and antibiotics, and incubated at 37°C, 5% CO2. Cells were stimulated with 40 ng/ml phorbol 12-myristate 13-acetate (PMA) and 1 μM ionomycin (Sigma, St. Louis, MO) (P+I). Actinomycin D (A.G. Scientific, San Diego, CA) was added at 10 ug/ml 30 minutes prior to activation with PMA + I.
Preparation and purification of RNA
Isolation of nuclei and preparation of nuclear run-on (NRO) RNA were as described [6, 21] with some modifications. Briefly, ~108 cells per sample were pelleted, washed with PBS, resuspended in ice-cold lysis buffer (20 mM Tris-HCl [pH 7.4], 20 mM NaCl, 5 mM MgCl2, 0.25% [v/v] NP-40), and incubated on ice for 5 min. Nuclei were spun down, resuspended in storage buffer (50 mM Tris-HCl [pH 8.3], 5 mM MgCl2, 0.1 mM EDTA-NaOH [pH 8.0], 40% [v/v] glycerol) and stored at -80°C until use. Thawed nuclei (200-μl aliquots) were mixed with 200 μl of reaction buffer (5 mM Tris-HCl [pH 8.0], 0.15 mM KCl, 2.5 mM MgCl2, 2.5 mM DTT, 0.5 mM of each ATP, UTP, GTP) plus 500 μCi of [α-33P]UTP (3,000 Ci/mmol, 10 mCi/ml; ICN), and incubated for 30 min at 30°C with shaking. Samples were then incubated with DNase I (100U, RNase-free; Roche Diagnostics) for 20 min at 37°C, and with proteinase K (1 μg/μl) for 1 h at 37°C. Finally, nascent RNA was purified by Sephadex G-50 column filtration. For the preparation of polyA mRNA cellular RNA, resting or stimulated Jurkat cells were lysed in STAT60™ (TEL-TEST, INC., Friendswood, TX) as described . RNA concentration and quality were assessed spectrophotometrically and by agarose gel electrophoresis. NRO RNA was stored at -80°C until use.
PolyA mRNA samples were radiolabeled and hybridized as previously described . In brief, 5 μg of total mRNA for each sample was annealed, in 16 μl H2O, with 1 μg of 24-mer poly(dT) primer (Research Genetics, Alabama), by heating at 65°C for 10 min and cooling on ice for 2 min. The RT reaction was performed by adding 8 μl of 5X first strand RT buffer (Life Technologies, Rockville, MD), 4 μl of 20 mM dNTPs minus dCTP) (Pharmacia, Piscataway, NJ), 4 μl of 0.1 M DTT, 40 U of RNAseOUT (Life Technologies), 6 μl of 3000 Ci/mmol α-33P dCTP (ICN Biomedicals, Costa Mesa, CA) to the RNA/primer mixture to a final volume of 40 μl. Two μl (400 U) of Superscript II reverse transcriptase (Life Technologies) was then added, and the sample was incubated for 60 min at 42°C. The reaction was stopped by the addition of five μl of 0.5 M EDTA. The samples were incubated at 65°C for 30 min after addition of 10 μl of 0.1 M NaOH in order to hydrolyze and remove RNA. The samples were pH neutralized by the addition of 25 μl of 0.5 M Tris, pH 8.0, and purified using Bio-Rad 6 purification columns (Hercules, CA). An aliquout of each labled sample was quantitated by liquid scintillation counting using a Beckman LS 6500. The entire remaining sample was stored at -20°C until use.
Array construction and hybridization
Microarray construction and hybridization were previously described . Briefly, NIA Human Focused Arrays consisting of a set of 4600 spotted cDNAs, arrayed in duplicate, representing a set of 2742 non-redundant genes (enriched in genes involved in immune function and signal transduction), were printed on Nytran + Supercharge nylon membranes (Schleicher & Schuell, Keene, NH), and hybridized with [α-33P]dCTP-labeled cDNA or [α-33P]UTP RNA probes overnight at 50°C as previously described , protocols available at http://www.grc.nia.nih.gov/branches/rrb/dna/dna.htm. Hybridized arrays were rinsed in 2 X SSC and 0.1% SDS twice at 55°C followed by washes in 2 X SSC and 0.1% SDS at 55°C. Microarrays were exposed for 1–3 d and scanned using a PhosphorImager (Molecular Dynamics, Sunnyvale, CA) at a 50-μm pixel resolution. ArrayPro software (MediaCybernetics, Silver Spring, MD) was used to convert the hybridization signals into raw intensity values; data generated were transferred into Microsoft Excel for further analysis.
Single endpoint PCR was carried out on polyA mRNA for gene detection and relative quantitative comparison between baseline (0 time) and I hour activation (with PMA +I). In brief, GS320 libraries for both control and activated polyA mRNA were prepared as described , normalized using a panel of ribosomal protein genes, and specific genes were amplified in duplicate using predefined GS320 primers.
Array data analysis
RNA samples (typically n = 3 to n = 5) were prepared from multiple experiments, each of which consisted of a consecutive series of time points. Equal amounts of total RNA (5 μg) or NRO RNA (108 cell-equivalents) were used in each hybridization. Raw intensity data for each experiment was transformed to log10, then used for the calculation of Z scores as described  [see Additional file 1]. Significant changes in gene expression were calculated in the form of Z ratios and/or Z test values (), using Z score values in all calculations. Z ratios constitute a measure of the change in gene expression of a given gene from its baseline value (in this case – time 0), expressed in units of standard deviation from the average change of all genes for that comparison. Z ratios are a direct measure of the likelihood that an observed change is an outlier in an otherwise normal distribution and, as such, are independent from underlying intensity values. Since the contents of the population of nuclear run-on and polyA mRNA are different in both complexity and number, care was taken not to compare Z score normalized intensities directly. Comparisons between Z ratios, however, test for equivalence of significant changes between the transcriptional and steady-state changes in gene expression each relative to its own population. All gene expression changes were assessed through comparison with untreated cells (time 0). A Z ratio value of ± 1.50 and/or a Z test value p < 0.0001 were the significance thresholds used in this study.
Hierarchical clustering was performed using the Cluster and TreeView software programs, developed at Stanford University . The clustering algorithm was set to complete linkage clustering using the uncentered Pearson correlation.
Mapping of transcription factor binding sites in selected genes was performed using software from Genomatix Software Gmbh, Munich, Germany http://www.genomatix.de. See Additional file 2 for a complete description of the analysis.
PMA plus Ionomycin
We are grateful to R.L. Wange (NIA, NIH) for kindly providing Jurkat cells.
- Legen J, Kemp S, Krause K, Profanter B, Herrmann RG, Maier RM: Comparative analysis of plastid transcription profiles of entire plastid chromosomes from tobacco attributed to wild-type and PEP- deficient transcription machineries. Plant J. 2002, 31: 171-188. 10.1046/j.1365-313X.2002.01349.x.PubMedView ArticleGoogle Scholar
- Schena M, Shalon D, Davis RW, Brown PO: Quantitative monitoring of gene expression patterns with a complementary DNA microarray. Science. 1995, 270: 467-470.PubMedView ArticleGoogle Scholar
- Eberwine J: Amplification of mRNA populations using aRNA generated from immobilized oligo(dT)-T7 primed cDNA. Biotechniques. 1996, 20: 584-591.PubMedGoogle Scholar
- Meininghaus M, Chapman RD, Horndasch M, Eick D: Conditional expression of RNA polymerase II in mammalian cells. Deletion of the carboxyl-terminal domain of the large subunit affects early steps in transcription. J Biol Chem. 2000, 275: 24375-24382. 10.1074/jbc.M001883200.PubMedView ArticleGoogle Scholar
- Schuhmacher M, Kohlhuber F, Holzel M, Kaiser C, Burtscher H, Jarsch M, Bornkamm GW, Laux G, Polack A, Weidle UH, Eick D: The transcriptional program of a human B cell line in response to Myc. Nucleic Acids Res. 2001, 29: 397-406. 10.1093/nar/29.2.397.PubMedPubMed CentralView ArticleGoogle Scholar
- Fan J, Yang X, Wang W, Wood WH, Becker KG, Gorospe M: Global analysis of stress-regulated mRNA turnover by using cDNA arrays. Proc Natl Acad Sci U S A. 2002, 99: 10611-10616. 10.1073/pnas.162212399.PubMedPubMed CentralView ArticleGoogle Scholar
- Ellisen LW, Palmer RE, Maki RG, Truong VB, Tamayo P, Oliner JD, Haber DA: Cascades of transcriptional induction during human lymphocyte activation. Eur J Cell Biol. 2001, 80: 321-328.PubMedView ArticleGoogle Scholar
- Cristillo AD, Bierer BE: Identification of novel targets of immunosuppressive agents by cDNA- based microarray analysis. J Biol Chem. 2002, 277: 4465-4476. 10.1074/jbc.M108598200.PubMedView ArticleGoogle Scholar
- Marrack P, Mitchell T, Hildeman D, Kedl R, Teague TK, Bender J, Rees W, Schaefer BC, Kappler J: Genomic-scale analysis of gene expression in resting and activated T cells. Curr Opin Immunol. 2000, 12: 206-209. 10.1016/S0952-7915(99)00075-8.PubMedView ArticleGoogle Scholar
- Rogge L, Bianchi E, Biffi M, Bono E, Chang SY, Alexander H, Santini C, Ferrari G, Sinigaglia L, Seiler M, Neeb M, Mous J, Sinigaglia F, Certa U: Transcript imaging of the development of human T helper cells using oligonucleotide arrays. Nat Genet. 2000, 25: 96-101. 10.1038/75671.PubMedView ArticleGoogle Scholar
- Feske S, Giltnane J, Dolmetsch R, Staudt LM, Rao A: Gene regulation mediated by calcium signals in T lymphocytes. Nat Immunol. 2001, 2: 316-324. 10.1038/86318.PubMedView ArticleGoogle Scholar
- Raghavan A, Ogilvie RL, Reilly C, Abelson ML, Raghavan S, Vasdewani J, Krathwohl M, Bohjanen PR: Genome-wide analysis of mRNA decay in resting and activated primary human T lymphocytes. Nucleic Acids Res. 2002, 30: 5529-5538. 10.1093/nar/gkf682.PubMedPubMed CentralView ArticleGoogle Scholar
- Grolleau A, Bowman J, Pradet-Balade B, Puravs E, Hanash S, Garcia-Sanz JA, Beretta L: Global and specific translational control by rapamycin in T cells uncovered by microarrays and proteomics. J Biol Chem. 2002, 277: 22175-22184. 10.1074/jbc.M202014200.PubMedView ArticleGoogle Scholar
- Crabtree GR: Contingent genetic regulatory events in T lymphocyte activation. Science. 1989, 243: 355-361.PubMedView ArticleGoogle Scholar
- Garcia-Martinez J, Aranda A, Perez-Ortin JE: Genomic run-on evaluates transcription rates for all yeast genes and identifies gene regulatory mechanisms. Mol Cell. 2004, 15: 303-313. 10.1016/j.molcel.2004.06.004.PubMedView ArticleGoogle Scholar
- Shimizu N, Ohta M, Fujiwara C, Sagara J, Mochizuki N, Oda T, Utiyama H: Expression of a novel immediate early gene during 12-O-tetradecanoylphorbol-13-acetate-induced macrophagic differentiation of HL-60 cells. J Biol Chem. 1991, 266: 12157-12161.PubMedGoogle Scholar
- Schena M, Shalon D, Heller R, Chai A, Brown PO, Davis RW: Parallel human genome analysis: microarray-based expression monitoring of 1000 genes. Proc Natl Acad Sci U S A. 1996, 93: 10614-10619. 10.1073/pnas.93.20.10614.PubMedPubMed CentralView ArticleGoogle Scholar
- Kizaki H, Tadakuma T, Odaka C, Muramatsu J, Ishimura Y: Activation of a suicide process of thymocytes through DNA fragmentation by calcium ionophores and phorbol esters. J Immunol. 1989, 143: 1790-1794.PubMedGoogle Scholar
- Ruiz-Ruiz MC, Oliver FJ, Izquierdo M, Lopez-Rivas A: Activation-induced apoptosis in Jurkat cells through a myc-independent mechanism. Mol Immunol. 1995, 32: 947-955. 10.1016/0161-5890(95)00059-N.PubMedView ArticleGoogle Scholar
- Raghavan A, Bohjanen PR: Microarray-based analyses of mRNA decay in the regulation of mammalian gene expression. Brief Funct Genomic Proteomic. 2004, 3: 112-124.PubMedView ArticleGoogle Scholar
- Gorospe M, Wang X, Holbrook NJ: p53-dependent elevation of p21Waf1 expression by UV light is mediated through mRNA stabilization and involves a vanadate-sensitive regulatory system. Mol Cell Biol. 1998, 18: 1400-1407.PubMedPubMed CentralView ArticleGoogle Scholar
- Cheadle C, Vawter M, Freed WJ, Becker KG: Analysis of microarray data using Z score transformation. Journal of Molecular Diagnostics. 2003, 5: 73-81.PubMedPubMed CentralView ArticleGoogle Scholar
- Barrett T, Cheadle C, Wood WB, Teichberg D, Donovan DM, Freed WJ, Becker KG, Vawter MP: Assembly and use of a broadly applicable neural cDNA microarray. Restor Neurol Neurosci. 2001, 18: 127-135.PubMedGoogle Scholar
- Wang A, Pierce A, Judson-Kremer K, Gaddis S, Aldaz CM, Johnson DG, MacLeod MC: Rapid analysis of gene expression (RAGE) facilitates universal expression profiling. Nucleic Acids Res. 1999, 27: 4609-4618. 10.1093/nar/27.23.4609.PubMedPubMed CentralView ArticleGoogle Scholar
- Nadon R, Woody E, Shi P, Rghei N, Hubschle H, Susko E, Ramm P: Statistical inference in array genomics. Microarrays for the Neurosciences. Edited by: Geschwind D and Gregg J. 2002, Cambridge, MIT PressGoogle Scholar
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.PubMedPubMed CentralView ArticleGoogle Scholar
- Dahlquist KD, Salomonis N, Vranizan K, Lawlor SC, Conklin BR: GenMAPP, a new tool for viewing and analyzing microarray data on biological pathways. Nat Genet. 2002, 31: 19-20. 10.1038/ng0502-19.PubMedView ArticleGoogle Scholar