Sample Preparation Approach Influences PAM50 Risk of Recurrence Score in Early Breast Cancer.

The PAM50 gene expression subtypes and the associated risk of recurrence (ROR) score are used to predict the risk of recurrence and the benefits of adjuvant therapy in early-stage breast cancer. The Prosigna assay includes the PAM50 subtypes along with their clinicopathological features, and is approved for treatment recommendations for adjuvant hormonal therapy and chemotherapy in hormone-receptor-positive early breast cancer. The Prosigna test utilizes RNA extracted from macrodissected tumor cells obtained from formalin-fixed, paraffin-embedded (FFPE) tissue sections. However, RNA extracted from fresh-frozen (FF) bulk tissue without macrodissection is widely used for research purposes, and yields high-quality RNA for downstream analyses. To investigate the impact of the sample preparation approach on ROR scores, we analyzed 94 breast carcinomas included in an observational study that had available gene expression data from macrodissected FFPE tissue and FF bulk tumor tissue, along with the clinically approved Prosigna scores for the node-negative, hormone-receptor-positive, HER2-negative cases (n = 54). ROR scores were calculated in R; the resulting two sets of scores from FFPE and FF samples were compared, and treatment recommendations were evaluated. Overall, ROR scores calculated based on the macrodissected FFPE tissue were consistent with the Prosigna scores. However, analyses from bulk tissue yielded a higher proportion of cases classified as normal-like; these were samples with relatively low tumor cellularity, leading to lower ROR scores. When comparing ROR scores (low, intermediate, and high), discordant cases between the two preparation approaches were revealed among the luminal tumors; the recommended treatment would have changed in a minority of cases.

Tailored graphical lasso for data integration in gene network reconstruction.

BACKGROUND: Identifying gene interactions is a topic of great importance in genomics, and approaches based on network models provide a powerful tool for studying these. Assuming a Gaussian graphical model, a gene association network may be estimated from multiomic data based on the non-zero entries of the inverse covariance matrix. Inferring such biological networks is challenging because of the high dimensionality of the problem, making traditional estimators unsuitable. The graphical lasso is constructed for the estimation of sparse inverse covariance matrices in such situations, using [Formula: see text]-penalization on the matrix entries. The weighted graphical lasso is an extension in which prior biological information from other sources is integrated into the model. There are however issues with this approach, as it naïvely forces the prior information into the network estimation, even if it is misleading or does not agree with the data at hand. Further, if an associated network based on other data is used as the prior, the method often fails to utilize the information effectively. RESULTS: We propose a novel graphical lasso approach, the tailored graphical lasso, that aims to handle prior information of unknown accuracy more effectively. We provide an R package implementing the method, tailoredGlasso. Applying the method to both simulated and real multiomic data sets, we find that it outperforms the unweighted and weighted graphical lasso in terms of all performance measures we consider. In fact, the graphical lasso and weighted graphical lasso can be considered special cases of the tailored graphical lasso, and a parameter determined by the data measures the usefulness of the prior information. We also find that among a larger set of methods, the tailored graphical is the most suitable for network inference from high-dimensional data with prior information of unknown accuracy. With our method, mRNA data are demonstrated to provide highly useful prior information for protein-protein interaction networks. CONCLUSIONS: The method we introduce utilizes useful prior information more effectively without involving any risk of loss of accuracy should the prior information be misleading.

Contrasting DCIS and invasive breast cancer by subtype suggests basal-like DCIS as distinct lesions.

Ductal carcinoma in situ (DCIS) is a non-invasive type of breast cancer with highly variable potential of becoming invasive and affecting mortality. Currently, many patients with DCIS are overtreated due to the lack of specific biomarkers that distinguish low risk lesions from those with a higher risk of progression. In this study, we analyzed 57 pure DCIS and 313 invasive breast cancers (IBC) from different patients. Three levels of genomic data were obtained; gene expression, DNA methylation, and DNA copy number. We performed subtype stratified analyses and identified key differences between DCIS and IBC that suggest subtype specific progression. Prominent differences were found in tumors of the basal-like subtype: Basal-like DCIS were less proliferative and showed a higher degree of differentiation than basal-like IBC. Also, core basal tumors (characterized by high correlation to the basal-like centroid) were not identified amongst DCIS as opposed to IBC. At the copy number level, basal-like DCIS exhibited fewer copy number aberrations compared with basal-like IBC. An intriguing finding through analysis of the methylome was hypermethylation of multiple protocadherin genes in basal-like IBC compared with basal-like DCIS and normal tissue, possibly caused by long range epigenetic silencing. This points to silencing of cell adhesion-related genes specifically in IBC of the basal-like subtype. Our work confirms that subtype stratification is essential when studying progression from DCIS to IBC, and we provide evidence that basal-like DCIS show less aggressive characteristics and question the assumption that basal-like DCIS is a direct precursor of basal-like invasive breast cancer.

Prognostic value of PAM50 and risk of recurrence score in patients with early-stage breast cancer with long-term follow-up.

BACKGROUND: The aim of this study was to investigate the prognostic value of the PAM50 intrinsic subtypes and risk of recurrence (ROR) score in patients with early breast cancer and long-term follow-up. A special focus was placed on hormone receptor-positive/human epidermal growth factor receptor 2-negative (HR+/HER2-) pN0 patients not treated with chemotherapy. METHODS: Patients with early breast cancer (n = 653) enrolled in the observational Oslo1 study (1995-1998) were followed for distant recurrence and breast cancer death. Clinicopathological parameters were collected from hospital records. The primary tumors were analyzed using the Prosigna® PAM50 assay to determine the prognostic value of the intrinsic subtypes and ROR score in comparison with pathological characteristics. The primary endpoints were distant disease-free survival (DDFS) and breast cancer-specific survival (BCSS). RESULTS: Of 653 tumors, 52.2% were classified as luminal A, 26.5% as luminal B, 10.6% as HER2-enriched, and 10.7% as basal-like. Among the HR+/HER2- patients (n = 476), 37.8% were categorized as low risk by ROR score, 22.7% as intermediate risk, and 39.5% as high risk. Median follow-up durations for BCSS and DDFS were 16.6 and 7.1 years, respectively. Multivariate analysis showed that intrinsic subtypes (all patients) and ROR risk classification (HR+/HER2- patients) yielded strong prognostic information. Among the HR+/HER2- pN0 patients with no adjuvant treatment (n = 231), 53.7% of patients had a low ROR, and their prognosis at 15 years was excellent (15-year BCSS 96.3%). Patients with intermediate risk had reduced survival compared with those with low risk (p = 0.005). In contrast, no difference in survival between the low- and intermediate-risk groups was seen for HR+/HER2- pN0 patients who received tamoxifen only. Ki-67 protein, grade, and ROR score were analyzed in the unselected, untreated pT1pN0 HR+/HER2- population (n = 171). In multivariate analysis, ROR score outperformed both Ki-67 and grade. Furthermore, 55% of patients who according to the PREDICT tool ( http://www.predict.nhs.uk/ ) would be considered chemotherapy candidates were ROR low risk (33%) or luminal A ROR intermediate risk (22%). CONCLUSIONS: The PAM50 intrinsic subtype classification and ROR score improve classification of patients with breast cancer into prognostic groups, allowing for a more precise identification of future recurrence risk and providing an improved basis for adjuvant treatment decisions. Node-negative patients with low ROR scores had an excellent outcome at 15 years even in the absence of adjuvant therapy.

Better prediction by use of co-data: adaptive group-regularized ridge regression.

For many high-dimensional studies, additional information on the variables, like (genomic) annotation or external p-values, is available. In the context of binary and continuous prediction, we develop a method for adaptive group-regularized (logistic) ridge regression, which makes structural use of such 'co-data'. Here, 'groups' refer to a partition of the variables according to the co-data. We derive empirical Bayes estimates of group-specific penalties, which possess several nice properties: (i) They are analytical. (ii) They adapt to the informativeness of the co-data for the data at hand. (iii) Only one global penalty parameter requires tuning by cross-validation. In addition, the method allows use of multiple types of co-data at little extra computational effort. We show that the group-specific penalties may lead to a larger distinction between 'near-zero' and relatively large regression parameters, which facilitates post hoc variable selection. The method, termed GRridge, is implemented in an easy-to-use R-package. It is demonstrated on two cancer genomics studies, which both concern the discrimination of precancerous cervical lesions from normal cervix tissues using methylation microarray data. For both examples, GRridge clearly improves the predictive performances of ordinary logistic ridge regression and the group lasso. In addition, we show that for the second study, the relatively good predictive performance is maintained when selecting only 42 variables.

Handling realistic assumptions in hypothesis testing of 3D co-localization of genomic elements.

The study of chromatin 3D structure has recently gained much focus owing to novel techniques for detecting genome-wide chromatin contacts using next-generation sequencing. A deeper understanding of the architecture of the DNA inside the nucleus is crucial for gaining insight into fundamental processes such as transcriptional regulation, genome dynamics and genome stability. Chromatin conformation capture-based methods, such as Hi-C and ChIA-PET, are now paving the way for routine genome-wide studies of chromatin 3D structure in a range of organisms and tissues. However, appropriate methods for analyzing such data are lacking. Here, we propose a hypothesis test and an enrichment score of 3D co-localization of genomic elements that handles intra- or interchromosomal interactions, both separately and jointly, and that adjusts for biases caused by structural dependencies in the 3D data. We show that maintaining structural properties during resampling is essential to obtain valid estimation of P-values. We apply the method on chromatin states and a set of mutated regions in leukemia cells, and find significant co-localization of these elements, with varying enrichment scores, supporting the role of chromatin 3D structure in shaping the landscape of somatic mutations in cancer.