Martin Lauss

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.

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.

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.

Gene expression profiling of primary male breast cancers reveals two unique subgroups and identifies N-acetyltransferase-1 (NAT1) as a novel prognostic biomarker

IntroductionMale breast cancer (MBC) is a rare and inadequately characterized disease. The aim of the present study was to characterize MBC tumors transcriptionally, to classify them into comprehensive subgroups, and to compare them with female breast cancer (FBC).MethodsA total of 66 clinicopathologically well-annotated fresh frozen MBC tumors were analyzed using Illumina Human HT-12 bead arrays, and a tissue microarray with 220 MBC tumors was constructed for validation using immunohistochemistry. Two external gene expression datasets were used for comparison purposes: 37 MBCs and 359 FBCs.ResultsUsing an unsupervised approach, we classified the MBC tumors into two subgroups, luminal M1 and luminal M2, respectively, with differences in tumor biological features and outcome, and which differed from the intrinsic subgroups described in FBC. The two subgroups were recapitulated in the external MBC dataset. Luminal M2 tumors were characterized by high expression of immune response genes and genes associated with estrogen receptor (ER) signaling. Luminal M1 tumors, on the other hand, despite being ER positive by immunohistochemistry showed a lower correlation to genes associated with ER signaling and displayed a more aggressive phenotype and worse prognosis. Validation of two of the most differentially expressed genes, class 1 human leukocyte antigen (HLA) and the metabolizing gene N-acetyltransferase-1 (NAT1), respectively, revealed significantly better survival associated with high expression of both markers (HLA, hazard ratio (HR) 3.6, P = 0.002; NAT1, HR 2.5, P = 0.033). Importantly, NAT1 remained significant in a multivariate analysis (HR 2.8, P = 0.040) and may thus be a novel prognostic marker in MBC.ConclusionsWe have detected two unique and stable subgroups of MBC with differences in tumor biological features and outcome. They differ from the widely acknowledged intrinsic subgroups of FBC. As such, they may constitute two novel subgroups of breast cancer, occurring exclusively in men, and which may consequently require novel treatment approaches. Finally, we identified NAT1 as a possible prognostic biomarker for MBC, as suggested by NAT1 positivity corresponding to better outcome.

Molecular Profiling Reveals Low- and High-Grade Forms of Primary Melanoma

Purpose: For primary melanomas, tumor thickness, mitotic rate, and ulceration are well-laid cornerstones of prognostication. However, a molecular exposition of melanoma aggressiveness is critically missing. We recently uncovered a four-class structure in metastatic melanoma, which predicts outcome and informs biology. This raises the possibility that a molecular structure exists even in the early stages of melanoma and that molecular determinants could underlie histophenotype and eventual patient outcome. Experimental Design: We subjected 223 archival primary melanomas to a horizontally integrated analysis of RNA expression, oncogenic mutations at 238 lesions, histomorphometry, and survival data. Results: Our previously described four-class structure that was elucidated in metastatic lesions was evident within the expression space of primary melanomas. Because these subclasses converged into two larger prognostic and phenotypic groups, we used the metastatic lesions to develop a binary subtype-based signature capable of distinguishing between “high” and “low” grade forms of the disease. The two-grade signature was subsequently applied to the primary melanomas. Compared with low-grade tumors, high-grade primary melanomas were significantly associated with increased tumor thickness, mitotic rate, ulceration (all P < 0.01), and poorer relapse-free (HR = 4.94; 95% CI, 2.84–8.59), and overall (HR = 3.66; 95% CI, 2.40–5.58) survival. High-grade melanomas exhibited elevated levels of proliferation and BRCA1/DNA damage signaling genes, whereas low-grade lesions harbored higher expression of immune genes. Importantly, the molecular-grade signature was validated in two external gene expression data sets. Conclusions: We provide evidence for a molecular organization within melanomas, which is preserved across all stages of disease. Clin Cancer Res; 18(15); 4026–36. ©2012 AACR.

Efficient Arsenic Metabolism - The AS3MT Haplotype Is Associated with DNA Methylation and Expression of Multiple Genes Around AS3MT.

Arsenic is a very potent toxicant. One major susceptibility factor for arsenic-related toxicity is the efficiency of arsenic metabolism. The efficiency, in turn, is associated with non-coding single nucleotide polymorphisms (SNPs) in the arsenic methyltransferase AS3MT on chromosome 10q24. However, the mechanism of action for these SNPs is not yet clarified. Here, we assessed the influence of genetic variation in AS3MT on DNA methylation and gene expression within 10q24, in people exposed to arsenic in drinking water. DNA was extracted from peripheral blood from women in the Argentinean Andes (N = 103) and from cord blood from new-borns in Bangladesh (N = 127). AS3MT SNPs were analyzed with Sequenom or Taqman assays. Whole genome epigenetic analysis with Infinium HumanMethylation450 BeadChip was performed on bisulphite-treated DNA. Whole genome gene expression analysis was performed with Illumina DirectHyb HumanHT-12 v4.0 on RNA from peripheral blood. Arsenic exposure was assessed by HPLC-ICPMS. In the Argentinean women, the major AS3MT haplotype, associated with more efficient arsenic metabolism, showed increased methylation of AS3MT (p = 10−6) and also differential methylation of several other genes within about 800 kilobasepairs: CNNM2 (p<10−16), NT5C2 (p<10−16), C10orf26 (p = 10−8), USMG5 (p = 10−5), TRIM8 (p = 10−4), and CALHM2 (p = 0.038) (adjusted for multiple comparisons). Similar, but weaker, associations between AS3MT haplotype and DNA methylation in 10q24 were observed in cord blood (Bangladesh). The haplotype-associated altered CpG methylation was correlated with reduced expression of AS3MT and CNNM2 (rs = −0.22 to −0.54), and with increased expression of NT5C2 and USMG5 (rs = 0.25 to 0.58). Taking other possibly influential variables into account in multivariable linear models did only to a minor extent alter the strength of the associations. In conclusion, the AS3MT haplotype status strongly predicted DNA methylation and gene expression of AS3MT as well as several genes in 10q24. This raises the possibility that several genes in this region are important for arsenic metabolism.

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.

Prediction of Stage, Grade, and Survival in Bladder Cancer Using Genome-wide Expression Data: A Validation Study

Purpose: To evaluate performances of published gene signatures for the assessment of urothelial carcinoma. Experimental Design: We evaluated 28 published gene signatures designed for diagnostic and prognostic purposes of urothelial cancer. The investigated signatures include eight signatures for stage, five for grade, four for progression, and six for survival. We used two algorithms for classification, nearest centroid classification and support vector machine, and Cox regression to evaluate signature performance in four independent data sets. Results: The overlap of genes among the signatures was low, ranging from 11% among stage signatures to 0.6% among survival signatures. The published signatures predicted muscle-invasive and high-grade tumors with accuracies in the range of 70% to 90%. The performance for a given signature varied considerably with the validation data set used, and interestingly, some of the best performing signatures were not designed for the tested classification problem. In addition, several nonbladder-derived gene signatures performed equally well. Large randomly selected gene signatures performed better than the published signatures, and by systematically increasing signature size, we show that signatures with >150 genes are needed to obtain robust performance in independent validation data sets. None of the published survival signatures performed better than random assignments when applied to independent validation data. Conclusion: We conclude that gene expression signatures with >150 genes predict muscle-invasive growth and high-grade tumors with robust accuracies. Special considerations have to be taken when designing gene signatures for outcome in bladder cancer. Clin Cancer Res; 16(17); 4421–33. ©2010 AACR.