Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Research article

How to get the most from microarray data: advice from reverse genomics

Ivan P Gorlov1*, Ji-Yeon Yang2, Jinyoung Byun3, Christopher Logothetis1, Olga Y Gorlova4, Kim-Anh Do5 and Christopher Amos3

Author Affiliations

1 Department of Genitourinary Medical Oncology, Unit 1374, The University of Texas MD Anderson Cancer Center, 1155 Pressler Street, Houston, TX 77030-3721, USA

2 Department of Bioinformatics and Computational Biology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA

3 Department of Community and Family Medicine, Geisel School of Medicine, Dartmouth College, Lebanon, NH, USA

4 Department of Eidemiology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA

5 Department of Biostatistics, The University of Texas MD Anderson Cancer Center, Houston, TX, USA

For all author emails, please log on.

BMC Genomics 2014, 15:223  doi:10.1186/1471-2164-15-223


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/15/223


Received:26 September 2012
Accepted:10 March 2014
Published:21 March 2014

© 2014 Gorlov et al.; licensee 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 credited.

Abstract

Background

Whole-genome profiling of gene expression is a powerful tool for identifying cancer-associated genes. Genes differentially expressed between normal and tumorous tissues are usually considered to be cancer associated. We recently demonstrated that the analysis of interindividual variation in gene expression can be useful for identifying cancer associated genes. The goal of this study was to identify the best microarray data–derived predictor of known cancer associated genes.

Results

We found that the traditional approach of identifying cancer genes—identifying differentially expressed genes—is not very efficient. The analysis of interindividual variation of gene expression in tumor samples identifies cancer-associated genes more effectively. The results were consistent across 4 major types of cancer: breast, colorectal, lung, and prostate. We used recently reported cancer-associated genes (2011–2012) for validation and found that novel cancer-associated genes can be best identified by elevated variance of the gene expression in tumor samples.

Conclusions

The observation that the high interindividual variation of gene expression in tumor tissues is the best predictor of cancer-associated genes is likely a result of tumor heterogeneity on gene level. Computer simulation demonstrates that in the case of heterogeneity, an assessment of variance in tumors provides a better identification of cancer genes than does the comparison of the expression in normal and tumor tissues. Our results thus challenge the current paradigm that comparing the mean expression between normal and tumorous tissues is the best approach to identifying cancer-associated genes; we found that the high interindividual variation in expression is a better approach, and that using variation would improve our chances of identifying cancer-associated genes.

Keywords:
Gene expression; Cancer genes; Interindividual variation in gene expression

Background

Global profiling of gene expression by microarray technology is widely used to study molecular mechanisms of cancer. Even though a number of more sophisticated methods have been developed [1,2] a typical approach to analyze gene expression data is to compare the expression level between normal and primary tumor tissues [3-6]; the genes showing the largest differences in expression are usually considered to be the top candidates as cancer genes.

Recently we hypothesized that high inter-tumor variation in gene expression may more effectively identify cancer-associated genes [7]. Tumors are heterogeneous at the molecular level: in different tumors, different subsets of cancer genes are drivers and therefore are upregulated or downregulated. This leads to a higher inter-tumor variation of the drivers and only slight differences in mean expression values between normal and tumorous tissues. The goal of this study was to comprehensively evaluate that hypothesis.

A previous census of cancer-associated genes identified 400 human cancer genes [8]. The real number of cancer-associated genes is likely to be higher because the cited studies used presence of recurrent somatic mutations in tumor tissue as the only criteria to define cancer-associated genes. We can use known cancer-associated genes to identify which microarray data–derived variables are the best predictors of known cancer genes. We call this a “reverse genomics approach”, and we used it to identify the best predictors of cancer genes for the 4 most common cancers: breast, colorectal, lung, and prostate.

Methods

Datasets and data processing

Figure 1 outlines the design of the study. Table 1 describes the datasets we used. We applied 2 criteria for selecting the datasets: the sample size should be large enough to allow reliable estimates of interindividual variance, and the gene expression data for both tumor and adjacent normal tissues should be available. We used only those probes that could be linked to a single gene. Because different datasets used different gene identifiers, we converted them to Entrez gene identification numbers by using bioDBnet [9]. GeWorkbench 2.3.0 was used to download Simple Omnibus Format in Text files [10]; the data were log2 transformed and normalized by variance-stabilizing algorithm [11] using VCN package in R from bioconductor http://bioconductor.org/ webcite.

thumbnailFigure 1. An outline of the study design.

Table 1. Brief description of the datasets used

Identification of known cancer genes

Our working hypothesis was that inter-individual variance for gene expression values are higher for the genes associated with cancer development because in different tumors different sets of genes are used to drive cancer development. To test this hypothesis we needed to identify genes known to be associated with a given type of cancer. We used text mining tools to identify cancer-associated genes. We evaluated 3 text mining tools: (1) Pathway Studio, (2) Ingenuity Systems, and (3) KnowledgeNet. Results were consistent across the 3 algorithms. Throughout this paper we use Pathway Studio [19] because it uses the most updated databases. The lists of known cancer-associated genes for the 4 cancers we analyzed can be found in Additional file 1. We excluded genes recently reported to be associated with cancer (defined as those reported from January 1, 2011, through July 25, 2012, the date we completed our analysis) because we used recently reported cancer genes for validation.

Additional file 1. Known cancer genes.

Format: XLS Size: 50KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Recently reported cancer genes

Our approach to retrieving recently reported cancer genes is exemplified by breast cancer genes. First we retrieved articles on breast cancer published during the period we defined as recent (from January 1, 2011, through July 25, 2012). In total, 23,550 papers were retrieved by using PubMed. We then scanned the abstracts by using MedScan [20]. One hundred forty-six associations between genes and breast cancer were identified. Seventy genes were novel and not previously reported (Additional file 2); 13 of them were microRNAs.

Additional file 2. Recently identified cancer genes.

Format: XLS Size: 144KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

The same approach was used to retrieve novel colorectal, lung, and prostate cancer genes. The list of the recently reported cancer genes and the corresponding publications is found in Additional file 2.

Microarray data–derived predictors of known cancer-associated genes

We used 6 microarray data–derived predictors: (1) mean gene expression in adjacent normal tissue (m(AN)); (2) mean gene expression in tumor tissue, (m(T)); (3) the degree of change (fold change, or FC) in expression level between tumorous and normal tissue; (4) –LOG(P), in which P is the type-I error identified by using Student’s t test to compare the mean expression between adjacent normal and tumor tissues; (5) the standard deviation of the gene expression in adjacent normal tissue (SD(AN)); and (6) the standard deviation of the gene expression in tumor tissue (SD(T)). Nonparametric Mann-Whitney U (MW) test was used to compare those 6 predictors between cancer-associated genes and all other genes in the human genome.

To estimate how efficiently microarray data–derived predictors predict cancer genes, we ranked the probes by the predictors’ values and estimated the percentage of the cancer-associated genes among the top 5% of the probes. An enrichment factor (EF) was used as a measure of identification efficacy: EF = P/0.05, in which P is the proportion of known cancer genes among the top 5% genes ranked by a given predictor. When EF equals 1, the predictors provide no advantage over random selection of cancer-associated genes; the higher the EF, the better the identification efficacy.

Analysis of outliers

We then identified gene expression outliers in tumor samples separately for cancer and noncancer genes. Outliers for a given gene were defined as tumors with expression level < <a onClick="popup('http://www.biomedcentral.com/1471-2164/15/223/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2164/15/223/mathml/M1">View MathML</a> or > <a onClick="popup('http://www.biomedcentral.com/1471-2164/15/223/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2164/15/223/mathml/M2">View MathML</a>, in which <a onClick="popup('http://www.biomedcentral.com/1471-2164/15/223/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2164/15/223/mathml/M3">View MathML</a> is the mean expression and SD(N) is the standard deviation of the gene expression in adjacent normal tissue. Known cancer genes are more likely to be differently expressed and therefore are more likely to be outliers. To account for the effect of differential gene expression on its probability to be an outlier we first i) sorted all genes according to –LOG(P) from largest to smallest, and for each known cancer gene, and then ii) took non-cancer genes from the list immediately above and below of a given known cancer genes. Those neighbouring non-cancer genes were used as a comparison group.

Computer simulation

We used computer simulation to compare the efficacy of identification of cancer genes by SD(T) and –LOG(P). We simulated the expression levels of 1,000 genes: 50 cancer genes and 950 noncancer genes. The total sample size was 40 tissues: 20 adjacent normal tissues and 20 tumor tissues. Expression values were sampled from the normal distribution with mean = 7 and SD = 0.6 {N(7.0,0.6)}, which are typical means and SDs for the datasets we used.

Two models were compared: the “shifting means model” and the “outlier’s model” (Figure 2). In the shifting means model, tumors are homogeneous: differences in mean expression levels between normal and tumor tissues are due to the shift of the distribution to the right (upregulation) or the left (downregulation). In the outlier’s scenario, a cancer gene is differently expressed in only a fraction of tumors. In different tumors, different cancer genes can be outliers.

thumbnailFigure 2. The “shifting means model” (upper panel) and the “outliers model” (lower panel) of gene expression in tumors. In the shifting means model, all tumors are similar in terms of gene expression. In the outlier’s model, the tumors are heterogeneous: a specific cancer gene is extremely upregulated or downregulated only in a small fraction of tumors in which this gene is a driver of tumorigenesis.

In the shifting means model, expression values were sampled from the normal distribution: N(7.0 + s, 0.6), in which s is the difference in mean expression values between tumor and adjacent normal tissues. Those mean and variance are typical for the studies we have used. Three different values of s (0.07, 0.35, and 0.7) were considered. The number of simulated outliers was defined to make the mean expression value in tumors similar to that of the shifting means model.

Binary logistic regression model

To explore whether the combination of several predictors can provide better identification of known cancer genes than any single predictor can, we ran a binary logistic regression model. A stepwise-forward likelihood ratio model was used to identify significant predictors in each dataset. This analysis was done for each cancer type separately. Known cancer genes were considered as outcome and –LOG(P), m(AN), m(T), FC, SD(AN), and SD(T) as predictors.

Raw and processed data

To estimate the effect of the data processing on the variance estimates, we computed variance in adjacent normal and tumor samples of (1) raw gene expression data, (2) log(2)-transformed data, and (3) data normalized by using a variance-stabilizing approach. We compared variance through these 3 levels using Kendall’s rank-correlation and estimated the efficacy of the variance-based identification of cancer genes for each level.

Results

SD(T) shows the strongest association with known cancer genes

The lists of the known breast, colorectal, lung, and prostate cancer genes are shown in Additional file 2. To check how known cancer genes differ from other genes, we used nonparametric MW testing; Table 2 shows the results of the comparison. SD(T) consistently revealed the strongest differences between cancer and noncancer genes: in 6 of 7 datasets, SD(T) had the largest MW statistics.

Table 2. Differences between cancer and all other genes for 4 cancer gene predictors in 7 datasets

We also ranked the probes according to the predicting variables and estimated percentage of the known cancer genes among the top 5% of the probes. Under the null hypothesis, one can expect 5% of the known cancer genes to be among the top 5% of the ranked probes. We found that all predictors identified more known cancer genes than what one could expect by chance. Figure 3 shows the results of that analysis.

thumbnailFigure 3. The enrichment factor (EF) for known cancer genes among the top 5% of the probes ranked on the basis of the predicting variables. The horizontal lines show the expected proportion of cancer genes under the null hypothesis. Left panel shows individual studies, right panel shows averages across the studies.

Outliers

In all datasets, the percentage of outliers for cancer genes was higher than it was for the comparison group (Table 3). In 3 datasets, the differences were significant. In the overall analysis (nonparametric MW test), the percentage of outliers was overall higher for cancer genes than it was for the paired controls: 3.3 ± 0.2 vs. 2.6 ± 0.1; MW test Z = 3.5; P = 0.0004.

Table 3. Percentage of outliers in tumor samples

Recently reported cancer-associated genes

We next compared the newly reported cancer genes (i.e., reported from January 1, 2011, through July 25, 2012) with all other genes in the human genome (Table 4). SD(T) was the most significant variable for recently identified cancer genes.

Table 4. Differences between recently reported cancer genes and all other genes in the human genome

Computer simulation

Figure 4 illustrates the results of the computer simulation. We found that for the shifting means model, –LOG(P) performed better than SD(T) did; however, for the outliers model, the identification efficacy was better for SD(T).

thumbnailFigure 4. The proportions of correctly predicted cancer genes for the shifting means (left panel) and outliers (right panel) models. The prediction based on – LOG(P) is shown in blue, and that based on SD(T) is shown in red.

Binary logistic regression model

The results of applying the binary logistic regression model to individual datasets are shown in Table 5. We found that SD(T) was the most significant predictor of known cancer genes in all models; in 3 of the 7 models, SD(T) was the single significant predictor. For the other predictors, the results were inconsistent across studies.

Table 5. Results of applying the binary logistic regression model to the 7 datasets

Raw and processed data

We found a strong correlation between variances of log(2)-transformed and normalized data: Kendall’s rank-correlation coefficient varied from 0.94 to 0.98. The prediction efficacy was essentially the same for the 2 types of data. For the raw data, SD(T) was the most significant predictor of known PCa genes. The average EF was slightly lower for the raw data than was it was for the processed data: 3.5 ± 0.3 vs. 3.8 ± 0.2, even though the correlation between the raw-data variance and the processed-data variance was relatively low: average, 0.52; range, 0.31 to 0.62.

Discussion

Our study showed that assessment of the interindividual variation of gene expression is better at predicting cancer-associated genes than is the traditional comparison of mean gene expression in adjacent normal and tumor tissues. The results were consistent for 4 major cancer types: breast, colorectal, lung, and prostate cancer. Additionally, we checked bladder cancer data (dataset GDS1479) [21] and found that SD(T) was also the best predictor of known bladder cancer genes (data not shown).

Why does interindividual variation in gene expression provide a more effective approach to identifying cancer genes than do differences in the mean expression? We believe that tumor heterogeneity at the genetic level is the most likely reason. Although technical errors can contribute to interindividual heterogeneity of gene expression, their effect should not be specific to cancer genes and usually it is much lower than that of gene expression heterogeneity resulting from biologic differences in gene expression [22-24].

It is well recognized that different tumors are driven by different sets of genes [25-28]. Variation among tumors will lead to a substantial increase in interindividual variation of cancer-associated genes and only slight differences in means. The results of our computer simulation support this explanation. When differences in expression levels between normal and tumor tissues are driven by outliers, SD(T) provides a superior prediction of cancer genes than –LOG(P) does. Consistent with these results, we found that cancer genes have a higher probability of being outliers than do other genes with comparable levels of differential expression between tumorous and adjacent normal tissue.

SD-based prediction of cancer genes seems to be robust for data processing. Regardless of whether raw or processed data were used for predicting cancer genes, the average enrichment factor was highest for SD(T). This is likely related to the fact that in extreme cases (i.e., very low or very high variance), it does not really matter whether we are using raw or processed data because the most variable genes tend to keep the same ranking across the different levels of data processing.

The binary logistic regression model identified SD(T) as the most significant and often the only predictor of known cancer genes, whereas the results for the other predictors we tested were inconsistent. Thus, our preliminary analysis does not support the idea that combining several predictors could be better than the cancer genes identification based on SD(T) only.

We noted that interindividual variation in the expression levels of cancer genes was higher not only for tumor samples but also for adjacent normal tissues. This may be a result of selection: adjacent normal tissue is not the same as normal tissue from healthy individuals. People differ by expression levels of cancer genes in normal target tissue, and those differences can contribute to the risk of developing cancer. Somatic alterations, such as methylation or loss or gain of chromosomal regions, may further modulate the expression of cancer genes [29-34], which may explain the better prediction of cancer genes by SD(T) than by SD(AN). The best way to validate this selection hypothesis would be to compare the gene expression in normal tissue (i.e., free of any pathologic changes) with that in adjacent normal tissues from cancer patients. Unfortunately, data on gene expression in normal tissues are usually not available. The only available dataset we found was GSE6919. For that dataset, we found that in normal prostate tissue from healthy individuals, the mean SD for cancer genes was 0.63 ± 0.03, and the mean SD for all other genes was 0.61 ± 0.01. That difference was not statistically significant (MW test, Z = 1.01, P = 0.31). However, in “normal” tissue adjacent to tumor, the genes associated with prostate cancer showed larger interindividual variation in expression than all other genes have: MW test, Z = 4.57, P = 0.000005. This suggests that “normal” tissue from cancer patients is different from normal prostate tissue from healthy individuals. It also suggests that the population of prostate cancer patients is heterogeneous in terms of the expression of prostate cancer genes: different sets of prostate cancer genes are upregulated or downregulated in different patients, leading to a greater interindividual variation in expression.

If a higher interindividual variation in the expression of cancer-associated genes results from genetic heterogeneity, so that different tumors use different sets of cancer genes, one can expect that genetically homogeneous cancers would not show high interindividual variation in expression. Clear cell renal cell carcinoma (CCRCC) is believed to be one of the least heterogeneous cancers [35] with only two major subtypes (ccA and ccB) identified by expression profiling [36]. For CCRCC we used GSE781 GEO dataset [37]. This dataset was selected because it was generated using Affymetrix U133A platform which makes the results comparable with the results on breast, colorectal and lung cancers.

Interindividual variation in expression values of CCRCC-associated genes was not different from interindividual variation in expression values for all other genes, neither in tumor: 0.24 ± 0.01 vs 0.23 ± 0.01; Mann-Whitney U Test Z-adjusted = 0.7, P = 0.48, nor in adjacent normal tissues: 0.23 ± 0.01 vs 0.22 ± 0.01; Mann-Whitney U Test Z-adjusted = 0.3, P = 0.74. Therefore the results of this analysis support the idea that a lower level of interindividual variation in the expression of cancer-associated genes is a result of tumor genetic homogeneity in this particular cancer type. However, in this specific analysis the lack of difference can be due to smaller sample size (9 tumor and 9 adjacent normal samples). This sample size is much smaller compared to the sample sizes used for other types of cancer (Table 1). To address this issue we randomly sampled 9 tumors from each of the datasets used for the analysis of the other types of cancer (Table 1). Twenty random samplings were performed for each datasets, 140 runs in total. For 134 of them or 96%, SD for cancer-associated genes was significantly higher compared to SD for all other genes. This suggests that a smaller sample size is unlikely to explain the lack of differences in SD between cancer-associated and all other genes in CCRCC sample. Therefore, the results of the analysis of CCRCC support the idea that genetic heterogeneity contributes to the higher interindividual variation in the expression of cancer-associated genes, but a larger study is needed to be definitive on this point.

Conclusion

In conclusion, we found that interindividual variation in gene expression more effectively identifies known cancer genes than does the difference in mean expression levels between adjacent normal and tumor tissues. The variation in gene expression levels was more effective at identifying know cancer genes than were differences in mean levels or p-values. Thus, if we use SD(T) instead of the traditional –LOG(P), we would increase our chances of identifying cancer-associated genes. Overall, our results suggest that it would be beneficial to analyze interindividual variation in gene expression.

Abbreviations

EF: Enrichment factor; FC: Fold change in gene expression between tumorous and normal tissue; –LOG(P): P is the type-I error from Student’s t test, comparing mean expression in adjacent normal and tumor tissues; m(AN): Mean gene expression in adjacent normal tissue; m(T): Mean gene expression in tumor tissue; MW: Mann-Whitney U testing; SD(AN): Standard deviation of the gene expression in adjacent normal tissue; SD(T): Standard deviation of the gene expression in tumor tissue.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

IG conceived the study and drafted the manuscript. JY participated in the bioinformatics and statistical analysis, JB participated in the statistical analysis, CL helped draft the manuscript, OG participated in the design of the study and helped to draft the manuscript, KD participated in the design of the study, CA participated in the design of the study and analyses. All authors read and approved the final manuscript.

Acknowledgements

This work was supported in part by the David H. Koch Center for Applied Research of Genitourinary Cancers, the National Institutes of Health Prostate SPORE grant 5 P50 CA140388, the National Institutes of Health U19 CA148127 grant, and the National Institutes of Health Cancer Center Support Grant 5 P30 CA016672. We thank Karen F. Phillips, ELS(D), from the MD Anderson Department of Genitourinary Medical Oncology, for editing the manuscript.

References

  1. Chon HS, Lancaster JM: Microarray-based gene expression studies in ovarian cancer.

    Cancer Control 2011, 18(1):8-15. PubMed Abstract | Publisher Full Text OpenURL

  2. Sorensen KD, Orntoft TF: Discovery of prostate cancer biomarkers by microarray gene expression profiling.

    Expert Rev Mol Diagn 2010, 10(1):49-64. PubMed Abstract | Publisher Full Text OpenURL

  3. Kauffmann A, Huber W: Microarray data quality control improves the detection of differentially expressed genes.

    Genomics 2010, 95(3):138-142. PubMed Abstract | Publisher Full Text OpenURL

  4. Owzar K, Barry WT, Jung SH: Statistical considerations for analysis of microarray experiments.

    Clin Transl Sci 2011, 4(6):466-477. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Suarez E, Burguete A, McLachlan GJ: Microarray data analysis for differential expression: a tutorial.

    P R Health Sci J 2009, 28(2):89-104. PubMed Abstract OpenURL

  6. Zhang Y, Szustakowski J, Schinke M: Bioinformatics analysis of microarray data.

    Methods Mol Biol 2009, 573:259-284. PubMed Abstract | Publisher Full Text OpenURL

  7. Gorlov BJ, Zhao H, Logothetis C, Gorlova O: Beyond comparing means: the usefulness of analyzing interindividual variation in gene expression for identifying genes associated with cancer development.

    J Bioinform Comput Biol 2012, 10(2):16-24. OpenURL

  8. Futreal PA, Coin L, Marshall M, Down T, Hubbard T, Wooster R, Rahman N, Stratton MR: A census of human cancer genes.

    Nat Rev Cancer 2004, 4(3):177-183. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Mudunuri U, Che A, Yi M, Stephens RM: bioDBnet: the biological database network.

    Bioinformatics 2009, 25(4):555-556. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Floratos A, Smith K, Ji Z, Watkinson J, Califano A: geWorkbench: an open source platform for integrative genomics.

    Bioinformatics 2010, 26(14):1779-1780. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Kreil DP, Karp NA, Lilley KS: DNA microarray normalization methods can remove bias from differential protein expression analysis of 2D difference gel electrophoresis results.

    Bioinformatics 2004, 20(13):2026-2034. PubMed Abstract | Publisher Full Text OpenURL

  12. Chen DT, Nasir A, Culhane A, Venkataramu C, Fulp W, Rubio R, Wang T, Agrawal D, McCarthy SM, Gruidl M, Bloom G, Anderson T, White J, Quackenbush J, Yeatman T: Proliferative genes dominate malignancy-risk gene signature in histologically-normal breast tissue.

    Breast Cancer Res Treat 2010, 119(2):335-346. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Graham K, de las Morenas A, Tripathi A, King C, Kavanah M, Mendez J, Stone M, Slama J, Miller M, Antoine G, Willers H, Sebastiani P, Rosenberg CL: Gene expression in histologically normal epithelium from breast cancer patients and from cancer-free prophylactic mastectomy patients shares a similar profile.

    Br J Cancer 2010, 102(8):1284-1293. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Loo LW, Cheng I, Tiirikainen M, Lum-Jones A, Seifried A, Dunklee LM, Church JM, Gryfe R, Weisenberger DJ, Haile RW, Gallinger S, Duggan DJ, Thibodeau SN, Casey G, Le Marchand L: cis-Expression QTL analysis of established colorectal cancer risk variants in colon tumors and adjacent normal tissue.

    PLoS One 2012, 7(2):e30477. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Hou J, Aerts J, den Hamer B, van Ijcken W, den Bakker M, Riegman P, van der Leest C, van der Spek P, Foekens JA, Hoogsteden HC, Grosveld F, Philipsen S: Gene expression-based classification of non-small cell lung carcinomas and survival prediction.

    PLoS One 2010, 5(4):e10312. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Sanchez-Palencia A, Gomez-Morales M, Gomez-Capilla JA, Pedraza V, Boyero L, Rosell R, Farez-Vidal ME: Gene expression profiling reveals novel biomarkers in nonsmall cell lung cancer.

    Int J Cancer 2011, 129(2):355-364. PubMed Abstract | Publisher Full Text OpenURL

  17. Taylor BS, Schultz N, Hieronymus H, Gopalan A, Xiao Y, Carver BS, Arora VK, Kaushik P, Cerami E, Reva B, Antipin Y, Mitsiades N, Landers T, Dolgalev I, Major JE, Wilson M, Socci ND, Lash AE, Heguy A, Eastham JA, Scher HI, Reuter VE, Scardino PT, Sander C, Sawyers CL, Gerald WL: Integrative genomic profiling of human prostate cancer.

    Cancer Cell 2010, 18(1):11-22. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Chandran UR, Ma C, Dhir R, Bisceglia M, Lyons-Weiler M, Liang W, Michalopoulos G, Becich M, Monzon FA: Gene expression profiles of prostate cancer reveal involvement of multiple molecular pathways in the metastatic process.

    BMC Cancer 2007, 7:64. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  19. Yue P, Melamud E, Moult J: SNPs3D: candidate gene and SNP selection for association studies.

    BMC Bioinformatics 2006, 7:166. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  20. Nikitin A, Egorov S, Daraselia N, Mazo I: Pathway studio–the analysis and navigation of molecular networks.

    Bioinformatics 2003, 19(16):2155-2157. PubMed Abstract | Publisher Full Text OpenURL

  21. Dyrskjot L, Kruhoffer M, Thykjaer T, Marcussen N, Jensen JL, Moller K, Orntoft TF: Gene expression in the urinary bladder: a common carcinoma in situ gene expression signature exists disregarding histopathological classification.

    Cancer Res 2004, 64(11):4040-4048. PubMed Abstract | Publisher Full Text OpenURL

  22. Anderson K, Hess KR, Kapoor M, Tirrell S, Courtemanche J, Wang B, Wu Y, Gong Y, Hortobagyi GN, Symmans WF, Pusztai L: Reproducibility of gene expression signature-based predictions in replicate experiments.

    Clin Cancer Res 2006, 12(6):1721-1727. PubMed Abstract | Publisher Full Text OpenURL

  23. Barry WT, Kernagis DN, Dressman HK, Griffis RJ, Hunter JD, Olson JA, Marks JR, Ginsburg GS, Marcom PK, Nevins JR, Geradts J, Datto MB: Intratumor heterogeneity and precision of microarray-based predictors of breast cancer biology and clinical outcome.

    J Clin Oncol 2010, 28(13):2198-2206. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Shi L, Reid LH, Jones WD, Shippy R, Warrington JA, Baker SC, Collins PJ, de Longueville F, Kawasaki ES, Lee KY, Luo Y, Sun YA, Willey JC, Setterquist RA, Fischer GM, Tong W, Dragan YP, Dix DJ, Frueh FW, Goodsaid FM, Herman D, Jensen RV, Johnson CD, Lobenhofer EK, Puri RK, Schrf U, Thierry-Mieg J, Wang C, Wilson M, Wolber PK, et al.: The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements.

    Nat Biotechnol 2006, 24(9):1151-1161. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Grunewald TG, Herbst SM, Heinze J, Burdach S: Understanding tumor heterogeneity as functional compartments–superorganisms revisited.

    J Transl Med 2011, 9:79. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  26. Huland H, Graefen M, Haese A, Hammerer PG, Palisaar J, Pichlmeier U, Henke RP, Erbersdobler A, Huland E, Lilja H: Prediction of tumor heterogeneity in localized prostate cancer.

    Urol Clin North Am 2002, 29(1):213-222. PubMed Abstract | Publisher Full Text OpenURL

  27. Johann DJ, Mukherjee S, Prieto DA, Veenstra TD, Blonder J: Profiling solid tumor heterogeneity by LCM and biological MS of fresh-frozen tissue sections.

    Methods Mol Biol 2011, 755:95-106. PubMed Abstract | Publisher Full Text OpenURL

  28. Marusyk A, Polyak K: Tumor heterogeneity: causes and consequences.

    Biochim Biophys Acta 2010, 1805(1):105-117. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Goering W, Kloth M, Schulz WA: DNA methylation changes in prostate cancer.

    Methods Mol Biol 2012, 863:47-66. PubMed Abstract | Publisher Full Text OpenURL

  30. Lu F, Zhang HT: DNA methylation and nonsmall cell lung cancer.

    Anat Rec (Hoboken) 2011, 294(11):1787-1795. PubMed Abstract | Publisher Full Text OpenURL

  31. Olivier M, Taniere P: Somatic mutations in cancer prognosis and prediction: lessons from TP53 and EGFR genes.

    Curr Opin Oncol 2011, 23(1):88-92. PubMed Abstract | Publisher Full Text OpenURL

  32. Ueno T, Emi M, Sato H, Ito N, Muta M, Kuroi K, Toi M: Genome-wide copy number analysis in primary breast cancer.

    Expert Opin Ther Targets 2012, 16(Suppl 1):S31-S35. PubMed Abstract | Publisher Full Text OpenURL

  33. Watanabe Y, Maekawa M: Methylation of DNA in cancer.

    Adv Clin Chem 2010, 52:145-167. PubMed Abstract OpenURL

  34. Yeang CH, McCormick F, Levine A: Combinatorial patterns of somatic gene mutations in cancer.

    FASEB J 2008, 22(8):2605-2622. PubMed Abstract | Publisher Full Text OpenURL

  35. Crumley SM, Divatia M, Truong L, Shen S, Ayala AG, Ro JY: Renal cell carcinoma: evolving and emerging subtypes.

    World J Clin Cases 2013, 1(9):262-275. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Brannon AR, Reddy A, Seiler M, Arreola A, Moore DT, Pruthi RS, Wallen EM, Nielsen ME, Liu H, Nathanson KL, Ljungberg B, Zhao H, Brooks JD, Ganesan S, Bhanot G, Rathmell WK: Molecular stratification of clear cell renal cell carcinoma by consensus clustering reveals distinct subtypes and survival patterns.

    Genes Cancer 2010, 1(2):152-163. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Lenburg ME, Liou LS, Gerry NP, Frampton GM, Cohen HT, Christman MF: Previously unidentified changes in renal cell carcinoma gene expression identified by parametric analysis of microarray data.

    BMC Cancer 2003, 3:31. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL