Abstract

We used whole-organ mapping to study loco-geographic molecular changes in evolution of human bladder cancer from mucosal field effects. The integrative multi-platform analyses based on genome-wide RNA sequencing, methylation, copy number variations, and whole exome sequencing identified over 100 dysregulated canonical pathways involving immunity, tissue differentiation and transformation as initiators of bladder carcinogenesis. Widespread dysregulation of interleukin signaling was the dominant change signifying the important role of inflammation and immunity in the incipient phases of urothelial carcinogenesis. The analyses of mutational patterns identified three types of mutations based on their geographic distribution and variant allele frequencies. The most common were low frequency subclonal mutations restricted to individual mucosal samples which were the progeny of their respective uroprogenitor cells. The two additional types of mutations were associated with clonal expansion and involved large areas of bladder mucosa. The first group referred to as mutations, showed a low mutational frequency across the mucosa. The second group referred to as {beta} mutations increased in their frequencies with disease progression and a large proportion of them represented mutated transcriptional regulators controlling proliferation. Time modeling revealed that bladder carcinogenesis is spanning 10-15 years and can be divided into dormant and progressive phases. The progressive phase lasted 1-2 years and was primarily driven by {beta} mutations with high proliferative advantage. This is the first detailed molecular characterization of mucosal field effects initiating bladder carcinogenesis on the whole-organ scale. It provides novel insights into incipient phases of bladder carcinogenesis and biomarkers for early detection of bladder cancer as well as targets for preventive therapies.