Statistical modeling for selecting housekeeper genes

There is a need for statistical methods to identify genes that have minimal variation in expression across a variety of experimental conditions. These 'housekeeper' genes are widely employed as controls for quantification of test genes using gel analysis and real-time RT-PCR. Using real-time quantitative RT-PCR, we analyzed 80 primary breast tumors for variation in expression of six putative housekeeper genes (MRPL19 (mitochondrial ribosomal protein L19), PSMC4 (proteasome (prosome, macropain) 26S subunit, ATPase, 4), SF3A1 (splicing factor 3a, subunit 1, 120 kDa), PUM1 (pumilio homolog 1 (Drosophila)), ACTB (actin, beta) and GAPD (glyceraldehyde-3-phosphate dehydrogenase)). We present appropriate models for selecting the best housekeepers to normalize quantitative data within a given tissue type (for example, breast cancer) and across different types of tissue samples

Correction: Statistical modeling for selecting housekeeper genes

A correction to Statistical modeling for selecting housekeeper genes by Aniko Szabo, Charles M Perou, Mehmet Karaca, Laurent Perreard, John F Quackenbush, and Philip S Bernard. Genome Biology 2004, 5:R5

Classification and risk stratification of invasive breast carcinomas using a real-time quantitative RT-PCR assay

INTRODUCTION: Predicting the clinical course of breast cancer is often difficult because it is a diverse disease comprised of many biological subtypes. Gene expression profiling by microarray analysis has identified breast cancer signatures that are important for prognosis and treatment. In the current article, we use microarray analysis and a real-time quantitative reverse-transcription (qRT)-PCR assay to risk-stratify breast cancers based on biological 'intrinsic' subtypes and proliferation. METHODS: Gene sets were selected from microarray data to assess proliferation and to classify breast cancers into four different molecular subtypes, designated Luminal, Normal-like, HER2+/ER-, and Basal-like. One-hundred and twenty-three breast samples (117 invasive carcinomas, one fibroadenoma and five normal tissues) and three breast cancer cell lines were prospectively analyzed using a microarray (Agilent) and a qRT-PCR assay comprised of 53 genes. Biological subtypes were assigned from the microarray and qRT-PCR data by hierarchical clustering. A proliferation signature was used as a single meta-gene (log(2 )average of 14 genes) to predict outcome within the context of estrogen receptor status and biological 'intrinsic' subtype. RESULTS: We found that the qRT-PCR assay could determine the intrinsic subtype (93% concordance with microarray-based assignments) and that the intrinsic subtypes were predictive of outcome. The proliferation meta-gene provided additional prognostic information for patients with the Luminal subtype (P = 0.0012), and for patients with estrogen receptor-positive tumors (P = 3.4 × 10(-6)). High proliferation in the Luminal subtype conferred a 19-fold relative risk of relapse (confidence interval = 95%) compared with Luminal tumors with low proliferation. CONCLUSION: A real-time qRT-PCR assay can recapitulate microarray classifications of breast cancer and can risk-stratify patients using the intrinsic subtype and proliferation. The proliferation meta-gene offers an objective and quantitative measurement for grade and adds significant prognostic information to the biological subtypes

The molecular portraits of breast tumors are conserved acress microarray platforms

Background Validation of a novel gene expression signature in independent data sets is a critical step in the development of a clinically useful test for cancer patient risk-stratification. However, validation is often unconvincing because the size of the test set is typically small. To overcome this problem we used publicly available breast cancer gene expression data sets and a novel approach to data fusion, in order to validate a new breast tumor intrinsic list. Results A 105-tumor training set containing 26 sample pairs was used to derive a new breast tumor intrinsic gene list. This intrinsic list contained 1300 genes and a proliferation signature that was not present in previous breast intrinsic gene sets. We tested this list as a survival predictor on a data set of 311 tumors compiled from three independent microarray studies that were fused into a single data set using Distance Weighted Discrimination. When the new intrinsic gene set was used to hierarchically cluster this combined test set, tumors were grouped into LumA, LumB, Basal-like, HER2+/ER-, and Normal Breast-like tumor subtypes that we demonstrated in previous datasets. These subtypes were associated with significant differences in Relapse-Free and Overall Survival. Multivariate Cox analysis of the combined test set showed that the intrinsic subtype classifications added significant prognostic information that was independent of standard clinical predictors. From the combined test set, we developed an objective and unchanging classifier based upon five intrinsic subtype mean expression profiles (i.e. centroids), which is designed for single sample predictions (SSP). The SSP approach was applied to two additional independent data sets and consistently predicted survival in both systemically treated and untreated patient groups. Conclusion This study validates the breast tumor intrinsic subtype classification as an objective means of tumor classification that should be translated into a clinical assay for further retrospective and prospective validation. In addition, our method of combining existing data sets can be used to robustly validate the potential clinical value of any new gene expression profile

Basal keratin expression in breast cancer by quantification of mRNA and by immunohistochemistry

Definitions of basal-like breast cancer phenotype vary, and microarray-based expression profiling analysis remains the gold standard for the identification of these tumors. Immunohistochemical identification of basal-like carcinomas is hindered with a fact, that on microarray level not all of them express basal-type cytokeratin 5/6, 14 and 17. We compared expression of cytokeratin 5, 14 and 17 in 115 patients with operable breast cancer estimated by real-time RT-PCR and immunohistochemistry

Expression of centromere protein F (CENP-F) associated with higher FDG uptake on PET/CT, detected by cDNA microarray, predicts high-risk patients with primary breast cancer

<p>Abstract</p> <p>Background</p> <p>Higher standardized uptake value (SUV) detected by 18F-fluorodeoxyglucose positron emission tomography/computed tomography (FDG PET/CT) correlates with proliferation of primary breast cancer. The purpose of this study is to identify specific molecules upregulated in primary breast cancers with a high SUV and to examine their clinical significance.</p> <p>Methods</p> <p>We compared mRNA expression profiles between 14 tumors with low SUVs and 24 tumors with high SUVs by cDNA microarray. We identified centromere protein F (CENP-F) and CDC6 were upregulated in tumors with high SUVs. RT-PCR and immunohistochemical analyses were performed to validate these data. Clinical implication of CENP-F and CDC6 was examined for 253 archival breast cancers by the tissue microarray.</p> <p>Results</p> <p>The relative ratios of CENP-F and CDC6 expression levels to β-actin were confirmed to be significantly higher in high SUV tumors than in low SUV tumors (<it>p </it>= 0.027 and 0.025, respectively) by RT-PCR. In immunohistochemical analysis of 47 node-negative tumors, the CENP-F expression was significantly higher in the high SUV tumors (74%) than the low SUV tumors (45%) (<it>p </it>= 0.04), but membranous and cytoplasmic CDC6 expressions did not significantly differ between both groups (<it>p </it>= 0.9 each). By the tissue microarray, CENP-F (HR = 2.94) as well as tumor size (HR = 4.49), nodal positivity (HR = 4.1), and Ki67 (HR = 2.05) showed independent impact on the patients' prognosis.</p> <p>Conclusion</p> <p>High CENP-F expression, correlated with high SUV, was the prognostic indicators of primary breast cancer. Tumoral SUV levels may serve as a pretherapeutic indicator of aggressiveness of breast cancer.</p

A new molecular breast cancer subclass defined from a large scale real-time quantitative RT-PCR study

BACKGROUND: Current histo-pathological prognostic factors are not very helpful in predicting the clinical outcome of breast cancer due to the disease's heterogeneity. Molecular profiling using a large panel of genes could help to classify breast tumours and to define signatures which are predictive of their clinical behaviour. METHODS: To this aim, quantitative RT-PCR amplification was used to study the RNA expression levels of 47 genes in 199 primary breast tumours and 6 normal breast tissues. Genes were selected on the basis of their potential implication in hormonal sensitivity of breast tumours. Normalized RT-PCR data were analysed in an unsupervised manner by pairwise hierarchical clustering, and the statistical relevance of the defined subclasses was assessed by Chi2 analysis. The robustness of the selected subgroups was evaluated by classifying an external and independent set of tumours using these Chi2-defined molecular signatures. RESULTS: Hierarchical clustering of gene expression data allowed us to define a series of tumour subgroups that were either reminiscent of previously reported classifications, or represented putative new subtypes. The Chi2 analysis of these subgroups allowed us to define specific molecular signatures for some of them whose reliability was further demonstrated by using the validation data set. A new breast cancer subclass, called subgroup 7, that we defined in that way, was particularly interesting as it gathered tumours with specific bioclinical features including a low rate of recurrence during a 5 year follow-up. CONCLUSION: The analysis of the expression of 47 genes in 199 primary breast tumours allowed classifying them into a series of molecular subgroups. The subgroup 7, which has been highlighted by our study, was remarkable as it gathered tumours with specific bioclinical features including a low rate of recurrence. Although this finding should be confirmed by using a larger tumour cohort, it suggests that gene expression profiling using a minimal set of genes may allow the discovery of new subclasses of breast cancer that are characterized by specific molecular signatures and exhibit specific bioclinical features

The molecular portraits of breast tumors are conserved across microarray platforms

BACKGROUND: Validation of a novel gene expression signature in independent data sets is a critical step in the development of a clinically useful test for cancer patient risk-stratification. However, validation is often unconvincing because the size of the test set is typically small. To overcome this problem we used publicly available breast cancer gene expression data sets and a novel approach to data fusion, in order to validate a new breast tumor intrinsic list. RESULTS: A 105-tumor training set containing 26 sample pairs was used to derive a new breast tumor intrinsic gene list. This intrinsic list contained 1300 genes and a proliferation signature that was not present in previous breast intrinsic gene sets. We tested this list as a survival predictor on a data set of 311 tumors compiled from three independent microarray studies that were fused into a single data set using Distance Weighted Discrimination. When the new intrinsic gene set was used to hierarchically cluster this combined test set, tumors were grouped into LumA, LumB, Basal-like, HER2+/ER-, and Normal Breast-like tumor subtypes that we demonstrated in previous datasets. These subtypes were associated with significant differences in Relapse-Free and Overall Survival. Multivariate Cox analysis of the combined test set showed that the intrinsic subtype classifications added significant prognostic information that was independent of standard clinical predictors. From the combined test set, we developed an objective and unchanging classifier based upon five intrinsic subtype mean expression profiles (i.e. centroids), which is designed for single sample predictions (SSP). The SSP approach was applied to two additional independent data sets and consistently predicted survival in both systemically treated and untreated patient groups. CONCLUSION: This study validates the "breast tumor intrinsic" subtype classification as an objective means of tumor classification that should be translated into a clinical assay for further retrospective and prospective validation. In addition, our method of combining existing data sets can be used to robustly validate the potential clinical value of any new gene expression profile

Can Survival Prediction Be Improved By Merging Gene Expression Data Sets?

BACKGROUND:High-throughput gene expression profiling technologies generating a wealth of data, are increasingly used for characterization of tumor biopsies for clinical trials. By applying machine learning algorithms to such clinically documented data sets, one hopes to improve tumor diagnosis, prognosis, as well as prediction of treatment response. However, the limited number of patients enrolled in a single trial study limits the power of machine learning approaches due to over-fitting. One could partially overcome this limitation by merging data from different studies. Nevertheless, such data sets differ from each other with regard to technical biases, patient selection criteria and follow-up treatment. It is therefore not clear at all whether the advantage of increased sample size outweighs the disadvantage of higher heterogeneity of merged data sets. Here, we present a systematic study to answer this question specifically for breast cancer data sets. We use survival prediction based on Cox regression as an assay to measure the added value of merged data sets. RESULTS:Using time-dependent Receiver Operating Characteristic-Area Under the Curve (ROC-AUC) and hazard ratio as performance measures, we see in overall no significant improvement or deterioration of survival prediction with merged data sets as compared to individual data sets. This apparently was due to the fact that a few genes with strong prognostic power were not available on all microarray platforms and thus were not retained in the merged data sets. Surprisingly, we found that the overall best performance was achieved with a single-gene predictor consisting of CYB5D1. CONCLUSIONS:Merging did not deteriorate performance on average despite (a) The diversity of microarray platforms used. (b) The heterogeneity of patients cohorts. (c) The heterogeneity of breast cancer disease. (d) Substantial variation of time to death or relapse. (e) The reduced number of genes in the merged data sets. Predictors derived from the merged data sets were more robust, consistent and reproducible across microarray platforms. Moreover, merging data sets from different studies helps to better understand the biases of individual studies and can lead to the identification of strong survival factors like CYB5D1 expression

The functional loss of the retinoblastoma tumour suppressor is a common event in basal-like and luminal B breast carcinomas

Introduction Breast cancers can be classified using whole genome expression into distinct subtypes that show differences in prognosis. One of these groups, the basal-like subtype, is poorly differentiated, highly metastatic, genomically unstable, and contains specific genetic alterations such as the loss of tumour protein 53 (TP53). The loss of the retinoblastoma tumour suppressor encoded by the RB1 locus is a well-characterised occurrence in many tumour types; however, its role in breast cancer is less clear with many reports demonstrating a loss of heterozygosity that does not correlate with a loss of RB1 protein expression. Methods We used gene expression analysis for tumour subtyping and polymorphic markers located at the RB1 locus to assess the frequency of loss of heterozygosity in 88 primary human breast carcinomas and their normal tissue genomic DNA samples. Results RB1 loss of heterozygosity was observed at an overall frequency of 39%, with a high frequency in basal-like (72%) and luminal B (62%) tumours. These tumours also concurrently showed low expression of RB1 mRNA. p16INK4a was highly expressed in basal-like tumours, presumably due to a previously reported feedback loop caused by RB1 loss. An RB1 loss of heterozygosity signature was developed and shown to be highly prognostic, and was potentially a predictive marker of response to neoadjuvant chemotherapy. Conclusions These results suggest that the functional loss of RB1 is common in basal-like tumours, which may play a key role in dictating their aggressive biology and unique therapeutic responses