Maximizing biomarker discovery by minimizing gene signatures

Prognostic value of a 92-probe signature in breast cancer

Clinical applications of gene expression signatures in breast cancer prognosis still remain limited due to poor predictive strength of single training datasets and appropriate invariable platforms. We proposed a gene expression signature by reducing baseline differences and analyzing common probes among three recent Affymetrix U133 plus 2 microarray data sets. Using a newly developed supervised method, a 92-probe signature found in this study was associated with overall survival. It was robustly validated in four independent data sets and then repeated on three subgroups by incorporating 17 breast cancer microarray datasets. The signature was an independent predictor of patients' survival in univariate analysis [(HR) 1.927, 95% CI (1.237–3.002); p < 0.01] as well as multivariate analysis after adjustment of clinical variables [(HR) 7.125, 95% CI (2.462–20.618); p < 0.001]. Consistent predictive performance was found in different multivariate models in increased patient population (p = 0.002). The survival signature predicted a late metastatic feature through 5-year disease free survival (p = 0.006). We identified subtypes within the lymph node positive (p < 0.001) and ER positive (p = 0.01) patients that best reflected the invasive breast cancer biology. In conclusion using the Common Probe Approach, we present a novel prognostic signature as a predictor in breast cancer late recurrences.

A Bayesian approach for inducing sparsity in generalized linear models with multi-category response

BACKGROUND

The dimension and complexity of high-throughput gene expression data create many challenges for downstream analysis. Several approaches exist to reduce the number of variables with respect to small sample sizes. In this study, we utilized the Generalized Double Pareto (GDP) prior to induce sparsity in a Bayesian Generalized Linear Model (GLM) setting. The approach was evaluated using a publicly available microarray dataset containing 99 samples corresponding to four different prostate cancer subtypes.

RESULTS

A hierarchical Sparse Bayesian GLM using GDP prior (SBGG) was developed to take into account the progressive nature of the response variable. We obtained an average overall classification accuracy between 82.5% and 94%, which was higher than Support Vector Machine, Random Forest or a Sparse Bayesian GLM using double exponential priors. Additionally, SBGG outperforms the other 3 methods in correctly identifying pre-metastatic stages of cancer progression, which can prove extremely valuable for therapeutic and diagnostic purposes. Importantly, using Geneset Cohesion Analysis Tool, we found that the top 100 genes produced by SBGG had an average functional cohesion p-value of 2.0E-4 compared to 0.007 to 0.131 produced by the other methods.

CONCLUSIONS

Using GDP in a Bayesian GLM model applied to cancer progression data results in better subclass prediction. In particular, the method identifies pre-metastatic stages of prostate cancer with substantially better accuracy and produces more functionally relevant gene sets.

Microarray Meta-Analysis and Cross-Platform Normalization: Integrative Genomics for Robust Biomarker Discovery

The diagnostic and prognostic potential of the vast quantity of publicly-available microarray data has driven the development of methods for integrating the data from different microarray platforms. Cross-platform integration, when appropriately implemented, has been shown to improve reproducibility and robustness of gene signature biomarkers. Microarray platform integration can be conceptually divided into approaches that perform early stage integration (cross-platform normalization) versus late stage data integration (meta-analysis). A growing number of statistical methods and associated software for platform integration are available to the user, however an understanding of their comparative performance and potential pitfalls is critical for best implementation. In this review we provide evidence-based, practical guidance to researchers performing cross-platform integration, particularly with an objective to discover biomarkers.

Testing an aflatoxin B1 gene signature in rat archival tissues

Archival tissues from laboratory studies represent a unique opportunity to explore the relationship between genomic changes and agent-induced disease. In this study, we evaluated the applicability of qPCR for detecting genomic changes in formalin-fixed, paraffin-embedded (FFPE) tissues by determining if a subset of 14 genes from a 90-gene signature derived from microarray data and associated with eventual tumor development could be detected in archival liver, kidney, and lung of rats exposed to aflatoxin B1 (AFB1) for 90 days in feed at 1 ppm. These tissues originated from the same rats used in the microarray study. The 14 genes evaluated were Adam8, Cdh13, Ddit4l, Mybl2, Akr7a3, Akr7a2, Fhit, Wwox, Abcb1b, Abcc3, Cxcl1, Gsta5, Grin2c, and the C8orf46 homologue. The qPCR FFPE liver results were compared to the original liver microarray data and to qPCR results using RNA from fresh frozen liver. Archival liver paraffin blocks yielded 30 to 50 μg of degraded RNA that ranged in size from 0.1 to 4 kB. qPCR results from FFPE and fresh frozen liver samples were positively correlated (p ≤ 0.05) by regression analysis and showed good agreement in direction and proportion of change with microarray data for 11 of 14 genes. All 14 transcripts could be amplified from FFPE kidney RNA except the glutamate receptor gene Grin2c; however, only Abcb1b was significantly upregulated from control. Abundant constitutive transcripts, S18 and β-actin, could be amplified from lung FFPE samples, but the narrow RNA size range (25-500 bp length) prevented consistent detection of target transcripts. Overall, a discrete gene signature derived from prior transcript profiling and representing cell cycle progression, DNA damage response, and xenosensor and detoxication pathways was successfully applied to archival liver and kidney by qPCR and indicated that gene expression changes in response to subchronic AFB1 exposure occurred predominantly in the liver, the primary target for AFB1-induced tumors. We conclude that an evaluation of gene signatures in archival tissues can be an important toxicological tool for evaluating critical molecular events associated with chemical exposures.

Genomic signatures and gene networking: challenges and promises

This is an editorial report of the supplement to BMC Genomics that includes 15 papers selected from the BIOCOMP'10 - The 2010 International Conference on Bioinformatics & Computational Biology as well as other sources with a focus on genomics studies.

BIOCOMP'10 was held on July 12-15 in Las Vegas, Nevada. The congress covered a large variety of research areas, and genomics was one of the major focuses because of the fast development in this field. We set out to launch a supplement to BMC Genomics with manuscripts selected from this congress and invited submissions. With a rigorous peer review process, we selected 15 manuscripts that showed work in cutting-edge genomics fields and proposed innovative methodology. We hope this supplement presents the current computational and statistical challenges faced in genomics studies, and shows the enormous promises and opportunities in the genomic future.