Image-guided Coring for Large-scale Studies in Molecular Pathology.

Sampling of formalin-fixed paraffin-embedded (FFPE) tissue blocks is a critical initial step in molecular pathology. Image-guided coring (IGC) is a new method for using digital pathology images to guide tissue block coring for molecular analyses. The goal of our study is to evaluate the use of IGC for both tissue-based and nucleic acid-based projects in molecular pathology. First, we used IGC to construct a tissue microarray (TMA); second, we used IGC for FFPE block sampling followed by RNA extraction; and third, we assessed the correlation between nuclear counts quantitated from the IGC images and RNA yields. We used IGC to construct a TMA containing 198 normal and breast cancer cores. Histopathologic analysis showed high accuracy for obtaining tumor and normal breast tissue. Next, we used IGC to obtain normal and tumor breast samples before RNA extraction. We selected a random subset of tumor and normal samples to perform computational image analysis to quantify nuclear density, and we built regression models to estimate RNA yields from nuclear count, age of the block, and core diameter. Number of nuclei and core diameter were the strongest predictors of RNA yields in both normal and tumor tissue. IGC is an effective method for sampling FFPE tissue blocks for TMA construction and nucleic acid extraction. We identify significant associations between quantitative nuclear counts obtained from IGC images and RNA yields, suggesting that the integration of computational image analysis with IGC may be an effective approach for tumor sampling in large-scale molecular studies.

Extensive rewiring of epithelial-stromal co-expression networks in breast cancer.

BACKGROUND: Epithelial-stromal crosstalk plays a critical role in invasive breast cancer pathogenesis; however, little is known on a systems level about how epithelial-stromal interactions evolve during carcinogenesis. RESULTS: We develop a framework for building genome-wide epithelial-stromal co-expression networks composed of pairwise co-expression relationships between mRNA levels of genes expressed in the epithelium and stroma across a population of patients. We apply this method to laser capture micro-dissection expression profiling datasets in the setting of breast carcinogenesis. Our analysis shows that epithelial-stromal co-expression networks undergo extensive rewiring during carcinogenesis, with the emergence of distinct network hubs in normal breast, and estrogen receptor-positive and estrogen receptor-negative invasive breast cancer, and the emergence of distinct patterns of functional network enrichment. In contrast to normal breast, the strongest epithelial-stromal co-expression relationships in invasive breast cancer mostly represent self-loops, in which the same gene is co-expressed in epithelial and stromal regions. We validate this observation using an independent laser capture micro-dissection dataset and confirm that self-loop interactions are significantly increased in cancer by performing computational image analysis of epithelial and stromal protein expression using images from the Human Protein Atlas. CONCLUSIONS: Epithelial-stromal co-expression network analysis represents a new approach for systems-level analyses of spatially localized transcriptomic data. The analysis provides new biological insights into the rewiring of epithelial-stromal co-expression networks and the emergence of epithelial-stromal co-expression self-loops in breast cancer. The approach may facilitate the development of new diagnostics and therapeutics targeting epithelial-stromal interactions in cancer.

Computational pathology to discriminate benign from malignant intraductal proliferations of the breast.

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.

Genome-wide association study identifies multiple loci associated with both mammographic density and breast cancer risk.

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 and SGSM3/MKL1), non-dense area (8p11.23) and percent density (PRDM6, 8p11.23 and 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.

The adherens junction protein afadin is an AKT substrate that regulates breast cancer cell migration.

UNLABELLED: The PI3K-AKT signaling pathway regulates all phenotypes that contribute to progression of human cancers, including breast cancer. AKT mediates signal relay by phosphorylating numerous substrates, which are causally implicated in biologic responses such as cell growth, survival, metabolic reprogramming, migration, and invasion. Here a new AKT substrate is identified, the adherens junction protein Afadin, which is phosphorylated by AKT at Ser1718. Importantly, under conditions of physiologic IGF-1 signaling and oncogenic PI3K and AKT, Afadin is phosphorylated by all AKT isoforms, and this phosphorylation elicits a relocalization of Afadin from adherens junctions to the nucleus. Also, phosphorylation of Afadin increased breast cancer cell migration that was dependent on Ser1718 phosphorylation. Finally, nuclear localization of Afadin was observed in clinical breast cancer specimens, indicating that regulation of Afadin by the PI3K-AKT pathway has pathophysiologic significance. IMPLICATIONS: Phosphorylation of the adhesion protein Afadin by AKT downstream of the PI3K pathway, leads to redistribution of Afadin and controls cancer cell migration.

Estrogen receptor negative/progesterone receptor positive breast cancer is not a reproducible subtype.

INTRODUCTION: Estrogen receptor (ER) and progesterone receptor (PR) testing are performed in the evaluation of breast cancer. While the clinical utility of ER as a predictive biomarker to identify patients likely to benefit from hormonal therapy is well-established, the added value of PR is less well-defined. The primary goals of our study were to assess the distribution, inter-assay reproducibility, and prognostic significance of breast cancer subtypes defined by patterns of ER and PR expression. METHODS: We integrated gene expression microarray (GEM) and clinico-pathologic data from 20 published studies to determine the frequency (n = 4,111) and inter-assay reproducibility (n = 1,752) of ER/PR subtypes (ER+/PR+, ER+/PR-, ER-/PR-, ER-/PR+). To extend our findings, we utilized a cohort of patients from the Nurses' Health Study (NHS) with ER/PR data recorded in the medical record and assessed on tissue microarrays (n = 2,011). In both datasets, we assessed the association of ER and PR expression with survival. RESULTS: In a genome-wide analysis, progesterone receptor was among the least variable genes in ER- breast cancer. The ER-/PR+ subtype was rare (approximately 1 to 4%) and showed no significant reproducibility (Kappa = 0.02 and 0.06, in the GEM and NHS datasets, respectively). The vast majority of patients classified as ER-/PR+ in the medical record (97% and 94%, in the GEM and NHS datasets) were re-classified by a second method. In the GEM dataset (n = 2,731), progesterone receptor mRNA expression was associated with prognosis in ER+ breast cancer (adjusted P <0.001), but not in ER- breast cancer (adjusted P = 0.21). PR protein expression did not contribute significant prognostic information to multivariate models considering ER and other standard clinico-pathologic features in the GEM or NHS datasets. CONCLUSION: ER-/PR+ breast cancer is not a reproducible subtype. PR expression is not associated with prognosis in ER- breast cancer, and PR does not contribute significant independent prognostic information to multivariate models considering ER and other standard clinico-pathologic factors. Given that PR provides no clinically actionable information in ER+ breast cancer, these findings question the utility of routine PR testing in breast cancer.

Functional variants at the 11q13 risk locus for breast cancer regulate cyclin D1 expression through long-range enhancers.

Analysis of 4,405 variants in 89,050 European subjects from 41 case-control studies identified three independent association signals for estrogen-receptor-positive tumors at 11q13. The strongest signal maps to a transcriptional enhancer element in which the G allele of the best candidate causative variant rs554219 increases risk of breast cancer, reduces both binding of ELK4 transcription factor and luciferase activity in reporter assays, and may be associated with low cyclin D1 protein levels in tumors. Another candidate variant, rs78540526, lies in the same enhancer element. Risk association signal 2, rs75915166, creates a GATA3 binding site within a silencer element. Chromatin conformation studies demonstrate that these enhancer and silencer elements interact with each other and with their likely target gene, CCND1.

Significance analysis of prognostic signatures.

A major goal in translational cancer research is to identify biological signatures driving cancer progression and metastasis. A common technique applied in genomics research is to cluster patients using gene expression data from a candidate prognostic gene set, and if the resulting clusters show statistically significant outcome stratification, to associate the gene set with prognosis, suggesting its biological and clinical importance. Recent work has questioned the validity of this approach by showing in several breast cancer data sets that "random" gene sets tend to cluster patients into prognostically variable subgroups. This work suggests that new rigorous statistical methods are needed to identify biologically informative prognostic gene sets. To address this problem, we developed Significance Analysis of Prognostic Signatures (SAPS) which integrates standard prognostic tests with a new prognostic significance test based on stratifying patients into prognostic subtypes with random gene sets. SAPS ensures that a significant gene set is not only able to stratify patients into prognostically variable groups, but is also enriched for genes showing strong univariate associations with patient prognosis, and performs significantly better than random gene sets. We use SAPS to perform a large meta-analysis (the largest completed to date) of prognostic pathways in breast and ovarian cancer and their molecular subtypes. Our analyses show that only a small subset of the gene sets found statistically significant using standard measures achieve significance by SAPS. We identify new prognostic signatures in breast and ovarian cancer and their corresponding molecular subtypes, and we show that prognostic signatures in ER negative breast cancer are more similar to prognostic signatures in ovarian cancer than to prognostic signatures in ER positive breast cancer. SAPS is a powerful new method for deriving robust prognostic biological signatures from clinically annotated genomic datasets.

Sex differences in social interaction behavior following social defeat stress in the monogamous California mouse (Peromyscus californicus).

Stressful life experiences are known to be a precipitating factor for many mental disorders. The social defeat model induces behavioral responses in rodents (e.g. reduced social interaction) that are similar to behavioral patterns associated with mood disorders. The model has contributed to the discovery of novel mechanisms regulating behavioral responses to stress, but its utility has been largely limited to males. This is disadvantageous because most mood disorders have a higher incidence in women versus men. Male and female California mice (Peromyscus californicus) aggressively defend territories, which allowed us to observe the effects of social defeat in both sexes. In two experiments, mice were exposed to three social defeat or control episodes. Mice were then behaviorally phenotyped, and indirect markers of brain activity and corticosterone responses to a novel social stimulus were assessed. Sex differences in behavioral responses to social stress were long lasting (4 wks). Social defeat reduced social interaction responses in females but not males. In females, social defeat induced an increase in the number of phosphorylated CREB positive cells in the nucleus accumbens shell after exposure to a novel social stimulus. This effect of defeat was not observed in males. The effects of defeat in females were limited to social contexts, as there were no differences in exploratory behavior in the open field or light-dark box test. These data suggest that California mice could be a useful model for studying sex differences in behavioral responses to stress, particularly in neurobiological mechanisms that are involved with the regulation of social behavior.