Molecular Classification and Tumor Microenvironment Characterization of Gallbladder Cancer by Comprehensive Genomic and Transcriptomic Analysis

SIMPLE SUMMARY: Gallbladder cancer (GBC) is a rare but lethal cancer. Molecular characterization of GBC is insufficient so far, and a comprehensive molecular portrait is warranted to uncover new targets and classify GBC. Clustering analysis of RNA expression revealed two subclasses of 36 GBCs, which reflects the status of the tumor microenvironment (TME) and poor prognosis of GBC, including epithelial–mesenchymal transition (EMT), immune suppression, and the TGF-β signaling pathway. The knockout of miR125B1 in GBC cell lines decreased its invasion ability and altered the EMT pathway. Mutations of the genes related to the TGF-β signaling pathway were enriched in the poor-prognosis/TME-rich cluster of GBCs. This comprehensive molecular analysis provides a new classification of GBCs based on the TME activity, which is involved with EMT and immune suppression for poor prognosis of GBC. This information may be useful for GBC prognosis and therapeutic decision-making. ABSTRACT: Gallbladder cancer (GBC), a rare but lethal disease, is often diagnosed at advanced stages. So far, molecular characterization of GBC is insufficient, and a comprehensive molecular portrait is warranted to uncover new targets and classify GBC. We performed a transcriptome analysis of both coding and non-coding RNAs from 36 GBC fresh-frozen samples. The results were integrated with those of comprehensive mutation profiling based on whole-genome or exome sequencing. The clustering analysis of RNA-seq data facilitated the classification of GBCs into two subclasses, characterized by high or low expression levels of TME (tumor microenvironment) genes. A correlation was observed between gene expression and pathological immunostaining. TME-rich tumors showed significantly poor prognosis and higher recurrence rate than TME-poor tumors. TME-rich tumors showed overexpression of genes involved in epithelial-to-mesenchymal transition (EMT) and inflammation or immune suppression, which was validated by immunostaining. One non-coding RNA, miR125B1, exhibited elevated expression in stroma-rich tumors, and miR125B1 knockout in GBC cell lines decreased its invasion ability and altered the EMT pathway. Mutation profiles revealed TP53 (47%) as the most commonly mutated gene, followed by ELF3 (13%) and ARID1A (11%). Mutations of ARID1A, ERBB3, and the genes related to the TGF-β signaling pathway were enriched in TME-rich tumors. This comprehensive analysis demonstrated that TME, EMT, and TGF-β pathway alterations are the main drivers of GBC and provides a new classification of GBCs that may be useful for therapeutic decision-making.