Quantitative comparison of microarray experiments with published leukemia related gene expression signatures

Background  

Multiple gene expression signatures derived from microarray experiments have been published in the field of leukemia research. A comparison of these signatures with results from new experiments is useful for verification as well as for interpretation of the results obtained. Currently, the percentage of overlapping genes is frequently used to compare published gene signatures against a signature derived from a new experiment. However, it has been shown that the percentage of overlapping genes is of limited use for comparing two experiments due to the variability of gene signatures caused by different array platforms or assay-specific influencing parameters. Here, we present a robust approach for a systematic and quantitative comparison of published gene expression signatures with an exemplary query dataset.     

Genomic signatures to guide the use of chemotherapeutics

Using in vitro drug sensitivity data coupled with Affymetrix microarray data, we developed gene expression signatures that predict sensitivity to individual chemotherapeutic drugs. Each signature was validated with response data from an independent set of cell line studies. We further show that many of these signatures can accurately predict clinical response in individuals treated with these drugs. Notably, signatures developed to predict response to individual agents, when combined, could also predict response to multidrug regimens. Finally, we integrated the chemotherapy response signatures with signatures of oncogenic pathway deregulation to identify new therapeutic strategies that make use of all available drugs. The development of gene expression profiles that can predict response to commonly used cytotoxic agents provides opportunities to better use these drugs, including using them in combination with existing targeted therapies.     

Network based consensus gene signatures for biomarker discovery in breast cancer

Diagnostic and prognostic biomarkers for cancer based on gene expression profiles are viewed as a major step towards a better personalized medicine. Many studies using various computational approaches have been published in this direction during the last decade. However, when comparing different gene signatures for related clinical questions often only a small overlap is observed. This can have various reasons, such as technical differences of platforms, differences in biological samples or their treatment in lab, or statistical reasons because of the high dimensionality of the data combined with small sample size, leading to unstable selection of genes. In conclusion retrieved gene signatures are often hard to interpret from a biological point of view. We here demonstrate that it is possible to construct a consensus signature from a set of seemingly different gene signatures by mapping them on a protein interaction network. Common upstream proteins of close gene products, which we identified via our developed algorithm, show a very clear and significant functional interpretation in terms of overrepresented KEGG pathways, disease associated genes and known drug targets. Moreover, we show that such a consensus signature can serve as prior knowledge for predictive biomarker discovery in breast cancer. Evaluation on different datasets shows that signatures derived from the consensus signature reveal a much higher stability than signatures learned from all probesets on a microarray, while at the same time being at least as predictive. Furthermore, they are clearly interpretable in terms of enriched pathways, disease associated genes and known drug targets. In summary we thus believe that network based consensus signatures are not only a way to relate seemingly different gene signatures to each other in a functional manner, but also to establish prior knowledge for highly stable and interpretable predictive biomarkers.     

From Description to Causality: Mechanisms of Gene Expression Signatures in Cancer

Global gene expression profiles of thousands of cancer samples have been completed, giving rise to hundreds of gene expression signatures (GES). Although many expression signatures show promise in predicting patient prognosis or response to therapies, the usefulness of the signatures in understanding the underlying mechanisms of cancer has not been fully exploited. While “reverse genomic” methods can test specific hypotheses of gene regulation, they fare less well in deciphering novel or combinatorial mechanisms of gene regulation. Recently we described SLAMS (stepwise linkage analysis of microarray signatures), a novel method that can prospectively identify genetic regulators of gene expression signatures in cancer. Applying SLAMS on a poor-prognosis wound signature in human breast cancer, we identified CSN5-mediated ubiquitination of MYC as a novel mechanism to activate a biological program favoring metastasis.     

A simple but highly effective approach to evaluate the prognostic performance of gene expression signatures

Background

Highly parallel analysis of gene expression has recently been used to identify gene sets or ‘signatures’ to improve patient diagnosis and risk stratification. Once a signature is generated, traditional statistical testing is used to evaluate its prognostic performance. However, due to the dimensionality of microarrays, this can lead to false interpretation of these signatures.

Principal Findings

A method was developed to test batches of a user-specified number of randomly chosen signatures in patient microarray datasets. The percentage of random generated signatures yielding prognostic value was assessed using ROC analysis by calculating the area under the curve (AUC) in six public available cancer patient microarray datasets. We found that a signature consisting of randomly selected genes has an average 10% chance of reaching significance when assessed in a single dataset, but can range from 1% to ∼40% depending on the dataset in question. Increasing the number of validation datasets markedly reduces this number.

Conclusions

We have shown that the use of an arbitrary cut-off value for evaluation of signature significance is not suitable for this type of research, but should be defined for each dataset separately. Our method can be used to establish and evaluate signature performance of any derived gene signature in a dataset by comparing its performance to thousands of randomly generated signatures. It will be of most interest for cases where few data are available and testing in multiple datasets is limited.     

Tumor Profiling Turmoil

Tumor profiling studies aim to determine gene expression signatures that can discriminate betweendifferent sub-types of tumors. We have recently discovered a signature that can reliably detectwhich primary head-neck tumors have metastasized to local lymph nodes. This signature has greatpotential for clinical application and also offers unique insights into how metastasis occurs. Despitethese obvious advances, discussed here alongside several other findings, such tumor profilingstudies are currently receiving harsh criticism. We make clear that such evaluations can themselvesalso be critically evaluated. The separation between the two factions actually shows that it is tooearly to either dismiss of exalt tumor profiling studies. Final judgment requires waiting for theresults of larger prospective studies carried out in parallel with current clinical practice.     

Induction of Wnt-Inducible Signaling Protein-1 Correlates with Invasive Breast Cancer Oncogenesis and Reduced Type 1 Cell-Mediated Cytotoxic Immunity: A Retrospective Study

Innate and type 1 cell-mediated cytotoxic immunity function as important extracellular control mechanisms that maintain cellular homeostasis. Interleukin-12 (IL12) is an important cytokine that links innate immunity with type 1 cell-mediated cytotoxic immunity. We recently observed in vitro that tumor-derived Wnt-inducible signaling protein-1 (WISP1) exerts paracrine action to suppress IL12 signaling. The objective of this retrospective study was three fold: 1) to determine whether a gene signature associated with type 1 cell-mediated cytotoxic immunity was correlated with overall survival, 2) to determine whether WISP1 expression is increased in invasive breast cancer, and 3) to determine whether a gene signature consistent with inhibition of IL12 signaling correlates with WISP1 expression. Clinical information and mRNA expression for genes associated with anti-tumor immunity were obtained from the invasive breast cancer arm of the Cancer Genome Atlas study. Patient cohorts were identified using hierarchical clustering. The immune signatures associated with the patient cohorts were interpreted using model-based inference of immune polarization. Reverse phase protein array, tissue microarray, and quantitative flow cytometry in breast cancer cell lines were used to validate observed differences in gene expression. We found that type 1 cell-mediated cytotoxic immunity was correlated with increased survival in patients with invasive breast cancer, especially in patients with invasive triple negative breast cancer. Oncogenic transformation in invasive breast cancer was associated with an increase in WISP1. The gene expression signature in invasive breast cancer was consistent with WISP1 as a paracrine inhibitor of type 1 cell-mediated immunity through inhibiting IL12 signaling and promoting type 2 immunity. Moreover, model-based inference helped identify appropriate immune signatures that can be used as design constraints in genetically engineering better pre-clinical models of breast cancer.     

Development and validation of a prognostic gene-expression signature for lung adenocarcinoma

Although several prognostic signatures have been developed in lung cancer, their application in clinical practice has been limited because they have not been validated in multiple independent data sets. Moreover, the lack of common genes between the signatures makes it difficult to know what biological process may be reflected or measured by the signature. By using classical data exploration approach with gene expression data from patients with lung adenocarcinoma (n = 186), we uncovered two distinct subgroups of lung adenocarcinoma and identified prognostic 193-gene gene expression signature associated with two subgroups. The signature was validated in 4 independent lung adenocarcinoma cohorts, including 556 patients. In multivariate analysis, the signature was an independent predictor of overall survival (hazard ratio, 2.4; 95% confidence interval, 1.2 to 4.8; p = 0.01). An integrated analysis of the signature revealed that E2F1 plays key roles in regulating genes in the signature. Subset analysis demonstrated that the gene signature could identify high-risk patients in early stage (stage I disease), and patients who would have benefit of adjuvant chemotherapy. Thus, our study provided evidence for molecular basis of clinically relevant two distinct two subtypes of lung adenocarcinoma.     

A network-based gene expression signature informs prognosis and treatment for colorectal cancer patients

Background

Several studies have reported gene expression signatures that predict recurrence risk in stage II and III colorectal cancer (CRC) patients with minimal gene membership overlap and undefined biological relevance. The goal of this study was to investigate biological themes underlying these signatures, to infer genes of potential mechanistic importance to the CRC recurrence phenotype and to test whether accurate prognostic models can be developed using mechanistically important genes.

Methods and Findings

We investigated eight published CRC gene expression signatures and found no functional convergence in Gene Ontology enrichment analysis. Using a random walk-based approach, we integrated these signatures and publicly available somatic mutation data on a protein-protein interaction network and inferred 487 genes that were plausible candidate molecular underpinnings for the CRC recurrence phenotype. We named the list of 487 genes a NEM signature because it integrated information from Network, Expression, and Mutation. The signature showed significant enrichment in four biological processes closely related to cancer pathophysiology and provided good coverage of known oncogenes, tumor suppressors, and CRC-related signaling pathways. A NEM signature-based Survival Support Vector Machine prognostic model was trained using a microarray gene expression dataset and tested on an independent dataset. The model-based scores showed a 75.7% concordance with the real survival data and separated patients into two groups with significantly different relapse-free survival (p = 0.002). Similar results were obtained with reversed training and testing datasets (p = 0.007). Furthermore, adjuvant chemotherapy was significantly associated with prolonged survival of the high-risk patients (p = 0.006), but not beneficial to the low-risk patients (p = 0.491).

Conclusions

The NEM signature not only reflects CRC biology but also informs patient prognosis and treatment response. Thus, the network-based data integration method provides a convergence between biological relevance and clinical usefulness in gene signature development.     

Cross-Platform Prediction of Gene Expression Signatures

Gene expression signatures can predict the activation of oncogenic pathways and other phenotypes of interest via quantitative models that combine the expression levels of multiple genes. However, as the number of platforms to measure genome-wide gene expression proliferates, there is an increasing need to develop models that can be ported across diverse platforms. Because of the range of technologies that measure gene expression, the resulting signal values can vary greatly. To understand how this variation can affect the prediction of gene expression signatures, we have investigated the ability of gene expression signatures to predict pathway activation across Affymetrix and Illumina microarrays. We hybridized the same RNA samples to both platforms and compared the resultant gene expression readings, as well as the signature predictions. Using a new approach to map probes across platforms, we found that the genes in the signatures from the two platforms were highly similar, and that the predictions they generated were also strongly correlated. This demonstrates that our method can map probes from Affymetrix and Illumina microarrays, and that this mapping can be used to predict gene expression signatures across platforms.     

Prediction of in vivo radiation dose status in radiotherapy patients using ex vivo and in vivo gene expression signatures

After a large-scale nuclear accident or an attack with an improvised nuclear device, rapid biodosimetry would be needed for triage. As a possible means to address this need, we previously defined a gene expression signature in human peripheral white blood cells irradiated ex vivo that predicts the level of radiation exposure with high accuracy. We now demonstrate this principle in vivo using blood from patients receiving total-body irradiation (TBI). Whole genome microarray analysis has identified genes responding significantly to in vivo radiation exposure in peripheral blood. A 3-nearest neighbor classifier built from the TBI patient data correctly predicted samples as exposed to 0, 1.25 or 3.75 Gy with 94% accuracy (P < 0.001) even when samples from healthy donor controls were included. The same samples were classified with 98% accuracy using a signature previously defined from ex vivo irradiation data. The samples could also be classified as exposed or not exposed with 100% accuracy. The demonstration that ex vivo irradiation is an appropriate model that can provide meaningful prediction of in vivo exposure levels, and that the signatures are robust across diverse disease states and independent sample sets, is an important advance in the application of gene expression for biodosimetry.     

ArrayVigil: a methodology for statistical comparison of gene signatures using segregated-one-tailed (SOT) Wilcoxon's signed-rank test

Due to versatile diagnostic and prognostic fidelity molecular signatures or fingerprints are anticipated as the most powerful tools for cancer management in the near future. Notwithstanding the experimental advancements in microarray technology, methods for analyzing either whole arrays or gene signatures have not been firmly established. Recently, an algorithm, ArraySolver has been reported by Khan for two-group comparison of microarray gene expression data using two-tailed Wilcoxon signed-rank test. Most of the molecular signatures are composed of two sets of genes (hybrid signatures) wherein up-regulation of one set and down-regulation of the other set collectively define the purpose of a gene signature. Since the direction of a selected gene's expression (positive or negative) with respect to a particular disease condition is known, application of one-tailed statistics could be a more relevant choice. A novel method, ArrayVigil, is described for comparing hybrid signatures using segregated-one-tailed (SOT) Wilcoxon signed-rank test and the results compared with integrated-two-tailed (ITT) procedures (SPSS and ArraySolver). ArrayVigil resulted in lower P values than those obtained from ITT statistics while comparing real data from four signatures.     

A signature inferred from Drosophila mitotic genes predicts survival of breast cancer patients

Introduction

The classification of breast cancer patients into risk groups provides a powerful tool for the identification of patients who will benefit from aggressive systemic therapy. The analysis of microarray data has generated several gene expression signatures that improve diagnosis and allow risk assessment. There is also evidence that cell proliferation-related genes have a high predictive power within these signatures.

Methods

We thus constructed a gene expression signature (the DM signature) using the human orthologues of 108 Drosophila melanogaster genes required for either the maintenance of chromosome integrity (36 genes) or mitotic division (72 genes).

Results

The DM signature has minimal overlap with the extant signatures and is highly predictive of survival in 5 large breast cancer datasets. In addition, we show that the DM signature outperforms many widely used breast cancer signatures in predictive power, and performs comparably to other proliferation-based signatures. For most genes of the DM signature, an increased expression is negatively correlated with patient survival. The genes that provide the highest contribution to the predictive power of the DM signature are those involved in cytokinesis.

Conclusion

This finding highlights cytokinesis as an important marker in breast cancer prognosis and as a possible target for antimitotic therapies.     

Improved breast cancer prognosis through the combination of clinical and genetic markers

MOTIVATION: Accurate prognosis of breast cancer can spare a significant number of breast cancer patients from receiving unnecessary adjuvant systemic treatment and its related expensive medical costs. Recent studies have demonstrated the potential value of gene expression signatures in assessing the risk of post-surgical disease recurrence. However, these studies all attempt to develop genetic marker-based prognostic systems to replace the existing clinical criteria, while ignoring the rich information contained in established clinical markers. Given the complexity of breast cancer prognosis, a more practical strategy would be to utilize both clinical and genetic marker information that may be complementary. METHODS: A computational study is performed on publicly available microarray data, which has spawned a 70-gene prognostic signature. The recently proposed I-RELIEF algorithm is used to identify a hybrid signature through the combination of both genetic and clinical markers. A rigorous experimental protocol is used to estimate the prognostic performance of the hybrid signature and other prognostic approaches. Survival data analyses is performed to compare different prognostic approaches. RESULTS: The hybrid signature performs significantly better than other methods, including the 70-gene signature, clinical makers alone and the St. Gallen consensus criterion. At the 90% sensitivity level, the hybrid signature achieves 67% specificity, as compared to 47% for the 70-gene signature and 48% for the clinical makers. The odds ratio of the hybrid signature for developing distant metastases within five years between the patients with a good prognosis signature and the patients with a bad prognosis is 21.0 (95% CI:6.5-68.3), far higher than either genetic or clinical markers alone. AVAILABILITY: The breast cancer dataset is available at www.nature.com and Matlab codes are available upon request.     

Strategy for encoding and comparison of gene expression signatures

EXALT (EXpression signature AnaLysis Tool) is a computational system enabling comparisons of microarray data across experimental platforms and different laboratories . An essential feature of EXALT is a database holding thousands of gene expression signatures extracted from the Gene Expression Omnibus, and encoded in a searchable format. This novel approach to performing global comparisons of shared microarray data may have enormous value when coupled directly with a shared data repository.     

Feasibility of using gene expression analysis to study canine soft tissue sarcomas

The prognosis given for canine soft tissue sarcomas (STSs) is based primarily on histopathologic grade. The decision to administer adjuvant chemotherapy is difficult since less than half of patients with high-grade STSs develop metastatic disease. We hypothesize that there is a gene signature that will improve our ability to predict development of metastatic disease in STS patients. The objective of this study was to determine the feasibility of using cDNA microarray and quantitative real-time PCR (qRT-PCR) analysis to determine gene expression patterns in metastatic versus nonmetastatic canine STSs, given the inherent heterogeneity of this group of tumors. Five STSs from dogs with metastatic disease were evaluated in comparison to eight STSs from dogs without metastasis. Tumor RNA was extracted, processed, and labeled for application to the Affymetrix Canine Genechip 2.0 Array. Array fluorescence was normalized using D-Chip software and data analysis was performed with JMP/Genomics. Differential gene expression was validated using qRT-PCR. Over 200 genes were differentially expressed at a false discovery rate of 5%. Differential gene expression was validated for five genes upregulated in metastatic tumors. Quantitative RT-PCR confirmed increased relative expression of all five genes of interest in the metastatic STSs. Our results demonstrate that microarray and qRT-PCR are feasible methods for comparing gene signatures in canine STSs. Further evaluation of the differences between gene expression in metastatic STSs and in nonmetastatic STSs is likely to identify genes that are important in the development of metastatic disease and improve our ability to prognosticate for individual patients.