Email updates

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

This article is part of the supplement: Proceedings of the 2009 AMIA Summit on Translational Bioinformatics

Open Access Introduction

Transcriptional network classifiers

Hsun-Hsien Chang* and Marco F Ramoni

Author Affiliations

Childrens' Hospital Informatics Program, Harvard-MIT Division of Health Sciences and Technology, Harvard Medical School, Boston, Massachusetts, USA

For all author emails, please log on.

BMC Bioinformatics 2009, 10(Suppl 9):S1-S17  doi:10.1186/1471-2105-10-S9-S1

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2105/10/S9/S1


Published:17 September 2009

© 2009 Chang and Ramoni; 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 cited.

Abstract

Background

Gene interactions play a central role in transcriptional networks. Many studies have performed genome-wide expression analysis to reconstruct regulatory networks to investigate disease processes. Since biological processes are outcomes of regulatory gene interactions, this paper develops a system biology approach to infer function-dependent transcriptional networks modulating phenotypic traits, which serve as a classifier to identify tissue states. Due to gene interactions taken into account in the analysis, we can achieve higher classification accuracy than existing methods.

Results

Our system biology approach is carried out by the Bayesian networks framework. The algorithm consists of two steps: gene filtering by Bayes factor followed by collinearity elimination via network learning. We validate our approach with two clinical data. In the study of lung cancer subtypes discrimination, we obtain a 25-gene classifier from 111 training samples, and the test on 422 independent samples achieves 95% classification accuracy. In the study of thoracic aortic aneurysm (TAA) diagnosis, 61 samples determine a 34-gene classifier, whose diagnosis accuracy on 33 independent samples achieves 82%. The performance comparisons with three other popular methods, PCA/LDA, PAM, and Weighted Voting, confirm that our approach yields superior classification accuracy and a more compact signature.

Conclusions

The system biology approach presented in this paper is able to infer function-dependent transcriptional networks, which in turn can classify biological samples with high accuracy. The validation of our classifier using clinical data demonstrates the promising value of our proposed approach for disease diagnosis.

Background

Genome-wide expression analysis has revolutionized disease diagnostic models through the identification of molecular signatures [1], which are selected from high ranked genes determined by statistical measures, such as fold change [2], t statistic [3], signal-to-noise ratio [4], or subnetwork scores [5]. Over the last decade, system biology researchers also exploited the comprehensive transcriptional landscape offered by microarrays to identify the transcriptional networks that unravel regulatory gene interactions and explain how diseases progress [6-8]. Although these two analysis approaches seem antithetic, they can be unified to create transcriptional network classifiers to enhance disease diagnosis accuracy. We can regard the transcriptional networks underpinning disease development as perturbed by the presence of diseases. The phenotype is treated as a binary perturbation of the overall transcriptional network. To reconstruct the classifier, our task is just to infer from expression profiles the function-dependent transcriptional network that modulates phenotypic traits.

Gene interactions play a central role in transcriptional networks. Abnormal interactions between gene transcripts will give rise to disease incursion [8,9]. To develop transcriptional network classifiers, we consider a system biology approach to capture gene interactions through the measurements of expression collinearity between genes. Our approach is carried out by the Bayesian networks framework, which is a powerful instrument to delineate dependence networks among variables. Bayesian networks have been extensively applied to analyze several types of genomic data, including gene regulation [10], protein-protein interactions [11], single-nucleotide polymorphisms [12] and pedigrees [13]. A Bayesian network is a directed acyclic graph in which nodes represent random variables and arcs define directed dependencies quantified by probability distributions. This study considers a mixed Bayesian network, where the tissue type is represented by a discrete variable and gene expression levels are modelled by continuous, log normal, distributions. Figure 1 illustrates a Bayesian network, where a node represents a gene or a phenotype, and a directed arc linking a pair of nodes records the conditional probability of the child (target) node on the parent (source) node.

thumbnailFigure 1. An example Bayesian network. A node represents a variable, and a directed arc linking a pair of nodes records the conditional probability of the child (target) node on the parent (source) node. In this network, genes 1, 2, 3 are under the phenotype's Markov blanket, so they form a signature for phenotype classification.

Both the graphical structure of a Bayesian network and the parameters of the conditional probabilities can be learned from the available database. Nevertheless, learning a network is computationally intensive because ideally the dependent relations of all pairs of variables must be evaluated. We circumvent the demanding computations by a two-stage learning process. Our algorithm begins with the use of Bayes factor to select the genes that are functionally dependent on the phenotype, since only function-dependent genes have potential to play a role in tissue discrimination. Then, we explore the detailed dependencies between the selected genes to reconstruct a transcriptional network. After the transcriptional network is learned, it can be exploited for tissue classification, again formulated in the Bayesian networks framework. In the learned network, the phenotype's Markov blanket is the set of nodes composed of the phenotype's parents, its children, and its children's parents. Given the genes under the Markov blanket, the phenotype is independent of the genes not covered by the Markov blanket. Hence, only the genes under the Markov blanket contribute to phenotype classification, and they assemble a signature. With reference to Figure 1, genes 1, 2, 3 are those under the phenotype's Markov blanket, consisting of a signature for tissue classification.

Results

We validate our approach by two clinical studies: discrimination of lung cancer subtypes and diagnosis of thoracic aortic aneurysm.

Discrimination of lung cancer subtypes

Lung adenocarcinoma (AC) and squamous cell carcinoma (SCC) are the most common subtypes of lung cancer. They are heterogeneous in many clinical aspects, such as responses to chemotherapy [14], tendency to metastases [15,16], and mortality rates [17,18]. Unfortunately, the current gold standard is histology which is subjective [19] and may fail when tumors are small [20] or when patients suffer from multiple types of primary lung carcinomas [21]. Gene expression profiling will avoid these problems and perform automatic discrimination of lung cancer subtypes. The classifier is trained by a Duke University data [22], which is available on Gene Expression Omnibus with accession number GSE3141, in a total of 58 ACs and 53 SCCs. The lung specimens are assayed by Affymetrix HG-U133A. Figure 2 shows the function-dependent transcriptional network inferred from the data. Of the 22,283 gene probes in the microarray, seventy seven probes are dependent, directly or indirectly, on the carcinoma subtypes. Of these 77 genes, 25 are under the phenotype Markov blanket, so they per se assemble a signature. Enrichment study shows that there are 23 unique genes in this signature, summarized in Table 1.

thumbnailFigure 2. The functional dependence network for lung cancer subtypes characterization. There are 77 genes dependent on the lung cancer subtypes and they are selected to build up this network, where 25 genes (in green) are under the phenotype's Markov blanket to assemble a signature.

Table 1. The signature of 25 genes for characterizing lung cancer subtypes. Enrichment shows that there are 23 unique genes in the signature.

The performance of 10-fold cross validation achieves 98.5% accuracy. We further test the classification accuracy of the network on seven independent study populations with Gene Expression Omnibus accession numbers GSE10072, GSE7670, GSE12667, GSE4824, GSE2109, GSE4573, and GSE6253, for a total of 422 samples, 232 AC and 190 SCC, from subjects of Caucasian, Asian and African descent representing 84.6%, 6.9%, and 2.8% of the data, respectively. On these independent samples, our transcriptional network classifier achieves an accuracy of 95.2%.

The 25-gene signature identified by the classifier is unique to discriminate AC and SCC with high accuracy. Furthermore, most of these genes have been reported their specificity to lung cancer. ABCC3, CLDN3, DPP4, MUC3B, MUC5B, NTRK2, SPINK1, TJP3 are specific markers of lung AC [23-29]. KRT6A, KRT6B, KRT6C, KRT17, RHCG, SPRR1A, and VSNL1 are unique to lung SCC [30-33]. BICD2, CDA, NMNAT2, SERPINB13, and TOX3 have no specificity to either AC or SCC but to lung cancer [34-38].

Diagnosis of thoracic aortic aneurysm

Thoracic aortic aneurysm (TAA) is usually asymptomatic and associated with high mortality. Identification of at-risk individuals is a challenging task. Gene expression patterns in peripheral blood cells are expected to assist the diagnosis of TAA. The data used to derive the classifier is publicly available on Gene Expression Omnibus with accession number GSE9106 [39], which involves 36 cases and 25 controls for training purpose. Peripheral blood samples were collected at Yale-New Haven Hospital. Gene expression experiments were carried out by Applied Biosystems Human Genome Survey Microarray v2.0, which is equipped with 32,878 probes. The utilization of Bayes factor in our algorithm first filters out 346 genes that are dependent on the phenotype. Bayesian network learning results in the functional dependence network shown in Figure 3. There are 34 genes under the phenotype's Markov blanket, and they form a signature for TAA diagnosis. Table 2 summarizes the annotations of the signature, where the nameless genes are provided with their probe identifies only. The genes ABCG4, ARNT2, BCOR, CABP2, CSTF2, DNTTIP1, FGG, IGF2BP1, MAL2, MMP11, RBM16, TM4SF1, ZBTB4, ZNF394 are involved in connective tissue disorders and inflammatory disease, which are prerequisite to TAA. The 10-fold cross validation of the classifier yields 97% accuracy. We further examine the classifier on the independent samples, 24 cases and 9 controls, also included in the Yale data GSE9106. The accuracy on these independent samples achieves 82%, demonstrating good performance of our approach.

thumbnailFigure 3. The functional dependence network for TAA diagnosis. There are 346 genes selected to reconstruct this network, because of their distinct expression patterns between TAA and normal samples. The signature consists of the 34 genes (in green) under the phenotype’s Markov blanket.

Table 2. The signature of 34 genes for diagnosing TAA.

Comparisons with other methods

We contrast our proposed system biology approach with other popular algorithms that do not take into account regulatory gene interactions:

1) Principal Component Analysis with Linear Discriminant Analysis (PCA/LDA): The PCA/LDA method begins with reducing the number of genes to a small number of principal genes and then searches for a discriminative linear function on expression values to separate tissues.

2) Prediction Analysis for Microarray (PAM) [40]: PAM utilizes signal to noise ratios to pick up a signature and uses the ratios to determine the tissue types of testing samples.

3) Weighted Voting [1]: This method ranks genes by the fold change of the means of the expression values. The classification is determined by how close to the high rank genes the testing data is.

Unlike our approach, the above methods neglect dependencies among genes, so they yield worse performance than our TNC. Table 3 and Table 4 summarize the comparisons of our approach with these methods on the lung cancer and TAA studies, respectively. The results show that our approach is superior to other algorithms. On the other hand, our approach leads to more compact signatures because collinearity elimination is addressed after gene selection. The differences between our approach and other schemes are statistically significant (p<0.005), except that weighted voting performs close to ours in the lung cancer study. Although weighted voting reaches high classification accuracy on the lung cancer data, it requires a large number of genes in the signature, giving rise to overfitting problem.

Table 3. Performance comparisons with other methods on the lung cancer data.

Table 4. Performance comparisons with other methods on the TAA data.

Discussion

The clinical application confirms improved accuracy of our proposed system biology approach. Literature survey on the functions of the signature genes also validates the capability of our approach to extract biologically reasonable signatures. Furthermore, the large-scale independent test on seven cohorts in the lung cancer study shows robustness of our classifier across platforms and populations. The two studies also demonstrate the capability of our method to analyze data assayed by microarrays manufactured by different makers.

Unlike existing methods that require the operator to specify a cutoff of statistical measures to select high ranked genes, our method is threshold free for signature selection, because the signature genes are determined once the transcriptional network is modelled. For phenotype classification, we need to keep the network merely composing of the signature genes, and the remaining network can be discarded; this way can save storage resources in clinical usage. Another feature of our transcriptional network classifier is its visualization of molecular dependence network, which will provide biologists a clue for gene causality investigation.

A recent work proposes to use prior knowledge of known pathway information to select gene subnetworks as features for tissue classification [5]. However, this method will discard a major portion of the data, because a large number of genes have not been discovered their functional pathways. Dissimilar to this method, our approach fully utilizes the entire data to screen the function-dependent genes and to reconstruct the network.

Conclusions

This paper uses a system biology approach to develop transcriptional network classifiers. The classifier can be thought of as a gene network perturbed by the presence of the phenotypic traits. We adopt Bayesian network framework to model the classifier. The algorithm uses Bayes factor for gene filtering, followed by collinearity elimination via network learning. The clinical applications of our approach to lung cancer subtypes classification and TAA diagnosis demonstrate high classification accuracy of the network based classifiers. The biological validation of the signatures further confirms the ability of the transcriptional network classifier to extract meaningful signatures.

Methods

Let Y1,Y2,...,YN be Gaussian random variables representing the expression levels of genes 1,...,N, and C be a multinomial random variable indicating tissue conditions. We use uppercase to denote random variables and lowercase to denote their values. Our algorithm first uses Bayes factor to filter function-dependent genes and then exploit Bayesian network learning to eliminate collinearity among these selected genes.

Gene filtering by Bayes factor

The genes functionally dependent on the phenotype are filtered in the beginning. The filtering can be realized by Bayes factor, which evaluates for each gene the ratio of its likelihood of being dependent on the phenotype to its likelihood of being independent of the phenotype. When the Bayes factor is greater than one, the gene is selected because it is more likely to be dependent on than to be independent of the phenotype.

Collinearity elimination via network learning

Without loss of generality, we assume that the first G out of N genes were selected by the preceding step. The gene expression data under consideration now is D = {y1,y2,...,YG,c}. When a gene Yi is collinearly expressed with another gene Yj, the dependence of gene Yi on the phenotype is mediated by gene Yj. In other words, our goal is to search which gene modulates gene Yi with the highest likelihood. When we find out for every gene its best upstream variable, the network is achieved. In the framework of Bayesian network, our objective is to learn from a set of candidate network models the optimal network fit best to the data D. Equivalently, we look for the highest posterior probability . Applying Bayes' theorem to results in , where is the prior probability of model and is the marginal likelihood. The computation of is to average out from the likelihood function , where is the values of the random vector parameterizing the distribution of y1,y2,...,YG,C conditional on . We can exploit the local Markov properties encoded by the network to rewrite the joint probability as

,

where pa(x) denotes the values of the parents Pa(X) of random variables X, and is the subset of parameters used to describe the dependence of variable X on its parents.

In this paper, we model a gene Yg to be dependent on either the phenotype C or another single geneYa, and the phenotype C is the root in the network without parents. We further can assume the J samples in the database are independent. The likelihood function becomes

,

where the subscripts j indicate the jth sample. The first term can be estimated by sample frequencies, and the second term can be derived using linear Gaussian model [41]. The marginal likelihood function is the solution of the integral

,

Due to limited space, we in this paper do not present the detailed computation, which can be derived from [41]. Finally, the determination of the best Bayesian network model is .

Sample classification

The phenotype classification of a sample is to find the maximum probability of the tissue class that the sample belongs to, conditional on the expression values of the sample. The formulation for the classification is as follows:

,

The application of Bayes' theorem leads to

,

where the second equality holds because the denominator p(y1,y2,...,yG) in the first line is not a function of c. Since only genes directly dependent on the phenotype variable C matter in the maximization, the tissue classification becomes

,

where H denotes the set of genes that are the children of the phenotype C in the network and assemble a signature. Equivalently, the set H of genes corresponds to the genes under the phenotype's Markov blanket.

Competing interests

The authors declare that they have no competing interest.

Authors' contributions

HHC designed the method and conducted the analysis; MFR directed the study; both authors prepared the manuscript.

Acknowledgements

This research is supported in part by NIH/NHGRI (R01HG003354).

This article has been published as part of BMC Bioinformatics Volume 10 Supplement 9, 2009: Proceedings of the 2009 AMIA Summit on Translational Bioinformatics. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/10?issue=S9.

References

  1. Golub TR, Slonim DK, Tamayo P, Huard C, Gaasenbeek M, Mesirov JP, Coller H, Loh ML, Downing JR, Caligiuri MA, et al.: Molecular classification of cancer: class discovery and class prediction by gene expression monitoring.

    Science 1999, 286(5439):531-537. PubMed Abstract | Publisher Full Text OpenURL

  2. Chen Y, Dougherty ER, Bittner ML: Ratio-based decisions and the quantitative analysis of cDNA microarray images.

    J Biomedical Optics 1997, 2(4):364-374. Publisher Full Text OpenURL

  3. Reich M, Ohm K, Angelo M, Tamayo P, Mesirov JP: GeneCluster 2.0: an advanced toolset for bioarray analysis.

    Bioinformatics 2004, 20(11):1797-1798. PubMed Abstract | Publisher Full Text OpenURL

  4. Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response.

    Proc Natl Acad Sci U S A 2001, 98(9):5116-5121. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Lee E, Chuang HY, Kim JW, Ideker T, Lee D: Inferring pathway activity toward precise disease classification.

    PLoS Comput Biol 2008, 4(11): e1000217. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Huang E, Ishida S, Pittman J, Dressman H, Bild A, Kloos M, D'Amico M, Pestell RG, West M, Nevins JR: Gene expression phenotypic models that predict the activity of oncogenic pathways.

    Nature genetics 2003, 34(2):226-230. PubMed Abstract | Publisher Full Text OpenURL

  7. Lamb J, Ramaswamy S, Ford HL, Contreras B, Martinez RV, Kittrell FS, Zahnow CA, Patterson N, Golub TR, Ewen ME: A mechanism of cyclin D1 action encoded in the patterns of gene expression in human cancer.

    Cell 2003, 114(3):323-334. PubMed Abstract | Publisher Full Text OpenURL

  8. Rhodes DR, Yu J, Shanker K, Deshpande N, Varambally R, Ghosh D, Barrette T, Pandey A, Chinnaiyan AM: Large-scale meta-analysis of cancer microarray data identifies common transcriptional profiles of neoplastic transformation and progression.

    Proceedings of the National Academy of Sciences of the United States of America 2004, 101(25):9309-9314. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Abdollahi A, Schwager C, Kleeff J, Esposito I, Domhan S, Peschke P, Hauser K, Hahnfeldt P, Hlatky L, Debus J, et al.: Transcriptional network governing the angiogenic switch in human pancreatic cancer.

    Proceedings of the National Academy of Sciences of the United States of America 2007, 104(31):12890-12895. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Friedman N: Inferring cellular networks using probabilistic graphical models.

    Science 2004, 303:799-805. PubMed Abstract | Publisher Full Text OpenURL

  11. Jansen R, Yu H, Greenbaum D, Kluger Y, Krogan NJ, Chung S, Emili A, Snyder M, Greenblatt JF, Gerstein M: A Bayesian networks approach for predicting protein-protein interactions from genomic data.

    Science 2003, 302:449-453. PubMed Abstract | Publisher Full Text OpenURL

  12. Sebastiani P, Ramoni MF, Nolan V, Baldwin CT, Steinberg MH: Genetic dissection and prognostic modeling of overt stroke in sickle cell anemia.

    Nature genetics 2005, 37(4):435-440. PubMed Abstract | Publisher Full Text OpenURL

  13. Lauritzen SL, Sheehan NA: Graphical models for genetic analysis.

    Statist Sci 2004, 18(4):489-514. OpenURL

  14. Kato H, Ichinose Y, Ohta M, Hata E, Tsubota N, Tada H, Watanabe Y, Wada H, Tsuboi M, Hamajima N: A randomized trial of adjuvant chemotherapy with uracil-tegafur for adenocarcinoma of the lung.

    N Engl J Med 2004, 350(17):1713-1721. PubMed Abstract | Publisher Full Text OpenURL

  15. Thomas P, Khokha R, Shepherd FA, Feld R, Tsao MS: Differential expression of matrix metalloproteinases and their inhibitors in non-small cell lung cancer.

    J Pathol 2000, 190(2):150-156. PubMed Abstract | Publisher Full Text OpenURL

  16. Yu CJ, Shih JY, Lee YC, Shun CT, Yuan A, Yang PC: Sialyl Lewis antigens: association with MUC5AC protein and correlation with post-operative recurrence of non-small cell lung cancer.

    Lung Cancer 2005, 47(1):59-67. PubMed Abstract | Publisher Full Text OpenURL

  17. Nesbitt JC, Putnam JB Jr., Walsh GL, Roth JA, Mountain CF: Survival in early-stage non-small cell lung cancer.

    Ann Thorac Surg 1995, 60(2):466-472. PubMed Abstract | Publisher Full Text OpenURL

  18. Okamoto T, Maruyama R, Suemitsu R, Aoki Y, Wataya H, Kojo M, Ichinose Y: Prognostic value of the histological subtype in completely resected non-small cell lung cancer.

    Interact Cardiovasc Thorac Surg 2006, 5(4):362-366. PubMed Abstract | Publisher Full Text OpenURL

  19. Jamieson LA, Carey FA: Pathology of lung tumours.

    SURGERY 2005, 23(11):389-393. OpenURL

  20. Wistuba II, Gazdar AF: Lung cancer preneoplasia.

    Annu Rev Pathol 2006, 1:331-348. PubMed Abstract | Publisher Full Text OpenURL

  21. Nonami Y, Ohtuki Y, Sasaguri S: Study of the diagnostic difference between the clinical diagnostic criteria and results of immunohistochemical staining of multiple primary lung cancers.

    J Cardiovasc Surg (Torino) 2003, 44(5):661-665. PubMed Abstract OpenURL

  22. Bild A, Yao G, Chang J, Wang Q, Potti A, Chasse D, Joshi M, Harpole D, Lancaster J, Berchuck A, et al.: Oncogenic pathway signatures in human cancers as a guide to targeted therapies.

    Nature 2006, 439(7074):353-357. PubMed Abstract | Publisher Full Text OpenURL

  23. Hanada S, Maeshima A, Matsuno Y, Ohta T, Ohki M, Yoshida T, Hayashi Y, Yoshizawa Y, Hirohashi S, Sakamoto M: Expression profile of early lung adenocarcinoma: identification of MRP3 as a molecular marker for early progression.

    J Pathol 2008, 216(1):75-82. PubMed Abstract | Publisher Full Text OpenURL

  24. Kuner R, Muley T, Meister M, Ruschhaupt M, Buness A, Xu EC, Schnabel P, Warth A, Poustka A, Sultmann H, et al.: Global gene expression analysis reveals specific patterns of cell junctions in non-small cell lung cancer subtypes.

    Lung Cancer 2009, 63(1):32-38. PubMed Abstract | Publisher Full Text OpenURL

  25. Wesley UV, Tiwari S, Houghton AN: Role for dipeptidyl peptidase IV in tumor suppression of human non small cell lung carcinoma cells.

    Int J Cancer 2004, 109(6):855-866. PubMed Abstract | Publisher Full Text OpenURL

  26. Nguyen PL, Niehans GA, Cherwitz DL, Kim YS, Ho SB: Membrane-bound (MUC1) and secretory (MUC2, MUC3, and MUC4) mucin gene expression in human lung cancer.

    Tumour Biol 1996, 17(3):176-192. PubMed Abstract | Publisher Full Text OpenURL

  27. Copin M, Buisine M, Leteurtre E, Marquette C, Porte H, Aubert J, Gosselin B, Porchet N: Mucinous bronchioloalveolar carcinomas display a specific pattern of mucin gene expression among primary lung adenocarcinomas.

    Hum Pathol 2001, 32(3):274-281. PubMed Abstract | Publisher Full Text OpenURL

  28. Ding L, Getz G, Wheeler DA, Mardis ER, McLellan MD, Cibulskis K, Sougnez C, Greulich H, Muzny DM, Morgan MB, et al.: Somatic mutations affect key pathways in lung adenocarcinoma.

    Nature 2008, 455(7216):1069-1075. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Borczuk A, Kim H, Yegen H, Friedman R, Powell C: Lung adenocarcinoma global profiling identifies type II transforming growth factor-beta receptor as a repressor of invasiveness.

    Am J Respir Crit Care Med 2005, 172(6):729-737. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Hawthorn L, Stein L, Panzarella J, Loewen G, Baumann H: Characterization of cell-type specific profiles in tissues and isolated cells from squamous cell carcinomas of the lung.

    Lung Cancer 2006, 53(2):129-142. PubMed Abstract | Publisher Full Text OpenURL

  31. Fujii T, Dracheva T, Player A, Chacko S, Clifford R, Strausberg R, Buetow K, Azumi N, Travis W, Jen J: A preliminary transcriptome map of non-small cell lung cancer.

    Cancer Res 2002, 62(12): 3340-3346. PubMed Abstract | Publisher Full Text OpenURL

  32. Chen BS, Xu ZX, Xu X, Cai Y, Han YL, Wang J, Xia SH, Hu H, Wei F, Wu M, et al.: RhCG is downregulated in oesophageal squamous cell carcinomas, but expressed in multiple squamous epithelia.

    Eur J Cancer 2002, 38(14):1927-1936. PubMed Abstract | Publisher Full Text OpenURL

  33. Fu J, Fong K, Bellacosa A, Ross E, Apostolou S, Bassi DE, Jin F, Zhang J, Cairns P, Ibanez de Caceres I, et al.: VILIP-1 downregulation in non-small cell lung carcinomas: mechanisms and prediction of survival.

    PLoS ONE 2008, 3(2):e1698. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Guha U, Chaerkady R, Marimuthu A, Patterson AS, Kashyap MK, Harsha HC, Sato M, Bader JS, Lash AE, Minna JD, et al.: Comparisons of tyrosine phosphorylated proteins in cells expressing lung cancer-specific alleles of EGFR and KRAS.

    Proc Natl Acad Sci U S A 2008, 105(37):14112-14117. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Tibaldi C, Giovannetti E, Vasile E, Mey V, Laan AC, Nannizzi S, Di Marsico R, Antonuzzo A, Orlandini C, Ricciardi S, et al.: Correlation of CDA, ERCC1, and XPD polymorphisms with response and survival in gemcitabine/cisplatin-treated advanced non-small cell lung cancer patients.

    Clin Cancer Res 2008, 14(6):1797-1803. PubMed Abstract | Publisher Full Text OpenURL

  36. Chari R, Lonergan KM, Ng RT, MacAulay C, Lam WL, Lam S: Effect of active smoking on the human bronchial epithelium transcriptome.

    BMC Genomics 2007, 8:297. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  37. Heighway J, Knapp T, Boyce L, Brennand S, Field JK, Betticher DC, Ratschiller D, Gugger M, Donovan M, Lasek A, et al.: Expression profiling of primary non-small cell lung cancer for target identification.

    Oncogene 2002, 21(50):7749-7763. PubMed Abstract | Publisher Full Text OpenURL

  38. Hu Z, Chen J, Tian T, Zhou X, Gu H, Xu L, Zeng Y, Miao R, Jin G, Ma H, et al.: Genetic variants of miRNA sequences and non-small cell lung cancer survival.

    J Clin Invest 2008, 118(7):2600-2608. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Wang Y, Barbacioru CC, Shiffman D, Balasubramanian S, Iakoubova O, Tranquilli M, Albornoz G, Blake J, Mehmet NN, Ngadimo D, et al.: Gene expression signature in peripheral blood detects thoracic aortic aneurysm.

    PLoS ONE 2007, 2(10):e1050. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Tibshirani R, Hastie T, Narasimhan B, Chu G: Diagnosis of multiple cancer types by shrunken centroids of gene expression.

    Proc Natl Acad Sci U S A 2002, 99(10):6567-6572. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  41. Ferrazzi F, Sebastiani P, Ramoni MF, Bellazzi R: Bayesian approaches to reverse engineer cellular systems: a simulation study on nonlinear Gaussian networks.

    BMC Bioinformatics 2007, 8(Suppl 5):S2. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL