Gottfrid Sjödahl

A Molecular Taxonomy for Urothelial Carcinoma

Purpose: Even though urothelial cancer is the fourth most common tumor type among males, progress in treatment has been scarce. A problem in day-to-day clinical practice is that precise assessment of individual tumors is still fairly uncertain; consequently efforts have been undertaken to complement tumor evaluation with molecular biomarkers. An extension of this approach would be to base tumor classification primarily on molecular features. Here, we present a molecular taxonomy for urothelial carcinoma based on integrated genomics. Experimental Design: We use gene expression profiles from 308 tumor cases to define five major urothelial carcinoma subtypes: urobasal A, genomically unstable, urobasal B, squamous cell carcinoma like, and an infiltrated class of tumors. Tumor subtypes were validated in three independent publically available data sets. The expression of 11 key genes was validated at the protein level by immunohistochemistry. Results: The subtypes show distinct clinical outcomes and differ with respect to expression of cell-cycle genes, receptor tyrosine kinases particularly FGFR3, ERBB2, and EGFR, cytokeratins, and cell adhesion genes, as well as with respect to FGFR3, PIK3CA, and TP53 mutation frequency. The molecular subtypes cut across pathologic classification, and class-defining gene signatures show coordinated expression irrespective of pathologic stage and grade, suggesting the molecular phenotypes as intrinsic properties of the tumors. Available data indicate that susceptibility to specific drugs is more likely to be associated with the molecular stratification than with pathologic classification. Conclusions: We anticipate that the molecular taxonomy will be useful in future clinical investigations. Clin Cancer Res; 18(12); 3377–86. ©2012 AACR.

Combined Gene Expression and Genomic Profiling Define Two Intrinsic Molecular Subtypes of Urothelial Carcinoma and Gene Signatures for Molecular Grading and Outcome

In the present investigation, we sought to refine the classification of urothelial carcinoma by combining information on gene expression, genomic, and gene mutation levels. For these purposes, we performed gene expression analysis of 144 carcinomas, and whole genome array-CGH analysis and mutation analyses of FGFR3, PIK3CA, KRAS, HRAS, NRAS, TP53, CDKN2A, and TSC1 in 103 of these cases. Hierarchical cluster analysis identified two intrinsic molecular subtypes, MS1 and MS2, which were validated and defined by the same set of genes in three independent bladder cancer data sets. The two subtypes differed with respect to gene expression and mutation profiles, as well as with the level of genomic instability. The data show that genomic instability was the most distinguishing genomic feature of MS2 tumors, and that this trait was not dependent on TP53/MDM2 alterations. By combining molecular and pathologic data, it was possible to distinguish two molecular subtypes of T(a) and T(1) tumors, respectively. In addition, we define gene signatures validated in two independent data sets that classify urothelial carcinoma into low-grade (G(1)/G(2)) and high-grade (G(3)) tumors as well as non-muscle and muscle-invasive tumors with high precisions and sensitivities, suggesting molecular grading as a relevant complement to standard pathologic grading. We also present a gene expression signature with independent prognostic effect on metastasis and disease-specific survival. We conclude that the combination of molecular and histopathologic classification systems might provide a strong improvement for bladder cancer classification and produce new insights into the development of this tumor type.

Integrated genomic and gene expression profiling identifies two major genomic circuits in urothelial carcinoma.

Similar to other malignancies, urothelial carcinoma (UC) is characterized by specific recurrent chromosomal aberrations and gene mutations. However, the interconnection between specific genomic alterations, and how patterns of chromosomal alterations adhere to different molecular subgroups of UC, is less clear. We applied tiling resolution array CGH to 146 cases of UC and identified a number of regions harboring recurrent focal genomic amplifications and deletions. Several potential oncogenes were included in the amplified regions, including known oncogenes like E2F3, CCND1, and CCNE1, as well as new candidate genes, such as SETDB1 (1q21), and BCL2L1 (20q11). We next combined genome profiling with global gene expression, gene mutation, and protein expression data and identified two major genomic circuits operating in urothelial carcinoma. The first circuit was characterized by FGFR3 alterations, overexpression of CCND1, and 9q and CDKN2A deletions. The second circuit was defined by E3F3 amplifications and RB1 deletions, as well as gains of 5p, deletions at PTEN and 2q36, 16q, 20q, and elevated CDKN2A levels. TP53/MDM2 alterations were common for advanced tumors within the two circuits. Our data also suggest a possible RAS/RAF circuit. The tumors with worst prognosis showed a gene expression profile that indicated a keratinized phenotype. Taken together, our integrative approach revealed at least two separate networks of genomic alterations linked to the molecular diversity seen in UC, and that these circuits may reflect distinct pathways of tumor development.

Impact of Molecular Subtypes in Muscle-invasive Bladder Cancer on Predicting Response and Survival after Neoadjuvant Chemotherapy

BACKGROUND An early report on the molecular subtyping of muscle-invasive bladder cancer (MIBC) by gene expression suggested that response to neoadjuvant chemotherapy (NAC) varies by subtype. OBJECTIVE To investigate the ability of molecular subtypes to predict pathological downstaging and survival after NAC. DESIGN, SETTING, AND PARTICIPANTS Whole transcriptome profiling was performed on pre-NAC transurethral resection specimens from 343 patients with MIBC. Samples were classified according to four published molecular subtyping methods. We developed a single-sample genomic subtyping classifier (GSC) to predict consensus subtypes (claudin-low, basal, luminal-infiltrated and luminal) with highest clinical impact in the context of NAC. Overall survival (OS) according to subtype was analyzed and compared with OS in 476 non-NAC cases (published datasets). INTERVENTION Gene expression analysis was used to assign subtypes. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS Receiver-operating characteristics were used to determine the accuracy of GSC. The effect of GSC on survival was estimated by Cox proportional hazard regression models. RESULTS AND LIMITATIONS The models generated subtype calls in expected ratios with high concordance across subtyping methods. GSC was able to predict four consensus molecular subtypes with high accuracy (73%), and clinical significance of the predicted consensus subtypes could be validated in independent NAC and non-NAC datasets. Luminal tumors had the best OS with and without NAC. Claudin-low tumors were associated with poor OS irrespective of treatment regimen. Basal tumors showed the most improvement in OS with NAC compared with surgery alone. The main limitations of our study are its retrospective design and comparison across datasets. CONCLUSIONS Molecular subtyping may have an impact on patient benefit to NAC. If validated in additional studies, our results suggest that patients with basal tumors should be prioritized for NAC. We discovered the first single-sample classifier to subtype MIBC, which may be suitable for integration into routine clinical practice. PATIENT SUMMARY Different molecular subtypes can be identified in muscle-invasive bladder cancer. Although cisplatin-based neoadjuvant chemotherapy improves patient outcomes, we identified that the benefit is highest in patients with basal tumors. Our newly discovered classifier can identify these molecular subtypes in a single patient and could be integrated into routine clinical practice after further validation.

A Systematic Study of Gene Mutations in Urothelial Carcinoma; Inactivating Mutations in TSC2 and PIK3R1.

Background Urothelial carcinoma (UC) is characterized by frequent gene mutations of which activating mutations in FGFR3 are the most frequent. Several downstream targets of FGFR3 are also mutated in UC, e.g., PIK3CA, AKT1, and RAS. Most mutation studies of UCs have been focused on single or a few genes at the time or been performed on small sample series. This has limited the possibility to investigate co-occurrence of mutations. Methodology/Principal Findings We performed mutation analyses of 16 genes, FGFR3, PIK3CA, PIK3R1 PTEN, AKT1, KRAS, HRAS, NRAS, BRAF, ARAF, RAF1, TSC1, TSC2, APC, CTNNB1, and TP53, in 145 cases of UC. We show that FGFR3 and PIK3CA mutations are positively associated. In addition, we identified PIK3R1 as a target for mutations. We demonstrate a negative association at borderline significance between FGFR3 and RAS mutations, and show that these mutations are not strictly mutually exclusive. We show that mutations in BRAF, ARAF, RAF1 rarely occurs in UC. Our data emphasize the possible importance of APC signaling as 6% of the investigated tumors either showed inactivating APC or activating CTNNB1 mutations. TSC1, as well as TSC2, that constitute the mTOR regulatory tuberous sclerosis complex were found to be mutated at a combined frequency of 15%. Conclusions/Significance Our data demonstrate a significant association between FGFR3 and PIK3CA mutations in UC. Moreover, the identification of mutations in PIK3R1 further emphasizes the importance of the PI3-kinase pathway in UC. The presence of TSC2 mutations, in addition to TSC1 mutations, underlines the involvement of mTOR signaling in UC.

Toward a Molecular Pathologic Classification of Urothelial Carcinoma

We recently defined molecular subtypes of urothelial carcinomas according to whole genome gene expression. Herein we describe molecular pathologic characterization of the subtypes using 20 genes and IHC of 237 tumors. In addition to differences in expression levels, the subtypes show important differences in stratification of protein expression. The selected genes included biological features central to bladder cancer biology, eg, cell cycle activity, cellular architecture, cell-cell interactions, and key receptor tyrosine kinases. We show that the urobasal (Uro) A subtype shares features with normal urothelium such as keratin 5 (KRT5), P-cadherin (P-Cad), and epidermal growth factor receptor (EGFR) expression confined to basal cells, and cell cycle activity (CCNB1) restricted to the tumor-stroma interface. In contrast, the squamous cell cancer-like (SCCL) subtype uniformly expresses KRT5, P-Cad, EGFR, KRT14, and cell cycle genes throughout the tumor parenchyma. The genomically unstable subtype shows proliferation throughout the tumor parenchyma and high ERBB2 and E-Cad expression but absence of KRT5, P-Cad, and EGFR expression. UroB tumors demonstrate features shared by both UroA and SCCL subtypes. A major transition in tumor progression seems to be loss of dependency of stromal interaction for proliferation. We present a simple IHC/histology-based classifier that is easy to implement as a standard pathologic evaluation to differentiate the three major subtypes: urobasal, genomically unstable, and SCCL. These three major subtypes exhibit important prognostic differences.

A Molecular Pathologic Framework for Risk Stratification of Stage T1 Urothelial Carcinoma

BACKGROUND One third of patients with stage T1 urothelial carcinoma (UC) progress to muscle-invasive disease requiring radical surgery. Thus, reliable tools are needed for risk stratification of stage T1 UC. OBJECTIVE To investigate the extent to which stratification of stage T1 tumours into previously described molecular pathologic UC subtypes can provide improved information on tumour progression. DESIGN, SETTING, AND PARTICIPANTS A population-based cohort of 167 primary stage T1 UCs was characterised by immunohistochemistry and classified into the molecular subtypes urobasal (Uro, 32%), genomically unstable (GU, 58%), and squamous-cell-carcinoma-like (SCCL, 10%). OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS Progression-free survival using univariate and multivariate models. RESULTS AND LIMITATIONS Subtype classification was validated using nine additional markers with known subtype-specific expression. Analysis of mRNA expression of progression biomarkers revealed a strong association with molecular subtype. Kaplan-Meier analyses showed that the risk of progression was low for Uro tumours and high for GU/SCCL tumours. High progression risk scores were found only for GU/SCCL tumours. Clinical risk factors such as multifocality, concomitant carcinoma in situ, invasion depth, lymphovascular invasion, and high CD3(+) lymphocyte infiltration were observed almost exclusively in GU/SCCL cases. CONCLUSIONS Molecular subtypes Uro, GU, and SCCL were identified in an independent population-based cohort of stage T1 UCs. Biomarkers and clinical risk factors for progression were associated with molecular subtype. Rapidly progressing T1 tumours were of subtype GU or SCCL and had either a high progression risk score or an elevated CD3(+) cell count. PATIENT SUMMARY We show that classification of stage T1 urothelial carcinoma into molecular subtypes can improve the identification of patients with progressing tumours.

DNA methylation analyses of urothelial carcinoma reveal distinct epigenetic subtypes and an association between gene copy number and methylation status

We assessed DNA methylation and copy number status of 27,000 CpGs in 149 urothelial carcinomas and integrated the findings with gene expression and mutation data. Methylation was associated with gene expression for 1,332 CpGs, of which 26% showed positive correlation with expression, i.e., high methylation and high gene expression levels. These positively correlated CpGs were part of specific transcription factor binding sites, such as sites for MYC and CREBP1, or located in gene bodies. Furthermore, we found genes with copy number gains, low expression and high methylation levels, revealing an association between methylation and copy number levels. This phenomenon was typically observed for developmental genes, such as HOX genes, and tumor suppressor genes. In contrast, we also identified genes with copy number gains, high expression and low methylation levels. This was for instance observed for some keratin genes. Tumor cases could be grouped into four subgroups, termed epitypes, by their DNA methylation profiles. One epitype was influenced by the presence of infiltrating immune cells, two epitypes were mainly composed of non-muscle invasive tumors, and the remaining epitype of muscle invasive tumors. The polycomb complex protein EZH2 that blocks differentiation in embryonic stem cells showed increased expression both at the mRNA and protein levels in the muscle invasive epitype, together with methylation of polycomb target genes and HOX genes. Our data highlights HOX gene silencing and EZH2 expression as mechanisms to promote a more undifferentiated and aggressive state in UC.

Molecular classification of urothelial carcinoma: global mRNA classification versus tumour-cell phenotype classification

Global mRNA expression analysis is efficient for phenotypic profiling of tumours, and has been used to define molecular subtypes for almost every major tumour type. A key limitation is that most tumours are communities of both tumour and non‐tumour cells. This problem is particularly pertinent for analysis of advanced invasive tumours, which are known to induce major changes and responses in both the tumour and the surrounding tissue. To identify bladder cancer tumour‐cell phenotypes and compare classification by tumour‐cell phenotype with classification by global gene expression analysis, we analysed 307 advanced bladder cancers (cystectomized) both by genome gene expression analysis and by immunohistochemistry with antibodies for 28 proteins. According to systematic analysis of gene and protein expression data, focusing on key molecular processes, we describe five tumour‐cell phenotypes of advanced urothelial carcinoma: urothelial‐like, genomically unstable, basal/SCC‐like, mesenchymal‐like, and small‐cell/neuroendocrine‐like. We provide molecular pathological definitions for each subtype. Tumours expressing urothelial differentiation factors show inconsistent and abnormal protein expression of terminal differentiation markers, suggesting pseudo‐differentiation. Cancers with different tumour‐cell phenotypes may co‐cluster (converge), and cases with identical tumour‐cell phenotypes may cluster apart (diverge), in global mRNA analyses. This divergence/convergence suggests that broad global commonalities related to the invasive process may exist between muscle‐invasive tumours regardless of specific tumour‐cell phenotype. Hence, there is a systematic disagreement in subtype classification determined by global mRNA profiling and by immunohistochemical profiling at the tumour‐cell level. We suggest that a combination of molecular pathology (tumour‐cell phenotype) and global mRNA profiling (context) is required for adequate subtype classification of muscle‐invasive bladder cancer.

Biological determinants of bladder cancer gene expression subtypes

Molecular stratification of tumors by gene expression profiling has been applied to a large number of human malignancies and holds great promise for personalized treatment. Comprehensive classification schemes for urothelial carcinoma have been proposed by three separate groups but have not previously been evaluated simultaneously in independent data. Here we map the interrelations between the proposed molecular subtypes onto the intrinsic structure of a rich independent dataset and show that subtype stratification within each scheme can be explained in terms of a set of common underlying biological processes. We highlight novel biological and genomic drivers of urothelial carcinoma molecular subtypes and show that tumors carrying genomic aberrations characteristic of distinct molecular pathways converge on a common top level phenotype corresponding to the two major molecular subtypes of non-muscle invasive disease.