Person:

Beck, Andrew

Leiomyosarcoma (LMS) is a malignant, soft-tissue tumor for which few effective therapies exist. Previously, we showed that there are three molecular subtypes of LMS. Here, we analyzed genes differentially expressed in each of the three LMS subtypes as compared to benign leiomyomas and then used the Connectivity Map (cmap) to calculate enrichment scores for the 1309 cmap drugs in order to identify candidate molecules with the potential to induce a benign, leiomyoma-like phenotype in LMS cells. 11 drugs were selected and tested for their ability to inhibit the growth of three human LMS cell lines. We identified two drugs with in vitro efficacy against LMS, one of which had a strongly negative enrichment score (Cantharidin) and the other of which had a strongly positive enrichment score (MG-132). Given MG-132's strong inhibitory effect on LMS cell viability, we hypothesized that LMS cells may be sensitive to treatment with other proteasome inhibitors and demonstrated that bortezomib, a clinically-approved proteasome inhibitor not included in the original cmap screen, potently inhibited the viability of the LMS cell lines. These findings suggest that systematically linking LMS subtype-specific expression signatures with drug-associated expression profiles represents a promising approach for the identification of new drugs for LMS.

Introduction: Screening mammography has contributed to a significant increase in the diagnosis of ductal carcinoma in situ (DCIS), raising concerns about overdiagnosis and overtreatment. Building on prior observations from lineage evolution analysis, we examined whether measuring genomic features of DCIS would predict association with invasive breast carcinoma (IBC). The long-term goal is to enhance standard clinicopathologic measures of low- versus high-risk DCIS and to enable risk-appropriate treatment. Methods: We studied three common chromosomal copy number alterations (CNA) in IBC and designed fluorescence in situ hybridization-based assay to measure copy number at these loci in DCIS samples. Clinicopathologic data were extracted from the electronic medical records of Stanford Cancer Institute and linked to demographic data from the population-based California Cancer Registry; results were integrated with data from tissue microarrays of specimens containing DCIS that did not develop IBC versus DCIS with concurrent IBC. Multivariable logistic regression analysis was performed to describe associations of CNAs with these two groups of DCIS. Results: We examined 271 patients with DCIS (120 that did not develop IBC and 151 with concurrent IBC) for the presence of 1q, 8q24 and 11q13 copy number gains. Compared to DCIS-only patients, patients with concurrent IBC had higher frequencies of CNAs in their DCIS samples. On multivariable analysis with conventional clinicopathologic features, the copy number gains were significantly associated with concurrent IBC. The state of two of the three copy number gains in DCIS was associated with a risk of IBC that was 9.07 times that of no copy number gains, and the presence of gains at all three genomic loci in DCIS was associated with a more than 17-fold risk (P = 0.0013). Conclusions: CNAs have the potential to improve the identification of high-risk DCIS, defined by presence of concurrent IBC. Expanding and validating this approach in both additional cross-sectional and longitudinal cohorts may enable improved risk stratification and risk-appropriate treatment in DCIS. Electronic supplementary material The online version of this article (doi:10.1186/s13058-015-0623-y) contains supplementary material, which is available to authorized users.

The categorization of intraductal proliferative lesions of the breast based on routine light microscopic examination of histopathologic sections is in many cases challenging, even for experienced pathologists. The development of computational tools to aid pathologists in the characterization of these lesions would have great diagnostic and clinical value. As a first step to address this issue, we evaluated the ability of computational image analysis to accurately classify DCIS and UDH and to stratify nuclear grade within DCIS. Using 116 breast biopsies diagnosed as DCIS or UDH from the Massachusetts General Hospital (MGH), we developed a computational method to extract 392 features corresponding to the mean and standard deviation in nuclear size and shape, intensity, and texture across 8 color channels. We used L1-regularized logistic regression to build classification models to discriminate DCIS from UDH. The top-performing model contained 22 active features and achieved an AUC of 0.95 in cross-validation on the MGH data-set. We applied this model to an external validation set of 51 breast biopsies diagnosed as DCIS or UDH from the Beth Israel Deaconess Medical Center, and the model achieved an AUC of 0.86. The top-performing model contained active features from all color-spaces and from the three classes of features (morphology, intensity, and texture), suggesting the value of each for prediction. We built models to stratify grade within DCIS and obtained strong performance for stratifying low nuclear grade vs. high nuclear grade DCIS (AUC = 0.98 in cross-validation) with only moderate performance for discriminating low nuclear grade vs. intermediate nuclear grade and intermediate nuclear grade vs. high nuclear grade DCIS (AUC = 0.83 and 0.69, respectively). These data show that computational pathology models can robustly discriminate benign from malignant intraductal proliferative lesions of the breast and may aid pathologists in the diagnosis and classification of these lesions.

Gene set enrichment analysis (GSEA) associates gene sets and phenotypes, its use is predicated on the choice of a pre-defined collection of sets. The defacto standard implementation of GSEA provides seven collections yet there are no guidelines for the choice of collections and the impact of such choice, if any, is unknown. Here we compare each of the standard gene set collections in the context of a large dataset of drug response in human cancer cell lines. We define and test a new collection based on gene co-expression in cancer cell lines to compare the performance of the standard collections to an externally derived cell line based collection. The results show that GSEA findings vary significantly depending on the collection chosen for analysis. Henceforth, collections should be carefully selected and reported in studies that leverage GSEA.

Background: Molecular characterization of tumors has been critical for identifying important genes in cancer biology and for improving tumor classification and diagnosis. Long non-coding RNAs, as a new, relatively unstudied class of transcripts, provide a rich opportunity to identify both functional drivers and cancer-type-specific biomarkers. However, despite the potential importance of long non-coding RNAs to the cancer field, no comprehensive survey of long non-coding RNA expression across various cancers has been reported. Results: We performed a sequencing-based transcriptional survey of both known long non-coding RNAs and novel intergenic transcripts across a panel of 64 archival tumor samples comprising 17 diagnostic subtypes of adenocarcinomas, squamous cell carcinomas and sarcomas. We identified hundreds of transcripts from among the known 1,065 long non-coding RNAs surveyed that showed variability in transcript levels between the tumor types and are therefore potential biomarker candidates. We discovered 1,071 novel intergenic transcribed regions and demonstrate that these show similar patterns of variability between tumor types. We found that many of these differentially expressed cancer transcripts are also expressed in normal tissues. One such novel transcript specifically expressed in breast tissue was further evaluated using RNA in situ hybridization on a panel of breast tumors. It was shown to correlate with low tumor grade and estrogen receptor expression, thereby representing a potentially important new breast cancer biomarker. Conclusions: This study provides the first large survey of long non-coding RNA expression within a panel of solid cancers and also identifies a number of novel transcribed regions differentially expressed across distinct cancer types that represent candidate biomarkers for future research.

Advancements in molecular biology have unveiled multiple breast cancer promoting pathways and potential therapeutic targets. Large randomized clinical trials remain the ultimate means of validating therapeutic efficacy, but they require large cohorts of patients and are lengthy and costly. A useful approach is to conduct a window of opportunity study in which patients are exposed to a drug pre-surgically during the interval between the core needle biopsy and the definitive surgery. These are non-therapeutic studies and the end point is not clinical or pathological response but rather evaluation of molecular changes in the tumor specimens that can predict response. However, since the end points of the non-therapeutic studies are biologic, it is critical to first define the biologic changes that occur in the absence of treatment. In this study, we compared the molecular profiles of breast cancer tumors at the time of the diagnostic biopsy versus the definitive surgery in the absence of any intervention using the Nanostring nCounter platform. We found that while the majority of the transcripts did not vary between the two biopsies, there was evidence of activation of immune related genes in response to the first biopsy and further investigations of the immune changes after a biopsy in early breast cancer seem warranted.

Mammographic density reflects the amount of stromal and epithelial tissues in relation to adipose tissue in the breast and is a strong risk factor for breast cancer. Here we report the results from meta-analysis of genome-wide association studies (GWAS) of three mammographic density phenotypes: dense area, non-dense area and percent density in up to 7,916 women in stage 1 and an additional 10,379 women in stage 2. We identify genome-wide significant (P<5×10−8) loci for dense area (AREG, ESR1, ZNF365, LSP1/TNNT3, IGF1, TMEM184B, SGSM3/MKL1), non-dense area (8p11.23) and percent density (PRDM6, 8p11.23, TMEM184B). Four of these regions are known breast cancer susceptibility loci, and four additional regions were found to be associated with breast cancer (P<0.05) in a large meta-analysis. These results provide further evidence of a shared genetic basis between mammographic density and breast cancer and illustrate the power of studying intermediate quantitative phenotypes to identify putative disease susceptibility loci.

We investigate 71 single nucleotide polymorphisms (SNPs) identified in meta-analytic studies of genome-wide association studies (GWAS) of breast cancer, the majority of which are located in intergenic or intronic regions. To explore regulatory impacts of these variants we conducted expression quantitative loci (eQTL) analyses on tissue samples from 376 invasive postmenopausal breast cancer cases in the Nurses’ Health Study (NHS) diagnosed from 1990–2004. Expression analysis was conducted on all formalin-fixed paraffin-embedded (FFPE) tissue samples (and on 264 adjacent normal samples) using the Affymetrix Human Transcriptome Array. Significance and ranking of associations between tumor receptor status and expression variation was preserved between NHS FFPE and TCGA fresh-frozen sample sets (Spearman r = 0.85, p<10^-10 for 17 of the 21 Oncotype DX recurrence signature genes). At an FDR threshold of 10%, we identified 27 trans-eQTLs associated with expression variation in 217 distinct genes. SNP-gene associations can be explored using an open-source interactive browser distributed in a Bioconductor package. Using a new a procedure for testing hypotheses relating SNP content to expression patterns in gene sets, defined as molecular function pathways, we find that loci on 6q14 and 6q25 affect various gene sets and molecular pathways (FDR < 10%). Although the ultimate biological interpretation of the GWAS-identified variants remains to be uncovered, this study validates the utility of expression analysis of this FFPE expression set for more detailed integrative analyses.

Long non-coding RNAs (lncRNAs) have been implicated in normal cellular homeostasis as well as pathophysiological conditions, including cancer. Here we performed global gene expression profiling of mammary epithelial cells transformed by oncogenic v-Src, and identified a large subset of uncharacterized lncRNAs potentially involved in breast cancer development. Specifically, our analysis revealed a novel lncRNA, LINC00520 that is upregulated upon ectopic expression of oncogenic v-Src, in a manner that is dependent on the transcription factor STAT3. Similarly, LINC00520 is also increased in mammary epithelial cells transformed by oncogenic PI3K and its expression is decreased upon knockdown of mutant PIK3CA. Additional expression profiling highlight that LINC00520 is elevated in a subset of human breast carcinomas, with preferential enrichment in the basal-like molecular subtype. ShRNA-mediated depletion of LINC00520 results in decreased cell migration and loss of invasive structures in 3D. RNA sequencing analysis uncovers several genes that are differentially expressed upon ectopic expression of LINC00520, a significant subset of which are also induced in v-Src-transformed MCF10A cells. Together, these findings characterize LINC00520 as a lncRNA that is regulated by oncogenic Src, PIK3CA and STAT3, and which may contribute to the molecular etiology of breast cancer.

Background: Enhancer of zeste homolog 2 (EZH2) is a polycomb-group protein that is involved in stem cell renewal and carcinogenesis. In breast cancer, increased EZH2 expression is associated with aggressiveness and has been suggested to identify normal breast epithelium at increased risk of breast cancer development. However, the association between EZH2 expression in benign breast tissue and breast cancer risk has not previously been evaluated in a large prospective cohort. Methods: We examined the association between EZH2 protein expression and subsequent breast cancer risk using logistic regression in a nested case-control study of benign breast disease (BBD) and breast cancer within the Nurses’ Health Studies. EZH2 immunohistochemical expression in normal breast epithelium and stroma was evaluated by computational image analysis and its association with breast cancer risk was analyzed after adjusting for matching factors between cases and controls, the concomitant BBD diagnosis, and the Ki67 proliferation index. Results: Women with a breast biopsy in which more than 20% of normal epithelial cells expressed EZH2 had a significantly increased risk of developing breast cancer (odds ratio (OR) 2.95, 95% confidence interval (CI) 1.11–7.84) compared to women with less than 10% EZH2 epithelial expression. The risk of developing breast cancer increased for each 5% increase in EZH2 expression (OR 1.22, 95% CI 1.02–1.46, p value 0.026). Additionally, women with high EZH2 expression and low estrogen receptor (ER) expression had a 4-fold higher risk of breast cancer compared to women with low EZH2 and low ER expression (OR 4.02, 95% CI 1.29–12.59). Conclusions: These results provide further evidence that EZH2 expression in the normal breast epithelium is independently associated with breast cancer risk and might be used to assist in risk stratification for women with benign breast biopsies. Electronic supplementary material The online version of this article (doi:10.1186/s13058-017-0817-6) contains supplementary material, which is available to authorized users.