Uncovering hidden cancer self-dependencies through analysis of shRNA-level dependency scores

Large-scale short hairpin RNA (shRNA) screens on well-characterized human cancer cell lines have been widely used to identify novel cancer dependencies. However, the off-target effects of shRNA reagents pose a significant challenge in the analysis of these screens. To mitigate these off-target effects, various approaches have been proposed that aggregate different shRNA viability scores targeting a gene into a single gene-level viability score. Most computational methods for discovering cancer dependencies rely on these gene-level scores. In this paper, we propose a computational method, named NBDep, to find cancer self-dependencies by directly analyzing shRNA-level dependency scores instead of gene-level scores. The NBDep algorithm begins by removing known batch effects of the shRNAs and selecting a subset of concordant shRNAs for each gene. It then uses negative binomial random effects models to statistically assess the dependency between genetic alterations and the viabilities of cell lines by incorporating all shRNA dependency scores of each gene into the model. We applied NBDep to the shRNA dependency scores available at Project DRIVE, which covers 26 different types of cancer. The proposed method identified more well-known and putative cancer genes compared to alternative gene-level approaches in pan-cancer and cancer-specific analyses. Additionally, we demonstrated that NBDep controls type-I error and outperforms statistical tests based on gene-level scores in simulation studies.

Discovering cancer dependencies, including synthetic lethality and self-dependencies, is of crucial significance for identifying new cancer treatment options.Synthetic lethality refers to the interaction of two genes, where the simultaneous loss of function through either genetic events or inhibition results in cell death, but the loss of function of either gene alone does not.Several computational methods have been previously proposed for identifying synthetic lethality interactions using loss-of-function RNAi and CRISPR screens [15][16][17][18][19] .Cancer self-dependency refers to a dependency type in which the loss of function of a specific gene leads to cell death preferentially in cells with specific molecular characteristics, such as mutations in the same gene.Studies have previously examined self-dependency in relation to missense and damaging mutations, as well as copy-number amplification, and have identified novel putative cancer genes through these investigations 9,20,21 .
To identify cancer dependency through gene perturbation, several approaches first aggregate scores at the gene level and then perform statistical tests at this level 9,11,12,22,23 .To the best of our knowledge, in analyzing RNAi screening data, the method proposed in this study is the only one that directly performs statistical tests at the shRNA level to identify cancer dependencies.However, for CRISPR analysis, MAGeCK method performs hypothetical testing at the sgRNA level 24 .Similarly, in protein expression data analysis, it has been shown that analyzing data at the peptide level offers greater sensitivity and specificity [25][26][27] .Additionally, there are methods tested at the shRNA level, such as siMEM, which is a time series-based design in RNAi screening, and Screen-Beam, a meta-analysis method that uses both shRNAs and sgRNAs.However, these methods are not comparable to this research 28,29 .
In this research, we investigated the potential of using shRNA-level viability scores instead of gene-level data to enhance the statistical power for identifying cancer self-dependencies.To tackle the challenge of off-target effects, we developed a statistical method for analyzing dependency scores from perturbation screens at the shRNA level.Our hypothesis was that analyzing at the shRNA level would result in greater statistical power for uncovering hidden cancer dependencies.To evaluate the performance of our proposed method, we analyzed pan-cancer and cancer-specific analyses on the Project DRIVE dataset, an RNAi screening project that used deep coverage shRNA lentiviral libraries to target genes across 398 cell lines, providing a valuable resource for exploring and discovering cancer dependencies.In addition, we conducted simulation studies to evaluate the type-I error and statistical power of our approach across various sample sizes.

Materials and methods shRNA viability data from Project DRIVE
In this research, we used the Project DRIVE data to find novel cancer dependencies.Project DRIVE conducted knockdown experiments in three different pools namely poolA, poolB, and BGPD using 158,114 shRNAs on 9850 genes in 398 cell lines across 39 cancer types by using a median of 20 pooled shRNA per gene.Pools poolA, poolB, and BGPD included 3492, 3577, and 4178 genes, respectively (Supplementary Fig. S1b).The shRNA viability scores were defined using next-generation sequencing counts as log fold change of shRNA read counts 14 days after the onset of the screen compared to shRNA initial abundance in the input library.For cancerspecific analyses, we only considered 26 cancer types that included at least four cell lines, following completion of all preprocessing steps.

Gene-level viability data from Project DRIVE
We used ATARiS and DEMETER2 methods for aggregating shRNA viability scores into gene-level dependency scores.ATARiS uses a subset of shRNAs with a consistent pattern of viability scores across cell lines and provides gene-level dependency scores relative to screened cell lines due to using median-centered shRNA viability scores.The ATARiS scores of Project DRIVE are available for 6557 genes in 398 cell lines.On the contrary, DEMETER2 scores are absolute dependency scores computed using a hierarchical Bayesian inference scheme by explicitly modeling the off-target effects associated with the seed sequence of each shRNA.The DEMETER2 scores of Project DRIVE are available for 7975 genes in 397 cell lines.

Copy number and mutation data
We downloaded the mutation data from the DepMap website and the copy number GISTIC2 data from the CBio-Portal website 30 .We considered three different classes of genetic alterations: non-missense mutations, missense mutations, and copy-number amplification.A gene has a non-missense mutation in a cell line if it harbors any of start codon deletion, stop codon deletion, start codon insertion, start codon insertion, frameshift deletion, frameshift insertion, in-frame deletion, in-frame insertion, nonsense mutation, and splice site.Missense mutations are simply defined as mutations annotated with missense mutation.We used the GISTIC2 to determine the copy-number status of a gene.We specifically used GISTIC2 values 2 and − 2 representing copy number amplification and deep deletion, respectively.Finally, a gene is considered wild-type in a cell line if it does not harbor missense, non-missense, and copy-number amplification and deep deletion.

TargetScan data
We used thermodynamic stability of seed pairing in shRNAs in RNAi screening to incorporate the off-target effects of seed sequences.The thermodynamic stability values of 7-mer seeds of the DRIVE data were extracted from TargetScan 31 .The total number of seeds in the TargetScan data and Project DRIVE are 16,384 ( 4 7 ) and 13,071, respectively.

IntOGen data
We used the IntOGen data (version 2020.02.01) for validation of our method.It consists of 568 genes involved in cancer among which 422 and 474 genes were also available in the ATARiS and DEMETER2 scores, respectively 32 .www.nature.com/scientificreports/ Step 2: batch-effect removal and ranking We then used the removeBatchEffect function in limma package in R 33 to remove two batch effects, namely pool and thermodynamic stability of seed sequences, from the shRNA data.Subsequently, we ranked all the corrected shRNA values per cell line such that shRNAs with lower ranks represented higher depletions.
Step 3: obtaining consistent shRNAs Similar to the ATARiS algorithm, though with a different approach, we found a subset of shRNAs for each gene which exhibited consistent behavior across all cell lines.To achieve this, we calculated Pearson correlation coefficients between all shRNA viability ranks of a gene and the average viability ranks across all shRNAs for all cell lines in a given gene and retained half of the shRNAs with the highest Pearson correlation coefficients.The final preprocessed data to perform NBDep algorithm contains information on 76,495 shRNAs targeting 7241 genes 339 cell lines.
Step 4: statistical testing using a negative binomial mixed effects model Finally, a negative binomial mixed effects model was performed to discover gene drivers in the refined data obtained in the previous steps.To incorporate the heterogeneity of cell lines into the model, they were entered as a random effect covariate in the model.In this model, we assumed that the viability rank for i th shRNA and k th cell line, denoted by R ik , follows a negative binomial distribution for i = 1, 2, . . ., n; k = 1, 2, . . ., l .Parameters µ ik and θ are mean and shape parameter of negative binomial dis- tribution.We selected the negative binomial distribution because it is a two-parameter distribution that offers the necessary flexibility to model non-negative integer values, making it suitable for the ranks.The probability mass function of R ik is given by where Ŵ(.) is the gamma function.The parameter µ ik is defined as: The mutational status of i th shRNA in k th cell line is a fixed effect, denoted as the variable X ik where 0 and 1 indicate wild-type and mutant cell lines, respectively.The variable u k represents the random effect for the kth cell line.Under the assumption of this model, the vector u = (u 1 , u 2 , . . ., u l ) T follows a multivariate normal with mean 0 and covariance matrix (i.e., u ∼ MVN(0, �) ).The maximum likelihood estimates of the model parameters, (θ, β, �) , were obtained using the glmmTMB package in R. To test for the dependency between mutational status and viability score of a gene, we performed one-sided hypothesis testing with H 0 : β = 0 and H 1 : β < 0 .Under the null hypothesis, β follows a normal distribution and consequently the p-value is straightforward to compute.

Simulation study
In order to evaluate the performance of NBDep and gene-level methods, we conducted two simulation studies to assess type I error and statistical power.

Simulation study: type I error analysis
In this simulation study, we selected a set of genes for type I error analysis, including known cancer genes KRAS, NRAS, BRAF, TP53, CDK4, AXIN1, DHX9, TOP2A, RRAS2, COL1A1 as well as a few genes that are not known to be involved in cancer namely UGT8, TRPV3, FUS, BRIX1, NAN60, RPS6KA4, PTBP1, SPRR1A, MYLIP, LRCH4, DUSP26, CSF2RB, CBX4, DLD, MCRS1, PYGB, ATP11B, and BTNL2.We assumed that there were m mutated and w wild-type cell lines with respect to a given gene.To generate 5000 random datasets, we performed the following steps: (1) We randomly selected n m cell lines from the set of all cell lines and assigned them to the mutant group.We generated n m from a uniform distribution between 2 and m. (2) We randomly selected n w cell lines from the remaining cell lines and assigned them to the wild-type group.
We generated n w from a uniform distribution between 2 and w.This ensured that both groups had the same distribution, and the generated data was suitable for assessment of type I error.(3) We then tested for a relationship between mutation status and viability using various methods and calculated type I error as the fraction of rejected tests over the total 5000 tests for each method.

Simulation study: statistical power analysis
In this section, we focused on known cancer genes having more than five mutant cell lines, identified through all methods including APSiC on ATARiS and DEMETER2 scores, as well as NBDep.These genes, namely KRAS, NRAS, PIK3CA, TP53, CTNNB1, and BRAF, were used in analyzing statistical power.To this end, we first simulated datasets in which the distributions of two groups were different.We then subsampled n m mutated cell lines and n w wild-type cell lines, where n m was between 2 and 10, and n w = qn m for each gene.The parameter q represents the ratio of wild-type to mutant cell lines, and we considered values of q = 1, 2, 3, 4. For each gene, www.nature.com/scientificreports/n m and q , we generated 800 datasets.The statistical power was then calculated as the ratio of rejected tests to the total number of tests for each setting.

Statistical testing on gene-level data
We applied the Wilcoxon rank sum test and APSiC, a ranked-based statistical approach based on Irwin-Hall and Bates distributions, to identify self-dependencies from perturbation screens 21 .These tests were conducted on gene-level scores obtained from ATARiS or DEMETER2 dependency scores, which are commonly used for aggregating shRNA viability scores targeting a gene into a single gene-level score.Since APSiC employs a onesided hypothesis testing, we used a one-sided hypothesis testing in NBDep too to facilitate direct comparison with APSiC results on gene-level ATARiS and DEMETER2 scores.We discovered that the statistical power of the Wilcoxon rank sum test was significantly lower than that of APSiC, and as a result, we only included the results of APSiC in our analyses.

PPI enrichment
We conducted an enrichment analysis on the set of genes identified through our pan-cancer analysis using data from the STRING database, which contains information on protein-protein interactions.The goal of this analysis was to determine if the identified genes had a greater number of interactions with each other compared to a set of randomly selected genes 34 .

Multiple testing correction
For multiple testing correction, we controlled the false discovery rate using Benjamini and Hochberg approach 35 .The significance level was set to 0.05.

Simulated studies
We first assessed the type I error of our method on synthetic data (see "Methods").At the significance level of 5%, the estimated type I error of NBDep was 0.067, slightly higher than the expected 5% error.Next, we compared the statistical power of our method with those of APSiC statistical test on ATARiS and DEMETER2 gene-level dependency scores on three well-known missense driver genes (KRAS, NRAS, and PIK3CA) using subsampling approaches.To determine whether our method is more robust to low sample size than other methods, we generated 800 datasets for each sample size resulting from various numbers of mutant cell lines and values of q (see "Materials and methods" section).In general, the statistical power of the NBDep algorithm is higher than other methods for small sample sizes with q = 1 (Fig. 2) and other values of q in the analyzed genes (Supplementary Figs.S2, S3).

Influence of denoising and ranking on NBDep's performance
In this section, we investigated the reason why the NBDep algorithm has enhanced statistical power.In particular, we considered whether it is due to denoising or ranking.To this end, we compared the following approaches: (1) a linear mixed-effects model on raw logFC of shRNA readouts, (2) a linear mixed-effects model on denoised logFC values, (3) a negative binomial mixed-effects model on rank data obtained from raw logFC values, (4) a negative binomial mixed-effects model on rank data obtained from denoised logFC values (i.e., NBDep algorithm).We performed a similar simulation as the previous section on KRAS, NRAS, PIK3CA.The Supplementary Fig. S4 demonstrates that the main reason for the enhanced statistical power is using ranks of logFC values.However, denoising can also slightly improve the statistical power.In order to ensure that the ranking does not result in elevated false positives, we calculated the type one error for the third model, which was 0.058 and plausible.

Pan-cancer analysis
NBDep method identified 89 missense mutational driver genes, 41 amplification driver genes, and 125 non-missense mutational driver genes in pan-cancer analyses.Top identified amplification, missense, and non-missense mutational genes by the NBDep method in pan-cancer are depicted in Fig. 3c.In analysis of missense driver genes, NBDep identified sixteen well-recognized genes reported in IntOGen namely KRAS, NRAS, BRAF, PIK3CA, CTNNB1, TP53, SMAD4, BCL2, SIX1, FBXO11, COL1A1, MAP2K1, TOP2A, DHX9, HRAS (Fig. 3a).In addition, unreported genes in IntOGen such as CANT1, IL20RA, BRX1, BUB1, ATIC, Twist2, and CSF2RB were identified by the NBDep method as putative missense mutational cancer genes www.nature.com/scientificreports/but these genes were not detected by using APSiC on ATARiS and DEMETER2 gene-level scores.Recently, the important role of these genes as driver genes in different stages of cancer has been presented [36][37][38][39][40][41][42] .Box plots depicting the standardized, denoised ranks of four example genes identified by the NBDep algorithm are displayed in Fig. 3b.As depicted in Supplementary Fig. S5, there was a significant PPI enrichment observed for both the identified missense driver genes and the identified genes not reported in IntOGen, with p-values of 0.000214 and 0.00767, respectively.One of the top missense driver genes identified by NBDep method is ANAPC1 which has not been reported in IntOGen and the literature research and it can be introduced as a novel cancer gene.Figure 4b indicates that ANAPC1 has significantly lower ranks in missense cell lines than wild-type ones even before taking account for cell line effect in the negative binomial mixed effects model.On the other hand, as depicted in Fig. 4a, nine identified amplification driver genes by NBDep are in the IntOGen gene list, namely ERBB2, CDK4, FLT3, MAP2K1, KRAS, NRAS, GNAS, RUNX1 as oncogenes and LATS2 as a tumor suppressor.Additionally, other amplification driver genes identified by NBDep, such as ITGAE, ALKBH3, SHPK, and GJB3 have been previously reported to play a significant role in cancer [43][44][45][46] .NBDep also identified 23 non-missense driver genes that were previously reported in IntOGen, namely SIN3A, NTRK1, FN1, MGA, MAP3K1, PTEN, DROSHA, WRN, RELA, CDX2, KDR, CDKN1B, POLD1, PIK3R1, NUP214, PML, GATA3, TOP2A, APC, TP53, www.nature.com/scientificreports/RAD21, GNAI2, and NXF1.Recent studies have also highlighted the role of some of the top non-missense driver genes that were not previously identified by IntOGen, such as ALOX12 47 .Furthermore, NBDep identified two putative non-missense driver genes MEPCE and TMEM41A that were not previously known to be cancer driver genes.
In the pan-cancer analysis of all three alterations, several genes were identified as both non-missense and missense driver genes, including ATIC, MCRS1, TBL3, CACNA1A, TOP2A, and TP53.Moreover, NRAS, KRAS, and MAP2K1 were found to be amplification and missense driver genes.Additionally, KCNG1 was identified as a driver gene in cancer by NBDep but has not been previously reported in both non-missense and amplification driver genes.Figure 4a shows that the NBDep algorithm identified more common genes across all three alterations compared to using the APSiC method on ATARiS and DEMETER2 scores.Figure 4c displays the top 30 genes identified by the NBDep algorithm in all three alterations, along with their p-values computed using the hypergeometric test for finding the statistical significance of overlaps between identified genes and the IntOGen gene list.It is worth noting that the number of identified driver genes by NBDep was higher in cancer types with more cell lines such as colorectal, skin, and pancreatic carcinoma.

Cancer-specific analysis
We employed the NBDep algorithm to identify driver genes across 26 cancer types only on genes with at least two mutant and two wild-type cell lines.Our analysis revealed 87 amplification, 177 missense, and 45 non-missense driver genes across cancer types.The most frequent identified amplification driver genes belonged to pancreatic and ovarian carcinoma.The identified missense mutational genes were found in colorectal carcinoma and ovarian carcinoma, and non-missense mutational genes were observed in colorectal carcinoma.
Among the amplification driver genes, NRAS (in breast carcinoma and skin melanoma), ERBB2 (in gastric carcinoma, esophageal carcinoma), and NFATC2 (in ovarian carcinoma and pancreatic carcinoma) were identified in more than one cancer type.Four genes, namely KRAS (in breast carcinoma, non-small cell lung cancer (NSCLC), lymphoma multiple myeloma carcinoma, colorectal carcinoma, and leukemia), NRAS (in pancreatic carcinoma, bladder carcinoma, primitive neuro-ectodermal tumors (PNET)), PIK3CA (in breast carcinoma and ovarian carcinoma), and ARHGAP31 (in colorectal carcinoma and pancreatic carcinoma) were identified as missense driver genes in more than one cancer type.Notably, the number of identified driver genes by NBDep was higher in cancer types with more cell lines, such as colorectal, skin, and pancreatic carcinoma.
Among the 89 identified missense driver genes in pan-cancer, 21 genes were detected in multiple cancer types, including known cancer genes such as TP53, KRAS, NRAS, and PIK3CA.Although most of these genes were identified in colorectal cancer, PIK3CA, TP53 and CTNNB1 were detected in other cancer types too (breast carcinoma, central nervous system glioma high grade, and gastric carcinoma).KRAS, IL2ORA, CANT1, NRAS and TP53 were the top genes in different cancers.A recent study indicates IL2ORA is an important regulator of oncogenic and immune pathways in colorectal carcinoma 36 (Supplementary Table 1).
NRAS (breast carcinoma and skin melanoma), CDK4 (NSCLC), and ERBB2 (gastric carcinoma and esophageal carcinoma) were identified as the top identified amplification driver genes in different cancer types (Supplementary Table 1).Also, we found three top genes which have not been reported in IntOGen, TNFRSF19, GJA3, and GJB2, recently found by researchers to have main role in colorectal cancer 48,49 (Supplementary Table 1).
Our method introduced SPATA13 a guanine-factor as a novel amplification driver gene in colorectal carcinoma, required for MMP9 up-regulation via the JNK signaling pathway in colorectal tumor cells.Also, NBDep method identified amplification WHSC1 in breast carcinoma as a novel gene.
For the cancer types with small number of cell lines, the NBDep method is able to identify well-known genes.For leukemia, lymphoma NH B-cell, lymphoma multiple myeloma cancer, and NSCLC that have very low number of cell lines, KRAS was detected as the most important missense gene driver.NBDep found MYC as an amplification driver gene in lymphoma multiple myeloma cancer while gene-level approaches did not identified this gene.Our finding was confirmed by a recent study where amplified MYC was shown to be effective in myeloma cancer 50 .In addition, the recent research approves the identified ATP8A1 by NBDep as an important gene in NSCLC 51 .The NBDep method also identified 43 missense driver genes in ovarian cancer having a significant protein-protein network (p-value = 2.12e−06) where well-known genes such as RHOA, PIK3CA, and ATM are among these genes [52][53][54] .Our method proposed PCDH8 and BAZ2A as novel missense driver genes in ovarian carcinoma.
In summary, the NBDep algorithm is able to identify well-known driver genes in cancer-specific analyses with large and small number of cell lines and to introduce novel putative driver genes (Fig. 4).In addition, NBDep algorithm identified more curated genes regarding to IntOGen in cancer-specific than ATARiS and DEMETER2 (Supplementary Fig. S6).

Discussion
Perturbation screens including RNAi screening has become increasingly popular in the field of cancer genomics over the past decade.One major limitation of RNAi screening is the off-targets issue where it poses a major challenge to infer actual gene effects in these screens.Different computational methods, including ATARiS and DEMETER2, have been developed to handle off-target effects of shRNAs, leading to gene-level scores known as dependency scores.ATARiS achieves an aggregate score for each gene by discarding shRNAs with non-consistent In this research, we presented a statistical framework aimed at handling off-target effects at shRNA level to identify driver genes associated to missense, amplification, and non-missense alterations.We applied our approach to 26 cancer-specific types as well as pan-cancer data using the Project DRIVE data.We coped with off target effects of shRNAs by incorporating thermodynamic stability of 7-mer seed shRNAs at the shRNA level as a main batch effect proposed by TargetScan.After calculating the denoised shRNA from the original shRNA scores, we then employed 50% of consistent shRNAs of each gene across cell lines.We subsequently performed a negative binomial mixed effects model to investigate association of gene perturbation and alteration statuses of genes.We additionally showed that NBDep algorithm is robust in small sample sizes and can detect driver genes more effectively than using the APSiC method on ATARiS and DEMETER2 gene-level scores.
Having compared the numbers of shRNAs designed for genes in IntOGen to other genes, we recognized that IntOGen genes have more shRNAs than non-well-known genes.This finding advocates that it may be more efficient to design more reagents for all genes when performing RNAi screening.
Our method was capable of identifying well-known genes in pan-cancer and 26 cancer-specific types including KRAS, NRAS, TP53, BRAF, and PIK3CA in missense alteration, ERBB2, CDK4, and FLT3 as amplification driver genes, and well-recognized non-missense driver genes such as TP53, TOP2A, APC, and GATA3.Additionally, in cancer-specific analyses, our method was able to identify well-known genes such as KRAS and PIK3CA as missense driver genes, NRAS and EGFR as amplification driver genes in breast carcinoma.NBDep method could not identify any driver gene in liver carcinoma.NBDep identified CCKAR and WHSC as novel missense and amplification driver genes, respectively, in breast carcinoma.CCKAR was reported as a driver gene in gallbladder and biliary tract cancer and WHSC as TP53 binding protein 55 .Moreover, PCDH8, a protein coding that acts as a cell adhesion molecule, and BAZ2A having DNA binding activity, were suggested as missense novel genes in ovarian cancer.
In summary, our method utilizes more information of shRNAs in RNAi screening and is capable to identify above-mentioned self-dependencies in pan-cancer and 26 cancer-specific types while handling off-target effects.Our approach also provides denoised shRNA ranks through which it is possible to explore other types of cancer dependencies such as synthetic lethality.

Figure 1 .
Figure 1.The flowchart of the NBDep algorithm.NBDep employs Project DRIVE, CCLE, GISTIC, and thermodynamic stability of shRNA's seeds obtained from TargetScan to identify missense, non-missense, and amplification driver genes.After eliminating batch effects arising from shRNA thermodynamics and pooling, all denoised logFCs in each cell line were ranked.Then, 50% of the most consistent shRNAs of each gene were selected based on their Pearson correlation with the average shRNA profile.Finally, a negative binomial mixed effects model was applied to the shRNA ranks in three alteration classes: missense, non-missense, and amplification alterations.The random effect in this model is associated with the cell line and the fixed effect is a binary variable indicating mutational status.

Figure 2 .
Figure 2. Comparison of statistical power between APSiC method applied to ATARiS and DEMETER2 scores and the NBDep algorithm for KRAS, NRAS, and PIK3CA across various sample sizes and q = 1.

Figure 3 .
Figure 3. (a) The number of identified genes common with IntOGen in pan-cancer analysis resulted from three methods: NBDep as well as applying the APSiC method on ATARiS and DEMETER gene-level scores in three alterations: amplification, missense, and non-missense.(b) Box plots of standardized denoised ranks of one novel missense-driver gene, ANAPC1, one experimentally curated missense driver gene, KRAS, two genes previously reported in literatures, RAB19 and IL20RA, identified by NBDep.(c) Top identified amplification, missense, and non-missense mutational genes by NBDep method in pan-cancer.Red color indicates IntOGen genes.Venn diagrams indicates the number of common genes with IntOGen in each of three alteration groups along with their enrichment p-values calculated by the hypergeometric test.

Figure 4 .
Figure 4.The identified amplification, missense, and non-missense driver genes by NBDep in cancer-specific analyses.The genes highlighted in red are in the IntOGen list.