Individualized discovery of rare cancer drivers in global network context

  1. Iurii Petrov
  2. Andrey Alexeyenko  Is a corresponding author
  1. Karolinska Institutet, Sweden


Late advances in genome sequencing expanded the space of known cancer driver genes several-fold. However, most of this surge was based on computational analysis of somatic mutation frequencies and/or their impact on the protein function. On the contrary, experimental research necessarily accounted for functional context of mutations interacting with other genes and conferring cancer phenotypes. Eventually, just such results become 'hard currency' of cancer biology. The new method, NEAdriver employs knowledge accumulated thus far in the form of global interaction network and functionally annotated pathways in order to recover known and predict novel driver genes. The driver discovery was individualized by accounting for mutations' co-occurrence in each tumour genome - as an alternative to summarizing information over the whole cancer patient cohorts. For each somatic genome change, probabilistic estimates from two lanes of network analysis were combined into joint likelihoods of being a driver. Thus, ability to detect previously unnoticed candidate driver events emerged from combining individual genomic context with network perspective. The procedure was applied to ten largest cancer cohorts followed by evaluating error rates against previous cancer gene sets. The discovered driver combinations were shown to be informative on cancer outcome. This revealed driver genes with individually sparse mutation patterns that would not be detectable by other computational methods and related to cancer biology domains poorly covered by previous analyses. In particular, recurrent mutations of collagen, laminin, and integrin genes were observed in the adenocarcinoma and glioblastoma cancers. Considering constellation patterns of candidate drivers in individual cancer genomes opens a novel avenue for personalized cancer medicine.

Data availability

All data generated or analysed during this study are included in the manuscript and supporting files.

The following previously published data sets were used

Article and author information

Author details

  1. Iurii Petrov

    Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden
    Competing interests
    The authors declare that no competing interests exist.
  2. Andrey Alexeyenko

    Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden
    For correspondence
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8812-6481


Vetenskapsrådet (2016-04940)

  • Iurii Petrov

Vetenskapsrådet (2016-04940)

  • Andrey Alexeyenko

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Reviewing Editor

  1. C Daniela Robles-Espinoza, International Laboratory for Human Genome Research, Mexico

Publication history

  1. Received: September 17, 2021
  2. Preprint posted: October 5, 2021 (view preprint)
  3. Accepted: May 20, 2022
  4. Accepted Manuscript published: May 20, 2022 (version 1)
  5. Version of Record published: June 1, 2022 (version 2)


© 2022, Petrov & Alexeyenko

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.


  • 636
    Page views
  • 133
  • 0

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Iurii Petrov
  2. Andrey Alexeyenko
Individualized discovery of rare cancer drivers in global network context
eLife 11:e74010.

Further reading

    1. Cancer Biology
    2. Microbiology and Infectious Disease
    Changkun Hu, Taylor Bugbee ... Nicholas Wallace
    Research Article

    Double strand breaks (DSBs) are one of the most lethal DNA lesions in cells. The E6 protein of beta-human papillomavirus (HPV8 E6) impairs two critical DSB repair pathways; homologous recombination (HR) and non-homologous end-joining (NHEJ). However, HPV8 E6 only delays DSB repair. How DSBs are repaired in cells with HPV8 E6 remains to be studied. We hypothesize that HPV8 E6 promotes a less commonly used DSB repair pathway, alternative end-joining (Alt-EJ). Using CAS9 based Alt-EJ reporters, we show that HPV8 E6 promotes Alt-EJ. Further, using small molecule inhibitors, CRISPR/CAS9 gene knockout, and HPV8 E6 mutant, we find that HPV8 E6 promotes Alt-EJ by binding p300, an acetyltransferase that facilitates DSB repair by HR and NHEJ. At least some of this repair occurs through a subset of Alt-EJ known as polymerase theta dependent end joining. Finally, whole genome sequencing analysis showed HPV8 E6 caused an increased frequency of deletions bearing the microhomology signatures of Alt-EJ. This study fills the knowledge gap of how DSB is repaired in cells with HPV8 E6 and the mutagenic consequences of HPV8 E6 mediated p300 destabilization. Broadly, this study supports the hypothesis that beta-HPV promotes cancer formation by increasing genomic instability.

    1. Cancer Biology
    Anthony Lozano, Francois-Régis Souche ... Damien Grégoire
    Research Article

    Quantitative differences in signal transduction are to date an understudied feature of tumour heterogeneity. The MAPK Erk pathway, which is activated in a large proportion of human tumours, is a prototypic example of distinct cell fates being driven by signal intensity. We have used primary hepatocyte precursors transformed with different dosages of an oncogenic form of Ras to model subclonal variations in MAPK signalling. Orthotopic allografts of Ras-transformed cells in immunocompromised mice gave rise to fast-growing aggressive tumours, both at the primary location and in the peritoneal cavity. Fluorescent labelling of cells expressing different oncogene levels, and consequently varying levels of MAPK Erk activation, highlighted the selection processes operating at the two sites of tumour growth. Indeed, significantly higher Ras expression was observed in primary as compared to secondary, metastatic sites, despite the apparent evolutionary trade-off of increased apoptotic death in the liver that correlated with high Ras dosage. Analysis of the immune tumour microenvironment at the two locations suggests that fast peritoneal tumour growth in the immunocompromised setting is abrogated in immunocompetent animals due to efficient antigen presentation by peritoneal dendritic cells. Furthermore, our data indicate that, in contrast to the metastatic-like outgrowth, strong MAPK signalling is required in the primary liver tumours to resist elimination by NK cells. Overall, this study describes a quantitative aspect of tumour heterogeneity and points to a potential vulnerability of a subtype of hepatocellular carcinoma as a function of MAPK Erk signalling intensity.