Search

Chromatin accessibility states that influence gene expression and other nuclear processes can be altered in disease. The constellation of transcription factors and chromatin regulatory complexes in cells results in characteristic patterns of chromatin accessibility. The study of these patterns in tissues has been limited because existing chromatin accessibility assays are ineffective for archival formalin-fixed, paraffin-embedded (FFPE) tissues. We have developed a method to efficiently extract intact chromatin from archival tissue via enhanced cavitation with a nanodroplet reagent consisting of a lipid shell with a liquid perfluorocarbon core. Inclusion of nanodroplets during the extraction of chromatin from FFPE tissues enhances the recovery of intact accessible and nucleosome-bound chromatin. We show that the addition of nanodroplets to the chromatin accessibility assay formaldehyde-assisted isolation of regulatory elements (FAIRE), does not affect the accessible chromatin signal. Applying the technique to FFPE human tumor xenografts, we identified tumor-relevant regions of accessible chromatin shared with those identified in primary tumors. Further, we deconvoluted non-tumor signal to identify cellular components of the tumor microenvironment. Incorporation of this method of enhanced cavitation into FAIRE offers the potential for extending chromatin accessibility to clinical diagnosis and personalized medicine, while also enabling the exploration of gene regulatory mechanisms in archival samples.

Inhibition of the HER2/ERBB2 receptor is a keystone to treating HER2-positive malignancies, particularly breast cancer, but a significant fraction of HER2-positive (HER2+) breast cancers recur or fail to respond. Anti-HER2 monoclonal antibodies, like trastuzumab or pertuzumab, and ATP active site inhibitors like lapatinib, commonly lack durability because of adaptive changes in the tumor leading to resistance. HER2+ cell line responses to inhibition with lapatinib were analyzed by RNAseq and ChIPseq to characterize transcriptional and epigenetic changes. Motif analysis of lapatinib-responsive genomic regions implicated the pioneer transcription factor FOXA1 as a mediator of adaptive responses. Lapatinib in combination with FOXA1 depletion led to dysregulation of enhancers, impaired adaptive upregulation of HER3, and decreased proliferation. HER2-directed therapy using clinically relevant drugs (trastuzumab with or without lapatinib or pertuzumab) in a 7-day clinical trial designed to examine early pharmacodynamic response to antibody-based anti-HER2 therapy showed reduced FOXA1 expression was coincident with decreased HER2 and HER3 levels, decreased proliferation gene signatures, and increased immune gene signatures. This highlights the importance of the immune response to anti-HER2 antibodies and suggests that inhibiting FOXA1-mediated adaptive responses in combination with HER2 targeting is a potential therapeutic strategy.

Background Human papillomavirus-positive (HPV+) squamous cell carcinoma of the oropharynx (OPSCC) is the most prevalent HPV-associated malignancy in the United States. Favorable treatment outcomes have led to increased interest in treatment de-escalation to reduce treatment morbidity as well as the development of prognostic markers to identify appropriately low-risk patients. Intratumoral genomic heterogeneity and copy number alteration burden have been demonstrated to be predictive of poor outcomes in many other cancers; therefore, we sought to determine whether intratumor heterogeneity and genomic instability are associated with poor outcomes in HPV+ OPSCC. Methods Tumor heterogeneity estimates were made based on targeted exome sequencing of 45 patients with HPV+ OPSCC tumors. Analysis of an additional cohort of HPV+ OPSCC tumors lacking matched normal sequencing allowed copy number analysis of 99 patient tumors. Results High intratumorally genomic heterogeneity and high numbers of copy number alterations were strongly associated with worse recurrence-free survival. Tumors with higher heterogeneity and frequent copy number alterations were associated with loss of distal 11q, which encodes key genes related to double-strand break repair, including ATM and MRE11A. Conclusions Both intratumor genomic heterogeneity and high-burden copy number alterations are strongly associated with poor recurrence-free survival in patients with HPV+ OPSCC. The drivers of genomic instability and heterogeneity in these tumors remains to be elucidated. However, 11q loss and defective DNA double-strand break repair have been associated with genomic instability in other solid tumors. Copy number alteration burden and intratumoral heterogeneity represent promising avenues for risk stratification of patients with HPV+OPSCC.

Background ONC201 is a dopamine receptor D2 (DRD2) antagonist that inhibits tumor growth in preclinical models through ClpP activation to induce integrated stress response pathway and mitochondrial events related to inhibition of cell growth, which is being explored in clinical trials for solid tumors and hematological malignancies. In this study, we investigated the anti-tumorigenic effect of ONC201 in endometrial cancer cell lines and a genetically engineered mouse model of endometrial cancer. Methods Cell proliferation was assessed by MTT and colony formation assays. Cell cycle and apoptosis were evaluated by Cellometer. Invasion capacity was tested using adhesion, transwell and wound healing assays. LKB1fl/flp53fl/fl mouse model of endometrial cancer were fed a control low fat diet versus a high fat diet to mimic diet-induced obesity. Following tumor onset, mice were treated with placebo or ONC201. Metabolomics and lipidomics were used to identify the obesity-dependent effects of ONC201 in the mouse endometrial tumors. DRD2 expression was analyzed by immunohistochemistry in human endometrioid and serous carcinoma specimens. DRD2 mRNA expression from the Cancer Genome Atlas (TCGA) database was compared between the four molecular subtypes of endometrial cancer. Results Increasing DRD2 expression in endometrial cancer was significantly associated with grade, serous histology and stage, as well as worse progression free survival and overall survival. Higher expression of DRD2 mRNA was found for the Copy Number High (CNH) subtype when compared to the other subtypes. ONC201 inhibited cell proliferation, induced cell cycle G1 arrest, caused cellular stress and apoptosis and reduced invasion in endometrial cancer cells. Diet-induced obesity promoted endometrial tumor growth while ONC201 exhibited anti-tumorigenic efficacy in the obese and lean LKB1fl/fl/p53fl/fl mice. Metabolomic analysis demonstrated that ONC201 reversed the obesity-driven upregulation of lipid biosynthesis and reduced protein biosynthesis in obese and lean mice. Conclusion ONC201 has anti-tumorigenic effects in endometrial cancer cells and a transgenic mouse model of endometrial cancer, and DRD2 expression was documented in both human serous and endometrioid endometrial cancer. These studies support DRD2 antagonism via ONC201 as a promising therapeutic strategy for endometrial cancer that has already demonstrated pharmacodynamic activity and clinical benefit in both serous and endometrioid endometrial cancer patients.

PURPOSE CALGB 40601 assessed whether dual versus single human epidermal growth factor receptor 2 (HER2) –targeting drugs added to neoadjuvant chemotherapy increased pathologic complete response (pCR). Here, we report relapse-free survival (RFS), overall survival (OS), and gene expression signatures that predict pCR and survival. PATIENTS AND METHODS Three hundred five women with untreated stage II and III HER2-positive breast cancer were randomly assigned to receive weekly paclitaxel combined with trastuzumab plus lapatinib (THL), trastuzumab (TH), or lapatinib (TL). The primary end point was pCR, and secondary end points included RFS, OS, and gene expression analyses. mRNA sequencing was performed on 264 pretreatment samples. RESULTS One hundred eighteen patients were randomly allocated to THL, 120 to TH, and 67 to TL. At more than 7 years of follow-up, THL had significantly better RFS and OS than did TH (RFS hazard ratio, 0.32; 95% CI, 0.14 to 0.71; P = .005; OS hazard ratio, 0.34; 95% CI, 0.12 to 0.94; P = .037), with no difference between TH and TL. Of 688 previously described gene expression signatures, significant associations were found in 215 with pCR, 45 with RFS, and only 22 with both pCR and RFS (3.2%). Specifically, eight immune signatures were significantly correlated with a higher pCR rate and better RFS. Among patients with residual disease, the immunoglobulin G signature was an independent, good prognostic factor, whereas the HER2-enriched signature, which was associated with a higher pCR rate, showed a significantly shorter RFS. CONCLUSION In CALGB 40601, dual HER2-targeting resulted in significant RFS and OS benefits. Integration of intrinsic subtype and immune signatures allowed for the prediction of pCR and RFS, both overall and within the residual disease group. These approaches may provide means for rational escalation and de-escalation treatment strategies in HER2-positive breast cancer.

Purpose Claudin-low molecular subtypes have been identified in breast and bladder cancers and are characterized by low expression of claudins, enrichment for epithelial-to-mesenchymal transition (EMT), and tumor-initiating cell (TIC) features. We evaluated whether the claudin-low subtype also exists in gastric cancer. Materials and Methods Four hundred fifteen tumors from The Cancer Genome Atlas (TCGA) gastric cancer mRNA data set were clustered on the claudin, EMT, and TIC gene sets to identify claudin-low tumors. We derived a 24-gene predictor that classifies gastric cancer into claudin-low and non–claudin-low subtypes. This predictor was validated with the Asian Cancer Research Group (ACRG) data set. We characterized molecular and clinical features of claudin-low tumors. Results We identified 46 tumors that had consensus enrichment for claudin-low features in TCGA data set. Claudin-low tumors were most commonly diffuse histologic type (82%) and originally classified as TCGA genomically stable (GS) subtype (78%). Compared with GS subtype, claudin-low subtype had significant activation in Rho family of GTPases signaling, which appears to play a key role in its EMT and TIC properties. In the ACRG data set, 28 of 300 samples were classified as claudin-low tumors by the 24-gene predictor and were phenotypically similar to the initially derived claudin-low tumors. Clinically, claudin-low subtype had the worst overall survival. Of note, the hazard ratios that compared claudin-low versus GS subtype were 2.10 (95% CI, 1.07 to 4.11) in TCGA and 2.32 (95% CI, 1.18 to 4.55) in the ACRG cohorts, with adjustment for age and pathologic stage. Conclusion We identified a gastric claudin-low subtype that carries a poor prognosis likely related to therapeutic resistance as a result of its EMT and TIC phenotypes.

Mechanisms driving tumor progression from less aggressive subtypes to more aggressive states represent key targets for therapy. We identified a subset of luminal A primary breast tumors that give rise to HER2-enriched (HER2E) subtype metastases, but remain clinically HER2 negative (cHER2(–)). By testing the unique genetic and transcriptomic features of these cases, we developed the hypothesis that FGFR4 likely participates in this subtype switching. To evaluate this, we developed 2 FGFR4 genomic signatures using a patient-derived xenograft (PDX) model treated with an FGFR4 inhibitor, which inhibited PDX growth in vivo. Bulk tumor gene expression analysis and single-cell RNA sequencing demonstrated that the inhibition of FGFR4 signaling caused molecular switching. In the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) breast cancer cohort, FGFR4-induced and FGFR4-repressed signatures each predicted overall survival. Additionally, the FGFR4-induced signature was an independent prognostic factor beyond subtype and stage. Supervised analysis of 77 primary tumors with paired metastases revealed that the FGFR4-induced signature was significantly higher in luminal/ER(+) tumor metastases compared with their primaries. Finally, multivariate analysis demonstrated that the FGFR4-induced signature also predicted site-specific metastasis for lung, liver, and brain, but not for bone or lymph nodes. These data identify a link between FGFR4-regulated genes and metastasis, suggesting treatment options for FGFR4-positive patients, whose high expression is not caused by mutation or amplification.

We hypothesized that human genes differ by their sensitivity to ultraviolet (UV) exposure. We used somatic mutations detected by genome-wide screens in melanoma and reported in the Catalog Of Somatic Mutations In Cancer. As a measure of UV sensitivity, we used the number of silent mutations generated by C>T transitions in pyrimidine dimers of a given transcript divided by the number of potential sites for this type of mutations in the transcript. We found that human genes varied by UV sensitivity by two orders of magnitude. We noted that the melanoma-associated tumor suppressor gene CDKN2A was among the top five most UV-sensitive genes in the human genome. Melanoma driver genes have a higher UV-sensitivity compared with other genes in the human genome. The difference was more prominent for tumor suppressors compared with oncogene. The results of this study suggest that differential sensitivity of human transcripts to UV light may explain melanoma specificity of some driver genes. Practical significance of the study relates to the fact that differences in UV sensitivity among human genes need to be taken into consideration whereas predicting melanoma-associated genes by the number of somatic mutations detected in a given gene.

BackgroundPan-cancer genomic analyses based on the magnitude of pathway activity are currently lacking. Focusing on the cell cycle, we examined the DNA mutations and chromosome arm-level aneuploidy within tumours with low, intermediate and high cell cycle activity.Patients and methodsMatching mRNA, gene mutational status, chromosomal arm-level aberrations and clinico-pathological data was assembled from pan-cancer studies of 9,515 patients with 32 different cancers. Cell cycle activity was estimated from mRNA data using the cell cycle score (CCS) signature. Barplots were used to visualise mutation and chromosomal aberration frequency within CCS subgroups. Kaplan-Meier and multivariable Cox-regression analyses were used to determine survival differences between CCS subgroups.ResultsCell cycle activity varied broadly across and within all cancers. TP53 and PIK3CA mutations were common in all CCS subgroups but with increasing frequency as cell cycle activity levels increased (P < 0.001). Mutations in BRAF and gains in 16p were less frequent CCS high tumours (P < 0.001). In Kaplan-Meier analysis, patients whose tumours were CCS Low had a longer PFI relative to intermediate or high (P < 0.001) and this significance remained in multivariable analysis (CCS intermediate: HR = 1.37; 95% CI 1.17 – 1.60, CCS high: 1.54; 1.29 – 1.84, CCS Low = Ref).ConclusionsCell cycle activity varies across and within cancers and whilst similar DNA alterations can be found at all activity levels, some notable exceptions exist. These data also demonstrate that independent prognostic information can be derived on a pan-cancer level from a simple measure of cell cycle activity.

Background: It has become increasingly clear that effective treatment of metastatic breast cancer (MBC) requires an in-depth understanding of the molecular differences between primary tumors and metastases. The AURORA US Network was established to collect primary breast cancer-metastasis pairs for multi-platform genomic profiling in order to identify the molecular drivers of metastatic disease. AURORA US has both a retrospective and prospective phase. This is the first report of the retrospective phase. Methods: Archived tissue samples from the primary tumor and at least one distant metastasis were retrospectively collected from 83 MBC patients. Following internal quality assessment, samples from 55 pts, including 105 distinct metastatic lesions, were subject to DNA low pass whole genome and exome sequencing, DNA methylation arrays, and RNA sequencing. Early analyses of these multi-platform data include: DNA methylation, tumor gene expression and microenvironmental signatures, somatic and germline variants, DNA copy number changes, and structural variants between breast primaries and matched metastases. Results: Median age at diagnosis was 49 years (25-76); 32 (58%) were Stage I or II at presentation, 27 (49%) had a family history of breast cancer, and 20 (36%) had a second breast primary. Median disease-free interval before developing MBC was 2 years (range 0-36, 5 patients presented with Stage IV). Median overall survival from initial presentation was 4 years (range 0-37). Median survival after developing MBC was 1 year (range 0-13), with a median of three treatments. Primary tissue samples were banked from 1977-2017 and metastases were banked from 1999-2017. Clinical phenotypes of the primaries included 27 HR+ (49%), 15 triple negative (TNBC, 27%), and 11 HER2+ (20%, 12 missing HER2 status). Intrinsic subtype distribution of the primaries included 17 Basal-like (31%), 17 Luminal A (31%), 7 Normal-like (13%), 5 HER2-enriched (9%), and 1 Luminal B, with 8 pending. All metastases from the Basal-like cases remained Basal-like, while metastases from luminal primaries tended to gain HER2-Enriched subtype features (5/18, p = 0.01). Overall, we identified significant metastasis-enriched alterations in metabolism pathways, an increase in proliferation, and the loss of differentiation signatures and immune infiltrates with progression; the latter being the most pronounced in brain metastases. The most frequent somatic mutations in this cohort were in TP53, NCOR1, and RUNX1. Interestingly, ERBB2, EGFR, and ATM were also mutated in ≥10% of the tumors sequenced. In almost all cases, CpG island hypermethylation was clonally present in the primary tumor and persisted stably in the majority of metastatic lesions. Promoter CpG island hypermethylation was also identified in some metastatic lesions at JAM3, an important cellular adhesion molecule,and this was accompanied by reduced mRNA expression. CONCLUSIONS: Collection of banked primary and metastatic tissue pairs identified a young MBC cohort with a high frequency of breast cancer family history and second breast primaries. Molecular characterization of luminal tumor pairs highlighted acquisition of aggressive traits including increased proliferation and loss of differentiation in the metastases. In contrast, basal-like pairs remained relatively unchanged, except for the loss of immune activation. Ongoing analyses to be presented include clonal heterogeneity and phylogeny, novel metastasis signature discovery, gene fusion, and endogenous retrovirus detection. Citation Format: Tari A King, Minetta C Liu, Marni B McClure, Toshinori Hinoue, Benjamin J Kelly, Chad J Creighton, Jay Bowen, Kristen Leraas, Robyn T Burns, Sara Coppens, Salma Rezk, Amy L Garrett, Justin M Balko, Joel S Parker, Ben H Park, Ian Krop, Carey Anders, Katherine A Hoadley, Julie Gastier-Foster, Mothaffar F Rimawi, Rita Nanda, Nancy U Lin, Claudine Isaacs, P. Kelly Marcom, Anna Maria Storniolo, Fergus J Couch, Elaine R Mardis, Adrian V Lee, Uma Chandran, Peter W Laird, Susan G Hilsenbeck, Larry Norton, Andrea L Richardson, W. Fraser Symmans, Lisa A Carey, Antonio C Wolff, Nancy E Davidson, Charles M Perou, the AURORA US Network. Multiplatform analysis of matched primary and metastatic breast tumors from the AURORA US Network [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr GS3-08.

Aberrant activation of the androgen receptor (AR) may play a critical role in castration resistant prostate cancer. After ligand binding, AR is recruited to the androgen responsive element (ARE) sequences on the DNA where AR interaction with coactivators and corepressors modulates transcription. We demonstrated that phosphorylation of AR at Tyr-267 by Ack1/TNK2 tyrosine kinase results in nuclear translocation, DNA binding, and androgen-dependent gene transcription in a low androgen environment. In order to dissect downstream mechanisms, we searched for proteins whose interaction with AR was regulated by Ack1. SLIRP (SRA stem-loop interacting RNA binding protein) was identified as a candidate protein. Interaction between AR and SLIRP was disrupted by Ack1 kinase activity as well as androgen or heregulin treatment. The noncoding RNA, SRA, was required for AR-SLIRP interaction. SLIRP was bound to ARE’s of AR target genes in the absence of androgen. Treatment with androgen or heregulin led to dissociation of SLIRP from the ARE. Whole transcriptome analysis of SLIRP knockdown in androgen responsive LNCaP cells showed that SLIRP affects a significant subset of androgen-regulated genes. Our data suggest that Ack1 kinase and androgen regulate interaction between AR and SLIRP and that SLIRP functions as a coregulator of AR with properties of a corepressor in a context-dependent manner.

Cutaneous melanoma can be lethal even if detected at an early stage. Epigenetic profiling may facilitate the identification of aggressive primary melanomas with unfavorable outcomes. We performed clustering of whole-genome methylation data to identify subclasses that were then assessed for survival, clinical features, methylation patterns, and biological pathways. Among 89 cutaneous primary invasive melanomas, we identified three methylation subclasses exhibiting low methylation, intermediate methylation, or hypermethylation of CpG islands, known as the CpG island methylator phenotype (CIMP). CIMP melanomas occurred as early as tumor stage 1b and, compared with low-methylation melanomas, were associated with age at diagnosis ≥65 years, lentigo maligna melanoma histologic subtype, presence of ulceration, higher American Joint Committee on Cancer stage and tumor stage, and lower tumor-infiltrating lymphocyte grade (all P0.05). Patients with CIMP melanomas had worse melanoma-specific survival (hazard ratio = 11.84; confidence interval = 4.65‒30.20) than those with low-methylation melanomas, adjusted for age, sex, American Joint Committee on Cancer stage, and tumor-infiltrating lymphocyte grade. Genes hypermethylated in CIMP compared with those in low-methylation melanomas included PTEN, VDR, PD-L1, TET2, and gene sets related to development/differentiation, the extracellular matrix, and immunity. CIMP melanomas exhibited hypermethylation of genes important in melanoma progression and tumor immunity, and although present in some early melanomas, CIMP was associated with worse survival independent of known prognostic factors.

Interactions between WD40 repeat domain protein 5 (WDR5) and its various partners such as mixed lineage leukemia (MLL) and c-MYC are essential for sustaining oncogenesis in human cancers. However, inhibitors that block protein-protein interactions (PPIs) between WDR5 and its binding partners exhibit modest cancer cell killing effects and lack in vivo efficacy. Here, we present pharmacological degradation of WDR5 as a promising therapeutic strategy for treating WDR5-dependent tumors and report two high-resolution crystal structures of WDR5-degrader-E3 ligase ternary complexes. We identified an effective WDR5 degrader via structure-based design and demonstrated its in vitro and in vivo antitumor activities. On the basis of the crystal structure of an initial WDR5 degrader in complex with WDR5 and the E3 ligase von Hippel–Lindau (VHL), we designed a WDR5 degrader, MS67, and demonstrated the high cooperativity of MS67 binding to WDR5 and VHL by another ternary complex structure and biophysical characterization. MS67 potently and selectively depleted WDR5 and was more effective than WDR5 PPI inhibitors in suppressing transcription of WDR5-regulated genes, decreasing the chromatin-bound fraction of MLL complex components and c-MYC, and inhibiting the proliferation of cancer cells. In addition, MS67 suppressed malignant growth of MLL-rearranged acute myeloid leukemia patient cells in vitro and in vivo and was well tolerated in vivo. Collectively, our results demonstrate that structure-based design can be an effective strategy to identify highly active degraders and suggest that pharmacological degradation of WDR5 might be a promising treatment for WDR5-dependent cancers.

Chromosome 11q13-14 amplification is a defining feature of high-risk hormone receptor-positive (HR+) breast cancer; however, the mechanism(s) by which this amplicon contributes to breast tumorigenesis remains unclear. In the current study, proteogenomic analyses of >3000 breast tumors from the TCGA, METABRIC and CPTAC studies demonstrated that carnitine palmitoyltransferase 1A (CPT1A), which is localized to this amplicon, is overexpressed at the mRNA and protein level in aggressive luminal tumors, strongly associated with indicators of tumor proliferation and a predictor of poor prognosis. In vitro genetic studies demonstrated that CPT1A is required for and can promote luminal breast cancer proliferation, survival, as well as colony and mammosphere formation. Since CPT1A is the rate-limiting enzyme during fatty acid oxidation (FAO), our data indicate that FAO may be essential for these tumors. Pharmacologic inhibition of FAO prevented in vitro and in vivo tumor growth and cell proliferation as well as promoted apoptosis in luminal breast cancer cells and orthotopic xenograft tumor models. Collectively, our data establish an oncogenic role for CPT1A and FAO in HR+ luminal tumors and provide preclinical evidence and rationale supporting further investigation of FAO as a potential therapeutic opportunity for the treatment of HR+ breast cancer.

Lymph node (LN) metastases correspond with a worse prognosis in nearly all cancers, yet the occurrence of cancer spreading from LNs remains controversial. Additionally, the mechanisms explaining how cancers survive and exit LNs are largely unknown. Here, we show that breast cancer patients frequently have LN metastases that closely resemble distant metastases. In addition, using a microsurgical model, we show how LN metastasis development and dissemination is regulated by the expression of a chromatin modifier, histone deacetylase 11 (HDAC11). Genetic and pharmacologic blockade of HDAC11 decreases LN tumor growth, yet substantially increases migration and distant metastasis formation. Collectively, we reveal a mechanism explaining how HDAC11 plasticity promotes breast cancer growth as well as dissemination from LNs and suggest caution with the use of HDAC inhibitors., The prognosis of cancer patients with lymph node (LN) metastasis is worse than those without. Here, the authors report that while histone deacetylase 11 (HDAC11) inhibition suppresses tumor growth within the LN, it also promotes cancer cell migration out of the LN to form distant metastasis, and therefore suggest caution with HDAC inhibitors.

Background Measures of the adaptive immune response have prognostic and predictive associations in melanoma and other cancer types. Specifically, intratumoral T cell density and function have considerable prognostic and predictive value in skin cutaneous melanoma (SKCM). Less is known about the significance of tumor-infiltrating B cells in SKCM. Our goal was to understand the prognostic and predictive value of B cell phenotypic subsets in SKCM using RNA sequencing. Methods We used our previously published algorithm, V’DJer, to assemble B cell receptor (BCR) repertoires and estimate diversity from short-read RNA sequencing (RNA-seq). We applied machine learning-based cellular phenotype classifiers to measure relative similarity of bulk tumor sample gene expression profiles and different B cell phenotypes. We assessed these aspects of B cell biology in 473 SKCM from the Cancer Genome Atlas Project (TCGA) as well as in RNA-seq data corresponding to tumor samples procured from patients who received CTLA-4 and PD-1 inhibitors for metastatic SKCM. Results We found that the BCR repertoire was associated with different clinical factors, such as tumor tissue site and sex. However, increased clonality of the BCR repertoire was favorably prognostic in SKCM and was prognostic even after first conditioning on various clinical factors. Mutation burden was not correlated with any BCR measurement, and no specific mutation had an altered BCR repertoire. Lack of an assembled BCR in pre-treatment tumor tissues was associated with a lack of anti-tumor response to a CTLA-4 inhibitor in metastatic SKCM. Conclusions These findings suggest an important prognostic and predictive role for B cell characteristics in SKCM. This has implications for melanoma immunobiology and potential development of immunogenomics features to predict survival and response to immunotherapy. Electronic supplementary material The online version of this article (10.1186/s13073-019-0647-5) contains supplementary material, which is available to authorized users.

In tumors, extravascular fibrin forms provisional scaffolds for endothelial cell (EC) growth and motility during angiogenesis. We report that fibrin-mediated angiogenesis was inhibited and tumor growth delayed following postnatal deletion of Tgfbr2 in the endothelium of Cdh5-Cre(ERT2) Tgfbr2(fl/fl) mice (Tgfbr2(iECKO) mice). ECs from Tgfbr2(iECKO) mice failed to upregulate the fibrinolysis inhibitor plasminogen activator inhibitor 1 (Serpine1, also known as PAI-1), due in part to uncoupled TGF-β–mediated suppression of miR-30c. Bypassing TGF-β signaling with vascular tropic nanoparticles that deliver miR-30c antagomiRs promoted PAI-1–dependent tumor growth and increased fibrin abundance, whereas miR-30c mimics inhibited tumor growth and promoted vascular-directed fibrinolysis in vivo. Using single-cell RNA-Seq and a NanoString miRNA array, we also found that subtypes of ECs in tumors showed spectrums of Serpine1 and miR-30c expression levels, suggesting functional diversity in ECs at the level of individual cells; indeed, fresh EC isolates from lung and mammary tumor models had differential abilities to degrade fibrin and launch new vessel sprouts, a finding that was linked to their inverse expression patterns of miR-30c and Serpine1 (i.e., miR-30c(hi) Serpine1(lo) ECs were poorly angiogenic and miR-30c(lo) Serpine1(hi) ECs were highly angiogenic). Thus, by balancing Serpine1 expression in ECs downstream of TGF-β, miR-30c functions as a tumor suppressor in the tumor microenvironment through its ability to promote fibrin degradation and inhibit blood vessel formation.

Immune checkpoint inhibitor (ICI) therapies have led to remarkable clinical responses in cancers such as melanoma and non-small cell lung cancer. In breast cancer, current immunotherapy trials have placed an emphasis on triple negative breast cancers (TNBC), where early results suggest response rates of 10-20%. Thus, it is critical to identify predictive biomarkers to enhance patient selection for immunotherapy. With this goal in mind, we simulated a clinical trial employing anti-PD1 and anti-CTLA therapies in immune-intact genetically engineered mouse models (GEMMs) of TNBC. Testing of ICI therapies on 8 different GEMMs demonstrated that each model was resistant. Whole exome sequencing showed that each model also harbored a low mutation burden. Given that mutation load is predictive of immunotherapy response in other cancer types, and that Apobec3B activity is associated with higher tumor mutation burden (TMB) in breast cancer, we created two different tumor lines with overexpression of murine Apobec3. In contrast to the parental lines, the Apobec3 overexpressing lines showed an elevated tumor mutation burden and new mutations were consistent with the Apobec mutation signature. These TNBC lines with new mutations resulting from Apobec3 activity were exquisitely sensitive to anti-PD1/anti-CTLA4 combination therapy; as assessed by reduction in tumor volume and extended overall survival. To identify features that predict response, we examined resistant and sensitive tumors at pretreatment, at 1 week of treatment, and at end stage by flow cytometry and mRNA-seq. Gene expression profiling identified multiple immune signatures as predictive of response to ICI therapy; specifically CD8+ T-effector memory cells, CD4+ T-cells, and activated B-Cells. Similarly, gene expression analysis showed that these cell types increased at 1 week of therapy in sensitive models but not in resistant models. Flow cytometry confirmed these predictions. Next, we used an antibody based approach to separately deplete CD4+ T-Cells, CD8+ T-cells, or B-cells in Apobec3 mutagenized murine tumors receiving aPD1/aCTLA4 combination therapy. In each case, depletion of these populations significantly reduced the therapeutic response. However, mice receiving combination immunotherapy and depleted for CD8+ T-cells still exhibited a significant extension in overall survival compared to non-treated controls. In contrast, the CD4+ T-cell depleted mice and B-cell depleted mice exhibited no ICI therapeutic benefit. Together, these data point to key immune biomarkers of response to anti-PD1/anti-CTLA4 therapy; we have further developed a genomic predictor of ICI response using our murine models and will test this on a human TNBC data set. Lastly, this GEMM system provides a rich RNA-seq resource, and new immune-activated models for TNBC, which uncovered a key role for B-cells and CD4+ T-cells in response to ICI therapies. Citation Format: Hollern DP, Xu N, Mott KR, He X, Carey-Ewend K, Marron DS, Ford J, Parker JS, Vincent BG, Serody JS, Perou CM. Apobec3 induced mutagenesis sensitizes murine models of triple negative breast cancer to immunotherapy by activating B-cells and CD4+ T-cells [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr GS1-05.

Background: The expression of ER, PR and HER2 in breast cancer is determined by routine pathology assessment to indicate systemic treatment. Here, we investigated the ability of a standardized mRNA-based assay to predict ER, PR and HER2 status and treatment response. Methods: ESR1, PGR and ERBB2 expression was analyzed from the standardized nCounter-based PAM50 assay in 1,544 FFPE breast tumors obtained from 13 independent studies (NeoEribulin, GEICAM 2012-09, TBCRC023/006, LPT109096, EGF117165, PAMELA, PerELISA, GEICAM 2003-11, VENTANA, IBIMA and HCB). All immunohistochemical and in situ hybridization analyses followed ASCO/CAP criteria and were performed in central labs except for Neoeribulin and VENTANA studies. To explore the best cutoff for each gene, we used Monte Carlo cross validation (repeated 1000 times, 2/3 training, 1/3 testing) to achieve the highest kappa values when gene expression was compared to pathology assessment. Receiver operating characteristic analysis and area under the ROC curve (AUC) was used to evaluate the performance of each gene to predict ER, PR and HER2 status. Finally, the association of each gene (using the pre-established cutoffs) with pathological complete response (pCR) was evaluated using univariate logistic regression analyses in 2 neoadjuvant cohorts: 1) A combined HER2+ cohort using 191 tumor samples from PAMELA/PerELISA phase II trials, where patients received neoadjuvant dual HER2 blockade without chemotherapy for 15-18 weeks and 2) a combined cohort of hormone receptor-positive/HER2-negative using 205 tumor samples from IBIMA/HCB consecutive series, where patients received anthracycline/taxane-based chemotherapy. Results: Concordance between ESR1 and ER was 95.3% (95% confidence interval [95CI] 94.0-96.4%; AUC=0.98). ER+ and ER- cases were classified as ESR1- and ESR1+ in 5.5% and 4.8% of the cases, respectively. Concordance between PR and PGR was 85.2% (95CI 83.1-87.1%; AUC=0.95), and between ERBB2 and HER2 was 92.8% (95CI 91.4-94.1%; AUC=0.95). HER2+ and HER2- cases were classified ERBB2- and ERBB2+ in 16.7% and 2.9% of the cases, respectively. In the neoadjuvant HER2+ cohort, the pCR rates were 44.6% in ESR1- and 18.8% in ESR1+ (odds ratio [OR] 3.9; 95CI 1.9-8.4; p Conclusion: ESR1, PGR and ERBB2 standardized mRNA levels show high concordance with ER, PR and HER2, and might provide a more objective and quantitative prediction of response to chemotherapy and anti-HER2 therapy. The level of concordance between GE and pathology-based assessments is similar to 2 FDA 510(k) cleared pathology-based ER assays (e.g. SP1 versus 6F11 antibody clones), local versus central HER2 testing, and ER or PR status when determined by the ligand-binding assay versus immunohistochemistry. Citation Format: Pascual T, Tsai YS, Martín M, Lluch A, Cortes J, Llombart A, Conte P, Guarneri V, Rimawi MF, Alba E, Ruiz-Borrego M, Rojo F, de la Haba J, Schiff R, Adamo B, Vidal M, Paré L, Chic N, Muñoz M, Galvan P, Gonzalez-Farre B, Brauer HA, Sullivan A, Nuciforo P, Parker JS, Prat A. Standardized nCounter-based determination of estrogen receptor (ER), progesterone receptor (PR) and HER2 receptor status in breast cancer [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P2-08-06.

This abstract was not presented at the conference. Citation Format: Hardeman A, Grushko T, Clayton W, Hurley I, Khramtsova G, Parker J, Perou C, Olopade O. Not presented [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P3-07-05.

The activation of cellular senescence throughout the lifespan promotes tumor suppression, whereas the persistence of senescent cells contributes to aspects of aging. This theory has been limited, however, by an inability to identify and isolate individual senescent cells within an intact organism. Toward that end, we generated a murine reporter strain by “knocking-in” a fluorochrome, tandem-dimer Tomato (tdTom), into exon 1α of the p16INK4a locus. We used this allele (p16tdTom) for the enumeration, isolation, and characterization of individual p16INK4a-expressing cells (tdTom+). The half-life of the knocked-in transcript was shorter than that of the endogenous p16INK4a mRNA, and therefore reporter expression better correlated with p16INK4a promoter activation than p16INK4a transcript abundance. The frequency of tdTom+ cells increased with serial passage in cultured murine embryo fibroblasts from p16tdTom/+ mice. In adult mice, tdTom+ cells could be readily detected at low frequency in many tissues, and the frequency of these cells increased with aging. Using an in vivo model of peritoneal inflammation, we compared the phenotype of cells with or without activation of p16INK4a and found that tdTom+ macrophages exhibited some features of senescence, including reduced proliferation, senescence-associated β-galactosidase (SA-β-gal) activation, and increased mRNA expression of a subset of transcripts encoding factors involved in SA-secretory phenotype (SASP). These results indicate that cells harboring activation of the p16INK4a promoter accumulate with aging and inflammation in vivo, and display characteristics of senescence.

Purpose: In this era of precision-based medicine, for optimal patient care, results reported from commercial next-generation sequencing (NGS) assays should adequately reflect the burden of somatic mutations in the tumor being sequenced. Here, we sought to determine the prevalence of clonal hematopoiesis leading to possible misattribution of tumor mutation calls on unpaired Foundation Medicine NGS assays. Experimental Design: This was a retrospective cohort study of individuals undergoing NGS of solid tumors from two large cancer centers. We identified and quantified mutations in genes known to be frequently altered in clonal hematopoiesis (DNMT3A, TET2, ASXL1, TP53, ATM, CHEK2, SF3B1, CBL, JAK2) that were returned to physicians on clinical Foundation Medicine reports. For a subset of patients, we explored the frequency of true clonal hematopoiesis by comparing mutations on Foundation Medicine reports with matched blood sequencing. Results: Mutations in genes that are frequently altered in clonal hematopoiesis were identified in 65% (1,139/1,757) of patients undergoing NGS. When excluding TP53, which is often mutated in solid tumors, these events were still seen in 35% (619/1,757) of patients. Utilizing paired blood specimens, we were able to confirm that 8% (18/226) of mutations reported in these genes were true clonal hematopoiesis events. The majority of DNMT3A mutations (64%, 7/11) and minority of TP53 mutations (4%, 2/50) were clonal hematopoiesis. Conclusions: Clonal hematopoiesis mutations are commonly reported on unpaired NGS testing. It is important to recognize clonal hematopoiesis as a possible cause of misattribution of mutation origin when applying NGS findings to a patient's care. See related commentary by Pollyea, p. 5790

Purpose: Response to a complex trastuzumab-based regimen is affected by multiple features of the tumor and its microenvironment. Developing a predictive algorithm is key to optimizing HER2-targeting therapy. Experimental Design: We analyzed 137 pretreatment tumors with mRNA-seq and DNA exome sequencing from CALGB 40601, a neoadjuvant phase III trial of paclitaxel plus trastuzumab with or without lapatinib in stage II to III HER2-positive breast cancer. We adopted an Elastic Net regularized regression approach that controls for covarying features within high-dimensional data. First, we applied 517 known gene expression signatures to develop an Elastic Net model to predict pCR, which we validated on 143 samples from four independent trials. Next, we performed integrative analyses incorporating clinicopathologic information with somatic mutation status, DNA copy number alterations (CNA), and gene signatures. Results: The Elastic Net model using only gene signatures predicted pCR in the validation sets (AUC = 0.76). Integrative analyses showed that models containing gene signatures, clinical features, and DNA information were better pCR predictors than models containing a single data type. Frequently selected variables from the multiplatform models included amplifications of chromosome 6p, TP53 mutation, HER2-enriched subtype, and immune signatures. Variables predicting resistance included Luminal/ER+ features. Conclusions: Models using RNA only, as well as integrated RNA and DNA models, can predict pCR with improved accuracy over clinical variables. Somatic DNA alterations (mutation, CNAs), tumor molecular subtype (HER2E, Luminal), and the microenvironment (immune cells) were independent predictors of response to trastuzumab and paclitaxel-based regimens. This highlights the complexity of predicting response in HER2-positive breast cancer. Clin Cancer Res; 24(21); 5292–304. ©2018 AACR.

Subependymal giant-cell astrocytomas (SEGAs) are slow-growing brain tumors that are a hallmark feature seen in 5–10% of patients with Tuberous Sclerosis Complex (TSC). Though histologically benign, they can cause serious neurologic symptoms, leading to death if untreated. SEGAs consistently show biallelic loss of TSC1 or TSC2. Herein, we aimed to define other somatic events beyond TSC1/TSC2 loss and identify potential transcriptional drivers that contribute to SEGA formation. Paired tumor-normal whole-exome sequencing was performed on 21 resected SEGAs from 20 TSC patients. Pathogenic variants in TSC1/TSC2 were identified in 19/21 (90%) SEGAs. Copy neutral loss of heterozygosity (size range: 2.2–46 Mb) was seen in 76% (16/21) of SEGAs (44% chr9q and 56% chr16p). An average of 1.4 other somatic variants (range 0–7) per tumor were identified, unlikely of pathogenic significance. Whole transcriptome RNA-sequencing analyses revealed 190 common differentially expressed genes in SEGA (n = 16, 13 from a prior study) in pairwise comparison to each of: low grade diffuse gliomas (n = 530) and glioblastoma (n = 171) from The Cancer Genome Atlas (TCGA) consortium, ganglioglioma (n = 10), TSC cortical tubers (n = 15), and multiple normal tissues. Among these, homeobox transcription factors (TFs) HMX3, HMX2, VAX1, SIX3; and TFs IRF6 and EOMES were all expressed >12-fold higher in SEGAs (FDR/q-value < 0.05). Immunohistochemistry supported the specificity of IRF6, VAX1, SIX3 for SEGAs in comparison to other tumor entities and normal brain. We conclude that SEGAs have an extremely low somatic mutation rate, suggesting that TSC1/TSC2 loss is sufficient to drive tumor growth. The unique and highly expressed SEGA-specific TFs likely reflect the neuroepithelial cell of origin, and may also contribute to the transcriptional and epigenetic state that enables SEGA growth following two-hit loss of TSC1 or TSC2 and mTORC1 activation.

The E3 ubiquitin ligase Rad18 promotes a damage-tolerant and error-prone mode of DNA replication termed trans-lesion synthesis that is pathologically activated in cancer. However, the impact of vertebrate Rad18 on cancer genomes is not known. To determine how Rad18 affects mutagenesis in vivo, we have developed and implemented a novel computational pipeline to analyze genomes of carcinogen (7, 12-Dimethylbenz[a]anthracene, DMBA)-induced skin tumors from Rad18+/+ and Rad18−/− mice. We show that Rad18 mediates specific mutational signatures characterized by high levels of A(T)>T(A) single nucleotide variations (SNVs). In Rad18−/- tumors, an alternative mutation pattern arises, which is characterized by increased numbers of deletions >4 bp. Comparison with annotated human mutational signatures shows that COSMIC signature 22 predominates in Rad18+/+ tumors whereas Rad18−/− tumors are characterized by increased contribution of COSMIC signature 3 (a hallmark of BRCA-mutant tumors). Analysis of The Cancer Genome Atlas shows that RAD18 expression is strongly associated with high SNV burdens, suggesting RAD18 also promotes mutagenesis in human cancers. Taken together, our results show Rad18 promotes mutagenesis in vivo, modulates DNA repair pathway choice in neoplastic cells, and mediates specific mutational signatures that are present in human tumors.

Genomics of radiation-induced damage The potential adverse effects of exposures to radioactivity from nuclear accidents can include acute consequences such as radiation sickness, as well as long-term sequelae such as increased risk of cancer. There have been a few studies examining transgenerational risks of radiation exposure but the results have been inconclusive. Morton et al. analyzed papillary thyroid tumors, normal thyroid tissue, and blood from hundreds of survivors of the Chernobyl nuclear accident and compared them against those of unexposed patients. The findings offer insight into the process of radiation-induced carcinogenesis and characteristic patterns of DNA damage associated with environmental radiation exposure. In a separate study, Yeager et al. analyzed the genomes of 130 children and parents from families in which one or both parents had experienced gonadal radiation exposure related to the Chernobyl accident and the children were conceived between 1987 and 2002. Reassuringly, the authors did not find an increase in new germline mutations in this population. Science , this issue p. eabg2538 , p. 725

BACKGROUND Surgery remains the frontline therapy for patients with localized clear cell renal cell carcinoma (ccRCC); however, 20%–40% recur. Angiogenesis inhibitors have improved survival in metastatic patients and may result in responses in the neoadjuvant setting. The impact of these agents on the tumor genetic heterogeneity or the immune milieu is largely unknown. This phase II study was designed to evaluate safety, response, and effect on tumor tissue of neoadjuvant pazopanib. METHODS ccRCC patients with localized disease received pazopanib (800 mg daily; median 8 weeks), followed by nephrectomy. Five tumors were examined for mutations by whole exome sequencing from samples collected before therapy and at nephrectomy. These samples underwent RNA sequencing; 17 samples were available for posttreatment assessment. RESULTS Twenty-one patients were enrolled. The overall response rate was 8 of 21 (38%). No patients with progressive disease. At 1-year, response-free survival and overall survival was 83% and 89%, respectively. The most frequent grade 3 toxicity was hypertension (33%, 7 of 21). Sequencing revealed strong concordance between pre- and posttreatment samples within individual tumors, suggesting tumors harbor stable core profiles. However, a reduction in private mutations followed treatment, suggesting a selective process favoring enrichment of driver mutations. CONCLUSION Neoadjuvant pazopanib is safe and active in ccRCC. Future genomic analyses may enable the segregation of driver and passenger mutations. Furthermore, tumor infiltrating immune cells persist during therapy, suggesting that pazopanib can be combined with immune checkpoint inhibitors without dampening the immune response. FUNDING Support was provided by Novartis and GlaxoSmithKline as part of an investigator-initiated study., Neoadjuvant treatment with pazopanib results in tumor response in clear cell renal cell carcinoma, and lends insight into genetic and immunological effects of these inhibitors.

Castration-resistant prostate cancer (CRPC) is a terminal disease, demanding a better understanding of its pathogenesis. Targeted therapy needs to be developed for CRPC due to its heterogeneity and resistance to current treatments. Here, through cistrome study of YY1, a transcription factor significantly overexpressed during prostate cancer progression, we identify a YY1-PFKP axis to be essential for CRPC tumorigenesis. Depletion of YY1 in independent CRPC models dramatically reduced tumor cell growth in vitro and delayed oncogenic progression in vivo. Importantly, YY1 functions as a master regulator of prostate tumor metabolism including the Warburg effect and mitochondria respiration. Loss-of-function and rescue studies further reveals a mechanistic underpinning in which YY1 directly binds and trans-activates PFKP, a gene encoding the rate-limiting enzyme for glycolysis, significantly contributing to the YY1-enforced oncogenic phenotypes such as enhanced tumor cell glycolysis and malignant growth. Additionally, a vast majority of gene-regulatory element in advanced prostate cancer cells are bound by YY1, lending a support for its role as a master regulator of prostate cancer progression. YY1 interactome studies point to bromodomain-containing coactivators in prostate cancer, which act as functional partners of YY1 to potentiate YY1-related target gene activation. Altogether, this study unveils an unexplored YY1:BRD4-PFKP oncogenic axis operating in advanced prostate cancer with implications for therapy.

Purpose Urachal adenocarcinoma is a rare type of primary bladder adenocarcinoma that comprises less than 1% of all bladder cancers. The low incidence of urachal adenocarcinomas does not allow for an evidence-based approach to therapy. Transcriptome profiling of urachal adenocarcinomas has not been previously reported. We hypothesized that an in-depth molecular understanding of urachal adenocarcinoma would uncover rational therapeutic strategies. Patients and Methods We performed targeted exon sequencing and global transcriptome profiling of 12 urachal tumors to generate a comprehensive molecular portrait of urachal adenocarcinoma. A single patient with an MSH6 mutation was treated with the anti–programmed death-ligand 1 antibody, atezolizumab. Results Urachal adenocarcinoma closely resembles colorectal cancer at the level of RNA expression, which extends previous observations that urachal tumors harbor genomic alterations that are found in colorectal adenocarcinoma. A subset of tumors was found to have alterations in genes that are associated with microsatellite instability ( MSH2 and MSH6) and hypermutation ( POLE). A patient with an MSH6 mutation was treated with immune checkpoint blockade, which resulted in stable disease. Conclusion Because clinical trials are next to impossible for patients with rare tumors, precision oncology may be an important adjunct for treatment decisions. Our findings demonstrate that urachal adenocarcinomas molecularly resemble colorectal adenocarcinomas at the level of RNA expression, are the first report, to our knowledge, of MSH2 and MSH6 mutations in this disease, and support the consideration of immune checkpoint blockade as a rational therapeutic treatment of this exceedingly rare tumor.

Introduction: Studies on the ubiquitous airway protein bactericidal/permeability-increasing protein fold-containing family member A1 (BPIFA1) have revealed important roles in reducing eosinophilic airway inflammation, relaxing airway smooth muscle, and defending against microbial infection. Aims and Objectives: The objective of this study was to test the hypothesis that BPIFA1 expression is reduced in asthmatics. Methods: BPIFA1 protein levels were measured by ELISA from bronchoalveolar lavage (BAL) supernatants from twenty mild asthmatics and twenty normal volunteers who underwent research bronchoscopy at our center. BPIFA1 gene expression in airway epithelia was analyzed from seven transcriptomic studies involving 355 asthmatics and 193 controls (Tsai et al. ERJ. 2018;51). Results: The concentration of BPIFA1 was significantly lower in BAL samples from asthmatics compared to normal volunteers (9 + 3 vs. 38 + 8 ng/mL, p Conclusions: BPIFA1 protein and BPIFA1 mRNA expression is reduced in the airways and airway epithelia of asthmatics, respectively, suggesting that reduced protein levels are due to down-regulation of expression rather than degradation. Further research characterizing BPIFA1 levels in larger numbers of asthmatics with well-defined phenotypes will help determine if BPIFA1-deficiency contributes to asthma pathogenesis.

Colonization of the lung by Yersinia pestis is a critical first step in establishing infection during primary pneumonic plague, a disease characterized by high lethality. However, the mechanisms by which Y. pestis adheres in the lung after inhalation remain elusive. Here, we used Tn-seq to identify Y. pestis genes important for adherence early during primary pneumonic plague. Our mutant enrichment strategy resulted in the identification of genes important for regulation and assembly of genes and proteins rather than adhesin genes themselves. These results reveal that there may be multiple Y. pestis adhesins or redundancy among adhesins. Identifying the adhesins regulated by the genes identified in our enrichment screen may reveal novel therapeutic targets for preventing Y. pestis adherence and the subsequent development of pneumonic plague., Following inhalation, Yersinia pestis rapidly colonizes the lung to establish infection during primary pneumonic plague. Although several adhesins have been identified in Yersinia spp., the factors mediating early Y. pestis adherence in the lung remain unknown. To identify genes important for Y. pestis adherence during primary pneumonic plague, we used transposon insertion sequencing (Tn-seq). Wild-type and capsule mutant (Δcaf1) Y. pestis transposon mutant libraries were serially passaged in vivo to enrich for nonadherent mutants in the lung using a mouse model of primary pneumonic plague. Sequencing of the passaged libraries revealed six mutants that were significantly enriched in both the wild-type and Δcaf1 Y. pestis backgrounds. The enriched mutants had insertions in genes that encode transcriptional regulators, chaperones, an endoribonuclease, and YPO3903, a hypothetical protein. Using single-strain infections and a transcriptional analysis, we identified a significant role for YPO3903 in Y. pestis adherence in the lung and showed that YPO3903 regulated transcript levels of psaA, which encodes a fimbria previously implicated in Y. pestis adherence in vitro. Deletion of psaA had a minor effect on Y. pestis adherence in the lung, suggesting that YPO3903 regulates other adhesins in addition to psaA. By enriching for mutations in genes that regulate the expression or assembly of multiple genes or proteins, we obtained screen results indicating that there may be not just one dominant adhesin but rather several factors that contribute to early Y. pestis adherence during primary pneumonic plague. IMPORTANCE Colonization of the lung by Yersinia pestis is a critical first step in establishing infection during primary pneumonic plague, a disease characterized by high lethality. However, the mechanisms by which Y. pestis adheres in the lung after inhalation remain elusive. Here, we used Tn-seq to identify Y. pestis genes important for adherence early during primary pneumonic plague. Our mutant enrichment strategy resulted in the identification of genes important for regulation and assembly of genes and proteins rather than adhesin genes themselves. These results reveal that there may be multiple Y. pestis adhesins or redundancy among adhesins. Identifying the adhesins regulated by the genes identified in our enrichment screen may reveal novel therapeutic targets for preventing Y. pestis adherence and the subsequent development of pneumonic plague.

SUMMARY Mucoepidermoid carcinoma (MEC) is a life-threatening salivary gland cancer that is driven primarily by a transcriptional coactivator fusion composed of cyclic AMP-regulated transcriptional coactivator 1 (CRTC1) and mastermind-like 2 (MAML2). The mechanisms by which the chimeric CRTC1/MAML2 (C1/M2) oncoprotein rewires gene expression programs that promote tumorigenesis remain poorly understood. Here, we show that C1/M2 induces transcriptional activation of the non-canonical peroxisome proliferator-activated receptor gamma coactivator-1 alpha (PGC-1α) splice variant PGC-1α4, which regulates peroxisome proliferator-activated receptor gamma (PPARγ)-mediated insulin-like growth factor 1 (IGF-1) expression. This mitogenic transcriptional circuitry is consistent across cell lines and primary tumors. C1/M2-positive tumors exhibit IGF-1 pathway activation, and small-molecule drug screens reveal that tumor cells harboring the fusion gene are selectively sensitive to IGF-1 receptor (IGF-1R) inhibition. Furthermore, this dependence on autocrine regulation of IGF-1 transcription renders MEC cells susceptible to PPARγ inhibition with inverse agonists. These results yield insights into the aberrant coregulatory functions of C1/M2 and identify a specific vulnerability that can be exploited for precision therapy., Graphical Abstract, In brief Musicant et al. demonstrate that the CRTC1-MAML2 gene fusion aberrantly regulates expression of the minor splice variant PGC-1α4, which in turn coactivates PPARγ to induce IGF-1 expression. This synthetic signal circuit establishes pro-growth and pro-survival signaling in mucoepidermoid carcinomas, sensitizing tumors to drugs that disable IGF-1 signaling by targeting PPARγ.

SUMMARYMucoepidermoid carcinoma (MEC) is a life-threatening salivary gland cancer that is driven primarily by the transcriptional co-activator fusion CRTC1-MAML2. The mechanisms by which the chimeric CRTC1-MAML2 oncoprotein rewires gene expression programs that promote tumorigenesis remain poorly understood. Here, we show that CRTC1-MAML2 induces transcriptional activation of the non-canonical PGC-1α splice variant PGC-1α4, which regulates PPARγ-dependent IGF-1 expression. This mitogenic transcriptional circuitry is consistent across cell lines and primary tumors. CRTC1-MAML2 positive tumors are dominated by IGF-1 pathway activation and small molecule drug screens reveal that tumor cells harboring the fusion gene are selectively sensitive to IGF-1R inhibition. Furthermore, this dependence on autocrine regulation of IGF-1 transcription renders MEC cells susceptible to PPARγinhibition with inverse agonists. These results yield insights into the aberrant co-regulatory functions of CRTC1-MAML2 and identify a specific vulnerability that can be exploited for precision therapy.

Purpose: Irinotecan/5-fluorouracil (5-FU; FOLFIRI) or oxaliplatin/5-FU (FOLFOX), combined with bevacizumab or cetuximab, are approved, first-line treatments for metastatic colorectal cancer (mCRC). We aimed at identifying germline variants associated with survival in patients with mCRC treated with these regimens in Cancer and Leukemia Group B/SWOG 80405. Experimental Design: Patients with mCRC receiving either FOLFOX or FOLFIRI were randomized to either cetuximab or bevacizumab. DNA from peripheral blood was genotyped for approximately 700,000 SNPs. The association between SNPs and overall survival (OS) was tested in 613 patients of genetically estimated European ancestry using Cox proportional hazards models. Results: The four most significant SNPs associated with OS were three haplotypic SNPs between microsomal glutathione S-transferase 1 (MGST1) and LIM domain only 3 (LMO3, representative HR, 1.56; P = 1.30 × 10−6), and rs11644916 in AXIN1 (HR, 1.39, P = 4.26 × 10−6). AXIN1 is a well-established tumor suppressor gene in colorectal cancer, and rs11644916 (G>A) conferred shorter OS. Median OS for patients with the AA, AG, or GG genotypes was 18.4, 25.6, or 36.4 months, respectively. In 90 patients with stage IV colorectal cancer from The Cancer Genome Atlas (TCGA), rs11649255 in AXIN1 [in almost complete linkage disequilibrium (LD) with rs11644916], was associated with shorter OS (HR, 2.24, P = 0.0096). Using rs11648673 in AXIN1 (in very high LD with rs11644916 and with functional evidence), luciferase activity in three colorectal cancer cell lines was reduced. Conclusions: This is the first large genome-wide association study ever conducted in patients with mCRC treated with first-line standard treatment in a randomized phase III trial. A common SNP in AXIN1 conferred worse OS and the effect was replicated in TCGA. Further studies in colorectal cancer experimental models are required.

Small cell carcinoma of the ovary, hypercalcemic type (SCCOHT) is a rare and aggressive form of ovarian cancer. SCCOHT tumors have inactivating mutations in SMARCA4 (BRG1), one of the two mutually exclusive ATPases of the SWI/SNF chromatin remodeling complex. To address the role that BRG1 loss plays in SCCOHT tumorigenesis, we performed integrative multi-omic analyses in SCCOHT cell lines +/- BRG1 reexpression. BRG1 reexpression induced a gene and protein signature similar to an epithelial cell and gained chromatin accessibility sites correlated with other epithelial originating TCGA tumors. Gained chromatin accessibility and BRG1 recruited sites were strongly enriched for transcription-factor-binding motifs of AP-1 family members. Furthermore, AP-1 motifs were enriched at the promoters of highly upregulated epithelial genes. Using a dominant-negative AP-1 cell line, we found that both AP-1 DNA-binding activity and BRG1 reexpression are necessary for the gene and protein expression of epithelial genes. Our study demonstrates that BRG1 reexpression drives an epithelial-like gene and protein signature in SCCOHT cells that depends upon by AP-1 activity.

Background Recent molecular characterization of urothelial cancer (UC) has suggested potential pathways in which to direct treatment, leading to a host of targeted therapies in development for UC. In parallel, gene expression profiling has demonstrated that high-grade UC is a heterogeneous disease. Prognostic basal-like and luminal-like subtypes have been identified and an accurate transcriptome BASE47 classifier has been developed. However, these phenotypes cannot be broadly investigated due to the lack of a clinically viable diagnostic assay. We sought to develop and evaluate a diagnostic classifier of UC subtype with the goal of accurate classification from clinically available specimens. Methods Tumor samples from 52 patients with high-grade UC were profiled for BASE47 genes concurrently by RNAseq as well as NanoString. After design and technical validation of a BASE47 NanoString probeset, results from the RNAseq and NanoString were used to translate diagnostic criteria to the Nanostring platform. Evaluation of repeatability and accuracy was performed to derive a final Nanostring based classifier. Diagnostic classification resulting from the NanoString BASE47 classifier was validated on an independent dataset (n = 30). The training and validation datasets accurately classified 87% and 93% of samples, respectively. Results Here we have derived a NanoString-platform BASE47 classifier that accurately predicts basal-like and luminal-like subtypes in high grade urothelial cancer. We have further validated our new NanoString BASE47 classifier on an independent dataset and confirmed high accuracy when compared with our original Transcriptome BASE47 classifier. Conclusions The NanoString BASE47 classifier provides a faster turnaround time, a lower cost per sample to process, and maintains the accuracy of the original subtype classifier for better clinical implementation.

T-cell receptor (TCR) repertoire profiling has emerged as a powerful tool for biological discovery and biomarker development in cancer immunology and immunotherapy. A key statistic derived from repertoire profiling data is diversity, which summarizes the frequency distribution of TCRs within a mixed population. Despite the growing use of TCR diversity metrics in clinical trial correlative studies in oncology, their accuracy has not been validated using published ground-truth datasets. Here, we reported the performance characteristics of methods for TCR repertoire profiling from RNA-sequencing data, showed undersampling as a prominent source of bias in diversity estimates, and derived a model via statistical learning that attenuates bias to produce corrected diversity estimates. This modeled diversity improved discrimination in The Cancer Genome Atlas data and associated with survival and treatment response in patients with melanoma treated with anti–PD-1 therapy, where the commonly used diversity normalizations did not. These findings have the potential to increase our understanding of the tumor immune microenvironment and improve the accuracy of predictions of patient responses to immunotherapy.

Background In early-stage HER2-positive breast cancer, escalation or de-escalation of systemic therapy is a controversial topic. As an aid to treatment decisions, we aimed to develop a prognostic assay that integrates multiple data types for predicting survival outcome in patients with newly diagnosed HER2-positive breast cancer. Methods We derived a combined prognostic model using retrospective clinical-pathological data on stromal tumourinfiltrating lymphocytes, PAM50 subtypes, and expression of 55 genes obtained from patients who participated in the Short-HER phase 3 trial. The trial enrolled patients with newly diagnosed, node-positive, HER2-positive breast cancer or, if node negative, with at least one risk factor (ie, tumour size >2 cm, histological grade 3, lymphovascular invasion, Ki67 >20%, age

Background Identification of HER2-positive breast cancers with high anti-HER2 sensitivity could help de-escalate chemotherapy. Here, we tested a clinically applicable RNA-based assay that combines ERBB2 and the HER2-enriched (HER2-E) intrinsic subtype in HER2-positive disease treated with dual HER2-blockade without chemotherapy. Methods A research-based PAM50 assay was applied in 422 HER2-positive tumors from five II–III clinical trials (SOLTI-PAMELA, TBCRC023, TBCRC006, PER-ELISA, EGF104090). In SOLTI-PAMELA, TBCRC023, TBCRC006, and PER-ELISA, all patients had early disease and were treated with neoadjuvant lapatinib or pertuzumab plus trastuzumab for 12–24 weeks. Primary outcome was pathological complete response (pCR). In EGF104900, 296 women with advanced disease were randomized to receive either lapatinib alone or lapatinib plus trastuzumab. Progression-free survival (PFS), overall response rate (ORR), and overall survival (OS) were evaluated. Results A total of 305 patients with early and 117 patients with advanced HER2-positive disease were analyzed. In early disease, HER2-E represented 83.8% and 44.7% of ERBB2-high and ERBB2-low tumors, respectively. Following lapatinib and trastuzumab, the HER2-E and ERBB2 (HER2-E/ERBB2)-high group showed a higher pCR rate compared to the rest (44.5%, 95% confidence interval [CI] = 35.4% to 53.9% vs 11.6%, 95% CI = 6.9% to 18.0%; adjusted odds ratio [OR] = 6.05, 95% CI = 3.10 to 11.80, P Conclusions Combining HER2-E subtype and ERBB2 mRNA into a single assay identifies tumors with high responsiveness to HER2-targeted therapy. This biomarker could help de-escalate chemotherapy in approximately 40% of patients with HER2-positive breast cancer.

Male breast cancer (BC) is rare, representing approximately 1% of cancers that occur in men and approximately 1% of all BCs worldwide. Because male BC is rare, not much is known about the disease, and treatment recommendations are typically extrapolated from data available from clinical trials enrolling female BC patients. We review the epidemiology, risk factors, prognosis, and the varied molecular and clinicopathologic features that characterize male BC. In addition, we summarize the available data for the use of systemic therapy in the treatment of male BC and explore the ongoing development of targeted therapeutic agents for the treatment of this subgroup of BCs. There are important biological differences between male and female BC. Male BC is almost exclusively hormone receptor positive (+), including the androgen receptor (AR), and is associated with an increased prevalence of BRCA2 germline mutations, especially in men with increased risk for developing high-risk BC. Additional research is warranted to better characterize male BC. To accomplish this, a multi-national consortium approach, such as the International Male Breast Cancer Program, is needed in response to the scarcity of patients. This approach allows the pooling of information from a large number of men with BC and the creation of registries for future therapeutic-focused clinical trials. Given the unique biology of BC in men, promising new therapeutic targets are currently under investigation, including the use of poly-ADP-ribose polymerase inhibitors or AR-targeted agents either as monotherapy or in combination with other agents.

The ability to characterize and predict tumor phenotypes is crucial to precision medicine. In this study, we present an integrative computational approach using a genome-wide association analysis and an Elastic Net prediction method to analyze the relationship between DNA copy number alterations and an archive of gene expression signatures. Across breast cancers, we are able to quantitatively predict many gene signatures levels within individual tumors with high accuracy based upon DNA copy number features alone, including proliferation status and Estrogen-signaling pathway activity. We can also predict many other key phenotypes, including intrinsic molecular subtypes, estrogen receptor status, and TP53 mutation. This approach is also applied to TCGA Pan-Cancer, which identify repeatedly predictable signatures across tumor types including immune features in lung squamous and basal-like breast cancers. These Elastic Net DNA predictors could also be called from DNA-based gene panels, thus facilitating their use as biomarkers to guide therapeutic decision making., Effective precision medicine strategies rely on the ability to predict tumour behaviour based on molecular characteristics. Here, the authors build models to predict multiple distinct gene expression patterns using DNA copy number alterations

Previous studies have shown that common variants of the gene coding for FK506-binding protein 51 (FKBP5), a critical regulator of glucocorticoid sensitivity, affect vulnerability to stress-related disorders. In a previous report,FKBP5rs1360780 was identified as a functional variant because of its effect on gene methylation. Here we report evidence for a novel functionalFKBP5allele, rs3800373. This study assessed the association between rs3800373 and post-traumatic chronic pain in 1607 women and men from two ethnically diverse human cohorts. The molecular mechanism through which rs3800373 affects adverse outcomes was established viain silico,in vivo, andin vitroanalyses. The rs3800373 minor allele predicted worse adverse outcomes after trauma exposure, such that individuals with the minor (risk) allele developed more severe post-traumatic chronic musculoskeletal pain. Among these individuals, peritraumatic circulatingFKBP5expression levels increased as cortisol and glucocorticoid receptor (NR3C1) mRNA levels increased, consistent with increased glucocorticoid resistance. Bioinformatic,in vitro, and mutational analyses indicate that the rs3800373 minor allele reduces the binding of a stress- and pain-associated microRNA, miR-320a, toFKBP5via altering theFKBP5mRNA 3′UTR secondary structure (i.e., is a riboSNitch). This results in relatively greaterFKBP5translation, unchecked by miR-320a. Overall, these results identify an important gene–miRNA interaction influencing chronic pain risk in vulnerable individuals and suggest that exogenous methods to achieve targeted reduction in poststressFKBP5mRNA expression may constitute useful therapeutic strategies.SIGNIFICANCE STATEMENTFKBP5is a critical regulator of the stress response. Previous studies have shown that dysregulation of the expression of this gene plays a role in the pathogenesis of chronic pain development as well as a number of comorbid neuropsychiatric disorders. In the current study, we identified a functional allele (rs3800373) in the 3′UTR ofFKBP5that influences vulnerability to chronic post-traumatic pain in two ethnic cohorts. Using multiple complementary experimental approaches, we show that theFKBP5rs3800373 minor allele alters the secondary structure ofFKBP5mRNA, decreasing the binding of a stress- and pain-associated microRNA, miR-320a. This results in relatively greaterFKBP5translation, unchecked by miR-320a, increasing glucocorticoid resistance and increasing vulnerability to post-traumatic pain.