Search

Molecular pathway data are essential in current computational and systems biology research. While there are many primary and integrated pathway databases, several challenges remain, including low proteome coverage (57%), low overlap across different databases, unavailability of direct information about underlying physical connectivity of pathway members, and high fraction of protein-coding genes without any pathway annotations, i.e. 'pathway orphans'. In order to address all these challenges, we developed pathDIP, which integrates data from 20 source pathway databases, 'core pathways', with physical protein-protein interactions to predict biologically relevant protein-pathway associations, referred to as 'extended pathways'. Cross-validation determined 71% recovery rate of our predictions. Data integration and predictions increase coverage of pathway annotations for protein-coding genes to 86%, and provide novel annotations for 5732 pathway orphans. PathDIP (http://ophid.utoronto.ca/pathdip) annotates 17 070 protein-coding genes with 4678 pathways, and provides multiple query, analysis and output options., (© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Introduction: Colorectal cancer is the third most common cancer in the world, a small fraction of which is represented by locally advanced rectal cancer (LARC). If not medically contraindicated, preoperative chemoradiotherapy, represent the standard of care for LARC patients. Unfortunately, patients shows a wide range of response rates in which approximately 20% has a complete pathological response, whereas in 20 to 40% the response is poor or absent., Results: The following specific gene signature, able to discriminate responders' patients from non-responders, were founded: AKR1C3, CXCL11, CXCL10, IDO1, CXCL9, MMP12 and HLA-DRA. These genes are mainly involved in immune system pathways and interact with drugs traditionally used in the adjuvant treatment of rectal cancer., Discussion: The present study suggests that new ideas for therapy could be found not only limited to studying genes differentially expressed between the two groups of patients but deepening the mechanisms, associated to response, in which they are involved., Methods: Gene expression studies performed by: Agostini et al., Rimkus et al. and Kim et al. have been merged through a meta-analysis of the raw data. Gene expression data-sets have been processed using A-MADMAN. Common differentially expressed gene (DEG) were identified through SAM analysis. To further characterize the identified DEG we deeply investigated its biological role using an integrative computational biology approach.

Repurposing FDA-approved drugs with the aid of gene signatures of disease can accelerate the development of new therapeutics. A major challenge to developing reliable drug predictions is heterogeneity. Different gene signatures of the same disease or drug treatment often show poor overlap across studies, as a consequence of both biological and technical variability, and this can affect the quality and reproducibility of computational drug predictions. Existing algorithms for signature-based drug repurposing use only individual signatures as input. But for many diseases, there are dozens of signatures in the public domain. Methods that exploit all available transcriptional knowledge on a disease should produce improved drug predictions. Here, we adapt an established meta-analysis framework to address the problem of drug repurposing using an ensemble of disease signatures. Our computational pipeline takes as input a collection of disease signatures, and outputs a list of drugs predicted to consistently reverse pathological gene changes. We apply our method to conduct the largest and most systematic repurposing study on lung cancer transcriptomes, using 21 signatures. We show that scaling up transcriptional knowledge significantly increases the reproducibility of top drug hits, from 44% to 78%. We extensively characterize drug hits in silico, demonstrating that they slow growth significantly in nine lung cancer cell lines from the NCI-60 collection, and identify CALM1 and PLA2G4A as promising drug targets for lung cancer. Our meta-analysis pipeline is general, and applicable to any disease context; it can be applied to improve the results of signature-based drug repurposing by leveraging the large number of disease signatures in the public domain.

Background: Clinical prognostic groupings for localised prostate cancers are imprecise, with 30-50% of patients recurring after image-guided radiotherapy or radical prostatectomy. We aimed to test combined genomic and microenvironmental indices in prostate cancer to improve risk stratification and complement clinical prognostic factors., Methods: We used DNA-based indices alone or in combination with intra-prostatic hypoxia measurements to develop four prognostic indices in 126 low-risk to intermediate-risk patients (Toronto cohort) who will receive image-guided radiotherapy. We validated these indices in two independent cohorts of 154 (Memorial Sloan Kettering Cancer Center cohort [MSKCC] cohort) and 117 (Cambridge cohort) radical prostatectomy specimens from low-risk to high-risk patients. We applied unsupervised and supervised machine learning techniques to the copy-number profiles of 126 pre-image-guided radiotherapy diagnostic biopsies to develop prognostic signatures. Our primary endpoint was the development of a set of prognostic measures capable of stratifying patients for risk of biochemical relapse 5 years after primary treatment., Findings: Biochemical relapse was associated with indices of tumour hypoxia, genomic instability, and genomic subtypes based on multivariate analyses. We identified four genomic subtypes for prostate cancer, which had different 5-year biochemical relapse-free survival. Genomic instability is prognostic for relapse in both image-guided radiotherapy (multivariate analysis hazard ratio [HR] 4·5 [95% CI 2·1-9·8]; p=0·00013; area under the receiver operator curve [AUC] 0·70 [95% CI 0·65-0·76]) and radical prostatectomy (4·0 [1·6-9·7]; p=0·0024; AUC 0·57 [0·52-0·61]) patients with prostate cancer, and its effect is magnified by intratumoral hypoxia (3·8 [1·2-12]; p=0·019; AUC 0·67 [0·61-0·73]). A novel 100-loci DNA signature accurately classified treatment outcome in the MSKCC low-risk to intermediate-risk cohort (multivariate analysis HR 6·1 [95% CI 2·0-19]; p=0·0015; AUC 0·74 [95% CI 0·65-0·83]). In the independent MSKCC and Cambridge cohorts, this signature identified low-risk to high-risk patients who were most likely to fail treatment within 18 months (combined cohorts multivariate analysis HR 2·9 [95% CI 1·4-6·0]; p=0·0039; AUC 0·68 [95% CI 0·63-0·73]), and was better at predicting biochemical relapse than 23 previously published RNA signatures., Interpretation: This is the first study of cancer outcome to integrate DNA-based and microenvironment-based failure indices to predict patient outcome. Patients exhibiting these aggressive features after biopsy should be entered into treatment intensification trials., Funding: Movember Foundation, Prostate Cancer Canada, Ontario Institute for Cancer Research, Canadian Institute for Health Research, NIHR Cambridge Biomedical Research Centre, The University of Cambridge, Cancer Research UK, Cambridge Cancer Charity, Prostate Cancer UK, Hutchison Whampoa Limited, Terry Fox Research Institute, Princess Margaret Cancer Centre Foundation, PMH-Radiation Medicine Program Academic Enrichment Fund, Motorcycle Ride for Dad (Durham), Canadian Cancer Society., (Copyright © 2014 Lalonde et al. Open Access article distributed under the terms of CC BY. Published by Elsevier Ltd. All rights reserved.)

Background: Patients with early-stage non-small-cell lung carcinoma (NSCLC) may benefit from treatments based on more accurate prognosis. A 15-gene prognostic classifier for NSCLC was identified from mRNA expression profiling of tumor samples from the NCIC CTG JBR.10 trial. In this study, we assessed its value in an independent set of cases., Methods: Expression profiling was performed on RNA from frozen, resected tumor tissues corresponding to 181 stage I and II NSCLC cases collected at University Health Network (UHN181). Kaplan-Meier methodology was used to estimate 5-year overall survival probabilities, and the prognostic effect of the classifier was assessed using the log-rank test. A Cox proportional hazards model evaluated the signature's effect adjusting for clinical prognostic factors., Results: Expression data of the 15-gene classifier stratified UHN181 cases into high- and low-risk subgroups with significantly different overall survival (hazard ratio [HR] = 1.92; 95% confidence interval [CI], 1.15-3.23; p = 0.012). In a subgroup analysis, this classifier predicted survival in 127 stage I patients (HR = 2.17; 95% CI, 1.12-4.20; p = 0.018) and the smaller subgroup of 48 stage IA patients (HR = 5.61; 95% CI, 1.19-26.45; p = 0.014). The signature was prognostic for both adenocarcinoma and squamous cell carcinoma cases (HR = 1.76, p = 0.058; HR = 4.19, p = 0.045, respectively)., Conclusion: The prognostic accuracy of a 15-gene classifier was validated in an independent cohort of 181 early-stage NSCLC samples including stage IA cases and in different NSCLC histologic subtypes.

Background: Uterine Leiomyomas (ULs) are the most common benign tumours affecting women of reproductive age. ULs represent a major problem in public health, as they are the main indication for hysterectomy. Approximately 40-50% of ULs have non-random cytogenetic abnormalities, and half of ULs may have copy number alterations (CNAs). Gene expression microarrays studies have demonstrated that cell proliferation genes act in response to growth factors and steroids. However, only a few genes mapping to CNAs regions were found to be associated with ULs., Methodology: We applied an integrative analysis using genomic and transcriptomic data to identify the pathways and molecular markers associated with ULs. Fifty-one fresh frozen specimens were evaluated by array CGH (JISTIC) and gene expression microarrays (SAM). The CONEXIC algorithm was applied to integrate the data., Principal Findings: The integrated analysis identified the top 30 significant genes (P<0.01), which comprised genes associated with cancer, whereas the protein-protein interaction analysis indicated a strong association between FANCA and BRCA1. Functional in silico analysis revealed target molecules for drugs involved in cell proliferation, including FGFR1 and IGFBP5. Transcriptional and protein analyses showed that FGFR1 (P = 0.006 and P<0.01, respectively) and IGFBP5 (P = 0.0002 and P = 0.006, respectively) were up-regulated in the tumours when compared with the adjacent normal myometrium., Conclusions: The integrative genomic and transcriptomic approach indicated that FGFR1 and IGFBP5 amplification, as well as the consequent up-regulation of the protein products, plays an important role in the aetiology of ULs and thus provides data for potential drug therapies development to target genes associated with cellular proliferation in ULs.

Motivation: Penalized regression methods have been adopted widely for high-dimensional feature selection and prediction in many bioinformatic and biostatistical contexts. While their theoretical properties are well-understood, specific methodology for their optimal application to genomic data has not been determined., Results: Through simulation of contrasting scenarios of correlated high-dimensional survival data, we compared the LASSO, Ridge and Elastic Net penalties for prediction and variable selection. We found that a 2D tuning of the Elastic Net penalties was necessary to avoid mimicking the performance of LASSO or Ridge regression. Furthermore, we found that in a simulated scenario favoring the LASSO penalty, a univariate pre-filter made the Elastic Net behave more like Ridge regression, which was detrimental to prediction performance. We demonstrate the real-life application of these methods to predicting the survival of cancer patients from microarray data, and to classification of obese and lean individuals from metagenomic data. Based on these results, we provide an optimized set of guidelines for the application of penalized regression for reproducible class comparison and prediction with genomic data., Availability and Implementation: A parallelized implementation of the methods presented for regression and for simulation of synthetic data is provided as the pensim R package, available at http://cran.r-project.org/web/packages/pensim/index.html., Contact: chuttenh@hsph.harvard.edu; juris@ai.utoronto.ca, Supplementary Information: Supplementary data are available at Bioinformatics online.

Background: Oral Squamous Cell Carcinoma (OSCC) is a major cause of cancer death worldwide, which is mainly due to recurrence leading to treatment failure and patient death. Histological status of surgical margins is a currently available assessment for recurrence risk in OSCC; however histological status does not predict recurrence, even in patients with histologically negative margins. Therefore, molecular analysis of histologically normal resection margins and the corresponding OSCC may aid in identifying a gene signature predictive of recurrence., Methods: We used a meta-analysis of 199 samples (OSCCs and normal oral tissues) from five public microarray datasets, in addition to our microarray analysis of 96 OSCCs and histologically normal margins from 24 patients, to train a gene signature for recurrence. Validation was performed by quantitative real-time PCR using 136 samples from an independent cohort of 30 patients., Results: We identified 138 significantly over-expressed genes (> 2-fold, false discovery rate of 0.01) in OSCC. By penalized likelihood Cox regression, we identified a 4-gene signature with prognostic value for recurrence in our training set. This signature comprised the invasion-related genes MMP1, COL4A1, P4HA2, and THBS2. Over-expression of this 4-gene signature in histologically normal margins was associated with recurrence in our training cohort (p = 0.0003, logrank test) and in our independent validation cohort (p = 0.04, HR = 6.8, logrank test)., Conclusion: Gene expression alterations occur in histologically normal margins in OSCC. Over-expression of the 4-gene signature in histologically normal surgical margins was validated and highly predictive of recurrence in an independent patient cohort. Our findings may be applied to develop a molecular test, which would be clinically useful to help predict which patients are at a higher risk of local recurrence.

Purpose: The JBR.10 trial demonstrated benefit from adjuvant cisplatin/vinorelbine (ACT) in early-stage non-small-cell lung cancer (NSCLC). We hypothesized that expression profiling may identify stage-independent subgroups who might benefit from ACT., Patients and Methods: Gene expression profiling was conducted on mRNA from 133 frozen JBR.10 tumor samples (62 observation [OBS], 71 ACT). The minimum gene set that was selected for the greatest separation of good and poor prognosis patient subgroups in OBS patients was identified. The prognostic value of this gene signature was tested in four independent published microarray data sets and by quantitative reverse-transcriptase polymerase chain reaction (RT-qPCR)., Results: A 15-gene signature separated OBS patients into high-risk and low-risk subgroups with significantly different survival (hazard ratio [HR], 15.02; 95% CI, 5.12 to 44.04; P < .001; stage I HR, 13.31; P < .001; stage II HR, 13.47; P < .001). The prognostic effect was verified in the same 62 OBS patients where gene expression was assessed by qPCR. Furthermore, it was validated consistently in four separate microarray data sets (total 356 stage IB to II patients without adjuvant treatment) and additional JBR.10 OBS patients by qPCR (n = 19). The signature was also predictive of improved survival after ACT in JBR.10 high-risk patients (HR, 0.33; 95% CI, 0.17 to 0.63; P = .0005), but not in low-risk patients (HR, 3.67; 95% CI, 1.22 to 11.06; P = .0133; interaction P < .001). Significant interaction between risk groups and ACT was verified by qPCR., Conclusion: This 15-gene expression signature is an independent prognostic marker in early-stage, completely resected NSCLC, and to our knowledge, is the first signature that has demonstrated the potential to select patients with stage IB to II NSCLC most likely to benefit from adjuvant chemotherapy with cisplatin/vinorelbine.

Background: MicroRNAs (miRs) are non-coding RNA molecules involved in post-transcriptional regulation, with diverse functions in tissue development, differentiation, cell proliferation and apoptosis. miRs may be less prone to degradation during formalin fixation, facilitating miR expression studies in formalin-fixed paraffin-embedded (FFPE) tissue., Results: Our study demonstrates that the TaqMan Human MicroRNA Array v1.0 (Early Access) platform is suitable for miR expression analysis in FFPE tissue with a high reproducibility (correlation coefficients of 0.95 between duplicates, p < 0.00001) and outlines the optimal performance conditions of this platform using clinical FFPE samples. We also outline a method of data analysis looking at differences in miR abundance between FFPE and fresh-frozen samples. By dividing the profiled miR into abundance strata of high (Ct<30), medium (30 < or = Ct < or = 35), and low (Ct>35), we show that reproducibility between technical replicates, equivalent dilutions, and FFPE vs. frozen samples is best in the high abundance stratum. We also demonstrate that the miR expression profiles of FFPE samples are comparable to those of fresh-frozen samples, with a correlation of up to 0.87 (p < 0.001), when examining all miRs, regardless of RNA extraction method used. Examining correlation coefficients between FFPE and fresh-frozen samples in terms of miR abundance reveals correlation coefficients of up to 0.32 (low abundance), 0.70 (medium abundance) and up to 0.97 (high abundance)., Conclusion: Our study thus demonstrates the utility, reproducibility, and optimization steps needed in miR expression studies using FFPE samples on a high-throughput quantitative PCR-based miR platform, opening up a realm of research possibilities for retrospective studies.

In non-small-cell lung cancer (NSCLC), molecular profiling of tumors has led to the identification of gene expression patterns that are associated with specific phenotypes and prognosis. Such correlations could identify early-stage patients who are at increased risk of disease recurrence and death after complete surgical resection and who might benefit from adjuvant therapy. Profiling may also identify aberrant molecular pathways that might lead to specific molecularly targeted therapies. The technology behind the capturing and correlating of molecular profiles with clinical and biologic endpoints have evolved rapidly since microarrays were first developed a decade ago. In this review, we discuss multiple methods that have been used to derive prognostic gene expression signatures in NSCLC. Despite the diversity in the approaches used, 3 main steps are followed. First, the expression levels of several hundred to tens of thousands of genes are quantified by microarray or quantitative polymerase chain reaction techniques; the data are then preprocessed, normalized, and possibly filtered. In the second step, expression data are combined and grouped by clustering, risk score generation, or other means, to generate a gene signature that correlates with a clinical outcome, usually survival. Finally, the signature is validated in datasets of independent cohorts. This review discusses the concepts and methodologies involved in these analytical steps, primarily to facilitate the understanding of reports on large dataset gene expression studies that focus on prognostic signatures in NSCLC.

Resectable non-small-cell lung cancer (NSCLC) patients have poor prognosis, with 30-50% relapsing within 5 years. Current staging criteria do not fully capture the complexity of this disease. Survival could be improved by identification of those early-stage patients who are most likely to benefit from adjuvant therapy. Molecular classification by using mRNA expression profiles has led to multiple, poorly overlapping signatures. We hypothesized that differing statistical methodologies contribute to this lack of overlap. To test this hypothesis, we analyzed our previously published quantitative RT-PCR dataset with a semisupervised method. A 6-gene signature was identified and validated in 4 independent public microarray datasets that represent a range of tumor histologies and stages. This result demonstrated that at least 2 prognostic signatures can be derived from this single dataset. We next estimated the total number of prognostic signatures in this dataset with a 10-million-signature permutation study. Our 6-gene signature was among the top 0.02% of signatures with maximum verifiability, reaffirming its efficacy. Importantly, this analysis identified 1,789 unique signatures, implying that our dataset contains >500,000 verifiable prognostic signatures for NSCLC. This result appears to rationalize the observed lack of overlap among reported NSCLC prognostic signatures.

Although prognostic gene expression signatures for survival in early-stage lung cancer have been proposed, for clinical application, it is critical to establish their performance across different subject populations and in different laboratories. Here we report a large, training-testing, multi-site, blinded validation study to characterize the performance of several prognostic models based on gene expression for 442 lung adenocarcinomas. The hypotheses proposed examined whether microarray measurements of gene expression either alone or combined with basic clinical covariates (stage, age, sex) could be used to predict overall survival in lung cancer subjects. Several models examined produced risk scores that substantially correlated with actual subject outcome. Most methods performed better with clinical data, supporting the combined use of clinical and molecular information when building prognostic models for early-stage lung cancer. This study also provides the largest available set of microarray data with extensive pathological and clinical annotation for lung adenocarcinomas.

Although microarray analysis has provided information regarding the dynamics of gene expression during development of the mouse lung, no extensive correlations have been made to the levels of corresponding protein products. Here, we present a global survey of protein expression during mouse lung organogenesis from embryonic day E13.5 until adulthood using gel-free two-dimensional liquid chromatography coupled to shotgun tandem mass spectrometry (MudPIT). Mathematical modeling of the proteomic profiles with parallel DNA microarray data identified large groups of gene products with statistically significant correlation or divergence in coregulation of protein and transcript levels during lung development. We also present an integrative analysis of mRNA and protein expression in Nmyc loss- and gain-of-function mutants. This revealed a set of 90 positively and negatively regulated putative target genes. These targets are evidence that Nmyc is a regulator of genes involved in mRNA processing and a repressor of the imprinted gene Igf2r in the developing lung.

Gene expression profiling was conducted using primary human nasopharyngeal carcinoma (NPC) biopsy samples to improve the understanding of the molecular pathways defining NPC and to identify novel potential therapeutic targets. RNA samples were extracted from 36 patients suspected to have NPC and hybridized onto the Affymetrix U133A chip. NPC was diagnosed in 19 patients, 11 had lymphoid hyperplasia (LH), and 6 were "normal" biopsies. Clinical stages for these NPC patients ranged from I-IV, including one M1. All NPC patients (except the M1) were treated with curative intent, which included radiotherapy alone (4 patients), or combined with chemotherapy (14 patients). Unsupervised clustering demonstrated a distinct NPC expression pattern, compared to normal biopsies. Subsequent Significance Analysis of Microarrays (SAM) derived from 14 NPC and 6 normal samples discovered 1,089 differentially regulated genes. Pathway analyses revealed novel insights into the mechanisms leading to NPC, whereby upregulation of NFkappaB2 and survivin play central roles in increasing resistance to apoptosis, and changes in integrin and WNT/beta-catenin signaling leading to uncontrolled proliferation. The role of survivin in resisting apoptosis in NPC was confirmed by RNA interference. Our data provide novel insights into the development and progression of NPC, and suggest survivin as a novel therapeutic target for NPC.

Purpose: To explore intratumor heterogeneity in gene expression profiles from patients with cervical cancer., Experimental Design: A total of 33 biopsies were obtained from 11 patients, sampling between two and five different areas for each tumor. The extracted RNA was hybridized onto the Affymetrix U133 Plus 2.0 oligonucleotide chip. The variance of expression within a patient (W), between patients (B) and the total variance (T = W + B) were calculated for each ProbeSet, and the ratio W/T was used as a measure of intratumor heterogeneity. Gene Ontology functional analysis was done to assess the function of genes that had high W/T (top 10%) and low W/T (bottom 10%) values., Results: In total, 448 ProbeSets (2.2% of the total) had W/T < 0.10, indicating low intratumor heterogeneity, and 537 ProbeSets (2.7% of the total) had W/T > 0.90, indicating high intratumor heterogeneity. In total 14,473 ProbeSets (72.4%) had higher intertumor than intratumor heterogeneity (W/T < 0.5). Genes with low intratumor heterogeneity were characterized by a statistically significant enrichment of immune-related functions (P < 0.0001). Genes with high intratumor heterogeneity were characterized by a significant tendency towards nuclear localization and nucleic acid binding (both P < 0.0001). For genes with W/T > 0.5, more than six biopsies would be required to minimize the intratumoral heterogeneity to <0.15; if W/T is 0.3 to 0.4, four biopsies are required; and for low W/T of 0.16 to 0.3, only two to three biopsies would be needed., Conclusion: Although the intratumor heterogeneity was low for the majority of the tested ProbeSets, for many genes, multiple biopsies are required to obtain a reliable estimate of gene expression.

Background: Oxygen plays a central role in human placental pathologies including preeclampsia, a leading cause of fetal and maternal death and morbidity. Insufficient uteroplacental oxygenation in preeclampsia is believed to be responsible for the molecular events leading to the clinical manifestations of this disease., Design: Using high-throughput functional genomics, we determined the global gene expression profiles of placentae from high altitude pregnancies, a natural in vivo model of chronic hypoxia, as well as that of first-trimester explants under 3 and 20% oxygen, an in vitro organ culture model. We next compared the genomic profile from these two models with that obtained from pregnancies complicated by preeclampsia. Microarray data were analyzed using the binary tree-structured vector quantization algorithm, which generates global gene expression maps., Results: Our results highlight a striking global gene expression similarity between 3% O(2)-treated explants, high-altitude placentae, and importantly placentae from preeclamptic pregnancies. We demonstrate herein the utility of explant culture and high-altitude placenta as biologically relevant and powerful models for studying the oxygen-mediated events in preeclampsia., Conclusion: Our results provide molecular evidence that aberrant global placental gene expression changes in preeclampsia may be due to reduced oxygenation and that these events can successfully be mimicked by in vivo and in vitro models of placental hypoxia.

Motivation: High-throughput experiments are being performed at an ever-increasing rate to systematically elucidate protein-protein interaction (PPI) networks for model organisms, while the complexities of higher eukaryotes have prevented these experiments for humans., Results: The Online Predicted Human Interaction Database (OPHID) is a web-based database of predicted interactions between human proteins. It combines the literature-derived human PPI from BIND, HPRD and MINT, with predictions made from Saccharomyces cerevisiae, Caenorhabditis elegans, Drosophila melanogaster and Mus musculus. The 23,889 predicted interactions currently listed in OPHID are evaluated using protein domains, gene co-expression and Gene Ontology terms. OPHID can be queried using single or multiple IDs and results can be visualized using our custom graph visualization program., Availability: Freely available to academic users at http://ophid.utoronto.ca, both in tab-delimited and PSI-MI formats. Commercial users, please contact I.J., Contact: juris@ai.utoronto.ca, Supplementary Information: http://ophid.utoronto.ca/supplInfo.pdf.

We previously reported that our cDNA microarray analysis of primary non-small cell lung carcinoma (NSCLC) could predict for patients at increased risk of cancer recurrence. From the result of this analysis, we selected 11 genes that were considered candidate prognostic marker genes and used the realtime reverse transcription polymerase chain reaction (RT-PCR) to investigate their expression in the same set of NSCLC cases used in the microarray study. Cluster analysis of the realtime RT-PCR data separated these patients into two groups with significantly different disease-free survivals (log-rank test, P < 0.017). In contrast, cluster analysis failed to confirm the prognostic significance of the realtime RT-PCR results for these 11 genes in a validation series of 92 NSCLC cases. In univariate analysis, hypoxia inducible factor 1alpha, Rho-GDP dissociation inhibitor (GDI) alpha (RhoGDI) and Citron/rho-interacting serine-threonine kinase 21 (Citron K21) were significant prognostic factors for disease-free survival in the entire cohort of 130 NSCLC patients, but none were significant in multivariate analysis. The results demonstrate that the prognostic significance of microarray (SAM) results can be partially validated using realtime RT-PCR, but secondary validation using larger and independent series of tumors is necessary to identify true prognostic marker genes.

One of the major concerns in microarray profiling studies of clinical samples is the effect of tissue sampling and RNA extraction on data. We analyzed gene expression in lung cancer specimens that were serially harvested from tumor mass and snap-frozen at several intervals up to 120 minutes after surgical resection. Global gene expression was profiled on cDNA microarrays, and selected stress and hypoxia-activated genes were evaluated using real-time reverse transcription polymerase chain reaction (RT-PCR). Remarkably, similar gene expression profiles were obtained for the majority of samples regardless of the time that had elapsed between resection and freezing. Real-time RT-PCR studies showed significant heterogeneity in the expression levels of stress and hypoxia-activated genes in samples obtained from different areas of a tumor specimen at one time point after resection. The variations between multiple samplings were significantly greater than those of elapsed time between sampling/freezing. Overall samples snap-frozen within 30 to 60 minutes of surgical resection are acceptable for gene expression studies, thus making sampling and snap-freezing of tumor samples in a routine surgical pathology laboratory setting feasible. However, sampling and pooling from multiple sites of each tumor may be necessary for expression profiling studies to overcome the molecular heterogeneity present in tumor specimens.

Motivation: With the increasing number of gene expression databases, the need for more powerful analysis and visualization tools is growing. Many techniques have successfully been applied to unravel latent similarities among genes and/or experiments. Most of the current systems for microarray data analysis use statistical methods, hierarchical clustering, self-organizing maps, support vector machines, or k-means clustering to organize genes or experiments into 'meaningful' groups. Without prior explicit bias almost all of these clustering methods applied to gene expression data not only produce different results, but may also produce clusters with little or no biological relevance. Of these methods, agglomerative hierarchical clustering has been the most widely applied, although many limitations have been identified., Results: Starting with a systematic comparison of the underlying theories behind clustering approaches, we have devised a technique that combines tree-structured vector quantization and partitive k-means clustering (BTSVQ). This hybrid technique has revealed clinically relevant clusters in three large publicly available data sets. In contrast to existing systems, our approach is less sensitive to data preprocessing and data normalization. In addition, the clustering results produced by the technique have strong similarities to those of self-organizing maps (SOMs). We discuss the advantages and the mathematical reasoning behind our approach.

Microarrays of mouse genes are now available from several sources, and they have so far given new insights into gene expression in embryonic development, regions of the brain and during apoptosis. Microarray data posted on the internet can be reanalyzed to study a range of questions.