The SORL1 p.Y1816C variant causes impaired endosomal dimerization and autosomal dominant Alzheimer's disease.

Truncating genetic variants of SORL1, encoding the endosome recycling receptor SORLA, have been accepted as causal of Alzheimer's disease (AD). However, most genetic variants observed in SORL1 are missense variants, for which it is complicated to determine the pathogenicity level because carriers come from pedigrees too small to be informative for penetrance estimations. Here, we describe three unrelated families in which the SORL1 coding missense variant rs772677709, that leads to a p.Y1816C substitution, segregates with Alzheimer's disease. Further, we investigate the effect of SORLA p.Y1816C on receptor maturation, cellular localization, and trafficking in cell-based assays. Under physiological circumstances, SORLA dimerizes within the endosome, allowing retromer-dependent trafficking from the endosome to the cell surface, where the luminal part is shed into the extracellular space (sSORLA). Our results showed that the p.Y1816C mutant impairs SORLA homodimerization in the endosome, leading to decreased trafficking to the cell surface and less sSORLA shedding. These trafficking defects of the mutant receptor can be rescued by the expression of the SORLA 3Fn-minireceptor. Finally, we find that iPSC-derived neurons with the engineered p.Y1816C mutation have enlarged endosomes, a defining cytopathology of AD. Our studies provide genetic as well as functional evidence that the SORL1 p.Y1816C variant is causal for AD. The partial penetrance of the mutation suggests this mutation should be considered in clinical genetic screening of multiplex early-onset AD families.

Connecting dementia risk loci to the CSF proteome identifies pathophysiological leads for dementia.

Genome-wide association studies have successfully identified many genetic risk loci for dementia, but exact biological mechanisms through which genetic risk factors contribute to dementia remains unclear. Integrating CSF proteomic data with dementia risk loci could reveal intermediate molecular pathways connecting genetic variance to the development of dementia. We tested to what extent effects of known dementia risk loci can be observed in CSF levels of 665 proteins [proximity extension-based (PEA) immunoassays] in a deeply-phenotyped mixed memory clinic cohort [n = 502, mean age (standard deviation, SD) = 64.1 (8.7) years, 181 female (35.4%)], including patients with Alzheimer's disease (AD, n = 213), dementia with Lewy bodies (DLB, n = 50) and frontotemporal dementia (FTD, n = 93), and controls (n = 146). Validation was assessed in independent cohorts (n = 99 PEA platform, n = 198, mass reaction monitoring-targeted mass spectroscopy and multiplex assay). We performed additional analyses stratified according to diagnostic status (AD, DLB, FTD and controls separately), to explore whether associations between CSF proteins and genetic variants were specific to disease or not. We identified four AD risk loci as protein quantitative trait loci (pQTL): CR1-CR2 (rs3818361, P = 1.65 × 10-8), ZCWPW1-PILRB (rs1476679, P = 2.73 × 10-32), CTSH-CTSH (rs3784539, P = 2.88 × 10-24) and HESX1-RETN (rs186108507, P = 8.39 × 10-8), of which the first three pQTLs showed direct replication in the independent cohorts. We identified one AD-specific association between a rare genetic variant of TREM2 and CSF IL6 levels (rs75932628, P = 3.90 × 10-7). DLB risk locus GBA showed positive trans effects on seven inter-related CSF levels in DLB patients only. No pQTLs were identified for FTD loci, either for the total sample as for analyses performed within FTD only. Protein QTL variants were involved in the immune system, highlighting the importance of this system in the pathophysiology of dementia. We further identified pQTLs in stratified analyses for AD and DLB, hinting at disease-specific pQTLs in dementia. Dissecting the contribution of risk loci to neurobiological processes aids in understanding disease mechanisms underlying dementia.

Cognitively healthy centenarians are genetically protected against Alzheimer's disease.

BACKGROUND: Alzheimer's disease (AD) prevalence increases with age, yet a small fraction of the population reaches ages > 100 years without cognitive decline. We studied the genetic factors associated with such resilience against AD. METHODS: Genome-wide association studies identified 86 single nucleotide polymorphisms (SNPs) associated with AD risk. We estimated SNP frequency in 2281 AD cases, 3165 age-matched controls, and 346 cognitively healthy centenarians. We calculated a polygenic risk score (PRS) for each individual and investigated the functional properties of SNPs enriched/depleted in centenarians. RESULTS: Cognitively healthy centenarians were enriched with the protective alleles of the SNPs associated with AD risk. The protective effect concentrated on the alleles in/near ANKH, GRN, TMEM106B, SORT1, PLCG2, RIN3, and APOE genes. This translated to >5-fold lower PRS in centenarians compared to AD cases (P = 7.69 × 10-71), and 2-fold lower compared to age-matched controls (P = 5.83 × 10-17). DISCUSSION: Maintaining cognitive health until extreme ages requires complex genetic protection against AD, which concentrates on the genes associated with the endolysosomal and immune systems. HIGHLIGHTS: Cognitively healthy cent enarians are enriched with the protective alleles of genetic variants associated with Alzheimer's disease (AD). The protective effect is concentrated on variants involved in the immune and endolysosomal systems. Combining variants into a polygenic risk score (PRS) translated to > 5-fold lower PRS in centenarians compared to AD cases, and ≈ 2-fold lower compared to middle-aged healthy controls.

snpXplorer: a web application to explore human SNP-associations and annotate SNP-sets.

Genetic association studies are frequently used to study the genetic basis of numerous human phenotypes. However, the rapid interrogation of how well a certain genomic region associates across traits as well as the interpretation of genetic associations is often complex and requires the integration of multiple sources of annotation, which involves advanced bioinformatic skills. We developed snpXplorer, an easy-to-use web-server application for exploring Single Nucleotide Polymorphisms (SNP) association statistics and to functionally annotate sets of SNPs. snpXplorer can superimpose association statistics from multiple studies, and displays regional information including SNP associations, structural variations, recombination rates, eQTL, linkage disequilibrium patterns, genes and gene-expressions per tissue. By overlaying multiple GWAS studies, snpXplorer can be used to compare levels of association across different traits, which may help the interpretation of variant consequences. Given a list of SNPs, snpXplorer can also be used to perform variant-to-gene mapping and gene-set enrichment analysis to identify molecular pathways that are overrepresented in the list of input SNPs. snpXplorer is freely available at https://snpxplorer.net. Source code, documentation, example files and tutorial videos are available within the Help section of snpXplorer and at https://github.com/TesiNicco/snpXplorer.

Resilience and resistance to the accumulation of amyloid plaques and neurofibrillary tangles in centenarians: An age-continuous perspective.

INTRODUCTION: With increasing age, neuropathological substrates associated with Alzheimer's disease (AD) accumulate in brains of cognitively healthy individuals-are they resilient, or resistant to AD-associated neuropathologies? METHODS: In 85 centenarian brains, we correlated NIA (amyloid) stages, Braak (neurofibrillary tangle) stages, and CERAD (neuritic plaque) scores with cognitive performance close to death as determined by Mini-Mental State Examination (MMSE) scores. We assessed centenarian brains against 2131 brains from AD patients, non-AD demented, and non-demented individuals in an age continuum ranging from 16 to 100+ years. RESULTS: With age, brains from non-demented individuals reached the NIA and Braak stages observed in AD patients, while CERAD scores remained lower. In centenarians, NIA stages varied (22.4% were the highest stage 3), Braak stages rarely exceeded stage IV (5.9% were V), and CERAD scores rarely exceeded 2 (4.7% were 3); within these distributions, we observed no correlation with the MMSE (NIA: P = 0.60; Braak: P = 0.08; CERAD: P = 0.16). DISCUSSION: Cognitive health can be maintained despite the accumulation of high levels of AD-related neuropathological substrates. HIGHLIGHTS: Cognitively healthy elderly have AD neuropathology levels similar to AD patients. AD neuropathology loads do not correlate with cognitive performance in centenarians. Some centenarians are resilient to the highest levels of AD neuropathology.

The correlation between neuropathology levels and cognitive performance in centenarians.

INTRODUCTION: Neuropathological substrates associated with neurodegeneration occur in brains of the oldest old. How does this affect cognitive performance? METHODS: The 100-plus Study is an ongoing longitudinal cohort study of centenarians who self-report to be cognitively healthy; post mortem brain donation is optional. In 85 centenarian brains, we explored the correlations between the levels of 11 neuropathological substrates with ante mortem performance on 12 neuropsychological tests. RESULTS: Levels of neuropathological substrates varied: we observed levels up to Thal-amyloid beta phase 5, Braak-neurofibrillary tangle (NFT) stage V, Consortium to Establish a Registry for Alzheimer's Disease (CERAD)-neuritic plaque score 3, Thal-cerebral amyloid angiopathy stage 3, Tar-DNA binding protein 43 (TDP-43) stage 3, hippocampal sclerosis stage 1, Braak-Lewy bodies stage 6, atherosclerosis stage 3, cerebral infarcts stage 1, and cerebral atrophy stage 2. Granulovacuolar degeneration occurred in all centenarians. Some high performers had the highest neuropathology scores. DISCUSSION: Only Braak-NFT stage and limbic-predominant age-related TDP-43 encephalopathy (LATE) pathology associated significantly with performance across multiple cognitive domains. Of all cognitive tests, the clock-drawing test was particularly sensitive to levels of multiple neuropathologies.

Cognition, function, and prevalent dementia in centenarians and near-centenarians: An individual participant data (IPD) meta-analysis of 18 studies.

INTRODUCTION: There are limited data on prevalence of dementia in centenarians and near-centenarians (C/NC), its determinants, and whether the risk of dementia continues to rise beyond 100. METHODS: Participant-level data were obtained from 18 community-based studies (N = 4427) in 11 countries that included individuals ≥95 years. A harmonization protocol was applied to cognitive and functional impairments, and a meta-analysis was performed. RESULTS: The mean age was 98.3 years (SD = 2.67); 79% were women. After adjusting for age, sex, and education, dementia prevalence was 53.2% in women and 45.5% in men, with risk continuing to increase with age. Education (OR 0.95;0.92-0.98) was protective, as was hypertension (odds ratio [OR] 0.51;0.35-0.74) in five studies. Dementia was not associated with diabetes, vision and hearing impairments, smoking, and body mass index (BMI). DISCUSSION: Among the exceptional old, dementia prevalence remains higher in the older participants. Education was protective against dementia, but other factors for dementia-free survival in C/NC remain to be understood.

PLCG2 protective variant p.P522R modulates tau pathology and disease progression in patients with mild cognitive impairment.

A rare coding variant (rs72824905, p.P522R) conferring protection against Alzheimer's disease (AD) was identified in the gene encoding the enzyme phospholipase-C-γ2 (PLCG2) that is highly expressed in microglia. To explore the protective nature of this variant, we employed latent process linear mixed models to examine the association of p.P522R with longitudinal cognitive decline in 3595 MCI patients, and in 10,097 individuals from population-based studies. Furthermore, association with CSF levels of pTau181, total tau, and Aß1-42 was assessed in 1261 MCI patients. We found that MCI patients who carried the p.P522R variant showed a slower rate of cognitive decline compared to non-carriers and that this effect was mediated by lower pTau181 levels in CSF. The effect size of the association of p.P522R with the cognitive decline and pTau181 was similar to that of APOE-ε4, the strongest genetic risk factor for AD. Interestingly, the protective effect of p.P522R was more pronounced in MCI patients with low Aß1-42 levels suggesting a role of PLCG2 in the response to amyloid pathology. In line with this hypothesis, we observed no protective effect of the PLCG2 variant on the cognitive decline in population-based studies probably due to the lower prevalence of amyloid positivity in these samples compared to MCI patients. Concerning the potential biological underpinnings, we identified a network of co-expressed proteins connecting PLCG2 to APOE and TREM2 using unsupervised co-regulatory network analysis. The network was highly enriched for the complement cascade and genes differentially expressed in disease-associated microglia. Our data show that p.P522R in PLCG2 reduces AD disease progression by mitigating tau pathology in the presence of amyloid pathology and, as a consequence, maintains cognitive function. Targeting the enzyme PLCG2 might provide a new therapeutic approach for treating AD.

Correction to: A nonsynonymous mutation in PLCG2 reduces the risk of Alzheimer's disease, dementia with Lewy bodies and frontotemporal dementia, and increases the likelihood of longevity.

The IPDGC (The International Parkinson Disease Genomics Consortium) and EADB (Alzheimer Disease European DNA biobank) are listed correctly as an author to the article, however, they were incorrectly listed more than once.

A nonsynonymous mutation in PLCG2 reduces the risk of Alzheimer's disease, dementia with Lewy bodies and frontotemporal dementia, and increases the likelihood of longevity.

The genetic variant rs72824905-G (minor allele) in the PLCG2 gene was previously associated with a reduced Alzheimer's disease risk (AD). The role of PLCG2 in immune system signaling suggests it may also protect against other neurodegenerative diseases and possibly associates with longevity. We studied the effect of the rs72824905-G on seven neurodegenerative diseases and longevity, using 53,627 patients, 3,516 long-lived individuals and 149,290 study-matched controls. We replicated the association of rs72824905-G with reduced AD risk and we found an association with reduced risk of dementia with Lewy bodies (DLB) and frontotemporal dementia (FTD). We did not find evidence for an effect on Parkinson's disease (PD), amyotrophic lateral sclerosis (ALS) and multiple sclerosis (MS) risks, despite adequate sample sizes. Conversely, the rs72824905-G allele was associated with increased likelihood of longevity. By-proxy analyses in the UK Biobank supported the associations with both dementia and longevity. Concluding, rs72824905-G has a protective effect against multiple neurodegenerative diseases indicating shared aspects of disease etiology. Our findings merit studying the PLCγ2 pathway as drug-target.

The 100-plus Study of cognitively healthy centenarians: rationale, design and cohort description.

Although the incidence of dementia increases exponentially with age, some individuals reach more than 100 years with fully retained cognitive abilities. To identify the characteristics associated with the escape or delay of cognitive decline, we initiated the 100-plus Study ( www.100plus.nl ). The 100-plus Study is an on-going prospective cohort study of Dutch centenarians who self-reported to be cognitively healthy, their first-degree family members and their respective partners. We collect demographics, life history, medical history, genealogy, neuropsychological data and blood samples. Centenarians are followed annually until death. PET-MRI scans and feces donation are optional. Almost 30% of the centenarians agreed to post-mortem brain donation. To date (September 2018), 332 centenarians were included in the study. We analyzed demographic statistics of the first 300 centenarians (25% males) included in the cohort. Centenarians came from higher socio-economic classes and had higher levels of education compared to their birth cohort; alcohol consumption of centenarians was similar, and most males smoked during their lifetime. At baseline, the centenarians had a median MMSE score of 25 points (IQR 22.0-27.5); most centenarians lived independently, retained hearing and vision abilities and were independently mobile. Mortality was associated with cognitive functioning: centenarians with a baseline MMSE score ≥ 26 points had a mortality percentage of 17% per annual year  in the second year after baseline, while centenarians with a baseline MMSE score < 26 points had a mortality of  42% per annual year (p = 0.003). The cohort was 2.1-fold enriched with the neuroprotective APOE-ε2 allele relative to 60-80 year-old population controls (p = 4.8 × 10-7), APOE-ε3 was unchanged and the APOE-ε4 allele was 2.3-fold depleted (p = 6.3 × 10-7). Comprehensive characterization of the 100-plus cohort of cognitively healthy centenarians might reveal protective factors that explain the physiology of long-term preserved cognitive health.

Somatic mutations found in the healthy blood compartment of a 115-yr-old woman demonstrate oligoclonal hematopoiesis.

The somatic mutation burden in healthy white blood cells (WBCs) is not well known. Based on deep whole-genome sequencing, we estimate that approximately 450 somatic mutations accumulated in the nonrepetitive genome within the healthy blood compartment of a 115-yr-old woman. The detected mutations appear to have been harmless passenger mutations: They were enriched in noncoding, AT-rich regions that are not evolutionarily conserved, and they were depleted for genomic elements where mutations might have favorable or adverse effects on cellular fitness, such as regions with actively transcribed genes. The distribution of variant allele frequencies of these mutations suggests that the majority of the peripheral white blood cells were offspring of two related hematopoietic stem cell (HSC) clones. Moreover, telomere lengths of the WBCs were significantly shorter than telomere lengths from other tissues. Together, this suggests that the finite lifespan of HSCs, rather than somatic mutation effects, may lead to hematopoietic clonal evolution at extreme ages.

Unbiased Quantitative Models of Protein Translation Derived from Ribosome Profiling Data.

Translation of RNA to protein is a core process for any living organism. While for some steps of this process the effect on protein production is understood, a holistic understanding of translation still remains elusive. In silico modelling is a promising approach for elucidating the process of protein synthesis. Although a number of computational models of the process have been proposed, their application is limited by the assumptions they make. Ribosome profiling (RP), a relatively new sequencing-based technique capable of recording snapshots of the locations of actively translating ribosomes, is a promising source of information for deriving unbiased data-driven translation models. However, quantitative analysis of RP data is challenging due to high measurement variance and the inability to discriminate between the number of ribosomes measured on a gene and their speed of translation. We propose a solution in the form of a novel multi-scale interpretation of RP data that allows for deriving models with translation dynamics extracted from the snapshots. We demonstrate the usefulness of this approach by simultaneously determining for the first time per-codon translation elongation and per-gene translation initiation rates of Saccharomyces cerevisiae from RP data for two versions of the Totally Asymmetric Exclusion Process (TASEP) model of translation. We do this in an unbiased fashion, by fitting the models using only RP data with a novel optimization scheme based on Monte Carlo simulation to keep the problem tractable. The fitted models match the data significantly better than existing models and their predictions show better agreement with several independent protein abundance datasets than existing models. Results additionally indicate that the tRNA pool adaptation hypothesis is incomplete, with evidence suggesting that tRNA post-transcriptional modifications and codon context may play a role in determining codon elongation rates.

Hi-C Chromatin Interaction Networks Predict Co-expression in the Mouse Cortex.

The three dimensional conformation of the genome in the cell nucleus influences important biological processes such as gene expression regulation. Recent studies have shown a strong correlation between chromatin interactions and gene co-expression. However, predicting gene co-expression from frequent long-range chromatin interactions remains challenging. We address this by characterizing the topology of the cortical chromatin interaction network using scale-aware topological measures. We demonstrate that based on these characterizations it is possible to accurately predict spatial co-expression between genes in the mouse cortex. Consistent with previous findings, we find that the chromatin interaction profile of a gene-pair is a good predictor of their spatial co-expression. However, the accuracy of the prediction can be substantially improved when chromatin interactions are described using scale-aware topological measures of the multi-resolution chromatin interaction network. We conclude that, for co-expression prediction, it is necessary to take into account different levels of chromatin interactions ranging from direct interaction between genes (i.e. small-scale) to chromatin compartment interactions (i.e. large-scale).

Scale-space measures for graph topology link protein network architecture to function.

MOTIVATION: The network architecture of physical protein interactions is an important determinant for the molecular functions that are carried out within each cell. To study this relation, the network architecture can be characterized by graph topological characteristics such as shortest paths and network hubs. These characteristics have an important shortcoming: they do not take into account that interactions occur across different scales. This is important because some cellular functions may involve a single direct protein interaction (small scale), whereas others require more and/or indirect interactions, such as protein complexes (medium scale) and interactions between large modules of proteins (large scale). RESULTS: In this work, we derive generalized scale-aware versions of known graph topological measures based on diffusion kernels. We apply these to characterize the topology of networks across all scales simultaneously, generating a so-called graph topological scale-space. The comprehensive physical interaction network in yeast is used to show that scale-space based measures consistently give superior performance when distinguishing protein functional categories and three major types of functional interactions-genetic interaction, co-expression and perturbation interactions. Moreover, we demonstrate that graph topological scale spaces capture biologically meaningful features that provide new insights into the link between function and protein network architecture. AVAILABILITY AND IMPLEMENTATION: Matlab(TM) code to calculate the scale-aware topological measures (STMs) is available at http://bioinformatics.tudelft.nl/TSSA

FluG affects secretion in colonies of Aspergillus niger.

Colonies of Aspergillus niger are characterized by zonal heterogeneity in growth, sporulation, gene expression and secretion. For instance, the glucoamylase gene glaA is more highly expressed at the periphery of colonies when compared to the center. As a consequence, its encoded protein GlaA is mainly secreted at the outer part of the colony. Here, multiple copies of amyR were introduced in A. niger. Most transformants over-expressing this regulatory gene of amylolytic genes still displayed heterogeneous glaA expression and GlaA secretion. However, heterogeneity was abolished in transformant UU-A001.13 by expressing glaA and secreting GlaA throughout the mycelium. Sequencing the genome of UU-A001.13 revealed that transformation had been accompanied by deletion of part of the fluG gene and disrupting its 3' end by integration of a transformation vector. Inactivation of fluG in the wild-type background of A. niger also resulted in breakdown of starch under the whole colony. Asexual development of the ∆fluG strain was not affected, unlike what was previously shown in Aspergillus nidulans. Genes encoding proteins with a signal sequence for secretion, including part of the amylolytic genes, were more often downregulated in the central zone of maltose-grown ∆fluG colonies and upregulated in the intermediate part and periphery when compared to the wild-type. Together, these data indicate that FluG of A. niger is a repressor of secretion.

An algorithm-based topographical biomaterials library to instruct cell fate.

It is increasingly recognized that material surface topography is able to evoke specific cellular responses, endowing materials with instructive properties that were formerly reserved for growth factors. This opens the window to improve upon, in a cost-effective manner, biological performance of any surface used in the human body. Unfortunately, the interplay between surface topographies and cell behavior is complex and still incompletely understood. Rational approaches to search for bioactive surfaces will therefore omit previously unperceived interactions. Hence, in the present study, we use mathematical algorithms to design nonbiased, random surface features and produce chips of poly(lactic acid) with 2,176 different topographies. With human mesenchymal stromal cells (hMSCs) grown on the chips and using high-content imaging, we reveal unique, formerly unknown, surface topographies that are able to induce MSC proliferation or osteogenic differentiation. Moreover, we correlate parameters of the mathematical algorithms to cellular responses, which yield novel design criteria for these particular parameters. In conclusion, we demonstrate that randomized libraries of surface topographies can be broadly applied to unravel the interplay between cells and surface topography and to find improved material surfaces.

Association of Polygenic Risk Score With Lifetime Risk of Developing Multiple Sclerosis in a Population-Based Birth-Year Cohort.

BACKGROUND AND OBJECTIVES: More than 200 genetic variants have been associated with multiple sclerosis (MS) susceptibility. However, it is unclear to what extent genetic factors influence lifetime risk of MS. Using a population-based birth-year cohort, we investigate the effect of genetics on lifetime risk of MS. METHODS: In the Project Y study, we tracked down almost all persons with MS (pwMS) from birth year 1966 in the Netherlands. As control participants, we included non-MS participants from the Project Y cohort (born 1965-1967 in the Netherlands) and non-MS participants from the Amsterdam Dementia Cohort born between 1963 and 1969. Genetic variants associated with MS were determined in pwMS and control participants using genotyping or imputation methods. Polygenic risk scores (PRSs) based on variants and weights from the largest genetic study in MS were calculated for each participant and assigned into deciles based on the PRS distribution in the control participants. We examined the lifetime risk for each decile and the association between PRS and MS disease variables, including age at onset and time to secondary progression. RESULTS: MS-PRS was calculated for 285 pwMS (mean age 53.0 ± 0.9 years, 72.3% female) and 267 control participants (mean age 51.8 ± 3.2 years, 58.1% female). Based on the lifetime risk estimation, we observed that 1:2,739 of the women with the lowest 30% genetic risk developed MS, whereas 1:92 of the women with the top 10% highest risk developed MS. For men, only 1:7,900 developed MS in the lowest 30% genetic risk group, compared with 1:293 men with the top 10% genetic risk. The PRS was not significantly associated with age at onset and time to secondary progression in both sexes. DISCUSSION: Our results show that the lifetime risk of MS is strongly influenced by genetic factors. Our findings have the potential to support diagnostic certainty in individuals with suspected MS: a high PRS could strengthen a diagnosis, but especially a PRS from the lowest tail of the PRS distribution should be considered a red flag and could prevent misdiagnosing conditions that mimic MS.

Detecting recurrent gene mutation in interaction network context using multi-scale graph diffusion.

BACKGROUND: Delineating the molecular drivers of cancer, i.e. determining cancer genes and the pathways which they deregulate, is an important challenge in cancer research. In this study, we aim to identify pathways of frequently mutated genes by exploiting their network neighborhood encoded in the protein-protein interaction network. To this end, we introduce a multi-scale diffusion kernel and apply it to a large collection of murine retroviral insertional mutagenesis data. The diffusion strength plays the role of scale parameter, determining the size of the network neighborhood that is taken into account. As a result, in addition to detecting genes with frequent mutations in their genomic vicinity, we find genes that harbor frequent mutations in their interaction network context. RESULTS: We identify densely connected components of known and putatively novel cancer genes and demonstrate that they are strongly enriched for cancer related pathways across the diffusion scales. Moreover, the mutations in the clusters exhibit a significant pattern of mutual exclusion, supporting the conjecture that such genes are functionally linked. Using multi-scale diffusion kernel, various infrequently mutated genes are found to harbor significant numbers of mutations in their interaction network neighborhood. Many of them are well-known cancer genes. CONCLUSIONS: The results demonstrate the importance of defining recurrent mutations while taking into account the interaction network context. Importantly, the putative cancer genes and networks detected in this study are found to be significant at different diffusion scales, confirming the necessity of a multi-scale analysis.