Comprehensive Molecular Portraits of Invasive Lobular Breast Cancer.

Invasive lobular carcinoma (ILC) is the second most prevalent histologic subtype of invasive breast cancer. Here, we comprehensively profiled 817 breast tumors, including 127 ILC, 490 ductal (IDC), and 88 mixed IDC/ILC. Besides E-cadherin loss, the best known ILC genetic hallmark, we identified mutations targeting PTEN, TBX3, and FOXA1 as ILC enriched features. PTEN loss associated with increased AKT phosphorylation, which was highest in ILC among all breast cancer subtypes. Spatially clustered FOXA1 mutations correlated with increased FOXA1 expression and activity. Conversely, GATA3 mutations and high expression characterized luminal A IDC, suggesting differential modulation of ER activity in ILC and IDC. Proliferation and immune-related signatures determined three ILC transcriptional subtypes associated with survival differences. Mixed IDC/ILC cases were molecularly classified as ILC-like and IDC-like revealing no true hybrid features. This multidimensional molecular atlas sheds new light on the genetic bases of ILC and provides potential clinical options.

Immune Landscape of Viral- and Carcinogen-Driven Head and Neck Cancer.

Head and neck squamous cell carcinoma (HNSCC) arises through exposure to environmental carcinogens or malignant transformation by human papillomavirus (HPV). Here, we assessed the transcriptional profiles of 131,224 single cells from peripheral and intra-tumoral immune populations from patients with HPV- and HPV+ HNSCC and healthy donors. Immune cells within tumors of HPV- and HPV+ HNSCC displayed a spectrum of transcriptional signatures, with helper CD4+ T cells and B cells being relatively divergent and CD8+ T cells and CD4+ regulatory T cells being relatively similar. Transcriptional results were contextualized through multispectral immunofluorescence analyses and evaluating putative cell-cell communication based on spatial proximity. These analyses defined a gene expression signature associated with CD4+ T follicular helper cells that is associated with longer progression-free survival in HNSCC patients. The datasets and analytical approaches herein provide a resource for the further study of the impact of immune cells on viral- and carcinogen-induced cancers.

Transcriptomic congruence analysis for evaluating model organisms.

Model organisms are instrumental substitutes for human studies to expedite basic, translational, and clinical research. Despite their indispensable role in mechanistic investigation and drug development, molecular congruence of animal models to humans has long been questioned and debated. Little effort has been made for an objective quantification and mechanistic exploration of a model organism's resemblance to humans in terms of molecular response under disease or drug treatment. We hereby propose a framework, namely Congruence Analysis for Model Organisms (CAMO), for transcriptomic response analysis by developing threshold-free differential expression analysis, quantitative concordance/discordance scores incorporating data variabilities, pathway-centric downstream investigation, knowledge retrieval by text mining, and topological gene module detection for hypothesis generation. Instead of a genome-wide vague and dichotomous answer of "poorly" or "greatly" mimicking humans, CAMO assists researchers to numerically quantify congruence, to dissect true cross-species differences from unwanted biological or cohort variabilities, and to visually identify molecular mechanisms and pathway subnetworks that are best or least mimicked by model organisms, which altogether provides foundations for hypothesis generation and subsequent translational decisions.

Systems approach for congruence and selection of cancer models towards precision medicine.

Cancer models are instrumental as a substitute for human studies and to expedite basic, translational, and clinical cancer research. For a given cancer type, a wide selection of models, such as cell lines, patient-derived xenografts, organoids and genetically modified murine models, are often available to researchers. However, how to quantify their congruence to human tumors and to select the most appropriate cancer model is a largely unsolved issue. Here, we present Congruence Analysis and Selection of CAncer Models (CASCAM), a statistical and machine learning framework for authenticating and selecting the most representative cancer models in a pathway-specific manner using transcriptomic data. CASCAM provides harmonization between human tumor and cancer model omics data, systematic congruence quantification, and pathway-based topological visualization to determine the most appropriate cancer model selection. The systems approach is presented using invasive lobular breast carcinoma (ILC) subtype and suggesting CAMA1 followed by UACC3133 as the most representative cell lines for ILC research. Two additional case studies for triple negative breast cancer (TNBC) and patient-derived xenograft/organoid (PDX/PDO) are further investigated. CASCAM is generalizable to any cancer subtype and will authenticate cancer models for faithful non-human preclinical research towards precision medicine.

Research autopsy programmes in oncology: shared experience from 14 centres across the world.

While there is a great clinical need to understand the biology of metastatic cancer in order to treat it more effectively, research is hampered by limited sample availability. Research autopsy programmes can crucially advance the field through synchronous, extensive, and high-volume sample collection. However, it remains an underused strategy in translational research. Via an extensive questionnaire, we collected information on the study design, enrolment strategy, study conduct, sample and data management, and challenges and opportunities of research autopsy programmes in oncology worldwide. Fourteen programmes participated in this study. Eight programmes operated 24 h/7 days, resulting in a lower median postmortem interval (time between death and start of the autopsy, 4 h) compared with those operating during working hours (9 h). Most programmes (n = 10) succeeded in collecting all samples within a median of 12 h after death. A large number of tumour sites were sampled during each autopsy (median 15.5 per patient). The median number of samples collected per patient was 58, including different processing methods for tumour samples but also non-tumour tissues and liquid biopsies. Unique biological insights derived from these samples included metastatic progression, treatment resistance, disease heterogeneity, tumour dormancy, interactions with the tumour micro-environment, and tumour representation in liquid biopsies. Tumour patient-derived xenograft (PDX) or organoid (PDO) models were additionally established, allowing for drug discovery and treatment sensitivity assays. Apart from the opportunities and achievements, we also present the challenges related with postmortem sample collections and strategies to overcome them, based on the shared experience of these 14 programmes. Through this work, we hope to increase the transparency of postmortem tissue donation, to encourage and aid the creation of new programmes, and to foster collaborations on these unique sample collections. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Liver tropism of ER mutant breast cancer is characterized by unique molecular changes and immune infiltration.

PURPOSE: Hotspot estrogen receptor alpha (ER/ESR1) mutations are recognized as the driver for both endocrine resistance and metastasis in advanced ER-positive (ER+) breast cancer, but their contributions to metastatic organ tropism remain insufficiently understood. In this study, we aim to comprehensively profile the organotropic metastatic pattern for ESR1 mutant breast cancer. METHODS: The organ-specific metastatic pattern of ESR1 mutant breast cancer was delineated using multi-omics data from multiple publicly available cohorts of ER+ metastatic breast cancer patients. Gene mutation/copy number variation (CNV) and differential gene expression analyses were performed to identify the genomic and transcriptomic alterations uniquely associated with ESR1 mutant liver metastasis. Upstream regulator, downstream pathway, and immune infiltration analysis were conducted for subsequent mechanistic investigations. RESULTS: ESR1 mutation-driven liver tropism was revealed by significant differences, encompassing a higher prevalence of liver metastasis in patients with ESR1 mutant breast cancer and an enrichment of mutations in liver metastatic samples. The significant enrichment of AGO2 copy number amplifications (CNAs) and multiple gene expression changes were revealed uniquely in ESR1 mutant liver metastasis. We also unveiled alterations in downstream signaling pathways and immune infiltration, particularly an enrichment of neutrophils, suggesting potential therapeutic vulnerabilities. CONCLUSION: Our data provide a comprehensive characterization of the behaviors and mechanisms of ESR1 mutant liver metastasis, paving the way for the development of personalized therapy to target liver metastasis for patients with ESR1 mutant breast cancer.

RET signaling in breast cancer therapeutic resistance and metastasis.

RET, a single-pass receptor tyrosine kinase encoded on human chromosome 10, is well known to the field of developmental biology for its role in the ontogenesis of the central and enteric nervous systems and the kidney. In adults, RET alterations have been characterized as drivers of non-small cell lung cancer and multiple neuroendocrine neoplasms. In breast cancer, RET signaling networks have been shown to influence diverse functions including tumor development, metastasis, and therapeutic resistance. While RET is known to drive the development and progression of multiple solid tumors, therapeutic agents selectively targeting RET are relatively new, though multiple multi-kinase inhibitors have shown promise as RET inhibitors in the past; further, RET has been historically neglected as a potential therapeutic co-target in endocrine-refractory breast cancers despite mounting evidence for a key pathologic role and repeated description of a bi-directional relationship with the estrogen receptor, the principal driver of most breast tumors. Additionally, the recent discovery of RET enrichment in breast cancer brain metastases suggests a role for RET inhibition specific to advanced disease. This review assesses the status of research on RET in breast cancer and evaluates the therapeutic potential of RET-selective kinase inhibitors across major breast cancer subtypes.

Mixed invasive ductal lobular carcinoma is clinically and pathologically more similar to invasive lobular than ductal carcinoma.

BACKGROUND: Mixed invasive ductal lobular carcinoma (mDLC) remains a poorly understood subtype of breast cancer composed of coexisting ductal and lobular components. METHODS: We sought to describe clinicopathologic characteristics and determine whether mDLC is clinically more similar to invasive ductal carcinoma (IDC) or invasive lobular carcinoma (ILC), using data from patients seen at the University of Pittsburgh Medical Center. RESULTS: We observed a higher concordance in clinicopathologic characteristics between mDLC and ILC, compared to IDC. There is a trend for higher rates of successful breast-conserving surgery after neoadjuvant chemotherapy in patients with mDLC compared to patients with ILC, in which it is known to be lower than in those with IDC. Metastatic patterns of mDLC demonstrate a propensity to develop in sites characteristic of both IDC and ILC. A meta-analysis evaluating mDLC showed shared features with both ILC and IDC with significantly more ER-positive and fewer high grades in mDLC compared to IDC, although mDLCs were significantly smaller and included fewer late-stage tumours compared to ILC. CONCLUSIONS: These findings support clinicopathologic characteristics of mDLC driven by individual ductal vs lobular components and given the dominance of lobular pathology, mDLC features are often more similar to ILC than IDC. This study exemplifies the complexity of mixed disease.

Semi-deconvolution of bulk and single-cell RNA-seq data with application to metastatic progression in breast cancer.

MOTIVATION: Identifying cell types and their abundances and how these evolve during tumor progression is critical to understanding the mechanisms of metastasis and identifying predictors of metastatic potential that can guide the development of new diagnostics or therapeutics. Single-cell RNA sequencing (scRNA-seq) has been especially promising in resolving heterogeneity of expression programs at the single-cell level, but is not always feasible, e.g. for large cohort studies or longitudinal analysis of archived samples. In such cases, clonal subpopulations may still be inferred via genomic deconvolution, but deconvolution methods have limited ability to resolve fine clonal structure and may require reference cell type profiles that are missing or imprecise. Prior methods can eliminate the need for reference profiles but show unstable performance when few bulk samples are available. RESULTS: In this work, we develop a new method using reference scRNA-seq to interpret sample collections for which only bulk RNA-seq is available for some samples, e.g. clonally resolving archived primary tissues using scRNA-seq from metastases. By integrating such information in a Quadratic Programming framework, our method can recover more accurate cell types and corresponding cell type abundances in bulk samples. Application to a breast tumor bone metastases dataset confirms the power of scRNA-seq data to improve cell type inference and quantification in same-patient bulk samples. AVAILABILITY AND IMPLEMENTATION: Source code is available on Github at https://github.com/CMUSchwartzLab/RADs.

Fusion-associated carcinomas of the breast: Diagnostic, prognostic, and therapeutic significance.

Recurrent gene fusions comprise a class of viable genetic targets in solid tumors that have culminated several recent breakthrough cancer therapies. Their role in breast cancer, however, remains largely underappreciated due to the complexity of genomic rearrangements in breast malignancy. Just recently, we and others have identified several recurrent gene fusions in breast cancer with important clinical and biological implications. Examples of the most significant recurrent gene fusions to date include (1) ESR1::CCDC170 gene fusions in luminal B and endocrine-resistant breast cancer that exert oncogenic function via modulating the HER2/HER3/SRC Proto-Oncogene (SRC) complex, (2) ESR1 exon 6 fusions in metastatic disease that drive estrogen-independent estrogen-receptor transcriptional activity, (3) BCL2L14::ETV6 fusions in a more aggressive form of the triple-negative subtype that prime epithelial-mesenchymal transition and endow paclitaxel resistance, (4) the ETV6::NTRK3 fusion in secretory breast carcinoma that constitutively activates NTRK3 kinase, (5) the oncogenic MYB-NFIB fusion as a genetic driver underpinning adenoid cystic carcinomas of the breast that activates MYB Proto-Oncogene (MYB) pathway, and (6) the NOTCH/microtubule-associated serine-threonine (MAST) kinase gene fusions that activate NOTCH and MAST signaling. Importantly, these fusions are enriched in more aggressive and lethal breast cancer presentations and appear to confer therapeutic resistance. Thus, these gene fusions could be utilized as genetic biomarkers to identify patients who require more intensive treatment and surveillance. In addition, kinase fusions are currently being evaluated in breast cancer clinical trials and ongoing mechanistic investigation is exposing therapeutic vulnerabilities in patients with fusion-positive disease.

Results of a worldwide survey on the currently used histopathological diagnostic criteria for invasive lobular breast cancer.

Invasive lobular carcinoma (ILC) represents the second most common subtype of breast cancer (BC), accounting for up to 15% of all invasive BC. Loss of cell adhesion due to functional inactivation of E-cadherin is the hallmark of ILC. Although the current world health organization (WHO) classification for diagnosing ILC requires the recognition of the dispersed or linear non-cohesive growth pattern, it is not mandatory to demonstrate E-cadherin loss by immunohistochemistry (IHC). Recent results of central pathology review of two large randomized clinical trials have demonstrated relative overdiagnosis of ILC, as only ~60% of the locally diagnosed ILCs were confirmed by central pathology. To understand the possible underlying reasons of this discrepancy, we undertook a worldwide survey on the current practice of diagnosing BC as ILC. A survey was drafted by a panel of pathologists and researchers from the European lobular breast cancer consortium (ELBCC) using the online tool SurveyMonkey®. Various parameters such as indications for IHC staining, IHC clones, and IHC staining procedures were questioned. Finally, systematic reporting of non-classical ILC variants were also interrogated. This survey was sent out to pathologists worldwide and circulated from December 14, 2020 until July, 1 2021. The results demonstrate that approximately half of the institutions use E-cadherin expression loss by IHC as an ancillary test to diagnose ILC and that there is a great variability in immunostaining protocols. This might cause different staining results and discordant interpretations. As ILC-specific therapeutic and diagnostic avenues are currently explored in the context of clinical trials, it is of importance to improve standardization of histopathologic diagnosis of ILC diagnosis.

Acquired mutations and transcriptional remodeling in long-term estrogen-deprived locoregional breast cancer recurrences.

BACKGROUND: Endocrine therapy resistance is a hallmark of advanced estrogen receptor (ER)-positive breast cancer. In this study, we aimed to determine acquired genomic changes in endocrine-resistant disease. METHODS: We performed DNA/RNA hybrid-capture sequencing on 12 locoregional recurrences after long-term estrogen deprivation and identified acquired genomic changes versus each tumor's matched primary. RESULTS: Despite being up to 7 years removed from the primary lesion, most recurrences harbored similar intrinsic transcriptional and copy number profiles. Only two genes, AKAP9 and KMT2C, were found to have single nucleotide variant (SNV) enrichments in more than one recurrence. Enriched mutations in single cases included SNVs within transcriptional regulators such as ARID1A, TP53, FOXO1, BRD1, NCOA1, and NCOR2 with one local recurrence gaining three PIK3CA mutations. In contrast to DNA-level changes, we discovered recurrent outlier mRNA expression alterations were common-including outlier gains in TP63 (n = 5 cases [42%]), NTRK3 (n = 5 [42%]), NTRK2 (n = 4 [33%]), PAX3 (n = 4 [33%]), FGFR4 (n = 3 [25%]), and TERT (n = 3 [25%]). Recurrent losses involved ESR1 (n = 5 [42%]), RELN (n = 5 [42%]), SFRP4 (n = 4 [33%]), and FOSB (n = 4 [33%]). ESR1-depleted recurrences harbored shared transcriptional remodeling events including upregulation of PROM1 and other basal cancer markers. CONCLUSIONS: Taken together, this study defines acquired genomic changes in long-term, estrogen-deprived disease; highlights the importance of longitudinal RNA profiling; and identifies a common ESR1-depleted endocrine-resistant breast cancer subtype with basal-like transcriptional reprogramming.

An immune-humanized patient-derived xenograft model of estrogen-independent, hormone receptor positive metastatic breast cancer.

BACKGROUND: Metastatic breast cancer (MBC) is incurable, with a 5-year survival rate of 28%. In the USA, more than 42,000 patients die from MBC every year. The most common type of breast cancer is estrogen receptor-positive (ER+), and more patients die from ER+ breast cancer than from any other subtype. ER+ tumors can be successfully treated with hormone therapy, but many tumors acquire endocrine resistance, at which point treatment options are limited. There is an urgent need for model systems that better represent human ER+ MBC in vivo, where tumors can metastasize. Patient-derived xenografts (PDX) made from MBC spontaneously metastasize, but the immunodeficient host is a caveat, given the known role of the immune system in tumor progression and response to therapy. Thus, we attempted to develop an immune-humanized PDX model of ER+ MBC. METHODS: NSG-SGM3 mice were immune-humanized with CD34+ hematopoietic stem cells, followed by engraftment of human ER+ endocrine resistant MBC tumor fragments. Strategies for exogenous estrogen supplementation were compared, and immune-humanization in blood, bone marrow, spleen, and tumors was assessed by flow cytometry and tissue immunostaining. Characterization of the new model includes assessment of the human tumor microenvironment performed by immunostaining. RESULTS: We describe the development of an immune-humanized PDX model of estrogen-independent endocrine resistant ER+ MBC. Importantly, our model harbors a naturally occurring ESR1 mutation, and immune-humanization recapitulates the lymphocyte-excluded and myeloid-rich tumor microenvironment of human ER+ breast tumors. CONCLUSION: This model sets the stage for development of other clinically relevant models of human breast cancer and should allow future studies on mechanisms of endocrine resistance and tumor-immune interactions in an immune-humanized in vivo setting.

Patient treatment and outcome after breast cancer orbital and periorbital metastases: a comprehensive case series including analysis of lobular versus ductal tumor histology.

BACKGROUND: Breast cancer is the most common malignancy to spread to the orbit and periorbit, and the invasive lobular carcinoma (ILC) histologic subtype of breast cancer has been reported to form these ophthalmic metastases (OM) more frequently than invasive ductal carcinomas (IDC). We herein report our single academic institution experience with breast cancer OM with respect to anatomical presentation, histology (lobular vs. ductal), treatment, and survival. METHODS: We employed the natural language processing platform, TIES (Text Information Extraction System), to search 2.3 million de-identified patient pathology and radiology records at our institution in order to identify patients with OM secondary to breast cancer. We then compared the resultant cohort, the "OM cohort," to two other representative metastatic breast cancer patient (MBC) databases from our institution. Histological analysis of selected patients was performed. RESULTS: Our TIES search and manual refinement ultimately identified 28 patients who were diagnosed with breast cancer between 1995 and 2016 that subsequently developed OM. Median age at diagnosis was 54 (range 28-77) years of age. ER, PR, and HER2 status from the 28 patients with OM did not differ from other patients with MBC from our institution. The relative proportion of patients with ILC was significantly higher in the OM cohort (32.1%) than in other MBC patients in our institution (11.3%, p = 0.007). Median time to first OM in the OM cohort was 46.7 months, and OM were the second most frequent first metastases after bony metastases. After diagnosis of the first distant metastasis of any kind, median survival of patients with ILC (21.4 months) was significantly shorter than that of patients with IDC (55.3 months, p = 0.03). Nine patients developed bilateral OM. We observed a significant co-occurrence of OM and central nervous system metastases (p = 0.0053). The histological analysis revealed an interesting case in which the primary tumor was of a mixed ILC/IDC subtype, while only ILC was present in the OM. CONCLUSIONS: OM from breast cancer are illustrative of the difference in metastatic behavior of ILC versus IDC and should be considered when treating patients with ILC, especially in those with complaints of visual acuity changes.

Targetable ERBB2 mutation status is an independent marker of adverse prognosis in estrogen receptor positive, ERBB2 non-amplified primary lobular breast carcinoma: a retrospective in silico analysis of public datasets.

BACKGROUND: Invasive lobular carcinoma (ILC) accounts for 10-15% of primary breast cancers and is typically estrogen receptor alpha positive (ER+) and ERBB2 non-amplified. Somatic mutations in ERBB2/3 are emerging as a tractable mechanism underlying enhanced human epidermal growth factor 2 (HER2) activity. We tested the hypothesis that therapeutically targetable ERBB2/3 mutations in primary ILC of the breast associate with poor survival outcome in large public datasets. METHODS: We performed in silico comparison of ERBB2 non-amplified cases of ER+ stage I-III primary ILC (N = 279) and invasive ductal carcinoma (IDC, N = 1301) using METABRIC, TCGA, and MSK-IMPACT information. Activating mutations amenable to HER2-directed therapy with neratinib were identified using existing functional data from in vitro cell line and xenograft experiments. Multivariate analysis of 10-year overall survival (OS) with tumor size, grade, and lymph node status was performed using a Cox regression model. Differential gene expression analyses by ERBB2 mutation and amplification status was performed using weighted average differences and an in silico model of response to neratinib derived from breast cancer cell lines. RESULTS: ILC tumors comprised 17.7% of all cases in the dataset but accounted for 47.1% of ERBB2-mutated cases. Mutations in ERBB2 were enriched in ILC vs. IDC cases (5.7%, N = 16 vs. 1.4%, N = 18, p < 0.0001) and clustered in the tyrosine kinase domain of HER2. ERBB3 mutations were not enriched in ILC (1.1%, N = 3 vs. 1.8%, N = 23; p = 0.604). Median OS for patients with ERBB2-mutant ILC tumors was 66 months vs. 211 months for ERBB2 wild-type (p = 0.0001), and 159 vs. 166 months (p = 0.733) for IDC tumors. Targetable ERBB2 mutational status was an independent prognostic marker of 10-year OS-but only in ILC (hazard ratio, HR = 3.7, 95% CI 1.2-11.0; p = 0.021). Findings were validated using a novel ERBB2 mutation gene enrichment score (HR for 10-year OS in ILC = 2.3, 95% CI 1.04-5.05; p = 0.040). CONCLUSIONS: Targetable ERBB2 mutations are enriched in primary ILC and their detection represents an actionable strategy with the potential to improve patient outcomes. Biomarker-led clinical trials of adjuvant HER-targeted therapy are warranted for patients with ERBB2-mutated primary ILC.

ADAM22/LGI1 complex as a new actionable target for breast cancer brain metastasis.

BACKGROUND: Metastatic breast cancer is a major cause of cancer-related deaths in woman. Brain metastasis is a common and devastating site of relapse for several breast cancer molecular subtypes, including oestrogen receptor-positive disease, with life expectancy of less than a year. While efforts have been devoted to developing therapeutics for extra-cranial metastasis, drug penetration of blood-brain barrier (BBB) remains a major clinical challenge. Defining molecular alterations in breast cancer brain metastasis enables the identification of novel actionable targets. METHODS: Global transcriptomic analysis of matched primary and metastatic patient tumours (n = 35 patients, 70 tumour samples) identified a putative new actionable target for advanced breast cancer which was further validated in vivo and in breast cancer patient tumour tissue (n = 843 patients). A peptide mimetic of the target's natural ligand was designed in silico and its efficacy assessed in in vitro, ex vivo and in vivo models of breast cancer metastasis. RESULTS: Bioinformatic analysis of over-represented pathways in metastatic breast cancer identified ADAM22 as a top ranked member of the ECM-related druggable genome specific to brain metastases. ADAM22 was validated as an actionable target in in vitro, ex vivo and in patient tumour tissue (n = 843 patients). A peptide mimetic of the ADAM22 ligand LGI1, LGI1MIM, was designed in silico. The efficacy of LGI1MIM and its ability to penetrate the BBB were assessed in vitro, ex vivo and in brain metastasis BBB 3D biometric biohybrid models, respectively. Treatment with LGI1MIM in vivo inhibited disease progression, in particular the development of brain metastasis. CONCLUSION: ADAM22 expression in advanced breast cancer supports development of breast cancer brain metastasis. Targeting ADAM22 with a peptide mimetic LGI1MIM represents a new therapeutic option to treat metastatic brain disease.

MetaOmics: analysis pipeline and browser-based software suite for transcriptomic meta-analysis.

SUMMARY: The rapid advances of omics technologies have generated abundant genomic data in public repositories and effective analytical approaches are critical to fully decipher biological knowledge inside these data. Meta-analysis combines multiple studies of a related hypothesis to improve statistical power, accuracy and reproducibility beyond individual study analysis. To date, many transcriptomic meta-analysis methods have been developed, yet few thoughtful guidelines exist. Here, we introduce a comprehensive analytical pipeline and browser-based software suite, called MetaOmics, to meta-analyze multiple transcriptomic studies for various biological purposes, including quality control, differential expression analysis, pathway enrichment analysis, differential co-expression network analysis, prediction, clustering and dimension reduction. The pipeline includes many public as well as >10 in-house transcriptomic meta-analytic methods with data-driven and biological-aim-driven strategies, hands-on protocols, an intuitive user interface and step-by-step instructions. AVAILABILITY AND IMPLEMENTATION: MetaOmics is freely available at https://github.com/metaOmics/metaOmics. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Network-guided prediction of aromatase inhibitor response in breast cancer.

Prediction of response to specific cancer treatments is complicated by significant heterogeneity between tumors in terms of mutational profiles, gene expression, and clinical measures. Here we focus on the response of Estrogen Receptor (ER)+ post-menopausal breast cancer tumors to aromatase inhibitors (AI). We use a network smoothing algorithm to learn novel features that integrate several types of high throughput data and new cell line experiments. These features greatly improve the ability to predict response to AI when compared to prior methods. For a subset of the patients, for which we obtained more detailed clinical information, we can further predict response to a specific AI drug.

Targeted mutation detection in breast cancer using MammaSeq™.

BACKGROUND: Breast cancer is the most common invasive cancer among women worldwide. Next-generation sequencing (NGS) has revolutionized the study of cancer across research labs around the globe; however, genomic testing in clinical settings remains limited. Advances in sequencing reliability, pipeline analysis, accumulation of relevant data, and the reduction of costs are rapidly increasing the feasibility of NGS-based clinical decision making. METHODS: We report the development of MammaSeq, a breast cancer-specific NGS panel, targeting 79 genes and 1369 mutations, optimized for use in primary and metastatic breast cancer. To validate the panel, 46 solid tumors and 14 plasma circulating tumor DNA (ctDNA) samples were sequenced to a mean depth of 2311× and 1820×, respectively. Variants were called using Ion Torrent Suite 4.0 and annotated with cravat CHASM. CNVKit was used to call copy number variants in the solid tumor cohort. The oncoKB Precision Oncology Database was used to identify clinically actionable variants. Droplet digital PCR was used to validate select ctDNA mutations. RESULTS: In cohorts of 46 solid tumors and 14 ctDNA samples from patients with advanced breast cancer, we identified 592 and 43 protein-coding mutations. Mutations per sample in the solid tumor cohort ranged from 1 to 128 (median 3), and the ctDNA cohort ranged from 0 to 26 (median 2.5). Copy number analysis in the solid tumor cohort identified 46 amplifications and 35 deletions. We identified 26 clinically actionable variants (levels 1-3) annotated by OncoKB, distributed across 20 out of 46 cases (40%), in the solid tumor cohort. Allele frequencies of ESR1 and FOXA1 mutations correlated with CA.27.29 levels in patient-matched blood draws. CONCLUSIONS: In solid tumor biopsies and ctDNA, MammaSeq detects clinically actionable mutations (OncoKB levels 1-3) in 22/46 (48%) solid tumors and in 4/14 (29%) of ctDNA samples. MammaSeq is a targeted panel suitable for clinically actionable mutation detection in breast cancer.

SNAIL is induced by tamoxifen and leads to growth inhibition in invasive lobular breast carcinoma.

PURPOSE: Invasive lobular carcinoma (ILC) is a histological subtype of breast cancer that is predominantly estrogen receptor alpha (ER)-positive (+) and is thus treated with endocrine therapies. Herein, we sought to understand the molecular underpinnings of the 4-hydroxytamoxifen (4OHT) resistance in ILC by assessing the potential role of the epithelial-to-mesenchymal transition transcription factor (EMT-TF) SNAIL (SNAI1). METHODS: Using a series of breast cancer cell lines, we measured the basal, estrogen and 4OHT-induced expression of SNAIL and other EMT-TF family members by quantitative reverse transcription-polymerase chain reaction and immunoblotting. Chromatin immunoprecipitation experiments were performed to assess ER binding to the SNAIL promoter. Cell proliferation, cell cycle and apoptosis were assessed in 2D cultures. 3D growth was assessed in Matrigel and Collagen I cultures. RESULTS: Estrogen and 4OHT induced SNAIL expression, but not that of the other EMT-TF family members SLUG (SNAI2) and SMUC (SNAI3), with the 4OHT effect being specific to the lobular but not the ductal subtype. We observed estrogen and 4OHT-induced ER recruitment to the SNAI1 promoter and high endogenous basal levels of SNAIL and several EMT-TFs in ILC cell lines. While SNAIL knockdown had a minor impact on the 4OHT partial agonism in estrogen-depleted conditions, it led to a surprising increase in cell proliferation in full serum. In complementary experiments, inducible SNAI1 overexpression caused decreased proliferation, associated with a cell cycle arrest in G0/G1. Additionally, apoptosis was observed in BCK4 cells. CONCLUSION: These data suggest a previously unrecognized role for SNAIL in ILC, substantiating a context-dependent behavior for this EMT-TF.