Comprehensive quality control utilizing the prehybridization third-dye image leads to accurate gene expression measurements by cDNA microarrays
© Wang et al; licensee BioMed Central Ltd. 2006
Received: 04 November 2005
Accepted: 14 August 2006
Published: 14 August 2006
Gene expression profiling using microarrays has become an important genetic tool. Spotted arrays prepared in academic labs have the advantage of low cost and high design and content flexibility, but are often limited by their susceptibility to quality control (QC) issues. Previously, we have reported a novel 3-color microarray technology that enabled array fabrication QC. In this report we further investigated its advantage in spot-level data QC.
We found that inadequate amount of bound probes available for hybridization led to significant, gene-specific compression in ratio measurements, increased data variability, and printing pin dependent heterogeneities. The impact of such problems can be captured through the definition of quality scores, and efficiently controlled through quality-dependent filtering and normalization. We compared gene expression measurements derived using our data processing pipeline with the known input ratios of spiked in control clones, and with the measurements by quantitative real time RT-PCR. In each case, highly linear relationships (R2>0.94) were observed, with modest compression in the microarray measurements (correction factor<1.17).
Our microarray analytical and technical advancements enabled a better dissection of the sources of data variability and hence a more efficient QC. With that highly accurate gene expression measurements can be achieved using the cDNA microarray technology.
Microarray technology allows a comprehensive examination of gene expression profiles and the regulations of their changes at a whole genome level [1–3]. It has great potential in the study of complex human diseases . However, the technology is prone to noise and low reproducibility . Correlations with other platforms including RT-PCR [4, 5], and between different microarray platforms are often unsatisfactory [6–9]. On the other hand, many disease processes involve subtle gene perturbations that require highly accurate gene expression measurements. The noise in microarrays if not adequately reduced, can obscure the true biological variations and presents an obstacle for data-mining tools to distinguish biology from artifacts. For this reason rigorous QC standards are needed for the microarrays . This in turn requires a clear understanding of the sources of data variability, so that the contributing factors can be appropriately dissected and most efficiently controlled.
Because of the lack of consistent quality control (QC) standards, spotted arrays fabricated in academic laboratories are usually more susceptible to QC problems than commercial arrays [6–8]. Their advantages include much reduced fabrication cost and higher content flexibility than commercial arrays. For example, they can be designed to target specific genetic pathways, or to perform promoter analysis for genes of interest . Therefore developing a generalizable, efficient microarray QC scheme for spotted arrays is important.
We have previously reported a microarray image processing software Matarray, which specializes in quantitative QC of data acquisition [12, 13]. Using it we have shown that several major sources of data variability are readily identifiable from the post-hybridization image, including high or non-uniform noise profiles, low or saturated signal intensities, and irregular spot sizes and shapes. Their resultant effect on data reliability can be well characterized through the definition of a set of individual quality scores each measuring the impact of a corresponding factor, and a composite score q com , which gives an overall assessment of the data quality acquired from each spot on the array . Through numerous experiments we have demonstrated the advantages of utilizing the ratio-q com plot for data filtering and normalization [12, 13]. Nevertheless, there are sources of variability that cannot be directly or quantitatively evaluated from the post-hybridization image. One important example is the quality of array fabrication. The generation of microarray slides involves coating of the glass slides, printing up to tens of thousands of amplified cDNA or oligonucleotide "probes" and fixing/blocking of the slide. During this process, variable amounts of material can be deposited and/or retained on the activated glass surface depending a number of variables. When the amount of immobilized probe is inadequate the measurements made on such arrays can be unreliable [14–18]. Noise and artifacts introduced to the arrays at this stage will also directly affect the quality of hybridization. Until recently, such problems have been difficult to quantitatively evaluate and control for each and every array, since the array is typically "invisible" prior to hybridization [14, 16, 17]. To overcome this difficulty, we have made a significant technological advancement in microarray QC by conceiving and developing a three-color microarray platform [15–17], which we termed third dye array visualization (TDAV) technology . The approach labels the cDNA probes printed on the array slides with a cyanine dye-compatible third dye (TD) fluorescein , and makes prehybridization quantitative assessment of array quality possible, so that precious samples as well as laboratory and analytical efforts will not be wasted over poor-quality slides. In the last several years, the quantitative third dye threshold for slide QC has been extensively investigated by us [15–17].
In this report, we further investigate the advantage of TDAV in better dissecting the sources for data variability at the spot level, and develop a data filtering and normalization procedure that incorporates the information from the TD image. We utilize data from four different microarray experiments to validate our procedure. We evaluate the accuracy of our microarray measurements by comparing them with the known input ratios for spiked in control clones, and with the measurements by quantitative real time RT-PCR.
Results and discussion
TD intensity and gene-dependent bias in ratio measurements
Compression in microarray measurements. The compression depends on the fold-of-change, and exhibits complex characteristics when the spot TD intensity falls below 5,000 RFU/pixel. Data presented are from the spiked-in Arabidopsis in vitro transcripts in experiment1. At each input ratio, the mean compression factor (ratio between the measured and input ratio) is determined for all spots with TD intensity above 5,000 RFU/pixel, and a linear regression is performed between the compression and the spot TD intensity for all spots below 5,000 RFU/pixel, in the same fashion as described in figure 1
Input ratio (Cy5:Cy3)
Compression factor for spots > 5,000 RFU/pixel
Slope of linear regression (< 5,000 RFU/pixel)
0.78 (R2~0.82, p < 0.0001)
0.42 (R2~0.62, p < 0.0001)
0.27 (R2~0.46, p < 0.0001)
0.10 (R2~0.30, p < 0.03)
0.05 (R2~0.04, p < 0.06)
0.14 (R2~0.18, p < 0.0001)
0.27 (R2~0.77, p < 0.0001)
We have also examined the data variability in relation to the TD spot intensity. At each spiked in ratio, we determined the standard deviation (SD) and the coefficient of variation (CV) in ratio measurements for every data point using its 20 nearest neighbors on the ratio-TD intensity plot. We find that above 5,000 RFU/pixel both SD and CV are essentially constants. With deficient TD intensities below 5,000 RFU/pixel SD increases initially followed by a drop at very low TD intensities due to the severe data compression in ratio measurements, whilst CV increases monotonically with decreasing TD intensity. The CV results for the data points corresponding to the spiked in ratio of 30:1 are given in figure 1B. In a real experiment where the transcript abundance for different genes spans a wide range and their folds of change vary, gene-dependent artifacts in measurements will occur. These results revealed that inadequate probe amount is an important major source of data variability that could cause complex features in data compression.
TD intensity and the spatial-dependent bias
Incorporating TD information in data filtering and normalization
Results in the preceding sections suggest that the TD intensity is a major factor that causes spot-level data reliability. In addition, other artifacts on the TD image can also influence the accuracy of expression measurements from post-hybridization images, including noise, spot size and shape irregularities [16, 17]. Based on these observations we formulated a quality measure for every spot from the TD image by defining
qTD = qint *qcom, TD (1)
where qcom, TD is the composite TD image quality score, defined according to signal-to-noise ratio, spot size, and background levels and variation, similarly as given in the equation (7) of . qint is given by:
In Matarray the default threshold = 5,000 RFU/pixel.
Incorporating our new quality score definition for the TD image, we have developed the following filtering and normalization procedure:
3. Evaluate the log R - qcom plot and decide on a filtering threshold for qcom . The quality score for all spots below this threshold will be reset to qcom = 0, as described in figure 4A.
4. Perform a local q com -dependent LOWESS normalization for all data points with qcom > 0. The LOWESS fit for SD will also be determined, and the Z-score will be calculated for normalized log R every spot by: . Figure 4B gives the log ratio after normalization for the same data set given in figure 4A.
After all normalization a final quality score will be defined
Qf = q TD •q com (3)
Only spots with Qf > 0 will be retained for further data mining and modeling.
The efficiency of our filtering and normalization procedure
The accuracy of gene expression measurements
In figure 6B the microarray measurements are compared to the RT-PCR results for the 22 genes in the rat liver experiments, and an overall good agreement is observed. The RT-PCR platform is generally considered more quantitative and accurate than the microarrays . 7 (open circles) of the 22 genes were identified as poor-quality data points as their Qf = 0 on at least one array. Excluding these genes a highly linear relationship (R2 ~ 0.94, p < 0.001) existed for the remaining 15. If we were to calibrate microarray measurements using the RT-PCR results as a standard, we would have Rcorrected = (Rmeasured) q , with a correction factor . This is much better than a q ~1.88 previously reported . If the 7 poor-quality data points were to be included, the agreement between the two platforms drops to R2 ~ 0.90 with a correction factor q ~ 1.27. This result strongly indicates that with stringent, efficient QC protocols, cDNA microarrays are able to generate high quality, quantitative gene expression measurements comparable to that by RT-PCR.
The impact on the biology
How much impact on the biological interpretation can such improvement in data quality bring about? To answer this question we have performed ontological analysis for the DE genes in each experiments using OntoExpress  and EASE . We found no significant difference in the pure number of DE gene predictions between data processed by either Z-norm or MA-norm (p > 0.5). However, at ontology level, a general trend appeared suggestive of Z-norm being able to lead to more focused, local biological themes, usually with more significant p-values (data not shown). For example, in experiment 3 at 6 hr after drug treatment, the apoptosis progression has been established with at least 40% cells were apoptotic according to the Annexin V/PI double staining . EASE analysis of the DE genes after Z-norm predicted enhanced presence of genes belonging to regulation of cell cycle, cell proliferation/death, lysosome, lytic vacuole, nucleus, etc, most of which were closely related to apoptosis. Using DE genes predicted after MA-norm, about one third of these categories were not detected. The interpretation of ontological analysis is a complex issue, with many unresolved problems . For example, since most of the ontological categories are not independent, it is still an open question on how to best recap the findings and evaluate significance . Therefore a quantitative evaluation of our findings awaits further methodology development in the field of ontological analysis.
In this report we have shown that when the probe amount is inadequate, severe compression in gene expression measurements occur with complex, gene-specific characteristics. Likewise, the normal variation in the amount of probes printed and immobilized is a major source for printing pin-dependent artifacts in the data. Utilizing TDAV, these problems rooted from array fabrication can be effectively dissected from hybridization problems, and efficiently controlled through the definition of a quality score qTD. In addition, we have developed a comprehensive two-step data filtering and normalization procedure based on the log R - qTD and log R - qcom plots, which was found to be more efficient than the commonly used MA-LOWESS approach. By confirming our microarray data with the known input ratio of spiked in controls clones, and with RT-PCR, we demonstrated that acquiring accurate measurements using cDNA microarrays is achievable with our TDAV technology and our data QC procedure. Furthermore, in a recent study where we compared measurements from our cDNA microarrays with those from Affymetrix and Agilent oligonucleotide array platforms, we observed a high correlation among the three, with no significant differences in terms of data quality. Specifically, using ANOVA we have found that the different platforms did not cause more variation than the replicate hybridizations within each individual platform . This is in contrast to the recent reports of poor performance by academic arrays which did not use our scheme [6–8]. Thus we wish to emphasize that academic labs can still perform high quality, inexpensive microarray experiments with our platform.
Our algorithms, although initially developed for cDNA arrays, are potentially applicable to spotted oligonucleotide arrays as well. Recently we have developed the three-color oligonucleotide array platform by introducing a third-dye labeled universal tracking oligonucleotide into the printing buffer, thus the quality of array fabrication can be quantitatively evaluated through the measurements of the tracking oligonucleotide . A high quality microarray platform will allow lab investigators to focus on their biological questions instead of the technical issues of the data, and will allow statisticians and bioinformatics investigators to develop more powerful complex analysis approaches.
Microarray slide fabrication
All the microarrays used in this report were fabricated in house with TDAV technology using an OmniGrid arrayer equipped with a stealth print head with 32 SMP3 MicroQuill pins (Telechem International, Sunnyvale, CA). The University of Iowa rat library, consisting of 35,040 clones, was used as a source of probe for the cDNA microarrays, and was printed over two poly-L-lysine coated glass slides, each possessing 18,432 features. Control clones were also printed on each array including β-actin, GAPDH, and 9 Arabidopsis clones, as well as negative controls including PolyA and spotting buffer. Specifically, the Arabidopsis clones were printed in a 2 fold serial dilution from 200 ng/ul to 6.25 ng/ul (6 dilutions in each series), 4 replicate series on each slide. To ensure consistent prehybridization TD image collection, we have implemented a confocal laser scanner calibration method utilizing FluorIS (CLONDIAG, Jena, Germany), a non-bleaching, reusable, calibration/standardization tool . This enabled us to set up universal, quantitative QC rules according to TD signal information. All the arrays used in our microarray experiments have passed our QC as previously described [16, 17].
All microarray data analyses were performed using our in-house Matarray package [12, 13]. Information from the prehybridization TD image as well as the two post hybridization cyanine images were acquired. Data from four different microarray experiments are utilized in this report. They comprise: (1) Profiling of BioBreeding (BB) rat thymus. Gene expressions were compared between the thymus of diabetes prone DRlyp/lyp (referred to as DP in this report) and diabetic resistant DR+/+ (DR) rats at day 40, and between day 65 and day 40 DP rats. This analysis utilized 4 animal pairs for each comparison, and 4 replicate hybridizations were performed for each pair, with 2 hybridizations reverse labeled to control for dye bias, totaling 32 hybridizations. For the 16 hybridizations (8 each way of dye labeling) that compared day 65 and day 40 DP rats, the labeling reactions of total thymus RNA were spiked with 4 Arabidopsis in vitro transcripts (cellulose synthase, chlorophyll a/b binding protein, ribulose-1,5-bisphosphate and triosphosphate isomerase) at known input ratios (30:1, 10:1, 1:1, and 1:0, respectively), giving rise to a total of 1216 data points. These clones enabled an evaluation of microarray measurement accuracy through the comparison of measured output ratios to the known RNA input ratios. (2) Gene expression profiling of the kidney from an end-stage renal failure (ESRD) rat model. Three pairs of 22 week-old Fawn Hooded Hypertensive (FHH) rats and control August Copenhagen Irish (ACI) rats were compared. For each animal pair, 2 replicate hybridizations were performed, with 1 reverse labeled, totaling 6 hybridizations. (3) Time course profiling of apoptosis progression in pancreatic islet β cells. Cells from a rat β cell line RIN-m5F were treated with a protein kinase C inhibitor staurosporine  at a high dose of 1 μM, and a low dose of 1 nM for 2, 4, and 6 hours, and were compared for differential gene expressions. At each time point, 6 replicate hybridizations were performed, with 3 of them reverse labeled, totaling 18 hybridizations. Cell apoptosis status were confirmed using Annexin V/PI double staining method as described in . (4) Profiling and comparison of liver gene expressions from day 65 BB-DR and Wistar-Furth (WF) rats. In this experiment, 4 animals from each strain were sacrificed and equal amounts of purified total RNA from the animals of the same strain were pooled. The two pools were then compared in a total of 6 replicate hybridizations, with 3 of them reverse labeled.
Real time quantitative RT-PCR of rat liver samples
In the last experiment that profiled the BB-DR and WF rat liver, expression of 22 genes that were deemed of biological interest according to our previous experiments had also been measured by quantitative real time RT-PCR. Two sets of primers were designed for each locus with Oligo 6.66 (Molecular Biology Insights, Cascade, CO) and were blasted against the whole genome to ensure gene specificity. The primers were synthesized by Sigma Genosys (The Woodlands, TX). A Rotor-Gene 3000 (Corbett Research, Morelake, Australia), was used to conduct monoplex quantitative real-time RT-PCR as previously described . Relative gene expression data (fold of change) between samples was determined using the mathematical model described in .
third dye array visualization
coefficient of variation
locally weighted scatter plot smoothing
This work is supported in part by a NIH/NIBIB grant (EB001421) awarded to MJH and XW, a NIH/NIAID grant (U19-AI62627) awarded to MJH, and by a special fund from the Children's Hospital of Wisconsin Foundation (2201377). We thank Rhonda Geoffrey, Sanaa Muheisen for assisting the data preparation, and Pawjai Khampang for careful reading of the manuscript.
- Brown PO, Botstein D: Exploring the new world of the genome with DNA microarrays. Nat Genet 1999, 21(1 Suppl):33–37. 10.1038/4462View ArticlePubMedGoogle Scholar
- Segal E, Shapira M, Regev A, Pe'er D, Botstein D, Koller D, Friedman N: Module networks: identifying regulatory modules and their condition-specific regulators from gene expression data. Nat Genet 2003, 34(2):166–176.View ArticlePubMedGoogle Scholar
- Ideker T, Thorsson V, Ranish JA, Christmas R, Buhler J, Eng JK, Bumgarner R, Goodlett DR, Aebersold R, Hood L: Integrated genomic and proteomic analyses of a systematically perturbed metabolic network. Science 2001, 292(5518):929–934. 10.1126/science.292.5518.929View ArticlePubMedGoogle Scholar
- Miklos GL, Maleszka R: Microarray reality checks in the context of a complex disease. Nat Biotechnol 2004, 22(5):615–621. 10.1038/nbt965View ArticlePubMedGoogle Scholar
- Chuaqui RF, Bonner RF, Best CJ, Gillespie JW, Flaig MJ, Hewitt SM, Phillips JL, Krizman DB, Tangrea MA, Ahram M, Linehan WM, Knezevic V, Emmert-Buck MR: Post-analysis follow-up and validation of microarray experiments. Nat Genet 2002, 32(Suppl):509–514. 10.1038/ng1034View ArticlePubMedGoogle Scholar
- Bammler T, Beyer RP, Bhattacharya S, Boorman GA, Boyles A, Bradford BU, Bumgarner RE, Bushel PR, Chaturvedi K, Choi D, Cunningham ML, Deng S, Dressman HK, Fannin RD, Farin FM, Freedman JH, Fry RC, Harper A, Humble MC, Hurban P, Kavanagh TJ, Kaufmann WK, Kerr KF, Jing L, Lapidus JA, Lasarev MR, Li J, Li YJ, Lobenhofer EK, Lu X, et al.: Standardizing global gene expression analysis between laboratories and across platforms. Nat Methods 2005, 2(5):351–356. 10.1038/nmeth754View ArticlePubMedGoogle Scholar
- Irizarry RA, Warren D, Spencer F, Kim IF, Biswal S, Frank BC, Gabrielson E, Garcia JG, Geoghegan J, Germino G, Griffin C, Hilmer SC, Hoffman E, Jedlicka AE, Kawasaki E, Martinez-Murillo F, Morsberger L, Lee H, Petersen D, Quackenbush J, Scott A, Wilson M, Yang Y, Ye SQ, Yu W: Multiple-laboratory comparison of microarray platforms. Nat Methods 2005, 2(5):345–350. 10.1038/nmeth756View ArticlePubMedGoogle Scholar
- Larkin JE, Frank BC, Gavras H, Sultana R, Quackenbush J: Independence and reproducibility across microarray platforms. Nat Methods 2005, 2(5):337–344. 10.1038/nmeth757View ArticlePubMedGoogle Scholar
- Yauk CL, Berndt ML, Williams A, Douglas GR: Comprehensive comparison of six microarray technologies. Nucleic Acids Res 2004, 32(15):e124. 10.1093/nar/gnh123PubMed CentralView ArticlePubMedGoogle Scholar
- Johnson K, Lin S: QA/QC as a pressing need for microarray analysis: meeting report from CAMDA'02. Biotechniques 2003, (Suppl):62–63.Google Scholar
- Odom DT, Zizlsperger N, Gordon DB, Bell GW, Rinaldi NJ, Murray HL, Volkert TL, Schreiber J, Rolfe PA, Gifford DK, Fraenkel E, Bell GI, Young RA: Control of pancreas and liver gene expression by HNF transcription factors. Science 2004, 303(5662):1378–1381. 10.1126/science.1089769PubMed CentralView ArticlePubMedGoogle Scholar
- Wang X, Ghosh S, Guo S-W: Quantitative quality control in microarray image processing and data acquisition. Nucleic Acids Research 2001, 29: E75–82. 10.1093/nar/29.15.e75PubMed CentralView ArticlePubMedGoogle Scholar
- Wang X, Hessner MJ, Wu Y, Pati N, Ghosh S: Quantitative quality control in microarray experiments and the application in data filtering, normalization and false positive rate prediction. Bioinformatics 2003, 19: 1341–1347. 10.1093/bioinformatics/btg154View ArticlePubMedGoogle Scholar
- Yue H, Eastman PS, Wang BB, Minor J, Doctolero MH, Nuttall RL, Stack R, Becker JW, Montgomery JR, Vainer M, Johnston R: An evaluation of the performance of cDNA microarrays for detecting changes in global mRNA expression. Nucleic Acids Res 2001, 29(8):E41–41. 10.1093/nar/29.8.e41PubMed CentralView ArticlePubMedGoogle Scholar
- Hessner MJ, Meyer L, Tackes J, Muheisen S, Wang X: Immobilized support-bound probe and glass surface chemistry as variables in microarray fabrication. BMC Genomics 2004, 5: 53. 10.1186/1471-2164-5-53PubMed CentralView ArticlePubMedGoogle Scholar
- Hessner M, Wang X, Hulse K, Meyer L, Wu Y, Nye S, Guo S-W, Ghosh S: Three color cDNA microarrays: quantitative assessment through the use of Fluorescein-Labeled Probes. Nucl Acids Res 2003, 31: e14. 10.1093/nar/gng014PubMed CentralView ArticlePubMedGoogle Scholar
- Hessner MJ, Wang X, Khan S, Meyer L, Schlicht M, Tackes J, Datta M, Jacob HJ, Ghosh S: Use of a three-color cDNA microarray platform to measure and control support-bound probe for improved data quality and reproducibility. Nucl Acids Res 2003, 31: e60. 10.1093/nar/gng059PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Y, Wang X, Guo SW, Ghosh S: Conditions to ensure competitive hybridization in two-color microarray: a theoretical and experimental analysis. Biotechniques 2002, 32(6):1342–1346.PubMedGoogle Scholar
- Wang X, Jiang N, Feng X, Xie Y, Tonellato PJ, Ghosh S, Hessner MJ: A Novel Approach For High Quality Microarray Processing Using Third-Dye Array Visualization Technology. IEEE Trans NanoBioscience 2003, 2(4):193–201. 10.1109/TNB.2003.816233View ArticlePubMedGoogle Scholar
- Rajeevan MS, Vernon SD, Taysavang N, Unger ER: Validation of array-based gene expression profiles by real-time (kinetic) RT-PCR. J Mol Diagn 2001, 3(1):26–31.PubMed CentralView ArticlePubMedGoogle Scholar
- Cleveland WS, Devlin SJ: Locally weighted regression: an approach to regression analysis by local fitting. journal of the American statistical Association 1988, 83(403):596–610. 10.2307/2289282View ArticleGoogle Scholar
- Bacon D, Watts D: Estimating the transition between two intersecting straight lines. Biometrika 1971, 58: 525–535. 10.2307/2334387View ArticleGoogle Scholar
- Colantuoni C, Henry G, Zeger S, Pevsner J: Local mean normalization of microarray element signal intensities across an array surface: quality control and correction of spatially systematic artifacts. Biotechniques 2002, 32(6):1316–1320.PubMedGoogle Scholar
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res 2002, 30(4):e15. 10.1093/nar/30.4.e15PubMed CentralView ArticlePubMedGoogle Scholar
- Wilson DL, Buckley MJ, Helliwell CA, Wilson IW: New normalization methods for cDNA microarray data. Bioinformatics 2003, 19(11):1325–1332. 10.1093/bioinformatics/btg146View ArticlePubMedGoogle Scholar
- Yuen T, Wurmbach E, Pfeffer RL, Ebersole BJ, Sealfon SC: Accuracy and calibration of commercial oligonucleotide and custom cDNA microarrays. Nucleic Acids Res 2002, 30(10):e48. 10.1093/nar/30.10.e48PubMed CentralView ArticlePubMedGoogle Scholar
- Draghici S, Khatri P, Shah A, Tainsky MA: Assessing the functional bias of commercial microarrays using the onto-compare database. Biotechniques 2003, (Suppl):55–61.
- Hosack DA, Dennis G, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome Biology 2003, 4(6):p4. 10.1186/gb-2003-4-6-p4View ArticleGoogle Scholar
- Wang X, Becker FF, Gascoyne PR: Membrane dielectric changes indicate induced apoptosis in HL-60 cells more sensitively than surface phosphatidylserine expression or DNA fragmentation. Biochim Biophys Acta 2002, 1564(2):412–420. 10.1016/S0005-2736(02)00495-9PubMed CentralView ArticlePubMedGoogle Scholar
- Khatri P, Draghici S: Ontological analysis of gene expression data: current tools, limitations, and open problems. Bioinformatics 2005.Google Scholar
- Hessner MJ, Xiang B, Jia S, Geoffrey R, Holmes S, Meyer L, Muheisen S, Wang X: Three-color cDNA microarrays with prehybridization quality control yield gene expression data comparable to that of commercial platforms. Physiol Genomics 2006, 25(1):166–178. 10.1152/physiolgenomics.00243.2005View ArticlePubMedGoogle Scholar
- Hessner MJ, Singh VK, Wang X, Khan K, Tschannen MR, Zahrt TC: Visualization and Quality Control of Spotted 70-mer Arrays Using a Labeled Tracking Oligonucleotide. BMC Genomics 2004, 5: 12. 10.1186/1471-2164-5-12PubMed CentralView ArticlePubMedGoogle Scholar
- Sanchez-Margalet V, Lucas M, Solano F, Goberna R: Sensitivity of insulin-secreting RIN m5F cells to undergoing apoptosis by the protein kinase C inhibitor staurosporine. Exp Cell Res 1993, 209(1):160–163. 10.1006/excr.1993.1297View ArticlePubMedGoogle Scholar
- Hessner M, Wang X, Meyer L, Geoffrey R, Jia S, Fuller J, Lernmark A, Ghosh S: Involvement of eotaxin, eosinophils, and pancreatic predisposition in development of type 1 diabetes mellitus in the BioBressding rat. journal of Immunology 2004, 173(11):6993–7002.View ArticleGoogle Scholar
- Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res 2001, 29(9):e45. 10.1093/nar/29.9.e45PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.