Meta-Analysis of the Luminal and Basal Subtypes of Bladder Cancer and the Identification of Signature Immunohistochemical Markers for Clinical Use

AbstractBackgroundIt has been suggested that bladder cancer can be divided into two molecular subtypes referred to as luminal and basal with distinct clinical behaviors and sensitivities to chemotherapy. We aimed to validate these subtypes in several clinical cohorts and identify signature immunohistochemical markers that would permit simple and cost-effective classification of the disease in primary care centers.MethodsWe analyzed genomic expression profiles of bladder cancer in three cohorts of fresh frozen tumor samples: MD Anderson (n=132), Lund (n=308), and The Cancer Genome Atlas (TCGA) (n=408) to validate the expression signatures of luminal and basal subtypes and relate them to clinical follow-up data. We also used an MD Anderson cohort of archival bladder tumor samples (n=89) and a parallel tissue microarray to identify immunohistochemical markers that permitted the molecular classification of bladder cancer.FindingsBladder cancers could be assigned to two candidate intrinsic molecular subtypes referred to here as luminal and basal in all of the datasets analyzed. Luminal tumors were characterized by the expression signature similar to the intermediate/superficial layers of normal urothelium. They showed the upregulation of PPARγ target genes and the enrichment for FGFR3, ELF3, CDKN1A, and TSC1 mutations. In addition, luminal tumors were characterized by the overexpression of E-Cadherin, HER2/3, Rab-25, and Src. Basal tumors showed the expression signature similar to the basal layer of normal urothelium. They showed the upregulation of p63 target genes, the enrichment for TP53 and RB1 mutations, and overexpression of CD49, Cyclin B1, and EGFR. Survival analyses showed that the muscle-invasive basal bladder cancers were more aggressive when compared to luminal cancers. The immunohistochemical expressions of only two markers, luminal (GATA3) and basal (KRT5/6), were sufficient to identify the molecular subtypes of bladder cancer with over 90% accuracy.InterpretationThe molecular subtypes of bladder cancer have distinct clinical behaviors and sensitivities to chemotherapy, and a simple two-marker immunohistochemical classifier can be used for prognostic and therapeutic stratification.FundingU.S. National Cancer Institute and National Institute of Health

Pan-Cancer Analysis of lncRNA Regulation Supports Their Targeting of Cancer Genes in Each Tumor Context

Long noncoding RNAs (lncRNAs) are commonly dys-regulated in tumors, but only a handful are known toplay pathophysiological roles in cancer. We inferredlncRNAs that dysregulate cancer pathways, onco-genes, and tumor suppressors (cancer genes) bymodeling their effects on the activity of transcriptionfactors, RNA-binding proteins, and microRNAs in5,185 TCGA tumors and 1,019 ENCODE assays.Our predictions included hundreds of candidateonco- and tumor-suppressor lncRNAs (cancerlncRNAs) whose somatic alterations account for thedysregulation of dozens of cancer genes and path-ways in each of 14 tumor contexts. To demonstrateproof of concept, we showed that perturbations tar-geting OIP5-AS1 (an inferred tumor suppressor) andTUG1 and WT1-AS (inferred onco-lncRNAs) dysre-gulated cancer genes and altered proliferation ofbreast and gynecologic cancer cells. Our analysis in-dicates that, although most lncRNAs are dysregu-lated in a tumor-specific manner, some, includingOIP5-AS1, TUG1, NEAT1, MEG3, and TSIX, synergis-tically dysregulate cancer pathways in multiple tumorcontexts

Pan-cancer Alterations of the MYC Oncogene and Its Proximal Network across the Cancer Genome Atlas

Although theMYConcogene has been implicated incancer, a systematic assessment of alterations ofMYC, related transcription factors, and co-regulatoryproteins, forming the proximal MYC network (PMN),across human cancers is lacking. Using computa-tional approaches, we define genomic and proteo-mic features associated with MYC and the PMNacross the 33 cancers of The Cancer Genome Atlas.Pan-cancer, 28% of all samples had at least one ofthe MYC paralogs amplified. In contrast, the MYCantagonists MGA and MNT were the most frequentlymutated or deleted members, proposing a roleas tumor suppressors.MYCalterations were mutu-ally exclusive withPIK3CA,PTEN,APC,orBRAFalterations, suggesting that MYC is a distinct onco-genic driver. Expression analysis revealed MYC-associated pathways in tumor subtypes, such asimmune response and growth factor signaling; chro-matin, translation, and DNA replication/repair wereconserved pan-cancer. This analysis reveals insightsinto MYC biology and is a reference for biomarkersand therapeutics for cancers with alterations ofMYC or the PMN

Genomic, Pathway Network, and Immunologic Features Distinguishing Squamous Carcinomas

This integrated, multiplatform PanCancer Atlas study co-mapped and identified distinguishing molecular features of squamous cell carcinomas (SCCs) from five sites associated with smokin

Spatial Organization and Molecular Correlation of Tumor-Infiltrating Lymphocytes Using Deep Learning on Pathology Images

Beyond sample curation and basic pathologic characterization, the digitized H&E-stained images of TCGA samples remain underutilized. To highlight this resource, we present mappings of tumorinfiltrating lymphocytes (TILs) based on H&E images from 13 TCGA tumor types. These TIL maps are derived through computational staining using a convolutional neural network trained to classify patches of images. Affinity propagation revealed local spatial structure in TIL patterns and correlation with overall survival. TIL map structural patterns were grouped using standard histopathological parameters. These patterns are enriched in particular T cell subpopulations derived from molecular measures. TIL densities and spatial structure were differentially enriched among tumor types, immune subtypes, and tumor molecular subtypes, implying that spatial infiltrate state could reflect particular tumor cell aberration states. Obtaining spatial lymphocytic patterns linked to the rich genomic characterization of TCGA samples demonstrates one use for the TCGA image archives with insights into the tumor-immune microenvironment

Integrated Genomic Analysis of the Ubiquitin Pathway across Cancer Types

Protein ubiquitination is a dynamic and reversibleprocess of adding single ubiquitin molecules orvarious ubiquitin chains to target proteins. Here,using multidimensional omic data of 9,125 tumorsamples across 33 cancer types from The CancerGenome Atlas, we perform comprehensive molecu-lar characterization of 929 ubiquitin-related genesand 95 deubiquitinase genes. Among them, we sys-tematically identify top somatic driver candidates,including mutatedFBXW7with cancer-type-specificpatterns and amplifiedMDM2showing a mutuallyexclusive pattern withBRAFmutations. Ubiquitinpathway genes tend to be upregulated in cancermediated by diverse mechanisms. By integratingpan-cancer multiomic data, we identify a group oftumor samples that exhibit worse prognosis. Thesesamples are consistently associated with the upre-gulation of cell-cycle and DNA repair pathways, char-acterized by mutatedTP53,MYC/TERTamplifica-tion, andAPC/PTENdeletion. Our analysishighlights the importance of the ubiquitin pathwayin cancer development and lays a foundation fordeveloping relevant therapeutic strategies

Machine Learning Identifies Stemness Features Associated with Oncogenic Dedifferentiation.

Cancer progression involves the gradual loss of a differentiated phenotype and acquisition of progenitor and stem-cell-like features. Here, we provide novel stemness indices for assessing the degree of oncogenic dedifferentiation. We used an innovative one-class logistic regression (OCLR) machine-learning algorithm to extract transcriptomic and epigenetic feature sets derived from non-transformed pluripotent stem cells and their differentiated progeny. Using OCLR, we were able to identify previously undiscovered biological mechanisms associated with the dedifferentiated oncogenic state. Analyses of the tumor microenvironment revealed unanticipated correlation of cancer stemness with immune checkpoint expression and infiltrating immune cells. We found that the dedifferentiated oncogenic phenotype was generally most prominent in metastatic tumors. Application of our stemness indices to single-cell data revealed patterns of intra-tumor molecular heterogeneity. Finally, the indices allowed for the identification of novel targets and possible targeted therapies aimed at tumor differentiation

Molecular characterization and clinical relevance of metabolic expression subtypes in human cancers.

Metabolic reprogramming provides critical information for clinical oncology. Using molecular data of 9,125 patient samples from The Cancer Genome Atlas, we identified tumor subtypes in 33 cancer types based on mRNA expression patterns of seven major metabolic processes and assessed their clinical relevance. Our metabolic expression subtypes correlated extensively with clinical outcome: subtypes with upregulated carbohydrate, nucleotide, and vitamin/cofactor metabolism most consistently correlated with worse prognosis, whereas subtypes with upregulated lipid metabolism showed the opposite. Metabolic subtypes correlated with diverse somatic drivers but exhibited effects convergent on cancer hallmark pathways and were modulated by highly recurrent master regulators across cancer types. As a proof-of-concept example, we demonstrated that knockdown of SNAI1 or RUNX1—master regulators of carbohydrate metabolic subtypes-modulates metabolic activity and drug sensitivity. Our study provides a system-level view of metabolic heterogeneity within and across cancer types and identifies pathway cross-talk, suggesting related prognostic, therapeutic, and predictive utility

Detection of bladder cancer using proteomic profiling of urine sediments.

We used protein expression profiles to develop a classification rule for the detection and prognostic assessment of bladder cancer in voided urine samples. Using the Ciphergen PBS II ProteinChip Reader, we analyzed the protein profiles of 18 pairs of samples of bladder tumor and adjacent urothelium tissue, a training set of 85 voided urine samples (32 controls and 53 bladder cancer), and a blinded testing set of 68 voided urine samples (33 controls and 35 bladder cancer). Using t-tests, we identified 473 peaks showing significant differential expression across different categories of paired bladder tumor and adjacent urothelial samples compared to normal urothelium. Then the intensities of those 473 peaks were examined in a training set of voided urine samples. Using this approach, we identified 41 protein peaks that were differentially expressed in both sets of samples. The expression pattern of the 41 protein peaks was used to classify the voided urine samples as malignant or benign. This approach yielded a sensitivity and specificity of 59% and 90%, respectively, on the training set and 80% and 100%, respectively, on the testing set. The proteomic classification rule performed with similar accuracy in low- and high-grade bladder carcinomas. In addition, we used hierarchical clustering with all 473 protein peaks on 65 benign voided urine samples, 88 samples from patients with clinically evident bladder cancer, and 127 samples from patients with a history of bladder cancer to classify the samples into Cluster A or B. The tumors in Cluster B were characterized by clinically aggressive behavior with significantly shorter metastasis-free and disease-specific survival