Network modeling links kidney developmental programs and the cancer type-specificity of VHL mutations.

Elucidating the molecular dependencies behind the cancer-type specificity of driver mutations may reveal new therapeutic opportunities. We hypothesized that developmental programs would impact the transduction of oncogenic signaling activated by a driver mutation and shape its cancer-type specificity. Therefore, we designed a computational analysis framework by combining single-cell gene expression profiles during fetal organ development, latent factor discovery, and information theory-based differential network analysis to systematically identify transcription factors that selectively respond to driver mutations under the influence of organ-specific developmental programs. After applying this approach to VHL mutations, which are highly specific to clear cell renal cell carcinoma (ccRCC), we revealed important regulators downstream of VHL mutations in ccRCC and used their activities to cluster patients with ccRCC into three subtypes. This classification revealed a more significant difference in prognosis than the previous mRNA profile-based method and was validated in an independent cohort. Moreover, we found that EP300, a key epigenetic factor maintaining the regulatory network of the subtype with the worst prognosis, can be targeted by a small inhibitor, suggesting a potential treatment option for a subset of patients with ccRCC. This work demonstrated an intimate relationship between organ development and oncogenesis from the perspective of systems biology, and the method can be generalized to study the influence of other biological processes on cancer driver mutations.