In silico APC/C substrate discovery reveals cell cycle-dependent degradation of UHRF1 and other chromatin regulators.

The anaphase-promoting complex/cyclosome (APC/C) is an E3 ubiquitin ligase and critical regulator of cell cycle progression. Despite its vital role, it has remained challenging to globally map APC/C substrates. By combining orthogonal features of known substrates, we predicted APC/C substrates in silico. This analysis identified many known substrates and suggested numerous candidates. Unexpectedly, chromatin regulatory proteins are enriched among putative substrates, and we show experimentally that several chromatin proteins bind APC/C, oscillate during the cell cycle, and are degraded following APC/C activation, consistent with being direct APC/C substrates. Additional analysis revealed detailed mechanisms of ubiquitylation for UHRF1, a key chromatin regulator involved in histone ubiquitylation and DNA methylation maintenance. Disrupting UHRF1 degradation at mitotic exit accelerates G1-phase cell cycle progression and perturbs global DNA methylation patterning in the genome. We conclude that APC/C coordinates crosstalk between cell cycle and chromatin regulatory proteins. This has potential consequences in normal cell physiology, where the chromatin environment changes depending on proliferative state, as well as in disease.

bcSeq: an R package for fast sequence mapping in high-throughput shRNA and CRISPR screens.

Summary: CRISPR-Cas9 and shRNA high-throughput sequencing screens have abundant applications for basic and translational research. Methods and tools for the analysis of these screens must properly account for sequencing error, resolve ambiguous mappings among similar sequences in the barcode library in a statistically principled manner, and be computationally efficient. Herein we present bcSeq, an open source R package that implements a fast and parallelized algorithm for mapping high-throughput sequencing reads to a barcode library while tolerating sequencing error. The algorithm uses a Trie data structure for speed and resolves ambiguous mappings by using a statistical sequencing error model based on Phred scores for each read. Availability and implementation: The package source code and an accompanying tutorial are available at http://bioconductor.org/packages/bcSeq/. Supplementary information: Supplementary data are available at Bioinformatics online.

Intrinsic subtypes of high-grade bladder cancer reflect the hallmarks of breast cancer biology.

We sought to define whether there are intrinsic molecular subtypes of high-grade bladder cancer. Consensus clustering performed on gene expression data from a meta-dataset of high-grade, muscle-invasive bladder tumors identified two intrinsic, molecular subsets of high-grade bladder cancer, termed "luminal" and "basal-like," which have characteristics of different stages of urothelial differentiation, reflect the luminal and basal-like molecular subtypes of breast cancer, and have clinically meaningful differences in outcome. A gene set predictor, bladder cancer analysis of subtypes by gene expression (BASE47) was defined by prediction analysis of microarrays (PAM) and accurately classifies the subtypes. Our data demonstrate that there are at least two molecularly and clinically distinct subtypes of high-grade bladder cancer and validate the BASE47 as a subtype predictor. Future studies exploring the predictive value of the BASE47 subtypes for standard of care bladder cancer therapies, as well as novel subtype-specific therapies, will be of interest.

Spatial Relationships in the Tumor Microenvironment Demonstrate Association with Pathologic Response to Neoadjuvant Chemoimmunotherapy in Muscle-invasive Bladder Cancer.

BACKGROUND: Platinum-based neoadjuvant chemotherapy (NAC) is standard for patients with muscle-invasive bladder cancer (MIBC). Pathologic response (complete: ypT0N0 and partial:

Poor Prognosis among Radiation-Associated Bladder Cancer Is Defined by Clinicogenomic Features.

Radiotherapy (RT) for prostate cancer has been associated with an increased risk for the development of bladder cancer. We aimed to integrate clinical and genomic data to better understand the development of RT-associated bladder cancer. A retrospective analysis was performed to identify control patients (CTRL; n = 41) and patients with RT-associated bladder cancer (n = 41). RT- and CTRL-specific features were then identified through integration and analysis of the genomic sequencing data and clinical variables. RT-associated bladder tumors were significantly enriched for alterations in KDM6A and ATM, whereas CTRL tumors were enriched for CDKN2A mutation. Globally, there were an increased number of variants within RT tumors, albeit at a lower variant allele frequency. Mutational signature analysis revealed three predominate motif patterns, with similarity to SBS2/13 (APOBEC3A), SBS5 (ERCC2/smoking), and SBS6/15 (MMR). Poor prognostic factors in the RT cohort include a short tumor latency, smoking status, the presence of the smoking and X-ray therapy mutational signatures, and CDKN2A copy number loss. Based on the clinical and genomic findings, we suggest at least two potential pathways leading to RT-associated bladder cancer: The first occurs in the setting of field cancerization related to smoking or preexisting genetic alterations and leads to the development of more aggressive bladder tumors, and the second involves RT initiating the oncogenic process in otherwise healthy urothelium, leading to a longer latency and less aggressive disease. SIGNIFICANCE: Clinicogenomic analysis of radiation-associated bladder cancer uncovered mutational signatures that, in addition to a short tumor latency, smoking, and CDKN2A loss, are associated with a poor outcome. These clinical and genomic features provide a potential method to identify patients with prostate cancer who are at an increased risk for the development of aggressive bladder cancer following prostate RT.

Determination of permissive and restraining cancer-associated fibroblast (DeCAF) subtypes.

Cancer-associated fibroblast (CAF) subpopulations in pancreatic ductal adenocarcinoma (PDAC) have been identified using single-cell RNA sequencing (scRNAseq) with divergent characteristics, but their clinical relevance remains unclear. We translate scRNAseq-derived CAF cell-subpopulation-specific marker genes to bulk RNAseq data, and develop a single- sample classifier, DeCAF, for the classification of clinically rest raining and perm issive CAF subtypes. We validate DeCAF in 19 independent bulk transcriptomic datasets across four tumor types (PDAC, mesothelioma, bladder and renal cell carcinoma). DeCAF subtypes have distinct histology features, immune landscapes, and are prognostic and predict response to therapy across cancer types. We demonstrate that DeCAF is clinically replicable and robust for the classification of CAF subtypes in patients for multiple tumor types, providing a better framework for the future development and translation of therapies against permissive CAF subtypes and preservation of restraining CAF subtypes. Significance: We introduce a replicable and robust classifier, DeCAF, that delineates the significance of the role of permissive and restraining CAF subtypes in cancer patients. DeCAF is clinically tractable, prognostic and predictive of treatment response in multiple cancer types and lays the translational groundwork for the preclinical and clinical development of CAF subtype specific therapies.

Immune features are associated with response to neoadjuvant chemo-immunotherapy for muscle-invasive bladder cancer.

Neoadjuvant cisplatin-based chemotherapy is standard of care for muscle-invasive bladder cancer (MIBC). Immune checkpoint inhibition (ICI) alone, and ICI in combination with chemotherapy, have demonstrated promising pathologic response (

FGFR inhibition augments anti-PD-1 efficacy in murine FGFR3-mutant bladder cancer by abrogating immunosuppression.

The combination of targeted therapy with immune checkpoint inhibition (ICI) is an area of intense interest. We studied the interaction of fibroblast growth factor receptor (FGFR) inhibition with ICI in urothelial carcinoma (UC) of the bladder, in which FGFR3 is altered in 50% of cases. Using an FGFR3-driven, Trp53-mutant genetically engineered murine model (UPFL), we demonstrate that UPFL tumors recapitulate the histology and molecular subtype of their FGFR3-altered human counterparts. Additionally, UPFL1 allografts exhibit hyperprogression to ICI associated with an expansion of T regulatory cells (Tregs). Erdafitinib blocked Treg proliferation in vitro, while in vivo ICI-induced Treg expansion was fully abrogated by FGFR inhibition. Combined erdafitinib and ICI resulted in high therapeutic efficacy. In aggregate, our work establishes that, in mice, co-alteration of FGFR3 and Trp53 results in high-grade, non-muscle-invasive UC and presents a previously underappreciated role for FGFR inhibition in blocking ICI-induced Treg expansion.

DNA Damage Repair Classifier Defines Distinct Groups in Hepatocellular Carcinoma.

DNA repair pathways have been associated with variability in hepatocellular carcinoma (HCC) clinical outcomes, but the mechanism through which DNA repair varies as a function of liver regeneration and other HCC characteristics is poorly understood. We curated a panel of 199 genes representing 15 DNA repair pathways to identify DNA repair expression classes and evaluate their associations with liver features and clinicopathologic variables in The Cancer Genome Atlas (TCGA) HCC study. We identified two groups in HCC, defined by low or high expression across all DNA repair pathways. The low-repair group had lower grade and retained the expression of classical liver markers, whereas the high-repair group had more clinically aggressive features, increased p53 mutant-like gene expression, and high liver regenerative gene expression. These pronounced features overshadowed the variation in the low-repair subset, but when considered separately, the low-repair samples included three subgroups: L1, L2, and L3. L3 had high DNA repair expression with worse progression-free (HR 1.24, 95% CI 0.81-1.91) and overall (HR 1.63, 95% CI 0.98-2.71) survival. High-repair outcomes were also significantly worse compared with the L1 and L2 groups. HCCs vary in DNA repair expression, and a subset of tumors with high regeneration profoundly disrupts liver biology and poor prognosis.

Activation of the CREB Coactivator CRTC2 by Aberrant Mitogen Signaling promotes oncogenic functions in HPV16 positive head and neck cancer.

Head and neck squamous cell carcinoma (HNSCC) is the 6th most common cancer worldwide and incidence rates are continuing to rise globally. Patients often present with locally advanced disease and a staggering 50% chance of relapse following treatment. Aberrant activation of adaptive response signaling pathways, such as the cAMP/PKA pathway, induce an array of genes associated with known cancer pathways that promote tumorigenesis and drug resistance. We identified the cAMP Regulated Transcription Coactivator 2 (CRTC2) to be overexpressed and constitutively activated in HNSCCs and this confers poor prognosis. CRTCs are regulated through their subcellular localization and we show that CRTC2 is exclusively nuclear in HPV(+) HNSCC, thus constitutively active, due to non-canonical Mitogen-Activated Kinase Kinase 1 (MEKK1)-mediated activation via a MEKK1-p38 signaling axis. Loss-of-function and pharmacologic inhibition experiments decreased CRTC2/CREB transcriptional activity by reducing nuclear CRTC2 via nuclear import inhibition and/or by eviction of CRTC2 from the nucleus. This shift in localization was associated with decreased proliferation, migration, and invasion. Our results suggest that small molecules that inhibit nuclear CRTC2 and p38 activity may provide therapeutic benefit to patients with HPV(+) HNSCC.

Collaborative study from the Bladder Cancer Advocacy Network for the genomic analysis of metastatic urothelial cancer.

Urothelial Cancer - Genomic Analysis to Improve Patient Outcomes and Research (NCT02643043), UC-GENOME, is a genomic analysis and biospecimen repository study in 218 patients with metastatic urothelial carcinoma. Here we report on the primary outcome of the UC-GENOME-the proportion of subjects who received next generation sequencing (NGS) with treatment options-and present the initial genomic analyses and clinical correlates. 69.3% of subjects had potential treatment options, however only 5.0% received therapy based on NGS. We found an increased frequency of TP53E285K mutations as compared to non-metastatic cohorts and identified features associated with benefit to chemotherapy and immune checkpoint inhibition, including: Ba/Sq and Stroma-rich subtypes, APOBEC mutational signature (SBS13), and inflamed tumor immune phenotype. Finally, we derive a computational model incorporating both genomic and clinical features predictive of immune checkpoint inhibitor response. Future work will utilize the biospecimens alongside these foundational analyses toward a better understanding of urothelial carcinoma biology.

Genomic characterization of rare molecular subclasses of hepatocellular carcinoma.

Primary liver cancer, consisting of both cholangiocarcinoma (CCA) and hepatocellular carcinoma (HCC), is the second leading cause of cancer deaths worldwide. Our goal is to genomically characterize rare HCC subclasses to provide insight into disease biology. Leveraging The Cancer Genome Atlas (TCGA) to perform a combined analysis of CCA (n = 36) and HCC (n = 275), we integrated multiple genomic platforms, to assess transcriptional profiles, mutational signatures, and copy number patterns to uncover underlying etiology and linage specific patterns. We identified two molecular classes distinct from prototypical HCC tumors. The first, CCA-Like, although histologically indistinguishable from HCC, had enrichment of CCA mutations (IDH1, BAP1), mutational signatures, and transcriptional patterns (SOX9, KRT19). CCA-Like, however, retained a copy number landscape similar to HCC, suggesting a hepatocellular linage. The second, Blast-Like, is enriched in TP53 mutations, HBV infection, exposure related mutational signatures and transcriptionally similar to hepatoblasts. Although these subclasses are molecularly distinct, they both have a worse progression-free survival compared to classical HCC tumors, yet are clinically treated the same. The identification of and characterization of CCA-Like and Blast-Like subclasses advance our knowledge of HCC as well as represents an urgent need for the identification of class specific biomarkers and targeted therapy.

Analytical protocol to identify local ancestry-associated molecular features in cancer.

People of different ancestries vary in cancer risk and outcome, and their molecular differences may indicate sources of these variations. Determining the "local" ancestry composition at each genetic locus across ancestry-admixed populations can suggest causal associations. We present a protocol to identify local ancestry and detect the associated molecular changes, using data from the Cancer Genome Atlas. This workflow can be applied to cancer cohorts with matched tumor and normal data from admixed patients to examine germline contributions to cancer. For complete details on the use and execution of this protocol, please refer to Carrot-Zhang et al. (2020).

Integrative modeling identifies genetic ancestry-associated molecular correlates in human cancer.

Cellular and molecular aberrations contribute to the disparity of human cancer incidence and etiology between ancestry groups. Multiomics profiling in The Cancer Genome Atlas (TCGA) allows for querying of the molecular underpinnings of ancestry-specific discrepancies in human cancer. Here, we provide a protocol for integrative associative analysis of ancestry with molecular correlates, including somatic mutations, DNA methylation, mRNA transcription, miRNA transcription, and pathway activity, using TCGA data. This protocol can be generalized to analyze other cancer cohorts and human diseases. For complete details on the use and execution of this protocol, please refer to Carrot-Zhang et al. (2020).

Identification of a Novel Inflamed Tumor Microenvironment Signature as a Predictive Biomarker of Bacillus Calmette-Guérin Immunotherapy in Non-Muscle-Invasive Bladder Cancer.

PURPOSE: Improved risk stratification and predictive biomarkers of treatment response are needed for non-muscle-invasive bladder cancer (NMIBC). Here we assessed the clinical utility of targeted RNA and DNA molecular profiling in NMIBC. EXPERIMENTAL DESIGN: Gene expression in NMIBC samples was profiled by NanoString nCounter, an RNA quantification platform, from two independent cohorts (n = 28, n = 50); targeted panel sequencing was performed in a subgroup (n = 50). Gene signatures were externally validated using two RNA sequencing datasets of NMIBC tumors (n = 438, n = 73). Established molecular subtype classifiers and novel gene expression signatures were assessed for associations with clinicopathologic characteristics, somatic tumor mutations, and treatment outcomes. RESULTS: Molecular subtypes distinguished between low-grade Ta tumors with FGFR3 mutations and overexpression (UROMOL-class 1) and tumors with more aggressive clinicopathologic characteristics (UROMOL-classes 2 and 3), which were significantly enriched with TERT promoter mutations. However, UROMOL subclasses were not associated with recurrence after bacillus Calmette-Guérin (BCG) immunotherapy in two independent cohorts. In contrast, a novel expression signature of an inflamed tumor microenvironment (TME) was associated with improved recurrence-free survival after BCG. Expression of immune checkpoint genes (PD-L1/PD-1/CTLA-4) was associated with an inflamed TME, but not with higher recurrence rates after BCG. FGFR3 mutations and overexpression were both associated with low immune signatures. CONCLUSIONS: Assessment of the immune TME, rather than molecular subtypes, is a promising predictive biomarker of BCG response. Modulating the TME in an immunologically "cold" tumor warrants further investigation. Integrated transcriptomic and exome sequencing should improve treatment selection in NMIBC.

Whole-genome characterization of lung adenocarcinomas lacking the RTK/RAS/RAF pathway.

RTK/RAS/RAF pathway alterations (RPAs) are a hallmark of lung adenocarcinoma (LUAD). In this study, we use whole-genome sequencing (WGS) of 85 cases found to be RPA(-) by previous studies from The Cancer Genome Atlas (TCGA) to characterize the minority of LUADs lacking apparent alterations in this pathway. We show that WGS analysis uncovers RPA(+) in 28 (33%) of the 85 samples. Among the remaining 57 cases, we observe focal deletions targeting the promoter or transcription start site of STK11 (n = 7) or KEAP1 (n = 3), and promoter mutations associated with the increased expression of ILF2 (n = 6). We also identify complex structural variations associated with high-level copy number amplifications. Moreover, an enrichment of focal deletions is found in TP53 mutant cases. Our results indicate that RPA(-) cases demonstrate tumor suppressor deletions and genome instability, but lack unique or recurrent genetic lesions compensating for the lack of RPAs. Larger WGS studies of RPA(-) cases are required to understand this important LUAD subset.

The IkappaB kinases IKKalpha and IKKbeta are necessary and sufficient for skeletal muscle atrophy.

Nuclear factor-kappaB (NF-kappaB) signaling is necessary for many types of muscle atrophy, yet only some of the required components have been identified. Gene transfer of a dominant negative (d.n.) IKKbeta into rat soleus muscles showed complete inhibition of 7-day disuse-induced activation of a kappaB reporter gene, while overexpression of wild-type (w.t.) IKKbeta did not. Overexpression of a d.n. IKKbeta-EGFP fusion protein showed that atrophy was inhibited by 50%, indicating that IKKbeta is required for the atrophy process. Overexpression of constitutively active (c.a.) IKKbeta-EGFP showed a marked increase in NF-kappaB activity and a decrease in fiber size of weight-bearing soleus muscles, while muscles overexpressing w.t. IKKbeta-HA had no effect. The same results were found for IKKalpha; overexpression of a d.n. form of the protein decreased unloading-induced NF-kappaB activation and inhibited atrophy by 50%, while overexpression of the w.t. protein had no effect. Overexpression of a c.a. IKKalpha-EGFP fusion protein showed that IKKalpha was sufficient to activate NF-kappaB activity and induce fiber atrophy in muscle. Overexpression of d.n. IKKbeta plus d.n. IKKalpha showed an additive effect on the inhibition of disuse atrophy (70%), suggesting that both kinases of the IKK complex are required for muscle atrophy. These data show that both IKKalpha and IKKbeta are necessary and sufficient for physiological muscle atrophy.

Development and validation of a NanoString BASE47 bladder cancer gene classifier.

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.

Adaptation and selection shape clonal evolution of tumors during residual disease and recurrence.

The survival and recurrence of residual tumor cells following therapy constitutes one of the biggest obstacles to obtaining cures in breast cancer, but it remains unclear how the clonal composition of tumors changes during relapse. We use cellular barcoding to monitor clonal dynamics during tumor recurrence in vivo. We find that clonal diversity decreases during tumor regression, residual disease, and recurrence. The recurrence of dormant residual cells follows several distinct routes. Approximately half of the recurrent tumors exhibit clonal dominance with a small number of subclones comprising the vast majority of the tumor; these clonal recurrences are frequently dependent upon Met gene amplification. A second group of recurrent tumors comprises thousands of subclones, has a clonal architecture similar to primary tumors, and is dependent upon the Jak/Stat pathway. Thus the regrowth of dormant tumors proceeds via multiple routes, producing recurrent tumors with distinct clonal composition, genetic alterations, and drug sensitivities.

Comprehensive Analysis of Genetic Ancestry and Its Molecular Correlates in Cancer.

We evaluated ancestry effects on mutation rates, DNA methylation, and mRNA and miRNA expression among 10,678 patients across 33 cancer types from The Cancer Genome Atlas. We demonstrated that cancer subtypes and ancestry-related technical artifacts are important confounders that have been insufficiently accounted for. Once accounted for, ancestry-associated differences spanned all molecular features and hundreds of genes. Biologically significant differences were usually tissue specific but not specific to cancer. However, admixture and pathway analyses suggested some of these differences are causally related to cancer. Specific findings included increased FBXW7 mutations in patients of African origin, decreased VHL and PBRM1 mutations in renal cancer patients of African origin, and decreased immune activity in bladder cancer patients of East Asian origin.